Initial-value ordinary differential equation solution via variable order Adams method (SIVA/DIVA) package is collection of subroutines for solution of nonstiff ordinary differential equations. There are versions for single-precision and double-precision arithmetic. Requires fewer evaluations of derivatives than other variable-order Adams predictor/ corrector methods. Option for direct integration of second-order equations makes integration of trajectory problems significantly more efficient. Written in FORTRAN 77.

Solving Ordinary Differential Equations

Numerical Initial Value Problems in Ordinary Differential Equations

Nucleosynthesis, light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive stars with solar metallicity and masses from 9.0 to 120 solar masses. The full evolution is followed using an adaptive reaction network of up to 2000 nuclei. A novel aspect of the survey is the use of a one-dimensional neutrino transport model for the explosion. This explosion model has been calibrated to give the observed energy for SN 1987A, using several standard progenitors, and for the Crab supernova using a 9.6 solar mass progenitor. As a result of using a calibrated central engine, the final kinetic energy of the supernova is variable and sensitive to the structure of the presupernova star. Many progenitors with extended core structures do not explode, but become black holes, and the masses of exploding stars do not form a simply connected set. The resulting nucleosynthesis agrees reasonably well with the sun provided that a reasonable contribution from Type Ia supernovae is also allowed, but with a deficiency of light s-process isotopes. The resulting neutron star IMF has a mean gravitational mass near 1.4 solar masses. The average black hole mass is about 9 solar masses if only the helium core implodes, and 14 solar masses if the entire presupernova star collapses. Only ~10% of supernovae come from stars over 20 solar masses and some of these are Type Ib or Ic. Some useful systematics of Type IIp light curves are explored.

https://iopscience.iop.org/article/10.3847/0004-637X/821/1/38/pdf

Core-Collapse Supernovae from 9 to 120 Solar Masses Based on Neutrino-powered Explosions

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  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  radiation transfer instrument, creates new avenues for exploring Type II supernova properties. These capabilities are exhibited with exploratory models of pair-instability supernovae, pulsational pair-instability supernovae, and the formation of stellar-mass black holes. The applicability of MESA is now widened by the capability to import multidimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, as well as four new software tools—  ,  -Docker,  , and mesastar.org—to enhance MESA's education and research impact.

https://iopscience.iop.org/article/10.3847/1538-4365/aaa5a8/pdf

Modules for Experiments in Stellar Astrophysics (MESA) : Convective Boundaries, Element Diffusion, and Massive Star Explosions

The well-studied type IIP SN 1987A, produced by the explosion of a blue supergiant (BSG) star, is a touchstone for massive-star evolution, simulations of neutrino-driven explosions, and modeling of light curves and spectra. In the framework of the neutrino-driven mechanism, we study the dependence of explosion properties on the structure of four different BSGs and compare the corresponding light curves with observations of SN 1987A. We perform 3D simulations with the PROMETHEUS code until about one day and map the results to the 1D code CRAB for the light curve calculations. All of our 3D models with explosion energies compatible with SN 1987A produce 56Ni in rough agreement with the amount deduced from fitting the radioactively powered light-curve tail. One of the progenitors yields maximum velocities of ~3000 km/s for the bulk of ejected 56Ni, consistent with observations. In all of our models inward mixing of hydrogen during the 3D evolution leads to minimum H-velocities below 100 km/s, in good agreement with spectral observations. The considered BSG models, 3D explosion simulations, and light-curve calculations can thus explain basic observational features of SN 1987A. However, all progenitors have too large pre-SN radii to reproduce the narrow initial luminosity peak, and the structure of their outer layers is not suitable to match the observed light curve during the first 30-40 days. Only one stellar model has a structure of the He core and the He/H composition interface that enables sufficient outward mixing of 56Ni and inward mixing of hydrogen to produce a good match of the dome-like shape of the observed light-curve maximum. But this model falls short of the He-core mass of 6 Msun inferred from the absolute luminosity of the pre-SN star. The lack of an adequate pre-SN model for SN 1987A is a pressing challenge for the theory of massive-star evolution. (Abridged)

/pdf/supernova-1987a-neutrino-driven-explosions-in-three-29dxrv7z52.pdf

Supernova 1987A: neutrino-driven explosions in three dimensions and light curves

Type II-plateau supernovae (SNe IIP) are the most numerous subclass of core-collapse SNe originating from massive stars. In the framework of the neutrino-driven explosion mechanism, we study the SN outburst properties for a red supergiant progenitor model and compare the corresponding light curves with observations of the ordinary Type IIP SN 1999em. Three-dimensional (3D) simulations of (parametrically triggered) neutrino-driven explosions are performed with the (explicit, finite-volume, Eulerian, multifluid hydrodynamics) code PROMETHEUS, using a presupernova model of a 15 Msun star as initial data. At approaching homologous expansion, the hydrodynamical and composition variables of the 3D models are mapped to a spherically symmetric configuration, and the simulations are continued with the (implicit, Lagrangian radiation-hydrodynamics) code CRAB to follow the blast-wave evolution during the SN outburst. Our 3D neutrino-driven explosion model with an explosion energy of about 0.5x10^51 erg produces Ni-56 in rough agreement with the amount deduced from fitting the radioactively powered light-curve tail of SN 1999em. The considered presupernova model, 3D explosion simulations, and light-curve calculations can explain the basic observational features of SN 1999em, except for those connected to the presupernova structure of the outer stellar layers. Our 3D simulations show that the distribution of Ni-rich matter in velocity space is asymmetric with a strong dipole component that is consistent with the observations of SN 1999em. The monotonic luminosity decline from the plateau to the radioactive tail in ordinary SNe IIP is a manifestation of the intense turbulent mixing at the He/H composition interface.

Light curve analysis of ordinary type IIP supernovae based on neutrino-driven explosion simulations in three dimensions

With the same method as used previously, we investigate neutrino-driven explosions of a larger sample of blue supergiant models. The larger sample includes three new presupernova stars. The results are compared with light-curve observations of the peculiar type IIP SN 1987A. The explosions were modeled in 3D with the neutrino-hydrodynamics code PROMETHEUS-HOTB, and light-curve calculations were performed in spherical symmetry with the radiation-hydrodynamics code CRAB. Our results confirm the basic findings of the previous work: 3D neutrino-driven explosions with SN 1987A-like energies synthesize an amount of Ni-56 that is consistent with the radioactive tail of the light curve. Moreover, the models mix hydrogen inward to minimum velocities below 400 km/s as required by spectral observations. Hydrodynamic simulations with the new progenitor models, which possess smaller radii than the older ones, show much better agreement between calculated and observed light curves in the initial luminosity peak and during the first 20 days. A set of explosions with similar energies demonstrated that a high growth factor of Rayleigh-Taylor instabilities at the (C+O)/He composition interface combined with a weak interaction of fast Rayleigh-Taylor plumes, where the reverse shock occurs below the He/H interface, provides a sufficient condition for efficient outward mixing of Ni-56 into the hydrogen envelope. This condition is realized to the required extent only in one of the older stellar models, which yielded a maximum velocity of around 3000 km/s for the bulk of ejected Ni-56, but failed to reproduce the helium-core mass of 6 Msun inferred from the absolute luminosity of the presupernova star. We conclude that none of the single-star progenitor models proposed for SN 1987A to date satisfies all constraints set by observations. (Abridged)

/pdf/three-dimensional-mixing-and-light-curves-constraints-on-the-10o2wvfh95.pdf

Three-dimensional mixing and light curves: constraints on the progenitor of supernova 1987A

Six binary-merger progenitors of Supernova 1987A (SN 1987A) with properties close to those of the blue supergiant Sanduleak -69 202 are exploded by neutrino heating and evolved until long after shock breakout in three dimensions (3D), and continued for light-curve calculations in spherical symmetry. Our results confirm previous findings for single-star progenitors: (1) 3D neutrino-driven explosions with SN 1987A-like energies synthesize Ni-56 masses consistent with the radioactive light-curve tail; (2) hydrodynamic models mix hydrogen inward to minimum velocities below 40 km/s compatible with spectral observations of SN 1987A; and (3) for given explosion energy the efficiency of outward radioactive Ni-56 mixing depends mainly on high growth factors of Rayleigh-Taylor instabilities at the (C+O)/He and He/H composition interfaces and a weak interaction of fast plumes with the reverse shock occurring below the He/H interface. All binary-merger models possess presupernova radii matching the photometric radius of Sanduleak -69 202 and a structure of the outer layers allowing them to reproduce the observed initial luminosity peak in the first about 7 days. Models that mix about 0.5 Msun of hydrogen into the He-shell and exhibit strong outward mixing of Ni-56 with maximum velocities exceeding the 3000 km/s observed for the bulk of ejected Ni-56 have light-curve shapes in good agreement with the dome of the SN 1987A light curve. A comparative analysis of the best representatives of our 3D neutrino-driven explosion models of SN 1987A based on single-star and binary-merger progenitors reveals that only one binary model fulfills all observational constraints, except one.

V. P. Utrobin

Papers

Supernova 1987A: neutrino-driven explosions in three dimensions and light curves

Light curve analysis of ordinary type IIP supernovae based on neutrino-driven explosion simulations in three dimensions

Three-dimensional mixing and light curves: constraints on the progenitor of supernova 1987A

Supernova 1987A: 3D Mixing and light curves for explosion models based on binary-merger progenitors.