[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
## [ESPResSo] Question about Tabulated Interaction

**From**: |
Lorenzo Isella |

**Subject**: |
[ESPResSo] Question about Tabulated Interaction |

**Date**: |
Tue, 1 Jul 2008 11:22:04 +0200 |

Dear All,
I would like to run some simulations to investigate the dynamics of
nanoparticles which are supposed to interact via a potential derived
from the integration of a Lennard-Jones potential over two spheres
(each particle is much larger than a molecule).
Now, I have the analytical expression for both the potential and the
force, so I am able to generate my own tabulated interaction. I simply
wonder how I should choose the minimum separation between these two
particles.
Namely, the potential blows up becoming infinitely repulsive when the
surface-to-surface distance, d, becomes zero and exhibits a minimum
above that distance before decaying quickly to zero.
So, physically, there is no possibility of an interpenetration of the
two particles, but I do not know how I should treat it numerically. It
seems that I have to specify the tabulated potential as a function of
the centre-of-mass-to-centre-of-mass separation between the two
particles, if I understood correctly.
But strictly speaking this is not defined for a
centre-of-mass-to-centre-of-mass distance smaller than 2r_p, where r_p
is the particle hard-core radius.
I think that a similar problem should arise also in the standard
Lennard-Jones potential treated for particle separations below r_off.
So, how is that dealt with in Espresso?
I am concerned that, during the evolution, a particle may end up in a
region where it cannot possibly have access to and that it would not
feel any repulsion at all.
What about introducing an artificial repulsive force for small separations?
I know this has been done in similar situations, but first I'd like to
know what is the standard treatment of such cases in Espresso.
Many thanks
Lorenzo

**[ESPResSo] Question about Tabulated Interaction**,
*Lorenzo Isella* **<=**