Cococubed.com MESA
Home Astronomy research   Software Infrastructure:      MESA      FLASH      STARLIB      MESA-Web      starkiller-astro      My instruments   White dwarf supernova:      Stable nickel production      Remnant metallicities      Colliding white dwarfs      Merging white dwarfs      Ignition conditions      Metallicity effects      Central density effects      Detonation density effects      Tracer particle burning      Subsonic burning fronts      Supersonic burning fronts      W7 profiles   Massive star supernova:      Rotating progenitors      3D evolution      26Al & 60Fe      44Ti, 60Co & 56Ni      Yields of radionuclides      Effects of 12C +12C      SN 1987A light curve      Constraints on Ni/Fe ratios      An r-process   Neutron Stars and Black Holes:      Black Hole Mass Gap      Compact object IMF   Stars:      Variable white dwarfs      Pop III with JWST      Micronovae      Neutrino HR diagram      Monte Carlo massive stars      Pre-supernovaneutrinos      Pre-supernova variations      Monte Carlo white dwarfs      SAGB stars      Nugrid Yields I      Classical novae      He shell convection      Presolar grains      He burn on neutron stars      BBFH at 40 years   Chemical Evolution:      Iron Pseudocarbynes      Radionuclides in the 2020s      Hypatia catalog      Zone models H to Zn      Mixing ejecta      γ-rays within 100 Mpc   Thermodynamics & Networks      Stellar EOS      12C(α,γ)16O Rate      Proton-rich NSE      Reaction networks      Bayesian reaction rates   Verification Problems:      Validating an astro code      Su-Olson      Cog8      Mader      RMTV      Sedov      Noh Software instruments Presentations Illustrations cococubed YouTube Bicycle adventures Public Outreach Education materials 2022 ASU Solar Systems Astronomy 2022 ASU Energy in Everyday Life AAS Journals AAS YouTube 2022 Earendel, A Highly Magnified Star 2022 TV Columbae, Micronova 2022 White Dwarfs and 12C(α,γ)16O 2022 MESA in Don't Look Up 2022 MESA Marketplace 2022 MESA Summer School 2022 MESA Classroom 2021 Bill Paxton, Tinsley Prize Contact: F.X.Timmes my one page vitae, full vitae, research statement, and teaching statement.	Modules for Experiments in Stellar Astrophysics (MESA ⅤⅠ): Time-Dependent Convection, Energy Conservation, Automatic Differentiation, and Infrastructure - 2022 In this ApJ Supplement submission, we update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The new auto_diff module implements automatic differentiation in MESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite difference approximations. We significantly enhance the treatment of the growth and decay of convection in MESA with a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron degenerate ignition events. We strengthen MESA's implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars in MESA, we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator split nuclear burning mode. We close by discussing major updates to MESA's software infrastructure that enhance source code development and community engagement. TDC schematic EOS regions Tests Submitted MESA appears in the Netflix movie "Don't Look Up" - 2022 First noted by Jared Goldberg on his Twitter feed, a lovely ApJ Supplement figure makes an appearance at about the 1 hr 11 min mark: Here is the relevant page in the movie's mashup "report": How and Why a MESA III figure appears in this movie is a story itself. Thank you Amy! Tinsley Prize - 2021 Bill Paxton was awarded the 2021 Beatrice M. Tinsley Prize "for his inspired work on providing, maintaining, and supporting the use of open-source stellar-evolution codes that have seeped into the foundation of research and education efforts and have given rise to an immense amount of new research across multiple subfields of astrophysics." The Tinsley Prize, awarded every two years by the American Astronomical Society, recognizes an outstanding research contribution to astronomy or astrophysics of an exceptionally creative or innovative character. Thank you Bill! Modules for Experiments in Stellar Astrophysics (MESA Ⅴ): Pulsating Variable Stars, Rotation, Convective Boundaries, and Energy Conservation - 2019 In this paper, 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 nonlinear 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, as well as 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. Nonlinear pulsations He flash and gold tolerances Equipotentials of Roche potential Grad profiles with CPM Modules for Experiments in Stellar Astrophysics (MESA Ⅳ): Convective Boundaries, Element Diffusion, and Massive Star Explosions - 2018 In this paper, 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 Msun 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. Sedov Entropy Trail H shell flash on He WD Modules for Experiments in Stellar Astrophysics (MESA Ⅲ): Binaries, Pulsations, and Explosions - 2015 In this paper, 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. 16 M☉ + 14 M☉ binary Instability strips 1.57×1051 erg explosion Modules for Experiments in Stellar Astrophysics (MESA Ⅱ): Giant Planets, Oscillations, Rotation, and Massive Stars - 2013 In this paper, we substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESAstar. Improvements in MESAstar's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESAstar. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESAstar calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 Msun stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESAstar solves the fully coupled stellar structure and composition equations, and we show how this has improved MESA's performance scaling on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. Giant planet evolution Echelle diagram of HD43985. Magnetic field and angular momentum Modules For Experiments In Stellar Astrophysics (MESA Ⅰ) - 2011 In this paper, we describe Modules for Experiments in Stellar Astrophysics (MESA), a suite of open source, robust, efficient, thread-safe libraries for a wide range of applications in computational stellar astrophysics. A one-dimensional stellar evolution module, MESAstar, combines many of the numerical and physics modules for simulations of a wide range of stellar evolution scenarios ranging from very low mass to massive stars, including advanced evolutionary phases. MESAstar solves the fully coupled structure and composition equations simultaneously. It uses adaptive mesh refinement and sophisticated timestep controls, and supports shared memory parallelism based on OpenMP. State-of-the-art modules provide equation of state, opacity, nuclear reaction rates, element diffusion data, and atmosphere boundary conditions. Each module is constructed as a separate Fortran 95 library with its own explicitly defined public interface to facilitate independent development. Several detailed examples indicate the extensive verification and testing that is continuously performed and demonstrate the wide range of capabilities that MESA possesses. MESA can be downloaded from the project Web site https://docs.mesastar.org Equation of State Opacity Massive stars in ρ-T plane