找回密码
 立即注册
搜索

ReaxFF力场格式规范

2024-6-12 21:17| 发布者: admin| 查看: 206| 评论: 0

摘要: General parameters(常规参数)本节的标题以 npar 开头,npar 是力场文件中存在的常规参数的数量。ehader 后面跟着 npar 行,每个行都包含一个参数值,后跟一个注释,例如:39 ! Number of general parameters50.0 ...
General parameters(常规参数)
本节的标题以 npar 开头,npar 是力场文件中存在的常规参数的数量。ehader 后面跟着 npar 行,每个行都包含一个参数值,后跟一个注释,例如:
 39     ! Number of general parameters
50.0000 !Overcoordination parameter

Atoms(原子)
原子参数部分从力场中存在的原子类型数开始,然后是三个额外的标题行和参数块,每个原子类型一个块。每个块由 4 行组成,从包含原子名称和 8 个参数值的行开始,格式为 (1x,a2,8f9.4),然后是三行,每行有 8 个参数值,格式为 (3x,8f9.4),例如:
3    ! Nr of atoms;cov.r;valency;a.m;Rvdw;Evdw;gammaEEM;cov.r2;#
       alfa;gammavdW;valency;Eunder;Eover;chiEEM;etaEEM;n.u.
       cov r3;Elp;Heat inc.;n.u.;n.u.;n.u.;n.u.
       ov/un;val1;n.u.;val3,vval4
 C    1.3817   4.0000  12.0000   1.8903   0.1838   0.   9000   1.1341   4.0000
      9.7559   2.1346   4.0000  34.9350  79.5548   5.   9666   7.0000   0.0000
      1.2114   0.0000 202.5551   8.9539  34.9289  13.   5366   0.8563   0.0000
     -2.8983   2.5000   1.0564   4.0000   2.9663   0.   0000   0.0000   0.0000
 H    0.7853   1.0000   1.0080   1.5904   0.0419   1.0206  -0.1000   1.0000
      9.3557   5.0518   1.0000   0.0000 121.1250   5.3200   7.4366   1.0000
     -0.1000   0.0000  62.4879   1.9771   3.3517   0.7571   1.0698   0.0000
    -15.7683   2.1488   1.0338   1.0000   2.8793   0.0000   0.0000   0.0000
... one more block ...

Bonds, angles, etc(键长,键角及其他)
在其余部分中,块键由两个或多个整数组成,每个整数都引用具有此索引的原子块。关键词中的整数取决于块类型(两个用于键,三个用于价角等)。 例如,下面的键参数块对应于上面所示的原子块的 C-H 键。
电荷参数部分以电荷类型的数量开头,后跟一个附加注释行。块的第一行采用 (2i3,8f9.4) 格式,第二行采用(6x,8f9.4):
   6      ! Nr of bonds;Edis1; LPpen;n.u.;pbe1;pbo5;13corr;pbo6
          pbe2;pbo3;pbo4;n.u.;pbo1;pbo2;ovc  orr
  1  1 156.5953 100.0397  80.0000  -0.8157  -0.4591   1.0000  37.7369   0.4235
         0.4527  -0.1000   9.2605   1.0000  -0.0750   6.8316   1.0000   0.0000
... five more blocks ...
对于其余部分,格式与键部分相似,只是它们没有额外的标题行,并且 Fortran 格式可能略有不同:(2i3,6f9.4) 表示非对角线,(3i3,7f9.4) 表示价角,(4i3,7f9.4) 表示扭转角,(3i3,4f9.4) 表示氢键。


鲜花

握手

雷人

路过

鸡蛋

最新评论

Archiver|手机版|小黑屋|凯算计算

GMT+8, 2025-4-2 22:20 , Processed in 0.102980 second(s), 19 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

返回顶部