Showing 1–50 of 128 results for author: Galloway, M

XENONnT WIMP Search: Signal & Background Modeling and Statistical Inference

Authors: XENON Collaboration, E. Aprile, J. Aalbers, K. Abe, S. Ahmed Maouloud, L. Althueser, B. Andrieu, E. Angelino, D. Antón Martin, F. Arneodo, L. Baudis, M. Bazyk, L. Bellagamba, R. Biondi, A. Bismark, K. Boese, A. Brown, G. Bruno, R. Budnik, J. M. R. Cardoso, A. P. Cimental Chávez, A. P. Colijn, J. Conrad, J. J. Cuenca-García, V. D'Andrea , et al. (139 additional authors not shown)

Abstract: The XENONnT experiment searches for weakly-interacting massive particle (WIMP) dark matter scattering off a xenon nucleus. In particular, XENONnT uses a dual-phase time projection chamber with a 5.9-tonne liquid xenon target, detecting both scintillation and ionization signals to reconstruct the energy, position, and type of recoil. A blind search for nuclear recoil WIMPs with an exposure of 1.1 t… ▽ More The XENONnT experiment searches for weakly-interacting massive particle (WIMP) dark matter scattering off a xenon nucleus. In particular, XENONnT uses a dual-phase time projection chamber with a 5.9-tonne liquid xenon target, detecting both scintillation and ionization signals to reconstruct the energy, position, and type of recoil. A blind search for nuclear recoil WIMPs with an exposure of 1.1 tonne-years yielded no signal excess over background expectations, from which competitive exclusion limits were derived on WIMP-nucleon elastic scatter cross sections, for WIMP masses ranging from 6 GeV/$c^2$ up to the TeV/$c^2$ scale. This work details the modeling and statistical methods employed in this search. By means of calibration data, we model the detector response, which is then used to derive background and signal models. The construction and validation of these models is discussed, alongside additional purely data-driven backgrounds. We also describe the statistical inference framework, including the definition of the likelihood function and the construction of confidence intervals. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: 20 pages, 10 figures

The LiteBIRD mission to explore cosmic inflation

Authors: T. Ghigna, A. Adler, K. Aizawa, H. Akamatsu, R. Akizawa, E. Allys, A. Anand, J. Aumont, J. Austermann, S. Azzoni, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, A. Basyrov, S. Beckman, M. Bersanelli, M. Bortolami, F. Bouchet, T. Brinckmann, P. Campeti, E. Carinos, A. Carones , et al. (134 additional authors not shown)

Abstract: LiteBIRD, the next-generation cosmic microwave background (CMB) experiment, aims for a launch in Japan's fiscal year 2032, marking a major advancement in the exploration of primordial cosmology and fundamental physics. Orbiting the Sun-Earth Lagrangian point L2, this JAXA-led strategic L-class mission will conduct a comprehensive mapping of the CMB polarization across the entire sky. During its 3-… ▽ More LiteBIRD, the next-generation cosmic microwave background (CMB) experiment, aims for a launch in Japan's fiscal year 2032, marking a major advancement in the exploration of primordial cosmology and fundamental physics. Orbiting the Sun-Earth Lagrangian point L2, this JAXA-led strategic L-class mission will conduct a comprehensive mapping of the CMB polarization across the entire sky. During its 3-year mission, LiteBIRD will employ three telescopes within 15 unique frequency bands (ranging from 34 through 448 GHz), targeting a sensitivity of 2.2\,$μ$K-arcmin and a resolution of 0.5$^\circ$ at 100\,GHz. Its primary goal is to measure the tensor-to-scalar ratio $r$ with an uncertainty $δr = 0.001$, including systematic errors and margin. If $r \geq 0.01$, LiteBIRD expects to achieve a $>5σ$ detection in the $\ell=$2-10 and $\ell=$11-200 ranges separately, providing crucial insight into the early Universe. We describe LiteBIRD's scientific objectives, the application of systems engineering to mission requirements, the anticipated scientific impact, and the operations and scanning strategies vital to minimizing systematic effects. We will also highlight LiteBIRD's synergies with concurrent CMB projects. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 23 pages, 9 figures, 1 table, SPIE Astronomical Telescopes + Instrumentation 2024

Cosmoglobe DR2. II. CIB monopole measurements from COBE-DIRBE through global Bayesian analysis

Authors: D. J. Watts, M. Galloway, E. Gjerløw, M. San, R. Aurlien, A. Basyrov, M. Brilenkov, H. K. Eriksen, U. Fuskeland, D. Herman, H. T. Ihle, J. G. S. Lunde, S. K. Næss, N. -O. Stutzer, H. Thommesen, I. K. Wehus

Abstract: We derive new constraints on the CIB monopole spectrum from reprocessed COBE-DIRBE sky maps with lower instrumental and astrophysical contamination than the legacy DIRBE maps. These maps have been generated through a global Bayesian analysis framework that simultaneously fits cosmological, astrophysical, and instrumental parameters, described in a series of papers referred to as Cosmoglobe DR2. We… ▽ More We derive new constraints on the CIB monopole spectrum from reprocessed COBE-DIRBE sky maps with lower instrumental and astrophysical contamination than the legacy DIRBE maps. These maps have been generated through a global Bayesian analysis framework that simultaneously fits cosmological, astrophysical, and instrumental parameters, described in a series of papers referred to as Cosmoglobe DR2. We have applied this method to the DIRBE Calibrated Individual Observations, complemented by selected HFI and FIRAS sky maps to break key astrophysical degeneracies, as well as WISE and Gaia compact object catalogs. In this paper, we focus on CIB constraints that result from this work. We report detections of an isotropic signal in five out of the ten DIRBE bands (1.25, 2.2, 3.5, 140, and 240 $μ$m). For the 2.2 $μ$m channel, we find an amplitude of $13\pm3\,\mathrm{nW\,m^{-2}\,sr^{-1}}$, 36% lower than reported from the official DIRBE release. For the 240 $μ$m channel, we find $10\pm3\,\mathrm{nW\,m^{-2}\,sr^{-1}}$, 26% lower than the official DIRBE release. We interpret these lower values as resulting from improved zodiacal light and Galactic foreground modeling. For the bands between 4.9 and 100 $μ$m, the presence of excess radiation in solar-centric coordinates reported in a companion paper precludes the definition of lower limits. However, the analysis still provides well-defined upper limits. For the 12 $μ$m channel, we find an upper 95% confidence limit of 103 $\mathrm{nW\,m^{-2}\,sr^{-1}}$, a factor of four lower than the corresponding legacy result of 468 $\mathrm{nW\,m^{-2}\,sr^{-1}}$. The results presented in this paper redefines the state-of-the-art CIB monopole constraints from COBE-DIRBE, and it provides a real-world illustration of the power of global end-to-end analysis of multiple complementary data sets which is the foundational idea of the Cosmoglobe project. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 14 pages, 7 figures, submitted to A&A

LiteBIRD Science Goals and Forecasts: Primordial Magnetic Fields

Authors: D. Paoletti, J. Rubino-Martin, M. Shiraishi, D. Molinari, J. Chluba, F. Finelli, C. Baccigalupi, J. Errard, A. Gruppuso, A. I. Lonappan, A. Tartari, E. Allys, A. Anand, J. Aumont, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, M. Bersanelli, M. Bortolami, T. Brinckmann, E. Calabrese, P. Campeti, A. Carones, F. J. Casas , et al. (75 additional authors not shown)

Abstract: We present detailed forecasts for the constraints on primordial magnetic fields (PMFs) that will be obtained with the LiteBIRD satellite. The constraints are driven by the effects of PMFs on the CMB anisotropies: the gravitational effects of magnetically-induced perturbations; the effects on the thermal and ionization history of the Universe; the Faraday rotation imprint on the CMB polarization; a… ▽ More We present detailed forecasts for the constraints on primordial magnetic fields (PMFs) that will be obtained with the LiteBIRD satellite. The constraints are driven by the effects of PMFs on the CMB anisotropies: the gravitational effects of magnetically-induced perturbations; the effects on the thermal and ionization history of the Universe; the Faraday rotation imprint on the CMB polarization; and the non-Gaussianities induced in polarization anisotropies. LiteBIRD represents a sensitive probe for PMFs and by exploiting all the physical effects, it will be able to improve the current limit coming from Planck. In particular, thanks to its accurate $B$-mode polarization measurement, LiteBIRD will improve the constraints on infrared configurations for the gravitational effect, giving $B_{\rm 1\,Mpc}^{n_{\rm B} =-2.9} < 0.8$ nG at 95% C.L., potentially opening the possibility to detect nanogauss fields with high significance. We also observe a significant improvement in the limits when marginalized over the spectral index, $B_{1\,{\rm Mpc}}^{\rm marg}< 2.2$ nG at 95% C.L. From the thermal history effect, which relies mainly on $E$-mode polarization data, we obtain a significant improvement for all PMF configurations, with the marginalized case, $\sqrt{\langle B^2\rangle}^{\rm marg}<0.50$ nG at 95% C.L. Faraday rotation constraints will take advantage of the wide frequency coverage of LiteBIRD and the high sensitivity in $B$ modes, improving the limits by orders of magnitude with respect to current results, $B_{1\,{\rm Mpc}}^{n_{\rm B} =-2.9} < 3.2$ nG at 95% C.L. Finally, non-Gaussianities of the $B$-mode polarization can probe PMFs at the level of 1 nG, again significantly improving the current bounds from Planck. Altogether our forecasts represent a broad collection of complementary probes, providing conservative limits on PMF characteristics that will be achieved with LiteBIRD. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 51 pages, 24 figures, abstract shortened

Offline tagging of radon-induced backgrounds in XENON1T and applicability to other liquid xenon detectors

Authors: E. Aprile, J. Aalbers, K. Abe, S. Ahmed Maouloud, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, D. Antón Martin, F. Arneodo, L. Baudis, A. L. Baxter, M. Bazyk, L. Bellagamba, R. Biondi, A. Bismark, E. J. Brookes, A. Brown, G. Bruno, R. Budnik, T. K. Bui, J. M. R. Cardoso, A. P. Cimental Chavez, A. P. Colijn, J. Conrad , et al. (142 additional authors not shown)

Abstract: This paper details the first application of a software tagging algorithm to reduce radon-induced backgrounds in liquid noble element time projection chambers, such as XENON1T and XENONnT. The convection velocity field in XENON1T was mapped out using $^{222}\text{Rn}$ and $^{218}\text{Po}$ events, and the root-mean-square convection speed was measured to be $0.30 \pm 0.01$ cm/s. Given this velocity… ▽ More This paper details the first application of a software tagging algorithm to reduce radon-induced backgrounds in liquid noble element time projection chambers, such as XENON1T and XENONnT. The convection velocity field in XENON1T was mapped out using $^{222}\text{Rn}$ and $^{218}\text{Po}$ events, and the root-mean-square convection speed was measured to be $0.30 \pm 0.01$ cm/s. Given this velocity field, $^{214}\text{Pb}$ background events can be tagged when they are followed by $^{214}\text{Bi}$ and $^{214}\text{Po}$ decays, or preceded by $^{218}\text{Po}$ decays. This was achieved by evolving a point cloud in the direction of a measured convection velocity field, and searching for $^{214}\text{Bi}$ and $^{214}\text{Po}$ decays or $^{218}\text{Po}$ decays within a volume defined by the point cloud. In XENON1T, this tagging system achieved a $^{214}\text{Pb}$ background reduction of $6.2^{+0.4}_{-0.9}\%$ with an exposure loss of $1.8\pm 0.2 \%$, despite the timescales of convection being smaller than the relevant decay times. We show that the performance can be improved in XENONnT, and that the performance of such a software-tagging approach can be expected to be further improved in a diffusion-limited scenario. Finally, a similar method might be useful to tag the cosmogenic $^{137}\text{Xe}$ background, which is relevant to the search for neutrinoless double-beta decay. △ Less

Submitted 19 June, 2024; v1 submitted 21 March, 2024; originally announced March 2024.

Comments: 17 pages, 19 figures

The XENONnT Dark Matter Experiment

Authors: XENON Collaboration, E. Aprile, J. Aalbers, K. Abe, S. Ahmed Maouloud, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, V. C. Antochi, D. Antón Martin, F. Arneodo, M. Balata, L. Baudis, A. L. Baxter, M. Bazyk, L. Bellagamba, R. Biondi, A. Bismark, E. J. Brookes, A. Brown, S. Bruenner, G. Bruno, R. Budnik, T. K. Bui , et al. (170 additional authors not shown)

Abstract: The multi-staged XENON program at INFN Laboratori Nazionali del Gran Sasso aims to detect dark matter with two-phase liquid xenon time projection chambers of increasing size and sensitivity. The XENONnT experiment is the latest detector in the program, planned to be an upgrade of its predecessor XENON1T. It features an active target of 5.9 tonnes of cryogenic liquid xenon (8.5 tonnes total mass in… ▽ More The multi-staged XENON program at INFN Laboratori Nazionali del Gran Sasso aims to detect dark matter with two-phase liquid xenon time projection chambers of increasing size and sensitivity. The XENONnT experiment is the latest detector in the program, planned to be an upgrade of its predecessor XENON1T. It features an active target of 5.9 tonnes of cryogenic liquid xenon (8.5 tonnes total mass in cryostat). The experiment is expected to extend the sensitivity to WIMP dark matter by more than an order of magnitude compared to XENON1T, thanks to the larger active mass and the significantly reduced background, improved by novel systems such as a radon removal plant and a neutron veto. This article describes the XENONnT experiment and its sub-systems in detail and reports on the detector performance during the first science run. △ Less

Submitted 15 February, 2024; originally announced February 2024.

Comments: 32 pages, 19 figures

Impact of beam far side-lobe knowledge in the presence of foregrounds for LiteBIRD

Authors: C. Leloup, G. Patanchon, J. Errard, C. Franceschet, J. E. Gudmundsson, S. Henrot-Versillé, H. Imada, H. Ishino, T. Matsumura, G. Puglisi, W. Wang, A. Adler, J. Aumont, R. Aurlien, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, A. Basyrov, M. Bersanelli, D. Blinov, M. Bortolami, T. Brinckmann, P. Campeti , et al. (86 additional authors not shown)

Abstract: We present a study of the impact of an uncertainty in the beam far side-lobe knowledge on the measurement of the Cosmic Microwave Background $B$-mode signal at large scale. It is expected to be one of the main source of systematic effects in future CMB observations. Because it is crucial for all-sky survey missions to take into account the interplays between beam systematic effects and all the dat… ▽ More We present a study of the impact of an uncertainty in the beam far side-lobe knowledge on the measurement of the Cosmic Microwave Background $B$-mode signal at large scale. It is expected to be one of the main source of systematic effects in future CMB observations. Because it is crucial for all-sky survey missions to take into account the interplays between beam systematic effects and all the data analysis steps, the primary goal of this paper is to provide the methodology to carry out the end-to-end study of their effect for a space-borne CMB polarization experiment, up to the cosmological results in the form of a bias $δr$ on the tensor-to-scalar ratio $r$. LiteBIRD is dedicated to target the measurement of CMB primordial $B$ modes by reaching a sensitivity of $σ\left( r \right) \leq 10^{-3}$ assuming $r=0$. As a demonstration of our framework, we derive the relationship between the knowledge of the beam far side-lobes and the tentatively allocated error budget under given assumptions on design, simulation and component separation method. We assume no mitigation of the far side-lobes effect at any stage of the analysis pipeline. We show that $δr$ is mostly due to the integrated fractional power difference between the estimated beams and the true beams in the far side-lobes region, with little dependence on the actual shape of the beams, for low enough $δr$. Under our set of assumptions, in particular considering the specific foreground cleaning method we used, we find that the integrated fractional power in the far side-lobes should be known at a level as tight as $\sim 10^{-4}$, to achieve the required limit on the bias $δr < 1.9 \times 10^{-5}$. The framework and tools developed for this study can be easily adapted to provide requirements under different design, data analysis frameworks and for other future space-borne experiments beyond LiteBIRD. △ Less

Submitted 14 December, 2023; originally announced December 2023.

LiteBIRD Science Goals and Forecasts: Improving Sensitivity to Inflationary Gravitational Waves with Multitracer Delensing

Authors: T. Namikawa, A. I. Lonappan, C. Baccigalupi, N. Bartolo, D. Beck, K. Benabed, A. Challinor, P. Diego-Palazuelos, J. Errard, S. Farrens, A. Gruppuso, N. Krachmalnicoff, M. Migliaccio, E. Martínez-González, V. Pettorino, G. Piccirilli, M. Ruiz-Granda, B. Sherwin, J. Starck, P. Vielva, R. Akizawa, A. Anand, J. Aumont, R. Aurlien, S. Azzoni , et al. (97 additional authors not shown)

Abstract: We estimate the efficiency of mitigating the lensing $B$-mode polarization, the so-called delensing, for the $LiteBIRD$ experiment with multiple external data sets of lensing-mass tracers. The current best bound on the tensor-to-scalar ratio, $r$, is limited by lensing rather than Galactic foregrounds. Delensing will be a critical step to improve sensitivity to $r$ as measurements of $r$ become mo… ▽ More We estimate the efficiency of mitigating the lensing $B$-mode polarization, the so-called delensing, for the $LiteBIRD$ experiment with multiple external data sets of lensing-mass tracers. The current best bound on the tensor-to-scalar ratio, $r$, is limited by lensing rather than Galactic foregrounds. Delensing will be a critical step to improve sensitivity to $r$ as measurements of $r$ become more and more limited by lensing. In this paper, we extend the analysis of the recent $LiteBIRD$ forecast paper to include multiple mass tracers, i.e., the CMB lensing maps from $LiteBIRD$ and CMB-S4-like experiment, cosmic infrared background, and galaxy number density from $Euclid$- and LSST-like survey. We find that multi-tracer delensing will further improve the constraint on $r$ by about $20\%$. In $LiteBIRD$, the residual Galactic foregrounds also significantly contribute to uncertainties of the $B$-modes, and delensing becomes more important if the residual foregrounds are further reduced by an improved component separation method. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: 21 pages, 7 figures

LiteBIRD Science Goals and Forecasts: A full-sky measurement of gravitational lensing of the CMB

Authors: A. I. Lonappan, T. Namikawa, G. Piccirilli, P. Diego-Palazuelos, M. Ruiz-Granda, M. Migliaccio, C. Baccigalupi, N. Bartolo, D. Beck, K. Benabed, A. Challinor, J. Errard, S. Farrens, A. Gruppuso, N. Krachmalnicoff, E. Martínez-González, V. Pettorino, B. Sherwin, J. Starck, P. Vielva, R. Akizawa, A. Anand, J. Aumont, R. Aurlien, S. Azzoni , et al. (97 additional authors not shown)

Abstract: We explore the capability of measuring lensing signals in $LiteBIRD$ full-sky polarization maps. With a $30$ arcmin beam width and an impressively low polarization noise of $2.16\,μ$K-arcmin, $LiteBIRD$ will be able to measure the full-sky polarization of the cosmic microwave background (CMB) very precisely. This unique sensitivity also enables the reconstruction of a nearly full-sky lensing map u… ▽ More We explore the capability of measuring lensing signals in $LiteBIRD$ full-sky polarization maps. With a $30$ arcmin beam width and an impressively low polarization noise of $2.16\,μ$K-arcmin, $LiteBIRD$ will be able to measure the full-sky polarization of the cosmic microwave background (CMB) very precisely. This unique sensitivity also enables the reconstruction of a nearly full-sky lensing map using only polarization data, even considering its limited capability to capture small-scale CMB anisotropies. In this paper, we investigate the ability to construct a full-sky lensing measurement in the presence of Galactic foregrounds, finding that several possible biases from Galactic foregrounds should be negligible after component separation by harmonic-space internal linear combination. We find that the signal-to-noise ratio of the lensing is approximately $40$ using only polarization data measured over $90\%$ of the sky. This achievement is comparable to $Planck$'s recent lensing measurement with both temperature and polarization and represents a four-fold improvement over $Planck$'s polarization-only lensing measurement. The $LiteBIRD$ lensing map will complement the $Planck$ lensing map and provide several opportunities for cross-correlation science, especially in the northern hemisphere. △ Less

Submitted 8 December, 2023; originally announced December 2023.

LiteBIRD Science Goals and Forecasts. A Case Study of the Origin of Primordial Gravitational Waves using Large-Scale CMB Polarization

Authors: P. Campeti, E. Komatsu, C. Baccigalupi, M. Ballardini, N. Bartolo, A. Carones, J. Errard, F. Finelli, R. Flauger, S. Galli, G. Galloni, S. Giardiello, M. Hazumi, S. Henrot-Versillé, L. T. Hergt, K. Kohri, C. Leloup, J. Lesgourgues, J. Macias-Perez, E. Martínez-González, S. Matarrese, T. Matsumura, L. Montier, T. Namikawa, D. Paoletti , et al. (85 additional authors not shown)

Abstract: We study the possibility of using the $LiteBIRD$ satellite $B$-mode survey to constrain models of inflation producing specific features in CMB angular power spectra. We explore a particular model example, i.e. spectator axion-SU(2) gauge field inflation. This model can source parity-violating gravitational waves from the amplification of gauge field fluctuations driven by a pseudoscalar "axionlike… ▽ More We study the possibility of using the $LiteBIRD$ satellite $B$-mode survey to constrain models of inflation producing specific features in CMB angular power spectra. We explore a particular model example, i.e. spectator axion-SU(2) gauge field inflation. This model can source parity-violating gravitational waves from the amplification of gauge field fluctuations driven by a pseudoscalar "axionlike" field, rolling for a few e-folds during inflation. The sourced gravitational waves can exceed the vacuum contribution at reionization bump scales by about an order of magnitude and can be comparable to the vacuum contribution at recombination bump scales. We argue that a satellite mission with full sky coverage and access to the reionization bump scales is necessary to understand the origin of the primordial gravitational wave signal and distinguish among two production mechanisms: quantum vacuum fluctuations of spacetime and matter sources during inflation. We present the expected constraints on model parameters from $LiteBIRD$ satellite simulations, which complement and expand previous studies in the literature. We find that $LiteBIRD$ will be able to exclude with high significance standard single-field slow-roll models, such as the Starobinsky model, if the true model is the axion-SU(2) model with a feature at CMB scales. We further investigate the possibility of using the parity-violating signature of the model, such as the $TB$ and $EB$ angular power spectra, to disentangle it from the standard single-field slow-roll scenario. We find that most of the discriminating power of $LiteBIRD$ will reside in $BB$ angular power spectra rather than in $TB$ and $EB$ correlations. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Comments: 22 pages, 13 figures. Submitted to JCAP

Cosmoglobe DR1. III. First full-sky model of polarized synchrotron emission from all WMAP and Planck LFI data

Authors: D. J. Watts, U. Fuskeland, R. Aurlien, A. Basyrov, L. A. Bianchi, M. Brilenkov, H. K. Eriksen, K. S. F. Fornazier, M. Galloway, E. Gjerløw, B. Hensley, L. T. Hergt, D. Herman, H. Ihle, K. Lee, J. G. S. Lunde, S. K. Nerval, M. San, N. O. Stutzer, H. Thommesen, I. K. Wehus

Abstract: We present the first model of full-sky polarized synchrotron emission that is derived from all WMAP and Planck LFI frequency maps. The basis of this analysis is the set of end-to-end reprocessed Cosmoglobe Data Release 1 sky maps presented in a companion paper, which have significantly lower instrumental systematics than the legacy products from each experiment. We find that the resulting polarize… ▽ More We present the first model of full-sky polarized synchrotron emission that is derived from all WMAP and Planck LFI frequency maps. The basis of this analysis is the set of end-to-end reprocessed Cosmoglobe Data Release 1 sky maps presented in a companion paper, which have significantly lower instrumental systematics than the legacy products from each experiment. We find that the resulting polarized synchrotron amplitude map has an average noise rms of $3.2\,\mathrm{μK}$ at 30 GHz and $2^{\circ}$ FWHM, which is 30% lower than the recently released BeyondPlanck model that included only LFI+WMAP Ka-V data, and 29% lower than the WMAP K-band map alone. The mean $B$-to-$E$ power spectrum ratio is $0.40\pm0.02$, with amplitudes consistent with those measured previously by Planck and QUIJOTE. Assuming a power law model for the synchrotron spectral energy distribution, and using the $T$--$T$ plot method, we find a full-sky inverse noise-variance weighted mean of $β_{\mathrm{s}}=-3.07\pm0.07$ between Cosmoglobe DR1 K-band and 30 GHz, in good agreement with previous estimates. In summary, the novel Cosmoglobe DR1 synchrotron model is both more sensitive and systematically cleaner than similar previous models, and it has a more complete error description that is defined by a set of Monte Carlo posterior samples. We believe that these products are preferable over previous Planck and WMAP products for all synchrotron-related scientific applications, including simulation, forecasting and component separation. △ Less

Submitted 20 October, 2023; originally announced October 2023.

Comments: 15 pages, 15 figures, submitted to A&A

Journal ref: A&A 686, A297 (2024)

Design and performance of the field cage for the XENONnT experiment

Authors: E. Aprile, K. Abe, S. Ahmed Maouloud, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, V. C. Antochi, D. Antón Martin, F. Arneodo, L. Baudis, A. L. Baxter, M. Bazyk, L. Bellagamba, R. Biondi, A. Bismark, E. J. Brookes, A. Brown, S. Bruenner, G. Bruno, R. Budnik, T. K. Bui, C. Cai, J. M. R. Cardoso, D. Cichon , et al. (139 additional authors not shown)

Abstract: The precision in reconstructing events detected in a dual-phase time projection chamber depends on an homogeneous and well understood electric field within the liquid target. In the XENONnT TPC the field homogeneity is achieved through a double-array field cage, consisting of two nested arrays of field shaping rings connected by an easily accessible resistor chain. Rather than being connected to t… ▽ More The precision in reconstructing events detected in a dual-phase time projection chamber depends on an homogeneous and well understood electric field within the liquid target. In the XENONnT TPC the field homogeneity is achieved through a double-array field cage, consisting of two nested arrays of field shaping rings connected by an easily accessible resistor chain. Rather than being connected to the gate electrode, the topmost field shaping ring is independently biased, adding a degree of freedom to tune the electric field during operation. Two-dimensional finite element simulations were used to optimize the field cage, as well as its operation. Simulation results were compared to ${}^{83m}\mathrm{Kr}$ calibration data. This comparison indicates an accumulation of charge on the panels of the TPC which is constant over time, as no evolution of the reconstructed position distribution of events is observed. The simulated electric field was then used to correct the charge signal for the field dependence of the charge yield. This correction resolves the inconsistent measurement of the drift electron lifetime when using different calibrations sources and different field cage tuning voltages. △ Less

Submitted 21 September, 2023; originally announced September 2023.

Journal ref: Eur. Phys. J. C 84, 138 (2024)

Cosmogenic background simulations for the DARWIN observatory at different underground locations

Authors: M. Adrover, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, B. Antunovic, E. Aprile, M. Babicz, D. Bajpai, E. Barberio, L. Baudis, M. Bazyk, N. Bell, L. Bellagamba, R. Biondi, Y. Biondi, A. Bismark, C. Boehm, A. Breskin, E. J. Brookes, A. Brown, G. Bruno, R. Budnik, C. Capelli, J. M. R. Cardoso , et al. (158 additional authors not shown)

Abstract: Xenon dual-phase time projections chambers (TPCs) have proven to be a successful technology in studying physical phenomena that require low-background conditions. With 40t of liquid xenon (LXe) in the TPC baseline design, DARWIN will have a high sensitivity for the detection of particle dark matter, neutrinoless double beta decay ($0νββ$), and axion-like particles (ALPs). Although cosmic muons are… ▽ More Xenon dual-phase time projections chambers (TPCs) have proven to be a successful technology in studying physical phenomena that require low-background conditions. With 40t of liquid xenon (LXe) in the TPC baseline design, DARWIN will have a high sensitivity for the detection of particle dark matter, neutrinoless double beta decay ($0νββ$), and axion-like particles (ALPs). Although cosmic muons are a source of background that cannot be entirely eliminated, they may be greatly diminished by placing the detector deep underground. In this study, we used Monte Carlo simulations to model the cosmogenic background expected for the DARWIN observatory at four underground laboratories: Laboratori Nazionali del Gran Sasso (LNGS), Sanford Underground Research Facility (SURF), Laboratoire Souterrain de Modane (LSM) and SNOLAB. We determine the production rates of unstable xenon isotopes and tritium due to muon-included neutron fluxes and muon-induced spallation. These are expected to represent the dominant contributions to cosmogenic backgrounds and thus the most relevant for site selection. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Cosmoglobe: Towards end-to-end CMB cosmological parameter estimation without likelihood approximations

Authors: J. R. Eskilt, K. Lee, D. J. Watts, V. Anshul, R. Aurlien, A. Basyrov, M. Bersanelli, L. P. L. Colombo, H. K. Eriksen, K. S. F. Fornazier, U. Fuskeland, M. Galloway, E. Gjerløw, L. T. Hergt, H. T. Ihle, J. G. S. Lunde, A. Marins, S. K. Nerval, S. Paradiso, F. Rahman, M. San, N. -O. Stutzer, I. K. Wehus

Abstract: We implement support for a cosmological parameter estimation algorithm as proposed by Racine et al. (2016) in Commander, and quantify its computational efficiency and cost. For a semi-realistic simulation similar to Planck LFI 70 GHz, we find that the computational cost of producing one single sample is about 20 CPU-hours and that the typical Markov chain correlation length is $\sim$100 samples. T… ▽ More We implement support for a cosmological parameter estimation algorithm as proposed by Racine et al. (2016) in Commander, and quantify its computational efficiency and cost. For a semi-realistic simulation similar to Planck LFI 70 GHz, we find that the computational cost of producing one single sample is about 20 CPU-hours and that the typical Markov chain correlation length is $\sim$100 samples. The net effective cost per independent sample is $\sim$2 000 CPU-hours, in comparison with all low-level processing costs of 812 CPU-hours for Planck LFI and WMAP in Cosmoglobe Data Release 1. Thus, although technically possible to run already in its current state, future work should aim to reduce the effective cost per independent sample by one order of magnitude to avoid excessive runtimes, for instance through multi-grid preconditioners and/or derivative-based Markov chain sampling schemes. This work demonstrates the computational feasibility of true Bayesian cosmological parameter estimation with end-to-end error propagation for high-precision CMB experiments without likelihood approximations, but it also highlights the need for additional optimizations before it is ready for full production-level analysis. △ Less

Submitted 31 October, 2023; v1 submitted 27 June, 2023; originally announced June 2023.

Comments: 10 pages, 8 figures. Accepted in A&A

Journal ref: A&A 678, A169 (2023)

Search for events in XENON1T associated with Gravitational Waves

Authors: XENON Collaboration, E. Aprile, K. Abe, S. Ahmed Maouloud, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, V. C. Antochi, D. Antoń Martin, F. Arneodo, L. Baudis, A. L. Baxter, M. Bazyk, L. Bellagamba, R. Biondi, A. Bismark, E. J. Brookes, A. Brown, S. Bruenner, G. Bruno, R. Budnik, T. K. Bui, C. Cai, J. M. R. Cardoso , et al. (138 additional authors not shown)

Abstract: We perform a blind search for particle signals in the XENON1T dark matter detector that occur close in time to gravitational wave signals in the LIGO and Virgo observatories. No particle signal is observed in the nuclear recoil, electronic recoil, CE$ν$NS, and S2-only channels within $\pm$ 500 seconds of observations of the gravitational wave signals GW170104, GW170729, GW170817, GW170818, and GW1… ▽ More We perform a blind search for particle signals in the XENON1T dark matter detector that occur close in time to gravitational wave signals in the LIGO and Virgo observatories. No particle signal is observed in the nuclear recoil, electronic recoil, CE$ν$NS, and S2-only channels within $\pm$ 500 seconds of observations of the gravitational wave signals GW170104, GW170729, GW170817, GW170818, and GW170823. We use this null result to constrain mono-energetic neutrinos and Beyond Standard Model particles emitted in the closest coalescence GW170817, a binary neutron star merger. We set new upper limits on the fluence (time-integrated flux) of coincident neutrinos down to 17 keV at 90% confidence level. Furthermore, we constrain the product of coincident fluence and cross section of Beyond Standard Model particles to be less than $10^{-29}$ cm$^2$/cm$^2$ in the [5.5-210] keV energy range at 90% confidence level. △ Less

Submitted 27 October, 2023; v1 submitted 20 June, 2023; originally announced June 2023.

Cosmoglobe DR1 results. II. Constraints on isotropic cosmic birefringence from reprocessed WMAP and Planck LFI data

Authors: J. R. Eskilt, D. J. Watts, R. Aurlien, A. Basyrov, M. Bersanelli, M. Brilenkov, L. P. L. Colombo, H. K. Eriksen, K. S. F. Fornazier, C. Franceschet, U. Fuskeland, M. Galloway, E. Gjerløw, B. Hensley, L. T. Hergt, D. Herman, H. T. Ihle, K. Lee, J. G. S. Lunde, S. K. Nerval, S. Paradiso, S. K. Patel, F. Rahman, M. Regnier, M. San , et al. (6 additional authors not shown)

Abstract: Cosmic birefringence is a parity-violating effect that might have rotated the plane of linearly polarized light of the cosmic microwave background (CMB) by an angle $β$ since its emission. This has recently been measured to be non-zero at a statistical significance of $3.6σ$ in the official Planck PR4 and 9-year WMAP data. In this work, we constrain $β$ using the reprocessed BeyondPlanck LFI and C… ▽ More Cosmic birefringence is a parity-violating effect that might have rotated the plane of linearly polarized light of the cosmic microwave background (CMB) by an angle $β$ since its emission. This has recently been measured to be non-zero at a statistical significance of $3.6σ$ in the official Planck PR4 and 9-year WMAP data. In this work, we constrain $β$ using the reprocessed BeyondPlanck LFI and Cosmoglobe DR1 WMAP polarization maps. These novel maps have both lower systematic residuals and a more complete error description than the corresponding official products. Foreground $EB$ correlations could bias measurements of $β$, and while thermal dust $EB$ emission has been argued to be statistically non-zero, no evidence for synchrotron $EB$ power has been reported. Unlike the dust-dominated Planck HFI maps, the majority of the LFI and WMAP polarization maps are instead dominated by synchrotron emission. Simultaneously constraining $β$ and the polarization miscalibration angle, $α$, of each channel, we find a best-fit value of $β=0.35^{\circ}\pm0.70^{\circ}$ with LFI and WMAP data only. When including the Planck HFI PR4 maps, but fitting $β$ separately for dust-dominated, $β_{>70\,\mathrm{GHz}}$, and synchrotron-dominated channels, $β_{\leq 70\,\mathrm{GHz}}$, we find $β_{\leq 70\,\mathrm{GHz}}=0.53^{\circ}\pm0.28^\circ$. This differs from zero with a statistical significance of $1.9σ$, and the main contribution to this value comes from the LFI 70 GHz channel. While the statistical significances of these results are low on their own, the measurement derived from the LFI and WMAP synchrotron-dominated maps agrees with the previously reported HFI-dominated constraints, despite the very different astrophysical and instrumental systematics involved in all these experiments. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Comments: 10 pages, 7 figures, 2 tables. Submitted to A&A

Journal ref: A&A 679, A144 (2023)

Searching for Heavy Dark Matter near the Planck Mass with XENON1T

Authors: E. Aprile, K. Abe, S. Ahmed Maouloud, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, V. C. Antochi, D. Antón Martin, F. Arneodo, L. Baudis, A. L. Baxter, M. Bazyk, L. Bellagamba, R. Biondi, A. Bismark, E. J. Brookes, A. Brown, S. Bruenner, G. Bruno, R. Budnik, T. K. Bui, C. Cai, J. M. R. Cardoso, D. Cichon , et al. (142 additional authors not shown)

Abstract: Multiple viable theoretical models predict heavy dark matter particles with a mass close to the Planck mass, a range relatively unexplored by current experimental measurements. We use 219.4 days of data collected with the XENON1T experiment to conduct a blind search for signals from Multiply-Interacting Massive Particles (MIMPs). Their unique track signature allows a targeted analysis with only 0.… ▽ More Multiple viable theoretical models predict heavy dark matter particles with a mass close to the Planck mass, a range relatively unexplored by current experimental measurements. We use 219.4 days of data collected with the XENON1T experiment to conduct a blind search for signals from Multiply-Interacting Massive Particles (MIMPs). Their unique track signature allows a targeted analysis with only 0.05 expected background events from muons. Following unblinding, we observe no signal candidate events. This work places strong constraints on spin-independent interactions of dark matter particles with a mass between 1$\times$10$^{12}\,$GeV/c$^2$ and 2$\times$10$^{17}\,$GeV/c$^2$. In addition, we present the first exclusion limits on spin-dependent MIMP-neutron and MIMP-proton cross-sections for dark matter particles with masses close to the Planck scale. △ Less

Submitted 21 April, 2023; originally announced April 2023.

Comments: 7 pages, 6 figures

Journal ref: Phys. Rev. Lett. 130, 261002 (2023)

Detector signal characterization with a Bayesian network in XENONnT

Authors: XENON Collaboration, E. Aprile, K. Abe, S. Ahmed Maouloud, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, V. C. Antochi, D. Antón Martin, F. Arneodo, L. Baudis, A. L. Baxter, M. Bazyk, L. Bellagamba, R. Biondi, A. Bismark, E. J. Brookes, A. Brown, S. Bruenner, G. Bruno, R. Budnik, T. K. Bui, C. Cai, J. M. R. Cardoso , et al. (142 additional authors not shown)

Abstract: We developed a detector signal characterization model based on a Bayesian network trained on the waveform attributes generated by a dual-phase xenon time projection chamber. By performing inference on the model, we produced a quantitative metric of signal characterization and demonstrate that this metric can be used to determine whether a detector signal is sourced from a scintillation or an ioniz… ▽ More We developed a detector signal characterization model based on a Bayesian network trained on the waveform attributes generated by a dual-phase xenon time projection chamber. By performing inference on the model, we produced a quantitative metric of signal characterization and demonstrate that this metric can be used to determine whether a detector signal is sourced from a scintillation or an ionization process. We describe the method and its performance on electronic-recoil (ER) data taken during the first science run of the XENONnT dark matter experiment. We demonstrate the first use of a Bayesian network in a waveform-based analysis of detector signals. This method resulted in a 3% increase in ER event-selection efficiency with a simultaneously effective rejection of events outside of the region of interest. The findings of this analysis are consistent with the previous analysis from XENONnT, namely a background-only fit of the ER data. △ Less

Submitted 26 July, 2023; v1 submitted 11 April, 2023; originally announced April 2023.

Comments: 11 pages, 8 figures

Journal ref: Phys. Rev. D 108, 012016 (2023)

First Dark Matter Search with Nuclear Recoils from the XENONnT Experiment

Authors: XENON Collaboration, E. Aprile, K. Abe, F. Agostini, S. Ahmed Maouloud, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, V. C. Antochi, D. Antón Martin, F. Arneodo, L. Baudis, A. L. Baxter, M. Bazyk, L. Bellagamba, R. Biondi, A. Bismark, E. J. Brookes, A. Brown, S. Bruenner, G. Bruno, R. Budnik, T. K. Bui, C. Cai , et al. (141 additional authors not shown)

Abstract: We report on the first search for nuclear recoils from dark matter in the form of weakly interacting massive particles (WIMPs) with the XENONnT experiment which is based on a two-phase time projection chamber with a sensitive liquid xenon mass of $5.9$ t. During the approximately 1.1 tonne-year exposure used for this search, the intrinsic $^{85}$Kr and $^{222}$Rn concentrations in the liquid targe… ▽ More We report on the first search for nuclear recoils from dark matter in the form of weakly interacting massive particles (WIMPs) with the XENONnT experiment which is based on a two-phase time projection chamber with a sensitive liquid xenon mass of $5.9$ t. During the approximately 1.1 tonne-year exposure used for this search, the intrinsic $^{85}$Kr and $^{222}$Rn concentrations in the liquid target were reduced to unprecedentedly low levels, giving an electronic recoil background rate of $(15.8\pm1.3)~\mathrm{events}/(\mathrm{t\cdot y \cdot keV})$ in the region of interest. A blind analysis of nuclear recoil events with energies between $3.3$ keV and $60.5$ keV finds no significant excess. This leads to a minimum upper limit on the spin-independent WIMP-nucleon cross section of $2.58\times 10^{-47}~\mathrm{cm}^2$ for a WIMP mass of $28~\mathrm{GeV}/c^2$ at $90\%$ confidence level. Limits for spin-dependent interactions are also provided. Both the limit and the sensitivity for the full range of WIMP masses analyzed here improve on previous results obtained with the XENON1T experiment for the same exposure. △ Less

Submitted 5 August, 2023; v1 submitted 26 March, 2023; originally announced March 2023.

Comments: Limit points are included in the submission file

Journal ref: Phys. Rev. Lett. 131, 041003 (2023)

Electron transport measurements in liquid xenon with Xenoscope, a large-scale DARWIN demonstrator

Authors: L. Baudis, Y. Biondi, A. Bismark, A. P. Cimental Chavez, J. J. Cuenca-Garcia, J. Franchi, M. Galloway, F. Girard, R. Peres, D. Ramirez Garcia, P. Sanchez-Lucas, K. Thieme, C. Wittweg

Abstract: There is a compelling physics case for a large, xenon-based underground detector devoted to dark matter and other rare-event searches. A two-phase time projection chamber as inner detector allows for a good energy resolution, a three-dimensional position determination of the interaction site and particle discrimination. To study challenges related to the construction and operation of a multi-tonne… ▽ More There is a compelling physics case for a large, xenon-based underground detector devoted to dark matter and other rare-event searches. A two-phase time projection chamber as inner detector allows for a good energy resolution, a three-dimensional position determination of the interaction site and particle discrimination. To study challenges related to the construction and operation of a multi-tonne scale detector, we have designed and constructed a vertical, full-scale demonstrator for the DARWIN experiment at the University of Zurich. Here we present first results from a several-months run with 343 kg of xenon and electron drift lifetime and transport measurements with a 53 cm tall purity monitor immersed in the cryogenic liquid. After 88 days of continuous purification, the electron lifetime reached a value of 664(23) microseconds. We measured the drift velocity of electrons for electric fields in the range (25--75) V/cm, and found values consistent with previous measurements. We also calculated the longitudinal diffusion constant of the electron cloud in the same field range, and compared with previous data, as well as with predictions from an empirical model. △ Less

Submitted 12 December, 2023; v1 submitted 24 March, 2023; originally announced March 2023.

Journal ref: Eur. Phys. J. C 83, 717 (2023)

Cosmoglobe DR1 results. I. Improved Wilkinson Microwave Anisotropy Probe maps through Bayesian end-to-end analysis

Authors: D. J. Watts, A. Basyrov, J. R. Eskilt, M. Galloway, L. T. Hergt, D. Herman, H. T. Ihle, S. Paradiso, F. Rahman, H. Thommesen, R. Aurlien, M. Bersanelli, L. A. Bianchi, M. Brilenkov, L. P. L. Colombo, H. K. Eriksen, C. Franceschet, U. Fuskeland, E. Gjerløw, B. Hensley, G. A. Hoerning, K. Lee, J. G. S. Lunde, A. Marins, S. K. Nerval , et al. (8 additional authors not shown)

Abstract: We present Cosmoglobe Data Release 1, which implements the first joint analysis of WMAP and Planck LFI time-ordered data, processed within a single Bayesian end-to-end framework. This framework builds directly on a similar analysis of the LFI measurements by the BeyondPlanck collaboration, and approaches the CMB analysis challenge through Gibbs sampling of a global posterior distribution, simultan… ▽ More We present Cosmoglobe Data Release 1, which implements the first joint analysis of WMAP and Planck LFI time-ordered data, processed within a single Bayesian end-to-end framework. This framework builds directly on a similar analysis of the LFI measurements by the BeyondPlanck collaboration, and approaches the CMB analysis challenge through Gibbs sampling of a global posterior distribution, simultaneously accounting for calibration, mapmaking, and component separation. The computational cost of producing one complete WMAP+LFI Gibbs sample is 812 CPU-hr, of which 603 CPU-hrs are spent on WMAP low-level processing; this demonstrates that end-to-end Bayesian analysis of the WMAP data is computationally feasible. We find that our WMAP posterior mean temperature sky maps and CMB temperature power spectrum are largely consistent with the official WMAP9 results. Perhaps the most notable difference is that our CMB dipole amplitude is $3366.2 \pm 1.4\ \mathrm{μK}$, which is $11\ \mathrm{μK}$ higher than the WMAP9 estimate and $2.5\ σ$ higher than BeyondPlanck; however, it is in perfect agreement with the HFI-dominated Planck PR4 result. In contrast, our WMAP polarization maps differ more notably from the WMAP9 results, and in general exhibit significantly lower large-scale residuals. We attribute this to a better constrained gain and transmission imbalance model. It is particularly noteworthy that the W-band polarization sky map, which was excluded from the official WMAP cosmological analysis, for the first time appears visually consistent with the V-band sky map. Similarly, the long standing discrepancy between the WMAP K-band and LFI 30 GHz maps is finally resolved, and the difference between the two maps appears consistent with instrumental noise at high Galactic latitudes. All maps and the associated code are made publicly available through the Cosmoglobe web page. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: 65 pages, 61 figures. Data available at cosmoglobe.uio.no. Submitted to A&A

Journal ref: A&A 679, A143 (2023)

Time-resolved Optical Polarization Monitoring of the Most Variable Brown Dwarf

Authors: Elena Manjavacas, Paulo A. Miles-Paez, Theodora Karalidi, Johanna M. Vos, Max L. Galloway, Julien H. Girard

Abstract: Recent atmospheric models for brown dwarfs suggest that the existence of clouds in substellar objects is not needed to reproduce their spectra, nor their rotationally-induced photometric variability, believed to be due to the heterogeneous cloud coverage of brown dwarf atmospheres. Cloud-free atmospheric models also predict that their flux should not be polarized, as polarization is produced by th… ▽ More Recent atmospheric models for brown dwarfs suggest that the existence of clouds in substellar objects is not needed to reproduce their spectra, nor their rotationally-induced photometric variability, believed to be due to the heterogeneous cloud coverage of brown dwarf atmospheres. Cloud-free atmospheric models also predict that their flux should not be polarized, as polarization is produced by the light-scattering of particles in the inhomogeneous cloud layers of brown dwarf atmospheres. To shed light on this dichotomy, we monitored the linear polarization and photometric variability of the most variable brown dwarf, 2MASS J21392676+0220226. We used FORS2 at the UT1 telescope to monitor the object in the z-band for six hours, split on two consecutive nights, covering one-third of its rotation period. We obtained the Stokes parameters, and we derived its time-resolved linear polarization, for which we did not find significant linear polarization (P = 0.14+\-0.07 %). We modeled the linear polarimetric signal expected assuming a map with one or two spot-like features and two bands using a polarization-enabled radiative-transfer code. We obtained values compatible with the time-resolved polarimetry obtained for 2MASS J21392676+0220226. The lack of significant polarization might be due to photometric variability produced mostly by banded structures or small-scale vortices, which cancel out the polarimetric signal from different regions of the dwarf's disk. Alternatively, the lack of clouds in 2MASS J21392676+0220226 would also explain the lack of polarization. Further linear polarimetric monitoring of 2MASS J21392676+0220226, during at least one full rotational period, would help to confirm or discard the existence of clouds in its atmosphere. △ Less

Submitted 13 March, 2023; v1 submitted 10 March, 2023; originally announced March 2023.

Comments: Accepted for publication Feb 28

Tensor-to-scalar ratio forecasts for extended LiteBIRD frequency configurations

Authors: U. Fuskeland, J. Aumont, R. Aurlien, C. Baccigalupi, A. J. Banday, H. K. Eriksen, J. Errard, R. T. Génova-Santos, T. Hasebe, J. Hubmayr, H. Imada, N. Krachmalnicoff, L. Lamagna, G. Pisano, D. Poletti, M. Remazeilles, K. L. Thompson, L. Vacher, I. K. Wehus, S. Azzoni, M. Ballardini, R. B. Barreiro, N. Bartolo, A. Basyrov, D. Beck , et al. (92 additional authors not shown)

Abstract: LiteBIRD is a planned JAXA-led CMB B-mode satellite experiment aiming for launch in the late 2020s, with a primary goal of detecting the imprint of primordial inflationary gravitational waves. Its current baseline focal-plane configuration includes 15 frequency bands between 40 and 402 GHz, fulfilling the mission requirements to detect the amplitude of gravitational waves with the total uncertaint… ▽ More LiteBIRD is a planned JAXA-led CMB B-mode satellite experiment aiming for launch in the late 2020s, with a primary goal of detecting the imprint of primordial inflationary gravitational waves. Its current baseline focal-plane configuration includes 15 frequency bands between 40 and 402 GHz, fulfilling the mission requirements to detect the amplitude of gravitational waves with the total uncertainty on the tensor-to-scalar ratio, $δr$, down to $δr<0.001$. A key aspect of this performance is accurate astrophysical component separation, and the ability to remove polarized thermal dust emission is particularly important. In this paper we note that the CMB frequency spectrum falls off nearly exponentially above 300 GHz relative to the thermal dust SED, and a relatively minor high frequency extension can therefore result in even lower uncertainties and better model reconstructions. Specifically, we compare the baseline design with five extended configurations, while varying the underlying dust modeling, in each of which the HFT (High-Frequency Telescope) frequency range is shifted logarithmically towards higher frequencies, with an upper cutoff ranging between 400 and 600 GHz. In each case, we measure the tensor-to-scalar ratio $r$ uncertainty and bias using both parametric and minimum-variance component-separation algorithms. When the thermal dust sky model includes a spatially varying spectral index and temperature, we find that the statistical uncertainty on $r$ after foreground cleaning may be reduced by as much as 30--50 % by extending the upper limit of the frequency range from 400 to 600 GHz, with most of the improvement already gained at 500 GHz. We also note that a broader frequency range leads to better ability to discriminate between models through higher $χ^2$ sensitivity. (abridged) △ Less

Submitted 15 August, 2023; v1 submitted 10 February, 2023; originally announced February 2023.

Comments: 18 pages, 13 figures. Published in A&A

Journal ref: A&A 676, A42 (2023)

The Triggerless Data Acquisition System of the XENONnT Experiment

Authors: E. Aprile, J. Aalbers, K. Abe, F. Agostini, S. Ahmed Maouloud, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, V. C. Antochi, D. Antón Martin, F. Arneodo, L. Baudis, A. L. Baxter, L. Bellagamba, R. Biondi, A. Bismark, E. J. Brookes, A. Brown, S. Bruenner, G. Bruno, R. Budnik, T. K. Bui, C. Cai, J. M. R. Cardoso , et al. (140 additional authors not shown)

Abstract: The XENONnT detector uses the latest and largest liquid xenon-based time projection chamber (TPC) operated by the XENON Collaboration, aimed at detecting Weakly Interacting Massive Particles and conducting other rare event searches. The XENONnT data acquisition (DAQ) system constitutes an upgraded and expanded version of the XENON1T DAQ system. For its operation, it relies predominantly on commerc… ▽ More The XENONnT detector uses the latest and largest liquid xenon-based time projection chamber (TPC) operated by the XENON Collaboration, aimed at detecting Weakly Interacting Massive Particles and conducting other rare event searches. The XENONnT data acquisition (DAQ) system constitutes an upgraded and expanded version of the XENON1T DAQ system. For its operation, it relies predominantly on commercially available hardware accompanied by open-source and custom-developed software. The three constituent subsystems of the XENONnT detector, the TPC (main detector), muon veto, and the newly introduced neutron veto, are integrated into a single DAQ, and can be operated both independently and as a unified system. In total, the DAQ digitizes the signals of 698 photomultiplier tubes (PMTs), of which 253 from the top PMT array of the TPC are digitized twice, at $\times10$ and $\times0.5$ gain. The DAQ for the most part is a triggerless system, reading out and storing every signal that exceeds the digitization thresholds. Custom-developed software is used to process the acquired data, making it available within $\mathcal{O}\left(10\text{ s}\right)$ for live data quality monitoring and online analyses. The entire system with all the three subsystems was successfully commissioned and has been operating continuously, comfortably withstanding readout rates that exceed $\sim500$ MB/s during calibration. Livetime during normal operation exceeds $99\%$ and is $\sim90\%$ during most high-rate calibrations. The combined DAQ system has collected more than 2 PB of both calibration and science data during the commissioning of XENONnT and the first science run. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Foreground Separation and Constraints on Primordial Gravitational Waves with the PICO Space Mission

Authors: Ragnhild Aurlien, Mathieu Remazeilles, Sebastian Belkner, Julien Carron, Jacques Delabrouille, Hans Kristian Eriksen, Raphael Flauger, Unni Fuskeland, Mathew Galloway, Krzysztof M. Gorski, Shaul Hanany, Brandon S. Hensley, J. Colin Hill, Charles R. Lawrence, Alexander van Engelen, Ingunn Kathrine Wehus

Abstract: PICO is a concept for a NASA probe-scale mission aiming to detect or constrain the tensor to scalar ratio $r$, a parameter that quantifies the amplitude of inflationary gravity waves. We carry out map-based component separation on simulations with five foreground models and input $r$ values $r_{in}=0$ and $r_{in} = 0.003$. We forecast $r$ determinations using a Gaussian likelihood assuming either… ▽ More PICO is a concept for a NASA probe-scale mission aiming to detect or constrain the tensor to scalar ratio $r$, a parameter that quantifies the amplitude of inflationary gravity waves. We carry out map-based component separation on simulations with five foreground models and input $r$ values $r_{in}=0$ and $r_{in} = 0.003$. We forecast $r$ determinations using a Gaussian likelihood assuming either no delensing or a residual lensing factor $A_{\rm lens}$ = 27%. By implementing the first full-sky, post component-separation, map-domain delensing, we show that PICO should be able to achieve $A_{\rm lens}$ = 22% - 24%. For four of the five foreground models we find that PICO would be able to set the constraints $r < 1.3 \times 10^{-4} \,\, \mbox{to} \,\, r <2.7 \times 10^{-4}\, (95\%)$ if $r_{in}=0$, the strongest constraints of any foreseeable instrument. For these models, $r=0.003$ is recovered with confidence levels between $18σ$ and $27σ$. We find weaker, and in some cases significantly biased, upper limits when removing few low or high frequency bands. The fifth model gives a $3σ$ detection when $r_{in}=0$ and a $3σ$ bias with $r_{in} = 0.003$. However, by correlating $r$ determinations from many small 2.5% sky areas with the mission's 555 GHz data we identify and mitigate the bias. This analysis underscores the importance of large sky coverage. We show that when only low multipoles $\ell \leq 12$ are used, the non-Gaussian shape of the true likelihood gives uncertainties that are on average 30% larger than a Gaussian approximation. △ Less

Submitted 16 June, 2023; v1 submitted 25 November, 2022; originally announced November 2022.

Comments: 34 pages, 13 figures, published in JCAP

Journal ref: JCAP06(2023)034

Low-energy Calibration of XENON1T with an Internal $^{37}$Ar Source

Authors: E. Aprile, K. Abe, F. Agostini, S. Ahmed Maouloud, M. Alfonsi, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, V. C. Antochi, D. Antón Martin, F. Arneodo, L. Baudis, A. L. Baxter, L. Bellagamba, R. Biondi, A. Bismark, A. Brown, S. Bruenner, G. Bruno, R. Budnik, T. K. Bui, C. Cai, C. Capelli, J. M. R. Cardoso , et al. (139 additional authors not shown)

Abstract: A low-energy electronic recoil calibration of XENON1T, a dual-phase xenon time projection chamber, with an internal $^{37}$Ar source was performed. This calibration source features a 35-day half-life and provides two mono-energetic lines at 2.82 keV and 0.27 keV. The photon yield and electron yield at 2.82 keV are measured to be (32.3$\pm$0.3) photons/keV and (40.6$\pm$0.5) electrons/keV, respecti… ▽ More A low-energy electronic recoil calibration of XENON1T, a dual-phase xenon time projection chamber, with an internal $^{37}$Ar source was performed. This calibration source features a 35-day half-life and provides two mono-energetic lines at 2.82 keV and 0.27 keV. The photon yield and electron yield at 2.82 keV are measured to be (32.3$\pm$0.3) photons/keV and (40.6$\pm$0.5) electrons/keV, respectively, in agreement with other measurements and with NEST predictions. The electron yield at 0.27 keV is also measured and it is (68.0$^{+6.3}_{-3.7}$) electrons/keV. The $^{37}$Ar calibration confirms that the detector is well-understood in the energy region close to the detection threshold, with the 2.82 keV line reconstructed at (2.83$\pm$0.02) keV, which further validates the model used to interpret the low-energy electronic recoil excess previously reported by XENON1T. The ability to efficiently remove argon with cryogenic distillation after the calibration proves that $^{37}$Ar can be considered as a regular calibration source for multi-tonne xenon detectors. △ Less

Submitted 21 March, 2023; v1 submitted 25 November, 2022; originally announced November 2022.

Weak lensing in the blue: a counter-intuitive strategy for stratospheric observations

Authors: Mohamed M. Shaaban, Ajay S. Gill, Jacqueline McCleary, Richard J. Massey, Steven J. Benton, Anthony M. Brown, Christopher J. Damaren, Tim Eifler, Aurelien A. Fraisse, Spencer Everett, Mathew N. Galloway, Michael Henderson, Bradley Holder, Eric M. Huff, Mathilde Jauzac, William C. Jones, David Lagattuta, Jason Leung, Lun Li, Thuy Vy T. Luu Johanna M. Nagy, C. Barth Netterfield, Susan F. Redmond, Jason D. Rhodes, Andrew Robertson, Jurgen Schmoll , et al. (2 additional authors not shown)

Abstract: The statistical power of weak lensing measurements is principally driven by the number of high redshift galaxies whose shapes are resolved. Conventional wisdom and physical intuition suggest this is optimised by deep imaging at long (red or near IR) wavelengths, to avoid losing redshifted Balmer break and Lyman break galaxies. We use the synthetic Emission Line EL-COSMOS catalogue to simulate lens… ▽ More The statistical power of weak lensing measurements is principally driven by the number of high redshift galaxies whose shapes are resolved. Conventional wisdom and physical intuition suggest this is optimised by deep imaging at long (red or near IR) wavelengths, to avoid losing redshifted Balmer break and Lyman break galaxies. We use the synthetic Emission Line EL-COSMOS catalogue to simulate lensing observations using different filters, from various altitudes. Here were predict the number of exposures to achieve a target z > 0.3 source density, using off-the-shelf and custom filters. Ground-based observations are easily better at red wavelengths, as (more narrowly) are space-based observations. However, we find that SuperBIT, a diffraction-limited observatory operating in the stratosphere, should instead perform its lensing-quality observations at blue wavelengths. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Effective Field Theory and Inelastic Dark Matter Results from XENON1T

Authors: E. Aprile, K. Abe, F. Agostini, S. Ahmed Maouloud, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, V. C. Antochi, D. Antón Martin, F. Arneodo, L. Baudis, A. L. Baxter, L. Bellagamba, R. Biondi, A. Bismark, A. Brown, S. Bruenner, G. Bruno, R. Budnik, C. Cai, C. Capelli, J. M. R. Cardoso, D. Cichon, M. Clark , et al. (135 additional authors not shown)

Abstract: In this work, we expand on the XENON1T nuclear recoil searches to study the individual signals of dark matter interactions from operators up to dimension-eight in a Chiral Effective Field Theory (ChEFT) and a model of inelastic dark matter (iDM). We analyze data from two science runs of the XENON1T detector totaling 1\,tonne$\times$year exposure. For these analyses, we extended the region of inter… ▽ More In this work, we expand on the XENON1T nuclear recoil searches to study the individual signals of dark matter interactions from operators up to dimension-eight in a Chiral Effective Field Theory (ChEFT) and a model of inelastic dark matter (iDM). We analyze data from two science runs of the XENON1T detector totaling 1\,tonne$\times$year exposure. For these analyses, we extended the region of interest from [4.9, 40.9]$\,$keV$_{\text{NR}}$ to [4.9, 54.4]$\,$keV$_{\text{NR}}$ to enhance our sensitivity for signals that peak at nonzero energies. We show that the data is consistent with the background-only hypothesis, with a small background over-fluctuation observed peaking between 20 and 50$\,$keV$_{\text{NR}}$, resulting in a maximum local discovery significance of 1.7\,$σ$ for the Vector$\otimes$Vector$_{\text{strange}}$ ($VV_s$) ChEFT channel for a dark matter particle of 70$\,$GeV/c$^2$, and $1.8\,σ$ for an iDM particle of 50$\,$GeV/c$^2$ with a mass splitting of 100$\,$keV/c$^2$. For each model, we report 90\,\% confidence level (CL) upper limits. We also report upper limits on three benchmark models of dark matter interaction using ChEFT where we investigate the effect of isospin-breaking interactions. We observe rate-driven cancellations in regions of the isospin-breaking couplings, leading to up to 6 orders of magnitude weaker upper limits with respect to the isospin-conserving case. △ Less

Submitted 17 October, 2022; v1 submitted 14 October, 2022; originally announced October 2022.

An approximate likelihood for nuclear recoil searches with XENON1T data

Authors: E. Aprile, K. Abe, F. Agostini, S. Ahmed Maouloud, M. Alfonsi, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, V. C. Antochi, D. Antón Martin, F. Arneodo, L. Baudis, A. L. Baxter, L. Bellagamba, R. Biondi, A. Bismark, A. Brown, S. Bruenner, G. Bruno, R. Budnik, C. Capelli, J. M. R. Cardoso, D. Cichon, B. Cimmino , et al. (129 additional authors not shown)

Abstract: The XENON collaboration has published stringent limits on specific dark matter -nucleon recoil spectra from dark matter recoiling on the liquid xenon detector target. In this paper, we present an approximate likelihood for the XENON1T 1 tonne-year nuclear recoil search applicable to any nuclear recoil spectrum. Alongside this paper, we publish data and code to compute upper limits using the method… ▽ More The XENON collaboration has published stringent limits on specific dark matter -nucleon recoil spectra from dark matter recoiling on the liquid xenon detector target. In this paper, we present an approximate likelihood for the XENON1T 1 tonne-year nuclear recoil search applicable to any nuclear recoil spectrum. Alongside this paper, we publish data and code to compute upper limits using the method we present. The approximate likelihood is constructed in bins of reconstructed energy, profiled along the signal expectation in each bin. This approach can be used to compute an approximate likelihood and therefore most statistical results for any nuclear recoil spectrum. Computing approximate results with this method is approximately three orders of magnitude faster than the likelihood used in the original publications of XENON1T, where limits were set for specific families of recoil spectra. Using this same method, we include toy Monte Carlo simulation-derived binwise likelihoods for the upcoming XENONnT experiment that can similarly be used to assess the sensitivity to arbitrary nuclear recoil signatures in its eventual 20 tonne-year exposure. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: Accepted by European Physical Journal C

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.

Journal ref: A&A 675, A4 (2023)

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.

Comments: 32 pages, 23 figures, data available from https://www.cosmoglobe.uio.no/

Journal ref: A&A 675, A10 (2023)

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.

Comments: 26 pages, 24 figures. Submitted to A&A. Part of the BeyondPlanck paper suite

Journal ref: A&A 675, A11 (2023)

Search for New Physics in Electronic Recoil Data from XENONnT

Authors: E. Aprile, K. Abe, F. Agostini, S. Ahmed Maouloud, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, V. C. Antochi, D. Antón Martin, F. Arneodo, L. Baudis, A. L. Baxter, L. Bellagamba, R. Biondi, A. Bismark, A. Brown, S. Bruenner, G. Bruno, R. Budnik, T. K. Bui, C. Cai, C. Capelli, J. M. R. Cardoso, D. Cichon , et al. (141 additional authors not shown)

Abstract: We report on a blinded analysis of low-energy electronic-recoil data from the first science run of the XENONnT dark matter experiment. Novel subsystems and the increased 5.9 tonne liquid xenon target reduced the background in the (1, 30) keV search region to $(15.8 \pm 1.3)$ events/(tonne$\times$year$\times$keV), the lowest ever achieved in a dark matter detector and $\sim$5 times lower than in XE… ▽ More We report on a blinded analysis of low-energy electronic-recoil data from the first science run of the XENONnT dark matter experiment. Novel subsystems and the increased 5.9 tonne liquid xenon target reduced the background in the (1, 30) keV search region to $(15.8 \pm 1.3)$ events/(tonne$\times$year$\times$keV), the lowest ever achieved in a dark matter detector and $\sim$5 times lower than in XENON1T. With an exposure of 1.16 tonne-years, we observe no excess above background and set stringent new limits on solar axions, an enhanced neutrino magnetic moment, and bosonic dark matter. △ Less

Submitted 15 November, 2022; v1 submitted 22 July, 2022; originally announced July 2022.

Comments: 10 pages, 6 figures. Data points are available in https://zenodo.org/record/7311940

Journal ref: Phys. Rev. Lett. 129, 161805 (2022)

Cosmoglobe: Simulating zodiacal emission with ZodiPy

Authors: M. San, D. Herman, G. B. Erikstad, M. Galloway, D. Watts

Abstract: We present ZodiPy, a modern and easy-to-use Python package for modeling the zodiacal emission seen by an arbitrary Solar System observer, which can be used for the removal of both thermal emission and scattered sunlight from interplanetary dust in astrophysical data. The code implements the COBE Diffuse Infrared Background Experiment (DIRBE) interplanetary dust model and the Planck extension, whic… ▽ More We present ZodiPy, a modern and easy-to-use Python package for modeling the zodiacal emission seen by an arbitrary Solar System observer, which can be used for the removal of both thermal emission and scattered sunlight from interplanetary dust in astrophysical data. The code implements the COBE Diffuse Infrared Background Experiment (DIRBE) interplanetary dust model and the Planck extension, which allows for zodiacal emission predictions at infrared wavelengths in the 1.25-240 $μ$m range and at microwave frequencies in the 30-857 GHz range. The predicted zodiacal emission may be extrapolated to frequencies and wavelengths not covered by the built-in models to produce forecasts for future experiments. ZodiPy attempts to enable the development of new interplanetary dust models by providing the community with an easy-to-use interface for testing both current and future models. We demonstrate how the software can be used by creating simulated zodiacal emission timestreams for the DIRBE experiment and show that these agree with corresponding timestreams produced with the DIRBE Zodiacal Light Prediction Software. We also make binned maps of the zodiacal emission as predicted to be observed by DIRBE and compare them with the DIRBE calibrated individual observations (CIO). △ Less

Submitted 30 August, 2022; v1 submitted 25 May, 2022; originally announced May 2022.

Comments: Accepted by A&A, 10 pages, 7 figures. Made edits in response to comments from the referee and the language editor. This work is part of the Cosmoglobe effort, and the code is published under an open-source license at https://github.com/Cosmoglobe/zodipy

Journal ref: A&A 666, A107 (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.

Comments: 13 pages, 10 figures

MSC Class: 83F05

Journal ref: A&A 675, A12 (2023)

Double-Weak Decays of $^{124}$Xe and $^{136}$Xe in the XENON1T and XENONnT Experiments

Authors: E. Aprile, K. Abe, F. Agostini, S. Ahmed Maouloud, M. Alfonsi, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, V. C. Antochi, D. Antón Martin, F. Arneodo, L. Baudis, A. L. Baxter, L. Bellagamba, R. Biondi, A. Bismark, A. Brown, S. Bruenner, G. Bruno, R. Budnik, C. Cai, C. Capelli, J. M. R. Cardoso, D. Cichon , et al. (135 additional authors not shown)

Abstract: We present results on the search for double-electron capture ($2ν\text{ECEC}$) of $^{124}$Xe and neutrinoless double-$β$ decay ($0νββ$) of $^{136}$Xe in XENON1T. We consider captures from the K- up to the N-shell in the $2ν\text{ECEC}$ signal model and measure a total half-life of $T_{1/2}^{2ν\text{ECEC}}=(1.1\pm0.2_\text{stat}\pm0.1_\text{sys})\times 10^{22}\;\text{yr}$ with a… ▽ More We present results on the search for double-electron capture ($2ν\text{ECEC}$) of $^{124}$Xe and neutrinoless double-$β$ decay ($0νββ$) of $^{136}$Xe in XENON1T. We consider captures from the K- up to the N-shell in the $2ν\text{ECEC}$ signal model and measure a total half-life of $T_{1/2}^{2ν\text{ECEC}}=(1.1\pm0.2_\text{stat}\pm0.1_\text{sys})\times 10^{22}\;\text{yr}$ with a $0.87\;\text{kg}\times\text{yr}$ isotope exposure. The statistical significance of the signal is $7.0\,σ$. We use XENON1T data with $36.16\;\text{kg}\times\text{yr}$ of $^{136}$Xe exposure to search for $0νββ$. We find no evidence of a signal and set a lower limit on the half-life of $T_{1/2}^{0νββ} > 1.2 \times 10^{24}\;\text{yr}\; \text{at}\; 90\,\%\;\text{CL}$. This is the best result from a dark matter detector without an enriched target to date. We also report projections on the sensitivity of XENONnT to $0νββ$. Assuming a $275\;\text{kg}\times\text{yr}$ $^{136}$Xe exposure, the expected sensitivity is $T_{1/2}^{0νββ} > 2.1 \times 10^{25}\;\text{yr}\; \text{at}\; 90\,\%\;\text{CL}$, corresponding to an effective Majorana mass range of $\langle m_{ββ} \rangle < (0.19 - 0.59)\;\text{eV/c}^2$. △ Less

Submitted 6 September, 2022; v1 submitted 9 May, 2022; originally announced May 2022.

Comments: 23 pages, 14 figures, 9 tables, version accepted for publication in Phys. Rev. C

Journal ref: Phys. Rev. C 106, 024328 (2022)

The upgraded low-background germanium counting facility Gator for high-sensitivity $γ$-ray spectrometry

Authors: G. R. Araujo, L. Baudis, Y. Biondi, A. Bismark, M. Galloway

Abstract: We describe the upgrade and performance of the high-purity germanium counting facility Gator, which is dedicated to low-background $γ$-ray spectrometry. Gator is operated at the Gran Sasso Underground Laboratory in Italy, at an average depth of 3600 meter water equivalent, and employed for material screening and selection in ultra-low background, rare-event search experiments in astroparticle phys… ▽ More We describe the upgrade and performance of the high-purity germanium counting facility Gator, which is dedicated to low-background $γ$-ray spectrometry. Gator is operated at the Gran Sasso Underground Laboratory in Italy, at an average depth of 3600 meter water equivalent, and employed for material screening and selection in ultra-low background, rare-event search experiments in astroparticle physics. The detector is equipped with a passive shield made of layers of copper, lead and polyethylene, and the sample cavity is purged with gaseous nitrogen maintained at positive pressure for radon suppression. After upgrading its enclosure, the background rate is (82.0$\pm$0.7) counts/(kg$\cdot$day) in the energy region 100 keV to 2700 keV, a 20% reduction compared to the previously reported rate. We show the stability of various operation parameters as a function of time. We also summarize the sample analysis procedure, and demonstrate Gator's sensitivity by examining one material sample, a candidate photosensor for the DARWIN experiment. △ Less

Submitted 5 July, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Journal ref: JINST 17 (2022) P08010

GPU-based optical simulation of the DARWIN detector

Authors: L. Althueser, B. Antunović, E. Aprile, D. Bajpai, L. Baudis, D. Baur, A. L. Baxter, L. Bellagamba, R. Biondi, Y. Biondi, A. Bismark, A. Brown, R. Budnik, A. Chauvin, A. P. Colijn, J. J. Cuenca-García, V. D'Andrea, P. Di Gangi, J. Dierle, S. Diglio, M. Doerenkamp, K. Eitel, S. Farrell, A. D. Ferella, C. Ferrari , et al. (55 additional authors not shown)

Abstract: Understanding propagation of scintillation light is critical for maximizing the discovery potential of next-generation liquid xenon detectors that use dual-phase time projection chamber technology. This work describes a detailed optical simulation of the DARWIN detector implemented using Chroma, a GPU-based photon tracking framework. To evaluate the framework and to explore ways of maximizing effi… ▽ More Understanding propagation of scintillation light is critical for maximizing the discovery potential of next-generation liquid xenon detectors that use dual-phase time projection chamber technology. This work describes a detailed optical simulation of the DARWIN detector implemented using Chroma, a GPU-based photon tracking framework. To evaluate the framework and to explore ways of maximizing efficiency and minimizing the time of light collection, we simulate several variations of the conventional detector design. Results of these selected studies are presented. More generally, we conclude that the approach used in this work allows one to investigate alternative designs faster and in more detail than using conventional Geant4 optical simulations, making it an attractive tool to guide the development of the ultimate liquid xenon observatory. △ Less

Submitted 11 July, 2022; v1 submitted 27 March, 2022; originally announced March 2022.

Comments: Updated to address the referees' comments, add few more authors. Journal reference added

Journal ref: JINST 17 (2022) P07018

BeyondPlanck V. Minimal ADC Corrections for Planck LFI

Authors: D. Herman, R. A. Watson, K. J. Andersen, R. Aurlien, R. Banjeri, 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, M. Iacobellis, M. Ieronymaki, H. T. Ihle, J. B. Jewell, A. Karakci, E. Keihänen , et al. (14 additional authors not shown)

Abstract: We describe the correction procedure for Analog-to-Digital Converter (ADC) differential non-linearities (DNL) adopted in the Bayesian end-to-end BeyondPlanck analysis framework. This method is nearly identical to that developed for the official LFI Data Processing Center (DPC) analysis, and relies on the binned rms noise profile of each detector data stream. However, rather than building the corre… ▽ More We describe the correction procedure for Analog-to-Digital Converter (ADC) differential non-linearities (DNL) adopted in the Bayesian end-to-end BeyondPlanck analysis framework. This method is nearly identical to that developed for the official LFI Data Processing Center (DPC) analysis, and relies on the binned rms noise profile of each detector data stream. However, rather than building the correction profile directly from the raw rms profile, we first fit a Gaussian to each significant ADC-induced rms decrement, and then derive the corresponding correction model from this smooth model. The main advange of this approach is that only samples which are significantly affected by ADC DNLs are corrected. The new corrections are only applied to data for which there is a clear detection of the non-linearities, and for which they perform at least comparably with the DPC corrections. Out of a total of 88 LFI data streams (sky and reference load for each of the 44 detectors) we apply the new minimal ADC corrections in 25 cases, and maintain the DPC corrections in 8 cases. All these correctsion are applited to 44 or 70 GHz channels, while, as in previous analyses, none of the 30 GHz ADCs show significant evidence of non-linearity. By comparing the BeyondPlanck and DPC ADC correction methods, we estimate that the residual ADC uncertainty is about two orders of magnitude below the total noise of both the 44 and 70 GHz channels, and their impact on current cosmological parameter estimation is small. However, we also show that non-idealities in the ADC corrections can generate sharp stripes in the final frequency maps, and these could be important for future joint analyses with HFI, WMAP, or other datasets. We therefore conclude that, although the existing corrections are adequate for LFI-based cosmological parameter analysis, further work on LFI ADC corrections is still warranted. △ Less

Submitted 6 April, 2022; v1 submitted 25 March, 2022; originally announced March 2022.

Comments: 9 pages, 8 figures, mild language edits

Journal ref: A&A 675, A5 (2023)

Passive low energy nuclear recoil detection with color centers -- PALEOCCENE

Authors: Krystal Alfonso, Gabriela R. Araujo, Laura Baudis, Nathaniel Bowden, Bernadette K. Cogswell, Anna Erickson, Michelle Galloway, Adam A. Hecht, Rathsara R. H. Herath Mudiyanselage, Patrick Huber, Igor Jovanovic, Giti A. Khodaparast, Brenden A. Magill, Thomas O'Donnell, Nicholas W. G. Smith, Xianyi Zhang

Abstract: The PALEOCCENE concept offers the potential for room-temperature, passive and robust detectors in the gram to kilogram range for the detection of low-energy nuclear recoil events. Nuclear recoil events can be caused by neutron scattering, coherent elastic neutrino nucleus scattering (CEvNS) or dark matter scattering and therefore, PALEOCCENE could find applications in all three areas. In this whit… ▽ More The PALEOCCENE concept offers the potential for room-temperature, passive and robust detectors in the gram to kilogram range for the detection of low-energy nuclear recoil events. Nuclear recoil events can be caused by neutron scattering, coherent elastic neutrino nucleus scattering (CEvNS) or dark matter scattering and therefore, PALEOCCENE could find applications in all three areas. In this white paper we present current and planned R&D efforts to study the feasibility of this technique. △ Less

Submitted 10 March, 2022; originally announced March 2022.

Comments: 22 pages, 6 figures, comments welcome. White paper contribution to Snowmass 2021

A Next-Generation Liquid Xenon Observatory for Dark Matter and Neutrino Physics

Authors: J. Aalbers, K. Abe, V. Aerne, F. Agostini, S. Ahmed Maouloud, D. S. Akerib, D. Yu. Akimov, J. Akshat, A. K. Al Musalhi, F. Alder, S. K. Alsum, L. Althueser, C. S. Amarasinghe, F. D. Amaro, A. Ames, T. J. Anderson, B. Andrieu, N. Angelides, E. Angelino, J. Angevaare, V. C. Antochi, D. Antón Martin, B. Antunovic, E. Aprile, H. M. Araújo , et al. (572 additional authors not shown)

Abstract: The nature of dark matter and properties of neutrinos are among the most pressing issues in contemporary particle physics. The dual-phase xenon time-projection chamber is the leading technology to cover the available parameter space for Weakly Interacting Massive Particles (WIMPs), while featuring extensive sensitivity to many alternative dark matter candidates. These detectors can also study neut… ▽ More The nature of dark matter and properties of neutrinos are among the most pressing issues in contemporary particle physics. The dual-phase xenon time-projection chamber is the leading technology to cover the available parameter space for Weakly Interacting Massive Particles (WIMPs), while featuring extensive sensitivity to many alternative dark matter candidates. These detectors can also study neutrinos through neutrinoless double-beta decay and through a variety of astrophysical sources. A next-generation xenon-based detector will therefore be a true multi-purpose observatory to significantly advance particle physics, nuclear physics, astrophysics, solar physics, and cosmology. This review article presents the science cases for such a detector. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 77 pages, 40 figures, 1262 references

Report number: INT-PUB-22-003

Journal ref: J. Phys. G: Nucl. Part. Phys. 50 (2023) 013001

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.

Comments: 11 figures, submitted to A&A. Includes updated instrument model and changes addressing referee comments

Journal ref: A&A 675, A16 (2023)

Understanding planetary context to enable life detection on exoplanets and test the Copernican principle

Authors: Joshua Krissansen-Totton, Maggie Thompson, Max L. Galloway, Jonathan J. Fortney

Abstract: The search for life on exoplanets is motivated by the universal ways in which life could modify its planetary environment. Atmospheric gases such as oxygen and methane are promising candidates for such environmental modification due to the evolutionary benefits their production would confer. However, confirming that these gases are produced by life, rather than by geochemical or astrophysical proc… ▽ More The search for life on exoplanets is motivated by the universal ways in which life could modify its planetary environment. Atmospheric gases such as oxygen and methane are promising candidates for such environmental modification due to the evolutionary benefits their production would confer. However, confirming that these gases are produced by life, rather than by geochemical or astrophysical processes, will require a thorough understanding of planetary context, including the expected counterfactual atmospheric evolution for lifeless planets. Here, we evaluate current understanding of planetary context for several candidate biosignatures and their upcoming observability. We review the contextual framework for oxygen and describe how conjectured abiotic oxygen scenarios may be testable. In contrast to oxygen, current understanding of how planetary context controls non-biological methane (CH$_4$) production is limited, even though CH$_4$ biosignatures in anoxic atmospheres may be readily detectable with the James Webb Space Telescope. We assess environmental context for CH$_4$ biosignatures and conclude that abundant atmospheric CH$_4$ coexisting with CO$_2$, and CO:CH$_4$ << 1 is suggestive of biological production, although precise thresholds are dependent on stellar context and sparsely characterized abiotic CH$_4$ scenarios. A planetary context framework is also considered for alternative or agnostic biosignatures. Whatever the distribution of life in the Universe, observations of terrestrial exoplanets in coming decades will provide a quantitative understanding of the atmospheric evolution of lifeless worlds. This knowledge will inform future instrument requirements to either corroborate the presence of life elsewhere or confirm its apparent absence. △ Less

Submitted 1 March, 2022; v1 submitted 21 February, 2022; originally announced February 2022.

Comments: 16 pages, 5 figures, updated with minor formatting corrections

Journal ref: Nat Astron 6, 189-198 (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.

Comments: 155 pages, accepted for publication in PTEP

Polarization angle requirements for CMB B-mode experiments. Application to the LiteBIRD satellite

Authors: P. Vielva, E. Martínez-González, F. J. Casas, T. Matsumura, S. Henrot-Versillé, E. Komatsu, J. Aumont, R. Aurlien, C. Baccigalupi, A. J. Banday, R. B. Barreiro, N. Bartolo, E. Calabrese, K. Cheung, F. Columbro, A. Coppolecchia, P. de Bernardis, T. de Haan, E. de la Hoz, M. De Petris, S. Della Torre, P. Diego-Palazuelos, H. K. Eriksen, J. Errard, F. Finelli , et al. (46 additional authors not shown)

Abstract: A methodology to provide the polarization angle requirements for different sets of detectors, at a given frequency of a CMB polarization experiment, is presented. The uncertainties in the polarization angle of each detector set are related to a given bias on the tensor-to-scalar ratio $r$ parameter. The approach is grounded in using a linear combination of the detector sets to obtain the CMB polar… ▽ More A methodology to provide the polarization angle requirements for different sets of detectors, at a given frequency of a CMB polarization experiment, is presented. The uncertainties in the polarization angle of each detector set are related to a given bias on the tensor-to-scalar ratio $r$ parameter. The approach is grounded in using a linear combination of the detector sets to obtain the CMB polarization signal. In addition, assuming that the uncertainties on the polarization angle are in the small angle limit (lower than a few degrees), it is possible to derive analytic expressions to establish the requirements. The methodology also accounts for possible correlations among detectors, that may originate from the optics, wafers, etc. The approach is applied to the LiteBIRD space mission. We show that, for the most restrictive case (i.e., full correlation of the polarization angle systematics among detector sets), the requirements on the polarization angle uncertainties are of around 1 arcmin at the most sensitive frequency bands (i.e., $\approx 150$ GHz) and of few tens of arcmin at the lowest (i.e., $\approx 40$ GHz) and highest (i.e., $\approx 400$ GHz) observational bands. Conversely, for the least restrictive case (i.e., no correlation of the polarization angle systematics among detector sets), the requirements are $\approx 5$ times less restrictive than for the previous scenario. At the global and the telescope levels, polarization angle knowledge of a few arcmins is sufficient for correlated global systematic errors and can be relaxed by a factor of two for fully uncorrelated errors in detector polarization angle. The reported uncertainty levels are needed in order to have the bias on $r$ due to systematics below the limit established by the LiteBIRD collaboration. △ Less

Submitted 18 April, 2022; v1 submitted 2 February, 2022; originally announced February 2022.

Comments: 26 pages, 6 figures. Minor changes to match the final version in JCAP

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.

Comments: 27 pages, 23 figures, 3 tables, part 13 of the BeyondPlanck release. All BeyondPlanck products and software will be released publicly at http://beyondplanck.science. Paper version as accepted by A&A

Journal ref: A&A 675, A13 (2023)

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.

Comments: 10 pages, 9 figures, part 15 of the BeyondPlanck release. All BeyondPlanck products and software will be released publicly at http://beyondplanck.science. Refereed, edited, and accepted by A&A. Added references, added arxiv IDs to BeyondPlanck papers. Figure update. Expanded discussion on AME models. Moved title from BeyondPlanck XV to BeyondPlanck XVI due to internal reordering

Journal ref: A&A 675, A15 (2023)

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.

Comments: 16 Pages, 7 Figures. Part of the BeyondPlanck paper suite

Journal ref: A&A 675, A3 (2023)

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.

Comments: 9 pages, 8 figures. Part of the BeyondPlanck paper suite

Journal ref: A&A 675, A8 (2023)