24小时热门版块排行榜    

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

broken1999

木虫 (正式写手)

[求助] NAMD软件在代码层面上是如何计算蛋白质势能的?

NAMD软件
CHARMM力场

蛋白质势能的计算是进行能量最小化的基础
想深入软件代码层面,完全掌握它的计算过程

查过资料,知道势能由以下几部分组成:键伸缩能bond stretching、键角弯曲能angle bending、二面角扭曲能torsion rotation、异常扭曲improper torsion、离平面弯曲Out-of –plane Bending、交叉项cross term、范德华相互作用能van der Waals interactions、静电相互作用能electrostatic contributions几部分组成

但是每个部分里面都涉及到好多参数,这些参数又取决于CHARMM力场
看过namdenergy.tcl的源代码,惘然而无所得,

几个问题各位大牛给予指导
1 NAMD是如何挑选合适的参数进行对势计算的
2 NAMD是否对运算过程进行过优化或近似?难道它就穷举所有原子对势?
3 有没有哪本书,或者资料对NAMD势能计算过程进行过讲解? 求指点
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yaozhq

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
broken1999: 金币+10, ★★★很有帮助 2013-06-06 20:07:47
NAMD的全名如果知道的话 你就不会问第二个问题了 对于求解大规模的动力学模拟 优化是必须的 而且程度非常高 最起码并行化就需要调整通讯和计算 任务分配 调度汇总等等问题 而近似是动力学天生的特性 在尽量保证精度和速度下 其实做了非常非常多的近似和简化处理 要理解这个就去多看看力场和算法的书吧 保证颠覆你的现有完美主义想法
6楼2013-05-23 09:25:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

fhh2626

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
broken1999: 金币+10, ★★★很有帮助 2013-06-06 20:07:16
1、参数自然是包括在CHARMM力场当中的,你要先知道CHARMM力场的格式
2、NAMD只计算一定范围内的原子(如10A),长程静电作用可以用一定的方法近似(如PME)
3、基本上任何关于分子模拟的书都会花几页讲一下力场的。。。比如陈正隆的那本

» 本帖已获得的红花(最新10朵)

2楼2013-05-21 14:29:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

broken1999

木虫 (正式写手)

引用回帖:
2楼: Originally posted by fhh2626 at 2013-05-21 14:29:43
1、参数自然是包括在CHARMM力场当中的,你要先知道CHARMM力场的格式
2、NAMD只计算一定范围内的原子(如10A),长程静电作用可以用一定的方法近似(如PME)
3、基本上任何关于分子模拟的书都会花几页讲一下力场的 ...

感谢您的回复

1 软件层面上,是通过什么东西识别原子类型,挑选对应参数的?有没有哪些资料对NAMD的实现细节进行过阐述?

2 陈正隆的书我看过,确实有讲,但是很简略,也没提NAMD的实现细节
3楼2013-05-21 19:25:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fhh2626

木虫 (正式写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by broken1999 at 2013-05-21 19:25:18
感谢您的回复

1 软件层面上,是通过什么东西识别原子类型,挑选对应参数的?有没有哪些资料对NAMD的实现细节进行过阐述?

2 陈正隆的书我看过,确实有讲,但是很简略,也没提NAMD的实现细节...

1、Top文件,你把CHARMM力场的格式搞明白就不会要什么具体实现细节了
2、同1
4楼2013-05-22 10:34:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见