Researcher profile

Kareem El-Badry

Kareem El-Badry contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
31works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2025arXiv

Millions of Main-Sequence Binary Stars from Gaia BP/RP Spectra

We present the main-sequence binary (MSMS) Catalog derived from Gaia Data Release 3 BP/RP (XP) spectra. Leveraging the vast sample of low-resolution Gaia XP spectra, we develop a forward modeling approach that maps stellar mass and photometric metallicity to XP spectra using a neural network. Our methodology identifies binary systems through statistical comparison of single- and binary-star model fits, enabling detection of binaries with mass ratios between 0.4 and 1.0 and flux ratios larger than 0.1. From an initial sample of 35 million stars within 1 kpc, we identify 14 million binary candidates and define a high-confidence "golden sample" of 1 million binary systems. This large, homogeneous sample enables detailed statistical analysis of binary properties across diverse Galactic environments, providing new insights into binary star formation and evolution. In addition, the $χ^2$ comparison allows us to distinguish stars with luminous companions from single stars or binaries with dark companions, such as white dwarfs, neutron stars and black hole candidates, improving our understanding of compact object populations.

preprint2023arXiv

Science with the Ultraviolet Explorer (UVEX)

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.

preprint2023arXiv

Wide post-common envelope binaries containing ultramassive white dwarfs: evidence for efficient envelope ejection in massive AGB stars

Post-common-envelope binaries (PCEBs) containing a white dwarf (WD) and a main-sequence (MS) star can constrain the physics of common envelope evolution and calibrate binary evolution models. Most PCEBs studied to date have short orbital periods ($P_{\rm orb} \lesssim 1\,$d), implying relatively inefficient harnessing of binaries' orbital energy for envelope expulsion. Here, we present follow-up observations of five binaries from {\it Gaia} DR3 containing solar-type MS stars and probable ultramassive WDs ($M\gtrsim 1.2\,M_{\odot}$) with significantly wider orbits than previously known PCEBs, $P_{\rm orb} = 18-49\,$d. The WD masses are much higher than expected for systems formed via stable mass transfer at these periods, and their near-circular orbits suggest partial tidal circularization when the WD progenitors were giants. These properties strongly suggest that the binaries are PCEBs. Forming PCEBs at such wide separations requires highly efficient envelope ejection, and we find that the observed periods can only be explained if a significant fraction of the energy released when the envelope recombines goes into ejecting it. Our 1D stellar models including recombination energy confirm prior predictions that a wide range of PCEB orbital periods, extending up to months or years, can potentially result from Roche lobe overflow of a luminous AGB star. This evolutionary scenario may also explain the formation of several wide WD+MS binaries discovered via self-lensing, as well as a significant fraction of post-AGB binaries and barium stars.

preprint2022arXiv

A Panchromatic Study of Massive Stars in the Extremely Metal-Poor Local Group Dwarf Galaxy Leo A

We characterize massive stars (M>8 M_sun) in the nearby (D~0.8 Mpc) extremely metal-poor (Z~5% Z_sun) galaxy Leo A using Hubble Space Telescope ultra-violet (UV), optical, and near-infrared (NIR) imaging along with Keck/LRIS and MMT/Binospec optical spectroscopy for 18 main sequence OB stars. We find that: (a) 12 of our 18 stars show emission lines, despite not being associated with an H II region, suggestive of stellar activity (e.g., mass loss, accretion, binary star interaction), which is consistent with previous predictions of enhanced activity at low metallicity; (b) 6 are Be stars, which are the first to be spectroscopically studied at such low metallicity -- these Be stars have unusual panchromatic SEDs; (c) for stars well-fit by the TLUSTY non-local thermodynamic equilibrium (non-LTE) models, the photometric and spectroscopic values of T_eff and log(g) agree to within ~0.01 dex and ~0.18 dex, respectively, indicating that NUV/optical/NIR imaging can be used to reliably characterize massive (M ~ 8-30 M_sun) main sequence star properties relative to optical spectroscopy; (d) the properties of the most massive stars in H II regions are consistent with constraints from previous nebular emission line studies; and (e) 13 stars with M>8 M_sun are >40 pc from a known star cluster or H II region. Our sample comprises ~50% of all known massive stars at Z < 10% Z_sun with derived stellar parameters, high-quality optical spectra, and panchromatic photometry.

preprint2022arXiv

An X-ray quiet black hole born with a negligible kick in a massive binary within the Large Magellanic Cloud

Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5-sigma confidence level. The minimum companion mass implies that it is a black hole. No other X-ray quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the Cosmos.

preprint2022arXiv

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

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 &#34;rotating cooling flow&#34; 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.

preprint2022arXiv

Improving white dwarfs as chronometers with Gaia parallaxes and spectroscopic metallicities

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

The gravitational redshift of solar-type stars from Gaia DR3 wide binaries

Light escaping from a gravitational potential suffers a redshift with magnitude proportional to the depth of the potential. This &#34;gravitational redshift&#34; is easily measurable in dense stars such as white dwarfs, but is much weaker and has evaded unambiguous detection in main-sequence stars. I show that the effect is directly measurable in the Gaia DR3 radial velocities (RVs) of the components of wide binary stars. In a sample of $\sim$500 wide binaries containing a solar-type main-sequence star and a red giant or red clump companion, the apparent RV of the giant is on average $0.49 \pm 0.02 \,\, \rm km\,s^{-1}$ lower than that of the main-sequence star. This owes primarily to the giants&#39; weaker gravitational fields and is in reasonably good agreement with the value expected from general relativity.

preprint2022arXiv

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

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 &#34;non-disky&#34; 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&#39;s intrinsic shape -- which is a direct quantification of the distribution of the galaxy&#39;s stellar content -- to interrogate the feedback implementation in simulated galaxies.

preprint2022arXiv

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

We analyze two binary systems containing giant stars, V723 Mon (&#34;the Unicorn&#34;) and 2M04123153+6738486 (&#34;the Giraffe&#34;). 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.

preprint2022arXiv

What are the spectroscopic binaries with high mass functions near the Gaia DR3 main sequence?

The 3rd data release of the Gaia mission includes orbital solutions for $> 10^5$ single-lined spectroscopic binaries, representing more than an order of magnitude increase in sample size over all previous studies. This dataset is a treasure trove for searches for quiescent black hole + normal star binaries. We investigate one population of black hole candidate binaries highlighted in the data release: sources near the main sequence in the color-magnitude diagram (CMD) with dynamically-inferred companion masses $M_2$ larger than the CMD-inferred mass of the luminous star. We model light curves, spectral energy distributions, and archival spectra of the 14 such objects in DR3 with high-significance orbital solutions and inferred $M_2 > 3\,M_{\odot}$. We find that 100\% of these sources are mass-transfer binaries containing a highly stripped lower giant donor ($0.2 \lesssim M/M_{\odot} \lesssim 0.4$) and much more massive ($2 \lesssim M/M_{\odot} \lesssim 2.5$) main-sequence accretor. The Gaia orbital solutions are for the donors, which contribute about half the light in the Gaia RVS bandpass but only $\lesssim 20\%$ in the $g-$band. The accretors&#39; broad spectral features likely prevented the sources from being classified as double-lined. The donors are all close to Roche lobe-filling ($R/R_{\rm Roche\,lobe}>0.8$), but modeling suggests that a majority are detached ($R/R_{\rm Roche\,lobe}<1$). Binary evolution models predict that these systems will soon become detached helium white dwarf + main sequence &#34;EL CVn&#34; binaries. Our investigation highlights both the power of Gaia data for selecting interesting sub-populations of binaries and the ways in which binary evolution can bamboozle standard CMD-based stellar mass estimates.

preprint2022arXiv

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

Due to the different environments in the Milky Way&#39;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.

preprint2022arXiv

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

The Gaia mission recently revealed an excess population of equal-mass &#34;twin&#34; 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&#39; 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.

preprint2021arXiv

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

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&#39; 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.

preprint2021arXiv

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

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&#39;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&#39;s motion. A stripped-star model reproduces the observed luminosity of the system, while a normal star with the B star&#39;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.

preprint2021arXiv

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

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&#39;s observed radius and the orbital period -- density relation for Roche lobe-filling stars: the star&#39;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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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 &#34;core-Einasto&#34; 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&#39;s dark matter halo is consistent with this expectation.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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 &#34;puffing&#34; 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 &#34;puffing&#34; 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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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 &#39;field stars&#39;. 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 &#39;dissolution&#39; of stellar birth associations into the field.

preprint2020arXiv

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

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&#39;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 &#34;normal&#34; 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.

preprint2020arXiv

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

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 &#34;formation&#34; 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.

preprint2019arXiv

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

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.

preprint2019arXiv

The Separation Distribution of Ultra-Wide Binaries across Galactic Populations

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&#39; 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&#39; 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&#39; 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&#39; birth clusters.