The recent discovery of two detached black hole–star (BH–star) binaries from Gaia’s third data release has sparked interest in understanding the formation mechanisms of these systems. We investigate the formation of these systems by dynamical processes in young star clusters (SCs) and via isolated binary (IB) evolution, using a combination of direct
New mass estimates and cumulative mass profiles with Bayesian credible regions for the Milky Way (MW) are found using the Galactic Mass Estimator (GME) code and dwarf galaxy (DG) kinematic data from multiple sources. GME takes a hierarchical Bayesian approach to simultaneously estimate the true positions and velocities of the DGs, their velocity anisotropy, and the model parameters for the Galaxy’s total gravitational potential. In this study, we incorporate meaningful prior information from past studies and simulations. The prior distributions for the physical model are informed by the results of Eadie & Jurić, who used globular clusters instead of DGs, as well as by the subhalo distributions of the Ananke Gaia-like surveys from Feedback in Realistic Environments-2 cosmological simulations (see Sanderson et al.). Using DGs beyond 45 kpc, we report median and 95% credible region estimates for
- Award ID(s):
- 2045928
- NSF-PAR ID:
- 10361746
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 924
- Issue:
- 2
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 131
- Size(s):
- ["Article No. 131"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract N -body and population synthesis simulations. We find that dynamical formation in SCs is nearly 50 times more efficient per unit of star formation at producing BH–star binaries than IB evolution. We expand this analysis to the full Milky Way (MW) using a FIRE-2 hydrodynamical simulation of an MW-mass galaxy. Even assuming that only 10% of star formation goes into SCs, we find that approximately four out of every five BH–star systems are formed dynamically, and that the MW contains a total of ∼2 × 105BH–star systems. Many of these dynamically formed systems have longer orbital periods, greater eccentricities, and greater black hole masses than their isolated counterparts. For binaries older than 100 Myr, we show that any detectable system withe ≳ 0.5 orM BH≳ 10M ⊙canonly be formed through dynamical processes. Our MW model predicts between 64 and 215 such detections from the complete DR4 Gaia catalog, with the majority of systems being dynamically formed in massive and metal-rich SCs. Finally, we compare our populations to the recently discovered Gaia BH1 and Gaia BH2, and conclude that the dynamical scenario is the most favorable formation pathway for both systems. -
more » « less
Abstract M35 is a young open cluster and home to an extensive binary population. Using Gaia Data Release 3, Pan-STARRS, and Two Micron All Sky Survey photometry with the Bayesian statistical software, BASE-9, we derive precise cluster parameters, identify single and binary cluster members, and extract their masses. We identify 571 binaries down to Gaia
G = 20.3 and a lower limit on the binary frequency off b = 0.41 ± 0.02. We extend the binary demographics by many magnitudes faint-ward of previous (radial-velocity) studies of this cluster and further away from the cluster center (1.°78, roughly 10 core radii). We find the binary stars to be more centrally concentrated than the single stars in the cluster. Furthermore, we find strong evidence for mass segregation within the binary population itself, with progressively more-massive binary samples becoming more and more centrally concentrated. For the single stars, we find weaker evidence for mass segregation; only the most massive single stars (>2.5M ⊙) appear more centrally concentrated. Given the cluster age of ∼200 Myr, and our derived half-mass relaxation time for the cluster of 230 ± 84 Myr, we estimate ∼47% of the binary stars and ∼12% of single stars in the cluster have had time to become dynamically mass segregated. Importantly, when we investigate only stars with mass segregation timescales greater than the cluster age, we still find the binaries to be more centrally concentrated than the singles, suggesting the binaries may have formed with a primordially different spatial distribution from that of the single stars. -
Abstract We report the kinematic, orbital, and chemical properties of 12 stellar streams with no evident progenitors using line-of-sight velocities and metallicities from the Southern Stellar Stream Spectroscopic Survey ( S 5 ), proper motions from Gaia EDR3, and distances derived from distance tracers or the literature. This data set provides the largest homogeneously analyzed set of streams with full 6D kinematics and metallicities. All streams have heliocentric distances between ∼10 and 50 kpc. The velocity and metallicity dispersions show that half of the stream progenitors were disrupted dwarf galaxies (DGs), while the other half originated from disrupted globular clusters (GCs), hereafter referred to as DG and GC streams. Based on the mean metallicities of the streams and the mass–metallicity relation, the luminosities of the progenitors of the DG streams range between those of Carina and Ursa Major I (−9.5 ≲ M V ≲ −5.5). Four of the six GC streams have mean metallicities of [Fe/H] < −2, more metal poor than typical Milky Way (MW) GCs at similar distances. Interestingly, the 300S and Jet GC streams are the only streams on retrograde orbits in our dozen-stream sample. Finally, we compare the orbital properties of the streams with known DGs and GCs in the MW, finding several possible associations. Some streams appear to have been accreted with the recently discovered Gaia–Enceladus–Sausage system, and others suggest that GCs were formed in and accreted together with the progenitors of DG streams whose stellar masses are similar to those of Draco to Carina (∼10 5 –10 6 M ⊙ ).more » « less
-
more » « less
ABSTRACT We identify Local Group (LG) analogues in the IllustrisTNG cosmological simulation, and use these to study two-mass estimators for the LG: One based on the timing argument (TA) and one based on the virial theorem (VT). Including updated measurements of the Milky Way-M31 tangential velocity and the cosmological constant, we show that the TA mass estimator slightly overestimates the true median LG-mass, though the ratio of the TA to the true mass is consistent at the approximate 90 per cent c.l.. These are in broad agreement with previous results using dark matter-only simulations. We show that the VT estimator better estimates the true LG-mass, though there is a larger scatter in the virial mass to true mass ratio relative to the corresponding ratio for the TA. We attribute the broader scatter in the VT estimator to several factors, including the predominantly radial orbits for LG satellite galaxies, which differs from the VT assumption of isotropic orbits. With the systematic uncertainties we derive, the updated measurements of the LG mass at 90 per cent c.l. are $4.75_{-2.41}^{+2.22} \times 10^{12}$ M⊙ from the TA and $2.0_{-1.5}^{+2.1} \times 10^{12}$ M⊙ from the VT. We consider the LMC’s effect on the TA and VT LG mass estimates, and do not find exact LMC–MW–M31 analogues in the Illustris simulations. However, in LG simulations with satellite companions as massive as the LMC, we find that the effect on the TA and VT estimators is small, though we need further studies on a larger sample of LMC–MW–M31 systems to confirm these results.
-
more » « less
Abstract The shape and orientation of dark matter (DM) halos are sensitive to the microphysics of the DM particles, yet in many mass models, the symmetry axes of the Milky Way’s DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well motivated for the inner DM halo, but not for the outer halo. We use zoomed-in cosmological baryonic simulations from the Latte suite of FIRE-2 Milky Way–mass galaxies to explore the evolution of the DM halo’s orientation with radius and time, with or without a major merger with a Large Magellanic Cloud analog, and when varying the DM model. In three of the four cold DM halos we examine, the orientation of the halo minor axis diverges from the stellar disk vector by more than 20° beyond about 30 galactocentric kpc, reaching a maximum of 30°–90°, depending on the individual halo’s formation history. In identical simulations using a model of self-interacting DM with
σ = 1 cm2g−1, the halo remains aligned with the stellar disk out to ∼200–400 kpc. Interactions with massive satellites (M ≳ 4 × 1010M ⊙at pericenter;M ≳ 3.3 × 1010M ⊙at infall) affect the orientation of the halo significantly, aligning the halo’s major axis with the satellite galaxy from the disk to the virial radius. The relative orientation of the halo and disk beyond 30 kpc is a potential diagnostic of self-interacting DM, if the effects of massive satellites can be accounted for.
