CMB likelihood approximation by a Gaussianized Blackwell-Rao estimator

We introduce a new CMB temperature likelihood approximation called the Gaussianized Blackwell-Rao (GBR) estimator. This estimator is derived by transforming the observed marginal power spectrum distributions obtained by the CMB Gibbs sampler into standard univariate Gaussians, and then approximate their joint transformed distribution by a multivariate Gaussian. The method is exact for full-sky coverage and uniform noise, and an excellent approximation for sky cuts and scanning patterns relevant for modern satellite experiments such as WMAP and Planck. A single evaluation of this estimator between l=2 and 200 takes ~0.2 CPU milliseconds, while for comparison, a single pixel space likelihood evaluation between l=2 and 30 for a map with ~2500 pixels requires ~20 seconds. We apply this tool to the 5-year WMAP temperature data, and re-estimate the angular temperature power spectrum, $C_{\ell}$, and likelihood, L(C_l), for l<=200, and derive new cosmological parameters for the standard six-parameter LambdaCDM model. Our spectrum is in excellent agreement with the official WMAP spectrum, but we find slight differences in the derived cosmological parameters. Most importantly, the spectral index of scalar perturbations is n_s=0.973 +/- 0.014, 1.9 sigma away from unity and 0.6 sigma higher than the official WMAP result, n_s = 0.965 +/- 0.014. This suggests that an exact likelihood treatment is required to higher l's than previously believed, reinforcing and extending our conclusions from the 3-year WMAP analysis. In that case, we found that the sub-optimal likelihood approximation adopted between l=12 and 30 by the WMAP team biased n_s low by 0.4 sigma, while here we find that the same approximation between l=30 and 200 introduces a bias of 0.6 sigma in n_s.Comment: 10 pages, 7 figures, submitted to Ap

Frequentist comparison of CMB local extrema statistics in the five-year WMAP data with two anisotropic cosmological models

We present local extrema studies of two models that introduce a preferred direction into the observed cosmic microwave background (CMB) temperature field. In particular, we make a frequentist comparison of the one- and two-point statistics for the dipole modulation and ACW models with data from the five-year Wilkinson Microwave Anisotropy Probe (WMAP). This analysis is motivated by previously revealed anomalies in the WMAP data, and particularly the difference in the statistical nature of the temperature anisotropies when analysed in hemispherical partitions. The analysis of the one-point statistics indicates that the previously determined hemispherical variance difficulties can be apparently overcome by a dipole modulation field, but new inconsistencies arise if the mean and the l-dependence of the statistics are considered. The two-point correlation functions of the local extrema, the temperature pair product and the point-point spatial pair-count, demonstrate that the impact of such a modulation is to over-`asymmetrise' the temperature field on smaller scales than the wave-length of the dipole or quadrupole, and this is disfavored by the observed data.The results from the ACW model predictions, however, are consistent with the standard isotropic hypothesis. The two-point analysis confirms that the impact of this type of violation of isotropy on the temperature extrema statistics is relatively weak. From this work, we conclude that a model with more spatial structure than the dipole modulated or rotational-invariance breaking models are required to fully explain the observed large-scale anomalies in the WMAP data.Comment: 10 pages, 4 figures, 3 tables, accepted for publication in MNRA

Systematic Distortion in Cosmic Microwave Background Maps

To minimize instrumentally induced systematic errors, cosmic microwave background (CMB) anisotropy experiments measure temperature differences across the sky using paires of horn antennas, temperature map is recovered from temperature differences obtained in sky survey through a map-making procedure. To inspect and calibrate residual systematic errors in recovered temperature maps is important as most previous studies of cosmology are based on these maps. By analyzing pixel-ring couping and latitude dependence of CMB temperatures, we find notable systematic deviation from CMB Gaussianity in released Wilkinson Microwave Anisotropy Probe (WMAP) maps. The detected deviation is hard to explain by any process in the early universe and can not be ignored for a precision cosmology study.Comment: accepted for publication in Sci China G-Phy Mech Astro

Slagtesvin på græs - produktionsmæssige muligheder og miljømæssige risici

Slagtesvin på friland er et reelt alternativ til den nuværende mest almindelige indendørsproduktion ud fra et produktions- og kødkvalitetsmæssigt synspunkt. Men det er dog også klart, at produktionen indebærer en betydelig risiko for miljømæssige ulemper i form af et stort N-tab til omgivelserne

Real space tests of the statistical isotropy and Gaussianity of the WMAP CMB data

ABRIDGED: We introduce and analyze a method for testing statistical isotropy and Gaussianity and apply it to the WMAP CMB foreground reduced, temperature maps, and cross-channel difference maps. We divide the sky into regions of varying size and shape and measure the first four moments of the one-point distribution within these regions, and using their simulated spatial distributions we test the statistical isotropy and Gaussianity hypotheses. By randomly varying orientations of these regions, we sample the underlying CMB field in a new manner, that offers a richer exploration of the data content, and avoids possible biasing due to a single choice of sky division. The statistical significance is assessed via comparison with realistic Monte-Carlo simulations. We find the three-year WMAP maps to agree well with the isotropic, Gaussian random field simulations as probed by regions corresponding to the angular scales ranging from 6 deg to 30 deg at 68% confidence level. We report a strong, anomalous (99.8% CL) dipole ``excess'' in the V band of the three-year WMAP data and also in the V band of the WMAP five-year data (99.3% CL). We notice the large scale hemispherical power asymmetry, and find that it is not highly statistically significant in the WMAP three-year data (<~ 97%) at scales l <= 40. The significance is even smaller if multipoles up to l=1024 are considered (~90% CL). We give constraints on the amplitude of the previously-proposed CMB dipole modulation field parameter. We easily detect the residual foregrounds in cross-band difference maps at rms level <~ 7 \mu K (at scales >~ 6 deg) and limit the systematical uncertainties to <~ 1.7 \mu K (at scales >~ 30 deg).Comment: 20 pages, 20 figures; more tests added; updated to match the version to be published in JCA

The Flux Auto- and Cross-Correlation of the Lyman-alpha Forest. II. Modelling Anisotropies with Cosmological Hydrodynamic Simulations

The isotropy of the Lyman-alpha forest in real-space uniquely provides a measurement of cosmic geometry at z > 2. The angular diameter distance for which the correlation function along the line of sight and in the transverse direction agree corresponds to the correct cosmological model. However, the Lyman-alpha forest is observed in redshift-space where distortions due to Hubble expansion, bulk flows, and thermal broadening introduce anisotropy. Similarly, a spectrograph's line spread function affects the autocorrelation and cross-correlation differently. In this the second paper of a series on using the Lyman-alpha forest observed in pairs of QSOs for a new application of the Alcock-Paczynski (AP) test, these anisotropies and related sources of potential systematic error are investigated with cosmological hydrodynamic simulations. Three prescriptions for galactic outflow were compared and found to have only a marginal effect on the Lyman-alpha flux correlation (which changed by at most 7% with use of the currently favored variable-momentum wind model vs. no winds at all). An approximate solution for obtaining the zero-lag cross-correlation corresponding to arbitrary spectral resolution directly from the zero-lag cross-correlation computed at full-resolution (good to within 2% at the scales of interest) is presented. Uncertainty in the observationally determined mean flux decrement of the Lyman-alpha forest was found to be the dominant source of systematic error; however, this is reduced significantly when considering correlation ratios. We describe a simple scheme for implementing our results, while mitigating systematic errors, in the context of a future application of the AP test.Comment: 20 page

Computation of local exchange coefficients in strongly interacting one-dimensional few-body systems: local density approximation and exact results

One-dimensional multi-component Fermi or Bose systems with strong zero-range interactions can be described in terms of local exchange coefficients and mapping the problem into a spin model is thus possible. For arbitrary external confining potentials the local exchanges are given by highly non-trivial geometric factors that depend solely on the geometry of the confinement through the single-particle eigenstates of the external potential. To obtain accurate effective Hamiltonians to describe such systems one needs to be able to compute these geometric factors with high precision which is difficult due to the computational complexity of the high-dimensional integrals involved. An approach using the local density approximation would therefore be a most welcome approximation due to its simplicity. Here we assess the accuracy of the local density approximation by going beyond the simple harmonic oscillator that has been the focus of previous studies and consider some double-wells of current experimental interest. We find that the local density approximation works quite well as long as the potentials resemble harmonic wells but break down for larger barriers. In order to explore the consequences of applying the local density approximation in a concrete setup we consider quantum state transfer in the effective spin models that one obtains. Here we find that even minute deviations in the local exchange coefficients between the exact and the local density approximation can induce large deviations in the fidelity of state transfer for four, five, and six particles.Comment: 12 pages, 7 figures, 1 table, final versio

Power Asymmetries in the Cosmic Microwave Background Temperature and Polarization patterns

We test the asymmetry of the Cosmic Microwave Background anisotropy jointly in temperature and polarization. We study the hemispherical asymmetry, previously found only in the temperature field, with respect to the axis identified by Hansen et al. (2009). To this extent, we make use of the low resolution WMAP 5 year temperature and polarization Nside=16 maps and the optimal angular power spectrum estimator BolPol (Gruppuso et al. 2009). We consider two simple estimators for the power asymmetry and we compare our findings with Monte Carlo simulations which take into account the full noise covariance matrix. We confirm an excess of power in temperature angular power spectrum in the Southern hemisphere at a significant level, between 3 sigma and 4 sigma depending on the exact range of multipoles considered. We do not find significant power asymmetry in the gradient (curl) component EE (BB) of polarized angular spectra. Also cross-correlation power spectra, i.e. TE, TB, EB, show no significant hemispherical asymmetry. We also show that the Cold Spot found by Vielva et al. (2004) in the Southern Galactic hemisphere does not alter the significance of the hemispherical asymmetries on multipoles which can be probed by maps at resolution Nside=16. Although the origin of the hemispherical asymmetry in temperature remains unclear, the study of the polarization patter could add useful information on its explanation. We therefore forecast by Monte Carlo the Planck capabilities in probing polarization asymmetries.Comment: 6 pages, 7 figure