<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 would expect nearly any cell to deform to a very elongated shape,<br>unless it contains a solid (crystal). In most cases an elongated<br>shape has the lowest (electrostatic) energy. Therefore anisotropic<br>coupling is only useful for crystal simulations.<br><br>Berk<br><br>&gt; Date: Mon, 1 Mar 2010 10:09:08 +0100<br>&gt; From: spoel@xray.bmc.uu.se<br>&gt; To: gmx-users@gromacs.org<br>&gt; Subject: Re: [gmx-users] Constraining periodic cell<br>&gt; <br>&gt; On 2010-03-01 09.49, David Chalmers wrote:<br>&gt; &gt; David,<br>&gt; &gt;<br>&gt; &gt; I get a different simulation result if I use isotropic or anisotropic cells.<br>&gt; &gt; I am not sure yet which is (more) correct, but I would like to be able to<br>&gt; &gt; run both types to see.  I understand that a constrained anisotropic cell is<br>&gt; &gt; not completely anisotropic, but a fully anistropic one clearly doesn't work.<br>&gt; &gt;<br>&gt; &gt; David<br>&gt; <br>&gt; You did not say which algorithm you use either, but making the time <br>&gt; constant longer might help. However, depending on your system, you may <br>&gt; need substantial surface tension to keep the system neutral if you do <br>&gt; use anisotropic coupling. We do not completely understand these effects <br>&gt; yet, see e.g http://bugzilla.gromacs.org/show_bug.cgi?id=165 for <br>&gt; additional discussions.<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;&gt;&gt; Dear All,<br>&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt; We are running some simulations using an anisotropic periodic cell.  We are<br>&gt; &gt;&gt;&gt; getting 'cell runaway' with the cell becoming very long and thin.  This<br>&gt; &gt;&gt;&gt; appears to being driven by the electrostatics of the system.  The simulation<br>&gt; &gt;&gt;&gt; then dies because the smallest cell dimension is less than 2*cutoff.<br>&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt; Is there a way that we can apply some constraints to the periodic cell?<br>&gt; &gt;&gt; Why not use isotropic pressure coupling?<br>&gt; &gt;&gt;<br>&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt; Regards<br>&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt; David<br>&gt; &gt;<br>&gt; &gt; ________________________________________________________________________<br>&gt; &gt;<br>&gt; &gt; David Chalmers                                      Lab: 9903 9110<br>&gt; &gt; Faculty of Pharmacy, Monash University<br>&gt; &gt; 381 Royal Pde, Parkville, Vic 3053. Australia<br>&gt; &gt; David.Chalmers@pharm.monash.edu.au<br>&gt; &gt; ________________________________________________________________________<br>&gt; &gt;<br>&gt; &gt;<br>&gt; <br>&gt; <br>&gt; -- <br>&gt; David van der Spoel, Ph.D., Professor of Biology<br>&gt; Dept. of Cell &amp; Molec. Biol., Uppsala University.<br>&gt; Box 596, 75124 Uppsala, Sweden. Phone:        +46184714205.<br>&gt; spoel@xray.bmc.uu.se    http://folding.bmc.uu.se<br>&gt; -- <br>&gt; gmx-users mailing list    gmx-users@gromacs.org<br>&gt; http://lists.gromacs.org/mailman/listinfo/gmx-users<br>&gt; Please search the archive at http://www.gromacs.org/search before posting!<br>&gt; Please don't post (un)subscribe requests to the list. Use the <br>&gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br>                                               <br /><hr />New Windows 7: Find the right PC for you. <a href='http://windows.microsoft.com/shop' target='_new'>Learn more.</a></body>
</html>