Counts-in-cells (CIC) measurements contain a wealth of cosmological information yet are seldom used to constrain theories. Although we can predict the shape of the distribution for a given cosmology, to fit a model to the observed CIC probabilities requires the covariance matrix -- both the variance of counts in one probability bin and the covariance between counts in different bins. To date, there have been no general expressions for these variances. Here we show that correlations of particular levels, or "slices," of the density field determine the variance and covariance of CIC probabilities. We derive explicit formulae that accurately predict the variance and covariance among subvolumes of a simulated galaxy catalog, opening the door to the use of CIC measurements for cosmological parameter estimation.

The counts-in-cells (CIC) galaxy probability distribution depends on both the dark matter clustering amplitude $\sigma_8$ and the galaxy bias $b$. We present a theory for the CIC distribution based on a previous prescription of the underlying dark matter distribution and a linear volume transformation to redshift space. We show that, unlike the power spectrum, the CIC distribution breaks the degeneracy between $\sigma_8$ and $b$ on scales large enough that both bias and redshift distortions are still linear; thus we obtain a simultaneous fit for both parameters. We first validate the technique on the Millennium Simulation and then apply it to the SDSS Main Galaxy Sample. We find $\sigma_8 = 0.94^{+.11}_{-.10}$ and $b = 1.36^{+.14}_{-.11}$, consistent with previous complementary results from redshift distortions and from Planck.

Repp and Szapudi (2019) present a physically-motivated galaxy bias model which remains physical in low-density regions and which also provides a better fit to simulation data than do typical survey-analysis bias models. Given plausible simplifying assumptions, the physics of this model (surprisingly) proves to be analogous to the Ising model of statistical mechanics. In the present work we present a method of testing this Ising bias model against empirical galaxy survey data. Using this method, we compare our model (as well as three reference models -- linear, quadratic, and logarithmic) to SDSS, 6dFGS, and COSMOS2015 results, finding that for spectroscopic redshift surveys, the Ising bias model provides a superior fit compared to the reference models. Photometric redshifts, on the other hand, introduce enough error into the radial coordinate that none of the models yields a good fit. A physically meaningful galaxy bias model is necessary for optimal extraction of cosmological information from dense galaxy surveys such as Euclid and WFIRST.

We perform an analysis of two-point galaxy clustering and galaxy bias using Subaru Hyper-Suprime Cam (HSC) data taken jointly by the Subaru Strategic Program and the University of Hawaii in the COSMOS field. The depth of the data is similar to the ongoing Hawaii Two-0 (H20) optical galaxy survey, thus the results are indicative of future constraints from tenfold area. We measure the angular auto-power spectra of the galaxy overdensity in three redshift bins, defined by dropouts from the g-, r- and i-bands, and compare them to the theoretical expectation from concordance cosmology with linear galaxy bias. We determine the redshift distribution of each bin using a standard template-based photometric redshift method, coupled with a self-organizing map (SOM) to quantify colour space coverage. We also investigate sources of systematic errors to inform the methodology and requirements for Hawaii Two-0. The linear galaxy bias fit results are $b_{\mathrm{gal,g}} = 3.90 \pm 0.33 (\mathrm{stat}) \substack{ +0.64 \\ -0.24 } (\mathrm{sys})$ at redshift $z \simeq 3.7$, $b_{\mathrm{gal,r}} = 8.44 \pm 0.63 (\mathrm{stat}) \substack{ +1.42 \\ -0.72 } (\mathrm{sys})$ at $z \simeq 4.7$, and $b_{\mathrm{gal,i}} = 11.94 \pm 2.24 (\mathrm{stat}) \substack{ +1.82 \\ -1.27 } (\mathrm{sys})$ at $z \simeq 5.9$.

A reliable model of galaxy bias is necessary for interpreting data from future dense galaxy surveys. Conventional bias models are inaccurate, in that they can yield unphysical results ($\delta_g < -1$) for voids that might contain half of the available cosmological information. For this reason, we present a physically-motivated bias model based on an analogy with the Ising model. With only two free parameters, the model produces sensible results for both high- and low-density regions. We also test the model using a catalog of Millennium Simulation galaxies in cubical survey pixels with side lengths from $2h^{-1}$--$31h^{-1}$Mpc, at redshifts from 0 to 2. We find the Ising model markedly superior to linear and quadratic bias models on scales smaller than $10h^{-1}$Mpc, while those conventional models fare better on scales larger than $30h^{-1}$Mpc. While the largest scale where the Ising model is applicable might vary for a specific galaxy catalog, it should be superior on any scale with a non-negligible fraction of cells devoid of galaxies.

Future galaxy surveys hope to realize significantly tighter constraints on various cosmological parameters. The higher number densities achieved by these surveys will allow them to probe the smaller scales affected by non-linear clustering. However, in these regimes, the standard power spectrum can extract only a portion of such surveys' cosmological information. In contrast, the alternate statistic $A^*$ has the potential to double these surveys' information return, provided one can predict the $A^*$-power spectrum for a given cosmology. Thus, in this work we provide a prescription for this power spectrum $P_{A^*}(k)$, finding that the prescription is typically accurate to about 5 per cent for near-concordance cosmologies. This prescription will thus allow us to multiply the information gained from surveys such as Euclid and WFIRST.

Evaluation of gravitational theories by means of cosmological data suffers from the fact that galaxies are biased tracers of dark matter. Current bias models focus primarily on high-density regions, whereas low-density regions carry significant amounts of information relevant to the constraint of dark energy and alternative gravity theories. Thus, proper treatment of both high and low densities is important for future surveys. Accordingly, we here present an interactionless Ising model for this bias, and we demonstrate that it exhibits a remarkably good fit to both Millennium Simulation and Sloan Digital Sky Survey data, at both density extremes. The quality of the fit indicates that galaxy formation is (to first order) an essentially local process determined by initial conditions.

The distribution of matter in the universe is, to first order, lognormal. Improving this approximation requires characterization of the third moment (skewness) of the log density field. Thus, using Millennium Simulation phenomenology and building on previous work, we present analytic fits for the mean, variance, and skewness of the log density field $A$. We further show that a Generalized Extreme Value (GEV) distribution accurately models $A$; we submit that this GEV behavior is the result of strong intrapixel correlations, without which the smoothed distribution would tend (by the Central Limit Theorem) toward a Gaussian. Our GEV model yields cumulative distribution functions accurate to within 1.7 per cent for near-concordance cosmologies, over a range of redshifts and smoothing scales.

A primary goal of galaxy surveys is to tighten constraints on cosmological parameters, and the power spectrum $P(k)$ is the standard means of doing so. However, at translinear scales $P(k)$ is blind to much of these surveys' information---information which the log density power spectrum recovers. For discrete fields (such as the galaxy density), $A^*$ denotes the statistic analogous to the log density: $A^*$ is a "sufficient statistic" in that its power spectrum (and mean) capture virtually all of a discrete survey's information. However, the power spectrum of $A^*$ is biased with respect to the corresponding log spectrum for continuous fields, and to use $P_{A^*}(k)$ to constrain the values of cosmological parameters, we require some means of predicting this bias. Here we present a prescription for doing so; for Euclid-like surveys (with cubical cells 16$h^{-1}$ Mpc across) our bias prescription's error is less than 3 per cent. This prediction will facilitate optimal utilization of the information in future galaxy surveys.

Hubble Space Telescope SNAPshot surveys of 86 X-ray selected galaxy clusters at $0.3 < z < 0.5$ from the MACS sample have proven invaluable for the exploration of a wide range of astronomical research topics. We here present an overview of the four MACS SNAPshot surveys conducted from Cycle 14 to Cycle 20 as part of a long-term effort aimed at identifying exceptional cluster targets for in-depth follow up by the extragalactic community. We also release redshifts and X-ray luminosities of all clusters observed as part of this initiative. To illustrate the power of SNAPshot observations of MACS clusters, we explore several aspects of galaxy evolution illuminated by the images obtained for these programmes. We confirm the high lensing efficiency of X-ray selected clusters at $z>0.3$. Examining the evolution of the slope of the cluster red sequence, we observe at best a slight decrease with redshift, indicating minimal age contribution since $z\sim 1$. Congruent to previous studies' findings, we note that the two BCGs which are significantly bluer ($\geq 5\sigma$) than their clusters' red sequences reside in relaxed clusters and exhibit pronounced internal structure. Thanks to our targets' high X-ray luminosity, the subset of our sample observed with Chandra adds valuable leverage to the X-ray luminosity--optical richness relation, which, albeit with substantial scatter, is now clearly established from groups to extremely massive clusters of galaxies. We conclude that SNAPshot observations of MACS clusters stand to continue to play a vital pathfinder role for astrophysical investigations across the entire electromagnetic spectrum.

We present the results of a 135-arcmin$^2$ search for high-redshift galaxies lensed by 29 clusters from the MAssive Cluster and extended MAssive Cluster Surveys (MACS and eMACS). We use relatively shallow images obtained with the Hubble Space Telescope in four passbands, namely, F606W, F814W, F110W, and F140W. We identify 130 F814W dropouts as candidates for galaxies at $z \le 6$. In order to fit the available broad-band photometry to galaxy spectral energy distribution (SED) templates, we develop a prior for the level of dust extinction at various redshifts. We also investigate the systematic biases incurred by the use of SED-fit software. The fits we obtain yield an estimate of 20 Lyman-break galaxies with photometric redshifts from $z \sim 7$ to 9. In addition, our survey has identified over 100 candidates with a significant probability of being lower-redshift ($z \sim 2$) interlopers. We conclude that even as few as four broad-band filters -- when combined with fitting the SEDs -- are capable of isolating promising objects. Such surveys thus allow one both to probe the bright end ($M_{1500} \le -19$) of the high-redshift UV luminosity function and to identify candidate massive evolved galaxies at lower redshifts.

At translinear scales, the log power spectrum captures significantly more cosmological information than the standard power spectrum. At high wavenumbers $k$, the Fisher information in the standard power spectrum $P(k)$ fails to increase in proportion to $k$ in part due to correlations between large- and small-scale modes. As a result, $P(k)$ suffers from an information plateau on these translinear scales, so that analysis with the standard power spectrum cannot access the information contained in these small-scale modes. The log power spectrum $P_A(k)$, on the other hand, captures the majority of this otherwise lost information. Until now there has been no means of predicting the amplitude of the log power spectrum apart from cataloging the results of simulations. We here present a cosmology-independent prescription for the log power spectrum; this prescription displays accuracy comparable to that of Smith et al. (2003), over a range of redshifts and smoothing scales, and for wavenumbers up to $1.5h$ Mpc$^{-1}$.

The primary science driver for 3D galaxy surveys is their potential to constrain cosmological parameters. Forecasts of these surveys' effectiveness typically assume Gaussian statistics for the underlying matter density, despite the fact that the actual distribution is decidedly non-Gaussian. To quantify the effect of this assumption, we employ an analytic expression for the power spectrum covariance matrix to calculate the Fisher information for BAO-type model surveys. We find that for typical number densities, at $k_\mathrm{max} = 0.5 h$ Mpc$^{-1}$, Gaussian assumptions significantly overestimate the information on all parameters considered, in some cases by up to an order of magnitude. However, after marginalizing over a six-parameter set, the form of the covariance matrix (dictated by $N$-body simulations) causes the majority of the effect to shift to the "amplitude-like" parameters, leaving the others virtually unaffected. We find that Gaussian assumptions at such wavenumbers can underestimate the dark energy parameter errors by well over 50 per cent, producing dark energy figures of merit almost 3 times too large. Thus, for 3D galaxy surveys probing the non-linear regime, proper consideration of non-Gaussian effects is essential.