Softened potentials and the multipolar expansion

 

F. C. Wachlin¹ and D. D. Carpintero²

 

¹ Facultad de Ciencias Astronómicas y Geofísicas, Universidad Nacional de La Plata, and Instituto de Astrofísica de La Plata, CONICET, Paseo del Bosque S/N, 1900 La Plata, Prov. Buenos Aires, Argentina (fcw@fcaglp.unlp.edu.ar).

² Facultad de Ciencias Astronómicas y Geofísicas, Universidad Nacional de La Plata, and Instituto de Astrofísica de La Plata, CONICET, Paseo del Bosque S/N, 1900 La Plata, Prov. Buenos Aires, Argentina (ddc@fcaglp.unlp.edu.ar).

 

Received 2005 November 3
Aaccepted 2006 May 18

 

RESUMEN

Cuando se calcula el desarrollo multipolar del potencial gravitatorio, los distintos multipolos quedan bien definidos, correspondiendo cada uno a una suma finita de términos de la serie. Sin embargo, al usar el potencial gravitatorio en simulaciones numéricas, suele desarrollarse en serie una versión suavizada del mismo. Ocurre que, en estos casos, el desarrollo multipolar estándar que suele utilizarse ya no aísla los multipolos, sino que cada uno de ellos queda distribuido en infinitos términos. En este artículo se muestra cómo recuperar los multipolos completos en estos casos. Afortunadamente, la diferencia entre usar multipolos incompletos y completos es despreciable en los casos de interés, por ejemplo, en su uso en códigos árbol.

 

ABSTRACT

When the gravitational potential is developed in a multipolar series, each multipole is well defined and corresponds to a finite sum of terms in the series. In order to use the gravitational potential in numerical simulations, however, a multipolar expansion is usually applied to a softened Newtonian potential. It turns out that the commonly used multipolar expansion in this case no longer isolates each multipole as in the former case; instead, each multipole is spilled over an infinity of terms. In this paper we show how to recover the complete multipoles. Fortunately, the overall effect of using incomplete multipoles instead of complete ones turns out to be negligible in the cases of interest, for example, in its use in tree codes.

Key Words: GALAXY: KINEMATICS AND DYNAMICS - METHODS: N-BODY SIMULATIONS - METHODS: NUMERICAL.

 

DESCARGAR ARTÍCULO EN FORMATO PDF

 

REFERENCES

Aguilar, L. A., & White, S. D. M. 1985, ApJ, 295, 374        [ Links ]

Ascasibar, Y., & Binney, J. 2005, MNRAS, 356, 872        [ Links ]

Barnes, J. E., & Hut, P. 1986, Nat, 324, 446        [ Links ]

----------. 1989, ApJS, 70, 389        [ Links ]

Dehnen, W. 2001, MNRAS, 324, 273        [ Links ]

Dyer, C. C, & Ip, P. S. S. 1993, ApJ, 409, 60        [ Links ]

Gradshteyn, I. S., & Ryzhik, I. M. 1980, Table of Integrals, Series and Products (New York: Academic Press)        [ Links ]

Hernquist, L., & Barnes, J. E. 1990, ApJ, 349, 562        [ Links ]

Kellogg, O. D., 1954, Foundations of Potential Theory (New York: Dover)        [ Links ]

King, I. R. 1966, AJ, 71, 64        [ Links ]

Merritt, D., & Tremblay, B. 1994, AJ, 108(2), 514        [ Links ]

Silverman, B. W. 1998, Density Estimation for Statistics and Data Analysis (Boca Raton: Chapman & Hall/CRC)        [ Links ]

van Albada, T. S. 1982, MNRAS, 201, 939        [ Links ]

Villumsen, J. V. 1982, MNRAS, 199, 493        [ Links ]

White, S. D. M. 1983, ApJ, 274, 53        [ Links ]