计算得到的吸附自由能比吸附能弱，请问是否正常？

a9471163

体系：A,B两个分子在水溶液中的稳定吸附体系
吸附能计算方法：gromacs中rerun定义两个能量组，分别是A，B分子，然后通过gmx energy命令，得到静电相互作用能和非静电相互作用能，总能量约为-100kJ/mol
吸附自由能计算方法：gromacs伞形采样后计算PMF曲线，得到横坐标是质心距离，纵坐标是能量的图，最后用PMF曲线最高点的y值，减去最低点的y值，计算得到的值是40kJ/mol，明显小于100kJ/mol
我个人觉得结果不正常，更换体系后（换了两个不同的分子，在水溶液中的吸附体系）再次计算，仍然是吸附自由能比吸附能弱，不清楚问题出在哪，还请老师指教，谢谢各位~
伞形采样时用到的umbrella.mdp参数如下
title       = Umbrella  simulation
define      = -DPOSRES_B
; Run parameters
integrator  = md
dt          = 0.002
tinit       = 0
nsteps      = 5000000   ; 10 ns
nstcomm     = 10
; Output parameters
nstxout     = 50000     ; every 100 ps
nstvout     = 50000
nstfout     = 5000
nstxtcout   = 5000      ; every 10 ps
nstenergy   = 5000
; Bond parameters
constraint_algorithm    = lincs
constraints             = all-bonds
continuation            = yes
; Single-range cutoff scheme
nstlist     = 5
ns_type     = grid
rlist       = 1.4
rcoulomb    = 1.4
rvdw        = 1.4
; PME electrostatics parameters
coulombtype     = PME
fourierspacing  = 0.12
fourier_nx      = 0
fourier_ny      = 0
fourier_nz      = 0
pme_order       = 4
ewald_rtol      = 1e-5
optimize_fft    = yes
; Berendsen temperature coupling is on in two groups
Tcoupl      = Nose-Hoover
tc_grps     = Other Water_and_ions
tau_t       = 0.5       0.5
ref_t       = 298.15       298.15
; Pressure coupling is on
Pcoupl          = Parrinello-Rahman
pcoupltype      = isotropic
tau_p           = 1.0      
compressibility = 4.5e-5
ref_p           = 1.0
refcoord_scaling = com
; Generate velocities is off
gen_vel     = no
; Periodic boundary conditions are on in all directions
pbc     = xyz
; Long-range dispersion correction
DispCorr    = EnerPres
; Pull code            
pull                    = yes      
pull_ngroups            = 2
pull_ncoords            = 1
pull_group1_name        = MOL
pull_group2_name        = SUB
pull_coord1_type        = umbrella      ; harmonic biasing force
pull_coord1_geometry    = distance      ; simple distance increase
pull_coord1_groups      = 1 2
pull_coord1_dim         = Y N N
pull_coord1_rate        = 0.0
pull_coord1_k           = 1000         
pull-coord1-start              = yes           ; define initial COM distance > 0

Tifo

k64_cc 发表于 2022-1-12 13:17
Variance大是一码事，和自由能曲线对不上是另一码事。

3 ns模拟，variance大是很正常的。一般来说，平 ...

感谢老师回复，我还有个问题想请教一下。就是我现在以一个离子为探针，把它从吸附稳定位置拉到体相位置，那等这个探针离子到体相位置的时候，其他离子难道不会重新吸附到之前探针离子所在的位点上吗，那这时候这个系统的总势能为什么还会和我拉那个探针之前的时候不一样？

k64_cc

Tifo 发表于 2022-1-10 21:10
老师您好，我也用和您说的一样的方法计算系统的总势能，就是在多个吸附离子中选择一个使用restraint将其 ...

Variance大是一码事，和自由能曲线对不上是另一码事。

3 ns模拟，variance大是很正常的。一般来说，平均势能的收敛需要更长的时间。我之前做水的PMF，已经很好弛豫了，还是需要几十ns的轨迹才能确保它收敛。如果收敛之后势能曲线还是和自由能对不上，那就说明是熵驱动的，不是焓驱动。这俩玩意也不一定非要对上。

Tifo

k64_cc 发表于 2020-11-5 11:58
不行，能量收敛比PMF慢，别想rerun了。而且G=H-TS的H是总势能，不是分解的，不用定义组，看总的Potential ...

老师您好，我也用和您说的一样的方法计算系统的总势能，就是在多个吸附离子中选择一个使用restraint将其控制在CV路径上的某一点3ns，再将获得的系统势能平均，但这样最后获得的势能随距离的变化曲线波动很大，与该离子的自由能变化对不上，请问是还有其他措施我没做吗？

k64_cc

a9471163 发表于 2020-11-6 15:36
我想用您说的这个方法计算熵变，请问有没有用过这种方法的参考文献可以引用

随便整一篇我老板从前的工作：

https://pubs.acs.org/doi/abs/10.1021/jacs.5b07232

a9471163

k64_cc 发表于 2020-11-5 11:58
不行，能量收敛比PMF慢，别想rerun了。而且G=H-TS的H是总势能，不是分解的，不用定义组，看总的Potential ...

我想用您说的这个方法计算熵变，请问有没有用过这种方法的参考文献可以引用

a9471163

k64_cc 发表于 2020-11-5 11:58
不行，能量收敛比PMF慢，别想rerun了。而且G=H-TS的H是总势能，不是分解的，不用定义组，看总的Potential ...

我明白了，您讲得很清楚，非常感谢老师

k64_cc

a9471163 发表于 2020-11-4 18:44
谢谢您的提醒~
请问PMF曲线上两个state之间的差值如何求出来，我目前想的办法是把state1和state2的分子 ...

不行，能量收敛比PMF慢，别想rerun了。而且G=H-TS的H是总势能，不是分解的，不用定义组，看总的Potential Energy就行。另外，水和溶质的相互作用能也会变化，中和一部分AB间作用能的变化，所以只分析这AB间相互作用也不是很充分。

饭得一口一口吃，你得先确定H主导了自由能变化，再分析什么因素主导了H变化。对于<H>的计算，你就在CV曲线上选两个点，用非常强的restraint把CV控制在那个点附近，然后跑三组长时间的模拟，算每组的平均势能，并用三组数据算variance。
Umbrella Sampling的原理是打格子数sample。你可以用一个模拟给多个格子提供数据。但是算PMF上某点平均势能的时候，你需要要尽量让sample处于你选定的那个格子里，所以要用一个很强的restraint来避免CV的波动。也是由于同个原因，US时的数据不能拿来直接用。

a9471163

k64_cc 发表于 2020-11-4 12:23
两个Energy Group之间的相互作用能不是dU，PMF曲线上两个state之间的差值才是能和dG对比分析的dU。能量分 ...

谢谢您的提醒~
请问PMF曲线上两个state之间<U>的差值如何求出来，我目前想的办法是把state1和state2的分子构型单独拿出来，用gromacs定义能量组rerun再求一遍静电能+非静电能，然后两个state的能量作差，不知道这种算法是不是合适啊

a9471163

sobereva 发表于 2020-11-4 14:34
看是对谁来说，对溶质的话，从结果上是如此
溶质之间的静电作用体现在结合的焓变的层面上，不是熵的层面 ...

明白了，谢谢老师

sobereva

a9471163 发表于 2020-11-4 09:20
假如物质疏水作用很强，物质表面的水会减少，所以能不能说，疏水作用强更有利于熵减小？

看是对谁来说，对溶质的话，从结果上是如此
溶质之间的静电作用体现在结合的焓变的层面上，不是熵的层面上

k64_cc

a9471163 发表于 2020-11-3 09:20
你好，我是这样理解这个问题的：
ΔG=ΔU+PΔV-TΔS，目前ΔG=-40，ΔU=-100，PΔV项可以忽略不计，因为 ...

两个Energy Group之间的相互作用能不是dU，PMF曲线上两个state之间<U>的差值才是能和dG对比分析的dU。能量分解得到的相互作用能可能对自由能变化起主导作用，但是不能拿来和自由能曲线一起做定量分析。

a9471163

sobereva 发表于 2020-11-4 04:17
对于溶质来说，熵减小是普遍现象。吸附前单体可以自由运动，吸附后运动受限，因此熵减小才是常见的。
虽 ...

假如物质疏水作用很强，物质表面的水会减少，所以能不能说，疏水作用强更有利于熵减小？

a9471163

sobereva 发表于 2020-11-4 04:17
对于溶质来说，熵减小是普遍现象。吸附前单体可以自由运动，吸附后运动受限，因此熵减小才是常见的。
虽 ...

谢谢sob老师，请问静电作用一般如何对熵作出贡献？
我的理解是，假如两个物质之间因为有疏水作用、静电作用而吸附在一起了，物质A和B表面的水化膜会解缚成水分子回到溶液中，这部分熵是增加的，而两个分子的吸附是熵减小的，所以也很难说疏水作用、静电作用对熵的总体贡献是增大熵还是减小熵，不知是否正确？

a9471163

wzkchem5 发表于 2020-11-3 10:59
多的是，没准比熵增大的例子还多。
给个最简单的例子，DNA单链结合成双链，除了端基的AT碱基对以外，所 ...

非常感谢老师，我明白了