| 查看: 1717 | 回复: 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
» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有10人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有11人回复
-
烦躁
已经有5人回复
-
寻找合作团队
已经有4人回复
|
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