24小时热门版块排行榜    

查看: 10029  |  回复: 210

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★ ★
lei0736(金币+5,VIP+0):谢谢 7-1 20:45
回【20】
你已经安装成功,
那个错误 No simulation config file specified on command lin,是说没有指定conf 文件。
namd 运行时候,所有的参数,比如温度,压力,步长,总时间等都要在conf 文件中指定。pdb/psf 文件存储分子构型,这2个文件也要在conf 文件中指定。你只是输入namd,它怎么会知道这些参数呢。namd 运行时候没有图形界面,要适应一下哦。

至于会写psf 和 conf文件了,你就学会NAMD 了。
这个要花几个月时间,不是看几个帖子能做到的。
我在这个帖子1楼列出几个链接,里面东西很多,可以下载一些看看。

[ Last edited by bay__gulf on 2009-7-1 at 14:25 ]
21楼2009-07-01 14:23:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wsyg1986

铜虫 (正式写手)

谢谢啦


小木虫(金币+0.5):给个红包,谢谢回帖交流
谢谢啦,我也下拉一个vmd,那个可以打开pdb/psf 文件,可以看到图形,是不是vmd构建分子结构,然后需要什么条件模拟就用conf文件写好,连接就可以完成任务啦。
   非常感谢您!
22楼2009-07-01 16:08:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★
lei0736(金币+3,VIP+0):谢谢 7-1 20:45
引用回帖:
Originally posted by wsyg1986 at 2009-7-1 16:08:
谢谢啦,我也下拉一个vmd,那个可以打开pdb/psf 文件,可以看到图形,是不是vmd构建分子结构,然后需要什么条件模拟就用conf文件写好,连接就可以完成任务啦。
   非常感谢您!

是的,namd只是一个计算模块,建模同数据处理都要在vmd中完成。
所以要先学vmd,建好模型后再用namd 跑动力学,数据出来后再用vmd分析。
namd同vmd 师出同门,现在已经实现无缝链接。
当然,因为文件格式兼容,namd 不只用vmd 建模分析,vmd 也不只是用于namd。

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

yfw7210

金虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
为什么自动产生的参数文件没有内容呢?
像下边:
!>>>>>> Combined CHARMM All-Hydrogen Parameter File for <<<<<<<<<
!>>>>>>>>> CHARMM22 Proteins and CHARMM27 Lipids <<<<<<<<<<
!from
!>>>> CHARMM22 All-Hydrogen Parameter File for Proteins <<<<<<<<<<
!>>>>>>>>>>>>>>>>>>>>>> August 1999 <<<<<<<<<<<<<<<<<<<<<<<<<<<<<
!>>>>>>> Direct comments to Alexander D. MacKerell Jr. <<<<<<<<<
!>>>>>> 410-706-7442 or email: alex,mmiris.ab.umd.edu  <<<<<<<<<
!and !  \\\ CHARMM27 All-Hydrogen Lipid Parameter File ///////
!  \\\\\\\\\ Developmental /////////////////////////
!              Alexander D. MacKerell Jr.
!                     August 1999
! All comments to ADM jr.  email:alex,mmiris.ab.umd.edu
!              telephone: 410-706-7442
!



! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
! + Residue topologies/parameters generated by PARATOOL.                      +
! +                                                                           +
! + PARATOOL is an interactive tool for generation of force field parameters. +
! + and can be obtained free of charge from                                   +
! + http://bioinf.charite.de/biophys/paratool                                 +
! +                                                                           +
! + Author:                                                                   +
! + Jan Saam                                                                  +
! + Institute of Biochemistry                                                 +
! + Charite Berlin                                                            +
! + Germany                                                                   +
! + saam@charite.de                                                           +
! +                                                                           +
! + CITATION:                                                                 +
! + If you use PARATOOL or topologies/parameters generated by this program,   +
! + please cite the following work:                                           +
! + Saam, et al. (2006), ...                                                  +
! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++


! Parameters for components GLY THR
! To be used in combination with:
! F://20090703/par_all27_prot_lipid.inp


BONDS
!
!V(bond) = Kb(b - b0)**2
!
!Kb: kcal/mole/A**2
!b0: A
!
!atom type   Kb         b0
!



ANGLES
!
!V(angle) = Ktheta(Theta - Theta0)**2
!
!V(Urey-Bradley) = Kub(S - S0)**2
!
!Ktheta: kcal/mole/rad**2
!Theta0: degrees
!Kub: kcal/mole/A**2 (Urey-Bradley)
!S0: A
!
!atom types        Ktheta    Theta0    Kub      S0
!



DIHEDRALS
!
!V(dihedral) = Kchi(1 + cos(n(chi) - delta))
!
!Kchi: kcal/mole
!n: multiplicity
!delta: degrees
!
!atom types             Kchi     n    delta
!


IMPROPER
!
!V(improper) = Kpsi(psi - psi0)**2
!
!Kpsi: kcal/mole/rad**2
!psi0: degrees
!note that the second column of numbers (0) is ignored
!
!atom types             Kpsi        0  psi0
!



NONBONDED nbxmod  5 atom cdiel shift vatom vdistance vswitch -
cutnb 14.0 ctofnb 12.0 ctonnb 10.0 eps 1.0 e14fac 1.0 wmin 1.5
                !adm jr., 5/08/91, suggested cutoff scheme
!
!V(Lennard-Jones) = Eps(i,j)*[(Rmin(i,j)/r(i,j))**12 - 2(Rmin(i,j)/r(i,j))**6]
!
!epsilon [kcal/mole]: Eps(i,j) = sqrt(eps(i) * eps(j))
!Rmin/2  [A]:         Rmin(i,j) = Rmin/2(i) + Rmin/2(j)
!
!atom  ignored   epsilon     Rmin/2  ignored   eps,1-4     Rmin/2,1-4
!

END
24楼2009-07-03 15:14:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★ ★ ★
gwdavid(金币+3,VIP+0):谢谢! 7-7 20:12
gwdavid(金币+3,VIP+0):谢谢! 7-7 20:13
回【24】
那些参数不是用某个程序就能得到的。
可以到 http://mackerell.umaryland.edu/CHARMM_ff_params.html 下载。
25楼2009-07-03 16:26:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yfw7210

金虫 (小有名气)

有关VMD插件Paratool


小木虫(金币+0.5):给个红包,谢谢回帖交流
版主,我在用VMD中的QMtool菜单,优化一个两个氨基酸组成的小分子,想得到它的最优化结构,我按照paratool 指南中操作,在写出了高斯输入文件后运行时,没有产生log文件,这是什么问题呢?
26楼2009-07-07 11:05:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★
gwdavid(金币+3,VIP+0):谢谢! 7-7 20:13
【回26楼】
vmd 不是专门的处理gaussian输入输出文件的软件,那个plugin也只是好事之人的习作,
不能指望它做的多么完美,
还是用gview 吧。
27楼2009-07-07 17:26:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

involute

新虫 (初入文坛)

请教NAMD对金电极的处理


小木虫(金币+0.5):给个红包,谢谢回帖交流
请教达人,我正在试图用NAMD做一个纳米金电极和DNA系统的分子动力学模拟,打算采取以下流程:

第一,在DNA的pdb文件中,手工加入一个Gold的residue,把组成纳米电极的全部金原子坐标写入。

第二,修改CHARMM的topology文件,加入对Gold这个residue的topology描述。

第三,在VMD界面下,用修改后的topology文件,和包含DNA、Gold的pdb文件,生成结构描述文件psf。

第四,修改CHARMM的parameter文件,加入对金原子和其他原子相互作用的描述。

第五,把整个系统放入water sphere,并做minimize和equilibrate。

想请教达人,是否可以这样做?如果可以,对Gold的topology和parameter如何描述?这些参数在哪儿找?

谢谢!
28楼2009-07-16 10:31:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★
lei0736(金币+3,VIP+0):谢谢 7-16 12:25
【回28】
步骤正确,但charmm力场对金属原子支持不强,par文件很难写。
这方面好像还没有好的工作。
要么减少一些关键的参数,要么换个力场吧
29楼2009-07-16 12:10:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

involute

新虫 (初入文坛)

回29楼


小木虫(金币+0.5):给个红包,谢谢回帖交流
回29楼:
谢谢达人指点!

      我看的参考文献上,对DNA游动通过纳米电极,他们采取的是CHARMM27力场,所以,我想他们是不是找到了跟金原子相关的参数,自己加进去了?

      另外,我本人是学微电子的,这是第一次做分子动力学模拟,是看了NAMD的tutorial,结合这个问题,估摸着这么一套方案[见27楼],心里完全没底,请教达人,确定是这么做吗?

      谢谢!
30楼2009-07-16 14:29:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bay__gulf 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 工材01送了吗? +10 xiaopang8958 2024-05-25 18/900 2024-05-31 10:15 by niuyunsong
[基金申请] 博后特别资助状态变化 +19 随梦而飞2017 2024-05-30 30/1500 2024-05-31 10:02 by 喵呜呜_c
[论文投稿] 纠结选哪一个期刊,电化学领域 50+7 Freya163 2024-05-28 8/400 2024-05-31 09:20 by 莱茵润色
[基金申请] B口人才项目 +6 WOWO159357 2024-05-29 9/450 2024-05-31 07:04 by dshy
[硕博家园] 每到中夜,情难自抑 +29 sioc-sunj 2024-05-28 53/2650 2024-05-30 22:09 by meier8
[基金申请] 化学口B0109(高分子合成),拿青年基金一般需要怎样的文章水平? +10 salmon95 2024-05-30 17/850 2024-05-30 21:58 by 3115321
[教师之家] 被惯着的学生终究要吃大亏 +31 535743368 2024-05-24 34/1700 2024-05-30 18:18 by kanyechris
[文学芳草园] 物是人非 +4 myrtle 2024-05-30 4/200 2024-05-30 15:05 by mapenggao
[基金申请] 基金上会 +15 mrKiller 2024-05-25 24/1200 2024-05-30 13:06 by mrKiller
[电化学] 有锂金属软包电池代工公司代做软包嘛 +5 yyywj995 2024-05-26 6/300 2024-05-30 08:33 by 那年的旧城
[硕博家园] 论大家对6070后普通教授导师的看法 +5 SNaiL1995 2024-05-28 9/450 2024-05-29 21:54 by SNaiL1995
[考博] 24年博士招生 +7 abinit432 2024-05-27 9/450 2024-05-29 17:16 by 中国银河
[论文投稿] materials letter +6 烟雨盛世 2024-05-24 6/300 2024-05-29 15:33 by ¥笑傲江湖¥
[考研] 研0 +3 徐小七七 2024-05-25 5/250 2024-05-29 13:49 by 1158057902
[论文投稿] 真急着毕业,CPB主编终审17天了,邮件催稿了两次,就是一点动静没有 5+3 kkkk夏 2024-05-28 6/300 2024-05-29 11:18 by hitsdu
[找工作] 找工作如此之难 +6 探123 2024-05-25 6/300 2024-05-29 08:11 by SNaiL1995
[有机交流] D-阿拉伯糖-1,4-内酯的合成 2+4 Leeu55 2024-05-24 12/600 2024-05-28 21:17 by Leeu55
[基金申请] E05现在涵评快结束了吧 +7 shrz98 2024-05-24 8/400 2024-05-28 10:35 by xuel2011
[有机交流] 奇怪的物质 100+4 桃桃PXS 2024-05-27 7/350 2024-05-28 10:22 by 091602
[硕博家园] 我是很理想化一人 +6 hahamyid 2024-05-26 6/300 2024-05-27 18:13 by 大飞鱼鱼鱼
信息提示
请填处理意见