<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content=text/html;charset=iso-8859-1 http-equiv=Content-Type>
<META name=GENERATOR content="MSHTML 8.00.6001.18702"></HEAD>
<BODY style="PADDING-LEFT: 10px; PADDING-RIGHT: 10px; PADDING-TOP: 15px" 
id=MailContainerBody leftMargin=0 topMargin=0 CanvasTabStop="true" 
name="Compose message area">
<DIV><FONT face=Calibri>I don't know the conditions you use with your system, 
but sometimes there are problems when you use pressure coupling in a system with 
positionally restrained or fixed atoms.</FONT></DIV>
<DIV><FONT face=Calibri>Cheers.</FONT></DIV>
<DIV><FONT face=Calibri>&nbsp;&nbsp;&nbsp;&nbsp; Lucio.</FONT></DIV>
<DIV style="FONT: 10pt Tahoma">
<DIV><BR></DIV>
<DIV style="BACKGROUND: #f5f5f5">
<DIV style="font-color: black"><B>From:</B> <A title=ilchorny@gmail.com 
href="mailto:ilchorny@gmail.com">Joe Joe</A> </DIV>
<DIV><B>Sent:</B> Friday, April 10, 2009 12:12 PM</DIV>
<DIV><B>To:</B> <A title=gmx-users@gromacs.org 
href="mailto:gmx-users@gromacs.org">Discussion list for GROMACS users</A> </DIV>
<DIV><B>Subject:</B> Re: [gmx-users] Pressure Coupling Problem</DIV></DIV></DIV>
<DIV><BR></DIV>I finally figured it out. I went through every parameter step by 
step and it turns our I had epsilon-r set to 80. Not sure why I had that. Wish 
gromacs would have given me a warning (hint hint). That explains why my P.E. was 
10^-5 instead of 10^-6. Thanks everyone for trying!!!. Sometime the most obvious 
mistake takes some time to figure out. 
<DIV><BR></DIV>
<DIV>Cheers,<BR>
<DIV><BR></DIV>
<DIV>Ilya</DIV>
<DIV><BR>
<DIV><BR>
<DIV class=gmail_quote>On Fri, Apr 10, 2009 at 8:05 AM, <SPAN dir=ltr>&lt;<A 
title="mailto:chris.neale@utoronto.ca&#10;CTRL + clic para seguir el vínculo" 
href="mailto:chris.neale@utoronto.ca">chris.neale@utoronto.ca</A>&gt;</SPAN> 
wrote:<BR>
<BLOCKQUOTE 
style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" 
class=gmail_quote>Alright, sorry that I wasn't able to help. I'm confused by 
  some apparent contradictions in your posts and I'm not sure that I'm going to 
  be useful to you here.<BR><BR>Quoting <A 
  title="http://www.gromacs.org/pipermail/gmx-users/2009-April/041173.html&#10;CTRL + clic para seguir el vínculo" 
  href="http://www.gromacs.org/pipermail/gmx-users/2009-April/041173.html" 
  target=_blank>http://www.gromacs.org/pipermail/gmx-users/2009-April/041173.html</A>: 
  "No matter how much minimization I do the volume of the box expands when run 
  it using berendsen pressure coupling with a tau_p -f .1."<BR><BR>Quoting <A 
  href="http://www.gromacs.org/pipermail/gmx-users/2009-April/041159.html" 
  target=_blank>http://www.gromacs.org/pipermail/gmx-users/2009-April/041159.html</A>: 
  "So I got my small water box (800 waters) to behave stably with pressure 
  coupling after more minimization ..."<BR><BR>Good luck.<BR>Chris.<BR><BR>-- 
  original message --<BR><BR><BR>Nope not an A/nm problem.<BR>As a simple test I 
  take spc.gro from share/top.<BR>I reconfigure the box (i.e. editconf -f 
  spc.gro -d 1.0 -c -bt cubic -o<BR>water_center.<BR><BR>I then solvate with 
  genbox, minimize and run using the mdp file I provided<BR>earlier.<BR><BR>No 
  matter how much minimization I do the volume of the box expands when run<BR>it 
  using berendsen pressure coupling with a tau_p -f 
  .1.<BR><BR>Ilya<BR><BR><BR>On Thu, Apr 9, 2009 at 11:22 AM, Chris Neale &lt;<A 
  href="mailto:chris.neale@utoronto.ca" 
  target=_blank>chris.neale@utoronto.ca</A>&gt;wrote:<BR><BR>[Hide Quoted 
  Text]<BR>So your problem with the small water box was solved simply by adding 
  more<BR>minimization? I then suspect that all of your problems are simply 
  related to<BR>a bad starting structure -- and by the sound of it is really is 
  very bad.<BR>Are you sure that you don't have an angstrom / nm problem 
  here?<BR><BR>Chris.<BR><BR>-- original message --<BR><BR>So I got my small 
  water box (800 waters) to behave stably with pressure<BR>coupling after more 
  minimization but I still can't get my large system to<BR>work with pressure 
  coupling. I tried minimizing but I can never get the<BR>Fmax to be less 102, 
  which is pretty normal for protein/water simulations of<BR>large proteins, at 
  least from my experience. &nbsp;I have since run 400 ps NVT<BR>as 
  ...<BR><BR>_______________________________________________<BR>gmx-users 
  mailing list &nbsp; &nbsp;<A href="mailto:gmx-users@gromacs.org" 
  target=_blank>gmx-users@gromacs.org</A><BR><A 
  href="http://www.gromacs.org/mailman/listinfo/gmx-users" 
  target=_blank>http://www.gromacs.org/mailman/listinfo/gmx-users</A><BR>Please 
  search the archive at <A href="http://www.gromacs.org/search" 
  target=_blank>http://www.gromacs.org/search</A> before posting!<BR>Please 
  don't post (un)subscribe requests to the list. Use thewww interface or send it 
  to <A href="mailto:gmx-users-request@gromacs.org" 
  target=_blank>gmx-users-request@gromacs.org</A>.<BR>Can't post? Read <A 
  href="http://www.gromacs.org/mailing_lists/users.php" 
  target=_blank>http://www.gromacs.org/mailing_lists/users.php</A><BR></BLOCKQUOTE></DIV><BR></DIV></DIV></DIV>
<P>
<HR>

<P></P>_______________________________________________<BR>gmx-users mailing 
list&nbsp;&nbsp;&nbsp; 
gmx-users@gromacs.org<BR>http://www.gromacs.org/mailman/listinfo/gmx-users<BR>Please 
search the archive at http://www.gromacs.org/search before posting!<BR>Please 
don't post (un)subscribe requests to the list. Use the <BR>www interface or send 
it to gmx-users-request@gromacs.org.<BR>Can't post? Read 
http://www.gromacs.org/mailing_lists/users.php</BODY></HTML>