The investigation of the metal-poor tail in the Galactic bulge provides unique information on the early Milky Way assembly and evolution. A chemo-dynamical analysis of 17 very metal-poor stars (VMP, [Fe/H < – 2.0]) selected from the Pristine Inner Galaxy Survey was carried out based on Gemini/GRACES spectra. The chemistry suggests that the majority of our stars are very similar to metal-poor stars in the Galactic halo. Orbits calculated from Gaia EDR3 imply these stars are brought into the bulge during the earliest Galactic assembly. Most of our stars have large [Na,Ca/Mg] abundances, and thus show little evidence of enrichment by pair-instability supernovae. Two of our stars (P171457, P184700) have chemical abundances compatible with second-generation globular cluster stars, suggestive of the presence of ancient and now dissolved globular clusters in the inner Galaxy. One of them (P171457) is extremely metal-poor ([Fe/H < – 3.0]) and well below the metallicity floor of globular clusters, which supports the growing evidence for the existence of lower-metallicity globular clusters in the early Universe. A third star (P180956, [Fe/H]∼ – 2) has low [Na,Ca/Mg] and very low [Ba/Fe] for its metallicity, which are consistent with formation in a system polluted by only one or a few low-mass supernovae. Interestingly, its orbit is confined to the Galactic plane, like other very metal-poor stars found in the literature, which have been associated with the earliest building blocks of the Milky Way.

The chemo-dynamics of the stellar populations in the Galactic Bulge inform and constrain the Milky Way’s formation and evolution. The metal-poor population is particularly important in light of cosmological simulations, which predict that some of the oldest stars in the Galaxy now reside in its center. The metal-poor bulge appears to consist of multiple stellar populations that require chemo-dynamical analyses to disentangle. In this paper, I describe the detailed chemo-dynamical study of the metal-poor stars in the inner Galaxy, named The COMBS Survey which uses VLT/FLAMES spectra of 350 metal-poor stars. I discuss the results and the implications for early Milky Way formation and chemical evolution. In addition, I preview results from an ongoing survey of carbon-enhanced metal-poor stars, which are thought to be solely enriched by the first generation of stars.

The Bulge Asymmetries and Dynamical Evolution (BAaDE) survey aims to explore the complex structure of the inner Galaxy and Galactic Bulge, by using the 43 GHz receivers at the Karl G. Jansky Very Large Array (VLA) and the 86 GHz receivers at the Atacama Large Millimeter/submillimeter Array (ALMA) to observe SiO maser lines in red giant stars. The goal is to construct a sample of stellar point-mass probes that can be used to test models of the gravitational potential, and the final sample is expected to provide at least 20,000 line-of-sight velocities and positions. A possible bias between the VLA and the ALMA SiO maser lines is explored, and the 86 GHz SiO line-peak velocities agree using either of the four sampled lines. Additionally, the SiO maser velocities agree with the OH maser derived velocities.

In the efforts to map the Milky Way structure, the central regions have remained very difficult to probe. The VISTA Variables in the Vía Láctea Survey (VVV) is a near-IR variability Survey that scans 560 sq.deg. across the Milky Way bulge and an adjacent section of the southern mid-plane. The main goal of the VVV Survey is to build a 3D map of the structure of the inner Galaxy and characterize its stellar populations. This survey has discovered different kinds of objects, such as globular clusters, Microlensing events, RR Lyrae stars, Cepheids, WITs, among others. The extension of the Survey (VVVX) is observing until 2020, tripling the areal coverage, and complementing the variability studies done by the VVV Survey.

We have derived absolute proper motions of stars in the Galactic bulge region combining the VVV InfraRed Astrometric Catalogue (VIRAC) and Gaia. We use the proper motions to study the kinematic structure of the bulge both integrated along the line-of-sight and in magnitude intervals using red clump stars as standard candles. In parallel we compare to a made-to-measure barred dynamical model, folding in the VIRAC selection function, to understand and interpret the structures that we observe. The barred dynamical model, which contains a boxy/peanut bulge, and has a pattern speed of 37.5 kms−1 kpc−1, is able to reproduce all structures impressively well.

In order to study the most reddened areas of the Milky Way we used near-IR data from the VVV Survey. For the first time, the VISTA telescope allows us to observe the mid-plane through the Galactic bulge and study the disk in the other side of the Milky Way. Motivated by the detection of hundreds of microlensing events in the inner regions of the Galaxy, we propose three new configurations of microlensing events, placing the sources in the far-disk and the lenses in the far-disk/bulge/near-disk. These new configurations will change the usual way to interpret the timescale distributions due to the different populations along the line of sight, that exhibit varied transverse velocities and relative distances.

The Bulge Asymmetries and Dynamical Evolution (BAaDE) survey aims to use circumstellar SiO maser line-of-sight velocities as probes for the Galactic gravitational potential and dynamical structure. The SiO masers are detected at a high rate in specific color-selected MSX infrared sources. Furthermore, the SiO maser properties and line ratios, in combination with infrared spectral energy distributions and location in the Galaxy, will statistically yield detailed information on population and evolution of low- to intermediate-mass evolved stars in the Galaxy.

Recent analyses of Lee et al. (2018, 2019) have confirmed that Galactic bulge consists of stellar populations originated from Milky Way globular clusters (MWGCs). Motivated by this, here we present the evolutionary population synthesis (EPS) for the Galactic bulge and early-type galaxies (ETGs) with the realistic treatment of individual variations in light elements observed in the MWGCs. We have utilized our model with GC-origin populations to explain the CN spread observed in ETGs, and the results show remarkable matches with the observations. We further employ our model to estimate the age of ETGs, which are considered as good analogs for the MW bulge. We find that, without the effect of our new treatments, EPS models will almost always underestimate the true age of ETGs. Our analysis indicates that the EPS with GC-origin populations is an essential constraint in determining the ETG formation epoch and is closely related to understanding the evolution of the Universe.

We report on the successful search for CO (2-1) and (3-2) emission associated with OH/IR stars in the Galactic Bulge. We observed a sample of eight extremely red AGB stars with the APEX telescope and detected seven. The sources were selected at sufficient high Galactic latitude to avoid interference by interstellar CO, which hampered previous studies of inner galaxy stars. We also collected photometric data and Spitzer IRS spectroscopy to construct the SEDs, which were analysed through radiative transfer modelling. We derived variability periods of our stars from the VVV and WISE surveys. Through dynamical modelling we then retrieve the total mass loss rates (MLR) and the gas-to-dust ratios. The luminosities range between approximately 4,000 and 5,500 L⊙ and periods are below 700 days. The total MLR ranges between 10−5 and 10−4 M⊙ yr−1. The results are presented in Blommaert et al. 2018 and summarized below.

We examine the metallicity trends in the Milky Way (MW) bulge – using APOGEE DR13 data – and explore their origin by comparing two N-body models of isolated galaxies which develop a bar and a boxy/peanut (b/p) bulge. Both models have been proposed as scenarios for reconciling a disc origin of the MW bulge with a negative vertical metallicity gradient. The first is a superposition of co-spatial disc populations, different scaleheights and metallicities (with flat gradients) where the thick, metal-poor populations contribute significantly to the stellar mass budget in the inner galaxy. The second model is a single disc with an initial steep radial metallicity gradient which gets mapped by the bar into the b/p bulge in such a way that the vertical metallicity gradient of the MW bulge is reproduced – as shown already in previous works in the literature. As we show here, the latter model does not reproduce the positive longitudinal metallicity gradient of the inner disc, nor the metal-poor innermost regions seen in the data. The model with co-spatial thin and thick disc populations reproduces all the aforementioned trends. We therefore see that it is possible to reconcile a (primarily) disc origin for the MW bulge with the observed trends in metallicity by mapping the inner thin and thick discs of the MW into a b/p.

Radio astrometric campaigns using VLBI have provided distances and proper motions for masers associated with young massive stars (BeSSeL survey). The ongoing BAaDE project plans to obtain astrometric information of SiO maser stars located in the inner Galaxy. These stars are associated with evolved, mass-losing stars. By overlapping optical (Gaia), infrared (2MASS, MSX and WISE) and radio (BAaDE) sources, we expect to obtain important clues on the intrinsic properties and population distribution of late-type stars. Moreover, a comparison of the Galactic parameters obtained with Gaia and VLBI can be done using radio observations on different targets: young massive stars (BeSSeL) and evolved stars (BAaDE).

A major avenue in the study of the Galaxy is the investigation of stellar populations and Galactic chemical evolution by stellar spectroscopy. Due to the dust obscuration, stars in the centre of the Galaxy can only be observed in the near-IR wavelength region. However, existing line lists in this wavelength region are demonstratively not of good enough quality for use in stellar spectroscopy. In response to this, we have developed an empirical astrophysical line list in the K-band based on modelling against the Sun and testing against Arcturus. Of ca. 700 identified interesting lines about 570 lines have been assigned empirically determined values.

Using the VLBA, the BeSSeL survey has provided distances and proper motions of young massive stars, allowing an accurate measure of the Galactic spiral structure. By the same technique, we are planning to map the inner Galaxy using positions and velocities of evolved stars (provided by the BAaDE survey). These radio astrometric measurements (BeSSeL and BAaDE) will be complementary to Gaia results and the overlap will provide important clues on the intrinsic properties and population distribution of the stars in the bulge.

Owing to their extreme crowding and high and variable extinction, stars in the Galactic Bulge, within ±2° of the Galactic plane, and especially those in the Nuclear Star Cluster, have only rarely been targeted for an analyses of their detailed abundances. There is also some disagreement about the high end of the abundance scale for these stars. It is now possible to obtain high dispersion, high S/N spectra in the infrared K band (~2.0 − 2.4 µm) for these giants; we report our progress at Keck and VLT in using these spectra to infer the composition of this stellar population.

I present examples of how chemo-dynamical N-body simulations can help understanding the structure and evolution of the inner Galaxy. Such simulations reproduce the observed links between kinematics, morphology and chemistry in the bar/bulge region and explain them by the self-consistent cohabitation of a number of components. Galactic archaeology, applied to simulation snapshots, explains the sequence in which the stars of the various components were formed. The thick disc stars form earlier than those of the thin disc and in a much shorter time scale. The bar in the thick disc is horizontally thicker than that of the thin disc and has a different vertical morphology. The Galaxy’s inner disc scalelength is much smaller than what is expected from nearby galaxies of similar stellar mass.

Milky Way globular clusters are excellent laboratories for stellar population detailed analysis that can be applied to extragalactic environments with the advent of the 40m-class telescopes like the ELT. The globular cluster population traces the early evolution of the Milky Way which is the field of Galactic archaeology. We present our GlObular clusTer Homogeneous Abundance Measurement (GOTHAM) survey. We derived radial velocities, Teff, log(g), [Fe/H], [Mg/Fe] for red giant stars in one third of all Galactic globular clusters that represent well the Milky Way globular cluster system in terms of metallicity, mass, reddening, and distance. Our method is based on low-resolution spectroscopy and is intrinsically reddening free and efficient even for faint stars. Our [Fe/H] determinations agree with high-resolution results to within 0.08 dex. The GOTHAM survey provides a new metallicity scale for Galactic globular clusters with a significant update of metallicities higher than [Fe/H] > -0.7. We show that the trend of [Mg/Fe] with metallicity is not constant as previously found, because now we have more metal-rich clusters. Moreover, peculiar clusters whose [Mg/Fe] does not match Galactic stars for a given metallicity are discussed. We also measured the CaII triplet index for all stars and we show that the different chemical evolution of Milky Way open clusters, field stars, and globular clusters implies different calibrations of calcium triplet to metallicity.

Recently the search for the oldest stars have started to focus on the Bulge region. The Galactic bulge hosts extremely old stars, with ages compatible with the ages of the oldest halo stars. The data coming from these recent observations present new chemical signatures and therefore provide complementary constraints to those already found in the halo. So, the study of the oldest bulge stars can improve dramatically the constraints on the nature of first stars and how they polluted the pristine ISM of our Galaxy. We present our first results regarding the light elements (CNO) and the neutron capture elements. Our findings in the oldest bulge stars support the scenario where the first stellar generations have been fast rotators.

Surveys of pulsating stars in the inner Galaxy have been very limited, but recent large-scale surveys are rapidly bringing us new samples of various kinds of variable stars and new insights into stellar populations therein. Because of the severe interstellar extinction along the Galactic disc, the stellar populations in the inner Galaxy are more easily observed in the infrared, but even in the infrared the interstellar extinction may cause a serious problem in revealing their accurate characteristics. Here we review recent discoveries of Cepheids and Miras, two kinds of luminous pulsating stars with period–luminosity relation, in the inner Galaxy.

4MOST is a new wide-field, high-multiplex spectroscopic survey facility for the VISTA telescope of ESO. Starting in 2022, 4MOST will deploy 2400 fibres in a 4.1 square degree field-of-view using a positioner based on the tilting spine principle. In this contribution we give an outline of the major science goals we wish to achieve with 4MOST in the area of Galactic Archeology. The 4MOST Galactic Archeology surveys have been designed to address long-standing and far-reaching problems in Galactic science. They are focused on four major themes: 1) Near-field cosmology tests, 2) Chemo-dynamical characterisation of the major Milky Way stellar components, 3) The Galactic Halo and beyond, and 4) Discovery and characterisation of extremely metal-poor stars. In addition to a top-level description of the Galactic surveys we provide information about how the community will be able to join 4MOST via a call for Public Spectroscopic Surveys that ESO will launch.

Some highlights are given of the IAU Symposium 334, Rediscovering our Galaxy, held in Potsdam, in July 2017: from the first stars fossil records found in the halo, the carbon-enhanced metal poor CEMP-no, to the cosmological simulations presenting possible scenarios for the Milky Way formation, passing through the chemo-dynamical models of the various components, thin and thick disks, box/peanut bulge, halo, etc. The domain is experiencing (or will be in the near future) huge improvements with precise and accurate stellar ages, provided by astero-seismology, precise stellar distances and kinematics (parallaxes and proper motions from GAIA), and the big data resulting from large surveys are treated with deep learning algorithms.