Microscopic Equations of State, Thermodynamical Properties, and System Collapse: Application of the Dynamical Equation of Period Vectors in Crystal Structure Prediction under External Temperature and Stress/Pressure

In crystal structure prediction, Newton's second law can always be applied to particles (atoms or ions) in a cell to determine their positions. However the values of the period vectors (crystal cell edge vectors) should be determined as well. Here we applied the dynamical equations of period vectors derived recently based on Newton's laws (doi:/10.1139/cjp-2014-0518) for that purpose, where the period vectors are driven by the imbalance between the internal and external pressures/stresses. For equilibrium states, they became equations of state, which essentially turn out the equilibrium conditions of crystals from mechanical point of view. Additionally for external pressure, they became Mie-Gruneisen equation supposing the Gruneisen constant is 1/3, which means phonon frequency is inversely proportional to the cell length. Since the internal stress has both a full kinetic energy term and a full interaction term, the influences of both external temperature and stress/pressure on crystal structures can be calculated, then thermodynamical properties and processes were presented. Contrary to usual ideas, the equations show that pure harmonic vibrating phonons can result in thermal expansion in crystals when the external temperature is changed. Finally, crystal system collapse due to temperature and/or stress/pressure change was discussed.