Theoretical astrophysics and plasma physics at RPC

Abstract:

ABSTRACT A theory of cosmic-ray transport in multiphase diffusive media is developed, with the specific application to cases in which the cosmic-ray diffusion coefficient has large spatial fluctuations that may be inherently multiscale. We demonstrate that the resulting transport of cosmic rays is diffusive in the long-time limit, with an average diffusion coefficient equal to the harmonic mean of the spatially varying diffusion coefficient. Thus, cosmic-ray transport is dominated by areas of low diffusion even if these areas occupy a relatively small, but not infinitesimal, fraction of the volume. On intermediate time-scales, the cosmic rays experience transient effective subdiffusion, as a result of low-diffusion regions interrupting long flights through high-diffusion regions. In the simplified case of a two-phase medium, we show that the extent and extremity of the subdiffusivity of cosmic-ray transport is controlled by the spectral exponent of the distribution of patch sizes of each of the phases. We finally show that, despite strongly influencing the confinement times, the multiphase medium is only capable of altering the energy dependence of cosmic-ray transport when there is a moderate (but not excessive) level of perpendicular diffusion across magnetic-field lines.

Abstract:

We introduce a new approach to measure the magnetic pitch angle profile in tokamak plasmas with Doppler backscattering (DBS), a technique traditionally used for measuring flows and density fluctuations. The DBS signal is maximised when its probe beam’s wavevector is perpendicular to the magnetic field at the cutoff location, independent of the density fluctuations (Hillesheim et al 2015 Nucl. Fusion 55 073024). Hence, if one could isolate this effect, DBS would then yield information about the magnetic pitch angle. By varying the toroidal launch angle, the DBS beam reaches cutoff with different angles with respect to the magnetic field, but with other properties remaining similar. Hence, the toroidal launch angle which gives maximum backscattered power is thus that which is matched to the pitch angle at the cutoff location, enabling inference of the magnetic pitch angle. We performed systematic scans of the DBS toroidal launch angle for repeated DIII-D tokamak discharges. Experimental DBS data from this scan were analysed and combined with Gaussian beam-tracing simulations using the Scotty code (Hall-Chen et al 2022 Plasma Phys. Control. Fusion 64 095002). The pitch-angle inferred from DBS is consistent with that from magnetics-only and motional-Stark-effect-constrained (MSE) equilibrium reconstruction in the edge. In the core, the pitch angles from DBS and magnetics-only reconstructions differ by one to two degrees, while simultaneous MSE measurements were not available. The uncertainty in these measurements was under a degree; we show that this uncertainty is primarily due to the error in toroidal steering, the number of toroidally separated measurements, and shot-to-shot repeatability. We find that the error of pitch-angle measurements can be reduced by optimising the poloidal launch angle and initial beam properties. Since DBS has high spatial and temporal resolutions, is non-perturbative, does not require neutral beams, and is likely robust to neutron damage of and debris on the first mirrors, using DBS to measure the pitch angle in future fusion energy systems is especially appealing.

Abstract:

The generation of dense electron-positron pair beams in the laboratory can enable direct tests of theoretical models of γ-ray bursts and active galactic nuclei. We have successfully achieved this using ultrarelativistic protons accelerated by the Super Proton Synchrotron at (CERN). In the first application of this experimental platform, the stability of the pair beam is studied as it propagates through a meter-length plasma, analogous to TeV γ-ray-induced pair cascades in the intergalactic medium. It has been argued that pair beam instabilities disrupt the cascade, thus accounting for the observed lack of reprocessed GeV emission from TeV blazars. If true, this would remove the need for a moderate strength intergalactic magnetic field to explain the observations. We find that the pair beam instability is suppressed if the beam is not perfectly collimated or monochromatic, hence the lower limit to the intergalactic magnetic field inferred from γ-ray observations of blazars is robust.

Abstract:

Abstract Merging black hole (BH) binaries in AGN discs formed through two-body scatterings via the “gas-capture” process may explain a significant fraction of BH mergers in AGN and a non-negligible contribution to the observed rate from LIGO-VIRGO-KAGRA. We perform Monte Carlo simulations of binary BH formation, evolution and mergers across the observed AGN mass function using a novel physically motivated treatment for the gas-capture process derived from hydrodynamical simulations of BH-BH encounters in AGN. Our models suggest that gas-captured binaries could result in merger rates of 0.73 − 7.1Gpc−3yr−1. Mergers from AGN are dominated by AGN with supermassive BH masses of ∼107M⊙, with 90 % of mergers occurring in the range ∼106M⊙ − 108M⊙. The merging mass distribution is flatter than the initial BH mass power law by a factor Δξ = 1.1 − 1.2, as larger BHs align with the disc and form binaries more efficiently. Similarly, the merging mass ratio distribution is flatter, therefore the AGN channel could explain high mass and unequal mass ratio detections such as GW190521 and GW190814. Using a simpler dynamical friction treatment for the binary formation process, the results are similar, where the primary bottleneck is the alignment time with the disc. The most influential parameters are the anticipated number of BHs and their mass function. Given the many uncertainties that remain in the AGN channel, we expect the true uncertainty extends beyond our predicted rates. Nonetheless, we conclude that AGN remain an important channel for consideration, particularly for gravitational wave detections involving one or two high mass BHs.

Abstract:

Gaia observations have revealed over a million stellar binary candidates within ∼1 kpc of the Sun, predominantly characterized by orbital separations >103 au and eccentricities >0.7. The prevalence of such wide, eccentric binaries has proven challenging to explain through canonical binary formation channels. However, recent advances in our understanding of three-body binary formation (3BBF)—new binary assembly by the gravitational scattering of three unbound bodies (3UB)—have shown that 3BBF in star clusters can efficiently generate wide, highly eccentric binaries. We further explore this possibility by constructing a semi-analytic model of the Galactic binary population in the solar neighborhood, originating from 3BBF in star clusters and subsequently migrating to the solar neighborhood within a Hubble time. The model relies on 3BBF scattering experiments to determine how the 3BBF rate and resulting binary properties scale with local stellar density, velocity dispersion, and physically motivated limits to 3UB encounters within a clusters’ tidal field. The Galactic star cluster population is modeled by incorporating up-to-date prescriptions for the Galaxy’s star formation history as well as the birth properties and internal evolution of its star clusters. Finally, we account for binary disruption induced by perturbations from stellar interactions before cluster dissolution and the subsequent changes and disruption of binary orbital elements induced by dynamical interactions in the Galactic field. Without any explicit fine-tuning, our model closely reproduces the total number of Gaia’s wide binaries and the separation and eccentricity distributions, suggesting that 3BBF may be an important formation channel for these enigmatic systems.