Theoretical astrophysics and plasma physics at RPC

Abstract:

A new algorithm for toroidal flow shear in a linearly implicit, local δf gyrokinetic code is described. Unlike the current approach followed by a number of codes, it treats flow shear continuously in time. In the linear gyrokinetic equation, time-dependences arising from the presence of flow shear are decomposed in such a way that they can be treated explicitly in time with no stringent constraint on the time step. Flow shear related time dependences in the nonlinear term are taken into account exactly, and time dependences in the quasineutrality equation are interpolated. Test cases validating the continuous-in-time implementation in the code GS2 are presented. Lastly, nonlinear gyrokinetic simulations of a JET discharge illustrate the differences observed in turbulent transport compared with the usual, discrete-in-time approach. The continuous-in-time approach is shown, in some cases, to produce fluxes that converge to a different value than with the discrete approach. The new approach can also lead to substantial computational savings by requiring radially narrower boxes. At fixed box size, the continuous implementation is only modestly slower than the previous, discrete approach.

Abstract:

We have used the local-δf gyrokinetic code GS2 to perform studies of the effect of flux-surface shaping on two highly-shaped, low- and high-β JT-60SA-relevant equilibria, including a successful benchmark with the GKV code. We find that for a high-performance plasma, i.e. one with high plasma beta and steep pressure gradients, the turbulent outwards radial fluxes may be reduced by minimizing the elongation. We explain the results as a competition between the local magnetic shear and finite-Larmor-radius (FLR) stabilization. Electromagnetic studies indicate that kinetic ballooning modes are stabilized by increased shaping due to an increased sensitivity to FLR effects, relative to the ion-temperature-gradient instability. Nevertheless, at high enough β, increased elongation degrades the local magnetic shear stabilization that enables access to the region of ballooning second-stability.

Abstract:

We extend the relativistic time-dependent thin-disc TDE model to describe non-thermal (2-10 keV) X-ray emission produced by the Compton up-scattering of thermal disc photons by a compact electron corona, developing analytical and numerical models of the evolving non-thermal X-ray light curves. In the simplest cases, these X-ray light curves follow power-law profiles in time. We suggest that TDE discs act in many respects as scaled-up versions of XRB discs, and that such discs should undergo state transitions into harder accretion states. XRB state transitions typically occur when the disc luminosity becomes roughly one per cent of its Eddington value. We show that if the same is true for TDE discs then this, in turn, implies that TDEs with non-thermal X-ray spectra should come preferentially from large-mass black holes. The characteristic hard-state transition mass is MHS ≃ 2 × 107M⊙. Hence, subpopulations of thermal and non-thermal X-ray TDEs should come from systematically different black hole masses. We demonstrate that the known populations of thermal and non-thermal X-ray TDEs do indeed come from different distributions of black hole masses. The null-hypothesis of identical black hole mass distributions is rejected by a two-sample Anderson-Darling test with a p-value <0.01. Finally, we present a model for the X-ray rebrightening of TDEs at late times as they transition into the hard state. These models of evolving TDE light curves are the first to join both thermal and non-thermal X-ray components in a unified scenario.