Large curvature perturbations can source an observable amount of stochastic gravitational wave background (SGWB). We consider several scenarios where small-scale curvature perturbations are naturally enhanced due to the presence of additional spectator fields during inflation. The same spectator fields can lead to a period of early matter domination (EMD) after inflation. We compute the inflationary spectrum of curvature perturbation and determine its evolution at later times, taking into account both the onset and the end of the EMD epoch, and also the impact of relative velocity perturbation between matter and radiation. The feature that the same field is responsible for both enhanced perturbations and the EMD era, leads to a predictive framework within which the full frequency dependence of SGWB can be computed. The SGWB can be observed in several detectors, including those focused on the nano-Hz regime. Our numerical framework can also be used to study other non-standard cosmological histories.

Axion strings are horizon-size topological defects that may be produced in the early Universe. Ultra-light axion-like particles may form strings that persist to temperatures below that of big bang nucleosynthesis. Such strings have been considered previously as sources of gravitational waves and cosmic microwave background (CMB) polarization rotation. In this work we show, through analytic arguments and dedicated adaptive mesh refinement cosmological simulations, that axion strings deposit a sub-dominant fraction of their energy into high-energy Standard Model (SM) final states, for example, by the direct production of heavy radial modes that subsequently decay to SM particles. This high-energy SM radiation is absorbed by the primordial plasma, leading to novel signatures in precision big bang nucleosynthesis, the CMB power spectrum, and gamma-ray surveys. In particular, we show that CMB power spectrum data constrains axion strings with decay constants $f_a \lesssim 10^{12}$ GeV, up to model dependence on the ultraviolet completion, for axion masses $m_a \lesssim 10^{-29}$ eV; future CMB surveys could find striking evidence of axion strings with lower decay constants.

Effective field theories (EFTs) of heavy particles coupled to the inflaton are rife with operator redundancies, frequently obscured by sensitivity to both boundary terms and field redefinitions. We initiate a systematic study of these redundancies by establishing a minimal operator basis for an archetypal example, the abelian gauge-Higgs-inflaton EFT. Working up to dimension 9, we show that certain low-dimensional operators are entirely redundant and identify new non-redundant operators with potentially interesting cosmological collider signals. Our methods generalize straightforwardly to other EFTs of heavy particles coupled to the inflaton.

Multiple fields can become dynamical during the inflationary epoch. We consider an example where a light field acquires isocurvature fluctuations during inflation and contributes to the dark matter abundance at late times. Interactions between the light field and the adiabatic sector contribute to mixed adiabatic-isocurvature non-Gaussianity (NG). We show the resulting form of NG has a different kinematic dependence than the 'local shape' commonly considered, and highlight the parameter space where a dedicated search is expected to significantly improve the current Planck sensitivity. We interpret our results in the context of the QCD axion and illustrate how the proposed NG searches can improve upon the existing searches for isocurvature power spectrum and bispectrum.

We consider first order cosmological phase transitions (PT) happening at late times, below Standard Model (SM) temperatures $T_{\rm PT} \lesssim$ GeV. The inherently stochastic nature of bubble nucleation and the finite number of bubbles associated with a late-time PT lead to superhorizon fluctuations in the PT completion time. We compute how such fluctuations eventually source curvature fluctuations with universal properties, independent of the microphysics of the PT dynamics. Using Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) measurements, we constrain the energy released in a dark-sector PT. For 0.1 eV $\lesssim T_{\rm PT} \lesssim$ keV this constraint is stronger than both the current bound from additional neutrino species $\Delta N_{\rm eff}$, and in some cases, even CMB-S4 projections. Future measurements of CMB spectral distortions and pulsar timing arrays will also provide competitive sensitivity for keV $\lesssim T_{\rm PT} \lesssim$ GeV.

We present a novel mechanism for gravitational wave generation in the early Universe. Light spectator scalar fields during inflation can acquire a blue-tilted power spectrum due to stochastic effects. We show that this effect can lead to large curvature perturbations at small scales (induced by the spectator field fluctuations) while maintaining the observed, slightly red-tilted curvature perturbations at large cosmological scales (induced by the inflaton fluctuations). Along with other observational signatures, such as enhanced dark matter substructure, large curvature perturbations can induce a stochastic gravitational wave background (SGWB). The predicted strength of SGWB in our scenario, $\Omega_{\rm GW}h^2 \simeq 10^{-20} - 10^{-15}$, can be observed with future detectors, operating between $10^{-5}$ Hz and 10 Hz. We note that, in order to accommodate the newly reported NANOGrav observation, one could consider the same class of spectator models. At the same time, one would need to go beyond the simple benchmark considered here and consider a regime in which a misalignment contribution is also important.

Non-adiabatic production of massive particles is a generic feature of many inflationary mechanisms. If sufficiently massive, these particles can leave features in the cosmic microwave background (CMB) that are not well-captured by traditional correlation function analyses. We consider a scenario in which particle production occurs only in a narrow time-interval during inflation, eventually leading to CMB hot- or coldspots with characteristic shapes and sizes. Searching for such features in CMB data is analogous to searching for late-Universe hot- or coldspots, such as those due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Exploiting this data-analysis parallel, we perform a search for particle-production hotspots in the Planck PR4 temperature dataset, which we implement via a matched-filter analysis. Our pipeline is validated on synthetic observations and found to yield unbiased constraints on sufficiently large hotspots across ≈ 60% of the sky. After removing point sources and tSZ clusters, we find no evidence for new physics and place novel bounds on the coupling between the inflaton and massive particles. These bounds are strongest for larger hotspots, produced early in inflation, whilst sensitivity to smaller hotspots is limited by noise and beam effects. Through such methods we can constrain particles with masses 𝒪(100) times larger than the inflationary Hubble scale, which represents possibly the highest energies ever directly probed with observational data.

The quantum chromodynamics (QCD) axion may solve the strong CP problem and explain the dark matter (DM) abundance of our Universe. The axion was originally proposed to arise as the pseudo-Nambu Goldstone boson of global $\mathrm{U}(1)_{\rm PQ}$ Peccei-Quinn (PQ) symmetry breaking, but axions also arise generically in string theory as zero modes of higher-dimensional gauge fields. In this work we show that string theory axions behave fundamentally differently from field theory axions in the early Universe. Field theory axions may form axion strings if the PQ phase transition takes place after inflation. In contrast, we show that string theory axions do not generically form axion strings. In special inflationary paradigms, such as D-brane inflation, string theory axion strings may form; however, their tension is parametrically larger than that of field theory axion strings. We then show that such QCD axion strings overproduce the DM abundance for all allowed QCD axion masses and are thus ruled out, except in scenarios with large warping. A loop-hole to this conclusion arises in the axiverse, where an axion string could be composed of multiple different axion mass eigenstates; a heavier eigenstate could collapse the network earlier, allowing for the QCD axion to produce the correct DM abundance and also generating observable gravitational wave signals.

We present the results of a search for heavy QCD axions performed by the ArgoNeuT experiment at Fermilab. We search for heavy axions produced in the NuMI neutrino beam target and absorber decaying into dimuon pairs, which can be identified using the unique capabilities of ArgoNeuT and the MINOS near detector. This decay channel is motivated by a broad class of heavy QCD axion models that address the strong CP and axion quality problems with axion masses above the dimuon threshold. We obtain new constraints at a 95% confidence level for heavy axions in the previously unexplored mass range of 0.2-0.9 GeV, for axion decay constants around tens of TeV.

We describe a new mechanism that gives rise to dissipation during cosmic inflation. In the simplest implementation, the mechanism requires the presence of a massive scalar field with a softly-broken global U(1) symmetry, along with the inflaton field. Particle production in this scenario takes place on parametrically sub-horizon scales, at variance with the case of dissipation into gauge fields. Consequently, the backreaction of the produced particles on the inflationary dynamics can be treated in a local manner, allowing us to compute their effects analytically. We determine the parametric dependence of the power spectrum which deviates from the usual slow-roll expression. Non-Gaussianities are always sizeable whenever perturbations are generated by the noise induced by dissipation: f NL eq ≳ O (10).

Positivity bounds represent nontrivial limitations on effective field theories (EFTs) if those EFTs are to be completed into a Lorentz-invariant, causal, local, and unitary framework. While such positivity bounds have been applied in a wide array of physical contexts to obtain useful constraints, their application to inflationary EFTs is subtle since Lorentz invariance is spontaneously broken during cosmic inflation. One path forward is to employ a Breit parameterization to ensure a crossing-symmetric and analytic S-matrix in theories with broken boosts. We extend this approach to a theory with multiple fields, and uncover a fundamental obstruction that arises unless all fields obey a dispersion relation that is approximately lightlike. We then apply the formalism to various classes of inflationary EFTs, with and without isocurvature perturbations, and employ this parameterization to derive new positivity bounds on such EFTs. For multifield inflation, we also consider bounds originating from the generalized optical theorem and demonstrate how these can give rise to stronger constraints on EFTs compared to constraints from traditional elastic positivity bounds alone. We compute various shapes of non-Gaussianity (NG), involving both adiabatic and isocurvature perturbations, and show how the observational parameter space controlling the strength of NG can be constrained by our bounds.

We study the prospect of measuring the electroweak quantum numbers of beyond the Standard Model (SM) color sextet particles that decay into same-sign top quark pairs. Among these particles, the color sextet scalars give rise to top quarks with the same chirality, while the top quarks coming from the color sextet vector would have opposite chirality. This difference gets encoded in the angular distributions of the bottom quarks and leptons originating from the decays of the top quarks. We utilize this feature and the energy distributions of the final state jets and leptons to distinguish among the three possible color sextet resonances, taking into account various SM background processes at the $13$ TeV LHC.

Identifying the anisotropies in a cosmologically sourced stochastic gravitational wave background (SGWB) would be of significance in shedding light on the nature of primordial inhomogeneities. For example, if SGWB carries isocurvature fluctuations, it would provide evidence for a multi-field inflationary origin of these inhomogeneities. However, this is challenging in practice due to finite detector sensitivity and also the presence of the astrophysical foregrounds that can compete with the cosmological signal. In this work, we explore the prospects for measuring cosmological SGWB anisotropies in the presence of an astrophysical counterpart and detector noise. To illustrate the main idea, we perform a Fisher analysis using a well-motivated cosmological SGWB template corresponding to a first order phase transition, and an astrophysical SGWB template corresponding to extra-galactic binary mergers, and compute the uncertainty with which various parameters characterizing the isotropic and anisotropic components can be extracted. We also discuss some subtleties and caveats involving shot noise in the astrophysical foreground. Overall, we show that upcoming experiments, e.g., LISA, Taiji, Einstein Telescope, Cosmic Explorer, and BBO, can all be effective in discovering plausible anisotropic cosmological SGWBs.

Heavy QCD axions are well-motivated extensions of the QCD axion that address the quality problem while still solving the strong CP problem. Owing to the gluon coupling, critical for solving the strong CP problem, these axions can be produced in significant numbers in beam dump and collider environments for axion decay constants as large as PeV, relevant for addressing the axion quality problem. In addition, if these axions have leptonic couplings, they can give rise to long-lived decay into lepton pairs, in particular, dominantly into muons above the dimuon threshold and below the GeV scale in a broad class of axion models. Considering existing constraints, primarily from rare meson decays, we demonstrate that current and future neutrino facilities and long-lived particle searches have the potential to probe significant parts of the heavy QCD axion parameter space via dimuon final states.

The quantum chromodynamics axion with a decay constant near the Grand Unification (GUT) scale has an ultralight mass near a neV. We show, however, that axionlike particles with masses near the keV-PeV range with GUT-scale decay constants are also well motivated in that they naturally arise from axiverse theories with dark nonabelian gauge groups. We demonstrate that the correct dark matter abundance may be achieved by the heavy axions in these models through the misalignment mechanism in combination with a period of early matter domination from the long-lived dark glueballs of the same gauge group. Heavy axion dark matter may decay to two photons, yielding mono-energetic electromagnetic signatures that may be detectable by current or nextgeneration space-based telescopes. We project the sensitivity of next-generation telescopes including Athena, AMEGO, and e-ASTROGAM to such decaying axion dark matter. If the dark sector contains multiple confining gauge groups, then the observed primordial baryon asymmetry may also be achieved in this scenario through spontaneous baryogenesis. We present explicit orbifold constructions where the dark gauge groups unify with the SM at the GUT scale and axions emerge as the fifth components of dark gauge fields with bulk Chern-Simons terms.

Features in the inflationary landscape can inject extra energies to inflation models and produce on-shell particles with masses much larger than the Hubble scale of inflation. This possibility extends the energy reach of the program of cosmological collider physics, in which signals associated with these particles are generically Boltzmann-suppressed. We study the mechanisms of this classical cosmological collider in two categories of primordial features. In the first category, the primordial feature is classical oscillation, which includes the case of coherent oscillation of a massive field and the case of oscillatory features in the inflationary potential. The second category includes any sharp feature in the inflation model. All these classical features can excite unsuppressed quantum modes of other heavy fields which leave observational signatures in primordial non-Gaussianities, including the information about the particle spectra of these heavy degrees of freedom.

We investigate the physics case for a dedicated trigger-level analysis for very low mass diphoton resonances at ATLAS and CMS and introduce a new photon isolation criterion. This greatly increases the acceptance for light particles at the expense of a somewhat larger Standard Model background. We show how such an analysis would likely yield new experimental coverage for axion-like particles for masses below 15 GeV, bridging the gap with the region covered by flavor factories.

Dark radiation (DR) appears as a new physics candidate in various scenarios beyond the Standard Model. While it is often assumed that perturbations in DR are adiabatic, they can easily have an isocurvature component if more than one field was present during inflation, and whose decay products did not all thermalize with each other. By implementing the appropriate isocurvature initial conditions (IC), we derive the constraints on both uncorrelated and correlated DR density isocurvature perturbations from the full Planck 2018 data alone, and also in combination with other cosmological data sets. Our study on free-streaming DR (FDR) updates and generalizes the existing bound on neutrino density isocurvature perturbations by including a varying number of relativistic degrees of freedom, and for coupled DR (CDR) isocurvature, we derive the first bound. We also show that for CDR qualitatively new physical effects arise compared to FDR. One such effect is that for isocurvature IC, FDR gives rise to larger CMB anisotropies compared to CDR -- contrary to the adiabatic case. More generally, we find that a blue-tilt of DR isocurvature spectrum is preferred. This gives rise to a larger value of the Hubble constant $H_0$ compared to the standard $\Lambda$CDM+$\Delta N_{\rm eff}$ cosmology with adiabatic spectra and relaxes the $H_0$ tension.

A bstract Heavy particles with masses much bigger than the inflationary Hubble scale H * , can get non-adiabatically pair produced during inflation through their couplings to the inflaton. If such couplings give rise to time-dependent masses for the heavy particles, then following their production, the heavy particles modify the curvature perturbation around their locations in a time-dependent and scale non-invariant manner. This results into a non-trivial spatial profile of the curvature perturbation that is preserved on superhorizon scales and eventually generates localized hot or cold spots on the CMB. We explore this phenomenon by studying the inflationary production of heavy scalars and derive the final temperature profile of the spots on the CMB by taking into account the subhorizon evolution, focusing in particular on the parameter space where pairwise hot spots (PHS) arise. When the heavy scalar has an $$ \mathcal{O} $$ O (1) coupling to the inflaton, we show that for an idealized situation where the dominant background to the PHS signal comes from the standard CMB fluctuations themselves, a simple position space search based on applying a temperature cut, can be sensitive to heavy particle masses M 0 /H * ∼ $$ \mathcal{O} $$ O (100). The corresponding PHS signal also modifies the CMB power spectra and bispectra, although the corrections are below (outside) the sensitivity of current measurements (searches).

While the QCD axion is often considered to be necessarily light (less than or similar to eV), recent work has opened a viable and interesting parameter space for heavy axions, which solve both the strong CP and the axion quality problems. These well-motivated heavy axions, as well as the generic axionlike particles, call for explorations in the GeV mass realm at collider and beam dump environments. The primary upcoming neutrino experiment, Deep Underground Neutrino Experiment, is simultaneously also a powerful beam dump experiment, enabled by its multipurpose near detector complex. In this study, we show with detailed analyses that the Deep Underground Neutrino Experiment near detector has a unique sensitivity to heavy axions for masses between 20 MeV and 2 GeV, complementary to other future experiments.