LAMMPS降(升)温的几个思路

LAMMPS模拟降温有两种方法:连续降温与准静态模拟。连续降温通过设置`fix nvt temp`实现,而准静态模拟则通过循环升温或降温并确保长时间平衡。在循环中,首先在最低温度下进行平衡,然后通过`npt`修正进行温度变化,并记录数据。此过程需确保足够的模拟步数以达到稳定。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

模拟降温,有大致两种办法,一个是连续(升)降温,也就是dt间隔下不停的温度下降(升);第二种方法是准静态模拟降温,每降低一个温度,跑很长时间;这两种方法大家可以根据自己情况进行模拟。

对于第一种方法,只需要进行fix nvt temp 1000 100 …的设置就行了,降温的快慢,由你run N 的N来决定。

对于第二种方法提供如下思路,也就是使用循环的方法(升温,降温改小参数,+改成-就行啦):

1.先在最低温度下,跑平衡输出restart.equ低温稳定的结构

units metal
boundary p p p
atom_style atomic

neighbor 0.3 bin
neigh_modify delay 10

lattice bcc 2.87
region box block 0 10 0 10 0 10
create_box 1 box
create_atoms 1 box

pair_style eam/fs
pair_coeff * * fe.eam Fe
mass 1 56
fix 1 all npt temp 373.0 373.0 100.0 iso 1.0 1.0 1000.0thermo 10
timestep 0.001

run 1000000000
write_restart restart.equ

  1. 用循环的方法进行准平衡升温抽样

units metal
boundary p p p
atom_style atomic

label loopa
variable i loop 25
variable t equal 273.0+100.0*$i
variable v equal vol

read_restart restart.equ

pair_style eam/fs
pair_coeff * * fe.eam Fe

fix 1 all npt temp $t $t 100.0 iso 1.0 1.0 1000.0

thermo 10
timestep 0.001

run 10000000(足够长)
write_restart restart.equ

print “$t $v” append gmm.txt

clear

loop

next i

jump in.script loopa

例子只提供思路,不保证每个命令的语法正确

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值