In a galactic halo like the Milky Way, bosonic dark matter particles lighter than about $100$ eV have a de Broglie wavelength larger than the average inter-particle separation and are therefore well described as a set of classical waves. This applies to, for instance, the QCD axion as well as to lighter axion-like particles such as fuzzy dark matter. We show that the interference of waves inside a halo inevitably leads to vortices, locations where chance destructive interference takes the density to zero. The phase of the wavefunction has non-trivial winding around these points. This can be interpreted as a non-zero velocity circulation, so that vortices are sites where the fluid velocity has a non-vanishing curl. Using analytic arguments and numerical simulations, we study the properties of vortices and show they have a number of universal features: (1) In three spatial dimensions, the generic defects take the form of vortex rings. (2) On average there is about one vortex ring per de Broglie volume and (3) generically only single winding ($\pm 1$) vortices are found in a realistic halo. (4) The density near a vortex scales as $r^2$ while the velocity goes as $1/r$, where $r$ is the distance to vortex. (5) A vortex segment moves at a velocity inversely proportional to its curvature scale so that smaller vortex rings move faster, allowing momentary motion exceeding escape velocity. We discuss observational/experimental signatures from vortices and, more broadly, wave interference. In the ultra-light regime, gravitational lensing by interference substructures leads to flux anomalies of $5-10 \%$ in strongly lensed systems. For QCD axions, vortices lead to a diminished signal in some detection experiments but not in others. We advocate the measurement of correlation functions by axion detection experiments as a way to probe and capitalize on the expected interference substructures.

It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with $m \sim 10^{-22}$ eV and often called "fuzzy" dark matter, provides intriguing solutions to some challenges at sub-Galactic scales for the standard cold dark matter model. As shown by Khmelnitsky and Rubakov, such a scalar field in the Galaxy would produce an oscillating gravitational potential with nanohertz frequencies, resulting in periodic variations in the times of arrival of radio pulses from pulsars. The Parkes Pulsar Timing Array (PPTA) has been monitoring 20 millisecond pulsars at two to three weeks intervals for more than a decade. In addition to the detection of nanohertz gravitational waves, PPTA offers the opportunity for direct searches for fuzzy dark matter in an astrophysically feasible range of masses. We analyze the latest PPTA data set which includes timing observations for 26 pulsars made between 2004 and 2016. We perform a search in this data set for evidence of ultralight dark matter in the Galaxy using Bayesian and Frequentist methods. No statistically significant detection has been made. We therefore place upper limits on the local dark matter density. Our limits, improving on previous searches by a factor of two to five, constrain the dark matter density of ultralight bosons with $m \leq 10^{-23}$ eV to be below $6\,\text{GeV}\,\text{cm}^{-3}$ with 95\% confidence in the Earth neighborhood. Finally, we discuss the prospect of probing the astrophysically favored mass range $m \gtrsim 10^{-22}$ eV with next-generation pulsar timing facilities.

Quasinormal modes describe the return to equilibrium of a perturbed system, in particular the ringdown phase of a black hole merger. But as globally-defined quantities, the quasinormal spectrum can be highly sensitive to global structure, including distant small perturbations to the potential. In what sense are quasinormal modes a property of the resulting black hole? We explore this question for the linearized perturbation equation with two potentials having disjoint bounded support. We give a composition law for the Wronskian that determines the quasinormal frequencies of the combined system. We show that over short time scales the evolution is governed by the quasinormal frequencies of the individual potentials, while the sensitivity to global structure can be understood in terms of echoes. We introduce an echo expansion of the Green's function and show that, as expected on general grounds, at any finite time causality limits the number of echoes that can contribute. We illustrate our results with the soluble example of a pair of $\delta$-function potentials. We explicate the causal structure of the Green's function, demonstrating under what conditions two very different quasinormal spectra give rise to very similar ringdown waveforms.

In this paper, we verify the large scale structure consistency relations using N-body simulations, including modes in the highly non-linear regime. These relations (pointed out by Kehagias & Riotto and Peloso & Pietroni) follow from the symmetry of the dynamics under a shift of the Newtonian potential by a constant and a linear gradient, and predict the absence of certain poles in the ratio between the (equal time) squeezed bispectrum and power spectrum. The consistency relations, as symmetry statements, are exact, but have not been previously checked beyond the perturbative regime. Our test using N-body simulations not only offers a non-perturbative check, but also serves as a warm-up exercise for applications to observational data. A number of subtleties arise when taking the squeezed limit of the bispectrum--we show how to circumvent or address them. An interesting by-product of our investigation is an explicit demonstration that the linear-gradient symmetry is unaffected by the periodic boundary condition of the simulations. Lastly, we verify using simulations that the consistency relations are violated when the initial conditions are non-gaussian (of the local fNL type). The methodology developed here paves the way for constraining primordial non-gaussianity using large scale structure data, including (numerous) highly non-linear modes that are otherwise hard to interpret and utilize.

We show that a black hole surrounded by scalar dark matter develops scalar hair. This is the generalization of a phenomenon pointed out by Jacobson, that a minimally coupled scalar with a non-trivial time dependence far away from the black hole would endow the black hole with hair. In our case, the time dependence arises from the oscillation of a scalar field with a non-zero mass. We systematically explore the scalar profile around the black hole for different scalar masses. In the small mass limit, the scalar field has a $1/r$ component at large radius $r$, consistent with Jacobson's result. In the large mass limit (with the Compton wavelength of order of the horizon or smaller), the scalar field has a $1/r^{3/4}$ profile yielding a pile-up close to the horizon, while distinctive nodes occur for intermediate masses. Thus, the dark matter profile around a black hole, while challenging to measure, contains information about the dark matter particle mass. As an application, we consider the case of the supermassive black hole at the center of M87, recently imaged by the Event Horizon Telescope. Its horizon size is roughly the Compton wavelength of a scalar particle of mass $10^{-20}$ eV. We consider the implications of the expected scalar pile-up close to the horizon, for fuzzy dark matter at a mass of $10^{-20}$ eV or below.

We reconsider the derivation of soft theorems associated with nonlinearly-realized symmetries in cosmology. Utilizing the path integral, we derive a generalized consistency relation that relates a squeezed $(N+1)$-point correlation function to an $N$-point function, where the relevant soft mode is at early rather than late time. This generalized (early-late-time) version has wider applicability than the standard consistency relation where all modes are evaluated at late times. We elucidate the conditions under which the latter follows from the former. A key ingredient is the physical mode condition: that the nonlinear part of the symmetry transformation must match the time dependence of the dominant, long wavelength physical mode. This is closely related to, but distinct from, the adiabatic mode condition. Our derivation sheds light on a number of otherwise puzzling features of the standard consistency relation: (1) the underlying nonlinearly-realized symmetries (such as dilation and special conformal transformation SCT) originate as residual gauge redundancies, yet the consistency relation has physical content---for instance, it can be violated; (2) the standard consistency relation is known to fail in ultra-slow-roll inflation, but since dilation and SCT remain good symmetries, there should be a replacement for the standard relation; (3) in large scale structure applications, it is known that the standard consistency relation breaks down if the long wavelength power spectrum is too blue. The early-late-time consistency relation helps address these puzzles. We introduce a toy model where explicit checks of this generalized consistency relation are simple to carry out. Our methodology can be adapted to cases where violations of the standard consistency relation involve additional light degrees of freedom beyond the inflaton.

The final ringdown phase in a coalescence process is a valuable laboratory to test General Relativity and potentially constrain additional degrees of freedom in the gravitational sector. We introduce here an effective description for perturbations around spherically symmetric spacetimes in the context of scalar-tensor theories, which we apply to study quasi-normal modes for black holes with scalar hair. We derive the equations of motion governing the dynamics of both the polar and the axial modes in terms of the coefficients of the effective theory. Assuming the deviation of the background from Schwarzschild is small, we use the WKB method to introduce the notion of "light ring expansion". This approximation is analogous to the slow-roll expansion used for inflation, and it allows us to express the quasinormal mode spectrum in terms of a small number of parameters. This work is a first step in describing, in a model independent way, how the scalar hair can affect the ringdown stage and leave signatures on the emitted gravitational wave signal. Potential signatures include the shifting of the quasi-normal spectrum, the breaking of isospectrality between polar and axial modes, and the existence of scalar radiation.

We investigate nonlinear structure formation in the fuzzy dark matter (FDM) model using both numerical and perturbative techniques. On the numerical side, we examine the virtues and limitations of a Schrodinger-Poisson solver (wave formulation) versus a fluid dynamics solver (Madelung formulation). We also carry out a perturbative computation of the one-loop mass power spectrum. We find that (1) in many cases, the fluid dynamics solver is capable of producing the expected interference patterns, but it fails where destructive interference causes the density to vanish which generally occurs in the nonlinear regime. (2) The Schrodinger-Poisson solver works well in all test cases, but it is demanding in resolution: one must resolve the small de Broglie scale to obtain the correct dynamics on large scales. (3) We compare the mass power spectrum from perturbation theory against that from the Schrodinger-Poisson solver, and find good agreement in the mildly nonlinear regime. Compared with fluid perturbation theory, wave perturbation theory has a more limited range of validity. (4) As an application, we compare the Lyman-alpha forest flux power spectrum obtained from the Schrodinger-Poisson solver versus one from an N-body simulation (which is often used as an approximate method to make predictions for FDM). At redshift 5, the two, starting from the same initial condition, agree to better than 10 % on observationally relevant scales as long as the FDM mass exceeds $2 \times 10^{-23}$ eV.

The EDGES experiment recently announced evidence for a broad absorption feature in the sky-averaged radio spectrum around 78 MHz, as may result from absorption in the 21 cm line by neutral hydrogen at z~15-20. If confirmed, one implication is that the spin temperature of the 21 cm line is coupled to the gas temperature by z=20. The known mechanism for accomplishing this is the Wouthuysen-Field effect, whereby Lyman-alpha photons scatter in the intergalactic medium (IGM) and impact the hyperfine level populations. This suggests that early star formation had already produced a copious Lyman-alpha background by z=20, and strongly constrains models in which the linear matter power spectrum is suppressed on small-scales, since halo and star formation are delayed in such scenarios. Here we consider the case that the dark matter consists of ultra-light axions with macroscopic de Broglie wavelengths (fuzzy dark matter, FDM). We assume that star formation tracks halo formation and adopt two simple models from the current literature for the halo mass function in FDM. We further suppose that the fraction of halo baryons which form stars is less than a conservative upper limit of $f_\star \leq 0.05$, and that ~10^4 Lyman-alpha to Lyman-limit photons are produced per stellar baryon. We find that the requirement that the 21 cm spin temperature is coupled to the gas temperature by $z=20$ places a lower-limit on the FDM particle mass of $m_a \geq 5 \times 10^{-21} {\rm eV}$. The constraint is insensitive to the precise minimum mass of halos where stars form. As the global 21 cm measurements are refined, the coupling redshift could change and we quantify how the FDM constraint would be modified. A rough translation of the FDM mass bound to a thermal relic warm dark matter (WDM) mass bound is also provided.

Scalar-tensor theories of gravity generally violate the strong equivalence principle, namely compact objects have a suppressed coupling to the scalar force, causing them to fall slower. A black hole is the extreme example where such a coupling vanishes, i.e. black hole has no scalar hair. Following earlier work, we explore observational scenarios for detecting strong equivalence principle violation, focusing on galileon gravity as an example. For galaxies in-falling towards galaxy clusters, the supermassive black hole can be offset from the galaxy center away from the direction of the cluster. Hence, well resolved images of galaxies around nearby clusters can be used to identify the displaced black hole via the star cluster bound to it. We show that this signal is accessible with imaging surveys, both ongoing ones such as the Dark Energy Survey, and future ground and space based surveys. Already, the observation of the central black hole in M~87 places new constraints on the galileon parameters, which we present here. $\mathcal{O}(1)$ matter couplings are disfavored for a large region of the parameter space. We also find a novel phenomenon whereby the black hole can escape the galaxy completely in less than one billion years.

An intriguing alternative to cold dark matter (CDM) is that the dark matter is a light ( $m \sim 10^{-22}$ eV) boson having a de Broglie wavelength $\lambda \sim 1$ kpc, often called fuzzy dark matter (FDM). We describe the arguments from particle physics that motivate FDM, review previous work on its astrophysical signatures, and analyze several unexplored aspects of its behavior. In particular, (i) FDM halos smaller than about $10^7 (m/10^{-22} {\rm eV})^{-3/2} M_\odot$ do not form. (ii) FDM halos are comprised of a core that is a stationary, minimum-energy configuration called a "soliton", surrounded by an envelope that resembles a CDM halo. (iii) The transition between soliton and envelope is determined by a relaxation process analogous to two-body relaxation in gravitating systems, which proceeds as if the halo were composed of particles with mass $\sim \rho\lambda^3$ where $\rho$ is the halo density. (iv) Relaxation may have substantial effects on the stellar disk and bulge in the inner parts of disk galaxies. (v) Relaxation can produce FDM disks but an FDM disk in the solar neighborhood must have a half-thickness of at least $300 (m/10^{-22} {\rm eV})^{-2/3}$ pc. (vi) Solitonic FDM sub-halos evaporate by tunneling through the tidal radius and this limits the minimum sub-halo mass inside 30 kpc of the Milky Way to roughly $10^8 (m/10^{-22} {\rm eV})^{-3/2} M_\odot$. (vii) If the dark matter in the Fornax dwarf galaxy is composed of CDM, most of the globular clusters observed in that galaxy should have long ago spiraled to its center, and this problem is resolved if the dark matter is FDM.

It has been shown recently that relativistic distortions generate a dipolar modulation in the two-point correlation function of galaxies. To measure this relativistic dipole it is necessary to cross-correlate different populations of galaxies with for example different luminosities or colours. In this paper, we construct an optimal estimator to measure the dipole with multiple populations. We show that this estimator increases the signal-to-noise of the dipole by up to 35 percent. Using 6 populations of galaxies, in a survey with halos and number densities similar to those of the millennium simulation, we forecast a cumulative signal-to-noise of 4.4. For the main galaxy sample of SDSS at low redshift z<0.2 our optimal estimator predicts a cumulative signal-to-noise of 2.4. Finally we forecast a cumulative signal-to-noise of 7.4 in the upcoming DESI survey. These forecasts indicate that with the appropriate choice of estimator the relativistic dipole should be detectable in current and future surveys.

It is usually assumed that in the linear regime the two-point correlation function of galaxies contains only a monopole, quadrupole and hexadecapole. Looking at cross-correlations between different populations of galaxies, this turns out not to be the case. In particular, the cross-correlations between a bright and a faint population of galaxies contain also a dipole. In this paper we present the first attempt to measure this dipole. We discuss the four types of effects that contribute to the dipole: relativistic distortions, evolution effect, wide-angle effect and large-angle effect. We show that the first three contributions are intrinsic anti-symmetric contributions that do not depend on the choice of angle used to measure the dipole. On the other hand the large-angle effect appears only if the angle chosen to extract the dipole breaks the symmetry of the problem. We show that the relativistic distortions, the evolution effect and the wide-angle effect are too small to be detected in the LOWz and CMASS sample of the BOSS survey. On the other hand with a specific combination of angles we are able to measure the large-angle effect with high significance. We emphasise that this large-angle dipole does not contain new physical information, since it is just a geometrical combination of the monopole and the quadrupole. However this measurement, which is in excellent agreement with theoretical predictions, validates our method for extracting the dipole from the two-point correlation function and it opens the way to the detection of relativistic effects in future surveys like e.g. DESI.

By examining the locations of central black holes in two elliptical galaxies, M\,32 and M\,87, we derive constraints on the violation of the strong equivalence principle for purely gravitational objects, i.e. black holes, of less than about two-thirds, $\eta_N<0.68$ from the gravitational interaction of M\,87 with its neighbours in the Virgo cluster. Although M\,32 appears to be a good candidate for this technique, the high concentration of stars near its centre substantially weakens the constraints. On the other hand, if a central black hole is found in NGC 205 or one of the other satellite ellipticals of M\,31, substantially better constraints could be obtained. In all cases the constraints could improve dramatically with better astrometry.

Lensing peaks have been proposed as a useful statistic, containing cosmological information from non-Gaussianities that is inaccessible from traditional two-point statistics such as the power spectrum or two-point correlation functions. Here we examine constraints on cosmological parameters from weak lensing peak counts, using the publicly available data from the 154 deg$^2$ CFHTLenS survey. We utilize a new suite of ray-tracing N-body simulations on a grid of 91 cosmological models, covering broad ranges of the three parameters $\Omega_m$, $\sigma_8$, and $w$, and replicating the Galaxy sky positions, redshifts, and shape noise in the CFHTLenS observations. We then build an emulator that interpolates the power spectrum and the peak counts to an accuracy of $\leq 5\%$, and compute the likelihood in the three-dimensional parameter space ($\Omega_m$, $\sigma_8$, $w$) from both observables. We find that constraints from peak counts are comparable to those from the power spectrum, and somewhat tighter when different smoothing scales are combined. Neither observable can constrain $w$ without external data. When the power spectrum and peak counts are combined, the area of the error "banana'' in the ($\Omega_m$, $\sigma_8$) plane reduces by a factor of $\approx2$, compared to using the power spectrum alone. For a flat $\Lambda$ cold dark matter model, combining both statistics, we obtain the constraint $\sigma_8(\Omega_m/0.27)^{0.63}=0.85\substack{+0.03 \\ -0.03}$.

We present a systematic exploration of dark energy and modified gravity models containing a single scalar field non-minimally coupled to the metric. Even though the parameter space is large, by exploiting an effective field theory (EFT) formulation and by imposing simple physical constraints such as stability conditions and (sub-)luminal propagation of perturbations, we arrive at a number of generic predictions. (1) The linear growth rate of matter density fluctuations is generally suppressed compared to $\Lambda$CDM at intermediate redshifts ($0.5 \lesssim z \lesssim 1$), despite the introduction of an attractive long-range scalar force. This is due to the fact that, in self-accelerating models, the background gravitational coupling weakens at intermediate redshifts, over-compensating the effect of the attractive scalar force. (2) At higher redshifts, the opposite happens; we identify a period of super-growth when the linear growth rate is larger than that predicted by $\Lambda$CDM. (3) The gravitational slip parameter $\eta$ - the ratio of the space part of the metric perturbation to the time part - is bounded from above. For Brans-Dicke-type theories $\eta$ is at most unity. For more general theories, $\eta$ can exceed unity at intermediate redshifts, but not more than about $1.5$ if, at the same time, the linear growth rate is to be compatible with current observational constraints. We caution against phenomenological parametrization of data that do not correspond to predictions from viable physical theories. We advocate the EFT approach as a way to constrain new physics from future large-scale-structure data.

Weak gravitational lensing is a powerful cosmological probe, with non--Gaussian features potentially containing the majority of the information. We examine constraints on the parameter triplet $(\Omega_m,w,\sigma_8)$ from non-Gaussian features of the weak lensing convergence field, including a set of moments (up to $4^{\rm th}$ order) and Minkowski functionals, using publicly available data from the 154deg$^2$ CFHTLenS survey. We utilize a suite of ray--tracing N-body simulations spanning 91 points in $(\Omega_m,w,\sigma_8)$ parameter space, replicating the galaxy sky positions, redshifts and shape noise in the CFHTLenS catalogs. We then build an emulator that interpolates the simulated descriptors as a function of $(\Omega_m,w,\sigma_8)$, and use it to compute the likelihood function and parameter constraints. We employ a principal component analysis to reduce dimensionality and to help stabilize the constraints with respect to the number of bins used to construct each statistic. Using the full set of statistics, we find $\Sigma_8\equiv\sigma_8(\Omega_m/0.27)^{0.55}=0.75\pm0.04$ (68% C.L.), in agreement with previous values. We find that constraints on the $(\Omega_m,\sigma_8)$ doublet from the Minkowski functionals suffer a strong bias. However, high-order moments break the $(\Omega_m,\sigma_8)$ degeneracy and provide a tight constraint on these parameters with no apparent bias. The main contribution comes from quartic moments of derivatives.

Consistency relations, which relate the squeezed limit of an (N+1)-point correlation function to an N-point function, are non-perturbative symmetry statements that hold even if the associated high momentum modes are deep in the nonlinear regime and astrophysically complex. Recently, Kehagias & Riotto and Peloso & Pietroni discovered a consistency relation applicable to large scale structure. We show that this can be recast into a simple physical statement in Lagrangian space: that the squeezed correlation function (suitably normalized) vanishes. This holds regardless of whether the correlation observables are at the same time or not, and regardless of whether multiple-streaming is present. The simplicity of this statement suggests that an analytic understanding of large scale structure in the nonlinear regime may be particularly promising in Lagrangian space.

The recently derived consistency relations for Large Scale Structure do not hold if the Equivalence Principle (EP) is violated. We show it explicitly in a toy model with two fluids, one of which is coupled to a fifth force. We explore the constraints that galaxy surveys can set on EP violation looking at the squeezed limit of the 3-point function involving two populations of objects. We find that one can explore EP violations of order 10^{-3} - 10^{-4} on cosmological scales. Chameleon models are already very constrained by the requirement of screening within the Solar System and only a very tiny region of the parameter space can be explored with this method. We show that no violation of the consistency relations is expected in Galileon models.

Consistency relations -- which relate an N-point function to a squeezed (N+1)-point function -- are useful in large scale structure (LSS) because of their non-perturbative nature: they hold even if the N-point function is deep in the nonlinear regime, and even if they involve astrophysically messy galaxy observables. The non-perturbative nature of the consistency relations is guaranteed by the fact that they are symmetry statements, in which the velocity plays the role of the soft pion. In this paper, we address two issues: (1) how to derive the relations systematically using the residual coordinate freedom in the Newtonian gauge, and relate them to known results in $\zeta$-gauge (often used in studies of inflation); (2) under what conditions the consistency relations are violated. In the non-relativistic limit, our derivation reproduces the Newtonian consistency relation discovered by Kehagias \& Riotto and Peloso & Pietroni. More generally, there is an infinite set of consistency relations, as is known in $\zeta$-gauge. There is a one-to-one correspondence between symmetries in the two gauges; in particular, the Newtonian consistency relation follows from the dilation and special conformal symmetries in $\zeta$-gauge. We probe the robustness of the consistency relations by studying models of galaxy dynamics and biasing. We give a systematic list of conditions under which the consistency relations are violated; violations occur if the galaxy bias is non-local in an infrared divergent way. We emphasize the relevance of the adiabatic mode condition, as distinct from symmetry considerations. As a by-product of our investigation, we discuss a simple fluid Lagrangian for LSS.

We study the two-point cross-correlation function between two populations of galaxies: for instance a bright population and a faint population. We show that this cross-correlation is asymmetric under the exchange of the line-of-sight coordinate of the galaxies, i.e. that the correlation is different if the bright galaxy is in front of, or behind, the faint galaxy. We give an intuitive, quasi-Newtonian derivation of all the effects that contribute to such an asymmetry in large-scale structure: gravitational redshift, Doppler shift, lensing, light-cone, evolution and Alcock-Paczynski effects - interestingly, the gravitational redshift term is exactly canceled by some of the others, assuming geodesic motion. Most of these effects are captured by previous calculations of general relativistic corrections to the observed galaxy density fluctuation; the asymmetry arises from terms that are suppressed by the ratio H/k - H is the Hubble constant and k is the wavenumber - which are more readily observable than the terms suppressed by (H/k)^2. Some of the contributions to the asymmetry, however, arise from terms that are generally considered 'Newtonian' - the lensing and evolution - and thus represent a contaminant in the search for general relativistic corrections. We propose methods to disentangle these different contributions. A simple method reduces the contamination to a level of < 10% for redshifts z<1. We also clarify the relation to recent work on measuring gravitational redshifts by stacking clusters.

We show that a static spherically symmetric black hole, in a generic theory of gravity with generic matter fields, has a two-dimensional Lorentz symmetry.

We show that the correlation functions of any single-field cosmological model with constant growing-modes are constrained by an infinite number of novel consistency relations, which relate (N+1)-point correlation functions with a soft-momentum scalar or tensor mode to a symmetry transformation on N-point correlation functions of hard-momentum modes. We derive these consistency relations from Ward identities for an infinite tower of non-linearly realized global symmetries governing scalar and tensor perturbations. These symmetries can be labeled by an integer n. At each order n, the consistency relations constrain - completely for n=0,1, and partially for n>= 2 - the q^n behavior of the soft limits. The identities at n=0 recover Maldacena's original consistency relations for a soft scalar and tensor mode, n=1 gives the recently-discovered conformal consistency relations, and the identities for n>= 2 are new. As a check, we verify directly that the n=2 identity is satisfied by known correlation functions in slow-roll inflation.

We compare the efficiency of moments and Minkowski functionals (MFs) in constraining the subset of cosmological parameters (Omega_m,w,sigma_8) using simulated weak lensing convergence maps. We study an analytic perturbative expansion of the MFs in terms of the moments of the convergence field and of its spatial derivatives. We show that this perturbation series breaks down on smoothing scales below 5', while it shows a good degree of convergence on larger scales (15'). Most of the cosmological distinguishing power is lost when the maps are smoothed on these larger scales. We also show that, on scales comparable to 1', where the perturbation series does not converge, cosmological constraints obtained from the MFs are approximately 1.5-2 times better than the ones obtained from the first few moments of the convergence distribution --- provided that the latter include spatial information, either from moments of gradients, or by combining multiple smoothing scales. Including either a set of these moments or the MFs can significantly tighten constraints on cosmological parameters, compared to the conventional method of using the power spectrum alone.

Magnetar magnetospheres are believed to be strongly twisted, due to shearing of the stellar crust by internal magnetic stresses. We present time-dependent axisymmetric simulations showing in detail the evolution of relativistic force-free magnetospheres subjected to slow twisting through large angles. When the twist amplitude is small, the magnetosphere moves quasi-statically through a sequence of equilibria of increasing free energy. At some twist amplitude the magnetosphere becomes tearing-mode unstable to forming a resistive current sheet, initiating large-scale magnetic reconnection in which a significant fraction of the magnetic free energy can be dissipated. This "critical" twist angle is insensitive to the resistive length scale. Rapid shearing temporarily stabilizes the magnetosphere beyond the critical angle, allowing the magnetosphere of a rapidly differentially rotating star to store, and dissipate, more free energy. In addition to these effects, shearing the surface of a rotating star increases the spindown torque applied to the star. If shearing is much slower than rotation, the resulting spikes in spindown rate can occur on timescales anywhere from the long twisting timescale to the stellar spin period or shorter, depending both on the stellar shear distribution and the existing distribution of magnetospheric twists. A model in which energy is stored in the magnetosphere, and released by a magnetospheric instability, therefore predicts large changes in the measured spindown rate before SGR giant flares.

The weak lensing power spectrum is a powerful tool to probe cosmological parameters. Additionally, lensing peak counts contain cosmological information beyond the power spectrum. Both of these statistics can be affected by the preferential selection of source galaxies in patches of the sky with high magnification, as well as by the dilution in the source galaxy surface density in such regions. If not accounted for, these biases introduce systematic errors for cosmological measurements. Here we quantify these systematic errors, using convergence maps from a suite of ray-tracing N-body simulations. At the cut-off magnitude m of on-going and planned major weak lensing surveys, the logarithmic slope of the cumulative number counts s = dlog[n(>m)]/dlog(m) is in the range 0.1 < s < 0.5. At s = 0.2, expected in the I band for LSST, the inferred values of Omega_m, w and sigma_8 are biased by many sigma (where sigma denotes the marginalized error) and therefore the biases will need to be carefully modeled. We also find that the parameters are biased differently in the (Omega_m, w, sigma_8) parameter space when the power spectrum and when the peak counts are used. In particular, w derived from the power spectrum is less affected than w derived from peak counts, while the opposite is true for the best-constrained combination of [sigma_8 Omega_m^gamma] (with gamma=0.62 from the power spectrum and gamma = 0.48 from peak counts). This suggests that the combination of the power spectrum and peak counts can help mitigate the impact of magnification and size biases.

The discovery of cosmic acceleration has stimulated theorists to consider dark energy or modifications to Einstein's General Relativity as possible explanations. The last decade has seen advances in theories that go beyond smooth dark energy -- modified gravity and interactions of dark energy. While the theoretical terrain is being actively explored, the generic presence of fifth forces and dark sector couplings suggests a set of distinct observational signatures. This report focuses on observations that differ from the conventional probes that map the expansion history or large-scale structure. Examples of such novel probes are: detection of scalar fields via lab experiments, tests of modified gravity using stars and galaxies in the nearby universe, comparison of lensing and dynamical masses of galaxies and clusters, and the measurements of fundamental constants at high redshift. The observational expertise involved is very broad as it spans laboratory experiments, high resolution astronomical imaging and spectroscopy and radio observations. In the coming decade, searches for these effects have the potential for discovering fundamental new physics. We discuss how the searches can be carried out using experiments that are already under way or with modest adaptations of existing telescopes or planned experiments. The accompanying paper on the Growth of Cosmic Structure describes complementary tests of gravity with observations of large-scale structure.

A theory of the quasidilaton is an extension of massive gravity by a scalar field, nonlinearly realizing a certain new global symmetry of the Lagrangian. It has been shown that unlike pure massive gravity, this theory does admit homogeneous and isotropic spatially flat solutions. Among the latter, selfaccelerated solutions attract a special attention. Previous studies of perturbations, performed in the decoupling limit, revealed one healthy scalar mode, while the second relevant scalar was not captured in that limit. Here we study full cosmological perturbations above the simplest selfaccelerated background. We show that the fluctuations of a mixed state of the quasidilaton and the helicity-0 graviton necessarily have a negative kinetic term at short distances, making this background unphysical. In addition, these cosmologies exhibit an order one sensitivity to higher dimensional terms suppressed by an energy scale that is parametrically higher than the strong coupling scale of the quasidilaton effective theory: such terms include Galileons, Goldstone-like selfinteractions and derivatives of the quasidilaton coupled to curvature, none of which introduce extra Ostrogradsky states. As one consequence, cosmology at the Hubble distances for this particular class of solutions depends on an unknown extension of the quasidilaton below its strong coupling distance scale. We note that non-FRW solutions that are similar to those of pure massive gravity should not necessarily suffer from these pathologies.

We derive consistency relations for correlators of scalar cosmological perturbations which hold in the "squeezed limit" in which one or more of the external momenta become soft. Our results are formulated as relations between suitably defined one-particle irreducible N-point and (N-1)-point functions that follow from residual spatial conformal diffeomorphisms of the unitary gauge Lagrangian. As such, some of these relations are exact to all orders in perturbation theory, and do not rely on approximate deSitter invariance or other dynamical assumptions (e.g., properties of the operator product expansion or the behavior of modes at horizon crossing). The consistency relations apply model-independently to cosmological scenarios where the time evolution is driven by a single scalar field. Besides reproducing the known results for single-field inflation in the slow roll limit, we verify that our consistency relations hold more generally, for instance in ghost condensate models in flat space. We comment on possible extensions of our results to multi-field models.

Gravitational waves at suitable frequencies can resonantly interact with a binary system, inducing changes to its orbit. A stochastic gravitational-wave background causes the orbital elements of the binary to execute a classic random walk, with the variance of orbital elements growing with time. The lack of such a random walk in binaries that have been monitored with high precision over long time-scales can thus be used to place an upper bound on the gravitational-wave background. Using periastron time data from the Hulse-Taylor binary pulsar spanning ~30 years, we obtain a bound of h_c < 7.9*10^(-14) at ~10^(-4) Hz, where h_c is the strain amplitude per logarithmic frequency interval. Our constraint complements those from pulsar timing arrays, which probe much lower frequencies, and ground-based gravitational-wave observations, which probe much higher frequencies. Interesting sources in our frequency band, which overlaps the lower sensitive frequencies of proposed space-based observatories, include white-dwarf/supermassive black-hole binaries in the early/late stages of inspiral, and TeV scale preheating or phase transitions. The bound improves as (time span)^(-2) and (sampling rate)^(-1/2). The Hulse-Taylor constraint can be improved to ~3.8*10^(-15) with a suitable observational campaign over the next decade. Our approach can also be applied to other binaries, including (with suitable care) the Earth-Moon system, to obtain constraints at different frequencies. The observation of additional binary pulsars with the SKA could reach a sensitivity of h_c ~ 3*10^(-17).

Non-Gaussianity in the inflationary perturbations can couple observable scales to modes of much longer wavelength (even superhorizon), leaving as a signature a large-angle modulation of the observed cosmic microwave background (CMB) power spectrum. This provides an alternative origin for a power asymmetry which is otherwise often ascribed to a breaking of statistical isotropy. The non-Gaussian modulation effect can be significant even for typical ~10^{-5} perturbations, while respecting current constraints on non-Gaussianity, if the squeezed limit of the bispectrum is sufficiently infrared divergent. Just such a strongly infrared-divergent bispectrum has been claimed for inflation models with a non-Bunch-Davies initial state, for instance. Upper limits on the observed CMB power asymmetry place stringent constraints on the duration of inflation in such models.

General Relativity (GR), with or without matter fields, admits a natural extension to a scale invariant theory that requires a dilaton. Here we show that the recently formulated massive GR, minimally coupled to matter, possesses a new global symmetry related to scaling of the reference coordinates w.r.t. the physical ones. The field enforcing this symmetry, dubbed here quasi-dilaton, coincides with an ordinary dilaton if only pure gravity is considered, but differs from it when the matter Lagrangian is present. We study: (1) Theoretical consistency of massive GR with the quasi-dilaton; (2) Consistency with observations for spherically symmetric sources on (nearly) flat backgrounds; (3) Cosmological implications of this theory. We find that: (I) The theory with the quasi-dilaton is as consistent as massive GR is. (II) The Vainshtein mechanism is generically retained, owing to the fact that in the decoupling limit there is an enhanced symmetry, which turns the quasi-dilaton into a second galileon, consistently coupled to a tensor field. (III) Unlike in massive GR, there exist flat FRW solutions. In particular, we find self-accelerated solutions and discuss their quadratic perturbations. These solutions are testable by virtue of the different effective Newton's constants that govern the Hubble expansion and structure growth.

We present the first simulations of evolving, strongly twisted magnetar magnetospheres. Slow shearing of the magnetar crust is seen to lead to a series of magnetospheric expansion and reconnection events, corresponding to X-ray flares and bursts. The axisymmetric simulations include rotation of the neutron star and the magnetic wind through the light cylinder. We study how the increasing twist affects the spindown rate of the star, finding that a dramatic increase in spindown occurs. Particularly spectacular are explosive events caused by the sudden opening of large amounts of overtwisted magnetic flux, which may be associated with the observed giant flares. These events are accompanied by a short period of ultra-strong spindown, resulting in an abrupt increase in spin period, such as was observed in the giant flare of SGR 1900+14.

The chameleon, or generalizations thereof, is a light scalar that couple to matter with gravitational strength, but whose manifestation depends on the ambient matter density. A key feature is that the screening mechanism suppressing its effects in high-density environments is determined by the local scalar field value. Under very general conditions, we prove two theorems limiting its cosmological impact: i) the Compton wavelength of such a scalar can be at most Mpc at present cosmic density, which restricts its impact to non-linear scales; ii) the conformal factor relating Einstein- and Jordan-frame scale factors is essentially constant over the last Hubble time, which precludes the possibility of self-acceleration. These results imply that chameleon-like scalar fields have a negligible effect on the linear-scale growth history; theories that invoke a chameleon-like scalar to explain cosmic acceleration rely on a form of dark energy rather than a genuine modified gravity effect. Our analysis applies to a broad class of chameleon, symmetron and dilaton theories.

We describe a new scheme for evolving the equations of force-free electrodynamics, the vanishing-inertia limit of magnetohydrodynamics. This pseudospectral code uses global orthogonal basis function expansions to take accurate spatial derivatives, allowing the use of an unstaggered mesh and the complete force-free current density. The method has low numerical dissipation and diffusion outside of singular current sheets. We present a range of one- and two-dimensional tests, and demonstrate convergence to both smooth and discontinuous analytic solutions. As a first application, we revisit the aligned rotator problem, obtaining a steady solution with resistivity localised in the equatorial current sheet outside the light cylinder.

We remark on the existence of non-linearly realized conformal symmetries for scalar adiabatic perturbations in cosmology. These conformal symmetries are present for any cosmological background, beyond any slow-roll or quasi-de Sitter approximation. The dilatation transformation shifts the curvature perturbation by a constant, and corresponds to the well-known symmetry under spatial rescaling. We argue that the scalar sector is also invariant under special conformal transformations, which shift the curvature perturbation by a term linear in the spatial coordinates. We discuss whether these conformal symmetries can be extended to include tensor perturbations. Tensor modes introduce their own set of non-linearly realized symmetries. We identify an infinite set of large gauge transformations which maintain the transverse, traceless gauge condition, while shifting the tensor mode non-trivially.

We study the scale-dependence of halo bias in generic (non-local) primordial non-Gaussian (PNG) initial conditions of the type motivated by inflation, parametrized by an arbitrary quadratic kernel. We first show how to generate non-local PNG initial conditions with minimal overhead compared to local PNG models for a general class of primordial bispectra that can be written as linear combinations of separable templates. We run cosmological simulations for the local, and non-local equilateral and orthogonal models and present results on the scale-dependence of halo bias. We also derive a general formula for the Fourier-space bias using the peak-background split (PBS) in the context of the excursion set approach to halos and discuss the difference and similarities with the known corresponding result from local bias models. Our PBS bias formula generalizes previous results in the literature to include non-Markovian effects and non-universality of the mass function and are in better agreement with measurements in numerical simulations than previous results for a variety of halo masses, redshifts and halo definitions. We also derive for the first time quadratic bias results for arbitrary non-local PNG, and show that non-linear bias loops give small corrections at large-scales. The resulting well-behaved perturbation theory paves the way to constrain non-local PNG from measurements of the power spectrum and bispectrum in galaxy redshift surveys.

We consider a galileon field coupled to gravity. The standard no-hair theorems do not apply, because of the galileon's peculiar derivative interactions. We prove that, nonetheless, static spherically symmetric black holes cannot sustain non-trivial galileon profiles. Our theorem holds regardless of whether there are non-minimal couplings between the galileon and gravity of the covariant galileon type.

Modified gravity theories capable of genuine self-acceleration typically invoke a galileon scalar which mediates a long range force, but is screened by the Vainshtein mechanism on small scales. In such theories, non-relativistic stars carry the full scalar charge (proportional to their mass), while black holes carry none. Thus, for a galaxy free-falling in some external gravitational field, its central massive black hole is expected to lag behind the stars. To look for this effect, and to distinguish it from other astrophysical effects, one can correlate the gravitational pull from the surrounding structure with the offset between the stellar center and the black hole. The expected offset depends on the central density of the galaxy, and ranges up to ~0.1 kpc for small galaxies. The observed offset in M87 cannot be explained by this effect unless the scalar force is significantly stronger than gravity. We also discuss the systematic offset of compact objects from the galactic plane as another possible signature.

We present an analysis of peculiar velocities and their effect on supernova cosmology. In particular, we study (a) the corrections due to our own motion, (b) the effects of correlations in peculiar velocities induced by large-scale structure, and (c) uncertainties arising from a possible local under- or over-density. For all of these effects we present a case study of their impact on the cosmology derived by the Sloan Digital Sky Survey-II Supernova Survey (SDSS-II SN Survey). Correcting supernova redshifts for the CMB dipole slightly over-corrects nearby supernovae that share some of our local motion. We show that while neglecting the CMB dipole would cause a shift in the derived equation of state of Delta w ~ 0.04 (at fixed matter density) the additional local-motion correction is currently negligible (Delta w<0.01). We use a covariance-matrix approach to statistically account for correlated peculiar velocities. This down-weights nearby supernovae and effectively acts as a graduated version of the usual sharp low-redshift cut. Neglecting coherent velocities in the current sample causes a systematic shift of ~2% in the preferred value of w and will therefore have to be considered carefully when future surveys aim for percent-level accuracy. Finally, we perform n-body simulations to estimate the likely magnitude of any local density fluctuation (monopole) and estimate the impact as a function of the low-redshift cutoff. We see that for this aspect the low-z cutoff of z=0.02 is well-justified theoretically, but that living in a putative local density fluctuation leaves an indelible imprint on the magnitude-redshift relation.

We consider the low-energy effective field theory describing the infrared dynamics of non-dissipative fluids. We extend previous work to accommodate conserved charges, and we clarify the matching between field theory variables and thermodynamical ones. We discuss the systematics of the derivative expansion, for which field theory offers a conceptually clear and technically neat scheme. As an example, we compute the correction to the sound-wave dispersion relation coming from a sample second-order term. This formalism forms the basis for a study of anomalies in hydrodynamics via effective field theory, which is initiated in a companion paper.

We develop the formalism that incorporates quantum anomalies in the effective field theory of non-dissipative fluids. We consider the effect of adding a Wess-Zumino-like term to the low-energy effective action to account for anomalies. In this paper we restrict to two spacetime dimensions. We find modifications to the constitutive relations for the current and the stress-energy tensor, and, more interestingly, half a new propagating mode (one-and-a-halfth sound): a left- or right-moving wave with propagation speed that goes to zero with the anomaly coefficient. Unlike for the chiral magnetic wave in four dimensions, this mode propagates even in the absence of external fields. We check our results against a more standard, purely hydrodynamical derivation. Unitarity of the effective field theory suggests an upper bound on the anomaly coefficient in hydrodynamics.

Similar to the magnification of the galaxies' fluxes by gravitational lensing, the extinction of the fluxes by comic dust, whose existence is recently detected by Menard et al (2009), also modify the distribution of a flux-selected galaxy sample. We study the anisotropic distortion by dust extinction to the 3D galaxy correlation function, including magnification bias and redshift distortion at the same time. We find the extinction distortion is most significant along the line of sight and at large separations, similar to that by magnification bias. The correction from dust extinction is negative except at sufficiently large transverse separations, which is almost always opposite to that from magnification bias (we consider a number count slope s > 0.4). Hence, the distortions from these two effects tend to reduce each other. At low z (~<1), the distortion by extinction is stronger than that by magnification bias, but at high z, the reverse holds. We also study how dust extinction affects probes in real space of the baryon acoustic oscillations (BAO) and the linear redshift distortion parameter beta. We find its effect on BAO is negligible. However, it introduces a positive scale-dependent correction to beta that can be as large as a few percent. At the same time, we also find a negative scale-dependent correction from magnification bias, which is up to percent level at low z, but to ~40% at high z. These corrections are non-negligible for precision cosmology, and should be considered when testing General Relativity through the scale-dependence of beta.

We derive a general expression for the large-scale halo bias, in theories with a scale-dependent linear growth, using the excursion set formalism. Such theories include modified gravity models, and models in which the dark energy clustering is non-negligible. A scale dependence is imprinted in both the formation and evolved biases by the scale-dependent growth. Mergers are accounted for in our derivation, which thus extends earlier work which focused on passive evolution. There is a simple analytic form for the bias for those theories in which the nonlinear collapse of perturbations is approximately the same as in general relativity. As an illustration, we apply our results to a simple Yukawa modification of gravity, and use SDSS measurements of the clustering of luminous red galaxies to constrain the theory's parameters.

Theories that attempt to explain cosmic acceleration by modifying gravity typically introduces a long-range scalar force that needs to be screened on small scales. One common screening mechanism is the chameleon, where the scalar force is screened in environments with a sufficiently deep gravitational potential, but acts unimpeded in regions with a shallow gravitational potential. This leads to a variation in the overall gravitational G with environment. We show such a variation can occur within a star itself, significantly affecting its evolution and structure, provided that the host galaxy is unscreened. The effect is most pronounced for red giants, which would be smaller by a factor of tens of percent and thus hotter by 100's of K, depending on the parameters of the underlying scalar-tensor theory. Careful measurements of these stars in suitable environments (nearby dwarf galaxies not associated with groups or clusters) would provide constraints on the chameleon mechanism that are four orders of magnitude better than current large scale structure limits, and two orders of magnitude better than present solar system tests.

In Lyman-alpha forest measurements it is generally assumed that quasars are mere background light sources which are uncorrelated with the forest. Gravitational lensing of the quasars violates this assumption. This effect leads to a measurement bias, but more interestingly it provides a valuable signal. The lensing signal can be extracted by correlating quasar magnitudes with the flux power spectrum and with the flux decrement. These correlations will be challenging to measure but their detection provides a direct measure of how features in the Lyman-alpha forest trace the underlying mass density field. Observing them will test the fundamental hypothesis that fluctuations in the forest are predominantly driven by fluctuations in mass, rather than in the ionizing background, helium reionization or winds. We discuss ways to disentangle the lensing signal from other sources of such correlations, including dust, continuum and background residuals. The lensing-induced measurement bias arises from sample selection: one preferentially collects spectra of magnified quasars which are behind overdense regions. This measurement bias is ~0.1-1% for the flux power spectrum, optical depth and the flux probability distribution. Since the effect is systematic, quantities such as the amplitude of the flux power spectrum averaged across scales should be interpreted with care.

The equivalence of inertial and gravitational masses is a defining feature of general relativity. Here, we clarify the status of the equivalence principle for interactions mediated by a universally coupled scalar, motivated partly by recent attempts to modify gravity at cosmological distances. Although a universal scalar-matter coupling is not mandatory, once postulated, it is stable against classical and quantum renormalizations in the matter sector. The coupling strength itself is subject to renormalization of course. The scalar equivalence principle is violated only for objects for which either the graviton self-interaction or the scalar self-interaction is important---the first applies to black holes, while the second type of violation is avoided if the scalar is Galilean-symmetric.

It has recently been shown in high resolution numerical simulations that relativistic collisions of bubbles in the context of a multi-vacua potential may lead to the creation of bubbles in a new vacuum. In this paper, we show that scalar fields with only potential interactions behave like free fields during high-speed collisions; the kick received by them in a collision can be deduced simply by a linear superposition of the bubble wall profiles. This process is equivalent to the scattering of solitons in 1+1 dimensions. We deduce an expression for the field excursion (shortly after a collision), which is related simply to the field difference between the parent and bubble vacua, i.e. contrary to expectations, the excursion cannot be made arbitrarily large by raising the collision energy. There is however a minimum energy threshold for this excursion to be realized. We verify these predictions using a number of 3+1 and 1+1 numerical simulations. A rich phenomenology follows from these collision induced excursions - they provide a new mechanism for scanning the landscape, they might end/begin inflation, and they might constitute our very own big bang, leaving behind a potentially observable anisotropy.

Scalar field theories with derivative interactions are known to possess solitonic excitations, but such solitons are generally unsatisfactory because the effective theory fails precisely where nonlinearities responsible for the solitons are important. A new class of theories possessing (internal) galilean invariance can in principle bypass this difficulty. Here, we show that these galileon theories do not possess stable solitonic solutions. As a by-product, we show that no stable solitons exist for a different class of derivatively coupled theories, describing for instance the infrared dynamics of superfluids, fluids, solids and some k-essence models.

Theories that attempt to explain the observed cosmic acceleration by modifying general relativity all introduce a new scalar degree of freedom that is active on large scales, but is screened on small scales to match experiments. We show that if such screening occurrs via the chameleon mechanism such as in f(R), it is possible to have order one violation of the equivalence principle, despite the absence of explicit violation in the microscopic action. Namely, extended objects such as galaxies or constituents thereof do not all fall at the same rate. The chameleon mechanism can screen the scalar charge for large objects but not for small ones (large/small is defined by the gravitational potential and controlled by the scalar coupling). This leads to order one fluctuations in the inertial to gravitational mass ratio. In Jordan frame, it is no longer true that all objects move on geodesics. In contrast, if the scalar screening occurrs via strong coupling, such as in the DGP braneworld model, equivalence principle violation occurrs at a much reduced level. We propose several observational tests of the chameleon mechanism: 1. small galaxies should fall faster than large galaxies, even when dynamical friction is negligible; 2. voids defined by small galaxies would be larger compared to standard expectations; 3. stars and diffuse gas in small galaxies should have different velocities, even on the same orbits; 4. lensing and dynamical mass estimates should agree for large galaxies but disagree for small ones. We discuss possible pitfalls in some of these tests. The cleanest is the third one where mass estimate from HI rotational velocity could exceed that from stars by 30 % or more. To avoid blanket screening of all objects, the most promising place to look is in voids.

Only galaxies bright enough and large enough to be unambiguously identified and measured are included in galaxy surveys used to estimate cosmic shear. We demonstrate that because gravitational lensing can scatter galaxies across the brightness and size thresholds, cosmic shear experiments suffer from lensing bias. We calculate the effect on the shear power spectrum and show that - unless corrected for - it will lead analysts to cosmological parameters estimates that are biased at the 2-3\sigma level in DETF Stage III experiments, such as the Dark Energy Survey.

Only certain galaxies are included in surveys: those bright and large enough to be detectable as extended sources. Because gravitational lensing can make galaxies appear both brighter and larger, the presence of foreground inhomogeneities can scatter galaxies across not only magnitude cuts but also size cuts, changing the statistical properties of the resulting catalog. Here we explore this size bias, and how it combines with magnification bias to affect galaxy statistics. We demonstrate that photometric galaxy samples from current and upcoming surveys can be even more affected by size bias than by magnification bias.

Given a scalar field with metastable minima, bubbles nucleate quantum mechanically. When bubbles collide, energy stored in the bubble walls is converted into kinetic energy of the field. This kinetic energy can facilitate the classical nucleation of new bubbles in minima that lie below those of the "parent" bubbles. This process is efficient and classical, and changes the dynamics and statistics of bubble formation in models with multiple vacua, relative to that derived from quantum tunneling.

We study the clustering of LRG galaxies in the latest spectroscopic SDSS data releases, DR6 and DR7, which sample over 1 Gpc^3/h^3 to z=0.47. The 2-point correlation function $\xisp$ is estimated as a function of perpendicular $\sigma$ and line-of-sight $\pi$ (radial) directions. We find a significant detection of a peak at $r\simeq 110$Mpc/h, which shows as a circular ring in the $\sigma-\pi$ plane. There is also significant evidence for a peak along the radial direction whose shape is consistent with its originating from the recombination-epoch baryon acoustic oscillations (BAO). A $\xisp$ model with no radial BAO peak is disfavored at $3.2\sigma$, whereas a model with no magnification bias is disfavored at $2\sigma$. The radial data enable, for the first time, a direct measurement of the Hubble parameter $H(z)$ as a function of redshift. This is independent from earlier BAO measurements which used the spherically averaged (monopole) correlation to constrain an integral of $H(z)$. Using the BAO peak position as a standard ruler in the radial direction, we find: $H(z=0.24)= 79.69 \pm 2.32 (\pm 1.29)$ km/s/Mpc for z=0.15-0.30 and $H(z=0.43)= 86.45 \pm 3.27 (\pm 1.69)$ km/s/Mpc for $z=0.40-0.47$. The first error is a model independent statistical estimation and the second accounts for systematics both in the measurements and in the model. For the full sample, $z=0.15-0.47$, we find $H(z=0.34)= 83.80 \pm 2.96 (\pm 1.59)$ km/s/Mpc.

Galaxy-galaxy or galaxy-quasar lensing can provide important information on the mass distribution in the Universe. It consists of correlating the lensing signal (either shear or magnification) of a background galaxy/quasar sample with the number density of a foreground galaxy sample. However, the foreground galaxy density is inevitably altered by the magnification bias due to the mass between the foreground and the observer, leading to a correction to the observed galaxy-lensing signal. The aim of this paper is to quantify this correction. The single most important determining factor is the foreground redshift z: the correction is small if the foreground galaxies are at low redshifts but can become non-negligible for sufficiently high redshifts. For instance, we find that for the multipole l=1000, the correction is above 1%*(5s-2)/b for z<0.37, and above 5%*(5s-2)/b for z<0.67, where s is the number count slope of the foreground sample, and b its galaxy bias. These considerations are particularly important for geometrical measures, such as the Jain and Taylor ratio or its generalization by Zhang et al. Assuming (5s-2)/b=1, we find that the foreground redshift should be limited to z<0.45 in order to avoid biasing the inferred dark energy equation of state w by more than 5%, and that even for a low foreground redshift (< 0.45), the background samples must be well separated from the foreground to avoid incurring a bias of similar magnitude. Lastly, we briefly comment on the possibility of obtaining these geometrical measures without using galaxy shapes, using instead magnification bias itself.

We conduct a Markov Chain Monte Carlo study of the Dvali-Gabadadze-Porrati (DGP) self-accelerating braneworld scenario given the cosmic microwave background (CMB) anisotropy, supernovae and Hubble constant data by implementing an effective dark energy prescription for modified gravity into a standard Einstein-Boltzmann code. We find no way to alleviate the tension between distance measures and horizon scale growth in this model. Growth alterations due to perturbations propagating into the bulk appear as excess CMB anisotropy at the lowest multipoles. In a flat cosmology, the maximum likelihood DGP model is nominally a 5.3 sigma poorer fit than Lambda CDM. Curvature can reduce the tension between distance measures but only at the expense of exacerbating the problem with growth leading to a 4.8 sigma result that is dominated by the low multipole CMB temperature spectrum. While changing the initial conditions to reduce large scale power can flatten the temperature spectrum, this also suppresses the large angle polarization spectrum in violation of recent results from WMAP5. The failure of this model highlights the power of combining growth and distance measures in cosmology as a test of gravity on the largest scales.

We discuss models that can account for today's dark energy. The underlying cosmological constant may be Planck scale but starts as a redundant coupling which can be eliminated by a field redefinition. The observed vacuum energy arises when the redundancy is explicitly broken, say by a non-minimal coupling to curvature. We give a recipe for constructing models, including R + 1/R type models, that realize this mechanism and satisfy all solar system constraints on gravity. A similar model, based on Gauss-Bonnet gravity, provides a technically natural explanation for dark energy and exhibits an interesting see-saw behavior: a large underlying cosmological constant gives rise to both low and high curvature solutions. Such models could be statistically favored in the string landscape.

In paper I of this series we discuss how magnification bias distorts the 3D correlation function by enhancing the observed correlation in the line-of-sight (LOS) orientation, especially on large scales. This lensing anisotropy is distinctive, making it possible to separately measure the galaxy-galaxy, galaxy-magnification {\it and} magnification-magnification correlations. Here we extend the discussion to the power spectrum and also to redshift space. In real space, pairs oriented close to the LOS direction are not protected against nonlinearity even if the pair separation is large; this is because nonlinear fluctuations can enter through gravitational lensing at a small transverse separation (or i.e. impact parameter). The situation in Fourier space is different: by focusing on a small wavenumber $k$, as is usually done, linearity is guaranteed because both the LOS and transverse wavenumbers must be small. This is why magnification distortion of the galaxy correlation appears less severe in Fourier space. Nonetheless, the effect is non-negligible, especially for the transverse Fourier modes, and should be taken into account in interpreting precision measurements of the galaxy power spectrum, for instance those that focus on the baryon oscillations. The lensing induced anisotropy of the power spectrum has a shape that is distinct from the more well known redshift space anisotropies due to peculiar motions and the Alcock-Paczynski effect. The lensing anisotropy is highly localized in Fourier space while redshift space distortions are more spread out. This means that one could separate the magnification bias component in real observations, implying that potentially it is possible to perform a gravitational lensing measurement without measuring galaxy shapes.

Fry (1996) showed that galaxy bias has the tendency to evolve towards unity, i.e. in the long run, the galaxy distribution tends to trace that of matter. Generalizing slightly Fry's reasoning, we show that his conclusion remains valid in theories of modified gravity (or equivalently, complex clustered dark energy). This is not surprising: as long as both galaxies and matter are subject to the same force, dynamics would drive them towards tracing each other. This holds, for instance, in theories where both galaxies and matter move on geodesics. This relaxation of bias towards unity is tempered by cosmic acceleration, however: the bias tends towards unity but does not quite make it, unless the formation bias were close to unity. Our argument is extended in a straightforward manner to the case of a stochastic or nonlinear bias. An important corollary is that dynamical evolution could imprint a scale dependence on the large scale galaxy bias. This is especially pronounced if non-standard gravity introduces new scales to the problem: the bias at different scales relaxes at different rates, the larger scales generally more slowly and retaining a longer memory of the initial bias. A consistency test of the current (general relativity + uniform dark energy) paradigm is therefore to look for departure from a scale independent bias on large scales. A simple way is to measure the relative bias of different populations of galaxies which are at different stages of bias relaxation. Lastly, we comment on the possibility of directly testing the Poisson equation on cosmological scales, as opposed to indirectly through the growth factor.

We develop an analytic approach to study inhomogeneous reionization on large scales by solving the equations of ionization balance and radiative transfer to first order in perturbations. Given the spatial distribution and spectrum of the ionizing sources, our formalism can be used to predict the large scale power spectra of fluctuations in the abundances of HII, HI and radiation. Our approach avoids common approximations/assumptions in existing analytic methods -- for instance, we do not assume a specific ionization topology from the outset; nor do we make a step-function bubble-like approximation to the HII distribution. Applying our formalism to sources biased according to the Press-Schechter prescription, we find: 1. reionization always proceeds ``inside--out'', with dense regions more highly ionized, at least on large scales; 2. on sufficiently large scales, HII, HI and radiation exhibits a scale independent bias relative to dark matter; 3. the bias is suppressed on scales comparable to or smaller than the mean free path of the ionizing photons; 4. if the ionizing source spectrum is sufficiently soft, the HII bias closely tracks the source bias for most of the reionization process but drops precipitously after percolation; 5. if the ionizing source spectrum is sufficiently hard, the HII bias drops in a more steady fashion throughout the reionization process. The tools developed here will be useful for interpreting future 21 cm, CMB and Lyman-alpha forest observations, both to learn about the reionization astrophysics (such as the hardness of the source spectrum and therefore the nature of the ionizing sources) and to possibly extract interesting cosmological information.

It is well known that magnification bias, the modulation of galaxy or quasar source counts by gravitational lensing, can change the observed angular correlation function. We investigate magnification-induced changes to the shape of the observed correlation function w(\theta) and the angular power spectrum C_{\ell}, paying special attention to the matter-radiation equality peak and the baryon wiggles. Lensing mixes the correlation function of the source galaxies with the matter correlation at the lower redshifts of the lenses. Since the lenses probe structure nearer to the observer, the angular scale dependence of the lensing terms is different from that of the sources, thus the observed correlation function is distorted. We quantify how the lensing corrections depend on the width of the selection function, the galaxy bias b, and the number count slope s. The correction increases with redshift and larger corrections are present for sources with steep number count slopes and/or broad redshift distributions. The most drastic changes to C_{\ell} occur for measurements at z >~1.5 and \ell <~ 100. For the source distributions we consider, magnification bias can shift the matter-radiation equality scale by 1-6% at z ~ 1.5 and by z ~ 3.5 the shift can be as large as 30%. The baryon bump in \theta^2w(\theta) is shifted by <~ 1% and the width is typically increased by ~10%. Shifts of >~ 0.5% and broadening of >~ 20% occur only for very broad selection functions and/or galaxies with (5s-2)/b>~2. However, near the baryon bump the magnification correction is not constant but a gently varying function which depends on the source population. Depending on how the w(\theta) data is fitted, this correction may need to be accounted for when using the baryon acoustic scale for precision cosmology.

We apply the technique of parameter-splitting to existing cosmological data sets, to check for a generic failure of dark energy models. Given a dark energy parameter, such as the energy density Omega_Lambda or equation of state w, we split it into two meta-parameters with one controlling geometrical distances, and the other controlling the growth of structure. Observational data spanning Type Ia Supernovae, the cosmic microwave background (CMB), galaxy clustering, and weak gravitational lensing statistics are fit without requiring the two meta-parameters to be equal. This technique checks for inconsistency between different data sets, as well as for internal inconsistency within any one data set (e.g., CMB or lensing statistics) that is sensitive to both geometry and growth. We find that the cosmological constant model is consistent with current data. Theories of modified gravity generally predict a relation between growth and geometry that is different from that of general relativity. Parameter-splitting can be viewed as a crude way to parametrize the space of such theories. Our analysis of current data already appears to put sharp limits on these theories: assuming a flat universe, current data constrain the difference Omega_Lambda(geom) - Omega_Lambda(grow) to be -0.0044 +/- 0.0058 (68% C.L.); allowing the equation of state w to vary, the difference w(geom) - w(grow) is constrained to be 0.37 +/- 0.37 (68% C.L.). Interestingly, the region w(grow) > w(geom), which should be generically favored by theories that slow structure formation relative to general relativity, is quite restricted by data already. We find w(grow) < -0.80 at 2 sigma. As an example, the best-fit flat Dvali-Gabadadze-Porrati (DGP) model approximated by our parametrization lies beyond the 3 sigma contour for constraints from all the data sets.

It is well known gravitational lensing, mainly via magnification bias, modifies the observed galaxy/quasar clustering. Such discussions have largely focused on the 2D angular correlation. Here and in a companion paper (Paper II) we explore how magnification bias distorts the 3D correlation function and power spectrum, as first considered by Matsubara. The interesting point is: the distortion is anisotropic. Magnification bias preferentially enhances the observed correlation in the line-of-sight (LOS) orientation, especially on large scales. For example at LOS separation of ~100 Mpc/h, where the intrinsic galaxy-galaxy correlation is rather weak, the observed correlation can be enhanced by lensing by a factor of a few, even at a modest redshift of z ~ 0.35. The opportunity: this lensing anisotropy is distinctive, making it possible to separately measure the galaxy-galaxy, galaxy-magnification and magnification-magnification correlations, without measuring galaxy shapes. The anisotropy is distinguishable from the well known distortion due to peculiar motions, as will be discussed in Paper II. The challenge: the magnification distortion of the galaxy correlation must be accounted for in interpreting data as precision improves. For instance, the ~100 Mpc/h baryon acoustic oscillation scale in the correlation function is shifted by up to ~3% in the LOS orientation, and up to ~0.6% in the monopole, depending on the galaxy bias, redshift and number count slope. The corresponding shifts in the inferred Hubble parameter and angular diameter distance, if ignored, could significantly bias measurements of the dark energy equation of state. Lastly, magnification distortion offers a plausible explanation for the well known excess correlations seen in pencil beam surveys.

Integrated Sachs-Wolfe (ISW) measurements, which involve cross-correlating the CMB with the foreground large-scale structure (e.g. galaxies/quasars), have proven to be an interesting probe of dark energy. We show that magnification bias, which is the inevitable modulation of the foreground number counts by gravitational lensing, alters both the shape and amplitude of the observed ISW signal. This is true especially at high redshifts because (1) the intrinsic galaxy-temperature signal diminishes greatly back in the matter dominated era, (2) the lensing efficiency increases with redshift and (3) the number count slope generally steepens with redshift in a magnitude limited sample. At z >~ 2, the magnification-temperature correlation dominates over the intrinsic galaxy-temperature correlation and causes the observed ISW signal to increase with z, despite dark energy subdominance -- a result of the fact that magnification probes structures between the observer and the sources. Ignoring magnification bias can then lead to erroneous conclusions about dark energy. While the lensing modulation opens up an interesting high z window for ISW measurements, high z measurements are not expected to add much new information to low z ones if dark energy is the cosmological constant. This is because lensing introduces significant covariance across redshifts. The most compelling reason to pursue high z ISW measurements is to look for a potential surprise such as early dark energy domination or the signature of modified gravity. We conclude with a discussion of existing measurements, the highest z of which is at the margin of being sensitive to magnification bias. We also develop a formalism which might be of general interest: to predict biases in estimating parameters when certain physical effects are ignored in interpreting data.

Large scale structure introduces two different kinds of errors in the luminosity distance estimates from standardizable candles such as supernovae Ia (SNe) - a Poissonian scatter for each SN and a coherent component due to correlated fluctuations between different SNe. Increasing the number of SNe helps reduce the first type of error but not the second. The coherent component has been largely ignored in forecasts of dark energy parameter estimation from upcoming SN surveys. For instance it is commonly thought, based on Poissonian considerations, that peculiar motion is unimportant, even for a low redshift SN survey such as the Nearby Supernova Factory (SNfactory; z = 0.03 - 0.08), which provides a useful anchor for future high redshift surveys by determining the SN zero-point. We show that ignoring coherent peculiar motion leads to an underestimate of the zero-point error by about a factor of 2, despite the fact that SNfactory covers almost half of the sky. More generally, there are four types of fluctuations: peculiar motion, gravitational lensing, gravitational redshift and what is akin to the integrated Sachs-Wolfe effect. Peculiar motion and lensing dominates at low and high redshifts respectively. Taking into account all significant luminosity distance fluctuations due to large scale structure leads to a degradation of up to 60% in the determination of the dark energy equation of state from upcoming high redshift SN surveys, when used in conjunction with a low redshift anchor such as the SNfactory. The most relevant fluctuations are the coherent ones due to peculiar motion and the Poissonian ones due to lensing, with peculiar motion playing the dominant role. We also discuss to what extent the noise here can be viewed as a useful signal, and whether corrections can be made to reduce the degradation.

Random walks with a general, nonlinear barrier have found recent applications ranging from reionization topology to refinements in the excursion set theory of halos. Here, we derive the first-crossing distribution of random walks with a moving barrier of an arbitrary shape. Such a distribution is shown to satisfy an integral equation that can be solved by a simple matrix inversion, without the need for Monte Carlo simulations, making this useful for exploring a large parameter space. We discuss examples in which common analytic approximations fail, a failure which can be remedied using the method described here.

The analysis of the Lyman-alpha forest of absorption lines in quasar spectra has emerged as a potentially powerful technique to constrain the linear matter power spectrum. In most previous work, the amplitude of the ionizing background was fixed by calibrating simulations to match the observed mean transmitted flux in the Lyman-alpha forest. This procedure is undesirable in principle as it requires the estimation of the unabsorbed quasar continuum level, a difficult undertaking subject to various sources of systematic error and bias. We suggest an alternative approach based on measuring the one-point probability distribution function (pdf) of the fluctuations in the flux about the mean, relative to the mean, i.e. the pdf of delta_f = (f - <f>)/<f>. This statistic, while sensitive to the amplitude of the ionizing background, has the virtue that its measurement does not require an estimate of the unabsorbed continuum level. We present a measurement of the pdf of delta_f from seven Keck HIRES spectra, spanning a redshift range of z = 2.2 - 4.4. To illustrate that our method is useful, we compare our measurements of the pdf of delta_f, and measurements of the flux power spectrum from Croft et al. (2002) at z = 2.72, with cosmological simulations. From this comparison, we obtain constraints on the mean transmission in the Lyman-alpha forest, the slope of the temperature-density relation, as well as the amplitude and slope of the mass power spectrum. Our methodology will be useful for obtaining more precise constraints with larger data samples from the Sloan Digital Sky Survey (SDSS).

The Supernova / Acceleration Probe (SNAP) is a proposed space-based experiment designed to study the dark energy and alternative explanations of the acceleration of the Universe's expansion by performing a series of complementary systematics-controlled measurements. We describe a self-consistent reference mission design for building a Type Ia supernova Hubble diagram and for performing a wide-area weak gravitational lensing study. A 2-m wide-field telescope feeds a focal plane consisting of a 0.7 square-degree imager tiled with equal areas of optical CCDs and near infrared sensors, and a high-efficiency low-resolution integral field spectrograph. The SNAP mission will obtain high-signal-to-noise calibrated light-curves and spectra for several thousand supernovae at redshifts between z=0.1 and 1.7. A wide-field survey covering one thousand square degrees resolves ~100 galaxies per square arcminute. If we assume we live in a cosmological-constant-dominated Universe, the matter density, dark energy density, and flatness of space can all be measured with SNAP supernova and weak-lensing measurements to a systematics-limited accuracy of 1%. For a flat universe, the density-to-pressure ratio of dark energy can be similarly measured to 5% for the present value w0 and ~0.1 for the time variation w'. The large survey area, depth, spatial resolution, time-sampling, and nine-band optical to NIR photometry will support additional independent and/or complementary dark-energy measurement approaches as well as a broad range of auxiliary science programs. (Abridged)

Absorption between the rest-frame wavelengths of 973 and 1026 Angstroms in quasar spectra arises from two sources (apart from occasional metals): one is due to Lyman-alpha (Lya) absorption by materials at a low redshift, and the other from Lyman-beta (Lyb) at a higher redshift. These two sources of absorption are to a good approximation uncorrelated because of their wide physical separation. Therefore, the two-point correlation of absorption in this region of quasar spectra neatly factorizes into two pieces: the Lyb correlation at high z, and the Lya correlation at low z. The latter can be independently measured from quasar spectra at lower redshifts using current techniques. A simple division then offers a way to statistically separate out the Lyb two-point correlation from the Lya correlation. Several applications of this technique are discussed. First, since the Lyb absorption cross-section is lower than Lya by about a factor of 5, the Lyb forest is a better probe of the intergalactic medium (IGM) at higher redshifts where Lya absorption is often saturated. Second, for the same reason, the Lyb forest allows a better measurement of the equation of state of the IGM at higher overdensities, yielding stronger constraints on its slope when used in conjunction with the Lya forest. Third, models of the Lya forest based on gravitational instability make unique predictions for the Lyb forest, which can be tested against observations. We briefly point out that feedback processes that affect higher density regions but leave low density structure intact may be better constrained by the Lyb forest.

We measure cosmological parameters using the three-dimensional power spectrum P(k) from over 200,000 galaxies in the Sloan Digital Sky Survey (SDSS) in combination with WMAP and other data. Our results are consistent with a ``vanilla'' flat adiabatic Lambda-CDM model without tilt (n=1), running tilt, tensor modes or massive neutrinos. Adding SDSS information more than halves the WMAP-only error bars on some parameters, tightening 1 sigma constraints on the Hubble parameter from h~0.74+0.18-0.07 to h~0.70+0.04-0.03, on the matter density from Omega_m~0.25+/-0.10 to Omega_m~0.30+/-0.04 (1 sigma) and on neutrino masses from <11 eV to <0.6 eV (95%). SDSS helps even more when dropping prior assumptions about curvature, neutrinos, tensor modes and the equation of state. Our results are in substantial agreement with the joint analysis of WMAP and the 2dF Galaxy Redshift Survey, which is an impressive consistency check with independent redshift survey data and analysis techniques. In this paper, we place particular emphasis on clarifying the physical origin of the constraints, i.e., what we do and do not know when using different data sets and prior assumptions. For instance, dropping the assumption that space is perfectly flat, the WMAP-only constraint on the measured age of the Universe tightens from t0~16.3+2.3-1.8 Gyr to t0~14.1+1.0-0.9 Gyr by adding SDSS and SN Ia data. Including tensors, running tilt, neutrino mass and equation of state in the list of free parameters, many constraints are still quite weak, but future cosmological measurements from SDSS and other sources should allow these to be substantially tightened.

Given a foreground galaxy-density field or shear field, its cross-correlation with the shear field from a background population of source galaxies scales with the source redshift in a way that is specific to lensing. Such a source-scaling can be exploited to effectively measure geometrical distances as a function of redshift and thereby constrain dark energy properties, free of any assumptions about the galaxy-mass/mass power spectrum (its shape, amplitude or growth). Such a geometrical method can yield a ~ 0.03 - 0.07 f_{sky}^{-1/2} measurement on the dark energy abundance and equation of state, for a photometric redshift accuracy of dz ~ 0.01 - 0.05 and a survey with median redshift of ~ 1. While these constraints are weaker than conventional weak lensing methods, they provide an important consistency check because the geometrical method carries less theoretical baggage: there is no need to assume any structure formation model (e.g. CDM). The geometrical method is at the most conservative end of a whole spectrum of methods which obtain smaller errorbars by making more restrictive assumptions -- we discuss some examples. Our geometrical approach differs from previous investigations along similar lines in three respects. First, the source-scaling we propose to use is less demanding on the photometric redshift accuracy. Second, the scaling works for both galaxy-shear and shear-shear correlations. Third, we find that previous studies underestimate the statistical errors associated with similar geometrical methods, the origin of which is discussed.

The microlensing interpretation of the optical afterglow of GRB 000301C seems naively surprising, since a simple estimate of the stellar microlensing rate gives less than one in four hundred for a flat Omega_Lambda=0.7 cosmology, whereas one event was seen in about thirty afterglows. Considering baryonic MACHOs making up half of the baryons in the universe, the microlensing probability per burst can be roughly 5% for a GRB at redshift z=2. We explore two effects that may enhance the probability of observing microlensed gamma-ray burst afterglows: binary lenses and double magnification bias. We find that the consideration of binary lenses can increase the rate only at the ~15% level. On the other hand, because gamma-ray bursts for which afterglow observations exist are typically selected based on fluxes at widely separated wavebands which are not necessarily well correlated (e.g. localization in X-ray, afterglow in optical/infrared), magnification bias can operate at an enhanced level compared to the usual single-bias case. We find that existing estimates of the slope of the luminosity function of gamma-ray bursts, while as yet quite uncertain, point to enhancement factors of more than three above the simple estimates of the microlensing rate. We find that the probability to observe at least one microlensing event in the sample of 27 measured afterglows can be 3-4% for stellar lenses, or as much as 25 Omega_lens for baryonic MACHOs. We note that the probability to observe at least one event over the available sample of afterglows is significant only if a large fraction of the baryons in the universe are condensed in stellar-mass objects. (ABRIDGED)

We measure the large-scale real-space power spectrum P(k) using a sample of 205,443 galaxies from the Sloan Digital Sky Survey, covering 2417 square degrees with mean redshift z~0.1. We employ a matrix-based method using pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 22 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.02 h/Mpc < k < 0.3h/Mpc. We pay particular attention to modeling, quantifying and correcting for potential systematic errors, nonlinear redshift distortions and the artificial red-tilt caused by luminosity-dependent bias. Our final result is a measurement of the real-space matter power spectrum P(k) up to an unknown overall multiplicative bias factor. Our calculations suggest that this bias factor is independent of scale to better than a few percent for k<0.1h/Mpc, thereby making our results useful for precision measurements of cosmological parameters in conjunction with data from other experiments such as the WMAP satellite. As a simple characterization of the data, our measurements are well fit by a flat scale-invariant adiabatic cosmological model with h Omega_m =0.201+/- 0.017 and L* galaxy sigma_8=0.89 +/- 0.02 when fixing the baryon fraction Omega_b/Omega_m=0.17 and the Hubble parameter h=0.72; cosmological interpretation is given in a companion paper.

A comparison between the line of sight power spectrum of absorption in the Lyman-alpha forest and the cross power spectrum between the absorption in neighboring lines of sight offers an evolution-free means to constrain the cosmological constant, or dark energy. Using cosmological simulations, we consider a maximum likelihood method to obtain constraints from this comparison. In our method, measurements of the auto and cross spectra from observations are compared with those from a multi-parameter grid of simulated models of the intergalactic medium. We then marginalize over nuisance parameters to obtain constraints on the cosmological constant. Redshift space distortions due to peculiar velocities and thermal broadening, a potential difficulty in applying this test, are explicitly modeled in our simulations. To illustrate our method, we measure the cross spectrum from a new sample of five close quasar pairs, with separations of 0.5 to 3 arcmin. We attempt to obtain a constraint on Omega_Lambda, but find only weak constraints. An Einstein-de-Sitter cosmology is, however, disfavored by the data at a ~ 2 sigma confidence level. We consider the power of future observations, paying particular attention to the effects of spectral resolution and shot-noise. We find that ~ 50 moderate resolution, FWHM ~ 150 km/s, close separation (~ 30-120 arcsec) pairs should allow a (2 sigma) constraint on Omega_Lambda at the level of 15%, if other modeling parameters are determined through other means. We find that there is a sizeable gain from observing very close, ~ 30 arcsec, separation pairs provided they are observed with high spectral resolution. A sample of ~ 10 such pairs gives similar constraints to the 50 moderate resolution pairs mentioned above.

We demonstrate that the gravity wave background amplitude implies a robust upper bound on the ratio: \lambda / H^{-1} < e^60, where \lambda is the proper wavelength of fluctuations of interest and H^{-1} is the horizon at the end of inflation. The bound holds as long as the energy density of the universe does not drop faster than radiation subsequent to inflation. This limit implies that the amount of expansion between the time the scales of interest leave the horizon and the end of inflation, denoted by e^N, is also bounded from above, by about e^60 times a factor that involves an integral over the first slow-roll parameter. In other words, the bound on N is model dependent -- we show that for vast classes of slow-roll models, N < 67. The quantities, \lambda / H^{-1} or N, play an important role in determining the nature of inflationary scalar and tensor fluctuations. We suggest ways to incorporate the above bounds when confronting inflation models with observations. As an example, this bound solidifies the tension between observations of cosmic microwave background (CMB) anisotropies and chaotic inflation with a \phi^4 potential by closing the escape hatch of large N (< 62).

The recent measurement by WMAP of a large electron scattering optical depth tau_e = 0.17 +- 0.04 is consistent with a simple model of reionization in which the intergalactic medium (IGM) is ionized at redshift z ~ 15, and remains highly ionized thereafter. Here, we show that existing measurements of the IGM temperature from the Lyman-alpha forest at z ~ 2 - 4 rule out this ``vanilla'' model. Under reasonable assumptions about the ionizing spectrum, as long as the universe is reionized before z = 10, and remains highly ionized thereafter, the IGM reaches an asymptotic thermal state which is too cold compared to observations. To simultaneously satisfy the CMB and forest constraints, the reionization history must be complex: reionization begins early at z >~ 15, but there must have been significant (order unity) changes in fractions of neutral hydrogen and/or helium at 6 < z < 10, and/or singly ionized helium at 4 < z < 10. We describe a physically motivated reionization model that satisfies all current observations. We also explore the impact of a stochastic reionization history and show that a late epoch of (HeII --> HeIII) reionization induces a significant scatter in the IGM temperature, but the scatter diminishes with time quickly. Finally, we provide an analytic formula for the thermal asymptote, and discuss possible additional heating mechanisms that might evade our constraints.

The flux/transmission power spectrum has become a popular statistical tool in studies of the high redshift ($z > 2$) Lyman-alpha forest. At low redshifts, where the forest has thinned out into a series of well-isolated absorption lines, the motivation for flux statistics is less obvious. Here, we study the relative merits of flux versus line correlations, and derive a simple condition under which one is favored over the other on purely statistical grounds. Systematic errors probably play an important role in this discussion, and they are outlined as well.

We discuss the determination of the linear power spectrum of dark matter fluctuations from measurements of the transmitted flux Ly-alpha Forest power spectrum. We show that at most scales probed by current measurements, the flux power spectrum is significantly affected by non-linear corrections to the linear dark matter power spectrum due to gravitational clustering. Inferring the linear dark matter power spectrum shape is therefore difficult due to non-linear corrections driving the non-linear power spectrum to a k^{-1.4} shape nearly independent of initial conditions. We argue that some methods used in previous estimates underestimate the uncertainties in the shape of the linear dark matter power spectrum.

XENON is a novel liquid xenon experiment concept for a sensitive dark matter search using a 1-tonne active target, distributed in an array of ten independent time projection chambers. The design relies on the simultaneous detection of ionization and scintillation signals in liquid xenon, with the goal of extracting as much information as possible on an event-by-event basis, while maintaining most of the target active. XENON is expected to have effective and redundant background identification and discrimination power, higher than 99.5%, and to achieve a very low threshold, on the order of 4 keV visible recoil energy. Based on this expectation and the 1-tonne mass of active xenon, we project a sensitivity of 0.0001 events/kg/day, after 3 yr operation in an appropriate underground location. The XENON experiment has been recently proposed to the National Science Foundation (NSF) for an initial development phase leading to the development of the 100 kg unit module.

Recent observations of high redshift quasar spectra reveal long gaps with little flux. A small or no detectable flux does not by itself imply the intergalactic medium (IGM) is neutral. Inferring the average neutral fraction from the observed absorption requires assumptions about clustering of the IGM, which the gravitational instability model supplies. Our most stringent constraint on the neutral fraction at z ~ 6 is derived from the mean Lyman-beta transmission measured from the z=6.28 SDSS quasar of Becker et al. -- the neutral hydrogen fraction at mean density has to be larger than 4.7 times 10^{-4}. This is substantially higher than the neutral fraction of ~ 3-5 times 10^{-5} at z = 4.5 - 5.7, suggesting that dramatic changes take place around or just before z ~ 6, even though current constraints are still consistent with a fairly ionized IGM at z ~ 6. An interesting alternative method to constrain the neutral fraction is to consider the probability of having many consecutive pixels with little flux, which is small unless the neutral fraction is high. This constraint is slightly weaker than the one obtained from the mean transmission. We show that while the derived neutral fraction at a given redshift is sensitive to the power spectrum normalization, the size of the jump around z ~ 6 is not. We caution that systematic uncertainties include spatial fluctuations in the ionizing background, and the continuum placement. Tests are proposed. In particular, the sightline to sightline dispersion in mean transmission might provide a useful diagnostic. We express the dispersion in terms of the transmission power spectrum, and develop a method to calculate the dispersion for spectra that are longer than the typical simulation box.

We review the current status of cosmological parameters, dark energy and large-scale structure, from a theoretical and observational perspective. We first present the basic cosmological parameters and discuss how they are measured with different observational techniques. We then describe the recent evidence for dark energy from Type Ia supernovae. Dynamical models of the dark energy, quintessence, are then described, as well as how they relate to theories of gravity and particle physics. The basic theory of structure formation via gravitational instability is then reviewed. Finally, we describe new observational probes of the large-structure of the universe, and how they constrain cosmological parameters.

The angular distribution of galaxies encodes a wealth of information about large scale structure. Ultimately, the Sloan Digital Sky Survey (SDSS) will record the angular positions of order 10^8 galaxies in five bands, adding significantly to the cosmological constraints. This is the first in a series of papers analyzing a rectangular stripe 2.5x90 degrees from early SDSS data. We present the angular correlation function for galaxies in four separate magnitude bins on angular scales ranging from 0.003 degrees to 15 degrees. Much of the focus of this paper is on potential systematic effects. We show that the final galaxy catalog -- with the mask accounting for regions of poor seeing, reddening, bright stars, etc. -- is free from external and internal systematic effects for galaxies brighter than r* = 22. Our estimator of the angular correlation function includes the effects of the integral constraint and the mask. The full covariance matrix of errors in these estimates is derived using mock catalogs with further estimates using a number of other methods.

We compute both extrinsic (lensing) and intrinsic contributions to the (galaxy-)density-ellipticity correlation function, the latter done using current analytic theories of tidal alignment. The gravitational lensing contribution has two components: one is analogous to galaxy-galaxy lensing and the other arises from magnification bias -- that gravitational lensing induces a modulation of the galaxy density as well as ellipticity. On the other hand, the intrinsic alignment contribution vanishes, even after taking into account source clustering corrections, which suggests the density-ellipticity correlation might be an interesting diagnostic in differentiating between intrinsic and extrinsic alignments. {\it However}, an important assumption, commonly adopted by current analytic alignment theories, is the Gaussianity of the tidal field. Inevitable non-Gaussian fluctuations from gravitational instability induces a non-zero intrinsic density-ellipticity correlation, which we estimate. We also argue that non-Gaussian contributions to the intrinsic {\it ellipticity-ellipticity} correlation are often non-negligible. This leads to a linear rather than, as is commonly assumed, quadratic scaling with the power spectrum on sufficiently large scales. Finally, we estimate the contribution of intrinsic alignment to low redshift galaxy-galaxy lensing measurements (e.g. SDSS), due to the partial overlap between foreground and background galaxies: the intrinsic contamination is about 10 - 30 % at 10'. Uncertainties in this estimate are discussed.

Conventional Cold Dark Matter cosmological models predict small scale structures, such as cuspy halos, which are in apparent conflict with observations. Several alternative scenarios based on modifying fundamental properties of the dark matter have been proposed. We show that general principles of quantum mechanics, in particular unitarity, imply interesting constraints on two proposals: collisional dark matter proposed by Spergel & Steinhardt, and strongly annihilating dark matter proposed by Kaplinghat, Knox & Turner. Efficient scattering required in both implies m < 12 GeV and m < 25 GeV respectively. The same arguments show that the strong annihilation in the second scenario implies the presence of significant elastic scattering, particularly for large enough masses. Recently, a variant of the collisional scenario has been advocated to satisfy simultaneously constraints from dwarf galaxies to clusters, with a cross section that scales inversely with velocity. We show that this scenario likely involves super-elastic processes, and the associated kinetic energy change must be taken into account when making predictions. Exceptions and implications for experimental searches are discussed.

Recent discussions suggest the possibility that short distance physics can significantly modify the behavior of quantum fluctuations in the inflationary universe, and alter the standard large scale structure predictions. Such modifications can be viewed as due to a different choice of the vacuum state. We show that such changes generally lead to violations of the well-known consistency relation between the scalar to tensor ratio and the tensor spectral index. Vacuum effects can introduce an observable modulation to the usual predictions for the scalar and tensor power spectra.

(Abridged) The solar system gas giant planets are oblate due to their rapid rotation. A measurement of the planet's projected oblateness would constrain the planet's rotational period. Planets that are synchronously rotating with their orbital revolution will be rotating too slowly to be significantly oblate; these include planets with orbital semi-major axes <~ 0.2AU (for M_P ~ M_J and M_* ~ M_sun). Jupiter-like planets in the range of orbital semi-major axes 0.1 AU to 0.2 AU will tidally evolve to synchronous rotation on a timescale similar to main sequence stars' lifetimes. In this case an oblateness detection will help constrain the planet's tidal Q value. The projected oblateness of a transiting extrasolar giant planet is measurable from a very high-photometric-precision transit light curve. For a sun-sized star and a Jupiter-sized planet the normalized flux difference in the transit ingress/egress light curve between a spherical and an oblate planet is a few to 15 x 10^{-5} for oblateness similar to Jupiter and Saturn respectively. The transit ingress and egress are asymmetric for an oblate planet with an orbital inclination different from 90 degrees and a non-zero projected obliquity. A photometric precision of 10^{-4} has been reached by HST observations of the known transiting extrasolar planet HD209458b. In addition several planned space missions (including MOST, MONS, Corot, and Kepler) will study or find transiting giant planets with photometric precision <~ 10^{-5}. In particular Kepler is expected to find ~40 transiting giant planets < 1 AU with ~20 of those > 0.2 AU.

We compute the angular power spectrum C_l from 1.5 million galaxies in early SDSS data on large angular scales, l<600. The data set covers about 160 square degrees, with a characteristic depth of order 1 Gpc/h in the faintest (21<r<22) of our four magnitude bins. Cosmological interpretations of these results are presented in a companion paper by Dodelson et al (2001). The data in all four magnitude bins are consistent with a simple flat ``concordance'' model with nonlinear evolution and linear bias factors of order unity. Nonlinear evolution is particularly evident for the brightest galaxies. A series of tests suggest that systematic errors related to seeing, reddening, etc., are negligible, which bodes well for the sixtyfold larger sample that the SDSS is currently collecting. Uncorrelated error bars and well-behaved window functions make our measurements a convenient starting point for cosmological model fitting.

We present measurements of parameters of the 3-dimensional power spectrum of galaxy clustering from 222 square degrees of early imaging data in the Sloan Digital Sky Survey. The projected galaxy distribution on the sky is expanded over a set of Karhunen-Loeve eigenfunctions, which optimize the signal-to-noise ratio in our analysis. A maximum likelihood analysis is used to estimate parameters that set the shape and amplitude of the 3-dimensional power spectrum. Our best estimates are Gamma=0.188 +/- 0.04 and sigma_8L = 0.915 +/- 0.06 (statistical errors only), for a flat Universe with a cosmological constant. We demonstrate that our measurements contain signal from scales at or beyond the peak of the 3D power spectrum. We discuss how the results scale with systematic uncertainties, like the radial selection function. We find that the central values satisfy the analytically estimated scaling relation. We have also explored the effects of evolutionary corrections, various truncations of the KL basis, seeing, sample size and limiting magnitude. We find that the impact of most of these uncertainties stay within the 2-sigma uncertainties of our fiducial result.

The Sloan Digital Sky Survey is one of the first multicolor photometric and spectroscopic surveys designed to measure the statistical properties of galaxies within the local Universe. In this Letter we present some of the initial results on the angular 2-point correlation function measured from the early SDSS galaxy data. The form of the correlation function, over the magnitude interval 18<r*<22, is shown to be consistent with results from existing wide-field, photographic-based surveys and narrower CCD galaxy surveys. On scales between 1 arcminute and 1 degree the correlation function is well described by a power-law with an exponent of ~ -0.7. The amplitude of the correlation function, within this angular interval, decreases with fainter magnitudes in good agreement with analyses from existing galaxy surveys. There is a characteristic break in the correlation function on scales of approximately 1-2 degrees. On small scales, < 1', the SDSS correlation function does not appear to be consistent with the power-law form fitted to the 1'< theta <0.5 deg data. With a data set that is less than 2% of the full SDSS survey area, we have obtained high precision measurements of the power-law angular correlation function on angular scales 1' < theta < 1 deg, which are robust to systematic uncertainties. Because of the limited area and the highly correlated nature of the error covariance matrix, these initial results do not yet provide a definitive characterization of departures from the power-law form at smaller and larger angles. In the near future, however, the area of the SDSS imaging survey will be sufficient to allow detailed analysis of the small and large scale regimes, measurements of higher-order correlations, and studies of angular clustering as a function of redshift and galaxy type.

Early photometric data from the Sloan Digital Sky Survey (SDSS) contain angular positions for 1.5 million galaxies. In companion papers, the angular correlation function $w(\theta)$ and 2D power spectrum $C_l$ of these galaxies are presented. Here we invert Limber's equation to extract the 3D power spectrum from the angular results. We accomplish this using an estimate of $dn/dz$, the redshift distribution of galaxies in four different magnitude slices in the SDSS photometric catalog. The resulting 3D power spectrum estimates from $w(\theta)$ and $C_l$ agree with each other and with previous estimates over a range in wavenumbers $0.03 < k/{\rm h Mpc}^{-1} < 1$. The galaxies in the faintest magnitude bin ($21 < \rstar < 22$, which have median redshift $z_m=0.43$) are less clustered than the galaxies in the brightest magnitude bin ($18 < \rstar < 19$ with $z_m=0.17$), especially on scales where nonlinearities are important. The derived power spectrum agrees with that of Szalay et al. (2001) who go directly from the raw data to a parametric estimate of the power spectrum. The strongest constraints on the shape parameter $\Gamma$ come from the faintest galaxies (in the magnitude bin $21 < \rstar < 22$), from which we infer $\Gamma = 0.14^{+0.11}_{-0.06}$ (95% C.L.).

We reassess constraints on the cosmological baryon density from observations of the mean decrement and power spectrum of the Lyman-alpha forest, taking into account uncertainties in all free parameters in the simplest gravitational instability model. The uncertainty is dominated by that of the photoionizing background, but incomplete knowledge of the thermal state of the intergalactic medium also contributes significantly to the error-budget. While current estimates of the baryon fraction from the forest do prefer values that are somewhat higher than the big bang nucleosynthesis value of \Omega_b h^2 = 0.02 \pm 0.001, the discrepancy is at best about 3 \sigma. For instance, assuming the highest estimate of the ionizing background, as indicated by recent measurements of a large escape fraction from Lyman-break galaxies by Steidel, Pettini & Adelberger, we find \Omega_b h^2 = 0.045 \pm 0.008. A recent measurement of the ionizing background from the proximity effect by Scott et al., on the other hand, implies \Omega_b h^2 = 0.03 \pm 0.01. We provide an expression from which future likelihoods for \Omega_b h^2 can be derived as measurements of the ionizing background improve -- consistency among constraints from the forest, nucleosynthesis and the microwave background will provide a powerful test of the gravitational instability model. We also develop a formalism which treats lower bounds on the baryon density in a statistical manner, which is appropriate if only a lower bound on the ionizing background is known. Finally, we discuss the implications of the escape fraction measurement for the age, structure and stellar content of Lyman-break galaxies.

Future high precision photometric measurements of transiting extrasolar planets promise to tell us much about the nature of these systems. We examine how atmospheric lensing and (projected) planet oblateness/ellipticity modify transit light curves. The large density gradients expected in planet atmospheres can offset the unfavorably large observer-lens to source-lens distance ratio, and allow the existence of caustics. Under such conditions, starlight from all points in the planet's shadow is refracted into view producing a characteristic slowing down of the dimming at ingress (vice versa for egress). A search over several parameters, the planet radius etc., cannot produce a nonlensed transit light curve that can mimic a lensed light curve. The fractional change in the diminution of starlight is approximately the ratio of atmospheric scale height to planet radius, expected to be 1 % or less. Planet oblateness induces an asymmetry to the transit light curve about the point of minimum flux, depending on the planet orientation. The fractional asymmetry is around 0.5 % for a projected ellipticity of 10 %, independent of lensing. For favorable ratios of planet to star radius, the above effects are potentially observable with future space-based missions. This will allow us to constrain the planet shape, and its atmospheric scale height, density and refractive coefficient, providing information on its rotation, temperature and composition. For HD 209458b, the only currently known transiting extrasolar planet, caustics are absent because of the very small lens-source separation. Finally, we provide estimates of other variations to transit light curves that could be of comparable importance, including rings, satellites, stellar oscillations, star spots, and weather.

The lifetime t_Q of the luminous phase of quasars is constrained by current observations to be between 10^6 and 10^8 years, but is otherwise unkown. We model the quasar luminosity function in detail in the optical and X-ray bands using the Press-Schechter formalism, and show that the expected clustering of quasars depends strongly on their assumed lifetime. We quantify this dependence, and find that existing measurements of the correlation length of quasars are consistent with the range 10^6 < t_Q < 10^8 years. We then show that future measurements of the power spectrum of quasars out to z=3, from the 2dF or Sloan Digital Sky Survey, can significantly improve this constraint, and in principle allow a precise determination of t_Q. We estimate the systematic errors introduced by uncertainties in the modeling of the quasar-halo relationship, as well as by the possible existence of obscured quasars.

We investigate correlations between the long wavelength fluctuations and the small scale power in the Lyman alpha forest. We show that such correlations can be used to discriminate between fluctuations induced by large scale structure and those produced by non-gravitational processes such as fluctuations in the continuum of the quasar or in the temperature of the gas. The correlations observed in Q1422+231 are in agreement with the predictions of numerical simulations, indicating that non-gravitational fluctuations on large scales have to be small compared those induced by the large scale structure of the universe, contributing less than 10 % (95 % confidence) of the observed power. We also show the sensitivity of such statistics to the equation of state of the gas and its mean temperature.

We use published measurements of the transmission power spectrum of the Lyman alpha forest to constrain several parameters that describe cosmology and thermal properties of the intergalactic medium (IGM). A 6 parameter grid is constructed using Particle-Mesh dark matter simulations together with scaling relations to make predictions for the gas properties. We fit for all parameters simultaneously and identify several degeneracies. We find that the temperature of the IGM can be well determined from the fall-off of the power spectrum at small scales. We find a temperature around 20.000 K, dependent on the slope of the gas equation of state. We see no evidence for evolution in the IGM temperature. We place constraints on the amplitude of the dark matter fluctuations. However, contrary to previous results, the slope of the dark matter power spectrum is poorly constrained. This is due to uncertainty in the effective Jeans smoothing scale, which depends on the temperature as well as the thermal history of the gas.

A simple model of how objects of different masses stream towards each other as they cluster gravitationally is described. The model shows how the mean streaming velocity of dark matter particles is related to the motions of the parent dark matter haloes. It also provides a reasonably accurate description of how the pairwise velocity dispersion of dark matter particles differs from that of the parent haloes. The analysis is then extended to describe the streaming motions of galaxies. This shows explicitly that the streaming motions measured in a given galaxy sample depend on how the sample was selected, and shows how to account for this dependence on sample selection. In addition,we show that the pairwise dispersion should also depend on sample type. Our model predicts that, on small scales, redshift space distortions should affect red galaxies more strongly than blue.

We write the correlation function of dark matter particles, xi(r), as the sum of two terms - one which accounts for nonlinear evolution, and dominates on small scales, and another which is essentially the term from linear theory, and dominates on large scales. We use models of the number and spatial distribution of haloes and halo density profiles to describe the nonlinear term and its evolution. The result provides a good description of the evolution of xi(r) in simulations. We then use this decomposition to provide simple and accurate models of how the single particle velocity dispersion evolves with time, and how the first and second moments of the pairwise velocity distribution depend on scale. The key idea is to use the simple physics of linear theory on large scales, the simple physics of the virial theorem on small scales, and our model for the correlation function to tell us how to weight the two types of contributions (linear and nonlinear) to the pairwise velocity statistics. When incorporated into the streaming model, our results will allow a simple accurate description of redshift-space distortions over the entire range of linear to highly nonlinear regimes.

Superheavy ($M>10^{10}$ GeV) particles produced during inflation may be the dark matter, independent of their interaction strength. Strongly interacting superheavy particles will be captured by the sun, and their annihilation in the center of the sun will produce a flux of energetic neutrinos that should be detectable by neutrino telescopes. Depending on the particle mass, event rates in a cubic-kilometer detector range from several per hour to several per year. The signature of the process is a predominance of tau neutrinos, with a relatively flat energy spectrum of events ranging from 50 GeV to many TeV, and with the mean energy of detected tau neutrinos about 3 TeV.

We study galaxy clustering using halo models, where gravitational clustering is described in terms of dark matter halos. At small scales, clustering statistics are dominated by halo density profiles, whereas at large scales, correlations are the result of combining non-linear perturbation theory with halo biasing. Galaxies are assumed to follow the dark matter profiles, and galaxy formation efficiency is given by the number of galaxies as a function of halo mass. This approach leads to generic predictions: the galaxy power spectrum shows a power-law behavior even though the dark matter does not, and the galaxy higher-order correlations show smaller amplitudes at small scales than their dark matter counterparts, as observed in galaxy catalogs. We find that requiring to fit both the second and third order moments of the APM galaxies provides a strong constraint on galaxy formation models. The data at large scales require that galaxy formation be relatively efficient at small masses, m =10^10 Msun/h, whereas data at smaller scales require that the number of galaxies in a halo scale as the mass to the 0.8th power in the high-mass limit. These constraints are independent of those derived from the luminosity function or Tully-Fisher relation. We also predict the power spectrum, bispectrum, and higher-order moments of the mass density field. Although halo models agree well with measurements of the mass power spectrum and the higher order Sp parameters in N-body simulations, the model assumption that halos are spherical leads to disagreement in the configuration dependence of the bispectrum at small scales. We stress the importance of finite volume effects in higher-order statistics and show how they can be estimated in this approach.

We present new methods to minimize the systematic and random errors for measuring the transmission power spectrum from the Lyman-alpha forest. Sources of systematic errors explored include metal line contamination and continuum-fitting. We advocate the technique of trend-removal in place of traditional continuum-fitting -- here, a spectrum is normalized by its (smoothly varying) mean rather than its continuum -- this method is easily automated and removes biases introduced by continuum-fitting. Trend- removal can be easily applied to spectra where continuum-fitting is difficult, such as when the resolution or signal-to-noise is low, or for spectra at high redshifts. Furthermore, a measurement of the continuum power spectrum using trend-removal, from either low redshift quasar spectra or the red-side of Lyman-alpha, allows in principle the removal of spurious power introduced by the continuum and thereby expanding scales probed to larger ones. We also derive expressions for the shot-noise bias and variance of the power spectrum estimate, taking into account the non-Poissonian nature of the shot-noise and the non-Gaussianity of the cosmic fluctuations. An appropriate minimum variance weighting of the data is given. Finally, we give practical suggestions on observing strategy: the desired resolution and S/N for different purposes, and how to distribute one's finite observing time among quasar targets. Also discussed is the quasar spectroscopic study of the Sloan Digital Sky Survey, which has the potential to measure the power spectrum at z ~ 2-4 accurate to better than 1 % per mode -- the techniques presented here will be useful for tackling the anticipated issues of shot-noise and continuum contamination.

It is quite plausible that the mass of the dark matter particle increases significantly after its freeze-out, due to a scalar field rolling to large values. We describe a realization of this scenario in the context of thermal inflation which naturally gives a cold dark matter particle with the correct cosmological abundance and a mass around $10^{10} {\rm GeV}$, evading the conventional upper bound of $10^5 {\rm GeV}$. We also discuss another realization which could produce a cosmologically interesting abundance of near Planck mass, possibly electromagnetically charged, particles. The detection and observational consequences of superheavy cold dark matter or WIMPZILLAs are briefly examined.

Gravitational clustering is an intrinsically non-linear process that generates significant non-Gaussian signatures in the density field. We consider how these affect power spectrum determinations from galaxy and weak-lensing surveys. Non-Gaussian effects not only increase the individual error bars compared to the Gaussian case but, most importantly, lead to non-trivial cross-correlations between different band-powers. We calculate the power-spectrum covariance matrix in non-linear perturbation theory (weakly non-linear regime), in the hierarchical model (strongly non-linear regime), and from numerical simulations in real and redshift space. We discuss the impact of these results on parameter estimation from power spectrum measurements and their dependence on the size of the survey and the choice of band-powers. We show that the non-Gaussian terms in the covariance matrix become dominant for scales smaller than the non-linear scale, depending somewhat on power normalization. Furthermore, we find that cross-correlations mostly deteriorate the determination of the amplitude of a rescaled power spectrum, whereas its shape is less affected. In weak lensing surveys the projection tends to reduce the importance of non-Gaussian effects. Even so, for background galaxies at redshift z=1, the non-Gaussian contribution rises significantly around l=1000, and could become comparable to the Gaussian terms depending upon the power spectrum normalization and cosmology. The projection has another interesting effect: the ratio between non-Gaussian and Gaussian contributions saturates and can even decrease at small enough angular scales if the power spectrum of the 3D field falls faster than 1/k^2.

Bernardeau et al. (1997), using perturbation theory, showed that the skewness of the large-scale lensing-convergence, or projected mass density, could be used to constrain $\Omega_m$, the matter content of the universe. On the other hand, deep weak-lensing field surveys in the near future will likely measure the convergence on small angular scales (< 10 arcmin.), where the signal will be dominated by highly nonlinear fluctuations. We develop a method to compute the small-scale convergence skewness, using a prescription for the highly nonlinear three-point function developed by Scoccimarro and Frieman (1998). This method gives predictions that agree well with existing results from ray-tracing N-body simulations, but is significantly faster, allowing the exploration of a large number of models. We demonstrate that the small-scale convergence skewness is insensitive to the shape and normalization of the primordial (CDM-type) power spectrum, making it dependent almost entirely on the cosmological energy contents, through their influence on the global geometrical distances and fluctuation growth rate. Moreover, nonlinear clustering appears to enhance the differences between predictions of the convergence skewness for a range of models. Hence, in addition to constraining $\Omega_m$, the small-scale convergence skewness from future deep several- degree-wide surveys can be used to differentiate between curvature dominated and cosmological constant ($\Lambda$) dominated models, as well as to constrain the equation of state of a quintessence component, thereby distinguishing $\Lambda$ from quintessence as well. Finally, our method can be easily generalized to other measures such as aperture mass statistics.

We examine several aspects of redshift distortions by expressing the redshift-space density in terms of the eigenvalues and orientation of the local Lagrangian deformation tensor. We explore the importance of multiple-streaming using the Zel'dovich approximation (ZA), and compute the average number of streams in real and redshift-space. It is found that multiple-streaming can be significant in redshift-space but negligible in real-space, even at moderate values of the linear fluctuation amplitude ($\sigma < 1$). Moreover, unlike their real-space counter-parts, redshift-space multiple-streams can flow past each other with minimal interactions. Such nonlinear redshift-space effects, which operate even when the real-space density field is quite linear, could suppress the classic compression of redshift-structures predicted by linear theory (Kaiser 1987). We also compute using the ZA the probability distribution function (PDF) of density, as well as $S_3$, in real and redshift-space, and compare it with the PDF measured from N-body simulations. The role of caustics in defining the character of the high density tail is examined. It is found that (non-Lagrangian) smoothing, due to both finite resolution or discreteness and small-scale velocity dispersions, is very effective in erasing caustic structures, unless the initial power spectrum is sufficiently truncated.

In recent years, a remarkably simple physical picture of the Lyman-alpha forest has emerged, which allows detailed predictions to be made and turns the forest into a powerful probe of cosmology. We point out ways in which such a picture can be tested observationally, and explore three areas in which the Lyman-alpha forest can yield valuable constraints: the reionization history, the primordial mass power spectrum and the cosmological constant or its variants. The possibility of combining with other high redshift observations, such as the Lyman-break galaxy surveys, to provide consistency checks and complementary information is also discussed.

Nonlinear combinations of direct observables are often used to estimate quantities of theoretical interest. Without sufficient caution, this could lead to biased estimations. An example of great interest is the skewness $S_3$ of the galaxy distribution, defined as the ratio of the third moment $\xibar_3$ and the variance squared $\xibar_2^2$. Suppose one is given unbiased estimators for $\xibar_3$ and $\xibar_2^2$ respectively, taking a ratio of the two does not necessarily result in an unbiased estimator of $S_3$. Exactly such an estimation-bias affects most existing measurements of $S_3$. Furthermore, common estimators for $\xibar_3$ and $\xibar_2$ suffer also from this kind of estimation-bias themselves: for $\xibar_2$, it is equivalent to what is commonly known as the integral constraint. We present a unifying treatment allowing all these estimation-biases to be calculated analytically. They are in general negative, and decrease in significance as the survey volume increases, for a given smoothing scale. We present a re-analysis of some existing measurements of the variance and skewness and show that most of the well-known systematic discrepancies between surveys with similar selection criteria, but different sizes, can be attributed to the volume-dependent estimation-biases. This affects the inference of the galaxy-bias(es) from these surveys. Our methodology can be adapted to measurements of analogous quantities in quasar spectra and weak-lensing maps. We suggest methods to reduce the above estimation-biases, and point out other examples in LSS studies which might suffer from the same type of a nonlinear-estimation-bias.

In this Letter we explore a version of the test of cosmological geometry proposed by Alcock and Paczynski (1979), using observations of the Lyman-alpha forest in the spectra of close quasar pairs. By comparing the correlations in absorption in one quasar spectrum with correlations between the spectra of neighboring quasars one can determine the relation of the redshift distance scale to the angle distance scale at the redshift of the absorbers, $z \sim 2 - 4$. Since this relationship depends on the parameters of the cosmological model, these parameters may be determined using the Lyman-alpha forest. While this test is relatively insensitive to the density parameter $\Omega_m$ in a dust-dominated universe, it is more sensitive to the presence of a matter component with large negative pressure (such as a cosmological constant $\Lambda$) and its equation of state. With only 25 pairs of quasar spectra at angular separations $0.5' - 2'$, one can discriminate between an $\Omega_m = 0.3$ open universe ($\Lambda=0$) and an $\Omega_m = 0.3$ flat ($\Lambda$-dominated) universe at the $4-\sigma$ level. The S/N can be enhanced by considering quasar pairs at smaller angular separations, but requires proper modeling of nonlinear redshift space distortions. Here the correlations and redshift space distortions are modeled using linear theory.

We propose a method for recovering the shape of the mass power spectrum on large scales from the transmission fluctuations of the Lyman-alpha forest, which takes into account directly redshift-space distortions. The procedure, in discretized form, involves the inversion of a triangular matrix which projects the mass power spectrum in 3-D real-space to the transmission power spectrum in 1-D redshift-space. We illustrate the method by performing a linear calculation relating the two. A method that does not take into account redshift-space anisotropy tends to underestimate the steepness of the mass power spectrum, in the case of linear distortions. The issue of the effective bias-factor for the linear distortion kernel is discussed.

One of the major goals of cosmological observations is to test theories of structure formation. The most straightforward way to carry out such tests is to compute the likelihood function L, the probability of getting the data given the theory. We write down this function for a general galaxy survey. The full likelihood function is very complex, depending on all of the $n$-point functions of the theory under consideration. Even in the simplest case, where only the two point function is non-vanishing (Gaussian perturbations), L cannot be calculated exactly, primarily because of the Poisson nature of the galaxy distribution. Here we expand L about the (trivial) zero correlation limit. As a first application, we take the binned values of the two point function as free parameters and show that L peaks at $(DD - DR + RR)/DD$. Using Monte Carlo techniques, we compare this estimator with the traditional $DD/DR$ and Landy & Szalay estimators. More generally, the success of this expansion should pave the way for further applications of the likelihood function.

We review the arguments for how the slope of the column density distribution of the Lyman-alpha forest should depend on the matter power spectrum. The latest progress, presented by various groups in this conference and elsewhere, is summarized.

A theory is developed which relates the observed b-parameter of a Lyman-alpha absorption line to the velocity-curvature of the corresponding peak in the optical depth fluctuation. Its relation to the traditional interpretation of b as the thermal broadening width is discussed. It is demonstrated that, independent of the details of the cosmological model, the differential b distribution has a high b asymptote of $dN/db \propto b^{-m}$, where $m \geq 5$, when we make the reasonable assumption that low-curvature fluctuations are statistically favored over high-curvature ones. There in general always exist absorption lines much broader than the thermal width. We then develop a linear perturbative analysis of the optical depth fluctuation, which yields a single-parameter prediction for the full b distribution: in addition to exhibiting the general high velocity tail, it explains the observed sharp low b cut-off. The dependence of the b distribution on cosmological parameters, such as $\Omega$ and the power spectrum, and reionization history as well as observation/simulation resolution is derived and discussed.

We develop a method to calculate the column density distribution of the Lyman-alpha forest for column densities in the range $10^{12.5} - 10^{14.5} cm^{-2}$. The Zel'dovich approximation, with appropriate smoothing, is used to compute the density and peculiar velocity fields. The effect of the latter on absorption profiles is discussed and it is shown to have little effect on the column density distribution. An approximation is introduced in which the column density distribution is related to a statistic of density peaks (involving its height and first and second derivatives along the line of sight) in real space. We show that the slope of the column density distribution is determined by the temperature-density relation as well as the power spectrum on scales $2 h Mpc^{-1} < k < 20 h Mpc^{-1}$. An expression relating the three is given. We find very good agreement between the column density distribution obtained by applying the Voigt-profile-fitting technique to the output of a full hydrodynamic simulation and that obtained using our approximate method for a test model. This formalism then is applied to study a group of CDM as well as CHDM models. We show that the amplitude of the column density distribution depends on the combination of parameters $(\Omega_b h^2)^2 T_0^{-0.7} J_{HI}^{-1}$, which is not well-constrained by independent observations. The slope of the distribution, on the other hand, can be used to distinguish between different models: those with a smaller amplitude and a steeper slope of the power spectrum on small scales give rise to steeper distributions, for the range of column densities we study. Comparison with high resolution Keck data is made.

We introduce an efficient and accurate alternative to full hydrodynamic simulations, Hydro-PM (HPM), for the study of the low column density Lyman-alpha forest. It consists of a Particle-Mesh solver, modified to compute, in addition to the gravitational potential, an effective potential due to the gas pressure. Such an effective potential can be computed from the density field because of a tight correlation between density and pressure in the low density limit, which can be calculated for any photo-reionization history by a method outlined in Hui & Gnedin (1997). Such a correlation exists, in part, because of minimal shock-heating in the low density limit. We compare carefully the density and velocity fields as well as absorption spectra, computed using HPM versus hydrodynamic simulations, and find good agreement. We show that HPM is capable of reproducing measurable quantities, such as the column density distribution, computed from full hydrodynamic simulations, to a precision comparable to that of observations. We discuss how, by virtue of its speed and accuracy, HPM can enable us to use the Lyman-alpha forest as a cosmological probe. We also discuss in detail the smoothing of the gas (or baryon) fluctuation relative to that of the dark matter on small scales due to finite gas pressure. It is shown the conventional wisdom that the linear gas fluctuation is smoothed on the Jeans scale is incorrect for general reionization (or reheating) history; the correct linear filtering scale is in general smaller than the Jeans scale after reheating, but larger prior to it. (abridged)

Absorption lines in the Lyman-alpha forest can be thought of as peaks in neutral hydrogen density along lines of sight. The column density distribution (the number density of absorption lines as a function of column density) is then a statistic of density peaks, which contains information about the underlying power spectrum. In particular, we show that the slope of the distribution provides a measure of power on scales smaller than those probed by studies of present-day large scale structure.

We develop an efficient method to study the effects of reionization history on the temperature-density relation of the intergalactic medium in the low density limit (overdensity less than 5). It is applied to the study of photo-reionization models in which the amplitude, spectrum and onset epoch of the ionizing flux, as well as the cosmology, are systematically varied. We find that the mean temperature-density relation at z=2-4 is well approximated by a power-law equation of state for uniform reionization models. We derive analytical expressions for its evolution and exhibit its asymptotic behavior: it is found that for sufficiently early reionization, imprints of reionization history prior to z=10 on the temperature-density relation are washed out. In this limit the temperature at cosmic mean density is proportional to (\Omega_b h/\sqrt\Omega_0)^{1/1.7}. While the amplitude of the radiation flux at the ionizing frequency of HI is found to have a negligible effect on the temperature-density relation as long as the universe reionizes before z=5, the spectrum can change the overall temperature by about 20%, through variations in the abundances of helium species. However the slope of the mean equation of state is found to lie within a narrow range for all reionization models we study, where reionization takes place before z=5. We discuss the implications of these findings for the observational properties of the Lyman-alpha forest. In particular, uncertainties in the temperature of the intergalactic medium, due to the uncertain reionization history of our universe, introduces a 30% scaling in the amplitude of the column density distribution while the the slope of the distribution is only affected by about 5%. Finally, we discuss how a fluctuating ionizing field affects the above results. We argue that under

We develop an analytical method based on the lognormal approximation to compute the column density distribution of the Lyman-alpha forest in the low column density limit. We compute the column density distributions for six different cosmological models and found that the standard, COBE-normalized CDM model cannot fit the observations of the Lyman-alpha forest at z=3. The amplitude of the fluctuations in that model has to be lowered by a factor of almost 3 to match observations. However, the currently viable cosmological models like the lightly tilted COBE-normalized CDM+Lambda model, the CHDM model with 20% neutrinos, and the low-amplitude Standard CDM model are all in agreement with observations, to within the accuracy of our approximation, for the value of the cosmological baryon density at or higher than the old Standard Bing Bang Nucleosynthesis value of 0.0125 for the currently favored value of the ionizing radiation intensity. With the low value for the baryon density inferred by Hogan & Rugers (1996), the models can only marginally match observations.

The observed 2B BATSE distribution is consistent with the faintest GRBs in our sample originating from a redshift of Zmax ~ 0.8-3.0 (90\%), with the most likely values in the range of 1.0-2.2, and is largely insensitive to Omega for models with no evolution. To constrain the model parameter Omega to the range 0.1-1.0 using only Log N -- Log P distributions, more than 4000 GRBs, with a most likely value of ~ 9,000 GRBs to BATSE sensitivity. This requires a live integration time of >6 years with BATSE. Detectors sensitive to much lower limits (~ 70-400 in sensitivity) require ~ 200 GRBs, with <0.03 year 4pi ster coverage. We place limits on the amount of frequency density and, separately, peak luminosity evolution in the sample of GRBs. We find that frequency density evolution models can place the faintest GRBs at Zmax ~ 10-200, without conflicting with the observations of relative time dilation of ~ 2.

We investigate three different local approximations for nonlinear gravitational instability in the framework of cosmological Lagrangian fluid dynamics of cold dust. They include the Zel'dovich approximation (ZA), the ``non-magnetic'' approximation of Bertschinger \& Jain (1994, NMA), and a new ``local tidal'' approximation (LTA). The LTA is exact for any perturbations whose gravitational and velocity equipotentials have the same constant shape with time, including spherical, cylindrical, and plane-parallel perturbations. We tested all three local approximations with the collapse of a homogeneous triaxial ellipsoid, for which an exact solution exists for an ellipsoid embedded in empty space and an excellent approximation is known in the cosmological context. We find that the LTA is significantly more accurate in general than the ZA and the NMA. Like the ZA, but unlike the NMA, the LTA generically leads to pancake collapse. For a randomly chosen mass element in an Einstein-de Sitter universe, assuming a Gaussian random field of initial density fluctuations, the LTA predicts that at least 78\% of initially underdense regions collapse owing to nonlinear effects of shear and tides.