ABSTRACT
White dwarf binaries with orbital periods below 1 h will be the most numerous sources for the space-based gravitational wave detector Laser Interferometer Space Antenna (LISA). Based on ...thousands of individually resolved systems, we will be able to constrain binary evolution and provide a new map of the Milky Way and its close surroundings. In this paper we predict the main properties of populations of different types of detached white dwarf binaries detected by LISA over time. For the first time, we combine a high-resolution cosmological simulation of a Milky Way-mass galaxy (taken from the FIRE project) with a binary population synthesis model for low- and intermediate-mass stars. Our Galaxy model therefore provides a cosmologically realistic star formation and metallicity history for the Galaxy and naturally produces its different components such as the thin and thick disc, the bulge, the stellar halo, and satellite galaxies and streams. Thanks to the simulation, we show how different Galactic components contribute differently to the gravitational wave signal, mostly due to their typical age and distance distributions. We find that the dominant LISA sources will be He–He double white dwarfs (DWDs) and He–CO DWDs with important contributions from the thick disc and bulge. The resulting sky map of the sources is different from previous models, with important consequences for the searches for electromagnetic counterparts and data analysis. We also emphasize that much of the science-enabling information regarding white dwarf binaries, such as the chirp mass and the sky localization, becomes increasingly rich with long observations, including an extended mission up to 8 yr.
Abstract Short-period double white dwarf (DWD) binaries will be the most prolific source of gravitational waves (GWs) for the Laser Interferometer Space Antenna (LISA). DWDs with GW frequencies below ...∼1 mHz will be the dominant contributor to a stochastic foreground caused by overlapping GW signals. Population modeling of Galactic DWDs typically assumes a binary fraction of 50% and a log-uniform zero-age main sequence (ZAMS) orbital period distribution. However, recent observations have shown that the binary fraction of close, solar-type stars exhibits a strong anticorrelation with metallicity that modulates the ZAMS orbital period distribution below 10 4 days. In this study, we perform the first simulation of the Galactic DWD population observable by LISA that incorporates an empirically derived metallicity-dependent binary fraction, using the binary population synthesis suite COSMIC and a metallicity-dependent star formation history. We compare two models: one that assumes a metallicity-dependent binary fraction, and one with a binary fraction of 50%. We repeat our analysis for three different assumptions for Roche-lobe overflow interactions. We find that while metallicity impacts the evolution and intrinsic properties of our simulated DWD progenitor binaries, the LISA-resolvable populations of the two models remain roughly indistinguishable. However, the size of the total Galactic DWD population orbiting in the LISA frequency band is reduced by more than half when accounting for a metallicity-dependent binary fraction for two of our four variations, which also lowers the effective foreground. The LISA population remains unchanged in number for two variations, highlighting the sensitivity of the population to binary evolution prescriptions.
Stellar streams result from the tidal disruption of satellites and star clusters as they orbit a host galaxy, and can be very sensitive probes of the gravitational potential of the host system. We ...select and study narrow stellar streams formed in a Milky-Way-like dark matter halo of the Aquarius suite of cosmological simulations, to determine if these streams can be used to constrain the present day characteristic parameters of the halo's gravitational potential. We find that orbits integrated in both spherical and triaxial static Navarro-Frenk-White potentials reproduce the locations and kinematics of the various streams reasonably well. To quantify this further, we determine the best-fit potential parameters by maximizing the amount of clustering of the stream stars in the space of their actions. We show that using our set of Aquarius streams, we recover a mass profile that is consistent with the spherically averaged dark matter profile of the host halo, although we ignored both triaxiality and time evolution in the fit. This gives us confidence that such methods can be applied to the many streams that will be discovered by the Gaia mission to determine the gravitational potential of our Galaxy.
Abstract
Detection of black holes (BHs) with detached luminous companions (LCs) can be instrumental in connecting the BH properties with their progenitors since the latter can be inferred from the ...observable properties of the LC. Past studies showed the promise of Gaia astrometry in detecting BH–LC binaries. We build on these studies by (1) initializing the zero-age binary properties based on realistic, metallicity-dependent star formation history in the Milky Way (MW); (2) evolving these binaries to current epoch to generate realistic MW populations of BH–LC binaries; (3) distributing these binaries in the MW, preserving the complex age–metallicity-Galactic position correlations; (4) accounting for extinction and reddening using three-dimensional dust maps; and (5) examining the extended Gaia mission’s ability to resolve BH–LC binaries. We restrict ourselves to detached BH–LC binaries with orbital period
P
orb
≤ 10 yr such that Gaia can observe at least one full orbit. We find that (1) the extended Gaia mission can astrometrically resolve ∼30–300 detached BH–LC binaries depending on our assumptions of supernova physics and astrometric detection threshold; (2) Gaia’s astrometry alone can indicate BH candidates for ∼10–100 BH–LC binaries by constraining the dark
primary
mass ≥3
M
⊙
; and (3) distributions of observables, including orbital periods, eccentricities, and component masses, are sensitive to the adopted binary evolution model and hence can directly inform binary evolution models. Finally, we comment on the potential to further characterize these BH binaries through radial velocity measurements and observation of X-ray counterparts.
With Gaia Data Release 2, the astronomical community is entering a new era of multidimensional surveys of the Milky Way. This new phase-space view of our Galaxy demands new tools for comparing ...observations to simulations of Milky Way-mass galaxies in a cosmological context, to test the physics of both dark matter and galaxy formation. We present ananke, a framework for generating synthetic phase-space surveys from high-resolution baryonic simulations, and use it to generate a suite of synthetic surveys resembling Gaia DR2 in data structure, magnitude limits, and observational errors. We use three cosmological simulations of Milky Way-mass galaxies from the Latte suite of the Feedback In Realistic Environments project, which feature self-consistent clustering of star formation in dense molecular clouds and thin stellar/gaseous disks in live cosmological halos with satellite dwarf galaxies and stellar halos. We select three solar viewpoints from each simulation to generate nine synthetic Gaia-like surveys. We sample synthetic stars by assuming each star particle (of mass 7070 M ) represents a single stellar population. At each viewpoint, we compute dust extinction from the simulated gas metallicity distribution and apply a simple error model to produce a synthetic Gaia-like survey that includes both observational properties and a pointer to the generating star particle. We provide the complete simulation snapshot at z = 0 for each simulated galaxy. We describe data access points, the data model, and plans for future upgrades. These synthetic surveys provide a tool for the scientific community to test analysis methods and interpret Gaia data.
Abstract
We analyse stellar streams in action-angle coordinates combined with recent local direct acceleration measurements to provide joint constraints on the potential of our galaxy. Our stream ...analysis uses the Kullback–Leibler divergence with a likelihood analysis based on the two-point correlation function. We provide joint constraints from pulsar accelerations and stellar streams for local and global parameters that describe the potential of the Milky Way (MW). Our goal is to build an “acceleration ladder,” where direct acceleration measurements that are currently limited in dynamic range are combined with indirect techniques that can access a much larger volume of the MW. To constrain the MW potential with stellar streams, we consider the Palomar 5, Orphan, Nyx, Helmi, and GD1 streams. Of the potential models that we have considered here, the preferred potential for the streams is a two-component Staeckel potential. We also compare the vertical accelerations from stellar streams and pulsar timing, defining the function
f
(
z
)
=
α
1
pulsar
z
−
∂
Φ
∂
z
, where Φ is the MW potential determined from stellar streams and
α
1 pulsar
z
is the vertical acceleration determined from pulsar timing observations. Our analysis indicates that the Oort limit determined from streams is consistently (regardless of the choice of potential) lower than that determined from pulsar timing observations. The calibration we have derived here may be used to correct the estimate of the acceleration from stellar streams.
ABSTRACT
The orbits of satellite galaxies encode rich information about their histories. We investigate the orbital dynamics and histories of satellite galaxies around Milky Way (MW)-mass host ...galaxies using the FIRE-2 cosmological simulations, which, as previous works have shown, produce satellite mass functions and spatial distributions that broadly agree with observations. We first examine trends in orbital dynamics at z = 0, including total velocity, specific angular momentum, and specific total energy: the time of infall into the MW-mass halo primarily determines these orbital properties. We then examine orbital histories, focusing on the lookback time of first infall into a host halo and pericentre distances, times, and counts. Roughly 37 per cent of galaxies with $M_{\rm star}\lesssim 10^7\, {\rm M}_{\odot }$ were ‘pre-processed’ as a satellite in a lower-mass group, typically $\approx 2.7\, {\rm Gyr}$ before falling into the MW-mass halo. Half of all satellites at z = 0 experienced multiple pericentres about their MW-mass host. Remarkably, for most (67 per cent) of these satellites, their most recent pericentre was not their minimum pericentre: the minimum typically was ∼40 per cent smaller and occurred $\sim 6\, {\rm Gyr}$ earlier. These satellites with growing pericentres appear to have multiple origins: for about half, their specific angular momentum gradually increased over time, while for the other half, most rapidly increased near their first apocentre, suggesting that a combination of a time-dependent MW-mass halo potential and dynamical perturbations in the outer halo caused these satellites’ pericentres to grow. Our results highlight the limitations of idealized, static orbit modelling, especially for pericentre histories.
The class of tidal features around galaxies known as 'shells' or 'umbrellas' comprises debris that has arisen from high-mass-ratio mergers with low-impact parameter; the nearly radial orbits of the ...debris give rise to a unique morphology, a universal density profile and a tight correlation between positions and velocities of the material. As such they are accessible to analytical treatment, and can provide a relatively clean system for probing the gravitational potential of the host galaxy. In this work, we present a simple analytical model that describes the density profile, phase-space distribution, and geometry of a shell and whose parameters are directly related to physical characteristics of the interacting galaxies. The model makes three assumptions: the orbit of the interacting galaxies is radial, the potential of the host galaxy at the shell radius is spherical and the satellite galaxy's initial velocity distribution is Maxwellian. We quantify the error introduced by the first two assumptions and show that selecting shells by their appearance on the sky is a sufficient basis to assume that these simplifications are valid. We further demonstrate that (1) given only an image of a shell, the radial gravitational force at the shell edge and the phase-space density of the satellite are jointly constrained, (2) combining the image with measurements of either point line-of-sight velocities or integrated-light spectra will yield an independent estimate of the gravitational force at a shell and (3) an independent measurement of this force is obtained for each shell observed around a given galaxy, potentially enabling a determination of the galactic mass distribution.
Abstract
We present a new set of high-resolution hydrodynamic cosmological zoom-in simulations that apply the Feedback In Realistic Environments physics to both Local Group (LG)-like and isolated ...Milky Way (MW)-like volumes (10 host systems in total with a baryonic particle mass ${\simeq } 3500\!-\!7000\,{\rm M}_\odot$). We study the stellar mass functions, circular velocity or mass profiles, and velocity dispersions of the dwarf galaxy populations. The simulations reproduce the stellar mass function and central densities of MW satellite dwarfs for $M_\ast \ge 10^{5.5}\,{\rm M}_\odot$ and predict the existence of ∼3 unidentified galaxies with $M_\ast \sim 10^5\,{\rm M}_\odot$ within 300 kpc of the MW. Overall, we find no evidence for the classical missing satellites or too-big-to-fail (TBTF) problems for satellite galaxies in our sample. Among the satellites, TBTF is resolved primarily by subhalo disruption and overall mass-loss; central density profiles of subhaloes are of secondary importance. For non-satellite galaxies, our LG-like simulations predict as many as ∼10 as-of-yet unseen galaxies at distances $0.3\!-\!1\, \mathrm{Mpc}$ from both hosts, with $M_\ast \simeq 10^{5-6}\,{\rm M}_\odot$ (in haloes with Vmax ∼ 20 km s−1), albeit with large halo-to-halo variance. None of our simulations produces a compact, baryon-dominated, high-density dwarf elliptical-type galaxy (with Vcirc ≳ 35 km s−1 at $r\lt 1\,$ kpc), of which six may appear in the LG (but none in the MW). It may therefore remain a challenge to reproduce the full diversity of the dwarf population, including both the highest and lowest density systems.
We present a new method for constraining the Milky Way halo gravitational potential by simultaneously fitting multiple tidal streams. This method requires three-dimensional positions and velocities ...for all stars to be fit, but does not require identification of any specific stream or determination of stream membership for any star. We exploit the principle that the action distribution of stream stars is most clustered when the potential used to calculate the actions is closest to the true potential. Clustering is quantified with the Kullback-Leibler Divergence (KLD), which also provides conditional uncertainties for our parameter estimates. We show, for toy Gaia-like data in a spherical isochrone potential, that maximizing the KLD of the action distribution relative to a smoother distribution recovers the input potential. The precision depends on the observational errors and number of streams; using K III giants as tracers, we measure the enclosed mass at the average radius of the sample stars accurate to 3% and precise to 20%-40%. Recovery of the scale radius is precise to 25%, biased 50% high by the small galactocentric distance range of stars in our mock sample (1-25 kpc, or about three scale radii, with mean 6.5 kpc). 20-25 streams with at least 100 stars each are required for a stable confidence interval. With radial velocities (RVs) to 100 kpc, all parameters are determined with ~10% accuracy and 20% precision (1.3% accuracy for the enclosed mass), underlining the need to complete the RV catalog for faint halo stars observed by Gaia.