Binary Neutron Star Mergers: Mass Ejection, Electromagnetic Counterparts, and Nucleosynthesis

We release dynamical ejecta data from binary neutron star merger simulations.  The outflows are extracted at a fixed coordinate sphere with radius 300 G/c^2 Msun (= 443 km). Only material unbound according to the geodesic criterion is considered to be part of the dynamical ejecta. See [1] for more details.

Included data:

• `Table2.txt`: Table 2 of the paper in machine readable format
• `tabulated_nucsyn.h5`: nucleosynthesis yields from pre-computed parametrized trajectories. The first three indices of each dataset are Ye, entropy, and expansion timescale tau. For example `Y_final[iYe, ientr, itau, iiso]` gives the final abundance of isotope `iiso` with `A[iiso]` and `Z[iiso]` for a trajectory with initial Ye = `Ye[iYe]`, initial entropy `s[ientr]`, and expansion timescale `tau[itau]`.
• `tabulated_rho.h5`: gives the density scale corresponding to the Ye, entropy, and expansion timescale  in `tabulated_nucsyn.h5`.
• `[model].tar`: ejecta data for individual simulations. The naming convention is the same as in the paper.

For each model we provide:

• `outflow.txt`: angle integrated outflow rate and cumulated ejecta mass. Data are given in units with Msun = G = c = 1 (eg, the conversion factor for time to seconds is 4.9258e-6).
• `hist_entropy.dat`: histogram of the ejecta as a function of the entropy (in kb)
• `hist_vinf.dat`: histogram of the ejecta as a function of the asymptotic velocity (in units of c)
• `hist_ye.dat`: histogram of the ejecta as a function of the electron fraction Ye.
• `profile.txt`: time integrated ejecta profiles as a function of the polar angle.
• `hist_vinf_theta.h5`: histograms of the ejecta as a function of the asymptotic velocity and the polar angle.
• `hist_ye_theta.h5`: histograms of the ejecta as a function of the asymptotic velocity and the polar angle.
• `hist_ye_entropy_tau.h5`: histograms of the ejecta as a function of Ye, entropy, and expansion timescale tau.

Additionally we distribute:

• Initial data generated with LORENE and associated EOS tables.
• EOS tables used for the evolution
• Parameter file used for each simulation

For the multidimensional histograms the indices are ordered as specified in the file name, ie the file `hist_ye_theta.h5` tabulates the ejecta mass as a function of Ye (first index) and polar angle theta (second index).

[1] D. Radice, A. Perego, K. Hotokezaka, S. A. Fromm, S. Bernuzzi, and L. F. Roberts, Binary Neutron Star Mergers: Mass Ejection, Electromagnetic Counterparts, and Nucleosynthesis, ApJ 869:130 (2018), arXiv:1809.11161