Problem with including two potentials inside input script

Thank you Jocabi. But here I am using two different potentials with same pair style, reax/c. When I remove hybrid I get another error: ERROR: Incorrect args for pair coefficients (…/pair_reaxc.cpp:288)
Last command: pair_coeff * * ffield.reax.cho C C C C C C H H H H H H H H H H H H O O O O O O
The input file is:

REAX potential for ZnOH2 system

units real

atom_style charge
read_data out_vmd.data

pair_style reax/c NULL
pair_coeff * * ffield.reax.cho C C C C C C H H H H H H H &
H H H H H O O O O O O
pair_coeff * * ffield.reax.ZnOH O Zn

neighbor 2 bin
neigh_modify every 2 delay 10 check yes page 100000

dump 1 all xyz 10 file2.xyz
#dump_modify 1 element C C C C C C H H H H H H H H H H H H O O O O O O

fix 2 all qeq/reax 1 0.0 10.0 1e-6 reax/c
minimize 1.0e-10 1.0e-10 10000 10000

fix 1 all nve
fix 2 all qeq/reax 1 0.0 10.0 1e-6 reax/c
fix 3 all temp/berendsen 500.0 500.0 100.0

reset_timestep 0
timestep 0.25
run 3000

print “All done!”

And the masses in data file are arranged in this way:
Masses

1 12.011000 # C1
2 12.011000 # C2
3 12.011000 # C3
4 12.011000 # C4
5 12.011000 # C5
6 12.011000 # C6
7 1.008000 # H1
8 1.008000 # H1o
9 1.008000 # H2
10 1.008000 # H2o
11 1.008000 # H3
12 1.008000 # H3o
13 1.008000 # H4
14 1.008000 # H4o
15 1.008000 # H5
16 1.008000 # H6R
17 1.008000 # H6S
18 1.008000 # H6o
19 15.999000 # O1h
20 15.999000 # O2h
21 15.999000 # O3h
22 15.999000 # O4h
23 15.999000 # O6h
24 15.999000 # OR
25 15.999000 # OZ
26 12.000000 # ZN

Please help me.

Pair_style reax/c doesn’t work this way. You cannot combine two ReaxFF parameter files with pair_style hybrid – full stop.

Ray

Thank you Shan for the clarification. Could you please tell me how do I combine these two reaxxff?