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.

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.

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)$.

A significant fraction of compact-object mergers in galactic nuclei are expected to be eccentric in the Laser-Interferometer-Space-Antenna (LISA) frequency sensitivity range, $10^{-4} - 10^{-1}\ \rm Hz$. Several compact binaries detected by the LIGO-Virgo-KAGRA Collaboration may retain hints of residual eccentricity at $10$~Hz, suggesting dynamical or triple origins for a significant fraction of the gravitational-wave-observable population. In triple systems, von-Zeipel-Lidov-Kozai oscillations perturb both the eccentricity and the argument of pericentre, $\omega$, of the inner black hole binary. The latter could be fully \textit{circulating}, where $\omega$ cycles through $2\pi$, or may \textit{librate}, with $\omega$ ranges about a fixed value with small or large variation. We use \texttt{TSUNAMI}, a regularised N-body code with up to 3.5 post-Newtonian (PN) term corrections, to identify four different families of orbits: (i) circulating, (ii) small and (iii) large amplitude librating, and (iv) merging orbits. We develop and demonstrate a new method to construct gravitational waveforms using the quadrupole formula utilising the instantaneous {\it total} acceleration of each binary component in \texttt{TSUNAMI}. We show that the four orbital families have distinct waveform phenomenologies, enabling them to be distinguished if observed in LISA. The orbits are also distinguishable from an isolated binary or from a binary perturbed by a different tertiary orbit, even if the secular timescale is the same. Future burst timing models will be able to distinguish the different orbital configurations. For efficient binary formation, about $\sim 50$ binaries can have librating orbits in the Galactic Centre.

One of the challenges of high granularity calorimeters, such as that to be built to cover the endcap region in the CMS Phase-2 Upgrade for HL-LHC, is that the large number of channels causes a surge in the computing load when clustering numerous digitised energy deposits (hits) in the reconstruction stage. In this article, we propose a fast and fully-parallelizable density-based clustering algorithm, optimized for high occupancy scenarios, where the number of clusters is much larger than the average number of hits in a cluster. The algorithm uses a grid spatial index for fast querying of neighbours and its timing scales linearly with the number of hits within the range considered. We also show a comparison of the performance on CPU and GPU implementations, demonstrating the power of algorithmic parallelization in the coming era of heterogeneous computing in high energy physics.

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 investigate how AGN disk turbulence affects the orbital dynamics of a stellar-mass black hole (BH) initially located at a migration trap, focusing on the long-term behavior of eccentricity and inclination in the quasi-embedded regime. We develop a semi-analytical framework in which turbulence is modeled as a stochastic velocity field acting through a modified drag force. We integrate the resulting stochastic differential equations both in Cartesian coordinates and in orbital elements using a linearized perturbative approach, and compare these results with full numerical simulations. Eccentricity and inclination evolve toward steady-state Rayleigh distributions, with variances determined by the local disk properties and the ratio of the gas damping rate to the orbital frequency. The analytical predictions agree well with the numerical simulations. We provide closed-form expressions for the variances in both the fast and slow damping regimes. These results are directly applicable to Monte Carlo population models and can serve as physically motivated initial conditions for hydrodynamical simulations. Turbulent forcing prevents full circularization and alignment of BH orbits in AGN disks, even in the presence of strong gas drag. This has important implications for BH merger and binary formation rates, which are sensitive to the residual eccentricity and inclination. Our results highlight the need to account for turbulence-induced stochastic heating when modeling the dynamical evolution of compact objects in AGN environments.

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)$.

We have applied an operator-overloading forward-mode algorithmic differentiation tool to the Monte-Carlo particle simulation toolkit Geant4. Our differentiated version of Geant4 allows computing mean pathwise derivatives of user-defined outputs of Geant4 applications with respect to user-defined inputs. This constitutes a major step towards enabling gradient-based optimization techniques in high-energy physics, as well as other application domains of Geant4. This is a preliminary report on the technical aspects of applying operator-overloading AD to Geant4, as well as a first analysis of some results obtained by our differentiated Geant4 prototype. We plan to follow up with a more refined analysis.

In the weak-coupling limit approach to open quantum systems, the presence of the bath is eliminated and accounted for by a master equation that introduces dissipative contributions to the system reduced dynamics. Within this framework, there are no bath entropy contributions to the entropy balance. We show that, as a consequence, the entropy production fails to be positive for a class of physically legitimate (i.e., completely positive and trace preserving) non-Markovian dynamical maps. Moreover, in the absence of the semigroup property, if the reduced dynamics has a thermal asymptotic state, this need not be stationary. In this case, even the integrated entropy production becomes negative. These observations imply that, when the conditions leading to reduced dynamics of semigroup type are relaxed, a consistent formulation of the second law of thermodynamics requires that the environment contribution to the entropy balance be explicitly taken into account.

We present the first measurement of cross-correlation between the lensing potential, reconstructed from cosmic microwave background (CMB) {\it polarization} data, and the cosmic shear field from galaxy shapes. This measurement is made using data from the POLARBEAR CMB experiment and the Subaru Hyper Suprime-Cam (HSC) survey. By analyzing an 11~deg$^2$ overlapping region, we reject the null hypothesis at 3.5$\sigma$\ and constrain the amplitude of the {\bf cross-spectrum} to $\widehat{A}_{\rm lens}=1.70\pm 0.48$, where $\widehat{A}_{\rm lens}$ is the amplitude normalized with respect to the Planck~2018{} prediction, based on the flat $\Lambda$ cold dark matter cosmology. The first measurement of this {\bf cross-spectrum} without relying on CMB temperature measurements is possible due to the deep POLARBEAR map with a noise level of ${\sim}$6\,$\mu$K-arcmin, as well as the deep HSC data with a high galaxy number density of $n_g=23\,{\rm arcmin^{-2}}$. We present a detailed study of the systematics budget to show that residual systematics in our results are negligibly small, which demonstrates the future potential of this cross-correlation technique.

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.

Two-dimensional materials on metallic surfaces or stacked one on top of the other can form a variety of moir\'e superstructures depending on the possible parameter and symmetry mismatch and misorientation angle. In most cases, such as incommensurate lattices or identical lattices but with a small twist angle, the common periodicity may be very large, thus making numerical simulations prohibitive. We propose here a general procedure to determine the minimal simulation cell which approximates, within a certain tolerance and a certain size, the primitive cell of the common superlattice, given the two interfacing lattices and the relative orientation angle. As case studies to validate our procedure, we report two applications of particular interest: the case of misaligned hexagonal/hexagonal identical lattices, describing a twisted graphene bilayer or a graphene monolayer grown on Ni(111), and the case of hexagonal/square lattices, describing for instance a graphene monolayer grown on Ni(100) surface. The first one, which has also analytic solutions, constitutes a solid benchmark for the algorithm; the second one shows that a very nice description of the experimental observations can be obtained also using the resulting relatively small coincidence cells.

The analysis of photometric large-scale structure data is often complicated by the need to account for many observational and astrophysical systematics. The elaborate models needed to describe them often introduce many ``nuisance parameters'', which can be a major inhibitor of an efficient parameter inference. In this paper we introduce an approximate method to analytically marginalise over a large number of nuisance parameters based on the Laplace approximation. We discuss the mathematics of the method, its relation to concepts such as volume effects and profile likelihood, and show that it can be further simplified for calibratable systematics by linearising the dependence of the theory on the associated parameters. We quantify the accuracy of this approach by comparing it with traditional sampling methods in the context of existing data from the Dark Energy Survey, as well as futuristic Stage-IV photometric data. The linearised version of the method is able to obtain parameter constraints that are virtually equivalent to those found by exploring the full parameter space for a large number of calibratable nuisance parameters, while reducing the computation time by a factor 3-10. Furthermore, the non-linearised approach is able to analytically marginalise over a large number of parameters, returning constraints that are virtually indistinguishable from the brute-force method in most cases, accurately reproducing both the marginalised uncertainty on cosmological parameters, and the impact of volume effects associated with this marginalisation. We provide simple recipes to diagnose when the approximations made by the method fail and one should thus resort to traditional methods. The gains in sampling efficiency associated with this method enable the joint analysis of multiple surveys, typically hindered by the large number of nuisance parameters needed to describe them.

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.

The High Energy Modular Array (HEMA) is one of three instruments that compose the STROBE-X mission concept. The HEMA is a large-area, high-throughput non-imaging pointed instrument based on the Large Area Detector developed as part of the LOFT mission concept. It is designed for spectral timing measurements of a broad range of sources and provides a transformative increase in sensitivity to X-rays in the energy range of 2--30 keV compared to previous instruments, with an effective area of 3.4 m$^{2}$ at 8.5 keV and an energy resolution of better than 300 eV at 6 keV in its nominal field of regard.

From a study of the subleading structure of the asymptotic equations of motion in Einstein-Yang-Mills theory, we construct charges that are conserved up to quadratic order in non-radiative vacuum. We then show that these higher spin charges obey the celestial $sw_{1+\infty}$ symmetry algebra found earlier from the OPE of positive-helicity conformally soft gluons and gravitons.

We investigate the time evolution of a non-resonant dressed-atom qubit in an XZ original configuration. It is composed of two electromagnetic fields, one oscillating parallel and the other orthogonal to the quantisation magnetic static field. The experiments are performed in rubidium and caesium atomic magnetometers, confined in a magneto-optical trap and in a vapour cell, respectively. Static fields in the $\mu$T range and kHz oscillating fields with large Rabi frequencies are applied. This dual-dressing configuration is an extension of the Landau-Zener multipassage interferometry in the presence of an additional dressing field controlling the tunneling process by its amplitude and phase. Our measurement of the qubit coherence introduces additional features to the transition probability readout of standard interferometry. The coherence time evolution is characterized by oscillations at several frequencies, each of them produced by a different quantum contribution. Such frequency description introduces a new picture of the qubit multipassage evolution. Because the present low-frequency dressing operation does not fall within the standard Floquet engineering paradigm based on the high-frequency expansion, we develop an ad-hoc dressing perturbation treatment. Numerical simulations support the adiabatic and non-adiabatic qubit evolution.

Gamma rays are produced by cosmic ray (CR) protons interacting with the particles at solar photosphere and by cosmic ray electrons and positrons (CRes) via inverse Compton scattering of solar photons. The former come from the solar disk while the latter extend beyond the disk. Evaluation of these emissions requires the flux and spectrum of CRs in the vicinity of the Sun, while most observations provide flux and spectra near the Earth, at around 1 AU from the Sun. Past estimates of the quiet Sun gamma-ray emission use phenomenological modulation procedures to estimate spectra near the Sun (see review by Orlando and Strong 2021 and references therein). We show that CRe transport in the inner heliosphere requires a kinetic approach and use a novel approximation to determine the variation of CRe flux and spectrum from 1 AU to the Sun including effects of (1) the structure of large scale magnetic field, (2) small scale turbulence in the solar wind from several in situ measurements, in particular, those by Parker Solar Probe that extend this information to 0.1 AU, and (3) most importantly, energy losses due to synchrotron and inverse Compton processes. We present results on the flux and spectrum variation of CRes from 1 AU to the Sun for several transport models. In forthcoming papers we will use these results for a more accurate estimate of quiet Sun inverse Compton gamma-ray spectra, and, for the first time, the spectrum of extreme ultraviolet to hard X-ray photons produced by synchrotron emission. These can be compared with the quiet Sun gamma-ray observation by Fermi (see, e.g.~Fermi-LAT Collaboration, 2011) and X-ray upper limits set by RHESSI (Hannah et al., 2010).