Journal
JOURNAL OF COMPUTATIONAL CHEMISTRY
Volume 25, Issue 6, Pages 888-897Publisher
WILEY-BLACKWELL
DOI: 10.1002/jcc.20019
Keywords
CRYSTAL code; vibrational frequencies; crystalline compounds; a-quartz
Categories
Ask authors/readers for more resources
The problem of numerical accuracy in the calculation of vibrational frequencies of crystalline compounds from the hessian matrix is discussed with reference to alpha-quartz (SiO2) as a case study and to the specific implementation in the CRYSTAL code. The Hessian matrix is obtained by numerical differentiation of the analytical gradient of the energy with respect to the atomic positions. The process of calculating vibrational frequencies involves two steps: the determination of the equilibrium geometry, and the calculation of the frequencies themselves. The parameters controlling the truncation of the Coulomb and exchange series in Hartree-Fock, the quality of the grid used for the numerical integration of the Exchange-correlation potential in Density Functional Theory, the SCF convergence criteria, the parameters controlling the convergence of the optimization process as well as those controlling the accuracy of the numerical calculation of the Hessian matrix can influence the obtained vibrational frequencies to some extent. The effect of all these parameters is discussed and documented. It is concluded that with relatively economical computational conditions the uncertainty related to these parameters is smaller than 2-4 cm(-1). In the case of the Local Density Approximation scheme, comparison is possible with recent calculations performed with a Density Functional Perturbation Theory method and a plane-wave basis set. (C) 2004 Wiley Periodicals, Inc.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available