Showing 1–50 of 70 results for author: Nataf, D M

Accurate, Precise, and Physically Self-consistent Ages and Metallicities for 400,000 Solar Neighborhood Subgiant Branch Stars

Authors: David M. Nataf, Kevin C. Schlaufman, Henrique Reggiani, Isabel Hahn

Abstract: Age is the most difficult fundamental stellar parameter to infer for isolated stars. While isochrone-based ages are in general imprecise for both main sequence dwarfs and red giants, precise isochrone-based ages can be obtained for stars on the subgiant branch transitioning from core to shell hydrogen burning. We synthesize Gaia DR3-based distance inferences, multiwavelength photometry from the ul… ▽ More Age is the most difficult fundamental stellar parameter to infer for isolated stars. While isochrone-based ages are in general imprecise for both main sequence dwarfs and red giants, precise isochrone-based ages can be obtained for stars on the subgiant branch transitioning from core to shell hydrogen burning. We synthesize Gaia DR3-based distance inferences, multiwavelength photometry from the ultraviolet to the mid infrared, and three-dimensional extinction maps to construct a sample of 289,759 solar-metallicity stars amenable to accurate, precise, and physically self-consistent age inferences. Using subgiants in the solar-metallicity open clusters NGC 2682 (i.e., M 67) and NGC 188, we show that our approach yields accurate and physically self-consistent ages and metallicities with median statistical precisions of 8\% and 0.06 dex. The inclusion of systematic uncertainties resulting from non-single or variable stars results in age and metallicity precisions of 9\% and 0.12 dex. We supplement this solar-metallicity sample with an additional 112,062 metal-poor subgiants, including over 3,000 stars with $[\text{Fe/H}]\lesssim-1.50$, 7\% age precisions, and apparent Gaia $G$-band magnitudes $G<14$. We further demonstrate that our inferred metallicities agree with those produced by multiplexed spectroscopic surveys. As an example of the scientific potential of this catalog, we show that the solar neighborhood star-formation history has three components at $([\text{Fe/H}],τ/\text{Gyr}) \approx (+0.0,4)$, $(+0.2,7)$, and a roughly linear sequence in age--metallicity space beginning at $([\text{Fe/H}],τ/\text{Gyr})\approx(+0.2,7)$ and extending to $(-0.5,13)$. Our analyses indicate that the solar neighborhood includes stars on disk-like orbits even at the oldest ages and lowest metallicities accessible by our samples. △ Less

Submitted 25 July, 2024; originally announced July 2024.

Comments: 31 pages, 13 figures, accepted for publication in the AAS Journals. Our input data tables, output tables of derived parameters, and an animation of the derived relationship between the dynamics and the star-formation history of solar neighbourhood subgiants, are available at: https://github.com/DavidMoiseNataf/Subgiants/

The Galactic Bulge exploration III.: Calcium Triplet Metallicities for RR Lyrae Stars

Authors: Andrea Kunder, Zdenek Prudil, Claire Skaggs, Henrique Reggiani, David M. Nataf, Joanne Hughes, Kevin R. Covey, Kathryn Devine

Abstract: RR Lyrae stars are excellent tracers of stellar populations for old, metal-poor components in the Milky Way Galaxy and the Local Group. Their luminosities have a metallicity-dependence, but determining spectroscopic [Fe/H] metallicities for RR Lyrae stars, especially at distances outside the solar neighbourhood, is challenging. Using 40 RRLs with metallicities derived from both Fe(II) and Fe(I) ab… ▽ More RR Lyrae stars are excellent tracers of stellar populations for old, metal-poor components in the Milky Way Galaxy and the Local Group. Their luminosities have a metallicity-dependence, but determining spectroscopic [Fe/H] metallicities for RR Lyrae stars, especially at distances outside the solar neighbourhood, is challenging. Using 40 RRLs with metallicities derived from both Fe(II) and Fe(I) abundances, we verify the calibration between the [Fe/H] of RR Lyrae stars from the Calcium triplet. Our calibration is applied to all RR Lyrae stars with Gaia RVS spectra in Gaia DR3 as well as to 80 stars in the inner Galaxy from the BRAVA-RR survey. The co-added Gaia RVS RR Lyrae spectra provide RR Lyrae metallicities with an uncertainty of 0.25~dex, which is a factor of two improvement over the Gaia photometric RR Lyrae metallicities. Within our Galactic bulge RR Lyrae star sample, we find a dominant fraction with low energies without a prominent rotating component. Due to the large fraction of such stars, we interpret these stars as belonging to the $in-situ$ metal-poor Galactic bulge component, although we can not rule out that a fraction of these belong to an ancient accretion event such as Kraken/Heracles. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: accepted for publication in The Astronomical Journal

Is Terzan 5 the remnant of a building block of the Galactic bulge? Evidence from APOGEE

Authors: Dominic J. Taylor, Andrew C. Mason, Ricardo P. Schiavon, Danny Horta, David M. Nataf, Doug Geisler, Shobhit Kisku, Siân G. Phillips, Roger E. Cohen, José G. Fernández-Trincado, Timothy C. Beers, Dmitry Bizyaev, Domingo Aníbal García-Hernández, Richard R. Lane, Penélope Longa-Peña, Dante Minniti, Cesar Muñoz, Kaike Pan, Sandro Villanova

Abstract: It has been proposed that the globular cluster-like system Terzan 5 is the surviving remnant of a primordial building block of the Milky Way bulge, mainly due to the age/metallicity spread and the distribution of its stars in the $α$-Fe plane. We employ Sloan Digital Sky Survey (SDSS-IV) data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE 2) to test this hypothesis. Adopti… ▽ More It has been proposed that the globular cluster-like system Terzan 5 is the surviving remnant of a primordial building block of the Milky Way bulge, mainly due to the age/metallicity spread and the distribution of its stars in the $α$-Fe plane. We employ Sloan Digital Sky Survey (SDSS-IV) data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE 2) to test this hypothesis. Adopting a random sampling technique, we contrast the abundances of 10 elements in Terzan 5 stars with those of their bulge field counterparts with comparable atmospheric parameters, finding that they differ at statistically significant levels. Abundances between the two groups differ by more than 1$σ$ in Ca, Mn, C, O, and Al, and more than 2$σ$ in Si and Mg. Terzan 5 stars have lower [$α$/Fe] and higher [Mn/Fe] than their bulge counterparts. Given those differences, we conclude that Terzan 5 is not the remnant of a $major$ building block of the bulge. We also estimate the stellar mass of the Terzan 5 progenitor based on predictions by the Evolution and Assembly of GaLaxies and their Environments (EAGLE) suite of cosmological numerical simulations, concluding that it may have been as low as $\sim3\times10^8$ M$_\odot$ so that it was likely unable to significantly influence the mean chemistry of the bulge/inner disk, which is significantly more massive ($\sim10^{10}$ M$_\odot$). We briefly discuss existing scenarios for the nature of Terzan 5 and propose an observational test that may help elucidate its origin. △ Less

Submitted 4 April, 2022; originally announced April 2022.

Comments: 15 pages, 10 figures, accepted for publication by MNRAS

Evidence for a high-energy tail in the gamma-ray spectra of globular clusters

Authors: Deheng Song, Oscar Macias, Shunsaku Horiuchi, Roland M. Crocker, David M. Nataf

Abstract: Millisecond pulsars are very likely the main source of gamma-ray emission from globular clusters. However, the relative contributions of two separate emission processes--curvature radiation from millisecond pulsar magnetospheres vs. inverse Compton emission from relativistic pairs launched into the globular cluster environment by millisecond pulsars--have long been unclear. To address this, we sea… ▽ More Millisecond pulsars are very likely the main source of gamma-ray emission from globular clusters. However, the relative contributions of two separate emission processes--curvature radiation from millisecond pulsar magnetospheres vs. inverse Compton emission from relativistic pairs launched into the globular cluster environment by millisecond pulsars--have long been unclear. To address this, we search for evidence of inverse Compton emission in 8-year $\textit{Fermi}$-LAT data from the directions of 157 Milky Way globular clusters. We find a mildly statistically significant (3.8$σ$) correlation between the measured globular cluster gamma-ray luminosities and their photon field energy densities. However, this may also be explained by a hidden correlation between the photon field densities and the stellar encounter rates of globular clusters. Analysed $\textit{in toto}$, we demonstrate that the gamma-ray emission of globular clusters can be resolved spectrally into two components: i) an exponentially cut-off power law and ii) a pure power law. The latter component--which we uncover at a significance of 8.2$σ$--has a power index of 2.79 $\pm$ 0.25. It is most naturally interpreted as inverse Compton emission by cosmic-ray electrons and positrons injected by millisecond pulsars. We find the luminosity of this power-law component is comparable to, or slightly smaller than, the luminosity of the curved component, suggesting the fraction of millisecond pulsar spin-down luminosity into relativistic leptons is similar to the fraction of the spin-down luminosity into prompt magnetospheric radiation. △ Less

Submitted 26 September, 2021; v1 submitted 29 January, 2021; originally announced February 2021.

Comments: 16+9 pages, 13+4 figures, 4+1 tables. match the published version in MNRAS

Journal ref: MNRAS 507, 5161-5176 (2021)

Buckling bars in nearly face-on galaxies observed with MaNGA

Authors: K. M. Xiang, D. M. Nataf, E. Athanassoula, N. L. Zakamska, K. Rowlands, K. Masters, A. Fraser-McKelvie, N. Drory, K. Kraljic

Abstract: Over half of disk galaxies are barred, yet the mechanisms for bar formation and the life-time of bar buckling remain poorly understood. In simulations, a thin bar undergoes a rapid (<1 Gyr) event called "buckling," during which the inner part of the bar is asymmetrically bent out of the galaxy plane and eventually thickens, developing a peanut/X-shaped profile when viewed side-on. Through analyzin… ▽ More Over half of disk galaxies are barred, yet the mechanisms for bar formation and the life-time of bar buckling remain poorly understood. In simulations, a thin bar undergoes a rapid (<1 Gyr) event called "buckling," during which the inner part of the bar is asymmetrically bent out of the galaxy plane and eventually thickens, developing a peanut/X-shaped profile when viewed side-on. Through analyzing stellar kinematics of N-body model snapshots of a galaxy before, during, and after the buckling phase, we confirm a distinct quadrupolar pattern of out-of-plane stellar velocities in nearly face-on galaxies. This kinematic signature of buckling allows us to identify five candidates of currently buckling bars among 434 barred galaxies in the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) Survey, an integral field unit (IFU) spectroscopic survey that measures the composition and kinematic structure of nearby galaxies. The frequency of buckling events detected is consistent with the 0.5-1 Gyr timescale predicted by simulations. The five candidates we present more than double the total number of candidate buckling bars, and are the only ones found using the kinematic signature. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: 22 pages, 12 figures. Accepted for publication in ApJ

The GALAH+ Survey: Third Data Release

Authors: Sven Buder, Sanjib Sharma, Janez Kos, Anish M. Amarsi, Thomas Nordlander, Karin Lind, Sarah L. Martell, Martin Asplund, Joss Bland-Hawthorn, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Michael R. Hayden, Geraint F. Lewis, Jane Lin, Katharine J. Schlesinger, Jeffrey D. Simpson, Dennis Stello, Daniel B. Zucker, Tomaz Zwitter, Kevin L. Beeson, Tobias Buck, Luca Casagrande, Jake T. Clark , et al. (22 additional authors not shown)

Abstract: The ensemble of chemical element abundance measurements for stars, along with precision distances and orbit properties, provides high-dimensional data to study the evolution of the Milky Way. With this third data release of the Galactic Archaeology with HERMES (GALAH) survey, we publish 678 423 spectra for 588 571 mostly nearby stars (81.2% of stars are within <2 kpc), observed with the HERMES spe… ▽ More The ensemble of chemical element abundance measurements for stars, along with precision distances and orbit properties, provides high-dimensional data to study the evolution of the Milky Way. With this third data release of the Galactic Archaeology with HERMES (GALAH) survey, we publish 678 423 spectra for 588 571 mostly nearby stars (81.2% of stars are within <2 kpc), observed with the HERMES spectrograph at the Anglo-Australian Telescope. This release (hereafter GALAH+ DR3) includes all observations from GALAH Phase 1 (bright, main, and faint survey, 70%), K2-HERMES (17%), TESS-HERMES (5%), and a subset of ancillary observations (8%) including the bulge and >75 stellar clusters. We derive stellar parameters $T_\text{eff}$, $\log g$, [Fe/H], $v_\text{mic}$, $v_\text{broad}$ & $v_\text{rad}$ using our modified version of the spectrum synthesis code Spectroscopy Made Easy (SME) and 1D MARCS model atmospheres. We break spectroscopic degeneracies in our spectrum analysis with astrometry from $Gaia$ DR2 and photometry from 2MASS. We report abundance ratios [X/Fe] for 30 different elements (11 of which are based on non-LTE computations) covering five nucleosynthetic pathways. We describe validations for accuracy and precision, flagging of peculiar stars/measurements and recommendations for using our results. Our catalogue comprises 65% dwarfs, 34% giants, and 1% other/unclassified stars. Based on unflagged chemical composition and age, we find 62% young low-$α$, 9% young high-$α$, 27% old high-$α$, and 2% stars with $\mathrm{[Fe/H]} \leq -1$. Based on kinematics, 4% are halo stars. Several Value-Added-Catalogues, including stellar ages and dynamics, updated after $Gaia$ eDR3, accompany this release and allow chrono-chemodynamic analyses, as we showcase. △ Less

Submitted 28 April, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: 48 pages, 27 figures, 5 pages appendix, accepted for publication in MNRAS, comments welcome, GALAH+ DR3 can be found at https://docs.datacentral.org.au/galah/dr3/overview/, Catalogs can be downloaded from https://cloud.datacentral.org.au/teamdata/GALAH/public/GALAH_DR3/

The Similarity of Abundance Ratio Trends and Nucleosynthetic Patterns in the Milky Way Disk and Bulge

Authors: Emily Griffith, David H. Weinberg, Jennifer A. Johnson, Rachael Beaton, D. A. García-Hernández, Sten Hasselquist, Jon Holtzman, James W. Johnson, Henrik Jönsson, Richard R. Lane, David M. Nataf, Alexandre Roman-Lopes

Abstract: We compare abundance ratio trends in a sample of $\sim 11,000$ Milky Way bulge stars ($R_{\rm GC} < 3$ kpc) from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) to those of APOGEE stars in the Galactic disk ($5$ kpc $< R_{\rm GC} < 11$ kpc). We divide each sample into low-Ia (high-[Mg/Fe]) and high-Ia (low-[Mg/Fe]) populations, and in each population we examine the median trend… ▽ More We compare abundance ratio trends in a sample of $\sim 11,000$ Milky Way bulge stars ($R_{\rm GC} < 3$ kpc) from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) to those of APOGEE stars in the Galactic disk ($5$ kpc $< R_{\rm GC} < 11$ kpc). We divide each sample into low-Ia (high-[Mg/Fe]) and high-Ia (low-[Mg/Fe]) populations, and in each population we examine the median trends of [X/Mg] vs. [Mg/H] for elements X = Fe, O, Na, Al, Si, P, S, K, Ca, V, Cr, Mn, Co, Ni, Cu, and Ce. To remove small systematic trends of APOGEE abundances with stellar $\log(g)$, we resample the disk stars to match the $\log(g)$ distributions of the bulge data. After doing so, we find nearly identical median trends for low-Ia disk and bulge stars for all elements. High-Ia trends are similar for most elements, with noticeable (0.05-0.1 dex) differences for Mn, Na, and Co. The close agreement of abundance trends (with typical differences $\lesssim 0.03$ dex) implies that similar nucleosynthetic processes enriched bulge and disk stars despite the different star formation histories and physical conditions of these regions. For example, we infer that differences in the high mass slope of the stellar initial mass function (IMF) between disk and bulge must have been $\lesssim 0.30$. This agreement, and the generally small scatter about the median sequences, means that one can predict all of a bulge star's APOGEE abundances with good accuracy knowing only its measured [Mg/Fe] and [Mg/H] and the observed trends of disk stars. △ Less

Submitted 25 October, 2021; v1 submitted 10 September, 2020; originally announced September 2020.

Comments: 32 pages, 22 figures, 5 tables

Exploring the stellar age distribution of the Milky Way Bulge using APOGEE

Authors: Sten Hasselquist, Gail Zasowski, Diane K. Feuillet, Mathias Schultheis, David M. Nataf, Borja Anguiano, Rachael L. Beaton, Timothy C. Beers, Roger E. Cohen, Katia Cunha, José G. Fernández-Trincado, D. A. García-Hernández, Doug Geisler, Jon A. Holtzman, Jennifer Johnson, Richard R. Lane, Steven R. Majewski, Christian Moni Bidin, Christian Nitschelm, Alexandre Roman-Lopes, Ricardo Schiavon, Verne V. Smith, Jennifer Sobeck

Abstract: We present stellar age distributions of the Milky Way (MW) bulge region using ages for $\sim$6,000 high-luminosity ($\log(g) < 2.0$), metal-rich ($\rm [Fe/H] \ge -0.5$) bulge stars observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). Ages are derived using {\it The Cannon} label-transfer method, trained on a sample of nearby luminous giants with precise parallaxes for w… ▽ More We present stellar age distributions of the Milky Way (MW) bulge region using ages for $\sim$6,000 high-luminosity ($\log(g) < 2.0$), metal-rich ($\rm [Fe/H] \ge -0.5$) bulge stars observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). Ages are derived using {\it The Cannon} label-transfer method, trained on a sample of nearby luminous giants with precise parallaxes for which we obtain ages using a Bayesian isochrone-matching technique. We find that the metal-rich bulge is predominantly composed of old stars ($>$8 Gyr). We find evidence that the planar region of the bulge ($|Z_{\rm GC}| \le 0.25$ kpc) enriched in metallicity, $Z$, at a faster rate ($dZ/dt \sim$ 0.0034 ${\rm Gyr^{-1}}$) than regions farther from the plane ($dZ/dt \sim$ 0.0013 ${\rm Gyr^{-1}}$ at $|Z_{\rm GC}| > 1.00$ kpc). We identify a non-negligible fraction of younger stars (age $\sim$ 2--5 Gyr) at metallicities of $\rm +0.2 < [Fe/H] < +0.4$. These stars are preferentially found in the plane ($|Z_{\rm GC}| \le 0.25$ kpc) and between $R_{\rm cy} \approx 2-3$ kpc, with kinematics that are more consistent with rotation than are the kinematics of older stars at the same metallicities. We do not measure a significant age difference between stars found in and outside of the bar. These findings show that the bulge experienced an initial starburst that was more intense close to the plane than far from the plane. Then, star formation continued at super-solar metallicities in a thin disk at 2 kpc $\lesssim R_{\rm cy} \lesssim$ 3 kpc until $\sim$2 Gyr ago. △ Less

Submitted 8 August, 2020; originally announced August 2020.

Comments: 21 pages, 15 figures, accepted for publication in ApJ

The contribution of N-rich stars to the Galactic stellar halo using APOGEE red giants

Authors: Danny Horta, J. Ted Mackereth, Ricardo P. Schiavon, Sten Hasselquist, Jo Bovy, Carlos Allende Prieto, Timothy C. Beers, Katia Cunha, D. A. García-Hernández, Shobhit S. Kisku, Richard R. Lane, Steven R. Majewski, Andrew C. Mason, David M. Nataf, Alexandre Roman-Lopes, Mathias Schultheis

Abstract: The contribution of dissolved globular clusters (GCs) to the stellar content of the Galactic halo is a key constraint on models for GC formation and destruction, and the mass assembly history of the Milky Way. Earlier results from APOGEE pointed to a large contribution of destroyed GCs to the stellar content of the inner halo, by as much as 25$\%$, which is an order of magnitude larger than previo… ▽ More The contribution of dissolved globular clusters (GCs) to the stellar content of the Galactic halo is a key constraint on models for GC formation and destruction, and the mass assembly history of the Milky Way. Earlier results from APOGEE pointed to a large contribution of destroyed GCs to the stellar content of the inner halo, by as much as 25$\%$, which is an order of magnitude larger than previous estimates for more distant regions of the halo. We set out to measure the ratio between N-rich and normal halo field stars, as a function of distance, by performing density modelling of halo field populations in APOGEE DR16. Our results show that at 1.5 kpc from the Galactic Centre, N-rich stars contribute a much higher 16.8$^{+10.0}_{-7.0}$$\%$ fraction to the total stellar halo mass budget than the 2.7$^{+1.0}_{-0.8}$$\%$ ratio contributed at 10 kpc. Under the assumption that N-rich stars are former GC members that now reside in the stellar halo field, and assuming the ratio between first-and second-population GC stars being 1:2, we estimate a total contribution from disrupted GC stars of the order of 27.5$^{+15.4}_{-11.5}$$\%$ at r = 1.5 kpc and 4.2$^{+1.5}_{-1.3}$$\%$ at r = 10 kpc. Furthermore, since our methodology requires fitting a density model to the stellar halo, we integrate such density within a spherical shell from 1.5-15 kpc in radius, and find a total stellar mass arising from dissolved and/or evaporated GCs of $M_{\mathrm{GC,total}}$ = 9.6$^{+4.0}_{-2.6}$ $\times$ 10$^{7}$ M$\odot$. △ Less

Submitted 12 October, 2022; v1 submitted 3 August, 2020; originally announced August 2020.

Comments: Paper accepted for Publication in MNRAS, correction of Fig 7

The Milky Way's bulge star formation history as constrained from its bimodal chemical abundance distribution

Authors: Jianhui Lian, Gail Zasowski, Sten Hasselquist, David M. Nataf, Daniel Thomas, Christian Moni Bidin, José G. Fernández-Trincado, D. A. Garcia-Hernandez, Richard R. Lane, Steven R. Majewski, Alexandre Roman-Lopes, Mathias Schultheis

Abstract: We conduct a quantitative analysis of the star formation history (SFH) of the Milky Way's bulge by exploiting the constraining power of its stellar [Fe/H] and [Mg/Fe] distribution functions. Using APOGEE data, we confirm the previously-established bimodal [Mg/Fe]--[Fe/H] distribution within 3 kpc of the inner Galaxy. Compared to that in the solar vicinity, the high-$α$ population in the bulge peak… ▽ More We conduct a quantitative analysis of the star formation history (SFH) of the Milky Way's bulge by exploiting the constraining power of its stellar [Fe/H] and [Mg/Fe] distribution functions. Using APOGEE data, we confirm the previously-established bimodal [Mg/Fe]--[Fe/H] distribution within 3 kpc of the inner Galaxy. Compared to that in the solar vicinity, the high-$α$ population in the bulge peaks at a lower [Fe/H]. To fit these observations, we use a simple but flexible star formation framework, which assumes two distinct stages of gas accretion and star formation, and systematically evaluate a wide multi-dimensional parameter space. We find that the data favor a three-phase SFH that consists of an initial starburst, followed by a rapid star formation quenching episode and a lengthy, quiescent secular evolution phase. The metal-poor, high-$α$ bulge stars ([Fe/H]<0.0 and [Mg/Fe]>0.15) are formed rapidly (<2 Gyr) during the early starburst. The density gap between the high- and low-$α$ sequences is due to the quenching process. The metal-rich, low-$α$ population ([Fe/H]>0.0 and [Mg/Fe]<0.15) then accumulates gradually through inefficient star formation during the secular phase. This is qualitatively consistent with the early SFH of the inner disk. Given this scenario, a notable fraction of young stars (age<5 Gyr) is expected to persist in the bulge. Combined with extragalactic observations, these results suggest that a rapid star formation quenching process is responsible for bimodal distributions in both the Milky Way's stellar populations and in the general galaxy population and thus plays a critical role in galaxy evolution. △ Less

Submitted 23 July, 2020; originally announced July 2020.

Comments: 16 pages, 12 figures. MNRAS in press

The Magellanic Edges Survey I. Description and First Results

Authors: L. R. Cullinane, A. D. Mackey, G. S. Da Costa, S. E. Koposov, V. Belokurov, D. Erkal, A. Koch, A. Kunder, D. M. Nataf

Abstract: We present an overview of, and first science results from, the Magellanic Edges Survey (MagES), an ongoing spectroscopic survey mapping the kinematics of red clump and red giant branch stars in the highly substructured periphery of the Magellanic Clouds. In conjunction with Gaia astrometry, MagES yields a sample of ~7000 stars with individual 3D velocities that probes larger galactocentric radii t… ▽ More We present an overview of, and first science results from, the Magellanic Edges Survey (MagES), an ongoing spectroscopic survey mapping the kinematics of red clump and red giant branch stars in the highly substructured periphery of the Magellanic Clouds. In conjunction with Gaia astrometry, MagES yields a sample of ~7000 stars with individual 3D velocities that probes larger galactocentric radii than most previous studies. We outline our target selection, observation strategy, data reduction and analysis procedures, and present results for two fields in the northern outskirts ($>10^{\circ}$ on-sky from the centre) of the Large Magellanic Cloud (LMC). One field, located in the vicinity of an arm-like overdensity, displays apparent signatures of perturbation away from an equilibrium disk model. This includes a large radial velocity dispersion in the LMC disk plane, and an asymmetric line-of-sight velocity distribution indicative of motions vertically out of the disk plane for some stars. The second field reveals 3D kinematics consistent with an equilibrium disk, and yields $V_{\text{circ}}=87.7\pm8.0$km s$^{-1}$ at a radial distance of ~10.5kpc from the LMC centre. This leads to an enclosed mass estimate for the LMC at this radius of $(1.8\pm0.3)\times10^{10}\text{M}_{\odot}$. △ Less

Submitted 13 July, 2020; originally announced July 2020.

Comments: 23 pages, 10 figures. Accepted by MNRAS

On the Color-Metallicity Relation of the Red Clump and the Reddening Toward the Magellanic Clouds

Authors: David M. Nataf, Santi Cassisi, Luca Casagrande, Wenlong Yuan, Adam G. Riess

Abstract: 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 widely-used metric to predict LMC… ▽ More 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 widely-used metric to predict LMC extinction, the intrinsic color of the red clump, $(V-I)_{RC,0}$, for the inner $\sim$3 degrees of that galaxy. The first approach was to empirically calibrate the color zeropoints of the BaSTI isochrones over a wide metallicity range of $Δ\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. With these three methods, we estimated the intrinsic color of the red clump of the LMC to be $(V-I)_{RC,0,\rm{LMC}} = \{ \approx 0.93,0.91 \pm 0.02,0.89 \pm 0.02\}$ respectively, and similarly using the first and third method we estimated $ (V-I)_{RC,0,\rm{SMC}} = \{\approx 0.85,0.84 \pm 0.02 \}$ respectively for the Small Magellanic Cloud. We estimate the luminosities to be $M_{I,RC,\rm{LMC}}=-0.26$ and $M_{I,RC,\rm{SMC}}=-0.37$. 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$. △ Less

Submitted 10 February, 2021; v1 submitted 5 June, 2020; originally announced June 2020.

Comments: 20 pages, 5 figures, 4 tables. Accepted for Publication in The Astrophysical Journal as of February 9th, 2021

The Predicted Properties of Helium-Enriched Globular Cluster Progenitors at High Redshift

Authors: David M. Nataf, Shunsaku Horiuchi, Guglielmo Costa, Rosemary F. G. Wyse, Yuan-Sen Ting, Roland Crocker, Christoph Federrath, Yang Chen

Abstract: 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 solut… ▽ More 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 ($Δ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 $λ\approx 2.0\, μm$, and to be 0.30 mag redder in the \textit{JWST}-NIRCam colour $(F070W-F200W)$, and to actually be fainter for $λ\lesssim 0.50 \, μm$. Separately, we find that the time-integrated shift in ionizing radiation is a negligible $\sim 5\%$, though we show that the Lyman-$α$ escape fraction could end up higher for helium-enriched stars. △ Less

Submitted 17 August, 2020; v1 submitted 13 May, 2020; originally announced May 2020.

Comments: Accepted for publication in MNRAS on 13 May, 2020 Replaced to update some bibliographic information

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 496, Issue 3, pp.3222-3234, 2020

Fundamental relations for the velocity dispersion of stars in the Milky Way

Authors: Sanjib Sharma, Michael R. Hayden, Joss Bland-Hawthorn, Dennis Stello, Sven Buder, Joel C. Zinn, Thomas Kallinger, Martin Asplund, Gayandhi M. De Silva, Valentina Dorazi, Ken Freeman, Janez Kos, Geraint F. Lewis, Jane Lin, Karin Lind, Sarah Martell, Jeffrey D. Simpson, Rob A. Wittenmyer, Daniel B. Zucker, Tomaz Zwitter, Boquan Chen, Klemen Cotar, James Esdaile, Marc Hon, Jonathan Horner , et al. (12 additional authors not shown)

Abstract: 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… ▽ More 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 $σ_z$ and $σ_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 $σ_z$ increases (implying flaring in the outer disc) while $σ_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-[$α$/Fe] stars follow the same relations as the low-[$α$/Fe] stars. △ Less

Submitted 14 April, 2020; originally announced April 2020.

Comments: 19 pages, 15 figures, submitted to MNRAS

The Chemical Compositions of Accreted and {\it in situ} Galactic Globular Clusters According to SDSS/APOGEE

Authors: Danny Horta, Ricardo P. Schiavon, J. Ted Mackereth, Timothy C. Beers, José G. Fernández-Trincado, Peter M. Frinchaboy, D. A. García-Hernández, Doug Geisler, Sten Hasselquist, Henrik Jönsson, Richard R. Lane, Steven R. Majewski, Szabolcs Mészáros, Christian Moni Bidin, David M. Nataf, Alexandre Roman-Lopes, Christian Nitschelm, Jaime Vargas-González, Gail Zasowski

Abstract: 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 s… ▽ More 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 $α$ (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. △ Less

Submitted 14 February, 2020; v1 submitted 9 January, 2020; originally announced January 2020.

Comments: 18 pages, 10 figures, 3 tables. Accepted for publication in MNRAS

The K2-HERMES Survey: Age and Metallicity of the Thick Disc

Authors: Sanjib Sharma, Dennis Stello, Joss Bland-Hawthorn, Michael R. Hayden, Joel C. Zinn, Thomas Kallinger, Marc Hon, Martin Asplund, Sven Buder, Gayandhi M. De Silva, Valentina D'Orazi, Ken Freeman, Janez Kos, Geraint F. Lewis, Jane Lin, Karin Lind, Sarah Martell, Jeffrey D. Simpson, Rob A. Wittenmyer, Daniel B. Zucker, Tomaz Zwitter, Timothy R. Bedding, Boquan Chen, Klemen Cotar, James Esdaile , et al. (12 additional authors not shown)

Abstract: 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 m… ▽ More 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 [$α$/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 [$α$/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 $α$-enhanced thick disc. △ Less

Submitted 29 April, 2019; originally announced April 2019.

Comments: 21 pages, submitted to MNRAS

The Relationship Between Globular Cluster Mass, Metallicity, and Light Element Abundance Variations

Authors: David M. Nataf, Rosemary Wyse, Ricardo P. Schiavon, Yuan-Sen Ting, Dante Minniti, Roger E. Cohen, José G. Fernández-Trincado, Douglas Geisler, Christian Nitschelm, Peter M. Frinchaboy

Abstract: 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 de… ▽ More 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. △ Less

Submitted 16 April, 2019; originally announced April 2019.

Comments: 39 pages, 24 figures, 8 tables, accepted for publication in The Astronomical Journal

HERBS II: Detailed chemical compositions of Galactic bulge stars

Authors: L. Duong, M. Asplund, D. M. Nataf, K. C. Freeman, M. Ness

Abstract: 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 obse… ▽ 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. △ Less

Submitted 3 April, 2019; originally announced April 2019.

Comments: 14 pages, 9 figures. Submitted to MNRAS July 2018; revised version. Online data access

HERBS I: Metallicity and alpha enhancement along the Galactic bulge minor axis

Authors: L. Duong, M. Asplund, D. M. Nataf, K. C. Freeman, M. Ness, L. M. Howes

Abstract: 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… ▽ More 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 $α$-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-$α$ 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. △ Less

Submitted 22 March, 2019; originally announced March 2019.

Comments: 18 pages, 14 figures. Submitted to MNRAS July 2018, revised version. Online data access

"Auxiliary" Science with the WFIRST Microlensing Survey

Authors: B. Scott Gaudi, Rachel Akeson, Jay Anderson, Etienne Bachelet, David P. Bennett, Aparna Bhattacharya, Valerio Bozza, Sebastiano Calchi Novati, Calen B. Henderson, Samson A. Johnson, Jeffrey Kruk, Jessica R. Lu, Shude Mao, Benjamin T. Montet, David M. Nataf, Matthew T. Penny, Radoslaw Poleski, Clément Ranc, Kailash Sahu, Yossi Shvartzvald, David N. Spergel, Daisuke Suzuki, Keivan G. Stassun, Rachel A. Street

Abstract: The Wide Field Infrared Survey Telescope (WFIRST) will monitor $\sim 2$ deg$^2$ toward the Galactic bulge in a wide ($\sim 1-2~μ$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, reachi… ▽ More The Wide Field Infrared Survey Telescope (WFIRST) will monitor $\sim 2$ deg$^2$ toward the Galactic bulge in a wide ($\sim 1-2~μ$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. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: 9 pages, 3 figures, submitted to the Astro2020 Science White Paper call. This white paper draws heavily from Chapter 2.5 of Spergel et al. 2015 (arXiv:1503.03757), which was written by the lead author of this white paper

The Scientific Context of WFIRST Microlensing in the 2020s

Authors: Jennifer Yee, Rachel Akeson, Jay Anderson, Etienne Bachelet, Charles Beichman, Andrea Bellini, David P. Bennett, Aparna Bhattacharya, Valerio Bozza, Geoffrey Bryden, Sebastiano Calchi Novati, B. Scott Gaudi, Andrew Gould, Calen B. Henderson, Savannah R. Jacklin, Samson A. Johnson, Naoki Koshimoto, Shude Mao, David M. Nataf, Matthew Penny, Radoslaw Poleski, Clément Ranc, Kailash Sahu, Yossi Shvartzvald, Keivan G. Stassun , et al. (1 additional authors not shown)

Abstract: [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… ▽ More [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. △ Less

Submitted 19 March, 2019; originally announced March 2019.

Comments: Science White Paper submitted to the Astro2020 Decadal Survey, 8 pages, 2 figures

High-Dimensional Dust Mapping

Authors: G. Zasowski, D. P. Finkbeiner, G. M. Green, J. A. Kollmeier, D. M. Nataf, J. E. G. Peek, E. F. Schlafly, V. Silva Aguirre, J. S. Speagle, K. Tchernyshyov, J. D. Trujillo, C. Zucker

Abstract: 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)… ▽ More 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. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 call for science white papers. 5 pages plus references, 1 figure

Discovery of a nitrogen-enhanced mildly metal-poor binary system: Possible evidence for pollution from an extinct AGB Star

Authors: José G. Fernández-Trincado, Ronald Mennickent, Mauricio Cabezas, Olga Zamora, Sarah L. Martell, Timothy C. Beers, Vinicius M. Placco, David M. Nataf, Szabolcs Mészáros, Dante Minniti, Dominik R. G. Schleicher, Baitian Tang, Angeles Pérez-Villegas, Annie C. Robin, Céline Reylé, Mario Ortigoza-Urdaneta

Abstract: 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 cl… ▽ More 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. △ Less

Submitted 3 September, 2019; v1 submitted 27 February, 2019; originally announced February 2019.

Comments: Accepted for publication in A&A, 13 pages, 8 figures, 7 tables

Journal ref: A&A 631, A97 (2019)

Strong Evidence that the Galactic Bulge is Shining in Gamma Rays

Authors: Oscar Macias, Shunsaku Horiuchi, Manoj Kaplinghat, Chris Gordon, Roland M. Crocker, David M. Nataf

Abstract: 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… ▽ More 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σ$ 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. △ Less

Submitted 20 September, 2019; v1 submitted 12 January, 2019; originally announced January 2019.

Comments: 21 pages, 13 figures, V2: Minor changes to match submitted version, V3: matches JCAP published version

Journal ref: JCAP09(2019)042

The Mira-based distance to the Galactic centre

Authors: Wenzer Qin, David M. Nataf, Nadia Zakamska, Peter R. Wood, Luca Casagrande

Abstract: 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 compone… ▽ More 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 $ΔM_{Ks} \approx -0.109(Δ\rm{[Fe/H]}) + 0.033( Δt/\rm{Gyr}) + 0.021 (Δ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. △ Less

Submitted 3 August, 2018; originally announced August 2018.

Comments: Accepted for publication in The Astrophysical Journal. 16 pages, 8 figures, 6 tables

Journal ref: 2018, ApJ, 865, 47

The GALAH survey: a catalogue of carbon-enhanced stars and CEMP candidates

Authors: Klemen Čotar, Tomaž Zwitter, Janez Kos, Ulisse Munari, Sarah L. Martell, Martin Asplund, Joss Bland-Hawthorn, Sven Buder, Gayandhi M. De Silva, Kenneth C. Freeman, Sanjib Sharma, Borja Anguiano, Daniela Carollo, Jonathan Horner, Geraint F. Lewis, David M. Nataf, Thomas Nordlander, Denis Stello, Yuan-Sen Ting, Chris Tinney, Gregor Traven, Rob A. Wittenmyer

Abstract: 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 ca… ▽ More 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. △ Less

Submitted 20 November, 2018; v1 submitted 20 July, 2018; originally announced July 2018.

Comments: 13+5 pages, 13 figures, 1 catalog, accepted to MNRAS

The GALAH Survey: Velocity fluctuations in the Milky Way using red clump giants

Authors: Shourya Khanna, Sanjib Sharma, Joss Bland-Hawthorn, Michael Hayden, David M. Nataf, Yuan-Sen Ting, Janez Kos, Sarah Martell, Tomaz Zwitter, Gayandhi De Silva, Martin Asplund, Sven Buder, Ly Duong, Jane Lin, Jeffrey D. Simpson, Borja Anguiano, Jonathan Horner, Prajwal R. Kafle, Geraint F. Lewis, Thomas Nordlander, Rosemary F. G. Wyse, Robert A. Wittenmyer, Daniel B. Zucker

Abstract: 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 eve… ▽ More 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. △ Less

Submitted 3 November, 2018; v1 submitted 19 April, 2018; originally announced April 2018.

Comments: 20 pages, 16 figures, accepted for publication in MNRAS

The GALAH Survey: Verifying abundance trends in the open cluster M67 using non-LTE spectroscopy

Authors: Xudong Gao, Karin Lind, Anish M. Amarsi, Sven Buder, Aaron Dotter, Thomas Nordlander, Martin Asplund, Joss Bland-Hawthorn, Gayandhi M. De Silva, Valentina D'Orazi, Ly Duong, Ken C. Freeman, Janez Kos, Geraint F. Lewis, Jane Lin, Sarah L. Martell, Katharine. J. Schlesinger, Sanjib Sharma, Jeffrey D. Simpson, Daniel B. Zucker, Tomaz Zwitter, Gary Da Costa, Borja Anguiano, Jonathan Horner, Elaina A. Hyde , et al. (6 additional authors not shown)

Abstract: 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 cl… ▽ More 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. △ Less

Submitted 23 August, 2018; v1 submitted 17 April, 2018; originally announced April 2018.

Comments: 14 pages, 10 figures; resubmitted to MNRAS

The GALAH Survey: Accurate Radial Velocities and Library of Observed Stellar Template Spectra

Authors: Tomaž Zwitter, Janez Kos, Andrea Chiavassa, Sven Buder, Gregor Traven, Klemen Čotar, Jane Lin, Martin Asplund, Joss Bland-Hawthorn, Andrew R. Casey, Gayandhi De Silva, Ly Duong, Kenneth C. Freeman, Karin Lind, Sarah Martell, Valentina D'Orazi, Katharine J. Schlesinger, Jeffrey D. Simpson, Sanjib Sharma, Daniel B. Zucker, Borja Anguiano, Luca Casagrande, Remo Collet, Jonathan Horner, Michael J. Ireland , et al. (12 additional authors not shown)

Abstract: 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,… ▽ More 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. △ Less

Submitted 6 September, 2018; v1 submitted 17 April, 2018; originally announced April 2018.

Comments: Accepted to MNRAS, compared to the initial submission it now reports radial velocities for more stars (336,215 instead of 199,628) and all radial velocities have been recalculated

Journal ref: 2018MNRAS.481..645Z

The GALAH Survey: Second Data Release

Authors: S. Buder, M. Asplund, L. Duong, J. Kos, K. Lind, M. K. Ness, S. Sharma, J. Bland-Hawthorn, A. R. Casey, G. M. De Silva, V. D'Orazi, K. C. Freeman, G. F. Lewis, J. Lin, S. L. Martell, K. J. Schlesinger, J. D. Simpson, D. B. Zucker, T. Zwitter, A. M. Amarsi, B. Anguiano, D. Carollo, K. Cotar, P. L. Cottrell, G. Da Costa , et al. (17 additional authors not shown)

Abstract: 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 abundance… ▽ More 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. △ Less

Submitted 17 April, 2018; originally announced April 2018.

Comments: 34+7 pages, 31 Figures, 5 tables, 1 catalog, submitted to MNRAS

The GALAH and TESS-HERMES surveys: high-resolution spectroscopy of luminous supergiants in the Magellanic Clouds and Bridge

Authors: Jeffrey D. Simpson, Dennis Stello, Sanjib Sharma, Yuan-Sen Ting, David M. Nataf, Gary Da Costa, Robert A. Wittenmyer, Jonathan Horner, Sarah L. Martell, Geraint F. Lewis, Gayandhi M. De Silva, Peter L. Cottrell, Martin Asplund, Joss Bland-Hawthorn, Sven Buder, Valentina D'Orazi, Ly Duong, Ken C. Freeman, Janez Kos, Jane Lin, Karin Lind, Katharine. J. Schlesinger, Daniel B. Zucker, Tomaž Zwitter, Prajwal R. Kafle , et al. (2 additional authors not shown)

Abstract: 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 ha… ▽ More 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. △ Less

Submitted 19 April, 2018; v1 submitted 16 April, 2018; originally announced April 2018.

Comments: 16 pages, 9 figures, submitted to MNRAS (updated to included missed co-author)

The GALAH survey: Co-orbiting stars and chemical tagging

Authors: Jeffrey D. Simpson, Sarah L. Martell, Gary Da Costa, Andrew R. Casey, Ken C. Freeman, Jonathan Horner, Yuan-Sen Ting, David M. Nataf, Geraint F. Lewis, Melissa K. Ness, Daniel B. Zucker, Peter L. Cottrell, Klemen Čotar, Martin Asplund, Joss Bland-Hawthorn, Sven Buder, Valentina D'Orazi, Gayandhi M. De Silva, Ly Duong, Janez Kos, Jane Lin, Karin Lind, Katharine J. Schlesinger, Sanjib Sharma, Tomaž Zwitter , et al. (2 additional authors not shown)

Abstract: 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-… ▽ More 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. △ Less

Submitted 19 June, 2018; v1 submitted 16 April, 2018; originally announced April 2018.

Comments: 14 pages, 9 figures, submitted to MNRAS. Updated for Gaia DR2 values

The GALAH survey: An abundance, age, and kinematic inventory of the solar neighbourhood made with TGAS

Authors: S. Buder, K. Lind, M. K. Ness, M. Asplund, L. Duong, J. Lin, J. Kos, L. Casagrande, A. R. Casey, J. Bland-Hawthorn, G. M. De Silva, V. D'Orazi, K. C. Freeman, S. L. Martell, K. J. Schlesinger, S. Sharma, J. D. Simpson, D. B. Zucker, T. Zwitter, K. Cotar, A. Dotter, M. R. Hayden, E. A. Hyde, P. R. Kafle, G. F. Lewis , et al. (9 additional authors not shown)

Abstract: 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… ▽ More 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-$α$ sequence are typically older than stars in the low-$α$ sequence, the latter spanning $-0.7<$[Fe/H]$<+0.5$. While these two sequences become indistinguishable in [$α$/Fe] vs. [Fe/H] at the metal-rich regime, we find that age can be used to separate stars from the extended high-$α$ & the low-$α$ 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-$α$-sequence up to super-solar [Fe/H] rather than a gap, which has been interpreted as a separate "high-$α$ metal-rich" population. Stars in our sample that are younger than 10 Gyr, are mainly found on the low $α$-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. △ Less

Submitted 1 December, 2018; v1 submitted 16 April, 2018; originally announced April 2018.

Comments: Revised paper resubmitted to A&A, 20 pages (+10 pages appendix), 25 figures, 2 tables, 1 catalog

Journal ref: A&A 624, A19 (2019)

Holistic spectroscopy: Complete reconstruction of a wide-field, multi-object spectroscopic image using a photonic comb

Authors: Janez Kos, Joss Bland-Hawthorn, Christopher H. Betters, Sergio Leon-Saval, Martin Asplund, Sven Buder, Andrew R. Casey, Valentina D'Orazi, Gayandhi de Silva, Ken Freeman, Geraint Lewis, Jane Lin, Sarah L. Martell, Katharine Schlesinger, Sanjib Sharma, Jeffrey D. Simpson, Daniel Zucker, Tomaž Zwitter, Michael Hayden, David M. Nataf, Yuan-Sen Ting

Abstract: 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… ▽ More 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. △ Less

Submitted 16 April, 2018; originally announced April 2018.

Comments: Submitted to MNRAS

White Paper: Exoplanetary Microlensing from the Ground in the 2020s

Authors: Jennifer C. Yee, Jay Anderson, Rachel Akeson, Etienne Bachelet, Charles Beichman, Andrea Bellini, David Bennett, Aparna Bhattacharya, Valerio Bozza, Sebastiano Calchi Novati, Will Clarkson, David R. Ciardi, Andrew Gould, Calen B. Henderson, Savannah R. Jacklin, Somayeh Khakpash, Shude Mao, Bertrand Mennesson, David M. Nataf, Matthew Penny, Joshua Pepper, Radek Poleski, Clement Ranc, Kailash Sahu, Y. Shvartzvald , et al. (3 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 21 March, 2018; originally announced March 2018.

Comments: White paper submitted to the National Academy of Sciences Committee on Exoplanet Science Strategy; 6 pages

Radial velocities of RR Lyrae stars in and around NGC 6441

Authors: Andrea Kunder, Arthur Mills, Joseph Edgecomb, Mathew Thomas, Levi Schilter, Craig Boyle, Stephen Parker, Gordon Bellevue, R. Michael Rich, Andreas Koch, Christian I. Johnson, David M. Nataf

Abstract: 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 presu… ▽ More 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. △ Less

Submitted 26 February, 2018; originally announced February 2018.

Comments: accepted to AJ

UKIRT-2017-BLG-001Lb: A giant planet detected through the dust

Authors: Y. Shvartzvald, S. Calchi Novati, B. S. Gaudi, G. Bryden, D. M. Nataf, M. T. Penny, C. Beichman, C. B. Henderson, S. Jacklin, E. F. Schlafly, M. J. Huston

Abstract: 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-inf… ▽ More 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. △ Less

Submitted 27 March, 2018; v1 submitted 19 February, 2018; originally announced February 2018.

Comments: 27 pages, 8 figures, 2 tables. Accepted for publication in ApJ Letters

The GALAH Survey: Stellar streams and how stellar velocity distributions vary with Galactic longitude, hemisphere and metallicity

Authors: Alice C. Quillen, Gayandhi De Silva, Sanjib Sharma, Michael Hayden, Ken Freeman, Joss Bland-Hawthorn, Maruša Žerjal, Martin Asplund, Sven Buder, Valentina D'Orazi, Ly Duong, Janez Kos, Jane Lin, Karin Lind, Sarah Martell, Katharine Schlesinger, Jeffrey D. Simpson, Daniel B. Zucker, Tomaz Zwitter, Borja Anguiano, Daniela Carollo, Luca Casagrande, Klemen Cotar, Peter L. Cottrell, Michael Ireland , et al. (8 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 4 April, 2018; v1 submitted 8 February, 2018; originally announced February 2018.

Comments: accepted for publication in MNRAS

The GALAH survey: properties of the Galactic disk(s) in the solar neighbourhood

Authors: L. Duong, K. C. Freeman, M. Asplund, L. Casagrande, S. Buder, K. Lind, M. Ness, J. Bland-Hawthorn, G. M. De Silva, V. D'Orazi, J. Kos, G. F. Lewis, J. Lin, S. L. Martell, K. Schlesinger, S. Sharma, J. D. Simpson, D. B. Zucker, T. Zwitter, B. Anguiano, G. S. Da Costa, E. Hyde, J. Horner, P. R. Kafle, D. M. Nataf , et al. (4 additional authors not shown)

Abstract: 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$). Thi… ▽ More 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-$α$ 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 $α$-abundance profile is almost flat, with a gradient of d[$α$/M]/d$z$ = $0.008 \pm 0.002$ dex kpc$^{-1}$. The steep vertical metallicity gradient of the low-$α$ population is in agreement with models where radial migration has a major role in the evolution of the thin disk. The thick disk (high-$α$ population) has a weaker vertical metallicity gradient d[M/H]/d$z = -0.058 \pm 0.003$ dex kpc$^{-1}$. The $α$-abundance of the thick disk is nearly constant with height, d[$α$/M]/d$z$ = $0.007 \pm 0.002$ dex kpc$^{-1}$. The negative gradient in metallicity and the small gradient in [$α$/M] indicate that the high-$α$ population experienced a settling phase, but also formed prior to the onset of major SNIa enrichment. We explore the implications of the distinct $α$-enrichments and narrow [$α$/M] range of the sub-populations in the context of thick disk formation. △ Less

Submitted 26 February, 2018; v1 submitted 4 January, 2018; originally announced January 2018.

Comments: 18 pages, 15 figures, 3 tables. Accepted for publication in MNRAS

The GALAH survey: Chemical Tagging of Star Clusters and New Members in the Pleiades

Authors: Janez Kos, Joss Bland-Hawthorn, Ken Freeman, Sven Buder, Gregor Traven, Gayandhi M. De Silva, Sanjib Sharma, Martin Asplund, Ly Duong, Jane Lin, Karin Lind, Sarah Martell, Jeffrey D. Simpson, Dennis Stello, Daniel B. Zucker, Tomaž Zwitter, Borja Anguiano, Gary Da Costa, Jonathan Horner, Prajwal R. Kafle, Geraint Lewis, Ulisse Munari, David M. Nataf, Melissa Ness, Warren Reid , et al. (3 additional authors not shown)

Abstract: 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 no… ▽ More 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. △ Less

Submitted 3 September, 2017; originally announced September 2017.

Was the Milky Way Bulge Formed From The Buckling Disk Instability, Hierarchical Collapse, Accretion of Clumps, or All of the Above?

Authors: David M. Nataf

Abstract: 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… ▽ More 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. △ Less

Submitted 3 August, 2017; originally announced August 2017.

Comments: 13 pages, 5 figures, accepted for publication in PASA

Diffuse Galactic antimatter from faint thermonuclear supernovae in old stellar populations

Authors: Roland M. Crocker, Ashley J. Ruiter, Ivo R. Seitenzahl, Fiona H. Panther, Stuart Sim, Holger Baumgardt, Anais Moller, David M. Nataf, Lilia Ferrario, J. J. Eldridge, Martin White, Brad E. Tucker, Felix Aharonian

Abstract: 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… ▽ More 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. △ Less

Submitted 23 May, 2017; v1 submitted 12 July, 2016; originally announced July 2016.

Comments: 28 pages main text with 4 figures in preprint style; 26 pages of Supplementary Information

Journal ref: Nature Astronomy 1, Article number: 0135 (2017)

The EMBLA Survey -- Metal-poor stars in the Galactic bulge

Authors: Louise M. Howes, Martin Asplund, Stefan C. Keller, Andrew R. Casey, David Yong, Karin Lind, Anna Frebel, Austin Hays, Alan Alves-Brito, Michael S. Bessell, Luca Casagrande, Anna F. Marino, David M. Nataf, Christopher I. Owen, Gary S. Da Costa, Brian P. Schmidt, Patrick Tisserand

Abstract: 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 th… ▽ More 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. △ Less

Submitted 26 April, 2016; originally announced April 2016.

Comments: 20 pages, 16 figures. Accepted to Monthly Notices of the Royal Astronomical Society. Online tables can be found at http://www.astro.lu.se/~louise/data

The GALAH Survey: Relative throughputs of the 2dF fibre positioner and the HERMES spectrograph from stellar targets

Authors: Jeffrey D. Simpson, G. M. De Silva, J. Bland-Hawthorn, K. C. Freeman, S. L. Martell, Katharine J. Schlesinger, Sanjib Sharma, D. B. Zucker, T. Zwitter, J. Kos, Borja Anguiano, David M. Nataf, Warren Reid, Robert A. Wittenmyer

Abstract: 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 throughpu… ▽ More 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. △ Less

Submitted 29 March, 2016; originally announced March 2016.

Comments: 14 pages, 10 figures. Accepted for publication in MNRAS

The Interstellar Extinction Toward the Milky Way Bulge with Planetary Nebulae, Red Clump, and RR Lyrae stars

Authors: David M. Nataf

Abstract: 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 exti… ▽ More 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. △ Less

Submitted 22 March, 2016; originally announced March 2016.

Comments: 12 pages, 6 figures, 2 Tables, invited article to PASA special review issue on the Galactic Bulge, edited by Beatriz Barbuy. This version includes minor modifications requested by referee prior to being accepted for publication

Before the Bar: Kinematic Detection of A Spheroidal Metal-Poor Bulge Component

Authors: Andrea Kunder, R. M. Rich, J. Storm, D. M. Nataf, R. De Propris, A. R. Walker, G. Bono, C. I. Johnson, J. Shen, Z. Y. Li

Abstract: 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… ▽ More 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. △ Less

Submitted 21 March, 2016; originally announced March 2016.

Comments: accepted to ApJL

Extremely metal-poor stars from the cosmic dawn in the bulge of the Milky Way

Authors: L. M. Howes, A. R. Casey, M. Asplund, S. C. Keller, D. Yong, D. M. Nataf, R. Poleski, K. Lind, C. Kobayashi, C. I. Owen, M. Ness, M. S. Bessell, G. S. Da Costa, B. P. Schmidt, P. Tisserand, A. Udalski, M. K. Szymański, I. Soszyński, G. Pietrzyński, K. Ulaczyk, Ł. Wyrzykowski, P. Pietrukowicz, J. Skowron, S. Kozłowski, P. Mróz

Abstract: 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 i… ▽ More 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. △ Less

Submitted 12 November, 2015; originally announced November 2015.

Comments: Published in Nature on 11/11/15, see http://dx.doi.org/10.1038/nature15747

Interstellar Extinction Curve Variations Toward the Inner Milky Way: A Challenge to Observational Cosmology

Authors: David M. Nataf, Oscar A. Gonzalez, Luca Casagrande, Gail Zasowski, Christopher Wegg, Christian Wolf, Andrea Kunder, Javier Alonso-Garcia, Dante Minniti, Marina Rejkuba, Roberto K. Saito, Elena Valenti, Manuela Zoccali, Radoslaw Poleski, Grzegorz Pietrzynski, Jan Skowron, Igor Soszynski, Michal K. Szymanski, Andrzej Udalski, Krzystof Ulaczyk, Lukasz Wyrzykowski

Abstract: 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 accura… ▽ More 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. △ Less

Submitted 18 December, 2015; v1 submitted 5 October, 2015; originally announced October 2015.

Comments: 16 pages, 13 figures, 3 tables, accepted to MNRAS, small revisions due to editor's comments and suggestions from the international astronomical community

The Controversial Star-Formation History and Helium Enrichment of the Milky Way Bulge

Authors: David M. Nataf

Abstract: 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]… ▽ More 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. △ Less

Submitted 31 August, 2015; originally announced September 2015.

Comments: 9 pages, 2 figures, contributed article to PASA special review issue on the Galactic Bulge, edited by Beatriz Barbuy, accepted for publication in PASA

A high-velocity bulge RR Lyrae variable on a halo-like orbit

Authors: Andrea Kunder, R. Michael Rich, Keith Hawkins, Radek Poleski, Jesper Storm, Christian I. Johnson, Juntai Shen, Zhao-Yu Li, Maria Jose Cordero, David M. Nataf, Giuseppe Bono, Alistair R. Walker, Andreas Koch, Roberto De Propris, Andrzej Udalski, Michal K. Szymanski, Igor Soszynski, Grzegorz Pietrzynski, Krzysztof Ulaczyk, Lukasz Wyrzykowski, Pawel Pietrukowicz, Jan Skowron, Szymon Kozlowski, Przemyslaw Mroz

Abstract: 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… ▽ More 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. △ Less

Submitted 8 June, 2015; originally announced June 2015.

Comments: accepted to ApJ Letters