24小时热门版块排行榜    

查看: 12717  |  回复: 210

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★ ★
lei0736(金币+5,VIP+0):谢谢 7-20 09:22
【回30】
你的思路非常正确,最大的问题是力场参数。

namd软件或是charmm力场做固体也大都是一些非金属材料形成的原子晶体。
看到过关于DNA穿越纳米孔道的也只是Si 的化合物,如Si3N4, SiO2 等,几乎没有看到过纯金属的。

charmm 中有一些金属的力场,那些全都是水化金属离子的。金属单质的势函数极为复杂,可能charmm 力场本身就表达不了。更不要说金属同有机物的相互作用了。

做金属单质的有不少用lammps 的,你可以关注一下

[ Last edited by bay__gulf on 2009-7-16 at 15:11 ]
31楼2009-07-16 15:05:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

involute

新虫 (初入文坛)

回31


小木虫(金币+0.5):给个红包,谢谢回帖交流
非常感谢您的回答!

想继续请教几个问题:
第一,你说的力场问题,我仔细查了参考文献,有这样一篇:Rappe, A. K.; Casewit, C. J.; Colwell, K. S. Goddard, W. A., III; Skiff, W. M. J. Am. Chem. Soc. 1992, 114, 10024. 它里面提出了一个Uniform force field (UFF),给出了各种势场的公式、参数,原则上可以计算CHARMM力场参数。

第二,一开始的时候,往DNA的pdb文件里,写入构成电极的全部金原子坐标位置,我打算利用materials studio的图形界面,画出组成电极的原子构型,然后保存成pdb文件,最后把这个pdb文件和DNA的pdb文件合二为一。你看可不可行?

第三,在topology文件里,把全部金原子写进一个residue,是不是太麻烦了?有无可能选取一个超晶胞(supercell)写residue,然后在生成psf时,其它所有金原子都会按照这个超晶胞,生成晶格结构?
32楼2009-07-19 19:20:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★ ★
lei0736(金币+5,VIP+0):谢谢 7-20 09:23
1
UFF 是求全型力场,几乎可以计算任何体系,但精度很差。
charmm 力场只包含键、角、二面角、交叉项(限制肽平面)、vdW 和静电力,这么力场形式可能太简单处理不了金属单质。
只是建议,我没有做过这些,你可以按自己的思路试一下。

2
可以,分别制作电极和DNA的pdb/psf,然后在psfgen 里面就能合并在一起

3
对于单个原子的,像水溶液中的离子,都是每个原子一个residue 的。
如果你用charmm形式的力场,可以每个residue 只有一个原子,这样top 文件就好写了。

你的很大的一个问题是力场。即要精确可信,又要兼容namd。剩下的问题,比如如何写top 文件等,都要等到力场确定后才能考虑。
33楼2009-07-19 20:50:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★
lei0736(金币+4,VIP+0):谢谢 7-21 10:47
找到一篇处理金-蛋白的文献,有vdW力和静电力的处理方法,
http://www.ks.uiuc.edu/Research/gbp/  or
Rosemary Braun, Mehmet Sarikaya, and Klaus Schulten.
Genetically engineered gold-binding polypeptides:
Structure prediction and molecular dynamics.
Journal of Biomaterials Science, 13:747-758, 2002.

不过他们用的是一个很大的金板,
如果你是针尖的话可能要变通一下才能用
34楼2009-07-21 08:48:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tangjikede

新虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
用vmd查看的坐标文件,如果是同一个个结构的连续变型,譬如说,碳纳米管的拉伸,要看到拉伸的整个过程,我已将xyz坐标文件的开头写上这根碳纳米管的原子个数,然后空一行,然后第一列是大写的C,第一步就是没拉伸的碳纳米管结构,第二步的结构就很乱,有些键应该断掉,但是还连着。是否在VMD中有断键的设置?
就是说两个原子超过多大距离以后就没有键了,不知道怎样能解决。
35楼2009-07-26 10:53:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★
lei0736(金币+2,VIP+0):谢谢 7-27 10:31
vmd 根据psf 或top 文件判断键连信息,在representation 框中dynamic bonds 试试看
36楼2009-07-26 11:26:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tangjikede

新虫 (小有名气)

Thanks.
37楼2009-08-02 15:16:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tangjikede

新虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
请问怎样把连续变化的结构变形,制作成动画,现在只能通过手动拉动滑条看到。我想用VMD 里extensions -->Visualization ---> Movie maker. 但是点击 Make movie 后只自动生成一些bmp 文件,会弹出一条信息,Program not found  could not locate 'videomach.exe' ......, 是否缺了一些东西,
不知道怎样能制作一个视频频文件,然后就可以自动播放卡看到整个过程。
38楼2009-08-02 15:29:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★
lei0736(金币+3,VIP+0):谢谢 8-4 20:15
VideoMach 是一款可以把图片做成视频的软件,vmd 本身不带。
这个很容易下载。

linux 下可以用bmp2avi 生成视频。
39楼2009-08-02 19:06:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wld0328_17

新虫 (小有名气)

新手提问 ,在线等


小木虫(金币+0.5):给个红包,谢谢回帖交流
我们学校有个SGI Altix XE 1300 的集群, 应该下载哪个版本的NAMD?
我试了几个Linux的,当输入命令 namd2 sim_short.conf 时总是提示“namd2 cannot found ”, 这是为什么?
脚踏实地。
40楼2009-08-05 09:46:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bay__gulf 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见