Finite size effects in neutron star and nuclear matter simulations

Autores
Giménez Molinelli, Pedro Agustín; Dorso, Claudio Oscar
Año de publicación
2015
Idioma
inglés
Tipo de recurso
artículo
Estado
versión publicada
Descripción
In this work we study molecular dynamics simulations of symmetric nuclear and neutron star matter using a semi-classical nucleon interaction model. Our aim is to gain insight on the nature of the so-called “finite size effects”, unavoidable in this kind of simulations, and to understand what they actually affect. To do so, we explore different geometries for the periodic boundary conditions imposed on the simulation cell: cube, hexagonal prism and truncated octahedron. For nuclear matter simulations we show that, at sub-saturation densities and low temperatures, the solutions are non-homogeneous structures reminiscent of the “nuclear pasta” phases expected in neutron star matter simulations, but only one structure per cell and shaped by specific artificial aspects of the simulations—for the same physical conditions (i.e. number density and temperature) different cells yield different solutions. The particular shape of the solution at low enough temperature and a given density can be predicted analytically by surface minimization. We also show that even if this behavior is due to the imposition of periodic boundary conditions on finite systems, this does not mean that it vanishes for very large systems, and it is actually independent of the system size. We conclude that, for nuclear matter simulations, the cells' size sets the only characteristic length scale for the inhomogeneities, and the geometry of the periodic cell determines the shape of those inhomogeneities. To model neutron star matter we add a screened Coulomb interaction between protons, and perform simulations in the three cell geometries. Our simulations indeed produce the well known nuclear pasta, with (in most cases) several structures per cell. However, we find that for systems not too large results are affected by finite size in different ways depending on the geometry of the cell. In particular, at the same certain physical conditions and system size, the hexagonal prism yields a single structure per cell while the cubic and truncated octahedron show consistent results, with more than one structure per cell. For systems of the size studied in this work these effects are still noticeable, but we find evidence to support that the dependence of the results on the cell geometry becomes smaller as the system size is increased. When the Coulomb interaction is present, the competition between opposing interactions of different range results in a proper, physically meaningful length scale that is independent of the system size and periodic cell of choice. Only under these conditions “finite size effects” will vanish for large enough systems (i.e. cells much larger than this characteristic length). Larger simulations are in order, but our computational capabilities forbid it for the time being.
Fil: Giménez Molinelli, Pedro Agustín. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; Argentina
Fil: Dorso, Claudio Oscar. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; Argentina
Materia
Nuclear Matter
Neutron Star Matter
Nuclear Pasta
Simulations
Nivel de accesibilidad
acceso abierto
Condiciones de uso
https://creativecommons.org/licenses/by-nc-sa/2.5/ar/
Repositorio
CONICET Digital (CONICET)
Institución
Consejo Nacional de Investigaciones Científicas y Técnicas
OAI Identificador
oai:ri.conicet.gov.ar:11336/18255

id CONICETDig_3be4cc33c363806f7bcb55dbb8d3ad15
oai_identifier_str oai:ri.conicet.gov.ar:11336/18255
network_acronym_str CONICETDig
repository_id_str 3498
network_name_str CONICET Digital (CONICET)
spelling Finite size effects in neutron star and nuclear matter simulationsGiménez Molinelli, Pedro AgustínDorso, Claudio OscarNuclear MatterNeutron Star MatterNuclear PastaSimulationshttps://purl.org/becyt/ford/1.3https://purl.org/becyt/ford/1In this work we study molecular dynamics simulations of symmetric nuclear and neutron star matter using a semi-classical nucleon interaction model. Our aim is to gain insight on the nature of the so-called “finite size effects”, unavoidable in this kind of simulations, and to understand what they actually affect. To do so, we explore different geometries for the periodic boundary conditions imposed on the simulation cell: cube, hexagonal prism and truncated octahedron. For nuclear matter simulations we show that, at sub-saturation densities and low temperatures, the solutions are non-homogeneous structures reminiscent of the “nuclear pasta” phases expected in neutron star matter simulations, but only one structure per cell and shaped by specific artificial aspects of the simulations—for the same physical conditions (i.e. number density and temperature) different cells yield different solutions. The particular shape of the solution at low enough temperature and a given density can be predicted analytically by surface minimization. We also show that even if this behavior is due to the imposition of periodic boundary conditions on finite systems, this does not mean that it vanishes for very large systems, and it is actually independent of the system size. We conclude that, for nuclear matter simulations, the cells' size sets the only characteristic length scale for the inhomogeneities, and the geometry of the periodic cell determines the shape of those inhomogeneities. To model neutron star matter we add a screened Coulomb interaction between protons, and perform simulations in the three cell geometries. Our simulations indeed produce the well known nuclear pasta, with (in most cases) several structures per cell. However, we find that for systems not too large results are affected by finite size in different ways depending on the geometry of the cell. In particular, at the same certain physical conditions and system size, the hexagonal prism yields a single structure per cell while the cubic and truncated octahedron show consistent results, with more than one structure per cell. For systems of the size studied in this work these effects are still noticeable, but we find evidence to support that the dependence of the results on the cell geometry becomes smaller as the system size is increased. When the Coulomb interaction is present, the competition between opposing interactions of different range results in a proper, physically meaningful length scale that is independent of the system size and periodic cell of choice. Only under these conditions “finite size effects” will vanish for large enough systems (i.e. cells much larger than this characteristic length). Larger simulations are in order, but our computational capabilities forbid it for the time being.Fil: Giménez Molinelli, Pedro Agustín. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; ArgentinaFil: Dorso, Claudio Oscar. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; ArgentinaElsevier Science2015-01info:eu-repo/semantics/articleinfo:eu-repo/semantics/publishedVersionhttp://purl.org/coar/resource_type/c_6501info:ar-repo/semantics/articuloapplication/pdfapplication/zipapplication/pdfhttp://hdl.handle.net/11336/18255Giménez Molinelli, Pedro Agustín; Dorso, Claudio Oscar; Finite size effects in neutron star and nuclear matter simulations; Elsevier Science; Nuclear Physics A; 933; 1-2015; 306-3240375-9474CONICET DigitalCONICETenginfo:eu-repo/semantics/altIdentifier/doi/10.1016/j.nuclphysa.2014.11.005info:eu-repo/semantics/altIdentifier/url/http://www.sciencedirect.com/science/article/pii/S0375947414005508info:eu-repo/semantics/altIdentifier/url/https://arxiv.org/abs/1403.5777info:eu-repo/semantics/openAccesshttps://creativecommons.org/licenses/by-nc-sa/2.5/ar/reponame:CONICET Digital (CONICET)instname:Consejo Nacional de Investigaciones Científicas y Técnicas2025-09-29T09:49:46Zoai:ri.conicet.gov.ar:11336/18255instacron:CONICETInstitucionalhttp://ri.conicet.gov.ar/Organismo científico-tecnológicoNo correspondehttp://ri.conicet.gov.ar/oai/requestdasensio@conicet.gov.ar; lcarlino@conicet.gov.arArgentinaNo correspondeNo correspondeNo correspondeopendoar:34982025-09-29 09:49:47.068CONICET Digital (CONICET) - Consejo Nacional de Investigaciones Científicas y Técnicasfalse
dc.title.none.fl_str_mv Finite size effects in neutron star and nuclear matter simulations
title Finite size effects in neutron star and nuclear matter simulations
spellingShingle Finite size effects in neutron star and nuclear matter simulations
Giménez Molinelli, Pedro Agustín
Nuclear Matter
Neutron Star Matter
Nuclear Pasta
Simulations
title_short Finite size effects in neutron star and nuclear matter simulations
title_full Finite size effects in neutron star and nuclear matter simulations
title_fullStr Finite size effects in neutron star and nuclear matter simulations
title_full_unstemmed Finite size effects in neutron star and nuclear matter simulations
title_sort Finite size effects in neutron star and nuclear matter simulations
dc.creator.none.fl_str_mv Giménez Molinelli, Pedro Agustín
Dorso, Claudio Oscar
author Giménez Molinelli, Pedro Agustín
author_facet Giménez Molinelli, Pedro Agustín
Dorso, Claudio Oscar
author_role author
author2 Dorso, Claudio Oscar
author2_role author
dc.subject.none.fl_str_mv Nuclear Matter
Neutron Star Matter
Nuclear Pasta
Simulations
topic Nuclear Matter
Neutron Star Matter
Nuclear Pasta
Simulations
purl_subject.fl_str_mv https://purl.org/becyt/ford/1.3
https://purl.org/becyt/ford/1
dc.description.none.fl_txt_mv In this work we study molecular dynamics simulations of symmetric nuclear and neutron star matter using a semi-classical nucleon interaction model. Our aim is to gain insight on the nature of the so-called “finite size effects”, unavoidable in this kind of simulations, and to understand what they actually affect. To do so, we explore different geometries for the periodic boundary conditions imposed on the simulation cell: cube, hexagonal prism and truncated octahedron. For nuclear matter simulations we show that, at sub-saturation densities and low temperatures, the solutions are non-homogeneous structures reminiscent of the “nuclear pasta” phases expected in neutron star matter simulations, but only one structure per cell and shaped by specific artificial aspects of the simulations—for the same physical conditions (i.e. number density and temperature) different cells yield different solutions. The particular shape of the solution at low enough temperature and a given density can be predicted analytically by surface minimization. We also show that even if this behavior is due to the imposition of periodic boundary conditions on finite systems, this does not mean that it vanishes for very large systems, and it is actually independent of the system size. We conclude that, for nuclear matter simulations, the cells' size sets the only characteristic length scale for the inhomogeneities, and the geometry of the periodic cell determines the shape of those inhomogeneities. To model neutron star matter we add a screened Coulomb interaction between protons, and perform simulations in the three cell geometries. Our simulations indeed produce the well known nuclear pasta, with (in most cases) several structures per cell. However, we find that for systems not too large results are affected by finite size in different ways depending on the geometry of the cell. In particular, at the same certain physical conditions and system size, the hexagonal prism yields a single structure per cell while the cubic and truncated octahedron show consistent results, with more than one structure per cell. For systems of the size studied in this work these effects are still noticeable, but we find evidence to support that the dependence of the results on the cell geometry becomes smaller as the system size is increased. When the Coulomb interaction is present, the competition between opposing interactions of different range results in a proper, physically meaningful length scale that is independent of the system size and periodic cell of choice. Only under these conditions “finite size effects” will vanish for large enough systems (i.e. cells much larger than this characteristic length). Larger simulations are in order, but our computational capabilities forbid it for the time being.
Fil: Giménez Molinelli, Pedro Agustín. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; Argentina
Fil: Dorso, Claudio Oscar. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; Argentina
description In this work we study molecular dynamics simulations of symmetric nuclear and neutron star matter using a semi-classical nucleon interaction model. Our aim is to gain insight on the nature of the so-called “finite size effects”, unavoidable in this kind of simulations, and to understand what they actually affect. To do so, we explore different geometries for the periodic boundary conditions imposed on the simulation cell: cube, hexagonal prism and truncated octahedron. For nuclear matter simulations we show that, at sub-saturation densities and low temperatures, the solutions are non-homogeneous structures reminiscent of the “nuclear pasta” phases expected in neutron star matter simulations, but only one structure per cell and shaped by specific artificial aspects of the simulations—for the same physical conditions (i.e. number density and temperature) different cells yield different solutions. The particular shape of the solution at low enough temperature and a given density can be predicted analytically by surface minimization. We also show that even if this behavior is due to the imposition of periodic boundary conditions on finite systems, this does not mean that it vanishes for very large systems, and it is actually independent of the system size. We conclude that, for nuclear matter simulations, the cells' size sets the only characteristic length scale for the inhomogeneities, and the geometry of the periodic cell determines the shape of those inhomogeneities. To model neutron star matter we add a screened Coulomb interaction between protons, and perform simulations in the three cell geometries. Our simulations indeed produce the well known nuclear pasta, with (in most cases) several structures per cell. However, we find that for systems not too large results are affected by finite size in different ways depending on the geometry of the cell. In particular, at the same certain physical conditions and system size, the hexagonal prism yields a single structure per cell while the cubic and truncated octahedron show consistent results, with more than one structure per cell. For systems of the size studied in this work these effects are still noticeable, but we find evidence to support that the dependence of the results on the cell geometry becomes smaller as the system size is increased. When the Coulomb interaction is present, the competition between opposing interactions of different range results in a proper, physically meaningful length scale that is independent of the system size and periodic cell of choice. Only under these conditions “finite size effects” will vanish for large enough systems (i.e. cells much larger than this characteristic length). Larger simulations are in order, but our computational capabilities forbid it for the time being.
publishDate 2015
dc.date.none.fl_str_mv 2015-01
dc.type.none.fl_str_mv info:eu-repo/semantics/article
info:eu-repo/semantics/publishedVersion
http://purl.org/coar/resource_type/c_6501
info:ar-repo/semantics/articulo
format article
status_str publishedVersion
dc.identifier.none.fl_str_mv http://hdl.handle.net/11336/18255
Giménez Molinelli, Pedro Agustín; Dorso, Claudio Oscar; Finite size effects in neutron star and nuclear matter simulations; Elsevier Science; Nuclear Physics A; 933; 1-2015; 306-324
0375-9474
CONICET Digital
CONICET
url http://hdl.handle.net/11336/18255
identifier_str_mv Giménez Molinelli, Pedro Agustín; Dorso, Claudio Oscar; Finite size effects in neutron star and nuclear matter simulations; Elsevier Science; Nuclear Physics A; 933; 1-2015; 306-324
0375-9474
CONICET Digital
CONICET
dc.language.none.fl_str_mv eng
language eng
dc.relation.none.fl_str_mv info:eu-repo/semantics/altIdentifier/doi/10.1016/j.nuclphysa.2014.11.005
info:eu-repo/semantics/altIdentifier/url/http://www.sciencedirect.com/science/article/pii/S0375947414005508
info:eu-repo/semantics/altIdentifier/url/https://arxiv.org/abs/1403.5777
dc.rights.none.fl_str_mv info:eu-repo/semantics/openAccess
https://creativecommons.org/licenses/by-nc-sa/2.5/ar/
eu_rights_str_mv openAccess
rights_invalid_str_mv https://creativecommons.org/licenses/by-nc-sa/2.5/ar/
dc.format.none.fl_str_mv application/pdf
application/zip
application/pdf
dc.publisher.none.fl_str_mv Elsevier Science
publisher.none.fl_str_mv Elsevier Science
dc.source.none.fl_str_mv reponame:CONICET Digital (CONICET)
instname:Consejo Nacional de Investigaciones Científicas y Técnicas
reponame_str CONICET Digital (CONICET)
collection CONICET Digital (CONICET)
instname_str Consejo Nacional de Investigaciones Científicas y Técnicas
repository.name.fl_str_mv CONICET Digital (CONICET) - Consejo Nacional de Investigaciones Científicas y Técnicas
repository.mail.fl_str_mv dasensio@conicet.gov.ar; lcarlino@conicet.gov.ar
_version_ 1844613538327298048
score 13.070432