任何之前的pair_coeff 和 pair_modify 命令设置将被保留。唯一例外的是,如果 pair_style 命令中全局截断距离被修改,它将覆盖之前任何 pair_modify 命令 中对应的截断距离。 不遵循此规则的两个对势样式是 pair_style table 和hybird 命令。新的 pair_style 命令将清除所有之前指定的 pair_coeff 和 pair_modify ...
lj/cut/coul/cut的截断半径有两个:普通LJ截断半径和库仑力截断半径。 如果只写一个截断半径,则默认cutoff2=cutoff。 #只设置一个截断半径pair_stylelj/cut/coul/cut10.0pair_coeff**100.03.0#设置两个截断半径pair_stylelj/cut/coul/cut10.08.0#原子1之间使用截断半径9.08.0pair_coeff11100.03.59.0#原子2之间使...
1. lj/cut在in文件中使用方法 在in文件中使用lj/cut,需要明确的指定相互作用的原子类型,如: pair_style lj/cut 10 pair_coeff * * 0.02 3.12 pair_coeff 1 3 0.03 3.22 8.5 上面两句指定全部原子之间的lj/cut力场,截断半径为10,epsilon=0.02、sigma=3.12。 1和3原子截断半径为8.5,epsilon=0.03、sigma=3.22。
lj/cut/coul/long力场公式与lj/cut/coul/cut类似,但是在计算库仑力时,不仅计算近邻原子之间的库伦力,也计算近邻原子之外的原子对中心原子的库伦作用力。 因此,需要配合kspace_style pppm或者kspace_style ewald命令使用。 lj/cut/coul/long使用方法与lj/cut/coul/long类似: kspace_style pppm 1.0e-5 pair_style ...
pair_style lj/cut 2.5 pair_style eam/alloy pair_style hybrid lj/charmm/coul/long 10.0 eam pair_style table linear 1000 pair_style none 说明 该命令用于设置每对原子间的相互作用的计算公式。在LAMMPS中,对势函数定义了截断半径内原子对的相互作用,随着原子间距离的变化,原子对间相互作用也不断变化。尽管...
atom_style atomic 2. lattice命令: 定义晶格结构(如果使用周期性边界条件)。 例如: lattice fcc 3.52 3. create_atoms命令: 创建原子,包括原子类型、位置等。 例如: create_atoms 1 box 4. pair_style命令: 选择力场模型。 例如: pair_style lj/cut 2.5 ...
1.单独使用eam势函数 pair style 和pair coeff 的写法: pair_style eam pair_coeff * * Cu.eam #cu-cu 2.混合势函数的写法: pair_style hybrid lj/cut/coul/long 10.0 lj/cut 10.0 eam pair_coeff 2 2 lj/cut/coul/long 0.006998 3.16435
3、entials(指定原子作用势)pair_stylelj/cut 2.5#指定lj势,截断半径为2.5pair_coeff * * 1.0 1.0 2.5#指定lj势参数# define groups(便于加载)Region 1 block INF INF INF 1.25 INF INFGroup lower region 1#定义lower组(便于施加外加速度)Region 2 block INF INF 38.75 INF INF INFGroup upper region 2...
我事先测试过 pair_style lj/cut 的 table potential 的形式，是可以正常运行的。希望有经验的朋友可以...
pair_style lj/cut/coul/cut command 是带电荷的L-J势，manual里面解释Qi and Qj are the charges ...