High Order Ab Initio Valence Force Field with Chemical Pattern Based Parameter Assignment.

Area of Science:

• Computational Chemistry
• Molecular Modeling
• Quantum Chemistry

Background:

• Molecular mechanics force fields (FFs) are crucial for simulating molecular structures and dynamics.
• Classical FFs often use simplified harmonic models for bonded interactions, neglecting anharmonicity and cross-coupling effects.
• These simplifications can limit the accuracy of vibrational properties and molecular configurations in simulations.

Purpose of the Study:

• To develop and validate improved bonded interaction parameters for molecular mechanics force fields.
• To incorporate vibrational anharmonicity and cross-coupling effects into valence models.
• To enhance the accuracy of molecular dynamics (MD) simulations by refining force field parameters.

Main Methods:

• Development of valence interaction parameters based on quantum mechanics (QM) calculations.
• Utilized approximately 270 atomic types defined by SMARTS strings.
• Employed an MM3(MM4)-style bonded interaction model, including out-of-plane bending for sp2-hybridized atoms.

Main Results:

• The developed valence parameters yield accurate vibrational frequencies with an RMSD of less than ~36.6 cm⁻¹ compared to QM.
• This represents a significant improvement over harmonic models, which typically show RMS errors greater than 60 cm⁻¹.
• The model accurately captures the potential energy surface for out-of-plane bending.

Conclusions:

• The new valence parameters offer enhanced accuracy for molecular mechanics force fields, particularly for vibrational properties.
• The model is applicable to the AMOEBA family and other MM3(MM4)-based FFs.
• This advancement can lead to more reliable molecular dynamics simulations across diverse chemical systems.