The Gaia-ESO Survey: radial metallicity gradients and age-metallicity relation of stars in the Milky Way disk

We study the relationship between age, metallicity, and α-enhancement of FGK stars in the Galactic disk. The results are based upon the analysis of high-resolution UVES spectra from the Gaia-ESO large stellar survey. We explore the limitations of the observed dataset, i.e. the accuracy of stellar parameters and the selection effects that are caused by the photometric target preselection. We find that the colour and magnitude cuts in the survey suppress old metal-rich stars and young metal-poor stars. This suppression may be as high as 97% in some regions of the age-metallicity relationship. The dataset consists of 144 stars with a wide range of ages from 0.5 Gyr to 13.5 Gyr, Galactocentric distances from 6 kpcto 9.5 kpc, and vertical distances from the plane 0 < |Z| < 1.5 kpc. On this basis, we find that i) the observed age-metallicity relation is nearly flat in the range of ages between 0 Gyr and 8 Gyr; ii) at ages older than 9 Gyr, we see a decrease in [Fe/H] and a clear absence of metal-rich stars; this cannot be explained by the survey selection functions; iii) there is a significant scatter of [Fe/H] at any age; and iv) [Mg/Fe] increases with age, but the dispersion of [Mg/Fe] at ages >9 Gyr is not as small as advocated by some other studies. In agreement with earlier work, we find that radial abundance gradients change as a function of vertical distance from the plane. The [Mg/Fe] gradient steepens and becomes negative. In addition, we show that the inner disk is not only more α-rich compared to the outer disk, but also older, as traced independently by the ages and Mg abundances of stars.

PLANETARY NEBULAE IN THE INNER MILKY WAY II: THE BULGE-DISK TRANSITION

In this work, a sample of planetary nebulae located in the inner-disk and bulge of the Galaxy is used in order to find the galactocentric distance which better separates these two populations, from the point of view of abundances. Statistical distance scales are used to study the distribution of abundances across the disk-bulge interface. A Kolmogorov-Smirnov test is used to find the distance at which the chemical properties of these regions better separate. The results of the statistical analysis indicate that, on the average, the inner population has lower abundances than the outer. Additionally, for the a-element abundances, the inner population does not follow the disk radial gradient towards the galactic center. Based on our results, we suggest a bulge-disk interface at 1.5 kpc, marking the transition between the bulge and inner-disk of the Galaxy as defined by the intermediate mass population.

The Milky Way's Dark Matter Halo Appears To Be Lopsided

The atomic hydrogen gas (H I) disk in the outer region (beyond similar to 10 kpc from the center) of Milky Way can provide valuable information about the structure of the dark matter halo. The recent three-dimensional thickness map of the outer H I disk from the all sky 21 cm line Leiden/Argentine/Bonn survey, gives us a unique opportunity to investigate the structure of the dark matter halo of Milky Way in great detail. A striking feature of this new survey is the north-south (N-S) asymmetry in the thickness map of the atomic hydrogen gas. Assuming vertical hydrostatic equilibrium under the total potential of the Galaxy, we derive the model thickness map of the H I gas. We show that simple axisymmetric halo models, such as softened isothermal halo (producing a flat rotation curve with V-c similar to 220 km s(-1)) or any halo with density falling faster than the isothermal one, are not able to explain the observed radial variation of the gas thickness. We also show that such axisymmetric halos along with different H I velocity dispersion in the two halves, cannot explain the observed asymmetry in the thickness map. Amongst the nonaxisymmetric models, it is shown that a purely lopsided (m = 1, first harmonic) dark matter halo with reasonable H I velocity dispersion fails to explain the N-S asymmetry satisfactorily. However, we show that by superposing a second harmonic (m = 2) out of phase onto a purely lopsided halo, e. g., our best fit and more acceptable model A (with parameters epsilon(1)(h) = 0.2, epsilon(2)(h) = 0.18, and sigma(H I) = 8.5 km s(-1)) can provide an excellent fit to the observation and reproduce the N-S asymmetry naturally. The emerging picture of the asymmetric dark matter halo is supported by the. cold dark matter halos formed in the cosmological N-body simulation.

VISTA Variables in the Via Lactea (VVV): The public ESO near-IR variability survey of the Milky Way

We describe the public ESO near-IR variability survey (VVV) scanning the Milky Way bulge and an adjacent section of the mid-plane where star formation activity is high. The survey will take 1929 h of observations with the 4-m VISTA telescope during 5 years (2010-2014), covering similar to 10(9) point sources across an area of 520 deg(2), including 33 known globular clusters and similar to 350 open clusters. The final product will be a deep near-IR atlas in five passbands (0.9-2.5 mu m) and a catalogue of more than 106 variable point sources. Unlike single-epoch surveys that, in most cases, only produce 2-D maps, the VVV variable star survey will enable the construction of a 3-D map of the surveyed region using well-understood distance indicators such as RR Lyrae stars, and Cepheids. It will yield important information on the ages of the populations. The observations will be combined with data from MACHO, OGLE, EROS, VST, Spitzer, HST, Chandra, INTEGRAL, WISE, Fermi LAT, XMM-Newton, GAIA and ALMA for a complete understanding of the variable sources in the inner Milky Way. This public survey will provide data available to the whole community and therefore will enable further studies of the history of the Milky Way, its globular cluster evolution, and the population census of the Galactic Bulge and center, as well as the investigations of the star forming regions in the disk. The combined variable star catalogues will have important implications for theoretical investigations of pulsation properties of stars. (C) 2009 Elsevier B.V. All rights reserved.

PLANETARY NEBULAE IN THE INNER MILKY WAY: NEW ABUNDANCES

The study of planetary nebulae in the inner-disk and bulge gives important information on the chemical abundances of elements such as He, N, O, Ar, Ne, and on the evolution of these abundances, which is associated with the evolution of intermediate-mass stars and the chemical evolution of time Galaxy. We present accurate abundances of the elements He, N, 5, 0, Ar, and Ne for a sample of 54 planetary nebulae located towards the bulge of the Galaxy, for 33 of which the abundances are derived here for the first time. The abundances are obtained based on observations in the optical domain made at the National Laboratory for Astrophysics (LNA, Brazil). The data show a good agreement; with other results in the literature, in the sense that the distribution of the abundances is similar to that of those works.

VVV DR1: The first data release of the Milky Way bulge and southern plane from the near-infrared ESO public survey VISTA variables in the Via Lactea

Context. The ESO public survey VISTA variables in the Via Lactea (VVV) started in 2010. VVV targets 562 sq. deg in the Galactic bulge and an adjacent plane region and is expected to run for about five years. Aims. We describe the progress of the survey observations in the first observing season, the observing strategy, and quality of the data obtained. Methods. The observations are carried out on the 4-m VISTA telescope in the ZYJHK(s) filters. In addition to the multi-band imaging the variability monitoring campaign in the K-s filter has started. Data reduction is carried out using the pipeline at the Cambridge Astronomical Survey Unit. The photometric and astrometric calibration is performed via the numerous 2MASS sources observed in each pointing. Results. The first data release contains the aperture photometry and astrometric catalogues for 348 individual pointings in the ZYJHK(s) filters taken in the 2010 observing season. The typical image quality is similar to 0 ''.9-1 ''.0. The stringent photometric and image quality requirements of the survey are satisfied in 100% of the JHK(s) images in the disk area and 90% of the JHK(s) images in the bulge area. The completeness in the Z and Y images is 84% in the disk, and 40% in the bulge. The first season catalogues contain 1.28 x 10(8) stellar sources in the bulge and 1.68 x 10(8) in the disk area detected in at least one of the photometric bands. The combined, multi-band catalogues contain more than 1.63 x 10(8) stellar sources. About 10% of these are double detections because of overlapping adjacent pointings. These overlapping multiple detections are used to characterise the quality of the data. The images in the JHK(s) bands extend typically similar to 4 mag deeper than 2MASS. The magnitude limit and photometric quality depend strongly on crowding in the inner Galactic regions. The astrometry for K-s = 15-18 mag has rms similar to 35-175 mas. Conclusions. The VVV Survey data products offer a unique dataset to map the stellar populations in the Galactic bulge and the adjacent plane and provide an exciting new tool for the study of the structure, content, and star-formation history of our Galaxy, as well as for investigations of the newly discovered star clusters, star-forming regions in the disk, high proper motion stars, asteroids, planetary nebulae, and other interesting objects.

The star formation rate in the inner Milky Way Galaxy

The present star formation rate (SFR) in the inner Galaxy is puzzling for the chemical evolution models (CEM). No static CEM is able to reproduce the peak of the SFR in the 4 kpc ring. The main reason is probably a shortage of gas, which could be due to the dynamical effects produced by the galactic bar, not considered by these models. We developed a CEM that includes radial gas flows in order to mimic the effects of the galactic bar in the first 5 kpc of the galactic disk. In this model, the star formation (SF) is a two-step process: first, the diffuse gas forms molecular clouds. Then, stars form from cloud-cloud collisions or by the interaction between massive stars and the molecular gas. The former is called spontaneous and the latter induced SF. The mass in the different phases of each region changes by the processes associated with the stellar formation and death by: the SF due to spontaneous fragmentation of gas in the halo; formation of gas clouds in the disk from the diffuse gas; induced SF in the disk due to the interaction between massive stars and gas clouds; and finally, the restitution of the diffuse gas associated to these process of cloud and star formation. In the halo, the star formation rate for the diffuse gas follows a Schmidt law with a power n = 1.5. In the disk, the stars form in two steps: first, molecular clouds are formed from the diffuse gas also following a Schmidt law with n=1.5 and a proportionality factor. Including a specific pattern of radial gas flows, the CEM is able to reproduce with success the peak in the SFR at 4 kpc (fig. 1).

Circumstellar emission in Be/X-ray binaries of the Magellanic Clouds and the Milky Way

Aims. We study the optical and near-infrared colour excesses produced by circumstellar emission in a sample of Be/X-ray binaries. Our main goals are exploring whether previously published relations, valid for isolated Be stars, are applicable to Be/X-ray binaries and computing the distance to these systems after correcting for the effects of the circumstellar contamination. Methods. Simultaneous UBVRI photometry and spectra in the 3500−7000 Å spectral range were obtained for 11 optical counterparts to Be/X-ray binaries in the LMC, 5 in the SMC and 12 in the Milky Way. As a measure of the amount of circumstellar emission we used the Hα equivalent width corrected for photospheric absorption. Results. We find a linear relationship between the strength of the Hα emission line and the component of E(B − V) originating from the circumstellar disk. This relationship is valid for stars with emission lines weaker than EW ≈ −15   Å. Beyond this point, the circumstellar contribution to E(B − V) saturates at a value ≈0.17   mag. A similar relationship is found for the (V − I) near infrared colour excess, albeit with a steeper slope and saturation level. The circumstellar excess in (B − V) is found to be about five times higher for Be/X-ray binaries than for isolated Be stars with the same equivalent width EW(Hα), implying significant differences in the physical properties of their circumstellar envelopes. The distance to Be/X-ray binaries (with non-shell Be star companions) can only be correctly estimated by taking into account the excess emission in the V band produced by free-free and free-bound transitions in the circumstellar envelope. We provide a simple method to determine the distances that includes this effect.

Minor mergers and their impact on the kinematics of old and young stellar populations in disk galaxies

By means of N-body simulations we investigate the impact of minor mergers on the angular momentum and dynamical properties of the merger remnant. Our simulations cover a range of initial orbital characteristics and gas-to-stellar mass fractions (from 0 to 20%), and include star formation and supernova feedback. We confirm and extend previous results by showing that the specific angular momentum of the stellar component always decreases independently of the orbital parameters or morphology of the satellite, and that the decrease in the rotation velocity of the primary galaxy is accompanied by a change in the anisotropy of the orbits. However, the decrease affects only the old stellar population, and not the new population formed from gas during the merging process. This means that the merging process induces an increasing difference in the rotational support of the old and young stellar components, with the old one lagging with respect to the new. Even if our models are not intended specifically to reproduce the Milky Way and its accretion history, we find that, under certain conditions, the modeled rotational lag found is compatible with that observed in the Milky Way disk, thus indicating that minor mergers can be a viable way to produce it. The lag can increase with the vertical distance from the disk midplane, but only if the satellite is accreted along a direct orbit, and in all cases the main contribution to the lag comes from stars originally in the primary disk rather than from stars in the satellite galaxy. We also discuss the possibility of creating counter-rotating stars in the remnant disk, their fraction as a function of the vertical distance from the galaxy midplane, and the cumulative effect of multiple mergers on their creation.

Evidence for dark matter in the inner Milky Way

The ubiquitous presence of dark matter in the Universe is today a central tenet in modern cosmology and astrophysics(1). Throughout the Universe, the evidence for dark matter is compelling in dwarfs, spiral galaxies, galaxy clusters as well as at cosmological scales. However, it has been historically difficult to pin down the dark matter contribution to the total mass density in the Milky Way, particularly in the innermost regions of the Galaxy and in the solar neighbourhood(2). Here we present an up-to-date compilation of Milky Way rotation curve measurements(3-13), and compare it with state-of-the-art baryonic mass distribution models(14-26). We show that current data strongly disfavour baryons as the sole contribution to the Galactic mass budget, even inside the solar circle. Our findings demonstrate the existence of dark matter in the inner Galaxy without making any assumptions about its distribution. We anticipate that this result will compel new model-independent constraints on the dark matter local density and profile, thus reducing uncertainties on direct and indirect dark matter searches, and will help reveal the structure and evolution of the Galaxy.

The dark matter profile of the Milky Way: a non-parametric reconstruction

We present the results of a new, non-parametric method to reconstruct the Galactic dark matter profile directly from observations. Using the latest kinematic data to track the total gravitational potential and the observed distribution of stars and gas to set the baryonic component, we infer the dark matter contribution to the circular velocity across the Galaxy. The radial derivative of this dynamical contribution is then estimated to extract the dark matter profile. The innovative feature of our approach is that it makes no assumption on the functional form or shape of the profile, thus allowing for a clean determination with no theoretical bias. We illustrate the power of the method by constraining the spherical dark matter profile between 2.5 and 25 kpc away from the Galactic center. The results show that the proposed method, free of widely used assumptions, can already be applied to pinpoint the dark matter distribution in the Milky Way with competitive accuracy, and paves the way for future developments.

Milky Way demographics with the VVV survey I. The 84-million star colour-magnitude diagram of the Galactic bulge

Context. The Milky Way (MW) bulge is a fundamental Galactic component for understanding the formation and evolution of galaxies, in particular our own. The ESO Public Survey VISTA Variables in the Via Lactea is a deep near-IR survey mapping the Galactic bulge and southern plane. Particularly for the bulge area, VVV is covering similar to 315 deg(2). Data taken during 2010 and 2011 covered the entire bulge area in the JHKs bands. Aims. We used VVV data for the whole bulge area as a single and homogeneous data set to build for the first time a single colour-magnitude diagram (CMD) for the entire Galactic bulge. Methods. Photometric data in the JHK(s) bands were combined to produce a single and huge data set containing 173 150 467 sources in the three bands, for the similar to 315 deg(2) covered by VVV in the bulge. Selecting only the data points flagged as stellar, the total number of sources is 84 095 284. Results. We built the largest colour-magnitude diagrams published up to date, containing 173.1+ million sources for all data points, and more than 84.0 million sources accounting for the stellar sources only. The CMD has a complex shape, mostly owing to the complexity of the stellar population and the effects of extinction and reddening towards the Galactic centre. The red clump (RC) giants are seen double in magnitude at b similar to -8 degrees-10 degrees, while in the inner part (b similar to -3 degrees) they appear to be spreading in colour, or even splitting into a secondary peak. Stellar population models show the predominance of main-sequence and giant stars. The analysis of the outermost bulge area reveals a well-defined sequence of late K and M dwarfs, seen at (J - K-s) similar to 0.7-0.9 mag and K-s greater than or similar to 14 mag. Conclusions. The interpretation of the CMD yields important information about the MW bulge, showing the fingerprint of its structure and content. We report a well-defined red dwarf sequence in the outermost bulge, which is important for the planetary transit searches of VVV. The double RC in magnitude seen in the outer bulge is the signature of the X-shaped MW bulge, while the spreading of the RC in colour, and even its splitting into a secondary peak, are caused by reddening effects. The region around the Galactic centre is harder to interpret because it is strongly affected by reddening and extinction.

The gaseous halo of the Milky Way

In the last decade, sensitive observations have revealed that disc galaxies are surrounded by multiphase gaseous halos produced by the circulation of gas from the discs to the environment and vice-versa. This Thesis is a study of the gaseous halo of the Milky Way carried out via the modelling of the HI emission and the available absorption-line data. We fitted simple kinematical models to the HI LAB Survey and found that the Galaxy has a massive (~3x10^8 Mo) HI halo extending a few kiloparsecs above the plane. This layer rotates more slowly than the disc and shows a global inflow motion, a kinematics similar to that observed in the HI halos of nearby galaxies. We built a dynamical model of the galactic fountain to reproduce the properties of this layer. In this model, fountain clouds are ejected from the disc by SN feedback and - as suggested by hydrodynamical simulations - triggers the cooling of coronal gas, which is entrained by the cloud wakes and accretes onto the disc when the clouds fall back. For a proper choice of the parameters, the model reproduces well the HI data and predicts an accretion of coronal gas onto the disc at a rate of 2 Mo/yr. We extended this model to the warm-hot component of the halo, showing that most of the ion absorption features observed towards background sources are consistent with being produced in the turbulent wakes that lag behind the fountain clouds. Specifically, the column densities, positions, and velocities of the absorbers are well reproduced by our model. Finally, we studied the gas content of galaxies extracted from a cosmological N-body+SPH simulation, and found that an HI halo with the forementioned properties is not observed, probably due ti the relatively low resolution of the simulations.

Resolved stellar populations in milky way and andromeda's new satellites

This thesis concerns the study of the variable stars and resolved stellar populations in four recently discovered dSphs, namely, Hercules and Ursa Major I (UMa I), which are UFD satellites of the MW; Andromeda XIX (And XIX) and Andromeda XXI (And XXI), which are satellites of M31. The main aim is to obtain detailed informations on the properties (age, metallicity, distance, and Oosterhoff type) of the stellar populations in these galaxies, to compare them with those of other satellites around the MW and M31, both ''classical'' dSphs and UFDs. The observables used to achieve these goals are the pulsating variables, especially the RR Lyrae stars, and the color magnitude diagram (CMD) of the resolved stellar populations. In particular, for UMa I, we combined B, V time-series observations from four different ground-based telescopes (Cassini, TLS, TT1 and Subaru) and for Hercules, we used archival data acquired with the Advanced Camera for Surveys (ACS) on board the HST. We used, instead B and V times-series photometry obtained with the Large Binocular Telescope (LBT) for And XIX and And XXI .