| 查看: 817 | 回复: 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怎么回事？不用具体的原子怎么写代码呢？ 新手跪谢！！！！ |
回复此楼» 猜你喜欢
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有6人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有257人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有24人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有22人回复