原子速度在分子动力学模拟中是一个重要的参数,不仅决定了当前时刻体系的温度,也决定了原子在下一个的位置。

速度的设置是分子动力学模拟必不可少的一个步骤。

在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命令三种使用方法相关推荐

  1. xp系统查不到共享计算机名,xp系统查找不到网络打印机怎么办?xp系统找不到网络打印机的三种解决方法图文教程...

    网络打印机相对共享打印机而言,实现了更加灵活.高效的打印任务,有时候在使用过程中也会遇到各种问题.这不一位xp纯净版系统用户说找不到网络打印机,导致无法正常添加网络打印机.有什么办法能够解决此问题呢? ...

  2. python调用cmd命令释放端口_详解python调用cmd命令三种方法

    目前我使用到的python中执行cmd的方式有三种 使用os.system("cmd") 该方法在调用完shell脚本后,返回一个16位的二进制数,低位为杀死所调用脚本的信号号码, ...

  3. 计算机系统盘涨满怎么办,win7系统C盘突然就满了怎么办?win7系统盘突然爆满了三种解决方法图文教程...

    最近,有win7系统用户反映,电脑的系统盘突然就爆满了,这该怎么办呢?大家可以释放C盘空间,下面脚本之家的小编和大家讲解win7系统盘突然就满了三种解决方法,一起来看看吧. win7系统盘突然就满了解 ...

  4. win7 计算机管理的命令,Win7运行cmd命令提示“命令提示符已被系统管理员停用”三种解决方法...

    在运行win7操作系统中,经常使用运行命令对系统进行一些设置.最近有用户在Win7系统下运行CDM命令时,提示"命令提示符已被系统管理员停用"(如图),导致无法正常使用命令功能,遇 ...

  5. 将python中的小数直接进位的函数_python保留小数位的三种实现方法

    python保留小数位的三种实现方法 前言 保留小数位是我们经常会碰到的问题,尤其是刷题过程中.那么在python中保留小数位的方法也非常多,但是笔者的原则就是什么简单用什么,因此这里介绍几种比较简单 ...

  6. 应用程序错误电脑黑屏_电脑运行程序出现APPCRASH错误的三种解决方法

    好多小伙伴说电脑在运行程序时出现appcrash错误,导致软件无法正常运行,其实电脑是很复杂的程序代码设计的,遇到一些奇怪的问题也不奇怪,系统运行程序出现APPCRASH错误,很多新手用户都不知道怎么 ...

  7. Win11系统错误powershell找不到mscoree.dll的三种解决方法

    近期有用户在更新微软Win11最新补丁KB5013943后,出现了严重 Bug,其中包含了powershell.exe - 系统错误,并提示找不到mscoree.dll,这该怎么解决呢?下面小编就为大 ...

  8. Using的三种使用方法

    Using的三种使用方法 一.作为指令使用:用于未命名空间创建别名或导入其他命名空间中定义的类型(https://img-blog.csdnimg.cn/5ac458dbe91a4fba892df35 ...

  9. 桌面计算机未响应怎么办,电脑桌面假死(点击无反应)怎么解决?试试这三种处理方法吧...

    原标题:电脑桌面假死(点击无反应)怎么解决?试试这三种处理方法吧 不少朋友提到遇到过电脑桌面假死,桌面图标怎么点击都没反应,最后无奈只能强制关机才能解决问题.有什么方法可以不用强制关机就能解决电脑假死 ...

  10. 基因注释 InterProScan的三种使用方法

    InterProScan的三种使用方法 Interproscan,通过蛋白质结构域和功能位点数据库预测蛋白质功能.是EBI开发的一个集成了蛋白质家族.结构域和功能位点的非冗余数据库.Interpros ...

最新文章

  1. 结对编程——单元测试
  2. 补一点COM的知识,有关_com_ptr_t
  3. 如何发布一个npm包?
  4. php while mysql_我怎么能避免在PHP的While循环中使用MySQL查询
  5. Vim设置支持鼠标,光标定位,滚轮上下移动
  6. 智能机维修暴利大起底:触摸屏成本30维修300元
  7. AMD 芯片被曝大量安全漏洞,Linux 之父怒评!
  8. win10怎么设置默认输入法_电脑输入法怎么设置成搜狗双拼?
  9. 机器学习--线性回归4(线性拟合、局部线性拟合实战)
  10. 用好Windows 7自带文件加密工具
  11. -XX:NewRatio 命令
  12. sqli-labs靶场Less-20~22 Cookie注入
  13. C++/C输出素数,满某个数就换行
  14. WIN10 EXCEL 快捷键
  15. 猎豹移动推出直播产品Live.me风靡美国
  16. Windows用户的分类
  17. 警告关于测试人员的职场生存,千万要避开这5个坑(不看后悔)
  18. 关于全志R11如何驱动4线ST7789芯片的屏幕
  19. 尝试反汇编 FLUENT ADDON 之 SOFC
  20. c++中set的介绍及用法

热门文章

  1. 我用Anylogic写了一个遗传算法
  2. Java教程(一)---JDK和Maven安装配置
  3. 什么录播软件好用?超级好用的录屏软件在这里
  4. mysql修改视图字段长度_SQL Server 数据库创建视图时修改字段长度
  5. app测试用例注意点
  6. 适合新手练习的二十套Java项目源码大分享
  7. c语言万能源代码,经典C语言源代码
  8. 怎么检查计算机和打印机是否连接网络,检查电脑是否正确连接网络打印机
  9. TransE算法的整理
  10. H5+springboot(集成ffmpeg)实现前端视频录制以及webm格式转mp4