| 8 | 1/1 | 返回列表 |
| 查看: 2177 | 回復(fù): 7 | ||
| 本帖產(chǎn)生 1 個(gè) 模擬EPI ,點(diǎn)擊這里進(jìn)行查看 | ||
plafalines捐助貴賓 (小有名氣)
|
[求助]
用tinker的amoeba09力場(chǎng)模擬~~總是提示參數(shù)未定義~~~怎么辦啊~熱切求助
|
|
|
用tinker的amoeba09模擬一個(gè)NaCl水溶液體系 總是提示: Random Number Generator Initialized with SEED : 123456789 MECHANIC -- Some Required Potential Energy Parameters are Undefined TINKER is Unable to Continue; Terminating the Current Calculation xyz文件和key文件見附件 求助各位高手 [ Last edited by plafalines on 2011-5-3 at 22:55 ] |
木蟲 (小有名氣)
捐助貴賓 (小有名氣)
木蟲 (初入文壇)
|
你的key文件里只有力場(chǎng)參數(shù),沒有動(dòng)力學(xué)計(jì)算需要的參數(shù),當(dāng)然會(huì)提示有些參數(shù)沒有定義啦~ 這里有我測(cè)試用的key文件。你用時(shí)把parameters這行后面的改成你自己裝的Tinker中amoeba09力場(chǎng)參數(shù)文件的位置。其他參數(shù)你改成適合你體系的就行了。 parameters /opt/tinker/params/amoeba09 verbose save-velocity archive polar-eps 0.01 cutoff 9.0 neighbor-list list-buffer 0.5 ewald ewald-cutoff 9.0 pme-grid 19 19 19 pme-order 5 fft-package fftw a-axis 19.8647 integrate verlet thermostat berendsen tau-temperature 0.2 barostat berendsen tau-pressure 1.0 compress 0.000046 |

木蟲 (正式寫手)
|
你這key文件定義了PME,但是4.2版手冊(cè)67頁上說: EWALD SUMMATION METHODS TINKER contains a versions of the Ewald summation technique for inclusion of long range electrostatic interactions via periodic boundaries. The particle mesh Ewald (PME) method is available for simple charge-charge potentials, while regular Ewald is provided for polarizable atomic multipole interactions. 即TINKER中PME方法不用于極化力場(chǎng)的靜電計(jì)算,它只適于固定電荷模型。 |
木蟲 (初入文壇)
|
這個(gè)設(shè)置PME參數(shù),計(jì)算時(shí)程序也會(huì)根據(jù)體系實(shí)際的力場(chǎng)形式選擇計(jì)算方法的,如果對(duì)于多偶極勢(shì)采用Ewald算法,那設(shè)置PME參數(shù)的參數(shù)應(yīng)該會(huì)自動(dòng)忽略掉~ 另外我用的tinker版本是5.0的。我看了下源碼,mpole的Ewald算法,其虛空間也是用PME算法處理了,不知道是不是新版本改了~可以看下源碼里的 empole1.f 和 pmepol.f |

木蟲 (初入文壇)

木蟲 (正式寫手)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085601材料工程專碩求調(diào)劑 +5 | 慕寒mio 2026-03-16 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +4 | zbzihdhd 2026-03-15 | 5/250 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
[考研] 302求調(diào)劑 +9 | 負(fù)心者當(dāng)誅 2026-03-11 | 9/450 |
|
|
[考研] 211本,11408一志愿中科院277分,曾在中科院自動(dòng)化所實(shí)習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研]
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[基金申請(qǐng)]
有必要更換申報(bào)口嗎
20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
[考研] 材料工程調(diào)劑 +4 | 咪咪空空 2026-03-11 | 4/200 |
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|