lammps教程:velocity命令三种使用方法
原子速度在分子动力学模拟中是一个重要的参数,不仅决定了当前时刻体系的温度,也决定了原子在下一个的位置。
速度的设置是分子动力学模拟必不可少的一个步骤。
在lammps中,使用velocity命令设置原子的速度,命令格式为:
velocity group-ID style args keyword value ...
其中,group-ID指定哪一组原子将要被设置速度,style决定了速度设置的类型。
常见的速度设置方式有以下三种:
(1)velocity create进行初始化
分子动力学模拟需要有一个初始的状态,因此,当原子模型建立完成,进行minimize能量最小化之后,需要进行速度初始化,体系的温度是由原子的温度计算得到,因此,这个过程也称为温度初始化。
使用的命令为velocity create,例如:
velocity all create 300 4928459
这一条语句的作用是对所有原子进行速度初始化,总体系的温度为300k,4928459为随机数,每个原子的速度被随机设置,默认的情况下,原子的速度符合正态分布。
也可以设置dist参数,使原子的速度符合高斯分布:
velocity all create 300 4928459 dist gaussian
当然,也可以对部分原子进行速度初始化,只需要改变原子组即可:
velocity mobile create 298 4928459
(2)velocity set设置原子在特定方向上的速度
在拉伸、剪切等动力学模拟中,需要这种一部分原子沿着某一个方向移动,可以使用velocity set命令设定特定方向的速度。
例如在沿着Z轴进行拉伸模拟时,需要把底部原子固定住,使上部(top)原子沿着Z轴施加一个速度,可以写为:
velocity top set NULL NULL 2 sum yes units box
上句代码中 NULL NULL 2 表示不设置x和y方向的速度,仅设置z方向速度,sum yes的意思是将速度2加到原子当前时刻速度分量上,如果不写sum yes,则表示忽略原来的速度,直接将z方向速度直接设为2。
默认的单位为晶格单位(lattice),units box关键字可将单位设置为实际速度单位,具体单位与体系units有关。
多个方向同时设置:
velocity top set 2 2 0 sum yes
将x、y方向移动速度设为2(晶格单位),z方向速度为0,覆盖掉原来的速度值。设置之后相当于top组原子沿着x和y方向的对角线方向移动。
(3)velocity ramp设置速度场
lammps允许在某一个方向上按照不同的坐标位置设置不同的速度,例如:
velocity flow ramp vx 0.0 5.0 y 5 15
这句代码是将y坐标从5到15的原子在x方向的速度设为0到5,y坐标小于5的原子x方向速度为0,y坐标大于15的原子x方向的速度为5,中间的原子x方向速度按照比例关系进行设定,例如y坐标为10的原子x方向上的速度为2.5。
掌握以上三种velocity用法,可以应付绝大部分的模拟场景,如果需要更复杂的速度设置,可以查lammps手册。
更多lammps教程请关注微信公众号:lammps加油站
lammps教程:velocity命令三种使用方法相关推荐
- xp系统查不到共享计算机名,xp系统查找不到网络打印机怎么办?xp系统找不到网络打印机的三种解决方法图文教程...
网络打印机相对共享打印机而言,实现了更加灵活.高效的打印任务,有时候在使用过程中也会遇到各种问题.这不一位xp纯净版系统用户说找不到网络打印机,导致无法正常添加网络打印机.有什么办法能够解决此问题呢? ...
- python调用cmd命令释放端口_详解python调用cmd命令三种方法
目前我使用到的python中执行cmd的方式有三种 使用os.system("cmd") 该方法在调用完shell脚本后,返回一个16位的二进制数,低位为杀死所调用脚本的信号号码, ...
- 计算机系统盘涨满怎么办,win7系统C盘突然就满了怎么办?win7系统盘突然爆满了三种解决方法图文教程...
最近,有win7系统用户反映,电脑的系统盘突然就爆满了,这该怎么办呢?大家可以释放C盘空间,下面脚本之家的小编和大家讲解win7系统盘突然就满了三种解决方法,一起来看看吧. win7系统盘突然就满了解 ...
- win7 计算机管理的命令,Win7运行cmd命令提示“命令提示符已被系统管理员停用”三种解决方法...
在运行win7操作系统中,经常使用运行命令对系统进行一些设置.最近有用户在Win7系统下运行CDM命令时,提示"命令提示符已被系统管理员停用"(如图),导致无法正常使用命令功能,遇 ...
- 将python中的小数直接进位的函数_python保留小数位的三种实现方法
python保留小数位的三种实现方法 前言 保留小数位是我们经常会碰到的问题,尤其是刷题过程中.那么在python中保留小数位的方法也非常多,但是笔者的原则就是什么简单用什么,因此这里介绍几种比较简单 ...
- 应用程序错误电脑黑屏_电脑运行程序出现APPCRASH错误的三种解决方法
好多小伙伴说电脑在运行程序时出现appcrash错误,导致软件无法正常运行,其实电脑是很复杂的程序代码设计的,遇到一些奇怪的问题也不奇怪,系统运行程序出现APPCRASH错误,很多新手用户都不知道怎么 ...
- Win11系统错误powershell找不到mscoree.dll的三种解决方法
近期有用户在更新微软Win11最新补丁KB5013943后,出现了严重 Bug,其中包含了powershell.exe - 系统错误,并提示找不到mscoree.dll,这该怎么解决呢?下面小编就为大 ...
- Using的三种使用方法
Using的三种使用方法 一.作为指令使用:用于未命名空间创建别名或导入其他命名空间中定义的类型(https://img-blog.csdnimg.cn/5ac458dbe91a4fba892df35 ...
- 桌面计算机未响应怎么办,电脑桌面假死(点击无反应)怎么解决?试试这三种处理方法吧...
原标题:电脑桌面假死(点击无反应)怎么解决?试试这三种处理方法吧 不少朋友提到遇到过电脑桌面假死,桌面图标怎么点击都没反应,最后无奈只能强制关机才能解决问题.有什么方法可以不用强制关机就能解决电脑假死 ...
- 基因注释 InterProScan的三种使用方法
InterProScan的三种使用方法 Interproscan,通过蛋白质结构域和功能位点数据库预测蛋白质功能.是EBI开发的一个集成了蛋白质家族.结构域和功能位点的非冗余数据库.Interpros ...
最新文章
- 结对编程——单元测试
- 补一点COM的知识,有关_com_ptr_t
- 如何发布一个npm包?
- php while mysql_我怎么能避免在PHP的While循环中使用MySQL查询
- Vim设置支持鼠标,光标定位,滚轮上下移动
- 智能机维修暴利大起底:触摸屏成本30维修300元
- AMD 芯片被曝大量安全漏洞,Linux 之父怒评!
- win10怎么设置默认输入法_电脑输入法怎么设置成搜狗双拼?
- 机器学习--线性回归4(线性拟合、局部线性拟合实战)
- 用好Windows 7自带文件加密工具
- -XX:NewRatio 命令
- sqli-labs靶场Less-20~22 Cookie注入
- C++/C输出素数,满某个数就换行
- WIN10 EXCEL 快捷键
- 猎豹移动推出直播产品Live.me风靡美国
- Windows用户的分类
- 警告关于测试人员的职场生存,千万要避开这5个坑(不看后悔)
- 关于全志R11如何驱动4线ST7789芯片的屏幕
- 尝试反汇编 FLUENT ADDON 之 SOFC
- c++中set的介绍及用法