Theoretical astrophysics and plasma physics at RPC

Abstract:

We comprehensively model the X-ray luminosity emergent from time-dependent relativistic accretion discs, developing analytical models of the X-ray luminosity of thermal disc systems as a function of black hole mass M, disc mass Md, and disc α-parameter. The X-ray properties of these solutions will be directly relevant for understanding tidal disruption event (TDE) observations. We demonstrate an extremely strong suppression of thermal X-ray luminosity from large mass black holes, LX ∼ exp (− m7/6), where m is a dimensionless mass, roughly the black hole mass in unity of 106M⊙. This strong suppression results in upper observable black hole mass limits, which we demonstrate to be of order Mlim ≃ 3 × 107M⊙, above which thermal X-ray emission will not be observable. This upper observable black hole mass limit is a function of the remaining disc parameters, and the full dependence can be described analytically (equation 82). We demonstrate that the current population of observed X-ray TDEs is indeed consistent with an upper black hole mass limit of order M ∼ 107M⊙, consistent with our analysis.

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.