Search

The [CII] 157.74 $\mu$m fine structure transition is one of the brightest and most well-studied emission lines in the far-infrared (FIR), produced in the interstellar medium (ISM) of galaxies. We study its properties in sub-pc resolution hydrodynamical simulations for an ISM patch with gas surface density of $\Sigma_{\rm{g}}=10\;M_{\odot}\;\rm{pc}^{-2}$, coupled with time-dependent chemistry, far-ultraviolet (FUV) dust and gas shielding, star formation, photoionization and supernova (SN) feedback, and full line-radiative transfer. We find a [CII]-to-H$_2$ conversion factor that scales weakly with metallicity $X_{\rm{[CII]}}=6.31\times 10^{19} \;Z^{\prime\;0.17}\; \rm{cm}^{-2}\;(\rm{K}\;\rm{km}\;\rm{s}^{-1})^{-1}$, where $Z^{\prime}$ is the normalized metallicity relative to solar. {The majority of [CII] originates from atomic gas with hydrogen number density $n\sim 10~{\rm cm^{-3}}$.} The [CII] line intensity positively correlates with the star formation rate (SFR), with a normalization factor that scales linearly with metallicity. We find that this is broadly consistent with $z\sim0$ observations. As such, [CII] is a good SFR tracer even in metal-poor environments where molecular lines might be undetectable. Resolving the clumpy structure of the dense ($n=10-10^3\;\rm{cm}^{-3}$) interstellar medium (ISM) is important as it dominates [CII] 157.74 $\mu$m emission. We compare our full radiative transfer computation with the optically-thin limit and find that the [CII] line becomes marginally optically thick only at super-solar metallicity for our assumed gas surface density., Comment: Accepted for publication in ApJ

In the current era of high-z galaxy discovery with JWST and ALMA, our ability to study the stellar populations and ISM conditions in a diverse range of galaxies at Cosmic Dawn has rapidly improved. At the same time, the need to understand the current limitations in modeling galaxy formation processes and physical properties in order to interpret these observations is critical. Here, we study the challenges in modeling galaxy dust temperatures, both in the context of forward modeling galaxy spectral properties from a hydrodynamical simulation and via backwards modeling galaxy physical properties from mock observations of far-infrared dust emission. Using the Simba model for galaxy formation combined with Powderday radiative transfer, we can accurately predict the evolution of dust at high redshift, though several aspects of the model are essentially free parameters (dust composition, sub-resolution dust in star-forming regions) that dull the predictive power of the model dust temperature distributions. We also highlight the uncertainties in the backwards modeling methods, where we find the commonly used models and assumptions to fit FIR SEDs and infer dust temperatures (e.g., single temperature, optically thin modified blackbody) largely fail to capture the complexity of high-z dusty galaxies. We caution that conclusions inferred from both simulations -- limited by resolution and post-processing techniques -- and observations -- limited by sparse data and simplistic model parameterizations -- are susceptible to unique and nuanced uncertainties that can limit the usefulness of current high-z dust measurements., Comment: 20 pages main text, 12 figures. Accepted to ApJ, newest version

We present a new methodology for simulating mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs) in galaxy evolution simulations. To do this, we combine theoretical models of PAH emission features as they respond to varying interstellar radiation fields, grain size distributions, and ionization states with a new on-the-fly model for dust evolution in hydrodynamic galaxy simulations. We apply these models to 3 idealized arepo galaxy evolution simulations within the smuggle physics framework. We use these simulations to develop numerical experiments investigating the buildup of PAH masses and luminosities in galaxies in idealized analogs of the Milky Way, a dwarf galaxy, and starburst disk. Our main results follow. Galaxies with high specific star formation rates have increased feedback energy per unit mass, and are able to efficiently shatter dust grains, driving up the fraction of ultra small grains. At the same time, in our model large radiation fields per unit gas density convert aliphatic grains into aromatics. The fraction of dust grains in the form of PAHs (q_PAH) can be understood as a consequence of these processes, and in our model PAHs form primarily from interstellar processing (shattering) of larger grains rather than from the growth of smaller grains. We find that the hardness of the radiation field plays a larger role than variations in the grain size distribution in setting the total integrated PAH luminosities, though cosmological simulations are necessary to fully investigate the complex interplay of processes that drive PAH band luminosities in galaxies. Finally, we highlight feature PAH strength variations, cautioning against the usage of emission templates with constant feature strength ratios., Comment: Submitted to ApJ; comments welcome

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., Comment: Published MNRAS, 27 pages, 17 figures; updated for consistency with journal version

Nearby dwarf irregular galaxies are ideal laboratories for studying the interstellar medium (ISM) at low metallicity, which is expected to be common for galaxies at very high redshift being observed by the James Webb Space Telescope. We present the first high-resolution ($\sim 0.2$~pc) hydrodynamical simulations of an isolated low-metallicity ($0.1~Z_\odot$) dwarf galaxy coupled with a time-dependent chemistry network and a dust evolution model where dust is locally produced and destroyed by various processes. To accurately model carbon monoxide (CO), we post-process the simulations with a detailed chemistry network including the time-dependent effect of molecular hydrogen (H$_2$). Our model successfully reproduces the observed star formation rate and CO(1-0) luminosity ($L_{\rm CO}$). We find that dust growth in dense gas is required to reproduce the observed $L_{\rm CO}$ as otherwise CO would be completely photodissociated. In contrast, the H$_2$ abundance is extremely small and is insensitive to dust growth, leading to a CO-to-H$_2$ conversion factor that is only slightly higher than the Milky Way value despite the low metallicity. An observationally inferred dust-to-gas ratio is thus underestimated if adopting the metallicity-dependent CO-to-H$_2$ conversion factor. The newly-produced dust in dense gas mixes with the ISM through supernova feedback without being completely destroyed by sputtering, which leads to galactic outflows 20% - 50% dustier than the ISM, providing a possible source for intergalactic dust., Comment: ApJ accepted. 19 pages, 10 figures

We present new GRIFFIN project hydrodynamical simulations that model the formation of galactic star cluster populations in low-metallicity ($Z=0.00021$) dwarf galaxies, including radiation, supernova and stellar wind feedback of individual massive stars. In the simulations, stars are sampled from the stellar initial mass function (IMF) down to the hydrogen burning limit of $0.08$ M$_\odot$. Mass conservation is enforced within a radius of $1$ pc for the formation of massive stars. We find that massive stars are preferentially found in star clusters and follow a correlation set at birth between the highest initial stellar mass and the star cluster mass that differs from pure stochastic IMF sampling. With a fully sampled IMF, star clusters lose mass in the galactic tidal field according to mass-loss rates observed in nearby galaxies. Of the released stellar feedback, $60\%$ of the supernova material and up to $35\%$ of the wind material reside either in the hot interstellar medium (ISM) or in gaseous, metal enriched outflows. While stellar winds (instantaneously) and supernovae (delayed) start enriching the ISM right after the first massive stars form, the formation of supernova-enriched stars and star clusters is significantly delayed (by $>50$ Myr) compared to the formation of stars and star clusters enriched by stellar winds. Overall, supernova ejecta dominate the enrichment by mass, while the number of enriched stars is determined by continuous stellar winds. These results present a concept for the formation of chemically distinct populations of stars in bound star clusters, reminiscent of multiple populations in globular clusters., Comment: 26 pages, 23 figures. Accepted for publication in MNRAS

Determining the atomic and molecular content of the interstellar medium (ISM) as a function of environmental parameters is of fundamental importance to understand the star-formation process across the epochs. Although there exist various three-dimensional hydro-chemical codes modelling the ISM at different scales and redshifts, they are computationally expensive and inefficient for studies over a large parameter space. Building on our earlier approach, we present PDFchem, a novel algorithm that models the cold ISM at moderate and large scales using functions connecting the quantities of the local ($A_{\rm V,eff}$) and the observed ($A_{\rm V,obs}$) visual extinctions, and the local number density, $n_{\rm H}$, with probability density functions (PDF) of $A_{\rm V,obs}$ on cloud scales typically tens-to-hundreds of pc as an input. For any given $A_{\rm V,obs}$-PDF, sampled with thousands of clouds, the algorithm instantly computes the average abundances of the most important species (HI, H$_2$, CII, CI, CO, OH, OH$^+$, H$_2$O$^+$, CH, HCO$^+$) and performs radiative transfer calculations to estimate the average emission of the most commonly observed lines ([CII]~$158\mu$m, both [CI] fine-structure lines and the first five rotational transitions of $^{12}$CO). We examine two $A_{\rm V,obs}$-PDFs corresponding to a non star-forming and a star-forming ISM region, under a variety of environmental parameters combinations. These cover FUV intensities in the range of $\chi/\chi_0=10^{-1}-10^3$, cosmic-ray ionization rates in the range of $\zeta_{\rm CR}=10^{-17}-10^{-13}\,{\rm s}^{-1}$ and metallicities in the range of $Z=0.1-2\,{\rm Z}_{\odot}$. PDFchem is fast, easy to use, reproduces the PDR quantities of the time-consuming hydrodynamical models and can be used directly with observed data to understand the evolution of the cold ISM chemistry., Comment: 33 pages, 27 figures. Accepted in MNRAS. The algorithm can be found in: https://github.com/tbisbas/PDFchem. Comments welcome!

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., Comment: accepted version by ApJ (including a new simulation in Appendix B suggested by the referee)

We present [CII] synthetic observations of smoothed particle hydrodynamics (SPH) simulations of a dwarf galaxy merger. The merging process varies the star-formation rate by more than three orders of magnitude. Several star clusters are formed, the feedback of which disperses and unbinds the dense gas through expanding HII regions and supernova (SN) explosions. For galaxies with properties similar to the modelled ones, we find that the [CII] emission remains optically thin throughout the merging process. We identify the Warm Neutral Medium ($3<\log T_{\rm gas}<4$ with $\chi_{\rm HI}>2\chi_{\rm H2}$) to be the primary source of [CII] emission ($\sim58\%$ contribution), although at stages when the HII regions are young and dense (during star cluster formation or SNe in the form of ionized bubbles) they can contribute $\gtrsim50\%$ to the total [CII] emission. We find that the [CII]/FIR ratio decreases due to thermal saturation of the [CII] emission caused by strong FUV radiation fields emitted by the massive star clusters, leading to a [CII]-deficit medium. We investigate the [CII]-SFR relation and find an approximately linear correlation which agrees well with observations, particularly those from the Dwarf Galaxy Survey. Our simulation reproduces the observed trends of [CII]/FIR versus $\Sigma_{\rm SFR}$ and $\Sigma_{\rm FIR}$, and it agrees well with the Kennicutt relation of SFR-FIR luminosity. We propose that local peaks of [CII] in resolved observations may provide evidence for ongoing massive cluster formation., Comment: To appear in ApJ. 24 pages, 17 figures. Comments welcome!

We study the CO(1-0)-to-H$_2$ conversion factor ($X_{\rm CO}$) and the line ratio of CO(2-1)-to-CO(1-0) ($R_{21}$) across a wide range of metallicity ($0.1 \leq Z/Z_\odot \leq 3$) in high-resolution (~0.2 pc) hydrodynamical simulations of a self-regulated multiphase interstellar medium. We construct synthetic CO emission maps via radiative transfer and systematically vary the "observational" beam size to quantify the scale dependence. We find that the kpc-scale $X_{\rm CO}$ can be over-estimated at low $Z$ if assuming steady-state chemistry or assuming that the star-forming gas is H$_2$-dominated. On parsec scales, $X_{\rm CO}$ varies by orders of magnitude from place to place, primarily driven by the transition from atomic carbon to CO. The pc-scale $X_{\rm CO}$ drops to the Milky Way value of $2\times 10^{20}\ {\rm cm^{-2}~(K~km~s^{-1})^{-1}}$ once dust shielding becomes effective, independent of $Z$. The CO lines become increasingly optically thin at lower $Z$, leading to a higher $R_{21}$. Most cloud area is filled by diffuse gas with high $X_{\rm CO}$ and low $R_{21}$, while most CO emission originates from dense gas with low $X_{\rm CO}$ and high $R_{21}$. Adopting a constant $X_{\rm CO}$ strongly over- (under-)estimates H$_2$ in dense (diffuse) gas. The line intensity negatively (positively) correlates with $X_{\rm CO}$ ($R_{21}$) as it is a proxy of column density (volume density). On large scales, $X_{\rm CO}$ and $R_{21}$ are dictated by beam averaging, and they are naturally biased towards values in dense gas. Our predicted $X_{\rm CO}$ is a multivariate function of $Z$, line intensity, and beam size, which can be used to more accurately infer the H$_2$ mass., Comment: Published in ApJ (typos in the abstract fixed). Code for interpolating Lagrangian (particle) data onto an adaptive mesh (and auxiliary scripts for RADMC-3D) available at https://github.com/huchiayu/ParticleGridMapper.jl

We study the metallicity dependence of the H/H$_2$ and C$^+$/C/CO distributions in a self-regulated interstellar medium (ISM) across a broad range of metallicities ($0.1 < Z/Z_\odot < 3$). To this end, we conduct high-resolution (particle mass of $1\ {\rm M_\odot}$) hydrodynamical simulations coupled with a time-dependent H$_2$ chemistry network. The results are then post-processed with an accurate chemistry network to model the associated C$^+$/C/CO abundances, based on the time-dependent non-steady-state (``non-equilibrium'') H$_2$ abundances. We find that the time-averaged star formation rate and the ISM structure are insensitive to metallicity. The column densities relevant for molecular shielding appear correlated with the volume densities in gravitationally unstable gas. As metallicity decreases, H$_2$ progressively deviates from steady state (``equilibrium'') and shows shallow abundance profiles until they sharply truncate at the photodissociation fronts. In contrast, the CO profile is sharp and controlled by photodissociation as CO quickly reaches steady state. We construct effective one-dimensional cloud models that successfully capture the time-averaged chemical distributions in simulations. At low metallicities, the steady-state model significantly overestimates the abundance of H$_2$ in the diffuse medium. The overestimated H$_2$, however, has little impact on CO. Consequently, the mass fraction of CO-dark H$_2$ gas is significantly lower than what a fully steady-state model predicts. The mass ratios of H$_2$/C$^+$ and H$_2$/C both show a weaker dependence on $Z^{\prime}$ than H$_2$/CO, which potentially indicates that C$^+$ and C could be alternative tracers for H$_2$ at low $Z^{\prime}$ in terms of mass budget. Our chemistry code for post-processing is publicly available., Comment: ApJ accepted, 32 pages (main text 24 pages), chemistry code available at https://github.com/huchiayu/AstroChemistry.jl

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., Comment: Published in MNRAS, 37 pages (including appendices), 23 figures. Updated for consistency with published version

We analyze the three-dimensional shapes and kinematics of the young star cluster population forming in a high-resolution GRIFFIN project simulation of a metal-poor dwarf galaxy starburst. The star clusters, which follow a power-law mass distribution, form from the cold ISM phase with an IMF sampled with individual stars down to 4 solar masses at sub-parsec spatial resolution. Massive stars and their important feedback mechanisms are modelled in detail. The simulated clusters follow a surprisingly tight relation between the specific angular momentum and mass with indications of two sub-populations. Massive clusters ($M_\mathrm{cl}\gtrsim 3\times 10^4 M_{\odot})$ have the highest specific angular momenta at low ellipticities ($\epsilon\sim 0.2$) and show alignment between their shapes and rotation. Lower mass clusters have lower specific angular momenta with larger scatter, show a broader range of elongations, and are typically misaligned indicating that they are not shaped by rotation. The most massive clusters $(M \gtrsim 10^5\,M_{\odot})$ accrete gas and proto-clusters from a $ \lesssim 100\,\rm pc$ scale local galactic environment on a $t \lesssim 10\,\rm Myr$ timescale, inheriting the ambient angular momentum properties. Their two-dimensional kinematic maps show ordered rotation at formation, up to $v \sim 8.5\,\rm km s^{-1}$, consistent with observed young massive clusters and old globular clusters, which they might evolve into. The massive clusters have angular momentum parameters $\lambda_R\lesssim 0.5$ and show Gauss-Hermite coefficients $h_3$ that are anti-correlated with the velocity, indicating asymmetric line-of-sight velocity distributions as a signature of a dissipative formation process., Comment: 21 pages, accepted for publication in ApJ

We perform a-priori validation tests of subgrid-scale (SGS) models for the turbulent transport of momentum, energy and passive scalars. To this end, we conduct two sets of high-resolution hydrodynamical simulations with a Lagrangian code: an isothermal turbulent box with rms Mach number of 0.3, 2 and 8, and the classical wind tunnel where a cold cloud traveling through a hot medium gradually dissolves due to fluid instabilities. Two SGS models are examined: the eddy diffusivity (ED) model wildly adopted in astrophysical simulations and the "gradient model" due to Clark et al. (1979). We find that both models predict the magnitude of the SGS terms equally well (correlation coefficient > 0.8). However, the gradient model provides excellent predictions on the orientation and shape of the SGS terms while the ED model predicts poorly on both, indicating that isotropic diffusion is a poor approximation of the instantaneous turbulent transport. The best-fit coefficient of the gradient model is in the range of [0.16, 0.21] for the momentum transport, and the turbulent Schmidt number and Prandtl number are both close to unity, in the range of [0.92, 1.15]., Comment: ApJ accepted; analysis code available at https://github.com/huchiayu/Lapriori.jl

We describe a population of young star clusters (SCs) formed in a hydrodynamical simulation of a gas-rich dwarf galaxy merger resolved with individual massive stars at sub-parsec spatial resolution. The simulation is part of the GRIFFIN (Galaxy Realizations Including Feedback From INdividual massive stars) project. The star formation environment during the simulation spans seven orders of magnitude in gas surface density and thermal pressure, and the global star formation rate surface density ($\Sigma_\mathrm{SFR}$) varies by more than three orders of magnitude during the simulation. Young SCs more massive than $M_{\mathrm{*,cl}}\sim 10^{2.5}\,M_{\odot}$ form along a mass function with a power-law index $\alpha\sim-1.7$ ($\alpha\sim-2$ for $M_{\mathrm{*,cl}}\gtrsim10^{3}\,M_{\odot}$) at all merger phases, while the normalization and the highest SC masses (up to $\sim 10^6 M_{\odot}$) correlate with $\Sigma_\mathrm{SFR}$. The cluster formation efficiency varies from $\Gamma\sim20\%$ in early merger phases to $\Gamma\sim80\%$ at the peak of the starburst and is compared to observations and model predictions. The massive SCs ($\gtrsim10^4\,M_{\odot}$) have sizes and mean surface densities similar to observed young massive SCs. Simulated lower mass clusters appear slightly more concentrated than observed. All SCs form on timescales of a few Myr and lose their gas rapidly resulting in typical stellar age spreads between $\sigma\sim0.1-2$ Myr ($1\sigma$), consistent with observations. The age spreads increase with cluster mass, with the most massive cluster ($\sim10^6\, M_{\odot}$) reaching a spread of $5\, \mathrm{Myr}$ once its hierarchical formation finishes. Our study shows that it is now feasible to investigate the SC population of entire galaxies with novel high-resolution numerical simulations., Comment: 21 pages, 10 figures, accepted for publication in ApJ

Supernovae (SN) generate hot gas in the interstellar medium (ISM), help setting the ISM structure and support the driving of outflows. It is important to resolve the hot gas generation for galaxy formation simulations at solar mass and sub-parsec resolution which realise individual supernova (SN) explosions with ambient densities varying by several orders of magnitude in a realistic multi-phase ISM. We test resolution requirements by simulating SN blast waves at three metallicities ($Z = 0.01, 0.1$ and $1 Z_{\odot}$), six densities and their respective equilibrium chemical compositions ($n=0.001$ cm$^{-3}$ - $100$ cm$^{-3}$), and four mass resolutions ($0.1$ - $100$ M$_{\odot}$), in three dimensions. We include non-equilibrium cooling and chemistry, a homogenous interstellar radiation field, and shielding with a modern pressure-energy smoothed particle hydrodynamics (SPH) method including isotropic thermal conduction and a meshless-finite-mass (MFM) solver. We find stronger resolution requirements for chemistry and hot phase generation than for momentum generation. While at $10$ M$_{\odot}$ the radial momenta at the end of the Sedov phase start converging, the hot phase generation and chemistry require higher resolutions to represent the neutral to ionised hydrogen fraction at the end of the Sedov phase correctly. Thermal conduction typically reduces the hot phase by $0.2$ dex and has little impact on the chemical composition. In general, our $1$, and $0.1$ M$_{\odot}$ results agree well with previous numerical and analytic estimates. We conclude that for the thermal energy injection SN model presented here resolutions higher than $10$ M$_{\odot}$ are required to model the chemistry, momentum and hot phase generation in a multi-phase ISM., Comment: 28 pages, 24 figures (2 in the appendix), accepted for publication in MNRAS, comments welcome

We present a hydro-dynamical simulation at sub-parsec and few solar mass resolution of a merger between two gas-rich dwarf galaxies. Our simulation includes a detailed model for the multi-phase interstellar medium (ISM) and is able to follow the entire formation history of spatially resolved star clusters, including feedback from individual massive stars. Shortly after the merger we find a population of $\sim 900$ stellar clusters with masses above $10^{2.5}\; \rm{M_\odot}$ and a cluster mass function (CMF), which is well fitted with a power-law with a slope of $\alpha=-1.70\pm0.08$. We describe here in detail the formation of the three most massive clusters (M$_{*} \gtrsim 10^5$ M$_\odot$), which populate the high-mass end of the CMF. The simulated clusters form rapidly on a timescale of $6$-$8$ Myr in converging flows of dense gas. The embedded merger phase has extremely high star formation rate surface densities of $\Sigma_\mathrm{SFR}>10\; \mathrm{M}_\odot\; \mathrm{yr}^{-1}\; \mathrm{kpc}^{-2}$ and thermal gas pressures in excess of $P_{\rm th}\sim10^7 \; \mathrm{k}_{\rm B}\;(\rm K\;\mathrm{cm}^{-3})^{-1}$. The formation process is terminated by rapid gas expulsion driven by the first generation of supernovae, after which the cluster centers relax and both their structure and kinematics become indistinguishable from observed local globular clusters. The simulation presented here provides a general model for the formation of metal-poor globular clusters in chemically unevolved starbursting environments of low-mass dwarf galaxies, which are common at high redshifts., Comment: Accepted for publication in ApJL. 8 pages, 5 figures

We study the destruction of interstellar dust via sputtering in supernova (SN) shocks using three-dimensional hydrodynamical simulations. With a novel numerical framework, we follow both sputtering and dust dynamics governed by direct collisions, plasma drag and betatron acceleration. Grain-grain collisions are not included and the grain-size distribution is assumed to be fixed. The amount of dust destroyed per SN is quantified for a broad range of ambient densities and fitting formulae are provided. Integrated over the grain-size distribution, nonthermal (inertial) sputtering dominates over thermal sputtering for typical ambient densities. We present the first simulations that explicitly follow dust sputtering within a turbulent multiphase interstellar medium. We find that the dust destruction timescales $\tau$ are 0.35 Gyr for silicate dust and 0.44 Gyr for carbon dust in solar neighborhood conditions. The SN environment has an important impact on $\tau$. SNe that occur in preexisting bubbles destroy less dust as the destruction is limited by the amount of dust in the shocked gas. This makes $\tau$ about 2.5 times longer than the estimate based on results from a single SN explosion. We investigate the evolution of the dust-to-gas mass ratio (DGR), and find that a spatial inhomogeneity of $\sim$ 14\% develops for scales below 10 pc. It locally correlates positively with gas density but negatively with gas temperature even in the exterior of the bubbles due to incomplete gas mixing. This leads to a $\sim$ 30\% lower DGR in the volume filling warm gas compared to that in the dense clouds., Comment: 20 pages, 16 figures, accepted version

We investigate galactic winds driven by supernova (SN) explosions in an isolated dwarf galaxy using high-resolution (particle mass $m_{\rm gas} = 1{\rm M_\odot}$, number of neighbor $N_{\rm ngb} = 100$) smoothed-particle hydrodynamics simulations that include non-equilibrium cooling and chemistry, individual star formation, stellar feedback and metal enrichment. Clustered SNe lead to the formation of superbubbles which break out of the disk and vent out hot gas, launching the winds. We find much weaker winds than what cosmological simulations typically adopt at this mass scale. At the virial radius, the time-averaged loading factors of mass, momentum and energy are 3, 1 and 0.05, respectively, and the metal enrichment factor is 1.5. Winds that escape the halo consist of two populations that differ in their launching temperatures. Hot gas acquires enough kinetic energy to escape when launched while warm gas does not. However, warm gas can be further accelerated by the ram pressure of the subsequently launched hot gas and eventually escape. The strong interactions between different temperature phases highlight the caveat of extrapolating properties of warm gas to large distances based on its local conditions (e.g. the Bernoulli parameter). Our convergence study finds that wind properties converge when the cooling masses of individual SNe are resolved, which corresponds to $m_{\rm gas}=5 {\rm M_\odot}$ with an injection mass of $500 {\rm M_\odot}$. The winds weaken dramatically once the SNe become unresolved. We demonstrate that injecting the terminal momentum of SNe, a popular sub-grid model in the literature, fails to capture SN winds irrespective of the inclusion of residual thermal energy., Comment: Version accepted for publication, including discussion on effects of residual thermal energy and injection mass (appendix A & B)

The cloud-scale physics of star formation and feedback represent the main uncertainty in galaxy formation studies. Progress is hampered by the limited empirical constraints outside the restricted environment of the Local Group. In particular, the poorly-quantified time evolution of the molecular cloud lifecycle, star formation, and feedback obstructs robust predictions on the scales smaller than the disc scale height that are resolved in modern galaxy formation simulations. We present a new statistical method to derive the evolutionary timeline of molecular clouds and star-forming regions. By quantifying the excess or deficit of the gas-to-stellar flux ratio around peaks of gas or star formation tracer emission, we directly measure the relative rarity of these peaks, which allows us to derive their lifetimes. We present a step-by-step, quantitative description of the method and demonstrate its practical application. The method's accuracy is tested in nearly 300 experiments using simulated galaxy maps, showing that it is capable of constraining the molecular cloud lifetime and feedback time-scale to $<0.1$ dex precision. Access to the evolutionary timeline provides a variety of additional physical quantities, such as the cloud-scale star formation efficiency, the feedback outflow velocity, the mass loading factor, and the feedback energy or momentum coupling efficiencies to the ambient medium. We show that the results are robust for a wide variety of gas and star formation tracers, spatial resolutions, galaxy inclinations, and galaxy sizes. Finally, we demonstrate that our method can be applied out to high redshift ($z\lesssim4$) with a feasible time investment on current large-scale observatories. This is a major shift from previous studies that constrained the physics of star formation and feedback in the immediate vicinity of the Sun., Comment: 86 pages (including appendices), 30 figures, 16 tables; accepted by MNRAS (April 27, 2018). The method is summarised in Figure 2 and 3, with the main results shown in Figures 18, 26, and 27. For a cursory overview of the paper (of about 10 pages), please consult Sections 1, 2, 3.1, 4.4, 6, and 8. Figure 27 replaced with authors' original version

We present high-resolution hydrodynamical simulations of isolated dwarf galaxies including self-gravity, non-equilibrium cooling and chemistry, interstellar radiation fields (IRSF) and shielding, star formation, and stellar feedback. This includes spatially and temporally varying photoelectric (PE) heating, photoionization, resolved supernova (SN) blast waves and metal enrichment. A new flexible method to sample the stellar initial mass function allows us to follow the contribution to the ISRF, the metal output and the SN delay times of individual massive stars. We find that SNe play the dominant role in regulating the global star formation rate, shaping the multi-phase interstellar medium (ISM) and driving galactic outflows. Outflow rates (with mass-loading factors of a few) and hot gas fractions of the ISM increase with the number of SNe exploding in low-density environments where radiative energy losses are low. While PE heating alone can suppress star formation slightly more (a factor of a few) than SNe alone can do, it is unable to drive outflows and reproduce the multi-phase ISM that emerges naturally when SNe are included. These results are in conflict with recent results of Forbes et al. who concluded that PE heating is the dominant process suppressing star formation in dwarfs, about an order of magnitude more efficient than SNe. Potential origins for this discrepancy are discussed. In the absence of SNe and photoionization (mechanisms to disperse dense clouds), the impact of PE heating is highly overestimated owing to the (unrealistic) proximity of dense gas to the radiation sources. This leads to a substantial boost of the infrared continuum emission from the UV-irradiated dust and a far infrared line-to-continuum ratio too low compared to observations. Though sub-dominant in regulating star formation, the ISRF controls the abundance of molecular hydrogen via photodissociation., Comment: submitted, comments welcome

Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies and the overall inefficient conversion of gas into stars. Some models, however, require resolution dependent scalings. We present a sub-resolution model representing the three major phases of supernova blast wave evolution $-$free expansion, energy conserving Sedov-Taylor, and momentum conserving snowplow$-$ with energy scalings adopted from high-resolution interstellar-medium simulations in both uniform and multiphase media. We allow for the effects of significantly enhanced SN remnant propagation in a multiphase medium with the cooling radius scaling with the hot volume fraction, $f_{\mathrm{hot}}$, as $(1 - f_{\mathrm{hot}})^{-4/5}$. We also include winds from young massive stars and AGB stars, Str\"omgren sphere gas heating by massive stars, and a gas cooling limiting mechanism driven by radiative recombination of dense HII regions. We present initial tests for isolated Milky-Way like systems simulated with the GADGET based code SPHgal with improved SPH prescription. Compared to pure thermal SN input, the model significantly suppresses star formation at early epochs, with star formation extended both in time and space in better accord with observations. Compared to models with pure thermal SN feedback, the age at which half the stellar mass is assembled increases by a factor of 2.4, and the mass loading parameter and gas outflow rate from the galactic disk increase by a factor of 2. Simulation results are converged for a two order of magnitude variation in particle mass in the range (1.3$-$130)$\times 10^4$ solar masses., Comment: Accepted for publication in the ApJ

We ask how the inclusion of various physical heating processes due to the metal content of gas affect the evolution of massive galaxies and compute a suite of cosmological hydrodynamical simulations that follow these systems and their supermassive black holes. We use a smoothed particle hydrodynamics code with a pressure-entropy formulation and a more accurate treatment of the metal production, turbulent diffusion and cooling rate based on individual element abundances. The feedback models include (1) AGN feedback via high velocity BAL winds and Compton/photoionization heating, (2) explicit stellar feedback from multiple processes including powerful winds from supernova events, stellar winds from young massive stars and AGB stars as well as radiative heating within Stromgren spheres around massive stars, and (3) additional heating effects due to the presence of metals including grain photoelectric heating, metallicity dependent X-ray heating by nearby accreting black holes and from the cosmic X-ray background, which are the major improvement in our feedback model. With a suite of zoom-in simulations of 30 halos with $M_{vir} \sim 10^{12-13.4}$, we show that energy and momentum budget from all feedback effects generate realistic galaxy properties. We explore the detailed role of each feedback model with three additional sets of simulations with varying input physics. We show that the metal induced heating mechanisms reduce the fraction of accreted stellar material by mainly suppressing the growth of diffuse small stellar systems at high redshift but overall have a relatively minor effect on the final stellar and gas properties of massive galaxies. The inclusion of AGN feedback significantly improves the ability of our cosmological simulations to yield realistic gas and stellar properties of massive galaxies with reasonable fraction of the final stellar mass accreted from other galaxies., Comment: 16 pages, 11 figures, accepted for publication in ApJ

We study the connection of star formation to atomic (HI) and molecular hydrogen (H$_2$) in isolated, low metallicity dwarf galaxies with high-resolution ($m_{\rm gas}$ = 4 M$_\odot$, $N_{\rm ngb}$ = 100) SPH simulations. The model includes self-gravity, non-equilibrium cooling, shielding from an interstellar radiation field, the chemistry of H$_2$ formation, H$_2$-independent star formation, supernova feedback and metal enrichment. We find that the H$_2$ mass fraction is sensitive to the adopted dust-to-gas ratio and the strength of the interstellar radiation field, while the star formation rate is not. Star formation is regulated by stellar feedback, keeping the gas out of thermal equilibrium for densities $n <$ 1 cm$^{-3}$. Because of the long chemical timescales, the H$_2$ mass remains out of chemical equilibrium throughout the simulation. Star formation is well-correlated with cold ( T $\leqslant$ 100 K ) gas, but this dense and cold gas - the reservoir for star formation - is dominated by HI, not H$_2$. In addition, a significant fraction of H$_2$ resides in a diffuse, warm phase, which is not star-forming. The ISM is dominated by warm gas (100 K $<$ T $\leqslant 3\times 10^4$ K) both in mass and in volume. The scale height of the gaseous disc increases with radius while the cold gas is always confined to a thin layer in the mid-plane. The cold gas fraction is regulated by feedback at small radii and by the assumed radiation field at large radii. The decreasing cold gas fractions result in a rapid increase in depletion time (up to 100 Gyrs) for total gas surface densities $\Sigma_{\rm HI+H_2} \lesssim$ 10 M$_\odot$pc$^{-2}$, in agreement with observations of dwarf galaxies in the Kennicutt-Schmidt plane., Comment: Accepted for publication in MNRAS. Changes (including a pamameter study in Appendix C) highlighted

We present the smoothed-particle hydrodynamics implementation SPHGal, which combines some recently proposed improvements in GADGET. This includes a pressure-entropy formulation with a Wendland kernel, a higher order estimate of velocity gradients, a modified artificial viscosity switch with a modified strong limiter, and artificial conduction of thermal energy. With a series of idealized hydrodynamic tests we show that the pressure-entropy formulation is ideal for resolving fluid mixing at contact discontinuities but performs conspicuously worse at strong shocks due to the large entropy discontinuities. Including artificial conduction at shocks greatly improves the results. In simulations of Milky Way like disk galaxies a feedback-induced instability develops if too much artificial viscosity is introduced. Our modified artificial viscosity scheme prevents this instability and shows efficient shock capturing capability. We also investigate the star formation rate and the galactic outflow. The star formation rates vary slightly for different SPH schemes while the mass loading is sensitive to the SPH scheme and significantly reduced in our favored implementation. We compare the accretion behavior of the hot halo gas. The formation of cold blobs, an artifact of simple SPH implementations, can be eliminated efficiently with proper fluid mixing, either by conduction and/or by using a pressure-entropy formulation., Comment: Replaced with the version accepted by MNRAS

The radio approach for detecting the ultra-high energy cosmic neutrinos has become a mature field. The Cherenkov pulse in radio detection originates from the charge excess of particle showers due to Askaryan effect. The conventional way of calculating the Cherenkov pulse by making far- field approximation fails when the size of elongated showers become comparable with detection distance. We investigate the Cherenkov pulse in near-field by a numerical code based on the finite- difference time-domain (FDTD) method. Our study shows that the near-field radiation exhibits very different behaviors from the far-field one and therefore can be easily recognized. For ground array neutrino detectors, the near-field radiation would provide a unique signature for ultra high energy electromagnetic showers induced by the electron neutrino charge-current interaction. This can be useful in neutrino flavor identification., Comment: 14 pages, 13 figures

The radio approach for detecting the ultra-high energy cosmic neutrinos has become a mature field. The Cherenkov signals in radio detection are originated from the charge excess of particle showers due to Askaryan effect. The conventional way of calculating the Cherenkov pulses by making Fraunhofer approximation fails when the sizes of the elongated showers become comparable with the detection distances. We present a calculation method of Cherenkov pulses based on the finite-difference time-domain (FDTD) method, and attain a satisfying effeciency via the GPU- acceleration. Our method provides a straightforward way of the near field calculation, which would be important for ultra high energy particle showers, especailly the electromagnetic showers induced by the high energy leptons produced in the neutrino charge current interactions., Comment: 8 pages, 9 figures, presented at the CosPA Symposium, Melbourne, Australia, November 2009