We study the gravitational collapse of axion dark matter fluctuations in the post-inflationary scenario, so-called axion miniclusters, with N-body simulations. Largely confirming theoretical expectations, overdensities begin to collapse in the radiation-dominated epoch and form an early distribution of miniclusters with masses up to $10^{-12}\,M_\odot$. After matter-radiation equality, ongoing mergers give rise to a steep power-law distribution of minicluster halo masses. The density profiles of well-resolved halos are NFW-like to good approximation. The fraction of axion DM in these bound structures is $\sim 0.75$ at redshift $z=100$.

We study the formation and the subsequent mass growth of axion stars inside axion miniclusters. Numerically solving the Schr\"odinger-Poisson equations with realistic initial conditions we find that the axion stars exhibit similar properties to solitonic cores in ultralight bosonic dark matter halos in terms of their radial density profiles and large-amplitude oscillations. A merger of two axion stars confirms a previously found empirical law for the mass of the merged axion star. Monitoring the axion star masses over time, we observe a mass growth consistent with the mass increase of Bose stars in the kinetic regime reported by Levkov et al., confirming that the mass evolution of axion stars can be understood in terms of wave condensation. Based on this result, we predict a saturation of mass growth in relation to the minicluster mass consistent with the core-halo mass relation previously found for ultra-light bosonic dark matter halos.