lammps教程：均方位移MSD计算方法介绍

大家好，我是小马老师。
在lammps扩散模拟中，大多需要计算MSD，今天分享一段MSD的代码：

compute 1 all msd com yes
variable msdx equal c_1[1]
variable msdy equal c_1[2]
variable msdz equal c_1[3]
variable msd equal  c_1[4]
variable istep equal step
fix 2 all print 1 "${istep} ${msdx} ${msdy} ${msdz} ${msd}" screen no file msd.dat

第一行代码使用compute msd 命令对所有原子的均方位移进行计算。
计算完成后会产生4个结果，保存到数组c_1[]中。
其中，c_1[1]、c_1[2]、c_1[3]分别对应x、y、z三个方向的msd，c_1[4]为平均后的msd。
最后一行使用fix print命令把步数和4个msd保存到msd.dat文件中。
在需要计算msd的时候，直接把代码段插入到in文件中即可。

更多lammps教程，请关注微信公众号：
lammps加油站

lammps教程：均方位移MSD计算方法介绍相关推荐

5种方法计算均方位移MSD
5种方法计算均方位移MSD (1)lammps计算 (2)VMD计算 (3)OVITO计算 (4)ISAACS软件计算 (5)自编程计算-附MATLAB代码总结在统计力学中,均方位移(MSD,均方 ...
lammps教程：实例讲解npt、nvt系综的选择
对于lammps新手来说,系综的选择令人困惑,到底是应该用nvt系综还是npt系综? 我在初学lammps时就被这个问题困扰了很久.如果去问别人,可能会告诉你一个答案:根据模拟需要确定. 其实和没说一 ...
lammps教程：平均值输出fix ave/time命令详解
lammps模拟输出方式有很多,前面几篇文章也介绍过,对于模拟总体系的值一般使用thermo命令输出,而单个原子的变量值使用dump命令输出. 今天主要讨论一下体系量值的输出,使用thermo命令有两 ...
lammps教程：径向分布函数g(r)的计算与输出方法
大家好,我是小马老师. 本文介绍如何使用lammps计算径向分布函数g®. 径向分布函数(Radial distribution function)是指给定某个粒子的坐标,其他粒子在空间的分布几率. ...
lammps教程：Cu孪晶建模方法
上一次文章介绍了Atomsk软件建立孪晶的方法,有朋友咨询使用lammps自带命令建立孪晶的方法. 本文以Cu为例介绍lammps孪晶建模方法. 使用lammps建模,首先要确定xyz坐标轴对应的晶向 ...
lammps教程：lj/cut/coul力场参数设置详解
本文介绍lj/cut/coul/cut以及lj/cut/coul/long力场的区别和用法. 前面一篇文章介绍了LJ势的用法,普通的LJ势只考虑原子之间的吸引力与排斥力,没有考虑原子之间电荷作用. 如 ...
lammps教程：推荐一个不同单位下压强自动转换网站
大家好,我是小马老师. 在进行lammps模拟时,不同units下同一个物理量的单位有所不同. 例如,在metal单位下,压强的单位为bars. 在real单位下,压强的单位为atmospheres. ...
python爬虫scrapy框架教程_Python爬虫教程-30-Scrapy 爬虫框架介绍
从本篇开始学习 Scrapy 爬虫框架 Python爬虫教程-30-Scrapy 爬虫框架介绍框架:框架就是对于相同的相似的部分,代码做到不出错,而我们就可以将注意力放到我们自己的部分了常见爬虫框 ...
Arduino可穿戴开发入门教程Arduino开发环境介绍
Arduino可穿戴开发入门教程Arduino开发环境介绍 Arduino开发环境介绍 Arduino不像我们使用的PC端操作系统一样,可以直接在操作系统中安装软件为操作系统编程.Arduino的软件 ...

lammps教程：均方位移MSD计算方法介绍

lammps教程：均方位移MSD计算方法介绍相关推荐

最新文章

热门文章