The zero point of the reddening toward the Large Magellanic Cloud (LMC) has been the subject of some dispute. Its uncertainty propagates as a systematic error for methods which measure the extragalactic distance scale through knowledge of the absolute extinction of LMC stars. In an effort to resolve this issue, we used three different methods to calibrate the most commonly used metric to predict LMC extinction, the intrinsic color of the red clump, $(V-I)_{RC,0}$. The first approach was to empirically calibrate the color zeropoints of the BaSTI isochrones over a wide metallicity range of ${\Delta}\rm{[Fe/H]} \approx 1.10$ using measurements of red clump stars in 47 Tuc, the Solar Neighborhood, and NGC 6791. From these efforts we also measure these properties of the Solar Neighborhood red clump, ($V-I$, $G_{BP}-K_{s}$, $G-K_{s}$, $G_{RP}-K_{s}$, $J-K_{s}$, $H-K_{s}$, $M_{I}$, $M_{Ks}$)$_{RC,0} =$ (1.02, 2.75, 2.18, 1.52, 0.64, 0.15, $-$0.23, $-$1.63). The second and third methods were to compare the observed colors of the red clump to those of Cepheids and RR Lyrae in the LMC which are standardizable crayons based on their periods and metallicities. With these three methods, we estimated the intrinsic color of the red clump of the LMC to be $(V-I)_{RC,0,\rm{LMC}} = \{0.90,0.93,0.93\}$ respectively, and similarly using the first and third method we estimated $ (V-I)_{RC,0,\rm{SMC}} = \{0.85,0.88 \}$ respectively for the Small Magellanic Cloud. We estimate the luminosities to be $M_{I,RC,\rm{LMC}}=-0.31$ and $M_{I,RC,\rm{SMC}}=-0.39$. We show that this has important implications for recent calibrations of the tip of the red giant branch in the Magellanic Clouds used to measure $H_0$.

Globular cluster progenitors may have been detected by \textit{HST}, and are predicted to be observable with \textit{JWST} and ground-based extremely-large telescopes with adaptive optics. This has the potential to elucidate the issue of globular cluster formation and the origins of significantly helium-enriched subpopulations, a problem in Galactic astronomy with no satisfactory theoretical solution. Given this context, we use model stellar tracks and isochrones to investigate the predicted observational properties of helium-enriched stellar populations in globular cluster progenitors. We find that, relative to helium-normal populations, helium-enriched (${\Delta}Y=+0.12$) stellar populations similar to those inferred in the most massive globular clusters, are expected, modulo some rapid fluctuations in the first $\sim$30 Myr, to be brighter and redder in the rest frame. At fixed age, stellar mass, and metallicity, a helium-enriched population is predicted to converge to being $\sim$0.40 mag brighter at $\lambda \approx 2.0\, {\mu}m$, and to be 0.30 mag redder in the \textit{JWST}-NIRCam colour $(F070W-F200W)$, and to actually be fainter for $\lambda \lesssim 0.50 \, {\mu}m$. Separately, we find that the time-integrated shift in ionizing radiation is a negligible $\sim 5\%$, though we show that the Lyman-$\alpha$ escape fraction could end up higher for helium-enriched stars.

We investigate the Fe, C, N, O, Mg, Al, Si, K, Ca, Ce and Nd abundances of 2283 red giant stars in 31 globular clusters from high-resolution spectra observed in both the northern and southern hemisphere by the SDSS-IV APOGEE-2 survey. This unprecedented homogeneous dataset, largest to date, allows us to discuss the intrinsic Fe spread, the shape and statistics of Al-Mg and N-C anticorrelations as a function of cluster mass, luminosity, age and metallicity for all 31 clusters. We find that the Fe spread does not depend on these parameters within our uncertainties including cluster metallicity, contradicting earlier observations. We do not confirm the metallicity variations previously observed in M22 and NGC 1851. Some clusters show a bimodal Al distribution, while others exhibit a continuous distribution as has been previously reported in the literature. We confirm more than 2 populations in $\omega$ Cen and NGC 6752, and find new ones in M79. We discuss the scatter of Al by implementing a correction to the standard chemical evolution of Al in the Milky Way. After correction, its dependence on cluster mass is increased suggesting that the extent of Al enrichment as a function of mass was suppressed before the correction. We observe a turnover in the Mg-Al anticorrelation at very low Mg in $\omega$ Cen, similar to the pattern previously reported in M15 and M92. $\omega$ Cen may also have a weak K-Mg anticorrelation, and if confirmed, it would be only the third cluster known to show such a pattern.

Radial velocities of 2768 fundamental mode RR Lyrae stars (RRLs) toward the Southern Galactic bulge are presented, spanning the southern bulge from -8 < l < +8 and -3 < b <-6. Distances derived from the pulsation properties of the RRLs are combined with Gaia proper motions to give constraints on the orbital motions of 1389 RRLs. The majority (~75%) of the bulge RRLs have orbits consistent with these stars being permanently bound to <3.5 kpc from the Galactic Center, similar to the bar. However, unlike the bulge giants, the RRLs exhibit slower rotation and a higher velocity dispersion. The higher velocity dispersion arises almost exclusively from halo interlopers passing through the inner Galaxy. We present 82 stars with space velocities > 500 km/s and find that the majority of these high-velocity stars are halo interlopers; it is unclear if a sub-sample of these stars with similar space velocities have a common origin. Once the 25% of the sample represented by halo interlopers is cleaned, we can clearly discern two populations of bulge RRLs in the inner Galaxy. One population of RRLs is not as tightly bound to the Galaxy (but is still confined to the inner ~3.5 kpc), and is both spatially and kinematically consistent with the barred bulge. The second population is more centrally concentrated and does not trace the bar. One possible interpretation is that this population was born prior to bar formation, as its spatial location, kinematics and pulsation properties suggest, possibly from an accretion event at high redshift.

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

We explore the fundamental relations governing the radial and vertical velocity dispersions of stars in the Milky Way, from combined studies of complementary surveys including GALAH, LAMOST, APOGEE, the NASA $Kepler$ and K2 missions, and $Gaia$ DR2. We find that different stellar samples, even though they target different tracer populations and employ a variety of age estimation techniques, follow the same set of fundamental relations. We provide the clearest evidence to date that, in addition to the well-known dependence on stellar age, the velocity dispersions of stars depend on orbital angular momentum $L_z$, metallicity and height above the plane $|z|$, and are well described by a multiplicatively separable functional form. The dispersions have a power-law dependence on age with exponents of 0.441$\pm 0.007$ and 0.251$\pm 0.006$ for $\sigma_z$ and $\sigma_R$ respectively, and the power law is valid even for the oldest stars. For the solar neighborhood stars, the apparent break in the power law for older stars, as seen in previous studies, is due to the anti-correlation of $L_z$ with age. The dispersions decrease with increasing $L_z$ until we reach the Sun's orbital angular momentum, after which $\sigma_z$ increases (implying flaring in the outer disc) while $\sigma_R$ flattens. The dispersions increase with decreasing metallicity, suggesting that the dispersions increase with birth radius. The dispersions also increase linearly with $|z|$. The same set of relations that work in the solar neighborhood also work for stars between $3<R/{\rm kpc}<20$. Finally, the high-[$\alpha$/Fe] stars follow the same relations as the low-[$\alpha$/Fe] stars.

Studies of the kinematics and chemical compositions of Galactic globular clusters (GCs) enable the reconstruction of the history of star formation, chemical evolution, and mass assembly of the Galaxy. Using the latest data release (DR16) of the SDSS/APOGEE survey, we identify 3,090 stars associated with 46 GCs. Using a previously defined kinematic association, we break the sample down into eight separate groups and examine how the kinematics-based classification maps into chemical composition space, considering only $\alpha$ (mostly Si and Mg) elements and Fe. Our results show that: (i) The loci of both in situ and accreted subgroups in chemical space match those of their field counterparts; (ii) GCs from different individual accreted subgroups occupy the same locus in chemical space. This could either mean that they share a similar origin or that they are associated with distinct satellites which underwent similar chemical enrichment histories; (iii) The chemical compositions of the GCs associated with the low orbital energy subgroup defined by Massari and collaborators is broadly consistent with an in situ origin. However, at the low metallicity end, the distinction between accreted and in situ populations is blurred; (iv) Regarding the status of GCs whose origin is ambiguous, we conclude the following: the position in Si-Fe plane suggests an in situ}origin for Liller 1 and a likely accreted origin for NGC 5904 and NGC 6388. The case of NGC 288 is unclear, as its orbital properties suggest an accretion origin, its chemical composition suggests it may have formed in situ.

There is growing evidence that the Galactic Center Excess identified in the $\textit{Fermi}$-LAT gamma-ray data arises from a population of faint astrophysical sources. We provide compelling supporting evidence by showing that the morphology of the excess traces the stellar over-density of the Galactic bulge. By adopting a template of the bulge stars obtained from a triaxial 3D fit to the diffuse near-infrared emission, we show that it is detected at high significance. The significance deteriorates when either the position or the orientation of the template is artificially shifted, supporting the correlation of the gamma-ray data with the Galactic bulge. In deriving these results, we have used more sophisticated templates at low-latitudes for the $\textit{Fermi}$ bubbles compared to previous work and the three-dimensional Inverse Compton (IC) maps recently released by the ${\tt GALPROP}$ team. Our results provide strong constraints on Millisecond Pulsar (MSP) formation scenarios proposed to explain the excess. We find that an $\textit{admixture formation}$ scenario, in which some of the relevant binaries are $\textit{primordial}$ and the rest are formed $\textit{dynamically}$, is preferred over a primordial-only formation scenario at $7.6\sigma$ confidence level. Our detailed morphological analysis also disfavors models of the disrupted globular clusters scenario that predict a spherically symmetric distribution of MSPs in the Galactic bulge. For the first time, we report evidence of a high energy tail in the nuclear bulge spectrum that could be the result of IC emission from electrons and positrons injected by a population of MSPs and star formation activity from the same site.

We report the serendipitous discovery of a nitrogen-rich, mildly metal-poor ([Fe/H]=-1.08) giant star in a single-lined spectroscopic binary system found in the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) survey, Data Release 14 (DR14). Previous work has assumed that the two percent of halo giants with unusual elemental abundances have been evaporated from globular clusters, but other origins for their abundance signatures, including binary mass transfer, must also be explored. We present the results of an abundance re-analysis of the APOGEE-2 high-resolution near-infrared spectrum of 2M12451043+1217401 with the Brussels Automatic Stellar Parameter (BACCHUS) automated spectral analysis code, and re-derive manually the main element families, namely the light elements (C, N), elements (O, Mg, Si), iron-peak element (Fe), \textit{s}-process element (Ce), and the light odd-Z element (Al). Our analysis confirm the N-rich nature of 2M12451043+1217401, which has a [N/Fe] ratio of $+0.69$, and shows that the abundances of C and Al are slightly discrepant from that of a typical mildly metal-poor RGB star, but exhibit Mg, Si, O and \textit{s}-process abundances (Ce) of typical field stars. We also detect a particularly large variability in its radial velocity over the period of the APOGEE-2 observations, and the most likely orbit fit to the radial velocity data has a period of $730.89\pm106.86$ days, a velocity semi-amplitude of $9.92 \pm 0.14$ km s$^{-1}$, and an eccentricity of $\sim 0.1276 \pm0.1174$, which support the hypothesis of a binary companion, and that has probably been polluted by a now-extinct AGB star.

The next decade affords tremendous opportunity to achieve the goals of Galactic archaeology. That is, to reconstruct the evolutionary narrative of the Milky Way, based on the empirical data that describes its current morphological, dynamical, temporal and chemical structures. Here, we describe a path to achieving this goal. The critical observational objective is a Galaxy-scale, contiguous, comprehensive mapping of the disk's phase space, tracing where the majority of the stellar mass resides. An ensemble of recent, ongoing, and imminent surveys are working to deliver such a transformative stellar map. Once this empirical description of the dust-obscured disk is assembled, we will no longer be operationally limited by the observational data. The primary and significant challenge within stellar astronomy and Galactic archaeology will then be in fully utilizing these data. We outline the next-decade framework for obtaining and then realizing the potential of the data to chart the Galactic disk via its stars. One way to support the investment in the massive data assemblage will be to establish a Galactic Archaeology Consortium across the ensemble of stellar missions. This would reflect a long-term commitment to build and support a network of personnel in a dedicated effort to aggregate, engineer, and transform stellar measurements into a comprehensive perspective of our Galaxy.

The problem of the origin of the elements is a fundamental one in astronomy and one that has many open questions. Prominent examples include (1) the nature of Type Ia supernovae and the timescale of their contributions; (2) the observational identification of elements such as titanium and potassium with the $\alpha$-elements in conflict with core-collapse supernova predictions; (3) the number and relative importance of r-process sites; (4) the origin of carbon and nitrogen and the influence of mixing and mass loss in winds; and (5) the origin of the intermediate elements, such as Cu, Ge, As, and Se, that bridge the region between charged-particle and neutron-capture reactions. The next decade will bring to maturity many of the new tools that have recently made their mark, such as large-scale chemical cartography of the Milky Way and its satellites, the addition of astrometric and asteroseismic information, the detection and characterization of gravitational wave events, 3-D simulations of convection and model atmospheres, and improved laboratory measurements for transition probabilities and nuclear masses. All of these areas are key for continued improvement, and such improvement will benefit many areas of astrophysics.

Radial velocities of fifteen double-mode bulge RR Lyrae (RR01) stars are presented, six of which belong to a compact group of RR01 stars in pulsation space, with the ratio of first-overtone period to fundamental-mode period, P_{fo}/P_{f}~0.74, and P_{f}~0.44. It has been suggested that these pulsationally clumped RR01 stars are a relic of a disrupted dwarf galaxy or stellar cluster, as they also appear to be spatially coherent in a vertical strip across the bulge. However, the radial velocities of the stars presented here, along with proper motions from Gaia DR2, show a large range of radial velocities, proper motions and distances for the bulge RR01 stars in the pulsation clump, much larger than the RR01 stars in the Sagittarius dwarf galaxy (Sgr). Therefore, in contrast to the kinematics of the RRL stars belonging to Sgr, and those in and surrounding the bulge globular cluster NGC~6441, there is no obvious kinematic signature within the pulsationally clumped RR01 stars. If the pulsationally clumped RR01 stars belonged to the same system in the past and were accreted, their accretion in the inner Galaxy was not recent, as the kinematic signature of this group has been lost (i.e., these stars are now well-mixed within the inner Galaxy). We show that the apparent spatial coherence reported for these stars could have been caused by small number statistics. The orbits of the RR01 stars in the inner Galaxy suggest they are confined to the innermost ~4~kpc of the Milky Way.

Asteroseismology is a promising tool to study Galactic structure and evolution because it can probe the ages of stars. Earlier attempts comparing seismic data from the {\it Kepler} satellite with predictions from Galaxy models found that the models predicted more low-mass stars compared to the observed distribution of masses. It was unclear if the mismatch was due to inaccuracies in the Galactic models, or the unknown aspects of the selection function of the stars. Using new data from the K2 mission, which has a well-defined selection function, we find that an old metal-poor thick disc, as used in previous Galactic models, is incompatible with the asteroseismic information. We show that spectroscopic measurements of [Fe/H] and [$\alpha$/Fe] elemental abundances from the GALAH survey indicate a mean metallicity of $\log (Z/Z_{\odot})=-0.16$ for the thick disc. Here $Z$ is the effective solar-scaled metallicity, which is a function of [Fe/H] and [$\alpha$/Fe]. With the revised disc metallicities, for the first time, the theoretically predicted distribution of seismic masses show excellent agreement with the observed distribution of masses. This provides an indirect verification of the asteroseismic mass scaling relation is good to within five percent. Using an importance-sampling framework that takes the selection function into account, we fit a population synthesis model of the Galaxy to the observed seismic and spectroscopic data. Assuming the asteroseismic scaling relations are correct, we estimate the mean age of the thick disc to be about 10 Gyr, in agreement with the traditional idea of an old $\alpha$-enhanced thick disc.

We investigate aluminum abundance variations in the stellar populations of globular clusters using both literature measurements of sodium and aluminum and APOGEE measurements of nitrogen and aluminum abundances. For the latter, we show that the Payne is the most suitable of the five available abundance pipelines for our purposes. Our combined sample of 42 globular clusters spans approximately 2 dex in [Fe/H] and 1.5 dex in $\log{M_{GC}/M_{\odot}}$. We find no fewer than five globular clusters with significant internal variations in nitrogen and/or sodium with little-to-no corresponding variation in aluminum, and that the minimum present-day cluster mass for aluminum enrichment in metal-rich systems is $\log{M_{GC}/M_{\odot}} \approx 4.50 + 2.17(\rm{[Fe/H]}+1.30)$. We demonstrate that the slopes of the [Al/Fe] vs [Na/Fe] and [Al/Fe] vs [N/Fe] relations for stars without field-like abundances are approximately log-linearly dependent on both the metallicity and the stellar mass of the globular clusters. In contrast, the relationship between [Na/Fe] and [N/Fe] shows no evidence of such dependencies. This suggests that there were (at least) two classes of non-supernovae chemical polluters that were common in the early universe, and that their relative contributions within globular clusters somehow scaled with the metallicity and mass of globular clusters. The first of these classes is predominantly responsible for the CNO and NeNa abundance variations, and likewise the second for the MgAl abundance variations. 47 Tuc and M4 are particularly striking examples of this dichotomy. As an auxiliary finding, we argue that abundance variations among Terzan 5 stars are consistent with it being a normal globular cluster.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

This work explores the detailed chemistry of the Milky Way bulge using the HERMES spectrograph on the Anglo-Australian Telescope. Here we present the abundance ratios of 13 elements for 832 red giant branch and clump stars along the minor bulge axis at latitudes $b=-10^{\circ}, -7.5$ and $-5^{\circ}$. Our results show that none of the abundance ratios vary significantly with latitude. We also observe {\color{red}disk-like} [Na/Fe] abundance ratios, which indicates the bulge does not contain helium-enhanced populations as observed in some globular clusters. Helium enhancement is therefore not the likely explanation for the double red-clump observed in the bulge. We confirm that bulge stars mostly follow abundance trends observed in the disk. However, this similarity is not confirmed across for all elements and metallicity regimes. The more metal-poor bulge population at [Fe/H] $\lesssim -0.8$ is enhanced in the elements associated with core collapse supernovae (SNeII). In addition, the [La/Eu] abundance ratio suggests higher $r$-process contribution, and likely higher star formation in the bulge compared to the disk. This highlights the complex evolution in the bulge, which should be investigated further, both in terms of modelling; and with additional observations of the inner Galaxy.

To better understand the origin and evolution of the Milky Way bulge, we have conducted a survey of bulge red giant branch and clump stars using the HERMES spectrograph on the Anglo-Australian Telescope. We targeted ARGOS survey stars with pre-determined bulge memberships, covering the full metallicity distribution function. The spectra have signal-to-noise ratios comparable to, and were analysed using the same methods as the GALAH survey. In this work, we present the survey design, stellar parameters, distribution of metallicity and alpha-element abundances along the minor bulge axis at latitudes $b$ = $-10^{\circ}, -7.5^{\circ}$ and $-5^{\circ}$. Our analysis of ARGOS stars indicates that the centroids of ARGOS metallicity components should be located $\approx$0.09 dex closer together. The vertical distribution of $\alpha$-element abundances is consistent with the varying contributions of the different metallicity components. Closer to the plane, alpha abundance ratios are lower as the metal-rich population dominates. At higher latitudes, the alpha abundance ratios increase as the number of metal-poor stars increases. However, we find that the trend of alpha-enrichment with respect to metallicity is independent of latitude. Comparison of our results with those of GALAH DR2 revealed that for [Fe/H] $\approx -0.8$, the bulge shares the same abundance trend as the high-$\alpha$ disk population. However, the metal-poor bulge population ([Fe/H] $\lesssim -0.8$) show enhanced alpha abundance ratios compared to the disk/halo. These observations point to fairly rapid chemical evolution in the bulge, and that the metal-poor bulge population does not share the same similarity with the disk as the more metal-rich populations.

The Wide Field Infrared Survey Telescope (WFIRST) will monitor $\sim 2$ deg$^2$ toward the Galactic bulge in a wide ($\sim 1-2~\mu$m) W149 filter at 15-minute cadence with exposure times of $\sim$50s for 6 seasons of 72 days each, for a total $\sim$41,000 exposures taken over $\sim$432 days, spread over the 5-year prime mission. This will be one of the deepest exposures of the sky ever taken, reaching a photon-noise photometric precision of 0.01 mag per exposure and collecting a total of $\sim 10^9$ photons over the course of the survey for a W149$_{\rm AB}\sim 21$ star. Of order $4 \times 10^7$ stars will be monitored with W149$_{\rm AB}$<21, and 10$^8$ stars with W145$_{\rm AB}$<23. The WFIRST microlensing survey will detect $\sim$54,000 microlensing events, of which roughly 1% ($\sim$500) will be due to isolated black holes, and $\sim$3% ($\sim$1600) will be due to isolated neutron stars. It will be sensitive to (effectively) isolated compact objects with masses as low as the mass of Pluto, thereby enabling a measurement of the compact object mass function over 10 orders of magnitude. Assuming photon-noise limited precision, it will detect $\sim 10^5$ transiting planets with sizes as small as $\sim 2~R_\oplus$, perform asteroseismology of $\sim 10^6$ giant stars, measure the proper motions to $\sim 0.3\%$ and parallaxes to $\sim 10\%$ for the $\sim 6 \times 10^6$ disk and bulge stars in the survey area, and directly detect $\sim 5 \times 10^3$ Trans-Neptunian objects (TNOs) with diameters down to $\sim 10$ km, as well as detect $\sim 10^3$ occulations of stars by TNOs during the survey. All of this science will completely serendipitous, i.e., it will not require modifications of the WFIRST optimal microlensing survey design. Allowing for some minor deviation from the optimal design, such as monitoring the Galactic center, would enable an even broader range of transformational science.

The Kepler, K2 and TESS transit surveys are revolutionizing our understanding of planets orbiting close to their host stars and our understanding of exoplanet systems in general, but there remains a gap in our understanding of wide-orbit planets. This gap in our understanding must be filled if we are to understand planet formation and how it affects exoplanet habitability. We summarize current and planned exoplanet detection programs using a variety of methods: microlensing (including WFIRST), radial velocities, Gaia astrometry, and direct imaging. Finally, we discuss the prospects for joint analyses using results from multiple methods and obstacles that could hinder such analyses. We endorse the findings and recommendations published in the 2018 National Academy report on Exoplanet Science Strategy. This white paper extends and complements the material presented therein.

[abridged] WFIRST is uniquely capable of finding planets with masses as small as Mars at separations comparable to Jupiter, i.e., beyond the current ice lines of their stars. These planets fall between the close-in planets found by Kepler and the wide separation gas giants seen by direct imaging and ice giants inferred from ALMA observations. Furthermore, the smallest planets WFIRST can detect are smaller than the planets probed by RV and Gaia at comparable separations. Interpreting planet populations to infer the underlying formation and evolutionary processes requires combining results from multiple detection methods to measure the full variation of planets as a function of planet size, orbital separation, and host star mass. Microlensing is the only way to find planets from 0.5 to 5M_E at 1 to 5au. The case for a microlensing survey from space has not changed in the past 20 yrs: space allows wide-field diffraction-limited observations that resolve main-sequence stars in the bulge, which allows the detection and characterization of the smallest signals including those from planets with masses at least as small as Mars. What has changed is that ground-based (GB) microlensing is reaching its limits, underscoring the scientific necessity for a space-based survey. GB microlensing has found a break in the mass-ratio distribution at about a Neptune, implying that these are the most common microlensing planet and that planets smaller than this are rare. However, GB microlensing reaches its detection limits only slightly below the observed break. WFIRST will measure the shape of the mass-ratio function below the break by finding numerous smaller planets: 500 Neptunes, 600 gas giants, 200 Earths, and planets as small as 0.1M_E. Because it will also measure host masses and distances, WFIRST will also track the behavior of the planet distribution as a function of separation and host star mass.

Galactic interstellar dust has a profound impact not only on our observations of objects throughout the Universe, but also on the morphology, star formation, and chemical evolution of the Galaxy. The advent of massive imaging and spectroscopic surveys (particularly in the infrared) places us on the threshold of being able to map the properties and dynamics of dust and the interstellar medium (ISM) in three dimensions throughout the Milky Way disk and bulge. These developments will enable a fundamentally new understanding of dust properties, including how grains respond to their local environment and how those environments affect dust attenuation of background objects of interest. Distance-resolved maps of dust motion also hold great promise for tracing the flow of interstellar material throughout the Galaxy on a variety of scales, from bar-streaming motions to the collapse and dissolution of individual molecular clouds. These advances require optical and infrared imaging of stars throughout the Galactic midplane, stretching many kiloparsecs from the Sun, matched with very dense spectroscopic coverage to probe the ISM's fine-grained structure.

The overlap between the spectroscopic Galactic Archaeology with HERMES (GALAH) survey & $Gaia$ provides a high-dimensional chemodynamical space of unprecedented size. We present a first analysis of a subset of this overlap, of 7066 dwarf, turn-off, & sub-giant stars. [...] We investigate correlations between chemical compositions, ages, & kinematics for this sample. Stellar parameters & elemental abundances are derived from the GALAH spectra with the spectral synthesis code SME. [...] We report Li, C, O, Na, Mg, Al, Si, K, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Y, as well as Ba & we note that we employ non-LTE calculations for Li, O, Al, & Fe. We show that the use of astrometric & photometric data improves the accuracy of the derived spectroscopic parameters, especially $\log g$. [...] we recover the result that stars of the high-$\alpha$ sequence are typically older than stars in the low-$\alpha$ sequence, the latter spanning $-0.7<$[Fe/H]$<+0.5$. While these two sequences become indistinguishable in [$\alpha$/Fe] vs. [Fe/H] at the metal-rich regime, we find that age can be used to separate stars from the extended high-$\alpha$ & the low-$\alpha$ sequence even in this regime. [...] we find that the old stars ($>8$ Gyr have lower angular momenta $L_z$ than the Sun, which implies that they are on eccentric orbits & originate from the inner disk. Contrary to some previous smaller scale studies we find a continuous evolution in the high-$\alpha$-sequence up to super-solar [Fe/H] rather than a gap, which has been interpreted as a separate "high-$\alpha$ metal-rich" population. Stars in our sample that are younger than 10 Gyr, are mainly found on the low $\alpha$-sequence & show a gradient in $L_z$ from low [Fe/H] ($L_z>L_{z,\odot}$) towards higher [Fe/H] ($L_z<L_{z,\odot}$), which implies that the stars at the ends of this sequence are likely not originating from the close solar vicinity.

Swan bands - characteristic molecular absorption features of the C$_2$ molecule - are a spectroscopic signature of carbon-enhanced stars. They can also be used to identify carbon-enhanced metal-poor (CEMP) stars. The GALAH (GALactic Archaeology with Hermes) is a magnitude-limited survey of stars producing high-resolution, high signal-to-noise spectra. We used 627,708 GALAH spectra to search for carbon-enhanced stars with a supervised and unsupervised classification algorithm, relying on the imprint of the Swan bands. We identified 918 carbon-enhanced stars, including 12 already described in the literature. An unbiased selection function of the GALAH survey allows us to perform a population study of carbon-enhanced stars. Most of them are giants, out of which we find 28 CEMP candidates. A large fraction of our carbon-enhanced stars with repeated observations show variation in radial velocity, hinting that there is a large fraction of variables among them. 32 of the detected stars also show strong Lithium enhancement in their spectra.

If the Galaxy is axisymmetric and in dynamical equilibrium, we expect negligible fluctuations in the residual line-of-sight velocity field. Recent results using the \apg{} survey find significant fluctuations in velocity for stars in the midplane ($|z|<$0.25 kpc) out to 5 kpc, suggesting that the dynamical influence of non-axisymmetric features i.e., the Milky Way's bar, spiral arms and merger events extends out to the Solar neighborhood. Their measured power spectrum has a characteristic amplitude of 11 \kms{} on a scale of 2.5 kpc. The existence of such large-scale streaming motions has important implications for determining the Sun's motion about the Galactic Centre. Using Red Clump stars from \glh{} and \apg{}, we map the line-of-sight velocities around the Sun (d$<$5 kpc), and $|z|<$1.25 kpc from the midplane. By subtracting a smooth axisymmetric model for the velocity field, we study the residual fluctuations and compare our findings with mock survey generated by \glx{}. We find negligible large-scale fluctuations away from the plane. In the mid-plane, we reproduce the earlier \apg{} power spectrum but with 20\% smaller amplitude (9.3 \kms{}) after taking into account a few systematics (e.g., volume completeness). Using a flexible axisymmetric model the power-amplitude is further reduced to 6.3 \kms{}. Additionally, our simulations show that, in the plane, distances are underestimated for high-mass Red Clump stars which can lead to spurious power-amplitude of about 5.2 \kms{}. Taking this into account, we estimate the amplitude of real fluctuations to be $<$4.6 \kms{}, about a factor of three less than the \apg{} result.

GALAH is a large-scale magnitude-limited southern stellar spectroscopic survey. Its second data release (GALAH DR2) provides values of stellar parameters and abundances of 23 elements for 342,682 stars (Buder et al.). Here we add a description of the public release of radial velocities with a typical accuracy of 0.1 km/s for 336,215 of these stars, achievable due to the large wavelength coverage, high resolving power and good signal to noise ratio of the observed spectra, but also because convective motions in stellar atmosphere and gravitational redshift from the star to the observer are taken into account. In the process we derive medians of observed spectra which are nearly noiseless, as they are obtained from between 100 and 1116 observed spectra belonging to the same bin with a width of 50 K in temperature, 0.2 dex in gravity, and 0.1 dex in metallicity. Publicly released 1181 median spectra have a resolving power of 28,000 and trace the well-populated stellar types with metallicities between -0.6 and +0.3. Note that radial velocities from GALAH are an excellent match to the accuracy of velocity components along the sky plane derived by Gaia for the same stars. The level of accuracy achieved here is adequate for studies of dynamics within stellar clusters, associations and streams in the Galaxy. So it may be relevant for studies of the distribution of dark matter.

Open cluster members are coeval and share the same initial bulk chemical compositions. Consequently, differences in surface abundances between members of a cluster that are at different evolutionary stages can be used to study the effects of mixing and internal chemical processing. We carry out an abundance analysis of seven elements (Li, O, Na, Mg, Al, Si, Fe) in 66 stars belonging to the open cluster M67, based on high resolution GALAH spectra, 1D MARCS model atmospheres, and, for the first time for a large sample of M67 stars, non-local thermodynamic equilibrium (non-LTE) radiative transfer. From the non-LTE analysis, we find a typical star-to-star scatter in the abundance ratios of around 0.05 dex; this scatter is slightly but systematically larger when LTE is assumed instead. We find trends in the abundance ratios with effective temperature, indicating systematic differences in the surface abundances between turn-off and giant stars; these trends are more pronounced when LTE is assumed. However, in the non-LTE analysis, most of the element trends have been flattened. Two species are exceptions to this behaviour, namely Al and Si, which both clearly display remaining trends in the non-LTE analysis. We comment on the possible origin of these trends, by comparing them with recent stellar models that include atomic diffusion.

Mira variables are useful distance indicators, due to their high luminosities and well-defined period-luminosity relation. We select 1863 Miras from SAAO and MACHO observations to examine their use as distance estimators in the Milky Way. We measure a distance to the Galactic centre of $R_0 = 7.9 \pm 0.3$ kpc, which is in good agreement with other literature values. The uncertainty has two components of $\sim$0.2 kpc each: the first is from our analysis and predominantly due to interstellar extinction, the second is due to zero-point uncertainties extrinsic to our investigation, such as the distance to the Large Magellanic Cloud (LMC). In an attempt to improve existing period-luminosity calibrations, we use theoretical models of Miras to determine the dependence of the period-luminosity relation on age, metallicity, and helium abundance, under the assumption that Miras trace the bulk stellar population. We find that at a fixed period of $\log P = 2.4$, changes in the predicted $K_s$ magnitudes can be approximated by $\Delta M_{Ks} \approx -0.109(\Delta \rm{[Fe/H]}) + 0.033( {\Delta}t/\rm{Gyr}) + 0.021 ({\Delta}Y/0.01)$, and these coefficients are nearly independent of period. The expected overestimate in the Galactic centre distance from using an LMC-calibrated relation is $\sim$0.3 kpc. This prediction is not validated by our analysis; a few possible reasons are discussed. We separately show that while the predicted color-color diagrams of solar-neighbourhood Miras work well in the near-infrared, though there are offsets from the model predictions in the optical and mid-infrared.

We present a study using the second data release of the GALAH survey of stellar parameters and elemental abundances of 15 pairs of stars identified by Oh et al 2017. They identified these pairs as potentially co-moving pairs using proper motions and parallaxes from Gaia DR1. We find that 11 very wide (>1.7 pc) pairs of stars do in fact have similar Galactic orbits, while a further four claimed co-moving pairs are not truly co-orbiting. Eight of the 11 co-orbiting pairs have reliable stellar parameters and abundances, and we find that three of those are quite similar in their abundance patterns, while five have significant [Fe/H] differences. For the latter, this indicates that they could be co-orbiting because of the general dynamical coldness of the thin disc, or perhaps resonances induced by the Galaxy, rather than a shared formation site. Stars such as these, wide binaries, debris of past star formation episodes, and coincidental co-orbiters, are crucial for exploring the limits of chemical tagging in the Milky Way.

We report the serendipitous observations of 571 luminous supergiants in the Magellanic Clouds by the spectroscopic GALAH and TESS-HERMES surveys: 434 stars in the Large Magellanic Cloud and 137 in the Small Magellanic Cloud. We also find one star that appears associated with structured star formation in the Magellanic Bridge. Both of these surveys are aimed at the local volume of the Galaxy but have simple, magnitude-limited selection functions that mean they include some observations of luminous extra-Galactic stars. The surveys determine stellar parameter and abundances using The Cannon, a data-driven generative modelling approach. In this work, we explore the results from The Cannon when it is fed the spectra of these intrinsically luminous supergiants in the Magellanic Clouds, which are well outside the normal bounds of The Cannon's training set. We find that, although the parameters are astrophysically incorrect, the $v\sin i$ and the abundances of lithium, barium, and magnesium are excellent discriminants of these stars. It shows that in the future, with an expanded training set, it should be possible to determine accurate values for these types of stars.

The Galactic Archaeology with HERMES (GALAH) survey is a large-scale stellar spectroscopic survey of the Milky Way and designed to deliver chemical information complementary to a large number of stars covered by the $Gaia$ mission. We present the GALAH second public data release (GALAH DR2) containing 342,682 stars. For these stars, the GALAH collaboration provides stellar parameters and abundances for up to 23 elements to the community. Here we present the target selection, observation, data reduction and detailed explanation of how the spectra were analysed to estimate stellar parameters and element abundances. For the stellar analysis, we have used a multi-step approach. We use the physics-driven spectrum synthesis of Spectroscopy Made Easy (SME) to derive stellar labels ($T_\mathrm{eff}$, $\log g$, $\mathrm{[Fe/H]}$, $\mathrm{[X/Fe]}$, $v_\mathrm{mic}$, $v \sin i$, $A_{K_S}$) for a representative training set of stars. This information is then propagated to the whole survey with the data-driven method of $The~Cannon$. Special care has been exercised in the spectral synthesis to only consider spectral lines that have reliable atomic input data and are little affected by blending lines. Departures from local thermodynamic equilibrium (LTE) are considered for several key elements, including Li, O, Na, Mg, Al, Si, and Fe, using 1D MARCS stellar atmosphere models. Validation tests including repeat observations, Gaia benchmark stars, open and globular clusters, and K2 asteroseismic targets lend confidence in our methods and results. Combining the GALAH DR2 catalogue with the kinematic information from $Gaia$ will enable a wide range of Galactic Archaeology studies, with unprecedented detail, dimensionality, and scope.

The primary goal of Galactic archaeology is to learn about the origin of the Milky Way from the detailed chemistry and kinematics of millions of stars. Wide-field multi-fibre spectrographs are increasingly used to obtain spectral information for huge samples of stars. Some surveys (e.g. GALAH) are attempting to measure up to 30 separate elements per star. Stellar abundance spectroscopy is a subtle art that requires a very high degree of spectral uniformity across each of the fibres. However wide-field spectrographs are notoriously non-uniform due to the fast output optics necessary to image many fibre outputs onto the detector. We show that precise spectroscopy is possible with such instruments across all fibres by employing a photonic comb -- a device that produces uniformly spaced spots of light on the CCD to precisely map complex aberrations. Aberrations are parametrized by a set of orthogonal moments with $\sim100$ independent parameters. We then reproduce the observed image by convolving high resolution spectral templates with measured aberrations as opposed to extracting the spectra from the observed image. Such a forward modeling approach also trivializes some spectroscopic reduction problems like fibre cross-talk, and reliably extracts spectra with a resolution $\sim2.3$ times above the nominal resolution of the instrument. Our rigorous treatment of optical aberrations also encourages a less conservative spectrograph design in the future.

Using GALAH survey data of nearby stars, we look at how structure in the planar (u,v) velocity distribution depends on metallicity and on viewing direction within the Galaxy. In nearby stars, with distance d < 1 kpc, the Hercules stream is most strongly seen in higher metallicity stars [Fe/H] > 0.2. The Hercules stream peak v value depends on viewed galactic longitude, which we interpret as due to the gap between the stellar stream and more circular orbits being associated with a specific angular momentum value of about 1640 km/s kpc. The association of the gap with a particular angular momentum value supports a bar resonant model for the Hercules stream. Moving groups previously identified in Hipparcos observations are easiest to see in stars nearer than 250 pc, and their visibility and peak velocities in the velocity distributions depends on both viewing direction (galactic longitude and hemisphere) and metallicity. We infer that there is fine structure in local velocity distributions that varies over distances of a few hundred pc in the Galaxy.

We report the discovery of a giant planet in event UKIRT-2017-BLG-001, detected by the UKIRT microlensing survey. The mass ratio between the planet and its host is $q=1.50_{-0.14}^{+0.17}\times10^{-3}$, about 1.5 times the Jupiter/Sun mass ratio. The event lies 0.35$^{\circ}$ from the Galactic center and suffers from high extinction of $A_K=1.68$. Therefore, it could be detected only by a near-infrared survey. The field also suffers from large spatial differential extinction, which makes it difficult to estimate the source properties required to derive the angular Einstein radius. Nevertheless, we find evidence suggesting that the source is located in the far disk. If correct, this would be the first source star of a microlensing event to be identified as belonging to the far disk. We estimate the lens mass and distance using a Bayesian analysis to find that the planet's mass is $1.28^{+0.37}_{-0.44}\,M_{J}$, and it orbits a $0.81^{+0.21}_{-0.27}\,M_{\odot}$ star at an instantaneous projected separation of $4.18^{+0.96}_{-0.88}$ AU. The system is at a distance of $6.3^{+1.6}_{-2.1}$ kpc, and so likely resides in the Galactic bulge. In addition, we find a non-standard extinction curve in this field, in agreement with previous results toward high-extinction fields near the Galactic center.

The Wide Field Infrared Survey Telescope (WFIRST) was the top ranked large space mission in the 2010 New Worlds, New Horizons decadal survey, and it was formed by merging the science programs of 3 different mission concepts, including the Microlensing Planet Finder (MPF) concept (Bennett \etal\ 2010). The WFIRST science program (Spergel \etal\ 2015) consists of a general observer program, a wavefront controlled technology program, and two targeted science programs: a program to study dark energy, and a statistical census of exoplanets with a microlensing survey, which uses nearly one quarter of WFIRST's observing time in the current design reference mission. The New Worlds, New Horizons (decadal survey) midterm assessment summarizes the science case for the WFIRST exoplanet microlensing survey with this statement: "WFIRST's microlensing census of planets beyond 1 AU will perfectly complement Kepler's census of compact systems, and WFIRST will also be able to detect free-floating planets unbound from their parent stars\rlap."

Microlensing can access planet populations that no other method can probe: cold wide-orbit planets beyond the snow line, planets in both the Galactic bulge and disk, and free floating planets (FFPs). The demographics of each population will provide unique constraints on planet formation. Over the past 5 years, U.S. microlensing campaigns with Spitzer and UKIRT have provided a powerful complement to international ground-based microlensing surveys, with major breakthroughs in parallax measurements and probing new regions of the Galaxy. The scientific vitality of these projects has also promoted the development of the U.S. microlensing community. In the 2020s, the U.S. can continue to play a major role in ground-based microlensing by leveraging U.S. assets to complement ongoing ground-based international surveys. LSST and UKIRT microlensing surveys would probe vast regions of the Galaxy, where planets form under drastically different conditions. Moreover, while ground-based surveys will measure the planet mass-ratio function beyond the snow line, adaptive optics (AO) observations with ELTs would turn all of these mass ratios into masses and also distinguish between very wide-orbit planets and genuine FFPs. To the extent possible, cooperation of U.S. scientists with international surveys should also be encouraged and supported.

Using data from the GALAH pilot survey, we determine properties of the Galactic thin and thick disks near the solar neighbourhood. The data cover a small range of Galactocentric radius ($7.9 \leq R_\mathrm{GC} \leq 9.5$ kpc), but extend up to 4 kpc in height from the Galactic plane, and several kpc in the direction of Galactic anti-rotation (at longitude $260 ^\circ \leq \ell \leq 280^\circ$). This allows us to reliably measure the vertical density and abundance profiles of the chemically and kinematically defined `thick' and `thin' disks of the Galaxy. The thin disk (low-$\alpha$ population) exhibits a steep negative vertical metallicity gradient, at d[M/H]/d$z=-0.18 \pm 0.01$ dex kpc$^{-1}$, which is broadly consistent with previous studies. In contrast, its vertical $\alpha$-abundance profile is almost flat, with a gradient of d[$\alpha$/M]/d$z$ = $0.008 \pm 0.002$ dex kpc$^{-1}$. The steep vertical metallicity gradient of the low-$\alpha$ population is in agreement with models where radial migration has a major role in the evolution of the thin disk. The thick disk (high-$\alpha$ population) has a weaker vertical metallicity gradient d[M/H]/d$z = -0.058 \pm 0.003$ dex kpc$^{-1}$. The $\alpha$-abundance of the thick disk is nearly constant with height, d[$\alpha$/M]/d$z$ = $0.007 \pm 0.002$ dex kpc$^{-1}$. The negative gradient in metallicity and the small gradient in [$\alpha$/M] indicate that the high-$\alpha$ population experienced a settling phase, but also formed prior to the onset of major SNIa enrichment. We explore the implications of the distinct $\alpha$-enrichments and narrow [$\alpha$/M] range of the sub-populations in the context of thick disk formation.

Detailed elemental abundance patterns of metal-poor ([Fe/H] ~ -1~dex) stars in the Galactic bulge indicate that a number of them are consistent with globular cluster (GC) stars and may be former members of dissolved GCs. This would indicate that a few per cent of the Galactic bulge was built up from destruction and/or evaporation of globular clusters. Here an attempt is made to identify such presumptive stripped stars originating from the massive, inner Galaxy globular cluster NGC~6441 using its rich RR Lyrae variable star (RRL) population. We present radial velocities of forty RRLs centered on the globular cluster NGC~6441. All of the 13 RRLs observed within the cluster tidal radius have velocities consistent with cluster membership, with an average radial velocity of 24 +- 5~km/s and a star-to-star scatter of 11~km/s. This includes two new RRLs that were previously not associated with the cluster. Eight RRLs with radial velocities consistent with cluster membership but up to three time the distance from the tidal radius are also reported. These potential extra-tidal RRLs also have exceptionally long periods, which is a curious characteristic of the NGC~6441 RRL population that hosts RRLs with periods longer than seen anywhere else in the Milky Way. As expected of stripped cluster stars, most are inline with the cluster's orbit. Therefore, either the tidal radius of NGC~6441 is underestimated and/or we are seeing dissolving cluster stars stemming from NGC~6441 that are building up the old spheroidal bulge.

The technique of chemical tagging uses the elemental abundances of stellar atmospheres to `reconstruct' chemically homogeneous star clusters that have long since dispersed. The GALAH spectroscopic survey --which aims to observe one million stars using the Anglo-Australian Telescope -- allows us to measure up to 30 elements or dimensions in the stellar chemical abundance space, many of which are not independent. How to find clustering reliably in a noisy high-dimensional space is a difficult problem that remains largely unsolved. Here we explore t-distributed stochastic neighbour embedding (t-SNE) -- which identifies an optimal mapping of a high-dimensional space into fewer dimensions -- whilst conserving the original clustering information. Typically, the projection is made to a 2D space to aid recognition of clusters by eye. We show that this method is a reliable tool for chemical tagging because it can: (i) resolve clustering in chemical space alone, (ii) recover known open and globular clusters with high efficiency and low contamination, and (iii) relate field stars to known clusters. t-SNE also provides a useful visualization of a high-dimensional space. We demonstrate the method on a dataset of 13 abundances measured in the spectra of 187,000 stars by the GALAH survey. We recover 7 of the 9 observed clusters (6 globular and 3 open clusters) in chemical space with minimal contamination from field stars and low numbers of outliers. With chemical tagging, we also identify two Pleiades supercluster members (which we confirm kinematically), one as far as 6$^\circ$ -- one tidal radius away from the cluster centre.

The assembly of the Milky Way bulge is an old topic in astronomy, one now in a period of renewed and rapid development. The dominant scenario for bulge formation is that of the Milky Way as a nearly pure disk galaxy, with the inner disk having formed a bar and buckled. This can potentially explain virtually all bulge stars with [Fe/H] <~ -1.0, comprising 95% of the stellar population. The evidence is the incredible success in N-body models of this type in making meaningful predictions, such as the rotation curve and velocity dispersion measured from radial velocities, and the spatial morphologies of the peanut/X-shape and the long bar. The classical bulge scenario remains viable for stars with [Fe/H] <~ -1.0 and potentially a minority of the other stars. A classical bulge is expected from Lambda-CDM cosmological simulations, can accentuate the properties of an existing bar in a hybrid system, and is most consistent with the bulge abundance trends such as [Mg/Fe], which are elevated relative to both the thin and thick disks. Finally, the clumpy-galaxy scenario is considered, as it is the correct description of most Milky Way precursors given observations of high-redshift galaxies. Simulations predict that these star-forming clumps will sometimes migrate to the centres of galaxies where they may form a bulge, and galaxies often include a bulge clump as well. They will possibly form a bar with properties consistent with those of the Milky Way, such as the exponential profile and metallicity gradient. Given the relative successes of these scenarios, the Milky Way bulge is plausibly of composite origin, with a classical bulge and/or inner halo numerically dominant for stars with [Fe/H] <~ -1.0, a buckling thick disk for stars with -1.0 <~ [Fe/H] <~ -0.50 perhaps descended from the clumpy galaxy phase, and a buckling thin disk for stars with [Fe/H] >~ -0.50.

Our Galaxy hosts the annihilation of a few $\times 10^{43}$ low-energy positrons every second. Radioactive isotopes capable of supplying such positrons are synthesised in stars, stellar remnants, and supernovae. For decades, however, there has been no positive identification of a main stellar positron source leading to suggestions that many positrons originate from exotic sources like the Galaxy's central super-massive black hole or dark matter annihilation. %, but such sources would not explain the recently-detected positron signal from the extended Galactic disk. Here we show that a single type of transient source, deriving from stellar populations of age 3-6 Gyr and yielding ~0.03 $M_\odot$ of the positron emitter $^{44}$Ti, can simultaneously explain the strength and morphology of the Galactic positron annihilation signal and the solar system abundance of the $^{44}$Ti decay product $^{44}$Ca. This transient is likely the merger of two low-mass white dwarfs, observed in external galaxies as the sub-luminous, thermonuclear supernova known as SN1991bg-like.

Galah is an ongoing high-resolution spectroscopic survey with the goal of disentangling the formation history of the Milky Way, using the fossil remnants of disrupted star formation sites which are now dispersed around the Galaxy. It is targeting a randomly selected, magnitude limited ($V \leq 14$) sample of stars, with the goal of observing one million objects. To date, 300,000 spectra have been obtained. Not all of them are correctly processed by parameter estimation pipelines and we need to know about them. We present a semi-automated classification scheme which identifies different types of peculiar spectral morphologies, in an effort to discover and flag potentially problematic spectra and thus help to preserve the integrity of the survey's results. To this end we employ a recently developed dimensionality reduction technique t-SNE (t-distributed Stochastic Neighbour Embedding), which enables us to represent the complex spectral morphology in a two-dimensional projection map while still preserving the properties of the local neighbourhoods of spectra. We find that the majority (178,483) of the 209,533 Galah spectra considered in this study represents normal single stars, whereas 31,050 peculiar and problematic spectra with very diverse spectral features pertaining to 28,579 stars are distributed into 10 classification categories: Hot stars, Cool metal-poor giants, Molecular absorption bands, Binary stars, H$\alpha$/H$\beta$ emission, H$\alpha$/H$\beta$ emission superimposed on absorption, H$\alpha$/H$\beta$ P-Cygni, H$\alpha$/H$\beta$ inverted P-Cygni, Lithium absorption, and Problematic. Classified spectra with supplementary information are presented in the catalogue, indicating candidates for follow-up observations and population studies of the short-lived phases of stellar evolution.

The Galactic Archaeology with HERMES (GALAH) Survey is a massive observational project to trace the Milky Way's history of star formation, chemical enrichment, stellar migration and minor mergers. Using high-resolution (R$\simeq$28,000) spectra taken with the High Efficiency and Resolution Multi-Element Spectrograph (HERMES) instrument at the Anglo-Australian Telescope (AAT), GALAH will determine stellar parameters and abundances of up to 29 elements for up to one million stars. Selecting targets from a colour-unbiased catalogue built from 2MASS, APASS and UCAC4 data, we expect to observe dwarfs at 0.3 to 3 kpc and giants at 1 to 10 kpc. This enables a thorough local chemical inventory of the Galactic thin and thick disks, and also captures smaller samples of the bulge and halo. In this paper we present the plan, process and progress as of early 2016 for GALAH survey observations. In our first two years of survey observing we have accumulated the largest high-quality spectroscopic data set at this resolution, over 200,000 stars. We also present the first public GALAH data catalogue: stellar parameters (Teff, log(g), [Fe/H], [alpha/Fe]), radial velocity, distance modulus and reddening for 10680 observations of 9860 Tycho-2 stars that may be included in the first Gaia data release.

We present the data reduction procedures being used by the GALAH survey, carried out with the HERMES fibre-fed, multi-object spectrograph on the 3.9~m Anglo-Australian Telescope. GALAH is a unique survey, targeting 1 million stars brighter than magnitude V=14 at a resolution of 28,000 with a goal to measure the abundances of 29 elements. Such a large number of high resolution spectra necessitates the development of a reduction pipeline optimized for speed, accuracy, and consistency. We outline the design and structure of the Iraf-based reduction pipeline that we developed, specifically for GALAH, to produce fully calibrated spectra aimed for subsequent stellar atmospheric parameter estimation. The pipeline takes advantage of existing Iraf routines and other readily available software so as to be simple to maintain, testable and reliable. A radial velocity and stellar atmospheric parameter estimator code is also presented, which is used for further data analysis and yields a useful verification of the reduction quality. We have used this estimator to quantify the data quality of GALAH for fibre cross-talk level ($\lesssim0.5$%) and scattered light ($\sim5$ counts in a typical 20 minutes exposure), resolution across the field, sky spectrum properties, wavelength solution reliability (better than $1$ $\mathrm{km\ s^{-1}}$ accuracy) and radial velocity precision.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

Cosmological models predict the oldest stars in the Galaxy should be found closest to the centre of the potential well, in the bulge. The EMBLA Survey successfully searched for these old, metal-poor stars by making use of the distinctive SkyMapper photometric filters to discover candidate metal-poor stars in the bulge. Their metal-poor nature was then confirmed using the AAOmega spectrograph on the AAT. Here we present an abundance analysis of 10 bulge stars with -2.8<[Fe/H]<-1.7 from MIKE/Magellan observations, in total determining the abundances of 22 elements. Combining these results with our previous high-resolution data taken as part of the Gaia-ESO Survey, we have started to put together a picture of the chemical and kinematic nature of the most metal-poor stars in the bulge. The currently available kinematic data is consistent with the stars belonging to the bulge, although more accurate measurements are needed to constrain the stars' orbits. The chemistry of these bulge stars deviates from that found in halo stars of the same metallicity. Two notable differences are the absence of carbon-enhanced metal-poor bulge stars, and the alpha-element abundances exhibit a large intrinsic scatter and include stars which are underabundant in these typically enhanced elements.

We present an analysis of the relative throughputs of the 3.9-metre Anglo-Australian Telescope's 2dF/HERMES system, based upon spectra acquired during the first two years of the GALAH survey. Averaged spectral fluxes of stars were compared to their photometry to determine the relative throughputs of fibres for a range of fibre position and atmospheric conditions. We find that overall the throughputs of the 771 usable fibres have been stable over the first two years of its operation. About 2.5 per cent of fibres have throughputs much lower than the average. There are also a number of yet unexplained variations between the HERMES bandpasses, and mechanically & optically linked fibre groups known as retractors or slitlets related to regions of the focal plane. These findings do not impact the science that HERMES will produce.

I review the literature covering the issue of interstellar extinction toward the Milky Way bulge, with emphasis placed on findings from planetary nebulae, RR Lyrae, and red clump stars. I also report on observations from HI gas and globular clusters. I show that there has been substantial progress in this field in recent decades, most particularly from red clump stars. The spatial coverage of extinction maps has increased by a factor $\sim 100 \times$ in the past twenty years, and the total-to-selective extinction ratios reported have shifted by $\sim$20-25\%, indicative of the improved accuracy and separately, of a steeper-than-standard extinction curve. Problems remain in modelling differential extinction, explaining anomalies involving the planetary nebulae, and understanding the difference between bulge extinction coefficients and "standard" literature values.

We present 947 radial velocities of RR Lyrae variable stars in four fields located toward the Galactic bulge, observed within the data from the ongoing Bulge RR Lyrae Radial Velocity Assay (BRAVA-RR). We show that these RR Lyrae stars exhibit hot kinematics and null or negligible rotation and are therefore members of a separate population from the bar/pseudobulge that currently dominates the mass and luminosity of the inner Galaxy. Our RR Lyrae stars predate these structures, and have metallicities, kinematics, and spatial distribution that are consistent with a "classical" bulge, although we cannot yet completely rule out the possibility that they are the metal-poor tail of a more metal rich ([Fe/H] ~ -1 dex) halo-bulge population. The complete catalog of radial velocities for the BRAVA-RR stars is also published electronically.

$K2$'s Campaign 9 ($K2$C9) will conduct a $\sim$3.7 deg$^{2}$ survey toward the Galactic bulge from 7/April through 1/July of 2016 that will leverage the spatial separation between $K2$ and the Earth to facilitate measurement of the microlens parallax $\pi_{\rm E}$ for $\gtrsim$127 microlensing events. These will include several that are planetary in nature as well as many short-timescale microlensing events, which are potentially indicative of free-floating planets (FFPs). These satellite parallax measurements will in turn allow for the direct measurement of the masses of and distances to the lensing systems. In this white paper we provide an overview of the $K2$C9 space- and ground-based microlensing survey. Specifically, we detail the demographic questions that can be addressed by this program, including the frequency of FFPs and the Galactic distribution of exoplanets, the observational parameters of $K2$C9, and the array of resources dedicated to concurrent observations. Finally, we outline the avenues through which the larger community can become involved, and generally encourage participation in $K2$C9, which constitutes an important pathfinding mission and community exercise in anticipation of $WFIRST$.

We investigate interstellar extinction curve variations toward $\sim$4 deg$^{2}$ of the inner Milky Way in $VIJK_{s}$ photometry from the OGLE-III and $VVV$ surveys, with supporting evidence from diffuse interstellar bands and $F435W,F625W$ photometry. We obtain independent measurements toward $\sim$2,000 sightlines of $A_{I}$, $E(V-I)$, $E(I-J)$, and $E(J-K_{s})$, with median precision and accuracy of 2%. We find that the variations in the extinction ratios $A_{I}/E(V-I)$, $E(I-J)/E(V-I)$ and $E(J-K_{s})/E(V-I)$ are large (exceeding 20%), significant, and positively correlated, as expected. However, both the mean values and the trends in these extinction ratios are drastically shifted from the predictions of Cardelli and Fitzpatrick, regardless of how $R_{V}$ is varied. Furthermore, we demonstrate that variations in the shape of the extinction curve has at least two degrees of freedom, and not one (e.g. $R_{V}$), which we conform with a principal component analysis. We derive a median value of $<A_{V}/A_{Ks}>=13.44$, which is $\sim$60% higher than the "standard" value. We show that the Wesenheit magnitude $W_{I}=I-1.61(I-J)$ is relatively impervious to extinction curve variations. Given that these extinction curves are linchpins of observational cosmology, and that it is generally assumed that $R_{V}$ variations correctly capture variations in the extinction curve, we argue that systematic errors in the distance ladder from studies of type Ia supernovae and Cepheids may have been underestimated. Moreover, the reddening maps from the Planck experiment are shown to systematically overestimate dust extinction by $\sim$100%, and lack sensitivity to extinction curve variations.

The first stars are predicted to have formed within 200 million years after the Big Bang, initiating the cosmic dawn. A true first star has not yet been discovered, although stars with tiny amounts of elements heavier than helium ('metals') have been found in the outer regions ('halo') of the Milky Way. The first stars and their immediate successors should, however, preferentially be found today in the central regions ('bulges') of galaxies, because they formed in the largest over-densities that grew gravitationally with time. The Milky Way bulge underwent a rapid chemical enrichment during the first 1-2 billion years, leading to a dearth of early, metal-poor stars. Here we report observations of extremely metal-poor stars in the Milky Way bulge, including one star with an iron abundance about 10,000 times lower than the solar value without noticeable carbon enhancement. We confirm that the most metal-poor bulge stars are on tight orbits around the Galactic Centre, rather than being halo stars passing through the bulge, as expected for stars formed at redshifts greater than 15. Their chemical compositions are in general similar to typical halo stars of the same metallicity although intriguing differences exist, including lower abundances of carbon.

For all exoplanet candidates, the reliability of a claimed detection needs to be assessed through a careful study of systematic errors in the data to minimize the false positives rate. We present a method to investigate such systematics in microlensing datasets using the microlensing event OGLE-2013-BLG-0446 as a case study. The event was observed from multiple sites around the world and its high magnification (A_{max} \sim 3000) allowed us to investigate the effects of terrestrial and annual parallax. Real-time modeling of the event while it was still ongoing suggested the presence of an extremely low-mass companion (\sim 3M_\oplus ) to the lensing star, leading to substantial follow-up coverage of the light curve. We test and compare different models for the light curve and conclude that the data do not favour the planetary interpretation when systematic errors are taken into account.

The stellar population of the Milky Way bulge is thoroughly studied, with a plethora of measurements from virtually the full suite of instruments available to astronomers. It is thus perhaps surprising that alongside well-established results lies some substantial uncertainty in its star-formation history. Cosmological models predict the bulge to host the Galaxy's oldest stars for [Fe/H]$\lesssim -1$, and this is demonstrated by RR Lyrae stars and globular cluster observations. There is consensus that bulge stars with [Fe/H]$\lesssim0$ are older than $t \approx10$ Gyr. However, at super-solar metallicity, there is a substantial unresolved discrepancy. Data from spectroscopic measurements of the main-sequence turnoff and subgiant branch, the abundances of asymptotic giant branch stars, the period distribution of Mira variables, the chemistry and central-star masses of planetary nebulae, all suggest a substantial intermediate-age population ($t \approx 3$ Gyr). This is in conflict with predictions from cosmologically-motivated chemical evolution models and photometric studies of the main-sequence turnoff region, which both suggest virtually no stars younger than $t \approx 8$ Gyr. A possible resolution to this conflict is enhanced helium-enrichment, as this would shift nearly all of the age estimates in the direction of decreasing discrepancy. Enhanced helium-enrichment is also arguably suggested by measurements of the red giant branch bump and the R-parameter.

We report on the RR Lyrae variable star, MACHO 176.18833.411, located toward the Galactic bulge and observed within the data from the ongoing Bulge RR Lyrae Radial Velocity Assay (BRAVA-RR), which has the unusual radial velocity of -372 +- 8 km/s and true space velocity of -482 +- 22 km/s relative to the Galactic rest frame. Located less than 1 kpc from the Galactic center and toward a field at (l,b)=(3,-2.5), this pulsating star has properties suggesting it belongs to the bulge RR Lyrae star population yet a velocity indicating it is abnormal, at least with respect to bulge giants and red clump stars. We show that this star is most likely a halo interloper and therefore suggest that halo contamination is not insignificant when studying metal-poor stars found within the bulge area, even for stars within 1 kpc of the Galactic center. We discuss the possibility that MACHO 176.18833.411 is on the extreme edge of the bulge RR Lyrae radial velocity distribution, and also consider a more exotic scenario in which it is a runaway star moving through the Galaxy.

Space-based microlens parallax measurements are a powerful tool for understanding planet populations, especially their distribution throughout the Galaxy. However, if space-based observations of the microlensing events must be specifically targeted, it is crucial that microlensing events enter the parallax sample without reference to the known presence or absence of planets. Hence, it is vital to define objective criteria for selecting events where possible and to carefully consider and minimize the selection biases where not possible so that the final sample represents a controlled experiment. We present objective criteria for initiating observations and determining their cadence for a subset of events, and we define procedures for isolating subjective decision making from information about detected planets for the remainder of events. We also define procedures to resolve conflicts between subjective and objective selections. These procedures maximize planet sensitivity of the sample as a whole by allowing for planet detections even if they occur before satellite observations for objectively-selected events and by helping to trigger fruitful follow-up observations for subjectively-chosen events. This paper represents our public commitment to these procedures, which is a necessary component of enforcing objectivity on the experimental protocol.

The GALAH survey is a large high-resolution spectroscopic survey using the newly commissioned HERMES spectrograph on the Anglo-Australian Telescope. The HERMES spectrograph provides high-resolution (R ~28,000) spectra in four passbands for 392 stars simultaneously over a 2 degree field of view. The goal of the survey is to unravel the formation and evolutionary history of the Milky Way, using fossil remnants of ancient star formation events which have been disrupted and are now dispersed throughout the Galaxy. Chemical tagging seeks to identify such dispersed remnants solely from their common and unique chemical signatures; these groups are unidentifiable from their spatial, photometric or kinematic properties. To carry out chemical tagging, the GALAH survey will acquire spectra for a million stars down to V~14. The HERMES spectra of FGK stars contain absorption lines from 29 elements including light proton-capture elements, alpha-elements, odd-Z elements, iron-peak elements and n-capture elements from the light and heavy s-process and the r-process. This paper describes the motivation and planned execution of the GALAH survey, and presents some results on the first-light performance of HERMES.

I revisit the Cepheid-distance determination to the nearby spiral galaxy M101 (Pinwheel Galaxy) of Shappee & Stanek (2011), in light of several recent investigations questioning the shape of the interstellar extinction curve at $\lambda \approx 8,000$ \AA (i.e. I-band). I find that the relatively steep extinction ratio $A_{I}/E(V-I)=1.1450$ (Fitzpatrick & Massa 2007) is slightly favoured relative to $A_{I}/E(V-I)=1.2899$ (Fitzpatrick 1999) and significantly favoured relative the historically canonical value of $A_{I}/E(V-I)=1.4695$ (Cardelli et al. 1989). The steeper extinction curves, with lower values of $A_{I}/E(V-I)$, yield fits with reduced scatter, metallicity-dependences to the dereddened Cepheid luminosities that are closer to values inferred in the local group, and that are less sensitive to the choice of reddening cut imposed in the sample selection. The increase in distance modulus to M101 when using the preferred extinction curve is ${\Delta}{\mu} \sim 0.06$ mag, resulting in an estimate of the distance modulus to M101 relative to the LMC of $ {\Delta}\mu_{\rm{LMC}} \approx 10.72 \pm 0.03$ (stat). The best-fit metallicity-dependence is $dM_{I}/d\rm{[O/H]} \approx (-0.38 \pm 0.14$ (stat)) mag dex$^{-1}$.

We model the split red clump of the Galactic bulge in OGLE-III photometry, and compare the results to predictions from two N-body models. Our analysis yields precise maps of the brightness of the two red clumps, the fraction of stars in the more distant peak, and their combined surface density. We compare the observations to predictions from two N-body models previously used in the literature. Both models correctly predict several features as long as one assumes an angle $\alpha_{\rm{Bar}} \approx 30^{\circ}$ between the Galactic bar's major axis and the line of sight to the Galactic centre. In particular that the fraction of stars in the faint red clump should decrease with increasing longitude. The biggest discrepancies between models and data are in the rate of decline of the combined surface density of red clump stars toward negative longitudes and of the brightness difference between the two red clumps toward positive longitudes, with neither discrepancy exceeding $\sim$25% in amplitude. Our analysis of the red giant luminosity function also yields an estimate of the red giant branch bump parameters toward these high-latitude fields, and evidence for a high rate ($\sim$25%) of disk contamination in the bulge at the colour and magnitude of the red clump, with the disk contamination rate increasing toward sightlines further distant from the plane.

We compare model predictions to observations of star counts in the red giant branch bump (RGBB) relative to the number density of first-ascent red giant branch at the magnitude of the RGBB, $EW_{RGBB}$. The predictions are shown to exceed the data by $(5.2 \pm 4.3)$% for the BaSTI models and by $(17.1 \pm 4.3)$% for the Dartmouth models, where the listed errors are purely statistical. These two offsets are brought to zero if the Galactic globular cluster metallicity scale is assumed to be overestimated by a linear shift of $\sim 0.11$ dex and $\sim 0.36$ dex respectively. This inference based on RGBB star counts goes in the opposite direction to the increase in metallicities of ${\Delta}$[M/H]$\approx$0.20 dex that would be required to fix the offset between predicted and observed RGBB luminosities. This comparison is a constraint on "deep mixing" models of stellar interiors, which predict decreased rather than increased RGBB star counts. We tabulate the predictions for RGBB star counts as a function of [Fe/H], [$\alpha$/Fe], CNONa, initial helium abundance, and age. Though our study suggests a small zero-point calibration issue, RGBB star counts should nonetheless be an actionable parameter with which to constrain stellar populations in the differential sense. The most significant outliers are toward the clusters NGC 5025 (M53), NGC 6723, and NGC 7089 (M2), each of which shows a $\sim 2 \sigma$ deficit in their RGBB star counts.

We present the first results of the EMBLA survey (Extremely Metal-poor BuLge stars with AAOmega), aimed at finding metal-poor stars in the Milky Way bulge, where the oldest stars should now preferentially reside. EMBLA utilises SkyMapper photometry to pre-select metal-poor candidates, which are subsequently confirmed using AAOmega spectroscopy. We describe the discovery and analysis of four bulge giants with -2.72<=[Fe/H]<=-2.48, the lowest metallicity bulge stars studied with high-resolution spectroscopy to date. Using FLAMES/UVES spectra through the Gaia-ESO Survey we have derived abundances of twelve elements. Given the uncertainties, we find a chemical similarity between these bulge stars and halo stars of the same metallicity, although the abundance scatter may be larger, with some of the stars showing unusual [{\alpha}/Fe] ratios.

We demonstrate that failure to properly account for stellar evolution can bias results in determinations of the spatial morphology of Galactic bulge stars, focusing on the question of whether or not the X-shape is more pronounced among the more metal-rich stars than among the metal-poor stars. We argue that this trend, a result recently claimed by three separate groups, may have been overestimated as it is relatively easier to detect a bimodality in the distance distribution function at higher metallicities. This is due to three factors. First, the intrinsic colour of red clump and red giant stars vary with metallicity, at the level d(V-I)_{RC}/d\[M/H] ~ 0.25 mag dex^{-1}, and thus the ratio of red clump to red giant stars within a spectroscopic sample will depend on the photometric selection of any investigation. Second, the duration of ascent of the red giant branch goes down and the red clump lifetime goes up as metallicity increases, which has the effect of increasing the ratio of red clump to red giant stars by as much as ~33% over the range of the bulge metallicity-distribution function. Finally, over the same metallicity interval, the effective number of red giant branch bump stars is predicted to increase by ~200%, and their presence becomes degenerate with the observational parameters of the two red clumps, creating an illusory increase in signal-to-noise for a second peak in the distance modulus distribution.

The Korean Microlensing Telescope Network (KMTNet) will consist of three 1.6m telescopes each with a 4 deg^{2} field of view (FoV) and will be dedicated to monitoring the Galactic Bulge to detect exoplanets via gravitational microlensing. KMTNet's combination of aperture size, FoV, cadence, and longitudinal coverage will provide a unique opportunity to probe exoplanet demographics in an unbiased way. Here we present simulations that optimize the observing strategy for, and predict the planetary yields of, KMTNet. We find preferences for four target fields located in the central Bulge and an exposure time of t_{exp} = 120s, leading to the detection of ~2,200 microlensing events per year. We estimate the planet detection rates for planets with mass and separation across the ranges 0.1 <= M_{p}/M_{Earth} <= 1000 and 0.4 <= a/AU <= 16, respectively. Normalizing these rates to the cool-planet mass function of Cassan (2012), we predict KMTNet will be approximately uniformly sensitive to planets with mass 5 <= M_{p}/M_{Earth} <= 1000 and will detect ~20 planets per year per dex in mass across that range. For lower-mass planets with mass 0.1 <= M_{p}/M_{Earth} < 5, we predict KMTNet will detect ~10 planets per year. We also compute the yields KMTNet will obtain for free-floating planets (FFPs) and predict KMTNet will detect ~1 Earth-mass FFP per year, assuming an underlying population of one such planet per star in the Galaxy. Lastly, we investigate the dependence of these detection rates on the number of observatories, the photometric precision limit, and optimistic assumptions regarding seeing, throughput, and flux measurement uncertainties.

We investigate the effect of helium enrichment on the evolution and nucleosynthesis of low-mass asymptotic giant branch (AGB) stars of 1.7Msun and 2.36Msun with a metallicity of Z=0.0006 ([Fe/H] = -1.4). We calculate evolutionary sequences with the primordial helium abundance (Y = 0.24) and with helium-enriched compositions (Y = 0.30, 0.35, 0.40). For comparison we calculate models of the same mass but at a lower metallicity Z=0.0003 ([Fe/H] = -1.8) with Y=0.24. Post-processing nucleosynthesis calculations are performed on each of the evolutionary sequences to determine the production of elements from hydrogen through to bismuth. Elemental surface abundance predictions and stellar yields are presented for each model. The models with enriched helium have shorter main sequence and AGB lifetimes, and enter the AGB with a more massive hydrogen exhausted core than the primordial helium model. The main consequences are 1) low-mass AGB models with enhanced helium will evolve more than twice as fast, giving them the chance to contribute sooner to the chemical evolution of the forming globular clusters, and 2) the stellar yields will be strongly reduced relative to their primordial helium counterparts. An increase of Delta Y = 0.10 at a given mass decreases the yields of carbon by up to ~60%, of fluorine by up to 80%, and decreases the yields of the s-process elements barium and lanthanum by about 45%. While the yields of first s-process peak elements strontium, yttrium and zirconium decrease by up to 50%, the yields of rubidium either do not change or increase.

I discuss the age and helium abundance of the Galactic bulge stellar population. I present examples as to how age and helium abundance can be degenerate observationally, and thus, how unstated assumptions of the helium-metallicity relation can result in biased age-metallicity determinations. I summarise efforts to tackle this degeneracy using the red giant branch bump, forward modelling of microlensing selection effects, and direct mass and metallicity measurements of stars along the red giant branch. I make the case that the degeneracy between these two parameters can be resolved in the near future.

We present a study of the luminosity density distribution of the Galactic bar using number counts of red clump giants (RCGs) from the OGLE-III survey. The data were recently published by Nataf et al. (2013) for 9019 fields towards the bulge and have $2.94\times 10^6$ RC stars over a viewing area of $90.25 \,\textrm{deg}^2$. The data include the number counts, mean distance modulus ($\mu$), dispersion in $\mu$ and full error matrix, from which we fit the data with several tri-axial parametric models. We use the Markov Chain Monte Carlo (MCMC) method to explore the parameter space and find that the best-fit model is the $E_3$ model, with the distance to the GC is 8.13 kpc, the ratio of semi-major and semi-minor bar axis scale lengths in the Galactic plane $x_{0},y_{0}$, and vertical bar scale length $z_0$, is $x_0:y_0:z_0 \approx 1.00:0.43:0.40$ (close to being prolate). The scale length of the stellar density profile along the bar's major axis is $\sim$ 0.67 kpc and has an angle of $29.4^\circ$, slightly larger than the value obtained from a similar study based on OGLE-II data. The number of estimated RC stars within the field of view is $2.78 \times 10^6$, which is systematically lower than the observed value. We subtract the smooth parametric model from the observed counts and find that the residuals are consistent with the presence of an X-shaped structure in the Galactic centre, the excess to the estimated mass content is $\sim 5.8%$. We estimate the total mass of the bar is $\sim 1.8 \times 10^{10} M_\odot$. Our results can be used as a key ingredient to construct new density models of the Milky Way and will have implications on the predictions of the optical depth to gravitational microlensing and the patterns of hydrodynamical gas flow in the Milky Way.

We combine VI photometry from OGLE-III with VVV and 2MASS measurements of E(J-K_{s}) to resolve the longstanding problem of the non-standard optical extinction toward the Galactic bulge. We show that the extinction is well-fit by the relation A_{I} = 0.7465*E(V-I) + 1.3700*E(J-K_{s}), or, equivalently, A_{I} = 1.217*E(V-I)(1+1.126*(E(J-K_{s})/E(V-I)-0.3433)). The optical and near-IR reddening law toward the inner Galaxy approximately follows an R_{V} \approx 2.5 extinction curve with a dispersion {\sigma}_{R_{V}} \approx 0.2, consistent with extragalactic investigations of the hosts of type Ia SNe. Differential reddening is shown to be significant on scales as small as as our mean field size of 6', with the 1{\sigma} dispersion in reddening averaging 9% of total reddening for our fields. The intrinsic luminosity parameters of the Galactic bulge red clump (RC) are derived to be (M_{I,RC}, \sigma_{I,RC,0}, (V-I)_{RC,0}, \sigma_{(V-I)_{RC}}, (J-K_{s})_{RC,0}) = (-0.12, 0.09, 1.06, 0.121, 0.66). Our measurements of the RC brightness, brightness dispersion and number counts allow us to estimate several Galactic bulge structural parameters. We estimate a distance to the Galactic center of 8.20 kpc, resolving previous discrepancies in distance determinations to the bulge based on I-band observations. We measure an upper bound on the tilt {\alpha} \approx 40{\deg}. between the bar's major axis and the Sun-Galactic center line of sight, though our brightness peaks are consistent with predictions of an N-body model oriented at {\alpha} \approx 25{\deg}. The number of RC stars suggests a total stellar mass for the Galactic bulge of 2.0*10^{10} M_{\odot}, if one assumes a Salpeter IMF.

We present the broadest and most precise empirical investigation of red giant branch bump (RGBB) brightness and number counts ever conducted. We implement a new method and use data from two \textit{Hubble Space Telescope (HST)} globular cluster (GC) surveys to measure the brightness and star counts of the RGBB in 72 GCs. The brightness is measured to a precision better than 0.01 mag while the precision in number counts reaches 10%. The position of the main-sequence turnoff (MSTO) and the number of horizontal branch (HB) stars are used as comparisons where appropriate. Several independent scientific conclusions are newly possible with our parametrization of the RGBB. Both brightness and number counts are shown to have second parameters in addition to their strong dependence on metallicity. The RGBBs are found to be anomalous in the GCs NGC 2808, 5286, 6388 and 6441, likely due to the presence of multiple populations. Finally, we use our empirical calibration to predict the properties of the Galactic bulge RGBB if the assumption of similar stellar physics for the bulge and Galactic GC system holds. The RGBB properties for the bulge are shown to differ from those of the Galactic GC system, with the former having lower number counts, a lower brightness dispersion and a brighter peak luminosity than would be expected from the latter. This discrepancy is well explained by the Galactic bulge having a higher helium abundance than expected from GCs, ${\Delta}$Y$\sim +$0.06 at the median metallicity.

The presence of two stellar populations in the Milky Way bulge has been reported recently. We aim at studying the abundances and kinematics of stars in the outer bulge, thereby providing additional constraints on models of its formation. Spectra of 401 red giant stars in a field at (l,b)=(0{\deg},-10{\deg}) were obtained with FLAMES at the VLT. Stars of luminosities down to below the two bulge red clumps (RCs) are included. From these spectra we measure general metallicities, abundances of Fe and the alpha-elements, and radial velocities (RV) of the stars. These measurements as well as photometric data are compared to simulations with the Besancon and TRILEGAL models of the Galaxy. We confirm the presence of two populations among our sample stars: i) a metal-rich one at [M/H] ~+0.3, comprising about 30% of the sample, with low RV dispersion and low alpha-abundance, and ii) a metal-poor population at [M/H] ~-0.6 with high RV dispersion and high alpha-abundance. The metal-rich population could be connected to the Galactic bar. We identify this population as the carrier of the double RC feature. We do not find a significant difference in metallicity or RV between the two RCs, a small difference in metallicity being probably due to a selection effect. The RV dispersion agrees well with predictions of the Besancon Galaxy model, but the metallicity of the "thick bulge" model component should be shifted to lower metallicity by 0.2 to 0.3dex to well reproduce the observations. We present evidence that the metallicity distribution function depends on the evolutionary state of the sample stars, suggesting that enhanced mass loss preferentially removes metal-rich stars. We also confirm the decrease of \alpha-element over-abundance with increasing metallicity.

We identify 34 highly-probable detached, red giant eclipsing binary pairs among 315 candidates in Devor's catalog of $\sim$10,000 OGLE-II eclipsing binaries. We estimate that there should be at least 200 such systems in OGLE-III. We show that spectroscopic measurements of the metallicities and radial-velocity-derived masses of these systems would independently constrain both the age-metallicity and helium-metallicity relations of the Galactic Bulge, potentially breaking the age-helium degeneracy that currently limits our ability to characterize the Bulge stellar population. Mass and metallicity measurements alone would be sufficient to immediately validate or falsify recent claims about the age and helium abundance of the Bulge. A spectroscopic survey of these systems would constrain models of Milky Way assembly, as well as provide significant auxiliary science on research questions such as mass loss on the red giant branch. We discuss the theoretical uncertainties in stellar evolution models that would need to be accounted for to maximize the scientific yield.

We show that the factor $\sim$2 discrepancy between spectroscopic and photometric age determinations of the Galactic bulge main-sequence turnoff can be naturally explained by positing an elevated helium enrichment for the bulge relative to that assumed by standard isochrones. We obtain an upper bound on the helium enrichment parameter of the bulge $({\Delta}Y/{\Delta}Z)_{\rm{Bulge}} \lesssim 5.0$ given the requirement that the spectroscopic and photometric ages be consistent and the limiting condition of instantaneous star formation. The corresponding mean age for the bulge is $t_{\rm{Bulge}} \approx 10$ Gyr. We discuss phenomenological evidence that the bulge may have had a chemical evolution that is distinct from the solar neighborhood in this manner, and we make several testable predictions. Should this emerging picture of the bulge as helium-enhanced hold, it will require the development of new isochrones, new model atmospheres, and modified analysis and cosmological interpretation of the integrated light of other bulges and elliptical galaxies.

Microlensing detections of cool planets are important for the construction of an unbiased sample to estimate the frequency of planets beyond the snow line, which is where giant planets are thought to form according to the core accretion theory of planet formation. In this paper, we report the discovery of a giant planet detected from the analysis of the light curve of a high-magnification microlensing event MOA-2010-BLG-477. The measured planet-star mass ratio is $q=(2.181\pm0.004)\times 10^{-3}$ and the projected separation is $s=1.1228\pm0.0006$ in units of the Einstein radius. The angular Einstein radius is unusually large $\theta_{\rm E}=1.38\pm 0.11$ mas. Combining this measurement with constraints on the "microlens parallax" and the lens flux, we can only limit the host mass to the range $0.13<M/M_\odot<1.0$. In this particular case, the strong degeneracy between microlensing parallax and planet orbital motion prevents us from measuring more accurate host and planet masses. However, we find that adding Bayesian priors from two effects (Galactic model and Keplerian orbit) each independently favors the upper end of this mass range, yielding star and planet masses of $M_*=0.67^{+0.33}_{-0.13}\ M_\odot$ and $m_p=1.5^{+0.8}_{-0.3}\ M_{\rm JUP}$ at a distance of $D=2.3\pm0.6$ kpc, and with a semi-major axis of $a=2^{+3}_{-1}$ AU. Finally, we show that the lens mass can be determined from future high-resolution near-IR adaptive optics observations independently from two effects, photometric and astrometric.

We have analyzed the data on 16,836 RR Lyrae (RR Lyr) variables observed toward the Galactic bulge during the third phase of the Optical Gravitational Lensing Experiment (OGLE-III), which took place in 2001-2009. Using these standard candles, we show that the ratio of total to selective extinction toward the bulge is given by R_I=A_I/E(V-I)=1.080+/-0.007 and is independent of color. We demonstrate that the bulge RR Lyr stars form a metal-uniform population, slightly elongated in its inner part. The photometrically derived metallicity distribution is sharply peaked at [Fe/H]=-1.02+/-0.18, with a dispersion of 0.25 dex. In the inner regions (|l|<3, |b|<4) the RR Lyr tend to follow the barred distribution of the bulge red clump giants. The distance to the Milky Way center inferred from the bulge RR Lyr is R_0=8.54+/-0.42 kpc. We report a break in the mean density distribution at a distance of ~0.5 kpc from the center indicating its likely flattening. Using the OGLE-III data, we assess that (4-7)x10^4 type ab RR Lyr variables should be detected toward the bulge area of the on going near-IR VISTA Variables in the Via Lactea (VVV) survey, where the uncertainty partially results from the unknown RR Lyr spatial density distribution within 0.2 kpc from the Galactic center.

We measure the intrinsic magnitude dispersion of the Galactic bulge red clump (RC) using two different methods and arrive at an estimate of \sigma_{I} ~ 0.17 mag. We first estimate the width of the RC by analyzing in detail a sample of RC measurements toward double RC sightlines, which are toward regions of the bulge that are geometrically thin. We obtain \sigma_{I} ~ 0.17 mag. We then estimate the value by comparing the luminosity function for bulge RC and red giant (RG) stars in V to that in I for a sightline near the plane, and we obtain a value of \sigma_{I} ~ 0.14-0.17 mag. This result has structural and evolutionary repercussions. It constrains model estimates of the bar's orientation angle derived from RC studies as well as the inferred characteristic length of the bar for a given characteristic height. Moreover, the value is too large to account for using the predicted evolutionary effects of the known metallicity distribution, implying the presence of other population effects.

We combine ground and space-based photometry of the Galactic globular cluster 47 Tuc to measure four independent lines of evidence for a helium gradient in the cluster, whereby stars in the cluster outskirts would have a lower initial helium abundance than stars in and near the cluster core. First and second, we show that the red giant branch bump (RGBB) stars exhibit gradients in their number counts and brightness. With increased separation from the cluster center, they become more numerous relative to the other red giant (RG) stars. They also become fainter. For our third and fourth lines of evidence, we show that the horizontal branch (HB) of the cluster becomes both fainter and redder for sightlines farther from the cluster center. These four results are respectively detected at the 2.3$\sigma$, 3.6$\sigma$, 7.7$\sigma$ and 4.1$\sigma$ levels. Each of these independent lines of evidence is found to be significant in the cluster-outskirts; closer in, the data are more compatible with uniform mixing. Our radial profile is qualitatively consistent with but quantitatively tighter than previous results based on CN absorption. These observations are qualitatively consistent with a scenario wherein a second generation of stars with modestly enhanced helium and CNO abundance formed deep within the gravitational potential of a cluster of previous generation stars having more canonical abundances.

We measure the red giant branch bump (RGBB) of the Galactic bulge, the most metal-rich RGBB ever detected. The RGBB luminosity functions peaks at the expected brightness, but its number density is very low relative to Galactic globular cluster calibrations, implying the Galactic bulge has a higher helium enrichment parameter {\Delta}Y/{\Delta}Z $\ge$ 4.0 for Y~0.35 rather than the standard 2.0 with Y=0.27. The RGBB is (0.71 +/- 0.02) mag fainter than the red clump (RC) in I toward the densest stellar regions imaged by the OGLE-III Galactic bulge photometric survey, (|l| $\le$ 4, 2 <~ |b| $\le$ 4). The number density of RGBB stars is (12.7 +/- 2.0)% that of RC stars. The brightness dispersion of the RGBB is significantly lower than that of the RC, a result that is difficult to explain as the luminosity of the RGBB is known to significantly vary with metallicity. Sightlines toward the Galactic bulge that have two RCs have two RGBBs with similar properties to one another, an expected outcome if the Milky Way's bulge is X-shaped. We also find preliminary evidence of the Galactic bulge asymptotic giant branch bump, at a brightness of ~1.1 mag brighter than the RC in I and with a number density ~1.5% that of the RC. Accounting for the RGBB has a small effect on the best-fit parameters of the RC, shifting its best-fit peak brightness and reducing its brightness dispersion by ~0.015 mag each.

We report on photometric I-band observations of 147 bright (8<I<13) periodic variables toward the Galactic bulge including 76 new discoveries. We used one of the HATnet telescopes to obtain 151 exposures spanning 88 nights in 2005 of an 8.4 x 8.4 deg^2 field of view (FOV) approximately centered on (l,b) = (1.73, -4.68). We observed the galactic bulge in 2005 as part of a microlensing feasibility study (Nataf et al. 2009), and here we discuss the periodic variables we found in our data. Among our discoveries we count 52 new eclipsing binaries and 24 other periodic variables.

The red clump is found to be split into two components along several sightlines toward the Galactic Bulge. This split is detected with high significance toward the areas (-3.5<l<1,b<-5) and (l,b)=(0,+5.2), i.e., along the Bulge minor axis and at least 5 degrees off the plane. The fainter (hereafter ``main'') component is the one that more closely follows the distance-longitude relation of the Bulge red clump. The main component is ~0.5 magnitudes fainter than the secondary component and with an overall approximately equal population. For sightlines further from the plane, the difference in brightness increases, and more stars are found in the secondary component than in the main component. The two components have very nearly equal (V-I) color.

Following the suggestion of Gould and Depoy (1998) we investigate the feasibility of studying the brightest microlensing events towards the Galactic bulge using a small aperture (~10 cm) telescope. We used one of the HAT telescopes to obtain 151 exposures spanning 88 nights in 2005 of an 8.4x8.4 square degree FOV centered on (l,b) = (2.85, -5.00). We reduced the data using image subtraction software. We find that such a search method can effectively contribute to monitoring bright microlensing events, as was advocated. Comparing this search method to the existing ones we find a dedicated bulge photometric survey of this nature would fulfill a significant niche at excellent performance and rather low cost. We obtain matches to 7 microlensing events listed in the 2005 OGLE archives. We find several other light curves whose fits closely resemble microlensing events. Unsurprisingly, many periodic variables and miscellaneous variables are also detected in our data, and we estimate approximately 50% of these are new discoveries. We conclude by briefly proposing Small Aperture Microlensing Survey, which would monitor the Galactic bulge around the clock to provide dense coverage of the highest magnification microlensing events.