Clayff力场（转载）

本文转载自：http://blog.sina.com.cn/s/blog_b48a7ac30102w4km.html

注意：lammps中并没有完整的clayff势函数，需要将一下三者结合使用。

angle_style harmonic

bond_style harmonic

pair_style lj/cut/coul/long 12.5

使用过程中data文件很重要，键类型，键角，电荷都需要制定完备。键系数可以放在in文件中，也可以放在data文件中。键系数需要查阅相关文献或者和MS中GULP模块的clayff.frc势文件view看到。

V=Ar^(-12)-Br^(-6)

V=DRo^12/r^12-2DRo^6/r^6

DRo^12=A

2DRo^6=B

D=B^2/(4A) ~D单位eV

Ro^6=2A/B ~ 单位埃

1 KCal/mol=0.04336 eV

The functional form of the energy takes into account both bonded i.e. atoms linked by covalent bonds (Ebonded) and non-bonded interaction terms (Enon−bonded).

The total energy in general form can then be written as（一个系统的势能由化学键势能和非化学键势能两个部分组成）

Etot = Ebonded + Enon−bonded

Bond stretching describes the change in energy which occurs due to the

change in bond length from its equilibrium value
Angle bending describes the change in energy due to change in the angle
between two sequential covalent bonds from the equilibrium value
Torsion describes the change in energy of three bonds connected as IJ, JK
and KL due to change in the dihedral (or torsional) angle between the planes
IJK and JKL from the equilibrium value
Inversion describes the energy of three atoms bonded to one central atom
in the same plane due to an out of plane configuration

Ebonded = Ebond stretch + Eangle bend + Etorsion + Einversion

The non-bonded energy takes into account the electrostatic (coulombic) and van der Waals interactions and sometimes also the effect of hydrogen bonding，

Clayff力场（转载）相关推荐

VMDNAMD命令规则（转载）
输出体系的整个带电量: measure sumweights $all weight charge 给PDB文件设置周期边界条件: pbc set {54 54 24 } -all 将此晶胞内原子脱除 ...
[转载]从WiMAX看自主创新
大家都知道3G有WCDMA.CDMA2000.TD-SCDMA三个国际标准,却不知道还有第四个标准WiMAX:大家都以为TD很短寿,却不知道WiMAX支撑的时间更短:有的人觉得TD是个坑,却不知WiM ...
【转载】游戏引擎发展史漫谈(资料整理)
游戏引擎发展史漫谈导读游戏引擎(Game Engine)是什么?大多数人给出的解释都是基于engine的英文翻译,将游戏引擎等同于汽车中的引擎(发动机),再好一些的解释成动力,这些解释都对,但是动 ...
如何提高效率By Aaron Swartz （转载）
天才的逝去总是让世人扼腕不息,今天从微博看到Aaron Swartz commits suicide,生命的脆弱让人唏嘘.转载一篇Aaron Swartz关于如何提高效率的文章,以此向这位天才致敬. ...
生信在线工具（转载）
Web服务器名称网址简要描述;简介 agriGO v2 http://systemsbiology.cau.edu.cn/agriGOv2/ GO分析农业物种 AMMOS2 http://drug ...
在史蒂夫传中寻找现实扭曲力场1
Hello,小杨写一下读史蒂夫传的一些感受. 记得在2011年的某一天,我在电视上看到一位叫史蒂夫&乔布斯的人逝世了,当时只知道他是IT界的一位巨头,进入2011年苹果就风靡中国了的大江南北, ...
十大事件与Java相关的四十个名字(转载)
1990-1994:Java缘起文/孟岩 Larry Wall说,优秀程序员应有的三个特点:懒惰.急躁和傲慢.Java就是诞生在一群懒惰.急躁而傲慢的程序天才之中. 1990年12月,Sun的工程 ...
Jzoj4877 力场护盾(红警系列)
ZMiG成功粉碎了707的基因突变计划,为了人类的安全,他决定向707的科学实验室发起进攻!707并没有想到有人敢攻击她的实验室,一时间不知所措,决定牺牲电力来换取自己实验室的平安. 在实验室 ...
关于towhee 回答zhang_jaj的问题-转载小木虫
以下是zhang_jaj 的问题: 你好,我想用towhee做NVT下argon气在surface面上nucleation过程,希望你给一些指导,首先我想把potential选为LJ,surface与 ...

Clayff力场（转载）

Clayff力场（转载）相关推荐

最新文章

热门文章