| 9 | 1/1 | 返回列表 |
| 查看: 1864 | 回復(fù): 8 | ||
broken1999木蟲 (正式寫手)
|
[求助]
NAMD軟件在代碼層面上是如何計(jì)算蛋白質(zhì)勢能的?
|
|
NAMD軟件 CHARMM力場 蛋白質(zhì)勢能的計(jì)算是進(jìn)行能量最小化的基礎(chǔ) 想深入軟件代碼層面,完全掌握它的計(jì)算過程 查過資料,知道勢能由以下幾部分組成:鍵伸縮能bond stretching、鍵角彎曲能angle bending、二面角扭曲能torsion rotation、異常扭曲improper torsion、離平面彎曲Out-of –plane Bending、交叉項(xiàng)cross term、范德華相互作用能van der Waals interactions、靜電相互作用能electrostatic contributions幾部分組成 但是每個(gè)部分里面都涉及到好多參數(shù),這些參數(shù)又取決于CHARMM力場 看過namdenergy.tcl的源代碼,惘然而無所得, 幾個(gè)問題各位大牛給予指導(dǎo) 1 NAMD是如何挑選合適的參數(shù)進(jìn)行對(duì)勢計(jì)算的 2 NAMD是否對(duì)運(yùn)算過程進(jìn)行過優(yōu)化或近似?難道它就窮舉所有原子對(duì)勢? 3 有沒有哪本書,或者資料對(duì)NAMD勢能計(jì)算過程進(jìn)行過講解? 求指點(diǎn) |
木蟲 (正式寫手)
|
1、參數(shù)自然是包括在CHARMM力場當(dāng)中的,你要先知道CHARMM力場的格式 2、NAMD只計(jì)算一定范圍內(nèi)的原子(如10A),長程靜電作用可以用一定的方法近似(如PME) 3、基本上任何關(guān)于分子模擬的書都會(huì)花幾頁講一下力場的。。。比如陳正隆的那本 |
木蟲 (正式寫手)
木蟲 (正式寫手)
金蟲 (小有名氣)
|
通過你給的原子類型定義 因?yàn)榱鍪峭ㄟ^定義一個(gè)原子類型 然后擬合其相關(guān)的各種參數(shù)來實(shí)現(xiàn)計(jì)算的 很簡單 就類似于一個(gè)表格 你在算一對(duì)兒原子作用的時(shí)候 分別查看他們的原子類型和他們的成鍵關(guān)系 然后加載力場參數(shù)去帶公式計(jì)算就好了 幾乎所有動(dòng)力學(xué)軟件都用拓?fù)湮募䜩肀4嬖宇愋秃统涉I關(guān)系的 軟件實(shí)現(xiàn)的細(xì)節(jié)會(huì)遠(yuǎn)比你想象的復(fù)雜 因?yàn)槠陂g有各種不同級(jí)別的算法的存在 包括并行化的代碼 直接看會(huì)不大容易理清其中的結(jié)構(gòu)的 |
金蟲 (小有名氣)
| NAMD的全名如果知道的話 你就不會(huì)問第二個(gè)問題了 對(duì)于求解大規(guī)模的動(dòng)力學(xué)模擬 優(yōu)化是必須的 而且程度非常高 最起碼并行化就需要調(diào)整通訊和計(jì)算 任務(wù)分配 調(diào)度匯總等等問題 而近似是動(dòng)力學(xué)天生的特性 在盡量保證精度和速度下 其實(shí)做了非常非常多的近似和簡化處理 要理解這個(gè)就去多看看力場和算法的書吧 保證顛覆你的現(xiàn)有完美主義想法 |
木蟲 (正式寫手)
木蟲 (正式寫手)
新蟲 (初入文壇)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|