Cosmological Simulations with TreeSPH

We describe numerical methods for incorporating gas dynamics into cosmological simulations and present illustrative applications to the cold dark matter (CDM) scenario. Our evolution code, a version of TreeSPH (Hernquist \& Katz 1989) generalized to handle comoving coordinates and periodic boundary conditions, combines smoothed--particle hydrodynamics (SPH) with the hierarchical tree method for computing gravitational forces. The Lagrangian hydrodynamics approach and individual time steps for gas particles give the algorithm a large dynamic range, which is essential for studies of galaxy formation in a cosmological context. The code incorporates radiative cooling for an optically thin, primordial composition gas in ionization equilibrium with a user-specified ultraviolet background. We adopt a phenomenological prescription for star formation that gradually turns cold, dense, Jeans-unstable gas into collisionless stars, returning supernova feedback energy to the surrounding medium. In CDM simulations, some of the baryons that fall into dark matter potential wells dissipate their acquired thermal energy and condense into clumps with roughly galactic masses. The resulting galaxy population is insensitive to assumptions about star formation; we obtain similar baryonic mass functions and galaxy correlation functions from simulations with star formation and from simulations without star formation in which we identify galaxies directly from the cold, dense gas.Comment: compressed postscript, 38 pages including 6 out of 7 embedded figures. Submitted to ApJ Supplements. Version with all 7 figures available from ftp://bessel.mps.ohio-state.edu/pub/dhw/Preprint

Photoionization, Numerical Resolution, and Galaxy Formation

Using cosmological simulations that incorporate gas dynamics and gravitational forces, we investigate the influence of photoionization by a UV radiation background on the formation of galaxies. In our highest resolution simulations, we find that photoionization has essentially no effect on the baryonic mass function of galaxies at $z=2$, down to our resolution limit of 5e9 M_\sun. We do, however, find a strong interplay between the mass resolution of a simulation and the microphysics included in the computation of heating and cooling rates. At low resolution, a photoionizing background can appear to suppress the formation of even relatively massive galaxies. However, when the same initial conditions are evolved with a factor of eight better mass resolution, this effect disappears. Our results demonstrate the need for care in interpreting the results of cosmological simulations that incorporate hydrodynamics and radiation physics. For example, we conclude that a simulation with limited resolution may yield more realistic results if it ignores some relevant physical processes, such as photoionization. At higher resolution, the simulated population of massive galaxies is insensitive to the treatment of photoionization and star formation, but it does depend significantly on the amplitude of the initial density fluctuations. By $z=2$, an $\Omega=1$ cold dark matter model normalized to produce the observed masses of present-day clusters has already formed galaxies with baryon masses exceeding 1e11 M_\sun.Comment: 25 pages, w/ embedded figures. Submitted to ApJ. Also available at http://www-astronomy.mps.ohio-state.edu/~dhw/Docs/preprints.htm

Simulating Cosmic Structure Formation

We describe cosmological simulation techniques and their application to studies of cosmic structure formation, with particular attention to recent hydrodynamic simulations of structure in the high redshift universe. Collisionless N-body simulations with Gaussian initial conditions produce a pattern of sheets, filaments, tunnels, and voids that resembles the observed large scale galaxy distribution. Simulations that incorporate gas dynamics and dissipation form dense clumps of cold gas with sizes and masses similar to the luminous parts of galaxies. Models based on inflation and cold dark matter predict a healthy population of high redshift galaxies, including systems with star formation rates of 20 M_{\sun}/year at z=6. At z~3, most of the baryons in these models reside in the low density intergalactic medium, which produces fluctuating Lyman-alpha absorption in the spectra of background quasars. The physical description of this ``Lyman-alpha forest'' is particularly simple if the absorption spectrum is viewed as a 1-dimensional map of a continuous medium instead of a collection of lines. The combination of superb observational data and robust numerical predictions makes the Lyman-alpha forest a promising tool for testing cosmological models.Comment: Latex w/ paspconf.sty, 25 pages, 8 ps figs. To appear in Origins, eds. J. M. Shull, C. E. Woodward, & H. Thronson (ASP Conference Series

Baryon Dynamics, Dark Matter Substructure, and Galaxies

By comparing a collisionless cosmological N-body simulation (DM) to an SPH simulation with the same initial conditions, we investigate the correspondence between the dark matter subhalos produced by collisionless dynamics and the galaxies produced by dissipative gas dynamics in a dark matter background. When galaxies in the SPH simulation become satellites in larger groups, they retain local dark matter concentrations (SPH subhalos) whose mass is typically five times their baryonic mass. The more massive subhalos of the SPH simulation have corresponding subhalos of similar mass and position in the DM simulation; at lower masses, there is fairly good correspondence, but some DM subhalos are in different spatial positions and some suffer tidal stripping or disruption. The halo occupation statistics of DM subhalos -- the mean number of subhalos, pairs, and triples as a function of host halo mass -- are very similar to those of SPH subhalos and SPH galaxies. Gravity of the dissipative baryon component amplifies the density contrast of subhalos in the SPH simulation, making them more resistant to tidal disruption. Relative to SPH galaxies and SPH subhalos, the DM subhalo population is depleted in the densest regions of the most massive halos. The good agreement of halo occupation statistics between the DM subhalo and SPH galaxy populations leads to good agreement of their two-point correlation functions and higher order moments on large scales. The depletion of DM subhalos in dense regions depresses their clustering at R<1 Mpc/h. In these simulations, the "conversation" between dark matter and baryons is mostly one-way, with dark matter dynamics telling galaxies where to form and how to cluster, but the "back talk" of the baryons influences small scale clustering by enhancing the survival of substructure in the densest environments.Comment: 32 pages including 16 figs. Submitted to ApJ. PDF file with higher quality versions of Figs 2 and 3 available at http://www.astronomy.ohio-state.edu/~dhw/Preprints/subhalo.pd

How Do Galaxies Get Their Gas?

We examine the temperature history of gas accreted by forming galaxies in SPH simulations. About half the gas shock heats to roughly the virial temperature of the galaxy potential well before cooling, condensing, and forming stars, but the other half radiates its acquired gravitational energy at much lower temperatures, typically T<10^5 K, and the histogram of maximum gas temperatures is clearly bimodal. The "cold mode" of gas accretion dominates for low mass galaxies (M_baryon < 10^{10.3}Msun or M_halo < 10^{11.4}Msun), while the conventional "hot mode" dominates the growth of high mass systems. Cold accretion is often directed along filaments, allowing galaxies to efficiently draw gas from large distances, while hot accretion is quasi-spherical. The galaxy and halo mass dependence leads to redshift and environment dependence of cold and hot accretion rates, with cold mode dominating at high redshift and in low density regions today, and hot mode dominating in group and cluster environments at low redshift. Star formation rates closely track accretion rates, and we discuss the physics behind the observed environment and redshift dependence of galactic scale star formation. If we allowed hot accretion to be suppressed by conduction or AGN feedback, then the simulation predictions would change in interesting ways, perhaps resolving conflicts with the colors of ellipticals and the cutoff of the galaxy luminosity function. The transition between cold and hot accretion at M_h ~ 10^{11.4}Msun is similar to that found by Birnboim & Dekel (2003) using 1-d simulations and analytic arguments. The corresponding baryonic mass is tantalizingly close to the scale at which Kauffmann et al. (2003) find a marked shift in galaxy properties. We speculate on connections between these theoretical and observational transitions.Comment: 1 figure added, Appendix discussing SAMs added, some text changes. Matches the version accepted by MNRAS. 31 pages (MNRAS style), 21 figures,For high resolution version of the paper (highly recommended) follow http://www.astro.umass.edu/~keres/paper/ms2.ps.g

The Lyman-alpha Forest as a Cosmological Tool

We review recent developments in the theory of the Lyman-alpha forest and their implications for the role of the forest as a test of cosmological models. Simulations predict a relatively tight correlation between the local Lya optical depth and the local gas or dark matter density. Statistical properties of the transmitted flux can constrain the amplitude and shape of the matter power spectrum at high redshift, test the assumption of Gaussian initial conditions, and probe the evolution of dark energy by measuring the Hubble parameter H(z). Simulations predict increased Lya absorption in the vicinity of galaxies, but observations show a Lya deficit within Delta_r ~ 0.5 Mpc/h (comoving). We investigate idealized models of "winds" and find that they must eliminate neutral hydrogen out to comoving radii ~1.5 Mpc/h to marginally explain the data. Winds of this magnitude suppress the flux power spectrum by \~0.1 dex but have little effect on the distribution function or threshold crossing frequency. In light of the stringent demands on winds, we consider the alternative possibility that extended Lya emission from target galaxies replaces absorbed flux, but we conclude that this explanation is unlikely. Taking full advantage of the data coming from large telescopes and from the Sloan Digital Sky Survey will require more complete understanding of the galaxy proximity effect, careful attention to continuum determination, and more accurate numerical predictions, with the goal of reaching 5-10% precision on key cosmological quantities.Comment: 13 pages, to appear in "The Emergence of Cosmic Structure," Proceedings of the 13th Annual Astrophysics Conference in Maryland, eds. S. Holt and C. Reynolds, AIP Press, 200

The Galaxy Proximity Effect in the Lyman-alpha Forest

Hydrodynamic cosmological simulations predict that the average opacity of the Ly-alpha forest should increase in the neighborhood of galaxies because galaxies form in dense environments. Recent observations (Adelberger et al. 2002) confirm this expectation at large scales, but they show a decrease of absorption at comoving separations Delta_r <~ 1 Mpc/h. We show that this discrepancy is statistically significant, especially for the innermost data point at Delta_r <= 0.5 Mpc/h, even though this data point rests on three galaxy-quasar pairs. Galaxy redshift errors of the expected magnitude are insufficient to resolve the conflict. Peculiar velocities allow gas at comoving distances >~ 1 Mpc/h to produce saturated absorption at the galaxy redshift, putting stringent requirements on any ``feedback'' solution. Local photoionization is insufficient, even if we allow for recurrent AGN activity that keeps the neutral hydrogen fraction below its equilibrium value. A simple ``wind'' model that eliminates all neutral hydrogen in spheres around the observed galaxies can marginally explain the data, but only if the winds extend to comoving radii ~1.5 Mpc/h.Comment: 4 pages, 1 figure; To appear in proceedings of the 13th Annual Astrophysics Conference in College Park, Maryland, The Emergence of Cosmic Structure, eds. S.Holt and C. Reynolds, (AIP

Trends in U. S. Wage Inequality: Re-Assessing the Revisionists

A large literature documents a substantial rise in U. S. wage inequality and educational wage differentials over the past several decades and finds that these trends can be primarily accounted for by shifts in the supply of and demand for skills reinforced by the erosion of labor market institutions affecting the wages of low- and middle-wage workers. Drawing on an additional decade of data, a number of recent contributions reject this consensus to conclude that (1) the rise in wage inequality was an “episodic” event of the first-half of the 1980s rather than a "secular” phenomenon, (2) this rise was largely caused by a falling minimum wage rather than by supply and demand factors; and (3) rising residual wage inequality since the mid-1980s is explained by confounding effects of labor force composition rather than true increases in inequality within detailed demographic groups. We reexamine these claims using detailed data from the Current Population Survey and find only limited support. Although the growth of overall inequality in the U. S. slowed in the 1990s, upper tail inequality rose almost as rapidly during the 1990s as during the 1980s. A decomposition applied to the CPS data reveals large and persistent rise in within-group earnings inequality over the past several decades, controlling for changes in labor force composition. While changes in the minimum wage can potentially account for much of the movement in lower tail earnings inequality, strong time series correlations of the evolution of the real minimum wage and upper tail wage inequality raise questions concerning the causal interpretation of such relationships. We also find that changes in the college/high school wage premium appear to be well captured by standard models emphasizing rapid secular growth in the relative demand for skills and fluctuations in the rate of growth of the relative supply of college workers – though these models do not accurately predict the slowdown in the growth of the college/high-school gap during the 1990s. We conclude that these patterns are not adequately explained by either a ‘unicausal’ skill-biased technical change explanation or a revisionist hypothesis focused primarily on minimum wages and mechanical labor force compositional effects. We speculate that these puzzles can be partially reconciled by a modified version of the skill-biased technical change hypothesis that generates a polarization of skill demands.

Rising Wage Inequality: The Role of Composition and Prices

During the early 1980s, earnings inequality in the U.S. labor market rose relatively uniformly throughout the wage distribution. But this uniformity gave way to a significant divergence starting in 1987, with upper-tail (90/50) inequality rising steadily and lower tail (50/10) inequality either flattening or compressing for the next 16 years (1987 to 2003). This paper applies and extends a quantile decomposition technique proposed by Machado and Mata (2005) to evaluate the role of changing labor force composition (in terms of education and experience) and changing labor market prices to the expansion and subsequent divergence of upper- and lower-tail inequality over the last three decades We show that the extended Machado-Mata quantile decomposition corrects shortcomings of the original Juhn-Murphy-Pierce (1993) full distribution accounting method and nests the kernel reweighting approach proposed by DiNardo, Fortin and Lemieux (1996). Our analysis reveals that shifts in labor force composition have positively impacted earnings inequality during the 1990s. But these compositional shifts have primarily operated on the lower half of the earnings distribution by muting a contemporaneous, countervailing lower-tail price compression. The steady rise of upper tail inequality since the late 1970s appears almost entirely explained by ongoing between-group price changes (particularly increasing wage differentials by education) and residual price changes.

The Polarization of the U.S. Labor Market

This paper analyzes a marked change in the evolution of the U.S. wage structure over the past fifteen years: divergent trends in upper-tail (90/50) and lower-tail (50/10) wage inequality. We document that wage inequality in the top half of distribution has displayed an unchecked and rather smooth secular rise for the last 25 years (since 1980). Wage inequality in the bottom half of the distribution also grew rapidly from 1979 to 1987, but it has ceased growing (and for some measures actually narrowed) since the late 1980s. Furthermore we find that occupational employment growth shifted from monotonically increasing in wages (education) in the 1980s to a pattern of more rapid growth in jobs at the top and bottom relative to the middles of the wage (education) distribution in the 1990s. We characterize these patterns as the %u201Cpolarization%u201D of the U.S. labor market, with employment polarizing into high-wage and low-wage jobs at the expense of middle-wage work. We show how a model of computerization in which computers most strongly complement the non-routine (abstract) cognitive tasks of high-wage jobs, directly substitute for the routine tasks found in many traditional middle-wage jobs, and may have little direct impact on non-routine manual tasks in relatively low-wage jobs can help explain the observed polarization of the U.S. labor market.