[CP2K:7606] Re: how to keep the atom's trajectory within the box

王胜杰 theancientg... at gmail.com
Thu Mar 24 15:37:58 UTC 2016

Dear Huan,

I didn't display the trajectory, because the trajectory file is really big, 
thus it can't be dragged to VMD.  I wrote a little program to identify it. 
First I found the min and max coordinates in the initial configuration, in 
ensemble NVT, when the atom coordinates less than min or larger than max, 
identify it crossing the boundary. And print out which atom and during 
which step.
Sure I can split the trajectory into several parts, but splitting costs too 
much time, and VMD always crash.

Best wishing,
Shengjie Wang 
在 2016年3月24日星期四 UTC+8下午6:07:21,Huan Wang写道:
> Hi Shengjie,
> An additional question is how do you display the trajectory? 
> If you use VMD, then you need to type some codes in *Extensions --> Tk 
> console* to let the VMD know the size of your simulation box.
> The codes are:
> pbc set {xxx xxx xxx} -all 
> # Here the "xxx" is the length of your (cubic) box. And "-all" means this 
> operation controls all frames in the current trajectroy.
> pbc box
> pbc wrap -all 
> # The wrap commend will show the atoms on the other side of the box when 
> atoms move out of the box.
> Good luck,
> Huan
> .
> On Thu, Mar 24, 2016 at 3:39 AM, 王胜杰 <theanci... at gmail.com 
> <javascript:>> wrote:
>> Hi Huan,
>> Thanks for your attention.
>> Yes, I did. I noticed and added the keyword the PERIODIC, and I tried the 
>> ensemble NPT_I and NVT, but cp2k still printed a real trajectory.
>> Regards,
>> Shengjie Wang
>> 在 2016年3月23日星期三 UTC+8上午4:45:33,hua... at mail.huji.ac.il写道:
>>> Hi Shengjie,
>>> Did you specify the keyword PERIODIC in the &SUBSYS &CELL section?
>>> Please see the CP2K manual
>>> https://manual.cp2k.org/cp2k-2_6-branch/CP2K_INPUT/FORCE_EVAL/SUBSYS/CELL.html#PERIODIC
>>> Regards,
>>> Huan
>>> .
>>> On Friday, March 18, 2016 at 6:41:41 PM UTC+2, 王胜杰 wrote:
>>>> Dear all
>>>> When I was dealing with the trajectory of MD, I found the atoms always 
>>>> move beyond the boundary, making some atoms' coordinate larger than the box 
>>>> length,which reflect the true trajectories. 
>>>> But I don't want a true trajectory, I want a periodic boundary.When the 
>>>> atoms cross the boundary, the coordinate will minus the box length 
>>>> automatically, which each direction of  atoms' coordinate is less than the 
>>>> box length. By adding what command line can I fix that?
>>>> Thanks in advance.
>>>> Shengjie Wang 
>>> -- 
>> You received this message because you are subscribed to a topic in the 
>> Google Groups "cp2k" group.
>> To unsubscribe from this topic, visit 
>> https://groups.google.com/d/topic/cp2k/A083JgzAgLo/unsubscribe.
>> To unsubscribe from this group and all its topics, send an email to 
>> cp2k+... at googlegroups.com <javascript:>.
>> To post to this group, send email to cp... at googlegroups.com <javascript:>
>> .
>> Visit this group at https://groups.google.com/group/cp2k.
>> For more options, visit https://groups.google.com/d/optout.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20160324/5f26b170/attachment.htm>

More information about the CP2K-user mailing list