These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.
Preprints are manuscripts made publicly available before they have been submitted for formal peer review and publication. They might contain new research findings or data. Preprints can be a draft or final version of an author's research but must not have been accepted for publication at the time of submission.
revised on 01.07.2020 and posted on 02.07.2020by Tomas Bucko, Monika Gešvandtnerová, Dario Rocca
While free energies are fundamental thermodynamic quantities to characterize chemical reactions, their calculation based on ab initio theory is usually limited by the high computational cost. This is particularly true if multiple levels of theory have to be tested to establish their relative accuracy, if highly expensive quantum mechanical approximations are of interest, and also if several different temperatures have to be considered. We present an ab initio approach that effectively couples perturbation theory and machine learning to make ab initio free energy calculations more affordable. Starting from results based on a certain production ab initio theory, perturbation theory is applied to obtain free energies. The large number of single point calculations required by a brute force application of this approach are here significantly decreased by applying machine learning techniques. Importantly, the training of the machine learning model requires only a small amount of data and does not need to be performed again when the temperature is decreased. The accuracy and efficiency of this method is demonstrated by computing the free energy of activation of the proton exchange reaction in the zeolite chabazite. Starting from an ab initio calculation based on a semilocal approximation of density functional theory, free energies based on significantly more expensive non-local van der Waals and hybrid functionals are obtained with only a few tens of additional single point calculations. In this way this work paves the route to quick free energy calculations using different levels of theory or approximations that would be too computationally expensive to be directly employed in molecular dynamics or Monte Carlo simulations.