skip to main content

Title: Modelling magnetically channeled winds in 3D – I. Isothermal simulations of a magnetic O supergiant

In this paper we present the first set of 3D magnetohydrodynamic (MHD) simulations performed with the riemann geomesh code. We study the dynamics of the magnetically channeled winds of magnetic massive stars in full three dimensions using a code that is uniquely suited to spherical problems. Specifically, we perform isothermal simulations of a smooth wind on a rotating star with a tilted, initially dipolar field. We compare the mass-loss, angular momentum loss, and magnetospheric dynamics of a template star (with the properties that are reminiscent of the O4 supergiant ζ Pup) over a range of rotation rates, magnetic field strengths, and magnetic tilt angles. The simulations are run up to a quasi-steady state and the results are observed to be consistent with the existing literature, showing the episodic centrifugal breakout events of the mass outflow, confined by the magnetic field loops that form the closed magnetosphere of the star. The catalogued results provide perspective on how angular-momentum loss varies for different configurations of rotation rate, magnetic field strength, and large magnetic tilt angles. In agreement with previous 2D MHD studies, we find that high magnetic confinement reduces the overall mass-loss rate, and higher rotation increases the mass-loss rate. This and future studies will be used to estimate the angular-momentum evolution, spin-down time, and mass-loss evolution of magnetic massive stars as a function of magnetic field strength, rotation rate, and dipole tilt.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
p. 237-255
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    We consider the general problem of a Parker-type non-relativistic isothermal wind from a rotating and magnetic star. Using the magnetohydrodynamics code athena++, we construct an array of simulations in the stellar rotation rate Ω* and the isothermal sound speed cT, and calculate the mass, angular momentum, and energy loss rates across this parameter space. We also briefly consider the 3D case, with misaligned magnetic and rotation axes. We discuss applications of our results to the spin-down of normal stars, highly irradiated exoplanets, and to nascent highly magnetic and rapidly rotating neutron stars born in massive star core-collapse.

    more » « less

    Magnetic fields are observed in massive Ap/Bp stars and are presumably present in the radiative zone of solar-like stars. To date, there is no clear understanding of the dynamics of the magnetic field in stably stratified layers. A purely toroidal magnetic field configuration is known to be unstable, developing mainly non-axisymmetric modes. Rotation and a poloidal field component may lead to stabilization. Here we perform global MHD simulations with the EULAG-MHD code to explore the evolution of a toroidal magnetic field located in a layer whose Brunt-Väisälä frequency resembles the lower solar tachocline. Our numerical experiments allow us to explore the initial unstable phase as well as the long-term evolution of such field. During the first Alfven cycles, we observe the development of the Tayler instability with the prominent longitudinal wavenumber, m = 1. Rotation decreases the growth rate of the instability and eventually suppresses it. However, after a stable phase, energy surges lead to the development of higher-order modes even for fast rotation. These modes extract energy from the initial toroidal field. Nevertheless, our results show that sufficiently fast rotation leads to a lower saturation energy of the unstable modes, resulting in a magnetic topology with only a small fraction of poloidal field, which remains steady for several hundreds of Alfven traveltimes. The system then becomes turbulent and the field is prone to turbulent diffusion. The final toroidal–poloidal configuration of the magnetic field may represent an important aspect of the field generation and evolution in stably stratified layers.

    more » « less

    Solar-type stars, which shed angular momentum via magnetized stellar winds, enter the main sequence with a wide range of rotational periods Prot. This initially wide range of rotational periods contracts and has mostly vanished by a stellar age $t\sim {0.6}\, {\rm Gyr}$, after which Solar-type stars spin according to the Skumanich relation $P_\text{rot}\propto \sqrt{t}$. Magnetohydrodynamic stellar wind models can improve our understanding of this convergence of rotation periods. We present wind models of 15 young Solar-type stars aged ∼24 Myr to ∼0.13 Gyr. With our previous wind models of stars aged ∼0.26 and ∼0.6 Gyr we obtain 30 consistent three-dimensional wind models of stars mapped with Zeeman–Doppler imaging – the largest such set to date. The models provide good cover of the pre-Skumanich phase of stellar spin-down in terms of rotation, magnetic field, and age. We find the mass-loss rate $\dot{M}\propto \Phi ^{{0.9\pm 0.1}}$ with a residual spread of ∼150 per cent and the wind angular momentum loss rate $\dot{J}\propto {}P_\text{rot}^{-1} \Phi ^{1.3\pm 0.2}$ with a residual spread of ∼500 per cent where Φ is the unsigned surface magnetic flux. When comparing different magnetic field scalings for each single star we find a gradual reduction in the power-law exponent with increasing magnetic field strength.

    more » « less
  4. ABSTRACT The time evolution of angular momentum and surface rotation of massive stars are strongly influenced by fossil magnetic fields via magnetic braking. We present a new module containing a simple, comprehensive implementation of such a field at the surface of a massive star within the Modules for Experiments in Stellar Astrophysics (mesa) software instrument. We test two limiting scenarios for magnetic braking: distributing the angular momentum loss throughout the star in the first case, and restricting the angular momentum loss to a surface reservoir in the second case. We perform a systematic investigation of the rotational evolution using a grid of OB star models with surface magnetic fields (M⋆ = 5–60 M⊙, Ω/Ωcrit = 0.2–1.0, Bp = 1–20 kG). We then employ a representative grid of B-type star models (M⋆ = 5, 10, 15 M⊙, Ω/Ωcrit = 0.2, 0.5, 0.8, Bp = 1, 3, 10, 30 kG) to compare to the results of a recent self-consistent analysis of the sample of known magnetic B-type stars. We infer that magnetic massive stars arrive at the zero-age main sequence (ZAMS) with a range of rotation rates, rather than with one common value. In particular, some stars are required to have close-to-critical rotation at the ZAMS. However, magnetic braking yields surface rotation rates converging to a common low value, making it difficult to infer the initial rotation rates of evolved, slowly rotating stars. 
    more » « less

    Stars form from the gravitational collapse of turbulent, magnetized molecular cloud cores. Our non-ideal MHD simulations reveal that the intrinsically anisotropic magnetic resistance to gravity during the core collapse naturally generates dense gravomagneto sheetlets within inner protostellar envelopes – disrupted versions of classical sheet-like pseudo-discs. They are embedded in a magnetically dominant background, where less dense materials flow along the local magnetic field lines and accumulate in the dense sheetlets. The sheetlets, which feed the disc predominantly through its upper and lower surfaces, are the primary channels for mass and angular momentum transfer from the envelope to the disc. The protostellar disc inherits a small fraction (up to 10 per cent) of the magnetic flux from the envelope, resulting in a disc-averaged net vertical field strength of 1–10 mG and a somewhat stronger toroidal field, potentially detectable through ALMA Zeeman observations. The inherited magnetic field from the envelope plays a dominant role in disc angular momentum evolution, enabling the formation of gravitationally stable discs in cases where the disc field is relatively well-coupled to the gas. Its influence remains significant even in marginally gravitationally unstable discs formed in the more magnetically diffusive cases, removing angular momentum at a rate comparable to or greater than that caused by spiral arms. The magnetically driven disc evolution is consistent with the apparent scarcity of prominent spirals capable of driving rapid accretion in deeply embedded protostellar discs. The dense gravomagneto sheetlets observed in our simulations may correspond to the ‘accretion streamers’ increasingly detected around protostars.

    more » « less