150 research outputs found
A Monte Carlo comparison between template-based and Wiener-filter CMB dipole estimators
We review and compare two different CMB dipole estimators discussed in the
literature, and assess their performances through Monte Carlo simulations. The
first method amounts to simple template regression with partial sky data, while
the second method is an optimal Wiener filter (or Gibbs sampling)
implementation. The main difference between the two methods is that the latter
approach takes into account correlations with higher-order CMB temperature
fluctuations that arise from non-orthogonal spherical harmonics on an
incomplete sky, which for recent CMB data sets (such as Planck) is the dominant
source of uncertainty. For an accepted sky fraction of 81% and an angular CMB
power spectrum corresponding to the best-fit Planck 2018 CDM model, we
find that the uncertainty on the recovered dipole amplitude is about six times
smaller for the Wiener filter approach than for the template approach,
corresponding to 0.5 and 3K, respectively. Similar relative differences
are found for the corresponding directional parameters and other sky fractions.
We note that the Wiener filter algorithm is generally applicable to any dipole
estimation problem on an incomplete sky, as long as a statistical and
computationally tractable model is available for the unmasked higher-order
fluctuations. The methodology described in this paper forms the numerical basis
for the most recent determination of the CMB solar dipole from Planck, as
summarized by arXiv:2007.04997.Comment: 8 pages, 10 figures, submitted to A&
Cosmoglobe DR1. III. First full-sky model of polarized synchrotron emission from all WMAP and Planck LFI data
We present the first model of full-sky polarized synchrotron emission that is
derived from all WMAP and Planck LFI frequency maps. The basis of this analysis
is the set of end-to-end reprocessed Cosmoglobe Data Release 1 sky maps
presented in a companion paper, which have significantly lower instrumental
systematics than the legacy products from each experiment. We find that the
resulting polarized synchrotron amplitude map has an average noise rms of
at 30 GHz and FWHM, which is 30% lower than
the recently released BeyondPlanck model that included only LFI+WMAP Ka-V data,
and 29% lower than the WMAP K-band map alone. The mean -to- power
spectrum ratio is , with amplitudes consistent with those measured
previously by Planck and QUIJOTE. Assuming a power law model for the
synchrotron spectral energy distribution, and using the -- plot method,
we find a full-sky inverse noise-variance weighted mean of
between Cosmoglobe DR1 K-band and 30 GHz, in
good agreement with previous estimates. In summary, the novel Cosmoglobe DR1
synchrotron model is both more sensitive and systematically cleaner than
similar previous models, and it has a more complete error description that is
defined by a set of Monte Carlo posterior samples. We believe that these
products are preferable over previous Planck and WMAP products for all
synchrotron-related scientific applications, including simulation, forecasting
and component separation.Comment: 15 pages, 15 figures, submitted to A&
Cosmoglobe DR1 results. II. Constraints on isotropic cosmic birefringence from reprocessed WMAP and Planck LFI data
Cosmic birefringence is a parity-violating effect that might have rotated the
plane of linearly polarized light of the cosmic microwave background (CMB) by
an angle since its emission. This has recently been measured to be
non-zero at a statistical significance of in the official Planck
PR4 and 9-year WMAP data. In this work, we constrain using the
reprocessed BeyondPlanck LFI and Cosmoglobe DR1 WMAP polarization maps. These
novel maps have both lower systematic residuals and a more complete error
description than the corresponding official products. Foreground
correlations could bias measurements of , and while thermal dust
emission has been argued to be statistically non-zero, no evidence for
synchrotron power has been reported. Unlike the dust-dominated Planck HFI
maps, the majority of the LFI and WMAP polarization maps are instead dominated
by synchrotron emission. Simultaneously constraining and the
polarization miscalibration angle, , of each channel, we find a
best-fit value of with LFI and WMAP data
only. When including the Planck HFI PR4 maps, but fitting separately
for dust-dominated, , and synchrotron-dominated
channels, , we find . This differs from zero with a
statistical significance of , and the main contribution to this
value comes from the LFI 70 GHz channel. While the statistical significances of
these results are low on their own, the measurement derived from the LFI and
WMAP synchrotron-dominated maps agrees with the previously reported
HFI-dominated constraints, despite the very different astrophysical and
instrumental systematics involved in all these experiments.Comment: 10 pages, 7 figures, 2 tables. Submitted to A&
Cosmoglobe DR1 results. I. Improved Wilkinson Microwave Anisotropy Probe maps through Bayesian end-to-end analysis
We present Cosmoglobe Data Release 1, which implements the first joint
analysis of WMAP and Planck LFI time-ordered data, processed within a single
Bayesian end-to-end framework. This framework builds directly on a similar
analysis of the LFI measurements by the BeyondPlanck collaboration, and
approaches the CMB analysis challenge through Gibbs sampling of a global
posterior distribution, simultaneously accounting for calibration, mapmaking,
and component separation. The computational cost of producing one complete
WMAP+LFI Gibbs sample is 812 CPU-hr, of which 603 CPU-hrs are spent on WMAP
low-level processing; this demonstrates that end-to-end Bayesian analysis of
the WMAP data is computationally feasible. We find that our WMAP posterior mean
temperature sky maps and CMB temperature power spectrum are largely consistent
with the official WMAP9 results. Perhaps the most notable difference is that
our CMB dipole amplitude is , which is $11\
\mathrm{\mu K}2.5\ {\sigma}$ higher than
BeyondPlanck; however, it is in perfect agreement with the HFI-dominated Planck
PR4 result. In contrast, our WMAP polarization maps differ more notably from
the WMAP9 results, and in general exhibit significantly lower large-scale
residuals. We attribute this to a better constrained gain and transmission
imbalance model. It is particularly noteworthy that the W-band polarization sky
map, which was excluded from the official WMAP cosmological analysis, for the
first time appears visually consistent with the V-band sky map. Similarly, the
long standing discrepancy between the WMAP K-band and LFI 30 GHz maps is
finally resolved, and the difference between the two maps appears consistent
with instrumental noise at high Galactic latitudes. All maps and the associated
code are made publicly available through the Cosmoglobe web page.Comment: 65 pages, 61 figures. Data available at cosmoglobe.uio.no. Submitted
to A&
Molecular and functional interactions between tumor necrosis factor-alpha receptors and the glutamatergic system in the mouse hippocampus : implications for seizure susceptibility
Tumor necrosis factor (TNF)-alpha is a proinflammatory cytokine acting on two distinct receptor subtypes, namely p55 and p75 receptors. TNF-alpha p55 and p75 receptor knockout mice were previously shown to display a decreased or enhanced susceptibility to seizures, respectively, suggesting intrinsic modifications in neuronal excitability. We investigated whether alterations in glutamate system function occur in these naive knockout mice with perturbed cytokine signaling that could explain their different propensity to develop seizures. Using Western blot analysis of hippocampal homogenates, we found that p55(-/-) mice have decreased levels of membrane GluR3 and NR1 glutamate receptor subunits while GluR1, GluR2, GluR6/7 and NR2A/B were unchanged as compared to wild-type mice. In p75(-/-) mice, GluR2, GluR3, GluR6/7 and NR2A/B glutamate receptor subunits were increased in the hippocampus while GluR1 and NR1 did not change. Extracellular single-cell recordings of the electrical activity of hippocampal neurons were carried out in anesthetized mice by standard electrophysiological techniques. Microiontophoretic application of glutamate increased the basal firing rate of hippocampal neurons in p75(-/-) mice versus wild-type mice, and this effect was blocked by 2-amino-5-phosphopentanoic acid and 6-nitro-7-sulfamoyl-benzo(f)quinoxaline-2,3-dione denoting the involvement of N-methyl-D-aspartic acid and AMPA receptors. In p55(-/-) mice, hippocampal neurons responses to glutamate were similar to wild-type mice. Spontaneous glutamate release measured by in vivo hippocampal microdialysis was significantly decreased only in p55(-/-) mice. No changes were observed in KCl-induced glutamate release in both receptor knockout mice strains versus wild-type mice. These findings highlight specific molecular and functional interactions between p55 and p75 receptor-mediated signaling and the glutamate system. These interactions may be relevant for controlling neuronal excitability in physiological and pathological conditions.peer-reviewe
BeyondPlanck IV. On end-to-end simulations in CMB analysis -- Bayesian versus frequentist statistics
End-to-end simulations play a key role in the analysis of any
high-sensitivity CMB experiment, providing high-fidelity systematic error
propagation capabilities unmatched by any other means. In this paper, we
address an important issue regarding such simulations, namely how to define the
inputs in terms of sky model and instrument parameters. These may either be
taken as a constrained realization derived from the data, or as a random
realization independent from the data. We refer to these as Bayesian and
frequentist simulations, respectively. We show that the two options lead to
significantly different correlation structures, as frequentist simulations,
contrary to Bayesian simulations, effectively include cosmic variance, but
exclude realization-specific correlations from non-linear degeneracies.
Consequently, they quantify fundamentally different types of uncertainties, and
we argue that they therefore also have different and complementary scientific
uses, even if this dichotomy is not absolute. Before BeyondPlanck, most
pipelines have used a mix of constrained and random inputs, and used the same
hybrid simulations for all applications, even though the statistical
justification for this is not always evident. BeyondPlanck represents the first
end-to-end CMB simulation framework that is able to generate both types of
simulations, and these new capabilities have brought this topic to the
forefront. The Bayesian BeyondPlanck simulations and their uses are described
extensively in a suite of companion papers. In this paper we consider one
important applications of the corresponding frequentist simulations, namely
code validation. That is, we generate a set of 1-year LFI 30 GHz frequentist
simulations with known inputs, and use these to validate the core low-level
BeyondPlanck algorithms; gain estimation, correlated noise estimation, and
mapmaking
BeyondPlanck II. CMB map-making through Gibbs sampling
We present a Gibbs sampling solution to the map-making problem for CMB
measurements, building on existing destriping methodology. Gibbs sampling
breaks the computationally heavy destriping problem into two separate steps;
noise filtering and map binning. Considered as two separate steps, both are
computationally much cheaper than solving the combined problem. This provides a
huge performance benefit as compared to traditional methods, and allows us for
the first time to bring the destriping baseline length to a single sample. We
apply the Gibbs procedure to simulated Planck 30 GHz data. We find that gaps in
the time-ordered data are handled efficiently by filling them with simulated
noise as part of the Gibbs process. The Gibbs procedure yields a chain of map
samples, from which we may compute the posterior mean as a best-estimate map.
The variation in the chain provides information on the correlated residual
noise, without need to construct a full noise covariance matrix. However, if
only a single maximum-likelihood frequency map estimate is required, we find
that traditional conjugate gradient solvers converge much faster than a Gibbs
sampler in terms of total number of iterations. The conceptual advantages of
the Gibbs sampling approach lies in statistically well-defined error
propagation and systematic error correction, and this methodology forms the
conceptual basis for the map-making algorithm employed in the BeyondPlanck
framework, which implements the first end-to-end Bayesian analysis pipeline for
CMB observations.Comment: 11 pages, 10 figures. All BeyondPlanck products and software will be
released publicly at http://beyondplanck.science during the online release
conference (November 18-20, 2020). Connection details will be made available
at the same website. Registration is mandatory for the online tutorial, but
optional for the conferenc
BeyondPlanck VII. Bayesian estimation of gain and absolute calibration for CMB experiments
We present a Bayesian calibration algorithm for CMB observations as
implemented within the global end-to-end BeyondPlanck (BP) framework, and apply
this to the Planck Low Frequency Instrument (LFI) data. Following the most
recent Planck analysis, we decompose the full time-dependent gain into a sum of
three orthogonal components: One absolute calibration term, common to all
detectors; one time-independent term that can vary between detectors; and one
time-dependent component that is allowed to vary between one-hour pointing
periods. Each term is then sampled conditionally on all other parameters in the
global signal model through Gibbs sampling. The absolute calibration is sampled
using only the orbital dipole as a reference source, while the two relative
gain components are sampled using the full sky signal, including the orbital
and Solar CMB dipoles, CMB fluctuations, and foreground contributions. We
discuss various aspects of the data that influence gain estimation, including
the dipole/polarization quadrupole degeneracy and anomalous jumps in the
instrumental gain. Comparing our solution to previous pipelines, we find good
agreement in general, with relative deviations of -0.84% (-0.67%) for 30 GHz,
-0.14% (0.02%) for 44 GHz and -0.69% (-0.08%) for 70 GHz, compared to Planck
2018 (NPIPE). The deviations we find are within expected error bounds, and we
attribute them to differences in data usage and general approach between the
pipelines. In particular, the BP calibration is performed globally, resulting
in better inter-frequency consistency. Additionally, WMAP observations are used
actively in the BP analysis, which breaks degeneracies in the Planck data set
and results in better agreement with WMAP. Although our presentation and
algorithm are currently oriented toward LFI processing, the procedure is fully
generalizable to other experiments.Comment: 18 pages, 15 figures. All BeyondPlanck products and software will be
released publicly at http://beyondplanck.science during the online release
conference (November 18-20, 2020). Connection details will be made available
at the same website. Registration is mandatory for the online tutorial, but
optional for the conferenc
BeyondPlanck XII. Cosmological parameter constraints with end-to-end error propagation
We present cosmological parameter constraints as estimated using the Bayesian
BeyondPlanck (BP) analysis framework. This method supports seamless end-to-end
error propagation from raw time-ordered data to final cosmological parameters.
As a first demonstration of the method, we analyze time-ordered Planck LFI
observations, combined with selected external data (WMAP 33-61GHz, Planck HFI
DR4 353 and 857GHz, and Haslam 408MHz) in the form of pixelized maps which are
used to break critical astrophysical degeneracies. Overall, all results are
generally in good agreement with previously reported values from Planck 2018
and WMAP, with the largest relative difference for any parameter of about 1
sigma when considering only temperature multipoles between 29<l<601. In cases
where there are differences, we note that the BP results are generally slightly
closer to the high-l HFI-dominated Planck 2018 results than previous analyses,
suggesting slightly less tension between low and high multipoles. Using low-l
polarization information from LFI and WMAP, we find a best-fit value of
tau=0.066 +/- 0.013, which is higher than the low value of tau=0.051 +/- 0.006
derived from Planck 2018 and slightly lower than the value of 0.069 +/- 0.011
derived from joint analysis of official LFI and WMAP products. Most
importantly, however, we find that the uncertainty derived in the BP processing
is about 30% larger than when analyzing the official products, after taking
into account the different sky coverage. We argue that this is due to
marginalizing over a more complete model of instrumental and astrophysical
parameters, and this results in both more reliable and more rigorously defined
uncertainties. We find that about 2000 Monte Carlo samples are required to
achieve robust convergence for low-resolution CMB covariance matrix with 225
independent modes.Comment: 13 pages, 10 figure
BeyondPlanck XIV. Polarized foreground emission between 30 and 70GHz
We constrain polarized foreground emission between 30 and 70GHz with the
Planck Low Frequency Instrument (LFI) and WMAP data within the framework of
BeyondPlanck global Bayesian analysis. We combine for the first time
full-resolution Planck LFI time-ordered data with low-resolution WMAP sky maps
at 33, 40 and 61GHz. Spectral parameters are fit with a likelihood defined at
the native resolution of each frequency channel. This analysis represents the
first implementation of true multi-resolution component separation applied to
CMB observations for both amplitude and spectral energy distribution (SED)
parameters. For synchrotron emission, we approximate the SED as a power-law in
frequency and find that the low signal-to-noise ratio of the data set strongly
limits the number of free parameters that may be robustly constrained. We
partition the sky into four large disjoint regions (High Latitude; Galactic
Spur; Galactic Plane; and Galactic Center), each associated with its own
power-law index. We find that the High Latitude region is prior-dominated,
while the Galactic Center region is contaminated by residual instrumental
systematics. The two remaining regions appear to be both signal-dominated and
clean of systematics, and for these we derive spectral indices of
and . This agrees qualitatively with the WMAP-only
polarization constraints presented by Dunkley et al. (2009), but contrasts with
several temperature-based analyses. For thermal dust emission we assume a
modified blackbody model and we fit the power-law index across the full sky. We
find , which is slightly steeper than that
previously reported from Planck HFI data, but still statistically consistent at
a 2 confidence level.Comment: 17 pages, 14 figures. All BeyondPlanck products and software will be
released publicly at http://beyondplanck.science during the online release
conference (November 18-20, 2020). Connection details will be made available
at the same website. Registration is mandatory for the online tutorial, but
optional for the conferenc
- …