Stratified Global MHD Models of Accretion Disks in Semi-Detached Binaries

Abstract

We present results of the first global magnetohydrodynamic (MHD) simulations of accretion disks fed by Roche lobe overflow, including vertical stratification, in order to investigate the roles of spiral shocks, magnetorotational instability (MRI), and the accretion stream on disk structure and evolution. Our models include a simple treatment of gas thermodynamics, with orbital Mach numbers at the inner edge of the disk $M_{\mathrm{in}}$ of 5 and 10. We find mass accretion rates to vary considerably on all time scales, with only the Mach 5 model reaching a clear quasi-stationary state. For Mach 10, the model undergoes an outside-in magnetically-driven accretion event occurring on a time scale of $\sim 10$ orbital periods of the binary. Both models exhibit spiral shocks inclined with respect to the binary plane, with their position and inclination changing rapidly. However, the time-averaged location of these shocks in the equatorial plane is well-fit by simple linear models. MRI turbulence in the disk generates toroidal magnetic field patterns (butterfly diagrams) that are in some cases irregular, perhaps due to interaction with spiral structure. While many of our results are in good agreement with local studies, we find some features (most notably those related to spiral shocks) can only be captured in global models such as studied here. Thus, while global studies remain computationally expensive -- even as idealized models -- they are essential (along with more sophisticated treatment of radiation transport and disk thermodynamics) for furthering our understanding of accretion in binary systems.

30 pages, 13 figures; Accepted for publication in ApJ

References in corpus (18)

• The NumPy array: a structure for efficient numerical computation
• Mayavi: a package for 3D visualization of scientific data
• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• An Unsplit Godunov Method for Ideal MHD via Constrained Transport in Three Dimensions
• MHD simulations of the magnetorotational instability in a shearing box with zero net flux. II. The effect of transport coefficients
• 3D Radiation Non-ideal Magnetohydrodynamical Simulations Of The Inner Rim In Protoplanetary Disks
• On the Magnetic Prandtl Number Behavior of Accretion Disks
• The Comptonisation of accretion disc X-ray emission: Consequences for X-ray reflection and the geometry of AGN coronae
• Convection Causes Enhanced Magnetic Turbulence in Accretion Disks in Outburst
• Dwarf Nova Outbursts with Magnetorotational Turbulence
• Magnetic Turbulence and Thermodynamics in the Inner Region of Protoplanetary Discs
• Impact of convection and resistivity on angular momentum transport in dwarf novae
• Convective Quenching of Field Reversals in Accretion Disc Dynamos
• 3D modelling of accretion disc in eclipsing binary system V1239 Her
• Spiral structures and temperature distribution in the quiescent accretion disc of the cataclysmic binary V2051 Ophiuchi
• L1 Stream Deflection and Ballistic Launching at the Disk Bow Shock: An Absorption-line Velocity Analysis in Semi-detached Binaries
• Analysis of spiral structure in the accretion disc of the novalike Cataclysmic Variable EC21178--5417
• A spectroscopic analysis of the eclipsing nova-like EC21178-5417 -- discovery of spiral density structures

Cited by in corpus (7)

• Disc tearing: implications for black hole accretion and AGN variability
• Modified models of radiation pressure instability in application to 10, 10${}^5$, and 10${}^7$ $M_{\odot }$ accreting black holes
• Investigating lack of accretion disk eccentricity growth in a global 3D MHD simulation of a superhump system
• Two-dimensional simulations of disks in close binaries: Simulating outburst cycles in cataclysmic variables
• Spiral-wave-driven accretion in quiescent dwarf novæ
• Three-dimensional simulations of accretion disks in pre-CE systems
• Gravitoturbulent dynamo in global simulations of gaseous disks