The standard description of cosmological observables is incomplete, because it does not take into account the correct angular parametrization of the sky, i.e. the one determined by the observer frame. The corresponding corrections must be taken into account for reliable results at non-linear orders. This can be accomplished by introducing an orthonormal basis, or "tetrad", at the observer point, representing the frame with respect to which observations are performed. In this work we consider the tetrad formulation of General Relativity, thus associating tetrads to sources as well, and develop a new formalism for describing cosmological observables associated with localized sources. It is based on a manifold which we call the "observer space-time", whose coordinates are the proper time, redshift and angles an observer uses to parametrize measurements, and on which the rest of the observables are defined. This manifold does not have to be diffeomorphic to the true space-time and allows us to resolve caustics in the latter, in contrast to similar coordinate-based formalisms. As a concrete example, we work out the definitions and equations for the angular diameter distance, weak lensing and number count observables. As for the observables associated to the CMB, they lie inside the phase space distribution of the photon fluid, so we also revisit the construction of general-relativistic matrix kinetic theory from the tetrad formalism viewpoint. Here too the latter appears as the natural approach for relating the macroscopic dynamics to the microscopic quantum field theory, and therefore for constructing the matrix Boltzmann equations, without any approximation on the gravitational side. (... more in the manuscript)

This paper focuses on two aspects of the statistics of cosmological observables that are important for the next stages of precision cosmology. First, we note that the theory of reduced angular $N$-point spectra has only been developed in detail up to the trispectrum case and in a fashion that makes it difficult to go beyond. To fill this gap, here we present a constructive approach that provides a systematic description of reduced angular $N$-point spectra and their covariance matrices, for arbitrary $N$. Second, we focus on the common practice in the literature on cosmological observables, which consists in simply discarding a part of the expression, namely, the terms containing fields evaluated at the observer position. We point out that this is not justified beyond linear order in perturbation theory, as these terms contribute to all the multipoles of the corresponding spectra and with a magnitude that is of the same order as the rest of the non-linear corrections. We consider the possibility that the reason for neglecting these terms is a conceptual discomfort when using ensemble averages, which originates in an apparent tension between the ergodic hypothesis and the privileged position of the observer on the light-cone. We clarify this subtle issue by performing a careful derivation of the relation between the theoretical statistical predictions and the observational estimators for all $N$. We conclude that there is no inconsistency whatsoever in ensemble-averaging fields at and near the observer position, thus clearing the way for consistent and robust high-precision calculations.

We present the galaxy power spectrum in general relativity. Using a novel approach to computing the power spectrum, we derive the galaxy power spectrum taking into account all the relativistic effects in observations. We first demonstrate that the presence of an infrared divergence in previous calculations of the galaxy power spectrum is inconsistent with the equivalence principle, and we show that this unphysical behavior on large scales vanishes if we consider all the relativistic effects consistently. As a consequence, the relativistic effects in the galaxy power spectrum yield no corrections to the measurement of the primordial non-Gaussianity. We present our numerical computation of the full galaxy power spectrum, with particular emphasis on the deviations from the standard redshift-space power spectrum on large scales. We conclude that, as relativistic effects significantly alter the galaxy power spectrum at $k\lesssim k_{eq}$, they need to be taken into account in the analysis of large-scale data.

We present the first statistical study of magnetic structures and associated energy dissipation observed during a single period of turbulent magnetic reconnection, by using the \textit{in-situ} measurements of the Magnetospheric Multiscale mission in the Earth's magnetotail on July 26, 2017. The structures are selected by identifying a bipolar signature in the magnetic field and categorized as plasmoids or current sheets via an automated algorithm which examines current density and plasma flow. The size of the plasmoids forms a decaying exponential distribution ranging from sub-electron up to ion scales. The presence of substantial number of current sheets is consistent with a physical picture of dynamic production and merging of plasmoids during turbulent reconnection. The magnetic structures are locations of significant energy dissipation via electric field parallel to the local magnetic field, while dissipation via perpendicular electric field dominates outside of the structures. Significant energy also returns from particles to fields.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

Type-Ia supernova observations yield estimates of the luminosity distance, which includes not only the background luminosity distance, but also the fluctuation due to inhomogeneities in the Universe. In particular, the spatial correlation of the host galaxies is a dominant source of the fluctuation in the luminosity distance measurements. Utilizing the recent theoretical framework that accurately quantifies the information contents accounting for the three-dimensional correlation of the observables on the past-light cone, we compute the maximum cosmological information obtainable from idealized supernova surveys as a function of maximum redshift $z_m$. Here we consider two cosmological parameters $\Omega_m$ and $w_0$ and show that these parameters can be constrained at maximum 1% levels in an idealized survey with $z_m=1$. We discuss how these fundamental limits set by cosmic variance can be overcome.

The background photon temperature $\bar T$ is one of the fundamental cosmological parameters. Despite its significance, $\bar T$ has never been allowed to vary in the data analysis, owing to the precise measurement of the comic microwave background (CMB) temperature by COBE FIRAS. However, even in future CMB experiments, $\bar T$ will remain unknown due to the unknown monopole contribution $\Theta_0$ at our position to the observed (angle-averaged) temperature $\langle T\rangle^{\rm obs}$. By fixing $\bar T\equiv\langle T\rangle^{\rm obs}$, the standard analysis underestimates the error bars on cosmological parameters, and the best-fit parameters obtained in the analysis are biased in proportion to the unknown amplitude of $\Theta_0$. Using the Fisher formalism, we find that these systematic errors are smaller than the error bars from the $Planck$ satellite. However, with $\bar T\equiv\langle T\rangle^{\rm obs}$, these systematic errors will always be present and irreducible, and future cosmological surveys might misinterpret the measurements.

We develop a theoretical framework to describe the cosmological observables on the past light cone such as the luminosity distance, weak lensing, galaxy clustering, and the cosmic microwave background anisotropies. We consider that all the cosmological observables include not only the background quantity, but also the perturbation quantity, and they are subject to cosmic variance, which sets the fundamental limits on the cosmological information that can be derived from such observables, even in an idealized survey with an infinite number of observations. To quantify the maximum cosmological information content, we apply the Fisher information matrix formalism and spherical harmonic analysis to cosmological observations, in which the angular and the radial positions of the observables on the light cone carry different information. We discuss the maximum cosmological information that can be derived from five different observables: (1) type Ia supernovae, (2) cosmic microwave background anisotropies, (3) weak gravitational lensing, (4) local baryon density, and (5) galaxy clustering. We compare our results with the cosmic variance obtained in the standard approaches, which treat the light cone volume as a cubic box of simultaneity. We discuss implications of our formalism and ways to overcome the fundamental limit.

We perform theoretical and numerical studies of the full relativistic two-point galaxy correlation function, considering the linear-order scalar and tensor perturbation contributions and the wide-angle effects. Using the gauge-invariant relativistic description of galaxy clustering and accounting for the contributions at the observer position, we demonstrate that the complete theoretical expression is devoid of any long-mode contributions from scalar or tensor perturbations and it lacks the infrared divergences in agreement with the equivalence principle. By showing that the gravitational potential contribution to the correlation function converges in the infrared, our study justifies an IR cut-off $(k_{\text{IR}} \leq H_0)$ in computing the gravitational potential contribution. Using the full gauge-invariant expression, we numerically compute the galaxy two-point correlation function and study the individual contributions in the conformal Newtonian gauge. We find that the terms at the observer position such as the coordinate lapses and the observer velocity (missing in the standard formalism) dominate over the other relativistic contributions in the conformal Newtonian gauge such as the source velocity, the gravitational potential, the integrated Sachs-Wolf effect, the Shapiro time-delay and the lensing convergence. Compared to the standard Newtonian theoretical predictions that consider only the density fluctuation and redshift-space distortions, the relativistic effects in galaxy clustering result in a few percent-level systematic errors beyond the scale of the baryonic acoustic oscillation. Our theoretical and numerical study provides a comprehensive understanding of the relativistic effects in the galaxy two-point correlation function, as it proves the validity of the theoretical prediction and accounts for effects that are often neglected in its numerical evaluation.

We present the second-order expression for the observed redshift, accounting for all the relativistic effects from the light propagation and from the frame change at the observer and the source positions. We derive the generic gauge-transformation law that any observable quantities should satisfy, and we verify our second-order expression for the observed redshift by explicitly checking its gauge transformation property. This is the first time an explicit verification is made for the second-order calculations of observable quantities. We present our results in popular gauge choices for easy use and discuss the origin of disagreements in previous calculations.

Cosmological weak lensing has been a highly successful and rapidly developing research field since the first detection of cosmic shear in 2000. However, it has recently been pointed out in Yoo et al. that the standard weak lensing formalism yields gauge-dependent results and, hence, does not meet the level of accuracy demanded by the next generation of weak lensing surveys. Here, we show that the Jacobi mapping formalism provides a solid alternative to the standard formalism, as it accurately describes all the relativistic effects contributing to the weak lensing observables. We calculate gauge-invariant expressions for the distortion in the luminosity distance, the cosmic shear components and the lensing rotation to linear order including scalar, vector and tensor perturbations. In particular, the Jacobi mapping formalism proves that the rotation is fully vanishing to linear order. Furthermore, the cosmic shear components contain an additional term in tensor modes which is absent in the results obtained with the standard formalism. Our work provides further support and confirmation of the gauge-invariant lensing formalism needed in the era of precision cosmology.

We present the gauge-invariant formalism of cosmological weak lensing, accounting for all the relativistic effects due to the scalar, vector, and tensor perturbations at the linear order. While the light propagation is fully described by the geodesic equation, the relation of the photon wavevector to the physical quantities requires the specification of the frames, where they are defined. By constructing the local tetrad bases at the observer and the source positions, we clarify the relation of the weak lensing observables such as the convergence, the shear, and the rotation to the physical size and shape defined in the source rest-frame and the observed angle and redshift measured in the observer rest-frame. Compared to the standard lensing formalism, additional relativistic effects contribute to all the lensing observables. We explicitly verify the gauge-invariance of the lensing observables and compare our results to previous work. In particular, we demonstrate that even in the presence of the vector and tensor perturbations, the physical rotation of the lensing observables vanishes at the linear order, while the tetrad basis rotates along the light propagation compared to a FRW coordinate. Though the latter is often used as a probe of primordial gravitational waves, the rotation of the tetrad basis is indeed not a physical observable. We further clarify its relation to the E-B decomposition in weak lensing. Our formalism provides a transparent and comprehensive perspective of cosmological weak lensing.

Baryon Acoustic Oscillations (BAOs) in the early Universe are predicted to leave an as yet undetected signature on the relative clustering of total mass versus luminous matter. This signature, a modulation of the relative large-scale clustering of baryons and dark matter, offers a new angle to compare the large scale distribution of light versus mass. A detection of this effect would provide an important confirmation of the standard cosmological paradigm and constrain alternatives to dark matter as well as non-standard fluctuations such as Compensated Isocurvature Perturbations (CIPs). The first attempt to measure this effect in the SDSS-III BOSS Data Release 10 CMASS sample remained inconclusive but allowed to develop a method, which we detail here and use to conduct the second observational search. When using the same model as in our previous study and including CIPs in the model, the DR12 data are consistent with a null-detection, a result in tension with the strong evidence previously measured with the DR10 data. This tension remains when we use a more realistic model taking into account our knowledge of the survey flux limit, as the data then privilege a zero effect. In the absence of CIPs, we obtain a null detection consistent with both the absence of the effect and the amplitude predicted in previous theoretical studies. This shows the necessity of more accurate data in order to prove or disprove the theoretical predictions.

We present a complete set of exact and fully non-linear equations describing all three types of cosmological perturbations -- scalar, vector and tensor perturbations. We derive the equations in a thoroughly gauge-ready manner, so that any spatial and temporal gauge conditions can be employed. The equations are completely general without any physical restriction except that we assume a flat homogeneous and isotropic universe as a background. We also comment briefly on the application of our formulation to the non-expanding Minkowski background.

We present the correlation function of the luminosity distances in a flat $\Lambda$CDM universe. Decomposing the luminosity distance fluctuation into the velocity, the gravitational potential, and the lensing contributions, we study their individual contributions to the correlation function. The lensing contribution is important at large redshift ($z\gtrsim 0.5$) but only for small angular separation ($\theta \lesssim 3^\circ$), while the velocity contribution dominates over the other contributions at low redshift or at larger separation. However, the gravitational potential contribution is always subdominant at all scale, if the correct gauge-invariant expression is used. The correlation function of the luminosity distances depends significantly on the matter content, especially for the lensing contribution, thus providing a novel tool of estimating cosmological parameters.

Theoretical descriptions of observable quantities in cosmological perturbation theory should be independent of coordinate systems. This statement is often referred to as gauge-invariance of observable quantities, and the sanity of their theoretical description is verified by checking its gauge-invariance. We argue that cosmological observables are invariant scalars under diffeomorphisms and as a consequence their theoretical description is gauge-invariant, only at linear order in perturbations. Beyond linear order, they are usually not gauge-invariant, and we provide the general law for the gauge-transformation that the perturbation part of an observable does obey. We apply this finding to derive the second-order expression for the observational light-cone average in cosmology and demonstrate that our expression is indeed invariant under diffeomorphisms.

The remarkable properties of the geodesic light-cone (GLC) coordinates allow analytic expressions for the light-cone observables, providing a new non-perturbative way for calculating the effects of inhomogeneities in our Universe. However, the gauge-invariance of these expressions in the GLC formalism has not been shown explicitly. Here we provide this missing part of the GLC formalism by proving the gauge-invariance of the GLC expressions for the light-cone observables, such as the observed redshift, the luminosity distance, and the physical area and volume of the observed sources. Our study provides a new insight on the properties of the GLC coordinates and it complements the previous work by the GLC collaboration, leading to a comprehensive description of light propagation in the GLC representation.

Measurements of the luminosity distance have played a key role in discovering the late-time cosmic acceleration. However, when accounting for inhomogeneities in the Universe, its interpretation has been plagued with infrared divergences in its theoretical predictions, which are in some cases used to explain the cosmic acceleration without dark energy. The divergences in most calculations are artificially removed by imposing an infrared cut-off scale. For the first time, we show that a gauge-invariant calculation of the luminosity distance is devoid of such divergences and consistent with the equivalence principle, eliminating the need to impose a cut-off scale. We present proper numerical calculations of the luminosity distance using the gauge-invariant expression and demonstrate that the numerical results with an ad hoc cut-off scale in previous calculations have negligible systematic errors as long as the cut-off scale is larger than the horizon scale. We discuss the origin of infrared divergences and their cancellation in the luminosity distance.

Comparing the luminosity distance measurements to its theoretical predictions is one of the cornerstones in establishing the modern cosmology. However, as shown in Biern & Yoo, its theoretical predictions in literature are often plagued with infrared divergences and gauge-dependences. This trend calls into question the sanity of the methods used to derive the luminosity distance. Here we critically investigate four different methods --- the geometric approach, the Sachs approach, the Jacobi mapping approach, and the geodesic light cone (GLC) approach to modeling the luminosity distance, and we present a unified treatment of such methods, facilitating the comparison among the methods and checking their sanity. All of these four methods, if exercised properly, can be used to reproduce the correct description of the luminosity distance.

We derive the exact third-order analytic solution of the matter density fluctuation in the proper-time hypersurface in a $\Lambda$CDM universe, accounting for the explicit time-dependence and clarifying the relation to the initial condition. Furthermore, we compare our analytic solution to the previous calculation in the comoving gauge, and to the standard Newtonian perturbation theory by providing Fourier kernels for the relativistic effects. Our results provide an essential ingredient for a complete description of galaxy bias in the relativistic context.

Liquid argon is a bright scintillator with potent particle identification properties, making it an attractive target for direct-detection dark matter searches. The DarkSide-50 dark matter search here reports the first WIMP search results obtained using a target of low-radioactivity argon. DarkSide-50 is a dark matter detector, using two-phase liquid argon time projection chamber, located at the Laboratori Nazionali del Gran Sasso. The underground argon is shown to contain Ar-39 at a level reduced by a factor (1.4 +- 0.2) x 10^3 relative to atmospheric argon. We report a background-free null result from (2616 +- 43) kg d of data, accumulated over 70.9 live-days. When combined with our previous search using an atmospheric argon, the 90 % C.L. upper limit on the WIMP-nucleon spin-independent cross section based on zero events found in the WIMP search regions, is 2.0 x 10^-44 cm^2 (8.6 x 10^-44 cm^2, 8.0 x 10^-43 cm^2) for a WIMP mass of 100 GeV/c^2 (1 TeV/c^2 , 10 TeV/c^2).

We present the third-order analytic solution of the matter density fluctuation in the proper-time hypersurface of nonrelativistic matter flows by solving the nonlinear general relativistic equations. The proper-time hypersurface provides a coordinate system that a local observer can set up without knowledge beyond its neighborhood, along with physical connections to the local Newtonian descriptions in the relativistic context. The initial condition of our analytic solution is set up by the curvature perturbation in the comoving gauge, clarifying its impact on the nonlinear evolution. We compute the effective non-Gaussian parameters due to the nonlinearity in the relativistic equations. With proper coordinate rescaling, we show that the equivalence principle is respected and the relativistic effect vanishes in the large-scale limit.

We report the first results of DarkSide-50, a direct search for dark matter operating in the underground Laboratori Nazionali del Gran Sasso (LNGS) and searching for the rare nuclear recoils possibly induced by weakly interacting massive particles (WIMPs). The dark matter detector is a Liquid Argon Time Projection Chamber with a (46.4+-0.7) kg active mass, operated inside a 30 t organic liquid scintillator neutron veto, which is in turn installed at the center of a 1 kt water Cherenkov veto for the residual flux of cosmic rays. We report here the null results of a dark matter search for a (1422+-67) kg d exposure with an atmospheric argon fill. This is the most sensitive dark matter search performed with an argon target, corresponding to a 90% CL upper limit on the WIMP-nucleon spin-independent cross section of 6.1x10^-44 cm^2 for a WIMP mass of 100 GeV/c^2.

It is generally inferred from astronomical measurements that Dark Matter (DM) comprises approximately 27\% of the energy-density of the universe. If DM is a subatomic particle, a possible candidate is a Weakly Interacting Massive Particle (WIMP), and the DarkSide-50 (DS) experiment is a direct search for evidence of WIMP-nuclear collisions. DS is located underground at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy, and consists of three active, embedded components; an outer water veto (CTF), a liquid scintillator veto (LSV), and a liquid argon (LAr) time projection chamber (TPC). This paper describes the data acquisition and electronic systems of the DS detectors, designed to detect the residual ionization from such collisions.

Current and future galaxy surveys cover a large fraction of the entire sky with a significant redshift range, and the recent theoretical development shows that general relativistic effects are present in galaxy clustering on very large scales. This trend has renewed interest in the wide angle effect in galaxy clustering measurements, in which the distant-observer approximation is often adopted. Using the full wide-angle formula for computing the redshift-space correlation function, we show that compared to the sample variance, the deviation in the redshift-space correlation function from the simple Kaiser formula with the distant-observer approximation is negligible in galaxy surveys such as the SDSS, Euclid and the BigBOSS, if the theoretical prediction from the Kaiser formula is properly averaged over the survey volume. We also find corrections to the wide-angle formula and clarify the confusion in literature between the wide angle effect and the velocity contribution in galaxy clustering. However, when the FKP method is applied, substantial deviations can be present in the power spectrum analysis in future surveys, due to the non-uniform distribution of galaxy pairs.

We compute the second-order density fluctuation in the proper-time hypersurface of non-relativistic matter flows and relate it to the galaxy number density fluctuation in general relativity. At the linear order, it is equivalent to the density fluctuation in the comoving synchronous gauge, in which two separate gauge conditions coincide. However, at the second order, the density fluctuations in these gauge conditions differ, while both gauge conditions represent the proper-time hypersurface. Compared to the density fluctuation in the temporal comoving and the spatial C-gauge conditions, the density fluctuation in the commonly used gauge condition ($N=1$ and $N^\alpha=0$) violates the mass conservation at the second order. We provide their physical interpretations in each gauge condition by solving the geodesic equation and the nonlinear evolution equations of non-relativistic matter. We apply this finding to the second-order galaxy biasing in general relativity, which complements the second-order relativistic description of galaxy clustering in Yoo & Zaldarriaga (2014).

The general relativistic description of galaxy clustering provides a complete and unified treatment of all the effects in galaxy clustering such as the redshift-space distortion, gravitational lensing, Sachs-Wolfe effects, and their relativistic effects. In particular, the relativistic description resolves the gauge issues in the standard Newtonian description of galaxy clustering by providing the gauge-invariant expression for the observed galaxy number density. The relativistic effect in galaxy clustering is significant on large scales, in which dark energy models or alternative theories of modified gravity deviate from general relativity. In this paper, we review the relativistic effect in galaxy clustering by providing a pedagogical derivation of the relativistic formula and by computing the observed galaxy two-point statistics. The relativistic description of galaxy clustering is an essential tool for testing general relativity and probing the early Universe on large scales in the era of precision cosmology.

We present the second-order general relativistic description of the observed galaxy number density in a cosmological framework. The observed galaxy number density is affected by the volume and the source effects, both of which arise due to the mismatch between physical and observationally inferred quantities such as the redshift, the angular position, the volume, and the luminosity of the observed galaxies. These effects are computed to the second order in metric perturbations without choosing a gauge condition or adopting any restrictions on vector and tensor perturbations, extending the previous linear-order calculations. Paying particular attention to the second-order gauge transformation, we explicitly isolate unphysical gauge modes and construct second-order gauge-invariant variables. Moreover, by constructing second-order tetrads in the observer's rest frame, we clarify the relation between the physical and the parametrized photon wavevectors. Our second-order relativistic description will provide an essential tool for going beyond the power spectrum in the era of precision measurements of galaxy clustering. We discuss potential applications and extensions of the second-order relativistic description of galaxy clustering.

We study the effect of the supersonic relative velocity between dark matter and baryons on large-scale galaxy clustering and derive the constraint on the relative velocity bias parameter from the Baryonic Oscillation Spectroscopic Survey (BOSS) power spectrum measurements. Recent work has shown that the relative velocity effect may have a dramatic impact on the star formation at high redshifts, if first stars are formed in minihalos around z~20, or if the effect propagates through secondary effects to stars formed at later redshifts. The relative velocity effect has particularly strong signatures in the large scale clustering of these sources, including the BAO position. Assuming that a small fraction of stars in low-redshift massive galaxies retain the memory of the primordial relative velocity effect, galaxy clustering measurements can be used to constrain the signatures of the first stars. Luminous red galaxies contain some of the oldest stars in the Universe and are ideally suited to search for this effect. Using the BOSS power spectrum measurements from the Sloan Data Release 9, in combination with Planck, we derive the upper limit on the fraction of the stars sensitive to relative velocity effect f_star<3.3% at the 95% confidence level in the CMASS galaxy sample. If additional galaxy sample not sensitive to the effect is available in a given survey, a joint multi-tracer analysis can be applied to construct a sample-variance cancelling combination, providing a model-independent way to verify the presence of the relative velocity effect in the galaxy power spectrum on large scales. Such a multi-tracer analysis in future galaxy surveys can greatly improve the current constraint, achieving a 0.1% level in f_star.

In this Report we discuss the four complementary searches for the identity of dark matter: direct detection experiments that look for dark matter interacting in the lab, indirect detection experiments that connect lab signals to dark matter in our own and other galaxies, collider experiments that elucidate the particle properties of dark matter, and astrophysical probes sensitive to non-gravitational interactions of dark matter. The complementarity among the different dark matter searches is discussed qualitatively and illustrated quantitatively in several theoretical scenarios. Our primary conclusion is that the diversity of possible dark matter candidates requires a balanced program based on all four of those approaches.

We explore the signatures of quintessence and modified gravity theories in the relativistic description of galaxy clustering within a parametrized post-Friedmann framework. For this purpose, we develop a calibration method to consistently account for horizon-scale effects in the linear parametrized Post-Friedmann perturbations of minimally and nonminimally coupled scalar-tensor theories and test it against the full model-specific fluctuations. We further study the relativistic effects in galaxy clustering for the normal and self-accelerating branches of the Dvali-Gabadadze-Porrati braneworld model as well as for phenomenological modifications of gravity. We quantify the impact of modified gravity and dark energy models on galaxy clustering by computing the velocity-to-matter density ratio F, the velocity contribution R, and the potential contribution P and give an estimate of their detectability in future galaxy surveys. Our results show that, in general, the relativistic correction contains additional information on gravity and dark energy, which needs to be taken into account in consistent horizon-scale tests of departures from LCDM using the galaxy-density field.

We report results of a search for Weakly Interacting Massive Particles (WIMPS) with the silicon detectors of the CDMS II experiment. This blind analysis of 140.2 kg-days of data taken between July 2007 and September 2008 revealed three WIMP-candidate events with a surface-event background estimate of 0.41^{+0.20}_{-0.08}(stat.)^{+0.28}_{-0.24}(syst.). Other known backgrounds from neutrons and 206Pb are limited to < 0.13 and <0.08 events at the 90% confidence level, respectively. The exposure of this analysis is equivalent to 23.4 kg-days for a recoil energy range of 7-100 keV for a WIMP of mass 10 GeV/c2. The probability that the known backgrounds would produce three or more events in the signal region is 5.4%. A profile likelihood ratio test of the three events that includes the measured recoil energies gives a 0.19% probability for the known-background-only hypothesis when tested against the alternative WIMP+background hypothesis. The highest likelihood occurs for a WIMP mass of 8.6 GeV/c2 and WIMP-nucleon cross section of 1.9e-41 cm2.

We report results of a search for Weakly Interacting Massive Particles (WIMPs) with the Si detectors of the CDMS II experiment. This report describes a blind analysis of the first data taken with CDMS II's full complement of detectors in 2006-2007; results from this exposure using the Ge detectors have already been presented. We observed no candidate WIMP-scattering events in an exposure of 55.9 kg-days before analysis cuts, with an expected background of ~1.1 events. The exposure of this analysis is equivalent to 10.3 kg-days over a recoil energy range of 7-100 keV for an ideal Si detector and a WIMP mass of 10 GeV/c2. These data set an upper limit of 1.7x10-41 cm2 on the WIMP-nucleon spin-independent cross section of a 10 GeV/c2 WIMP. These data exclude parameter space for spin-independent WIMP-nucleon elastic scattering that is relevant to recent searches for low-mass WIMPs.

We perform an all-sky analysis of the general relativistic galaxy power spectrum using the well-developed spherical Fourier decomposition. Spherical Fourier analysis expresses the observed galaxy fluctuation in terms of the spherical harmonics and spherical Bessel functions that are angular and radial eigenfunctions of the Helmholtz equation, providing a natural orthogonal basis for all-sky analysis of the large-scale mode measurements. Accounting for all the relativistic effects in galaxy clustering, we compute the spherical power spectrum and its covariance matrix and compare it to the standard three-dimensional power spectrum to establish a connection. The spherical power spectrum recovers the three-dimensional power spectrum at each wavenumber k with its angular dependence mu_k encoded in angular multipole l, and the contributions of the line-of-sight projection to galaxy clustering such as the gravitational lensing effect can be readily accommodated in the spherical Fourier analysis. A complete list of formulas for computing the relativistic spherical galaxy power spectrum is also presented.

Mounting observational data confirm that about 73% of the energy density consists of dark energy which is responsible for the current accelerated expansion of the Universe. We present observational evidences and dark energy projects. We then review various theoretical ideas that have been proposed to explain the origin of dark energy; they contain the cosmological constant, modified matter models, modified gravity models and the inhomogeneous model. The cosmological constant suffers from two major problems: one regarding fine-tuning and the other regarding coincidence. To solve them there arose modified matter models such as quintessence, k-essence, coupled dark energy, and unified dark energy. We compare those models by presenting attractive aspects, new rising problems and possible solutions. Furthermore we review modified gravity models that lead to late-time accelerated expansion without invoking a new form of dark energy; they contain f(R) gravity and the Dvali-Gabadadze-Porrati model. We also discuss observational constraints on those models and on future modified gravity theories. Finally we review the inhomogeneous Lemaitre-Tolman-Bondi model that drops an assumption of the spatial homogeneity of the Universe. We also present basics of cosmology and scalar field theory, which are useful especially for students and novices to understand dark energy models.

The proper general relativistic description of the observed galaxy power spectrum is substantially different from the standard Newtonian description on large scales, providing a unique opportunity to test general relativity on horizon scales. Using the Einstein equations, the general relativistic effects can be classified as two new terms that represent the velocity and the gravitational potential, coupling to the time evolution of galaxy number density and Hubble parameter. Compared to the dominant density and velocity redshift-space distortion terms, the former scales as H/k and correlates the real and imaginary parts of the Fourier modes, while the latter scales as (H/k)^2, where k is the comoving wave number and H is the conformal Hubble parameter. We use the recently developed methods to reduce the sampling variance and shot noise to show that in an all sky galaxy redshift survey at low redshift the velocity term can be measured at 10-sigma confidence level, if one can utilize halos of mass M>10^{10} Msun, while the gravitational potential term itself can only be marginally detected. We also demonstrate that the general relativistic effect is not degenerate with the primordial non-Gaussian signature in galaxy bias, and the ability to detect the primordial non-Gaussianity is little compromised.

We explore three different methods based on weak lensing to extract cosmological constraints from the large-scale structure. In the first approach (method I), small-scale galaxy lensing measurements of their halo mass provide a constraint on the halo bias, which can be combined with the large-scale galaxy clustering to measure the dark matter clustering. In the second approach (method II), large-scale galaxy clustering and large-scale galaxy-galaxy lensing can be combined into a direct measurement of the dark matter clustering. These two methods can be combined into one method I+II to make use of lensing measurements on all scales. In the third approach (method III), we add abundance information to the method I. We explore the statistical power of these three approaches as a function of galaxy luminosity to investigate the optimal mass range for each method and their cosmological constraining power. In the case of the SDSS, we find that the three methods give comparable constraints, but not in the same mass range: the method II works best for halos of M~10^13 Msun, and the methods I and III work best for halos of M~10^14 Msun. We discuss the robustness of each method against various systematics. Furthermore, we extend the analysis to the future large-scale galaxy surveys and find that the cluster abundance method is not superior to the combined method I+II, both in terms of statistical power and robustness against systematic errors. The cosmic shear-shear correlation analysis in the future surveys yields constraints as strong as the combined method, but suffer from additional systematic effects. We thus advocate the combined analysis of clustering and lensing (method I+II) as a powerful alternative to other large-scale probes. Our analysis provides a guidance to observers planning large-scale galaxy surveys such as the DES, Euclid, and the LSST.

We report limits on annual modulation of the low-energy event rate from the Cryogenic Dark Matter Search (CDMS II) experiment at the Soudan Underground Laboratory. Such a modulation could be produced by interactions from Weakly Interacting Massive Particles (WIMPs) with masses ~10 GeV/c^2. We find no evidence for annual modulation in the event rate of veto-anticoincident single-detector interactions consistent with nuclear recoils, and constrain the magnitude of any modulation to <0.06 event [keVnr kg day]^-1 in the 5-11.9 keVnr energy range at the 99% confidence level. These results disfavor an explanation for the reported modulation in the 1.2-3.2 keVee energy range in CoGeNT in terms of nuclear recoils resulting from elastic scattering of WIMPs at >98% confidence. For events consistent with electron recoils, no significant modulation is observed for either single- or multiple-detector interactions in the 3.0-7.4 keVee range.

Kaiser redshift-space distortion formula describes well the clustering of galaxies in redshift surveys on small scales, but there are numerous additional terms that arise on large scales. Some of these terms can be described using Newtonian dynamics and have been discussed in the literature, while the others require proper general relativistic description that was only recently developed. Accounting for these terms in galaxy clustering is the first step toward tests of general relativity on horizon scales. The effects can be classified as two terms that represent the velocity and the gravitational potential contributions. Their amplitude is determined by effects such as the volume and luminosity distance fluctuation effects and the time evolution of galaxy number density and Hubble parameter. We compare the Newtonian approximation often used in the redshift-space distortion literature to the fully general relativistic equation, and show that Newtonian approximation accounts for most of the terms contributing to velocity effect. We perform a Fisher matrix analysis of detectability of these terms and show that in a single tracer survey they are completely undetectable. To detect these terms one must resort to the recently developed methods to reduce sampling variance and shot noise. We show that in an all-sky galaxy redshift survey at low redshift the velocity term can be measured at a few sigma if one can utilize halos of mass M>10^12 Msun (this can increase to 10-sigma or more in some more optimistic scenarios), while the gravitational potential term itself can only be marginally detected. We also demonstrate that the general relativistic effect is not degenerate with the primordial non-Gaussian signature in galaxy bias, and the ability to detect primordial non-Gaussianity is little compromised.

The CDMS and EDELWEISS collaborations have combined the results of their direct searches for dark matter using cryogenic germanium detectors. The total data set represents 614 kg.d equivalent exposure. A straightforward method of combination was chosen for its simplicity before data were exchanged between experiments. The results are interpreted in terms of limits on spin-independent WIMP-nucleon cross-section. For a WIMP mass of 90 GeV/c^2, where this analysis is most sensitive, a cross-section of 3.3 x 10^{-44} cm^2 is excluded at 90% CL. At higher WIMP masses, the combination improves the individual limits, by a factor 1.6 above 700 GeV/c^2. Alternative methods of combining the data provide stronger constraints for some ranges of WIMP masses and weaker constraints for others.

Inhomogeneous reionization gives rise to angular fluctuations in the Cosmic Microwave Background (CMB) optical depth tau(n) to the last scattering surface, correlating different spherical harmonic modes and imprinting characteristic non-Gaussianity on CMB maps. Recently the minimum variance quadratic estimator has been derived using this mode-coupling signal, and found that the optical depth fluctuations could be detected with (S/N)^2 ~ 100 in futuristic experiments like CMBPol. We first demonstrate that the non-Gaussian signal from gravitational lensing of CMB is the dominant source of contamination for reconstructing inhomogeneous reionization signals, even with 98% of its contribution removed by delensing. We then construct unbiased estimators that simultaneously reconstruct inhomogeneous reionization signals tau(n) and gravitational lensing potential phi(n). We apply our new unbiased estimators to future CMB experiment to assess the detectability of inhomogeneous reionization signals. With more physically motivated simulations of inhomogenuous reionizations that predict an order of magnitude smaller optical depth power spectrum than previous studies, we show that a CMBPol-like experiment could achieve a marginal detection of inhomogeneous reionization,(S/N)^2 ~ O(1) with this quadratic estimator to ~O(10) with the analogous maximum likelihood estimator.

Results are presented from a reanalysis of the entire five-tower data set acquired with the Cryogenic Dark Matter Search (CDMS II) experiment at the Soudan Underground Laboratory, with an exposure of 969 kg-days. The analysis window was extended to a recoil energy of 150 keV, and an improved surface-event background-rejection cut was defined to increase the sensitivity of the experiment to the inelastic dark matter (iDM) model. Three dark matter candidates were found between 25 keV and 150 keV. The probability to observe three or more background events in this energy range is 11%. Because of the occurrence of these events the constraints on the iDM parameter space are slightly less stringent than those from our previous analysis, which used an energy window of 10-100 keV.

We investigate the effect of supersonic relative velocities between baryons and dark matter, recently shown to arise generically at high redshift, on baryonic acoustic oscillation (BAO) measurements at low redshift. The amplitude of the relative velocity effect at low redshift is model-dependent, but can be parameterized by using an unknown bias. We find that if unaccounted, the relative velocity effect can shift the BAO peak position and bias estimates of the dark energy equation-of-state due to its non-smooth, out-of-phase oscillation structure around the BAO scale. Fortunately, the relative velocity effect can be easily modeled in constraining cosmological parameters without substantially inflating the error budget. We also demonstrate that the presence of the relative velocity effect gives rise to a unique signature in the galaxy bispectrum, which can be utilized to isolate this effect. Future dark energy surveys can accurately measure the relative velocity effect and subtract it from the power spectrum analysis to constrain dark energy models with high precision.

We report results from a reanalysis of data from the Cryogenic Dark Matter Search (CDMS II) experiment at the Soudan Underground Laboratory. Data taken between October 2006 and September 2008 using eight germanium detectors are reanalyzed with a lowered, 2 keV recoil-energy threshold, to give increased sensitivity to interactions from Weakly Interacting Massive Particles (WIMPs) with masses below ~10 GeV/c^2. This analysis provides stronger constraints than previous CDMS II results for WIMP masses below 9 GeV/c^2 and excludes parameter space associated with possible low-mass WIMP signals from the DAMA/LIBRA and CoGeNT experiments.

Data taken during the final shallow-site run of the first tower of the Cryogenic Dark Matter Search (CDMS II) detectors have been reanalyzed with improved sensitivity to small energy depositions. Four ~224 g germanium and two ~105 g silicon detectors were operated at the Stanford Underground Facility (SUF) between December 2001 and June 2002, yielding 118 live days of raw exposure. Three of the germanium and both silicon detectors were analyzed with a new low-threshold technique, making it possible to lower the germanium and silicon analysis thresholds down to the actual trigger thresholds of ~1 keV and ~2 keV, respectively. Limits on the spin-independent cross section for weakly interacting massive particles (WIMPs) to elastically scatter from nuclei based on these data exclude interesting parameter space for WIMPs with masses below 9 GeV/c^2. Under standard halo assumptions, these data partially exclude parameter space favored by interpretations of the DAMA/LIBRA and CoGeNT experiments' data as WIMP signals, and exclude new parameter space for WIMP masses between 3 GeV/c^2 and 4 GeV/c^2.

We extend the general relativistic description of galaxy clustering developed in Yoo, Fitzpatrick, and Zaldarriaga (2009). For the first time we provide a fully general relativistic description of the observed matter power spectrum and the observed galaxy power spectrum with the linear bias ansatz. It is significantly different from the standard Newtonian description on large scales and especially its measurements on large scales can be misinterpreted as the detection of the primordial non-Gaussianity even in the absence thereof. The key difference in the observed galaxy power spectrum arises from the real-space matter fluctuation defined as the matter fluctuation at the hypersurface of the observed redshift. As opposed to the standard description, the shape of the observed galaxy power spectrum evolves in redshift, providing additional cosmological information. While the systematic errors in the standard Newtonian description are negligible in the current galaxy surveys at low redshift, correct general relativistic description is essential for understanding the galaxy power spectrum measurements on large scales in future surveys with redshift depth z>3. We discuss ways to improve the detection significance in the current galaxy surveys and comment on applications of our general relativistic formalism in future surveys.

Measurements of the baryonic acoustic oscillation (BAO) peak in the redshift-space correlation function yield the angular diameter distance D_A(z) and the Hubble parameter H(z) as a function of redshift, constraining the properties of dark energy and space curvature. We discuss the perturbations introduced in the galaxy correlation function by gravitational lensing through the effect of magnification bias and its cross-correlation with the galaxy density. At the BAO scale, gravitational lensing adds a small and slowly varying component to the galaxy correlation function and does not change its shape significantly, through which the BAO peak is measured. The relative shift in the position of the BAO peak caused by gravitational lensing in the angle-averaged correlation function is 10^-4 at z=1, rising to 10^-3 at z=2.5. Lensing effects are stronger near the line-of-sight, however the relative peak shift increases only to 10^-3.3 and 10^-2.4 at z=1 and z=2.5, when the galaxy correlation is averaged within 5 degrees of the line-of-sight (containing only 0.4% of the galaxy pairs in a survey). Furthermore, the lensing contribution can be measured separately and subtracted from the observed correlation at the BAO scale.

We extend our maximum likelihood method for reconstructing the cluster-mass cross-correlation from cosmic microwave background (CMB) temperature anisotropies and develop new estimators that utilize six different quadratic combinations of CMB temperature and polarization fields. Our maximum likelihood estimators are constructed with delensed CMB temperature and polarization fields by using an assumed model of the convergence field and they can be iteratively applied to a set of clusters, approaching to the optimal condition for the lensing reconstruction as the assumed initial model is refined. Using smoothed particle hydrodynamics simulations, we create a catalog of realistic clusters obtainable from the current Sunyaev-Zel'dovich (SZ) surveys, and we demonstrate the ability of the maximum likelihood estimators to reconstruct the cluster-mass cross-correlation from the massive clusters. The iTT temperature estimator provides a signal-to-noise ratio of a factor 3 larger than the iEB polarization estimator, unless the detector noise for measuring polarization anisotropies is controlled under 3 microK.

We present a general relativistic description of galaxy clustering in a FLRW universe. The observed redshift and position of galaxies are affected by the matter fluctuations and the gravity waves between the source galaxies and the observer, and the volume element constructed by using the observables differs from the physical volume occupied by the observed galaxies. Therefore, the observed galaxy fluctuation field contains additional contributions arising from the distortion in observable quantities and these include tensor contributions as well as numerous scalar contributions. We generalize the linear bias approximation to relate the observed galaxy fluctuation field to the underlying matter distribution in a gauge-invariant way. Our full formalism is essential for the consistency of theoretical predictions. As our first application, we compute the angular auto correlation of large-scale structure and its cross correlation with CMB temperature anisotropies. We comment on the possibility of detecting primordial gravity waves using galaxy clustering and discuss further applications of our formalism.

We report on the analysis of the low-energy electron-recoil spectrum from the CDMS II experiment using data with an exposure of 443.2 kg-days. The analysis provides details on the observed counting rate and possible background sources in the energy range of 2 - 8.5 keV. We find no significant excess in the counting rate above background, and compare this observation to the recent DAMA results. In the framework of a conversion of a dark matter particle into electromagnetic energy, our 90% confidence level upper limit of 0.246 events/kg/day at 3.15 keV is lower than the total rate above background observed by DAMA by 8.9$\sigma$. In absence of any specific particle physics model to provide the scaling in cross section between NaI and Ge, we assume a Z^2 scaling. With this assumption the observed rate in DAMA differs from the upper limit in CDMS by 6.8$\sigma$. Under the conservative assumption that the modulation amplitude is 6% of the total rate we obtain upper limits on the modulation amplitude a factor of ~2 less than observed by DAMA, constraining some possible interpretations of this modulation.

We revisit the well known discrepancy between the observed number of Milky Way (MW) dwarf satellite companions and the predicted population of cold dark matter (CDM) sub-halos, in light of the dozen new low luminosity satellites found in SDSS imaging data and our recent calibration of the SDSS satellite detection efficiency, which implies a total population far larger than these dozen discoveries. We combine a dynamical model for the CDM sub-halo population with simple, physically motivated prescriptions for assigning stellar content to each sub-halo, then apply observational selection effects and compare to the current observational census. As expected, models in which the stellar mass is a constant fraction F(Omega_b/Omega_m) of the sub-halo mass M_sat at the time it becomes a satellite fail for any choice of F. However, previously advocated models that invoke suppression of gas accretion after reionization in halos with circular velocity v_c <~ 35 km/s can reproduce the observed satellite counts for -15 < M_V < 0, with F ~ 10^{-3}. Successful models also require strong suppression of star formation BEFORE reionization in halos with v_c <~ 10 km/s; models without pre-reionization suppression predict far too many satellites with -5 < M_V < 0. Our models also reproduce the observed stellar velocity dispersions ~ 5-10 km/s of the SDSS dwarfs given the observed sizes of their stellar distributions, and model satellites have M(<300 pc) ~ 10^7 M_sun as observed even though their present day total halo masses span more than two orders of magnitude. Our modeling shows that natural physical mechanisms acting within the CDM framework can quantitatively explain the properties of the MW satellite population as it is presently known, thus providing a convincing solution to the `missing satellite' problem.

We report the first results from the GammeV search for chameleon particles, which may be created via photon-photon interactions within a strong magnetic field. Chameleons are hypothesized scalar fields that could explain the dark energy problem. We implement a novel technique to create and trap the reflective particles within a jar and to detect them later via their afterglow as they slowly convert back into photons. These measurements provide the first experimental constraints on the couplings of chameleons to photons.

We present a coherent theoretical framework for computing gravitational lensing effects and redshift-space distortions in an inhomogeneous universe and investigate their impacts on galaxy two-point statistics. Adopting the linearized FRW metric, we derive the gravitational lensing and the generalized Sachs-Wolfe effects that include the weak lensing distortion, magnification, and time delay effects, and the redshift-space distortion, Sachs-Wolfe, and integrated Sachs-Wolfe effects, respectively. Based on this framework, we first compute their effects on observed source fluctuations, separating them as two physically distinct origins: the volume effect that involves the change of volume and is always present in galaxy two-point statistics, and the source effect that depends on the intrinsic properties of source populations. Then we identify several terms that are ignored in the standard method, and we compute the observed galaxy two-point statistics, an ensemble average of all the combinations of the intrinsic source fluctuations and the additional contributions from the gravitational lensing and the generalized Sachs-Wolfe effects. This unified treatment of galaxy two-point statistics clarifies the relation of the gravitational lensing and the generalized Sachs-Wolfe effects to the metric perturbations and the underlying matter fluctuations. For near future dark energy surveys, we compute additional contributions to the observed galaxy two-point statistics and analyze their impact on the anisotropic structure. Thorough theoretical modeling of galaxy two-point statistics would be not only necessary to analyze precision measurements from upcoming dark energy surveys, but also provide further discriminatory power in understanding the underlying physical mechanisms.

We develop a new method for reconstructing cluster mass profiles and large-scale structure from the cosmic microwave background (CMB). By analyzing the likelihood of CMB lensing, we analytically prove that standard quadratic estimators for CMB lensing are unbiased and achieve the optimal condition only in the limit of no lensing; they become progressively biased and sub-optimal, when the lensing effect is large, especially for clusters that can be found by ongoing Sunyaev-Zel'dovich surveys. Adopting an alternative approach to the CMB likelihood, we construct a new maximum likelihood estimator that utilizes delensed CMB temperature fields based on an assumed model. We analytically show that this estimator asymptotically approaches the optimal condition as our assumed model is refined, and we numerically show that our estimator quickly converges to the true model as we iteratively apply it to CMB maps. For realistic CMB experiments, we demonstrate the applicability of the maximum likelihood estimator with tests against numerical simulations in the presence of CMB secondary contaminants. With significant improvement on the signal-to-noise ratio, our new maximum likelihood estimator can be used to measure the cluster-mass cross-correlation functions at different redshifts, probing the evolution of dark energy.

The shape of the primordial matter power spectrum Plin(k) encodes critical information on cosmological parameters. At large scales, the observable galaxy power spectrum Pobs(k) is expected to follow the shape of Plin(k), but on smaller scales the effects of nonlinearity and galaxy bias make the ratio Pobs(k)/Plin(k) scale-dependent. We develop a method that can extend the dynamic range of the Plin(k) recovery by incorporating constraints on the galaxy halo occupation distribution (HOD) from the projected galaxy correlation function wp. We devise an analytic model to calculate Pobs(k) in real-space and redshift-space. Once HOD parameters are determined by matching wp for a given cosmological model, galaxy bias is completely specified, and our analytic model predicts both the shape and normalization of Pobs(k). Applying our method to SDSS main galaxy samples, we find that the real-space Pobs(k) follows the shape of the nonlinear matter power spectrum at the 1-2% level up to k=0.2 h/Mpc. When we apply our method to SDSS LRG samples, the linear bias approximation is accurate to 5% at k<0.08 h/Mpc, but the scale-dependence of LRG bias prevents the use of linear theory at k>0.08 h/Mpc. Our HOD model prediction is in good agreement with the recent SDSS LRG Pobs(k) measurements at all measured scales (k<0.2 h/Mpc), naturally explaining the shape of Pobs(k). The "Q-model" prescription is a poor description of galaxy bias for the LRG samples, and it can lead to biased cosmological parameter estimates when measurements at k>0.1 h/Mpc are included in the analysis. We quantify the potential bias and constraints on cosmological parameters that arise from applying linear theory and Q-model fitting, and we demonstrate the utility of HOD modeling of future high precision measurements of Pobs(k) on quasi-linear scales.

We explore the phenomenology of Kaluza-Klein (KK) dark matter in very general models with universal extra dimensions (UEDs), emphasizing the complementarity between high-energy colliders and dark matter direct detection experiments. In models with relatively small mass splittings between the dark matter candidate and the rest of the (colored) spectrum, the collider sensitivity is diminished, but direct detection rates are enhanced. UEDs provide a natural framework for such mass degeneracies. We consider both 5-dimensional and 6-dimensional non-minimal UED models, and discuss the detection prospects for various KK dark matter candidates: the KK photon $\gamma_1$, the KK $Z$-boson $Z_1$, the KK Higgs boson $H_1$ and the spinless KK photon $\gamma_H$. We combine collider limits such as electroweak precision data and expected LHC reach, with cosmological constraints from WMAP, and the sensitivity of current or planned direct detection experiments. Allowing for general mass splittings, we show that neither colliders, nor direct detection experiments by themselves can explore all of the relevant KK dark matter parameter space. Nevertheless, they probe different parameter space regions, and the combination of the two types of constraints can be quite powerful. For example, in the case of $\gamma_1$ in 5D UEDs the relevant parameter space will be almost completely covered by the combined LHC and direct detection sensitivities expected in the near future.

We report the first results of the GammeV experiment, a search for milli-eV mass particles with axion-like couplings to two photons. The search is performed using a "light shining through a wall" technique where incident photons oscillate into new weakly interacting particles that are able to pass through the wall and subsequently regenerate back into detectable photons. The oscillation baseline of the apparatus is variable, thus allowing probes of different values of particle mass. We find no excess of events above background and are able to constrain the two-photon couplings of possible new scalar (pseudoscalar) particles to be less than 3.1x10^{-7} GeV^{-1} (3.5x10^{-7} GeV^{-1}) in the limit of massless particles.

Recent observations support the idea that nuclear black holes grew by gas accretion while shining as luminous quasars at high redshift, and they establish a relation of the black hole mass with the host galaxy's spheroidal stellar system. We develop an analytic model to calculate the expected impact of mergers on the masses of black holes in massive clusters of galaxies. We use the extended Press-Schechter formalism to generate Monte Carlo merger histories of halos with a mass 10^{15} h^{-1} Msun. We assume that the black hole mass function at z=2 is similar to that inferred from observations at z=0 (since quasar activity declines markedly at z<2), and we assign black holes to the progenitor halos assuming a monotonic relation between halo mass and black hole mass. We follow the dynamical evolution of subhalos within larger halos, allowing for tidal stripping, the loss of orbital energy by dynamical friction, and random orbital perturbations in gravitational encounters with subhalos, and we assume that mergers of subhalos are followed by mergers of their central black holes. Our analytic model reproduces numerical estimates of the subhalo mass function. We find that the most massive black holes in massive clusters typically grow by a factor ~ 2 by mergers after gas accretion has stopped. In our ten realizations of 10^{15} h^{-1} Msun clusters, the highest initial (z=2) black hole masses are 5-7 x 10^9 Msun, but four of the clusters contain black holes in the range 1-1.5 x 10^{10} Msun at z=0. Satellite galaxies may host black holes whose mass is comparable to, or even greater than, that of the central galaxy. Thus, black hole mergers can significantly extend the very high end of the black hole mass function.

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 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.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

Galaxy-galaxy lensing measures the mean excess surface density DS(r) around a sample of lensing galaxies. We develop a method for combining DS(r) with the galaxy correlation function xi_gg(r) to constrain Omega_m and sigma_8, going beyond the linear bias model to reach the level of accuracy demanded by current and future measurements. We adopt the halo occupation distribution (HOD) framework, and we test its applicability to this problem by examining the effects of replacing satellite galaxies in the halos of an SPH simulation with randomly selected dark matter particles from the same halos. The difference between dark matter and satellite galaxy radial profiles has a ~10% effect on DS(r) at r<1 Mpc/h. However, if radial profiles are matched, the remaining impact of individual subhalos around satellite galaxies and environmental dependence of the HOD at fixed halo mass is <5% in DS(r) for 0.1<r<15 Mpc/h. We develop an analytic approximation for DS(r) that incorporates halo exclusion and scale-dependent halo bias, and we demonstrate its accuracy with tests against a suite of populated N-body simulations. We use the analytic model to investigate the dependence of DS(r) and the galaxy-matter correlation function xi_gm(r) on Omega_m and sigma_8, once HOD parameters for a given cosmological model are pinned down by matching xi_gg(r). The linear bias prediction is accurate for r>2 Mpc/h, but it fails at the 30-50% level on smaller scales. The scaling of DS(r) ~ Omega_m^a(r) sigma_8^b(r) approaches the linear bias expectation a=b=1 at r>10 Mpc/h, but a(r) and b(r) vary from 0.8 to 1.6 at smaller r. We calculate a fiducial DS(r) and scaling indices a(r) and b(r) for two SDSS galaxy samples; galaxy-galaxy lensing measurements for these samples can be combined with our predictions to constrain Omega_m and sigma_8.

We explore the halo structure of four gravitational lenses with well-observed, thin Einstein rings. We find that the gravitational potentials are well described by ellipsoidal density distributions in the sense that the best-fit nonellipsoidal models have parameters consistent with their ellipsoidal counterparts. We find upper limits on the standard parameters for the deviation from an ellipse of |a_3/a_0|<0.023, 0.019, 0.037, and 0.035, and |a_4/a_0|<0.034, 0.041, 0.051, and 0.064 for SDSS J0924+0219, HE0435-1223, B1938+666, and PG1115+080, respectively. We find that the lens galaxies are at the centers of their dark matter halos, and obtain upper limits for the offset of each center of mass from the center of light of |Delta x|<0.004, 0.005, 0.009, and 0.005 arcsec, corresponding to 22, 29, 70, and 23 pc. These limits also exclude the possibility of any significant lopsidednessof the dark matter halos and set an upper limit of f_sat<sqrt(N)% on the mass fraction of massive substructures inside the Einstein ring if they are divided over N satellites. We also explore the properties of galaxies as substructures in groups for the lens PG1115+080, finding evidence for dark matter halos associated with the galaxies but no evidence for a clear distinction between satellite and central galaxies.

Presently the CDMS-II collaboration's Weakly Interacting Massive Particle (WIMP) search at the Soudan Underground Laboratory sets the most stringent exclusion limits of any WIMP cold dark matter direct-detection experiment. To extend our reach further, to WIMP-nucleon cross sections in the range $10^{-46} - 10^{-44}$cm$^2$, we propose SuperCDMS, which would take advantage of a very deep site. One promising site is the recently approved SNOLab facility in Canada. In this paper we will present our overall plan, identify primary issues, and set the goals that need to be met prior to embarking upon each phase of SuperCDMS.

We use the structure of the Einstein ring image of the quasar host galaxy in the four-image quasar lens PG1115+080 to determine the angular structure of the gravitational potential of the lens galaxy. We find that it is well described as an ellipsoid and that the best fit non-ellipsoidal models are consistent with the ellipsoidal model. We find upper limits on the standard parameters for the m=3 and m=4 deviations from an ellipse of <0.035 and <0.064, respectively. We also find that the position of the center of mass is consistent with the center of light, with an upper limit of 0.005 arcsec on the offset between them. Neither the ellipsoidal nor the non-ellipsoidal models can reproduce the observed image flux ratios while simultaneously maintaining a reasonable fit to the Einstein ring, so the anomalous flux ratio of the A_1 and A_2 quasar images must be due to substructure in the gravitational potential such as compact satellite galaxies or stellar microlenses rather than odd angular structure in the lens galaxy.

The recent discovery of luminous quasars up to a redshift z=6.43 has renewed interest in the formation of black holes massive enough to power quasars. If black holes grow by Eddington-limited gas accretion with a radiative efficiency of at least 10%, the time required to grow from a stellar black hole to ~10^9 msun is ~10^9 years, close to the age of the universe at z=6. Black hole mergers may accelerate the rate of mass growth, but can also completely eject black holes from halo centers owing to the gravitational wave recoil effect. Recently, Haiman concluded that black hole ejections likely do not allow black holes to grow to ~10^9 msun by z=6.43. We reexamine this problem and show that, by using a different halo escape velocity, accounting for the dependence of the recoil velocity on the black hole binary mass ratio and spins, and allowing seed black holes to form in all halos down to virial temperatures of 2000 K, black hole masses may reach ~10^9 msun as early as z=9 starting from stellar seeds, without super-Eddington accretion. In this particular case, we find that these massive black holes form from the merger of ~10^4 stellar black holes formed in low-mass halos at z~20, which must all grow close to the maximum Eddington rate over most of the time available from their birth to z~6. The alternative is that black holes can grow more rapidly by super-Eddington accretion.

We develop a new method of modeling microlensing events based on a Monte Carlo simulation that incorporates both a Galactic model and the constraints imposed by the observed characteristics of the event. The method provides an unbiased way to analyze the event especially when parameters are poorly constrained by the observed lightcurve. We apply this method to search for planetary companions of the lens in OGLE-2003-BLG-423, whose maximum magnification A_max=256+-43 (or A_max=400+-115 from the lightcurve data alone) is the highest among single-lens events ever recorded. The method permits us, for the first time, to place constraints directly in the planet-mass/projected-physical-separation plane rather than in the mass-ratio/Einstein-radius plane as was done previously. For example, Jupiter-mass companions of main-sequence stars at 2.5 AU are excluded with 80% efficiency.

We analyze OGLE-2003-BLG-262, a relatively short, t_E=12.5+-0.1day, microlensing event generated by a point-mass lens transiting the face of a K giant source in the Galactic bulge. We use the resulting finite-source effects to measure the angular Einstein radius, theta_E=195+-17muas, and so constrain the lens mass to the full-width half-maximum interval 0.08 < M/M_sun < 0.54. The lens-source relative proper motion is mu_rel = 27+-2 km/s/kpc. Both values are typical of what is expected for lenses detected toward the bulge. Despite the short duration of the event, we detect marginal evidence for a "parallax asymmetry", but argue that this is more likely to be induced by acceleration of the source, a binary lens, or possibly by statistical fluctuations. Although OGLE-2003-BLG-262 is only the second published event to date in which the lens transits the source, such events will become more common with the new OGLE-III survey in place. We therefore give a detailed account of the analysis of this event to facilitate the study of future events of this type.

We simulate the evolution of halo wide binaries in the presence of MAssive Compact Halo Objects (MACHOs) and compare our results to the sample of wide binaries of Chaname & Gould (2003). The observed distribution is well fit by a single power law for angular separation, 3.5" < theta < 900", whereas the simulated distributions show a break in the power law whose location depends on the MACHO mass and density. This allows us to place upper limits on MACHO density as a function of their assumed mass. At the 95% confidence level, we exclude MACHOs with mass M > 43Msun at the standard local halo density rho_H. This all but removes the last permitted window for a full MACHO halo for masses M > 10^{-7.5}Msun.

We derive constraints on the time variation of the Higgs vacuum expectation value $<\phi>$ through the effects on Big Bang nucleosynthesis (BBN) and the cosmic microwave background (CMB). In the former case, we include the (previously-neglected) effect of the change in the deuteron binding energy, which alters both the $^4$He and deuterium abundances significantly. We find that the current BBN limits on the relative change in \higgs are $-(0.6 - 0.7) \times 10^{-2} < Delta<\phi>/<\phi> < (1.5 - 2.0) \times 10^{-2}$, where the exact limits depend on the model we choose for the dependence of the deuteron binding energy on \higgs.The limits from the current CMB data are much weaker.

In many symbiotic stars there appear broad wings around H alpha, of which the formation mechanism proposed thus far includes a fast outflow, the inner accretion disc motion, electron scattering and Raman scattering of Ly beta. We adopt a Monte Carlo technique to simulate the Raman scattering of UV photons that are converted into optical photons around H alpha forming broad wings, and compute its polarisation. Noting that many symbiotic stars exhibit a bipolar nebular morphology and polarisation flip in the red wing part of the Raman scattered O {\tiny{VI}} features, we assume that the neutral scattering region is composed of the two components. The first component is a static cylindrical shell with finite thickness and the other component is a finite slab that is moving away with velocity v_p = 100 km s^{-1} along the symmetry axis of the first component. The cylindrical shell component yields polarisation in the direction parallel to the cylinder axis. The polarisation near the line-centre is weaker than in the far wing regions because of the large Rayleigh scattering numbers due to the large scattering cross sections near the line centre. The receding polar scattering component produces strong polarisation in the direction perpendicular to the cylinder axis. When the both scattering components coexist, the polarisation is characterised by weak parallel polarisation near the line-centre and strong perpendicular polarisation in the red part. We discuss the observational implications of our computation.

We compute the profiles of resonance doublet lines formed in bipolar winds with velocity greater than the doublet separation in symbiotic stars. Particular attention has been paid on the doublet line ratio, where the conversion of the short wavelength component arising from the $2S_{1/2}-2P_{3/2}$ transition into the long wavelengh component for the transition $2S_{1/2}-2P_{1/2}$ plays an essential role in determining the flux line ratio of the doublet. We adopted a Monte Carlo technique and the Sobolev approximation. Our bipolar winds take the form of a cone and are characterised by the terminal wind velocity, mass loss rate and the opening angle of the cone. When the observer is in the polar direction and the Sobolev optical depth is moderate $\tau_{Sob}\simeq 1$, then we mainly obtain profiles with inverted flux line ratios, where the short wavelength component is weaker than the long wavelength component. When the observer is in the equatorial direction, we find that the profiles are characterised by two broad components, where the long wavelength component is the broader and stronger of the two. We conclude that the profiles obtained in our model provide qualitative understanding of broad and inverted intensity ratio of the doublets in symbiotic stars and that the bipolarity of the stellar winds can be inferred from the broader and stronger long wavelength components in the doublet.