 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 2058 | 回復: 6 | ||
[求助]
gromacs做全原子模擬時候，突然到MD這一步出錯
|
 MD的控制文件是這樣的：title = MD simulation ;Preprocessor cpp = /lib/cpp ;Directories to include in the topology format include = -I../top ;Run control: A leap-frog algorithm for integrating Newton's equations. integrator = md ;Total simulation time: 1000 ps :time step in femtoseconds dt = 0.002 ;number of steps nsteps = 10000000 ;frequency to write coordinates to output trajectory file nstxout = 5000 ;frequency to write velocities to output trajectory file nstvout = 5000 ;frequency to write energies to log file nstlog = 1000 ;frequency to write energies to energy file nstenergy = 1000 ;frequency to write coordinates to xtc trajectory nstxtcout = 1000 ;group(s) to write to xtc trajectory xtc_grps = OCT ;group(s) to write to energy file energygrps = OCT ;Frequency to update the neighbor list (and the long-range forces, ;when using twin-range cut-off's). nstlist = 5 ;Make a grid in the box and only check atoms in neighboring grid cells ;when constructing a new neighbor list every nstlist steps. ns_type = grid ;cut-off distance for the short-range neighbor list rlist = 1.0 ;treatment of electrostatic interactions coulombtype = PME rcoulomb = 1.0 fourierspacing =0.12 pme_order =4 ewald_rtol =1e-5 optimize_fft =yes ;treatment of van der waals interactions rvdw = 1.4 ; Periodic boudary conditions in all the directions pbc = xyz ;Temperature coupling tcoupl = berendsen tc-grps = OCT tau_t = 0.1 ref_t = 290 ;Pressure coupling ;Pcoupl = berendsen ;Pcoupltype = isotropic ;tau_p = 1.0 ;compressibility = 4.5e-5 ;ref_p = 1.0 ;Velocity generation gen_vel = yes gen_temp = 290 gen_seed = 173529 ; Type of constraint algorithm constraint-algorithm = Lincs lincs-iter = 4 ;Constrain bonds constraints = none 出現(xiàn)的錯誤是： creating statusfile for 1 node... ' for variable integrator, using 'md' Next time use one of: 'md' 'steep' 'cg' 'bd' 'sd' 'nm' 'l-bfgs' 'tpi' ' for variable ns-type, using 'Grid' Next time use one of: 'Grid' 'Simple' ' for variable pbc, using 'xyz' Next time use one of: 'xyz' 'no' 'full' ' for variable coulombtype, using 'Cut-off' Next time use one of: 'Cut-off' 'Reaction-Field' 'Generalized-Reaction-Field' 'PME' 'Ewald' 'PPPM' 'Poisson' 'Switch' 'Shift' 'User' 'Generalized-Born' 'Reaction-Field-nec' 'Encad-shift' 'PME-User' ' for variable optimize_fft, using 'no' Next time use one of: 'no' 'yes' ' for variable tcoupl, using 'No' Next time use one of: 'No' 'Berendsen' 'Nose-Hoover' 'yes' 'Andersen' 'Andersen-interval' ' for variable gen-vel, using 'no' Next time use one of: 'no' 'yes' ' for variable constraints, using 'none' Next time use one of: 'none' 'h-bonds' 'all-bonds' 'h-angles' 'all-angles' ' for variable constraint-algorithm, using 'Lincs' Next time use one of: 'Lincs' 'Shake' Fatal error: number of coordinates in coordinate file (box.pdb, 3900) does not match topology (octanebox.top, 0) top文件在能量最小化的那一步?jīng)]有任何問題，但是在MD這一步就不行了。 |
回復此樓» 猜你喜歡
0703化學336分求調(diào)劑
已經(jīng)有6人回復
-
304求調(diào)劑
已經(jīng)有9人回復
-
293求調(diào)劑
已經(jīng)有13人回復
-
281求調(diào)劑（0805）
已經(jīng)有8人回復
-
環(huán)境領域全國重點實驗室招收博士1-2名
已經(jīng)有3人回復
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復
-
301求調(diào)劑
已經(jīng)有6人回復
-
一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復
-
302求調(diào)劑
已經(jīng)有6人回復
-
26博士申請
已經(jīng)有3人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 在gromacs中運行mdrun命令時出錯
已經(jīng)有10人回復
- Gromacs中加入離子
已經(jīng)有5人回復
- 【轉(zhuǎn)帖】GROMACS分析程序討論，在gromacs中添加Amber力場
已經(jīng)有8人回復
- 【求助】gromacs-4 mdrun死掉
已經(jīng)有20人回復
- 【求助】這是我安裝gromacs過程，但是出錯阿，這個怎么辦呢？謝謝
已經(jīng)有11人回復
- 【求助】gromacs的又一個問題
已經(jīng)有22人回復
- 【求助】GROMACS分析二面角隨時間的變化
已經(jīng)有4人回復
- 【原創(chuàng)】在gromacs的gro中插入小分子
已經(jīng)有12人回復
- 【求助】gromacs出現(xiàn)了一個問題
已經(jīng)有10人回復
2樓2011-04-19 07:14:17
|
; ; File 'octanebox.top' was generated ; By user: public (501) ; On host: cluster.hpc.org ; At date: Thu Apr 14 17:21:22 2011 ; ; This is your topology file ; PRODRG COORDS ; ; Include forcefield parameters #include "ffoplsaa.itp" ;include needed parameters #include "octane.itp" ; Include surfacetant topology #include "sds.itp" ;include water parameters #include "spce.itp" #ifdef POSRES_WATER ; Position restraint for each water oxygen [ position_restraints ] ; i funct fcx fcy fcz 1 1 1000 1000 1000 #endif ; Include generic topology for ions #include "ions.itp" [ system ] ; Name octane box [ molecules ] ; Compound #mols OCT 150 |
3樓2011-04-19 13:56:05
|
; ; .itp file for octane in OPLS force field ; Generated by HYF ; ; [ moleculetype ] ; Name nrexcl OCT 3 [ atoms ] ; nr type resnr resid atom cgnr charge mass 1 opls_135 1 OCT C 1 -0.18 12.011 2 opls_140 1 OCT H 1 0.06 1.008 3 opls_140 1 OCT H 1 0.06 1.008 4 opls_140 1 OCT H 1 0.06 1.008 5 opls_136 1 OCT C 2 -0.12 12.011 6 opls_140 1 OCT H 2 0.06 1.008 7 opls_140 1 OCT H 2 0.06 1.008 8 opls_136 1 OCT C 3 -0.12 12.011 9 opls_140 1 OCT H 3 0.06 1.008 10 opls_140 1 OCT H 3 0.06 1.008 11 opls_136 1 OCT C 4 -0.12 12.011 12 opls_140 1 OCT H 4 0.06 1.008 13 opls_140 1 OCT H 4 0.06 1.008 14 opls_136 1 OCT C 5 -0.12 12.011 15 opls_140 1 OCT H 5 0.06 1.008 16 opls_140 1 OCT H 5 0.06 1.008 17 opls_136 1 OCT C 6 -0.12 12.011 18 opls_140 1 OCT H 6 0.06 1.008 19 opls_140 1 OCT H 6 0.06 1.008 20 opls_136 1 OCT C 7 -0.12 12.011 21 opls_140 1 OCT H 7 0.06 1.008 22 opls_140 1 OCT H 7 0.06 1.008 23 opls_135 1 OCT C 8 -0.18 12.011 24 opls_140 1 OCT H 8 0.06 1.008 25 opls_140 1 OCT H 8 0.06 1.008 26 opls_140 1 OCT H 8 0.06 1.008 [ bonds ] ; ai aj fu c0, c1, ... 1 2 1 1 3 1 1 4 1 1 5 1 5 6 1 5 7 1 5 8 1 8 9 1 8 10 1 8 11 1 11 12 1 11 13 1 11 14 1 14 15 1 14 16 1 14 17 1 17 18 1 17 19 1 17 20 1 20 21 1 20 22 1 20 23 1 23 24 1 23 25 1 23 26 1 [ pairs ] ; ai aj fu c0, c1, ... 2 6 1 2 7 1 2 8 1 3 6 1 3 7 1 3 8 1 4 6 1 4 7 1 4 8 1 1 9 1 1 10 1 1 11 1 6 9 1 6 10 1 6 11 1 7 9 1 7 10 1 7 11 1 5 12 1 5 13 1 5 14 1 9 12 1 9 13 1 9 14 1 10 12 1 10 13 1 10 14 1 8 15 1 8 16 1 8 17 1 12 15 1 12 16 1 12 17 1 13 15 1 13 16 1 13 17 1 11 18 1 11 19 1 11 20 1 15 18 1 15 19 1 15 20 1 16 18 1 16 19 1 16 20 1 14 21 1 14 22 1 14 23 1 18 21 1 18 22 1 18 23 1 19 21 1 19 22 1 19 23 1 17 24 1 17 25 1 17 26 1 21 24 1 21 25 1 21 26 1 22 24 1 22 25 1 22 26 1 [ angles ] ; ai aj ak fu c0, c1, ... 2 1 3 1 2 1 4 1 2 1 5 1 3 1 4 1 3 1 5 1 4 1 5 1 1 5 6 1 1 5 7 1 1 5 8 1 6 5 7 1 6 5 8 1 7 5 8 1 5 8 9 1 5 8 10 1 5 8 11 1 9 8 10 1 9 8 11 1 10 8 11 1 8 11 12 1 8 11 13 1 8 11 14 1 12 11 13 1 12 11 14 1 13 11 14 1 11 14 15 1 11 14 16 1 11 14 17 1 15 14 16 1 15 14 17 1 16 14 17 1 14 17 18 1 14 17 19 1 14 17 20 1 18 17 19 1 18 17 20 1 19 17 20 1 17 20 21 1 17 20 22 1 17 20 23 1 21 20 22 1 21 20 23 1 22 20 23 1 20 23 24 1 20 23 25 1 20 23 26 1 24 23 25 1 24 23 26 1 25 23 26 1 [ dihedrals ] ; ai aj ak al fu c0, c1, m, ... 2 1 5 6 3 2 1 5 7 3 2 1 5 8 3 3 1 5 6 3 3 1 5 7 3 3 1 5 8 3 4 1 5 6 3 4 1 5 7 3 4 1 5 8 3 1 5 8 9 3 1 5 8 10 3 1 5 8 11 3 6 5 8 9 3 6 5 8 10 3 6 5 8 11 3 7 5 8 9 3 7 5 8 10 3 7 5 8 11 3 5 8 11 12 3 5 8 11 13 3 5 8 11 14 3 9 8 11 12 3 9 8 11 13 3 9 8 11 14 3 10 8 11 12 3 10 8 11 13 3 10 8 11 14 3 8 11 14 15 3 8 11 14 16 3 8 11 14 17 3 12 11 14 15 3 12 11 14 16 3 12 11 14 17 3 13 11 14 15 3 13 11 14 16 3 13 11 14 17 3 11 14 17 18 3 11 14 17 19 3 11 14 17 20 3 15 14 17 18 3 15 14 17 19 3 15 14 17 20 3 16 14 17 18 3 16 14 17 19 3 16 14 17 20 3 14 17 20 21 3 14 17 20 22 3 14 17 20 23 3 18 17 20 21 3 18 17 20 22 3 18 17 20 23 3 19 17 20 21 3 19 17 20 22 3 19 17 20 23 3 17 20 23 24 3 17 20 23 25 3 17 20 23 26 3 21 20 23 24 3 21 20 23 25 3 21 20 23 26 3 22 20 23 24 3 22 20 23 25 3 22 20 23 26 3 |
4樓2011-04-19 13:56:30
5樓2012-12-06 22:47:30
明月照清渠
金蟲 (小有名氣)
- 應助: 10 (幼兒園)
- 金幣: 3550.4
- 散金: 1
- 紅花: 3
- 帖子: 95
- 在線: 69.3小時
- 蟲號: 1101813
- 注冊: 2010-09-17
- 性別: MM
- 專業(yè): 理論和計算化學
6樓2012-12-07 10:15:02
yuanbei1987
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 96.5
- 帖子: 12
- 在線: 22.9小時
- 蟲號: 1484155
- 注冊: 2011-11-09
- 性別: MM
- 專業(yè): 蛋白質(zhì)組學
7樓2012-12-07 15:24:11
相關版塊跳轉(zhuǎn)
我要訂閱樓主 zyj8119 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學336分求調(diào)劑 +5 | zbzihdhd 2026-03-15 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +8 | zjl的號 2026-03-16 | 13/650 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +5 | 慕寒mio 2026-03-16 | 5/250 |
|
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
|
[考研] 化學工程321分求調(diào)劑 +11 | 大米飯！ 2026-03-15 | 14/700 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +4 | 梨花珞晚風 2026-03-17 | 4/200 |
|
|
|
[考研] 一志愿南京大學，080500材料科學與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[碩博家園] 深圳大學碩士招生（2026秋，傳感器方向，僅錄取第一志愿） +4 | xujiaoszu 2026-03-11 | 9/450 |
|
|
|
[考研] 285化工學碩求調(diào)劑（081700） +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
|
[考研] 080500，材料學碩302分求調(diào)劑學校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 0703化學調(diào)劑 +4 | 快樂的香蕉 2026-03-11 | 4/200 |
|
|
|
[考研] 四川大學085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +5 | 一定有學上- 2026-03-12 | 5/250 |
|
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|