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.

We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with M < 8 M become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools − MESA-Web, MESA-Docker, pyMESA, and mesastar.org − to enhance MESA's education and research impact.

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001%. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen-and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development.

No abstract

With recent advances in gravitational-wave astronomy, the direct detection of gravitational waves from the merger of two stellarmass compact objects has become a realistic prospect. Evolutionary scenarios towards mergers of various double compact objects generally invoke so-called common-envelope evolution, which is poorly understood and leads to large uncertainties in the predicted merger rates. Here we explore, as an alternative, the scenario of massive overcontact binary (MOB) evolution, which involves two very massive stars in a very tight binary that remain fully mixed as a result of their tidally induced high spin. While many of these systems merge early on, we find many MOBs that swap mass several times, but survive as a close binary until the stars collapse. The simplicity of the MOB scenario allows us to use the efficient public stellar-evolution code MESA to explore it systematically by means of detailed numerical calculations. We find that, at low metallicity, MOBs produce double-black-hole (BH+BH) systems that will merge within a Hubble time with mass-ratios close to one, in two mass ranges, about 25 . . . 60 M and > ∼ 130 M , with pairinstability supernovae (PISNe) being produced at intermediate masses. Our models are also able to reproduce counterparts of various stages in the MOB scenario in the local Universe, providing direct support for the scenario. We map the initial binary parameter space that produces BH+BH mergers, determine the expected chirp mass distribution, merger times, and expected Kerr parameters, and predict event rates. We find typically one BH+BH merger event for ∼1000 core-collapse supernovae for Z < ∼ Z /10 . The advanced LIGO (aLIGO) detection rate is more uncertain and depends on the cosmic metallicity evolution. From deriving upper and lower limits from a local and a global approximation for the metallicity distribution of massive stars, we estimate aLIGO detection rates (at the aLIGO design limit) of ∼19−550 yr −1 for BH-BH mergers below the PISN gap and of ∼2.1−370 yr −1 above the PISN gap. Even with conservative assumptions, we find that aLIGO will probably soon detect BH+BH mergers from the MOB scenario. These could be the dominant source for aLIGO detections.