计算化学公社

 找回密码 Forget password
 注册 Register
Views: 18662|回复 Reply: 37
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 使用现成的甲烷水合物跑MD,在能量最小化这一步,原子名就对不上,无法继续

[复制链接 Copy URL]

1478

帖子

0

威望

4537

eV
积分
6015

Level 6 (一方通行)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 牧生 于 2021-8-18 22:13 编辑

一、尝试做甲烷水合物的模拟,使用CrystalMaker直接生成一个甲烷水合物的pdb晶胞,


二、使用ligpargen在线得到甲烷的itp文件,把[ atomtypes ]内容剪切到ffnonbonded.itp,用opls力场,把forcefield.itp,ffnonbonded.itp等相关文件都复制到当前目录下
三、根据碳原子和氧原子的数量,手搓top文件为
#include "forcefield.itp"
;#include "tip4p-ice.itp"
#include "spce.itp"                ;这里用的三点水模型,计划在能量最小化以后,再转为四点水模型
#include "jiawan.itp"            ;其次才出现甲烷的itp
[ system ]
jiawan in water
[ molecules ]
; Compound        nmols
SOL           2484    ;先出现水
jiawan      432      ;再出现甲烷


四、对这个水合物晶胞进行能量最小化,就报错:

Setting the LD random seed to 1566708099
Generated 337431 of the 337431 non-bonded parameter combinations
Generating 1-4 interactions: fudge = 0.5
Generated 337431 of the 337431 1-4 parameter combinations
Excluding 2 bonded neighbours molecule type 'SOL'
Excluding 3 bonded neighbours molecule type 'jiawan'
Warning: atom name 1 in topol.top and 333.gro does not match (OW - C)
Warning: atom name 2 in topol.top and 333.gro does not match (HW1 - C)
Warning: atom name 3 in topol.top and 333.gro does not match (HW2 - C)
Warning: atom name 4 in topol.top and 333.gro does not match (OW - C)
Warning: atom name 5 in topol.top and 333.gro does not match (HW1 - C)
Warning: atom name 6 in topol.top and 333.gro does not match (HW2 - C)
Warning: atom name 7 in topol.top and 333.gro does not match (OW - C)
Warning: atom name 8 in topol.top and 333.gro does not match (HW1 - C)
Warning: atom name 9 in topol.top and 333.gro does not match (HW2 - C)
Warning: atom name 10 in topol.top and 333.gro does not match (OW - C)
Warning: atom name 11 in topol.top and 333.gro does not match (HW1 - C)
Warning: atom name 12 in topol.top and 333.gro does not match (HW2 - C)
Warning: atom name 13 in topol.top and 333.gro does not match (OW - C)
Warning: atom name 14 in topol.top and 333.gro does not match (HW1 - C)
Warning: atom name 15 in topol.top and 333.gro does not match (HW2 - C)
Warning: atom name 16 in topol.top and 333.gro does not match (OW - C)
Warning: atom name 17 in topol.top and 333.gro does not match (HW1 - C)
Warning: atom name 18 in topol.top and 333.gro does not match (HW2 - C)
Warning: atom name 19 in topol.top and 333.gro does not match (OW - C)
Warning: atom name 20 in topol.top and 333.gro does not match (HW1 - C)
(more than 20 non-matching atom names)

WARNING 1 [file topol.top, line 17]:
  9612 non-matching atom names
  atom names from topol.top will be used
  atom names from 333.gro will be ignored

Removing all charge groups because cutoff-scheme=Verlet
Analysing residue names:
There are:  2484      Water residues
There are:   432      Other residues
Analysing residues not classified as Protein/DNA/RNA/Water and splitting into groups...
Number of degrees of freedom in T-Coupling group rest is 28833.00
Calculating fourier grid dimensions for X Y Z
Using a fourier grid of 60x32x32, spacing 0.118 0.111 0.111
Estimate for the relative computational load of the PME mesh part: 0.18
This run will generate roughly 57 Mb of data

There was 1 warning

-------------------------------------------------------
Program:     gmx grompp, version 2019.6
Source file: src\gromacs\gmxpreprocess\grompp.cpp (line 2315)

Fatal error:
Too many warnings (1).
If you are sure all warnings are harmless, use the -maxwarn option.
但是看起来,top中,甲烷和水的顺序是对了的,为何原子名也是对不上的?




五、如果,我把top里面甲烷和水的位置换一下,会有如下报错(当然,我知道这个是因为顺序不对而报错)。
Setting the LD random seed to 368889344
Generated 337431 of the 337431 non-bonded parameter combinations
Generating 1-4 interactions: fudge = 0.5
Generated 337431 of the 337431 1-4 parameter combinations
Excluding 3 bonded neighbours molecule type 'jiawan'
Excluding 2 bonded neighbours molecule type 'SOL'
Warning: atom name 1 in topol.top and 333.gro does not match (C00 - C)
Warning: atom name 2 in topol.top and 333.gro does not match (H01 - C)
Warning: atom name 3 in topol.top and 333.gro does not match (H02 - C)
Warning: atom name 4 in topol.top and 333.gro does not match (H03 - C)
Warning: atom name 5 in topol.top and 333.gro does not match (H04 - C)
Warning: atom name 6 in topol.top and 333.gro does not match (C00 - C)
Warning: atom name 7 in topol.top and 333.gro does not match (H01 - C)
Warning: atom name 8 in topol.top and 333.gro does not match (H02 - C)
Warning: atom name 9 in topol.top and 333.gro does not match (H03 - C)
Warning: atom name 10 in topol.top and 333.gro does not match (H04 - C)
Warning: atom name 11 in topol.top and 333.gro does not match (C00 - C)
Warning: atom name 12 in topol.top and 333.gro does not match (H01 - C)
Warning: atom name 13 in topol.top and 333.gro does not match (H02 - C)
Warning: atom name 14 in topol.top and 333.gro does not match (H03 - C)
Warning: atom name 15 in topol.top and 333.gro does not match (H04 - C)
Warning: atom name 16 in topol.top and 333.gro does not match (C00 - C)
Warning: atom name 17 in topol.top and 333.gro does not match (H01 - C)
Warning: atom name 18 in topol.top and 333.gro does not match (H02 - C)
Warning: atom name 19 in topol.top and 333.gro does not match (H03 - C)
Warning: atom name 20 in topol.top and 333.gro does not match (H04 - C)
(more than 20 non-matching atom names)

WARNING 1 [file topol.top, line 16]:
  9612 non-matching atom names
  atom names from topol.top will be used
  atom names from 333.gro will be ignored


请大佬帮忙看看,为何这个pdb晶胞,在进行能量最小化这一步,原子名都对不上呢。如果强行忽略这个报错,那么能量最小化以后,晶体结构就散了。

附上所用的文件

em.mdp (388 Bytes, 下载次数 Times of downloads: 50)
jiawan.itp (1.36 KB, 下载次数 Times of downloads: 78)
topol.top (200 Bytes, 下载次数 Times of downloads: 60) ffnonbonded.itp (62.01 KB, 下载次数 Times of downloads: 39)
ffbonded.itp (204.7 KB, 下载次数 Times of downloads: 31)
forcefield.itp (862 Bytes, 下载次数 Times of downloads: 30)


pdb文件略大,传到了蓝奏云,将后缀改为pdb即可https://wwa.lanzoui.com/i9fATstho8b







又菜又爱玩

308

帖子

2

威望

3557

eV
积分
3905

Level 5 (御坂)

2#
发表于 Post on 2021-8-18 23:27:04 | 只看该作者 Only view this author
pdb中甲烷分子,水分子内部原子顺序都不是严格按照C-H-H-H-H和O-H-H的顺序出现的,能模拟成功就见鬼了。

1478

帖子

0

威望

4537

eV
积分
6015

Level 6 (一方通行)

3#
 楼主 Author| 发表于 Post on 2021-8-19 08:11:01 | 只看该作者 Only view this author
lyj714 发表于 2021-8-18 23:27
pdb中甲烷分子,水分子内部原子顺序都不是严格按照C-H-H-H-H和O-H-H的顺序出现的,能模拟成功就见鬼了。

请大佬能不能授人以渔啊,如何得到一个排序正确的甲烷水合物pdb
又菜又爱玩

308

帖子

2

威望

3557

eV
积分
3905

Level 5 (御坂)

4#
发表于 Post on 2021-8-21 11:54:56 | 只看该作者 Only view this author
牧生 发表于 2021-8-19 08:11
请大佬能不能授人以渔啊,如何得到一个排序正确的甲烷水合物pdb

我的建议是自己编写程序,对每一个C和O氧原子周围的H做一个遍历循环,计算找出C-H和O-H距离符合实际键长的H,按照C-H-H-H和O-H-H的顺序输出坐标即可。这就是授人以渔,但实际上如果你连最基本的编程都不会,那我的渔对你一点用处都没有,所以鱼比渔要更直接。

1478

帖子

0

威望

4537

eV
积分
6015

Level 6 (一方通行)

5#
 楼主 Author| 发表于 Post on 2021-8-23 11:20:40 | 只看该作者 Only view this author
lyj714 发表于 2021-8-21 11:54
我的建议是自己编写程序,对每一个C和O氧原子周围的H做一个遍历循环,计算找出C-H和O-H距离符合实际键长 ...

感谢大佬对我的帮助。
sob大佬已经给予了鱼和渔了。但我仍然遇到了一个看似简单,我却一直无法解决的问题。请帮忙看一下问题出在哪里。

分别都试过使用amber力场,opls力场,都是一样的情况。

第一步:CrystalMaker建立了甲烷水合物晶胞,用Multiwfn理顺了pdb中的原子排序。确保了pdb中原子顺序严格按照C-H-H-H-H和O-H-H出现,晶胞大小为2.4×2.4×2.4

第二步:使用packmol建立了甲烷水合物和水的一个盒子,大小为2.4×2.4×4.0

第三步:能量最小以后,将得到的gro文件转为四点水模型。修改四点水模型后的gro文件中原子总数,并理顺gro中的原子编号。

第四步:低温下(260K)进行NPT模拟,可以看到水合物的部分基本保持原状,水的部分有缓慢结冰的迹象。

图形如下:    


问题出在下面的操作中,一旦用packmol加入的是甲烷和水分子,进行操作就出错了,具体如下:

第一步和之前一样:CrystalMaker建立了甲烷水合物晶胞,用Multiwfn理顺了pdb中的原子排序。确保了pdb中原子顺序严格按照C-H-H-H-H和O-H-H出现,晶胞大小为2.4×2.4×2.4

第二步,仅仅多加入了几个甲烷:使用packmol建立了甲烷水合物和水和10个甲烷分子的一个盒子,大小为2.4×2.4×4.0      (也尝试过只加一个甲烷,也会有问题。也试过建立层状的水合物——水——甲烷的体系,也有问题)

第三步:能量最小这一步就错了。虽然表面上看来进行完成了,但实际上,vmd打开得到的gro文件就有问题,
Warning) Unusual bond between residues:  2 (none) and 391 (waters)
Warning) Unusual bond between residues:  8 (none) and 389 (waters)
Warning) Unusual bond between residues:  11 (none) and 348 (waters)
Warning) Unusual bond between residues:  12 (none) and 229 (waters)
Warning) Unusual bond between residues:  15 (none) and 352 (waters)
Warning) Unusual bond between residues:  18 (none) and 203 (waters)
Warning) Unusual bond between residues:  18 (none) and 203 (waters)
Warning) Unusual bond between residues:  20 (none) and 261 (waters)
Warning) Unusual bond between residues:  22 (none) and 406 (waters)
Warning) Unusual bond between residues:  24 (none) and 209 (waters)
Warning) Unusual bond between residues:  24 (none) and 265 (waters)
Warning) Unusual bond between residues:  24 (none) and 209 (waters)
Warning) Unusual bond between residues:  30 (none) and 126 (waters)
Warning) Unusual bond between residues:  30 (none) and 226 (waters)
Warning) Unusual bond between residues:  35 (none) and 340 (waters)
Warning) Unusual bond between residues:  39 (none) and 344 (waters)
Warning) Unusual bond between residues:  40 (none) and 273 (waters)
Warning) Unusual bond between residues:  43 (none) and 298 (waters)
Warning) Unusual bond between residues:  43 (none) and 298 (waters)
Warning) Unusual bond between residues:  46 (none) and 367 (waters)
Warning) Unusual bond between residues:  46 (none) and 367 (waters)
Warning) Unusual bond between residues:  64 (none) and 270 (waters)
Warning) Unusual bond between residues:  66 (none) and 130 (waters)
Warning) Unusual bond between residues:  66 (none) and 166 (waters)
Warning) Unusual bond between residues:  66 (none) and 387 (waters)
Warning) Unusual bond between residues:  66 (none) and 387 (waters)
Warning) Unusual bond between residues:  67 (none) and 178 (waters)
Warning) Unusual bond between residues:  67 (none) and 225 (waters)
Warning) Unusual bond between residues:  67 (none) and 368 (waters)
Warning) Unusual bond between residues:  67 (none) and 384 (waters)
Warning) Unusual bond between residues:  68 (none) and 158 (waters)
Warning) Unusual bond between residues:  71 (none) and 397 (waters)
Warning) Unusual bond between residues:  71 (none) and 268 (waters)
Warning) Unusual bond between residues:  73 (none) and 171 (waters)
Info)    Segments: 1
Info)    Fragments: 1402   Protein: 0   Nucleic: 0
The alascan package could not be loaded:

图形如下:


使用的文件如下:


222.pdb (111.37 KB, 下载次数 Times of downloads: 39)
em.mdp (387 Bytes, 下载次数 Times of downloads: 24)
ffbonded.itp (35.13 KB, 下载次数 Times of downloads: 23)
ffnonbonded.itp (5.04 KB, 下载次数 Times of downloads: 20)
forcefield.itp (929 Bytes, 下载次数 Times of downloads: 21)
mix.gro (120.17 KB, 下载次数 Times of downloads: 51)
topol.top (202 Bytes, 下载次数 Times of downloads: 51)

运行以下命令后,盒子中的水分子就乱了,使用的em.mdp文件都是同一个
gmx grompp -f em.mdp -c mix.gro -p topol.top -o em-sol.tpr -maxwarn 80
gmx mdrun -deffnm em-sol -v


请教一下怎么办
又菜又爱玩

308

帖子

2

威望

3557

eV
积分
3905

Level 5 (御坂)

6#
发表于 Post on 2021-8-23 14:01:00 | 只看该作者 Only view this author
牧生 发表于 2021-8-23 11:20
感谢大佬对我的帮助。
sob大佬已经给予了鱼和渔了。但我仍然遇到了一个看似简单,我却一直无法解决的问 ...

top中的分子不能合并这一点我就先不说了,你这mix.gro中的甲烷水合物中的H2O顺序不对啊,323行你开始看,明显的H-H-O的顺序。

1478

帖子

0

威望

4537

eV
积分
6015

Level 6 (一方通行)

7#
 楼主 Author| 发表于 Post on 2021-8-23 17:32:52 | 只看该作者 Only view this author
lyj714 发表于 2021-8-23 14:01
top中的分子不能合并这一点我就先不说了,你这mix.gro中的甲烷水合物中的H2O顺序不对啊,323行你开始看, ...

感谢大佬的帮助。。

我重新调了顺序,使用packmol建立了甲烷水合物+200水分子+10个甲烷分子的盒子,大小2.4×2.4×4,并检查看,甲烷全部是C-H-H-H-H,水分子都是H-H-O    (这里水分子用了培训班的pdb文件)

但是能量最小化,还是一样的提示,

Warning: atom name 1 in topol.top and mix.gro does not match (C00 - C)
Warning: atom name 2 in topol.top and mix.gro does not match (H01 - H)
Warning: atom name 3 in topol.top and mix.gro does not match (H02 - H)
Warning: atom name 4 in topol.top and mix.gro does not match (H03 - H)
Warning: atom name 5 in topol.top and mix.gro does not match (H04 - H)
Warning: atom name 6 in topol.top and mix.gro does not match (C00 - C)
Warning: atom name 7 in topol.top and mix.gro does not match (H01 - H)
Warning: atom name 8 in topol.top and mix.gro does not match (H02 - H)
Warning: atom name 9 in topol.top and mix.gro does not match (H03 - H)
Warning: atom name 10 in topol.top and mix.gro does not match (H04 - H)
Warning: atom name 11 in topol.top and mix.gro does not match (C00 - C)
Warning: atom name 12 in topol.top and mix.gro does not match (H01 - H)
Warning: atom name 13 in topol.top and mix.gro does not match (H02 - H)
Warning: atom name 14 in topol.top and mix.gro does not match (H03 - H)
Warning: atom name 15 in topol.top and mix.gro does not match (H04 - H)
Warning: atom name 16 in topol.top and mix.gro does not match (C00 - C)
Warning: atom name 17 in topol.top and mix.gro does not match (H01 - H)
Warning: atom name 18 in topol.top and mix.gro does not match (H02 - H)
Warning: atom name 19 in topol.top and mix.gro does not match (H03 - H)
Warning: atom name 20 in topol.top and mix.gro does not match (H04 - H)
(more than 20 non-matching atom names)



这些错误只是原子名对不上而已,我觉得是可以忽略的。-maxwarn 忽略以后,得到了gro文件,VMD打开这个gro时,有提示

Warning) Unusual bond between residues:  2 (none) and 315 (waters)
Warning) Unusual bond between residues:  3 (none) and 579 (waters)
Warning) Unusual bond between residues:  4 (none) and 109 (waters)
Warning) Unusual bond between residues:  4 (none) and 109 (waters)
Warning) Unusual bond between residues:  5 (none) and 275 (waters)
Warning) Unusual bond between residues:  7 (none) and 126 (waters)
Warning) Unusual bond between residues:  8 (none) and 362 (waters)
Warning) Unusual bond between residues:  9 (none) and 302 (waters)
Warning) Unusual bond between residues:  10 (none) and 388 (waters)
Warning) Unusual bond between residues:  10 (none) and 227 (waters)
Warning) Unusual bond between residues:  11 (none) and 165 (waters)
Warning) Unusual bond between residues:  12 (none) and 94 (waters)
Warning) Unusual bond between residues:  12 (none) and 294 (waters)
Warning) Unusual bond between residues:  13 (none) and 360 (waters)
Warning) Unusual bond between residues:  13 (none) and 360 (waters)
Warning) Unusual bond between residues:  13 (none) and 375 (waters)
Warning) Unusual bond between residues:  14 (none) and 426 (waters)
Warning) Unusual bond between residues:  14 (none) and 426 (waters)
Warning) Unusual bond between residues:  14 (none) and 426 (waters)
Warning) Unusual bond between residues:  14 (none) and 86 (waters)
Warning) Unusual bond between residues:  15 (none) and 485 (waters)
Warning) Unusual bond between residues:  15 (none) and 485 (waters)
Warning) Unusual bond between residues:  19 (none) and 169 (waters)
Warning) Unusual bond between residues:  20 (none) and 253 (waters)
Warning) Unusual bond between residues:  20 (none) and 165 (waters)
Warning) Unusual bond between residues:  20 (none) and 213 (waters)
Warning) Unusual bond between residues:  21 (none) and 294 (waters)
Warning) Unusual bond between residues:  28 (none) and 530 (waters)
Warning) Unusual bond between residues:  29 (none) and 362 (waters)
Warning) Unusual bond between residues:  29 (none) and 362 (waters)
Warning) Unusual bond between residues:  34 (none) and 267 (waters)
Warning) Unusual bond between residues:  34 (none) and 267 (waters)
Warning) Unusual bond between residues:  38 (none) and 256 (waters)
Warning) Unusual bond between residues:  39 (none) and 382 (waters)
Warning) Unusual bond between residues:  40 (none) and 141 (waters)
Warning) Unusual bond between residues:  44 (none) and 312 (waters)
Warning) Unusual bond between residues:  45 (none) and 77 (waters)
Warning) Unusual bond between residues:  46 (none) and 254 (waters)
Warning) Unusual bond between residues:  47 (none) and 170 (waters)
Warning) Unusual bond between residues:  47 (none) and 384 (waters)
Warning) Unusual bond between residues:  47 (none) and 238 (waters)
Warning) Unusual bond between residues:  48 (none) and 342 (waters)
Warning) Unusual bond between residues:  48 (none) and 348 (waters)
Warning) Unusual bond between residues:  48 (none) and 394 (waters)
Warning) Unusual bond between residues:  48 (none) and 346 (waters)
Warning) Unusual bond between residues:  54 (none) and 416 (waters)
Warning) Unusual bond between residues:  57 (none) and 412 (waters)
Warning) Unusual bond between residues:  60 (none) and 197 (waters)
Warning) Unusual bond between residues:  61 (none) and 82 (waters)
Warning) Unusual bond between residues:  62 (none) and 199 (waters)
Warning) Unusual bond between residues:  63 (none) and 327 (waters)
Warning) Unusual bond between residues:  65 (none) and 197 (waters)
Warning) Unusual bond between residues:  65 (none) and 282 (waters)
Warning) Unusual bond between residues:  65 (none) and 293 (waters)
Warning) Unusual bond between residues:  65 (none) and 405 (waters)
Warning) Unusual bond between residues:  67 (none) and 127 (waters)
Warning) Unusual bond between residues:  72 (none) and 383 (waters)
Warning) Unusual bond between residues:  72 (none) and 126 (waters)
Warning) Unusual bond between residues:  72 (none) and 383 (waters)
Warning) Unusual bond between residues:  72 (none) and 398 (waters)



既然顺序对了,那么可能的错误就是您提到的top中的分子不能合并这一点我就先不说了”,请能否再次帮我一下,如何在这里写一个正确的top。

topol.top (201 Bytes, 下载次数 Times of downloads: 45)

mix.gro (93.23 KB, 下载次数 Times of downloads: 43)


又菜又爱玩

308

帖子

2

威望

3557

eV
积分
3905

Level 5 (御坂)

8#
发表于 Post on 2021-8-23 19:11:16 | 只看该作者 Only view this author
本帖最后由 lyj714 于 2021-8-23 19:16 编辑

后面回复

308

帖子

2

威望

3557

eV
积分
3905

Level 5 (御坂)

9#
发表于 Post on 2021-8-23 19:15:57 | 只看该作者 Only view this author
lyj714 发表于 2021-8-23 19:11
哪里来的自信说顺序对了,H-H-O你认为个spc.itp中的[ atoms ]字段出现的顺序一样???你仔细看看吧。我 ...

哪里来的自信说顺序对了,你这个gro中的是H-H-O,但是你spc.itp中的atoms字段可是OHH的顺序啊,除非你把spc.itp中所有东西,包括原子顺序,成键原子弄对,否则就错误的。另外我说的不能合并指的是top你的两部分甲烷不能合并到一起,分子字段这一部分展开以后必须和你的gro匹配才行,也就是说你只要另外加入了其他分子,不管之前有没有相同的分子,都不要合并,按照你的这个gro,后面的几个甲烷必须写到SOL的后面,

1478

帖子

0

威望

4537

eV
积分
6015

Level 6 (一方通行)

10#
 楼主 Author| 发表于 Post on 2021-9-1 17:33:30 | 只看该作者 Only view this author
本帖最后由 牧生 于 2021-9-1 17:36 编辑
lyj714 发表于 2021-8-23 19:15
哪里来的自信说顺序对了,你这个gro中的是H-H-O,但是你spc.itp中的atoms字段可是OHH的顺序啊,除非你把s ...
感谢大佬的回答。做的甲烷水合物,有点像那么回事了。
有一个小疑问,看过一些文献,好多都是在纳秒就结晶了,我都跑了几微秒,才有一点结晶的迹象,这让我有点困惑。






现在有个新出现的问题,我做二氧化碳水合物,将二氧化碳水合物单胞扩展为2×2×2的超胞,然后使用packmol将二氧化碳水合物晶体和水,二氧化碳建立在一个盒子里面,然后将mix.pdb转为mix.gro,并在末尾加上盒子的尺寸,能量最小化以后,二氧化碳水合物晶体就瓦解了,证明我的参数不对。
提示报错为:
Warning: atom name 140 in topol.top and mix.gro does not match (O1 - O)
Warning: atom name 141 in topol.top and mix.gro does not match (O2 - O)
……
Warning: atom name 306 in topol.top and mix.gro does not match (O2 - O)

(more than 20 non-matching atom names)

top中也是按照gro中出现的水-二氧化碳-水-二氧化碳-水-二氧化碳的顺序来编排的。
我检查了一下,发现CO2中的氧原子被识别为水中的氧原子,请问一下如何解决。(原谅我又菜又爱玩)

CO2.itp (511 Bytes, 下载次数 Times of downloads: 39)
mix.gro (118.92 KB, 下载次数 Times of downloads: 58)
topol.top (374 Bytes, 下载次数 Times of downloads: 44)
ffnonbonded.itp (4.78 KB, 下载次数 Times of downloads: 40)
em.mdp (393 Bytes, 下载次数 Times of downloads: 32)
forcefield.itp (929 Bytes, 下载次数 Times of downloads: 30)




又菜又爱玩

1478

帖子

0

威望

4537

eV
积分
6015

Level 6 (一方通行)

11#
 楼主 Author| 发表于 Post on 2021-9-5 16:19:53 | 只看该作者 Only view this author
本帖最后由 牧生 于 2021-9-5 16:22 编辑

问题已经解决,我犯了个很基础的错误。

将[ atomtypes ]写进ffnonbonded.itp,且分开写O1和O2,就能和itp中的O1和O2对应
C1         6     12.01     0.0000    A       0.2757   0.23388
O1         8     16.00     0.0000    A       0.3033   0.66937
O2         8     16.00     0.0000    A       0.3033   0.66937

我真是又菜又爱玩
又菜又爱玩

2

帖子

0

威望

87

eV
积分
89

Level 2 能力者

12#
发表于 Post on 2021-12-9 17:55:40 | 只看该作者 Only view this author
请问楼主,如何用Multiwfn理顺pdb文件中的甲烷和水的原子顺序呀

1478

帖子

0

威望

4537

eV
积分
6015

Level 6 (一方通行)

13#
 楼主 Author| 发表于 Post on 2021-12-9 19:00:48 | 只看该作者 Only view this author
foodie 发表于 2021-12-9 17:55
请问楼主,如何用Multiwfn理顺pdb文件中的甲烷和水的原子顺序呀

http://bbs.keinsci.com/thread-21676-2-1.html

22楼和25楼
又菜又爱玩

2

帖子

0

威望

87

eV
积分
89

Level 2 能力者

14#
发表于 Post on 2021-12-9 21:31:07 | 只看该作者 Only view this author
牧生 发表于 2021-12-9 19:00
http://bbs.keinsci.com/thread-21676-2-1.html

22楼和25楼

感谢 成功解决~

86

帖子

0

威望

1040

eV
积分
1126

Level 4 (黑子)

15#
发表于 Post on 2022-3-15 22:58:22 | 只看该作者 Only view this author
牧生 发表于 2021-8-23 11:20
感谢大佬对我的帮助。
sob大佬已经给予了鱼和渔了。但我仍然遇到了一个看似简单,我却一直无法解决的问 ...

楼主好,我也遇到了同样的问题,想请教一下sob大佬说的方法是什么?

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-13 03:04 , Processed in 0.270941 second(s), 29 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list