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 at T = 6 GK corresponding to the Ye, entropy, and expansion timescale used 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</p