A 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 LIOUVILLE equation. The expression for the HELMHOLTZ free 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 heats Cv and Cp are calculated using the pair potentials of BARKER‐POMPE, of BARKER‐BOBETIC, and of LENNARD‐JONES, together with the AXILROD‐TELLER three‐body potential. The obtained results are compared with the experimental data and discussed. Copyright © 1974 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim