RR Lyrae Variables as Tracers of the Galactic Bulge Kinematic Structure

: RR Lyrae stars are recognized as some of the oldest stars in the Universe. In addition, they are some of the few old celestial objects for which distances can be reliably inferred. As such, these stars are excellent tracers of the oldest structures that exist in the inner Galaxy. Although the inner Galaxy is where the oldest structures in the Milky Way are thought to be hidden, it is also a region notoriously difﬁcult to study due to high extinction and crowding. Here, I will summarize how RR Lyrae stars have been used to obtain a more complete picture of the inner Galaxy. In particular, recently, a large sample of RR Lyrae star motions through space have been obtained and compared to younger, more metal-rich stars in the bulge/bar. It is seen that the inner Galaxy RR Lyrae star kinematics are complicated by a mix of a variety of Galactic components. After isolating only those RR Lyrae stars that are conﬁned to the bulge, a subsample of these stars have slower rotation and are less barred than the dominant bar/bulge. Curiously, there is no discernible metallicity [ Fe/H ] difference between these two subsamples. Old, metal-poor stars in the inner Galaxy need to be properly accounted for when discussing processes that gave rise to the formation of the inner Galaxy and the Galactic bar/bulge.


Introduction
Variable stars known as RR Lyrae stars (RRLs, hereafter) are a stellar tracer that has been trusted for decades to help piece together the size and structure of the Milky Way (MW) Galaxy [1][2][3][4][5].RRLs are variables that fall under radially pulsating stars with periods ranging from ∼0.2 to 1.0 days.These stars begin their lives on the Hertzsprung-Russell (HR) diagram with masses and radii similar to that of our Sun [6].They have since evolved off the main-sequence, and after helium flash, reside now on the horizontal branch, undergoing core helium burning.In particular, they fall within the instability strip that crosses the horizontal branch (see Figure 6 in Catelan 2007), [6].Therefore, they have a relatively low mass (∼0.6-0.8M ), temperatures between 6000 K-7250 K, and radii typically between 4-6 R [7,8].Due to this stellar evolution path, the existence of RRLs in a stellar population signifies a system older than about 10 Gyr [9].RRLs therefore tell us about conditions in the Galaxy at the time they were formed [10].
The oscillatory phenomenon seen in RRLs is attributed to the propagation of sound waves throughout their stellar interior [11].One of the first papers that systematically categorized variable stars was from Bailey (1902), who studied stars that change their brightness in the Globular cluster ω Centauri (NGC 5139) [12].Bailey divided the RRLs into three subclasses, a, b, and c, which were distinguished by the shape of their light curve as well as their pulsation amplitudes and periods.It is now known that RRLs of "Bailey" types of a or b are RRLs that pulsate in the fundamental mode with periods typically longer than ∼0.45 days, and RRLs of "Bailey" type of c are RRLs that pulsate in the first overtone with periods typically shorter than ∼0.45 days [13].Since then, RRLs that pulsate in a combination of the first and fundamental overtone-RRd-have been reported [14][15][16], as well as those that pulsate in the second overtone-RRe [17,18].
The physics of these stellar pulsations in RRLs is an active field.However, RRLs are also important standard candles for measuring the distances to very old stellar populations, such as globular clusters, as was realized already by Shapley (1918) [19].In fact, RRLs were some of the first celestial objects that provided a distance to the center of the Galaxy [2,20].In these studies, RRL distances were obtained using V-magnitudes of RR Lyrae stars.Indeed, until relatively recently, the absolute magnitude of an RRL was investigated mostly using the optical passbands, particularly the V-band.
RRLs have long been known to exist in the bulge [20].Microlensing surveys such as the MACHO project [21] and Optical Gravitational Lensing Experiment (OGLE) [22] have uncovered thousands of RRLs in the central regions of the Galaxy.RRLs can be easily identified on the basis of their colors (bluer than typical bulge stars) and distinctive sawtooth-shaped light curves (with periods of about 1/2 day), thus enabling the selection of pure bulge samples.As they are well-understood standard candles (with a precision of <10%), their presence within the bulge can be securely established.RRLs are the most significant and identifiable population of metal-poor stars in the bulge.
The oldest stellar populations in the bulge can reveal the early history of formation and assembly for our Galaxy and can constrain models for the formation and assembly of spiral galaxies in a cosmological context.Here, we discuss how RRLs are contributing to the knowledge of the formation of the bulge and inner Galaxy.In Section 2, both spectroscopic and photometric metallicities of RRLs in the inner Galaxy are discussed.The RRL spatial distribution of the RRLs is reviewed in Section 3 particularly with respect to different RRL surveys.In Section 4 the motions of these stars are reviewed, as well as how velocities have helped constrain scenarios explaining the formation of the inner Galaxy.

Metallicities of Inner Galaxy RR Lyrae Stars
A fundamental characteristic of any stellar system is its metallicity distribution function (MDF).The [Fe/H] distribution is related to the stellar system's star formation and also provides constraints on models of chemical evolution.The majority of the bulge stars have [Fe/H] metalllicites that are metal-rich (approximately solar), and their α-elements are enhanced [23], suggesting rapid early star formation.The bulge seems to show a radial abundance gradient in the sense that more metal-rich stars dominate closer to the Galactic plane and more sub-solar metallicity stars increase further from the Galactic plane [24,25].For RR Lyrae stars in particular, [Fe/H] abundances are also desirable because the absolute magnitude of RRLs are dependent on metallicity, as described in Section 3.1 below.

Spectroscopic Metallicities of Inner Galaxy RR Lyrae Stars
Unfortunately, there are very few spectroscopic abundances of bulge RRLs.This is partly due to the complexity of measuring metallicities for RRLs at the distance of the bulge.The commonly used lines in the blue region of the spectrum (Ca H&K, the Mg complex) are difficult to use, firstly because spectrographs are relatively inefficient in the blue, and secondarily because these moderately hot stars are faint at these wavelengths (I ∼16).Furthermore, in order to obtain accurate metallicities, spectra can only be taken when the atmospheres are most stable, at minimum light, as the expansion and contraction during the pulsations create shocks and alter abundances and chemistry in difficult to predict ways.
For example, Pancino et al. analyzed a sample of 18 RRLs in high-resolution (R∼30,000), high signal-to-noise (SNR > 30) and measured the atmospheric parameter variations due to RR Lyrae pulsations [26].They found that pulsations can reach ∆T eff ≈ 1500 Kelvin and ∆ log g ≈ 0.70 dex.To overcome this, the optimal region for spectroscopic analysis is near minimum light, with a duration of approximately 10% of the pulsational period [26].
A proper survey of metal abundances for bulge RRL has not been carried out since the early work by Walker & Terndrup (1991) [27].They analyzed RRLs in Baade's window in the 3750-5300 Å range, comparing the hydrogen spectral type, using the equivalent widths of the H γ and H δ lines, to the Ca II K line spectral type, using the equivalent width of the Ca II K line.From 59 RRLs, they found that the bulge RRL has a sharply peaked [Fe/H] distribution with a dispersion of only 0.16 dex, far smaller than other bulge stellar populations [27].
More recently, Savino et al. (2020) used the calcium II triplet (CaT) region, covering a wavelength range of about 8300-8800 Å, with a resolution of R∼10,000, to determine [Fe/H] metallicities [28].Because the SNR of the spectra was ∼20, their individual [Fe/H] metallicities are not as reliable as looking at the mean of the MDF.Their bulge RRL MDF with ∼900 RRL measurements covers a wide range of metallicities, with a single prominent peak at [Fe/H] = −1.39.This RRL metallicity distribution is shown in Figure 1, and compared to the [Fe/H] distribution of ∼2 million bulge red clump stars from the Blanco DECam Bulge Survey (BDBS) [25].Both BDBS and the Savino et al. (2020) RRL sample focus on stars located in the Southern bulge, and cover a similar range in Galactic latitude.However, BDBS covers a larger, more homogenous area of the sky, and extends to a significantly larger range in Galactic longitude.The uncertainty distribution of the RRL metallicity measurements almost entirely accounts for the width of the observed MDF peak, which suggests that the true underlying MDF has a very narrow, dominant, component, superimposed to wide metal-rich and metal-poor tails.
The bulge RRLs therefore appear to have a metallicity distribution function centered around [Fe/H] ∼ −1.0-1.4 [27,28].This is at least 10× lower in metallicity than the dominant bulge population [29,30].Their lower MDF indicates that they either trace a more ancient population as compared to the bulk of bulge stars or that their progenitors come preferentially from the more metal-poor stars (more metal-rich stars are less likely to evolve on to the instability strip).This is important because separating bulge stars into older/younger stars, or more metal-rich/metal-poor stars is challenging; the MDF of the bulge RRLs clearly indicates that targeting RRLs allows one to probe the older and metal-poor stellar population in the bulge.2020) [25].The distribution of RRLs is considerably more metal-poor than the dominant bulge population.

Photometric Metallicities of Inner Galaxy RR Lyrae Stars
Besides deriving [Fe/H] metallicities from spectra, there is also the possibility to estimate [Fe/H] from the photometric light curves of RRLs.For photometric metallicity calibrations of RRLs, [Fe/H] is estimated from the pulsation properties of the star (e.g., period and/or amplitude) and their Fourier coefficients.Five notable photometric metallicity calibrations of RRLs include those from Jurcsik and Kovacs (1996), Sandage (2004) Each of these calibrations is widely used in the literature, and each is calibrated to be consistent with the scale of Jurcsik & Kovacs (1996) [31,35], which is itself compiled from a heterogeneous list of 31 different references, and includes only two stars with [Fe/H] ≤ −2.0.It is worth noting that none of the above calibrations accounts for the effects of variable [α/Fe] abundances.This may be a problem when using photometric calibrations to determine [Fe/H] metallicities for RRLs in the bulge, since, for bulge RRLs, a shift of 0.4 dex in [α/Fe] is equivalent to a 0.3 dex shift in [Fe/H] [36].Even in the solar neighborhood, Prudil et al. (2020) showed that a number of local RRLs originated from the Gaia-Enceladus system [37][38][39], and these have [α/Fe] values that are ≈0.20 dex lower than the RRL in the halo stars [40].For photometric metallicities accurate at the ∼0.1-0.2 dex level, accounting for [α/Fe] is likely important.
Recently, Dekany et al. (2021) revisited the photometric metallicity prediction in the Cousins I passband, which is used in the OGLE survey [41].They used a calibration sample of 80 RRab stars ranging from −2.6 to +0.2 dex, 60 of which had multiple independent metallicity determinations.They found that all previously used formulae overestimate the modes of the MDFs of the Galactic bulge, disk, and the Magellanic Clouds by ∼0.35-0.4dex and cause distortions in their shapes.As such, at this time, it is not completely clear how well photometric metallicities predict the underlying MDF of the bulge RRLs-if the scatter in the [Fe/H] metallicity distribution is dominated by intrinsic scatter or by difficulties associated in the calibration and methodology in estimating photometric metallicities.
The MACHO and OGLE surveys have used the extensive photometric RRL measurements in their surveys to investigate the bulge MDF using photometric metallicities.Results using 2435 MACHO RRLs [42], using 16,836 OGLE-III RRL [43], and using 27,258 OGLE-IV RRL [44] all suggest an inner Galaxy RRL bulge MDF with <[Fe/H]> = −1.25 (on the Zinn and West 1984 scale [45]), or <[Fe/H]> = −1.0(on the Jurcsik 1995 scale [35]), with a broad metallicity range from [Fe/H] = −2.26 to −0.15 dex, and a dispersion of 0.25-0.30dex.These studies also indicate that there is a slight gradient in metallicity as a function of distance from the Galactic center, with the central-most RRLs being more metal-rich than those as further Galactocentric distances.The most recent photometric metallicity analysis of bulge stars, by Dekany et al. (2021), finds that the mode of the bulge RRL MDF is approximately [Fe/H] ∼ −1.4 dex, and they also see a shift in the bulge RRL MDF to slightly toward lower metallicity with increasing Galactocentric distance [41].Pietrukowicz et al. (2020) suggest that the inner Galaxy MDF is better approximated by three peaks, and that each [Fe/H] peak represents a different Galactic component that is superimposed within the inner Galaxy-the halo, bulge and disk [46].There is no doubt that a complexity in determining the true MDF of bulge RRLs is the reality that all Galactic components overlap in the inner Galaxy, and, especially at the metal-poor end ([Fe/H] ∼ −1.0 dex), there will be a significant contamination from halo stars.It is now possible to clean the inner Galaxy RRLs from outliers, as has recently been done thanks to Gaia proper motions, and ongoing radial velocity surveys, as described in Section 4.4.
For photometric metallicities, both a homogeneously derived set of chemical abundances for a large sample of RR Lyrae stars with high-resolution, high signal-to-noise ratio spectra as well as inclusion of more RR Lyrae stars in the calibrating samples at both low ([Fe/H] ≤ −2.0) and high ([Fe/H] ≥ −0.3) metallicities would more robustly anchor photometric parameters and the trends with metallicity-in addition to being able to distinguish between the effects of [Fe/H] and [α/Fe] would be important when moving toward RRL in the bulge.For spectroscopic metallicities, it would be valuable to have more spectroscopic [Fe/H] studies, and, especially, being able to measure elemental abundances of bulge RRLs spectroscopically would be exciting.

Age-Metallicity Relation
The dominant stellar populations of the bulge appear to be largely old (∼11 Gyr) [47,48], with little evidence of younger stars, although there have been also claims of a substantial intermediate-age population, especially from microlensing observations [49,50].
A landmark theoretical calculation by Lee (1992) [10] estimated that the RR Lyrae in the inner Milky Way were ∼1.5 Gyr older than typical Halo stars with [Fe/H] ∼ −1.50.It is still widely accepted that the bulge RRLs are some of the oldest stars in the MW.RRLs can be used as age indicators because RRLs are horizontal branch (HB) stars, and the HB morphology changes dramatically as a function of metallicity.Therefore, the metallicity of RRLs yields constraints on this population's age [10,28].
For example, using HB population models from the BaSTI alpha-enhanced isochrones, complemented with the BaSTI alpha-enhanced HB tracks [51], stars with an age of 10 Gyr and a [Fe/H] metallicity of −1.62 dex will end up as red horizontal branch star-too red to populate the instability strip and become an RRL.However, stars with an age of 10 Gyr and a [Fe/H] metallicity of −2.14 dex, will evolve to become a horizontal branch star blue enough to fall within the instability strip and will therefore pulsate as an RRL.At a given age, the more metal-poor a star, the bluer the star will be on the HB.
By increasing the age of the [Fe/H] = −1.62 dex stars to 13 Gyr, a stellar evolutionary path results that takes the star through the instability strip.Whereas different HB population models as well as the specific contribution of helium and heavy elements will factor into the exact position a star falls on the HB with age, metal-rich stars need to be older in age to become blue enough to fall within the instability strip (see, e.g., Figure 6 in Savino et al. 2020) [28].
In Savino et al. (2020), a set of synthetic stellar population models was built to match the modern metallicity distribution of the bulge RRLs [28].They concluded that the bulge RRL population has an extremely ancient age of 13.41 ± 0.54 Gyr, which is the oldest among the Milky Way stellar populations to date.
The ages of the bulge RRLs are compared in Figure 2 to bulge field stars determined from Hasselquist et al. (2020) [52].In Hasselquist et al. (2020), ages of ∼6000 bulge stars with [Fe/H] > −0.5 stars from the APOGEE survey were calculated using the Cannon labeltransfer method [52].They find that nearly all stars have ages consistent with being formed 8-10 Gyr ago, with a slight vertical age gradient in a sense that stars closer to the plane are younger than those further from the plane.Although ages of stars are notorious for having large uncertainties of ∼25-50%, these ages are similar to most studies of bulge stars that indicate that all, or nearly all, of the stars appear to be old->9-10 Gyr [48,[53][54][55].However, there is also evidence for significant fractions of bulge stars with ages <8 Gyr [50, [56][57][58].The Hasselquist et al. (2020) analysis explains younger stars in the bulge as the vertical age gradient suggests that younger stars will preferentially be found at distances closer than 0.5 kpc from the plane and with [Fe/H] > 0.2.

Spatial Distribution of RR Lyrae Stars in the Inner Galaxy
The appearance of the inner Galaxy is dominated by a bulge of stars, which is asymmetric and elongated, so that the bulge is a peanut-shaped bar pointing toward us [59,60].In this section, we discuss how RRLs are used as distance indicators to probe the inner Galaxy and the bar/bulge.

Distances Derived from RR Lyrae Stars
It is typically believed that, in the V-band, the absolute magnitude of an RRL is a function of metallicity [61,62], and traditionally a linear relationship between [Fe/H] and M V (RR) is assumed, e.g., M V (RR) = α[Fe/H] + β.Because this relation of RRLs has been arguably the most widely discussed RRL luminosity relation in the literature, it is briefly reviewed here.Theoretical models do not predict that a simple linear relationship between M V (RR) and [Fe/H] exists [63][64][65][66].For example, using computations of convective pulsation models, it was shown in Caputo et al. (2000) that the slope of the M V (RR) and [Fe/H] relation will increase as RRLs become more metal-rich.Using RRL in globular clusters, they show that observationally this effect can be seen as well.Despite these limitations in a M V (RR) -[Fe/H] relation, the slope (α) and zero-point (β) of this relationship are the subject of much discussion.
The Hipparcos satellite was only able to provide a parallax with an accuracy of <20% for one RRL (RR Lyrae itself), so different techniques, such as Statistical Parallax [67], Baade-Wesselink analysis [68], and theoretical HB models [64], have been used to calibrate the M V (RR).Values for the slope α range from 0.30-0.37[69,70] to 0.13 [71].This means that the M V (RR) at [Fe/H] = −1.6 ranges from 0.45 to 0.77 mag, depending on the technique used.One popular value for α is ∼0.21-0.25 mag/dex, which is not only the best estimate of the M V (RR) calibration from a number of different techniques [72] but is also in line with the value found from a photometric and spectroscopic study of about a hundred RRLs in the bar of the LMC from Clementini et al. (2003) and Gratton et al. (2004) [73,74].
It is exciting that today Gaia parallaxes of a statistical sample of RRLs can be used to investigate the absolute magnitudes of RRLs.Dambis et al. (2013) collected and homogenized literature values of period, pulsation mode and extinction in the visual passband (A V ), metal abundance ([Fe/H]) and intensity-averaged magnitudes in the Johnson V, 2MASS K s and W ISE W1 passbands for 403 field RRLs [75].The majority of these have Gaia DR2 parallaxes with uncertainties <20%.Using Gaia DR2 parallaxes and 381 RRLs from the Dambis et al. (2013) compilation, a surprisingly large value of α = 0.34 ± 0.03 mag/dex is found [76].When restricting the sample to 23 RRLs with metallicity estimated from high-resolution spectroscopy, and which are closer to us and span a narrower range of the distances, the slope of the M V (RR) -[Fe/H] decreases to α = 0.25 ± 0.05 mag/dex, more in line with previous estimates.It should be noted that this analysis takes into account the correction for the zero point parallax offset of Gaia DR2 data.It is therefore still not settled what the best relation is between M V (RR) and [Fe/H], but the situation will continually be revisited as Gaia parallaxes become more and more precise [77].
In recent years, there has been a growing body of literature exploring RRL luminosity relations in other passbands.In particular, there are advantages of using longer wavelengths of light when using RRLs as standard candles.Not only is longer-wavelength light less sensitive to reddening, but passbands such as the infrared K-band are shown to have less dependence on metallicity and also less dependence on evolutionary effects.Instead, an RRLs absolute magnitude M K (RR) depends mainly on the RRLs period [78,79].Such a period-luminosity (PL) relation is expected at nearly all wavelengths, with the exception of the V-band, as shown by theoretical models [64,80,81] Therefore, currently RRL PL relations are well characterized from optical (R and I) to mid-IR wavelengths [82][83][84][85][86][87][88].Combining period-luminosity relations with a metallicity term in general decreases the intrinsic scatter in finding M K (RR); with Gaia DR2, the scatter in period-luminosity-[Fe/H] relations is on the order of 0.2 mag, still dominated by uncertainties in the parallaxes [89,90].In the PLZ relation with the form M K s (RR) = a + b logP + c [Fe/H], the slopes in the literature in log P range from b = −2.4 to −2.9 and the slopes in [Fe/H] range from c = 0.03 to 0.18.Theoretical PLZ relations predict the slope in log P to range from b = −2.2 to −2.4 and the slope in [Fe/H] to be c ∼ 0.18 [64,81].As the PLZ relation becomes tighter constrained with forthcoming, improved Gaia parallaxes, it will be interesting to see if the slopes found observationally (b = −2.4 to −2.9) are consistently more negative than what the theoretical models predict (b = −2.2 to −2.4).
Regarding RRL in the inner Galaxy, the OGLE survey is the largest and most comprehensive publicly available RRL photometric catalog.Since OGLE is carried out in the V and I passbands, luminosity relations in M V and M I are needed for the many investigations that use the OGLE RRL as their foundation.Distances of RRL toward the Galactic center are typically found using the theoretical relations by Catelan et al. (2004b), given by (1) where P is the stars period, and Z is the RRL metal abundance [64].Distances for the RRL from the VVV survey, which uses the JHK S passbands, use a variety of different M J and M K S relations, mainly including those from Catelan et al.

Inner Galaxy RR Lyrae Stars Spatial Distribution
Gaia parallaxes of RRLs alone are not sufficient for RRLs at the distance of the bulge.Distances to these stars are therefore determined adopting absolute magnitudes from the relations discussed above (Section 3.1).All-sky reddening maps, such as those from the VVV [98,99], can be used to enable corrections of the large and variable extinctions present in the inner Galaxy.Still, reddening variations down to a few arcseconds are clearly seen in maps of specific windows in the inner Galaxy [98][99][100], and interstellar extinction curve variations can exceed 20 percent [101].It is therefore difficult for bulge RRL distances to be established with uncertainties of <10%, although, for some inner Galaxy GCs with well-studied extinctions and/or small differential reddening, RRL distances with a 5-8% precision have been reported [96].
From a sample of ∼4000 MACHO RRLs, it was reported that the inner Galaxy RRLs represent primarily the inner extension of the Galactic halo and not the bulge/bar [42,102].This is because it was found that the spatial distribution of the inner Galaxy RRLs did not appear to trace out the bar in the bulge, whereas the more metal-rich and more populous red clump stars show a strong barred signature [42,102,103].
In contrast, the OGLE survey, using tens of thousands of inner Galaxy RRLs, reported that the bulge RRLs do trace the barred bulge, with a stronger signature for the stars closest to the Galactic plane [43,44].These results favor that the inner Galaxy RRLs lie in the same gravitational potential as the bar.The OGLE photometry is more accurate for variable star studies than the MACHO photometry, which may be a factor in this discrepancy [104,105].However, the VISTA Variables in the Vía Láctea (VVV) survey, using ∼8000 OGLE RRLs, have also reported a spheroidal, barely barred distribution [91,106].The VVV survey has the advantage of being less sensitive to reddening than those in the optical passbands, such as OGLE and MACHO.
The completeness of the OGLE survey has been reported to be about 99% [107] based on the cross-identification with the MACHO RRLs.Therefore, this comparison is valid for RRL in the southern half of the bulge and those a few degrees from the Galactic plane (b < −3 • ).We note that OGLE and MACHO are both ground-based photometric surveys carried out in the optical passbands, so both will suffer from the same biases regarding crowding and reddening.The completeness of the OGLE sample falls with decreasing latitude, as at |b| <2 • , the apparent magnitudes of RRLs are heavily extincted and hence too faint for detection within OGLE and MACHO.
The VVV survey has by and large taken advantage of the RRLs identified in OGLE, and has supplemented these with their JHK s photometry [91].For areas not able to be probed by OGLE, VVV has been successful in identifying thousands of RRLs, although the RRL along the slight-lines of dust clouds in the close proximity of the mid-plane still remain uncovered [106,108,109].
To date, there is no consensus on the exact 3D structure of the inner Galaxy RRLs.Recently, Kunder et al. (2020) argue that there are two populations of bulge RRLs, one population that follows the bar and is dominant at Galactocentric distances larger than ∼1 kpc, and another that is more centrally concentrated that does not follow the bar [110].These results are based on ∼1000 RRLs with known radial velocities, proper motions, and extinction estimates.
That there are two populations of bulge RRLs-one population that is part of the bar, and one that does not have a bar-like signature-would suggest that the different RRL spatial distributions seen by OGLE and VVV are both correct, and that the differences arise due to both surveys focussing on RRLs in slightly different locations within the bulge.In particular, the results from the OGLE papers would be drawn from distances of RRL dominating further from the plane, and the VVV papers would be incorporating more distances of RRL closer to the plane.Because the longer passbands from the VVV photometry are less sensitive to dust, VVV has the advantage of being less hampered by reddening, especially for stars with a closer proximity to the plane.This could in turn give greater accuracy in the VVV distances, especially for RRL closer to the plane.
The BRAVA-RR survey is not thought to be affected by any significant biases.Their sample is selected from OGLE, which has a high completeness in the Southern fields that were selected by BRAVA-RR, and the full RRL period-amplitude plane is sampled.The BRAVA-RR fields have RRLs with apparent magnitudes that not only reach the bulge, but also to that of the Sagittarius dwarf Galaxy at ∼25 kpc.
Figure 3 shows the spatial distribution of ∼4000 inner bulge RRLs as a function of Galactocentric distance.The top panel highlights the spatial location of these RRLs, as this is an area where contiguous celestial coverage can not be provided by OGLE; instead, the first statistical sample of RRLs in this region was uncovered using the VVV by Dekany et al. (2020) [109].These authors analyzed time series data from the the VVV survey to identify RRLs also in some of the more extinct regions of the Galaxy-an area that is too faint for optical surveys to probe.To find distances to these stars, the theoretical calibration of the RRL period-luminosity relation from the Catelan et al. (2004b) [64] on the VIRCAM/VISTA filter system is used: where P is the period of the RRL, and Z is its metallicity [111].
To correct for reddening, the Bulge Extinction And Metallicity (BEAM) extinction values [98] is used with the extinction law from Cardelli et al. (1989) [112], which in the VVV passbands is A Ks = 0.73E(J − K s ) [92].It is then straightforward to calculate the distance to the individual RRLs, where D the individual distance in parsecs.The spatial distribution of the VVV inner bulge RRL sample does not clearly appear to trace the bar/bulge.In Figure 3, the x and y distances of VVV inner Galaxy RRL are compared to the spatial distribution of APOGEE giants with distances from StarHorse [113] and VVV red clump stars studied by Gonzalez et al. (2012) [98] (which use the same BEAM reddening map).The bulk of the inner Galaxy RRL stars do not appear to trace the bar as clearly as the giants or red clump giants.Combined with the metal-poor nature of these RRLs (see Figure 1), it is apparent that the inner Galaxy RRLs warrant further investigations on how they fit into the bulge population, or if they are even a bulge population [114].The spatial properties of the inner Galaxy RRLs in combination with their kinematics is useful for a more detailed look as to their origin.

The First Inner Galaxy Kinematic Surveys
One of the first indications of rotation of the bulge came from radial velocity measurements of planetary nebulae (PNe) [115].The velocities of the PNe toward the Galactic center showed a large dispersion with a trend of increasing velocity with increasing longitude.Twenty-five years later, 26 bulge Miras were observed to show a rotation gradient, whereas disk Miras did not show such rotation [116].Then, ∼1000 K-giants in two off-axis fields at (l, b) = (8, 7) and (l, b) = (12, 3) were used as metal-rich bulge tracers to show a similar result [117].These observations therefore suggested rotation of the bulge/bar.
The first dynamical survey that covered a wide area of the bulge was carried out by Beaulieu et al. (2000) [118].They combined PNe from a new survey with those already known from existing catalogs to obtain a sample of ∼400 PNe that cover both the Northern and Southern portions of the bulge (−20 • < l < 20 • and 3.3 • < |b| < 15 • ).The velocities of the PNe in this survey clearly indicated rotation of the bulge.It was further found that N-body models, in which the bar/bulge grows from the disk via bar-forming instabilities, give a good representation of these bulge PNe kinematics [118].
PNe are short-lived and are therefore rare, and their population memberships (and especially distances) are uncertain.A more extensive survey of the bulge using red giant stars, which are the most common type of luminous evolved stars in the bulge/bar, was undertaken by Rich (2007) [119].In this Bulge Radial Velocity Assay (BRAVA), eventually 10,000 stars in 84 different fields across the southern bulge were observed [120].There it was demonstrated that the bulge follows pure cylindrical rotation [121], similar to boxy/peanut bulges in external galaxies [122].It was also shown in Shen et al. (2010) that the BRAVA kinematics (Figure 2 in Shen et al. 2010) [123] is consistent with most of the mass being in the stellar bar, with <8% (of the mass of the disk) allowed within a slow/non-rotating pressure-supported classical component.
The Abundances and Radial velocity Galactic Origins Survey (ARGOS) Survey used ∼28,000 red clump giant stars (RCGs) over 28 fields to probe the chemistry and dynamics of the bulge [124,125].RCGs are representative for most of the bulge stars, tracing old stellar populations except for low metallicities [126].In the ARGOS survey, ∼95% of the sample contains stars with with metallicities [Fe/H] > −1.0.The near-cylindrical rotation is seen for RGC stars of all metallicities up to [Fe/H] ∼ −1.
ARGOS and BRAVA by and large avoided the plane of the Galaxy; the plane is plagued by high and variable reddening, as well as stellar blending and crowding issues.To fill in this gap, the Giraffe Inner Bulge Survey (GIBS) [29,127] targeted ∼5000 red clump giants in 24 bulge fields, to probe both elemental abundances and velocities of stars also close to the plane of the Galaxy (b = −1).GIBS confirmed cylindrical rotation also along the plane of the Galaxy and showed that differences may exist between super-solar and sub-solar bulge stars.In particular, they find that the metal poor population ( [Fe/H] ∼ −0.3) is more centrally concentrated than the metal rich one ( [Fe/H] ∼ +0.1), and with a more axisymmetric spatial distribution.The metal rich population, on the other hand, is arranged in a boxy distribution, consistent with an edge-on bar.
The APOGEE survey [128,129] also is collecting spectroscopic observations of stars along the plane of the Galaxy.From ∼10,000 giants at low-latitudes, the kinematic profile of the bulge is found with stars with [Fe/H] > −1.0 having cylindrical rotation.
From the Gaia-ESO Survey (GES) [130,131], ∼2500 red clump stars in 11 bulge fields were compared to ∼6300 disk stars to study the chemo-dynamics of the bulge.They found differences in dynamics of the metal-rich stars and the metal-poor subpopulation.
Figure 4 shows the extent of the bulge covered by these surveys and the metal-poor surveys described in the following section.Because the Northern bulge tends to be more crowded than the Southern bulge, this area of the sky is not as well studied.There is no indication that the Northern bulge has different global kinematic properties than the Southern bulge, and it is generally assumed that the bulge exhibits symmetry about the Galactic plane [132,133].

Kinematics of Stars with [Fe/H] > −1
Whereas there is disagreement as to exact metallicity distribution of the bulge stars from different bulge surveys (see, e.g., Figure 4 in Barbuy et al. 2018), there is a remarkable agreement regarding the distribution of radial velocities of bulge stars [55].Figure 5 (left panel) shows the Galactocentric velocities as a function of longitude for the BRAVA, GIBS, APOGEE and ARGOS fields within |b| < 5 • .All surveys show that the bulge is rotating, as the radial velocity is positive (stars are, on average, going away from the Galactic center) at positive longitudes and negative (stars are on average going toward the Galactic center) at negative longitudes. Figure 5 (right panel) shows only BRAVA velocities for all Galactic latitudes.It is evident that the Galactocentric velocities do not change as a function of latitude or height above the disk plane-this same result is found for other bulge kinematic surveys as well.This is what is expected for the bulge rotating in a cylindrical manner [121].If the bulge was rotating spherically (e.g., similar to the Sun or the Earth), the velocities would be faster near the Galactic plane (the mid-point) and slower closer to the poles.Cylindrical rotation of bulges is also seen in most other boxy bulges of other galaxies, and so this rotation is strong evidence that the MW bulge is an edge-on bar [134].
The velocity dispersion is not the same at all latitudes.Instead, the velocity dispersion drops with latitude, indicating that the stellar populations near the Galactic plane are kinematically hotter than the stellar populations further from the Galactic plane.This is expected from N-body models in which the Galaxy self-consistently develops a bar [123].
N-body models of bars also rotate cylindrically when viewed edge-on [135][136][137].Therefore, all of these surveys agree that the dynamics are consistent with a boxy bulge being formed via dynamical instabilities in massive early disks [123,138].Although a radial metallicity gradient observed in bulges was at first considered a signature of a classical bulge, it has now been shown that metallicity gradients may be reproduced in these models [139].This suggests that our galaxy has not experienced a significant merger since at least the epoch of disk formation (z ∼3), since otherwise one expects a classical bulge to result [124].
Using observations from the GIBS survey, Zoccali et al. (2017) showed that, upon dividing bulge red giant stars by [Fe/H] ∼0, those that have [Fe/H] <0 (i.e., the subsolar population) are more centrally concentrated than those with [Fe/H] >0 (i.e., the metal rich one), and with a more axisymmetric spatial distribution [29].The more metal rich population has a more boxy distribution, consistent with an edge-on bar.Curiously, the kinematics of the two populations show no statistical difference, so it is unclear if these spatial differences arise from having different origins.Figure 6 shows the different [Fe/H] metallicity components found in the GIBS survey, but that these do not have different rotation signatures.Rojas-Arriagada et al. (2019) argue that the apparent spatial dichotomy for subsolar and solar metallicity bulge stars does not translate into a dichotomy in α-abundance space [30].At least within the solar and subsolar bulge stellar population, there is a limited possibility for different origin scenarios, and hence no strong indication that accretion played a relevant role in the formation of the bulge.
Recent discovery of accretion events that undoubtedly extend to the inner Galaxy include the Inner Galaxy Structure [140] and the Kraken [141], accretion events that occurred around z = 2.2, and which may be the same event [140], Gaia-Enceladus [39,142], an accretion event that occurred around z = 1.4,and the Sagittarius dwarf galaxy [143], an accretion event that occurred around z = 0.8.These have all been shown to be significant building blocks of the Galactic halo, contributing to the properties and build-up of the halo.The Milky Way halo also extends to the inner Galaxy but is more difficult to study there due to the high extinction in the inner Galaxy as well as the complications of disentangling inner halo stars from the foreground stars and bulge stars.Because the halo does have a connection to the formation of the disk [142], it is not implausible that there may be some signatures of these accretion events within bulge stars.To date, however, there is no connection reported between bulge stars and those from the halo accretion events listed above.

Kinematics of Stars with [Fe/H] < −1
Most of the recent spectroscopic surveys do not specifically target metal-poor stars: BRAVA and APOGEE use K and M giants as tracers and these largely consist of the dominant metal-rich stellar population as metal-poor stars do not become M giants.It is therefore not surprising that APOGEE has found very few stars with [Fe/H] < −1.ARGOS and GIBS use red clump giants, i.e., metal-rich core Helium burning stars.ARGOS identified but 16 stars with [Fe/H] < −2.0 dex out of 14,000, while, in the APOGEE survey, 10 stars with −1.8 < [Fe/H] < −0.8 are highlighted [144], still a considerably larger population than previously observed [145,146].The EMBLA survey, dedicated to searching for metal-poor stars in the bulge, has uncovered ∼40 metal-poor stars [147][148][149], including a handful with [Fe/H] < −3.0.This sample is too small for a large scale kinematic investigation, but they do show that at least some of these metal-poor stars are confined to the bulge and on tight orbits.
The Pristine Inner Galaxy Survey (PIGS) first data release [150], on the other hand, did target metal-poor stars and was successful in uncovering ∼600 stars in the inner Galaxy with [Fe/H] < −2.0 dex.With especially the PIGS sample, it seems that the bulge rotation curves are similar for different metallicity bins until [Fe/H] ∼ −1.5 dex.As the [Fe/H] becomes more metal-poor, the velocity dispersion increases.The rotation curve begins to slow with inclusion of stars more metal-poor than [Fe/H] ∼ −1.5 dex, and the velocity dispersion continues to increase.Figure 7 shows this trend from the PIGS survey.The reason for the slower rotation and higher velocity in the more metal-poor inner Galaxy stars could be due to a number of different factors.First, as the stellar population becomes more metal-poor, greater contamination from, e.g., the halo will enter the sample.Halo stars would naturally cause the bulge/bar rotation curve to slow and the velocity dispersion to be higher.Second, we could be seeing a change from the bulge/bar being dominated by a rotation-supported population of a bar or disk at higher metallicities ([Fe/H] > −1.0) to being fully dominated by the pressure-supported component, like a classical bulge at the lowest metallicities.Third, the metal-poor stars could be a few Gyr older than the more metal-rich stars, and therefore have been mapped into the boxy/peanut bulge in different ways because of their different initial velocity distributions at the time of the bar formation [151,152].This is commonly referred to as kinematic fractionation.
Whether the PIGS metal-poor stars are actually confined to the inner bulge could be tested through an orbital analysis.In this way, halo outliers could be removed, and a "clean" sample of bulge stars could be obtained.Unfortunately, whereas most of the PIGS stars will have proper motions form Gaia, uncertain distances have prevented an orbital analysis from being carried out.At the moment, therefore, it is not clear if these metal-poor stars are actually confined to the bulge, or an inner halo or halo interlopers.
Whether the PIGS metal-poor stars are a rotation-supported classical bulge could be investigated by exploring the differences in elemental abundances between the most metalpoor stars and the "normal" bulge stars.Recently, Arentsen et al. ( 2021) have compared the fraction of their carbon-enhanced metal-poor (CEMP) stars in their sample to other Galactic components [153].The PIGS inner Galaxy CEMP fraction for [Fe/H] < −3 is consistent with the halo fraction found in the literature.In contrast, at higher metallicities ([Fe/H] between −1.5 and −2.5), the PIGS fraction is substantially lower.This may indicate a fast chemical enrichment in the early building blocks of the inner Galaxy and/or a larger contribution of the low-metallicity bulge from former globular cluster stars.More highresolution follow-up of the PIGS CEMP stars, as well as an orbital analysis of PIGS stars, can make progress on this front.

Kinematics of RR Lyrae Stars
The first study of the kinematics of RRLs toward the inner Galaxy came from Kunder et al. (2016) [154].Using ∼1000 RRLs, it was found that the RRL kinematics do not follow the rotation pattern of any of the other bulge populations.They trace a slow or non-rotating component, possibly corresponding to the classical (pressure-supported) bulge.Du et al. (2020) use Gaia DR2 proper motions of ∼15,000 OGLE RRL to probe the rotation of the Galactic bulge [155].They confirm slower rotation of the inner Galaxy RRL, as first reported in the BRAVA-RR survey.Splitting their sample of RRLs into those that are more metal-poor than [Fe/H] < −1.0 and those more metal-rich than [Fe/H] > −0.5, they find that the more metal-rich RRLs ([Fe/H] > −0.5) rotate similarly to the bar, and also with a similar velocity dispersion.The metal-poor RRLs are those that rotate slowest.The Du et al. (2020) sample of RRLs do not have radial velocities, and the Kunder et al. (2016) sample of RRLs did not have proper motions, so neither survey was able to clean for halo stars, which will be more present in a metal-poor population, and will of course inflate any velocity dispersion and decrease any rotation signature.
Missing, therefore, are the full space velocities of the RRLs.This is a problem because not all RRLs in the direction of the Galactic bulge are actually confined to the bulge, as shown explicitly for especially the high-velocity RRLs [156].Other metal-poor stars located toward the Galactic bulge have also been shown to be interlopers from other Galactic components [148].Obtaining a "pure" sample of bulge stars based on photometry/distance and spatial position alone, which is the usual approach when selecting bulge targets, is not statistically probable for more metal-poor stellar populations.Prudil et al. (2019) were the first to derive orbits for a significant sample of bulge RRLs [37].They reported that the RRLs in the central region of the Galactic bulge are consistent with random motions expected for a classical bulge component.They further searched for differences in the bulge RRLs based on their position in a period-amplitude diagram, or Oosterhoff group [37].From ∼450 RRLs, they did not find significant kinematic differences between a RRL and its position in the period-amplitude diagram.
Using a sample of ∼1400 stars with radial velocities and Gaia DR2 proper motions, Kunder et al. (2020) carried out an orbital analysis to determine which RRLs actually belong to the bulge [110].Now that Gaia eDR3 proper motions are available, this analysis can be repeated using more precise proper motions from eDR3.There is furthermore a ∼40% greater overlap between a BRAVA-RR star and one with an eDR3 Gaia proper motion.Therefore, the sample of BRAVA-RR stars for which orbits can be computed is almost doubled-2628 RRLs.
Orbits of the BRAVA-RR + Gaia eDR3 sample are calculated employing a non-axisymmetric model for the Galactic gravitational potential, as in Kunder et al. (2020).The model has an axisymmetric background made by an exponential disk built from the superposition of three Miyamoto-Nagai potentials [157], and a Navarro-Frenk-White density profile [158] to model the dark matter halo, which has a circular velocity V 0 = 241 km s −1 at R 0 = 8.2 kpc.The RRLs are separated into two groups, based on their apocenter distances, R apo , and the maximum distance from the Galactic plane that the stars reach during their orbit, z max .
Figure 8 shows the rotation curve of the BRAVA-RR stars in the inner bulge, where the inner bulge is defined as those stars having R apo <3.5 kpc and z max <2.5 kpc.It is believed that the bulge/bar likely extends out to ∼5 kpc [159], but the focus is on the inner bulge to ensure the RRL sample is truly confined to the bulge.There is a striking difference in the rotation profile of the stars on the most centrally concentrated orbits (R apo <1.8 kpc) compared to those of the bar/bulge RRL that have apocenter distances between 1.8 kpc and 3.5 kpc.The central bulge RRLs show almost no rotation, even out to |l| ∼7 • .In contrast, the inner bulge RRLs show rotation, although the amplitude of rotation is smaller.The velocity dispersion of the two populations is also different, with the central bulge RRLs having a smaller velocity dispersion than those in the inner bulge.Kunder et al. (2020) illustrate that it is not only the kinematics that differ between the more centrally concentrated RRLs and those that are less tightly bound, but that these populations also have different spatial distributions.Figure 9 shows that, although neither of these populations of RRL show a strong bar/bar spatial distribution, the less tightly bound RRL follow the bar/bulge weakly.In contrast, the more central RRLs show almost no trace of a bar.Both of these RRL groups have similar [Fe/H] distributions.Therefore, within the metal-poor stars, two distinct spatial and kinematic groups can be seen based on how centrally confined they are within the bulge.As such, it may be more appropriate to separate bulge stars based on their orbits, rather than their [Fe/H] metallicity.Theoretically, it is also predicted that stellar populations with the oldest ages and origins will have centrally confined orbits compared to other bulge populations [160].

Kinematics of RR Lyrae Stars in the Context of Milky Way Formation Scenarios
Martinez-Valpuesta and Gerhard (2011) show that, when the disk evolves to form a boxy bulge through a buckling instability, the stellar disk can also extend to a longer in-plane bar [161].Therefore, secular evolution of a disk galaxy, a natural mechanism for galaxies like ours, can be responsible for a variety of different stellar features.Using a similar N-boday simulation, Shen et al. (2010) find that the velocities of the boxy bulge agree well with what is found from bulge radial velocity surveys [123], leaving little doubt that the boxy bulge was formed from secular evolution of the disk.
Regarding the chemical gradients observed in the bar/bulge, Martinez-Valpuesta and Gerhard (2013) further show that chemical gradients in the initial disk can transform into gradients in the boxy bulge.Therefore, depending on when the bulge formed and on the properties of the precursor disk at that time, boxy bulges can have stars with differences in [Fe/H] metallicities depending on their spatial location within the bulge.It is therefore expected that a mix of stars can co-exist in the boxy bulge, all resulting from secular evolution of the MW disk.Whether this is the case also for the old, metal-poor RRLs is still under debate.
One promising explanation for a barred population of stars co-existing with a population of stars that is either less-barred or un-barred is that the bulge stars formed from disk buckling and ended up with different spatial distributions due to having slightly different ages and initial velocity distributions, i.e., kinematic fractionation [152].In this scenario, the different stellar populations originate from the Galaxy's disk.They would have slightly different ages (difference of ∼1.5-2 Gyr) and different initial velocity distributions.As the Galaxy evolves and the disk buckles, they would end up tracing the barred bulge in different ways-the hotter and older initial population would weakly trace the bar, whereas the cooler and younger population would follow a stronger barred spatial distribution.
The issue with kinematic fractionation being an explanation for the RRL populations is that the weakly-barred, central RRLs have a larger velocity dispersion compared to the barred RRLs (see Figure 8).Another issue is that since, in general, [Fe/H] metallicity and age of an RRL are related (see Section 2.3), and there is no difference in the mean [Fe/H] distributions of the two populations (see Figure 9), it is difficult to conclude these two populations have age differences on the order of 1.5-2 Gyr.However, until the α abundances are sufficiently known, in particular, if the RRLs have a spread in [α/Fe], relative ages will have uncertainties at the ∼1 Gyr level.
That the ages of the RRLs are >10 Gyr also renders the cosmological hydrodynamical simulation of Buck et al. (2018) unable to explain the bulge RRL properties.In Buck et al. (2018), stars with ages >10 Gyr form a roughly spheroidal and nonbarred distribution of stars in the central few kiloparsecs of the Galaxy, and only stars with ages ∼8 Gyr or younger contribute to the bar [162].However, the inner Galaxy RRLs do show that there is a population of bulge RRLs that contribute to the bar, despite their old ages.
There has been indication from a number of studies that the metal-rich and metalpoor bulge stars show different spatial distribution [29], which could suggest that the metal-poor and metal-rich bulge stars trace different Galactic components.However, the photometric [Fe/H] distribution of the two groups of RRLs do not show a difference in mean photometric metallicity-both groups have an average [Fe/H] of −1.23 dex with a dispersion in [Fe/H] of 0.30 dex.Therefore, formation scenarios indicating metal-rich and metal-poor stars trace different bulge populations are not complete and/or too limited.Within the old, metal-poor bulge RRL (R apo < 3.5 kpc) different bulge populations are seen-one that traces out the bar with a similar spatial distribution and kinematic signature as seen in more metal-rich populations, and one that does not trace out the bar, with slower rotation.The bulge RRLs suggest that it is not the case that bulge metal-poor stars arose from one common origin.Saha et al. (2016) predict that a massive classical bulge will gain some angular momentum from the bar, and so classical bulge rotation will be small in the center, but significant at larger radii [163].Rotation at the larger longitudes that is absent in the center of the bulge could be interpreted as a population of stars belonging to a relatively massive classical bulge.Low-mass classical bulges would likely follow the cylindrical rotation even at the center [164].The RRLs do show some suggestion of rotation at larger longitudes compared to the center, but the sample size of the RRLs is still small.The Milky Way harboring a massive classical bulge would be surprising, especially given the paucity of studies, suggesting that a classical bulge exists in the MW.

Conclusions
The oldest stars in the Bulge that likely witnessed the earliest episodes of star formation and assembly of our Galaxy are key to our understanding of its mechanisms.For example, within CDM-based simulations, the oldest stars in the Galaxy, those formed in the primordial dark matter halos at z ∼10-15, are expected to be tightly bound in the center of the Milky Way as it forms from the inside out.These stars lie in the bulge, but do not share its abundance patterns or kinematics [160].The greatest chance of uncovering a classical component, formed by the early (z >10) collapse of baryons within dark matter halos, would therefore be within the old, metal-poor stars in the inner few kpc of our Galaxy.
Unfortunately, it has been difficult for studies to specifically target old stars, although the stellar population of the bulge appears to be predominantly old and approximately solar in metal abundance.For example, the K and M giants and red clump giants, used as tracers in the past and current Bulge surveys (BRAVA, ARGOS, GIBS and APOGEE) are a mixed population of stars of different ages and metal abundances.Surveys such as EMBLA and PIGS are now beginning to uncover larger samples of metal-poor stars in the inner Galaxy.Still, today there are only a few low metallicity stars known to reside in the bulge, and, from these samples, it is not possible to have a complete picture into the structure and properties of the (likely) oldest stars in the bulge of our Galaxy.
This review discusses a number of advantages in using the bulge RRL population to trace the old, metal-poor bulge.These stars are at least ∼10x more metal-poor than the subsolar metallicity stars probed by many other bulge surveys, (Section 2) they are distance indicators (Section 3), and it is well-known from stellar evolution theory that RRLs are old (>10 Gyr) core Helium burning stars.
There is still no consensus as to if the bulge RRLs spatially trace out the bar, or if they are more consistent with a more spheroidal bulge population; different RRL photometric surveys in different wavelengths (i.e., OGLE in the optical and VVV in the infrared) present differing results.There is a broad consensus from the RRL proper motions and the RRL radial velocities that the kinematics of the bulge RRLs rotate slower than the bar/bulge, although there does appear to be a subset of bulge RRLs that rotate in a similar manner to the bar [110,155,165].
Possibly the most intriguing aspect of the RRLs is that there are two different metalpoor, bulge populations that coexist together in the bulge [110].The population of RRLs on the tightest orbits (R apo < 1.8 kpc) have both a spatial distribution and kinematic signature unlike the bar.Instead, this population has an average Galactocentric velocity consistent with no rotation, and a small velocity dispersion.The inner Galaxy RRLs (R apo between 1.8 kpc and 3.5 kpc) are more consistent with following the bar, with a signature of rotation and an elongated spatial distribution consistent with the bar.
These observations can be compared in the context of bulge formation scenarios.In Saha et al. (2016), a classical bulge coexists with a pseudobulge [163]; in Debattista et al. (2017), kinematic fractionation naturally gives rise to bulge populations with different spatial and kinematic signatures [152], and in Buck et al. (2018), the oldest (>10 Gyr) stars have different properties than younger (∼8 Gyr and younger) stars [162].Although these scenarios naturally give rise to a bulge population that is barred and one that is not because the barred and un-barred RRL populations are so similar in [Fe/H] and therefore age, and because the velocity dispersion of the unbarred RRL is smaller than that of the barred RRL, none of these scenarios adequately interpret the RRL observations.It is therefore still not clear as to the origin of the bulge RRLs.One promising avenue to make progress on this front is to combine the RRL kinematics with reliable spectroscopic metallicities of bulge RRLs, in particular elemental abundances.In this way, the bulge RRLs may be able to be conclusively linked to the bar, or to inner Galaxy merger events (e.g., the IGS [140]), or to the dissolution of a globular cluster or dwarf Galaxy stars [166,167].

Figure 1 .
Figure 1.The normalized metallicity distribution for ∼900 BRAVA-RR inner Galaxy RRLs (shaded histogram) in comparison with the distribution of ∼2 million inner Galaxy red giants (solid thick line) selected from the Bulge DECam Bulge Survey catalog of Johnson et al. (2020)[25].The distribution of RRLs is considerably more metal-poor than the dominant bulge population.

Figure 2 .
Figure 2. Age distributions of the bulge field stars as compared to the ages of bulge RR Lyrae stars.

Figure 3 .
Figure 3.A comparison of the spatial properties of the RRLs in the inner Galaxy (red) compared to the giants in the inner Galaxy (blue).The top panel shows the location in Galactic coordinates of the APOGEE giants (blue) and the VVV RRLs from Dekany and Grebel (2020) [109] (red).The bottom panel shows the mean spatial distribution of these two inner Galaxy stellar traces.The red-clump stars studied in Gonzalez et al. (2012) [98] to trace the bar/bulge are shown as triangles.

Figure 5 .
Figure 5. (Top) The mean velocity profiles for the stars within |b| < 5 • in the ARGOS (magenta), APOGEE (light blue), BRAVA (green) and GIBS (grey) surveys as a function of Galactic longitude.The lines represent the model by Shen et al. (2010) with a pure bar and no classical component.(Bottom) The velocity dispersion of the same surveys as a function of Galactic longitude.The surveys all agree well, and are consistent with the bulge being in pure cylindrical rotation.There are no obvious signs of cold substructure, indicating the lack of recent mergers.

Figure 6 .
Figure 6.(Left) Metallicity versus Galactocentric velocity for the GIBS sample.A division between sub-solar ([Fe/H] < −0.1 dex) and solar ([Fe/H] > −0.1 dex) samples can be seen.(Right)The Galactocentric velocity versus longitude distribution at fixed latitudes, split up into sub-solar (grey) and solar (black) samples.There is no statistical difference in kinematics between these two metallicity groups.

Figure 7 .
Figure 7. (Top) Galactic longitude versus mean Galactocentric velocity for different metallicity ranges.(Middle row) the same, but separately for each observed field with colors indicating the latitude of the field.(Bottom row) Galactic longitude versus standard deviation of the velocities, figure from Arentsen et al. (2020) [150].Reproduced with permission from Arentsen, A., Monthly Notices of the Royal Astronomical Society Letters; published by Oxford University Press, 2020.

Figure 8 .
Figure 8. (Left) The mean velocity for RR Lyrae stars in the BRAVA-RR survey as a function of Galactic longitude.(Right) The velocity dispersion profile for the same stars.The lines represent the model by Shen et al. (2010) with a pure bar and no classical component.

Figure 9 .
Figure 9. (Left) The mean spatial distribution of the most centrally concentrated RRLs (R apo < 1.8 kpc) as compared to bulge RRLs with apocenter distances between 1.8 and 3.5 kpc.The spatial distribution of red clump stars from Gonzalez et al. (2012) that trace out the bar/bulge is also shown.(Right) The photometric [Fe/H] metallicity distribution for the RRLs in the inner bulge (green, filled histogram) as well as the RRLs in the central bulge (black, solid line histogram).There is no obvious [Fe/H] metallicity difference between these two RRL populations.