Analytic second nuclear derivatives of Hartree-Fock and DFT using multi-resolution analysis.

Area of Science:

• Computational chemistry
• Theoretical physics
• Molecular modeling

Background:

• Accurate computation of molecular properties is crucial for understanding chemical reactions and material properties.
• Existing methods for calculating second derivatives, essential for vibrational analysis, can be computationally intensive and face numerical challenges.

Purpose of the Study:

• To develop and implement a novel computational approach for calculating molecular second derivatives using multi-resolution analysis.
• To assess the accuracy and numerical stability of the proposed method for determining vibrational frequencies.

Main Methods:

• Developed a formalism for computing second derivatives of molecular systems within a multi-resolution analysis framework.
• Employed partial regularization of singular nuclear potentials to enhance numerical stability.
• Validated the method against large basis set Linear Combination of Atomic Orbitals (LCAO) calculations.

Main Results:

• Vibrational frequencies were reproduced with a Root Mean Square (RMS) deviation of a few cm^-1 compared to LCAO calculations.
• Identified that intermolecular modes, hindered rotations, and heavy atoms can impact precision.
• Demonstrated the necessity of tight precision thresholds for achieving numerically stable results.

Conclusions:

• The multi-resolution analysis approach provides a stable and accurate method for computing molecular second derivatives and vibrational frequencies.
• Careful selection of precision thresholds is essential for reliable results, particularly for complex molecular systems.
• This method offers a promising alternative for molecular dynamics and spectroscopic studies.