| 5 | 1/1 | 返回列表 |
| 查看: 3634 | 回復: 9 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
彬哥木蟲 (小有名氣)
|
[求助]
gromacs運行問題 已有2人參與
|
||
|
各位大神,我用一個小程序對膜進行剪切處理后,進行能量最小化時出現(xiàn)下列問題: Steepest Descents: Tolerance (Fmax) = 1.00000e+03 Number of steps = 500000 Step Time Lambda 0 0.00000 0.00000 DD step 0 load imb.: force 12.9% Step Time Lambda 1 1.00000 0.00000 Step Time Lambda 2 2.00000 0.00000 Step Time Lambda 3 3.00000 0.00000 Step Time Lambda 4 4.00000 0.00000 Step Time Lambda 5 5.00000 0.00000 Step Time Lambda 6 6.00000 0.00000 Step Time Lambda 7 7.00000 0.00000 Step Time Lambda 8 8.00000 0.00000 Step Time Lambda 9 9.00000 0.00000 Step Time Lambda 10 10.00000 0.00000 Step Time Lambda 11 11.00000 0.00000 Step Time Lambda 12 12.00000 0.00000 Step Time Lambda 13 13.00000 0.00000 Step Time Lambda 14 14.00000 0.00000 Stepsize too small, or no change in energy. Converged to machine precision, but not to the requested precision Fmax < 1000 Double precision normally gives you higher accuracy. Steepest Descents converged to machine precision in 15 steps, but did not reach the requested Fmax < 1000. Potential Energy = 1.5981487e+19 Maximum force = inf on atom 5867 Norm of force = inf 這是嘛回事呢,咋解決? |
木蟲 (小有名氣)
|
Initializing Domain Decomposition on 8 nodes Dynamic load balancing: auto Will sort the charge groups at every domain (re)decomposition Initial maximum inter charge-group distances: two-body bonded interactions: 0.661 nm, LJ-14, atoms 4431 4452 multi-body bonded interactions: 0.661 nm, Proper Dih., atoms 4431 4452 Minimum cell size due to bonded interactions: 0.727 nm Maximum distance for 5 constraints, at 120 deg. angles, all-trans: 0.771 nm Estimated maximum distance required for P-LINCS: 0.771 nm This distance will limit the DD cell size, you can override this with -rcon Using 0 separate PME nodes Scaling the initial minimum size with 1/0.8 (option -dds) = 1.25 Optimizing the DD grid for 8 cells with a minimum initial size of 0.964 nm The maximum allowed number of cells is: X 11 Y 11 Z 6 Domain decomposition grid 8 x 1 x 1, separate PME nodes 0 PME domain decomposition: 8 x 1 x 1 Domain decomposition nodeid 0, coordinates 0 0 0 Table routines are used for coulomb: TRUE Table routines are used for vdw: FALSE Will do PME sum in reciprocal space. ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++ U. Essman, L. Perela, M. L. Berkowitz, T. Darden, H. Lee and L. G. Pedersen A smooth particle mesh Ewald method J. Chem. Phys. 103 (1995) pp. 8577-8592 -------- -------- --- Thank You --- -------- -------- Will do ordinary reciprocal space Ewald sum. Using a Gaussian width (1/beta) of 0.384195 nm for Ewald Cut-off's: NS: 1.2 Coulomb: 1.2 LJ: 1.2 Long Range LJ corr.: <C6> 6.3119e-03 System total charge: -0.000 Generated table with 1100 data points for Ewald. Tabscale = 500 points/nm Generated table with 1100 data points for LJ6. Tabscale = 500 points/nm Generated table with 1100 data points for LJ12. Tabscale = 500 points/nm Generated table with 1100 data points for 1-4 COUL. Tabscale = 500 points/nm Generated table with 1100 data points for 1-4 LJ6. Tabscale = 500 points/nm Generated table with 1100 data points for 1-4 LJ12. Tabscale = 500 points/nm Configuring nonbonded kernels... Configuring standard C nonbonded kernels... Testing x86_64 SSE2 support... present. Removing pbc first time Initializing Parallel LINear Constraint Solver ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++ B. Hess P-LINCS: A Parallel Linear Constraint Solver for molecular simulation J. Chem. Theory Comput. 4 (2008) pp. 116-122 -------- -------- --- Thank You --- -------- -------- The number of constraints is 18564 There are inter charge-group constraints, will communicate selected coordinates each lincs iteration Linking all bonded interactions to atoms There are 45500 inter charge-group exclusions, will use an extra communication step for exclusion forces for PME The initial number of communication pulses is: X 1 The initial domain decomposition cell size is: X 1.44 nm The maximum allowed distance for charge groups involved in interactions is: non-bonded interactions 1.200 nm (the following are initial values, they could change due to box deformation) two-body bonded interactions (-rdd) 1.200 nm multi-body bonded interactions (-rdd) 1.200 nm atoms separated by up to 5 constraints (-rcon) 1.438 nm When dynamic load balancing gets turned on, these settings will change to: The maximum number of communication pulses is: X 2 The minimum size for domain decomposition cells is 1.049 nm The requested allowed shrink of DD cells (option -dds) is: 0.80 The allowed shrink of domain decomposition cells is: X 0.73 The maximum allowed distance for charge groups involved in interactions is: non-bonded interactions 1.200 nm two-body bonded interactions (-rdd) 1.200 nm multi-body bonded interactions (-rdd) 1.049 nm atoms separated by up to 5 constraints (-rcon) 1.049 nm Making 1D domain decomposition grid 8 x 1 x 1, home cell index 0 0 0 Center of mass motion removal mode is Linear We have the following groups for center of mass motion removal: 0: POPC There are: 18928 Atoms Charge group distribution at step 0: 1547 1628 1762 1544 1654 1623 1699 1647 Grid: 4 x 13 x 7 cells Constraining the starting coordinates (step 0) Constraining the coordinates at t0-dt (step 0) RMS relative constraint deviation after constraining: 1.41e-06 Initial temperature: 0.00294003 K Started mdrun on node 0 Sat Apr 12 23:30:51 2014 Step Time Lambda 0 0.00000 0.00000 Energies (kJ/mol) G96Angle Proper Dih. Improper Dih. LJ-14 Coulomb-14 1.79277e+04 2.19389e+04 1.21721e+03 7.28810e+02 3.52807e+04 LJ (SR) Disper. corr. Coulomb (SR) Coul. recip. Potential -7.88621e+04 -3.55517e+03 -7.88547e+04 -3.15743e+05 -3.99921e+05 Kinetic En. Total Energy Temperature Pres. DC (bar) Pressure (bar) 5.89721e+00 -3.99915e+05 3.71179e-02 -1.53150e+02 -5.35095e+02 Constr. rmsd 我也不清楚那里有問題,然后命令停止,輸出幾個被切割開的pdb文件,這是神馬節(jié)奏? |
木蟲 (正式寫手)
木蟲 (小有名氣)
木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 346求調劑 +4 | 鄭誠樂 2026-04-02 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿華東理工大學,080500學碩,317分,求調劑 +12 | s1145 2026-03-31 | 12/600 |
|
|
[考研] 0856材料與化工調劑,339 +13 | 10213207 2026-03-31 | 13/650 |
|
|
[考研] 348求調劑 +6 | 小懶蟲不懶了 2026-03-27 | 7/350 |
|
|
[考研] 321求調劑 +9 | y-yh 2026-04-01 | 10/500 |
|
|
[考研] 考研調劑0702 +3 | 梅子味晚霞 2026-03-28 | 3/150 |
|
|
[考研] 302求調劑一志愿北航070300,本科鄭大化學 +8 | 圣日耳曼條 2026-04-01 | 11/550 |
|
|
[考研] 0856,269分求調劑 +8 | 有學上就行求求?/a> 2026-03-30 | 11/550 |
|
|
[考研] 一志愿北交材料工程總分358 +5 | cs0106 2026-04-01 | 7/350 |
|
|
[考研] 326求調劑 +4 | 崽崽仔 2026-03-31 | 4/200 |
|
|
[考研] 080200學碩,機械工程專業(yè)277分,求帶走! +4 | 瓶子PZ 2026-03-31 | 4/200 |
|
|
[考研] 考研調劑求助 +7 | 13287130938 2026-03-31 | 7/350 |
|
|
[考研] 343求調劑 +8 | 愛羈絆 2026-03-28 | 8/400 |
|
|
[考研] 求收留 +8 | 1943443204 2026-03-28 | 8/400 |
|
|
[考研] 276求調劑 +3 | 趙久華 2026-03-29 | 3/150 |
|
|
[考研] 08工科求調劑286 +5 | tgs_001 2026-03-28 | 5/250 |
|
|
[考研] 081200-11408-276學碩求調劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
[考研] 308求調劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
[考研] 298調劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|