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.