Search

Galaxy cluster mass halos ("clusters") in a dark matter simulation are matched to nodes in several different cosmic webs found using the Disperse cosmic web finder. The webs have different simulation smoothings and Disperse parameter choices; for each, 4 methods are considered for matching Disperse nodes to clusters. For most of the webs, Disperse nodes outnumber clusters, but not every cluster has a Disperse node match (and sometimes $>1$ cluster matches to the same Disperse node). The clusters frequently lacking a matching Disperse node have a different distribution of local shear trends and perhaps merger histories. It might be interesting to see in what other ways, e.g., observational properties, these clusters differ. For the webs with smoothing $\leq$ 2.5 $Mpc/h$, and all but the most restrictive matching criterion, $\sim$3/4 of the clusters always have a Disperse node counterpart. The nearest cluster to a given Disperse node and vice versa, within twice the smoothing length, obey a cluster mass-Disperse node density relation. Cluster pairs where both clusters match Disperse nodes can also be assigned the filaments between those nodes, but as the web and matching methods are varied, most such filaments do not remain. There is an enhancement of subhalo counts and halo mass between cluster pairs, averaging over cluster pairs assigned Disperse filaments increases the enhancement. The approach here also lends itself to comparing nodes across many cosmic web constructions, using the fixed underlying cluster distribution to make a correspondence., Comment: 15 pages, 15 Figures. Revised version incorporating many helpful suggestions from the referee; more information about nodes matched or not to clusters, cluster pair profile quantities changed

One emerging application of machine learning methods is the inference of galaxy cluster masses. In this note, machine learning is used to directly combine five simulated multiwavelength measurements in order to find cluster masses. This is in contrast to finding mass estimates for each observable, normally by using a scaling relation, and then combining these scaling law based mass estimates using a likelihood. We also illustrate how the contributions of each observable to the accuracy of the resulting mass measurement can be compared via model-agnostic Importance Permutation values. Thirdly, as machine learning relies upon the accuracy of the training set in capturing observables, their correlations, and the observational selection function, and as the machine learning training set originates from simulations, two tests of whether a simulation's correlations are consistent with observations are suggested and explored as well., Comment: 12 pages, final version to appear in MNRAS. Expanded explanations, and added examples

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space., Comment: Science white paper submitted to the Astro2020 Decadal Survey

Upcoming imaging surveys, such as LSST, will provide an unprecedented view of the Universe, but with limited resolution along the line-of-sight. Common ways to increase resolution in the third dimension, and reduce misclassifications, include observing a wider wavelength range and/or combining the broad-band imaging with higher spectral resolution data. The challenge with these approaches is matching the depth of these ancillary data with the original imaging survey. However, while a full 3D map is required for some science, there are many situations where only the statistical distribution of objects (dN/dz) in the line-of-sight direction is needed. In such situations, there is no need to measure the fluxes of individual objects in all of the surveys. Rather a stacking procedure can be used to perform an `ensemble photo-z'. We show how a shallow, higher spectral resolution survey can be used to measure dN/dz for stacks of galaxies which coincide in a deeper, lower resolution survey. The galaxies in the deeper survey do not even need to appear individually in the shallow survey. We give a toy model example to illustrate tradeoffs and considerations for applying this method. This approach will allow deep imaging surveys to leverage the high resolution of spectroscopic and narrow/medium band surveys underway, even when the latter do not have the same reach to high redshift., Comment: 5 pages, 3 figures

Galaxies evolve via a complex interaction of numerous different physical processes, scales and components. In spite of this, overall trends often appear. Simplified models for galaxy histories can be used to search for and capture such emergent trends, and thus to interpret and compare results of galaxy formation models to each other and to nature. Here, two approximations are applied to galaxy integrated star formation rate histories, drawn from a semi-analytic model grafted onto a dark matter simulation. Both a lognormal functional form and principal component analysis (PCA) approximate the integrated star formation rate histories fairly well. Machine learning, based upon simplified galaxy halo histories, is somewhat successful at recovering both fits. The fits to the histories give fixed time star formation rates which have notable scatter from their true fixed time rates at final time, especially for quiescent and "green valley" galaxies, and more so for the PCA fit. For classifying galaxies into subfamilies sharing similar integrated histories, both approximations are better than using final stellar mass or specific star formation rate. Several subsamples from the simulation illustrate how these simple parameterizations can provide points of contact for comparisons between different galaxy formation samples, or more generally, models. As a side result, the halo masses of simulated galaxies with early peak star formation rate (according to the lognormal fit) lie on one of two tracks. The small fraction of galaxies with a lower halo mass at peak star formation rate appear to stall in their halo growth, even though they are central in their host halos., Comment: Final version, to appear in the MNRAS. Helpful suggestions and tests from referee now included. Reference to PCA work on star formation rate histories now correctly attributed to Sparre et al (2015)

Dark matter simulations can serve as a basis for creating galaxy histories via the galaxy-dark matter connection. Here, one such model by Becker (2015) is implemented with several variations on three different dark matter simulations. Stellar mass and star formation rates are assigned to all simulation subhalos at all times, using subhalo mass gain to determine stellar mass gain. The observational properties of the resulting galaxy distributions are compared to each other and observations for a range of redshifts from 0-2. Although many of the galaxy distributions seem reasonable, there are noticeable differences as simulations, subhalo mass gain definitions, or subhalo mass definitions are altered, suggesting that the model should change as these properties are varied. Agreement with observations may improve by including redshift dependence in the added-by-hand random contribution to star formation rate. There appears to be an excess of faint quiescent galaxies as well (perhaps due in part to differing definitions of quiescence). The ensemble of galaxy formation histories for these models tend to have more scatter around their average histories (for a fixed final stellar mass) than the two more predictive and elaborate semi-analytic models of Guo et al (2013) and Henriques et al (2015), and require more basis fluctuations (using PCA) to capture 90 percent of the scatter around their average histories. The codes to plot model predictions (in some cases alongside observational data) are publicly available to test other mock catalogues at https://github.com/jdcphysics/validation . Information on how to use these codes is in the appendix., Comment: Paper also serves as documentation for public code validating mock catalogues which have M*, star formation rates, stellar mass to halo mass information (redshifts 0 to ~ 4, most data sets for 0 to 2)

Acoustic waves traveling through the early Universe imprint a characteristic scale in the clustering of galaxies, QSOs and inter-galactic gas. This scale can be used as a standard ruler to map the expansion history of the Universe, a technique known as Baryon Acoustic Oscillations (BAO). BAO offer a high-precision, low-systematics means of constraining our cosmological model. The statistical power of BAO measurements can be improved if the `smearing' of the acoustic feature by non-linear structure formation is undone in a process known as reconstruction. In this paper we use low-order Lagrangian perturbation theory to study the ability of $21\,$cm experiments to perform reconstruction and how augmenting these surveys with galaxy redshift surveys at relatively low number densities can improve performance. We find that the critical number density which must be achieved in order to benefit $21\,$cm surveys is set by the linear theory power spectrum near its peak, and corresponds to densities achievable by upcoming surveys of emission line galaxies such as eBOSS and DESI. As part of this work we analyze reconstruction within the framework of Lagrangian perturbation theory with local Lagrangian bias, redshift-space distortions, ${\bf k}$-dependent noise and anisotropic filtering schemes., Comment: 10 pages, final version to appear in MNRAS, helpful suggestions from referee and others included

Halos, filaments, sheets and voids in the cosmic web can be defined in terms of the eigenvalues of the smoothed shear tensor and a threshold $\lambda_{\rm th}$. Using analytic methods, we construct mean maps centered on these types of structures for Gaussian random fields corresponding to cosmological initial conditions. Each map also requires a choice of shear at the origin; we consider three possibilities. We find characteristic sizes, shapes and other properties of the central objects in these mean maps and explore how these properties change with varying the threshold and smoothing scale, i.e. varying the separation of the cosmic web into different kinds of components. The mean maps become increasingly complex as the threshold $\lambda_{\rm th}$ decreases to zero. We also describe scatter around these mean maps, subtleties which can arise in their construction, and some comparisons between halos in the maps and collapsed halos at final times., Comment: 16 pages, 8 figures. Added explanations and examples in response to comments from referee and others. Final version to appear in MNRAS

Cosmological galaxy formation simulations can now produce rich and diverse ensembles of galaxy histories. These simulated galaxy histories, taken all together, provide an answer to the question "how do galaxies form?" for the models used to construct them. We characterize such galaxy history ensembles both to understand their properties and to identify points of comparison for histories within a given galaxy formation model or between different galaxy formation models and simulations. We focus primarily on stellar mass histories of galaxies with the same final stellar mass, for six final stellar mass values and for three different simulated galaxy formation models (a semi-analytic model built upon the dark matter Millennium simulation and two models from the hydrodynamical OverWhelmingly Large Simulations project). Using principal component analysis (PCA) to classify scatter around the average stellar mass history, we find that one fluctuation dominates for all sets of histories we consider, although its shape and contribution can vary between samples. We correlate the PCA characterization with several $z=0$ galaxy properties (to connect with survey observables) and also compare it to some other galaxy history properties. We then explore separating galaxy stellar mass histories into classes, using the largest PCA contribution, k-means clustering, and simple Gaussian mixture models. For three component models, these different methods often gave similar results. These history classification methods provide a succinct and often quick way to characterize changes in the full ensemble of histories of a simulated population as physical assumptions are varied, to compare histories of different simulated populations to each other, and to assess the relation of simulated histories to fixed time observations., Comment: 17 pages, 16 figures. Final version to appear in MNRAS; expanded introduction and comparisons of additional measurements

One of the most striking properties of galaxies is the bimodality in their star-formation rates. A major puzzle is why any given galaxy is star-forming or quiescent, and a wide range of physical mechanisms have been proposed as solutions. We consider how observations, such as might be available in upcoming large galaxy surveys, might distinguish different galaxy quenching scenarios. To do this, we combine an N-body simulation and multiple prescriptions from the literature to create several quiescent galaxy mock catalogues. Each prescription uses a different set of galaxy properties (such as history, environment, centrality) to assign individual simulation galaxies as quiescent. We find how and how much the resulting quiescent galaxy distributions differ from each other, both intrinsically and observationally. In addition to tracing observational consequences of different quenching mechanisms, our results indicate which sorts of quenching models might be most readily disentangled by upcoming observations and which combinations of observational quantities might provide the most discriminatory power. Our observational measures are auto, cross, and marked correlation functions, pro- jected density distributions, and group multiplicity functions, which rely upon galaxy positions, stellar masses and of course quiescence. Although degeneracies between models are present for individual observations, using multiple observations in concert allows us to distinguish between all ten models we consider. In addition to identifying intrinsic and observational consequences of quiescence prescriptions and testing these quiescence models against each other and observations, these methods can also be used to validate colors (or other history and environment dependent properties) in simulated mock catalogues., Comment: 20 pages, 14 figures. Version to appear in MNRAS, incorporating helpful suggestions from referee and others

The challenge of obtaining galaxy cluster masses is increasingly being addressed by multiwavelength measurements. As scatters in measured cluster masses are often sourced by properties of or around the clusters themselves, correlations between mass scatters are frequent and can be significant, with consequences for errors on mass estimates obtained both directly and via stacking. Using a high resolution 250 Mpc/h side N-body simulation, combined with proxies for observational cluster mass measurements, we obtain mass scatter correlations and covariances for 243 individual clusters along ~96 lines of sight each, both separately and together. Many of these scatters are quite large and highly correlated. We use principal component analysis (PCA) to characterize scatter trends and variations between clusters. PCA identifies combinations of scatters, or variations more generally, which are uncorrelated or non-covariant. The PCA combination of mass measurement techniques which dominates the mass scatter is similar for many clusters, and this combination is often present in a large amount when viewing the cluster along its long axis. We also correlate cluster mass scatter, environmental and intrinsic properties, and use PCA to find shared trends between these. For example, if the average measured richness, velocity dispersion and Compton decrement mass for a cluster along many lines of sight are high relative to its true mass, in our simulation the cluster's mass measurement scatters around this average are also high, its sphericity is high, and its triaxiality is low. Our analysis is based upon estimated mass distributions for fixed true mass. Extensions to observational data would require further calibration from numerical simulations, tuned to specific observational survey selection functions and systematics., Comment: 18 pages, 12 figures, final version to appear in MNRAS, helpful changes from referee and others incorporated

Weak gravitational lensing is a valuable probe of galaxy formation and cosmology. Here we quantify the effects of using photometric redshifts (photo-z) in galaxy-galaxy lensing, for both sources and lenses, both for the immediate goal of using galaxies with photo-z as lenses in the Sloan Digital Sky Survey (SDSS) and as a demonstration of methodology for large, upcoming weak lensing surveys that will by necessity be dominated by lens samples with photo-z. We calculate the bias in the lensing mass calibration as well as consequences for absolute magnitude (i.e., k-corrections) and stellar mass estimates, for a large sample of SDSS Data Release 8 (DR8) galaxies. The redshifts are obtained with the template based photo-z code ZEBRA on the SDSS DR8 ugriz photometry. We assemble and characterise the calibration samples (~9k spectroscopic redshifts from four surveys) to obtain photometric redshift errors and lensing biases corresponding to our full SDSS DR8 lens and source catalogues. Our tests of the calibration sample also highlight the impact of observing conditions in the imaging survey when the spectroscopic calibration covers a small fraction of its footprint; atypical imaging conditions in calibration fields can lead to incorrect conclusions regarding the photo-z of the full survey. For the SDSS DR8 catalogue, we find sigma_z/(1+z)=0.096 and 0.113 for the lens and source catalogues, with flux limits of r=21 and r=21.8, respectively. We also explore the systematic uncertainty in the lensing signal calibration when using source photo-z, and both lens and source photo-z; given the size of existing training samples, we can constrain the lensing signal calibration (and therefore the normalization of the surface mass density) to within 2 and 4 per cent, respectively. [ABRIDGED], Comment: 26 pages, 28 figures, minor changes, accepted for publication in MNRAS

Galaxies which fall into clusters as part of the same infall halo can retain correlations due to their shared origin. N-body simulations are used to study properties of such galaxy subgroups within clusters, including their richnesses and prevalence. The sizes, densities and velocity dispersions of all subgroups with >= 8 galaxies are found and compared to those of the host clusters. The largest galaxy subgroup provides a preferred direction in the cluster and is compared to other preferred directions in the cluster. Scatter in cluster mass measurements (via five observables), along ~ 96 lines of sight, is compared to the relation of the line of sight to this preferred direction: scatter in cluster velocity dispersion measurements show the strongest correlation. The Dressler-Shectman test (Dressler & Shectman 1988), is applied to these clusters, to see whether the substructure it identifies is related to these subgroups. The results for any specific line of sight seem noisy; however, clusters with large subgroups tend to have a higher fraction of lines of sight where the test detects substructure., Comment: 12 pages, final version for publication with helpful comments from referee and others included

We construct a filament catalogue using an extension of the halo based filament finder of Zhang et al.(2009), in a 250 Mpc/h side N-body simulation, and study the properties of filaments ending upon or surrounding galaxy clusters (within 10 Mpc/h). In this region, the majority of filamentary mass, halo mass, and galaxy richness centered upon the cluster tends to lie in sheets, which are not always coincident. Fixing a sheet width of 3 Mpc/h for definiteness, we find the sheet orientations and (connected) filamentary mass, halo mass and richness fractions relative to the surrounding sphere. Filaments usually have one or more endpoints outside the sheet determined by filament or halo mass or richness, with at least one having a large probability to be aligned with the perpendicular of the plane. Scatter in mock cluster mass measurements, for several observables, is often correlated with the observational direction relative to these local sheets, most often for richness and weak lensing, somewhat less for Compton decrement, and least often for velocity dispersions. The long axis of the cluster also tends to lie in the sheet andits orientiation relative to line of sight also correlates with mass scatter., Comment: References and explanations added as suggested, revised results for cluster inertia tensor (some correlations stronger than reported previously)

We use a high-resolution N-body simulation to study how the influence of large-scale structure in and around clusters causes correlated signals in different physical probes and discuss some implications this has for multi-physics probes of clusters. We pay particular attention to velocity dispersions, matching galaxies to subhalos which are explicitly tracked in the simulation. We find that not only do halos persist as subhalos when they fall into a larger host, groups of subhalos retain their identity for long periods within larger host halos. The highly anisotropic nature of infall into massive clusters, and their triaxiality, translates into an anisotropic velocity ellipsoid: line-of-sight galaxy velocity dispersions for any individual halo show large variance depending on viewing angle. The orientation of the velocity ellipsoid is correlated with the large-scale structure, and thus velocity outliers correlate with outliers caused by projection in other probes. We quantify this orientation uncertainty and give illustrative examples. Such a large variance suggests that velocity dispersion estimators will work better in an ensemble sense than for any individual cluster, which may inform strategies for obtaining redshifts of cluster members. We similarly find that the ability of substructure indicators to find kinematic substructures is highly viewing angle dependent. While groups of subhalos which merge with a larger host halo can retain their identity for many Gyr, they are only sporadically picked up by substructure indicators. (Abridged), Comment: 19 pages, 14 figures, published version (revised according to suggestions from referee and others)

Recently Eisenstein and collaborators introduced a method to `reconstruct' the linear power spectrum from a non-linearly evolved galaxy distribution in order to improve precision in measurements of baryon acoustic oscillations. We reformulate this method within the Lagrangian picture of structure formation, to better understand what such a method does, and what the resulting power spectra are. We show that reconstruction does not reproduce the linear density field, at second order. We however show that it does reduce the damping of the oscillations due to non-linear structure formation, explaining the improvements seen in simulations. Our results suggest that the reconstructed power spectrum is potentially better modeled as the sum of three different power spectra, each dominating over different wavelength ranges and with different non-linear damping terms. Finally, we also show that reconstruction reduces the mode-coupling term in the power spectrum, explaining why mis-calibrations of the acoustic scale are reduced when one considers the reconstructed power spectrum., Comment: 6 pages, 3 figs, PRD submitted. (v2) typo fixed in Eq. 5. (v3) Matches version accepted for publication

High-resolution simulations of cosmological structure formation indicate that dark matter substructure in dense environments, like groups and clusters, may survive for a long time. These dark matter subhalos are the likely hosts of galaxies. We examine the small-scale spatial clustering of subhalo major mergers at high redshift using high-resolution N-body simulations of cosmological volumes. Recently merged, massive subhalos exhibit enhanced clustering on scales ~100-300 kpc/h, relative to all subhalos of the same infall mass, for a short time after a major merger (< 500 Myr). The small-scale clustering enhancement is smaller for lower mass subhalos, which also show a deficit on scales just beyond the excess. Halos hosting recent subhalo mergers tend to have more subhalos; for massive subhalos the excess is stronger and it tends to increase for the most massive host halos. The subhalo merger fraction is independent of halo mass for the scales we probe. In terms of satellite and central subhalos, the merger increase in small-scale clustering for massive subhalos arises from recently merged massive central subhalos having an enhanced satellite population. Our mergers are defined via their parent infall mass ratios. Subhalos experiencing major mass gains also exhibit a small-scale clustering enhancement, but these correspond to two-body interactions leading to two final subhalos, rather than subhalo coalescence., Comment: 10 pages, 5 figures, accepted in MNRAS, version matches published one

Galaxies are believed to be in one-to-one correspondence with simulated dark matter subhalos. We use high-resolution N-body simulations of cosmological volumes to calculate the statistical properties of subhalo (galaxy) major mergers at high redshift (z=0.6-5). We measure the evolution of the galaxy merger rate, finding that it is much shallower than the merger rate of dark matter host halos at z>2.5, but roughly parallels that of halos at z<1.6. We also track the detailed merger histories of individual galaxies and measure the likelihood of multiple mergers per halo or subhalo. We examine satellite merger statistics in detail: 15%-35% of all recently merged galaxies are satellites and satellites are twice as likely as centrals to have had a recent major merger. Finally, we show how the differing evolution of the merger rates of halos and galaxies leads to the evolution of the average satellite occupation per halo, noting that for a fixed halo mass, the satellite halo occupation peaks at z~2.5., Comment: 17 pages, 12 figures, accepted in MNRAS, version matches published one

We apply simple optical and SZ cluster finders to mock galaxy catalogues and SZ flux maps created from dark matter halos in a (1 Gpc/h)^3 dark matter simulation, at redshifts 0.5 and 0.9. At each redshift, the two catalogues are then combined to assess how well they can improve each other, and compared to several variants of catalogues made using SZ flux and galaxy information simultaneously. We use several different criteria to compare the catalogues, and illustrate some of the tradeoffs which arise in tuning the galaxy cluster finders with respect to these criteria. We detail many of the resulting improvements and issues which arise in comparing and combining these two types of data sets., Comment: 14 pages, added information thanks to helpful suggestions from referee

Recent measurements of high-redshift QSO clustering from the Sloan Digital Sky Survey indicate that QSOs at z~4 have a bias b~14. We find that this extremely high clustering amplitude, combined with the corresponding space density, constrains the dispersion in the L-Mhalo relation to be less than 50% at 99% confidence for the most conservative case of a 100% duty cycle. This upper limit to the intrinsic dispersion provides as strong a constraint as current upper limits to the intrinsic dispersion in the local M_BH-sigma relation and the ratio of bolometric to Eddington luminosity of luminous QSOs., Comment: 6 pages, 4 figures. Minor revisions to address questions from referee. References updated

We investigate halo mass selection properties of red-sequence cluster finders using galaxy populations of the Millennium Simulation (MS). A clear red sequence exists for MS galaxies in massive halos at redshifts z < 1, and we use this knowledge to inform a cluster-finding algorithm applied to 500 Mpc/h projections of the simulated volume. At low redshift (z=0.4), we find that 90% of the clusters found have galaxy membership dominated by a single, real-space halo, and that 10% are blended systems for which no single halo contributes a majority of a cluster's membership. At z=1, the fraction of blends increases to 22%, as weaker redshift evolution in observed color extends the comoving length probed by a fixed range of color. Other factors contributing to the increased blending at high-z include broadening of the red sequence and confusion from a larger number of intermediate mass halos hosting bright red galaxies of magnitude similar to those in higher mass halos. Our method produces catalogs of cluster candidates whose halo mass selection function, p(M|\Ngal,z), is characterized by a bimodal log-normal model with a dominant component that reproduces well the real-space distribution, and a redshift-dependent tail that is broader and displaced by a factor ~2 lower in mass. We discuss implications for X-ray properties of optically selected clusters and offer ideas for improving both mock catalogs and cluster-finding in future surveys., Comment: final version to appear in MNRAS. Appendix added on purity and completeness, small shift in red sequence due to correcting an error in finding it

We use a suite of high-resolution N-body simulations to study the properties, abundance and clustering of high mass halos at high redshift, including their mass assembly histories and mergers. We find that the analytic form which best fits the abundance of halos depends sensitively on the assumed definition of halo mass, with common definitions of halo mass differing by a factor of two for these low concentration, massive halos. A significant number of massive halos are undergoing rapid mass accretion, with major merger activity being common. We compare the mergers and mass accretion histories to the extended Press-Schechter formalism. We consider how major merger induced star formation or black hole accretion may change the distribution of photon production from collapsed halos, and hence reionization, using some simplified examples. In all of these, the photon distribution for a halo of a given mass acquires a large scatter. If rare, high mass halos contribute significantly to the photon production rates, the scatter in photon production rate can translate into additional scatter in the sizes of ionized bubbles., Comment: 9 pages, 9 figures, resubmitted to MNRAS, added some references and convergence tests to better constrain bias

The clustering properties of dark matter halos are a firm prediction of modern theories of structure formation. We use two large volume, high-resolution N-body simulations to study how the correlation function of massive dark matter halos depends upon their mass and formation history. We find that halos with the lowest concentrations are presently more clustered than those of higher concentration, the size of the effect increasing with halo mass; this agrees with trends found in studies of lower mass halos. The clustering dependence on other characterizations of the full mass accretion history appears weaker than the effect with concentration. Using the integrated correlation function, marked correlation functions, and a power-law fit to the correlation function, we find evidence that halos which have recently undergone a major merger or a large mass gain have slightly enhanced clustering relative to a randomly chosen population with the same mass distribution., Comment: 10 pages, 8 figures; text improved, references and one figure added; accepted for publication in ApJ

The recently introduced models of reionization bubbles based on extended Press-Schechter theory (Furlanetto, Zaldarriaga & Hernquist 2004) are generalized to include mergers of ionization sources. Sources with a recent major merger are taken to have enhanced photon production due to star formation, and accretion onto a central black hole if a black hole is present. This produces a scatter in the number of ionized photons corresponding to a halo of a given mass and a change in photon production over time for any given halo mass. Photon production histories, bubble distributions, and ionization histories are computed for several different parameter and recombination assumptions; the resulting distributions interpolate between previously calculated limiting cases., Comment: 44 pages, 11 figures, version to appear in MNRAS. Some discussion of case with WMAP parameters and expanded explanations

Galaxy cluster counts at high redshift, binned into spatial pixels and binned into ranges in an observable proxy for mass, contain a wealth of information on both the dark energy equation of state and the mass selection function required to extract it. The likelihood of the number counts follows a Poisson distribution whose mean fluctuates with the large-scale structure of the universe. We develop a joint likelihood method that accounts for these distributions. Maximization of the likelihood over a theoretical model that includes both the cosmology and the observable-mass relations allows for a joint extraction of dark energy and cluster structural parameters., Comment: 4 pages

A sample of hundreds of simulated galaxy clusters is used to study the statistical properties of galaxy cluster formation. Individual assembly histories are discussed, the degree of virialization is demonstrated and various commonly used formation times are measured and inter-compared. In addition, the fraction of clusters which have ``recently'' undergone a major merger or significant mass jump is calculated as a function of lookback time and interval. The fraction of three- and four-body mergers is also studied., Comment: 24 pages, final version to appear in Astroparticle Physics

Poisson distributed shot noise is normally considered in the Gaussian limit in cosmology. However, if the shot noise is large enough and the correlation function/power spectrum conspires, the Gaussian approximation mis-estimates the errors and their covariance significantly. The power spectrum, even for initially Gaussian densities,acquires cross correlations which can be large, while the change in the correlation function error matrix is diagonal except at zero separation. Two and three dimensional power law correlation function and power spectrum examples are given. These corrections appear to have a large effect when applied to galaxy clusters, e.g. for SZ selected galaxy clusters in 2 dimensions. This can increase the error estimates for cosmological parameter estimation and consequently affect survey strategies, as the corrections are minimized for surveys which are deep and narrow rather than wide and shallow. In addition, a rewriting of the error matrix for the power spectrum/correlation function is given which eliminates most of the Bessel function dependence (in two dimensions) and all of it (in three dimensions), which makes the calculation of the error matrix more tractable. This applies even when the shot noise is in the (usual) Gaussian limit., Comment: 22 pages, 4 figures, 3 equations corrected/figures updated, results unchanged

Large SZ selected galaxy cluster surveys are beginning imminently. We compare the dependence of the galaxy cluster angular power spectrum on cosmological parameters, different modeling assumptions and statistical observational errors. We quantify the degeneracies between theoretical assumptions such as the mass function and cosmological parameters such as $\sigma_8$. We also identify a rough scaling behavior of this angular power spectrum with $\sigma_8$ alone., Comment: Final version to appear in ApJ. Expanded the discussion on the degeneracy of sigma_8 and Omega_m

Surveys for gravitational lens systems have typically found a significantly larger fraction of lenses with four (or more) images than are predicted by standard ellipsoidal lens models (50% versus 25-30%). We show that including the effects of smaller satellite galaxies, with an abundance normalized by the observations, significantly increases the expected number of systems with more than two images and largely explains the discrepancy. The effect is dominated by satellites with ~20% the luminosity of the primary lens, in rough agreement with the typical luminosities of the observed satellites. We find that the lens systems with satellites cannot, however, be dropped from estimates of the cosmological model based on gravitational lens statistics without significantly biasing the results., Comment: 23 pages, 7 figures, more discussion of sis vs sie and inclusion of uncorrelated contributions

This Resource Letter provides a guide to the literature on the theory of anisotropies in the cosmic microwave background. Journal articles, web pages, and books are cited for the following topics: discovery, cosmological origin, early work, recombination, general CMB anisotropy references, primary CMB anisotropies (numerical, analytical work), secondary effects, Sunyaev-Zel'dovich effect(s), lensing, reionization, polarization, gravity waves, defects, topology, origin of fluctuations, development of fluctuations, inflation and other ties to particle physics, parameter estimation, recent constraints, web resources, foregrounds, observations and observational issues, and gaussianity., Comment: AJP/AAPT Bibliographic Resource letter published Feb. 2002, 24 pages (9 of text), 1 figure

We use large volume, high resolution, N-body simulations of 3 different $\Lambda$CDM models, with different clustering strengths, to generate dark matter halo merging histories. Over the reliable range of halo masses, roughly galaxy groups to rich clusters of galaxies, we quantify the number density of major mergers for two different time intervals and compare with analytic predictions based on the extended Press-Schechter (1974) theory., Comment: 24 pages, 18 figures. Comments welcome. Final version to appear MNRAS, helpful suggestions from referee and others incorporated

We consider a wide range of parametric mass models for B 1933+503, a ten-image radio lens, and identify shared properties of the models with the best fits. The approximate rotation curves varies by less than 8.5% from the average value between the innermost and the outermost image (1.5h^{-1} kpc to 4.1h^{-1} kpc) for models within 1 \sigma of the best fit, and the radial dependence of the shear strength and angle also have common behavior for the best models. The time delay between images 1 and 6, the longest delay between the radio cores, is \Delta t = (10.6^{+2.4}_{-1.1})h^{-1} days (\Omega_0=0.3, \lambda_0=0.7) including all the modeling uncertainties. Deeper infrared observations, to more precisely register the lens galaxy with the radio images and to measure the properties of the Einstein ring image of the radio source's host galaxy, would significantly improve the model constraints and further reduce the uncertainties in the mass distribution and time delay., Comment: 24 pages, 10 figures, final version to appear in ApJ. Some minor corrections (e.g. constraint on central unseen image was stronger than intended earlier, now agrees with text, conventions on angles fixed in text/plots). Resulting model fits have some change in chi squareds and best parameters (e.g. cores, flatness of rotation curve) have some changes. Properties of model families and trends for best fitting models very close to earlier results; general conclusions the same

Non-abelian discrete gauge symmetries can provide the inflaton with a flat potential even when one takes into account gravitational strength effects. The discreteness of the symmetries also provide special field values where inflation can end via a hybrid type mechanism. An interesting feature of this method is that it can naturally lead to extremely flat potentials and so, in principle, to inflation at unusually low energy scales. Two examples of effective field theories with this mechanism are given, one with a hybrid exit and one with a mutated hybrid exit. They include an explicit example in which the single field consistency condition is violated., Comment: 24 pages, uses revtex.sty, submitted to PRD (Nov. 1999) Final version to appear in PRD. Background information on supergravity expanded

Obtaining a potential flat enough to provide slow roll inflation is often difficult when gravitational effects are included. Non-abelian discrete gauge symmetries can guarantee the flatness of the inflaton potential in this case, and also provide special field values where inflation can end., Comment: 9 pages, to appear in Physics Letters B

A short pedagogical note about the consequences of a nonzero cosmological constant in physical cosmology., Comment: 16 pages, uses epsfig,prd, aps style files, 5 figures included in text, final version to appear as invited review in Astrophysics and Space Science

Scalar and gravity wave subcurvature fluctuations are calculated for a background approximating the Hawking and Turok open universe model. The gravity wave cosmic microwave background contribution is finite and it appears that a normalizable scalar supercurvature mode is possible in some regions of parameter space., Comment: 10 pages, latex, uses revtex, prd and aps style files, final version to appear in Phys. Rev. D

In the hyperbolic slicing of de Sitter space appropriate for open universe models, a curvature scale is present and supercurvature fluctuations are possible. In some cases, the expansion of a scalar field in the Bunch-Davies vacuum includes supercurvature modes, as shown by Sasaki, Tanaka and Yamamoto. We express the normalizable vacuum supercurvature modes for a massless scalar field in terms of the basis modes for the spatially-flat slicing of de Sitter space., Comment: 23 pages, 3 .eps figures included, uses prd,aps, epsfig style files final version to appear in Phys Rev D, expanded concluding section

An overview of some new results in open inflation over the past year, including the calculation of gravity wave contributions to the Cosmic Microwave Background. Proceedings for COSMO-97, held in Ambleside, U.K., hence very short due to space limitations (the version for the web is more than 3 pages because archive numbers have been added for the references)., Comment: LaTeX, uses sprocl.tex, 4 pages, corrected reference typo

We calculate the spectrum of gravitational waves generated during inflation in open $(\Omega _0<1)$ inflationary models. In such models an initial epoch of old inflation solves the horizon and flatness problems, and during this first epoch of inflation the quantum state of the graviton field rapidly approaches the Bunch-Davies vacuum. Then old inflation ends by the nucleation of a single bubble, inside of which there is a shortened epoch of slow-roll inflation giving $\Omega _0<1$ today. In this paper we re-express the Bunch-Davies vacuum for the graviton field in terms of the hyperbolic modes inside the bubble and propagate these modes forward in time into the present era. We derive the expression for the contribution from these gravity waves to the cosmic microwave background anisotropy including the effect of a finite energy difference across the bubble wall., Comment: 40 pages, TEX with phyzzx macro, 5 figures

An introduction to models of open universes originating from bubbles, including a summary of recent theoretical results for the power spectrum. To appear in the proceedings of the XXXIth Moriond meeting, "Microwave Background Anisotropies.", Comment: LaTeX file, uses epsf.tex, 3 figures, 8 pages, minor wording changes to clarify discussion of coordinate systems

The interior of a vacuum bubble in de Sitter space may give an open universe with sufficient homogeneity to agree with observations. Here, previous work by Bucher, Goldhaber and Turok is extended to describe a thin bubble wall with nonzero radius and energy difference across the wall. The vacuum modes present before formation of the bubble propagate into the interior of open universe and the power spectrum of the resulting gauge invariant gravitational potential is calculated. It is found to become scale invariant on small scales, with onset at about the same scale as that found in the zero radius case. There is sensitivity to the radius and energy difference at large scales, but it is expected that they cannot be strongly constrained because of cosmic variance. As the prediction of a scale invariant spectrum is robust with respect to variation of these parameters at small scales, it is a generic feature of the contribution of these modes for these thin wall models., Comment: 15 pages, LaTeX and revtex macros, final version to appear in Physical Review D, typos and minor wording changes

A strong-weak coupling duality symmetry of the string equations of motion has been suggested in the literature. This symmetry implies that vacua occur in pairs. Since the coupling constant is a dynamical variable in string theory, tunneling solutions between strong and weak coupling vacua may exist. Such solutions would naturally lead to nonperturbative effects with anomalous coupling dependence. A highly simplified example is given., Comment: 7 pages, iassns-hep-93-24, fermilab 93/087-T, nsf-itp-93-45 (A simple example of the mechanism suggested has been added.)

This TASI lecture covers the material in hep-th/9205026. It reviewed the theory of effective strings, with particular emphasis on the manner in which Lorentz invariance is represented. The quantum properties of an example of an effective string are derived from the underlying field theory. A comparison is made with what one would expect if one assumed that quantum effective strings were governed by fundamental string actions such as the Nambu-Goto or the Polyakov actions. It is shown that the requirements on dimensions for consistent quantizations of fundamental strings imply no contradictions for effective strings., Comment: 11 pages(phyzzx), Fermilab 92/293-C, iassns-hep-92-73

Starting from a Poincar\'e invariant field theory of a real scalar field with interactions governed by a double-well potential in 2+1 dimensions, the Lorentz representation induced on the collective coordinates describing low-energy excitations about an effective string background is derived. In this representation, Lorentz transformations are given in terms of an infinite series, in powers of derivatives along the worldsheet. Transformations that act on the direction transverse to the string worldsheet involve a universal dimension $-1$ term. As a consequence, Lorentz invariance holds in this theory of long effective strings due to cancellations in the action between irrelevant terms and the dimension two term that describes free massless scalar fields in two dimensions. (in plain tex, no macropackages necessary), Comment: 9 pages

We derive the supersymmetric collective field theory for the Marinari-Parisi model. For a specific choice of the superpotential, to leading order we find a one parameter family of ground states which can be connected via instantons. At this level of analysis the instanton size implied by the underlying matrix model does not appear., Comment: 14 pages