Dear all，
Thank you for your reply that help me a lot！But the convergence problem didn’t solve，The qeq parameters is according to the example for ZnOH2 file. I think that the parameters for qeq is correct. I don’t know what’s wrong?
###########The input file#####################################
units metal
atom_style charge
boundary p p p
variable a index 3.35
variable b equal $asqrt(3)0.5
variable c equal -$a0.5
variable d equal $a1.63
lattice custom 1 &
a1 $b $c 0.0 &
a2 0.0 $a 0.0 &
a3 0.0 0.0 $d &
basis 0.3333333333333333 0.6666666666666666 0.3825999999999995 &
basis 0.6666666666666666 0.3333333333333333 0.8825999999999994 &
basis 0.3333333333333333 0.6666666666666666 0.0000000000000000 &
basis 0.6666666666666666 0.3333333333333333 0.5000000000000000 &
region box block 0 10 0 5 0.25 5.25
create_box 2 box
create_atoms 1 box basis 1 1 basis 2 1 basis 3 2 basis 4 2
mass 1 15.999400
mass 2 65.408997
pair_style hybrid/overlay buck 10.0 morse 10.0
pair_coeff 1 1 buck 234.7246659 0.626654164 0.0
pair_coeff 2 2 buck 146.2713521 0.574294176 0.0
pair_coeff 1 2 morse 5.0955 1.89 1.98438
fix 1 all npt temp 300 300 0.1 z 0.0 0.0 0.1
fix 2 all qeq/reax 1 0.0 10.0 1e-6 param.qeq
timestep 0.001
dump 1 all atom 30 dump.reax.znoh
run 3000
###########The param.qeq file#####################################
1 8.5000 8.3122 1.0898
2 2.0219 5.7915 0.4828
