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.

We study the diffusion of helium and other heavy elements in the solar interior by solving exactly the set of flow equations developed by Burgers for a multi-component fluid, including the residual heat-flow terms. No approximation is made concerning the relative concentrations and no restriction is placed on the number of elements considered. We give improved diffusion velocities for hydrogen, helium, oxygen and iron, in the analytic form derived previously by Bahcall and Loeb. These expressions for the diffusion velocities are simple to program in stellar evolution codes and are expected to be accurate to ∼ 15%. We find that the inclusion of the residual heat flow terms leads to an increase in the hydrogen diffusion velocity. We compare our numerical results with those obtained analytically by Bahcall and Loeb using a simplified treatment, as well as with those derived numerically by Michaud and Proffitt. We find that for conditions characteristic of the sun, the results of Bahcall and Loeb for the hydrogen diffusion velocity are smaller than our more accurate numerical results by ∼ 30%, except very near the center where the error becomes larger. The Michaud and Proffitt results differ from the numerical results derived here by < ∼ 15%. Our complete treatment of element diffusion can be directly incorporated in a standard stellar evolution code by means of an exportable subroutine, but, for convenience, we also give simple analytical fits to our numerical results.

Photoionization by the high-redshift ultraviolet radiation background heats low density gas before it falls into dark matter potential wells, and it eliminates the neutral hydrogen and singly ionized helium that dominate cooling of primordial gas at temperatures of 10 4 10 5 K. We investigate the in uence of photoionization on galaxy formation using high-resolution simulations with a 1-dimensional, spherically symmetric, Lagrangian hydrodynamics/gravity code. We nd that the presence of a photoionizing background suppresses the formation of galaxies with circular velocities v circ < 30 km s 1 and substantially reduces the mass of cooled baryons in systems with circular velocities up to v circ 50 km s 1 . Above v circ 75 km s 1 , photoionization has no signi cant e ect. Photoionization exerts its in uence primarily by heating gas before collapse; the elimination of line cooling processes is less important. We discuss the implications of these results for hierarchical theories of galaxy formation.

We investigate the early formation of bound objects with masses comparable to the cosmological Jeans mass (10^5 solar masses). We follow the growth of isolated spherically symmetric density peaks starting from the linear perturbative regime. The initial parameters correspond to density peaks of various widths and heights in a Cold Dark Matter cosmology. We use a one-dimensional spherical Lagrangian hydrodynamics code to follow the dynamical, thermal, and non-equilibrium chemical evolution of the gas. The system includes a collisionless dark matter component and a baryonic component composed of the nine species H, H^-, H^+, He, He^+, He^{++}, H_2, H_2^+, and e^-. All relevant chemical reactions between these species and their cooling mechanisms are included in the calculations. We find that radiative cooling by H_2 affects the collapse dynamics of the gas only after it has already virialized and become part of the bound object. Therefore, radiative cooling is unlikely to have triggered the initial collapse of perturbations at redshifts z>10. Nevertheless, objects with baryonic masses well below the linear-theory Jeans mass (<10^3 solar masses) collapse due to shell crossing by the dark matter. Such objects could be the progenitors of a primordial population of high-mass stars in the intergalactic medium.Comment: 40 pages, uuencoded compressed Postscript, 14 figures included as three separate file