Theoretical astrophysics and plasma physics at RPC

Abstract:

Multiple space and time scales arise in plasma turbulence in magnetic confinement fusion devices because of the smallness of the square root of the electron-to-ion mass ratio ${\left({m}_{{\rm{e}}}/{m}_{{\rm{i}}}\right)}^{1/2}$ and the consequent disparity of the ion and electron thermal gyroradii and thermal speeds. Direct simulations of this turbulence that include both ion and electron space–time scales indicate that there can be significant interactions between the two scales. The extreme computational expense and complexity of these direct simulations motivates the desire for reduced treatment. By exploiting the scale-separation between ion scales (IS) and electron scales (ES), and expanding the gyrokinetic equations for the turbulence in ${\left({m}_{{\rm{e}}}/{m}_{{\rm{i}}}\right)}^{1/2}$, we derive such a reduced system of gyrokinetic equations that describes cross-scale interactions. The coupled gyrokinetic equations contain novel terms which provide candidate mechanisms for the observed cross-scale interaction. The ES turbulence experiences a modified drive due to gradients in the IS distribution function, and is advected by the IS $E\times B$ drift, which varies in the direction parallel to the magnetic field line. The largest possible cross-scale term in the IS equations is sub-dominant in our ${\left({m}_{{\rm{e}}}/{m}_{{\rm{i}}}\right)}^{1/2}$ expansion. Hence, in our model the IS turbulence evolves independently of the ES turbulence. To complete the scale-separated approach, we provide and justify a parallel boundary condition for the coupled gyrokinetic equations in axisymmetric equilibria based on the standard 'twist-and-shift' boundary condition. This approach allows one to simulate multi-scale turbulence using ES flux tubes nested within an IS flux tube.

Abstract:

The growing number of stellar-mass binary black hole mergers discovered by Advanced LIGO and Advanced Virgo are starting to constrain the binaries' origin and environment. However, we still lack sufficiently accurate modeling of binary formation channels to obtain strong constraints, or to identify sub-populations. One promising formation mechanism that could result in different black hole properties is binaries merging within the accretion disks of Active Galactic Nuclei (AGN). Here we show that the black holes' orbital alignment with the AGN disks preferentially selects heavier black holes. We carry out Monte Carlo simulations of orbital alignment with AGN disks, and find that AGNs harden the initial black hole mass function. Assuming an initial power law mass distribution $M_{\rm bh}^{-\beta}$, we find that the power law index changes by $\Delta \beta\sim1.3$, resulting in a more top-heavy population of merging black holes. This change is independent of the mass of, and accretion rate onto, the supermassive black hole in the center of the AGN. Our simulations predict an AGN-assisted merger rate of $\sim4$Gpc$^{-3}$yr$^{-1}$. With its hardened mass spectra, the AGN channel could be responsible for $10-50$% of gravitational-wave detections.

Abstract:

With the hundreds of merging binary black hole (BH) signals expected to be detected by LIGO/Virgo, LISA and other instruments in the next few years, the modeling of astrophysical channels that lead to the formation of compact-object binaries has become of fundamental importance. In this paper, we carry out a systematic statistical study of quadruple BHs consisting of two binaries in orbit around their center of mass, by means of high-precision direct $N$-body simulations including Post-Newtonian (PN) terms up to 2.5PN order. We found that most merging systems have high initial inclinations and the distributions peak at $\sim 90^\circ$ as for triples, but with a more prominent broad distribution tail. We show that BHs merging through this channel have a significant eccentricity in the LIGO band, typically much larger than BHs merging in isolated binaries and in binaries ejected from star clusters, but comparable to that of merging binaries formed via the GW capture scenario in clusters, mergers in hierarchical triples, or BH binaries orbiting intermediate-mass black holes in star clusters. We show that the merger fraction can be up to $\sim 3$--$4\times$ higher for quadruples than for triples. Thus even if the number of quadruples is $20\%$--$25\%$ of the number of triples, the quadruple scenario can represent an important contribution to the events observed by LIGO/VIRGO.

Abstract:

Stellar-mass black hole binaries (BHBs) near supermassive black holes (SMBH) in galactic nuclei undergo eccentricity oscillations due to gravitational perturbations from the SMBH. Previous works have shown that this channel can contribute to the overall BHB merger rate detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo Interferometer. Significantly, the SMBH gravitational perturbations on the binary's orbit may produce eccentric BHBs which are expected to be visible using the upcoming Laser Interferometer Space Antenna (LISA) for a large fraction of their lifetime before they merge in the LIGO/Virgo band. For a proof-of-concept, we show that the eccentricity oscillations of these binaries can be detected with LISA for BHBs in the local universe up to a few Mpcs, with observation periods shorter than the mission lifetime, thereby disentangling this merger channel from others. The approach presented here is straightforward to apply to a wide variety of compact object binaries with a tertiary companion.