author: L. Pasquini; C.Melo; C. Chavero ; D. Dravins; H.-G. Ludwig; P.Bonifacio; R. De La Reza
date: Received / Accepted
title: Gravitational redshifts in main-sequence and giant stars [^1]

Precise analyses of stellar radial velocities reveal intrinsic effects causing wavelength shifts of spectral lines (other than Doppler shifts due to radial motion), such as gravitational redshifts and convective blueshifts. Gravitational redshifts in solar-type main-sequence stars are expected to be some 500 m s$^{-1}$ greater than those in giants. Such a signature is searched for between groups of open-cluster stars which share the same average space motion and thus have the same average Doppler shift. 144 main-sequence stars and cool giants were observed in the M67 open cluster using the ESO FEROS spectrograph, obtaining radial velocities by cross correlation with a spectral template. Binaries and doubtful members were removed, and averaged low-noise spectra calculated for different classes of stars. M67 dwarf and giant radial-velocity distributions are well represented by Gaussian functions, sharing the same apparent average radial velocity within $\simeq$ 100 m s$^{-1}$. In addition, dwarfs in M67 appear to be dynamically hotter ($\sigma$ = 0.90 km s$^{-1}$) than giants ($\sigma$ = 0.68 km s$^{-1}$).

Explanations for the lack of an expected signal are sought: a likely cause is the differential wavelength shifts produced by different hydrodynamics in dwarf and giant atmospheres. Radial-velocity differences measured between unblended lines in low-noise averaged spectra vary with line-strenth: stronger lines are more blushifted in dwarfs than in giants, apparently 'compensating' for the gravitational redshift.

Synthetic high-resolution spectra are computed from 3-dimensional hydrodynamic model atmospheres for both giants and dwarfs, and synthetic wavelength shifts obtained. In agreement with observations, 3D models predict substantially smaller wavelength-shift differences than expected from gravitational redshift only. The procedures developed could be used to test 3D models for different classes of stars, but will ultimately require high-fidelity spectra for measurements of wavelength shifts in individual spectral lines.

# Introduction

## Gravitational redshift

Spectroscopic measurements of stellar radial velocities have reached levels of precision which not only enable new types of studies, but also require a better understanding of effects other than a star's barycentric motion that may cause displacements of spectral-line wavelengths. One such effect is the gravitational redshift originating from the propagation of light between different gravitational potentials at the source and at the observer.

Any theory involving conservation of energy requires that radiation leaving a gravitational field loses energy, and with the early 20$^\mathrm{th}$ century understanding of the photon and the relation between its energy and wavelength, predictions could be made of the wavelength displacement expected for radiation leaving various stars, and these effects have been sought observationally.

While a determination of an absolute wavelength shift requires the true stellar center-of-mass motion (and its induced Doppler shift) to be known, *relative shifts* can be determined between stars sharing a common space motion. This applies to the orbit-averaged velocity of binary stars (where each component shares the same systemic velocity), and for stars in \[open\] clusters which share the same space velocity vector.

The first claimed confirmation of the predicted gravitational redshift came from the measurement of the apparent radial velocity of Sirius B, the white-dwarf companion to Sirius (Adams 1925), identified as the difference from the motion expected in its 50-year orbit. However, this faint white dwarf is overwhelmed by scattered light from the much brighter primary and the validity (and even the honesty!) of these early measurements has been questioned (Hetherington 1980; Greenstein et al. 1985).

Gravitational redshifts have since been more reliably determined for numerous white dwarfs in binary systems and in nearby open clusters such as the Hyades and the Pleiades, finding values typically on order 30–40 km s$^{-1}$, as measured against wavelength positions in the spectra of ordinary stars. Also, very much larger gravitational redshifts have been identified from spectra of X-ray bursts on neutron-star surfaces.

The gravitational redshift is proportional to the gravitational potential at the stellar surface and scales proportional to the stellar mass, and inversely as its radius $R$:

$v_{\rm grav} = GM/Rc$,

where $G$ is the gravitational constant. The expected solar gravitational redshift for light escaping from the solar photospheric surface to infinity is 636.486 $\pm$  0.024 m s$^{-1}$ (Lindegren & Dravins 2003) and – since the shift diminishes with distance from the stellar center as $r^{-1}$ – becomes 633.5 m s$^{-1}$ for light intercepted at the Earth's mean distance from the Sun ($r = 215R_\odot$). A solar spectral line instead formed at chromospheric heights (30 Mm, say; $r = 1.04R_\odot$) will have this redshift decreased by some 20 m s$^{-1}$, and a coronal line by perhaps 100 m s$^{-1}$.

For other stars, the shift scales as $(M/M_\odot)(R/R_\odot)^{-1}$, or as $(g/g_\odot)(R/R_\odot)$, where $g$ is the surface gravity. For distant stars, additional lineshifts may originate from large-scale gravitational fields in the Milky Way Galaxy. For instance, spectra of stars in the central Galactic bulge ($R_\ast \simeq 1$ kpc) may be gravitationally redshifted by 300–400 m s$^{-1}$, while stars in the Magellanic Clouds ($R_\ast \simeq 55$ kpc) might be blueshifted by a similar amount, as seen by an observer near the solar galactocentric position at $R_{\rm obs}\simeq 8.5$ kpc (Lindegren & Dravins 2003).

Across the Hertzsprung-Russell diagram, the gravitational redshifts are not expected to vary much along the main sequence between A5 V and K0 V (being around 650 m s$^{-1}$), but to reach twice that value for more massive early-B stars with $M\simeq$ 10$M_\odot$. For giants the shifts decrease substantially and for red supergiants it should be only some 20 or 30 m s$^{-1}$, three orders of magnitude smaller than for white dwarfs (Griffin 1982; Dravins 2005). Thus, the difference between the shifts in dwarf and giant stars is expected to be some 0.5 km s$^{-1}$, orders of magnitude greater than current measurement precision. Such values are significant, both to study the precise dynamics in binary- and multiple-star systems (Griffin 1982, Pourbaix et al. 2002), or to disentangle wavelength shifts caused by stellar atmospheric dynamics by convective motions (e.g., Dravins 1982; Asplund 2005; Dravins et al. 2005). Together with photospheric line asymmetries, such intrinsic wavelength shifts make up significant diagnostic tools for probing the 3-dimensional structure of stellar atmospheres (Dravins 2008; Nagendra et al. 2009), a field where additional diagnostics are needed to reliably segregate among hydrodynamic models predicting significantly different chemical compositions, oscillation properties, and so on.

## Searching for shifts in ordinary stars

At first sight, the Sun would appear to be a promising source for determining absolute gravitational redshifts, since the radial Sun-Earth motion is both very small and well known. However, already long ago, spectral lines in the Sun were found to be shifted by different but comparable amounts, including variation of absorption-line wavelengths across the solar disk from its center towards the limb that was initially a mystery. This is now well understood as convective lineshifts, originating from a summation of contributions from hot and rising (thus locally blueshifted) elements, and cool and sinking ones, and various effects of observational perspective across the disk. Many efforts were made to measure the solar gravitational redshift (and several publications claimed its detection) although, in hindsight, almost all such studies were constrained by the limited understanding of the interplay with other mechanisms causing subtle wavelength shifts (besides effects from atmospheric motions, also the limited accuracy for laboratory wavelengths, and even effects of pressure shifts come to play). Possibly, the most convincing demonstration of a gravitational redshift in the solar spectrum was by Beckers (1977). The magnetic fields in the cores of sunspots inhibit convective motions and greatly reduce the amount of both convective wavelength shifts, and their center-to-limb variation, enabling rather precise wavelength shift measurements in magnetically insensitive lines. From such lines, Beckers deduced a shift of 613 m s$^{-1}$, consistent with the expected magnitude.

Just as for white dwarfs, systematic wavelength displacements between ordinary stars may be identified in binaries sharing a common systemic velocity, and for stars in open clusters, moving in space with sufficiently small velocity dispersion. The earliest such attempts seem to have been by E. Finlay-Freundlich who, starting around 1914, tried to find statistical evidence of gravitational redshifts (for a lucid review, see Hentschel 1994). Somewhat later, Trumpler (1935), recognized that massive O stars should be excellent targets, because they have the largest gravitational redshift while they often appeared in young clusters where they, being massive, by equipartition of energy should have the smallest velocity dispersion. In some clusters, shifts of some 10 km s$^{-1}$ were found relative to comparison stars, suggesting enormous stellar masses. A later understanding of these stars showed that shifts should be much less (Conti et al. 1977), and while effects of stellar variability, line asymmetries, and differences of apparent velocities in different spectral features did begin to get appreciated, also a certain literature grew on suggested 'anomalous' redshifts in binary high-luminosity stars, and other objects.

Only after measuring precisions for individual stars had reached levels of around 100 m s$^{-1}$, did it become practical to search for systematic differences in apparent velocities among ordinary stars, e.g., between main-sequence dwarfs and giants. One hint came from a study of the open cluster NGC 3680 by Nordström et al. (1997), who found the giants (though the sample contained only six such stars) to be blueshifted on average 0.4 km s$^{-1}$ relative to the dwarfs, as expected from gravity.

For certain groups of stars, including those in undispersed moving clusters, *absolute* stellar radial motion can be determined from precise astrometry, not invoking any spectroscopy. The difference to spectroscopic lineshifts then reveals shifts intrinsic to stellar atmospheres. In the Hyades (Madsen et al. 2002), such signatures were seen for some classes of hotter main-sequence stars, but the spectral displacements of the (only three) Hyades K-giants did not show any significant deviation from that of main-sequence K-type stars.

Different gravitational redshifts in dwarf and giant stars have thus been sought for almost a century. The value and usefulness of such measurements has been pointed out from time to time (e.g., von Hippel 1996) but conclusive or statistically significant observations are still lacking, which was a motivation for the present project.

# Radial velocities in open clusters

The aim is to observe apparent radial velocities of main-sequence and giant stars in some stellar cluster(s) with sufficient precision to identify expected differences in gravitational redshift between these groups. The ideal cluster should obviously be populous and have stars covering a large range of $M/R$ ratios.

Not all clusters are suitable. To assure good measuring precision for spectroscopic radial velocity, the stars must have spectral lines that are both narrow and numerous, which excludes young clusters where main-sequence stars still are rapidly rotating and produce broadened lines. Also, higher-temperature giants are unfit due to their few spectral lines and often complex atmospheric dynamics. Thus, the cluster must be sufficiently old, with main-sequence rotation slower than some 5–10 km s$^{-1}$, and must have a well developed red-giant branch with many cool giants that are good targets for radial-velocity measurements (e.g., Setiawan et al. 2004).

The relatively old cluster M67 (NGC 2682) at $\sim$`<!-- -->`{=html}850 pc distance was chosen as a nearly ideal object. Besides having a well-developed red-giant branch, advantages include that its composition is quite similar to solar (Tautvaišienė et al. 2000; Randich et al. 2006; Pace et al. 2008) and the age somewhat comparable: 2.6 Gy (WEBDA 2009). It has been well studied in various aspects (Pasquini et al. 2008; Sandquist 2004), including a census of binaries (Latham et al. 1992). Finally, we had already observed a homogeneous group of stars in M67 (Melo et al. 2001; hereafter Paper I), verifying the modest rotational velocities, and could use these data as a starting point.

From previous radial-velocity data (e.g., Mathieu et al. 1986), the internal velocity dispersion of the space motions within M67 was found to be $\approx$ 0.5–0.8 km s$^{-1}$, similar to that deduced from proper motions (Girard et al. 1989). Other studies suggest the expansion and rotation of the cluster core to be below 1 km s$^{-1}$ (Zhao et al. 1996; Loktin 2005).

Since these numbers are comparable to the expected difference in gravitational redshifts between dwarfs and giants, averages over many stars will be required, and attention must be paid to calibrating radial-velocity measures between stars of different spectral types, identifying stellar multiplicity, possibly different velocities for stars mass-segregated into different parts of the cluster, and the physics of line formation in dynamic stellar atmospheres. In any case, since the internal measuring precision of modern radial-velocity instruments is very much better than this dispersion, significant progress should be possible.

# Observations and data reduction

An observing program was carried out with FEROS (*Fiberfed Extended Range Optical Spectrograph*; Kaufer et al. 1999) in three different runs. The first took place during FEROS commissioning in October and November 1998, observing 28 main-sequence, turnoff, and giant stars. The radial velocities and $V\sin i$ were published in Paper I and are listed in the online Appendix for sake of completeness. Additionally 68, mostly main-sequence and turnoff stars, were observed over four nights during the first week of February 2002 (Run II) in parallel with other projects during Brazilian time. Finally in Run III, 66 additional M67 members were observed in late February 2008. For the first two runs, the spectrograph was attached to the ESO 1.5-m telescope, and during the third to the MPG/ESO 2.2-m telescope. In total 162 observations for 144 different stars were obtained.

FEROS is a bench-mounted, thermally controlled, prism-cross-dispersed échelle spectrograph. It offers rather high resolution (R $\simeq$ 48,000), and in a single recording with 39 échelle orders provides almost complete spectral coverage from $\approx$ 350–920 nm. On bright stars, FEROS has been demonstrated to be precise to better than 20 m s$^{-1}$ (Setiawan et al. 2004). Our targets are substantially fainter, but an exceedingly high precision is not required. The M67 velocity dispersion is on the order of 1 km s$^{-1}$, so a precision of a few hundred meters in a single exposure is sufficient. Exposure times were set to have a S/N around 30 or better. With this photon noise, the radial velocity (RV) precision limit is of the order 20 m s$^{-1}$. The total error budget is discussed in Section 6.1.

We aimed to observe many stars with different $M/R$ ratios, to be able to precisely determine mean velocities for different groups of stars. Known binaries were excluded (M67 is rich in those), although some single-lined spectroscopic binaries likely are still present. We also observed some stars with deviant velocities which appear to be non-members. The list for the whole sample, including data from Paper I, is given in the Appendix, available online.

# Radial velocities and membership

Radial velocities were computed as the Gaussian fit to the cross-correlation function (CCF) between the observed stellar spectrum and a spectral template mask. This binary mask was based on a stellar spectrum, following the recipe in Baranne et al. (1979). As in Paper I, a K0 III synthetic mask (Baranne et al. 1996) was used for all stars, assuring the same reference for all spectra.

The mask is composed mostly of metallic lines, whose number is sufficiently high (about 2000 lines in the 500–700 nm wavelength interval used) that we expect that line blends, which vary according to the stellar effective temperature and gravity, are substantially symmetric and do not change the average shift of the spectrum. The cross-correlation mask is built by selecting the signal from spectral lines below a given residual flux. This cut is made in order to compromise between the number of points producing the cross-correlation signal and the width of the lines: too high a threshold (including points close to the continuum) might produce too wide CCFs, decreasing the measuring precision (see Baranne et al. 1996). For our mask a value of 0.7 was selected, i.e., only wavelength points with a residual intensity of less than 0.7 in the K0 III reference spectrum contribute to the digital mask. This implies that the selected lines have relatively large equivalent widths in the M67 giants.

Only a few stars have repeated observations, and the comparison between multiple measurements is shown in Figure . The average radial velocity of these stars in the first observing runs, made at the ESO 1.5-m, is 33.93 km s$^{-1}$, which compares very well with 33.96 km s$^{-1}$ for the same objects observed at MPG/ESO 2.2-m, although the dispersion is quite large, around 0.45 km s$^{-1}$. This number is, however, an upper limit to the measurement error, because among the 14 stars observed at the two telescopes, some probably are binaries. Three stars (S1016, S1075, S1305) differ by more than 0.8 km s$^{-1}$, and once these are removed from the sample, the dispersion drops to 0.3 km s$^{-1}$. We consider this as a good estimate of the real uncertainty of a single measurement, though it might still be an upper limit. This value is similar to the dispersion in measurements of our RV standard star, HD 10700 (see section 5.1 and 6.1.2).

Figure  shows the radial velocity distribution of all the observed stars. The M67 members are clearly visible, and their RV distribution can be well represented by a Gaussian distribution centered on RV=33.73, with $\sigma$=0.83 km s$^{-1}$. We made a 2.5 $\sigma$ clipping and considered only the stars left as bona fide M67 single members. Using this selection, a total of 33 stars were discarded. We note that S1011, 1064, 1221, 1250, 1314, all classified as binaries in Paper I, appear in the list of binary stars. Out of the Paper I and Latham et al. (1992) lists, only S986 was not recovered as a binary. The final sample consists of 110 stars and their color-magnitude diagram is in Figure ; bona fide single members are marked as filled squares, stars discarded as empty squares. The targets are well distributed in the color-magnitude diagram, and numeorus evolved stars were observed to ensure good statistics also in that less populated group.

M67 is a well studied cluster, and among previous radial velocity studies, the most appropriate comparison is probably with the pioneering work by Mathieu et al. (1986), which is similar in size to ours (104 likely single members), with comparable (but slightly worse) RV precision. A difference is that their sample was limited to brighter magnitudes, since it included all proper motion candidates brighter than V=12.8, which corresponds to the region around the turnoff. Those authors used (as we do) a cross correlation technique, although their spectra were limited to a smaller spectral range, centered on 520 nm. We analyzed the Mathieu et al. (1986) radial velocities exactly in the same way as ours, finding a mean velocity 33.57 m s$^{-1}$($\sigma$= 0.78 km s$^{-1}$) and, anticipating the results of the next section, also their data do not show evidence of gravitational redshift.

# Analysis

After cleaning the sample, we performed a gross split between 'main sequence' and 'giants', defining as giants all stars along the subgiant and red-giant branches redder than $B-V$=0.7, and all the other stars as belonging to the main sequence. The giants' radii are so large that gravitational redshift is expected to be negligible in these, and therefore their average RV should be bluer than the main sequence stars by several hundred m s$^{-1}$. This simple split of the sample however does not produce the expected result, as evident from Figure :

1.  Contrary to what is expected from gravitational redshift predictions, the two distributions are centered on almost the same velocity (33.75 km s$^{-1}$)

2.  Main sequence stars and giants have substantially different velocity dispersions: $\sim$ 0.90 and 0.68 km s$^{-1}$ respectively.

In proceeding with a more complete analysis, masses and radii were computed for all stars. While the mass range is small, the radii vary considerably. These were computed using stellar bolometric magnitude and effective temperature, calculated from the $B-V$ vs. T$_\mathrm{eff}$ relation and the bolometric magnitudes of Alonso et al. (1996, 1999). A reddening of 0.041 (An et al. 2007; Taylor 2007) and a distance modulus of 9.61 (Pasquini et al. 2008) for the cluster were assumed. For the masses, a fit between mass and luminosity was derived from evolutionary isochrones (Girardi et al. 2000) assuming solar metallicity and solar age. For the evolved stars, a flat initial mass of 1.32 M$_{\odot}$ was adopted (however, small imprecisions in masses or radii do not influence our conclusions).

The result is in Figure : RV as a function of the mass-to-radius ($M/R$) ratio. If gravitational redshift were the only mechanism acting, one would expect a linear relationship, which is not present: the overall slope of the observed data can actually be fitted with a coefficient consistent with zero.

A closer analysis suggests that a trend among the dwarfs (e.g., $M/R >$`<!-- -->`{=html}0.5) might be present, fitted by a RV = 0.99($M/R$) law. However, given the small range in $M/R$ (a factor 2 only) covered by the main-sequence stars, the uncertainty in this slope is large, and the slope is not statistically significant. Giants, if analyzed alone, show instead a negative slope of RV vs. $M/R$, also with a large error.

## Maximum likelihood estimation of the kinematic cluster parameters

Following the approach laid out by Pryor & Meylan (1993) we performed a maximum likelihood estimation (MLE) of the kinematic cluster parameters assuming a simple model including the effect of the gravitational redshift, and separating the two stellar groups of giants and dwarfs. We are interested in the cluster systemic velocity $v_\mathrm{c}$, the velocity dispersion of giants and dwarfs, $\sigma_\mathrm{g}$ and $\sigma_\mathrm{d}$, respectively, as well as the dependence of the spectroscopically measured velocity on (M/R) – presumably due to the gravitational redshift. We assumed that each individual measurement $i$ was drawn from a probability density $$l_i(\ensuremath{v_\mathrm{c}},\ensuremath{\alpha},\ensuremath{\sigma_\mathrm{\!x}^2}) = \frac{1}{\sqrt{2\pi(\ensuremath{\sigma_\mathrm{\!x}^2}+\ensuremath{\sigma_\mathrm{\!m}^2})}}\,
\exp\left(-\frac{(\ensuremath{v_i}-\ensuremath{v_\mathrm{c}}-\ensuremath{\alpha}\ensuremath{P_i})^2}{2(\ensuremath{\sigma_\mathrm{\!x}^2}+\ensuremath{\sigma_\mathrm{\!m}^2})}\right).$$ $v_i$ is the spectroscopically measured radial velocity, $\sigma_\mathrm{x}$ the velocity dispersion of the group of stars the object $i$ belongs to, $\sigma_\mathrm{m}$ the measurement error on a single star for which we assumed it is the same for all stars and amounts to 0.3 km s$^{-1}$. For brevity, we introduced $\ensuremath{P_i}\equiv(M/R)_i$. $\alpha$ is a constant of proportionality which should come out as 0.636 km s$^{-1}$ per unit $M/R$ (in solar units) if the expected gravitational redshift was present[^2]. We assumed that there is no measurement uncertainty on the $P_i$. This approximation is justified since the uncertainties are dominated by the spectroscopic ones which correspond to an uncertainty in $P_i$ of $\sim 0.5$, much worse than the precision expected for $P_i$-values from isochrone fitting.

We maximised the overall likelihood $$L(\ensuremath{v_\mathrm{c}},\ensuremath{\alpha},\ensuremath{\sigma_\mathrm{\!g}^2},\ensuremath{\sigma_\mathrm{\!d}^2}) = \Pi_{i\in\{\mathrm{giants}\}} l_i(\ensuremath{v_\mathrm{c}},\ensuremath{\alpha},\ensuremath{\sigma_\mathrm{\!g}^2}) \,
\Pi_{j\in\{\mathrm{dwarfs}\}} l_j(\ensuremath{v_\mathrm{c}},\ensuremath{\alpha},\ensuremath{\sigma_\mathrm{\!d}^2}),$$ i.e., we assumed that giants and dwarfs share the same systemic velocity, and exhibit an intrinsic line shift proportional to their mass to radius ratio. Giants and dwarfs can potentially have different velocity dispersions – as is expected from the previous analysis. The ansatz for the likelihood above emphasises the role of intrinsic line shifts in the determination of the kinematic parameters of the cluster, and the need for correcting for intrinsic shifts if one aims at a precise determination of the kinematic properties of cluster members.

We followed the standard procedures for determining the most likely parameters and their covariance matrix, except that we used the variances instead of standard deviations as actual parameters in the minimisation process. This had mostly practical reasons since some algebraic expressions became much more compact this way. For a single group of stars with velocity dispersion $\sigma_\mathrm{x}$ the best fitting parameters could be found analytically as

$$\ensuremath{\alpha}= \frac{\ensuremath{\left\langle\ensuremath{v_i}\ensuremath{P_i}\right\rangle}-\ensuremath{\left\langle\ensuremath{v_i}\right\rangle}\ensuremath{\left\langle\ensuremath{P_i}\right\rangle}}{\ensuremath{\left\langle\ensuremath{P_i}^2\right\rangle}-\ensuremath{\left\langle\ensuremath{P_i}\right\rangle}^2},$$ $$\ensuremath{v_\mathrm{c}}= \ensuremath{\left\langle\ensuremath{v_i}\right\rangle} - \ensuremath{\alpha}\ensuremath{\left\langle\ensuremath{P_i}\right\rangle},$$ and $$\ensuremath{\sigma_\mathrm{\!x}^2}= \ensuremath{\left\langle(\ensuremath{v_i}-\ensuremath{v_\mathrm{c}}-\ensuremath{\alpha}\ensuremath{P_i})^2\right\rangle} - \ensuremath{\sigma_\mathrm{\!m}^2}.$$ $\left\langle.\right\rangle$ denotes the (ensemble) average. We did not find an analytical solution for the case of two different velocity dispersions. For this case we solved the minimisation problem numerically by direct iteration on $\sigma_\mathrm{\!g}^2$ and $\sigma_\mathrm{\!d}^2$ as suggested by Pryor & Meylan (1993). With the best fitting parameters we evaluated the covariance matrix $(-\partial^2 L
/\partial\theta_i\partial\theta_j)^{-1}$ where $\mathbf{\theta}$ is the vector of parameters $\mathbf{\theta}=(\ensuremath{v_\mathrm{c}},\ensuremath{\alpha},\ensuremath{\sigma_\mathrm{\!g}^2},\ensuremath{\sigma_\mathrm{\!d}^2})$. We obtained the best fitting parameters and their $1\,\sigma$ uncertainties as

$$\begin{array}{llrll}
\ensuremath{v_\mathrm{c}}& = &  33.824 & \pm\,0.121 &\mbox{km\,s$^{-1}$}\\
\ensuremath{\alpha}& = & -0.134  & \pm\,0.230 &\mbox{km\,s$^{-1}$}\\
\ensuremath{\sigma_\mathrm{g}}& = &  0.566  & \pm\,0.068 &\mbox{km\,s$^{-1}$}\\
\ensuremath{\sigma_\mathrm{d}}& = &  0.849  & \pm\,0.078 &\mbox{km\,s$^{-1}$}.
\end{array}$$

The errors of the standard deviations were calculated from the errors on the variances applying the standard error propagation formula. From the covariance matrix we calculated the correlation coefficients among parameter pairs: $-$`<!-- -->`{=html}0.840 ($v_\mathrm{c}$, $\alpha$); $-$`<!-- -->`{=html}0.067 ($v_\mathrm{c}$, $\sigma_\mathrm{\!g}^2$); 0.059 ($v_\mathrm{c}$, $\sigma_\mathrm{\!d}^2$); 0.093 ($\alpha$, $\sigma_\mathrm{\!g}^2$); $-$`<!-- -->`{=html}0.081 ($\alpha$, $\sigma_\mathrm{\!d}^2$); and $-$`<!-- -->`{=html}0.008 ($\sigma_\mathrm{\!g}^2$, $\sigma_\mathrm{\!d}^2$). The only significant and actually strong correlation among the parameters is the one between $v_\mathrm{c}$ and $\alpha$. This is understandable since changes in the systemic velocity can to some extent be traded against changes of the intrinsic line shifts. The velocity dispersions in each group of stars are largely insensitive to the systemic velocity and the intrinsic line shifts. All in all, the MLE analysis confirms our findings that there is a clear difference ($3.6\,\sigma$) between the velocity dispersions of giants and dwarfs, and a dependence of the intrinsic line shifts on $(M/R)$ which deviates significantly ($3.3\,\sigma$) from the expected behaviour for a signal solely due to a gravitational redshift.

# Discussion

Our data show two main results: that no obvious signature of gravitational redshift is seen in the data, and that the dispersion of the radial velocities is higher in un-evolved stars than in evolved ones. We discuss these two results separately:

## Missing gravitational redshift

As discussed above, a radial velocity shift of up to 600 m s$^{-1}$ is expected between main sequence stars and giants. Such a difference is not present in our data as illustrated in Figures  and  . Could we have missed the gravitational redshift because of some systematic effect?

### Template effects

Given that the mask used is of spectral type K0 III, a better match between template and target star is expected for the sub-giants and giants than for main sequence stars. One clear effect of using the K0 III mask with G-type turn-off stars is that the resulting CCF is shallower than for the giants and sub-giants (Figure 3 of Paper I). This is both because of the mismatch between star and template, and because of the generally weaker lines in dwarf stars; the depth of the CCF is thus related to the luminosity.

Could the mismatch between luminosity class and mask have caused a velocity shift in the turnoff stars? To investigate this, we computed new cross-correlation functions to all stars using a G2 V mask instead. The mean RV difference was found to be a well defined value (30 m s$^{-1}$) due to template zero-point mismatch. The RV difference distribution has an RMS of 120 m s$^{-1}$. A detailed inspection of spectra showing differences larger than 100 m s$^{-1}$ reveals that low S/N along with a not optimum cosmic ray cleaning are responsible. In addition, most of the discrepant spectra were recorded during the last season, when a simultaneous Th-Ar lamp was used, contaminating faint spectra through its emission lines. If only data from the first two seasons are considered, the RMS drops to 90 m s$^{-1}$, as in Fig. .

While these tests show some limits in the precision, the use of a different mask also shows that there is no offset introduced by the choice of the mask differentially affecting various groups of stars. However, an extra error of about 100 m s$^{-1}$ due to poor cosmic-ray cleaning and ThAr contamination has to be added to the formal error of the cross-correlation fit.

We also tested different mask characteristics, such as dependence of the measured RV upon the wavelength range, finding that values retrieved from the 500-700 nm mask are within a few tens of m s$^{-1}$ from those found for the 'full' 400-1000 nm one, for different types of stars. We therefore conclude that the mask and wavelength range used do not induce systematic effects which could cancel the gravitational redshift signal.

### Spectrograph drifts

Variations in temperature and atmospheric pressure during the observing night can cause spectral drifts of a few hundred m s$^{-1}$ (e.g., Platais et al. 2007). This effect is seen in our data collected in the third season, where the simultaneous calibration fiber was used. Unfortunately, the first two datasets did not make use of the simultaneous calibration fiber and the spectrograph drifts could not be removed.

In order to estimate the amplitude of the drifts caused by changes in atmospheric conditions we used 52 observations of $\tau$ Cet (HD 10700) taken from the SACY survey (Torres et al. 2006) and from our own observations. These spectra were reduced and analyzed in a similar way to the present M67 stars. They cover a time span of 3000 days from Oct. 1999 to Feb. 2008. The mean RV is $-$`<!-- -->`{=html}16.37 km s$^{-1}$ with an RMS of 0.3 km s$^{-1}$. The best fit, excluding 6 extreme values, gives a similar mean, and a sigma of 0.18 km s$^{-1}$. Thus we believe that the drift error is about 200 m s$^{-1}$.

Simulations of the effect of photon noise in the RV computation show that already with a S/N of 10–30, radial velocity errors are as low as a few tens of m s$^{-1}$ (Paper I; Bouchy et al. 2001). However those simulations only consider white (photon) noise, and (as discussed above) there are additional sources of errors associated with low S/N. In Figure  we plot the scatter diagram as a function of stellar magnitude and indeed the behavior is that expected, if the noise was not negligible.

We investigate if one can reproduce the missing gravitational trend and the scatter seen in Figure  accounting for the different sources of spread in the VR. First, the intrisic velocity dispersion of the cluster of 0.8 km s$^{-1}$ was considered. A Gaussian drift error drawn from a distribution with an RMS of 200 m s$^{-1}$ and a photon noise of 150 m s$^{-1}$ for a V=14 star were added. This noise was scaled as the square root of the flux for a given magnitude (i.e., $1/SN$, assuming that photon noise is the limiting factor for faint targets).

We then computed the difference of the center of the RV distribution of main sequence and giant stars: this is always below 100 m s$^{-1}$. Moreover, the RMS of both distributions is clearly dominated by the cluster velocity dispersion. An additional gravitational shift of 600 m s$^{-1}$ added to the distribution of sub-giant and giant stars can be retrieved in all cases in spite of their limited number.

We cannot so far identify any instrumental or observational effect which could cancel out the expected gravitational redshift signal, and we conclude that our observations show the absence of evidence for a gravitational redshift trend.

Nordström et al. (1997) briefly discuss a similar case, because they found a difference of 0.4 km s$^{-1}$ (dwarfs redder than giants) between the average RV of giants and main sequence stars in the open cluster NGC 3680, almost exactly what is predicted from gravitational redshift. NGC 3680 is an almost completely evaporated cluster, with only a few bona-fide members left. The paucity of the sample precludes certain conclusions, but could suggest that our results about M67 may not be generalized to other clusters. On the other hand, the fact that no gravitational redshift is found by analyzing the published data by Mathieu et al. (1986) shows that this result is not limited to our observations only, but it is a characteristic of this cluster.

###  Cluster dynamics

Even if open clusters are not known to rotate or to have complex dynamical topologies, we shall analyze whether the missing gravitational signature could be due to a combination of complex cluster dynamics and bias selection, namely that the observed giants and dwarfs are situated in different spatial regions of the cluster, perhaps affected by different global velocities. In this case the stars would not share the same systemic velocity and this bias could conspire to cancel the gravitational redshift signature. For M67, there is no report of complex dynamics, with the exception of a claim by Loktin et al. (2005) of a possible expansion of the cluster core. No evidence of rotation is found in the analysis of the Mathieu et al. (1986) data by the same authors. In order to investigate effects related to cluster dynamics we have analyzed the velocity distribution of the observed stars as a function of R (distance from the cluster center) and $\theta$ (angular position) for giants and dwarfs, finding no evidence of RV variations with these parameters (Figure ). We note in passing that we do not find evidence for a larger dispersion at large radii. We therefore conclude that complex cluster dynamics coupled to a selection bias cannot explain the non-detection of gravitational redshift.

### Convective shifts

Dynamic effects in stellar atmospheres may produce net shifts in photospheric absorption lines, and we cannot exclude that lines in main-sequence stars could be more blueshifted than in giants, thus (partially) canceling the expected gravitational redshift signature.

Such a differential shift need not be present in all spectral lines, but only (in a statistical sense) in the deeper ones. In fact, for a given metallicity, typical lines (metallic neutral ions) are generally deeper in cool giants than in dwarfs. By using the same mask for all stars, we use lines with different equivalent widths in the two classes of stars, and relatively deeper in the giant spectra.

Since single spectra of each star have a rather low S/N ratio, we have constructed two 'average' low-noise spectra: one from 64 selected dwarfs (called 'dwarf') and one for 20 selected giants ('giant'). The spectra were first normalized in flux at 600 nm, and then averaged by calculating the medians. Given the low velocity dispersion of the cluster stars, these co-added spectra are not appreciably broadened. The S/N ratio of these spectra is about 200 in the region of interest for the 'dwarf' spectrum and about 140 for the 'giant' one. The co-added spectra still show – within the uncertainties – the same radial velocity, for 'giant' and 'dwarf' alike. We selected a number of lines used for chemical analysis of solar metallicity stars (Pasquini et al. 2004), inspecting them visually in the co-added spectra. Out of the original list, 52 lines passed this quality check, and were fitted with a Gaussian function, using the ARES package (Sousa et al. 2008). We first tested with the IRAF task Splot the consistency between the automatic program and the manual fitting. The central wavelengths in the 'dwarf' and 'giant' spectra were recorded, and the DV = V$_\mathrm{giant}$ – V$_\mathrm{dwarf}$ is plotted in Figure as function of line equivalent widths. For 'giant' lines weaker than about 4 pm (40 mÅ), the measures were not very reliable, and such weak lines are not used. On the other hand, lines with equivalent widths above 15 pm (150 mÅ) were discarded to minimize possible contamination by blends. In both plots of Figure a trend is present: the stronger the line, the redder the shift in the giant spectrum. The observed relative shifts are up to $\pm$`<!-- -->`{=html}500 m s$^{-1}$, comparable to the expected gravitational redshift. We looked for, but did not find any dependence with other quantities (e.g., wavelength).

Formal fits to the trends give: DV=2.1\*EW$_\mathrm{dwarfs}-$`<!-- -->`{=html}350 m$s^{-1}$ and DV=5.0\*EW$_\mathrm{giants}-$`<!-- -->`{=html}418 ms$^{-1}$. The trends are only marginally significant (correlation coefficients 0.25 and 0.35). The formal error in measurements of single lines is about 50 ms$^{-1}$, and we believe that the observed spread is real. The lines used are mostly from neutral elements (dominated by Fe) with different characteristics, but also a few . With higher S/N spectra for single stars, it would be possible to further investigate trends between different chemical elements and line parameters, but we believe that for our current data this would be an over-interpretation.

Wavelength shifts induced by convective motions have in particular been studied for lines in the Sun and in Procyon, where it has been possible both to measure and to model intrinsic lineshapes under very high spectral resolution (Dravins 1982; Allende Prieto et al. 2002a); a successively increased blueshift for weaker lines is a signature of the three-dimensional photospheric hydrodynamics (e.g., Asplund et al. 2000). Also, in other cool stars, Allende Prieto et al. (2002b) and Ramı́rez et al. (2008) observed how lines from neutral species get more blueshifted as they become fainter.

Since we are able to compare dwarfs and giants of the same cluster, some information can be added: the trend seems to differ between main sequence and evolved stars. In the M67 evolved stars, the blueshift is less than in turnoff ones, and the expected gravitational redshift is seen only for a few of the faintest lines. The blueshift variation with equivalent width is larger in giants than in dwarfs, as suggested by the trend in Figure . A caveat, however, is that further parameters may enter, in particular line broadening by stellar rotation which differs among dwarfs and giants: rotational broadening changes line asymmetries and may affect the velocity signal measured by cross correlation (Dravins & Nordlund 1990).

This discussion suggests a plausible explanation for our non-detection: gravitational redshift is largely canceled because our correlation mask is dominated by rather deep lines, which are more blueshifted in main sequence stars than in giants because of photospheric convective motions.

## Convective line shifts in 3D model atmospheres

In order to investigate the apparent absence of gravitational redshift, synthetic spectral lines were computed from 3-dimensional and time-dependent model atmospheres. We synthesized 15 artificial lines, spanning the range of excitation potential and equivalent width of the observed ones. All lines were placed at a wavelength of 620 nm, which is the average of the observed lines in the range 586–672.7 nm. The 3D model atmospheres providing the thermal and kinematic information for the spectral synthesis were taken from the CFIST 3D model atmosphere grid (Ludwig et al. 2009). For representing the 'dwarf' and 'giant', we chose for each a single 3D model atmosphere. Table  summarizes their basic properties.

```latex
\begin{table*}[t]\caption{3-dimensional radiation-hydrodynamics model atmospheres.
\label{tab:3dmodels}}
\begin{center}
\begin{tabular}[t]{lllllllllllll}
\hline\noalign{\smallskip}
Model   &  \ensuremath{T_{\mathrm{eff}}}& \ensuremath{\log g}& \ensuremath{[\mathrm{Fe/H}]}& \ensuremath{l_\mathrm{x}},\ensuremath{l_\mathrm{y}}& \ensuremath{l_\mathrm{z}}& \ensuremath{N_\mathrm{x}},\,\ensuremath{N_\mathrm{y}}& \ensuremath{N_\mathrm{z}}& \ensuremath{\mathrm{N}_\mathrm{obm}}& $T$ & \ensuremath{\delta I_\mathrm{rms}/I}& $v_\mathrm{rms}$ & Modelcode\\
\mbox{} & [K]    & [$\mathrm{cm/s^2}$]& & [Mm]& [Mm]&         &     &       & [h] &   $\%$    & [km\,s$^{-1}$]\\
\noalign{\smallskip}
\hline\noalign{\smallskip}
dwarf     & 5\,782  & 4.44 & 0.0 & 5.6 & 2.27  &  140 & 150 & 12  & 1.2 & 14.4 & 2.24 & d3gt57g44n58\\
giant     & 4\,968  & 2.50 & 0.0 & 573 & 243   &  160 & 200 &  5  & 152 & 18.9 & 3.04 & d3t50g25mm00n01\\
\noalign{\smallskip}
\hline
\end{tabular}
\end{center}
\tablefoot{ ``Model'' is the
    model's name used in this paper, \ensuremath{T_{\mathrm{eff}}}\ the effective temperature,
    \ensuremath{\log g}\ the gravitational acceleration, \ensuremath{[\mathrm{Fe/H}]}\ the metallicity, $\ensuremath{l_\mathrm{x}}=\ensuremath{l_\mathrm{y}}$ the linear horizontal size of
    the square-shaped computational box, \ensuremath{l_\mathrm{z}}\ its vertical extent,
    $\ensuremath{N_\mathrm{x}}=\ensuremath{N_\mathrm{y}}$ the number of grid points in the horizontal directions, \ensuremath{N_\mathrm{z}}\ the number in vertical direction,
    $T$ the duration of the simulated time series,
    \ensuremath{\mathrm{N}_\mathrm{obm}}\ the number of equivalent frequency points considered in the solution of the radiative     transfer equation,
    \ensuremath{\delta I_\mathrm{rms}/I}\ the relative spatial white light intensity contrast at stellar disk center,
     $v_\mathrm{rms}$ the temporal and spatial average of the vertical velocity at Rosseland optical depth unity,
    ``Modelcode'' is an internal identifier of the model sequence.}
\end{table*}
```

The wavelength shifts of the synthesized lines were measured applying the same procedure as for the observed ones, i.e., by fitting a Gaussian profile. In contrast to the observations, the artificial lines were available at very high spectral resolution ($\approx 2\times 10^6$) and signal-to-noise ratio so that their shifts with respect to the laboratory wavelength could be calculated with high fidelity. For the 32 observed lines we calculated convective line shifts by interpolating in the grid of synthetic lines according to the observed equivalent width and excitation potential of the line. The typical statistical uncertainty to which the wavelength position of an individual line is predicted is (from previous experience) $~30\,\mbox{m\,s$^{-1}$}$. This means that for the difference between wavelength shift of the same line between two models we expect a typical precision of $~42\,\mbox{m\,s$^{-1}$}$. The statistical uncertainty is related to issues like how well the acoustic oscillations present in the computational box are averaged out. The example illustrates that we further expect that the line shifts among lines of the same model are correlated, hence we do not expect that the uncertainty of the predicted shifts among giant and dwarf drops as the square-root of the number of lines used to characterise it. Thus, we expect a modeling uncertainty of about $\pm 40\,\mbox{m\,s$^{-1}$}$ for the line shift difference between the giant and dwarf – not a large value but also not completely insignificant.

Figure  shows the predicted convective shifts in dwarfs and giants. Before determining the line shift, the synthetic spectral line profiles were broadened with a Gaussian profile of $\sigma$ 3.9 km s$^{-1}$ which corresponds to the spectrograph resolution, combined with the macroturbulence-like velocity dispersion stemming from the motions of individual cluster stars. In addition, it was assumed that the dwarf rotates with a projected rotational velocity of 5 km s$^{-1}$ while the giant exihibits no rotation. We find for the difference of the shifts of lines between the giant and dwarf a minimum of 62  a maximum of 211  and an average of $146\,\mbox{m\,s$^{-1}$}$. We note that the correlation between the strengths of lines in the giant and dwarf is not very tight. However, as seen in Fig. , we find a trend of the convective shift with equivalent width which is continuous between dwarfs and giants: on average the stronger lines are less blueshifted in the giant.

Figure  depicts the calculated differences of the shift of the same line between dwarf and giant, comparing them to observations. For the synthetic lines, assumptions about spectrograph resolution and stellar rotation were the same as before. An average gravitational redshift was now included, so that the giant synthetic points were blueshifted by 382  (applying an average M/R of 0.73 for dwarfs and 0.13 for giants, as for the stars used to compute the average spectra). The agreement with the observations is now quite good, with a difference of only $\sim 140\,\mbox{m\,s$^{-1}$}$ between the synthetic and the observed average shifts.

However, the observed spectra show effects not predicted by the modeling: dependence of the shift on the equivalent width, and a considerable spread in the measurements. Nevertheless, we conclude that convective line shifts are the most likely explanation as to why the expected gravitational redshift signature was not observed.

Among possible causes of remaining discrepancies between models and observations, we note that the chosen 3D models may not be the optimal representatives of the 'dwarf' and 'giant' spectra. Their choice was dictated by availability. For example, the iron abundance determined from the dwarf model is approximately solar but for the giant it turns out to be 0.5 dex greater. This may indicate that the giant model is hotter than the typical observed giant. For the dwarfs, solar temperature is appropriate for our fainter stars, while higher temperatures are expected for the stars closer to the turnoff. Indeed a similar Fe abundance is obtained by adopting 6000 and 4600 K for dwarfs and giants respectively (using 1D models). However, we do not expect the line shifts to sensitively depend on stellar parameters (as long as the line strength is matched) so we argue that our approach should give a reasonably accurate estimate of the expected line shift differences.

Another possibility is the presence of noise and biases in our combined spectra; in order to test that we repeated the analysis using the Kitt Peak solar flux spectrum (Kurucz et al. 1984) and the Arcturus spectrum (Bagnulo et al. 2003). Although here we can only test the slope and the spread around the relationship, this gives qualitatively similar results to M67, although the slope is shallower ($\sim$ 2 m s$^{-1}$/mÅ instead of 5) and the scatter is smaller (100 m s$^{-1}$ instead of 240 m s$^{-1}$).

Being the first attempt of such a test, we believe it is encouraging. In presence of higher S/N data it will be possible to test other effects, such as the dependence of the line shifts on other parameters (element, excitation potential, ionization, wavelength). One could also apply the same technique to clusters of different metallicity and chemical composition, and to work out comparison grids for theoretical models.

## Dynamically hot dwarfs

The observation that dwarfs are dynamically hotter than giants might seem counter-intuitive, because giants are known to be intrinsic RV variables, with variations of up to a few hundreds m s$^{-1}$ (e.g., Setiawan et al. 2004), and their intrinsic RV variability will add to their RV dispersion. In addition, main sequence stars have lower masses, and in case of mass segregation the giants should be (as they are in M67) mostly concentrated in the middle of the cluster, near the center of the potential well. On the other hand, in an equipartition regime, giants should be dynamically cooler, because they are more massive.

The proper-motion dispersion of the bright M67 stars as measured by Girard et al. (1989) is almost 0.2 mas/yr. At a distance of 836 pc (Pasquini et al. 2008) it corresponds to a radial velocity dispersion of 0.8 km s$^{-1}$, comparing very well with the dispersion of RV=0.83 km s$^{-1}$ measured from all stars and also with the 0.68 km s$^{-1}$ dispersion of the giants only. The comparison of the radial velocity distributions for both populations in Figure  suggests that the dwarfs are dynamically hotter than the giants in M67, i.e., they have a larger velocity dispersion.

In order to check whether this might originate by some observational bias, we performed the following test: we again simulated the radial velocity distribution for the dwarfs using the same error parameters and a population of 100 stars as in Figure . The RMS of this artificial distribution was then computed, and the procedure repeated 5000 times. The result is a distribution of RMSes, centered at the velocity dispersion of the cluster (800 m s$^{-1}$) with values between 0.6 and 1.0 km s$^{-1}$. If in the simulations we now turn off the spectrometer drift error (Sect. 6.1.2), while maintaining all other parameters, we end up with similar figures. This already indicates that the observed velocity dispersion does not depend on measurement errors, but rather on the number of objects.

As a further test, we set the drift error to its initial value and generated the RMS distributions for two populations: one with N=100 members (corresponding to our real case) and another with 1000 members. The result is in Figure . With 1000 members, the velocity dispersion is well defined with a center value of 0.8 km s$^{-1}$ and 0.73 to 0.88 km s$^{-1}$ peak-to-peak. Therefore, the observation that the M67 dwarfs have a larger velocity dispersion than the giants seems not caused by our measurement errors but possibly is an effect of low number statistics.

We note that Nordström et al. (1997) found a similar result for NGC 3680, with the dwarfs showing a dispersion of 0.67 km s$^{-1}$ compared to 0.36 of the giants. Also, analysis of the Mathieu et al. (1986) data (which, however, contain a limited number of main sequence stars) is in perfect agreement with our findings, giving $\sigma(RV)$= 0.89 and 0.66 km s$^{-1}$ for dwarfs and giants respectively.

Since we cannot exclude that the effect is real, it is worth mentioning that a similar result was found by Zhao et al. (1996) and Loktin et al. (2005) from analyzing proper motion data. Zhao et al. (1996), analyzing a sample of about 200 stars, found that the dispersion of proper motions of M67 bright stars (V$<13.5$ ) is significantly smaller than that of fainter objects (13.5 $<$ V $<$ 14.5), by a factor $\sim$`<!-- -->`{=html}1.3. They also found that the brighter (more massive ) stars are more concentrated towards the center of the cluster. Also Loktin et al. (2005) found a similar result, with a much larger sample. The larger sample allows them to subdivide the color-magnitude diagram into some sub-sample with giants, and another with turnoff and upper main sequence stars. These two subsamples closely correspond to our division between giants and dwarfs. For dwarfs they find dispersion velocities on average almost 1.3 times higher than for giants (their Table 3).

We also note that the tangential and radial velocity dispersions derived by Loktin et al. (2005) are substantially larger than our radial-velocity dispersions: 0.84 and 0.77 km s$^{-1}$ for the giants; 1.01 and 1.04 km s$^{-1}$ for dwarfs ('radial' in their work refers to polar coordinates in the sky plane). Loktin et al. (2005) adopt a distance to the cluster of 904 parsecs. Using our preferred distance of 831 pc, the radial and tangential dispersions become (0.77, 0.71) and (0.93, 0.96) km s$^{-1}$, respectively for giants and dwarfs, fully compatible with the observed radial velocity dispersions of 0.68 (giants) and 0.9 km s$^{-1}$ (dwarfs). This result indirectly confirms the shorter distance scale for this cluster that emerges from recent works (cf. Pasquini et al. 2008).

It is remarkable that in both our data and those by Mathieu et al. (1986), the ratio between the radial velocity dispersions of main sequence and evolved stars is 1.3, the same as found by Zhao et al. (1996) and Lotkin et al. (2005) for the dispersion between fainter and brighter stars. Zhao et al. (1996) interpret these results claiming that the dynamical energy of the more massive stars has been passed in the time evolution of the cluster to the less massive objects, which are now escaping the cluster.

While a full dynamical analysis of the cluster is beyond the scope of this work, the possibility of joining precise radial velocity and proper motion data seems attractive for detailed future investigations.

M67 is by now the best studied open cluster when combining both, the number of stars observed and the precision of the measurements. Substancially larger samples of radial velocity measurements have been gathered in globular clusters (see e.g. Sommariva et al. 2009), but they are dynamically much more complex objects and with a larger velocity dispersion.

# Conclusions

A search was carried out to find the expected differences in gravitational redshift between the spectra of dwarf and giant stars in the M67 open cluster, whose stars should share the same systemic velocity, and where the average radial motion for different types of stars thus should produce the same Doppler shift. However, despite adequate measuring precision, the expected signature was not seen. Investigations of the reasons for absence of signal – in particular following high-resolution spectral synthesis from 3-dimensional hydrodynamic model atmospheres – indicate that likely reasons are different amounts of convective blushifts in different types of stars. An enhanced convective blueshift in, particularly, those stronger spectral lines in dwarf stars that carry a larger weight in the cross-correlation techniques used for radial-velocity measurement, appear adequate to cancel the expected gravitational redshift.

This study thus also illustrates the limitations of cross-correlation techniques for high-accuracy lineshift measurements. To understand the origin of wavelength shifts that are only a small fraction of any spectral-line width, a cross-correlation measure over an extended spectral range may not be adequate, but rather one needs to examine (groups of) individually selected line profiles measured with sufficient spectral fidelity to reveal also their intrinsic asymmetries and amount of rotational broadening. These can then be compared to synthetic spectral lines from 3-D atmospheres complete with convective blueshifts and gravitational redshifts. Not many current spectrometers have adequate performance for such tasks, but work is in progress towards realizing such high-fidelity spectroscopy with an accurate wavelength calibration at the largest telescopes (Dravins 2010; Pepe et al. 2010; Pasquini et al. 2010). With these instruments it will be possible to extend the present study to less explored stars, for instance in the metal-poor and in the metal-rich regimes.

[^1]: Based on observations collected at ESO, La Silla, Chile, during the agreement between the Observatorio Nacional at Rio de Janeiro and ESO.

[^2]: Here, we neglected the 3 m s$^{-1}$ redshift in the gravitational potential of the Sun at Earth orbit.


abstract: Spam costs US corporations upwards of \$8.9 billion a year, and comprises as much as 40% of all email received . Solutions exist to reduce the amount of spam seen by end users, but cannot withstand sophisticated attacks. Worse yet, many will occasionally misclassify and silently drop legitimate email. Spammers take advantage of the near-zero cost of sending email to flood the network, knowing that success even a tiny fraction of the time means a profit. End users, however, have proven unwilling to pay money to send email to friends and family.
  .
  We show that it is feasible to extend the existing mail system to reduce the amount of unwanted email, without misclassifying email, and without charging well-behaved users. We require that bulk email senders accurately classify each email message they send as an advertisement with an area of interest or else be charged a small negative incentive per message delivered. Recipients are able to filter out email outside their scope of interest, while senders are able to focus their sendings to the appropriate audience.
author: Evan P. Greenberg and David R. Cheriton
Stanford University
{evang,firstname.lastname@example.com
title: Enforcing Bulk Mail Classification

# Introduction

Unsolicited email has become a real problem. Companies around the globe are losing billions of dollars each year in the form of lost productivity. And the spammers are profiting; 8% of Internet users bought something from spam in 2003 .

Unfortunately, the numbers in this game favor the spammers. A marketer can send billions of emails at negligible cost; to realize a profit, only a handful need result in a sale. Every email which makes it through our arsenal of spam filters carries with it a very real cost to the recipient, who must spend time to identify and discard the message. Spammers have strong incentives to evade or subvert any mechanism we deploy, while the recipients bear the cost of dealing with any email that leaks through.

Antispam solutions generally attempt either to reduce the leakage by improving the quality of filtration or to reduce the profit potential by imposing a cost for sending spam. Approaches which attempt the former are in the unenviable position of not being able to safely err on either side in their classification, as they risk either being ineffectual or dropping desired email. In this paper, we will primarily take the latter approach, of introducing penalties for sending spam, and we will show that in our system:

-   Legitimate email is never misclassified or dropped.

-   Well-behaved users do not pay to send email.

-   Legitimate mass mailing is supported.

Clearly it is also necessary to provide attack resistance, assuming sophisticated spammers who will optimize their strategies for any solution we deploy.

This is accomplished by requiring bulk email senders to accurately classify each email message they send as an advertisement and its area of interest or else be charged a small negative incentive per message delivered to a recipient, a micro-penalty. This micro-penalty multiplied by the scale required for spam makes it uneconomic to lie. Reliable classification makes it feasible for email recipients to filter out advertising emails that are outside their scope of interest. We show that it is feasible to log all email so that violations can be acted on quickly and securely while resisting attacks, supporting opt-in mailing lists as well as direct person-to-person email, and requiring only modest changes to existing email servers while being economic to deploy to Internet scale.

# Tagging Bulk Mail

We require senders to annotate all bulk email with a new header field, X-Bulk-Mail. At this time we define two types of bulk mail: advertising and mailing lists (described further in Section ). Advertising emails are tagged with an ADV leader followed by a comma-separated list of interest groups. Interest groups are defined hierarchically, with terms separated by dots, and can be defined organically as the system evolves. For instance, an advertisement for suntan lotion might include the header

X-Bulk-Mail:

:   ADV: rec.sports.swimming,
    rec.sports.sailing

A user can easily establish filters to discard mail for categories he isn't interested in. Ideally this would be provided as a service by his ISP.

When email received is considered to have been misclassified, whether because it is lacking the appropriate header or because it has been designated with an inappropriate interest group, the user will forward the email to his ISP using his email client's bounce or redirect command. This preserves the full headers of the email, allowing the penalty phase of the system to take effect.

In our system, we define personal email to be any email which does not carry an X-Bulk-Mail tag. It is treated in the same way as bulk mail, and is presumed to be of interest to all recipients. Email which is not of interest, such as advertising, can be objected to on the basis of misclassification.

# Penalizing Misbehavior

Each email server is extended to log relevant header information from the mail it forwards over a window of two weeks. Specifically, we store the (cryptographically secure) hash of the Date, To, From, and Received fields.

When an ISP receives a spam complaint, it can perform a lookup to validate that the complaint references a genuine email. It then forwards the complaint on to the previous relay in the path. It also requires the upstream relay pay it the required micro-penalty. That relay will recover the cost when it passes the complaint on.

The result of this process is as if the micro-penalty were levied against the ISP of the sender by the ISP of the recipient. The sender's ISP can choose how to deal with the complaint, but should at a minimum include forwarding the complaint to the user and imposing some form of sanction.

For this system of penalties to work, we must also establish a few procedural elements. First, as stated before, we have a history of forwarded email to work from to prevent fraudulent use of the penalty system. We also refuse to forward email that is more than a week old. An enforceable penalty process must be set up between SMTP peers, and relays must refuse unknown connections. In section we discuss how an ISP might redirect unknown relays to more permissive ingress points. Finally, ISPs must specify some form of ingress rate limit, to bound their potential outstanding liability. The monthly service fee paid by end users could be considered a bond against which the user borrows to send email.

The impact of these constraints on typical clients is negligible. Consider an ISP which charges \$30/mo for Internet access, and a micro-penalty fee of \$0.10. The ISP limits customers to sending 100 emails per week, and terminates accounts on receipt of 10 spam complaints over a 3-month period. In this scenario, the ISP has a maximum outstanding liability of \$20.90 per client, which is safely less than the monthly service fee. The ISP could also choose to allow clients to reset their spam counter in exchange for a \$1 fee. Clients would need to send more than 100 emails in a one week or send 10 offensive emails before noticing the constraints of the system. Upgraded accounts with an explicit surety could be provided to address these issues for the handful of customers who require greater flexibility.

To the spammer, however, even a \$0.10 fee is profound. As a reference point consider DoubleClick, which provides legitimate marketers with the tools needed to maintain high-quality lists of interested users. For these lists in Q1 of 2004 , DoubleClick customers saw a revenue of \$0.23 per email. Few enterprises could afford to spend 40% of their revenue on marketing, and few spammers will have returns anywhere close to that seen by the carefully culled lists maintained by professional marketing companies.

# Attack Resistance

For clarity of presentation, we have so far ignored the possibility of malicious entities. It is critical, however, to demonstrate that attackers cannot evade the system to send email without being held accountable, and that they cannot subvert the system to charge innocent users.

## Assigning Responsibility

Suppose Alice receives spam, for which she forwards a complaint to her ISP. There are only three possible attackers: the sender, some host on the path, and the receiver (her ISP). There is no way to send email in our record route environment without being on-path (unless there is an on-path accomplice, which is itself an attacker).

Assuming that there are no misbehaving hosts on the path between Alice and the sender, it is clear the spammer will be charged. Using a fake return address will not help, since the complaint will use the unforgeable portion of the Received header path to reach his ISP, which has a clear economic incentive to correctly identify the true sender as well as the means to do so, since it controls where and how email enters its network. The complaint cannot be forwarded by the sender to another host because it will not be present in the mail history for that node.

Suppose there is some host on the path between Alice and the sender which is misbehaving. It cannot refuse to validate the spam complaint without jeopardizing its relationship with its downstream peer. The complaint will, in turn, only be accepted by the upstream ISP towards the true sender. Spam forged by this host would result in the host receiving complaints which could not be forwarded to any "upstream" peers.

This misbehaving host could hijack legitimate mail traffic by replacing the body of the email. This would be detected by the endpoints of the communication channel, and eventually the culprit uncloaked. This more pernicious attack also requires the malignant host be on the path between Alice and some other user she wants to talk with, which is already unlikely given the short length of typical SMTP relay paths and the relative trustworthiness of genuinely on-path hosts.

Finally, Alice's ISP could itself send her spam, or refuse to take action against spammers within the local network. The free market suggests that such ISPs will not last very long, as users who find this behavior distasteful will simply take their business elsewhere.

## Protecting Innocent Users

We also prevent malicious nodes from abusing the spam reporting system itself. In this case, the attacker could be an endpoint, some ISP on the path between two nodes, or an off-path entity.

A client can only complain about an email she actually receives because ISPs validate all complaints against their histories. This also automatically excludes reverse path forgery and makes it easy to prevent recipients from filing more than one complaint per spam.

A malicious on-path node is prevented from forging complaints for the same reasons as for an endpoint. It would be possible for the node to generate spam complaints about email which it was only supposed to forward. However, because the sender is informed when complaints are lodged, this type of misbehavior can be easily detected.

A true off-path attacker would need to somehow derive the headers of emails traversing the network in order to file spam complaints. Guessing these would not be feasible, but the attacker might try to snoop email sessions or have an on-path accomplice forward the relevant data. However, it still would not be possible for the attacker to introduce a spam complaint into the network at any of the well-behaved relays, as these nodes would know the attacker was not a suitable next hop for the original email.

## Other Issues

It is conceivable that users might accidentally file spam complaints against desirable email. Anecdotal evidence from the SPAM-L mailing list indicates that some of the spam reports generated in AOL's Feedback Loop are a result of just this effect. In general, lower rates of spam are likely to lead to less accidental complaints. Simple interface improvements, such as requesting confirmation or password entry may also help to further reduce this effect. Finally, our future work on conversational indemnity should also help reduce the impact of accidental reports.

The issue of zombie hosts is also of concern, because it violates the assumption that email sent by a host is intentional by the user. However, the amount of damage a single host can do is severely constrained due to the requisite ingress rate limiting. Furthermore, the threat of losing email privileges gives users a real incentive to notice and promptly clean up infections. Finally, consumer ISPs could provide as a service to customers some level of early notification and/or prevention through signature analysis or other monitoring technique. In the end, we must hold a user accountable for the behavior of his computer.

# Opt-in Bulk Mail

Mailing lists are special cases of bulk mailers, because although they may need to send a number of messages each day to thousands of subscribers, they have explicitly been allowed to do so by the process of subscription. In light of this, it does not seem practical for the mailing list to assume liability for every subscriber. Nor does it seem fair to charge the original sender for every ultimate recipient, when the sender cannot necessarily determine in advance who those individuals are. However, it is also clear that mailing lists are a widely used, and therefore desirable, feature of the current email system.

To solve this we introduce recipient-side whitelisting. When an opt-in bulk sender (such as a mailing list) sends email, it flags it with a special X-Bulk-Mail header, which consists of a LIST leader followed by an identifier for the list. The identifier could either be derived from the list name or could be a random nonce. For instance,

X-Bulk-Mail:

:   LIST: freefood.348290

SMTP relays along the path do not include entries in their histories for such email. The final relay, however, will drop any email marked as list mail which is not explicitly whitelisted by that recipient. Whitelisting is done based on the list identifier and, optionally, the reverse path to the mailing list's remailer.

This mechanism works in exactly the way we expect mailing lists to operate. The act of subscription adds the user's email address to the list of addresses on the list and also adds the list to the local whitelist. Unsubscription is simply removing the mailing list from the local whitelist. Additionally, we gain the ability to force withdrawal from a list by removing the local whitelist entry.

From the viewpoint of the list operator, opt-in lists in our system carry the benefit of guaranteed protection from liability, because users cede their ability to complain by virtue of subscribing to the list. List operators could abuse this power by filling the list with spam, but the ability to force unsubscription means that users only pay temporarily for their mistakes.

# Implementation

As a reference implementation we have produced an SMTP proxy intended to be placed in front of a site's existing infrastructure. From this vantage point it is straightforward to inspect all messages entering the network from an SMTP peer and to log the requisite hash. Similarly, messages from internal users destined for the outside world can be logged and rate controlled.

This is also a suitable location to deal with spam complaints. When a user receives an undesirable message, he simply redirects the message (this preserves the original headers) to spamsink@localdomain, and the proxy intercepts this message. It strips off the downstream headers, validates the message against its logs, and then forwards the complaint upstream. In particular, note that it is not necessary to log at every SMTP host; within a single domain of trust, only one log point is necessary.

To avoid duplication of configuration information, our proxy is tightly coupled to a slave SMTP server. Commands are passed directly through to the downstream server, except where needed to preserve the semantics of the SMTP protocol. The proxy can then determine whether a command was successful by inspecting the server's response, and takes that into account when updating its state. The upshot of this is that the proxy server can depend upon the real server to enforce any policies the site may have in place (for instance, the set of addressable destinations) without needing any separate logic or configuration. The only notable exception to this is user authentication, which the proxy needs to know about in order to perform rate limiting. But it only makes sense that the proxy should provide this, as existing servers lack the functionality.

The exact performance characteristics of the proxy will, naturally, depend upon the load characteristics of the site it serves. We assume an ISP such as Princeton University (which publishes its SMTP traffic statistics online). A generous interpretation of the statistics available shows that in the month of February, 2005, Princeton saw at most a few thousand SMTP messages per minute. This figure is roughly in line with CAIDA measurements of 1.6 million SMTP sessions per day .

If we assume a constant rate of 3000 messages per minute, or 60 million messages over the two-week window, our solution requires about 3 gigabytes of storage for logging. Our experiments show that storing the log entry for each email costs approximately 3.1msec on a single-proc P4 3.2GHz with 10k RPM SATA drive; for the suggested load this means the SMTP relay (or a dedicated logging machine) will spend 15% of its time dealing with hashes. Checking a spam complaint against the database of log entries costs just 0.77 msec. Computing the hash for an email costs just 39 usec.

Full performance results, including analysis of performance under load bursts, are pending. However, the results shown here do not include any attempts to optimize the logging mechanism; each write results in synchronous disk I/O. Note that it should be possible to maintain the entire 3GB of history in memory for reasonable cost, with periodic writes to disk just to prevent data loss in the event of a crash. Writing the whole hash table to disk takes less than 5 minutes, and could of course be done in stages.

# Incremental Deployment

Although we present this work in a context which assumes global deployment, it is also incrementally deployable, by which we mean that there is a benefit realized by ISPs who deploy this system, even if all of their neighbors do not. Specifically, an ISP who deploys our proposal can guarantee that it will never drop email for its customers, while also moderating the risk that this will cause an increase in spam load.

Assume some ISP chooses to deploy our system. It will then have two categories of neighbors: those with whom it can negotiate contracts, and those with whom it cannot. When an ISP can set up the contractual arrangement we specified earlier, it need not care if the peer ISP has also deployed the system or not, because that ISP has already agreed to be accountable for all the email it sends.

When an ISP cannot set up a contractual relationship with another ISP, it must instead treat that neighbor as a potential spammer. Specifically, it should keep track of the amount of spam received from that specific ISP so that exceeding some threshold (say 100 objectionable emails in a month) will cause the peer to be cut off. Future email from that peer can then be bounced, with an indication that the ISP has been identified as a source of spam and a few remedies suggested (such as switching to a compliant ISP or asking your ISP to become compliant). Our deploying ISP does assume financial responsibility for any email accepted from this peer ISP, particularly if it is not destined for one of its direct customers. It is important to realize that there is no requirement that one must accept email from anyone on the Internet; as such these peering relationships can be pre-screened to eliminate known spammers.

Additionally, during this transitional phase, ISPs can continue to use traditional antispam techniques such as content-based filtration or blacklisting. These techniques are not particularly useful in a world where everyone has deployed our system, but they can be applied against these "unaccountable" peer ISPs. However, it should be understood that to avoid silently dropping email, an indicative error response must be sent back if the traditional methods cause the email to be classified as spam.

In this paper we treat the SMTP reverse path is if there could be any number of ISPs on the path between sender and recipient. In practice, of course, the sender ISP will contact the recipient's ISP directly (by performing a lookup of an MX record for the destination domain). But it is impractical to assume that every pair of ISPs will be able to negotiate (and enforce) the contractual relationship we require. Rather, deployment of our system will result in baby ISPs signing contracts with larger for the ability to participate in the global exchange of email. Conveniently, there is enough flexibility in the system of MX records to allow automatic traversal of the tree (or several trees). When an unknown peer attempts to send email, the recipient's server will simply terminate the session, at which point the sender will happily try the next entry in the list. In this manner, email ingress gradually moves from less permissive to more permissive, eventually settling on the appropriate entry point.

# Related Work

A number of different solutions to the spam problem have been suggested; some have even seen widespread deployment. Unfortunately, the existing literature has not accurately gauged the sophistication of the attackers, leading to solutions which are not effective enough to stem the ongoing flood of spam.

Content-based filtering solutions, particularly those based on machine learning, have been the most effective anti-spam treatment to date. But they will never be able to guarantee that they will not decide to randomly drop legitimate email. This is an anathema in a world where we are trying to push to five-nine reliability and beyond. Worse yet, because these solutions will always leak some amount of traffic, they fail to address the core issue that spamming will continue to be profitable so long as a trickle of bandwidth exists. In general, not only is content-based filtering dangerous because it leads to false positives, it continually becomes less effective as spammers learn how to construct emails to get around the filtering algorithms.

More than a decade ago, Dwork and Naor introduced the idea of using computational stamps to price email . This idea has been respun in a variety of different clothes, including a recent paper by Balakrishnan and Karger . Unfortunately, such solutions are unlikely to ever succeed, because users are simply not willing to pay for a service they view as being free. And if the fee is purely computational (as suggested by the original work), then it is too easy to farm the computation off to unknowing users through the use of a botnet, cooperative spyware, or even a background javascript task on a popular webpage. Advances in reverse Turing tests such as CAPTCHAs , though appealing in concept as a method of postage stamping will suffer from the same popularity problems. Additionally, spammers have apparently found ways to offload this processing by using the puzzle images as conditions for admission to porn sites .

There is also a continual effort from within the Internet community to maintain various blacklists, which are then used to drop spam. In theory these should never produce false positives, as only mail validated as spam is added to the list. But they cannot prevent letting some spam through, and thus are not a true deterrent. Additionally, lists driven by Internet telescopy will only work so long as spammers are not aware of the set of trigger addresses. Lists driven by voluntary reporting are even scarier; eventually spammers will find a way to pollute these lists with false reports, at which point claims of never dropping desirable mail go out the window. Explicit whitelist solutions such as ChoiceMail provide an interesting twist. They do provide the additional benefit of requiring email to have a legitimate return address. However, the whitelist request mechanism can itself be leveraged to spam. Reverse Turing tests, though better suited to this application than email stamping, still suffer the same issues of being sharable.

The sheer scale of the spam problem on the Internet today, coupled with the technological savvy exhibited by spammers, leaves us with little choice but to introduce accountability into the network, and then establishing a system of economic disincentives for abusers. Our solution does this by introducing logs within the network. The Bonded Sender and SHRED systems suggest achieving this goal through the use of what are effectively prepurchased cryptographic stamps belonging to some trusted third party. The Bonded Sender solution, though, because it only tries to target (legitimate) bulk mailers, is only able to solve part of the problem. The SHRED proposal, recent research by Krishnamurthy and Blackmond, does seem more promising. It is our belief that by providing validation from within the network, we are able to achieve better attack resistance. In particular, we demonstrate in this paper that the economic disincentives of our system will always work against spammers; solutions requiring trusted parties are likely to find it difficult to provide the same level of assurance. Additionally, by providing spammers with the ability to achieve legitimacy by classifying their email, we simultaneously reduce the load on the penalty system and give users better automated visibility into the email they receive.

# Future Work

One way in which this system could readily be extended would be to provide conversational indemnity. That is, include the ability to negate a spam complaint by presenting evidence (in the form of a prior email) that the other party initiated communication. This could be validated by the ISPs at each end of the communication chain using the same hash history table which we use to validate spam complaints.

This sort of conversation-based indemnity could also be used to reduce the effect of accidental complaints. Additionally, this mechanism could be used to make it safe to set up an email autoresponder (along the lines of the unix *vacation* program) without turning your computer into an attack reflector. It might also be possible for ISPs to use this mechanism to relax ingress limits on users, since a response to an email would constitute an acknowledgment that the original email was not spam.

The challenge here is in ensuring that an attacker cannot leverage the indemnity mechanism to evade penalty (for instance, by playing with the 2-week history window or reusing the same refutation multiple times). Additionally, the ISP should be able to automatically refute a complaint on behalf of the user, but we do not want the ISP to maintain copies of all of a user's email.

# Conclusion

We describe a system in which bulk email senders are required to accurately classify each email message they send or be charged a small negative incentive for each recipient. A history of forwarded email receipts maintained by relays along the path provides automated verification for spam complaints by recipients, which are forwarded towards senders to provide a robust, attack-resistant penalty mechanism.

We show that this system allows person-to-person and opt-in bulk email for free, without ever misclassifying or dropping email, needs only modest changes to the existing SMTP infrastructure, and remains scalable to Internet scale. Novel contributions include the use of service fees as sureties for typical consumers, the separation of and improved protection for inherently opt-in mailing list bulk mail versus opt-out lists, and the enforceability of voluntary classification by advertisers.