Simulations Of Accretion Disks Research Articles

Exploring Viscosity Effects in Accretion Disk Simulations: A Code Comparison between FLASH and Athena

This project investigates the impact of viscosity on the evolution of black hole accretion disks through a comparative analysis between FLASH and Athena simulations. Our focus lies in understanding the role of viscosity in overall disk evolution. We present mass accretion rate plots that reveal near identical patterns in disk evolution when varying viscosity. We emphasize the need for additional code-comparison research, including exploration of a wider variety of codes and physics included to enhance our understanding and confidence in simulated models.

Inner dusty regions of protoplanetary discs - III. The role of non-radial radiation pressure in dust dynamics

Abstract We explore dynamical behaviour of dust particles that populate the surface of inner optically thick protoplanetary discs. This is a disc region with the hottest dust and of a great importance for planet formation and dust evolution, but we still struggle to understand all the forces that shape this environment. In our approach we combine results from two separate numerical studies - one is the wind velocity and density distributions obtained from magnetohydrodynamical simulations of accretion discs, and the other is a high-resolution multigrain dust radiation transfer. In our previous paper in the series, we described the methodology for utilising these results as an environmental input for the integration of dust trajectories driven by gravity, gas drag, and radiation pressure. Now we have two improvements - we incorporate time changes in the wind density and velocity, and we implement the non-radial radiation pressure force. We applied our analysis on Herbig Ae and T Tau stars. We confirm that the radiation pressure force can lead to dust outflow, especially in the case of more luminous stars. Additionally, it opposes dust accretion at the inner disc edge and reduces dust settling. These effects are enhanced by the disc wind, especially in the zone where the stellar and the disc magnetic fields meet. Our results suggest that dust grains can stay in the hottest disc region for an extended period and then can end up ejected into the outer disc regions.

Magnetorotational dynamo can generate large-scale vertical magnetic fields in 3D GRMHD simulations of accreting black holes

ABSTRACT Jetted astrophysical phenomena with black hole engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed three-dimensional global general relativistic magnetohydrodynamical simulations of accretion discs to quantify, for the first time, a dynamo mechanism that generates large-scale magnetic fields. This dynamo mechanism primarily arises from the non-linear evolution of the magnetorotational instability (MRI). In this mechanism, large non-axisymmetric MRI-amplified shearing wave modes, mediated by the axisymmetric azimuthal magnetic field, generate and sustain the large-scale vertical magnetic field through their non-linear interactions. We identify the advection of magnetic loops as a crucial feature, transporting the large-scale vertical magnetic field from the outer regions to the inner regions of the accretion disc. This leads to a larger characteristic size of the, now advected, magnetic field when compared to the local disc height. We characterize the complete dynamo mechanism with two time-scales: one for the local magnetic field generation, $t_{\rm gen}$, and one for the large-scale scale advection, $t_{\rm adv}$. Whereas the dynamo we describe is non-linear, we explore the potential of linear mean field models to replicate its core features. Our findings indicate that traditional $\alpha$-dynamo models, often computed in stratified shearing box simulations, are inadequate and that the effective large-scale dynamics is better described by the shear current effects or stochastic $\alpha$-dynamos.

Introducing two improved methods for approximating radiative cooling in hydrodynamical simulations of accretion discs

ABSTRACT The evolution of many astrophysical systems depends strongly on the balance between heating and cooling, in particular star formation in giant molecular clouds and the evolution of young protostellar systems. Protostellar discs are susceptible to the gravitational instability, which can play a key role in their evolution and in planet formation. The strength of the instability depends on the rate at which the system loses thermal energy. To study the evolution of these systems, we require radiative cooling approximations because full radiative transfer is generally too expensive to be coupled to hydrodynamical models. Here, we present two new approximate methods for computing radiative cooling that make use of the polytropic cooling approximation. This approach invokes the assumption that each parcel of gas is located within a spherical pseudo-cloud, which can then be used to approximate the optical depth. The first method combines the methods introduced by Stamatellos et al. and Lombardi et al. to overcome the limitations of each method at low and high optical depths, respectively. The second method, the ‘modified Lombardi’ method, is specifically tailored for self-gravitating discs. This modifies the scale height estimate from the method of Lombardi et al. using the analytical scale height for a self-gravitating disc. We show that the modified Lombardi method provides an excellent approximation for the column density in a fragmenting disc, a regime in which the existing methods fail to recover the clumps and spiral structures. We therefore recommend this improved radiative cooling method for more realistic simulations of self-gravitating discs.

Plasmoid identification and statistics in two-dimensional Harris sheet and GRMHD simulations

ABSTRACT Magnetic reconnection is a ubiquitous phenomenon for magnetized plasma and leads to the rapid reconfiguration of magnetic field lines. During reconnection events, plasma is heated and accelerated until the magnetic field lines enclose and capture the plasma within a circular configuration. These so-called plasmoids could therefore observationally manifest themselves as hotspots, which are associated with flaring behaviour in supermassive black hole systems, such as Sagittarius A*. We have developed a novel algorithm for identifying plasmoid structures, which incorporates watershed and custom closed contouring steps. From the identified structures, we determine the plasmoids’ plasma characteristics and energetics in magnetohydrodynamical simulations. The algorithm’s performance is showcased for a high-resolution suite of axisymmetric ideal and resistive magnetohydrodynamical simulations of turbulent accretion discs surrounding a supermassive black hole. For validation purposes, we also evaluate several Harris current sheets that are well-investigated in the literature. We recover the characteristic power-law distribution of plasmoid sizes for both the black hole and Harris sheet simulations. This indicates that while the dynamics are vastly different, with different dominant plasma instabilities, the plasmoid creation characteristics are similar. Plasmoid formation rates for resistive general relativistic magnetohydrodynamical simulations are significantly higher than for their ideal counterpart. Moreover, the largest identified plasmoids are consistent with sizes typically assumed for semi-analytical interpretation of observations. We recover a positive correlation between the plasmoid formation rate and decreases in black-hole-horizon-penetrating magnetic flux, during which the accretion flow is temporarily halted. These results demonstrate the efficacy of the newly developed algorithm which has enabled an extensive quantitative analysis of plasmoid formation in black hole accretion simulations.

Emergent Nucleosynthesis from a 1.2 s Long Simulation of a Black Hole Accretion Disk

We simulate a black hole accretion disk system with full-transport general relativistic neutrino radiation magnetohydrodynamics for 1.2 s. This system is likely to form after the merger of two compact objects and is thought to be a robust site of r-process nucleosynthesis. We consider the case of a black hole accretion disk arising from the merger of two neutron stars. Our simulation time coincides with the nucleosynthesis timescale of the r-process (∼1 s). Because these simulations are time-consuming, it is common practice to run for a “short” duration of approximately 0.1–0.3 s. We analyze the nucleosynthetic outflow from this system and compare the results of stopping at 0.12 and 1.2 s. We find that the addition of mass ejected in the longer simulation as well as more favorable thermodynamic conditions from emergent viscous ejecta greatly impacts the nucleosynthetic outcome. We quantify the error in nucleosynthetic outcomes between short and long cuts.

Evolution of the Magnetic Field in High- and Low-β Disks with Initially Toroidal Fields

Abstract We present the results from a pair of high-resolution, long-timescale (∼105 GM/c 3), global, three-dimensional magnetohydrodynamical accretion disk simulations with differing initial magnetic plasma β in order to study the effects of the initial toroidal field strength on the production of a large-scale poloidal field. We initialize our disks in approximate equilibrium with purely toroidal magnetic fields of strength β 0 = 5 and 200. We also perform a limited resolution study. We find that simulations of differing field strengths diverge early in their evolution and remain distinct over the time studied, indicating that the initial magnetic conditions leave a persistent imprint in our simulations. Neither simulation enters the magnetically arrested disk regime. Both simulations are able to produce poloidal fields from initially toroidal fields, with the β 0 = 5 simulation evolving clear signs of a large-scale poloidal field. We make a cautionary note that computational artifacts in the form of large-scale vortices may be introduced in the combination of initially weak field and disk-internal mesh refinement boundaries, as evidenced by the production of an m = 1 mode overdensity in the weak field simulation. Our results demonstrate that the initial toroidal field strength plays a vital role in the simulated disk evolution for the models studied.

General relativistic magnetohydrodynamic simulations of accretion disks around tilted binary black holes of unequal mass

General relativistic magnetohydrodynamic simulations of accretion disks around tilted binary black holes of unequal mass

Multifrequency General Relativistic Radiation Magnetohydrodynamic Simulations of Thin Disks

We present a set of six general relativistic, multifrequency, radiation magnetohydrodynamic simulations of thin accretion disks with different target mass accretion rates around black holes with spins ranging from nonrotating to rapidly spinning. The simulations use the M 1 closure scheme with 12 independent frequency (or energy) bins ranging logarithmically from 5 × 10−3 keV to 5 × 103 keV. The multifrequency capability allows us to generate crude spectra and energy-dependent light curves directly from the simulations without a need for special postprocessing. While we generally find roughly thermal spectra with peaks around 1–4 keV, our high-spin cases showed harder-than-expected tails for the soft or thermally dominant state. This leads to radiative efficiencies that are up to five times higher than expected for a Novikov–Thorne disk at the same spin. We attribute these high efficiencies to the high-energy, coronal emission. These coronae mostly occupy the effectively optically thin regions near the inner edges of the disks and also cover or sandwich the inner ∼15GM/c 2 of the disks.

Chaotic Type I migration in turbulent discs

ABSTRACT By performing global hydrodynamical simulations of accretion discs with driven turbulence models, we demonstrate that elevated levels of turbulence induce highly stochastic migration torques on low-mass companions embedded in these discs. This scenario applies to planets migrating within gravito-turbulent regions of protoplanetary discs as well as stars and black holes embedded in the outskirts of active galactic nucleus (AGN) accretion discs. When the turbulence level is low, linear Lindblad torque persists in the background of stochastic forces and its accumulative effect can still dominate over relatively long time-scales. However, in the presence of very stronger turbulence, classical flow patterns around the companion embedded in the disc are disrupted, leading to significant deviations from the expectations of classical Type I migration theory over arbitrarily long time-scales. Our findings suggest that the stochastic nature of turbulent migration can prevent low-mass companions from monotonically settling into universal migration traps within the traditional laminar disc framework, thus reducing the frequency of three-body interactions and hierarchical mergers compared to previously expected. We propose a scaling for the transition mass ratio from classical to chaotic migration q ∝ αR, where αR is the Reynolds viscosity stress parameter, which can be further tested and refined by conducting extensive simulations over the relevant parameter space.

Radiation Transport Two-temperature GRMHD Simulations of Warped Accretion Disks

In many black hole (BH) systems, the accretion disk is expected to be misaligned with respect to the BH spin axis. If the scale height of the disk is much smaller than the misalignment angle, the spin of the BH can tear the disk into multiple, independently precessing “sub-disks.” This is most likely to happen during outbursts in black hole X-Ray binaries (BHXRBs) and in active galactic nuclei (AGNs) accreting above a few percent of the Eddington limit, because the disk becomes razor-thin. Disk tearing has the potential to explain variability phenomena including quasi-periodic oscillations in BHXRBs and changing-look phenomena in AGNs. Here, we present the first radiative two-temperature general relativistic magnetohydrodynamic (GRMHD) simulation of a strongly tilted (65°) accretion disk around an M BH = 10 M ⊙ BH, which tears and precesses. This leads to luminosity swings between a few percent and 50% of the Eddington limit on sub-viscous timescales. Surprisingly, even where the disk is radiation-pressure-dominated, the accretion disk is thermally stable over t ≳ 14,000 r g /c. This suggests warps play an important role in stabilizing the disk against thermal collapse. The disk forms two nozzle shocks perpendicular to the line of nodes where the scale height of the disk decreases tenfold and the electron temperature reaches T e ∼ 108–109 K. In addition, optically thin gas crossing the tear between the inner and outer disk gets heated to T e ∼ 108 K. This suggests that warped disks may emit a Comptonized spectrum that deviates substantially from idealized models.

CuHARM: A New GPU-accelerated GRMHD Code and Its Application to ADAF Disks

We introduce a new GPU-accelerated general relativistic magnetohydrodynamic code based on HARM, which we call cuHARM. The code is written in CUDA-C and uses OpenMP to parallelize multi-GPU setups. Our code allows us to run high-resolution simulations of accretion disks and the formation and structure of jets without a need for multinode supercomputer infrastructure. A 2563 simulation is well within the reach of an Nvidia DGX-V100 server, with the computation being about 10 times faster than if only the CPU is used. We use this code to examine several disk structures all in the SANE state. We find that (i) increasing the magnetic field in the SANE state does not affect the mass accretion rate; (ii) simultaneously increasing the disk size and the magnetic field, while keeping the ratio of energies fixed, leads to the destruction of the jet once the magnetic flux through the horizon decreases below a certain limit (this demonstrates that the existence of a jet is a function of the initial mass distribution, and not of the initial intensity of the magnetic field, since the magnetorotational instability dictates the evolution of the magnetic field); and (iii) the structure of the jet is a weak function of the adiabatic index of the gas, with relativistic gas tending to have a wider jet.

GRRMHD simulations of MAD accretion discs declining from super-Eddington to sub-Eddington accretion rates

ABSTRACT We present two general relativistic radiation magnetohydrodynamics (GRRMHD) simulations of magnetically arrested discs (MADs) around non-spinning (a* = 0) and spinning (a* = 0.9) supermassive black holes (BHs). In each simulation, the mass accretion rate is decreased with time such that we sample Eddington-scaled rates over the range $3 \gtrsim \dot{M}/\dot{M}_{\rm {Edd}}\gtrsim 0.3$. For the non-spinning BH model, the total and radiative efficiencies increase as the accretion rate decreases, varying over the range $\eta _{\rm {tot}}\sim 9\!-\!16{{\ \rm per\ cent}}$ and $\eta _{\rm {rad}}\sim 6{-}12{{\ \rm per\ cent}}$, respectively. This model shows very little jet activity. In contrast, the spinning BH model has a strong relativistic jet powered by spin energy extracted from the BH. The jet power declines with accretion rate such that $\eta _{\rm {jet}}\sim 18{-}39{{\ \rm per\ cent}}$ while the total and radiative efficiencies are $\eta _{\rm {tot}}\sim 64{-}100{{\ \rm per\ cent}}$ and $\eta _{\rm {rad}}\sim 45{-}79{{\ \rm per\ cent}}$, respectively. We confirm that mildly sub-Eddington discs can extract substantial power from a spinning BH, provided they are in the MAD state. The jet profile out to $100\, GM/c^2$ is roughly parabolic with a power-law index of k ≈ 0.43−0.53 during the sub-Eddington evolution. Both models show significant variability in the outgoing radiation which is likely associated with episodes of magnetic flux eruptions. The a* = 0.9 model shows semiregular variations with a period of $\sim 2000\, GM/c^3$ over the final $\sim 10\, 000\, GM/c^3$ of the simulation, which suggests that magnetic flux eruptions may be an important source of quasi-periodic variability. For the simulated accretion rates, the a* = 0 model is spinning up while the a* = 0.9 model is spinning down. Spinup–spindown equilibrium of the BH will likely be achieved at 0.5 < a*, eq < 0.6, assuming continuous accretion in the MAD state.

Fast flavor instability in hypermassive neutron star disk outflows

We examine the effect of neutrino flavor transformation by the fast flavor instability (FFI) on long-term mass ejection from accretion disks formed after neutron star mergers. Neutrino emission and absorption in the disk set the composition of the disk ejecta, which subsequently undergoes $r$-process nucleosynthesis upon expansion and cooling. Here we perform 28 time-dependent, axisymmetric, viscous-hydrodynamic simulations of accretion disks around hypermassive neutron stars (HMNSs) of variable lifetime, using a 3-species neutrino leakage scheme for emission and an annular-lightbulb scheme for absorption. We include neutrino flavor transformation due the FFI in a parametric way, by modifying the absorbed neutrino fluxes and temperatures, allowing for flavor mixing at various levels of flavor equilibration, and also in a way that aims to respect the lepton-number preserving symmetry of the neutrino self-interaction Hamiltonian. We find that for a promptly-formed black hole (BH), the FFI lowers the average electron fraction of the disk outflow due to a decrease in neutrino absorption, driven primarily by a drop in electron neutrino/antineutrino flux upon flavor mixing. For a long-lived HMNS, the disk emits more heavy lepton neutrinos and reabsorbs more electron neutrinos than for a BH, with a smaller drop in flux compensated by a higher neutrino temperature upon flavor mixing. The resulting outflow has a broader electron fraction distribution, a more proton-rich peak, and undergoes stronger radiative driving. Disks with intermediate HMNS lifetimes show results that fall in between these two limits. In most cases, the impact of the FFI on the outflow is moderate, with changes in mass ejection, average velocity, and average electron fraction of order $\sim 10\%$, and changes in the lanthanide/actinide mass fraction of up to a factor $\sim 2$.

Dynamical Unification of Tidal Disruption Events

The ∼100 tidal disruption events (TDEs) observed so far exhibit a wide range of emission properties both at peak and over their lifetimes. Some TDEs radiate predominantly at X-ray energies, while others radiate chiefly at UV and optical wavelengths. While the peak luminosities across TDEs show distinct properties, the evolutionary behavior can also vary between TDEs with similar peak emission properties. In particular, for optical TDEs, while their UV and optical emissions decline somewhat following the fallback pattern, some events can greatly rebrighten in X-rays at late time. In this Letter, we conduct three-dimensional general relativistic radiation magnetohydrodynamics simulations of TDE accretion disks at varying accretion rates in the regime of super-Eddington accretion. We make use of Monte Carlo radiative transfer simulations to calculate the reprocessed spectra at various inclinations and at different evolutionary stages. We confirm the unified model proposed by Dai et al., which predicts that the observed emission largely depends on the viewing angle of the observer with respect to the disk orientation. Furthermore, we find that disks with higher accretion rates have elevated wind and disk densities, which increases the reprocessing of the high-energy radiation and thus generally augments the optical-to-X-ray flux ratio along a particular viewing angle. This implies that at later times, as the accretion level declines, we expect that more X-rays will leak out along intermediate viewing angles. Such dynamical model for TDEs can provide a natural explanation for the diversity in the emission properties observed in TDEs at peak and along their temporal evolution.

Extending GRMHD for thin disks to non-Kerr spacetimes

Black hole based tests of general relativity have proliferated in recent times with new and improved detectors and telescopes. Modeling of the black hole neighborhood, where most of the radiation carrying strong-field signature originates, is of utmost importance for robust and accurate constraints on possible violations of general relativity. As a first step, this paper presents the extension of general relativistic magnetohydrodynamic simulations of thin accretion disks to parametrically deformed black holes that generalize the Kerr solution. The extension is based on harmpi, a publicly available member of the harm family of codes, and uses a phenomenological metric to study parametric deviations away from Kerr. The extended model is used to study the disk structure, stability, and radiative efficiency. We also compute the Fe $\mathrm{K}\ensuremath{\alpha}$ profiles in simplified scenarios and present an outlook for the future.

Observational properties of puffy discs: radiative GRMHD spectra of mildly sub-Eddington accretion

ABSTRACT Numerical general relativistic radiative magnetohydrodynamic simulations of accretion discs around a stellar-mass black hole with a luminosity above 0.5 of the Eddington value reveal their stratified, elevated vertical structure. We refer to these thermally stable numerical solutions as puffy discs. Above a dense and geometrically thin core of dimensionless thickness h/r ∼ 0.1, crudely resembling a classic thin accretion disc, a puffed-up, geometrically thick layer of lower density is formed. This puffy layer corresponds to h/r ∼ 1.0, with a very limited dependence of the dimensionless thickness on the mass accretion rate. We discuss the observational properties of puffy discs, particularly the geometrical obscuration of the inner disc by the elevated puffy region at higher observing inclinations, and collimation of the radiation along the accretion disc spin axis, which may explain the apparent super-Eddington luminosity of some X-ray objects. We also present synthetic spectra of puffy discs, and show that they are qualitatively similar to those of a Comptonized thin disc. We demonstrate that the existing xspec spectral fitting models provide good fits to synthetic observations of puffy discs, but cannot correctly recover the input black hole spin. The puffy region remains optically thick to scattering; in its spectral properties, the puffy disc roughly resembles that of a warm corona sandwiching the disc core. We suggest that puffy discs may correspond to X-ray binary systems of luminosities above 0.3 of the Eddington luminosity in the intermediate spectral states.

Fast neutrino conversion in hydrodynamic simulations of neutrino-cooled accretion disks

The outflows from neutrino-cooled black-hole (BH) accretion disks formed in neutron-star mergers or cores of collapsing stars are expected to be neutron-rich enough to explain a large fraction of elements created by the rapid neutron-capture (r-) process, but their precise chemical composition remains elusive. Here, we investigate the role of fast neutrino flavor conversion, motivated by the findings of our post-processing analysis that shows evidence of electron-neutrino lepton-number (ELN) crossings deep inside the disk, hence suggesting possibly non-trivial effects due to neutrino flavor mixing. We implement a parametric, dynamically self-consistent treatment of fast conversion in time-dependent simulations and examine the impact on the disk and its outflows. By activating the, otherwise inefficient, emission of heavy-lepton neutrinos, fast conversions enhance the disk cooling rates and reduce the absorption rates of electron-type neutrinos, causing a reduction of the electron fraction in the disk by 0.03-0.06 and in the ejected material by 0.01-0.03. The r-process yields are enhanced by typically no more than a factor of two, rendering the overall impact of fast conversions modest. The kilonova is prolonged as a net result of increased lanthanide opacities and enhanced radioactive heating rates. We observe only mild sensitivity to the disk mass, the condition for the onset of flavor conversion, and to the considered cases of flavor mixing. Remarkably, parametric models of flavor mixing that conserve the lepton numbers per family result in an overall smaller impact than models invoking three-flavor equipartition, often assumed in previous works.

Long-term 3D MHD simulations of black hole accretion discs formed in neutron star mergers

ABSTRACT We examine the long-term evolution of accretion tori around black hole (BH) remnants of compact object mergers involving at least one neutron star, to better understand their contribution to kilonovae and the synthesis of r-process elements. To this end, we modify the unsplit magnetohydrodynamic (MHD) solver in flash 4.5 to work in non-uniform three-dimensional spherical coordinates, enabling more efficient coverage of a large dynamic range in length scales while exploiting symmetries in the system. This modified code is used to perform BH accretion disc simulations that vary the initial magnetic field geometry and disc compactness, utilizing a physical equation of state, a neutrino leakage scheme for emission and absorption, and modelling the BH’s gravity with a pseudo-Newtonian potential. Simulations run for long enough to achieve a radiatively inefficient state in the disc. We find robust mass ejection with both poloidal and toroidal initial field geometries, and suppressed outflow at high disc compactness. With the included physics, we obtain bimodal velocity distributions that trace back to mass ejection by magnetic stresses at early times, and to thermal processes in the radiatively inefficient state at late times. The electron fraction distribution of the disc outflow is broad in all models, and the ejecta geometry follows a characteristic hourglass shape. We test the effect of removing neutrino absorption or nuclear recombination with axisymmetric models, finding ${\sim} 50{{\ \rm per\ cent}}$ less mass ejection and more neutron-rich composition without neutrino absorption, and a subdominant contribution from nuclear recombination. Tests of the MHD and neutrino leakage implementations are included.

General Relativistic Radiation Magnetohydrodynamics Simulations of Black Hole Accretion Disks: Comparison of Methods Based on Variable Eddington Tensor and Based on M1 Closure

We perform 2.5-dimensional general relativistic radiation magnetohydrodynamics simulations of black hole accretion disks and disk winds in the range of mass accretion rate from 0.1 to 104.5 times the Eddington limit. In this paper, we compare the results of the INAZUMA code, in which the frequency-integrated time-dependent radiation transfer equation is solved in order to evaluate the Eddington tensor, with those of the first momentum (M1) approximation method. In both methods, accretion disks and disk winds appear, and there is no remarkable difference in accretion rate, outflow rate, or luminosity. However, the significant difference in the radiation field appears around the rotation axis. In the M1 method, the radial component of the radiation flux tends to be amplified owing to unphysical radiation collisions. Such an enhancement of the outward radiation flux does not appear in INAZUMA. Also, the problem of radiation not reaching the rotation axis occurs with M1, but not with INAZUMA. Our results indicate that the radiation transfer equation should be solved to obtain the accurate radiation field in the optically thin region around the rotation axis.