Re: [ESPResSo-users] track a particle in each time step

[Narges Nikoofard]

Dear Axel and Stefan,

Many thanks for the answers.

In more details, in my problem, a single polymer is confined eg. inside a sphere. I let the polymer to fluctuate in the sphere for a long time. A small volume is assumed inside this confining sphere. I need the simulation time, when one end of the polymer enters this volume. The simulation continues until the end monomers enters the specified volume, 100 times.

I read Sec. 9.1 of the userguide. I do not know which of the opinions below is better.

1. Write the positions of the two end monomers in an output file, then analyze this file with a separate C program. I see in the userguide that observable "particle_positions" is defined and specification of the monomers is possible. However, my simulations are commonly very long, and the output file will be very huge.

2. Write the required simulation times (that end monomers enter the specified volume) in an output file. I do not know if it is possible with the observable command.

Also, I would be grateful if I can have an example on using the command "observable".

Many thanks again,

Narges

On Tue, Mar 18, 2014 at 2:21 AM, Stefan Kesselheim <address@hidden> wrote:

Hi Narges,
as Axel pointed out, the observable concept might solve your problem. It however depends on what exactly you want to do with the information if certain monomers are in certain places. You of course need to do the postprocessing on the c-level. If you tell me more details, I might be able to suggest a solution.
Cheers
Stefan

On Mar 16, 2014, at 9:59 AM, Narges Nikoofard <address@hidden> wrote:

> Dear Espresso users,
>
> I am simulating a confined polymer. I should count the number of times the end monomers enter a pre-defined region in the space. So, I need to track the end monomers' position. I need to know the end monomers' position in each time step of the simulation.
>
> As a result, I have used "integrate 1" in my code. So, I can get the monomers position in each time step and check if they are in the specified region. This makes the code very slow. Is there any way to solve this problem?
>
> Many thanks in advance,
> Narges
>
>