AbstractA crystal is described by a probability density which is not symmetric regarding the interchange of phase‐space coordinates between two molecules. The set of approximate equations for one‐particle distribution functions with many‐body forces is derived from the LIOUVILLEequation. The expression for the HELMHOLTZfree energy and the equation of state are obtained. The developed method is applied to the determination of the thermodynamic properties of crystalline argon, krypton, and xenon. The nearest‐neighbor distances, internal energies, isothermal and adiabatic compressibilities, linear thermal expansion coefficients, specific heatsCvandCpare calculated using the pair potentials of BARKER‐POMPE, of BARKER‐BOBETIC, and of LENNARD‐JONES, together with the AXILROD‐TELLERthree‐body potential. The obtained results are compared with the experimental d