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) 表示氢键。 |
Powered by Discuz! X3.5
© 2001-2024 Discuz! Team.