OPLS全原子力场在lammps中的使用

opls_aa力场定义

OPLS是optimized potentials for liquid simulations,适用于液体体系。和其他势能函数一样,opls_aa势能函数包括键结作用和非键相互作用。键结作用包括harmonics键长和键角势能,Fourier扭矩势能,非键相互作用包括LJ势和Coulomb势[1]。力场表达式如下:


其中,LJ对相互作用采用混合法则采用标准混合法则,比如LB。lammps计算时注意扭矩表达式(3)中的系数1/2不包含在系数K中

lammps中opls_aa力场的应用

lammps中harmonics键长/键角势能,Fourier扭矩势能,LJ势和Coulomb势均有相应的命令与之对应。所以,opls_aa在lammps中采用如下指令处理:

atom_style full
bond_style harmonic
angle_style harmonic
dihedral_style opls
pair_style lj/cut/coul/long
pair_modify mix geometric
special_bonds lj/coul 0.0 0.0 0.5

采用以上命令组合时填入合适的系数即可。注意,全原子力场中1-2 1-3非键相互作用相对键长键角势能可忽略不计,1-4非键相互作用与二面角势能在一个数量级,这里考虑将之乘以0.5。不同力场的special_bonds设置不同。improper看所给的力场中是什么形式,harmonic或者cvff均有可能。

参考文献

[1] Watkins and Jorgensen, J Phys Chem A, 105, 4118-4125 (2001)

OPLS全原子力场在lammps中的使用相关推荐

  1. 全原子力场中各参数的含义及其在LAMMPS中的应用

    全原子力场中各参数的含义及其在LAMMPS中的应用 OPLSAA力场参数是分子动力学中广泛应用的力场.但有些同学对力场中的参数含义不是很了解.下面就OPLSAA力场参数如何应用于LAMMPS做一个简单 ...

  2. 了解lammps中NVE/NVT/NPT三种系综的特性和区别

    了解lammps中NVE/NVT/NPT三种系综的特性和区别: 系综是lammps模拟中的一个重要环节,系综种类及对应参数选择的直接决定了模拟结果程度的好坏.然而,由于各种体系.特定条件.系综种类和参 ...

  3. lammps中fix deform变形命令详细解释

    学习目标:理解和掌握lammps中fix deform变形命令 提示:不讲多余背景,精简版直接上手 主要内容: 一个番茄钟理解deform变形 一个番茄钟使用deform变形 理解deform: de ...

  4. 低秩矩阵补全算法matlab实现,推荐系统中的矩阵补全算法

    最基本的问题,以用户电影评分为例,也就是这个用户-电影矩阵. 表中是用户多电影的评分,但评分有缺失,因为用户不可能对所有电影作出评价. 那么推荐问题就是给用户合理推荐一个没看过的电影,合理是指,预测用 ...

  5. 电影下载全攻略 [初、中、高级]——老猫

    电影下载全攻略 [初.中.高级]--老猫 初级篇 出于保护版权的目的,RM和ASF这两种格式文件的播放器--RealPlay和MediaPlayer都没有提供SAVE AS选项.如果希望看完节目后在w ...

  6. 超级十全大补贴!!打造您心目中完美的TREO!!『上』

    超级十全大补贴!!打造您心目中完美的TREO!!『上』 我是去年8月用600,再到10月分入手650,一直到现在,不知不觉已经过了那么久了.快乐的时间总是过得快啊,怪不得王羲之说:快然自足,曾不知老之 ...

  7. 英文写一篇文章,名称为磺化聚醚醚酮/沸石分子筛(ZSM-5)杂化膜的制备及在全钒液流电池中的应用...

    磺化聚醚醚酮/沸石分子筛(ZSM-5)杂化膜的制备及在全钒液流电池中的应用 磺化聚醚醚酮/沸石分子筛(ZSM-5)是一种经过改性的分子筛材料,具有优异的电化学性能和良好的电解质分离能力.本文将介绍如何 ...

  8. 一、在全系1000名学生中,征集慈善募捐, 当总数达到10万元时就结束,统计此时捐款的人数,以及平均每人捐款的数目。

    #include <stdio.h> #include <stdbool.h> /*一.在全系1000名学生中,征集慈善募捐, 当总数达到10万元时就结束,统计此时捐款的人数, ...

  9. Lammps中常用的势函数和晶体库资源收集

    关注 Mr.material,\color{Violet} \rm Mr.material\ ,Mr.material , 更\color{red}{更}更多\color{blue}{多}多精\col ...

最新文章

  1. Android --- 很好用的时间选择器
  2. 让“云”无处不在-Citrix Xenserver之三 license server
  3. 防注入php 留言板代码,简单的 php 防注入代码
  4. js声明变量三种方式
  5. 卸载iis express后80端口仍然被占用的解决方法
  6. sql脚本比较大,sqlserver 无法导入,就用cmd命令执行
  7. 掐头去尾,只吃中间一段
  8. miniusbpcb封装_MINI MICRO A型B型D型USB多类型Altium AD元件库 PCB封装库的电路方案详解...
  9. 机器学习(Machine Learning)与深度学习(Deep Learning)资料汇总
  10. Tomcat9的新特性和Tomcat8的区别
  11. 如何构建高效可信的持续交付能力,华为云有绝活!
  12. 《PS技巧精华全部在这里!》
  13. Mysql数据库管理系统原理及基本操作
  14. 让你一分钟认识电子身份验证系统EID
  15. 从业了八年程序员觉得已经是非常的枯燥了
  16. 『HTML5梦幻之旅』-滚动播放的幻灯片效果
  17. 去掉Echarts饼状图的引导线
  18. 修改字符串指定范围的颜色
  19. 在区块链世界中的token到底是什么?
  20. 吉林省二级计算机vb c语言 access vf 考哪个有用,计算机二级考试选择

热门文章

  1. 大学物理·第2章牛顿定律
  2. 如何对待焦虑与拖延症
  3. 练习-Java类和对象之类的声明之复数运算
  4. [解决]尝试安装时出现“无法安装 Office(64 位或 32 位)”错误
  5. 怎样通过几何画板绘制相似三角形
  6. 【原理篇】一文读懂Faster RCNN
  7. mysql group by 后求和_MySQL分组求和GROUP BY
  8. touchGFX安装与初识
  9. 机房安全监控设备连接图及指令发送
  10. 补丁工具Quilt与Git