Abstract
We present full-volume cosmological simulations, using the moving-mesh code arepo to study the coevolution of dust and galaxies. We extend the dust model in arepo to include thermal ...sputtering of grains and investigate the evolution of the dust mass function, the cosmic distribution of dust beyond the interstellar medium and the dependence of dust-to-stellar mass ratio on galactic properties. The simulated dust mass function is well described by a Schechter fit and lies closest to observations at z = 0. The radial scaling of projected dust surface density out to distances of 10 Mpc around galaxies with magnitudes 17 < i < 21 is similar to that seen in Sloan Digital Sky Survey data, albeit with a lower normalization. At z = 0, the predicted dust density of Ωdust ≈ 1.3 × 10−6 lies in the range of Ωdust values seen in low-redshift observations. We find that the dust-to-stellar mass ratio anticorrelates with stellar mass for galaxies living along the star formation main sequence. Moreover, we estimate the 850 μm number density functions for simulated galaxies and analyse the relation between dust-to-stellar flux and mass ratios at z = 0. At high redshift, our model fails to produce enough dust-rich galaxies, and this tension is not alleviated by adopting a top-heavy initial mass function. We do not capture a decline in Ωdust from z = 2 to 0, which suggests that dust production mechanisms more strongly dependent on star formation may help to produce the observed number of dusty galaxies near the peak of cosmic star formation.
We present an analytic model for how momentum deposition from stellar feedback simultaneously regulates star formation and drives outflows in a turbulent interstellar medium (ISM). Because the ISM is ...turbulent, a given patch of ISM exhibits sub-patches with a range of surface densities. The high-density patches are 'pushed' by feedback, thereby driving turbulence and self-regulating local star formation. Sufficiently low-density patches, however, are accelerated to above the escape velocity before the region can self-adjust and are thus vented as outflows. When the gas fraction is ... 0.3, the ratio of the turbulent velocity dispersion to the circular velocity is sufficiently high that at any given time, of the order of half of the ISM has surface density less than the critical value and thus can be blown out on a dynamical time. The resulting outflows have a mass-loading factor (...M sub( out)/M*) that is inversely proportional to the gas fraction times the circular velocity. At low gas fractions, the star formation rate needed for local self-regulation, and corresponding turbulent Mach number, declines rapidly; the ISM is 'smoother', and it is actually more difficult to drive winds with large mass-loading factors. Crucially, our model predicts that stellar-feedback-driven outflows should be suppressed at z ... 1 in M* ... 10 super( 10) M... galaxies. This mechanism allows massive galaxies to exhibit violent outflows at high redshifts and then 'shut down' those outflows at late times, thereby enabling the formation of a smooth, extended thin stellar disc. We provide simple fitting functions for ... that should be useful for sub-resolution and semi-analytic models. (ProQuest: ... denotes formulae/symbols omitted.)
Understanding the physical processes that drive star formation is a key challenge for galaxy formation models. In this paper, we study the tight correlation between the star formation rate (SFR) and ...stellar mass of galaxies at a given redshift, how halo growth influences star formation, and star formation histories of individual galaxies. We study these topics using Illustris, a state-of-the-art cosmological hydrodynamical simulation of galaxy formation. Illustris reproduces the observed relation (the star formation main sequence, SFMS) between SFR and stellar mass at redshifts z = 0 and 4, but at intermediate redshifts of z ... 1-2, the simulated SFMS has a significantly lower normalization than reported by observations. The scatter in the relation is consistent with the observed scatter. However, the fraction of outliers above the SFR-stellar mass relation in Illustris is less than that observed. Galaxies with halo masses of ~10... M... dominate the SFR density of the Universe, in agreement with the results of abundance matching. Furthermore, more-massive galaxies tend to form the bulk of their stars at high redshift, which indicates that 'downsizing' occurs in Illustris. We also studied the star formation histories of individual galaxies, including the use of a principal component analysis decomposition. We find that for fixed stellar mass, galaxies that form earlier have more-massive black holes at z = 0, indicating that star formation and black hole growth are tightly linked processes in Illustris. While many of the properties of normal star-forming galaxies are well reproduced in the Illustris simulation, forming a realistic population of starbursts will likely require higher resolution and probably a more sophisticated treatment of star formation and feedback from stars and black holes. (ProQuest: ... denotes formulae/symbols omitted.)
Abstract
Galaxy formation models are now able to reproduce observed relations such as the relation between galaxies’ star formation rates (SFRs) and stellar masses (M
*) and the ...stellar-mass–halo-mass relation. We demonstrate that comparisons of the short-time-scale variability in galaxy SFRs with observational data provide an additional useful constraint on the physics of galaxy formation feedback. We apply SFR indicators with different sensitivity time-scales to galaxies from the Feedback in Realistic Environments (FIRE) simulations. We find that the SFR–M
* relation has a significantly greater scatter when the Hα-derived SFR is considered compared with when the far-ultraviolet (FUV)-based SFR is used. This difference is a direct consequence of bursty star formation because the FIRE galaxies exhibit order-of-magnitude SFR variations over time-scales of a few Myr. We show that the difference in the scatter between the simulated Hα- and FUV-derived SFR–M
* relations at z = 2 is consistent with observational constraints. We also find that the Hα/FUV ratios predicted by the simulations at z = 0 are similar to those observed for local galaxies except for a population of low-mass (M
* ≲ 109.5 M⊙) simulated galaxies with lower Hα/FUV ratios than observed. We suggest that future cosmological simulations should compare the Hα/FUV ratios of their galaxies with observations to constrain the feedback models employed.
Abstract
We present a suite of 34 high-resolution cosmological zoom-in simulations consisting of thousands of haloes up to $M_{\rm halo}\sim 10^{12}\, \mathrm{M}_{\odot }$ ($M_{\ast }\sim 10^{10.5}\, ...\mathrm{M}_{\odot }$) at z ≥ 5 from the Feedback in Realistic Environments project. We post-process our simulations with a three-dimensional Monte Carlo dust radiative transfer code to study dust attenuation, dust emission, and dust temperature within these simulated z ≥ 5 galaxies. Our sample forms a tight correlation between infrared excess (IRX ≡ FIR/FUV) and ultraviolet (UV)-continuum slope (βUV), despite the patchy, clumpy dust geometry shown in our simulations. We find that the IRX–βUV relation is mainly determined by the shape of the attenuation law and is independent of its normalization (set by the dust-to-gas ratio). The bolometric IR luminosity (LIR) correlates with the intrinsic UV luminosity and the star formation rate (SFR) averaged over the past 10 Myr. We predict that at a given LIR, the peak wavelength of the dust spectral energy distributions for z ≥ 5 galaxies is smaller by a factor of 2 (due to higher dust temperatures on average) than at z = 0. The higher dust temperatures are driven by higher specific SFRs and SFR surface densities with increasing redshift. We derive the galaxy UV luminosity functions (UVLFs) at z = 5–10 from our simulations and confirm that a heavy attenuation is required to reproduce the observed bright-end UVLFs. We also predict the IR luminosity functions (IRLFs) and UV luminosity densities at z = 5–10. We discuss the implications of our results on current and future observations probing dust attenuation and emission in z ≥ 5 galaxies.
Large-scale outflows in star-forming galaxies are observed to be ubiquitous and are a key aspect of theoretical modeling of galactic evolution, the focus of the Simulating Multiscale Astrophysics to ...Understand Galaxies (SMAUG) project. Gas blown out from galactic disks, similar to gas within galaxies, consists of multiple phases with large contrasts of density, temperature, and other properties. To study multiphase outflows as emergent phenomena, we run a suite of rougly parsec-resolution local galactic disk simulations using the TIGRESS framework. Explicit modeling of the interstellar medium (ISM), including star formation and self-consistent radiative heating plus supernova feedback, regulates ISM properties and drives the outflow. We investigate the scaling of outflow mass, momentum, energy, and metal loading factors with galactic disk properties, including star formation rate (SFR) surface density ( SFR ∼ 10−4 − 1 M kpc−2 yr−1), gas surface density ( ), and total midplane pressure (or weight; ). The main components of outflowing gas are mass-delivering cool gas (T ∼ 104 K) and energy/metal-delivering hot gas (T 106 K). Cool mass outflow rates measured at outflow launch points (one or two scale heights ) are 1-100 times the SFR (decreasing with SFR), although in massive galaxies most mass falls back owing to insufficient outflow velocity. The hot galactic outflow carries mass comparable to 10% of the SFR, together with 10%-20% of the energy and 30%-60% of the metal mass injected by SN feedback. Importantly, our analysis demonstrates that in any physically motivated cosmological wind model it is crucial to include at least two distinct thermal wind components.
ABSTRACT
Machine learning is becoming a popular tool to quantify galaxy morphologies and identify mergers. However, this technique relies on using an appropriate set of training data to be ...successful. By combining hydrodynamical simulations, synthetic observations, and convolutional neural networks (CNNs), we quantitatively assess how realistic simulated galaxy images must be in order to reliably classify mergers. Specifically, we compare the performance of CNNs trained with two types of galaxy images, stellar maps and dust-inclusive radiatively transferred images, each with three levels of observational realism: (1) no observational effects (idealized images), (2) realistic sky and point spread function (semirealistic images), and (3) insertion into a real sky image (fully realistic images). We find that networks trained on either idealized or semireal images have poor performance when applied to survey-realistic images. In contrast, networks trained on fully realistic images achieve 87.1 per cent classification performance. Importantly, the level of realism in the training images is much more important than whether the images included radiative transfer, or simply used the stellar maps ($87.1{{\ \rm per\ cent}}$ compared to $79.6{{\ \rm per\ cent}}$ accuracy, respectively). Therefore, one can avoid the large computational and storage cost of running radiative transfer with a relatively modest compromise in classification performance. Making photometry-based networks insensitive to colour incurs a very mild penalty to performance with survey-realistic data ($86.0{{\ \rm per\ cent}}$ with r-only compared to $87.1{{\ \rm per\ cent}}$ with gri). This result demonstrates that while colour can be exploited by colour-sensitive networks, it is not necessary to achieve high accuracy and so can be avoided if desired. We provide the public release of our statistical observational realism suite, RealSim, as a companion to this paper.
On the dust temperatures of high-redshift galaxies Liang, Lichen; Feldmann, Robert; Kereš, Dušan ...
Monthly notices of the Royal Astronomical Society,
10/2019, Letnik:
489, Številka:
1
Journal Article
Recenzirano
Odprti dostop
Abstract
Dust temperature is an important property of the interstellar medium (ISM) of galaxies. It is required when converting (sub)millimetre broad-band flux to total infrared luminosity (LIR), and ...hence star formation rate, in high-redshift galaxies. However, different definitions of dust temperatures have been used in the literature, leading to different physical interpretations of how ISM conditions change with, e.g. redshift and star formation rate. In this paper, we analyse the dust temperatures of massive ($M_{\rm star} \gt 10^{10}\, \mathrm{M}_{\odot }$) $z$ = 2–6 galaxies with the help of high-resolution cosmological simulations from the Feedback in Realistic Environments (fire) project. At $z$ ∼ 2, our simulations successfully predict dust temperatures in good agreement with observations. We find that dust temperatures based on the peak emission wavelength increase with redshift, in line with the higher star formation activity at higher redshift, and are strongly correlated with the specific star formation rate. In contrast, the mass-weighted dust temperature, which is required to accurately estimate the total dust mass, does not strongly evolve with redshift over $z$ = 2–6 at fixed IR luminosity but is tightly correlated with LIR at fixed $z$. We also analyse an ‘equivalent’ dust temperature for converting (sub)millimetre flux density to total IR luminosity, and provide a fitting formula as a function of redshift and dust-to-metal ratio. We find that galaxies of higher equivalent (or higher peak) dust temperature (‘warmer dust’) do not necessarily have higher mass-weighted temperatures. A ‘two-phase’ picture for interstellar dust can explain the different scaling relations of the various dust temperatures.
PDF
