The proper treatment in this scenario is to solve numerically the equation of motion of the PQ field in an expanding Universe, in order to capture all features coming from the misalignment mechanism, including the contribution from topological defects like "axionic" strings and domain walls. An axion mass estimate between 0.05 and 1.50 meV was reported by Borsanyi et al. (2016). The result was calculated by simulating the formation of axions during the post-inflation period on a supercomputer.