Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThree-dimensional, global, radiative GRMHD simulations of a thermally unstable disc
188
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present results of a set of three-dimensional, general relativistic radiation magnetohydrodynamics simulations of thin accretion discs around a non-rotating black hole to test their thermal stability. We consider two cases, one that is initially radiation pressure dominated and expected to be thermally unstable and another that is initially gas-pressure dominated and expected to remain stable. Indeed, we find that cooling dominates over heating in the radiation pressure dominated model, causing the disc to collapse vertically on roughly the local cooling timescale. We also find that heating and cooling within the disc have a different dependence on the mid-plane pressure, a prerequisite of thermal instability. Comparison of our data with the relevant thin-disc thermal equilibrium curve suggests that our disc may be headed for the thermally stable, gas-pressure-dominated branch. However, because the disc collapses to the point that we are no longer able to resolve it, we had to terminate the simulation. On the other hand, the gas pressure dominated model, which was run for twice as long as the radiation pressure dominated one, remains stable, with heating and cooling roughly in balance. Finally, the radiation pressure dominated simulation shows some evidence of viscous instability. The strongest evidence is in plots of surface density, which show the disc breaking up into rings.
rate research
Read More
One of the main dissipation processes acting on all scales in relativistic jets is thought to be governed by magnetic reconnection. Such dissipation processes have been studied in idealized environments, such as reconnection layers, which evolve in merging islands and lead to the production of plasmoids, ultimately resulting in efficient particle acceleration. In accretion flows onto black holes, reconnection layers can be developed and destroyed rapidly during the turbulent evolution of the flow. We present a series of two-dimensional general-relativistic magnetohydrodynamic simulations of tori accreting onto rotating black holes focusing our attention on the formation and evolution of current sheets. Initially, the tori are endowed with a poloidal magnetic field having a multi-loop structure along the radial direction and with an alternating polarity. During reconnection processes, plasmoids and plasmoid chains are developed leading to a flaring activity and hence to a variable electromagnetic luminosity. We describe the methods developed to track automatically the plasmoids that are generated and ejected during the simulation, contrasting the behaviour of multi-loop initial data with that encountered in typical simulations of accreting black holes having initial dipolar field composed of one loop only. Finally, we discuss the implications that our results have on the variability to be expected in accreting supermassive black holes.
The candidate supermassive black hole in the Galactic Centre, Sagittarius A* (Sgr A*), is known to be fed by a radiatively inefficient accretion flow (RIAF), inferred by its low accretion rate. Consequently, radiative cooling has in general been overlooked in the study of Sgr A*. However, the radiative properties of the plasma in RIAFs are poorly understood. In this work, using full 3D general-relativistic magneto-hydrodynamical simulations, we study the impact of radiative cooling on the dynamical evolution of the accreting plasma, presenting spectral energy distributions and synthetic sub-millimeter images generated from the accretion flow around Sgr A*. These simulations solve the approximated equations for radiative cooling processes self-consistently, including synchrotron, bremsstrahlung, and inverse Compton processes. We find that radiative cooling plays an increasingly important role in the dynamics of the accretion flow as the accretion rate increases: the mid-plane density grows and the infalling gas is less turbulent as cooling becomes stronger. The changes in the dynamical evolution become important when the accretion rate is larger than $10^{-8},M_{odot}~{rm yr}^{-1}$ ($gtrsim 10^{-7} dot{M}_{rm Edd}$, where $dot{M}_{rm Edd}$ is the Eddington accretion rate). The resulting spectra in the cooled models also differ from those in the non-cooled models: the overall flux, including the peak values at the sub-mm and the far-UV, is slightly lower as a consequence of a decrease in the electron temperature. Our results suggest that radiative cooling should be carefully taken into account in modelling Sgr A* and other low-luminosity active galactic nuclei that have a mass accretion rate of $dot{M} > 10^{-7},dot{M}_{rm Edd}$.
Luminous active galactic nuclei (AGN) and X-Ray binaries (XRBs) tend to be surrounded by geometrically thin, radiatively cooled accretion discs. According to both theory and observations, these are -- in many cases -- highly misaligned with the black hole spin axis. In this work we present the first general relativistic magnetohydrodynamic simulations of very thin ($h/r sim 0.015-0.05$) accretion discs around rapidly spinning ($a sim 0.9$) black holes and tilted by 45-65 degrees. We show that the inner regions of the discs with $h/r lesssim 0.03$ align with the black hole equator, though at smaller radii than predicted by theoretical work. The inner aligned and outer misaligned disc regions are separated by a sharp break in tilt angle accompanied by a sharp drop in density. We find that frame-dragging by the spinning black hole overpowers the disc viscosity, which is self-consistently produced by magnetized turbulence, tearing the disc apart and forming a rapidly precessing inner sub-disc surrounded by a slowly precessing outer sub-disc. We find that at all tilt values the system produces a pair of relativistic jets. At small distances the jets precess rapidly together with the inner sub-disc, whereas at large distances they partially align with the outer sub-disc and precess more slowly. If the tearing radius can be modeled accurately in future work, emission model independent measurements of black hole spin based on precession-driven quasi-periodic oscillations may become possible.
The structure and evolution of protoplanetary disks (PPDs) are largely governed by disk angular momentum transport, mediated by magnetic fields. In the most observable outer disk, PPD gas dynamics is primarily controlled by ambipolar diffusion as the dominant non-ideal magnetohydrodynamic (MHD) effect. In this work, we study the gas dynamics in outer PPDs by conducting a set of global 3D non-ideal MHD simulations with ambipolar diffusion and net poloidal magnetic flux, using the Athena++ MHD code, with resolution comparable to local simulations. Our simulations demonstrate the co-existence of magnetized disk wind and turbulence driven by the magneto-rotational instability (MRI). While MHD wind dominates disk angular momentum transport, the MRI turbulence also contributes significantly. We observe that magnetic flux spontaneously concentrate into axisymmetric flux sheets, leading to radial variations in turbulence levels, stresses, and accretion rates. Annular substructures arise as a natural consequence of magnetic flux concentration. The flux concentration phenomena show diverse properties with different levels of disk magnetization and ambipolar diffusion. The disk generally loses magnetic flux over time, though flux sheets could prevent the leak of magnetic flux in some cases. Our results demonstrate the ubiquity of disk annular substructures in weakly MRI turbulent outer PPDs, and imply a stochastic nature of disk evolution.
Essentially all low-mass X-ray binaries (LMXBs) in the soft state appear to drive powerful equatorial disc winds. A simple mechanism for driving such outflows involves X-ray heating of the top of the disc atmosphere to the Compton temperature. Beyond the Compton radius, the thermal speed exceeds the escape velocity, and mass loss is inevitable. Here, we present the first coupled radiation-hydrodynamic simulation of such thermally-driven disc winds. The main advance over previous modelling efforts is that the frequency-dependent attenuation of the irradiating SED is taken into account. We can therefore relax the approximation that the wind is optically thin throughout which is unlikely to hold in the crucial acceleration zone of the flow. The main remaining limitations of our simulations are connected to our treatment of optically thick regions. Adopting parameters representative of the wind-driving LMXB GRO~J1655-40, our radiation-hydrodynamic model yields a mass-loss rate that is $simeq5times$ lower than that suggested by pure hydrodynamic, optically thin models. This outflow rate still represents more than twice the accretion rate and agrees well with the mass-loss rate inferred from Chandra/HETG observations of GRO~J1655-40 at a time when the system had a similar luminosity to that adopted in our simulations. The Fe XXV and Fe XXVI Lyman $rm{alpha}~$ absorption line profiles observed in this state are slightly stronger than those predicted by our simulations but the qualitative agreement between observed and simulated outflow properties means that thermal driving is a viable mechanism for powering the disc winds seen in soft-state LMXBs.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
