| 查看: 833 | 回复: 0 | ||
[求助]
lammps 中的lj势
|
|
请问lj势用的时候可以不涉及到具体的原子么？两种类型的原子该怎么写呢？我先用了具体的原子这样写的 # LJ 12-6 system with npt ensemble units lj atom_style atomic boundary p p p lattice fcc 1.073 lattice fcc 0.8 region box block 0 10 0 5 0 5 create_box 2 box create_atoms 1 box create_atoms 2 box mass 1 1.0 mass 2 1.0 velocity all create 0.01 872877 timestep 0.01 dump 1 all xyz 1000 lj.xyz pair_style lj/cut 2.5 pair_coeff 1 1 1.0 1.0 2.5 pair_coeff 2 2 1.6 1.1765 2.5 pair_coeff 1 2 0.2 1.0882 2.5 #pair_modify shift yes mix geometric neighbor 0.3 bin neigh_modify every 10 delay 0 check yes thermo 1000 fix 1 all npt temp 0.01 0.01 1.00 iso 0.0024 0.0024 100.0 run 50000 unfix 1 thermo 1000 运行结果总出现nan怎么回事？不用具体的原子怎么写代码呢？ 新手跪谢！！！！ |
回复此楼» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有143人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有6人回复
-
散金币祈福
已经有32人回复
