Search

We introduce and apply a new approach to probe the response of galactic stellar haloes to the interplay between cosmological merger histories and galaxy formation physics. We perform dark matter-only, zoomed simulations of two Milky Way-mass hosts and make targeted, controlled changes to their cosmological histories using the genetic modification technique. Populating each history's stellar halo with a semi-empirical, particle tagging approach then enables a controlled study, with all instances converging to the same large-scale structure, dynamical and stellar mass at z  = 0 as their reference. These related merger scenarios alone generate an extended spread in stellar halo mass fractions (1.5 dex) comparable to the observed population, with the largest scatter achieved by growing late (z ≤ 1) major mergers that spread out existing stars to create massive, in - situ dominated stellar haloes. Increasing a last major merger at z ∼ 2 brings more accreted stars into the inner regions, resulting in smaller scatter in the outskirts which are predominantly built by subsequent minor events. Exploiting the flexibility of our semi-empirical approach, we show that the diversity of stellar halo masses across scenarios is reduced by allowing shallower slopes in the stellar mass–halo mass relation for dwarf galaxies, while it remains conserved when central stars are born with hotter kinematics across cosmic time. The merger-dependent diversity of stellar haloes thus responds distinctly to assumptions in modelling the central and dwarf galaxies respectively, opening exciting prospects to constrain star formation and feedback at different galactic mass-scales with the coming generation of deep, photometric observatories. [ABSTRACT FROM AUTHOR]

We present the empirical dust attenuation (EDA) framework—a flexible prescription for assigning realistic dust attenuation to simulated galaxies based on their physical properties. We use the EDA to forward model synthetic observations for three state-of-the-art large-scale cosmological hydrodynamical simulations: SIMBA, IllustrisTNG, and EAGLE. We then compare the optical and UV color–magnitude relations, (g − r) − Mr and (far-UV −near-UV) − Mr, of the simulations to a Mr < − 20 and UV complete Sloan Digital Sky Survey galaxy sample using likelihood-free inference. Without dust, none of the simulations match observations, as expected. With the EDA, however, we can reproduce the observed color–magnitude with all three simulations. Furthermore, the attenuation curves predicted by our dust prescription are in good agreement with the observed attenuation–slope relations and attenuation curves of star-forming galaxies. However, the EDA does not predict star-forming galaxies with low AV since simulated star-forming galaxies are intrinsically much brighter than observations. Additionally, the EDA provides, for the first time, predictions on the attenuation curves of quiescent galaxies, which are challenging to measure observationally. Simulated quiescent galaxies require shallower attenuation curves with lower amplitude than star-forming galaxies. The EDA, combined with forward modeling, provides an effective approach for shedding light on dust in galaxies and probing hydrodynamical simulations. This work also illustrates a major limitation in comparing galaxy formation models: by adjusting dust attenuation, simulations that predict significantly different galaxy populations can reproduce the same UV and optical observations. [ABSTRACT FROM AUTHOR]

We present the empirical dust attenuation (EDA) frameworkâ€"a flexible prescription for assigning realistic dust attenuation to simulated galaxies based on their physical properties. We use the EDA to forward model synthetic observations for three state-of-the-art large-scale cosmological hydrodynamical simulations: SIMBA, IllustrisTNG, and EAGLE. We then compare the optical and UV colorâ€"magnitude relations, (g âˆ' r) âˆ' M r and (far-UV âˆ'near-UV) âˆ' M r , of the simulations to a M r < âˆ' 20 and UV complete Sloan Digital Sky Survey galaxy sample using likelihood-free inference. Without dust, none of the simulations match observations, as expected. With the EDA, however, we can reproduce the observed colorâ€"magnitude with all three simulations. Furthermore, the attenuation curves predicted by our dust prescription are in good agreement with the observed attenuationâ€"slope relations and attenuation curves of star-forming galaxies. However, the EDA does not predict star-forming galaxies with low A V since simulated star-forming galaxies are intrinsically much brighter than observations. Additionally, the EDA provides, for the first time, predictions on the attenuation curves of quiescent galaxies, which are challenging to measure observationally. Simulated quiescent galaxies require shallower attenuation curves with lower amplitude than star-forming galaxies. The EDA, combined with forward modeling, provides an effective approach for shedding light on dust in galaxies and probing hydrodynamical simulations. This work also illustrates a major limitation in comparing galaxy formation models: by adjusting dust attenuation, simulations that predict significantly different galaxy populations can reproduce the same UV and optical observations. [ABSTRACT FROM AUTHOR]

We compare three major large-scale hydrodynamical galaxy simulations (EAGLE, Illustris-TNG, and SIMBA) by forward modeling simulated galaxies into observational space and computing the fraction of isolated and quiescent low-mass galaxies as a function of stellar mass. Using SDSS as our observational template, we create mock surveys and synthetic spectroscopic and photometric observations of each simulation, adding realistic noise and observational limits. All three simulations show a decrease in the number of quiescent, isolated galaxies in the mass range M* = 109−10 M⊙, in broad agreement with observations. However, even after accounting for observational and selection biases, none of the simulations reproduce the observed absence of quiescent field galaxies below M* = 109 M⊙. We find that the low-mass quiescent populations selected via synthetic observations have consistent quenching timescales, despite an apparent variation in the late-time star formation histories. The effect of increased numerical resolution is not uniform across simulations and cannot fully mitigate the differences between the simulations and the observations. The framework presented here demonstrates a path toward more robust and accurate comparisons between theoretical simulations and galaxy survey observations, while the quenching threshold serves as a sensitive probe of feedback implementations. [ABSTRACT FROM AUTHOR]

The observed breakBRD ("break bulges in red disks") galaxies are a nearby sample of face-on disk galaxies with particularly centrally concentrated star formation: they have red disks but recent star formation in their centers as measured by the Dn4000 spectral index. In this paper, we search for breakBRD analogs in the IllustrisTNG simulation and describe their history and future. We find that a small fraction (∼4% at z = 0; ∼1% at z = 0.5) of galaxies fulfill the breakBRD criteria, in agreement with observations. In comparison with the mass-weighted parent IllustrisTNG sample, these galaxies tend to consist of a higher fraction of satellite and splashback galaxies. However, the central, non-splashback breakBRD galaxies show similar environments, black hole masses, and merger rates, indicating that there is not a single formation trigger for inner star formation and outer quenching. We determine that breakBRD analog galaxies as a whole are in the process of quenching. The breakBRD state, with its highly centrally concentrated star formation, is uncommon in the history of either currently quiescent or star-forming galaxies; however, approximately 10% of 1010 < M*/M⊙ < 1011 quiescent galaxies at z = 0 have experienced SFR concentrations comparable to those of the breakBRDs in their past. Additionally, the breakBRD state is short lived, lasting a few hundred Myr up to ∼2 Gyr. The observed breakBRD galaxies may therefore be a unique sample of outside-in quenching galaxies. [ABSTRACT FROM AUTHOR]

We study the relationship between supermassive black hole (BH) feedback, BH luminosity, and the kinematics of stars and gas for galaxies in IllustrisTNG. We use galaxies with mock MaNGA observations to identify kinematic misalignment at z  = 0 (difference in rotation of stars and gas), for which we follow the evolutionary history of BH activity and gas properties over the last 8 Gyr. Misaligned low-mass galaxies (⁠|$\mathrm{\mathit{ M}_{stel} \lt 10^{10.2}\, \mathrm{M}_{\odot }}$|⁠) typically have boosted BH luminosity and BH growth, and have had more energy injected by BHs into the gas over the last 8 Gyr in comparison to low-mass aligned galaxies. These properties likely lead to outflows and gas removal, in agreement with active low mass galaxies in observations. Splitting on BH luminosity at z  = 0 produces consistent distributions of kinematic misalignment at z  = 0; however, splitting on the maximum BH luminosity over the last 8 Gyr produces statistically significant different distributions. While instantaneous correlation at z  = 0 is difficult due to misalignment persisting on longer time-scales, the relationship between BH activity and misalignment is clear. High-mass quenched galaxies (⁠|$\mathrm{\mathit{ M}_{stel} \gt 10^{10.2}\, \mathrm{M}_{\odot }}$|⁠) with misalignment typically have similar BH luminosities, show lower gas fractions, and have typically lower gas phase metallicity over the last 8 Gyr in comparison to the high mass aligned. [ABSTRACT FROM AUTHOR]

We present a new suite of cosmological zoom-in hydrodynamical (⁠|$\approx 20\, \mathrm{pc}$| spatial resolution) simulations of Milky-Way mass galaxies to study how a varying mass ratio for a Gaia-Sausage-Enceladus (GSE) progenitor impacts the z  = 0 chemodynamics of halo stars. Using the genetic modification approach, we create five cosmological histories for a Milky-Way-mass dark matter halo (⁠|$M_{200}\approx 10^{12} \, \mbox{M}_\mathrm{\odot }$|⁠), incrementally increasing the stellar mass ratio of a z ≈ 2 merger from 1:25 to 1:2, while fixing the galaxy's final dynamical, stellar mass, and large-scale environment. We find markedly different morphologies at z  = 0 following this change in early history, with a growing merger resulting in increasingly compact and bulge-dominated galaxies. Despite this structural diversity, all galaxies show a radially biased population of inner halo stars like the Milky-Way's GSE which, surprisingly, has a similar magnitude, age, |$\rm [Fe/H]$|⁠ , and |$\rm [\alpha /Fe]$| distribution whether the z ≈ 2 merger is more minor or major. This arises because a smaller ex-situ population at z ≈ 2 is compensated by a larger population formed in an earlier merger-driven starburst whose contribution to the GES can grow dynamically over time, and with both populations strongly overlapping in the |$\rm [Fe/H]-\rm [\alpha /Fe]$| plane. Our study demonstrates that multiple high-redshift histories can lead to similar z  = 0 chemodynamical features in the halo, highlighting the need for additional constraints to distinguish them, and the importance of considering the full spectrum of progenitors when interpreting z  = 0 data to reconstruct our Galaxy's past. [ABSTRACT FROM AUTHOR]

Stellar streams from globular clusters (GCs) offer constraints on the nature of dark matter and have been used to explore the dark matter halo structure and substructure of our Galaxy. Detection of GC streams in other galaxies would broaden this endeavor to a cosmological context, yet no such streams have been detected to date. To enable such exploration, we develop the Hough Stream Spotter (HSS), and apply it to the Pan-Andromeda Archaeological Survey (PAndAS) photometric data of resolved stars in M31's stellar halo. We first demonstrate that our code can re-discover known dwarf streams in M31. We then use the HSS to blindly identify 27 linear GC stream-like structures in the PAndAS data. For each HSS GC stream candidate, we investigate the morphologies of the streams and the colors and magnitudes of all stars in the candidate streams. We find that the five most significant detections show a stronger signal along the red giant branch in color–magnitude diagrams than spurious non-stream detections. Lastly, we demonstrate that the HSS will easily detect globular cluster streams in future Nancy Grace Roman Space Telescope data of nearby galaxies. This has the potential to open up a new discovery space for GC stream studies, GC stream gap searches, and for GC stream-based constraints on the nature of dark matter. [ABSTRACT FROM AUTHOR]

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

Understanding the variability of galaxy star formation histories (SFHs) across a range of time-scales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at z  = 0 from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIRE-2, g14, and Marvel/Justice League), semi-analytic models (Santa Cruz SAM) and empirical models (UniverseMachine), and quantify the variability of these SFHs on different time-scales using the power spectral density (PSD) formalism. We find that the PSDs are well described by broken power laws, and variability on long time-scales (≳1 Gyr) accounts for most of the power in galaxy SFHs. Most hydrodynamical models show increased variability on shorter time-scales (≲300 Myr) with decreasing stellar mass. Quenching can induce ∼0.4−1 dex of additional power on time-scales >1 Gyr. The dark matter accretion histories of galaxies have remarkably self-similar PSDs and are coherent with the in situ star formation on time-scales >3 Gyr. There is considerable diversity among the different models in their (i) power due to star formation rate variability at a given time-scale, (ii) amount of correlation with adjacent time-scales (PSD slope), (iii) evolution of median PSDs with stellar mass, and (iv) presence and locations of breaks in the PSDs. The PSD framework is a useful space to study the SFHs of galaxies since model predictions vary widely. Observational constraints in this space will help constrain the relative strengths of the physical processes responsible for this variability. [ABSTRACT FROM AUTHOR]

Aims. Several kinematic and chemical substructures have been recently found amongst Milky Way halo stars with retrograde motions. It is currently unclear how these various structures are related to each other. This Letter aims to shed light on this issue. Methods. We explore the retrograde halo with an augmented version of the Gaia DR2 RVS sample, extended with data from three large spectroscopic surveys, namely RAVE, APOGEE, and LAMOST. In this dataset, we identify several structures using the HDBSCAN clustering algorithm. We discuss their properties and possible links using all the available chemical and dynamical information. Results. In concordance with previous work, we find that stars with [Fe/H] < −1 have more retrograde motions than those with [Fe/H] > −1. The retrograde halo contains a mixture of debris from objects like Gaia-Enceladus, Sequoia, and even the chemically defined thick disc. We find that the Sequoia has a smaller range in orbital energies than previously suggested and is confined to high energy. Sequoia could be a small galaxy in itself, but since it overlaps both in integrals-of-motion space and chemical abundance space with the less bound debris of Gaia-Enceladus, its nature cannot yet be fully settled. In the low-energy part of the halo, we find evidence for at least one more distinct structure: Thamnos. Stars in Thamnos are on low-inclination, mildly eccentric retrograde orbits, moving at vϕ ≈ −150 km s−1, and are chemically distinct from the other structures. Conclusions. Even with the excellent Gaia DR2 data, piecing together all the fragments found in the retrograde halo remains challenging. At this point, we are very much in need of large datasets with high-quality high-resolution spectra and tailored high-resolution hydrodynamical simulations of galaxy mergers. [ABSTRACT FROM AUTHOR]