AMBER Archive (2004)

Subject: Re: AMBER: solute leaving PBC box

From: opitz_at_che.udel.edu
Date: Fri Oct 08 2004 - 08:43:06 CDT

('binary' encoding is not supported, stored as-is) So I have to run ptraj to image the system, even though iwrap was turned
on? And even though the solute is sticking outside of the iwarped box,
the molecule is still interacting with the solvent, thereby the
simulation is still correct?

Yes, you are right dt=0.001, I just typed the input on my windows machine
instead of copying it over from the unix workstation.

==============Original message text===============
On Thu, 07 Oct 2004 17:43:49 EST Carlos Simmerling wrote:

with iwrap, the molecule will be translated, but part of the molecule
may still move outside the box before the whole molecule gets
translated. this translation happens when the center of the molecule
moves out of the box (see subroutine wrap_molecules in ew_box.f).

in the future, you should send actual input files, since I doubt that
you really
ran with dt=0.01 (10fs), and there is no "tmp0", etc. it's very hard
to help people when what they say they did is not what they really did.

carlos

>Sorry, I should have made my question clearer.
>If iwarp is turned on, shouldn't the molecule be imaged back into the box
>(even, if that means "breaking" up the molecule and have it come into the
>box from different sides?). And based on that is the behavior of the
>molecule affected by this or is it still ok?
>
>Best Regards,
>
>Armin
>
>==============Original message text===============
>On Thu, 07 Oct 2004 16:48:12 EST Carlos Simmerling wrote:
>
>this is fine. you can use the image command in ptraj to
>make it easier to view (center on the molecule of interest
>and then image). check the ptraj documentation for details.
>
>
>
>>Dear Amber Community,
>>
>>I am running MD in a periodic box, with the following input:
>>&cntrl
>>imin=0, irest=1, ntx=5,
>>iwrap=1,
>>ntb=2, ntp=1, pres0=1, taup=1.0
>>ntf=1,
>>igb=0,
>>ntr=0, cut=12,
>>ntt=1, tempi=300, tmp0=300, tautp=1
>>nstlim=2600000, dt=0.01, t=400,
>>ntpr=1000, ntwx=1000, ntwr=1000
>>&end
>>At the end of this I took a look at my simulation, and my solute (the
>>molecule of interest), has migrated through the box to one of its
>>corners. Parts of the molecule are sticking out of the box, and aren't
>>imaged back into the box.
>>Does this affect the molecule behavior? If so how could i prevent it from
>>leaving the box?
>>Please help.
>>
>>Best Regards,
>>
>>Armin

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu