TY - JOUR
T1 - Pasta nucleosynthesis
T2 - Molecular dynamics simulations of nuclear statistical equilibrium
AU - Caplan, M. E.
AU - Schneider, A. S.
AU - Horowitz, C. J.
AU - Berry, D. K.
N1 - Publisher Copyright:
© 2015 American Physical Society.
PY - 2015/6/22
Y1 - 2015/6/22
N2 - Background: Exotic nonspherical nuclear pasta shapes are expected in nuclear matter at just below saturation density because of competition between short-range nuclear attraction and long-range Coulomb repulsion. Purpose: We explore the impact nuclear pasta may have on nucleosynthesis during neutron star mergers when cold dense nuclear matter is ejected and decompressed. Methods: We use a hybrid CPU/GPU molecular dynamics (MD) code to perform decompression simulations of cold dense matter with 51200 and 409600 nucleons from 0.080fm-3 down to 0.00125fm-3. Simulations are run for proton fractions YP= 0.05, 0.10, 0.20, 0.30, and 0.40 at temperatures T= 0.5, 0.75, and 1.0 MeV. The final composition of each simulation is obtained using a cluster algorithm and compared to a constant density run. Results: Size of nuclei in the final state of decompression runs are in good agreement with nuclear statistical equilibrium (NSE) models for temperatures of 1MeV while constant density runs produce nuclei smaller than the ones obtained with NSE. Our MD simulations produces unphysical results with large rod-like nuclei in the final state of T=0.5MeV runs. Conclusions: Our MD model is valid at higher densities than simple nuclear statistical equilibrium models and may help determine the initial temperatures and proton fractions of matter ejected in mergers.
AB - Background: Exotic nonspherical nuclear pasta shapes are expected in nuclear matter at just below saturation density because of competition between short-range nuclear attraction and long-range Coulomb repulsion. Purpose: We explore the impact nuclear pasta may have on nucleosynthesis during neutron star mergers when cold dense nuclear matter is ejected and decompressed. Methods: We use a hybrid CPU/GPU molecular dynamics (MD) code to perform decompression simulations of cold dense matter with 51200 and 409600 nucleons from 0.080fm-3 down to 0.00125fm-3. Simulations are run for proton fractions YP= 0.05, 0.10, 0.20, 0.30, and 0.40 at temperatures T= 0.5, 0.75, and 1.0 MeV. The final composition of each simulation is obtained using a cluster algorithm and compared to a constant density run. Results: Size of nuclei in the final state of decompression runs are in good agreement with nuclear statistical equilibrium (NSE) models for temperatures of 1MeV while constant density runs produce nuclei smaller than the ones obtained with NSE. Our MD simulations produces unphysical results with large rod-like nuclei in the final state of T=0.5MeV runs. Conclusions: Our MD model is valid at higher densities than simple nuclear statistical equilibrium models and may help determine the initial temperatures and proton fractions of matter ejected in mergers.
UR - http://www.scopus.com/inward/record.url?scp=84936804610&partnerID=8YFLogxK
U2 - 10.1103/PhysRevC.91.065802
DO - 10.1103/PhysRevC.91.065802
M3 - ???researchoutput.researchoutputtypes.contributiontojournal.article???
AN - SCOPUS:84936804610
SN - 0556-2813
VL - 91
JO - Physical Review C - Nuclear Physics
JF - Physical Review C - Nuclear Physics
IS - 6
M1 - 065802
ER -