We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development.

https://iopscience.iop.org/article/10.1088/0067-0049/220/1/15/pdf

Modules for experiments in stellar astrophysics (mesa): binaries, pulsations, and explosions

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

The presence of a nearby companion alters the evolution of massive stars in binary systems, leading to phenomena such as stellar mergers, x-ray binaries, and gamma-ray bursts. Unambiguous constraints on the fraction of massive stars affected by binary interaction were lacking. We simultaneously measured all relevant binary characteristics in a sample of Galactic massive O stars and quantified the frequency and nature of binary interactions. More than 70% of all massive stars will exchange mass with a companion, leading to a binary merger in one-third of the cases. These numbers greatly exceed previous estimates and imply that binary interaction dominates the evolution of massive stars, with implications for populations of massive stars and their supernovae.

/pdf/binary-interaction-dominates-the-evolution-of-massive-stars-37fgh83sbu.pdf

Binary Interaction Dominates the Evolution of Massive Stars

We present a rapid binary-evolution algorithm that enables modelling of even the most complex binary systems. In addition to all aspects of single-star evolution, features such as mass transfer, mass accretion, common-envelope evolution, collisions, supernova kicks and angular momentum loss mechanisms are included. In particular, circularization and synchronization of the orbit by tidal interactions are calculated for convective, radiative and degenerate damping mechanisms. We use this algorithm to study the formation and evolution of various binary systems. We also investigate the effect that tidal friction has on the outcome of binary evolution. Using the rapid binary code, we generate a series of large binary populations and evaluate the formation rate of interesting individual species and events. By comparing the results for populations with and without tidal friction, we quantify the hitherto ignored systematic effect of tides and show that modelling of tidal evolution in binary systems is necessary in order to draw accurate conclusions from population synthesis work. Tidal synchronism is important but, because orbits generally circularize before Roche lobe overflow, the outcome of the interactions of systems with the same semilatus rectum is almost independent of eccentricity. It is not necessary to include a distribution of eccentricities in population synthesis of interacting binaries; however, the initial separations should be distributed according to the observed distribution of semilatera recta rather than periods or semimajor axes.

https://researchbank.swinburne.edu.au/file/ab4c7d3f-f3a0-4cb3-bb4f-9790d3f33b19/1/PDF (Accepted manuscript).pdf

Evolution of binary stars and the effect of tides on binary populations

The temporal evolution of the optical spectra of various types of supernovae (SNe) is illustrated, in part to aid observers classifying supernova candidates. Type II SNe are defined by the presence of hydrogen, and they exhibit a very wide variety of photometric and spectroscopic properties. Among hydrogen-deficient SNe (Type I), three subclasses are now known: those whose early-time spectra show strong Si II (Ia), prominent He I (Ib), or neither Si II nor He I (Ic). The late-time spectra of SNe Ia consist of a multitude of blended emission lines of iron-group elements; in sharp contrast, those of SNe Ib and SNe Ic (which are similar to each other) are dominated by several relatively unblended lines of intermediatemass elements. Although SNe Ia, which result from the thermonuclear runaway of white dwarfs, constitute a rather homogeneous subclass, important variations in their photometric and spectroscopic properties are undeniably present. SNe Ib/Ic probably result from core collapse in massive stars largely stripped of their hydrogen (Ib) and helium (Ic) envelopes, and hence they are physically related to SNe II. Indeed, the progenitors of some SNe II seem to have only a low-mass skin of hydrogen; their spectra gradually evolve to resemble those of SNe Ib. In addition to the two well-known photometric subclasses (linear and plateau) of SNe II, which may exhibit minor spectroscopic differences, there is a new subclass (SNe IIn) distinguished by relatively narrow emission lines with little or no P Cygni absorption component and slowly declining light curves. These objects probably have unusually dense circumstellar gas with which the ejecta interact.

/pdf/optical-spectra-of-supernovae-2pjopknob7.pdf

Optical spectra of supernovae

The way in which binary interaction affects the presupernova evolution of massive close binaries and the resulting supernova explosions is investigated systematically by means of a Henyey-type stellar evolution code that was modified to allow its application to binary stellar evolution calculations. The code makes it possible to trace the effects of mass and angular momentum loss from the binary, as well as mass transfer within the binary system. It is found that a large number of binary scenarios can be distinguished, depending on the type of binary interaction and the evolutionary stage of the supernova progenitor at the time of the interaction. Monte Carlo simulations are performed to estimate the frequencies of the occurrence of various scenarios. It is found that, because of a previous binary interaction, 15-30 percent of all massive stars (with initial masses greater than about 8 solar masses) become helium stars, and another 5 percent of all massive stars end their lives as blue supergiants rather than as red supergiants.

Presupernova Evolution in Massive Interacting Binaries

The development of appropriate computer programs has made it possible to conduct studies of stellar evolution which are more detailed and accurate than the investigations previously feasible. However, the use of such programs can also entail some serious drawbacks which are related to the time and expense required for the work. One approach for overcoming these drawbacks involves the employment of simplified stellar evolution codes which incorporate the essential physics of the problem of interest without attempting either great generality or maximal accuracy. Rappaport et al. (1982) have developed a simplified code to study the evolution of close binary stellar systems composed of a collapsed object and a low-mass secondary. The present investigation is concerned with a more general, but still simplified, technique for calculating the evolution of close binary systems with collapsed binaries and mass-losing secondaries.

A new technique for calculations of binary stellar evolution, with application to magnetic braking

A new theoretical treatment of the evolution of highly compact binary systems is presented. The evolution is calculated until almost the entire mass of the secondary has been transferred to the primary or lost from the system. It is assumed that gravitational radiation from the system is the cause of mass transfer. It is found that the structure of the mass-losing star can be approximated by an n = 3/2 polytrope, and as a result a relatively large number of different cases can be explored and some general conclusions drawn. An explanation is found for the existence of a cutoff in the orbital period distribution among the cataclysmic variables and light is shed upon the possible generic relationships among cataclysmic variables, the low-mass X-ray binaries, and the spectrally soft transient X-ray sources.

The evolution of highly compact binary stellar systems

The relation between white dwarf mass and orbital period in wide binary radio pulsars

The early history of the development of our understanding of neutron stars has been described by others, and we summarize it only briefly here. Following an early speculation by Baade & Zwicky ( 1934) soon after the discovery of the neutron in 1932, Oppenheimer & Volkoff(1939) presented a classic analysis of the theoretical viability of neutron stars, based on the general theory of relativity and the nuclear physics that was then known. A number of further theoretical studies were published during the next three decades, but it seemed that neutron stars might be virtually unobservable and thus remain little more than a theorist's curiosity piece. In a prophetic paper, Pacini (1967) noted that a rotating neutron star with a magnetic field whose axis was misaligned with the rotation axis should emit intense magnetic dipole radiation. Almost immediately thereafter (but quite independently), Hewish et al. (1968a,b) announced the discovery of radio pulsars. For a brief period, several theoretical models for the pulsar phenomenon were in contention. However, the discovery of the pulsar NP 0532 in the Crab Nebula (Staelin & Reifenstein i968, Comella et al. 1 969), followed by a notably straightforward set of observations and theoretical arguments, soon clinched the case for neutron stars (Gold 1969a).1 (For a review of these remarkable developments, see Gold 1969b.)

Paul C. Joss

Papers

Presupernova Evolution in Massive Interacting Binaries

A new technique for calculations of binary stellar evolution, with application to magnetic braking

The evolution of highly compact binary stellar systems

The relation between white dwarf mass and orbital period in wide binary radio pulsars

Neutron stars in interacting binary systems