Showing 51–100 of 122 results for author: El-Badry, K

Testing White Dwarf Age Estimates using Wide Double White Dwarf Binaries from Gaia EDR3

Authors: Tyler M. Heintz, J. J. Hermes, Kareem El-Badry, Charlie Walsh, Jennifer L. van Saders, C. E. Fields, Detlev Koester

Abstract: White dwarf (WD) stars evolve simply and predictably, making them reliable age indicators. However, self-consistent validation of the methods for determining WD total ages has yet to be widely performed. This work uses 1565 wide ( > 100 au) WD+WD binaries and 24 new triples containing at least two WDs to test the accuracy and validity of WD total age determinations. For these 1589 wide double-WD b… ▽ More White dwarf (WD) stars evolve simply and predictably, making them reliable age indicators. However, self-consistent validation of the methods for determining WD total ages has yet to be widely performed. This work uses 1565 wide ( > 100 au) WD+WD binaries and 24 new triples containing at least two WDs to test the accuracy and validity of WD total age determinations. For these 1589 wide double-WD binaries and triples, we derive total ages of each WD using photometric data from all-sky surveys, in conjunction with Gaia parallaxes and current hydrogen-atmosphere WD models. Ignoring initial-to-final-mass relations and considering only WD cooling ages, we find that roughly 21-36% of the more massive WDs in a system have a shorter cooling age. Since more massive WDs should be born as more massive main-sequence stars, we attribute this unphysical disagreement as evidence of prior mergers or the presence of an unresolved companion, suggesting that roughly 21-36% of wide WD+WD binaries were once triples. Among the 423 wide WD+WD pairs that pass high-fidelity cuts, we find that 25% total age uncertainties are generally appropriate for WDs with masses > 0.63 Msun and temperatures < 12,000 K, and provide suggested inflation factors for age uncertainties for higher-mass WDs. Overall, WDs return reliable stellar ages, but we detail cases where total ages are least reliable, especially for WDs < 0.63 Msun. △ Less

Submitted 31 May, 2022; originally announced June 2022.

Comments: 24 pages, 10 figures. Submitted to ApJ. Catalog available at https://sites.bu.edu/buwd/files/2022/05/Table_A1.csv

Wide twin binaries are extremely eccentric: evidence of twin binary formation in circumbinary disks

Authors: Hsiang-Chih Hwang, Kareem El-Badry, Hans-Walter Rix, Chris Hamilton, Yuan-Sen Ting, Nadia Zakamska

Abstract: The Gaia mission recently revealed an excess population of equal-mass "twin" wide binaries, with mass ratio $q\gtrsim 0.95$, extending to separations of at least 1000 AU. The origin of this population is an enigma: twin binaries are thought to form via correlated accretion in circumbinary disks, but the typical observed protostellar disks have radii of $\sim100$ AU, far smaller than the separation… ▽ More The Gaia mission recently revealed an excess population of equal-mass "twin" wide binaries, with mass ratio $q\gtrsim 0.95$, extending to separations of at least 1000 AU. The origin of this population is an enigma: twin binaries are thought to form via correlated accretion in circumbinary disks, but the typical observed protostellar disks have radii of $\sim100$ AU, far smaller than the separations of the widest twins. Here, we infer the eccentricity distribution of wide twins from the distribution of their $v$-$r$ angles, i.e., the angle between the components' separation and relative velocity vectors. We find that wide twins must be on extremely eccentric orbits. For the excess-twin population at 400-1000 AU, we infer a near-delta function excess of high-eccentricity system, with eccentricity $0.95 \lesssim e \leq 1$. These high eccentricities for wide twins imply pericenter distances of order $10$ AU and suggest that their orbits were scattered via dynamical interactions in their birth environments, consistent with a scenario in which twins are born in circumbinary disks and subsequently widened. These results further establish twin wide binaries as a distinct population and imply that wide twins can be used as a probe of the dynamical history of stellar populations. △ Less

Submitted 7 July, 2022; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: accepted for publication in The Astrophysical Journal Letter

Formation of proto-globular cluster candidates in cosmological simulations of dwarf galaxies at $z>4$

Authors: Omid Sameie, Michael Boylan-Kolchin, Philip F. Hopkins, Andrew Wetzel, Xiangcheng Ma, James S. Bullock, Kareem El-Badry, Eliot Quataert, Jenna Samuel, Anna T. P. Schauer, Daniel R. Weisz

Abstract: We perform cosmological hydrodynamical simulations to study the formation of proto-globular cluster candidates in progenitors of present-day dwarf galaxies $(M_{\rm vir} \approx 10^{10}\, {\rm M}_\odot$ at $z=0$) as part of the "Feedback in Realistic Environment" (FIRE) project. Compact ($r_{1/2}<30$ pc), relatively massive ($0.5 \times 10^5 \lesssim M_{\star}/{\rm M}_\odot \lesssim 5\times10^5$),… ▽ More We perform cosmological hydrodynamical simulations to study the formation of proto-globular cluster candidates in progenitors of present-day dwarf galaxies $(M_{\rm vir} \approx 10^{10}\, {\rm M}_\odot$ at $z=0$) as part of the "Feedback in Realistic Environment" (FIRE) project. Compact ($r_{1/2}<30$ pc), relatively massive ($0.5 \times 10^5 \lesssim M_{\star}/{\rm M}_\odot \lesssim 5\times10^5$), self-bound stellar clusters form at $11\gtrsim z \gtrsim 5$ in progenitors with $M_{\rm vir} \approx 10^9\,{\rm M}_\odot$. Cluster formation is triggered when at least $10^7\,{\rm M}_\odot$ of dense, turbulent gas reaches $Σ_{\rm gas} \approx 10^4\, {\rm M}_\odot\, {\rm pc}^{-2}$ as a result of the compressive effects of supernova feedback or from cloud-cloud collisions. The clusters can survive for $2-3\,{\rm Gyr}$; absent numerical effects, they would likely survive substantially longer, perhaps to $z=0$. The longest-lived clusters are those that form at significant distance -- several hundreds of pc -- from their host galaxy. We therefore predict that globular clusters forming in progenitors of present-day dwarf galaxies will be offset from any pre-existing stars within their host dark matter halos as opposed to deeply embedded within a well-defined galaxy. Properties of the nascent clusters are consistent with observations of some of the faintest and most compact high-redshift sources in \textit{Hubble Space Telescope} lensing fields and are at the edge of what will be detectable as point sources in deep imaging of non-lensed fields with the \textit{James Webb Space Telescope}. By contrast, the star clusters' host galaxies will remain undetectable. △ Less

Submitted 1 April, 2022; originally announced April 2022.

Comments: 14 pages, 5 figures, submitted to MNRAS

Improving white dwarfs as chronometers with Gaia parallaxes and spectroscopic metallicities

Authors: Adam Moss, Ted von Hippel, Elliot Robinson, Kareem El-Badry, David C. Stenning, David van Dyk, Morgan Fouesneau, Coryn A. L. Bailer-Jones, Elizabeth Jeffery, Jimmy Sargent, Isabelle Kloc, Natalię Moticska

Abstract: White dwarfs (WDs) offer unrealized potential in solving two problems in astrophysics: stellar age accuracy and precision. WD cooling ages can be inferred from surface temperatures and radii, which can be constrained with precision by high quality photometry and parallaxes. Accurate and precise Gaia parallaxes along with photometric surveys provide information to derive cooling and total ages for… ▽ More White dwarfs (WDs) offer unrealized potential in solving two problems in astrophysics: stellar age accuracy and precision. WD cooling ages can be inferred from surface temperatures and radii, which can be constrained with precision by high quality photometry and parallaxes. Accurate and precise Gaia parallaxes along with photometric surveys provide information to derive cooling and total ages for vast numbers of WDs. Here we analyse 1372 WDs found in wide binaries with MS companions, and report on the cooling and total age precision attainable in these WD+MS systems. The total age of a WD can be further constrained if its original metallicity is known, because the main-sequence (MS) progenitor lifetime depends on metallicity at fixed mass, yet metallicity is unavailable via spectroscopy of the WD. We show that incorporating spectroscopic metallicity constraints from 38 wide binary MS companions substantially decreases internal uncertainties in WD total ages compared to a uniform constraint. Averaged over the 38 stars in our sample, the total (internal) age uncertainty improves from 21.04% to 16.77% when incorporating the spectroscopic constraint. Higher mass WDs yield better total age precision; for 8 WDs with zero age main sequence masses $\geq$ 2.0 solar masses, the mean uncertainty in total ages improves from 8.61% to 4.54% when incorporating spectroscopic metallicities. We find that it is often possible to achieve 5% total age precision for WDs with progenitor masses above 2.0 solar masses if parallaxes with $\leq$ 1% precision and Pan-STARRS $g, r$, and $i$ photometry with $\leq$ 0.01 mag precision are available. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Comments: 14 pages, 9 figure, accepted for publication in The Astrophysical Journal

Unicorns and Giraffes in the binary zoo: stripped giants with subgiant companions

Authors: Kareem El-Badry, Rhys Seeburger, Tharindu Jayasinghe, Hans-Walter Rix, Silvia Almada, Charlie Conroy, Adrian M. Price-Whelan, Kevin Burdge

Abstract: We analyze two binary systems containing giant stars, V723 Mon ("the Unicorn") and 2M04123153+6738486 ("the Giraffe"). Both giants orbit more massive but less luminous companions, previously proposed to be mass-gap black holes. Spectral disentangling reveals luminous companions with star-like spectra in both systems. Joint modeling of the spectra, light curves, and spectral energy distributions ro… ▽ More We analyze two binary systems containing giant stars, V723 Mon ("the Unicorn") and 2M04123153+6738486 ("the Giraffe"). Both giants orbit more massive but less luminous companions, previously proposed to be mass-gap black holes. Spectral disentangling reveals luminous companions with star-like spectra in both systems. Joint modeling of the spectra, light curves, and spectral energy distributions robustly constrains the masses, temperatures, and radii of both components: the primaries are luminous, cool giants ($T_{\rm eff,\,giant} = 3,800\,\rm K$ and $4,000\,\rm K$, $R_{\rm giant}= 22.5\,R_{\odot}$ and $25\,R_{\odot}$) with exceptionally low masses ($M_{\rm giant} \approx 0.4\,M_{\odot}$) that likely fill their Roche lobes. The secondaries are only slightly warmer subgiants ($T_{\rm eff,\,2} = 5,800\,\rm K$ and $5,150\,\rm K$, $R_2= 8.3\,R_{\odot}$ and $9\,R_{\odot}$) and thus are consistent with observed UV limits that would rule out main-sequence stars with similar masses ($M_2 \approx 2.8\,M_{\odot}$ and $\approx 1.8\,M_{\odot}$). In the Unicorn, rapid rotation blurs the spectral lines of the subgiant, making it challenging to detect even at wavelengths where it dominates the total light. Both giants have surface abundances indicative of CNO processing and subsequent envelope stripping. The properties of both systems can be reproduced by binary evolution models in which a $1-2\,M_{\odot}$ primary is stripped by a companion as it ascends the giant branch. The fact that the companions are also evolved implies either that the initial mass ratio was very near unity, or that the companions are temporarily inflated due to rapid accretion. The Unicorn and Giraffe offer a window into into a rarely-observed phase of binary evolution preceding the formation of wide-orbit helium white dwarfs, and eventually, compact binaries containing two helium white dwarfs. △ Less

Submitted 21 March, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: 20 pages, 16 figures. Accepted to MNRAS

Public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation

Authors: Andrew Wetzel, Christopher C. Hayward, Robyn E. Sanderson, Xiangcheng Ma, Daniel Angles-Alcazar, Robert Feldmann, T. K Chan, Kareem El-Badry, Coral Wheeler, Shea Garrison-Kimmel, Farnik Nikakhtar, Nondh Panithanpaisal, Arpit Arora, Alexander B. Gurvich, Jenna Samuel, Omid Sameie, Viraj Pandya, Zachary Hafen, Cameron Hummels, Sarah Loebman, Michael Boylan-Kolchin, James S. Bullock, Claude-Andre Faucher-Giguere, Dusan Keres, Eliot Quataert , et al. (1 additional authors not shown)

Abstract: We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation, available at http://flathub.flatironinstitute.org/fire, from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multi-phase interstellar medium while implementing direct models for stellar evolution and feedback, includ… ▽ More We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation, available at http://flathub.flatironinstitute.org/fire, from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multi-phase interstellar medium while implementing direct models for stellar evolution and feedback, including stellar winds, core-collapse and Ia supernovae, radiation pressure, photoionization, and photoelectric heating. We release complete snapshots from 3 suites of simulations. The first comprises 20 simulations that zoom in on 14 Milky Way-mass galaxies, 5 SMC/LMC-mass galaxies, and 4 lower-mass galaxies including 1 ultra-faint; we release 39 snapshots across z = 0 - 10. The second comprises 4 massive galaxies, with 19 snapshots across z = 1 - 10. Finally, a high-redshift suite comprises 22 simulations, with 11 snapshots across z = 5 - 10. Each simulation also includes dozens of resolved lower-mass (satellite) galaxies in its zoom-in region. Snapshots include all stored properties for all dark matter, gas, and star particles, including 11 elemental abundances for stars and gas, and formation times (ages) of star particles. We also release accompanying (sub)halo catalogs, which include galaxy properties and member star particles. For the simulations to z = 0, including all Milky Way-mass galaxies, we release the formation coordinates and an "ex-situ" flag for all star particles, pointers to track particles across snapshots, catalogs of stellar streams, and multipole basis expansions for the halo mass distributions. We describe publicly available python packages for reading and analyzing these simulations. △ Less

Submitted 29 March, 2023; v1 submitted 14 February, 2022; originally announced February 2022.

Comments: 22 pages. Accepted for publication in ApJS, matches the published version, with updated references. Now includes up to 39 snapshots per simulation and formation coordinates for all star particles. Data available at http://flathub.flatironinstitute.org/fire

The "Giraffe": Discovery of a stripped red giant in an interacting binary with a ${\sim}2~M_\odot$ lower giant

Authors: T. Jayasinghe, Todd A. Thompson, C. S. Kochanek, K. Z. Stanek, D. M. Rowan, D. V. Martin, Kareem El-Badry, P. J. Vallely, J. T. Hinkle, D. Huber, H. Isaacson, J. Tayar, K. Auchettl, I. Ilyin, A. W. Howard, C. Badenes

Abstract: We report the discovery of a stripped giant + lower giant binary, 2M04123153+6738486 (2M0412), identified during a search for non-interacting compact object-star binaries. 2M0412 is an evolved ($T_{\rm eff, giant}\simeq4000$ K), luminous ($L_{\rm giant}\simeq150~L_\odot$) red giant in a circular $P=81.2$ day binary. 2M0412 is a known variable star previously classified as a semi-regular variable.… ▽ More We report the discovery of a stripped giant + lower giant binary, 2M04123153+6738486 (2M0412), identified during a search for non-interacting compact object-star binaries. 2M0412 is an evolved ($T_{\rm eff, giant}\simeq4000$ K), luminous ($L_{\rm giant}\simeq150~L_\odot$) red giant in a circular $P=81.2$ day binary. 2M0412 is a known variable star previously classified as a semi-regular variable. The cross-correlation functions of follow-up Keck/HIRES and LBT/PEPSI spectra show an RV-variable second component with implied mass ratio $q=M_{\rm giant}/M_{\rm comp}\simeq0.20\pm0.01$. The ASAS-SN, ATLAS, TESS and ZTF light curves show that the giant is a Roche lobe filling ellipsoidal variable with an inclination of $49.4^\circ{}\pm{0.3^{\circ}}$, and a giant mass of $M_{\rm giant}=0.38\pm0.01~ M_\odot$ for a distance of $\simeq3.7$ kpc. The mass of the giant indicates that its envelope has been stripped. The giant companion on the lower red giant branch has a mass of $M_{\rm comp}=1.91\pm0.03~M_\odot$ with $T_{\rm eff, comp}\simeq5000$ K, $L_{\rm comp}\simeq60~L_\odot$ and $R_{\rm comp}\simeq11~R_\odot$. We also identify an orbital phase dependent, broad $\rm Hα$ emission line which could indicate ongoing accretion from the stripped red giant onto the companion. △ Less

Submitted 15 April, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: Revised version of the paper following comments from the referee, submitted to MNRAS. arXiv admin note: text overlap with arXiv:2101.02212

Hot-mode accretion and the physics of thin-disk galaxy formation

Authors: Zachary Hafen, Jonathan Stern, James Bullock, Alex B. Gurvich, Sijie Yu, Claude-Andre Faucher-Giguere, Drummond B. Fielding, Daniel Angles-Alcazar, Eliot Quataert, Andrew Wetzel, Tjitske Starkenburg, Michael Boylan-Kolchin, Jorge Moreno, Robert Feldmann, Kareem El-Badry, T. K. Chan, Cameron Trapp, Dusan Keres, Philip F. Hopkins

Abstract: We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation *prior* to joining the galaxy. Among galaxies with a high fraction (>70%) of their young stars in a thin disk (h/R~0.1) we find that:… ▽ More We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation *prior* to joining the galaxy. Among galaxies with a high fraction (>70%) of their young stars in a thin disk (h/R~0.1) we find that: (i) hot, virial-temperature gas dominates the inflowing gas mass on halo scales (>~20 kpc), with radiative losses offset by compression heating; (ii) this hot accretion proceeds until angular momentum support slows inward motion, at which point the gas cools to T~10^4 K or less; (iii) prior to cooling, the accreting gas develops an angular momentum distribution that is aligned with the galaxy disk, and while cooling transitions from a quasi-spherical spatial configuration to a more flattened, disk-like configuration. We show that the existence of this "rotating cooling flow" accretion mode is strongly correlated with the fraction of stars forming in a thin disk among a sample of 17 z~0 galaxies spanning a halo mass range of 10^10.5 solar masses to 10^12 solar masses, or a stellar mass range 10^8 solar masses to 10^11 solar masses. Notably, galaxies with a thick disk or irregular morphology do not undergo significant angular momentum alignment of gas prior to accretion and show no correspondence between halo gas cooling and flattening. Our results suggest that rotating cooling flows (or, more generally, rotating subsonic flows) that become coherent and angular momentum-supported prior to accretion onto the galaxy are likely a necessary condition for the formation of thin, star-forming disk galaxies in a LambdaCDM universe. △ Less

Submitted 6 June, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 17 pages, 11 figures. Comments welcome. Accompanying interactive visualization at zhafen.github.io/rotating-cooling-flows

Birth of a Be star: an APOGEE search for Be stars forming through binary mass transfer

Authors: Kareem El-Badry, Charlie Conroy, Eliot Quataert, Hans-Walter Rix, Jonathan Labadie-Bartz, Tharindu Jayasinghe, Todd Thompson, Phillip Cargile, Keivan G. Stassun, Ilya Ilyin

Abstract: Motivated by recent suggestions that many Be stars form through binary mass transfer, we searched the APOGEE survey for Be stars with bloated, stripped companions. From a well-defined parent sample of 297 Be stars, we identified one mass-transfer binary, HD 15124. The object consists of a main-sequence Be star ($M_{\rm Be}=5.3\pm 0.6 \,M_{\odot}$) with a low-mass (… ▽ More Motivated by recent suggestions that many Be stars form through binary mass transfer, we searched the APOGEE survey for Be stars with bloated, stripped companions. From a well-defined parent sample of 297 Be stars, we identified one mass-transfer binary, HD 15124. The object consists of a main-sequence Be star ($M_{\rm Be}=5.3\pm 0.6 \,M_{\odot}$) with a low-mass ($M_{\rm donor}=0.92\pm 0.22\,M_{\odot}$), subgiant companion on a 5.47-day orbit. The emission lines originate in an accretion disk caused by ongoing mass transfer, not from a decretion disk as in classical Be stars. Both stars have surface abundances bearing imprint of CNO processing in the donor's core: the surface helium fraction is $Y_{\rm He}\approx 0.6$, and the nitrogen-to-carbon ratio is 1000 times the solar value. The system's properties are well-matched by binary evolution models in which mass transfer begins while a $3-5\,M_{\odot}$ donor leaves the main sequence, with the secondary becoming the Be star. These models predict that the system will soon become a detached Be + stripped star binary like HR 6819 and LB-1, with the stripped donor eventually contracting to become a core helium-burning sdOB star. Discovery of one object in this short-lived ($\sim$1 Myr) evolutionary phase implies the existence of many more that have already passed through it and are now Be + sdOB binaries. We infer that $(28_{-16}^{+27})\,\%$ of Be stars have stripped companions, most of which are faint. Together with the dearth of main-sequence companions to Be stars and recent discovery of numerous Be + sdOB binaries in the UV, our results imply that binarity plays an important role in the formation of Be stars. △ Less

Submitted 24 August, 2022; v1 submitted 14 January, 2022; originally announced January 2022.

Comments: 25 pages, 19 figures, accepted to MNRAS

NGC 2004 #115: a black hole imposter containing three luminous stars

Authors: Kareem El-Badry, Kevin B. Burdge, Przemek Mróz

Abstract: NGC 2004 #115 is a recently identified black hole (BH) candidate in the Large Magellanic Cloud (LMC) containing a B star orbiting an unseen companion in a 2.9 day orbit and Be star tertiary. We show that the unseen companion is not a $25\,M_{\odot}$ BH, but a $(2-3)\,M_{\odot}$ luminous star. Analyzing the OGLE and MACHO light curves of the system, we detect ellipsoidal variability with amplitude… ▽ More NGC 2004 #115 is a recently identified black hole (BH) candidate in the Large Magellanic Cloud (LMC) containing a B star orbiting an unseen companion in a 2.9 day orbit and Be star tertiary. We show that the unseen companion is not a $25\,M_{\odot}$ BH, but a $(2-3)\,M_{\odot}$ luminous star. Analyzing the OGLE and MACHO light curves of the system, we detect ellipsoidal variability with amplitude 10 times larger than would be expected if the companion were a $25\,M_{\odot}$ BH, ruling out the low inclination required for a massive companion. The light curve also shows a clear reflection effect that is well-modeled with a $2.5\,M_{\odot}$ main-sequence secondary, ruling out a lower-mass BH or neutron star companion. We consider and reject models in which the system is a binary containing a stripped star orbiting the Be star: only a triple model with an outer Be star can explain both the observed light curve and radial velocities. Our results imply that the B star, whose slow projected rotation velocity and presumed tidal synchronization were interpreted as evidence for a low inclination (and thus a high companion mass), is far from being tidally synchronized: despite being in a 2.9 day orbit that is fully or nearly circularized ($e < 0.04$), its surface rotation period appears to be at least 20 days. We offer cautionary notes on the interpretation of dormant BH candidates in binaries. △ Less

Submitted 27 January, 2022; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: Accepted to MNRAS

Science with the Ultraviolet Explorer (UVEX)

Authors: S. R. Kulkarni, Fiona A. Harrison, Brian W. Grefenstette, Hannah P. Earnshaw, Igor Andreoni, Danielle A. Berg, Joshua S. Bloom, S. Bradley Cenko, Ryan Chornock, Jessie L. Christiansen, Michael W. Coughlin, Alexander Wuollet Criswell, Behnam Darvish, Kaustav K. Das, Kishalay De, Luc Dessart, Don Dixon, Bas Dorsman, Kareem El-Badry, Christopher Evans, K. E. Saavik Ford, Christoffer Fremling, Boris T. Gansicke, Suvi Gezari, Y. Goetberg , et al. (31 additional authors not shown)

Abstract: UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky… ▽ More UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky survey 50/100 times deeper than GALEX in the NUV/FUV, cadenced surveys of the Large and Small Magellanic Clouds, rapid target of opportunity followup, as well as spectroscopic followup of samples of stars and galaxies. The science program is built around three pillars. First, UVEX will explore the low-mass, low-metallicity galaxy frontier through imaging and spectroscopic surveys that will probe key aspects of the evolution of galaxies by understanding how star formation and stellar evolution at low metallicities affect the growth and evolution of low-metallicity, low-mass galaxies in the local universe. Such galaxies contain half the mass in the local universe, and are analogs for the first galaxies, but observed at distances that make them accessible to detailed study. Second, UVEX will explore the dynamic universe through time-domain surveys and prompt spectroscopic followup capability will probe the environments, energetics, and emission processes in the early aftermaths of gravitational wave-discovered compact object mergers, discover hot, fast UV transients, and diagnose the early stages of stellar explosions. Finally, UVEX will become a key community resource by leaving a large all-sky legacy data set, enabling a wide range of scientific studies and filling a gap in the new generation of wide-field, sensitive optical and infrared surveys provided by the Rubin, Euclid, and Roman observatories. This paper discusses the scientific potential of UVEX, and the broad scientific program. △ Less

Submitted 17 January, 2023; v1 submitted 30 November, 2021; originally announced November 2021.

Comments: 69 pages, 43 figures

NGC 1850 BH1 is another stripped-star binary masquerading as a black hole

Authors: Kareem El-Badry, Kevin Burdge

Abstract: We show that the radial velocity-variable star in the black hole candidate NGC 1850 BH1 cannot be a normal $\approx 5\,M_{\odot}$ subgiant, as was proposed, but is an overluminous stripped-envelope star with mass $\approx 1 M_{\odot}$. The result follows directly from the star's observed radius and the orbital period -- density relation for Roche lobe-filling stars: the star's density, as constrai… ▽ More We show that the radial velocity-variable star in the black hole candidate NGC 1850 BH1 cannot be a normal $\approx 5\,M_{\odot}$ subgiant, as was proposed, but is an overluminous stripped-envelope star with mass $\approx 1 M_{\odot}$. The result follows directly from the star's observed radius and the orbital period -- density relation for Roche lobe-filling stars: the star's density, as constrained by the observed ellipsoidal variability, is too low for its mass to exceed $\approx 1.5\,M_{\odot}$. This lower mass significantly reduces the implied mass of the unseen companion and qualitative interpretation of the system, such that a normal main-sequence companion with mass $(2.5-5)\,M_{\odot}$ is fully consistent with the data. We explore evolutionary scenarios that could produce the binary using MESA and find that its properties can be matched by models in which a $\sim5\,M_{\odot}$ primary loses most of its envelope to a companion and is observed in a bloated state before contracting to become a core helium burning sdOB star. This is similar to the scenario proposed to explain the binaries LB-1 and HR 6819. Though it likely does not contain a black hole, NGC 1850 BH1 provides an interesting test case for binary evolution models, particularly given its membership in a cluster of known age. △ Less

Submitted 29 December, 2021; v1 submitted 15 November, 2021; originally announced November 2021.

Comments: 5 pages, 3 figures. Accepted to MNRAS Letters

RR Lyrae-based Distances for 39 Nearby Dwarf Galaxies Calibrated to Gaia eDR3

Authors: Pranav Nagarajan, Daniel R. Weisz, Kareem El-Badry

Abstract: We provide uniform RR Lyrae-based distances to 39 dwarf galaxies in and around the Local Group. We determine distances based on a Bayesian hierarchical model that uses periods and magnitudes of published RR Lyrae in dwarf galaxies and is anchored to well-studied Milky Way (MW) RR Lyrae with spectroscopic metallicities and Gaia eDR3 parallaxes. Gaia eDR3 parallaxes for the anchor sample are a facto… ▽ More We provide uniform RR Lyrae-based distances to 39 dwarf galaxies in and around the Local Group. We determine distances based on a Bayesian hierarchical model that uses periods and magnitudes of published RR Lyrae in dwarf galaxies and is anchored to well-studied Milky Way (MW) RR Lyrae with spectroscopic metallicities and Gaia eDR3 parallaxes. Gaia eDR3 parallaxes for the anchor sample are a factor of 2, on average, more precise than DR2 parallaxes, and allow for a much better constrained period-luminosity-metallicity relation. While $\sim75$% of our distances are within 1-$σ$ of recent literature RR Lyrae distances, our distances are also $\sim2$-$3$ times more precise than literature distances, on average. On average, our distances are $\sim0.05$ mag closer than literature distances, as well as $\sim0.06$ mag closer than distances derived using a theoretical period-luminosity-metallicity relation. These discrepancies are largely due to our eDR3 parallax anchor. We show that the Hipparcos-anchored RR Lyrae distance scale of Carretta et al. (2000) over-predicts distances to MW RR Lyrae by $\sim0.05$ mag. The largest uncertainties in our distances are (i) the lack of direct metallicity measurements for RR Lyrae and (ii) the heterogeneity of published RR Lyrae photometry. We provide simple formulae to place new dwarf galaxies with RR Lyrae on a common distance scale with this work. We provide public code that can easily incorporate additional galaxies and data from future surveys, providing a versatile framework for cartography of the local Universe with RR Lyrae. △ Less

Submitted 16 June, 2022; v1 submitted 12 November, 2021; originally announced November 2021.

Comments: 34 pages, 15 figures, Published in ApJ

Wide binaries from the H3 survey: the thick disk and halo have similar wide binary fractions

Authors: Hsiang-Chih Hwang, Yuan-Sen Ting, Charlie Conroy, Nadia L. Zakamska, Kareem El-Badry, Phillip Cargile, Dennis Zaritsky, Vedant Chandra, Jiwon Jesse Han, Joshua S. Speagle, Ana Bonaca

Abstract: Due to the different environments in the Milky Way's disk and halo, comparing wide binaries in the disk and halo is key to understanding wide binary formation and evolution. By using Gaia Early Data Release 3, we search for resolved wide binary companions in the H3 survey, a spectroscopic survey that has compiled $\sim$150,000 spectra for thick-disk and halo stars to date. We identify 800 high-con… ▽ More Due to the different environments in the Milky Way's disk and halo, comparing wide binaries in the disk and halo is key to understanding wide binary formation and evolution. By using Gaia Early Data Release 3, we search for resolved wide binary companions in the H3 survey, a spectroscopic survey that has compiled $\sim$150,000 spectra for thick-disk and halo stars to date. We identify 800 high-confidence (a contamination rate of 4%) wide binaries and two resolved triples, with binary separations mostly between $10^3$-$10^5$ AU and a lowest [Fe/H] of $-2.7$. Based on their Galactic kinematics, 33 of them are halo wide binaries, and most of those are associated with the accreted Gaia-Sausage-Enceladus galaxy. The wide binary fraction in the thick disk decreases toward the low metallicity end, consistent with the previous findings for the thin disk. Our key finding is that the halo wide binary fraction is consistent with the thick-disk stars at a fixed [Fe/H]. There is no significant dependence of the wide binary fraction on the $α$-captured abundance. Therefore, the wide binary fraction is mainly determined by the iron abundance, not their disk or halo origin nor the $α$-captured abundance. Our results suggest that the formation environments play a major role for the wide binary fraction, instead of other processes like radial migration that only apply to disk stars. △ Less

Submitted 7 March, 2022; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS. Fig. 7 is the key plot. Online data are available at: https://drive.google.com/file/d/1TXKEwKBHcP4j9tMeJpa2hWDFHsEc1Qd3/view?usp=sharing. Comments are welcome!

Predictions for Complex Distributions of Stellar Elemental Abundances in Low-Mass Galaxies

Authors: Preet B. Patel, Sarah R. Loebman, Andrew Wetzel, Claude-André Faucher-Giguère, Kareem El-Badry, Jeremy Bailin

Abstract: We investigate stellar elemental abundance patterns at z = 0 in 8 low-mass (M_* = 10^6 - 10^9 M_sun) galaxies in the Feedback in Realistic Environments (FIRE-2) cosmological simulations. Using magnesium (Mg) as a representative alpha-element, we explore stellar abundance patterns in [Mg/Fe] versus [Fe/H], which follow an overall monotonic trend that evolved slowly over time. Beyond this, we explor… ▽ More We investigate stellar elemental abundance patterns at z = 0 in 8 low-mass (M_* = 10^6 - 10^9 M_sun) galaxies in the Feedback in Realistic Environments (FIRE-2) cosmological simulations. Using magnesium (Mg) as a representative alpha-element, we explore stellar abundance patterns in [Mg/Fe] versus [Fe/H], which follow an overall monotonic trend that evolved slowly over time. Beyond this, we explore 3 notable secondary features in enrichment (found in three different case-study galaxies) that arise from a galaxy merger or bursty star formation. First, we observe a secondary track with a lower [Mg/Fe] than the main trend. At z = 0, stars from this track are predominantly found within 2-6 kpc of the center; they were accreted in a 1:3 total-mass-ratio merger ~ 0.4 Gyr ago. Second, we find a distinct elemental bi-modality that forms following a strong burst in star formation in a galaxy at t_lookback ~ 10 Gyr. This burst quenched star formation for ~ 0.66 Gyr, allowing Ia supernovae to enrich the system with iron before star formation resumed. Third, we examine stripes in enrichment that run roughly orthogonal to the dominant [Mg/Fe] versus [Fe/H] trend; these stripes correspond to short bursts of star formation during which core-collapse supernovae enrich the surrounding medium with Mg (and Fe) on short timescales. If observed, these features would substantiate the utility of elemental abundances in revealing the assembly and star formation histories of dwarf galaxies. We explore the observability of these features for upcoming spectroscopic studies. Our results show that precise measurements of elemental abundance patterns can reveal critical events in the formation histories of low-mass galaxies. △ Less

Submitted 23 March, 2022; v1 submitted 15 October, 2021; originally announced October 2021.

Comments: v1: 14 pages, 8 figures. Submitted to MNRAS on Oct. 18th, 2021 | v2: 15 pages, 8 figures. Accepted by MNRAS on Mar. 17th, 2022

The In-situ Origins of Dwarf Stellar Outskirts in FIRE-2

Authors: Erin Kado-Fong, Robyn E. Sanderson, Jenny E. Greene, Emily C. Cunningham, Coral Wheeler, T. K. Chan, Kareem El-Badry, Philip F. Hopkins, Andrew Wetzel, Michael Boylan-Kolchin, Claude-André Faucher-Giguère, Song Huang, Eliot Quataert, Tjitske Starkenburg

Abstract: Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies ($10^{8.5}<M_\star/M_\odot<10^{9.6}$) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosm… ▽ More Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies ($10^{8.5}<M_\star/M_\odot<10^{9.6}$) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations are capable of producing dwarf galaxies with realistic structure, including both a thick disk and round stellar halo. Crucially, these stellar halos are formed in-situ, largely via the outward migration of disk stars. However, there also exists a large population of "non-disky" dwarfs in FIRE-2 that lack a well-defined disk/halo and do not resemble the observed dwarf population. These non-disky dwarfs tend to be either more gas poor or to have burstier recent star formation histories than the disky dwarfs, suggesting that star formation feedback may be preventing disk formation. Both classes of dwarfs underscore the power of a galaxy's intrinsic shape -- which is a direct quantification of the distribution of the galaxy's stellar content -- to interrogate the feedback implementation in simulated galaxies. △ Less

Submitted 2 May, 2022; v1 submitted 10 September, 2021; originally announced September 2021.

Comments: 28 pages, 15 figures. Accepted to ApJ

Testing the Relationship Between Bursty Star Formation and Size Fluctuations of Local Dwarf Galaxies

Authors: Najmeh Emami, Brian Siana, Kareem El-Badry, David Cook, Xiangcheng Ma, Daniel Weisz, Joobin Gharibshah, Sara Alaee, Claudia Scarlata, Evan Skillman

Abstract: Stellar feedback in dwarf galaxies plays a critical role in regulating star formation via galaxy-scale winds. Recent hydrodynamical zoom simulations of dwarf galaxies predict that the periodic outward flow of gas can change the gravitational potential sufficiently to cause radial migration of stars. To test the effect of bursty star formation on stellar migration, we examine star formation observa… ▽ More Stellar feedback in dwarf galaxies plays a critical role in regulating star formation via galaxy-scale winds. Recent hydrodynamical zoom simulations of dwarf galaxies predict that the periodic outward flow of gas can change the gravitational potential sufficiently to cause radial migration of stars. To test the effect of bursty star formation on stellar migration, we examine star formation observables and sizes of 86 local dwarf galaxies. We find a correlation between the R-band half-light radius (R$_e$) and far-UV luminosity (L$_{FUV}$) for stellar masses below 10$^8 $ M$_{\odot}$ and a weak correlation between the R$_e$ and H$α$ luminosity (L$_{Hα}$). We produce mock observations of eight low-mass galaxies from the FIRE-2 cosmological simulations and measure the similarity of the time sequences of R$_e$ and a number of star formation indicators with different timescales. Major episodes of R$_e$ time sequence align very well with the major episodes of star formation, with a delay of $\sim$ 50 Myrs. This correlation decreases towards SFR indicators of shorter timescales such that $R_e$ is weakly correlated with L$_{FUV}$ (10-100 Myr timescale) and is completely uncorrelated with L$_{Hα}$ (a few Myr timescale), in agreement with the observations. Our findings based on FIRE-2 suggest that the R-band size of a galaxy reacts to star formation variations on a $\sim50$ Myr timescale. With the advent of a new generation of large space telescopes (e.g., JWST), this effect can be examined explicitly in galaxies at higher redshifts where bursty star formation is more prominent. △ Less

Submitted 19 August, 2021; originally announced August 2021.

Comments: 17 pages, 11 figures, Accepted to ApJ

Birth of the ELMs: a ZTF survey for evolved cataclysmic variables turning into extremely low-mass white dwarfs

Authors: Kareem El-Badry, Hans-Walter Rix, Eliot Quataert, Thomas Kupfer, Ken J. Shen

Abstract: We present a systematic survey for mass-transferring and recently-detached cataclysmic variables (CVs) with evolved secondaries, which are progenitors of extremely low-mass white dwarfs (ELM WDs), AM CVn systems, and detached ultracompact binaries. We select targets below the main sequence in the Gaia color-magnitude diagram with ZTF light curves showing large-amplitude ellipsoidal variability and… ▽ More We present a systematic survey for mass-transferring and recently-detached cataclysmic variables (CVs) with evolved secondaries, which are progenitors of extremely low-mass white dwarfs (ELM WDs), AM CVn systems, and detached ultracompact binaries. We select targets below the main sequence in the Gaia color-magnitude diagram with ZTF light curves showing large-amplitude ellipsoidal variability and orbital period $P_{\rm orb} < 6$ hr. This yields 51 candidates brighter than G=18, of which we have obtained many-epoch spectra for 21. We confirm all 21 to be completely -- or nearly -- Roche lobe filling close binaries. 13 show evidence of ongoing mass transfer, which has likely just ceased in the other 8. Most of the secondaries are hotter than any previously known CV donors, with temperatures $4700<T_{\rm eff}/{\rm K}<8000$. Remarkably, all secondaries with $T_{\rm eff} \gtrsim 7000\,\rm K$ appear to be detached, while all cooler secondaries are still mass-transferring. This transition likely marks the temperature where magnetic braking becomes inefficient due to loss of the donor's convective envelope. Most of the proto-WD secondaries have masses near $0.15\,M_{\odot}$; their companions have masses near $0.8\,M_{\odot}$. We infer a space density of $\sim 60\,\rm kpc^{-3}$, roughly 80 times lower than that of normal CVs and three times lower than that of ELM WDs. The implied Galactic birth rate, $\mathcal{R}\sim 60\,\rm Myr^{-1}$, is half that of AM CVn binaries. Most systems are well-described by MESA models for CVs in which mass transfer begins only as the donor leaves the main sequence. All are predicted to reach minimum periods $5\lesssim P_{\rm orb}/{\rm min}\lesssim30$ within a Hubble time, where they will become AM CVn binaries or merge. This sample triples the known evolved CV population and offers broad opportunities for improving understanding of the compact binary population. △ Less

Submitted 30 November, 2021; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: 26 pages, 16 figures, plus appendices. Magnitudes in Table 1 fixed; machine-readable table included in source files. Accepted to MNRAS

The effects of LMC-mass environments on their dwarf satellite galaxies in the FIRE simulations

Authors: Ethan D. Jahn, Laura V. Sales, Andrew Wetzel, Jenna Samuel, Kareem El-Badry, Michael Boylan-Kolchin, James S. Bullock

Abstract: Characterizing the predicted environments of dwarf galaxies like the Large Magellanic Cloud (LMC) is becoming increasingly important as next generation surveys push sensitivity limits into this low-mass regime at cosmological distances. We study the environmental effects of LMC-mass halos ($M_{200m} \sim 10^{11}$ M$_\odot$) on their populations of satellites ($M_\star \geq 10^4$ M$_\odot$) using a… ▽ More Characterizing the predicted environments of dwarf galaxies like the Large Magellanic Cloud (LMC) is becoming increasingly important as next generation surveys push sensitivity limits into this low-mass regime at cosmological distances. We study the environmental effects of LMC-mass halos ($M_{200m} \sim 10^{11}$ M$_\odot$) on their populations of satellites ($M_\star \geq 10^4$ M$_\odot$) using a suite of zoom-in simulations from the Feedback In Realistic Environments (FIRE) project. Our simulations predict significant hot coronas with $T\sim10^5$ K and $M_\text{gas}\sim10^{9.5}$ M$_\odot$. We identify signatures of environmental quenching in dwarf satellite galaxies, particularly for satellites with intermediate mass ($M_\star = 10^{6-7}$ M$_\odot$). The gas content of such objects indicates ram-pressure as the likely quenching mechanism, sometimes aided by star formation feedback. Satellites of LMC-mass hosts replicate the stellar mass dependence of the quiescent fraction found in satellites of MW mass hosts (i.e. that the quiescent fraction increases as stellar mass decreases). Satellites of LMC-mass hosts have a wider variety of quenching times when compared to the strongly bi-modal distribution of quenching times of nearby centrals. Finally, we identify significant tidal stellar structures around four of our six LMC-analogs, suggesting that stellar streams may be common. These tidal features originated from satellites on close orbits, extend to $\sim$80 kpc from the central galaxy, and contain $\sim10^{6-7}$ M$_\odot$~of stars. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: 14 pages, 10 figures, submitted to MNRAS

The Early Phases of Supernova 2020pni: Shock-Ionization of the Nitrogen-Enriched Circumstellar Material

Authors: G. Terreran, W. V. Jacobson-Galan, J. H. Groh, R. Margutti, D. L. Coppejans, G. Dimitriadis, C. D. Kilpatrick, D. J. Matthews, M. R. Siebert, C. R. Angus, T. G. Brink, A. V. Filippenko, R. J. Foley, D. O. Jones, S. Tinyanont, C. Gall, H. Pfister, Y. Zenati, Z. Ansari, K. Auchettl, K. El-Badry, E. A. Magnier, W. Zheng

Abstract: We present multiwavelength observations of the Type II SN 2020pni. Classified at $\sim 1.3$ days after explosion, the object showed narrow (FWHM $<250\,\textrm{km}\,\textrm{s}^{-1}$) recombination lines of ionized helium, nitrogen, and carbon, as typically seen in flash-spectroscopy events. Using the non-LTE radiative transfer code CMFGEN to model our first high resolution spectrum, we infer a pro… ▽ More We present multiwavelength observations of the Type II SN 2020pni. Classified at $\sim 1.3$ days after explosion, the object showed narrow (FWHM $<250\,\textrm{km}\,\textrm{s}^{-1}$) recombination lines of ionized helium, nitrogen, and carbon, as typically seen in flash-spectroscopy events. Using the non-LTE radiative transfer code CMFGEN to model our first high resolution spectrum, we infer a progenitor mass-loss rate of $\dot{M}=(3.5-5.3)\times10^{-3}$ M$_{\odot}$ yr$^{-1}$ (assuming a wind velocity of $v_w=200\,\textrm{km}\,\textrm{s}^{-1}$), estimated at a radius of $R_{\rm in}=2.5\times10^{14}\,\rm{cm}$. In addition, we find that the progenitor of SN 2020pni was enriched in helium and nitrogen (relative abundances in mass fractions of 0.30$-$0.40, and $8.2\times10^{-3}$, respectively). Radio upper limits are also consistent with a dense CSM, and a mass-loss rate of $\dot M>5 \times 10^{-4}\,\rm{M_{\odot}\,yr^{-1}}$. During the first 4 days after first light, we also observe an increase in velocity of the hydrogen lines (from $\sim 250\,\textrm{km}\,\textrm{s}^{-1}$ to $\sim 1000\,\textrm{km}\,\textrm{s}^{-1}$), which suggests a complex CSM. The presence of dense and confined CSM, as well as its inhomogeneous structure, suggest a phase of enhanced mass loss of the progenitor of SN 2020pni during the last year before explosion. Finally, we compare SN 2020pni to a sample of other shock-photoionization events. We find no evidence of correlations among the physical parameters of the explosions and the characteristics of the CSM surrounding the progenitors of these events. This favors the idea that the mass-loss experienced by massive stars during their final years could be governed by stochastic phenomena, and that, at the same time, the physical mechanisms responsible for this mass-loss must be common to a variety of different progenitors. △ Less

Submitted 1 March, 2022; v1 submitted 25 May, 2021; originally announced May 2021.

Comments: 19 pages, 14 figures

Neutral CGM as damped Lyα absorbers at high redshift

Authors: Jonathan Stern, Amiel Sternberg, Claude-André Faucher-Giguère, Zachary Hafen, Drummond Fielding, Eliot Quataert, Andrew Wetzel, Daniel Anglés-Alcázar, Kareem El-Badry, Dušan Kereš, Philip F. Hopkins

Abstract: Recent searches for the hosts of high-redshift ($z \sim 4$) damped Ly$α$ absorbers (DLAs) have detected bright galaxies at distances of tens of kpc from the DLA. Using the FIRE-2 cosmological zoom simulations, we argue that these relatively large distances are due to a predominantly cool and neutral inner circumgalactic medium (CGM) surrounding high-redshift galaxies. The inner CGM is cool because… ▽ More Recent searches for the hosts of high-redshift ($z \sim 4$) damped Ly$α$ absorbers (DLAs) have detected bright galaxies at distances of tens of kpc from the DLA. Using the FIRE-2 cosmological zoom simulations, we argue that these relatively large distances are due to a predominantly cool and neutral inner circumgalactic medium (CGM) surrounding high-redshift galaxies. The inner CGM is cool because of the short cooling time of hot gas in $\lesssim10^{12}$ Msun halos, which implies that accretion and feedback energy are radiated quickly, while it is neutral due to the high volume densities and column densities at high redshift which shield cool gas from photoionization. Our analysis predicts large DLA covering factors ($\gtrsim50\%$) out to impact parameters $\sim0.3((1 + z)/5)^{3/2}\ R_{\rm vir}$ from the central galaxies at $z > 1$, equivalent to a physical distance of $\sim 21 M_{12}^{1/3} ((1 + z)/5)^{1/2}$ kpc ($R_{\rm vir}$ and $M_{12}$ are the halo virial radius and mass in units of $10^{12}$ Msun, respectively). This implies that DLA covering factors at $z \sim 4$ may be comparable to unity out to a distance $\sim 10$ times larger than stellar half-mass radii. A predominantly neutral inner CGM in the early universe suggests that its mass and metallicity can be directly constrained by CGM absorption surveys, without resorting to large ionization corrections as required for ionized CGM. △ Less

Submitted 30 July, 2021; v1 submitted 13 May, 2021; originally announced May 2021.

Comments: 16 pages, 13 figures, accepted to MNRAS

Distant Relatives: The Chemical Homogeneity of Comoving Pairs Identified in Gaia

Authors: Tyler Nelson, Yuan-Sen Ting, Keith Hawkins, Alexander Ji, Harshil Kamdar, Kareem El-Badry

Abstract: Comoving pairs, even at the separations of $\mathcal{O}(10^6)\,$AU, are a predicted reservoir of conatal stars. We present detailed chemical abundances of 62 stars in 31 comoving pairs with separations of $10^2 - 10^7\,$AU and 3D velocity differences $< 2 \mathrm{\ km \ s^{-1}}$. This sample includes both bound comoving pairs/wide binaries and unbound comoving pairs. Observations were taken using… ▽ More Comoving pairs, even at the separations of $\mathcal{O}(10^6)\,$AU, are a predicted reservoir of conatal stars. We present detailed chemical abundances of 62 stars in 31 comoving pairs with separations of $10^2 - 10^7\,$AU and 3D velocity differences $< 2 \mathrm{\ km \ s^{-1}}$. This sample includes both bound comoving pairs/wide binaries and unbound comoving pairs. Observations were taken using the MIKE spectrograph on the Magellan/Clay Telescope at high resolution ($\mathrm{R} \sim 45,000$) with a typical signal-to-noise ratio of 150 per pixel. With these spectra, we measure surface abundances for 24 elements, including Li, C, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Sr, Y, Zr, Ba, La, Nd, Eu. Taking iron as the representative element, our sample of wide binaries is chemically homogeneous at the level of $0.05$ dex, which agrees with prior studies on wide binaries. Importantly, even systems at separations $2\times10^5-10^7\,$AU are homogeneous to $0.09$ dex, as opposed to the random pairs which have a dispersion of $0.23\,$dex. Assuming a mixture model of the wide binaries and random pairs, we find that $73 \pm 22\%$ of the comoving pairs at separations $2\times10^5-10^7\,$AU are conatal. Our results imply that a much larger parameter space of phase space may be used to find conatal stars, to study M-dwarfs, star cluster evolution, exoplanets, chemical tagging, and beyond. △ Less

Submitted 15 June, 2021; v1 submitted 26 April, 2021; originally announced April 2021.

Comments: 16 pages, 8 figures, submitted to ApJ

LAMOST J0140355+392651: An evolved cataclysmic variable donor transitioning to become an extremely low mass white dwarf

Authors: Kareem El-Badry, Eliot Quataert, Hans-Walter Rix, Daniel R. Weisz, Thomas Kupfer, Ken J. Shen, Maosheng Xiang, Yong Yang, Xiaowei Liu

Abstract: We present LAMOST J0140355+392651 (hereafter J0140), a close ($P_{\rm orb} = 3.81$ hours) binary containing a bloated, low-mass ($M \approx 0.15 M_{\odot}$) proto-white dwarf (WD) and a massive ($M\approx 0.95\,M_{\odot}$) WD companion. The system's optical light curve is dominated by large-amplitude ellipsoidal variability but also exhibits additional scatter, likely driven by pulsations. The pro… ▽ More We present LAMOST J0140355+392651 (hereafter J0140), a close ($P_{\rm orb} = 3.81$ hours) binary containing a bloated, low-mass ($M \approx 0.15 M_{\odot}$) proto-white dwarf (WD) and a massive ($M\approx 0.95\,M_{\odot}$) WD companion. The system's optical light curve is dominated by large-amplitude ellipsoidal variability but also exhibits additional scatter, likely driven by pulsations. The proto-WD is cooler ($T_{\rm eff} = 6800\pm 100$ K) and more puffy ($\log\left[g/\left({\rm cm\,s^{-2}}\right)\right]=4.74\pm0.07$) than any known extremely low mass (ELM) WD, but hotter than any known cataclysmic variable (CV) donor. It either completely or very nearly fills its Roche lobe ($R/R_{\rm Roche\,lobe}=0.99\pm0.01$), suggesting ongoing or recently terminated mass transfer. No dwarf nova-like outbursts have been observed. The spectrum is dominated by the proto-WD but shows tentative hints of H$α$ emission, perhaps due to accretion onto the massive WD. The properties of the system are well-matched by MESA binary evolution models of CVs with donors that underwent significant nuclear evolution before the onset of mass transfer. In these models, the bloated proto-WD is either still losing mass via stable Roche lobe overflow or was doing so until very recently. In either case, it is evolving toward higher temperatures at near-constant luminosity to become an ELM WD. If the system is detached, mass transfer likely ended when the donor became too hot for magnetic braking to remain efficient. Evolutionary models predict that the binary will shrink to $P_{\rm orb}\lesssim 10$ minutes within a few Gyr, when it will either merge or become an AM CVn binary. J0140 provides an observational link between the formation channels of CVs, ELM WDs, detached ultracompact WD binaries, and AM CVn systems. △ Less

Submitted 6 May, 2021; v1 submitted 14 April, 2021; originally announced April 2021.

Comments: 19 pages, 14 figures. Accepted to MNRAS

The origin of metal-poor stars on prograde disk orbits in FIRE simulations of Milky Way-mass galaxies

Authors: Isaiah B. Santistevan, Andrew Wetzel, Robyn E. Sanderson, Kareem El-Badry, Jenna Samuel, Claude-André Faucher-Giguère

Abstract: In hierarchical structure formation, metal-poor stars in and around the Milky Way (MW) originate primarily from mergers of lower-mass galaxies. A common expectation is therefore that metal-poor stars should have isotropic, dispersion-dominated orbits that do not correlate strongly with the MW disk. However, recent observations of stars in the MW show that metal-poor ([Fe/H] < -2) stars are prefere… ▽ More In hierarchical structure formation, metal-poor stars in and around the Milky Way (MW) originate primarily from mergers of lower-mass galaxies. A common expectation is therefore that metal-poor stars should have isotropic, dispersion-dominated orbits that do not correlate strongly with the MW disk. However, recent observations of stars in the MW show that metal-poor ([Fe/H] < -2) stars are preferentially on prograde orbits with respect to the disk. Using the FIRE-2 suite of cosmological zoom-in simulations of MW/M31-mass galaxies, we investigate the prevalence and origin of prograde metal-poor stars. Almost all (11 of 12) of our simulations have metal-poor stars on preferentially prograde orbits today and throughout most of their history: we thus predict that this is a generic feature of MW/M31-mass galaxies. The typical prograde-to-retrograde ratio is ~2:1, which depends weakly on stellar metallicity at Fe/H] < -1. These trends predicted by our simulations agree well with MW observations. Prograde metal-poor stars originate largely from a single LMC/SMC-mass gas-rich merger 7 - 12.5 Gyr ago, which deposited existing metal-poor stars and significant gas on an orbital vector that sparked the formation of and/or shaped the orientation of a long-lived stellar disk, giving rise to a prograde bias for all low-metallicity stars. We find sub-dominant contributions from in-situ stars formed in the host galaxy before this merger, and in some cases, additional massive mergers. We find few clear correlations between any properties of our MW/M31-mass galaxies at z = 0 and the degree of this prograde bias as a result of diverse merger scenarios. △ Less

Submitted 7 May, 2021; v1 submitted 5 February, 2021; originally announced February 2021.

Comments: 18 pages, 8 figures. Accepted for publication in MNRAS

A classifier for spurious astrometric solutions in Gaia EDR3

Authors: Jan Rybizki, Gregory Green, Hans-Walter Rix, Kareem El-Badry, Markus Demleitner, Eleonora Zari, Andrzej Udalski, Richard L. Smart, Andrew Gould

Abstract: The Gaia early Data Release 3 has delivered exquisite astrometric data for 1.47 billion sources, which is revolutionizing many fields in astronomy. For a small fraction of these sources, the astrometric solutions are poor, and the reported values and uncertainties may not apply. Before any analysis, it is important to recognize and excise these spurious results - this is commonly done by means of… ▽ More The Gaia early Data Release 3 has delivered exquisite astrometric data for 1.47 billion sources, which is revolutionizing many fields in astronomy. For a small fraction of these sources, the astrometric solutions are poor, and the reported values and uncertainties may not apply. Before any analysis, it is important to recognize and excise these spurious results - this is commonly done by means of quality flags in the Gaia catalog. Here, we devise a means of separating 'good' from 'bad' astrometric solutions that is an order of magnitude cleaner than any single flag: 99.3% pure and 97.3% complete, as validated on our test data. We devise an extensive sample of manifestly bad astrometric solutions, with parallax that is negative at > 4.5 sigma; and a corresponding sample of presumably good solutions, including sources in HEALPix pixels on the sky that do not contain such negative parallaxes, and sources that fall on the main sequence in a color-absolute magnitude diagram. We then train a neural network that uses 17 pertinent Gaia catalog entries and information about nearby sources to discriminate between these two samples, captured in a single 'astrometric fidelity' parameter. A diverse set of verification tests shows that our approach works very cleanly, including for sources with positive parallaxes. The main limitations of our approach are in the very low-SNR and the crowded regime. Our astrometric fidelities for all of eDR3 can be queried via the Virtual Observatory, our code and data are public. △ Less

Submitted 7 December, 2021; v1 submitted 27 January, 2021; originally announced January 2021.

Comments: 21 pages, 20 figures, accepted with MNRAS. Fidelities available via Virtual Observatory (GAVO). Catalog information: https://dc.zah.uni-heidelberg.de/tableinfo/gedr3spur.main -- Colab python notebook: https://colab.research.google.com/drive/1lPzhGSSIjx2nQ7XM2v8bQZtkf0Atrk0z?usp=sharing

A million binaries from Gaia eDR3: sample selection and validation of Gaia parallax uncertainties

Authors: Kareem El-Badry, Hans-Walter Rix, Tyler M. Heintz

Abstract: We construct from Gaia eDR3 an extensive catalog of spatially resolved binary stars within $\approx$ 1 kpc of the Sun, with projected separations ranging from a few au to 1 pc. We estimate the probability that each pair is a chance alignment empirically, using the Gaia catalog itself to calculate the rate of chance alignments as a function of observables. The catalog contains 1.3 (1.1) million bin… ▽ More We construct from Gaia eDR3 an extensive catalog of spatially resolved binary stars within $\approx$ 1 kpc of the Sun, with projected separations ranging from a few au to 1 pc. We estimate the probability that each pair is a chance alignment empirically, using the Gaia catalog itself to calculate the rate of chance alignments as a function of observables. The catalog contains 1.3 (1.1) million binaries with >90% (>99%) probability of being bound, including 16,000 white dwarf -- main sequence (WD+MS) binaries and 1,400 WD+WD binaries. We make the full catalog publicly available, as well as the queries and code to produce it. We then use this sample to calibrate the published Gaia DR3 parallax uncertainties, making use of the binary components' near-identical parallaxes. We show that these uncertainties are generally reliable for faint stars ($G\gtrsim 18$), but are underestimated significantly for brighter stars. The underestimates are generally $\le 30\%$ for isolated sources with well-behaved astrometry, but are larger (up to 80%) for apparently well-behaved sources with a companion within $\lesssim 4$ arcsec, and much larger for sources with poor astrometric fits. We provide an empirical fitting function to inflate published $σ_{\varpi}$ values for isolated sources. The public catalog offers wide ranging follow-up opportunities: from calibrating spectroscopic surveys, to precisely constraining ages of field stars, to the masses and the initial-final mass relation of white dwarfs, to dynamically probing the Galactic tidal field. △ Less

Submitted 2 February, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: Accepted to MNRAS. Catalog available at https://zenodo.org/record/4435257

TOI-1259Ab -- a gas giant planet with 2.7% deep transits and a bound white dwarf companion

Authors: David V. Martin, Kareem El-Badry, Vedad Kunovac Hodžić, Amaury H. M. J. Triaud, Ruth Angus, Jessica Birky, Daniel Foreman-Mackey, Christina Hedges, Benjamin T. Montet, Simon J. Murphy, Alexandre Santerne, Keivan G. Stassun, Alexander P. Stephan, Ji Wang, Paul Benni, Vadim Krushinsky, Nikita Chazov, Nikolay Mishevskiy, Carl Ziegler, Abderahmane Soubkiou, Zouhair Benkhaldoun, Douglas A. Caldwell, Karen Collins, Christopher E. Henze, Natalia M. Guerrero , et al. (11 additional authors not shown)

Abstract: We present TOI-1259Ab, a 1.0 Rjup gas giant planet transiting a 0.71 Rsun K-dwarf on a 3.48 day orbit. The system also contains a bound white dwarf companion TOI-1259B with a projected distance of approximately 1600 AU from the planet host. Transits are observed in nine TESS sector and are 2.7 per cent deep - among the deepest known - making TOI-1259Ab a promising target for atmospheric characteri… ▽ More We present TOI-1259Ab, a 1.0 Rjup gas giant planet transiting a 0.71 Rsun K-dwarf on a 3.48 day orbit. The system also contains a bound white dwarf companion TOI-1259B with a projected distance of approximately 1600 AU from the planet host. Transits are observed in nine TESS sector and are 2.7 per cent deep - among the deepest known - making TOI-1259Ab a promising target for atmospheric characterization. Our follow-up radial velocity measurements indicate a variability of semiamplitude K = 71 m/s, implying a planet mass of 0.44 Mjup. By fitting the spectral energy distribution of the white dwarf we derive a total age of 4.08 (+1.21 -0.53) Gyr for the system. The K-dwarf's light curve reveals a rotational variability with a period of 28 days, which implies a gyrochronology age broadly consistent with the white dwarf's total age. △ Less

Submitted 1 July, 2021; v1 submitted 7 January, 2021; originally announced January 2021.

Comments: Accepted to MNRAS. Some structural changes from first arXiv version but no significant changes to results. One figure got a bit prettier

Probing the CGM of Low-redshift Dwarf Galaxies Using FIRE Simulations

Authors: Fei Li, Mubdi Rahman, Norman Murray, Zachary Hafen, Claude-André Faucher-Giguère, Jonathan Stern, Cameron B. Hummels, Philip F. Hopkins, Kareem El-Badry, Dušan Kereš

Abstract: Observations of UV metal absorption lines have provided insight into the structure and composition of the circumgalactic medium (CGM) around galaxies. We compare these observations with the low-redshift ($z \leq 0.3$) CGM around dwarf galaxies in high-resolution cosmological zoom-in runs in the FIRE-2 simulation suite. We select simulated galaxies that match the halo mass, stellar mass, and redshi… ▽ More Observations of UV metal absorption lines have provided insight into the structure and composition of the circumgalactic medium (CGM) around galaxies. We compare these observations with the low-redshift ($z \leq 0.3$) CGM around dwarf galaxies in high-resolution cosmological zoom-in runs in the FIRE-2 simulation suite. We select simulated galaxies that match the halo mass, stellar mass, and redshift of the observed samples. We produce absorption measurements using Trident for UV transitions of C IV, O VI, Mg II and Si III. The FIRE equivalent width (EW) distributions and covering fractions for the C IV ion are broadly consistent with observations inside $0.5 R_{vir}$, but are under-predicted for O VI, Mg II, and Si III. The absorption strengths of the ions in the CGM are moderately correlated with the masses and star formation activity of the galaxies. The correlation strengths increase with the ionization potential of the ions. The structure and composition of the gas from the simulations exhibit three zones around dwarf galaxies characterized by distinct ion column densities: the disky ISM, the inner CGM (the wind-dominated regime), and the outer CGM (the IGM accretion-dominated regime). We find that the outer CGM in the simulations is nearly but not quite supported by thermal pressure, so it is not in hydrostatic equilibrium (HSE), resulting in halo-scale bulk inflow and outflow motions. The net gas inflow rates are comparable to the SFR of the galaxy, but the bulk inflow and outflow rates are greater by an order of magnitude, with velocities comparable to the virial velocity of the halo. These roughly virial velocities (${\sim} 100 km s^{-1}$) produce large EWs in the simulations. This supports a picture for dwarf galaxies in which the dynamics of the CGM at large scales are coupled to the small-scale star formation activity near the centre of their halos. △ Less

Submitted 3 November, 2020; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: To appear in MNRAS

A Relationship Between Stellar Metallicity Gradients and Galaxy Age in Dwarf Galaxies

Authors: Francisco J. Mercado, James S. Bullock, Michael Boylan-Kolchin, Jorge Moreno, Andrew Wetzel, Kareem El-Badry, Andrew S. Graus, Alex Fitts, Philip F. Hopkins, Claude-André Faucher-Giguère

Abstract: We explore the origin of stellar metallicity gradients in simulated and observed dwarf galaxies. We use FIRE-2 cosmological baryonic zoom-in simulations of 26 isolated galaxies as well as existing observational data for 10 Local Group dwarf galaxies. Our simulated galaxies have stellar masses between $10^{5.5}$ and $10^{8.6} \msun$. Whilst gas-phase metallicty gradients are generally weak in our s… ▽ More We explore the origin of stellar metallicity gradients in simulated and observed dwarf galaxies. We use FIRE-2 cosmological baryonic zoom-in simulations of 26 isolated galaxies as well as existing observational data for 10 Local Group dwarf galaxies. Our simulated galaxies have stellar masses between $10^{5.5}$ and $10^{8.6} \msun$. Whilst gas-phase metallicty gradients are generally weak in our simulated galaxies, we find that stellar metallicity gradients are common, with central regions tending to be more metal-rich than the outer parts. The strength of the gradient is correlated with galaxy-wide median stellar age, such that galaxies with younger stellar populations have flatter gradients. Stellar metallicty gradients are set by two competing processes: (1) the steady "puffing" of old, metal-poor stars by feedback-driven potential fluctuations, and (2) the accretion of extended, metal-rich gas at late times, which fuels late-time metal-rich star formation. If recent star formation dominates, then extended, metal-rich star formation washes out pre-existing gradients from the "puffing" process. We use published results from ten Local Group dwarf galaxies to show that a similar relationship between age and stellar metallicity-gradient strength exists among real dwarfs. This suggests that observed stellar metallicity gradients may be driven largely by the baryon/feedback cycle rather than by external environmental effects. △ Less

Submitted 2 December, 2020; v1 submitted 2 September, 2020; originally announced September 2020.

Comments: 14 pages, 13 figures, submitted to MNRAS

Data-driven spectroscopic estimates of absolute magnitude, distance and binarity -- method and catalog of 16,002 O- and B-type stars from LAMOST

Authors: Mao-Sheng Xiang, Hans-Walter Rix, Yuan-Sen Ting, Eleonora Zari, Kareem El-Badry, Hai-Bo Yuan, Wen-Yuan Cui

Abstract: We present a data-driven method to estimate absolute magnitudes for O- and B-type stars from the LAMOST spectra, which we combine with {\it Gaia} parallaxes to infer distance and binarity. The method applies a neural network model trained on stars with precise {\it Gaia} parallax to the spectra and predicts $K_{\rm s}$-band absolute magnitudes $M_{Ks}$ with a precision of 0.25\,mag, which correspo… ▽ More We present a data-driven method to estimate absolute magnitudes for O- and B-type stars from the LAMOST spectra, which we combine with {\it Gaia} parallaxes to infer distance and binarity. The method applies a neural network model trained on stars with precise {\it Gaia} parallax to the spectra and predicts $K_{\rm s}$-band absolute magnitudes $M_{Ks}$ with a precision of 0.25\,mag, which corresponds to a precision of 12\% in spectroscopic distance. For distant stars (e.g. $>5$\,kpc), the inclusion of constraints from spectroscopic $M_{Ks}$ significantly improves the distance estimates compared to inferences from {\it Gaia} parallax alone. Our method accommodates for emission line stars by first identifying them via PCA reconstructions and then treating them separately for the $M_{Ks}$ estimation. We also take into account unresolved binary/multiple stars, which we identify through deviations in the spectroscopic $M_{Ks}$ from the geometric $M_{Ks}$ inferred from {\it Gaia} parallax. This method of binary identification is particularly efficient for unresolved binaries with near equal-mass components and thus provides an useful supplementary way to identify unresolved binary or multiple-star systems. We present a catalog of spectroscopic $M_{Ks}$, extinction, distance, flags for emission lines, and binary classification for 16,002 OB stars from LAMOST DR5. As an illustration of the method, we determine the $M_{Ks}$ and distance to the enigmatic LB-1 system, where Liu et al. (2019) had argued for the presence of a black hole and incorrect parallax measurement, and we do not find evidence for errorneous {\it Gaia} parallax. △ Less

Submitted 24 August, 2020; originally announced August 2020.

Comments: 21 pages, 16 figures, submitted to ApJ Supplement. Comments are welcome

The time-scales probed by star formation rate indicators for realistic, bursty star formation histories from the FIRE simulations

Authors: José A. Flores Velázquez, Alexander B. Gurvich, Claude-André Faucher-Giguère, James S. Bullock, Tjitske K. Starkenburg, Jorge Moreno, Alexandres Lazar, Francisco J. Mercado, Jonathan Stern, Martin Sparre, Christopher C. Hayward, Andrew Wetzel, Kareem El-Badry

Abstract: Understanding the rate at which stars form is central to studies of galaxy formation. Observationally, the star formation rates (SFRs) of galaxies are measured using the luminosity in different frequency bands, often under the assumption of a time-steady SFR in the recent past. We use star formation histories (SFHs) extracted from cosmological simulations of star-forming galaxies from the FIRE pro… ▽ More Understanding the rate at which stars form is central to studies of galaxy formation. Observationally, the star formation rates (SFRs) of galaxies are measured using the luminosity in different frequency bands, often under the assumption of a time-steady SFR in the recent past. We use star formation histories (SFHs) extracted from cosmological simulations of star-forming galaxies from the FIRE project to analyze the time-scales to which the H$α$ and far-ultraviolet (FUV) continuum SFR indicators are sensitive. In these simulations, the SFRs are highly time variable for all galaxies at high redshift, and continue to be bursty to z=0 in dwarf galaxies. When FIRE SFHs are partitioned into their bursty and time-steady phases, the best-fitting FUV time-scale fluctuates from its ~10 Myr value when the SFR is time-steady to >~100 Myr immediately following particularly extreme bursts of star formation during the bursty phase. On the other hand, the best-fitting averaging time-scale for H$α$ is generally insensitive to the SFR variability in the FIRE simulations and remains ~5 Myr at all times. These time-scales are shorter than the 100 Myr and 10 Myr time-scales sometimes assumed in the literature for FUV and H$α$, respectively, because while the FUV emission persists for stellar populations older than 100 Myr, the time-dependent luminosities are strongly dominated by younger stars. Our results confirm that the ratio of SFRs inferred using H$α$ vs. FUV can be used to probe the burstiness of star formation in galaxies. △ Less

Submitted 15 December, 2020; v1 submitted 19 August, 2020; originally announced August 2020.

Comments: 14 pages, 10 figures, accepted to MNRAS

Spatial and Kinematic Clustering of Stars in the Galactic Disk

Authors: Harshil Kamdar, Charlie Conroy, Yuan-Sen Ting, Kareem El-Badry

Abstract: The Galactic disk is expected to be spatially, kinematically, and chemically clustered on many scales due to both star formation and non-axisymmetries in the Galactic potential. In this work we calculate the spatial and kinematic two-point correlation functions using a sample of $1.7 \times 10^6$ stars within 1 kpc of the Sun with 6D phase space information available from \textit{Gaia} DR2. Cluste… ▽ More The Galactic disk is expected to be spatially, kinematically, and chemically clustered on many scales due to both star formation and non-axisymmetries in the Galactic potential. In this work we calculate the spatial and kinematic two-point correlation functions using a sample of $1.7 \times 10^6$ stars within 1 kpc of the Sun with 6D phase space information available from \textit{Gaia} DR2. Clustering is detected on spatial scales of 1-300 pc and velocity scales of at least 15 km s$^{-1}$. With bound structures included, the data have a power-law index ($ξ(Δr) \propto Δr^γ$) of $γ\approx-2$ at most spatial scales, which is in line with theoretical predictions. After removing bound structures, the data have a power-law index of $γ\approx-1$ for $<100$ pc and $γ\lesssim -1.5$ for $>100$ pc. We interpret these results with the aid of a novel star-by-star simulation of the Galaxy in which stars are born in clusters orbiting in a realistic potential that includes spiral arms, a bar, and GMCs. We find that the simulation largely agrees with the observations (within a factor of 2-3) at all spatial and kinematic scales. In detail, the correlation function in the simulation is shallower than the data at $< 20$ pc scales, and steeper than the data at $> 30$ pc scales. We also find a persistent clustering signal in the kinematic correlation function for the data at large $Δv$ ($>5$ km s$^{-1}$) not present in the simulations. We speculate that this mismatch between observations and simulations may be due to two processes not included in the simulation: hierarchical star formation and transient spiral arms. We also use the simulations to predict the clustering signal as a function of pair-wise metallicity and age separations. Ages and metallicities measured with a precision of $50\%$ and $0.05$ dex are required in order to enhance the clustering signal. △ Less

Submitted 21 July, 2020; originally announced July 2020.

Comments: 21 pages, 13 figures. Submitted to ApJ

A proto-helium white dwarf stripped by a substellar companion via common-envelope ejection -- Uncovering the true nature of a candidate hypervelocity B-type star

Authors: Andreas Irrgang, Stephan Geier, Ulrich Heber, Thomas Kupfer, Kareem El-Badry, Steven Bloemen

Abstract: In the past, SDSS J160429.12+100002.2 was spectroscopically classified as a blue horizontal branch (BHB) star. Assuming a luminosity that is characteristic of BHB stars, the object's radial velocity and proper motions from Gaia Early Data Release 3 would imply that its Galactic rest-frame velocity exceeds its local escape velocity. Consequently, the object would be considered a hypervelocity star,… ▽ More In the past, SDSS J160429.12+100002.2 was spectroscopically classified as a blue horizontal branch (BHB) star. Assuming a luminosity that is characteristic of BHB stars, the object's radial velocity and proper motions from Gaia Early Data Release 3 would imply that its Galactic rest-frame velocity exceeds its local escape velocity. Consequently, the object would be considered a hypervelocity star, which would prove particularly interesting because its Galactic trajectory points in our direction. However, based on the spectroscopic analysis of follow-up observations, we show that the object is actually a short-period ($P \approx 3.4$ h) single-lined spectroscopic binary system with a visible B-type star (effective temperature $T_{\mathrm{eff}} = 15\,840\pm160$ K and surface gravity $\log(g) = 4.86\pm0.04$) that is less luminous than typical BHB stars. Accordingly, the distance of the system is lower than originally thought, which renders its Galactic orbit bound to the Galaxy. Nevertheless, it is still an extreme halo object on a highly retrograde orbit. The abundances of He, C, N, O, Ne, Mg, Al, Si, S, and Ca are subsolar by factors from 3 to more than 100, while Fe is enriched by a factor of about 6. This peculiar chemical composition pattern is most likely caused by atomic diffusion processes. Combining constraints from astrometry, orbital motion, photometry, and spectroscopy, we conclude that the visible component is an unevolved proto-helium white dwarf with a thin hydrogen envelope that was stripped by a substellar companion through common-envelope ejection. Its unique configuration renders the binary system an interesting test bed for stellar binary evolution in general and common-envelope evolution in particular. △ Less

Submitted 19 June, 2021; v1 submitted 7 July, 2020; originally announced July 2020.

Comments: Accepted for publication in A&A (Astronomy and Astrophysics)

Journal ref: A&A 650, A102 (2021)

Virialization of the inner CGM in the FIRE simulations and implications for galaxy discs, star formation and feedback

Authors: Jonathan Stern, Claude-André Faucher-Giguère, Drummond Fielding, Eliot Quataert, Zachary Hafen, Alexander B. Gurvich, Xiangcheng Ma, Lindsey Byrne, Kareem El-Badry, Daniel Anglés-Alcázar, T. K. Chan, Robert Feldmann, Dušan Kereš, Andrew Wetzel, Norman Murray, Philip F. Hopkins

Abstract: We use the FIRE-2 cosmological simulations to study the formation of a quasi-static, virial-temperature gas phase in the circumgalactic medium (CGM) at redshifts 0<z<5, and how the formation of this virialized phase affects the evolution of galactic discs. We demonstrate that when the halo mass crosses ~10^12 M_sun, the cooling time of shocked gas in the inner CGM (~0.1 R_vir, where R_vir is the v… ▽ More We use the FIRE-2 cosmological simulations to study the formation of a quasi-static, virial-temperature gas phase in the circumgalactic medium (CGM) at redshifts 0<z<5, and how the formation of this virialized phase affects the evolution of galactic discs. We demonstrate that when the halo mass crosses ~10^12 M_sun, the cooling time of shocked gas in the inner CGM (~0.1 R_vir, where R_vir is the virial radius) exceeds the local free-fall time. The inner CGM then experiences a transition from on average sub-virial temperatures (T<<T_vir), large pressure fluctuations and supersonic inflow/outflow velocities, to virial temperatures (T~T_vir), uniform pressures and subsonic velocities. This transition occurs when the outer CGM (~0.5 R_vir) is already subsonic and has a temperature ~T_vir, indicating that the longer cooling times at large radii allow the outer CGM to virialize at lower halo masses than the inner CGM. This outside-in CGM virialization scenario is in contrast with inside-out scenarios commonly envisioned based on more idealized simulations. We demonstrate that inner CGM virialization coincides with abrupt changes in the central galaxy and its stellar feedback: the galaxy settles into a stable rotating disc, star formation transitions from `bursty' to `steady,' and stellar-driven galaxy-scale outflows are suppressed. Our results thus suggest that CGM virialization is initially associated with the formation of rotation-dominated thin galactic discs, rather than with the quenching of star formation as often assumed. △ Less

Submitted 29 December, 2020; v1 submitted 24 June, 2020; originally announced June 2020.

Comments: Accepted for publication in ApJ. Online versions of figures 13 and 15 are available at https://sites.northwestern.edu/jonathanstern/online-figures/

A stripped-companion origin for Be stars: clues from the putative black holes HR 6819 and LB-1

Authors: Kareem El-Badry, Eliot Quataert

Abstract: HR 6819 is a bright ($V=5.36$), blue star recently proposed to be a triple containing a detached black hole (BH). We show that the system is a binary and does not contain a BH. Using spectral decomposition, we disentangle the observed composite spectra into two components: a rapidly rotating Be star and a slowly rotating B star with low surface gravity $(\log g \approx 2.75)$. Both stars show peri… ▽ More HR 6819 is a bright ($V=5.36$), blue star recently proposed to be a triple containing a detached black hole (BH). We show that the system is a binary and does not contain a BH. Using spectral decomposition, we disentangle the observed composite spectra into two components: a rapidly rotating Be star and a slowly rotating B star with low surface gravity $(\log g \approx 2.75)$. Both stars show periodic radial velocity (RV) variability, but the RV semi-amplitude of the B star's orbit is $K_{\rm B}= (62.7 \pm 1)\,\rm km\,s^{-1}$, while that of the Be star is only $K_{\rm Be} = (4.5\pm 2)\,\rm km\,s^{-1}$. This implies that the B star is less massive by at least a factor of 10. The surface abundances of the B star bear imprints of CNO burning. We argue that the B star is a bloated, recently stripped helium star with mass $\approx 0.5\,M_{\odot}$ that is currently contracting to become a hot subdwarf. The orbital motion of the Be star obviates the need for a BH to explain the B star's motion. A stripped-star model reproduces the observed luminosity of the system, while a normal star with the B star's temperature and gravity would be more than 10 times too luminous. HR 6819 and the binary LB-1 probably formed through similar channels. We use MESA models to investigate their evolutionary history, finding that they likely formed from intermediate-mass ($3-7\,M_{\odot}$) primaries stripped by slightly lower-mass secondaries and are progenitors to Be + sdOB binaries such as $φ$ Persei. The lifetime of their current evolutionary phase is on average $2\times 10^5$ years, of order half a percent of the total lifetime of the Be phase. This implies that many Be stars have hot subdwarf and white dwarf companions, and that a substantial fraction ($20-100\%$) of field Be stars form through accretion of material from a binary companion. △ Less

Submitted 29 January, 2021; v1 submitted 21 June, 2020; originally announced June 2020.

Comments: Accepted to MNRAS. Figure E1 is new

A dark matter profile to model diverse feedback-induced core sizes of $Λ$CDM haloes

Authors: Alexandres Lazar, James S. Bullock, Michael Boylan-Kolchin, T. K. Chan, Philip F. Hopkins, Andrew S. Graus, Andrew Wetzel, Kareem El-Badry, Coral Wheeler, Maria C. Straight, Dušan Kereš, Claude-André Faucher-Giguère, Alex Fitts, Shea Garrison-Kimmel

Abstract: We analyze the cold dark matter density profiles of 54 galaxy halos simulated with FIRE-2 galaxy formation physics, each resolved within $0.5\%$ of the halo virial radius. These halos contain galaxies with masses that range from ultra-faint dwarfs ($M_\star \simeq 10^{4.5} M_{\odot}$) to the largest spirals ($M_\star \simeq 10^{11} M_{\odot}$) and have density profiles that are both cored and cusp… ▽ More We analyze the cold dark matter density profiles of 54 galaxy halos simulated with FIRE-2 galaxy formation physics, each resolved within $0.5\%$ of the halo virial radius. These halos contain galaxies with masses that range from ultra-faint dwarfs ($M_\star \simeq 10^{4.5} M_{\odot}$) to the largest spirals ($M_\star \simeq 10^{11} M_{\odot}$) and have density profiles that are both cored and cuspy. We characterize our results using a new analytic density profile that extends the standard Einasto form to allow for a pronounced constant-density core in the resolved innermost radius. With one additional core-radius parameter, $r_{c}$, this "core-Einasto" profile is able to characterize the shape and normalization of our feedback-impacted dark matter halos. In order to enable comparisons with observations, we provide fitting functions for $r_{c}$ and other profile parameters as a function of both $M_\star$ and $M_{\star}/M_{\rm halo}$. In agreement with similar studies done in the literature, we find that dark matter core formation is most efficient at the characteristic stellar-mass to halo-mass ratio $M_\star/M_{\rm halo} \simeq 5 \times 10^{-3}$, or $M_{\star} \sim 10^9 \, M_{\odot}$, with cores that are roughly the size of the galaxy half-light radius, $r_{c} \simeq 1-5$ kpc. Furthermore, we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with $M_{\star}/M_{\rm halo} < 5\times 10^{-4}$ or $M_\star \lesssim 10^6 \, M_{\odot}$. For Milky Way-size galaxies, baryonic contraction often makes halos significantly more concentrated and dense at the stellar half-light radius than dark matter only runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of $\simeq 0.5-2$ kpc in size. Recent evidence for a ${\sim} 2$ kpc core in the Milky Way's dark matter halo is consistent with this expectation. △ Less

Submitted 8 July, 2020; v1 submitted 22 April, 2020; originally announced April 2020.

Comments: 27 pages; 19 figures; Accepted by MNRAS

From birth associations to field stars: mapping the small-scale orbit distribution in the Galactic disc

Authors: Johanna Coronado, Hans-Walter Rix, Wilma H. Trick, Kareem El-Badry, Jan Rybizki, Maosheng Xiang

Abstract: Stars born at the same time in the same place should have formed from gas of the same element composition. But most stars subsequently disperse from their birth siblings, in orbit and orbital phase, becoming 'field stars'. Here we explore and provide direct observational evidence for this process in the Milky Way disc, by quantifying the probability that orbit-similarity among stars implies indist… ▽ More Stars born at the same time in the same place should have formed from gas of the same element composition. But most stars subsequently disperse from their birth siblings, in orbit and orbital phase, becoming 'field stars'. Here we explore and provide direct observational evidence for this process in the Milky Way disc, by quantifying the probability that orbit-similarity among stars implies indistinguishable metallicity. We define the orbit similarity among stars through their distance in action-angle space, $Δ(J,θ)$, and their abundance similarity simply by $Δ$[Fe/H]. Analyzing a sample of main sequence stars from Gaia DR2 and LAMOST, we find an excess of pairs with the same metallicity ($Δ\mathrm{[Fe/H]}<0.1$) that extends to remarkably large separations in $Δ(J,θ)$ that correspond to nearly 1 kpc distances. We assess the significance of this effect through a mock sample, drawn from a smooth and phase-mixed orbit distribution. Through grouping such star pairs into associations with a friend-of-friends algorithm linked by $Δ(J,θ)$, we find 100s of mono-abundance groups with $\ge 3$ (to $\gtrsim 20$) members; these groups -- some clusters, some spread across the sky -- are over an order-of-magnitude more abundant than expected for a smooth phase-space distribution, suggesting that we are witnessing the 'dissolution' of stellar birth associations into the field. △ Less

Submitted 18 May, 2020; v1 submitted 21 February, 2020; originally announced February 2020.

Comments: 16 pages, 13 figures. Accepted for publication in MNRAS

The Formation Times and Building Blocks of Milky Way-mass Galaxies in the FIRE Simulations

Authors: Isaiah B. Santistevan, Andrew Wetzel, Kareem El-Badry, Joss Bland-Hawthorn, Michael Boylan-Kolchin, Jeremy Bailin, Claude-Andre Faucher-Giguere, Samantha Benincasa

Abstract: Surveys of the Milky Way (MW) and M31 enable detailed studies of stellar populations across ages and metallicities, with the goal of reconstructing formation histories across cosmic time. These surveys motivate key questions for galactic archaeology in a cosmological context: when did the main progenitor of a MW/M31-mass galaxy form, and what were the galactic building blocks that formed it? We in… ▽ More Surveys of the Milky Way (MW) and M31 enable detailed studies of stellar populations across ages and metallicities, with the goal of reconstructing formation histories across cosmic time. These surveys motivate key questions for galactic archaeology in a cosmological context: when did the main progenitor of a MW/M31-mass galaxy form, and what were the galactic building blocks that formed it? We investigate the formation times and progenitor galaxies of MW/M31-mass galaxies using the FIRE-2 cosmological simulations, including 6 isolated MW/M31-mass galaxies and 6 galaxies in Local Group (LG)-like pairs at z = 0. We examine main progenitor "formation" based on two metrics: (1) transition from primarily ex-situ to in-situ stellar mass growth and (2) mass dominance compared to other progenitors. We find that the main progenitor of a MW/M31-mass galaxy emerged typically at z ~ 3-4 (11.6-12.2 Gyr ago), while stars in the bulge region (inner 2 kpc) at z = 0 formed primarily in a single main progenitor at z < 5 (< 12.6 Gyr ago). Compared with isolated hosts, the main progenitors of LG-like paired hosts emerged significantly earlier (Δz ~ 2, Δt ~ 1.6 Gyr), with ~ 4x higher stellar mass at all z > 4 (> 12.2 Gyr ago). This highlights the importance of environment in MW/M31-mass galaxy formation, especially at early times. Overall, about 100 galaxies with M_star > 10^5 M_sun formed a typical MW/M31-mass system. Thus, surviving satellites represent a highly incomplete census (by ~ 5x) of the progenitor population. △ Less

Submitted 28 July, 2020; v1 submitted 9 January, 2020; originally announced January 2020.

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

Effects of Stellar Feedback on Stellar and Gas Kinematics of Star-Forming Galaxies at 0.6<z<1.0

Authors: Debora Pelliccia, Bahram Mobasher, Behnam Darvish, Brian C. Lemaux, Lori M. Lubin, Jessie Hirtenstein, Lu Shen, Po-Feng Wu, Kareem El-Badry, Andrew Wetzel, Tucker Jones

Abstract: Recent zoom-in cosmological simulations have shown that stellar feedback can flatten the inner density profile of the dark matter halo in low-mass galaxies. A correlation between the stellar/gas velocity dispersion ($σ_{star}$, $σ_{gas}$) and the specific star formation rate (sSFR) is predicted as an observational test of the role of stellar feedback in re-shaping the dark matter density profile.… ▽ More Recent zoom-in cosmological simulations have shown that stellar feedback can flatten the inner density profile of the dark matter halo in low-mass galaxies. A correlation between the stellar/gas velocity dispersion ($σ_{star}$, $σ_{gas}$) and the specific star formation rate (sSFR) is predicted as an observational test of the role of stellar feedback in re-shaping the dark matter density profile. In this work we test the validity of this prediction by studying a sample of star-forming galaxies at $0.6<z<1.0$ from the LEGA-C survey, which provides high signal-to-noise measurements of stellar and gas kinematics. We find that a weak but significant correlation between $σ_{star}$ (and $σ_{gas}$) and sSFR indeed exists for galaxies in the lowest mass bin (M$_\ast\sim10^{10}\,$M$_\odot$). This correlation, albeit with a $\sim$35% scatter, holds for different tracers of star formation, and becomes stronger with redshift. This result generally agrees with the picture that at higher redshifts star formation rate was generally higher, and galaxies at M$_\ast\lesssim10^{10}\,$M$_\odot$ have not yet settled into a disk. As a consequence, they have shallower gravitational potentials more easily perturbed by stellar feedback. The observed correlation between $σ_{star}$ (and $σ_{gas}$) and sSFR supports the scenario predicted by cosmological simulations, in which feedback-driven outflows cause fluctuations in the gravitation potential that flatten the density profiles of low-mass galaxies. △ Less

Submitted 12 June, 2020; v1 submitted 2 January, 2020; originally announced January 2020.

Comments: 8 pages, 2 figures, 3 tables. Revised version. ApJL in press

Identical or fraternal twins? : The chemical homogeneity of wide binaries from Gaia DR2

Authors: Keith Hawkins, Madeline Lucey, Yuan-Sen Ting, Alexander Ji, Dustin Katzberg, Megan Thompson, Kareem El-Badry, Johanna Teske, Tyler Nelson, Andreia Carrillo

Abstract: One of the high-level goals of Galactic archaeology is chemical tagging of stars across the Milky Way to piece together its assembly history. For this to work, stars born together must be uniquely chemically homogeneous. Wide binary systems are an important laboratory to test this underlying assumption. Here we present the detailed chemical abundance patterns of 50 stars across 25 wide binary syst… ▽ More One of the high-level goals of Galactic archaeology is chemical tagging of stars across the Milky Way to piece together its assembly history. For this to work, stars born together must be uniquely chemically homogeneous. Wide binary systems are an important laboratory to test this underlying assumption. Here we present the detailed chemical abundance patterns of 50 stars across 25 wide binary systems comprised of main-sequence stars of similar spectral type identified in Gaia DR2 with the aim of quantifying their level of chemical homogeneity. Using high-resolution spectra obtained with McDonald Observatory, we derive stellar atmospheric parameters and precise detailed chemical abundances for light/odd-Z (Li, C, Na, Al, Sc, V, Cu), $α$ (Mg, Si, Ca), Fe-peak (Ti, Cr, Mn, Fe, Co, Ni, Zn), and neutron capture (Sr, Y, Zr, Ba, La, Nd, Eu) elements. Results indicate that 80% (20 pairs) of the systems are homogeneous in [Fe/H] at levels below 0.02 dex. These systems are also chemically homogeneous in all elemental abundances studied, with offsets and dispersions consistent with measurement uncertainties. We also find that wide binary systems are far more chemically homogeneous than random pairings of field stars of similar spectral type. These results indicate that wide binary systems tend to be chemically homogeneous but in some cases they can differ in their detailed elemental abundances at a level of [X/H] ~ 0.10 dex, overall implying chemical tagging in broad strokes can work. △ Less

Submitted 18 December, 2019; originally announced December 2019.

Comments: 20 pages, 9 Figures, 5 Tables, Accepted for Publication in MNRAS

Not so fast: LB-1 is unlikely to contain a 70 $M_{\odot}$ black hole

Authors: Kareem El-Badry, Eliot Quataert

Abstract: The recently discovered binary LB-1 has been reported to contain a $\sim$\,$70\,M_{\odot}$ black hole (BH). The evidence for the unprecedentedly high mass of the unseen companion comes from reported radial velocity (RV) variability of the H$α$ emission line, which has been proposed to originate from an accretion disk around a BH. We show that there is in fact no evidence for RV variability of the… ▽ More The recently discovered binary LB-1 has been reported to contain a $\sim$\,$70\,M_{\odot}$ black hole (BH). The evidence for the unprecedentedly high mass of the unseen companion comes from reported radial velocity (RV) variability of the H$α$ emission line, which has been proposed to originate from an accretion disk around a BH. We show that there is in fact no evidence for RV variability of the H$α$ emission line, and that its apparent shifts instead originate from shifts in the luminous star's H$α$ absorption line. If not accounted for, such shifts will cause a stationary emission line to appear to shift in anti-phase with the luminous star. We show that once the template spectrum of a B star is subtracted from the observed Keck/HIRES spectra of LB-1, evidence for RV variability vanishes. Indeed, the data rule out periodic variability of the line with velocity semi-amplitude $K_{\rm Hα} > 1.3\,\rm km\,s^{-1}$. This strongly suggests that the observed H$α$ emission does not originate primarily from an accretion disk around a BH, and thus that the mass ratio cannot be constrained from the relative velocity amplitudes of the emission and absorption lines. The nature of the unseen companion remains uncertain, but a "normal" stellar-mass BH with mass $5\lesssim M/M_{\odot}\lesssim 20 $ seems most plausible. The H$α$ emission likely originates primarily from circumbinary material, not from either component of the binary. △ Less

Submitted 3 January, 2020; v1 submitted 9 December, 2019; originally announced December 2019.

Comments: 5 pages, 4 figures. Accepted to MNRAS letters

The Fates of the Circumgalactic Medium in the FIRE Simulations

Authors: Z. Hafen, C. -A. Faucher-Giguere, D. Angles-Alcazar, J. Stern, D. Keres, C. Esmerian, A. Wetzel, K. El-Badry, T. K. Chan, N. Murray

Abstract: We analyze the different fates of the circumgalactic medium (CGM) in FIRE-2 cosmological simulations, focusing on the redshifts z=0.25 and z=2 representative of recent surveys. Our analysis includes 21 zoom-in simulations covering the halo mass range Mh(z=0) ~ 10^10 - 10^12 Msun. We analyze both where the gas ends up after first leaving the CGM (its "proximate" fate), as well as its location at z=… ▽ More We analyze the different fates of the circumgalactic medium (CGM) in FIRE-2 cosmological simulations, focusing on the redshifts z=0.25 and z=2 representative of recent surveys. Our analysis includes 21 zoom-in simulations covering the halo mass range Mh(z=0) ~ 10^10 - 10^12 Msun. We analyze both where the gas ends up after first leaving the CGM (its "proximate" fate), as well as its location at z=0 (its "ultimate" fate). Of the CGM at z=2, about half is found in the ISM or stars of the central galaxy by z=0 in Mh(z=2) ~ 5e11 Msun halos, but most of the CGM in lower-mass halos ends up in the IGM. This is so even though most of the CGM in M_h(z=2) ~ 5e10 Msun halos first accretes onto the central galaxy before being ejected into the IGM. On the other hand, most of the CGM mass at z=0.25 remains in the CGM by z=0 at all halo masses analyzed. Of the CGM gas that subsequently accretes onto the central galaxy in the progenitors of Mh(z=0) ~10^12 Msun halos, most of it is cool (T~10^4 K) at z=2 but hot (~Tvir) at z=0.25, consistent with the expected transition from cold mode to hot mode accretion. Despite the transition in accretion mode, at both z=0.25 and z=2 >~80% of the cool gas in Mh >~ 10^11 Msun halos will accrete onto a galaxy. We find that the metallicity of CGM gas is typically a poor predictor of both its proximate and ultimate fates. This is because there is in general little correlation between the origin of CGM gas and its fate owing to substantial mixing while in the CGM. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Comments: 16 pages, 10 figures

The Separation Distribution of Ultra-Wide Binaries across Galactic Populations

Authors: Hai-Jun Tian, Kareem El-Badry, Hans-Walter Rix, Andrew Gould

Abstract: We present an extensive sample of ultra-wide binary stars in the solar neighborhood, focusing on separations of $0.01<s/{pc}<1$. Using data from Gaia DR2, we define kinematic sub-populations via the systems' tangential velocities, i.e., disk-like ($v_{\perp, tot}<40$km/s), intermediate ($v_{\perp,tot}=(40-85)$km/s), and halo-like ($v_{\perp,tot}>85$km/s) samples, presuming that these velocity cuts… ▽ More We present an extensive sample of ultra-wide binary stars in the solar neighborhood, focusing on separations of $0.01<s/{pc}<1$. Using data from Gaia DR2, we define kinematic sub-populations via the systems' tangential velocities, i.e., disk-like ($v_{\perp, tot}<40$km/s), intermediate ($v_{\perp,tot}=(40-85)$km/s), and halo-like ($v_{\perp,tot}>85$km/s) samples, presuming that these velocity cuts represent a rough ordering in the binaries' age and metallicity. Through stringent cuts on astrometric precision, we can obtain pure binary samples with thousands of binaries in each sample. For all three populations, the distribution of binary separations extends smoothly to 1pc, displaying neither strong truncation nor bimodality. Fitting a smoothly-broken power law for the separation distribution, we find that its slope at separations $s=10^{2.5-4}$AU is the same for all sub-populations, $p(s)\propto s^γ$, with $γ\sim-1.54$. However, the logarithmic slope of $p(s)$ steepens at $s \gtrsim 10^4$AU. We find some evidences that the degree of steepening increases with the binaries' age, with a slope-change of only $Δγ\sim0.5$ for disk-like stars, but $Δγ>1$ for halo-like stars. This trend is contrary to what might be expected if steepening at wide separations were due to gravitational perturbations by molecular clouds or stars, which would preferentially disrupt disk binaries. If we were to interpret steepening at $s\gtrsim 10^4$AU as a consequence of disruption by massive halo objects, we would have to invoke a MACHO population inconsistent with other constraints. As a more plausible alternative, we propose a simple model to predict the separation distribution of wide binaries formed in dissolving star clusters. This model generically predicts $γ\sim-1.5$ as observed, with steepening at larger separations due to the finite size of binaries' birth clusters. △ Less

Submitted 6 November, 2019; v1 submitted 10 September, 2019; originally announced September 2019.

Comments: 16 pages, 9 figures, and 2 tables, Accepted to ApJS

Dark and luminous satellites of LMC-mass galaxies in the FIRE simulations

Authors: Ethan D. Jahn, Laura V. Sales, Andrew Wetzel, Michael Boylan-Kolchin, T. K. Chan, Kareem El-Badry, Alexandres Lazar, James S. Bullock

Abstract: Within $Λ$CDM, dwarf galaxies like the Large Magellanic Cloud (LMC) are expected to host numerous dark matter subhalos, several of which should host faint dwarf companions. Recent Gaia proper motions confirm new members of the LMC-system in addition to the previously known SMC, including two classical dwarf galaxies ($M_\ast$ > $10^5$ M$_{\odot}$; Carina and Fornax) as well as several ultra-faint… ▽ More Within $Λ$CDM, dwarf galaxies like the Large Magellanic Cloud (LMC) are expected to host numerous dark matter subhalos, several of which should host faint dwarf companions. Recent Gaia proper motions confirm new members of the LMC-system in addition to the previously known SMC, including two classical dwarf galaxies ($M_\ast$ > $10^5$ M$_{\odot}$; Carina and Fornax) as well as several ultra-faint dwarfs (Car2, Car3, Hor1, and Hyd1). We use the Feedback In Realistic Environments (FIRE) simulations to study the dark and luminous (down to ultrafaint masses, $M_\ast$ ~ $6\times10^3$ M$_{\odot}$) substructure population of isolated LMC-mass hosts ($M_\text{200}$ = $1-3\times10^{11}$ M$_{\odot}$) and place the Gaia + DES results in a cosmological context. By comparing number counts of subhalos in simulations with and without baryons, we find that, within 0.2 $r_\text{200}$, LMC-mass hosts deplete ~30% of their substructure, significantly lower than the ~70% of substructure depleted by Milky Way (MW) mass hosts. For our highest resolution runs ($m_\text{bary}$ = 880 M$_{\odot}$), ~5-10 subhalos form galaxies with $M_\ast$ > $10^4$ M$_{\odot}$, in agreement with the 7 observationally inferred pre-infall LMC companions. However, we find steeper simulated luminosity functions than observed, hinting at observation incompleteness at the faint end. The predicted DM content for classical satellites in FIRE agrees with observed estimates for Carina and Fornax, supporting the case for an LMC association. We predict that tidal stripping within the LMC potential lowers the inner dark matter density of ultra faint companions of the LMC. Thus, in addition to their orbital consistency, the low densities of dwarfs Car2, Hyd1, and Hyd2 reinforce their likelihood of Magellanic association. △ Less

Submitted 31 August, 2019; v1 submitted 5 July, 2019; originally announced July 2019.

Comments: 19 pages, 9 figures, 4 tables. Accepted for publication in MNRAS

Discovery of an equal-mass "twin" binary population reaching 1000+ AU separations

Authors: Kareem El-Badry, Hans-Walter Rix, Haijun Tian, Gaspard Duchêne, Maxwell Moe

Abstract: We use a homogeneous catalog of 42,000 main-sequence wide binaries identified by Gaia to measure the mass ratio distribution, p(q), of binaries with primary masses $0.1<M_1/M_{\odot}<2.5$, mass ratios $0.1 \lesssim q<1$, and separations $50<s/{\rm AU}<50,000$. A well-understood selection function allows us to constrain p(q) in 35 independent bins of primary mass and separation, with hundreds to th… ▽ More We use a homogeneous catalog of 42,000 main-sequence wide binaries identified by Gaia to measure the mass ratio distribution, p(q), of binaries with primary masses $0.1<M_1/M_{\odot}<2.5$, mass ratios $0.1 \lesssim q<1$, and separations $50<s/{\rm AU}<50,000$. A well-understood selection function allows us to constrain p(q) in 35 independent bins of primary mass and separation, with hundreds to thousands of binaries in each bin. Our investigation reveals a sharp excess of equal-mass "twin" binaries that is statistically significant out to separations of 1,000 to 10,000 AU, depending on primary mass. The excess is narrow: a steep increase in p(q) at $0.95 \lesssim q<1$, with no significant excess at $q\lesssim 0.95$. A range of tests confirm the signal is real, not a data artifact or selection effect. Combining the Gaia constraints with those from close binaries, we show that the twin excess decreases with increasing separation, but its width ($q\gtrsim 0.95$) is constant over $0.01<a/{\rm AU}<10,000$. The wide twin population would be difficult to explain if the components of all wide binaries formed via core fragmentation, which is not expected to produce strongly correlated component masses. We conjecture that wide twins formed at closer separations ($a \lesssim 100$ AU), likely via accretion from circumbinary disks, and were subsequently widened by dynamical interactions in their birth environments. The separation-dependence of the twin excess then constrains the efficiency of dynamical widening and disruption of binaries in young clusters. We also constrain p(q) across $0.1 \lesssim q<1$. Besides changes in the twin fraction, p(q) is independent of separation at fixed primary mass over $100 \lesssim s/{\rm AU} < 50,000$. It is flatter than expected for random pairings from the IMF but more bottom-heavy for wide binaries than for binaries with $a\lesssim$100 AU. △ Less

Submitted 3 September, 2019; v1 submitted 24 June, 2019; originally announced June 2019.

Comments: 24 pages, 17 figures, plus appendices. Accepted to MNRAS. Fig 17 and Appendix F are new since v1

A profile in FIRE: resolving the radial distributions of satellite galaxies in the Local Group with simulations

Authors: Jenna Samuel, Andrew Wetzel, Erik Tollerud, Shea Garrison-Kimmel, Sarah Loebman, Kareem El-Badry, Philip F. Hopkins, Michael Boylan-Kolchin, Claude-Andre Faucher-Giguere, James S. Bullock, Samantha Benincasa, Jeremy Bailin

Abstract: While many tensions between Local Group (LG) satellite galaxies and LCDM cosmology have been alleviated through recent cosmological simulations, the spatial distribution of satellites remains an important test of physical models and physical versus numerical disruption in simulations. Using the FIRE-2 cosmological zoom-in baryonic simulations, we examine the radial distributions of satellites with… ▽ More While many tensions between Local Group (LG) satellite galaxies and LCDM cosmology have been alleviated through recent cosmological simulations, the spatial distribution of satellites remains an important test of physical models and physical versus numerical disruption in simulations. Using the FIRE-2 cosmological zoom-in baryonic simulations, we examine the radial distributions of satellites with Mstar > 10^5 Msun around 8 isolated Milky Way- (MW) mass host galaxies and 4 hosts in LG-like pairs. We demonstrate that these simulations resolve the survival and physical destruction of satellites with Mstar >~ 10^5 Msun. The simulations broadly agree with LG observations, spanning the radial profiles around the MW and M31. This agreement does not depend strongly on satellite mass, even at distances <~ 100 kpc. Host-to-host variation dominates the scatter in satellite counts within 300 kpc of the hosts, while time variation dominates scatter within 50 kpc. More massive host galaxies within our sample have fewer satellites at small distances, likely because of enhanced tidal destruction of satellites via the baryonic disks of host galaxies. Furthermore, we quantify and provide fits to the tidal depletion of subhalos in baryonic relative to dark matter-only simulations as a function of distance. Our simulated profiles imply observational incompleteness in the LG even at Mstar >~ 10^5 Msun: we predict 2-10 such satellites to be discovered around the MW and possibly 6-9 around M31. To provide cosmological context, we compare our results with the radial profiles of satellites around MW analogs in the SAGA survey, finding that our simulations are broadly consistent with most SAGA systems. △ Less

Submitted 17 March, 2020; v1 submitted 25 April, 2019; originally announced April 2019.

Comments: 18 pages, 10 figures, plus appendices. Main results in figures 2, 3, and 4. Accepted version

Star formation histories of dwarf galaxies in the FIRE simulations: dependence on mass and Local Group environment

Authors: Shea Garrison-Kimmel, Andrew Wetzel, Phillip F. Hopkins, Robyn Sanderson, Kareem El-Badry, Andrew Graus, T. K. Chan, Robert Feldmann, Michael Boylan-Kolchin, Christopher Hayward, James S. Bullock, Alex Fitts, Jenna Samuel, Coral Wheeler, Dusan Keres, Claude-Andre Faucher-Giguere

Abstract: We study star formation histories (SFHs) of $\simeq500$ dwarf galaxies (stellar mass $M_\ast = 10^5 - 10^9\,M_\odot$) from FIRE-2 cosmological zoom-in simulations. We compare dwarfs around individual Milky Way (MW)-mass galaxies, dwarfs in Local Group (LG)-like environments, and true field (i.e. isolated) dwarf galaxies. We reproduce observed trends wherein higher-mass dwarfs quench later (if at a… ▽ More We study star formation histories (SFHs) of $\simeq500$ dwarf galaxies (stellar mass $M_\ast = 10^5 - 10^9\,M_\odot$) from FIRE-2 cosmological zoom-in simulations. We compare dwarfs around individual Milky Way (MW)-mass galaxies, dwarfs in Local Group (LG)-like environments, and true field (i.e. isolated) dwarf galaxies. We reproduce observed trends wherein higher-mass dwarfs quench later (if at all), regardless of environment. We also identify differences between the environments, both in terms of "satellite vs. central" and "LG vs. individual MWvs. isolated dwarf central." Around the individual MW-mass hosts, we recover the result expected from environmental quenching: central galaxies in the "near field" have more extended SFHs than their satellite counterparts, with the former more closely resemble isolated ("true field") dwarfs (though near-field centrals are still somewhat earlier forming). However, this difference is muted in the LG-like environments, where both near-field centrals and satellites have similar SFHs, which resemble satellites of single MW-mass hosts. This distinction is strongest for $M_\ast = 10^6 - 10^7\,M_\odot$ but exists at other masses. Our results suggest that the paired halo nature of the LG may regulate star formation in dwarf galaxies even beyond the virial radii of the MW and Andromeda. Caution is needed when comparing zoom-in simulations targeting isolated dwarf galaxies against observed dwarf galaxies in the LG. △ Less

Submitted 25 March, 2019; originally announced March 2019.

Comments: Main text: 11 pages, 8 figures; appendices: 4 pages, 4 figures. Submitted to MNRAS; comments welcome

Evolution of supernovae-driven superbubbles with conduction and cooling

Authors: Kareem El-Badry, Eve C. Ostriker, Chang-Goo Kim, Eliot Quataert, Daniel R. Weisz

Abstract: We use spherically symmetric hydrodynamic simulations to study the dynamical evolution and internal structure of superbubbles (SBs) driven by clustered supernovae (SNe), focusing on the effects of thermal conduction and cooling in the interface between the hot bubble interior and cooled shell. Our simulations employ an effective diffusivity to account for turbulent mixing from nonlinear instabilit… ▽ More We use spherically symmetric hydrodynamic simulations to study the dynamical evolution and internal structure of superbubbles (SBs) driven by clustered supernovae (SNe), focusing on the effects of thermal conduction and cooling in the interface between the hot bubble interior and cooled shell. Our simulations employ an effective diffusivity to account for turbulent mixing from nonlinear instabilities that are not captured in 1D. The conductive heat flux into the shell is balanced by a combination of cooling in the interface and evaporation of shell gas into the bubble interior. This evaporation increases the density, and decreases the temperature, of the SB interior by more than an order of magnitude relative to simulations without conduction. However, most of the energy conducted into the interface is immediately lost to cooling, reducing the evaporative mass flux required to balance conduction. As a result, the evaporation rate is typically a factor of $\sim$3-30 lower than predicted by the classical similarity solution of Weaver et al. (1977), which neglects cooling. Blast waves from the first $\sim$30 SNe remain supersonic in the SB interior because reduced evaporation from the interface lowers the mass they sweep up in the hot interior. Updating the Weaver solution to include cooling, we construct a new analytic model to predict the cooling rate, evaporation rate, and temporal evolution of SBs. The cooling rate, and hence the hot gas mass, momentum, and energy delivered by SBs, is set by the ambient ISM density and the efficiency of nonlinear mixing at the bubble/shell interface. △ Less

Submitted 1 October, 2019; v1 submitted 25 February, 2019; originally announced February 2019.

Comments: 26 pages, 14 figures, plus appendices. Accepted to MNRAS

AGN All the Way Down? AGN-like Line Ratios are Common In the Lowest-Mass Isolated Quiescent Galaxies

Authors: C. Dickey, M. Geha, A. Wetzel, K. El-Badry

Abstract: We investigate the lowest-mass quiescent galaxies known to exist in isolated environments ($\mathrm{M^* = 10^{9.0-9.5} \ M_\odot}$; 1.5 Mpc from a more massive galaxy). This population may represent the lowest stellar mass galaxies in which internal feedback quenches galaxy-wide star formation. We present Keck/ESI long-slit spectroscopy for 27 isolated galaxies in this regime: 20 quiescent galaxie… ▽ More We investigate the lowest-mass quiescent galaxies known to exist in isolated environments ($\mathrm{M^* = 10^{9.0-9.5} \ M_\odot}$; 1.5 Mpc from a more massive galaxy). This population may represent the lowest stellar mass galaxies in which internal feedback quenches galaxy-wide star formation. We present Keck/ESI long-slit spectroscopy for 27 isolated galaxies in this regime: 20 quiescent galaxies and 7 star-forming galaxies. We measure emission line strengths as a function of radius and place galaxies on the Baldwin Phillips Terlevich (BPT) diagram. Remarkably, 16 of 20 quiescent galaxies in our sample host central AGN-like line ratios. Only 5 of these quiescent galaxies were identified as AGN-like in SDSS due to lower spatial resolution and signal-to-noise. We find that many of the quiescent galaxies in our sample have spatially-extended emission across the non-SF regions of BPT-space. When considering only the central 1$^{\prime\prime}$, we identify a tight relationship between distance from the BPT star-forming sequence and host galaxy stellar age as traced by $\mathrm{D_n4000}$, such that older stellar ages are associated with larger distances from the star-forming locus. Our results suggest that the presence of hard ionizing radiation (AGN-like line ratios) is intrinsically tied to the quenching of what may be the lowest-mass self-quenched galaxies. △ Less

Submitted 4 February, 2019; originally announced February 2019.

Comments: 16 pages, 9 figures, submitted to ApJ

Dwarf Galaxies in CDM, WDM, and SIDM: Disentangling Baryons and Dark Matter Physics

Authors: Alex Fitts, Michael Boylan-Kolchin, Brandon Bozek, James S. Bullock, Andrew Graus, Victor Robles, Philip F. Hopkins, Kareem El-Badry, Shea Garrison-Kimmel, Claude-André Faucher-Giguère, Andrew Wetzel, Dušan Kereš

Abstract: We present a suite of FIRE-2 cosmological zoom-in simulations of isolated field dwarf galaxies, all with masses of $M_\mathrm{halo} \approx 10^{10}\,$M$_\odot$ at $z=0$, across a range of dark matter models. For the first time, we compare how both self-interacting dark matter (SIDM) and/or warm dark matter (WDM) models affect the assembly histories as well as the central density structure in fully… ▽ More We present a suite of FIRE-2 cosmological zoom-in simulations of isolated field dwarf galaxies, all with masses of $M_\mathrm{halo} \approx 10^{10}\,$M$_\odot$ at $z=0$, across a range of dark matter models. For the first time, we compare how both self-interacting dark matter (SIDM) and/or warm dark matter (WDM) models affect the assembly histories as well as the central density structure in fully hydrodynamical simulations of dwarfs. Dwarfs with smaller stellar half-mass radii (r$_{1/2}<500$ pc) have lower $σ_\star/V_\mathrm{max}$ ratios, reinforcing the idea that smaller dwarfs may reside in halos that are more massive than is naively expected. The majority of dwarfs simulated with self-interactions actually experience contraction of their inner density profiles with the addition of baryons relative to the cores produced in dark-matter-only runs, though the simulated dwarfs are always less centrally dense than in $Λ$CDM. The V$_{1/2}-$r$_{1/2}$ relation across all simulations is generally consistent with observations of Local Field dwarfs, though compact objects such as Tucana provide a unique challenge. Spatially-resolved rotation curves in the central regions ($<400$ pc) of small dwarfs could provide a way to distinguish between CDM, WDM, and SIDM, however: at the masses probed in this simulation suite, cored density profiles in dwarfs with small r$_{1/2}$ values can only originate from dark matter self-interactions. △ Less

Submitted 24 September, 2019; v1 submitted 28 November, 2018; originally announced November 2018.

Comments: 16 pages, 12 figures. V2: matches version accepted by MNRAS