版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

vip_zzc

銅蟲 (小有名氣)

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

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

我用amber力場（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è)力場ff99SB，ff12SB，及 AmberTool1.2、 AmberTool15兩個(gè)版本，用來生成.prmtop （leap 顯示存在 IMPROPER）
結(jié)果NAMD讀出的IMPROPER都是0 ，雖然NAMD可以正常運(yùn)行！

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

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

求解求解～

回復(fù)此樓

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

科研之寶，所向披靡

» 猜你喜歡

310求調(diào)劑已經(jīng)有5人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有17人回復(fù)
材料工程專碩求調(diào)劑已經(jīng)有3人回復(fù)
337求調(diào)劑已經(jīng)有5人回復(fù)
材料與化工328分調(diào)劑已經(jīng)有9人回復(fù)
318一志愿吉林大學(xué)生物與醫(yī)藥求調(diào)劑已經(jīng)有5人回復(fù)
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有6人回復(fù)
085600材料與化工調(diào)劑已經(jīng)有4人回復(fù)
327求調(diào)劑已經(jīng)有4人回復(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的回帖

查看全部 5 個(gè)回答

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力場時(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力場時(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力場參數(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的回帖

簡單回復(fù)

kiss200855372樓

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

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

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

查看全部 5 個(gè)回答

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

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

[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +4	兔兔兔111223314 2026-03-29	4/200	2026-03-30 06:55 by ilovexiaobin
[考研] 309求調(diào)劑 +8	誰不是少年 2026-03-29	8/400	2026-03-30 00:55 by 我是小康
[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +5	Shern—- 2026-03-28	5/250	2026-03-29 23:32 by 飛行日記西
[考研] 南京大學(xué)化學(xué)調(diào)劑 +9	景隨風(fēng) 2026-03-29	14/700	2026-03-29 21:22 by chixmc
[考研] 材料專碩調(diào)劑 +5	椰椰。 2026-03-29	5/250	2026-03-29 20:50 by 唐沐兒
[考研] 考研調(diào)劑 +7	小蠟新筆 2026-03-29	7/350	2026-03-29 19:00 by 學(xué)員8dgXkO
[考研] 277跪求調(diào)劑 +6	1915668 2026-03-27	10/500	2026-03-29 16:03 by 王亮_大連醫(yī)科大
[考研] 315求調(diào)劑 +4	akie... 2026-03-28	5/250	2026-03-28 21:05 by zhq0425
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +5	niko- 2026-03-27	5/250	2026-03-28 14:25 by 唐沐兒
[考研] 0703一志愿9，初試成績：338，四六級(jí)已過，有科研經(jīng)歷，求調(diào)劑！ +4	Zuhui0306 2026-03-25	4/200	2026-03-28 13:07 by 唐沐兒
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8	邱gl 2026-03-27	8/400	2026-03-28 12:42 by 唐沐兒
[考研] 286求調(diào)劑 +4	丟掉懶惰 2026-03-27	7/350	2026-03-28 08:07 by baoball
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4	Ripcord上岸 2026-03-27	4/200	2026-03-27 17:00 by zhyzzh
[考研] 求調(diào)劑323材料與化工 +7	1124361 2026-03-24	7/350	2026-03-27 10:22 by wangjy2002
[考研] 求調(diào)劑 +6	林之夕 2026-03-24	6/300	2026-03-27 08:38 by hypershenger
[考研] 284求調(diào)劑 +11	junqihahaha 2026-03-26	12/600	2026-03-27 04:37 by wxiongid
[考研] 342求調(diào)劑 +3	加油a李zs 2026-03-26	3/150	2026-03-27 00:29 by wxiongid
[考研] 321求調(diào)劑 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚的貓
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5	星沉uu 2026-03-26	6/300	2026-03-26 19:02 by macy2011
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4	鋼鐵大炮 2026-03-24	4/200	2026-03-26 16:28 by dick_runner

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

vip_zzc

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

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

» 猜你喜歡

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

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

smutao

vip_zzc

vip_zzc

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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

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

查看全部散金貼