Assessing the Partial Hessian Approximation in QM/MM-Based Vibrational Analysis.

Area of Science:

• Computational Chemistry
• Molecular Mechanics
• Quantum Mechanics

Background:

• Calculating the full Hessian matrix for vibrational analysis in QM/MM systems is computationally expensive.
• The partial Hessian approximation focuses on the QM subsystem, aligning with QM/MM's core principles.
• This approximation is commonly used to study local vibrational modes within the QM subsystem.

Purpose of the Study:

• To evaluate the accuracy and applicability of the partial Hessian vibrational analysis (PHVA) as typically implemented in QM/MM.
• To differentiate and quantify errors arising from the partial Hessian approximation versus QM/MM embedding.
• To investigate the treatment of pseudotranslational and pseudorotational modes in PHVA.

Main Methods:

• Compared PHVA with full Hessian calculations at the QM level for normal modes, frequencies, and IR/Raman intensities.
• Quantified QM/MM embedding errors by comparing PHVA results from QM/MM-type embedding with QM-level PHVA.
• Analyzed and proposed methods for handling pseudotranslational and pseudorotational modes.

Main Results:

• The partial Hessian approximation shows good accuracy for local vibrational modes in QM/MM systems.
• QM/MM embedding introduces errors that need to be quantified separately from the partial Hessian approximation errors.
• Removing pseudotranslational and pseudorotational modes by projection can negatively impact other normal modes; identification and removal are preferred.

Conclusions:

• PHVA is a viable approach for vibrational analysis in QM/MM systems, particularly for local modes.
• Careful consideration of QM/MM embedding errors and proper handling of spurious modes are crucial for accurate results.
• The study provides insights into optimizing PHVA for reliable vibrational analysis in complex molecular systems.