the formation mechanism of the thick disk is one of the most noticeable grey points of our current understanding of the galactic formation and evolution process .
this uncertainty is particularly unfortunate because , being the thick disk a ubiquitous feature among spiral galaxies @xcite , its formation must be a common stage in the early evolution of late - type galaxies . during the nearly three decades since its discovery @xcite ,
many models have been proposed .
the merging scenario , in which the early merging of a small satellite galaxy heats a primordial galactic disk producing an old , thick , and kinematically hot disk population @xcite , has been the most favored model in the last decade , following the evidence that the thin and thick disks are chemically distinct @xcite .
nevertheless , even this scenario is not free of problems ( see , for example , * ? ? ?
* ) , and alternative models have recently been drawn to attention ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
it is clear that the general properties of the galactic thick disk , such as its mean metallicity or mean kinematics , are not enough to distinguish between the models proposed for its formation .
moreover , the merging scenario has many free parameters , such as the mass of the merging satellite and the inclination of its orbit with respect to the galactic disk , and the observations must constrain them if the quality of the model is to be finally assessed . in the last few years
, theoretical simulations have started to cast predictions of the stellar distribution of stars in the multi - dimensional spatial - kinematical - chemical space ( e.g. , * ? ? ?
* ; * ? ? ?
for example , @xcite have shown that , within the merging scenario , the time of impact , and the mass and orbit of the infalling satellite can be deduced from the distribution of the expected kinematical clumps in the @xmath7@xmath8 plane , while @xcite , @xcite , and @xcite find that the @xmath4/@xmath5 ratio and its variation with galactocentric distance are excellent indicators of the inclination angle of the merging event . @xcite and @xcite have also shown that the observed decrease of rotation velocity with distance from the plane points to a low - latitude merging .
at the same time , the observations are gathering information about the spatial variations of the chemical composition and velocity distribution ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) .
detailed results of this kind are strongly needed because , through comparison with the expectations from the different models , they can be key to solve the puzzle of the galactic thick disk formation .
we are undertaking an extensive survey aimed to reveal the kinematical and chemical vertical structure of the galactic thick disk @xcite .
preliminary kinematical results were presented by @xcite and @xcite , while the sample was also used to investigate the galactic dark disk @xcite , and lithium - rich stars in the galactic thick disk @xcite . in this paper , we focus on the kinematical results , studying the trend of kinematical quantities as a function of distance from the galactic plane . in later papers of this series , the collected spectra will be used to measure the metallicity and chemical abundances of the sample , to study the variation of the thick disk chemistry with galactic height .
our investigation is based on the sample of @xmath61200 red giants defined by @xcite , vertically distributed with respect to the galactic plane in a cone of fifteen degrees radius centered on the south galactic pole .
the sample covers a large range of galactic heights , from the plane to @xmath95  kpc , and it is volume - complete up to @xmath10=3  kpc .
all the objects have 2mass photometry @xcite , and absolute proper motions from the spm3 catalog @xcite .
the sample was defined by the color cut 0.7@xmath111.1 , to isolate intermediate - metallicity stars .
main sequence ( ms ) dwarfs were excluded both by a sloped cut at fainter magnitudes , parallel to the ms , which excludes all but the nearest ( d@xmath1263  pc ) dwarfs , and by conservative kinematical criteria imposing a stellar velocity lower than the local escape velocity ( see * ? ? ? * for more details ) .
l l c l c run & date & n. nights & instrument & n. stars + coralie1 & 2005 , september & 4 & coralie@euler & 39 + feros1 & 2005 , september & 6 & feros@2.2 & 183 + dupont1 & 2005 , october & 6 & echelle@dupont & 168 + coralie2 & 2005 , october & 3 & coralie@euler & 25 + feros2 & 2006 , august & 6 & feros@2.2 & 161 + dupont2 & 2006 , september & 6 & echelle@dupont & 119 + mike1 & 2006 , aug.nov . &
5x0.5 & mike@magellan & 77 + mike2 & 2007 , october & 2 & mike@magellan & 52 +    we collected high - resolution echelle spectra for two - thirds of the girard et  al.s sample ( 824 stars ) .
all the stars in the brighter half of the sample were spectroscopically observed , while the completeness decreases with magnitude for fainter objects .
we found a high contamination ( @xmath635% ) by nearby dwarf in the faintest end of the distribution , hence stars with @xmath13 were given lower priority after the first observations .
the distribution of proper motions and colors of the observed sub - sample is shown in figure  [ f_sample ] .
the comparison with the complete sample reveals that no selection effect was introduced : the slight overabundance of red stars in the observed sample is due to the higher completeness at brighter magnitudes , where stars are on average redder .    as discussed later
, we will restrict our investigation to stars with galactic height @xmath141.3  kpc , to avoid a strong thin disk contamination .
this reduced the sample under study to 514 stars .
we visually inspected all the spectra , identifying 46 probable dwarf stars .
as expected , they show on average low radial velocities but high u and v components , and only one is found at @xmath159 .
we also found 22 stars with [ fe / h]@xmath161.5 , as deduced from comparisons with synthetic spectra .
they were considered probable halo contaminants and were excluded from further analysis .
we note that the metal - poor tail of the thick disk extends to much lower metallicities @xcite , but below this limit thick disk stars are outnumbered by halo members by a factor of nearly ten @xcite . finally , we also excluded from analysis 34 stars that , in the galactic cylindrical coordinate system , had velocity components ( defined in section  [ s_disp ] ) outside the range @xmath17u@xmath18300  km s@xmath3 , @xmath19v@xmath12300  km s@xmath3 , and @xmath17w@xmath18150  km s@xmath3 .
the cut in @xmath20 was aimed to remove the residual halo contamination , while the other components were restricted to exclude stars with wrong distance or bad proper motion .
the cut in @xmath8 was offset toward negative values to take into account the vertical shear ( see section  [ s_resmean ] ) .
our final sample thus comprised 412 stars .
the interstellar reddening e(@xmath21 ) was derived for each target from the @xcite maps , and the @xmath22 magnitude and ( @xmath23 ) color were de - reddened by means of the transformations of @xcite .
the distance to each star was then estimated through a color - absolute magnitude relation calibrated on 47tucanae , a disk globular cluster @xcite whose stellar population is very similar to the galactic thick disk both in age and metallicity @xcite . the fit of the cluster red giant branch yields the relation ( l. miller 2008 , priv .
comm . ) : @xmath24 which were translated into absolute magnitude and de - reddened color assuming @xmath25=13.50@xmath260.08 and e(@xmath21)=0.024@xmath260.004 for the cluster distance modulus and reddening , respectively @xcite .    inspecting the 2mass photometric data of 47tuc
used to derive equation  ( [ eq_dist ] ) , we found that the uncertainty on @xmath27 is of the order of @xmath60.2 magnitudes .
this is only a marginal contribution to the final error in distance , because the main source of uncertainty is the relatively wide range of age and metallicity covered by thick disk stars .
indeed , the 2mass photometric errors have only negligible impact , because they do not exceed 0.03  magnitudes for our fainter targets ( @xmath2810.6 ) .
we estimated the effect of the age and metallicity distribution on the derived absolute magnitudes by means of yale - yonsei isochrones @xcite .
we assumed a scatter of 0.3  dex in metallicity , that should include the bulk of thick disk stars @xcite when excluding the scarcely - populated low- and high- metallicity tails @xcite , and a scatter of 2  gyr in age @xcite .
they were considered uncorrelated , because the age - metallicity relation for the thick disk is very weak @xcite .
we finally estimated the error on distance to be @xmath620% , quadratically summing all the relevant contributions .
thin disk stars do not follow the age and metallicity distribution assumed to derive equation  ( [ eq_dist ] ) , and their distances should be systematically biased . indeed , younger , more metal - rich red giants are intrinsically fainter than our estimate . comparing the absolute magnitudes
calculated from equation  ( [ eq_dist ] ) with yale - yonsei isochrones between 2 and 8  gyr , and metallicity following the age - metallicity relation of @xcite and @xcite , we found that the distance of thin disk stars would be overestimated by 1020% .
this systematic error is small , but it has a relevant consequence on the contamination of the sample , which is artificially stretched to larger heights from the plane .
our sample is contaminated by the thin disk , because the targets were selected through a color cut that efficiently eliminates low - metallicity stars , but which excludes only a small fraction of metal - rich objects , as it can be deduced from figure  1 of @xcite . in figure  [ f_contam ]
we plot an estimate of the fraction of thin disk stars in the sample , assuming 300 and 900  pc as thin and thick disk scale height , respectively , and a thick disk local normalization of 12% @xcite .
we also calculated the contamination in the case of a distance overestimate of 15% for thin disk objects .
the curves can be considered an upper limit , because the color cut , unaccounted for in the calculation , could have removed part of the contamination . in the same figure
we indicate the results of a rough empirical estimate , obtained fitting a double gaussian to the distribution of the vertical velocity component ( @xmath20 ) in five bins of 50  stars .
the dispersions of the two gaussian components were kept fixed , and the only free parameter of the fit was the number of thin disk stars in the bin .
the error bars show the results when varying @xmath5 in the range 14 - 16 and 35 - 40  km  s@xmath3 for the thin and the thick disk , respectively .
this test was performed only for @xmath291 - 2  kpc , where both populations contribute by more than 20% of the objects , i.e. 10 stars in each bin .
the empirical errorbars agree well with the curve obtained when the distance bias is considered , although the observed contamination fall - off with @xmath1 seems steeper than the modeled one .    in conclusion ,
our sample is contaminated by thin disk stars , whose distances are overestimated .
this affects even the estimate of their kinematics , which is also biased .
therefore , we will not analyze the results for the thin disk , and we will limit our study to @xmath141.3  kpc , to avoid the strong contamination ( @xmath3050% ) of the nearest subsample .
@xcite estimated that the halo contamination in the sample should be about 8% , hence we would expect 41 halo members among the 514 targets with @xmath141.3  kpc .
the cut in @xmath20 removed 10 probable halo contaminants and , according to the statistics of @xcite , 18 of the excluded low - metallicity targets are expected to be halo stars . the cuts in @xmath7 and @xmath8 could also have removed some halo objects , but the effect of this cut is harder to quantify , because an unknown fraction of the outliers are probably objects with problematic proper motion or wrong distance . the residual halo contamination in our sample of 412 targets should therefore be lower than @xmath63% ( 13 stars ) , and it can be safely neglected .
the spectra were collected during six observing runs between 2005 and 2007 , at la silla and las campanas observatories .
the details of the observations and data reduction were presented in @xcite , and we will outline here only the most relevant points .
four different telescopes were used , because the stars span a wide range in magnitude ( @xmath8=516 ) , but the instruments had similar characteristics .
the basic information of the runs is given in table  [ t_obs ] .
the spectra covered the whole visual range from the atmospheric cutoff on the blue side to @xmath69000   , except for coralie data which only reached 6800    on the red edge .
we did not reduce the spectra collected with the mike blue arm because of too low stellar flux , and mike data were thus limited to 4850    on the blue end . the spectral resolution varied between 32000 and 50000 , depending on the spectrograph .
in particular , the 0@xmath317 slit was always used for mike runs ( r=32000 ) , while at dupont both the 0@xmath3175 and 1@xmath32 slits were used ( r=40000 and 30000 , respectively ) , depending on weather conditions . during each run
we secured the spectra of up to seventeen bright stars with accurate parameters from the literature ( radial and rotational velocities , temperature , gravity , metallicity ) .
exposure times were chosen according to weather conditions , to reach s / n=70100 for all the targets .
comparison lamp arcs were acquired only during daytime calibrations for the fiber spectrographs coralie and feros .
the second fiber of these two instruments was allocated to the sky background , because the spectra were usually collected next to full moon .
lamp arcs were collected at intervals of about 2 hours during dupont and mike runs , and each spectrum was calibrated with the average of the two lamps temporally enclosing it , linearly weighted by the temporal distance between science and calibration frames .
spectra were reduced with standard iraf[multiblock footnote omitted ] tasks , and we kept the procedures as much as possible identical for all the data .
we a posteriori verified that the reduced feros and coralie spectra were of the same quality as the products of their dedicated pipelines .
we analyzed the dark frames collected for all the runs , but we always found a negligible instrumental dark current and no clear 2d pattern , hence no dark correction was applied .
the wavelength calibration lamp spectra were extracted in the same position on the ccd as science targets , to avoid systematics introduced by the curvature of the orders .
the spectrum of a bright , well - exposed star was used to trace the orders on the ccd in each observing night , allowing only for rigid shifts among the frames .
then , the spectra were extracted with an optimum extraction algorithm @xcite , normalized , and finally merged .
the radial velocity ( rv ) for all the targets was measured with a cross - correlation ( cc ) technique @xcite as implemented in the iraf _ fxcor _ task .
the procedure was described in detail in @xcite , and we give only a brief summary here .
the spectrum of each star was cross - correlated with three template stars observed in the same run , namely hd180540 , hd223559 , and hd213893 ( this last replaced by hd224834 for coralie data ) , encompassing the temperature range of the targets .
their rv was fixed by a cc with a synthetic spectrum generated with the spectrum code @xcite , because we found poor agreement between the available literature sources .
the analysis of the seventeen standard stars acquired during observations , and of the solar spectra collected each night , revealed that the rv zero - point thus defined was biased by 0.3 - 0.7  km  s@xmath3 , depending on the instrument , and this offset was removed .
the three measurements were averaged , although they never differed by more than 0.2  km  s@xmath3 .
the cc was restricted to the range 4600 - 6800    ( 5000 - 6800    for mike data ) , to avoid the poor - quality blue section and the telluric bands at longer wavelengths .
all rvs were corrected to heliocentric velocities , then the position of the telluric band at 6875    was used to correct the rvs for instrumental effects , mainly caused by a displacement on ccd between the lamp and science spectra , and an off - center position of the star inside the slit ( see , for example , the analysis of * ? ? ?
corrections of up to 2  km  s@xmath3 were applied , but with little scatter ( @xmath60.5  km  s@xmath3 ) within each observing night .
the final rv errors were estimated as the quadratic sum of the five relevant uncertainties introduced in the procedure : the cc and wavelength calibration error , the uncertainty on the zero - point and its offset , and the error on the estimate of the correction for instrumental effects .
the resulting errors were typically in the range 0.4 - 0.7  km  s@xmath3 .
the final rvs of all the 824 stars will be made available at the cds website .
l c c reference & n. stars & @xmath33 + & & km s@xmath3 + @xcite & 145 & 0.8@xmath262.5 + @xcite & 9 & @xmath341.5 + @xcite & 9 & 0.4@xmath261.8 + @xcite & 8 & @xmath352.1 + @xcite & 6 & @xmath363.8 + @xcite & 6 & @xmath370.9 + @xcite & 6 & 1.8@xmath264.2 +    we found 211 previous rv measurements for the stars in our sample .
our results agree excellently with literature sources : the mean difference ( in the sense ours@xmath38literature ) is 0.4@xmath262.7  km  s@xmath3 , where the uncertainty is the rms of the differences . the comparison with the works that share with us more than five stars in common is given in table  [ t_lit ] .
the mean rv difference is always of the order of 1  km  s@xmath3 , except for @xcite , whose rvs are higher than ours by a non - negligible amount ( @xmath65  km  s@xmath3 ) .
however , the number of stars in common is too small to conclude that this offset is significant .
the proper motion , radial velocity , and distance of each target were combined to derive its ( @xmath39 ) velocity components in the galactic cylindrical reference frame , where @xmath7 is positive toward the galactic center , @xmath8 is directed in the sense of galactic rotation , and @xmath20 points toward the north galactic pole .
the error on these velocities was derived propagating the uncertainty on the observed quantities .
we assumed a proper motion error of 3  mas  yr@xmath3 for all the stars , as this value is more realistic than the uncertainties quoted in the spm3 catalog ( t. girard , 2009 , priv .
see also @xcite ) .
the velocities were corrected subtracting the solar peculiar motion ( @xmath40)=(10.0 , 5.1 , 7.2 )  km  s@xmath3 @xcite .
@xcite recently proposed the revised values ( @xmath40)=(11.0 , 12.2 , 7.3 )  km  s@xmath3 , but we preferred to adopt the older ones , for sake of continuity with previous works . in any case , the definition of the solar motion does not affect the velocity dispersions nor the off - diagonal terms of the dispersion matrix ( equation  [ e_cross ] ) , while the effects on the mean values are discussed in section  [ s_resmean ] .
the sample was then divided into several bins , in which we calculated the mean velocities , the dispersions ( @xmath41 ) , and the non - diagonal terms of the dispersion matrix @xmath42 where the sum is extended to all the stars in the bin , and @xmath43 .
the results are given in table  [ t_res ] .
the bins were defined by the requirement that their centers were equally spaced by 0.1  kpc from @xmath1=1.5 to 4.5  kpc , to uniformly sample the variation with @xmath1 of the kinematical quantities , while the width was imposed by the fixed number of stars per bin ( see below ) .
this implied a partial overlap of the bins , increasing with distance from the plane due to the decreasing density of observed stars .
however , @xcite have shown that the binning definition does not alter the results , and in fact our results are very similar to that work , despite the very different binning schemes .
the bins with @xmath442.4  kpc comprised 50  stars each , while at lower heights , where the number of observed stars is larger , the bin size was increased to 80  targets for 2.1@xmath452.4  kpc , and 100 targets for @xmath462.1  kpc .
we thus ensured a good statistic of thick disk stars in bins contaminated by the thin disk .    in each bin , the mean value and dispersion of each velocity component were measured by means of the analysis of the corresponding probability plot @xcite .
this powerful tool is highly insensitive to outliers , and it can be employed even in case of poorly populated bins .
the data are first ordered with increasing value of the velocity .
each point @xmath47 is then assigned a value @xmath48 , equal to the deviation from the mean expected for the @xmath49 point of a normal distribution , in units of the standard deviation . when the underlying distribution is gaussian
, the data follow a linear relation in the @xmath48-velocity plane , whose intercept and slope are given by the mean value and the dispersion , respectively .
the probability plots were fitted with a straight line , weighting each data point by its error , two examples are shown in figure  [ f_probplot ] .
the intrinsic velocity dispersions were then obtained quadratically subtracting the mean velocity error in the bin .
the wings of the distribution showing any deviation from linearity were excluded from the fit , suspected to be affected by objects with wrong distance or problematic proper motion , or by residual halo members .
the analysis of artificial data samples , generated through monte - carlo simulations and analyzed as the real data , indicated that the velocity dispersion is underestimated by 1 - 2  km  s@xmath3 when the thin disk contamination approaches 10% .
the contamination was therefore neglected for @xmath442.5  kpc , where it is expected to be lower than 5% ( i.e. two stars per bin , see figure  [ f_contam ] ) . following @xcite , in the presence of a non - negligible thin disk contamination ( @xmath462.5  kpc ) we derived the mean velocity and the dispersion of thick disk stars fitting only the wings of the probability plot , and only the negative one for @xmath8 .
the cut excluded the ranges @xmath17u@xmath1850  km  s@xmath3 , @xmath17w@xmath1830  km  s@xmath3 , and @xmath5060  km  s@xmath3 , approximatively 1.5  times the expected thin disk dispersion , thus removing about 90% of the contaminants .
the two wings of the probability plot should have the same slope but a different intercept , and the mean velocity was obtained averaging the two intercepts obtained from the fit . although the mean values thus derived agree well with the trends observed in the contamination - free bins at higher @xmath1 , we regard them as less reliable
.    the formal errors of the least - square procedure used to obtain the intercept and the slope in the probability plots mainly indicate the deviations from gaussianity of the observed distribution .
hence , they are not a good estimate of the real uncertainties .
the errors were therefore derived by means of monte - carlo simulations . for each bin
, we generated one thousand artificial samples of 50 stars , changing the proper motion , distance , and rv of each star assuming a gaussian random noise with dispersion equal to the estimated errors ( 3  mas  yr@xmath3 , 20% , and 0.4 - 0.7  km  s@xmath3 , respectively ) .
the measurements were repeated in each artificial bin as done for the real data , and the dispersion of these artificial measurements was assumed as an estimate of the observational errors .
the final errors are given in table  [ t_res ] .
the non - diagonal terms of the dispersion matrix , @xmath51 , @xmath52 , and @xmath53 , were calculated by means of equation  [ e_cross ] .
the errors were estimated from monte - carlo simulations , as done for the velocity dispersions . for each bin ,
one thousand artificial data sets of 50  values were simulated , drawn from a distribution characterized by the same dispersion matrix as the real data . in each simulation , a gaussian random noise with a dispersion equal to the observational errors was added , and the non - diagonal terms were estimated by means of equation  [ e_cross ] .
the errors were defined by the scatter of these simulated measurements . with the same simulations we also evaluated the systematic error introduced by the observational uncertainties ( e.g. * ? ? ?
* ) , comparing the results when the random noise was considered or not .
we found that the expected systematic thus introduced is negligible , being about one order of magnitude lower than the uncertainties on the measured values .
therefore , we did not correct the observed non - diagonal terms for observational errors .
c r r r r r r r r r @xmath1 & @xmath54 & @xmath55 & @xmath56 & @xmath4 & @xmath57 & @xmath5 & @xmath52 & @xmath51 & @xmath53 + kpc & km  s@xmath3 & km  s@xmath3 & km  s@xmath3 & km  s@xmath3 & km  s@xmath3 & km  s@xmath3 & km@xmath58  s@xmath59 & km@xmath58  s@xmath59 & km@xmath58  s@xmath59 + 1.5 & @xmath3831@xmath264 & @xmath3846@xmath264 & 17@xmath261 & 81@xmath265 & 57@xmath265 & 38@xmath261 & 570@xmath26310 & 400@xmath26470 & 10@xmath26210 + 1.6 & @xmath3820@xmath264 & @xmath3855@xmath264 & 5@xmath261 & 77@xmath265 & 59@xmath265 & 37@xmath261 & 840@xmath26290 & @xmath38180@xmath26460 & @xmath3820@xmath26220 + 1.7 & @xmath3828@xmath264 & @xmath3861@xmath265 & 3@xmath261 & 79@xmath265 & 64@xmath265 & 38@xmath261 & 930@xmath26310 & @xmath38130@xmath26510 & @xmath3840@xmath26250 + 1.8 & @xmath3824@xmath264 & @xmath3869@xmath265 & 8@xmath261 & 79@xmath265 & 61@xmath266 & 40@xmath261 & 810@xmath26330 & @xmath38760@xmath26500 & @xmath3880@xmath26240 + 1.9 & @xmath3816@xmath264 & @xmath3877@xmath265 & 3@xmath261 & 78@xmath265 & 60@xmath265 & 40@xmath261 & 850@xmath26330 & @xmath381060@xmath26500 & @xmath38120@xmath26240 + 2.0 & @xmath3810@xmath265 & @xmath3885@xmath264 & 4@xmath261 & 83@xmath266 & 55@xmath266 & 39@xmath261 & 940@xmath26340 & 30@xmath26460 & @xmath3890@xmath26210 + 2.1 & @xmath3822@xmath265 & @xmath3878@xmath265 & 3@xmath261 & 80@xmath266 & 58@xmath266 & 38@xmath261 & 1060@xmath26320 & 470@xmath26470 & @xmath38130@xmath26220 + 2.2 & @xmath3816@xmath265 & @xmath3876@xmath267 & @xmath3813@xmath261 & 77@xmath266 & 63@xmath267 & 42@xmath261 & 690@xmath26330 & 0@xmath26550 & 230@xmath26260 + 2.3 & @xmath381@xmath266 & @xmath3881@xmath265 & @xmath389@xmath261 & 81@xmath267 & 58@xmath266 & 40@xmath261 & 750@xmath26340 & @xmath38500@xmath26470 & 10@xmath26230 + 2.4 & @xmath3812@xmath265 & @xmath3885@xmath266 & @xmath3812@xmath261 & 80@xmath266 & 59@xmath266 & 40@xmath261 & 410@xmath26330 & @xmath38710@xmath26480 & 180@xmath26240 + 2.5 & @xmath3818@xmath266 & @xmath3885@xmath266 & @xmath382@xmath261 & 78@xmath267 & 63@xmath267 & 42@xmath261 & 630@xmath26330 & @xmath381870@xmath26600 & @xmath3840@xmath26270 + 2.6 & @xmath3825@xmath266 & @xmath3890@xmath266 & 1@xmath261 & 71@xmath267 & 66@xmath267 & 42@xmath261 & 530@xmath26310 & @xmath382390@xmath26660 & @xmath38460@xmath26290 + 2.7 & @xmath3823@xmath266 & @xmath3898@xmath266 & 7@xmath261 & 72@xmath267 & 62@xmath267 & 39@xmath261 & 870@xmath26290 & @xmath381150@xmath26520 & @xmath381360@xmath26300 + 2.8 & @xmath3829@xmath266 & @xmath3895@xmath267 & @xmath382@xmath261 & 76@xmath267 & 62@xmath267 & 41@xmath261 & 880@xmath26320 & @xmath381740@xmath26580 & @xmath381310@xmath26300 + 2.9 & @xmath3817@xmath266 & @xmath38115@xmath267 & 5@xmath261 & 83@xmath267 & 68@xmath267 & 40@xmath261 & 840@xmath26340 & @xmath38820@xmath26590 & @xmath381230@xmath26300 + 3.0 & @xmath387@xmath266 & @xmath38126@xmath267 & 8@xmath261 & 90@xmath267 & 67@xmath268 & 42@xmath261 & 1440@xmath26410 & @xmath3850@xmath26610 & @xmath38460@xmath26280 + 3.1 & @xmath386@xmath267 & @xmath38129@xmath268 & 1@xmath261 & 101@xmath268 & 67@xmath268 & 43@xmath261 & 1610@xmath26450 & @xmath381150@xmath26690 & @xmath38650@xmath26290 + 3.2 & 2@xmath266 & @xmath38131@xmath268 & @xmath381@xmath261 & 99@xmath268 & 63@xmath268 & 42@xmath261 & 1470@xmath26440 & @xmath381820@xmath26650 & @xmath38440@xmath26280 + 3.3 & 8@xmath267 & @xmath38140@xmath268 & 0@xmath261 & 101@xmath268 & 66@xmath268 & 44@xmath262 & 2180@xmath26490 & @xmath382610@xmath26740 & @xmath38310@xmath26300 + 3.4 & 12@xmath267 & @xmath38140@xmath268 & 1@xmath261 & 98@xmath268 & 63@xmath269 & 43@xmath262 & 2100@xmath26480 & @xmath381880@xmath26660 & 260@xmath26280 + 3.5 & 18@xmath267 & @xmath38135@xmath269 & 5@xmath261 & 95@xmath269 & 63@xmath269 & 44@xmath262 & 1950@xmath26460 & @xmath383390@xmath26740 & 260@xmath26280 + 3.6 & 29@xmath267 & @xmath38137@xmath268 & 9@xmath261 & 101@xmath269 & 64@xmath269 & 44@xmath262 & 2290@xmath26510 & @xmath382170@xmath26700 & 360@xmath26290 + 3.7 & 7@xmath267 & @xmath38132@xmath269 & 0@xmath262 & 91@xmath269 & 61@xmath2610 & 44@xmath262 & 2180@xmath26460 & @xmath382150@xmath26610 & 610@xmath26280 + 3.8 & 15@xmath267 & @xmath38134@xmath2610 & 1@xmath262 & 92@xmath269 & 68@xmath2610 & 43@xmath262 & 1990@xmath26450 & @xmath382120@xmath26700 & @xmath3830@xmath26300 + 3.9 & 5@xmath268 & @xmath38139@xmath269 & @xmath386@xmath262 & 94@xmath2610 & 66@xmath2610 & 43@xmath262 & 1540@xmath26430 & @xmath38520@xmath26620 & 210@xmath26280 + 4.0 & 7@xmath268 & @xmath38143@xmath2610 & @xmath381@xmath262 & 93@xmath269 & 66@xmath2611 & 41@xmath262 & 1910@xmath26420 & @xmath38100@xmath26610 & 10@xmath26280 + 4.1 & 17@xmath268 & @xmath38145@xmath2610 & 1@xmath262 & 92@xmath2610 & 68@xmath2610 & 44@xmath262 & 1570@xmath26440 & 720@xmath26640 & @xmath3830@xmath26300 + 4.2 & 22@xmath268 & @xmath38142@xmath2611 & 2@xmath262 & 95@xmath2610 & 72@xmath2611 & 45@xmath262 & 1480@xmath26460 & 880@xmath26710 & 10@xmath26320 + 4.3 & 12@xmath269 & @xmath38148@xmath2610 & @xmath382@xmath262 & 94@xmath2610 & 69@xmath2611 & 48@xmath262 & 2320@xmath26510 & 1170@xmath26680 & 180@xmath26330 + 4.4 & 17@xmath269 & @xmath38152@xmath2611 & 0@xmath262 & 96@xmath2610 & 72@xmath2611 & 46@xmath262 & 2330@xmath26500 & 1370@xmath26710 & @xmath3850@xmath26330 + 4.5 & 19@xmath2610 & @xmath38158@xmath2611 & @xmath382@xmath262 & 93@xmath2611 & 76@xmath2611 & 46@xmath262 & 2180@xmath26480 & 3960@xmath26930 & 240@xmath26360 +
the trend of the mean velocities with distance from the plane is shown in figure  [ f_resmean ] .
@xmath56 is consistent with no significant departure from zero in the whole range .
some deviations are observed for @xmath602.5  kpc , but in this range the results are affected by large errors and , as already commented , not very trustworthy .
this result agrees with previous studies , that found no net vertical motion up to various kpc from the galactic plane @xcite . on the contrary
, @xmath61 has a more complex behavior : between @xmath1=1.5 and 3  kpc we detect a non - null mean velocity directed toward the galactic anticenter , with an average value @xmath54=@xmath623  km  s@xmath3 . beyond 3  kpc from the plane
@xmath54 abruptly increases and changes sign , and the net velocity between 3.5 and 4.5  kpc is 15@xmath262  km  s@xmath3 toward the galactic center .
another possible interpretation of the results is that @xmath54 linearly increases with @xmath1 , and a s - shaped feature is superimposed to this trend between 2.5 and 3.5  kpc , as observed for @xmath4 and other kinematical quantities ( see section  [ s_comove ] ) . after the exclusion of this range , the fit returns a vertical increase of 15@xmath63  km  s@xmath3  kpc@xmath3 .
adopting the new values of @xcite for the solar peculiar motion causes not a relevant change on the results , as @xmath54 would be higher by about 1  km  s@xmath3 only .
previous studies have shown evidence that stars more distant than 1  kpc from the plane show a non - null net radial motion of the order of @xmath69  km  s@xmath3 toward the galactic anticenter @xcite .
our results between @xmath1=1.5 and 3  kpc agree with their conclusion , although our mean value is larger by about a factor of two .
the sign flip observed at larger @xmath1 , on the contrary , has never been reported in the literature .
nevertheless , @xcite analyze halo stars , @xcite do not reach these heights above the plane , and @xcite have too few stars in this range . very interestingly ,
@xcite detect a small positive mean radial component out to @xmath64  kpc , and a negative value in their more distant bins , although the authors conclude that these results are consistent , within errors , with a null net motion .
@xcite proposed that an inward motion of the local standard of rest ( lsr ) of the order of 10  km  s@xmath3 is responsible for the observed non - null mean value of @xmath54 , as suggested by the fact that this is observed among both disk and halo stars more distant than 1  kpc from the sun .
this hypothesis can not account for any vertical trend of @xmath54 other than a constant non - zero value at any @xmath1 .
thus , while it is not contradicted by our results , at least another effect must be invoked .
for example , an outward motion of @xmath625  km  s@xmath3 of the stars between 3.5 and 4.5  kpc , coupled with the mentioned lsr motion , could explain our results .
this clumpy kinematical distribution would not be surprising , because perturbations produced by the bar and the spiral arms are expected to form groups of stars with velocity offset as large as 50  km  s@xmath3 @xcite .
@xcite studied the proper motions of our same stellar sample , and their results for @xmath61 are overplotted to ours in the upper panel of figure  [ f_resmean ] , after correcting for the solar peculiar motion and changing the sign of @xmath7 to account for the different direction of the reference axis .
the agreement with our results is poor : while their nearly flat profile roughly coincides with our mean value ( @xmath387.7  km  s@xmath3 ) , they did not detect any change of sign , nor a steep positive gradient .
the different approach to the same data must have caused the different results .
for example the features observed by us could have been masked in @xcite by their smoothed , statistical distance determination , or by dwarf stars and halo contaminants , removed in our work .    in the middle panel of figure  [ f_resmean ]
we show the measured trend of @xmath65 , compared to other results from the literature .
the mean rotational velocity of thick disk stars decreases with distance from the plane , due to the well - known vertical shear first detected by @xcite .
the fit of the data , after the exclusion of the less reliable bins at @xmath602.5  kpc , yields @xmath66  ( km  s@xmath3 ) .
as can be seen in figure  [ f_resmean ] , this solution excellently matches the results of @xcite , and even the data points of @xcite , at @xmath67  kpc , are well aligned with the derived linear relation .
both these investigations measure a vertical shear of @xmath3830  km  s@xmath3  kpc@xmath3 , as also recently found by @xcite , and @xcite .
had we fitted all the data points down to @xmath1=1.5  kpc , we would have found a steeper slope ( @xmath68  km  s@xmath3  kpc@xmath3 ) , at the upper limit of the range spanned by literature values , which vary from @xmath694 @xcite to @xmath701  km  s@xmath3 @xcite . the revised values for the solar peculiar motion proposed by @xcite cause an upward revision of the results by 7.1  km  s@xmath3 .
@xcite and @xcite measured a shallower slope ( @xmath712.1 and @xmath622  km  s@xmath3  kpc@xmath3 , respectively ) between 0.7 and 2.8  kpc .
while our results are compatible with theirs in our nearest bins , we find a lower mean rotational velocity beyond @xmath1=2  kpc .
nevertheless , the different vertical rotational gradient found by these studies is not necessarily a disagreement , because they sample a different @xmath1-range , and the underlying shear is not required to be strictly linear .
in fact , @xcite have proposed the non - linear relation @xmath72  km  s@xmath3 from the analysis of sdss data .
the combined data points of @xcite , @xcite , @xcite , and of the present work , closely follow this equation ( see figure  [ f_vall ] ) , and a fit in the form @xmath73  km  s@xmath3 returns a very similar solution , with @xmath74  km  s@xmath3 , @xmath75  km  s@xmath3  kpc@xmath3 , and @xmath761.23 .
very noticeably , the results of four surveys finding a different linear relation are all well reproduced by a single non - linear curve proposed by an independent work .
the underlying vertical trend of the thick disk rotational velocity is therefore most probably close to but not exactly linear , and better represented by a low - order power law .
the vertical profile of the velocity dispersions is shown in figure  [ f_disp ] , compared to other works from the literature .
we excluded from comparison the results of @xcite , because they do not distinguish between different disk components , thus finding steeper gradients as a result of the mix of thin and thick disk stars .
caution must also be taken when comparing the results obtained in different ranges of @xmath1 , because the underlying trend of the dispersions is not necessarily linear , and the gradient can assume different values .
for example , the models of @xcite indicate that the vertical profile should be progressively flatter at increasing @xmath1 .
the linear fit of the trends with @xmath1 yields : @xmath77 @xmath78 @xmath79 where @xmath1 is in kpc .
the quoted errors were obtained restricting the fit to a subset of eight non - overlapping bins , to avoid the underestimate induced by the correlation between the data points .
the results are identical to those of @xcite , despite the different binning scheme , except for @xmath57 for which we derive a slightly smaller gradient , but the difference is not significant ( 0.4  km  s@xmath3  kpc@xmath3 ) .
the gradients are small , and the three dispersions increase by only @xmath67% between 2.5 and 3.5  kpc .
this explains why the change of the thick disk kinematics with distance from the plane has not been clearly detected for about two decades after its discovery .
the results for @xmath5 are the most precise and reliable , because @xmath690% of the information on @xmath20 comes from rvs , whose errors are an order of magnitude smaller than those in proper motions . @xmath80 shows a small but clear and constant increase , with small scatter around the derived linear relation . the vertical gradient is smaller than the one found by @xcite by about a factor of two , but the results are consistent at the 1@xmath48 level .
@xcite propose an even steeper profile ( @xmath610  km  s@xmath3  kpc@xmath3 ) in the range @xmath1=02  kpc , but the same authors suspect that this could be due to contamination by halo stars increasing with galactic height .
the dispersion of the rotational velocity component , @xmath57 , also shows a clear increase with @xmath1 , although the data points are affected by larger errors and are more scattered . the vertical gradient is smaller than the results of both @xcite and @xcite , by a factor of three and two , respectively , and the difference is at the 2@xmath48 level in both cases .
the vertical profile of @xmath4 shows a very peculiar behavior , with large deviations from linearity between 2.5 and 3.5  kpc .
this feature will be discussed in section  [ s_comove ] . outside this range ,
the data points closely follow a linear relation , and the solution given in equation  [ e_su ] was obtained excluding this interval from the fit .
the derived vertical gradient of @xmath4 is very similar to that found by @xcite , but their solution is offset toward higher values by about 10  km  s@xmath3 , while the gradient measured by @xcite is 4@xmath48 times higher than ours .    in conclusion , we confirm that the velocity dispersions of the galactic thick disk increase with distance from the plane , as suggested by previous investigations @xcite , but we derive vertical gradients that are in general smaller than other studies
. the differences with @xcite can be at least in part due to the aforementioned expected decrease of the gradient with distance from the plane , because the results for @xmath5 are consistent , and the solutions proposed for the other two components overlap in the @xmath0-range in common ( 1.52.2  kpc ) .
on the contrary , @xcite studied our same sample , and the different data analysis must be the cause of the disagreement .
since they measured both higher dispersions and steeper vertical gradients than us , their results could have been affected by the dwarf stars and halo contaminators removed by us , and/or the thin disk contamination , that they did not take into account .
the three dispersions increase with distance from the plane proportionally , and the anisotropy is approximatively constant at all @xmath1 . in figure  [ f_ratio ]
the vertical profile of the ratios @xmath4/@xmath5 and @xmath4/@xmath57 is shown , where the errors come from propagation of the uncertainties on the dispersions .
the mean value of @xmath4/@xmath5 is 2.08@xmath260.06 , where the error was estimated from the statistical error on the mean for eight uncorrelated bins , as before .
the linear fit of the data points returns a negligible gradient ( 0.09@xmath260.07  kpc@xmath3 ) , and the data are consistent with a flat profile .
literature values for @xmath4/@xmath5 span a wide range from 1.1 @xcite to 1.9 ( e.g. , * ? ? ?
* ) , and our results are at the upper end of this distribution . as shown in the upper panel of figure  [ f_ratio ] , the results of @xcite agree with our measurements in the z - range in common , but they deviate from our measurements if extrapolated to higher distance from the plane
.    the results for @xmath4/@xmath57 also show no significant gradient ( 0.06@xmath260.05  kpc@xmath3 ) and a mean value of 1.36@xmath260.05 .
this is lower than the value predicted by the epicyclic approximation ( 1.49 , cfr .
equation 3 - 76 of * ? ? ?
* ) , indicating that the population under analysis can not be assumed in nearly circular orbits .
our results are well within the range spanned by literature values , which vary from @xmath811 ( e.g. , * ? ? ?
* ; * ? ? ?
* ) , to @xmath811.6 @xcite .
as in the case of @xmath4/@xmath5 , the results of @xcite overlap with our data points in the range in common , but their extrapolation to higher @xmath1 do not .
the analytical expressions given by @xcite for @xmath4 and @xmath57 , on the contrary , return a value of their ratio much lower than our data points at any galactic height .
this is probably due to their much higher @xmath57 , as shown in the middle panel of figure  [ f_disp ] .
the non - diagonal term of the dispersion matrix @xmath82 gives a measurement of the correlation between the @xmath47-th and @xmath83-th velocity components .
in fact , if two principal axis of the dispersion ellipsoid are aligned with the @xmath47-th and @xmath83-th axis , we have @xmath82=0 , and the two velocity components are uncorrelated .
the orientation of the dispersion ellipsoid of an old , dynamically relaxed population is related to the shape of the galactic potential @xcite , and can indicate the presence of non - axisymmetric structures in the disk @xcite .
moreover , @xcite and @xcite have shown that the vertical tilt of the ellipsoid allows to estimate the flattening of the dark halo . however , calculating the expectation of a specific halo model through integration of orbits is beyond the scope of the paper , and here we will only present the observational results , for use of future works .    our results for the three cross - terms are shown in figure  [ f_cross ] .
the profile of @xmath51 closely follows a decreasing linear relation up to @xmath1=4  kpc , but abruptly jumps to positive values at larger heights .
this is probably due to the sign flip of @xmath54 analyzed in section  [ s_comove ] . in any case , the results indicate that the @xmath7 and @xmath8 velocity components are correlated , and the velocity ellipsoid is titled in the radial - longitudinal plane .
the rotation angle , i.e. the vertex deviation @xmath84 , was calculated in each bin by means of the relation : @xmath85 @xcite .
the results are shown in the upper panel of figure  [ f_angles ] .
we measure a non - negligible vertex deviation , increasing from nearly zero at @xmath1=1.5  kpc to @xmath86 at 3.5  kpc .
the linear fit in the range 1.54  kpc yields the relation @xmath87 .
our results agree with @xcite , who measured @xmath88 at @xmath0=1.1  kpc .
previous investigations evidenced that the vertex deviation decreases from @xmath89 for young populations to nearly zero for old , metal - poor disk stars @xcite .
our results are not at variance with this conclusion if the vertex deviation increases with @xmath1 , as suggested by our observations , because previous studies were limited to small galactic heights .
@xmath53 shows no significant deviation from zero in the whole range of @xmath1 , except between 2.5 and 3.5  kpc , as discussed in section  [ s_comove ] .
this indicates that the @xmath8 and @xmath20 velocity components are not correlated . on the contrary , @xmath52 is significantly different from zero , and steadily increases with @xmath1 . proposing a reliable linear expression for @xmath90 is not straightforward , because it shows irregular fluctuations . as expected by @xcite and @xcite , @xmath52
is always bound between zero and @xmath91 , the value assumed when the velocity ellipsoid is aligned with the spherical coordinate system , because we find that @xmath92  km@xmath93  s@xmath59 at any galactic height .
we also find that the expression proposed by @xcite , obtained under the assumption that the ellipsoid points toward the galactic center and its axis ratio is constant in the spherical coordinate system , is a relatively good approximation of the measured value up to about 3  kpc from the plane , but the predicted slope is too shallow and the agreement with observations degrades with @xmath1 ( see figure  [ f_cross ] ) .
the tilt angle in the @xmath7-@xmath20 plane can be calculated , analogously to the vertex deviation , from the equation : @xmath94 and the results are shown in figure  [ f_angles ] .
we measure a slight increase of @xmath95 with @xmath1 , and the linear fit of the data , after the exclusion of two deviating points , yields the relation @xmath96 $ ] .
the results agree with previous works at @xmath1=11.5  kpc @xcite , but not with the measurements of @xcite at @xmath971  kpc , whose sample is most probably dominated by thin disk stars .
the dashed curve in figure  [ f_angles ] indicates the value of @xmath95 when the dispersion ellipsoid points toward the galactic center , assuming r@xmath98=8  kpc .
the tilt angle is constantly lower than this in the range 1.54.5  kpc , and the ellipsoid thus is directed toward a point located behind the galactic center , at a galactocentric distance increasing with @xmath1 from r@xmath99=2.3 to 9  kpc .
orbit integration studies indeed predict this result @xcite , but the inclination is noticeably higher than the expectations , because theoretical calculations return r@xmath99=510  kpc at @xmath1=1.1  kpc .
this could indicate that the galactic potential used in these studies needs to be refined to match the observations . on the contrary ,
the increase of @xmath100 modeled by @xcite is too steep , and their expectation @xmath95(3.5  kpc)=26@xmath101 is at variance with our results .
@xcite showed that the vertical trend of the tilt angle is an excellent observational signature of the underlying gravity law , and it can be used as a test for the modified newtonian dynamics theory ( mond , * ? ? ?
* ) , because its expectation diverges from that of the newtonian dynamics with distance from the plane .
unfortunately , the calculations of @xcite are limited to @xmath972  kpc , and the range of overlap with our data is very narrow .
the value of @xmath95 expected by @xcite at @xmath1=2  kpc is 12@xmath101 for newtonian dynamics and
10@xmath101 for mond .
we find @xmath95(2  kpc)=@xmath102 and the mean of the five measurements in the range 1.8 - 2.2  kpc is @xmath103 , in excellent agreement with mond expectations .
the linear fit presented in section  [ s_resangles ] also implies @xmath95(2  kpc)=@xmath104 .
however , we derive a vertical gradient ( @xmath105  kpc@xmath3 ) much shallower than the predictions of @xcite for both mond ( 5@xmath101  kpc@xmath3 ) and newtonian dynamics ( 6@xmath101  kpc@xmath3 ) . in conclusion ,
our observations agree better with the models that @xcite derived for mond than with those from newtonian dynamics , but the extension of @xcite calculations to higher @xmath1 , and to specific initial conditions for the thick disk , are required to perform a reliable test of the gravitational law .
the vertical profile of @xmath106 shows a puzzling behavior between 2.5 and 3.5  kpc from the galactic plane ( top panel , figure  [ f_disp ] ) , where the dispersion deviates from the linear trend .
a similar feature is clearly observable in the trend of other quantities related to the radial velocity component @xmath7 , like @xmath54 and @xmath52 ( as already noted by * ? ? ?
* ) , but deviations from the linear trend in the range @xmath1=2.53  kpc could be present even in other profiles , such as @xmath65 , @xmath107 , @xmath80 , and @xmath108 .
the origin of this behavior is unclear .
it is possible that a group of comoving stars , forming a sub - structure in the galactic thick disk , is affecting the measured kinematics between 2.5 and 3.5  kpc from the plane . in particular , the similarity of the profile of two totally unrelated quantity such as @xmath108 and @xmath106 is instructive , because it excludes the possibility that this behavior is only due to some bad measurements . in any case , the existence of a kinematical substructure among our stars can not be claimed on the basis of these results only .
our targets are distributed in a very narrow range of @xmath109 , and the data do not provide any direct information about the radial behavior of the kinematical quantities .
nevertheless , some indirect indication can be derived . manipulating the radial jeans equation in cylindrical coordinates in steady state , with the radial component of the force expressed as f@xmath110
, we obtain : @xmath111 where we also assumed that @xmath112 , @xmath113/@xmath114=0 , and that the density decays exponentially with both @xmath109 and @xmath0 , with scale length @xmath115 and scale height @xmath116 . inserting the observed quantities in the right hand side of equation  [ e_jeans ]
we can thus estimate , as an exercise , the radial behavior of @xmath4 .
the large uncertainties involved prevent a precise measurement , but the comparison with the expectations provides a consistency check for the kinematical data here presented , because an unphysical result could indicate a problem with them .
the results are shown in figure  [ f_jeans ] , where we assumed @xmath115=3.6  kpc and @xmath116=0.9  kpc @xcite , r@xmath98=8  kpc , and @xmath117=220  km  s@xmath3 .
the data points are scattered around the mean value @xmath118  km  s@xmath3  kpc@xmath3 , following the behavior of @xmath106 discussed in section  [ s_comove ] , with no clear vertical trend . in the nearest bins the results are partially consistent , but more negative , than what found by @xcite and ( * ? ? ?
* @xmath119  km  s@xmath3  kpc@xmath3 ) , although neese & yoss refer to the thin disk .
the observations show @xcite that @xmath5 should exponentially decay in the radial direction , with the same scale length of the mass density .
this is often considered valid even for @xmath4 , under the assumption of a constant anisotropy , and @xcite found that this should be the best approximation at the solar position .
the results also roughly agree with the expectation of this model , as shown in figure  [ f_jeans ] . in conclusion ,
the presented kinematics implies a radial behavior of @xmath4 consistent with both previous observations and with the expectations of theoretical predictions .      the presence of kinematical gradients in the galactic thick disk is a powerful diagnostic to discriminate between the various models of its formation @xcite , but unfortunately the model expectations concerning the vertical gradient have not been investigated so far
. it would be interesting to compare , in the near future , the predictions of the various models with the observed trend with @xmath0 of the velocity dispersions .
on the contrary , the radial gradient of the dispersions has been modeled in the context of thick disk formation through disk heating by the merging of minor satellites @xcite .
the observations of @xcite agree with the results of these simulations , provided that a low inclination orbit of the merging satellite is assumed .
a similar conclusion is drawn comparing our measurements of @xmath4/@xmath5 with the theoretical expectations of the merging scenario .
in fact , the simulations showed that this ratio is strongly linked to the inclination angle of the orbit @xcite , and our result ( @xmath4/@xmath5=2.08 ) favors a small inclination angle , @xmath120030@xmath101 ( see figure  15 of * ? ? ? * ) .
a low inclination orbit of the infalling satellite is also required to reproduce the rather large vertical gradient of the rotational velocity , found in this work as in previous investigations @xcite . on the contrary ,
the radial migration model predicts a much shallower gradient ( @xmath3817  km  s@xmath3  kpc@xmath3 , * ? ? ?
* ) , incompatible with the observations .    in conclusion
, the models of thick disk formation have so far provided only fragmentary predictions about its kinematical properties , and our observations can not still be used to fully discriminate between them .
however , we find that our results are consistent with the scenario where the thick disk formed through dynamical heating of a pre - existing galactic disk , induced by the merging of a minor satellite .
moreover , all the kinematical evidence shows that , if this is the correct model , a low - latitude ( @xmath121 ) merging event is strongly preferred .
we have analyzed a sample of @xmath6400 thick disk stars , measuring the variation of their kinematical properties as a function of distance from the galactic plane , from @xmath1=1.5 to 4.5  kpc .
our results can be summarized as follows :    * while the mean vertical velocity component @xmath56 shows no significant deviation from zero in the whole range , between 1.5 and 3  kpc , we find a net radial motion of about 20  km  s@xmath3 directed toward the galactic anticenter .
other authors have recently found evidence for a similar behavior , proposing an inward motion of the lsr , although our mean velocity is larger than their proposed value by a factor of two .
however , we find that @xmath54 changes sign for @xmath1223  kpc , and that a radial motion of the lsr , although not excluded , can not alone explain this behavior .
* the mean rotational velocity of the thick disk decreases with distance from the galactic plane , as found by many previous investigations .
the linear fit of our data returns a gradient of @xmath2  km  s@xmath3  kpc@xmath3 , our data points and the results of three other previous works are globally better represented by a power - law with index 1.25 , very similar to what has recently been proposed by @xcite . *
all the velocity dispersions steadily increase with distance from the galactic plane , closely following a linear relation .
the gradients we found are , however , smaller than those proposed by previous works .
* while the velocity dispersions increase with @xmath1 , the ratios @xmath123 and @xmath124 show no significant vertical trend . the observations thus indicate a substantial constancy with @xmath1 of the anisotropy . *
we find a non - negligible vertex deviation , increasing with @xmath1 from values close to zero to @xmath86 at @xmath1=3.5  kpc .
this is consistent with previous investigations , which found a very small vertex deviation of old stellar population close to the galactic plane .
* the tilt angle steadily increases with distance from the galactic plane .
as expected , the orientation of the velocity ellipsoid in the @xmath7-@xmath20 plane results , at any @xmath1 , intermediate between alignment with the cylindrical and spherical coordinate systems . according to calculations by @xcite , the tilt angle at @xmath1=2  kpc coincides with the expectation of mond , although the extension of their models to higher galactic heights is required to perform a conclusive test of the underlying gravitational law . *
the vertical trend of many kinematical quantities show deviations from linearity between 2.5 and 3.5  kpc .
the origin of these features is unknown , but it could indicate the presence of a sub - structure at this galactic height , such as a comoving group of stars . *
the results are fully consistent with the model of thick disk formation through dynamical heating of a pre - existing galactic disk . if this is the correct scenario , a low inclination angle of the merging event is strongly preferred .
however , not all the models proposed so far could be tested by our observations , and more simulations are required to obtain a detailed comparison able to discriminate between them .
c.m.b . and r.a.m .
acknowledge support from the chilean centro de astrofsica fondap no . 15010003 , and the chilean centro de excelencia en astrofsica y tecnologas afines ( cata ) basal pfb/06 . c.m.b .
also thanks f. mauro for useful discussions .
all authors acknowledge partial support from the yale university / universidad de chile collaboration . the spm3 catalog was funded in part by grants from the us national science foundation , yale university and the universidad nacional de san juan , argentina .
we warmly thank w.  f. van altena , v.  i. korchagin , t.  m. girard , and d.  i. casetti - dinescu for their help and suggestions .
allende prieto , c. , beers , t.  c. , wilhelm , r. , et al .
2006 , apj , 636 , 804 amendt , p. , & cuddeford , p. 1991 , apj , 368 , 79 assmann , p. , fellhauer , m. , kroupa , p. , brns , r.  c. , & smith , r. 2011 , mnras , 415 , 1280 beers , t. c. , & sommer - larsen , j. 1995 , apjs , 96 , 175 bekki , k. , & tsujimoto , t. 2011 , apj , 738 , 4 bensby , t. , feltzing , s. , & lundstrm , i. 2003 , a&a , 410 , 527 bensby , t. , zenn , a. r. , oey , m. s. , & feltzing , s. , 2007 , apj , 663 , l13 bienaym , o. 1999 , a&a , 341 , 86 bienaym , o. 2009 , a&a , 500 , 781 bienaym , o. , famaey , b. , wu , x. , zhao , h.  s. , & aubert , d. 2009 , a&a , 500 , 801 binney , j. , & merrifield , m. 1998 , galactic astronomy ( princeton , nj : princeton university press ) , 1998 binney , j. , & spergel , d. 1983 , the nearby stars and the stellar luminosity , ed .
philip , & u. schenectady ( n.y .
: davis press ) , iau col .
, 76 , 259 bochanski , j.  j. , munn , j.  a. , hawley , s.  l. , et al .
2007 , aj , 134 , 2418 bond , n.  a. , ivezi , z. , sesar , b. , et  al .
2010 , apj , 716 , 1 bournaud , f. , elmegreen , b.  g. , & martig , m. 2009 , apj , 707 , l1 brown , w.  r. , beers , t.  c. , wilhelm , r. , et  al .
2008 , aj , 135 , 564 cardelli , j. a. , geoffrey , c. c. , & mathis , j. s. 1989 , apj , 345 , 245 carney , b.  w. , latham , d.  w. , & laird , j.  b. 1989 , aj , 97 , 423 carollo , d. , beers , t.  c. , chiba , m. , et  al .
2010 , apj , 712 , 692 carraro , g. , van altena , w.  f. , moni bidin , c. , et  al .
2005 , aas , 37 , 1378 casetti - dinescu , d.  i. , girard , t.  m. , korchagin , v.  i. , & van altena , w.  f. 2011 , apj , 728 , 7 chiba , m. , & beers , t. c. 2000 , aj , 119 , 2843 cuddeford , p. , & amendt , p. 1992 , mnras , 256 , 166 dalcanton , j.  j. , & bernstein , r.  a. 2002 , aj , 124 , 1328 dehnen , w. , & binney , j. 1998 , mnras , 298 , 387 duflot , m. , figon , p. , & meyssonnier , n. 1995 , a&as , 114 , 269 evans , l.  t. 1970 , mnras , 149 , 311 feltzing , s. , bensby , t. , & lundstrm , i. 2003 , a&a , 397 , l1 fuchs , b. , dettbarn , c. , rix , h .- w .
, et al . 2009 , aj , 137 , 4149 fuhrmann , k. 1998 , a&a , 338 , 161 flynn , c. , & freeman , k.  c. 1993 , a&as , 97 , 835 gilmore , g. , & reid , n. 1983 , mnras , 202 , 1025 girard , t.  m. , dinescu , d.  i. , van altena , w.  f. , et al .
2004 , aj , 127 , 3060 girard , t.  m. , korchagin , v.  i. , casetti - dinescu , d.  i. , et al .
2006 aj , 132 , 1768 gmez , f.  a. , minchev , i. , villalobos , a. , oshea , b.  w. , & williams , m.  e.  k. 2011 ( arxiv:1105.4231 ) gratton , r.  g. , bragaglia , a. , carretta , e. , et al .
2003 , a&a , 408 , 529 gratton , r.  g. , carretta , e. , matteucci , f. , & sneden , c. 2001 , a&a , 358 , 671 gray , r.  o. , & corbally , c.  j. 1994 , aj , 107 , 742 hamaker , h.  c. 1978 , appl .
stat . , 27 , 7 haywood , m. 2001 , mnras , 325 , 1365 hori , l. , & liu , t. 1963 , pasj , 15 , 100 horne , k. 1986 , pasp , 98 , 609 ivezi , z. , sesar , b. , juri , m. , et  al .
2008 , apj , 684 , 287 jones , d.  h.  p. 1972 ,
apj , 178 , 467 juri , m. , ivezi , z. , brooks , a. , et  al .
2008 , apj , 673 , 864 kent , s. , & de zeew , t. 1991 , aj , 102 , 1994 kobayashi , c. , & nakasato , n. 2011 , apj , 729 , 16 kuijken , k. , & gilmore , g. 1989 , mnras , 239 , 571 kuijken , k. , & tremaine , s. 1991 , in dynamics of disc galaxies , varberg castle , sweden , 71 loebman , s.  r. , roskar , r. , debattista , v.  p. , et  al .
2011 , apj , 737 , 8 lutz , t. e. , & hanson , r. b. 1992 , asp conf .
series , 25 , 257 lynden - bell , d. 1962 , mnras , 124 , 95 majewski , s.  r. 1992 , apjs , 78 , 87 majewski , s.  r. 1993 , ara&a , 31 , 575 majewski , s.  r. , kunkel , w.  e. , law , d.  r. , et  al .
2004 , aj , 128 , 245 milgrom , m. 1983 , apj , 270 , 371 monaco , l. , villanova , s. , moni bidin , c. , et  al .
2011 , a&a , 529 , a90 moni bidin , c. , moehler , s. , piotto , g. , et  al .
2006 , a&a , 451 , 499 moni bidin , c. , girard , t.  m. , carraro , g. , et  al .
2009 , rmxac , 35 , 109 moni bidin , c. 2009 , phd thesis , univ . de chile , 1 moni bidin , c. , carraro , g. , mndez , r.  a. , & van altena , w.  f. 2010 , apj , 724 , l122 neese , c.  l. , & yoss , k.  m. 1988 , aj , 95 , 463 ollongren , a. 1962 , baas , 16 , 561 purcell , c.  w. , bullock , j.  s. , & kaplinghat , m. 2009 , apj , 703 , 2275 quillen , a.  c. , doughery , j. , bagley , m.  b. , minchev , i. , & comparetta , j. 2010 , mnras , submitted ( arxiv:1010.5745 ) quinn , p.  j. , hernquist , l. , & fullagar , d.  p. 1993 , apj , 403 , 74 reddy , b.  e. , lambert , d.  l. , & allende prieto , c. 2006 , mnras , 367 , 1329 rocha - pinto , h.  j. , maciel , w.  j. , scalo , j. & flynn , c. 2000 , a&a , 358 , 850 roskar , r. , debattista , v.  p. , stinson , g.  s. , et  al .
2008 , apj , 675 , l65 rybka , s.  p. , & yatsenko , a.  i. 2010 , kpcb , 25 , 309 schlegel , d.  j. , finkbeiner , d.  p. , & davis , m. 1998 , apj , 500 , 525 schnrich , r. , & binney , j. 2009 , mnras , 396 , 203 schnrich , r. , binney , j. , & dehnen , w. 2010 , mnras , 403 , 1829 seth , a.  c. , dalcanton , j.  j. , & de jong , r.  s. 2005 , aj , 130 , 1574 shapiro , k. , gerssen , j. , & van der marel , r. 2003 , apj , 126 , 2707 siebert , a. , bienaym , o. , binney , j. , et  al .
2008 , mnras , 391 , 793 soubiran , c. , bienaym , o. , & siebert , a. 2003 , a&a , 398 , 141 spagna , a. , lattanzi , m.  g. , re  fiorentin , p. , & smart , r.  l. 2010 , a&a , 510 , 4 skrutskie , m.  f. , cutri , r.  m. , stiening , r. , et  al .
2006 , aj , 131 , 1163 smith , m.  c. , evans , n.  w. , belokurov , v. et  al .
2009 , mnras , 399 , 1223 tonry , j. & davis , m. 1979 , aj , 84 , 1511 vallenari , a. , pasetto , s. , bertelli , g. , et  al .
2006 , a&a , 451 , 125 van der kruit , p.  c. , & searle , l. 1981 , a&a , 95 , 105 van der kruit , p.  c. , & searle , l. 1982 , a&a , 110 , 61 veltz , l. , bienaym , o. , freeman , k.  c. , et  al .
2008 , a&a , 2008 , 480 , 753 villalobos , a. , & helmi , a. 2008 , mnras , 391 , 1806 villalobos , a. , & helmi , a. 2009 , mnras , 399 , 166 villalobos , a. , kazantzidis , s. , & helmi , a. 2010 , apj , 718 , 314 wyse , r.  f.  g. & gilmore , g. 2005 , astro - ph/0510025 yi , s. k. , kim , y .- c .
, & demarque , p. 2003 , apjs , 144 , 259 yoss , k.  m. , neese , c.  l. , & hartkopf , w.  i. 1987 , aj , 94 , 1600 yoss , k.  m. , & griffin , r.  f. 1997 , japa , 18 , 161 walker , i.  r. , mihos , j.  c. , & hernquist , l. 1996 , apj , 460 , 121 zinn , r. 1985 , apj , 293 , 424 zwitter , t. , siebert , a. , munari , u. , et  al .
2008 , aj , 136 , 421