Showing 1–50 of 245 results for author: Perna, R

Time-dependent models of AGN disks with radiation from embedded stellar-mass black holes

Authors: Marguerite Epstein-Martin, Hiromichi Tagawa, Zoltan Haiman, Rosalba Perna

Abstract: The brightest steady sources of radiation in the universe, active galactic nuclei (AGN), are powered by gas accretion onto a central supermassive black hole (SMBH). The large sizes and accretion rates implicated in AGN accretion disks are expected to lead to gravitational instability and fragmentation, effectively cutting off mass inflow to the SMBH. Radiative feedback from disk-embedded stars has… ▽ More The brightest steady sources of radiation in the universe, active galactic nuclei (AGN), are powered by gas accretion onto a central supermassive black hole (SMBH). The large sizes and accretion rates implicated in AGN accretion disks are expected to lead to gravitational instability and fragmentation, effectively cutting off mass inflow to the SMBH. Radiative feedback from disk-embedded stars has been invoked to yield marginally stable, steady-state solutions in the outer disks. Here, we examine the consequences of this star formation with a semi-analytical model in which stellar-mass black hole (sBH) remnants in the disk provide an additional source of stabilizing radiative feedback. Assuming star formation seeds the embedded sBH population, we model the time-evolving feedback from both stars and the growing population of accreting sBHs. We find that in the outer disk, the luminosity of the sBHs quickly dominates that of their parent stars. However, because sBHs consume less gas than stars to stabilize the disk, the presence of the sBHs enhances the mass flux to the inner disk. As a result, star formation persists over the lifetime of the AGN, damped in the outer disk, but amplified in a narrow ring in the inner disk. Heating from the embedded sBHs significantly modifies the disk's temperature profile and hardens its spectral energy distribution, and direct emission from the sBHs adds a new hard X ray component, resembling a Compton reflection "hump". △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 21 pages, 15 figures

Ephemeral Flaring Transients Following Supernova Explosions in Black-Hole Binary Systems

Authors: Davide Lazzati, Rosalba Perna, Taeho Ryu

Abstract: The majority of massive stars are members of binary systems, a fraction of which is expected to remain bound after the first star goes off as a supernova. When the second star also explodes, the SN ejecta are bound to interact with the compact object companion. We explore the consequences of this interaction in the case of a black hole (BH) companion. We show that accretion of the SN ejecta by the… ▽ More The majority of massive stars are members of binary systems, a fraction of which is expected to remain bound after the first star goes off as a supernova. When the second star also explodes, the SN ejecta are bound to interact with the compact object companion. We explore the consequences of this interaction in the case of a black hole (BH) companion. We show that accretion of the SN ejecta by the BH generally occurs with highly super-Eddington rates, which can be strongly modulated if the ejecta are clumpy, as typically observed in supernova remnants. This late accretion produces transient flaring activity with a time delay from the SN explosion of days to months, depending on the orbital separation and the velocity of the ejecta. The flares are expected to have a non-thermal, broad band spectrum, but their high-frequency emission (UV and X-rays) would be absorbed within the remnant. Observed flares should therefore be in the optical and near-infrared. We propose this model as an explanation to the late-time flaring activity observed in the fast blue optical transient AT2022tsd. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: Submitted to ApJL, comments most welcome

The role of Ohmic dissipation of internal currents on Hot Jupiter radii

Authors: Taner Akgün, Clàudia Soriano-Guerrero, Albert Elias-López, Daniele Viganò, Rosalba Perna, Fabio Del Sordo

Abstract: The inflated radii observed in hundreds of Hot Jupiters represent a long-standing open issue. The observed correlation between radii and irradiation strength, and the occasional extreme cases, nearly double the size of Jupiter, remain without a comprehensive quantitative explanation. In this investigation, we delve into this issue within the framework of Ohmic dissipation, one of the most promisin… ▽ More The inflated radii observed in hundreds of Hot Jupiters represent a long-standing open issue. The observed correlation between radii and irradiation strength, and the occasional extreme cases, nearly double the size of Jupiter, remain without a comprehensive quantitative explanation. In this investigation, we delve into this issue within the framework of Ohmic dissipation, one of the most promising mechanisms for explaining the radius anomaly. Using the evolutionary code MESA, we simulate the evolution of irradiated giant planets, spanning the range 1 to 8 Jupiter masses, incorporating an internal source of Ohmic dissipation located beneath the radiative-convective boundary. Our modeling is based on physical parameters, and accounts for the approximated conductivity and the evolution of the magnetic fields, utilizing widely-used scaling laws. We compute the radius evolution across a spectrum of masses and equilibrium temperatures, considering varying amounts of Ohmic dissipation, calculated with the internal conductivity profile and an effective parametrization of the currents, based on the typical radius of curvature of the field lines. Our analysis reveals that this internal Ohmic dissipation can broadly reproduce the range of observed radii using values of radius of curvature up to about one order of magnitude lower than what we estimate from the Juno measurements of the Jovian magnetosphere and from MHD dynamo simulations presented herein. The observed trend with equilibrium temperature can be explained if the highly-irradiated planets have more intense and more small-scale magnetic fields. This suggests the possibility of an interplay between atmospherically induced currents and the interior, via turbulence, in agreement with recent box simulations of turbulent MHD in atmospheric columns. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: 13 pages, 14 figures

Tidal Disruption Encores

Authors: Taeho Ryu, Rosalba Perna, Matteo Cantiello

Abstract: Nuclear star clusters (NSCs), made up of a dense concentrations of stars and the compact objects they leave behind, are ubiquitous in the central regions of galaxies, surrounding the central supermassive black hole (SMBH). Close interactions between stars and stellar-mass black holes (sBH) lead to tidal disruption events (TDEs). We uncover an interesting new phenomenon: For a subset of these, the… ▽ More Nuclear star clusters (NSCs), made up of a dense concentrations of stars and the compact objects they leave behind, are ubiquitous in the central regions of galaxies, surrounding the central supermassive black hole (SMBH). Close interactions between stars and stellar-mass black holes (sBH) lead to tidal disruption events (TDEs). We uncover an interesting new phenomenon: For a subset of these, the unbound debris (to the sBH) remain bound to the SMBH, accreting at a later time, and thus giving rise to a second flare. We compute the rate of such events, and find them ranging within $10^{-6}$ -$10^{-3}$ yr$^{-1}$gal$^{-1}$ for SMBH mass $\simeq 10^{6}-10^{9}M_\odot$. Time delays between the two flares spread over a wide range, from less than a year to hundreds of years. The temporal evolution of the light curves of the second flare can vary between the standard $t^{-5/3}$ power-law to much steeper decays, providing a natural explanation for observed light curves in tension with the classical TDE model. Our predictions have implications for learning about NSC properties and calibrating its sBH population. Some double flares may be electromagnetic counterparts to LISA Extreme-Mass-Ratio-Inspiral (EMRI) sources. Another important implication is the possible existence of TDE-like events in very massive SMBHs, where TDEs are not expected. Such flares can affect spin measurements relying on TDEs in the upper SMBH range. △ Less

Submitted 27 March, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJL

3D code for MAgneto-Thermal evolution in Isolated Neutron Stars, MATINS: thermal evolution and lightcurves

Authors: Stefano Ascenzi, Daniele Viganò, Clara Dehman, José A. Pons, Nanda Rea, Rosalba Perna

Abstract: The thermal evolution of isolated neutron stars is a key element in unraveling their internal structure and composition and establishing evolutionary connections among different observational subclasses. Previous studies have predominantly focused on one-dimensional or axisymmetric two-dimensional models. In this study, we present the thermal evolution component of the novel three-dimensional magn… ▽ More The thermal evolution of isolated neutron stars is a key element in unraveling their internal structure and composition and establishing evolutionary connections among different observational subclasses. Previous studies have predominantly focused on one-dimensional or axisymmetric two-dimensional models. In this study, we present the thermal evolution component of the novel three-dimensional magnetothermal code MATINS (MAgneto-Thermal evolution of Isolated Neutron Star). MATINS employs a finite volume scheme and integrates a realistic background structure, along with state-of-the-art microphysical calculations for the conductivities, neutrino emissivities, heat capacity, and superfluid gap models. This paper outlines the methodology employed to solve the thermal evolution equations in MATINS, along with the microphysical implementation which is essential for the thermal component. We test the accuracy of the code and present simulations with non-evolving magnetic fields of different topologies to produce temperature maps of the neutron star surface. Additionally, for a specific magnetic field topology, we show one fully coupled evolution of magnetic field and temperature. Subsequently, we use a ray-tracing code to link the neutron star surface temperature maps obtained by MATINS with the phase-resolved spectra and pulsed profiles that would be detected by distant observers. This study, together with our previous article focused on the magnetic formalism, presents in detail the most advanced evolutionary code for isolated neutron stars, with the aim of comparison with their timing properties, thermal luminosities and the associated X-ray light curves. △ Less

Submitted 28 January, 2024; originally announced January 2024.

Comments: 22 pages, 13 figures; submitted to MNRAS

Shock cooling and breakout emission for optical flares associated with gravitational wave events

Authors: Hiromichi Tagawa, Shigeo S. Kimura, Zoltán Haiman, Rosalba Perna, Imre Bartos

Abstract: The astrophysical origin of stellar-mass black hole (BH) mergers discovered through gravitational waves (GWs) is widely debated. Mergers in the disks of active galactic nuclei (AGN) represent promising environments for at least a fraction of these events, with possible observational clues in the GW data. An additional clue to unveil AGN merger environments is provided by possible electromagnetic e… ▽ More The astrophysical origin of stellar-mass black hole (BH) mergers discovered through gravitational waves (GWs) is widely debated. Mergers in the disks of active galactic nuclei (AGN) represent promising environments for at least a fraction of these events, with possible observational clues in the GW data. An additional clue to unveil AGN merger environments is provided by possible electromagnetic emission from post-merger accreting BHs. Associated with BH mergers in AGN disks, emission from shocks emerging around jets launched by accreting merger remnants is expected. In this paper we compute the properties of the emission produced during breakout and the subsequent adiabatic expansion phase of the shocks, and we then apply this model to optical flares suggested to be possibly associated with GW events. We find that the majority of the reported flares can be explained by the breakout and the shock cooling emission. If these events are real, then the merging locations of binaries are constrained depending on the emission processes. If the optical flares are produced by shock cooling emission, they would display moderate color evolution, possibly color variations among different events, a positive correlation between the delay time and the duration of flares, and accompanying breakout emission in X-ray bands before the optical flares. If the breakout emission dominates the observed lightcurve, it is expected that the color is distributed in a narrow range in the optical band, and the delay time from GW to electromagnetic emission is longer than $\sim 2$ days. Hence, further explorations of the distributions of delay times, color evolution of the flares, and associated X-ray emission will be useful to test the proposed emission model for the observed flares. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: 15 pages, 3 figures

Floating binary planets from ejections during close stellar encounters

Authors: Yihan Wang, Rosalba Perna, Zhaohuan Zhu

Abstract: The discovery of planetary systems beyond our solar system has challenged established theories of planetary formation. Planetary orbits display a variety of unexpected architectures, and free-floating planets appear ubiquitous. The recent detection of candidate Jupiter Mass Binary Objects (JuMBOs) by the James Webb Space Telescope (JWST) has added another puzzling layer. Here, through direct few-b… ▽ More The discovery of planetary systems beyond our solar system has challenged established theories of planetary formation. Planetary orbits display a variety of unexpected architectures, and free-floating planets appear ubiquitous. The recent detection of candidate Jupiter Mass Binary Objects (JuMBOs) by the James Webb Space Telescope (JWST) has added another puzzling layer. Here, through direct few-body simulations, we demonstrate that JuMBOs could arise from the ejection of double giant planets following a close encounter with a passing star, if the two planets are nearly aligned at the closest approach. These ejected JuMBOs typically possess an average semi-major axis approximately three times the orbital separation within their original planetary system and a high eccentricity, characterized by a superthermal distribution that sets them apart from those formed primordially. We estimate the JuMBO formation rate per planetary system in typical and densely populated clusters, revealing a significant environmental dependence. In dense clusters, this formation rate can reach a few percent for wide planetary systems. Comparative analysis of JuMBO rates and properties with current and forthcoming JWST observations across various environments promises insights into the conditions under which these giant planets formed in protoplanetary disks, thereby imposing constraints on theories of giant planet formation. △ Less

Submitted 2 May, 2024; v1 submitted 9 October, 2023; originally announced October 2023.

Comments: Matches version published in Nature Astronomy; Videos are available at https://yihanwangastro.github.io/#posts2

Tidal Disruption Events from three-body scatterings and eccentricity pumping in the disks of Active Galactic Nuclei

Authors: Chaitanya Prasad, Yihan Wang, Rosalba Perna, K. E. Saavik Ford, Barry McKernan

Abstract: Tidal Disruption Events (TDEs) are routinely observed in quiescent galaxies, as stars from the nuclear star cluster are scattered into the loss cone of the central supermassive black hole (SMBH). TDEs are also expected to occur in Active Galactic Nuclei (AGN), due to scattering or orbital eccentricity pumping of stars embedded in the innermost regions of the AGN accretion disk. Encounters with emb… ▽ More Tidal Disruption Events (TDEs) are routinely observed in quiescent galaxies, as stars from the nuclear star cluster are scattered into the loss cone of the central supermassive black hole (SMBH). TDEs are also expected to occur in Active Galactic Nuclei (AGN), due to scattering or orbital eccentricity pumping of stars embedded in the innermost regions of the AGN accretion disk. Encounters with embedded stellar-mass black holes (BH) can result in AGN $μ$TDEs. AGN TDEs and $μ$TDEs could therefore account for a fraction of observed AGN variability. Here, by performing scattering experiments with the few-body code {\tt SpaceHub}, we compute the probability of AGN TDEs and $μ$TDEs as a result of 3-body interactions between stars and binary BHs. We find that AGN TDEs are more probable during the early life of the AGNs, when rates are $\sim (6\times 10^{-5}-5 \times 10^{-2}) (f_\bullet/0.01)$ $~\rm{AGN}^{-1}$~yr$^{-1}$ (where $f_\bullet$ is the ratio between the number density of BHs and stars), generally higher than in quiescent galactic nuclei. By contrast, $μ$TDEs should occur throughout the AGN lifetime at a rate of $\sim (1\times 10^{-4} - 4\times 10^{-2} (f_\bullet/0.01)$ $~\rm{AGN}^{-1}$~yr$^{-1}$. Detection and characterization of AGN TDEs and $μ$ AGN TDEs with future surveys using {\em Rubin} and {\em Roman} will help constrain the populations of stars and compact objects embedded in AGN disks, a key input for the LVK AGN channel. △ Less

Submitted 23 May, 2024; v1 submitted 26 September, 2023; originally announced October 2023.

Comments: Final version accepted by MNRAS

Close Encounters of Star - Black Hole Binaries with Single Stars

Authors: Taeho Ryu, Selma de Mink, Rob Farmer, Ruediger Pakmor, Rosalba Perna, Volker Springel

Abstract: Multi-body dynamical interactions of binaries with other objects are one of the main driving mechanisms for the evolution of star clusters. It is thus important to bring our understanding of three-body interactions beyond the commonly employed point-particle approximation. To this end we here investigate the hydrodynamics of three-body encounters between star-black hole (BH) binaries and single st… ▽ More Multi-body dynamical interactions of binaries with other objects are one of the main driving mechanisms for the evolution of star clusters. It is thus important to bring our understanding of three-body interactions beyond the commonly employed point-particle approximation. To this end we here investigate the hydrodynamics of three-body encounters between star-black hole (BH) binaries and single stars, focusing on the identification of final outcomes and their long-term evolution and observational properties, using the moving-mesh hydrodynamics code AREPO. This type of encounters produces five types of outcomes: stellar disruption, stellar collision, weak perturbation of the original binary, binary member exchange, and triple formation. The two decisive parameters are the binary phase angle, which determines which two objects meet at the first closest approach, and the impact parameter, which sets the boundary between violent and non-violent interactions. When the impact parameter is smaller than the semimajor axis of the binary, tidal disruptions and star-BH collisions frequently occur when the BH and the incoming star first meet, while the two stars mostly merge when the two stars meet first instead. In both cases, the BHs accrete from an accretion disk at super-Eddington rates, possibly generating flares luminous enough to be observed. The stellar collision products either form a binary with the BH or remain unbound to the BH. Upon collision, the merged stars are hotter and larger than main sequence stars of the same mass at similar age. Even after recovering their thermal equilibrium state, stellar collision products, if isolated, would remain hotter and brighter than main sequence stars until becoming giants. △ Less

Submitted 10 July, 2023; v1 submitted 6 July, 2023; originally announced July 2023.

Comments: 16 pages, 10 figures, 2 tables. Submitted to MNRAS. Movie here: https://studio.youtube.com/playlist/PLxLK3qI02cQf_x96DSKQWVS251Jog7IOE/videos

Large-scale Evolution of Seconds-long Relativistic Jets from Black Hole-Neutron Star Mergers

Authors: Ore Gottlieb, Danat Issa, Jonatan Jacquemin-Ide, Matthew Liska, Francois Foucart, Alexander Tchekhovskoy, Brian D. Metzger, Eliot Quataert, Rosalba Perna, Daniel Kasen, Matthew D. Duez, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: We present the first numerical simulations that track the evolution of a black hole-neutron star (BH-NS) merger from pre-merger to $r\gtrsim10^{11}\,{\rm cm}$. The disk that forms after a merger of mass ratio $q=2$ ejects massive disk winds ($3-5\times10^{-2}\,M_{\odot}$). We introduce various post-merger magnetic configurations and find that initial poloidal fields lead to jet launching shortly a… ▽ More We present the first numerical simulations that track the evolution of a black hole-neutron star (BH-NS) merger from pre-merger to $r\gtrsim10^{11}\,{\rm cm}$. The disk that forms after a merger of mass ratio $q=2$ ejects massive disk winds ($3-5\times10^{-2}\,M_{\odot}$). We introduce various post-merger magnetic configurations and find that initial poloidal fields lead to jet launching shortly after the merger. The jet maintains a constant power due to the constancy of the large-scale BH magnetic flux until the disk becomes magnetically arrested (MAD), where the jet power falls off as $L_j\sim t^{-2}$. All jets inevitably exhibit either excessive luminosity due to rapid MAD activation when the accretion rate is high or excessive duration due to delayed MAD activation compared to typical short gamma-ray bursts (sGRBs). This provides a natural explanation for long sGRBs such as GRB 211211A but also raises a fundamental challenge to our understanding of jet formation in binary mergers. One possible implication is the necessity of higher binary mass ratios or moderate BH spins to launch typical sGRB jets. For post-merger disks with a toroidal magnetic field, dynamo processes delay jet launching such that the jets break out of the disk winds after several seconds. We show for the first time that sGRB jets with initial magnetization $σ_0>100$ retain significant magnetization ($σ\gg1$) at $r>10^{10}\,{\rm cm}$, emphasizing the importance of magnetic processes in the prompt emission. The jet-wind interaction leads to a power-law angular energy distribution by inflating an energetic cocoon whose emission is studied in a companion paper. △ Less

Submitted 18 August, 2023; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: For movies of the simulations, see https://oregottlieb.com/bhns.html

Hours-long Near-UV/Optical Emission from Mildly Relativistic Outflows in Black Hole-Neutron Star Mergers

Authors: Ore Gottlieb, Danat Issa, Jonatan Jacquemin-Ide, Matthew Liska, Alexander Tchekhovskoy, Francois Foucart, Daniel Kasen, Rosalba Perna, Eliot Quataert, Brian D. Metzger

Abstract: The ongoing LIGO-Virgo-KAGRA observing run O4 provides an opportunity to discover new multi-messenger events, including binary neutron star (BNS) mergers such as GW170817, and the highly anticipated first detection of a multi-messenger black hole-neutron star (BH-NS) merger. While BNS mergers were predicted to exhibit early optical emission from mildly relativistic outflows, it has remained uncert… ▽ More The ongoing LIGO-Virgo-KAGRA observing run O4 provides an opportunity to discover new multi-messenger events, including binary neutron star (BNS) mergers such as GW170817, and the highly anticipated first detection of a multi-messenger black hole-neutron star (BH-NS) merger. While BNS mergers were predicted to exhibit early optical emission from mildly relativistic outflows, it has remained uncertain whether the BH-NS merger ejecta provides the conditions for similar signals to emerge. We present the first modeling of early near-ultraviolet/optical emission from mildly relativistic outflows in BH-NS mergers. Adopting optimal binary properties: a mass ratio of $q=2$ and a rapidly rotating BH, we utilize numerical relativity and general relativistic magnetohydrodynamic (GRMHD) simulations to follow the binary's evolution from pre-merger to homologous expansion. We use an M1 neutrino transport GRMHD simulation to self-consistently estimate the opacity distribution in the outflows and find a bright near-ultraviolet/optical signal that emerges due to jet-powered cocoon cooling emission, outshining the kilonova emission at early time. The signal peaks at an absolute magnitude of $\sim -15$ a few hours after the merger, longer than previous estimates, which did not consider the first principles-based jet launching. By late 2024, the Rubin Observatory will have the capability to track the entire signal evolution or detect its peak up to distances of $\gtrsim1$ Gpc. In 2026, ULTRASAT will conduct all-sky surveys within minutes, detecting some of these events within $\sim 200$ Mpc. The BH-NS mergers with higher mass ratios or lower BH spins would produce shorter and fainter signals. △ Less

Submitted 8 August, 2023; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: For movies of the simulations, see https://oregottlieb.com/bhns.html

Magnetar central engines in gamma-ray busts follow the universal relation of accreting magnetic stars

Authors: Simone Dall'Osso, Giulia Stratta, Rosalba Perna, Giovanni de Cesare, Luigi Stella

Abstract: Gamma-ray bursts (GRBs), both long and short, are explosive events whose inner engine is generally expected to be a black hole or a highly magnetic neutron star (magnetar) accreting high density matter. Recognizing the nature of GRB central engines, and in particular the formation of neutron stars (NSs), is of high astrophysical significance. A possible signature of NSs in GRBs is the presence of… ▽ More Gamma-ray bursts (GRBs), both long and short, are explosive events whose inner engine is generally expected to be a black hole or a highly magnetic neutron star (magnetar) accreting high density matter. Recognizing the nature of GRB central engines, and in particular the formation of neutron stars (NSs), is of high astrophysical significance. A possible signature of NSs in GRBs is the presence of a plateau in the early X-ray afterglow. Here we carefully select a subset of long and short GRBs with a clear plateau, and look for an additional NS signature in their prompt emission, namely a transition between accretion and propeller in analogy with accreting, magnetic compact objects in other astrophysical sources. We estimate from the prompt emission the minimum accretion luminosity below which the propeller mechanism sets in, and the NS magnetic field and spin period from the plateau. We demonstrate that these three quantities obey the same universal relation in GRBs as in other accreting compact objects switching from accretion to propeller. This relation provides also an estimate of the radiative efficiency of GRBs, which we find to be several times lower than radiatively efficient accretion in X-ray binaries and in agreement with theoretical expectations. These results provide additional support to the idea that at least some GRBs are powered by magnetars surrounded by an accretion disc. △ Less

Submitted 28 April, 2023; originally announced May 2023.

Comments: 15 pages, 5 figures, accepted for publication in The Astrophysical Journal Letters

Magnetic winding and turbulence in ultra-hot Jupiters

Authors: Clàudia Soriano-Guerrero, Daniele Viganò, Rosalba Perna, Taner Akgün, Carlos Palenzuela

Abstract: While magnetism in exoplanets remains largely unknown, Hot Jupiters have been considered as natural candidates to harbour intense magnetic fields, both due to their large masses and their high energy budgets coming from irradiation as a consequence of their vicinity to their host stars. In this work we perform MHD simulations of a narrow day-side atmospheric column of ultra-hot Jupiters, suitable… ▽ More While magnetism in exoplanets remains largely unknown, Hot Jupiters have been considered as natural candidates to harbour intense magnetic fields, both due to their large masses and their high energy budgets coming from irradiation as a consequence of their vicinity to their host stars. In this work we perform MHD simulations of a narrow day-side atmospheric column of ultra-hot Jupiters, suitable for very high local temperatures (T > 3000 K). Since the conductivity in this regime is very high, the dominant effect is winding due to the intense zonal winds. By including a forcing that mimics the wind profiles obtained in global circulation models, the shear layer induces a strong toroidal magnetic field (locally reaching hundreds of gauss), supported by meridional currents. Such fields and the sustaining currents don$'$t depend on the internally generated field, but are all confined in the thin (less than a scale-height) shear layer around 1 bar. Additionally, we add random perturbations that induce turbulent motions, which lead to further (but much smaller) magnetic field generation to a broader range of depths. These results allow an evaluation of the currents induced by the atmospheric dynamo. Although here we use ideal MHD and the only resistivity comes from the numerical scheme, we estimate a-posteriori the amount of Ohmic heat deposited in the outer layers, which could be employed in evolutionary models for Hot Jupiters' inflated radii. △ Less

Submitted 27 July, 2023; v1 submitted 14 April, 2023; originally announced April 2023.

Comments: Accepted for publication in MNRAS

Close encounters of black hole - star binaries with stellar-mass black holes

Authors: Taeho Ryu, Ruggero Valli, Rudiger Pakmor, Rosalba Perna, Selma E. de Mink, Volker Springel

Abstract: Dynamical interactions involving binaries play a crucial role in the evolution of star clusters and galaxies. We continue our investigation of the hydrodynamics of three-body encounters, focusing on binary black hole (BBH) formation, stellar disruption, and electromagnetic (EM) emission in dynamical interactions between a BH-star binary and a stellar-mass BH, using the moving-mesh hydrodynamics co… ▽ More Dynamical interactions involving binaries play a crucial role in the evolution of star clusters and galaxies. We continue our investigation of the hydrodynamics of three-body encounters, focusing on binary black hole (BBH) formation, stellar disruption, and electromagnetic (EM) emission in dynamical interactions between a BH-star binary and a stellar-mass BH, using the moving-mesh hydrodynamics code {\small AREPO}. This type of encounters can be divided into two classes depending on whether the final outcome includes BBHs. This outcome is primarily determined by which two objects meet at the first closest approach. BBHs are more likely to form when the star and the incoming BH encounter first with an impact parameter smaller than the binary's semimajor axis. In this case, the star is frequently disrupted. On the other hand, when the two BHs encounter first, frequent consequences are an orbit perturbation of the original binary or a binary member exchange. For the parameters chosen in this study, BBH formation, accompanied by stellar disruption, happens in roughly 1 out of 4 encounters. The close correlation between BBH formation and stellar disruption has possible implications for EM counterparts at the binary's merger. The BH that disrupts the star is promptly surrounded by an optically and geometrically thick disk with accretion rates exceeding the Eddington limit. If the debris disk cools fast enough to become long-lived, EM counterparts can be produced at the time of the BBH merger. △ Less

Submitted 10 July, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Comments: 16 pages, 11 figures, 2 table. Accepted for publication in MNRAS, movie here: https://studio.youtube.com/playlist/PLxLK3qI02cQfXhvhRZp9gkvts_E46t0kH/videos

GRB 191019A: a short gamma-ray burst in disguise from the disk of an active galactic nucleus

Authors: Davide Lazzati, Rosalba Perna, Benjamin Gompertz, Andrew Levan

Abstract: Long and short gamma-ray bursts (GRBs), canonically separated at around 2 seconds duration, are associated with different progenitors: the collapse of a massive star and the merger of two compact objects, respectively. GRB 191019A was a long GRB ($T_{90}\sim64$ s). Despite the relatively small redshift z=0.248 and HST followup observations, an accompanying supernova was not detected. In addition,… ▽ More Long and short gamma-ray bursts (GRBs), canonically separated at around 2 seconds duration, are associated with different progenitors: the collapse of a massive star and the merger of two compact objects, respectively. GRB 191019A was a long GRB ($T_{90}\sim64$ s). Despite the relatively small redshift z=0.248 and HST followup observations, an accompanying supernova was not detected. In addition, the host galaxy did not have significant star formation activity. Here we propose that GRB 191019A was produced by a binary compact merger, whose prompt emission was stretched in time by the interaction with a dense external medium. This would be expected if the burst progenitor was located in the disk of an active galactic nucleus, as supported by the burst localization close to the center of its host galaxy. We show that the light curve of GRB 191019A can be well modeled by a burst of intrinsic duration t=1.1 s and of energy $E_{\rm{iso}}=10^{51}$ erg seen moderately off-axis, exploding in a medium of density $10^7-10^8$ cm$^{-3}$. The double-peaked light curve carries the telltale features predicted for GRBs in high-density media, where the first peak is produced by the photosphere, and the second by the overlap of reverse shocks that take place before the internal shocks could happen. This would make GRB 191019A the first confirmed stellar explosion from within an accretion disk, with important implications for the formation and evolution of stars in accretion flows and for gravitational waves source populations. △ Less

Submitted 26 April, 2023; v1 submitted 22 March, 2023; originally announced March 2023.

Comments: 6 pages, 2 figures. Accepted for publication in ApJL

Tidal Peeling Events: low-eccentricity tidal disruption of a star by a stellar-mass black hole

Authors: Chengcheng Xin, Zoltan Haiman, Rosalba Perna, Yihan Wang, Taeho Ryu

Abstract: Close encounters between stellar-mass black holes (BHs) and stars occur frequently in dense star clusters and in the disks of active galactic nuclei (AGNs). Recent studies have shown that in highly eccentric close encounters, the star can be tidally disrupted by the BH (micro-tidal disruption event, or micro-TDE), resulting in rapid mass accretion and possibly bright electromagnetic signatures. He… ▽ More Close encounters between stellar-mass black holes (BHs) and stars occur frequently in dense star clusters and in the disks of active galactic nuclei (AGNs). Recent studies have shown that in highly eccentric close encounters, the star can be tidally disrupted by the BH (micro-tidal disruption event, or micro-TDE), resulting in rapid mass accretion and possibly bright electromagnetic signatures. Here we consider a scenario in which the star might approach the stellar-mass BH in a gradual, nearly circular inspiral, under the influence of dynamical friction on a circum-binary gas disk or three-body interactions in a star cluster. We perform hydro-dynamical simulations of this scenario using the smoothed particle hydrodynamics code PHANTOM. We find that the mass of the star is slowly stripped away by the BH. We call this gradual tidal disruption a "tidal-peeling event", or a TPE. Depending on the initial distance and eccentricity of the encounter, TPEs might exhibit significant accretion rates and orbital evolution distinct from those of a typical (eccentric) micro-TDE. △ Less

Submitted 22 March, 2023; originally announced March 2023.

On the Tidal Capture of White Dwarfs by Intermediate-mass Black Holes in Dense Stellar Environments

Authors: Claire S. Ye, Giacomo Fragione, Rosalba Perna

Abstract: Intermediate-mass black holes (IMBHs) are the missing link between stellar-mass and supermassive black holes, widely believed to reside in at least some dense star clusters, but not yet observed directly. Tidal disruptions of white dwarfs (WDs) are luminous only for black holes less massive than $\sim 10^5\,M_{\odot}$, therefore providing a unique smoking gun that could finally prove the existence… ▽ More Intermediate-mass black holes (IMBHs) are the missing link between stellar-mass and supermassive black holes, widely believed to reside in at least some dense star clusters, but not yet observed directly. Tidal disruptions of white dwarfs (WDs) are luminous only for black holes less massive than $\sim 10^5\,M_{\odot}$, therefore providing a unique smoking gun that could finally prove the existence of IMBHs beyond any reasonable doubt. Here, we investigate the tidal captures of WDs by IMBHs in dense star clusters, and estimate a typical rate of $\sim 1\,{\rm Myr}^{-1}$ for galactic nuclei and $\sim 0.01\,{\rm Myr}^{-1}$ for globular clusters. Following the capture, the WD inspirals onto the IMBH producing gravitational waves detectable out to $\sim100$ Mpc by LISA for $\sim 10^4\,M_{\odot}$ IMBHs. The subsequent tidal stripping/disruption of the WD can also release bright X-ray and gamma-ray emission with luminosities of at least $\gtrsim10^{40}\,\rm{erg\,s^{-1}}$, detectable by \textit{Chandra}, \textit{Swift}, and upcoming telescopes, such as the \textit{Einstein Probe}. △ Less

Submitted 12 August, 2023; v1 submitted 13 March, 2023; originally announced March 2023.

Comments: 13 pages, 8 figures. Published at ApJ

Effective two-body scatterings around a massive object

Authors: Yihan Wang, Saavik Ford, Rosalba Perna, Barry McKernan, Zhaohuan Zhu, Bing Zhang

Abstract: Two-body scatterings under the potential of a massive object are very common in astrophysics. If the massive body is far enough away that the two small bodies are in their own gravitational sphere of influence, the gravity of the massive body can be temporarily ignored. However, this requires the scattering process to be fast enough that the small objects do not spend too much time at distances ne… ▽ More Two-body scatterings under the potential of a massive object are very common in astrophysics. If the massive body is far enough away that the two small bodies are in their own gravitational sphere of influence, the gravity of the massive body can be temporarily ignored. However, this requires the scattering process to be fast enough that the small objects do not spend too much time at distances near the surface of the sphere of influence. In this paper, we derive the validation criteria for effective two-body scattering and establish a simple analytical solution for this process, which we verify through numerical scattering experiments. We use this solution to study star-black hole scatterings in the disks of Active Galactic Nuclei and planet-planet scatterings in planetary systems, and calculate their one-dimensional cross-section analytically. Our solution will be valuable in reducing computational time when treating two-body scatterings under the potential of a much more massive third body, provided that the problem settings are in the valid parameter space region identified by our study. △ Less

Submitted 8 March, 2023; originally announced March 2023.

Comments: Comments are Welcome

Observable signatures of stellar-mass black holes in active galactic nuclei

Authors: Hiromichi Tagawa, Shigeo S. Kimura, Zoltán Haiman, Rosalba Perna, Imre Bartos

Abstract: Stellar-mass black holes (BHs) are predicted to be embedded in the disks of active galactic nuclei (AGN) due to gravitational drag and in-situ star formation. However, clear evidence for AGN disk-embedded BHs is currently lacking. Here, as possible electromagnetic signatures of these BHs, we investigate breakout emission from shocks emerging around Blandford-Znajek jets launched from accreting BHs… ▽ More Stellar-mass black holes (BHs) are predicted to be embedded in the disks of active galactic nuclei (AGN) due to gravitational drag and in-situ star formation. However, clear evidence for AGN disk-embedded BHs is currently lacking. Here, as possible electromagnetic signatures of these BHs, we investigate breakout emission from shocks emerging around Blandford-Znajek jets launched from accreting BHs in AGN disks. We assume that the majority of the highly super-Eddington flow reaches the BH, produces a strong jet, and the jet produces feedback that shuts off accretion and thus leads to episodic flaring. While these assumptions are highly uncertain at present, they predict a breakout emission characterized by luminous thermal emission in the X-ray bands, and bright, broadband non-thermal emission from the infrared to the gamma-ray bands. The flare duration depends on the BH's distance $r$ from the central supermassive BH, varying between $10^3-10^6$ s for $r \sim 0.01-1$ pc. This emission can be discovered by current and future infrared, optical, and X-ray wide-field surveys and monitoring campaigns of nearby AGNs. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Comments: 10 pages, 4 figures, accepted in ApJL

Observable signature of merging stellar-mass black holes in active galactic nuclei

Authors: Hiromichi Tagawa, Shigeo S. Kimura, Zoltán Haiman, Rosalba Perna, Imre Bartos

Abstract: The origin of stellar-mass black hole mergers discovered through gravitational waves is being widely debated. Mergers in the disks of active galactic nuclei (AGN) represent a promising source of origin, with possible observational clues in the gravitational wave data. Beyond gravitational waves, a unique signature of AGN-assisted mergers is electromagnetic emission from the accreting black holes.… ▽ More The origin of stellar-mass black hole mergers discovered through gravitational waves is being widely debated. Mergers in the disks of active galactic nuclei (AGN) represent a promising source of origin, with possible observational clues in the gravitational wave data. Beyond gravitational waves, a unique signature of AGN-assisted mergers is electromagnetic emission from the accreting black holes. Here we show that jets launched by accreting black holes merging in an AGN disk can be detected as peculiar transients by infrared, optical, and X-ray observatories We further show that this emission mechanism can explain the possible associations between gravitational wave events and the optical transient ZTF19abanrhr and the proposed gamma-ray counterparts GW150914-GBM and LVT151012-GBM. We demonstrate how these associations, if genuine, can be used to reconstruct the properties of these events' environments. Searching for infrared and X-ray counterparts to similar electromagnetic transients in the future, once host galaxies are localized by optical observations, could provide a smoking gun signature of the mergers' AGN origin. △ Less

Submitted 18 March, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: 21 pages, 6 figures, accepted in ApJ

Deep X-ray and radio observations of the first outburst of the young magnetar Swift J1818.0-1607

Authors: A. Y. Ibrahim, A. Borghese, N. Rea, F. Coti Zelati, E. Parent, T. D. Russell, S. Ascenzi, R. Sathyaprakash, D. Gotz, S. Mereghetti, M. Topinka, M. Rigoselli, V. Savchenko, S. Campana, G. L. Israel, A. Tiengo, R. Perna, R. Turolla, S. Zane, P. Esposito, G. A. Rodrıguez Castillo, V. Graber, A. Possenti, C. Dehman, M. Ronchi , et al. (1 additional authors not shown)

Abstract: Swift J1818.0-1607 is a radio-loud magnetar with a spin period of 1.36 s and a dipolar magnetic field strength of B~3E14 G, which is very young compared to the Galactic pulsar population. We report here on the long-term X-ray monitoring campaign of this young magnetar using XMM-Newton, NuSTAR, and Swift from the activation of its first outburst in March 2020 until October 2021, as well as INTEGRAL… ▽ More Swift J1818.0-1607 is a radio-loud magnetar with a spin period of 1.36 s and a dipolar magnetic field strength of B~3E14 G, which is very young compared to the Galactic pulsar population. We report here on the long-term X-ray monitoring campaign of this young magnetar using XMM-Newton, NuSTAR, and Swift from the activation of its first outburst in March 2020 until October 2021, as well as INTEGRAL upper limits on its hard X-ray emission. The 1-10 keV magnetar spectrum is well modeled by an absorbed blackbody with a temperature of kT_BB~1.1 keV, and apparent reduction in the radius of the emitting region from ~0.6 to ~0.2 km. We also confirm the bright diffuse X-ray emission around the source extending between ~50'' and ~110''. A timing analysis revealed large torque variability, with an average spin-down rate nudot~-2.3E-11 Hz^2 that appears to decrease in magnitude over time. We also observed Swift J1818.0-1607 with the Karl G. Jansky Very Large Array (VLA) on 2021 March 22. We detected the radio counterpart to Swift J1818.0-1607 measuring a flux density of S_v = 4.38+/-0.05 mJy at 3 GHz, and a half ring-like structure of bright diffuse radio emission located at ~90'' to the west of the magnetar. We tentatively suggest that the diffuse X-ray emission is due to a dust scattering halo and that the radio structure may be associated with the supernova remnant of this young pulsar, based on its morphology. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Comments: 19 pages, 8 figures, accepted for publication on ApJ

Close Encounters of Tight Binary Stars with Stellar-mass Black Holes

Authors: Taeho Ryu, Rosalba Perna, Ruediger Pakmor, Jing-Ze Ma, Rob Farmer, Selma E. de Mink

Abstract: Strong dynamical interactions among stars and compact objects are expected in a variety of astrophysical settings, such as star clusters and the disks of active galactic nuclei. Via a suite of 3D hydrodynamics simulations using the moving-mesh code AREPO, we investigate the formation of transient phenomena and their properties in close encounters between an $2M_{\odot}$ or $20M_{\odot}$ equal-mass… ▽ More Strong dynamical interactions among stars and compact objects are expected in a variety of astrophysical settings, such as star clusters and the disks of active galactic nuclei. Via a suite of 3D hydrodynamics simulations using the moving-mesh code AREPO, we investigate the formation of transient phenomena and their properties in close encounters between an $2M_{\odot}$ or $20M_{\odot}$ equal-mass circular binary star and single $20M_{\odot}$ black hole (BH). Stars can be disrupted by the BH during dynamical interactions, naturally producing electromagnetic transient phenomena. Encounters with impact parameters smaller than the semimajor axis of the initial binary frequently lead to a variety of transients whose electromagnetic signatures are qualitatively different from those of ordinary disruption events involving just two bodies. These include the simultaneous or successive disruptions of both stars and one full disruption of one star accompanied by successive partial disruptions of the other star. On the other hand, when the impact parameter is larger than the semimajor axis of the initial binary, the binary is either simply tidally perturbed or dissociated into bound and unbound single stars ("micro-Hills" mechanism). The dissociation of $20M_{\odot}$ binaries can produce a runaway star and an active BH moving away from one another. Also, the binary dissociation can either produce an interacting binary with the BH, or a non-interacting, hard binary; both could be candidates of BH high- and low-mass X-ray binaries. Hence our simulations especially confirm that strong encounters can lead to the formation of the (generally difficult to form) BH low-mass X-ray binaries. △ Less

Submitted 5 January, 2023; v1 submitted 4 November, 2022; originally announced November 2022.

Comments: 13 pages, 8 figures, 2 tables. Accepted for publication in MNRAS. Movies: https://www.youtube.com/watch?v=DGz2WokfZLA&list=PLxLK3qI02cQeEc3ktlYAOLY4RHVX1V_nZ

The effects of Time-Variable Absorption due to Gamma-Ray Bursts In Active Galactic Nuclei Accretion Disks

Authors: Michael Ray, Davide Lazzati, Rosalba Perna

Abstract: Both long and short gamma-ray bursts (GRBs) are expected to occur in the dense environments of active galactic nuclei (AGN) accretion disks. As these bursts propagate through the disks they live in, they photoionize the medium causing time-dependent opacity that results in transients with unique spectral evolution. In this paper we use a line-of-sight radiation transfer code coupling metal and dus… ▽ More Both long and short gamma-ray bursts (GRBs) are expected to occur in the dense environments of active galactic nuclei (AGN) accretion disks. As these bursts propagate through the disks they live in, they photoionize the medium causing time-dependent opacity that results in transients with unique spectral evolution. In this paper we use a line-of-sight radiation transfer code coupling metal and dust evolution to simulate the time-dependent absorption that occurs in the case of both long and short GRBs. Through these simulations, we investigate the parameter space in which dense environments leave a potentially observable imprint on the bursts. Our numerical investigation reveals that time dependent spectral evolution is expected for central supermassive black hole masses between $10^5$ and $5 \times 10^7$ solar masses in the case of long GRBs, and between $10^4$ and $10^7$ solar masses in the case of short GRBs. Our findings can lead to the identification of bursts exploding in AGN disk environments through their unique spectral evolution coupled with a central location. In addition, the study of the time-dependent evolution would allow for studying the disk structure, once the identification with an AGN has been established. Finally, our findings lead to insight into whether GRBs contribute to the AGN emission, and which kind, thus helping to answer the question of whether GRBs can be the cause of some of the as-of-yet unexplained AGN time variability. △ Less

Submitted 17 March, 2023; v1 submitted 31 October, 2022; originally announced October 2022.

Disks of Stars in the Galactic center triggered by Tidal Disruption Events

Authors: Rosalba Perna, Evgeni Grishin

Abstract: In addition to a supermassive black hole (SMBH), the central parsec of the Milky Way hosts over a hundred of massive, high velocity young stars whose existence, and organisation of a subset of them in one, or possibly two, mis-aligned disks, is puzzling. Due to a combination of low medium density and strong tidal forces in the vicinity of Sgr A*, stars are not expected to form. Here we propose a n… ▽ More In addition to a supermassive black hole (SMBH), the central parsec of the Milky Way hosts over a hundred of massive, high velocity young stars whose existence, and organisation of a subset of them in one, or possibly two, mis-aligned disks, is puzzling. Due to a combination of low medium density and strong tidal forces in the vicinity of Sgr A*, stars are not expected to form. Here we propose a novel scenario for their in-situ formation: a jetted tidal disruption event (TDE) from an older wandering star triggers an episode of positive feedback of star formation in the plane perpendicular to the jet, as demonstrated via numerical simulations in the context of jet-induced feedback in galactic outflows. An over-pressured cocoon surrounding the jet shock-compresses clumps to densities high enough to resist the SMBH tidal field. The TDE rate of $10^{-5}-10^{-4}$ yr$^{-1}$ per galaxy, out of which a few percent events are jetted, implies a jetted TDE event per galaxy to occur every few million years. This timescale is interestingly of the same order of the age of the disk stars. The mass function predicted by our mechanism is top-heavy. Additionally, since TDEs are isotropic, our model predicts a random orientation for the disk of stars with respect to the plane of the galaxy and, due to the relatively high TDE rate, it can account for multiple disks of stars with uncorrelated orientations. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: ApJL in press

Prompt Emission of Gamma-Ray Bursts in the High-density Environment of Active Galactic Nuclei Accretion Disks

Authors: Davide Lazzati, Gustavo Soares, Rosalba Perna

Abstract: Long and short gamma-ray bursts are traditionally associated with galactic environments, where circumburst densities are small or moderate (few to hundreds of protons per cubic cm). However, both are also expected to occur in the disks of Active Galactic Nuclei, where the ambient medium density can be much larger. In this work we study, via semi-analytical methods, the propagation of the GRB outfl… ▽ More Long and short gamma-ray bursts are traditionally associated with galactic environments, where circumburst densities are small or moderate (few to hundreds of protons per cubic cm). However, both are also expected to occur in the disks of Active Galactic Nuclei, where the ambient medium density can be much larger. In this work we study, via semi-analytical methods, the propagation of the GRB outflow, its interaction with the external material, and the ensuing prompt radiation. In particular, we focus on the case in which the external shock develops early in the evolution, at a radius that is smaller than the internal shock one. We find that bursts in such high density environments are likely characterized by a single, long emission episode that is due to the superposition of individual pulses, with a characteristic hard to soft evolution irrespective of the light curve luminosity. While multi-pulse light curves are not impossible, they would require the central engine to go dormant for a long time before re-igniting. In addition, short GRB engines would produce bursts with prompt duration that would exceed the canonical 2 s separation threshold and would likely be incorrectly classified as long events, even though they would not be accompanied by a simultaneous supernova. Finally, these events have a large dynamical efficiency which would produce a bright prompt emission followed by a somewhat dim afterglow. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 9 pages, 5 figures, under review for AAS Journals

The emergence of diffused Gamma-Ray Burst afterglows from the disks of Active Galactic Nuclei

Authors: Yihan Wang, Davide Lazzati, Rosalba Perna

Abstract: The disks of Active Galactic Nuclei (AGNs) have emerged as rich environments for the production and capture of stars and the compact objects that they leave behind. These stars produce long Gamma-Ray Bursts (LGRBs) at their deaths, while frequent interactions among compact objects form binary neutron stars and neutron star-black hole binaries, leading to short GRBs (SGRBs) upon their merger. Predi… ▽ More The disks of Active Galactic Nuclei (AGNs) have emerged as rich environments for the production and capture of stars and the compact objects that they leave behind. These stars produce long Gamma-Ray Bursts (LGRBs) at their deaths, while frequent interactions among compact objects form binary neutron stars and neutron star-black hole binaries, leading to short GRBs (SGRBs) upon their merger. Predicting the properties of these transients as they emerge from the dense environments of AGN disks is key to their proper identification and to better constrain the star and compact object population in AGN disks. Some of these transients would appear unusual because they take place in much higher densities than the interstellar medium. Others, which are the subject of this paper, would additionally be modified by radiation diffusion since they are generated within optically thick regions of the accretion disks. Here we compute the GRB afterglow light curves for diffused GRB sources for a representative variety of central black-hole masses and disk locations. We find that the radiation from radio to UV and soft X-rays can be strongly suppressed by synchrotron self-absorption in the dense medium of the AGN disk. In addition, photon diffusion can significantly delay the emergence of the emission peak, turning a beamed, fast transient into a slow, isotropic, and dimmer one. These would appear as broadband-correlated AGN variability with dominance at the higher frequencies. Their properties can constrain both the stellar populations within AGN disks as well as the disk structure. △ Less

Submitted 11 July, 2022; originally announced July 2022.

Comments: Accepted for publication in MNRAS

Close Encounters of Stars with Stellar-mass Black Hole Binaries

Authors: Taeho Ryu, Rosalba Perna, Yihan Wang

Abstract: Many astrophysical environments, from star clusters and globular clusters to the disks of Active Galactic Nuclei, are characterized by frequent interactions between stars and the compact objects that they leave behind. Here, using a suite of $3-D$ hydrodynamics simulations, we explore the outcome of close interactions between $1M_{\odot}$ stars and binary black holes (BBHs) in the gravitational wa… ▽ More Many astrophysical environments, from star clusters and globular clusters to the disks of Active Galactic Nuclei, are characterized by frequent interactions between stars and the compact objects that they leave behind. Here, using a suite of $3-D$ hydrodynamics simulations, we explore the outcome of close interactions between $1M_{\odot}$ stars and binary black holes (BBHs) in the gravitational wave regime, resulting in a tidal disruption event (TDE) or a pure scattering, focusing on the accretion rates, the back reaction on the BH binary orbital parameters and the increase in the binary BH effective spin. We find that TDEs can make a significant impact on the binary orbit, which is often different from that of pure scattering. Binaries experiencing a prograde (retrograde) TDE tend to be widened (hardened) by up to $\simeq 20\%$. Initially circular binaries become more eccentric by $\lesssim 10\%$ by a prograde or retrograde TDE, whereas the eccentricity of initially eccentric binaries increases (decreases) by a retrograde (prograde) TDE by $\lesssim 5\%$. Overall a single TDE can generally result in changes of the gravitational wave-driven merger time scale by order unity. The accretion rates of both black holes are very highly super-Eddington, showing modulations (preferentially for retrograde TDEs) on a time scale of the orbital period, which can be a characteristic feature of BBH-driven TDEs. Prograde TDEs result in the effective spin parameter $χ$ to vary by $\lesssim 0.02$ while $χ\gtrsim -0.005$ for retrograde TDEs. △ Less

Submitted 2 August, 2022; v1 submitted 1 June, 2022; originally announced June 2022.

Comments: 14 pages, 12 figures, 2 tables. Accepted for publication in MNRAS

Searching for Ultra-light Bosons and Constraining Black Hole Spin Distributions with Stellar Tidal Disruption Events

Authors: Peizhi Du, Daniel Egana-Ugrinovic, Rouven Essig, Giacomo Fragione, Rosalba Perna

Abstract: Stars that pass close to the supermassive black holes located in the center of galaxies can be violently disrupted by tidal forces, leading to flares that are observed as bright transient events in sky surveys. The rate for these events to occur depends on the black hole spins, which in turn can be affected by ultra-light bosons due to superradiance. We perform a detailed analysis of these effects… ▽ More Stars that pass close to the supermassive black holes located in the center of galaxies can be violently disrupted by tidal forces, leading to flares that are observed as bright transient events in sky surveys. The rate for these events to occur depends on the black hole spins, which in turn can be affected by ultra-light bosons due to superradiance. We perform a detailed analysis of these effects and show that searches for stellar tidal disruptions have a significant potential to uncover the existence of ultra-light bosons. In particular, we find that upcoming stellar tidal disruption rate measurements by the Vera Rubin Observatory's Legacy Survey of Space and Time can be used to either discover or rule out bosons with masses ranging from $10^{-20}$ to $10^{-18}$ eV. Our analysis also indicates that these measurements may be used to constrain a variety of supermassive black hole spin distributions and determine if close-to maximal spins are preferred. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: 7 pages, 5 figures

Measuring the Non-Axially-Symmetric Surface Temperature Distribution of the Central Compact Object in Puppis A

Authors: J. A. J. Alford, E. V. Gotthelf, R. Perna, J. P. Halpern

Abstract: The surface temperature distributions of central compact objects (CCOs) are powerful probes of their crustal magnetic field strengths and geometries. Here we model the surface temperature distribution of RX J0822$-$4300, the CCO in the Puppis A supernova remnant (SNR), using $471$ ks of XMM-Newton data. We compute the energy-dependent pulse profiles in sixteen energy bands, fully including the gen… ▽ More The surface temperature distributions of central compact objects (CCOs) are powerful probes of their crustal magnetic field strengths and geometries. Here we model the surface temperature distribution of RX J0822$-$4300, the CCO in the Puppis A supernova remnant (SNR), using $471$ ks of XMM-Newton data. We compute the energy-dependent pulse profiles in sixteen energy bands, fully including the general relativistic effects of gravitational redshift and light bending, to accurately model the two heated surface regions of different temperatures and areas, in addition to constraining the viewing geometry. This results in precise measurements of the two temperatures: $kT_{\rm warm} = (1+z) \times 0.222_{-0.019}^{+0.018}$ keV and $kT_{\rm hot} = (1+z) \times 0.411\pm0.011$ keV. For the first time, we are able to measure a deviation from a pure antipodal hot-spot geometry, with a minimum value of $1.\!^{\circ}1 \pm 0.\!^{\circ}2$, and an expectation value of $9.\!^{\circ}35 \pm 0.\!^{\circ}17$ among the most probable geometries. The discovery of this asymmetry, along with the factor of $\approx2$ temperature difference between the two emitting regions, may indicate that RX J0822$-$4300 was born with a strong, tangled crustal magnetic field. △ Less

Submitted 22 December, 2021; originally announced December 2021.

Comments: 19 pages, 10 figures, Submitted to ApJ

Host galaxies and electromagnetic counterparts to binary neutron star mergers across the cosmic time: Detectability of GW170817-like events

Authors: Rosalba Perna, M. Celeste Artale, Yi-Han Wang, Michela Mapelli, Davide Lazzati, Cecilia Sgalletta, Filippo Santoliquido

Abstract: The detection of electromagnetic radiation (EM) accompanying the gravitational wave (GW) signal from the binary neutron star (BNS) merger GW170817 has revealed that these systems constitute at least a fraction of the progenitors of short gamma-ray bursts (SGRBs). As gravitational wave detectors keep pushing their detection horizons, it is important to assess coupled GW/EM probabilities, and how to… ▽ More The detection of electromagnetic radiation (EM) accompanying the gravitational wave (GW) signal from the binary neutron star (BNS) merger GW170817 has revealed that these systems constitute at least a fraction of the progenitors of short gamma-ray bursts (SGRBs). As gravitational wave detectors keep pushing their detection horizons, it is important to assess coupled GW/EM probabilities, and how to maximize observational prospects. Here we perform population synthesis calculations of BNS evolution with the code MOBSE, and seed the binaries in galaxies at three representative redshifts (z=0.01,0.1,1) of the Illustris TNG50 simulation. The binaries are evolved and their locations numerically tracked in the host galactic potentials until merger. Adopting the astrophysical parameters of GRB170817A as a prototype, we numerically compute the broadband lightcurves of jets from BNS mergers, with the afterglow brightness depending on the local medium density at the merger sites. We perform Monte Carlo simulations of the resulting EM population assuming either a random viewing angle with respect to the jet, or a jet aligned with the orbital angular momentum of the binary, which biases the viewing angle probability for GW-triggered events. We find that ~70-80% of BNSs from z=0.01 should be detectable in gamma-rays. The afterglow detection probabilities of GW-triggered BNS mergers vary between ~0.3-0.7%, with higher values for jets aligned with the BNS angular momentum, and are comparable across the high and low-energy bands, unlike gamma-ray-triggered events (cosmological SGRBs) which are significantly brighter at higher energies. We further quantify observational biases with respect to host galaxy masses. △ Less

Submitted 9 December, 2021; originally announced December 2021.

Comments: 15 pages, 12 figures, submitted to MNRAS

Can stellar-mass black hole growth disrupt disks of active galactic nuclei? The role of mechanical feedback

Authors: Hiromichi Tagawa, Shigeo S. Kimura, Zoltán Haiman, Rosalba Perna, Hidekazu Tanaka, Imre Bartos

Abstract: Stellar-mass BHs (sBHs) are predicted to be embedded in active galactic nuclei (AGN) disks due to gravitational drag and in-situ star formation. However, we find that due to a high gas density in an AGN disk environment, compact objects may rapidly grow to intermediate-mass BHs and deplete matter from the AGN disk unless accretion is suppressed by some feedback process(es). These consequences are… ▽ More Stellar-mass BHs (sBHs) are predicted to be embedded in active galactic nuclei (AGN) disks due to gravitational drag and in-situ star formation. However, we find that due to a high gas density in an AGN disk environment, compact objects may rapidly grow to intermediate-mass BHs and deplete matter from the AGN disk unless accretion is suppressed by some feedback process(es). These consequences are inconsistent with AGN observations and the dynamics of the Galactic center. Here we consider mechanical feedback mechanisms for the reduction of gas accretion. Rapidly accreting sBHs launch winds and/or jets via the Blandford-Znajek mechanism, which produce high-pressure shocks and cocoons. Such a shock and cocoon can spread laterally in the plane of the disk, eject the outer regions of a circum-sBH disk (CsBD) and puncture a hole in the AGN disk with horizontal size comparable to the disk scale-height. Since the depletion timescale of the bound CsBD is much shorter than the resupply timescale of gas to the sBH, the time-averaged accretion rate onto sBHs is reduced by this process by a factor of $\sim 10$--$100$. This feedback mechanism can therefore help alleviate the sBH over-growth and AGN-disk depletion problems. On the other hand, we find that cocoons of jets can unbind a large fraction of the gas accreting in the disks of less massive SMBHs, which may help explain the dearth of high-Eddington ratio AGNs with SMBH mass $\lesssim10^5{\rm M_\odot}$. △ Less

Submitted 2 December, 2021; originally announced December 2021.

Comments: 19 pages, 9 figures, revised version after the first referee report

Symmetry Breaking in Dynamical Encounters in the Disks of Active Galactic Nuclei

Authors: Yihan Wang, Barry McKernan, Saavik Ford, Rosalba Perna, Nathan W. C. Leigh, Mordecai-Mark Mac Low

Abstract: Active galactic nucleus (AGN) disks may be important sites of binary black hole (BBH) mergers. Here we show via numerical experiments with the high-accuracy, high precision code {\tt SpaceHub} that broken symmetry in dynamical encounters in AGN disks can lead to an asymmetry between prograde and retrograde BBH mergers. The direction of the hardening asymmetry depends on the initial binary semi-maj… ▽ More Active galactic nucleus (AGN) disks may be important sites of binary black hole (BBH) mergers. Here we show via numerical experiments with the high-accuracy, high precision code {\tt SpaceHub} that broken symmetry in dynamical encounters in AGN disks can lead to an asymmetry between prograde and retrograde BBH mergers. The direction of the hardening asymmetry depends on the initial binary semi-major axis. An asymmetric distribution of mass-weighted projected spin $χ_{\rm eff}$ should therefore be expected in LIGO-Virgo detections of BBH mergers from AGN disks. This channel further predicts that negative $χ_{\rm eff}$ BBH mergers are most likely for massive binaries. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: submitted to ApJL. Comments are welcome

Hot Jupiter formation in dense clusters: secular chaos in multi-planetary systems

Authors: Yihan Wang, Rosalba Perna, Nathan W. C. Leigh, Michael M. Shara

Abstract: Exoplanetary observations reveal that the occurrence rate of hot Jupiters is correlated with star clustering. In star clusters, interactions between planetary systems and close fly-by stars can significantly change the architecture of primordially coplanar, circular planetary systems. Flybys in dense clusters have a significant impact on hot Jupiter formation via activation of high eccentricity ex… ▽ More Exoplanetary observations reveal that the occurrence rate of hot Jupiters is correlated with star clustering. In star clusters, interactions between planetary systems and close fly-by stars can significantly change the architecture of primordially coplanar, circular planetary systems. Flybys in dense clusters have a significant impact on hot Jupiter formation via activation of high eccentricity excitation mechanisms such as the Zeipel-Lidov-Kozai (ZLK) effect and planet-planet scattering. Previous studies have shown that if there are two giant planets in the planetary system, close flybys can efficiently activate the ZLK mechanism, thus triggering high eccentricity tidal migration and ultimately form hot Jupiters in star clusters. Here we extend our previous study with a multi-planet (triple) system. We perform high precision, high-accuracy few-body simulations of stellar flybys and subsequent planetary migration within the perturbed planetary systems using the code {\tt SpaceHub}. Our simulations demonstrate that a single close flyby on a multi-planet system in a cluster can activate secular chaos and ultimately lead to hot Jupiter formation via high eccentricity migration. We find that the hot Jupiter formation rate per system increases with both the size of the planetary system as well as with the mass of the outer planet, and we quantify the relative formation fractions for a range of parameters. Hot Jupiters formed via secular chaos are expected to be accompanied by massive companions with very long periods. Our study further shows that this flyby-induced secular chaos is preferred in low-density clusters where multi-planet systems are more likely to survive, and that it contributes a significant fraction of the hot Jupiter formation in star clusters compared to the flyby-induced ZLK mechanism. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: submitted to MNRAS

Two steps forward and one step sideways: the propagation of relativistic jets in realistic binary neutron star merger ejecta

Authors: Davide Lazzati, Rosalba Perna, Riccardo Ciolfi, Bruno Giacomazzo, Diego Lopez-Camara, Brian Morsony

Abstract: The association of GRB170817A with GW170817 has confirmed the long-standing hypothesis that binary neutron star (BNS) mergers are the progenitors of at least some short gamma-ray bursts (SGRBs). This connection has ushered in an era in which broadband observations of SGRBs, together with measurements of the time delay between the gravitational waves and the electromagnetic radiation, allow to prob… ▽ More The association of GRB170817A with GW170817 has confirmed the long-standing hypothesis that binary neutron star (BNS) mergers are the progenitors of at least some short gamma-ray bursts (SGRBs). This connection has ushered in an era in which broadband observations of SGRBs, together with measurements of the time delay between the gravitational waves and the electromagnetic radiation, allow to probe the properties of the emitting outflow and its engine to an unprecedented detail. Since the structure of the radiating outflow is molded by the interaction of a relativistic jet with the binary ejecta, it is of paramount importance to study the system in a realistic setting. Here we present a three-dimensional hydrodynamic simulation of a relativistic jet propagating in the ejecta of a BNS merger, which were computed with a general relativistic magnetohydrodynamic simulation. We find that the jet's centroid oscillates around the axis of the system, due to inhomogeneities encountered in the propagation. These oscillations allow the jet to find the path of least resistance and travel faster than an identical jet in smooth ejecta. In our setup the breakout time is ~0.6 sec, comparable to the expected central engine duration in SGRBs and possibly a non-negligible fraction of the total delay between the gravitational and gamma-ray signals. Our simulation also shows that energy is carried in roughly equal amounts by the jet and by the cocoon, and that about 20 per cent of the injected energy is transferred to the ejecta via mechanical work. △ Less

Submitted 16 July, 2021; originally announced July 2021.

Comments: Under review for publication in ApJL

Micro Tidal Disruption Events in Active Galactic Nuclei

Authors: Y. Yang, I. Bartos, G. Fragione, Z. Haiman, M. Kowalski, S. Marka, R. Perna, H. Tagawa

Abstract: Active galactic nuclei (AGNs) can act as black hole assembly lines, funneling some of the stellar-mass black holes from the vicinity of the galactic center into the inner plane of the AGN disk where the black holes can merge through dynamical friction and gravitational wave emission. Here, we show that stars near the galactic center are also brought into the AGN disk, where they can be tidally dis… ▽ More Active galactic nuclei (AGNs) can act as black hole assembly lines, funneling some of the stellar-mass black holes from the vicinity of the galactic center into the inner plane of the AGN disk where the black holes can merge through dynamical friction and gravitational wave emission. Here, we show that stars near the galactic center are also brought into the AGN disk, where they can be tidally disrupted by the stellar-mass black holes in the disk. Such micro-tidal disruption events (micro-TDEs) could be useful probe of stellar interaction with the AGN disk. We find that micro-TDEs in AGNs occur at a rate of $\sim170$ Gpc$^{-3}$yr$^{-1}$. Their cleanest observational probe may be the detection of tidal disruption in AGNs by heavy supermassive black holes ($M_{\bullet}\gtrsim10^{8}$ M$_{\odot}$) so that cannot tidally disrupt solar-type stars. We discuss two such TDE candidates observed to date (ASASSN-15lh and ZTF19aailpwl). △ Less

Submitted 5 May, 2021; originally announced May 2021.

Comments: 7 pages, 1 figure

The X-ray evolution and geometry of the 2018 outburst of XTE J1810-197

Authors: A. Borghese, N. Rea, R. Turolla, M. Rigoselli, J. A. J. Alford, E. V. Gotthelf, M. Burgay, A. Possenti, S. Zane, F. Coti Zelati, R. Perna, P. Esposito, S. Mereghetti, D. Viganó, A. Tiengo, D. Götz, A. Ibrahim, G. L. Israel, J. Pons, R. Sathyaprakash

Abstract: After 15 years, in late 2018, the magnetar XTE J1810-197 underwent a second recorded X-ray outburst event and reactivated as a radio pulsar. We initiated an X-ray monitoring campaign to follow the timing and spectral evolution of the magnetar as its flux decays using Swift, XMM-Newton, NuSTAR, and NICER observations. During the year-long campaign, the magnetar reproduced similar behaviour to that… ▽ More After 15 years, in late 2018, the magnetar XTE J1810-197 underwent a second recorded X-ray outburst event and reactivated as a radio pulsar. We initiated an X-ray monitoring campaign to follow the timing and spectral evolution of the magnetar as its flux decays using Swift, XMM-Newton, NuSTAR, and NICER observations. During the year-long campaign, the magnetar reproduced similar behaviour to that found for the first outburst, with a factor of two change in its spin-down rate from $\sim7.2\times10^{-12}$ s s$^{-1}$ to $\sim1.5\times10^{-11}$ s s$^{-1}$ after two months. Unique to this outburst, we confirm the peculiar energy-dependent phase shift of the pulse profile. Following the initial outburst, the spectrum of XTE J1810-197 is well-modelled by multiple blackbody components corresponding to a pair of non-concentric, hot thermal caps surrounded by a cooler one, superposed to the colder star surface. We model the energy-dependent pulse profile to constrain the viewing and surface emission geometry and find that the overall geometry of XTE J1810-197 has likely evolved relative to that found for the 2003 event. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Comments: 14 pages, 12 figures, 4 tables. Accepted for publications in MNRAS

SpaceHub: A high-performance gravity integration toolkit for few-body problems in astrophysics

Authors: Yi-Han Wang, Nathan Leigh, Bin Liu, Rosalba Perna

Abstract: We present the open source few-body gravity integration toolkit {\tt SpaceHub}. {\tt SpaceHub} offers a variety of algorithmic methods, including the unique algorithms AR-Radau, AR-Sym6, AR-ABITS and AR-chain$^+$ which we show out-perform other methods in the literature and allow for fast, precise and accurate computations to deal with few-body problems ranging from interacting black holes to plan… ▽ More We present the open source few-body gravity integration toolkit {\tt SpaceHub}. {\tt SpaceHub} offers a variety of algorithmic methods, including the unique algorithms AR-Radau, AR-Sym6, AR-ABITS and AR-chain$^+$ which we show out-perform other methods in the literature and allow for fast, precise and accurate computations to deal with few-body problems ranging from interacting black holes to planetary dynamics. We show that AR-Sym6 and AR-chain$^+$, with algorithmic regularization, chain algorithm, active round-off error compensation and a symplectic kernel implementation, are the fastest and most accurate algorithms to treat black hole dynamics with extreme mass ratios, extreme eccentricities and very close encounters. AR-Radau, the first regularized Radau integrator with round off error control down to 64 bits floating point machine precision, has the ability to handle extremely eccentric orbits and close approaches in long-term integrations. AR-ABITS, a bit efficient arbitrary precision method, achieves any precision with the least CPU cost compared to other open source arbitrary precision few-body codes. With the implementation of deep numerical and code optimization, these new algorithms in {\tt SpaceHub} prove superior to other popular high precision few-body codes in terms of performance, accuracy and speed. △ Less

Submitted 29 April, 2021; v1 submitted 13 April, 2021; originally announced April 2021.

Comments: Submitted to MNRAS. Comments are welcome

Accretion-Induced Collapse of Neutron Stars in the Disks of Active Galactic Nuclei

Authors: Rosalba Perna, Hiromichi Tagawa, Zoltan Haiman, Imre Bartos

Abstract: The disks of active galactic nuclei (AGNs) have emerged as a rich environment for the evolution of stars and their compact remnants. The very dense medium favors rapid accretion, while torques and migration traps enhance binary formation and mergers. Both long and short gamma-ray bursts (GRBs) are hence expected. We show that AGN disks constitute an ideal environment for another interesting phenom… ▽ More The disks of active galactic nuclei (AGNs) have emerged as a rich environment for the evolution of stars and their compact remnants. The very dense medium favors rapid accretion, while torques and migration traps enhance binary formation and mergers. Both long and short gamma-ray bursts (GRBs) are hence expected. We show that AGN disks constitute an ideal environment for another interesting phenomenon: the accretion induced collapse (AIC) of neutron stars (NSs) to black holes (BHs). Rapid accretion in the dense disks can cause NSs to grow to the point of exceeding the maximum mass allowed by their equation of state. General relativistic magnetohydrodynamical simulations have shown that electromagnetic signatures are expected if the NS is surrounded by a mini-disk prior to collapse, which then rapidly accretes onto the BH, and/or if the NS is highly magnetized, from reconnection of the magnetosphere during collapse. Here we compute the rates of AICs and their locations within the disks for both isolated NSs, and for (initially stable) NSs formed from NS-NS mergers. We find that the global AIC rates are $\sim 0.07-20$~Gpc$^{-3}$~yr$^{-1}$, and we discuss their observable prospects and signatures as they emerge from the dense disk environments. △ Less

Submitted 19 March, 2021; originally announced March 2021.

Comments: Submitted to AAS Journals. Comments are welcome

Partial tidal disruption events by stellar mass black holes: gravitational instability of stream and impact from remnant core

Authors: Yi-Han Wang, Rosalba Perna, Philip J. Armitage

Abstract: In dense star clusters, such as globular and open clusters, dynamical interactions between stars and black holes (BHs) can be extremely frequent, leading to various astrophysical transients. Close encounters between a star and a stellar mass BH make it possible for the star to be tidally disrupted by the BH. Due to the relative low mass of the BH and the small cross section of the tidal disruption… ▽ More In dense star clusters, such as globular and open clusters, dynamical interactions between stars and black holes (BHs) can be extremely frequent, leading to various astrophysical transients. Close encounters between a star and a stellar mass BH make it possible for the star to be tidally disrupted by the BH. Due to the relative low mass of the BH and the small cross section of the tidal disruption event (TDE) for cases with high penetration, disruptions caused by close encounters are usually partial disruptions. The existence of the remnant stellar core and its non-negligible mass compared to the stellar mass BH alters the accretion process significantly. We study this problem with SPH simulations using the code {\tt Phantom}, with the inclusion of radiation pressure, which is important for small mass BHs. Additionally, we develop a new, more general method of computing the fallback rate which does not rely on any approximation. Our study shows that the powerlaw slope of the fallback rate has a strong dependence on the mass of the BH in the stellar mass BH regime. Furthermore, in this regime, self-gravity of the fallback stream and local instabilities become more significant, and cause the disrupted material to collapse into small clumps before returning to the BH. This results in an abrupt increase of the fallback rate, which can significantly deviate from a powerlaw. Our results will help in the identification of TDEs by stellar mass BHs in dense clusters. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: Accepted for publication in MNRAS

The variable absorption in the X-ray spectrum of GRB 190114C

Authors: Sergio Campana, Davide Lazzati, Rosalba Perna, Maria Grazia Bernardini, Lara Nava

Abstract: GRB 190114C was a bright burst that occurred in the local Universe (z=0.425). It was the first gamma-ray burst (GRB) ever detected at TeV energies, thanks to MAGIC. We characterize the ambient medium properties of the host galaxy through the study of the absorbing X-ray column density. Joining Swift, XMM-Newton, and NuSTAR observations, we find that the GRB X-ray spectrum is characterized by a hig… ▽ More GRB 190114C was a bright burst that occurred in the local Universe (z=0.425). It was the first gamma-ray burst (GRB) ever detected at TeV energies, thanks to MAGIC. We characterize the ambient medium properties of the host galaxy through the study of the absorbing X-ray column density. Joining Swift, XMM-Newton, and NuSTAR observations, we find that the GRB X-ray spectrum is characterized by a high column density that is well in excess of the expected Milky Way value and decreases, by a factor of ~2, around ~$10^5$ s. Such a variability is not common in GRBs. The most straightforward interpretation of the variability in terms of photoionization of the ambient medium is not able to account for the decrease at such late times, when the source flux is less intense. Instead, we interpret the decrease as due to a clumped absorber, denser along the line of sight and surrounded by lower-density gas. After the detection at TeV energies of GRB 190114C, two other GRBs were promptly detected. They share a high value of the intrinsic column density and there are hints for a decrease of the column density, too. We speculate that a high local column density might be a common ingredient for TeV-detected GRBs. △ Less

Submitted 9 March, 2021; originally announced March 2021.

Comments: 6 pages, 4 figures. Accepted for publication in A&A

Journal ref: A&A 649, A135 (2021)

Stellar Evolution in the Disks of Active Galactic Nuclei Produces Rapidly Rotating Massive Stars

Authors: Adam S. Jermyn, Alexander J. Dittmann, Matteo Cantiello, Rosalba Perna

Abstract: Stars can either be formed in or captured by the accretion disks in Active Galactic Nuclei (AGN). These AGN stars are irradiated and subject to extreme levels of accretion, which can turn even low-mass stars into very massive ones ($M > 100 {\rm M}_\odot$) whose evolution may result in the formation of massive compact objects ($M > 10 {\rm M}_\odot$). Here we explore the spins of these AGN stars a… ▽ More Stars can either be formed in or captured by the accretion disks in Active Galactic Nuclei (AGN). These AGN stars are irradiated and subject to extreme levels of accretion, which can turn even low-mass stars into very massive ones ($M > 100 {\rm M}_\odot$) whose evolution may result in the formation of massive compact objects ($M > 10 {\rm M}_\odot$). Here we explore the spins of these AGN stars and the remnants they leave behind. We find that AGN stars rapidly spin up via accretion, eventually reaching near-critical rotation rates. They further maintain near-critical rotation even as they shed their envelopes, become compact, and undergo late stages of burning. This makes them good candidates to produce high-spin massive black holes, such as the ones seen by LIGO-Virgo in GW190521g, as well as long Gamma Ray Bursts (GRBs) and the associated chemical pollution of the AGN disk. △ Less

Submitted 25 February, 2021; originally announced February 2021.

Comments: 19 pages, 7 figures. Submitted to ApJ, comments welcome!

Electromagnetic Signatures of Relativistic Explosions in AGN Disks

Authors: Rosalba Perna, Davide Lazzati, Matteo Cantiello

Abstract: The disks of Active Galactic Nuclei (AGNs), traditionally studied as the feeders of the supermassive black holes (SMBHs) at their centers, have recently triggered a lot of interest also as hosts to massive stars and hence their neutron star (NS) and black hole (BH) remnants. Migration traps and gas torques in these disks favor binary formation and enhance the rate of compact object mergers. In the… ▽ More The disks of Active Galactic Nuclei (AGNs), traditionally studied as the feeders of the supermassive black holes (SMBHs) at their centers, have recently triggered a lot of interest also as hosts to massive stars and hence their neutron star (NS) and black hole (BH) remnants. Migration traps and gas torques in these disks favor binary formation and enhance the rate of compact object mergers. In these environments both long gamma-ray bursts (GRBs) from the death of massive stars and short GRBs from NS-NS and NS-BH mergers are expected. However, their properties in the unique environments of AGN disks have never been studied. Here we show that GRBs in AGNs can display unique features, owing to the unusual relative position of the shocks that characterize the burst evolution and the Thomson photosphere of the AGN disk. In dense environments, for example, the external shock develops before the internal shocks, leading to prompt emission powered by a relativistic reverse shock. The transient's time evolution is also compressed, yielding afterglow emission that is much brighter and peaks much earlier than for GRBs in the interstellar medium. Additionally, in regions of the disk that are sufficiently dense and extended, the light curves are dominated by diffusion, since the fireball is trapped inside the disc photosphere. These diffusion-dominated transients emerge on timescales of days in disks around SMBHs of $\sim 10^6 M_\odot$ to years for SMBHs of $\sim 10^8 M_\odot$. Finally, a large fraction of events, especially in AGNs with SMBHs $\lesssim 10^7 M_\odot$, display time-variable absorption in the X-ray band. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: Submitted to ApJ Letters

Hot Jupiter and ultra-cold Saturn formation in dense star clusters

Authors: Yi-Han Wang, Nathan W. C. Leigh, Rosalba Perna, Michael M. Shara

Abstract: The discovery of high incidence of hot Jupiters in dense clusters challenges the field-based hot Jupiter formation theory. In dense clusters, interactions between planetary systems and flyby stars are relatively common. This has a significant impact on planetary systems, dominating hot Jupiter formation. In this paper, we perform high precision, few-body simulations of stellar flybys and subsequen… ▽ More The discovery of high incidence of hot Jupiters in dense clusters challenges the field-based hot Jupiter formation theory. In dense clusters, interactions between planetary systems and flyby stars are relatively common. This has a significant impact on planetary systems, dominating hot Jupiter formation. In this paper, we perform high precision, few-body simulations of stellar flybys and subsequent planet migration in clusters. A large parameter space exploration demonstrates that close flybys that change the architecture of the planetary system can activate high eccentricity migration mechanisms: Lidov-Kozai and planet-planet scattering, leading to high hot Jupiter formation rate in dense clusters. Our simulations predict that many of the hot Jupiters are accompanied by "ultra-cold Saturns", expelled to apastra of thousands of AU. This increase is particularly remarkable for planetary systems originally hosting two giant planets with semi-major axis ratios $\sim$ 4 and the flyby star approaching nearly perpendicular to the planetary orbital plane. The estimated lower limit to the hot Jupiter formation rate of a virialized cluster is $\sim 1.6\times10^{-4}(σ/{\rm 1kms^{-1}})^5({a_{\rm p}}/{\rm 20 AU})({M_{\rm c}}/{\rm 1000M_\odot})^{-2}$Gyr$^{-1}$ per star, where $σ$ is the cluster velocity dispersion, $a_{\rm p}$ is the size of the planetary system and $M_{\rm c}$ is the mass of the cluster. Our simulations yield a hot Jupiter abundance which is $\sim$ 50 times smaller than that observed in the old open cluster M67. We expect that interactions involving binary stars, as well as a third or more giant planets, will close the discrepancy. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: Accepted for publication in ApJS

On the rate of crustal failures in young magnetars

Authors: Clara Dehman, Daniele Viganò, Nanda Rea, Jose A. Pons, Rosalba Perna, Alberto Gracía-Gracía

Abstract: The activity of magnetars is powered by their intense and dynamic magnetic fields and has been proposed as the trigger to extragalactic Fast Radio Bursts. Here we estimate the frequency of crustal failures in young magnetars, by computing the magnetic stresses in detailed magneto-thermal simulations including Hall drift and Ohmic dissipation. The initial internal topology at birth is poorly known… ▽ More The activity of magnetars is powered by their intense and dynamic magnetic fields and has been proposed as the trigger to extragalactic Fast Radio Bursts. Here we estimate the frequency of crustal failures in young magnetars, by computing the magnetic stresses in detailed magneto-thermal simulations including Hall drift and Ohmic dissipation. The initial internal topology at birth is poorly known but is likely to be much more complex than a dipole. Thus, we explore a wide range of initial configurations, finding that the expected rate of crustal failures varies by orders of magnitude depending on the initial magnetic configuration. Our results show that this rate scales with the crustal magnetic energy, rather than with the often used surface value of the dipolar component related to the spin-down torque. The estimated frequency of crustal failures for a given dipolar component can vary by orders of magnitude for different initial conditions, depending on how much magnetic energy is distributed in the crustal non-dipolar components, likely dominant in newborn magnetars. The quantitative reliability of the expected event rate could be improved by a better treatment of the magnetic evolution in the core and the elastic/plastic crustal response, here not included. Regardless of that, our results are useful inputs in modelling the outburst rate of young Galactic magnetars, and their relation with the Fast Radio Bursts in our and other galaxies. △ Less

Submitted 1 October, 2020; originally announced October 2020.

Comments: 9 pages, 4 figures, ApJ Letter in press

Calibrating the binary black hole population in nuclear star clusters through tidal disruption events

Authors: Giacomo Fragione, Rosalba Perna, Abraham Loeb

Abstract: As the sensitivity of gravitational wave (GW) instruments improves and new networks start operating, hundreds of merging stellar-mass black holes (SBHs) and intermediate-mass black holes (IMBHs) are expected to be observed in the next few years. The origin and distribution of SBH and IMBH binaries in various dynamical environments is a fundamental scientific question in GW astronomy. In this paper… ▽ More As the sensitivity of gravitational wave (GW) instruments improves and new networks start operating, hundreds of merging stellar-mass black holes (SBHs) and intermediate-mass black holes (IMBHs) are expected to be observed in the next few years. The origin and distribution of SBH and IMBH binaries in various dynamical environments is a fundamental scientific question in GW astronomy. In this paper, we discuss ways tidal disruption events (TDEs) may provide a unique electromagnetic window into the assembly and merger of binary SBHs and IMBHs in nuclear star clusters (NSCs). We discuss how the host NSC mass and density and the slope of the black-hole mass function set the orbital properties and the masses of the binaries that undergo a TDE. For typical NSC properties, we predict a TDE rate of $\sim 10^{-6}$--$10^{-7}\ {\rm yr}^{-1}$ per galaxy. The lightcurve of TDEs in NSCs could be interrupted and modulated by the companion black hole on the orbital period of the binary. These should be readily detectable by optical transient surveys such as the Zwicky Transient Facility and LSST. △ Less

Submitted 5 November, 2020; v1 submitted 25 June, 2020; originally announced June 2020.

Comments: 13 pages, 7 figures, 1 table, accepted by MNRAS

Mergers of Equal-Mass Binaries with Compact Object Companions from Mass Transfer in Triple Star Systems

Authors: N. W. C. Leigh, S. Toonen, S. F. Portegies Zwart, R. Perna

Abstract: In this paper, we consider triple systems composed of main-sequence (MS) stars, and their internal evolution due to stellar and binary evolution. Our focus is on triples that produce white dwarfs (WDs), where Roche lobe overflow of an evolving tertiary triggers accretion onto the inner binary via a circumbinary disk (CBD) driving it toward a mass ratio of unity. We present a combination of analyti… ▽ More In this paper, we consider triple systems composed of main-sequence (MS) stars, and their internal evolution due to stellar and binary evolution. Our focus is on triples that produce white dwarfs (WDs), where Roche lobe overflow of an evolving tertiary triggers accretion onto the inner binary via a circumbinary disk (CBD) driving it toward a mass ratio of unity. We present a combination of analytic- and population synthesis-based calculations performed using the \texttt{SeBa} code to constrain the expected frequency of such systems, given a realistic initial population of MS triples, and provide the predicted distributions of orbital periods. We identify the parameter space for triples that can accommodate a CBD, to inform future numerical simulations of suitable initial conditions. We find that $\lesssim$ 10\% of all MS triples should be able to accommodate a CBD around the inner binary, and compute lower limits for the production rates. This scenario broadly predicts mergers of near equal-mass binaries, producing blue stragglers (BSs), Type Ia supernovae, gamma ray bursts and gravitational wave-induced mergers, along with the presence of an outer WD tertiary companion. We compare our predicted distributions to a sample of field BS binaries, and argue that our proposed mechanism explains the observed range of orbital periods. Finally, the mechanism considered here could produce hypervelocity MS stars, WDs and even millisecond pulsars with masses close to the Chandrasekhar mass limit, and be used to constrain the maximum remnant masses at the time of any supernova explosion. △ Less

Submitted 23 June, 2020; v1 submitted 13 May, 2020; originally announced May 2020.

Comments: 15 pages, 7 figures; accepted for publication in MNRAS

Intrinsic properties of the engine and jet that powered the short gamma-ray burst associated with GW170817

Authors: Davide Lazzati, Riccardo Ciolfi, Rosalba Perna

Abstract: GRB 170817A was a subluminous short gamma-ray burst detected about 1.74 s after the gravitational wave signal GW170817 from a binary neutron star (BNS) merger. It is now understood as an off-axis event powered by the cocoon of a relativistic jet pointing 15 to 30 degrees away from the direction of observation. The cocoon was energized by the interaction of the incipient jet with the non-relativist… ▽ More GRB 170817A was a subluminous short gamma-ray burst detected about 1.74 s after the gravitational wave signal GW170817 from a binary neutron star (BNS) merger. It is now understood as an off-axis event powered by the cocoon of a relativistic jet pointing 15 to 30 degrees away from the direction of observation. The cocoon was energized by the interaction of the incipient jet with the non-relativistic baryon wind from the merger remnant, resulting in a structured outflow with a narrow core and broad wings. In this paper, we couple the observational constraints on the structured outflow with a model for the jet-wind interaction to constrain the intrinsic properties with which the jet was launched by the central engine, including its time delay from the merger event. Using wind prescriptions inspired by magnetized BNS merger simulations, we find that the jet was launched within about 0.4 s from the merger, implying that the 1.74 s observed delay was dominated by the fireball propagation up to the photospheric radius. We also constrain, for the first time for any gamma-ray burst, the jet opening angle at injection and set a lower limit to its asymptotic Lorentz factor. These findings suggest an initially Poynting-flux dominated jet, launched via electromagnetic processes. If the jet was powered by an accreting black hole, they also provide a significant constraint on the survival time of the metastable neutron star remnant. △ Less

Submitted 9 June, 2020; v1 submitted 21 April, 2020; originally announced April 2020.

Comments: Accepted for publication in ApJ

A very young radio-loud magnetar

Authors: P. Esposito, N. Rea, A. Borghese, F. Coti Zelati, D. Viganò, G. L. Israel, A. Tiengo, A. Ridolfi, A. Possenti, M. Burgay, D. Götz, F. Pintore, L. Stella, C. Dehman, M. Ronchi, S. Campana, A. Garcia-Garcia, V. Graber, S. Mereghetti, R. Perna, G. A. Rodríguez Castillo, R. Turolla, S. Zane

Abstract: The magnetar Swift ,J1818.0-1607 was discovered in March 2020 when Swift detected a 9 ms hard X-ray burst and a long-lived outburst. Prompt X-ray observations revealed a spin period of 1.36 s, soon confirmed by the discovery of radio pulsations. We report here on the analysis of the Swift burst and follow-up X-ray and radio observations. The burst average luminosity was… ▽ More The magnetar Swift ,J1818.0-1607 was discovered in March 2020 when Swift detected a 9 ms hard X-ray burst and a long-lived outburst. Prompt X-ray observations revealed a spin period of 1.36 s, soon confirmed by the discovery of radio pulsations. We report here on the analysis of the Swift burst and follow-up X-ray and radio observations. The burst average luminosity was $L_{\rm burst} \sim2\times 10^{39}$ erg/s (at 4.8 kpc). Simultaneous observations with XMM-Newton and NuSTAR three days after the burst provided a source spectrum well fit by an absorbed blackbody ($N_{\rm H} = (1.13\pm0.03) \times 10^{23}$ cm$^{-2}$ and $kT = 1.16\pm0.03$ keV) plus a power-law ($Γ=0.0\pm1.3$) in the 1-20 keV band, with a luminosity of $\sim$$8\times10^{34}$ erg/s, dominated by the blackbody emission. From our timing analysis, we derive a dipolar magnetic field $B \sim 7\times10^{14}$ G, spin-down luminosity $\dot{E}_{\rm rot} \sim 1.4\times10^{36}$ erg/s and characteristic age of 240 yr, the shortest currently known. Archival observations led to an upper limit on the quiescent luminosity $<$$5.5\times10^{33}$ erg/s, lower than the value expected from magnetar cooling models at the source characteristic age. A 1 hr radio observation with the Sardinia Radio Telescope taken about 1 week after the X-ray burst detected a number of strong and short radio pulses at 1.5 GHz, in addition to regular pulsed emission; they were emitted at an average rate 0.9 min$^{-1}$ and accounted for $\sim$50% of the total pulsed radio fluence. We conclude that Swift ,J1818.0-1607 is a peculiar magnetar belonging to the small, diverse group of young neutron stars with properties straddling those of rotationally and magnetically powered pulsars. Future observations will make a better estimation of the age possible by measuring the spin-down rate in quiescence. △ Less

Submitted 25 May, 2020; v1 submitted 8 April, 2020; originally announced April 2020.

Comments: 10 pages, 4 figures, 1 table; revised version accepted for publication in The Astrophysical Journal Letters

Journal ref: The Astrophysical Journal Letters, Volume 896, L30 (2020)

The X-ray outburst of the Galactic Center magnetar over six years of Chandra observations

Authors: N. Rea, F. Coti Zelati, D. Vigano', A. Papitto, F. Baganoff, A. Borghese, S. Campana, P. Esposito, D. Haggard, G. L. Israel, S. Mereghetti, R. Mignani, R. Perna, J. A. Pons, G. Ponti, L. Stella, D. F. Torres, R. Turolla, S. Zane

Abstract: The magnetar SGR J1745-2900 discovered at parsecs distance from the Milky Way central black hole, Sagittarius A*, represents the closest pulsar to a supermassive black hole ever detected. Furthermore, its intriguing radio emission has been used to study the environment of the black hole, as well as to derive a precise position and proper motion for this object. The discovery of SGR J1745-2900 has… ▽ More The magnetar SGR J1745-2900 discovered at parsecs distance from the Milky Way central black hole, Sagittarius A*, represents the closest pulsar to a supermassive black hole ever detected. Furthermore, its intriguing radio emission has been used to study the environment of the black hole, as well as to derive a precise position and proper motion for this object. The discovery of SGR J1745-2900 has opened interesting debates about the number, age and nature of pulsars expected in the Galactic center region. In this work, we present extensive X-ray monitoring of the outburst of SGR J1745-2900 using the Chandra X-ray Observatory, the only instrument with the spatial resolution to distinguish the magnetar from the supermassive black hole (2.4" angular distance). It was monitored from its outburst onset in April 2013 until August 2019, collecting more than fifty Chandra observations for a total of more than 2.3 Ms of data. Soon after the outburst onset, the magnetar emission settled onto a purely thermal emission state that cooled from a temperature of about 0.9 to 0.6 keV over 6 years. The pulsar timing properties showed at least two changes in the period derivative, increasing by a factor of about 4 during the outburst decay. We find that the long-term properties of this outburst challenge current models for the magnetar outbursts. △ Less

Submitted 24 March, 2020; v1 submitted 16 March, 2020; originally announced March 2020.

Comments: 11 pages, 6 figures. Accepted by ApJ

Giant planet swaps during close stellar encounters

Authors: Yi-Han Wang, Rosalba Perna, Nathan W. C. Leigh

Abstract: The discovery of planetary systems outside of the solar system has challenged some of the tenets of planetary formation. Among the difficult-to-explain observations, are systems with a giant planet orbiting a very-low mass star, such as the recently discovered GJ~3512b planetary system, where a Jupiter-like planet orbits an $M$-star in a tight and eccentric orbit. Systems such as this one are not… ▽ More The discovery of planetary systems outside of the solar system has challenged some of the tenets of planetary formation. Among the difficult-to-explain observations, are systems with a giant planet orbiting a very-low mass star, such as the recently discovered GJ~3512b planetary system, where a Jupiter-like planet orbits an $M$-star in a tight and eccentric orbit. Systems such as this one are not predicted by the core accretion theory of planet formation. Here we suggest a novel mechanism, in which the giant planet is born around a more typical Sun-like star ($M_{*,1}$), but is subsequently exchanged during a dynamical interaction with a flyby low-mass star ($M_{*,2}$). We perform state-of-the-art $N$-body simulations with $M_{*,1}=1M_\odot$ and $M_{*,2}=0.1M_\odot$ to study the statistical outcomes of this interaction, and show that exchanges result in high eccentricities for the new orbit around the low-mass star, while about half of the outcomes result in tighter orbits than the planet had around its birth star. We numerically compute the cross section for planet exchange, and show that an upper limit for the probability per planetary system to have undergone such an event is $Γ\sim 4.4(M_{\rm c}/100M_\odot)^{-2}(a_{\rm p}/{\rm AU}) (σ/1\,{\rm km}\,{\rm s}^{-1})^{5}$Gyr$^{-1}$, where $a_{\rm p}$ is the planet semi-major axis around the birth star, $σ$ the velocity dispersion of the star cluster, and $M_{\rm c}$ the total mass of the star cluster. Hence these planet exchanges could be relatively common for stars born in open clusters and groups, should already be observed in the exoplanet database, and provide new avenues to create unexpected planetary architectures. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: Accepted to ApJ Letters