Search for: All records

Abstract ΛCDM cosmology predicts the hierarchical formation of galaxies, which build up mass by merger events and accreting smaller systems. The stellar halo of the Milky Way (MW) has proven to be useful a tool for tracing this accretion history. However, most of this work has focused on the outer halo where dynamical times are large and the dynamical properties of accreted systems are preserved. In this work, we investigate the inner galaxy regime, where dynamical times are relatively small and systems are generally completely phase mixed. Using the FIRE-2 and Auriga cosmological zoom-in simulation suites of MW-mass galaxies, we find the stellar density profiles along the minor axis (perpendicular to the galactic disk) within the Navarro–Frenk–White scale radii (R ≈ 15 kpc) are best described as an exponential disk with scale height < 0.3 kpc and a power-law component with slopeα ≈ −4. The stellar density amplitude and slope for the power-law component are not significantly correlated with metrics of the galaxy’s accretion history. Instead, we find the stellar profiles strongly correlate with the dark matter profile. Across simulation suites, the galaxies studied in this work have a stellar-to-dark-matter mass ratio that decreases as 1/r²along the minor axis. 

Abstract We utilize the cosmological volume simulation FIREbox to investigate how a galaxy’s environment influences its size and dark matter content. Our study focuses on approximately 1200 galaxies (886 central and 332 satellite halos) in the low-mass regime, with stellar masses between 10⁶and 10⁹M_⊙. We analyze the size–mass relation (r₅₀–M_⋆), the inner dark matter mass–stellar mass ( $M_{\mathrm{DM}}^{50}$–M_⋆) relation, and the halo mass–stellar mass (M_halo–M_⋆) relation. At fixed stellar mass, we find that galaxies experiencing stronger tidal influences, indicated by higher Perturbation Indices (PI > 1) are generally larger and have lower halo masses relative to their counterparts with lower Perturbation Indices (PI < 1). Applying a Random Forest regression model, we show that both the environment (PI) and halo mass (M_halo) are significant predictors of a galaxy’s relative size and dark matter content. Notably, becauseM_halois also strongly affected by the environment, our findings indicate that environmental conditions not only influence galactic sizes and relative inner dark matter content directly, but also indirectly, through their impact on halo mass. Our results highlight a critical interplay between environmental factors and halo mass in shaping galaxy properties, affirming the environment as a fundamental driver in galaxy formation and evolution. 

Context.The detection of supermassive black holes (SMBHs) in high-redshift luminous quasars may require a phase of rapid accretion, and as a precondition, substantial gas influx toward seed black holes (BHs) from kiloparsec or parsec scales. Our previous research demonstrated the plausibility of such gas supply for BH seeds within star-forming giant molecular clouds (GMCs) with high surface density (∼10⁴ M_⊙ pc⁻²), facilitating “hyper-Eddington” accretion via efficient feeding by dense clumps, which are driven by turbulence and stellar feedback. Aims.This article presents an investigation of the impacts of feedback from accreting BHs on this process, including radiation, mechanical jets, and highly relativistic cosmic rays. Methods.We ran a suite of numerical simulations to explore diverse parameter spaces of BH feedback, including the subgrid accretion model, feedback energy efficiency, mass loading factor, and initial metallicity. Results.Using radiative feedback models inferred from the slim disk, we find that hyper-Eddington accretion is still achievable, yielding BH bolometric luminosities of as high as 10⁴¹ − 10⁴⁴ erg/s, depending on the GMC properties and specific feedback model assumed. We find that the maximum possible mass growth of seed BHs (ΔM^max_BH) is regulated by the momentum-deposition rate from BH feedback,ṗ_feedback/(Ṁ_BHc), which leads to an analytic scaling that agrees well with simulations. This scenario predicts the rapid formation of ∼10⁴M_⊙intermediate-massive BHs (IMBHs) from stellar-mass BHs within ∼1 Myr. Furthermore, we examine the impacts of subgrid accretion models and how BH feedback may influence star formation within these cloud complexes. 

Abstract The physical mechanisms responsible for bar formation and destruction in galaxies remain a subject of debate. While we have gained valuable insight into how bars form and evolve from isolated idealized simulations, in the cosmological domain, galactic bars evolve in complex environments, with mergers and gas accretion events occurring in the presence of the turbulent interstellar medium with multiple star formation episodes, in addition to coupling with their host galaxies’ dark matter halos. We investigate the bar formation in 13 Milky Way–mass galaxies from the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulations. 8 of the 13 simulated galaxies form bars at some point during their history: three from tidal interactions and five from internal evolution of the disk. The bars in FIRE-2 are generally shorter than the corotation radius (mean bar radius ∼1.53 kpc), have a wide range of pattern speeds (36–97 km s⁻¹kpc⁻¹), and live for a wide range of dynamical times (2–160 bar rotations). We find that the bar formation in FIRE-2 galaxies is influenced by satellite interactions and the stellar-to-dark-matter mass ratio in the inner galaxy, but neither is a sufficient condition for bar formation. Bar formation is more likely to occur, with the bars formed being stronger and longer-lived, if the disks are kinematically cold; galaxies with high central gas fractions and/or vigorous star formation, on the other hand, tend to form weaker bars. In the case of the FIRE-2 galaxies, these properties combine to produce ellipsoidal bars with strengthsA₂/A₀∼ 0.1–0.2. 

Abstract Stellar feedback influences the star formation rate (SFR) and the interstellar medium of galaxies in ways that are difficult to quantify numerically, because feedback is an essential ingredient of realistic simulations. To overcome this, we conduct a feedback-halting experiment starting with a Milky Way–mass galaxy in the second-generation Feedback In Realistic Environments (FIRE-2) simulation framework. By terminating feedback, and comparing to a simulation in which feedback is maintained, we monitor how the runs diverge. We find that without feedback, the interstellar turbulent velocities decay. There is a marked increase of dense material, while the SFR increases by over an order of magnitude. Importantly, this SFR boost is a factor of ∼15–20 larger than is accounted for by the increased freefall rate caused by higher densities. This implies that feedback moderates the star formation efficiency per freefall time more directly than simply through the density distribution. To probe changes at the scale of giant molecular clouds (GMCs), we identify GMCs using density and virial parameter thresholds, tracking clouds as the galaxy evolves. Halting feedback stimulates rapid changes, including a proliferation of new bound clouds, a decrease of turbulent support in loosely bound clouds, an overall increase in cloud densities, and a surge of internal star formation. Computing the cloud-integrated SFR using several theories of turbulence regulation, we show that these theories underpredict the surge in SFR by at least a factor of 3. We conclude that galactic star formation is essentially feedback regulated on scales that include GMCs, and that stellar feedback affects GMCs in multiple ways. 

Abstract Feedback from supermassive black holes is believed to be a critical driver of the observed color bimodality of galaxies above the Milky Way mass scale. Active galactic nuclei (AGN) feedback has been modeled in many galaxy formation simulations, but most implementations have involved simplified prescriptions or a coarse-grained interstellar medium (ISM). We present the first set of Feedback In Realistic Environments (FIRE)-3 cosmological zoom-in simulations with AGN feedback evolved toz∼ 0, examining the impact of AGN feedback on a set of galaxies with halos in the mass range 10¹²–10¹³M_⊙. These simulations combine detailed stellar and ISM physics with multichannel AGN feedback including radiative feedback, mechanical outflows, and, in some simulations, cosmic rays (CRs). We find that massive (>L*) galaxies in these simulations can match local scaling relations including the stellar mass–halo mass relation and theM_BH–σrelation; in the stronger model with CRs, they also match the size–mass relation and the Faber–Jackson relation. Many of the massive galaxies in the simulations with AGN feedback have quenched star formation and elliptical morphologies, in qualitative agreement with observations. In contrast, simulations at the massive end without AGN feedback produce galaxies that are too massive and form stars too rapidly, are order-of-magnitude too compact, and have velocity dispersions well above Faber–Jackson. Despite these successes, the AGN models analyzed do not produce uniformly realistic galaxies when the feedback parameters are held constant: While the stronger model produces the most realistic massive galaxies, it tends to overquench the lower-mass galaxies. This indicates that further refinements of the AGN modeling are needed. 

Abstract We analyze the first 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 modeled as a constant fractional energy dissipation (f_diss= 0.75). In this paper, we examine the properties of dwarf galaxies withM_*∼ 10⁵–10⁹M_⊙in both isolation and within Milky Way–mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of disk formation in dSIDM with (σ/m) ≤ 1 cm²g⁻¹. On the contrary, models with (σ/m) = 10 cm²g⁻¹produce puffier stellar distributions that are in tension with the observed size–mass relation. In addition, owing to the steeper central density profiles, the subkiloparsec circular velocities of isolated dwarfs when (σ/m) ≥ 0.1 cm²g⁻¹are enhanced by about a factor of 2, which are still consistent with the kinematic measurements of Local Group dwarfs but in tension with the Hirotation curves of more massive field dwarfs. Meanwhile, for satellites of Milky Way–mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help promote the structural diversity of dwarf satellites. 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 (σ/m) ≳ 0.1 cm²g⁻¹,f_diss= 0.75 are in tension in massive dwarfs (M_halo∼ 10¹¹M_⊙) due to circular velocity constraints. However, models with lower effective cross sections (at this halo mass/velocity scale) are still viable and can produce nontrivial observable signatures. 

Abstract Planet formation is expected to be severely limited in disks of low metallicity, owing to both the small solid mass reservoir and the low-opacity accelerating the disk gas dissipation. While previous studies have found a weak correlation between the occurrence rates of small planets (≲4R_⊕) and stellar metallicity, so far no studies have probed below the metallicity limit beyond which planet formation is predicted to be suppressed. Here, we constructed a large catalog of ∼110,000 metal-poor stars observed by the TESS mission with spectroscopically derived metallicities, and systematically probed planet formation within the metal-poor regime ([Fe/H] ≤−0.5) for the first time. Extrapolating known higher-metallicity trends for small, short-period planets predicts the discovery of ∼68 super-Earths around these stars (∼85,000 stars) after accounting for survey completeness; however, we detect none. As a result, we have placed the most stringent upper limit on super-Earth occurrence rates around metal-poor stars (−0.75 < [Fe/H] ≤ −0.5) to date, ≤ 1.67%, a statistically significant (p-value = 0.000685) deviation from the prediction of metallicity trends derived with Kepler and K2. We find a clear host star metallicity cliff for super-Earths that could indicate the threshold below which planets are unable to grow beyond an Earth-mass at short orbital periods. This finding provides a crucial input to planet-formation theories, and has implications for the small planet inventory of the Galaxy and the galactic epoch at which the formation of small planets started. 

ABSTRACT We present predictions for the high-redshift halo–galaxy–supermassive black hole (SMBH) connection from the Trinity model. Matching a comprehensive compilation of galaxy (0 ≤ z ≤ 13) and SMBH data sets (0 ≤ z ≤ 6.5), Trinity finds: (1) The number of SMBHs with M• > 109 M⊙ in the observable Universe increases by five orders of magnitude from z ∼ 10 to z ∼ 2, and by another factor of ∼3 from z ∼ 2 to z = 0; (2) The M• > 109 and 1010 M⊙ SMBHs at z ∼ 6 live in haloes with ∼(2 − 3) and (3 − 5) × 1012 M⊙; (3) the newly discovered JWST AGN candidates at 7 ≲ z ≲ 11 are overmassive compared to the intrinsic SMBH mass–galaxy mass relation from Trinity, but they are still broadly consistent with Trinity predictions for flux limited AGN samples with Lauer bias. This bias favours the detection for overmassive SMBHs due to higher luminosities at a fixed Eddington ratio. However UHZ1’s M•/M* ratio is still some 1 dex higher than Trinity AGNs, indicating a discrepancy; (4) Trinity underpredicts the number densities of GN-z11 and CEERS_1019 analogues. But given the strong constraints from existing data in Trinity, the extra constraint from GN-z11 and CEERS_1019 does not significantly change trinity model results. (5) z = 6–10 quasar luminosity functions will reduce uncertainties in the trinity prediction of the z = 6–10 SMBH mass–galaxy mass relation by up to ∼0.5 dex. These luminosity functions will be available with future telescopes, such as Roman and Euclid. 

Abstract We make an in-depth analysis of different active galactic nuclei (AGN) jet models’ signatures, inducing quiescence in galaxies with a halo mass of 10¹²M_⊙. Three jet models, including cosmic-ray-dominant, hot thermal, and precessing kinetic jets, are studied at two energy flux levels each, compared to a jet-free, stellar feedback-only simulation. Each of our simulations is idealized isolated galaxy simulations with AGN jet powers that are constant in time and generated using GIZMO and with FIRE stellar feedback. We examine the distribution of Mgii, Ovi, and Oviiiions, alongside gas temperature and density profiles. Low-energy ions, like Mgii, concentrate in the interstellar medium (ISM), while higher energy ions, e.g., Oviii, prevail at the AGN jet cocoon’s edge. High-energy flux jets display an isotropic ion distribution with lower overall density. High-energy thermal or cosmic-ray jets pressurize at smaller radii, significantly suppressing core density. The cosmic-ray jet provides extra pressure support, extending cool and warm gas distribution. A break in the ion-to-mass ratio slope in Oviand Oviiiis demonstrated in the ISM-to-circumgalactic medium (CGM) transition (between 10 and 30 kpc), growing smoothly toward the CGM at greater distances.