Hi Xie Yinghong,<br>
Since one of your vectors is exactly perpendicular to the other two,
lets just focus on the base. And remember that what goes in 2D counts
in 3D too. Your base is a parallelogram with angles of 120 and 60
degrees. Now copy that several times. Then take one of the corner
points in the centre. Draw lines to mark the borders where any point in
the plane lies closer to that point than to any of the other corner
points. You will get a marked shape, which is generally called the
Voronoi region. In your case, you'll find that the Voronoi region of
your system is a perfect hexagon. With your third vector, you'll have a
perfectly regular hexagonal prism.<br>
Now if you do it for the other points two, you'll find that the same
hexagon will be located at each corner point in the grid created by the
copied parallelograms. Thus, you can equally well look at your system
as an infinite replicated system of hexagons. These will be connected
by the edges of the original unit cells. Imagine the first vector
aligned with the x-axis, and the central corner point at the origin
(thus with one hexagon centred at the origin) and you'll have one
neighbouring hexagon at +1x. Now you defined your second vector to have
an angle of 120 degrees. So in your case you'll take the second hexagon
at -0.5x, +0.866y. But if you look at the copied system, you'll find
that there's also a hexagon at +0.5x, +0.866y. That's equally valid to
take as a second copy, and if you draw a line connecting from the
origin to that, it will have an angle of 60 degrees with the positive
x-axis. In other words, with an angle of 60 degrees (and equal vector
lengths) you have defined the exact same infinite simulation system.<br>
I encourage you to draw the things I explained above. Also google for
'lattice', 'voronoi', 'triclinic' and do have a look at the paper of
Bekker I suggested before.<br>
I hope this helps,<br>
Tsjerk<br><br><div><span class="gmail_quote">On 10/14/05, <b class="gmail_sendername">Yinghong</b> &lt;<a href="mailto:xieyh@hkusua.hku.hk">xieyh@hkusua.hku.hk</a>&gt; wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">

<p>Dear Tsjerk:</p>
<p>Thanks for your kind help. And, I need furthur confirmation about the 
following things.</p>
<p>My box is actually based on hexagonal prism, but its cross-section is not a 
hexagon. In fact, it is only a monoclinic cell totoally with 6 faces in the 
whole box, instead of 8 faces in hexagonal prism. Accordingly, is my previous 
definition correct?</p>
<p>So, in&nbsp;such case,&nbsp;I used &quot;editconf -bt tric&quot; to build the box, is 
that ok?????&nbsp;As told by you, all the box types are&nbsp;the same for 
simulation, to this point, I&nbsp;need your&nbsp;confirmation, thanks.</p>
<p>By the way, what is the meaning of the third angle equal to 60 degree? I 
searched many previous achieves, 120 degree was always adopted.</p>
<p>Xie Yinghong</p>
<p>HongKong University</p><div><span class="e" id="q_106ece24c62cfd7b_1">
<p>&gt;Hi Yinghong,</p>
<p>&gt; I believe you want to set up a hexagonal prism according to your 
previous<br>&gt; mails, and the procedure you outline would be consistent with 
that. For<br>&gt;visualizing the system represented by a triclinic box the 
command for<br>&gt;trjconv is also correct. You might also want to try compact 
and rect for -ur<br>&gt;and compare the results, keeping in mind that it's all 
the same system, but<br>&gt;just different representations (interconversion of 
box types).</p>
<p>&gt;Do make sure that the box dimensions are large enough to contain 
your<br>&gt;protein when setting up the box with editconf. By the way, coming to 
think<br>&gt;of it, you may as well set the third angle to 60 degrees. That will 
give the<br>&gt;same box :p (that will describe the same lattice), but is more 
likely to<br>&gt;comply to Gromacs prerequisites for boxes.</p>
<p><br>On 10/13/05, Rodrigo Reston &lt;<a href="mailto:rodrigoreston@yahoo.com.br" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">rodrigoreston@yahoo.com.br</a>&gt; 
wrote:<br>&gt;<br>&gt; I'm not sure if your problem is this one... Have 
you<br>&gt; visualized what the triclinic shape looks like? If you<br>&gt; 
haven't, search for the word &quot;triclinic&quot; on Google<br>&gt; Images to see one and 
maybe then it'll help you design<br>&gt; the best approach to encase your system 
(by the way,<br>&gt; do that to the other shapes available in editconf).<br>&gt; 
Also, you might not need to specify the dimensions of<br>&gt; the triclinic box 
- perhaps, just using the option -d<br>&gt; and typing the distance (like [0.5] 
nm) from your<br>&gt; molecule to anyone of the triclinic box faces will 
do.<br>&gt; Rodrigo S. Reston, BSc.<br>&gt; UFMG, Brazil.<br>&gt;<br>&gt; --- 
Yinghong &lt;<a href="mailto:xieyh@hkusua.hku.hk" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">xieyh@hkusua.hku.hk</a>&gt; 
wrote:<br>&gt;<br>&gt; &gt; Dear users:<br>&gt; &gt;<br>&gt; &gt; I am very 
strange with the concept of triclinic, so<br>&gt; &gt; I need your instruction 
for the choice of my<br>&gt; &gt; simulation box.<br>&gt; &gt;<br>&gt; &gt; In 
my current system, the lengths of a b c are<br>&gt; &gt; respectively 5.0, 5.0, 
10.0nm, and the angles<br>&gt; &gt; between bc, ac, ab are 90, 90 and 120 
degrees.<br>&gt; &gt; Because of my misunderstanding of triclinic, 
maybe,<br>&gt; &gt; the following steps are somewhat wrong, can you help<br>&gt; 
&gt; me to point them out?<br>&gt; &gt;<br>&gt; &gt; 1. editconf -f mole.gro -o 
out.gro -bt tric -box 5<br>&gt; &gt; 5 10 -angles 90 90 120 -c<br>&gt; &gt; 
Should I use the option &quot;-bt tric&quot; here? If not<br>&gt; &gt; triclinic, 
possibly, &quot;-bt tric&quot; is wrong.<br>&gt; &gt;<br>&gt; &gt; 2. genbox -cp out -cs 
-p topol.top -o b4em.gro<br>&gt; &gt;<br>&gt; &gt; 3. grompp -v -f em -c b4em -p 
topol -o box.tpr<br>&gt; &gt;<br>&gt; &gt; 4.trjconv -f b4em.gro -s box.tpr -o 
b4em2.gro -ur<br>&gt; &gt; tric -pbc inbox ---- for visualization only.<br>&gt; 
&gt; What should my choice for options &quot;-ur &amp; -pbc&quot; in<br>&gt; &gt; this 
case?<br>&gt; &gt;<br>&gt; &gt; Urgent equiry and any help will be 
greatly<br>&gt; &gt; apprieciated.<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt; Xie 
Yinghong<br>&gt; &gt; Hong Kong University</p>

</span></div><br>_______________________________________________<br>gmx-users mailing list<br><a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">
http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please don't post (un)subscribe requests to the list. Use the<br>www interface or send it to <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users-request@gromacs.org">
gmx-users-request@gromacs.org</a>.<br><br></blockquote></div><br><br clear="all"><br>-- <br><br>Tsjerk A. Wassenaar, M.Sc.<br>Groningen Biomolecular Sciences and Biotechnology Institute (GBB)<br>Dept. of Biophysical Chemistry
<br>University of Groningen<br>Nijenborgh 4<br>9747AG Groningen, The Netherlands<br>+31 50 363 4336<br>