Potentialdensity pairs and vertical tilt of the stellar velocity ellipsoid
O. Bienaymé
Université de Strasbourg, CNRS, Observatoire Astronomique, France
Received 16 January 2008 / Accepted 9 March 2009
Abstract
We define new potentialdensity pairs and examine the
impact of the potential flattening on the vertical velocity ellipsoid
tilt, .
By means of numerical integrations and analytical
calculations, we estimate
in a variety of galactic
potentials. We show that at 1 kpc above the Galactic plane at the
solar radius,
can differ by 5 degrees, depending on whether
the dark matter halo is flat or spherical. This result excludes the
possibility of an extremely flattened Galactic dark halo.
Key words: gravitation  Galaxy: kinematics and dynamics  Galaxy: structure  stars: kinematics
1 Introduction
The velocity ellipsoid shape can be determined from the stellar kinematics in our Galaxy, see for instance Chereul et al. (1998) from Hipparcos data, Soubiran et al. (2008) from distant red giants or Siebert et al. (2008) from Rave data. It is known that the vertical tilt of the velocity ellipsoid is related to the galactic gravitational potential and that, in the peculiar case of the Stäckel potentials, it only depends on the gravitational potential. Probably because of the lack of data, very few works have examined this question since the pioneering work of Ollongren (1962). We can note the work by Kuijken & Gilmore (1989) that proposes an approximate determination of the inclination at the solar position above the galactic plane and that points out the necessity to consider this inclination to determine the K_{z} force out of the galactic plane. There are also the analytical studies by Cuddeford & Amendt (1991) of the Jeans equations to the fourth order and their analytical approximation of the vertical tilt. Many other studies concern Stäckel potentials for which the ellipsoid orientation is easily obtained.
Thus for realistic galactic potentials and from orbit integration, both Kuijken & Gilmore (1989) and Binney & Spergel (1983) show that at 1 kpc out of the galactic plane above the Sun, the velocity ellipsoid points towards a point located at 5 to 10 kpc behind the galactic center (see also Kent & de Zeeuw 1991; and Shapiro et al. 2003).
In Sect. 2, we define a new potentialdensity pair to measure (Sect. 3) the vertical tilt of the velocity ellipsoid numerically within a spheroidal potential by varying its flattening. In Sect. 4, we generalize the potentialdensity pair to a threeparameter family with properties similar to the Miyamoto & Nagai potentials, but with flat circular velocity curves at large radii. In Sect. 5, we define a simple expression to estimate the tilt and examine its accuracy.
2 A new potentialdensity pair
We make the definition of the Mestel disk more general by utilizing
the scalefree axisymmetric potential
It generates, within the plane of symmetry z=0, a constant circular velocity curve of amplitude . The resulting density is always positive:
Figures 1, 2 show the isodensity and isopotential contours for some values of q. The axis ratio of isodensity contours is q^{3/2} and for the isopotentials it is (1+q)/2, never less than 1/2. It can be shown (Toomre 1982) that any scalefree potential at the spherical limit (obtained here with q=1) is the potential of the singular isothermal sphere, while the flat disk limit (q=0) is the Mestel flat disk. With q >1, the Eq. (1) potential is prolate and at the limit the potential and density distributions do not depend on the z coordinate (cylindrical distribution).
3 Velocity ellipsoid tilt
Since the Eq. (1) potential is scalefree we simply need to determine the velocity ellipsoid orientation versus z at some given radius R. We build the velocity ellipsoid numerically from orbit integrations ( orbits with at z=0). The velocity ellipsoid orientation for these distribution functions is shown (Fig. 3), for R=1, versus z within the meridional plane for potentials with flattenings q=0.3, 0.7. Other approaches could be tried: see for instance, the analysis by Evans et al. (1997) from the Jeans equations for scalefree potentials with a flat rotation curve by applying an approximate third integral (see also for scalefree potentials: Qian et al. 1995; Evans & de Zeeuw 1992; Qian 1992; Hunter et al. 1984; de Zeeuw et al. 1996; Evans et al. 1997).
Figure 1: Contours of isopotentials with q=0.8 and 0.2. 

Open with DEXTER 
Figure 2: Contours of isodensities with q=0.8 and 0.2. 

Open with DEXTER 
Figure 4 shows the variation in the tilt angle at R=8 and z=1 above the midplane versus the coefficient q. varies by 5 degrees from 2.2 to 7.1 degrees depending on whether the density distribution is flat or spherical. This may be compared with the recent determination by Siebert et al. (2008) of the inclination of the velocity ellipsoid at 1 kpc below the Galactic plane. From RAVE data (Zwitter et al. 2008), they find an inclination degrees towards the Galactic plane. They show it implies a moderate flattening of the dark halo by exploring the velocity ellipsoid tilt within various galactic potentials (Dehnen & Binney 1998). Within the potential of the Besançon galactic model (Robin et al. 2003) (probably the best current Galactic baryonic mass model) that includes a spherical dark halo, we also determine the velocity ellispoid tilt from numerical orbit integration and obtain degrees (Bienaymé et al. 2009). To summarize, within a spherical Galactic potential, the tilt at 1 kpc above the plane at the solar radius is about 7 degrees, and an extreme flattening of the dark matter component is excluded, since the tilt would be close to 2 or 3 degrees (Fig. 4).
4 Other potentialdensity pairs
The Mestel disk potential can be generalized by introducing a core radius
a. It gives (Evans & Collett 1993; and quoted references
Rybicki 1974; and Zang 1976)
with the resulting surface density
and the rotation curve , flat at large R, is given by
We generalize these Mestel disks with a core by applying the Miyamoto & Nagai (1975) method to generalize the Plummer spheroid and Kuzmin disk (see also Baes 2009). Thus we define the potential
It can be verified that this potential generates a positive spheroidal density distribution. With b=0 it may generate both a flat disk and a spheroid. The sum (a+b) constrains the size of the core, and bthe thickness of the disky part of the spheroid.
Figure 3: Vertical tilt of the velocity ellipsoid in degrees versus z at R=1 for models with q=0.3 (left) and q=0.7(right). Black line: analytical estimate (see Sect. 5), red line: result of numerical orbit integration, dotted line: tilt within a spherical potential. 

Open with DEXTER 
Figure 4: Velocity ellipsoid inclination in degrees versus potential flattening q at R=8 and z=1. Dotted line: the expected inclination from Sect. 5, crosses: numerical determinations from orbit integration. 

Open with DEXTER 
Figure 5: Isolevels for the density distribution resulting from the potential given by Eq. (4) with (a,b,q)=(5, 2, 1.3), (5, 2, 0.7), (0, 5, 0.5) from left to right. 

Open with DEXTER 
Potentials given by Eqs. (1) and (3) are peculiar cases (q=0 or a=b=0) of the more general potential
and the corresponding density is (if )
with
(6) 
If a=0, it may be shortly written as
(7) 
When b=0 and , it is the combination of a spheroid and a flat disk
(8) 
(9) 
With q=1, this is the potentialdensity pair defined by Brada & Milgrom (1995).
More generally, we obtain a spheroid potential producing a rising velocity curve, flat at large radii, the core radius depending only on q and on a+b. The parameter q modifies the flattening at large R and b the thickness of the disky part of the density distribution.
Thus, this generalization introduces a new family of potentialdensity pairs with general properties similar to the Miyamoto & Nagai potentials, but with flat circular velocity curves at large radii. Figures 5 illustrate isodensity contours for three such potentials. To conclude, we recall another threeparameter family of potentialdensity pairs obtained by Zhao (1996).
5 Anaytical estimate of the vertical velocity ellipsoid tilt
We examine the reliability of the Cuddeford & Amendt formulae (1991, Eqs. (90), (91) and (E10)) and Amendt & Cuddeford (1991, Eq. (104)) that predict the vertical tilt of the velocity ellipsoid within an axisymmetric potential, and we propose another formula that is more accurate at greater distances from the plane of symmetry z=0. For various potentials, we compare the expected tilt from these two expressions with the tilt obtained from numerical orbit integrations.
Cuddeford & Amendt (1991) analyze the consecutive moments of the Boltzmann equation up to the 4th order by expanding these moments in Taylor series assuming that is small. Combining these moment equations, they derive expressions for the velocity moments and obtain an approximate expression for the tilt angle, , which only depends on the gravitational potential. They show that their expression is exact in the case of axisymmetric Stäckel potentials, peculiar potentials for which the tilt only depends on the potentials but not on the exact distribution function. According to their initial hypothesis, the formula must be generally valid when the considered distribution of stellar orbits covers a sufficiently small domain (expecting that the covered domain can be approximated with a Stäckel potential). This should be the case of nearly circular orbits with small vertical or radial oscillations.
The simplicity of the result obtained by Cuddeford & Amendt (1991) makes it extremely attractive. For instance Vallenari et al. (2006) use it to improve their recent galactic model of star counts and kinematics, which is dynamically selfconsistent locally.
Inspired by their work, we propose a more direct expression derivated
from the generic equation (Ollongren 1962, Eq. (7.2)) that
defines axisymmetric Stäckel potentials:
with a positive sign for prolate spheroidal coordinates, and a negative one for oblate ones (see de Zeeuw 1985, for a detailed description).
For a given potential ,
we substitute z_{0} from Eq. (10) in
the expression for the velocity ellipsoid tilt
given by Hori
& Liu (1963), which is exact in the case of Stäckel
potentials:
Combining Eqs. (10) and (11) gives an estimate of the velocity ellipsoid inclination within axisymmetric potentials. We note that the Amendt & Cuddeford (1991) formula could have been obtained in the same way, just by differentiating the Ollongren equation with respect to z and replacing it within Eq. (11).
A noticeable difference between their expression and ours is that their formula is defined at z=0 and is exact for the Stäckel potential only at z=0, while our expression remains exact (for Stäckel potentials) at any z. This may be why our formula is in better agreement with the results of the numerical explorations described below.
Now, to verify the reliability of the estimated tilt angle given by Eqs. (10), (11), and the one given by Amendt & Cuddeford (1991) formulae, we measure the tilt numerically within a series of potentials. For this purpose, we build stationary distribution functions with a library of orbits for each model using a 7thorder RungeKutta (Fehlberg 1968) with a relative accuracy . Each orbit is integrated for 80 rotations, the initial conditions being drawn from a Shu distribution function (we also tried a Dehnen disk distribution function) with and (we also consider the 0.20.1 and 0.60.3 pairs). One point is selected per orbit in the last 40 rotations.
Figures 6 show the inclination versus the distance z from the plane of symmetry in the case of the Brada & Milgrom (1995) disk+halo potential with a core radius a=3. The tilt is plotted versus z at four different galactic radii R=1.5, 4.5, 7.5, and 10.5. The dark line is the analytical estimate from Eqs. (10), (11), while the red line is the result of numerical orbit integrations. The agreement between numerical measures and the analytical estimate is satisfying in the range of considered positions R and z. The estimated vertical tilt is also shown in Fig. 4: at position R=8 and z=1, is plotted versus the flattening q of the Eq. (1) potential. The agreement is satisfying for q between 0 and 1.
Figure 6: Velocity ellipsoid tilt angles in degrees versus z within the Brada & Milgrom (1995) Disk+Halo potential at four galactic radii R=1.5, 4.5, 7.5 and 10.5 (from top left to bottom right). Red irregular line: numerical determination. Dark continuous line: analytical estimate. Dotted line: tilt within a spherical potential. 

Open with DEXTER 
We explored other potentials: a Kuzmin disk (a Stäckel potential to test our numerical process), the Mestel disk, the Brada & Milgrom potential, a logarithmic flattened potential, combined Kuzmin+logarithmic halo potentials. For these potentials our analytical estimate give excellent results in agreement with numerical determinations at a precision better than 1 degree at z/R=0.15.
With the slightly flattened spherical logarithmic potentials, the Amendt & Cuddeford (1991) formula is also an excellent estimate for the tilt (1 percent difference at z/R=0.15 with a potential axis ratio q=0.9). However, with the Mestel potential, it fails at any z (a null tilt is predicted), while it works correctly with the Brada and Milgrom potential, but only at radius close to or smaller than the core radius.
In conclusion, our analytical estimate of the vertical tilt of the stellar velocity ellipsoid is accurate for various potentials and vertical distances less than 0.15R. It also shows that at 1 kpc above the Galactic plane at the solar radius, varies by 5 degrees depending on whether the potential is spherical or the dark matter component is flat.
However, for an exponential disk, this estimate fails in the range R<2 or >4, in units of the scale length. The situation is similar for the potential given by Eq. (1) for ``flattening'' of q=2and neighboring values. An explanation for this deficiency is that the denominator of the r.h.s. of Eq. (11) is zero. In that case, the resulting value for z_{0} varies strongly with R or z. Since Eq. (10) is not a Stäckel fit, its domain of validity is determined by the condition that the estimate for z_{0} from Eq. (11) does not vary strongly over the considered (R,z)domain. For deviations from this limit, more sophisticated methods of fitting potentials with Stäckel potentials are required (see for instance de Bruyne et al. 2000) to obtain more reliable estimates of the vertical tilt angle.
References
 Amendt, P., & Cuddeford, P. 1991, ApJ, 368, 79 [NASA ADS] [CrossRef] (In the text)
 Baes, M. 2009, MNRAS, 392, 1503 [NASA ADS] [CrossRef] (In the text)
 Bienaymé, O., Famaey, B., Wu, X., Zhao, H. S., & Aubert, D. 2009, A&A, 500, 801 [CrossRef] [EDP Sciences] (In the text)
 Binney, J., & Spergel, D. 1983, The nearby stars and the stellar luminosity, ed. A.G.D. Philip, & U. Schenectady (N.Y.: Davis Press), IAU Col., 76, 259 (In the text)
 Brada, R., & Milgrom, M. 1995, MNRAS, 276, 453 [NASA ADS] (In the text)
 Cuddeford, P., & Amendt, P. 1991, MNRAS, 253, 427 [NASA ADS] (In the text)
 Chereul, E., Crézé, M., & Bienaymé, O. 1998, A&A, 340, 384 [NASA ADS] (In the text)
 de Bruyne, V., Leeuwin, F., & Dejonghe, H. 2000, MNRAS, 311, 297 [NASA ADS] [CrossRef] (In the text)
 Dehnen, W., & Binney, J. 1998, MNRAS, 294, 429 [NASA ADS] [CrossRef] (In the text)
 de Zeeuw, P. T. 1985, MNRAS, 216, 273 [NASA ADS] (In the text)
 de Zeeuw, P. T., Evans, N. W., & Schwarzschild, M. 1996, MNRAS, 280, 903 [NASA ADS] (In the text)
 Evans, N. W., & Collett, J. L. 1993, MNRAS, 264 353 (In the text)
 Evans, N. W., & de Zeeuw, P. T. 1992, MNRAS, 257, 152 [NASA ADS] (In the text)
 Evans, N. W., van Häfner, R. M., & de Zeeuw, P. T. 1997, MNRAS, 286, 315 [NASA ADS] (In the text)
 Fehlberg E. 1968, NASA Technical Report, R287 (In the text)
 Hori, L., & Liu, T. 1963, PASJ, 15, 100 [NASA ADS] (In the text)
 Hunter, J. H., Ball, R., & Gottesman, S. T. 1984, MNRAS, 208, 1 [NASA ADS] (In the text)
 Kent, S., & de Zeew, T. 1991, AJ, 102, 1994 [NASA ADS] [CrossRef] (In the text)
 Kuijken, K., & Gilmore, G. 1989, MNRAS, 239, 571 [NASA ADS] (In the text)
 Miyamoto, M., & Nagai, R. 1975, PASJ, 27, 533 [NASA ADS] (In the text)
 Ollongren, A. 1962, B.A.N., 16, 241 [NASA ADS] (In the text)
 Qian, E. 1992, MNRAS, 257, 581 [NASA ADS] (In the text)
 Qian, E., de Zeeuw, P. T., van der Marel, R. P., & Hunter, C. 1995, MNRAS, 274, 602 [NASA ADS] (In the text)
 Robin, A. C., Reylé, C., Derrière, S., & Picaud, S. 2003, A&A, 409, 523 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Shapiro, K., Gerssen, J., & van der Marel, R. 2003, ApJ, 126, 2707 [NASA ADS] (In the text)
 Siebert, A., Bienaymé, O., Binney, J., et al. 2008, MNRAS, 391, 793 [NASA ADS] [CrossRef] (In the text)
 Soubiran, C., Bienaymé, O., Mishenina, T. V., & Kovtyukh, V. V. 2008, A&A, 480, 91 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Toomre, A. 1982, ApJ, 259, 535 [NASA ADS] [CrossRef] (In the text)
 Vallenari, A., Pasetto, S., Bertelli, G., et al. 2006, A&A, 451, 125 [NASA ADS] [CrossRef] [EDP Sciences]
 Zang, T. A. 1976, Ph.D. Thesis, Massachusetts Institute of Technology (In the text)
 Zhao, H. S. 1996, MNRAS, 278, 488 [NASA ADS] (In the text)
 Zwitter, T., Siebert, A., Munari, U., et al. 2008, AJ, 136, 421 [NASA ADS] [CrossRef] (In the text)
All Figures
Figure 1: Contours of isopotentials with q=0.8 and 0.2. 

Open with DEXTER  
In the text 
Figure 2: Contours of isodensities with q=0.8 and 0.2. 

Open with DEXTER  
In the text 
Figure 3: Vertical tilt of the velocity ellipsoid in degrees versus z at R=1 for models with q=0.3 (left) and q=0.7(right). Black line: analytical estimate (see Sect. 5), red line: result of numerical orbit integration, dotted line: tilt within a spherical potential. 

Open with DEXTER  
In the text 
Figure 4: Velocity ellipsoid inclination in degrees versus potential flattening q at R=8 and z=1. Dotted line: the expected inclination from Sect. 5, crosses: numerical determinations from orbit integration. 

Open with DEXTER  
In the text 
Figure 5: Isolevels for the density distribution resulting from the potential given by Eq. (4) with (a,b,q)=(5, 2, 1.3), (5, 2, 0.7), (0, 5, 0.5) from left to right. 

Open with DEXTER  
In the text 
Figure 6: Velocity ellipsoid tilt angles in degrees versus z within the Brada & Milgrom (1995) Disk+Halo potential at four galactic radii R=1.5, 4.5, 7.5 and 10.5 (from top left to bottom right). Red irregular line: numerical determination. Dark continuous line: analytical estimate. Dotted line: tilt within a spherical potential. 

Open with DEXTER  
In the text 
Copyright ESO 2009