Cosmic Microwave Background Polarization Research Articles

Measuring the reionization optical depth without large-scale CMB polarization

We study the possibility of measuring the optical depth at reionization τ without relying on large-scale cosmic microwave background (CMB) polarization. Our analysis is driven by the need to obtain competitive measurements that can validate the state-of-the-art constraints on this parameter, widely based on E-mode polarization measurements at ℓ≤30. This need is partially motivated by the typical concerns regarding anomalies observed in the Planck large-scale CMB data as well as by the remarkable fact that, excluding these latter, τ consistently exhibits correlations with anomalous parameters, such as Alens and Ωk, suggesting that slightly higher values of the optical depth at reionization could significantly alleviate or even eliminate anomalies. Within the Λ cold dark matter model, our most constraining result is τ=0.080±0.012, obtained by combining Planck temperature and polarization data at ℓ>30, the Atacama Cosmology Telescope (ACT) and Planck measurements of the lensing potential, baryon acoustic oscillations (BAOs), and type-Ia supernova data from the Pantheon+catalog. Notably, using only ACT temperature, polarization, and lensing data in combination with BAOs and supernovae, we obtain τ=0.076±0.015, which is entirely independent of Planck. The relative precision of these results is approaching the constraints based on large-scale CMB polarization (τ=0.054±0.008). Despite the overall agreement, we report a slight 1.8σ shift toward larger values of τ. We also test how these results change by extending the cosmological model. While in many extensions they remain robust, in general, obtaining precise measurements of τ may become significantly more challenging. Published by the American Physical Society 2024

Deployment of POLARBEAR-2b

POLARBEAR-2b (PB-2b) is the second receiver in the Simons Array, a cosmic microwave background (CMB) polarization experiment. The Simons Array uses dichroic polarization sensitive lenslet-coupled sinuous antennas and transition-edge sensor (TES) bolometers made of superconducting films. These bolometers are read out with frequency multiplexing electronics. PB-2b contains ∼\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim$$\\end{document} 7500 detectors in two bands at 90 and 150 GHz with arcminute resolution. The polarization of these detectors is modulated by a cryogenic continuously rotating half-wave plate. PB-2b was installed on its telescope in 2022 in the Atacama Desert at an altitude of 5.2 km. This paper will detail initial readout commissioning, test of a new loopgain monitoring method, and focusing the optics. Work is ongoing to commission the remaining ambient temperature readout electronics, measure detector time constants, and observe with the cryogenic half-wave plate spinning

High-accuracy emulators for observables in ΛCDM, Neff, Σmν, and w cosmologies

ABSTRACT We use the emulation framework CosmoPower to construct and publicly release neural network emulators of cosmological observables, including the cosmic microwave background (CMB) temperature and polarization power spectra, matter power spectrum, distance-redshift relation, baryon acoustic oscillation (BAO) and redshift-space distortion (RSD) observables, and derived parameters. We train our emulators on Einstein–Boltzmann calculations obtained with high-precision numerical convergence settings, for a wide range of cosmological models including ΛCDM, wCDM, ΛCDM + Neff, and ΛCDM + Σmν. Our CMB emulators are accurate to better than 0.5 per cent out to ℓ = 104, which is sufficient for Stage-IV data analysis, and our P(k) emulators reach the same accuracy level out to $k=50 \, \, \mathrm{Mpc}^{-1}$, which is sufficient for Stage-III data analysis. We release the emulators via an online repository (CosmoPower Organisation), which will be continually updated with additional extended cosmological models. Our emulators accelerate cosmological data analysis by orders of magnitude, enabling cosmological parameter extraction analyses, using current survey data, to be performed on a laptop. We validate our emulators by comparing them to class and camb and by reproducing cosmological parameter constraints derived from Planck TT, TE, EE, and CMB lensing data, as well as from the Atacama Cosmology Telescope Data Release 4 CMB data, Dark Energy Survey Year-1 galaxy lensing and clustering data, and Baryon Oscillation Spectroscopic Survey Data Release 12 BAO and RSD data.

Implementation and testing of a 100 mK ADR backed by a 4He/3He sorption fridge for CMB-S4 project

CMB-S4 is a ground-based experiment that will use large-format bolometer arrays to map Cosmic Microwave Background (CMB) polarization with unprecedented sensitivity, including the search for a B-mode polarization pattern associated with cosmic inflation. The focal planes on which the superconducting detectors are mounted weigh about 20 kg each and must be cooled at 100 mK. To reach this temperature for a first receiver fielding prototype hardware, we have developed an Adiabatic Demagnetization Refrigerator (ADR) used in conjunction with an existing 4He/3He sorption refrigerator.The ADR unit described here is based on a Chromium Potassium Alum (CPA) salt pill. It is backed by a 3-stage (4He/3He/3He) sorption fridge, coupled through a 3He gas-gap heat switch. This heat switch is connected to a 3He stage of the sorption fridge to dissipate the heat generated by magnetization. The ADR is cycled with a 1.3 T magnetic field. The parasitic field at the detectors’ location during observations is kept below 2.5 µT thanks to a ferromagnetic shield. The pill is supported by Kevlar lines, which are doubly intercepted by the 4He and 3He stages of the sorption fridge to reduce conductive heat losses. With a total recycling time of 12 hours, this cryogenic assembly provides 2 µW of cooling power during 48 hours at 100 mK. The ADR has been designed by CEA/DSBT in France and is now implemented in a cryostat at Caltech.

B-mode polarization forecasts for GreenPol

We present tensor-to-scalar ratio forecasts for GreenPol, a hypothetical ground-based B-mode experiment aiming to survey the cleanest regions of the Northern Galactic Hemisphere at five frequencies between 10 and 44 GHz. Its primary science goal would be to measure large-scale cosmic microwave background (CMB) polarization fluctuations at multipoles ℓ ≲ 500, and thereby constrain the primordial tensor-to-scalar ratio r. The observations for the suggested experiment would take place at the Summit Station (72 ° 34N, 38 ° 27W) on Greenland, at an altitude of 3216 m above sea level. For this paper we simulated various experimental setups, and derived limits on the tensor-to-scalar ratio after CMB component separation using a Bayesian component separation implementation called Commander. When combining the proposed experiment with Planck HFI observations for constraining polarized thermal dust emission, we found a projected limit of r < 0.02 at 95% confidence for the baseline configuration. This limit is very robust with respect to a range of important experimental parameters, including sky coverage, detector weighting, foreground priors, among others. Overall, GreenPol would have the possibility to provide deep CMB polarization measurements of the Northern Galactic Hemisphere at low frequencies.

THEMIS 2.0: A self-consistent model for dust extinction, emission, and polarisation

Context. Recent observational constraints in emission, extinction, and polarisation have at least partially invalidated most of the astronomical standard grain models for the diffuse interstellar medium. Moreover, laboratory measurements on interstellar silicate analogues have shown quite significant differences with the optical properties used in these standard models. Aims. To address these issues, our objective is twofold: (i) to update the optical properties of silicates and (ii) to develop The Heterogeneous dust Evolution Model for Interstellar Solids (THEMIS) to allow the calculation of polarised extinction and emission. Methods. Based on optical constants measured in the laboratory from 5 µm to 1 mm for amorphous silicates and on observational constraints in mid-IR extinction and X-ray scattering, we defined new optical constants for the THEMIS silicates. Absorption and scattering efficiencies for spheroidal grains using these properties were subsequently derived with the discrete dipole approximation. Results. These new optical properties make it possible to explain the dust emission and extinction, both total and polarised. It is noteworthy that the model is not yet pushed to its limits since it does not require the perfect alignment of all grains to explain the observations and it therefore has the potential to accommodate the highest polarisation levels inferred from extinction measurements. Moreover, the dispersion of the optical properties of the different silicates measured in the laboratory naturally explain the variations in both the total and polarised emission and extinction observed in the diffuse interstellar medium. Conclusions. A single, invariant model calibrated on one single set of observations is obsolete for explaining contemporary observations. We are proposing a completely flexible dust model based entirely on laboratory measurements that has the potential to make major advances in understanding the exact nature of interstellar grains and how they evolve as a function of their radiative and dynamic environment. Even if challenging, this is also relevant for future cosmic microwave background (CMB) missions that will aim to perform precise measurements of the CMB spectral distortions and polarisation.

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

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 the polarization of starlight in a region of the sky of about four square degrees. We propose a Bayesian method to decompose the stellar-polarization source field along the distance to invert the three-dimensional (3D) volume occupied by the observed stars. We used this method 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 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) ≈ (103.3°, 22.3°), 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 and which we identified as being an intersection of the wall of the Local Bubble and the Cepheus Flare. Farther clouds, with a distance of up to 2 kpc, were similarly detected. Some of them likely correspond to intermediate-velocity clouds seen in H I spectra in this region of the sky. We found that the orientation of the plane-of-sky 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 plane-of-sky component of the magnetic field in dusty regions. 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.

CLASS Angular Power Spectra and Map-component Analysis for 40 GHz Observations through 2022

Measurement of the largest angular scale (ℓ < 30) features of the cosmic microwave background (CMB) polarization is a powerful way to constrain the optical depth to reionization and search for the signature of inflation through the detection of primordial B-modes. We present an analysis of maps covering 73.6% of the sky made from the 40 GHz channel of the Cosmology Large Angular Scale Surveyor (CLASS) from 2016 August to 2022 May. Taking advantage of the measurement stability enabled by front-end polarization modulation and excellent conditions from the Atacama Desert, we show this channel achieves higher sensitivity than the analogous frequencies from satellite measurements in the range 10 < ℓ < 100. Simulations show the CLASS linear (circular) polarization maps have a white noise level of 125(130)μKarcmin . We measure the Galaxy-masked EE and BB spectra of diffuse synchrotron radiation and compare to space-based measurements at similar frequencies. In combination with external data, we expand measurements of the spatial variations of the synchrotron spectral energy density (SED) to include new sky regions and measure the diffuse SED in the harmonic domain. We place a new upper limit on a background of circular polarization in the range 5 < ℓ < 125 with the first bin showing D ℓ < 0.023 μKCMB2 at 95% confidence. These results establish a new standard for recovery of the largest-scale CMB polarization from the ground and signal exciting possibilities when the higher sensitivity and higher-frequency CLASS channels are included in the analysis.

Testing multi-field inflation with LiteBIRD

We investigate expected constraints on the primordial tensor power spectrum from the future cosmic microwave background polarization experiment LiteBIRD as a test of multi-field inflation, where we specifically consider spectator models as representative examples. We argue that the measurements of the tensor-to-scalar ratio and the tensor spectral index, in combination with the constraints on the scalar spectral index from the Planck observation, are useful in testing multi-field inflation models. We also discuss implications for multi-field inflationary model building.

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

We present new measurements of cosmic microwave background (CMB) lensing over 9400 deg2 of the sky. These lensing measurements are derived from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) CMB data set, 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 that our results are robust, we analyze an extensive set of null tests, consistency tests, and systematic error estimates and employ a blinded analysis framework. Our CMB lensing power spectrum measurement provides constraints on the amplitude of cosmic structure that do not depend on Planck or galaxy survey data, thus giving independent information about large-scale structure growth and potential tensions in structure measurements. The baseline spectrum is well fit by a lensing amplitude of A lens = 1.013 ± 0.023 relative to the Planck 2018 CMB power spectra best-fit ΛCDM model and A lens = 1.005 ± 0.023 relative to the ACT DR4 + WMAP best-fit model. From our lensing power spectrum measurement, we derive constraints on the parameter combination S8CMBL≡σ8Ωm/0.30.25 of S8CMBL=0.818±0.022 from ACT DR6 CMB lensing alone and S8CMBL=0.813±0.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 ACT DR4 + WMAP CMB power spectrum measurements. Our lensing measurements from redshifts z ∼ 0.5–5 are thus fully consistent with ΛCDM structure growth predictions based on CMB anisotropies probing primarily z ∼ 1100. We find no evidence for a suppression of the amplitude of cosmic structure at low redshifts.

The Simons Observatory: Cryogenic half wave plate rotation mechanism for the small aperture telescopes.

We present the requirements, design, and evaluation of the cryogenic continuously rotating half-wave plate (CHWP) for the Simons Observatory (SO). SO is a cosmic microwave background polarization experiment at Parque Astronómico de Atacama in northern Chile that covers a wide range of angular scales using both small (⌀0.42m) and large (⌀6m) aperture telescopes. In particular, the small aperture telescopes (SATs) focus on large angular scales for primordial B-mode polarization. To this end, the SATs employ a CHWP to modulate the polarization of the incident light at 8Hz, suppressing atmospheric 1/f noise and mitigating systematic uncertainties that would otherwise arise due to the differential response of detectors sensitive to orthogonal polarizations. The CHWP consists of a 505mm diameter achromatic sapphire HWP and a cryogenic rotation mechanism, both of which are cooled down to ∼50K to reduce detector thermal loading. Under normal operation, the HWP is suspended by a superconducting magnetic bearing and rotates with a constant 2Hz frequency, controlled by an electromagnetic synchronous motor. We find that the number of superconductors and the number of magnets that make up the superconducting magnetic bearing are important design parameters, especially for the rotation mechanism's vibration performance. The rotation angle is detected through an angular encoder with a noise level of 0.07 μrad s. During a cooldown process, the rotor is held in place by a grip-and-release mechanism that serves as both an alignment device and a thermal path. In this paper, we provide an overview of the SO SAT CHWP: its requirements, hardware design, and laboratory performance.

Patchy screening of the CMB from dark photons

We study anisotropic (patchy) screening induced by the resonant conversion of cosmic microwave background (CMB) photons into dark-sector massive vector bosons (dark photons) as they cross non-linear large scale structure (LSS). Resonant conversion takes place through the kinetic mixing of the photon with the dark photon, one of the simplest low energy extensions to the Standard Model. In the early Universe, resonant conversion can occur when the photon plasma mass, obtained as the photon propagates through the ionized interstellar and intergalactic media, matches the dark photon mass. After the epoch of reionization, resonant conversion occurs mainly in the ionized gas that occupies virialized dark matter halos, for a range of dark photon masses between 10-13 eV ≲ m A' ≲ 10-11 eV. This leads to new CMB anisotropies that are correlated with LSS, which we refer to as patchy dark screening, in analogy with anisotropies from Thomson screening. Its unique frequency dependence allows it to be distinguished from the blackbody CMB. In this paper, we use a halo model approach to predict the imprint of dark screening on the CMB temperature and polarization anisotropies, as well as their correlation with LSS. We then examine the two- and three-point correlation functions of the dark-screened CMB, as well as correlation functions between CMB and LSS observables, to project the sensitivity of future measurements to the kinetic mixing parameter and dark photon mass. We demonstrate that an analysis with existing CMB data can improve upon current constraints on the kinetic mixing parameter by two orders of magnitude with the two-point correlation functions, while data from upcoming CMB experiments and LSS surveys can further improve the reach by another order of magnitude with two- and three-point correlation functions.

Minkowski Functionals in 𝖲𝖮(3) for the spin-2 CMB polarisation field

The study of the angular power spectrum of Cosmic Microwave Background (CMB) anisotropies, both in intensity and in polarisation, has led to the tightest constraints on cosmological parameters. However, this statistical quantity is not sensitive to any deviation from Gaussianity and statistical isotropy in the CMB data. Minkowski Functionals (MFs) have been adopted as one of the most powerful statistical tools to study such deviations, since they characterise the topology and geometry of the field of interest. In this paper, we extend the application of MFs to CMB polarisation data by introducing a new formalism, where we lift the spin 2 polarisation field to a scalar function in a higher-dimensional manifold: the group of rotations of the sphere, SO(3). Such a function is defined as f = Q cos(2ζ) - U sin(2ζ). We analytically obtain the expected values for the MFs of f in the case of Gaussian isotropic polarisation maps. Furthermore, we present a new pipeline which estimates these MFs from input HEALPix polarisation maps. We apply it to CMB simulations in order to validate the theoretical results and the methodology. The pipeline is to be included in the publicly available Python package Pynkowski .

The Simons Observatory: A fully remote controlled calibration system with a sparse wire grid for cosmic microwave background telescopes.

For cosmic microwave background (CMB) polarization observations, calibration of detector polarization angles is essential. We have developed a fully remote controlled calibration system with a sparse wire grid that reflects linearly polarized light along the wire direction. The new feature is a remote-controlled system for regular calibration, which has not been possible in sparse wire grid calibrators in past experiments. The remote control can be achieved by two electric linear actuators that load or unload the sparse wire grid into a position centered on the optical axis of a telescope between the calibration time and CMB observation. Furthermore, the sparse wire grid can be rotated by using a motor. A rotary encoder and a gravity sensor are installed on the sparse wire grid to monitor the wire direction. They allow us to achieve detector polarization angle calibration with an expected systematic error of 0.08°. The calibration system will be installed in small-aperture telescopes at Simons Observatory.

Tracking the multifield dynamics with cosmological data: a Monte Carlo approach

We introduce a numerical method specifically designed for investigating generic multifield models of inflation where a number of scalar fields ϕ K are minimally coupled to gravity and live in a field space with a non-trivial metric 𝒢> IJ (ϕ K ). Our algorithm consists of three main parts. Firstly, we solve the field equations through the entire inflationary period, deriving predictions for observable quantities such as the spectrum of scalar perturbations, primordial gravitational waves, and isocurvature modes. We also incorporate the transfer matrix formalism to track the behavior of adiabatic and isocurvature modes on super-horizon scales and the transfer of entropy to scalar modes after the horizon crossing. Secondly, we interface our algorithm with Boltzmann integrator codes to compute the subsequent full cosmology, including the cosmic microwave background anisotropies and polarization angular power spectra. Finally, we develop a novel sampling algorithm able to efficiently explore a large volume of the parameter space and identify a sub-region where theoretical predictions agree with observations. In this way, sampling over the initial conditions of the fields and the free parameters of the models, we enable Monte Carlo analysis of multifield scenarios. We test all the features of our approach by analyzing a specific model and deriving constraints on its free parameters. Our methodology provides a robust framework for studying multifield inflation, opening new avenues for future research in the field.

Hemispherical power asymmetry in intensity and polarization for Planck PR4 data

One of the foundations of the Standard Model of Cosmology is statistical isotropy, which can be tested, among other probes, through the study of the Cosmic Microwave Background (CMB). However, a hemispherical power asymmetry on large scales has been reported for WMAP and Planck data by different works.The statistical significance is above 3σ for temperature, suggesting a directional dependence of the local power spectrum, and thus a feature beyond the ΛCDM model. With the third release of the Planck data (PR3), a new analysis was performed including the E-mode polarization maps, finding an asymmetry at a modest level of significance. In this work, we perform an asymmetry analysis in intensity and polarization maps for the latest Planck processing pipeline (PR4). We obtain similar results to those obtained with PR3, with a slightly lower significance (2.8% for the Sevem method)for the amplitude of the E-mode local variance dipole as well as a significant variability with the considered mask.In addition, a hint of a possible T-E alignment between the asymmetry axes is found at the level of ∼ 5%. For the analysis, we have implemented an alternative inpainting approach in order to get an accurate reconstruction of the E-modes. More sensitive all-sky CMB polarization data, such as those expected from the future LiteBIRD experiment, are needed to reach a more robust conclusion on the possible existence of deviations from statistical isotropy in the form of a hemispherical power asymmetry.

Spectro-spatial evolution of the CMB. Part I. Discretisation of the thermalisation Green's function

Spectral distortions of the cosmic microwave background (CMB) have been recognized as an important future probe of the early Universe. Existing theoretical studies primarily focused on describing the evolution and creation of average distortions, ignoring spatial perturbations in the plasma. One of the main reasons for this choice is that a treatment of the spectro-spatial evolution of the photon field deep into the primordial Universe requires solving a radiative transfer problem for the distortion signals, which in full detail is computationally challenging. Here we provide the first crucial step towards tackling this problem by formulating a new spectral discretisation of the underlying average thermalisation Green's function. Our approach allows us to convert the high-dimensional partial differential equation system (≃ 103–104 equations) into and set of ordinary differential equations of much lower dimension (≃ 10 equations). We demonstrate the precision of the approach and highlight how it may be further improved in the future. We also clarify the link of the observable spectral distortion parameters (e.g., μ and y) to the computational spectral basis that we use in our frequency discretisation. This reveals how several basis-dependent ambiguities can be interpreted in future CMB analysis. Even if not exact, the new Green's function discretisation can be used to formulate a generalised photon Boltzmann-hierarchy, which can then be solved with methods that are familiar from theoretical studies of the CMB temperature and polarisation anisotropies. We will carry this program out in a series of companion papers, thereby opening the path to full spectro-spatial exploration of the CMB with future CMB imagers and spectrometers.

Photonic Integrated Circuits for Microwave Astronomy

One of the main objectives of cosmology is the measurement of the Cosmic Microwave Background (CMB) polarization and, in particular, the so-called B-modes, which could demonstrate the existence of the primordial gravitational waves generated at the early stages of the Universe. For this reason, a new integrated version of the heterodyne polarimeter photonic part shown on a previously proposed scheme that can be used in both direct imaging and interferometric instruments is presented. We have tested its properties using specific commercial software, obtaining promising results. Working as a direct imaging instrument, it is shown that the proposed polarimeter can provide sufficiently low polarization angle and polarization efficiency systematic errors, which are considered acceptable for current ground-based CMB polarization experiments dedicated to the characterization of the foreground signals affecting the lowest part of the frequency spectrum.

Disentangling patchy reionization signatures from primordial gravitational waves using CMB E-mode and B-mode polarization

ABSTRACT The detection of large angular scale B mode in the cosmic microwave background (CMB) polarization signal will open a direct window into not only the primary CMB anisotropies caused by the primordial gravitational waves (PGW) originating in the epoch of inflation, but also the secondary anisotropies imprinted during the epoch of cosmic reionization. The existence of patchiness in the electron density during reionization produces a unique distortion in the CMB B-mode polarization, which can be distinguished from the PGW signal with the aid of spatial frequency modes. In this work, we employ an EB estimator by combining E-mode and B-mode polarization for the τ power spectrum signal generated in a photon-conserving seminumerical reionization model called SCRIPT. We developed a Bayesian framework for the joint detection of the PGW and reionization signal from CMB observations and show the efficacy of this technique for upcoming CMB experiments. We find that, for our model, the τ power spectrum signal effectively tracks the inhomogeneous electron density field, allowing for robust constraints on the patchy B-mode signal. Further, our results indicate that employing the EB estimator for the τ signal will facilitate ground-based CMB-S4 to detect the patchy B-mode signal at approximately ≥2σ confidence level, while observations with space-based PICO will improve this detection to ≥3σ going as high as ≥7σ for extreme reionization models. These findings not only highlight the future potential of these experiments to provide an improved picture of the reionization process but also have important implications towards an unbiased measurement of r.

Minkowski functionals of CMB polarization intensity with pynkowski: theory and application to Planck and future data

ABSTRACT The angular power spectrum of the cosmic microwave background (CMB) anisotropies is a key tool to study the Universe. However, it is blind to the presence of non-Gaussianities and deviations from statistical isotropy, which can be detected with other statistics such as Minkowski functionals (MFs). These tools have been applied to CMB temperature and E-mode anisotropies with no detection of deviations from Gaussianity and isotropy. In this work, we extend the MF formalism to the CMB polarization intensity, P2 = Q2 + U2. We use the Gaussian kinematic formula to derive the theoretical predictions of MFs for Gaussian isotropic fields. We develop a software that computes MFs on P2healpix maps and apply it to simulations to verify the robustness of both theory and methodology. We then estimate MFs of P2 maps from Planck, both in pixel space and needlet domain, comparing them with realistic simulations that include CMB and instrumental noise residuals. We find no significant deviations from Gaussianity or isotropy in Planck CMB polarization intensity. However, MFs could play an important role in the analysis of CMB polarization measurements from upcoming experiments with improved sensitivity. Therefore, we forecast the ability of MFs applied to P2 maps to detect much fainter non-Gaussian anisotropic signals than with Planck data for two future complementary experiments: the LiteBIRD satellite and the ground-based Simons Observatory. We publicly release the software to compute MFs in arbitrary scalar healpix maps as a fully documented python package called pynkowski.