计算化学公社

标题: 求助：计算溶剂化自由能在能量最小化步骤报错 [打印本页]

作者
Author: suxiao 时间: 2022-3-4 10:49
标题: 求助：计算溶剂化自由能在能量最小化步骤报错
我在计算喜树碱分子在离子液体中的溶剂化自由能，考虑了范德华力和库仑力的转变，但是多次更改参数文件，仍旧报相同的错误，请各位老师指导

报错信息如图所示：
Fatal error:
Function type Fourier Dih. not implemented in ip_pert

mdp文件参数如下：
integrator             = sd
dt                      = 0.002    ; 2 fs
nsteps                = 50000    ; 100 ps

nstenergy             = 200
nstlog                = 2000
nstxout-compressed    = 10000

gen-vel                = yes
gen-temp                = 298.15

constraint-algorithm    = lincs
constraints             = h-bonds

cutoff-scheme          = Verlet

coulombtype             = PME
rcoulomb                = 1.0

vdwtype                = Cut-off
rvdw                   = 1.0
DispCorr                = EnerPres

tcoupl                = Nose-Hoover
tc-grps                = System
tau-t                   = 2.0
ref-t                   = 298.15
nhchainlength          = 1

free-energy             = yes
init-lambda-state       = 0
calc-lambda-neighbors = -1
vdw-lambdas             = 0.00 0.20 0.40 0.60 0.80 1.00 1.00 1.00 1.00 1.00
coul-lambdas          = 0.00 0.00 0.00 0.00 0.00 0.00 0.25 0.50 0.75 1.00
couple-moltype          = CPT
couple-lambda0          = vdw-q
couple-lambda1          = none
couple-intramol       = no
nstdhdl                = 100
sc-alpha                = 0.5
sc-coul                = no
sc-power                = 1
sc-sigma                = 0.3

作者
Author: sobereva 时间: 2022-3-4 18:43
报错写了，此版本当前任务不支持这种扭转势函数的情况
尝试最新版本看看能否跑，要么换个力场

作者
Author: ccccc123 时间: 2024-11-14 20:29
同样用gmx2021遇到了这个问题，装了gmx2022.6就可以了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)