Search

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 26 galaxies at $0.5 \lesssim z \lesssim 1.7$ from the MSA-3D 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.05 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: 23 pages, 6 figures, 1 table; submitted to 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: Submitted 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

Die studies are fundamental to quantifying ancient monetary production, providing insights into the relationship between coinage, politics, and history. The process requires tedious manual work, which limits the size of the corpora that can be studied. Few works have attempted to automate this task, and none have been properly released and evaluated from a computer vision perspective. We propose a fully automatic approach that introduces several innovations compared to previous methods. We rely on fast and robust local descriptors matching that is set automatically. Second, the core of our proposal is a clustering-based approach that uses an intrinsic metric (that does not need the ground truth labels) to determine its critical hyper-parameters. We validate the approach on two corpora of Greek coins, propose an automatic implementation and evaluation of previous baselines, and show that our approach significantly outperforms them., Comment: code: https://cea-list-lasti.github.io/projects/studies/studies.html

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: 26 pages, 14 figures, 5 tables, submitted to MNRAS, comments welcome

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

We use simulated galaxy observations from the NIHAO-SKIRT-Catalog to test the accuracy of Spectral Energy Distribution (SED) modeling techniques. SED modeling is an essential tool for inferring star-formation histories from nearby galaxy observations, but is fraught with difficulty due to our incomplete understanding of stellar populations, chemical enrichment processes, and the nonlinear, geometry-dependent effects of dust on our observations. The NIHAO-SKIRT-Catalog uses hydrodynamic simulations and radiative transfer to produce SEDs from the ultraviolet (UV) through the infrared (IR), accounting for the effects of dust. We use the commonly used Prospector software to perform inference on these SEDs, and compare the inferred stellar masses and star-formation rates (SFRs) to the known values in the simulation. We match the stellar population models to isolate the effects of differences in the star-formation history, the chemical evolution history, and the dust. We find that the combined effect of model mismatches for high mass ($> 10^{9.5} M_{\odot}$) galaxies leads to inferred SFRs that are on average underestimated by a factor of 2 when fit to UV through IR photometry, and a factor of 3 when fit to UV through optical photometry. These biases lead to significant inaccuracies in the resulting sSFR-mass relations, with UV through optical fits showing particularly strong deviations from the true relation of the simulated galaxies. In the context of massive existing and upcoming photometric surveys, these results highlight that star-formation history inference from photometry remains imprecise and inaccurate, and that there is a pressing need for more realistic testing of existing techniques.

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

We present the characterization of ultrashort laser pulses by using the plasma-induced frequency resolved optical switching (PI-FROSt) technique, implemented in ambient air. This recently developed method allows for a temporal reconstruction of a pulse at its focal spot by utilizing a moderately intense pump laser pulse for generating an ionization-induced ultrafast defocusing lens. When propagating through the produced plasma lens, the probe beam to characterize experiences an increase of its size in the far field. The spectrum of the defocused probe field, measured as a function of the pump-probe delay, allows for a comprehensive characterization of the temporal and spectral attributes of the pulse. We report herein the ability of this technique, initially designed for use in rare gases, to operate in ambient air conditions with similar performance. The method is remarkably straightforward to implement and requires no additional optical component other than a focusing mirror, while delivering laser pulse reconstructions of high reliability.

Cosmological galaxy formation simulations are powerful tools to understand the complex processes that govern the formation and evolution of galaxies. However, evaluating the realism of these simulations remains a challenge. The two common approaches for evaluating galaxy simulations is either through scaling relations based on a few key physical galaxy properties, or through a set of pre-defined morphological parameters based on galaxy images. This paper proposes a novel image-based method for evaluating the quality of galaxy simulations using unsupervised deep learning anomaly detection techniques. By comparing full galaxy images, our approach can identify and quantify discrepancies between simulated and observed galaxies. As a demonstration, we apply this method to SDSS imaging and NIHAO simulations with different physics models, parameters, and resolution. We further compare the metric of our method to scaling relations as well as morphological parameters. We show that anomaly detection is able to capture similarities and differences between real and simulated objects that scaling relations and morphological parameters are unable to cover, thus indeed providing a new point of view to validate and calibrate cosmological simulations against observed data., Comment: 13 pages, 16 figures, published in MNRAS

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 work is concerned with implementing the hybridizable discontinuous Galerkin (HDG) method to solve the linear anisotropic elastic equation in the frequency domain. First-order formulation with the compliance tensor and Voigt notation are employed to provide a compact description of the discretized problem and flexibility with highly heterogeneous media. We further focus on the question of optimal choices of stabilization in the definition of HDG numerical traces. For this purpose, we construct a hybridized Godunov-upwind flux for anisotropic elastic media possessing three distinct wavespeeds. This stabilization removes the need to choose a scaling factor, contrary to the identity and Kelvin-Christoffel based stabilizations which are popular choices in the literature. We carry out comparisons among these families for isotropic and anisotropic material, with constant background and highly heterogeneous ones, in two and three dimensions. These experiments establish the optimality of the Godunov stabilization which can be used as a reference choice for a generic material in which different types of waves propagate., Comment: 34 pages, 9 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

Solar oscillations can be modeled by Galbrun's equation which describes Lagrangian wave displacement in a self-gravitating stratified medium. For spherically symmetric backgrounds, we construct an algorithm to compute efficiently and accurately the coefficients of the Green's tensor of the time-harmonic equation in vector spherical harmonic basis. With only two resolutions, our algorithm provides values of the kernels for all heights of source and receiver, and prescribes analytically the singularities of the kernels. We also derive absorbing boundary conditions (ABC) to model wave propagation in the atmosphere above the cut-off frequency. The construction of ABC, which contains varying gravity terms, is rendered difficult by the complex behavior of the solar potential in low atmosphere and for frequencies below the Lamb frequency. We carry out extensive numerical investigations to compare and evaluate the efficiency of the ABCs in capturing outgoing solutions. Finally, as an application towards helioseismology, we compute synthetic solar power spectra that contain pressure modes as well as internal-gravity (g-) and surface-gravity (f-) ridges which are missing in simpler approximations of the wave equation. For purpose of validation, the location of the ridges in the synthetic power spectra are compared with observed solar modes., Comment: 34 pages; 12 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

We use simulated attenuation curves from the NIHAO-SKIRT-Catalog to test the flexibility of commonly used dust attenuation models in the face of the variations expected from realistic star-dust geometries. Motivated by lack of flexibility in these existing models, we propose a novel dust attenuation model with three free parameters that can accurately recover the simulated attenuation curves as well as the best-fitting curves from the commonly used models. This new model is fully analytic and treats all starlight equally, in contrast to two-component dust attenuation models. We use the parametrization to investigate the relationship between the overall attenuation law shape and the strength of the 2175 \AA\ bump. Our results indicate variation in star-dust geometry leads these features to correlate tightly, with grayer attenuation curves having weaker bumps.

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

With increasing interest in the use of glassy carbon (GC) for a wide variety of application areas, the need for developing fundamental understanding of its mechanical properties has come to the forefront. Further, recent theoretical and modeling works that shed some light on the synthesis of GC through the process of pyrolysis of polymer precursors have highlighted the possibilities of a revisit to investigation of its mechanical properties at a fundamental level. While there are isolated reports on the experimental determination of its elastic modulus, insights into stress-strain behavior of GC material under tension and compression obtained through simulation, either at molecular level or for the bulk material is missing. This current study fills the gap at the molecular level and investigates the mechanical properties of GC using molecular dynamics (MD) simulations which model the atomistic level formation and breaking of bonds using bond-order based reactive force field formulations. The molecular model considered for this simulation has a characteristics 3D cagey structure of 5-, 6-, and 7-membered carbon rings and graphitic domain of a flat graphene-like structure. The GC molecular model was subjected to loading under varying strain rates (0.4/ns, 0.6/ns, 1.25/ns, and 2.5/ns) and varying temperatures (300 - 800 K) in each of the three axes x, y, and z. The simulation showed that GC molecule has distinct stress-strain curves under tension and compression. In tension, MD modeling predicted mean elastic modulus of 5.71 GPa for a single GC molecule with some dependency on strain rates and temperature, while in compression, the elastic modulus was also found to depend on the strain rates as well as temperature and was predicted to have a mean value of 35 GPa

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

We present simulated galaxy spectral energy distributions (SEDs) from the far ultraviolet through the far infrared, created using hydrodynamic simulations and radiative transfer calculations, suitable for the validation of SED modeling techniques. SED modeling is an essential tool for inferring star formation histories from nearby galaxy observations, but is fraught with difficulty due to our incomplete understanding of stellar populations, chemical enrichment processes, and the non-linear, geometry dependent effects of dust on our observations. Our simulated SEDs will allow us to assess the accuracy of these inferences against galaxies with known ground truth. To create the SEDs, we use simulated galaxies from the NIHAO suite and the radiative transfer code SKIRT. We explore different sub-grid post-processing recipes, using color distributions and their dependence on axis ratio of galaxies in the nearby universe to tune and validate them. We find that sub-grid post-processing recipes that mitigate limitations in the temporal and spatial resolution of the simulations are required for producing FUV to FIR photometry that statistically reproduce the colors of galaxies in the nearby universe. With this paper we release resolved photometry and spatially integrated spectra for our sample galaxies, each from a range of different viewing angles. Our simulations predict that there is a large variation in attenuation laws among galaxies, and that from any particular viewing angle that energy balance between dust attenuation and reemission can be violated by up to a factor of 3. These features are likely to affect SED modeling accuracy.

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

Without AGN feedback, simulated massive, star-forming galaxies become too compact relative to observed galaxies at z<2. In this paper, we perform high-resolution re-simulations of a massive (M_star~10^11 M_sol) galaxy at z~2.3, drawn from the Feedback in Realistic Environments (FIRE) project. In the simulation without AGN feedback, the galaxy experiences a rapid starburst and shrinking of its half-mass radius. We experiment with driving mechanical AGN winds, using a state-of-the-art hyper-Lagrangian refinement technique to increase particle resolution. These winds reduce the gas surface density in the inner regions of the galaxy, suppressing the compact starburst and maintaining an approximately constant half-mass radius. Using radiative transfer, we study the impact of AGN feedback on the magnitude and extent of the multi-wavelength continuum emission. When AGN winds are included, the suppression of the compact, dusty starburst results in lowered flux at FIR wavelengths (due to decreased star formation) but increased flux at optical-to-near-IR wavelengths (due to decreased dust attenuation, in spite of the lowered star formation rate), relative to the case without AGN winds. The FIR half-light radius decreases from ~1 kpc to ~0.1 kpc in <40 Myr when AGN winds are not included, but increases to ~2 kpc when they are. Interestingly, the half-light radius at optical-NIR wavelengths remains approximately constant over 35 Myr, for simulations with and without AGN winds. In the case without winds, this occurs despite the rapid compaction, and is due to heavy dust obscuration in the inner regions of the galaxy. This work highlights the importance of forward-modelling when comparing simulated and observed galaxy populations., Comment: 14 pages, 10 figures. Accepted for publication in MNRAS

Despite their remarkable success in approximating a wide range of operators defined by PDEs, existing neural operators (NOs) do not necessarily perform well for all physics problems. We focus here on high-frequency waves to highlight possible shortcomings. To resolve these, we propose a subfamily of NOs enabling an enhanced empirical approximation of the nonlinear operator mapping wave speed to solution, or boundary values for the Helmholtz equation on a bounded domain. The latter operator is commonly referred to as the ''forward'' operator in the study of inverse problems. Our methodology draws inspiration from transformers and techniques such as stochastic depth. Our experiments reveal certain surprises in the generalization and the relevance of introducing stochastic depth. Our NOs show superior performance as compared with standard NOs, not only for testing within the training distribution but also for out-of-distribution scenarios. To delve into this observation, we offer an in-depth analysis of the Rademacher complexity associated with our modified models and prove an upper bound tied to their stochastic depth that existing NOs do not satisfy. Furthermore, we obtain a novel out-of-distribution risk bound tailored to Gaussian measures on Banach spaces, again relating stochastic depth with the bound. We conclude by proposing a hypernetwork version of the subfamily of NOs as a surrogate model for the mentioned forward operator.

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