Fowles Daniel J, Palmer David S
Department of Pure and Applied Chemistry, University of Strathclyde, Thomas Graham Building, 295 Cathedral Street, Glasgow, Scotland G1 1XL, UK.
Phys Chem Chem Phys. 2023 Mar 1;25(9):6944-6954. doi: 10.1039/d3cp00199g.
Simultaneous calculation of entropies, enthalpies and free energies has been a long-standing challenge in computational chemistry, partly because of the difficulty in obtaining estimates of all three properties from a single consistent simulation methodology. This has been particularly true for methods from the Integral Equation Theory of Molecular Liquids such as the Reference Interaction Site Model which have traditionally given large errors in solvation thermodynamics. Recently, we presented pyRISM-CNN, a combination of the 1 Dimensional Reference Interaction Site Model (1D-RISM) solver, pyRISM, with a deep learning based free energy functional, as a method of predicting solvation free energy (SFE). With this approach, a 40-fold improvement in prediction accuracy was delivered for a multi-solvent, multi-temperature dataset when compared to the standard 1D-RISM theory [Fowles , 2023, , 177-188]. Here, we report three further developments to the pyRISM-CNN methodology. Firstly, solvation free energies have been introduced for organic molecular ions in methanol or water solvent systems at 298 K, with errors below 4 kcal mol obtained without the need for corrections or additional descriptors. Secondly, the number of solvents in the training data has been expanded from carbon tetrachloride, water and chloroform to now also include methanol. For neutral solutes, prediction errors nearing or below 1 kcal mol are obtained for each organic solvent system at 298 K and water solvent systems at 273-373 K. Lastly, pyRISM-CNN was successfully applied to the simultaneous prediction of solvation enthalpy, entropy and free energy through a multi-task learning approach, with errors of 1.04, 0.98 and 0.47 kcal mol, respectively, for water solvent systems at 298 K.
同时计算熵、焓和自由能一直是计算化学领域长期存在的挑战,部分原因在于难以从单一一致的模拟方法中获得这三种性质的估计值。对于分子液体积分方程理论中的方法,如参考相互作用位点模型,情况尤其如此,该模型在溶剂化热力学中传统上会产生较大误差。最近,我们提出了pyRISM-CNN,它将一维参考相互作用位点模型(1D-RISM)求解器pyRISM与基于深度学习的自由能泛函相结合,作为预测溶剂化自由能(SFE)的一种方法。与标准的1D-RISM理论相比,这种方法在多溶剂、多温度数据集上的预测精度提高了40倍[福尔斯,2023,,177 - 188]。在此,我们报告了pyRISM-CNN方法的三个进一步发展。首先,引入了298 K下甲醇或水溶剂体系中有机分子离子的溶剂化自由能,在无需校正或额外描述符的情况下,误差低于4 kcal/mol。其次,训练数据中的溶剂数量已从四氯化碳、水和氯仿扩展到现在还包括甲醇。对于中性溶质,在298 K时每个有机溶剂体系以及在273 - 373 K时水溶剂体系的预测误差接近或低于1 kcal/mol。最后,pyRISM-CNN通过多任务学习方法成功应用于溶剂化焓、熵和自由能的同时预测,对于298 K的水溶剂体系,误差分别为1.04、0.98和0.47 kcal/mol。