First Principles Calculation of Terahertz Spectroscopy of Molecular Crystals: the crucial role of periodic boundary conditions benchmarked with experimental L-ascorbic acid spectra
Ying Wang,Huifang Ma,Yanzhao Yang,Jiantao Qi,Guiming Zhang,Hao Ren,Wenyue Guo
DOI: https://doi.org/10.1039/d2cp05688g
IF: 3.3
2023-04-05
Physical Chemistry Chemical Physics
Abstract:The terahertz (THz) region vibration spectral signatures of molecular crystals can usually be ascribed to the low-frequency vibrational modes related to weak intermolecular interactions, e.g. van der Waals (vdW) interactions or hydrogen bonding. These interactions dictate the compositional units deviate from their equilibrium configurations collectively. The collective movement are intrinsically long-range, hence the boundary condition used for theoretical calculation would affect the corresponding potential energy gradients and alter the vibrational features. In this work, we constructed a series of finite-sized cluster models with varying sizes and an extended periodic crystal model for the L-ascorbic acid (L-AA) crystal. Density functionals with both semi-local contributions and nonlocal vdW terms, implemented with either atom-centered Gaussian basis or plane waves, were tested. By comparing first principles calculations with experimental time-domain spectra (TDS), we found that the non-local vdW functional opt-B88 combined with a periodic boundary condition is capable to assign all the experimental features in the 0.2~1.6 THz region. Calculations with cluster models failed in this task. Even worse, the deficiency of the cluster models varies upon the cluster sizes, and did not converge as the cluster size grows. Our results demonstrate that an appropriate periodic boundary condition is essential to correctly assign and analyze THz vibration spectra of molecular crystals.
chemistry, physical,physics, atomic, molecular & chemical