24小时热门版块排行榜    

查看: 815  |  回复: 0

影月beatrice

铁虫 (初入文坛)

[求助] 关于锆合金到辐照计算

我想对锆合金辐照进行模拟,用 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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 影月beatrice 的主题更新
信息提示
请填处理意见