vip_zzc

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 105
帖子: 87
在線: 49小時(shí)
蟲號(hào): 758728

[交流] AMBER力場(chǎng) improper 參數(shù) 在NAMD中不能被讀入（0 IMPROPER）

我用amber力場(chǎng)（ff12SB, ff99SB兩個(gè)都試過了）在NAMD中跑MD，在NAMD輸出log文件中，顯示 0 IMPROPER，如下：

PARM file in AMBER 7 format
Warning: Skipping ATOMIC_NUMBER in parm file while seeking MASS.
Warning: Encounter 10-12 H-bond term
Warning: Found 87876 H-H bonds.
Info: LONG-RANGE LJ: APPLYING ANALYTICAL CORRECTIONS TO ENERGY AND PRESSURE
Info: LONG-RANGE LJ: AVERAGE A AND B COEFFICIENTS -1.95755e+07 AND -18374.1
Info: SUMMARY OF PARAMETERS:
Info: 104 BONDS
Info: 219 ANGLES
Info: 139 DIHEDRAL
Info: 0 IMPROPER
Info: 0 CROSSTERM
Info: 0 VDW
Info: 325 VDW_PAIRS
Info: 0 NBTHOLE_PAIRS
Info: READING EXTRA BONDS FILE restraints/step6.1_equilibration.dihe
Info: READ 524 EXTRA DIHEDRALS
Info: READ 524 EXTRA IMPROPERS
Info: TIME FOR READING PDB FILE: 9.53674e-07
Info:
Info: ****************************
Info: STRUCTURE SUMMARY:
Info: 381982 ATOMS
Info: 381366 BONDS
Info: 221289 ANGLES
Info: 455628 DIHEDRALS
Info: 0 IMPROPERS
Info: 0 CROSSTERMS
Info: 914447 EXCLUSIONS
Info: 330613 RIGID BONDS
Info: 815333 DEGREES OF FREEDOM
Info: 139245 HYDROGEN GROUPS
Info: 4 ATOMS IN LARGEST HYDROGEN GROUP
Info: 139245 MIGRATION GROUPS
Info: 4 ATOMS IN LARGEST MIGRATION GROUP
Info: TOTAL MASS = 2.33708e+06 amu
Info: TOTAL CHARGE = -0.00574125 e
Info: MASS DENSITY = 0.994389 g/cm^3
Info: ATOM DENSITY = 0.0978737 atoms/A^3
Info: *****************************

如上所示，在parameter summary 和 structure summary中，IMPROPER 一項(xiàng)為 0 ��！
我試過兩個(gè)力場(chǎng)ff99SB，ff12SB，及 AmberTool1.2、 AmberTool15兩個(gè)版本，用來生成.prmtop （leap 顯示存在 IMPROPER）
結(jié)果NAMD讀出的IMPROPER都是0 ，雖然NAMD可以正常運(yùn)行！

經(jīng)測(cè)試，發(fā)現(xiàn)：若用CHARMM力場(chǎng)在NAMD跑，確實(shí)給出了相應(yīng)的IMPROPER數(shù)
請(qǐng)問有沒有人注意到這個(gè)奇葩的現(xiàn)象呢？

1. 難道NAMD讀AMBER力場(chǎng)時(shí)已將IMPROPER納入到DIHEDRAL了？
2. 還是壓根就是程序本身的bug，果真如此，NAMD跑出的結(jié)果就是不對(duì)了？

求解求解～

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

科研之寶，所向披靡

» 猜你喜歡

309求調(diào)劑已經(jīng)有8人回復(fù)
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有8人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
求化學(xué)調(diào)劑已經(jīng)有5人回復(fù)
085600，材料與化工321分求調(diào)劑已經(jīng)有10人回復(fù)
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分已經(jīng)有5人回復(fù)
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑已經(jīng)有3人回復(fù)
085600，專業(yè)課化工原理，320分求調(diào)劑已經(jīng)有4人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有16人回復(fù)
290求調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2014-02-18 15:46:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

smutao

禁蟲 (著名寫手)

★ ★ ★
vip_zzc(金幣+2): 謝謝參與
月只藍(lán): 金幣+2, 鼓勵(lì)交流！ 2014-02-18 18:27:42

本帖內(nèi)容被屏蔽

3樓2014-02-18 17:28:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

vip_zzc

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 105
帖子: 87
在線: 49小時(shí)
蟲號(hào): 758728

引用回帖:

3樓: Originally posted by smutao at 2014-02-18 17:28:47
請(qǐng)樓主參見
RUN NAMD with Amber Force Field
http://ambermd.org/namd/namd_amber.html
里面好像有講到樓主的那個(gè)問題

這個(gè)之前就看過了，主要說明是用AMBER力場(chǎng)時(shí)，在NAMD配置文件應(yīng)該如何對(duì)相應(yīng)選項(xiàng)進(jìn)行設(shè)置的問題。
但并未提到我所描述的問題。。
不過，仍然感謝你的回復(fù)~~

贊一下

回復(fù)此樓

4樓2014-02-18 20:31:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

vip_zzc

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 105
帖子: 87
在線: 49小時(shí)
蟲號(hào): 758728

★
月只藍(lán): 金幣+1, 鼓勵(lì)交流！ 2014-02-19 11:01:55

引用回帖:

4樓: Originally posted by vip_zzc at 2014-02-18 20:31:52
這個(gè)之前就看過了，主要說明是用AMBER力場(chǎng)時(shí)，在NAMD配置文件應(yīng)該如何對(duì)相應(yīng)選項(xiàng)進(jìn)行設(shè)置的問題。
但并未提到我所描述的問題。。
不過，仍然感謝你的回復(fù)~~

...

grep FLAG .prmtop 后
%FLAG TITLE
%FLAG POINTERS
%FLAG ATOM_NAME
%FLAG CHARGE
%FLAG ATOMIC_NUMBER
%FLAG MASS
%FLAG ATOM_TYPE_INDEX
%FLAG NUMBER_EXCLUDED_ATOMS
%FLAG NONBONDED_PARM_INDEX
%FLAG RESIDUE_LABEL
%FLAG RESIDUE_POINTER
%FLAG BOND_FORCE_CONSTANT
%FLAG BOND_EQUIL_VALUE
%FLAG ANGLE_FORCE_CONSTANT
%FLAG ANGLE_EQUIL_VALUE
%FLAG DIHEDRAL_FORCE_CONSTANT
%FLAG DIHEDRAL_PERIODICITY
%FLAG DIHEDRAL_PHASE
%FLAG SCEE_SCALE_FACTOR
%FLAG SCNB_SCALE_FACTOR
%FLAG SOLTY
%FLAG LENNARD_JONES_ACOEF
%FLAG LENNARD_JONES_BCOEF
%FLAG BONDS_INC_HYDROGEN
%FLAG BONDS_WITHOUT_HYDROGEN
%FLAG ANGLES_INC_HYDROGEN
%FLAG ANGLES_WITHOUT_HYDROGEN
%FLAG DIHEDRALS_INC_HYDROGEN
%FLAG DIHEDRALS_WITHOUT_HYDROGEN
%FLAG EXCLUDED_ATOMS_LIST
%FLAG HBOND_ACOEF
%FLAG HBOND_BCOEF
%FLAG HBCUT
%FLAG AMBER_ATOM_TYPE
%FLAG TREE_CHAIN_CLASSIFICATION
%FLAG JOIN_ARRAY
%FLAG IROTAT
%FLAG RADIUS_SET
%FLAG RADII
%FLAG SCREEN
%FLAG IPOL
我們會(huì)發(fā)現(xiàn)在amber .prmtop文件中，并沒有IMPROPER這一標(biāo)識(shí)，說明AMBER應(yīng)該是將
IMPROPER 和 PROPER 整合在一起了，NAMD根據(jù)FLAG讀取amber力場(chǎng)參數(shù)，故只能得到
DIHEDRAL，我認(rèn)為 NAMD .log 輸出中給出的 DIHEDRAL = PROPER + IMPROPER
（當(dāng)然，這是對(duì)amber FF作為輸入而言）

轉(zhuǎn) http://simulation.haotui.com/thread-45212-1-1.html

贊一下(1人)

回復(fù)此樓

5樓2014-02-18 21:22:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

簡(jiǎn)單回復(fù)

kiss200855372樓

2014-02-18 15:55 回復(fù)

vip_zzc(金幣+2): 謝謝參與

1 [ 發(fā)自小木蟲客戶端 ]

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 327求調(diào)劑 +3	小卡不卡. 2026-03-29	3/150	2026-03-29 22:29 by 無際的草原
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5	hanamiko 2026-03-29	5/250	2026-03-29 16:46 by 學(xué)員8dgXkO
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12	c297914 2026-03-28	12/600	2026-03-29 16:21 by Serene1974
[考研] 學(xué)碩274求調(diào)劑 +9	Li李魚 2026-03-26	9/450	2026-03-28 21:42 by bymhappy
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4	小張zxy張 2026-03-26	8/400	2026-03-28 19:25 by lbsjt
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +3	nhj_ 2026-03-25	3/150	2026-03-28 13:24 by Iveryant
[考研] 299求調(diào)劑 +7	嗯嗯嗯嗯2 2026-03-27	7/350	2026-03-28 13:09 by 唐沐兒
[考研] 085701環(huán)境工程，267求調(diào)劑 +16	minht 2026-03-26	16/800	2026-03-28 12:16 by zllcz
[考研] 292求調(diào)劑 +14	鵝鵝鵝額額額額?/a> 2026-03-25	15/750	2026-03-28 08:45 by WYUMater
[考研] 286求調(diào)劑 +4	丟掉懶惰 2026-03-27	7/350	2026-03-28 08:07 by baoball
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級(jí)，六級(jí)521，求調(diào)劑 +4	美味可樂雞翅 2026-03-26	4/200	2026-03-27 15:27 by 星空星月
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5	XQTJZ 2026-03-26	5/250	2026-03-27 14:49 by 狂炫麥當(dāng)當(dāng)
[考研] 化學(xué)308分求調(diào)劑 +8	你好明天你好 2026-03-23	9/450	2026-03-27 14:01 by 楊光于青云
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6	kongweilin 2026-03-26	8/400	2026-03-27 10:18 by FF_16
[考研] 求調(diào)劑，一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +4	@taotao 2026-03-26	5/250	2026-03-27 08:10 by hypershenger
[考研] 281求調(diào)劑 +6	Koxui 2026-03-24	7/350	2026-03-26 15:37 by 無際的草原
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3	江往賣魚 2026-03-26	3/150	2026-03-26 09:42 by 王小欠i
[考研] 290分調(diào)劑求助 +3	吉祥止止陳 2026-03-25	3/150	2026-03-25 19:58 by barlinike
[考研] 0854AI CV方向招收調(diào)劑 +4	章小魚567 2026-03-23	4/200	2026-03-25 17:04 by CoderLoser
[考研] 293求調(diào)劑 +7	加一一九 2026-03-24	7/350	2026-03-25 12:02 by userper