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

Optical and X-ray pulsations from Hercules X-1 were observed during the 1977 July/August high state of the source. Some of the optical pulastions apparently originated near the line of sight from the neutron star to the Earth. The X-ray pulse timing observations were sufficiently accurate and adjacent in time to the optical observations that a direct comparison of the optical and X-ray pulse arrival times is possible. If these optical pulsations originate within approx.10/sup 9/ cm of the neutron star, then they are emitted in phase with the X-ray pulsations, to within th measurement errors. However, it seems more likely that the optical pulsations originate in a region of enhanced density near the interaction of the accretion disk with an accretion stream from the companion star, HZ Herculis; this region should have been approximately along the line of sight from the neutron star to the Earth at the time of the observations. We use the X-ray timing data to refine the measured spin-up rate of the neutron star and the orbital parameters of the binary system.

Coordinated observations of optical and X-ray pulsations from Hercules X-1

X-ray pulsations from Her X-1 with energies between 0.1 and 30 keV were observed for four days with the SAS-3 satellite, during the 1975 July-August ON state of the source. The existence of a strong flux between 0.1 and 0.4 keV, with pulsations that are out of phase with those above 1 keV, is confirmed. A pulsed flux in the 19-30 keV band was discovered. The average fractional rate of change in pulse period between 1972 and 1975 was about 3 x 10 to the -6th/yr, and the absolute value of the average fractional rate of change in orbital period during the same interval was not greater than 5 x 10 to the -7th/yr

Observations of Hercules X-1 with SAS-3 during 1975 July.

A search for X-rays from Nova Cygni 1975 was carried out before, during, and after the time of optical maximum. No X-rays were detected over the spectral range 0.1-50 keV. On the basis of these results a strong upper limit of .0001 has been placed on the ratio of X-ray to optical luminosity for this nova, consistent with effective temperatures of about 10,000 K. If Nova Cygni 1975 is a virgin nova, its low mass exchange rate would imply that any associated X-ray emission would not be detectable by present techniques.

Search for X-ray emission from Nova Cygni 1975

A model for the light curve of HZ Herculis

We have calculated the secular evolution of a highly compact binary, composed of a degenerate-dwarf primary and a low-mass main-sequence secondary, taking into account the combined effects of (1) the gradual breakdown of quasi-equilibrium of the /sup 3/He abundance in the interior of the secondary, and (2) the progressive mixture of fresh /sup 3/He into the core of the secondary due to the increasing depth of the stellar surface convection zone. We find that these effects may provide a natural explanation for the apparent sharp decrease in the discovery probability of a cataclysmic variable as its orbital period decreases through approx.3 hours. However, the effects that we consider do not account for the relatively large number of observed cataclysmic variables with orbital periods between approx.80 minutes and approx.2 hours, so that the overall distribution of orbital periods among cataclysmic variables is not yet fully understood.

The breakdown of nuclear quasi-equilibrium in highly compact binaries and the origin of the 2-3 hour gap in the orbital period distribution of cataclysmic variables

Paul C. Joss

Papers

Coordinated observations of optical and X-ray pulsations from Hercules X-1

Observations of Hercules X-1 with SAS-3 during 1975 July.

Search for X-ray emission from Nova Cygni 1975

A model for the light curve of HZ Herculis

The breakdown of nuclear quasi-equilibrium in highly compact binaries and the origin of the 2-3 hour gap in the orbital period distribution of cataclysmic variables