计算化学公社

标题: xtc轨迹文件输出的PDB异常 [打印本页]

作者
Author: 1154975925 时间: 2024-11-8 21:32
标题: xtc轨迹文件输出的PDB异常
我将gromacs模拟输出的gro文件载入VMD后，再load xtc 文件，之后执行下面的代码输出所有帧的PDB，下图是其中某一帧的PDB，结构是异常的，不知道为啥，下面是我执行的代码
set all [atomselect top "all"]
set nf [molinfo top get numframes]
for { set i 0 } { $i < $nf } { incr i } {
$all frame $i
$all writepdb $i.pdb
}

(, 下载次数 Times of downloads: 5)

作者
Author: sobereva 时间: 2024-11-8 21:54
弄清楚是不是跨盒子导致键的判断问题导致的

作者
Author: 1154975925 时间: 2024-11-9 12:44

sobereva 发表于 2024-11-8 21:54
弄清楚是不是跨盒子导致键的判断问题导致的

可是已经PBC了呀

作者
Author: sobereva 时间: 2024-11-9 18:08

1154975925 发表于 2024-11-9 12:44
可是已经PBC了呀

莫名其妙
正因为用了PBC，才会出现跨盒子两头的键

作者
Author: 1154975925 时间: 2024-11-10 18:41

sobereva 发表于 2024-11-9 18:08
莫名其妙
正因为用了PBC，才会出现跨盒子两头的键

我试了一下,没有PBC 也是这种情况，会不会是周期性边界条件没设置好，比如gmx editconf中的d，我这里设置的是1，或许需要调大一些？

作者
Author: student0618 时间: 2024-11-10 19:06
虽然不知为什么要pdb，但trjconv也可以输出pdb哦。

作者
Author: 1154975925 时间: 2024-11-11 10:26

student0618 发表于 2024-11-10 19:06
虽然不知为什么要pdb，但trjconv也可以输出pdb哦。

和代码效果一样同样存在这个问题

作者
Author: Seyilaxa 时间: 2024-11-11 11:00
是不是结构中有不同的chain没区分开来

作者
Author: student0618 时间: 2024-11-11 13:55
本帖最后由 student0618 于 2024-11-11 14:01 编辑

也可以试试删掉所有CONECT 行让PyMOL根据距离判断成键。

作者
Author: 1154975925 时间: 2024-11-11 16:18

student0618 发表于 2024-11-11 13:55
也可以试试删掉所有CONECT 行让PyMOL根据距离判断成键。

里面没有conect

作者
Author: 1154975925 时间: 2024-11-11 16:23

Seyilaxa 发表于 2024-11-11 11:00
是不是结构中有不同的chain没区分开来

说的很对里面的A B链都变成的X 很奇怪

作者
Author: Seyilaxa 时间: 2024-11-12 13:14

1154975925 发表于 2024-11-11 16:23
说的很对里面的A B链都变成的X 很奇怪

之前我也遇到过，手动把链标识改过来就行了

作者
Author: 1154975925 时间: 2024-11-12 19:39

Seyilaxa 发表于 2024-11-12 13:14
之前我也遇到过，手动把链标识改过来就行了

感谢

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