No Arabic abstract
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 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).
Observational evidence accumulated over the past decade indicates that accretion discs in X-ray binaries are viscously stable unless they accrete very close to the Eddington limit. This is at odds with the most basic standard accretion disc theory, but could be explained by either having the discs to be much cooler whereby they are not radiation pressure dominated, or by a more sophisticated viscosity law. Here we argue that the latter is taking place in practice, on the basis of a stability analysis that assumes that the magneto-rotational-instability (MRI) responsible for generating the turbulent stresses inside the discs is also the source for a magnetically dominated corona. We show that observations of stable discs in the high/soft states of black hole binaries, on the one hand, and of the strongly variable microquasar GRS 1915+105 on the other, can all be explained if the magnetic turbulent stresses inside the disc scale proportionally to the geometric mean of gas and total pressure with a constant of proportionality (viscosity parameter) having a value of a few times 10^{-2}. Implications for bright AGN are also briefly discussed
Origin of hydrodynamical instability and turbulence in the Keplerian accretion disc as well as similar laboratory shear flows, e.g. plane Couette flow, is a long standing puzzle. These flows are linearly stable. Here we explore the evolution of perturbation in such flows in the presence of an additional force. Such a force, which is expected to be stochastic in nature hence behaving as noise, could be result of thermal fluctuations (however small be), Brownian ratchet, grain-fluid interactions and feedback from outflows in astrophysical discs etc. We essentially establish the evolution of nonlinear perturbation in the presence of Coriolis and external forces, which is modified Landau equation. We show that even in the linear regime, under suitable forcing and Reynolds number, the otherwise least stable perturbation evolves to a very large saturated amplitude, leading to nonlinearity and plausible turbulence. Hence, forcing essentially leads a linear stable mode to unstable. We further show that nonlinear perturbation diverges at a shorter timescale in the presence of force, leading to a fast transition to turbulence. Interestingly, emergence of nonlinearity depends only on the force but not on the initial amplitude of perturbation, unlike original Landau equation based solution.
The rather elusive High-Frequency Quasi-Periodic Oscillations(HFQPO) observed in the X-ray lightcurve of black holes have been seen in a wide range of frequencies, even within one source. It is also notable to have been detected in pairs of HFQPOs with a close to integer ratio between the frequencies. The aim of this paper is to investigate some of the possible observable that we could obtain from having the Rossby Wave Instability (RWI) active in the accretion disc surrounding the compact object. Using the newly developed GR-AMRVAC code able to follow the evolution of the RWI in a full general relativistic framework, we explore how RWI can reproduce observed HFQPO frequencies ratios and if it is compatible with the observations. In order to model the emission coming from the disc we have linked our general relativistic simulations to the general relativistic ray-tracing GYOTO code and delivered synthetic observables that can be confronted to actual data from binary systems hosting HFQPOs.}{We have demonstrated in our study that some changes in the physical conditions prevailing in the part of the disc where RWI can be triggered leads to various dominant RWI modes whose ratio recovers frequency ratios observed in various X-ray binary systems. In addition, we have also highlighted that when RWI is triggered near the last stable orbit of a spinning black hole, the amplitude of the X-ray modulation increases with the spin of the black hole. Revisiting published data on X-ray binary systems, we show that this type of relationship actually exists in the five systems where an indirect measure of the spin of the black hole is available.
The magneto-rotational instability (MRI) is the most likely mechanism for transportation of angular momentum and dissipation of energy within hot, ionized accretion discs. This instability is produced through the interactions of a differentially rotating plasma with an embedded magnetic field. Like all substances in nature, the plasma in an accretion disc has the potential to become magnetically polarized when it interacts with the magnetic field. In this paper, we study the effect of this magnetic susceptibility, parameterized by $chi_m$, on the MRI, specifically within the context of black hole accretion. We find from a linear analysis within the Newtonian limit that the minimum wavelength of the first unstable mode and the wavelength of the fastest growing mode are shorter in paramagnetic ($chi_m>0$) than in diamagnetic ($chi_m<0$) discs, all other parameters being equal. Furthermore, the magnetization parameter (ratio of gas to magnetic pressure) in the saturated state should be smaller when the magnetic susceptibility is positive than when it is negative. We confirm this latter prediction through a set of numerical simulations of magnetically polarized black hole accretion discs. We additionally find that the vertically integrated stress and mass accretion rate are somewhat larger when the disc is paramagnetic than when it is diamagnetic. If astrophysical discs are able to become magnetically polarized to any significant degree, then our results would be relevant to properly interpreting observations.