The Kullback-Leibler Divergence as an Estimator of the Statistical Properties of CMB Maps

The identification of unsubtracted foreground residuals in the cosmic microwave background maps on large scales is of crucial importance for the analysis of polarization signals. These residuals add a non-Gaussian contribution to the data. We propose the Kullback-Leibler (KL) divergence as an effective, non-parametric test on the one-point probability distribution function of the data. With motivation in information theory, the KL divergence takes into account the entire range of the distribution and is highly non-local. We demonstrate its use by analyzing the large scales of the Planck 2013 SMICA temperature fluctuation map and find it consistent with the expected distribution at a level of 6%. Comparing the results to those obtained using the more popular Kolmogorov-Smirnov test, we find the two methods to be in general agreement.Comment: 12 pages, 5 figures, minor change, as published in JCA

On the Derivation of Vector Radiative Transfer Equation for Polarized Radiative Transport in Graded Index Media

Light transport in graded index media follows a curved trajectory determined by the Fermat's principle. Besides the effect of variation of the refractive index on the transport of radiative intensity, the curved ray trajectory will induce geometrical effects on the transport of polarization ellipse. This paper presents a complete derivation of vector radiative transfer equation for polarized radiation transport in absorption, emission and scattering graded index media. The derivation is based on the analysis of the conserved quantities for polarized light transport along curved trajectory and a novel approach. The obtained transfer equation can be considered as a generalization of the classic vector radiative transfer equation that is only valid for uniform refractive index media. Several variant forms of the transport equation are also presented, which include the form for Stokes parameters defined with a fixed reference and the Eulerian forms in the ray coordinate and in several common orthogonal coordinate systems.Comment: This paper has been submitted to JQSR

Quantifying jet transport properties via large pTp_T hadron production

Nuclear modification factor $R_{AA}$ for large $p_T$ single hadron is studied in a next-to-leading order (NLO) perturbative QCD (pQCD) parton model with medium-modified fragmentation functions (mFFs) due to jet quenching in high-energy heavy-ion collisions. The energy loss of the hard partons in the QGP is incorporated in the mFFs which utilize two most important parameters to characterize the transport properties of the hard parton jets: the jet transport parameter $\hat q_{0}$ and the mean free path $\lambda_{0}$, both at the initial time $\tau_0$. A phenomenological study of the experimental data for $R_{AA}(p_{T})$ is performed to constrain the two parameters with simultaneous $\chi^2/{\rm d.o.f}$ fits to RHIC as well as LHC data. We obtain for energetic quarks $\hat q_{0}\approx 1.1 \pm 0.2$ GeV$^2$/fm and $\lambda_{0}\approx 0.4 \pm 0.03$ fm in central $Au+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV, while $\hat q_{0}\approx 1.7 \pm 0.3$ GeV$^2$/fm, and $\lambda_{0}\approx 0.5 \pm 0.05$ fm in central $Pb+Pb$ collisions at $\sqrt{s_{NN}}=2.76$ TeV. Numerical analysis shows that the best fit favors a multiple scattering picture for the energetic jets propagating through the bulk medium, with a moderate averaged number of gluon emissions. Based on the best constraints for $\lambda_{0}$ and $\tau_0$, the estimated value for the mean-squared transverse momentum broadening is moderate which implies that the hard jets go through the medium with small reflection.Comment: 8 pages, 6 figures, revised versio