To deepen our understanding of Quantum Gravity and its connections with black holes and cosmology, building a common language and exchanging ideas across different approaches is crucial. The Nordita Program "Quantum Gravity: from gravitational effective field theories to ultraviolet complete approaches" created a platform for extensive discussions, aimed at pinpointing both common grounds and sources of disagreements, with the hope of generating ideas and driving progress in the field. This contribution summarizes the twelve topical discussions held during the program and collects individual thoughts of speakers and panelists on the future of the field in light of these discussions.

Cosmic inflation provides a window to the highest energy densities accessible in nature, far beyond those achievable in any realistic terrestrial experiment. Theoretical insights into the inflationary era and its observational probes may therefore shed unique light on the physical laws underlying our universe. This white paper describes our current theoretical understanding of the inflationary era, with a focus on the statistical properties of primordial fluctuations. In particular, we survey observational targets for three important signatures of inflation: primordial gravitational waves, primordial non-Gaussianity and primordial features. With the requisite advancements in analysis techniques, the tremendous increase in the raw sensitivities of upcoming and planned surveys will translate to leaps in our understanding of the inflationary paradigm and could open new frontiers for cosmology and particle physics. The combination of future theoretical and observational developments therefore offer the potential for a dramatic discovery about the nature of cosmic acceleration in the very early universe and physics on the smallest scales.

Axion as one of the promising dark matter candidates can be detected through narrow radio lines emitted from the magnetic white dwarf stars. Due to the existence of the strong magnetic field, the axion may resonantly convert into the radio photon (Primakoff effect) when it passes through a narrow region in the corona of the magnetic white dwarf, where the plasma frequency is equal to the axion mass. We show that for the magnetic white dwarf WD 2010+310, the future experiment SKA phase 1 with 100 hours of observation can effectively probe the parameter space of the axion-photon coupling $g_{a\gamma}$ up to $\sim 10^{-12}~ \text{GeV}^{-1}$ for the axion mass range of $0.2 \sim 3.7~ \mu$eV. Note that in the low mass region ($m_a \lesssim 1.5 ~\mu\text{eV}$), the WD 2010+310 could give greater sensitivity than the neutron star RX J0806.4-4123.

We present a measurement of the $B$-mode polarization power spectrum of the cosmic microwave background (CMB) using taken from July 2014 to December 2016 with the POLARBEAR experiment. The CMB power spectra are measured using observations at 150 GHz with an instantaneous array sensitivity of $\mathrm{NET}_\mathrm{array}=23\, \mu \mathrm{K} \sqrt{\mathrm{s}}$ on a 670 square degree patch of sky centered at (RA, Dec)=($+0^\mathrm{h}12^\mathrm{m}0^\mathrm{s},-59^\circ18^\prime$). A continuously rotating half-wave plate is used to modulate polarization and to suppress low-frequency noise. We achieve $32\,\mu\mathrm{K}$-$\mathrm{arcmin}$ effective polarization map noise with a knee in sensitivity of $\ell = 90$, where the inflationary gravitational wave signal is expected to peak. The measured $B$-mode power spectrum is consistent with a $\Lambda$CDM lensing and single dust component foreground model over a range of multipoles $50 \leq \ell \leq 600$. The data disfavor zero$C_\ell^{BB}$at$2.2\sigma$ using this $\ell$ range of POLARBEAR data alone. We cross-correlate our data with Planck high frequency maps and find the low-$\ell$ $B$-mode power in the combined dataset to be consistent with thermal dust emission. We place an upper limit on the tensor-to-scalar ratio r < 0.90 at 95% confidence level after marginalizing over foregrounds.

Very light pseudoscalar fields, often referred to as axions, are compelling dark matter candidates and can potentially be detected through their coupling to the electromagnetic field. Recently a novel detection technique using the cosmic microwave background (CMB) was proposed, which relies on the fact that the axion field oscillates at a frequency equal to its mass in appropriate units, leading to a time-dependent birefringence. For appropriate oscillation periods this allows the axion field at the telescope to be detected via the induced sinusoidal oscillation of the CMB linear polarization. We search for this effect in two years of POLARBEAR data. We do not detect a signal, and place a median $95 \%$ upper limit of $0.65 ^\circ$ on the sinusoid amplitude for oscillation frequencies between $0.02\,\text{days}^{-1}$ and $0.45\,\text{days}^{-1}$, which corresponds to axion masses between $9.6 \times 10^{-22} \, \text{eV}$and$2.2\times 10^{-20} \,\text{eV}$. Under the assumptions that 1) the axion constitutes all the dark matter and 2) the axion field amplitude is a Rayleigh-distributed stochastic variable, this translates to a limit on the axion-photon coupling $g_{\phi \gamma} < 2.4 \times 10^{-11} \,\text{GeV}^{-1} \times ({m_\phi}/{10^{-21} \, \text{eV}})$.

The upcoming Simons Observatory Small Aperture Telescopes aim at achieving a constraint on the primordial tensor-to-scalar ratio $r$ at the level of $\sigma(r=0)\lesssim0.003$, observing the polarized CMB in the presence of partial sky coverage, cosmic variance, inhomogeneous non-white noise, and Galactic foregrounds. We present three different analysis pipelines able to constrain $r$ given the latest available instrument performance, and compare their predictions on a set of sky simulations that allow us to explore a number of Galactic foreground models and elements of instrumental noise, relevant for the Simons Observatory. The three pipelines employ different combinations of parametric and non-parametric component separation at the map and power spectrum levels, and use B-mode purification to estimate the CMB B-mode power spectrum. We applied them to a common set of simulated realistic frequency maps, and compared and validated them with focus on their ability to extract robust constraints on the tensor-to-scalar ratio $r$. We evaluated their performance in terms of bias and statistical uncertainty on this parameter. In most of the scenarios the three methodologies achieve similar performance. Nevertheless, several simulations with complex foreground signals lead to a $>2\sigma$ bias on $r$ if analyzed with the default versions of these pipelines, highlighting the need for more sophisticated pipeline components that marginalize over foreground residuals. We show two such extensions, using power-spectrum-based and map-based methods, that are able to fully reduce the bias on $r$ below the statistical uncertainties in all foreground models explored, at a moderate cost in terms of $\sigma(r)$.

The detection of primordial $B$-mode polarisation of the Cosmic Microwave Background (CMB) is a major observational goal in modern Cosmology, offering a potential window into inflationary physics through the measurement of the tensor-to-scalar ratio $r$. However, the presence of Galactic foregrounds poses significant challenges, possibly biasing the $r$ estimate. In this study we explore the viability of using Minkowski functionals (MFs) as a robustness test to validate a potential $r$ detection by identifying non-Gaussian features associated with foregrounds contamination. To do so, we simulate sky maps as observed by a LiteBIRD-like CMB experiment, with realistic instrumental and foregrounds modelling. The CMB $B$-mode signal is recovered through blind component separation algorithms, and the obtained (biased) value of $r$ is used to generate Gaussian realisation of CMB signal. Their MFs are then compared with those computed on maps contaminated by foreground residual left by component separation, looking for a detection of non-Gaussianity. Our results demonstrate that, with the experimental configuration considered here, MFs can not be reliably adopted as a robustness test of an eventual $r$ detection, as we find that in the majority of the cases MFs are not able to raise significant warnings about the non-Gaussianity induced by the presence of foreground residuals. In the most realistic and refined scenario we adopted, the test is able to flag non-Gaussianity in $\sim 26\%$ of the simulations, meaning that there is no warning on the biased tensor-to-scalar ratio in $\sim 74\%$ of cases. These results suggest that more advanced statistics than MFs must be considered to look for non-Gaussian signatures of foregrounds, in order to be able to perform reliable null tests in future CMB missions.

Thermal dust is the major polarized foreground hindering the detection of primordial CMB B-modes. Its signal presents a complex behavior in frequency space inherited from the joint variation in our Galaxy of the orientation of magnetic fields and of the spectral properties of dust grains aligned with magnetic field lines. In this work, we present a new framework for analyzing the thermal dust signal using polarized microwave data. We introduce residual maps, represented as complex quantities, that capture deviations of the local polarized SED from the mean complex SED averaged over the sky mask. We make simple predictions that relate the values of the statistical correlation and covariances between the residual maps to the physical properties of the emitting aligned grains. The validation or falsification of these predictions provide valuable information on the nature of the dust signal. We confront our predictions with Planck data over a 97% mask excluding the inner Galactic plane. Despite its simplicity, our model captures a significant part of the statistical properties of the data. The spectral dependence of the covariances between residual maps is, for the Sroll2 version of the data, compatible with a dust model assuming only temperature variations rather than spectral index variations, while it is incompatible with either of them for the PR4 Planck official release. Our methodology can be used to analyze future high precision polarization data and build more accurate dust models to be used by the CMB community.

The Crab Nebula, also known as Tau A, is a polarized astronomical source at millimeter wavelengths. It has been used as a stable light source for polarization angle calibration in millimeter-wave astronomy. However, it is known that its intensity and polarization vary as a function of time at a variety of wavelengths. Thus, it is of interest to verify the stability of the millimeter-wave polarization. If detected, polarization variability may be used to better understand the dynamics of Tau~A, and for understanding the validity of Tau~A as a calibrator. One intriguing application of such observation is to use it for the search of axion-like particles (ALPs). Ultralight ALPs couple to photons through a Chern-Simons term, and induce a temporal oscillation in the polarization angle of linearly polarized sources. After assessing a number of systematic errors and testing for internal consistency, we evaluate the variability of the polarization angle of the Crab Nebula using 2015 and 2016 observations with the 150 GHz POLARBEAR instrument. We place a median 95% upper bound of polarization oscillation amplitude A &lt; 0.065^\circ over the oscillation frequencies from $0.75~\mathrm{year}^{-1}$ to $0.66~\mathrm{hour}^{-1}$. Assuming that no sources other than ALP are causing Tau A's polarization angle variation, that the ALP constitutes all the dark matter, and that the ALP field is a stochastic Gaussian field, this bound translates into a median 95% upper bound of ALP-photon coupling $g_{a\gamma\gamma} < 2.16\times10^{-12}\,\mathrm{GeV}^{-1}\times(m_a/10^{-21} \mathrm{eV})$in the mass range from$9.9\times10^{-23} \mathrm{eV}$ to $7.7\times10^{-19} \mathrm{eV}$. This demonstrates that this type of analysis using bright polarized sources is as competitive as those using the polarization of cosmic microwave background in constraining ALPs.

We report a measurement of the E-mode polarization power spectrum of the cosmic microwave background (CMB) using 150 GHz data taken from July 2014 to December 2016 with the POLARBEAR experiment. We reach an effective polarization map noise level of $32\,\mu\mathrm{K}$-$\mathrm{arcmin}$ across an observation area of 670 square degrees. We measure the EE power spectrum over the angular multipole range 500 \leq \ell &lt;3000, tracing the third to seventh acoustic peaks with high sensitivity. The statistical uncertainty on E-mode bandpowers is $\sim 2.3 \mu {\rm K}^2$ at $\ell \sim 1000$ with a systematic uncertainty of 0.5$\mu {\rm K}^2$. The data are consistent with the standard $\Lambda$CDM cosmological model with a probability-to-exceed of 0.38. We combine recent CMB E-mode measurements and make inferences about cosmological parameters in $\Lambda$CDM as well as in extensions to $\Lambda$CDM. Adding the ground-based CMB polarization measurements to the Planck dataset reduces the uncertainty on the Hubble constant by a factor of 1.2 to $H_0 = 67.20 \pm 0.57 {\rm km\,s^{-1} \,Mpc^{-1}}$. When allowing the number of relativistic species ($N_{eff}$) to vary, we find $N_{eff} = 2.94 \pm 0.16$, which is in good agreement with the standard value of 3.046. Instead allowing the primordial helium abundance ($Y_{He}$) to vary, the data favor $Y_{He} = 0.248 \pm 0.012$. This is very close to the expectation of 0.2467 from Big Bang Nucleosynthesis. When varying both $Y_{He}$ and $N_{eff}$, we find $N_{eff} = 2.70 \pm 0.26$ and $Y_{He} = 0.262 \pm 0.015$.

The upcoming Simons Observatory (SO) Small Aperture Telescopes aim at observing the degree-scale anisotropies of the polarized CMB to constrain the primordial tensor-to-scalar ratio $r$ at the level of $\sigma(r=0)\lesssim0.003$ to probe models of the very early Universe. We present three complementary $r$ inference pipelines and compare their results on a set of sky simulations that allow us to explore a number of Galactic foreground and instrumental noise models, relevant for SO. In most scenarios, the pipelines retrieve consistent and unbiased results. However, several complex foreground scenarios lead to a $>2\sigma$ bias on $r$ if analyzed with the default versions of these pipelines, highlighting the need for more sophisticated pipeline components that marginalize over foreground residuals. We present two such extensions, using power-spectrum-based and map-based methods, and show that they fully reduce the bias on $r$ to sub-sigma level in all scenarios, and at a moderate cost in terms of $\sigma(r)$.

Polarized foreground emission from the Galaxy is one of the biggest challenges facing current and upcoming cosmic microwave background (CMB) polarization experiments. We develop new models of polarized Galactic dust and synchrotron emission at CMB frequencies that draw on the latest observational constraints, that employ the ``polarization fraction tensor'' framework to couple intensity and polarization in a physically motivated way, and that allow for stochastic realizations of small-scale structure at sub-arcminute angular scales currently unconstrained by full-sky data. We implement these models into the publicly available Python Sky Model (PySM) software and additionally provide PySM interfaces to select models of dust and CO emission from the literature. We characterize the behavior of each model by quantitatively comparing it to observational constraints in both maps and power spectra, demonstrating an overall improvement over previous PySM models. Finally, we synthesize models of the various Galactic foreground components into a coherent suite of three plausible microwave skies that span a range of astrophysical complexity allowed by current data.

The detection of primordial $B$-mode polarisation of the Cosmic Microwave Background (CMB) is a major observational goal in modern Cosmology, offering a potential window into inflationary physics through the measurement of the tensor-to-scalar ratio $r$. However, the presence of Galactic foregrounds poses significant challenges, possibly biasing the $r$ estimate. In this study we explore the viability of using Minkowski functionals (MFs) as a robustness test to validate a potential $r$ detection by identifying non-Gaussian features associated with foregrounds contamination. To do so, we simulate sky maps as observed by a LiteBIRD-like CMB experiment, with realistic instrumental and foregrounds modelling. The CMB $B$-mode signal is recovered through blind component separation algorithms, and the obtained (biased) value of $r$ is used to generate Gaussian realisation of CMB signal. Their MFs are then compared with those computed on maps contaminated by foreground residual left by component separation, looking for a detection of non-Gaussianity. Our results demonstrate that, with the experimental configuration considered here, MFs can not be reliably adopted as a robustness test of an eventual $r$ detection, as we find that in the majority of the cases MFs are not able to raise significant warnings about the non-Gaussianity induced by the presence of foreground residuals. In the most realistic and refined scenario we adopted, the test is able to flag non-Gaussianity in $\sim 26\%$ of the simulations, meaning that there is no warning on the biased tensor-to-scalar ratio in $\sim 74\%$ of cases. These results suggest that more advanced statistics than MFs must be considered to look for non-Gaussian signatures of foregrounds, in order to be able to perform reliable null tests in future CMB missions.

We quantify the calibration requirements for systematic uncertainties for next-generation ground-based observatories targeting the large-angle $B$-mode polarization of the Cosmic Microwave Background, with a focus on the Simons Observatory (SO). We explore uncertainties on gain calibration, bandpass center frequencies, and polarization angles, including the frequency variation of the latter across the bandpass. We find that gain calibration and bandpass center frequencies must be known to percent levels or less to avoid biases on the tensor-to-scalar ratio $r$ on the order of $\Delta r\sim10^{-3}$, in line with previous findings. Polarization angles must be calibrated to the level of a few tenths of a degree, while their frequency variation between the edges of the band must be known to ${\cal O}(10)$ degrees. Given the tightness of these calibration requirements, we explore the level to which residual uncertainties on these systematics would affect the final constraints on $r$ if included in the data model and marginalized over. We find that the additional parameter freedom does not degrade the final constraints on $r$ significantly, broadening the error bar by ${\cal O}(10\%)$ at most. We validate these results by reanalyzing the latest publicly available data from the BICEP2/Keck collaboration within an extended parameter space covering both cosmological, foreground and systematic parameters. Finally, our results are discussed in light of the instrument design and calibration studies carried out within SO.

We report a 4.8$\sigma$ measurement of the cross-correlation signal between the cosmic microwave background (CMB) lensing convergence reconstructed from measurements of the CMB polarization made by the POLARBEAR experiment and the infrared-selected galaxies of the Herschel-ATLAS survey. This is the first measurement of its kind. We infer a best-fit galaxy bias of $b = 5.76 \pm 1.25$, corresponding to a host halo mass of$\log_{10}(M_h/M_\odot) =13.5^{+0.2}_{-0.3}$at an effective redshift of$z \sim 2$ from the cross-correlation power spectrum. Residual uncertainties in the redshift distribution of the sub-mm galaxies are subdominant with respect to the statistical precision. We perform a suite of systematic tests, finding that instrumental and astrophysical contaminations are small compared to the statistical error. This cross-correlation measurement only relies on CMB polarization information that, differently from CMB temperature maps, is less contaminated by galactic and extra-galactic foregrounds, providing a clearer view of the projected matter distribution. This result demonstrates the feasibility and robustness of this approach for future high-sensitivity CMB polarization experiments.

Cosmic Microwave Background (CMB) observations are used to constrain reheating to Standard Model (SM) particles after a period of inflation. As a light spectator field, the SM Higgs boson acquires large field values from its quantum fluctuations during inflation, gives masses to SM particles that vary from one Hubble patch to another, and thereby produces large density fluctuations. We consider both perturbative and resonant decay of the inflaton to SM particles. For the case of perturbative decay from coherent oscillations of the inflaton after high scale inflation, we find strong constraints on the reheat temperature for the inflaton decay into heavy SM particles. For the case of resonant particle production (preheating) to (Higgsed) SM gauge bosons, we find temperature fluctuations larger than observed in the CMB for a range of gauge coupling that includes those found in the SM and conclude that such preheating cannot be the main source of reheating the Universe after inflation.

The optical depth to reionization, $\tau$, is the least constrained parameter of the cosmological $\Lambda$CDM model. To date, its most precise value is inferred from large-scale polarized CMB power spectra from the ${\it Planck}$ High-Frequency Instrument (HFI). These maps are known to contain significant contamination by residual non-Gaussian systematic effects, which are hard to model analytically. Therefore, robust constraints on $\tau$ are currently obtained through an empirical cross-spectrum likelihood built from simulations. In this paper, we present a likelihood-free inference of $\tau$ from polarized ${\it Planck}$ HFI maps which, for the first time, is fully based on neural networks (NNs). NNs have the advantage of not requiring an analytical description of the data and can be trained on state-of-the-art simulations, combining information from multiple channels. By using Gaussian sky simulations and ${\it Planck}$ ${\tt SRoll2}$ simulations, including CMB, noise, and residual instrumental systematic effects, we train, test and validate NN models considering different setups. We infer the value of $\tau$ directly from $Q$ and $U$ maps at $\sim 4^\circ$ pixel resolution, without computing power spectra. On ${\it Planck}$ data, we obtain $\tau_{NN}=0.058 \pm 0.008$, compatible with current EE cross-spectrum results but with a $\sim30\%$ larger uncertainty, which can be assigned to the inherent non-optimality of our estimator and to the retraining procedure applied to avoid biases. While this paper does not improve on current cosmological constraints, our analysis represents a first robust application of NN-based inference on real data and highlights its potential as a promising tool for complementary analysis of near-future CMB experiments, also in view of the ongoing challenge to achieve a detection of primordial gravitational waves.

Paleo-detectors are a proposed experimental technique to search for dark matter (DM). In lieu of the conventional approach of operating a tonne-scale real-time detector to search for DM-induced nuclear recoils, paleo-detectors take advantage of small samples of naturally occurring rocks on Earth that have been deep underground ($\gtrsim 5$ km), accumulating nuclear damage tracks from recoiling nuclei for $\mathcal{O}(1)$ Gyr. Modern microscopy techniques promise the capability to read out nuclear damage tracks with nanometer resolution in macroscopic samples. Thanks to their $\mathcal{O}(1)$ Gyr integration times, paleo-detectors could constitute nuclear recoil detectors with keV recoil energy thresholds and 100 kilotonne-yr exposures. This combination would allow paleo-detectors to probe DM-nucleon cross sections orders of magnitude below existing upper limits from conventional direct detection experiments. In this article, we use improved background modeling and a new spectral analysis technique to update the sensitivity forecast for paleo-detectors. We demonstrate the robustness of the sensitivity forecast to the (lack of) ancillary measurements of the age of the samples and the parameters controlling the backgrounds, systematic mismodeling of the spectral shape of the backgrounds, and the radiopurity of the mineral samples. Specifically, we demonstrate that even if the uranium concentration in paleo-detector samples is $10^{-8}$ (per weight), many orders of magnitude larger than what we expect in the most radiopure samples obtained from ultra basic rock or marine evaporite deposits, paleo-detectors could still probe DM-nucleon cross sections below current limits. For DM masses $\lesssim 10$ GeV/$c^2$, the sensitivity of paleo-detectors could still reach down all the way to the conventional neutrino floor in a Xe-based direct detection experiment.

Phonons in Bose-Einstein condensates propagate as massless scalar particles on top of an emergent acoustic metric. This hydrodynamics/gravity analogy can be exploited to realize acoustic black holes, featuring an event horizon that traps phonons. We show that by an appropriate external potential, gravitational wave-like perturbations of the acoustic metric can be produced. Such perturbations can be used to excite an acoustic black hole, which should then relax by phonon emission.