Theoretical astrophysics and plasma physics at RPC

Abstract:

Here we present details of an operator-split, implicit-explicit numerical scheme for the solution of the gyrokinetic-Poisson system of equations in the local limit. This scheme has been implemented in a new code called stella, which is capable of evolving electrostatic fluctuations with full kinetic electron effects and an arbitrary number of ion species in general magnetic geometry. We demonstrate the advantages of this mixed approach over a fully explicit treatment and provide linear and nonlinear benchmark comparisons for both axisymmetric and non-axisymmetric magnetic equilibria.

Abstract:

© 2018 IOP Publishing Ltd. Differential rotation is induced in tokamak plasmas when an underlying symmetry of the governing gyrokinetic-Maxwell system of equations is broken. One such symmetry-breaking mechanism is considered here: the turbulent acceleration of particles along the mean magnetic field. This effect, often referred to as the 'parallel nonlinearity', has been implemented in the δf gyrokinetic code stella and used to study the dependence of turbulent momentum transport on the plasma size and on the strength of the turbulence drive. For JET-like parameters with a wide range of driving temperature gradients, the momentum transport induced by the inclusion of turbulent acceleration is similar to or smaller than the ratio of the ion Larmor radius to the plasma minor radius. This low level of momentum transport is explained by demonstrating an additional symmetry that prohibits momentum transport when the turbulence is driven far above marginal stability.

Abstract:

In a recent paper, we determined the measurement accuracy of physical parameters for eccentric, precessing, non-spinning, inspiraling, stellar-mass black hole - black hole (BH-BH) binaries for the upcoming second-generation LIGO/VIRGO/KAGRA detector network at design sensitivity using the Fisher matrix method. Here we extend that study to a wide range of binary masses including neutron star - neutron star (NS-NS), NS-BH, and BH-BH binaries with BH masses up to $110 \, M_{\odot}$. The measurement error of eccentricity $e_{10 \,\rm Hz}$ at a gravitational-wave (GW) frequency of $10 \, {\rm Hz}$ is in the range $(10^{-4}-10^{-3}) \times (D_L/ 100\,\rm Mpc)$ for NS-NS, NS-BH, and BH-BH binaries at a luminosity distance of $D_L$ if $e_{10 \,\rm Hz} \gtrsim 0.1 $. For events with masses and distances similar to the detected 10 GW transients, we show that nonzero orbital eccentricities may be detected if $0.081 \lesssim e_{10 \,\rm Hz}$. Consequently, the LIGO/VIRGO/KAGRA detector network at design sensitivity will have the capability to distinguish between eccentric waveforms and circular waveforms. In comparison to circular inspirals, we find that the chirp mass measurement precision can improve by up to a factor of $\sim 20$ and $\sim 50-100$ for NS-NS and NS-BH binaries with BH companion masses $\lesssim 40 \, M_{\odot}$, respectively. The identification of eccentric sources may give information on their astrophysical origin; it would indicate merging binaries in triple or higher multiplicity systems or dynamically formed binaries in dense stellar systems such as globular clusters or galactic nuclei.

Abstract:

We use the second data releases of the European Space AgencyGaia astrometric survey and the high-resolution Galactic Archaeology with HERMES (GALAH) spectroscopic survey to analyse the structure of our Galaxy’s disc components. With GALAH, we separate the α-rich and α-poor discs (with respect to Fe), which are superposed in both position and velocity space, and examine their distributions in action space. We study the distribution of stars in the zVz phase plane, for both Vϕ and VR, and recover the remarkable ‘phase spiral’ discovered by Gaia. We identify the anticipated quadrupole signature in zVz of a tilted velocity ellipsoid for stars above and below the Galactic plane. By connecting our work with earlier studies, we show that the phase spiral is likely to extend well beyond the narrow solar neighbourhood cylinder in which it was found. The phase spiral is a signature of corrugated waves that propagate through the disc, and the associated non-equilibrium phase mixing. The radially asymmetric distribution of stars involved in the phase spiral reveals that the corrugation, which is mostly confined to the α-poor disc, grows in z-amplitude with increasing radius. We present new simulations of tidal disturbance of the Galactic disc by the Sagittarius (Sgr) dwarf. The effect on the zVz phase plane lasts ≳2Gyr⁠, but a subsequent disc crossing wipes out the coherent structure. We find that the phase spiral was excited ≲0.5Gyr ago by an object like Sgr with total mass ∼3 × 1010 M⊙ (stripped down from ∼5 × 1010 M⊙ when it first entered the halo) passing through the plane.