A&A 456, 517-522 (2006)
DOI: 10.1051/0004-6361:20064960

Absolute motions of globular clusters

II. HST astrometry and VLT radial velocities in NGC 6397[*]

A. P. Milone1 - S. Villanova1 - L. R. Bedin2 - G. Piotto1 - G. Carraro1,3 - J. Anderson4 - I. R. King5 - S. Zaggia6


1 - Dip. di Astronomia, Univ. degli studi di Padova, vic. Osservatorio 2, 35122 Padova, Italy
2 - European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching, Germany
3 - Andes Fellow, Departamento de Astronómia, Universidad de Chile, Casilla 36-d, Santiago, Chile
4 - Dept. of Physics and Astronomy, Mail Stop 108, Rice University, 6100 Main Street, Houston, TX 77005, USA
5 - Dept. of Astronomy, Univ. of Washington, Box 351580, Seattle, WA 98195-1580, USA
6 - INAF - Osservatorio Astronomico di Trieste, via Tiepolo 11, 34131 Trieste, Italy

Received 3 February 2006 / Accepted 27 April 2006

Abstract
In this paper we present a new, accurate determination of the three components of the absolute space velocity of the Galactic globular cluster NGC 6397 ( $l\simeq338^\circ$, $b\simeq-12^\circ$). We used three HST/WFPC2 fields with multi-epoch observations to obtain astrometric measurements of objects in three different fields in this cluster. The identification of 33 background galaxies with sharp nuclei allowed us to determine an absolute reference point and measure the absolute proper motion of the cluster. The third component was obtained from radial velocities measured from spectra from the multi-fiber spectrograph FLAMES at UT2-VLT. We find ( $\mu_\alpha \cos{\delta}$, $\mu_\delta$) $_{\rm J2000.0}$ =( $+3.39 \pm 0.15$, $-17.55 \pm 0.15$) mas yr-1 and $V_{\rm rad} =
+18.36 \pm 0.09$ ($\pm$0.10) km s-1. Assuming a Galactic potential, we calculate the cluster orbit for various assumed distances and briefly discuss the implications.

Key words: astrometry - globular clusters: individual: NGC 6397 - techniques: spectroscopic - Galaxy: kinematics and dynamics

1 Introduction

After 15 years of activity, the Hubble Space Telescope (HST) is producing large volumes of accurate astrometry from the measurements that the high spatial resolution of WFPC2 and ACS images allows. In the archive are numerous fields with multi-epoch observations that are well separated in time, and first-epoch observations exist for even more fields. These are a gold mine for proper-motion measurements in a large variety of Galactic and extragalactic fields, including the crowded regions of globular clusters (GCs).

Globular clusters are important Galactic tracers; their orbital motion offers insight into the kinematics and dynamics of the halo, and their motion provides important clues to the distribution of mass in the Galaxy. Moreover, knowledge of GC orbits enables us to study their interaction with the Galactic disk and halo and is important to dynamical studies of the relationship of a cluster present-day mass function to the initial one.

In the first paper of this series (Bedin et al. 2003, Paper I), we presented the absolute proper motion of the globular cluster M 4, using a background QSO as reference. In this paper we use the background galaxies in the field of another nearby GC, NGC 6397. We made use of archival material, as well as of new images taken as a second epoch (GO-8656). We also used ground-based spectra from the multi-fiber spectrograph FLAMES (Pasquini et al. 2002) on UT2-VLT to measure the radial velocity of the cluster. Finally, with the new measurements of proper motion and radial velocity, and an assumed Galactic potential, we explore the orbit of the cluster.


  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{4960fig1.eps} \end{figure} Figure 1: Finding chart of the WFPC2 fields considered for this work. The region is $20'\times 20'$ wide and is centered on the center of NGC 6397. The fields measured in this work are labeled "B'',"C'', and "D''. For completeness we also show field "A'', which we did not use here (see text) but will use in forthcoming papers. Fields A, B, C, D are centered 5', 4 $.\mkern -4mu^\prime $5, 8', and 8' from the cluster center, respectively.
Open with DEXTER

2 Observations, data reduction

2.1 WFPC2/HST: proper motions

In Fig. 1 we use equatorial coordinates to show the footprints of the three WFPC2 fields of NGC 6397 that we used in the present work. The center of the figure coincides with the cluster center. In the following we refer to each field as B, C, and D, as labeled in the figure. Field A was not used in the present work because the two epochs were oriented differently, and this would have required a different approach to the reduction, since the chips can no longer be treated independently.

Table 1: Description of the data set used in this work.

Table 1 lists for our three fields the exposures that we used. Almost all of the image sets are dithered well, with both large and fractional pixel offsets, allowing precise astrometric measurements and a more accurate assessment of the errors.


  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{4960fig2.eps} \end{figure} Figure 2: Vector-point diagram of the proper motions in equatorial coordinates, combining all the fields. Open circles with error bars represent the selected reference galaxies. Dots represent the stars; a circle at 3 mas/yr separates cluster members from field objects. The arrow is the adopted absolute proper motion of NGC 6397.
Open with DEXTER

Positions were measured for each filter, chip, and epoch independently by using the effective-point-spread-function (ePSF) algorithms described by Anderson & King (2000). As in Paper I, for proper-motion measurements we used local transformations, and the best distortion corrections (Anderson & King 2003), after removing the 34th-row error (Anderson & King 1999). The photometry was calibrated to the Vega-mag flight system, following the recipe of Holtzman et al. (1995). For stars with >100 DN in their brightest pixel, we were able to measure positions in a single image with an error of <0.04 pixel in each coordinate.

In order to avoid introducing systematic errors, we measured the proper motions from two epochs in the same filter. When two epochs were available for each of two filters (fields B and C), the proper motions were obtained from the weighted mean of the two measurements. For field B, we linked F814W (2001) with F814W (1994) and also F606W (2001) with F606W (1996); in field C we combined F814W (2001) with F814W (1996), and also F606W (2000) with F606W (1996). In field D we matched F814W (2001) with F814W (1995).

A visual inspection of the images reveals many background galaxies. Since 33 of these have a point-like nucleus, we used our ePSF-fitting procedure to measure their positions. As was to be expected, the errors found from multiple measurements of galaxy positions are much larger than those that are typical of stars, and they are a strong function of the galaxy morphology. Nevertheless, those galaxies that give good positions provide an ideal reference point for the measurement of the absolute proper motions of stars.

Figure 2 shows the vector point diagrams for the three fields B, C, and D combined. Small points show the proper motions of the stars, and open circles show the motions measured for the selected galaxies. For each galaxy we also show the proper-motion error bar. The axes are parallel to right ascension and declination. Cluster members were used to compute the transformations and to derive the relative proper motions (see Paper I for details). For the pixel scale we adopted the number given by Anderson & King (2003), and for the orientations we used information in the HST image headers. Cluster members constitute the concentrated distribution of points at the lower left. The separation between field stars and cluster stars is very well defined. Finally, the zero point of the figure was placed at the weighted mean motion of the galaxies; this converts all the motions to absolute. The centroid proper motions of the cluster stars are in good agreement for the three fields, as is shown in Table 2.

Taking the weighted mean of the positions of the whole sample of 33 galaxies as the zero point of the motion, we find for NGC 6397 an absolute proper motion, in the J2000.0 system, of

\begin{displaymath}\mu_\alpha \cos\delta = +3.39 \pm 0.15\; {\rm mas\; yr^{-1}},
\end{displaymath}


\begin{displaymath}\mu_\delta = -17.55 \pm 0.15\; {\rm mas\; yr^{-1}}.
\end{displaymath}

In Galactic coordinates this corresponds to

\begin{displaymath}\mu_\ell \cos b = -13.52 \pm 0.15\; {\rm mas\; yr^{-1}},
\end{displaymath}


\begin{displaymath}\mu_b = -11.70 \pm 0.15\; {\rm mas\; yr^{-1}.}
\end{displaymath}

The dominant source of error is the uncertainty in the zero point given by the galaxies.

We note that the uncertainties of the results for the individual fields are quite different. This in turn rests on the wide range of the measurement errors for the individual galaxies, as indicated by the lengths of the error bars in Fig. 2. Nearly all of the weight comes from the few galaxies that have the best measurements; the majority of these are in Field B. This shows the importance of finding well-exposed galaxies with sharp centers.

These values are in marginal agreement with those found by Cudworth & Hanson (1993), which give a proper motion that is $\sim$13% smaller. However, their values were based on a correction from relative to absolute proper motion that was obtained from a Galactic model, and their procedure seems to provide an analogous discrepancy ($\sim$20% smaller) in the case of M 4 (Cudworth & Rees 1990), when compared with more recent direct determinations of its absolute proper motion (Dinescu et al. 1999; and Paper I).

Table 2: Summary of absolute motion (in [mas/yr]) for each field. The number of background galaxies used to compute the reference point is also indicated.

2.2 FLAMES+GIRAFFE/VLT: radial velocities

With the aim of measuring the dispersion of radial velocities in NGC 6397, in order to estimate the distance of the cluster by comparing with the dispersion of proper motions (see Bedin et al. 2003b for a brief description of the method), we collected spectra of 1508 stars at the UT2-VLT telescope in May 2003 with the FLAMES/GIRAFFE spectrograph (proposal No. 71.D-0076). In the present paper we use these spectra to measure the mean radial velocity of the cluster. Observations were performed in the MEDUSA mode with simultaneous calibration, which allows simultaneously obtaining about 130 spectra of objects over a field of view of 25 arcmin in diameter. The high-resolution set-up H9 (now called H9B) was used; this set-up allows a spectral resolution R=26 000 in the 5140-5350 Å range. This spectral region is well-suited for radial-velocity measurements because of the large number of metallic lines and the absence of H features. It also includes the Mg triplet, which is particularly important in radial-velocity determinations for metal-poor stars like those of NGC 6397 ( $\rm [Fe/H] \sim-2.0$), because it is the strongest feature also clearly visible in low S/N spectra.

The target stars for radial velocities were primarily selected on the basis of their position in the color magnitude diagram (CMD). The primary criterion used was to select stars along the red giant branch (RGB), the subgiant branch (SGB), and the main sequence (MS) down to V=17.5, on a narrow stripe 0.08 mag wide in color. This allowed us to drastically limit the contamination from field stars. Around the turn-off the stars of NGC 6397 have a color that is 0.2 mag bluer than that of the contaminating stars of the Galactic disk. Moreover, the evolved stars in NGC 6397 occupy a position in the CMD that is well above the location of the nearby disk stars. By analyzing in an outer field, far from the cluster center, a stripe in the CMD of a similar size to our selection strip, we estimated that on average only 12 stars should be non-members, i.e., a smaller contamination than 1%. (Furthermore, our radial velocities showed that stars in this part of the CMD are non-members.) This left virtually no contamination in the sample.

Data were automatically pre-reduced using the GIRAFFE pipeline GIRBLDRS[*], in which the spectra have been de-biased, flat-field corrected, extracted, and wavelength-calibrated, using both prior and simultaneous calibration-lamp spectra. The resulting spectra have a dispersion of 0.05 Å/pixel. A sky correction was applied to each stellar spectrum by using the IRAF sarith subroutine for subtraction of the average of the sky spectra that were observed simultaneously (same FLAMES plate) with the stars.

Radial velocities were obtained from the IRAF fxcor cross-correlation subroutine. Stellar spectra were cross-correlated with a synthetic template calculated by SPECTRUM[*], for the mean temperature, $\log{g}$, and metallicity appropriate for each observed star. The accuracy of this template then was tested on a solar spectrum obtained with the UVES spectrograph; the systematic error is less than 50 m s-1.

From multiple observations of 81 stars it was also possible to make a direct estimate of the errors. For stars of magnitude V=14, 15, 16, 17, we find errors of 0.2, 0.3, 0.7, and 1.4 km s-1, respectively. We also searched for possible systematic errors due to misalignment of the fibers or to the difference in the illumination of the simultaneous calibration fibers and the fibers used for the stars[*]. To this end we reduced the daytime wavelength-calibration plates (the plates where all fibers are pointed at the wavelength-calibration ThAr lamp) using the same procedure that we used for our stellar plates. Then we measured the velocity shift of the fibers, using as a template spectrum the same ThAr template that was used for the wavelength calibration. As expected, the average radial velocity obtained from the spectra from the five simultaneus calibration fibers is zero, within the errors. By contrast, the radial velocities from the spectra from the remaining fibers show a systematic shift, different for each fiber, probably due to their incorrect alignment at the entrance of the spectrograph. For the purposes of the present paper, it is sufficient to consider the mean systematic shift in radial velocity. We obtained for the two positioner plates Medusa 1 and Medusa 2:

\begin{eqnarray*}\Delta V_{\rm rad}~ [{\rm Medusa\ 1, NGC~6397}] &=& -0.175 \pm ...
...rm Medusa\ 2, NGC~6397}] &=& -0.225 \pm 0.008\ {\rm km~s^{-1}},
\end{eqnarray*}


and the mean value is

\begin{eqnarray*}\Delta V_{\rm rad} = -0.200 \pm 0.025\ {\rm km~s^{-1}},
\end{eqnarray*}


which must be subtracted from the cluster mean radial velocity.

As a further test we measured the systematic error affecting similar observations of another cluster, M 4, which were collected a month later. In this way we can test the time stability of the spectrograph by seeing whether these numbers remain constant. We obtained in an analogous way for the data-set of M 4 the mean systematic shift in radial velocity, for the two positioner plates Medusa 1 and Medusa 2:

\begin{eqnarray*}\Delta V_{\rm rad}~ [{\rm Medusa\ 1, M~4}] &=& -0.176 \pm 0.008...
... [{\rm Medusa\ 2, M~4}] &=& -0.236 \pm 0.008\ {\rm
km~s^{-1}}.
\end{eqnarray*}


These values agree very well with the previous ones. From these results we can assert that GIRAFFE is stable, at least for a period of about a month.

The spatial distribution of our stars, and their heliocentric radial velocities as a function of visual magnitude, are shown in Figs. 3 and 4.

  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{4960fig3.eps} \end{figure} Figure 3: Spatial distribution of stars with radial velocity measurement. The circles indicate the radial distances of 5 and 10 arc min from the center of the cluster.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{4960fig4.eps} \end{figure} Figure 4: Radial velocities of the observed sample of stars as a function of magnitude.
Open with DEXTER

The not-weighted average radial velocity[*], corrected for the systematic shift discussed previously and clipped for outliers more than 15 km s-1 from the mean, is

\begin{displaymath}V_{\rm rad} = 18.36 \pm 0.09 (\pm 0.03)\ {\rm km ~ s^{-1}}.
\end{displaymath}

The two errors that are given are first the random error and then the uncertainty in the correction for the systematic shift. This mean comes from 1486 stars (22 were rejected).

However, there are other effects of the same size that afflict spectroscopic measurements of radial velocities, such as the blue-shift of convective bubbles in giant stars or gravitational redshift in dwarf stars (Dravins et al. 1999). Since we have observed both giant and dwarf stars, we add a global systematic error of 0.10 km s-1 to our uncertainties. Finally,

\begin{displaymath}V_{\rm rad} = 18.36 \pm 0.09 (\pm 0.10)\ {\rm km ~ s^{-1}}.
\end{displaymath}

[When we analyze the cluster internal dynamics, we will carefully correct for all these effects, and also each fiber will be corrected individually for the calibration bias described above.]

Our heliocentric average radial velocity for NGC 6397 is within ${\sim}1.3\sigma$ of the Da Costa et al. (1977) value of $V_{\rm
rad}=19.8\pm1.1$ km s-1, within ${\sim}2.1\sigma$ of the Dubath et al. (1997) value of $V_{\rm rad}=15.1\pm1.6$ km s-1, but there is a large disagreement ( ${\sim}5\sigma$) with the $V_{\rm
rad}=18.80\pm0.10$ km s-1 given by Meylan & Mayor (1991).

Dubath et al. measured the radial velocity on an integrated-light echelle spectrum of the cluster. Kinematical measurements of this kind of spectra are known to be affected by strong observational biases due to small-number statistics and strong luminosity weighting (of a few red giants). As clearly demonstrated by Zaggia et al. (1992, 1993), these biases influence both velocity dispersion and radial velocity: systematic offsets can be as large as $1\times$ the velocity dispersion itself. Therefore it is not surprising to find a discrepancy between the integrated-light radial velocity and the average radial velocity of individual stars.

It is more problematic to explain the discrepancy with the Meylan & Mayor (1991) radial velocity, considering that errors in their single velocities are similar to ours, even though their sample is only ${\sim}1/12$ the size of our. Their measure is based on CORAVEL spectrometer observations of 127 RGB and SGB stars. A possible source of systematic error could possibly be in a spectral mismatch of the CORAVEL mask with the low-metallicity spectra of NGC 6397 stars, as Meylan & Mayor indeed mention in their paper. Another source of systematic error could be the presence of unrecognized field stars which could bias the mean radial velocity, considering the low velocity of NGC 6397. Meylan & Mayor observed mainly RGB stars in NGC 6397, while the majority of our targets are turn-off and MS stars, which should be relatively free of field contamination. Unfortunately a direct comparison of the radial velocities of the two datasets is not possible.

Considering the work done to limit all possible systematics in the GIRAFFE spectra, we believe that with a final systematic error <0.1 km s-1 and the large-number statistics of nearly 1600 radial velocities, our mean radial velocity is a significant improvement over previous measurements.

Table 3: Input conditions for orbit calculation. Space velocities in km s-1.


  \begin{figure}
\par\includegraphics[width=8cm,clip]{4960fig5.eps}\protect\end{figure} Figure 5: Calculated orbits for assumed distances 1.8, 2.0, 2.2, and 2.4 kpc.
Open with DEXTER

3 Calculation of the orbital parameters

Our cluster radial velocity and absolute proper motion allow us to derive the three velocity components of NGC 6397 and estimate its Galactic orbit. The new orbit should help us to trace back the dynamical history of the cluster, clarify how close it goes to the Galactic center, and assess the possible impact of the motion on its internal dynamics and its mass function.

Orbit integration for a globular cluster requires adopting a model of the potential of the Milky Way. We chose that of Allen & Santillan (1991), which was constructed to fit an assumed rotation curve and a distance and velocity of the Sun with respect to the Galactic center. The potential is time-independent, so clearly a crude approximation, because a significant variation of the Galactic potential is expected over a typical globular-cluster lifetime. Nevertheless, it is reasonable to believe that the Galactic potential has not changed much in the last few Gyr, so that the derived parameters for the orbit of NGC 6397, such as the apo- and perigalacticon, can be considered to be reasonable estimates.

The Galactic potential of Allen & Santillan (1991) assumes density distributions in the bulge, disk, and halo components. This gravitational potential is time-independent, axisymmetric, fully analytic, and mathematically very simple. It has already been used to derive the Galactic orbits of nearby stars (Allen & Martos 1986), open clusters (Carraro & Chiosi 1994), and disk and halo GCs (Odenkirchen & Brosche 1992). The initial conditions are given in Table 3 for different cluster distances, assuming $R_{\circ} = 8.5$ kpc, and the value of $V_{\circ} / R_{\circ}$ derived in Paper I.

The integration routine is a fifteenth-order symmetric, simplectic Runge-Kutta, using the Radau scheme (Everhart 1985); it guarantees energy and momentum conservation at the levels of 10-12 and 10-9, respectively, over the whole orbit integration. The orbits have been integrated back in time for 5 Gyr and are shown in Fig. 5.

The orbital parameters are summarized in Table 4, where Col. (1) lists the adopted cluster heliocentric distance, Col. (2) the z-component of the angular momentum, Col. (3) the total energy, Col. (4) the orbital period (for radial motion and for vertical motion), Cols. (5) and (6) the apo- and peri-center of the orbit, Col. (7) the maximum vertical distance the cluster reaches and Col. (8) the eccentricity, defined as $(R_{\rm a}-R_{\rm p})/(R_{\rm a}+R_{\rm p})$. The orbit parameters that we found do not differ much from those published by Dauphole et al. (1996) and by Dinescu et al. (1999), which means basically that the Galactic model we adopt in this paper is consistent with those used by the previous authors. Indeed, Dauphole et al. (1996) used an earlier version (Allen & Martos 1986) of the Galactic model we use here (Allen & Santillan 1991). In any case, our newly derived space velocities make our results an improvement over previous investigations. In particular, this is the first time that proper motions calibrated directly against extragalactic sources have been used in studying the orbit of NGC 6397.

Independent of the initial conditions chosen, the orbit of NGC 6397 has a boxy nature. The cluster oscillates rapidly through a dense part of the Galactic plane, only a few kiloparsecs from the Galactic center. This kind of orbit makes NGC 6397 very vulnerable to tidal shocks, with a consequent preferential loss of low-mass stars, which, as first recognized by Piotto et al. (1997, see their detailed discussion), is the likely cause of its anomalously flat mass function.

As the heliocentric distance increases (from bottom to top panel) the cluster tends to have a longer period, to reach greater heights above the Galactic plane and to show a larger epicyclic amplitude, dipping closer to the Galactic center. The assumed distance has a complicated effect on tidal shocks, which work over time to flatten the mass function of the cluster. Tidal shocks occur three times in each orbital period: at perigalactic passage and at each passage through the Galactic plane. For a greater assumed distance the cluster is shocked more frequently, and the shocks due to the Galactic center are stronger, but the shocks due to the plane occur at higher cluster velocity and are thus less effective.

Table 4: Orbit parameters, for the four different distances. Unit description: d[kpc], $L_{\rm z}$[kpc km s-1], $E_{\rm tot}$[ ${\rm 10 \times km^2~ s^{-2}}$], P[Myr], $R_{\rm a}$[kpc], $R_{\rm p}$[kpc], $z_{\rm max}$[kpc], e[pure number].

Note that the major source of uncertainty in the initial conditions is our ignorance of the cluster's distance. A more accurate distance determination is in progress, based on the same observing material that is presented in this paper. It will be a fundamental step in our understanding of the motion of NGC 6397.

Acknowledgements
G.P. and A.M. acknowledge the support by the MIUR under the program PRIN2003. I.R.K. and J.A. were supported by STScI grants GO-8656 and AR-8736. S.Z. was supported by the MIUR program PRIN2004. We thank the referee, Dr. C. Pryor, for a careful reading of the manuscript and for many useful comments.

References

 

Copyright ESO 2006