Thermophysical properties of the Lennard-Jones substance, i.e. particles interacting with the Lennard-Jones potential can be obtained using statistical mechanics. Some properties can be computed analytically, i.e. with machine precision, whereas most properties can only be obtained by performing molecular simulations. The latter will in general be superimposed by both statistical and systematic uncertainties. The virial coefficients can for example be computed directly from the Lennard-potential using algebraic expressions and reported data has therefore no uncertainty. Molecular simulation results, e.g. the pressure at a given temperature and density has both statistical and systematic uncertainties. Molecular simulations of the Lennard-Jones potential can in general be performed using either molecular dynamics (MD) simulations or Monte Carlo (MC) simulation. For MC simulations, the Lennard-Jones potential V L J ( r ) is directly used, whereas MD simulations are always based on the derivative of the potential, i.e. the force F = d V / d r . These differences in combination with differences in the treatment of the long-range interactions (see below) can influence computed thermophysical properties.