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
- Institución
- Consejo Nacional de Investigaciones Científicas y Técnicas
- OAI Identificador
- oai:ri.conicet.gov.ar:11336/18255
Ver los metadatos del registro completo
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 |