Thanks and happy new year! Considering the classical molecular dynamics simulation of copper or aluminum crystal, I notice in most simulation papers, thermodynamic ensemble is defined and thermodynamic constraints are applied. I am just wondering what these constraints physically mean. Does it mean a crystal with right density of vacancies and interstitial sites atoms and other excitation?