We examine whether the Hubble tension, the mismatch between early and late measurements of $H_0$, can be alleviated by ultralight scalar fields in the early universe, and we assess its plausibility within UV physics. Since their energy density needs to rapidly redshift away, we explore decays to massless fields around the era of matter-radiation equality. We highlight a concrete implementation of ultralight pseudo-scalars, axions, that decay to an abelian dark sector. This scenario circumvents major problems of other popular realizations of early universe scalar models in that it uses a regular scalar potential that is quadratic around the minimum, instead of the extreme fine-tuning of many existing models. The idea is that the scalar is initially frozen in its potential until $H\sim m$, then efficient energy transfer from the scalar to the massless field can occur shortly after the beginning of oscillations due to resonance. We introduce an effective fluid model which captures the transition from the frozen scalar phase to the radiation dark sector phase. We perform a fit to a combined Planck 2018, BAO, SH$_0$ES and Pantheon supernovae dataset and find that the model gives $H_0=69.69_{-0.83}^{+0.81}$ km/s/Mpc with $\Delta\chi^2 \approx -8$ compared to $\Lambda$CDM; while inclusions of other data sets may worsen the fit. Importantly, we find that large values of the coupling between fields is required for sufficiently rapid decay: For axion-gauge field models $\phi F\tilde{F}/\Lambda$ it requires $\Lambda\lesssim f/80$, where $2\pi f$ is the field range. We find related conclusions for scalar-scalar models $\sim\phi\,\chi^2$ and for models that utilize perturbative decays. We conclude that these sorts of ultralight scalar models that purport to alleviate the Hubble tension, while being reasonable effective field theories, require features that are difficult to embed within UV physics.

We present a scenario that ameliorates the tuning problems present in models of primordial black hole dark matter from inflation. Our setup employs the advantages of gravitational collapse in a long epoch of early matter domination with reheating temperature $\lesssim 10^{6}~\text{GeV}$. Furthermore, we make use of a string-inspired class of models where the inflaton is identified with a non-compact axion field. In this framework, the presence of multiple local minima in the inflaton potential can be traced back to an approximate discrete shift symmetry. This scenario allows the formation of primordial black holes in the observationally viable range of masses ($M_{\text{PBH}}\sim 10^{-13}M_{\odot}-10^{-16}M_{\odot}$) accounting for all dark matter, and in excellent agreement with the CMB. Crucially, we find a significant reduction in the required tuning of the parameters of the inflationary potential, in contrast to the standard case of primordial black hole formation during radiation domination. However, abundant formation of light PBHs during an early phase of matter domination can be more easily in conflict with evaporation bounds. We discuss how these can be avoided under mild assumptions on the collapsing energy density fraction.

The LIGO observatories can potentially detect stochastic gravitational waves arising from phase transitions which happened in the early universe at temperatures around $T\sim 10^{8}$ GeV. This provides an extraordinary opportunity for discovering the phase transition associated with the breaking of the Peccei-Quinn symmetry, required in QCD axion models. Here we consider the simplest Peccei-Quinn models and study under which conditions a strong first-order phase transition can occur, analyzing its associated gravitational wave signal. To be detectable at LIGO, we show that some supercooling is needed, which can arise either in Coleman-Weinberg-type symmetry breaking or in strongly-coupled models. We also investigate phase transitions that interestingly proceed by first breaking the electroweak symmetry at large scales before tunneling to the Peccei-Quinn breaking vacuum. In this case, the associated gravitational wave signal is more likely to be probed at the proposed Einstein Telescope.

We investigate whether the $4.4\sigma$ tension on $H_0$ between SH$_{0}$ES 2019 and Planck 2018 can be alleviated by a variation of Newton's constant $G_N$ between the early and the late Universe. This changes the expansion rate before recombination, similarly to the addition of $\Delta N_{\rm eff}$ extra relativistic degrees of freedom . We implement a varying $G_N$ in a scalar-tensor theory of gravity, with a non-minimal coupling of the form $(M^2+\beta \phi^2)R$. If the scalar $\phi$ starts in the radiation era at an initial value $\phi_I \approx 0.3 \, M_{Pl}$ and with $\beta\approx-0.8$, a dynamical transition occurs naturally around the epoch of matter-radiation equality and the field evolves towards zero at late times. As a consequence the $H_0$ tension between SH$_{0}$ES (2019) and Planck 2018+BAO decreases, as in $\Delta N_{\rm eff}$ models. However, mostly due to late-time constraints from Post-Newtonian (PN) local gravity, the tension is reduced only to 3.5$\sigma$ level. When including also the SH$_{0}$ES data in the fit, the varying $G_N$ model has $H_0=69.2_{-0.75}^{+0.62}$ and an improvement of $\Delta\chi^2=-3.6$ compared to $\Lambda$CDM, at the cost of 2 extra parameters. This corresponds to a decrease of $7_{-6}^{+3}$ percent in the value of $G_N$ from the radiation era to the present time. For comparison, we update the fit of the $\Delta N_{\rm eff}$ model to the same dataset. We find that the $\Delta N_{\rm eff}$ model performs better than the simplest varying $G_N$ scenario, with $H_0=70_{-0.95}^{+0.93}$ and $\Delta\chi^2=-5.5$. The $\Lambda$CDM limit of the $\Delta N_{\rm eff}$ model is disfavored at slightly more than 2$\sigma$, since $\Delta N_{\rm eff}=0.316_{-0.15}^{+0.15}$.

Oscillons are bound states sustained by self-interactions that appear in rather generic scalar models. They can be extremely long-lived and in the context of cosmology they have a built-in formation mechanism - parametric resonance instability. These features suggest that oscillons can affect the standard picture of scalar ultra-light dark matter (ULDM) models. We explore this idea along two directions. First, we investigate numerically oscillon lifetimes and their dependence on the shape of the potential. We find that scalar potentials that occur in well motivated axion-like models can lead to oscillons that live up to $10^8$ cycles or more. Second, we discuss the observational constraints on the ULDM models once the presence of oscillons is taken into account. For a wide range of axion masses, oscillons decay around or after matter-radiation equality and can thus act as early seeds for structure formation. We also discuss the possibility that oscillons survive up to today. In this case they can most easily play the role of dark matter.

Strongly-coupled theories at the TeV can naturally drive a long period of supercooling in the early universe. Trapped into the deconfined phase, the universe could inflate and cool down till the temperature reaches the QCD strong scale. We show how at these low temperatures QCD effects are important and could trigger the exit from the long supercooling era. We also study the implications on relic abundances. In particular, the latent heat released at the end of supercooling could be the reason for the similarities between dark matter and baryon energy densities. The axion abundance could also be significantly affected, allowing for larger values of the axion decay constant. Finally, we discuss how a long supercooling epoch could lead to an enhanced gravitational wave signal.

We propose a mechanism to generate Primordial Black Holes (PBHs) which is independent of cosmological inflation and occurs slightly below the QCD phase transition. Our setup relies on the collapse of long-lived string-domain wall networks and is naturally realized in QCD axion models with domain wall number $N_{DW}>1$ and Peccei-Quinn symmetry broken after inflation. In our framework, dark matter is mostly composed of axions in the meV mass range along with a small fraction, $\Omega_{\text{PBH}} \gtrsim 10^{-6} \Omega_{\text{CDM}} $ of heavy $M \sim 10^4-10^7 M_\odot$ PBHs. The latter could play a role in alleviating some of the shortcomings of the $\Lambda$CDM model on sub-galactic scales. The scenario has distinct signatures in ongoing axion searches as well as gravitational wave observatories.

We propose a framework in which the QCD axion has an exponentially large coupling to photons, relying on the "clockwork" mechanism. We discuss the impact of present and future axion experiments on the parameter space of the model. In addition to the axion, the model predicts a large number of pseudo-scalars which can be light and observable at the LHC. In the most favorable scenario, axion Dark Matter will give a signal in multiple axion detection experiments and the pseudo-scalars will be discovered at the LHC, allowing to determine most of the parameters of the model.

Axions with broken discrete shift symmetry (axion monodromy) have recently played a central role both in the discussion of inflation and the `relaxion' approach to the hierarchy problem. We suggest a very minimalist way to constrain such models by the weak gravity conjecture for domain walls: While the electric side of the conjecture is always satisfied if the cosine-oscillations of the axion potential are sufficiently small, the magnetic side imposes a cutoff, $\Lambda^3 \sim m f M_{pl}$, independent of the height of these `wiggles'. We compare our approach with the recent related proposal by Ibanez, Montero, Uranga and Valenzuela. We also discuss the non-trivial question which version, if any, of the weak gravity conjecture for domain walls should hold. In particular, we show that string compactifications with branes of different dimensions wrapped on different cycles lead to a `geometric weak gravity conjecture' relating volumes of cycles, norms of corresponding forms and the volume of the compact space. Imposing this `geometric conjecture', e.g.~on the basis of the more widely accepted weak gravity conjecture for particles, provides at least some support for the (electric and magnetic) conjecture for domain walls.

Large field inflation is arguably the simplest and most natural variant of slow-roll inflation. Axion monodromy may be the most promising framework for realising this scenario. As one of its defining features, the long-range polynomial potential possesses short-range, instantonic modulations. These can give rise to a series of local minima in the post-inflationary region of the potential. We show that for certain parameter choices the inflaton populates more than one of these vacua inside a single Hubble patch. This corresponds to a dynamical phase decomposition, analogously to what happens in the course of thermal first-order phase transitions. In the subsequent process of bubble wall collisions, the lowest-lying axionic minimum eventually takes over all space. Our main result is that this violent process sources gravitational waves, very much like in the case of a first-order phase transition. We compute the energy density and peak frequency of the signal, which can lie anywhere in the mHz-GHz range, possibly within reach of next-generation interferometers. We also note that this "dynamical phase decomposition" phenomenon and its gravitational wave signal are more general and may apply to other inflationary or reheating scenarios with axions and modulated potentials.

We continue the development of axion monodromy inflation, focussing in particular on the backreaction of complex structure moduli. In our setting, the shift symmetry comes from a partial large complex structure limit of the underlying type IIB orientifold or F-theory fourfold. The coefficient of the inflaton term in the superpotential has to be tuned small to avoid conflict with Kahler moduli stabilisation. To allow such a tuning, this coefficient necessarily depends on further complex structure moduli. At large values of the inflaton field, these moduli are then in danger of backreacting too strongly. To avoid this, further tunings are necessary. In weakly coupled type IIB theory at the orientifold point, implementing these tunings appears to be difficult if not impossible. However, fourfolds or models with mobile D7-branes provide enough structural freedom. We calculate the resulting inflaton potential and study the feasibility of the overall tuning given the limited freedom of the flux landscape. Our preliminary investigations suggest that, even imposing all tuning conditions, the remaining choice of flux vacua can still be large enough for such models to provide a promising path to large-field inflation in string theory.

We present a new model of large field inflation along a winding trajectory in the field space of two axionic fields, where the 'axions' originate from the complex structure moduli sector of a Calabi-Yau 3-fold at large complex structure. The winding trajectory arises from fixing one combination of axions by bulk fluxes and allows for a transplanckian effective field range. The inflaton potential arises from small 'instantonic' corrections to the geometry and realises natural inflation. By working in a regime of large complex structure for two complex structure moduli the inflaton potential can be made subdominant without severe tuning. We also discuss the impact of the recent 'no-go theorems' for transplanckian axion periodicities on our work. Interestingly, our setup seems to realise a loophole pointed out in arXiv:1503.00795 and arXiv:1503.04783: our construction is a candidate for a string theory model of large field inflation which is consistent with the mild form of the weak gravity conjecture for axions.