Theoretical astrophysics and plasma physics at RPC

Abstract:

Supermassive black hole (SMBH) binary systems are an unavoidable outcome of galaxy mergers. Their dynamics encode valuable information about their formation and growth, the composition of their host galactic nuclei, the evolution of galaxies, and the nature of gravity. Many SMBH binaries with separations pc-kpc have been found, but closer (subparsec) binaries remain to be confirmed. Identifying these systems may elucidate how binaries evolve past the “final parsec” until gravitational radiation drives them to coalescence. Methods to discover and characterize SMBH binaries can shed light on these important questions and potentially open new multimessenger channels. Here we show that SMBH binaries in nonactive galactic nuclei can be identified and characterized by the gravitational lensing of individual bright stars, located behind them in the host galaxy. The rotation of “caustics”—regions where sources are hugely magnified due to the SMBH binary’s orbit and inspiral—leads to quasiperiodic lensing of starlight (QPLS). The extreme lensing magnification of individual bright stars produces a significant variation in the host galaxies’ luminosity; their lightcurve traces the orbit of the SMBH binary and its evolution, analogous to the waveforms recorded by gravitational-wave (GW) detectors. QPLS probes the population of sources observable by pulsar timing arrays and space detectors (LISA, TianQin), offering advance warning triggers for merging SMBHs for coincident or follow-up GW detections. SMBH population models predict 1–50 $[190–5000](n_{\star }/{\mathrm{pc}}^{-3})$ QPLS binaries with period less than 10[40] yr with comparable masses and redshift $z<0.3$ , where $n_{\star }$ is the stellar number density. Additionally, stellar and orbital motion will lead to frequent instances of single or double flares caused by SMBHBs with longer periods. This novel signature can be searched for in a wealth of existing and upcoming time-domain photometric data: identifying quasiperiodic variability in galactic lightcurves will reveal an ensemble of binary systems and illuminate outstanding questions around them.

Abstract:

The centre of the Galaxy harbours a supermassive black hole, Sgr A*, which is surrounded by a massive star cluster known as the S-cluster. The most extensively studied star in this cluster is the B-type main-sequence S2 star (also known as S0-2). These types of stars are commonly found in binary systems in the Galactic field, but observations do not seem to detect a companion to S2. This absence may be attributed to observational biases or to a dynamically hostile environment caused by phenomena such as tidal disruption or mergers. Using a N-body code with first-order post-Newtonian corrections, we investigate whether S2 can host a stellar or planetary companion. We perform 10⁵ simulations adopting uniform distributions for the orbital elements of the companion. Our results show that companions may exist for orbital periods shorter than 100 days, eccentricities below 0.8, and across the full range of mutual inclinations. The number of surviving companions increases with shorter orbital periods, lower eccentricities, and nearly coplanar orbits. We also find that the disruption mechanisms include mergers driven by Lidov–Kozai cycles and breakups that occur when the companion surpasses the Hill radius of its orbit. Finally, we find that the presence of a companion would alter S2’s astrometric signal by no more than 5 μas. Current radial-velocity detection limits constrain viable stellar binary configurations to approximately 4.4% of the simulated cases. Including astrometric limits reduces to 4.3%. Imposing an additional constraint that any companion must have a mass ≲2 M ⊙ (otherwise it would be visible) narrows the fraction of undetectable stellar binaries to just 3.0%.

Abstract:

Relativistic magnetic reconnection has been proposed as an important non-thermal particle acceleration (NTPA) mechanism that generates power-law spectra and high-energy emissions. Power-law particle spectra are in general characterised by three parameters: the power-law index, the high-energy cutoff and the low-energy cutoff (i.e. the injection energy). Particle injection into the non-thermal power law, despite also being a critical step in the NTPA chain, has received considerably less attention than the subsequent acceleration to high energies. Open questions on particle injection that are important for both physical understanding and astronomical observations include how the upstream magnetisation influences the injection energy and the contributions of the known injection mechanisms (i.e. direct acceleration by the reconnection electric field, Fermi kicks and pickup acceleration) to the injected particle population. Using fully kinetic particle-in-cell simulations, we uncover these relationships by systematically measuring the injection energy and calculating the contributions of each acceleration mechanism to the total injected particle population. We also present a theoretical model to explain these results. Additionally, we compare two- and three-dimensional simulations to assess the impact of the flux-rope kink and drift-kink instability on particle injection. We conclude with comparisons with previous work and outlook for future work.

Abstract:

Compact objects evolving in an astrophysical environment experience a gravitational drag force known as dynamical friction. We present a multipole-frequency decomposition to evaluate the orbit-averaged energy and angular momentum dissipation experienced by point masses on periodic orbits within a homogeneous, fluidlike background. Our focus is on eccentric Keplerian trajectories. Although our approach is currently restricted to linear response theory, it is fully consistent within that framework. We validate our theoretical expressions for the specific case of an ideal fluid, using semi-numerical simulations of the linear response acoustic wake. We demonstrate that, for a finite-time perturbation switched on at t=0, a steady dissipation state is reached after a time bounded by twice the sound crossing time of the apocenter distance. We apply our results to model the secular evolution of compact eccentric binaries in a gaseous medium, assuming low-density conditions where the orbital elements evolve adiabatically. For unequal-mass systems with moderate initial eccentricity, the late-time eccentricity growth is significantly delayed compared to the equal-mass case, due to the binary components becoming transonic at different times along their orbital trajectory. Our approach offers a computationally efficient alternative to full simulations of the linear response wake.