Examination of structural stability and phase transitions in constant-pressure first-principles molecular dynamics simulations

Constant-pressure first-principles molecular dynamics (FPMD) simulation is a powerful tool for investigations of structures in crystals. However, it needs enourmous computations so that highly accurate calculations for electronic states cannot be employed at present. In this report, we examined the reliability and applicability of constant-pressure FPMD in the study of structural properties under this limitation. Crystalline silicon was employed as a benchmark to perform constant-pressure FPMD simulations (with a deformable simulation cell). It is found that, in high pressure (metallic) phases, crystalline symmetry is broken with the present simulation conditions. Several structural transformations were realized by compression and decompression, but they are not entirely consistent with experiment. We discuss this discrepancy and conclude that the number of k point sampling in the Brillouin zone is crucial. It is recommended that constant-pressure FPMD is employed to explore candidate structures for unknown solid phases at present computational resources.