24小时热门版块排行榜    

查看: 2586  |  回复: 5
【悬赏金币】回答本帖问题,作者ccc阿鑫将赠送您 10 个金币

ccc阿鑫

新虫 (小有名气)

[求助] VMD生成的石墨烯pbd文件如何正确转化为psf 已有1人参与

在VMD中使用工具可以得到石墨烯,通过保存坐标可以轻松得到pdb文件。但是psf文件怎么得到呢?
我试了VMD的autopsf工具,会自动闪退,没办法得到psf。
然后我在Tcl中,使用一下命令可以得到psf
package require topotools
topo retypebonds
topo guessangles
topo guessdihedrals mol reanalyze 0
animate write psf name.psf
生成的psf看起来都没有问题,但是在仿真时候会报错,"FATAL ERROR: UNABLE TO FIND ANGLE PARAMETERS FOR C C C."

有没有大神知道如何得到正确的psf,十分感谢

@smutao @oxox6085 发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

先点鼠标:
Extensions --  Modeling --  Nanotube Builder  进入石墨烯和碳纳米管的建模界面
输入你想要的尺寸,然后点击最下面的Generate Sheet就可以建好石墨烯模型
进入 VMD TKconsole 输入命令:
set all [atomselect top all]
$all writepdb graphene.pdb
$all writepsf graphene.psf

既可把石墨烯模型保存在你的电脑里,位置在vmd的默认路径,你可以通过  pwd  命令查看
2楼2024-03-12 09:24:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

"FATAL ERROR: UNABLE TO FIND ANGLE PARAMETERS FOR C C C."

这个错误貌似不是PSF的问题,是你在力场参数文件里没有正确设置的问题
3楼2024-03-12 09:25:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ccc阿鑫

新虫 (小有名气)

引用回帖:
2楼: Originally posted by jiaoyixiong at 2024-03-12 09:24:00
先点鼠标:
Extensions --  Modeling --  Nanotube Builder  进入石墨烯和碳纳米管的建模界面
输入你想要的尺寸,然后点击最下面的Generate Sheet就可以建好石墨烯模型
进入 VMD TKconsole 输入命令:
set all  ...

谢谢解答! 但是这样只能保存完整的石墨烯膜,我想做石墨烯孔,类似DNA测序那种,我是通过保存时候,选中部分原子做的。这样的话也可以用这个方式保存吗

发自小木虫Android客户端
4楼2024-03-28 00:07:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ccc阿鑫

新虫 (小有名气)

引用回帖:
3楼: Originally posted by jiaoyixiong at 2024-03-12 09:25:16
"FATAL ERROR: UNABLE TO FIND ANGLE PARAMETERS FOR C C C."

这个错误貌似不是PSF的问题,是你在力场参数文件里没有正确设置的问题

谢谢,我检查一下,c  c  c这种力场文件没有,我可以自己改力场文件加上这部分数据嘛

发自小木虫Android客户端
5楼2024-03-28 00:08:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

引用回帖:
4楼: Originally posted by ccc阿鑫 at 2024-03-28 00:07:34
谢谢解答! 但是这样只能保存完整的石墨烯膜,我想做石墨烯孔,类似DNA测序那种,我是通过保存时候,选中部分原子做的。这样的话也可以用这个方式保存吗
...

可以的,你输入
set all [atomselect top “sqrt(x*x+y*y)>10”]
就可以保持纳米孔直径为10艾米的石墨烯片。

“想做石墨烯孔,类似DNA测序那种”这个在NAMD 官网上有这个具体的例子教程,你可以直接去下载
6楼2024-04-17 16:12:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ccc阿鑫 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-10-01 21:59 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 CfXuS1rDhLYN 2026-09-28 4/200 2026-10-01 09:23 by Equinoxhua
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:45 by Z9YWQ5EAO3qp
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见