Path Integral Monte Carlo calculation of momentum distribution in solid \⁴He

We perform calculations of the momentum distribution n(k) in solid \^4\He by means of path integral Monte Carlo methods. We see that, in perfect crystal, n(k) does not depend on temperature T and that is different from the classical Gaussian shape of Maxwell-Boltzmann distribution, even though these discrepancies decrease when the density of the system increases. In crystal presenting vacancies, we see that for T \ge 0.75K, n(k) presents the same behavior as in the perfect crystal, but, at lower T, it presents a peak when k\to0.