author: C. Boeche, A. Siebert, T. Piffl, A. Just, M. Steinmetz, S. Sharma, G. Kordopatis, G. Gilmore, C. Chiappini, M. Williams, E. K. Grebel, J. Bland-Hawthorn, B. K. Gibson, U. Munari, A. Siviero, O. Bienaymé, J.F. Navarro, Q. A. Parker, W. Reid, G. M. Seabroke, F. G. Watson, R. F. G. Wyse, T. Zwitter
title: Chemical gradients in the Milky Way from the RAVE data. I. Dwarf stars.

We aim at measuring the chemical gradients of the elements Mg, Al, Si, and Fe along the Galactic radius to provide new constraints on the chemical evolution models of the Galaxy and Galaxy models such as the Besançon model. Thanks to the large number of stars of our RAVE sample we can study how the gradients vary as function of the distance from the Galactic plane. We analysed three different samples selected from three independent datasets: a sample of 19 962 dwarf stars selected from the RAVE database, a sample of 10 616 dwarf stars selected from the Geneva-Copenhagen Survey (GCS) dataset, and a mock sample (equivalent to the RAVE sample) created by using the GALAXIA code, which is based on the Besançon model. We integrated the Galactic orbits and obtained the guiding radii ($R_\mathrm{g}$) and the maximum distances from the Galactic plane reached by the stars along their orbits ($Z_\mathrm{max}$). We measured the chemical gradients as functions of $R_\mathrm{g}$ at different $Z_\mathrm{max}$. The chemical gradients of the RAVE and GCS samples are negative and show consistent trends, although they are not equal: at $Z_\mathrm{max}$$<$`<!-- -->`{=html}0.4 kpc and 4.5$<$$R_\mathrm{g}$(kpc)$<$`<!-- -->`{=html}9.5, the iron gradient for the RAVE sample is $d[Fe/H]/dR_\mathrm{g}=-0.065$ dex kpc$^{-1}$, whereas for the GCS sample it is $d[Fe/H]/dR_\mathrm{g}=-0.043$ dex kpc$^{-1}$ with internal errors of $\pm0.002$ and $\pm0.004$ dex kpc$^{-1}$, respectively. The gradients of the RAVE and GCS samples become flatter at larger $Z_\mathrm{max}$. Conversely, the mock sample has a positive iron gradient of $d[Fe/H]/dR_\mathrm{g}=+0.053\pm0.003$ dex kpc$^{-1}$ at $Z_\mathrm{max}$$<$`<!-- -->`{=html}0.4 kpc and remains positive at any $Z_\mathrm{max}$. These positive and unrealistic values originate from the lack of correlation between metallicity and tangential velocity in the Besançon model. The discrepancies between the observational samples and the mock sample can be reduced by i) decreasing the density, ii) decreasing the vertical velocity, and iii) increasing the metallicity of the thick disc in the Besançon model.

# Introduction

The spatial distribution of the chemical abundances in the Milky Way provides important constraints for our understanding of galaxy formation. The distribution of the chemical species is a tracer of the way the galaxies formed because the abundances in the stellar atmospheres of a stellar population do not change during its lifetime[^1], in contrast to its spatial distribution and kinematics. Stars inherit the metallicity of the interstellar medium from which they form. Its metallicity depends in a complex way on the infall and outflow of gas and the enrichment by older stellar populations. Additionally, the element abundance ratios depend on the speed of enrichment. Fast enrichment results in an $\alpha$-enhancement because $\alpha$-elements are mainly produced by supernovae (SNe) type II on a short timescale, whereas most of the iron is produced by SN type Ia with a significant time delay. Therefore the abundances of stellar populations are tightly connected to the gas infall and star formation history. As a consequence, two locations that experienced different star formation histories exhibit a difference (gradient) in chemical abundances.
Today the existence of a negative chemical gradient along the Galactic radius is well established thanks to several dedicated observational works (see Maciel & Costa and references therein). Chemical gradients were measured by means of several tracers (such as Cepheids, planetary nebulae, open clusters) that cover a large range of Galactocentric distances ($R$). The chemical gradients found by previous investigations can span quite a wide range (from $-0.029$ dex kpc$^{-1}$ by Andrievsky et al., using Cepheids, to $-0.17$ dex kpc$^{-1}$ by Sestito et al., using open clusters). Most of the works seem to converge on a value close to $\sim-0.06$ dex kpc$^{-1}$ by using different tracers (Cepheids by Luck et al. , Luck & Lambert , Yong et al. , Lemasle et al. among others; open clusters by Friel et al. , Sestito et al. , Pancino et al. among others; planetary nebulae by Pasquali & Perinotto , Maciel & Quireza ; main-sequence turn-off stars by Cheng et al. , Coşkunoǧlu et al. ). It has been shown that the radial gradient becomes flatter with increasing distance from the Galactic plane (Pasquali & Perinotto ; Boeche ; Cheng et al. ) . If this flattening is caused by a zero radial gradient in the thick disc (; Ruchti et al. ; Bilir et al. ) or an effect of a flattening agent such as radial mixing (as dynamical and chemo-dynamical models would suggest, see Sellwood & Binney ; Schönrich & Binney ; Sánchez-Blázquez et al. ; Minchev et al. , ) and energy feedback (Gibson et al. ) is still unclear.
To explain the gradients observed in the Milky Way, an inside-out formation scenario of the Milky Way has been proposed (Matteucci & François, ) in which the inner parts of the Galaxy experienced a higher star formation rate than the outer parts. Galactochemical models (Chiappini et al. , Cescutti et al., ) and hydrodynamical simulations (Pilkington et al. ) with inside-out formation show radial gradients consistent with the observed one. The study by Schönrich & Binney claims that there is no need to postulate an inside-out star formation to obtain negative gradients. In their model the whole disc forms simultaneously, and a metallicity gradient of $\sim-0.11$ dex kpc$^{-1}$ is finally obtained because of the advection inwards of heavier elements released by SNe. Observational data (see Sestito et al. , Luck et al. , Carraro et al. , Yong et al. among others) suggest that there are variations of the gradient along $R$: a steeper gradient in the inner disc ($R$$\lesssim5$ kpc) and a flatter gradient in the outer disc ($R$$\gtrsim10$ kpc).
In this paper we make use of the RAVE data to investigate the radial abundance gradients of the Milky Way. RAVE is a large spectroscopic survey that collected 574 630 spectra of 483 330 stars of the Milky Way. The spectra are centred on the near-infrared triplet window (8410-8795Å). The resolution of R$\simeq$`<!-- -->`{=html}7500 allows measuring radial velocities (Steinmetz et al., ), stellar parameters (Zwitter et al. ; Siebert et al. ; Kordopatis et al. ), and chemical abundances (Boeche et al. ; Kordopatis et al. ). Distance estimates (Breddels et al. ; Zwitter et al. ; Burnett et al. ; Binney et al. ) enable one to locate the stars in three-dimensional space. Proper motions are listed in a variety of catalogues: Tycho2 (Høg et al., ), the PPM-Extended catalogues PPMX and PPMXL (Roeser et al., , ), and the second and third U.S. Naval Observatory CCD Astrograph Catalog UCAC2 and UCAC3 (Zacharias et al., ). Chemical abundances for six elements are provided by the RAVE chemical pipeline (Boeche et al. ; Kordopatis et al. ), which delivers abundances with errors of $\sim0.10-0.35$ dex, depending on the element, stellar parameters, and signal-to-noise ratio (S/N). The full set of kinematic and chemical data is a rich source of information for studying the Galactic disc and has been recently employed for several studies. By using the RAVE space and velocity data set, Siebert et al. found a gradient in the mean Galactocentric radial velocity of stars in the extended solar neighbourhood, which has been confirmed by a recent work of Williams et al. who, using red-clump giant stars, found that it is more marked in the southern part of the disc, and furthermore highlights wave-like structures in the disc from the velocity distributions of the sample. The structure of the disk was studied by Pasetto et al. (, ), Antoja et al. , Wilson et al. , and Veltz et al. . The full set of RAVE chemo-kinematical data was recently used by Boeche et al. to study the relation between abundances and kinematics of the Galactic disc populations. The radial metallicity gradient was previously studied in recent works by using RAVE data (Karataş & Klement , Coşkunoǧlu et al. ), but none of them investigated the gradients for the individual elemental abundances. Recently, Bilir et al. split a sample of red-clump RAVE stars into thin- and thick-disc subsamples by labelling the stars as a function of their kinematic characteristics, and measured the metallicity gradients in these subsamples.
In the present work we investigate the chemical gradients of four elements over the Galactocentric distance range of 4.5–9.5 kpc by using the chemical and kinematic data of dwarf stars in the RAVE database. We also investigate how the gradients change as function of the distance from the Galactic plane. We applied the same analysis to the Geneva-Copenhagen Survey data set (Nordström et al. ) with the recent data refinement by Casagrande et al. and with a mock sample constructed with the code GALAXIA (Sharma et al., ) to test the robustness of our results and compare them with Galactic models. We discuss how different stellar populations in the sample can generate a bias and affect the gradients measurements. This paper will be followed by a second paper in which we study the chemical gradients by using RAVE giant stars.
The paper is structured as follows: Section 2 describes the data and the selection of the star samples used to measure the gradients, in Section 3 we explain the analysis method applied to the samples, in Section 4 and 5 we report the chemical gradients obtained, in Section 6 we discuss biases and results and we conclude in Section 7.

# Data

## RAVE sample

We selected the sample from the RAVE internal data archive, which contains stellar atmospheric parameters ($T_{\fontsize{6}{6}\selectfont \mbox{eff}}$, $\log g$ and \[M/H\]) of 451 413 spectra, and chemical elemental abundances for 450 641 spectra. The stellar atmospheric parameters are measured with the recently adopted RAVE pipeline (Kordopatis et al. , and Kordopatis at al. ). The chemical abundances (\[X/H\]) for Fe, Al, Mg, Si, Ti, and Ni are derived by the RAVE chemical pipeline (Boeche et al. , Kordopatis et al. ). We adopted the parallaxes by Binney et al. (submitted) computed with an improved version of the recipe described in Burnett et al. .
With these data we determined Galactocentric positions and absolute velocities. Galactic orbital parameters such as the guiding radius ($R_\mathrm{g}$) and maximum distance from the Galactic plane ($Z_\mathrm{max}$) were considered and computed as explained in Sec. . In this first paper we limit our investigation to dwarf stars with good-quality RAVE spectra (i.e. good stellar parameters and chemical abundances) and reliable distances. The RAVE chemical catalogue provides abundances with errors $\sigma\sim0.2$ dex for spectra with an S/N$>$`<!-- -->`{=html}40. These errors hold for the elements Fe, Mg, Al, and Si (Ti and Ni measurements have low precision for dwarf stars). Therefore, our selection criteria were i) stars with spectra with an S/N$>$`<!-- -->`{=html}40, ii) dwarf stars with gravity $\log g$$>$`<!-- -->`{=html}3.8 dex and effective temperatures $T_{\fontsize{6}{6}\selectfont \mbox{eff}}$ from 5250 to 7000 K, and iii) stars with distance uncertainties smaller than 30%. Moreover, we considered only the four elements Fe, Mg, Al, and Si because of their good precision ($\sigma\sim0.15$ dex) within the selected $T_{\fontsize{6}{6}\selectfont \mbox{eff}}$ and S/N ranges. We took care that no peculiar stars were included in the sample by selecting only spectra flagged as normal stars by Matijevič et al. .
Among these stars, we furthermore selected only spectra for which the RAVE pipeline converged to a single point of the parameter space ($Algo\_Conv$ parameter equal to zero). These spectra are well fit by the reconstructed spectra of the RAVE chemical pipeline ($\chi^2<2000$), and with minor or no continuum defects[^2] ($frac>0.7$). Since RAVE re-observed $\sim$`<!-- -->`{=html}10% of the objects, some stars have two or more spectroscopic observations. To avoid multiple observations in our sample, we chose only the spectrum with the highest S/N for re-observed stars. This selection yielded 19 962 stars, whose spatial distribution is shown in Figure . Most of the stars lie within 0.3 kpc from the Galactic plane and in the Galactic radius interval $R$=7.6–8.3 kpc. This is a small interval where we can estimate the local gradient, although its small spatial coverage does not allow precise measurements. We extended the gradient measurements to the Galactic radius range $R\sim$`<!-- -->`{=html}4.5-9.5 kpc by using the guiding radius $R_\mathrm{g}$.

## Geneva-Copenhagen Survey sample

The GCS provides information such as the temperature, metallicity, kinematics, distance, and age of 16 682 F and G dwarf stars in the solar neighbourhood (Nordström et al., ). We here refer to the more recent values reported by Casagrande et al. . We selected stars by using the same constraints in $T_{\fontsize{6}{6}\selectfont \mbox{eff}}$, $\log g$, and guiding radii $R_\mathrm{g}$ as were used for RAVE stars to keep the characteristics of the GCS and RAVE samples as similar as possible. With these constraints we selected 10 616 stars. Still, there are some differences. In fact, GCS stars are spherically distributed in space and most of them do not lie farther away than 0.2 kpc from the Galactic plane, whereas the RAVE stars have a cone-shaped distribution in space and extend up to 0.5 kpc from the Galactic plane (Figure ). Because most of the RAVE sample is inside the scale height of the thin disc ($\sim$`<!-- -->`{=html}0.3 kpc), it is reasonable to assume that the difference in stellar populations between the two samples is small. We expect a slight underrepresentation of the young, dynamically cold sub population in the RAVE sample compared with the GCS sample because of the cone structure of the RAVE volume. In Figure  we compare the $T_{\fontsize{6}{6}\selectfont \mbox{eff}}$ and \[Fe/H\] distributions of the two samples. The GCS sample has a temperature distribution favouring higher $T_{\fontsize{6}{6}\selectfont \mbox{eff}}$ than the RAVE sample, and the \[Fe/H\] distribution is shifted by $\sim$`<!-- -->`{=html}0.1 dex towards higher metallicities. Despite the highlighted differences, the GCS sample shows slightly flatter chemical gradients (but still in reasonable agreement) than the RAVE sample, as we will see below.

## RAVE mock sample

The magnitude and spatial distributions of the RAVE sample may generate observational biases that can lead to misinterpretations. To avoid this, we created an equivalent mock sample with the stellar population synthesis code GALAXIA , which uses analytical density profiles based on the Besançon model . The mock sample reproduces the detailed RAVE selection function in I magnitude and target density, both as functions of angular position on the sky[^3]. Only RAVE spectra with S/N$>40$ and flagged as normal stars by Matijevič et al. were considered.
The resulting mock catalog agrees well with the original RAVE catalog for the stellar parameters and the distance and radial velocity distributions. From this sample of 97 485 entries, we imposed the same selection criteria in $T_{\fontsize{6}{6}\selectfont \mbox{eff}}$ and $\log g$ to mimic our RAVE sample, obtaining a mock sample of 26 198 entries. The higher number of stars in the mock sample than the RAVE sample is due to other selection criteria, that are applied to the RAVE sample (such as $Algo\_Conv$ parameter convergence, $\chi^2$, $frac$ parameters and a distance estimate) and not to the underlying Galaxy model. The selected cut in S/N (S/N$>$`<!-- -->`{=html}40) helps minimising these differences, which should not affect the results on the abundance gradients significantly. In fact, although the spectrum quality parameters change as a function of the distance (because RAVE is a magnitude-limited survey with constant integration time, more distant stars have on average lower S/N spectra), this will equally affect stars at any Galactocentric distance. However, to stay on the safe side, we verified that the results were consistent with those obtained by applying the cut when no cuts in distance error were applied to the RAVE data.

# Method and error estimation

## Galactic orbit integration

Because stars are oscillating radially and vertically along their orbits, stellar samples in the solar neighbourhood contain information for a much larger volume. To take advantage of this, we integrated the Galactic orbits for a time of 40 Gyr (corresponding to more than 50 periods for most of the stars) with the code NEMO . As Galactic potential we adopted the potential model n.2 by Dehnen & Binney , which assumes R$_\odot$=8.0 kpc, and best matches the observed properties of the Galaxy. We used the same potential for all samples, including the mock sample based on the Besancon model. In axisymmetric equilibrium models the z-component of the angular momentum $L\mathrm{z}$ and the vertical action $J\mathrm{z}$ are constants of motion. With the help of the rotation curve $v_\mathrm{c}(R_\mathrm{g})$ the angular momentum can be converted into the guiding radius $R_\mathrm{g}$. The commonly used mean radius $R_\mathrm{m}$ of apocentre $R_\mathrm{a}$ and pericentre $R_\mathrm{p}$ is a proxy for the guiding radius only for low eccentricities, therefore $R_\mathrm{g}$ has to be preferred. With these quantities we can study the gradients at Galactic radii beyond the range currently occupied by the stars. Vertically, the longest distance to the Galactic plane $Z_\mathrm{max}$ is a measure of the vertical action.
The values $R_\mathrm{g}$ and $Z_\mathrm{max}$ were computed with the method just described for the RAVE, the GCS, and the mock samples.

## Gradient error estimation

We measured the abundance gradients by fitting the distribution of the stars in the ($R_\mathrm{g}$,\[X/H\]) plane with a linear regression. We also investigated more sophisticated methods such as modelling the density distribution of the stars with a two-dimensional function such as a Gaussian or lognormal function, and performing the fitting procedure with the maximum-likelihood method. The different methods provided consistent results. The linear regression was ultimately preferred because it is simpler and faster, in particular during the confidence interval evaluation. The gradient confidence intervals reported in this paper were evaluated by using a resampling method such as bootstrapping. Means and 1$\sigma$ errors (68% confidence) were computed over 1000 realizations for every gradient reported. These errors represent the internal errors, as explained in the following section.

## Bias in the gradients and gradient errors

One needs to be aware that the use of the guiding radius $R_\mathrm{g}$ introduces a kinematic bias for volume-limited samples such as the RAVE and GCS samples. Stars with small $R_\mathrm{g}$ can reach the solar neighbourhood only if they move in eccentric (kinematically hot) orbits. Such stars are on average more metal-poor than stars at the same $R_\mathrm{g}$ and in circular orbits, which will never reach the solar neighbourhood. This causes a lack of metal-rich stars at small $R_\mathrm{g}$ in volume-limited samples, and for these samples the measured gradient is less negative. The lack of metal-rich stars occurs at large $R_\mathrm{g}$ as well, with the difference that in the outer Galaxy the stellar density is lower, therefore the effect is weaker. Although we cannot precisely quantify this bias, we discuss the impact of this bias on the measured gradients.

# Analysis and results

In this section we split the samples into three ranges of $Z_\mathrm{max}$, discuss the \[Fe/H\] distributions, and derive the radial gradients in \[Fe/H\] and in \[$\alpha$/Fe\]. In Sec.  we refine the method for the RAVE sample to improve the vertical resolution and extend the analysis to the elements Fe, Al, Mg, and Si and the relative abundance gradients \[X/Fe\].

## Gradients from the RAVE sample

The resulting gradient for iron at the present position of the stars is $\frac{d [Fe/H]}{d R}= -0.029\pm0.011$ dex kpc$^{-1}$. The small range in $R$ covered by the sample (7.6$<$$R$(kpc)$<$`<!-- -->`{=html}8.3) does not allow a precise measurement. The derived metallicity gradient is significantly flatter than expected, but this may be due to the conical volume, where a negative vertical gradient biases the radial gradient to lower values. Moreover, $R$ may not be representative of the distance at which the stars were born, since many of them do not have circular orbits. From the probability point of view, the guiding radius $R_\mathrm{g}$ better represents the Galactocentric distance of the star during one revolution around the Galactic centre. By using $R_\mathrm{g}$ and the constraints 4.5$<R_\mathrm{g}$(kpc)$<$`<!-- -->`{=html}9.5 (to avoid the few outlier stars that can affect the estimation) we obtained the result
$\frac{d [Fe/H]}{d R_\mathrm{g}}= -0.059\pm0.002$ dex kpc$^{-1}$.
We found in a previous work (Boeche, ) that different gradients are obtained at different $Z_\mathrm{max}$. We divided the full sample into three subsamples with different $Z_\mathrm{max}$ranges: $0.0<Z_\mathrm{max}$ (kpc)$\leq 0.4$ (16 456 stars, named $Z_\mathrm{0.0}^\mathrm{RAVE}$ sample) $0.4<Z_\mathrm{max}$ (kpc)$\leq 0.8$ (3032 stars, the $Z_\mathrm{0.4}^\mathrm{RAVE}$ sample), and $Z_\mathrm{max}$ (kpc)$>0.8$ (399 stars, the $Z_\mathrm{0.8}^\mathrm{RAVE}$ sample). For these samples we found
$\frac{d [Fe/H]}{d R_\mathrm{g}}
\mbox{($Z_\mathrm{0.0}^\mathrm{RAVE}$)}=-0.065\pm0.002$ dex kpc$^{-1}$
$\frac{d [Fe/H]}{d R_\mathrm{g}}\mbox{($Z_\mathrm{0.4}^\mathrm{RAVE}$)}=-0.059\pm0.005$ dex kpc$^{-1}$
$\frac{d [Fe/H]}{d R_\mathrm{g}}\mbox{($Z_\mathrm{0.8}^\mathrm{RAVE}$)}=0.005\pm0.012$ dex kpc$^{-1}$.
In Figure  the gradients are represented by the slope of the straight line (left panels). We also report the gradients for \[$\alpha$/Fe\][^4]:
$\frac{d[\alpha/Fe]}{dR_\mathrm{g}}\mbox{($Z_\mathrm{0.0}^\mathrm{RAVE}$)}
=-0.004\pm0.001$ dex kpc$^{-1}$,
$\frac{d[\alpha/Fe]}{dR_\mathrm{g}}\mbox{($Z_\mathrm{0.4}^\mathrm{RAVE}$)} =
-0.005\pm0.002$ dex kpc$^{-1}$,
$\frac{d[\alpha/Fe]}{dR_\mathrm{g}}\mbox{($Z_\mathrm{0.8}^\mathrm{RAVE}$)} =
-0.020\pm0.005$ dex kpc$^{-1}$.
These results confirm that the radial gradients in the RAVE sample vary with increasing $Z_\mathrm{max}$. The complete set of measured gradients of abundances \[X/H\] and relative abundance \[X/Fe\] for the RAVE elements are given in Tables  and .

## Gradients from the Geneva-Copenhagen survey sample

To evaluate the robustness of our results we repeated the same procedure for the GCS data. With the actual Galactocentric radius $R$ of all stars, the gradient is $\frac{d [Fe/H]}{d R}=-0.004\pm0.047$ dex kpc$^{-1}$, which is unreliable because the very small radial extent of the sample.
With the guiding radius in the range $R_\mathrm{g}$$\sim$`<!-- -->`{=html}4.5-9.5 kpc we obtained a gradient $\frac{d [Fe/H]}{d R_\mathrm{g}}=-0.027\pm0.004$ dex kpc$^{-1}$ .
We followed the procedure used for the RAVE sample and selected three subsamples in the following $Z_\mathrm{max}$ ranges: $0.0<Z_\mathrm{max}$ (kpc)$<0.4$ (8130 stars, named the $Z_\mathrm{0.0}^\mathrm{GCS}$ sample) $0.4<Z_\mathrm{max}$ (kpc)$<0.8$ (621 stars, the $Z_\mathrm{0.4}^\mathrm{GCS}$ sample) and $Z_\mathrm{max}$ (kpc)$>0.8$ (197 stars, the $Z_\mathrm{0.8}^\mathrm{GCS}$ sample). The resulting gradients are (see also Figure , middle)
$\frac{d[Fe/H]}{dR_\mathrm{g}}\mbox{($Z_\mathrm{0.0}^\mathrm{GCS}$)} =-0.043\pm0.004$ dex kpc$^{-1}$,
$\frac{d[Fe/H]}{dR_\mathrm{g}}\mbox{($Z_\mathrm{0.4}^\mathrm{GCS}$)} =
-0.008\pm0.011$ dex kpc$^{-1}$,
$\frac{d[Fe/H]}{dR_\mathrm{g}}\mbox{($Z_\mathrm{0.8}^\mathrm{GCS}$)} =
+0.056\pm0.019$ dex kpc$^{-1}$.
The metallicity distribution shows a more densely populated low-metallicity tail at large $Z_\mathrm{max}$ than the RAVE sample. The gradients $d[Fe/H]/dR_\mathrm{g}$ are in general less negative than the RAVE sample. The trend of a flatter gradient with $Z_\mathrm{max}$ agrees with the RAVE data, with the exception of the $Z_\mathrm{0.8}^\mathrm{GCS}$ subsample (which has a positive gradient), which is particularly affected by the bias discussed in Sec. . This gradient is probably not reliable, and in Sec.  we provide a more realistic value.
We tentatively measured the gradients of the $\alpha$-elements given by Casagrande et al. , with the warning that these abundances were estimated from photometry and must be considered only a proxy of \[$\alpha$/Fe\]. For \[$\alpha$/Fe\] we obtained
$\frac{d[\alpha/Fe]}{dR_\mathrm{g}}\mbox{($Z_\mathrm{0.0}^\mathrm{GCS}$)}
=+0.010\pm0.002$ dex kpc$^{-1}$,
$\frac{d[\alpha/Fe]}{dR_\mathrm{g}}\mbox{($Z_\mathrm{0.4}^\mathrm{GCS}$)} =
-0.006\pm0.005$ dex kpc$^{-1}$,
$\frac{d[\alpha/Fe]}{dR_\mathrm{g}}\mbox{($Z_\mathrm{0.8}^\mathrm{GCS}$)} =
-0.023\pm0.007$ dex kpc$^{-1}$.
Using only stars with more reliable temperature estimates (treated with the infrared flux method of Casagrande et al. ), the results do not change significantly (differences smaller than 1$\sigma$).

## Gradients from the RAVE mock sample

The mock sample allows us to separate the influence of the different components (thin disc, thick disc, stellar halo) on the properties of the RAVE sample. In Figure  (bottom panel) the contribution of the three components are separately shown in every $Z_\mathrm{max}$ slice. The halo contribution is negligible, and the thick disc is detected as a prominent low-metallicity peak at $Z_\mathrm{max}$$>0.8$ kpc. The Besançon model assumes a radial metallicity gradient of d\[Fe/H\]/dR= –0.07 dex kpc$^{-1}$ for the thin disc and no gradient in the thick disc. The gradients we found deviate significantly from these values for reasons that we analyse in Sec. . Using the actual galactocentric radius $R$, the gradient for the mock sample is $\frac{d [Fe/H]}{d R}=-0.156\pm0.010$ dex kpc$^{-1}$, significantly more negative than expected. This is probably caused by the uneven spatial coverage in RAVE.
Using the guiding radius $R_\mathrm{g}$, we obtained a gradient $\frac{d [Fe/H]}{d R_\mathrm{g}}=+0.060\pm0.003$ dex kpc$^{-1}$ .
Selecting three subsamples in the three different $Z_\mathrm{max}$ ranges as we did before, the resulting gradients are
$\frac{d[Fe/H]}{dR_\mathrm{g}}\mbox{($Z_\mathrm{0.0}^\mathrm{mock}$)}
=+0.053\pm0.003$ dex kpc$^{-1}$,
$\frac{d[Fe/H]}{dR_\mathrm{g}}\mbox{($Z_\mathrm{0.4}^\mathrm{mock}$)} =
+0.065\pm0.005$ dex kpc$^{-1}$,
$\frac{d[Fe/H]}{dR_\mathrm{g}}\mbox{($Z_\mathrm{0.8}^\mathrm{mock}$)} =
+0.063\pm0.008$ dex kpc$^{-1}$.
We were unable measure the \[$\alpha$/Fe\] gradient for the mock sample because the Besançon model does not provide this observable.

```latex
\begin{table*}[t]
\caption[]{Radial abundance gradients measured in the RAVE sample for Fe, Mg, Al, and Si
expressed as dex kpc$^{-1}$ for three
ranges of $Z_\mathrm{max}$. Uncertainties of 68\% confidence are
obtained with the bootstrap method.
}
\label{tab_XH_grad_Rg}
\vskip 0.3cm
\centering
\begin{tabular}{l|ccccc}
\hline
\noalign{\smallskip}
        &  $\frac{d[Fe/H]}{dR_\mathrm{g}}$  & $\frac{d[Mg/H]}{dR_\mathrm{g}}$  & $\frac{d[Al/H]}{dR_\mathrm{g}}$
& $\frac{d[Si/H]}{dR_\mathrm{g}}$ \\
\noalign{\smallskip}
\hline
\noalign{\smallskip}
0.0$\leq Z_\mathrm{max}\mathrm{(kpc)}<$0.4           &$-0.065\pm0.003$&$-0.073\pm0.002$&$-0.085\pm0.003$&$-0.064\pm0.003$\\
%\noalign{\smallskip}

0.4$\leq Z_\mathrm{max}\mathrm{(kpc)}<$0.8           &$-0.059\pm0.006$ &$-0.061\pm0.005$&$-0.076\pm0.007$ &$-0.063\pm0.005$\\
%\noalign{\smallskip}

\phantom{0.4$\leq$}$Z_\mathrm{max}\mathrm{(kpc)}>$0.8& $+0.006\pm0.015$&$-0.026\pm0.011$&$-0.017\pm0.016$ & $-0.008\pm0.013$\\
%\noalign{\smallskip}
\noalign{\smallskip}
\hline
\end{tabular}

\end{table*}
```
```latex
\begin{table*}[t]
\caption[]{As Table~\ref{tab_XH_grad_Rg}, but for relative abundances [X/Fe].}
\label{tab_XFe_grad_Rg}
\vskip 0.3cm
\centering
\begin{tabular}{l|ccccc}
\hline
\noalign{\smallskip}
        &  & $\frac{d[Mg/Fe]}{dR_\mathrm{g}}$ &
$\frac{d[Al/Fe]}{dR_\mathrm{g}}$ & $\frac{d[Si/Fe]}{dR_\mathrm{g}}$ \\
\noalign{\smallskip}
\hline
\noalign{\smallskip}
0.0$\leq Z_\mathrm{max}\mathrm{(kpc)}<$0.4            &       &$-0.009\pm0.002$&$-0.022\pm0.002$&$+0.001\pm0.001$\\
%\noalign{\smallskip}

0.4$\leq Z_\mathrm{max}\mathrm{(kpc)} <$0.8            &       & $-0.004\pm0.003$&$-0.018\pm0.004$ & $-0.004\pm0.003$\\
%\noalign{\smallskip}

\phantom{0.4$\leq$}$ Z_\mathrm{max}\mathrm{(kpc)} >$0.8 &       &$-0.027\pm0.007$&$-0.024\pm0.009$ &$-0.013\pm0.006$\\
%\noalign{\smallskip}
\noalign{\smallskip}
\hline
\end{tabular}
%\normalsize

\end{table*}
```

# Gradient estimates with moving box car

The previous results show that the gradient flattens when moving to higher $Z_\mathrm{max}$. Whether this change occurs smoothly or with a sudden transition can be investigated by measuring the gradient in an interval of constant width, centred on progressively larger distances from the plane.
We measured the gradient $\frac{d[X/H]}{dR_\mathrm{g}}$ of a sample of stars lying in the interval \[$Z_\mathrm{max}$-0.1,$Z_\mathrm{max}$+0.1\] and moved this sampling interval by 0.1 kpc steps from 0.1 kpc towards larger $Z_\mathrm{max}$. To always have a statistically significant sample in these intervals, we imposed the condition that at every step the sample must contain no fewer than 200 stars. If there were fewer stars, the $Z_\mathrm{max}$ interval was extended until the condition was matched. Undersampling occurs in particular at $Z_\mathrm{max}$$>$`<!-- -->`{=html}0.7 kpc, where the number of stars is small. We present the result in Figure , where the gradient $\frac{d[X/H]}{dR_\mathrm{g}}$ is shown as a function of $Z_\mathrm{max}$. We found that the gradients remained constant for all elements until $Z_\mathrm{max}$$\sim 0.7$ kpc, then they increase to shallower (or positive) gradients with a hint of a step occurring between 0.6 and 0.9 kpc. For \[Fe/H\] the gradient of the RAVE sample is –0.06 dex kpc$^{-1}$ and $\sim0$ at small and large $Z_\mathrm{max}$, respectively, slightly more negative than the GCS sample, where it amounts to –0.05 and +0.02 dex kpc$^{-1}$, respectively. The mock sample shows positive gradients at any $Z_\mathrm{max}$. The value of +0.06 dex kpc$^{-1}$ at low $Z_\mathrm{max}$ decreases to +0.02 at large $Z_\mathrm{max}$, where the thick disc dominates.
In Figure  we show the abundance distributions \[X/H\] as a function of $Z_\mathrm{max}$ for the four elements of the RAVE sample and for \[Fe/H\] in the GCS and the mock sample. In this figure we allowed a minimum of 20 stars in the interval to see local changes in the mean abundance. There is a clear negative trend visible in all cases. However, the trend of \[Fe/H\] for the GCS sample is slightly smoother but covers the same range. The mock sample behaves quite differently: although the mean iron abundance at small and large $Z_\mathrm{max}$ is similar to that of the RAVE and GCS samples ($\sim0.6$ dex difference between $Z_\mathrm{max}$=0 and $Z_\mathrm{max}$=2 kpc), its evolution as a function of $Z_\mathrm{max}$ is quite different. While for RAVE and GCS the decrease of the mean Fe abundance in the range 0.0$<$$Z_\mathrm{max}$$<$`<!-- -->`{=html}1.6 kpc appears to be approximately linear, for the mock sample it has a negative second derivative up to $Z_\mathrm{max}$$\sim$`<!-- -->`{=html}0.6 kpc and becomes positive afterwards. This difference suggests that the transition between the thin and thick disc in RAVE and GCS differs from the prediction of the mock sample.

# Discussion

An important goal of investigating the chemical gradients is to distinguish the effect of biases caused by the selection function, the impact of the thin disc/thick disc dichotomy, and the features of the intrinsic evolution of the thin disc by an inside-out growth (Matteucci & François, ). The three samples analysed here show some similarities, but also striking differences.

## Abundance gradients in the GALAXIA/Besançon mock sample

In the local sample of any Milky Way model the spatial structure of components or sub populations are represented only indirectly and only if the model is set up in a dynamically consistent way. The radial scalelength of the density in the disc is represented by a proper choice of the asymmetric drift (the shift of the tangential velocity distribution with respect to the circular speed). Vertically, the thickness corresponds to the assigned velocity dispersion. In addition to the density distributions and kinematics, stellar metallicities are assigned to each sub population according to their current position, in the Besançon model, regardless of kinematics. The vertical metallicity distribution is set up in a consistent way, because each sub population is assigned a spatially homogeneous metallicity distribution. Radially, an uncorrelated local metallicity distribution would also require a homogeneous metallicity distribution. The additionally assigned spatial metallicity gradient of d\[Fe/H\]/dR=–0.07 dex kpc$^{-1}$ for the thin-disc sub populations would require the correlation of the metallicity with the tangential velocity or angular momentum to be dynamically consistent. This means that the local sample with no correlation of tangential velocity and metallicity looks like a sample of components with no radial metallicity gradient. This caveat could easily be overcome by assigning the metallicity distribution at each location as a function of the angular momentum (or tangential velocity), leading to a radial gradient with respect to $R_\mathrm{g}$. The radial gradient as a function of the current position $R$ would then be very similar (see e.g. Binney for a dynamical model).

Nevertheless, we can learn about the impact of a distinct thick disc with lower metallicity and higher velocity dispersion and stronger asymmetric drift on the mock sample. The thick disc with a mean metallicity of –0.78 dex and 0.3 dex internal dispersion has an asymmetric drift of 53 km s$^{-1}$ and a tangential velocity dispersion of 51 km s$^{-1}$ compared with 15 and 28 km s$^{-1}$ for the old thin disc. This shifts the thick-disc stars to the lower left in the ($R_\mathrm{g}$,\[Fe/H\]) diagrams (Figure  bottom) with respect to the thin-disc stars, generating an apparent decreasing mean \[Fe/H\] with decreasing $R_\mathrm{g}$, or in other words, a fictitious positive gradient. This trend is stronger for larger $Z_\mathrm{max}$, because the fraction of thick-disc stars increases with a still dominating thin disc at large $R_\mathrm{g}$. For the analysis of thin disc properties the contamination by thick disc stars in the RAVE sample results in a bias towards less negative or even positive radial metallicity gradients with respect to $R_\mathrm{g}$.

An additional striking difference of the mock sample compared with the RAVE and the GCS samples is the large number of low-metallicity stars for $Z_\mathrm{max}$$>$`<!-- -->`{=html}0.8 kpc (see Figure ). The reasons for this deviation can be the i) thick-disc density, ii) its vertical velocity dispersion, and iii) its mean metallicity adopted by the Besançon model. A lower local density in the model (at present assumed to be $\sim$`<!-- -->`{=html}10%) would reduce the peak proportionally, a lower vertical velocity dispersion (at present assumed to be 42 km s$^{-1}$) would shift a significant part to the lower $Z_\mathrm{max}$ bin and slightly lift the low-metallicity tail, whereas a higher mean metallicity (at present assumed –0.78dex) would merge the peak with the thin-disc distribution.

Previous authors have already shown that the thick disc might have a higher mean metallicity than the one assumed in the Besançon model. A thick-disc mean metallicity of $\sim-0.5$ dex was found by some authors (; ; ; ) with a gradient that decreases the mean metallicity $\sim-0.8$ dex beyond 3 kpc (; ). The lack of a prominent low-metallicity tail in the RAVE and GCS samples compared with the mock sample supports these studies.

Despite the unrealistic gradients exhibited, the mock sample gives very precious information for data interpretation. It suggests that short scale-length metal-poor components (thick disc and halo stars) may accumulate at low $R_\mathrm{g}$, and, together with the thin-disc stars, cause a fictitious flatter or positive gradient that depends on the relative fraction of the stellar populations in the sample. This might be one of the possible explanations of the positive metallicity gradients in thick-disc-like stars found in some previous works (Nordström et al. ; Carrell et al. ). In observational data the population mix depends on the selection function of the sample under analysis, and only a complete (ideal) sample can deliver unbiased gradients. We must therefore be careful when interpreting the gradient measurement of any sample available today.

## Abundance gradients of the observational samples

In the previous section we saw how the mix of thin- and thick-disc stars in the Besançon model causes a positive gradient when the kinematically dependent quantity $R_\mathrm{g}$ is used. This effect is enhanced in our samples by the observational bias discussed in Sec. : metal-rich stars in (nearly) circular orbits at small $R_\mathrm{g}$ are missing because they cannot reach our volume-limited samples located in the solar neighbourhood, and the gradient becomes flatter than it should be. This bias is clearly visible in the GSC sample (Figure , middle panel) where stars at 4.5$\leq$$R_\mathrm{g}$(kpc)$\leq$`<!-- -->`{=html}5.5 are clearly \[Fe/H\] poorer than stars at higher $R_\mathrm{g}$. Excluding this $R_\mathrm{g}$ interval from the measurement, the gradients for the GCS data are $d [Fe/H]/d R_\mathrm{g}=-0.054\pm0.003$ dex kpc$^{-1}$, $-0.028\pm0.012$ dex kpc$^{-1}$, and $+0.004\pm0.025$ dex kpc$^{-1}$ for the subsamples $Z_\mathrm{0.0}^\mathrm{GCS}$, $Z_\mathrm{0.4}^\mathrm{GCS}$, and $Z_\mathrm{0.8}^\mathrm{GCS}$, respectively, whereas for the RAVE data they are $d [Fe/H]/d R_\mathrm{g}=-0.066\pm0.003$ dex kpc$^{-1}$, $-0.063\pm0.006$ dex kpc$^{-1}$, and $-0.002\pm0.016$ dex kpc$^{-1}$ for the subsamples $Z_\mathrm{0.0}^\mathrm{RAVE}$, $Z_\mathrm{0.4}^\mathrm{RAVE}$, and $Z_\mathrm{0.8}^\mathrm{RAVE}$, respectively. The new GCS results are closer to those found with the RAVE data and, in particular, the peculiar positive gradient found with the $Z_\mathrm{0.8}^\mathrm{GCS}$ sample disappears. We verified that adopting other cuts at low or high $R_\mathrm{g}$ causes changes in gradient of the order of 2-3 $\sigma$ for both the RAVE and the GCS samples (that means $\sim0.01$ dex kpc$^{-1}$). These results suggest that the errors obtained with the bootstrap method (explained in Sec. ) represent the internal errors, in which observational biases are neglected. Local inhomogeneities in the stellar population in the solar neighbourhood (such as disrupted open clusters and moving groups of kinematic origin, see Famaey et al. ; Antoja et al. ) can also explain the differences in gradient observed between the RAVE and the GCS samples and between different $R_\mathrm{g}$ intervals considered. Our conclusion is that these observational biases can affect the measured gradients for a quantity of about $\sim0.01$ dex kpc$^{-1}$, which is larger than our internal errors, but it does not affect the interpretation of our results.

## Radial chemical gradients of the RAVE sample as a function of $Z_\mathrm{max}$

In Section  we began the analysis by dividing the sample into three subsamples $Z_\mathrm{0.0}^\mathrm{RAVE}$, $Z_\mathrm{0.4}^\mathrm{RAVE}$, and $Z_\mathrm{0.8}^\mathrm{RAVE}$. In the $R_\mathrm{g}$ interval 4.5-9.5 kpc, the $Z_\mathrm{0.0}^\mathrm{RAVE}$ subsample exhibits a gradient of $\frac{d[Fe/H]}{dR_m}\sim-0.065$ dex, which agrees with some of the works cited in Sec. . Moving to the subsamples $Z_\mathrm{0.4}^\mathrm{RAVE}$ and $Z_\mathrm{0.8}^\mathrm{RAVE}$, the gradients become less negative and then flatten, as also found in previous works (Cheng et al. ; Boeche ; Pasquali & Perinotto ), but with these analyses it is not clear if the gradient changes smoothly or abruptly. Thanks to the large number of RAVE stars, we can use a moving box car to measure the radial gradients along $Z_\mathrm{max}$ (Figure ), and show that the gradients (of four different elements) are negative and remain approximately constant up to $Z_\mathrm{max}$$\sim0.6$ kpc, then they progressively increase to zero or positive values. In the light of the previous discussion of the Besançon model, the reasons for this change may be that i) the kinematically hot component(s) have no radial gradients, ii) the fraction of kinematically hot and metal poor component(s) increase and become significant at $Z_\mathrm{max}$$>0.6$ kpc, generating a zero or positive gradient, or iii) both explanations are valid.
The information available up to this point does not allow us to decide whether the change of gradient at $Z_\mathrm{max}$$\sim0.6$ kpc indicates the transition between one population (thin disc?) to another (thick disc?) or whether it is a mere artefact. If it were caused by a change from a population with a gradient (thin disc) to another one without a gradient (old thin disc or thick disc?) we would also expect a change in mean abundances at the same $Z_\mathrm{max}$. Figure  shows that the abundances gently diminish up to $Z_\mathrm{max}$$\sim1.5$ kpc from the plane. Then, the slope becomes steeper until $Z_\mathrm{max}$$\sim1.8$ kpc, where the abundances become constant and where one would expect the thick disc to become dominant (but be aware of the poor statistic between 1.5 and 1.8 kpc). The GCS sample (Figures  and ) exhibits a similar pattern but at different $Z_\mathrm{max}$: the rise to a flatter gradient occurs between 0.5 and 0.6 kpc and a change of slope in iron abundance at $Z_\mathrm{max}$$>$ 2.0 kpc. The change of slope at $Z_\mathrm{max}$$\sim0.6$ kpc for the gradients and at $Z_\mathrm{max}$$\sim1.5$ kpc for the mean abundance is imitated by the abundance enhancement \[X/Fe\] shown in Figures  and . The stars at $Z_\mathrm{max}$$<$`<!-- -->`{=html}0.6 kpc and $Z_\mathrm{max}$$>$`<!-- -->`{=html}1.5 kpc can be identified as thin- and thick-disc stars, respectively, while the stars in the middle region might be both the kinematically hot (old) part of the thin disc or a mix of thin- and thick-disc stars. The absence of a clear step in abundances at $Z_\mathrm{max}$$\sim0.6$ kpc excludes the transition between two populations at that height and supports the hypothesis that the change in gradients at $Z_\mathrm{max}$$\sim0.6$ kpc is caused by metal-poorer, inner disc (thick disc?) stars, which creates the fictitious flatter or positive gradient discussed in Sec. . This does not exclude the presence of an old thin disk population or the absence of a gradient in the thick-disk population reported by some authors (Allende Prieto et al. ; Bilir et al. ).

## Observed radial chemical gradients of the thin disc (the $Z_\mathrm{0.0}^\mathrm{RAVE}$ sample)

In Sec.  we discussed that a mix of different stellar populations can bias the measured gradients. In our data, the stars close to the Galactic plane are the less affected by this bias. This is confirmed in Figure , which shows that the gradients remain constant up to $Z_\mathrm{max}$$\sim0.6$ kpc because the number of thin-disc stars is larger than any other stellar populations that might be present in this region. The $Z_\mathrm{0.0}^\mathrm{RAVE}$ sample is therefore more suitable for analysing and evaluating the measured chemical gradients in the framework of the star formation history of the Galactic disc.
In the Galactic radius interval 4.5-9.5 kpc, the $Z_\mathrm{0.0}^\mathrm{RAVE}$ sample exhibits iron gradients of d\[Fe/H\]/d$R_g=-0.065$ dex kpc$^{-1}$, which lies in the range covered by many observational studies cited in Sec. . The alpha elements Mg and Si exhibit similar gradients. Aluminium has a particularly negative gradient that distinguishes it from the other elements, and we briefly discuss it below.
Under the assumption of negligible gas infall, the abundance \[X/H\] is a monotonically growing function of the number of SNe; this means that a negative abundance gradient implies a higher number of SNe that exploded in the inner disc than the outer disc in relation to the gas reservoir. According to this simple remark, chemical evolution models of the Galaxy with a process of inside-out formation of the disc were developed (Chiappini et al. , ; Cescutti el al. ; Gibson et al. ). Colavitti et al. showed that an inside-out formation of the disc is a necessary condition to explain the observed metallicity gradient. They predict negative abundance gradients that are slightly flatter than our values. For instance, Cescutti et al. found gradients of –0.052, –0.039 and –0.035 dex kpc$^{-1}$ for Fe, Mg, and Si, respectively. For comparison, ours are –0.065, –0.073 and –0.064 (see Table ). Chiappini et al. reported a slightly negative gradient in \[$\alpha$/Fe\] in the range 4-8 kpc, which turns slighly positive at 8-10 kpc. Gibson et al. predicted negative or slightly negative \[O/Fe\] gradients at the solar radius. Both works give predicted gradients close to our findings.
The model by Schönrich & Binney differs from the previous ones by construction. This model assumes that the disc forms simultaneously at all radii (no inside-out formation), with a higher infall rate at smaller radii. This causes a higher star formation rate at small Galactic radii, which generates a negative radial gradient, enhanced at late times because the inward flow of Fe-enriched gas increases the iron abundance (and decreases the ratio \[$\alpha$/Fe\]). In addition to these mechanisms, radial migration constantly acts as flattening agent. This model predicts an iron gradient of $\simeq$–0.11dex/kpc and a positive gradient for \[$\alpha$/Fe\]. This last has the opposite sign to the sign predicted by the inside-out formation models and is also opposite to our observational data, which render a shallow \[Mg/Fe\] and \[Si/Fe\] gradient (–0.01 and 0.00, respectively). This seems to favour the inside-out formation scenario, but it does not exclude the flattening action of radial migration on the gradients.
The comparison with other recent chemo-dynamical models are not conclusive. The predicted metallicity gradients range from –0.02 dex kpc$^{-1}$ (Kobayashi & Nakasato ) to –0.06 dex kpc$^{-1}$ (Minchev et al. ). This last gives a Fe gradient that closer agrees with our data, whereas the comparison with the \[Mg/Fe\] and \[Si/Fe\] gradient is not possible. Indeed, it is not correct to directly compare the gradients obtained in the present work to the ones in Figure 5 of Minchev et al. The reason is that in the RAVE analysis our sample was subdivided in three subsamples as a function of Zmax, hence each of them composed of stars with a different age mix. As an example, although in the Figure 5 (bottom) of Minchev et al. the predicted \[O/Fe\] gradients are mostly positive, to compute the gradient for a sample composed by stars with different ages one needs to take in account the density of stars of different age bins. A detailed comparison of the \[Mg/Fe\] and \[Si/Fe\] gradients with the predictions of the chemodynamical models of Minchev et al. will be carried out in a forthcoming paper.
The aluminium abundance gradient deserves a separate discussion. It is much more negative ($\sim-0.85$ dex kpc$^{-1}$) than the other elements. The aluminium abundance suffers from non-LTE effects, which depend on the stellar parameters and metallicity (Baumüller & Gehren ). Although the RAVE chemical abundances are based on LTE spectral analysis, non-LTE effects are expected to cause mere systematic errors in our Al abundance estimates based on the Al doublet at $\sim$`<!-- -->`{=html}8773Å of metal-rich (\[m/H\]$>-1.0$ dex) dwarf stars. This has no effects on the gradient determination. Aluminium is also reported to be a useful element to distinguish different stellar populations, such as disc and halo stars (Gehren et al. ). Its peculiar gradient may serve as discriminant between different chemical model of the Galaxy. Unfortunately, these element have not had enough attention in the cited chemical evolution studies, and a direct comparison with our findings is not possible.

# Conclusions

We have analysed 19 962 dwarf stars selected from the RAVE chemical catalogue and 10 616 dwarf stars selected from the Geneva-Copenhagen survey (GCS) making use of chemical, kinematic, and Galactic orbit parameters such as chemical abundances, guiding radii ($R_\mathrm{g}$), and maximum distance from the Galactic plane reached by the stars along their orbits ($Z_\mathrm{max}$). The same analyses were also applied to a RAVE mock sample created using the stellar population synthesis code GALAXIA (Sharma et al. ), which is based on the Besançon model (Robin et al. ). Using the RAVE and GCS sample, we found that the Galactic disc has a negative iron abundance gradient of $-0.065$ dex kpc$^{-1}$ and $-0.043$ dex kpc$^{-1}$ for the RAVE and GCS samples, respectively. We found that for the GCS sample the gradient is more negative ($-0.054$ dex kpc$^{-1}$) when stars at small $R_\mathrm{g}$ ($<5.5$ kpc) are excluded from the measurement. This is due to an observational bias that affects volume-limited samples such as the GCS and RAVE: stars from the inner disc can reach the Sun's neighbourhood only if they move on eccentric orbits. Such kinematically hot stars are on average more metal-poor than stars with the same $R_\mathrm{g}$, but are on more circular orbits, which will be missing from our sample. This creates a bias against high-metallicity and low-$R_\mathrm{g}$ stars, which contributes to make the gradient flatter than it really is. This effect is evident in the GCS sample and barely visible in the RAVE one. To assess the impact of this bias, we measured the gradients with different cuts at low and high $R_\mathrm{g}$, with differences in gradients no larger than 0.01 dex kpc$^{-1}$. This uncertainty in the chemical gradients is small enough not to affect our data interpretation. The presence of inhomogeneities in the volumes probed by the RAVE and GCS samples (such as disrupted open clusters and moving groups) cannot be excluded to explain the different obtained gradients. Despite the differences just highlighted, for both RAVE and GCS samples the gradient shows a common trend: it becomes flatter at increasing $Z_\mathrm{max}$ (as found by previous studies, e.g. Cheng et al. ; Boeche ; Pasquali & Perinotto ). The RAVE sample also allows gradient measurements for the elements Mg, Al, and Si, for which we found the same trend.

The GALAXIA/Besonçon mock sample shows striking differences compared with the RAVE and GCS samples, because the iron gradients are positive at any $Z_\mathrm{max}$. This surprising and unrealistic result is caused by the lack of correlation between metallicity and tangential velocity (or angular momentum) in the Besançon model. In fact, although the Besançon model assumes a radial gradient of –0.07 dex kpc$^{-1}$ for the thin disc (and zero for the other populations), the metallicities are assigned regardless of the kinematics. This means that stars with small $R_\mathrm{g}$ have the same probability to be metal-rich as stars with large $R_\mathrm{g}$, while in the real Galaxy stars coming from the inner disc are more likely to be metal-rich than stars coming from the outer disc. The absence of a correlation between chemistry and kinematics in the Besançon model removes the metallicity gradient in the mock sample when the kinematically dependent $R_\mathrm{g}$ quantity is used. Moreover, the thick-disc stars, with a mean metallicity of $-0.78$ dex and stronger asymmetric drift than the thin-disc stars, shift to lower $R_\mathrm{g}$ and \[Fe/H\], creating a positive metallicity gradient in the full sample even when the individual components have a negative or no radial gradient. The mock sample also shows an excess of thick disc stars compared with the RAVE (and GCS) sample, with a mean metallicity peak at $-0.78$ dex. This appears too low compared with the RAVE and GCS samples, where a value of $-0.5$ dex seems more realistic (also suggested by ; ; ; ). The discrepancies between the observational samples and the model can be reduced by i) decreasing the local density, ii) decreasing the vertical velocity dispersion, and iii) raising the mean metallicity of the thick disc.
Our analysis of the gradients as function of $Z_\mathrm{max}$ with a moving box car revealed that the abundance gradients remain approximately constant up to $\sim$`<!-- -->`{=html}0.6 kpc, and then increase to zero. The change in gradients at $\sim$`<!-- -->`{=html}0.6 kpc does not correspond to a change in average abundance, suggesting that the gradients become flatter because the fraction of kinematically hot, metal-poor (thick-disc) stars become significant at that height and leads to a flatter gradient, as discussed before. Another explanation is that this population of stars has no radial gradient. Both can be true, and our analysis is not able to distinguish between these possibilities. The constancy of the radial gradients up to $Z_\mathrm{max}$$\sim$`<!-- -->`{=html}0.6 kpc reveals that a stellar population with gradient zero (probably the thick disc) does not affect the thin disc properties at these heights and the values obtained here are reliable. For a sample of RAVE stars at $Z_\mathrm{max}$$<$`<!-- -->`{=html}0.4 kpc we measured gradients that are steeper (d\[Fe/H\]/d$R_\mathrm{g}$=–0.065 dex kpc$^{-1}$) than predicted by some chemical models (for instance –0.040 dex kpc$^{-1}$ by Chiappini et al. , –0.052 dex kpc$^{-1}$ by Cescutti et al. ) and flatter than others (–0.11 dex kpc$^{-1}$ by Schönrich & Binney ). We also found slightly negative (or zero) gradients for abundances relative to iron (i.e. d\[Mg/Fe\]/d$R_\mathrm{g}$=–0.009 dex kpc$^{-1}$), which supports the models of Chiappini et al. and Gibson et al. , and it is in contrast with the model by Schönrich & Binney which predict positive gradients.
The last model considers radial migration to be an important flattening agent for the radial gradients. Although it is common sense to expect that radial mixing is acting in our Galaxy, it is difficult to estimate the impact of its action on the observed chemical gradients, because migrated stars changed their original $R_\mathrm{g}$ (Roškar et al. ; Sánchez-Blázquez et al. ) and have become kinematically indistinguishable from the local sample. Detailed chemical analyses are therefore necessary to distinguish migrated from the locally born stars. High-resolution and high signal-to-noise spectroscopic surveys such as the Gaia-ESO survey , the GALactic Archaeology with HERMES (GALAH) survey , and the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey may be able to identify and quantify migrated stars.

This paper will be followed by a second one in which we measure the radial chemical gradients by using RAVE giant stars. Because RAVE is a magnitude-limited survey, and because giants stars are luminous objects, we are able to cover a larger volume, use the actual Galactocentric distance $R$ in addition to the guiding radius $R_\mathrm{g}$, and probe the Galactic disc at larger heights. This analysis will provide additional constraints on the chemical and kinematical characteristics of the Milky Way.

[^1]: An exception to this occurs in the last stages of the stellar evolution, when the chemical abundances of the stellar atmosphere can be changed by dredge-ups, which bring the products of the nucleosynthesis to the surface.

[^2]: The $frac$ parameter described in Boeche et al. gives the fraction of pixels that are non-defective and represents the goodness of the continuum fitting.

[^3]: The statistics on the celestial sphere were derived with the HEALPix package .

[^4]: We recall that in RAVE the $\alpha$-enhancement is computed as \[$\alpha$/Fe\]=$\frac{[Mg/Fe]+[Si/Fe]}{2}$.


abstract: In this paper, we study the variations of group galaxy properties according to the assembly history in SDSS-DR6 selected groups. Using mock SDSS group catalogues, we find two suitable indicators of group formation time: i) the isolation of the group, defined as the distance to the nearest neighbor in terms of its virial radius, and ii) the concentration, measured as the group inner density calculated using the fifth nearest bright galaxy to the group centre. Groups within narrow ranges of mass in the mock catalogue show increasing group age with isolation and concentration. However, in the observational data the stellar age, as indicated by the spectral type, only shows a correlation with concentration.
  .
  We study groups of similar mass and different assembly history, finding important differences in their galaxy population. Particularly, in high mass SDSS groups, the number of members, mass-to-light ratios, red galaxy fractions and the magnitude difference between the brightest and second brightest group galaxies, show different trends as a function of isolation and concentration, even when it is expected that the latter two quantities correlate with group age. Conversely, low mass SDSS groups appear to be less sensitive to their assembly history.
  .
  The correlations detected in the SDSS are not consistent with the trends measured in the mock catalogues. However, discrepancies can be explained in terms of the disagreement found in the age-isolation trends, suggesting that the model might be overestimating the effects of environment. We discuss how the modeling of the cold gas in satellite galaxies could be responsible for this problem. These results can be used to improve our understanding of the evolution of galaxies in high-density environments.
author: Tatiana Zapata$^{1,2}$[^1], Josefa Perez$^{3}$, Nelson Padilla$^{1}$, Patricia Tissera$^{3}$
$^{1}$Departamento de Astronomı́a y Astrofı́sica, Pontificia Universidad Católica de Chile, Santiago, Chile.
$^{2}$Instituto de Astronomia, Geofı́sica e Ciências Atmosféricas, Universidade de São Paulo, Brasil.
$^{3}$Instituto de Astronomı́a y Fı́sica del Espacio,Conicet-UBA, CC67, Suc.28,Ciudad de Buenos Aires, Argentina.
date: Accepted ????. Received ????; in original form ????
title: The influence of halo assembly on galaxies and galaxy groups

# Introduction

Up to only a few years back the mass of a dark-matter halo was the only parameter considered to affect its clustering properties in an important or detectable way (see for instance Padilla & Baugh 2002; Sheth, Mo & Tormen 1998; Mo & White 1996) as well as the galaxy population inside haloes. The latter has been studied in the framework of the Halo Model (Cooray & Sheth 2002; Cooray, 2005, 2006, and references therein), whereby the number of galaxies per dark-matter halo, as well as central and satellite luminosities and possibly colours, are proposed to depend only on the host-halo mass (Wang et al., 2008); by assuming this, most of the global fundamental statistics of the galaxy population can be then reproduced.

Recent studies on the clustering of dark-matter haloes found that the age and assembly history of haloes of similar mass influence their clustering amplitude. The dependence on halo age was first reported by Gao, Springel & White (2005), and has been studied into more detail finding a general dependency on the way in which a halo is assembled by Gao & White (2007), Jing, Suto & Mo (2007), Croton, Gao & White (2007), and Wechsler et al., (2007). The assembly bias has also been detected in observations by Wang et al. (2008); their studies only concentrate on the clustering amplitude of different samples selected according to group colour which they associate with group age.

Our aim is to detect the effects of halo assembly on the galaxy population. Note that even though there have been several attempts to study the environmental dependence of galaxy properties in the Sloan Digital Sky Survey (as for instance in Blanton & Berlind, 2007, and O'Mill, Padilla & Lambas, 2008), and other galaxy catalogues (in the 2-degree Field Galaxy Redshift Survey by González et al., 2005), none of these works relate their results to the influence of halo assembly. In addition to using group colour as an indicator of halo age we will explore other parameters as indicators such as the degree of dynamical relaxation of a group given by its apparent projected shape (rounder objects tend to be older, see Paz et al., 2006) or concentration (Wechsler et al. 2007, Gao & White 2007), luminosity difference between the brightest and second brightest galaxies in the group (used to identify old, fossil groups, Mendes de Oliveira & Carrasco 2007; von Benda-Beckmann et al. 2007), and group isolation.

This paper is organised as follows, Section 2 presents the sets of observational data used in this work as well as the mock catalogues constructed from semi-analytic galaxies; in Section 3 we explore different indicators of group age, Section 4 shows the variations of galaxy properties as a function of group isolation and concentration. Section 5 presents a discussion and Section 6 summarises the main results of this work.

# Data

## Observations

We use the Sloan Digital Sky Survey-Data Release 6 (SDSS-DR6, Adelman-McCarthy et al., 2007) to identify groups following the prescription presented by Merchán & Zandivarez (2005), which consists on a friends-of-friends algorithm with varying projected linking length $\sigma$, with $\sigma_0=0.239$ h$^{-1}$Mpc and fixed radial linking length $\Delta v=450 \ {\rm km s^{-1}}$. These parameters correspond to the values found by Merchán & Zandivarez to produce a reasonably complete sample with low contamination ($95\%$ and $\la 8\%$, respectively). The minimum number of galaxies per group in our sample is set to $10$.

Even though there are previous compilations of groups in the SDSS such as the catalogues compiled from the SDSS-DR3 (Merchán & Zandivarez, 2005) and DR4 (Yang et al., 2007), and by using percolation algorithms (Berlind et al., 2006) on the NYU-VAGC (Blanton et al. 2005), we need to perform our own identification in order to control the measurements of group centre, mass and membership. Furthermore, given the low statistical significance and high susceptibility to observational biases, we also need mock catalogues on which to perform our analyses to test for these effects. Therefore we apply exactly the same group detection algorithm to a SDSS mock catalogue which we will describe in the following subsection.

We use the virial theorem to compute the virial mass of groups which is given by $$M_{vir}=\frac{3\sigma_{v}^2 R_{vir}}{G}$$ where $\sigma_{v}$ is the line-of-sight velocity dispersion and $R_{vir}$ is estimated as in Merchán & Zandivarez, (2005), $$R_{vir}=\frac{\pi}{2} \frac{N_{g}(N_{g}-1)}{\Sigma_{i>j}R_{ij}^{-1}}$$ where $N_{g}$ is the number of galaxy members and $R_{ij}$ are the galaxy relative projected distances.

Throughout this paper the quoted errors in measurements involving SDSS groups will correspond to errors in the mean obtained using the jacknife technique with $10$ subsamples.

## Mock catalogues

We construct a mock SDSS catalogue using the semi-analytic model (SAM) from De Lucia et al. (2006), who use the merger trees from the Millennium Simulation (Springel et al., 2005) to follow different physical processes that lead to the final $z=0$ galaxy population in a $\Lambda$CDM cosmology, with cosmological parameters determined from the combined analysis of the 2-degree Field Galaxy Redshift Survey (2dFGRS, Colless et al. 2001) and the first year Wilkinson Microwave Anisotropy Probe (WMAP) data (Spergel et al. 2003): $\Omega_{m}=0.25$, $\Omega_{b}=0.045$, $\Omega_{\Lambda}=0.75$, $h=0.73$, $n=1$ and $\sigma_{8}=0.9$. The Millennium simulation consists of $2160^3$ particles, in a box of $500$h$^{-1}$Mpc on a side, with a particle mass resolution of $8.6\times 10^8$h$^{-1}$M$_{\odot}$. The resulting galaxies conform a complete sample down to a magnitude $M_r=-16$.

In order to mimic the SDSS, we place four observers on different locations in the $z=0$ output from the numerical simulation so that the volumes surveyed by each observer are as independent as possible; this is achieved by placing observers in the centre of four opposite octants in a diagonal direction, with two of the observers (on a long diagonal) facing in a direction offset by 90 degrees (in the $x-y$ and $y-z$ planes) with respect to the other two observers. The minimum distance between observers is $\simeq 350$h$^{-1}$Mpc, and they sample different regions of the simulation in different directions. We apply the same angular mask affecting the spectroscopic main sample, and apply a magnitude limit cut of $r_{\rm lim}=17.77$. For each galaxy we store its observed properties, such as redshift (which includes the peculiar motion), angular position in the sky, and apparent magnitudes in several bands. However, the main advantage of the mock catalogue is that we also store several intrinsic properties such as luminosity in different bands and properties of the host halo including its mass.

We apply the adaptive FOF algorithm to each mock SDSS catalogue, using exactly the same search parameters applied to the real data. In this way, we obtain four mock group catalogues which are studied the same way as the real SDSS groups; in particular, we use the same minimum number of galaxies per group, as well as the same method to estimate group masses.

In the remainder of this paper the results obtained from the mocks will correspond to the average over our four mock SDSS catalogues; the uncertainties represented by the errorbars will correspond to the errors in the mean obtained using the jacknife technique with $10$ subsamples, averaged over the mock catalogues.

# Searching for observational indicators of group assembly

In this section, we analise the variation of galaxy group properties with the aim of finding an indicator of group age which can be related to the assembly history of the group host halo. We will analyse several candidates of observational parameters that will allow us to rank our observational and simulated group samples according to their assembly history.

In the case of the mock catalogues, the group age is directly estimated by the mean stellar-mass weighted age of the galaxies hosted by a given group. The star formation histories of galaxies in groups are determined by the rate of gas consumption which is modulated by the assembly history of the haloes (i.e. collapse, infall, mergers) and feedback mechanisms which in effect are regulated by the size of the haloes. A recent work by Li, Mo & Gao (2008) shows that a halo age related to the formation history of its member galaxies is the most appropriate for detecting the assembly effects on the bias parameter; in this case the effect is seen for all halo masses, including the ranges we will study in this paper, $12.9<\log_{10}(M_{\rm vir}$ h$^{-1}$M$_{\odot})<14.0$. This justifies our choice of halo age indicator.

Observationally, this parameter is more complicated to estimate. In order to find the best observational indicator which could be related to group age, we first search for possible correlations between the mock group ages and other group parameters.

Although the first candidate for age that comes to mind is the group mass as larger groups started to collapse earlier (see for instance Paz et al. 2006), our aim is to explore any possible variation of group properties with age for fixed values of group mass. Hence, for fixed group masses, we test a series of observable parameters:

-   the projected group shape measured by the ratio between the minor and major semi-axes of an ellipse such that the eigenvalues of the anisotropic part of its inertia tensor are the same as obtained from the projected positions of the group members,

-   the difference in magnitudes between the first and second brightest galaxies in the group ($\Delta M_{12}$),

-   the fraction of red galaxies in the group ($u-r \ge 2.0$),

-   the degree of isolation of a group defined as the distance $D_{\rm neighbor}/R_{\rm vir}$ (measured in terms of the group virial radius) to the nearest group neighbour,

-   the concentration measured by the group inner density ($\Sigma_5$), calculated by using the projected distance to the 5$^{th}$ galaxy, brighter than M$_{\rm r} = -20.5$, and

-   the mass-to-light ratio which we obtain using luminosities measured in the r-band.

We find that the group shape is not a good indicator of group age due to the lack of statistical significance of a relation between them. Conversely, the fraction of red galaxies, the $M/L$ ratio, $D_{\rm neighbor}/R_{\rm vir}$ and $\Sigma_5$ show clear correlations with mean stellar-mass weighted ages as it can be appreciated in Fig. . In all cases we divide the group catalogues into subsamples defined by two narrow (0.4 dex wide) ranges of halo mass: high mass groups selected within $13.6<\log_{10}(M_{\rm vir}$h$^{-1}$M$_{\odot})<14.0$ and low mass groups, within $12.9<\log_{10}(M_{\rm vir}$h$^{-1}$M$_{\odot})<13.3$ (from this point on we will use the same halo mass ranges to study galaxy group properties as a function of different age indicators).

Even though there is also a clear trend of age with the fraction of red galaxies and the mass-to-light ratio, it has been shown that semi-analytic models overestimate the fraction of red galaxies (e.g. Weinmann et al. 2006). Conversely, the trends with group concentration and isolation are very well defined, showing that for the model, the most isolated and concentrated groups tend to be older. Since these parameters are more directly related to the structure of the dark matter, we adopt both estimators as good indicators of group age and assembly history. As it can be seen from the lower panels of Fig. , a semi-analytic galaxy group with no neighbour groups out to ten times its virial radius, or with an internal density larger than $80 \  h^2 {\rm Mpc^{-2}}$, is likely to be up to $800$ Myr older than a group living in a more populated area or with a lower concentration.

This analysis also shows that both group age indicators (isolation and concentration) depend on halo mass. In fact, in Fig.  we plot halo mass as a function of concentration (left panel) and isolation (right panel) for the mock (solid symbols) and the SDSS (open symbols). Both the observations and model show a trend with both parameters, however the dependence of mass on concentration is very weak, indicating that the concentration may be better suited to study assembly effects without introducing biases. Nevertheless, given that the isolation of a group can store information on its large-scale neighbourhood, we will also include it in our analysis.

In order to detect dependencies of different group properties on these two age indicators, we will separate groups according to the proposed age parameter in two bins corresponding to young and old objects. Then for each bin, we calculate the mean group properties and errors using the jacknife technique, and characterise the statistical significance of a trend by quoting the resulting difference between the galaxy population property in old and young groups ($\chi^2$). Notice that for one degree of freedom, $\chi^2=n$ indicates a $n\sigma$ detection. In the case of concentration, we will consider as old objects those characterised by $\Sigma_5\geq70$h$^2$Mpc$^{-2}$, and as young objects those with $\Sigma_5\leq30$h$^2$Mpc$^{-2}$. The limits for the case of isolation are set to $D_{\rm neighbor}/R_{vir}\geq8$ for old objects and $D_{\rm neighbor}/R_{vir}\leq7$ for young objects. One important advantage of this approach is that correlations between different isolation or concentration bins are reduced to a minimum.

A possible equivalent observational estimator of stellar mass age is the galaxy spectral type (Yip et al. 2004); lower values of spectral type correspond to earlier galaxy types and, therefore, to older stellar populations. In a similar way, spectral types correlate with colours, corresponding to a value of -0.08 for the limit spectral type which separates the red and blue sequences. In order to test the validity of the previous theoretical indicators of age, Fig. shows the mean spectral type of SDSS groups as a function of concentration (left panel) and isolation (right panel). For massive groups, we find different trends with these parameters. While higher concentration massive SDSS groups are populated with systematically lower spectral type (i.e., older stellar populations, $4.85\sigma$ effect, see Table 1), the more isolated massive groups contain somewhat higher spectral type galaxies (i.e., younger, $1.39\sigma$ effect, see Table 1). These correlations show that SDSS groups confirm the predicted correlation between concentration and age found in mock catalogues (Fig. ). However a disagreement in the SDSS and mock isolation-age/spectral type trends reveals that the model may be overestimating the effects of environment. We will come back to this point in the following sections. On the other hand, low mass SDSS groups show no significant trend (at most $1\sigma$). with either isolation or concentration, again in contradiction to the findings obtained from the mock catalogues.

In the following section, we turn to the study of the variation of galaxy group properties as a function of assembly epoch as indicated by group concentration and isolation. When using isolation as an age indicator we bear in mind the differences between observations and models.

# galaxy group properties as a function of group age

```latex
\begin{table*}
 \centering
 \begin{minipage}{160mm}
  \caption{SDSS results: significance of trends in spectral type (column 3), mass-to-light
ratios (column 4), occupation numbers (column 5), fraction of red galaxies (column 6), and magnitude difference
between brightest and second brightest group member (column 7), as a function of concentration (top rows) and isolation
(bottom rows).  All quantities are rounded to $\sim1$ significant figure in the error, estimated using the jacknife technique.}
   \label{table1}
  \begin{tabular}{@{}ccccccc@{}}
  \hline
  Mass range & quantity & Spec. Type & $M_{vir}/L_r/(M_{\odot}$h$^{-1}L_{\odot}^{-1})$ & $Log_{10}(N_{gal})$ & fr. red gals. & $\Delta M_{12}$\\
  \hline
  high mass& low $\Sigma_5$                      &$-0.097\pm0.002$  &$572\pm13$    & $1.178\pm0.032$ &$0.784\pm0.004$  &$0.49\pm0.04$\\
           & high $\Sigma_5$                     &$-0.108\pm0.001$  &$443\pm8$     & $1.454\pm0.053$ &$0.809\pm0.005$  &$0.54\pm0.04$\\
           & $\chi^2$ (trend)                    &$4.85$ (decr.)    &$8.37$ (decr.)& $4.47$ (incr.)  &$4.02$ (incr.)   &$0.9$ (none)\\
  low mass & low $\Sigma_5$                      &$-0.08\pm0.002$   &$180\pm3$     & $1.065\pm0.017$ &$0.77\pm0.01$    &$0.49\pm0.06$\\
           & high $\Sigma_5$                     &$-0.077\pm0.002$  &$158\pm5$     & $1.064\pm0.013$ &$0.815\pm0.006$  &$0.64\pm0.04$\\
           & $\chi^2$ (trend)                    &$1.07$ (incr.)    &$3.94$ (decr.)& $0.006$(none)   &$3.84$ (incr.)   &$1.99$ (incr.)\\
  \hline
  high mass& low $D_{\rm neighbor}/R_{\rm vir}$  &$-0.099\pm0.001$  &$498\pm13$    &$1.302\pm0.021$ &$0.809\pm0.005$   &$0.45\pm0.01$\\
           & high $D_{\rm neighbor}/R_{\rm vir}$ &$-0.091\pm0.006$  &$704\pm40$    &$1.121\pm0.038$ &$0.763\pm0.016$   &$0.51\pm0.05$\\
           & $\chi^2$ (trend)                    &$1.39$ (incr.)    &$4.92$ (incr.)&$7.9$ (decr.)   &$2.7$ (decr.)     &$1.18$ (incr.)\\
  low mass & low $D_{\rm neighbor}/R_{\rm vir}$  &$-0.0735\pm0.001$ &$227\pm6$     &$1.086\pm0.014$ &$0.797\pm0.004$   &$0.44\pm0.02$\\
           & high $D_{\rm neighbor}/R_{\rm vir}$ &$-0.0733\pm0.003$ &$144\pm3$     &$1.069\pm0.027$ &$0.749\pm0.007$   &$0.52\pm0.02$\\
           & $\chi^2$ (trend)                    &$0.06$ (none)     &$11.51$ (decr.)&$8.02$ (decr.)  &$5.33$ (decr.)   &$2.53$ (incr.)\\
  \hline
  \end{tabular}
  \end{minipage}
\end{table*}
```
```latex
\begin{table*}
 \centering
 \begin{minipage}{160mm}
  \caption{Mock results: columns and rows are as in Table 1.
}
   \label{table2}
  \begin{tabular}{@{}ccccccc@{}}
  \hline
  Mass range & quantity & Spec. Type & $M_{vir}/L_r/(M_{\odot}$h$^{-1}L_{\odot}^{-1})$ & $Log_{10}(N_{gal})$ & fr. red gals. & $\Delta M_{12}$\\
  \hline
  high mass& low $\Sigma_5$                      &N/A &$202\pm2$     & $1.231\pm0.006$ &$0.863\pm0.002$  &$0.70\pm0.01$\\
           & high $\Sigma_5$                     &N/A &$177\pm2$     & $1.379\pm0.013$ &$0.913\pm0.002$  &$0.77\pm0.03$\\
           & $\chi^2$ (trend)                    &N/A &$10.3$ (decr.)& $10.0$ (incr.)  &$17.6$ (incr.)   &$2.23$ (incr.)\\
  low mass & low $\Sigma_5$                      &N/A &$60\pm1.6$    & $1.088\pm0.008$ &$0.854\pm0.005$  &$0.61\pm0.02$\\
           & high $\Sigma_5$                     &N/A &$56\pm1.7$    & $1.087\pm0.014$ &$0.938\pm0.006$  &$0.67\pm0.03$\\
           & $\chi^2$ (trend)                    &N/A &$1.9$ (decr.) & $0.09$ (none)   &$11.4$ (incr.)   &$1.45$ (incr.)\\
  \hline
  high mass& low $D_{\rm neighbor}/R_{\rm vir}$  &N/A& $198\pm2$     & $1.313\pm0.006$ &$0.881\pm0.002$  &$0.75\pm0.01$\\
           & high $D_{\rm neighbor}/R_{\rm vir}$ &N/A& $177\pm2$     & $1.227\pm0.009$ &$0.921\pm0.002$  &$0.73\pm0.02$\\
           & $\chi^2$ (trend)                    &N/A& $7.2$ (decr.) & $8.1$ (decr.)   &$11.28$ (incr.)  &$0.86$ (none)\\
  low mass & low $D_{\rm neighbor}/R_{\rm vir}$  &N/A& $64\pm1.1$    & $1.062\pm0.006$ &$0.852\pm0.01$   &$0.64\pm0.02$\\
           & high $D_{\rm neighbor}/R_{\rm vir}$ &N/A& $53\pm0.7$    & $1.07\pm0.01$   &$0.916\pm0.002$  &$0.64\pm0.02$\\
           & $\chi^2$ (trend)                    &N/A& $8.7$ (decr.) & $0.67$ (none)   &$8.4$ (incr.)    &$0.08$ (none)\\
  \hline
  \end{tabular}
  \end{minipage}
\end{table*}
```

In principle, halo models (cf. Cooray & Sheth, 2002) and the conditional luminosity function (Cooray 2005, 2006) are able to establish the properties of galaxies in dark-matter haloes using only the halo mass. This technique has been very successful in reproducing several galaxy statistics such as the correlation function of SDSS galaxies (Zehavi et al., 2004) and numerical simulations (Zheng et al., 2005), and the evolution of the luminosity function of galaxies up to $z=6$ (Cooray 2005). However, as we will show in this section, the properties of galaxies in groups of similar masses are also dependent on the concentration and isolation of the group, parameters related to the assembly history of the system.

We start analysing the variations of mass-to-light ratios, $M/L$, of low and high mass groups for different degrees of isolation and concentration. Fig.   shows the $M/L$ ratios as a function of group concentration (left panel) and isolation (right panel). In consistency with our previous result for spectral types, high mass SDSS groups show two clear opposite trends of decreasing ($8.37\sigma$, Table 1) and increasing $M/L$ ($4.92\sigma$, Table 1) with concentration and isolation, respectively, of almost a $50\%$ change. As can be seen, both trends are comparatively significant. A more detailed analysis shows that $(90\pm10)\%$ of the $M/L$ ratios exhibited by groups isolated by more than $8$ virial radii lies above the median relation between $M/L$ ratio to $\Sigma_5$; however, the values of $\Sigma_5$ for these isolated groups span almost the entire range of concentrations shown in the figure. This tells us that more isolated groups at a given concentration show higher $M/L$ ratios. Note also that since the analysis is carried out in a narrow mass range, the increase or decrease in $M/L$ can be directly related to a decrease or increase, respectively in the r-band luminosity. Therefore, the opposite behaviours in the $M/L$ ratios observed for the SDSS high mass groups with concentration and isolation can also be associated to opposing trends in their r-band luminosities. This can be explained in terms of their spectral types; while high mass SDSS groups tend to be older with increasing concentration (and consequently brighter in the r-band), they tend to be younger with increasing isolation (and fainter in the same band).

On the other hand, Fig.   also shows that low mass SDSS groups exhibit compatible trends of decreasing $M/L$ with both concentration and isolation, with a high statistical significance (Table 1). Interestingly, groups in the mock catalogues show this trend for both high and low mass groups (the mass-to-light ratios decreases as a function of isolation and concentration). This discrepancy between observational and mock catalogues for high mass groups is an indication of possible problems in the modeling of the physical processes associated to the baryonic component in the SAM, since the mass-to-light ratio, $M/L$, is strongly dependent on the galaxy population.

In order to investigate further the origin of these behaviours, we also estimate the number of group members as a function of concentration and isolation, for both SDSS and mock groups. As it can be seen in Fig. , the number of members for high mass SDSS groups decreases with isolation ($7.9\sigma$, Table 1) and increases with concentration ($4.47\sigma$, Table 1). This result is also consistent with the $M/L$ and r-band luminosity tendencies since a decreasing number of members produces lower total luminosities as the group isolation increases, and an increasing number of galaxy members results in higher luminosities as the concentration increases. Note that, although high mass groups in the mock catalogues exhibit similar trends than the observations (see Table 2) for the number of galaxy members as a function of concentration and isolation (Fig ), they totally disagree in reproducing the $M/L$ behavior. This suggests that the intrinsic luminosity of the galaxy members vary differently in the model and observations as a function of isolation and concentration. This discrepancy between mock and SDSS groups may be an indication that the model fails in describing the physical processes associated to the baryonic component to some level. Several facts could be involved: there could be an offset in the gas accretion rates onto the halo, the star formation timescales could be too short, the threshold cold gas surface density for star formation could be too low, the suppression of star formation by AGN heating could be overestimated, the merging timescales could be too short, environmental effects could be exacerbated, and so on. Because the most significant discrepancy is observed for the most massive haloes in low density environments (large isolation), we can rule out some of these sources of disagreement. In the perspective of the current hierarchical model of structure formation, where massive haloes assemble rather late through the merging of smaller systems, gas cooling processes are the main contributors to the galaxy mass-growth in low mass haloes. On the other hand, galaxy-galaxy mergers are more efficient within small haloes with low velocity dispersion. As a consequence, both processes would not explain the observed disagreement. AGN feedback acts regulating the star formation activity in massive galaxies and has been proved to be reasonably well suited to reproduce several observational statistics, so it is not likely to be responsible for the present problem. Thus, the remaining possible causes could lie in the star formation prescription or in environmental processes. Recently, analytic and hydrodynamic simulations have shown that the environment in groups and cluster is less aggressive than previously assumed. It has been argued that the excess of red satellite owes to an oversimplified treatment of ram pressure stripping of their hot gaseous halos. In general, SAMs assume that the entire hot gas reservoir of a galaxy is instantaneously stripped at the moment of accretion into a larger halo (when it becomes a satellite galaxy). Recent improvements in SAMs allow a better agreement with observations for the fractions of red satellites once the following considerations are implemented; (i) a decreased stripping efficiency, and (ii) a significant loss of low-mass satellites tidally disrupted during the accretion by a central host galaxy (Font et al 2008; Kang et al. 2008). Thus, satellite galaxies are able to retain a significant fraction of their hot gas for several Gyr., replenishing their reservoirs of cold, star forming gas and remaining blue for relatively longer periods of time.

Given that the fraction of red galaxies is a good indicator of the evolutionary state of a group (Bell et al. 2004), we also investigate the variations of this parameter with concentration and isolation in the mock and SDSS groups (Fig.). We find that both SDSS and mock groups show a similar trend of increasing red galaxy fraction as a function of concentration (left panel of the figure, significance levels greater than $4\sigma$, Tables 1 and 2). This result is consistent with the expected behaviour since concentration correlates with age for both SDSS and mock groups. However, it is noticeable the different fractions present in the model and in the SDSS; the former are $\sim 10\%$ larger, in agreement with previous studies that indicate that semi-analytic colours in the Croton et al. (2006) model are systematically redder than SDSS galaxies (Weinmann et al. 2006)[^2]. On the other hand, as can be seen in the right panel of this figure, high mass SDSS groups show a decreasing fraction of red galaxies as groups become more isolated, in agreement with the larger fraction of late type galaxies found to populate them (Fig. ). We calculate how often extremely isolated SDSS massive groups with distances to their closest neighbors larger than $>8R_{\rm vir}$, show red galaxy fractions below the median as a function of $\Sigma_5$, and find that this occurs with a frequency of a $(70\pm17)\%$, again with $\Sigma_5$ values over the whole range shown in the figure. This is in agreement with the results found for the $M/L$ dependence on isolation and concentration.

In contrast, as we have already discussed, mock catalogues show that high mass groups have a larger fraction of red galaxies as the isolation increases. In fact, this disagreement suggests that the semi-analytic model fails to reproduce the evolution of galaxy population in the lowest density environments, overestimating the level of aging. This suggests that the model should include some physical process to stop the reddening of the group/cluster galaxy population in more isolated environments as was discussed previously.

So far, our analysis of SDSS groups proves that at a given concentration, the degree of isolation for high mass groups may play a role at not only determining the number of group members but also at modulating the star formation histories (i.e. high mass, isolated groups have a smaller fraction of red galaxies). Again, even though the mock catalogues are able to reproduce the dependence of number of members on isolation, they fail to show the observed trend for the red galaxy fraction.

With the purpose of confirming that the galaxy population is intrinsically different for groups of fixed mass and different degrees of concentration and isolation, we estimate the magnitude difference between the brightest and second brightest group galaxies, $\Delta M_{12}$. As it can be seen in Fig.  , both the observed and simulated groups show a clear albeit weak tendency for higher values of $\Delta M_{12}$ for more concentrated groups (greater than $1\sigma$ detections, Tables 1 and 2), indicative of older group ages (Mendes de Oliveira & Carrasco, 2007). In the case of the isolation, as shown in the right panel of Fig.   the dependence of $\Delta M_{12}$ is similar, although not that clear for massive SDSS groups. These findings suggest that at a given mass, groups show slightly different galaxy populations depending on their concentration and isolation. It is important to recall that for the more massive systems, concentration and isolation are not correlated. Notice the higher typical values of $\Delta M_{12}$ in the mock catalogues suggesting a larger difference between central and satellite galaxies in the model than those in the SDSS. This should be considered a possible problem in the galaxy formation modeling which suggests again an overestimation of the effects of environment in the semi-analytic scheme.

Finally, we studied the variation of group projected shapes on concentration and therefore group age. We were not able to detect any significant variations in this sense for either the mock or SDSS group samples, even though it is expected that older groups or dark-matter haloes show rounder intrinsic and apparent shapes (Paz et al., 2006). Note that the analysis was carried out using the same narrow ($0.4$ dex) ranges of group mass and different concentrations, that according to our mock catalogues, should have included groups with ages differing by up to $800$Myr.

# Discussion

The results found in the previous sections on the effects of isolation and concentration on the properties of groups and galaxy groups, can be used to shed light on several aspects of the galaxy formation process, in particular in view of several discrepancies in the behaviour of real and simulated galaxies.

We start by reviewing the results obtained for the mock catalogues. In the mocks, more isolated and concentrated groups, independently of their masses, show almost constant mass-to-light ratios, slightly higher $\Delta M_{12}$ and higher red galaxy fractions. With respect to the number of members, mock catalogues show different behaviours; low mass groups have constant galaxy numbers independently of their concentration or isolation, and massive groups show opposing trends of increasing and decreasing number of members as the concentration and isolation increase (i.e. as the age increases), respectively.

On the other hand, high mass SDSS-DR6 groups show different trends depending on their concentration and isolation. We now analyse high- and low-mass group results separately. i) For high masses, the most isolated groups are characterised by higher mass-to-light ratios, lower number of members and lower red galaxy fractions as a consequence of being more populated by late type galaxies. However, as a function of increasing concentration we find lower mass-to-light ratios and higher number of members and red galaxy fractions, in agreement with their older stellar populations as determined by their spectral types. Considering that the group age is proportional to the concentration and anti-correlates with the isolation, the discrepant trends of $M/L$ and red galaxy fractions in SDSS-DR6 groups suggest different assembly processes indicated by different degrees of concentration and isolation; in particular, the conditions placed by the large-scale structure around groups seems to be important in shaping a group galaxy population such that more isolated groups present lower numbers of galaxies, bluer colours and consequently larger mass-to-light ratios. ii) For low mass SDSS groups, trends are more difficult to detect with either concentration or isolation, except for a consistent decrease of $M/L$ with concentration and isolation (and consistent with the results from the mocks).

As can be seen, the characteristics of galaxy groups in the SAM only reproduce the observations to some level; in particular, massive groups show a dependence of $M/L$ and red galaxy fraction on isolation and concentration in stark contrast to the observational results. Also, the values of $\Delta M_{12}$ and red galaxy fractions are higher than in the SDSS regardless of group mass. These problems may be understood by considering that more isolated groups are older in the model, in contradiction with the results from the observations. This suggests an overestimation of the environmental effects on the evolution of systems, with an excess of old galaxies in the lowest density environments. Particularly, the treatment of satellites in the SAM, which lose their hot gas component and, consequently, their sources of gas accretion when entering into a new dark-matter halo, could be responsible for most of the trends observed in the mock groups, including the high values of red galaxy fraction and $\Delta M_{12}$, as well as the lower $M/L$ values.

The process thought to be responsible for the removal of cold gas from satellites is ram-pressure (Haynes & Giovanelli, 1986, Solanes et al., 2001, and more recently Brueggen & De Lucia, 2007, and references therein), which some of the current models of semi-analytic galaxy formation (Cole et al., 2000, Baugh et al., 2005, De Lucia et al., 2006, Croton et al., 2006, Bower et al., 2006, Lagos, Cora & Padilla, 2008) consider to produce an instantaneous result. Our observational results show several details that could be used to improve the modeling of processes that drive the evolution of galaxies. For instance, if satellites in isolated groups or clusters were allowed to form stars after being accreted, the fractions of red galaxies and $\Delta M_{12}$ values would decrease and the agreement with the observations would improve.

# Conclusions

We have searched for possible dependencies of the galaxy population in groups on their assembly history, focusing on possible observational indicators of group age, with the aim of applying the study to the SDSS-DR6 dataset.

In order to find indicators of group age, we use mock SDSS group catalogues and find two suitable candidates: i) the isolation of the group defined as the distance to the nearest neighbor in terms of its virial radius, and ii) the concentration, measured as the density calculated using the fifth closest, bright galaxy to the group centre. Groups within narrow ranges of group mass in the mock catalogue show an increase of group age with isolation and concentration. We also find that the group age is correlated with the luminosity and colour of group galaxies, and with the magnitude difference between the brightest and second brightest group galaxies.

In order to test if these two theoretical age indicators can also be used in SDSS groups, we analyse the dependence of group spectral types on isolation and concentration finding consistent results only for the latter

Using both, the isolation and concentration of galaxy groups as indicators of their assembly history, we study properties of group galaxies in SDSS and mock groups of similar mass. We find important variations in mass-to-light ratios, red galaxy fractions, and $\Delta M_{12}$. These correlations or relations are not compatible between model and observations, a problem that can be explained in terms of the differing trends found for age and spectral type as a function of isolation. We discuss how the modeling of the cold gas in satellite galaxies in the SAM could be responsible for these discrepancies.

In this paper we have confirmed that in addition to a difference in clustering amplitude, groups of similar mass and different assembly histories also show important differences in the characteristics of their member galaxies, an aspect that can be added to current halo model and conditional luminosity functions, and that can also be used to improve the modeling and understanding of the evolutionary processes that shape galaxies in the observed Universe.

# Acknowledgments

We thank the Referee for his/her useful comments which have contributed to improve this paper. This work was supported in part by the FONDAP "Centro de Astrofı́sica", Fundación Andes, Consejo Nacional de Ciencia y Tecnologı́a (PIP 6446) y Agencia Nacional de Promoción Cientı́fica y Técnica (PICT 32342). NP was supported by a Proyecto Fondecyt Regular No. 1071006. This work was supported by the European Commission's ALFA-II programme through its funding of the Latin-american European Network for Astrophysics and Cosmology (LENAC).

[^1]: E-mail: email@example.com

[^2]: Croton & Farrar (2008) study galaxy colours in voids in the 2-degree Field Galaxy Redshift Survey (Colles et al., 2001), and find that the semi-analytic model by Croton et al. (2006) produces similar galaxy colours in such low-density environments. Notice that this is not to be compared to our results since in our case, even in the most isolated environment galaxies are embedded in virialised groups of at least $10$ galaxy members, i.e., not in galaxy voids.