The second data release of the Gaia mission has revealed, in stellar velocity and action space, multiple ridges, the exact origin of which is still debated. Recently, we demonstrated that a large Galactic bar with pattern speed 39 km/s/kpc does create most of the observed ridges. Among those ridges, the Hercules moving group would then be associated to orbits trapped at the co-rotation resonance of the bar. Here we show that a distinctive prediction of such a model is that the angular momentum of Hercules at the Sun's radius must significantly decrease with increasing Galactocentric azimuth, i.e. when getting closer to the major axis of the bar. We show that such a dependence of the angular momentum of trapped orbits on the azimuth would on the other hand not happen close to the outer Lindblad resonance of a faster bar, unless the orbital distribution is still far from phase-mixed, namely for a bar perturbation younger than ~ 2 Gyr. Using Gaia DR2 and Bayesian distances from the StarHorse code, and tracing the average Galactocentric radial velocity as a function of angular momentum and azimuth, we show that the Hercules angular momentum changes significantly with azimuth as expected for the co-rotation resonance of a dynamically old large bar.

We test the performance of the semi-analytic self-consistent Just-Jahrei{\ss} disc model (JJ model) with the astrometric data from the Tycho-Gaia Astrometric Solution (TGAS) sub-catalogue of the first Gaia data release (Gaia DR1), as well as the radial velocities from the fifth data release of the Radial Velocity Experiment survey (RAVE DR5). We use a sample of 19,746 thin disc stars from the TGAS$\times$RAVE cross-match selected in the local solar cylinder of 300 pc radius and 1 kpc height below the Galactic plane and simulate this sample via the forward modelling technique. First, we convert the predicted vertical density laws of the thin disc populations into a mock sample. Then the obtained mock populations are reddened with a 3D dust map and are subjected to the selection criteria corresponding to the RAVE and TGAS observational limitations as well as to additional cuts applied to the data sample. We calculate the quantities of interest separately at different heights above the Galactic plane taking into account the distance error effects separately in horizontal and vertical directions. We investigate the simulated sample in terms of the vertical number density profiles, Hess diagrams and velocity distribution functions. Basing on a good agreement of our simulations with the data, we conclude that our fiducial disc model confidently reproduces the vertical trends in the thin disc stellar population properties. Thus, it can serve as a starting point for the future extension of the JJ model to other Galactocentric distances.

Context. Dynamically self-consistent galactic models are necessary for analysing and interpreting star counts, stellar density distributions, and stellar kinematics in order to understand the formation and the evolution of our Galaxy. Aims. We modify and improve the dynamical self-consistency of the Besan\c{c}on Galaxy model in the case of a stationary and axisymmetric gravitational potential. Methods. Each stellar orbit is modelled by determining a Staeckel approximate integral of motion. Generalised Shu distribution functions with three integrals of motion are used to model the stellar distribution functions. Results. This new version of the Besan\c{c}on model is compared with the previous axisymmetric BGM2014 version and we find that the two versions have similar densities for each stellar component. The dynamically self-consistency is improved and can be tested by recovering the forces and the potential through the Jeans equations applied to each stellar distribution function. Forces are recovered with an accuracy better than one per cent over most of the volume of the Galaxy.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

We explore the connections between stellar age, chemistry, and kinematics across a Galactocentric distance of $7.5 < R\,(\mathrm{kpc}) < 9.0$, using a sample of $\sim 12\,000$ intermediate-mass (FGK) turnoff stars observed with the RAdial Velocity Experiment (RAVE) survey. The kinematics of this sample are determined using radial velocity measurements from RAVE, and parallax and proper motion measurements from the Tycho-Gaia Astrometric Solution (TGAS). In addition, ages for RAVE stars are determined using a Bayesian method, taking TGAS parallaxes as a prior. We divide our sample into young ($0 < \tau < 3$ Gyr) and old ($8 < \tau < 13$ Gyr) populations, and then consider different metallicity bins for each of these age groups. We find significant differences in kinematic trends of young and old, metal-poor and metal-rich, stellar populations. In particular, we find a strong metallicity dependence in the mean Galactocentric radial velocity as a function of radius ($\partial {V_{\rm R}}/\partial R$) for young stars, with metal-rich stars having a much steeper gradient than metal-poor stars. For $\partial {V_{\phi}}/\partial R$, young, metal-rich stars significantly lag the LSR with a slightly positive gradient, while metal-poor stars show a negative gradient above the LSR. We interpret these findings as correlations between metallicity and the relative contributions of the non-axisymmetries in the Galactic gravitational potential (the spiral arms and the bar) to perturb stellar orbits.

We use data from the Radial Velocity Experiment (RAVE) and the Tycho-Gaia astrometric solution catalogue (TGAS) to compute the velocity fields yielded by the radial (VR), azimuthal (Vphi) and vertical (Vz) components of associated Galactocentric velocity. We search in particular for variation in all three velocity components with distance above and below the disc midplane, as well as how each component of Vz (line-of-sight and tangential velocity projections) modifies the obtained vertical structure. To study the dependence of velocity on proper motion and distance we use two main samples: a RAVE sample including proper motions from the Tycho-2, PPMXL and UCAC4 catalogues, and a RAVE-TGAS sample with inferred distances and proper motions from the TGAS and UCAC5 catalogues. In both samples, we identify asymmetries in VR and Vz. Below the plane we find the largest radial gradient to be dVR / dR = -7.01+- 0.61 km\s kpc, in agreement with recent studies. Above the plane we find a similar gradient with dVR / dR= -9.42+- 1.77 km\s kpc. By comparing our results with previous studies, we find that the structure in Vz is strongly dependent on the adopted proper motions. Using the Galaxia Milky Way model, we demonstrate that distance uncertainties can create artificial wave-like patterns. In contrast to previous suggestions of a breathing mode seen in RAVE data, our results support a combination of bending and breathing modes, likely generated by a combination of external or internal and external mechanisms.

We study the kinematics of stellar disc populations within the solar neighbourhood in order to find the imprints of the Galactic bar. We carry out the analysis by developing a numerical resolution of the 2D2V Collisionless Boltzmann Equation (CBE) and modelling the stellar motions within the plane of the Galaxy within the solar neighbourhood. We recover similar results to these obtained by different authors using N-body simulations, but we can also numerically identify faint structures thanks to the canceling of the Poisson noise. We find that the ratio of the bar pattern speed to the local circular frequency is in the range $\Omega_B/\Omega$ = 1.77 to 1.91. If the Galactic bar angle orientation is within the range from 24 to 45 degrees, the bar pattern speed is between 46 and 49 km/s/kpc.

We attempt to constrain the kinematics of the thin and thick disks using the Besancon population synthesis model together with RAVE DR4 and Gaia first data release (TGAS). The RAVE fields were simulated by applying a detailed target selection function and the kinematics was computed using velocity ellipsoids depending on age in order to study the secular evolution. We accounted for the asymmetric drift computed from fitting a Stackel potential to orbits. Model parameters such as velocity dispersions, mean motions, and velocity gradients were adjusted using an ABC-MCMC method. We made use of the metallicity to enhance the separation between thin and thick disks. We show that this model is able to reproduce the kinematics of the local disks in great detail. The disk follows the expected secular evolution, in very good agreement with previous studies of the thin disk. The new asymmetric drift formula, fitted to our previously described St\"ackel potential, fairly well reproduces the velocity distribution in a wide solar neighborhood. The U and W components of the solar motion determined with this method agree well with previous studies. However, we find a smaller V component than previously thought, essentially because we include the variation of the asymmetric drift with distance to the plane. The thick disk is represented by a long period of formation (at least 2 Gyr), during which, as we show, the mean velocity increases with time while the scale height and scale length decrease, very consistently with a collapse phase with conservation of angular momentum. This new Galactic dynamical model is able to reproduce the observed velocities in a wide solar neighborhood at the quality level of the TGAS-RAVE sample, allowing us to constrain the thin and thick disk dynamical evolution, as well as determining the solar motion.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

We present a novel analysis of the metal-poor star sample in the complete Radial Velocity Experiment (RAVE) Data Release 5 catalog with the goal of identifying and characterizing all very metal-poor stars observed by the survey. Using a three-stage method, we first identified the candidate stars using only their spectra as input information. We employed an algorithm called t-SNE to construct a low-dimensional projection of the spectrum space and isolate the region containing metal-poor stars. Following this step, we measured the equivalent widths of the near-infrared CaII triplet lines with a method based on flexible Gaussian processes to model the correlated noise present in the spectra. In the last step, we constructed a calibration relation that converts the measured equivalent widths and the color information coming from the 2MASS and WISE surveys into metallicity and temperature estimates. We identified 877 stars with at least a 50% probability of being very metal-poor $(\rm [Fe/H] < -2\,\rm dex)$, out of which 43 are likely extremely metal-poor $(\rm [Fe/H] < -3\,\rm dex )$. The comparison of the derived values to a small subsample of stars with literature metallicity values shows that our method works reliably and correctly estimates the uncertainties, which typically have values $\sigma_{\rm [Fe/H]} \approx 0.2\,\mathrm{dex}$. In addition, when compared to the metallicity results derived using the RAVE DR5 pipeline, it is evident that we achieve better accuracy than the pipeline and therefore more reliably evaluate the very metal-poor subsample. Based on the repeated observations of the same stars, our method gives very consistent results. The method used in this work can also easily be extended to other large-scale data sets, including to the data from the Gaia mission and the upcoming 4MOST survey.

We characterize the selection function of RAVE using 2MASS as our underlying population, which we assume represents all stars which could have potentially been observed. We evaluate the completeness fraction as a function of position, magnitude, and color in two ways: first, on a field-by-field basis, and second, in equal-size areas on the sky. Then, we consider the effect of the RAVE stellar parameter pipeline on the final resulting catalogue, which in principle limits the parameter space over which our selection function is valid. Our final selection function is the product of the completeness fraction and the selection function of the pipeline. We then test if the application of the selection function introduces biases in the derived parameters. To do this, we compare a parent mock catalogue generated using Galaxia with a mock-RAVE catalogue where the selection function of RAVE has been applied. We conclude that for stars brighter than I = 12, between $4000 \rm K < T_{\rm eff} < 8000 \rm K$ and $0.5 < \rm{log}\,g < 5.0$, RAVE is kinematically and chemically unbiased with respect to expectations from Galaxia.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

A large sample of over 38,000 chromospherically active candidate solar-like stars and cooler dwarfs from the RAVE survey is addressed in this paper. An improved activity identification with respect to the previous study was introduced to build a catalog of field stars in the Solar neighborhood with an excess emission flux in the calcium infrared triplet wavelength region. The central result of this work is the calibration of the age--activity relation for the main sequence dwarfs in a range from a few $10 \; \mathrm{Myr}$ up to a few Gyr. It enabled an order of magnitude age estimation of the entire active sample. Almost 15,000 stars are shown to be younger than $1\;\mathrm{Gyr}$ and $\sim$2000 younger than $100\;\mathrm{Myr}$. The young age of the most active stars is confirmed by their position off the main sequence in the $J-K$ versus $N_{UV}-V$ diagram showing strong ultraviolet excess, mid-infrared excess in the $J-K$ versus $W_1-W_2$ diagram and very cool temperatures ($J-K>0.7$). They overlap with the reference pre-main sequence RAVE stars often displaying X-ray emission. The activity level increasing with the color reveals their different nature from the solar-like stars and probably represents an underlying dynamo generating magnetic fields in cool stars. 50\% of the RAVE objects from DR5 are found in the TGAS catalog and supplemented with accurate parallaxes and proper motions by Gaia. This makes the database of a large number of young stars in a combination with RAVE's radial velocities directly useful as a tracer of the very recent large-scale star formation history in the Solar neighborhood. The data are available online in the Vizier database.

Both the three-dimensional density of red clump giants and the gas kinematics in the inner Galaxy indicate that the pattern speed of the Galactic bar could be much lower than previously estimated. Here, we show that such slow bar models are unable to reproduce the bimodality observed in local stellar velocity space. We do so by computing the response of stars in the Solar neighbourhood to the gravitational potential of slow and fast bars, in terms of their perturbed distribution function in action-angle space up to second order, as well as by identifying resonantly trapped orbits. We also check that the bimodality is unlikely to be produced through perturbations from spiral arms, and conclude that, contrary to gas kinematics, local stellar kinematics still favour a fast bar in the Milky Way, with a pattern speed of the order of almost twice (and no less than 1.8 times) the circular frequency at the Sun's position. This leaves open the question of the nature of the long flat extension of the bar in the Milky Way.

The orbits, atmospheric parameters, chemical abundances, and ages of individual stars in the Milky Way provide the most comprehensive illustration of galaxy formation available. The Tycho-Gaia Astrometric Solution (TGAS) will deliver astrometric parameters for the largest ever sample of Milky Way stars, though its full potential cannot be realized without the addition of complementary spectroscopy. Among existing spectroscopic surveys, the RAdial Velocity Experiment (RAVE) has the largest overlap with TGAS ($\gtrsim$200,000 stars). We present a data-driven re-analysis of 520,781 RAVE spectra using The Cannon. For red giants, we build our model using high-fidelity APOGEE stellar parameters and abundances for stars that overlap with RAVE. For main-sequence and sub-giant stars, our model uses stellar parameters from the K2/EPIC. We derive and validate effective temperature $T_{\rm eff}$, surface gravity $\log{g}$, and chemical abundances of up to seven elements (O, Mg, Al, Si, Ca, Fe, Ni). We report a total of 1,685,851 elemental abundances with a typical precision of 0.07 dex, a substantial improvement over previous RAVE data releases. The synthesis of RAVE-on and TGAS is the most powerful data set for chemo-dynamic analyses of the Milky Way ever produced.

We present evidence from the RAdial Velocity Experiment (RAVE) survey of chemically separated, kinematically distinct disc components in the solar neighbourhood. We apply probabilistic chemical selection criteria to separate our sample into $\alpha$-low (`thin disc') and $\alpha$-high (`thick disc') sequences. Using newly derived distances, which will be utilized in the upcoming RAVE DR5, we explore the kinematic trends as a function of metallicity for each of the disc components. For our $\alpha$-low disc, we find a negative trend in the mean rotational velocity ($V_{\mathrm{\phi}}$) as a function of iron abundance ([Fe/H]). We measure a positive gradient $\partial V_{\mathrm{\phi}}$/$\partial$[Fe/H] for the $\alpha$-high disc, consistent with results from high-resolution surveys. We also find differences between the $\alpha$-low and $\alpha$-high discs in all three components of velocity dispersion. We discuss the implications of an $\alpha$-low, metal-rich population originating from the inner Galaxy, where the orbits of these stars have been significantly altered by radial mixing mechanisms in order to bring them into the solar neighbourhood. The probabilistic separation we propose can be extended to other data sets for which the accuracy in [$\alpha$/Fe] is not sufficient to disentangle the chemical disc components a priori. For such datasets which will also have significant overlap with Gaia DR1, we can therefore make full use of the improved parallax and proper motion data as it becomes available to investigate kinematic trends in these chemical disc components.

Using six-dimesional phase-space information from the Fourth Data release of the Radial Velocity Experiment (RAVE) over the range of Galactic longitude 240$^{\circ}< l <$ 360$^{\circ}$ and $V_{LSR} < -239$ kms$^{-1}$, we have computed orbits for 329 RAVE stars that were originally selected as chemically and kinematically related to $\omega$ Centauri. The orbits were integrated in a Milky-Way-like axisymmetric Galactic potential, ignoring the effects of the dynamical evolution of $\omega$ Centauri due to the tidal effects of the Galaxy disk on the cluster along time. We also ignored secular changes in the Milky Way potential over time. In a Monte Carlo scheme, and under the assumption that the stars may have been ejected with velocities greater than the escape velocity ($V_{rel}>V_{esc,0}$) from the cluster, we identified 15 stars as having close encounters with $\omega$ Centauri: (\textit{i}) 8 stars with relative velocities $V_{rel}< 200 $ kms$^{-1}$ may have been ejected $\sim$ 200 Myr ago from $\omega$ Centauri; (\textit{ii}) other group of 7 stars were identified with high relative velocity $V_{rel}> 200 $ kms$^{-1}$ during close encounters, and seems unlikely that they have been ejected from $\omega$ Centauri. We also confirm the link between J131340.4-484714 as potential member of $\omega$ Centauri, and probably ejected $\sim$ 2.0 Myr ago, with a relative velocity $V_{rel}\sim80$ kms$^{-1}$.

We present an estimate of the third integral of motion for axisymmetric three-dimensional potentials. This estimate is based on a Staeckel approximation and is explicitly written as a function of the potential. We tested this scheme for the Besancon Galactic model and two other disc-halo models and find that orbits of disc stars have an accurately conserved third quasi integral. The accuracy ranges from of 0.1% to 1% for heights varying from z = 0~kpc to z= 6 kpc and Galactocentric radii R from 5 to 15kpc. We also tested the usefulness of this quasi integral in analytic distribution functions of disc stellar populations: we show that the distribution function remains approximately stationary and that it allows to recover the potential and forces by applying Jeans equations to its moments.

We present the identification of potential members of nearby Galactic globular clusters using radial velocities from the RAdial Velocity Experiment Data Release 4 (RAVE-DR4) survey database. Our identifications are based on three globular clusters -- NGC 3201, NGC 5139 ($\omega$ Cen) and NGC 362 -- all of which are shown to have |RV|>100 km/s. The identification of globular cluster stars in RAVE DR4 data offers a unique opportunity to test the precision and accuracy of the stellar parameters determined with the currently available Stellar Parameter Pipelines (SPPs) used in the survey, as globular clusters are ideal testbeds for the validation of stellar atmospheric parameters, abundances, distances and ages. For both NGC 3201 and $\omega$ Cen, there is compelling evidence for numerous members (> 10) in the RAVE database; in the case of NGC 362 the evidence is more ambiguous, and there may be significant foreground and/or background contamination in our kinematically-selected sample. A comparison of the RAVE-derived stellar parameters and abundances with published values for each cluster and with BASTI isochrones for ages and metallicities from the literature reveals overall good agreement, with the exception of the apparent underestimation of surface gravities for giants, in particular for the most metal-poor stars. Moreover, if the selected members are part of the main body of each cluster our results would also suggest that the distances from Binney et al. 2013, where only isochrones more metal-rich than -0.9 dex were used, are typically underestimated by ~ 40% with respect to the published distances for the clusters, while the distances from Zwitter et al. 2010 show stars ranging from 1 to ~ 6.5 kpc -- with indications of a trend toward higher distances at lower metallicities -- for the three clusters analysed in this study.

This is a brief rebuttal to arXiv:1502.03821, which claims to provide the first observational proof of dark matter interior to the solar circle. We point out that this result is not new, and can be traced back at least a quarter century.

We study the kinematics of a local sample of stars, located within a cylinder of 500 pc radius centered on the Sun, in the RAVE dataset. We find clear asymmetries in the $v_R$-$v_\phi$ velocity distributions of thin and thick disk stars: here are more stars moving radially outwards for low azimuthal velocities and more radially inwards for high azimuthal velocities. Such asymmetries have been previously reported for the thin disk as being due to the Galactic bar, but this is the first time that the same type of structures are seen in the thick disk. Our findings imply that the velocities of thick disk stars should no longer be described by Schwarzschild's, multivariate Gaussian or purely axisymmetric distributions. Furthermore, the nature of previously reported substructures in the thick disk needs to be revisited as these could be associated with dynamical resonances rather than to accretion events. It is clear that dynamical models of the Galaxy must fit the 3D velocity distributions of the disks, rather than the projected 1D, if we are to understand the Galaxy fully.

We determine the Galactic potential in the solar neigbourhood from RAVE observations. We select red clump stars for which accurate distances, radial velocities, and metallicities have been measured. Combined with data from the 2MASS and UCAC catalogues, we build a sample of 4600 red clump stars within a cylinder of 500 pc radius oriented in the direction of the South Galactic Pole, in the range of 200 pc to 2000 pc distances. We deduce the vertical force and the total mass density distribution up to 2 kpc away from the Galactic plane by fitting a distribution function depending explicitly on three isolating integrals of the motion in a separable potential locally representing the Galactic one with four free parameters. Because of the deep extension of our sample, we can determine nearly independently the dark matter mass density and the baryonic disc surface mass density. We find (i) at 1kpc Kz/(2piG) = 68.5 pm 1.0 Msun/pc2, and (ii) at 2 kpc Kz/(2piG) = 96.9 pm 2.2 Msun/pc2. Assuming the solar Galactic radius at R0 = 8.5 kpc, we deduce the local dark matter density rhoDM (z=0) = 0.0143 pm 0.0011Msun pc3 = 0.542 pm 0.042 Gev/cm3 and the baryonic surface mass density Sigma = 44.4 pm 4.1 Msun/pc2 . Our results are in agreement with previously published Kz determinations up to 1 kpc, while the extension to 2 kpc shows some evidence for an unexpectedly large amount of dark matter. A flattening of the dark halo of order 0.8 can produce such a high local density in combination with a circular velocity of 240 km/s . Another explanation, allowing for a lower circular velocity, could be the presence of a secondary dark component, a very thick disc resulting either from the deposit of dark matter from the accretion of multiple small dwarf galaxies, or from the presence of an effective phantom thick disc in the context of effective galactic-scale modifications of gravity.

We use the kinematics of $\sim200\,000$ giant stars that lie within $\sim 1.5$ kpc of the plane to measure the vertical profile of mass density near the Sun. We find that the dark mass contained within the isodensity surface of the dark halo that passes through the Sun ($(6\pm0.9)\times10^{10}\,\mathrm{M_\odot}$), and the surface density within $0.9$ kpc of the plane ($(69\pm10)\,\mathrm{M_\odot\,pc^{-2}}$) are almost independent of the (oblate) halo's axis ratio $q$. If the halo is spherical, 46 per cent of the radial force on the Sun is provided by baryons, and only 4.3 per cent of the Galaxy's mass is baryonic. If the halo is flattened, the baryons contribute even less strongly to the local radial force and to the Galaxy's mass. The dark-matter density at the location of the Sun is $0.0126\,q^{-0.89}\,\mathrm{M_\odot\,pc^{-3}}=0.48\,q^{-0.89}\,\mathrm{GeV\,cm^{-3}}$. When combined with other literature results we find hints for a mildly oblate dark halo with $q \simeq 0.8$. Our value for the dark mass within the solar radius is larger than that predicted by cosmological dark-matter-only simulations but in good agreement with simulations once the effects of baryonic infall are taken into account. Our mass models consist of three double-exponential discs, an oblate bulge and a Navarro-Frenk-White dark-matter halo, and we model the dynamics of the RAVE stars in the corresponding gravitational fields by finding distribution functions $f(\mathbf{J})$ that depend on three action integrals. Statistical errors are completely swamped by systematic uncertainties, the most important of which are the distance to the stars in the photometric and spectroscopic samples and the solar distance to the Galactic centre. Systematics other than the flattening of the dark halo yield overall uncertainties $\sim 15$ per cent.

The diffuse interstellar bands (DIBs) are absorption lines observed in visual and near infrared spectra of stars. Understanding their origin in the interstellar medium (ISM) is one of the oldest problems in astronomical spectroscopy, as DIBs ahave been known since 1922. In a completely new approach to understanding DIBs, we combined information from nearly 500,000 stellar spectra obtained by the massive spectroscopic survey RAVE (Radial Velocity Experiment) to produce the first pseudo three-dimensional map of the strength of the DIB at 8620~{\AA}ngstroms covering the nearest 3 kiloparsecs from the Sun, and show that it follows our independently constructed spatial distribution of extinction by interstellar dust along the Galactic plane. Despite having a similar distribution in the Galactic plane, the DIB 8620 carrier has a significantly larger vertical scale height than the dust. Even if one DIB may not represent the general DIB population, our observations outline the future direction of DIB research.

We provide new constraints on the chemo-dynamical models of the Milky Way by measuring the radial and vertical chemical gradients for the elements Mg, Al, Si, Ti, and Fe in the Galactic disc and the gradient variations as a function of the distance from the Galactic plane ($Z$). We selected a sample of giant stars from the RAVE database using the gravity criterium 1.7$<$log g$<$2.8. We created a RAVE mock sample with the Galaxia code based on the Besan\c con model and selected a corresponding mock sample to compare the model with the observed data. We measured the radial gradients and the vertical gradients as a function of the distance from the Galactic plane $Z$ to study their variation across the Galactic disc. The RAVE sample exhibits a negative radial gradient of $d[Fe/H]/dR=-0.054$ dex kpc$^{-1}$ close to the Galactic plane ($|Z|<0.4$ kpc) that becomes flatter for larger $|Z|$. Other elements follow the same trend although with some variations from element to element. The mock sample has radial gradients in fair agreement with the observed data. The variation of the gradients with $Z$ shows that the Fe radial gradient of the RAVE sample has little change in the range $|Z|\lesssim0.6$ kpc and then flattens. The iron vertical gradient of the RAVE sample is slightly negative close to the Galactic plane and steepens with $|Z|$. The mock sample exhibits an iron vertical gradient that is always steeper than the RAVE sample. The mock sample also shows an excess of metal-poor stars in the [Fe/H] distributions with respect to the observed data. These discrepancies can be reduced by decreasing the number of thick disc stars and increasing their average metallicity in the Besan\c con model.

Aim: 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 Besancon model. Methods: 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 Besancon model. We measured the chemical gradients as functions of the guiding radius (Rg) at different distances from the Galactic plane reached by the stars along their orbit (Zmax). Results: The chemical gradients of the RAVE and GCS samples are negative and show consistent trends, although they are not equal: at Zmax<0.4 kpc and 4.5<Rg(kpc)<9.5, the iron gradient for the RAVE sample is d[Fe/H]/dRg=-0.065 dex kpc^{-1}, whereas for the GCS sample it is d[Fe/H]/dRg=-0.043 dex kpc^{-1} with internal errors +-0.002 and +-0.004 dex kpc^{-1}, respectively. The gradients of the RAVE and GCS samples become flatter at larger Zmax. Conversely, the mock sample has a positive iron gradient of d[Fe/H]/dRg=+0.053+-0.003 dex kpc^{-1} at Zmax<0.4 kpc and remains positive at any Zmax. These positive and unrealistic values originate from the lack of correlation between metallicity and tangential velocity in the Besancon 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 Besancon model.

RAVE, the unbiased magnitude limited survey of the southern sky stars, contained 456,676 medium-resolution spectra at the time of our analysis. Spectra cover the CaII IRT range which is a known indicator of chromospheric activity. Our previous work (Matijevi\v{c} et al. 2012) classified all spectra using locally linear embedding. It identified 53,347 cases with a suggested emission component in calcium lines. Here we use a spectral subtraction technique to measure the properties of this emission. Synthetic templates are replaced by the observed spectra of non-active stars to bypass the difficult computations of non-LTE profiles of the line cores and stellar parameter dependence. We derive both the equivalent width of the excess emission for each calcium line on a 5\AA\ wide interval and their sum EW_IRT for ~44,000 candidate active dwarf stars with S/N>20 and with no respect to the source of their emission flux. From these ~14,000 show a detectable chromospheric flux with at least 2\sigma\ confidence level. Our set of active stars vastly enlarges previously known samples. Atmospheric parameters and in some cases radial velocities of active stars derived from automatic pipeline suffer from systematic shifts due to their shallower calcium lines. We re-estimate the effective temperature, metallicity and radial velocities for candidate active stars. The overall distribution of activity levels shows a bimodal shape, with the first peak coinciding with non-active stars and the second with the pre main-sequence cases. The catalogue will be publicly available with the next RAVE public data releases.

Context. Galactic open clusters (OCs) mainly belong to the young stellar population in the Milky Way disk, but are there groups and complexes of OCs that possibly define an additional level in hierarchical star formation? Current compilations are too incomplete to address this question, especially regarding radial velocities (RVs) and metallicities ($[M/H]$). Aims. Here we provide and discuss newly obtained RV and $[M/H]$ data, which will enable us to reinvestigate potential groupings of open clusters and associations. Methods.We extracted additional RVs and $[M/H]$ from the RAdial Velocity Experiment (RAVE) via a cross-match with the Catalogue of Stars in Open Cluster Areas (CSOCA). For the identified OCs in RAVE we derived RV and $[M/H]$ from a cleaned working sample and compared the results with previous findings. Results. Although our RAVE sample does not show the same accuracy as the entire survey, we were able to derive reliable RV for 110 Galactic open clusters. For 37 OCs we publish RV for the first time. Moreover, we determined $[M/H]$ for 81 open clusters, extending the number of OCs with $[M/H]$ by 69.

The determination of the LSR is still a matter of debate. The classical value of the tangential peculiar motion of the Sun with respect to the LSR was challenged in recent years, claiming a significantly larger value. We show that the RAdial Velocity Experiment (RAVE) sample of dwarf stars is an excellent data set to derive tighter boundary conditions to chemodynamical evolution models of the extended solar neighbourhood. We present an improved Jeans analysis, which allows a better interpretation of the measured kinematics of stellar populations in the Milky Way disc. We propose an improved version of the Str\"omberg relation with the radial scalelengths as the only unknown. Binning RAVE stars in metallicity reveals a bigger asymmetric drift (corresponding to a smaller radial scalelength) for more metal-rich populations. With the standard assumption of velocity-dispersion independent radial scalelengths in each metallicity bin, we redetermine the LSR. The new Str\"omberg equation yields a joint LSR value of V_\sun=3.06 \pm 0.68 km/s, which is even smaller than the classical value based on Hipparcos data. The corresponding radial scalelength increases from 1.6 kpc for the metal-rich bin to 2.9 kpc for the metal-poor bin, with a trend of an even larger scalelength for young metal-poor stars. When adopting the recent Sch\"onrich value of V_\sun=12.24 km/s for the LSR, the new Str\"omberg equation yields much larger individual radial scalelengths of the RAVE subpopulations, which seem unphysical in part. The new Str\"omberg equation allows a cleaner interpretation of the kinematic data of disc stars in terms of radial scalelengths. Lifting the LSR value by a few km/s compared to the classical value results in strongly increased radial scalelengths with a trend of smaller values for larger velocity dispersions.

We determine approximate numerical integrals of motion of 2D symmetric Hamiltonian systems. We detail for a few gravitational potentials the conditions under which quasi-integrals are obtained as polynomial series. We show that each of these potentials has a wide range of regular orbits that are accurately modelled with a unique approximate integral of motion.

We performed a spectroscopic survey of nearly 700 stars probing the galactic thick disc far from the solar neighbourhood towards the galactic coordinates (l~277, b~47). The derived effective temperatures, surface gravities and overall metallicities were then combined with stellar evolution isochrones, radial velocities and proper motions to derive the distances, kinematics and orbital parameters of the sample stars. The targets belonging to each galactic component (thin disc, thick disc, halo) were selected either on their kinematics or according to their position above the galactic plane, and the vertical gradients were also estimated. We present here atmospheric parameters, distances and kinematics for this sample, and a comparison of our kinematic and metallicity distributions with the Besancon model of the Milky Way. The thick disc far from the solar neighbourhood is found to differ only slightly from the thick disc properties as derived in the solar vicinity. For regions where the thick disc dominates, we measured vertical velocity and metallicity trends of d(V_phi)/dZ = 19 +/- 8 km/s/kpc and d[M/H]/dZ = -0.14 +/- 0.05 dex/kpc, respectively. These trends can be explained as a smooth transition between the different galactic components, although intrinsic gradients could not be excluded. In addition, a correlation d(V_phi)/d[M/H] = -45 +/- 12 km/s/dex between the orbital velocity and the metallicity of the thick disc is detected. This gradient is inconsistent with the SDSS photometric survey analysis, which did not detect any such trend, and challenges radial migration models of thick disc formation. Estimations of the scale heights and scale lengths for different metallicity bins of the thick disc result in consistent values, with hR~3.4 \pm 0.7 kpc, and hZ~694 \pm 45 pc, showing no evidence of relics of destroyed massive satellites.

We present chemical elemental abundances for $36,561$ stars observed by the RAdial Velocity Experiment (RAVE), an ambitious spectroscopic survey of our Galaxy at Galactic latitudes $|$b$|>25^{\circ}$ and with magnitudes in the range 9$<I_{DENIS}<$13. RAVE spectra cover the Ca-triplet region at 8410--8795\AA\ with resolving power R$\sim$7500. This first data release of the RAVE chemical catalogue is complementary to the third RAVE data release of radial velocities and stellar parameters, and it contains chemical abundances for the elements Mg, Al, Si, Ca, Ti, Fe and Ni, with a mean error of $\sim$0.2 dex, as judged from accuracy tests performed on synthetic and real spectra. Abundances are estimated through a dedicated processing pipeline in which the curve of growth of individual lines is obtained from a library of absorption-line equivalent widths to construct a model spectrum that is then matched to the observed spectrum via a $\chi^2$-minimization technique. We plan to extend this pipeline to include estimates for other elements, such as oxygen and sulfur, in future data releases.

The Hyades stream has long been thought to be a dispersed vestige of the Hyades cluster. However, recent analyses of the parallax distribution, of the mass function, and of the action-space distribution of stream stars have shown it to be rather composed of orbits trapped at a resonance of a density disturbance. This resonant scenario should leave a clearly different signature in the element abundances of stream stars than the dispersed cluster scenario, since the Hyades cluster is chemically homogeneous. Here, we study the metalllicty as well as the element abundances of Li, Na, Mg, Fe, Zr, Ba, La, Ce, Nd, and Eu for a random sample of stars belonging to the Hyades stream, and compare them with those of stars from the Hyades cluster. From this analysis: (i) we independently confirm that the Hyades stream cannot be solely composed of stars originating in the Hyades cluster; (ii) we show that some stars from the Hyades stream nevertheless have abundances compatible with an origin in the cluster; (iii) we emphasize that the use of Li as a chemical tag of the cluster origin of main-sequence stars is very efficient in the range 5500 - 6200 K, since the Li sequence in the Hyades cluster is very tight; (iv) we show that, while this evaporated population has a metallicity excess of ~0.2 dex w.r.t. the local thin disk population, identical to that of the Hyades cluster, the remainder of the Hyades stream population has still a metallicity excess of ~0.06 to 0.15 dex, consistent with an origin in the inner Galaxy; (v) we show that the Hyades stream can be interpreted as an inner 4:1 resonance of the spiral pattern: this then also reproduces an orbital family compatible with the Sirius stream, and places the origin of the Hyades stream up to 1 kpc inwards from the solar radius, which might explain the observed metallicity excess of the stream population.

We identify a new, nearby (0.5 < d < 10 kpc) stream in data from the RAdial Velocity Experiment (RAVE). As the majority of stars in the stream lie in the constellation of Aquarius we name it the Aquarius Stream. We identify 15 members of the stream lying between 30 < l < 75 and -70< b <-50, with heliocentric line-of-sight velocities V_los~-200 km/s. The members are outliers in the radial velocity distribution, and the overdensity is statistically significant when compared to mock samples created with both the Besan\c{c}on Galaxy model and newly-developed code Galaxia. The metallicity distribution function and isochrone fit in the log g - T_eff plane suggest the stream consists of a 10 Gyr old population with [m/H]~-1.0. We explore relations to other streams and substructures, finding the stream cannot be identified with known structures: it is a new, nearby substructure in the Galaxy's halo. Using a simple dynamical model of a dissolving satellite galaxy we account for the localization of the stream. We find that the stream is dynamically young and therefore likely the debris of a recently disrupted dwarf galaxy or globular cluster. The Aquarius stream is thus a specimen of ongoing hierarchical Galaxy formation, rare for being right in the solar suburb.

The RAdial Velocity Experiment (RAVE) is a spectroscopic survey of the Milky Way. We use the subsample of spectra with spectroscopically determined values of stellar parameters to determine the distances to these stars. The list currently contains 235,064 high quality spectra which show no peculiarities and belong to 210,872 different stars. The numbers will grow as the RAVE survey progresses. The public version of the catalog will be made available through the CDS services along with the ongoing RAVE public data releases. The distances are determined with a method based on the work by Breddels et al.~(2010). Here we assume that the star undergoes a standard stellar evolution and that its spectrum shows no peculiarities. The refinements include: the use of either of the three isochrone sets, a better account of the stellar ages and masses, use of more realistic errors of stellar parameter values, and application to a larger dataset. The derived distances of both dwarfs and giants match within ~21% to the astrometric distances of Hipparcos stars and to the distances of observed members of open and globular clusters. Multiple observations of a fraction of RAVE stars show that repeatability of the derived distances is even better, with half of the objects showing a distance scatter of \simlt 11%. RAVE dwarfs are ~300 pc from the Sun, and giants are at distances of 1 to 2 kpc, and up to 10 kpc. This places the RAVE dataset between the more local Geneva-Copenhagen survey and the more distant and fainter SDSS sample. As such it is ideal to address some of the fundamental questions of Galactic structure and evolution in the pre-Gaia era. Individual applications are left to separate papers, here we show that the full 6-dimensional information on position and velocity is accurate enough to discuss the vertical structure and kinematic properties of the thin and thick disks.

We have chosen the name of GYES, one of the mythological giants with one hundred arms, offspring of Gaia and Uranus, for our instrument study of a multifibre spectrograph for the prime focus of the Canada-France-Hawaii Telescope. Such an instrument could provide an excellent ground-based complement for the Gaia mission and a northern complement to the HERMES project on the AAT. The CFHT is well known for providing a stable prime focus environment, with a large field of view, which has hosted several imaging instruments, but has never hosted a multifibre spectrograph. Building upon the experience gained at GEPI with FLAMES-Giraffe and X-Shooter, we are investigating the feasibility of a high multiplex spectrograph (about 500 fibres) over a field of view 1 degree in diameter. We are investigating an instrument with resolution in the range 15000 to 30000, which should provide accurate chemical abundances for stars down to 16th magnitude and radial velocities, accurate to 1 km/s for fainter stars. The study is led by GEPI-Observatoire de Paris with a contribution from Oxford for the study of the positioner. The financing for the study comes from INSU CSAA and Observatoire de Paris. The conceptual study will be delivered to CFHT for review by October 1st 2010.

Theories of thick disk formation can be differentiated by measurements of stellar elemental abundances. We have undertaken a study of metal-poor stars selected from the RAVE spectroscopic survey of bright stars to establish whether or not there is a significant population of metal-poor thick-disk stars ([Fe/H] <~ -1.0) and to measure their elemental abundances. In this paper, we present abundances of four alpha-elements (Mg, Si, Ca, Ti) and iron for a subsample of 212 RGB and 31 RC/HB stars from this study. We find that the [alpha/Fe] ratios are enhanced implying that enrichment proceeded by purely core-collapse supernovae. This requires that star formation in each star forming region had a short duration. The relative lack of scatter in the [alpha/Fe] ratios implies good mixing in the ISM prior to star formation. In addition, the ratios resemble that of the halo, indicating that the halo and thick disk share a similar massive star IMF. We conclude that the alpha enhancement of the metal-poor thick disk implies that direct accretion of stars from dwarf galaxies similar to surviving dwarf galaxies today did not play a major role in the formation of the thick disk.

We derive a simple analytical expression for the two-body force in a sub-class of MOND-like theories and make testable predictions in the modification to the two-body orbital period, shape, and precession rate, and escape speed etc. We demonstrate the applications of the modified Kepler's law in the timing of satellite orbits around the Milky Way, and checking the feasibility of MOND in the orbit of Large Magellanic Cloud, the M31 galaxy, and the merging Bullet Clusters. MOND appears to be consistent with satellite orbits although with a tight margin. Our results on two-bodies are also generalized to restricted three-body, many-body problems, rings and shells.

CONTEXT: The RAVE spectroscopic survey for galactic structure and evolution obtains 8400-8800 Ang spectra at 7500 resolving power at the UK Schmidt Telescope using the 6dF multi-fiber positioner. More than 300,000 9<I<12 and |b|>25 deg southern stars have been observed to date. AIMS: This paper presents the first intrinsic examination of stellar spectra from the RAVE survey, aimed at evaluating their diagnostic potential for peculiar stars and at contributing to the general understanding of Luminous Blue Variables (LBVs). METHODS: We used the multi-epoch spectra for all seven LBVs observed, between 2005 and 2008, in the Large Magellanic Cloud (LMC) by the RAVE survey. RESULTS: We demonstrate that RAVE spectra possess significant diagnostic potential when applied to peculiar stars and, in particular, LBVs. The behaviour of the radial velocities for both emission and absorption lines, and the spectral changes between outburst and quiescence states are described and found to agree with evidence gathered at more conventional wavelengths. The wind outflow signatures and their variability are investigated, with multi-components detected in S Doradus. Photoionisation modelling of the rich emission line spectrum of R 127 shows evidence of a massive detached ionised shell that was ejected during the 1982-2000 outburst. Surface inhomogeneities in the nuclear-processed material, brought to the surface by heavy mass loss, could have been observed in S Doradus, even if alternative explanations are possible. We also detect the transition from quiescence to outburst state in R 71. Finally, our spectrum of R 84 offers one of the clearest views of its cool companion.

We look for observational signatures that could discriminate between Newtonian and modified Newtonian (MOND) dynamics in the Milky Way, in view of the advent of large astrometric and spectroscopic surveys. Indeed, a typical signature of MOND is an apparent disk of "phantom" dark matter, which is uniquely correlated with the visible disk-density distribution. Due to this phantom dark disk, Newtonian models with a spherical halo have different signatures from MOND models close to the Galactic plane. The models can thus be differentiated by measuring dynamically (within Newtonian dynamics) the disk surface density at the solar radius, the radial mass gradient within the disk, or the velocity ellipsoid tilt angle above the Galactic plane. Using the most realistic possible baryonic mass model for the Milky Way, we predict that, if MOND applies, the local surface density measured by a Newtonist will be approximately 78 Msun/pc2 within 1.1 kpc of the Galactic plane, the dynamically measured disk scale-length will be enhanced by a factor of 1.25 with respect to the visible disk scale-length, and the local vertical tilt of the velocity ellipsoid at 1 kpc above the plane will be approximately 6 degrees. None of these tests can be conclusive for the present-day accuracy of Milky Way data, but they will be of prime interest with the advent of large surveys such as GAIA.

We define new potential-density pairs and examine the impact of the potential flattening on the vertical velocity ellipsoid tilt, $\delta$. By means of numerical integrations and analytical calculations, we estimate $\delta$ in a variety of galactic potentials. We show that at 1 kpc above the Galactic plane at the solar radius, $\delta$ can differ by 5 degrees, depending on whether the dark matter halo is flat or spherical. This result excludes the possibility of an extremely flattened Galactic dark halo.

We present a measure of the inclination of the velocity ellipsoid at 1 kpc below the Galactic plane using a sample of red clump giants from the RAVE DR2 release. We find that the velocity ellipsoid is tilted towards the Galactic plane with an inclination of 7.3 +/-1.8 degree. We compare this value to computed inclinations for two mass models of the Milky Way. We find that our measurement is consistent with a short scale length of the stellar disc (Rd ~2 kpc) if the dark halo is oblate or with a long scale length (Rd~3 kpc) if the dark halo is prolate. Once combined with independent constraints on the flattening of the halo, our measurement suggests that the scale length is approximately halfway between these two extreme values, with a preferred range [2.5-2.7] kpc for a nearly spherical halo. Nevertheless, no model can be clearly ruled out. With the continuation of the RAVE survey, it will be possible to provide a strong constraint on the mass distribution of the Milky Way using refined measurements of the orientation of the velocity ellipsoid.

We analyze the distribution of G and K type stars towards the Galactic poles using RAVE and ELODIE radial velocities, 2MASS photometric star counts, and UCAC2 proper motions. The combination of photometric and 3D kinematic data allows us to disentangle and describe the vertical distribution of dwarfs, sub-giants and giants and their kinematics. We identify discontinuities within the kinematics and magnitude counts that separate the thin disk, thick disk and a hotter component. The respective scale heights of the thin disk and thick disk are 225$\pm$10 pc and 1048$\pm$36 pc. We also constrain the luminosity function and the kinematic distribution function. The existence of a kinematic gap between the thin and thick disks is incompatible with the thick disk having formed from the thin disk by a continuous process, such as scattering of stars by spiral arms or molecular clouds. Other mechanisms of formation of the thick disk such as `created on the spot' or smoothly `accreted' remain compatible with our findings.

We used red clump stars to measure the surface mass density of the Galactic disk in the solar neighbourhood. High resolution spectra of red clump stars towards the NGP have been obtained with the ELODIE spectrograph at OHP for Tycho-2 selected stars, and nearby Hipparcos counterparts were also observed. We determined their distances, velocities, and metallicities to measure the gravitational force law perpendicular to the Galactic plane. As in most previous studies, we applied one-parameter models of the vertical gravitational potential. We obtained a disk surface mass density within 1.1kpc of the Galactic plane, Sigma_{1.1kpc}=64+/-5 Msun_pc^{-2}, with an excellent formal accuracy, however we found that such one-parameter models can underestimate the real uncertainties. Applying two-parameter models, we derived more realistic estimates of the total surface mass density within 800pc from the Galactic plane, Sigma_{0.8kpc}=57-66 Msun pc^{-2}, and within 1.1kpc, Sigma{1.1kpc}=57-79 Msun pc^{-2}. This can be compared to literature estimates of \sim40 Msun pc^{-2} in stars and to 13 Msun pc^{-2} in the less accurately measured ISM contribution. We conclude that there is no evidence of large amounts of dark matter in the disk and, furthermore, that the dark matter halo is round or not vey much flattened. A by-product of this study is the determination of the half period of oscillation by the Sun through the Galactic plane, 42+/-2Myr, which cannot be related to the possible period of large terrestrial impact craters \sim 33-37Myr.

Low-mass and brown dwarfs have recently been found as wide companions to many nearby stars, formerly believed to be single. Wide binaries are usually found as common proper motion pairs. Sometimes, more than two objects share the same large proper motion, identifying them as nearby systems. We have found a third, low-mass component to a known wide binary at a distance of about 21pc, consisting of a red and a white dwarf (LHS4039 and LHS4040; about 150AU separation). The new companion, APMPMJ2354-3316C separated by about 2200AU, was classified as M8.5 dwarf. In recent spectroscopic observations it shows a very strong $H_{\alpha}$ emission line and blue continuum. Comparing this event to flares in late-type M dwarfs, we find some similarity with LHS2397a, a nearby M8 dwarf which is so far the only known example of a low-mass star with a tight brown dwarf companion (separation of less than 4AU). The level of the activity as measured by $L_{H_{\alpha}}/L_{bol}$ is comparable to that of the M9.5 dwarf 2MASSWJ0149+29 both during the flare and in quiescence.

Nearly 400 Tycho-2 stars have been observed in a 720 square degree field in the direction of the North Galactic Pole with the high resolution echelle spectrograph ELODIE. Absolute magnitudes, effective temperatures, gravities and metallicities have been estimated, as well as distances and 3D velocities. Most of these stars are clump giants and span typical distances from 200pc to 800pc to the galactic mid-plane. This new sample, free of any kinematical and metallicity bias, is used to investigate the vertical distribution of disk stars. The old thin disk and thick disk populations are deconvolved from the velocity-metallicity distribution of the sample and their parameters are determined. The thick disk is found to have a moderate rotational lag of -51+-5 km/s with respect to the Sun with velocity ellipsoid (sigma_U, sigma_V, sigma_W)=(63+-6, 39+-4, 39+-4) km/s, mean metallicity of [Fe/H]=-0.48+-0.05 and a high local normalization of 15+-7%. Combining this NGP sample with a local sample of giant stars from the Hipparcos catalogue, the orientation of the velocity ellipsoid is investigated as a function of distance to the plane and metallicity. We find no vertex deviation for old stars, consistent with an axisymmetric Galaxy. Paper II is devoted to the dynamical analysis of the sample, puting new constraints on the vertical force perpendicular to the galactic plane and on the total mass density in the galactic plane.

The velocity distribution of neighbouring stars is deduced from the Hipparcos proper motions. We have used a classical Schwarzschild decomposition and also developed a dynamical model for quasi-exponential stellar discs. This model is a 3-D derivation of Shu's model in the framework of St\"ackel potentials with three integrals of motion. We determine the solar motion relative to the local standard of rest (LSR), the density and kinematic radial gradients, as well as the local slope of the velocity curve. From the stellar disc dynamical model, we determine explicitly the link between the tangential-vertical velocity (v_\theta, v_z) coupling and the local shape of the potential. Using a restricted sample of 3-D velocity data, we measure $z_o$, the focus of the spheroidal coordinate system defining the best fitted St\"ackel potential. The parameter $z_o$ is related to the tilt of the velocity ellipsoid and more fundamentally to the mass gradient in the galactic disc. This parameter is found to be 5.7\pm1.4 kpc. This implies that the galactic potential is not extremly flat and that the dark matter component is not confined in the galactic plane.

This paper presents the detailed results obtained in the search of density- velocity inhomogeneities in a volume limited and absolute magnitude limited sample of A-F type dwarfs within 125 parsecs of the Sun. A 3-D wavelet analysis is used to extract inhomogeneities, both in the density and velocity distributions. Having established a real picture of the phase space without assumption we come back to previously known observational facts regarding clusters and associations, superclusters. In the 3-D position space, well known open clusters (Hyades, Coma Berenices and Ursa Major), associations (parts of the Scorpio-Centaurus association) as well as the Hyades evaporation track are retrieved. Three new probably loose clusters are identified (Bootes, Pegasus 1 and 2). The sample is relatively well mixed in the position space since less than 7 per cent of the stars belong to structures with coherent kinematics, most likely gravitationally bound. In the velocity space, the majority of large scale velocity structures ($\sigma$ ~ 6.3 $km s-1$) are Eggen's superclusters (Pleiades SCl, Hyades SCl and Sirius SCl) with the whole Centaurus association. A new supercluster-like structure is found with a mean velocity between the Sun and Sirius SCl velocities. These structures are all characterized by a large age range which reflects the overall sample age distribution. Moreover, a few old streams of ~ 2 Gyr are also extracted at this scale with high U components. We show that all these large velocity dispersion structures represent 46% of the sample. Smaller scales (\sigma ~ 3.8 and 2.4 $km s-1$) reveal that superclusters are always substructured by 2 or more streams which generally exhibit a coherent age distribution. Percentages of stars in these streams are 38% and 18% respectively.

A volume limited and absolute magnitude limited sample of A-F type dwarfs within 125 parsecs of the Sun is searched for inhomogeneities in the density- velocity space, expecting signatures of the cluster evaporation, phase mixing and possible disc heating mechanisms. A 3-D wavelet analysis is used to extract inhomogeneities, both in the density and velocity distributions. Thus, a real picture of the phase space is produced. Not only are some clusters and streams detected, but the fraction of clumped stars can be measured. By estimating individual stellar ages one can relate the streams and clusters to the state of the interstellar medium at star formation time and provide a quantitative view of cluster evaporation and stream mixing. As a result, we propose a coherent interpretation of superclusters and derive some quantitative evolutionary tracers which we expect to serve in the understanding of the large scale evolution of the galactic disc. The sample is relatively well mixed in the position space since less than 7 per cent of the stars are proper motion confirmed cluster members. We also detect star evaporation out of the Hyades open cluster. In the velocity space, the typical scale of so-called Eggen's superclusters ($\sigma 6.3 km/s$) does not seem to correspond to any physical entity. The picture they form, their frequency and their divisions at smaller scales are well compatible with their creation by chance coincidence of physically homogeneous smaller scale structures ($\sigma 3.8 or 2.4 km/s$). The detailed analysis of the phase space structures provides a scenario of kinematical evolution in the solar neighbourhood.

Hipparcos data provide the first, volume limited and absolute magnitude limited homogeneous tracer of stellar density and velocity distributions in the solar neighbourhood. The density of A-type stars more luminous than $M_v=2.5$ can be accurately mapped within a sphere of 125 pc radius, while proper motions in galactic latitude provide the vertical velocity distribution near the galactic plane. The potential well across the galactic plane is traced practically hypothesis-free and model-free. The local dynamical density comes out as $\rho_{0}=0.076 \pm0.015~M_{\sun}~{pc}^{-3}$ a value well below all previous determinations leaving no room for any disk shaped component of dark matter.

In a recent Letter, Cowsik, Ratnam and Bhattacharjee (1996a) have built a dynamically self-consistent spatial distribution of particles of galactic dark matter. They have come up with the rather unorthodox conclusion that the mean velocity dispersion of dark matter particles ``should be $600 kms$ or larger''. Their letter triggered immediate comments (Evans 1996, Gates et al. 1996). Here we find, as did Cowsik et al. (1996a), that models of the dark matter halo can be made consistent with velocity dispersion much larger than that expected from a simple application of the virial theorem in the solar neighbourhood. But in contrast to their conclusions, we show that using their model, we also obtain solutions with smaller velocity dispersion down to $\sim 270 kms$. These more orthodox dispersions arise because of less constraining boundary conditions for the central density but do not rely on indirect or model dependent measurements of the large scale distance behaviour of the rotation curve.

Accurate characterization of thick disc properties from recent kinematic and photometric surveys provides converging evidences that this intermediate population is a sequel of the violent heating of early disc populations by a merging satellite galaxy. The thick disc population is revisited under the light of new data in a number of galactic sample fields. Various thick disc hypotheses are fitted to observational data through a maximum likelihood technique. The resulting characteristics of the thick disc are the following : a scale height of 760 $\pm$ 50 pc, with a local density of 5.6 $\pm$ 1 \% of the thin disc. The scale length is constrained to be 2.8 $\pm$ 0.8 kpc, well in agreement with the disc scale length (2.5 $\pm$ 0.3 kpc). The mean metallicity of the thick disc is found to be -0.7 $\pm$ 0.2 dex, with no {significant} metallicity gradients. These photometric constraints in combination with kinematic data give new constraints on the thick disc formation. We show that thick disc characteristics are {hardly compatible with a top-down formation scenario} but fully compatible with a violent merging event arising at the early thin disc life time as described by Quinn, Hernquist \& Fullagar (1993).