Publication:
An adaptive algorithm for N-body field expansions

Simple item page
dc.contributor.authorWeinberg, MD
dc.date2023-09-22T21:50:47.000
dc.date.accessioned2024-04-26T08:29:37Z
dc.date.available2024-04-26T08:29:37Z
dc.date.issued1999
dc.descriptionThis is a pre-published version which is collected from arXiv link. The published version is at http://iopscience.iop.org/1538-3881/117/1/629
dc.description.abstractAn expansion of a density field or particle distribution in basis functions that solve the Poisson equation both provides an easily parallelized N-body force algorithm and simplifies perturbation theories. The expansion converges quickly and provides the highest computational advantage if the lowest order potential-density pair in the basis looks like the unperturbed galaxy or stellar system. Unfortunately, there are only a handful of such bases in the literature that limit this advantage. This paper presents an algorithm for deriving these bases to match a wide variety of galaxy models. The method is based on efficient numerical solution of the Sturm-Liouville equation and can be used for any geometry with a separable Laplacian. Two cases are described in detail. First, for the spherical case, the lowest order basis function pair may be chosen to be exactly that of the underlying model. The profile may be (1) cuspy or have a core and (2) truncated or of infinite extent. Second, the method yields a three-dimensional cylindrical basis appropriate for studying galactic disks. In this case, the vertical and radial bases are coupled; the lowest order radial part of the basis function can be chosen to match the underlying profile only in the disk plane. Practically, this basis is still a very good match to the overall disk profile and converges in a small number of terms. The ease of combining several bases makes this force solver ideally suited to multicomponent simulations, such as those of disks embedded in halos.
dc.description.pages629-637
dc.identifier.doihttps://doi.org/10.1086/300669
dc.identifier.urihttps://hdl.handle.net/20.500.14394/3354
dc.relation.ispartofASTRONOMICAL JOURNAL
dc.relation.urlhttps://scholarworks.umass.edu/cgi/viewcontent.cgi?article=1069&context=astro_faculty_pubs&unstamped=1
dc.source.issue1
dc.source.issue117
dc.source.statuspublished
dc.subjectcelestial mechanics
dc.subjectstellar dynamics
dc.subjectgalaxies : structure
dc.subjectGalaxy : structure
dc.subjectmethods : numerical
dc.subjectAstrophysics and Astronomy
dc.titleAn adaptive algorithm for N-body field expansions
dc.typearticle
dc.typearticle
digcom.contributor.authorisAuthorOfPublication|email:weinberg@astro.umass.edu|Weinberg, MD
digcom.identifierastro_faculty_pubs/70
digcom.identifier.contextkey1604224
digcom.identifier.submissionpathastro_faculty_pubs/70
dspace.entity.typePublication
relation.isAuthorOfPublication6387ef0a-b08d-4ee2-b19d-775a09801631
relation.isAuthorOfPublication.latestForDiscovery6387ef0a-b08d-4ee2-b19d-775a09801631
Files
Original bundle
Now showing 1 - 1 of 1
Thumbnail Image
Name:
x_An_adaptive_algorithm_for_n_body_field_expansions.pdf
Size:
359.71 KB
Format:
Adobe Portable Document Format
Thumbnail Image
Name: x_An_adaptive_algorithm_for_n_body_field_expansions.pdf
Size: 359.71 KB
Kind: Adobe PDF
Collections
Astronomy Department Faculty Publication Series