Regarding Potential in a Cu-H-C system.

Hi,

In my system I have 6 types of atoms ( Cu,Cu,H and 3 layers of Graphene ©). I was trying to use rebo potential for the Graphene layers and lj/cut for the VDW between two graphene layers. And for the Cu-Cu-H interactions I was trying to use bond order potential.The potential portion of my code is shown below:

pair_style hybrid bop bop lj/cut 3.5 rebo rebo rebo rebo

pair_coeff * * CuH.bop.table Cu Cu H NULL NULL NULL

pair_coeff * * CCu_v2.bop.table Cu Cu NULL C C C

pair_coeff 4 5 lj/cut 0.00239 3.4 3.5

pair_coeff 5 6 lj/cut 0.00239 3.4 3.5

pair_coeff 4 6 lj/cut 0.00239 3.4 3.5

pair_coeff * * rebo 1 CH.airebo NULL NULL NULL C NULL NULL

pair_coeff * * rebo 2 CH.airebo NULL NULL NULL NULL C NULL

pair_coeff * * rebo 3 CH.airebo NULL NULL NULL NULL NULL C

pair_coeff * * rebo 4 CH.airebo NULL NULL H C C C.

But it is showing Pair coeff for hybrid has invalid style. Can I use bond order potential with hybrid pair style actually? Thanks in advance.

Your first 2 pair_coeff commands do not list which sub-style of hybrid
they are for (and the index if bop), hence the error.

Steve

Hi,

[...]

But it is showing Pair coeff for hybrid has invalid style. Can I use bond order potential with hybrid pair style actually? Thanks in advance.

yes, you can, but not in the way you seem to intend to do it. also
your choice of REBO cannot work. steve was telling you about the
syntax error in your input.

but there is also a grave conceptual problem: for manybody potentials,
atom types must be "wholly" contained. i.e. you can have a specific
atom type only mapped to one manybody potential at a time.
when you map two atom types, you also have all interactions between
those two atom types (pairs, triples, etc.) contained. you cannot
selectively override specific pairs or triples.
you can only selectively provide inter-type pairwise interactions with
pairwise additive potentials, e.g. how you use lj/cut (but the cutoff,
you are using there is odd. do you want to do a repulsive only
interaction? for a system that isn't?).

thus the first bob pair_coeff collides with the second and that
collides with all rebo coeff statements. the use of the first three
rebo potentials is correct, but the fourth is very wrong and also
collides with the second bop coeff statement.

in short, you have two syntax errors, but they are a minor issue
compared with a completely bogus choice of parameters and intentions
to mix and match interactions in an unphysical way.

axel.

Dear Axel,

Thanks a lot for your detailed reply. I got the point you are trying to indicate. But can you please suggest me what type of potential I should be using for Cu-graphene and H-graphene interactions? Thanks again.
Joy

Dear Axel,

Thanks a lot for your detailed reply. I got the point you are trying to indicate. But can you please suggest me what type of potential I should be using for Cu-graphene and H-graphene interactions? Thanks again.

this is not my area of research, so i don't know by heart. you have to
do what everybody has to do in this kind of situation and search the
literature.

axel.