No Arabic abstract
We use three dimensional radiation magneto-hydrodynamic simulations to study the effects of magnetic fields on the energy transport and structure of radiation pressure dominated main sequence massive star envelopes at the region of the iron opacity peak. We focus on the regime where the local thermal timescale is shorter than the dynamical timescale, corresponding to inefficient convective energy transport. We begin with initially weak magnetic fields relative to the thermal pressure, from 100-1000G in differing geometries. The unstable density inversion amplifies the magnetic field, increasing the magnetic energy density to values close to equipartition with the turbulent kinetic energy density. By providing pressure support, the magnetic fields presence significantly increases the density fluctuations in the turbulent envelope, thereby enhancing the radiative energy transport by allowing photons to diffuse out through low density regions. Magnetic buoyancy brings small scale magnetic fields to the photosphere and increases the vertical energy transport with the energy advection velocity proportional to the Alfven velocity, although in all cases we study photon diffusion still dominates the energy transport. The increased radiative and advective energy transport causes the stellar envelope to shrink by several scale heights. We also find larger turbulent velocity fluctuations compared to the purely hydrodynamic case, reaching $approx$ 100 km/s at the stellar photosphere. The photosphere also shows vertical oscillations with similar averaged velocities and periods of a few hours. The increased turbulent velocity and oscillations will have strong impacts on the line broadening and periodic signals in massive stars.
Stellar feedback in the form of radiation pressure and magnetically-driven collimated outflows may limit the maximum mass that a star can achieve and affect the star-formation efficiency of massive pre-stellar cores. Here we present a series of 3D adaptive mesh refinement radiation-magnetohydrodynamic simulations of the collapse of initially turbulent, massive pre-stellar cores. Our simulations include radiative feedback from both the direct stellar and dust-reprocessed radiation fields, and collimated outflow feedback from the accreting stars. We find that protostellar outflows punches holes in the dusty circumstellar gas along the stars polar directions, thereby increasing the size of optically thin regions through which radiation can escape. Precession of the outflows as the stars spin axis changes due to the turbulent accretion flow further broadens the outflow, and causes more material to be entrained. Additionally, the presence of magnetic fields in the entrained material leads to broader entrained outflows that escape the core. We compare the injected and entrained outflow properties and find that the entrained outflow mass is a factor of $sim$3 larger than the injected mass and the momentum and energy contained in the entrained material are $sim$25% and $sim$5% of the injected momentum and energy, respectively. As a result, we find that, when one includes both outflows and radiation pressure, the former are a much more effective and important feedback mechanism, even for massive stars with significant radiative outputs.
(Abridged) Stars more massive than $20-30M_{odot}$ are so luminous that the radiation force on the cooler, more opaque outer layers can balance or exceed the force of gravity. These near or super-Eddington outer envelopes represent a long standing challenge for calculating the evolution of massive stars in one dimension, a situation that limits our understanding of the stellar progenitors of some of the most exciting and energetic explosions in the universe. In particular, the proximity to the Eddington limit has been the suspected cause for the variability, large mass loss rate and giant eruptions of an enigmatic class of massive stars: the luminous blue variables (LBVs). When in quiescence, LBVs are usually found on the hot ($T_{eff} approx 2 - 4 times 10^4$ K) S Dor instability strip. While in outburst, most LBVs stay on the cold S Dor instability strip with a $T_{eff} approx 9000$ K. Here we show that physically realistic three dimensional global radiation hydrodynamic simulations of radiation dominated massive stars with the largest supercomputers in the world naturally reproduce many observed properties of LBVs, specifically their location in the Hertzsprung-Russell (HR) diagram and their episodic mass loss with rates of $10^{-7}-10^{-5} M_{odot}/yr$. The helium opacity peak is found to play an important role to determine these properties, which is not realized in the traditional one dimensional models of massive stars. The simulations also predict that convection causes irregular envelope oscillations and 10-30% brightness variations on a typical timescale of a few days. The variability is more prominent in our models that are on the cool part of the S Dor instability. These calculations pave the way to a quantitative understanding of the structure, stability and the dominant mode of mass loss of massive stars.
$tau$ Sco, a well-studied magnetic B-type star in the Upper Sco association, has a number of surprising characteristics. It rotates very slowly and shows nitrogen excess. Its surface magnetic field is much more complex than a purely dipolar configuration which is unusual for a magnetic massive star. We employ the CMFGEN radiative transfer code to determine the fundamental parameters and surface CNO and helium abundances. Then, we employ MESA and GENEC stellar evolution models accounting for the effects of surface magnetic fields. To reconcile $tau$ Scos properties with single-star models, an increase is necessary in the efficiency of rotational mixing by a factor of 3 to 10 and in the efficiency of magnetic braking by a factor of 10. The spin down could be explained by assuming a magnetic field decay scenario. However, the simultaneous chemical enrichment challenges the single-star scenario. Previous works indeed suggested a stellar merger origin for $tau$ Sco. However, the merger scenario also faces similar challenges as our magnetic single-star models to explain $tau$ Scos simultaneous slow rotation and nitrogen excess. In conclusion, the single-star channel seems less likely and versatile to explain these discrepancies, while the merger scenario and other potential binary-evolution channels still require further assessment as to whether they may self-consistently explain the observables of $tau$ Sco.
Surface magnetic fields have a strong impact on stellar mass loss and rotation and, as a consequence, on the evolution of massive stars. In this work we study the influence of an evolving dipolar surface fossil magnetic field with an initial field strength of 4 kG on the characteristics of 15 M$_{odot}$ solar metallicity models using the Geneva stellar evolution code. Non-rotating and rotating models considering two different scenarios for internal angular momentum transport are computed, including magnetic field evolution, mass-loss quenching, and magnetic braking. Magnetic field evolution results in weakening the initially strong magnetic field, however, in our models an observable magnetic field is still maintained as the star evolves towards the red supergiant phase. At the given initial mass of the models, mass-loss quenching is modest. Magnetic braking greatly enhances chemical element mixing if radial differential rotation is allowed for, on the other hand, the inclusion of surface magnetic fields yields a lower surface enrichment in the case of near solid-body rotation. Models including surface magnetic fields show notably different trends on the Hunter diagram (plotting nitrogen abundance vs $v sin i$) compared to those that do not. The magnetic models agree qualitatively with the anomalous `Group 2 stars, showing slow surface rotation and high surface nitrogen enhancement on the main sequence.
One of the defining processes which govern massive star evolution is their continuous mass loss via dense, supersonic line-driven winds. In the case of those OB stars which also host a surface magnetic field, the interaction between that field and the ionized outflow leads to complex circumstellar structures known as magnetospheres. In this contribution, we review recent developments in the field of massive star magnetospheres, including current efforts to characterize the largest magnetosphere surrounding an O star: that of NGC 1624-2. We also discuss the potential of the `analytic dynamical magnetosphere (ADM) model to interpret multi-wavelength observations. Finally, we examine the possible effects of -- heretofore undetected -- small-scale magnetic fields on massive star winds and compare their hypothetical consequences to existing, unexplained observations.