Re: [ESPResSo-users] Hard Dimers in 2D

[Salvador H-V]

Hi Rudolf,

Thanks a lot for your comments.

I am interested in the patch to set the Langevin Gamma for rotation, that you mentioned in your email. 

I would like to decrease the rotation of the dimers. In my current simulations the dumbbells rotate very fast; this is an unphysical scenario when particles are supposed to be immersed in a solvent.

Could you please send the patch to me. Also, some guidance of how to implement the patch would be great.

Thanks a lot for your help,

Salvador

On Thu, Jan 9, 2014 at 2:13 AM, Rudolf Weeber <address@hidden> wrote:

Hi,

On Wed, Jan 08, 2014 at 03:38:59PM -0600, Salvador H-V wrote:
> I want to simulate a two dimensional system of  hard tangent spheres
> (dimers).
>
> I am interested in a rigid dumbbell, so, I need to avoid sliding of one
> particles on the surface of the other particle.
>
> I was wondering what is the proper way to setup this system with ESPresSo.
>
> I am considering two options:
> i) rigid bond with length bond equal to the diameter of the particles  +
> exclusion command
>
> ii) Virtual sites to create a rigid object + ¿exclusion command?

If you only need to keep the distance between the two spheres forming the dumbbell constant, both will do the trick.
If you decide to use virtual sites, place the non-virtual particle in the center of the dumbbell, don't define any interaction for it, and vs_auto_relate the two spheres to it. Define LJ potentials for the two spheres.

The question of "sliding" is only relevant, if you couple something to the rotational degree of freedom of the individual spheres. Tis would be the case, e.g., if you assign the individual spheres forming the dumbbell a magnetic moment.
If that is, what you need, let me know.
Otherwise, you can just ignore the rotational degree of freedom of the individual spheres.

> Finally, I would like to know the reliability of the Langevin Dynamics
> thermostat, to calculate the dynamical properties of this system, such as
> the Mean Square Displacement of the center of mass of the dimers. What
> about rotational properties?

The rotational degrees of freedom follow a similar Langevin equation as the translational ones.
You cannot set the Langevin gamma for the rotation independently, right now. Should you need that, I have a (littlte tested) patch, which I can provide if need be.

The Langevin thermostat influences the dynamics of the system.
I think, at least in cases where the Langevin thermostat is used just for thermalization and not to model an implicit solvent, people equilibrate with the thermostat and then turn it off for the MSD measurement, but I have no personal experience with that.


Hope, that helps,
Regards, Rudolf

--
=o=o=o=o=o=o=o=o=o=o=o=o=o=o=o=o=o=o=o

Dr. Salvador Herrera-Velarde
Simulacion Computacional de Materia Condensada Blanda
Sub-Direccion de Investigacion y Posgrado
Instituto Tecnologico Superior de Xalapa (ITSX)
Instituto Tecnológico Superior de Xalapa
Sección 5A Reserva Territorial s/n, 

Col. Santa Bárbara, 

CP 91096
Xalapa, Ver. Mexico


=o=o=o=o=o=o=o=o=o=o=o=o=o=o=o=o=o=o=o

Este correo ha sido editado para evitar el uso de acentos y guardar compatibilidad entre
diferentes distribuciones