Showing 1–50 of 114 results for author: Burkhart, B

Any Way the Wind Blows: Quantifying Superbubbles and their Outflows in Simulated Galaxies across $z \approx 0-3$

Authors: Lori E. Porter, Matthew E. Orr, Blakesley Burkhart, Andrew Wetzel, Dušan Kereš, Claude-André Faucher-Giguère, Philip F. Hopkins

Abstract: 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 ob… ▽ More 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. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Objects May Be Closer Than They Appear: Significant Host Galaxy Dispersion Measures of Fast Radio Bursts in Zoom-in Simulations

Authors: Matthew E. Orr, Blakesley Burkhart, Wenbin Lu, Sam B. Ponnada, Cameron B. Hummels

Abstract: We investigate the contribution of host galaxies to the overall Dispersion Measures (DMs) for Fast Radio Bursts (FRBs) using the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulation suite. We calculate DMs from every star particle in the simulated L* galaxies by ray-tracing through their multi-phase interstellar medium (ISM), summing the line-of-sight free thermal electron c… ▽ More We investigate the contribution of host galaxies to the overall Dispersion Measures (DMs) for Fast Radio Bursts (FRBs) using the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulation suite. We calculate DMs from every star particle in the simulated L* galaxies by ray-tracing through their multi-phase interstellar medium (ISM), summing the line-of-sight free thermal electron column for all gas elements within $\pm$20 kpc of the galactic mid-plane. At $z=0$, we find average (median) host-galaxy DMs of 74 (43) and 210 (94) pc cm$^{-3}$ for older ($\gtrsim$10 Myr) and younger ($\lesssim$10 Myr) stellar populations, respectively. Inclination raises the median DM measured for older populations ($\gtrsim$10 Myr) in the simulations by a factor of $\sim$2, but generally does not affect the younger stars deeply embedded in H{\small II} regions except in extreme edge-on cases (inclination $\gtrsim 85^\circ$). In kinematically disturbed snapshots ($z = 1$ in FIRE), the average (median) host-galaxy DMs are higher: 80 (107) and 266 (795) pc cm$^{-3}$ for older ($\gtrsim$10 Myr) and younger ($\lesssim$10 Myr) stellar populations, respectively. FIRE galaxies tend to have higher DM values than cosmological simulations such as IllustrisTNG. As a result, FRB host galaxies may be closer (lower redshift) than previously inferred. Furthermore, constraining host-galaxy DM distributions may help significantly constrain FRB progenitor models. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: 17 pages, 12 figures, submitted to ApJ Letters

The Supersonic Project: Early Star Formation with the Streaming Velocity

Authors: William Lake, Claire E. Williams, Smadar Naoz, Federico Marinacci, Blakesley Burkhart, Mark Vogelsberger, Naoki Yoshida, Gen Chiaki, Avi Chen, Yeou S. Chiou

Abstract: At high redshifts ($z\gtrsim12$), the relative velocity between baryons and dark matter (the so-called streaming velocity) significantly affects star formation in low-mass objects. Streaming substantially reduces the abundance of low-mass gas objects while simultaneously allowing for the formation of supersonically-induced gas objects (SIGOs) and their associated star clusters outside of dark matt… ▽ More At high redshifts ($z\gtrsim12$), the relative velocity between baryons and dark matter (the so-called streaming velocity) significantly affects star formation in low-mass objects. Streaming substantially reduces the abundance of low-mass gas objects while simultaneously allowing for the formation of supersonically-induced gas objects (SIGOs) and their associated star clusters outside of dark matter halos. Here, we present a study of the population-level effects of streaming on star formation within both halos and SIGOs in a set of simulations with and without streaming. Notably, we find that streaming actually enhances star formation within individual halos of all masses at redshifts between $z=12$ and $z=20$. This is demonstrated both as an increased star formation rate per object as well as an enhancement of the Kennicutt-Schmidt relation for objects with streaming. We find that our simulations are consistent with some observations at high redshift, but on a population level, they continue to under-predict star formation relative to the majority of observations. However, simulations of overdense regions (both with and without streaming) agree with observations, suggesting a strategy for extracting information about the overdensity and streaming velocity in a given survey volume in future observations. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 16 pages, 8 figures

The Interplay between the IMF and Star Formation Efficiency through Radiative Feedback at High Stellar Surface Densities

Authors: Shyam H. Menon, Lachlan Lancaster, Blakesley Burkhart, Rachel S. Somerville, Avishai Dekel, Mark R. Krumholz

Abstract: The observed rest-UV luminosity function at cosmic dawn ($z \sim 8-14$) measured by JWST revealed an excess of UV-luminous galaxies relative to many pre-launch theoretical predictions. A high star-formation efficiency (SFE) and a top-heavy initial mass function (IMF) are among the mechanisms proposed for explaining this excess. Although a top-heavy IMF has been proposed for its ability to increase… ▽ More The observed rest-UV luminosity function at cosmic dawn ($z \sim 8-14$) measured by JWST revealed an excess of UV-luminous galaxies relative to many pre-launch theoretical predictions. A high star-formation efficiency (SFE) and a top-heavy initial mass function (IMF) are among the mechanisms proposed for explaining this excess. Although a top-heavy IMF has been proposed for its ability to increase the light-to-mass ratio (\(Ψ_{\mathrm{UV}}\)), the resulting enhanced radiative pressure from young stars could decrease the star formation efficiency (SFE), potentially driving galaxy luminosities back down. In this Letter, we use idealized radiation hydrodynamic simulations of star cluster formation to explore the effects of a top-heavy IMF on the SFE of clouds typical of the high pressure conditions found at these redshifts. We find that the SFE in star clusters with solar neighbourhood-like dust abundance decreases with increasingly top-heavy IMF's -- by $\sim 20 \%$ for an increase of factor 4 in $Ψ_{\mathrm{UV}}$, and by $50 \%$ for a factor $ \sim 10$ in $Ψ_{\mathrm{UV}}$. However, we find that an expected decrease in the dust-to-gas ratio ($\sim 0.01 \times \mathrm{Solar}$) at these redshifts can completely compensate for the enhanced light output. This leads to a (cloud-scale; $\sim 10 \, \mathrm{pc}$) SFE that is $\gtrsim 70\%$ even for a factor 10 increase in $Ψ_{\mathrm{UV}}$, implying that highly efficient star formation is unavoidable for high surface density and low metallicity conditions. Our results suggest that a top-heavy IMF, if present, likely coexists with efficient star formation in these galaxies. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 12 pages, 4 figures. Accepted to the Astrophysical Journal Letters

The Molecular Cloud Lifecycle II: Formation and Destruction of Molecular Clouds Diagnosed via H$_2$ Fluorescent Emission Emission

Authors: Blakesley Burkhart, Shmuel Bialy, Daniel Seifried, Stefanie Walch, Erika Hamden, Thomas J. Haworth, Keri Hoadley, Shuo Kong, Madisen Johnson, Sarah Jeffreson, Mark R. Krumholz, Min-Young Lee, Amiel Sternberg, Neal J. Turner

Abstract: Molecular hydrogen (H$_2$) formation and dissociation are key processes that drive the gas lifecycle in galaxies. Using the SImulating the LifeCycle of Molecular Clouds (SILCC) zoom-in simulation suite, we explore the utility of future observations of H$_2$ dissociation and formation for tracking the lifecycle of molecular clouds. The simulations used in this work include non-equilibrium H$_2$ for… ▽ More Molecular hydrogen (H$_2$) formation and dissociation are key processes that drive the gas lifecycle in galaxies. Using the SImulating the LifeCycle of Molecular Clouds (SILCC) zoom-in simulation suite, we explore the utility of future observations of H$_2$ dissociation and formation for tracking the lifecycle of molecular clouds. The simulations used in this work include non-equilibrium H$_2$ formation, stellar radiation, sink particles, and turbulence. We find that, at early times in the cloud evolution, H$_2$ formation rapidly outpaces dissociation and molecular clouds build their mass from the atomic reservoir in their environment. Rapid H$_2$ formation is also associated with a higher early star formation rate. For the clouds studied here, H$_2$ is strongly out of chemical equilibrium during the early stages of cloud formation but settles into a bursty chemical steady-state about 2 Myrs after the first stars form. At the latest stage of cloud evolution, dissociation outweighs formation and the clouds enter a dispersal phase. We discuss how theories for the molecular cloud lifecycle and the star formation efficiency may be distinguished with observational measurements of H$_2$ fluorescence with a space-based high-resolution FUV spectrometer, such as the proposed Hyperion and Eos NASA Explorer missions. Such missions would enable measurements of the H$_2$ dissociation and formation rates, which we demonstrate can be connected to different phases in a molecular cloud's star-forming life, including cloud building, rapidly star-forming, H$_2$ chemical equilibrium, and cloud destruction. △ Less

Submitted 5 February, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: Submitted to ApJ, comments welcome

Removing Dust from CMB Observations with Diffusion Models

Authors: David Heurtel-Depeiges, Blakesley Burkhart, Ruben Ohana, Bruno Régaldo-Saint Blancard

Abstract: In cosmology, the quest for primordial $B$-modes in cosmic microwave background (CMB) observations has highlighted the critical need for a refined model of the Galactic dust foreground. We investigate diffusion-based modeling of the dust foreground and its interest for component separation. Under the assumption of a Gaussian CMB with known cosmology (or covariance matrix), we show that diffusion m… ▽ More In cosmology, the quest for primordial $B$-modes in cosmic microwave background (CMB) observations has highlighted the critical need for a refined model of the Galactic dust foreground. We investigate diffusion-based modeling of the dust foreground and its interest for component separation. Under the assumption of a Gaussian CMB with known cosmology (or covariance matrix), we show that diffusion models can be trained on examples of dust emission maps such that their sampling process directly coincides with posterior sampling in the context of component separation. We illustrate this on simulated mixtures of dust emission and CMB. We show that common summary statistics (power spectrum, Minkowski functionals) of the components are well recovered by this process. We also introduce a model conditioned by the CMB cosmology that outperforms models trained using a single cosmology on component separation. Such a model will be used in future work for diffusion-based cosmological inference. △ Less

Submitted 11 December, 2023; v1 submitted 24 October, 2023; originally announced October 2023.

Comments: 5+6 pages, 2+3 figures, accepted at the NeurIPS 2023 workshop on "Machine Learning and the Physical Sciences" and selected for a spotlight talk

The Supersonic Project: Lighting up the faint end of the JWST UV luminosity function

Authors: Claire E. Williams, William Lake, Smadar Naoz, Blakesley Burkhart, Tommaso Treu, Federico Marinacci, Yurina Nakazato, Mark Vogelsberger, Naoki Yoshida, Gen Chiaki, Yeou S. Chiou, Avi Chen

Abstract: The James Webb Space Telescope (JWST) is capable of probing extremely early eras of our Universe when the supersonic relative motions between dark matter and baryonic overdensities modulate structure formation ($z>\sim 10$). We study low-mass galaxy formation including this "stream velocity" using high resolution AREPO hydrodynamics simulations, and present theoretical predictions of the UV lumino… ▽ More The James Webb Space Telescope (JWST) is capable of probing extremely early eras of our Universe when the supersonic relative motions between dark matter and baryonic overdensities modulate structure formation ($z>\sim 10$). We study low-mass galaxy formation including this "stream velocity" using high resolution AREPO hydrodynamics simulations, and present theoretical predictions of the UV luminosity function (UVLF) and galaxy stellar mass function (GSMF) down to extremely faint and low mass galaxies ($M_{UV}>\sim-15$, $10^4M_\odot<=M_*<=10^8 M_\odot)$. We show that, although the stream velocity suppresses early star formation overall, it induces a short period of rapid star formation in some larger dwarfs, leading to an enhancement in the faint-end of the UVLF at $z=12$. We demonstrate that JWST observations are close to this enhanced regime, and propose that the UVLF may constitute an important probe of the stream velocity at high redshift for JWST and future observatories. △ Less

Submitted 15 December, 2023; v1 submitted 5 October, 2023; originally announced October 2023.

Comments: 13 pages, 7 figures. Accepted ApJL

Probing the Conditions for the Hı-to-H$_{2}$ Transition in the Interstellar Medium

Authors: Gyueun Park, Min-Young Lee, Shmuel Bialy, Blakesley Burkhart, J. R. Dawson, Carl Heiles, Di Li, Claire Murray, Hiep Nguyen, Anita Hafner, Daniel R. Rybarczyk, Snežana Stanimirović

Abstract: In this paper, we investigate the conditions for the HI-to-H$_{2}$ transition in the solar neighborhood by analyzing HI emission and absorption measurements toward 58 Galactic lines of sight (LOSs) along with $^{12}$CO(1$-$0) (CO) and dust data. Based on the accurate column densities of the cold and warm neutral medium (CNM and WNM), we first perform a decomposition of gas into atomic and molecula… ▽ More In this paper, we investigate the conditions for the HI-to-H$_{2}$ transition in the solar neighborhood by analyzing HI emission and absorption measurements toward 58 Galactic lines of sight (LOSs) along with $^{12}$CO(1$-$0) (CO) and dust data. Based on the accurate column densities of the cold and warm neutral medium (CNM and WNM), we first perform a decomposition of gas into atomic and molecular phases and show that the observed LOSs are mostly HI-dominated. In addition, we find that the CO-dark H$_{2}$, not the optically thick HI, is a major ingredient of the dark gas in the solar neighborhood. To examine the conditions for the formation of CO-bright molecular gas, we analyze the kinematic association between HI and CO and find that the CNM is kinematically more closely associated with CO than the WNM. When CNM components within CO line widths are isolated, we find the following characteristics: spin temperature $<$ 200 K, peak optical depth $>$ 0.1, CNM fraction of $\sim$0.6, and $V$-band dust extinction $>$ 0.5 mag. These results suggest that CO-bright molecular gas preferentially forms in environments with high column densities where the CNM becomes colder and more abundant. Finally, we confront the observed CNM properties with the steady-state H$_{2}$ formation model of Sternberg et al. and infer that the CNM must be clumpy with a small volume filling factor. Another possibility would be that missing processes in the model, such as cosmic-rays and gas dynamics, play an important role in the HI-to-H$_{2}$ transition. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: 25 pages, 12 figures, 6 tables, Accepted for publication in ApJ

An Exploration of AGN and Stellar Feedback Effects in the Intergalactic Medium via the Low Redshift Lyman-$α$ Forest

Authors: Megan Taylor Tillman, Blakesley Burkhart, Stephanie Tonnesen, Simeon Bird, Greg L. Bryan, Daniel Anglés-Alcázar, Sultan Hassan, Rachel S. Somerville, Romeel Davé, Federico Marinacci, Lars Hernquist, Mark Vogelsberger

Abstract: We explore the role of galactic feedback on the low redshift Lyman-$α$ (Ly$α$)~forest ($z \lesssim 2$) statistics and its potential to alter the thermal state of the intergalactic medium. Using the Cosmology and Astrophysics with Machine Learning Simulations (CAMELS) suite, we explore variations of the AGN and stellar feedback models in the IllustrisTNG and Simba sub-grid models. We find that both… ▽ More We explore the role of galactic feedback on the low redshift Lyman-$α$ (Ly$α$)~forest ($z \lesssim 2$) statistics and its potential to alter the thermal state of the intergalactic medium. Using the Cosmology and Astrophysics with Machine Learning Simulations (CAMELS) suite, we explore variations of the AGN and stellar feedback models in the IllustrisTNG and Simba sub-grid models. We find that both AGN and stellar feedback in Simba play a role in setting the Ly$α$ forest column density distribution function (CDD) and the Doppler width ($b$-value) distribution. The Simba AGN jet feedback mode is able to efficiently transport energy out to the diffuse IGM causing changes in the shape and normalization of the CDD and a broadening of the $b$-value distribution. We find that stellar feedback plays a prominent role in regulating supermassive black hole growth and feedback, highlighting the importance of constraining stellar and AGN feedback simultaneously. In IllustrisTNG, the AGN feedback variations explored in CAMELS do not affect the Ly$α$ forest, but varying the stellar feedback model does produce subtle changes. Our results imply that the low-$z$ Ly$α$ forest can be sensitive to changes in the ultraviolet background (UVB), stellar and black hole feedback, and that AGN jet feedback in particular can have a strong effect on the thermal state of the IGM. △ Less

Submitted 1 November, 2023; v1 submitted 12 July, 2023; originally announced July 2023.

Comments: 27 pages, 12 figures, 2 tables, accepted to ApJ

Extending a Physics-Informed Machine Learning Network for Superresolution Studies of Rayleigh-Bénard Convection

Authors: Diane M. Salim, Blakesley Burkhart, David Sondak

Abstract: Advancing our understanding of astrophysical turbulence is bottlenecked by the limited resolution of numerical simulations that may not fully sample scales in the inertial range. Machine learning (ML) techniques have demonstrated promise in up-scaling resolution in both image analysis and numerical simulations (i.e., superresolution). Here we employ and further develop a physics-constrained convol… ▽ More Advancing our understanding of astrophysical turbulence is bottlenecked by the limited resolution of numerical simulations that may not fully sample scales in the inertial range. Machine learning (ML) techniques have demonstrated promise in up-scaling resolution in both image analysis and numerical simulations (i.e., superresolution). Here we employ and further develop a physics-constrained convolutional neural network (CNN) ML model called "MeshFreeFlowNet'' (MFFN) for superresolution studies of turbulent systems. The model is trained both on the simulation images as well as the evaluated PDEs, making it sensitive to the underlying physics of a particular fluid system. We develop a framework for 2D turbulent Rayleigh-Bénard convection (RBC) generated with the \textsc{Dedalus} code by modifying the MFFN architecture to include the full set of simulation PDEs and the boundary conditions. Our training set includes fully developed turbulence sampling Rayleigh numbers ($Ra$) of $Ra=10^6-10^{10}$. We evaluate the success of the learned simulations by comparing the power spectra of the direct \textsc{Dedalus} simulation to the predicted model output, and compare both ground truth and predicted power spectral inertial range scalings to theoretical predictions. We find that the updated network performs well at all $Ra$ studied here in recovering large-scale information, including the inertial range slopes. The superresolution prediction is overly dissipative at smaller scales than that of the inertial range in all cases, but the smaller-scales are better recovered in more turbulent, than laminar, regimes. This is likely because more turbulent systems have a rich variety of structures at many length scales compared to laminar flows. △ Less

Submitted 31 January, 2024; v1 submitted 5 July, 2023; originally announced July 2023.

Comments: 11 pages, 10 figures

The Supersonic Project: Star Formation in Early Star Clusters without Dark Matter

Authors: William Lake, Smadar Naoz, Federico Marinacci, Blakesley Burkhart, Mark Vogelsberger, Claire E. Williams, Yeou S. Chiou, Gen Chiaki, Yurina Nakazato, Naoki Yoshida

Abstract: The formation mechanism of globular clusters (GCs) has long been debated by astronomers. It was recently proposed that Supersonically Induced Gas Objects (SIGOs), which formed in the early Universe due to the supersonic relative motion of baryons and dark matter at recombination, could be the progenitors of early globular clusters. In order to become GCs, SIGOs must form stars relatively efficient… ▽ More The formation mechanism of globular clusters (GCs) has long been debated by astronomers. It was recently proposed that Supersonically Induced Gas Objects (SIGOs), which formed in the early Universe due to the supersonic relative motion of baryons and dark matter at recombination, could be the progenitors of early globular clusters. In order to become GCs, SIGOs must form stars relatively efficiently despite forming outside of dark matter halos. We investigate the potential for star formation in SIGOs using cosmological hydrodynamic simulations, including the aforementioned relative motions of baryons and dark matter, molecular hydrogen cooling in primordial gas clouds, and including explicit star formation. We find that SIGOs do form stars and that the nascent star clusters formed through this process are accreted by dark matter halos on short timescales (a few hundreds of Myr). Thus, SIGOs may be found as intact substructures within these halos, analogous to many present-day GCs. From this result, we conclude that SIGOs are capable of forming star clusters with similar properties to globular clusters in the early Universe and we discuss their detectablity by upcoming JWST surveys. △ Less

Submitted 18 September, 2023; v1 submitted 1 June, 2023; originally announced June 2023.

Comments: 12 pages, 5 figures

Fitting Probability Distribution Functions in Turbulent Star-Forming Molecular Clouds

Authors: Avery Kiihne, Sabrina M. Appel, Blakesley Burkhart, Vadim A. Semenov, Christoph Federrath

Abstract: We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback and magnetic fields, to investigate the effectiveness of different fitting methods for volume and column density probability distribution functions (PDFs). The first method fits a piecewise lognormal and power-law (PL) function to recover PDF parameters such as the PL slope and transition density.… ▽ More We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback and magnetic fields, to investigate the effectiveness of different fitting methods for volume and column density probability distribution functions (PDFs). The first method fits a piecewise lognormal and power-law (PL) function to recover PDF parameters such as the PL slope and transition density. The second method fits a polynomial spline function and examines the first and second derivatives of the spline to determine the PL slope and the functional transition density. We demonstrate that fitting a spline allows us to directly determine if the data has multiple PL slopes. The first PL (set by the transition between lognormal and PL function) can also be visualized in the derivatives directly. In general, the two methods produce fits that agree reasonably well for volume density but vary for column density, likely due to the increased statistical noise in column density maps as compared to volume density. We test a well-known conversion for estimating volume density PL slopes from column density slopes and find that the spline method produces a better match (\c{hi}2 of 2.38 vs \c{hi}2 of 5.92), albeit with a significant scatter. Ultimately, we recommend the use of both fitting methods on column density data to mitigate the effects of noise. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: 14 pages, 8 figures, Submitted to ApJ. Comments welcome

JWST constraints on the UV luminosity density at cosmic dawn: implications for 21-cm cosmology

Authors: Sultan Hassan, Christopher C. Lovell, Piero Madau, Marc Huertas-Company, Rachel S. Somerville, Blakesley Burkhart, Keri L. Dixon, Robert Feldmann, Tjitske K. Starkenburg, John F. Wu, Christian Kragh Jespersen, Joseph D. Gelfand, Ankita Bera

Abstract: An unprecedented array of new observational capabilities are starting to yield key constraints on models of the epoch of first light in the Universe. In this Letter we discuss the implications of the UV radiation background at cosmic dawn inferred by recent JWST observations for radio experiments aimed at detecting the redshifted 21-cm hyperfine transition of diffuse neutral hydrogen. Under the ba… ▽ More An unprecedented array of new observational capabilities are starting to yield key constraints on models of the epoch of first light in the Universe. In this Letter we discuss the implications of the UV radiation background at cosmic dawn inferred by recent JWST observations for radio experiments aimed at detecting the redshifted 21-cm hyperfine transition of diffuse neutral hydrogen. Under the basic assumption that the 21-cm signal is activated by the Ly$α$ photon field produced by metal-poor stellar systems, we show that a detection at the low frequencies of the EDGES and SARAS3 experiments may be expected from a simple extrapolation of the declining UV luminosity density inferred at $z\lesssim 14$ from JWST early galaxy data. Accounting for an early radiation excess above the CMB suggests a shallower or flat evolution to simultaneously reproduce low and high-$z$ current UV luminosity density constraints, which cannot be entirely ruled out, given the large uncertainties from cosmic variance and the faint-end slope of the galaxy luminosity function at cosmic dawn. Our findings raise the intriguing possibility that a high star formation efficiency at early times may trigger the onset of intense Ly$α$ emission at redshift $z\lesssim 20$ and produce a cosmic 21-cm absorption signal 200 Myr after the Big Bang. △ Less

Submitted 11 October, 2023; v1 submitted 4 May, 2023; originally announced May 2023.

Comments: 8 pages, 3 figures, accepted to ApJL

Molecular gas and star formation in nearby starburst galaxy mergers

Authors: Hao He, Connor Bottrell, Christine Wilson, Jorge Moreno, Blakesley Burkhart, Christopher C. Hayward, Lars Hernquist, Angela Twum

Abstract: We employ the Feedback In Realistic Environments (FIRE-2) physics model to study how the properties of giant molecular clouds (GMCs) evolve during galaxy mergers. We conduct a pixel-by-pixel analysis of molecular gas properties in both the simulated control galaxies and galaxy major mergers. The simulated GMC-pixels in the control galaxies follow a similar trend in a diagram of velocity dispersion… ▽ More We employ the Feedback In Realistic Environments (FIRE-2) physics model to study how the properties of giant molecular clouds (GMCs) evolve during galaxy mergers. We conduct a pixel-by-pixel analysis of molecular gas properties in both the simulated control galaxies and galaxy major mergers. The simulated GMC-pixels in the control galaxies follow a similar trend in a diagram of velocity dispersion ($σ_v$) versus gas surface density ($Σ_{\mathrm{mol}}$) to the one observed in local spiral galaxies in the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) survey. For GMC-pixels in simulated mergers, we see a significant increase of factor of 5 - 10 in both $Σ_{\mathrm{mol}}$ and $σ_v$, which puts these pixels above the trend of PHANGS galaxies in the $σ_v$ vs $Σ_{\mathrm{mol}}$ diagram. This deviation may indicate that GMCs in the simulated mergers are much less gravitationally bound compared with simulated control galaxies with virial parameter ($α_{\mathrm{vir}}$) reaching 10 - 100. Furthermore, we find that the increase in $α_{\mathrm{vir}}$ happens at the same time as the increase in global star formation rate (SFR), which suggests stellar feedback is responsible for dispersing the gas. We also find that the gas depletion time is significantly lower for high $α_{\mathrm{vir}}$ GMCs during a starburst event. This is in contrast to the simple physical picture that low $α_{\mathrm{vir}}$ GMCs are easier to collapse and form stars on shorter depletion times. This might suggest that some other physical mechanisms besides self-gravity are helping the GMCs in starbursting mergers collapse and form stars. △ Less

Submitted 31 March, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: 22 pages, 11 figures. Accepted to ApJ. Link to animation updated

What Sets the Star Formation Rate of Molecular Clouds? The Density Distribution as a Fingerprint of Compression and Expansion Rates

Authors: Sabrina M. Appel, Blakesley Burkhart, Vadim A. Semenov, Christoph Federrath, Anna L. Rosen, Jonathan C. Tan

Abstract: We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycl… ▽ More We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycling between high and low density states. We find that the inclusion of protostellar jets produces rapidly expanding and compressing low-density gas. We find that the gas mass flux peaks at the transition between the lognormal and power-law forms of the density probability distribution function (PDF). This is consistent with the transition density tracking the post-shock density, which promotes an enhancement of mass at this density (i.e., shock compression and filament formation). At high densities, the gas dynamics are dominated by self-gravity: the compression rate in all of our runs matches the rate of the run with only gravity, suggesting that processes other than self-gravity have little effect at these densities. The net gas mass flux becomes constant at a density below the sink formation threshold, where it equals the star formation rate. The density at which the net gas mass flux equals the star formation rate is one order of magnitude lower than our sink threshold density, corresponds to the formation of the second power-law tail in the density PDF, and sets the overall star formation rates of these simulations. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 22 pages, 13 figures. Submitted to ApJ. Comments welcome

Arkenstone I: A Novel Method for Robustly Capturing High Specific Energy Outflows In Cosmological Simulations

Authors: Matthew C. Smith, Drummond B. Fielding, Greg L. Bryan, Chang-Goo Kim, Eve C. Ostriker, Rachel S. Somerville, Jonathan Stern, Kung-Yi Su, Rainer Weinberger, Chia-Yu Hu, John C. Forbes, Lars Hernquist, Blakesley Burkhart, Yuan Li

Abstract: Arkenstone is a new model for multiphase, stellar feedback driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealised non-cosmological simulations of a galaxy with a realistic CGM, using the Arep… ▽ More Arkenstone is a new model for multiphase, stellar feedback driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealised non-cosmological simulations of a galaxy with a realistic CGM, using the Arepo code. Hot, fast gas phases with low mass loadings are predicted to dominate the energy content of multiphase outflows. In order to treat the huge dynamic range of spatial scales involved in cosmological galaxy formation at feasible computational expense, cosmological volume simulations typically employ a Lagrangian code or else use adaptive mesh refinement with a quasi-Lagrangian refinement strategy. However, it is difficult to inject a high specific energy wind in a Lagrangian scheme without incurring artificial burstiness. Additionally, the low densities inherent to this type of flow result in poor spatial resolution. Arkenstone addresses these issues with a novel scheme for coupling energy into the ISM/CGM transition region which also provides the necessary level of refinement at the base of the wind. In the absence of our improvements, we show that poor spatial resolution near the sonic point of a hot, fast outflow leads to an underestimation of gas acceleration as the wind propagates. We explore the different mechanisms by which low and high specific energy winds can regulate the SFR of galaxies. In future work, we will demonstrate other aspects of the Arkenstone model. △ Less

Submitted 7 March, 2024; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: Published MNRAS, 27 pages, 17 figures; updated for consistency with journal version

Journal ref: MNRAS, 527, 1216 (2024)

Hyperion: The origin of the stars A far-UV space telescope for high-resolution spectroscopy over wide fields

Authors: Erika Hamden, David Schiminovich, Shouleh Nikzad, Neal J. Turner, Blakesley Burkhart, Thomas J. Haworth, Keri Hoadley, Jinyoung Serena Kim, Shmuel Bialyh, Geoff Bryden, Haeun Chung, Nia Imara, Rob Kennicutt, Jorge Pineda, Shuo Konga, Yasuhiro Hasegawa, Ilaria Pascucci, Benjamin Godard, Mark Krumholz, Min-Young Lee, Daniel Seifried, Amiel Sternberg, Stefanie Walch, Miles Smith, Stephen C. Unwin , et al. (8 additional authors not shown)

Abstract: We present Hyperion, a mission concept recently proposed to the December 2021 NASA Medium Explorer announcement of opportunity. Hyperion explores the formation and destruction of molecular clouds and planet-forming disks in nearby star-forming regions of the Milky Way. It does this using long-slit, high-resolution spectroscopy of emission from fluorescing molecular hydrogen, which is a powerful fa… ▽ More We present Hyperion, a mission concept recently proposed to the December 2021 NASA Medium Explorer announcement of opportunity. Hyperion explores the formation and destruction of molecular clouds and planet-forming disks in nearby star-forming regions of the Milky Way. It does this using long-slit, high-resolution spectroscopy of emission from fluorescing molecular hydrogen, which is a powerful far-ultraviolet (FUV) diagnostic. Molecular hydrogen (H2) is the most abundant molecule in the universe and a key ingredient for star and planet formation, but is typically not observed directly because its symmetric atomic structure and lack of a dipole moment mean there are no spectral lines at visible wavelengths and few in the infrared. Hyperion uses molecular hydrogen's wealth of FUV emission lines to achieve three science objectives: (1) determining how star formation is related to molecular hydrogen formation and destruction at the boundaries of molecular clouds; (2) determining how quickly and by what process massive star feedback disperses molecular clouds; and (3) determining the mechanism driving the evolution of planet-forming disks around young solar-analog stars. Hyperion conducts this science using a straightforward, highly-efficient, single-channel instrument design. Hyperion's instrument consists of a 48 cm primary mirror, with an f/5 focal ratio. The spectrometer has two modes, both covering 138.5-161.5 nm bandpasses. A low resolution mode has a spectral resolution of R>10,000 with a slit length of 65 arcmin, while the high resolution mode has a spectral resolution of R>50,000 over a slit length of 5 armin. Hyperion occupies a 2 week long, high-earth, Lunar resonance TESS-like orbit, and conducts 2 weeks of planned observations per orbit, with time for downlinks and calibrations. Hyperion was reviewed as Category I, which is the highest rating possible, but was not selected. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: Accepted to JATIS, 9 Figures

X-ray Absorption Lines in the Warm-Hot Intergalactic Medium: Probing Chandra observations with the CAMEL simulations

Authors: Amanda Butler Contreras, Erwin T. Lau, Benjamin D. Oppenheimer, Ákos Bogdán, Megan Tillman, Daisuke Nagai, Orsolya E. Kovács, Blakesley Burkhart

Abstract: Known as the "Missing Baryon Problem", about one-third of baryons in the local universe remain unaccounted for. The missing baryons are thought to reside in the warm-hot intergalactic medium (WHIM) of the cosmic web filaments, which are challenging to detect. Recent Chandra X-ray observations used a novel stacking analysis and detected an OVII absorption line toward the sightline of a luminous qua… ▽ More Known as the "Missing Baryon Problem", about one-third of baryons in the local universe remain unaccounted for. The missing baryons are thought to reside in the warm-hot intergalactic medium (WHIM) of the cosmic web filaments, which are challenging to detect. Recent Chandra X-ray observations used a novel stacking analysis and detected an OVII absorption line toward the sightline of a luminous quasar, hinting that the missing baryons may reside in the WHIM. To explore how the properties of the OVII absorption line depend on feedback physics, we compare the observational results with predictions obtained from the Cosmology and Astrophysics with MachinE Learning (CAMEL) Simulation suite. CAMELS consists of cosmological simulations with state-of-the-art supernova (SN) and active galactic nuclei (AGN) feedback models from the IllustrisTNG and SIMBA simulations, with varying strengths. We find that the simulated OVII column densities are higher in the outskirts of galaxies than in the large-scale WHIM, but they are consistently lower than those obtained in the Chandra observations, for all feedback runs. We establish that the OVII distribution is primarily sensitive to changes in the SN feedback prescription, whereas changes in the AGN feedback prescription have minimal impact. We also find significant differences in the OVII column densities between the IllustrisTNG and SIMBA runs. We conclude that the tension between the observed and simulated OVII column densities cannot be explained by the wide range of feedback models implemented in CAMELS. △ Less

Submitted 12 January, 2023; v1 submitted 28 November, 2022; originally announced November 2022.

Comments: 11 pages, 9 figures. Matching to the published version on MNRAS

Journal ref: MNRAS, Volume 519, Issue 2, February 2023, Pages 2251-2261

Active galactic nucleus jet feedback in hydrostatic halos

Authors: Rainer Weinberger, Kung-Yi Su, Kristian Ehlert, Christoph Pfrommer, Lars Hernquist, Greg L. Bryan, Volker Springel, Yuan Li, Blakesley Burkhart, Ena Choi, Claude-André Faucher-Giguère

Abstract: 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 funct… ▽ More 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. △ Less

Submitted 5 May, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: 22 pages, 20 figures, accepted for publication in MNRAS, comments welcome

CARMA-NRO Orion Survey: unbiased survey of dense cores and core mass functions in Orion A

Authors: Hideaki Takemura, Fumitaka Nakamura, Héctor G. Arce, Nicola Schneider, Volker Ossenkopf-Okada, Shuo Kong, Shun Ishii, Kazuhito Dobashi, Tomomi Shimoikura, Patricio Sanhueza, Takashi Tsukagoshi, Paolo Padoan, Ralf S. Klessen, Paul. F. Goldsmith, Blakesley Burkhart, Dariusz C. Lis Álvaro Sánchez-Monge, Yoshito Shimajiri, Ryohei Kawabe

Abstract: The mass distribution of dense cores is a potential key to understand the process of star formation. Applying dendrogram analysis to the CARMA-NRO Orion C$^{18}$O ($J$=1--0) data, we identify 2342 dense cores, about 22 \% of which have virial ratios smaller than 2, and can be classified as gravitationally bound cores. The derived core mass function (CMF) for bound starless cores which are not asso… ▽ More The mass distribution of dense cores is a potential key to understand the process of star formation. Applying dendrogram analysis to the CARMA-NRO Orion C$^{18}$O ($J$=1--0) data, we identify 2342 dense cores, about 22 \% of which have virial ratios smaller than 2, and can be classified as gravitationally bound cores. The derived core mass function (CMF) for bound starless cores which are not associate with protostars has a slope similar to Salpeter's initial mass function (IMF) for the mass range above 1 $M_\odot$, with a peak at $\sim$ 0.1 $M_\odot$. We divide the cloud into four parts based on the declination, OMC-1/2/3, OMC-4/5, L1641N/V380 Ori, and L1641C, and derive the CMFs in these regions. We find that starless cores with masses greater than 10 $M_\odot$ exist only in OMC-1/2/3, whereas the CMFs in OMC-4/5, L1641N, and L1641C are truncated at around 5--10 $M_\odot$. From the number ratio of bound starless cores and Class II objects in each subregion, the lifetime of bound starless cores is estimated to be 5--30 free-fall times, consistent with previous studies for other regions. In addition, we discuss core growth by mass accretion from the surrounding cloud material to explain the coincidence of peak masses between IMFs and CMFs. The mass accretion rate required for doubling the core mass within a core lifetime is larger than that of Bondi-Hoyle accretion by a factor of order 2. This implies that more dynamical accretion processes are required to grow cores. △ Less

Submitted 21 November, 2022; v1 submitted 18 November, 2022; originally announced November 2022.

Comments: 58 pages, 33 figures, 21 tables, accepted by ApJS

The Supersonic Project: The eccentricity and rotational support of SIGOs and DM GHOSts

Authors: Claire E. Williams, Smadar Naoz, William Lake, Yeou S. Chiou, Blakesley Burkhart, Federico Marinacci, Mark Vogelsberger, Gen Chiaki, Yurina Nakazato, Naoki Yoshida

Abstract: A supersonic relative velocity between dark matter (DM) and baryons (the stream velocity) at the time of recombination induces the formation of low mass objects with anomalous properties in the early Universe. We widen the scope of the `Supersonic Project' paper series to include objects we term Dark Matter + Gas Halos Offset by Streaming (DM GHOSts)--diffuse, DM-enriched structures formed because… ▽ More A supersonic relative velocity between dark matter (DM) and baryons (the stream velocity) at the time of recombination induces the formation of low mass objects with anomalous properties in the early Universe. We widen the scope of the `Supersonic Project' paper series to include objects we term Dark Matter + Gas Halos Offset by Streaming (DM GHOSts)--diffuse, DM-enriched structures formed because of a physical offset between the centers of mass of DM and baryonic overdensities. We present an updated numerical investigation of DM GHOSts and Supersonically Induced Gas Objects (SIGOs), including the effects of molecular cooling, in high resolution hydrodynamic simulations using the AREPO code. Supplemented by an analytical understanding of their ellipsoidal gravitational potentials, we study the population-level properties of these objects, characterizing their morphology, spin, radial mass, and velocity distributions in comparison to classical structures in non-streaming regions. The stream velocity causes deviations from sphericity in both the gas and DM components and lends greater rotational support to the gas. Low mass ($<\sim 10^{5.5}$ M$_\odot$) objects in regions of streaming demonstrate core-like rotation and mass profiles. Anomalies in the rotation and morphology of DM GHOSts could represent an early Universe analogue to observed ultra-faint dwarf galaxies with variations in DM content and unusual rotation curves. △ Less

Submitted 13 February, 2023; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 27 pages, 20 figures. Accepted ApJ

Efficient long-range AGN feedback affects the low redshift Lyman-$α$ forest

Authors: Megan Taylor Tillman, Blakesley Burkhart, Stephanie Tonnesen, Simeon Bird, Greg L. Bryan, Daniel Anglés-Alcázar, Romeel Davé, Shy Genel

Abstract: Active galactic nuclei (AGN) feedback models are generally calibrated to reproduce galaxy observables such as the stellar mass function and the bimodality in galaxy colors. We use variations of the AGN feedback implementations in the IllustrisTNG (TNG) and Simba cosmological hydrodynamic simulations to show that the low redshift Lyman-$α$ forest can provide constraints on the impact of AGN feedbac… ▽ More Active galactic nuclei (AGN) feedback models are generally calibrated to reproduce galaxy observables such as the stellar mass function and the bimodality in galaxy colors. We use variations of the AGN feedback implementations in the IllustrisTNG (TNG) and Simba cosmological hydrodynamic simulations to show that the low redshift Lyman-$α$ forest can provide constraints on the impact of AGN feedback. We show that TNG over-predicts the number density of absorbers at column densities $N_{\rm HI} < 10^{14}$ cm$^{-2}$ compared to data from the Cosmic Origins Spectrograph (in agreement with previous work), and we demonstrate explicitly that its kinetic feedback mode, which is primarily responsible for galaxy quenching, has a negligible impact on the column density distribution (CDD) of absorbers. In contrast, we show that the fiducial Simba model which includes AGN jet feedback is the preferred fit to the observed CDD of the $z = 0.1$ Lyman-$α$ forest across five orders of magnitude in column density. We show that the Simba results with jets produce a quantitatively better fit to the observational data than the Simba results without jets, even when the UVB is left as a free parameter. AGN jets in Simba are high speed, collimated, weakly-interacting with the interstellar medium (via brief hydrodynamic decoupling) and heated to the halo virial temperature. Collectively these properties result in stronger long-range impacts on the IGM when compared to TNG's kinetic feedback mode, which drives isotropic winds with lower velocities at the galactic radius. Our results suggest that the low redshift Lyman-$α$ forest provides plausible evidence for long-range AGN jet feedback. △ Less

Submitted 7 March, 2023; v1 submitted 5 October, 2022; originally announced October 2022.

Comments: 15 pages, 5 figures, 1 table, accepted to ApJL

Journal ref: ApJL 945 L17 (2023)

Spiral Arms are Metal Freeways: Azimuthal Gas-Phase Metallicity Variations in Simulated Cosmological Zoom-in Flocculent Disks

Authors: Matthew E. Orr, Blakesley Burkhart, Andrew Wetzel, Philip F. Hopkins, Ivanna A. Escala, Allison L. Strom, Paul F. Goldsmith, Jorge L. Pineda, Christopher C. Hayward, Sarah R. Loebman

Abstract: We examine the azimuthal variations in gas-phase metallicity profiles in simulated Milky Way mass disk galaxies from the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulation suite, which includes a sub-grid turbulent metal mixing model. We produce spatially resolved maps of the disks at $z \approx 0$ with pixel sizes ranging from 250 to 750~pc, analogous to modern integral f… ▽ More We examine the azimuthal variations in gas-phase metallicity profiles in simulated Milky Way mass disk galaxies from the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulation suite, which includes a sub-grid turbulent metal mixing model. We produce spatially resolved maps of the disks at $z \approx 0$ with pixel sizes ranging from 250 to 750~pc, analogous to modern integral field unit (IFU) galaxy surveys, mapping the gas-phase metallicities in both the cold & dense gas and the ionized gas correlated with HII regions. We report that the spiral arms alternate in a pattern of metal rich and metal poor relative to the median metallicity on the order of $\lesssim 0.1$~dex, appearing generally in this sample of flocculent spirals. The pattern persists even in a simulation with different strengths of metal mixing, indicating that the pattern emerges from physics above the sub-grid scale. Local enrichment does not appear to be the dominant source of the azimuthal metallicity variations at $z \approx 0$: there is no correlation with local star formation on these spatial scales. Rather, the arms are moving inwards and outwards relative to each other, carrying their local metallicity gradients with them radially before mixing into the larger-scale interstellar medium. We propose that the arms act as freeways channeling relatively metal poor gas radially inwards, and relatively enriched gas radially outwards. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 16 Pages, 14 Figures, submitted to MNRAS

Code Comparison in Galaxy Scale Simulations with Resolved Supernova Feedback: Lagrangian vs. Eulerian Methods

Authors: Chia-Yu Hu, Matthew C. Smith, Romain Teyssier, Greg L. Bryan, Robbert Verbeke, Andrew Emerick, Rachel S. Somerville, Blakesley Burkhart, Yuan Li, John C. Forbes, Tjitske Starkenburg

Abstract: We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: {\sc Gizmo}, {\sc Arepo}, {\sc Gadget}, and {\sc Ramses}. All codes adopt the same physical model which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to sub-grid mo… ▽ More We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: {\sc Gizmo}, {\sc Arepo}, {\sc Gadget}, and {\sc Ramses}. All codes adopt the same physical model which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to sub-grid models, eliminating one of the major uncertainties in cosmological simulations. We find reasonable agreement on the time-averaged star formation rates as well as the joint density-temperature distributions between all codes. However, the Lagrangian codes show significantly burstier star formation, larger supernova-driven bubbles, and stronger galactic outflows compared to the Eulerian code. This is caused by the behavior in the dense, collapsing gas clouds when the Jeans length becomes unresolved: gas in Lagrangian codes collapses to much higher densities than in Eulerian codes, as the latter is stabilized by the minimal cell size. Therefore, more of the gas cloud is converted to stars and SNe are much more clustered in the Lagrangian models, amplifying their dynamical impact. The differences between Lagrangian and Eulerian codes can be reduced by adopting a higher star formation efficiency in Eulerian codes, which significantly enhances SN clustering in the latter. Adopting a zero SN delay time reduces burstiness in all codes, resulting in vanishing outflows as SN clustering is suppressed. △ Less

Submitted 3 May, 2023; v1 submitted 22 August, 2022; originally announced August 2022.

Comments: accepted version by ApJ (including a new simulation in Appendix B suggested by the referee)

The Supersonic Project: The Early Evolutionary Path of SIGOs

Authors: William Lake, Smadar Naoz, Blakesley Burkhart, Federico Marinacci, Mark Vogelsberger, Gen Chiaki, Yeou S. Chiou, Naoki Yoshida, Yurina Nakazato, Claire E. Williams

Abstract: Supersonically Induced Gas Objects (SIGOs) are a class of early Universe objects that have gained attention as a potential formation route for globular clusters. SIGOs have only recently begun to be studied in the context of molecular hydrogen cooling, which is key to characterizing their structure and evolution. Studying the population-level properties of SIGOs with molecular cooling is important… ▽ More Supersonically Induced Gas Objects (SIGOs) are a class of early Universe objects that have gained attention as a potential formation route for globular clusters. SIGOs have only recently begun to be studied in the context of molecular hydrogen cooling, which is key to characterizing their structure and evolution. Studying the population-level properties of SIGOs with molecular cooling is important for understanding their potential for collapse and star formation, and central for addressing whether SIGOs can survive to the present epoch. Here, we investigate the evolution of SIGOs before they form stars, using a combination of numerical and analytical analysis. For example, we study various timescales important to the evolution of SIGOs at a population level in the presence of molecular cooling. Revising the previous formulation for the critical density of collapse for SIGOs allows us to show that their prolateness tends to act as an inhibiting factor to collapse. We find that simulated SIGOs are limited by artificial two-body relaxation effects that tend to disperse them, an effect of their limited resolution. We expect that SIGOs in nature will be longer-lived compared to our simulations. Further, the fall-back timescale on which SIGOs fall into nearby dark matter halos, potentially producing a globular-cluster-like system, is frequently longer than their cooling timescale and the collapse timescale on which they shrink through gravity. Therefore, some SIGOs have time to cool and collapse outside of halos despite initially failing to exceed the critical density, even without considering metal line cooling. From this analysis we conclude that SIGOs should form stars outside of halos in non-negligible stream velocity patches in the Universe. △ Less

Submitted 9 January, 2023; v1 submitted 11 August, 2022; originally announced August 2022.

Comments: 18 pages, 11 figures

On the impact of runaway stars on dwarf galaxies with resolved interstellar medium

Authors: Ulrich P. Steinwandel, Greg L. Bryan, Rachel S. Somerville, Christopher C. Hayward, Blakesley Burkhart

Abstract: About ten to 20 percent of massive stars may be kicked out of their natal clusters before exploding as supernovae. These "runaway stars" might play a crucial role in driving galactic outflows and enriching the circumgalactic medium with metals. To study this effect, we carry out high resolution dwarf galaxy simulations that include velocity kicks to massive O/B stars above 8 M$_{\odot}$. We consid… ▽ More About ten to 20 percent of massive stars may be kicked out of their natal clusters before exploding as supernovae. These "runaway stars" might play a crucial role in driving galactic outflows and enriching the circumgalactic medium with metals. To study this effect, we carry out high resolution dwarf galaxy simulations that include velocity kicks to massive O/B stars above 8 M$_{\odot}$. We consider two scenarios, one that adopts a power law velocity distribution for kick velocities, resulting in more stars with high velocity kicks, and a more moderate scenario with a Maxwellian velocity distribution. We explicitly resolve the multi-phase interstellar medium (ISM), and include non-equilibrium cooling and chemistry channels. We adopt a resolved feedback scheme (\textsc{Griffin}) where we sample individual massive stars from an IMF. We follow the lifetime of these stars and add their photoionising radiation, their UV radiation field, and their photoelectric heating rate to the surrounding gas. At the end of their lifetime we explode the massive population as core collapse supernovae (CCSN). In the simulations with runaway massive stars, we add additional (natal) velocity kicks that mimic two and three body interactions that cannot be fully resolved in our simulations. We find that the inclusion of runaway or walkaway star scenarios has an impact on mass, metal, momentum and energy outflows as well as the respective loading factors. We find an increase in mass, metal and momentum loading by a factor of 2-3, whereas we find an increase in the mean energy loading by a factor of 5 in the runaway case and a factor of 3 in the walkaway case. However, we find that the peak values are increased by a factor of up to 10, independent of the adopted velocity kick model. We conclude that the inclusion of runaway stars could have a significant impact on the global outflow properties of dwarf galaxies. △ Less

Submitted 19 May, 2022; originally announced May 2022.

Comments: 18 pages, 13 Figures, 3 tables, submitted to MNRAS, comments welcome :-)

The low redshift Lyman-$α$ Forest as a constraint for models of AGN feedback

Authors: Blakesley Burkhart, Megan Tillman, Alexander B. Gurvich, Simeon Bird, Stephanie Tonnesen, Greg L. Bryan, Lars E. Hernquist, Rachel S. Somerville

Abstract: We study the sensitivity of the $z=0.1$ Lyman-$α$ Forest observables, such as the column density distribution function (CDD), flux PDF, flux power spectrum, and line width distribution, to sub-grid models of active galactic nuclei (AGN) feedback using the Illustris and IllustrisTNG (TNG) cosmological simulations. The two simulations share an identical Ultraviolet Background (UVB) prescription and… ▽ More We study the sensitivity of the $z=0.1$ Lyman-$α$ Forest observables, such as the column density distribution function (CDD), flux PDF, flux power spectrum, and line width distribution, to sub-grid models of active galactic nuclei (AGN) feedback using the Illustris and IllustrisTNG (TNG) cosmological simulations. The two simulations share an identical Ultraviolet Background (UVB) prescription and similar cosmological parameters, but TNG features an entirely reworked AGN feedback model. Due to changes in the AGN radio mode model, the original Illustris simulations have a factor of 2-3 fewer Lyman-$α$ absorbers than TNG at column densities $N_{\rm HI}< 10^{15.5}$ cm$^{-2}$. We compare the simulated forest statistics to UV data from the Cosmic Origins Spectrograph (COS) and find that neither simulation can reproduce the slope of the absorber distribution. Both Illustris and TNG also produce significantly smaller line width distributions than observed in the COS data. We show that TNG is in much better agreement with the observed $z=0.1$ flux power spectrum than Illustris. We explore which statistics can disentangle the effects of AGN feedback from alternative UVB models by rescaling the UVB of Illustris to produce a CDD match to TNG. While this UVB rescaling is degenerate with the effect of AGN feedback on the CDD, the amplitude and shape of the flux PDF and 1D flux power spectrum change in a way distinct from a scaling of the UVB. Our study suggests that the $z=0.1$ Lyman-$α$ forest observables can be used as a diagnostic of AGN feedback models. △ Less

Submitted 18 July, 2022; v1 submitted 20 April, 2022; originally announced April 2022.

Comments: ApJL

Spatially Resolved Gas-phase Metallicity in FIRE-2 Dwarfs: Late-Time Evolution of Metallicity Relations in Simulations with Feedback and Mergers

Authors: Lori E. Porter, Matthew E. Orr, Blakesley Burkhart, Andrew Wetzel, Xiangcheng Ma, Philip F. Hopkins, Andrew Emerick

Abstract: We present an analysis of spatially resolved gas-phase metallicity relations in five dwarf galaxies ($M_{halo} \approx 10^{11} M_\odot$, $M_\star \approx 10^{8.8}-10^{9.6} M_\odot$) from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, which include an explicit model for sub-grid turbulent mixing of metals in gas, near $z\approx 0$, over a period of 1.4 Gyrs,… ▽ More We present an analysis of spatially resolved gas-phase metallicity relations in five dwarf galaxies ($M_{halo} \approx 10^{11} M_\odot$, $M_\star \approx 10^{8.8}-10^{9.6} M_\odot$) from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, which include an explicit model for sub-grid turbulent mixing of metals in gas, near $z\approx 0$, over a period of 1.4 Gyrs, and compare our findings with observations. While these dwarf galaxies represent a diverse sample, we find that all simulated galaxies match the observed mass-metallicity (MZR) and mass-metallicity gradient (MZGR) relations. We note that in all five galaxies, the metallicities are effectively identical between phases of the interstellar medium (ISM), with 95$\%$ being within $\pm$0.1 dex between various ISM phases, including the cold and dense gas ($T < 500$ K and $n_{\rm H} > 1$ cm$^{-3}$), ionized gas (near the H$α$ $T \approx 10^4$ K ridge-line), and nebular regions (ionized gas where the 10 Myr-averaged star formation rate is non-zero). We find that most of the scatter in relative metallicity between cold and dense gas and ionized gas/nebular regions can be attributed to either local starburst events or metal-poor inflows. We also note the presence of a major merger in one of our galaxies, m11e, with a substantial impact on the metallicity distribution in the spatially resolved map, showing two strong metallicity peaks and triggering a starburst in the main galaxy. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: 14 pages, 10 figures, submitted to MNRAS

Gas Accretion Can Drive Turbulence in Galaxies

Authors: John C. Forbes, Razieh Emami, Rachel S. Somerville, Shy Genel, Dylan Nelson, Annalisa Pillepich, Blakesley Burkhart, Greg L. Bryan, Mark R. Krumholz, Lars Hernquist, Stephanie Tonnesen, Paul Torrey, Viraj Pandya, Christopher C. Hayward

Abstract: The driving of turbulence in galaxies is deeply connected with the physics of feedback, star formation, outflows, accretion, and radial transport in disks. The velocity dispersion of gas in galaxies therefore offers a promising observational window into these processes. However, the relative importance of each of these mechanisms remains controversial. In this work we revisit the possibility that… ▽ More The driving of turbulence in galaxies is deeply connected with the physics of feedback, star formation, outflows, accretion, and radial transport in disks. The velocity dispersion of gas in galaxies therefore offers a promising observational window into these processes. However, the relative importance of each of these mechanisms remains controversial. In this work we revisit the possibility that turbulence on galactic scales is driven by the direct impact of accreting gaseous material on the disk. We measure this effect in a disk-like star-forming galaxy in IllustrisTNG, using the high-resolution cosmological magnetohydrodynamical simulation TNG50. We employ Lagrangian tracer particles with a high time cadence of only a few Myr to identify accretion and other events, such as star formation, outflows, and movement within the disk. The energies of particles as they arrive in the disk are measured by stacking the events in bins of time before and after the event. The average effect of each event is measured on the galaxy by fitting explicit models for the kinetic and turbulent energies as a function of time in the disk. These measurements are corroborated by measuring the cross-correlation of the turbulent energy in the different annuli of the disk with other time series, and searching for signals of causality, i.e. asymmetries in the cross-correlation across zero time lag. We find that accretion contributes to the large-scale turbulent kinetic energy even if it is not the dominant driver of turbulence in this $\sim 5 \times 10^{9} M_\odot$ stellar mass galaxy. Extrapolating this finding to a range of galaxy masses, we find that there are regimes where energy from direct accretion may dominate the turbulent energy budget, particularly in disk outskirts, galaxies less massive than the Milky Way, and at redshift $\sim 2$. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: Submitted to AAS Journals, comments welcome

Studying Interstellar Turbulence Driving Scales using the Bispectrum

Authors: Michael J. O'Brien, Blakesley Burkhart, Michael J. Shelley

Abstract: We demonstrate the utility of the bispectrum, the Fourier three-point correlation function, for studying driving scales of magnetohydrodynamic (MHD) turbulence in the interstellar medium. We calculate the bispectrum by implementing a parallelized Monte Carlo direct measurement method, which we have made publicly available. In previous works, the bispectrum has been used to identify non-linear scal… ▽ More We demonstrate the utility of the bispectrum, the Fourier three-point correlation function, for studying driving scales of magnetohydrodynamic (MHD) turbulence in the interstellar medium. We calculate the bispectrum by implementing a parallelized Monte Carlo direct measurement method, which we have made publicly available. In previous works, the bispectrum has been used to identify non-linear scaling correlations and break degeneracies in lower-order statistics like the power spectrum. We find that the bicoherence, a related statistic which measures phase coupling of Fourier modes, identifies turbulence driving scales using density and column density fields. In particular, it shows that the driving scale is phase-coupled to scales present in the turbulent cascade. We also find that the presence of an ordered magnetic field at large-scales enhances phase coupling as compared to a pure hydrodynamic case. We therefore suggest the bispectrum and bicoherence as tools for searching for non-locality for wave interactions in MHD turbulence. △ Less

Submitted 24 March, 2022; originally announced March 2022.

Comments: 12 pages, 9 figures

The CAMELS project: public data release

Authors: Francisco Villaescusa-Navarro, Shy Genel, Daniel Anglés-Alcázar, Lucia A. Perez, Pablo Villanueva-Domingo, Digvijay Wadekar, Helen Shao, Faizan G. Mohammad, Sultan Hassan, Emily Moser, Erwin T. Lau, Luis Fernando Machado Poletti Valle, Andrina Nicola, Leander Thiele, Yongseok Jo, Oliver H. E. Philcox, Benjamin D. Oppenheimer, Megan Tillman, ChangHoon Hahn, Neerav Kaushal, Alice Pisani, Matthew Gebhardt, Ana Maria Delgado, Joyce Caliendo, Christina Kreisch , et al. (22 additional authors not shown)

Abstract: The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present… ▽ More The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present the CAMELS public data release, describing the characteristics of the CAMELS simulations and a variety of data products generated from them, including halo, subhalo, galaxy, and void catalogues, power spectra, bispectra, Lyman-$α$ spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over one thousand catalogues that contain billions of galaxies from CAMELS-SAM: a large collection of N-body simulations that have been combined with the Santa Cruz Semi-Analytic Model. We release all the data, comprising more than 350 terabytes and containing 143,922 snapshots, millions of halos, galaxies and summary statistics. We provide further technical details on how to access, download, read, and process the data at \url{https://camels.readthedocs.io}. △ Less

Submitted 4 January, 2022; originally announced January 2022.

Comments: 18 pages, 3 figures. More than 350 Tb of data from thousands of simulations publicly available at https://www.camel-simulations.org

The Simons Observatory: Galactic Science Goals and Forecasts

Authors: Brandon S. Hensley, Susan E. Clark, Valentina Fanfani, Nicoletta Krachmalnicoff, Giulio Fabbian, Davide Poletti, Giuseppe Puglisi, Gabriele Coppi, Jacob Nibauer, Roman Gerasimov, Nicholas Galitzki, Steve K. Choi, Peter C. Ashton, Carlo Baccigalupi, Eric Baxter, Blakesley Burkhart, Erminia Calabrese, Jens Chluba, Josquin Errard, Andrei V. Frolov, Carlos Hervías-Caimapo, Kevin M. Huffenberger, Bradley R. Johnson, Baptiste Jost, Brian Keating , et al. (9 additional authors not shown)

Abstract: Observing in six frequency bands from 27 to 280 GHz over a large sky area, the Simons Observatory (SO) is poised to address many questions in Galactic astrophysics in addition to its principal cosmological goals. In this work, we provide quantitative forecasts on astrophysical parameters of interest for a range of Galactic science cases. We find that SO can: constrain the frequency spectrum of pol… ▽ More Observing in six frequency bands from 27 to 280 GHz over a large sky area, the Simons Observatory (SO) is poised to address many questions in Galactic astrophysics in addition to its principal cosmological goals. In this work, we provide quantitative forecasts on astrophysical parameters of interest for a range of Galactic science cases. We find that SO can: constrain the frequency spectrum of polarized dust emission at a level of $Δβ_d \lesssim 0.01$ and thus test models of dust composition that predict that $β_d$ in polarization differs from that measured in total intensity; measure the correlation coefficient between polarized dust and synchrotron emission with a factor of two greater precision than current constraints; exclude the non-existence of exo-Oort clouds at roughly 2.9$σ$ if the true fraction is similar to the detection rate of giant planets; map more than 850 molecular clouds with at least 50 independent polarization measurements at 1 pc resolution; detect or place upper limits on the polarization fractions of CO(2-1) emission and anomalous microwave emission at the 0.1% level in select regions; and measure the correlation coefficient between optical starlight polarization and microwave polarized dust emission in $1^\circ$ patches for all lines of sight with $N_{\rm H} \gtrsim 2\times10^{20}$ cm$^{-2}$. The goals and forecasts outlined here provide a roadmap for other microwave polarization experiments to expand their scientific scope via Milky Way astrophysics. △ Less

Submitted 3 November, 2021; originally announced November 2021.

Comments: Submitted to AAS journals. 33 pages, 10 figures

Bursting Bubbles: Feedback from Clustered SNe and the Trade-off Between Turbulence and Outflows

Authors: Matthew E. Orr, Drummond B. Fielding, Christopher C. Hayward, Blakesley Burkhart

Abstract: We present an analytic model for clustered supernovae (SNe) feedback in galaxy disks, incorporating the dynamical evolution of superbubbles formed from spatially overlapping SNe remnants. We propose two realistic outcomes for the evolution of superbubbles in galactic disks: (1) the expansion velocity of the shock front falls below the turbulent velocity dispersion of the ISM in the galaxy disk, wh… ▽ More We present an analytic model for clustered supernovae (SNe) feedback in galaxy disks, incorporating the dynamical evolution of superbubbles formed from spatially overlapping SNe remnants. We propose two realistic outcomes for the evolution of superbubbles in galactic disks: (1) the expansion velocity of the shock front falls below the turbulent velocity dispersion of the ISM in the galaxy disk, whereupon the superbubble stalls and fragments, depositing its momentum entirely within the galaxy disk, or (2) the superbubble grows in size to reach the gas scale height, breaking out of the galaxy disk and driving galactic outflows/fountains. In either case, we find that superbubble breakup/breakout almost always occurs before the last Type-II SN ($\lesssim$40 Myr) in the recently formed star cluster, assuming a standard high-end IMF slope, and scalings between stellar lifetimes and masses. The threshold between these two cases implies a break in the effective strength of feedback in driving turbulence within galaxies, and a resulting change in the scalings of, e.g., star formation rates with gas surface density (the Kennicutt-Schmidt relation) and the star formation efficiency in galaxy disks. △ Less

Submitted 2 May, 2022; v1 submitted 29 September, 2021; originally announced September 2021.

Comments: 24 pages, 11 figures, accepted to ApJ

Bursting Bubbles: Clustered Supernova Feedback in Local and High-redshift Galaxies

Authors: Matthew E. Orr, Drummond B. Fielding, Christopher C. Hayward, Blakesley Burkhart

Abstract: We compare an analytic model for the evolution of supernova-driven superbubbles with observations of local and high-redshift galaxies, and the properties of intact HI shells in local star-forming galaxies. Our model correctly predicts the presence of superwinds in local star-forming galaxies (e.g., NGC 253) and the ubiquity of outflows near $z \sim 2$. We find that high-redshift galaxies may `capt… ▽ More We compare an analytic model for the evolution of supernova-driven superbubbles with observations of local and high-redshift galaxies, and the properties of intact HI shells in local star-forming galaxies. Our model correctly predicts the presence of superwinds in local star-forming galaxies (e.g., NGC 253) and the ubiquity of outflows near $z \sim 2$. We find that high-redshift galaxies may `capture' 20-50\% of their feedback momentum in the dense ISM (with the remainder escaping into the nearby CGM), whereas local galaxies may contain $\lesssim$10\% of their feedback momentum from the central starburst. Using azimuthally averaged galaxy properties, we predict that most superbubbles stall and fragment \emph{within} the ISM, and that this occurs at, or near, the gas scale height. We find a consistent interpretation in the observed HI bubble radii and velocities, and predict that most will fragment within the ISM, and that those able to break-out originate from short dynamical time regions (where the dynamical time is shorter than feedback timescales). Additionally, we demonstrate that models with constant star cluster formation efficiency per Toomre mass are inconsistent with the occurrence of outflows from high-$z$ starbursts and local circumnuclear regions. △ Less

Submitted 21 January, 2022; v1 submitted 29 September, 2021; originally announced September 2021.

Comments: 10 pages, 4 figures, accepted to ApJL

Journal ref: ApJL 924 L28 (2022)

The Effects of Magnetic Fields and Outflow Feedback on the Shape and Evolution of the Density PDF in Turbulent Star-Forming Clouds

Authors: Sabrina M. Appel, Blakesley Burkhart, Vadim A. Semenov, Christoph Federrath, Anna L. Rosen

Abstract: Using a suite of 3D hydrodynamical simulations of star-forming molecular clouds, we investigate how the density probability distribution function (PDF) changes when including gravity, turbulence, magnetic fields, and protostellar outflows and heating. We find that the density PDF is not lognormal when outflows and self-gravity are considered. Self-gravity produces a power-law tail at high densitie… ▽ More Using a suite of 3D hydrodynamical simulations of star-forming molecular clouds, we investigate how the density probability distribution function (PDF) changes when including gravity, turbulence, magnetic fields, and protostellar outflows and heating. We find that the density PDF is not lognormal when outflows and self-gravity are considered. Self-gravity produces a power-law tail at high densities and the inclusion of stellar feedback from protostellar outflows and heating produces significant time-varying deviations from a lognormal distribution at the low densities. The simulation with outflows has an excess of diffuse gas compared to the simulations without outflows, exhibits increased average sonic Mach number, and maintains a slower star formation rate over the entire duration of the run. We study the mass transfer between the diffuse gas in the lognormal peak of the PDF, the collapsing gas in the power-law tail, and the stars. We find that the mass fraction in the power-law tail is constant, such that the stars form out of the power-law gas at the same rate at which the gas from the lognormal part replenishes the power-law. We find that turbulence does not provide significant support in the dense gas associated with the power-law tail. When including outflows and magnetic fields in addition to driven turbulence, the rate of mass transfer from the lognormal to the power-law, and then to the stars, becomes significantly slower, resulting in slower star formation rates and longer depletion times. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 15 pages, 10 figures. Submitted to ApJ. Comments welcome

First extragalactic measurement of the turbulence driving parameter: ALMA observations of the star-forming region N159E in the Large Magellanic Cloud

Authors: Piyush Sharda, Shyam H. Menon, Christoph Federrath, Mark R. Krumholz, James R. Beattie, Katherine E. Jameson, Kazuki Tokuda, Blakesley Burkhart, Roland M. Crocker, Charles J. Law, Amit Seta, Terrance J. Gaetz, Nickolas M. Pingel, Ivo R. Seitenzahl, Hidetoshi Sano, Yasuo Fukui

Abstract: Studying the driving modes of turbulence is important for characterizing the impact of turbulence in various astrophysical environments. The driving mode of turbulence is parameterized by $b$, which relates the width of the gas density PDF to the turbulent Mach number; $b\approx 1/3$, $1$, and $0.4$ correspond to driving that is solenoidal, compressive, and a natural mixture of the two, respective… ▽ More Studying the driving modes of turbulence is important for characterizing the impact of turbulence in various astrophysical environments. The driving mode of turbulence is parameterized by $b$, which relates the width of the gas density PDF to the turbulent Mach number; $b\approx 1/3$, $1$, and $0.4$ correspond to driving that is solenoidal, compressive, and a natural mixture of the two, respectively. In this work, we use high-resolution (sub-pc) ALMA $^{12}$CO ($J$ = $2-1$), $^{13}$CO ($J$ = $2-1$), and C$^{18}$O ($J$ = $2-1$) observations of filamentary molecular clouds in the star-forming region N159E (the Papillon Nebula) in the Large Magellanic Cloud (LMC) to provide the first measurement of turbulence driving parameter in an extragalactic region. We use a non-local thermodynamic equilibrium (NLTE) analysis of the CO isotopologues to construct a gas density PDF, which we find to be largely log-normal in shape with some intermittent features indicating deviations from lognormality. We find that the width of the log-normal part of the density PDF is comparable to the supersonic turbulent Mach number, resulting in $b \approx 0.9$. This implies that the driving mode of turbulence in N159E is primarily compressive. We speculate that the compressive turbulence could have been powered by gravo-turbulent fragmentation of the molecular gas, or due to compression powered by H I flows that led to the development of the molecular filaments observed by ALMA in the region. Our analysis can be easily applied to study the nature of turbulence driving in resolved star-forming regions in the local as well as the high-redshift Universe. △ Less

Submitted 19 October, 2021; v1 submitted 8 September, 2021; originally announced September 2021.

Comments: 14 pages, 7 figures. MNRAS in press

Diagnosing Turbulence in the Neutral and Molecular Interstellar Medium of Galaxies

Authors: Blakesley Burkhart

Abstract: Magnetohydrodynamic (MHD) turbulence is a crucial component of the current paradigms of star formation, dynamo theory, particle transport, magnetic reconnection and evolution of structure in the interstellar medium (ISM) of galaxies. Despite the importance of turbulence to astrophysical fluids, a full theoretical framework based on solutions to the Navier-Stokes equations remains intractable. Obse… ▽ More Magnetohydrodynamic (MHD) turbulence is a crucial component of the current paradigms of star formation, dynamo theory, particle transport, magnetic reconnection and evolution of structure in the interstellar medium (ISM) of galaxies. Despite the importance of turbulence to astrophysical fluids, a full theoretical framework based on solutions to the Navier-Stokes equations remains intractable. Observations provide only limited line-of-sight information on densities, temperatures, velocities and magnetic field strengths and therefore directly measuring turbulence in the ISM is challenging. A statistical approach has been of great utility in allowing comparisons of observations, simulations and analytic predictions. In this review article we address the growing importance of MHD turbulence in many fields of astrophysics and review statistical diagnostics for studying interstellar and interplanetary turbulence. In particular, we will review statistical diagnostics and machine learning algorithms that have been developed for observational data sets in order to obtain information about the turbulence cascade, fluid compressibility (sonic Mach number), and magnetization of fluid (Alfvénic Mach number). These techniques have often been tested on numerical simulations of MHD turbulence, which may include the creation of synthetic observations, and are often formulated on theoretical expectations for compressible magnetized turbulence. We stress the use of multiple techniques, as this can provide a more accurate indication of the turbulence parameters of interest. We conclude by describing several open-source tools for the astrophysical community to use when dealing with turbulence. △ Less

Submitted 3 June, 2021; originally announced June 2021.

Comments: submitted as a review to Publications of the Astronomical Society of the Pacific; comments welcome

The Supersonic Project: SIGOs, a Proposed Progenitor to Globular Clusters, and their Connections to Gravitational Wave Anisotropies

Authors: William Lake, Smadar Naoz, Yeou S. Chiou, Blakesley Burkhart, Federico Marinacci, Mark Vogelsberger, Kyle Kremer

Abstract: Supersonically Induced Gas Objects (SIGOs), are structures with little to no dark matter component predicted to exist in regions of the Universe with large relative velocities between baryons and dark matter at the time of recombination. They have been suggested to be the progenitors of present-day globular clusters. Using simulations, SIGOs have been studied on small scales (around 2 Mpc), where… ▽ More Supersonically Induced Gas Objects (SIGOs), are structures with little to no dark matter component predicted to exist in regions of the Universe with large relative velocities between baryons and dark matter at the time of recombination. They have been suggested to be the progenitors of present-day globular clusters. Using simulations, SIGOs have been studied on small scales (around 2 Mpc), where these relative velocities are coherent. However, it is challenging to study SIGOs using simulations on large scales due to the varying relative velocities at scales larger than a few Mpc. Here, we study SIGO abundances semi-analytically: using perturbation theory, we predict the number density of SIGOs analytically, and compare these results to small-box numerical simulations. We use the agreement between the numerical and analytic calculations to extrapolate the large-scale variation of SIGO abundances over different stream velocities. As a result, we predict similar large-scale variations of objects with high gas densities before reionization that could possibly be observed by JWST. If indeed SIGOs are progenitors of globular clusters, then we expect a similar variation of globular cluster abundances over large scales. Significantly, we find that the expected number density of SIGOs is consistent with observed globular cluster number densities. As a proof-of-concept, and because globular clusters were proposed to be natural formation sites for gravitational wave sources from binary black hole (BBH) mergers, we show that SIGOs should imprint an anisotropy on the gravitational wave signal on the sky, consistent with SIGOs' distribution. △ Less

Submitted 24 August, 2021; v1 submitted 22 April, 2021; originally announced April 2021.

Comments: 17 pages, 8 figures, Accepted to ApJ

The Core Mass Function in the Orion Nebula Cluster Region: What Determines the Final Stellar Masses?

Authors: Hideaki Takemura, Fumitaka Nakamura, Shuo Kong, Héctor G. Arce, John M. Carpenter, Volker Ossenkopf-Okada, Ralf Klessen, Patricio Sanhueza, Yoshito Shimajiri, Takashi Tsukagoshi, Ryohei Kawabe, Shun Ishii, Kazuhito Dobashi, Tomomi Shimoikura, Paul F. Goldsmith, Álvaro Sánchez-Monge, Jens Kauffmann, Thushara Pillai, Paolo Padoan, Adam Ginsberg, Rowan J. Smith, John Bally, Steve Mairs, Jaime E. Pineda, Dariusz C. Lis , et al. (7 additional authors not shown)

Abstract: Applying dendrogram analysis to the CARMA-NRO C$^{18}$O ($J$=1--0) data having an angular resolution of $\sim$ 8", we identified 692 dense cores in the Orion Nebula Cluster (ONC) region. Using this core sample, we compare the core and initial stellar mass functions in the same area to quantify the step from cores to stars. About 22 \% of the identified cores are gravitationally bound. The derived… ▽ More Applying dendrogram analysis to the CARMA-NRO C$^{18}$O ($J$=1--0) data having an angular resolution of $\sim$ 8", we identified 692 dense cores in the Orion Nebula Cluster (ONC) region. Using this core sample, we compare the core and initial stellar mass functions in the same area to quantify the step from cores to stars. About 22 \% of the identified cores are gravitationally bound. The derived core mass function (CMF) for starless cores has a slope similar to Salpeter's stellar initial mass function (IMF) for the mass range above 1 $M_\odot$, consistent with previous studies. Our CMF has a peak at a subsolar mass of $\sim$ 0.1 $M_\odot$, which is comparable to the peak mass of the IMF derived in the same area. We also find that the current star formation rate is consistent with the picture in which stars are born only from self-gravitating starless cores. However, the cores must gain additional gas from the surroundings to reproduce the current IMF (e.g., its slope and peak mass), because the core mass cannot be accreted onto the star with a 100\% efficiency. Thus, the mass accretion from the surroundings may play a crucial role in determining the final stellar masses of stars. △ Less

Submitted 25 February, 2021; originally announced March 2021.

Comments: 13 pages, 5 figures, accepted by ApJL

SIGAME v3: Gas Fragmentation in Post-processing of Cosmological Simulations for More Accurate Infrared Line Emission Modeling

Authors: Karen Pardos Olsen, Blakesley Burkhart, Mordecai-Mark Mac Low, Robin G. Treß, Thomas R. Greve, David Vizgan, Jay Motka, Josh Borrow, Gergö Popping, Romeel Davé, Rowan J. Smith, Desika Narayanan

Abstract: We present an update to the framework called SImulator of GAlaxy Millimeter/submillimeter Emission (SÍGAME). SÍGAME derives line emission in the far-infrared (FIR) for galaxies in particle-based cosmological hydrodynamics simulations by applying radiative transfer and physics recipes via a post-processing step after completion of the simulation. In this version, a new technique is developed to mod… ▽ More We present an update to the framework called SImulator of GAlaxy Millimeter/submillimeter Emission (SÍGAME). SÍGAME derives line emission in the far-infrared (FIR) for galaxies in particle-based cosmological hydrodynamics simulations by applying radiative transfer and physics recipes via a post-processing step after completion of the simulation. In this version, a new technique is developed to model higher gas densities by parametrizing the gas density probability distribution function (PDF) in higher resolution simulations for use as a look-up table, allowing for more adaptive PDFs than in previous work. SÍGAME v3 is tested on redshift z = 0 galaxies drawn from the SIMBA cosmological simulation for eight FIR emission lines tracing vastly different interstellar medium phases. Including dust radiative transfer with SKIRT and high resolution photo-ionization models with Cloudy, this new method is able to self-consistently reproduce observed relations between line luminosity and star formation rate in all cases, except for [NII]122, [NII]205 and [OI]63, the luminosities of which are overestimated by median factors of 1.6, 1.2 and 1.2 dex, respectively. We attribute the remaining disagreement with observations to the lack of precise attenuation of the interstellar light on subgrid scales (<200 pc). △ Less

Submitted 3 December, 2021; v1 submitted 4 February, 2021; originally announced February 2021.

Comments: 31 pages, 15 figures. Accepted for publication in ApJ

Spatially Resolved Star Formation and Inside-out Quenching in the TNG50 Simulation and 3D-HST Observations

Authors: Erica J. Nelson, Sandro Tacchella, Benedikt Diemer, Joel Leja, Lars Hernquist, Katherine E. Whitaker, Rainer Weinberger, Annalisa Pillepich, Dylan Nelson, Bryan A. Terrazas, Rebecca Nevin, Gabriel B. Brammer, Blakesley Burkhart, Rachel Cochrane, Pieter van Dokkum, Benjamin D. Johnson, Lamiya Mowla, Rudiger Pakmor, Rosalind E. Skelton, Joshua Speagle, Volker Springel, Paul Torrey, Mark Vogelsberger, Stijn Wuyts

Abstract: We compare the star forming main sequence (SFMS) -- both integrated and resolved on 1kpc scales -- between the high-resolution TNG50 simulation of IllustrisTNG and observations from the 3D-HST slitless spectroscopic survey at z~1. Contrasting integrated star formation rates (SFRs), we find that the slope and normalization of the star-forming main sequence in TNG50 are quantitatively consistent wit… ▽ More We compare the star forming main sequence (SFMS) -- both integrated and resolved on 1kpc scales -- between the high-resolution TNG50 simulation of IllustrisTNG and observations from the 3D-HST slitless spectroscopic survey at z~1. Contrasting integrated star formation rates (SFRs), we find that the slope and normalization of the star-forming main sequence in TNG50 are quantitatively consistent with values derived by fitting observations from 3D-HST with the Prospector Bayesian inference framework. The previous offsets of 0.2-1dex between observed and simulated main sequence normalizations are resolved when using the updated masses and SFRs from Prospector. The scatter is generically smaller in TNG50 than in 3D-HST for more massive galaxies with M_*>10^10Msun, even after accounting for observational uncertainties. When comparing resolved star formation, we also find good agreement between TNG50 and 3D-HST: average specific star formation rate (sSFR) radial profiles of galaxies at all masses and radii below, on, and above the SFMS are similar in both normalization and shape. Most noteworthy, massive galaxies with M_*>10^10.5Msun, which have fallen below the SFMS due to ongoing quenching, exhibit a clear central SFR suppression, in both TNG50 and 3D-HST. In TNG this inside-out quenching is due to the supermassive black hole (SMBH) feedback model operating at low accretion rates. In contrast, the original Illustris simulation, without this same physical SMBH mechanism, does not reproduce the central SFR profile suppression seen in data. The observed sSFR profiles provide support for the TNG quenching mechanism and how it affects gas on kiloparsec scales in the centers of galaxies. △ Less

Submitted 28 January, 2021; originally announced January 2021.

Comments: Main results in Figs 1 and 3. The TNG50 data is now fully publicly available at https://www.tng-project.org/

Magnetic Field Orientation in Self-Gravitating Turbulent Molecular Clouds

Authors: Lucas Barreto-Mota, Elisabete M. de Gouveia Dal Pino, Blakesley Burkhart, Claudio Melioli, Reinaldo Santos-Lima, Luís H. S. Kadowaki

Abstract: Stars form inside molecular cloud filaments from the competition of gravitational forces with turbulence and magnetic fields. The exact orientation of these filaments with the magnetic fields depends on the strength of these fields, the gravitational potential, and the line-of-sight (LOS) relative to the mean field. To disentangle these effects we employ three-dimensional magnetohydrodynamical num… ▽ More Stars form inside molecular cloud filaments from the competition of gravitational forces with turbulence and magnetic fields. The exact orientation of these filaments with the magnetic fields depends on the strength of these fields, the gravitational potential, and the line-of-sight (LOS) relative to the mean field. To disentangle these effects we employ three-dimensional magnetohydrodynamical numerical simulations that explore a wide range of initial turbulent and magnetic states, i.e., sub-Alfvénic to super-Alfvénic turbulence, with and without gravity. We use histogram of relative orientation (HRO) and the associated projected Rayleigh statistics (PRS) to study the orientation of density and, in order to compare with observations, the integrated density relative to the magnetic field. We find that in sub-Alfvénic systems the initial coherence of the magnetic is maintained inside the cloud and filaments form perpendicular to the field. This trend is not observed in super-Alfvénic models, where the lines are dragged by gravity and turbulence and filaments are mainly aligned to the field. The PRS analysis of integrated maps shows that LOS effects are important only for sub-Alfvénic clouds. When the LOS is perpendicular to the initial field orientation most of the filaments are perpendicular to the projected magnetic field. The inclusion of gravity increases the number of dense structures perpendicular to the magnetic field, reflected as lower values of the PRS for denser regions, regardless of whether the model is sub- or super-Alfvénic. The comparison of our results with observed molecular clouds reveal that most are compatible with sub-Alfvénic models. △ Less

Submitted 12 March, 2021; v1 submitted 8 January, 2021; originally announced January 2021.

Comments: Paper accepted for publication in MNRAS

Turbulence and Particle Acceleration in Radiative Shock Waves in the Cygnus Loop II: Development of Postshock Turbulence

Authors: John C. Raymond, Jonathan D. Slavin, William P. Blair, Igor V. Chilingarian, Blakesley Burkhart, Ravi Sankrit

Abstract: Radiative shock waves in the Cygnus Loop and other supernova remnants show different morphologies in [O III] and Hα emission. We use HST spectra and narrowband images to study the development of turbulence in the cooling region behind a shock on the west limb of the Cygnus Loop. We refine our earlier estimates of shock parameters that were based upon ground-based spectra, including ram pressure, v… ▽ More Radiative shock waves in the Cygnus Loop and other supernova remnants show different morphologies in [O III] and Hα emission. We use HST spectra and narrowband images to study the development of turbulence in the cooling region behind a shock on the west limb of the Cygnus Loop. We refine our earlier estimates of shock parameters that were based upon ground-based spectra, including ram pressure, vorticity and magnetic field strength. We apply several techniques, including Fourier power spectra and the Rolling Hough Transform, to quantify the shape of the rippled shock front as viewed in different emission lines. We assess the relative importance of thermal instabilities, the thin shell instability, upstream density variations, and upstream magnetic field variations in producing the observed structure. △ Less

Submitted 27 October, 2020; v1 submitted 24 October, 2020; originally announced October 2020.

Comments: 40 pages, 22 Figures. Accepted for publication in ApJ

Classification of Magnetohydrodynamic Simulations using Wavelet Scattering Transforms

Authors: Andrew K. Saydjari, Stephen K. N. Portillo, Zachary Slepian, Sule Kahraman, Blakesley Burkhart, Douglas P. Finkbeiner

Abstract: The complex interplay of magnetohydrodynamics, gravity, and supersonic turbulence in the interstellar medium (ISM) introduces non-Gaussian structure that can complicate comparison between theory and observation. We show that the Wavelet Scattering Transform (WST), in combination with linear discriminant analysis (LDA), is sensitive to non-Gaussian structure in 2D ISM dust maps. WST-LDA classifies… ▽ More The complex interplay of magnetohydrodynamics, gravity, and supersonic turbulence in the interstellar medium (ISM) introduces non-Gaussian structure that can complicate comparison between theory and observation. We show that the Wavelet Scattering Transform (WST), in combination with linear discriminant analysis (LDA), is sensitive to non-Gaussian structure in 2D ISM dust maps. WST-LDA classifies magnetohydrodynamic (MHD) turbulence simulations with up to a 97\% true positive rate in our testbed of 8 simulations with varying sonic and Alfvénic Mach numbers. We present a side-by-side comparison with two other methods for non-Gaussian characterization, the Reduced Wavelet Scattering Transform (RWST) and the 3-Point Correlation Function (3PCF). We also demonstrate the 3D-WST-LDA and apply it to classification of density fields in position-position-velocity (PPV) space, where density correlations can be studied using velocity coherence as a proxy. WST-LDA is robust to common observational artifacts, such as striping and missing data, while also sensitive enough to extract the net magnetic field direction for sub-Alfvénic turbulent density fields. We include a brief analysis of the effect of point spread functions and image pixelization on 2D-WST-LDA applied to density fields, which informs the future goal of applying WST-LDA to 2D or 3D all-sky dust maps to extract hydrodynamic parameters of interest. △ Less

Submitted 22 October, 2020; originally announced October 2020.

Comments: 21 pages, 13 figures, submitted to ApJ

The Catalogue for Astrophysical Turbulence Simulations (CATS)

Authors: B. Burkhart, S. Appel, S. Bialy, J. Cho, A. J. Christensen, D. Collins, C. Federrath, D. Fielding, D. Finkbeiner, A. S. Hill, J. C. Ibanez-Mejia, M. R. Krumholz, A. Lazarian, M. Li, P. Mocz, M. -M. Mac Low, J. Naiman, S. K. N. Portillo, B. Shane, Z. Slepian, Y. Yuan

Abstract: Turbulence is a key process in many fields of astrophysics. Advances in numerical simulations of fluids over the last several decades have revolutionized our understanding of turbulence and related processes such as star formation and cosmic ray propagation. However, data from numerical simulations of astrophysical turbulence are often not made public. We introduce a new simulation-oriented databa… ▽ More Turbulence is a key process in many fields of astrophysics. Advances in numerical simulations of fluids over the last several decades have revolutionized our understanding of turbulence and related processes such as star formation and cosmic ray propagation. However, data from numerical simulations of astrophysical turbulence are often not made public. We introduce a new simulation-oriented database for the astronomical community: The Catalogue for Astrophysical Turbulence Simulations (CATS), located at www.mhdturbulence.com. CATS includes magnetohydrodynamic (MHD) turbulent box simulation data products generated by the public codes athena++, arepo, enzo, and flash. CATS also includes several synthetic observational data sets, such as turbulent HI data cubes. We also include measured power spectra and 3-point correlation functions from some of these data. We discuss the importance of open source statistical and visualization tools for the analysis of turbulence simulations such as those found in CATS. △ Less

Submitted 21 October, 2020; originally announced October 2020.

Comments: Accepted to ApJ

The CAMELS project: Cosmology and Astrophysics with MachinE Learning Simulations

Authors: Francisco Villaescusa-Navarro, Daniel Anglés-Alcázar, Shy Genel, David N. Spergel, Rachel S. Somerville, Romeel Dave, Annalisa Pillepich, Lars Hernquist, Dylan Nelson, Paul Torrey, Desika Narayanan, Yin Li, Oliver Philcox, Valentina La Torre, Ana Maria Delgado, Shirley Ho, Sultan Hassan, Blakesley Burkhart, Digvijay Wadekar, Nicholas Battaglia, Gabriella Contardo, Greg L. Bryan

Abstract: We present the Cosmology and Astrophysics with MachinE Learning Simulations --CAMELS-- project. CAMELS is a suite of 4,233 cosmological simulations of $(25~h^{-1}{\rm Mpc})^3$ volume each: 2,184 state-of-the-art (magneto-)hydrodynamic simulations run with the AREPO and GIZMO codes, employing the same baryonic subgrid physics as the IllustrisTNG and SIMBA simulations, and 2,049 N-body simulations.… ▽ More We present the Cosmology and Astrophysics with MachinE Learning Simulations --CAMELS-- project. CAMELS is a suite of 4,233 cosmological simulations of $(25~h^{-1}{\rm Mpc})^3$ volume each: 2,184 state-of-the-art (magneto-)hydrodynamic simulations run with the AREPO and GIZMO codes, employing the same baryonic subgrid physics as the IllustrisTNG and SIMBA simulations, and 2,049 N-body simulations. The goal of the CAMELS project is to provide theory predictions for different observables as a function of cosmology and astrophysics, and it is the largest suite of cosmological (magneto-)hydrodynamic simulations designed to train machine learning algorithms. CAMELS contains thousands of different cosmological and astrophysical models by way of varying $Ω_m$, $σ_8$, and four parameters controlling stellar and AGN feedback, following the evolution of more than 100 billion particles and fluid elements over a combined volume of $(400~h^{-1}{\rm Mpc})^3$. We describe the simulations in detail and characterize the large range of conditions represented in terms of the matter power spectrum, cosmic star formation rate density, galaxy stellar mass function, halo baryon fractions, and several galaxy scaling relations. We show that the IllustrisTNG and SIMBA suites produce roughly similar distributions of galaxy properties over the full parameter space but significantly different halo baryon fractions and baryonic effects on the matter power spectrum. This emphasizes the need for marginalizing over baryonic effects to extract the maximum amount of information from cosmological surveys. We illustrate the unique potential of CAMELS using several machine learning applications, including non-linear interpolation, parameter estimation, symbolic regression, data generation with Generative Adversarial Networks (GANs), dimensionality reduction, and anomaly detection. △ Less

Submitted 15 August, 2021; v1 submitted 1 October, 2020; originally announced October 2020.

Comments: 34 pages, 19 figures, Matches published version. CAMELS webpage at https://www.camel-simulations.org

Journal ref: ApJ, 915, 71, (2021)

Efficient early stellar feedback can suppress galactic outflows by reducing supernova clustering

Authors: Matthew C. Smith, Greg L. Bryan, Rachel S. Somerville, Chia-Yu Hu, Romain Teyssier, Blakesley Burkhart, Lars Hernquist

Abstract: We present a novel set of stellar feedback models, implemented in the moving-mesh code Arepo, designed for galaxy formation simulations with near-parsec (or better) resolution. These include explicit sampling of stars from the IMF, allowing feedback to be linked to individual massive stars, an improved method for the modelling of H II regions, photoelectric heating from a spatially varying FUV fie… ▽ More We present a novel set of stellar feedback models, implemented in the moving-mesh code Arepo, designed for galaxy formation simulations with near-parsec (or better) resolution. These include explicit sampling of stars from the IMF, allowing feedback to be linked to individual massive stars, an improved method for the modelling of H II regions, photoelectric heating from a spatially varying FUV field and supernova feedback. We perform a suite of 32 simulations of isolated $M_\mathrm{vir} = 10^{10}\,\mathrm{M_\odot}$ galaxies with a baryonic mass resolution of $20\,\mathrm{M_\odot}$ in order to study the non-linear coupling of the different feedback channels. We find that photoionization and supernova feedback are both independently capable of regulating star formation to the same level, while photoelectric heating is inefficient. Photoionization produces a considerably smoother star formation history than supernovae. When all feedback channels are combined, the additional suppression of star formation rates is minor. However, outflow rates are substantially reduced relative to the supernova only simulations. We show that this is directly caused by a suppression of supernova clustering by the photoionization feedback, disrupting star forming clouds prior to the first supernovae. We demonstrate that our results are robust to variations of our star formation prescription, feedback models and the gas fraction of the disk. Our results also imply that the burstiness of star formation and the mass loading of outflows may be overestimated if the adopted star particle mass is considerably larger than the mass of individual stars because this imposes a minimum cluster size. △ Less

Submitted 20 August, 2021; v1 submitted 23 September, 2020; originally announced September 2020.

Comments: Published in MNRAS, 37 pages (including appendices), 23 figures. Updated for consistency with published version

Journal ref: 2021, MNRAS, Volume 506, Issue 3, pp.3882-3915

The Supersonic Project: To cool or not to cool Supersonically Induced Gas Objects (SIGOs)?

Authors: Yeou S. Chiou, Smadar Naoz, Blakesley Burkhart, Federico Marinacci, Mark Vogelsberger

Abstract: Supersonically Induced Gas Objects (SIGOs) primarily form in the early Universe, outside of dark matter halos due to the presence of a relative stream velocity between baryons and dark matter. These structures may be the progenitors of globular clusters. Since SIGOs are made out of pristine gas, we investigate the effect of atomic cooling on their properties. We run a suite of simulations by using… ▽ More Supersonically Induced Gas Objects (SIGOs) primarily form in the early Universe, outside of dark matter halos due to the presence of a relative stream velocity between baryons and dark matter. These structures may be the progenitors of globular clusters. Since SIGOs are made out of pristine gas, we investigate the effect of atomic cooling on their properties. We run a suite of simulations by using the moving-mesh code {\sc arepo}, with and without baryon-dark matter relative velocity and with and without the effects of atomic cooling. We show that SIGO's density, temperature, and prolateness are determined by gravitational interactions rather than cooling. The cold gas fraction in SIGOs is much higher than that of dark matter halos. Specifically, we show that the SIGO's characteristic low temperature and extreme high gas density forges a nurturing ground for the earliest star formation sites. △ Less

Submitted 8 November, 2020; v1 submitted 6 August, 2020; originally announced August 2020.

Comments: 12 pages, 9 figures, accepted by ApJ, expanded motivation and updated figures

Understanding biases in measurements of molecular cloud kinematics using line emission

Authors: Yuxuan Yuan, Mark R. Krumholz, Blakesley Burkhart

Abstract: Molecular line observations using a variety of tracers are often used to investigate the kinematic structure of molecular clouds. However, measurements of cloud velocity dispersions with different lines, even in the same region, often yield inconsistent results. The reasons for this disagreement are not entirely clear since molecular line observations are subject to a number of biases. In this pap… ▽ More Molecular line observations using a variety of tracers are often used to investigate the kinematic structure of molecular clouds. However, measurements of cloud velocity dispersions with different lines, even in the same region, often yield inconsistent results. The reasons for this disagreement are not entirely clear since molecular line observations are subject to a number of biases. In this paper, we untangle and investigate various factors that drive linewidth measurement biases by constructing synthetic position-position-velocity cubes for a variety of tracers from a suite of self-gravitating magnetohydrodynamic simulations of molecular clouds. We compare linewidths derived from synthetic observations of these data cubes to the true values in the simulations. We find that differences in linewidth as measured by different tracers are driven by a combination of density-dependent excitation, whereby tracers that are sensitive to higher densities sample smaller regions with smaller velocity dispersions, opacity broadening, especially for highly optically thick tracers such as CO, and finite resolution and sensitivity, which suppress the wings of emission lines. We find that, at fixed signal to noise ratio, three commonly-used tracers, the J=4->3 line of CO, the J=1->0 line of C^{18}O and the (1,1) inversion transition of NH_3, generally offer the best compromise between these competing biases, and produce estimates of the velocity dispersion that reflect the true kinematics of a molecular cloud to an accuracy of about 10% regardless of the cloud magnetic field strengths, evolutionary state, or orientations of the line of sight relative to the magnetic field. △ Less

Submitted 27 July, 2020; originally announced July 2020.

Comments: 18 pages, 16 figures

First results from SMAUG: The need for preventative stellar feedback and improved baryon cycling in semi-analytic models of galaxy formation

Authors: Viraj Pandya, Rachel S. Somerville, Daniel Anglés-Alcázar, Christopher C. Hayward, Greg L. Bryan, Drummond B. Fielding, John C. Forbes, Blakesley Burkhart, Shy Genel, Lars Hernquist, Chang-Goo Kim, Stephanie Tonnesen, Tjitske Starkenburg

Abstract: Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same… ▽ More Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same halo merger trees, with an emphasis on the global mass flow cycle. Our study includes 13 halos spanning low-mass dwarfs (M_vir~10^10 M_sun at z=0), intermediate-mass dwarfs (M_vir~10^11 M_sun) and Milky Way-mass galaxies (M_vir~10^12 M_sun). The SAM and FIRE-2 predictions agree relatively well with each other in terms of stellar and interstellar mass, but differ dramatically on circumgalactic mass (the SAM is lower than FIRE-2 by ~3 orders of magnitude for dwarfs). Strikingly, the SAM predicts higher gas accretion rates for dwarfs compared to FIRE-2 by factors of ~10-100, and this is compensated for with higher mass outflow rates in the SAM. We argue that the most severe model discrepancies are caused by the lack of preventative stellar feedback and the assumptions for halo gas cooling and recycling in the SAM. As a first step towards resolving these model tensions, we present a simple yet promising new preventative stellar feedback model in which the energy carried by supernova-driven winds is allowed to heat some fraction of gas outside of halos to at least the virial temperature such that accretion is suppressed. △ Less

Submitted 30 October, 2020; v1 submitted 29 June, 2020; originally announced June 2020.

Comments: Accepted to ApJ with minor revisions. More information on the SMAUG project and the other first result papers can be found here: https://www.simonsfoundation.org/flatiron/center-for-computational-astrophysics/galaxy-formation/smaug/papersplash1