Theoretical astrophysics and plasma physics at RPC

Abstract:

The evidence of a large diffusive turbulent contribution to the radial impurity transport in Wendelstein 7-X (W7-X) plasmas has been experimentally inferred during the first campaigns and numerically confirmed by means of gyrokinetic simulations with the code stella. In general, the absence of strong impurity accumulation during the initial W7-X campaigns is attributed to this diffusive term. Given the large variety of possible stellarator configurations, in the present work the diffusive contribution is also calculated in other stellarator plasmas. In particular, a numerical cross-device comparison is presented, where the diffusion (D) and convection (V) coefficients of carbon and iron impurities produced by ion-temperature-gradient (ITG) turbulence are obtained. The simulations have been performed for the helias W7-X, the heliotron LHD, the heliac TJ-II and the quasi-axisymmetric stellarator NCSX at the radial position r/a = 0.75. The results show that, although the size of D and V can differ across the four devices, inward convection is found for all of them. For W7-X, TJ-II and NCSX the two coefficients are comparable and the turbulent peaking factor is surprisingly similar. In LHD, appreciably weaker diffusive and convective impurity transport and significantly larger turbulent peaking factor, in comparison with the other three stellarators, are predicted. All this suggests that ITG turbulence, although not strongly, would lead to negative impurity density gradients in stellarators. Then, considering mixed ITG/trapped electron mode (TEM) turbulence for the specific case of W7-X, it has been quantitatively assessed to what degree pellet fueled reduced turbulence scenarios feature reduced turbulent transport of impurities as well. The results for trace iron impurities show that, although their turbulent transport is not entirely suppressed, a significant reduction of the convection and a stronger decrease of the diffusion term are found. Although the diffusion is still above neoclassical levels, the neoclassical convection would gain under such conditions a greater specific weight on the dynamics of impurities in comparison with standard ECRH scenarios without pellet fueling.

Abstract:

We use the problem of dynamical friction within the periodic cube to illustrate the application of perturbation theory in stellar dynamics, testing its predictions against measurements from N-body simulations. Our development is based on the explicitly time-dependent Volterra integral equation for the cube’s linear response, which avoids the subtleties encountered in analyses based on complex frequency. We obtain an expression for the self-consistent response of the cube to steady stirring by an external perturber. From this, we show how to obtain the familiar Chandrasekhar dynamical friction formula and construct an elementary derivation of the Lenard–Balescu equation for the secular quasi-linear evolution of an isolated cube composed of N equal-mass stars. We present an alternative expression for the (real-frequency) van Kampen modes of the cube and show explicitly how to decompose any linear perturbation of the cube into a superposition of such modes.

Abstract:

Recent gravitational wave (GW) observations by LIGO/Virgo show evidence for hierarchical mergers, where the merging BHs are the remnants of previous BH merger events. These events may carry important clues about the astrophysical host environments of the GW sources. In this paper, we present the distributions of the effective spin parameter (χeff), the precession spin parameter (χp), and the chirp mass (mchirp) expected in hierarchical mergers. Under a wide range of assumptions, hierarchical mergers produce (i) a monotonic increase of the average of the typical total spin for merging binaries, which we characterize with χ¯typ≡(χ2eff+χ2p)1/2⁠, up to roughly the maximum mchirp among first-generation (1g) BHs, and (ii) a plateau at χ¯typ∼0.6 at higher mchirp. We suggest that the maximum mass and typical spin magnitudes for 1g BHs can be estimated from χ¯typ as a function of mchirp. The GW data observed in LIGO/Virgo O1–O3a prefers an increase in χ¯typ at low mchirp, which is consistent with the growth of the BH spin magnitude by hierarchical mergers at ∼2σ confidence. A Bayesian analysis using the χeff, χp, and mchirp distributions suggests that 1g BHs have the maximum mass of ∼15–30M⊙ if the majority of mergers are of high-generation BHs (not among 1g–1g BHs), which is consistent with mergers in active galactic nucleus discs and/or nuclear star clusters, while if mergers mainly originate from globular clusters, 1g BHs are favoured to have non-zero spin magnitudes of ∼0.3. We also forecast that signatures for hierarchical mergers in the χ¯typ distribution can be confidently recovered once the number of GW events increases to ≳ O(100).