A new method to find the potential center of N-body systems

Autores
Aguilar, L. A.; Cruz, F.; Carpintero, Daniel Diego
Año de publicación
2002
Idioma
inglés
Tipo de recurso
artículo
Estado
versión publicada
Descripción
We present a new and fast method to nd the potential center of an N-body distribution. The method uses an iterative algorithm which exploits the fact that the gradient of the potential is null at its center: it uses a smoothing radius to avoid getting trapped in secondary minima. We have tested this method on several random realizations of King models (in which the numerical computation of this center is rather dicult, due to the constant density within their cores), and com- pared its performance and accuracy against a more straightforward, but computer intensive method, based on cartesian meshes of increasing spatial resolution. In all cases, both methods converged to the same center, within the mesh resolution, but the new method is two orders of magnitude faster. We have also tested the method with one astronomical problem: the evolu- tion of a 105 particle King model orbiting around a xed potential that represents our Galaxy. We used a spherical harmonics expansion N-body code, in which the potential center determination is crucial for the correct force computation. We compared this simulation with another one in which a method previously used to determine the expansion center is employed (White 1983). Our routine gives better results in energy conservation and mass loss.
Fil: Aguilar, L. A.. Universidad Nacional Autonoma de Mexico. Instituto de Astronomia; México
Fil: Cruz, F.. Universidad Nacional Autonoma de Mexico. Instituto de Astronomia; México
Fil: Carpintero, Daniel Diego. Universidad Nacional de la Plata. Facultad de Ciencias Astronómicas y Geofísicas; Argentina
Materia
Globular clusters
Numerical methods
N Body 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/36956

id CONICETDig_be9c724efd3e97f457be3f370b50122d
oai_identifier_str oai:ri.conicet.gov.ar:11336/36956
network_acronym_str CONICETDig
repository_id_str 3498
network_name_str CONICET Digital (CONICET)
spelling A new method to find the potential center of N-body systemsAguilar, L. A.Cruz, F.Carpintero, Daniel DiegoGlobular clustersNumerical methodsN Body simulationsWe present a new and fast method to nd the potential center of an N-body distribution. The method uses an iterative algorithm which exploits the fact that the gradient of the potential is null at its center: it uses a smoothing radius to avoid getting trapped in secondary minima. We have tested this method on several random realizations of King models (in which the numerical computation of this center is rather dicult, due to the constant density within their cores), and com- pared its performance and accuracy against a more straightforward, but computer intensive method, based on cartesian meshes of increasing spatial resolution. In all cases, both methods converged to the same center, within the mesh resolution, but the new method is two orders of magnitude faster. We have also tested the method with one astronomical problem: the evolu- tion of a 105 particle King model orbiting around a xed potential that represents our Galaxy. We used a spherical harmonics expansion N-body code, in which the potential center determination is crucial for the correct force computation. We compared this simulation with another one in which a method previously used to determine the expansion center is employed (White 1983). Our routine gives better results in energy conservation and mass loss.Fil: Aguilar, L. A.. Universidad Nacional Autonoma de Mexico. Instituto de Astronomia; MéxicoFil: Cruz, F.. Universidad Nacional Autonoma de Mexico. Instituto de Astronomia; MéxicoFil: Carpintero, Daniel Diego. Universidad Nacional de la Plata. Facultad de Ciencias Astronómicas y Geofísicas; ArgentinaUniversidad Nacional Autónoma de México. Instituto de Astronomía2002-09info:eu-repo/semantics/articleinfo:eu-repo/semantics/publishedVersionhttp://purl.org/coar/resource_type/c_6501info:ar-repo/semantics/articuloapplication/pdfapplication/pdfhttp://hdl.handle.net/11336/36956Aguilar, L. A.; Cruz, F.; Carpintero, Daniel Diego; A new method to find the potential center of N-body systems; Universidad Nacional Autónoma de México. Instituto de Astronomía; Revista Mexicana de Astronomía y Astrofísica; 38; 9-2002; 225-2310185-1101CONICET DigitalCONICETenginfo:eu-repo/semantics/altIdentifier/url/http://adsabs.harvard.edu/abs/2002RMxAA..38..225Cinfo:eu-repo/semantics/altIdentifier/url/http://www.astrosmo.unam.mx/rmaa/RMxAA..38-2/PDF/RMxAA..38-2_cruz.pdfinfo: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-10-15T14:23:19Zoai:ri.conicet.gov.ar:11336/36956instacron: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-10-15 14:23:19.722CONICET Digital (CONICET) - Consejo Nacional de Investigaciones Científicas y Técnicasfalse
dc.title.none.fl_str_mv A new method to find the potential center of N-body systems
title A new method to find the potential center of N-body systems
spellingShingle A new method to find the potential center of N-body systems
Aguilar, L. A.
Globular clusters
Numerical methods
N Body simulations
title_short A new method to find the potential center of N-body systems
title_full A new method to find the potential center of N-body systems
title_fullStr A new method to find the potential center of N-body systems
title_full_unstemmed A new method to find the potential center of N-body systems
title_sort A new method to find the potential center of N-body systems
dc.creator.none.fl_str_mv Aguilar, L. A.
Cruz, F.
Carpintero, Daniel Diego
author Aguilar, L. A.
author_facet Aguilar, L. A.
Cruz, F.
Carpintero, Daniel Diego
author_role author
author2 Cruz, F.
Carpintero, Daniel Diego
author2_role author
author
dc.subject.none.fl_str_mv Globular clusters
Numerical methods
N Body simulations
topic Globular clusters
Numerical methods
N Body simulations
dc.description.none.fl_txt_mv We present a new and fast method to nd the potential center of an N-body distribution. The method uses an iterative algorithm which exploits the fact that the gradient of the potential is null at its center: it uses a smoothing radius to avoid getting trapped in secondary minima. We have tested this method on several random realizations of King models (in which the numerical computation of this center is rather dicult, due to the constant density within their cores), and com- pared its performance and accuracy against a more straightforward, but computer intensive method, based on cartesian meshes of increasing spatial resolution. In all cases, both methods converged to the same center, within the mesh resolution, but the new method is two orders of magnitude faster. We have also tested the method with one astronomical problem: the evolu- tion of a 105 particle King model orbiting around a xed potential that represents our Galaxy. We used a spherical harmonics expansion N-body code, in which the potential center determination is crucial for the correct force computation. We compared this simulation with another one in which a method previously used to determine the expansion center is employed (White 1983). Our routine gives better results in energy conservation and mass loss.
Fil: Aguilar, L. A.. Universidad Nacional Autonoma de Mexico. Instituto de Astronomia; México
Fil: Cruz, F.. Universidad Nacional Autonoma de Mexico. Instituto de Astronomia; México
Fil: Carpintero, Daniel Diego. Universidad Nacional de la Plata. Facultad de Ciencias Astronómicas y Geofísicas; Argentina
description We present a new and fast method to nd the potential center of an N-body distribution. The method uses an iterative algorithm which exploits the fact that the gradient of the potential is null at its center: it uses a smoothing radius to avoid getting trapped in secondary minima. We have tested this method on several random realizations of King models (in which the numerical computation of this center is rather dicult, due to the constant density within their cores), and com- pared its performance and accuracy against a more straightforward, but computer intensive method, based on cartesian meshes of increasing spatial resolution. In all cases, both methods converged to the same center, within the mesh resolution, but the new method is two orders of magnitude faster. We have also tested the method with one astronomical problem: the evolu- tion of a 105 particle King model orbiting around a xed potential that represents our Galaxy. We used a spherical harmonics expansion N-body code, in which the potential center determination is crucial for the correct force computation. We compared this simulation with another one in which a method previously used to determine the expansion center is employed (White 1983). Our routine gives better results in energy conservation and mass loss.
publishDate 2002
dc.date.none.fl_str_mv 2002-09
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/36956
Aguilar, L. A.; Cruz, F.; Carpintero, Daniel Diego; A new method to find the potential center of N-body systems; Universidad Nacional Autónoma de México. Instituto de Astronomía; Revista Mexicana de Astronomía y Astrofísica; 38; 9-2002; 225-231
0185-1101
CONICET Digital
CONICET
url http://hdl.handle.net/11336/36956
identifier_str_mv Aguilar, L. A.; Cruz, F.; Carpintero, Daniel Diego; A new method to find the potential center of N-body systems; Universidad Nacional Autónoma de México. Instituto de Astronomía; Revista Mexicana de Astronomía y Astrofísica; 38; 9-2002; 225-231
0185-1101
CONICET Digital
CONICET
dc.language.none.fl_str_mv eng
language eng
dc.relation.none.fl_str_mv info:eu-repo/semantics/altIdentifier/url/http://adsabs.harvard.edu/abs/2002RMxAA..38..225C
info:eu-repo/semantics/altIdentifier/url/http://www.astrosmo.unam.mx/rmaa/RMxAA..38-2/PDF/RMxAA..38-2_cruz.pdf
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/pdf
dc.publisher.none.fl_str_mv Universidad Nacional Autónoma de México. Instituto de Astronomía
publisher.none.fl_str_mv Universidad Nacional Autónoma de México. Instituto de Astronomía
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_ 1846082641464393728
score 13.22299