skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Higher-order gravitational lensing reconstruction using Feynman diagrams

Abstract

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 (φ{sup 4}) in the lensing potential φ. We consider both the diagonal noise TT TT, EB EB, etc. and, for the first time, the off-diagonal noise TT TE, TB EB, etc. The previously noted large  O (φ{sup 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 φ expansion. These improved estimators have almost no bias for the off-diagonal case involving only one B component of the CMB, such as EE EB.

Authors:
; ;  [1];  [2]
  1. Department of Physics, University of California at San Diego, La Jolla, CA 92093 (United States)
  2. Nikhef, Theory Group, Science Park 105, 1098 XG, Amsterdam (Netherlands)
Publication Date:
OSTI Identifier:
22375883
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Cosmology and Astroparticle Physics; Journal Volume: 2014; Journal Issue: 09; Other Information: Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
79 ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; FARADAY EFFECT; FEYNMAN DIAGRAM; GRAVITATIONAL LENSES; NOISE; POTENTIALS; RELICT RADIATION; SPECTRA

Citation Formats

Jenkins, Elizabeth E., Manohar, Aneesh V., Yadav, Amit P.S., and Waalewijn, Wouter J., E-mail: ejenkins@ucsd.edu, E-mail: amanohar@ucsd.edu, E-mail: W.J.Waalewijn@uva.nl, E-mail: ayadav@physics.ucsd.edu. Higher-order gravitational lensing reconstruction using Feynman diagrams. United States: N. p., 2014. Web. doi:10.1088/1475-7516/2014/09/024.
Jenkins, Elizabeth E., Manohar, Aneesh V., Yadav, Amit P.S., & Waalewijn, Wouter J., E-mail: ejenkins@ucsd.edu, E-mail: amanohar@ucsd.edu, E-mail: W.J.Waalewijn@uva.nl, E-mail: ayadav@physics.ucsd.edu. Higher-order gravitational lensing reconstruction using Feynman diagrams. United States. doi:10.1088/1475-7516/2014/09/024.
Jenkins, Elizabeth E., Manohar, Aneesh V., Yadav, Amit P.S., and Waalewijn, Wouter J., E-mail: ejenkins@ucsd.edu, E-mail: amanohar@ucsd.edu, E-mail: W.J.Waalewijn@uva.nl, E-mail: ayadav@physics.ucsd.edu. Mon . "Higher-order gravitational lensing reconstruction using Feynman diagrams". United States. doi:10.1088/1475-7516/2014/09/024.
@article{osti_22375883,
title = {Higher-order gravitational lensing reconstruction using Feynman diagrams},
author = {Jenkins, Elizabeth E. and Manohar, Aneesh V. and Yadav, Amit P.S. and Waalewijn, Wouter J., E-mail: ejenkins@ucsd.edu, E-mail: amanohar@ucsd.edu, E-mail: W.J.Waalewijn@uva.nl, E-mail: ayadav@physics.ucsd.edu},
abstractNote = {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 (φ{sup 4}) in the lensing potential φ. We consider both the diagonal noise TT TT, EB EB, etc. and, for the first time, the off-diagonal noise TT TE, TB EB, etc. The previously noted large  O (φ{sup 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 φ expansion. These improved estimators have almost no bias for the off-diagonal case involving only one B component of the CMB, such as EE EB.},
doi = {10.1088/1475-7516/2014/09/024},
journal = {Journal of Cosmology and Astroparticle Physics},
number = 09,
volume = 2014,
place = {United States},
year = {Mon Sep 01 00:00:00 EDT 2014},
month = {Mon Sep 01 00:00:00 EDT 2014}
}
  • We extend to higher order a recently published method for calculating the deflection angle of light in a general static and spherically symmetric metric. We have tested our method on the metric of Schwarzschild and Reissner-Nordstroem black holes, on the metric of a charged black hole coupled to Born-Infeld electrodynamics and on the metric of Weyl gravity. Since our method is geometrically convergent, as proved in our previous work, our analytical formulas obtained working to fourth order are sufficient to reach errors of few percents even in proximity of the photon sphere.
  • While the ML-EM (maximum-likelihood-expectation maximization) algorithm for reconstruction for emission tomography is unstable due to the ill-posed nature of the problem, Bayesian reconstruction methods overcome this instability by introducing prior information, often in the form of a spatial smoothness regularizer. More elaborate forms of smoothness constraints may be used to extend the role of the prior beyond that of a stabilizer in order to capture actual spatial information about the object. Previously proposed forms of such prior distributions were based on the assumption of a piecewise constant source distribution. Here, the authors propose an extension to a piecewise linear model--themore » weak plate--which is more expressive than the piecewise constant model. The weak plate prior not only preserves edges but also allows for piecewise ramplike regions in the reconstruction. Indeed, for the application in SPECT, such ramplike regions are observed in ground-truth source distributions in the form of primate autoradiographs of rCBF radionuclides. To incorporate the weak plate prior in a MAP approach, the authors model the prior as a Gibbs distribution and use a GEM formulation for the optimization. They compare quantitative performance of the ML-EM algorithm, a GEM algorithm with a prior favoring piecewise constant regions, and a GEM algorithm with the weak plate prior. Pointwise and regional bias and variance of ensemble image reconstructions are used as indications of image quality. The results show that the weak plate and membrane priors exhibit improved bias and variance relative to ML-EM techniques.« less
  • This article is devoted to the numerical discretisation of the hyperbolic two-phase flow model of Baer and Nunziato. A special attention is paid on the discretisation of intercell flux functions in the framework of Finite Volume and Discontinuous Galerkin approaches, where care has to be taken to efficiently approximate the non-conservative products inherent to the model equations. Various upwind approximate Riemann solvers have been tested on a bench of discontinuous test cases. New discretisation schemes are proposed in a Discontinuous Galerkin framework following the criterion of Abgrall and the path-conservative formalism. A stabilisation technique based on artificial viscosity is appliedmore » to the high-order Discontinuous Galerkin method and compared against classical TVD-MUSCL Finite Volume flux reconstruction.« less
  • PET images are affected by the presence of scattered photons. Incorrect scatter-correction may cause artifacts, particularly in 3D PET systems. Current scatter reconstruction methods do not distinguish between single and higher order scattered photons. A dual-scattered reconstruction method (GDS-MLEM) that is independent of the number of Compton scattering interactions and less sensitive to the need for high energy resolution detectors, is proposed. To avoid overcorrecting for scattered coincidences, the attenuation coefficient was calculated by integrating the differential Klein-Nishina cross-section over a restricted energy range, accounting only for scattered photons that were not detected. The optimum image can be selected bymore » choosing an energy threshold which is the upper energy limit for the calculation of the cross-section and the lower limit for scattered photons in the reconstruction. Data was simulated using the GATE platform. 500,000 multiple scattered photon coincidences with perfect energy resolution were reconstructed using various methods. The GDS-MLEM algorithm had the highest confidence (98%) in locating the annihilation position and was capable of reconstructing the two largest hot regions. 100,000 photon coincidences, with a scatter fraction of 40%, were used to test the energy resolution dependence of different algorithms. With a 350–650 keV energy window and the restricted attenuation correction model, the GDS-MLEM algorithm was able to improve contrast recovery and reduce the noise by 7.56%–13.24% and 12.4%–24.03%, respectively. This approach is less sensitive to the energy resolution and shows promise if detector energy resolutions of 12% can be achieved.« less