Dark matter comprised of axion-like particles (ALPs) generated by the realignment mechanism in the post-inflationary scenario leads to primordial isocurvature fluctuations. The power spectrum of these fluctuations is flat for small wave numbers, extending to scales accessible with cosmological surveys. We use the latest measurements of Cosmic Microwave Background (CMB) primary anisotropies together with CMB lensing, Baryonic Acoustic Oscillations (BAO) and Sunyaev Zel'dovich (SZ) cluster counts to measure the amplitude and tilt of the isocurvature component. We find preference for a white-noise isocurvature component in the CMB primary anisotropies; this conclusion is, however, weakened by current large-scale structure (LSS) data. Interpreting the result as a conservative upper limit on the isocurvature component, the combined bound on the ALP mass from all probes is $m_{a} \gtrsim 10^{-19}$ eV, with some dependence on how $m_{a}$ evolves with temperature. The expected sensitivity of cosmic shear and galaxy clustering from future LSS experiments and CMB lensing suggests improved bounds of $m_{a} \gtrsim 10^{-18}$-$10^{-13}$ eV, depending on scale cuts used to avoid non-linearities and the ALP mass-temperature dependence.

We assume that dark matter is comprised of axion-like particles (ALPs) generated by the realignment mechanism in the post-inflationary scenario. This leads to isocurvature fluctuations with an amplitude of order one for scales comparable to the horizon at the time when the ALP field starts oscillating. The power spectrum of these fluctuations is flat for small wave numbers, extending to scales relevant for cosmological observables. Denoting the relative isocurvature amplitude at $k_*$ = 0.05 Mpc$^{-1}$ by $f_{\rm iso}$, Planck observations of the cosmic microwave background (CMB) yield $f_{\rm iso}$ < 0.31 at the 2$\sigma$-level. This excludes the hypothesis of post-inflationary ALP dark matter with masses $m_{a}$ < 10$^{-20}-$10$^{-16}$ eV, where the range is due to details of the ALP mass-temperature dependence. Future CMB stage IV and 21-cm intensity mapping experiments may improve these limits by 1$-$2 orders of magnitude in $m_{a}$.

When Peccei-Quinn (PQ) symmetry breaking happens after inflation, the axion field takes random values in causally disconnected regions. This leads to fluctuations of order one in the axion energy density around the QCD epoch. These over-densities eventually decouple from the Hubble expansion and form so-called miniclusters. We present a semi-analytical method to calculate the average axion energy density, as well as the power spectrum, from the re-alignment mechanism in this scenario. Furthermore, we develop a modified Press & Schechter approach, suitable to describe the collapse of non-linear density fluctuations during radiation domination, which is relevant for the formation of axion miniclusters. It allows us to calculate the double differential distribution of gravitationally collapsed miniclusters as a function of their mass and size. For instance, assuming a PQ scale of $10^{11}$ GeV, minicluster masses range from about $5 \times 10^{-16}$ to $3 \times 10^{-13}$ solar masses and have sizes from about $4\times 10^4$ to $7\times 10^5$ km at the time they start to collapse.