| 查看: 855 | 回复: 0 | |||
[求助]
关于锆合金到辐照计算
|
|
我想对锆合金辐照进行模拟，用 fix heat 提高一个原子到能量去撞击其他原子，但是输出显示ERROR: Fix heat kinetic energy went negative ，对加热却使动能变为负值表示不理解，请各位大神教教如何处理这个问题。 units metal boundary p p p atom_style atomic atom_modify map array read_data data.Zr group 1 type 1 pair_style eam/fs pair_coeff * * Zr_mm.eam.fs Zr neighbor 2.0 bin neigh_modify delay 10 check yes velocity 1 create 300.0 1234 #fix 1 all nvt temp 11385 11385 100 fix 1 all nve timestep 0.001 thermo 100 thermo_style custom temp pe ke run 100 unfix 1 group 2 id 1 fix 3 2 heat 1 10.0 #run 10 #fix 4 2 heat 10 -10.0 #fix 5 all nve #dump 1 all atom 1 Zr.data #dump 2 all atom 1 Zr.lammpstrj thermo 100 thermo_style custom temp pe ke run 100 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有53人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
30