Analytic cubic and quartic force fields using density-functional theory
- Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Tromsø—The Arctic University of Norway, 9037 Tromsø (Norway)
- Theoretical Chemistry and Biology, School of Biotechnology, Royal Institute of Technology, AlbaNova University Center, S-10691 Stockholm, Sweden and PDC Center for High Performance Computing, Royal Institute of Technology, S-10044 Stockholm (Sweden)
- Center for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo (Norway)
We present the first analytic implementation of cubic and quartic force constants at the level of Kohn–Sham density-functional theory. The implementation is based on an open-ended formalism for the evaluation of energy derivatives in an atomic-orbital basis. The implementation relies on the availability of open-ended codes for evaluation of one- and two-electron integrals differentiated with respect to nuclear displacements as well as automatic differentiation of the exchange–correlation kernels. We use generalized second-order vibrational perturbation theory to calculate the fundamental frequencies of methane, ethane, benzene, and aniline, comparing B3LYP, BLYP, and Hartree–Fock results. The Hartree–Fock anharmonic corrections agree well with the B3LYP corrections when calculated at the B3LYP geometry and from B3LYP normal coordinates, suggesting that the inclusion of electron correlation is not essential for the reliable calculation of cubic and quartic force constants.
- OSTI ID:
- 22255276
- Journal Information:
- Journal of Chemical Physics, Vol. 140, Issue 3; Other Information: (c) 2014 Author(s); Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
A New, Self-Contained Asymptotic Correction Scheme To Exchange-Correlation Potentials For Time-Dependent Density Functional Theory
Evaluating interaction energies of weakly bonded systems using the Buckingham-Hirshfeld method