| 查看: 3678 | 回复: 9 | |||
[求助]
gromacs运行问题 已有2人参与
|
|
各位大神，我用一个小程序对膜进行剪切处理后，进行能量最小化时出现下列问题： 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 这是嘛回事呢，咋解决？ |
回复此楼» 猜你喜欢
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有181人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有17人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
此贴已删除
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- abinit软件安装步骤（亲测）
已经有16人回复
- GROMACS运行出错了，求助中
已经有15人回复
- gromacs NPT restrain 两个分子
已经有5人回复
- Gromacs能量最小化问题
已经有5人回复
- autodock与gromacs在Windows下的安装
已经有10人回复
- GROMACS运行完后的结果分析
已经有4人回复
- compute 命令出错，急求
已经有5人回复
- 这个错误警告怎么回事啊？不能运行了
已经有5人回复
- ubuntu运行vasp5.3.3出错
已经有6人回复
- gromacs编译安装
已经有11人回复
- Gromacs运行时itp文件语法出错
已经有5人回复
- gromacs安装出错求助
已经有4人回复
- Gromacs动力学运行出错
已经有4人回复
- grompp无法识别原子类型
已经有7人回复
- 请教GROMACS安装
已经有8人回复
- gromacs运行出错，求方案，谢谢！！！
已经有6人回复
- 急求fftw3.3.1安装问题，解答
已经有13人回复
- 在gromacs中运行mdrun命令时出错
已经有10人回复
- 求Gromacs中Glycam力场数据包。。非常感谢
已经有5人回复
- gromacs运行出错
已经有3人回复
- 【求助】安装后修改.bash_profile的问题
已经有5人回复
- 【求助】这是我安装gromacs过程，但是出错阿，这个怎么办呢？谢谢
已经有11人回复
- 【求助】gromacs的又一个问题
已经有22人回复
- 【求助】Gromacs运行出错
已经有4人回复
- 【求助】MPI2-1.0.8安装出错求助
已经有5人回复
softmatter
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 4439.2
- 帖子: 723
- 在线: 73.5小时
- 虫号: 298219
- 注册: 2006-11-19
- 专业: 凝聚态物性I:结构、力学和
2楼2014-04-12 15:35:51
3楼2014-04-12 23:16:02
softmatter
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 4439.2
- 帖子: 723
- 在线: 73.5小时
- 虫号: 298219
- 注册: 2006-11-19
- 专业: 凝聚态物性I:结构、力学和
4楼2014-04-12 23:51:27
★
fegg7502: 金币+1, 鼓励交流 2014-05-05 08:22:47
fegg7502: 金币+1, 鼓励交流 2014-05-05 08:22:47
|
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文件，这是神马节奏？ |
5楼2014-04-13 00:16:30
softmatter
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 4439.2
- 帖子: 723
- 在线: 73.5小时
- 虫号: 298219
- 注册: 2006-11-19
- 专业: 凝聚态物性I:结构、力学和
6楼2014-04-13 08:48:44
7楼2014-04-13 09:48:31
youngfi
木虫 (小有名气)
- 应助: 9 (幼儿园)
- 金币: 2419.9
- 散金: 15
- 红花: 1
- 帖子: 144
- 在线: 100.9小时
- 虫号: 1850824
- 注册: 2012-06-07
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
8楼2014-04-28 08:19:07
9楼2014-05-04 23:29:30
10楼2017-01-06 12:40:58
50