Computing pKa Values with a Mixing Hamiltonian Quantum Mechanical/Molecular Mechanical Approach

Area of Science:

• Computational Chemistry
• Physical Chemistry
• Chemical Physics

Background:

• Accurate computation of pKa values in solution is crucial but challenging.
• Existing methods may lack efficiency or comprehensive simulation capabilities for protonation/deprotonation processes.

Purpose of the Study:

• To develop a novel mixing quantum mechanical/molecular mechanical (QM/MM) Hamiltonian method for accurate pKa calculations.
• To enable efficient free-energy simulations of protonation/deprotonation in solution.

Main Methods:

• A new mixing QM/MM Hamiltonian method was developed, alchemically varying the nuclear charge of the transforming nucleus.
• The method simulates free-energy changes by varying the proton's charge fractionally from 0 to 1.
• Inspired by previous QM/MM free energy simulation methods and the linear combination of atomic potential approach.

Main Results:

• The pKa values of methanol and methanethiol in aqueous solution were calculated using the new method.
• The computational results demonstrated satisfactory agreement with experimental data.
• The method proves efficient for quantum mechanical free-energy simulations.

Conclusions:

• The developed mixing QM/MM Hamiltonian method provides an accurate and efficient approach for calculating pKa values.
• This method advances the simulation of protonation/deprotonation processes in solution.
• The technique offers a valuable tool for computational chemistry research.