| 查看: 776 | 回复: 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怎么回事？不用具体的原子怎么写代码呢？ 新手跪谢！！！！ |
回复此楼» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有202人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复