Search

Numerous stars exhibit surprisingly large variations in their refractory element abundances, often interpreted as signatures of planetary ingestion events. In this study, we propose that differences in the dust-to-gas ratio near stars during their formation can produce similar observational signals. We investigate this hypothesis using a suite of radiation-dust-magnetohydrodynamic STARFORGE simulations of star formation. Our results show that the distribution of refractory abundance variations ($\rm \Delta [X/H]$) has extended tails, with about 10-30% of all stars displaying variations around $\sim$0.1 dex. These variations are comparable to the accretion of $2-5 \rm M_\oplus$ of planetary material into the convective zones of Sun-like stars. The width of the distributions increases with the incorporation of more detailed dust physics, such as radiation pressure and back-reaction forces, as well as with larger dust grain sizes and finer resolutions. Furthermore, our simulations reveal no correlation between $\rm \Delta [X/H]$ and stellar separations, suggesting that dust-to-gas fluctuations likely occur on scales smaller than those of wide binaries. These findings highlight the importance of considering dust dynamics as a potential source of the observed chemical enrichment in stars., Comment: 10 pages, 3 figures. Submitted to ApJL

We present the first suite of cosmological hydrodynamical zoom-in simulations of isolated dwarf galaxies for a dark sector that consists of Cold Dark Matter and a strongly-dissipative sub-component. The simulations are implemented in GIZMO and include standard baryons following the FIRE-2 galaxy formation physics model. The dissipative dark matter is modeled as Atomic Dark Matter (aDM), which forms a dark hydrogen gas that cools in direct analogy to the Standard Model. Our suite includes seven different simulations of $\sim 10^{10} M_{\odot}$ systems that vary over the aDM microphysics and the dwarf's evolutionary history. We identify a region of aDM parameter space where the cooling rate is aggressive and the resulting halo density profile is universal. In this regime, the aDM gas cools rapidly at high redshifts and only a small fraction survives in the form of a central dark gas disk; the majority collapses centrally into collisionless dark "clumps", which are clusters of sub-resolution dark compact objects. These dark clumps rapidly equilibrate in the inner galaxy, resulting in an approximately isothermal distribution that can be modeled with a simple fitting function. Even when only a small fraction ($\sim 5\%$) of the total dark matter is strongly dissipative, the central densities of classical dwarf galaxies can be enhanced by over an order of magnitude, providing a sharp prediction for observations., Comment: 17 pages, 8 pages of appendices

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 ($\leq$4R$_\oplus$) 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] $\leq$ -0.5) for the first time. Extrapolating known higher-metallicity trends for small, short-period planets predicts the discovery of ~68 superEarths 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] $\leq$ -0.5) to date, $\leq$ 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., Comment: Accepted to AJ

The dust grain size distribution (GSD) likely varies significantly across different star-forming environments in the Universe, but the overall impact of this variation on star formation remains unclear. This ambiguity arises because the GSD interacts non-linearly with processes like heating/cooling, radiation, and chemistry, which have competing effects and different environmental dependencies. In this study, we investigate the effects of GSD variation on the thermochemistry and evolution of giant molecular clouds (GMCs). To achieve this, we conducted radiation-dust-magnetohydrodynamic simulations spanning a range of cloud masses and grain sizes, which explicitly incorporate the dynamics of dust grains within the full-physics framework of the STARFORGE project. We find that differences in grain size significantly alter the thermochemistry of GMCs. Specifically, we show that the leading-order effect is that larger grains, under fixed dust mass and dust-to-gas ratio conditions, result in lower dust opacities. This reduced opacity permits ISRF photons to penetrate more deeply and allows internal radiation field photons to permeate more extensively into the cloud, resulting in rapid gas heating and the inhibition of star formation. We find that star formation efficiency is highly sensitive to grain size, with an order of magnitude reduction in efficiency when grain size increases from 0.1 $\rm\mu m$ to 10 $\rm\mu m$. Additionally, we note that warmer gas suppresses the formation of low-mass stars. Moreover, as a consequence of the decreased opacities, we observe a greater proportion of gas residing in diffuse ionized structures., Comment: 13 pages, 5 figures, submitted to ApJ

Recent simulations have demonstrated the formation of 'flux-frozen' and hyper-magnetized disks, qualitatively distinct from both classical $\alpha$ disks and magnetically-arrested disks, as a natural consequence of fueling gas to supermassive black holes in galactic nuclei. We previously showed that the dynamical structure of said disks can be approximated by simple analytic similarity models. Here we study the thermal properties of these models over a wide range of physical scales and accretion rates. We show there are several characteristic zones: a dusty 'torus'-like region, a multi-phase neutral and then multi-phase ionized, broad line-emitting region interior to the sublimation radius, before finally a transition to a thermal accretion disk with a warm Comptonizing layer. The disks are strongly-flared with large scale heights, and reprocess and/or scatter an order-one fraction of the central disk emission. As a result, this simple accretion disk model predicts phenomena including the existence of a dusty torus and its covering factor, geometry, clumpiness, and dust temperatures; a broad-line-region (BLR) with its characteristic sizes and luminosities and ionization properties; extended scattering/reprocessing surfaces producing cooler disk continuum and apparently large observed disk sizes; and existence of warm Comptonizing layers and hard coronal gas. Remarkably, these properties emerge without our having to introduce new components or parameters: they are all part of the accretion flow if the disks are in the hyper-magnetized limit., Comment: Submitted to the Open Journal of Astrophysics. Comments welcome. 36 pages, 14 figures

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 FIRE-2 simulation framework. Terminating feedback, and comparing to a simulation in which feedback is maintained, we monitor how the runs diverge. We find that without feedback, 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 $\sim$15-20 larger than is accounted for by the increased free fall rate caused by higher densities. This implies that feedback moderates the star formation efficiency per free-fall 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 three. We conclude that galactic star formation is essentially feedback-regulated on scales that include GMCs, and that stellar feedback affects GMCs in multiple ways., Comment: 20 pages, 10 figures, under ongoing review at ApJ

Stars form within dense cores composed of both gas and dust within molecular clouds. However, despite the crucial role that dust plays in the star formation process, its dynamics is frequently overlooked, with the common assumption being a constant, spatially uniform dust-to-gas ratio and grain size spectrum. In this study, we introduce a set of radiation-dust-magnetohydrodynamic simulations of star forming molecular clouds from the {\small STARFORGE} project. These simulations expand upon the earlier radiation MHD models, which included cooling, individual star formation, and feedback. Notably, they explicitly address the dynamics of dust grains, considering radiation, drag, and Lorentz forces acting on a diverse size spectrum of live dust grains. We find that interactions between radiation and dust significantly influence the properties of gas surrounding and accreting onto massive stars. Specifically, we find that once stars exceed a certain mass threshold ($\sim 2 M_{\odot}$), their emitted radiation can evacuate dust grains from their vicinity, giving rise to a dust-suppressed zone of size $\sim 100$ AU. Commencing during the early accretion stages and preceding the Main-sequence phase, this process results in a mass-dependent depletion in the accreted dust-to-gas (ADG) mass ratio within both the circumstellar disc and the star. We predict massive stars ($\gtrsim 10 M_{\odot}$) would exhibit ADG ratios that are approximately one order of magnitude lower than that of their parent clouds. Consequently, stars, their discs, and circumstellar environments would display notable deviations in the abundances of elements commonly associated with dust grains, such as carbon and oxygen., Comment: 14 pages, 9 figures, submitted to ApJ

We present an investigation of clustered stellar feedback in the form of superbubbles identified within eleven galaxies from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, at both cosmic noon (1 < z < 3) and in the local Universe. We study the spatially-resolved multiphase outflows that these supernovae drive, comparing our findings with recent theory and observations. These simulations consist of five LMC-mass galaxies and six Milky Way-mass progenitors (with a minimum baryonic particle mass of $m_{b.min} = 7100 M_{\odot}$), for which we calculate the local mass and energy loading factors on 750~pc scales from the identified outflows. We also characterize the multiphase morphology and properties of the identified superbubbles, including the `shell' of cool ($T<10^5$ K) gas and break out of energetic hot ($T>10^5$ K) gas when the shell bursts. For all galaxies, the outflow mass, momentum, and energy fluxes appear to reach their peak during the identified superbubbles, and we investigate the effects on the interstellar medium (ISM), circumgalactic medium (CGM), and subsequent star formation rates. We find that these simulations, regardless of redshift, have mass-loading factors and momentum fluxes in the cool gas that largely agree with recent observations. Lastly, we also investigate how methodological choices in measuring outflows can affect loading factors for galactic winds.

The origins and mergers of supermassive black holes (BHs) remain a mystery. We describe a scenario from a novel multi-physics simulation featuring rapid ($\lesssim 1\,$Myr) hyper-Eddington gas capture by a $\sim 1000\,{\rm M}_{\odot}$ ``seed'' BH up to supermassive ($\gtrsim 10^{6}\,M_{\odot}$) masses, in a massive, dense molecular cloud complex typical of high-redshift starbursts. Due to the high cloud density, stellar feedback is inefficient and most of the gas turns into stars in star clusters which rapidly merge hierarchically, creating deep potential wells. Relatively low-mass BH seeds at random positions can be ``captured'' by merging sub-clusters and migrate to the center in $\sim1$ free-fall time (vastly faster than dynamical friction). This also efficiently produces a paired BH binary with $\sim 0.1$\,pc separation. The centrally-concentrated stellar density profile (akin to a ``proto-bulge'') allows the cluster as a whole to capture and retain gas and build up a large (pc-scale) circum-binary accretion disk with gas coherently funnelled to the central BH (even when the BH radius of influence is small). The disk is ``hyper-magnetized'' and ``flux-frozen'': dominated by a toroidal magnetic field with plasma $\beta \sim 10^{-3}$, with the fields amplified by flux-freezing. This drives hyper-Eddington inflow rates $\gtrsim 1\,\rm M_\odot yr^{-1}$, which also drive the two BHs to nearly-equal masses. The late-stage system appears remarkably similar to recently-observed high-redshift ``little red dots.'' This scenario can provide an explanation for rapid SMBH formation, growth and mergers in high-redshift galaxies., Comment: 9 pages, 5 figures, final version matching the publication

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 kilo-parsec 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 ($\sim 10^4\,{\rm {\rm M_\odot}\, pc}^{-2}$), facilitating ``hyper-Eddington'' accretion via efficient feeding by dense clumps which are driven by turbulence and stellar feedback. This article investigates the impacts of feedback from accreting BHs on this process, including radiation, mechanical jets, and highly relativistic cosmic rays. We run a suite of numerical simulations to explore diverse parameter spaces of BH feedback, including the sub-grid accretion model, feedback energy efficiency, mass loading factor, and initial metallicity. Utilizing radiative feedback models inferred from the slim disk, we find that hyper-Eddington accretion is still achievable, yielding BH bolometric luminosities as high as $10^{41}$ -- $10^{44}\,\rm erg/s$, depending on the GMC properties and specific feedback model assumed. We find the maximum possible mass growth of seed BHs ($\Delta M_{\rm BH}^{\rm max}$) is regulated by the momentum deposition rate from BH feedback, $\dot{p}_{\rm feedback}/(\dot{M}_{\rm BH} c)$, which leads to an analytic scaling that agrees well with simulations. This scenario predicts the rapid formation of $\sim 10^4\,\rm M_\odot$ intermediate-massive BHs (IMBHs) from stellar-mass BHs within $\sim \rm Myr$. Furthermore, we examine the impacts of sub-grid accretion models and how BH feedback may influence star formation within these cloud complexes., Comment: Comments welcome

The magnetohydrodynamics (MHD) equations plus 'non-ideal' (Ohmic, Hall, ambipolar) resistivities are widely used to model weakly-ionized astrophysical systems. We show that if gradients in the magnetic field become too steep, the implied charge drift speeds become much faster than microphysical signal speeds, invalidating the assumptions used to derive both the resistivities and MHD equations themselves. Generically this situation will excite microscale instabilities that suppress the drift and current. We show this could be relevant at low ionization fractions especially if Hall terms appear significant, external forces induce supersonic motions, or dust grains become a dominant charge carrier. Considering well-established treatments of super-thermal drifts in laboratory, terrestrial, and Solar plasmas as well as conduction and viscosity models, we generalize a simple prescription to rectify these issues, where the resistivities are multiplied by a correction factor that depends only on already-known macroscopic quantities. This is generalized for multi-species and weakly-ionized systems, and leaves the equations unchanged in the drift limits for which they are derived, but restores physical behavior (driving the system back towards slow drift by diffusing away small-scale gradients in the magnetic field) if the limits are violated. This has important consequences: restoring intuitive behaviors such as the system becoming hydrodynamic in the limit of zero ionization; suppressing magnetic structure on scales below a critical length which can comparable to circumstellar disk sizes; limiting the maximum magnetic amplification; and suppressing the effects of the Hall term in particular. This likely implies that the Hall term does not become dynamically important under most conditions of interest in these systems., Comment: 16 pages, 4 figures. Submitted to the Open Journal of Astrophysics. Comments welcome

Fueling star formation in large, discy galaxies requires a continuous supply of gas accreting into star-forming regions. Previously, we characterized this accretion in 4 Milky Way mass galaxies ($M_{\rm halo}\sim10^{12}M_{\odot}$) in the FIRE-2 cosmological zoom-in simulations. At $z\sim0$, we found that gas within the inner circumgalactic medium (iCGM) approaches the disc with comparable angular momentum (AM) to the disc edge, joining in the outer half of the gaseous disc. Within the disc, gas moves inward at velocities of $\sim$1-5~km~s$^{-1}$ while fully rotationally supported. In this study, we analyze the torques that drive these flows. In all cases studied, we find that the torques in discs enable gas accreted near the disc edge to transport inwards and fuel star formation in the central few kpc. The primary sources of torque come from gravity, hydrodynamical forces, and the sub-grid $P dV$ work done by supernova (SNe) remnants interacting with gas on $\lesssim$10 pc scales. These SNe remnant interactions induce negative torques within the inner disc and positive torques in the outer disc. The gas-gas gravitational, hydro, and "feedback" torques transfer AM outward to where accreting gas joins the disc, playing an important role in driving inflows and regulating disc structure. Gravitational torques from stars and dark matter provide an AM sink within the innermost regions of the disc and iCGM, respectively. Feedback torques are dominant within the disc, while gravitational and hydrodynamical torques have similar significance depending on the system/region. Torques from viscous shearing, magnetic forces, stellar winds, and radiative transfer are less significant., Comment: 20 pages, 18 figures. Accepted for publication in MNRAS. Added 1 figure to Appendix B, modified text

In galaxy formation simulations, it is increasingly common to represent supernovae (SNe) at finite resolution (when the Sedov-Taylor phase is unresolved) via hybrid energy-momentum coupling with some 'terminal momentum' $p_{\rm term}$ (depending weakly on ambient density and metallicity) that represents unresolved work from an energy-conserving phase. Numerical implementations can ensure momentum and energy conservation of such methods, but these require that couplings depend on the surrounding gas velocity field (radial velocity $\langle v_{r} \rangle$). This raises the question of whether $p_{\rm term}$ should also be velocity-dependent, which we explore analytically and in simulations. We show that for simple spherical models, the dependence of $p_{\rm term}$ on $\langle v_{r} \rangle$ introduces negligible corrections beyond those already imposed by energy conservation if $\langle v_{r} \rangle \ge 0$. However, for SNe in some net converging flow ($\langle v_{r} \rangle<0$), naively coupling the total momentum when a blastwave reaches the standard cooling/snowplow phase (or some effective cooling time/velocity/temperature criterion) leads to enormous $p_{\rm term}$ and potentially pathological behaviors. We propose an alternative $\Delta$-Momentum formulation which represents the differential SNe effect and show this leads to opposite behavior of $p_{\rm term}$ in this limit. We also consider a more conservative velocity-independent formulation. Testing in numerical simulations, these directly translate to large effects on predicted star formation histories and stellar masses of massive galaxies, explaining differences between some models and motivating further study in idealized simulations., Comment: 9 pages, 3 figures. Submitted to the Open Journal of Astrophysics. Comments welcome

While galaxy rotation curves provide one of the most powerful methods for measuring dark matter profiles in the inner regions of rotation-supported galaxies, at the dwarf scale there are factors that can complicate this analysis. Given the expectation of a universal profile in dark matter-only simulations, the diversity of observed rotation curves has become an often-discussed issue in Lambda Cold Dark Matter cosmology on galactic scales. We analyze a suite of Feedback in Realistic Environments (FIRE) simulations of $10^{10}-10^{12}$ $M_\odot$ halos with standard cold dark matter, and compare the true circular velocity to rotation curve reconstructions. We find that, for galaxies with well-ordered gaseous disks, the measured rotation curve may deviate from true circular velocity by at most 10% within the radius of the disk. However, non-equilibrium behavior, non-circular motions, and non-thermal and non-kinetic stresses may cause much larger discrepancies of 50% or more. Most rotation curve reconstructions underestimate the true circular velocity, while some reconstructions transiently over-estimate it in the central few kiloparsecs due to dynamical phenomena. We further demonstrate that the features that contribute to these failures are not always visibly obvious in HI observations. If such dwarf galaxies are included in galaxy catalogs, they may give rise to the appearance of "artificial" rotation curve diversity that does not reflect the true variation in underlying dark matter profiles., Comment: 20 pages, 10 figures

Recently, we demonstrated self-consistent formation of strongly-magnetized quasar accretion disks (QADs) from cosmological radiation-magnetohydrodynamic-thermochemical galaxy-star formation simulations, including the full STARFORGE physics shown previously to produce a reasonable IMF under typical ISM conditions. Here we study star formation and the stellar IMF in QADs, on scales from 100 au to 10 pc from the SMBH. We show it is critical to include physics often previously neglected, including magnetic fields, radiation, and (proto)stellar feedback. Closer to the SMBH, star formation is suppressed, but the (rare) stars that do form exhibit top-heavy IMFs. Stars can form only in special locations (e.g. magnetic field switches) in the outer QAD. Protostars accrete their natal cores rapidly but then dynamically decouple from the gas and wander, ceasing accretion on timescales ~100 yr. Their jets control initial core accretion, but the ejecta are swept up into the larger-scale QAD flow without much dynamical effect. The strong tidal environment strongly suppresses common-core multiplicity. The IMF shape depends sensitively on un-resolved dynamics of protostellar disks (PSDs), as the global dynamical times can become incredibly short ($\ll$ yr) and tidal fields are incredibly strong, so whether PSDs can efficiently transport angular momentum or fragment catastrophically at $\lesssim 10$ au scales requires novel PSD simulations to properly address. Most analytic IMF models and analogies with planet formation in PSDs fail qualitatively to explain the simulation IMFs, though we discuss a couple of viable models., Comment: 39 pages, 23 figures, replaced with version accepted to the Open Journal of Astrophysics. Additional images and movies from the simulations available at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html

The observationally-inferred size versus stellar-mass relationship (SMR) for low-mass galaxies provides an important test for galaxy formation models. However, the relationship relies on assumptions that relate observed luminosity profiles to underlying stellar mass profiles. Here we use the Feedback in Realistic Environments simulations of low-mass galaxies to explore how the predicted SMR changes depending on whether one uses star-particle counts directly or mock observations. We reproduce the SMR found in The Exploration of Local Volume Satellites survey remarkably well only when we infer stellar masses and sizes using mock observations. However, when we use star particles to directly infer stellar masses and half-mass radii, we find that our galaxies are too large and obey a SMR with too little scatter compared to observations. This discrepancy between the "true" galaxy size and mass and those derived in the mock observation approach is twofold. First, our simulated galaxies have higher and more varied MLRs at a fixed colour than those commonly-adopted, which tends to underestimate their stellar masses compared to their true, simulated values. Second, our galaxies have radially increasing MLR gradients therefore using a single MLR tends to under-predict the mass in the outer regions. Similarly, the true half-mass radius is larger than the half-light radius because the light is more concentrated than the mass. If our simulations are accurate representations of the real universe, then the relationship between galaxy size and stellar mass is even tighter for low-mass galaxies than is commonly inferred from observed relations., Comment: 11 pages, 7 figures, 1 table

We present a series of studies on the connection between galaxy morphology and the structure of host dark-matter (DM) haloes using cosmological simulations. In this work, we introduce a new kinematic decomposition scheme that features physical identification of morphological components, enabling robust separation of thin and thick discs; and measure a wide range of halo properties, including their locations in the cosmic web, internal structures, and assembly histories. Our analysis of the TNG50 simulation reveals that the orbital-circularity threshold for disc differentiation varies across galaxies, with systematic trends in mass and redshift, so the widely used decomposition method with constant circularity cuts is oversimplified and underestimates thin disc at JWST redshifts. The energy threshold between the stellar halo and the inner galaxy is also a function of mass and redshift, minimizing at the sub-Galactic halo mass, where the circularity threshold peaks. Revisiting the issue of galaxy size predictor, we show that disc sizes in TNG50 exhibit correlations with three structural parameters besides virial mass and redshift: 1) a positive correlation with halo spin $\lambda$ across redshifts -- stronger than previously reported for zoom-in simulations but still weaker than the simple $r_{1/2}/R_{\rm vir} \propto \lambda$ scaling; 2) an anti-correlation with DM concentration $c$ that is well described by $r_{1/2}/R_{\rm vir} \propto c^{-0.7}$ even when $c$ is measured in the DM only run; 3) more actively accreting haloes having slightly larger discs, as well as more significant stellar haloes and lower thin-to-thick ratio. Disc mass fraction is higher in rounder haloes and in cosmic knots and filaments, implying that disc development needs both stable halo conditions and continuous material supply. Our methodology is public and adaptable to other simulations., Comment: 20 pages, 17 figures, submitted to MNRAS

The cold ($\sim 10^{4}\,{\rm K}$) component of the circumgalactic medium (CGM) accounts for a significant fraction of all galactic baryons. However, using current galaxy-scale simulations to determine the origin and evolution of cold CGM gas poses a significant challenge, since it is computationally infeasible to directly simulate a galactic halo alongside the sub-pc scales that are crucial for understanding the interactions between cold CGM gas and the surrounding ''hot'' medium. In this work, we introduce a new approach: the Cold Gas Subgrid Model (CGSM), which models unresolved cold gas as a second fluid in addition to the standard ''normal'' gas fluid. The CGSM tracks the total mass density and bulk momentum of unresolved cold gas, deriving the properties of its unresolved cloudlets from the resolved gas phase. The interactions between the subgrid cold fluid and the resolved fluid are modeled by prescriptions from high-resolution simulations of ''cloud crushing'' and thermal instability. Through a series of idealized tests, we demonstrate the CGSM's ability to overcome the resolution limitations of traditional hydrodynamics simulations, successfully capturing the correct cold gas mass, its spatial distribution, and the timescales for cloud destruction and growth. We discuss the implications of using this model in cosmological simulations to more accurately represent the microphysics that govern the galactic baryon cycle., Comment: Submitted to MNRAS

Observations indicate dust populations vary between galaxies and within them, suggesting a complex life cycle and evolutionary history. Here we investigate the evolution of galactic dust populations across cosmic time using a suite of cosmological zoom-in simulations from the Feedback in Realistic Environments (FIRE) project, spanning $M_{\rm vir}=10^{9-12}M_{\odot};\,M_{*}=10^{6-11}\,M_{\odot}$. Our simulations incorporate a dust evolution model that accounts for the dominant sources of dust production, growth, and destruction and follows the evolution of specific dust species. All galactic dust populations in our suite exhibit similar evolutionary histories, with gas-dust accretion being the dominant producer of dust mass for all but the most metal-poor galaxies. Similar to previous works, we find the onset of efficient gas-dust accretion occurs above a `critical' metallicity threshold ($Z_{\rm crit}$). Due to this threshold, our simulations reproduce observed trends between galactic D/Z and metallicity and element depletion trends in the ISM. However, we find $Z_{\rm crit}$ varies between dust species due to differences in key element abundances, dust physical properties, and life cycle processes resulting in $Z_{\rm crit}\sim0.05Z_{\odot},\,0.2Z_{\odot},\,0.5Z_{\odot}$ for metallic iron, silicates, and carbonaceous dust, respectively. These variations could explain the lack of small carbonaceous grains observed in the Magellanic Clouds. We also find a delay between the onset of gas-dust accretion and when a dust population reaches equilibrium, which we call the equilibrium timescale ($\tau_{\rm eq}$). The relation between $\tau_{\rm eq}$ and the metal enrichment timescale of a galaxy, determined by its recent evolutionary history, can contribute to the scatter in the observed relation between galactic D/Z and metallicity., Comment: 24 pages, 18 figures, accepted for publication in MNRAS. Minor changes

We make an in-depth analysis of different AGN jet models' signatures, inducing quiescence in galaxies with a halo mass of $10^{12} M_\odot$. 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. We examine the distribution of Mg II, O VI, and O VIII ions, alongside gas temperature and density profiles. Low-energy ions, like Mg II, concentrate in the ISM, while higher energy ions, e.g., O VIII, 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 O VI and O VIII is demonstrated in the ISM-to-CGM transition (between 10-30 kpc), growing smoothly towards the CGM at greater distances., Comment: 17 pages, 13 figures

Carbon monoxide (CO) is a poor tracer of H$_{2}$ in the diffuse interstellar medium (ISM), where most of the carbon is not incorporated into CO molecules unlike the situation at higher extinctions. We present a novel, indirect method to constrain H$_{2}$ column densities ($N_{H_{2}}$) without employing CO observations. We show that previously-recognized nonlinearities in the relation between the extinction, $A_{V} ({H}_{2})$, derived from dust emission and the HI column density ($N_{HI}$) are due to the presence of molecular gas. We employ archival $N_{H_{2}}$ data, obtained from the UV spectra of stars, and calculate $A_{V} ({H}_{2})$ towards these sight lines using 3D extinction maps. We derive an empirical relation between $A_{V} ({H}_{2})$ and $N_{H_{2}}$, which we use to constrain $N_{H_{2}}$ in the diffuse ISM. We construct a $N_{H_{2}}$ map of our Galaxy and compare it to the CO integrated intensity ($W_{CO}$) distribution. We find that the average ratio ($X_{CO}$) between $N_{H_{2}}$ and $W_{CO}$ is approximately equal to $2 \times 10^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$, consistent with previous estimates. However, we find that the $X_{CO}$ factor varies by orders of magnitude on arcminute scales between the outer and the central portions of molecular clouds. For regions with $N_{H_{2}} \gtrsim 10^{20}$ cm$^{-2}$, we estimate that the average H$_{2}$ fractional abundance, $f_{{H}_{2}}$ = $2N_{H_{2}}$/(2$N_{H_{2}}$ + $N_{{HI}}$), is 0.25. Multiple (distinct) largely atomic clouds are likely found along high-extinction sightlines ($A_{V} \geq 1$ mag), hence limiting $f_{{H}_{2}}$ in these directions. More than $50 \%$ of the lines of sight with $N_{H_{2}} \geq 10^{20}$ cm$^{-2}$ are untraceable by CO with a $J$ = 1-0 sensitivity limit $W_{CO} = 1$ K km s$^{-1}$., Comment: Accepted for publication in A&A (edited author list)

Using cosmological hydrodynamical zoom-in simulations, we explore the properties of subhalos in Milky Way analogs that contain a sub-component of Atomic Dark Matter (ADM). ADM differs from Cold Dark Matter (CDM) due to the presence of self interactions that lead to energy dissipation and bound-state formation, analogous to Standard Model baryons. This model can arise in complex dark sectors that are natural and theoretically-motivated extensions to the Standard Model. The simulations used in this work were carried out using GIZMO and utilize the FIRE-2 galaxy formation physics in the Standard Model baryonic sector. For the parameter points we consider, the ADM gas cools efficiently, allowing it to collapse to the center of subhalos. This increases a subhalo's central density and affects its orbit, with more subhalos surviving small pericentric passages. The subset of subhalos that host visible satellite galaxies have cuspier density profiles and smaller stellar-half-mass radii relative to CDM. The entire population of dwarf galaxies produced in the ADM simulations is much more compact than those seen in CDM simulations, unable to reproduce the entire diversity of observed dwarf galaxy structures. Additionally, we also identify a population of highly compact subhalos that consist nearly entirely of ADM and form in the central region of the host, where they can leave distinctive imprints in the baryonic disk. This work presents the first detailed exploration of subhalo properties in a strongly dissipative dark matter scenario, providing intuition for how other regions of ADM parameter space, as well as other dark sector models, would impact galactic-scale observables., Comment: 18 pages, 10 figures, 8 pages of appendices

We investigate the formation of dense stellar clumps in a suite of high-resolution cosmological zoom-in simulations of a massive, star forming galaxy at $z \sim 2$ under the presence of strong quasar winds. Our simulations include multi-phase ISM physics from the Feedback In Realistic Environments (FIRE) project and a novel implementation of hyper-refined accretion disk winds. We show that powerful quasar winds can have a global negative impact on galaxy growth while in the strongest cases triggering the formation of an off-center clump with stellar mass ${\rm M}_{\star}\sim 10^{7}\,{\rm M}_{\odot}$, effective radius ${\rm R}_{\rm 1/2\,\rm Clump}\sim 20\,{\rm pc}$, and surface density $\Sigma_{\star} \sim 10^{4}\,{\rm M}_{\odot}\,{\rm pc}^{-2}$. The clump progenitor gas cloud is originally not star-forming, but strong ram pressure gradients driven by the quasar winds (orders of magnitude stronger than experienced in the absence of winds) lead to rapid compression and subsequent conversion of gas into stars at densities much higher than the average density of star-forming gas. The AGN-triggered star-forming clump reaches ${\rm SFR} \sim 50\,{\rm M}_{\odot}\,{\rm yr}^{-1}$ and $\Sigma_{\rm SFR} \sim 10^{4}\,{\rm M}_{\odot}\,{\rm yr}^{-1}\,{\rm kpc}^{-2}$, converting most of the progenitor gas cloud into stars in $\sim$2\,Myr, significantly faster than its initial free-fall time and with stellar feedback unable to stop star formation. In contrast, the same gas cloud in the absence of quasar winds forms stars over a much longer period of time ($\sim$35\,Myr), at lower densities, and losing spatial coherency. The presence of young, ultra-dense, gravitationally bound stellar clumps in recently quenched galaxies could thus indicate local positive feedback acting alongside the strong negative impact of powerful quasar winds, providing a plausible formation scenario for globular clusters., Comment: 14 pages, 12 figures. Accepted for publication in MNRAS; doi:10.1093/mnras/stae1021

In the absence of supplementary heat, the radiative cooling of halo gas around massive galaxies (Milky Way mass and above) leads to an excess of cold gas or stars beyond observed levels. AGN jet-induced heating is likely essential, but the specific properties of the jets remain unclear. Our previous work (Su et al. 2021) concludes from simulations of a halo with $10^{14} M_\odot$ that a successful jet model should have an energy flux comparable to the free-fall energy flux at the cooling radius and should inflate a sufficiently wide cocoon with a long enough cooling time. In this paper, we investigate three jet modes with constant fluxes satisfying the criteria, including high-temperature thermal jets, cosmic ray (CR)-dominant jets, and widely precessing kinetic jets in $10^{12}-10^{15}\,{\rm M}_{\odot}$ halos using high-resolution, non-cosmological MHD simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We find that scaling the jet energy according to the free-fall energy at the cooling radius can successfully suppress the cooling flows and quench galaxies without obviously violating observational constraints. We investigate an alternative scaling method in which we adjust the energy flux based on the total cooling rate within the cooling radius. However, we observe that the strong interstellar medium (ISM) cooling dominates the total cooling rate in this scaling approach, resulting in a jet flux that exceeds the amount needed to suppress the cooling flows. With the same energy flux, the CR-dominant jet is most effective in suppressing the cooling flow across all the surveyed halo masses due to the enhanced CR pressure support. We confirm that the criteria for a successful jet model, which we proposed in Su et al. (2021), work across a much wider range, encompassing halo masses of $10^{12}-10^{15} {\rm M_\odot}$., Comment: 15 pages, 10 figures, Submitted to MNRAS; comments welcome!

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. 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 FIRE-3 cosmological zoom-in simulations with AGN feedback evolved to $z\sim0$, examining the impact of AGN feedback on a set of galaxies with halos in the mass range $10^{12}-10^{13} M_{\odot}$. These simulations combine detailed stellar and ISM physics with multi-channel 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 the $M_{\rm BH}$-$\sigma$ 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 over-quench the lower-mass galaxies. This indicates that further refinements of the AGN modeling are needed., Comment: 18 pages, 8 figures, accepted by ApJ

Recent numerical cosmological radiation-magnetohydrodynamic-thermochemical-star formation simulations have resolved the formation of quasar accretion disks with Eddington or super-Eddington accretion rates onto supermassive black holes (SMBHs) down to a few hundred gravitational radii. These 'flux-frozen' and hyper-magnetized disks appear to be qualitatively distinct from classical $\alpha$ disks and magnetically-arrested disks: the midplane pressure is dominated by toroidal magnetic fields with plasma $\beta \ll 1$ powered by advection of magnetic flux from the interstellar medium (ISM), and they are super-sonically and trans-Alfvenically turbulent with cooling times short compared to dynamical times yet remain gravitationally stable owing to magnetic support. In this paper, we present a simple analytic similarity model for such disks. For reasonable assumptions, the model is entirely specified by the boundary conditions (inflow rate at the BH radius of influence [BHROI]). We show that the scalings from this model are robust to various detailed assumptions, agree remarkably well with the simulations (given their simplicity), and demonstrate the self-consistency and gravitational stability of such disks even in the outer accretion disk (approaching the BHROI) at hyper-Eddington accretion rates., Comment: 12 pages, 1 figure. Replaced with accepted version. Companion paper to 'FORGE'D IN FIRE II' (arXiv:2310.04506, part of a series with arXiv:2309.13115 -- animations of the simulations referred to here can be viewed at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html)

In a companion paper, we reported the self-consistent formation of quasar accretion disks with inflow rates $\sim 10\,{\rm M_{\odot}\,yr^{-1}}$ down to <300 Schwarzschild radii from cosmological radiation-magneto-thermochemical-hydrodynamical galaxy and star formation simulations. We see the formation of a well-defined, steady-state accretion disk which is stable against star formation at sub-pc scales. The disks are optically thick, with radiative cooling balancing accretion, but with properties that are distinct from those assumed in most previous accretion disk models. The pressure is strongly dominated by (primarily toroidal) magnetic fields, with a plasma $\beta \sim 10^{-4}$ even in the disk midplane. They are qualitatively distinct from magnetically elevated or arrested disks. The disks are strongly turbulent, with trans-Alfvenic and highly super-sonic turbulence, and balance this via a cooling time that is short compared to the disk dynamical time, and can sustain highly super-Eddington accretion rates. Their surface and 3D densities at $\sim 10^{3}-10^{5}$ gravitational radii are much lower than in a Shakura-Sunyaev disk, with important implications for their thermo-chemistry and stability. We show how the magnetic field strengths and geometries arise from rapid advection of flux with the inflow from much weaker galaxy-scale fields in these 'flux-frozen' disks, and how this stabilizes the disk and gives rise to efficient torques. Re-simulating without magnetic fields produces catastrophic fragmentation with a vastly smaller, lower-$\dot{M}$ Shakura-Sunyaev-like disk., Comment: 46 pages, 30 figures. Updated to match accepted version from The Open Journal of Astrophysics. Part of a series with arXiv:2309.13115 . Animations of the simulations can be viewed at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html

Understanding what shapes the cold gas component of galaxies, which both provides the fuel for star formation and is strongly affected by the subsequent stellar feedback, is a crucial step towards a better understanding of galaxy evolution. Here, we analyse the HI properties of a sample of 46 Milky Way halo-mass galaxies, drawn from cosmological simulations (EMP-Pathfinder and FIREbox). This set of simulations comprises galaxies evolved self-consistently across cosmic time with different baryonic sub-grid physics: three different star formation models [constant star formation efficiency (SFE) with different star formation eligibility criteria, and an environmentally-dependent, turbulence-based SFE] and two different feedback prescriptions, where only one sub-sample includes early stellar feedback. We use these simulations to assess the impact of different baryonic physics on the HI content of galaxies. We find that the galaxy-wide HI properties agree with each other and with observations. However, differences appear for small-scale properties. The thin HI discs observed in the local Universe are only reproduced with a turbulence-dependent SFE and/or early stellar feedback. Furthermore, we find that the morphology of HI discs is particularly sensitive to the different physics models: galaxies simulated with a turbulence-based SFE have discs that are smoother and more rotationally symmetric, compared to those simulated with a constant SFE; galaxies simulated with early stellar feedback have more regular discs than supernova-feedback-only galaxies. We find that the rotational asymmetry of the HI discs depends most strongly on the underlying physics model, making this a promising observable for understanding the physics responsible for shaping the interstellar medium of galaxies., Comment: 13 pages, 10 figures + appendices (7 pages, 10 figures); updated to version accepted by MNRAS

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, gas accretion events, in presence of turbulent Inter Stellar Medium (ISM) with multiple star formation episodes, in addition to coupling to their host galaxies' dark matter halos. We investigate bar formation in 13 Milky Way-mass galaxies from the FIRE-2 (Feedback in Realistic Environments) 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 $\sim 1.53$ kpc), have a wide range of pattern speeds (36--97 km s$^{-1}$kpc$^{-1}$), and live for a wide range of dynamical times (2--160 bar rotations). We find that 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, and the bars formed are 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 strengths $A_2/A_0 \sim$ 0.1--0.2., Comment: 36 pages (including Appendix), 22 figures, to be submitted to the Astrophysical Journal

Cosmic rays (CRs) may drive outflows and alter the phase structure of the circumgalactic medium, with potentially important implications on galaxy formation. However, these effects ultimately depend on the dominant mode of transport of CRs within and around galaxies, which remains highly uncertain. To explore potential observable constraints on CR transport, we investigate a set of cosmological FIRE-2 CR-MHD simulations of L$_{\ast}$ galaxies which evolve CRs with transport models motivated by self-confinement (SC) and extrinsic turbulence (ET) paradigms. To first order, the synchrotron properties diverge between SC and ET models due to a CR physics driven hysteresis. SC models show a higher tendency to undergo `ejective' feedback events due to a runaway buildup of CR pressure in dense gas due to the behavior of SC transport scalings at extremal CR energy densities. The corresponding CR wind-driven hysteresis results in brighter, smoother, and more extended synchrotron emission in SC runs relative to ET and constant diffusion runs. The differences in synchrotron arise from different morphology, ISM gas and \textbf{B} properties, potentially ruling out SC as the dominant mode of CR transport in typical star-forming L$_{\ast}$ galaxies, and indicating the potential for non-thermal radio continuum observations to constrain CR transport physics., Comment: 5 pages, 3 figures. Accepted version

It has recently become possible to zoom-in from cosmological to sub-pc scales in galaxy simulations to follow accretion onto supermassive black holes (SMBHs). However, at some point the approximations used on ISM scales (e.g. optically-thin cooling and stellar-population-integrated star formation [SF] and feedback [FB]) break down. We therefore present the first cosmological radiation-magnetohydrodynamic (RMHD) simulation which self-consistently combines the FIRE physics (relevant on galactic/ISM scales where SF/FB are ensemble-averaged) and STARFORGE physics (relevant on small scales where we track individual (proto)stellar formation and evolution), together with explicit RMHD (including non-ideal MHD and multi-band M1-RHD) which self-consistently treats both optically-thick and thin regimes. This allows us to span scales from ~100 Mpc down to <100 au (~300 Schwarzschild radii) around a SMBH at a time where it accretes as a bright quasar, in a single simulation. We show that accretion rates up to $\sim 10-100\,{\rm M_{\odot}\,yr^{-1}}$ can be sustained into the accretion disk at $\ll 10^{3}\,R_{\rm schw}$, with gravitational torques between stars and gas dominating on sub-kpc scales until star formation is shut down on sub-pc scales by a combination of optical depth to cooling and strong magnetic fields. There is an intermediate-scale, flux-frozen disk which is gravitoturbulent and stabilized by magnetic pressure sustaining strong turbulence and inflow with persistent spiral modes. In this paper we focus on how gas gets into the small-scale disk, and how star formation is efficiently suppressed., Comment: 38 pages, 18 figures. Replaced with accepted version

We present predictions for the high-redshift halo-galaxy-supermassive black hole (SMBH) connection from the TRINITY model. Constrained by a comprehensive compilation of galaxy ($0\leq z \leq 10$) and SMBH datasets ($0\leq z \leq 6.5$), TRINITY finds: 1) The number of SMBHs with $M_\bullet > 10^9 M_\odot$ in the observable Universe increases by six orders of magnitude from $z\sim10$ to $z\sim2$, and by another factor of $\sim 3$ from $z\sim2$ to $z=0$; 2) The $M_\bullet > 10^9/10^{10} M_\odot$ SMBHs at $z\sim 6$ live in haloes with $\sim (2-3)/(3-5) \times 10^{12} M_\odot$; 3) the new JWST AGNs at $7\lesssim z \lesssim 11$ are broadly consistent with the median SMBH mass-galaxy mass relation for AGNs from TRINITY; 4) Seeds from runaway mergers in nuclear star clusters are viable progenitors for the SMBHs in GN-z11 ($z=10.6$) and CEERS_1019 ($z=8.7$); 5) $z=6-10$ quasar luminosity functions from wide area surveys by, e.g., Roman and Euclid, will reduce uncertainties in the $z=6-10$ SMBH mass-galaxy mass relation by up to $\sim 0.5$ dex., Comment: 15 pages, 12 figures, submitted to MNRAS. Questions and comments are welcome!

Synchrotron emission is one of few observable tracers of galactic magnetic fields (\textbf{B}) and cosmic rays (CRs). Much of our understanding of \textbf{B} in galaxies comes from utilizing synchrotron observations in conjunction with several simplifying assumptions of equipartition models, however it remains unclear how well these assumptions hold, and what \textbf{B} these estimates physically represent. Using FIRE simulations which self consistently evolve CR proton, electron, and positron spectra from MeV to TeV energies, we present the first synthetic synchrotron emission predictions from simulated L$_{*}$ galaxies with "live" spectrally-resolved CR-MHD. We find that synchrotron emission can be dominated by relatively cool and dense gas, resulting in equipartition estimates of \textbf{B} with fiducial assumptions underestimating the "true" \textbf{B} in the gas that contributes the most emission by factors of 2-3 due to small volume filling factors. Motivated by our results, we present an analytic framework that expands upon equipartition models for estimating \textbf{B} in a multi-phase medium. Comparing our spectrally-resolved synchrotron predictions to simpler spectral assumptions used in galaxy simulations with CRs, we find that spectral evolution can be crucial for accurate synchrotron calculations towards galactic centers, where loss terms are large., Comment: 10 pages, 5 figures (with 3 additional figures in the Appendix). Accepted version (MNRAS December 21, 2023)

The opacity limit is an important concept in star formation: isothermal collapse cannot proceed without limit, because eventually cooling radiation is trapped and the temperature rises quasi-adiabatically, setting a minimum Jeans mass $M_{\rm J}^{\rm min}$. Various works have considered this scenario and derived expressions for $M_{\rm J}^{\rm min}$, generally $\sim 10^{-3}-10^{-2}M_\odot$ in normal star-forming conditions, but with conflicting results about the scaling with ambient conditions and material properties. We derive expressions for the thermal evolution of dust-cooled collapsing gas clumps in various limiting cases, given a general ambient radiation field ($u_{\rm rad}$, $T_{\rm rad}$) and a general power-law dust opacity law $\sigma_{\rm d} = A_{\rm d} T^{\beta}$. By accounting for temperature evolution self-consistently we rule out a previously-proposed regime in which the adiabatic transition occurs while the core is still optically-thin. If the radiation field is weak or dust opacity is small, $M_{\rm J}^{\rm min}$ is insensitive to dust properties/abundance ($\sim A_{\rm d}^{-\frac{1}{11}}-A_{\rm d}^{-\frac{1}{15}}$), but if the radiation field is strong and dust is abundant it scales $\propto A_{\rm d}^{1/3}$. This could make the IMF less bottom-heavy in dust-rich and/or radiation-dense environments, e.g. galactic centers, starburst galaxies, massive high-$z$ galaxies, and proto-star clusters that are already luminous., Comment: To be submitted to the Open Journal of Astrophysics. Comments welcome

Cosmic rays (CRs) with energies $\ll$ TeV comprise a significant component of the interstellar medium (ISM). Major uncertainties in CR behavior on observable scales (much larger than CR gyroradii) stem from how magnetic fluctuations scatter CRs in pitch angle. Traditional first-principles models, which assume these magnetic fluctuations are weak and uniformly scatter CRs in a homogeneous ISM, struggle to reproduce basic observables such as the dependence of CR residence times and scattering rates on rigidity. We therefore explore a new category of "patchy" CR scattering models, wherein CRs are predominantly scattered by intermittent strong scattering structures with small volume-filling factors. These models produce the observed rigidity dependence with a simple size distribution constraint, such that larger scattering structures are rarer but can scatter a wider range of CR energies. To reproduce the empirically-inferred CR scattering rates, the mean free path between scattering structures must be $\ell_{\rm mfp} \sim 10$ pc at GeV energies. We derive constraints on the sizes, internal properties, mass/volume-filling factors, and the number density any such structures would need to be both physically and observationally consistent. We consider a range of candidate structures, both large-scale (e.g. H II regions) and small-scale (e.g. intermittent turbulent structures, perhaps even associated with radio plasma scattering) and show that while many macroscopic candidates can be immediately ruled out as the primary CR scattering sites, many smaller structures remain viable and merit further theoretical study. We discuss future observational constraints that could test these models., Comment: 9 pages, 3 figures, accepted to MNRAS

The Radial Acceleration Relation (RAR) connects the total gravitational acceleration of a galaxy at a given radius, $a_{\rm tot}(r)$, with that accounted for by baryons at the same radius, $a_{\rm bar}(r)$. The shape and tightness of the RAR for rotationally-supported galaxies have characteristics in line with MOdified Newtonian Dynamics (MOND) and can also arise within the Cosmological Constant + Cold Dark Matter ($\Lambda$CDM) paradigm. We use zoom simulations of 20 galaxies with stellar masses of $M_{\star} \, \simeq 10^{7-11} \, M_{\odot}$ to study the RAR in the \texttt{FIRE-2} simulations. We highlight the existence of simulated galaxies with non-monotonic RAR tracks that ``hook'' down from the average relation. These hooks are challenging to explain in Modified Inertia theories of MOND, but naturally arise in all of our \lcdm-simulated galaxies that are dark-matter dominated at small radii and have feedback-induced cores in their dark matter haloes. We show, analytically and numerically, that downward hooks are expected in such cored haloes because they have non-monotonic acceleration profiles. We also extend the relation to accelerations below those traced by disc galaxy rotation curves. In this regime, our simulations exhibit ``bends'' off of the MOND-inspired extrapolation of the RAR, which, at large radii, approach $a_{\rm tot} \, \approx \, a_{\rm bar} \, /f_{\rm b}$, where $f_{\rm b}$ is the cosmic baryon fraction. Future efforts to search for these hooks and bends in real galaxies will provide interesting tests for MOND and $\Lambda$CDM., Comment: 15 pages, 9 figures; Accepted for publication in MNRAS

When a detailed model of a stellar population is unavailable, it is most common to assume that stellar masses are independently and identically distributed according to some distribution: the universal initial mass function (IMF). However, stellar masses resulting from causal, long-ranged physics cannot be truly random and independent, and the IMF may vary with environment. To compare stochastic sampling with a physical model, we run a suite of 100 STARFORGE radiation magnetohydrodynamics simulations of low-mass star cluster formation in $2000M_\odot$ clouds that form $\sim 200$ stars each on average. The stacked IMF from the simulated clouds has a sharp truncation at $\sim 28 M_\odot$, well below the typically-assumed maximum stellar mass $M_{\rm up} \sim 100-150M_\odot$ and the total cluster mass. The sequence of star formation is not totally random: massive stars tend to start accreting sooner and finish later than the average star. However, final cluster properties such as maximum stellar mass and total luminosity have a similar amount of cloud-to-cloud scatter to random sampling. Therefore stochastic sampling does not generally model the stellar demographics of a star cluster as it is forming, but may describe the end result fairly well, if the correct IMF -- and its environment-dependent upper cutoff -- are known., Comment: Accepted for publication in Open Journal of Astrophysics

We present the redshift evolution of quasar luminosity functions decomposed by halo mass, galaxy mass, supermassive black hole (SMBH) mass, and Eddington ratio, as well as SMBH kinetic/radiative energy output ratios from TRINITY, a flexible empirical model that self-consistently infers the halo--galaxy--SMBH connection that match observational data. Key findings include: 1) The normalization of QLF increases by ~3-4 dex from z~10 to z~4, due to the fast mass build-up of different SMBH populations; 2) From z~4 to z~1, less massive galaxies and SMBHs make up bigger and bigger fractions of QLFs, due to the AGN downsizing effect; 3) At z~0, massive haloes/galaxies/SMBHs are responsible for most bright quasars due to low Eddington ratios among all SMBHs; 4) The bright ends of quasar luminosity functions (QLFs) are dominated by SMBHs that are at least 0.3 dex over-massive relative to the median SMBH mass-galaxy mass relation; 5) QLFs at z~6-7 are dominated by SMBHs accreting at Eddington ratios 0.1 < $\eta_\mathrm{rad}$ < 1, but super-Eddington AGNs contribute more significantly to QLFs towards z~9-10., Comment: 18 pages, 14 figures. Accepted by MNRAS. Comments welcome!

Dark sector theories naturally lead to multi-component scenarios for dark matter where a sub-component can dissipate energy through self-interactions, allowing it to efficiently cool inside galaxies. We present the first cosmological hydrodynamical simulations of Milky Way analogues where the majority of dark matter is collisionless Cold Dark Matter (CDM), but a sub-component (6%) is strongly dissipative minimal Atomic Dark Matter (ADM). The simulations, implemented in GIZMO and utilizing FIRE-2 galaxy formation physics to model the standard baryonic sector, demonstrate that the addition of even a small fraction of dissipative dark matter can significantly impact galactic evolution despite being consistent with current cosmological constraints. We show that ADM gas with roughly Standard-Model-like masses and couplings can cool to form a rotating "dark disk" with angular momentum closely aligned with the visible stellar disk. The morphology of the disk depends sensitively on the parameters of the ADM model, which affect the cooling rates in the dark sector. The majority of the ADM gas gravitationally collapses into dark "clumps" (regions of black hole or mirror star formation), which form a prominent bulge and a rotating thick disk in the central galaxy. These clumps form early and quickly sink to the inner ~kpc of the galaxy, affecting the galaxy's star-formation history and present-day baryonic and CDM distributions., Comment: 9 pages, 4 figures, 12 pages of appendices and supplementary figures

We study the kinematics of stars both at their formation and today within 14 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. We quantify the relative importance of cosmological disk settling and post-formation dynamical heating. We identify three eras: a Pre-Disk Era (typically >8 Gyr ago), when stars formed on dispersion-dominated orbits; an Early-Disk Era (~ 8 - 4 Gyr ago), when stars started to form on rotation-dominated orbits but with high velocity dispersion, sigma_form; and a Late-Disk Era (< 4 Gyr ago), when stars formed with low sigma_form. sigma_form increased with time during the Pre-Disk Era, peaking ~ 8 Gyr ago, then decreased throughout the Early-Disk Era as the disk settled and remained low throughout the Late-Disk Era. By contrast, the velocity dispersion measured today, sigma_now, increases monotonically with age because of stronger post-formation heating for Pre-Disk stars. Importantly, most of sigma_now was in place at formation, not added post-formation, for stars younger than ~ 10 Gyr. We compare the evolution of the three velocity components: at all times, sigma_R,form > sigma_phi,form > sigma_Z,form. Post-formation heating primarily increased sigma_R at ages < 4 Gyr but acted nearly isotropically for older stars. The kinematics of young stars in FIRE-2 broadly agree with the range observed across the MW, M31, M33, and PHANGS-MUSE galaxies. The lookback time that the disk began to settle correlates with its dynamical state today: earlier-settling galaxies currently form colder disks. Including stellar cosmic-ray feedback does not significantly change disk rotational support at fixed stellar mass., Comment: 23 pages, accepted for publication in MNRAS

Using recent empirical constraints on the dark matter halo--galaxy--supermassive black hole (SMBH) connection from $z=0-7$, we infer how undermassive, typical, and overmassive SMBHs contribute to the quasar luminosity function (QLF) at $z=6$. We find that beyond $L_\mathrm{bol} = 5 \times 10^{46}$ erg/s, the $z=6$ QLF is dominated by SMBHs that are at least 0.3 dex above the $z=6$ median $M_\bullet-M_*$ relation. The QLF is dominated by typical SMBHs (i.e., within $\pm 0.3$ dex around the $M_\bullet-M_*$ relation) at $L_\mathrm{bol} \lesssim 10^{45}$ erg/s. At $z\sim 6$, the intrinsic $M_\bullet-M_*$ relation for all SMBHs is slightly steeper than the $z=0$ scaling, with a similar normalization at $M_* \sim 10^{11} M_\odot$. We also predict the $M_\bullet-M_*$ relation for $z=6$ bright quasars selected by different bolometric luminosity thresholds, finding very good agreement with observations. For quasars with $L_\mathrm{bol} > 3 \times 10^{46}$ ($10^{48}$) erg/s, the scaling relation is shifted upwards by $\sim0.35$ (1.0) dex for $10^{11} M_\odot$ galaxies. To accurately measure the intrinsic $M_\bullet-M_*$ relation, it is essential to include fainter quasars with $L_\mathrm{bol} \lesssim 10^{45}$ erg/s. At high redshifts, low-luminosity quasars are thus the best targets for understanding typical formation paths for SMBHs in galaxies., Comment: 5 pages, 3 figures. Submitted to MNRAS Letters. Comments welcome!

As they grow, galaxies can transition from irregular/spheroidal with 'bursty' star formation histories (SFHs), to disky with smooth SFHs. But even in simulations, the direct physical cause of such transitions remains unclear. We therefore explore this in a large suite of numerical experiments re-running portions of cosmological simulations with widely varied physics, further validated with existing FIRE simulations. We show that gas supply, cooling/thermodynamics, star formation model, Toomre scale, galaxy dynamical times, and feedback properties do not have a direct causal effect on these transitions. Rather, both the formation of disks and cessation of bursty star formation are driven by the gravitational potential, but in different ways. Disk formation is promoted when the mass profile becomes sufficiently centrally-concentrated in shape (relative to circularization radii): we show that this provides a well-defined dynamical center, ceases to support the global 'breathing modes' which can persist indefinitely in less-concentrated profiles and efficiently destroy disks, promotes orbit mixing to form a coherent angular momentum, and stabilizes the disk. Smooth SF is promoted by the potential or escape velocity (not circular velocity) becoming sufficiently large at the radii of star formation that cool, mass-loaded (momentum-conserving) outflows are trapped/confined near the galaxy, as opposed to escaping after bursts. We discuss the detailed physics, how these conditions arise in cosmological contexts, their relation to other correlated phenomena (e.g. inner halo virialization, vertical disk 'settling'), and observations., Comment: Updated to match MNRAS accepted version. 45 pages, 32 figures

Observations of local star-forming galaxies (SFGs) show a tight correlation between their singly ionized carbon line luminosity ($L_{\rm [C_{II}]}$) and star formation rate (SFR), suggesting that $L_{\rm [C_{II}]}$ may be a useful SFR tracer for galaxies. Some other galaxy populations, however, are found to have lower $L_{\rm [C_{II}]}{}/{}\rm SFR$ than the local SFGs, including the infrared-luminous, starburst galaxies at low and high redshifts, as well as some moderately star-forming galaxies at the epoch of re-ionization (EoR). The origin of this `$\rm [C_{II}]$ deficit' is unclear. In this work, we study the $L_{\rm [C_{II}]}$-SFR relation of galaxies using a sample of $z=0-8$ galaxies with $M_*\approx10^7-5\times10^{11}\,M_\odot$ extracted from cosmological volume and zoom-in simulations from the Feedback in Realistic Environments (FIRE) project. We find a simple analytic expression for $L_{\rm [C_{II}]}$/SFR of galaxies in terms of the following parameters: mass fraction of $\rm [C_{II}]$-emitting gas ($f_{\rm [C_{II}]}$), gas metallicity ($Z_{\rm gas}$), gas density ($n_{\rm gas}$) and gas depletion time ($t_{\rm dep}{}={}M_{\rm gas}{}/{}\rm SFR$). We find two distinct physical regimes, where $t_{\rm dep}$ ($Z_{\rm gas}$) is the main driver of the $\rm [C_{II}]$ deficit in $\rm H_2$-rich ($\rm H_2$-poor) galaxies. The observed $\rm [C_{II}]$ deficit of IR-luminous galaxies and early EoR galaxies, corresponding to the two different regimes, is due to short gas depletion time and low gas metallicity, respectively. Our result indicates that $\rm [C_{II}]$ deficit is a common phenomenon of galaxies, and caution needs to be taken when applying a constant $L_{\rm [C_{II}]}$-to-SFR conversion factor derived from local SFGs to estimate cosmic SFR density at high redshifts and interpret data from upcoming $\rm [C_{II}]$ line intensity mapping experiments., Comment: 44 page, 24 figures, 8 tables and 8 appendices, accepted for publication in MNRAS. Key figures: fig. 16 & 17. Comments are welcome

Negative feedback from accreting supermassive black holes is regarded as a key ingredient in suppressing star formation and quenching massive galaxies. However, several models and observations suggest that black hole feedback may have a positive effect, triggering star formation by compressing interstellar medium gas to higher densities. We investigate the dual role of black hole feedback using cosmological hydrodynamic simulations from the Feedback In Realistic Environments (FIRE) project, including a novel implementation of hyper-refined accretion-disc winds. Focusing on a massive, star-forming galaxy at $z \sim 2$ ($M_{\rm halo} \sim 10^{12.5} \, {\rm M}_{\odot}$), we show that strong quasar winds with kinetic power $\sim$10$^{46}$ erg/s acting for $>$20$\,$Myr drive the formation of a central gas cavity and can dramatically reduce the star formation rate surface density across the galaxy disc. The suppression of star formation is primarily driven by reducing the amount of gas that can become star-forming, compared to directly evacuating the pre-existing star-forming gas reservoir (preventive feedback dominates over ejective feedback). Despite the global negative impact of quasar winds, we identify several plausible signatures of local positive feedback, including: (1) spatial anti-correlation of wind-dominated regions and star-forming clumps, (2) higher local star formation efficiency in compressed gas near the edge of the cavity, and (3) increased local contribution of outflowing material to star formation. Stars forming under the presence of quasar winds tend to do so at larger radial distances. Our results suggest that positive and negative AGN feedback can coexist in galaxies, but local positive triggering of star formation plays a minor role in global galaxy growth., Comment: 17 pages, 12 figures. Accepted for publication in MNRAS

Modeling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve adaptively, in any high-dynamic range simulation, but this poses unique challenges of consistency, conservation, and accuracy, especially in multi-physics simulations where species with different softening laws may interact. We therefore derive a generalized form of the energy-and-momentum conserving gravitational equations of motion, applicable to arbitrary rules used to determine the force softening, together with consistent associated timestep criteria, interaction terms between species with different softening laws, and arbitrary maximum/minimum softenings. We also derive new methods to maintain better accuracy and conservation when symmetrizing forces between particles. We review and extend previously-discussed adaptive softening schemes based on the local neighbor particle density, and present several new schemes for scaling the softening with properties of the gravitational field, i.e. the potential or acceleration or tidal tensor. We show that the tidal softening scheme not only represents a physically-motivated, translation and Galilean invariant and equivalence-principle respecting (and therefore conservative) method, but imposes negligible timestep or other computational penalties, ensures that pairwise two-body scattering is small compared to smooth background forces, and can resolve outstanding challenges in properly capturing tidal disruption of substructures (minimizing artificial destruction) while also avoiding excessive N-body heating. We make all of this public in the GIZMO code., Comment: 25 pages, 13 figures. Updated to match MNRAS accepted version. Six appendices with additional details and tests expanded

The shape and orientation of dark matter (DM) halos are sensitive to the micro-physics of the DM particle, yet in many mass models, the symmetry axes of the Milky Way's DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well-motivated for the inner DM halo but not for the outer halo. We use zoomed cosmological-baryonic simulations from the Latte suite of FIRE-2 Milky Way-mass galaxies to explore the evolution of the DM halo's orientation with radius and time, with or without a major merger with a Large Magellanic Cloud (LMC) analog, and when varying the DM model. In three of the four CDM halos we examine, the orientation of the halo minor axis diverges from the stellar disk vector by more than 20 degrees beyond about 30 galactocentric kpc, reaching a maximum of 30--90 degrees depending on the individual halo's formation history. In identical simulations using a model of self-interacting DM with $\sigma = 1 \, \mathrm{cm}^2 \, \mathrm{g}^{-1}$, the halo remains aligned with the stellar disk out to $\sim$200--400 kpc. Interactions with massive satellites ($M \gtrsim 4 \times 10^{10} \, \rm{M_\odot}$ at pericenter; $M \gtrsim 3.3 \times 10^{10} \, \rm{M_\odot}$ at infall) affect the orientation of the halo significantly, aligning the halo's major axis with the satellite galaxy from the disk to the virial radius. The relative orientation of the halo and disk beyond 30 kpc is a potential diagnostic of SIDM if the effects of massive satellites can be accounted for., Comment: Submitted to ApJ, 19 pages, 12 figures, 3 tables

Many recent numerical studies have argued that cosmic rays (CRs) from supernovae (SNe) or active galactic nuclei (AGN) could play a crucial role in galaxy formation, in particular by establishing a CR-pressure dominated circum-galactic medium (CGM). But explicit CR-magneto-hydrodynamics (CR-MHD) remains computationally expensive, and it is not clear whether those results can be applied to simulations that do not explicitly treat magnetic fields or resolved ISM phase structure. We therefore present an intentionally extremely-simplified sub-grid model for CRs, which attempts to capture the key qualitative behaviors of greatest interest for those interested in simulations or semi-analytic models including some approximate CR effects on galactic (>kpc) scales, while imposing negligible computational overhead. The model is numerically akin to some recently-developed sub-grid models for radiative feedback, and allows for a simple constant parameterization of the CR diffusivity and/or streaming speed; it allows for an arbitrary distribution of sources (proportional to black hole accretion rates or star-particle SNe rates or gas/galaxy star formation rates), and interpolates between the limits where CRs escape the galaxies with negligible losses and those where CRs lose most of their energy catastrophically before escape (relevant in e.g. starburst galaxies). The numerical equations are solved trivially alongside gravity in most codes. We compare this to explicit CR-MHD simulations and discuss where the (many) sub-grid approximations break down, and what drives the major sources of uncertainty., Comment: 12 pages, 4 figures. Updated to match published MNRAS version

Recent theoretical studies predict that the circumgalactic medium (CGM) around low-redshift, $\sim L_*$ galaxies could have substantial nonthermal pressure support in the form of cosmic rays. However, these predictions are sensitive to the specific model of cosmic-ray transport employed, which is theoretically and observationally underconstrained. In this work, we propose a novel observational constraint for calculating the lower limit of the radially-averaged, effective cosmic-ray transport rate, $\kappa_{\rm min}^{\rm eff}$. Under a wide range of assumptions (so long as cosmic rays do not lose a significant fraction of their energy in the galactic disk, regardless of whether the cosmic-ray pressure is important or not in the CGM), we demonstrate a well-defined relationship between $\kappa_{\rm min}^{\rm eff}$ and three observable galaxy properties: the total hydrogen column density, the average star formation rate, and the gas circular velocity. We use a suite of FIRE-2 galaxy simulations with a variety of cosmic-ray transport physics to demonstrate that our analytic model of $\kappa_{\rm min}^{\rm eff}$ is a robust lower limit of the true cosmic-ray transport rate. We then apply our new model to calculate $\kappa_{\rm min}^{\rm eff}$ for galaxies in the COS-Halos sample, and confirm this already reveals strong evidence for an effective transport rate which rises rapidly away from the interstellar medium to values $\kappa_{\rm min}^{\rm eff}\gtrsim 10^{30-31}\,{\rm cm}^2\,{\rm s}^{-1}$ (corresponding to anisotropic streaming velocities of $v^{\rm stream}_{\rm eff} \gtrsim 1000\,{\rm km}\,{\rm s}^{-1}$) in the diffuse CGM, at impact parameters larger than $50-100$\,kpc. We discuss how future observations can provide qualitatively new constraints in our understanding of cosmic rays in the CGM and intergalactic medium., Comment: 9 pages, 2 figures, accepted to MNRAS

Partial dust obscuration in active galactic nuclei (AGN) has been proposed as a potential explanation for some cases of AGN variability. The dust-gas mixture present in AGN torii is accelerated by radiation pressure leading to the launching of an AGN wind. Dust under these conditions has been shown to be unstable to a generic class of fast growing resonant drag instabilities (RDIs). We present the first set of numerical simulations of radiation driven outflows that include explicit dust dynamics in conditions resembling AGN winds and discuss the implications of the RDIs on the morphology of the AGN torus, AGN variability, and the ability of the radiation to effectively launch a wind. We find that the RDIs rapidly develop reaching saturation at times much shorter than the global timescales of the outflows, resulting in the formation of filamentary structure on box-size scales with strong dust clumping and super-Alfv\'enic velocity dispersions on micro-scales. This results in 10-20% fluctuations in dust opacity and gas column density, integrated along mock observed lines-of-sight to the quasar accretion disk, over year to decade timescales with a red-noise power spectrum which is commonly observed for AGN. Additionally, all our simulations show that the radiation is sufficiently coupled to the dust-gas mixture launching highly super-sonic winds, which entrain 70-90% of gas, with a factor of $\lesssim 3$ photon momentum loss relative to the ideal case. Therefore, the RDIs could play an important role in generating the clumpy nature of AGN torii and driving AGN variability consistent with observations., Comment: 19 pages, 15 Figures

Several recent simulations of galaxy formation predict two main phases of supermassive black hole (BH) accretion: an early, highly intermittent phase (during which BHs are under-massive relative to local scaling relations), followed by a phase of accelerated growth. We investigate physical factors that drive the transition in BH accretion in cosmological zoom-in simulations from the FIRE project, ranging from dwarf galaxies to galaxies sufficiently massive to host luminous quasars. The simulations model multi-channel stellar feedback, but neglect AGN feedback. We show that multiple physical properties, including halo mass, galaxy stellar mass, and depth of the central gravitational potential correlate with accelerated BH fueling: constant thresholds in these properties are typically crossed within ~0.1 Hubble time of accelerated BH fueling. Black hole masses increase sharply when the stellar surface density in the inner 1 kpc crosses a threshold Sigma1 ~ 10^9.5 Msun/kpc^2, a characteristic value above which gravity prevents stellar feedback from ejecting gas, and similar to the value above which galaxies are observed to quench. We further show that accelerated BH growth correlates with the emergence of long-lived, thin gas disks, as well as with virialization of the inner circumgalactic medium. The halo mass Mh ~ 10^12 Msun and stellar mass Mstar ~ 10^10.5 Msun at which BH growth accelerates correspond to ~L* galaxies. The fact that stellar feedback becomes inefficient at ejecting gas from the nucleus above this mass scale may play an important role in explaining why AGN feedback appears to be most important in galaxies above ~L*., Comment: 21 pages, 9 figures, submitted to MNRAS