We examine the quenched fraction of central and satellite galaxies as a function of galaxy stellar mass, halo mass, and the matter density of their large-scale environment. Matter densities are ...inferred from our ELUCID simulation, a constrained simulation of the local universe sampled by SDSS, while halo masses and central/satellite classification are taken from the galaxy group catalog of Yang et al. The quenched fraction for the total population increases systematically with the three quantities. We find that the "environmental quenching efficiency," which quantifies the quenched fraction as a function of halo mass, is independent of stellar mass. And this independence is the origin of the stellar mass independence of density-based quenching efficiency found in previous studies. Considering centrals and satellites separately, we find that the two populations follow similar correlations of quenching efficiency with halo mass and stellar mass, suggesting that they have experienced similar quenching processes in their host halo. We demonstrate that satellite quenching alone cannot account for the environmental quenching efficiency of the total galaxy population, and that the difference between the two populations found previously arises mainly from the fact that centrals and satellites of the same stellar mass reside, on average, in halos of different mass. After removing these effects of halo mass and stellar mass, there remains a weak, but significant, residual dependence on environmental density, which is eliminated when halo assembly bias is taken into account. Our results therefore indicate that halo mass is the prime environmental parameter that regulates the quenching of both centrals and satellites.
Abstract We conducted an investigation of the Coma cluster of galaxies by running a series of constrained hydrodynamic simulations with GIZMO-SIMBA and GADGET-3 based on initial conditions ...reconstructed from the SDSS survey volume in the ELUCID project. We compared simulation predictions and observations for galaxies, intracluster medium (ICM) and intergalactic medium (IGM) in and around the Coma cluster to constrain galaxy formation physics. Our results demonstrate that this type of constrained investigation allows us to probe in more detail the implemented physical processes, because the comparison between simulations and observations is free of cosmic variance and hence can be conducted in a “one-to-one” manner. We found that an increase in the earlier star formation rate and the supernova feedback of the original GIZMO-SIMBA model is needed to match observational data on stellar, interstellar medium, and ICM metallicity. The simulations without active galactic nucleus (AGN) feedback can well reproduce the observational ICM electron density, temperature, and entropy profiles, ICM substructures, and the IGM temperature–density relation, while the ones with AGN feedback usually fail. However, one requires something like AGN feedback to reproduce a sufficiently large population of quiescent galaxies, particularly in low-density regions. The constrained simulations of the Coma cluster thus provide a test bed to understand processes that drive galaxy formation and evolution.
We have simulated the formation of a galaxy cluster in a Λ cold dark matter universe using 13 different codes modelling only gravity and non-radiative hydrodynamics (ramses, ART, arepo, hydra and ...nine incarnations of gadget). This range of codes includes particle-based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various gadget implementations span classic and modern smoothed particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at z = 0, global properties such as mass and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ramses, art and arepo form extended entropy cores in the gas with rising central gas temperatures. Those codes employing classic SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid-based methods.
Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 ...state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern smoothed particle hydrodynamics (SPH) approaches. For each code represented we have a dark-matter-only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency ⪅20 per cent) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R
200. However, we see larger differences for quantities within R
2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include active galactic nucleus feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers.
Generation of mock tidal streams Fardal, Mark A; Huang, Shuiyao; Weinberg, Martin D
Monthly Notices of the Royal Astronomical Society,
09/2015, Letnik:
452, Številka:
1
Journal Article
Recenzirano
Odprti dostop
In this paper, we discuss a method for the generation of mock tidal streams. Using an ensemble of simulations in an isochrone potential where the actions and frequencies are known, we derive an ...empirical recipe for the evolving satellite mass and the corresponding mass-loss rate, and the ejection conditions of the stream material. The resulting stream can then be quickly generated either with direct orbital integration, or by using the action-angle formalism. The model naturally produces streaky features within the stream. These are formed due to the radial oscillation of the progenitor and the bursts of stars emitted near pericentre, rather than clumping at particular oscillation phases as sometimes suggested. When detectable, these streaky features are a reliable diagnostic for the stream's direction of motion and encode other information on the progenitor and its orbit. We show several tests of the recipe in alternate potentials, including a case with a chaotic progenitor orbit which displays a marked effect on the width of the stream. Although the specific ejection recipe may need adjusting when elements such as the orbit or satellite density profile are changed significantly, our examples suggest that model tidal streams can be quickly and accurately generated by models of this general type for use in Bayesian sampling.
ABSTRACT
The propagation and evolution of cold galactic winds in galactic haloes is crucial to galaxy formation models. However, modelling of this process in hydrodynamic simulations of galaxy ...formation is oversimplified owing to a lack of numerical resolution and often neglects critical physical processes such as hydrodynamic instabilities and thermal conduction. We propose an analytic model, Physically Evolved Winds, that calculates the evolution of individual clouds moving supersonically through a uniform ambient medium. Our model reproduces predictions from very high resolution cloud-crushing simulations that include isotropic thermal conduction over a wide range of physical conditions. We discuss the implementation of this model into cosmological hydrodynamic simulations of galaxy formation as a subgrid prescription to model galactic winds more robustly both physically and numerically.
The sources that drove cosmological reionization left clues regarding their identity in the slope and inhomogeneity of the ultraviolet ionizing background (UVB): bright quasars (QSOs) generate a hard ...UVB with predominantly large-scale fluctuations while Population II stars generate a softer one with smaller scale fluctuations. Metal absorbers probe the UVB's slope because different ions are sensitive to different energies. Likewise, they probe spatial fluctuations because they originate in regions where a galaxy-driven UVB is harder and more intense. We take a first step towards studying the reionization-epoch UVB's slope and inhomogeneity by comparing observations of 12 metal absorbers at z ∼ 6 versus predictions from a cosmological hydrodynamic simulation using three different UVBs: a soft, spatially inhomogeneous ‘galaxies+QSOs’ UVB; a homogeneous ‘galaxies+QSOs’ UVB, and a ‘QSOs-only’ model. All UVBs reproduce the observed column density distributions of
$\mathrm{C\,\small {ii}}$
,
$\mathrm{Si\,\small {iv}}$
, and
$\mathrm{C\,\small {iv}}$
reasonably well although high-column, high-ionization absorbers are underproduced, reflecting numerical limitations. With upper limits treated as detections, only a soft, fluctuating UVB reproduces both the observed
$\mathrm{Si\,\small {iv}}/\mathrm{C\,\small {iv}}$
and
$\mathrm{C\,\small {ii}}/\mathrm{C\,\small {iv}}$
distributions. The QSOs-only UVB overpredicts both
$\mathrm{C\,\small {iv}}/\mathrm{C\,\small {ii}}$
and
$\mathrm{C\,\small {iv}}/\mathrm{Si\,\small {iv}}$
, indicating that it is too hard. The Haardt & Madau (2012) UVB underpredicts
$\mathrm{C\,\small {iv}}/\mathrm{Si\,\small {iv}}$
, suggesting that it lacks amplifications near galaxies. Hence current observations prefer a soft, fluctuating UVB as expected from a predominantly Population II background although they cannot rule out a harder one. Future observations probing a factor of 2 deeper in metal column density will distinguish between the soft, fluctuating and QSOs-only UVBs.
Abstract
Using reconstructed initial conditions in the Sloan Digital Sky Survey (SDSS) survey volume, we carry out constrained hydrodynamic simulations in three regions representing different types ...of the cosmic web: the Coma cluster of galaxies; the SDSS Great Wall; and a large low-density region at
z
∼ 0.05. These simulations, which include star formation and stellar feedback but no active galactic nucleus formation and feedback, are used to investigate the properties and evolution of intergalactic and intracluster media. About half of the warm-hot intergalactic gas is associated with filaments in the local cosmic web. Gas in the outskirts of massive filaments and halos can be heated significantly by accretion shocks generated by mergers of filaments and halos, respectively, and there is a tight correlation between the gas temperature and the strength of the local tidal field. The simulations also predict some discontinuities associated with shock fronts and contact edges, which can be tested using observations of the thermal Sunyaev–Zel’dovich effect and X-rays. A large fraction of the sky is covered by Ly
α
and O
vi
absorption systems, and most of the O
vi
systems and low-column-density H
i
systems are associated with filaments in the cosmic web. The constrained simulations, which follow the formation and heating history of the observed cosmic web, provide an important avenue to interpret observational data. With full information about the origin and location of the cosmic gas to be observed, such simulations can also be used to develop observational strategies.
ABSTRACT
Although galactic winds play a critical role in regulating galaxy formation, hydrodynamic cosmological simulations do not resolve the scales that govern the interaction between winds and the ...ambient circumgalactic medium (CGM). We implement the Physically Evolved Wind (PhEW) model of Huang et al. in the gizmo hydrodynamics code and perform test cosmological simulations with different choices of model parameters and numerical resolution. PhEW adopts an explicit subgrid model that treats each wind particle as a collection of clouds that exchange mass and metals with their surroundings and evaporate by conduction and hydrodynamic instabilities as calibrated on much higher resolution cloud scale simulations. In contrast to a conventional wind algorithm, we find that PhEW results are robust to numerical resolution and implementation details because the small scale interactions are defined by the model itself. Compared to our previous wind simulations with the same resolution, our PhEW simulations are in better agreement with low-redshift galactic stellar mass functions at M* < 1011M⊙ because PhEW particles shed mass to the CGM before escaping low mass haloes. PhEW radically alters the CGM metal distribution because PhEW particles disperse metals to the ambient medium as their clouds dissipate, producing a CGM metallicity distribution that is skewed but unimodal and is similar between cold and hot gas. While the temperature distributions and radial profiles of gaseous haloes are similar in simulations with PhEW and conventional winds, these changes in metal distribution will affect their predicted UV/X-ray properties in absorption and emission.
The reionization of carbon Finlator, Kristian; Thompson, Robert; Huang, Shuiyao ...
Monthly Notices of the Royal Astronomical Society,
03/2015, Letnik:
447, Številka:
3
Journal Article
Recenzirano
Odprti dostop
Observations suggest that C ii was more abundant than C iv in the intergalactic medium towards the end of the hydrogen reionization epoch (z ~ 6). This transition provides a unique opportunity to ...study the enrichment history of intergalactic gas and the growth of the ionizing ultraviolet background (UVB) at early times. We study how carbon absorption evolves from z = 10 to 5 using a cosmological hydrodynamic simulation that includes a self-consistent multifrequency UVB as well as a well-constrained model for galactic outflows to disperse metals. Our predicted UVB is within ~2-4 times of that from Haardt & Madau, which is fair agreement given the uncertainties. Nonetheless, we use a calibration in post-processing to account for Lyman ... forest measurements while preserving the predicted spectral slope and inhomogeneity. The UVB fluctuates spatially in such a way that it always exceeds the volume average in regions where metals are found. This implies both that a spatially uniform UVB is a poor approximation and that metal absorption is not sensitive to the epoch when H ii regions overlap globally even at column densities of 10... cm... We find, consistent with observations, that the C ii mass fraction drops to low redshift while C iv rises owing the combined effects of a growing UVB and continued addition of carbon in low-density regions. This is mimicked in absorption statistics, which broadly agree with observations at z = 6-3 while predicting that the absorber column density distributions rise steeply to the lowest observable columns. Our model reproduces the large observed scatter in the number of low-ionization absorbers per sightline, implying that the scatter does not indicate a partially neutral Universe at z ~ 6. (ProQuest: ... denotes formulae/symbols omitted.)