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

This program was supported by the the Kavli Foundation, Danish National Research Foundation, the Niels Bohr International Academy, and the DARK Cosmology Centre. The UCSC group is supported in part by NSF grant AST-1518052, the Gordon & Betty Moore Foundation, the Heising-Simons Foundation, generous donations from many individuals through a UCSC Giving Day grant, and from fellowships from the Alfred P. Sloan Foundation (R.J.F.), the David and Lucile Packard Foundation (R.J.F. and E.R.) and the Niels Bohr Professorship from the DNRF (E.R.). AMB acknowledges support from a UCMEXUS-CONACYT Doctoral Fellowship. Support for this work was provided by NASA through Hubble Fellowship grants HST-HF-51348.001 (B.J.S.) and HST-HF-51373.001 (M.R.D.) awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. This paper includes data gathered with the 1 meter Swope and 6.5 meter Magellan Telescopes located at Las Campanas Observatory, Chile.r (AGILE) The AGILE Team thanks the ASI management, the technical staff at the ASI Malindi ground station, the technical support team at the ASI Space Science Data Center, and the Fucino AGILE Mission Operation Center. AGILE is an ASI space mission developed with programmatic support by INAF and INFN. We acknowledge partial support through the ASI grant No. I/028/12/2. We also thank INAF, Italian Institute of Astrophysics, and ASI, Italian Space Agency.r (ANTARES) The ANTARES Collaboration acknowledges the financial support of: Centre National de la Recherche Scientifique (CNRS), Commissariat a l'energie atomique et aux energies alternatives (CEA), Commission Europeenne (FEDER fund and Marie Curie Program), Institut Universitaire de France (IUF), IdEx program and UnivEarthS Labex program at Sorbonne Paris Cite (ANR-10-LABX-0023 and ANR-11-IDEX-0005-02), Labex OCEVU (ANR-11-LABX-0060) and the A*MIDEX project (ANR-11-IDEX-0001-02), Region Ile-de-France (DIM-ACAV), Region Alsace (contrat CPER), Region Provence-Alpes-Cite d'Azur, Departement du Var and Ville de La Seyne-sur-Mer, France; Bundesministerium fur Bildung und Forschung (BMBF), Germany; Istituto Nazionale di Fisica Nucleare (INFN), Italy; Nederlandse organisatie voor Wetenschappelijk Onderzoek (NWO), the Netherlands; Council of the President of the Russian Federation for young scientists and leading scientific schools supporting grants, Russia; National Authority for Scientific Research (ANCS), Romania;...

/pdf/multi-messenger-observations-of-a-binary-neutron-star-12ssh2a28t.pdf

Multi-messenger Observations of a Binary Neutron Star

On 2019 April 25, the LIGO Livingston detector observed a compact binary coalescence with signal-to-noise ratio 12.9. The Virgo detector was also taking data that did not contribute to detection due to a low signal-to-noise ratio, but were used for subsequent parameter estimation. The 90% credible intervals for the component masses range from to if we restrict the dimensionless component spin magnitudes to be smaller than 0.05). These mass parameters are consistent with the individual binary components being neutron stars. However, both the source-frame chirp mass and the total mass of this system are significantly larger than those of any other known binary neutron star (BNS) system. The possibility that one or both binary components of the system are black holes cannot be ruled out from gravitational-wave data. We discuss possible origins of the system based on its inconsistency with the known Galactic BNS population. Under the assumption that the signal was produced by a BNS coalescence, the local rate of neutron star mergers is updated to 250-2810.

/pdf/gw190425-observation-of-a-compact-binary-coalescence-with-lyej4fjfjg.pdf

GW190425: Observation of a Compact Binary Coalescence with Total Mass ∼ 3.4 M O

We report the observation of a compact binary coalescence involving a 22.2 -
24.3 $M_{\odot}$ black hole and a compact object with a mass of 2.50 - 2.67
$M_{\odot}$ (all measurements quoted at the 90$\%$ credible level). The
gravitational-wave signal, GW190814, was observed during LIGO's and Virgo's
third observing run on August 14, 2019 at 21:10:39 UTC and has a
signal-to-noise ratio of 25 in the three-detector network. The source was
localized to 18.5 deg$^2$ at a distance of $241^{+41}_{-45}$ Mpc; no
electromagnetic counterpart has been confirmed to date. The source has the most
unequal mass ratio yet measured with gravitational waves,
$0.112^{+0.008}_{-0.009}$, and its secondary component is either the lightest
black hole or the heaviest neutron star ever discovered in a double
compact-object system. The dimensionless spin of the primary black hole is
tightly constrained to $\leq 0.07$. Tests of general relativity reveal no
measurable deviations from the theory, and its prediction of higher-multipole
emission is confirmed at high confidence. We estimate a merger rate density of
1-23 Gpc$^{-3}$ yr$^{-1}$ for the new class of binary coalescence sources that
GW190814 represents. Astrophysical models predict that binaries with mass
ratios similar to this event can form through several channels, but are
unlikely to have formed in globular clusters. However, the combination of mass
ratio, component masses, and the inferred merger rate for this event challenges
all current models for the formation and mass distribution of compact-object
binaries.

GW190814: Gravitational Waves from the Coalescence of a 23 M$_\odot$ Black Hole with a 2.6 M$_\odot$ Compact Object

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

Thus far, judging the fate of a massive star (either a neutron star (NS) or a black hole) solely by its structure prior to core collapse has been ambiguous. Our work and previous attempts find a non-monotonic variation of successful and failed supernovae with zero-age main-sequence mass, for which no single structural parameter can serve as a good predictive measure. However, we identify two parameters computed from the pre-collapse structure of the progenitor, which in combination allow for a clear separation of exploding and non-exploding cases with only few exceptions ( 1‐2.5%) in our set of 621 investigated stellar models. One parameter is M4, defining the enclosed mass for a dimensionless entropy per nucleon of s = 4, and the other is 4 dm=drjs=4, being the mass-derivative at this location. The two parameters 4 and M4 4 can be directly linked to the mass-infall rate, ˙ M, of the collapsing star and the electron-type neutrino luminosity of the accreting proto-NS, L e / Mns ˙ M, which play a crucial role in the “critical luminosity” concept for the theoretical description of neutrino-driven explosions as runaway phenomenon of the stalled accretion shock. All models were evolved employing the approach of Ugliano et al. for simulating neutrino-driven explosions in spherical symmetry. The neutrino emission of the accretion layer is approximated by a gray transport solver, while the uncertain neutrino emission of the 1.1 M proto-NS core is parametrized by an analytic model. The free parameters connected to the core-boundary prescription are calibrated to reproduce the observables of Supernova 1987A for five di erent progenitor models. Subject headings: supernovae: general — stars: massive — hydrodynamics — neutrinos

http://iopscience.iop.org/article/10.3847/0004-637X/818/2/124/pdf

A two-parameter criterion for classifying the explodability of massive stars by the neutrino-driven mechanism

Thus far, judging the fate of a massive star (either a neutron star (NS) or a black hole) solely by its structure prior to core collapse has been ambiguous. Our work and previous attempts find a non-monotonic variation of successful and failed supernovae with zero-age main-sequence mass, for which no single structural parameter can serve as a good predictive measure. However, we identify two parameters computed from the pre-collapse structure of the progenitor, which in combination allow for a clear separation of exploding and non-exploding cases with only few exceptions (~1-2.5%) in our set of 621 investigated stellar models. One parameter is M4, defining the normalized enclosed mass for a dimensionless entropy per nucleon of s=4, and the other is mu4 = d(m/M_sun)/d(r/1000 km) at s=4, being the normalized mass-derivative at this location. The two parameters mu4 and M4*mu4 can be directly linked to the mass-infall rate, Mdot, of the collapsing star and the electron-type neutrino luminosity of the accreting proto-NS, L_nue ~ M_ns*Mdot, which play a crucial role in the "critical luminosity" concept for the theoretical description of neutrino-driven explosions as runaway phenomenon of the stalled accretion shock. All models were evolved employing the approach of Ugliano et al. for simulating neutrino-driven explosions in spherical symmetry. The neutrino emission of the accretion layer is approximated by a gray transport solver, while the uncertain neutrino emission of the 1.1 M_sun proto-NS core is parametrized by an analytic model. The free parameters connected to the core-boundary prescription are calibrated to reproduce the observables of Supernova 1987A for five different progenitor models.

/pdf/a-two-parameter-criterion-for-classifying-the-explodability-3dysb7fzka.pdf

T. Ertl

Papers

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

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

A two-parameter criterion for classifying the explodability of massive stars by the neutrino-driven mechanism

A two-parameter criterion for classifying the explodability of massive stars by the neutrino-driven mechanism

The Explosion of Helium Stars Evolved with Mass Loss