We develop a Feynman diagram approach to calculating correlations of the Cosmic Microwave Background (CMB) in the presence of distortions. As one application, we focus on CMB distortions due to gravitational lensing by Large Scale Structure (LSS). We study the Hu-Okamoto quadratic estimator for extracting lensing from the CMB and derive the noise of the estimator up to ${\mathcal O}(\phi^4)$ in the lensing potential $\phi$. The previously noted large ${\mathcal O}(\phi^4)$ term can be significantly reduced by a reorganization of the $\phi$ expansion. Our approach makes it simple to obtain expressions for quadratic estimators based on any CMB channel. We briefly discuss other applications to cosmology of this diagrammatic approach, such as distortions of the CMB due to patchy reionization, or due to Faraday rotation from primordial axion fields.

We develop a method for calculating the correlation structure of the Cosmic Microwave Background (CMB) using Feynman diagrams, when the CMB has been modified by gravitational lensing, Faraday rotation, patchy reionization, or other distorting effects. This method is used to calculate the bias of the Hu-Okamoto quadratic estimator in reconstructing the lensing power spectrum up to O(\phi^4) in the lensing potential $\phi$. We consider both the diagonal noise TTTT, EBEB, etc. and, for the first time, the off-diagonal noise TTTE, TBEB, etc. The previously noted large O(\phi^4) term in the second order noise is identified to come from a particular class of diagrams. It can be significantly reduced by a reorganization of the $\phi$ expansion. These improved estimators have almost no bias for the off-diagonal case involving only one $B$ component of the CMB, such as EEEB.

We analyze seven year and nine year WMAP temperature maps for signatures of three finite flat topologies M_0=T^3, M_1=T^2 x R^1, and M_2=S^1 x R^2. We use Monte-Carlo simulations with the Feldman-Cousins method to obtain confidence intervals for the size of the topologies considered. We analyze the V, W, and Q frequency bands along with the ILC map and find no significant difference in the results. The 95.5% confidence level lower bound on the size of the topology is 1.5L_0 for M_0, 1.4L_0 for M_1, and 1.1L_0 for M_2, where L_0 is the radius of the last scattering surface. Our results agree very well with the recently released results from the Planck temperature data. We show that the likelihood function is not Gaussian in the size, and therefore simulations are important for obtaining accurate bounds on the size. We then introduce the formalism for including polarization data in the analysis. The improvement that we find from WMAP polarization maps is small because of the high level of instrumental noise, but our forecast for Planck maps shows a much better improvement on the lower bound for L. For the M_0 topology we expect an improvement on the lower bound of L from 1.7L_0 to 1.9L_0 at 95.5% confidence level. Using both polarization and temperature data is important because it tests the hypothesis that deviations in the TT spectrum at small l originate in the primordial perturbation spectrum.

We place limits on semiclassical fluctuations that might be present in the primordial perturbation spectrum. These can arise if some signatures of pre-inflationary features survive the expansion, or could be created by whatever mechanism ends inflation. We study two possible models for such remnant fluctuations, both of which break the isotropy of CMB on large scales. We first consider a semiclassical fluctuation in one Fourier mode of primordial perturbations. The second scenario we analyze is a semiclassical Gaussian bump somewhere in space. These models are tested with the seven-year WMAP data using a Markov Chain Monte Carlo Bayesian analysis, and we place limits on these fluctuations. The upper bound for the amplitude of a fluctuation in a single Fourier mode is a<=10^(-4), while for the Gaussian bump a<=10^(-3).

Constraints are placed on the spin-independent interaction cross section of dark matter with regular matter by refining two methods. First, dark matter--cosmic ray interactions are considered, wherein cosmic ray protons collide with dark matter to contribute to the gamma ray sky. This constraint is developed using the NFW and Moore dark matter density profiles and new data from the Fermi gamma ray space telescope. Second, the Earth capture scenario is considered, wherein particles that are captured self-annihilate at Earth's center, thus adding to its internal heat flow. The constraint presented here is developed based on analysis of the drift time of dark matter particles through Earth, modeled as a core composed of iron and a mantle composed of oxygen with linear density gradients between layers. An analysis of the cosmic ray constraint (which rules out dark matter--regular matter interaction cross sections greater than its value) shows that it overlaps significantly with the Earth drift time constraint (which rules out cross sections smaller than its value), closing the window on strongly interacting dark matter particles up to a mass of about 10^{17} GeV when combined with other exclusions.

We cross correlate the gravitational lensing map extracted from cosmic microwave background measurements by the Wilkinson Microwave Anisotropy Probe (WMAP) with the radio galaxy distribution from the NRAO VLA Sky Survey (NVSS) by using a quadratic estimator technique. We use the full covariance matrix to filter the data, and calculate the cross-power spectra for the lensing-galaxy correlation. We explore the impact of changing the values of cosmological parameters on the lensing reconstruction, and obtain statistical detection significances at $>3\sigma$. The results of all cross correlations pass the curl null test as well as a complementary diagnostic test using the NVSS data in equatorial coordinates. We forecast the potential for Planck and NVSS to constrain the lensing-galaxy cross correlation as well as the galaxy bias. The lensing-galaxy cross-power spectra are found to be Gaussian distributed.

We study the possibility that the universe has compact topologies T^3, T^2 x R^1, or S^1 x R^2 using the seven-year WMAP data. The maximum likelihood 95% confidence intervals for the size L of the compact direction are 1.7 < L/L_0 < 2.1, 1.8 < L/L_0 < 2.0, 1.2 < L/L_0 < 2.1 for the three cases, respectively, where L_0=14.4 Gpc is the distance to the last scattering surface. An infinite universe is compatible with the data at 4.3 sigma. We find using a Bayesian analysis that the most probable universe has topology T^2 x R^1, with L/L_0=1.9.

Dark energy models which alter the relative scaling behavior of dark energy and matter could provide a natural solution to the cosmic coincidence problem - why the densities of dark energy and dark matter are comparable today. A generalized class of dark energy models is introduced which allows non-canonical scaling of the ratio of dark matter and dark energy with the Robertson-Walker scale factor a(t). Upcoming observations, such as a high redshift supernova survey, application of the Alcock-Paczynski test to quasar pairs, and cluster evolution, will strongly constrain the relative scaling of dark matter and dark energy as well as the equation of state of the dark energy. Thus, whether there actually is a coincidence problem, and the extent of cosmic coincidence in the universe's recent past can be answered observationally in the near future. Determining whether today is a special time in the history of the universe will be a SNAP.