<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
</style>
</head>
<body class='hmmessage'>
Hi,<br><br>I tried to reproduce this, but I did not succeed.<br>mdrun always seems to write the correct z size in all output files.<br><br>Are you sure mdrun outputs a file (and which file, if so),<br>or is it another Gromacs program that produced the file with 0 z size?<br><br>You can simply manually correct the size in the file.<br>But please report back which program produced which incorrect file.<br><br>Berk<br><br><hr id="stopSpelling">Date: Wed, 6 May 2009 00:39:18 +0800<br>From: crickzhang1@gmail.com<br>To: gmx-users@gromacs.org<br>Subject: [gmx-users] help: z-direction component of box vector becomes 0.00000 when using pdc = xy after<br><br><div>Hi, all</div>
<div> </div>
<div>There are some things confusing and resulting in error when I run mdrun for energy minization and md.</div>
<div> </div>
<div>I use Gromacs 4.0.4 to simulate electroosmotic flow in slab geometry with only 2 dimensional </div>
<div>periodic boundary (x and y), the system is aqueous solution with some ions, eg. CL- and Na+,</div>
<div>flowing through a plate channel, with external electrical field applied along x-direction, the height of channel</div>
<div>is z-direction. </div>
<div> </div>
<div>Coulomb interactions are computed with PME, and use pbc=xy, after running EM, all ouput tpr box sizes </div>
<div>in z-direction becomes zero, eg. originally the box size is 4.50000,4.50000,3.90000, after EM, it became</div>
<div>4.50000,4.50000,0.00000. so I can't use the output gro of EM for the MD run.</div>
<div> </div>
<div>Did someone encounter this before and can give me some suggestion for my mdp setup?</div>
<div> </div>
<div>Thanks in advance.</div>
<div> </div>
<div>(sorry for bad English expressions)</div>
<div> </div>
<div>some options for mdp file of mine are here:</div>
<div> </div>
<div>pbc = xy</div>
<div>rlist = 1.1 </div>
<div>rcoulomb = 1.1 <br>rvdw = 1.1 </div>
<div> </div>
<div>; electrostatics<br>coulombtype = PME </div>
<div>fourierspacing = 0.12 </div>
<div>;fourier_nx = </div>
<div>;fourier_ny = <br>fourier_nz = 30<br>pme_order = 4 <br> dim</div>
<div>ewald_geometry = 3dc ; IMPORTANT</div>
<div>optimize_fft = yes ; </div>
<div> </div>
<div>nwall = 2 <br>wall_type = 10-4 <br>;wall_r_linpot = <br>wall_atomtype = SI SI <br>wall_density = 7.12 7.12 <br>
wall_ewald_zfac = 3 </div>
<div> </div>
<div>; NEMD<br>acc_grps = WA WAb WAu<br>accelerate = 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0<br>freezegrps = WA WAb WAu<br>freezedim = Y Y Y Y Y Y Y Y Y<br>
</div>
<div>; Electric fields <br>E_x = 1 -0.55 0</div>
<div></div>
<div> </div>
<div><br clear="all">Regards,<br>Yinhe Zhang<br></div><br /><hr />Express yourself instantly with MSN Messenger! <a href='http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/' target='_new'>MSN Messenger</a></body>
</html>