3D MHD Simulation of a Pulsationally-Driven MRI Decretion Disc

Abstract

We explore the pulsationally driven orbital mass ejection mechanism for Be star disc formation using isothermal, 3D magnetohydrodynamic (MHD) and hydrodynamic simulations. Non-radial pulsations are added to a star rotating at 95\% of critical as an inner boundary condition that feeds gas into the domain. In MHD, the initial magnetic field within the star is weak. The hydrodynamics simulation has limited angular momentum transport, resulting in repeating cycles of mass accumulation into a rotationally-supported disc at small radii followed by fall-back onto the star. The MHD simulation, conversely, has efficient (Maxwell ${\alpha }_{\mathrm{M}}$ $\sim$ 0.04) angular momentum transport provided by both of turbulent and coherent magnetic fields; a slowly decreting midplane driven by the magnetorotational instability and a supersonic wind on the surface of the disc driven by global magnetic torques. The angle and time-averaged properties near the midplane agree reasonably well with a 1D viscous decretion disc model with a modified $\tilde{\alpha }=0.5$, in which the gas transitions from a subsonic thin disc to a supersonic spherical wind at the critical point. 1D models, however, cannot capture the multi-phase decretion/angular structure seen in our simulations. Our results demonstrate that, at least under certain conditions, non-radial pulsations on the surface of a rapidly rotating, weakly magnetized star can drive a Keplerian disc with the basic properties of the viscous decretion disc paradigm, albeit coupled to a laminar wind away from the midplane. Future modeling of Be star discs should consider the possible existence of such a surface wind.

Accepted by MNRAS. 16 pages, 10 figures

References in corpus (10)

• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• Dynamical Simulations of Magnetically Channeled Line-Driven Stellar Winds: II. The Effects of Field-Aligned Rotation
• Dynamical Simulations of Magnetically Channeled Line-Driven Stellar Winds: III. Angular Momentum Loss and Rotational Spindown
• Giant Outbursts in Be/X-ray Binaries
• Non-LTE Monte Carlo Radiative Transfer: II. Non-Isothermal Solutions for Viscous Keplerian Disks
• The signature of the magnetorotational instability in the Reynolds and Maxwell stress tensors in accretion discs
• A Pulsational Mechanism for Producing Keplerian Disks around Be Stars
• Modelling decretion discs in Be/X-ray binaries
• Long-term evolution of accretion discs in Be/X-ray binaries
• Magnetorotational instability in decretion disks of critically rotating stars and the outer structure of Be and Be/X-ray disks

Cited by in corpus (1)

• A survey for variable young stars with small telescopes: VI -- Analysis of the outbursting Be stars NSW284, Gaia19eyy, and VES263