Thermal analysis by classical molecular dynamics simulations is discussed on hand of heat capacity of crystals of 9600 atoms. The differences between quantum mechanical and classical mechanical calculations are shown. Anharmonicity is proven to be an important factor. Finally, it is found that defects contribute to an increase in heat capacity before melting. The energy of conformational gauche defects within the crystal is only about 10% due to internal rotation. The other energy must be generated by cooperative strain. The conclusion is that the next generation of faster computers may permit wider use of molecular dynamics simulations in support of the interpretation of thermal analysis.