Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStrongly magnetized accretion disks: structure and accretion from global magnetohydrodynamic simulations
112
0
0.0
(
0
)
Added by
Bhupendra Mishra Dr.
Publication date
2019
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We use global magnetohydrodynamic simulations to study the influence of net vertical magnetic fields on the structure of geometrically thin ($H/r approx 0.05$) accretion disks in the Newtonian limit. We consider initial mid-plane gas to magnetic pressure ratios $beta_0 = 1000,, 300$ and $100$, spanning the transition between weakly and strongly magnetized accretion regimes. We find that magnetic pressure is important for the disks vertical structure in all three cases, with accretion occurring at $z/Rapprox 0.2$ in the two most strongly magnetized models. The disk midplane shows outflow rather than accretion. Accretion through the surface layers is driven mainly by stress due to coherent large scale magnetic field rather than by turbulent stress. Equivalent viscosity parameters measured from our simulations show similar dependencies on initial $beta_0$ to those seen in shearing box simulations, though the disk midplane is not magnetic pressure dominated even for the strongest magnetic field case. Winds are present but are not the dominant driver of disk evolution. Over the (limited) duration of our simulations, we find evidence that the net flux attains a quasi-steady state at levels that can stably maintain a strongly magnetized disk. We suggest that geometrically thin accretion disks in observed systems may commonly exist in a magnetically elevated state, characterized by non-zero but modest vertical magnetic fluxes, with potentially important implications for disk phenomenology in X-ray binaries (XRBs) and active galactic nuclei (AGN).
rate research
Read More
We examine the properties of strongly magnetized accretion discs in a global framework, with particular focus on the evolution of magnetohydrodynamic instabilities such as the magnetorotational instability (MRI). Work by Pessah and Psaltis showed that MRI is stabilized beyond a critical toroidal field in compressible, differentially rotating flows and, also, reported the appearance of two new instabilities beyond this field. Their results stemmed from considering geometric curvature effects due to the suprathermal background toroidal field, which had been previously ignored in weak-field studies. However, their calculations were performed under the local approximation, which poses the danger of introducing spurious behavior due to the introduction of global geometric terms in an otherwise local framework. In order to avoid this, we perform a global eigenvalue analysis of the linearized MHD equations in cylindrical geometry. We confirm that MRI indeed tends to be highly suppressed when the background toroidal field attains the Pessah-Psaltis limit. We also observe the appearance of two new instabilities that emerge in the presence of highly suprathermal toroidal fields. These results were additionally verified using numerical simulations in PLUTO. There are, however, certain differences between the the local and global results, especially in the vertical wavenumber occupancies of the various instabilities, which we discuss in detail. We also study the global eigenfunctions of the most unstable modes in the suprathermal regime, which are inaccessible in the local analysis. Overall, our findings emphasize the necessity of a global treatment for accurately modeling strongly magnetized accretion discs.
We carried out two dimensional high-resolution magnetohydrodynamic (MHD) simulations of an accretion disk around a weakly magnetized neutron star. General relativistic effects are taken into account by using pseudo-Newtonian potential of Pacz{y}nski, B., P. J. Witta (1980). When magnetic loops connect the neutron star and the accretion disk, the twist injection from the disk or from the rotating neutron star triggers expansion of the loops. Since the expanding magnetic loops prevent inflow toward the magnetic poles of the neutron star, disk matter accumulates on the boundary between the magnetosphere and the disk. Magnetic reconnection taking place in the loops creates a channel along which the disk matter can accrete and unloads the magnetosphere. This process produces quasi-periodic variation of the accretion flow in the innermost region of the disk. We found two kinds of oscillations. One is the magnetospheric oscillation regulated by magnetic reconnection. The other is the radial disk oscillation. The typical frequency of the oscillations is 100 Hz to 2 kHz. Furthermore, we predict that QPO sources inevitably accompany X-ray flares by magnetic reconnection and bipolar outflows of hot X-ray emitting plasma similar to the optical jets in protostars.
We model the vertical structure of magnetized accretion disks subject to viscous and resistive heating, and irradiation by the central star. We apply our formalism to the radial structure of magnetized accretion disks threaded by a poloidal magnetic field dragged during the process of star formation developed by Shu and coworkers. We consider disks around low mass protostars, T Tauri, and FU Orionis stars. We consider two levels of disk magnetization, $lambda_{sys} = 4$ (strongly magnetized disks), and $lambda_{sys} = 12$ (weakly magnetized disks). The rotation rates of strongly magnetized disks have large deviations from Keplerian rotation. In these models, resistive heating dominates the thermal structure for the FU Ori disk. The T Tauri disk is very thin and cold because it is strongly compressed by magnetic pressure; it may be too thin compared with observations. Instead, in the weakly magnetized disks, rotation velocities are close to Keplerian, and resistive heating is always less than 7% of the viscous heating. In these models, the T Tauri disk has a larger aspect ratio, consistent with that inferred from observations. All the disks have spatially extended hot atmospheres where the irradiation flux is absorbed, although most of the mass ($sim 90-95$ %) is in the disk midplane. With the advent of ALMA one expects direct measurements of magnetic fields and their morphology at disk scales. It will then be possible to determine the mass-to-flux ratio of magnetized accretion disks around young stars, an essential parameter for their structure and evolution. Our models contribute to the understanding of the vertical structure and emission of these disks.
We calculate the emission of protoplanetary disks threaded by a poloidal magnetic field and irradiated by the central star. The radial structure of these disks was studied by Shu and collaborators and the vertical structure was studied by Lizano and collaborators. We consider disks around low mass protostars, T Tauri stars, and FU Ori stars with different mass-to-flux ratios $lambda_{rm sys}$. We calculate the spectral energy distribution and the antenna temperature profiles at 1 mm and 7 mm convolved with the ALMA and VLA beams. We find that disks with weaker magnetization (high values of $lambda_{rm sys}$) emit more than disks with stronger magnetization (low values of $lambda_{rm sys}$). This happens because the former are denser, hotter and have larger aspect ratios, receiving more irradiation from the central star. The level of magnetization also affects the optical depth at millimeter wavelengths, being larger for disks with high $lambda_{rm sys}$. In general, disks around low mass protostars and T Tauri stars are optically thin at 7 mm while disks around FU Ori are optically thick. A qualitative comparison of the emission of these magnetized disks, including heating by an external envelope, with the observed millimeter antenna temperature profiles of HL Tau indicates that large cm grains are required to increase the optical depth and reproduce the observed 7 mm emission at large radii.
Global three dimensional magnetohydrodynamic (MHD) simulations of turbulent accretion disks are presented which start from fully equilibrium initial conditions in which the magnetic forces are accounted for and the induction equation is satisfied. The local linear theory of the magnetorotational instability (MRI) is used as a predictor of the growth of magnetic field perturbations in the global simulations. The linear growth estimates and global simulations diverge when non-linear motions - perhaps triggered by the onset of turbulence - upset the velocity perturbations used to excite the MRI. The saturated state is found to be independent of the initially excited MRI mode, showing that once the disk has expelled the initially net flux field and settled into quasi-periodic oscillations in the toroidal magnetic flux, the dynamo cycle regulates the global saturation stress level. Furthermore, time-averaged measures of converged turbulence, such as the ratio of magnetic energies, are found to be in agreement with previous works. In particular, the globally averaged stress normalized to the gas pressure, <alpha_{rm P}> = 0.034, with notably higher values achieved for simulations with higher azimuthal resolution. Supplementary tests are performed using different numerical algorithms and resolutions. Convergence with resolution during the initial linear MRI growth phase is found for 23-35 cells per scaleheight (in the vertical direction).
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
