Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNew Global 3D MHD Simulations of Black Hole Disk Accretion and Outflows
125
0
0.0
(
0
)
Authors
Peter B. Dobbie
Ask ChatGPT about the research
No Arabic abstract
It is widely accepted that quasars and other active galactic nuclei (AGN) are powered by accretion of matter onto a central supermassive black hole. While numerical simulations have demonstrated the importance of magnetic fields in generating the turbulence believed necessary for accretion, so far they have not produced the high mass accretion rates required to explain the most powerful sources. We describe new global 3D simulations we are developing to assess the importance of radiation and non-ideal MHD in generating magnetized outflows that can enhance the overall rates of angular momentum transport and mass accretion.
rate research
Read More
We perform a full 3D general relativistic magnetohydrodynamical (GRMHD) simulation of an equal-mass, spinning, binary black hole approaching merger, surrounded by a circumbinary disk and with a mini-disk around each black hole. For this purpose, we evolve the ideal GRMHD equations on top of an approximated spacetime for the binary that is valid in every position of space, including the black hole horizons, during the inspiral regime. We use relaxed initial data for the circumbinary disk from a previous long-term simulation, where the accretion is dominated by a $m=1$ overdensity called the lump. We compare our new spinning simulation with a previous non-spinning run, studying how spin influences the mini-disk properties. We analyze the accretion from the inner edge of the lump to the black hole, focusing on the angular momentum budget of the fluid around the mini-disks. We find that mini-disks in the spinning case have more mass over a cycle than the non-spinning case. However, in both cases, we find most of the mass received by the black holes is delivered by the direct plunging of material from the lump. We also analyze the morphology and variability of the electromagnetic fluxes and we find they share the same periodicities of the accretion rate. In the spinning case, we find that the outflows are $8$ times stronger than the non-spinning case. Our results will be useful to understand and produce realistic synthetic light curves and spectra, which can be used in future observations.
We present the results of a 3D global magnetohydrodynamic (MHD) simulation of an AM CVn system that was aimed at exploring eccentricity growth in the accretion disc self-consistently from a first principles treatment of the MHD turbulence. No significant eccentricity growth occurs in the simulation. In order to investigate the reasons why, we ran 2D alpha disc simulations with alpha values of 0.01, 0.1, and 0.2, and found that only the latter two exhibit significant eccentricity growth. We present an equation expressing global eccentricity evolution in terms of contributing forces and use it to analyze the simulations. As expected, we find that the dominant term contributing to the growth of eccentricity is the tidal gravity of the companion star. In the 2D simulations, the alpha viscosity directly contributes to eccentricity growth. In contrast, the overall magnetic forces in the 3D simulation damp eccentricity. We also analyzed the mode-coupling mechanism of Lubow, and confirmed that the spiral wave excited by the 3:1 resonance was the dominant contributor to eccentricity growth in the 2D $alpha=0.1$ simulations, but other waves also contribute significantly. We found that the $alpha=0.1$ and 0.2 simulations had more relative mass at larger radii compared to the $alpha=0.01$ and 3D MHD simulation, which also had an effective $alpha$ of 0.01. This suggests that in 3D MHD simulations without sufficient poloidal magnetic flux, MRI turbulence does not saturate at a high enough $alpha$ to spread the disc to large enough radii to reproduce the superhumps observed in real systems.
Gas falling into a black hole (BH) from large distances is unaware of BH spin direction, and misalignment between the accretion disc and BH spin is expected to be common. However, the physics of tilted discs (e.g., angular momentum transport and jet formation) is poorly understood. Using our new GPU-accelerated code H-AMR, we performed 3D general relativistic magnetohydrodynamic simulations of tilted thick accretion discs around rapidly spinning BHs, at the highest resolution to date. We explored the limit where disc thermal pressure dominates magnetic pressure, and showed for the first time that, for different magnetic field strengths on the BH, these flows launch magnetized relativistic jets propagating along the rotation axis of the tilted disc (rather than of the BH). If strong large-scale magnetic flux reaches the BH, it bends the inner few gravitational radii of the disc and jets into partial alignment with the BH spin. On longer time scales, the simulated disc-jet system as a whole undergoes Lense-Thirring precession and approaches alignment, demonstrating for the first time that jets can be used as probes of disc precession. When the disc turbulence is well-resolved, our isolated discs spread out, causing both the alignment and precession to slow down.
The Fermi Gamma-Ray Space Telescope observations of blazars show a strong correlation between the spectral index of their gamma-ray spectra and their synchrotron peak frequency $ u_{rm{pk}}^{rm{syn}}$; additionally, the rate of Compton Dominance of these sources also seems to be a function of $ u_{rm{pk}}^{rm{syn}}$. In this work, we adopt the assumption that the nonthermal emission of blazars is primarily due to radiation by a population of Fermi-accelerated electrons in a relativistic outflow (jet) along the symmetry axis of the blazars accretion disk. Furthermore, we assume that the Compton component is related to an external photon field of photons, which are scattered from particles of the magnetohydrodynamic (MHD) wind emanating from the accretion disk. Our results reproduce well the aforementioned basic observational trends of blazar classification by varying just one parameter, namely the mass accretion rate onto the central black hole.
We perform GR-MHD simulations of outflow launching from thin accretion disks. As in the non-relativistic case, resistivity is essential for the mass loading of the disk wind. We implemented resistivity in the ideal GR-MHD code HARM3D, extending previous works (Qian et al. 2017, 2018) for larger physical grids, higher spatial resolution, and longer simulation time. We consider an initially thin, resistive disk orbiting the black hole, threaded by a large-scale magnetic flux. As the system evolves, outflows are launched from the black hole magnetosphere and the disk surface. We mainly focus on disk outflows, investigating their MHD structure and energy output in comparison with the Poynting-dominated black hole jet. The disk wind encloses two components -- a fast component dominated by the toroidal magnetic field and a slower component dominated by the poloidal field. The disk wind transitions from sub to super-Alfvenic speed, reaching velocities $simeq 0.1c$. We provide parameter studies varying spin parameter and resistivity level, and measure the respective mass and energy fluxes. A higher spin strengthens the $B_{phi}$-dominated disk wind along the inner jet. We disentangle a critical resistivity level that leads to a maximum matter and energy output for both, resulting from the interplay between re-connection and diffusion, which in combination govern the magnetic flux and the mass loading. For counter-rotating black holes the outflow structure shows a magnetic field reversal. We estimate the opacity of the inner-most accretion stream and the outflow structure around it. This stream may be critically opaque for a lensed signal, while the axial jet funnel remains optically thin.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
