分子动力学理论篇（1）——牛顿力学、哈密顿动力学和相空间

Mr. Material

已于 2023-03-05 11:10:47 修改

阅读量2.3k

点赞数 5

CC 4.0 BY-SA版权

分类专栏：分子动力学理论文章标签： lammps

于 2020-10-06 20:06:31 首次发布

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

分子动力学理论专栏收录该内容

2 篇文章

订阅专栏

该系列教程详细介绍了LAMMPS分子动力学模拟的各个方面，包括安装、基本原理、小技巧、In文件详解、编程技巧、后处理程序和Matlab集成函数，以及使用Python进行SCI论文绘图和Ovito渲染案例。内容涵盖了从入门到进阶的多个层次，适合不同需求的读者学习。此外，还提供了定制后处理程序的咨询服务。

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

关注 $\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$

分子动力学理论-总结

在这里插入图片描述

接下来我们画出harmonic spring 在相空间的轨迹：在这里插入图片描述
MALAB 代码如下：

clc;clear;
%--------------------------------%
% initial setting 
step = 0.1;
X_initial  = -4;
V_initial = 0;
lambda    = 1;
mass      = 1;
X_current = X_initial;
V_current = V_initial;
t         = 0;run = 0;
%--------------------------------%
% Run
while (t<=100)
  t   = t+step;
  half_t    = 0.5*step;
  a = -X_current;
  V_current = V_current+a*step;
  X_current = X_current+V_current*step;
  momentum  = V_current;
 %---------------------------------------%
 % plot 
  hold on;  
  s = scatter(X_current,momentum);
  s.LineWidth = 0.6;
  s.MarkerEdgeColor = 'b';
  s.MarkerFaceColor = [0 0.5 0.5];
  run = run +1;  
 %---------------------------------------%
 % output 
  if(mod(run,100)==0)
      fprintf("The step %d is finished!\n",run);      
  end 
end
disp("--------------------------");
disp("--------ALL DONE!!!-------");
disp("--------------------------");