The kinetic Sunyaev Zel'dovich (kSZ) and polarized Sunyaev Zel'dovich (pSZ) effects are temperature and polarization anisotropies induced by the scattering of CMB photons from structure in the post-reionization Universe. In the case of the kSZ effect, small angular scale anisotropies in the optical depth are modulated by the cosmic microwave background (CMB) dipole field, i.e. the CMB dipole observed at each spacetime point, which is sourced by the primordial dipole and especially the local peculiar velocity. In the case of the pSZ effect, similar small-scale anisotropies are modulated by the CMB quadrupole field, which receives contributions from both scalar and tensor modes. Statistical anisotropies in the cross correlations of CMB temperature and polarization with tracers of the inhomogeneous distribution of electrons provide a means of isolating and reconstructing the dipole and quadrupole fields. In this paper, we present a set of unbiased minimum variance quadratic estimators for the reconstruction of the dipole and quadrupole fields, and forecast the ability of future CMB experiments and large scale structure surveys to perform this reconstruction. Consistent with previous work, we find that a high fidelity reconstruction of the dipole and quadrupole fields over a variety of scales is indeed possible, and demonstrate the sensitivity of the pSZ effect to primordial tensor modes. Using a principle component analysis, we estimate how many independent modes could be accessed in such a reconstruction. We also comment on a few first applications of a detection of the dipole and quadrupole fields, including a reconstruction of the primordial contribution to our locally observed CMB dipole, a test of statistical homogeneity on large scales from the first modes of the quadrupole field, and a reconstruction technique for the primordial potential on the largest scales.

Secondary CMB polarization is induced by the late-time scattering of CMB photons by free electrons on our past light cone. This polarized Sunyaev Zel'dovich (pSZ) effect is sensitive to the electrons' locally observed CMB quadrupole, which is sourced primarily by long wavelength inhomogeneities. By combining the remote quadrupoles measured by free electrons throughout the Universe after reionization, the pSZ effect allows us to obtain additional information about large scale modes beyond what can be learned from our own last scattering surface. Here we determine the power of pSZ tomography, in which the pSZ effect is cross-correlated with the density field binned at several redshifts, to provide information about the long wavelength Universe. The signal we explore here is a power asymmetry in the cross-correlation between $E$ or $B$ mode CMB polarization and the density field. We compare this to the cosmic variance limited noise: the random chance to get a power asymmetry in the absence of a large scale quadrupole field. By computing the necessary transfer functions and cross-correlations, we compute the signal-to-noise ratio attainable by idealized next generation CMB experiments and galaxy surveys. We find that a signal-to-noise ratio of $\sim 1-10$ is in principle attainable over a significant range of power multipoles, with the strongest signal coming from the first multipoles in the lowest redshift bins. These results prompt further assessment of realistically measuring the pSZ signal and the potential impact for constraining cosmology on large scales.

Due to cosmic variance we cannot learn any more about large-scale inhomogeneities from the primary cosmic microwave background (CMB) alone. More information on large scales is essential for resolving large angular scale anomalies in the CMB. Here we consider cross correlating the large-scale kinetic Sunyaev Zel'dovich (kSZ) effect and probes of large-scale structure, a technique known as kSZ tomography. The statistically anisotropic component of the cross correlation encodes the CMB dipole as seen by free electrons throughout the observable Universe, providing information about long wavelength inhomogeneities. We compute the large angular scale power asymmetry, constructing the appropriate transfer functions, and estimate the cosmic variance limited signal to noise for a variety of redshift bin configurations. The signal to noise is significant over a large range of power multipoles and numbers of bins. We present a simple mode counting argument indicating that kSZ tomography can be used to estimate more modes than the primary CMB on comparable scales. A basic forecast indicates that a first detection could be made with next-generation CMB experiments and galaxy surveys. This paper motivates a more systematic investigation of how close to the cosmic variance limit it will be possible to get with future observations.

Massive bigravity, a theoretically consistent modification of general relativity with an additional dynamical rank two tensor, successfully describes the observed accelerated expansion of the Universe without a cosmological constant. Recent analyses of perturbations around a cosmological background have revealed power law instabilities in both scalar and tensor perturbations, motivating an analysis of the initial conditions, evolution, and cosmological observables to determine the viability of these theories. In this paper we focus on the tensor sector, and study a primordial stochastic gravitational wave background in massive bigravity. The phenomenology can differ from standard General Relativity due to non-trivial mixing between the two linearized tensor fluctuations in the theory, only one of which couples to matter. We study perturbations about two classes of cosmological solutions in bigravity, computing the tensor contribution to the temperature anisotropies in the Cosmic Microwave Background radiation and the present stochastic gravitational wave background. The result is strongly dependent on the choice of cosmological background and initial conditions. One class of background solution generically displaying tremendous growth in the amplitude of large-wavelength gravitational waves, while the other remains observationally indistinguishable from standard General Relativity for a wide variety of initial conditions. We analyze the initial conditions for tensor modes expected in an inflationary cosmology, finding again that there is a strong dependence on the assumed background. For one choice of background, the semi-classical theory is beyond the perturbative regime. For the other choice, inflation generically yields initial conditions that, when evolved, give rise to a stochastic background observationally indistinguishable from standard General Relativity.

Scalar forces "screened" by the Vainshtein mechanism may hold the key to understanding the cosmological expansion of our universe, while predicting new and exciting features in the interaction between massive bodies. Here we explore the dynamics of the Vainshtein screening mechanism, focusing on the decoupling limit of the DGP braneworld scenario and dRGT massive gravity. We show that there is a vast set of initial conditions whose evolution is well defined and which are driven to the static, screening solutions of these theories. Screening solutions are stable and behave coherently under small fluctuations: they oscillate and eventually settle to an equilibrium configuration, the timescale for the oscillations and damping being dictated by the Vainshtein radius of the screening solutions. At very late times, a power-law decay ensues, in agreement with known analytical results. However, we also conjecture that physically interesting processes such as gravitational collapse of compact stars may not possess a well-posed initial value problem. Finally, we construct solutions with nontrivial multipolar structure describing the screening field of deformed, asymmetric bodies and show that higher multipoles are screened more efficiently than the monopole component.