field主命令

作用

设置场参数，如相互作用力场的计算类型（势函数类型）及相关参数、体系温度场等。

形式

field -estyle potential/helmholtz
      -vstyle virial
      set temperature value
      use name ...
      set name ...
      coeff name ...

说明

• field -estyle [potential|Helmholtz] : 设置HYMACS进行势能或自由能极小化，默认为自由能极小化。

• field -vstyle virial : 设置体系压力分量的计算形式，默认采用维里应力来计算。

• field temperature value : 设置体系的温度值为value。

• field use name ... : 设置HYMACS是否使用name类型的势函数，并进行相应的参数设置。目前支持的作用势类型有：lj/cut，morse，eam，buckingham。可以进行一步通过field set name和 field coeff name来设置势函数的相关选项以及势函数参数。

目前支持的势函数类型和设置形式有：

lj/cut势函数

• field use lj/cut
• field set lj/cut –offset [yes|no]
• field coeff lj/cut type1 type2 rcut epsilon sigma

其中type1、type2表示相互作用的两种原子类型，rcut表示截断半径，epsilon， sigma分别表示对应的势函数参数 和 。–offset选项表示是否计入截断半径之后的势能，保持势函数的连续性，一般选择yes。

morse势函数

• field use morse
• field set morse –offset [yes|no]
• field coeff morse type1 type2 rcut D0 alpha r0

其中type1、type2表示相互作用的两种原子类型，rcut表示截断半径，D0，alpha，r0分别表示对应的势函数参数。–offset选项表示是否计入截断半径之后的势能，保持势函数的连续性，一般选择yes。

eam势函数

field use eam file itype typeName … field coeff morse type1 type2 其中type1、type2表示相互作用的两种原子类型，file表示势函数数据文件（一般由开发者成生），itype表示HYMACS中的类型标号，typeName表示数据文件中的类型名称。

示例

field -estyle helmholtz
field use eam Cu.eam 1 Cu
field use morse
field coeff eam 1 1
field coeff morse 1 2 7.0 0.087 5.140 2.05 
field coeff morse 2 2 7.0 2.4230 2.5550 2.5220