 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2060 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
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這一步就不行了。 |
回復(fù)此樓» 猜你喜歡
280求調(diào)劑
已經(jīng)有7人回復(fù)
-
環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗室招收博士1-2名
已經(jīng)有5人回復(fù)
-
303求調(diào)劑
已經(jīng)有6人回復(fù)
-
290求調(diào)劑
已經(jīng)有10人回復(fù)
-
0703化學(xué)求調(diào)劑 總分331
已經(jīng)有3人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有8人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有3人回復(fù)
-
材料，紡織，生物（0856、0710），化學(xué)招生啦
已經(jīng)有9人回復(fù)
-
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有3人回復(fù)
-
0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 在gromacs中運(yùn)行mdrun命令時出錯
已經(jīng)有10人回復(fù)
- Gromacs中加入離子
已經(jīng)有5人回復(fù)
- 【轉(zhuǎn)帖】GROMACS分析程序討論，在gromacs中添加Amber力場
已經(jīng)有8人回復(fù)
- 【求助】gromacs-4 mdrun死掉
已經(jīng)有20人回復(fù)
- 【求助】這是我安裝gromacs過程，但是出錯阿，這個怎么辦呢？謝謝
已經(jīng)有11人回復(fù)
- 【求助】gromacs的又一個問題
已經(jīng)有22人回復(fù)
- 【求助】GROMACS分析二面角隨時間的變化
已經(jīng)有4人回復(fù)
- 【原創(chuàng)】在gromacs的gro中插入小分子
已經(jīng)有12人回復(fù)
- 【求助】gromacs出現(xiàn)了一個問題
已經(jīng)有10人回復(fù)
|
; ; .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
yazhiyang
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 20.4
- 帖子: 13
- 在線: 21.2小時
- 蟲號: 423619
- 注冊: 2007-07-21
- 專業(yè): 化學(xué)
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
5樓2012-12-06 22:47:30
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +3 | 1孫悟空 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 297求調(diào)劑 +4 | 戲精丹丹丹 2026-03-17 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實(shí)用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研]
機(jī)械專碩調(diào)劑
+3
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
|
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +5 | ZY-05 2026-03-13 | 6/300 |
|
|
|
[考研] 333求調(diào)劑 +3 | 球球古力 2026-03-11 | 3/150 |
|
|
|
[考研] 材料工程調(diào)劑 +4 | 咪咪空空 2026-03-11 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|