24小时热门版块排行榜    

CyRhmU.jpeg
查看: 11279  |  回复: 210
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

bay__gulf

金虫 (著名写手)

刘苏州

[交流] 【在线答疑】关于NAMD/VMD已有27人参与

看到这个板块的在线答疑方式办得很好,这里也尽自己的薄力给大家一些关于NAMD/VMD 原理及使用方面的帮助。

本帖旨在邀请各位朋友集中提问和回答NAMD/VMD的问题,欢迎每个虫子参与提问与回答。我将尽可能的抽出时间来尽我所能回答该帖相关问题。
注意:
1. 请回答问题的虫子勿必回帖时注明:【回复__楼】;
2. 提问的虫子们得到答案后应该重新编辑自己的提问帖,在下面附上:【已解决,答案见__楼】


namd/vmd 官网上的一些好东西,一个免费软件提供的服务比收费软件还要全面。
这里列几个,不保证全。有些链接很深,需要仔细点找。

http://www.ks.uiuc.edu/      首页
http://www.ks.uiuc.edu/Highlights/index.cgi    一些亮点文章 ,可以索要,但不如看看网页的介绍
http://www.ks.uiuc.edu/Research/Categories/   一些分类算例
http://www.ks.uiuc.edu/Research/vmd/  VMD
   http://www.ks.uiuc.edu/Research/vmd/current/docs.html   vmd文档,里面的链接最好都看看
   http://www.ks.uiuc.edu/Research/vmd/plugins/                  VMD Plugin Library
   http://www.ks.uiuc.edu/Research/vmd/script_library/         VMD Script Library
   http://www.ks.uiuc.edu/Research/vmd/mailing_list/            VMD-L Mailing List
   http://biocore.ks.uiuc.edu/bioco ... ips/tclscripts.html     一些脚本,或许对你有用
http://www.ks.uiuc.edu/Research/namd/    NAMD
   http://www.ks.uiuc.edu/Research/namd/wiki/                    NamdWiki
   http://www.ks.uiuc.edu/Research/namd/mailing_list/        NAMD Mailing List
   http://www.ks.uiuc.edu/Research/namd/documentation.html    NAMD Documentation,推荐看一下NAMD 1.5 Programmer's Guide ,我还没有看,谁看完了讲一下。
   http://www.ks.uiuc.edu/Research/namd/announcements.html   NAMD Announcements
   http://www.ks.uiuc.edu/Training/                                                Training
   http://www.ks.uiuc.edu/Training/CaseStudies/index.html           Case Studies
http://www.ks.uiuc.edu/Training/Tutorials/                Tutorials,那些pdf 和代码,常更新。就是说本来对你有用的可能下个月就找不到了

[ Last edited by lei0736 on 2009-12-21 at 18:04 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

引用回帖:
Originally posted by jianchaoyv at 2009-11-15 13:28:
>Main< (vmd) 3 % source ubq.pgn
couldn't read file "ubq.pgn": no such file or directory

请问怎样解决这个问题,初学VMD请高手指教!!

和当前路径有关,按tutorial 一步一步来.
补习pwd, cd, ls 等操作方法
105楼2009-11-15 14:26:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★
zeoliters(金币+2,VIP+0):多谢回帖交流! 12-25 16:44
引用回帖:
Originally posted by jianchaoyv at 2009-11-16 16:49:
------------- Processor 0 Exiting: Called CmiAbort ------------
Reason: FATAL ERROR: can't read "coorList(577)": no such element in array
    while executing
"lindex $coorList($i) ...

这是user define force , 属于namd 的高级用法, 涉及的内容很多,
在一个帖子中讲不完. 你还是先学一下基础的吧.
把这个帖子提到的东西看完再说
http://www.mdbbs.org/thread-13237-1-2.html
108楼2009-11-16 17:05:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★
zeoliters(金币+3,VIP+0):多谢回帖交流! 12-25 16:45
这是一篇教程的例子,模拟的是水分子在加压情况下通过SWNT的过程.
User-Defined Forces 就是用来模拟外压的.
http://www.ks.uiuc.edu/Training/ ... bes-html/index.html
但你的问题是什么
110楼2009-11-16 18:29:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★
lei0736(金币+2,VIP+0):谢谢答疑 11-18 21:16
引用回帖:
Originally posted by jianchaoyv at 2009-11-18 16:36:
我的问题运行.conf时,提示:
------------- Processor 0 Exiting: Called CmiAbort ------------
Reason: FATAL ERROR: can't read "coorList(577)": no such element in array
    while executin ...

你好
不是我不认真回复,
实在是这个问题涉及namd 的高级用法, 不是一两句话能说明白的.
我的建议是你从头开始一步一步学习namd, 有一定基础之后自然会明白.
科研道路没有捷径, 知识也只能靠自己去学, 毕竟MD 不是play game.

你是女生吧, 你留个QQ, 我给你开个单线.
112楼2009-11-18 17:37:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★
lei0736(金币+4,VIP+0):谢谢 12-12 11:35
buried   : resname ALA LEU VAL ILE PHE CYS MET TRP
surface : protein and not buried

hydrophobic : resname ALA LEU VAL ILE PRO PHE MET TRP
polor             : protein and not hydrophobic
115楼2009-12-06 11:44:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★
lei0736(金币+2,VIP+0):谢谢 12-12 11:35
引用回帖:
Originally posted by arthurii at 2009-12-6 14:03:
依照这个规定,那假设对象是核酸还有小分子的话,就需要另外自行定义polar和non polar了吧?

用别的方法定义吧, 比如resid
不要用surface, polor 了, 这个不适合核酸.
117楼2009-12-06 20:49:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★ ★
lei0736(金币+5,VIP+0):感谢分享 12-12 11:36
自己用的, 略作修改
将就看吧
CODE:
set numFrame [molinfo top get numframes]

for { set dt 1 } { $dt < $numFrame/2 } { incr dt } {
  set r2 0
  for { set fr0 1 } { $fr0 < $numFrame/2 } { incr fr0 } {
    set sel [atomselect top "water and name OH2 frame $fr0]
      foreach wti [$sel get index] {
         set wt0 [atomselect top "index $index" frame $fr0]
         set wtt [atomselect top "index $index" frame [expr $fr0 + $dt]]
         set cd0 [measure center $wt0]
         set cdt [measure center $wtt]
         set dr [veclength [vecsub $cdt $cd0] ]
         set r2 [expr $r2 + $dr*$dr]
         $wt0 delete
         $wtt delete
      }
      $sel delete
  }
  set dr2 [expr $r2/$n]
  puts "$n $dt\t$dr2"
}

[ Last edited by bay__gulf on 2009-12-9 at 17:46 ]
120楼2009-12-09 17:37:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lei0736(金币+2):谢谢 2010-01-29 10:04
lei0736(金币+20):应wanlichuan要求转移给你 感谢你赠送邀请码 2010-02-06 10:07
namd用psfgen生成psf文件。参考
http://www.ks.uiuc.edu/Research/namd/2.7b2/ug/node16.html
122楼2010-01-28 17:15:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★
lei0736(金币+4):谢谢 2010-02-13 10:45
引用回帖:
Originally posted by zh1987hs at 2010-02-10 13:36:23:
想请问一下,当处理分子中带有杂原子的时候,比如FAD,或是金属原子的时候,它们的参数如何获得呢?看到VMD中有PARATOOL可以建立缺失原子类型的常数,但是在用的时候总出错,想问下楼主知道PARATOOL可以处理以上两 ...

charmm力场可以处理有机小分子,但大部分没有现成的参数,水合离子除Na,K,Mg,Ca,Zn,Ce,Cl外也没有,需要自己制作,制作方法参考http://www.mdbbs.org/thread-16875-1-1.html

可以以两种方法加入namd计算,如果力场参数形式可以被charmm接受,则可以加到相应par文件里面;如果需要特殊的函数形式,则可以用tclForce,tclBC或External Program Forces等处理,但这类原子数目太多会影响计算速度.

[ Last edited by bay__gulf on 2010-2-10 at 20:29 ]
124楼2010-02-10 20:28:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★ ★
lei0736(金币+5):谢谢 2010-03-07 11:55
1 FEP计算方法参考
Research of Calculation and Simulation, 2010,Vol. 1,No.1,Page22-26.
注意不要到google 学术上去搜索

2 不知道你做什么体系, 关于alchemify
计算体系有两种
a 不涉及共价键断裂,如水中加一个甲烷
b 涉及共价键断裂,如氨基酸突变,改变的部分同不变的部分有bond相连.
无论在ug还是tutorial中都是分开来讨论的.那个alchemify 只是针对后者的工具

两种情况在namd 的算法上没有区别, 只是在建psf 时候, 后者需要注意几个问题, 那个alchemify 就是为了这个才做的. 前者的例子不需要alchemify, 就像平时那样就可以了,当原子消失时候,那些键/角/二面角会自然消失.后者经alchemify 处理,得到特殊的psf 文件,当一个基团还未消失,另一个已经出现时候,连在Ca上的五根 bonds 的相互关系不会出错.

3 关于运行
不清楚你的系统环境, 也不知道具体的报错信息
开启rsh服务进程, 设置无密码登陆
用带local的方法
把namd2 写成绝对路径
不行再发问,写清楚linux版本,环境设置和详细的保存信息
127楼2010-03-07 11:23:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bay__gulf 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见