Showing 1–50 of 98 results for author: Hensley, B

The TEMPO Survey II: Science Cases Leveraged from a Proposed 30-Day Time Domain Survey of the Orion Nebula with the Nancy Grace Roman Space Telescope

Authors: Melinda Soares-Furtado, Mary Anne Limbach, Andrew Vanderburg, John Bally, Juliette Becker, Anna L. Rosen, Luke G. Bouma, Johanna M. Vos, Steve B. Howell, Thomas G. Beatty, William M. J. Best, Anne Marie Cody, Adam Distler, Elena D'Onghia, René Heller, Brandon S. Hensley, Natalie R. Hinkel, Brian Jackson, Marina Kounkel, Adam Kraus, Andrew W. Mann, Nicholas T. Marston, Massimo Robberto, Joseph E. Rodriguez, Jason H. Steffen , et al. (4 additional authors not shown)

Abstract: The TEMPO (Transiting Exosatellites, Moons, and Planets in Orion) Survey is a proposed 30-day observational campaign using the Nancy Grace Roman Space Telescope. By providing deep, high-resolution, short-cadence infrared photometry of a dynamic star-forming region, TEMPO will investigate the demographics of exosatellites orbiting free-floating planets and brown dwarfs -- a largely unexplored disco… ▽ More The TEMPO (Transiting Exosatellites, Moons, and Planets in Orion) Survey is a proposed 30-day observational campaign using the Nancy Grace Roman Space Telescope. By providing deep, high-resolution, short-cadence infrared photometry of a dynamic star-forming region, TEMPO will investigate the demographics of exosatellites orbiting free-floating planets and brown dwarfs -- a largely unexplored discovery space. Here, we present the simulated detection yields of three populations: extrasolar moon analogs orbiting free-floating planets, exosatellites orbiting brown dwarfs, and exoplanets orbiting young stars. Additionally, we outline a comprehensive range of anticipated scientific outcomes accompanying such a survey. These science drivers include: obtaining observational constraints to test prevailing theories of moon, planet, and star formation; directly detecting widely separated exoplanets orbiting young stars; investigating the variability of young stars and brown dwarfs; constraining the low-mass end of the stellar initial mass function; constructing the distribution of dust in the Orion Nebula and mapping evolution in the near-infrared extinction law; mapping emission features that trace the shocked gas in the region; constructing a dynamical map of Orion members using proper motions; and searching for extragalactic sources and transients via deep extragalactic observations reaching a limiting magnitude of $m_{AB}=29.7$\,mag (F146 filter). △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 15 pages, 6 figures, submitted to OJAp

The Primordial Inflation Explorer (PIXIE): Mission Design and Science Goals

Authors: Alan Kogut, Eric Switzer, Dale Fixsen, Nabila Aghanim, Jens Chluba, Dave Chuss, Jacques Delabrouille, Cora Dvorkin, Brandon Hensley, Colin Hill, Bruno Maffei, Anthony Pullen, Aditya Rotti, Alina Sabyr, Leander Thiele, Ed Wollack, Ioana Zelko

Abstract: The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission concept to measure the energy spectrum and linear polarization of the cosmic microwave background (CMB). A single cryogenic Fourier transform spectrometer compares the sky to an external blackbody calibration target, measuring the Stokes I, Q, U parameters to levels ~200 Jy/sr in each 2.65 degree diameter beam over the full sky… ▽ More The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission concept to measure the energy spectrum and linear polarization of the cosmic microwave background (CMB). A single cryogenic Fourier transform spectrometer compares the sky to an external blackbody calibration target, measuring the Stokes I, Q, U parameters to levels ~200 Jy/sr in each 2.65 degree diameter beam over the full sky, in each of 300 frequency channels from 28 GHz to 6 THz. With sensitivity over 1000 times greater than COBE/FIRAS, PIXIE opens a broad discovery space for the origin, contents, and evolution of the universe. Measurements of small distortions from a CMB blackbody spectrum provide a robust determination of the mean electron pressure and temperature in the universe while constraining processes including dissipation of primordial density perturbations, black holes, and the decay or annihilation of dark matter. Full-sky maps of linear polarization measure the optical depth to reionization at nearly the cosmic variance limit and constrain models of primordial inflation. Spectra with sub-percent absolute calibration spanning microwave to far-IR wavelengths provide a legacy data set for analyses including line intensity mapping of extragalactic emission and the cosmic infrared background amplitude and anisotropy. We describe the PIXIE instrument sensitivity, foreground subtraction, and anticipated science return from both the baseline 2-year mission and a potential extended mission. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 44 pages with 18 figures and 4 tables, including full sensitivity curves. Sensitivity curves are available for download at https://lambda.gsfc.nasa.gov/product/pixie/pixie_baseline_noise_get.html

The Metallicity Dependence of PAH Emission in Galaxies I: Insights from Deep Radial Spitzer Spectroscopy

Authors: Cory M. Whitcomb, J. -D. T. Smith, Karin Sandstrom, Carl A. Starkey, Grant P. Donnelly, Bruce T. Draine, Evan D. Skillman, Daniel A. Dale, Lee Armus, Brandon S. Hensley, Thomas S. -Y. Lai, Robert C. Kennicutt

Abstract: We use deep Spitzer mid-infrared spectroscopic maps of radial strips across three nearby galaxies with well-studied metallicity gradients (M101, NGC 628, and NGC 2403) to explore the physical origins of the observed deficit of polycyclic aromatic hydrocarbons (PAHs) at sub-solar metallicity (i.e. the PAH-metallicity relation or PZR). These maps allow us to trace the evolution of all PAH features f… ▽ More We use deep Spitzer mid-infrared spectroscopic maps of radial strips across three nearby galaxies with well-studied metallicity gradients (M101, NGC 628, and NGC 2403) to explore the physical origins of the observed deficit of polycyclic aromatic hydrocarbons (PAHs) at sub-solar metallicity (i.e. the PAH-metallicity relation or PZR). These maps allow us to trace the evolution of all PAH features from 5-18 $μ$m as metallicity decreases continuously from solar ($Z_\odot$) to 0.2 $Z_\odot$. The total PAH to dust luminosity ratio remains relatively constant until reaching a threshold of $\sim$$\frac{2}{3}$$Z_\odot$, below which it declines smoothly but rapidly. The PZR has been attributed to preferential destruction of the smallest grains in the hard radiation environments found at low metallicity. In this scenario, a decrease in emission from the shortest wavelength PAH features is expected. In contrast, we find a steep decline in long wavelength power below $Z_\odot$, especially in the 17 $μ$m feature, with the shorter wavelength PAH bands carrying an increasingly large fraction of power at low metallicity. We use newly developed grain models to reproduce the observed PZR trends, including these variations in fractional PAH feature strengths. The model that best reproduces the data employs an evolving grain size distribution that shifts to smaller sizes as metallicity declines. We interpret this as a result of inhibited grain growth at low metallicity, suggesting continuous replenishment in the interstellar medium is the dominant process shaping the PAH grain population in galaxies. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 22 pages, 13 figures, submitted for publication in ApJ

COMAP Galactic Science I: Observations of Spinning Dust Emission at 30GHz in Dark Clouds Surrounding the λ-Orionis Hii Region

Authors: Stuart E. Harper, Clive Dickinson, Kieran A. Cleary, Brandon S. Hensley, Gabriel A. Hoerning, Roberta Paladini, Thomas J. Rennie, Roke Cepeda-Arroita, Delaney A. Dunne, Hans Kristian Eriksen, Joshua Ott Gundersen, Havard T. Ihle, Jonas G. S. Lunde, Roberto Ricci, Jeroen Stil, Nils-Ole Stutzer, A. R. Taylor, Ingunn Kathrine Wehus

Abstract: Anomalous Microwave Emission (AME) is a major component of Galactic emission in the frequency band 10 to 60 GHz and is commonly modelled as rapidly rotating spinning dust grains. The photodissociation region (PDR) at the boundary of the $λ$-Orionis Hii region has been identified by several recent analyses as one of the brightest spinning dust emitting sources in the sky. We investigate the Barnard… ▽ More Anomalous Microwave Emission (AME) is a major component of Galactic emission in the frequency band 10 to 60 GHz and is commonly modelled as rapidly rotating spinning dust grains. The photodissociation region (PDR) at the boundary of the $λ$-Orionis Hii region has been identified by several recent analyses as one of the brightest spinning dust emitting sources in the sky. We investigate the Barnard 30 dark cloud, a dark cloud embedded within the $λ$-Orionis PDR. We use total-power observations of Barnard 30 from the CO Mapping Array Project (COMAP) pathfinder instrument at 26 to 34GHz with a resolution of 4.5 arcminutes alongside existing data from Planck, WISE, IRAS, ACT, and the 1.447GHz GALFACTS survey. We use aperture photometry and template fitting to measure the spectral energy distribution of Barnard 30. We find that the spinning dust is the dominant emission component in the 26 to 34GHz range at the $7 σ$ level ($S_{30GHz} = 2.85\pm0.43$Jy). We find no evidence that polycyclic aromatic hydrocarbons are the preferred carrier for the spinning dust emission, suggesting that the spinning dust carriers are due to a mixed population of very small grains. Finally, we find evidence for variations in spinning dust emissivity and peak frequency within Barnard 30, and that these variations are possibly driven by changes in dust grain population and the total radiation field. Confirming the origin of the variations in the spinning dust spectrum will require both future COMAP observations at 15GHz combined with spectroscopic mid-infrared data of Barnard 30. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Comments: 21 pages, 13 figures. Submitted to MNRAS

Simulation of the Far-Infrared Polarimetry Approach Envisioned for the PRIMA Mission

Authors: C. Darren Dowell, Brandon S. Hensley, Marc Sauvage

Abstract: Interest in the study of magnetic fields and the properties of interstellar dust, explored through increasingly capable far-IR/submillimeter polarimetry, along with maturing detector technology, have set the stage for a transformative leap in polarization mapping capability using a cryogenic space telescope. We describe the approach pursued by the proposed Probe far-Infrared Mission for Astrophysi… ▽ More Interest in the study of magnetic fields and the properties of interstellar dust, explored through increasingly capable far-IR/submillimeter polarimetry, along with maturing detector technology, have set the stage for a transformative leap in polarization mapping capability using a cryogenic space telescope. We describe the approach pursued by the proposed Probe far-Infrared Mission for Astrophysics (PRIMA) to make ultra-deep maps of intensity and polarization in four bands in the 91-232 micron range. A simple, polarimetry-optimized PRIMA Polarimetric Imager (PPI) is designed for this purpose, consisting of arrays of single-polarization Kinetic Inductance Detectors oriented with three angles which allow measurement of Stokes I, Q, and U in single scans. In this study, we develop an end-to-end observation simulator to perform a realistic test of the approach for the case of mapping a nearby galaxy. The observations take advantage of a beam-steering mirror to perform efficient, two-dimensional, crossing scans. Map making is based on 'destriping' approaches demonstrated for Herschel/SPIRE and Planck. Taking worst-case assumptions for detector sensitivity including 1/f noise, we find excellent recovery of simulated input astrophysical maps, with I, Q, and U detected at near fundamental limits. We describe how PPI performs detector relative calibration and mitigates the key systematic effects to accomplish PRIMA polarization science goals. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: 18 pages, 10 figures, 1 table, to be submitted to PASP

The first degree-scale starlight-polarization-based tomography map of the magnetized interstellar medium

Authors: V. Pelgrims, N. Mandarakas, R. Skalidis, K. Tassis, G. V. Panopoulou, V. Pavlidou, D. Blinov, S. Kiehlmann, S. E. Clark, B. S. Hensley, S. Romanopoulos, A. Basyrov, H. K. Eriksen, M. Falalaki, T. Ghosh, E. Gjerløw, J. A. Kypriotakis, S. Maharana, A. Papadaki, T. J. Pearson, S. B. Potter, A. N. Ramaprakash, A. C. S. Readhead, I. K. Wehus

Abstract: We present the first degree-scale tomography map of the dusty magnetized interstellar medium (ISM) from stellar polarimetry and distance measurements. We used the RoboPol polarimeter at Skinakas Observatory to conduct a survey of starlight polarization in a region of the sky of 4 square degrees. We propose a Bayesian method to decompose the stellar-polarization source field along the distance to i… ▽ More We present the first degree-scale tomography map of the dusty magnetized interstellar medium (ISM) from stellar polarimetry and distance measurements. We used the RoboPol polarimeter at Skinakas Observatory to conduct a survey of starlight polarization in a region of the sky of 4 square degrees. We propose a Bayesian method to decompose the stellar-polarization source field along the distance to invert the 3D volume occupied by the observed stars. We used it to obtain the first 3D map of the dusty magnetized ISM. Specifically, we produced a tomography map of the orientation of the plane-of-sky (POS) component of the magnetic field threading the diffuse, dusty regions responsible for the stellar polarization. For the targeted region centered on Galactic coordinates $(l,b) \approx (103.3^\circ, 22.3^\circ)$, we identified several ISM clouds. Most of the lines of sight intersect more than one cloud. A very nearby component was detected in the foreground of a dominant component from which most of the polarization signal comes. Farther clouds, with a distance of up to 2~kpc, were similarly detected. Some of them likely correspond to intermediate-velocity clouds seen in HI spectra in this region of the sky. We found that the orientation of the POS component of the magnetic field changes along distance for most of the lines of sight. Our study demonstrates that starlight polarization data coupled to distance measures have the power to reveal the great complexity of the dusty magnetized ISM in 3D and, in particular, to provide local measurements of the POS component of the magnetic field. This demonstrates that the inversion of large data volumes, as expected from the PASIPHAE survey, will provide the necessary means to move forward in the modeling of the Galactic magnetic field and of the dusty magnetized ISM as a contaminant in observations of the cosmic microwave background polarization. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: Preprint of the accepted version at Astronomy & Astrophysics. The 3D map obtained in this paper can be visualized online at https://pasiphae.science/visualization and a video featuring it is accessible at https://www.youtube.com/watch?v=dB_6J1zhmPI

Confusion of extragalactic sources in the far infrared: a baseline assessment of the performance of PRIMAger in intensity and polarization

Authors: Matthieu Béthermin, Alberto D. Bolatto, François Boulanger, Charles M. Bradford, Denis Burgarella, Laure Ciesla, James Donnellan, Brandon S. Hensley, Jason Glenn, Guilaine Lagache, Enrique Lopez-Rodriguez, Seb Oliver, Alexandra Pope, Marc Sauvage

Abstract: Because of their limited angular resolution, far-infrared telescopes are usually affected by confusion phenomenon. Since several galaxies can be located in the same instrumental beam, only the brightest objects emerge from the fluctuations caused by fainter sources. The probe far-infrared mission for astrophysics imager (PRIMAger) will observe the mid- and far-infrared (25-235 $μ$m) sky both in in… ▽ More Because of their limited angular resolution, far-infrared telescopes are usually affected by confusion phenomenon. Since several galaxies can be located in the same instrumental beam, only the brightest objects emerge from the fluctuations caused by fainter sources. The probe far-infrared mission for astrophysics imager (PRIMAger) will observe the mid- and far-infrared (25-235 $μ$m) sky both in intensity and polarization. We aim to provide predictions of the confusion level and its consequences for future surveys. We produced simulated PRIMAger maps affected only by the confusion noise using the simulated infrared extragalactic sky (SIDES) semi-empirical simulation. We then estimated the confusion limit in these maps and extracted the sources using a basic blind extractor. By comparing the input galaxy catalog and the extracted source catalog, we derived various performance metrics as completeness, purity, and the accuracy of various measurements. In intensity, we predict that the confusion limit increases rapidly with increasing wavelength. The confusion limit in polarization is more than 100x lower. The measured flux density is dominated by the brightest galaxy in the beam, but other objects also contribute at longer wavelength (~30% at 235 $μ$m). We also show that galaxy clustering has a mild impact on confusion in intensity (up to 25%), while it is negligible in polarization. In intensity, a basic blind extraction will be sufficient to detect galaxies at the knee of the luminosity function up to z~3 and 10$^{11}$ M$_\odot$ main-sequence galaxies up to z~5. In polarization for a conservative sensitivity, we expect ~8 000 detections up to z=2.5 opening a totally new window on the high-z dust polarization. Finally, we show that intensity surveys at short wavelength and polarization surveys at long wavelength tend to reach confusion at similar depth. There is thus a strong synergy. △ Less

Submitted 28 April, 2024; v1 submitted 5 April, 2024; originally announced April 2024.

Comments: 21 pages, 15 figures, 4 tables, submitted to A&A, minor update on the names of the sensitivity scenarii (for consistency with other papers)

The BLAST Observatory: A Sensitivity Study for Far-IR Balloon-borne Polarimeters

Authors: The BLAST Observatory Collaboration, Gabriele Coppi, Simon Dicker, James E. Aguirre, Jason E. Austermann, James A. Beall, Susan E. Clark, Erin G. Cox, Mark J. Devlin, Laura M. Fissel, Nicholas Galitzki, Brandon S. Hensley, Johannes Hubmayr, Sergio Molinari, Federico Nati, Giles Novak, Eugenio Schisano, Juan D. Soler, Carole E. Tucker, Joel N. Ullom, Anna Vaskuri, Michael R. Vissers, Jordan D. Wheeler, Mario Zannoni

Abstract: Sensitive wide-field observations of polarized thermal emission from interstellar dust grains will allow astronomers to address key outstanding questions about the life cycle of matter and energy driving the formation of stars and the evolution of galaxies. Stratospheric balloon-borne telescopes can map this polarized emission at far-infrared wavelengths near the peak of the dust thermal spectrum… ▽ More Sensitive wide-field observations of polarized thermal emission from interstellar dust grains will allow astronomers to address key outstanding questions about the life cycle of matter and energy driving the formation of stars and the evolution of galaxies. Stratospheric balloon-borne telescopes can map this polarized emission at far-infrared wavelengths near the peak of the dust thermal spectrum - wavelengths that are inaccessible from the ground. In this paper we address the sensitivity achievable by a Super Pressure Balloon (SPB) polarimetry mission, using as an example the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) Observatory. By launching from Wanaka, New Zealand, BLAST Observatory can obtain a 30-day flight with excellent sky coverage - overcoming limitations of past experiments that suffered from short flight duration and/or launch sites with poor coverage of nearby star-forming regions. This proposed polarimetry mission will map large regions of the sky at sub-arcminute resolution, with simultaneous observations at 175, 250, and 350 $μm$, using a total of 8274 microwave kinetic inductance detectors. Here, we describe the scientific motivation for the BLAST Observatory, the proposed implementation, and the forecasting methods used to predict its sensitivity. We also compare our forecasted experiment sensitivity with other facilities. △ Less

Submitted 23 May, 2024; v1 submitted 25 January, 2024; originally announced January 2024.

Comments: Published in PASP

Journal ref: 2024 PASP 136 035003

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

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

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

Submitted 20 October, 2023; originally announced October 2023.

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

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

GOALS-JWST: Small neutral grains and enhanced 3.3 micron PAH emission in the Seyfert galaxy NGC 7469

Authors: Thomas S. -Y. Lai, Lee Armus, Marina Bianchin, Tanio Diaz-Santos, Sean T. Linden, George C. Privon, Hanae Inami, Vivian U, Thomas Bohn, Aaron S. Evans, Kirsten L. Larson, Brandon S. Hensley, J. -D. T. Smith, Matthew A. Malkan, Yiqing Song, Sabrina Stierwalt, Paul P. van der Werf, Jed McKinney, Susanne Aalto, Victorine A. Buiten, Jeff Rich, Vassilis Charmandaris, Philip Appleton, Loreto Barcos-Munoz, Torsten Boker , et al. (14 additional authors not shown)

Abstract: We present James Webb Space Telescope (JWST) Near Infrared Spectrograph (NIRSpec) integral-field spectroscopy of the nearby luminous infrared galaxy, NGC 7469. We take advantage of the high spatial/spectral resolution and wavelength coverage of JWST /NIRSpec to study the 3.3 um neutral polycyclic aromatic hydrocarbon (PAH) grain emission on ~60 pc scales. We find a clear change in the average grai… ▽ More We present James Webb Space Telescope (JWST) Near Infrared Spectrograph (NIRSpec) integral-field spectroscopy of the nearby luminous infrared galaxy, NGC 7469. We take advantage of the high spatial/spectral resolution and wavelength coverage of JWST /NIRSpec to study the 3.3 um neutral polycyclic aromatic hydrocarbon (PAH) grain emission on ~60 pc scales. We find a clear change in the average grain properties between the star-forming ring and the central AGN. Regions in the vicinity of the AGN, with [NeIII]/[NeII]>0.25, tend to have larger grain sizes and lower aliphatic-to-aromatic (3.4/3.3) ratios indicating that smaller grains are preferentially removed by photo-destruction in the vicinity of the AGN. We find an overall suppression of the total PAH emission relative to the ionized gas in the central 1 kpc region of the AGN in NGC 7469 compared to what has been observed with Spitzer on 3 kpc scales. However, the fractional 3.3 um to total PAH power is enhanced in the starburst ring, possibly due to a variety of physical effects on sub-kpc scales, including recurrent fluorescence of small grains or multiple photon absorption by large grains. Finally, the IFU data show that while the 3.3 um PAH-derived star formation rate (SFR) in the ring is 8% higher than that inferred from the [NeII] and [NeIII] emission lines, the integrated SFR derived from the 3.3 um feature would be underestimated by a factor of two due to the deficit of PAHs around the AGN, as might occur if a composite system like NGC 7469 were to be observed at high-redshift. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 14 pages, 5 figures, 2 tables, Submitted to ApJL

The Atacama Cosmology Telescope: Galactic Dust Structure and the Cosmic PAH Background in Cross-correlation with WISE

Authors: Rodrigo Córdova Rosado, Brandon S. Hensley, Susan E. Clark, Adriaan J. Duivenvoorden, Zachary Atkins, Elia Stefano Battistelli, Steve K. Choi, Jo Dunkley, Carlos Hervías-Caimapo, Zack Li, Thibaut Louis, Sigurd Naess, Lyman A. Page, Bruce Partridge, Cristóbal Sifón, Suzanne T. Staggs, Cristian Vargas, Edward J. Wollack

Abstract: We present a cross-correlation analysis between $1'$ resolution total intensity and polarization observations from the Atacama Cosmology Telescope (ACT) at 150 and 220 GHz and 15$''$ mid-infrared photometry from the Wide-field Infrared Survey Explorer (WISE) over 107 12.5$^\circ\times$12.5$^\circ$ patches of sky. We detect a spatially isotropic signal in the WISE$\times$ACT $TT$ cross power spectr… ▽ More We present a cross-correlation analysis between $1'$ resolution total intensity and polarization observations from the Atacama Cosmology Telescope (ACT) at 150 and 220 GHz and 15$''$ mid-infrared photometry from the Wide-field Infrared Survey Explorer (WISE) over 107 12.5$^\circ\times$12.5$^\circ$ patches of sky. We detect a spatially isotropic signal in the WISE$\times$ACT $TT$ cross power spectrum at 30$σ$ significance that we interpret as the correlation between the cosmic infrared background at ACT frequencies and polycyclic aromatic hydrocarbon (PAH) emission from galaxies in WISE, i.e., the cosmic PAH background. Within the Milky Way, the Galactic dust $TT$ spectra are generally well-described by power laws in $\ell$ over the range 10$^3 < \ell < $10$^4$, but there is evidence both for variability in the power law index and for non-power law behavior in some regions. We measure a positive correlation between WISE total intensity and ACT $E$-mode polarization at 1000$ < \ell \lesssim $6000 at $>$3$σ$ in each of 35 distinct $\sim$100 deg$^2$ regions of the sky, suggesting alignment between Galactic density structures and the local magnetic field persists to sub-parsec physical scales in these regions. The distribution of $TE$ amplitudes in this $\ell$ range across all 107 regions is biased to positive values, while there is no evidence for such a bias in the $TB$ spectra. This work constitutes the highest-$\ell$ measurements of the Galactic dust $TE$ spectrum to date and indicates that cross-correlation with high-resolution mid-infrared measurements of dust emission is a promising tool for constraining the spatial statistics of dust emission at millimeter wavelengths. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: 20 pages, 14 figures, submitted to ApJ

The Atacama Cosmology Telescope: High-resolution component-separated maps across one-third of the sky

Authors: William R. Coulton, Mathew S. Madhavacheril, Adriaan J. Duivenvoorden, J. Colin Hill, Irene Abril-Cabezas, Peter A. R. Ade, Simone Aiola, Tommy Alford, Mandana Amiri, Stefania Amodeo, Rui An, Zachary Atkins, Jason E. Austermann, Nicholas Battaglia, Elia Stefano Battistelli, James A. Beall, Rachel Bean, Benjamin Beringue, Tanay Bhandarkar, Emily Biermann, Boris Bolliet, J Richard Bond, Hongbo Cai, Erminia Calabrese, Victoria Calafut , et al. (129 additional authors not shown)

Abstract: Observations of the millimeter sky contain valuable information on a number of signals, including the blackbody cosmic microwave background (CMB), Galactic emissions, and the Compton-$y$ distortion due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Extracting new insight into cosmological and astrophysical questions often requires combining multi-wavelength observations to spectrally isolate one… ▽ More Observations of the millimeter sky contain valuable information on a number of signals, including the blackbody cosmic microwave background (CMB), Galactic emissions, and the Compton-$y$ distortion due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Extracting new insight into cosmological and astrophysical questions often requires combining multi-wavelength observations to spectrally isolate one component. In this work, we present a new arcminute-resolution Compton-$y$ map, which traces out the line-of-sight-integrated electron pressure, as well as maps of the CMB in intensity and E-mode polarization, across a third of the sky (around 13,000 sq.~deg.). We produce these through a joint analysis of data from the Atacama Cosmology Telescope (ACT) Data Release 4 and 6 at frequencies of roughly 93, 148, and 225 GHz, together with data from the \textit{Planck} satellite at frequencies between 30 GHz and 545 GHz. We present detailed verification of an internal linear combination pipeline implemented in a needlet frame that allows us to efficiently suppress Galactic contamination and account for spatial variations in the ACT instrument noise. These maps provide a significant advance, in noise levels and resolution, over the existing \textit{Planck} component-separated maps and will enable a host of science goals including studies of cluster and galaxy astrophysics, inferences of the cosmic velocity field, primordial non-Gaussianity searches, and gravitational lensing reconstruction of the CMB. △ Less

Submitted 3 July, 2023; originally announced July 2023.

Comments: The Compton-y map and associated products will be made publicly available upon publication of the paper. The CMB T and E mode maps will be made available when the DR6 maps are made public

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

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

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

Submitted 3 May, 2023; originally announced May 2023.

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

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

The Atacama Cosmology Telescope: DR6 Gravitational Lensing Map and Cosmological Parameters

Authors: Mathew S. Madhavacheril, Frank J. Qu, Blake D. Sherwin, Niall MacCrann, Yaqiong Li, Irene Abril-Cabezas, Peter A. R. Ade, Simone Aiola, Tommy Alford, Mandana Amiri, Stefania Amodeo, Rui An, Zachary Atkins, Jason E. Austermann, Nicholas Battaglia, Elia Stefano Battistelli, James A. Beall, Rachel Bean, Benjamin Beringue, Tanay Bhandarkar, Emily Biermann, Boris Bolliet, J Richard Bond, Hongbo Cai, Erminia Calabrese , et al. (134 additional authors not shown)

Abstract: We present cosmological constraints from a gravitational lensing mass map covering 9400 sq. deg. reconstructed from CMB measurements made by the Atacama Cosmology Telescope (ACT) from 2017 to 2021. In combination with BAO measurements (from SDSS and 6dF), we obtain the amplitude of matter fluctuations $σ_8 = 0.819 \pm 0.015$ at 1.8% precision, $S_8\equivσ_8({Ω_{\rm m}}/0.3)^{0.5}=0.840\pm0.028$ an… ▽ More We present cosmological constraints from a gravitational lensing mass map covering 9400 sq. deg. reconstructed from CMB measurements made by the Atacama Cosmology Telescope (ACT) from 2017 to 2021. In combination with BAO measurements (from SDSS and 6dF), we obtain the amplitude of matter fluctuations $σ_8 = 0.819 \pm 0.015$ at 1.8% precision, $S_8\equivσ_8({Ω_{\rm m}}/0.3)^{0.5}=0.840\pm0.028$ and the Hubble constant $H_0= (68.3 \pm 1.1)\, \text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}$ at 1.6% precision. A joint constraint with CMB lensing measured by the Planck satellite yields even more precise values: $σ_8 = 0.812 \pm 0.013$, $S_8\equivσ_8({Ω_{\rm m}}/0.3)^{0.5}=0.831\pm0.023$ and $H_0= (68.1 \pm 1.0)\, \text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}$. These measurements agree well with $Λ$CDM-model extrapolations from the CMB anisotropies measured by Planck. To compare these constraints to those from the KiDS, DES, and HSC galaxy surveys, we revisit those data sets with a uniform set of assumptions, and find $S_8$ from all three surveys are lower than that from ACT+Planck lensing by varying levels ranging from 1.7-2.1$σ$. These results motivate further measurements and comparison, not just between the CMB anisotropies and galaxy lensing, but also between CMB lensing probing $z\sim 0.5-5$ on mostly-linear scales and galaxy lensing at $z\sim 0.5$ on smaller scales. We combine our CMB lensing measurements with CMB anisotropies to constrain extensions of $Λ$CDM, limiting the sum of the neutrino masses to $\sum m_ν < 0.12$ eV (95% c.l.), for example. Our results provide independent confirmation that the universe is spatially flat, conforms with general relativity, and is described remarkably well by the $Λ$CDM model, while paving a promising path for neutrino physics with gravitational lensing from upcoming ground-based CMB surveys. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Comments: 30 pages, 16 figures, prepared for submission to ApJ. Cosmological likelihood data is here: https://lambda.gsfc.nasa.gov/product/act/actadv_prod_table.html ; likelihood software is here: https://github.com/ACTCollaboration/act_dr6_lenslike . Also see companion papers Qu et al and MacCrann et al. Mass maps will be released when papers are published

The Atacama Cosmology Telescope: A Measurement of the DR6 CMB Lensing Power Spectrum and its Implications for Structure Growth

Authors: Frank J. Qu, Blake D. Sherwin, Mathew S. Madhavacheril, Dongwon Han, Kevin T. Crowley, Irene Abril-Cabezas, Peter A. R. Ade, Simone Aiola, Tommy Alford, Mandana Amiri, Stefania Amodeo, Rui An, Zachary Atkins, Jason E. Austermann, Nicholas Battaglia, Elia Stefano Battistelli, James A. Beall, Rachel Bean, Benjamin Beringue, Tanay Bhandarkar, Emily Biermann, Boris Bolliet, J Richard Bond, Hongbo Cai, Erminia Calabrese , et al. (133 additional authors not shown)

Abstract: We present new measurements of cosmic microwave background (CMB) lensing over $9400$ sq. deg. of the sky. These lensing measurements are derived from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) CMB dataset, which consists of five seasons of ACT CMB temperature and polarization observations. We determine the amplitude of the CMB lensing power spectrum at $2.3\%$ precision ($43σ$ sign… ▽ More We present new measurements of cosmic microwave background (CMB) lensing over $9400$ sq. deg. of the sky. These lensing measurements are derived from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) CMB dataset, which consists of five seasons of ACT CMB temperature and polarization observations. We determine the amplitude of the CMB lensing power spectrum at $2.3\%$ precision ($43σ$ significance) using a novel pipeline that minimizes sensitivity to foregrounds and to noise properties. To ensure our results are robust, we analyze an extensive set of null tests, consistency tests, and systematic error estimates and employ a blinded analysis framework. The baseline spectrum is well fit by a lensing amplitude of $A_{\mathrm{lens}}=1.013\pm0.023$ relative to the Planck 2018 CMB power spectra best-fit $Λ$CDM model and $A_{\mathrm{lens}}=1.005\pm0.023$ relative to the $\text{ACT DR4} + \text{WMAP}$ best-fit model. From our lensing power spectrum measurement, we derive constraints on the parameter combination $S^{\mathrm{CMBL}}_8 \equiv σ_8 \left({Ω_m}/{0.3}\right)^{0.25}$ of $S^{\mathrm{CMBL}}_8= 0.818\pm0.022$ from ACT DR6 CMB lensing alone and $S^{\mathrm{CMBL}}_8= 0.813\pm0.018$ when combining ACT DR6 and Planck NPIPE CMB lensing power spectra. These results are in excellent agreement with $Λ$CDM model constraints from Planck or $\text{ACT DR4} + \text{WMAP}$ CMB power spectrum measurements. Our lensing measurements from redshifts $z\sim0.5$--$5$ are thus fully consistent with $Λ$CDM structure growth predictions based on CMB anisotropies probing primarily $z\sim1100$. We find no evidence for a suppression of the amplitude of cosmic structure at low redshifts △ Less

Submitted 28 May, 2024; v1 submitted 11 April, 2023; originally announced April 2023.

Comments: 45+22 pages, 50 figures. v2 matches with published version in ApJ. Cosmological likelihood data and lensing maps are here: https://lambda.gsfc.nasa.gov/product/act/actadv_prod_table.html ; likelihood software is here: https://github.com/ACTCollaboration/act_dr6_lenslike . Also see companion papers Madhavacheril et al and MacCrann et al

Report number: FERMILAB-PUB-23-237-PPD

Systematic error mitigation for the PIXIE Fourier transform spectrometer

Authors: A. Kogut, Dale Fixsen, Nabila Aghanim, Jens Chluba, David T. Chuss, Jacques Delabrouille, Brandon S. Hensley, J. Colin Hill, Bruno Maffei, Anthony R. Pullen, Aditya Rotti, Eric R. Switzer, Edward J. Wollack, Ioana Zelko

Abstract: The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission concept to measure the spectrum and polarization of the cosmic microwave background. Cosmological signals are small compared to the instantaneous instrument noise, requiring strict control of instrumental signals. The instrument design provides multiple levels of null operation, signal modulation, and signal differences, with o… ▽ More The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission concept to measure the spectrum and polarization of the cosmic microwave background. Cosmological signals are small compared to the instantaneous instrument noise, requiring strict control of instrumental signals. The instrument design provides multiple levels of null operation, signal modulation, and signal differences, with only few-percent systematic error suppression required at each level. Jackknife tests based on discrete instrument symmetries provide an independent means to identify, model, and remove remaining instrumental signals. We use detailed time-ordered simulations, including realistic performance and tolerance parameters, to evaluate the instrument response to broad classes of systematic errors for both spectral distortions and polarization. The largest systematic errors contribute additional white noise at the few-percent level compared to the dominant photon noise. Coherent instrumental effects which do not integrate down are smaller still, and remain several orders of magnitude below the targeted cosmological signals. △ Less

Submitted 31 March, 2023; originally announced April 2023.

Comments: 37 pages including 23 figures

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

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

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

Submitted 14 March, 2023; originally announced March 2023.

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

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

A Framework for Modeling Polycyclic Aromatic Hydrocarbon Emission in Galaxy Evolution Simulations

Authors: Desika Narayanan, J. D. Smith, Brandon Hensley, Qi Li, Chia-Yu Hu, Karin Sandstrom, Paul Torrey, Mark Vogelsberger, Federico Marinacci, Laura Sales

Abstract: We present a new methodology for simulating mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs) in galaxy evolution simulations. To do this, we combine theoretical models of PAH emission features as they respond to varying interstellar radiation fields, grain size distributions, and ionization states with a new on-the-fly model for dust evolution in hydrodynamic galaxy simulations.… ▽ More We present a new methodology for simulating mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs) in galaxy evolution simulations. To do this, we combine theoretical models of PAH emission features as they respond to varying interstellar radiation fields, grain size distributions, and ionization states with a new on-the-fly model for dust evolution in hydrodynamic galaxy simulations. We apply these models to 3 idealized arepo galaxy evolution simulations within the smuggle physics framework. We use these simulations to develop numerical experiments investigating the buildup of PAH masses and luminosities in galaxies in idealized analogs of the Milky Way, a dwarf galaxy, and starburst disk. Our main results follow. Galaxies with high specific star formation rates have increased feedback energy per unit mass, and are able to efficiently shatter dust grains, driving up the fraction of ultra small grains. At the same time, in our model large radiation fields per unit gas density convert aliphatic grains into aromatics. The fraction of dust grains in the form of PAHs (q_PAH) can be understood as a consequence of these processes, and in our model PAHs form primarily from interstellar processing (shattering) of larger grains rather than from the growth of smaller grains. We find that the hardness of the radiation field plays a larger role than variations in the grain size distribution in setting the total integrated PAH luminosities, though cosmological simulations are necessary to fully investigate the complex interplay of processes that drive PAH band luminosities in galaxies. Finally, we highlight feature PAH strength variations, cautioning against the usage of emission templates with constant feature strength ratios. △ Less

Submitted 17 January, 2023; originally announced January 2023.

Comments: Submitted to ApJ; comments welcome

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

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

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

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

Comments: 34 pages, 13 figures, published in JCAP

Journal ref: JCAP06(2023)034

The TEMPO Survey I: Predicting Yields of the Transiting Exosatellites, Moons, and Planets from a 30-day Survey of Orion with the Nancy Grace Roman Space Telescope

Authors: Mary Anne Limbach, Melinda Soares-Furtado, Andrew Vanderburg, William M. J. Best, Ann Marie Cody, Elena D'Onghia, René Heller, Brandon S. Hensley, Marina Kounkel, Adam Kraus, Andrew W. Mann, Massimo Robberto, Anna L. Rosen, Richard Townsend, Johanna M. Vos

Abstract: We present design considerations for the Transiting Exosatellites, Moons, and Planets in Orion (TEMPO) Survey with the Nancy Grace Roman Space Telescope. This proposed 30-day survey is designed to detect a population of transiting extrasolar satellites, moons, and planets in the Orion Nebula Cluster (ONC). The young (1-3 Myr), densely-populated ONC harbors about a thousand bright brown dwarfs (BDs… ▽ More We present design considerations for the Transiting Exosatellites, Moons, and Planets in Orion (TEMPO) Survey with the Nancy Grace Roman Space Telescope. This proposed 30-day survey is designed to detect a population of transiting extrasolar satellites, moons, and planets in the Orion Nebula Cluster (ONC). The young (1-3 Myr), densely-populated ONC harbors about a thousand bright brown dwarfs (BDs) and free-floating planetary-mass objects (FFPs). TEMPO offers sufficient photometric precision to monitor FFPs with ${\rm M}\geq1{\rm M}_{\rm J}$ for transiting satellites. The survey is also capable of detecting FFPs down to sub-Saturn masses via direct imaging, although follow-up confirmation will be challenging. TEMPO yield estimates include 14 (3-22) exomoons/satellites transiting FFPs and 54 (8-100) satellites transiting BDs. Of this population, approximately $50\%$ of companions would be "super-Titans" (Titan to Earth mass). Yield estimates also include approximately $150$ exoplanets transiting young Orion stars, of which $>50\%$ will orbit mid-to-late M dwarfs and approximately ten will be proto-habitable zone, terrestrial ($0.1{\rm M}_{\oplus} - 5{\rm M}_{\oplus}$) exoplanets. TEMPO would provide the first census demographics of small exosatellites orbiting FFPs and BDs, while simultaneously offering insights into exoplanet evolution at the earliest stages. This detected exosatellite population is likely to be markedly different from the current census of exoplanets with similar masses (e.g., Earth-mass exosatellites that still possess H/He envelopes). Although our yield estimates are highly uncertain, as there are no known exoplanets or exomoons analogous to these satellites, the TEMPO survey would test the prevailing theories of exosatellite formation and evolution, which limit the certainty surrounding detection yields. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: Submitted to PASP

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)

The Astrodust+PAH Model: A Unified Description of the Extinction, Emission, and Polarization from Dust in the Diffuse Interstellar Medium

Authors: Brandon S. Hensley, B. T. Draine

Abstract: We present a new model of interstellar dust in which large grains are a single composite material, ``astrodust,'' and nanoparticle-sized grains come in distinct varieties including polycyclic aromatic hydrocarbons (PAHs). We argue that a single-composition model for grains larger than $\sim$0.02 $μ$m most naturally explains the lack of frequency dependence in the far-infrared (FIR) polarization fr… ▽ More We present a new model of interstellar dust in which large grains are a single composite material, ``astrodust,'' and nanoparticle-sized grains come in distinct varieties including polycyclic aromatic hydrocarbons (PAHs). We argue that a single-composition model for grains larger than $\sim$0.02 $μ$m most naturally explains the lack of frequency dependence in the far-infrared (FIR) polarization fraction and the characteristic ratio of optical to FIR polarization. We derive a size distribution and alignment function for 1.4:1 oblate astrodust grains that, with PAHs, reproduce the mean wavelength dependence and polarization of Galactic extinction and emission from the diffuse interstellar medium while respecting constraints on solid phase abundances. All model data and Python-based interfaces are made publicly available. △ Less

Submitted 6 June, 2023; v1 submitted 25 August, 2022; originally announced August 2022.

Comments: 34 pages, 19 figures. Data available at https://dataverse.harvard.edu/dataverse/astrodust

Journal ref: ApJ, 948, 55 (2023)

SOFIA/HAWC+ observations of the Crab Nebula: dust properties from polarised emission

Authors: Jérémy Chastenet, Ilse De Looze, Brandon S. Hensley, Bert Vandenbroucke, Mike J. Barlow, Jeonghee Rho, Aravind P. Ravi, Haley L. Gomez, Florian Kirchschlager, Juan Macías-Pérez, Mikako Matsuura, Kate Pattle, Nicolas Ponthieu, Felix D. Priestley, Monica Relaño, Alessia Ritacco, Roger Wesson

Abstract: Supernova remnants (SNRs) are well-recognised dust producers, but their net dust production rate remains elusive due to uncertainties in grain properties that propagate into observed dust mass uncertainties, and determine how efficiently these grains are processed by reverse shocks. In this paper, we present a detection of polarised dust emission in the Crab pulsar wind nebula, the second SNR with… ▽ More Supernova remnants (SNRs) are well-recognised dust producers, but their net dust production rate remains elusive due to uncertainties in grain properties that propagate into observed dust mass uncertainties, and determine how efficiently these grains are processed by reverse shocks. In this paper, we present a detection of polarised dust emission in the Crab pulsar wind nebula, the second SNR with confirmed polarised dust emission after Cassiopeia A. We constrain the bulk composition of the dust with new SOFIA/HAWC+ polarimetric data in band C 89 um and band D 154 um. After correcting for synchrotron polarisation, we report dust polarisation fractions ranging between 3.7-9.6 per cent and 2.7-7.6 per cent in three individual dusty filaments at 89 and 154 um, respectively. The detected polarised signal suggests the presence of large (> 0.05-0.1 um) grains in the Crab Nebula. With the observed polarisation, and polarised and total fluxes, we constrain the temperatures and masses of carbonaceous and silicate grains. We find that the carbon-rich grain mass fraction varies between 12 and 70 per cent, demonstrating that carbonaceous and silicate grains co-exist in this SNR. Temperatures range from 40 K to 70 K and from 30 K to 50 K for carbonaceous and silicate grains, respectively. Dust masses range from 10^{-4} Msol to 10^{-2} Msol for carbonaceous grains and to 10^{-1} Msol for silicate grains, in three individual regions. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: Accepted for publication in MNRAS. 16 pages, 8 figures + 2 Appendices, 2 figures

Tracing PAH Emission in $λ$-Orionis Using COBE/DIRBE Data

Authors: David T. Chuss, Brandon S. Hensley, Alan J. Kogut, Jordan A. Guerra, Hayley C. Nofi, Javad Siah

Abstract: We use archival COBE/DIRBE data to construct a map of polycyclic aromatic hydrocarbon (PAH) emission in the $λ$-Orionis region. The presence of the 3.3 $μ$m PAH feature within the DIRBE 3.5 $μ$m band and the corresponding lack of significant PAH spectral features in the adjacent DIRBE bands (1.25, 2.2, and 4.9 $μ$m) enable estimation of the PAH contribution to the 3.5 $μ$m data. Having the shortes… ▽ More We use archival COBE/DIRBE data to construct a map of polycyclic aromatic hydrocarbon (PAH) emission in the $λ$-Orionis region. The presence of the 3.3 $μ$m PAH feature within the DIRBE 3.5 $μ$m band and the corresponding lack of significant PAH spectral features in the adjacent DIRBE bands (1.25, 2.2, and 4.9 $μ$m) enable estimation of the PAH contribution to the 3.5 $μ$m data. Having the shortest wavelength of known PAH features, the 3.3 $μ$m feature probes the smallest PAHs, which are also the leading candidates for carriers of anomalous microwave emission (AME). We use this map to investigate the association between the AME and the emission from PAH molecules. We find that the spatial correlation in $λ$-Orionis is higher between AME and far-infrared dust emission (as represented by the DIRBE 240 $μ$m map) than it is between our PAH map and AME. This finding, in agreement with previous studies using PAH features at longer wavelengths, is in tension with the hypothesis that AME is due to spinning PAHs. However, the expected correlation between mid-infrared and microwave emission could potentially be degraded by different sensitivities of each emission mechanism to local environmental conditions even if PAHs are the carriers of both. △ Less

Submitted 7 November, 2022; v1 submitted 18 August, 2022; originally announced August 2022.

Comments: 9 pages, 6 figures, accepted, ApJ

Characterizing line-of-sight variability of polarized dust emission with future CMB experiments

Authors: Lisa McBride, Philip Bull, Brandon S. Hensley

Abstract: Galactic dust emission is often accounted for in cosmic microwave background (CMB) analyses by fitting a simple two-parameter modified blackbody (MBB) model in each pixel, which nominally accounts for the temperature and opacity of the dust. While this may be a good approximation for individual dust clouds, typically a number of such clouds are found along each line of sight and within each angula… ▽ More Galactic dust emission is often accounted for in cosmic microwave background (CMB) analyses by fitting a simple two-parameter modified blackbody (MBB) model in each pixel, which nominally accounts for the temperature and opacity of the dust. While this may be a good approximation for individual dust clouds, typically a number of such clouds are found along each line of sight and within each angular pixel, resulting in a superposition of their spectra. In this paper, we study the effects of this superposition on pixel-based foreground fitting strategies by modelling the spectral energy distribution (SED) in each pixel as the integral of individual MBB spectra over various physically-motivated statistical distributions of dust cloud properties. We show that fitting these SEDs with the simple two-parameter MBB model generally results in unbiased estimates of the CMB Stokes Q and U amplitudes in each pixel, unless there are significant changes in both the dust SED and polarization angle along the line of sight, in which case significant ($ > 10σ$) biases are observed in an illustrative model. We also find that the best-fit values of the dust temperature, $T_d$, and spectral index, $β_d$, are significantly biased away from the mean/median of the corresponding statistical distributions when the distributions are broad, suggesting that MBB model fits can give an unrepresentative picture of the physical properties of the dust at microwave wavelengths if not interpreted carefully. Using a Fisher matrix analysis, we also determine the experimental sensitivity required to recover the parameters of the $T_d$ and $β_d$ distributions themselves by fitting a probabilistic MBB model, finding that only the parameters of broad distributions can be measured by SED fitting on a single line of sight. △ Less

Submitted 28 July, 2022; originally announced July 2022.

Comments: 15 pages, 10 figures; see Table 3 for main results

Constraints on the Optical Depth to Reionization from Balloon-Borne CMB Measurements

Authors: Josquin Errard, Mathieu Remazeilles, Jonathan Aumont, Jacques Delabrouille, Daniel Green, Shaul Hanany, Brandon S. Hensley, Alan Kogut

Abstract: We assess the uncertainty with which a balloon-borne experiment, nominally called Tau Surveyor ($τS$), can measure the optical depth to reionization $σ(τ)$ with given realistic constraints of instrument noise and foreground emissions. Using a $τS$ fiducial design with six frequency bands between 150 and 380 GHz with white and uniform map noise of 7 $μ$K arcmin, achievable with a single mid-latitud… ▽ More We assess the uncertainty with which a balloon-borne experiment, nominally called Tau Surveyor ($τS$), can measure the optical depth to reionization $σ(τ)$ with given realistic constraints of instrument noise and foreground emissions. Using a $τS$ fiducial design with six frequency bands between 150 and 380 GHz with white and uniform map noise of 7 $μ$K arcmin, achievable with a single mid-latitude flight, and including Planck's 30 and 44 GHz data we assess the error $σ(τ)$ obtained with three foreground models and as a function of sky fraction $f_{\rm sky}$ between 40% and 54%. We carry out the analysis using both parametric and blind foreground separation techniques. We compare $σ(τ)$ values to those obtained with low frequency and high frequency versions of the experiment called $τS$-lf and $τS$-hf that have only four and up to eight frequency bands with narrower and wider frequency coverage, respectively. We find that with $τS$ the lowest constraint is $σ(τ)=0.0034$, obtained for one of the foreground models with $f_{\rm sky}$=54%. $σ(τ)$ is larger, in some cases by more than a factor of 2, for smaller sky fractions, with $τS$-lf, or as a function of foreground model. The $τS$-hf configuration does not lead to significantly tighter constraints. Exclusion of the 30 and 44 GHz data, which give information about synchrotron emission, leads to significant $τ$ mis-estimates. Decreasing noise by an ambitious factor of 10 while keeping $f_{\rm sky}$=40% gives $σ(τ) =0.0031$. The combination of $σ(τ) =0.0034$, BAO data from DESI, and future CMB B-mode lensing data from CMB-S3/S4 experiments could give $σ(\sum m_ν) = 17$ meV. △ Less

Submitted 22 November, 2022; v1 submitted 7 June, 2022; originally announced June 2022.

Comments: 12 pages, 7 Figures, reflects version published in ApJ

Journal ref: ApJ 940:68, November 2022

From BeyondPlanck to Cosmoglobe: Open Science, Reproducibility, and Data Longevity

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

Abstract: The BeyondPlanck and Cosmoglobe collaborations have implemented the first integrated Bayesian end-to-end analysis pipeline for CMB experiments. The primary long-term motivation for this work is to develop a common analysis platform that supports efficient global joint analysis of complementary radio, microwave, and sub-millimeter experiments. A strict prerequisite for this to succeed is broad part… ▽ More The BeyondPlanck and Cosmoglobe collaborations have implemented the first integrated Bayesian end-to-end analysis pipeline for CMB experiments. The primary long-term motivation for this work is to develop a common analysis platform that supports efficient global joint analysis of complementary radio, microwave, and sub-millimeter experiments. A strict prerequisite for this to succeed is broad participation from the CMB community, and two foundational aspects of the program are therefore reproducibility and Open Science. In this paper, we discuss our efforts toward this aim. We also discuss measures toward facilitating easy code and data distribution, community-based code documentation, user-friendly compilation procedures, etc. This work represents the first publicly released end-to-end CMB analysis pipeline that includes raw data, source code, parameter files, and documentation. We argue that such a complete pipeline release should be a requirement for all major future and publicly-funded CMB experiments, noting that a full public release significantly increases data longevity by ensuring that the data quality can be improved whenever better processing techniques, complementary datasets, or more computing power become available, and thereby also taxpayers' value for money; providing only raw data and final products is not sufficient to guarantee full reproducibility in the future. △ Less

Submitted 14 March, 2023; v1 submitted 20 May, 2022; originally announced May 2022.

BeyondPlanck XII. Cosmological parameter constraints with end-to-end error propagation

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

Abstract: We present cosmological parameter constraints as estimated using the Bayesian BeyondPlanck (BP) analysis framework. This method supports seamless end-to-end error propagation from raw time-ordered data to final cosmological parameters. As a first demonstration of the method, we analyze time-ordered Planck LFI observations, combined with selected external data (WMAP 33-61GHz, Planck HFI DR4 353 and… ▽ More We present cosmological parameter constraints as estimated using the Bayesian BeyondPlanck (BP) analysis framework. This method supports seamless end-to-end error propagation from raw time-ordered data to final cosmological parameters. As a first demonstration of the method, we analyze time-ordered Planck LFI observations, combined with selected external data (WMAP 33-61GHz, Planck HFI DR4 353 and 857GHz, and Haslam 408MHz) in the form of pixelized maps which are used to break critical astrophysical degeneracies. Overall, all results are generally in good agreement with previously reported values from Planck 2018 and WMAP, with the largest relative difference for any parameter of about 1 sigma when considering only temperature multipoles between 29<l<601. In cases where there are differences, we note that the BP results are generally slightly closer to the high-l HFI-dominated Planck 2018 results than previous analyses, suggesting slightly less tension between low and high multipoles. Using low-l polarization information from LFI and WMAP, we find a best-fit value of tau=0.066 +/- 0.013, which is higher than the low value of tau=0.051 +/- 0.006 derived from Planck 2018 and slightly lower than the value of 0.069 +/- 0.011 derived from joint analysis of official LFI and WMAP products. Most importantly, however, we find that the uncertainty derived in the BP processing is about 30% larger than when analyzing the official products, after taking into account the different sky coverage. We argue that this is due to marginalizing over a more complete model of instrumental and astrophysical parameters, and this results in both more reliable and more rigorously defined uncertainties. We find that about 2000 Monte Carlo samples are required to achieve robust convergence for low-resolution CMB covariance matrix with 225 independent modes. △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: 13 pages, 10 figures

MSC Class: 83F05

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

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)

Snowmass 2021 CMB-S4 White Paper

Authors: Kevork Abazajian, Arwa Abdulghafour, Graeme E. Addison, Peter Adshead, Zeeshan Ahmed, Marco Ajello, Daniel Akerib, Steven W. Allen, David Alonso, Marcelo Alvarez, Mustafa A. Amin, Mandana Amiri, Adam Anderson, Behzad Ansarinejad, Melanie Archipley, Kam S. Arnold, Matt Ashby, Han Aung, Carlo Baccigalupi, Carina Baker, Abhishek Bakshi, Debbie Bard, Denis Barkats, Darcy Barron, Peter S. Barry , et al. (331 additional authors not shown)

Abstract: This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan. This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. arXiv admin note: substantial text overlap with arXiv:1908.01062, arXiv:1907.04473

Snowmass2021 Cosmic Frontier: Cosmic Microwave Background Measurements White Paper

Authors: Clarence L. Chang, Kevin M. Huffenberger, Bradford A. Benson, Federico Bianchini, Jens Chluba, Jacques Delabrouille, Raphael Flauger, Shaul Hanany, William C. Jones, Alan J. Kogut, Jeffrey J. McMahon, Joel Meyers, Neelima Sehgal, Sara M. Simon, Caterina Umilta, Kevork N. Abazajian, Zeeshan Ahmed, Yashar Akrami, Adam J. Anderson, Behzad Ansarinejad, Jason Austermann, Carlo Baccigalupi, Denis Barkats, Darcy Barron, Peter S. Barry , et al. (107 additional authors not shown)

Abstract: This is a solicited whitepaper for the Snowmass 2021 community planning exercise. The paper focuses on measurements and science with the Cosmic Microwave Background (CMB). The CMB is foundational to our understanding of modern physics and continues to be a powerful tool driving our understanding of cosmology and particle physics. In this paper, we outline the broad and unique impact of CMB science… ▽ More This is a solicited whitepaper for the Snowmass 2021 community planning exercise. The paper focuses on measurements and science with the Cosmic Microwave Background (CMB). The CMB is foundational to our understanding of modern physics and continues to be a powerful tool driving our understanding of cosmology and particle physics. In this paper, we outline the broad and unique impact of CMB science for the High Energy Cosmic Frontier in the upcoming decade. We also describe the progression of ground-based CMB experiments, which shows that the community is prepared to develop the key capabilities and facilities needed to achieve these transformative CMB measurements. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: contribution to Snowmass 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)

Probing Cosmic Inflation with the LiteBIRD Cosmic Microwave Background Polarization Survey

Authors: LiteBIRD Collaboration, E. Allys, K. Arnold, J. Aumont, R. Aurlien, S. Azzoni, C. Baccigalupi, A. J. Banday, R. Banerji, R. B. Barreiro, N. Bartolo, L. Bautista, D. Beck, S. Beckman, M. Bersanelli, F. Boulanger, M. Brilenkov, M. Bucher, E. Calabrese, P. Campeti, A. Carones, F. J. Casas, A. Catalano, V. Chan, K. Cheung , et al. (166 additional authors not shown)

Abstract: LiteBIRD, the Lite (Light) satellite for the study of B-mode polarization and Inflation from cosmic background Radiation Detection, is a space mission for primordial cosmology and fundamental physics. The Japan Aerospace Exploration Agency (JAXA) selected LiteBIRD in May 2019 as a strategic large-class (L-class) mission, with an expected launch in the late 2020s using JAXA's H3 rocket. LiteBIRD is… ▽ More LiteBIRD, the Lite (Light) satellite for the study of B-mode polarization and Inflation from cosmic background Radiation Detection, is a space mission for primordial cosmology and fundamental physics. The Japan Aerospace Exploration Agency (JAXA) selected LiteBIRD in May 2019 as a strategic large-class (L-class) mission, with an expected launch in the late 2020s using JAXA's H3 rocket. LiteBIRD is planned to orbit the Sun-Earth Lagrangian point L2, where it will map the cosmic microwave background (CMB) polarization over the entire sky for three years, with three telescopes in 15 frequency bands between 34 and 448 GHz, to achieve an unprecedented total sensitivity of 2.2$μ$K-arcmin, with a typical angular resolution of 0.5$^\circ$ at 100 GHz. The primary scientific objective of LiteBIRD is to search for the signal from cosmic inflation, either making a discovery or ruling out well-motivated inflationary models. The measurements of LiteBIRD will also provide us with insight into the quantum nature of gravity and other new physics beyond the standard models of particle physics and cosmology. We provide an overview of the LiteBIRD project, including scientific objectives, mission and system requirements, operation concept, spacecraft and payload module design, expected scientific outcomes, potential design extensions and synergies with other projects. △ Less

Submitted 27 March, 2023; v1 submitted 6 February, 2022; originally announced February 2022.

Comments: 155 pages, accepted for publication in PTEP

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)

An Optical Spectrum of the Diffuse Galactic Light from BOSS and IRIS

Authors: Blake Chellew, Timothy D. Brandt, Brandon S. Hensley, Bruce T. Draine, Eve Matthaey

Abstract: We present a spectrum of the diffuse Galactic light (DGL) between 3700 and 10,000 A, obtained by correlating optical sky intensity with far-infrared dust emission. We use nearly 250,000 blank-sky spectra from BOSS/SDSS-III together with IRIS-reprocessed maps from the IRAS satellite. The larger sample size compared to SDSS-II results in a factor-of-two increase in signal to noise. We combine these… ▽ More We present a spectrum of the diffuse Galactic light (DGL) between 3700 and 10,000 A, obtained by correlating optical sky intensity with far-infrared dust emission. We use nearly 250,000 blank-sky spectra from BOSS/SDSS-III together with IRIS-reprocessed maps from the IRAS satellite. The larger sample size compared to SDSS-II results in a factor-of-two increase in signal to noise. We combine these data sets with a model for the optical/far-infrared correlation that accounts for self-absorption by dust. The spectral features of the DGL agree remarkably well with features present in stellar spectra. There is evidence for a difference in the DGL continuum between the regions covered by BOSS in the northern and southern Galactic hemisphere. We interpret the difference at red wavelengths as the result of a difference in stellar populations, with mainly old stars in both regions but a higher fraction of young stars in the south. There is also a broad excess in the southern DGL spectrum over the prediction of a simple radiative transfer model, without a clear counterpart in the north. We interpret this excess, centered at ~6500 A, as evidence for luminescence in the form of extended red emission (ERE). The observed strength of the 4000 A break indicates that at most ~7% of the dust-correlated light at 4000 A can be due to blue luminescence. Our DGL spectrum provides constraints on dust scattering and luminescence independent of measurements of extinction. △ Less

Submitted 4 July, 2022; v1 submitted 4 January, 2022; originally announced January 2022.

Comments: 16 pages, 8 figures. Revised version after reviewer comments

Journal ref: The Astrophysical Journal, Volume 932, Number 2. June 22, 2022

COMAP Early Science: VI. A First Look at the COMAP Galactic Plane Survey

Authors: Thomas J. Rennie, Stuart E. Harper, Clive Dickinson, Liju Philip, Kieran A. Cleary, Richard J. Bond, Jowita Borowska, Patrick C. Breysse, Morgan Catha, Roke Cepeda-Arroita, Dongwoo T. Chung, Sarah E. Church, Delaney A. Dunne, Hans Kristian Eriksen, Marie Kristine Foss, Todd Gaier, Joshua Ott Gunderson, Andrew I. Harris, Brandon Hensley, Richard Hobbs, Håvard T. Ihle, James W. Lamb, Charles R. Lawrence, Jonas G. S. Lunde, Roberta Paladini , et al. (7 additional authors not shown)

Abstract: We present early results from the COMAP Galactic Plane Survey conducted between June 2019 and April 2021, spanning $20^\circ<\ell<40^\circ$ in Galactic longitude and $|b|<1.\!\!^{\circ}5$ in Galactic latitude with an angular resolution of $4.5^{\prime}$. The full survey will span $\ell \sim 20^{\circ}$- $220^{\circ}$ and will be the first large-scale radio continuum survey at $30$ GHz with sub-deg… ▽ More We present early results from the COMAP Galactic Plane Survey conducted between June 2019 and April 2021, spanning $20^\circ<\ell<40^\circ$ in Galactic longitude and $|b|<1.\!\!^{\circ}5$ in Galactic latitude with an angular resolution of $4.5^{\prime}$. The full survey will span $\ell \sim 20^{\circ}$- $220^{\circ}$ and will be the first large-scale radio continuum survey at $30$ GHz with sub-degree resolution. We present initial results from the first part of the survey, including diffuse emission and spectral energy distributions (SEDs) of HII regions and supernova remnants. Using low and high frequency surveys to constrain free-free and thermal dust emission contributions, we find evidence of excess flux density at $30\,$GHz in six regions that we interpret as anomalous microwave emission. Furthermore we model UCHII contributions using data from the $5\,$GHz CORNISH catalogue and reject this as the cause of the $30\,$GHz excess. Six known supernova remnants (SNR) are detected at $30\,$GHz, and we measure spectral indices consistent with the literature or show evidence of steepening. The flux density of the SNR W44 at $30\,$GHz is consistent with a power-law extrapolation from lower frequencies with no indication of spectral steepening in contrast with recent results from the Sardinia Radio Telescope. We also extract five hydrogen radio recombination lines to map the warm ionized gas, which can be used to estimate electron temperatures or to constrain continuum free-free emission. The full COMAP Galactic plane survey, to be released in 2023/2024, will be an invaluable resource for Galactic astrophysics. △ Less

Submitted 21 March, 2022; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 6 of 7 in series. 28 pages, 10 figures, submitted to ApJ

COMAP Early Science: V. Constraints and Forecasts at $z \sim 3$

Authors: Dongwoo T. Chung, Patrick C. Breysse, Kieran A. Cleary, Håvard T. Ihle, Hamsa Padmanabhan, Marta B. Silva, J. Richard Bond, Jowita Borowska, Morgan Catha, Sarah E. Church, Delaney A. Dunne, Hans Kristian Eriksen, Marie Kristine Foss, Todd Gaier, Joshua Ott Gundersen, Stuart E. Harper, Andrew I. Harris, Brandon Hensley, Richard Hobbs, Laura C. Keating, Junhan Kim, James W. Lamb, Charles R. Lawrence, Jonas Gahr Sturtzel Lunde, Norman Murray , et al. (12 additional authors not shown)

Abstract: We present the current state of models for the $z\sim3$ carbon monoxide (CO) line-intensity signal targeted by the CO Mapping Array Project (COMAP) Pathfinder in the context of its early science results. Our fiducial model, relating dark matter halo properties to CO luminosities, informs parameter priors with empirical models of the galaxy-halo connection and previous CO(1-0) observations. The Pat… ▽ More We present the current state of models for the $z\sim3$ carbon monoxide (CO) line-intensity signal targeted by the CO Mapping Array Project (COMAP) Pathfinder in the context of its early science results. Our fiducial model, relating dark matter halo properties to CO luminosities, informs parameter priors with empirical models of the galaxy-halo connection and previous CO(1-0) observations. The Pathfinder early science data spanning wavenumbers $k=0.051$-$0.62\,$Mpc$^{-1}$ represent the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum. Our 95% upper limit on the redshift-space clustering amplitude $A_{\rm clust}\lesssim70\,μ$K$^2$ greatly improves on the indirect upper limit of $420\,μ$K$^2$ reported from the CO Power Spectrum Survey (COPSS) measurement at $k\sim1\,$Mpc$^{-1}$. The COMAP limit excludes a subset of models from previous literature, and constrains interpretation of the COPSS results, demonstrating the complementary nature of COMAP and interferometric CO surveys. Using line bias expectations from our priors, we also constrain the squared mean line intensity-bias product, $\langle{Tb}\rangle^2\lesssim50\,μ$K$^2$, and the cosmic molecular gas density, $ρ_\text{H2}<2.5\times10^8\,M_\odot\,$Mpc$^{-3}$ (95% upper limits). Based on early instrument performance and our current CO signal estimates, we forecast that the five-year Pathfinder campaign will detect the CO power spectrum with overall signal-to-noise of 9-17. Between then and now, we also expect to detect the CO-galaxy cross-spectrum using overlapping galaxy survey data, enabling enhanced inferences of cosmic star-formation and galaxy-evolution history. △ Less

Submitted 4 March, 2022; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 5 of 7 in series. 17 pages + appendix and bibliography (30 pages total); 15 figures, 6 tables; accepted for publication in ApJ; v3 reflects the accepted version with minor changes and additions to text

Journal ref: ApJ, 933, 186 (2022)

COMAP Early Science: I. Overview

Authors: Kieran A. Cleary, Jowita Borowska, Patrick C. Breysse, Morgan Catha, Dongwoo T. Chung, Sarah E. Church, Clive Dickinson, Hans Kristian Eriksen, Marie Kristine Foss, Joshua Ott Gundersen, Stuart E. Harper, Andrew I. Harris, Richard Hobbs, Håvard, T. Ihle, Junhan Kim, Jonathon Kocz, James W. Lamb, Jonas G. S. Lunde, Hamsa Padmanabhan, Timothy J. Pearson, Liju Philip, Travis W. Powell, Maren Rasmussen, Anthony C. S. Readhead , et al. (18 additional authors not shown)

Abstract: The CO Mapping Array Project (COMAP) aims to use line intensity mapping of carbon monoxide (CO) to trace the distribution and global properties of galaxies over cosmic time, back to the Epoch of Reionization (EoR). To validate the technologies and techniques needed for this goal, a Pathfinder instrument has been constructed and fielded. Sensitive to CO(1-0) emission from $z=2.4$-$3.4$ and a fainte… ▽ More The CO Mapping Array Project (COMAP) aims to use line intensity mapping of carbon monoxide (CO) to trace the distribution and global properties of galaxies over cosmic time, back to the Epoch of Reionization (EoR). To validate the technologies and techniques needed for this goal, a Pathfinder instrument has been constructed and fielded. Sensitive to CO(1-0) emission from $z=2.4$-$3.4$ and a fainter contribution from CO(2-1) at $z=6$-8, the Pathfinder is surveying $12$ deg$^2$ in a 5-year observing campaign to detect the CO signal from $z\sim3$. Using data from the first 13 months of observing, we estimate $P_\mathrm{CO}(k) = -2.7 \pm 1.7 \times 10^4μ\mathrm{K}^2 \mathrm{Mpc}^3$ on scales $k=0.051-0.62 \mathrm{Mpc}^{-1}$ - the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum. Based on these observations alone, we obtain a constraint on the amplitude of the clustering component (the squared mean CO line temperature-bias product) of $\langle Tb\rangle^2<49$ $μ$K$^2$ - nearly an order-of-magnitude improvement on the previous best measurement. These constraints allow us to rule out two models from the literature. We forecast a detection of the power spectrum after 5 years with signal-to-noise ratio (S/N) 9-17. Cross-correlation with an overlapping galaxy survey will yield a detection of the CO-galaxy power spectrum with S/N of 19. We are also conducting a 30 GHz survey of the Galactic plane and present a preliminary map. Looking to the future of COMAP, we examine the prospects for future phases of the experiment to detect and characterize the CO signal from the EoR. △ Less

Submitted 29 November, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 1 of 7 in series. 18 pages, 16 figures, submitted to ApJ

Polycyclic Aromatic Hydrocarbons, the Anomalous Microwave Emission, and Their Connection to the Cold Neutral Medium

Authors: Brandon S. Hensley, Claire E. Murray, Mark Dodici

Abstract: Using new large area maps of the cold neutral medium (CNM) fraction, $f_{\rm CNM}$, we investigate the relationship between the CNM, the abundance of polycyclic aromatic hydrocarbons (PAHs), and the anomalous microwave emission (AME). We first present our $f_{\rm CNM}$ map based on full-sky HI4PI data, using a convolutional neural network to covert the spectroscopic HI data to $f_{\rm CNM}$. We de… ▽ More Using new large area maps of the cold neutral medium (CNM) fraction, $f_{\rm CNM}$, we investigate the relationship between the CNM, the abundance of polycyclic aromatic hydrocarbons (PAHs), and the anomalous microwave emission (AME). We first present our $f_{\rm CNM}$ map based on full-sky HI4PI data, using a convolutional neural network to covert the spectroscopic HI data to $f_{\rm CNM}$. We demonstrate that $f_{\rm CNM}$ is strongly correlated with the fraction of dust in PAHs as estimated from mid- and far-infrared dust emission. In contrast, we find no correlation between $f_{\rm CNM}$ and the amount of AME per dust emission, nor between $f_{\rm CNM}$ and the AME peak frequency. These results suggest PAHs preferentially reside in cold, relatively dense gas, perhaps owing to enhanced destruction in more diffuse media. The lack of correlation between $f_{\rm CNM}$ and AME peak frequency is in tension with expectations from theoretical models positing different spectral energy distributions of AME in the cold versus warm neutral medium. We suggest that different PAH abundances and emission physics in different interstellar environments may explain the weaker than expected correlation between 12$μ$m PAH emission and AME even if PAHs are the AME carriers. △ Less

Submitted 4 November, 2021; originally announced November 2021.

Comments: Submitted to ApJ. 12 pages, 4 figures, comments welcome

The Simons Observatory: Galactic Science Goals and Forecasts

Authors: Brandon S. Hensley, Susan E. Clark, Valentina Fanfani, Nicoletta Krachmalnicoff, Giulio Fabbian, Davide Poletti, Giuseppe Puglisi, Gabriele Coppi, Jacob Nibauer, Roman Gerasimov, Nicholas Galitzki, Steve K. Choi, Peter C. Ashton, Carlo Baccigalupi, Eric Baxter, Blakesley Burkhart, Erminia Calabrese, Jens Chluba, Josquin Errard, Andrei V. Frolov, Carlos Hervías-Caimapo, Kevin M. Huffenberger, Bradley R. Johnson, Baptiste Jost, Brian Keating , et al. (9 additional authors not shown)

Abstract: Observing in six frequency bands from 27 to 280 GHz over a large sky area, the Simons Observatory (SO) is poised to address many questions in Galactic astrophysics in addition to its principal cosmological goals. In this work, we provide quantitative forecasts on astrophysical parameters of interest for a range of Galactic science cases. We find that SO can: constrain the frequency spectrum of pol… ▽ More Observing in six frequency bands from 27 to 280 GHz over a large sky area, the Simons Observatory (SO) is poised to address many questions in Galactic astrophysics in addition to its principal cosmological goals. In this work, we provide quantitative forecasts on astrophysical parameters of interest for a range of Galactic science cases. We find that SO can: constrain the frequency spectrum of polarized dust emission at a level of $Δβ_d \lesssim 0.01$ and thus test models of dust composition that predict that $β_d$ in polarization differs from that measured in total intensity; measure the correlation coefficient between polarized dust and synchrotron emission with a factor of two greater precision than current constraints; exclude the non-existence of exo-Oort clouds at roughly 2.9$σ$ if the true fraction is similar to the detection rate of giant planets; map more than 850 molecular clouds with at least 50 independent polarization measurements at 1 pc resolution; detect or place upper limits on the polarization fractions of CO(2-1) emission and anomalous microwave emission at the 0.1% level in select regions; and measure the correlation coefficient between optical starlight polarization and microwave polarized dust emission in $1^\circ$ patches for all lines of sight with $N_{\rm H} \gtrsim 2\times10^{20}$ cm$^{-2}$. The goals and forecasts outlined here provide a roadmap for other microwave polarization experiments to expand their scientific scope via Milky Way astrophysics. △ Less

Submitted 3 November, 2021; originally announced November 2021.

Comments: Submitted to AAS journals. 33 pages, 10 figures

The Galaxy Evolution Probe

Authors: Jason Glenn, Charles M. Bradford, Erik Rosolowsky, Rashied Amini, Katherine Alatalo, Lee Armus, Andrew J. Benson, Tzu-Ching Chang, Jeremy Darling, Peter K. Day, Jeanette Domber, Duncan Farrah, Brandon Hensley, Sarah Lipscy, Bradley Moore, Seb Oliver, Joanna Perido, David Redding, Michael Rodgers, Raphael Shirley, Howard A. Smith, John B. Steeves, Carole Tucker, Jonas Zmuidzinas

Abstract: The Galaxy Evolution Probe (GEP) is a concept for a mid- and far-infrared space observatory to measure key properties of large samples of galaxies with large and unbiased surveys. GEP will attempt to achieve zodiacal light and Galactic dust emission photon background-limited observations by utilizing a 6 Kelvin, 2.0 meter primary mirror and sensitive arrays of kinetic inductance detectors. It will… ▽ More The Galaxy Evolution Probe (GEP) is a concept for a mid- and far-infrared space observatory to measure key properties of large samples of galaxies with large and unbiased surveys. GEP will attempt to achieve zodiacal light and Galactic dust emission photon background-limited observations by utilizing a 6 Kelvin, 2.0 meter primary mirror and sensitive arrays of kinetic inductance detectors. It will have two instrument modules: a 10 - 400 micron hyperspectral imager with spectral resolution R = 8 (GEP-I) and a 24 - 193 micron, R = 200 grating spectrometer (GEP-S). GEP-I surveys will identify star-forming galaxies via their thermal dust emission and simultaneously measure redshifts using polycyclic aromatic hydrocarbon emission lines. Galaxy luminosities derived from star formation and nuclear supermassive black hole accretion will be measured for each source, enabling the cosmic star formation history to be measured to much greater precision than previously possible. Using optically thin far-infrared fine-structure lines, surveys with GEP-S will measure the growth of metallicity in the hearts of galaxies over cosmic time and extraplanar gas will be mapped in spiral galaxies in the local universe to investigate feedback processes. The science case and mission architecture designed to meet the science requirements are described, and the kinetic inductance detector and readout electronics state of the art and needed developments are described. This paper supersedes the GEP concept study report cited in it by providing new content, including: a summary of recent mid-infrared KID development, a discussion of microlens array fabrication for mid-infrared KIDs, and additional context for galaxy surveys. The reader interested in more technical details may want to consult the concept study report. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: 84 pages, 24 figures, published in JATIS

Journal ref: J. Astron. Telesc. Instrum. Syst. 7(3), 034004 (2021)

CCAT-prime Collaboration: Science Goals and Forecasts with Prime-Cam on the Fred Young Submillimeter Telescope

Authors: CCAT-Prime collaboration, M. Aravena, J. E. Austermann, K. Basu, N. Battaglia, B. Beringue, F. Bertoldi, F. Bigiel, J. R. Bond, P. C. Breysse, C. Broughton, R. Bustos, S. C. Chapman, M. Charmetant, S. K. Choi, D. T. Chung, S. E. Clark, N. F. Cothard, A. T. Crites, A. Dev, K. Douglas, C. J. Duell, R. Dunner, H. Ebina, J. Erler , et al. (62 additional authors not shown)

Abstract: We present a detailed overview of the science goals and predictions for the Prime-Cam direct detection camera/spectrometer being constructed by the CCAT-prime collaboration for dedicated use on the Fred Young Submillimeter Telescope (FYST). The FYST is a wide-field, 6-m aperture submillimeter telescope being built (first light in mid-2024) by an international consortium of institutions led by Corn… ▽ More We present a detailed overview of the science goals and predictions for the Prime-Cam direct detection camera/spectrometer being constructed by the CCAT-prime collaboration for dedicated use on the Fred Young Submillimeter Telescope (FYST). The FYST is a wide-field, 6-m aperture submillimeter telescope being built (first light in mid-2024) by an international consortium of institutions led by Cornell University and sited at more than 5600 meters on Cerro Chajnantor in northern Chile. Prime-Cam is one of two instruments planned for FYST and will provide unprecedented spectroscopic and broadband measurement capabilities to address important astrophysical questions ranging from Big Bang cosmology through reionization and the formation of the first galaxies to star formation within our own Milky Way galaxy. Prime-Cam on the FYST will have a mapping speed that is over ten times greater than existing and near-term facilities for high-redshift science and broadband polarimetric imaging at frequencies above 300 GHz. We describe details of the science program enabled by this system and our preliminary survey strategies. △ Less

Submitted 8 August, 2022; v1 submitted 21 July, 2021; originally announced July 2021.

Comments: 61 pages, 16 figures. Resubmitted to ApJSS July 11, 2022

A study of 90 GHz dust emissivity on molecular cloud and filament scales

Authors: Ian Lowe, Brian Mason, Tanay Bhandarkar, S. E. Clark, Mark Devlin, Simon R. Dicker, Shannon M. Duff, Rachel Friesen, Alvaro Hacar, Brandon Hensley, Tony Mroczkowski, Sigurd Naess, Charles Romero, Sarah Sadavoy, Maria Salatino, Craig Sarazin, John Orlowski-Scherer, Alessandro Schillaci, Jonathan Sievers, Thomas Stanke, Amelia Stutz, Zhilei Xu

Abstract: Recent observations from the MUSTANG2 instrument on the Green Bank Telescope have revealed evidence of enhanced long-wavelength emission in the dust spectral energy distribution (SED) in the Orion Molecular Cloud (OMC) 2/3 filament on 25" (0.1 pc) scales. Here we present a measurement of the SED on larger spatial scales (map size 0.5-3 degrees or 3-20 pc), at somewhat lower resolution (120", corre… ▽ More Recent observations from the MUSTANG2 instrument on the Green Bank Telescope have revealed evidence of enhanced long-wavelength emission in the dust spectral energy distribution (SED) in the Orion Molecular Cloud (OMC) 2/3 filament on 25" (0.1 pc) scales. Here we present a measurement of the SED on larger spatial scales (map size 0.5-3 degrees or 3-20 pc), at somewhat lower resolution (120", corresponding to 0.25 pc at 400 pc) using data from the Herschel satellite and Atacama Cosmology Telescope (ACT). We then extend the 120"-scale investigation to other regions covered in the Herschel Gould Belt Survey (HGBS) specifically: the dense filaments in the southerly regions of Orion A; Orion B; and Serpens-S. Our dataset in aggregate covers approximately 10 square degrees, with continuum photometry spanning from 160um to 3mm. These OMC 2/3 data display excess emission at 3mm, though less (10.9% excess) than what is seen at higher resolution. Strikingly, we find that the enhancement is present even more strongly in the other filaments we targeted, with an average excess of 42.4% and 30/46 slices showing an inconsistency with the modified blackbody to at least 4σ. Applying this analysis to the other targeted regions, we lay the groundwork for future high-resolution analyses. Additionally, we also consider a two-component dust model motivated by Planck results and an amorphous grain dust model. While both of these have been proposed to explain deviations in emission from a generic modified blackbody (MBB), we find that they have significant drawbacks, requiring many spectral points or lacking experimental data coverage. △ Less

Submitted 14 July, 2022; v1 submitted 27 May, 2021; originally announced May 2021.

Journal ref: The Astrophysical Journal 929 (2022) 102

The Atacama Cosmology Telescope: Microwave Intensity and Polarization Maps of the Galactic Center

Authors: Yilun Guan, Susan E. Clark, Brandon S. Hensley, Patricio A. Gallardo, Sigurd Naess, Cody J. Duell, Simone Aiola, Zachary Atkins, Erminia Calabrese, Steve K. Choi, Nicholas F. Cothard, Mark Devlin, Adriaan J. Duivenvoorden, Jo Dunkley, Rolando Dünner, Simone Ferraro, Matthew Hasselfield, John P. Hughes, Brian J. Koopman, Arthur B. Kosowsky, Mathew S. Madhavacheril, Jeff McMahon, Federico Nati, Michael D. Niemack, Lyman A. Page , et al. (8 additional authors not shown)

Abstract: We present arcminute-resolution intensity and polarization maps of the Galactic center made with the Atacama Cosmology Telescope (ACT). The maps cover a 32 deg$^2$ field at 98, 150, and 224 GHz with $\vert l\vert\le4^\circ$, $\vert b\vert\le2^\circ$. We combine these data with Planck observations at similar frequencies to create coadded maps with increased sensitivity at large angular scales. With… ▽ More We present arcminute-resolution intensity and polarization maps of the Galactic center made with the Atacama Cosmology Telescope (ACT). The maps cover a 32 deg$^2$ field at 98, 150, and 224 GHz with $\vert l\vert\le4^\circ$, $\vert b\vert\le2^\circ$. We combine these data with Planck observations at similar frequencies to create coadded maps with increased sensitivity at large angular scales. With the coadded maps, we are able to resolve many known features of the Central Molecular Zone (CMZ) in both total intensity and polarization. We map the orientation of the plane-of-sky component of the Galactic magnetic field inferred from the polarization angle in the CMZ, finding significant changes in morphology in the three frequency bands as the underlying dominant emission mechanism changes from synchrotron to dust emission. Selected Galactic center sources, including Sgr A*, the Brick molecular cloud (G0.253+0.016), the Mouse pulsar wind nebula (G359.23-0.82), and the Tornado supernova remnant candidate (G357.7-0.1), are examined in detail. These data illustrate the potential for leveraging ground-based Cosmic Microwave Background polarization experiments for Galactic science. △ Less

Submitted 14 September, 2021; v1 submitted 11 May, 2021; originally announced May 2021.

Comments: 26 pages, 14 figures, accepted for publication in ApJ