espressomd-users
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [ESPResSo-users] Espressomd-users Digest, Vol 76, Issue 9


From: Nairhita Samanta
Subject: Re: [ESPResSo-users] Espressomd-users Digest, Vol 76, Issue 9
Date: Wed, 1 Mar 2017 01:48:11 +0530

Hello Pierre!

Thank you very much for the information.

Regards,
Nairhita

On Tue, Feb 28, 2017 at 10:30 PM, <address@hidden> wrote:
Send Espressomd-users mailing list submissions to
        address@hidden

To subscribe or unsubscribe via the World Wide Web, visit
        https://lists.nongnu.org/mailman/listinfo/espressomd-users
or, via email, send a message with subject or body 'help' to
        espressomd-users-request@nongnu.org

You can reach the person managing the list at
        address@hiddenorg

When replying, please edit your Subject line so it is more specific
than "Re: Contents of Espressomd-users digest..."


Today's Topics:

   1. Re: Fwd:  Active matter (Pierre de Buyl)


----------------------------------------------------------------------

Message: 1
Date: Mon, 27 Feb 2017 21:33:31 +0100
From: Pierre de Buyl <address@hidden>
To: address@hidden
Subject: Re: [ESPResSo-users] Fwd:  Active matter
Message-ID: <address@hiddenx230>
Content-Type: text/plain; charset=us-ascii

Hi,

The rotational component of the motion is defined with quaternions as follows:

1. The position of the particles relative to the rigid body center of mass
position, in the laboratory reference frame, is

r = q r_B q*

Where r_B is the body frame position (fixed in time by definition) and q depends
on time.

2. There is a equation of motion for q:

dq / dt = 1/2 omega q

where omega is the angular velocity.

You can look the algorithm up in Martys and Mountain
doi:10.1103/PhysRevE.59.3733 that is the reference used by ESPResSo (see
src/core/rotation.cpp).

The force for the swim mechanism comes relies on the orientation of the swimmer
that is defined by the quaternion (hence the explanation above).

For more information, look for "swimming" in the user guide.

Regards,

Pierre

On Fri, Feb 24, 2017 at 12:31:47PM +0530, Nairhita Samanta wrote:
> Dear Joost,
>
> Thank you for your reply. I have seen this paper earlier. However, I still
> exactly dont understand the mechanism used in Espresso for the swimming
> feature, that how this orientation (*u*) is updated. Does it follow a
> separate equation of motion like a Langevin equation for rotation, which
> gets updated at each step along with position and then the force (F) is
> applied in that direction? Otherwise it is difficult for Espresso to find
> the orientation of a symmetric single particle. I am trying to understand
> for a system without any Hydrodynamic interaction.
>
> Regards,
> Nairhita
>
> On Thu, Feb 23, 2017 at 9:47 PM, Joost de Graaf <address@hidden>
> wrote:
>
> > Dear Samantha,
> >
> > Perhaps reading the 2016 J. Chem. Phys. by De Graaf et al. (144, 134106)
> > on anisotropic swimmers will be helpful.
> >
> > Best Wishes,
> >
> > Joost
> >
> > On Thu, 23 Feb 2017 at 12:31, Nairhita Samanta <address@hidden>
> > wrote:
> >
> >> Hi all!
> >>
> >> I am using the v_swim/f_swim feature of ESPResSo. But I really want to
> >> understand the more clearly. So, can you please guide me to find the
> >> algorithm, or at least the articles which is followed for the equation of
> >> motion of the quaternions.
> >>
> >> Thanking you.
> >>
> >> With regards,
> >> Nairhita
> >>
> >>
> >>

--
-----------------------------------------------------------
Pierre de Buyl
KU Leuven - Institute for Theoretical Physics
T +32 16 3 27355
W http://pdebuyl.be/
-----------------------------------------------------------



------------------------------

Subject: Digest Footer

_______________________________________________
Espressomd-users mailing list
address@hidden
https://lists.nongnu.org/mailman/listinfo/espressomd-users


------------------------------

End of Espressomd-users Digest, Vol 76, Issue 9
***********************************************


reply via email to

[Prev in Thread] Current Thread [Next in Thread]