Showing 1–41 of 41 results for author: Bortolas, E

An effective model for the tidal disruption of satellites undergoing minor mergers with axisymmetric primaries

Authors: Ludovica Varisco, Massimo Dotti, Matteo Bonetti, Elisa Bortolas, Alessandro Lupi

Abstract: According to the hierarchical formation paradigm, galaxies form through mergers of smaller entities and massive black holes (MBHs), if lurking at their centers, migrate to the nucleus of the newly formed galaxy, where they form binary systems. The formation and evolution of MBH binaries, and in particular their coalescence timescale, is very relevant for current and future facilities aimed at dete… ▽ More According to the hierarchical formation paradigm, galaxies form through mergers of smaller entities and massive black holes (MBHs), if lurking at their centers, migrate to the nucleus of the newly formed galaxy, where they form binary systems. The formation and evolution of MBH binaries, and in particular their coalescence timescale, is very relevant for current and future facilities aimed at detecting the gravitational-wave signal produced by the MBH close to coalescence. While most of the studies targeting this process are based on hydrodynamic simulations, the high computational cost makes a complete parameter space exploration prohibitive. Semi-analytic approaches represent a valid alternative, but they require ad-hoc prescriptions for the mass loss of the merging galaxies in minor mergers due to tidal stripping, which is not commonly considered or at most modelled assuming very idealised geometries. In this work, we propose a novel, effective model for the tidal stripping in axisymmetric potentials, to be implemented in semi-analytic models. We validate our semi-analytic approach against N-body simulations considering different galaxy sizes, inclinations, and eccentricities, finding only a moderate dependence on the orbit eccentricity. In particular, we find that, for almost circular orbits, our model mildly overestimates the mass loss, and this is due to the adjustment of the stellar distribution after the mass is removed. Nonetheless, the model exhibits a very good agreement with simulations in all the considered conditions, and thus represents an extremely powerful addition to semi-analytic calculations. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 15 pages, 8 figures, 2 tables. Accepted for publication on A&A

Eccentricity evolution of PTA sources from cosmological initial conditions

Authors: Federica Fastidio, Alessia Gualandris, Alberto Sesana, Elisa Bortolas, Walter Dehnen

Abstract: Recent results from pulsar timing arrays (PTAs) show evidence for a gravitational wave background (GWB) consistent with a population of unresolved supermassive black hole (SMBH) binaries (BHBs). While the data do not yet constrain the slope of the spectrum, this appears to flatten at the lowest frequencies, deviating from the power-law shape expected for circular binaries evolving solely due to gr… ▽ More Recent results from pulsar timing arrays (PTAs) show evidence for a gravitational wave background (GWB) consistent with a population of unresolved supermassive black hole (SMBH) binaries (BHBs). While the data do not yet constrain the slope of the spectrum, this appears to flatten at the lowest frequencies, deviating from the power-law shape expected for circular binaries evolving solely due to gravitational wave (GW) emission. Interestingly, such flattening can be explained with a population of eccentric rather than circular binaries. The eccentricity of BHBs is notoriously difficult to predict based simply on the parameters of the host galaxies and the initial galactic orbit, as it is subject to stochastic effects. We study the evolution of the eccentricity of BHBs formed in galactic mergers with cosmological initial conditions from pairing to coalescence, with a focus on potential PTA sources. We select galactic mergers from the IllustrisTNG100-1 simulation and re-simulate them at high resolution with the N-body code Griffin down to binary separations of the order of a parsec. We then estimate coalescence timescales with a semi-analytical model of the evolution under the effects of GW emission and stellar hardening. We find that most mergers in IllustrisTNG100-1 occur on highly eccentric orbits, and that the eccentricity of BHBs at binary formation correlates with the initial eccentricity of the merger, if this is no larger than approximately 0.9. For extremely eccentric mergers, the binaries tend to form with modest eccentricities. We discuss the implications of these results on the interpretation of the observed GWB. △ Less

Submitted 11 June, 2024; v1 submitted 4 June, 2024; originally announced June 2024.

Comments: 10 pages, 11 figures, Accepted for publication in MNRAS

The Lunar Gravitational-wave Antenna: Mission Studies and Science Case

Authors: Parameswaran Ajith, Pau Amaro Seoane, Manuel Arca Sedda, Riccardo Arcodia, Francesca Badaracco, Enis Belgacem, Stefano Benetti, Alexey Bobrick, Alessandro Bonforte, Elisa Bortolas, Valentina Braito, Marica Branchesi, Adam Burrows, Enrico Cappellaro, Roberto Della Ceca, Chandrachur Chakraborty, Shreevathsa Chalathadka Subrahmanya, Michael W. Coughlin, Stefano Covino, Andrea Derdzinski, Aayushi Doshi, Maurizio Falanga, Stefano Foffa, Alessia Franchini, Alessandro Frigeri , et al. (58 additional authors not shown)

Abstract: The Lunar Gravitational-wave Antenna (LGWA) is a proposed array of next-generation inertial sensors to monitor the response of the Moon to gravitational waves (GWs). Given the size of the Moon and the expected noise produced by the lunar seismic background, the LGWA would be able to observe GWs from about 1 mHz to 1 Hz. This would make the LGWA the missing link between space-borne detectors like L… ▽ More The Lunar Gravitational-wave Antenna (LGWA) is a proposed array of next-generation inertial sensors to monitor the response of the Moon to gravitational waves (GWs). Given the size of the Moon and the expected noise produced by the lunar seismic background, the LGWA would be able to observe GWs from about 1 mHz to 1 Hz. This would make the LGWA the missing link between space-borne detectors like LISA with peak sensitivities around a few millihertz and proposed future terrestrial detectors like Einstein Telescope or Cosmic Explorer. In this article, we provide a first comprehensive analysis of the LGWA science case including its multi-messenger aspects and lunar science with LGWA data. We also describe the scientific analyses of the Moon required to plan the LGWA mission. △ Less

Submitted 14 April, 2024; originally announced April 2024.

Repeating partial disruptions and two-body relaxation

Authors: Luca Broggi, Nicholas C. Stone, Taeho Ryu, Elisa Bortolas, Massimo Dotti, Matteo Bonetti, Alberto Sesana

Abstract: Two-body relaxation may drive stars onto near-radial orbits around a massive black hole, resulting in a tidal disruption event (TDE). In some circumstances, stars are unlikely to undergo a single terminal disruption, but rather to have a sequence of many grazing encounters with the black hole. It has long been unclear what is the physical outcome of this sequence: each of these encounters can only… ▽ More Two-body relaxation may drive stars onto near-radial orbits around a massive black hole, resulting in a tidal disruption event (TDE). In some circumstances, stars are unlikely to undergo a single terminal disruption, but rather to have a sequence of many grazing encounters with the black hole. It has long been unclear what is the physical outcome of this sequence: each of these encounters can only liberate a small amount of stellar mass, but may significantly alter the orbit of the star. We study the phenomenon of repeating partial tidal disruptions (pTDEs) by building a semi-analytical model that accounts for mass loss and tidal excitation. In the empty loss cone regime, where two-body relaxation is weak, we estimate the number of consecutive partial disruptions that a star can undergo, on average, before being significantly affected by two-body encounters. We find that in this empty loss cone regime, a star will be destroyed in a sequence of weak pTDEs, possibly explaining the tension between the low observed TDE rate and its higher theoretical estimates. △ Less

Submitted 19 June, 2024; v1 submitted 8 April, 2024; originally announced April 2024.

Comments: Open Journal of Astrophysics, accepted for publication

Demographics of Tidal Disruption Events with L-Galaxies: I. Volumetric TDE rates and the abundance of Nuclear Star Clusters

Authors: M. Polkas, S. Bonoli, E. Bortolas, D. Izquierdo-Villalba, A. Sesana, L. Broggi, N. Hoyer, D. Spinoso

Abstract: Stars can be ripped apart by tidal forces in the vicinity of a massive black hole (MBH), causing luminous flares known as tidal disruption events (TDEs). These events could be contributing to the mass growth of intermediate-mass MBHs, and new samples from transient surveys can provide useful information on this growth channel. This work aims to study the demographics of TDEs by modeling the co-evo… ▽ More Stars can be ripped apart by tidal forces in the vicinity of a massive black hole (MBH), causing luminous flares known as tidal disruption events (TDEs). These events could be contributing to the mass growth of intermediate-mass MBHs, and new samples from transient surveys can provide useful information on this growth channel. This work aims to study the demographics of TDEs by modeling the co-evolution of MBHs and their galactic environments in a cosmological framework. We use the semi-analytic galaxy formation model \emph{L-Galaxies}BH, which follows the evolution of galaxies as well as of MBHs, including multiple scenarios for MBH seeds and growth, spin evolution, and binary MBH dynamics. Time-dependent TDE rates are associated with each MBH depending on the stellar environment, following the solutions to the 1-D Fokker Planck equation solved with \textsc{PhaseFlow}. Our model produces volumetric rates that are in agreement with the latest optical and previous X-ray samples. This agreement requires a high occupation fraction of nuclear star clusters with MBHs since these star reservoirs host the majority of TDEs at all mass regimes. We predict that TDE rates are an increasing function of MBH mass up to $\sim\, 10^{5.5}$M$_\odot$, beyond which the distribution flattens and eventually drops for $>\,10^{7}$M$_\odot$. In general, volumetric rates are predicted to be redshift-independent at $z\,{<}\,1$. We discuss how the spin distribution of MBHs around the event horizon suppression can be constrained via TDE rates and what is the average contribution of TDEs to the MBH growth. In our work, the majority of low-mass galaxies host nuclear star clusters that have their loss-cone depleted by $z\,=\,0$, explaining why TDEs are rare in these systems. This highlights that time-dependent TDE rates are essential for any model to be in good agreement with observations at all mass regimes. △ Less

Submitted 1 June, 2024; v1 submitted 20 December, 2023; originally announced December 2023.

Comments: 26 pages, 11 figures

Searching for Intermediate Mass Black Holes in Globular Clusters Through Tidal Disruption Events

Authors: Vivian L. Tang, Piero Madau, Elisa Bortolas, Eric W. Peng

Abstract: Intermediate mass black holes (IMBHs) may be the link between stellar mass holes and the supermassive variety in the nuclei of galaxies, and globular clusters (GCs) may be one of the most promising environments for their formation. Here we carry out a pilot study of the observability of tidal disruption events (TDEs) from 10^3 Msun < M_BH < 10^5 Msun IMBHs embedded in stellar cusps at the center o… ▽ More Intermediate mass black holes (IMBHs) may be the link between stellar mass holes and the supermassive variety in the nuclei of galaxies, and globular clusters (GCs) may be one of the most promising environments for their formation. Here we carry out a pilot study of the observability of tidal disruption events (TDEs) from 10^3 Msun < M_BH < 10^5 Msun IMBHs embedded in stellar cusps at the center of GCs. We model the long super-Eddington accretion phase and ensuing optical flare, and derive the disruption rate of main-sequence stars as a function of black hole mass and GC properties with the help of a 1D Fokker-Planck approach. The photospheric emission of the adiabatically expanding outflow dominates the observable radiation and peaks in the NUV/optical bands, outshining the brightness of the (old) stellar population of GCs in Virgo for a period of months to years. A search for TDE events in a sample of nearly 4,000 GCs observed at multiple epochs by the Next Generation Virgo Cluster Survey (NGVS) yields null results. Given our model predictions, this sample is too small to set stringent constraints on the present-day occupation fraction of GCs hosting IMBHs. Naturally, better simulations of the properties of the cluster central stellar distribution, TDE light curves and rates, together with larger surveys of GCs are all needed to gain deeper insights into the presence of IMBHs in GCs. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Comments: 12 pages, 9 figures, submitted for publication in The Astrophysical Journal

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

Authors: The International Pulsar Timing Array Collaboration, G. Agazie, J. Antoniadis, A. Anumarlapudi, A. M. Archibald, P. Arumugam, S. Arumugam, Z. Arzoumanian, J. Askew, S. Babak, M. Bagchi, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, A. Bathula, B. Bécsy, A. Berthereau, N. D. R. Bhat, L. Blecha, M. Bonetti, E. Bortolas, A. Brazier, P. R. Brook, M. Burgay , et al. (220 additional authors not shown)

Abstract: The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTA… ▽ More The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1σ$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

The second data release from the European Pulsar Timing Array: VI. Challenging the ultralight dark matter paradigm

Authors: Clemente Smarra, Boris Goncharov, Enrico Barausse, J. Antoniadis, S. Babak, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, G. Desvignes, M. Falxa, R. D. Ferdman, A. Franchini, J. R. Gair, E. Graikou, J. -M. Grie , et al. (46 additional authors not shown)

Abstract: Pulsar Timing Array experiments probe the presence of possible scalar or pseudoscalar ultralight dark matter particles through decade-long timing of an ensemble of galactic millisecond radio pulsars. With the second data release of the European Pulsar Timing Array, we focus on the most robust scenario, in which dark matter interacts only gravitationally with ordinary baryonic matter. Our results s… ▽ More Pulsar Timing Array experiments probe the presence of possible scalar or pseudoscalar ultralight dark matter particles through decade-long timing of an ensemble of galactic millisecond radio pulsars. With the second data release of the European Pulsar Timing Array, we focus on the most robust scenario, in which dark matter interacts only gravitationally with ordinary baryonic matter. Our results show that ultralight particles with masses $10^{-24.0}~\text{eV} \lesssim m \lesssim 10^{-23.3}~\text{eV}$ cannot constitute $100\%$ of the measured local dark matter density, but can have at most local density $ρ\lesssim 0.3$ GeV/cm$^3$. △ Less

Submitted 25 October, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 5 pages + acknowledgements + refs, 2 figures. Text and figures match the version published in PRL

The second data release from the European Pulsar Timing Array: IV. Implications for massive black holes, dark matter and the early Universe

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, P. Auclair, S. Babak, M. Bagchi, A. -S. Bak Nielsen, E. Barausse, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, C. Caprini, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, M. Crisostomi, S. Dandapat, D. Deb , et al. (89 additional authors not shown)

Abstract: The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling sup… ▽ More The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling supermassive black hole binaries (SMBHBs); inflation, phase transitions, cosmic strings and tensor mode generation by non-linear evolution of scalar perturbations in the early Universe; oscillations of the Galactic potential in the presence of ultra-light dark matter (ULDM). At the current stage of emerging evidence, it is impossible to discriminate among the different origins. Therefore, in this paper, we consider each process separately, and investigate the implications of the signal under the hypothesis that it is generated by that specific process. We find that the signal is consistent with a cosmic population of inspiralling SMBHBs, and its relatively high amplitude can be used to place constraints on binary merger timescales and the SMBH-host galaxy scaling relations. If this origin is confirmed, this is the first direct evidence that SMBHBs merge in nature, adding an important observational piece to the puzzle of structure formation and galaxy evolution. As for early Universe processes, the measurement would place tight constraints on the cosmic string tension and on the level of turbulence developed by first-order phase transitions. Other processes would require non-standard scenarios, such as a blue-tilted inflationary spectrum or an excess in the primordial spectrum of scalar perturbations at large wavenumbers. Finally, a ULDM origin of the detected signal is disfavoured, which leads to direct constraints on the abundance of ULDM in our Galaxy. △ Less

Submitted 15 May, 2024; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 23 figures, replaced to match the version published in Astronomy & Astrophysics, note the change in the numbering order in the series (now paper IV)

The second data release from the European Pulsar Timing Array V. Search for continuous gravitational wave signals

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, S. Babak, M. Bagchi, A. S. Bak Nielsen, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, S. Dandapat, D. Deb, S. Desai, G. Desvignes, N. Dhanda-Batra, C. Dwivedi , et al. (75 additional authors not shown)

Abstract: We present the results of a search for continuous gravitational wave signals (CGWs) in the second data release (DR2) of the European Pulsar Timing Array (EPTA) collaboration. The most significant candidate event from this search has a gravitational wave frequency of 4-5 nHz. Such a signal could be generated by a supermassive black hole binary (SMBHB) in the local Universe. We present the results o… ▽ More We present the results of a search for continuous gravitational wave signals (CGWs) in the second data release (DR2) of the European Pulsar Timing Array (EPTA) collaboration. The most significant candidate event from this search has a gravitational wave frequency of 4-5 nHz. Such a signal could be generated by a supermassive black hole binary (SMBHB) in the local Universe. We present the results of a follow-up analysis of this candidate using both Bayesian and frequentist methods. The Bayesian analysis gives a Bayes factor of 4 in favor of the presence of the CGW over a common uncorrelated noise process, while the frequentist analysis estimates the p-value of the candidate to be 1%, also assuming the presence of common uncorrelated red noise. However, comparing a model that includes both a CGW and a gravitational wave background (GWB) to a GWB only, the Bayes factor in favour of the CGW model is only 0.7. Therefore, we cannot conclusively determine the origin of the observed feature, but we cannot rule it out as a CGW source. We present results of simulations that demonstrate that data containing a weak gravitational wave background can be misinterpreted as data including a CGW and vice versa, providing two plausible explanations of the EPTA DR2 data. Further investigations combining data from all PTA collaborations will be needed to reveal the true origin of this feature. △ Less

Submitted 25 June, 2024; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 13 figures, 15 pages, accepted

The second data release from the European Pulsar Timing Array II. Customised pulsar noise models for spatially correlated gravitational waves

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, S. Babak, M. Bagchi, A. S. Bak Nielsen, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, S. Dandapat, D. Deb, S. Desai, G. Desvignes, N. Dhanda-Batra, C. Dwivedi , et al. (73 additional authors not shown)

Abstract: The nanohertz gravitational wave background (GWB) is expected to be an aggregate signal of an ensemble of gravitational waves emitted predominantly by a large population of coalescing supermassive black hole binaries in the centres of merging galaxies. Pulsar timing arrays, ensembles of extremely stable pulsars, are the most precise experiments capable of detecting this background. However, the su… ▽ More The nanohertz gravitational wave background (GWB) is expected to be an aggregate signal of an ensemble of gravitational waves emitted predominantly by a large population of coalescing supermassive black hole binaries in the centres of merging galaxies. Pulsar timing arrays, ensembles of extremely stable pulsars, are the most precise experiments capable of detecting this background. However, the subtle imprints that the GWB induces on pulsar timing data are obscured by many sources of noise. These must be carefully characterized to increase the sensitivity to the GWB. In this paper, we present a novel technique to estimate the optimal number of frequency coefficients for modelling achromatic and chromatic noise and perform model selection. We also incorporate a new model to fit for scattering variations in the pulsar timing package temponest and created realistic simulations of the European Pulsar Timing Array (EPTA) datasets that allowed us to test the efficacy of our noise modelling algorithms. We present an in-depth analysis of the noise properties of 25 millisecond pulsars (MSPs) that form the second data release (DR2) of the EPTA and investigate the effect of incorporating low-frequency data from the Indian PTA collaboration. We use enterprise and temponest packages to compare noise models with those reported with the EPTA DR1. We find that, while in some pulsars we can successfully disentangle chromatic from achromatic noise owing to the wider frequency coverage in DR2, in others the noise models evolve in a more complicated way. We also find evidence of long-term scattering variations in PSR J1600$-$3053. Through our simulations, we identify intrinsic biases in our current noise analysis techniques and discuss their effect on GWB searches. The results presented here directly help improve sensitivity to the GWB and are already being used as part of global PTA efforts. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 20 pages, 6 figures, 9 tables

Journal ref: A&A 678, A49 (2023)

The second data release from the European Pulsar Timing Array I. The dataset and timing analysis

Authors: J. Antoniadis, S. Babak, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, G. Desvignes, M. Falxa, R. D. Ferdman, A. Franchini, J. R. Gair, B. Goncharov, E. Graikou, J. -M. Grießmeier, L. Guillemot, Y. J. Guo , et al. (44 additional authors not shown)

Abstract: Pulsar timing arrays offer a probe of the low-frequency gravitational wave spectrum (1 - 100 nanohertz), which is intimately connected to a number of markers that can uniquely trace the formation and evolution of the Universe. We present the dataset and the results of the timing analysis from the second data release of the European Pulsar Timing Array (EPTA). The dataset contains high-precision pu… ▽ More Pulsar timing arrays offer a probe of the low-frequency gravitational wave spectrum (1 - 100 nanohertz), which is intimately connected to a number of markers that can uniquely trace the formation and evolution of the Universe. We present the dataset and the results of the timing analysis from the second data release of the European Pulsar Timing Array (EPTA). The dataset contains high-precision pulsar timing data from 25 millisecond pulsars collected with the five largest radio telescopes in Europe, as well as the Large European Array for Pulsars. The dataset forms the foundation for the search for gravitational waves by the EPTA, presented in associated papers. We describe the dataset and present the results of the frequentist and Bayesian pulsar timing analysis for individual millisecond pulsars that have been observed over the last ~25 years. We discuss the improvements to the individual pulsar parameter estimates, as well as new measurements of the physical properties of these pulsars and their companions. This data release extends the dataset from EPTA Data Release 1 up to the beginning of 2021, with individual pulsar datasets with timespans ranging from 14 to 25 years. These lead to improved constraints on annual parallaxes, secular variation of the orbital period, and Shapiro delay for a number of sources. Based on these results, we derived astrophysical parameters that include distances, transverse velocities, binary pulsar masses, and annual orbital parallaxes. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 29 pages, 9 figures, 13 tables, Astronomy & Astrophysics in press

The second data release from the European Pulsar Timing Array III. Search for gravitational wave signals

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, S. Babak, M. Bagchi, A. -S. Bak Nielsen, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, S. Dandapat, D. Deb, S. Desai, G. Desvignes, N. Dhanda-Batra, C. Dwivedi , et al. (73 additional authors not shown)

Abstract: We present the results of the search for an isotropic stochastic gravitational wave background (GWB) at nanohertz frequencies using the second data release of the European Pulsar Timing Array (EPTA) for 25 millisecond pulsars and a combination with the first data release of the Indian Pulsar Timing Array (InPTA). We analysed (i) the full 24.7-year EPTA data set, (ii) its 10.3-year subset based on… ▽ More We present the results of the search for an isotropic stochastic gravitational wave background (GWB) at nanohertz frequencies using the second data release of the European Pulsar Timing Array (EPTA) for 25 millisecond pulsars and a combination with the first data release of the Indian Pulsar Timing Array (InPTA). We analysed (i) the full 24.7-year EPTA data set, (ii) its 10.3-year subset based on modern observing systems, (iii) the combination of the full data set with the first data release of the InPTA for ten commonly timed millisecond pulsars, and (iv) the combination of the 10.3-year subset with the InPTA data. These combinations allowed us to probe the contributions of instrumental noise and interstellar propagation effects. With the full data set, we find marginal evidence for a GWB, with a Bayes factor of four and a false alarm probability of $4\%$. With the 10.3-year subset, we report evidence for a GWB, with a Bayes factor of $60$ and a false alarm probability of about $0.1\%$ ($\gtrsim 3σ$ significance). The addition of the InPTA data yields results that are broadly consistent with the EPTA-only data sets, with the benefit of better noise modelling. Analyses were performed with different data processing pipelines to test the consistency of the results from independent software packages. The inferred spectrum from the latest EPTA data from new generation observing systems is rather uncertain and in mild tension with the common signal measured in the full data set. However, if the spectral index is fixed at 13/3, the two data sets give a similar amplitude of ($2.5\pm0.7)\times10^{-15}$ at a reference frequency of $1\,{\rm yr}^{-1}$. By continuing our detection efforts as part of the International Pulsar Timing Array (IPTA), we expect to be able to improve the measurement of spatial correlations and better characterise this signal in the coming years. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 21 pages, 14 figures, 4 appendix figures, accepted for publication in A&A

Journal ref: A&A 678, A50 (2023)

Practical approaches to analyzing PTA data: Cosmic strings with six pulsars

Authors: Hippolyte Quelquejay Leclere, Pierre Auclair, Stanislav Babak, Aurélien Chalumeau, Danièle A. Steer, J. Antoniadis, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, G. Desvignes, M. Falxa, R. D. Ferdman, A. Franchini, J. R. Gair, B. Goncharov, E. Graikou , et al. (47 additional authors not shown)

Abstract: We search for a stochastic gravitational wave background (SGWB) generated by a network of cosmic strings using six millisecond pulsars from Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA). We perform a Bayesian analysis considering two models for the network of cosmic string loops, and compare it to a simple power-law model which is expected from the population of supermassive blac… ▽ More We search for a stochastic gravitational wave background (SGWB) generated by a network of cosmic strings using six millisecond pulsars from Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA). We perform a Bayesian analysis considering two models for the network of cosmic string loops, and compare it to a simple power-law model which is expected from the population of supermassive black hole binaries. Our main strong assumption is that the previously reported common red noise process is a SGWB. We find that the one-parameter cosmic string model is slightly favored over a power-law model thanks to its simplicity. If we assume a two-component stochastic signal in the data (supermassive black hole binary population and the signal from cosmic strings), we get a $95\%$ upper limit on the string tension of $\log_{10}(Gμ) < -9.9$ ($-10.5$) for the two cosmic string models we consider. In extended two-parameter string models, we were unable to constrain the number of kinks. We test two approximate and fast Bayesian data analysis methods against the most rigorous analysis and find consistent results. These two fast and efficient methods are applicable to all SGWBs, independent of their source, and will be crucial for analysis of extended data sets. △ Less

Submitted 3 May, 2024; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: 14 pages, 6 figures; typo corrected in (5)

Journal ref: Phys. Rev. D 108 (2023), 123527

Quasi-periodic eruptions from impacts between the secondary and a rigidly precessing accretion disc in an extreme mass-ratio inspiral system

Authors: Alessia Franchini, Matteo Bonetti, Alessandro Lupi, Giovanni Miniutti, Elisa Bortolas, Margherita Giustini, Massimo Dotti, Alberto Sesana, Riccardo Arcodia, Taeho Ryu

Abstract: X-ray quasi-periodic eruptions (QPEs) represent a recently discovered example of extreme X-ray variability associated with supermassive black holes. These are high-amplitude bursts recurring every few hours that are detected in the soft X-ray band from the nuclei of nearby galaxies whose optical spectra lack the broad emission lines typically observed in unobscured active galaxies. The physical or… ▽ More X-ray quasi-periodic eruptions (QPEs) represent a recently discovered example of extreme X-ray variability associated with supermassive black holes. These are high-amplitude bursts recurring every few hours that are detected in the soft X-ray band from the nuclei of nearby galaxies whose optical spectra lack the broad emission lines typically observed in unobscured active galaxies. The physical origin of this new X-ray variability phenomenon is still unknown and several theoretical models have been presented. However, no attempt has been made so far to account for the varying QPE recurrence time and luminosity in individual sources, nor for the diversity of the QPE phenomenology in the different known erupters. We present a semi-analytical model based on an extreme mass-ratio inspiral (EMRI) system where the secondary intersects, along its orbit, a rigidly precessing accretion disc surrounding the primary. We assume that QPEs result from emission from an adiabatically expanding, initially optically thick gas cloud expelled from the disc plane at each impact. We produced synthetic X-ray light curves, which we then compared with X-ray data from four QPE sources: GSN 069, eRO-QPE1, eRO-QPE2, and RX J1301.9+2747. Our model aptly reproduces the diversity of QPE properties between the considered objects and it is also able to naturally account for the varying QPE amplitudes and recurrence times in individual sources. Future implementations will enable us to refine the match with the data and to estimate the system parameters precisely, making additional use of multi-epoch QPE data. We briefly discuss the nature of the secondary object, as well as the possible implications of our findings for the EMRI population at large. △ Less

Submitted 14 June, 2023; v1 submitted 3 April, 2023; originally announced April 2023.

Comments: 11 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 675, A100 (2023)

Partial stellar tidal disruption events and their rates

Authors: Elisa Bortolas, Taeho Ryu, Luca Broggi, Alberto Sesana

Abstract: Tidal disruption events (TDEs) of stars operated by massive black holes (MBHs) will be detected in thousands by upcoming facilities such as the Vera Rubin Observatory. In this work, we assess the rates of standard total TDEs, destroying the entire star, and partial TDEs, in which a stellar remnant survives the interaction, by solving 1-D Fokker-Planck equations. Our rate estimates are based on a n… ▽ More Tidal disruption events (TDEs) of stars operated by massive black holes (MBHs) will be detected in thousands by upcoming facilities such as the Vera Rubin Observatory. In this work, we assess the rates of standard total TDEs, destroying the entire star, and partial TDEs, in which a stellar remnant survives the interaction, by solving 1-D Fokker-Planck equations. Our rate estimates are based on a novel definition of the loss cone whose size is commensurate to the largest radius at which partial disruptions can occur, as motivated by relativistic hydrodynamical simulations. Our novel approach unveils two important results. First, partial TDEs can be more abundant than total disruptions by a factor of a few to a few tens. Second, the rates of complete stellar disruptions can be overestimated by a factor of a few to a few tens if one neglects partial TDEs, as we find that many of the events classified at total disruptions in the standard framework are in fact partial TDEs. Accounting for partial TDEs is particularly relevant for galaxies harbouring a nuclear stellar cluster featuring many events coming from the empty loss cone. Based on these findings, we stress that partial disruptions should be considered when constraining the luminosity function of TDE flares; accounting for this may reconcile the theoretically estimated TDE rates with the observed ones. △ Less

Submitted 4 July, 2023; v1 submitted 6 March, 2023; originally announced March 2023.

Comments: 12 pages + Appendix, MNRAS, accepted

MSC Class: 85A05

Extreme mass ratio inspirals and tidal disruption events in nuclear clusters. I. Time dependent rates

Authors: Luca Broggi, Elisa Bortolas, Matteo Bonetti, Alberto Sesana, Massimo Dotti

Abstract: In this paper we develop a computationally efficient, two-population, time-dependent Fokker-Plank approach in the two dimensions of energy and angular momentum to study the rates of tidal disruption events (TDEs), extreme mass ratio inspirals (EMRIs) and direct plunges occurring around massive black holes (MBHs) in galactic nuclei. We test our code by exploring a wide range of the astrophysically… ▽ More In this paper we develop a computationally efficient, two-population, time-dependent Fokker-Plank approach in the two dimensions of energy and angular momentum to study the rates of tidal disruption events (TDEs), extreme mass ratio inspirals (EMRIs) and direct plunges occurring around massive black holes (MBHs) in galactic nuclei. We test our code by exploring a wide range of the astrophysically relevant parameter space, including MBH masses, galaxy central densities and inner density slopes. We find that mass segregation and, more in general, the time dependency of the distribution function regulate the event rate: TDEs always decline with time, whereas EMRIs and plunges reach a maximum and undergo a subsequent nearly exponential decay. Once suitably normalized, the rates associated to different choices of MBH mass and galaxy density overlap nearly perfectly. Based on this, we provide a simple scaling that allows to reproduce the time-dependent event rates for any MBH mass and underlying galactic nucleus. Although our peak rates are in general agreement with the literature relying on the steady-state (non-time dependent) assumption, those can be sustained on a timescale that strongly depends on the properties of the system. In particular this can be much shorter than a Gyr for relatively light MBHs residing in dense systems. This warns against using steady state models to compute global TDE, EMRI and plunge rates and calls for a more sophisticated, time dependent treatment of the problem. △ Less

Submitted 12 May, 2022; originally announced May 2022.

Comments: 14 pages, 10 figures, submitted to MNRAS

Optical follow-up of the tick-tock massive black hole binary candidate

Authors: Massimo Dotti, Matteo Bonetti, Fabio Rigamonti, Elisa Bortolas, Matteo Fossati, Roberto Decarli, Stefano Covino, Alessandro Lupi, Alessia Franchini, Alberto Sesana, Giorgio Calderone

Abstract: The observation of a population of massive black hole binaries (MBHBs) is key for our complete understanding of galaxy mergers and for the characterization of the expected gravitational waves (GWs) signal. However, MBHBs still remain elusive with only a few candidates proposed to date. Among these, SDSSJ143016.05+230344.4 ('tick-tock' hereafter) is the only candidate with a remarkably well sampled… ▽ More The observation of a population of massive black hole binaries (MBHBs) is key for our complete understanding of galaxy mergers and for the characterization of the expected gravitational waves (GWs) signal. However, MBHBs still remain elusive with only a few candidates proposed to date. Among these, SDSSJ143016.05+230344.4 ('tick-tock' hereafter) is the only candidate with a remarkably well sampled light curve showing a clear reduction of the modulation period and amplitude over three years of observations. This particular feature has been recently claimed to be the signature of a MBHB that is about to merge. In this paper, we provide an optical follow-up of the tick-tock source using the Rapid Eye Mount (REM) telescope. The decreasing luminosity observed in our follow up is hardly explained within the binary scenario. We speculate about an alternative scenario that might explain the observed light curve through relativistic Lense-Thirring precession of an accretion disc around a single massive black hole. △ Less

Submitted 14 November, 2022; v1 submitted 12 May, 2022; originally announced May 2022.

Comments: 9 pages, 2 figure, accepted for publication in MNRAS

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

Eccentricity evolution of massive black hole binaries from formation to coalescence

Authors: Alessia Gualandris, Fazeel Mahmood Khan, Elisa Bortolas, Matteo Bonetti, Alberto Sesana, Peter Berczik, Kelly Holley-Bockelmann

Abstract: Coalescing supermassive black hole binaries (BHBs) are expected to be the loudest sources of gravitational waves (GWs) in the Universe. Detection rates for ground or space-based detectors based on cosmological simulations and semi-analytic models are highly uncertain. A major difficulty stems from the necessity to model the BHB from the scale of the merger to that of inspiral. Of particular releva… ▽ More Coalescing supermassive black hole binaries (BHBs) are expected to be the loudest sources of gravitational waves (GWs) in the Universe. Detection rates for ground or space-based detectors based on cosmological simulations and semi-analytic models are highly uncertain. A major difficulty stems from the necessity to model the BHB from the scale of the merger to that of inspiral. Of particular relevance to the GW merger timescale is the binary eccentricity. Here we present a self-consistent numerical study of the eccentricity of BHBs formed in massive gas-free mergers from the early stages of the merger to the hardening phase, followed by a semi-analytical model down to coalescence. We find that the early eccentricity of the unbound black hole pair is largely determined by the initial orbit. It systematically decreases during the dynamical friction phase. The eccentricity at binary formation is affected by stochasticity and noise owing to encounters with stars, but preserves a strong correlation with the initial orbital eccentricity. Binding of the black holes is a phase characterised by strong perturbations, and we present a quantitative definition of the time of binary formation. During hardening the eccentricity increases in minor mergers, unless the binary is approximately circular, but remains largely unchanged in major mergers, in agreement with predictions from semi-analytical models based on isotropic scattering experiments. Coalescence times due to hardening and GW emission in gas-poor non-rotating ellipticals are <~0.5 Gyr for the large initial eccentricities (0.5 < e < 0.9) typical of galaxy mergers in cosmological simulations. △ Less

Submitted 21 January, 2022; originally announced January 2022.

Comments: 15 pages, 13 figures, accepted for publication in MNRAS

Tidal disruption events in post-starburst galaxies: the importance of a complete stellar mass function

Authors: Elisa Bortolas

Abstract: A tidal disruption event (TDE) occurs when a star is destroyed by the strong tidal shear of a massive black hole (MBH). The accumulation of TDE observations over the last years has revealed that post-starburst galaxies are significantly overrepresented in the sample of TDE hosts. Here we address the post-starburst preference by investigating the decline of TDE rates in a Milky-Way like nuclear ste… ▽ More A tidal disruption event (TDE) occurs when a star is destroyed by the strong tidal shear of a massive black hole (MBH). The accumulation of TDE observations over the last years has revealed that post-starburst galaxies are significantly overrepresented in the sample of TDE hosts. Here we address the post-starburst preference by investigating the decline of TDE rates in a Milky-Way like nuclear stellar cluster featuring either a monochromatic (1 M$\odot$) or a complete, evolved stellar mass function. In the former case, the decline of TDE rates with time is very mild, and generally up to a factor of a few in 10 Gyr. Conversely, if a complete mass function is considered, a strong TDE burst over the first 0.1-1 Gyr is followed by a considerable rate drop, by at least an order of magnitude over 10 Gyr. The decline starts after a mass segregation time-scale, and it is more pronounced assuming a more top-heavy initial mass function and/or an initially denser nucleus. Our results thus suggest that the post-starburst preference can be accounted for in realistic systems featuring a complete stellar mass function, even in moderately dense galactic nuclei. Overall, our findings support the idea that starbursting galactic nuclei are characterized by a top-heavy initial mass function; we speculate that accounting for this can reconcile the discrepancy between observed and theoretically predicted TDE rates even in quiescent galaxies. △ Less

Submitted 14 February, 2022; v1 submitted 20 January, 2022; originally announced January 2022.

Comments: 12 pages, 7 figures, 1 table. Accepted for publication in MNRAS. Corrected typo in Equation (1)

MSC Class: 85A05

The competing effect of gas and stars in massive black hole binaries evolution

Authors: Elisa Bortolas, Alessia Franchini, Matteo Bonetti, Alberto Sesana

Abstract: Massive black hole binaries are predicted to form during the hierarchical assembly of cosmic structures and will represent the loudest sources of low-frequency gravitational waves (GWs) detectable by present and forthcoming GW experiments. Before entering the GW-driven regime, their evolution is driven by the interaction with the surrounding stars and gas. While stellar interactions are found to a… ▽ More Massive black hole binaries are predicted to form during the hierarchical assembly of cosmic structures and will represent the loudest sources of low-frequency gravitational waves (GWs) detectable by present and forthcoming GW experiments. Before entering the GW-driven regime, their evolution is driven by the interaction with the surrounding stars and gas. While stellar interactions are found to always shrink the binary, recent studies predict the possibility of binary outspiral mediated by the presence of a gaseous disk, which could endlessly delay the coalescence and impact the merger rates of massive binaries. Here we implement a semi-analytical treatment that follows the binary evolution under the combined effect of stars and gas. We find that binaries may outspiral only if they accrete near or above their Eddington limit and only until their separation reaches the gaseous disk self-gravitating radius. Even in case of an outspiral, the binary eventually reaches a large enough mass for GW to take over and drive it to coalescence. The combined action of stellar hardening, mass growth and GW-driven inspiral brings binaries to coalescence in few hundreds Myr at most, implying that gas-driven expansion will not severely affect the detection prospects of upcoming GW facilities. △ Less

Submitted 30 August, 2021; originally announced August 2021.

Comments: 11 pages, 6 figures, Published in ApJ Letters

MSC Class: 85A05

Journal ref: ApJL 918 L15 (2021)

Revised event rates for extreme and extremely large mass-ratio inspirals

Authors: Verónica Vázquez-Aceves, Lorenz Zwick, Elisa Bortolas, Pedro R. Capelo, Pau Amaro-Seoane, Lucio Mayer, Xian Chen

Abstract: One of the main targets of the Laser Interferometer Space Antenna (LISA) is the detection of extreme mass-ratio inspirals (EMRIs) and extremely large mass-ratio inspirals (X-MRIs). Their orbits are expected to be highly eccentric and relativistic when entering the LISA band. Under these circumstances, the inspiral time-scale given by Peters' formula loses precision and the shift of the last-stable… ▽ More One of the main targets of the Laser Interferometer Space Antenna (LISA) is the detection of extreme mass-ratio inspirals (EMRIs) and extremely large mass-ratio inspirals (X-MRIs). Their orbits are expected to be highly eccentric and relativistic when entering the LISA band. Under these circumstances, the inspiral time-scale given by Peters' formula loses precision and the shift of the last-stable orbit (LSO) caused by the massive black hole spin could influence the event rates estimate. We re-derive EMRIs and X-MRIs event rates by implementing two different versions of a Kerr loss-cone angle that includes the shift in the LSO, and a corrected version of Peters' time-scale that accounts for eccentricity evolution, 1.5 post-Newtonian hereditary fluxes, and spin-orbit coupling. The main findings of our study are summarized as follows: (1) implementing a Kerr loss-cone changes the event rates by a factor ranging between 0.9 and 1.1; (2) the high-eccentricity limit of Peters' formula offers a reliable inspiral time-scale for EMRIs and X-MRIs, resulting in an event rate estimate that deviates by a factor of about 0.9 to 3 when compared to event rates computed with the corrected version of Peters' time-scale and the usual loss-cone definition. (3) Event rates estimates for systems with a wide range of eccentricities should be revised. Peters' formula overestimates the inspiral rates of highly eccentric systems by a factor of about 8 to 30 compared to the corrected values. Besides, for e$_0 \lesssim$0.8, implementing the corrected version of Peters' formula would be necessary to obtain accurate estimates. △ Less

Submitted 1 December, 2021; v1 submitted 30 July, 2021; originally announced August 2021.

Comments: 14 pages, 10 figures

Stellar hardening of massive black hole binaries: the impact of the host rotation

Authors: Ludovica Varisco, Elisa Bortolas, Massimo Dotti, Alberto Sesana

Abstract: Massive black hole binaries (MBHBs) with masses of ~ 10^4 to ~ 10^10 of solar masses are one of the main targets for currently operating and forthcoming space-borne gravitational wave observatories. In this paper, we explore the effect of the stellar host rotation on the bound binary hardening efficiency, driven by three-body stellar interactions. As seen in previous studies, we find that the cent… ▽ More Massive black hole binaries (MBHBs) with masses of ~ 10^4 to ~ 10^10 of solar masses are one of the main targets for currently operating and forthcoming space-borne gravitational wave observatories. In this paper, we explore the effect of the stellar host rotation on the bound binary hardening efficiency, driven by three-body stellar interactions. As seen in previous studies, we find that the centre of mass (CoM) of a prograde MBHB embedded in a rotating environment starts moving on a nearly circular orbit about the centre of the system shortly after the MBHB binding. In our runs, the oscillation radius is approximately 0.25 ( approximately 0.1) times the binary influence radius for equal mass MBHBs (MBHBs with mass ratio 1:4). Conversely, retrograde binaries remain anchored about the centre of the host. The binary shrinking rate is twice as fast when the binary CoM exhibits a net orbital motion, owing to a more efficient loss cone repopulation even in our spherical stellar systems. We develop a model that captures the CoM oscillations of prograde binaries; we argue that the CoM angular momentum gain per time unit scales with the internal binary angular momentum, so that most of the displacement is induced by stellar interactions occurring around the time of MBHB binding, while the subsequent angular momentum enhancement gets eventually quashed by the effect of dynamical friction. The effect of the background rotation on the MBHB evolution may be relevant for LISA sources, that are expected to form in significantly rotating stellar systems. △ Less

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

Comments: 10 pages, 9 figures, 2 tables. Accepted for publication on MNRAS

The role of bars on the dynamical-friction driven inspiral of massive perturbers

Authors: Elisa Bortolas, Matteo Bonetti, Massimo Dotti, Alessandro Lupi, Pedro R. Capelo, Lucio Mayer, Alberto Sesana

Abstract: In this paper, we systematically explore the impact of a galactic bar on the inspiral time-scale of a massive object (MO) within a Milky Way-like galaxy. We integrate the orbit of MOs in a multi-component galaxy model via a semi-analytical approach that accounts for dynamical friction generalized to rotationally supported backgrounds. We compare the MO evolution in a galaxy featuring a Milky Way-l… ▽ More In this paper, we systematically explore the impact of a galactic bar on the inspiral time-scale of a massive object (MO) within a Milky Way-like galaxy. We integrate the orbit of MOs in a multi-component galaxy model via a semi-analytical approach that accounts for dynamical friction generalized to rotationally supported backgrounds. We compare the MO evolution in a galaxy featuring a Milky Way-like rotating bar to the evolution within an analogous axisymmetric galaxy without the bar. In agreement with previous studies, we find that the bar presence may significantly affect the inspiral, sometimes making it shorter by a factor of a few, sometimes hindering it for a Hubble time. The erratic behaviour is mainly impacted by the relative phase at which the MO encounters the stronger bar-induced resonances. In particular, the effect of the bar is more prominent for initially in-plane, prograde MOs, especially those crossing the bar co-rotation radius or outer Lindblad resonance. In the barred galaxy, we find the sinking of the most massive MOs (>~10^7.5 Msun) approaching the galaxy from large separations (>~8 kpc) to be most efficiently hampered. Neglecting the effect of global torques associated with the non-symmetric mass distribution is thus not advisable even within an idealized, smooth galaxy model; we further note that spiral patterns are unlikely to affect the inspiral due to their transient and fluctuating nature. We speculate that the sinking efficiency of massive black holes involved in minor galaxy mergers may be hampered in barred galaxies, making them less likely to host a gravitational wave signal accessible to low-frequency detectors. △ Less

Submitted 7 March, 2022; v1 submitted 12 March, 2021; originally announced March 2021.

Comments: 19 pages, 12 figures, 2 tables. MNRAS, accepted

MSC Class: 85A05

Improved gravitational radiation time-scales II: spin-orbit contributions and environmental perturbations

Authors: Lorenz Zwick, Pedro R. Capelo, Elisa Bortolas, Veronica Vazquez-Aceves, Lucio Mayer, Pau Amaro-Seoane

Abstract: Peters' formula is an analytical estimate of the time-scale of gravitational wave (GW)-induced coalescence of binary systems. It is used in countless applications, where the convenience of a simple formula outweighs the need for precision. However, many promising sources of the Laser Interferometer Space Antenna (LISA), such as supermassive black hole binaries and extreme mass-ratio inspirals (EMR… ▽ More Peters' formula is an analytical estimate of the time-scale of gravitational wave (GW)-induced coalescence of binary systems. It is used in countless applications, where the convenience of a simple formula outweighs the need for precision. However, many promising sources of the Laser Interferometer Space Antenna (LISA), such as supermassive black hole binaries and extreme mass-ratio inspirals (EMRIs), are expected to enter the LISA band with highly eccentric ($e \gtrsim$ 0.9) and highly relativistic orbits. These are exactly the two limits in which Peters' estimate performs the worst. In this work, we expand upon previous results and give simple analytical fits to quantify how the inspiral time-scale is affected by the relative 1.5 post-Newtonian (PN) hereditary fluxes and spin-orbit couplings. We discuss several cases that demand a more accurate GW time-scale. We show how this can have a major influence on quantities that are relevant for LISA event-rate estimates, such as the EMRI critical semi-major axis. We further discuss two types of environmental perturbations that can play a role in the inspiral phase: the gravitational interaction with a third massive body and the energy loss due to dynamical friction and torques from a surrounding gas medium ubiquitous in galactic nuclei. With the aid of PN corrections to the time-scale in vacuum, we find simple analytical expressions for the regions of phase space in which environmental perturbations are of comparable strength to the effects of any particular PN order, being able to qualitatively reproduce the results of much more sophisticated analyses. △ Less

Submitted 24 June, 2021; v1 submitted 29 January, 2021; originally announced February 2021.

Comments: Accepted for publication in MNRAS. Comments welcome!

Generation of gravitational waves and tidal disruptions in clumpy galaxies

Authors: Boris Pestoni, Elisa Bortolas, Pedro R. Capelo, Lucio Mayer

Abstract: Obtaining a better understanding of intermediate-mass black holes (IMBHs) is crucial, as their properties could shed light on the origin and growth of their supermassive counterparts. Massive star-forming clumps, which are present in a large fraction of massive galaxies at $z \sim$ 1-3, are amongst the venues wherein IMBHs could reside. We perform a series of Fokker-Planck simulations to explore t… ▽ More Obtaining a better understanding of intermediate-mass black holes (IMBHs) is crucial, as their properties could shed light on the origin and growth of their supermassive counterparts. Massive star-forming clumps, which are present in a large fraction of massive galaxies at $z \sim$ 1-3, are amongst the venues wherein IMBHs could reside. We perform a series of Fokker-Planck simulations to explore the occurrence of tidal disruption (TD) and gravitational wave (GW) events about an IMBH in a massive star-forming clump, modelling the latter so that its mass ($10^8 \,{\rm M}_{\odot}$) and effective radius ($100$ pc) are consistent with the properties of both observed and simulated clumps. We find that the TD and GW event rates are in the ranges $10^{-6}$-$10^{-5}$ and $10^{-8}$-$10^{-7}$ yr$^{-1}$, respectively, depending on the assumptions for the initial inner density profile of the system ($ρ\propto r^{-2}$ or $\propto r^{-1}$) and the initial mass of the central IMBH ($10^5$ or $10^3\,{\rm M}_{\odot}$). By integrating the GW event rate over $z$ = 1-3, we expect that the Laser Interferometer Space Antenna will be able to detect $\sim$2 GW events per yr coming from these massive clumps; the intrinsic rate of TD events from these systems amounts instead to a few $10^3$ per yr, a fraction of which will be observable by, e.g. the Square Kilometre Array and the Advanced Telescope for High Energy Astrophysics. In conclusion, our results support the idea that the forthcoming GW and electromagnetic facilities may have the unprecedented opportunity of unveiling the lurking population of IMBHs. △ Less

Submitted 15 December, 2020; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: 13 pages, 4 figures, 4 tables, added references

Dynamical evolution of massive perturbers in realistic multi-component galaxy models I: implementation and validation

Authors: Matteo Bonetti, Elisa Bortolas, Alessandro Lupi, Massimo Dotti

Abstract: Galaxies are self-gravitating structures composed by several components encompassing spherical, axial and triaxial symmetry. Although real systems feature heterogeneous components whose properties are intimately connected, semi-analytical approaches often exploit the linearity of the Poisson's equation to represent the potential and mass distribution of a multi-component galaxy as the sum of the i… ▽ More Galaxies are self-gravitating structures composed by several components encompassing spherical, axial and triaxial symmetry. Although real systems feature heterogeneous components whose properties are intimately connected, semi-analytical approaches often exploit the linearity of the Poisson's equation to represent the potential and mass distribution of a multi-component galaxy as the sum of the individual components. In this work, we expand the semi-analytical framework developed in Bonetti et al. (2020) by including both a detailed implementation of the gravitational potential of exponential disc (modelled with a ${\rm sech}^2$ and an exponential vertical profile) and an accurate prescription for the dynamical friction experienced by massive perturbers in composite galaxy models featuring rotating disc structures. Such improvements allow us to evolve arbitrary orbits either within or outside the galactic disc plane. We validate the results obtained by our numerical model against public semi-analytical codes as well as full N-body simulations, finding that our model is in excellent agreement to the codes it is compared with. The ability to reproduce the relevant physical processes responsible for the evolution of massive perturber orbits and its computational efficiency make our framework perfectly suited for large parameter-space exploration studies. △ Less

Submitted 23 January, 2021; v1 submitted 16 October, 2020; originally announced October 2020.

Comments: 17 pages, 9 figures; accepted for publication in MNRAS

Global torques and stochasticity as the drivers of massive black hole pairing in the young Universe

Authors: Elisa Bortolas, Pedro R. Capelo, Tommaso Zana, Lucio Mayer, Matteo Bonetti, Massimo Dotti, Melvyn B. Davies, Piero Madau

Abstract: The forthcoming Laser Interferometer Space Antenna (LISA) will probe the population of coalescing massive black hole (MBH) binaries up to the onset of structure formation. Here we simulate the galactic-scale pairing of $\sim10^6 M_\odot$ MBHs in a typical, non-clumpy main-sequence galaxy embedded in a cosmological environment at $z = 7-6$. In order to increase our statistical sample, we adopt a st… ▽ More The forthcoming Laser Interferometer Space Antenna (LISA) will probe the population of coalescing massive black hole (MBH) binaries up to the onset of structure formation. Here we simulate the galactic-scale pairing of $\sim10^6 M_\odot$ MBHs in a typical, non-clumpy main-sequence galaxy embedded in a cosmological environment at $z = 7-6$. In order to increase our statistical sample, we adopt a strategy that allows us to follow the evolution of six secondary MBHs concomitantly. We find that the magnitude of the dynamical-friction induced torques is significantly smaller than that of the large-scale, stochastic gravitational torques arising from the perturbed and morphologically evolving galactic disc, suggesting that the standard dynamical friction treatment is inadequate for realistic galaxies at high redshift. The dynamical evolution of MBHs is very stochastic, and a variation in the initial orbital phase can lead to a drastically different time-scale for the inspiral. Most remarkably, the development of a galactic bar in the host system either significantly accelerates the inspiral by dragging a secondary MBH into the centre, or ultimately hinders the orbital decay by scattering the MBH in the galaxy outskirts. The latter occurs more rarely, suggesting that galactic bars overall promote MBH inspiral and binary coalescence. The orbital decay time can be an order of magnitude shorter than what would be predicted relying on dynamical friction alone. The stochasticity, and the important role of global torques, have crucial implications for the rates of MBH coalescences in the early Universe: both have to be accounted for when making predictions for the upcoming LISA observatory. △ Less

Submitted 25 August, 2020; v1 submitted 5 May, 2020; originally announced May 2020.

Comments: Accepted for publication in MNRAS; 15 pages, 10 Figures, 2 Tables

MSC Class: 85A05

The erratic path to coalescence of LISA massive black hole binaries in sub-pc resolution simulations of smooth circumnuclear gas disks

Authors: Rafael Souza Lima, Lucio Mayer, Pedro R. Capelo, Elisa Bortolas, Thomas R. Quinn

Abstract: We report on high-resolution simulations that explore the orbital decay of massive black hole (MBH) pairs with masses between $10^5$ and $10^7 M_{\odot}$ embedded in a circumnuclear gas disk (CND). An adiabatic equation of state is adopted, with a range of adiabatic indices, which maintains a smooth flow. Mergers between MBHs in this mass range would be detectable by the upcoming Laser Inteferomet… ▽ More We report on high-resolution simulations that explore the orbital decay of massive black hole (MBH) pairs with masses between $10^5$ and $10^7 M_{\odot}$ embedded in a circumnuclear gas disk (CND). An adiabatic equation of state is adopted, with a range of adiabatic indices, which maintains a smooth flow. Mergers between MBHs in this mass range would be detectable by the upcoming Laser Inteferometer Space Antenna (LISA). The orbital evolution is followed from the CND scale ($100$~pc) down to separations of $0.1$--$0.01$~pc at which a circumbinary disk (CBD) could form. The decay is erratic and strongly dependent on the gas flow within the disk, that ultimately determines the net torques experienced by the sinking MBH. Overall, we can identify three different evolutionary stages: (i) an initially slow decay that leads to no significant change in the orbital angular momentum, resulting in some circularization; (ii) a fast migration phase in which the orbital angular momentum decreases rapidly; and (iii) a final, very slow decay phase, in which orbital angular momentum can even increase, and a CBD can form. The fast migration phase owes to disk-driven torques originating primarily from the co-orbital region of the secondary MBH, at a distance of 1--3 Hill radii. We find strong analogies with fast Type III migration for massive planets in protoplanetary disks. The CBD forms only when the decay rate becomes small enough to allow it enough time to carve a cavity around the primary MBH, at scales $\lesssim 1$~pc; when this happens, the MBH separation nearly stalls in our higher-resolution run. We suggest an empirically modified gap-opening criterion that takes into account such timescale effects as well as other deviations from standard assumptions made in the literature. [Abriged] △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: 28 pages, 16 figures, 1 table, submitted to ApJ

Dynamical friction-driven orbital circularisation in rotating discs: a semi-analytical description

Authors: Matteo Bonetti, Elisa Bortolas, Alessandro Lupi, Massimo Dotti, Sandra I. Raimundo

Abstract: We present and validate a novel semi-analytical approach to study the effect of dynamical friction on the orbits of massive perturbers in rotating stellar discs. We find that dynamical friction efficiently circularises the orbit of co-rotating perturbers, while it constantly increases the eccentricity of counter-rotating ones until their angular momenta reverse, then once again promoting circulari… ▽ More We present and validate a novel semi-analytical approach to study the effect of dynamical friction on the orbits of massive perturbers in rotating stellar discs. We find that dynamical friction efficiently circularises the orbit of co-rotating perturbers, while it constantly increases the eccentricity of counter-rotating ones until their angular momenta reverse, then once again promoting circularisation. Such "drag toward circular corotation" could shape the distribution of orientations of kinematically decoupled cores in disc galaxies, naturally leading to the observed larger fraction of co-rotating cores. △ Less

Submitted 4 April, 2020; v1 submitted 11 February, 2020; originally announced February 2020.

Comments: 7 pages, 3 figures, accepted for publication in MNRAS

Improved gravitational radiation time-scales: significance for LISA and LIGO-Virgo sources

Authors: Lorenz Zwick, Pedro R. Capelo, Elisa Bortolas, Lucio Mayer, Pau Amaro-Seoane

Abstract: We present a revised version of Peters' (1964) time-scale for the gravitational-wave (GW) induced decay of two point masses. The new formula includes the effects of the first-order post-Newtonian perturbation and additionally provides a simple fit to account for the Newtonian self-consistent evolution of the eccentricity. The revised time-scale is found by multiplying Peters' estimate by two facto… ▽ More We present a revised version of Peters' (1964) time-scale for the gravitational-wave (GW) induced decay of two point masses. The new formula includes the effects of the first-order post-Newtonian perturbation and additionally provides a simple fit to account for the Newtonian self-consistent evolution of the eccentricity. The revised time-scale is found by multiplying Peters' estimate by two factors, $R(e_0)= 8^{1-\sqrt{1-e_0}}$ and $Q_{\rm f}(p_0) = \exp \left(2.5 (r_{\rm S}/p_0) \right)$, where $e_0$ and $p_0$ are the initial eccentricity and periapsis, respectively, and $r_{\rm S}$ the Schwarzschild radius of the system. Their use can correct errors of a factor of 1-10 that arise from using the original Peters' formula. We apply the revised time-scales to a set of typical sources for existing ground-based laser interferometers and for the future Laser Interferometer Space Antenna (LISA), at the onset of their GW driven decay. We argue that our more accurate model for the orbital evolution will affect current event- and detection-rate estimates for mergers of compact object binaries, with stronger deviations for eccentric LISA sources, such as extreme and intermediate mass-ratio inspirals. We propose the correction factors $R$ and $Q_{\rm f}$ as a simple prescription to quantify decay time-scales more accurately in future population synthesis models. We also suggest that the corrected time-scale may be used as a computationally efficient alternative to numerical integration in other applications that include the modelling of radiation reaction for eccentric sources. △ Less

Submitted 7 May, 2020; v1 submitted 14 November, 2019; originally announced November 2019.

Comments: Accepted for publication in MNRAS

Improved constraints from ultra-faint dwarf galaxies on primordial black holes as dark matter

Authors: Jakob Stegmann, Pedro R. Capelo, Elisa Bortolas, Lucio Mayer

Abstract: Soon after the recent first ever detection of gravitational waves from merging black holes it has been suggested that their origin is primordial. Appealingly, a sufficient number of primordial black holes (PBHs) could also partially or entirely constitute the dark matter (DM) in the Universe. However, recent studies on PBHs in ultra-faint dwarf galaxies (UFDGs) suggest that they would dynamically… ▽ More Soon after the recent first ever detection of gravitational waves from merging black holes it has been suggested that their origin is primordial. Appealingly, a sufficient number of primordial black holes (PBHs) could also partially or entirely constitute the dark matter (DM) in the Universe. However, recent studies on PBHs in ultra-faint dwarf galaxies (UFDGs) suggest that they would dynamically heat up the stellar component due to two-body relaxation processes. From the comparison with the observed stellar velocity dispersions and the stellar half-light radii it was claimed that only PBHs with masses $\lesssim10\,M_\odot$ can significantly contribute to the DM. In this work, we improve the latter constraints by considering the largest observational sample of UFDGs and by allowing the PBH masses to follow an extended (log-normal) distribution. By means of collisional Fokker-Planck simulations, we explore a wide parameter space of UFDGs containing PBHs. The analysis of the half-light radii and velocity dispersions resulting from the simulations leads to three general findings that exclude PBHs with masses $\sim\mathcal{O}(1$-$100)\,M_\odot$ from constituting all of the DM: (i) We identify a critical sub-sample of UFDGs that only allows for $\sim\mathcal{O}(1)\,M_\odot$ PBH masses; (ii) for any PBH mass, there is an UFDG in our sample that disfavours it; (iii) the spatial extensions of a majority of simulated UFDGs containing PBHs are too large to match the observed. △ Less

Submitted 20 January, 2020; v1 submitted 10 October, 2019; originally announced October 2019.

Can supernova kicks trigger EMRIs in the Galactic Centre?

Authors: Elisa Bortolas, Michela Mapelli

Abstract: One of the most promising gravitational wave (GW) sources detectable by the forthcoming LISA observatory are the so-called extreme-mass ratio inspirals (EMRIs), i.e. GW-driven inspirals of stellar-mass compact objects onto supermassive black holes (SMBHs). In this paper, we suggest that supernova (SN) kicks may trigger EMRIs in galactic nuclei by scattering newborn stellar black holes and neutron… ▽ More One of the most promising gravitational wave (GW) sources detectable by the forthcoming LISA observatory are the so-called extreme-mass ratio inspirals (EMRIs), i.e. GW-driven inspirals of stellar-mass compact objects onto supermassive black holes (SMBHs). In this paper, we suggest that supernova (SN) kicks may trigger EMRIs in galactic nuclei by scattering newborn stellar black holes and neutron stars on extremely eccentric orbits; as a consequence, the time-scale over which these compact objects are expected to inspiral onto the central SMBH via GW emission may become shorter than the time-scale for other orbital perturbations to occur. By applying this argument to the Galactic Centre, we show that the S-cluster and the clockwise disc are optimal regions for the generation of such events: one SN out of about 10 000 (100 000) occurring in the S-cluster (clockwise disc) is expected to induce an EMRI. If we assume that the natal kicks affecting stellar black holes are significantly slower than those experienced by neutron stars, we find that most SN-driven EMRIs involve neutron stars. We further estimate the time spanning from the SN to the final plunge onto the SMBH to be of the order of few Myr. Finally, we extrapolate the rate of SN-driven EMRIs per Milky Way to be up to 10^-8/yr, thus we expect that LISA will detect up to a few tens of SN-driven EMRIs every year. △ Less

Submitted 26 February, 2019; v1 submitted 12 February, 2019; originally announced February 2019.

Comments: 15 pages, 10 figures, 4 tables; MNRAS, in press. Changes: we corrected the e-mail address for correspondence, and we added a few references

MSC Class: 85A05

Star Cluster Disruption by a Supermassive Black Hole Binary

Authors: Elisa Bortolas, Michela Mapelli, Mario Spera

Abstract: Supermassive black hole binaries (BHBs) are expected to be one of the most powerful sources of low-frequency gravitational waves (GWs) for future space-borne detectors. Prior to the GW emission stage, BHBs evolving in gas-poor nuclei shrink primarily through the slingshot ejection of stars approaching the BHB from sufficiently close distances. Here we address the possibility that the BHB shrinking… ▽ More Supermassive black hole binaries (BHBs) are expected to be one of the most powerful sources of low-frequency gravitational waves (GWs) for future space-borne detectors. Prior to the GW emission stage, BHBs evolving in gas-poor nuclei shrink primarily through the slingshot ejection of stars approaching the BHB from sufficiently close distances. Here we address the possibility that the BHB shrinking rate is enhanced through the infall of a star cluster (SC) onto the BHB. We present the results of direct summation N-body simulations exploring different orbits for the SC infall, and we show that SCs reaching the BHB on non-zero angular momentum orbits (with eccentricity 0.75) fail to enhance the BHB hardening, while SCs approaching the BHB on radial orbits reduce the BHB separation by ~10% in less than 10 Myr, effectively shortening the BHB path towards GWs. △ Less

Submitted 17 May, 2018; originally announced May 2018.

Comments: 5 pages, 2 figures; to appear in Proceedings of Science; proceeding of the GRAvitational-waves Science & technology Symposium - GRASS2018, 1-2 March 2018, Palazzo Moroni, Padova (Italy)

MSC Class: 85A05

Star cluster disruption by a massive black hole binary

Authors: Elisa Bortolas, Michela Mapelli, Mario Spera

Abstract: Massive black hole binaries (BHBs) are expected to form as the result of galaxy mergers; they shrink via dynamical friction and stellar scatterings, until gravitational waves (GWs) bring them to the final coalescence. It has been argued that BHBs may stall at a parsec scale and never enter the GW stage if stars are not continuously supplied to the BHB loss cone. Here we perform several N-body expe… ▽ More Massive black hole binaries (BHBs) are expected to form as the result of galaxy mergers; they shrink via dynamical friction and stellar scatterings, until gravitational waves (GWs) bring them to the final coalescence. It has been argued that BHBs may stall at a parsec scale and never enter the GW stage if stars are not continuously supplied to the BHB loss cone. Here we perform several N-body experiments to study the effect of an 80,000 solar masses stellar cluster (SC) infalling on a parsec-scale BHB. We explore different orbital elements for the SC and we perform runs both with and without accounting for the influence of a rigid stellar cusp (modelled as a rigid Dehnen potential). We find that the semi-major axis of the BHB shrinks by more than 10 per cent if the SC is on a nearly radial orbit; the shrinking is more efficient when a Dehnen potential is included and the orbital plane of the SC coincides with that of the BHB. In contrast, if the SC orbit has non-zero angular momentum, only a few stars enter the BHB loss cone and the resulting BHB shrinking is negligible. Our results indicate that SC disruption might significantly contribute to the shrinking of a parsec-scale BHB only if the SC approaches the BHB on a nearly radial orbit. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Comments: 12 pages, 12 figures, 2 tables; accepted for publication in MNRAS

MSC Class: 85A05

The influence of Massive Black Hole Binaries on the Morphology of Merger Remnants

Authors: Elisa Bortolas, Alessia Gualandris, Massimo Dotti, Justin I. Read

Abstract: Massive black hole (MBH) binaries, formed as a result of galaxy mergers, are expected to harden by dynamical friction and three-body stellar scatterings, until emission of gravitational waves (GWs) leads to their final coalescence. According to recent simulations, MBH binaries can efficiently harden via stellar encounters only when the host geometry is triaxial, even if only modestly, as angular m… ▽ More Massive black hole (MBH) binaries, formed as a result of galaxy mergers, are expected to harden by dynamical friction and three-body stellar scatterings, until emission of gravitational waves (GWs) leads to their final coalescence. According to recent simulations, MBH binaries can efficiently harden via stellar encounters only when the host geometry is triaxial, even if only modestly, as angular momentum diffusion allows an efficient repopulation of the binary loss cone. In this paper, we carry out a suite of N-body simulations of equal-mass galaxy collisions, varying the initial orbits and density profiles for the merging galaxies and running simulations both with and without central MBHs. We find that the presence of an MBH binary in the remnant makes the system nearly oblate, aligned with the galaxy merger plane, within a radius enclosing 100 MBH masses. We never find binary hosts to be prolate on any scale. The decaying MBHs slightly enhance the tangential anisotropy in the centre of the remnant due to angular momentum injection and the slingshot ejection of stars on nearly radial orbits. This latter effect results in about 1% of the remnant stars being expelled from the galactic nucleus. Finally, we do not find any strong connection between the remnant morphology and the binary hardening rate, which depends only on the inner density slope of the remnant galaxy. Our results suggest that MBH binaries are able to coalesce within a few Gyr, even if the binary is found to partially erase the merger-induced triaxiality from the remnant. △ Less

Submitted 23 March, 2018; v1 submitted 12 October, 2017; originally announced October 2017.

Comments: 16 pages, 13 figures, 4 tables; accepted for publication in MNRAS

MSC Class: 85A05

Supernova Kicks and Dynamics of Compact Remnants in the Galactic Centre

Authors: E. Bortolas, M. Mapelli, M. Spera

Abstract: The Galactic Centre (GC) is a unique place to study the extreme dynamical processes occurring near a super-massive black hole (SMBH). Here we investigate the role of supernova (SN) explosions occurring in massive binary systems lying in a disc-like structure within the innermost parsec. We use a regularized algorithm to simulate 30,000 isolated three-body systems composed of a stellar binary orbit… ▽ More The Galactic Centre (GC) is a unique place to study the extreme dynamical processes occurring near a super-massive black hole (SMBH). Here we investigate the role of supernova (SN) explosions occurring in massive binary systems lying in a disc-like structure within the innermost parsec. We use a regularized algorithm to simulate 30,000 isolated three-body systems composed of a stellar binary orbiting the SMBH. We start the integration when the primary member undergoes a SN explosion, and we analyze the impact of SN kicks on the orbits of stars and compact remnants. We find that SN explosions scatter the lighter stars in the pair on completely different orbits, with higher eccentricity and inclination. In contrast, stellar-mass black holes (BHs) and massive stars retain memory of the orbit of their progenitor star. Our results suggest that SN kicks are not sufficient to eject BHs from the GC. We thus predict that all BHs that form in situ in the central parsec of our Galaxy remain in the GC, building up a cluster of dark remnants. In addition, the change of NS orbits induced by SNe may partially account for the observed dearth of NSs in the GC. About 40 per cent of remnants stay bound to the stellar companion after the kick; we expect up to 70 per cent of them might become X-ray binaries through Roche-lobe filling. Finally, the eccentricity of some light stars becomes > 0.7 as an effect of the SN kick, producing orbits similar to those of the G1 and G2 dusty objects. △ Less

Submitted 19 April, 2017; originally announced April 2017.

Comments: 12 pages, 9 figures, 4 tables; accepted for publication in MNRAS

MSC Class: 85A05

Collisionless loss-cone refilling: there is no final parsec problem

Authors: Alessia Gualandris, Justin I. Read, Walter Dehnen, Elisa Bortolas

Abstract: Coalescing massive black hole binaries, formed during galaxy mergers, are expected to be a primary source of low frequency gravitational waves. Yet in isolated gas-free spherical stellar systems, the hardening of the binary stalls at parsec-scale separations owing to the inefficiency of relaxation-driven loss-cone refilling. Repopulation via collisionless orbit diffusion in triaxial systems is mor… ▽ More Coalescing massive black hole binaries, formed during galaxy mergers, are expected to be a primary source of low frequency gravitational waves. Yet in isolated gas-free spherical stellar systems, the hardening of the binary stalls at parsec-scale separations owing to the inefficiency of relaxation-driven loss-cone refilling. Repopulation via collisionless orbit diffusion in triaxial systems is more efficient, but published simulation results are contradictory. While sustained hardening has been reported in simulations of galaxy mergers with $N \sim 10^6$ stars and in early simulations of rotating models, in isolated non-rotating triaxial models the hardening rate continues to fall with increasing N, a signature of spurious two-body relaxation. We present a novel approach for studying loss cone repopulation in galactic nuclei. Since loss cone repopulation in triaxial systems owes to orbit diffusion, it is a purely collisionless phenomenon and can be studied with an approximated force calculation technique, provided the force errors are well behaved and sufficiently small. We achieve this using an accurate fast multipole method and define a proxy for the hardening rate that depends only on stellar angular momenta. We find that the loss cone is efficiently replenished even in very mildly triaxial models (with axis ratios 1 : 0.9 : 0.8). Such triaxiality is unavoidable following galactic mergers and can drive binaries into the gravitational wave regime. We conclude that there is no final parsec problem. △ Less

Submitted 29 September, 2016; originally announced September 2016.

Comments: 11 pages, 10 figures, accepted for publication in MNRAS

Dynamics of supernova remnants in the Galactic Centre

Authors: Elisa Bortolas, Michela Mapelli, Mario Spera

Abstract: The Galactic centre (GC) is a unique place to study the extreme dynamical processes occurring near a super-massive black hole (SMBH). Here we simulate a large set of binaries orbiting the SMBH while the primary member undergoes a supernova (SN) explosion, in order to study the impact of SN kicks on the orbits of stars and dark remnants in the GC. We find that SN explosions are efficient in scatter… ▽ More The Galactic centre (GC) is a unique place to study the extreme dynamical processes occurring near a super-massive black hole (SMBH). Here we simulate a large set of binaries orbiting the SMBH while the primary member undergoes a supernova (SN) explosion, in order to study the impact of SN kicks on the orbits of stars and dark remnants in the GC. We find that SN explosions are efficient in scattering neutron stars and other light stars on new (mostly eccentric) orbits, while black holes (BHs) tend to retain memory of the orbit of their progenitor star. SN kicks are thus unable to eject BHs from the GC: a cusp of dark remnants may be lurking in the central parsec of our Galaxy. △ Less

Submitted 22 June, 2016; originally announced June 2016.

Comments: 2 pages, 2 figures; To appear in: "Cosmic-Lab: Star Clusters as Cosmic Laboratories for Astrophysics, Dynamics and Fundamental Physics", F.R. Ferraro & B. Lanzoni eds, Mem. SAIt, Vol 87

Brownian motion of massive black hole binaries and the final parsec problem

Authors: E. Bortolas, A. Gualandris, M. Dotti, M. Spera, M. Mapelli

Abstract: Massive black hole binaries (BHBs) are expected to be one of the most powerful sources of gravitational waves (GWs) in the frequency range of the pulsar timing array and of forthcoming space-borne detectors. They are believed to form in the final stages of galaxy mergers, and then harden by slingshot ejections of passing stars. However, evolution via the slingshot mechanism may be ineffective if t… ▽ More Massive black hole binaries (BHBs) are expected to be one of the most powerful sources of gravitational waves (GWs) in the frequency range of the pulsar timing array and of forthcoming space-borne detectors. They are believed to form in the final stages of galaxy mergers, and then harden by slingshot ejections of passing stars. However, evolution via the slingshot mechanism may be ineffective if the reservoir of interacting stars is not readily replenished, and the binary shrinking may come to a halt at roughly a parsec separation. Recent simulations suggest that the departure from spherical symmetry, naturally produced in merger remnants, leads to efficient loss cone refilling, preventing the binary from stalling. However, current N-body simulations able to accurately follow the evolution of BHBs are limited to very modest particle numbers. Brownian motion may artificially enhance the loss cone refilling rate in low-N simulations, where the binary encounters a larger population of stars due its random motion. Here we study the significance of Brownian motion of BHBs in merger remnants in the context of the final parsec problem. We simulate mergers with various particle numbers (from 8k to 1M) and with several density profiles. Moreover, we compare simulations where the BHB is fixed at the centre of the merger remnant with simulations where the BHB is free to random walk. We find that Brownian motion does not significantly affect the evolution of BHBs in simulations with particle numbers in excess of one million, and that the hardening measured in merger simulations is due to collisionless loss cone refilling. △ Less

Submitted 21 June, 2016; originally announced June 2016.

Comments: 9 pages, 6 figures, 3 tables in Monthly Notices of the Royal Astronomical Society 2016