Search

Simulations and observations suggest that galaxy interactions may enhance the star formation rate (SFR) in merging galaxies. One proposed mechanism is the torque exerted on the gas and stars in the larger galaxy by the smaller galaxy. We analyze the interaction torques and star formation activity on six galaxies from the FIRE-2 simulation suite with masses comparable to the Milky Way galaxy at redshift $z=0$. We trace the halos from $z = 3.6$ to $z=0$, calculating the torque exerted by the nearby galaxies on the gas in the central galaxy. We calculate the correlation between the torque and the SFR across the simulations for various mass ratios. For near-equal-stellar-mass-ratio interactions in the galaxy sample, occurring between $z=1.2-3.6$, there is a positive and statistically significant correlation between the torque from nearby galaxies on the gas of the central galaxies and the SFR. For all other samples, no statistically significant correlation is found between the torque and the SFR. Our analysis shows that some, but not all, major interactions cause starbursts in the simulated Milky Way-mass galaxies, and that most starbursts are not caused by galaxy interactions. The transition from `bursty' at high redshift ($z\gtrsim1$) to `steady' star-formation state at later times is independent of the interaction history of the galaxies, and most of the interactions do not leave significant imprints on the overall trend of the star formation history of the galaxies., Comment: Submitted to ApJ

The observed prevalence of galaxies exhibiting bursty star formation histories (SFHs) at $z\gtrsim6$ has created new challenges and opportunities for understanding their formation pathways. The degenerate effects of the efficiency and burstiness of star formation on the observed UV luminosity function are separable by galaxy clustering. However, quantifying the timescales of burstiness requires more than just the continuum UV measurements. Here we develop a flexible semi-analytic framework for modeling both the amplitude of star formation rate (SFR) variations and their temporal correlation, from which the luminosity function and clustering can be derived for SFR indicators tracing different characteristic timescales (e.g., UV continuum and H$\alpha$ luminosities). Based on this framework, we study the prospect of using galaxy summary statistics to distinguish models where SFR fluctuations are prescribed by different power spectral density (PSD) forms. Using the Fisher matrix approach, we forecast the constraints on parameters in our PSD-based model that can be extracted from mock JWST observations of the UV and H$\alpha$ luminosity functions and clustering bias factors at $z\sim6$. These constraints demonstrate the feasibility of constraining the burstiness of high-$z$ galaxies solely from their one-point and two-point statistics and underscore the importance of combining tracers of both long-term and short-term SFR variations. Our flexible framework can be readily extended to characterize the SFH of high-redshift galaxies with a wider range of observational diagnostics., Comment: 21 pages, 7 figures, submitted to JCAP

The circumgalactic medium (CGM) in $\gtrsim 10^{12}$ $\mathrm{M}_{\odot}$ halos is dominated by a hot phase ($T \gtrsim 10^{6}$ K). While many models exist for the hot gas structure, there is as yet no consensus. We compare cooling flow models, in which the hot CGM flows inward due to radiative cooling, to the CGM of $\sim 10^{12}-10^{13}$ $\mathrm{M}_{\odot}$ halos in galaxy formation simulations from the FIRE project at $z\sim0$. The simulations include realistic cosmological evolution and feedback from stars but neglect AGN feedback. At both mass scales, CGM inflows are typically dominated by the hot phase rather than by the `precipitation' of cold gas. Despite being highly idealized, we find that cooling flows describe $\sim 10^{13}$ $\mathrm{M}_{\odot}$ halos very well, with median agreement in the density and temperature profiles of $\sim 20\%$ and $\sim 10\%$, respectively. This indicates that stellar feedback has little impact on CGM scales in those halos. For $\sim 10^{12}$ $\mathrm{M}_{\odot}$ halos, the thermodynamic profiles are also accurately reproduced in the outer CGM. For some of these lower-mass halos, cooling flows significantly overpredict the hot gas density in the inner CGM. This could be due to multidimensional angular momentum effects not well captured by our 1D cooling flow models and/or to the larger cold gas fractions in these regions. Turbulence, which contributes $\sim 10-40\%$ of the total pressure, must be included to accurately reproduce the temperature profiles. Overall, cooling flows predict entropy profiles in better agreement with the FIRE simulations than other idealized models in the literature., Comment: 20 pages, 14 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

The radial gradient of gas-phase metallicity is a powerful probe of the chemical and structural evolution of star-forming galaxies, closely tied to disk formation and gas kinematics in the early universe. We present spatially resolved chemical and dynamical properties for a sample of 25 galaxies at $0.5 \lesssim z \lesssim 1.7$ from the \msasd survey. These innovative observations provide 3D spectroscopy of galaxies at a spatial resolution approaching JWST's diffraction limit and a high spectral resolution of $R\simeq2700$. The metallicity gradients measured in our galaxy sample range from $-$0.03 to 0.02 dex~kpc$^{-1}$. Most galaxies exhibit negative or flat radial gradients, indicating lower metallicity in the outskirts or uniform metallicity throughout the entire galaxy. We confirm a tight relationship between stellar mass and metallicity gradient at $z\sim1$ with small intrinsic scatter of 0.02 dex~kpc$^{-1}$. Our results indicate that metallicity gradients become increasingly negative as stellar mass increases, likely because the more massive galaxies tend to be more ``disky". This relationship is consistent with the predictions from cosmological hydrodynamic zoom-in simulations with strong stellar feedback. This work presents the effort to harness the multiplexing capability of JWST NIRSpec/MSA in slit-stepping mode to map the chemical and kinematic profiles of high-redshift galaxies in large samples and at high spatial and spectral resolution., Comment: 24 pages, 6 figures, 1 table; Accepted for publication in ApJL

The shallow potential wells of star-forming dwarf galaxies make their surrounding circumgalactic and intergalactic medium (CGM/IGM) sensitive laboratories for studying the inflows and outflows thought to regulate galaxy evolution. We present new absorption-line measurements in quasar sightlines probing within projected distances of $<300$ kpc from 91 star-forming field dwarf galaxies with a median stellar mass of $\log{M_\star/\rm{M_\odot}} \approx 8.3$ at $0.077 < z < 0.73$ from the Cosmic Ultraviolet Baryon Survey (CUBS). In this redshift range, the CUBS quasar spectra cover a suite of transitions including H I, low and intermediate metal ions (e.g., C II, Si II, C III, and Si III), and highly ionized O VI. This CUBS-Dwarfs survey enables constraints with samples 9$\times$ larger than past dwarf CGM/IGM studies with similar ionic coverage. We find that low and intermediate ionization metal absorption is rare around dwarf galaxies, consistent with previous surveys of local dwarfs. In contrast, highly ionized O VI is commonly observed in sightlines that pass within the virial radius of a dwarf, and O VI detection rates are non-negligible at projected distances of 1$-$2$\times$ the virial radius. Based on these measurements, we estimate that the O VI-bearing phase of the CGM/IGM accounts for a dominant share of the metal budget of dwarf galaxies. The absorption kinematics suggest that a relatively modest fraction of the O VI-bearing gas is formally unbound. Together, these results imply that low-mass systems at $z\lesssim 1$ effectively retain a substantial fraction of their metals within the nearby CGM and IGM., Comment: Accepted to ApJ. 18 pages, 6 figures, 2 tables

Integral field spectroscopy (IFS) is a powerful tool for understanding the formation of galaxies across cosmic history. We present the observing strategy and first results of MSA-3D, a novel JWST program using multi-object spectroscopy in a slit-stepping strategy to produce IFS data cubes. The program observed 43 normal star-forming galaxies at redshifts $0.5 \lesssim z \lesssim 1.5$, corresponding to the epoch when spiral thin-disk galaxies of the modern Hubble sequence are thought to emerge, obtaining kpc-scale maps of rest-frame optical nebular emission lines with spectral resolution $R\simeq2700$. Here we describe the multiplexed slit-stepping method which is $>15$ times more efficient than the NIRSpec IFS mode for our program. As an example of the data quality, we present a case study of an individual galaxy at $z=1.104$ (stellar mass $M_{*} = 10^{10.3}~M_{\odot}$, star formation rate~$=3~M_{\odot}$ yr$^{-1}$) with prominent face-on spiral structure. We show that the galaxy exhibits a rotationally supported disk with moderate velocity dispersion ($\sigma = 36^{+5}_{-4}$~\kms), a negative radial metallicity gradient ($-0.020\pm0.002$~dex\,kpc$^{-1}$), a dust attenuation gradient, and an exponential star formation rate density profile which closely matches the stellar continuum. These properties are characteristic of local spirals, indicating that mature galaxies are in place at $z\sim1$. We also describe the customized data reduction and original cube-building software pipelines which we have developed to exploit the powerful slit-stepping technique. Our results demonstrate the ability of JWST slit-stepping to study galaxy populations at intermediate to high redshifts, with data quality similar to current surveys of the $z\sim0.1$ universe., Comment: 24 pages, 12 figures, 1 table; submitted to ApJ, comments welcome. Our custom designed data reduction pipeline (MSA3D) and example processed data-cube will be publicly released before the JWST Cycle 4 deadline. MSA3D package on GitHub: https://github.com/barisiciv/msa3d.git

Recent observations with the James Webb Space Telescope (JWST) have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the Big Bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox-HR, a high-resolution, cosmological hydrodynamical simulation from the Feedback in Realistic Environments project, which offers insights into the SFE of galaxies during the first billion years of cosmic time. FIREbox-HR re-simulates the cosmic volume (L = 22.1 cMpc) of the original FIREbox run with eight times higher mass resolution (m_b ~ 7800 M_sun), but with identical physics, down to z ~ 6. FIREbox-HR predicts ultraviolet (UV) luminosity functions in good agreement with available observational data. The simulation also successfully reproduces the observed cosmic UV luminosity density at z ~ 6 - 14, demonstrating that relatively high star formation activity in the early Universe is a natural outcome of the baryonic processes encoded in the FIRE-2 model. According to FIREbox-HR, the SFE - halo mass relation for intermediate mass halos (M_halo ~ 10^9 - 10^11 M_sun) does not significantly evolve with redshift and is only weakly mass-dependent. These properties of the SFE - halo mass relation lead to a larger contribution from lower mass halos at higher z, driving the gradual evolution of the observed cosmic UV luminosity density. A theoretical model based on the SFE - halo mass relation inferred from FIREbox-HR allows us to explore implications for galaxy evolution. Future observations of UV faint galaxies at z > 12 will provide an opportunity to further test these predictions and deepen our understanding of star formation during Cosmic Dawn., Comment: 28 pages, 14 figures, 5 tables, revised to match version accepted by MNRAS, Appendix D added, UV LF shown in Fig. 1 provided as txt file

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.

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

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

Our ability to trace the star-forming molecular gas is important to our understanding of the Universe. We can trace this gas using CO emission, converting the observed CO intensity into the H$_2$ gas mass of the region using the CO-to-H$_2$ conversion factor (Xco). In this paper, we use simulations to study the conversion factor and the molecular gas within galaxies. We analysed a suite of simulations of isolated disc galaxies, ranging from dwarfs to Milky Way-mass galaxies, that were run using the FIRE-2 subgrid models coupled to the CHIMES non-equilibrium chemistry solver. We use the non-equilibrium abundances from the simulations, and we also compare to results using abundances assuming equilibrium, which we calculate from the simulation in post-processing. Our non-equilibrium simulations are able to reproduce the relation between CO and H$_2$ column densities, and the relation between Xco and metallicity, seen within observations of the Milky Way. We also compare to the xCOLD GASS survey, and find agreement with their data to our predicted CO luminosities at fixed star formation rate. We also find the multivariate function used by xCOLD GASS overpredicts the H$_2$ mass for our simulations, motivating us to suggest an alternative multivariate function of our fitting, though we caution that this fitting is uncertain due to the limited range of galaxy conditions covered by our simulations. We also find that the non-equilibrium chemistry has little effect on the conversion factor (<5\%) for our high-mass galaxies, though still affects the H$_2$ mass and Lco by $\approx$25\%., Comment: 19 pages, 16 figures (including appendices). Accepted for publication in MNRAS

The unprecedented infrared spectroscopic capabilities of JWST have provided high-quality interstellar medium (ISM) metallicity measurements and enabled characterization of the gas-phase mass-metallicity relation (MZR) for galaxies at $z \gtrsim 5$ for the first time. We analyze the gas-phase MZR and its evolution in a high-redshift suite of FIRE-2 cosmological zoom-in simulations at $z=5-12$ and for stellar masses $M_* \sim 10^6-10^{10} \rm{M}_\odot$. These simulations implement a multi-channel stellar feedback model and produce broadly realistic galaxy properties, including when evolved to $z=0$. The simulations predict very weak redshift evolution of the MZR over the redshift range studied, with the normalization of the MZR increasing by less than $0.01$ dex as redshift decreases from $z = 12$ to $z=5$. The median MZR in the simulations is well-approximated as a constant power-law relation across this redshift range given by $\log(Z/Z_\odot) = 0.37\log(M_*/\rm{M}_\odot) - 4.3$. We find good agreement between our best-fit model and recent observations made by JWST at high redshift. The weak evolution of the MZR at $z > 5$ contrasts with the evolution at $z \lesssim 3$, where increasing normalization of the MZR with decreasing redshift is observed and predicted by most models. The FIRE-2 simulations predict increasing scatter in the gas-phase MZR with decreasing stellar mass, in qualitative agreement with some observations., Comment: 13 pages, 6 figures

This paper presents a newly established sample of 103 unique galaxies or galaxy groups at $0.4\lesssim z\lesssim 0.7$ from the Cosmic Ultraviolet Baryon Survey (CUBS) for studying the warm-hot circumgalactic medium (CGM) probed by both O VI and Ne VIII absorption. The galaxies and associated neighbors are identified at $< 1$ physical Mpc from the sightlines toward 15 CUBS QSOs at $z_{\rm QSO}\gtrsim 0.8$. A total of 30 galaxies or galaxy groups exhibit associated O VI $\lambda\lambda$ 1031, 1037 doublet absorption within a line-of-sight velocity interval of $\pm250$ km/s, while the rest show no trace of O VI to a detection limit of $\log N_{\rm OVI}/{\rm cm^{-2}}\approx13.7$. Meanwhile, only five galaxies or galaxy groups exhibit the Ne VIII $\lambda\lambda$ 770,780 doublet absorption, down to a limiting column density of $\log N_{\rm NeVIII}/{\rm cm^{-2}}\approx14.0$. These O VI- and Ne VIII-bearing halos reside in different galaxy environments with stellar masses ranging from $\log M_{\rm star}/M_\odot \approx 8$ to $\approx11.5$. The warm-hot CGM around galaxies of different stellar masses and star formation rates exhibits different spatial profiles and kinematics. In particular, star-forming galaxies with $\log M_{\rm star}/M_\odot\approx9-11$ show a significant concentration of metal-enriched warm-hot CGM within the virial radius, while massive quiescent galaxies exhibit flatter radial profiles of both column densities and covering fractions. In addition, the velocity dispersion of O VI absorption is broad with $\sigma_v > 40$ km/s for galaxies of $\log M_{\rm star}/M_\odot>9$ within the virial radius, suggesting a more dynamic warm-hot halo around these galaxies. Finally, the warm-hot CGM probed by O VI and Ne VIII is suggested to be the dominant phase in sub-$L^*$ galaxies with $\log M_{\rm star}/M_\odot\approx9-10$ based on their high ionization fractions in the CGM., Comment: Submitted to ApJ after addressing the referee's comments; 28 pages, 16 figures

Observations show a tight correlation between the stellar mass of galaxies and their gas-phase metallicity (MZR). This relation evolves with redshift, with higher-redshift galaxies being characterized by lower metallicities. Understanding the physical origin of the slope and redshift evolution of the MZR may provide important insight into the physical processes underpinning it: star formation, feedback, and cosmological inflows. While theoretical models ascribe the shape of the MZR to the lower efficiency of galactic outflows in more massive galaxies, what drives its evolution remains an open question. In this letter, we analyze how the MZR evolves over $z=0-3$, combining results from the FIREbox cosmological volume simulation with analytical models. Contrary to a frequent assertion in the literature, we find that the evolution of the gas fraction does not contribute significantly to the redshift evolution of the MZR. Instead, we show that the latter is driven by the redshift-dependence of the inflow metallicity, outflow metallicity, and mass loading factor, whose relative importance depends on stellar mass. These findings also suggest that the evolution of the MZR is not explained by galaxies moving along a fixed surface in the space spanned by stellar mass, gas phase metallicity, and star formation rate., Comment: 6 pages, 4 figures. Submitted for publication on MNRAS

The properties of warm-hot gas around $\sim L_{*}$ galaxies can be studied with absorption lines from highly ionized metals. We predict Ne VIII column densities from cosmological zoom-in simulations of halos with masses in $\sim 10^{12}$ and $\sim 10^{13}\,\mathrm{M}_{\odot}$ from the FIRE project. Ne VIII traces the volume-filling, virial-temperature gas in $\sim 10^{12}\,\mathrm{M}_{\odot}$ halos. In $\sim 10^{13}\,\mathrm{M}_{\odot}$ halos the Ne VIII gas is clumpier, and biased towards the cooler part of the warm-hot phase. We compare the simulations to observations by the CASBaH and CUBS surveys. We show that when inferring halo masses from stellar masses to compare simulated and observed halos, it is important to account for the scatter in the stellar-mass-halo-mass relation, especially at $\,\mathrm{M}_{\star} \gtrsim 10^{10.5} \,\mathrm{M}_{\odot}$. Median Ne VIII columns in the fiducial FIRE-2 model are about as high as observed upper limits allow, while the simulations analyzed do not reproduce the highest observed columns. This suggests that the median Ne VIII profiles predicted by the simulations are consistent with observations, but that the simulations may underpredict the scatter. We find similar agreement with analytical models that assume a product of the halo gas fraction and metallicity (relative to solar) $\sim 0.1$-$0.3$, indicating that observations are consistent with plausible CGM temperatures, metallicities, and gas masses. Variants of the FIRE simulations with a modified supernova feedback model and/or AGN feedback included (as well as some other cosmological simulations from the literature) more systematically underpredict Ne VIII columns. The circumgalactic Ne VIII observations therefore provide valuable constraints on simulations that otherwise predict realistic galaxy properties., Comment: version accepted by ApJ

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

Atomic hydrogen (HI) serves a crucial role in connecting galactic-scale properties such as star formation with the large-scale structure of the Universe. While recent numerical simulations have successfully matched the observed covering fraction of HI near Lyman Break Galaxies (LBGs) and in the foreground of luminous quasars at redshifts $z \lesssim 3$, the low-mass end remains as-of-yet unexplored in observational and computational surveys. We employ a cosmological, hydrodynamical simulation (FIREbox) supplemented with zoom-in simulations (MassiveFIRE) from the Feedback In Realistic Environments (FIRE) project to investigate the HI covering fraction of Lyman Limit Systems ($N_{\mathrm{HI}} \gtrsim 10^{17.2}$ cm$^{-2}$) across a wide range of redshifts ($z=0-6$) and halo masses ($10^8-10^{13} M_{\odot}$ at $z=0$, $10^8-10^{11} M_{\odot}$ at $z=6$) in the absence of feedback from active galactic nuclei. We find that the covering fraction inside haloes exhibits a strong increase with redshift, with only a weak dependence on halo mass for higher-mass haloes. For massive haloes ($M_{\mathrm{vir}} \sim 10^{11}-10^{12} M_{\odot}$), the radial profiles showcase scale-invariance and remain independent of mass. The radial dependence is well-captured by a fitting function. The covering fractions in our simulations are in good agreement with measurements of the covering fraction in LBGs. Our comprehensive analysis unveils a complex dependence with redshift and halo mass for haloes with $M_{\mathrm{vir}} \lesssim 10^{10} M_{\odot}$ that future observations aim to constrain, providing key insights into the physics of structure formation and gas assembly.

Lightcurves of many classical novae deviate from the canonical "fast rise - smooth decline" pattern and display complex variability behavior. We present the first TESS-space-photometry-based investigation of this phenomenon. We use TESS Sector 41 full-frame images to extract a lightcurve of the slow Galactic nova V606 Vul that erupted nine days prior to the start of the TESS observations. The lightcurve covers the first of two major peaks of V606 Vul that was reached 19 days after the start of the eruption. The nova reached its brightest visual magnitude V=9.9 in its second peak 64 days after the eruption onset, following the completion of Sector 41 observations. To increase the confidence level of the extracted lightcurve, we performed the analysis using four different codes implementing the aperture photometry (Lightkurve, VaST) and image subtraction (TESSreduce, tequila_shots) and find good agreement between them. We performed ground-based photometric and spectroscopic monitoring to complement the TESS data. The TESS lightcurve reveals two features: periodic variations (0.12771d, 0.01mag average peak-to-peak amplitude) that disappeared when the source was within 1mag of peak optical brightness and a series of isolated mini-flares (with peak-to-peak amplitudes of up to 0.5mag) appearing at seemingly random times. We interpret the periodic variations as the result of azimuthal asymmetry of the photosphere engulfing the nova-hosting binary that was distorted by and rotating with the binary. Whereas we use spectra to associate the two major peaks in the nova lightcurve with distinct episodes of mass ejection, the origin of mini-flares remains elusive., Comment: 31 pages, 10 figures, submitted to ApJ; comments welcome

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

The origin of the cold phase in the CGM is a highly debated question. We investigate the contribution of satellite galaxies to the cold gas budget in the circumgalactic medium (CGM)of a Milky Way-like host galaxy. We perform controlled experiments with three different satellite mass distributions and identify several mechanisms by which satellites can add cold gas to the CGM, including ram pressure stripping and induced cooling in the mixing layer of the stripped cold gas. These two mechanisms contribute a comparable amount of cold gas to the host CGM. We find that the less massive satellites ($\leq 10^9 M_\odot$) not only lose all of their cold gas in a short period ($\sim$ 0.5-1 Gyr), but their stripped cold clouds also mix with the hot CGM gas and get heated up quickly. However, stellar feedback from these less massive satellites can hugely alter the fate of their stripped gas. Feedback speeds up the destruction of the stripped cold clouds from these satellites by making them more diffuse with more surface area. On the other hand, the more massive satellites (LMC or SMC-like $\sim 10^{10} M_\odot$) can add cold gas to the total gas budget of the host CGM for several Gyrs., Comment: 17 pages, 12 figures, 2 tables Accepted for publication in MNRAS on Oct 6, 2023

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

The shape of the low-mass (faint) end of the galaxy stellar mass function (SMF) or ultraviolet luminosity function (UVLF) at z > 6 is an open question for understanding which galaxies primarily drove cosmic reionisation. Resolved photometry of Local Group low-mass galaxies allows us to reconstruct their star formation histories, stellar masses, and UV luminosities at early times, and this fossil record provides a powerful `near-far' technique for studying the reionisation-era SMF/UVLF, probing orders of magnitude lower in mass than direct HST/JWST observations. Using 882 low-mass (Mstar < 10^9 Msun) galaxies across 11 Milky Way- and Local Group-analogue environments from the FIRE-2 cosmological baryonic zoom-in simulations, we characterise their progenitors at z ~ 6 - 9, the mergers/disruption of those progenitors over time, and how well their present-day fossil record traces the high-redshift SMF. A present-day galaxy with Mstar ~ 10^5 Msun (10^9 Msun) had ~1 (~30) progenitors at z ~ 7, and its main progenitor comprised ~100% (~50%) of the total stellar mass of all its progenitors at z ~ 7. We show that although only ~ 15% of the early population of low-mass galaxies survives to present day, the fossil record of surviving Local Group galaxies accurately traces the low-mass slope of the SMF at z ~ 6 - 9. We find no obvious mass dependence to the mergers and accretion, and show that applying this reconstruction technique to just the low-mass galaxies at z = 0 and not the MW/M31 hosts correctly recovers the slope of the SMF down to Mstar ~ 10^4.5 Msun at z > 6. Thus, we validate the `near-far' approach as an unbiased tool for probing low-mass reionisation-era galaxies., Comment: 16 pages, 6 figures, 3 appendix figures, and 1 table. Submitted to MNRAS

This paper reports the first measurement of the relationship between turbulent velocity and cloud size in the diffuse circumgalactic medium (CGM) in typical galaxy halos at redshift z~0.4-1. Through spectrally-resolved absorption profiles of a suite of ionic transitions paired with careful ionization analyses of individual components, cool clumps of size as small as l_cl~1 pc and density lower than nH = 0.001 cm^-3 are identified in galaxy halos. In addition, comparing the line widths between different elements for kinematically matched components provides robust empirical constraints on the thermal temperature T and the non-thermal motions bNT, independent of the ionization models. On average, bNT is found to increase with l_cl following bNT \propto l_cl^0.3 over three decades in spatial scale from l_cl~1 pc to l_cl~1 kpc. Attributing the observed bNT to turbulent motions internal to the clumps, the best-fit bNT-l_cl relation shows that the turbulence is consistent with Kolmogorov at <1 kpc with a roughly constant energy transfer rate per unit mass of epsilon~0.003 cm^2 s^-3 and a dissipation time scale of <~ 100 Myr. No significant difference is found between massive quiescent and star-forming halos in the sample on scales less than 1 kpc. While the inferred epsilon is comparable to what is found in CIV absorbers at high redshift, it is considerably smaller than observed in star-forming gas or in extended line-emitting nebulae around distant quasars. A brief discussion of possible sources to drive the observed turbulence in the cool CGM is presented., Comment: Accepted for publication in ApJ Letters; 10 pages, 4 figures

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)

Observational studies are finding stars believed to be relics of the earliest stages of hierarchical mass assembly of the Milky Way (i.e., proto-Galaxy). In this work, we contextualize these findings by studying the masses, ages, spatial distributions, morphology, kinematics, and chemical compositions of proto-galaxy populations from the 13 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. Our findings indicate that proto-Milky Way populations: i) can have a stellar mass range between $1\times10^{8}<\mathrm{M}_{\star}<2\times10^{10}[\mathrm{M}_{\odot}]$, a virial mass range between $3\times10^{10}<\mathrm{M}_{\star}<6\times10^{11}[\mathrm{M}_{\odot}]$, and be as young as $8 \lesssim \mathrm{Age} \lesssim 12.8$ [Gyr] ($1\lesssim z \lesssim 6$); ii) are predominantly centrally concentrated, with $\sim50\%$ of the stars contained within $5-10$ kpc; iii) on average show weak but systematic net rotation in the plane of the host's disc at $z=0$ (i.e., $0.25\lesssim\langle\kappa/\kappa_{\mathrm{disc}}\rangle\lesssim0.8$); iv) present [$\alpha$/Fe]-[Fe/H] compositions that overlap with the metal-poor tail of the host's old disc; v) tend to assemble slightly earlier in Local Group-like environments than in systems in isolation. Interestingly, we find that ~60% of the proto-Milky Way galaxies are comprised by 1 dominant system ($1/5\lesssim$M$_{\star}$/M$_{\star,\mathrm{proto-Milky Way}}$$\lesssim4/5$) and 4-5 lower mass systems (M$_{\star}$/M$_{\star,\mathrm{proto-Milky Way}}$$\lesssim1/10$); the other ~40% are comprised by 2 dominant systems and 3-4 lower mass systems. These massive/dominant proto-Milky Way fragments can be distinguished from the lower mass ones in chemical-kinematic samples, but appear (qualitatively) indistinguishable from one another. Our results could help observational studies disentangle if the Milky Way formed from one or two dominant systems., Comment: Accepted in MNRAS. 16 pages, 12 figures and 3 tables

Recent discoveries of a significant population of bright galaxies at cosmic dawn $\left(z \gtrsim 10\right)$ have enabled critical tests of cosmological galaxy formation models. In particular, the bright end of the galaxy UV luminosity function (UVLF) appears higher than predicted by many models. Using approximately 25,000 galaxy snapshots at $8 \leq z \leq 12$ in a suite of FIRE-2 cosmological "zoom-in'' simulations from the Feedback in Realistic Environments (FIRE) project, we show that the observed abundance of UV-bright galaxies at cosmic dawn is reproduced in these simulations with a multi-channel implementation of standard stellar feedback processes, without any fine-tuning. Notably, we find no need to invoke previously suggested modifications such as a non-standard cosmology, a top-heavy stellar initial mass function, or a strongly enhanced star formation efficiency. We contrast the UVLFs predicted by bursty star formation in these original simulations to those derived from star formation histories (SFHs) smoothed over prescribed timescales (e.g., 100 Myr). The comparison demonstrates that the strongly time-variable SFHs predicted by the FIRE simulations play a key role in correctly reproducing the observed, bright-end UVLFs at cosmic dawn: the bursty SFHs induce order-or-magnitude changes in the abundance of UV-bright ($M_\mathrm{UV} \lesssim -20$) galaxies at $z \gtrsim 10$. The predicted bright-end UVLFs are consistent with both the spectroscopically confirmed population and the photometrically selected candidates. We also find good agreement between the predicted and observationally inferred integrated UV luminosity densities, which evolve more weakly with redshift in FIRE than suggested by some other models., Comment: 12 pages, 4 figures + 1 table, ApJL accepted

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

This paper presents a new sample of 19 unique galaxies and galaxy groups at $z\approx1$ from the CUBS program, which is designated as the CUBSz1 sample. In this CUBSz1 sample, nine galaxies or galaxy groups show absorption features, while ten systems do not have detectable absorption with 2-$\sigma$ upper limits of log$N$(HeI)/cm$^{-2}\lesssim 13.5$ and log$N$(OV)/cm$^{-2}\lesssim 13.3$. Environmental properties of the galaxies, including galaxy overdensities, the total stellar mass and gravitational potential summed over all nearby neighbors, and the presence of local ionizing sources, are found to have a significant impact on the observed CGM absorption properties. Specifically, massive galaxies and galaxies in overdense regions exhibit a higher rate of incidence of absorption. At the same time, the observed CGM absorption properties in galaxy groups appear to be driven by the galaxy closest to the QSO sightline, rather than by the most massive galaxy or by mass-weighted properties. We introduce a total projected gravitational potential $\psi$, defined as $-\psi/G =\sum M_{{\rm halo}}/d_{{\rm proj}}$ summed over all group members, to characterize the overall galaxy environment. This projected gravitational potential correlates linearly with the maximum density detected in each sightline, consistent with higher-pressure gas being confined in deeper gravitational potential wells. In addition, we find that the radial profile of cool gas density exhibits a general decline from the inner regions to the outskirts, being in pressure balance with the hot halo. Finally, we note that the ionizing flux from nearby galaxies can generate an elevated $N$(HI)/$N$(HeI) ratio, which in turn provides a unique diagnostic of possible local sources contributing to the ionizing radiation field., Comment: 17 pages, 11 figures and 1 table in the main text; 40 figures and 1 table in the appendix; accepted for publication in MNRAS

Observed accretion rates onto the Milky-Way and other local spirals fall short of that required to sustain star formation for cosmological timescales. A potential avenue for this unseen accretion is an inflow in the volume-filling hot phase ($\sim10^6$ K) of the circumgalactic medium (CGM), as suggested by some cosmological simulations. Using hydrodynamic simulations and a new analytic solution valid in the slow-rotation limit, we show that a hot inflow spins up as it approaches the galaxy, while remaining hot, subsonic and quasi-spherical. At the radius of angular momentum support ($\approx15$ kpc for the Milky-Way) the hot flow flattens into a disk geometry and then cools from $\sim10^6$ K to $\sim10^4$ K at the disk-halo interface. Cooling affects all hot gas, rather than just a subset of individual gas clouds, implying that accretion via hot inflows does not rely on local thermal instability in contrast with 'precipitation' models for galaxy accretion. Prior to cooling and accretion the inflow completes $\sim t_{\rm cool}/t_{\rm ff}$ radians of rotation, where $t_{\rm cool}/t_{\rm ff}$ is the cooling time to free-fall time ratio in hot gas immediately outside the galaxy. The ratio $t_{\rm cool}/t_{\rm ff}$ may thus govern the development of turbulence and enhancement of magnetic fields in gas accreting onto low-redshift spirals. We argue that accretion via hot inflows can explain the observed truncation of nearby thin stellar disks at $\approx4$ disk radii. We also show that if rotating hot inflows are common in Milky-Way size disk galaxies, as predicted, then signatures should be observable with X-ray telescopes and FRB surveys., Comment: 21 pages, 12 figures, resubmitted to MNRAS following referee report

Understanding the star formation rate (SFR) variability and how it depends on physical properties of galaxies is important for developing and testing the theory of galaxy formation. We investigate how statistical measurements of the extragalactic background light (EBL) can shed light on this topic and complement traditional methods based on observations of individual galaxies. Using semi-empirical models of galaxy evolution and SFR indicators sensitive to different star formation timescales (e.g., H$\alpha$ and UV continuum luminosities), we show that the SFR variability, quantified by the joint probability distribution of the SFR indicators (i.e., the bivariate conditional luminosity function), can be characterized as a function of galaxy mass and redshift through the cross-correlation between deep, near-infrared maps of the EBL and galaxy distributions. As an example, we consider combining upcoming SPHEREx maps of the EBL with galaxy samples from Rubin/LSST. We demonstrate that their cross-correlation over a sky fraction of $f_\mathrm{sky}\sim0.5$ can constrain the joint SFR indicator distribution at high significance up to $z\sim2.5$ for mass-complete samples of galaxies down to $M_{*}\sim10^9\,M_{\odot}$. These constraints not only allow models of different SFR variability to be distinguished, but also provide unique opportunities to investigate physical mechanisms that require large number statistics such as environmental effects. The cross-correlations investigated illustrate the power of combining cosmological surveys to extract information inaccessible from each data set alone, while the large galaxy populations probed capture ensemble-averaged properties beyond the reach of targeted observations towards individual galaxies., Comment: 12 pages, 7 figures, MNRAS accepted

Both observations and simulations have shown strong evidence for highly time-variable star formation in low-mass and/or high-redshift galaxies, which has important observational implications because high-redshift galaxy samples are rest-UV selected and therefore particularly sensitive to the recent star formation. Using a suite of cosmological "zoom-in" simulations at $z>5$ from the Feedback in Realistic Environments (FIRE) project, we examine the implications of bursty star formation histories for observations of high-redshift galaxies with JWST. We characterize how the galaxy observability depends on the star formation history. We also investigate selection effects due to bursty star formation on the physical properties measured, such as the gas fraction, specific star formation rate, and metallicity. We find the observability to be highly time-dependent for galaxies near the survey's limiting flux due to the SFR variability: as the star formation rate fluctuates, the same galaxy oscillates in and out of the observable sample. The observable fraction $f_\mathrm{obs} = 50\%$ at $z \sim 7$ and $M_{\star} \sim 10^{8.5}$ to $10^{9}\,M_{\odot}$ for a JWST/NIRCam survey reaching a limiting magnitude of $m^\mathrm{lim}_\mathrm{AB} \sim 29$-$30$, representative of surveys such as JADES and CEERS. JWST-detectable galaxies near the survey limit tend to have properties characteristic of galaxies in the bursty phase: on average, they show approximately 2.5 times higher cold, dense gas fractions and 20 times higher specific star formation rates at a given stellar mass than galaxies below the rest-UV detection threshold. Our study represents a first step in quantifying selection effects and the associated biases due to bursty star formation in studying high-redshift galaxy properties., Comment: 8 pages, 4 figures, resubmitted after incorporating referee's comments; analysis expanded to include more galaxies and some quantitative results corrected

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

Spurred by rich, multi-wavelength observations and enabled by new simulations, ranging from cosmological to sub-pc scales, the last decade has seen major theoretical progress in our understanding of the circumgalactic medium. We review key physical processes in the CGM. Our conclusions include: (1) The properties of the CGM depend on a competition between gravity-driven infall and gas cooling. When cooling is slow relative to free fall, the gas is hot (roughly virial temperature) whereas the gas is cold (T~10^4 K) when cooling is rapid. (2) Gas inflows and outflows play crucial roles, as does the cosmological environment. Large-scale structure collimates cold streams and provides angular momentum. Satellite galaxies contribute to the CGM through winds and gas stripping. (3) In multiphase gas, the hot and cold phases continuously exchange mass, energy and momentum. The interaction between turbulent mixing and radiative cooling is critical. A broad spectrum of cold gas structures, going down to sub-pc scales, arises from fragmentation, coagulation, and condensation onto gas clouds. (4) Magnetic fields, thermal conduction and cosmic rays can substantially modify how the cold and hot phases interact, although microphysical uncertainties are presently large. Key open questions for future work include the mutual interplay between small-scale structure and large-scale dynamics, and how the CGM affects the evolution of galaxies., Comment: 69 pages, 13 figures. Accepted for publication in Annual Review of Astronomy and Astrophysics. Authors' draft. Edited version will appear in the next volume

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

Milky Way-mass galaxies in the FIRE-2 simulations demonstrate two main modes of star formation. At high redshifts star formation occurs in a series of short and intense bursts, while at low redshifts star formation proceeds at a steady rate with a transition from one mode to another at times ranging from 3 to 7 Gyr ago for different galaxies. We analyse how the mode of star formation affects iron and alpha-element abundance. We find that the early bursty regime imprints a measurable pattern in stellar elemental abundances in the form of a "sideways chevron" shape on the [Fe/H] - [O/Fe] plane and the scatter in [O/Fe] at a given stellar age is higher than when a galaxy is in the steady regime. That suggests that the evolution of [O/Fe] scatter with age provides an estimate of the end of the bursty phase. We investigate the feasibility of observing of this effect by adding mock observational errors to a simulated stellar survey and find that the transition between the bursty and steady phase should be detectable in the Milky Way, although larger observational uncertainties make the transition shallower. We apply our method to observations of the Milky Way from the Second APOKASC Catalog and estimate that the transition to steady star formation in the Milky Way happened 7-8 Gyrs ago, earlier than transition times measured in the simulations., Comment: 15 pages; Accepted 17 January 2023 in MNRAS

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

Low-mass galaxies are highly susceptible to environmental effects that can efficiently quench star formation. We explore the role of ram pressure in quenching low-mass galaxies ($M_{*}\sim10^{5-9}\,\rm{M}_{\odot}$) within 2 Mpc of Milky Way (MW) hosts using the FIRE-2 simulations. Ram pressure is highly variable across different environments, within individual MW haloes, and for individual low-mass galaxies over time. The impulsiveness of ram pressure -- the maximum ram pressure scaled to the integrated ram pressure prior to quenching -- correlates with whether a galaxy is quiescent or star-forming. The time-scale between maximum ram pressure and quenching is anticorrelated with impulsiveness, such that high impulsiveness corresponds to quenching time-scales $<1$ Gyr. Galaxies in low-mass groups ($M_\mathrm{*,host}\sim10^{7-9}\,\rm{M}_{\odot}$) outside of MW haloes experience typical ram pressure only slightly lower than ram pressure on MW satellites, helping to explain effective quenching via group pre-processing. Ram pressure on MW satellites rises sharply with decreasing distance to the host, and, at a fixed physical distance, more recent pericentre passages are typically associated with higher ram pressure because of greater gas density in the inner host halo at late times. Furthermore, the ram pressure and gas density in the inner regions of Local Group-like paired host haloes is higher at small angles off the host galaxy disc compared to isolated hosts. The quiescent fraction of satellites within these low-latitude regions is also elevated in the simulations and observations, signaling possible anisotropic quenching via ram pressure around MW-mass hosts., Comment: 13 pages, 10 figures, + appendices. Accepted for publication in MNRAS

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

Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated halos. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution function only require resolutions sufficient to keep the jet-inflated cavities stable. Comparing different model variations, as well as an independent jet model using a different hydrodynamics code, we show that the dominant uncertainties are the choices of jet properties within a given model. Independent of implementation, we find that light, thermal jets with low momentum flux tend to delay the onset of a cooling flow more efficiently on a $50$ Myr timescale than heavy, kinetic jets. The delay of the cooling flow originates from a displacement and boost in entropy of the central gas. If the jet kinetic luminosity depends on accretion rate, collimated, light, hydrodynamic jets are able to reduce cooling flows in halos, without a need for jet precession or wide opening angles. Comparing the jet feedback with a `kinetic wind' implementation shows that equal amounts of star formation rate reduction can be achieved by different interactions with the halo gas: the jet has a larger effect on the hot halo gas while leaving the denser, star forming phase in place, while the wind acts more locally on the star forming phase, which manifests itself in different time-variability properties., Comment: 22 pages, 20 figures, accepted for publication in MNRAS, comments welcome

The circumgalactic medium (CGM) plays a pivotal role in regulating gas flows around galaxies and thus shapes their evolution. However, the details of how galaxies and their CGM co-evolve remain poorly understood. We present a new time-dependent two-zone model that self-consistently tracks not just mass and metal flows between galaxies and their CGM but also the evolution of the global thermal and turbulent kinetic energy of the CGM. Our model accounts for heating and turbulence driven by both supernova winds and cosmic accretion as well as radiative cooling, turbulence dissipation, and halo outflows due to CGM overpressurization. We demonstrate that, depending on parameters, the CGM can undergo a phase transition (``thermalization'') from a cool, turbulence-supported phase to a virial-temperature, thermally-supported phase. This CGM phase transition is largely determined by the ability of radiative cooling to balance heating from supernova winds and turbulence dissipation. We perform an initial calibration of our model to the FIRE-2 cosmological hydrodynamical simulations and show that it can approximately reproduce the baryon cycles of the simulated halos. In particular, we find that, for these parameters, the phase transition occurs at high-redshift in ultrafaint progenitors and at low redshift in classical $M_{\rm vir}\sim10^{11}M_{\odot}$ dwarfs, while Milky Way-mass halos undergo the transition at $z\approx0.5$. We see a similar transition in the simulations though it is more gradual, likely reflecting radial dependence and multi-phase gas not captured by our model. We discuss these and other limitations of the model and possible future extensions., Comment: Updated to reflect version accepted by ApJ; minor text edits to clarify limitations and caveats of model

In the $\Lambda$-Cold Dark Matter model of the Universe, galaxies form in part through accreting satellite systems. Previous work have built an understanding of the signatures of these processes contained within galactic stellar halos. This work revisits that picture using seven Milky Way-like galaxies in the \textit{Latte} suite of FIRE-2 cosmological simulations. The resolution of these simulations allows a comparison of contributions from satellites above M$_{*}$$\gtrsim$10$\times$$^{7}$M$_{\odot}$, enabling the analysis of observable properties for disrupted satellites in a fully self-consistent and cosmological context. Our results show that, the time of accretion and the stellar mass of an accreted satellite are fundamental parameters that in partnership dictate the resulting spatial distribution, orbital energy, and [$\alpha$/Fe]-[Fe/H] compositions of the stellar debris of such mergers $at$ $present$ $day$. These parameters also govern the resulting dynamical state of an accreted galaxy at $z=0$, leading to the expectation that the inner regions of the stellar halo (R$_{\mathrm{GC}}$ $\lesssim$30 kpc) should contain fully phase-mixed debris from both lower and higher mass satellites. In addition, we find that a significant fraction of the lower mass satellites accreted at early times deposit debris in the outer halo (R$_{\mathrm{GC}}$ $>$50 kpc) that are $not$ fully phased-mixed, indicating that they could be identified in kinematic surveys. Our results suggest that, as future surveys become increasingly able to map the outer halo of our Galaxy, they may reveal the remnants of long-dead dwarf galaxies whose counterparts are too faint to be seen $in$ $situ$ in higher redshift surveys., Comment: Submitted for publication in ApJ

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