lammps模拟输出单个原子的能量

大家好，我是小马老师。
本文介绍lammps输出单个原子总能量的方法。
能量的计算与输出相对比较简单，本来没有写本文的计划，但最近有几个学员咨询到相关的问题，想想还是写出来吧，万一哪天能用的上呢。
在lammps模拟中，单个原子量的计算主要通过compute **/atom类的命令实现。
只要是带有atom关键词的compute命令，基本都是计算单个原子的量，总共41个计算命令，功能已经比较全面，但是没有单个原子总能量的计算命令。
虽然lammps没有提供总能量的计算方法，但是提供了势能和动能的计算命令。
只要计算出动能和势能，求和即可得到总能量。
原子的总能量=动能+势能
代码如下：

#计算单个原子的动能
compute Ke_atom all ke/atom
#计算单个原子的势能
compute Pe_atom all pe/atom
#求和得到总能量
variable E_total atom c_Ke_atom+c_Pe_atom

注意，variable命令中不能用equal，必须用atom关键词，否则会提示错误。
equal对应标量，atom对应矢量，因为每一个原子具有一个势能或者动能值，所以此处的能量均为矢量。
单个原子量的输出需要用dump命令：

dump 1 all custom 100 npt.xyz id type x y z v_E_total

原子能量值保存到轨迹文件中，用ovito就可以显示原子总能量云图：

公众号：lammps加油站

lammps模拟输出单个原子的能量相关推荐

lammps教程：平均值输出fix ave/time命令详解
lammps模拟输出方式有很多,前面几篇文章也介绍过,对于模拟总体系的值一般使用thermo命令输出,而单个原子的变量值使用dump命令输出. 今天主要讨论一下体系量值的输出,使用thermo命令有两 ...
Gromacs文件-Chapter1
Gromacs的文件非常的多,这是官方文档地址:http://manual.gromacs.org/online/files.html. 本文章部分内容来自以下网址https://zhuanlan.z ...
vasp 模拟退火_VASP 计算问题小结
本文转载自博主一个人就是一个叠加态,有部分删减修改,文中对相关概念方法做了详细的总结,留坑待填... 1. 第一原理计算的一些心得 1. 第一性原理第一性原理其实是包括基于密度泛函的从头算和基于Ha ...
在一个程序中什么算是亮点_VASP 计算问题小结
本文转载自博主一个人就是一个叠加态,有部分删减修改,文中对相关概念方法做了详细的总结,留坑待填... 1. 第一原理计算的一些心得 1. 第一性原理第一性原理其实是包括基于密度泛函的从头算和基于Ha ...
差分电荷密度 matlab,关注：电荷密度计算及图形化处理，电声耦合常数计算
1) 电荷密度计算及图形化处理: 2) 电声耦合常数计算: 3) Raman光谱计算与分析题记:好久没碰过计算了,一篇稿子搁置了一年,是时候重新拿起来了,希望两周内能完成:这段时间希望在朋友们的帮助 ...
关于LAMMPS能量最小化
本文来源:"分子动力学"公众号原文链接:https://mp.weixin.qq.com/s/f5xk0dqJtD4DET__j4m0Tw LAMMPS能量最小化有几种方法? 适 ...
lammps构建高熵合金模型＋结构优化初步筛选能量最小的结构
文章目录前言一.ubuntu并行安装Lammps 安装教程二.单个高熵合金建模参考资料操作步骤三.构建1000个高熵合金模型筛选能量最小的结构参考资料操作步骤 1.构建高熵合金模型 2 ...
lammps教程：如何计算单个原子的体积
大家好,我是小马老师. 本文介绍如何使用lammps命令计算单原子的体积. lammps提供了compute voronoi/atom命令计算单原子体积. 基本原理可参考voronoi算法,根据vor ...
lammps和atomsk中常用的建模方法
文章目录 1.Atomsk生成刃型位错 2.Atomsk生成任意角度晶界 3.Atomsk生成孪晶 4.Atomsk生成非晶 5.合金多晶结构建模 6.Atomsk建立管道类模型 7.Atomsk建立 ...

lammps模拟输出单个原子的能量

lammps模拟输出单个原子的能量相关推荐

最新文章

热门文章