Explore the vast range of titles available.

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.

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.

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.

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\\mu) < -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.

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 \\(\\rho\\lesssim 0.3\\) GeV/cm\\(^3\\).

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.

Ground-based gravitational-wave (GW) observatories have transformed our view of compact-object mergers, yet their reach still limits a comprehensive reconstruction of the processes that generate these systems. Only next-generation observatories, with order-of-magnitude improvements in sensitivity and access to lower frequencies, will be capable of radically extending this detection horizon. GW observations will make it possible to detect the complete population of binary black hole (BBH) mergers out to redshifts of \\(z \\simeq 100\\). This capability will deliver an unprecedented map of merger events across cosmic time and enable precise reconstruction of their mass and spin distributions, while for several thousand events the signal-to-noise ratio will surpass 100, enabling precision physics of BHs and neutron stars (NSs). The access to lower frequencies will also open the intermediate-mass window, detecting systems of order \\(\\sim 10^3 M_\\odot\\), potentially in coordination with multi-band observations from LISA. At higher redshifts, where Population III stars have so far remained beyond reach - even for the James Webb Space Telescope - GW observations by next-generation detectors will routinely provide observations of BH mergers thought to originate from these primordial stellar populations. Such measurements are expected to play a central role in clarifying the early assembly of supermassive black holes. A single detection of a binary BH system at \\(z \\gtrsim 30\\), or of a compact object with sub-solar mass and no tidal deformability, would constitute strong evidence for the existence of primordial black holes. Such a discovery would have profound consequences for our understanding of dark matter and the early Universe. Ultimately, the GW observations will become revolutionary for identifying the physical channels responsible for compact binary formation.

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.

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.

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\\sigma\\) 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.