Related Experiment Video
Updated: Dec 4, 2025

Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry
Published on: April 8, 2020
Analytical approach to phonon calculations in the SCC-DFTB framework.
Vladimir Bačić1, Thomas Heine2, Agnieszka Kuc1
1Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany.
A new analytical method accurately and efficiently calculates the Hessian for periodic systems using self-consistent-charge density-functional based tight-binding (SCC-DFTB). This approach outperforms traditional numerical methods for various materials.
Area of Science:
- Computational Materials Science
- Condensed Matter Physics
- Quantum Chemistry
Background:
- Calculating the Hessian is crucial for understanding vibrational properties and material stability.
- Traditional numerical force differentiation for Hessian calculation can be computationally expensive and less accurate for periodic systems.
- The self-consistent-charge density-functional based tight-binding (SCC-DFTB) method offers a computationally efficient approach to electronic structure calculations.
Purpose of the Study:
- To present a detailed derivation of an analytical, reciprocal-space approach for Hessian calculation within the SCC-DFTB framework.
- To establish an accurate and efficient method for obtaining the SCC-DFTB Hessian of periodic systems.
- To demonstrate the superiority of the analytical approach over traditional numerical methods.
Main Methods:
- Development of an analytical, reciprocal-space formalism for Hessian matrix computation.
- Implementation within the self-consistent-charge density-functional based tight-binding (SCC-DFTB) method.
- Validation against the traditional numerical force differentiation technique.
Main Results:
- The analytical, reciprocal-space approach provides accurate SCC-DFTB Hessian values for periodic systems.
- This method demonstrates significant efficiency gains compared to numerical differentiation.
- Superior performance was observed for diverse systems including doped graphene, graphene nanoribbons, boron-nitride nanotubes, and bulk zinc-oxide.
Conclusions:
- The presented analytical, reciprocal-space method is a highly accurate and efficient tool for SCC-DFTB Hessian calculations in periodic systems.
- This advancement offers a superior alternative to conventional numerical approaches for materials simulations.
- The method is broadly applicable to various periodic materials, facilitating advanced vibrational and stability analyses.
Related Concept Videos
IR Spectroscopy: Hooke's Law Approximation of Molecular Vibration
According to Hooke's law, the vibrational frequency is directly proportional to...
Molecular Orbital Theory I
Molecular Orbital Theory II
Crystal Field Theory - Octahedral Complexes
To explain the observed behavior of transition metal complexes (such as colors), a model involving electrostatic interactions between the electrons from the ligands and the electrons in the unhybridized d orbitals of the central metal atom has been developed. This electrostatic model is crystal field theory (CFT). It helps to understand, interpret, and predict the colors, magnetic behavior, and some structures of coordination compounds of transition metals.
CFT focuses on...
Van der Waals Equation
First, the attractive forces between molecules, which are stronger at higher densities and reduce the pressure, are considered by adding to the pressure a term equal to the square of the molar density multiplied by a positive coefficient a. Second, the volume...
Crystal Field Theory - Tetrahedral and Square Planar Complexes
Crystal field theory (CFT) is applicable to molecules in geometries other than octahedral. In octahedral complexes, the lobes of the dx2−y2 and dz2 orbitals point directly at the ligands. For tetrahedral complexes, the d orbitals remain in place, but with only four ligands located between the axes. None of the orbitals points directly at the tetrahedral ligands. However, the dx2−y2 and dz2 orbitals (along the Cartesian axes) overlap with the ligands less than the dxy,...

