| 查看: 1959 | 回复: 3 | |||
[交流]
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} # 产出步数 # 跑产出阶段 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
2026国自然函评费到账
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
2楼2020-05-12 18:06:31
lqpzh 
新虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 111.6
- 散金: 5
- 帖子: 90
- 在线: 5.8小时
- 虫号: 9150633
- 注册: 2018-06-21
- 专业: 金属材料的合金相、相变及
3楼2020-05-13 07:15:26
4楼2020-05-18 08:52:22
40