A QUICKSTEP-based quantum mechanics/molecular mechanics approach for silica

Area of Science:

• Computational materials science
• Quantum chemistry
• Solid-state physics

Background:

• Quantum mechanics/molecular mechanics (QM/MM) methods are vital for simulating localized properties in silica systems.
• Accurate modeling of defects like oxygen vacancies in silica is crucial for understanding material behavior.
• Existing QM/MM frameworks require efficient schemes for large-scale simulations.

Purpose of the Study:

• To present a novel QM/MM scheme tailored for silica systems.
• To integrate this scheme into a large-scale molecular dynamics simulation framework (QUICKSTEP).
• To validate the scheme by studying an oxygen vacancy in alpha-quartz.

Main Methods:

• Development of a QM/MM approach for silica, utilizing the QUICKSTEP framework.
• Application of the scheme to an oxygen vacancy in alpha-quartz.
• Calculation of structural properties (Si-Si bond length, formation energy) and dynamical properties (Si-Si bond frequency).

Main Results:

• Good convergence for Si-Si bond length and formation energy was achieved with a quantum cluster of only eight atoms.
• The QM/MM scheme demonstrated suitability for molecular dynamics simulations.
• The Si-Si bond frequency was successfully evaluated using the velocity-velocity correlation function.

Conclusions:

• The developed QM/MM scheme provides an efficient and accurate method for simulating defects in silica.
• A small quantum cluster size (eight atoms) is sufficient for reliable calculations of key properties.
• The method is well-suited for large-scale molecular dynamics simulations of silica-based materials.