What can I use instead of compute group/group?

Hello all

I want to print out the force (fx fy fz or overall force) of several atom groups during MD simulation. I know there is a good command for this “compute group/group”, but I can’t use this command since I use multi-body potential (reax/c) My atom groups consist of different types of atoms, and each groups have different size and number of atoms.

Is there any compute or fix command to calculate the fx, fy, and fz of group of atoms? I already have per-atom fx fy and fz by dump command, but this is not what I want to get.

Thanks

Best,

Joon

Take a simple example for a simple many-body potential, like EAM.

The force on one atom due to some small number of other atoms

is not easily defined, since those other atoms combine with

all the atoms to create terms in the potential, like the

embedding force on the one atom.

There is no code in LAMMPS that takes derivatives

of arbitrary many-body potential formulas to get the force

on atom I due to just atom J. So I don’t think you are asking for

something that is easy to calculate (or means anything useful

in most contexts).

Steve

In addition to that, you can look into this recent work: http://arxiv.org/pdf/1503.06565.pdf . Authors derive a universal expression of the pairwise force between two atoms, for a general classical many-body potential (see equation 57 in particular), validated numerically. However, that’s not implemented in LAMMPS and you may want to double check if it actually applies to ReaxFF.

Arthur

So, there is no way to calculate the interaction energy when using Reax FF?
compute group/group is actually a very useful command :frowning:

So, there is no way to calculate the interaction energy when using Reax FF?

not as it is currently implemented in LAMMPS.

compute group/group is actually a very useful command :frowning:

compute group/group is dependent on your pair style being a pair-wise
additive potential.

of course you are free to try and implement the feature from the
referenced pre-print...

axel.

Nice :slight_smile:

Is it possible to run the main simulation by Reax FF, and then rerun it by defining another FF (a pairwise additive potential) to calculate the interaction energy?? At least this can give an approximation of the energy.

Nice :slight_smile:

Is it possible to run the main simulation by Reax FF, and then rerun it by
defining another FF (a pairwise additive potential) to calculate the
interaction energy?? At least this can give an approximation of the energy.

yes, if there is a separation between the two groups, the
approximation can be quite good.
if there are bonds between the two groups, it may be tricky.

axel.

A couple of comments:

1. The reason why there is no group/group feature for non-pairwise
additive potentials is that the many-body interaction energy between
groups is not a well-defined mathematical quantity.
2. Since you started out asking about forces and later switched to
energy, you are clearly not looking for a specific mathematical
quantity. Rather, you seem to after the subjective concept of whether
two groups are interacting or not.
3. Interaction between two groups is strongly correlated with spatial
proximity of member atoms.
4. So what you really need is a way to measure whether or not two
groups have member atoms that are near each other, a kind of minimum
distance measure.
5. This is easily calculated in a post-processing step. There might be
a way to do it in LAMMPS on the fly, but I'm not sure it is worth the
trouble.