• 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 Problem (1610.03159)

    Demitri Muna, Michael Alexander, Alice Allen, Richard Ashley, Daniel Asmus, Ruyman Azzollini, Michele Bannister, Rachael Beaton, Andrew Benson, G. Bruce Berriman, Maciej Bilicki, Peter Boyce, Joanna Bridge, Jan Cami, Eryn Cangi, Xian Chen, Nicholas Christiny, Christopher Clark, Michelle Collins, Johan Comparat, Neil Cook, Darren Croton, Isak Delberth Davids, Éric Depagne, John Donor, Leonardo A. dos Santos, Stephanie Douglas, Alan Du, Meredith Durbin, Dawn Erb, Daniel Faes, J. G. Fernández-Trincado, Anthony Foley, Sotiria Fotopoulou, Søren Frimann, Peter Frinchaboy, Rafael Garcia-Dias, Artur Gawryszczak, Elizabeth George, Sebastian Gonzalez, Karl Gordon, Nicholas Gorgone, Catherine Gosmeyer, Katie Grasha, Perry Greenfield, Rebekka Grellmann, James Guillochon, Mark Gurwell, Marcel Haas, Alex Hagen, Daryl Haggard, Tim Haines, Patrick Hall, Wojciech Hellwing, Edmund Christian Herenz, Samuel Hinton, Renee Hlozek, John Hoffman, Derek Holman, Benne Willem Holwerda, Anthony Horton, Cameron Hummels, Daniel Jacobs, Jens Juel Jensen, David Jones, Arna Karick, Luke Kelley, Matthew Kenworthy, Ben Kitchener, Dominik Klaes, Saul Kohn, Piotr Konorski, Coleman Krawczyk, Kyler Kuehn, Teet Kuutma, Michael T. Lam, Richard Lane, Jochen Liske, Diego Lopez-Camara, Katherine Mack, Sam Mangham, Qingqing Mao, David J. E. Marsh, Cecilia Mateu, Loïc Maurin, James McCormac, Ivelina Momcheva, Hektor Monteiro, Michael Mueller, Roberto Munoz, Rohan Naidu, Nicholas Nelson, Christian Nitschelm, Chris North, Juan Nunez-Iglesias, Sara Ogaz, Russell Owen, John Parejko, Vera Patrício, Joshua Pepper, Marshall Perrin, Timothy Pickering, Jennifer Piscionere, Richard Pogge, Radek Poleski, Alkistis Pourtsidou, Adrian M. Price-Whelan, Meredith L. Rawls, Shaun Read, Glen Rees, Hanno Rein, Thomas Rice, Signe Riemer-Sørensen, Naum Rusomarov, Sebastian F. Sanchez, Miguel Santander-García, Gal Sarid, William Schoenell, Aleks Scholz, Robert L. Schuhmann, William Schuster, Peter Scicluna, Marja Seidel, Lijing Shao, Pranav Sharma, Aleksandar Shulevski, David Shupe, Cristóbal Sifón, Brooke Simmons, Manodeep Sinha, Ian Skillen, Bjoern Soergel, Thomas Spriggs, Sundar Srinivasan, Abigail Stevens, Ole Streicher, Eric Suchyta, Joshua Tan, O. Grace Telford, Romain Thomas, Chiara Tonini, Grant Tremblay, Sarah Tuttle, Tanya Urrutia, Sam Vaughan, Miguel Verdugo, Alexander Wagner, Josh Walawender, Andrew Wetzel, Kyle Willett, Peter K. G. Williams, Guang Yang, Guangtun Zhu, Andrea Zonca
    Oct. 11, 2016 physics.soc-ph, astro-ph.IM
    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.