ABSTRACT
We introduce the thesan project, a suite of large volume ($L_\mathrm{box} = 95.5 \, \mathrm{cMpc}$) radiation-magnetohydrodynamic simulations that simultaneously model the large-scale ...statistical properties of the intergalactic medium during reionization and the resolved characteristics of the galaxies responsible for it. The flagship simulation has dark matter and baryonic mass resolutions of $3.1 \times 10^6\, {\rm M_\odot }$ and $5.8 \times 10^5\, {\rm M_\odot }$, respectively. The gravitational forces are softened on scales of 2.2 ckpc with the smallest cell sizes reaching 10 pc at z = 5.5, enabling predictions down to the atomic cooling limit. The simulations use an efficient radiation hydrodynamics solver (arepo-rt) that precisely captures the interaction between ionizing photons and gas, coupled to well-tested galaxy formation (IllustrisTNG) and dust models to accurately predict the properties of galaxies. Through a complementary set of medium resolution simulations we investigate the changes to reionization introduced by different assumptions for ionizing escape fractions, varying dark matter models, and numerical convergence. The fiducial simulation and model variations are calibrated to produce realistic reionization histories that match the observed evolution of the global neutral hydrogen fraction and electron scattering optical depth to reionization. They also match a wealth of high-redshift observationally inferred data, including the stellar-to-halo-mass relation, galaxy stellar mass function, star formation rate density, and the mass–metallicity relation, despite the galaxy formation model being mainly calibrated at z = 0. We demonstrate that different reionization models give rise to varied bubble size distributions that imprint unique signatures on the 21 cm emission, especially on the slope of the power spectrum at large spatial scales, enabling current and upcoming 21 cm experiments to accurately characterize the sources that dominate the ionizing photon budget.
ABSTRACT
The diffuse medium in and around galaxies can exist in a multiphase state: small, cold gas clouds contributing significantly to the total mass embedded in pressure equilibrium with a hotter, ...more diffuse volume-filling component. Modelling this multiphase state in cosmological simulations poses a significant challenge due to the requirements to spatially resolve the clouds and consequently the interactions between the phases. In this paper, we present a novel method to model this gas state in cosmological hydrodynamical simulations. We solve the compressible two-fluid hydrodynamic equations using a moving-mesh finite-volume method and define mass, momentum, and energy exchange terms between the phases as operator-split source terms. Using a stratified flow model, our implementation is able to maintain volume fraction discontinuities in pressure equilibrium to machine precision, allowing for the treatment of both resolved and unresolved multiphase fluids. The solver remains second order accurate on smooth hydrodynamics problems. We use the source and sink terms of an existing two-phase model for the interstellar medium to demonstrate the value of this type of approach in simulations of galaxy formation, compare it to its effective equation of state implementation, and discuss its advantages in future large-scale simulations of galaxy formation.
We present a model for star formation and supernova feedback that describes the multiphase structure of star-forming gas on scales that are typically not resolved in cosmological simulations. Our ...approach includes radiative heating and cooling, the growth of cold clouds embedded in an ambient hot medium, star formation in these clouds, feedback from supernovae in the form of thermal heating and cloud evaporation, galactic winds and outflows, and metal enrichment. Implemented using smoothed particle hydrodynamics, our scheme is a significantly modified and extended version of the grid-based method of Yepes et al., and enables us to achieve a high dynamic range in simulations of structure formation. We discuss properties of the feedback model in detail and show that it predicts a self-regulated, quiescent mode of star formation, which, in particular, stabilizes the star-forming gaseous layers of disc galaxies. The parametrization of this mode can be reduced to a single free quantity that determines the overall time-scale for star formation. We fix this parameter numerically to match the observed rates of star formation in local disc galaxies. When normalized in this manner, cosmological simulations employing our model nevertheless overproduce the observed cosmic abundance of stellar material. We are thus motivated to extend our feedback model to include galactic winds associated with star formation. Using small-scale simulations of individual star-forming disc galaxies, we show that these winds produce either galactic fountains or outflows, depending on the depth of the gravitational potential. In low-mass haloes, winds can greatly suppress the overall efficiency of star formation. When incorporated into cosmological simulations, our combined model for star formation and winds predicts a cosmic star formation density that is consistent with observations, provided that the winds are sufficiently energetic. Moreover, outflows from galaxies in these simulations drive chemical enrichment of the intergalactic medium — in principle, accounting for the presence of metals in the Lyman α forest.
ABSTRACT We use a set of high-resolution N-body simulations of the Galactic disk to study its interactions with the population of cosmologically predicted satellites. One simulation illustrates that ...multiple passages of massive satellites with different velocities through the disk generate a wobble, which has the appearance of rings in face-on projections of the stellar disk. They also produce flares in the outer disk parts and gradually heat the disk through bending waves. A different numerical experiment shows that an individual satellite as massive as the Sagittarius dwarf galaxy passing through the disk will drive coupled horizontal and vertical oscillations of stars in underdense regions with small associated heating. This experiment shows that vertical excursions of stars in these low-density regions can exceed 1 kpc in the Solar neighborhood, resembling the recently locally detected coherent vertical oscillations. They can also induce non-zero vertical streaming motions as large as 10-20 km s−1, which is consistent with recent observations in the Galactic disk. This phenomenon appears as a local ring with modest associated disk heating.
Abstract
We discuss differences in simulation results that arise between the use of either the thermal energy or the entropy as an independent variable in smoothed particle hydrodynamics (SPH). In ...this context, we derive a new version of SPH that, when appropriate, manifestly conserves both energy and entropy if smoothing lengths are allowed to adapt freely to the local mass resolution. To test various formulations of SPH, we consider point-like energy injection, as in certain models of supernova feedback, and find that powerful explosions are well represented by SPH even when the energy is deposited into a single particle, provided that the entropy equation is integrated. If the thermal energy is instead used as an independent variable, unphysical solutions can be obtained for this problem.
We also examine the radiative cooling of gas spheres that collapse and virialize in isolation, and of haloes that form in cosmological simulations of structure formation. When applied to these problems, the thermal energy version of SPH leads to substantial overcooling in haloes that are resolved with up to a few thousand particles, while the entropy formulation is biased only moderately low for these haloes under the same circumstances. For objects resolved with much larger particle numbers, the two approaches yield consistent results. We trace the origin of the differences to systematic resolution effects in the outer parts of cooling flows. When the thermal energy equation is integrated and the resolution is low, the compressional heating of the gas in the inflow region is underestimated, violating entropy conservation and improperly accelerating cooling.
The cumulative effect of this overcooling can be significant. In cosmological simulations of moderate size, we find that the fraction of baryons which cool and condense can be reduced by up to a factor ∼2 if the entropy equation is employed rather than the thermal energy equation, partly explaining discrepancies with semi-analytic treatments of galaxy formation. We also demonstrate that the entropy method leads to a greatly reduced scatter in the density-temperature relation of the low-density Lyα forest relative to the thermal energy approach, in accord with theoretical expectations.
Previous simulations of the growth of cosmic structures have broadly reproduced the 'cosmic web' of galaxies that we see in the Universe, but failed to create a mixed population of elliptical and ...spiral galaxies, because of numerical inaccuracies and incomplete physical models. Moreover, they were unable to track the small-scale evolution of gas and stars to the present epoch within a representative portion of the Universe. Here we report a simulation that starts 12 million years after the Big Bang, and traces 13 billion years of cosmic evolution with 12 billion resolution elements in a cube of 106.5 megaparsecs a side. It yields a reasonable population of ellipticals and spirals, reproduces the observed distribution of galaxies in clusters and characteristics of hydrogen on large scales, and at the same time matches the 'metal' and hydrogen content of galaxies on small scales.
Celotno besedilo
Dostopno za:
DOBA, IJS, IZUM, KILJ, KISLJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
ABSTRACT
The high-redshift intergalactic medium (IGM) and the primeval galaxy population are rapidly becoming the new frontier of extragalactic astronomy. We investigate the IGM properties and their ...connection to galaxies at z ≥ 5.5 under different assumptions for the ionizing photon escape and the nature of dark matter, employing our novel thesan radiation-hydrodynamical simulation suite, designed to provide a comprehensive picture of the emergence of galaxies in a full reionization context. Our simulations have realistic ‘late’ reionization histories, match available constraints on global IGM properties, and reproduce the recently observed rapid evolution of the mean free path of ionizing photons. We additionally examine high-z Lyman-α transmission. The optical depth evolution is consistent with data, and its distribution suggests an even-later reionization than simulated, although with a strong sensitivity to the source model. We show that the effects of these two unknowns can be disentangled by characterizing the spectral shape and separation of Lyman-α transmission regions, opening up the possibility to observationally constrain both. For the first time in simulations, thesan reproduces the modulation of the Lyman-α flux as a function of galaxy distance, demonstrating the power of coupling a realistic galaxy formation model with proper radiation hydrodynamics. We find this feature to be extremely sensitive on the timing of reionization, while being relatively insensitive to the source model. Overall, thesan produces a realistic IGM and galaxy population, providing a robust framework for future analysis of the high-z Universe.
Context.
Early-type galaxies (ETGs) are found to follow a wide variety of merger and accretion histories in cosmological simulations.
Aims.
We characterize the photometric and kinematic properties of ...simulated ETG stellar halos, and compare them to the observations.
Methods.
We selected a sample of 1114 ETGs in the TNG100 simulation and 80 in the higher-resolution TNG50. These ETGs span a stellar mass range of 10
10.3
− 10
12
M
⊙
and they were selected within the range of
g
−
r
colour and
λ
-ellipticity diagram populated by observed ETGs. We determined photometric parameters, intrinsic shapes, and kinematic observables in their extended stellar halos. We compared the results with central IFU kinematics and ePN.S planetary nebula velocity fields at large radii, studying the variation in kinematics from center to halo, and connecting it to a change in the intrinsic shape of the galaxies.
Results.
We find that the simulated galaxy sample reproduces the diversity of kinematic properties observed in ETG halos. Simulated fast rotators (FRs) divide almost evenly in one third having flat
λ
profiles and high halo rotational support, a third with gently decreasing profiles, and another third with low halo rotation. However, the peak of rotation occurs at larger
R
than in observed ETG samples. Slow rotators (SRs) tend to have increased rotation in the outskirts, with half of them exceeding
λ
= 0.2. For
M
*
> 10
11.5
M
⊙
halo rotation is unimportant. A similar variety of properties is found for the stellar halo intrinsic shapes. Rotational support and shape are deeply related: the kinematic transition to lower rotational support is accompanied by a change towards rounder intrinsic shape. Triaxiality in the halos of FRs increases outwards and with stellar mass. Simulated SRs have relatively constant triaxiality profiles.
Conclusions.
Simulated stellar halos show a large variety of structural properties, with quantitative but no clear qualitative differences between FRs and SRs. At the same stellar mass, stellar halo properties show a more gradual transition and significant overlap between the two families, despite the clear bimodality in the central regions. This is in agreement with observations of extended photometry and kinematics.
ABSTRACT
The visibility of high-redshift Lyman-alpha emitting galaxies (LAEs) provides important constraints on galaxy formation processes and the Epoch of Reionization (EoR). However, predicting ...realistic and representative statistics for comparison with observations represents a significant challenge in the context of large-volume cosmological simulations. The thesan project offers a unique framework for addressing such limitations by combining state-of-the-art galaxy formation (IllustrisTNG) and dust models with the arepo-rt radiation-magnetohydrodynamics solver. In this initial study, we present Lyman-alpha centric analysis for the flagship simulation that resolves atomic cooling haloes throughout a $(95.5\, \text{cMpc})^3$ region of the Universe. To avoid numerical artefacts, we devise a novel method for accurate frequency-dependent line radiative transfer in the presence of continuous Hubble flow, transferable to broader astrophysical applications as well. Our scalable approach highlights the utility of LAEs and red damping-wing transmission as probes of reionization, which reveal nontrivial trends across different galaxies, sightlines, and frequency bands that can be modelled in the framework of covering fractions. In fact, after accounting for environmental factors influencing large-scale ionized bubble formation such as redshift and UV magnitude, the variation across galaxies and sightlines mainly depends on random processes including peculiar velocities and self-shielded systems that strongly impact unfortunate rays more than others. Throughout the EoR local and cosmological optical depths are often greater than or less than unity such that the exp (− τ) behaviour leads to anisotropic and bimodal transmissivity. Future surveys will benefit by targeting both rare bright objects and Goldilocks zone LAEs to infer the presence of these (un)predictable (dis)advantages.
Dust entrainment in galactic winds Kannan, R; Vogelsberger, M; Marinacci, F ...
Monthly notices of the Royal Astronomical Society,
05/2021, Letnik:
503, Številka:
1
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
Winds driven by stellar feedback are an essential part of the galactic ecosystem and are the main mechanism through which low-mass galaxies regulate their star formation. These winds are ...generally observed to be multiphase with detections of entrained neutral and molecular gas. They are also thought to enrich the circumgalactic medium around galaxies with metals and dust. This ejected dust encodes information about the integrated star formation and outflow history of the galaxy. Therefore it is important to understand how much dust is entrained and driven out of the disc by galactic winds. Here, we demonstrate that stellar feedback is efficient in driving dust-enriched winds and eject enough material to account for the amount of extraplanar dust observed in nearby galaxies. The amount of ejected dust depends on the sites from where they are launched, with dustier galaxies launching more dust-enriched outflows. Moreover, the outflowing cold and dense gas is significantly more dust enriched than the volume filling hot and tenuous material. These results provide an important new insight into the dynamics, structure, and composition of galactic winds and their role in determining the dust content of the extragalactic gas in galaxies.