Simulation of Material Properties Below the Debye Temperature: A Path-Integral Molecular Dynamics Case Study of Quartz

Abstract

Classical and path integral molecular dynamics (PIMD) simulations are used to study alpha-quartz and beta-quartz in a large range of temperatures at zero external stress. PIMD account for quantum fluctuations of atomic vibrations, which can modify material properties at temperatures below the Debye temperature. The difference between classical and quantum mechanical results for bond lengths, bond angles, elastic modulii, and some dynamical properties is calculated and comparison to experimental data is done. Only quantum mechanical simulations are able to reflect the correct thermomechanical properties below room temperature. It is discussed in how far classical and PIMD simulations can be helpful in constructing improved potential energy surfaces for silica.

8 pages, 9 figures, submitted to J. Chem. Phys

References in corpus (1)

• Path integral Monte Carlo simulations of silicates

Cited by in corpus (5)

• Molecular Dynamics Simulations
• Portable implementation of a quantum thermal bath for molecular dynamics simulations
• Elastic Constants of Quantum Solids by Path Integral Simulations
• Comparison of two non-primitive methods for path integral simulations: Higher-order corrections vs. an effective propagator approach
• Computing the thermal transport coefficient of neutral amorphous polymers using exact vibrational density of states: Comparison with experiments