We present a method for fast optimal estimation of the temperature angular power spectrum from observations of the cosmic microwave background. We employ a Hamiltonian Monte Carlo (HMC) sampler to ...obtain samples from the posterior probability distribution of all the power spectrum coefficients given a set of observations. We compare the properties of the HMC and the related Gibbs sampling approach on low-resolution simulations. The correlation lengths of the samples produced by Gibbs and HMC are generally comparable. At very high signal-to-noise (≳10) the Gibbs outperforms the HMC while for intermediate signal-to-noise (∼0.1–10) the HMC results in shorter correlation lengths than the Gibbs. At low signal-to-noise (≲0.1) both samplers appear to perform similarly. We also demonstrate the method on high-resolution data by applying it to simulated Wilkinson Microwave Anisotropy Probe (WMAP) data. Analysis of a single map from a WMAP-resolution data set is possible in around 80 h on a high-end desktop computer. HMC uses around a factor of 3–4 fewer spherical harmonic transforms than is reported for a Gibbs sampler (with pre-conditioning) and therefore produces noticeable performance gains. HMC imposes few conditions on the distribution to be sampled and obviates the requirement of having to be able to sample efficiently from the conditional densities, thus providing us with an extremely flexible approach upon which to build.
We present the Planck Sky Model (PSM), a parametric model for generating all-sky, few arcminute resolution maps of sky emission at submillimetre to centimetre wavelengths, in both intensity and ...polarisation. Several options are implemented to model the cosmic microwave background, Galactic diffuse emission (synchrotron, free-free, thermal and spinning dust, CO lines), Galactic H ii regions, extragalactic radio sources, dusty galaxies, and thermal and kinetic Sunyaev-Zeldovich signals from clusters of galaxies. Each component is simulated by means of educated interpolations/extrapolations of data sets available at the time of the launch of the Planck mission, complemented by state-of-the-art models of the emission. Distinctive features of the simulations are spatially varying spectral properties of synchrotron and dust; different spectral parameters for each point source; modelling of the clustering properties of extragalactic sources and of the power spectrum of fluctuations in the cosmic infrared background. The PSM enables the production of random realisations of the sky emission, constrained to match observational data within their uncertainties. It is implemented in a software package that is regularly updated with incoming information from observations. The model is expected to serve as a useful tool for optimising planned microwave and sub-millimetre surveys and testing data processing and analysis pipelines. It is, in particular, used to develop and validate data analysis pipelines within the Planck collaboration. A version of the software that can be used for simulating the observations for a variety of experiments is made available on a dedicated website.
Precise measurement of the angular power spectrum of the cosmic microwave background (CMB) temperature and polarization anisotropy can tightly constrain many cosmological models and parameters. ...However, accurate measurements can only be realized in practice provided all major systematic effects have been taken into account. Beam asymmetry, coupled with the scan strategy, is a major source of systematic error in scanning CMB experiments such as Planck, the focus of our current interest. We envision Monte Carlo methods to rigorously study and account for the systematic effect of beams in CMB analysis. Toward that goal, we have developed a fast pixel space convolution method that can simulate sky maps observed by a scanning instrument, taking into account real beam shapes and scan strategy. The essence is to pre-compute the 'effective beams' using a computer code, 'Fast Effective Beam Convolution in Pixel space' (FEBeCoP), that we have developed for the Planck mission. The code computes effective beams given the focal plane beam characteristics of the Planck instrument and the full history of actual satellite pointing, and performs very fast convolution of sky signals using the effective beams. In this paper, we describe the algorithm and the computational scheme that has been implemented. We also outline a few applications of the effective beams in the precision analysis of Planck data, for characterizing the CMB anisotropy and for detecting and measuring properties of point sources.
ABSTRACT
We present QUIJOTE intensity and polarization maps in four frequency bands centred around 11, 13, 17, and 19 GHz, and covering approximately 29 000 deg2, including most of the northern sky ...region. These maps result from 9000 h of observations taken between May 2013 and June 2018 with the first QUIJOTE multifrequency instrument (MFI), and have angular resolutions of around 1°, and sensitivities in polarization within the range 35–40 µK per 1° beam, being a factor ∼2–4 worse in intensity. We discuss the data processing pipeline employed, and the basic characteristics of the maps in terms of real space statistics and angular power spectra. A number of validation tests have been applied to characterize the accuracy of the calibration and the residual level of systematic effects, finding a conservative overall calibration uncertainty of 5 per cent. We also discuss flux densities for four bright celestial sources (Tau A, Cas A, Cyg A, and 3C274), which are often used as calibrators at microwave frequencies. The polarization signal in our maps is dominated by synchrotron emission. The distribution of spectral index values between the 11 GHz and WMAP 23 GHz map peaks at β = −3.09 with a standard deviation of 0.14. The measured BB/EE ratio at scales of ℓ = 80 is 0.26 ± 0.07 for a Galactic cut |b| > 10°. We find a positive TE correlation for 11 GHz at large angular scales (ℓ ≲ 50), while the EB and TB signals are consistent with zero in the multipole range 30 ≲ ℓ ≲ 150. The maps discussed in this paper are publicly available.
Planck intermediate results Ade, P A R; Aghanim, N; Arnaud, M ...
Astronomy and astrophysics (Berlin),
08/2015, Letnik:
580
Journal Article
Recenzirano
Odprti dostop
Any variation in the fundamental physical constants, more particularly in the fine structure constant, alpha, or in the mass of the electron, m sub(e), affects the recombination history of the ...Universe and cause an imprint on the cosmic microwave background angular power spectra. We show that the Planck data allow one to improve the constraint on the time variation of the fine structure constant at redshift z ~ 10 super(3) by about a factor of 5 compared to WMAP data, as well as to break the degeneracy with the Hubble constant, H sub(0). In addition to alpha, we can set a constraint on the variation in the mass of the electron, m sub(e), and in the simultaneous variation of the two constants. We examine in detail the degeneracies between fundamental constants and the cosmological parameters, in order to compare the limits obtained from Planck and WMAP and to determine the constraining power gained by including other cosmological probes. We conclude that independent time variations of the fine structure constant and of the mass of the electron are constrained by Planck to Deltaalpha/alpha = (3.6 + or - 3.7) x 10 super(-3) and Deltam sub(e)/m sub(e) = (4 + or - 11) x 10 super(-3) at the 68% confidence level. We also investigate the possibility of a spatial variation of the fine structure constant. The relative amplitude of a dipolar spatial variation in a (corresponding to a gradient across our Hubble volume) is constrained to be deltaalpha/alpha = (-2.4 + or - 3.7) x 10 super(-2).
ABSTRACT
Map-making is an important step for the data analysis of cosmic microwave background (CMB) experiments. It consists of converting the data, which are typically a long, complex, and noisy ...collection of measurements, into a map, which is an image of the observed sky. We present in this paper a new map-making code named picasso (Polarization and Intensity CArtographer for Scanned Sky Observations), which was implemented to construct intensity and polarization maps from the Multi Frequency Instrument (MFI) of the QUIJOTE (Q-U-I Joint TEnerife) CMB polarization experiment. picasso is based on the destriping algorithm, and is suited to address specific issues of ground-based microwave observations, with a technique that allows the fit of a template function in the time domain, during the map-making step. This paper describes the picasso code, validating it with simulations and assessing its performance. For this purpose, we produced realistic simulations of the QUIJOTE-MFI survey of the northern sky (approximately ∼20 000 deg2), and analysed the reconstructed maps with picasso, using real and harmonic space statistics. We show that, for this sky area, picasso is able to reconstruct, with high fidelity, the injected signal, recovering all the scales with ℓ > 10 in TT, EE, and BB. The signal error is better than 0.001 per cent at 20 < ℓ < 200. Finally, we validated some of the methods that will be applied to the real wide-survey data, like the detection of the CMB anisotropies via cross-correlation analyses. Despite that the implementation of picasso is specific for QUIJOTE-MFI data, it could be adapted to other experiments.
In this paper, we present Q-U-I JOint Tenerife Experiment (QUIJOTE) 10–20 GHz observations (194 h in total over ≈250 deg2) in intensity and polarisation of G159.6-18.5, one of the most widely studied ...regions harbouring anomalous microwave emission (AME). By combining with other publicly available intensity data, we achieve the most precise spectrum of the AME measured to date in an individual region, with 13 independent data points between 10 and 50 GHz being dominated by this emission. The four QUIJOTE data points provide the first independent confirmation of the downturn of the AME spectrum at low frequencies, initially unveiled by the COSMOlogical Structures On Medium Angular Scales experiment in this region. Our polarisation maps, which have an angular resolution of ≈1° and a sensitivity of ≈ 25 μK beam−1, are consistent with zero polarisation. We obtain upper limits on the polarisation fraction of Π < 6.3 and <2.8 per cent (95 per cent C.L.), respectively, at 12 and 18 GHz (ΠAME < 10.1 and <3.4 per cent with respect to the residual AME intensity), a frequency range where no AME polarisation observations have been reported to date. The combination of these constraints with those from other experiments confirm that all the magnetic dust models based on single-domain grains, and most of those considering randomly oriented magnetic inclusions, predict higher polarisation levels than is observed towards regions with AME. Also, neither of the two considered models of electric dipole emission seems to be compatible with all the observations together. More stringent constraints of the AME polarisation at 10–40 GHz are necessary to disentangle between different models, to which future QUIJOTE data will contribute.
ABSTRACT
The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be ...correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarize the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterize 233 245 neutral hydrogen (H i) sources in a simulated data product representing a 2000 h SKA-Mid spectral line observation from redshifts 0.25–0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, ‘reproducibility awards’ were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy – which combined predictions from two independent machine learning techniques to yield a 20 per cent improvement in overall performance – underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical data sets.
ABSTRACT
Anomalous Microwave Emission (AME) is a significant component of Galactic diffuse emission in the frequency range 10–$60\, \mathrm{GHz}$ and a new window into the properties of ...sub-nanometre-sized grains in the interstellar medium. We investigate the morphology of AME in the ≈10○ diameter λ Orionis ring by combining intensity data from the QUIJOTE experiment at 11, 13, 17, and $19\, \mathrm{GHz}$ and the C-Band All Sky Survey (C-BASS) at $4.76\, \mathrm{GHz}$, together with 19 ancillary data sets between 1.42 and $3000\, \mathrm{GHz}$. Maps of physical parameters at 1○ resolution are produced through Markov chain Monte Carlo (MCMC) fits of spectral energy distributions (SEDs), approximating the AME component with a lognormal distribution. AME is detected in excess of $20\, \sigma$ at degree-scales around the entirety of the ring along photodissociation regions (PDRs), with three primary bright regions containing dark clouds. A radial decrease is observed in the AME peak frequency from $\approx 35\, \mathrm{GHz}$ near the free–free region to $\approx 21\, \mathrm{GHz}$ in the outer regions of the ring, which is the first detection of AME spectral variations across a single region. A strong correlation between AME peak frequency, emission measure and dust temperature is an indication for the dependence of the AME peak frequency on the local radiation field. The AME amplitude normalized by the optical depth is also strongly correlated with the radiation field, giving an overall picture consistent with spinning dust where the local radiation field plays a key role.