BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

K. M. Gorski

K. M. Gorski contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-phastro-ph.COphysics.data-anastro-ph.IM

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2011arXiv

Search for non-Gaussianities in the WMAP data with the Scaling Index Method

In the recent years, non-Gaussianity and statistical isotropy of the Cosmic Microwave Background (CMB) was investigated with various statistical measures, first and foremost by means of the measurements of the WMAP satellite. In this Review, we focus on the analyses that were accomplished with a measure of local type, the so-called Scaling Index Method (SIM). The SIM is able to detect structural characteristics of a given data set, and has proven to be highly valuable in CMB analysis. It was used for comparing the data set with simulations as well as surrogates, which are full sky maps generated by randomisation of previously selected features of the original map. During these investigations, strong evidence for non-Gaussianities as well as asymmetries and local features could be detected. In combination with the surrogates approach, the SIM detected the highest significances for non-Gaussianity to date.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Measuring Planck beams with planets

Accurate measurement of the cosmic microwave background (CMB) anisotropy requires precise knowledge of the instrument beam. We explore how well the Planck beams will be determined from observations of planets, developing techniques that are also appropriate for other experiments. We simulate planet observations with a Planck-like scanning strategy, telescope beams, noise, and detector properties. Then we employ both parametric and non-parametric techniques, reconstructing beams directly from the time-ordered data. With a faithful parameterization of the beam shape, we can constrain certain detector properties, such as the time constants of the detectors, to high precision. Alternatively, we decompose the beam using an orthogonal basis. For both techniques, we characterize the errors in the beam reconstruction with Monte Carlo realizations. For a simplified scanning strategy, we study the impact on estimation of the CMB power spectrum. Finally, we explore the consequences for measuring cosmological parameters, focusing on the spectral index of primordial scalar perturbations, n_s. The quality of the power spectrum measurement will be significantly influenced by the optical modeling of the telescope. In our most conservative case, using no information about the optics except the measurement of planets, we find that a single transit of Jupiter across the focal plane will measure the beam window functions to better than 0.3% for the channels at 100-217 GHz that are the most sensitive to the CMB. Constraining the beam with optical modeling can lead to much higher quality reconstruction. Depending on the optical modeling, the beam errors may be a significant contribution to the measurement systematics for n_s.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Performance of XFaster likelihood in real CMB experiments

We assess the strengths and weaknesses of several likelihood formalisms, including the XFaster likelihood. We compare the performance of the XFaster likelihood to that of the Offset Lognormal Bandpower likelihood on simulated data for the Planck satellite. Parameters estimated with these two likelihoods are in good agreement. The advantages of the XFaster likelihood can therefore be realized without compromising performance.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Probing non-Gaussianities on Large Scales in WMAP5 and WMAP7 Data using Surrogates

Probing Gaussianity represents one of the key questions in modern cosmology, because it allows to discriminate between different models of inflation. We test for large-scale non-Gaussianities in the cosmic microwave background (CMB) in a model-independent way. To this end, so-called first and second order surrogates are generated by first shuffling the Fourier phases belonging to the scales not of interest and then shuffling the remaining phases for the length scales under study. Using scaling indices as test statistics we find highly significant signatures for both non-Gaussianities and asymmetries on large scales for the WMAP data of the CMB. We find remarkably similar results when analyzing different ILC-maps based on the WMAP five and seven year data. Such features being independent from the map-making procedure would disfavor the fundamental principle of isotropy as well as canonical single-field slow-roll inflation - unless there is some undiscovered systematic error in the collection or reduction of the CMB data or yet unknown foreground contributions.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Increasing evidence for hemispherical power asymmetry in the five-year WMAP data

(Abridged)Motivated by the recent results of Hansen et al. (2008) concerning a noticeable hemispherical power asymmetry in the WMAP data on small angular scales, we revisit the dipole modulated signal model introduced by Gordon et al. (2005). This model assumes that the true CMB signal consists of a Gaussian isotropic random field modulated by a dipole, and is characterized by an overall modulation amplitude, A, and a preferred direction, p. Previous analyses of this model has been restricted to very low resolution due to computational cost. In this paper, we double the angular resolution, and compute the full corresponding posterior distribution for the 5-year WMAP data. The results from our analysis are the following: The best-fit modulation amplitude for l <= 64 and the ILC data with the WMAP KQ85 sky cut is A=0.072 +/- 0.022, non-zero at 3.3sigma, and the preferred direction points toward Galactic coordinates (l,b) = (224 degree, -22 degree) +/- 24 degree. The corresponding results for l <~ 40 from earlier analyses was A = 0.11 +/- 0.04 and (l,b) = (225 degree,-27 degree). The statistical significance of a non-zero amplitude thus increases from 2.8sigma to 3.3sigma when increasing l_max from 40 to 64, and all results are consistent to within 1sigma. Similarly, the Bayesian log-evidence difference with respect to the isotropic model increases from Delta ln E = 1.8 to Delta ln E = 2.6, ranking as &#34;strong evidence&#34; on the Jeffreys&#39; scale. The raw best-fit log-likelihood difference increases from Delta ln L = 6.1 to Delta ln L = 7.3. Similar, and often slightly stronger, results are found for other data combinations. Thus, we find that the evidence for a dipole power distribution in the WMAP data increases with l in the 5-year WMAP data set, in agreement with the reports of Hansen et al. (2008).

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Residual noise covariance for Planck low-resolution data analysis

Aims: Develop and validate tools to estimate residual noise covariance in Planck frequency maps. Quantify signal error effects and compare different techniques to produce low-resolution maps. Methods: We derive analytical estimates of covariance of the residual noise contained in low-resolution maps produced using a number of map-making approaches. We test these analytical predictions using Monte Carlo simulations and their impact on angular power spectrum estimation. We use simulations to quantify the level of signal errors incurred in different resolution downgrading schemes considered in this work. Results: We find an excellent agreement between the optimal residual noise covariance matrices and Monte Carlo noise maps. For destriping map-makers, the extent of agreement is dictated by the knee frequency of the correlated noise component and the chosen baseline offset length. The significance of signal striping is shown to be insignificant when properly dealt with. In map resolution downgrading, we find that a carefully selected window function is required to reduce aliasing to the sub-percent level at multipoles, ell > 2Nside, where Nside is the HEALPix resolution parameter. We show that sufficient characterization of the residual noise is unavoidable if one is to draw reliable contraints on large scale anisotropy. Conclusions: We have described how to compute the low-resolution maps, with a controlled sky signal level, and a reliable estimate of covariance of the residual noise. We have also presented a method to smooth the residual noise covariance matrices to describe the noise correlations in smoothed, bandwidth limited maps.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

A Markov Chain Monte Carlo Algorithm for analysis of low signal-to-noise CMB data

We present a new Monte Carlo Markov Chain algorithm for CMB analysis in the low signal-to-noise regime. This method builds on and complements the previously described CMB Gibbs sampler, and effectively solves the low signal-to-noise inefficiency problem of the direct Gibbs sampler. The new algorithm is a simple Metropolis-Hastings sampler with a general proposal rule for the power spectrum, C_l, followed by a particular deterministic rescaling operation of the sky signal. The acceptance probability for this joint move depends on the sky map only through the difference of chi-squared between the original and proposed sky sample, which is close to unity in the low signal-to-noise regime. The algorithm is completed by alternating this move with a standard Gibbs move. Together, these two proposals constitute a computationally efficient algorithm for mapping out the full joint CMB posterior, both in the high and low signal-to-noise regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Bayesian analysis of the low-resolution polarized 3-year WMAP sky maps

We apply a previously developed Gibbs sampling framework to the foreground corrected 3-yr WMAP polarization data and compute the power spectrum and residual foreground template amplitude posterior distributions. We first analyze the co-added Q- and V-band data, and compare our results to the likelihood code published by the WMAP team. We find good agreement, and thus verify the numerics and data processing steps of both approaches. However, we also analyze the Q- and V-bands separately, allowing for non-zero EB cross-correlations and including two individual foreground template amplitudes tracing synchrotron and dust emission. In these analyses, we find tentative evidence of systematics: The foreground tracers correlate with each of the Q- and V-band sky maps individually, although not with the co-added QV map; there is a noticeable negative EB cross-correlation at l <~ 16 in the V-band map; and finally, when relaxing the constraints on EB and BB, noticeable differences are observed between the marginalized band powers in the Q- and V-bands. Further studies of these features are imperative, given the importance of the low-l EE spectrum on the optical depth of reionization tau and the spectral index of scalar perturbations n_s.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Joint Bayesian component separation and CMB power spectrum estimation

We describe and implement an exact, flexible, and computationally efficient algorithm for joint component separation and CMB power spectrum estimation, building on a Gibbs sampling framework. Two essential new features are 1) conditional sampling of foreground spectral parameters, and 2) joint sampling of all amplitude-type degrees of freedom (e.g., CMB, foreground pixel amplitudes, and global template amplitudes) given spectral parameters. Given a parametric model of the foreground signals, we estimate efficiently and accurately the exact joint foreground-CMB posterior distribution, and therefore all marginal distributions such as the CMB power spectrum or foreground spectral index posteriors. The main limitation of the current implementation is the requirement of identical beam responses at all frequencies, which restricts the analysis to the lowest resolution of a given experiment. We outline a future generalization to multi-resolution observations. To verify the method, we analyse simple models and compare the results to analytical predictions. We then analyze a realistic simulation with properties similar to the 3-yr WMAP data, downgraded to a common resolution of 3 degree FWHM. The results from the actual 3-yr WMAP temperature analysis are presented in a companion Letter.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The joint large-scale foreground-CMB posteriors of the 3-year WMAP data

Using a Gibbs sampling algorithm for joint CMB estimation and component separation, we compute the large-scale CMB and foreground posteriors of the 3-yr WMAP temperature data. Our parametric data model includes the cosmological CMB signal and instrumental noise, a single power law foreground component with free amplitude and spectral index for each pixel, a thermal dust template with a single free overall amplitude, and free monopoles and dipoles at each frequency. This simple model yields a surprisingly good fit to the data over the full frequency range from 23 to 94 GHz. We obtain a new estimate of the CMB sky signal and power spectrum, and a new foreground model, including a measurement of the effective spectral index over the high-latitude sky. A particularly significant result is the detection of a common spurious offset in all frequency bands of ~ -13muK, as well as a dipole in the V-band data. Correcting for these is essential when determining the effective spectral index of the foregrounds. We find that our new foreground model is in good agreement with template-based model presented by the WMAP team, but not with their MEM reconstruction. We believe the latter may be at least partially compromised by the residual offsets and dipoles in the data. Fortunately, the CMB power spectrum is not significantly affected by these issues, as our new spectrum is in excellent agreement with that published by the WMAP team. The corresponding cosmological parameters are also virtually unchanged.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere

HEALPix -- the Hierarchical Equal Area iso-Latitude Pixelization -- is a versatile data structure with an associated library of computational algorithms and visualization software that supports fast scientific applications executable directly on very large volumes of astronomical data and large area surveys in the form of discretized spherical maps. Originally developed to address the data processing and analysis needs of the present generation of cosmic microwave background (CMB) experiments (e.g. BOOMERanG, WMAP), HEALPix can be expanded to meet many of the profound challenges that will arise in confrontation with the observational output of future missions and experiments, including e.g. Planck, Herschel, SAFIR, and the Beyond Einstein CMB polarization probe. In this paper we consider the requirements and constraints to be met in order to implement a sufficient framework for the efficient discretization and fast analysis/synthesis of functions defined on the sphere, and summarise how they are satisfied by HEALPix.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Spherical Mexican Hat wavelet: an application to detect non-Gaussianity in the COBE-DMR maps

The spherical Mexican Hat wavelet is introduced in this paper, with the aim of testing the Gaussianity of the Cosmic Microwave Background temperature fluctuations. Using the information given by the wavelet coefficients at several scales, we have performed several statistical tests on the COBE-DMR maps to search for evidence of non-Gaussianity. Skewness, kurtosis, scale-scale correlations (for two and three scales) and Kolmogorov-Smirnov tests indicate that the COBE-DMR data are consistent with a Gaussian distribution. We have extended the analysis to compare temperature values provided by COBE-DMR data with distributions (obtained from Gaussian simulations) at each pixel and at each scale. The number of pixels with temperature values outside the 95% and the 99% is consistent with that obtained from Gaussian simulations, at all scales. Moreover, the extrema values for COBE-DMR data (maximum and minimum temperatures in the map) are also consistent with those obtained from Gaussian simulations.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

A Bayesian estimate of the skewness of the Cosmic Microwave Background

We propose a formalism for estimating the skewness and angular power spectrum of a general Cosmic Microwave Background data set. We use the Edgeworth Expansion to define a non-Gaussian likelihood function that takes into account the anisotropic nature of the noise and the incompleteness of the sky coverage. The formalism is then applied to estimate the skewness of the publicly available 4 year Cosmic Background Explorer (COBE) Differential Microwave Radiometer data. We find that the data is consistent with a Gaussian skewness, and with isotropy. Inclusion of non Gaussian degrees of freedom has essentially no effect on estimates of the power spectrum, if each $C_\ell$ is regarded as a separate parameter or if the angular power spectrum is parametrized in terms of an amplitude (Q) and spectral index (n). Fixing the value of the angular power spectrum at its maxiumum likelihood estimate, the best fit skewness is $S=6.5\pm6.0\times10^4(\muK)^3$; marginalizing over Q the estimate of the skewness is $S=6.5\pm8.4\times10^4(\muK)^3$ and marginalizing over n one has $S=6.5\pm8.5\times10^4(\muK)^3$.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

On the Non-Gaussianity Observed in the COBE-DMR Sky Maps

In this paper we pursue the origin of the non-Gaussianity determined by a bispectrum analysis of the COBE-DMR 4-year sky maps. The robustness of the statistic is demonstrated by the rebinning of the data into 12 coordinate systems. By computing the bispectrum statistic as a function of various data partitions - by channel, frequency, and time interval, we show that the observed non-Gaussian signal is driven by the 53 GHz data. This frequency dependence strongly rejects the hypothesis that the signal is cosmological in origin. A jack-knife analysis of the coadded 53 and 90 GHz sky maps reveals those sky pixels to which the bispectrum statistic is particularly sensitive. We find that by removing data from the 53 GHz sky maps for periods of time during which a known systematic effect perturbs the 31 GHz channels, the amplitudes of the bispectrum coefficients become completely consistent with that expected for a Gaussian sky. We conclude that the non-Gaussian signal detected by the normalised bispectrum statistic in the publicly available DMR sky maps is due to a systematic artifact. The impact of removing the affected data on estimates of the normalisation of simple models of cosmological anisotropy is negligible.

0 citations0 reviewsNo code linkedOpen access