pair_modify mix geometric # TIP3P水分子参数 # O-O: ε=0.1521 kcal/mol, σ=3.1589Å pair_coeff 1 1 0.1521 3.1589 # O-H: 禁用LJ作用 pair_coeff 1 2 0.0 0.0 # H-H: 禁用LJ作用 pair_coeff 2 2 0.0 0.0 # OPLS-AA甲烷参数 # C-C: ε=0.0661 kcal/mol, σ=3.50Å pair_coeff 3 ...
# --- Init Section --- units real boundary p p p atom_style full pair_style lj/cut/coul/cut 10.0 pair_modify mix arithmetic special_bonds lj 0.0 0.0 0.5 coul 0.0 0.0 0.8333 bond_style harmonic angle_style harmonic dihedral_style fourier # --- Atom Definition Section --- read_data "...
pair_style body:体类型粒子间的相互作用。 pair_style bop: 基于量子力学计算 sigma 和 pi 键的键序势。 pair_style comb: 优化电荷作用的多体势。 pair_style coul/cut:截断库伦势。 pair_style dpd:耗散粒子动力学。 pair_style dsmc:直接模拟蒙特卡洛方法计算粒子碰撞。 ...
势函数为混合势,具体写法如下: pair_stylehybrideamlj/class2/coul/long${cutoff}${charge_cutoff}#9.5bond_styleclass2angle_styleclass2dihedral_styleclass2improper_styleclass2pair_modifymixsixthpowertailyesspecial_bondslj/coul001#PairCoeffspair_coeff1111eamAl_jnp_real.eampair_coeff11lj/class2/coul/long...
pair_coeff 1 1 lj/cut 0.10492547 1.7155 pair_modify mix arithmetic kspace_style none group CA type 1 fix 1 CA setforce 0.0 0.0 0.0 thermo 200 thermo_style custom step temp pe etotal press vol ke enthalpy evdwl epair ebond eangle edihed ...
neigh_modify every 1 delay 0 check yes #读取模型 read_data co2.data group CO2 type 1 2 #势函数类型设置 pair_style lj/cut/coul/cut 12.0 bond_style harmonic angle_style harmonic pair_modify mix geometric #参数设置 pair_coeff 1 1 0.0537 2.80 ...
pair_modify mix arithmetic #CLAYFF力场中的参数组合形式 pair_modify shift yes bond_style harmonic angle_style harmonic kspace_style ewald 1.0e-6 read_data P0.0001GPaeq-zlength.data #MS导出的data文件 set type 1 charge 1.50133929 #ao 由1.575改为1.49507186 ...
pair_modify mix arithmetic tail yes bong_style harmonic #键类型弹簧势 angle_style harmon...
手册里有提到有种情况下，lammps会没有错误信息输出就停止计算，还有就是reax不支持pair_modify mix。
pair_coeff 2 2 1.4564 0.0252 2.5 pair_modify mix arithmetic 如果我的第一行改为：pair_co...