Search

We propose a novel dust battery mechanism for generating seed magnetic fields in the early universe, in which charged dust grains are radiatively accelerated, inducing strong electric currents that subsequently generate magnetic fields. Our analysis demonstrates that this process is effective even at very low metallicities (approximately $ \sim 10^{-5} Z_\odot$), and capable of producing seed fields with significant amplitudes of $B \sim \rm \mu G$ around luminous sources over timescales of years to Myr and across spatial scales ranging from au to kpc. Crucially, we find that this mechanism is generically $\sim10^8$ times more effective than the radiatively-driven electron battery or Biermann battery in relatively cool gas ($\ll 10^{5}\,$K), including both neutral and ionized gas. Furthermore, our results suggest that, to first order, dissipation effects do not appear to significantly impede this process, and that it can feasibly generate coherent seed fields on macroscopically large ISM scales (much larger than turbulent dissipation scales or electron mean-free-paths in the ISM). These seed fields could then be amplified by subsequent dynamo actions to the observed magnetic fields in galaxies. Additionally, we propose a sub-grid model for integration into cosmological simulations, and the required electric-field expressions for magnetohydrodynamic-particle-in-a-cell (MHD-PIC) simulations that explicitly model dust dynamics. Finally, we explore the broad applicability of this mechanism across different scales and conditions, emphasizing its robustness compared to other known battery mechanisms., Comment: 19 pages, 2 figures, 1 appendix. Submitted to ApJ

Coupling black hole (BH) feeding and feedback involves interactions across vast spatial and temporal scales that is computationally challenging. Tracking gas inflows and outflows from kilo-parsec scales to the event horizon for non-spinning BHs in the presence of strong magnetic fields, Cho et al. (2023, 2024) report strong suppression of accretion on horizon scales and low (2%) feedback efficiency. In this letter, we explore the impact of these findings for the supermassive BHs M87* and Sgr A*, using high-resolution, non-cosmological, magnetohydrodynamic (MHD) simulations with the Feedback In Realistic Environments (FIRE-2) model. With no feedback, we find rapid BH growth due to "cooling flows," and for 2% efficiency feedback, while accretion is suppressed, the rates still remain higher than constraints from Event Horizon Telescope (EHT) data (Event Horizon Telescope Collaboration et al. 2021, 2022) for M87* and Sgr A*. To match EHT observations of M87*, a feedback efficiency greater than 15% is required, suggesting the need to include enhanced feedback from BH spin. Similarly, a feedback efficiency of $>15\%$ is needed for Sgr A* to match the estimated observed star formation rate of $\lesssim 2 {\rm M_\odot}$ yr$^{-1}$. However, even with 100% feedback efficiency, the accretion rate onto Sgr A* matches with EHT data only on rare occasions in the simulations, suggesting that Sgr A* is likely in a temporary quiescent phase currently. Bridging accretion and feedback across scales, we conclude that higher feedback efficiency, possibly due to non-zero BH spin, is necessary to suppress "cooling flows" and match observed accretion and star formation rates in M87* and Sgr A*., Comment: 13 pages, 4 figures

$\Lambda$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 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 Milky Way-mass galaxies, we find the stellar density profiles along the minor axis (perpendicular to the galactic disk) within the NFW scale radii (R$\approx$15 kpc) are best described as an exponential disk with scale height <0.3 kpc and a power law component with slope $\alpha\approx$-4. The stellar density amplitude and slope for the power law component is 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^2$ along the minor axis., Comment: 29 pages, 13 figures, 1 table, Submitted to ApJ, comments welcome!

The far-infrared (FIR) - radio correlation (FRC) is one of the most promising empirical constraints on the role of cosmic-rays (CRs) and magnetic fields (\textbf{B}) in galaxy formation and evolution. While many theories have been proposed in order to explain the emergence and maintenance of the FRC across a gamut of galaxy properties and redshift, the non-linear physics at play remain unexplored in full complexity and cosmological context. We present the first reproduction of the $z \sim 0$ FRC using detailed synthetic observations of state-of-the-art cosmological zoom-in simulations from the FIRE-3 suite with explicitly-evolved CR proton and electron (CRe) spectra, for three models for CR transport and multi-channel AGN feedback. In doing so, we generally verify the predictions of `calorimeter' theories at high FIR luminosities (\Lsixty\, $\gtrsim$ 10$^{9.5}$) and at low FIR luminosities (\Lsixty\, $\lesssim$ 10$^{9.5}$) the so-called `conspiracy' of increasing ultraviolet radiation escape in tandem with increasing CRe escape, and find that the global FRC is insensitive to \textit{orders-of-magnitude} locally-variable CR transport coefficients. Importantly, the indirect effect of AGN feedback on emergent observables highlights novel interpretations of outliers in the FRC. In particular, we find that in many cases, `radio-excess' objects can be better understood as \textit{IR-dim} objects with longer-lived radio contributions at low $z$ from Type Ia SNe and intermittent black hole accretion in quenching galaxies, though this is sensitive to the interplay of CR transport and AGN feedback physics. This creates characteristic evolutionary tracks leading to the $z=0$ FRC, which shape the subsequent late-time behavior of each model., Comment: 21 pages, of which 16 are main text, with 6 Figures, 1 Table, and one additional Appendix Figure. The remaining 4 pages are references. Submitted to The Astrophysical Journal

Quasars are powered by supermassive black hole (SMBH) accretion disks, yet standard disk models are inconsistent with many quasar observations. Recently, Hopkins et al. (2024) simulated the formation of a quasar disk feeding a SMBH of mass $M=1.3\times10^7\,M_\odot$ in a host galaxy that evolved from cosmological initial conditions. The disk had surprisingly strong toroidal magnetic fields that supported it vertically from gravity and powered fast accretion. What radiation and feedback can such a system produce? To answer this, we must follow the gas to the event horizon. For this, we interpolated the accretion system onto the grid of the general-relativistic radiation magnetohydrodynamics code H-AMR and performed 3D simulations with BH spins $a=0$ and $a=0.9375$. This remapping generates spurious magnetic monopoles, which we erase using a novel divergence cleaning approach. Despite the toroidal magnetic field's dominance at large radii, vertical magnetic flux builds up at the event horizon. This causes a magnetic state transition within the inner $200$ gravitational radii of the disk, where net vertical magnetic flux begins dominating the accretion flow. This powers strong winds and, if the BH spins, relativistic jets that can spin-down the BH within $5-10\,{\rm Myrs}$. Sometimes, vertical magnetic fields of opposite sign reach the BH, causing polarity inversion events that briefly destroy the jets and, possibly, the X-ray corona. The disk powers accretion at rates $5\times$ the Eddington limit, which can double the BH mass in $5-10\,{\rm Myrs}$. When $a=0.9375$ ($a=0$), the energy in outflows and radiation equals about $60\%$ ($10\%$) and $100\%$ ($3\%$) of the accreted rest mass energy, respectively. Much of the light escapes in cool, extended $\gtrsim1300\,{\rm au}$ photospheres, consistent with quasar microlensing and the ``big blue bump'' seen in spectral energy distributions., Comment: 33 pages, 18 figures. Submitted to ApJ. Comments welcome!

We infer supermassive black hole (SMBH) accretion rates and Eddington ratios as a function of SMBH/host galaxy mass and redshift with the empirical TRINITY model of dark matter halo--galaxy--SMBH connection. The galaxy--SMBH mass and growth rate connection from TRINITY is constrained by galaxy observables from $06$: at these redshifts, dark matter halos experience $\sim$Eddington specific growth rates, driving $\sim$Eddington specific growth rates in both galaxies and SMBHs., Comment: 14 pages, 19 figures, submitted to MNRAS

Using the FIRE-2 cosmological zoom-in simulations, we investigate the temporal evolution of gas-phase metallicity radial gradients of Milky Way-mass progenitors in the redshift range of $0.4

We show that there exist two qualitatively different turbulent states of the zero-net-vertical-flux shearing box. The first, which has been studied in detail previously, is characterized by a weakly magnetized ($\beta\sim50$) midplane with slow periodic reversals of the mean azimuthal field (dynamo cycles). The second (the "low-$\beta$ state"), which is the main subject of this paper, is characterized by a strongly magnetized $\beta\sim1$ midplane dominated by a coherent azimuthal field with much stronger turbulence and much larger accretion stress $\alpha \sim 1$. The low-$\beta$ state is realized in simulations that begin with sufficiently strong azimuthal magnetic fields. The mean azimuthal field in the low-$\beta$ state is quasi steady (no cycles) and is sustained by a dynamo mechanism that compensates for the continued loss of magnetic flux through the vertical boundaries; we attribute the dynamo to the combination of differential rotation and the Parker instability, although many of its details remain unclear. Vertical force balance in the low-$\beta$ state is dominated by the mean magnetic pressure except at the midplane, where thermal pressure support is always important (this is true even when simulations are initialized at $\beta\ll1$, provided the thermal scale-height of the disk is well-resolved). The efficient angular momentum transport in the low-$\beta$ state may resolve long-standing tension between predictions of magnetorotational turbulence (at high $\beta$) and observations; likewise, the bifurcation in accretion states we find may be important for understanding the state transitions observed in dwarf novae, X-ray binaries, and changing-look AGN. We discuss directions for future work including the implications of our results for global accretion disk simulations.

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 star-forming environments in the Universe, but its impact 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 varying environmental dependencies. Processes such as grain coagulation, expected to be efficient in dense star-forming regions, reduce the abundance of small grains and increase that of larger grains. Motivated by this, we investigate the effects of similar GSD variations on the thermochemistry and evolution of giant molecular clouds (GMCs) using 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 \SF project. We find that grain size variations significantly alter GMC thermochemistry: with the leading-order effect is that larger grains, under fixed dust mass, GSD dynamic range, and dust-to-gas ratio, result in lower dust opacities. This reduced opacity permits ISRF and internal radiation photons to penetrate more deeply. This leads to rapid gas heating and inhibited star formation. Star formation efficiency is highly sensitive to grain size, with an order of magnitude reduction when grain size dynamic range increases from $10^{-3}$-0.1 $\rm\mu m$ to 0.1-10 $\rm\mu m$. Additionally, warmer gas suppresses low-mass star formation, and decreased opacities result in a greater proportion of gas 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 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. This removal of dust, which interacts with gas through cooling, chemistry, drag, and radiative transfer, alters the gas properties in the region. 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: Updated version; accepted to A&A

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: 39 pages (including Appendix), 25 figures, accepted in 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