-
Deep Learning as Ricci Flow
Authors:
Anthony Baptista,
Alessandro Barp,
Tapabrata Chakraborti,
Chris Harbron,
Ben D. MacArthur,
Christopher R. S. Banerji
Abstract:
Deep neural networks (DNNs) are powerful tools for approximating the distribution of complex data. It is known that data passing through a trained DNN classifier undergoes a series of geometric and topological simplifications. While some progress has been made toward understanding these transformations in neural networks with smooth activation functions, an understanding in the more general settin…
▽ More
Deep neural networks (DNNs) are powerful tools for approximating the distribution of complex data. It is known that data passing through a trained DNN classifier undergoes a series of geometric and topological simplifications. While some progress has been made toward understanding these transformations in neural networks with smooth activation functions, an understanding in the more general setting of non-smooth activation functions, such as the rectified linear unit (ReLU), which tend to perform better, is required. Here we propose that the geometric transformations performed by DNNs during classification tasks have parallels to those expected under Hamilton's Ricci flow - a tool from differential geometry that evolves a manifold by smoothing its curvature, in order to identify its topology. To illustrate this idea, we present a computational framework to quantify the geometric changes that occur as data passes through successive layers of a DNN, and use this framework to motivate a notion of `global Ricci network flow' that can be used to assess a DNN's ability to disentangle complex data geometries to solve classification problems. By training more than $1,500$ DNN classifiers of different widths and depths on synthetic and real-world data, we show that the strength of global Ricci network flow-like behaviour correlates with accuracy for well-trained DNNs, independently of depth, width and data set. Our findings motivate the use of tools from differential and discrete geometry to the problem of explainability in deep learning.
△ Less
Submitted 22 April, 2024;
originally announced April 2024.
-
BeyondPlanck IV. On end-to-end simulations in CMB analysis -- Bayesian versus frequentist statistics
Authors:
M. Brilenkov,
K. S. F. Fornazier,
L. T. Hergt,
G. A. Hoerning,
A. Marins,
T. Murokoshi,
F. Rahman,
N. -O. Stutzer,
Y. Zhou,
F. B. Abdalla,
K. J. Andersen,
R. Aurlien,
R. Banerji,
A. Basyrov,
A. Battista,
M. Bersanelli,
S. Bertocco,
S. Bollanos,
L. P. L. Colombo,
H. K. Eriksen,
J. R. Eskilt,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta
, et al. (26 additional authors not shown)
Abstract:
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…
▽ More
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.
△ Less
Submitted 9 September, 2022;
originally announced September 2022.
-
BeyondPlanck X. Planck LFI frequency maps with sample-based error propagation
Authors:
A. Basyrov,
A. -S. Suur-Uski,
L. P. L. Colombo,
J. R. Eskilt,
S. Paradiso,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle
, et al. (15 additional authors not shown)
Abstract:
We present Planck LFI frequency sky maps derived within the BeyondPlanck framework. This framework draws samples from a global posterior distribution that includes instrumental, astrophysical and cosmological parameters, and the main product is an entire ensemble of frequency sky map samples. This ensemble allows for computationally convenient end-to-end propagation of low-level instrumental uncer…
▽ More
We present Planck LFI frequency sky maps derived within the BeyondPlanck framework. This framework draws samples from a global posterior distribution that includes instrumental, astrophysical and cosmological parameters, and the main product is an entire ensemble of frequency sky map samples. This ensemble allows for computationally convenient end-to-end propagation of low-level instrumental uncertainties into higher-level science products. We show that the two dominant sources of LFI instrumental systematic uncertainties are correlated noise and gain fluctuations, and the products presented here support - for the first time - full Bayesian error propagation for these effects at full angular resolution. We compare our posterior mean maps with traditional frequency maps delivered by the Planck collaboration, and find generally good agreement. The most important quality improvement is due to significantly lower calibration uncertainties in the new processing, as we find a fractional absolute calibration uncertainty at 70 GHz of $δg_{0}/g_{0} =5 \cdot 10^{-5}$, which is nominally 40 times smaller than that reported by Planck 2018. However, the original Planck 2018 estimate has a non-trivial statistical interpretation, and this further illustrates the advantage of the new framework in terms of producing self-consistent and well-defined error estimates of all involved quantities without the need of ad hoc uncertainty contributions. We describe how low-resolution data products, including dense pixel-pixel covariance matrices, may be produced directly from the posterior samples without the need for computationally expensive analytic calculations or simulations. We conclude that posterior-based frequency map sampling provides unique capabilities in terms of low-level systematics modelling and error propagation, and may play an important role for future CMB B-mode experiments. (Abridged.)
△ Less
Submitted 30 August, 2022;
originally announced August 2022.
-
BeyondPlanck XI. Bayesian CMB analysis with sample-based end-to-end error propagation
Authors:
L. P. L. Colombo,
J. R. Eskilt,
S. Paradiso,
H. Thommesen,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell
, et al. (14 additional authors not shown)
Abstract:
We present posterior sample-based cosmic microwave background (CMB) constraints from Planck LFI and WMAP observations derived through global end-to-end Bayesian processing. We use these samples to study correlations between CMB, foreground, and instrumental parameters, and we identify a particularly strong degeneracy between CMB temperature fluctuations and free-free emission on intermediate angul…
▽ More
We present posterior sample-based cosmic microwave background (CMB) constraints from Planck LFI and WMAP observations derived through global end-to-end Bayesian processing. We use these samples to study correlations between CMB, foreground, and instrumental parameters, and we identify a particularly strong degeneracy between CMB temperature fluctuations and free-free emission on intermediate angular scales, which is mitigated through model reduction, masking, and resampling. We compare our posterior-based CMB results with previous Planck products, and find generally good agreement, but with higher noise due to exclusion of HFI data. We find a best-fit CMB dipole amplitude of $3362.7\pm1.4μK$, in excellent agreement with previous Planck results. The quoted uncertainty is derived directly from the sampled posterior distribution, and does not involve any ad hoc contribution for systematic effects. Similarly, we find a temperature quadrupole amplitude of $σ^{TT}_2=229\pm97μK^2$, in good agreement with previous results in terms of the amplitude, but the uncertainty is an order of magnitude larger than the diagonal Fisher uncertainty. Relatedly, we find lower evidence for a possible alignment between $\ell = 2$ and $\ell = 3$ than previously reported due to a much larger scatter in the individual quadrupole coefficients, caused both by marginalizing over a more complete set of systematic effects, and by our more conservative analysis mask. For higher multipoles, we find that the angular temperature power spectrum is generally in good agreement with both Planck and WMAP. This is the first time the sample-based asymptotically exact Blackwell-Rao estimator has been successfully established for multipoles up to $\ell\le600$, and it now accounts for the majority of the cosmologically important information. Cosmological parameter constraints are presented in a companion paper. (Abriged)
△ Less
Submitted 30 August, 2022;
originally announced August 2022.
-
From BeyondPlanck to Cosmoglobe: Open Science, Reproducibility, and Data Longevity
Authors:
S. Gerakakis,
M. Brilenkov,
M. Ieronymaki,
M. San,
D. J. Watts,
K. J. Andersen,
R. Aurlien,
R. Banerji,
A. Basyrov,
M. Bersanelli,
S. Bertocco,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
J. R. Eskilt,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
H. T. Ihle
, et al. (17 additional authors not shown)
Abstract:
The BeyondPlanck and Cosmoglobe collaborations have implemented the first integrated Bayesian end-to-end analysis pipeline for CMB experiments. The primary long-term motivation for this work is to develop a common analysis platform that supports efficient global joint analysis of complementary radio, microwave, and sub-millimeter experiments. A strict prerequisite for this to succeed is broad part…
▽ More
The BeyondPlanck and Cosmoglobe collaborations have implemented the first integrated Bayesian end-to-end analysis pipeline for CMB experiments. The primary long-term motivation for this work is to develop a common analysis platform that supports efficient global joint analysis of complementary radio, microwave, and sub-millimeter experiments. A strict prerequisite for this to succeed is broad participation from the CMB community, and two foundational aspects of the program are therefore reproducibility and Open Science. In this paper, we discuss our efforts toward this aim. We also discuss measures toward facilitating easy code and data distribution, community-based code documentation, user-friendly compilation procedures, etc. This work represents the first publicly released end-to-end CMB analysis pipeline that includes raw data, source code, parameter files, and documentation. We argue that such a complete pipeline release should be a requirement for all major future and publicly-funded CMB experiments, noting that a full public release significantly increases data longevity by ensuring that the data quality can be improved whenever better processing techniques, complementary datasets, or more computing power become available, and thereby also taxpayers' value for money; providing only raw data and final products is not sufficient to guarantee full reproducibility in the future.
△ Less
Submitted 14 March, 2023; v1 submitted 20 May, 2022;
originally announced May 2022.
-
BeyondPlanck XII. Cosmological parameter constraints with end-to-end error propagation
Authors:
S. Paradiso,
L. P. L. Colombo,
K. J. Andersen,
R. Aurlien,
R. Banerji,
A. Basyrov,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
H. K. Eriksen,
J. R. Eskilt,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell
, et al. (16 additional authors not shown)
Abstract:
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…
▽ More
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.
△ Less
Submitted 20 May, 2022;
originally announced May 2022.
-
From BeyondPlanck to Cosmoglobe: Preliminary $\mathit{WMAP}$ $\mathit Q$-band analysis
Authors:
D. J. Watts,
M. Galloway,
H. T. Ihle,
K. J. Andersen,
R. Aurlien,
R. Banerji,
A. Basyrov,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
J. R. Eskilt,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
J. B. Jewell
, et al. (18 additional authors not shown)
Abstract:
We present the first application of the Cosmoglobe analysis framework by analyzing 9-year $\mathit{WMAP}$ time-ordered observations using similar machinery as BeyondPlanck utilizes for $\mathit{Planck}$ LFI. We analyze only the $\mathit Q$-band (41 GHz) data and report on the low-level analysis process from uncalibrated time-ordered data to calibrated maps. Most of the existing BeyondPlanck pipeli…
▽ More
We present the first application of the Cosmoglobe analysis framework by analyzing 9-year $\mathit{WMAP}$ time-ordered observations using similar machinery as BeyondPlanck utilizes for $\mathit{Planck}$ LFI. We analyze only the $\mathit Q$-band (41 GHz) data and report on the low-level analysis process from uncalibrated time-ordered data to calibrated maps. Most of the existing BeyondPlanck pipeline may be reused for $\mathit{WMAP}$ analysis with minimal changes to the existing codebase. The main modification is the implementation of the same preconditioned biconjugate gradient mapmaker used by the $\mathit{WMAP}$ team. Producing a single $\mathit{WMAP}$ $\mathit Q$1-band sample requires 22 CPU-hrs, which is slightly more than the cost of a $\mathit{Planck}$ 44 GHz sample of 17 CPU-hrs; this demonstrates that full end-to-end Bayesian processing of the $\mathit{WMAP}$ data is computationally feasible. In general, our recovered maps are very similar to the maps released by the $\mathit{WMAP}$ team, although with two notable differences. In temperature we find a $\sim2\,\mathrm{μK}$ quadrupole difference that most likely is caused by different gain modeling, while in polarization we find a distinct $2.5\,\mathrm{μK}$ signal that has been previously called poorly-measured modes by the $\mathit{WMAP}$ team. In the Cosmoglobe processing, this pattern arises from temperature-to-polarization leakage from the coupling between the CMB Solar dipole, transmission imbalance, and sidelobes. No traces of this pattern are found in either the frequency map or TOD residual map, suggesting that the current processing has succeeded in modelling these poorly measured modes within the assumed parametric model by using $\mathit{Planck}$ information to break the sky-synchronous degeneracies inherent in the $\mathit{WMAP}$ scanning strategy.
△ Less
Submitted 23 May, 2022; v1 submitted 24 February, 2022;
originally announced February 2022.
-
Probing Cosmic Inflation with the LiteBIRD Cosmic Microwave Background Polarization Survey
Authors:
LiteBIRD Collaboration,
E. Allys,
K. Arnold,
J. Aumont,
R. Aurlien,
S. Azzoni,
C. Baccigalupi,
A. J. Banday,
R. Banerji,
R. B. Barreiro,
N. Bartolo,
L. Bautista,
D. Beck,
S. Beckman,
M. Bersanelli,
F. Boulanger,
M. Brilenkov,
M. Bucher,
E. Calabrese,
P. Campeti,
A. Carones,
F. J. Casas,
A. Catalano,
V. Chan,
K. Cheung
, et al. (166 additional authors not shown)
Abstract:
LiteBIRD, the Lite (Light) satellite for the study of B-mode polarization and Inflation from cosmic background Radiation Detection, is a space mission for primordial cosmology and fundamental physics. The Japan Aerospace Exploration Agency (JAXA) selected LiteBIRD in May 2019 as a strategic large-class (L-class) mission, with an expected launch in the late 2020s using JAXA's H3 rocket. LiteBIRD is…
▽ More
LiteBIRD, the Lite (Light) satellite for the study of B-mode polarization and Inflation from cosmic background Radiation Detection, is a space mission for primordial cosmology and fundamental physics. The Japan Aerospace Exploration Agency (JAXA) selected LiteBIRD in May 2019 as a strategic large-class (L-class) mission, with an expected launch in the late 2020s using JAXA's H3 rocket. LiteBIRD is planned to orbit the Sun-Earth Lagrangian point L2, where it will map the cosmic microwave background (CMB) polarization over the entire sky for three years, with three telescopes in 15 frequency bands between 34 and 448 GHz, to achieve an unprecedented total sensitivity of 2.2$μ$K-arcmin, with a typical angular resolution of 0.5$^\circ$ at 100 GHz. The primary scientific objective of LiteBIRD is to search for the signal from cosmic inflation, either making a discovery or ruling out well-motivated inflationary models. The measurements of LiteBIRD will also provide us with insight into the quantum nature of gravity and other new physics beyond the standard models of particle physics and cosmology. We provide an overview of the LiteBIRD project, including scientific objectives, mission and system requirements, operation concept, spacecraft and payload module design, expected scientific outcomes, potential design extensions and synergies with other projects.
△ Less
Submitted 27 March, 2023; v1 submitted 6 February, 2022;
originally announced February 2022.
-
Predictive Closed-Loop Service Automation in O-RAN based Network Slicing
Authors:
Joseph Thaliath,
Solmaz Niknam,
Sukhdeep Singh,
Rahul Banerji,
Navrati Saxena,
Harpreet S. Dhillon,
Jeffrey H. Reed,
Ali Kashif Bashir,
Avinash Bhat,
Abhishek Roy
Abstract:
Network slicing provides introduces customized and agile network deployment for managing different service types for various verticals under the same infrastructure. To cater to the dynamic service requirements of these verticals and meet the required quality-of-service (QoS) mentioned in the service-level agreement (SLA), network slices need to be isolated through dedicated elements and resources…
▽ More
Network slicing provides introduces customized and agile network deployment for managing different service types for various verticals under the same infrastructure. To cater to the dynamic service requirements of these verticals and meet the required quality-of-service (QoS) mentioned in the service-level agreement (SLA), network slices need to be isolated through dedicated elements and resources. Additionally, allocated resources to these slices need to be continuously monitored and intelligently managed. This enables immediate detection and correction of any SLA violation to support automated service assurance in a closed-loop fashion. By reducing human intervention, intelligent and closed-loop resource management reduces the cost of offering flexible services. Resource management in a network shared among verticals (potentially administered by different providers), would be further facilitated through open and standardized interfaces. Open radio access network (O-RAN) is perhaps the most promising RAN architecture that inherits all the aforementioned features, namely intelligence, open and standard interfaces, and closed control loop. Inspired by this, in this article we provide a closed-loop and intelligent resource provisioning scheme for O-RAN slicing to prevent SLA violations. In order to maintain realism, a real-world dataset of a large operator is used to train a learning solution for optimizing resource utilization in the proposed closed-loop service automation process. Moreover, the deployment architecture and the corresponding flow that are cognizant of the O-RAN requirements are also discussed.
△ Less
Submitted 3 February, 2022;
originally announced February 2022.
-
BeyondPlanck XIII. Intensity foreground sampling, degeneracies, and priors
Authors:
K. J. Andersen,
D. Herman,
R. Aurlien,
R. Banerji,
A. Basyrov,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
J. R. Eskilt,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci
, et al. (19 additional authors not shown)
Abstract:
We present the intensity foreground algorithms and model employed within the BeyondPlanck analysis framework. The BeyondPlanck analysis is aimed at integrating component separation and instrumental parameter sampling within a global framework, leading to complete end-to-end error propagation in the $Planck$ Low Frequency Instrument (LFI) data analysis. Given the scope of the BeyondPlanck analysis,…
▽ More
We present the intensity foreground algorithms and model employed within the BeyondPlanck analysis framework. The BeyondPlanck analysis is aimed at integrating component separation and instrumental parameter sampling within a global framework, leading to complete end-to-end error propagation in the $Planck$ Low Frequency Instrument (LFI) data analysis. Given the scope of the BeyondPlanck analysis, a limited set of data is included in the component separation process, leading to foreground parameter degeneracies. In order to properly constrain the Galactic foreground parameters, we improve upon the previous $\texttt{Commander}$ component separation implementation by adding a suite of algorithmic techniques. These algorithms are designed to improve the stability and computational efficiency for weakly constrained posterior distributions. These are: 1) joint foreground spectral parameter and amplitude sampling, building on ideas from Miramare; 2) component-based monopole determination; 3) joint spectral parameter and monopole sampling; and 4) application of informative spatial priors for component amplitude maps. We find that the only spectral parameter with a significant signal-to-noise ratio using the current BeyondPlanck data set is the peak frequency of the anomalous microwave emission component, for which we find $ν_{\mathrm{p}}=25.3\pm0.5$ GHz; all others must be constrained through external priors. Future works will be aimed at integrating many more data sets into this analysis, both map and time-ordered based, thereby gradually eliminating the currently observed degeneracies in a controlled manner with respect to both instrumental systematic effects and astrophysical degeneracies. When this happens, the simple LFI-oriented data model employed in the current work will need to be generalized to account for both a richer astrophysical model and additional instrumental effects.
△ Less
Submitted 1 October, 2022; v1 submitted 20 January, 2022;
originally announced January 2022.
-
BeyondPlanck XVI. Limits on Large-Scale Polarized Anomalous Microwave Emission from Planck LFI and WMAP
Authors:
D. Herman,
B. Hensley,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci,
E. Keihänen,
R. Keskitalo
, et al. (13 additional authors not shown)
Abstract:
We constrain the level of polarized anomalous microwave emission (AME) on large angular scales using $\textit{Planck}$ LFI and $\textit{WMAP}$ polarization data within a Bayesian CMB analysis framework. We model synchrotron emission with a power-law spectral energy distribution, and the sum of AME and thermal dust emission through linear regression with the $\textit{Planck}$ HFI 353 GHz data. This…
▽ More
We constrain the level of polarized anomalous microwave emission (AME) on large angular scales using $\textit{Planck}$ LFI and $\textit{WMAP}$ polarization data within a Bayesian CMB analysis framework. We model synchrotron emission with a power-law spectral energy distribution, and the sum of AME and thermal dust emission through linear regression with the $\textit{Planck}$ HFI 353 GHz data. This template-based dust emission model allows us to constrain the level of polarized AME while making minimal assumptions on its frequency dependence. We neglect cosmic microwave background fluctuations, but show through simulations that these have a minor impact on the results. We find that the resulting AME polarization fraction confidence limit is sensitive to the polarized synchrotron spectral index prior, and for priors steeper than $β_{\mathrm{s}} = -3.1\pm0.1$ we find an upper limit of $p_{\mathrm{AME}}^{\rm max}\lesssim 0.6\,\%$ ($95\,\%$ confidence). In contrast, for $β_{\mathrm{s}}=-3.0\pm0.1$, we find a nominal detection of $p_{\mathrm{AME}}=2.5\pm1.0\,\%$ ($95\,\%$ confidence). These data are thus not strong enough to simultaneously and robustly constrain both polarized synchrotron emission and AME, and our main result is therefore a constraint on the AME polarization fraction explicitly as a function of $β_\mathrm{s}$. Combining the current $\textit{Planck}$ and $\textit{WMAP}$ observations with measurements from high-sensitivity low-frequency experiments such as C-BASS and QUIJOTE will be critical to improve these limits further.
△ Less
Submitted 26 April, 2022; v1 submitted 10 January, 2022;
originally announced January 2022.
-
BeyondPlanck III. Commander3
Authors:
M. Galloway,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci,
E. Keihänen,
R. Keskitalo
, et al. (13 additional authors not shown)
Abstract:
We describe the computational infrastructure for end-to-end Bayesian CMB analysis implemented by the BeyondPlanck collaboration. This code is called commander3, and provides a statistically consistent framework for global analysis of CMB and microwave observations, and may be useful for a wide range of legacy, current, and future experiments. The paper has three main goals. Firstly, we provide a h…
▽ More
We describe the computational infrastructure for end-to-end Bayesian CMB analysis implemented by the BeyondPlanck collaboration. This code is called commander3, and provides a statistically consistent framework for global analysis of CMB and microwave observations, and may be useful for a wide range of legacy, current, and future experiments. The paper has three main goals. Firstly, we provide a high-level overview of the existing code base, aiming to guide readers who wish to extend and adapt the code according to their own needs, or to reimplement it from scratch in a different programming language. Secondly, we discuss some critical computational challenges that arise within any global CMB analysis framework, for instance in-memory compression of time-ordered data, FFT optimization, and parallelization and load-balancing. Thirdly, we quantify the CPU and RAM requirements for the current BeyondPlanck analysis, and find that a total of 1.5 TB of RAM is required for efficient analysis, and the total cost of a full Gibbs sample is 170 CPU-hrs, including both low-level processing and high-level component separation, which is well within the capabilities of current low-cost computing facilities. The existing code base is made publicly available under a GNU General Public Library (GPL) license.
△ Less
Submitted 10 January, 2022;
originally announced January 2022.
-
BeyondPlanck VIII. Efficient Sidelobe Convolution and Correction through Spin Harmonics
Authors:
M. Galloway,
M. Reinecke,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci,
E. Keihänen
, et al. (13 additional authors not shown)
Abstract:
We introduce a new formulation of the Conviqt convolution algorithm in terms of spin harmonics, and apply this to the problem of sidelobe correction for BeyondPlanck, the first end-to-end Bayesian Gibbs sampling framework for CMB analysis. We compare our implementation to the previous Planck LevelS implementation, and find good agreement between the two codes in terms of accuracy, but with a speed…
▽ More
We introduce a new formulation of the Conviqt convolution algorithm in terms of spin harmonics, and apply this to the problem of sidelobe correction for BeyondPlanck, the first end-to-end Bayesian Gibbs sampling framework for CMB analysis. We compare our implementation to the previous Planck LevelS implementation, and find good agreement between the two codes in terms of accuracy, but with a speed-up reaching a factor of 3--10, depending on the frequency bandlimits, $l_{\textrm{max}}$ and $m_{\textrm{max}}$. The new algorithm is significantly simpler to implement and maintain, since all low-level calculations are handled through an external spherical harmonic transform library. We find that our mean sidelobe estimates for Planck LFI agree well with previous efforts. Additionally, we present novel sidelobe rms maps that quantify the uncertainty in the sidelobe corrections due to variations in the sky model.
△ Less
Submitted 10 January, 2022;
originally announced January 2022.
-
BeyondPlanck X. Bandpass and beam leakage corrections
Authors:
T. L. Svalheim,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci,
E. Keihänen
, et al. (14 additional authors not shown)
Abstract:
We discuss the treatment of bandpass and beam leakage corrections in the Bayesian BeyondPlanck CMB analysis pipeline as applied to the Planck LFI measurements. As a preparatory step, we first apply three corrections to the nominal LFI bandpass profiles including removal of a known systematic effect in the ground measuring equipment at 61 GHz; smoothing of standing wave ripples; and edge regulariza…
▽ More
We discuss the treatment of bandpass and beam leakage corrections in the Bayesian BeyondPlanck CMB analysis pipeline as applied to the Planck LFI measurements. As a preparatory step, we first apply three corrections to the nominal LFI bandpass profiles including removal of a known systematic effect in the ground measuring equipment at 61 GHz; smoothing of standing wave ripples; and edge regularization. The main net impact of these modifications is an overall shift in the 70 GHz bandpass of +0.6 GHz; we argue that any analysis of LFI data products, either from Planck or BeyondPlanck, should use these new bandpasses. In addition, we fit a single free bandpass parameter for each radiometer of the form $Δ_i = Δ_0 + δ_i$, where $Δ_0$ represents an absolute frequency shift per frequency band and $δ_i$ is a relative shift per detector. The absolute correction is only fitted at 30 GHz with a full $χ^2$-based likelihood, resulting in a correction of $Δ_{30}=0.24\pm0.03\,$GHz. The relative corrections are fitted using a spurious map approach, fundamentally similar to the method pioneered by the WMAP team, but without introducing many additional degrees of freedom. All bandpass parameters are sampled using a standard Metropolis sampler within the main BeyondPlanck Gibbs chain, and bandpass uncertainties are thus propagated to all other data products in the analysis. In total, we find that our bandpass model significantly reduces leakage effects. For beam leakage corrections, we adopt the official Planck LFI beam estimates without additional degrees of freedom, and only marginalize over the underlying sky model. We note that this is the first time leakage from beam mismatch has been included for Planck LFI maps.
△ Less
Submitted 10 January, 2022;
originally announced January 2022.
-
Practical Adoption of Cloud Computing in Power Systems- Drivers, Challenges, Guidance, and Real-world Use Cases
Authors:
Song Zhang,
Amritanshu Pandey,
Xiaochuan Luo,
Maggy Powell,
Ranjan Banerji,
Lei Fan,
Abhineet Parchure,
Edgardo Luzcando
Abstract:
Motivated by The Federal Energy Regulatory Commission's (FERC) recent direction and ever-growing interest in cloud adoption by power utilities, a Task Force was established to assist power system practitioners with secure, reliable and cost-effective adoption of cloud technology to meet various business needs. This paper summarizes the business drivers, challenges, guidance, and best practices for…
▽ More
Motivated by The Federal Energy Regulatory Commission's (FERC) recent direction and ever-growing interest in cloud adoption by power utilities, a Task Force was established to assist power system practitioners with secure, reliable and cost-effective adoption of cloud technology to meet various business needs. This paper summarizes the business drivers, challenges, guidance, and best practices for cloud adoption in power systems from the Task Force's perspective, after extensive review and deliberation by its members, including grid operators, utility companies, software vendors, and cloud providers. The paper begins by enumerating various business drivers for cloud adoption in the power industry. It follows with the discussion of the challenges and risks of migrating power grid utility workloads to the cloud. Next, for each corresponding challenge or risk, the paper provides appropriate guidance. Notably, the guidance is directed toward power industry professionals who are considering cloud solutions and are yet hesitant about the practical execution. Finally, to tie all the sections together, the paper documents various real-world use cases of cloud technology in the power system domain, which both the power industry practitioners and software vendors can look toward to design and select their own future cloud solutions. We hope that the information in this paper will serve as helpful guidance for the development of NERC guidelines and standards relevant to cloud adoption in the industry.
△ Less
Submitted 2 February, 2022; v1 submitted 31 July, 2021;
originally announced August 2021.
-
Overview of the Medium and High Frequency Telescopes of the LiteBIRD satellite mission
Authors:
L. Montier,
B. Mot,
P. de Bernardis,
B. Maffei,
G. Pisano,
F. Columbro,
J. E. Gudmundsson,
S. Henrot-Versillé,
L. Lamagna,
J. Montgomery,
T. Prouvé,
M. Russell,
G. Savini,
S. Stever,
K. L. Thompson,
M. Tsujimoto,
C. Tucker,
B. Westbrook,
P. A. R. Ade,
A. Adler,
E. Allys,
K. Arnold,
D. Auguste,
J. Aumont,
R. Aurlien
, et al. (212 additional authors not shown)
Abstract:
LiteBIRD is a JAXA-led Strategic Large-Class mission designed to search for the existence of the primordial gravitational waves produced during the inflationary phase of the Universe, through the measurements of their imprint onto the polarization of the cosmic microwave background (CMB). These measurements, requiring unprecedented sensitivity, will be performed over the full sky, at large angular…
▽ More
LiteBIRD is a JAXA-led Strategic Large-Class mission designed to search for the existence of the primordial gravitational waves produced during the inflationary phase of the Universe, through the measurements of their imprint onto the polarization of the cosmic microwave background (CMB). These measurements, requiring unprecedented sensitivity, will be performed over the full sky, at large angular scales, and over 15 frequency bands from 34GHz to 448GHz. The LiteBIRD instruments consist of three telescopes, namely the Low-, Medium- and High-Frequency Telescope (respectively LFT, MFT and HFT). We present in this paper an overview of the design of the Medium-Frequency Telescope (89-224GHz) and the High-Frequency Telescope (166-448GHz), the so-called MHFT, under European responsibility, which are two cryogenic refractive telescopes cooled down to 5K. They include a continuous rotating half-wave plate as the first optical element, two high-density polyethylene (HDPE) lenses and more than three thousand transition-edge sensor (TES) detectors cooled to 100mK. We provide an overview of the concept design and the remaining specific challenges that we have to face in order to achieve the scientific goals of LiteBIRD.
△ Less
Submitted 1 February, 2021;
originally announced February 2021.
-
Concept Design of Low Frequency Telescope for CMB B-mode Polarization satellite LiteBIRD
Authors:
Y. Sekimoto,
P. A. R. Ade,
A. Adler,
E. Allys,
K. Arnold,
D. Auguste,
J. Aumont,
R. Aurlien,
J. Austermann,
C. Baccigalupi,
A. J. Banday,
R. Banerji,
R. B. Barreiro,
S. Basak,
J. Beall,
D. Beck,
S. Beckman,
J. Bermejo,
P. de Bernardis,
M. Bersanelli,
J. Bonis,
J. Borrill,
F. Boulanger,
S. Bounissou,
M. Brilenkov
, et al. (212 additional authors not shown)
Abstract:
LiteBIRD has been selected as JAXA's strategic large mission in the 2020s, to observe the cosmic microwave background (CMB) $B$-mode polarization over the full sky at large angular scales. The challenges of LiteBIRD are the wide field-of-view (FoV) and broadband capabilities of millimeter-wave polarization measurements, which are derived from the system requirements. The possible paths of stray li…
▽ More
LiteBIRD has been selected as JAXA's strategic large mission in the 2020s, to observe the cosmic microwave background (CMB) $B$-mode polarization over the full sky at large angular scales. The challenges of LiteBIRD are the wide field-of-view (FoV) and broadband capabilities of millimeter-wave polarization measurements, which are derived from the system requirements. The possible paths of stray light increase with a wider FoV and the far sidelobe knowledge of $-56$ dB is a challenging optical requirement. A crossed-Dragone configuration was chosen for the low frequency telescope (LFT : 34--161 GHz), one of LiteBIRD's onboard telescopes. It has a wide field-of-view ($18^\circ \times 9^\circ$) with an aperture of 400 mm in diameter, corresponding to an angular resolution of about 30 arcminutes around 100 GHz. The focal ratio f/3.0 and the crossing angle of the optical axes of 90$^\circ$ are chosen after an extensive study of the stray light. The primary and secondary reflectors have rectangular shapes with serrations to reduce the diffraction pattern from the edges of the mirrors. The reflectors and structure are made of aluminum to proportionally contract from warm down to the operating temperature at $5\,$K. A 1/4 scaled model of the LFT has been developed to validate the wide field-of-view design and to demonstrate the reduced far sidelobes. A polarization modulation unit (PMU), realized with a half-wave plate (HWP) is placed in front of the aperture stop, the entrance pupil of this system. A large focal plane with approximately 1000 AlMn TES detectors and frequency multiplexing SQUID amplifiers is cooled to 100 mK. The lens and sinuous antennas have broadband capability. Performance specifications of the LFT and an outline of the proposed verification plan are presented.
△ Less
Submitted 15 January, 2021;
originally announced January 2021.
-
BeyondPlanck XV. Polarized foreground emission between 30 and 70 GHz
Authors:
T. L. Svalheim,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci,
E. Keihänen
, et al. (13 additional authors not shown)
Abstract:
We constrain polarized foreground emission between 30 and 70 GHz with the Planck Low Frequency Instrument (LFI) and WMAP data within the global Bayesian BeyondPlanck framework. We combine for the first time full-resolution Planck LFI time-ordered data with low-resolution WMAP sky maps at 33, 40 and 61 GHz. Spectral parameters are fit with a likelihood defined at the native resolution of each frequ…
▽ More
We constrain polarized foreground emission between 30 and 70 GHz with the Planck Low Frequency Instrument (LFI) and WMAP data within the global Bayesian BeyondPlanck framework. We combine for the first time full-resolution Planck LFI time-ordered data with low-resolution WMAP sky maps at 33, 40 and 61 GHz. 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 current data 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 signal-dominated, and for these we derive spectral indices of $β_{\mathrm s}^{\mathrm{Spur}}=-3.17\pm0.06$ and $β_{\mathrm s}^{\mathrm{Plane}}=-3.03\pm0.07$, in good agreement with previous results. For thermal dust emission we assume a modified blackbody model and we fit a single power-law index across the full sky. We find $β_{\mathrm{d}}=1.64\pm0.03$, which is slightly steeper than reported from Planck HFI data, but still statistically consistent at the 2$σ$ confidence level.
△ Less
Submitted 18 January, 2022; v1 submitted 17 November, 2020;
originally announced November 2020.
-
BeyondPlanck VII. Bayesian estimation of gain and absolute calibration for CMB experiments
Authors:
E. Gjerløw,
H. T. Ihle,
S. Galeotta,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
M. Galloway,
S. Gerakakis,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
J. B. Jewell,
A. Karakci,
E. Keihänen,
R. Keskitalo
, et al. (13 additional authors not shown)
Abstract:
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-ind…
▽ More
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.67% (-0.84%) for 30 GHz, 0.12% (-0.04%) for 44 GHz and -0.03% (-0.64%) for 70 GHz, compared to Planck DR4 (Planck 2018). 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.
△ Less
Submitted 20 May, 2022; v1 submitted 16 November, 2020;
originally announced November 2020.
-
BeyondPlanck VI. Noise characterization and modelling
Authors:
H. T. Ihle,
M. Bersanelli,
C. Franceschet,
E. Gjerløw,
K. J. Andersen,
R. Aurlien,
R. Banerji,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
J. R. Eskilt,
M. K. Foss,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
J. B. Jewell,
A. Karakci,
E. Keihänen
, et al. (16 additional authors not shown)
Abstract:
We present a Bayesian method for estimating instrumental noise parameters and propagating noise uncertainties within the global BeyondPlanck Gibbs sampling framework, and apply this to Planck Low Frequency Instrument (LFI) time-ordered data. Following previous literature, we initially adopt a $1/f$ model for the noise power spectral density (PSD), but find the need for an additional lognormal comp…
▽ More
We present a Bayesian method for estimating instrumental noise parameters and propagating noise uncertainties within the global BeyondPlanck Gibbs sampling framework, and apply this to Planck Low Frequency Instrument (LFI) time-ordered data. Following previous literature, we initially adopt a $1/f$ model for the noise power spectral density (PSD), but find the need for an additional lognormal component in the noise model for the 30 and 44\,GHz bands. We implement an optimal Wiener-filter (or constrained realization) gap-filling procedure to account for masked data. We then use this procedure to both estimate the gapless correlated noise in the time-domain, $n_\mathrm{corr}$, and to sample the noise PSD parameters, $ξ^n = \{σ_0, f_\mathrm{knee}, α, A_\mathrm{p}\}$. In contrast to previous \textit{Planck} analyses, we assume piecewise stationary noise only within each pointing period (PID), not throughout the full mission, but we adopt the LFI Data Processing Center (DPC) results as priors on $α$ and $f_\mathrm{knee}$. On average, we find best-fit correlated noise parameters that are mostly consistent with previous results, with a few notable exceptions. However, a detailed inspection of the time-dependent results reveals many important findings. First and foremost, we find strong evidence for statistically significant temporal variations in all noise PSD parameters, many of which are directly correlated with satellite housekeeping data. Second, while the simple $1/f$ model appears to be an excellent fit for the LFI 70 GHz channel, there is evidence for additional correlated noise not described by a $1/f$ model in the 30 and 44 GHz channels, including within the primary science frequency range of 0.1--1 Hz. (Abridged)
△ Less
Submitted 19 July, 2022; v1 submitted 12 November, 2020;
originally announced November 2020.
-
BeyondPlanck II. CMB map-making through Gibbs sampling
Authors:
E. Keihänen,
A. -S. Suur-Uski,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci
, et al. (15 additional authors not shown)
Abstract:
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 benef…
▽ More
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.
△ Less
Submitted 11 November, 2020;
originally announced November 2020.
-
BeyondPlanck I. Global Bayesian analysis of the Planck Low Frequency Instrument data
Authors:
BeyondPlanck Collaboration,
K. J. Andersen,
R. Aurlien,
R. Banerji,
A. Basyrov,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
J. R. Eskilt,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell
, et al. (20 additional authors not shown)
Abstract:
We describe the BeyondPlanck project in terms of motivation, methodology and main products, and provide a guide to a set of companion papers that describe each result in fuller detail. We implement a complete end-to-end Bayesian analysis framework for the Planck LFI observations. The primary product is a full joint posterior distribution $P(ω|d)$, where $ω$ represents the set of all free instrumen…
▽ More
We describe the BeyondPlanck project in terms of motivation, methodology and main products, and provide a guide to a set of companion papers that describe each result in fuller detail. We implement a complete end-to-end Bayesian analysis framework for the Planck LFI observations. The primary product is a full joint posterior distribution $P(ω|d)$, where $ω$ represents the set of all free instrumental, astrophysical, and cosmological parameters. Notable advantages of this approach are seamless end-to-end propagation of uncertainties; accurate modeling of both astrophysical and instrumental effects in the most natural basis for each uncertain quantity; optimized computational costs with little or no need for intermediate human interaction between various analysis steps; and a complete overview of the entire analysis process within one single framework. We focus in particular on low-$\ell$ CMB polarization reconstruction with Planck LFI. We identify several important new effects that have not been accounted for in previous pipelines, including gain over-smoothing and time-variable and non-$1/f$ correlated noise in the 30 and 44 GHz channels. We find that all results are consistent with the $Λ$CDM model, and we constrain the reionization optical depth to $τ=0.066\pm0.013$, with a low-resolution $χ^2$ probability-to-exceed of 32%. This uncertainty is about 30% larger than the official pipelines, arising from taking into account a more complete instrumental model. The marginal CMB Solar dipole amplitude is $3362.7\pm1.4μ\mathrm{K}$, where the error bar is derived directly from the posterior distribution without the need of any ad-hoc instrumental corrections. We are currently not aware of any significant unmodelled systematic effects remaining in the Planck LFI data, and, for the first time, the 44 GHz channel is fully exploited. (Abridged.)
△ Less
Submitted 12 September, 2022; v1 submitted 11 November, 2020;
originally announced November 2020.
-
A Monte Carlo comparison between template-based and Wiener-filter CMB dipole estimators
Authors:
H. Thommesen,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Brilenkov,
H. K. Eriksen,
U. Fuskeland,
M. Galloway,
L. M. Mocanu,
T. L. Svalheim,
I. K. Wehus
Abstract:
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 a…
▽ More
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 3$~μ$K, 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.
△ Less
Submitted 13 July, 2020;
originally announced July 2020.
-
Intelligent O-RAN for Beyond 5G and 6G Wireless Networks
Authors:
Solmaz Niknam,
Abhishek Roy,
Harpreet S. Dhillon,
Sukhdeep Singh,
Rahul Banerji,
Jeffery H. Reed,
Navrati Saxena,
Seungil Yoon
Abstract:
Building on the principles of openness and intelligence, there has been a concerted global effort from the operators towards enhancing the radio access network (RAN) architecture. The objective is to build an operator-defined RAN architecture (and associated interfaces) on open hardware that provides intelligent radio control for beyond fifth generation (5G) as well as future sixth generation (6G)…
▽ More
Building on the principles of openness and intelligence, there has been a concerted global effort from the operators towards enhancing the radio access network (RAN) architecture. The objective is to build an operator-defined RAN architecture (and associated interfaces) on open hardware that provides intelligent radio control for beyond fifth generation (5G) as well as future sixth generation (6G) wireless networks. Specifically, the open-radio access network (O-RAN) alliance has been formed by merging xRAN forum and C-RAN alliance to formally define the requirements that would help achieve this objective. Owing to the importance of O-RAN in the current wireless landscape, this article provides an introduction to the concepts, principles, and requirements of the Open RAN as specified by the O-RAN alliance. In order to illustrate the role of intelligence in O-RAN, we propose an intelligent radio resource management scheme to handle traffic congestion and demonstrate its efficacy on a real-world dataset obtained from a large operator. A high-level architecture of this deployment scenario that is compliant with the O-RAN requirements is also discussed. The article concludes with key technical challenges and open problems for future research and development.
△ Less
Submitted 17 May, 2020;
originally announced May 2020.
-
Updated design of the CMB polarization experiment satellite LiteBIRD
Authors:
H. Sugai,
P. A. R. Ade,
Y. Akiba,
D. Alonso,
K. Arnold,
J. Aumont,
J. Austermann,
C. Baccigalupi,
A. J. Banday,
R. Banerji,
R. B. Barreiro,
S. Basak,
J. Beall,
S. Beckman,
M. Bersanelli,
J. Borrill,
F. Boulanger,
M. L. Brown,
M. Bucher,
A. Buzzelli,
E. Calabrese,
F. J. Casas,
A. Challinor,
V. Chan,
Y. Chinone
, et al. (196 additional authors not shown)
Abstract:
Recent developments of transition-edge sensors (TESs), based on extensive experience in ground-based experiments, have been making the sensor techniques mature enough for their application on future satellite CMB polarization experiments. LiteBIRD is in the most advanced phase among such future satellites, targeting its launch in Japanese Fiscal Year 2027 (2027FY) with JAXA's H3 rocket. It will ac…
▽ More
Recent developments of transition-edge sensors (TESs), based on extensive experience in ground-based experiments, have been making the sensor techniques mature enough for their application on future satellite CMB polarization experiments. LiteBIRD is in the most advanced phase among such future satellites, targeting its launch in Japanese Fiscal Year 2027 (2027FY) with JAXA's H3 rocket. It will accommodate more than 4000 TESs in focal planes of reflective low-frequency and refractive medium-and-high-frequency telescopes in order to detect a signature imprinted on the cosmic microwave background (CMB) by the primordial gravitational waves predicted in cosmic inflation. The total wide frequency coverage between 34GHz and 448GHz enables us to extract such weak spiral polarization patterns through the precise subtraction of our Galaxy's foreground emission by using spectral differences among CMB and foreground signals. Telescopes are cooled down to 5Kelvin for suppressing thermal noise and contain polarization modulators with transmissive half-wave plates at individual apertures for separating sky polarization signals from artificial polarization and for mitigating from instrumental 1/f noise. Passive cooling by using V-grooves supports active cooling with mechanical coolers as well as adiabatic demagnetization refrigerators. Sky observations from the second Sun-Earth Lagrangian point, L2, are planned for three years. An international collaboration between Japan, USA, Canada, and Europe is sharing various roles. In May 2019, the Institute of Space and Astronautical Science (ISAS), JAXA selected LiteBIRD as the strategic large mission No. 2.
△ Less
Submitted 6 January, 2020;
originally announced January 2020.
-
Constraints on the spectral index of polarized synchrotron emission from WMAP and Faraday-corrected S-PASS data
Authors:
U. Fuskeland,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Brilenkov,
H. K. Eriksen,
M. Galloway,
E. Gjerløw,
S. K. Næss,
T. L. Svalheim,
I. K. Wehus
Abstract:
We constrain the spectral index of polarized synchrotron emission, $β_s$, by correlating the recently released 2.3 GHz S-Band Polarization All Sky Survey (S-PASS) data with the 23 GHz 9-year Wilkinson Microwave Anisotropy Probe (WMAP) sky maps. We subdivide the S-PASS field, which covers the southern ecliptic hemisphere, into 95 $15^{\circ}\times15^{\circ}$ regions and estimate the spectral index…
▽ More
We constrain the spectral index of polarized synchrotron emission, $β_s$, by correlating the recently released 2.3 GHz S-Band Polarization All Sky Survey (S-PASS) data with the 23 GHz 9-year Wilkinson Microwave Anisotropy Probe (WMAP) sky maps. We subdivide the S-PASS field, which covers the southern ecliptic hemisphere, into 95 $15^{\circ}\times15^{\circ}$ regions and estimate the spectral index of polarized synchrotron emission within each region using a simple but robust T-T plot technique. Three different versions of the S-PASS data are considered, corresponding to: no correction for Faraday rotation; Faraday correction based on the rotation measure model presented by the S-PASS team; or Faraday correction based on a rotation measure model presented by Hutschenreuter and Enßlin. We find that the correlation between S-PASS and WMAP is strongest when applying the S-PASS model. Adopting this correction model, we find that the mean spectral index of polarized synchrotron emission gradually steepens from $β_s\approx-2.8$ at low Galactic latitudes to $β_s\approx-3.3$ at high Galactic latitudes, in good agreement with previously published results. Finally, we consider two special cases defined by the BICEP2 and SPIDER fields and obtain mean estimates of $β_{BICEP2}=-3.22\pm0.06$ and $β_{SPIDER}=-3.21\pm0.03$, respectively. A comparison with a similar analysis performed in the 23-33 GHz range suggests a flattening of about $Δβ_s \sim 0.1 \pm 0.2$ from low to higher frequencies, but with no statistical significance due to high uncertainties.
△ Less
Submitted 11 March, 2021; v1 submitted 12 September, 2019;
originally announced September 2019.
-
PICO: Probe of Inflation and Cosmic Origins
Authors:
S. Hanany,
M. Alvarez,
E. Artis,
P. Ashton,
J. Aumont,
R. Aurlien,
R. Banerji,
R. B. Barreiro,
J. G. Bartlett,
S. Basak,
N. Battaglia,
J. Bock,
K. K. Boddy,
M. Bonato,
J. Borrill,
F. Bouchet,
F. Boulanger,
B. Burkhart,
J. Chluba,
D. Chuss,
S. Clark,
J. Cooperrider,
B. P. Crill,
G. De Zotti,
J. Delabrouille
, et al. (57 additional authors not shown)
Abstract:
The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.
The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.
△ Less
Submitted 20 August, 2019;
originally announced August 2019.
-
PICO: Probe of Inflation and Cosmic Origins
Authors:
Shaul Hanany,
Marcelo Alvarez,
Emmanuel Artis,
Peter Ashton,
Jonathan Aumont,
Ragnhild Aurlien,
Ranajoy Banerji,
R. Belen Barreiro,
James G. Bartlett,
Soumen Basak,
Nick Battaglia,
Jamie Bock,
Kimberly K. Boddy,
Matteo Bonato,
Julian Borrill,
François Bouchet,
François Boulanger,
Blakesley Burkhart,
Jens Chluba,
David Chuss,
Susan E. Clark,
Joelle Cooperrider,
Brendan P. Crill,
Gianfranco De Zotti,
Jacques Delabrouille
, et al. (57 additional authors not shown)
Abstract:
The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $μ$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a…
▽ More
The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $μ$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $μ$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5σ)$, or will rule out with more than $5σ$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4σ$, constrain the effective number of light particle species with $ΔN_{\rm eff}<0.06~(2σ)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $σ_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.
△ Less
Submitted 5 March, 2019; v1 submitted 26 February, 2019;
originally announced February 2019.
-
Bandpass mismatch error for satellite CMB experiments II: Correcting for the spurious signal
Authors:
Ranajoy Banerji,
Guillaume Patanchon,
Jacques Delabrouille,
Masashi Hazumi,
Duc Thuong Hoang,
Hirokazu Ishino,
Tomotake Matsumura
Abstract:
Future Cosmic Microwave Background (CMB) satellite missions aim at using the B-mode polarisation signal to measure the tensor-to-scalar ratio $r$ with a sensitivity $σ(r)$ of the order of $\leq 10^{-3}$. Small uncertainties in the characterisation of instrument properties such as the spectral filters can lead to a leakage of the intensity signal to polarisation and can possibly bias any measuremen…
▽ More
Future Cosmic Microwave Background (CMB) satellite missions aim at using the B-mode polarisation signal to measure the tensor-to-scalar ratio $r$ with a sensitivity $σ(r)$ of the order of $\leq 10^{-3}$. Small uncertainties in the characterisation of instrument properties such as the spectral filters can lead to a leakage of the intensity signal to polarisation and can possibly bias any measurement of a primordial signal. In this paper we discuss methods for avoiding and correcting for the intensity to polarisation leakage due to bandpass mismatch among detector sets. We develop a template fitting map-maker to obtain an unbiased estimate of the leakage signal and subtract it out of the total signal. Using simulations we show how such a method can reduce the bias on the observed B-mode signal by up to $3$ orders of magnitude in power.
△ Less
Submitted 4 June, 2019; v1 submitted 1 February, 2019;
originally announced February 2019.
-
Exploring cosmic origins with CORE: mitigation of systematic effects
Authors:
P. Natoli,
M. Ashdown,
R. Banerji,
J. Borrill,
A. Buzzelli,
G. de Gasperis,
J. Delabrouille,
E. Hivon,
D. Molinari,
G. Patanchon,
L. Polastri,
M. Tomasi,
F. R. Bouchet,
S. Henrot-Versillé,
D. T. Hoang,
R. Keskitalo,
K. Kiiveri,
T. Kisner,
V. Lindholm,
D. McCarthy,
F. Piacentini,
O. Perdereau,
G. Polenta,
M. Tristram,
A. Achucarro
, et al. (101 additional authors not shown)
Abstract:
We present an analysis of the main systematic effects that could impact the measurement of CMB polarization with the proposed CORE space mission. We employ timeline-to-map simulations to verify that the CORE instrumental set-up and scanning strategy allow us to measure sky polarization to a level of accuracy adequate to the mission science goals. We also show how the CORE observations can be proce…
▽ More
We present an analysis of the main systematic effects that could impact the measurement of CMB polarization with the proposed CORE space mission. We employ timeline-to-map simulations to verify that the CORE instrumental set-up and scanning strategy allow us to measure sky polarization to a level of accuracy adequate to the mission science goals. We also show how the CORE observations can be processed to mitigate the level of contamination by potentially worrying systematics, including intensity-to-polarization leakage due to bandpass mismatch, asymmetric main beams, pointing errors and correlated noise. We use analysis techniques that are well validated on data from current missions such as Planck to demonstrate how the residual contamination of the measurements by these effects can be brought to a level low enough not to hamper the scientific capability of the mission, nor significantly increase the overall error budget. We also present a prototype of the CORE photometric calibration pipeline, based on that used for Planck, and discuss its robustness to systematics, showing how CORE can achieve its calibration requirements. While a fine-grained assessment of the impact of systematics requires a level of knowledge of the system that can only be achieved in a future study phase, the analysis presented here strongly suggests that the main areas of concern for the CORE mission can be addressed using existing knowledge, techniques and algorithms.
△ Less
Submitted 13 July, 2017;
originally announced July 2017.
-
Exploring cosmic origins with CORE: gravitational lensing of the CMB
Authors:
Anthony Challinor,
Rupert Allison,
Julien Carron,
Josquin Errard,
Stephen Feeney,
Thomas Kitching,
Julien Lesgourgues,
Antony Lewis,
Íñigo Zubeldía,
Ana Achucarro,
Peter Ade,
Mark Ashdown,
Mario Ballardini,
A. J. Banday,
Ranajoy Banerji,
James Bartlett,
Nicola Bartolo,
Soumen Basak,
Daniel Baumann,
Marco Bersanelli,
Anna Bonaldi,
Matteo Bonato,
Julian Borrill,
François Bouchet,
François Boulanger
, et al. (88 additional authors not shown)
Abstract:
Lensing of the CMB is now a well-developed probe of large-scale clustering over a broad range of redshifts. By exploiting the non-Gaussian imprints of lensing in the polarization of the CMB, the CORE mission can produce a clean map of the lensing deflections over nearly the full-sky. The number of high-S/N modes in this map will exceed current CMB lensing maps by a factor of 40, and the measuremen…
▽ More
Lensing of the CMB is now a well-developed probe of large-scale clustering over a broad range of redshifts. By exploiting the non-Gaussian imprints of lensing in the polarization of the CMB, the CORE mission can produce a clean map of the lensing deflections over nearly the full-sky. The number of high-S/N modes in this map will exceed current CMB lensing maps by a factor of 40, and the measurement will be sample-variance limited on all scales where linear theory is valid. Here, we summarise this mission product and discuss the science that it will enable. For example, the summed mass of neutrinos will be determined to an accuracy of 17 meV combining CORE lensing and CMB two-point information with contemporaneous BAO measurements, three times smaller than the minimum total mass allowed by neutrino oscillations. In the search for B-mode polarization from primordial gravitational waves with CORE, lens-induced B-modes will dominate over instrument noise, limiting constraints on the gravitational wave power spectrum amplitude. With lensing reconstructed by CORE, one can "delens" the observed polarization internally, reducing the lensing B-mode power by 60%. This improves to 70% by combining lensing and CIB measurements from CORE, reducing the error on the gravitational wave amplitude by 2.5 compared to no delensing (in the null hypothesis). Lensing measurements from CORE will allow calibration of the halo masses of the 40000 galaxy clusters that it will find, with constraints dominated by the clean polarization-based estimators. CORE can accurately remove Galactic emission from CMB maps with its 19 frequency channels. We present initial findings that show that residual Galactic foreground contamination will not be a significant source of bias for lensing power spectrum measurements with CORE. [abridged]
△ Less
Submitted 7 July, 2017;
originally announced July 2017.
-
Bandpass mismatch error for satellite CMB experiments I: Estimating the spurious signal
Authors:
Duc Thuong Hoang,
Guillaume Patanchon,
Martin Bucher,
Tomotake Matsumura,
Ranajoy Banerji,
Hirokazu Ishino,
Masashi Hazumi,
Jacques Delabrouille
Abstract:
Future Cosmic Microwave Background (CMB) satellite missions aim to use the $B$ mode polarization to measure the tensor-to-scalar ratio $r$ with a sensitivity of about $10^{-3}$. Achieving this goal will not only require sufficient detector array sensitivity but also unprecedented control of all systematic errors inherent to CMB polarization measurements. Since polarization measurements derive from…
▽ More
Future Cosmic Microwave Background (CMB) satellite missions aim to use the $B$ mode polarization to measure the tensor-to-scalar ratio $r$ with a sensitivity of about $10^{-3}$. Achieving this goal will not only require sufficient detector array sensitivity but also unprecedented control of all systematic errors inherent to CMB polarization measurements. Since polarization measurements derive from differences between observations at different times and from different sensors, detector response mismatches introduce leakages from intensity to polarization and thus lead to a spurious $B$ mode signal. Because the expected primordial $B$ mode polarization signal is dwarfed by the known unpolarized intensity signal, such leakages could contribute substantially to the final error budget for measuring $r.$ Using simulations we estimate the magnitude and angular spectrum of the spurious $B$ mode signal resulting from bandpass mismatch between different detectors. It is assumed here that the detectors are calibrated, for example using the CMB dipole, so that their sensitivity to the primordial CMB signal has been perfectly matched. Consequently the mismatch in the frequency bandpass shape between detectors introduces difference in the relative calibration of galactic emission components. We simulate using a range of scanning patterns being considered for future satellite missions. We find that the spurious contribution to $r$ from reionization bump on large angular scales ($\ell < 10$) is $\approx 10^{-3}$ assuming large detector arrays and 20 percent of the sky masked. We show how the amplitude of the leakage depends on the angular coverage per pixels that results from the scan pattern.
△ Less
Submitted 20 November, 2017; v1 submitted 28 June, 2017;
originally announced June 2017.
-
Exploring Cosmic Origins with CORE: Survey requirements and mission design
Authors:
J. Delabrouille,
P. de Bernardis,
F. R. Bouchet,
A. Achúcarro,
P. A. R. Ade,
R. Allison,
F. Arroja,
E. Artal,
M. Ashdown,
C. Baccigalupi,
M. Ballardini,
A. J. Banday,
R. Banerji,
D. Barbosa,
J. Bartlett,
N. Bartolo,
S. Basak,
J. J. A. Baselmans,
K. Basu,
E. S. Battistelli,
R. Battye,
D. Baumann,
A. Benoît,
M. Bersanelli,
A. Bideaud
, et al. (178 additional authors not shown)
Abstract:
Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium…
▽ More
Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.
△ Less
Submitted 14 June, 2017;
originally announced June 2017.
-
Exploring Cosmic Origins with CORE: The Instrument
Authors:
P. de Bernardis,
P. A. R. Ade,
J. J. A. Baselmans,
E. S. Battistelli,
A. Benoit,
M. Bersanelli,
A. Bideaud,
M. Calvo,
F. J. Casas,
G. Castellano,
A. Catalano,
I. Charles,
I. Colantoni,
F. Columbro,
A. Coppolecchia,
M. Crook,
G. D'Alessandro,
M. De Petris,
J. Delabrouille,
S. Doyle,
C. Franceschet,
A. Gomez,
J. Goupy,
S. Hanany,
M. Hills
, et al. (104 additional authors not shown)
Abstract:
We describe a space-borne, multi-band, multi-beam polarimeter aiming at a precise and accurate measurement of the polarization of the Cosmic Microwave Background. The instrument is optimized to be compatible with the strict budget requirements of a medium-size space mission within the Cosmic Vision Programme of the European Space Agency. The instrument has no moving parts, and uses arrays of diffr…
▽ More
We describe a space-borne, multi-band, multi-beam polarimeter aiming at a precise and accurate measurement of the polarization of the Cosmic Microwave Background. The instrument is optimized to be compatible with the strict budget requirements of a medium-size space mission within the Cosmic Vision Programme of the European Space Agency. The instrument has no moving parts, and uses arrays of diffraction-limited Kinetic Inductance Detectors to cover the frequency range from 60 GHz to 600 GHz in 19 wide bands, in the focal plane of a 1.2 m aperture telescope cooled at 40 K, allowing for an accurate extraction of the CMB signal from polarized foreground emission. The projected CMB polarization survey sensitivity of this instrument, after foregrounds removal, is 1.7 μK$\cdot$arcmin. The design is robust enough to allow, if needed, a downscoped version of the instrument covering the 100 GHz to 600 GHz range with a 0.8 m aperture telescope cooled at 85 K, with a projected CMB polarization survey sensitivity of 3.2 μK$\cdot$arcmin.
△ Less
Submitted 22 May, 2017; v1 submitted 5 May, 2017;
originally announced May 2017.
-
Exploring cosmic origins with CORE: effects of observer peculiar motion
Authors:
C. Burigana,
C. S. Carvalho,
T. Trombetti,
A. Notari,
M. Quartin,
G. De Gasperis,
A. Buzzelli,
N. Vittorio,
G. De Zotti,
P. de Bernardis,
J. Chluba,
M. Bilicki,
L. Danese,
J. Delabrouille,
L. Toffolatti,
A. Lapi,
M. Negrello,
P. Mazzotta,
D. Scott,
D. Contreras,
A. Achucarro,
P. Ade,
R. Allison,
M. Ashdown,
M. Ballardini
, et al. (94 additional authors not shown)
Abstract:
We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide freque…
▽ More
We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide frequency range, will provide a more accurate estimate of the CMB dipole. The extension of boosting effects to polarization and cross-correlations will enable a more robust determination of purely velocity-driven effects that are not degenerate with the intrinsic CMB dipole, allowing us to achieve a S/N ratio of 13; this improves on the Planck detection and essentially equals that of an ideal cosmic-variance-limited experiment up to a multipole l of 2000. Precise inter-frequency calibration will offer the opportunity to constrain or even detect CMB spectral distortions, particularly from the cosmological reionization, because of the frequency dependence of the dipole spectrum, without resorting to precise absolute calibration. The expected improvement with respect to COBE-FIRAS in the recovery of distortion parameters (in principle, a factor of several hundred for an ideal experiment with the CORE configuration) ranges from a factor of several up to about 50, depending on the quality of foreground removal and relative calibration. Even for 1% accuracy in both foreground removal and relative calibration at an angular scale of 1 deg, we find that dipole analyses for a mission like CORE will be able to improve the recovery of the CIB spectrum amplitude by a factor of 17 in comparison with current results based on FIRAS. In addition to the scientific potential of a mission like CORE for these analyses, synergies with other planned and ongoing projects are also discussed.
△ Less
Submitted 30 August, 2017; v1 submitted 19 April, 2017;
originally announced April 2017.
-
Exploring Cosmic Origins with CORE: B-mode Component Separation
Authors:
M. Remazeilles,
A. J. Banday,
C. Baccigalupi,
S. Basak,
A. Bonaldi,
G. De Zotti,
J. Delabrouille,
C. Dickinson,
H. K. Eriksen,
J. Errard,
R. Fernandez-Cobos,
U. Fuskeland,
C. Hervías-Caimapo,
M. López-Caniego,
E. Martinez-González,
M. Roman,
P. Vielva,
I. Wehus,
A. Achucarro,
P. Ade,
R. Allison,
M. Ashdown,
M. Ballardini,
R. Banerji,
N. Bartolo
, et al. (91 additional authors not shown)
Abstract:
We demonstrate that, for the baseline design of the CORE satellite mission, the polarized foregrounds can be controlled at the level required to allow the detection of the primordial cosmic microwave background (CMB) $B$-mode polarization with the desired accuracy at both reionization and recombination scales, for tensor-to-scalar ratio values of ${r\gtrsim 5\times 10^{-3}}$. We consider detailed…
▽ More
We demonstrate that, for the baseline design of the CORE satellite mission, the polarized foregrounds can be controlled at the level required to allow the detection of the primordial cosmic microwave background (CMB) $B$-mode polarization with the desired accuracy at both reionization and recombination scales, for tensor-to-scalar ratio values of ${r\gtrsim 5\times 10^{-3}}$. We consider detailed sky simulations based on state-of-the-art CMB observations that consist of CMB polarization with $τ=0.055$ and tensor-to-scalar values ranging from $r=10^{-2}$ to $10^{-3}$, Galactic synchrotron, and thermal dust polarization with variable spectral indices over the sky, polarized anomalous microwave emission, polarized infrared and radio sources, and gravitational lensing effects. Using both parametric and blind approaches, we perform full component separation and likelihood analysis of the simulations, allowing us to quantify both uncertainties and biases on the reconstructed primordial $B$-modes. Under the assumption of perfect control of lensing effects, CORE would measure an unbiased estimate of $r=\left(5 \pm 0.4\right)\times 10^{-3}$ after foreground cleaning. In the presence of both gravitational lensing effects and astrophysical foregrounds, the significance of the detection is lowered, with CORE achieving a $4σ$-measurement of $r=5\times 10^{-3}$ after foreground cleaning and $60$% delensing. For lower tensor-to-scalar ratios ($r=10^{-3}$) the overall uncertainty on $r$ is dominated by foreground residuals, not by the 40% residual of lensing cosmic variance. Moreover, the residual contribution of unprocessed polarized point-sources can be the dominant foreground contamination to primordial B-modes at this $r$ level, even on relatively large angular scales, $\ell \sim 50$. Finally, we report two sources of potential bias for the detection of the primordial $B$-modes.[abridged]
△ Less
Submitted 19 June, 2017; v1 submitted 14 April, 2017;
originally announced April 2017.
-
Exploring Cosmic Origins with CORE: Cluster Science
Authors:
J. -B. Melin,
A. Bonaldi,
M. Remazeilles,
S. Hagstotz,
J. M. Diego,
C. Hernández-Monteagudo,
R. T. Génova-Santos,
G. Luzzi,
C. J. A. P. Martins,
S. Grandis,
J. J. Mohr,
J. G. Bartlett,
J. Delabrouille,
S. Ferraro,
D. Tramonte,
J. A. Rubiño-Martín,
J. F. Macìas-Pérez,
A. Achúcarro,
P. Ade,
R. Allison,
M. Ashdown,
M. Ballardini,
A. J. Banday,
R. Banerji,
N. Bartolo
, et al. (96 additional authors not shown)
Abstract:
We examine the cosmological constraints that can be achieved with a galaxy cluster survey with the future CORE space mission. Using realistic simulations of the millimeter sky, produced with the latest version of the Planck Sky Model, we characterize the CORE cluster catalogues as a function of the main mission performance parameters. We pay particular attention to telescope size, key to improved…
▽ More
We examine the cosmological constraints that can be achieved with a galaxy cluster survey with the future CORE space mission. Using realistic simulations of the millimeter sky, produced with the latest version of the Planck Sky Model, we characterize the CORE cluster catalogues as a function of the main mission performance parameters. We pay particular attention to telescope size, key to improved angular resolution, and discuss the comparison and the complementarity of CORE with ambitious future ground-based CMB experiments that could be deployed in the next decade. A possible CORE mission concept with a 150 cm diameter primary mirror can detect of the order of 50,000 clusters through the thermal Sunyaev-Zeldovich effect (SZE). The total yield increases (decreases) by 25% when increasing (decreasing) the mirror diameter by 30 cm. The 150 cm telescope configuration will detect the most massive clusters ($>10^{14}\, M_\odot$) at redshift $z>1.5$ over the whole sky, although the exact number above this redshift is tied to the uncertain evolution of the cluster SZE flux-mass relation; assuming self-similar evolution, CORE will detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure individual cluster halo masses through lensing of the cosmic microwave background anisotropies with a 1-$σ$ sensitivity of $4\times10^{14} M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm one. [abridged]
△ Less
Submitted 30 March, 2017;
originally announced March 2017.
-
Exploring Cosmic Origins with CORE: Inflation
Authors:
CORE Collaboration,
Fabio Finelli,
Martin Bucher,
Ana Achúcarro,
Mario Ballardini,
Nicola Bartolo,
Daniel Baumann,
Sébastien Clesse,
Josquin Errard,
Will Handley,
Mark Hindmarsh,
Kimmo Kiiveri,
Martin Kunz,
Anthony Lasenby,
Michele Liguori,
Daniela Paoletti,
Christophe Ringeval,
Jussi Väliviita,
Bartjan van Tent,
Vincent Vennin,
Rupert Allison,
Frederico Arroja,
Marc Ashdown,
A. J. Banday,
Ranajoy Banerji
, et al. (107 additional authors not shown)
Abstract:
We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of…
▽ More
We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 μ$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.
△ Less
Submitted 5 April, 2017; v1 submitted 25 December, 2016;
originally announced December 2016.
-
Exploring Cosmic Origins with CORE: Cosmological Parameters
Authors:
Eleonora Di Valentino,
Thejs Brinckmann,
Martina Gerbino,
Vivian Poulin,
François R. Bouchet,
Julien Lesgourgues,
Alessandro Melchiorri,
Jens Chluba,
Sebastien Clesse,
Jacques Delabrouille,
Cora Dvorkin,
Francesco Forastieri,
Silvia Galli,
Deanna C. Hooper,
Massimiliano Lattanzi,
Carlos J. A. P. Martins,
Laura Salvati,
Giovanni Cabass,
Andrea Caputo,
Elena Giusarma,
Eric Hivon,
Paolo Natoli,
Luca Pagano,
Simone Paradiso,
Jose Alberto Rubino-Martin
, et al. (103 additional authors not shown)
Abstract:
We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particul…
▽ More
We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED)
△ Less
Submitted 5 April, 2017; v1 submitted 30 November, 2016;
originally announced December 2016.
-
Exploring Cosmic Origins with CORE: Extragalactic sources in Cosmic Microwave Background maps
Authors:
G. De Zotti,
J. Gonzalez-Nuevo,
M. Lopez-Caniego,
M. Negrello,
J. Greenslade,
C. Hernandez-Monteagudo,
J. Delabrouille,
Z. -Y. Cai,
M. Bonato,
A. Achucarro,
P. Ade,
R. Allison,
M. Ashdown,
M. Ballardini,
A. J. Banday,
R. Banerji,
J. G. Bartlett,
N. Bartolo,
S. Basak,
M. Bersanelli,
M. Biesiada,
M. Bilicki,
A. Bonaldi,
J. Borrill,
F. Bouchet
, et al. (99 additional authors not shown)
Abstract:
We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size…
▽ More
We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.
△ Less
Submitted 18 May, 2017; v1 submitted 23 September, 2016;
originally announced September 2016.
-
Increased signaling entropy in cancer requires the scale-free property of protein interaction networks
Authors:
Andrew E. Teschendorff,
Christopher R. S. Banerji,
Simone Severini,
Reimer Kuehn,
Peter Sollich
Abstract:
One of the key characteristics of cancer cells is an increased phenotypic plasticity, driven by underlying genetic and epigenetic perturbations. However, at a systems-level it is unclear how these perturbations give rise to the observed increased plasticity. Elucidating such systems-level principles is key for an improved understanding of cancer. Recently, it has been shown that signaling entropy,…
▽ More
One of the key characteristics of cancer cells is an increased phenotypic plasticity, driven by underlying genetic and epigenetic perturbations. However, at a systems-level it is unclear how these perturbations give rise to the observed increased plasticity. Elucidating such systems-level principles is key for an improved understanding of cancer. Recently, it has been shown that signaling entropy, an overall measure of signaling pathway promiscuity, and computable from integrating a sample's gene expression profile with a protein interaction network, correlates with phenotypic plasticity and is increased in cancer compared to normal tissue. Here we develop a computational framework for studying the effects of network perturbations on signaling entropy. We demonstrate that the increased signaling entropy of cancer is driven by two factors: (i) the scale-free (or near scale-free) topology of the interaction network, and (ii) a subtle positive correlation between differential gene expression and node connectivity. Indeed, we show that if protein interaction networks were random graphs, described by Poisson degree distributions, that cancer would generally not exhibit an increased signaling entropy. In summary, this work exposes a deep connection between cancer, signaling entropy and interaction network topology.
△ Less
Submitted 1 April, 2015;
originally announced April 2015.
-
Cellular network entropy as the energy potential in Waddington's differentiation landscape
Authors:
Christopher R. S. Banerji,
Diego Miranda-Saavedra,
Simone Severini,
Martin Widschwendter,
Tariq Enver,
Joseph X. Zhou,
Andrew E. Teschendorff
Abstract:
Differentiation is a key cellular process in normal tissue development that is significantly altered in cancer. Although molecular signatures characterising pluripotency and multipotency exist, there is, as yet, no single quantitative mark of a cellular sample's position in the global differentiation hierarchy. Here we adopt a systems view and consider the sample's network entropy, a measure of si…
▽ More
Differentiation is a key cellular process in normal tissue development that is significantly altered in cancer. Although molecular signatures characterising pluripotency and multipotency exist, there is, as yet, no single quantitative mark of a cellular sample's position in the global differentiation hierarchy. Here we adopt a systems view and consider the sample's network entropy, a measure of signaling pathway promiscuity, computable from a sample's genome-wide expression profile. We demonstrate that network entropy provides a quantitative, in-silico, readout of the average undifferentiated state of the profiled cells, recapitulating the known hierarchy of pluripotent, multipotent and differentiated cell types. Network entropy further exhibits dynamic changes in time course differentiation data, and in line with a sample's differentiation stage. In disease, network entropy predicts a higher level of cellular plasticity in cancer stem cell populations compared to ordinary cancer cells. Importantly, network entropy also allows identification of key differentiation pathways. Our results are consistent with the view that pluripotency is a statistical property defined at the cellular population level, correlating with intra-sample heterogeneity, and driven by the degree of signaling promiscuity in cells. In summary, network entropy provides a quantitative measure of a cell's undifferentiated state, defining its elevation in Waddington's landscape.
△ Less
Submitted 26 October, 2013;
originally announced October 2013.
-
A notion of graph likelihood and an infinite monkey theorem
Authors:
Christopher R. S. Banerji,
Toufik Mansour,
Simone Severini
Abstract:
We play with a graph-theoretic analogue of the folklore infinite monkey theorem. We define a notion of graph likelihood as the probability that a given graph is constructed by a monkey in a number of time steps equal to the number of vertices. We present an algorithm to compute this graph invariant and closed formulas for some infinite classes. We have to leave the computational complexity of the…
▽ More
We play with a graph-theoretic analogue of the folklore infinite monkey theorem. We define a notion of graph likelihood as the probability that a given graph is constructed by a monkey in a number of time steps equal to the number of vertices. We present an algorithm to compute this graph invariant and closed formulas for some infinite classes. We have to leave the computational complexity of the likelihood as an open problem.
△ Less
Submitted 12 April, 2013;
originally announced April 2013.
-
Network Transfer Entropy and Metric Space for Causality Inference
Authors:
Christopher R. S. Banerji,
Simone Severini,
Andrew E. Teschendorff
Abstract:
A measure is derived to quantify directed information transfer between pairs of vertices in a weighted network, over paths of a specified maximal length. Our approach employs a general, probabilistic model of network traffic, from which the informational distance between dynamics on two weighted networks can be naturally expressed as a Jensen Shannon Divergence (JSD). Our network transfer entropy…
▽ More
A measure is derived to quantify directed information transfer between pairs of vertices in a weighted network, over paths of a specified maximal length. Our approach employs a general, probabilistic model of network traffic, from which the informational distance between dynamics on two weighted networks can be naturally expressed as a Jensen Shannon Divergence (JSD). Our network transfer entropy measure is shown to be able to distinguish and quantify causal relationships between network elements, in applications to simple synthetic networks and a biological signalling network. We conclude with a theoretical extension of our framework, in which the square root of the JSD induces a metric on the space of dynamics on weighted networks. We prove a convergence criterion, demonstrating that a form of convergence in the structure of weighted networks in a family of matrix metric spaces implies convergence of their dynamics with respect to the square root JSD metric.
△ Less
Submitted 2 November, 2013; v1 submitted 1 March, 2013;
originally announced March 2013.
-
Dissipative Cylindrical Collapse of Charged Anisotropic Fluid
Authors:
Sarbari Guha,
Ranajoy Banerji
Abstract:
We have studied the dynamics of a cylindrical column of anisotropic, charged fluid which is experiencing dissipation in the form of heat flow, free-streaming radiation, and shearing viscosity, undergoing gravitational collapse. We calculate the Einstein-Maxwell field equations and, using the Darmois junction conditions, match the interior non-static cylindrically symmetric space-time with the exte…
▽ More
We have studied the dynamics of a cylindrical column of anisotropic, charged fluid which is experiencing dissipation in the form of heat flow, free-streaming radiation, and shearing viscosity, undergoing gravitational collapse. We calculate the Einstein-Maxwell field equations and, using the Darmois junction conditions, match the interior non-static cylindrically symmetric space-time with the exterior anisotropic, charged, cylindrically symmetric space-time. The behavior of the density, pressure and luminosity of the collapsing matter has been analyzed. From the dynamical equations, the effect of charge and dissipative quantities over the cylindrical collapse are studied. Finally, we have derived the solutions for the collapsing matter which is valid during the later stages of collapse and have discussed the significance from a physical standpoint.
△ Less
Submitted 13 December, 2013; v1 submitted 25 February, 2013;
originally announced February 2013.