Source author record

M. Schultheis

M. Schultheis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

51works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

51 published item(s)

preprint2026arXiv

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes in the underlying $N$-body simulation, as well as the relative positions of galaxies within their halo. Alignment strengths are moderated via stochastic misalignments, separately for central and satellite galaxies and conditional on the galaxy's redshift, luminosity, and rest-frame colour. The resulting model is calibrated against galaxy ellipticity statistics from the COSMOS Survey, selected alignment measurements based on Sloan Digital Sky Survey samples, and galaxy orientations extracted from the Horizon-AGN hydrodynamic simulation at redshift $z=1$. The best-fit model has a total of 12 alignment parameters and generally reproduces the calibration data sets well within the $1σ$ statistical uncertainties of the observations and the \flagship simulation, with notable exceptions for the most luminous sub-samples on small physical scales. The statistical power of the calibration data and the volume of the single \flagship realisation are still too small to provide informative prior ranges for intrinsic alignment amplitudes in relevant galaxy samples. As a first application, we predict that \Euclid end-of-mission tomographic weak gravitational lensing two-point statistics are modified by up to order $10\,\%$ due to intrinsic alignments.

preprint2026arXiv

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

We model intrinsic alignments (IA) in Euclid's Flagship simulation to investigate its impact on Euclid's weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Linear Alignment (NLA) and the Tidal Alignment and Tidal Torquing (TATT) models. We measure the amplitude of the simulated IA signal as a function of galaxy magnitude and colour in the redshift range $0.1<z<2.1$. We find that both NLA and TATT can accurately describe the IA signal in the simulation down to scales of $6$-$7 \,h^{-1}\,$Mpc. We measure alignment amplitudes for red galaxies comparable to those of the observations, with samples not used in the calibration procedure. For blue galaxies, our constraints are consistent with zero alignments in our first redshift bin $0.1 < z < 0.3$, but we detect a non-negligible signal at higher redshift, which is, however, consistent with the upper limits set by observational constraints. Additionally, several hydrodynamical simulations predict alignment for spiral galaxies, in agreement with our findings. Finally, the evolution of alignment with redshift is realistic and comparable to that determined in the observations. However, we find that the commonly adopted redshift power-law for IA fails to reproduce the simulation alignments above $z=1.1$. A significantly better agreement is obtained when a luminosity dependence is included, capturing the intrinsic luminosity evolution with redshift in magnitude-limited surveys. We conclude that the Flagship IA simulation is a useful tool for translating current IA constraints into predictions for IA contamination of Euclid-like samples.

preprint2025arXiv

Supergiant GCIRS 22 in the Milky Way Nuclear Star Cluster: Revised alpha abundances

Context: The chemical abundances of alpha-elements in Galactic Centre (GC) supergiants provide key insights into the chemical enrichment and star formation history of the Milky Way's Nuclear Star Cluster. Previous studies have reported enhanced alpha-element abundances, raising questions about the chemical evolution of this unique region. Aims: We aim to reassess the alpha-element abundances in the GC supergiant GCIRS 22 using updated spectral modelling and non-local thermodynamic equilibrium (NLTE) corrections to resolve discrepancies from earlier abundance analyses. Methods: High-resolution near-infrared spectra of GCIRS 22 were analysed using contemporary line lists and precise stellar parameters derived from scandium line diagnostics. We applied comprehensive NLTE corrections to accurately determine the abundances of silicon and calcium. Results: Our analysis reveals solar-scale alpha abundances ([Ca/Fe] = 0.06 $\pm$ 0.07; [Si/Fe] = $-$0.08 $\pm$ 0.20) for GCIRS 22, significantly lower than previous LTE-based findings. NLTE corrections reduce the calcium abundance by approximately 0.3,dex compared to LTE estimates, aligning our results with recent studies and highlighting the importance of accurate NLTE modelling. Conclusions: The solar-scale alpha-element abundances observed in GCIRS 22 suggest that recent star formation in the region has not been dominated by Type II supernovae, such as those expected from a recent starburst. Our findings support a scenario of episodic star formation, characterized by intermittent bursts separated by extended quiescent phases, or potentially driven by gas inflows from the inner disk, funnelled by the Galactic bar. Future comprehensive NLTE studies of additional GC stars will be essential for refining our understanding of the region's chemical evolution and star formation history.

preprint2022arXiv

Abundance analysis of APOGEE spectra for 58 metal-poor stars from the bulge spheroid

The central part of the Galaxy host a multitude of stellar populations, including the spheroidal bulge stars, stars moved to the bulge through secular evolution of the bar, inner halo, inner thick disk, inner thin disk, as well as debris from past accretion events. We identified a sample of 58 candidate stars belonging to the stellar population of the spheroidal bulge, and analyse their abundances. The present calculations of Mg, Ca, and Si lines are in agreement with the APOGEE-ASPCAP abundances, whereas abundances of C, N, O, and Ce are re-examined. We find normal $α$-element enhancements in oxygen, similar to magnesium, Si, and Ca abundances, which are typical of other bulge stars surveyed in the optical in Baade's Window. The enhancement of [O/Fe] in these stars suggests that they do not belong to accreted debris. No spread in N abundances is found, and none of the sample stars is N-rich, indicating that these stars are not second generation stars originated in globular clusters. Ce instead is enhanced in the sample stars, which points to an s-process origin such as due to enrichment from early generations of massive fast rotating stars, the so-called spinstars

preprint2022arXiv

Euclid preparation: XIII. Forecasts for galaxy morphology with the Euclid Survey using Deep Generative Models

We present a machine learning framework to simulate realistic galaxies for the Euclid Survey. The proposed method combines a control on galaxy shape parameters offered by analytic models with realistic surface brightness distributions learned from real Hubble Space Telescope observations by deep generative models. We simulate a galaxy field of $0.4\,\rm{deg}^2$ as it will be seen by the Euclid visible imager VIS and show that galaxy structural parameters are recovered with similar accuracy as for pure analytic Sérsic profiles. Based on these simulations, we estimate that the Euclid Wide Survey will be able to resolve the internal morphological structure of galaxies down to a surface brightness of $22.5\,\rm{mag}\,\rm{arcsec}^{-2}$, and $24.9\,\rm{mag}\,\rm{arcsec}^{-2}$ for the Euclid Deep Survey. This corresponds to approximately $250$ million galaxies at the end of the mission and a $50\,\%$ complete sample for stellar masses above $10^{10.6}\,\rm{M}_\odot$ (resp. $10^{9.6}\,\rm{M}_\odot$) at a redshift $z\sim0.5$ for the wide (resp. deep) survey. The approach presented in this work can contribute to improving the preparation of future high-precision cosmological imaging surveys by allowing simulations to incorporate more realistic galaxies.

preprint2022arXiv

Euclid Preparation: XIV. The Complete Calibration of the Color-Redshift Relation (C3R2) Survey: Data Release 3

The Complete Calibration of the Color-Redshift Relation (C3R2) survey is obtaining spectroscopic redshifts in order to map the relation between galaxy color and redshift to a depth of i ~ 24.5 (AB). The primary goal is to enable sufficiently accurate photometric redshifts for Stage IV dark energy projects, particularly Euclid and the Roman Space Telescope, which are designed to constrain cosmological parameters through weak lensing. We present 676 new high-confidence spectroscopic redshifts obtained by the C3R2 survey in the 2017B-2019B semesters using the DEIMOS, LRIS, and MOSFIRE multi-object spectrographs on the Keck telescopes. Combined with the 4454 redshifts previously published by this project, the C3R2 survey has now obtained and published 5130 high-quality galaxy spectra and redshifts. If we restrict consideration to only the 0.2 < z(phot) < 2.6 range of interest for the Euclid cosmological goals, then with the current data release C3R2 has increased the spectroscopic redshift coverage of the Euclid color space from 51% (as reported by Masters et al. 2015) to the current 91%. Once completed and combined with extensive data collected by other spectroscopic surveys, C3R2 should provide the spectroscopic calibration set needed to enable photometric redshifts to meet the cosmology requirements for Euclid, and make significant headway toward solving the problem for Roman.

preprint2022arXiv

Euclid preparation: XIX. Impact of magnification on photometric galaxy clustering

We investigate the importance of lensing magnification for estimates of galaxy clustering and its cross-correlation with shear for the photometric sample of Euclid. Using updated specifications, we study the impact of lensing magnification on the constraints and the shift in the estimation of the best fitting cosmological parameters that we expect if this effect is neglected. We follow the prescriptions of the official Euclid Fisher matrix forecast for the photometric galaxy clustering analysis and the combination of photometric clustering and cosmic shear. The slope of the luminosity function (local count slope), which regulates the amplitude of the lensing magnification, and the galaxy bias have been estimated from the Euclid Flagship simulation.We find that magnification significantly affects both the best-fit estimation of cosmological parameters and the constraints in the galaxy clustering analysis of the photometric sample. In particular, including magnification in the analysis reduces the 1$σ$ errors on $Ω_{\text{m},0}, w_{0}, w_a$ at the level of 20-35%, depending on how well we will be able to independently measure the local count slope. In addition, we find that neglecting magnification in the clustering analysis leads to shifts of up to 1.6$σ$ in the best-fit parameters. In the joint analysis of galaxy clustering, cosmic shear, and galaxy-galaxy lensing, magnification does not improve precision, but it leads to an up to 6$σ$ bias if neglected. Therefore, for all models considered in this work, magnification has to be included in the analysis of galaxy clustering and its cross-correlation with the shear signal ($3\times2$pt analysis) for an accurate parameter estimation.

preprint2022arXiv

Stellar ages, masses, extinctions and orbital parameters based on spectroscopic parameters of Gaia DR3

Gaia DR3 provides radial velocities for 33 million stars and spectroscopically derived atmospheric parameters for more than five million targets. When combined with the astrometric data, these allow us to derive orbital and stellar parameters that are key in order to understand the stellar populations of the Milky Way and perform galactic archaeology. We use the calibrated atmospheric parameters, 2MASS and Gaia-EDR3 photometry, and parallax-based distances to compute, via an isochrone fitting method, the ages, initial stellar masses and reddenings for the stars with spectroscopic parameters. We also derive the orbits (actions, eccentricities, apocentre, pericentre and Zmax) for all of the stars with measured radial velocities and astrometry, adopting two sets of line-of-sight distances from the literature and an axisymmetric potential of the Galaxy. Comparisons with reference catalogues of field and cluster stars suggest that reliable ages are obtained for stars younger than 9-10Gyr when the estimated relative age uncertainty is <50%. For older stars, ages tend to be under-estimated. The most reliable stellar type for age determination are turn-off stars, even when the input atmospheric parameters have large uncertainties. Ages for giants and main-sequence stars are retrieved with uncertainties of ~2Gyr when extinction towards the star's line-of sight is smaller than A_V<2.5mag. The full catalogue is made publicly available to be downloaded. With it, the full chemo-dynamical properties of the extended Solar neighbourhood unfold, and allow us to better identify the properties of the spiral arms, to parameterise the dynamical heating of the disc, or to thoroughly study the chemical enrichment of the Milky Way.

preprint2021arXiv

$Euclid$ preparation: XV. Forecasting cosmological constraints for the $Euclid$ and CMB joint analysis

The combination and cross-correlation of the upcoming $Euclid$ data with cosmic microwave background (CMB) measurements is a source of great expectation since it will provide the largest lever arm of epochs, ranging from recombination to structure formation across the entire past light cone. In this work, we present forecasts for the joint analysis of $Euclid$ and CMB data on the cosmological parameters of the standard cosmological model and some of its extensions. This work expands and complements the recently published forecasts based on $Euclid$-specific probes, namely galaxy clustering, weak lensing, and their cross-correlation. With some assumptions on the specifications of current and future CMB experiments, the predicted constraints are obtained from both a standard Fisher formalism and a posterior-fitting approach based on actual CMB data. Compared to a $Euclid$-only analysis, the addition of CMB data leads to a substantial impact on constraints for all cosmological parameters of the standard $Λ$-cold-dark-matter model, with improvements reaching up to a factor of ten. For the parameters of extended models, which include a redshift-dependent dark energy equation of state, non-zero curvature, and a phenomenological modification of gravity, improvements can be of the order of two to three, reaching higher than ten in some cases. The results highlight the crucial importance for cosmological constraints of the combination and cross-correlation of $Euclid$ probes with CMB data.

preprint2021arXiv

Euclid preparation: I. The Euclid Wide Survey

Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced straylight, producing very crisp images. This paper presents the building of the Euclid reference survey: the sequence of pointings of EWS, Deep fields, Auxiliary fields for calibrations, and spacecraft movements followed by Euclid as it operates in a step-and-stare mode from its orbit around the Lagrange point L2. Each EWS pointing has four dithered frames; we simulate the dither pattern at pixel level to analyse the effective coverage. We use up-to-date models for the sky background to define the Euclid region-of-interest (RoI). The building of the reference survey is highly constrained from calibration cadences, spacecraft constraints and background levels; synergies with ground-based coverage are also considered. Via purposely-built software optimized to prioritize best sky areas, produce a compact coverage, and ensure thermal stability, we generate a schedule for the Auxiliary and Deep fields observations and schedule the RoI with EWS transit observations. The resulting reference survey RSD_2021A fulfills all constraints and is a good proxy for the final solution. Its wide survey covers 14,500 square degrees. The limiting AB magnitudes ($5σ$ point-like source) achieved in its footprint are estimated to be 26.2 (visible) and 24.5 (near-infrared); for spectroscopy, the H$_α$ line flux limit is $2\times 10^{-16}$ erg cm$^{-2}$ s$^{-1}$ at 1600 nm; and for diffuse emission the surface brightness limits are 29.8 (visible) and 28.4 (near-infrared) mag arcsec$^{-2}$.

preprint2021arXiv

Euclid preparation: XVIII. Cosmic Dawn Survey. Spitzer observations of the Euclid deep fields and calibration fields

We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer's Infrared Array Camera (IRAC). We have combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11% of the total Spitzer mission time. The resulting mosaics cover a total of approximately 71.5deg$^2$ in the 3.6 and 4.5um bands, and approximately 21.8deg$^2$ in the 5.8 and 8um bands. They reach at least 24 AB magnitude (measured to sigma, in a 2.5 arcsec aperture) in the 3.6um band and up to ~ 5 mag deeper in the deepest regions. The astrometry is tied to the Gaia astrometric reference system, and the typical astrometric uncertainty for sources with 16<[3.6]<19 is <0.15 arcsec. The photometric calibration is in excellent agreement with previous WISE measurements. We have extracted source number counts from the 3.6um band mosaics and they are in excellent agreement with previous measurements. Given that the Spitzer Space Telescope has now been decommissioned these mosaics are likely to be the definitive reduction of these IRAC data. This survey therefore represents an essential first step in assembling multi-wavelength data on the Euclid deep fields which are set to become some of the premier fields for extragalactic astronomy in the 2020s.

preprint2020arXiv

Detailed abundances in the Galactic center: Evidence of a metal-rich alpha-enhanced stellar population

We present a detailed study of the composition of 20 M giants in the Galactic center with 15 of them confirmed to be in the Nuclear Star Cluster. As a control sample we have also observed 7 M giants in the Milky Way Disk with similar stellar parameters. All 27 stars are observed using the NIRSPEC spectograph on the KECK II telescope in the K-band at a resolving power of R=23,000. We report the first silicon abundance trends versus [Fe/H] for stars in the Galactic center. While finding a disk/bulge like trend at subsolar metallicities, we find that [Si/Fe] is enhanced at supersolar metallicities. We speculate on possible enrichment scenarios to explain such a trend. However, the sample size is modest and the result needs to be confirmed by additional measurements of silicon and other \textalpha-elements. We also derive a new distribution of [Fe/H] and find the most metal rich stars at [Fe/H]=+0.5 dex, confirming our earlier conclusions that the Galactic center hosts no stars with extreme chemical composition.

preprint2020arXiv

From the bulge to the outer disc: StarHorse stellar parameters, distances, and extinctions for stars in APOGEE DR16 and other spectroscopic surveys

We combine high-resolution spectroscopic data from APOGEE-2 Survey Data Release 16 (DR16) with broad-band photometric data from several sources, as well as parallaxes from {\it Gaia} Data Release 2 (DR2). Using the Bayesian isochrone-fitting code {\tt StarHorse}, we derive distances, extinctions and astrophysical parameters for around 388,815 APOGEE stars, achieving typical distance uncertainties of $\sim 6\%$ for APOGEE giants, $\sim 2\%$ for APOGEE dwarfs, as well as extinction uncertainties of $\sim 0.07$ mag when all photometric information is available, and $\sim 0.17$ mag if optical photometry is missing. {\tt StarHorse} uncertainties vary with the input spectroscopic catalogue, with the available photometry, and with the parallax uncertainties. To illustrate the impact of our results, we show that, thanks to {\it Gaia} DR2 and the now larger sky coverage of APOGEE-2 (including APOGEE-South), we obtain an extended map of the Galactic plane, providing unprecedented coverage of the disk close to the Galactic mid-plane ($|Z_{Gal}|<1$ kpc) from the Galactic Centre out to $R_{\rm Gal}\sim 20$ kpc. The improvements in statistics as well as distance and extinction uncertainties unveil the presence of the bar in stellar density, as well as the striking chemical duality in the innermost regions of the disk, now clearly extending to the inner bulge. We complement this paper with distances and extinctions for stars in other public released spectroscopic surveys: 324,999 in GALAH DR2, 4,928,715 in LAMOST DR5, 408,894 in RAVE DR6, and 6,095 in GES DR3

preprint2020arXiv

Heavy element evolution in the inner regions of the Milky Way

We present results for the evolution of the abundances of heavy elements (O, Mg, Al, Si, K, Ca, Cr, Mn, Ni and Fe) in the inner Galactic regions ($R_{GC} < 4$kpc). We adopt a detailed chemical evolution model already tested for the Galactic bulge and compare the results with APOGEE data. We start with a set of yields from the literature which are considered the best to reproduce the abundance patterns in the solar vicinity. We find that in general the predicted trends nicely reproduce the data but in some cases either the trend or the absolute values of the predicted abundances need to be corrected, even by large factors, in order to reach the best agreement. We suggest how the current stellar yields should be modified to reproduce the data and we discuss whether such corrections are reasonable in the light of the current knowledge of stellar nucleosynthesis. However, we also critically discuss the observations. Our results suggest that Si, Ca, Cr and Ni are the elements for which the required corrections are the smallest, while for Mg and Al moderate modifications are necessary. On the other hand, O and K need the largest corrections to reproduce the observed patterns, a conclusion already reached for solar vicinity abundance patterns, with the exception of oxygen. For Mn we apply corrections already suggested in previous works. \end{abstract}

preprint2020arXiv

MOONS Surveys of the Milky Way and its Satellites

The study of resolved stellar populations in the Milky Way and other Local Group galaxies can provide us with a fossil record of their chemo-dynamical and star-formation histories over timescales of many billions of years. In the galactic components and stellar systems of the Milky Way and its satellites, individual stars can be resolved. Therefore, they represent a unique laboratory in which to investigate the details of the processes behind the formation and evolution of the disc and dwarf/irregular galaxies. MOONS at the VLT represents a unique combination of an efficient infrared multi-object spectrograph and a large-aperture 8-m-class telescope which will sample the cool stellar populations of the dense central regions of the Milky Way and its satellites, delivering accurate radial velocities, metallicities, and other chemical abundances for several millions of stars over its lifetime (see Cirasuolo et al., this issue). MOONS will observe up to 1000 targets across a 25-arcminute field of view in the optical and near-infrared (0.6-1.8 micron) simultaneously. A high-resolution (R~19700) setting in the H band has been designed for the accurate determination of stellar abundances such as alpha, light, iron-peak and neutron-capture elements.

preprint2020arXiv

Optical interferometry and Gaia measurement uncertainties reveal the physics of asymptotic giant branch stars

Context. Asymptotic giant branch stars are cool luminous evolved stars that are well observable across the Galaxy and populating Gaia data. They have complex stellar surface dynamics Aims. On the AGB star CL Lac, it has been shown that the convection-related variability accounts for a substantial part of the Gaia DR2 parallax error. We observed this star with the MIRC-X beam combiner installed at the CHARA interferometer to detect the presence of stellar surface inhomogeneities. Methods. We performed the reconstruction of aperture synthesis images from the interferometric observations at different wavelengths. Then, we used 3D radiative hydrodynamics simulations of stellar convection with CO5BOLD and the post-processing radiative transfer code Optim3D to compute intensity maps in the spectral channels of MIRC-X observations. Then, we determined the stellar radius and compared the 3D synthetic maps to the reconstructed ones focusing on matching the intensity contrast, the morphology of stellar surface structures, and the photocentre position at two different spectral channels, 1.52 and 1.70 micron, simultaneously. Results. We measured the apparent diameter of CL Lac at two wavelengths and recovered the radius using a Gaia parallax. In addition to this, the reconstructed images are characterised by the presence of a brighter area that largely affects the position of the photocentre. The comparison with 3D simulation shows good agreement with the observations both in terms of contrast and surface structure morphology, meaning that our model is adequate for explaining the observed inhomogenities. Conclusions. This work confirms the presence of convection-related surface structures on an AGB star of Gaia DR2. Our result will help us to take a step forward in exploiting Gaia measurement uncertainties to extract the fundamental properties of AGB stars using appropriate RHD simulations.

preprint2020arXiv

The diffuse interstellar band around 8620 Å I. Methods and application to the GIBS data set

We developed a set of procedures to automatically detect and measure the DIB around 8620 Å (the Gaia DIB) for a wide range of temperatures. The DIB profile is fit with a Gaussian function. Specifically, the DIB feature is extracted from the spectra of late-type stars by subtracting the corresponding synthetic spectra. For early-type stars we applied a specific model based on the Gaussian process that needs no prior knowledge of the stellar parameters. The method was tested on $\sim$5000 spectra from the Giraffe Inner Bulge Survey (GIBS). After validation, we obtained 4194 reasonable fitting results from the GIBS database. An EW versus $E(J\,{-}\,K_{\rm S})$ relation is derived as $E(J\,{-}\,K_{\rm S})\,{=}\,1.875\,({\pm}\,0.152)\,{\times}\,{\rm EW}\,{-}\,0.011\,({\pm}\,0.048)$, according to $E(B\,{-}\,V)/{\rm EW}\,{=}\,2.721$, which is highly consistent with previous results toward similar sightlines. After a correction based on the VVV database for both EW and reddening, the coefficient derived from individual GIBS fields, $E(J\,{-}\,K_{\rm S})/{\rm EW}\,{=}\,1.884\,{\pm}\,0.225$, is also in perfect agreement with literature values. Based on a subsample of 1015 stars toward the Galactic center within $-3^{\circ}\,{<}\,b\,{<}\,3^{\circ}$ and $-6^{\circ}\,{<}\,l\,{<}\,3^{\circ}$, we determined a rest-frame wavelength of the Gaia DIB as 8620.55 Å. A Gaussian profile is proved to be a proper and stable assumption for the Gaia DIB as no intrinsic asymmetry is found.

preprint2016arXiv

Chemical trends in the Galactic Halo with APOGEE data

The galaxy formation process in the $Λ$-Cold Dark Matter scenario can be constrained from the analysis of stars in the Milky Way's halo system. We examine the variation of chemical abundances in distant halo stars observed by the Apache Point Galactic Evolution Experiment (APOGEE), as a function of distance from the Galactic center ($r$) and iron abundance ([M/H]), in the range 5 $\lesssim r \lesssim$ 30 kpc and $-2.5 <$ [M/H] $<$ 0.0. We perform a statistical analysis of the abundance ratios derived by the APOGEE pipeline (ASPCAP) and distances calculated by several approaches. Our analysis reveals signatures of a different chemical enrichment between the inner and outer regions of the halo, with a transition at about 15 kpc. The derived metallicity distribution function exhibits two peaks, at [M/H] $\sim -1.5$ and $\sim -2.1$, consistent with previously reported halo metallicity distributions. We obtain a difference of $\sim 0.1$ dex for $α$-element-to-iron ratios for stars at $r > 15$ kpc and [M/H] $> -1.1$ (larger in the case of O, Mg and S) with respect to the nearest halo stars. This result confirms previous claims for low-$α$ stars found at larger distances. Chemical differences in elements with other nucleosynthetic origins (Ni, K, Na, and Al) are also detected. C and N do not provide reliable information about the interstellar medium from which stars formed because our sample comprises RGB and AGB stars and can experience mixing of material to their surfaces.

preprint2016arXiv

Discovery of a Metal-Poor Field Giant with a Globular Cluster Second-Generation Abundance Pattern

We report on detection, from observations obtained with the APOGEE spectroscopic survey, of a metal-poor ([Fe/H] $= -1.3$ dex) field giant star with an extreme Mg-Al abundance ratio ([Mg/Fe] $= -0.31$ dex; [Al/Fe] $= 1.49$ dex). Such low Mg/Al ratios are seen only among the second-generation population of globular clusters, and are not present among Galactic disk field stars. The light element abundances of this star, 2M16011638-1201525, suggest that it could have been born in a globular cluster. We explore several origin scenarios, in particular studying the orbit of the star to check the probability of it being kinematically related to known globular clusters. We performed simple orbital integrations assuming the estimated distance of 2M16011638-1201525 and the available six-dimensional phase-space coordinates of 63 globular clusters, looking for close encounters in the past with a minimum distance approach within the tidal radius of each cluster. We found a very low probability that 2M16011638-1201525 was ejected from most globular clusters; however, we note that the best progenitor candidate to host this star is globular cluster $ω$ Centauri (NGC 5139). Our dynamical investigation demonstrates that 2M16011638-1201525 reaches a distance $|Z_{max}| < 3 $ kpc from the Galactic plane and a minimum and maximum approach to the Galactic center of $R_{min}<0.62$ kpc and $R_{max}<7.26$ kpc in an eccentric ($e\sim0.53$) and retrograde orbit. Since the extreme chemical anomaly of 2M16011638-1201525 has also been observed in halo field stars, this object could also be considered a halo contaminant, likely been ejected into the Milky Way disk from the halo. We conclude that, 2M16011638-20152 is also kinematically consistent with the disk but chemically consistent with halo field stars.

preprint2016arXiv

Galactic Archaeology with asteroseismology and spectroscopy: Red giants observed by CoRoT and APOGEE

With the advent of the space missions CoRoT and Kepler, it has become feasible to determine precise asteroseismic masses and ages for large samples of red-giant stars. In this paper, we present the CoRoGEE dataset -- obtained from CoRoT lightcurves for 606 red giant stars in two fields of the Galactic disc which have been co-observed for an ancillary project of APOGEE. We have used the Bayesian parameter estimation code PARAM to calculate distances, extinctions, masses, and ages for these stars in a homogeneous analysis, resulting in relative statistical uncertainties of $\sim2\%$ in distance, $\sim4\%$ in radius, $\sim9\%$ in mass and $\sim25\%$ in age. We also assess systematic age uncertainties due to different input physics and mass loss. We discuss the correlation between ages and chemical abundance patterns of field stars over a large radial range of the Milky Way's disc (5 kpc $<R_{\rm Gal}<$ 14 kpc), focussing on the [$α$/Fe]-[Fe/H]-age plane in five radial bins of the Galactic disc. We find an overall agreement with the expectations of chemical-evolution models computed before the present data were available, especially for the outer regions. However, our data also indicate that a significant fraction of stars now observed near and beyond the Solar Neighbourhood migrated from inner regions. Mock CoRoGEE observations of a chemo-dynamical Milky Way disc model show that the number of high-metallicity stars in the outer disc is too high to be accounted for even by the strong radial mixing present in the model. The mock observations also reveal that the age distribution of the [$α$/Fe]-enhanced sequence in the CoRoGEE inner-disc field is much broader than expected from a combination of radial mixing and observational errors. We suggest that a thick disc/bulge component that formed stars for more than 3 Gyr may account for these discrepancies.

preprint2016arXiv

Temperatures and metallicities of M giants in the galactic Bulge from low-resolution K-band spectra

With the existing and upcoming large multi-fibre low-resolution spectrographs, the question arises how precise stellar parameters such as Teff and [Fe/H] can be obtained from low-resolution K-band spectra with respect to traditional photometric temperature measurements. Until now, most of the effective temperatures in galactic Bulge studies come directly from photometric techniques. Uncertainties in interstellar reddening and in the assumed extinction law could lead to large systematic errors. We aim to obtain and calibrate the relation between Teff and the $\rm ^{12}CO$ first overtone bands for M giants in the galactic Bulge covering a wide range in metallicity. We use low-resolution spectra for 20 M giants with well-studied parameters from photometric measurements covering the temperature range 3200 < Teff < 4500 K and a metallicity range from 0.5 dex down to -1.2 dex and study the behaviour of Teff and [Fe/H] on the spectral indices. We find a tight relation between Teff and the $\rm ^{12}CO(2-0)$ band with a dispersion of 95 K as well as between Teff and the $\rm ^{12}CO(3-1)$ with a dispersion of 120 K. We do not find any dependence of these relations on the metallicity of the star, making them relation attractive for galactic Bulge studies. This relation is also not sensitive to the spectral resolution allowing to apply this relation in a more general way. We also found a correlation between the combination of the NaI, CaI and the $\rm ^{12}CO$ band with the metallicity of the star. However this relation is only valid for sub-solar metallicities. We show that low-resolution spectra provide a powerful tool to obtain effective temperatures of M giants. We show that this relation does not depend on the metallicity of the star within the investigated range and is also applicable to different spectral resolution.

preprint2015arXiv

Chemical evolution of the Galactic Center

In recent years, the Galactic Center (GC) region (200 pc in radius) has been studied in detail with spectroscopic stellar data as well as an estimate of the ongoing star formation rate. The aims of this paper are to study the chemical evolution of the GC region by means of a detailed chemical evolution model and to compare the results with high resolution spectroscopic data in order to impose constraints on the GC formation history.The chemical evolution model assumes that the GC region formed by fast infall of gas and then follows the evolution of alpha-elements and Fe. We test different initial mass functions (IMFs), efficiencies of star formation and gas infall timescales. To reproduce the currently observed star formation rate, we assume a late episode of star formation triggered by gas infall/accretion. We find that, in order to reproduce the [alpha/Fe] ratios as well as the metallicity distribution function observed in GC stars, the GC region should have experienced a main early strong burst of star formation, with a star formation efficiency as high as 25 Gyr^{-1}, occurring on a timescale in the range 0.1-0.7 Gyr, in agreement with previous models of the entire bulge. Although the small amount of data prevents us from drawing firm conclusions, we suggest that the best IMF should contain more massive stars than expected in the solar vicinity, and the last episode of star formation, which lasted several hundred million years, should have been triggered by a modest episode of gas infall/accretion, with a star formation efficiency similar to that of the previous main star formation episode. This last episode of star formation produces negligible effects on the abundance patterns and can be due to accretion of gas induced by the bar. Our results exclude an important infall event as a trigger for the last starburst.

preprint2015arXiv

Evidence for a metal-poor population in the inner Galactic Bulge

The inner Galactic Bulge has, until recently, been avoided in chemical evolution studies due to extreme extinction and stellar crowding. Large, near-IR spectroscopic surveys, such as APOGEE, allow for the first time the measurement of metallicities in the inner region of our Galaxy. We study metallicities of 33 K/M giants situated in the Galactic Center region from observations obtained with the APOGEE survey. We selected K/M giants with reliable stellar parameters from the APOGEE/ASPCAP pipeline. Distances, interstellar extinction values, and radial velocities were checked to confirm that these stars are indeed situated in the inner Galactic Bulge. We find a metal-rich population centered at [M/H] = +0.4 dex, in agreement with earlier studies of other bulge regions, but also a peak at low metallicity around $\rm [M/H] = -1.0\,dex$, suggesting the presence of a metal-poor population which has not previously been detected in the central region. Our results indicate a dominant metal-rich population with a metal-poor component that is enhanced in the $α$-elements. This metal-poor population may be associated with the classical bulge and a fast formation scenario.

preprint2015arXiv

The Gaia-ESO Survey: Separating disk chemical substructures with cluster models

(Abridged) Recent spectroscopic surveys have begun to explore the Galactic disk system outside the solar neighborhood on the basis of large data samples. In this way, they provide valuable information for testing spatial and temporal variations of disk structure kinematics and chemical evolution. We used a Gaussian mixture model algorithm, as a rigurous mathematical approach, to separate in the [Mg/Fe] vs. [Fe/H] plane a clean disk star subsample from the Gaia-ESO survey internal data release 2. We find that the sample is separated into five groups associated with major Galactic components; the metal-rich end of the halo, the thick disk, and three subgroups for the thin disk sequence. This is confirmed with a sample of red clump stars from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. The two metal-intermediate and metal-rich groups of the thin disk decomposition ([Fe/H]>-0.25 dex) highlight a change in the slope at solar metallicity. This holds true at different radial regions. The distribution of Galactocentric radial distances of the metal-poor part of the thin disk ([Fe/H]<-0.25 dex) is shifted to larger distances than those of the more metal-rich parts. Moreover, the metal-poor part of the thin disk presents indications of a scale height intermediate between those of the thick and the rest of the thin disk, and it displays higher azimuthal velocities than the latter. These stars might have formed and evolved in parallel and/or dissociated from the inside-out formation taking place in the internal thin disk. Their enhancement levels might be due to their origin from gas pre-enriched by outflows from the thick disk or the inner halo. The smooth trends of their properties (their spatial distribution with respect to the plane, in particular) with [Fe/H] and [Mg/Fe] suggested by the data indicates a quiet dynamical evolution, with no relevant merger events.

preprint2015arXiv

The Gaia-ESO Survey: Tracing interstellar extinction

Large spectroscopic surveys have enabled in the recent years the computation of three-dimensional interstellar extinction maps thanks to accurate stellar atmospheric parameters and line-of-sight distances. Such maps are complementary to 3D maps extracted from photometry, allowing a more thorough study of the dust properties. Our goal is to use the high-resolution spectroscopic survey Gaia-ESO in order to obtain with a good distance resolution the interstellar extinction and its dependency as a function of the environment and the Galactocentric position. We use the stellar atmospheric parameters of more than 5000 stars, obtained from the Gaia-ESO survey second internal data release, and combine them with optical (SDSS) and near-infrared (VISTA) photometry as well as different sets of theoretical stellar isochrones, in order to calculate line-of-sight extinction and distances. The extinction coefficients are then compared with the literature to discuss their dependancy on the stellar parameters and position in the Galaxy. Within the errors of our method, our work does not show that there is any dependence of the interstellar extinction coefficient on the atmospheric parameters of the stars. We do not find any evidence of the variation of E(J-H)/E(J-K) with the angle from the Galactic centre nor with Galactocentric distance. This suggests that we are dealing with a uniform extinction law in the SDSS ugriz bands and the near-IR JHKs bands. Therefore, extinction maps using mean colour-excesses and assuming a constant extinction coefficient can be used without introducing any systematic errors.

preprint2015arXiv

Understanding AGB evolution in Galactic bulge stars from high-resolution infrared spectroscopy

An analysis of high-resolution near-infrared spectra of a sample of 45 asymptotic giant branch (AGB) stars towards the Galactic bulge is presented. The sample consists of two subsamples, a larger one in the inner and intermediate bulge, and a smaller one in the outer bulge. The data are analysed with the help of hydrostatic model atmospheres and spectral synthesis. We derive the radial velocity of all stars, and the atmospheric chemical mix ([Fe/H], C/O, $^{12}$C/$^{13}$C, Al, Si, Ti, and Y) where possible. Our ability to model the spectra is mainly limited by the (in)completeness of atomic and molecular line lists, at least for temperatures down to $T_{\rm eff}\approx3100$ K. We find that the subsample in the inner and intermediate bulge is quite homogeneous, with a slightly sub-solar mean metallicity and only few stars with super-solar metallicity, in agreement with previous studies of non-variable M-type giants in the bulge. All sample stars are oxygen-rich, C/O$<$1.0. The C/O and carbon isotopic ratios suggest that third dredge-up (3DUP) is absent among the sample stars, except for two stars in the outer bulge that are known to contain technetium. These stars are also more metal-poor than the stars in the intermediate or inner bulge. Current stellar masses are determined from linear pulsation models. The masses, metallicities and 3DUP behaviour are compared to AGB evolutionary models. We conclude that these models are partly in conflict with our observations. Furthermore, we conclude that the stars in the inner and intermediate bulge belong to a more metal-rich population that follows bar-like kinematics, whereas the stars in the outer bulge belong to the metal-poor, spheroidal bulge population.

preprint2015arXiv

Young [$α$/Fe]-enhanced stars discovered by CoRoT and APOGEE: What is their origin?

We report the discovery of a group of apparently young CoRoT red-giant stars exhibiting enhanced [alpha/Fe] abundance ratios (as determined from APOGEE spectra) with respect to Solar values. Their existence is not explained by standard chemical evolution models of the Milky Way, and shows that the chemical-enrichment history of the Galactic disc is more complex. We find similar stars in previously published samples for which isochrone-ages could be robustly obtained, although in smaller relative numbers, which could explain why these stars have not received prior attention. The young [alpha/Fe]-rich stars are much more numerous in the CoRoT-APOGEE (CoRoGEE) inner-field sample than in any other high-resolution sample available at present, as only CoRoGEE can explore the inner-disc regions and provide ages for its field stars. The kinematic properties of the young [$α$/Fe]-rich stars are not clearly thick-disc like, despite their rather large distances from the Galactic mid-plane. Our tentative interpretation of these and previous intriguing observations in the Milky Way is that these stars were formed close to the end of the Galactic bar, near corotation -- a region where gas can be kept inert for longer times, compared to other regions shocked more frequently by the passage of spiral arms. Moreover, that is where the mass return from older inner-disc stellar generations should be maximal (according to an inside-out disc-formation scenario), further diluting the in-situ gas. Other possibilities to explain these observations (e.g., a recent gas-accretion event) are also discussed.

preprint2014arXiv

A high resolution spectroscopic atlas of M subdwarfs - Effective temperature and metallicity

Context. M subdwarfs are metal poor and cool stars. They are important probes of the old galactic populations. However, they remain elusive due to their low luminosity. Observational and modeling efforts are required to fully understand their physics and to investigate the effects of metallicity in their cool atmospheres. Aims. We perform a detailed study of a sample of subdwarfs to determine their stellar parameters and constrain the stat-of-the art atmospheric models. Methods. We present UVES/VLT high resolution spectra of three late-K subdwarfs and 18 M subdwarfs. Our atlas covers the optical region from 6400 AA up to the near infrared at 8900 AA. We show spectral details of cool atmospheres at very high resolution (R= 40 000) and compare with synthetic spectra computed from the recent BT-Settl atmosphere models. Results. Our comparison shows that molecular features (TiO, VO, CaH), and atomic features (Fe I, Ti I, Na I, K I) are well fitted by current models. We produce an effective temperature versus spectral type relation all over the subdwarf spectral sequence. Thanks to the high resolution of our spectra, we perform a detailed comparison of line profiles of individual elements such as Fe I, Ca II, Ti I, and are able to determine accurate metallicities of these stars. These determinations contribute to calibrate the relation between metallicity and molecular band strength indices from low-resolution spectra. Conclusions. This work shows that the new generation of models are able to reproduce various spectral features of M subdwarfs. Working with these high resolution spectra allowed us to disentangle the atmospheric parameters (effective temperature, gravity, metallicity), which is not possible when using low resolution spectroscopy or photometry.

preprint2014arXiv

A three dimensional extinction map of the Galactic Anticentre from multi-band photometry

We present a three dimensional extinction map in $r$ band. The map has a spatial angular resolution, depending on latitude, between 3 -- 9\,arcmin and covers the entire XSTPS-GAC survey area of over 6,000\,$\rm deg^2$ for Galactic longitude $\rm 140 \leq$ $l$ $ \leq 220°$ and latitude $\rm -40\leq$ $b$ $ \leq 40°$. By cross-matching the photometric catalog of the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC) with those of 2MASS and WISE, we have built a multi-band photometric stellar sample of about 30 million stars and applied spectral energy distribution (SED) fitting to the sample. By combining photometric data from the optical to the near-infrared, we are able to break the degeneracy between the intrinsic stellar colours and the amounts of extinction by dust grains for stars with high photometric accuracy, and trace the extinction as a function of distance for low Galactic latitude and thus highly extincted regions. This has allowed us to derive the best-fit extinction and distance information of more than 13 million stars, which are used to construct the three dimensional extinction map. We have also applied a Rayleigh-Jeans colour excess (RJCE) method to the data using the 2MASS and WISE colour $(H-W2)$. The resulting RJCE extinction map is consistent with the integrated two dimensional map deduced using the best-fit SED algorithm. However for individual stars, the amounts of extinction yielded by the RJCE method suffer from larger errors than those given by the best-fit SED algorithm.

preprint2014arXiv

Chemical Abundances of M giants in the Galactic Center: a Single Metal-Rich Population with Low [alpha/Fe]

Context. The formation and evolution of the Milky Way bulge is still largely an unanswered question. One of the most essential observables needed in its modelling are the metallicity distribution and the trends of the alpha elements as measured in stars. While Bulge regions beyond R > 50 pc of the centre has been targeted in several surveys, the central part has escaped detailed study due to the extreme extinction and crowding. The abundance gradients from the center are, however, of large diagnostic value. Aims. We aim at investigating the Galactic Centre environment by probing M giants in the field, avoiding supergiants and cluster members. Methods. For 9 field M-giants in the Galactic Centre region, we have obtained high- and low-resolution spectra observed simultaneously with CRIRES and ISAAC on UT1 and UT3 of the VLT. The low-resolution spectra provide a means of determining the effective temperatures, and the high-resolution spectra provide detailed abundances of Fe, Mg, Si, and Ca. Results. We find a metal-rich population at [Fe/H]=+0.11+-0.15 and a lack of the metal-poor population, found further out in the Bulge, corroborating earlier studies. Our [alpha/Fe] element trends, however, show low values, following the outer Bulge trends. A possible exception of the [Ca/Fe] trend is found and needs further investigation. Conclusions. The results of the analysed field M-giants in the Galactic Centre region, excludes a scenario with rapid formation, in which SNIIe played a dominated role in the chemical enrichment of the gas. The metal-rich metallicities together with low alpha-enhancement seems to indicate a bar-like population perhaps related to the nuclear bar.

preprint2014arXiv

Chemodynamics of the Milky Way. I. The first year of APOGEE data

We investigate the chemo-kinematic properties of the Milky Way disc by exploring the first year of data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE), and compare our results to smaller optical high-resolution samples in the literature, as well as results from lower resolution surveys such as GCS, SEGUE and RAVE. We start by selecting a high-quality sample in terms of chemistry ($\sim$ 20.000 stars) and, after computing distances and orbital parameters for this sample, we employ a number of useful subsets to formulate constraints on Galactic chemical and chemodynamical evolution processes in the Solar neighbourhood and beyond (e.g., metallicity distributions -- MDFs, [$α$/Fe] vs. [Fe/H] diagrams, and abundance gradients). Our red giant sample spans distances as large as 10 kpc from the Sun. We find remarkable agreement between the recently published local (d $<$ 100 pc) high-resolution high-S/N HARPS sample and our local HQ sample (d $<$ 1 kpc). The local MDF peaks slightly below solar metallicity, and exhibits an extended tail towards [Fe/H] $= -$1, whereas a sharper cut-off is seen at larger metallicities. The APOGEE data also confirm the existence of a gap in the [$α$/Fe] vs. [Fe/H] abundance diagram. When expanding our sample to cover three different Galactocentric distance bins, we find the high-[$α$/Fe] stars to be rare towards the outer zones, as previously suggested in the literature. For the gradients in [Fe/H] and [$α$/Fe], measured over a range of 6 $ < $ R $ <$ 11 kpc in Galactocentric distance, we find a good agreement with the gradients traced by the GCS and RAVE dwarf samples. For stars with 1.5 $<$ z $<$ 3 kpc, we find a positive metallicity gradient and a negative gradient in [$α$/Fe].

preprint2014arXiv

Extinction Maps towards the Milky Way Bulge: 2D and 3D Tests with APOGEE

Galactic interstellar extinction maps are powerful and necessary tools for Milky Way structure and stellar population analyses, particularly toward the heavily-reddened bulge and in the midplane. However, due to the difficulty of obtaining reliable extinction measures and distances for a large number of stars that are independent of these maps, tests of their accuracy and systematics have been limited. Our goal is to assess a variety of photometric stellar extinction estimates, including both 2D and 3D extinction maps, using independent extinction measures based on a large spectroscopic sample of stars towards the Milky Way bulge. We employ stellar atmospheric parameters derived from high-resolution $H$-band APOGEE spectra, combined with theoretical stellar isochrones, to calculate line-of-sight extinction and distances for a sample of more than 2400 giants towards the Milky Way bulge. We compare these extinction values to those predicted by individual near-IR and near+mid-IR stellar colors, 2D bulge extinction maps and 3D extinction maps. The long baseline, near+mid-IR stellar colors are, on average, the most accurate predictors of the APOGEE extinction estimates, and the 2D and 3D extinction maps derived from different stellar populations along different sightlines show varying degrees of reliability. We present the results of all of the comparisons and discuss reasons for the observed discrepancies. We also demonstrate how the particular stellar atmospheric models adopted can have a strong impact on this type of analysis, and discuss related caveats.

preprint2014arXiv

Mapping the Milky Way bulge at high resolution: the 3D dust extinction, CO, and X factor maps

Three dimensional interstellar extinction maps provide a powerful tool for stellar population analysis. We use data from the VISTA Variables in the Via Lactea survey together with the Besançon stellar population synthesis model of the Galaxy to determine interstellar extinction as a function of distance in the Galactic bulge covering $ -10 < l < 10$ and $-10 < b <5$. We adopted a recently developed method to calculate the colour excess. First we constructed the H-Ks vs. Ks and J-Ks vs. Ks colour-magnitude diagrams based on the VVV catalogues that matched 2MASS. Then, based on the temperature-colour relation for M giants and the distance-colour relations, we derived the extinction as a function of distance. The observed colours were shifted to match the intrinsic colours in the Besançon model as a function of distance iteratively. This created an extinction map with three dimensions: two spatial and one distance dimension along each line of sight towards the bulge. We present a 3D extinction map that covers the whole VVV area with a resolution of 6' x 6', using distance bins of 0.5 kpc. The high resolution and depth of the photometry allows us to derive extinction maps for a range of distances up to 10 kpc and up to 30 magnitudes of extinction in $A_{V}$. Integrated maps show the same dust features and consistent values as other 2D maps. We discuss the spatial distribution of dust features in the line of sight, which suggests that there is much material in front of the Galactic bar, specifically between 5-7 kpc. We compare our dust extinction map with high-resolution $\rm ^{12}CO$ maps towards the Galactic bulge, where we find a good correlation between $\rm ^{12}CO$ and $\rm A_{V}$. We determine the X factor by combining the CO map and our dust extinction map. Our derived average value is consistent with the canonical value of the Milky Way.

preprint2014arXiv

Metallicity and kinematics of the bar in-situ

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.

preprint2014arXiv

Milky Way demographics with the VVV Survey III. Evidence for a Great Dark Lane in the 157 Million Star Bulge Color-Magnitude Diagram

The new generation of IR surveys are revealing and quantifying Galactic features, providing an improved 3-D interpretation of our own Galaxy. We present an analysis of the global distribution of dust clouds in the bulge using the near-IR photometry of 157 million stars from the VVV Survey. We investigate the color magnitude diagram of the Milky Way bulge which shows a red giant clump of core He burning stars that is split in two color components, with a mean color difference of (Z-Ks)=0.55 magnitudes equivalent to A_V=2.0 magnitudes. We conclude that there is an optically thick dust lane at intermediate latitudes above and below the plane, that runs across several square degrees from l=-10 deg to l=+10 deg. We call this feature the "Great Dark Lane". Although its exact distance is uncertain, it is located in front of the bulge. The evidence for a large-scale great dark lane within the Galactic bulge is important in order to constrain models of the barred Milky Way bulge and to compare our galaxy with external barred galaxies, where these kinds of features are prominent. We discuss two other potential implications of the presence of the Great Dark Lane for microlensing and bulge stellar populations studies.

preprint2014arXiv

The Chemical Evolution of Fluorine in the Bulge - High-resolution K-band spectra of giants in three fields

Possible main formation sites of F in the Universe include AGB stars, the ν-process in Type II SNe, and/or W-R stars. The importance of the W-R stars has theoretically been questioned and they are probably not needed in the modelling of the chemical evolution of F in the solar neighborhood. It has, however, been suggested that W-R stars are indeed needed to explain the chemical evolution of F in the Bulge. The molecular spectral data of the often used HF-molecule has not been presented in a complete and consistent way and has recently been debated in the literature. In this article we determine the [F/O] vs. [O/H] trend in the Bulge to investigate the possible contribution from W-R stars. Additionally, we present here a HF line list for the K- and L-bands (including the often used 23358.33 Å line) and an accompanying partition function. The F abundances were determined using spectral fitting from hi-res NIR spectra of eight K giants recorded by the spectrograph CRIRES. We have also re-analyzed five previously published Bulge giants using our new HF molecular data. We find that the F-O abundance in the Bulge probably cannot be explained with chemical evolution models including only AGB-stars and the ν-process in SNe Type II, i.e. a significant amount of F production in W-R stars is likely needed to explain the F abundance in the Bulge. Concerning the HF line list, we find that a possible reason for the inconsistencies in the literature, with two different excitation energies being used, is two different definitions of the zero-point energy for the HF molecule and therefore also two accompanying different dissociation energies. Both line lists are correct, as long as the corresponding consistent partition function is used in the spectral synthesis. However, we suspect this has not been the case in several earlier works leading to F abundances 0.3 dex too high.

preprint2014arXiv

The Gaia-ESO Survey: metallicity and kinematic trends in the Milky Way bulge

(Abridged) We analyzed the stellar parameters and radial velocities of ~1200 stars in five bulge fields as determined from the Gaia-ESO survey data (iDR1). We use VISTA Variables in The Via Lactea (VVV) photometry to obtain reddening values by using a semi-empirical T_eff-color calibration. From a Gaussian decomposition of the metallicity distribution functions, we unveil a clear bimodality in all fields, with the relative size of components depending of the specific position on the sky. In agreement with some previous studies, we find a mild gradient along the minor axis (-0.05 dex/deg between b=-6 and b=-10) that arises from the varying proportion of metal-rich and metal-poor components. The number of metal-rich stars fades in favor of the metal-poor stars with increasing b. The K-magnitude distribution of the metal-rich population splits into two peaks for two of the analyzed fields that intersects the near and far branches of the X-shaped bulge structure. In addition, two lateral fields at (l,b)=(7,-9) and (l,b)=(-10,-8) present contrasting characteristics. In the former, the metallicity distribution is dominated by metal-rich stars, while in the latter it presents a mix of a metal-poor population and and a metal-intermediate one, of nearly equal sizes. Finally, we find systematic differences in the velocity dispersion between the metal-rich and the metal-poor components of each field. Our chemo-kinematical analysis is consistent with a varying field-to-field proportion of stars belonging to (i) a metal-rich boxy/peanut X-shaped component, with bar-like kinematics, and (ii) a metal-poor more extended rotating structure with a higher velocity dispersion that dominates far from the Galactic plane. These first GES data allow studying the detailed spatial dependence of the Galactic bulge populations, thanks to the analysis of individual fields with relatively high statistics.

preprint2013arXiv

High Resolution Spectroscopy of M subdwarfs

M subdwarfs are metal poor and cool stars. They are important probes of the old galactic populations. However, they remain elusive due to their low luminosity. Observational and modeling efforts are required to fully understand the physics and to investigate the effect of metallicity in their cool atmo- spheres. We perform a detail study of a sample of subdwarfs to determine their stellar parameters and constrain the atmosphere models. We present UVES/VLT high resolution spectra of 21 M subdwarfs. Our atlas covers the optical region from 6400°A up to the near infrared at 10000°A. We show spectral details of cool atmospheres at very high resolution (R = 40 000) and compare with synthetic spectra computed from the recent BT-Settl atmosphere models. Our comparison shows that molecular features (TiO, VO, CaH), and atomic features (Fe, Ti, Na, K) are well fitted by current models. We produce an effective temperature versus spectral type relation all over the subdwarf spectral sequence.

preprint2013arXiv

Target Selection for the Apache Point Observatory Galactic Evolution Experiment (APOGEE)

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) is a high-resolution infrared spectroscopic survey spanning all Galactic environments (i.e., bulge, disk, and halo), with the principal goal of constraining dynamical and chemical evolution models of the Milky Way. APOGEE takes advantage of the reduced effects of extinction at infrared wavelengths to observe the inner Galaxy and bulge at an unprecedented level of detail. The survey's broad spatial and wavelength coverage enables users of APOGEE data to address numerous Galactic structure and stellar populations issues. In this paper we describe the APOGEE targeting scheme and document its various target classes to provide the necessary background and reference information to analyze samples of APOGEE data with awareness of the imposed selection criteria and resulting sample properties. APOGEE's primary sample consists of ~100,000 red giant stars, selected to minimize observational biases in age and metallicity. We present the methodology and considerations that drive the selection of this sample and evaluate the accuracy, efficiency, and caveats of the selection and sampling algorithms. We also describe additional target classes that contribute to the APOGEE sample, including numerous ancillary science programs, and we outline the targeting data that will be included in the public data releases.

preprint2013arXiv

The Effective Temperature Scale of M dwarfs

Despite their large number in the Galaxy, M dwarfs remain elusive objects and the modeling of their photospheres has long remained a challenge (molecular opacities, dust cloud formation). Our objectives are to validate the BT-Settl model atmospheres, update the M dwarf T eff -spectral type relation, and find the atmospheric parameters of the stars in our sample. We compare two samples of optical spectra covering the whole M dwarf sequence with the most recent BT-Settl synthetic spectra and use a \c{hi}2 minimization technique to determine Teff . The first sample consists of 97 low-resolution spectra obtained with NTT at La Silla Observatory. The second sample contains 55 mid-resolution spectra obtained at the Siding Spring Observatory (SSO). The spectral typing is realized by comparison with already classified M dwarfs. We show that the BT-Settl synthetic spectra reproduce the slope of the spectral energy distribution and most of its features. Only the CaOH band at 5570Å and AlH and NaH hydrides in the blue part of the spectra are still missing in the models. The Teff-scale obtained with the higher resolved SSO 2.3 m spectra is consistent with that obtained with the NTT spectra. We compare our Teff scale with those of other authors and to published isochrones using the BT-Settl colors. We also present relations between effective temperature, spectral type and colors of the M dwarfs.

preprint2012arXiv

Reddening and metallicity maps of the Milky Way bulge from VVV and 2MASS II. The complete high resolution extinction map and implications for Bulge studies

We use the Vista Variables in the Via Lactea (VVV) ESO public survey data to measure extinction values in the complete area of the Galactic bulge covered by the survey at high resolution. We derive reddening values using the method described in Paper I. This is based on measuring the mean (J-Ks) color of red clump giants in small subfields of 2' to 6' in the following bulge area: -10.3<b<+5.1 and -10<l<+10.4. To determine the reddening values E(J-Ks) for each region, we measure the RC color and compare it to the (J-Ks) color of RC stars measured in Baade's window, for which we adopt E(B-V)=0.55. This allows us to construct a reddening map sensitive to small scale variations minimizing the problems arising from differential extinction. The significant reddening variations are clearly observed on spatial scales as small as 2'. We find a good agreement between our extinction measurements and Schlegel maps in the outer bulge, but, as already stated in the literature the Schlegel maps are not reliable for regions within |b| < 6. In the inner regions we compare our results with maps derived from DENIS and Spitzer surveys. While we find good agreement with other studies in the corresponding overlapping regions, our extinction map has better quality due to both higher resolution and a more complete spatial coverage in the Bulge. We investigate the importance of differential reddening and demonstrate the need for high resolution extinction maps for detailed studies of Bulge stellar populations and structure. The extinction variations on scales of up to 2'-6', must be taken into account when analysing the stellar populations of the Bulge.

preprint2012arXiv

The very low mass multiple system LHS\,1070 -- a testbed for model atmospheres for the lower end of the main sequence

LHS1070 is a nearby multiple system of low mass stars. It is an important source of information for probing the low mass end of the main sequence, down to the hydrogen-burning limit. The primary of the system is a mid-M dwarf and two components are late-M to early L dwarfs, at the star-brown dwarf transition. Hence LHS1070 is a valuable object to understand the onset of dust formation in cool stellar atmospheres.This work aims at determining the fundamental stellar parameters of LHS1070 and to test recent model atmospheres: BT-Dusty,BT-Settl, DRIFT, and MARCS models.Unlike in previous studies, we have performed a chi^2-minimization comparing well calibrated optical and infrared spectra with recent cool star synthetic spectra leading to the determination of the physical stellar parameters Teff, radius, and log g for each of the three components of LHS1070.

preprint2012arXiv

Three dimensional interstellar extinction map towards the Galactic Bulge

Studies of the properties of the inner Galactic Bulge depend strongly on the assumptions about the interstellar extinction. Most of the extinction maps available in the literature lack the information about the distance. We combine the observations with the Besancon model of the Galaxy to investigate the variations of extinction along different lines of sight towards the inner Galactic bulge as a function of distance. In addition we study the variations in the extinction law in the Bulge. We construct color-magnitude diagrams with the following sets of colors: H-Ks and J-Ks from the VVV catalogue as well as Ks-[3.6], Ks-[4.5], Ks-[5.8] and Ks-[8.0] from GLIMPSE-II catalogue matched with 2MASS. Using the newly derived temperature-color relation for M giants that match better the observed color-magnitude diagrams we then use the distance-color relations to derive the extinction as a function of distance. The observed colors are shifted to match the intrinsic colors in the Besançon model, as a function of distance, iteratively thereby creating an extinction map with three dimensions: two spatial and one distance dimension along each line of sight towards the bulge. Colour excess maps are presented at a resolution of 15' x 15' for 6 different combinations of colors, in distance bins of 1 kpc. The high resolution and depth of the photometry allows to derive extinction maps to 10 kpc distance and up to 35 magnitudes of extinction in Av (3.5 mag in Aks). Integrated maps show the same dust features and consistent values with the other 2D maps. Starting from the color excess between the observations and the model we investigate the extinction law in near-infrared and its variation along different lines of sight.

preprint2011arXiv

The Effective temperature scale of M dwarfs from spectral synthesis

We present a comparison of low-resolution spectra of 60 stars covering the whole M-dwarf sequence. Using the most recent PHOENIX BT-Settl stellar model atmospheres (see paper by F. Allard, in this book) we do a first quantitative compari- son to our observed spectra in the wavelength range 550-950 nm. We perform a first confrontation between models and observations and we assign an effective tempera- tures to the observed M-dwarfs. Teff-spectral type relations are then compared with the published ones. This comparison also aims at improving the models' opacities.

preprint2011arXiv

The metallicity distribution of bulge clump giants in Baade's Window

We seek to constrain the formation of the Galactic bulge by means of analysing the detailed chemical composition of a large sample of red clump stars in Baade's window. We measure [Fe/H] in a sample of 219 bulge red clump stars from R=20000 resolution spectra obtained with FLAMES/GIRAFFE at the VLT, using an automatic procedure, differentially to the metal-rich local reference star muLeo. For a subsample of 162 stars, we also derive [Mg/H] from spectral synthesis around the MgI triplet at 6319A. The Fe and Mg metallicity distributions are both asymmetric, with median values of +0.16 and +0.21 respectively. The iron distribution is clearly bimodal, as revealed both by a deconvolution (from observational errors) and a Gaussian decomposition. The decomposition of the observed Fe and Mg metallicity distributions into Gaussian components yields two populations of equal sizes (50% each): a metal-poor component centred around [Fe/H]=-0.30 and [Mg/H]=-0.06 with a large dispersion and a narrow metal-rich component centred around [Fe/H]=+0.32 and [Mg/H]=+0.35. The metal poor component shows high [Mg/Fe] ratios (around 0.3) whereas stars in the metal rich component are found to have near solar ratios. Babusiaux et al. (2010) also find kinematical differences between the two components: the metal poor component shows kinematics compatible with an old spheroid whereas the metal rich component is consistent with a population supporting a bar. In view of their chemical and kinematical properties, we suggest different formation scenarios for the two populations: a rapid formation timescale as an old spheroid for the metal poor component (old bulge) and for the metal rich component, a formation over a longer time scale driven by the evolution of the bar (pseudo-bulge).

preprint2010arXiv

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

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. 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, b=-4) and two other fields along the bulge minor axis (l=0, b=-6 and b=-12). We combine the radial velocity and the [Fe/H] measurements obtained with FLAMES/GIRAFFE at the VLT with a spectral resolution of R=20000, plus for the Baade's Window field the OGLE-II proper motions, and compare these with published N-body simulations of the Galactic bulge. We confirm the presence of two distinct populations in Baade's Window found in Hill et al. 2010: the metal-rich population presents bar-like kinematics while the metal-poor population shows kinematics corresponding to an old spheroid or a thick disc one. In this context the metallicity gradient along the bulge minor axis observed by Zoccali et al. (2008), 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. We show evidences that the two main scenarios for the bulge formation co-exist within the Milky Way bulge.

preprint2008arXiv

Photometric redshifts for the CFHTLS T0004 Deep and Wide fields

We compute photometric redshifts based on the template-fitting method in the fourth public release of the Canada-France-Hawaii Telescope Legacy Survey. This unique multi-colour catalogue comprises u*,g',r',i',z' photometry in four deep fields of 1 deg2 each and 35 deg2 distributed over three Wide fields. Our photometric redshifts are calibrated with and compared to 16,983 high-quality spectroscopic redshifts from several surveys. We find a dispersion of 0.028 and an outlier rate of 3.5% in the Deep field at i'AB < 24 and a dispersion of 0.036 and an outlier rate of 2.8% in the Wide field at i'AB < 22.5. Beyond i'AB = 22.5 in the Wide field the number of outliers rises from 5% to 10% at i'AB<23 and i'AB<24 respectively. For the Wide sample, we find the systematic redshift bias keeps below 1% to i'AB < 22.5, whereas we find no significant bias in the Deep field. We investigated the effect of tile-to-tile photometric variations and demonstrate that the accuracy of our photometric redshifts is reduced by at most 21%. We separate stars from galaxies using both the size and colour information, reducing the contamination by stars in our catalogues from 50% to 8% at i'AB < 22.5 in fields with the highest stellar density while keeping a complete galaxy sample. Our CFHTLS T0004 photometric redshifts are distributed to the community. Our release include 592,891 (i'AB < 22.5) and 244,701 (i'AB < 24) reliable galaxy photometric redshifts in the Wide and Deep fields, respectively.

preprint2006arXiv

Modelling the Galactic Interstellar Extinction Distribution in Three Dimensions

The Two Micron All Sky Survey, along with the Stellar Population Synthesis Model of the Galaxy, developed in Besancon, is used to calculate the extinction distribution along different lines of sight. By combining many lines of sight, the large scale distribution of interstellar material can be deduced. The Galaxy model is used to provide the intrinsic colour of stars and their probable distances, so that the near infrared colour excess, and hence the extinction, may be calculated and its distance evaluated. Such a technique is dependent on the model used, however we are able to show that moderate changes in the model parameters result in insignificant changes in the predicted extinction. This technique has now been applied to over 64000 lines of sight, each separated by 15 arcmin, in the inner Galaxy (|l|<=100 deg, |b|<=10 deg). Using our extinction map, we have derived the main characteristics of the large scale structure of the dust distribution: scale height and warp of the ISM disc as well as the angle of the dust in the Galactic Bar. This resulting extinction map will be useful for studies of the inner Galaxy and its stellar populations.

preprint2006arXiv

Recent star formation in the inner Galactic Bulge seen by ISOGAL. I - Classification of bright mid-IR sources in a test field

Context: The stellar populations in the central region of the Galaxy are poorly known because of the high visual extinction and very great source density in this direction. Aims: To use recent infrared surveys for studying the dusty stellar objects in this region. Methods: We analyse the content of a 20x20 arcmin^2 field centred at (l,b)=(-0.27,-0.06) observed at 7 and 15 microns as part of the ISOGAL survey. These ISO observations are more than an order of magnitude better in sensitivity and spatial resolution than the IRAS observations. The sources are cross-associated with other catalogues to identify various types of objects. We then derive criteria to distinguish young objects from post-main sequence stars. Results: We find that a sample of about 50 young stellar objects and ultra-compact HII regions emerges, out of a population of evolved AGB stars. We demonstrate that the sources colours and spatial extents, as they appear in the ISOGAL catalogue, possibly complemented with MSX photometry at 21 microns, can be used to determine whether the ISOGAL sources brighter than 300 mJy at 15 microns (or [15] < 4.5 mag) are young objects or late-type evolved stars.

preprint2005arXiv

Stellar populations in the CFHTLS. I. New constraints on the IMF at low mass

We present a stellar populations analysis of the first release of the CFHTLS (Canada-France-Hawai Telescope Legacy Survey) data. A detailed comparison between the Besançon model of the Galaxy and the first data release of the CFHTLS-Deep survey is performed by implementing the MEGACAM photometric system in this model using stellar atmosphere model libraries. The reliability of the theoretical libraries to reproduce the observed colours in the MEGACAM system is investigated. The locations of various stellar species like subdwarfs, white dwarfs, late-type and brown dwarfs, binary systems are identified. The contamination of the stellar sample by quasars and compact galaxies is quantified using spectroscopic data from the VIMOS-VLT Deep Survey (VVDS) as a function of $i'$ magnitude and $r'-i'$ colour. A comparison between simulated counts using the standard IMF at low masses show that the number of very low mass dwarfs may have been underestimated in previous studies. These observations favour a power law IMF following $d(n)/dm \propto m^{-α}$ with $α=2.5$ for $m < $ 0.25 \Msun or $α=3.0$ for $m < $ 0.2 \Msun for single stars. The resulting LF is in agreement with the local LF as measured from the 5 or 25 pc samples. It is in strong disagreement with the Zheng et al (2001) LF measured from deep HST data. We show that this discrepancy can be understood as an indication of a different IMF at low masses at early epochs of the Galaxy compared to the local thin disc IMF.

preprint1998arXiv

Field #3 of the Palomar-Groningen Survey II. Near-infrared photometry of semiregular variables

Near-infrared photometry (JHKL'M) was obtained for 78 semiregular variables (SRVs) in field #3 of the Palomar-Groningen survey (PG3, l=0, b=-10). Together with a sample of Miras in this field a comparison is made with a sample of field SRVs and Miras. The PG3 SRVs form a sequence (period-luminosity & period-colour) with the PG3 Miras, in which the SRVs are the short period extension to the Miras. The field and PG3 Miras follow the same P/(J--K)o relation, while this is not the case for the field and PG3 SRVs. Both the PG3 SRVs and Miras follow the SgrI period-luminosity relation adopted from Glass et al. (1995, MNRAS 273, 383). They are likely pulsating in the fundamental mode and have metallicities spanning the range from intermediate to approximately solar.