24小时热门版块排行榜    

查看: 2316  |  回复: 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的回帖

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的回帖
查看全部 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的回帖

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的回帖

yaozhq

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
3楼: Originally posted by broken1999 at 2013-05-21 19:25:18
感谢您的回复

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

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

通过你给的原子类型定义 因为力场是通过定义一个原子类型 然后拟合其相关的各种参数来实现计算的 很简单 就类似于一个表格 你在算一对儿原子作用的时候 分别查看他们的原子类型和他们的成键关系 然后加载力场参数去带公式计算就好了 几乎所有动力学软件都用拓扑文件来保存原子类型和成键关系的  软件实现的细节会远比你想象的复杂 因为期间有各种不同级别的算法的存在 包括并行化的代码 直接看会不大容易理清其中的结构的

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

5楼2013-05-23 09:17:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +3 Tide man 2026-08-14 4/200 2026-08-14 16:33 by Yeuchan
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +4 Ldrop2023 2026-08-13 4/200 2026-08-14 14:50 by Tide man
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[基金申请] 好奇怪的filecode +6 布布和一二 2026-08-08 7/350 2026-08-14 13:38 by Equinoxhua
[基金申请] filecode +5 cratir 2026-08-14 9/450 2026-08-14 12:25 by 医学老男孩
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 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
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
信息提示
请填处理意见