| 查看: 1523 | 回复: 1 | |||
| 【悬赏金币】回答本帖问题，作者六氟磷酸锂将赠送您 5 个金币 | |||
[求助]
Non-numeric positions报错怎么解决
|
|||
用lammps做al/rdx体系的恒温热浴模拟，力场选择为wang.et al等人2017修正的reaxff-lg力场[j. phys. chem. c 2017, 121, 14597−14610]，在做3500 k恒温热浴计算时大概20 ps后出现non-numeric positions-simulation unstable错误，具体in文件如下，请问各位虫友怎么解决这个问题？ 力场.jpg |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : in
2021-02-04 11:19:27, 1.5 K
» 猜你喜欢
085701环境工程专硕求调剂
已经有6人回复
-
0817化学工程319求调剂
已经有10人回复
-
伙伴们，祝我生日快乐吧
已经有5人回复
-
274求调剂0856材料化工
已经有10人回复
-
337一志愿华南理工0805材料求调剂
已经有6人回复
-
（081700）化学工程与技术-298分求调剂
已经有7人回复
-
材料与化工304求B区调剂
已经有5人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有9人回复
-
一志愿北京化工大学材料与化工296分求调剂
已经有17人回复
-
0703化学一志愿211 总分320求调剂
已经有3人回复
|
in 文件在这里 ######################################################################## # Energy in eV, time in fs, distance in angstroms. units real boundary p p p atom_style charge read_data RDX-Al.data variable pe equal pe variable ke equal ke variable p equal (pxx+pyy+pzz)/3 variable temp equal temp variable step equal step compute ke all ke/atom # Create region for wall. region RDX-Al-Sphere sphere 125 125 125 60 units box region O2-Wall sphere 125 125 125 62.5 side out units box #Define force field pair_style reax/c lmp_control checkqeq yes safezone 3 mincap 600 pair_coeff * * ffield.reax.CHONAl C H O N Al #CHON mass 1 12.0 mass 2 1.008 mass 3 15.999 mass 4 14.000 mass 5 26.982 neighbor 2 bin neigh_modify every 1 delay 0 check no fix 1 all qeq/reax 1 0.0 10.0 1e-6 reax/c ###### Static relaxation at 0 K ######### timestep 0.1 thermo 100 thermo_style custom step temp lx ly lz vol press pe ke etotal thermo_modify lost ignore min_style cg minimize 1.0e-6 1.0e-8 5000 100000 ######## Establish wall so that RDX-Al will not interact with O2 ########### fix wall all wall/region RDX-Al-Sphere lj126 1.0 1.0 2.5 ####### Dynamic relaxation at 1500 K ########### reset_timestep 0 restart 10000 *.restart fix bond all reax/c/bonds 1000 Decomposition.reax dump 1 all custom 1000 Decomposition.lammpstrj id type x y z q fix plot1 all ave/time 1 1000 1000 v_temp v_pe v_ke v_p file Decomposition.txt velocity all create 3500.0 9999 fix nvt all nvt temp 3500.0 3500.0 10.0 run 100000 |
2楼2021-02-04 11:28:39
