-
Cosmic acceleration is the most surprising cosmological discovery in many
decades. Testing and distinguishing among possible explanations requires
cosmological measurements of extremely high precision probing the full history
of cosmic expansion and structure growth and, ideally, compare and contrast
matter and relativistic tracers of the gravity potential. This program is one
of the defining objectives of the Wide-Field Infrared Survey Telescope
(WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010.
The WFIRST mission has the ability to improve these measurements by 1-2 orders
of magnitude compared to the current state of the art, while simultaneously
extending their redshift grasp, greatly improving control of systematic
effects, and taking a unified approach to multiple probes that provide
complementary physical information and cross-checks of cosmological results. We
describe in this annual report the activities of the Science Investigation Team
(SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017.
This team was selected by NASA in December 2015 in order to address the
stringent challenges of the WFIRST dark energy (DE) program through the
Project's formulation phase. This SIT has elected to jointly address Galaxy
Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the
fact that the imaging and spectroscopic elements of the HLS will be realized as
an integrated observing program, and they jointly impose requirements on
performance and operations. WFIRST is designed to be able to deliver a
definitive result on the origin of cosmic acceleration. It is not optimized for
Figure of Merit sensitivity but for control of systematic uncertainties and for
having multiple techniques each with multiple cross-checks. Our SIT work
focuses on understanding the potential systematics in the WFIRST DE
measurements.
-
Knowledge of the number density of H$\alpha$ emitting galaxies is vital for
assessing the scientific impact of the Euclid and WFIRST missions. In this work
we present predictions from a galaxy formation model, Galacticus, for the
cumulative number counts of H$\alpha$-emitting galaxies. We couple Galacticus
to three different dust attenuation methods and examine the counts using each
method. A $\chi^2$ minimisation approach is used to compare the model
predictions to observed galaxy counts and calibrate the dust parameters. We
find that weak dust attenuation is required for the Galacticus counts to be
broadly consistent with the observations, though the optimum dust parameters
return large values for $\chi^2$, suggesting that further calibration of
Galacticus is necessary. The model predictions are also consistent with
observed estimates for the optical depth and the H$\alpha$ luminosity function.
Finally we present forecasts for the redshift distributions and number counts
for two Euclid-like and one WFIRST-like survey. For a Euclid-like survey with
redshift range $0.9\leqslant z\leqslant 1.8$ and H$\alpha+{\rm [NII]}$ blended
flux limit of $2\times 10^{-16}{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^{-2}$ we
predict a number density between 3900--4800 galaxies per square degree. For a
WFIRST-like survey with redshift range $1\leqslant z\leqslant 2$ and blended
flux limit of $1\times 10^{-16}{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^{-2}$ we
predict a number density between 10400--15200 galaxies per square degree.
-
We present a comparison of nine galaxy formation models, eight
semi-analytical and one halo occupation distribution model, run on the same
underlying cold dark matter simulation (cosmological box of co-moving width
125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24\times 10^9 h^{-1}$
Msun) and the same merger trees. While their free parameters have been
calibrated to the same observational data sets using two approaches, they
nevertheless retain some 'memory' of any previous calibration that served as
the starting point (especially for the manually-tuned models). For the first
calibration, models reproduce the observed z = 0 galaxy stellar mass function
(SMF) within 3-{\sigma}. The second calibration extended the observational data
to include the z = 2 SMF alongside the z~0 star formation rate function, cold
gas mass and the black hole-bulge mass relation. Encapsulating the observed
evolution of the SMF from z = 2 to z = 0 is found to be very hard within the
context of the physics currently included in the models. We finally use our
calibrated models to study the evolution of the stellar-to-halo mass (SHM)
ratio. For all models we find that the peak value of the SHM relation decreases
with redshift. However, the trends seen for the evolution of the peak position
as well as the mean scatter in the SHM relation are rather weak and strongly
model dependent. Both the calibration data sets and model results are publicly
available.
-
We present the public release of the MultiDark-Galaxies: three distinct
galaxy catalogues derived from one of the Planck cosmology MultiDark
simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution
of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models
GALACTICUS, SAG, and SAGE to it. We compare the three models and their
conformity with observational data for a selection of fundamental properties of
galaxies like stellar mass function, star formation rate, cold gas fractions,
and metallicities - noting that they sometimes perform differently reflecting
model designs and calibrations. We have further selected galaxy subsamples of
the catalogues by number densities in stellar mass, cold gas mass, and star
formation rate in order to study the clustering statistics of galaxies. We show
that despite different treatment of orphan galaxies, i.e. galaxies that lost
their dark-matter host halo due to the finite mass resolution of the N-body
simulation or tidal stripping, the clustering signal is comparable, and
reproduces the observations in all three models - in particular when selecting
samples based upon stellar mass. Our catalogues provide a powerful tool to
study galaxy formation within a volume comparable to those probed by on-going
and future photometric and redshift surveys. All model data consisting of a
range of galaxy properties - including broad-band SDSS magnitudes - are
publicly available.
-
Dwarf galaxies are known to have remarkably low star formation efficiency due
to strong feedback. Adopting the dwarf galaxies of the Milky Way as a
laboratory, we explore a flexible semi-analytic galaxy formation model to
understand how the feedback processes shape the satellite galaxies of the Milky
Way. Using Markov-Chain Monte-Carlo, we exhaustively search a large parameter
space of the model and rigorously show that the general wisdom of strong
outflows as the primary feedback mechanism cannot simultaneously explain the
stellar mass function and the mass--metallicity relation of the Milky Way
satellites. An extended model that assumes that a fraction of baryons is
prevented from collapsing into low-mass halos in the first place can be
accurately constrained to simultaneously reproduce those observations. The
inference suggests that two different physical mechanisms are needed to explain
the two different data sets. In particular, moderate outflows with weak halo
mass dependence are needed to explain the mass--metallicity relation, and
prevention of baryons falling into shallow gravitational potentials of low-mass
halos (e.g. "pre-heating") is needed to explain the low stellar mass fraction
for a given subhalo mass.
-
Using a cosmological $N$-body simulation, we investigate the origin and
distribution of stars in the intracluster light (ICL) of a Fornax-like cluster.
In a dark matter only simulation we identify a halo which, at $z=0$, has $M_200
\simeq 4.1 \times 10^{13}M_{sun}$ and $r_{200} = 700kpc$, and replace infalling
subhalos with models that include spheroid and disc components. As they fall
into the cluster, the stars in some of these galaxies are stripped from their
hosts, and form the ICL. We consider the separate contributions to the ICL from
stars which originate in the haloes and the discs of the galaxies. We find that
disc ICL stars are more centrally concentrated than halo ICL stars. The
majority of the disc ICL stars are associated with one initially disc-dominated
galaxy that falls to the centre of the cluster and is heavily disrupted,
producing part of the cD galaxy. At radial distances greater than 200kpc, well
beyond the stellar envelope of the cD galaxy, stars formerly from the stellar
haloes of galaxies dominate the ICL. Therefore at large distances, the ICL
population is dominated by older stars.
-
We present a clustering comparison of 12 galaxy formation models (including
Semi-Analytic Models (SAMs) and Halo Occupation Distribution (HOD) models) all
run on halo catalogues and merger trees extracted from a single {\Lambda}CDM
N-body simulation. We compare the results of the measurements of the mean halo
occupation numbers, the radial distribution of galaxies in haloes and the
2-Point Correlation Functions (2PCF). We also study the implications of the
different treatments of orphan (galaxies not assigned to any dark matter
subhalo) and non-orphan galaxies in these measurements. Our main result is that
the galaxy formation models generally agree in their clustering predictions but
they disagree significantly between HOD and SAMs for the orphan satellites.
Although there is a very good agreement between the models on the 2PCF of
central galaxies, the scatter between the models when orphan satellites are
included can be larger than a factor of 2 for scales smaller than 1 Mpc/h. We
also show that galaxy formation models that do not include orphan satellite
galaxies have a significantly lower 2PCF on small scales, consistent with
previous studies. Finally, we show that the 2PCF of orphan satellites is
remarkably different between SAMs and HOD models. Orphan satellites in SAMs
present a higher clustering than in HOD models because they tend to occupy more
massive haloes. We conclude that orphan satellites have an important role on
galaxy clustering and they are the main cause of the differences in the
clustering between HOD models and SAMs.
-
The Astropy Project (http://astropy.org) is, in its own words, "a community
effort to develop a single core package for Astronomy in Python and foster
interoperability between Python astronomy packages." For five years this
project has been managed, written, and operated as a grassroots,
self-organized, almost entirely volunteer effort while the software is used by
the majority of the astronomical community. Despite this, the project has
always been and remains to this day effectively unfunded. Further, contributors
receive little or no formal recognition for creating and supporting what is now
critical software. This paper explores the problem in detail, outlines possible
solutions to correct this, and presents a few suggestions on how to address the
sustainability of general purpose astronomical software.
-
Many properties of the Milky Way's dark matter halo, including its mass
assembly history, concentration, and subhalo population, remain poorly
constrained. We explore the connection between these properties of the Milky
Way and its satellite galaxy population, especially the implication of the
presence of the Magellanic Clouds for the properties of the Milky Way halo.
Using a suite of high-resolution $N$-body simulations of Milky Way-mass halos
with a fixed final Mvir ~ 10^{12.1}Msun, we find that the presence of
Magellanic Cloud-like satellites strongly correlates with the assembly history,
concentration, and subhalo population of the host halo, such that Milky
Way-mass systems with Magellanic Clouds have lower concentration, more rapid
recent accretion, and more massive subhalos than typical halos of the same
mass. Using a flexible semi-analytic galaxy formation model that is tuned to
reproduce the stellar mass function of the classical dwarf galaxies of the
Milky Way with Markov-Chain Monte-Carlo, we show that adopting host halos with
different mass-assembly histories and concentrations can lead to different
best-fit models for galaxy-formation physics, especially for the strength of
feedback. These biases arise because the presence of the Magellanic Clouds
boosts the overall population of high-mass subhalos, thus requiring a different
stellar-mass-to-halo-mass ratio to match the data. These biases also lead to
significant differences in the mass--metallicity relation, the kinematics of
low-mass satellites, the number counts of small satellites associated with the
Magellanic Clouds, and the stellar mass of Milky Way itself. Observations of
these galaxy properties can thus provide useful constraints on the properties
of the Milky Way halo.
-
We develop a simple analytical model that tracks galactic metallicities
governed by star formation and feedback to gain insight from the observed
galaxy stellar mass-metallicity relations over a large range of stellar masses
and redshifts. The model reveals the following implications of star formation
and feedback processes in galaxy formation. First, the observed metallicity
relations provide a stringent upper limit for the averaged outflow mass-loading
factors of local galaxies, which is ~20 for M_*~10^9Msun galaxies and
monotonically decreases to ~1 for M_*~10^{11}Msun galaxies. Second, the
inferred upper-limit for the outflow mass-loading factor sensitively depends on
whether the outflow is metal-enriched with respect to the ISM metallicity. If
half of the metals ejected from SNe leave the galaxy in metal-enriched winds,
the outflow mass-loading factor for galaxies at any mass can barely be higher
than ~10, which puts strong constraints on galaxy formation models. Third, the
relatively lower stellar-phase to gas-phase metallicity ratio for lower-mass
galaxies indicate that low-mass galaxies are still rapidly enriching their
metallicities in recent times, while high-mass galaxies are more settled, which
seems to show a downsizing effect in the metallicity evolution of galaxies. The
analysis presented in the paper demonstrates the importance of accurate
measurements of galaxy metallicities and the cold gas fraction of galaxies at
different redshifts for constraining star formation and feedback processes, and
demonstrates the power of these relations in constraining the physics of galaxy
formation.
-
We present a comparison of 14 galaxy formation models: 12 different
semi-analytical models and 2 halo-occupation distribution models for galaxy
formation based upon the same cosmological simulation and merger tree
information derived from it. The participating codes have proven to be very
successful in their own right but they have all been calibrated independently
using various observational data sets, stellar models, and merger trees. In
this paper we apply them without recalibration and this leads to a wide variety
of predictions for the stellar mass function, specific star formation rates,
stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter
is much larger than seen in previous comparison studies primarily because the
codes have been used outside of their native environment within which they are
well tested and calibrated. The purpose of the `nIFTy comparison of galaxy
formation models' is to bring together as many different galaxy formation
modellers as possible and to investigate a common approach to model
calibration. This paper provides a unified description for all participating
models and presents the initial, uncalibrated comparison as a baseline for our
future studies where we will develop a common calibration framework and address
the extent to which that reduces the scatter in the model predictions seen
here.
-
We present a simulation setup for studying the dynamical and chemical
evolution of the intracluster medium (ICM) and analyze a sample of 12 galaxy
clusters that are diverse both kinetically (pre-merger, merging, virialized)
and in total mass (M vir = 1.17 x 10^14 - 1.06 x 10^15 M). We analyzed the
metal mass fraction in the ICM as a function of redshift and discuss radial
trends as well as projected 2D metallicity maps. The setup combines high mass
resolution N-body simulations with the semi-analytical galaxy formation model
Galacticus for consistent treatment of the subgrid physics (such as galactic
winds and ram-pressure stripping) in the cosmological hydrodynamical
simulations. The interface between Galacticus and the hydro simulation of the
ICM with FLASH is discussed with respect to observations of star formation rate
histories, radial star formation trends in galaxy clusters, and the metallicity
at different redshifts. As a test for the robustness of the wind model, we
compare three prescriptions from different approaches. For the wind model
directly taken from Galacticus, we find mean ICM metallicities between 0.2 -
0.8Z within the inner 1Mpc at z = 0. The main contribution to the metal mass
fraction comes from galactic winds. The outflows are efficiently mixed in the
ICM, leading to a steady homogenization of metallicities until ram-pressure
stripping becomes effective at low redshifts. We find a very peculiar and yet
common drop in metal mass fractions within the inner ~200kpc of the cool cores,
which is due to a combination of wind suppression by outer pressure within our
model and a lack of mixing after the formation of these dense regions.
-
We study the oscillations and stability of self-gravitating cylindrically
symmetric fluid systems and collisionless systems. This is done by studying
small perturbations to the equilibrium system and finding the normal modes,
using methods similar to those used in astroseismology. We find that there is a
single sequence of purely radial modes that become unstable if the adiabatic
exponent is less than 1. Nonradial modes can be divided into p modes, which are
stable and pressure-driven, and g modes, which are are gravity driven. The g
modes become unstable if the adiabatic exponent is greater than the polytrope
index. These modes are analogous to the modes of a spherical star, but their
behavior is somewhat different because a cylindrical geometry has less symmetry
than a spherical geometry. This implies that perturbations are classified by a
radial quantum number, an azimuthal quantum number, and wavelength in the z
direction, which can become arbitrarily large. We find that decreasing this
wavelength increases the frequency of stable modes and increases the growth
rate of unstable modes. We use use variational arguments to demonstrate that
filaments of collisionless matter with ergodic distribution functions are
stable to purely radial perturbations, and that filaments with ergodic
power-law distribution functions are stable to all perturbations.
-
Well-motivated particle physics theories predict the existence of particles
(such as sterile neutrinos) which acquire non-negligible thermal velocities in
the early universe. These particles could behave as warm dark matter (WDM) and
generate a small-scale cutoff in the linear density power spectrum which scales
approximately inversely with the particle mass. If this mass is of order a keV,
the cutoff occurs on the scale of dwarf galaxies. Thus, in WDM models the
abundance of small galaxies, such as the satellites that orbit in the halo of
the Milky Way, depends on the mass of the warm particle. The abundance also
scales with the mass of the host galactic halo. We use the \galform
semi-analytic model of galaxy formation to calculate the properties of galaxies
in universes in which the dark matter is warm. Using this method, we can
compare the predicted satellite luminosity functions to the observed data for
the Milky Way dwarf spheroidals, and determine a lower bound on the thermally
produced WDM particle mass. This depends strongly on the value of the Milky Way
halo mass and, to some extent, on the baryonic physics assumed; we examine both
of these dependencies. For our fiducial model we find that for a particle mass
of 3.3 keV (the 2$\sigma$ lower limit found by Viel et al. from a recent
analysis of the Lyman-$\alpha$ forest) the Milky Way halo mass is required to
be $> 1.4 \times 10^{12}$ \msun. For this same fiducial model, we also find
that all WDM particle masses are ruled out (at 95% confidence) if the halo of
the Milky Way has a mass smaller than $1.1 \times 10^{12}$ \msun, while if the
mass of the Galactic halo is greater than 1.8 $\times 10^{12}$ \msun, only WDM
particle masses larger than 2 keV are allowed.
-
We present a spectro-photometric survey of 2522 extragalactic globular
clusters (GCs) around twelve early-type galaxies, nine of which have not been
published previously. Combining space-based and multi-colour wide field
ground-based imaging, with spectra from the Keck DEIMOS instrument, we obtain
an average of 160 GC radial velocities per galaxy, with a high velocity
precision of 15 km/s per GC. After studying the photometric properties of the
GC systems, such as their spatial and colour distributions, we focus on the
kinematics of metal-poor (blue) and metal-rich (red) GC subpopulations to an
average distance of ~8 effective radii from the galaxy centre. Our results show
that for some systems the bimodality in GC colour is also present in GC
kinematics. The kinematics of the red GC subpopulations are strongly coupled
with the host galaxy stellar kinematics. The blue GC subpopulations are more
dominated by random motions, especially in the outer regions, and decoupled
from the red GCs. Peculiar GC kinematic profiles are seen in some galaxies: the
blue GCs in NGC 821 rotate along the galaxy minor axis, whereas the GC system
of the lenticular galaxy NGC 7457 appears to be strongly rotation supported in
the outer region. We supplement our galaxy sample with data from the literature
and carry out a number of tests to study the kinematic differences between the
two GC subpopulations. We confirm that the GC kinematics are coupled with the
host galaxy properties and find that the velocity kurtosis and the slope of
their velocity dispersion profiles is different between the two GC
subpopulations in more massive galaxies.
-
Using data from the COSMOS survey, we perform the first joint analysis of
galaxy-galaxy weak lensing, galaxy spatial clustering, and galaxy number
densities. Carefully accounting for sample variance and for scatter between
stellar and halo mass, we model all three observables simultaneously using a
novel and self-consistent theoretical framework. Our results provide strong
constraints on the shape and redshift evolution of the stellar-to-halo mass
relation (SHMR) from z=0.2 to z=1. At low stellar mass, we find that halo mass
scales as Mh M*^0.46 and that this scaling does not evolve significantly with
redshift to z=1. We show that the dark-to-stellar ratio, Mh/M*, varies from low
to high masses, reaching a minimum of Mh/M*~27 at M*=4.5x10^10 Msun and
Mh=1.2x10^12 Msun. This minimum is important for models of galaxy formation
because it marks the mass at which the accumulated stellar growth of the
central galaxy has been the most efficient. We describe the SHMR at this
minimum in terms of the "pivot stellar mass", M*piv, the "pivot halo mass",
Mhpiv, and the "pivot ratio", (Mh/M*)piv. Thanks to a homogeneous analysis of a
single data set, we report the first detection of mass downsizing trends for
both Mhpiv and M*piv. The pivot stellar mass decreases from
M*piv=5.75+-0.13x10^10 Msun at z=0.88 to M*piv=3.55+-0.17x10^10 Msun at z=0.37.
Intriguingly, however, the corresponding evolution of Mhpiv leaves the pivot
ratio constant with redshift at (Mh/M*)piv~27. We use simple arguments to show
how this result raises the possibility that star formation quenching may
ultimately depend on Mh/M* and not simply Mh, as is commonly assumed. We show
that simple models with such a dependence naturally lead to downsizing in the
sites of star formation. Finally, we discuss the implications of our results in
the context of popular quenching models, including disk instabilities and AGN
feedback.
-
While the dark matter content within the most massive giant and smallest
dwarf galaxies has been probed -- spanning a range of over one million in mass
-- an important observational gap remains for galaxies of intermediate mass.
This gap covers K band magnitudes of approximately -16 > M_K > -18 (for which
dwarf galaxies have B--K ~ 2). On the high mass side of the gap are dwarf
elliptical (dE) galaxies, that are dominated by stars in their inner regions.
While the low mass side includes dwarf spheroidal (dSph) galaxies that are dark
matter-dominated and ultra compact dwarf (UCD) objects that are star-dominated.
Evolutionary pathways across the gap have been suggested but remain largely
untested because the `gap' galaxies are faint, making dynamical measurements
very challenging. With long exposures on the Keck telescope using the ESI
instrument we have succeeded in bridging this gap by measuring the dynamical
mass for five dwarf galaxies with M_K ~ -17.5 (M_B ~ --15.5). With the
exception of our brightest dwarf galaxy, they possess relatively flat velocity
dispersion profiles of around 20 km/s. By examining their 2D scaling relations
and 3D fundamental manifold, we found that the sizes and velocity dispersions
of these gap galaxies reveal continuous trends from dE to dSph galaxies. We
conclude that low-luminosity dwarf elliptical galaxies are dominated by stars,
not by dark matter, within their half light radii. This finding can be
understood if internal feedback processes are operating most efficiently in gap
galaxies, gravitationally heating the centrally-located dark matter to larger
radii. Whereas external environmental processes, which can strip away stars,
have a greater influence on dSph galaxies resulting in their higher dark matter
fractions. Abridged.
-
Models of disk galaxy formation commonly predict the existence of an extended
reservoir of accreted hot gas surrounding massive spirals at low redshift. As a
test of these models, we use X-ray and H-alpha data of the two massive,
quiescent edge-on spirals NGC 5746 and NGC 5170 to investigate the amount and
origin of any hot gas in their halos. Contrary to our earlier claim, the
Chandra analysis of NGC 5746, employing more recent calibration data, does not
reveal any significant evidence for diffuse X-ray emission outside the optical
disk, with a 3-sigma upper limit to the halo X-ray luminosity of 4e39 erg/s. An
identical study of the less massive NGC 5170 also fails to detect any
extraplanar X-ray emission. By extracting hot halo properties of disk galaxies
formed in cosmological hydrodynamical simulations, we compare these results to
expectations for cosmological accretion of hot gas by spirals. For Milky
Way-sized galaxies, these high-resolution simulations predict hot halo X-ray
luminosities which are lower by a factor of ~2 compared to our earlier results
reported by Toft et al. (2002). We find the new simulation predictions to be
consistent with our observational constraints for both NGC 5746 and NGC 5170,
while also confirming that the hot gas detected so far around more actively
star-forming spirals is in general probably associated with stellar activity in
the disk. Observational results on quiescent disk galaxies at the high-mass end
are nevertheless providing powerful constraints on theoretical predictions, and
hence on the assumed input physics in numerical studies of disk galaxy
formation and evolution.
-
There is a vast menagerie of plausible candidates for the constituents of
dark matter, both within and beyond extensions of the Standard Model of
particle physics. Each of these candidates may have scattering (and other)
cross section properties that are consistent with the dark matter abundance,
BBN, and the most scales in the matter power spectrum; but which may have
vastly different behavior at sub-galactic "cutoff" scales, below which dark
matter density fluctuations are smoothed out. The only way to quantitatively
measure the power spectrum behavior at sub-galactic scales at distances beyond
the local universe, and indeed over cosmic time, is through probes available in
multiply imaged strong gravitational lenses. Gravitational potential
perturbations by dark matter substructure encode information in the observed
relative magnifications, positions, and time delays in a strong lens. Each of
these is sensitive to a different moment of the substructure mass function and
to different effective mass ranges of the substructure. The time delay
perturbations, in particular, are proving to be largely immune to the
degeneracies and systematic uncertainties that have impacted exploitation of
strong lenses for such studies. There is great potential for a coordinated
theoretical and observational effort to enable a sophisticated exploitation of
strong gravitational lenses as direct probes of dark matter properties. This
opportunity motivates this white paper, and drives the need for: a) strong
support of the theoretical work necessary to understand all astrophysical
consequences for different dark matter candidates; and b) tailored
observational campaigns, and even a fully dedicated mission, to obtain the
requisite data.
-
We present new measurements of the evolution in the Lyman break galaxy (LBG)
population between z~4 and z~6. By utilizing the extensive multiwavelength
datasets available in the GOODS fields, we identify 2443 B, 506 V, and 137
i'-band dropout galaxies likely to be at z~4, 5, and 6. With the goal of
understanding the duration of typical star formation episodes in galaxies at
z>4, we examine the distribution of stellar masses and ages as a function of
cosmic time. We find that at a fixed rest-UV luminosity, the average stellar
masses and ages of galaxies do not increase significantly between z~6 and 4. In
order to maintain this near equilibrium in the average properties of high
redshift LBGs, we argue that there must be a steady flux of young,
newly-luminous objects at each successive redshift. When considered along with
the short duty cycles inferred from clustering measurements, these results may
suggest that galaxies are undergoing star formation episodes lasting only
several hundred million years. In contrast to the unchanging relationship
between the average stellar mass and rest-UV luminosity, we find that the
number density of massive galaxies increases considerably with time over 4<z<6.
Given this rapid increase of UV luminous massive galaxies, we explore the
possibility that a significant fraction of massive (M*>1e11 Msun) z~2-3 distant
red galaxies (DRGs) were in part assembled in an LBG phase at earlier times.
Integrating the growth in the stellar mass function of actively forming LBGs
over 4<z<6 down to z~2, we find that z>3 LBGs could have contributed
significantly to the quiescent DRG population, indicating that the intense
star-forming systems probed by current sub-millimeter observations are not the
only route toward the assembly of DRGs at z~2.
-
Both preheating of the intergalactic medium and radiative cooling of low
entropy gas have been proposed to explain the deviation from self-similarity in
the cluster L_x-T_x relation and the observed entropy floor in these systems.
However, severe overcooling of gas in groups is necessary for radiative cooling
alone to explain the observations. Non-gravitational entropy injection must
therefore still be important in these systems. We point out that on scales of
groups and below, gas heated to the required entropy floor cannot cool in a
Hubble time, regardless of its subsequent adiabatic compression. Preheating
therefore shuts off the gas supply to galaxies, and should be an important
global feedback mechanism for galaxy formation. Constraints on global gas
cooling can be placed from the joint evolution of the comoving star formation
rate and neutral gas density. Preheating at high redshift can be ruled out;
however the data does not rule out passive gas consumption without inflow since
z~2. Since for preheated gas t_cool > t_dyn, we speculate that preheating could
play a role in determining the Hubble sequence: at a given mass scale, high
sigma peaks in the density field collapse early to form ellipticals, while low
sigma peaks collapse late and quiescently accrete preheated gas to form
spirals. The entropy produced by large scale shock-heating of the intergalatic
medium is significant only at late times, z<1, and cannot produce these
effects.
-
We develop a hybrid galaxy formation model which uses outputs from an N-body
simulation to follow the merger histories (or ``merger trees'') of dark matter
halos and treats baryonic processes, such as the cooling of gas within halos
and subsequent star formation, using the semi-analytic model of Cole et al. We
compare this hybrid model to an otherwise identical model which utilises merger
tree realisations generated by a Monte-Carlo algorithm and find that, apart
from the limited mass resolution imposed by the N-body particle mass, the only
significant differences between the models are due to the known discrepancy
between the distribution of halo progenitor masses predicted by extended Press
Schechter theory and that found in N-body simulations. We investigate the
effect of limited mass resolution on the hybrid model by comparing to a purely
semi-analytic model with greatly improved mass resolution. We find that the
mass resolution of the simulation we use, which has a particle mass of
1.4x10^10 Msolar/h, is insufficient to produce a reasonable luminosity function
for galaxies with magnitudes in the b_J band fainter than -17.
-
We compare the results of two techniques used to calculate the evolution of
cooling gas during galaxy formation: Smooth Particle Hydrodynamics (SPH)
simulations and semi-analytic modelling. We improve upon the earlier
statistical comparison of Benson et al. by taking halo merger histories from
the dark matter component of the SPH simulation, which allows us to compare the
evolution of galaxies on an object-by-object basis in the two treatments. We
use a ``stripped-down'' version of the semi-analytic model described by Helly
et al. which includes only shock heating and radiative cooling of gas and which
is adjusted to mimic the resolution and other parameters of a comparison SPH
simulation as closely as possible. We compare the total mass of gas that cools
in halos of different mass as a function of redshift as well as the masses and
spatial distribution of individual ``galaxies.'' At redshift z=0, the cooled
gas mass in well-resolved halos agrees remarkably well (to better than ~20%) in
the SPH simulation and stripped-down semi-analytic model. At high redshift,
resolution effects in the simulation become increasingly important and, as a
result, more gas tends to cool in low mass halos in the SPH simulation than in
the semi-analytic model. The cold gas mass function of individual galaxies in
the two treatments at z=0 also agrees very well and, when the effects of
mergers are accounted for, the masses of individual galaxies and their 2-point
correlation functions are also in excellent agreement in the two treatments.
Thus, our comparison confirms and extends the earlier conclusion of Benson et
al. that SPH simulations and semi-analytic models give consistent results for
the evolution of cooling galactic gas.
-
We present predictions for two statistical measures of the hydrogen
reionization process at high redshift. The first statistic is the number of
neutral segments identified in spectra of high redshift QSOs as a function of
their length. The second is the cross-correlation of neutral regions with
possible sources of ionizing radiation. These independent probes are sensitive
to the topology of the ionized regions. If reionization proceeded from high to
low density regions then the cross-correlation will be negative, while if voids
were ionized first then we expect a positive correlation and a relatively small
number of long neutral segments. We test the sensitivity of these statistics
for reionization by stars in high redshift galaxies. The flux of ionizing
radiation emitted from stars is estimated by identifying galaxies in an N-body
simulation using a semi-analytic galaxy formation model. The spatial
distribution of ionized gas is traced in various models for the propagation of
the ionization fronts. A model with ionization proceeding from high to low
density regions is consistent with the observations of Becker et al. (2001),
while models in which ionization begins in the lowest density regions appear to
be inconsistent with the present data.
-
The main ingredients of recent semi-analytic models of galaxy formation are
summarised. We present predictions for the galaxy clustering properties of a
well specified LCDM model whose parameters are constrained by observed local
galaxy properties. We present preliminary predictions for evolution of
clustering that can be probed with deep pencil beam surveys.