We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA star. Improvements in MESA star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 M ? stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA star solves the fully coupled stellar structure and composition equations, and we show how this has improved the scaling of MESA's calculational speed on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit that packages all the required components needed to form a unified, maintained, and well-validated build environment for MESA. We also highlight a few tools developed by the community for rapid visualization of MESA star results.

/pdf/modules-for-experiments-in-stellar-astrophysics-mesa-planets-29j324a9bq.pdf

Modules for Experiments in Stellar Astrophysics (MESA): Planets, Oscillations, Rotation, and Massive Stars

Using a two-dimensional hydrodynamics code (PROMETHEUS), we explore the continued evolution of rotating helium stars, Mα 10 M☉, in which iron-core collapse does not produce a successful outgoing shock but instead forms a black hole of 2-3 M☉. The model explored in greatest detail is the 14 M☉ helium core of a 35 M☉ main-sequence star. The outcome is sensitive to the angular momentum. For j16 ≡ j/(1016 cm2 s-1) 3, material falls into the black hole almost uninhibited. No outflows are expected. For j16 20, the infalling matter is halted by centrifugal force outside 1000 km where neutrino losses are negligible. The equatorial accretion rate is very low, and explosive oxygen burning may power a weak equatorial explosion. For 3 j16 20, however, a reasonable value for such stars, a compact disk forms at a radius at which the gravitational binding energy can be efficiently radiated as neutrinos or converted to beamed outflow by magnetohydrodynamical (MHD) processes. These are the best candidates for producing gamma-ray bursts (GRBs). Here we study the formation of such a disk, the associated flow patterns, and the accretion rate for disk viscosity parameter α ≈ 0.001 and 0.1. Infall along the rotational axis is initially uninhibited, and an evacuated channel opens during the first few seconds. Meanwhile the black hole is spun up by the accretion (to a ≈ 0.9), and energy is dissipated in the disk by MHD processes and radiated by neutrinos. For the α = 0.1 model, appreciable energetic outflows develop between polar angles of 30° and 45°. These outflows, powered by viscous dissipation in the disk, have an energy of up to a few times 1051 ergs and a mass ~1 M☉ and are rich in 56Ni. They constitute a supernova-like explosion by themselves. Meanwhile accretion through the disk is maintained for approximately 10-20 s but is time variable (±30%) because of hydrodynamical instabilities at the outer edge in a region where nuclei are experiencing photodisintegration. Because the efficiency of neutrino energy deposition is sensitive to the accretion rate, this instability leads to highly variable energy deposition in the polar regions. Some of this variability, which has significant power at 50 ms and overtones, may persist in the time structure of the burst. During the time followed, the average accretion rate for the standard α = 0.1 and j16 = 10 model is 0.07 M☉ s-1. The total energy deposited along the rotational axes by neutrino annihilation is (1-14) × 1051 ergs, depending upon the evolution of the Kerr parameter and uncertain neutrino efficiencies. Simulated deposition of energy in the polar regions, at a constant rate of 5 × 1050 ergs s-1 per pole, results in strong relativistic outflow jets beamed to about 1% of the sky. These jets may be additionally modulated by instabilities in the sides of the "nozzle" through which they flow. The jets blow aside the accreting material, remain highly focused, and are capable of penetrating the star in ~10 s. After the jet breaks through the surface of the star, highly relativistic flow can emerge. Because of the sensitivity of the mass ejection and jets to accretion rate, angular momentum, and disk viscosity, and the variation of observational consequences with viewing angle, a large range of outcomes is possible, ranging from bright GRBs like GRB 971214 to faint GRB-supernovae like SN 1998bw. X-ray precursors are also possible as the jet first breaks out of the star. While only a small fraction of supernovae make GRBs, we predict that collapsars will always make supernovae similar to SN 1998bw. However, hard, energetic GRBs shorter than a few seconds will be difficult to produce in this model and may require merging neutron stars and black holes for their explanation.

https://iopscience.iop.org/article/10.1086/307790/pdf

Collapsars: Gamma-ray bursts and explosions in 'failed supernovae'

How massive stars die-what sort of explosion and remnant each produces-depends chiefly on the masses of their helium cores and hydrogen envelopes at death. For single stars, stellar winds are the only means of mass loss, and these are a function of the metallicity of the star. We discuss how metallicity, and a simplified prescription for its effect on mass loss, affects the evolution and final fate of massive stars. We map, as a function of mass and metallicity, where black holes and neutron stars are likely to form and where different types of supernovae are produced. Integrating over an initial mass function, we derive the relative populations as a function of metallicity. Provided that single stars rotate rapidly enough at death, we speculate on stellar populations that might produce gamma-ray bursts and jet-driven supernovae.

/pdf/how-massive-single-stars-end-their-life-3v2i056c9w.pdf

How Massive Single Stars End Their Life

amount of energy, a tiny fraction of which is sufficient to explode the star as a supernova. The authors examine our current understanding of the lives and deaths of massive stars, with special attention to the relevant nuclear and stellar physics. Emphasis is placed upon their post-helium-burning evolution. Current views regarding the supernova explosion mechanism are reviewed, and the hydrodynamics of supernova shock propagation and ‘‘fallback’’ is discussed. The calculated neutron star masses, supernova light curves, and spectra from these model stars are shown to be consistent with observations. During all phases, particular attention is paid to the nucleosynthesis of heavy elements. Such stars are capable of producing, with few exceptions, the isotopes between mass 16 and 88 as well as a large fraction of still heavier elements made by the r and p processes.

/pdf/the-evolution-and-explosion-of-massive-stars-57jc6k6s1i.pdf

The evolution and explosion of massive stars

Aims. Many topical astrophysical research areas, such as the properties of planet host stars, the nature of the progenitors of different types of supernovae and gamma ray bursts, and the evolution of galaxies, require complete and homogeneous sets of stellar models at different metallicities in order to be studied during the whole of cosmic history. We present here a first set of models for solar metallicity, where the effects of rotation are accounted for in a homogeneous way.Methods. We computed a grid of 48 different stellar evolutionary tracks, both rotating and non-rotating, at Z  = 0.014, spanning a wide mass range from 0.8 to 120 M ⊙ . For each of the stellar masses considered, electronic tables provide data for 400 stages along the evolutionary track and at each stage, a set of 43 physical data are given. These grids thus provide an extensive and detailed data basis for comparisons with the observations. The rotating models start on the zero-age main sequence (ZAMS) with a rotation rate υ ini /υ crit  = 0.4. The evolution is computed until the end of the central carbon-burning phase, the early asymptotic giant branch (AGB) phase, or the core helium-flash for, respectively, the massive, intermediate, and both low and very low mass stars. The initial abundances are those deduced by Asplund and collaborators, which best fit the observed abundances of massive stars in the solar neighbourhood. We update both the opacities and nuclear reaction rates, and introduce new prescriptions for the mass-loss rates as stars approach the Eddington and/or the critical velocity. We account for both atomic diffusion and magnetic braking in our low-mass star models.Results. The present rotating models provide a good description of the average evolution of non-interacting stars. In particular, they reproduce the observed main-sequence width, the positions of the red giant and supergiant stars in the Hertzsprung-Russell (HR) diagram, the observed surface compositions and rotational velocities. Very interestingly, the enhancement of the mass loss during the red-supergiant stage, when the luminosity becomes supra-Eddington in some outer layers, help models above 15−20 M ⊙ to lose a significant part of their hydrogen envelope and evolve back into the blue part of the HR diagram. This result has interesting consequences for the blue to red supergiant ratio, the minimum mass for stars to become Wolf-Rayet stars, and the maximum initial mass of stars that explode as type II−P supernovae.

/pdf/grids-of-stellar-models-with-rotation-i-models-from-0-8-to-3bjuk0gnjq.pdf

Grids of stellar models with rotation - I. Models from 0.8 to 120 M⊙ at solar metallicity (Z = 0.014)

Magnetic fields can be created in stably stratified (non-convective) layers in a differentially rotating star. A magnetic instability in the toroidal field (wound up by differential rotation) replaces the role of convection in closing the field amplification loop. Tayler instability is likely to be the most relevant magnetic instability. A dynamo model is developed from these ingredients, and applied to the problem of angular momentum transport in stellar interiors. It produces a predominantly horizontal field. This dynamo process is found to be more effective in transporting angular momentum than the known hydrodynamic mechanisms. It might account for the observed pattern of rotation in the solar core.

/pdf/dynamo-action-by-differential-rotation-in-a-stably-a52fb1vx2e.pdf

Dynamo action by differential rotation in a stably stratified stellar interior

Some main-sequence stars of spectral type A are observed to have a strong (0.03–3 tesla), static, large-scale magnetic field, of a chiefly dipolar shape: they are known as ‘Ap stars’1,2,3,4, such as Alioth, the fifth star in the Big Dipper. Following the discovery of these fields, it was proposed that they are remnants of the star's formation, a ‘fossil’ field5,6. An alternative suggestion is that they could be generated by a dynamo process in the star's convective core7. The dynamo hypothesis, however, has difficulty explaining high field strengths and the observed lack of a correlation with rotation. The weakness of the fossil-field theory has been the absence of field configurations stable enough to survive in a star over its lifetime. Here we report numerical simulations that show that stable magnetic field configurations, with properties agreeing with those observed, can develop through evolution from arbitrary, unstable initial fields. The results are applicable equally to Ap stars, magnetic white dwarfs and some highly magnetized neutron stars known as magnetars. This establishes fossil fields as the natural, unifying explanation for the magnetism of all these stars.

https://arxiv.org/pdf/astro-ph/0502043.pdf

A fossil origin for the magnetic field in A stars and white dwarfs

Radio pulsars are thought to born with spin periods of 0.02–0.5 s
and space velocities of 100–1,000 kms^(-1), and they are inferred to
have initial dipole magnetic fields of 10^(11)–10^(13) G. The
average space velocity of their progenitor stars is less than 15 kms^(-1),
which means that pulsars must receive a substantial ‘kick’ at birth.
Here we propose that the birth characteristics of pulsars have a
simple physical connection with each other. Magnetic fields
maintained by differential rotation between the core and envelope
of the progenitor would keep the whole star in a state of
approximately uniform rotation until 10 years before the explosion.
Such a slowly rotating core has 1,000 times less angular
momentum than required to explain the rotation of pulsars. The
specific physical process that ‘kicks’ the neutron star at birth has
not been identified, but unless its force is exerted exactly head-on
it will also cause the neutron star to rotate. We identify this
process as the origin of the spin of pulsars. Such kicks may cause a
correlation between the velocity and spin vectors of pulsars. We
predict that many neutron stars are born with periods longer than
2 s, and never become radio pulsars.

/pdf/birth-kicks-as-the-origin-of-pulsar-rotation-577faeihvk.pdf

Birth kicks as the origin of pulsar rotation

Propagation speeds are derived for the wave modes of a thin magnetic tube in an otherwise homogeneous magnetized or unmagnetized fluid. These results generalize results obtained by previous authors. There are three types of wave, a (torsional) Alfven wave and two waves which are specific for the thin tube. These are named the longitudinal and transversal tube waves, according to their polarization properties. They can be camped by radiating an MHD or acoustic wave into the surroundings of the tube. Conditions for occurrence of this acoustic damping, and the damping rates, are derived. The behavior of the waves in the solar convection zone and corona is discussed.

https://link.springer.com/content/pdf/10.1007/BF00153456.pdf

Propagation Speeds and Acoustic Damping of Waves in Magnetic Flux Tubes

The Sun supports a rich spectrum of internal waves that are continuously excited by turbulent convection. The Global Oscillation Network Group (GONG) network and the SOHO/MDI (Solar and Heliospheric Observatory/Michelson Doppler Imager) space instrument provide an exceptional database of spatially resolved observations of solar oscillations, covering more than an entire sunspot cycle (11 years). Local helioseismology is a set of tools for probing the solar interior in three dimensions using measurements of wave travel times and local mode frequencies. Local helioseismology has discovered (a) near-surface vector flows associated with convection, (b) 250 m s−1 subsurface horizontal outflows around sunspots, (c) ∼50 m s−1 extended horizontal flows around active regions (converging near the surface and diverging below), (d) the effect of the Coriolis force on convective flows and active region flows, (e) the subsurface signature of the 15 m s−1 poleward meridional flow, (f) a ±5 m s−1 time-varying depth-depen...

http://inspirehep.net/record/841923/files/arXiv:1001.0930.pdf

Hendrik C. Spruit

Papers

Dynamo action by differential rotation in a stably stratified stellar interior

A fossil origin for the magnetic field in A stars and white dwarfs

Birth kicks as the origin of pulsar rotation

Propagation Speeds and Acoustic Damping of Waves in Magnetic Flux Tubes

Local Helioseismology: Three-Dimensional Imaging of the Solar Interior