Showing 1–50 of 79 results for author: Jewell, J

Efficient ultra-broadband low-resolution astrophotonic spectrographs

Authors: Pradip Gatkine, Greg Sercel, Nemanja Jovanovic, Ronald Broeke, Katarzyna Lawniczuk, Marco Passoni, Ashok Balakrishnan, Serge Bidnyk, Jielong Yin, Jeffrey Jewell, J. Kent Wallace, Dimitri Mawet

Abstract: Broadband low-resolution near-infrared spectrographs in a compact form are crucial for ground- and space-based astronomy and other fields of sensing. Astronomical spectroscopy poses stringent requirements including high efficiency, broad band operation ($>$ 300 nm), and in some cases, polarization insensitivity. We present and compare experimental results from the design, fabrication, and characte… ▽ More Broadband low-resolution near-infrared spectrographs in a compact form are crucial for ground- and space-based astronomy and other fields of sensing. Astronomical spectroscopy poses stringent requirements including high efficiency, broad band operation ($>$ 300 nm), and in some cases, polarization insensitivity. We present and compare experimental results from the design, fabrication, and characterization of broadband (1200 - 1650 nm) arrayed waveguide grating (AWG) spectrographs built using the two most promising low-loss platforms - Si$_3$N$_4$ (rectangular waveguides) and doped-SiO$_2$ (square waveguides). These AWGs have a resolving power ($λ/Δλ$) of ~200, a free spectral range of ~ 200-350 nm, and a small footprint of ~ 50-100 mm$^2$. The peak overall (fiber-chip-fiber) efficiency of the doped-SiO$_2$ AWG was ~ 79\% (1 dB), and it exhibited a negligible polarization-dependent shift compared to the channel spacing. For Si$_3$N$_4$ AWGs, the peak overall efficiency in TE mode was ~ 50\% (3 dB), and the main loss component was found to be fiber-to-chip coupling losses. These broadband AWGs are key to enabling compact integrations such as multi-object spectrographs or dispersion back-ends for other astrophotonic devices such as photonic lanterns or nulling interferometers. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: 16 pages, 8 Figures, accepted for publication in Optics Express

Thermal control of long delay lines in a high-resolution astrophotonic spectrograph

Authors: Gregory P. Sercel, Pradip R. Gatkine, Nemanja Jovanovic, Jeffrey B. Jewell, Luis Pereira da Costa, J. Kent Wallace, Dimitri P. Mawet

Abstract: High-resolution astronomical spectroscopy carried out with a photonic Fourier transform spectrograph (FTS) requires long asymmetrical optical delay lines that can be dynamically tuned. For example, to achieve a spectral resolution of R = 30,000, a delay line as long as 1.5 cm would be required. Such delays are inherently prone to phase errors caused by temperature fluctuations. This is due to the… ▽ More High-resolution astronomical spectroscopy carried out with a photonic Fourier transform spectrograph (FTS) requires long asymmetrical optical delay lines that can be dynamically tuned. For example, to achieve a spectral resolution of R = 30,000, a delay line as long as 1.5 cm would be required. Such delays are inherently prone to phase errors caused by temperature fluctuations. This is due to the relatively large thermo-optic coefficient and long lengths of the waveguides, in this case composed of SiN, resulting in thermally dependent changes to the optical path length. To minimize phase error to the order of 0.05 radians, thermal stability of the order of 0.05° C is necessary. A thermal control system capable of stability such as this would require a fast thermal response and minimal overshoot/undershoot. With a PID temperature control loop driven by a Peltier cooler and thermistor, we minimized interference fringe phase error to +/- 0.025 radians and achieved temperature stability on the order of 0.05° C. We present a practical system for precision temperature control of a foundry-fabricated and packaged FTS device on a SiN platform with delay lines ranging from 0.5 to 1.5 cm in length using inexpensive off-the-shelf components, including design details, control loop optimization, and considerations for thermal control of integrated photonics. △ Less

Submitted 2 April, 2024; v1 submitted 15 March, 2024; originally announced March 2024.

2023 Astrophotonics Roadmap: pathways to realizing multi-functional integrated astrophotonic instruments

Authors: Nemanja Jovanovic, Pradip Gatkine, Narsireddy Anugu, Rodrigo Amezcua-Correa, Ritoban Basu Thakur, Charles Beichman, Chad Bender, Jean-Philippe Berger, Azzurra Bigioli, Joss Bland-Hawthorn, Guillaume Bourdarot, Charles M. Bradford, Ronald Broeke, Julia Bryant, Kevin Bundy, Ross Cheriton, Nick Cvetojevic, Momen Diab, Scott A. Diddams, Aline N. Dinkelaker, Jeroen Duis, Stephen Eikenberry, Simon Ellis, Akira Endo, Donald F. Figer , et al. (55 additional authors not shown)

Abstract: Photonics offer numerous functionalities that can be used to realize astrophotonic instruments. The most spectacular example to date is the ESO Gravity instrument at the Very Large Telescope in Chile. Integrated astrophotonic devices stand to offer critical advantages for instrument development, including extreme miniaturization, as well as integration, superior thermal and mechanical stabilizatio… ▽ More Photonics offer numerous functionalities that can be used to realize astrophotonic instruments. The most spectacular example to date is the ESO Gravity instrument at the Very Large Telescope in Chile. Integrated astrophotonic devices stand to offer critical advantages for instrument development, including extreme miniaturization, as well as integration, superior thermal and mechanical stabilization owing to the small footprint, and high replicability offering cost savings. Numerous astrophotonic technologies have been developed to address shortcomings of conventional instruments to date, including for example the development of photonic lanterns, complex aperiodic fiber Bragg gratings, complex beam combiners to enable long baseline interferometry, and laser frequency combs for high precision spectral calibration of spectrometers. Despite these successes, the facility implementation of photonic solutions in astronomical instrumentation is currently limited because of (1) low throughputs from coupling to fibers, coupling fibers to chips, propagation and bend losses, device losses, etc, (2) difficulties with scaling to large channel count devices needed for large bandwidths and high resolutions, and (3) efficient integration of photonics with detectors, to name a few. In this roadmap, we identify 24 areas that need further development. We outline the challenges and advances needed across those areas covering design tools, simulation capabilities, fabrication processes, the need for entirely new components, integration and hybridization and the characterization of devices. To realize these advances the astrophotonics community will have to work cooperatively with industrial partners who have more advanced manufacturing capabilities. With the advances described herein, multi-functional instruments will be realized leading to novel observing capabilities for both ground and space platforms. △ Less

Submitted 1 November, 2023; originally announced November 2023.

Comments: 191 pages, 47 figures. This is the version of the article before peer review or editing, as submitted by an author to J. Phys. Photonics. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at https://iopscience.iop.org/article/10.1088/2515-7647/ace869/meta

Journal ref: J. Phys. Photonics 5 042501 (2023)

The path to detecting extraterrestrial life with astrophotonics

Authors: Nemanja Jovanovic, Yinzi Xin, Michael P. Fitzgerald, Olivier Guyon, Peter Tuthill, Barnaby Norris, Pradip Gatkine, Greg Sercel, Svarun Soda, Yoo Jung Kim, Jonathan Lin, Sergio Leon-Saval, Rodrigo Amezcua-Correa, Stephanos Yerolatsitis, Julien Lozi, Sebastien Vievard, Chris Betters, Steph Sallum, Daniel Levinstein, Dimitri Mawet, Jeffrey Jewell, J. Kent Wallace, Nick Cvetojevic

Abstract: Astrophysical research into exoplanets has delivered thousands of confirmed planets orbiting distant stars. These planets span a wide ranges of size and composition, with diversity also being the hallmark of system configurations, the great majority of which do not resemble our own solar system. Unfortunately, only a handful of the known planets have been characterized spectroscopically thus far,… ▽ More Astrophysical research into exoplanets has delivered thousands of confirmed planets orbiting distant stars. These planets span a wide ranges of size and composition, with diversity also being the hallmark of system configurations, the great majority of which do not resemble our own solar system. Unfortunately, only a handful of the known planets have been characterized spectroscopically thus far, leaving a gaping void in our understanding of planetary formation processes and planetary types. To make progress, astronomers studying exoplanets will need new and innovative technical solutions. Astrophotonics -- an emerging field focused on the application of photonic technologies to observational astronomy -- provides one promising avenue forward. In this paper we discuss various astrophotonic technologies that could aid in the detection and subsequent characterization of planets and in particular themes leading towards the detection of extraterrestrial life. △ Less

Submitted 15 September, 2023; originally announced September 2023.

Comments: 9 pages, 2 figures, SPIE Optics and Photonics conference

Report number: 12680-17

Dual Purpose Lyot Coronagraph Masks for Simultaneous High-Contrast Imaging and High-Resolution Wavefront Sensing

Authors: Garreth Ruane, J. Kent Wallace, A J Eldorado Riggs, Tobias Wenger, Mahmood Bagheri, Jeffrey Jewell, Nasrat Raouf, Gregory Allan, Camilo Mejia Prada, Matthew Noyes, Alex B. Walter

Abstract: Directly imaging Earth-sized exoplanets with a visible-light coronagraph instrument on a space telescope will require a system that can achieve $\sim10^{-10}$ raw contrast and maintain it for the duration of observations (on the order of hours or more). We are designing, manufacturing, and testing Dual Purpose Lyot coronagraph (DPLC) masks that allow for simultaneous wavefront sensing and control… ▽ More Directly imaging Earth-sized exoplanets with a visible-light coronagraph instrument on a space telescope will require a system that can achieve $\sim10^{-10}$ raw contrast and maintain it for the duration of observations (on the order of hours or more). We are designing, manufacturing, and testing Dual Purpose Lyot coronagraph (DPLC) masks that allow for simultaneous wavefront sensing and control using out-of-band light to maintain high contrast in the science focal plane. Our initial design uses a tiered metallic focal plane occulter to suppress starlight in the transmitted coronagraph channel and a dichroic-coated substrate to reflect out-of-band light to a wavefront sensing camera. The occulter design introduces a phase shift such that the reflected channel is a Zernike wavefront sensor. The dichroic coating allows higher-order wavefront errors to be detected which is especially critical for compensating for residual drifts from an actively-controlled segmented primary mirror. A second-generation design concept includes a metasurface to create polarization-dependent phase shifts in the reflected beam, which has several advantages including an extended dynamic range. We will present the focal plane mask designs, characterization, and initial testing at NASA's High Contrast Imaging Testbed (HCIT) facility. △ Less

Submitted 20 August, 2023; originally announced August 2023.

Comments: To appear in the Proceedings of the SPIE, Techniques and Instrumentation for Detection of Exoplanets XI

Flattening laser frequency comb spectra with a high dynamic range, broadband spectral shaper on-a-chip

Authors: Nemanja Jovanovic, Pradip Gatkine, Boqiang Shen, Maodong Gao, Nick Cvetojevic, Katarzyna Lawniczuk, Ronald Broeke, Charles Beichman, Stephanie Leifer, Jeffery Jewell, Gautam Vasisht, Dimitri Mawet

Abstract: Spectral shaping is critical to many fields of science. In astronomy for example, the detection of exoplanets via the Doppler effect hinges on the ability to calibrate a high resolution spectrograph. Laser frequency combs can be used for this, but the wildly varying intensity across the spectrum can make it impossible to optimally utilize the entire comb, leading to a reduced overall precision of… ▽ More Spectral shaping is critical to many fields of science. In astronomy for example, the detection of exoplanets via the Doppler effect hinges on the ability to calibrate a high resolution spectrograph. Laser frequency combs can be used for this, but the wildly varying intensity across the spectrum can make it impossible to optimally utilize the entire comb, leading to a reduced overall precision of calibration. To circumvent this, astronomical applications of laser frequency combs rely on a bulk optic setup which can flatten the output spectrum before sending it to the spectrograph. Such flatteners require complex and expensive optical elements like spatial light modulators and have non-negligible bench top footprints. Here we present an alternative in the form of an all-photonic spectral shaper that can be used to flatten the spectrum of a laser frequency comb. The device consists of a circuit etched into a silicon nitride wafer that supports an arrayed-waveguide grating to disperse the light over hundreds of nanometers in wavelength, followed by Mach-Zehnder interferometers to control the amplitude of each channel, thermo-optic phase modulators to phase the channels and a second arrayed-waveguide grating to recombine the spectrum. The demonstrator device operates from 1400 to 1800 nm (covering the astronomical H band), with twenty 20 nm wide channels. The device allows for nearly 40 dBs of dynamic modulation of the spectrum via the Mach-Zehnders , which is greater than that offered by most spatial light modulators. With a superluminescent diode, we reduced the static spectral variation to ~3 dB, limited by the properties of the components used in the circuit and on a laser frequency comb we managed to reduce the modulation to 5 dBs, sufficient for astronomical applications. △ Less

Submitted 3 October, 2022; originally announced October 2022.

Comments: 15 pages, 10 figures. arXiv admin note: substantial text overlap with arXiv:2209.09455

Journal ref: Opt. Express 30, 36745-36760 (2022)

An all-photonic, dynamic device for flattening the spectrum of a laser frequency comb for precise calibration of radial velocity measurements

Authors: Nemanja Jovanovic, Pradip Gatkine, Boqiang Shen, Maodong Gao, Nick Cvetojevic, Katarzyna Ławniczuk, Ronald Broeke, Charles Beichman, Stephanie Leifer, Jeffery Jewell, Gautam Vasisht, Dimitri Mawet

Abstract: Laser frequency combs are fast becoming critical to reaching the highest radial velocity precisions. One shortcoming is the highly variable brightness of the comb lines across the spectrum (up to 4-5 orders of magnitude). This can result in some lines saturating while others are at low signal and lost in the noise. Losing lines to either of these effects reduces the precision and hence effectivene… ▽ More Laser frequency combs are fast becoming critical to reaching the highest radial velocity precisions. One shortcoming is the highly variable brightness of the comb lines across the spectrum (up to 4-5 orders of magnitude). This can result in some lines saturating while others are at low signal and lost in the noise. Losing lines to either of these effects reduces the precision and hence effectiveness of the comb. In addition, the brightness of the comb lines can vary with time which could drive comb lines with initially reasonable SNR's into the two regimes described above. To mitigate these two effects, laser frequency combs use optical flattener's. Flattener's are typically bulk optic setups that disperse the comb light with a grating, and then use a spatial light modulator to control the amplitude across the spectrum before recombining the light into another single mode fiber and sending it to the spectrograph. These setups can be large (small bench top), expensive (several hundred thousand dollars) and have limited stability. To address these issues, we have developed an all-photonic spectrum flattener on a chip. The device is constructed from optical waveguides on a SiN chip. The light from the laser frequency comb's output optical fiber can be directly connected to the chip, where the light is first dispersed using an arrayed waveguide grating. To control the brightness of each channel, the light is passed through a Mach-Zehnder interferometer before being recombined with a second arrayed waveguide grating. Thermo-optic phase modulators are used in each channel before recombination to path length match the channels as needed. Here we present the results from our first generation prototype. The device operates from 1400-1800 nm (covering the H band), with 20, 20 nm wide channels. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: 7 pages, 5 figures, conference

Journal ref: Proc. SPIE 12188, Advances in Optical and Mechanical Technologies for Telescopes and Instrumentation V, 121885D (29 August 2022);

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)

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)

Axion Dark Matter

Authors: C. B. Adams, N. Aggarwal, A. Agrawal, R. Balafendiev, C. Bartram, M. Baryakhtar, H. Bekker, P. Belov, K. K. Berggren, A. Berlin, C. Boutan, D. Bowring, D. Budker, A. Caldwell, P. Carenza, G. Carosi, R. Cervantes, S. S. Chakrabarty, S. Chaudhuri, T. Y. Chen, S. Cheong, A. Chou, R. T. Co, J. Conrad, D. Croon , et al. (130 additional authors not shown)

Abstract: Axions are well-motivated dark matter candidates with simple cosmological production mechanisms. They were originally introduced to solve the strong CP problem, but also arise in a wide range of extensions to the Standard Model. This Snowmass white paper summarizes axion phenomenology and outlines next-generation laboratory experiments proposed to detect axion dark matter. There are vibrant synerg… ▽ More Axions are well-motivated dark matter candidates with simple cosmological production mechanisms. They were originally introduced to solve the strong CP problem, but also arise in a wide range of extensions to the Standard Model. This Snowmass white paper summarizes axion phenomenology and outlines next-generation laboratory experiments proposed to detect axion dark matter. There are vibrant synergies with astrophysical searches and advances in instrumentation including quantum-enabled readout, high-Q resonators and cavities and large high-field magnets. This white paper outlines a clear roadmap to discovery, and shows that the US is well-positioned to be at the forefront of the search for axion dark matter in the coming decade. △ Less

Submitted 29 March, 2023; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: restore and expand author list

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)

Simulating the Study of Exoplanets Using Photonic Spectrographs

Authors: Marcos Perez, Pradip Gatkine, Nemanja Jovanovic, Jeffrey Jewell, J. Kent Wallace, Dimitri Mawet

Abstract: Photonic spectrographs offer a highly miniaturized, flexible, and stable on-chip solution for astronomical spectroscopy and can be used for various science cases such as determining the atmospheric composition of exoplanets to understand their habitability, formation, and evolution. Arrayed Waveguide Gratings (AWGs) have shown the best promise to be used as an astrophotonic spectrograph. We develo… ▽ More Photonic spectrographs offer a highly miniaturized, flexible, and stable on-chip solution for astronomical spectroscopy and can be used for various science cases such as determining the atmospheric composition of exoplanets to understand their habitability, formation, and evolution. Arrayed Waveguide Gratings (AWGs) have shown the best promise to be used as an astrophotonic spectrograph. We developed a publically-available tool to conduct a preliminary examination of the capability of the AWGs in spectrally resolving exoplanet atmospheres. We derived the Line-Spread-Function (LSF) as a function of wavelength and the Full-Width-at-Half-Maximum (FWHM) of the LSF as a function of spectral line width to evaluate the response of a discretely- and continuously-sampled low-resolution AWG (R $\sim$ 1000). We observed that the LSF has minimal wavelength dependence ($\sim$5\%), irrespective of the offset with respect to the center-wavelengths of the AWG channels, contrary to the previous assumptions. We further confirmed that the observed FWHM scales linearly with the emission line width. Finally, we present simulated extraction of a sample molecular absorption spectrum with the discretely- and continuously-sampled low-resolution AWGs. From this, we show that while the discrete AWG matches its expected resolving power, the continuous AWG spectrograph can, in principle, achieve an effective resolution significantly greater ($\sim$ 2x) than the discrete AWG. This detailed examination of the AWGs will be foundational for future deployment of AWG spectrographs for astronomical science cases such as exoplanet atmospheres. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: 13 pages, 6 figures, Presented at and published in the proceedings of SPIE Photonics West 2022

Journal ref: Proc. SPIE 12008, Photonic Instrumentation Engineering IX, 120080C (Presented at SPIE OPTO: 5 March 2022; Published: 5 March 2022)

An on-chip astrophotonic spectrograph with a resolving power of 12,000

Authors: Pradip Gatkine, Nemanja Jovanovic, Jeffrey Jewell, J. Kent Wallace, Dimitri Mawet

Abstract: With the upcoming extremely large telescopes (ELTs), the volume, mass, and cost of the associated spectrographs will scale with the telescope diameter. Astrophotonics offers a unique solution to this problem in the form of single-mode fiber-fed diffraction-limited spectrographs on a chip. These highly miniaturized chips offer great flexibility in terms of coherent manipulation of photons. Such pho… ▽ More With the upcoming extremely large telescopes (ELTs), the volume, mass, and cost of the associated spectrographs will scale with the telescope diameter. Astrophotonics offers a unique solution to this problem in the form of single-mode fiber-fed diffraction-limited spectrographs on a chip. These highly miniaturized chips offer great flexibility in terms of coherent manipulation of photons. Such photonic spectrographs are well-suited to disperse the light from directly imaged planets (post-coronagraph, collected using a single-mode fiber) to characterize exoplanet atmospheres. Here we present the results from a proof-of-concept high-resolution astrophotonic spectrograph using the arrayed waveguide gratings (AWG) architecture. This chip uses the low-loss SiN platform (SiN core, SiO$_2$ cladding) with square waveguides (800 nm $\times$ 800 nm). The AWG has a measured resolving power ($λ/δλ$) of $\sim$ 12,000 and a free spectral range (FSR) of 2.8 nm. While the FSR is small, the chip operates over a broad band (1200 $-$ 1700 nm). The peak on-chip throughput (excluding the coupling efficiency) is $\sim$40\% (- 4 dB) and the overall throughput (including the coupling loss) is $\sim$ 11\% (- 9.6 dB) in the TE mode. Thanks to the high-confinement waveguide geometry, the chip is highly miniaturized with a size of only 7.4 mm $\times$ 2 mm. This demonstration highlights the utility of SiN platform for astrophotonics, particularly, the capability of commercial SiN foundries to fabricate ultra-small, high-resolution, high-throughput AWG spectrographs on a chip suitable for both ground- and space-based telescopes. △ Less

Submitted 9 March, 2022; originally announced March 2022.

Comments: 10 pages, 7 figures, Published in the proceedings of SPIE Optics + Photonics conference 2021

Journal ref: Proc. SPIE 11819, UV/Optical/IR Space Telescopes and Instruments: Innovative Technologies and Concepts X, 118190I (20 August 2021)

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)

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)

BeyondPlanck X. Bandpass and beam leakage corrections

Authors: T. L. Svalheim, K. J. Andersen, R. Aurlien, R. Banerji, M. Bersanelli, S. Bertocco, M. Brilenkov, M. Carbone, L. P. L. Colombo, H. K. Eriksen, M. K. Foss, C. Franceschet, U. Fuskeland, S. Galeotta, M. Galloway, S. Gerakakis, E. Gjerløw, B. Hensley, D. Herman, M. Iacobellis, M. Ieronymaki, H. T. Ihle, J. B. Jewell, A. Karakci, E. Keihänen , et al. (14 additional authors not shown)

Abstract: We discuss the treatment of bandpass and beam leakage corrections in the Bayesian BeyondPlanck CMB analysis pipeline as applied to the Planck LFI measurements. As a preparatory step, we first apply three corrections to the nominal LFI bandpass profiles including removal of a known systematic effect in the ground measuring equipment at 61 GHz; smoothing of standing wave ripples; and edge regulariza… ▽ More We discuss the treatment of bandpass and beam leakage corrections in the Bayesian BeyondPlanck CMB analysis pipeline as applied to the Planck LFI measurements. As a preparatory step, we first apply three corrections to the nominal LFI bandpass profiles including removal of a known systematic effect in the ground measuring equipment at 61 GHz; smoothing of standing wave ripples; and edge regularization. The main net impact of these modifications is an overall shift in the 70 GHz bandpass of +0.6 GHz; we argue that any analysis of LFI data products, either from Planck or BeyondPlanck, should use these new bandpasses. In addition, we fit a single free bandpass parameter for each radiometer of the form $Δ_i = Δ_0 + δ_i$, where $Δ_0$ represents an absolute frequency shift per frequency band and $δ_i$ is a relative shift per detector. The absolute correction is only fitted at 30 GHz with a full $χ^2$-based likelihood, resulting in a correction of $Δ_{30}=0.24\pm0.03\,$GHz. The relative corrections are fitted using a spurious map approach, fundamentally similar to the method pioneered by the WMAP team, but without introducing many additional degrees of freedom. All bandpass parameters are sampled using a standard Metropolis sampler within the main BeyondPlanck Gibbs chain, and bandpass uncertainties are thus propagated to all other data products in the analysis. In total, we find that our bandpass model significantly reduces leakage effects. For beam leakage corrections, we adopt the official Planck LFI beam estimates without additional degrees of freedom, and only marginalize over the underlying sky model. We note that this is the first time leakage from beam mismatch has been included for Planck LFI maps. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: 14 pages, 11 figures, this is part 10 in the BeyondPlanck release. All BeyondPlanck products and software will be released publicly at http://beyondplanck.science. Submitted to A&A

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

Potential of commercial SiN MPW platforms for developing mid/high-resolution integrated photonic spectrographs for astronomy

Authors: Pradip Gatkine, Nemanja Jovanovic, Christopher Hopgood, Simon Ellis, Ronald Broeke, Katarzyna Ławniczuk, Jeffrey Jewell, J. Kent Wallace, Dimitri Mawet

Abstract: Integrated photonic spectrographs offer an avenue to extreme miniaturization of astronomical instruments, which would greatly benefit extremely large telescopes and future space missions. These devices first require optimization for astronomical applications, which includes design, fabrication and field-testing. Given the high costs of photonic fabrication, Multi-Project Wafer (MPW) SiN offerings,… ▽ More Integrated photonic spectrographs offer an avenue to extreme miniaturization of astronomical instruments, which would greatly benefit extremely large telescopes and future space missions. These devices first require optimization for astronomical applications, which includes design, fabrication and field-testing. Given the high costs of photonic fabrication, Multi-Project Wafer (MPW) SiN offerings, where a user purchases a portion of a wafer, provide a convenient and affordable avenue to develop this technology. In this work we study the potential of two commonly used SiN waveguide geometries by MPW foundries, i.e. square and rectangular profiles to determine how they affect the performance of mid-high resolution arrayed waveguide grating spectrometers around 1.5 $μ$m. Specifically, we present results from detailed simulations on the mode sizes, shapes, and polarization properties, and on the impact of phase errors on the throughput and cross talk as well as some laboratory results of coupling and propagation losses. From the MPW-run tolerances and our phase-error study, we estimate that an AWG with R $\sim$ 10,000 can be developed with the MPW runs and even greater resolving power is achievable with more reliable, dedicated fabrication runs. Depending on the fabrication and design optimizations, it is possible to achieve throughputs $\sim 60\%$ using the SiN platform. Thus, we show that SiN MPW offerings are highly promising and will play a key role in integrated photonic spectrograph developments for astronomy. △ Less

Submitted 8 June, 2021; originally announced June 2021.

Comments: Accepted at Applied Optics, Special Issue on Astrophotonics (24 pages, 13 figures, 3 tables)

Journal ref: Applied Optics Vol. 60, Issue 19, pp. D15-D32 (2021)

BeyondPlanck XV. Polarized foreground emission between 30 and 70 GHz

Authors: T. L. Svalheim, K. J. Andersen, R. Aurlien, R. Banerji, M. Bersanelli, S. Bertocco, M. Brilenkov, M. Carbone, L. P. L. Colombo, H. K. Eriksen, M. K. Foss, C. Franceschet, U. Fuskeland, S. Galeotta, M. Galloway, S. Gerakakis, E. Gjerløw, B. Hensley, D. Herman, M. Iacobellis, M. Ieronymaki, H. T. Ihle, J. B. Jewell, A. Karakci, E. Keihänen , et al. (13 additional authors not shown)

Abstract: We constrain polarized foreground emission between 30 and 70 GHz with the Planck Low Frequency Instrument (LFI) and WMAP data within the global Bayesian BeyondPlanck framework. We combine for the first time full-resolution Planck LFI time-ordered data with low-resolution WMAP sky maps at 33, 40 and 61 GHz. Spectral parameters are fit with a likelihood defined at the native resolution of each frequ… ▽ More We constrain polarized foreground emission between 30 and 70 GHz with the Planck Low Frequency Instrument (LFI) and WMAP data within the global Bayesian BeyondPlanck framework. We combine for the first time full-resolution Planck LFI time-ordered data with low-resolution WMAP sky maps at 33, 40 and 61 GHz. Spectral parameters are fit with a likelihood defined at the native resolution of each frequency channel. This analysis represents the first implementation of true multi-resolution component separation applied to CMB observations for both amplitude and spectral energy distribution (SED) parameters. For synchrotron emission, we approximate the SED as a power-law in frequency and find that the low signal-to-noise ratio of the current data strongly limits the number of free parameters that may be robustly constrained. We partition the sky into four large disjoint regions (High Latitude; Galactic Spur; Galactic Plane; and Galactic Center), each associated with its own power-law index. We find that the High Latitude region is prior-dominated, while the Galactic Center region is contaminated by residual instrumental systematics. The two remaining regions appear to be signal-dominated, and for these we derive spectral indices of $β_{\mathrm s}^{\mathrm{Spur}}=-3.17\pm0.06$ and $β_{\mathrm s}^{\mathrm{Plane}}=-3.03\pm0.07$, in good agreement with previous results. For thermal dust emission we assume a modified blackbody model and we fit a single power-law index across the full sky. We find $β_{\mathrm{d}}=1.64\pm0.03$, which is slightly steeper than reported from Planck HFI data, but still statistically consistent at the 2$σ$ confidence level. △ Less

Submitted 18 January, 2022; v1 submitted 17 November, 2020; originally announced November 2020.

Comments: 18 pages, 15 figures. All BeyondPlanck products and software will be released publicly at http://beyondplanck.science. This paper has been submitted to A&A. Major revision, matches submitted version

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

BeyondPlanck VII. Bayesian estimation of gain and absolute calibration for CMB experiments

Authors: E. Gjerløw, H. T. Ihle, S. Galeotta, K. J. Andersen, R. Aurlien, R. Banerji, M. Bersanelli, S. Bertocco, M. Brilenkov, M. Carbone, L. P. L. Colombo, H. K. Eriksen, M. K. Foss, C. Franceschet, U. Fuskeland, M. Galloway, S. Gerakakis, B. Hensley, D. Herman, M. Iacobellis, M. Ieronymaki, J. B. Jewell, A. Karakci, E. Keihänen, R. Keskitalo , et al. (13 additional authors not shown)

Abstract: We present a Bayesian calibration algorithm for CMB observations as implemented within the global end-to-end BeyondPlanck (BP) framework, and apply this to the Planck Low Frequency Instrument (LFI) data. Following the most recent Planck analysis, we decompose the full time-dependent gain into a sum of three orthogonal components: One absolute calibration term, common to all detectors; one time-ind… ▽ More We present a Bayesian calibration algorithm for CMB observations as implemented within the global end-to-end BeyondPlanck (BP) framework, and apply this to the Planck Low Frequency Instrument (LFI) data. Following the most recent Planck analysis, we decompose the full time-dependent gain into a sum of three orthogonal components: One absolute calibration term, common to all detectors; one time-independent term that can vary between detectors; and one time-dependent component that is allowed to vary between one-hour pointing periods. Each term is then sampled conditionally on all other parameters in the global signal model through Gibbs sampling. The absolute calibration is sampled using only the orbital dipole as a reference source, while the two relative gain components are sampled using the full sky signal, including the orbital and Solar CMB dipoles, CMB fluctuations, and foreground contributions. We discuss various aspects of the data that influence gain estimation, including the dipole/polarization quadrupole degeneracy and anomalous jumps in the instrumental gain. Comparing our solution to previous pipelines, we find good agreement in general, with relative deviations of -0.67% (-0.84%) for 30 GHz, 0.12% (-0.04%) for 44 GHz and -0.03% (-0.64%) for 70 GHz, compared to Planck DR4 (Planck 2018). The deviations we find are within expected error bounds, and we attribute them to differences in data usage and general approach between the pipelines. In particular, the BP calibration is performed globally, resulting in better inter-frequency consistency. Additionally, WMAP observations are used actively in the BP analysis, which breaks degeneracies in the Planck data set and results in better agreement with WMAP. Although our presentation and algorithm are currently oriented toward LFI processing, the procedure is fully generalizable to other experiments. △ Less

Submitted 20 May, 2022; v1 submitted 16 November, 2020; originally announced November 2020.

Comments: 19 pages, 18 figures. All BeyondPlanck products and software will be released publicly at http://beyondplanck.science

BeyondPlanck VI. Noise characterization and modelling

Authors: H. T. Ihle, M. Bersanelli, C. Franceschet, E. Gjerløw, K. J. Andersen, R. Aurlien, R. Banerji, S. Bertocco, M. Brilenkov, M. Carbone, L. P. L. Colombo, H. K. Eriksen, J. R. Eskilt, M. K. Foss, U. Fuskeland, S. Galeotta, M. Galloway, S. Gerakakis, B. Hensley, D. Herman, M. Iacobellis, M. Ieronymaki, J. B. Jewell, A. Karakci, E. Keihänen , et al. (16 additional authors not shown)

Abstract: We present a Bayesian method for estimating instrumental noise parameters and propagating noise uncertainties within the global BeyondPlanck Gibbs sampling framework, and apply this to Planck Low Frequency Instrument (LFI) time-ordered data. Following previous literature, we initially adopt a $1/f$ model for the noise power spectral density (PSD), but find the need for an additional lognormal comp… ▽ More We present a Bayesian method for estimating instrumental noise parameters and propagating noise uncertainties within the global BeyondPlanck Gibbs sampling framework, and apply this to Planck Low Frequency Instrument (LFI) time-ordered data. Following previous literature, we initially adopt a $1/f$ model for the noise power spectral density (PSD), but find the need for an additional lognormal component in the noise model for the 30 and 44\,GHz bands. We implement an optimal Wiener-filter (or constrained realization) gap-filling procedure to account for masked data. We then use this procedure to both estimate the gapless correlated noise in the time-domain, $n_\mathrm{corr}$, and to sample the noise PSD parameters, $ξ^n = \{σ_0, f_\mathrm{knee}, α, A_\mathrm{p}\}$. In contrast to previous \textit{Planck} analyses, we assume piecewise stationary noise only within each pointing period (PID), not throughout the full mission, but we adopt the LFI Data Processing Center (DPC) results as priors on $α$ and $f_\mathrm{knee}$. On average, we find best-fit correlated noise parameters that are mostly consistent with previous results, with a few notable exceptions. However, a detailed inspection of the time-dependent results reveals many important findings. First and foremost, we find strong evidence for statistically significant temporal variations in all noise PSD parameters, many of which are directly correlated with satellite housekeeping data. Second, while the simple $1/f$ model appears to be an excellent fit for the LFI 70 GHz channel, there is evidence for additional correlated noise not described by a $1/f$ model in the 30 and 44 GHz channels, including within the primary science frequency range of 0.1--1 Hz. (Abridged) △ Less

Submitted 19 July, 2022; v1 submitted 12 November, 2020; originally announced November 2020.

Comments: 26 pages, 29 figures, as accepted in A&A

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

BeyondPlanck II. CMB map-making through Gibbs sampling

Authors: E. Keihänen, A. -S. Suur-Uski, K. J. Andersen, R. Aurlien, R. Banerji, M. Bersanelli, S. Bertocco, M. Brilenkov, M. Carbone, L. P. L. Colombo, H. K. Eriksen, M. K. Foss, C. Franceschet, U. Fuskeland, S. Galeotta, M. Galloway, S. Gerakakis, E. Gjerløw, B. Hensley, D. Herman, M. Iacobellis, M. Ieronymaki, H. T. Ihle, J. B. Jewell, A. Karakci , et al. (15 additional authors not shown)

Abstract: We present a Gibbs sampling solution to the map-making problem for CMB measurements, building on existing destriping methodology. Gibbs sampling breaks the computationally heavy destriping problem into two separate steps; noise filtering and map binning. Considered as two separate steps, both are computationally much cheaper than solving the combined problem. This provides a huge performance benef… ▽ More We present a Gibbs sampling solution to the map-making problem for CMB measurements, building on existing destriping methodology. Gibbs sampling breaks the computationally heavy destriping problem into two separate steps; noise filtering and map binning. Considered as two separate steps, both are computationally much cheaper than solving the combined problem. This provides a huge performance benefit as compared to traditional methods, and allows us for the first time to bring the destriping baseline length to a single sample. We apply the Gibbs procedure to simulated Planck 30 GHz data. We find that gaps in the time-ordered data are handled efficiently by filling them with simulated noise as part of the Gibbs process. The Gibbs procedure yields a chain of map samples, from which we may compute the posterior mean as a best-estimate map. The variation in the chain provides information on the correlated residual noise, without need to construct a full noise covariance matrix. However, if only a single maximum-likelihood frequency map estimate is required, we find that traditional conjugate gradient solvers converge much faster than a Gibbs sampler in terms of total number of iterations. The conceptual advantages of the Gibbs sampling approach lies in statistically well-defined error propagation and systematic error correction, and this methodology forms the conceptual basis for the map-making algorithm employed in the BeyondPlanck framework, which implements the first end-to-end Bayesian analysis pipeline for CMB observations. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Comments: 11 pages, 10 figures. All BeyondPlanck products and software will be released publicly at http://beyondplanck.science during the online release conference (November 18-20, 2020). Connection details will be made available at the same website. Registration is mandatory for the online tutorial, but optional for the conference

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

BeyondPlanck I. Global Bayesian analysis of the Planck Low Frequency Instrument data

Authors: BeyondPlanck Collaboration, K. J. Andersen, R. Aurlien, R. Banerji, A. Basyrov, M. Bersanelli, S. Bertocco, M. Brilenkov, M. Carbone, L. P. L. Colombo, H. K. Eriksen, J. R. Eskilt, M. K. Foss, C. Franceschet, U. Fuskeland, S. Galeotta, M. Galloway, S. Gerakakis, E. Gjerløw, B. Hensley, D. Herman, M. Iacobellis, M. Ieronymaki, H. T. Ihle, J. B. Jewell , et al. (20 additional authors not shown)

Abstract: We describe the BeyondPlanck project in terms of motivation, methodology and main products, and provide a guide to a set of companion papers that describe each result in fuller detail. We implement a complete end-to-end Bayesian analysis framework for the Planck LFI observations. The primary product is a full joint posterior distribution $P(ω|d)$, where $ω$ represents the set of all free instrumen… ▽ More We describe the BeyondPlanck project in terms of motivation, methodology and main products, and provide a guide to a set of companion papers that describe each result in fuller detail. We implement a complete end-to-end Bayesian analysis framework for the Planck LFI observations. The primary product is a full joint posterior distribution $P(ω|d)$, where $ω$ represents the set of all free instrumental, astrophysical, and cosmological parameters. Notable advantages of this approach are seamless end-to-end propagation of uncertainties; accurate modeling of both astrophysical and instrumental effects in the most natural basis for each uncertain quantity; optimized computational costs with little or no need for intermediate human interaction between various analysis steps; and a complete overview of the entire analysis process within one single framework. We focus in particular on low-$\ell$ CMB polarization reconstruction with Planck LFI. We identify several important new effects that have not been accounted for in previous pipelines, including gain over-smoothing and time-variable and non-$1/f$ correlated noise in the 30 and 44 GHz channels. We find that all results are consistent with the $Λ$CDM model, and we constrain the reionization optical depth to $τ=0.066\pm0.013$, with a low-resolution $χ^2$ probability-to-exceed of 32%. This uncertainty is about 30% larger than the official pipelines, arising from taking into account a more complete instrumental model. The marginal CMB Solar dipole amplitude is $3362.7\pm1.4μ\mathrm{K}$, where the error bar is derived directly from the posterior distribution without the need of any ad-hoc instrumental corrections. We are currently not aware of any significant unmodelled systematic effects remaining in the Planck LFI data, and, for the first time, the 44 GHz channel is fully exploited. (Abridged.) △ Less

Submitted 12 September, 2022; v1 submitted 11 November, 2020; originally announced November 2020.

Comments: 61 pages, 31 figures, submitted to A&A. Major revision: Previously reported unresolved effects are now understood and mitigated. This paper version concludes the BeyondPlanck project. All BeyondPlanck products and software are available at https://cosmoglobe.uio.no

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

Wavefront sensing and control in space-based coronagraph instruments using Zernike's phase-contrast method

Authors: Garreth Ruane, J. Kent Wallace, John Steeves, Camilo Mejia Prada, Byoung-Joon Seo, Eduardo Bendek, Carl Coker, Pin Chen, Brendan Crill, Jeff Jewell, Brian Kern, David Marx, Phillip K. Poon, David Redding, A J Eldorado Riggs, Nicholas Siegler, Robert Zimmer

Abstract: Future space telescopes with coronagraph instruments will use a wavefront sensor (WFS) to measure and correct for phase errors and stabilize the stellar intensity in high-contrast images. The HabEx and LUVOIR mission concepts baseline a Zernike wavefront sensor (ZWFS), which uses Zernike's phase contrast method to convert phase in the pupil into intensity at the WFS detector. In preparation for th… ▽ More Future space telescopes with coronagraph instruments will use a wavefront sensor (WFS) to measure and correct for phase errors and stabilize the stellar intensity in high-contrast images. The HabEx and LUVOIR mission concepts baseline a Zernike wavefront sensor (ZWFS), which uses Zernike's phase contrast method to convert phase in the pupil into intensity at the WFS detector. In preparation for these potential future missions, we experimentally demonstrate a ZWFS in a coronagraph instrument on the Decadal Survey Testbed in the High Contrast Imaging Testbed facility at NASA's Jet Propulsion Laboratory. We validate that the ZWFS can measure low- and mid-spatial frequency aberrations up to the control limit of the deformable mirror, with surface height sensitivity as small as 1 pm, using a configuration similar to the HabEx and LUVOIR concepts. Furthermore, we demonstrate closed-loop control, resolving an individual DM actuator, with residuals consistent with theoretical models. In addition, we predict the expected performance of a ZWFS on future space telescopes using natural starlight from a variety of spectral types. The most challenging scenarios require ~1 hr of integration time to achieve picometer sensitivity. This timescale may be drastically reduced by using internal or external laser sources for sensing purposes. The experimental results and theoretical predictions presented here advance the WFS technology in the context of the next generation of space telescopes with coronagraph instruments. △ Less

Submitted 3 November, 2020; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: Accepted for publication in JATIS

Journal ref: J. Astron. Telesc. Instrum. Syst. 6(4), 045005 (2020)

Astro2020 APC White Paper: The need for better tools to design future CMB experiments

Authors: G. Rocha, A. J. Banday, R. Belen Barreiro, A. Challinor, K. M. Górski, B. Hensley, T. Jaffe, J. Jewell, B. Keating, A. Kogut, C. Lawrence, G. Panopoulou, B. Partridge, T. Pearson, J. Silk, P. Steinhardt, I. Whehus, J. Bock, B. Crill, J. Delabrouille, O. Doré, R. Fernandez-Cobos, A. Ijjas, R. Keskitalo, A. Kritsuk , et al. (5 additional authors not shown)

Abstract: This white paper addresses key challenges for the design of next-decade Cosmic Microwave Background (CMB) experiments, and for assessing their capability to extract cosmological information from CMB polarization. We focus here on the challenges posed by foreground emission, CMB lensing, and instrumental systematics to detect the signal that arises from gravitational waves sourced by inflation and… ▽ More This white paper addresses key challenges for the design of next-decade Cosmic Microwave Background (CMB) experiments, and for assessing their capability to extract cosmological information from CMB polarization. We focus here on the challenges posed by foreground emission, CMB lensing, and instrumental systematics to detect the signal that arises from gravitational waves sourced by inflation and parameterized by $r$, at the level of $r \sim 10^{-3}$ or lower, as proposed for future observational efforts. We argue that more accurate and robust analysis and simulation tools are required for these experiments to realize their promise. We are optimistic that the capability to simulate the joint impact of foregrounds, CMB lensing, and systematics can be developed to the level necessary to support the design of a space mission at $r \sim 10^{-4}$ in a few years. We make the case here for supporting such work. Although ground-based efforts present additional challenges (e.g., atmosphere, ground pickup), which are not addressed here, they would also benefit from these improved simulation capabilities. △ Less

Submitted 5 August, 2019; originally announced August 2019.

Comments: 10 pages + references, 4 figures; Submitted to the Astro2020 call for APC white papers

Enabling the next generation of scientific discoveries by embracing photonic technologies

Authors: Nemanja Jovanovic, Charles Beichman, Cullen Blake, Michael Bottom, Jeffrey Chilcote, Carl Coker, Jonathan Crass, Justin R. Crepp, Nick Cvetojevic, Miguel Daal, Mario Dagenais, Kristina Davis, Richard Dekany, Don Figer, Michael P. Fitzgerald, Pradip Gatkine, Olivier Guyon, Sam Halverson, Robert J. Harris, Philip M. Hinz, David Hover, Andrew W. Howard, Rebecca Jensen-Clem, Jeffrey Jewell, Colby Jurgenson , et al. (24 additional authors not shown)

Abstract: The fields of Astronomy and Astrophysics are technology limited, where the advent and application of new technologies to astronomy usher in a flood of discoveries altering our understanding of the Universe (e.g., recent cases include LIGO and the GRAVITY instrument at the VLTI). Currently, the field of astronomical spectroscopy is rapidly approaching an impasse: the size and cost of instruments, e… ▽ More The fields of Astronomy and Astrophysics are technology limited, where the advent and application of new technologies to astronomy usher in a flood of discoveries altering our understanding of the Universe (e.g., recent cases include LIGO and the GRAVITY instrument at the VLTI). Currently, the field of astronomical spectroscopy is rapidly approaching an impasse: the size and cost of instruments, especially multi-object and integral field spectrographs for extremely large telescopes (ELTs), are pushing the limits of what is feasible, requiring optical components at the very edge of achievable size and performance. For these reasons, astronomers are increasingly looking for innovative solutions like photonic technologies that promote instrument miniaturization and simplification, while providing superior performance. Astronomers have long been aware of the potential of photonic technologies. The goal of this white paper is to draw attention to key photonic technologies and developments over the past two decades and demonstrate there is new momentum in this arena. We outline where the most critical efforts should be focused over the coming decade in order to move towards realizing a fully photonic instrument. A relatively small investment in this technology will advance astronomical photonics to a level where it can reliably be used to solve challenging instrument design limitations. For the benefit of both ground and space borne instruments alike, an endorsement from the National Academy of Sciences decadal survey will ensure that such solutions are set on a path to their full scientific exploitation, which may one day address a broad range of science cases outlined in the KSPs. △ Less

Submitted 13 August, 2019; v1 submitted 17 July, 2019; originally announced July 2019.

Comments: 11 pages, 3 figures, NAS astro2020 whitepaper

On the Use of Planetary Science Data for Studying Extrasolar Planets

Authors: Daniel J. Crichton, J. Steve Hughes, Gael Roudier, Robert A. West, Jeffrey Jewell, Geoffrey Bryden, Mark Swain, T. Joseph W. Lazio

Abstract: There is an opportunity to advance both solar system and extrasolar planetary studies that does not require the construction of new telescopes or new missions but better use and access to inter-disciplinary data sets. This approach leverages significant investment from NASA and international space agencies in exploring this solar system and using those discoveries as "ground truth" for the study o… ▽ More There is an opportunity to advance both solar system and extrasolar planetary studies that does not require the construction of new telescopes or new missions but better use and access to inter-disciplinary data sets. This approach leverages significant investment from NASA and international space agencies in exploring this solar system and using those discoveries as "ground truth" for the study of extrasolar planets. This white paper illustrates the potential, using phase curves and atmospheric modeling as specific examples. A key advance required to realize this potential is to enable seamless discovery and access within and between planetary science and astronomical data sets. Further, seamless data discovery and access also expands the availability of science, allowing researchers and students at a variety of institutions, equipped only with Internet access and a decent computer to conduct cutting-edge research. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

Characterization of the Hamamatsu VUV4 MPPCs for nEXO

Authors: G. Gallina, P. Giampa, F. Retiere, J. Kroeger, G. Zhang, M. Ward, P. Margetak, G. Lic, T. Tsang, L. Doria, S. Al Kharusi, M. Alfaris, G. Anton, I. J. Arnquist, I. Badhrees, P. S. Barbeau, D. Beck, V. Belov, T. Bhatta, J. Blatchford, J. P. Brodsky, E. Brown, T. Brunner, G. F. Cao, L. Cao , et al. (126 additional authors not shown)

Abstract: In this paper we report on the characterization of the Hamamatsu VUV4 (S/N: S13370-6152) Vacuum Ultra-Violet (VUV) sensitive Silicon Photo-Multipliers (SiPMs) as part of the development of a solution for the detection of liquid xenon scintillation light for the nEXO experiment. Various SiPM features, such as: dark noise, gain, correlated avalanches, direct crosstalk and Photon Detection Efficiency… ▽ More In this paper we report on the characterization of the Hamamatsu VUV4 (S/N: S13370-6152) Vacuum Ultra-Violet (VUV) sensitive Silicon Photo-Multipliers (SiPMs) as part of the development of a solution for the detection of liquid xenon scintillation light for the nEXO experiment. Various SiPM features, such as: dark noise, gain, correlated avalanches, direct crosstalk and Photon Detection Efficiency (PDE) were measured in a dedicated setup at TRIUMF. SiPMs were characterized in the range $163 \text{ } \text{K} \leq \text{T}\leq 233 \text{ } \text{K}$. At an over voltage of $3.1\pm0.2$ V and at $\text{T}=163 \text{ }\text{K}$ we report a number of Correlated Avalanches (CAs) per pulse in the $1 \upmu\text{s}$ interval following the trigger pulse of $0.161\pm0.005$. At the same settings the Dark-Noise (DN) rate is $0.137\pm0.002 \text{ Hz/mm}^{2}$. Both the number of CAs and the DN rate are within nEXO specifications. The PDE of the Hamamatsu VUV4 was measured for two different devices at $\text{T}=233 \text{ }\text{K}$ for a mean wavelength of $189\pm7\text{ nm}$. At $3.6\pm0.2$ V and $3.5\pm0.2$ V of over voltage we report a PDE of $13.4\pm2.6\text{ }\%$ and $11\pm2\%$, corresponding to a saturation PDE of $14.8\pm2.8\text{ }\%$ and $12.2\pm2.3\%$, respectively. Both values are well below the $24\text{ }\%$ saturation PDE advertised by Hamamatsu. More generally, the second device tested at $3.5\pm0.2$ V of over voltage is below the nEXO PDE requirement. The first one instead yields a PDE that is marginally close to meeting the nEXO specifications. This suggests that with modest improvements the Hamamatsu VUV4 MPPCs could be considered as an alternative to the FBK-LF SiPMs for the final design of the nEXO detector. △ Less

Submitted 7 June, 2019; v1 submitted 8 March, 2019; originally announced March 2019.

Review of high-contrast imaging systems for current and future ground-based and space-based telescopes III. Technology opportunities and pathways

Authors: Frans Snik, Olivier Absil, Pierre Baudoz, Mathilde Beaulieu, Eduardo Bendek, Eric Cady, Brunella Carlomagno, Alexis Carlotti, Nick Cvetojevic, David Doelman, Kevin Fogarty, Raphaël Galicher, Olivier Guyon, Sebastiaan Haffert, Elsa Huby, Jeffrey Jewell, Nemanja Jovanovic, Christoph Keller, Matthew Kenworthy, Justin Knight, Jonas Kühnn, Johan Mazoyer, Kelsey Miller, Mamadou N'Diaye, Barnaby Norris , et al. (8 additional authors not shown)

Abstract: The Optimal Optical CoronagraphWorkshop at the Lorentz Center in September 2017 in Leiden, the Netherlands gathered a diverse group of 25 researchers working on exoplanet instrumentation to stimulate the emergence and sharing of new ideas. This contribution is the final part of a series of three papers summarizing the outcomes of the workshop, and presents an overview of novel optical technologies… ▽ More The Optimal Optical CoronagraphWorkshop at the Lorentz Center in September 2017 in Leiden, the Netherlands gathered a diverse group of 25 researchers working on exoplanet instrumentation to stimulate the emergence and sharing of new ideas. This contribution is the final part of a series of three papers summarizing the outcomes of the workshop, and presents an overview of novel optical technologies and systems that are implemented or considered for high-contrast imaging instruments on both ground-based and space telescopes. The overall objective of high contrast instruments is to provide direct observations and characterizations of exoplanets at contrast levels as extreme as 10^-10. We list shortcomings of current technologies, and identify opportunities and development paths for new technologies that enable quantum leaps in performance. Specifically, we discuss the design and manufacturing of key components like advanced deformable mirrors and coronagraphic optics, and their amalgamation in "adaptive coronagraph" systems. Moreover, we discuss highly integrated system designs that combine contrast-enhancing techniques and characterization techniques (like high-resolution spectroscopy) while minimizing the overall complexity. Finally, we explore extreme implementations using all-photonics solutions for ground-based telescopes and dedicated huge apertures for space telescopes. △ Less

Submitted 18 July, 2018; originally announced July 2018.

Journal ref: Proc. SPIE 10706-91 (2018)

Review of high-contrast imaging systems for current and future ground-based and space-based telescopes II. Common path wavefront sensing/control and Coherent Differential Imaging

Authors: Nemanja Jovanovic, Olivier Absil, Pierre Baudoz, Mathilde Beaulieu, Michael Bottom, Eric Cady, Brunella Carlomagno, Alexis Carlotti, David Doelman, Kevin Fogarty, Raphael Galicher, Olivier Guyon, Sebastiaan Haffert, Elsa Huby, Jeffrey Jewell, Christoph Keller, Matthew A. Kenworthy, Justin Knight, Jonas Kuhn, Kelsey Miller, Johan Mazoyer, Mamadou N'Diaye, Emiel Por, Laurent Pueyo, A J Eldorado Riggs , et al. (6 additional authors not shown)

Abstract: The Optimal Optical Coronagraph (OOC) Workshop held at the Lorentz Center in September 2017 in Leiden, the Netherlands, gathered a diverse group of 25 researchers working on exoplanet instrumentation to stimulate the emergence and sharing of new ideas. In this second installment of a series of three papers summarizing the outcomes of the OOC workshop (see also~\citenum{ruane2018,snik2018}), we pre… ▽ More The Optimal Optical Coronagraph (OOC) Workshop held at the Lorentz Center in September 2017 in Leiden, the Netherlands, gathered a diverse group of 25 researchers working on exoplanet instrumentation to stimulate the emergence and sharing of new ideas. In this second installment of a series of three papers summarizing the outcomes of the OOC workshop (see also~\citenum{ruane2018,snik2018}), we present an overview of common path wavefront sensing/control and Coherent Differential Imaging techniques, highlight the latest results, and expose their relative strengths and weaknesses. We layout critical milestones for the field with the aim of enhancing future ground/space based high contrast imaging platforms. Techniques like these will help to bridge the daunting contrast gap required to image a terrestrial planet in the zone where it can retain liquid water, in reflected light around a G type star from space. △ Less

Submitted 18 July, 2018; originally announced July 2018.

Comments: 18 pages, 9 figures

Report number: 10703-67

Journal ref: SPIE, paper 10703-67, 2018

Review of high-contrast imaging systems for current and future ground- and space-based telescopes I. Coronagraph design methods and optical performance metrics

Authors: G. Ruane, A. Riggs, J. Mazoyer, E. H. Por, M. N'Diaye, E. Huby, P. Baudoz, R. Galicher, E. Douglas, J. Knight, B. Carlomagno, K. Fogarty, L. Pueyo, N. Zimmerman, O. Absil, M. Beaulieu, E. Cady, A. Carlotti, D. Doelman, O. Guyon, S. Haffert, J. Jewell, N. Jovanovic, C. Keller, M. A. Kenworthy , et al. (7 additional authors not shown)

Abstract: The Optimal Optical Coronagraph (OOC) Workshop at the Lorentz Center in September 2017 in Leiden, the Netherlands gathered a diverse group of 25 researchers working on exoplanet instrumentation to stimulate the emergence and sharing of new ideas. In this first installment of a series of three papers summarizing the outcomes of the OOC workshop, we present an overview of design methods and optical… ▽ More The Optimal Optical Coronagraph (OOC) Workshop at the Lorentz Center in September 2017 in Leiden, the Netherlands gathered a diverse group of 25 researchers working on exoplanet instrumentation to stimulate the emergence and sharing of new ideas. In this first installment of a series of three papers summarizing the outcomes of the OOC workshop, we present an overview of design methods and optical performance metrics developed for coronagraph instruments. The design and optimization of coronagraphs for future telescopes has progressed rapidly over the past several years in the context of space mission studies for Exo-C, WFIRST, HabEx, and LUVOIR as well as ground-based telescopes. Design tools have been developed at several institutions to optimize a variety of coronagraph mask types. We aim to give a broad overview of the approaches used, examples of their utility, and provide the optimization tools to the community. Though it is clear that the basic function of coronagraphs is to suppress starlight while maintaining light from off-axis sources, our community lacks a general set of standard performance metrics that apply to both detecting and characterizing exoplanets. The attendees of the OOC workshop agreed that it would benefit our community to clearly define quantities for comparing the performance of coronagraph designs and systems. Therefore, we also present a set of metrics that may be applied to theoretical designs, testbeds, and deployed instruments. We show how these quantities may be used to easily relate the basic properties of the optical instrument to the detection significance of the given point source in the presence of realistic noise. △ Less

Submitted 18 July, 2018; originally announced July 2018.

Comments: To appear in Proceedings of the SPIE, vol. 10698

Journal ref: Proc. SPIE 10698, 106982S (2018)

Characterization of microdot apodizers for imaging exoplanets with next-generation space telescopes

Authors: Manxuan Zhang, Garreth Ruane, Jacques-Robert Delorme, Dimitri Mawet, Nemanja Jovanavic, Jeffrey Jewell, Stuart Shaklan, J. Kent Wallace

Abstract: A major science goal of future, large-aperture, optical space telescopes is to directly image and spectroscopically analyze reflected light from potentially habitable exoplanets. To accomplish this, the optical system must suppress diffracted light from the star to reveal point sources approximately ten orders of magnitude fainter than the host star at small angular separation. Coronagraphs with m… ▽ More A major science goal of future, large-aperture, optical space telescopes is to directly image and spectroscopically analyze reflected light from potentially habitable exoplanets. To accomplish this, the optical system must suppress diffracted light from the star to reveal point sources approximately ten orders of magnitude fainter than the host star at small angular separation. Coronagraphs with microdot apodizers achieve the theoretical performance needed to image Earth-like planets with a range of possible telescope designs, including those with obscured and segmented pupils. A test microdot apodizer with various bulk patterns (step functions, gradients, and sinusoids) and 4 different dot sizes (3, 5, 7, and 10 $μ$m) made of small chrome squares on anti-reflective glass was characterized with microscopy, optical laser interferometry, as well as transmission and reflectance measurements at wavelengths of 600 and 800 nm. Microscopy revealed the microdots were fabricated to high precision. Results from laser interferometry showed that the phase shifts observed in reflection vary with the local microdot fill factor. Transmission measurements showed that microdot fill factor and transmission were linearly related for dot sizes >5 $μ$m. However, anomalously high transmittance was measured when the dot size is <5x the wavelength and the fill factor is approximately 50%, where the microdot pattern becomes periodic. The transmission excess is not as prominent in the case of larger dot sizes suggesting that it is likely to be caused by the interaction between the incident field and electronic resonances in the surface of the metallic microdots. We used our empirical models of the microdot apodizers to optimize a second generation of reflective apodizer designs and confirmed that the amplitude and phase of the reflected beam closely matches the ideal wavefront. △ Less

Submitted 18 July, 2018; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: Space Telescopes and Instrumentation 2018: Optical, Infrared, and Millimeter Wave

Journal ref: Proc. SPIE 10698, 106985X (2018)

Fast Linearized Coronagraph Optimizer (FALCO) IV. Coronagraph design survey for obstructed and segmented apertures

Authors: G. Ruane, A. Riggs, C. T. Coker, S. B. Shaklan, E. Sidick, D. Mawet, J. Jewell, K. Balasubramanian, C. C. Stark

Abstract: Coronagraph instruments on future space telescopes will enable the direct detection and characterization of Earth-like exoplanets around Sun-like stars for the first time. The quest for the optimal optical coronagraph designs has made rapid progress in recent years thanks to the Segmented Coronagraph Design and Analysis (SCDA) initiative led by the Exoplanet Exploration Program at NASA's Jet Propu… ▽ More Coronagraph instruments on future space telescopes will enable the direct detection and characterization of Earth-like exoplanets around Sun-like stars for the first time. The quest for the optimal optical coronagraph designs has made rapid progress in recent years thanks to the Segmented Coronagraph Design and Analysis (SCDA) initiative led by the Exoplanet Exploration Program at NASA's Jet Propulsion Laboratory. As a result, several types of high-performance designs have emerged that make use of dual deformable mirrors to (1) correct for optical aberrations and (2) suppress diffracted starlight from obstructions and discontinuities in the telescope pupil. However, the algorithms used to compute the optimal deformable mirror surface tend to be computationally intensive, prohibiting large scale design surveys. Here, we utilize the Fast Linearized Coronagraph Optimizer (FALCO), a tool that allows for rapid optimization of deformable mirror shapes, to explore trade-offs in coronagraph designs for obstructed and segmented space telescopes. We compare designs for representative shaped pupil Lyot and vortex coronagraphs, two of the most promising concepts for the LUVOIR space mission concept. We analyze the optical performance of each design, including their throughput and ability to passively suppress light from partially resolved stars in the presence of low-order aberrations. Our main result is that deformable mirror based apodization can sufficiently suppress diffraction from support struts and inter-segment gaps whose widths are on the order of $\sim$0.1% of the primary mirror diameter to detect Earth-sized planets within a few tens of milliarcseconds from the star. △ Less

Submitted 17 July, 2018; originally announced July 2018.

Comments: To appear in the Proceedings of the SPIE, vol. 10698

Journal ref: Proc. SPIE 10698, 106984U (2018)

Vortex coronagraphs for the Habitable Exoplanet Imaging Mission (HabEx) concept: theoretical performance and telescope requirements

Authors: Garreth Ruane, Dimitri Mawet, Bertrand Mennesson, Jeffrey Jewell, Stuart Shaklan

Abstract: The Habitable Exoplanet Imaging Mission (HabEx) concept requires an optical coronagraph that provides deep starlight suppression over a broad spectral bandwidth, high throughput for point sources at small angular separation, and insensitivity to temporally-varying, low-order aberrations. Vortex coronagraphs are a promising solution that perform optimally on off-axis, monolithic telescopes and may… ▽ More The Habitable Exoplanet Imaging Mission (HabEx) concept requires an optical coronagraph that provides deep starlight suppression over a broad spectral bandwidth, high throughput for point sources at small angular separation, and insensitivity to temporally-varying, low-order aberrations. Vortex coronagraphs are a promising solution that perform optimally on off-axis, monolithic telescopes and may also be designed for segmented telescopes with minor losses in performance. We describe the key advantages of vortex coronagraphs on off-axis telescopes: 1) Unwanted diffraction due to aberrations is passively rejected in several low-order Zernike modes relaxing the wavefront stability requirements for imaging Earth-like planets from <10 to >100 pm rms. 2) Stars with angular diameters >0.1 $λ/D$ may be sufficiently suppressed. 3) The absolute planet throughput is >10%, even for unfavorable telescope architectures. 4) Broadband solutions ($Δλ/λ>0.1$) are readily available for both monolithic and segmented apertures. The latter make use of grayscale apodizers in an upstream pupil plane to provide suppression of diffracted light from amplitude discontinuities in the telescope pupil without inducing additional stroke on the deformable mirrors. We set wavefront stability requirements on the telescope, based on a stellar irradiance threshold set at an angular separation of 3$\pm$0.5 $λ/D$ from the star, and discuss how some requirements may be relaxed by trading robustness to aberrations for planet throughput. △ Less

Submitted 11 March, 2018; originally announced March 2018.

Comments: accepted for publication in the Journal of Astronomical Telescopes, Instruments, and Systems. arXiv admin note: text overlap with arXiv:1708.05787

Journal ref: J. of Astronomical Telescopes, Instruments, and Systems, 4(1), 015004 (2018)

Segmented coronagraph design and analysis (SCDA): an initial design study of apodized vortex coronagraphs

Authors: Garreth Ruane, Jeffrey Jewell, Dimitri Mawet, Stuart Shaklan, Christopher C. Stark

Abstract: The segmented coronagraph design and analysis (SCDA) study is a coordinated effort, led by Stuart Shaklan (JPL) and supported by NASA's Exoplanet Exploration Program (ExEP), to provide efficient coronagraph design concepts for exoplanet imaging with future segmented aperture space telescopes. This document serves as an update on the apodized vortex coronagraph designs devised by the Caltech/JPL SC… ▽ More The segmented coronagraph design and analysis (SCDA) study is a coordinated effort, led by Stuart Shaklan (JPL) and supported by NASA's Exoplanet Exploration Program (ExEP), to provide efficient coronagraph design concepts for exoplanet imaging with future segmented aperture space telescopes. This document serves as an update on the apodized vortex coronagraph designs devised by the Caltech/JPL SCDA team. Apodized vortex coronagraphs come in two flavors, where the apodization is achieved either by use of 1) a gray-scale semi-transparent pupil mask or 2) a pair of deformable mirrors in series. Each approach has attractive benefits. This document presents a comprehensive review of the former type. Future theoretical investigations will further explore the use of deformable mirrors for apodization. △ Less

Submitted 6 December, 2017; originally announced December 2017.

Comments: White Paper (2016-2017)

Performance and sensitivity of vortex coronagraphs on segmented space telescopes

Authors: Garreth Ruane, Dimitri Mawet, Jeffrey Jewell, Stuart Shaklan

Abstract: The detection of molecular species in the atmospheres of earth-like exoplanets orbiting nearby stars requires an optical system that suppresses starlight and maximizes the sensitivity to the weak planet signals at small angular separations. Achieving sufficient contrast performance on a segmented aperture space telescope is particularly challenging due to unwanted diffraction within the telescope… ▽ More The detection of molecular species in the atmospheres of earth-like exoplanets orbiting nearby stars requires an optical system that suppresses starlight and maximizes the sensitivity to the weak planet signals at small angular separations. Achieving sufficient contrast performance on a segmented aperture space telescope is particularly challenging due to unwanted diffraction within the telescope from amplitude and phase discontinuities in the pupil. Apodized vortex coronagraphs are a promising solution that theoretically meet the performance needs for high contrast imaging with future segmented space telescopes. We investigate the sensitivity of apodized vortex coronagraphs to the expected aberrations, including segment co-phasing errors in piston and tip/tilt as well as other low-order and mid-spatial frequency aberrations. Coronagraph designs and their associated telescope requirements are identified for conceptual HabEx and LUVOIR telescope designs. △ Less

Submitted 18 August, 2017; originally announced August 2017.

Comments: To appear in the Proceedings of the SPIE: Techniques and Instrumentation for Detection of Exoplanets VIII

Apodized vortex coronagraph designs for segmented aperture telescopes

Authors: Garreth Ruane, Jeffrey Jewell, Dimitri Mawet, Laurent Pueyo, Stuart Shaklan

Abstract: Current state-of-the-art high contrast imaging instruments take advantage of a number of elegant coronagraph designs to suppress starlight and image nearby faint objects, such as exoplanets and circumstellar disks. The ideal performance and complexity of the optical systems depends strongly on the shape of the telescope aperture. Unfortunately, large primary mirrors tend to be segmented and have v… ▽ More Current state-of-the-art high contrast imaging instruments take advantage of a number of elegant coronagraph designs to suppress starlight and image nearby faint objects, such as exoplanets and circumstellar disks. The ideal performance and complexity of the optical systems depends strongly on the shape of the telescope aperture. Unfortunately, large primary mirrors tend to be segmented and have various obstructions, which limit the performance of most conventional coronagraph designs. We present a new family of vortex coronagraphs with numerically-optimized gray-scale apodizers that provide the sensitivity needed to directly image faint exoplanets with large, segmented aperture telescopes, including the Thirty Meter Telescope (TMT) as well as potential next-generation space telescopes. △ Less

Submitted 21 July, 2016; originally announced July 2016.

Comments: To appear in SPIE proceedings vol. 9912

Journal ref: Proc. SPIE 9912, 99122L (2016)

Cosmological Parameters from CMB Maps without Likelihood Approximation

Authors: Benjamin Racine, Jeffrey B. Jewell, Hans Kristian K. Eriksen, Ingunn K. Wehus

Abstract: We propose an efficient Bayesian MCMC algorithm for estimating cosmological parameters from CMB data without use of likelihood approximations. It builds on a previously developed Gibbs sampling framework that allows for exploration of the joint CMB sky signal and power spectrum posterior, P(s,Cl|d), and addresses a long-standing problem of efficient parameter estimation simultaneously in high and… ▽ More We propose an efficient Bayesian MCMC algorithm for estimating cosmological parameters from CMB data without use of likelihood approximations. It builds on a previously developed Gibbs sampling framework that allows for exploration of the joint CMB sky signal and power spectrum posterior, P(s,Cl|d), and addresses a long-standing problem of efficient parameter estimation simultaneously in high and low signal-to-noise regimes. To achieve this, our new algorithm introduces a joint Markov Chain move in which both the signal map and power spectrum are synchronously modified, by rescaling the map according to the proposed power spectrum before evaluating the Metropolis-Hastings accept probability. Such a move was already introduced by Jewell et al. (2009), who used it to explore low signal-to-noise posteriors. However, they also found that the same algorithm is inefficient in the high signal-to-noise regime, since a brute-force rescaling operation does not account for phase information. This problem is mitigated in the new algorithm by subtracting the Wiener filter mean field from the proposed map prior to rescaling, leaving high signal-to-noise information invariant in the joint step, and effectively only rescaling the low signal-to-noise component. To explore the full posterior, the new joint move is then interleaved with a standard conditional Gibbs sky map move. We apply our new algorithm to simplified simulations for which we can evaluate the exact posterior to study both its accuracy and performance, and find good agreement with the exact posterior; marginal means agree to less than 0.006 sigma, and standard deviations to better than 3%. The Markov Chain correlation length is of the same order of magnitude as those obtained by other standard samplers in the field. △ Less

Submitted 28 March, 2016; v1 submitted 21 December, 2015; originally announced December 2015.

Comments: 9 pages, 3 figures, Published in ApJ

Journal ref: B. Racine et al 2016 ApJ 820 31

Optimized Large-Scale CMB Likelihood And Quadratic Maximum Likelihood Power Spectrum Estimation

Authors: E. Gjerløw, L. P. L. Colombo, H. K. Eriksen, K. M. Górski, A. Gruppuso, J. B. Jewell, S. Plaszczynski, I. K. Wehus

Abstract: We revisit the problem of exact CMB likelihood and power spectrum estimation with the goal of minimizing computational cost through linear compression. This idea was originally proposed for CMB purposes by Tegmark et al.\ (1997), and here we develop it into a fully working computational framework for large-scale polarization analysis, adopting \WMAP\ as a worked example. We compare five different… ▽ More We revisit the problem of exact CMB likelihood and power spectrum estimation with the goal of minimizing computational cost through linear compression. This idea was originally proposed for CMB purposes by Tegmark et al.\ (1997), and here we develop it into a fully working computational framework for large-scale polarization analysis, adopting \WMAP\ as a worked example. We compare five different linear bases (pixel space, harmonic space, noise covariance eigenvectors, signal-to-noise covariance eigenvectors and signal-plus-noise covariance eigenvectors) in terms of compression efficiency, and find that the computationally most efficient basis is the signal-to-noise eigenvector basis, which is closely related to the Karhunen-Loeve and Principal Component transforms, in agreement with previous suggestions. For this basis, the information in 6836 unmasked \WMAP\ sky map pixels can be compressed into a smaller set of 3102 modes, with a maximum error increase of any single multipole of 3.8\% at $\ell\le32$, and a maximum shift in the mean values of a joint distribution of an amplitude--tilt model of 0.006$σ$. This compression reduces the computational cost of a single likelihood evaluation by a factor of 5, from 38 to 7.5 CPU seconds, and it also results in a more robust likelihood by implicitly regularizing nearly degenerate modes. Finally, we use the same compression framework to formulate a numerically stable and computationally efficient variation of the Quadratic Maximum Likelihood implementation that requires less than 3 GB of memory and 2 CPU minutes per iteration for $\ell \le 32$, rendering low-$\ell$ QML CMB power spectrum analysis fully tractable on a standard laptop. △ Less

Submitted 10 September, 2015; v1 submitted 13 June, 2015; originally announced June 2015.

Comments: 13 pages, 13 figures, accepted by ApJS

A Joint Analysis of BICEP2/Keck Array and Planck Data

Authors: BICEP2/Keck, Planck Collaborations, :, P. A. R. Ade, N. Aghanim, Z. Ahmed, R. W. Aikin, K. D. Alexander, M. Arnaud, J. Aumont, C. Baccigalupi, A. J. Banday, D. Barkats, R. B. Barreiro, J. G. Bartlett, N. Bartolo, E. Battaner, K. Benabed, A. Benoît, A. Benoit-Lévy, S. J. Benton, J. -P. Bernard, M. Bersanelli, P. Bielewicz, C. A. Bischoff , et al. (254 additional authors not shown)

Abstract: We report the results of a joint analysis of data from BICEP2/Keck Array and Planck. BICEP2 and Keck Array have observed the same approximately 400 deg$^2$ patch of sky centered on RA 0h, Dec. $-57.5°$. The combined maps reach a depth of 57 nK deg in Stokes $Q$ and $U$ in a band centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 to 353 GHz, but much… ▽ More We report the results of a joint analysis of data from BICEP2/Keck Array and Planck. BICEP2 and Keck Array have observed the same approximately 400 deg$^2$ patch of sky centered on RA 0h, Dec. $-57.5°$. The combined maps reach a depth of 57 nK deg in Stokes $Q$ and $U$ in a band centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 to 353 GHz, but much less deeply in any given region (1.2 $μ$K deg in $Q$ and $U$ at 143 GHz). We detect 150$\times$353 cross-correlation in $B$-modes at high significance. We fit the single- and cross-frequency power spectra at frequencies $\geq 150$ GHz to a lensed-$Λ$CDM model that includes dust and a possible contribution from inflationary gravitational waves (as parameterized by the tensor-to-scalar ratio $r$), using a prior on the frequency spectral behavior of polarized dust emission from previous \planck\ analysis of other regions of the sky. We find strong evidence for dust and no statistically significant evidence for tensor modes. We probe various model variations and extensions, including adding a synchrotron component in combination with lower frequency data, and find that these make little difference to the $r$ constraint. Finally we present an alternative analysis which is similar to a map-based cleaning of the dust contribution, and show that this gives similar constraints. The final result is expressed as a likelihood curve for $r$, and yields an upper limit $r_{0.05}<0.12$ at 95% confidence. Marginalizing over dust and $r$, lensing $B$-modes are detected at $7.0\,σ$ significance. △ Less

Submitted 14 April, 2015; v1 submitted 2 February, 2015; originally announced February 2015.

Comments: Version updated to match published. Data and figures available for download at http://bicepkeck.org and http://pla.esac.esa.int/pla

Journal ref: Phys. Rev. Lett. 114, 101301 (2015)

Planck intermediate results. XXXIV. The magnetic field structure in the Rosette Nebula

Authors: Planck Collaboration, N. Aghanim, M. I. R. Alves, M. Arnaud, D. Arzoumanian, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, N. Bartolo, E. Battaner, K. Benabed, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, A. Bracco, C. Burigana, E. Calabrese , et al. (175 additional authors not shown)

Abstract: Planck has mapped the polarized dust emission over the whole sky, making it possible to trace the Galactic magnetic field structure that pervades the interstellar medium (ISM). We combine polarization data from Planck with rotation measure (RM) observations towards a massive star-forming region, the Rosette Nebula in the Monoceros molecular cloud, to study its magnetic field structure and the impa… ▽ More Planck has mapped the polarized dust emission over the whole sky, making it possible to trace the Galactic magnetic field structure that pervades the interstellar medium (ISM). We combine polarization data from Planck with rotation measure (RM) observations towards a massive star-forming region, the Rosette Nebula in the Monoceros molecular cloud, to study its magnetic field structure and the impact of an expanding HII region on the morphology of the field. We derive an analytical solution for the magnetic field, assumed to evolve from an initially uniform configuration following the expansion of ionized gas and the formation of a shell of swept-up ISM. From the RM data we estimate a mean value of the line-of-sight component of the magnetic field of about 3microG (towards the observer) in the Rosette Nebula, for a uniform electron density of about 12cm-3. The dust shell that surrounds the Rosette HII region is clearly observed in the Planck intensity map at 353 GHz, with a polarization signal significantly different from that of the local background when considered as a whole. The Planck observations constrain the plane-of-the-sky orientation of the magnetic field in the Rosette's parent molecular cloud to be mostly aligned with the large-scale field along the Galactic plane. The Planck data are compared with the analytical model, which predicts the mean polarization properties of a spherical and uniform dust shell for a given orientation of the field. This comparison leads to an upper limit of about 45degr on the angle between the line of sight and the magnetic field in the Rosette complex, for an assumed intrinsic dust polarization fraction of 4%. This field direction can reproduce the RM values detected in the ionized region if the magnetic field strength in the Monoceros molecular cloud is in the range 6.5--9microG. △ Less

Submitted 15 September, 2015; v1 submitted 5 January, 2015; originally announced January 2015.

Comments: 16 pages, 11 figures. A&A accepted. Corresponding author: M.I.R. Alves

Journal ref: A&A 586, A137 (2016)

A CMB Gibbs sampler for localized secondary anisotropies

Authors: Philip Bull, Ingunn K. Wehus, Hans Kristian Eriksen, Pedro G. Ferreira, Unni Fuskeland, Krzysztof M. Gorski, Jeffrey B. Jewell

Abstract: As well as primary fluctuations, CMB temperature maps contain a wealth of additional information in the form of secondary anisotropies. Secondary effects that can be identified with individual objects, such as the thermal and kinetic Sunyaev-Zel'dovich (SZ) effects due to galaxy clusters, are difficult to unambiguously disentangle from foreground contamination and the primary CMB however. We devel… ▽ More As well as primary fluctuations, CMB temperature maps contain a wealth of additional information in the form of secondary anisotropies. Secondary effects that can be identified with individual objects, such as the thermal and kinetic Sunyaev-Zel'dovich (SZ) effects due to galaxy clusters, are difficult to unambiguously disentangle from foreground contamination and the primary CMB however. We develop a Bayesian formalism for rigorously characterising anisotropies that are localised on the sky, taking the TSZ and KSZ effects as an example. Using a Gibbs sampling scheme, we are able to efficiently sample from the joint posterior distribution for a multi-component model of the sky with many thousands of correlated physical parameters. The posterior can then be exactly marginalised to estimate properties of the secondary anisotropies, fully taking into account degeneracies with the other signals in the CMB map. We show that this method is computationally tractable using a simple implementation based on the existing Commander component separation code, and also discuss how other types of secondary anisotropy can be accommodated within our framework. △ Less

Submitted 15 October, 2015; v1 submitted 9 October, 2014; originally announced October 2014.

Comments: Updated to published version. Sections added on instrumental properties and effects of IR sources

Journal ref: ApJS, 219, 10 (2015)

Planck intermediate results. XXX. The angular power spectrum of polarized dust emission at intermediate and high Galactic latitudes

Authors: Planck Collaboration, R. Adam, P. A. R. Ade, N. Aghanim, M. Arnaud, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, N. Bartolo, E. Battaner, K. Benabed, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, A. Bracco, M. Bucher , et al. (208 additional authors not shown)

Abstract: The polarized thermal emission from Galactic dust is the main foreground present in measurements of the polarization of the cosmic microwave background (CMB) at frequencies above 100GHz. We exploit the Planck HFI polarization data from 100 to 353GHz to measure the dust angular power spectra $C_\ell^{EE,BB}$ over the range $40<\ell<600$ well away from the Galactic plane. These will bring new insigh… ▽ More The polarized thermal emission from Galactic dust is the main foreground present in measurements of the polarization of the cosmic microwave background (CMB) at frequencies above 100GHz. We exploit the Planck HFI polarization data from 100 to 353GHz to measure the dust angular power spectra $C_\ell^{EE,BB}$ over the range $40<\ell<600$ well away from the Galactic plane. These will bring new insights into interstellar dust physics and a precise determination of the level of contamination for CMB polarization experiments. We show that statistical properties of the emission can be characterized over large fractions of the sky using $C_\ell$. For the dust, they are well described by power laws in $\ell$ with exponents $α^{EE,BB}=-2.42\pm0.02$. The amplitudes of the polarization $C_\ell$ vary with the average brightness in a way similar to the intensity ones. The dust polarization frequency dependence is consistent with modified blackbody emission with $β_d=1.59$ and $T_d=19.6$K. We find a systematic ratio between the amplitudes of the Galactic $B$- and $E$-modes of 0.5. We show that even in the faintest dust-emitting regions there are no "clean" windows where primordial CMB $B$-mode polarization could be measured without subtraction of dust emission. Finally, we investigate the level of dust polarization in the BICEP2 experiment field. Extrapolation of the Planck 353GHz data to 150GHz gives a dust power $\ell(\ell+1)C_\ell^{BB}/(2π)$ of $1.32\times10^{-2}μ$K$_{CMB}^2$ over the $40<\ell<120$ range; the statistical uncertainty is $\pm0.29$ and there is an additional uncertainty (+0.28,-0.24) from the extrapolation, both in the same units. This is the same magnitude as reported by BICEP2 over this $\ell$ range, which highlights the need for assessment of the polarized dust signal even in the cleanest windows of the sky. △ Less

Submitted 8 December, 2014; v1 submitted 19 September, 2014; originally announced September 2014.

Comments: 25 pages, 25 figures, accepted in A&A

Journal ref: A&A 586, A133 (2016)

Planck 2013 results. XI. All-sky model of thermal dust emission

Authors: Planck Collaboration, A. Abergel, P. A. R. Ade, N. Aghanim, M. I. R. Alves, G. Aniano, C. Armitage-Caplan, M. Arnaud, M. Ashdown, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoît, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. Bobin, J. J. Bock, A. Bonaldi , et al. (222 additional authors not shown)

Abstract: This paper presents an all-sky model of dust emission from the Planck 857, 545 and 353 GHz, and IRAS 100 micron data. Using a modified black-body fit to the data we present all-sky maps of the dust optical depth, temperature, and spectral index over the 353-3000 GHz range. This model is a tight representation of the data at 5 arc min. It shows variations of the order of 30 % compared with the wide… ▽ More This paper presents an all-sky model of dust emission from the Planck 857, 545 and 353 GHz, and IRAS 100 micron data. Using a modified black-body fit to the data we present all-sky maps of the dust optical depth, temperature, and spectral index over the 353-3000 GHz range. This model is a tight representation of the data at 5 arc min. It shows variations of the order of 30 % compared with the widely-used model of Finkbeiner, Davis, and Schlegel. The Planck data allow us to estimate the dust temperature uniformly over the whole sky, providing an improved estimate of the dust optical depth compared to previous all-sky dust model, especially in high-contrast molecular regions. An increase of the dust opacity at 353 GHz, tau_353/N_H, from the diffuse to the denser interstellar medium (ISM) is reported. It is associated with a decrease in the observed dust temperature, T_obs, that could be due at least in part to the increased dust opacity. We also report an excess of dust emission at HI column densities lower than 10^20 cm^-2 that could be the signature of dust in the warm ionized medium. In the diffuse ISM at high Galactic latitude, we report an anti-correlation between tau_353/N_H and T_obs while the dust specific luminosity, i.e., the total dust emission integrated over frequency (the radiance) per hydrogen atom, stays about constant. The implication is that in the diffuse high-latitude ISM tau_353 is not as reliable a tracer of dust column density as we conclude it is in molecular clouds where the correlation of tau_353 with dust extinction estimated using colour excess measurements on stars is strong. To estimate Galactic E(B-V) in extragalactic fields at high latitude we develop a new method based on the thermal dust radiance, instead of the dust optical depth, calibrated to E(B-V) using reddening measurements of quasars deduced from Sloan Digital Sky Survey data. △ Less

Submitted 23 September, 2014; v1 submitted 4 December, 2013; originally announced December 2013.

Comments: 37 pages. accepted version. This paper is one of a set associated with the 2013 data release from the Planck mission

Journal ref: A&A 571, A11 (2014)

A multi-level solver for Gaussian constrained CMB realizations

Authors: D. S. Seljebotn, K. -A. Mardal, J. B. Jewell, H. K. Eriksen, P. Bull

Abstract: We present a multi-level solver for drawing constrained Gaussian realizations or finding the maximum likelihood estimate of the CMB sky, given noisy sky maps with partial sky coverage. The method converges substantially faster than existing Conjugate Gradient (CG) methods for the same problem. For instance, for the 143 GHz Planck frequency channel, only 3 multi-level W-cycles result in an absolute… ▽ More We present a multi-level solver for drawing constrained Gaussian realizations or finding the maximum likelihood estimate of the CMB sky, given noisy sky maps with partial sky coverage. The method converges substantially faster than existing Conjugate Gradient (CG) methods for the same problem. For instance, for the 143 GHz Planck frequency channel, only 3 multi-level W-cycles result in an absolute error smaller than 1 microKelvin in any pixel. Using 16 CPU cores, this translates to a computational expense of 6 minutes wall time per realization, plus 8 minutes wall time for a power spectrum-dependent precomputation. Each additional W-cycle reduces the error by more than an order of magnitude, at an additional computational cost of 2 minutes. For comparison, we have never been able to achieve similar absolute convergence with conventional CG methods for this high signal-to-noise data set, even after thousands of CG iterations and employing expensive preconditioners. The solver is part of the Commander 2 code, which is available with an open source license at http://commander.bitbucket.org/. △ Less

Submitted 24 February, 2014; v1 submitted 24 August, 2013; originally announced August 2013.

Comments: Updated to match published version (no major changes)

Journal ref: ApJS 210, 24 (2014)