Divergence in perturbative expansions is where interesting physics takes place. Particle production on time-dependent backgrounds, as one such example, is interpreted as transition from one vacuum to another. Vacuum is typically defined as an asymptotic state in which the WKB approximation is valid. The use of the WKB method, however, poses several conceptual and computational issues, as the WKB series is divergent in general, quantization is insensitive to higher orders in the series, and the global behavior of solutions cannot be captured. Exact WKB analysis is a powerful resummation technology that provides an analytical tool for a global structure of exact solutions to overcome these problems. In this paper, we establish quantization by fully employing the exact WKB solutions as mode functions and by defining the vacua with respect to them. We provide a self-contained exact WKB formulation to obtain evolution matrices without resorting to the use of known special functions and without approximations. We find that the quantity called Voros coefficient plays an important role to re-normalize the exact WKB solutions compatible with asymptotic states. We show that the ambiguity that coexists with nontrivial Voros coefficients is eliminated by requiring physical quantization conditions. Our formalism provides a conceptual as well as practical framework to upgrade our treatment of quantization and particle production. Combined with other approximating techniques, it can form a basis to tackle a broad class of problems that are beyond technical ability of the existing formulations.

Recent investigations of vacuum polarization in extremely compact stars suggest that, in such regimes, the effective matter content of spacetime may acquire a vacuum-energy equation of state with negative energy density, mimicking a negative cosmological constant. Motivated by this observation, we introduce a general algorithm to modify well-known spherically symmetric regular black hole metrics by replacing their usual de Sitter cores (dSC) with Anti-de Sitter cores (AdSC). Like their dSC counterparts, these AdSC solutions may exhibit two, one, or no horizons depending on the value of a regularization parameter l. We present explicit examples of AdSC-Bardeen and AdSC-Dymnikova metrics, analyze their main properties, and investigate some of their phenomenological signatures using test fields. In particular, we compare their fundamental quasinormal modes and echo signals with those of the dSC cases, highlighting potential avenues for distinguishing them observationally.

Cosmic microwave background (CMB) photons are deflected by large-scale structure through gravitational lensing. This secondary effect introduces higher-order correlations in CMB anisotropies, which are used to reconstruct lensing deflections. This allows mapping of the integrated matter distribution along the line of sight, probing the growth of structure, and recovering an undistorted view of the last-scattering surface. Gravitational lensing has been measured by previous CMB experiments, with $\textit{Planck}$'s $42\,\sigma$ detection being the current best full-sky lensing map. We present an enhanced $\textit{LiteBIRD}$ lensing map by extending the CMB multipole range and including the minimum-variance estimation, leading to a $49$ to $58\,\sigma$ detection over $80\,\%$ of the sky, depending on the final complexity of polarized Galactic emission. The combination of $\textit{Planck}$ and $\textit{LiteBIRD}$ will be the best full-sky lensing map in the 2030s, providing a $72$ to $78\,\sigma$ detection over $80\,\%$ of the sky, almost doubling $\textit{Planck}$'s sensitivity. Finally, we explore different applications of the lensing map, including cosmological parameter estimation using a lensing-only likelihood and internal delensing, showing that the combination of both experiments leads to improved constraints. The combination of $\textit{Planck}$ + $\textit{LiteBIRD}$ will improve the $S_8$ constraint by a factor of 2 compared to $\textit{Planck}$, and $\textit{Planck}$ + $\textit{LiteBIRD}$ internal delensing will improve $\textit{LiteBIRD}$'s tensor-to-scalar ratio constraint by $6\,\%$. We have tested the robustness of our results against foreground models of different complexity, showing that improvements remains even for the most complex foregrounds.

Composite axions offer a scenario where the axion emerges as a pion-like state, avoiding fine-tuning of elementary scalars and ameliorating the axion quality problem. Despite these advantages, their post-inflationary cosmology remains largely unexplored, with challenges including the domain wall problem and the presence of exotic relics. We propose two composite axion models with an effective domain wall number $N_\text{DW} = 1$ and study the dilution of relics via a short period of inflation. One model is based on an $SU(5)$ chiral gauge theory, while the other employs a ``gauged'' $U(1)$ Peccei-Quinn symmetry in vector-like $SU(N)$ gauge theories. We identify the viable parameter space in which axion strings re-enter the horizon before or even after the QCD transition and axion dark matter is dominantly produced from the decay of the string-wall network.

The Euclid mission aims to measure the positions, shapes, and redshifts of over a billion galaxies to provide unprecedented constraints on the nature of dark matter and dark energy. Achieving this goal requires a continuous reassessment of the mission's scientific performance, particularly in terms of its ability to constrain cosmological parameters, as our understanding of how to model large-scale structure observables improves. In this study, we present the first scientific forecasts using CLOE (Cosmology Likelihood for Observables in Euclid), a dedicated Euclid cosmological pipeline developed to support this endeavour. Using advanced Bayesian inference techniques applied to synthetic Euclid-like data, we sample the posterior distribution of cosmological and nuisance parameters across a variety of cosmological models and Euclid primary probes: cosmic shear, angular photometric galaxy clustering, galaxy-galaxy lensing, and spectroscopic galaxy clustering. We validate the capability of CLOE to produce reliable cosmological forecasts, showcasing Euclid's potential to achieve a figure of merit for the dark energy parameters $w_0$ and $w_a$ exceeding 400 when combining all primary probes. Furthermore, we illustrate the behaviour of the posterior probability distribution of the parameters of interest given different priors and scale cuts. Finally, we emphasise the importance of addressing computational challenges, proposing further exploration of innovative data science techniques to efficiently navigate the Euclid high-dimensional parameter space in upcoming cosmological data releases.

A methodology to provide the polarization angle requirements for different sets of detectors, at a given frequency of a CMB polarization experiment, is presented. The uncertainties in the polarization angle of each detector set are related to a given bias on the tensor-to-scalar ratio $r$ parameter. The approach is grounded in using a linear combination of the detector sets to obtain the CMB polarization signal. In addition, assuming that the uncertainties on the polarization angle are in the small angle limit (lower than a few degrees), it is possible to derive analytic expressions to establish the requirements. The methodology also accounts for possible correlations among detectors, that may originate from the optics, wafers, etc. The approach is applied to the LiteBIRD space mission. We show that, for the most restrictive case (i.e., full correlation of the polarization angle systematics among detector sets), the requirements on the polarization angle uncertainties are of around 1 arcmin at the most sensitive frequency bands (i.e., $\approx 150$ GHz) and of few tens of arcmin at the lowest (i.e., $\approx 40$ GHz) and highest (i.e., $\approx 400$ GHz) observational bands. Conversely, for the least restrictive case (i.e., no correlation of the polarization angle systematics among detector sets), the requirements are $\approx 5$ times less restrictive than for the previous scenario. At the global and the telescope levels, polarization angle knowledge of a few arcmins is sufficient for correlated global systematic errors and can be relaxed by a factor of two for fully uncorrelated errors in detector polarization angle. The reported uncertainty levels are needed in order to have the bias on $r$ due to systematics below the limit established by the LiteBIRD collaboration.

Various so-called anomalies have been found in both the WMAP and Planck cosmic microwave background (CMB) temperature data that exert a mild tension against the highly successful best-fit 6 parameter cosmological model, potentially providing hints of new physics to be explored. That these are real features on the sky is uncontested. However, given their modest significance, whether they are indicative of true departures from the standard cosmology or simply statistical excursions, due to a mildly unusual configuration of temperature anisotropies on the sky which we refer to as the "fluke hypothesis", cannot be addressed further without new information. No theoretical model of primordial perturbations has to date been constructed that can explain all of the temperature anomalies. Therefore, we focus in this paper on testing the fluke hypothesis, based on the partial correlation between the temperature and $E$-mode CMB polarisation signal. In particular, we compare the properties of specific statistics in polarisation, built from unconstrained realisations of the $\Lambda$CDM cosmological model as might be observed by the LiteBIRD satellite, with those determined from constrained simulations, where the part of the $E$-mode anisotropy correlated with temperature is constrained by observations of the latter. Specifically, we use inpainted Planck 2018 SMICA temperature data to constrain the $E$-mode realisations. Subsequent analysis makes use of masks defined to minimise the impact of the inpainting procedure on the $E$-mode map statistics. We find that statistical assessments of the $E$-mode data alone do not provide any evidence for or against the fluke hypothesis. However, tests based on cross-statistical measures determined from temperature and $E$ modes can allow this hypothesis to be rejected with a moderate level of probability.

The Python Sky Model (PySM) is a Python package used by Cosmic Microwave Background (CMB) experiments to simulate maps, in HEALPix pixelization, of the various diffuse astrophysical components of Galactic emission relevant at CMB frequencies (i.e. dust, synchrotron, free-free and Anomalous Microwave Emission), as well as the CMB itself. These maps may be integrated over a given instrument bandpass and smoothed with a given instrument beam. PySM 2, released in 2016, has become the de-facto standard for simulating Galactic emission, for example it is used by CMB-S4, Simons Observatory, LiteBird, PICO, CLASS, POLARBEAR and other CMB experiments, as shown by the 80+ citations of the PySM 2 publication. As the resolution of upcoming experiments increases, the PySM 2 software has started to show some limitations, the solution to these issues was to reimplement PySM from scratch focusing on these features: reimplement all the models with the numba Just-In-Time compiler for Python to reduce memory overhead and optimize performance; use MPI through mpi4py to coordinate execution of PySM 3 across multiple nodes and rely on libsharp, for distributed spherical harmonic transforms; employ the data utilities infrastructure provided by astropy to download the input templates and cache them when requested. At this stage we strive to maintain full compatibility with PySM 2, therefore we implement the exact same astrophysical emission models with the same naming scheme. In the extensive test suite we compare the output of each PySM 3 model with the results obtained by PySM 2.

We present ForSE+, a Python package that produces non-Gaussian diffuse Galactic thermal dust emission maps at arcminute angular scales and that has the capacity to generate random realizations of small scales. This represents an extension of the ForSE (Foreground Scale Extender) package, which was recently proposed to simulate non-Gaussian small scales of thermal dust emission using generative adversarial networks (GANs). With the input of the large-scale polarization maps from observations, ForSE+ has been trained to produce realistic polarized small scales at 3' following the statistical properties, mainly the non-Gaussianity, of observed intensity small scales, which are evaluated through Minkowski functionals. Furthermore, by adding different realizations of random components to the large-scale foregrounds, we show that ForSE+ is able to generate small scales in a stochastic way. In both cases, the output small scales have a similar level of non-Gaussianity compared with real observations and correct amplitude scaling as a power law. These realistic new maps will be useful, in the future, to understand the impact of non-Gaussian foregrounds on the measurements of the cosmic microwave background (CMB) signal, particularly on the lensing reconstruction, de-lensing, and the detection of cosmological gravitational waves in CMB polarization B-modes.

Boson stars have attracted much attention in recent decades as simple, self-consistent models of compact objects and also as self-gravitating structures formed in some dark-matter scenarios. Direct detection of these hypothetical objects through electromagnetic signatures would be unlikely because their bosonic constituents are not expected to interact significantly with ordinary matter and radiation. However, binary boson stars might form and coalesce emitting a detectable gravitational wave signal which might distinguish them from ordinary compact object binaries containing black holes and neutron stars. We study the merger of two boson stars by numerically evolving the fully relativistic Einstein-Klein-Gordon equations for a complex scalar field with a solitonic potential that generates very compact boson stars. Owing to the steep mass-radius diagram, we can study the dynamics and gravitational radiation from unequal-mass binary boson stars with mass ratios up to $q\approx23$ without the difficulties encountered when evolving binary black holes with large mass ratios. Similar to the previously-studied equal-mass case, our numerical evolutions of the merger produce either a nonspinning boson star or a spinning black hole, depending on the initial masses and on the binary angular momentum. We do not find any evidence of synchronized scalar clouds forming around either the remnant spinning black hole or around the remnant boson stars. Interestingly, in contrast to the equal-mass case, one of the mechanisms to dissipate angular momentum is now asymmetric, and leads to large kick velocities (up to a few $10^4\,{\rm km/s}$) which could produce wandering remnant boson stars. We also compare the gravitational wave signals predicted from boson star binaries with those from black hole binaries, and comment on the detectability of the differences with ground interferometers.

This white paper discusses the current landscape and prospects for experiments sensitive to particle dark matter processes producing photons and cosmic rays. Much of the gamma-ray sky remains unexplored on a level of sensitivity that would enable the discovery of a dark matter signal. Currently operating GeV-TeV observatories, such as Fermi-LAT, atmospheric Cherenkov telescopes, and water Cherenkov detector arrays continue to target several promising dark matter-rich environments within and beyond the Galaxy. Soon, several new experiments will continue to explore, with increased sensitivity, especially extended targets in the sky. This paper reviews the several near-term and longer-term plans for gamma-ray observatories, from MeV energies up to hundreds of TeV. Similarly, the X-ray sky has been and continues to be monitored by decade-old observatories. Upcoming telescopes will further bolster searches and allow new discovery space for lines from, e.g., sterile neutrinos and axion-photon conversion. Furthermore, this overview discusses currently operating cosmic-ray probes and the landscape of future experiments that will clarify existing persistent anomalies in cosmic radiation and spearhead possible new discoveries. Finally, the article closes with a discussion of necessary cross section measurements that need to be conducted at colliders to reduce substantial uncertainties in interpreting photon and cosmic-ray measurements in space.

We present ForSE (Foreground Scale Extender), a novel Python package which aims at overcoming the current limitations in the simulation of diffuse Galactic radiation, in the context of Cosmic Microwave Background experiments (CMB). ForSE exploits the ability of generative adversarial neural networks (GANs) to learn and reproduce complex features present in a set of images, with the goal of simulating realistic and non-Gaussian foreground radiation at sub-degree angular scales. This is of great importance in order to estimate the foreground contamination to lensing reconstruction, de-lensing and primordial B-modes, for future CMB experiments. We applied this algorithm to Galactic thermal dust emission in both total intensity and polarization. Our results show how ForSE is able to generate small scale features (at 12 arc-minutes) having as input the large scale ones (80 arc-minutes). The injected structures have statistical properties, evaluated by means of the Minkowski functionals, in good agreement with those of the real sky and which show the correct amplitude scaling as a function of the angular dimension. The obtained thermal dust Stokes Q and U full sky maps as well as the ForSE package are publicly available for download.

We present a comprehensive joint analysis of two distinct methodologies for measuring the mass of galaxy clusters: hydrostatic measurements and caustic techniques. We show that by including cluster-specific assumptions obtained from hydrostatic measurements in the caustic method, the potential mass bias between these approaches can be significantly reduced. Applying this approach to two well-observed massive galaxy clusters A2029 and A2142. We find no discernible mass bias, affirming the method's validity. We then extend the analysis to modified gravity models and draw a similar conclusion when applying our approach. Specifically, our implementation allows us to investigate Chameleon and Vainshtein screening mechanisms, tightening the posteriors and enhancing our understanding of these modified gravity scenarios.

Future cosmic microwave background (CMB) experiments are primarily targeting a detection of the primordial $B$-mode polarisation. The faintness of this signal requires exquisite control of systematic effects which may bias the measurements. In this work, we derive requirements on the relative calibration accuracy of the overall polarisation gain ($\Delta g_\nu$) for LiteBIRD experiment, through the application of the blind Needlet Internal Linear Combination (NILC) foreground-cleaning method. We find that minimum variance techniques, as NILC, are less affected by gain calibration uncertainties than a parametric approach, which requires a proper modelling of these instrumental effects. The tightest constraints are obtained for frequency channels where the CMB signal is relatively brighter (166 GHz channel, $\Delta {g}_\nu \approx 0.16 \%$), while, with a parametric approach, the strictest requirements were on foreground-dominated channels. We then propagate gain calibration uncertainties, corresponding to the derived requirements, into all frequency channels simultaneously. We find that the overall impact on the estimated $r$ is lower than the required budget for LiteBIRD by almost a factor $5$. The adopted procedure to derive requirements assumes a simple Galactic model. We therefore assess the robustness of obtained results against more realistic scenarios by injecting the gain calibration uncertainties, according to the requirements, into LiteBIRD simulated maps and assuming intermediate- and high-complexity sky models. In this case, we employ the so-called Multi-Clustering NILC (MC-NILC) foreground-cleaning pipeline and obtain that the impact of gain calibration uncertainties on $r$ is lower than the LiteBIRD gain systematics budget for the intermediate-complexity sky model. For the high-complexity case, instead, it would be necessary to tighten the requirements by a factor $1.8$.

We constrain cosmological parameters from a joint cosmic shear analysis of peak-counts and the two-point shear correlation functions, as measured from the Dark Energy Survey (DES-Y1). We find the structure growth parameter $S_8\equiv \sigma_8\sqrt{\Omega_{\rm m}/0.3} = 0.766^{+0.033}_{-0.038}$, which at 4.8% precision, provides one of the tightest constraints on $S_8$ from the DES-Y1 weak lensing data. In our simulation-based method we determine the expected DES-Y1 peak-count signal for a range of cosmologies sampled in four $w$CDM parameters ($\Omega_{\rm m}$, $\sigma_8$, $h$, $w_0$). We also determine the joint covariance matrix with over 1000 realisations at our fiducial cosmology. With mock DES-Y1 data we calibrate the impact of photometric redshift and shear calibration uncertainty on the peak-count, marginalising over these uncertainties in our cosmological analysis. Using dedicated training samples we show that our measurements are unaffected by mass resolution limits in the simulation, and that our constraints are robust against uncertainty in the effect of baryon feedback. Accurate modelling for the impact of intrinsic alignments on the tomographic peak-count remains a challenge, currently limiting our exploitation of cross-correlated peak counts between high and low redshift bins. We demonstrate that once calibrated, a fully tomographic joint peak-count and correlation functions analysis has the potential to reach a 3% precision on $S_8$ for DES-Y1. Our methodology can be adopted to model any statistic that is sensitive to the non-Gaussian information encoded in the shear field. In order to accelerate the development of these beyond-two-point cosmic shear studies, our simulations are made available to the community upon request.

The stochastic gravitational-wave background (SGWB) produced by merging neutron stars exhibits a peak in the kHz band. In this paper, we develop a theoretical framework to exploit this distinctive feature through a Markov Chain Monte Carlo analysis using a simulated dataset of SGWB measurements within this frequency range. The aim is to use the SGWB peak as an observable to constrain a set of astrophysical and cosmological parameters that accurately describe the sources of the SGWB. We examine how variations in these parameters impact the morphology of the SGWB and investigate the necessary sensitivity to effectively constrain them. Given our priors on astrophysical and cosmological parameters, and assuming a power-law integrated sensitivity curve of the order of $10^{-11}$ between 1 kHz and 5 kHz, we show that the values of the chirp mass and common envelope efficiency of the binary systems are retrieved with percent accuracy. Furthermore, the method allows for the reconstruction of the cosmological expansion history populated by these binaries, encompassing the Hubble constant, matter abundance, and the effective equation of state of dark energy.

We explore the capability of measuring lensing signals in $LiteBIRD$ full-sky polarization maps. With a $30$ arcmin beam width and an impressively low polarization noise of $2.16\,\mu$K-arcmin, $LiteBIRD$ will be able to measure the full-sky polarization of the cosmic microwave background (CMB) very precisely. This unique sensitivity also enables the reconstruction of a nearly full-sky lensing map using only polarization data, even considering its limited capability to capture small-scale CMB anisotropies. In this paper, we investigate the ability to construct a full-sky lensing measurement in the presence of Galactic foregrounds, finding that several possible biases from Galactic foregrounds should be negligible after component separation by harmonic-space internal linear combination. We find that the signal-to-noise ratio of the lensing is approximately $40$ using only polarization data measured over $90\%$ of the sky. This achievement is comparable to $Planck$'s recent lensing measurement with both temperature and polarization and represents a four-fold improvement over $Planck$'s polarization-only lensing measurement. The $LiteBIRD$ lensing map will complement the $Planck$ lensing map and provide several opportunities for cross-correlation science, especially in the northern hemisphere.

In this paper we discuss the impact of the s-process nucleosynthesis in Asymptotic Giant Branch stars on the enrichment of heavy elements. We review the main steps made on this subject in the last 40 years and discuss the importance of modelling the evolution of the abundances of such elements in our Milky Way. From the comparison between model results and observations, we can impose strong constraints on stellar nucleosynthesis as well as on the evolution of the Milky Way.

Time-dependent energy spectra of galactic cosmic rays (GCRs) carry fundamental information regarding their origin and propagation. When observed at the Earth, these spectra are significantly affected by the solar wind and the embedded solar magnetic field that permeates the heliosphere, changing significantly over an 11-year solar cycle. Energy spectra of GCRs measured during different epochs of solar activity provide crucial information for a thorough understanding of solar and heliospheric phenomena. The PAMELA experiment had collected data for almost ten years (15th June 2006 - 23rd January 2016), including the minimum phase of solar cycle 23 and the maximum phase of solar cycle 24. In this paper, we present new spectra for helium nuclei measured by the PAMELA instrument from January 2010 to September 2014 over a three Carrington rotation time basis. These data are compared to the PAMELA spectra measured during the previous solar minimum providing a picture of the time dependence of the helium nuclei fluxes over a nearly full solar cycle. Time and rigidity dependencies are observed in the proton-to-helium flux ratios. The force-field approximation of the solar modulation was used to relate these dependencies to the shapes of the local interstellar proton and helium-nuclei spectra.

We implement a sub-resolution prescription for the unresolved dynamical friction onto black holes (BHs) in the OpenGadget3 code. We carry out cosmological simulations of a volume of 16 cMpc3 and zoom-ins of a galaxy group and of a galaxy cluster. The advantages of our new technique are assessed in comparison to commonly adopted methods to hamper spurious BH displacements, i.e. repositioning onto a local minimum of the gravitational potential and ad-hoc boosting of the BH particle dynamical mass. The newly-introduced dynamical friction correction provides centering of BHs on host halos which is at least comparable with the other techniques. It predicts half as many merger events with respect to the repositioning prescription, with the advantage of being less prone to leave sub-structures without any central BH. Simulations featuring our dynamical friction prescription produce a smaller (by up to 50% with respect to repositioning) population of wandering BHs and final BH masses in good agreement with observations. As for individual BH-BH interactions, our dynamical friction model captures the gradual inspiraling of orbits before the merger occurs. By contrast, the repositioning scheme, in its most classical renditions considered, describes extremely fast mergers, while the dynamical mass misrepresents the BHs' dynamics, introducing numerical scattering between the orbiting BHs. Given its performances in describing the centering of BHs within host galaxies and the orbiting of BH pair before their merging, our dynamical friction correction opens interesting applications for an accurate description of the evolution of BH demography within cosmological simulations of galaxy formation at different cosmic epochs and within different environments.