We review the landscape of QCD axion models. Theoretical constructions that extend the window for the axion mass and couplings beyond conventional regions are highlighted and classified. Bounds from cosmology, astrophysics and experimental searches are reexamined and updated.

Recently, Tsupko et al. have put forward the very interesting proposal to use the shadows of high-redshift supermassive black holes (SMBHs) as standard rulers. This would in principle allow us to probe the expansion history within a redshift range which would otherwise be challenging to access. In this note, we critically examine this proposal, and identify a number of important issues which had been previously overlooked. These include difficulties in obtaining reliable SMBH mass estimates and reaching the required angular resolution, and an insufficient knowledge of the accretion dynamics of high-redshift SMBHs. While these issues currently appear to prevent high-redshift SMBH shadows from being used as robust standard rulers, we hope that our flagging them early will help in making this probe theoretically mature by the time it will be experimentally feasible.

We consider the possibility that dark energy and baryons might scatter off each other. The type of interaction we consider leads to a pure momentum exchange, and does not affect the background evolution of the expansion history. We parametrize this interaction in an effective way at the level of Boltzmann equations. We compute the effect of dark energy-baryon scattering on cosmological observables, focusing on the Cosmic Microwave Background (CMB) temperature anisotropy power spectrum and the matter power spectrum. Surprisingly, we find that even huge dark energy-baryon cross-sections $\sigma_{xb} \sim {\cal O}({\rm b})$, which are generically excluded by non-cosmological probes such as collider searches or precision gravity tests, only leave an insignificant imprint on the observables considered. In the case of the CMB temperature power spectrum, the only imprint consists in a sub-percent enhancement or depletion of power (depending whether or not the dark energy equation of state lies above or below $-1$) at very low multipoles, which is thus swamped by cosmic variance. These effects are explained in terms of differences in how gravitational potentials decay in the presence of a dark energy-baryon scattering, which ultimately lead to an increase or decrease in the late-time integrated Sachs-Wolfe power. Even smaller related effects are imprinted on the matter power spectrum. The imprints on the CMB are not expected to be degenerate with the effects due to altering the dark energy sound speed. We conclude that, while strongly appealing, the prospects for a direct detection of dark energy through cosmology do not seem feasible when considering realistic dark energy-baryon cross-sections. As a caveat, our results hold to linear order in perturbation theory.

If the symmetry breaking leading to the origin of the axion dark matter field occurs after the end of inflation and is never restored, then overdensities in the axion field collapse to form dense objects known in the literature as axion miniclusters. The estimates of the typical minicluster mass and radius strongly depend on the details of the cosmology at which the onset of axion oscillations begin. In this work we study the properties and phenomenology of miniclusters in alternative cosmological histories and find that they can change by many orders of magnitude. Our findings have direct implications on current and future experimental searches and, in the case of discovery, could be used to learn something about the universe expansion prior to Big-Bang-Nucleosynthesis.

The last decade witnessed an increasing interest in axions and axion-like particles with many theoretical works published and many new experimental proposals that started a real race towards their discovery. This paper is the Conceptual Design Report of the KLASH (KLoe magnet for Axion SearcH) experiment at the Laboratori Nazionali di Frascati (LNF). The idea of this experiment has been stimulated by the availability of the large volume superconducting magnet, with a moderate magnetic field of 0.6 T, used in the KLOE detector at the DAFNE collider. The main conclusion we draw from this report is the possibility to build and put in operation at LNF in 2-3 years a large haloscope with the sensitivity to KSVZ axions in the low mass range between 0.2 and 1 $\mu$eV, complementary to that of other experiments. Timeline and cost are competitive with respect to other proposals in the same mass region thanks to the availability of most of the infrastructure, in particular the superconducting magnet and the cryogenics plant.

The Standard Model Higgs boson, which has previously been shown to develop an effective vacuum expectation value during inflation, can give rise to large particle masses during inflation and reheating, leading to temporary blocking of the reheating process and a lower reheat temperature after inflation. We study the effects on the multiple stages of reheating: resonant particle production (preheating) as well as perturbative decays from coherent oscillations of the inflaton field. Specifically, we study both the cases of the inflaton coupling to Standard Model fermions through Yukawa interactions as well as to Abelian gauge fields through a Chern-Simons term. We find that, in the case of perturbative inflaton decay to SM fermions, reheating can be delayed due to Higgs blocking and the reheat temperature can decrease by up to an order of magnitude. In the case of gauge-reheating, Higgs-generated masses of the gauge fields can suppress preheating even for large inflaton-gauge couplings. In extreme cases, preheating can be shut down completely and must be substituted by perturbative decay as the dominant reheating channel. Finally, we discuss the distribution of reheat temperatures in different Hubble patches, arising from the stochastic nature of the Higgs VEV during inflation and its implications for the generation of both adiabatic and isocurvature fluctuations.

We consider a quintessence field which transitions from a matter-like to a cosmological constant behavior between recombination and the present time. We aim at easing the tension in the measurement of the present Hubble rate, and we assess the $\Lambda$CDM model properly enlarged to include our quintessence field against cosmological observations. The model does not address the scope we proposed. This result allows us to exclude a class of quintessential models as a solution to the tension in the Hubble constant measurements.

In order to accommodate the QCD axion as the dark matter (DM) in a model in which the Peccei-Quinn (PQ) symmetry is broken before the end of inflation, a relatively low scale of inflation has to be invoked in order to avoid bounds from DM isocurvature fluctuations, $H_* \lesssim \mathcal{O}(10^9)\,$GeV. We construct a simple model in which the Standard Model Higgs field is non-minimally coupled to Palatini gravity and acts as the inflaton, leading to a scale of inflation $H_* \sim 10^8\,$GeV. When the energy scale at which the PQ symmetry breaks is much larger than the scale of inflation, we find that in this scenario the required axion mass for which the axion constitutes all DM is $m_0 \lesssim 0.05{\rm \,\mu eV}$ for a quartic Higgs self-coupling $\lambda_\phi = 0.1$, which correspond to the PQ breaking scale $v_\sigma \gtrsim 10^{14}\,$GeV and tensor-to-scalar ratio $r \sim 10^{-12}$. Future experiments sensitive to the relevant QCD axion mass scale can therefore shed light on the physics of the Universe before the end of inflation.

The Event Horizon Telescope (EHT) collaboration has recently released the first image of a black hole (BH), opening a new window onto tests of general relativity in the strong field regime. In this paper, we derive constraints on the nature of M87* (the supermassive object at the centre of the galaxy M87), exploiting the fact that its shadow appears to be highly circular, and using measurements of its angular size. We first consider the simple case where M87* is assumed to be a Kerr BH. We find that the inferred circularity of M87* excludes Kerr BHs with observation angle $\theta_{\rm obs} \gtrsim 45^{\circ}$ for dimensionless rotational parameter $0.95 \lesssim a_* \leq 1$ whereas the observation angle is unbounded for $a_* \lesssim 0.9$. We then consider the possibility that M87* might be a superspinar, i.e. an object described by the Kerr solution and spinning so fast that it violates the Kerr bound by having $|a_*| > 1$. We find that, within certain regions of parameter space, the inferred circularity and size of the shadow of M87* do not exclude the possibility that this object might be a superspinar.

Persisting tensions between high-redshift and low-redshift cosmological observations suggest the dark energy sector of the Universe might be more complex than the positive cosmological constant of the $\Lambda$CDM model. Motivated by string theory, wherein symmetry considerations make consistent AdS backgrounds (\textit (i.e.) maximally symmetric spacetimes with a negative cosmological constant) ubiquitous, we explore a scenario where the dark energy sector consists of two components: a negative cosmological constant, with a dark energy component with equation of state $w_{\phi}$ on top. We test the consistency of the model against low-redshift Baryon Acoustic Oscillation and Type Ia Supernovae distance measurements, assessing two alternative choices of distance anchors: the sound horizon at baryon drag determined by the \textit{Planck} collaboration, and the Hubble constant determined by the SH0ES program. We find no evidence for a negative cosmological constant, and mild indications for an effective phantom dark energy component on top. A model comparison analysis reveals the $\Lambda$CDM model is favoured over our negative cosmological constant model. While our results are inconclusive, should low-redshift tensions persist with future data, it would be worth reconsidering and further refining our toy negative cosmological constant model by considering realistic string constructions.

The Event Horizon Telescope has recently provided the first image of the dark shadow around the supermassive black hole M87*. The observation of a highly circular shadow provides strong limits on deviations of M87*'s quadrupole moment from the Kerr value. We show that the absence of such a deviation can be used to constrain the physics of extra dimensions of spacetime. Focusing on the Randall-Sundrum AdS$_5$ brane-world scenario, we show that the observation of M87*'s dark shadow sets the limit $\ell \lesssim 170\,{\rm AU}$, where $\ell$ is the AdS$_5$ curvature radius. This limit is among the first quantitative constraints on exotic physics obtained from the extraordinary first ever image of the dark shadow of a black hole.

We show results for the expected reach of the network of experiments that is being set up globally with the aim of detecting the "invisible" axion, in light of a non-standard thermal history of the universe. Assuming that the axion is the dark matter, we discuss the reach of a successful detection by a given experimental setup in a particular axion mass window for different modifications of the cosmological background before primordial nucleosynthesis occurred. Results are presented both in the case where the present energy budget in cold axions is produced through the vacuum realignment mechanism alone, or in the case in which axionic strings also provide with additional contributions to the axion energy density. We also show that in some cosmological models, the spectrum of gravitational waves from the axionic string network would be within reach of the future network of detectors like LISA and DECIGO-BBO. We conclude that some scenarios describing the early universe can be probed jointly by the experimental efforts on axion detection and by gravity wave multi-messenger astronomy.

We consider single-field inflation in light of string-motivated "swampland" conjectures suggesting that effective scalar field theories with a consistent UV completion must have field excursion $\Delta \phi \lesssim M_{\rm Pl}$, in combination with a sufficiently steep potential, $M_{\rm Pl} V_\phi/V \gtrsim {\cal O}(1)$. Here, we show that the swampland conjectures are inconsistent with existing observational constraints on single-field inflation. Focusing on the observationally favoured class of concave potentials, we map the allowed swampland region onto the $n_S$-$r$ "zoo plot" of inflationary models, and find that consistency with the Planck satellite and BICEP2/Keck Array requires $M_{\rm Pl} V_\phi/V \lesssim 0.1$ and $-0.02 \lesssim M_{\rm Pl}^2 V_{\phi\phi}/V < 0$, in strong tension with swampland conjectures. Extension to non-canonical models such as DBI Inflation does not significantly weaken the bound.

In the simplest picture, the masses of string axions populating the axiverse depend on two parameters: the supersymmetry breaking scale $M_{\rm susy}$ and the action $S$ of the string instantons responsible for breaking the axion shift symmetry. In this work, we explore whether cosmological data can be used to probe these two parameters. Adopting string-inspired flat priors on $\log_{10}M_{\rm susy}$ and $S$, and imposing that $M_{\rm susy}$ be sub-Planckian, we find $S=198\pm28$. These bounds suggest that cosmological data complemented with string-inspired priors select a quite narrow axion mass range within the axiverse, $\log_{10}(m_a/{\rm eV}) = -21.5^{+1.3}_{-2.3}$. We find that $M_{\rm susy}$ remains unconstrained due to a fundamental parameter degeneracy with $S$. We explore the significant impact of other choices of priors on the results, and we comment on similar findings in recent previous literature.

Axion-like particles (ALPs) might constitute the totality of the cold dark matter (CDM) observed. The parameter space of ALPs depends on the mass of the particle $m$ and on the energy scale of inflation $H_I$ , the latter being bound by the non-detection of primordial gravitational waves. We show that the bound on HI implies the existence of a mass scale $m_\chi = 10 {\rm \,neV} {\div} 0.5 {\rm \,peV}$, depending on the ALP susceptibility $\chi$, such that the energy density of ALPs of mass smaller than $m_\chi$ is too low to explain the present CDM budget, if the ALP field has originated after the end of inflation. This bound affects Ultra-Light Axions (ULAs), which have recently regained popularity as CDM candidates. Light ($m < m_\chi$) ALPs can then be CDM candidates only if the ALP field has already originated during the inflationary period, in which case the parameter space is constrained by the non-detection of axion isocurvature fluctuations. We comment on the effects on these bounds from additional physics beyond the Standard Model, besides ALPs.

Understanding the nature of the Dark Matter (DM) is one of the current challenges in modern astrophysics and cosmology. Knowing the properties of the DM particle would shed light on physics beyond the Standard Model and even provide us with details of the early Universe. In fact, the detection of such a relic would bring us information from the pre-Big Bang Nucleosynthesis (BBN) period, an epoch from which we have no data, and could even hint at inflationary physics. In this work, we assume that the expansion rate of the Universe after inflationary is governed by the kinetic energy of a scalar field $\phi$, in the so-called "kination" model. We assume that the $\phi$ field decays into both radiation and DM particles, which we take to be Weakly Interacting Massive Particles (WIMPs). The present abundance of WIMPs is then fixed during the kination period through either a thermal "freeze-out" or "freeze-in" mechanism, or through a non-thermal process governed by the decay of $\phi$. We explore the parameter space of this theory with the requirement that the present WIMP abundance provides the correct DM relic budget. Requiring that BBN occurs during the standard cosmological scenario sets a limit on the temperature at which the kination period ends. Using this limit and assuming the WIMP has a mass $m_\chi = 100\,$GeV, we obtain that the thermally-averaged WIMP annihilation cross section has to satisfy the constraints $4 \times 10^{-16}{\rm \,GeV^{-2}} \lesssim \langle\sigma v\rangle \lesssim 2\times 10^{-5} {\rm \,GeV^{-2}}$ in order for having at least one of the production mechanism to yield the observed amount of DM. This result shows how the properties of the WIMP particle, if ever measured, can yield information on the pre-BBN content of the Universe.

We derive constraints on mixed dark-matter scenarios consisting of primordial black holes (PBHs) and weakly interacting massive particles (WIMPs). In these scenarios, we expect a density spike of the WIMPs that are gravitationally bound to the PBHs, which results in an enhanced annihilation rate and increased indirect detection prospects. We show that such scenarios provide strong constraints on the allowed fraction of PBHs that constitutes the dark matter, depending on the WIMP mass $m_{\chi}$ and the velocity-averaged annihilation cross-section $\langle \sigma v \rangle$. For the standard scenario with $m_{\chi} = 100\,{\rm GeV}$ and $\langle \sigma v \rangle = 3 \times 10^{-26}\,{\rm cm}^{3} / {\rm s}$, we derive bounds that are stronger than all existing bounds for PBHs with masses $10^{-12}\,M_{\odot} \lesssim M_{\rm BH} \lesssim 10^{4}\,M_{\odot}$, where $M_{\odot}$ is the solar mass, and mostly so by several orders of magnitude.

The search for extra dimensions is a challenging endeavor to probe physics beyond the Standard Model. The joint detection of gravitational waves (GW) and electromagnetic (EM) signals from the merging of a binary system of compact objects like neutron stars (NS), can help constrain the geometry of extra dimensions beyond our 3+1 spacetime ones. A theoretically well-motivated possibility is that our observable Universe is a 3+1-dimensional hypersurface, or brane, embedded in a higher 4+1-dimensional Anti-de Sitter (AdS$_5$) spacetime, in which gravity is the only force which propagates through the infinite bulk space, while other forces are confined to the brane. In these types of brane-world models, GW and EM signals between two points on the brane would, in general, travel different paths. This would result in a time-lag between the detection of GW and EM signals emitted simultaneously from the same source. We consider the recent near-simultaneous detection of the GW event GW170817 from the LIGO/Virgo collaboration, and its EM counterpart, the short gamma-ray burst GRB170817A detected by the Fermi Gamma-ray Burst Monitor and the INTEGRAL Anti-Coincidence Shield spectrometer. Assuming the standard $\Lambda$-Cold Dark Matter ($\Lambda$CDM) scenario and performing a likelihood analysis which takes into account astrophysical uncertainties associated to the measured time-lag, we set an upper limit of $\ell \lesssim 0.535\,$Mpc at $68\%$ confidence level on the AdS$_5$ radius of curvature $\ell$. Although the bound is not competitive with current Solar System constraints, it is the first time that data from a multi-messenger GW-EM measurement is used to constrain extra-dimensional models. Thus, our work provides a proof-of-principle for the possibility of using multi-messenger astronomy for probing the geometry of our space-time.

Axion stars are hypothetical objects formed of axions, obtained as localized and coherently oscillating solutions to their classical equation of motion. Depending on the value of the field amplitude at the core $|\theta_0| \equiv |\theta(r=0)|$, the equilibrium of the system arises from the balance of the kinetic pressure and either self-gravity or axion self-interactions. Starting from a general relativistic framework, we obtain the set of equations describing the configuration of the axion star, which we solve as a function of $|\theta_0|$. For small $|\theta_0| \lesssim 1$, we reproduce results previously obtained in the literature, and we provide arguments for the stability of such configurations in terms of first principles. We compare qualitative analytical results with a numerical calculation. For large amplitudes $|\theta_0| \gtrsim 1$, the axion field probes the full non-harmonic QCD chiral potential and the axion star enters the {\it dense} branch. Our numerical solutions show that in this latter regime the axions are relativistic, and that one should not use a single frequency approximation, as previously applied in the literature. We employ a multi-harmonic expansion to solve the relativistic equation for the axion field in the star, and demonstrate that higher modes cannot be neglected in the dense regime. We interpret the solutions in the dense regime as pseudo-breathers, and show that the life-time of such configurations is much smaller than any cosmological time scale.

Bulk viscosity leads to negative pressure, which is a key ingredient for successful inflation. We build an inflationary model where the slow-roll of the inflaton field is driven by a viscous component instead of the flat potential commonly used. Since viscosity does not contribute to the energy density of the Universe, the Hubble rate depends only on the kinetic energy of the inflaton. The power spectrum is almost scale-invariant, with a tilt depending on the slow-roll parameters defined within the model.

Weakly Interacting Massive Particles (WIMPs), which are among the best motivated dark matter (DM) candidates, could make up all or only a fraction of the total DM budget. We consider a scenario in which WIMPs are a sub-dominant DM component; such a scenario would affect both current direct and indirect bounds on the WIMP-nucleon scattering cross section. In this paper we focus on indirect searches for the neutrino flux produced by annihilation of sub-dominant WIMPs captured by the Sun or the Earth via either spin-dependent or spin-independent scattering. We derive the annihilation rate and the expected neutrino flux at neutrino observatories. In our computation, we include an updated chemical composition of the Earth with respect to the previous literature, leading to an increase of the Earth's capture rate for spin-dependent scattering by a factor three. Results are compared with current bounds from Super-Kamiokande and IceCube. We discuss the scaling of bounds from both direct and indirect detection methods with the WIMP abundance.

Warm inflation is, as of today, one of the best motivated mechanisms for explaining an early inflationary period. In this paper, we derive and analyze the current bounds on warm inflation with a monomial potential $U\propto \phi^p$, using the constraints from the PLANCK mission. In particular, we discuss the parameter space of the tensor-to-scalar ratio $r$ and the potential coupling $\lambda$ of the monomial warm inflation in terms of the number of e-folds. We obtain that the theoretical tensor-to-scalar ratio $r\sim 10^{-8}$ is much smaller than the current observational constrain $r \lesssim 0.12$, despite a relatively large value of the field excursion $\Delta \phi \sim 0.1\MP$. Warm inflation thus eludes the Lyth bound set on the tensor-to-scalar ratio by the field excursion.

We consider the steady-state regime describing the density profile of a dark matter halo, if dark matter is treated as a Bose-Einstein condensate. We first solve the fluid equation for "canonical" cold dark matter, obtaining a class of density profiles which includes the Navarro-Frenk-White profile, and which diverge at the halo core. We then solve numerically the equation obtained when an additional "quantum pressure" term is included in the computation of the density profile. The solution to this latter case is finite at the halo core, possibly avoiding the "cuspy halo problem" present in some cold dark matter theories. Within the model proposed, we predict the mass of the cold dark matter particle to be of the order of M_chi c2 = 10^-24 eV, which is of the same order of magnitude as that predicted in ultra-light scalar cold dark matter models. Finally, we derive the differential equation describing perturbations in the density and the pressure of the dark matter fluid.

We present a general expression for the values of the average kinetic energy and of the temperature of kinetic decoupling of a WIMP, valid for any cosmological model. We show an example of the usage of our solution when the Hubble rate has a power-law dependence on temperature.

Using a WKB approximation of the Dirac equation in a curved spacetime, we obtain the expression for the phase shift between two neutrino mass eigenstates in a generic gravitational field. We apply this expression to two specific space-time geometries, namely the Kerr-Newman metric describing a rotating and charged black hole, and the Friedmann-Robertson-Walker metric.

We present a general expression for the values of the average kinetic energy and of the temperature of kinetic decoupling of a WIMP, valid for any cosmological model. We show an example of the usage of our solution when the Hubble rate has a power-law dependence on temperature, and we show results for the specific cases of kination cosmology and low- temperature reheating cosmology.

Within the framework of the interacting fluid formalism, we provide the numerical solution to the Boltzmann equation describing the evolution of an inflaton field coupled to radiation. We study the behavior of the system during the slow-roll regime, in the case in which an additional stochastic source term is included in the set of equations, and we recover the expression for the cosmological perturbations previously obtained in the Warm inflation scenarios.

The subjects of this thesis are the invisible axion and the more general family of axion-like particles. The invisible axion is a hypothetical elementary particle and a cold dark matter candidate. I present an improved computation of the constraints on the parameter space of the cold dark matter axion in the standard cosmology, that includes the contributions from anharmonicities in the axion potential and from the decay of axionic strings. In this scenario, I update the value of the mass of the cold dark matter axion, finding the value $(67\pm17){\rm \mu eV}$, approximately one order of magnitude larger than previous computations. The effect of nonstandard cosmological scenarios on the parameter space of axion cold dark matter is studied for the first time. In particular, I consider the cases of low-temperature reheating and kination cosmologies, and I show that the mass of the cold dark matter axion can differ from the value in the standard cosmological scenario by orders of magnitude. Finally, I consider the family of axion-like particles, assuming that these particles serve as the inflaton in the context of warm inflation. I find that the axion energy scale $f$, which in the standard inflation scenario is of the order of the Planck mass, can be lowered to the much safer Grand Unification Theory scale $f \sim 10^{16}{\rm GeV}$. I also constrain the parameter space and the amount of gravitational waves from this model, using results from the Wilkinson Microwave Anisotropy Probe 7-year data.

We derive the requirements that a generic axion-like field has to satisfy in order to play the role of the inflaton field in the warm inflation scenario. Compared to the parameter space in ordinary natural inflation models, we find that the parameter space in our model is enlarged. In particular, we avoid the problem of having an axion decay constant $f$ that relates to the Planck scale, which is instead present in the ordinary natural inflation models; in fact, our model can easily accommodate values of the axion decay constant that lie well below the Planck scale.

We study the parameter space of cold dark matter axions in two cosmological scenarios with non-standard thermal histories before Big Bang nucleosynthesis: the Low Temperature Reheating (LTR) cosmology and the kination cosmology. If the Peccei-Quinn symmetry breaks during inflation, we find more allowed parameter space in the LTR cosmology than in the standard cosmology and less in the kination cosmology. On the contrary, if the Peccei-Quinn symmetry breaks after inflation, the Peccei-Quinn scale is orders of magnitude higher than standard in the LTR cosmology and lower in the kination cosmology. We show that the axion velocity dispersion may be used to distinguish some of these non-standard cosmologies. Thus, axion cold dark matter may be a good probe of the history of the Universe before Big Bang nucleosynthesis.

We study for what specific values of the theoretical parameters the axion can form the totality of cold dark matter. We examine the allowed axion parameter region in the light of recent data collected by the WMAP5 mission plus baryon acoustic oscillations and supernovae, and assume an inflationary scenario and standard cosmology. If the Peccei-Quinn symmetry is restored after inflation, we recover the usual relation between axion mass and density, so that an axion mass $m_a =67\pm2{\rm \mu eV}$ makes the axion 100% of the cold dark matter. If the Peccei-Quinn symmetry is broken during inflation, the axion can instead be 100% of the cold dark matter for $m_a < 15{\rm meV}$ provided a specific value of the initial misalignment angle $\theta_i$ is chosen in correspondence to a given value of its mass $m_a$. Large values of the Peccei-Quinn symmetry breaking scale correspond to small, perhaps uncomfortably small, values of the initial misalignment angle $\theta_i$.

We study for what specific values of the theoretical parameters the axion can form the totality of cold dark matter. We examine the allowed axion parameter region in the light of recent data collected by the WMAP5 mission plus baryon acoustic oscillations and supernovae \cite{komatsu}, and assume an inflationary scenario and standard cosmology. We also upgrade the treatment of anharmonicities in the axion potential, which we find important in certain cases. If the Peccei-Quinn symmetry is restored after inflation, we recover the usual relation between axion mass and density, so that an axion mass $m_a =(85 \pm 3){\rm \mu eV}$ makes the axion 100% of the cold dark matter. If the Peccei-Quinn symmetry is broken during inflation, the axion can instead be 100% of the cold dark matter for $m_a < 15{\rm meV}$ provided a specific value of the initial misalignment angle $\theta_i$ is chosen in correspondence to a given value of its mass $m_a$. Large values of the Peccei-Quinn symmetry breaking scale correspond to small, perhaps uncomfortably small, values of the initial misalignment angle $\theta_i$.