| 查看: 1140 | 回复: 2 | ||
[求助]
不知道为什么错，求指导
|
|
LAMMPS (1 Feb 2014) ### REAX potential for GO ### This script is prepared on Feb. 9, 2011 ### QPHLL ## Initialization units real dimension 3 boundary p p p atom_style charge ## Atom definition read_data data.CHO orthogonal box = (2.1226 -0.257523 -0.855987) to (77.1226 74.7425 24.704) 4 by 4 by 1 MPI processor grid reading atoms ... 3078 atoms #Force Field, in ReaxFF, 1-C, 2-H, 3-O, 4-N, 5-S pair_style reax pair_coeff * * ffield.reax 1 2 3 # Settings # mass 1 12.0107 # mass 2 1.00794 # mass 3 15.9994 neighbor 2.5 bin neigh_modify delay 0 every 2 check yes # TimeUnit=fs timestep 0.25 compute reax all pair reax variable eb equal c_reax[1] variable ea equal c_reax[2] variable elp equal c_reax[3] variable emol equal c_reax[4] variable ev equal c_reax[5] variable epen equal c_reax[6] variable ecoa equal c_reax[7] variable ehb equal c_reax[8] variable et equal c_reax[9] variable eco equal c_reax[10] variable ew equal c_reax[11] variable ep equal c_reax[12] variable efi equal c_reax[13] variable eqeq equal c_reax[14] ## Output frequency thermo 10 thermo_style custom step cpu press vol temp etotal pe ke v_eb v_ea v_elp v_emol v_ev v_epen v_ecoa v_ehb v_et v_eco v_ew v_ep v_efi v_eqeq thermo_modify line one restart 10000 restart.* # AtomEye, scaled to box size, xs, ys, zs will be [0,1] dump 1 all cfg 1000 dump.*.cfg mass type xs ys zs vx vy vz fx fy fz #dump_modify 1 element C H O dump_modify 1 sort id # VMD dump 2 all dcd 100 movie.dcd dump_modify 2 sort id # XYZ dump 3 all xyz 10000 structure.*.xyz dump_modify 3 sort id # Similar to ReaxFF_xmolout, further modification needed! dump 4 all custom 10 xmolout type x y z dump_modify 4 sort id # Charge information dump 5 all custom 10 charge.q id type q dump_modify 5 sort id # Charge informaiton with coordinates dump 6 all custom 10 coordinate.q id type x y z q dump_modify 6 sort id minimize 1.0e-6 1.0e-6 100000 100000 WARNING: Resetting reneighboring criteria during minimization (../min.cpp:173) |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 速求各位大神指点，小弟感激不尽！
已经有4人回复
- Matlab用fsolve求解非线性方程，卡了三四天了跪求帮助！
已经有35人回复
- Materials Express投稿问题，急求...
已经有24人回复
- 用Matlab求一个二阶微分方程，出错了，求指导（求助）
已经有13人回复
- comsol multiphysics中导入AutoCAD模型出现错误
已经有12人回复
- matlab求模型参数有错误，求指导
已经有5人回复
- 求助大神！DNA-PAGE一直出现问题，不知道什么原因，请求大神回复帮忙解决。
已经有15人回复
- 审稿人提了一个ORR分析的问题，不知道什么意思，求指导！解释清楚会再加分！！
已经有12人回复
- 现阶段不知道该做什么准备，求指导
已经有3人回复
- 导师的回复让我不知道后续要怎么办求指导？?急！！
已经有12人回复
- 我想做3D图，但用好几版ORIGIN都没做出来，尤其是转矩阵都不行
已经有11人回复
- 初学TD-DFT计算一直出错，求解
已经有8人回复
- 6.22考试，啥都不知道，开始准备没几天的求指导
已经有13人回复
- MATLAB初学者，这样的式子怎么就是画不了图，总是出错!求指导!!!!!!!!!!!!!
已经有8人回复
- 谈妥了了一个老板，结果发现才是assistant professor,不知道怎么样？？？求指导
已经有7人回复
- 【求助】BSSE of O2运行出错
已经有17人回复
- 【求助】formchk出错
已经有23人回复
- 【求助】MS测试时出错
已经有4人回复
差点儿是帅哥
金虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 1943.8
- 散金: 300
- 红花: 22
- 帖子: 1356
- 在线: 168.6小时
- 虫号: 2641553
- 注册: 2013-09-08
- 专业: 零件加工制造
2楼2015-05-14 10:38:17
3楼2017-12-25 20:49:48
5