Researcher profile

Michael Boylan-Kolchin

Michael Boylan-Kolchin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
31works
0followers
8topics
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)

preprint2024arXiv

An analytic surface density profile for $Λ$CDM halos and gravitational lensing studies

We introduce an analytic surface density profile for dark matter halos that accurately reproduces the structure of simulated halos of mass $M_{\rm vir} = 10^{7-11}\ M_\odot$, making it useful for modeling line-of-sight perturbers in strong gravitational lensing models. The two-parameter function has an analytic deflection potential and is more accurate than the projected Navarro, Frenk & White (NFW) profile commonly adopted at this mass scale for perturbers, especially at the small radii of most relevant for lensing perturbations. Using a characteristic radius, $R_{-1}$, where the log slope of surface density is equal to $-1$, and an associated surface density, $Σ_{-1}$, we can represent the expected lensing signal from line-of-sight halos statistically, for an ensemble of halo orientations, using a distribution of {\em projected concentration} parameters, $\mathcal{C}_{\rm vir} := r_{\rm vir}/ R_{-1}$. Though an individual halo can have a projected concentration that varies with orientation with respect to the observer, the range of projected concentrations correlates with the usual three-dimensional halo concentration in a way that enables ease of use.

preprint2023arXiv

On the Cosmic Web Elongation in Fuzzy Dark Matter Cosmologies: Effects on Density Profiles, Shapes and Alignments of Halos

The fuzzy dark matter (FDM) scenario has received increased attention in recent years due to the small-scale challenges of the vanilla Lambda cold dark matter ($Λ$CDM) cosmological model and the lack of any experimental evidence for any candidate particle. In this study, we use cosmological $N$-body simulations to investigate high-redshift dark matter halos and their responsiveness to an FDM-like power spectrum cutoff on small scales in the primordial density perturbations. We study halo density profiles, shapes and alignments in FDM-like cosmologies (the latter two for the first time) by providing fits and quantifying departures from $Λ$CDM as a function of the particle mass $m$. Compared to $Λ$CDM, the concentrations of FDM-like halos are lower, peaking at an $m$-dependent halo mass and thus breaking the approximate universality of density profiles in $Λ$CDM. The intermediate-to-major and minor-to-major shape parameter profiles are monotonically increasing with ellipsoidal radius in $N$-body simulations of $Λ$CDM. In FDM-like cosmologies, the monotonicity is broken, halos are more elongated around the virial radius than their $Λ$CDM counterparts and less elongated closer to the center. Finally, intrinsic alignment correlations, stemming from the deformation of initially spherically collapsing halos in an ambient gravitational tidal field, become stronger with decreasing $m$. At $z\sim 4$, we find a $6.4 σ$-significance in the fractional differences between the isotropised linear alignment magnitudes $D_{\text{iso}}$ in the $m=10^{-22}$ eV model and $Λ$CDM. Such FDM-like imprints on the internal properties of virialised halos are expected to be strikingly visible in the high-$z$ Universe.

preprint2023arXiv

The JWST Resolved Stellar Populations Early Release Science Program II. Survey Overview

We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color-magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen burning limit in M92 ($<0.08$ $M_{\odot}$; SNR $\sim5$ at $m_{F090W}\sim28.2$; $M_{F090W}\sim+13.6$), include the lowest-mass stars observed outside the Milky Way in Draco II (0.09 $M_{\odot}$; SNR $=10$ at $m_{F090W}\sim29$; $M_{F090W}\sim+12.1$), and reach $\sim1.5$ magnitudes below the oldest main sequence turnoff in WLM (SNR $=10$ at $m_{F090W}\sim29.5$; $M_{F090W}\sim+4.6$). The PARSEC stellar models provide a good qualitative match to the NIRCam CMDs, though are $\sim0.05$ mag too blue compared to M92 F090W$-$F150W data. The NIRCam ETC (v2.0) matches the SNRs based on photon noise from DOLPHOT stellar photometry in uncrowded fields, but the ETC may not be accurate in more crowded fields, similar to what is known for HST. We release beta versions of DOLPHOT NIRCam and NIRISS modules to the community. Results from this ERS program will establish JWST as the premier instrument for resolved stellar populations studies for decades to come.

preprint2022arXiv

Amplified J-factors in the Galactic Center for velocity-dependent darkmatter annihilation in FIRE simulations

We use FIRE-2 zoom cosmological simulations of Milky Way size galaxy halos to calculate astrophysical J-factors for dark matter annihilation and indirect detection studies. In addition to velocity-independent (s-wave) annihilation cross sections $σ_v$, we also calculate effective J-factors for velocity-dependent models, where the annihilation cross section is either either p-wave ($\propto v^2/c^2$) or d-wave ($\propto v^4/c^4$). We use 12 pairs of simulations, each run with dark-matter-only (DMO) physics and FIRE-2 physics. We observe FIRE runs produce central dark matter velocity dispersions that are systematically larger than in DMO runs by factors of $\sim 2.5-4$. They also have a larger range of central ($\sim 400$ pc) dark matter densities than the DMO runs ($ρ_{\rm FIRE}/ρ_{\rm DMO} \simeq 0.5 - 3$) owing to the competing effects of baryonic contraction and feedback. At 3 degrees from the Galactic Center, FIRE J-factors are $5-50$ (p-wave) and $15-500$ (d-wave) times higher than in the DMO runs. The change in s-wave signal at 3 degrees is more modest and can be higher or lower ($\sim 0.3-6$), though the shape of the emission profile is flatter (less peaked towards the Galactic Center) and more circular on the sky in FIRE runs. Our results for s-wave are broadly consistent with the range of assumptions in most indirect detection studies. We observe p-wave J-factors that are significantly enhanced compared to most past estimates. We find that thermal models with p-wave annihilation may be within range of detection in the near future.

preprint2022arXiv

Dissipative Dark Matter on FIRE: II. Observational signatures and constraints from local dwarf galaxies

We analyze the first set of cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modelled as a constant fractional energy dissipation ($f_{\rm diss}=0.5$). In this paper, we examine the properties of dwarf galaxies with $M_{\ast} \sim 10^{5}\operatorname{-}10^{9}\,{\rm M}_{\odot}$ in both isolation and within Milky Way-mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of stellar/neutral gas disk formation in dSIDM with $(σ/m)\leq 1\,{\rm cm^2\,g^{-1}}$ but they are still consistent with observed galaxy sizes and masses. In addition, as a result of the steeper central density profiles developed in dSIDM, the sub-kpc circular velocities of isolated dwarfs in models with $(σ/m)\geq 0.1\,{\rm cm^2\,g^{-1}}$ are enhanced by about a factor of two, which are still consistent with the measured stellar velocity dispersions of Local Group dwarfs but in tension with the HI rotation curves of more massive field dwarfs. Meanwhile, for satellites of the simulated Milky Way-mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help address the missing compact dwarf satellites in CDM. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with constant cross-section $(σ/m) \gtrsim 0.1\,{\rm cm^2\,g^{-1}}$ (assuming $f_{\rm diss}=0.5$) are effectively ruled out in bright dwarfs ($M_{\rm halo}\sim 10^{11}\,{\rm M}_{\odot}$) by circular velocity constraints. However, models with lower effective cross-sections (at this halo mass/velocity scale) are still viable and can give rise to non-trivial observable signatures.

preprint2022arXiv

Extinguishing the FIRE: environmental quenching of satellite galaxies around Milky Way-mass hosts in simulations

The star formation and gas content of satellite galaxies around the Milky Way (MW) and Andromeda (M31) are depleted relative to more isolated galaxies in the Local Group (LG) at fixed stellar mass. We explore the environmental regulation of gas content and quenching of star formation in $z=0$ galaxies at $M*=10^{5-10}\rm{M}_{\odot}$ around 14 MW-mass hosts from the FIRE-2 simulations. Lower-mass satellites ($M*\lesssim10^7\rm{M}_{\odot}$) are mostly quiescent and higher-mass satellites ($M*\gtrsim10^8\rm{M}_{\odot}$) are mostly star-forming, with intermediate-mass satellites ($M*\approx10^{7-8}\rm{M}_{\odot}$) split roughly equally between quiescent and star-forming. Hosts with more gas in their circumgalactic medium have a higher quiescent fraction of massive satellites ($M*=10^{8-9}\rm{M}_{\odot}$). We find no significant dependence on isolated versus paired (LG-like) host environments, and the quiescent fractions of satellites around MW-mass and LMC-mass hosts from the FIRE-2 simulations are remarkably similar. Environmental effects that lead to quenching can also occur as preprocessing in low-mass groups prior to MW infall. Lower-mass satellites typically quenched before MW infall as central galaxies or rapidly during infall into a low-mass group or a MW-mass galaxy. Most intermediate- to high-mass quiescent satellites have experienced $\geq1-2$ pericentre passages ($\approx2.5-5$ Gyr) within a MW-mass halo. Most galaxies with $M*\gtrsim10^{6.5}\rm{M}_{\odot}$ did not quench before falling into a host, indicating a possible upper mass limit for isolated quenching. The simulations reproduce the average trend in the LG quiescent fraction across the full range of satellite stellar masses. Though the simulations are consistent with the SAGA survey&#39;s quiescent fraction at $M*\gtrsim10^8\rm{M}_{\odot}$, they do not generally reproduce SAGA&#39;s turnover at lower masses.

preprint2022arXiv

Galaxies lacking dark matter produced by close encounters in a cosmological simulation

The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark matter and modified gravity. Apprehension grew after several cosmological simulations using the conventional model failed to form adequate numerical analogues with comparable internal characteristics (stellar masses, sizes, velocity dispersions and morphologies). Here we show that the standard paradigm naturally produces galaxies lacking dark matter with internal characteristics in agreement with observations. Using a state-of-the-art cosmological simulation and a meticulous galaxy-identification technique, we find that extreme close encounters with massive neighbours can be responsible for this. We predict that approximately 30 percent of massive central galaxies (with at least 1e11 solar masses in stars) harbour at least one dark-matter-deficient satellite (with 1e8 - 1e9 solar masses in stars). This distinctive class of galaxies provides an additional layer in our understanding of the role of interactions in shaping galactic properties. Future observations surveying galaxies in the aforementioned regime will provide a crucial test of this scenario.

preprint2022arXiv

Great Balls of FIRE I: The formation of star clusters across cosmic time in a Milky Way-mass galaxy

The properties of young star clusters formed within a galaxy are thought to vary in different interstellar medium (ISM) conditions, but the details of this mapping from galactic to cluster scales are poorly understood due to the large dynamic range involved in galaxy and star cluster formation. We introduce a new method for modeling cluster formation in galaxy simulations: mapping giant molecular clouds (GMCs) formed self-consistently in a FIRE-2 MHD galaxy simulation onto a cluster population according to a GMC-scale cluster formation model calibrated to higher-resolution simulations, obtaining detailed properties of the galaxy&#39;s star clusters in mass, metallicity, space, and time. We find $\sim 10\%$ of all stars formed in the galaxy originate in gravitationally-bound clusters overall, and this fraction increases in regions with elevated $Σ_{\rm gas}$ and $Σ_{\rm SFR}$, because such regions host denser GMCs with higher star formation efficiency. These quantities vary systematically over the history of the galaxy, driving variations in cluster formation. The mass function of bound clusters varies -- no single Schechter-like or power-law distribution applies at all times. In the most extreme episodes, clusters as massive as $7\times 10^6 M_\odot$ form in massive, dense clouds with high star formation efficiency. The initial mass-radius relation of young star clusters is consistent with an environmentally-dependent 3D density that increases with $Σ_{\rm gas}$ and $Σ_{\rm SFR}$. The model does not reproduce the age and metallicity statistics of old ($>11\rm Gyr$) globular clusters found in the Milky Way, possibly because it forms stars more slowly at $z>3$.

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

On the probability of the extremely lensed $z$=6.2 Earendel source being a Population~III star

The recent discovery of the extremely lensed Earendel object at $z=6.2$ is remarkable in that it is likely a single star or stellar multiple, observed within the first billion years of cosmic history. Depending on its mass, which is still uncertain but will soon be more tightly constrained with the James Webb Space Telescope, the Earendel star might even be a member of the first generation of stars, the so-called Population~III (Pop~III). By combining results from detailed cosmological simulations of the assembly of the first galaxies, including the enrichment of the pristine gas with heavy chemical elements, with assumptions on key stellar parameters, we quantify the probability that Earendel has indeed a Pop~III origin. We find that this probability is non-negligible throughout the mass range inferred for Earendel, specifically ranging from a few percent at the lower-mass end to near unity for some Pop~III initial mass function (IMF) models towards the high-mass end of the allowed range. For models that extend the metal-enriched IMF to $500$\,M$_\odot$, the likelihood of Earendel being a Pop~III star stays at the few to ten percent level. We discuss the implications of such a discovery for the overall endeavor to probe the hitherto so elusive first stars in the universe.

preprint2022arXiv

Shapes of Milky-Way-Mass Galaxies with Self-Interacting Dark Matter

Self-interacting dark matter (SIDM) models offer one way to reconcile inconsistencies between observations and predictions from collisionless cold dark matter (CDM) models on dwarf-galaxy scales. In order to incorporate the effects of both baryonic and SIDM interactions, we study a suite of cosmological-baryonic simulations of Milky-Way (MW)-mass galaxies from the Feedback in Realistic Environments (FIRE-2) project where we vary the SIDM self-interaction cross-section $σ/m$. We compare the shape of the main dark matter (DM) halo at redshift $z=0$ predicted by SIDM simulations (at $σ/m=0.1$, $1$, and $10$ cm$^2$ g$^{-1}$) with CDM simulations using the same initial conditions. In the presence of baryonic feedback effects, we find that SIDM models do not produce the large differences in the inner structure of MW-mass galaxies predicted by SIDM-only models. However, we do find that the radius where the shape of the total mass distribution begins to differ from that of the stellar mass distribution is dependent on $σ/m$. This transition could potentially be used to set limits on the SIDM cross-section in the MW.

preprint2022arXiv

Snowmass2021 CMB-HD White Paper

CMB-HD is a proposed millimeter-wave survey over half the sky that would be ultra-deep (0.5 uK-arcmin) and have unprecedented resolution (15 arcseconds at 150 GHz). Such a survey would answer many outstanding questions about the fundamental physics of the Universe. Major advances would be 1.) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10 h Mpc^(-1)), which probes dark matter particle properties. It will also allow 2.) measurements of the thermal and kinetic Sunyaev-Zel&#39;dovich effects on small scales to map the gas density and velocity, another probe of cosmic structure. In addition, CMB-HD would allow us to cross critical thresholds: 3.) ruling out or detecting any new, light (< 0.1 eV) particles that were in thermal equilibrium with known particles in the early Universe, 4.) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe, and 5.) ruling out or detecting inflationary magnetic fields. CMB-HD would also provide world-leading constraints on 6.) axion-like particles, 7.) cosmic birefringence, 8.) the sum of the neutrino masses, and 9.) the dark energy equation of state. The CMB-HD survey would be delivered in 7.5 years of observing 20,000 square degrees of sky, using two new 30-meter-class off-axis crossed Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

preprint2022arXiv

Snowmass2021 Cosmic Frontier: Synergies between dark matter searches and multiwavelength/multimessenger astrophysics

This whitepaper focuses on the astrophysical systematics which are encountered in dark matter searches. Oftentimes in indirect and also in direct dark matter searches, astrophysical systematics are a major limiting factor to sensitivity to dark matter. Just as there are many forms of dark matter searches, there are many forms of backgrounds. We attempt to cover the major systematics arising in dark matter searches using photons -- radio and gamma rays -- to cosmic rays, neutrinos and gravitational waves. Examples include astrophysical sources of cosmic messengers and their interactions which can mimic dark matter signatures. In turn, these depend on commensurate studies in understanding the cosmic environment -- gas distributions, magnetic field configurations -- as well as relevant nuclear astrophysics. We also cover the astrophysics governing celestial bodies and galaxies used to probe dark matter, from black holes to dwarf galaxies. Finally, we cover astrophysical backgrounds related to probing the dark matter distribution and kinematics, which impact a wide range of dark matter studies. In the future, the rise of multi-messenger astronomy, and novel analysis methods to exploit it for dark matter, will offer various strategic ways to continue to enhance our understanding of astrophysical backgrounds to deliver improved sensitivity to dark matter.

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.

preprint2021arXiv

HETDEX [OIII] Emitters I: A spectroscopically selected low-redshift population of low-mass, low-metallicity galaxies

We assemble a sample of 17 low metallicity (7.45 < log(O/H)+12 < 8.12) galaxies with z < 0.1 found spectroscopically, without photometric pre-selection, in early data from the Hobby Eberly Telescope Dark Energy Experiment (HETDEX). Star forming galaxies that occupy the lowest mass and metallicity end of the mass-metallicity relation tend to be under sampled in continuum-based surveys as their spectra are typically dominated by emission from newly forming stars. We search for galaxies with high [OIII]$λ$5007 / [OII]$λ$3727, implying highly ionized nebular emission often indicative of low metallicity systems. With the Second Generation Low Resolution Spectrograph on the Hobby Eberly Telescope we acquired follow-up spectra, with higher resolution and broader wavelength coverage, of each low-metallicity candidate in order to confirm the redshift, measure the H$α$ and [NII] line strengths and, in many cases, obtain deeper spectra of the blue lines. We find our galaxies are consistent with the mass-metallicity relation of typical low mass galaxies. However, galaxies in our sample tend to have similar specific star formation rates (sSFRs) as the incredibly rare &#34;blueberry&#34; galaxies found in (Yang et. al. 2017). We illustrate the power of spectroscopic surveys for finding low mass and metallicity galaxies and reveal that we find a sample of galaxies that are a hybrid between the properties of typical dwarf galaxies and the more extreme blueberry galaxies.

preprint2021arXiv

Metallicity Distribution Function of the Eridanus~II Ultra-Faint Dwarf Galaxy from Hubble Space Telescope Narrow-band Imaging

We use deep narrowband Ca H&K ($F395N$) imaging taken with the Hubble Space Telescope (HST) to construct the metallicity distribution function (MDF) of Local Group (LG) ultra-faint dwarf (UFD) galaxy Eridanus II (Eri II). When combined with archival $F475W$ and $F814W$ data, we measure metallicities for 60 resolved red giant branch stars as faint as $m_{F475W}\sim24$ mag, a factor of $\sim4$x more stars than current spectroscopic MDF determinations. We find that Eri II has a mean metallicity of [Fe/H]$=$-2.50$^{+0.07}_{-0.07}$ and a dispersion of $σ_{\mbox{[Fe/H]}}=0.42^{+0.06}_{-0.06}$, which are consistent with spectroscopic MDFs, though more precisely constrained owing to a larger sample. We identify a handful of extremely metal-poor star candidates (EMP; [Fe/H] $< -3$) that are marginally bright enough for spectroscopic follow up. Eri II&#39;s MDF appears well-described by a leaky box chemical evolution model. We also compute an updated orbital history for Eri II using Gaia eDR3 proper motions, and find that it is likely on first infall into the Milky Way. Our findings suggest that Eri II underwent an evolutionary history similar to that of an isolated galaxy. Compared to MDFs for select cosmological simulations of similar mass galaxies, we find that Eri II has a lower fraction of stars with [Fe/H] $< -3$, though such comparisons should currently be treated with caution due to a paucity of simulations, selection effects, and known limitations of Ca H&K for EMPs. This study demonstrates the power of deep HST CaHK imaging for measuring the MDFs of UFDs.

preprint2021arXiv

Out of sight, out of mind? The impact of correlated clustering in substructure lensing

A promising route for revealing the existence of dark matter structures on mass scales smaller than the faintest galaxies is through their effect on strong gravitational lenses. We examine the role of local, lens-proximate clustering in boosting the lensing probability relative to contributions from substructure and unclustered line-of-sight (LOS) halos. Using two cosmological simulations that can resolve halo masses of $M_{\rm halo} \simeq 10^{9}\ M_{\odot}$ (in a simulation box of length $L_{\rm box}{\sim}100\,{\rm Mpc}$) and $10^{7}\ M_{\odot}$ ($L_{\rm box}\sim20\,{\rm Mpc}$), we demonstrate that clustering in the vicinity of the lens host produces a clear enhancement relative to an assumption of unclustered halos that persists to $> 20\,R_{\rm vir}$. This enhancement exceeds estimates that use a two-halo term to account for clustering, particularly within $2-5\,R_{\rm vir}$. We provide an analytic expression for this excess, clustered contribution. We find that local clustering boosts the expected count of $10^9 \ M_\odot$ perturbing halos by ${\sim}35\%$ compared to substructure alone, a result that will significantly enhance expected signals for low-redshift ($z_l \simeq 0.2$) lenses, where substructure contributes substantially compared to LOS halos. We also find that the orientation of the lens with respect to the line of sight (e.g., whether the line of sight passes through the major axis of the lens) can also have a significant effect on the lensing signal, boosting counts by an additional $\sim 50\%$ compared to a random orientations. This could be important if discovered lenses are biased to be oriented along their principal axis.

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

CMB-HD: Astro2020 RFI Response

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

preprint2020arXiv

Evolution of giant molecular clouds across cosmic time

Giant molecular clouds (GMCs) are well-studied in the local Universe, however, exactly how their properties vary during galaxy evolution is poorly understood due to challenging resolution requirements, both observational and computational. We present the first time-dependent analysis of giant molecular clouds in a Milky Way-like galaxy and an LMC-like dwarf galaxy of the FIRE-2 (Feedback In Realistic Environments) simulation suite, which have sufficient resolution to predict the bulk properties of GMCs in cosmological galaxy formation self-consistently. We show explicitly that the majority of star formation outside the galactic center occurs within self-gravitating gas structures that have properties consistent with observed bound GMCs. We find that the typical cloud bulk properties such as mass and surface density do not vary more than a factor of 2 in any systematic way after the first Gyr of cosmic evolution within a given galaxy from its progenitor. While the median properties are constant, the tails of the distributions can briefly undergo drastic changes, which can produce very massive and dense self-gravitating gas clouds. Once the galaxy forms, we identify only two systematic trends in bulk properties over cosmic time: a steady increase in metallicity produced by previous stellar populations and a weak decrease in bulk cloud temperatures. With the exception of metallicity we find no significant differences in cloud properties between the Milky Way-like and dwarf galaxies. These results have important implications for cosmological star and star cluster formation and put especially strong constraints on theories relating the stellar initial mass function to cloud properties.

preprint2020arXiv

Self-consistent proto-globular cluster formation in cosmological simulations of high-redshift galaxies

We report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z>5 galaxies from the FIRE project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over 10^4 Msun pc^-2, where the cloud-scale star formation efficiency is near unity and young stars born in these regions are gravitationally bound at birth. These high-pressure clouds are compressed by feedback-driven winds and/or collisions of smaller clouds/gas streams in highly gas-rich, turbulent environments. The newly formed clusters follow a power-law mass function of dN/dM~M^-2. The cluster formation efficiency is similar across galaxies with stellar masses of ~10^7-10^10 Msun at z>5. The age spread of cluster stars is typically a few Myrs and increases with cluster mass. The metallicity dispersion of cluster members is ~0.08 dex in [Z/H] and does not depend on cluster mass significantly. Our findings support the scenario that present-day old globular clusters (GCs) were formed during relatively normal star formation in high-redshift galaxies. Simulations with a stricter/looser star formation model form a factor of a few more/fewer bound clusters per stellar mass formed, while the shape of the mass function is unchanged. Simulations with a lower local star formation efficiency form more stars in bound clusters. The simulated clusters are larger than observed GCs due to finite resolution. Our simulations are among the first cosmological simulations that form bound clusters self-consistently in a wide range of high-redshift galaxies.

preprint2020arXiv

Simulating the &#34;hidden giant&#34; in cold and self-interacting dark matter models

We perform a series of controlled N-body simulations to study realizations of the recently discovered Antlia 2 galaxy in cold dark matter (CDM) and self-interacting dark matter (SIDM) scenarios. Our simulations contain six benchmark models, where we vary the initial halo concentration and the self-scattering cross section. We adopt well-motivated initial stellar and halo masses, and our fiducial orbit has a small pericenter. After evolving in the Milky Way&#39;s tidal field, the simulated galaxies experience significant mass loss and their stellar distributions expand accordingly. These tidal effects are more prominent if the initial halo concentration is lower and if the self-scattering cross section is larger. Our results show that Antlia 2-like galaxies could be realized in CDM if the halo concentration is low and the stellar distribution is diffuse at the infall time, while these conditions could be relaxed in SIDM. We also find all the simulated galaxies predict approximately the same stellar velocity dispersion after imposing selection criteria for stellar particles. This has important implications for testing dark matter models using tidally disturbed systems.

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.

preprint2020arXiv

The Orbital Histories of Magellanic Satellites Using Gaia DR2 Proper Motions

With the release of Gaia DR2, it is now possible to measure the proper motions (PMs) of the lowest mass, ultra-faint satellite galaxies in the Milky Way&#39;s (MW) halo for the first time. Many of these faint satellites are posited to have been accreted as satellites of the Magellanic Clouds (MCs). Using their 6-dimensional phase space information, we calculate the orbital histories of 13 ultra-faint satellites and five classical dwarf spheroidals in a combined MW+LMC+SMC potential to determine which galaxies are dynamically associated with the MCs. These 18 galaxies are separated into four classes: i.) long-term Magellanic satellites that have been bound to the MCs for at least the last two consecutive orbits around the MCs (Carina 2, Carina 3, Horologium 1, Hydrus 1); ii.) Magellanic satellites that were recently captured by the MCs $<$ 1 Gyr ago (Reticulum 2, Phoenix 2); iii.) MW satellites that have interacted with the MCs (Sculptor 1, Tucana 3, Segue 1); and iv.) MW satellites (Aquarius 2, Canes Venatici 2, Crater 2, Draco 1, Draco 2, Hydra 2, Carina, Fornax, Ursa Minor). Results are reported for a range of MW and LMC masses. Contrary to previous work, we find no dynamical association between Carina, Fornax, and the MCs. Finally, we determine that the addition of the SMC&#39;s gravitational potential affects the longevity of satellites as members of the Magellanic system (long-term versus recently captured), but it does not change the total number of Magellanic satellites.

preprint2020arXiv

The universal acceleration scale from stellar feedback

It has been established for decades that rotation curves deviate from the Newtonian gravity expectation given baryons alone below a characteristic acceleration scale $g_{\dagger}\sim 10^{-8}\,\rm{cm\,s^{-2}}$, a scale promoted to a new fundamental constant in MOND. In recent years, theoretical and observational studies have shown that the star formation efficiency (SFE) of dense gas scales with surface density, SFE $\sim Σ/Σ_{\rm crit}$ with $Σ_{\rm crit} \sim \langle\dot{p}/m_{\ast}\rangle/(π\,G)\sim 1000\,\rm{M_{\odot}\,pc^{-2}}$ (where $\langle \dot{p}/m_{\ast}\rangle$ is the momentum flux output by stellar feedback per unit stellar mass in a young stellar population). We argue that the SFE, more generally, should scale with the local gravitational acceleration, i.e. that SFE $\sim g_{\rm tot}g_\mathrm{crit} \equiv (G\,M_{\rm tot}/R^{2}) / \langle\dot{p}/m_{\ast}\rangle$, where $M_{\rm tot}$ is the total gravitating mass and $g_\mathrm{crit}=\langle\dot{p}/m_{\ast}\rangle = π\,G\,Σ_{\rm crit} \approx 10^{-8}\,\rm{cm\,s^{-2}} \approx g_{\dagger}$. Hence the observed $g_\dagger$ may correspond to the characteristic acceleration scale above which stellar feedback cannot prevent efficient star formation, and baryons will eventually come to dominate. We further show how this may give rise to the observed acceleration scaling $g_{\rm obs}\sim(g_{\rm baryon}\,g_{\dagger})^{1/2}$ (where $g_{\rm baryon}$ is the acceleration due to baryons alone) and flat rotation curves. The derived characteristic acceleration $g_{\dagger}$ can be expressed in terms of fundamental constants (gravitational constant, proton mass, and Thomson cross section): $g_{\dagger}\sim 0.1\,G\,m_{p}/σ_{\rm T}$.

preprint2019arXiv

CMB-HD: An Ultra-Deep, High-Resolution Millimeter-Wave Survey Over Half the Sky

A millimeter-wave survey over half the sky, that spans frequencies in the range of 30 to 350 GHz, and that is both an order of magnitude deeper and of higher-resolution than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. By providing such a deep, high-resolution millimeter-wave survey (about 0.5 uK-arcmin noise and 15 arcsecond resolution at 150 GHz), CMB-HD will enable major advances. It will allow 1) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10/hMpc), which probes dark matter particle properties. It will also allow 2) measurements of the thermal and kinetic Sunyaev-Zel&#39;dovich effects on small scales to map the gas density and gas pressure profiles of halos over a wide field, which probes galaxy evolution and cluster astrophysics. In addition, CMB-HD would allow us to cross critical thresholds in fundamental physics: 3) ruling out or detecting any new, light (< 0.1eV), thermal particles, which could potentially be the dark matter, and 4) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe. Such a survey would also 5) monitor the transient sky by mapping the full observing region every few days, which opens a new window on gamma-ray bursts, novae, fast radio bursts, and variable active galactic nuclei. Moreover, CMB-HD would 6) provide a census of planets, dwarf planets, and asteroids in the outer Solar System, and 7) enable the detection of exo-Oort clouds around other solar systems, shedding light on planet formation. CMB-HD will deliver this survey in 5 years of observing half the sky, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. The telescopes will field about 2.4 million detectors (600,000 pixels) in total.

preprint2019arXiv

Galaxy Formation with BECDM -- II. Cosmic Filaments and First Galaxies

Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like (&#34;WDM&#34;) model often used as a proxy for BECDM. Our simulations confirm that &#34;WDM&#34; is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/&#34;WDM&#34; compared to the CDM case: in BECDM/&#34;WDM&#34; first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in &#34;WDM&#34;, is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM.

preprint2019arXiv

Stars made in outflows may populate the stellar halo of the Milky Way

We study stellar-halo formation using six Milky Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that $5-40\%$ of the outer ($50-300$ kpc) stellar halo in each system consists of $\textit{in-situ}$ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by super-bubble winds, which can be compressed, cool, and form co-moving stars. The majority of these stars remain bound to the halo and fall back with orbital properties similar to the rest of the stellar halo at $z=0$.In the outer halo, outflow stars are more spatially homogeneous, metal rich, and alpha-element-enhanced than the accreted stellar halo. At the solar location, up to $\sim 10 \%$ of our kinematically-identified halo stars were born in outflows; the fraction rises to as high as $\sim 40\%$ for the most metal-rich local halo stars ([Fe/H] $> -0.5$). We conclude that the Milky Way stellar halo could contain local counterparts to stars that are observed to form in molecular outflows in distant galaxies. Searches for such a population may provide a new, near-field approach to constraining feedback and outflow physics. A stellar halo contribution from outflows is a phase-reversal of the classic halo formation scenario of Eggen, Lynden-Bell $\&$ Sandange, who suggested that halo stars formed in rapidly $\textit{infalling}$ gas clouds. Stellar outflows may be observable in direct imaging of external galaxies and could provide a source for metal-rich, extreme velocity stars in the Milky Way.