Francesco D'Eramo

When the axion mass is generated during a first-order phase transition and becomes nonvanishing only inside expanding true-vacuum bubbles, the standard picture of misalignment production is qualitatively modified. Using lattice simulations in an expanding universe, we study dark matter production within such a framework and identify two distinct regimes. For rapid transitions, the onset of oscillations is delayed until bubble percolation, enhancing the relic abundance. For slower transitions, spatial gradients generated by expanding bubbles suppress the effective misalignment angle through the bubble misalignment mechanism. We derive a semianalytical expression for the relic density that provides a unified description of both regimes and accurately reproduces the simulation results. Finally, we show how this mechanism also modifies isocurvature perturbations and the small-scale matter power spectrum, with important implications for axion minicluster formation.

We present a framework in which three classes of dark matter number-changing processes can affect both the relic abundance via thermal freeze-out in the early universe and the generation of indirect cosmic-ray signals today. These processes are: (i) direct annihilations into Standard Model final states; (ii) annihilations into metastable on-shell mediators that subsequently decay into Standard Model particles; (iii) semi-annihilation processes featuring a dark matter particle in the final state, accompanied by a metastable mediator. A central element of our analysis is the systematic inclusion of semi-annihilation alongside the more commonly considered channels. This setup is largely model-independent, as we only assume the presence of one or more of these processes with unsuppressed s-wave contributions. We analyze representative benchmarks for the dominant decay modes of the mediator and show how the resulting injection spectra for γ rays, neutrinos, and cosmic-ray antimatter vary with the relative importance of the three classes of processes. As an application, we evaluate the observable γ-ray fluxes from dwarf spheroidal galaxies in the GeV–TeV window. Finally, we provide explicit model realizations in which multiple processes coexist, and discuss how their interplay shapes indirect detection signatures. Our results provide a consistent connection between early-universe dynamics and present-day observables, revealing distinctive features that arise when multiple dark matter processes contribute simultaneously.

The meV mass range has emerged as a focal point in axion physics, where advances in theory, cosmology, astrophysics, and experimental techniques converge. Axions in this mass range are theoretically well motivated, can arise in ultraviolet-complete models, and can have significant cosmological impacts as dark matter or dark radiation. In parallel, their efficient production in stellar and supernova environments provides powerful astrophysical probes. Here, we provide a comprehensive overview of meV axions across these domains, highlighting both established results and open questions. We discuss the theoretical underpinnings of meV axions, their cosmological and astrophysical signatures, and the diverse experimental strategies — ranging from helioscopes and haloscopes to quasiparticle systems and large-volume Cherenkov detectors — that aim to explore this regime. The convergence of these approaches emphasizes the pivotal role of the meV mass range for axion discovery in the coming years, identifying meV axions as a key probe for testing beyond-Standard-Model physics. This review document is the direct outcome of the discussions at the dedicated workshop “The meV Mass Axion Frontier: Challenges and Opportunities”, held at Laboratori Nazionali di Frascati (IT) on 27–28 October 2025, and organized by the EU funded COST Action “Cosmic WISPers in the Dark Universe: Theory, astrophysics, and experiments” (CA21106, https://www.cost.eu/actions/CA21106). Its aim is to provide an overview of current efforts in meV axion research, their motivations, and the research goals that animate the community involved in this search.

We introduce a novel class of bosonic dark matter candidates that we dub wallions, featuring boundaries in field space. The wallion mass is exponentially suppressed when the separation between boundaries far exceeds their intrinsic width and remains radiatively stable under self-interactions. We study the early universe evolution of wallions and the associated cosmological signatures. Finally, we show that instanton effects can dynamically generate field-space boundaries and discuss possible experimental probes once the wallion couples to standard model fields.

Quantifying the imprints of freeze-in dark matter (DM) on cosmological structures requires knowledge of its phase-space distribution. We investigate how variations in the cosmological history before nucleosynthesis, the “weather” of that epoch, give rise to distinct “seasons” in the DM momentum distribution that governs its warmness. Studying decay-driven production across diverse cosmological histories, we map how these conditions shape DM phase-space properties. Our study quantifies how the early Universe composition plays a key role in determining the mass bound on freeze-in DM.

The early universe grants access to energy scales far beyond those achievable in terrestrial experiments and allows unstable Standard Model particles to play an active dynamical role. In this contribution, we focus on recent studies aimed at quantifying the potential of the early universe to probe the properties and interactions of axions. The discussion is organized around four classes of axion scenarios, ordered from long to short lifetimes: (i) stable or long-lived axions contributing to dark radiation; (ii) stable or long-lived axions produced out-of-equilibrium and constituting dark matter; (iii) metastable axions whose decays inject energy into the primordial plasma and leave observable signatures in the global 21 cm signal; and (iv) very short-lived axions that act only as portals to additional degrees of freedom. Together, these scenarios highlight the interplay between axion phenomenology and early universe cosmology and demonstrate the potential of cosmological data to probe axions over a broad range of masses and lifetimes.

Observations of the hydrogen hyperfine transition through the 21 cm line near the end of the cosmic dark ages provide unique opportunities to probe new physics. In this work, we investigate the potential of the sky-averaged 21 cm signal to constrain metastable particles produced in the early universe that decay at later times, thereby modifying the thermal and ionization history of the intergalactic medium. The study begins by extending previous analyses of decaying dark matter (DM), incorporating back-reaction effects and tightening photon decay constraints down to DM masses as low as 20.4 eV. The focus then shifts to non-minimal dark sectors with multiple interacting components. The analysis covers two key scenarios: a hybrid setup comprising a stable cold DM component alongside a metastable sub-component, and a two-component dark sector of nearly degenerate states with a metastable heavier partner. A general parameterization based on effective mass spectra and fractional densities allows for a model-independent study. The final part presents two explicit realizations: an axion-like particle coupled to photons, and pseudo-Dirac DM interacting via vector portals or electromagnetic dipoles. These scenarios illustrate how 21 cm cosmology can set leading bounds and probe otherwise inaccessible regions of parameter space.

The suppression of cosmological structure at small scales is a key signature of dark matter (DM) produced via freeze-in in the low-mass regime. We present a comprehensive analysis of its impact, incorporating recent constraints from Milky Way satellite counts, strong gravitational lensing with James Webb Space Telescope (JWST) data, and the Lyman- α forest. We adopt a general strategy to translate existing warm dark matter (WDM) bounds into lower mass limits for a broad class of DM candidates characterized by quasithermal phase space distributions. The benefits of this approach include computational efficiency and the ability to explore a wide range of models. We derive model-independent bounds for DM produced via two-body decays, scatterings, and three-body decays, and apply the framework to concrete scenarios such as the Higgs portal, sterile neutrinos, axionlike particles, and the dark photon portal. Results from specific models confirm the validity of the model-independent analysis.

This report, summarising work achieved in the context of the LHC Dark Matter Working Group, investigates the phenomenology of t-channel dark matter models, spanning minimal setups with a single dark matter candidate and mediator to more complex constructions closer to UV-complete models. For each considered class of models, we examine collider, cosmological and astrophysical implications. In addition, we explore scenarios with either promptly decaying or long-lived particles, as well as featuring diverse dark matter production mechanisms in the early universe. By providing a unified analysis framework, numerical tools and guidelines, this work aims to support future experimental and theoretical efforts in exploring t-channel dark matter models at colliders and in cosmology.