In the local (redshift z~0) Universe, collisional ring galaxies make up only ~0.01% of galaxies and are formed by head-on galactic collisions that trigger radially propagating density waves. These striking systems provide key snapshots for dissecting galactic disks and are studied extensively in the local Universe. However, not much is known about distant (z>0.1) collisional rings. Here we present a detailed study of a ring galaxy at a look-back time of 10.8 Gyr (z=2.19). Compared with our Milky Way, this galaxy has a similar stellar mass, but has a stellar half-light radius that is 1.5-2.2 times larger and is forming stars 50 times faster. The large, diffuse stellar light outside the star-forming ring, combined with a radial velocity on the ring and an intruder galaxy nearby, provides evidence for this galaxy hosting a collisional ring. If the ring is secularly evolved, the implied large bar in a giant disk would be inconsistent with the current understanding of the earliest formation of barred spirals. Contrary to previous predictions, this work suggests that massive collisional rings were as rare 11 Gyr ago as they are today. Our discovery offers a unique pathway for studying density waves in young galaxies, as well as constraining the cosmic evolution of spiral disks and galaxy groups.

To understand how strong emission line galaxies (ELGs) contribute to the overall growth of galaxies and star formation history of the universe, we target Strong ELGs (SELGs) from the ZFOURGE imaging survey that have blended (Hb+[OIII]) rest-frame equivalent widths of >230A and 2.5<zphot<4.0. Using Keck/MOSFIRE, we measure 49 redshifts for galaxies brighter than Ks=25 mag as part of our Multi-Object Spectroscopic Emission Line (MOSEL) survey. Our spectroscopic success rate is ~53% and zphot uncertainty is sigma_z= [Delta(z)/(1+z)]=0.0135. We confirm 31 ELGs at 3<zspec<3.8 and show that Strong ELGs have spectroscopic rest-frame [OIII]5007A equivalent widths of 100-500A and tend to be lower mass systems [log(Mstar/Msun)~8.2-9.6] compared to more typical star-forming galaxies. The Strong ELGs lie ~0.9 dex above the star-forming main-sequence at z~3.5 and have high inferred gas fractions of fgas~>60%, i.e. the inferred gas masses can easily fuel a starburst to double stellar masses within ~10-100 Myr. Combined with recent results using ZFOURGE, our analysis indicates that 1) strong [OIII]5007A emission signals an early episode of intense stellar growth in low mass (Mstar<0.1M*) galaxies and 2) many, if not most, galaxies at z>3 go through this starburst phase. If true, low-mass galaxies with strong [OIII]5007A emission (EW_rest>200A) may be an increasingly important source of ionizing UV radiation at z>3.

We use K-band spectroscopic data from the Multi-Object Spectroscopic Emission Line (MOSEL) survey to analyze the kinematic properties of galaxies at z>3. Our sample consists of 34 galaxies at 3.0<zspec<3.8 between 9.0<log(M_star)<11.0. We find that galaxies with log(M_star) > 10.2 at z > 3 have 56 +/- 21 km/s lower integrated velocity dispersion compared to galaxies at z ~ 2 of similar stellar mass. Massive galaxies at z > 3 have either a flat or declining star formation history (SFH), whereas similar stellar mass galaxies at z~2.0 exhibit a slight peak in the past 500 Myrs. Comparing with the IllustrisTNG cosmological simulation, we find that (i) the dynamical mass of massive galaxies in simulations (log(M_star) > 10.0) increases by ~0.1 dex at a fixed stellar mass between z=2.0-3.0, and (ii) dynamical mass growth is coupled with a rapid rise in the ex-situ stellar mass fraction (stars accreted from other galaxies) for massive galaxies at z < 3.5. We speculate that the rising contribution of ex-situ stellar mass to the total stellar mass growth of massive galaxies is driving the higher integrated velocity dispersion and rising SFHs of massive galaxies at z~2.0 compared to galaxies of similar stellar masses at z > 3.

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.

The ZFIRE survey has spectroscopically confirmed two proto-clusters using the MOSFIRE instrument on Keck 1: one at z=2.095 in COSMOS and another at z=1.62 in UDS. Here we use an updated ZFIRE dataset to derive the properties of ionized gas regions of proto-cluster galaxies by extracting fluxes from emission lines H$\beta$ 4861\AA, [OIII] 5007\AA, H$\alpha$ 6563\AA, [NII] 6585\AA, and [SII] 6716,6731\AA. We measure gas-phase metallicity of members in both proto-clusters using two indicators, including a strong-line indicator relatively independent of ionization parameter and electron density. Proto-cluster and field galaxies in both UDS and COSMOS lie on the same Mass-Metallicity Relation with both metallicity indicators. We compare our results to recent IllustrisTNG results, which reports no significant gas-phase metallicity offset between proto-cluster and field galaxies until z=1.5. This is in agreement with our observed metallicities, where no offset is measured between proto-cluster and field populations. We measure tentative evidence from stacked spectra that UDS high mass proto-cluster and field galaxies have differing [OIII]/H$\beta$ ratios, however these results are dependent on the sample size of the high mass stacks.

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.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

The "cosmic dawn" refers to the period of the Universe's history when stars and black holes first formed and began heating and ionizing hydrogen in the intergalactic medium (IGM). Though exceedingly difficult to detect directly, the first stars and black holes can be constrained indirectly through measurements of the cosmic 21-cm background, which traces the ionization state and temperature of intergalactic hydrogen gas. In this white paper, we focus on the science case for such observations, in particular those targeting redshifts z $\gtrsim$ 10 when the IGM is expected to be mostly neutral. 21-cm observations provide a unique window into this epoch and are thus critical to advancing first star and black hole science in the next decade.

In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $\Lambda$CDM paradigm at high redshifts.

The Dark Ages are the period between the last scattering of the cosmic microwave background and the appearance of the first luminous sources, spanning approximately 1100 < z < 30. The only known way to measure fluctuations in this era is through the 21-cm line of neutral hydrogen. Such observations have enormous potential for cosmology, because they span a large volume while the fluctuations remain linear even on small scales. Observations of 21-cm fluctuations during this era can therefore constrain fundamental aspects of our Universe, including inflation and any exotic physics of dark matter. While the observational challenges to these low-frequency 21-cm observations are enormous, especially from the terrestrial environment, they represent an important goal for cosmology.

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

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.

Using the \prospector\ spectral energy distribution (SED) fitting code, we analyze the properties of 19 Extreme Emission Line Galaxies (EELGs) identified in the bluest composite SED in the \zfourge\ survey at $2.5 \leq z \leq 4$. \prospector\ includes a physical model for nebular emission and returns probability distributions for stellar mass, stellar metallicity, dust attenuation, and nonparametric star formation history (SFH). The EELGs show evidence for a starburst in the most recent 50 Myr, with the median EELG having a specific star formation rate (sSFR) of 4.6 Gyr$^{-1}$ and forming 15\% of its mass in this short time. For a sample of more typical star-forming galaxies (SFGs) at the same redshifts, the median SFG has a sSFR of 1.1 Gyr$^{-1}$ and forms only $4\%$ of its mass in the last 50 Myr. We find that virtually all of our EELGs have rising SFHs, while most of our SFGs do not. From our analysis, we hypothesize that many, if not most, star-forming galaxies at $z \geq 2.5$ undergo an extreme H$\beta$+$[\hbox{{\rm O}\kern 0.1em{\sc iii}}]$ emission line phase early in their lifetimes. In a companion paper, we obtain spectroscopic confirmation of the EELGs as part of our {\sc MOSEL} survey. In the future, explorations of uncertainties in modeling the UV slope for galaxies at $z>2$ are needed to better constrain their properties, e.g. stellar metallicities.

We investigate the properties of galaxies as they shut off star formation over the 4 billion years surrounding peak cosmic star formation. To do this we categorize $\sim7000$ galaxies from $1<z<4$ into $90$ groups based on the shape of their spectral energy distributions (SEDs) and build composite SEDs with $R\sim 50$ resolution. These composite SEDs show a variety of spectral shapes and also show trends in parameters such as color, mass, star formation rate, and emission line equivalent width. Using emission line equivalent widths and strength of the 4000\AA\ break, $D(4000)$, we categorize the composite SEDs into five classes: extreme emission line, star-forming, transitioning, post-starburst, and quiescent galaxies. The transitioning population of galaxies show modest H$\alpha$ emission ($EW_{\rm REST}\sim40$\AA) compared to more typical star-forming composite SEDs at $\log_{10}(M/M_\odot)\sim10.5$ ($EW_{\rm REST}\sim80$\AA). Together with their smaller sizes (3 kpc vs. 4 kpc) and higher S\'ersic indices (2.7 vs. 1.5), this indicates that morphological changes initiate before the cessation of star formation. The transitional group shows a strong increase of over one dex in number density from $z\sim3$ to $z\sim1$, similar to the growth in the quiescent population, while post-starburst galaxies become rarer at $z\lesssim1.5$. We calculate average quenching timescales of 1.6 Gyr at $z\sim1.5$ and 0.9 Gyr at $z\sim2.5$ and conclude that a fast quenching mechanism producing post-starbursts dominated the quenching of galaxies at early times, while a slower process has become more common since $z\sim2$.

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.

The recent discovery of gravitational waves (GWs) from stellar-mass binary black holes (BBHs) provided direct evidence of the existence of these systems. BBH lenses would have gravitational microlensing signatures that are distinct from single-lens signals. We apply Bayesian statistics to examine the distinguishability of BBH microlensing events from single-lens events under ideal observing conditions, using the photometric capabilities of the Korean Microlensing Telescope Network. Given one year of observations, a source star at the Galactic Centre, a GW150914-like BBH lens (total mass 65M$_\odot$, mass ratio 0.8) at half that distance, and an impact parameter of 0.4 Einstein radii, we find that binarity is detectable for BBHs with separations down to 0.0250 Einstein radii, which is nearly 3.5 times greater than the maximum separation for which such BBHs would merge within the age of the Universe. Microlensing searches are thus sensitive to more widely separated BBHs than GW searches, perhaps allowing the discovery of BBH populations produced in different channels of binary formation.

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.

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.

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.

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.

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.

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.

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.

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]

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)

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

The cold dark matter paradigm has been extremely successful for explaining a wide range of cosmological phenomena. Nevertheless, since evidence for non-baryonic dark matter remains indirect, all reasonable alternatives should be explored. One proposed idea, involving a fundamental acceleration scale a_0 = 1-2e-10 m/s^2, is called MOdified Newtonian Dynamics or MOND. MOND was suggested to explain the flat rotation curves of galaxies without the need for dark matter. Whether or not it can adequately fit the available data has been debated for almost 20 years (and we summarise many of these studies), but only recently have there been attempts to extend MOND to larger scale regimes. We discuss how the basic properties of MOND make it at best ambiguous to apply these ideas to cosmological scales. We emphasize the difficulties inherent in developing a full theory in which to embed MOND. Without such a theory there is no obviously consistent way to discuss the early Universe and the growth of perturbations. Recent claims that MONDian cosmology works very well are therefore not supportable. We also provide an argument for why a_0~cH_0 naturally, a coincidence which is often suggested as a motivation for taking MOND seriously. We discuss other alternative theories of gravity concluding, as others have, that no metric theory extensions appear workable for explaining rotation curves as well as other observed phenomena. The whole premise of many of these attempts is fatally flawed - galaxies are not pre-selected, discrete, isolated regions which formed monolithically and around which one can construct an axially-symmetric dynamical model in order to remove the need for dark matter. In the modern view, galaxies are part of a dynamic continuum of objects which collectively make up the evolving large-scale structure of the universe.

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.

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.

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.

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.

A short pedagogical note about the consequences of a nonzero cosmological constant in physical cosmology.

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.

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.

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).

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.

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."

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.