Lammps将当前体系温度作为输入，循环升温

Mr. Material

已于 2023-03-05 11:03:08 修改

阅读量4.2k

点赞数 5

分类专栏： LAMMPS小技巧文章标签： lammps 循环升温

于 2021-11-11 16:53:23 首次发布

本文链接：https://blog.csdn.net/qq_43689832/article/details/121270798

版权

LAMMPS小技巧专栏收录该内容

32 篇文章

订阅专栏

关注 $\color{Violet} \rm Mr.material\ ,$ $\color{red}{更}$ $\color{blue}{多}$ $\color{orange}{精}$ $\color{green}{彩}$ ！

主要专栏内容包括：
$\textbf{ \underline{\dag《LAMMPS小技巧》：}}$ 主要介绍采用分子动力学（ $L amm p s$ ）模拟相关安装教程、原理以及模拟小技巧（难度： $\bigstar$ ）
$\textbf{ \underline{\dag\dag《LAMMPS实例教程—In文件详解》：}}$ 主要介绍采用分子动力学（ $L amm p s$ ）模拟相关物理过程模拟。（包含：热导率计算、定压比热容计算，难度： $\bigstar$ $\bigstar$ $\bigstar$ ）
$\textbf{ \underline{\dag\dag\dag《Lammps编程技巧及后处理程序技巧》：}}$ 主要介绍针对分子模拟的动力学过程（轨迹文件）进行后相关的处理分析（需要一定编程能力。难度： $\bigstar$ $\bigstar$ $\bigstar$ $\bigstar$ $\bigstar$ ）。
$\textbf{ \underline{\dag\dag\dag\dag《分子动力学后处理集成函数—Matlab》：}}$ 主要介绍针对后处理过程中指定函数，进行包装，方便使用者直接调用（需要一定编程能力，难度： $\bigstar$ $\bigstar$ $\bigstar$ $\bigstar$ ）。
$\textbf{ \underline{\dag\dag\dag\dag\dag《SCI论文绘图—Python绘图常用模板及技巧》：}}$ 主要介绍针对处理后的数据可视化，并提供对应的绘图模板（需要一定编程能力，难度： $\bigstar$ $\bigstar$ $\bigstar$ $\bigstar$ ）。
$\textbf{ \underline{\dag\dag\dag\dag\dag\dag《分子模拟—Ovito渲染案例教程》：}}$ 主要采用 $\rm Ovito$ 软件，对 $\rm Lammps$ 生成的轨迹文件进行渲染（难度： $\bigstar$ $\bigstar$ ）。

$\color{red}{\textbf{ \underline{专栏说明（订阅后可浏览对应专栏全部博文）：}}}$
$\color{red} 注意：$ 如需只订阅某个单独博文，请联系博主邮箱咨询。 $\rm lammps\_materials@163.com$

$\spadesuit$ $\dag$ 开源后处理集成程序：请关注专栏《LAMMPS后处理——MATLAB子函数合集整理》
$\spadesuit$ $\dag$ $\dag$ 需要付费定制后处理程序请邮件联系： $\rm lammps\_materials@163.com$


#--------------------------------------------------------------#
#      500-1000 1000-1500 1500-2000 2000-2500 2500-3000 
#--------------------------------------------------------------#
variable         KB           equal  8.625e-5

variable         2ns          equal  2000/${TS}

variable         per_loop_run equal  ${2ns}/10
#--------------------------------------------------------------#
###################
#     run 2       #
###################

variable          temp_incr loop 5 
label             loop

variable          temp_now equal (ke(all)/count(all)/${KB}/1.5)

print             "######################"
print             "Current temperature is: ${temp_now}"
print             "######################"

variable          current_target_temp equal (${temp_incr}+1)*500

fix               1    all nvt temp ${temp_now} ${current_target_temp}  ${Tdamp} 
#fix               1   all temp/rescale 1 ${temp_now} ${current_target_temp} 0 1

dump              loop all custom 10000 dump.${temp_incr}               &
                  id type x y z vx vy vz fx fy fz c_ke c_pe c_2[*]

thermo_style      custom step temp press ke  pe  vol  density 

timestep          ${TS}
run              ${per_loop_run}

#run              1

undump            loop 

unfix             1

next  temp_incr
jump  SELF loop
#--------------------------------------------------------------#