24小时热门版块排行榜

南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知

返回列表

/name-?=p/

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 186.9
红花: 1
帖子: 49
在线: 9.4小时
虫号: 7570339
注册: 2017-12-05
专业: 物理学I

[交流] lammps在npt平衡后换成nve温度一直升高是怎么回事？已有2人参与

各位大神请教一个问题
用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} # 产出步数
# 跑产出阶段

回复此楼