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

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

We present analytic formulae that approximate the evolution of stars for a wide range of mass M and metallicity Z Stellar luminosity, radius and core mass are given as a function of age, M and Z, for all phases from the zero-age main sequence up to, and including, the remnant stages For the most part we find continuous formulae accurate to within 5 per cent of detailed models These formulae are useful for purposes such as population synthesis that require very rapid but accurate evaluation of stellar properties, and in particular for use in combination with N-body codes We describe a mass-loss prescription that can be used with these formulae, and investigate the resulting stellar remnant distribution

https://researchbank.swinburne.edu.au/file/9b350da3-a16d-47fc-927b-4cdb1a8aa58f/1/PDF (Accepted manuscript).pdf

Comprehensive analytic formulae for stellar evolution as a function of mass and metallicity

Modules for Experiments in Stellar Astrophysics (MESA): Binaries, Pulsations, and Explosions

The major uncertainties involved in the Chandrasekhar mass models for Type Ia supernovae (SNe Ia) are related to the companion star of their accreting white dwarf progenitor (which determines the accretion rate and consequently the carbon ignition density) and the flame speed after the carbon ignition. We calculate explosive nucleosynthesis in relatively slow deflagrations with a variety of deflagration speeds and ignition densities to put new constraints on the above key quantities. The abundance of the Fe group, in particular of neutron-rich species like 48Ca,50Ti,54Cr,54,58Fe, and 58Ni, is highly sensitive to the electron captures taking place in the central layers. The yields obtained from such a slow central deflagration, and from a fast deflagration or delayed detonation in the outer layers, are combined and put to comparison with solar isotopic abundances. To avoid excessively large ratios of 54Cr/56Fe and 50Ti/56Fe, the central density of the average white dwarf progenitor at ignition should be as low as 2 ? 109 g cm-3. To avoid the overproduction of 58Ni and 54Fe, either the flame speed should not exceed a few percent of the sound speed in the central low Ye layers or the metallicity of the average progenitors has to be lower than solar. Such low central densities can be realized by a rapid accretion as fast as -->img1.gif 1 ? 10-7 M? yr-1. In order to reproduce the solar abundance of 48Ca, one also needs progenitor systems that undergo ignition at higher densities. Even the smallest laminar flame speeds after the low-density ignitions would not produce sufficient amount of this isotope. We also found that the total amount of 56Ni, the Si-Ca/Fe ratio, and the abundance of some elements like Mn and Cr (originating from incomplete Si burning), depend on the density of the deflagration-detonation transition in delayed detonations. Our nucleosynthesis results favor transition densities slightly below 2.2 ? 107 g cm-3.

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

Nucleosynthesis in Chandrasekhar Mass Models for Type Ia Supernovae and Constraints on Progenitor Systems and Burning-Front Propagation

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

A model of bright, low-mass X-ray binaries is proposed which features a lower giant-branch star losing mass on a nuclear time scale to an accreting compact companion. Simple numerical models show that mass transfer rates > or =10/sup -9/ M/sub sun/ yr/sup -1/ are sustained at very nearly a constant rate until the envelope of the donor star is exhausted. The model predicts orbital periods in the range 1/sup d/-200/sup d/ and X-ray to optical luminosity ratios L/sub x//L/sub opt/roughly-equal200-1000 for these sources. It accounts in a natural way for the large fraction of the total galactic bulge luminosity emitted by a few bright (> or =10/sup 37/ ergs s/sup -1/) sources. It also accords very well with the observed X-ray and optical properties of the halo source Cyg X-2 and also with those of 2S 0921-63, provided this latter system contains a massive accreting white dwarf rather than a neutron star. Problems of the prior evolution of low-mass X-ray sources are also briefly delineated.

Evolutionary status of bright, low-mass x-ray sources/sup 1/

Le temps de croissance d'une instabilite de transfert de masse dans un systeme binaire serre de faible masse est calcule. La variation du rayon du lobe de Roche due au transfert de masse est aussi determinee. Les resultats obtenus sont appliques au cas de formation d'un pulsar milliseconde et a celui du transfert de masse naissant d'une geante vers une etoile a neutrons

Mass Transfer Instabilities Due to Angular Momentum Flows in Close Binaries

Luminous supersoft X-ray sources, with characteristic luminosities of approximately 10(exp 38) ergs/s and temperatures, kT, of approximately 35 eV, have been established as a new and distinct class of X-ray source through recent Roentgen Satellite (ROSAT) observations. Several possible physical models have been proposed for these sources. One promising scenario (van den Heuvel et al. 1992) involves mass transfer, which is unstable on a thermal timescale, from a main-sequence or subgiant donor star onto the surface of a white dwarf. For a narrow range of accretion rates, steady nuclear burning of the accreted matter can take place. This process can provide the high luminosities and the correct range of temperatures observed in the supersoft sources. However, given the limited range of mass transfer rates that are consistent with this phenomenon, it is far from obvious that a sufficient population of such systems exists in galaxies such as our own, M31, and the Magellanic Clouds, in order to account for the large number of supersoft sources which can be inferred from present observations. This work addresses the population question in detail, through a Monte Carlo simulation of the formation and evolution of such systems, which starts with zero-age primordial binaries. In order to evolve into close binary systems which contain a white dwarf component and a companion transferring mass at a rate within the requisite narrow range, a binary system must undergo a specific progression of evolutionary steps. We find that a sufficient subset of our initial binaries evolve to become systems with the requisite properties, so that they can account for the population of supersoft sources that is inferred from observations. In particular, we find that there should be more than 1000 systems in the Galaxy today with properties that very closely match those of the observed supersoft sources. From our models, we find expected luminosities, white dwarf effective temperatures, and orbital periods in the ranges of 10(exp 37) - 10(exp 38) ergs/s, (1-5) x 10(exp 5) K, and 8 hr-1.4 days, respectively. The masses of the white dwarf and donor star are expected to lie in the range of 0.7-1.05 solar mass and 1.3-2.7 solar mass. Finally, we discuss the role that supersoft X-ray sources may play as the progenitors of Type Ia supernovae, and estimate the rate of production via this channel.

Saul Rappaport

Papers

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

The evolution of highly compact binary stellar systems

Evolutionary status of bright, low-mass x-ray sources/sup 1/

Mass Transfer Instabilities Due to Angular Momentum Flows in Close Binaries

Formation and evolution of luminous supersoft X-ray sources