The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes.

Significant investments to upgrade or construct large-scale scientific facilities demand commensurate investments in R&D to design algorithms and computing approaches to enable scientific and engineering breakthroughs in the big data era. The remarkable success of Artificial Intelligence (AI) algorithms to turn big-data challenges in industry and technology into transformational digital solutions that drive a multi-billion dollar industry, which play an ever increasing role shaping human social patterns, has promoted AI as the most sought after signal processing tool in big-data research. As AI continues to evolve into a computing tool endowed with statistical and mathematical rigor, and which encodes domain expertise to inform and inspire AI architectures and optimization algorithms, it has become apparent that single-GPU solutions for training, validation, and testing are no longer sufficient. This realization has been driving the confluence of AI and high performance computing (HPC) to reduce time-to-insight and to produce robust, reliable, trustworthy, and computationally efficient AI solutions. In this white paper, we present a summary of recent developments in this field, and discuss avenues to accelerate and streamline the use of HPC platforms to design accelerated AI algorithms.

In this paper we introduce a new method for analysing Milky Way phase-space which allows us to reveal the imprint left by the Milky Way bar and spiral arms on the stars with full phase-space data in Gaia Data Release 2. The unprecedented quality and extended spatial coverage of these data enable us to discover six prominent stellar density structures in the disc to a distance of 5 kpc from the Sun. Four of these structures correspond to the spiral arms detected previously in the gas and young stars (Scutum-Centaurus, Sagittarius, Local and Perseus). The remaining two are associated with the main resonances of the Milky Way bar where corotation is placed at around 6.2 kpc and the outer Lindblad resonance beyond the Solar radius, at around 9 kpc. For the first time we provide evidence of the imprint left by spiral arms and resonances in the stellar densities not relying on a specific tracer, through enhancing the signatures left by these asymmetries. Our method offers new avenues for studying how the stellar populations in our Galaxy are shaped.

Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravitational wave sources and their electromagnetic and astroparticle counterparts, and make a number of recommendations to maximize their potential for scientific discovery. These recommendations refer to the design of scalable and computationally efficient machine learning algorithms; the cyber-infrastructure to numerically simulate astrophysical sources, and to process and interpret multi-messenger astrophysics data; the management of gravitational wave detections to trigger real-time alerts for electromagnetic and astroparticle follow-ups; a vision to harness future developments of machine learning and cyber-infrastructure resources to cope with the big-data requirements; and the need to build a community of experts to realize the goals of multi-messenger astrophysics.

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

Previous studies based on the analysis of Gaia DR2 data have revealed that accreted stars, possibly originating from a single progenitor satellite, are a significant component of the halo of our Galaxy, potentially constituting most of the halo stars at $\rm [Fe/H] < -1$ within a few kpc from the Sun and beyond. In this paper, we couple astrometric data from Gaia DR2 with elemental abundances from APOGEE DR14 to characterize the kinematics and chemistry of in-situ and accreted populations up to $\rm [Fe/H] \sim -2$. Accreted stars appear to significantly impact the Galactic chemo-kinematic relations, not only at $\rm [Fe/H] < -1$, but also at metallicities typical of the thick and metal-poor thin discs. They constitute about 60% of all stars at $\rm [Fe/H] < -1$, the remaining 40% being made of (metal-weak) thick disc stars. We find that the stellar kinematic fossil record shows the imprint left by this accretion event which heated the old Galactic disc. We are able to age-date this kinematic imprint, showing that the accretion occurred between 9 and 11 Gyr ago, and that it led to the last significant heating of the Galactic disc. An important fraction of stars with abundances typical of the (metal-rich) thick disc, and heated by this interaction, is now found in the Galactic halo. Indeed about half of the kinematically defined halo at few kpc from the Sun is composed of metal-rich thick disc stars. Moreover, we suggest that this metal-rich thick disc component dominates the stellar halo of the inner Galaxy. The new picture that emerges from this study is one where the standard non-rotating in-situ halo population, the collapsed halo, seems to be more elusive than ever.

The period 2035-50 considered in the ESA Voyage long-term plan will coincide with a series of foreseeable advances in the characterization of the stellar content of the Milky Way. The Gaia mission, combined with large-scale spectroscopic surveys, is helping to build an unprecedented census in terms of the astrometric, kinematic and chemical properties of Galactic stellar populations. Within a decade, precise measurements of such properties will be available for hundreds of millions of stars. Meanwhile, time-domain surveys initiated with CoRoT and Kepler/K2 and carried on by space missions such as TESS and PLATO or ground-based projects like the LSST, will have brought asteroseismology to a high level of maturity. The combination of precise ages from asteroseismology with astrometric and spectroscopic data, on large stellar samples, is allowing Galactic archaeologists to gain new insight into the assembly history of the Milky Way. Recent breakthroughs --- based on the detection of solar-like oscillations in tens of thousands of red-giant stars --- demonstrate the potential of such approach. Therefore, we are convinced that an all-sky, high-cadence, long-duration stellar variability survey will become a scientific priority in the 2035-50 period. The Chronos concept presented here consists in a time-domain extension to Gaia. It will allow for mass and age estimates for half a million red giants within 1.7 kpc from the Sun and hence shed a new light on our understanding of the Galactic dynamics and archaeology. In terms of the targeted pulsators, Chronos will bridge the gap between PLATO and the LSST by surveying stars all the way from the subgiant branch to the early AGB. Finally, it will surpass all previous surveys capable of conducting asteroseismology in terms of the combined sky coverage and duration of the observations (2 x 3.75 months over the whole sky and >5 years in the CVZ).

We explore the origin of phase-space substructures revealed by the second Gaia data release in the disc of the Milky Way, such as the ridges in the $V_{\phi}$-$r$ plane, the undulations in the $V_{\phi}$-$r$-$V_r$ space and the streams in the $V_{\phi}$-$V_r$ plane. We use a collisionless N-body simulation with co-spatial thin and thick discs, along with orbit integration, to study the orbital structure close to the Outer Lindblad Resonance (OLR) of the bar. We find that a prominent, long-lived ridge is formed in the $V_{\phi}$-$r$ plane due to the OLR which translates to streams in the $V_{\phi}$-$V_r$ plane and examine which closed periodic and trapped librating orbits are responsible for these features. We find that orbits which carry out small librations around the $x_1(1)$ family are preferentially found at negative $V_r$, giving rise to a `horn'-like feature, while orbits with larger libration amplitudes, trapped around the $x_1(2)$ and $x_1(1)$ families, constitute the positive $V_r$ substructure, i.e. the Hercules-like feature. This changing libration amplitude of orbits will translate to a changing ratio of thin/thick disc stars, which could have implications on the metallicity distribution in this plane. We find that a scenario in which the Sun is placed close to the OLR gives rise to a strong asymmetry in $V_r$ in the $V_{\phi}$-$V_r$ plane (i.e. Hercules vs. `the horn') and subsequently to undulations in the $V_{\phi}$-$r$-$V_r$ space. We also explore a scenario in which the Sun is placed closer to the bar corotation and find that the bar perturbation $alone$ cannot give rise to the these features.

The importance of studying the Gaia DR2 parallax zero-point by external means was underlined by Lindegren et al. (2018), and initiated by several works making use of Cepheids, eclipsing binaries, and asteroseismology. Despite a very efficient elimination of basic-angle variations, a small fluctuation remains and shows up as a small offset in the Gaia DR2 parallaxes. By combining astrometric, asteroseismic, spectroscopic, and photometric constraints, we undertake a new analysis of the Gaia parallax offset for nearly 3000 red-giant branch (RGB) and 2200 red clump (RC) stars observed by Kepler, as well as about 500 and 700 red giants (both RGB and RC) selected by the K2 Galactic Archaeology Program in campaigns 3 and 6. Engaging into a thorough comparison of the astrometric and asteroseismic parallaxes, we are able to highlight the influence of the asteroseismic method, and measure parallax offsets in the Kepler field that are compatible with independent estimates from literature and open clusters. Moreover, adding the K2 fields to our investigation allows us to retrieve a clear illustration of the positional dependence of the zero-point, in general agreement with the information provided by quasars. Lastly, we initiate a two-step methodology to make progress in the simultaneous calibration of the asteroseismic scaling relations and of the Gaia DR2 parallax offset, which will greatly benefit from the gain in precision with the third Data Release of Gaia.

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.

Context. Open clusters are very good tracers of the evolution of the Galactic disc. Thanks to Gaia, their kinematics can be investigated with an unprecedented precision and accuracy. Aims. The distribution of open clusters in the 6D phase space is revisited with Gaia DR2. Methods. The weighted mean radial velocity of open clusters was determined, using the most probable members available from a previous astrometric investigation that also provided mean parallaxes and proper motions. Those parameters, all derived from Gaia DR2 only, were combined to provide the 6D phase space information of 861 clusters. The velocity distribution of nearby clusters was investigated, as well as the spatial and velocity distributions of the whole sample as a function of age. A high quality subsample was used to investigate some possible pairs and groups of clusters sharing the same Galactic position and velocity. Results. For the high quality sample that has 406 clusters, the median uncertainty of the weighted mean radial velocity is 0.5 km/s. The accuracy, assessed by comparison to ground-based high resolution spectroscopy, is better than 1 km/s. Open clusters nicely follow the velocity distribution of field stars in the close Solar neighbourhood previously revealed by Gaia DR2. As expected, the vertical distribution of young clusters is very flat but the novelty is the high precision to which this can be seen. The dispersion of vertical velocities of young clusters is at the level of 5 km/s. Clusters older than 1 Gyr span distances to the Galactic plane up to 1 kpc with a vertical velocity dispersion of 14 km/s, typical of the thin disc. Five pairs of clusters and one group with five members are possibly physically related. Other binary candidates previously identified turn out to be chance alignment.

Gaia data are revolutionizing our knowledge of the evolutionary history of the Milky Way. 3D maps of the interstellar dust provide complementary information and are a tool for a wide range of uses. We aimed at building 3D maps of the dust in the Local arm and surrounding regions. To do so, Gaia DR2 photometric data were combined with 2MASS measurements to derive extinction towards stars that possess accurate photometry and relative uncertainties on DR2 parallaxes smaller than 20\%. We applied to the extinctions a new hierarchical inversion algorithm adapted to large datasets and to a inhomogeneous target distribution. Each step associates regularized Bayesian inversions along radial directions and a subsequent inversion in 3D of their results. Each inverted distribution serves as a prior for the subsequent step and the spatial resolution is progressively increased. We present the resulting 3D distribution of the dust in a 6 x 6 x 0.8 kpc3 volume around the Sun. Its main features are found to be elongated along different directions that vary from below to above the mid-plane: the outer part of Carina-Sagittarius, mainly located above the mid-plane, the Local arm/Cygnus Rift around and above the mid-plane and the fragmented Perseus arm are oriented close to the direction of circular motion. The long spur (nicknamed the split) that extends between the Local Arm and Carina-Sagittarius, the compact near side of Carina-Sagittarius and the Cygnus Rift below the Plane are oriented along l=40 to 55 deg. Dust density images in vertical planes reveal in some regions a wavy pattern and show that the solar neighborhood within 500 pc remains atypical by its extent above and below the Plane. We show several comparisons with the locations of molecular clouds, HII regions, O stars and masers. The link between the dust concentration and these tracers is markedly different from one region to the other.

The advent of experimental science facilities-instruments and observatories, such as the Large Hadron Collider, the Laser Interferometer Gravitational Wave Observatory, and the upcoming Large Synoptic Survey Telescope-has brought about challenging, large-scale computational and data processing requirements. Traditionally, the computing infrastructure to support these facility's requirements were organized into separate infrastructure that supported their high-throughput needs and those that supported their high-performance computing needs. We argue that to enable and accelerate scientific discovery at the scale and sophistication that is now needed, this separation between high-performance computing and high-throughput computing must be bridged and an integrated, unified infrastructure provided. In this paper, we discuss several case studies where such infrastructure has been implemented. These case studies span different science domains, software systems, and application requirements as well as levels of sustainability. A further aim of this paper is to provide a basis to determine the common characteristics and requirements of such infrastructure, as well as to begin a discussion of how best to support the computing requirements of existing and future experimental science facilities.

This report provides an overview of recent work that harnesses the Big Data Revolution and Large Scale Computing to address grand computational challenges in Multi-Messenger Astrophysics, with a particular emphasis on real-time discovery campaigns. Acknowledging the transdisciplinary nature of Multi-Messenger Astrophysics, this document has been prepared by members of the physics, astronomy, computer science, data science, software and cyberinfrastructure communities who attended the NSF-, DOE- and NVIDIA-funded "Deep Learning for Multi-Messenger Astrophysics: Real-time Discovery at Scale" workshop, hosted at the National Center for Supercomputing Applications, October 17-19, 2018. Highlights of this report include unanimous agreement that it is critical to accelerate the development and deployment of novel, signal-processing algorithms that use the synergy between artificial intelligence (AI) and high performance computing to maximize the potential for scientific discovery with Multi-Messenger Astrophysics. We discuss key aspects to realize this endeavor, namely (i) the design and exploitation of scalable and computationally efficient AI algorithms for Multi-Messenger Astrophysics; (ii) cyberinfrastructure requirements to numerically simulate astrophysical sources, and to process and interpret Multi-Messenger Astrophysics data; (iii) management of gravitational wave detections and triggers to enable electromagnetic and astro-particle follow-ups; (iv) a vision to harness future developments of machine and deep learning and cyberinfrastructure resources to cope with the scale of discovery in the Big Data Era; (v) and the need to build a community that brings domain experts together with data scientists on equal footing to maximize and accelerate discovery in the nascent field of Multi-Messenger Astrophysics.

Using a single N-body simulation ($N=0.14\times 10^9$) we explore the formation, evolution and spatial variation of the phase-space spirals similar to those recently discovered by Antoja et al. in the Milky Way disk, with Gaia DR2. For the first time in the literature, we use a self-consistent N-body simulation of an isolated Milky Way-type galaxy to show that the phase-space spirals develop naturally from vertical oscillations driven by the buckling of the stellar bar. We claim that the physical mechanism standing behind the observed incomplete phase-space mixing process can be internal and not necessarily due to the perturbation induced by a massive satellite. In our model, the bending oscillations propagate outwards and produce axisymmetric variations of the mean vertical coordinate and of the vertical velocity component. As a consequence, the phase-space wrapping results in the formation of patterns with various morphology across the disk, depending on the bar orientation, distance to the galactic center and time elapsed since the bar buckling. Once bending waves appear, they are supported for a long time via disk self-gravity. The underlying physical mechanism implies the link between in-plane and vertical motion that leads directly to phase-space structures whose amplitude and shape are in remarkable agreement with those of the phase-space spirals observed in the Milky Way disk. In our isolated galaxy simulation, phase-space spirals are still distinguishable, at the solar neighbourhood, 3 Gyr after the buckling phase. The long-lived character of the phase-space spirals generated by the bar buckling instability cast doubts on the timing argument used so far to get back at the time of the onset of the perturbation: phase-space spirals may have been caused by perturbations originated several Gyrs ago, and not as recent as suggested so far.

The evolution of the Milky Way disk, which contains most of the stars in the Galaxy, is affected by several phenomena. For example, the bar and the spiral arms of the Milky Way induce radial migration of stars and can trap or scatter stars close to orbital resonances. External perturbations from satellite galaxies can also have a role, causing dynamical heating of the Galaxy, ring-like structures in the disk and correlations between different components of the stellar velocity. These perturbations can also cause 'phase wrapping' signatures in the disk, such as arched velocity structures in the motions of stars in the Galactic plane. Some manifestations of these dynamical processes have already been detected, including kinematic substructure in samples of nearby stars, density asymmetries and velocities across the Galactic disk that differ from the axisymmetric and equilibrium expectations, especially in the vertical direction, and signatures of incomplete phase mixing in the disk. Here we report an analysis of the motions of six million stars in the Milky Way disk. We show that the phase-space distribution contains different substructures with various morphologies, such as snail shells and ridges, when spatial and velocity coordinates are combined. We infer that the disk must have been perturbed between 300 million and 900 million years ago, consistent with estimates of the previous pericentric passage of the Sagittarius dwarf galaxy. Our findings show that the Galactic disk is dynamically young and that modelling it as time-independent and axisymmetric is incorrect.

Using Gaia DR2 astrometry, we map the kinematic signature of the Galactic stellar warp out to a distance of 7 kpc from the Sun. Combining Gaia DR2 and 2MASS photometry, we identify, via a probabilistic approach, 599 494 upper main sequence stars and 12 616 068 giants without the need for individual extinction estimates. The spatial distribution of the upper main sequence stars clearly shows segments of the nearest spiral arms. The large-scale kinematics of both the upper main sequence and giant populations show a clear signature of the warp of the Milky Way, apparent as a gradient of 5-6 km/s in the vertical velocities from 8 to 14 kpc in Galactic radius. The presence of the signal in both samples, which have different typical ages, suggests that the warp is a gravitationally induced phenomenon.

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.

This paper presents the specification, design, and development of the Radial Velocity Spectrometer (RVS) on the European Space Agency's Gaia mission. Starting with the rationale for the full six dimensions of phase space in the dynamical modelling of the Galaxy, the scientific goals and derived top-level instrument requirements are discussed, leading to a brief description of the initial concepts for the instrument. The main part of the paper is a description of the flight RVS, considering the optical design, the focal plane, the detection and acquisition chain, and the as-built performance drivers and critical technical areas. After presenting the pre-launch performance predictions, the paper concludes with the post-launch developments and mitigation strategies, together with a summary of the in-flight performance at the end of commissioning.

The Gaia Data Release 2 contains the 1st release of radial velocities complementing the kinematic data of a sample of about 7 million relatively bright, late-type stars. Aims: This paper provides a detailed description of the Gaia spectroscopic data processing pipeline, and of the approach adopted to derive the radial velocities presented in DR2. Methods: The pipeline must perform four main tasks: (i) clean and reduce the spectra observed with the Radial Velocity Spectrometer (RVS); (ii) calibrate the RVS instrument, including wavelength, straylight, line-spread function, bias non-uniformity, and photometric zeropoint; (iii) extract the radial velocities; and (iv) verify the accuracy and precision of the results. The radial velocity of a star is obtained through a fit of the RVS spectrum relative to an appropriate synthetic template spectrum. An additional task of the spectroscopic pipeline was to provide 1st-order estimates of the stellar atmospheric parameters required to select such template spectra. We describe the pipeline features and present the detailed calibration algorithms and software solutions we used to produce the radial velocities published in DR2. Results: The spectroscopic processing pipeline produced median radial velocities for Gaia stars with narrow-band near-IR magnitude Grvs < 12 (i.e. brighter than V~13). Stars identified as double-lined spectroscopic binaries were removed from the pipeline, while variable stars, single-lined, and non-detected double-lined spectroscopic binaries were treated as single stars. The scatter in radial velocity among different observations of a same star, also published in DR2, provides information about radial velocity variability. For the hottest (Teff > 7000 K) and coolest (Teff < 3500 K) stars, the accuracy and precision of the stellar parameter estimates are not sufficient to allow selection of appropriate templates. [Abridged]

To illustrate the potential of GDR2, we provide a first look at the kinematics of the Milky Way disc, within a radius of several kiloparsecs around the Sun. We benefit for the first time from a sample of 6.4 million F-G-K stars with full 6D phase-space coordinates, precise parallaxes, and precise Galactic cylindrical velocities . From this sample, we extracted a sub-sample of 3.2 million giant stars to map the velocity field of the Galactic disc from $\sim$5~kpc to $\sim$13~kpc from the Galactic centre and up to 2~kpc above and below the plane. We also study the distribution of 0.3 million solar neighbourhood stars ($r < 200$~pc), with median velocity uncertainties of 0.4~km/s, in velocity space and use the full sample to examine how the over-densities evolve in more distant regions. GDR2 allows us to draw 3D maps of the Galactocentric median velocities and velocity dispersions with unprecedented accuracy, precision, and spatial resolution. The maps show the complexity and richness of the velocity field of the galactic disc. We observe streaming motions in all the components of the velocities as well as patterns in the velocity dispersions. For example, we confirm the previously reported negative and positive galactocentric radial velocity gradients in the inner and outer disc, respectively. Here, we see them as part of a non-axisymmetric kinematic oscillation, and we map its azimuthal and vertical behaviour. We also witness a new global arrangement of stars in the velocity plane of the solar neighbourhood and in distant regions in which stars are organised in thin substructures with the shape of circular arches that are oriented approximately along the horizontal direction in the $U-V$ plane. Moreover, in distant regions, we see variations in the velocity substructures more clearly than ever before, in particular, variations in the velocity of the Hercules stream. (abridged)

The European Space Agency Gaia satellite was launched into orbit around L2 in December 2013. This ambitious mission has strict requirements on residual systematic errors resulting from instrumental corrections in order to meet a design goal of sub-10 microarcsecond astrometry. During the design and build phase of the science instruments, various critical calibrations were studied in detail to ensure that this goal could be met in orbit. In particular, it was determined that the video-chain offsets on the analogue side of the analogue-to-digital conversion electronics exhibited instabilities that could not be mitigated fully by modifications to the flight hardware. We provide a detailed description of the behaviour of the electronic offset levels on microsecond timescales, identifying various systematic effects that are known collectively as offset non-uniformities. The effects manifest themselves as transient perturbations on the gross zero-point electronic offset level that is routinely monitored as part of the overall calibration process. Using in-orbit special calibration sequences along with simple parametric models, we show how the effects can be calibrated, and how these calibrations are applied to the science data. While the calibration part of the process is relatively straightforward, the application of the calibrations during science data processing requires a detailed on-ground reconstruction of the readout timing of each charge-coupled device (CCD) sample on each device in order to predict correctly the highly time-dependent nature of the corrections. We demonstrate the effectiveness of our offset non-uniformity models in mitigating the effects in Gaia data. We demonstrate for all CCDs and operating instrument and modes on board Gaia that the video-chain noise-limited performance is recovered in the vast majority of science samples.

Aims. The Radial Velocity Spectrometer (RVS) on board the ESA satellite mission Gaia has no calibration device. Therefore, the radial velocity zero point needs to be calibrated with stars that are proved to be stable at a level of 300 m/s during the Gaia observations. Methods. We compiled a dataset of ~71000 radial velocity measurements from five high-resolution spectrographs. A catalogue of 4813 stars was built by combining these individual measurements. The zero point was established using asteroids. Results. The resulting catalogue has seven observations per star on average on a typical time baseline of six years, with a median standard deviation of 15 m/s. A subset of the most stable stars fulfilling the RVS requirements was used to establish the radial velocity zero point provided in Gaia Data Release 2. The stars that were not used for calibration are used to validate the RVS data.

For Gaia DR2 (GDR2), 280 million spectra, collected by the RVS instrument on-board Gaia, were processed and median radial velocities were derived for 9.8 million sources brighter than Grvs = 12 mag. This paper describes the validation and properties of the median radial velocities published in GDR2. Quality tests and filters are applied to select, from the 9.8 million radial velocities, those with the quality to be published in GDR2. The accuracy of the selected sample is assessed with respect to ground-based catalogues. Its precision is estimated using both ground-based catalogues and the distribution of the Gaia radial velocity uncertainties. GDR2 contains median radial velocities for 7 224 631 stars, with Teff in the range [3550, 6900] K, which passed succesfully the quality tests. The published median radial velocities provide a full sky-coverage and have a completness with respect to the astrometric data of 77.2\% (for $G \leq 12.5$ mag). The median radial velocity residuals with respect to the ground-based surveys vary from one catalogue to another, but do not exceed a few 100s m/s. In addition, the Gaia radial velocities show a positive trend as a function of magnitude, which starts around Grvs $\sim 9$ mag and reaches about $+500$ m/s at Grvs $= 11.75$ mag. The overall precision, estimated from the median of the Gaia radial velocity uncertainties, is 1.05 km/s. The radial velocity precision is function of many parameters, in particular the magnitude and effective temperature. For bright stars, Grvs in [4, 8] mag, the precision is in the range 200-350 m/s, which is about 3 to 5 times more precise than the pre-launch specification of 1 km/s. At the faint end, Grvs = 11.75 mag, the precisions for Teff = 5000 K and 6500 K are respectively 1.4 km/s and 3.7 km/s.

The second Gaia data release (DR2) is scheduled for April 2018. While Gaia DR1 had increased the number of stars with parallaxes by a factor 20 with respect to the Hipparcos catalogue, Gaia DR2 will bring another factor 500 increase, with parallaxes (and proper motions) for more than a billion stars. In addition, Gaia DR2 will deliver improved accuracy and precision for the astrometric and photometric data, $G$, $G_\mathrm{BP}$, $G_\mathrm{RP}$ magnitudes, radial velocities, identification and characterisation of variable stars and asteroids as well as stellar parameters for stars down to $G = 17$ mag. On behalf of the teams of the Gaia-DPAC consortium, these proceedings give a foretaste of Gaia DR2, 6 months before the release.

Deciphering the assembly history of the Milky Way is a formidable task, which becomes possible only if one can produce high-resolution chrono-chemo-kinematical maps of the Galaxy. Data from large-scale astrometric and spectroscopic surveys will soon provide us with a well-defined view of the current chemo-kinematical structure of the Milky Way, but will only enable a blurred view on the temporal sequence that led to the present-day Galaxy. As demonstrated by the (ongoing) exploitation of data from the pioneering photometric missions CoRoT, Kepler, and K2, asteroseismology provides the way forward: solar-like oscillating giants are excellent evolutionary clocks thanks to the availability of seismic constraints on their mass and to the tight age-initial-mass relation they adhere to. In this paper we identify five key outstanding questions relating to the formation and evolution of the Milky Way that will need precise and accurate ages for large samples of stars to be addressed, and we identify the requirements in terms of number of targets and the precision on the stellar properties that are needed to tackle such questions. By quantifying the asteroseismic yields expected from PLATO for red-giant stars, we demonstrate that these requirements are within the capabilities of the current instrument design, provided that observations are sufficiently long to identify the evolutionary state and allow robust and precise determination of acoustic-mode frequencies. This will allow us to harvest data of sufficient quality to reach a 10% precision in age. This is a fundamental pre-requisite to then reach the more ambitious goal of a similar level of accuracy, which will only be possible if we have to hand a careful appraisal of systematic uncertainties on age deriving from our limited understanding of stellar physics, a goal which conveniently falls within the main aims of PLATO's core science.

We explore trends in the morphology and line-of-sight (los) velocity of stellar populations in the inner regions of disc galaxies, using N-body simulations with both a thin (kinematically cold) and a thick (kinematically hot) disc which form a bar and boxy/peanut (b/p) bulge. The bar in the thin disc component is $\sim$50\% stronger than the thick disc bar and is more elongated, with an axis ratio almost half that of the thick disc bar. The thin disc b/p bulge has a pronounced X-shape, while the thick disc b/p is weaker with a rather boxy shape. This leads to the signature of the b/p bulge in the thick disc to be weaker and further away from the plane than in the thin disc. Regarding the kinematics, we find that the los velocity of thick disc stars in the outer parts of the b/p bulge can be \emph{larger} than that of thin disc stars, by up to 40\% and 20\% for side-on and Milky Way-like orientations of the bar respectively. This is due to the different orbits followed by thin and thick disc stars in the bar-b/p region, which are affected by the fact that: i) thin disc stars are trapped more efficiently in the bar - b/p instability and thus lose more angular momentum than their thick disc counterparts and ii) thick disc stars have large radial excursions and therefore stars from large radii with high angular momenta can be found in the bar region. We also find that the difference between the los velocities of the thin and thick disc in the b/p bulge ($\Delta v_{los}$) correlates with the initial difference between the radial velocity dispersions of the two discs ($\Delta \sigma$) . We therefore conclude that stars in the bar - b/p bulge will have considerably different morphologies and kinematics depending on the kinematic properties of the disc population they originate from.

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.

Aims: A previous study claimed that the star HIP 21539 passed close to the Sun, at a distance of 1.9 pc, around 0.14 Myr ago. We show that this is not the case. Methods: We redetermined the trajectory of the star relative to the Sun using a new accurate radial velocity from the HARPS spectrograph combined with the recent Gaia-TGAS astrometry. Results: With this new data, the closest approach of HIP 21539 to the Sun is now 17 pc, instead of 1.9 pc. Conclusions: At this distance, the star has not perturbed the Oort cloud.

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.

Gaia is a space mission currently measuring the five astrometric parameters as well as spectrophotometry of at least 1 billion stars to G = 20.7 mag with unprecedented precision. The sixth parameter in phase space (radial velocity) is also measured thanks to medium-resolution spectroscopy being obtained for the 150 million brightest stars. During the commissioning phase, two fields, one around each ecliptic pole, have been repeatedly observed to assess and to improve the overall satellite performances as well as the associated reduction and analysis software. A ground-based photometric and spectroscopic survey was therefore initiated in 2007, and is still running in order to gather as much information as possible about the stars in these fields. This work is of particular interest to the validation of the Radial Velocity Spectrometer (RVS) outputs. The paper presents the radial velocity measurements performed for the Southern targets in the 12 - 17 R magnitude range on high- to mid-resolution spectra obtained with the GIRAFFE and UVES spectrographs.

Gaia's Radial Velocity Spectrometer (RVS) has been operating in routine phase for over one year since initial commissioning. RVS continues to work well but the higher than expected levels of straylight reduce the limiting magnitude. The end-of-mission radial-velocity (RV) performance requirement for G2V stars was 15 km/s at V = 16.5 mag. Instead, 15 km/s precision is achieved at 15 < V < 16 mag, consistent with simulations that predict a loss of 1.4 mag. Simulations also suggest that changes to Gaia's onboard software could recover ~0.14 mag of this loss. Consequently Gaia's onboard software was upgraded in April 2015. The status of this new commissioning period is presented, as well as the latest scientific performance of the on-ground processing of RVS spectra. We illustrate the implications of the RVS limiting magnitude on Gaia's view of the Milky Way's halo in 6D using the Gaia Universe Model Snapshot (GUMS).

Quenching, the cessation of star formation, is one of the most significant events in the life cycle of galaxies. We show here the first evidence that the Milky Way experienced a generalised quenching of its star formation at the end of its thick disk formation $\sim$9 Gyr ago. Elemental abundances of stars studied as part of the APOGEE survey reveal indeed that in less than $\sim$2 Gyr the star formation rate in our Galaxy dropped by an order-of-magnitude. Because of the tight correlation between age and alpha abundance, this event reflects in the dearth of stars along the inner disk sequence in the [Fe/H]-[$\alpha$/Fe] plane. Before this phase, which lasted about 1.5 Gyr, the Milky Way was actively forming stars. Afterwards, the star formation resumed at a much lower level to form the thin disk. These events are very well matched by the latest observation of MW-type progenitors at high redshifts. In late type galaxies, quenching is believed to be related to a long and secular exhaustion of gas. In our Galaxy, it occurred on a much shorter time scale, while the chemical continuity before and after the quenching indicates that it was not due to the exhaustion of the gas. While quenching is generally associated with spheroids, our results show that it also occurs in galaxies like the Milky Way, possibly when they are undergoing a morphological transition from thick to thin disks. Given the demographics of late type galaxies in the local universe, in which classical bulges are rare, we suggest further that this may hold true generally in galaxies with mass lower than or approximately $M^*$, where quenching could be directly a consequence of thick disk formation. We emphasize that the quenching phase in the Milky Way could be contemporaneous with, and related to, the formation of the bar. We sketch a scenario on how a strong bar may inhibit star formation.

(Abridged) By using an N-body simulation of a bulge that was formed via a bar instability mechanism, we analyse the imprints of the initial (i.e. before bar formation) location of stars on the bulge kinematics, in particular on the heliocentric radial velocity distribution of bulge stars. Four different latitudes were considered: $b=-4^\circ$, $-6^\circ$, $-8^\circ$, and $-10^\circ$, along the bulge minor axis as well as outside it, at $l=\pm5^\circ$ and $l=\pm10^\circ$. The bulge X-shaped structure comprises stars that formed in the disk at different locations. Stars formed in the outer disk, beyond the end of the bar, which are part of the boxy peanut-bulge structure may show peaks in the velocity distributions at positive and negative heliocentric radial velocities with high absolute values that can be larger than 100 $\rm km$ $\rm s^{-1}$, depending on the observed direction. In some cases the structure of the velocity field is more complex and several peaks are observed. Stars formed in the inner disk, the most numerous, contribute predominantly to the X-shaped structure and present different kinematic characteristics. Our results may enable us to interpret the cold high-velocity peak observed in the APOGEE commissioning data, as well as the excess of high-velocity stars in the near and far arms of the X-shaped structure at $l$=$0^\circ$ and $b$=$-6^\circ$. When compared with real data, the kinematic picture becomes more complex due to the possible presence in the observed samples of classical bulge and/or thick disk stars. Overall, our results point to the existence of complex patterns and structures in the bulge velocity fields, which are generated by the bar. This suggests that caution should be used when interpreting the bulge kinematics: the presence of substructures, peaks and clumps in the velocity fields is not necessarily a sign of past accretion events.

The sum of active neutrino masses is well constrained, $58$ meV $\leq m_\nu \lesssim 0.23$ eV, but the origin of this scale is not well understood. Here we investigate the possibility that it arises by environmental selection in a large landscape of vacua. Earlier work had noted the detrimental effects of neutrinos on large scale structure. However, using Boltzmann codes to compute the smoothed density contrast on Mpc scales, we find that dark matter halos form abundantly for $m_\nu \gtrsim 10$ eV. This finding rules out an anthropic origin of $m_\nu$, unless a different catastrophic boundary can be identified. Here we argue that galaxy formation becomes inefficient for $m_\nu \gtrsim 10$ eV. We show that in this regime, structure forms late and is dominated by cluster scales, as in a top-down scenario. This is catastrophic: baryonic gas will cool too slowly to form stars in an abundance comparable to our universe. With this novel cooling boundary, we find that the anthropic prediction for $m_\nu$ agrees at better than $2\sigma$ with current observational bounds. A degenerate hierarchy is mildly preferred.

We develop a chemical evolution model in order to study the star formation history of the Milky Way. Our model assumes that the Milky Way is formed from a closed box-like system in the inner regions, while the outer parts of the disc experience some accretion. Unlike the usual procedure, we do not fix the star formation prescription (e.g. Kennicutt law) in order to reproduce the chemical abundance trends. Instead, we fit the abundance trends with age in order to recover the star formation history of the Galaxy. Our method enables one to recover with unprecedented accuracy the star formation history of the Milky Way in the first Gyrs, in both the inner (R<7-8kpc) and outer (R>9-10kpc) discs as sampled in the solar vicinity. We show that, in the inner disc, half of the stellar mass formed during the thick disc phase, in the first 4-5 Gyr. This phase was followed by a significant dip in the star formation activity (at 8-9 Gyr) and a period of roughly constant lower level star formation for the remaining 8 Gyr. The thick disc phase has produced as many metals in 4 Gyr as the thin disc in the remaining 8 Gyr. Our results suggest that a closed box model is able to fit all the available constraints in the inner disc. A closed box system is qualitatively equivalent to a regime where the accretion rate, at high redshift, maintains a high gas fraction in the inner disc. In such conditions, the SFR is mainly governed by the high turbulence of the ISM. By z~1 it is possible that most of the accretion takes place in the outer disc, while the star formation activity in the inner disc is mostly sustained by the gas not consumed during the thick disc phase, and the continuous ejecta from earlier generations of stars. The outer disc follows a star formation history very similar to that of the inner disc, although initiated at z~2, about 2 Gyr before the onset of the thin disc formation in the inner disc.

By analyzing a N-body simulation of a bulge formed simply via a bar instability mechanism operating on a kinematically cold stellar disk, and by comparing the results of this analysis with the structural and kinematic properties of the main stellar populations of the Milky Way bulge, we conclude that the bulge of our Galaxy is not a pure stellar bar formed from a pre-existing thin stellar disk, as some studies have recently suggested. On the basis of several arguments emphasized in this paper, we propose that the bulge population which, in the Milky Way, is observed not to be part of the peanut structure corresponds to the old galactic thick disk, thus implying that the Milky Way is a pure thin+thick disk galaxy, with only a possible limited contribution of a classical bulge.

In 2010 Spergel introduced an alternative to the traditional Sersic form for galactic luminosity profiles based on modified Bessel functions of the second kind. His motivation was the desire for an accurate one-parameter profile form with a simple Fourier transform (in contrast to the Fourier transform of the Sersic profile which can't be written in closed form), but we have found that the Spergel profile almost universally makes integrals easier when it replaces the Sersic profile in the integrand. In the original paper on the subject Spergel noted that his profile seems to fit galaxies on average just as well as Sersic's. Here we make this observation quantitative by comparing the residuals from fitting Sersic and Spergel forms to data. We find that the Spergel profile actually fits better than the Sersic for a random sample of 16 galaxies.

(abridged) We have observed almost 1/3 of the globular clusters in the Milky Way, targeting distant and/or highly reddened objects, besides a few reference clusters. A large sample of red giant stars was observed with FORS2@VLT/ESO at R ~ 2,000. The method for derivation of stellar parameters is presented with application to six reference clusters. We aim at deriving the stellar parameters effective temperature, gravity, metallicity and alpha-element enhancement, as well as radial velocity, for membership confirmation of individual stars in each cluster. We analyse the spectra collected for the reference globular clusters NGC 6528, NGC 6553, M 71, NGC 6558, NGC 6426 and Terzan 8. They cover the full range of globular cluster metallicities, and are located in the bulge, disc and halo. Full spectrum fitting techniques are applied, by comparing each target spectrum with a stellar library in the optical region at 4560-5860 A. We employed the library of observed spectra MILES, and the synthetic library by Coelho et al. (2005). Validation of the method is achieved through recovery of the known atmospheric parameters for 49 well-studied stars that cover a wide range in the parameter space. We adopted as final stellar parameters (effective temperatures, gravities, metallicities) the average of results using MILES and Coelho et al. libraries. We identified 4 member stars in NGC 6528, 13 in NGC 6553, 10 in M 71, 5 in NGC 6558, 5 in NGC 6426 and 12 in Terzan 8. Radial velocities, Teff, log(g), [Fe/H] and alpha-element enhancements were derived. We derived abundances for NGC 6426 from spectroscopy for the first time. The method proved to be reliable for red giant stars observed with resolution R ~ 2,000, yielding results compatible with high-resolution spectroscopy. The derived alpha-element abundances show [A/Fe] vs. [Fe/H] consistent with that of field stars at the same metallicities.

The Gaia payload ensures maximum passive stability using a single material, SiC, for most of its elements. Dedicated metrology instruments are, however, required to carry out two functions: monitoring the basic angle and refocusing the telescope. Two interferometers fed by the same laser are used to measure the basic angle changes at the level of $\mu$as (prad, micropixel), which is the highest level ever achieved in space. Two Shack-Hartmann wavefront sensors, combined with an ad-hoc analysis of the scientific data are used to define and reach the overall best-focus. In this contribution, the systems, data analysis, procedures and performance achieved during commissioning are presented

The derivation of radial velocities from large numbers of spectra that typically result from survey work, requires automation. However, except for the classical cases of slowly rotating late-type spectra, existing methods of measuring Doppler shifts require fine-tuning to avoid a loss of accuracy due to the idiosyncrasies of individual spectra. The radial velocity spectrometer (RVS) on the Gaia mission, which will start operating very soon, prompted a new attempt at creating a measurement pipeline to handle a wide variety of spectral types. The present paper describes the theoretical background on which this software is based. However, apart from the assumption that only synthetic templates are used, we do not rely on any of the characteristics of this instrument, so our results should be relevant for most telescope-detector combinations. We propose an approach based on the simultaneous use of several alternative measurement methods, each having its own merits and drawbacks, and conveying the spectral information in a different way, leading to different values for the measurement. A comparison or a combination of the various results either leads to a "best estimate" or indicates to the user that the observed spectrum is problematic and should be analysed manually. We selected three methods and analysed the relationships and differences between them from a unified point of view; with each method an appropriate estimator for the individual random error is chosen. We also develop a procedure for tackling the problem of template mismatch in a systematic way. Furthermore, we propose several tests for studying and comparing the performance of the various methods as a function of the atmospheric parameters of the observed objects. Finally, we describe a procedure for obtaining a knowledge-based combination of the various Doppler-shift measurements.

By means of idealized, dissipationless N-body simulations which follow the formation and subsequent buckling of a stellar bar, we study the characteristics of boxy/peanut-shaped bulges and compare them with the properties of the stellar populations in the Milky Way bulge. The main results of our modeling, valid for the general family of boxy/peanut shaped bulges, are the following: (i) because of the spatial redistribution in the disk initiated at the epoch of bar formation, stars from the innermost regions to the outer Lindblad resonance of the stellar bar are mapped into a boxy bulge; (ii) the contribution of stars to the local bulge density depends on their birth radius: stars born in the innermost disk tend to dominate the innermost regions of the boxy bulge, while stars originating closer to the OLR are preferentially found in the outer regions of the boxy/peanut structure; (iii) stellar birth radii are imprinted in the bulge kinematics, the larger the birth radii of stars ending up in the bulge, the greater their rotational support and the higher their line-of- sight velocity dispersions (but note that this last trend depends on the bar viewing angle); (iv) the higher the classical bulge-over-disk ratio, the larger its fractional contribution of stars at large vertical distance from the galaxy mid-plane. (ABRIDGED) On the basis of their chemical and kinematic characteristics, the results of our modeling suggests that the populations A, B and C, as defined by the ARGOS survey, can be associated, respectively, with the inner thin disk, to the young thick and to the old thick disk, following the nomenclature recently suggested for stars in the solar neighborhood by Haywood et al. (2013).

Constraints on the Galactic bulge/bar structure and formation history from stellar kinematics and metallicities mainly come from relatively high-latitude fields (|b|>4) where a complex mix of stellar population is seen. We aim here to constrain the formation history of the Galactic bar by studying the radial velocity and metallicity distributions of stars in-situ (|b|<1). We observed red clump stars in four fields along the bar's major axis (l=10,6,-6 and b=0 plus a field at l=0,b=1) with low-resolution spectroscopy from VLT/FLAMES, observing around the CaII triplet. We developed robust methods for extracting radial velocity and metallicity estimates from these low signal-to-noise spectra. We derived distance probability distributions using Bayesian methods rigorously handling the extinction law. We present radial velocities and metallicity distributions, as well as radial velocity trends with distance. We observe an increase in the radial velocity dispersion near the Galactic plane. We detect the streaming motion of the stars induced by the bar in fields at l=+/-6, the highest velocity components of this bar stream being metal-rich ([Fe/H]~0.2 dex). Our data is consistent with a bar inclined at 26+/-3 from the Sun-Galactic centre line. We observe a significant fraction of metal-poor stars, in particular in the field at l=0,b=1. We confirm the flattening of the metallicity gradient along the minor axis when getting closer to the plane, with a hint that it could actually be inverted. Our stellar kinematics corresponds to the expected behaviour of a bar issued from the secular evolution of the Galactic disc. The mix of several populations, seen further away from the plane, is also seen in the bar in-situ since our metallicity distributions highlight a different spatial distribution between metal-poor and metal-rich stars, the more metal-poor stars being more centrally concentrated.

We report the first robust measurement of the Milky Way star formation history using the imprint left on chemical abundances of long-lived stars. The formation of the Galactic thick disc occurs during an intense star formation phase between 9.0 (z~1.5) and 12.5 Gyr (z~4.5) ago and is followed by a dip (at z~1.1) lasting about 1 Gyr. Our results imply that the thick disc is as massive as the Milky Way's thin disc, suggesting a fundamental role of this component in the genesis of our Galaxy, something that had been largely unrecognized. This new picture implies that huge quantities of gas necessary to feed the building of the thick disc must have been present at these epochs, in contradiction with the long-term infall assumed by chemical evolution models in the last two decades. These results allow us to fit the Milky Way within the emerging features of the evolution of disc galaxies in the early Universe.

We analyze high quality abundances data of solar neighborhood stars and show that there are two distinct regimes of [alpha/Fe] versus age which we identify as the epochs of the thick and thin disk formation. A tight correlation between metallicity and [alpha/Fe] versus age is clearly identifiable on thick disk stars, implying that this population formed from a well mixed ISM, over a time scale of 4-5 Gyr. Thick disk stars vertical velocity dispersion correlate with age, with the youngest objects having as small scale heights as those of thin disk stars. A natural consequence of these two results is that a vertical metallicity gradient is expected in this population. We suggest that the thick disk set the initial conditions for the formation of the inner thin disk. This provides also an explanation of the apparent coincidence between the step in metallicity at 7-10 kpc in the thin disk and the confinment of the thick disk at about R<10 kpc. We suggest that the outer thin disk developped outside the influence of the thick disk, but also that the high alpha-enrichment of the outer regions may originate from a primordial pollution by the gas expelled from the thick disk. Local metal-poor thin disk stars, whose properties are best explained by an origin in the outer disk, are shown to be as old as the youngest thick disk (9-10 Gyr), implying that the outer thin disk started to form while the thick disk formation was still on-going in the inner Galaxy. We point out that, given the tight age-abundance relations in the thick disk, an inside-out process would give rise to a radial gradient in abundances in this population which is not observed. Finally, we argue that the data discussed here leave little room for radial migration, either to have contaminated the solar vicinity, or, to have redistributed stars in significant proportion across the solar annulus.

The Radial Velocity Spectrograph (RVS) on board of Gaia needs to be calibrated using stable reference stars known in advance. The catalogue presented here has being built for that purpose. It includes 1420 radial velocity standard star candidates selected on strict criteria in order to fulfill the Gaia-RVS requirements. A large programme of ground based observations is underway since 2006 to monitor these stars and verify their stability which has to be better than 300 m/s over several years. The observations have been made on the echelle spectrographs ELODIE and SOPHIE on the 1.93-m telescope at Observatoire de Haute-Provence (OHP), NARVAL on the Telescope Bernard Lyot at Observatoire du Pic du Midi and CORALIE on the Euler-Swiss Telescope at La Silla. Data from the OHP and Geneva Observatory archives have also been retrieved as well as HARPS spectra from the ESO archive. We provide for each star a mean radial velocity in the SOPHIE scale, derived from the combination of velocities measured with those instruments, after having carefully estimated their differences of zero-points. In total 10214 radial velocity measurements have been obtained for the 1420 stars. With a mean time baseline of 6.35 years, 92.9% of the candidates fulfill a target stability criterion of 300 m/s. 343 stars are found to be constant at the level of 100 m/s over 10 years. Comparisons with the catalogues of Nidever et al.(2002) and Chubak et al.(2012) show an excellent agreement for FGK stars, with zero-point differences lower than 100 m/s and a remarkably low RMS scatter of 33 m/s in the case of Nidever et al.(2002), suggesting that the precision of the catalogue presented here is better than this value. This catalog will likely be useful for other large-scale spectroscopic surveys such as APOGEE, Gaia-ESO, HERMES and LAMOST.

Context. As increasingly more spectroscopic data are being delivered by medium- and high-resolving power multi-object spectrographs, more automatic stellar parameter determination softwares are being developed. The quality of the spectra collected also allows the determination of elemental abundances. Aims. SPADES is an automated software for determining: the radial velocity (Vr), the effective temperature (Teff), the surface gravity (log g), the metallicity ([Fe/H]), and most importantly, the individual abundances. In this first version it is targeted on the analysis of mid-F-G dwarfs, but is meant to evolve to analyze any type of single stars. Methods. SPADES relies on a line-by-line modeling to determine the stellar parameters. Results. The internal systematic and random errors of SPADES were assessed by Monte Carlo method simulations with synthetic spectra and the external systematic errors by analysing real ground-based observed spectra. For example, by simulating the Giraffe setups HR13 and HR14B with synthetic spectra for a dwarf with Teff = 5800 K, log g = 4.5, [Fe/H] = 0.0 dex and with a signal-tonoise ratio (S/N) of 100, the stellar parameters are recovered with no significant bias and with 1-{\sigma} precisions of 8 K for Teff, 0.05 for log g, 0.009 for [Fe/H], 0.003 for [Ti/Fe] and 0.01 for [Ni/Fe].

With the large amounts of spectroscopic data available today and the very large surveys to come (e.g. Gaia), the need for automatic data analysis software is unquestionable. We thus developed an automatic spectra analysis program for the determination of stellar parameters: radial velocity, effective temperature, surface gravity, micro-turbulence, metallicity and the elemental abundances of the elements present in the spectral range. Target stars for this software should include all types of stars. The analysis method relies on a line by line comparison of the spectrum of a target star to a library of synthetic spectra. The idea is built on the experience acquired in developing the TGMET (Katz et al. 1998 and Soubiran et al. 2003) ETOILE (Katz 2001) and Abbo (Bonifacio & Caffau 2003) softwares. The method is presented and the performances are illustrated with GIRAFFE-like simulated spectra with high resolution (R = 25000), with high and low signal to noise ratios (down to SNR= 30). These spectra should be close to what could be targeted by the Gaia-ESO Survey (GCDS).

The Radial Velocity Spectrometer (RVS) on board of Gaia will perform a large spectroscopic survey to determine the radial velocities of some 1.5x10^8 stars. We present the status of ground-based observations of a sample of 1420 candidate standard stars designed to calibrate the RVS. Each candidate star has to be observed several times before Gaia launch (and at least once during the mission) to ensure that its radial velocity remains stable during the whole mission. Observations are performed with the high-resolution spectrographs SOPHIE, NARVAL and CORALIE, completed with archival data of the ELODIE and HARPS instruments. The analysis shows that about 7% of the current catalogue exhibits variations larger than the adopted threshold of 300 m/s. Consequently, those stars should be rejected as reference targets, due to the expected accuracy of the Gaia RVS. Emphasis is also put here on our observations of bright asteroids to calibrate the ground-based velocities by a direct comparison with celestial mechanics. It is shown that the radial velocity zero points of SOPHIE, NARVAL and CORALIE are consistent with each other, within the uncertainties. Despite some scatter, their temporal variations remain small with respect to our adopted stability criterion.

This chapter describes how astronomical imaging survey data have become a vital part of modern astronomy, how these data are archived and then served to the astronomical community through on-line data access portals. The Virtual Observatory, now under development, aims to make all these data accessible through a uniform set of interfaces. This chapter also describes the scientific need for one common image processing task, that of composing individual images into large scale mosaics and introduces Montage as a tool for this task. Montage, as distributed, can be used in four ways: as a single thread/process on a single CPU, in parallel using MPI to distribute similar tasks across a parallel computer, in parallel using grid tools (Pegasus/DAGMan) to distributed tasks across a grid, or in parallel using a script-driven approach (Swift). An on-request web based Montage service is available for users who do not need to build a local version. We also introduce some work on a new scripted version of Montage, which offers ease of customization for users. Then, we discuss various ideas where Web 2.0 technologies can help the Montage community.

The calibration of the Radial Velocity Spectrometer (RVS) onboard the ESA Gaia satellite (to be launched in 2012) requires a list of standard stars with a radial velocity (RV) known with an accuracy of at least 300 m/s. The IAU Commission 30 lists of RV standard stars are too bright and not dense enough. We describe the selection criteria due to the RVS constraints for building an adequate full-sky list of at least 1000 RV standards from catalogues already published in the literature. A preliminary list of 1420 candidate standard stars is built and its properties are shown. An important re-observation programme has been set up in order to ensure within it the selection of objects with a good stability until the end of the Gaia mission (around 2018). The present list of candidate standards is available at CDS and usable for many other projects.

A new full-sky catalog of Radial Velocity standard stars is being built for the determination of the Radial Velocity Zero Point of the RVS on board of Gaia. After a careful selection of 1420 candidates matching well defined criteria, we are now observing all of them to verify that they are stable enough over several years to be qualified as reference stars. We present the status of this long-term observing programme on three spectrographs : SOPHIE, NARVAL and CORALIE, complemented by the ELODIE and HARPS archives. Because each instrument has its own zero-point, we observe intensively IAU RV standards and asteroids to homogenize the radial velocity measurements. We can already estimate that ~8% of the candidates have to be rejected because of variations larger than the requested level of 300 m/s.

The lack of radial velocity data in the Hipparcos catalogue was considered a significant deficiency, so when Gaia was conceived, a spectrometer was a core constituent of its payload. The Gaia Radial Velocity Spectrometer faced a number of design challenges, in particular set by the need to balance kinematic and astrophysical capability. We present an overview of the evolution of the instrument to its present form, identifying the competing technical, performance and programmatic factors which have shaped it.

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.

Montage is a portable software toolkit for constructing custom, science-grade mosaics by composing multiple astronomical images. The mosaics constructed by Montage preserve the astrometry (position) and photometry (intensity) of the sources in the input images. The mosaic to be constructed is specified by the user in terms of a set of parameters, including dataset and wavelength to be used, location and size on the sky, coordinate system and projection, and spatial sampling rate. Many astronomical datasets are massive, and are stored in distributed archives that are, in most cases, remote with respect to the available computational resources. Montage can be run on both single- and multi-processor computers, including clusters and grids. Standard grid tools are used to run Montage in the case where the data or computers used to construct a mosaic are located remotely on the Internet. This paper describes the architecture, algorithms, and usage of Montage as both a software toolkit and as a grid portal. Timing results are provided to show how Montage performance scales with number of processors on a cluster computer. In addition, we compare the performance of two methods of running Montage in parallel on a grid.

This document is divided in two parts. The first part deals with the radial velocities (RV) distributions for B-type stars and nebulosities observed with the VLT-GIRAFFE in the Large and Small Magellanic Clouds towards the open clusters NGC2004 and NGC330. Thanks to the resolution of GIRAFFE spectra, we found that the RV distribution for the nebulosities in the LMC is bi-modal. This bi-modality can be interpreted, in term of dynamics, by the expansion of the LMC4 superbubble. The second part deals with the GAIA space mission and the determination of the radial velocities by using Radial Velocity Spectrometer (RVS) spectra. The methods to determine the radial velocities are presented as well as preliminary results on simulated RVS spectra.

Gaia is an ambitious space astrometry mission of ESA with a main objective to map the sky in astrometry and photometry down to a magnitude 20 by the end of the next decade. While the mission is built and operated by ESA and an industrial consortium, the data processing is entrusted to a consortium formed by the scientific community, which was formed in 2006 and formally selected by ESA one year later. The satellite will downlink around 100 TB of raw telemetry data over a mission duration of 5 years from which a very complex iterative processing will lead to the final science output: astrometry with a final accuracy of a few tens of microarcseconds, epoch photometry in wide and narrow bands, radial velocity and spectra for the stars brighter than 17 mag. We discuss the general principles and main difficulties of this very large data processing and present the organisation of the European Consortium responsible for its design and implementation.

The Naval Research Laboratory and the National Radio Astronomy Observatory completed implementation of a low frequency capability on the VLA at 73.8 MHz in 1998. This frequency band offers unprecedented sensitivity (~25 mJy/beam) and resolution (~25 arcsec) for low-frequency observations. We review the hardware, the calibration and imaging strategies, comparing them to those at higher frequencies, including aspects of interference excision and wide-field imaging. Ionospheric phase fluctuations pose the major difficulty in calibrating the array. Over restricted fields of view or at times of extremely quiescent ionospheric ``weather'', an angle-invariant calibration strategy can be used. In this approach a single phase correction is devised for each antenna, typically via self-calibration. Over larger fields of view or at times of more normal ionospheric ``weather'' when the ionospheric isoplanatic patch size is smaller than the field of view, we adopt a field-based strategy in which the phase correction depends upon location within the field of view. This second calibration strategy was implemented by modeling the ionosphere above the array using Zernike polynomials. Images of 3C sources of moderate strength are provided as examples of routine, angle-invariant calibration and imaging. Flux density measurements indicate that the 74 MHz flux scale at the VLA is stable to a few percent, and tied to the Baars et al. value of Cygnus A at the 5 percent level. We also present an example of a wide-field image, devoid of bright objects and containing hundreds of weaker sources, constructed from the field-based calibration. We close with a summary of lessons the 74 MHz system offers as a model for new and developing low-frequency telescopes. (Abridged)

Gaia is ESA's ambitious space astrometry mission the main objective of which is to astrometrically and spectro-photometrically map 1000 Million celestial objects (mostly in our galaxy) with unprecedented accuracy. The announcement of opportunity for the data processing will be issued by ESA late in 2006. The Gaia Data Processing and Analysis Consortium (DPAC) has been formed recently and is preparing an answer. The satellite will downlink close to 100 TB of raw telemetry data over 5 years. To achieve its required accuracy of a few 10s of Microarcsecond astrometry, a highly involved processing of this data is required. In addition to the main astrometric instrument Gaia will host a Radial Velocity instrument, two low-resolution dispersers for multi-color photometry and two Star Mappers. Gaia is a flying Giga Pixel camera. The various instruments each require relatively complex processing while at the same time being interdependent. We describe the overall composition of the DPAC and the envisaged overall architecture of the Gaia data processing system. We shall delve further into the core processing - one of the nine, so-called, coordination units comprising the Gaia processing system.

The Gaia mission is designed as a Galaxy explorer, and will measure simultaneously, in a survey mode, the five or six phase space parameters of all stars brighter than 20th magnitude, as well as providing a description of their astrophysical characteristics. These measurements are obtained by combining an astrometric instrument with micro-arcsecond capabilities, a photometric system giving the magnitudes and colours in 15 bands and a medium resolution spectrograph named the Radial Velocity Spectrometer (RVS). The latter instrument will produce spectra in the 848 to 874 nm wavelength range, with a resolving power R = 11 500, from which radial velocities, rotational velocities, atmospheric parameters and abundances can be derived. A companion paper (Katz et al. 2004) presents the characteristics of the RVS and its performance. This paper details the outstanding scientific impact of this important part of the Gaia satellite on some key open questions in present day astrophysics. The unbiased and simultaneous acquisition of multi-epoch radial velocities and individual abundances of key elements in parallel with the astrometric parameters is essential for the determination of the dynamical state and formation history of our Galaxy. Moreover, for stars brighter than V=15, the resolving power of the RVS will give information about most of the effects which influence the position of a star in the Hertzsprung-Russell diagram, placing unprecedented constraints on the age, internal structure and evolution of stars of all types. Finally, the RVS multi-epoch observations are ideally suited to the identification, classification and characterisation of the many types of double, multiple and variable stars.

Matched resolution multi-frequency VLA observations of four radio galaxies are used to derive the asymptotic low energy slope of the relativistic electron distribution. Where available, low energy slopes are also determined for other sources in the literature. They provide information on the acceleration physics independent of radiative and other losses, which confuse measurements of the synchrotron spectra in most radio, optical and X-ray studies. We find a narrow range of inferred low energy electron energy slopes, n(E)=const*E^-2.1 for the currently small sample of lower luminosity sources classified as FRI (not classical doubles). This distribution is close to, but apparently inconsistent with, the test particle limit of n(E)=const*E^-2.0 expected from strong diffusive shock acceleration in the non-relativistic limit. Relativistic shocks or those modified by the back-pressure of efficiently accelerated cosmic rays are two alternatives to produce somewhat steeper spectra. We note for further study the possiblity of acceleration through shocks, turbulence or shear in the flaring/brightening regions in FRI jets as they move away from the nucleus. Jets on pc scales and the collimated jets and hot spots of FRII (classical double) sources would be governed by different acceleration sites and mechanisms; they appear to show a much wider range of spectra than for FRI sources.

The Grist project (http://grist.caltech.edu/) is developing a grid-technology based system as a research environment for astronomy with massive and complex datasets. This knowledge extraction system will consist of a library of distributed grid services controlled by a workflow system, compliant with standards emerging from the grid computing, web services, and virtual observatory communities. This new technology is being used to find high redshift quasars, study peculiar variable objects, search for transients in real time, and fit SDSS QSO spectra to measure black hole masses. Grist services are also a component of the ``hyperatlas'' project to serve high-resolution multi-wavelength imagery over the Internet. In support of these science and outreach objectives, the Grist framework will provide the enabling fabric to tie together distributed grid services in the areas of data access, federation, mining, subsetting, source extraction, image mosaicking, statistics, and visualization.

Stellar spectra obtained with Radial Velocity Spectrometer on-board Gaia will enable determination of projected rotational velocities (v sin i) from rotational broadening of spectral lines. To estimate the accuracy with which v sin i can be determined for different stellar types and as a function of magnitude, we perform simulations in which "observed" spectra (as obtained from Gaia RVS simulator) are fitted with Kurucz spectra from template library using the least square method. We compare results on the v sin i accuracy as obtained in the case of spectra differing only in v sin i and in the case of a more general spectral mismatch (in T_eff, log g, [Fe/H]).

The GAUDI database (Ground-based Asteroseismology Uniform Database Interface, http://sdc.laeff.esa.es/gaudi/) is a preparatory archive for the COROT (COnvection, ROtation and planetary Transits, http://www.astrsp-mrs.fr/projets/corot/) mission developed at LAEFF (Laboratory for Space Astrophysics and Theoretical Physics, http://www.laeff.esa.es). Its intention is to make the ground-based observations obtained in the preparation of the asteroseismology programme available in a simple and efficient way. It contains spectroscopic and photometric data together with inferred physical parameters for more than 1500 objects gathered since January 1998 in 6 years of observational campaigns. In this paper, the main functionalities and characteristics of the system are described. The observations have been collected at ESO-La Silla, Telescopio Nazionale Galileo, Observatoire de Haute-Provence, South African Astronomical Observatory, Tautenberg Observatory and Sierra Nevada Observatory.

The definition and optimisation studies for the Gaia satellite spectrograph, the Radial Velocity Spectrometer (RVS), converged in late 2002 with the adoption of the instrument baseline. This paper reviews the characteristics of the selected configuration and presents its expected performance. The RVS is a 2.0 by 1.6 degree integral field spectrograph, dispersing the light of all sources entering its field of view with a resolving power R=11 500 over the wavelength range [848, 874] nm. The RVS will continuously and repeatedly scan the sky during the 5 years of the Gaia mission. On average, each source will be observed 102 times over this period. The RVS will collect the spectra of about 100-150 million stars up to magnitude V~17-18. At the end of the mission, the RVS will provide radial velocities with precisions of ~2 km/s at V=15 and \~15-20 km/s at V=17, for a solar metallicity G5 dwarf. The RVS will also provide rotational velocities, with precisions (at the end of the mission) for late type stars of sigma_vsini ~5 km/s at V~15 as well as atmospheric parameters up to V~14-15. The individual abundances of elements such as Silicon and Magnesium, vital for the understanding of Galactic evolution, will be obtained up to V~12-13. Finally, the presence of the 862.0 nm Diffuse Interstellar Band (DIB) in the RVS wavelength range will make it possible to derive the three dimensional structure of the interstellar reddening.

One of the goals of the ground-based support program for the COROT and Roemer satellite missions is to characterize suitable target stars for the part of the missions dedicated to asteroseismology. We present the detailed abundance analysis of nine of the potential COROT main targets using the semi-automatic software VWA. For two additional COROT targets we could not perform the analysis due to the high rotational velocity of these stars. For five stars with low rotational velocity we have also performed abundance analysis by a classical equivalent width method in order to test the reliability of the VWA software. The agreement between the different methods is good. We find that it is necessary to measure abundances extracted from each line relative to the abundances found from a spectrum of the Sun in order to remove systematic errors. We have constrained the global atmospheric parameters Teff, log g, and [Fe/H] to within 70-100 K, 0.1-0.2 dex, and 0.1 dex for five stars which are slow rotators (v sin i < 15 km/s). For most of the stars stars we find good agreement with the parameters found from line depth ratios, H-alpha lines, Str\"omgren indices, previous spectroscopic studies, and also log g determined from the HIPPARCOS parallaxes. For the fast rotators v sin i > 60 km/s it is not possible to constrain the atmospheric parameters.

We report the determination from high-resolution spectra of the atmospheric parameters and abundances of 13 chemical species (among which lithium) in 8 single-lined active binaries. These data are combined with our previous results for 6 other RS CVn systems to examine a possible relationship between the photospheric abundance patterns and the stellar activity level. The stars analyzed are generally found to exhibit peculiar abundance ratios compared to inactive, galactic disk stars of similar metallicities. We argue that this behaviour is unlikely an artefact of errors in the determination of the atmospheric parameters or non-standard mixing processes along the red giant branch, but diagnoses instead the combined action of various physical processes related to activity. The most promising candidates are cool spot groups covering a very substantial fraction of the stellar photosphere or NLTE effects arising from nonthermal excitation. However, we cannot exclude the possibility that more general shortcomings in our understanding of K-type stars (e.g. inadequacies in the atmospheric models) also play a significant role. Lastly, we call attention to the unreliability of the (V-R) and (V-I) colour indices as temperature indicators in chromospherically active stars.

Photospheric parameters and abundances are presented for a sample of single-lined chromospherically active binaries from a differential LTE analysis of high-resolution spectra. Abundances have been derived for 13 chemical species, including several key elements such as Li, Mg, and Ca. Two methods have been used. The effective temperatures, surface gravities and microturbulent velocities were first derived from a fully self-consistent analysis of the spectra, whereby the temperature is determined from the excitation equilibrium of the Fe I lines. The second approach relies on temperatures derived from the (B-V) colour index. These two methods give broadly consistent results for the stars in our sample, suggesting that the neutral iron lines are formed under conditions close to LTE. We discuss the reliability in the context of chromospherically active stars of various colour indices used as temperature indicators, and conclude that the (V-R) and (V-I) colours are likely to be significantly affected by activity processes. Irrespective of the method used, our results indicate that the X-ray active binaries studied are not as metal poor as previously claimed, but are at most mildly iron-depleted relative to the Sun (-0.41<[Fe/H]<+0.11). A significant overabundance of several chemical species is observed (e.g., the alpha-synthezised elements). These abundance patterns are discussed in relation to stellar activity.

The high-resolution optical spectra of the two X-ray active binaries RS CVn stars HD 113816 (IS Vir) and HD 119285 (V851 Cen) are analysed and their Na, Mg, Al, Si, Ca, Sc, Ti, Co and Ni contents determined. The analysis of IS Vir and V851 Cen is performed with three different LTE methods. In the first one, abundances are derived using measured equivalent widths and Kurucz LTE model atmospheres as input for the MOOG software package. The input atmospheric parameters and abundances are iteratively modified until (i) the Fe I abundances exhibit no trend with excitation potential or equivalent width, (ii) Fe I and Fe II average abundances are the same and (iii) Fe and Alpha elements average abundances are consistent with the input values. The second method follows a similar approach, but relies on the B-V and V-I colour indices to determine the temperature. The third method relies on fitting the 6162 A Ca I line wing profiles to derive the surface gravity. The reliability of these methods is investigated in the context of single line RS CVn stars. It is shown that the V-I photometric index gives, on a broader sample of stars, significantly cooler estimates of the effective temperature than the B-V index. All other approaches give results in good agreement with each other. The analysis of IS Vir and V851 Cen results in both cases in their primaries being giant stars of near-solar metallicity. Their parameters as derived with the first method are respectively Teff = 4720 K, log g = 2.65, [Fe/H] = +0.04 and Teff = 4700 K, log g = 3.0 and [Fe/H] = -0.13. In the case of V851 Cen the derived iron content is significantly higher than a previous determination in the literature. Both stars exhibit relative overabundances of several elements (e.g. Ca) with respect to the solar mix.

We present several new sets of grids of model stellar atmospheres computed with modified versions of the ATLAS9 code. Each individual set consists of several grids of models with different metallicities ranging from [M/H] = -2.0 to +1.0 dex. The grids range from 4000 to 10000 K in T_eff and from 2.0 to 5.0 dex in logg. The individual sets differ from each other and from previous ones essentially in the physics used for the treatment of the convective energy transport, in the higher vertical resolution of the atmospheres and in a finer grid in the (T_eff, logg) plane. These improvements enable the computation of derivatives of color indices accurate enough for pulsation mode identification. In addition, we show that the chosen vertical resolution is necessary and sufficient for the purpose of stellar interior modelling. To explain the physical differences between the model grids we provide a description of the currently available modifications of ATLAS9 according to their treatment of convection. Our critical analysis of the dependence of the atmospheric structure and observable quantities on convection treatment, vertical resolution and metallicity reveals that spectroscopic and photometric observations are best represented when using an inefficient convection treatment. This conclusion holds whatever convection formulation investigated here is used, i.e. MLT(alpha=0.5), CM and CGM are equivalent. We also find that changing the convection treatment can lead to a change in the effective temperature estimated from Stroemgren color indices from 200 to 400 K.

We present a method estimating the atmospheric parameters Teff, log g, [Fe/H] for stars observed, even at low signal to noise ratio, with the echelle spectrograph ELODIE on the 193cm telescope at Observatoire de Haute-Provence. The method relies on the least-square comparison of the spectrum of a target star to a library of 211 spectra of reference stars for which the atmospheric parameters are well known and which were observed with the same instrument. In order to obtain a meaningful comparison between target and reference spectra, all features which are not intrinsic to the objects must be removed. This treatment involves the correction of the blaze efficiency for each order, cosmic rays hits and telluric line removal, convolution of the spectra to a common spectral resolution, wavelength scale and flux level adjustment. The library available at the present time covers the effective temperature range [4000K, 6300K], the metallicity range [-2.9, +0.35] and the gravities of both unevolved and evolved stars existing at these temperatures and metallicities. Tests performed with the actual library allow us to estimate the internal accuracy to be 86 K, 0.28 dex and 0.16 dex for Teff, log g, [Fe/H] for a target star with S/N = 100 and 102 K, 0.29 dex and 0.17 dex at S/N = 10. This accuracy will improve in the future as the number of reference stars in the library will increase. The software (named TGMET) has been installed at Observatoire de Haute-Provence for the on-line analysis of the high-resolution spectra of ELODIE, which was originally conceived for accurate radial velocity measurements.

A library of 211 echelle spectra taken with ELODIE at the Observatoire de Haute-Provence is presented. It provides a set of spectroscopic standards covering the full range of gravities and metallicities in the effective temperature interval [4000 K, 6300 K]. The spectra are straightened, wavelength calibrated, cleaned of cosmic ray hits, bad pixels and telluric lines. They cover the spectral range [440 nm, 680 nm] with an instrumental resolution of 42000. For each star, basic data were compiled from the Hipparcos catalogue and the Hipparcos Input Catalogue. Radial velocities with a precision better than 100 m/s are given. Atmospheric parameters (Teff, log g, [Fe/H]) from the literature are discussed. Because of scattered determinations in the bibliography, even for the most well-known stars, these parameters were adjusted by an iterative process which takes account of common or different spectral features between the standards, using our homogeneous set of spectra. Revised values of (Teff, log g, [Fe/H]) are proposed. They are still consistent with the literature, and also lead to the self-consistency of the library, in the sense that similar spectra have similar atmospheric parameters. This adjustment was performed by using step by step a method based on the least square comparison of carefully prepared spectra, which was originally developed for the on-line estimation of the atmospheric parameters of faint field stars (companion paper in the main journal). The spectra and corresponding data will only be available in electronic form at the CDS (ftp cdsarc.u-strasbg.fr or http://cdsweb.u-strasbg.fr/Abstract.html).

We show how previous works (Fuhrmann et al., van 't Veer-Menneret & Megessier) demonstrate the efficiency of the use of Balmer line profiles for effective temperature determination. In agreement with them, we insist on the physical interest of this method based on the behaviour of these lines with the variations of the parameters involved in the treatment of the convective transport. The comparison between Fuhrmann's results and ours, independently obtained, exhibits a quite good agreement. We show new results of effective temperature, gravity and metallicities for a few of our programme stars, ranging from solar to overabundant metallicities.

The extragalactic applications of NIR surveys are summarized with a focus on the ability to map the interstellar extinction of our Galaxy. Very preliminary extraction of galaxies on a set of 180 consecutive images is presented, and the results illustrate some of the pitfalls in attempting an homogeneous extraction of galaxies from these wide-angle and shallow surveys.