19 research outputs found
Generation of mock tidal streams
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 pericenter, 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.Comment: 20 pages, 11 figures, 1 table; submitted to MNRA
A New Model For Including Galactic Winds in Simulations of Galaxy Formation II: Implementation of PhEW in Cosmological Simulations
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. (2020) 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, metals,
and momentum 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 conventional wind simulations with the same resolution, our PhEW simulations
produce similar galaxy stellar mass functions at z ≥ 1 but are in better agreement with low-redshift observations
at M∗ < 1011M⊙ because PhEW particles shed mass to the CGM before escaping low mass halos. 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 halos 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
A New Model For Including Galactic Winds in Simulations of Galaxy Formation II: Implementation of PhEW in Cosmological Simulations
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. (2020) 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, metals, and momentum 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 conventional wind simulations with the same
resolution, our PhEW simulations produce similar galaxy stellar mass functions
at but are in better agreement with low-redshift observations at because PhEW particles shed mass to the CGM before escaping
low mass halos. 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 halos 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.Comment: 23 pages, 17 figures, MNRAS accepte
The Impact of Wind Scalings on Stellar Growth and the Baryon Cycle in Cosmological Simulations
Many phenomenologically successful cosmological simulations employ kinetic winds to model galactic outflows. Yet systematic studies of how variations in kinetic wind scalings might alter observable galaxy properties are rare. Here we employ GADGET-3 simulations to study how the baryon cycle, stellar mass function, and other galaxy and CGM predictions vary as a function of the assumed outflow speed and the scaling of the mass-loading factor with velocity dispersion. We design our fiducial model to reproduce the measured wind properties at 25 per cent of the virial radius from the Feedback In Realistic Environments simulations
ELUCID. VII. Using constrained hydro simulations to explore the gas component of the cosmic web
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 strategie
ELUCID IV: Galaxy Quenching and its Relation to Halo Mass, Environment, and Assembly Bias
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 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
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 the difference between the two
populations found previously mainly arises from the fact that centrals and
satellites of the same stellar mass reside, on average, in halos of different
mass. After removing these halo-mass and stellar-mass effects, 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.Comment: 21 pages, 16 figures, submitted to Ap
The robustness of cosmological hydrodynamic simulation predictions to changes in numerics and cooling physics
We test and improve the numerical schemes in our smoothed particle
hydrodynamics (SPH) code for cosmological simulations, including the
pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a
refined timestep criterion, and metal-line cooling that accounts for
photoionisation in the presence of a recently refined Haardt \& Madau (2012)
model of the ionising background. The PESPH algorithm effectively removes the
artificial surface tension present in the traditional SPH formulation, and in
our test simulations it produces better qualitative agreement with mesh-code
results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set
of cosmological simulations, we examine many of the quantities we have studied
in previous work. Results for galaxy stellar and HI mass functions, star
formation histories, galaxy scaling relations, and statistics of the Ly
forest are robust to the changes in numerics and microphysics. As in our
previous simulations, cold gas accretion dominates the growth of high-redshift
galaxies and of low mass galaxies at low redshift, and recycling of winds
dominates the growth of massive galaxies at low redshift. However, the PESPH
simulation removes spurious cold clumps seen in our earlier simulations, and
the accretion rate of hot gas increases by up to an order of magnitude at some
redshifts. The new numerical model also influences the distribution of metals
among gas phases, leading to considerable differences in the statistics of some
metal absorption lines, most notably NeVIII.Comment: 29 pages, 25 figures, accepted by MNRA
nIFTy galaxy cluster simulations – II. Radiative models
We have simulated the formation of a massive galaxy cluster (M = 1.110) in a CDM universe using
10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling
hydrodynamics with full radiative subgrid physics. These codes include
Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH
schemes, adaptive mesh and moving mesh codes. Our goal is to study the
consistency between simulated clusters modeled with different radiative
physical implementations - such as cooling, star formation and AGN feedback. We
compare images of the cluster at , global properties such as mass, and
radial profiles of various dynamical and thermodynamical quantities. We find
that, with respect to non-radiative simulations, dark matter is more centrally
concentrated, the extent not simply depending on the presence/absence of AGN
feedback. The scatter in global quantities is substantially higher than for
non-radiative runs. Intriguingly, adding radiative physics seems to have washed
away the marked code-based differences present in the entropy profile seen for
non-radiative simulations in Sembolini et al. (2015): radiative physics +
classic SPH can produce entropy cores. Furthermore, the inclusion/absence of
AGN feedback is not the dividing line -as in the case of describing the stellar
content- for whether a code produces an unrealistic temperature inversion and a
falling central entropy profile. However, AGN feedback does strongly affect the
overall stellar distribution, limiting the effect of overcooling and reducing
sensibly the stellar fraction.Comment: 20 pages, 13 figures, submitted to MNRA
nIFTy galaxy cluster simulations - IV. Quantifying the influence of baryons on halo properties
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-matteronly
(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 R200. However, we see larger differences for quantities within R2500,
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