各位大神请教一个问题
用lammps模拟计算3C-SiC热导率,时间步长1fs  在1000K下，先在NVT下模拟10万步，使系统达到平衡然后换上NVE系综后在运行125万步。但是当换上NVE系综不久后，温度就开始升高，到结束时温度都达到几千甚至上万K，这么高温度碳化硅固体早就没了。125万步不是设定步数，125万步时原子丢失了，在低温时 温度也会升高  想问这是哪里出问题了
这是in文件   晶格常数设置的读坐标文件  盒子边长大小是十个晶格常数也就是一共8000个原子
variable    T     equal 1000   # 温度为 400 K
#variable A equal 4.348        # 晶格常数大概为 6.35A
variable    DT    equal 0.001  # 积分步长为 0.3 fs

# 初始化坐标、势函  数、速度等数据

boundary           p p p

read_data        SiC-beta.txt

pair_style tersoff

pair_coeff    * * SiC.tersoff C Si

min_style     cg                                 
#minimize      1.0e-9 1.0e-9 10000 30000   
velocity all create ${T} 123456  dist gaussian         # 根据温度初始化速度
# 在NVT系综下平衡体系
fix    NVT all nvt temp ${T} ${T} 0.1        # Nose-Hoover 热浴
timestep     ${DT}                        # 设置积分步长
thermo_style custom step temp press       # 观察温度和压强的变化
thermo       1000                         # 100 步输出一次

# 根据以上设置跑若干步，目的是使系统达到热力学平衡态
run         100000                       # 跑 10000 步

# 在产出阶段将用NVE系综
unfix NVT         # 撤掉之前的 NVT 系综
fix NVE all nve         # 换上 NVE 系综

# 在产出阶段计算热流
compute KE all ke/atom         # 计算单原子动能
compute PE all pe/atom         # 计算单原子势能
compute V all stress/atom NULL virial         # 计算单原子位力
compute J all heat/flux KE PE V         # 计算单原子热流

# 在产出阶段计算热流自关联函数 <J(0)J(t)>
variable Ns equal 100         # 取样间隔
variable Nc equal 100         # 关联数据量（最大关联时间为 Ns*Nc*DT）
variable Np equal ${Ns}*${Nc}*10         # 产出步数（这是一个好的选择）
fix         HAC all ave/correlate ${Ns} ${Nc} ${Np} c_J[1] c_J[2] c_J[3] &     
        type auto file hac.txt                                       #！！！！！！！！！！！！！！！   

variable M equal 14 # 热流自关联函数的个数
variable N equal ${M}*${Np} # 产出步数
# 跑产出阶段 返回小木虫查看更多