JCenter 发表于 2023-7-10 18:39 没法跟相互作用能或者其中任何的成份比较，本身疏水作用就不体现在力场项里面 |
sobereva 发表于 2021-10-17 17:55 老师，您好。我的目标是对水中的两个分子（A和B）间计算弱相互作用和疏水效应。首先，关于弱相互作用，我是通过molcus+genmer搜索和优化AB团簇构型，对能量对低的团簇构型的波函数文件通过基于Multiwfn的EDA-FF方法可以得到总的弱相互作用能，以及分解出来的静电相互作用、色散相互作用的值（虽然目前已经有了sobEDA方法，但我之前已经通过EDA-FF算完了，本质是通过GAFF分子立场下对其进行分析的）。对于疏水效应，我想与EDA-FF方法分解出来的静电作用、色散进行横向对比分析。看了大量帖子后，有些感想和疑问，想请老师解疑答惑： 1.通过gmx模拟去掉非键相互作用，计算逐渐远离过程的PMF曲线，用PMF曲线收敛到平坦部分时的值减去最小点数值就可以体现疏水作用能，这种方式算出来的值可以和EDA-FF方法得出来的值进行对比吗，如果不能的话，用什么其他方法可以与EDA-FF方法得出来的静电、色散比较呢。 2.在gmx模拟做PMF曲线时。每种分子的原子电荷是否需要设置为0，不知道老师有无关于做PMF曲线的链接或者幻灯片呢。 3.关于另外一种疏水效应比较的方法，溶解自由能中的非极性部分，但并非完全是疏水效应带给非极性的贡献，这种计算感觉可以侧面仅仅比较疏水效应的强弱，但由于非极性部分包含的较多，如疏水效应、孔洞能等，应该不能与EDA-FF分解出来的各项作用进行讨论分析，不知道我是否理解的对。 4. 关于相互作用能，eda-ff方法和通过gmx模拟得出的静电作用能和范德华作用能，在精度上，这两种哪种更好呢，(sobEDA方法精度比这两种高） |
a9471163 发表于 2021-10-17 11:19 不同 这种方式得到的是两个分子间的疏水作用能，和某个分子的溶解自由能里的疏水效应对溶解自由能的贡献不是一码事 很少有专门算孔洞能的。Gaussian里支持一种算孔洞能的方式，原文是A Scaled Particle Theory of Aqueous and Nonaqueous Solutions，但我没仔细看，也不是什么流行的做法。 |
sobereva 发表于 2021-10-17 06:05 感谢社长，能否理解为“用PMF曲线收敛到平坦部分时的值减去最小点数值得到的疏水作用能”在数值上等于“疏水效应对溶解自由能的贡献”；有没有办法单独求出孔洞能 |
|
1、2 用PMF曲线收敛到平坦部分时的值减去最小点数值就可以体现疏水作用能 3 -odg那个给出的只不过是把SASA乘上个自定义的系数，如果靠这个得到溶解自由能的非极性部分，不能说这就是疏水效应对溶解自由能的贡献，本来非极性部分还有溶剂-溶质之间的范德华作用、孔洞能等物理成分 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-24 05:50 , Processed in 0.168951 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites