ABSTRACT The physics of cosmic rays (CRs) is a promising candidate for explaining the driving of galactic winds and outflows. Recent galaxy formation simulations have demonstrated the need for active ...CR transport either in the form of diffusion or streaming to successfully launch winds in galaxies. However, due to computational limitations, most previous simulations have modeled CR transport isotropically. Here, we discuss high-resolution simulations of isolated disk galaxies in a 1011 M halo with the moving-mesh code Arepo that include injection of CRs from supernovae, advective transport, CR cooling, and CR transport through isotropic or anisotropic diffusion. We show that either mode of diffusion leads to the formation of strong bipolar outflows. However, they develop significantly later in the simulation with anisotropic diffusion compared to the simulation with isotropic diffusion. Moreover, we find that isotropic diffusion allows most of the CRs to quickly diffuse out of the disk, while in the simulation with anisotropic diffusion, most CRs remain in the disk once the magnetic field becomes dominated by its azimuthal component, which occurs after ∼300 Myr. This has important consequences for the gas dynamics in the disk. In particular, we show that isotropic diffusion strongly suppresses the amplification of the magnetic field in the disk compared to anisotropic or no diffusion models. We therefore conclude that reliable simulations which include CR transport inevitably need to account for anisotropic diffusion.
Cosmological models in which dark matter consists of cold elementary particles predict that the dark halo population should extend to masses many orders of magnitude below those at which galaxies can ...form
. Here we report a cosmological simulation of the formation of present-day haloes over the full range of observed halo masses (20 orders of magnitude) when dark matter is assumed to be in the form of weakly interacting massive particles of mass approximately 100 gigaelectronvolts. The simulation has a full dynamic range of 30 orders of magnitude in mass and resolves the internal structure of hundreds of Earth-mass haloes in as much detail as it does for hundreds of rich galaxy clusters. We find that halo density profiles are universal over the entire mass range and are well described by simple two-parameter fitting formulae
. Halo mass and concentration are tightly related in a way that depends on cosmology and on the nature of the dark matter. For a fixed mass, the concentration is independent of the local environment for haloes less massive than those of typical galaxies. Haloes over the mass range of 10
to 10
solar masses contribute about equally (per logarithmic interval) to the luminosity produced by dark matter annihilation, which we find to be smaller than all previous estimates by factors ranging up to one thousand
.
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.
The progenitors of Type Ia supernovae (SNe Ia) are still unknown, despite significant progress during the past several years in theory and observations. Violent mergers of two carbon-oxygen (CO) ...white dwarfs (WDs) are a candidate scenario suggested to be responsible for at least a significant fraction of normal SNe Ia. Here, we simulate the merger of two CO WDs using a moving-mesh code that allows for the inclusion of thin helium (He) shells (0.01 M sub(middot in circle)) on top of the WDs at an unprecedented numerical resolution. The accretion of He onto the primary WD leads to the formation of a detonation in its He shell. This detonation propagates around the CO WD and sends a converging shock wave into its core, known to robustly trigger a second detonation, as in the well-known double-detonation scenario for He-accreting CO WDs. However, in contrast to that scenario where a massive He shell is required to form a detonation through thermal instability, here the He detonation is ignited dynamically. Accordingly the required He-shell mass is significantly smaller, and hence its burning products are unlikely to affect the optical display of the explosion. We show that this scenario, which works for CO primary WDs with CO- as well as He-WD companions, has the potential to explain the different brightness distributions, delay times, and relative rates of normal and fast declining SNe Ia. Finally, we discuss extensions to our unified merger model needed to obtain a comprehensive picture of the full observed diversity of SNe Ia.
The abundance and structure of dark matter subhaloes have been analysed extensively in recent studies of dark-matter-only simulations, but comparatively little is known about the impact of baryonic ...physics on halo substructures. We here extend the subfind algorithm for substructure identification such that it can be reliably applied to dissipative hydrodynamical simulations that include star formation. This allows, in particular, the identification of galaxies as substructures in simulations of clusters of galaxies and a determination of their content of gravitationally bound stars, dark matter and hot and cold gas. Using a large set of cosmological cluster simulations, we present a detailed analysis of halo substructures in hydrodynamical simulations of galaxy clusters, focusing in particular on the influence both of radiative and non-radiative gas physics and of non-standard physics such as thermal conduction and feedback by galactic outflows. We also examine the impact of numerical nuisance parameters such as artificial viscosity parameterizations. We find that diffuse hot gas is efficiently stripped from subhaloes when they enter the highly pressurized cluster atmosphere. This has the effect of decreasing the subhalo mass function relative to a corresponding dark-matter-only simulation. These effects are mitigated in radiative runs, where baryons condense in the central subhalo regions and form compact stellar cores. However, in all cases, only a very small fraction, of the order of one per cent, of subhaloes within the cluster virial radii preserve a gravitationally bound hot gaseous atmosphere. The fraction of mass contributed by gas in subhaloes is found to increase with the cluster-centric distance. Interestingly, this trend extends well beyond the virial radii, thus showing that galaxies feel the environment of the pressurized cluster gas over fairly large distances. The compact stellar cores (i.e. galaxies) are generally more resistant against tidal disruption than pure dark matter subhaloes. Still, the fraction of star-dominated substructures within our simulated clusters is only ∼10 per cent. We expect that the finite resolution in our simulations makes the galaxies overly susceptible to tidal disruption, hence the above fraction of star-dominated galaxies should represent a lower limit for the actual fraction of galaxies surviving the disruption of their host dark matter subhalo.
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 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.
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.
Simulating cosmic ray physics on a moving mesh Pfrommer, C; Pakmor, R; Schaal, K ...
Monthly notices of the Royal Astronomical Society,
03/2017, Letnik:
465, Številka:
4
Journal Article
Recenzirano
Odprti dostop
We discuss new methods to integrate the cosmic ray (CR) evolution equations coupled to magnetohydrodynamics on an unstructured moving mesh, as realized in the massively parallel arepo code for ...cosmological simulations. We account for diffusive shock acceleration of CRs at resolved shocks and at supernova remnants in the interstellar medium (ISM) and follow the advective CR transport within the magnetized plasma, as well as anisotropic diffusive transport of CRs along the local magnetic field. CR losses are included in terms of Coulomb and hadronic interactions with the thermal plasma. We demonstrate the accuracy of our formalism for CR acceleration at shocks through simulations of plane-parallel shock tubes that are compared to newly derived exact solutions of the Riemann shock-tube problem with CR acceleration. We find that the increased compressibility of the post-shock plasma due to the produced CRs decreases the shock speed. However, CR acceleration at spherically expanding blast waves does not significantly break the self-similarity of the Sedov-Taylor solution; the resulting modifications can be approximated by a suitably adjusted, but constant adiabatic index. In first applications of the new CR formalism to simulations of isolated galaxies and cosmic structure formation, we find that CRs add an important pressure component to the ISM that increases the vertical scaleheight of disc galaxies and thus reduces the star formation rate. Strong external structure formation shocks inject CRs into the gas, but the relative pressure of this component decreases towards halo centres as adiabatic compression favours the thermal over the CR pressure.