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.

Solar twins in M 67: evolutionary status and lithium abundance

Aims. We determine the age and mass of the three best solar twin candidates in open cluster M 67 through lithium evolutionary models. Methods. We computed a grid of evolutionary models with non-standard mixing at metallicity [Fe/H] = 0.01 with the Toulouse-Geneva evolution code for a range of stellar masses. We estimated the mass and age of 10 solar analogs belonging to the open cluster M 67. We made a detailed study of the three solar twins of the sample, YPB637, YPB1194, and YPB1787. Results. We obtained a very accurate estimation of the mass of our solar analogs in M 67 by interpolating in the grid of evolutionary models. The three solar twins allowed us to estimate the age of the open cluster, which is 3.87(-0.66)(+0.55) Gyr, which is better constrained than former estimates. Conclusions. Our results show that the 3 solar twin candidates have one solar mass within the errors and that M 67 has a solar age within the errors, validating its use as a solar proxy. M 67 is an important cluster when searching for solar twins.

Insights on the Milky Way bulge formation from the correlations between kinematics and metallicity

Context. Two main scenarios for the formation of the Galactic bulge are invoked, the first one through gravitational collapse or hierarchical merging of subclumps, the second through secular evolution of the Galactic disc. Aims. We aim to constrain the formation of the Galactic bulge through studies of the correlation between kinematics and metallicities in Baade's Window (l = 1 degrees, b = -4 degrees) and two other fields along the bulge minor axis (l = 0 degrees, b = -6 degrees and b = -12 degrees). Methods. We combine the radial velocity and the [Fe/H] measurements obtained with FLAMES/GIRAFFE at the VLT with a spectral resolution of R = 20 000, plus for the Baade's Window field the OGLE-II proper motions, and compare these with published N-body simulations of the Galactic bulge. Results. We confirm the presence of two distinct populations in Baade's Window found in Hill et al. (2010, A&A, submitted): the metal-rich population presents bar-like kinematics while the metal-poor population shows kinematics corresponding to an old spheroid or a thick disc. In this context the metallicity gradient along the bulge minor axis observed by Zoccali et al. (2008, A&A, 486, 177), visible also in the kinematics, can be related to a varying mix of these two populations as one moves away from the Galactic plane, alleviating the apparent contradiction between the kinematic evidence of a bar and the existence of a metallicity gradient. Conclusions. We show evidence that the two main scenarios for the bulge formation co-exist within the Milky Way bulge.

Beryllium abundances along the evolutionary sequence of the open cluster IC 4651-A new test for hydrodynamical stellar models

Context. Previous analyses of lithium abundances in main sequence and red giant stars have revealed the action of mixing mechanisms other than convection in stellar interiors. Beryllium abundances in stars with Li abundance determinations can offer valuable complementary information on the nature of these mechanisms. Aims. Our aim is to derive Be abundances along the whole evolutionary sequence of an open cluster. We focus on the well-studied open cluster IC 4651. These Be abundances are used with previously determined Li abundances, in the same sample stars, to investigate the mixing mechanisms in a range of stellar masses and evolutionary stages. Methods. Atmospheric parameters were adopted from a previous abundance analysis by the same authors. New Be abundances have been determined from high-resolution, high signal-to-noise UVES spectra using spectrum synthesis and model atmospheres. The careful synthetic modeling of the Be lines region is used to calculate reliable abundances in rapidly rotating stars. The observed behavior of Be and Li is compared to theoretical predictions from stellar models including rotation-induced mixing, internal gravity waves, atomic diffusion, and thermohaline mixing. Results. Beryllium is detected in all the main sequence and turn-off sample stars, both slow- and fast-rotating stars, including the Li-dip stars, but is not detected in the red giants. Confirming previous results, we find that the Li dip is also a Be dip, although the depletion of Be is more modest than for Li in the corresponding effective temperature range. For post-main-sequence stars, the Be dilution starts earlier within the Hertzsprung gap than expected from classical predictions, as does the Li dilution. A clear dispersion in the Be abundances is also observed. Theoretical stellar models including the hydrodynamical transport processes mentioned above are able to reproduce all the observed features well. These results show a good theoretical understanding of the Li and Be behavior along the color-magnitude diagram of this intermediate-age cluster for stars more massive than 1.2 M(circle dot).

CNO and F abundances in the globular cluster M22 (NGC 6656)

Context. Recent studies have confirmed the long standing suspicion that M 22 shares a metallicity spread and complex chemical enrichment history similar to that observed in omega Cen. M 22 is among the most massive Galactic globular clusters and its color-magnitude diagram and chemical abundances reveal the existence of sub-populations. Aims. To further constrain the chemical diversity of M 22, necessary to interpret its nucleosynthetic history, we seek to measure relative abundance ratios of key elements (carbon, nitrogen, oxygen, and fluorine) best studied, or only available, using high-resolution spectra at infrared wavelengths. Methods. High-resolution (R = 50 000) and high S/N infrared spectra were acquired of nine red giant stars with Phoenix at the Gemini-South telescope. Chemical abundances were calculated through a standard 1D local thermodynamic equilibrium analysis using Kurucz model atmospheres. Results. We derive [Fe/H] = -1.87 to -1.44, confirming at infrared wavelengths that M 22 does present a [Fe/H] spread. We also find large C and N abundance spreads, which confirm previous results in the literature but based on a smaller sample. Our results show a spread in A(C+N+O) of similar to 0.7 dex. Similar to mono-metallic globular clusters, M 22 presents a strong [Na/Fe]-[O/Fe] anticorrelation as derived from Na and CO lines in the K band. For the first time we recover F abundances in M 22 and find that it exhibits a 0.6 dex variation. We find tentative evidence for a flatter A(F)-A(O) relation compared to higher metallicity globular clusters. Conclusions. Our study confirms and expands upon the chemical diversity seen in this complex stellar system. All elements studied to date show large abundance spreads which require contributions from both massive and low mass stars.

Multiple stellar populations in globular clusters with photometry and low resolution spectroscopy

Our view of Globular Clusters has deeply changed in the last decade. Modern spectroscopic and photometric data have conclusively established that globulars are neither coeval nor monometallic, reopening the issue of the formation of such systems. Their formation is now schematized as a two-step process, during which the polluted matter from the more massive stars of a first generation gives birth, in the cluster innermost regions, to a second generation of stars with the characteristic signature of fully CNO-processed matter. To date, star-to-star variations in abundances of the light elements (C, N, O, Na) have been observed in stars of all evolutionary phases in all properly studied Galactic globular clusters. Multiple or broad evolutionary sequences have also been observed in nearly all the clusters that have been observed with good signal-to-noise in the appropriate photometric bands. The body of evidence suggests that spreads in light-element abundances can be fairly well traced by photometric indices including near ultraviolet passbands, as CNO abundance variations affect mainly wavelengths shorter than ~400 nm owing to the rise of some NH and CN molecular absorption bands. Here, we exploit this property of near ultraviolet photometry to trace internal chemical variations and combined it with low resolution spectroscopy aimed to derive carbon and nitrogen abundances in order to maximize the information on the multiple populations. This approach has been proven to be very effective in (i) detecting multiple population, (ii) characterizing their global properties (i.e., relative fraction of stars, location in the color-magnitude diagram, spatial distribution, and trends with cluster parameters) and (iii) precisely tagging their chemical properties (i.e., extension of the C-N anticorrelation, bimodalities in the N content).

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 .

The star formation and extinction coevolution of UV-selected galaxies over 0.05 < z < 1.2

We use a new stacking technique to obtain mean mid-IR and far-IR to far-UV flux ratios over the rest-frame near-UV, near-IR color-magnitude diagram. We employ COMBO-17 redshifts and COMBO-17 optical, GALEX far- and near-UV, and Spitzer IRAC and MIPS mid-IR photometry. This technique permits us to probe the infrared excess (IRX), the ratio of far-IR to far-UV luminosity, and the specific star formation rate (SSFR) and their coevolution over 2 orders of magnitude of stellar mass and over redshift 0.1 < z < 1.2. We find that the SSFR and the characteristic mass (Script M_0) above which the SSFR drops increase with redshift (downsizing). At any given epoch, the IRX is an increasing function of mass up to Script M_0. Above this mass the IRX falls, suggesting gas exhaustion. In a given mass bin below Script M_0, the IRX increases with time in a fashion consistent with enrichment. We interpret these trends using a simple model with a Schmidt-Kennicutt law and extinction that tracks gas density and enrichment. We find that the average IRX and SSFR follow a galaxy age parameter ξ, which is determined mainly by the galaxy mass and time since formation. We conclude that blue-sequence galaxies have properties which show simple, systematic trends with mass and time such as the steady buildup of heavy elements in the interstellar media of evolving galaxies and the exhaustion of gas in galaxies that are evolving off the blue sequence. The IRX represents a tool for selecting galaxies at various stages of evolution.

Resolved Stars' Insights into Galaxy Physics

Thesis (Ph.D.)--University of Washington, 2016-08

PFLEIDERER 2: IDENTIFICATION OF A NEW GLOBULAR CLUSTER IN THE GALAXY

We provide evidence that indicates the star cluster Pfleiderer 2, which is projected in a rich field, as a newly identified Galactic globular cluster. Since it is located in a crowded field, core extraction and decontamination tools were applied to reveal the cluster sequences in B, V, and I color-magnitude diagrams (CMDs). The main CMD features of Pfleiderer 2 are a tilted red giant branch and a red horizontal branch, indicating a high metallicity around solar. The reddening is E(B - V) = 1.01. The globular cluster is located at a distance of d(circle dot) = 16 +/- 2 kpc from the Sun. The cluster is located 2.7 kpc above the Galactic plane and at a distance of R(GC) = 9.7 kpc from the Galactic center, which is unusual for a metal-rich globular cluster.

Diagnostic tool to analyse colour-magnitude diagrams of poorly populated stellar concentrations

The dynamical processes that lead to open cluster disruption cause its mass to decrease. To investigate such processes from the observational point of view, it is important to identify open cluster remnants (OCRs), which are intrinsically poorly populated. Due to their nature, distinguishing them from field star fluctuations is still an unresolved issue. In this work, we developed a statistical diagnostic tool to distinguish poorly populated star concentrations from background field fluctuations. We use 2MASS photometry to explore one of the conditions required for a stellar group to be a physical group: to produce distinct sequences in a colour-magnitude diagram (CMD). We use automated tools to (i) derive the limiting radius; (ii) decontaminate the field and assign membership probabilities; (iii) fit isochrones; and (iv) compare object and field CMDs, considering the isochrone solution, in order to verify the similarity. If the object cannot be statistically considered as a field fluctuation, we derive its probable age, distance modulus, reddening and uncertainties in a self-consistent way. As a test, we apply the tool to open clusters and comparison fields. Finally, we study the OCR candidates DoDz 6, NGC 272, ESO 435 SC48 and ESO 325 SC15. The tool is optimized to treat these low-statistic objects and to separate the best OCR candidates for studies on kinematics and chemical composition. The study of the possible OCRs will certainly provide a deep understanding of OCR properties and constraints for theoretical models, including insights into the evolution of open clusters and dissolution rates.

VLT-FLAMES analysis of 8 giants in the bulge metal-poor globular cluster NGC 6522: oldest cluster in the Galaxy? Analysis of 8 giants in NGC 6522

Context. NGC 6522 has been the first metal-poor globular cluster identified in the bulge by Baade. Despite its importance, very few high-resolution abundance analyses of stars in this cluster are available. The bulge metal-poor clusters may be important tracers of the early chemical enrichment of the Galaxy. Aims. The main purpose of this study is to determine metallicity and elemental ratios in individual stars of NGC 6522. Methods. High-resolution spectra of 8 giants of the bulge's globular cluster NGC 6522 were obtained at the 8m VLT UT2-Kueyen telescope with the FLAMES+GIRAFFE spectrograph. Multiband V, I, J, K(s) photometry was used to derive effective temperatures as reference values. Spectroscopic parameters were derived from Fe I and Fe II lines, and adopted for the derivation of abundance ratios. Results. The present analysis provides a metallicity [Fe/H] = -1.0 +/- 0.2. The alpha-elements oxygen, magnesium and silicon show [O/Fe] = +0.4 +/- 0.3, [Mg/Fe] = [Si/Fe] = +0.25 +/- 0.15, whereas calcium and titanium show shallower ratios of [Ca/Fe] = [Ti/Fe] = +0.15 +/- 0.15. The neutron-capture r-process element europium appears to be overabundant by [Eu/Fe] = +0.4 +/- 0.4. The neutron-capture s-elements lanthanum and barium are enhanced by [La/Fe] = +0.35 +/- 0.2 and [Ba/Fe] = +0.5 +/- 0.5. The large internal errors, indicating the large star-to-star variation in the barium and europium abundances, are also discussed. Conclusions. The moderate metallicity combined to a blue horizontal branch (BHB), are characteristics similar to those of HP 1 and NGC 6558, pointing to a population of very old globular clusters in the Galactic bulge. Also, the abundance ratios in NGC 6522 resemble those in HP 1 and NGC 6558. The ultimate conclusion is that the bulge is old, and went through an early prompt chemical enrichment.

Recovery of the star formation history of the LMC from the VISTA survey of the Magellanic system

The VISTA near infrared survey of the Magellanic System (VMC) will provide deep YJK(s) photometry reaching stars in the oldest turn-off point throughout the Magellanic Clouds (MCs). As part of the preparation for the survey, we aim to access the accuracy in the star formation history (SFH) that can be expected from VMC data, in particular for the Large Magellanic Cloud (LMC). To this aim, we first simulate VMC images containing not only the LMC stellar populations but also the foreground Milky Way (MW) stars and background galaxies. The simulations cover the whole range of density of LMC field stars. We then perform aperture photometry over these simulated images, access the expected levels of photometric errors and incompleteness, and apply the classical technique of SFH-recovery based on the reconstruction of colour-magnitude diagrams (CMD) via the minimisation of a chi-squared-like statistics. We verify that the foreground MW stars are accurately recovered by the minimisation algorithms, whereas the background galaxies can be largely eliminated from the CMD analysis due to their particular colours and morphologies. We then evaluate the expected errors in the recovered star formation rate as a function of stellar age, SFR(t), starting from models with a known age-metallicity relation (AMR). It turns out that, for a given sky area, the random errors for ages older than similar to 0.4 Gyr seem to be independent of the crowding. This can be explained by a counterbalancing effect between the loss of stars from a decrease in the completeness and the gain of stars from an increase in the stellar density. For a spatial resolution of similar to 0.1 deg(2), the random errors in SFR(t) will be below 20% for this wide range of ages. On the other hand, due to the lower stellar statistics for stars younger than similar to 0.4 Gyr, the outer LMC regions will require larger areas to achieve the same level of accuracy in the SFR( t). If we consider the AMR as unknown, the SFH-recovery algorithm is able to accurately recover the input AMR, at the price of an increase of random errors in the SFR(t) by a factor of about 2.5. Experiments of SFH-recovery performed for varying distance modulus and reddening indicate that these parameters can be determined with (relative) accuracies of Delta(m-M)(0) similar to 0.02 mag and Delta E(B-V) similar to 0.01 mag, for each individual field over the LMC. The propagation of these errors in the SFR(t) implies systematic errors below 30%. This level of accuracy in the SFR(t) can reveal significant imprints in the dynamical evolution of this unique and nearby stellar system, as well as possible signatures of the past interaction between the MCs and the MW.

New Galactic star clusters discovered in the VVV survey

Context. VISTA Variables in the Via Lactea (VVV) is one of the six ESO Public Surveys operating on the new 4-m Visible and Infrared Survey Telescope for Astronomy (VISTA). VVV is scanning the Milky Way bulge and an adjacent section of the disk, where star formation activity is high. One of the principal goals of the VVV Survey is to find new star clusters of different ages. Aims. In order to trace the early epochs of star cluster formation we concentrated our search in the directions to those of known star formation regions, masers, radio, and infrared sources. Methods. The disk area covered by VVV was visually inspected using the pipeline processed and calibrated K(S)-band tile images for stellar over-densities. Subsequently, we examined the composite JHK(S) and ZJK(S) color images of each candidate. PSF photometry of 15 x 15 arcmin fields centered on the candidates was then performed on the Cambridge Astronomy Survey Unit reduced images. After statistical field-star decontamination, color-magnitude and color-color diagrams were constructed and analyzed. Results. We report the discovery of 96 new infrared open clusters and stellar groups. Most of the new cluster candidates are faint and compact (with small angular sizes), highly reddened, and younger than 5 Myr. For relatively well populated cluster candidates we derived their fundamental parameters such as reddening, distance, and age by fitting the solar-metallicity Padova isochrones to the color-magnitude diagrams.

Relevé spectroscopique et étude des propriétés physiques des étoiles naines blanches à moins de 40 parsecs du Soleil

Les étoiles naines blanches représentent la fin de l’évolution de 97% des étoiles de notre galaxie, dont notre Soleil. L’étude des propriétés globales de ces étoiles (distribution en température, distribution de masse, fonction de luminosité, etc.) requiert l’élaboration d’ensembles statistiquement complets et bien définis. Bien que plusieurs relevés d’étoiles naines blanches existent dans la littérature, la plupart de ceux-ci souffrent de biais statistiques importants pour ce genre d’analyse. L’échantillon le plus représentatif de la population d’étoiles naines blanches demeure à ce jour celui défini dans un volume complet, restreint à l’environnement immédiat du Soleil, soit à une distance de 20 pc (∼ 65 années-lumière) de celui-ci. Malheureusement, comme les naines blanches sont des étoiles intrinsèquement peu lumineuses, cet échantillon ne contient que ∼ 130 objets, compromettant ainsi toute étude statistique significative. Le but de notre étude est de recenser la population d’étoiles naines blanches dans le voisinage solaire a une distance de 40 pc, soit un volume huit fois plus grand. Nous avons ainsi entrepris de répertorier toutes les étoiles naines blanches à moins de 40 pc du Soleil à partir de SUPERBLINK, un vaste catalogue contenant le mouvement propre et les données photométriques de plus de 2 millions d’étoiles. Notre approche est basée sur la méthode des mouvements propres réduits qui permet d’isoler les étoiles naines blanches des autres populations stellaires. Les distances de toutes les candidates naines blanches sont estimées à l’aide de relations couleur-magnitude théoriques afin d’identifier les objets se situant à moins de 40 pc du Soleil, dans l’hémisphère nord. La confirmation spectroscopique du statut de naine blanche de nos ∼ 1100 candidates a ensuite requis 15 missions d’observations astronomiques sur trois grands télescopes à Kitt Peak en Arizona, ainsi qu’une soixantaine d’heures allouées sur les télescopes de 8 m des observatoires Gemini Nord et Sud. Nous avons ainsi découvert 322 nouvelles étoiles naines blanches de plusieurs types spectraux différents, dont 173 sont à moins de 40 pc, soit une augmentation de 40% du nombre de naines blanches connues à l’intérieur de ce volume. Parmi ces nouvelles naines blanches, 4 se trouvent probablement à moins de 20 pc du Soleil. De plus, nous démontrons que notre technique est très efficace pour identifier les étoiles naines blanches dans la région peuplée du plan de la Galaxie. Nous présentons ensuite une analyse spectroscopique et photométrique détaillée de notre échantillon à l’aide de modèles d’atmosphère afin de déterminer les propriétés physiques de ces étoiles, notamment la température, la gravité de surface et la composition chimique. Notre analyse statistique de ces propriétés, basée sur un échantillon presque trois fois plus grand que celui à 20 pc, révèle que nous avons identifié avec succès les étoiles les plus massives, et donc les moins lumineuses, de cette population qui sont souvent absentes de la plupart des relevés publiés. Nous avons également identifié plusieurs naines blanches très froides, et donc potentiellement très vieilles, qui nous permettent de mieux définir le côté froid de la fonction de luminosité, et éventuellement l’âge du disque de la Galaxie. Finalement, nous avons aussi découvert plusieurs objets d’intérêt astrophysique, dont deux nouvelles étoiles naines blanches variables de type ZZ Ceti, plusieurs naines blanches magnétiques, ainsi que de nombreux systèmes binaires non résolus.