<p>Hey :)</p>
<p>A little while ago I tried growing in molecules in a solution starting from shrunken coordinates. It would have been very useful if the box also had taken part in the energy minimization... That would have saved quite bit of trouble relaxing with NPT MD.</p>
<p>Just my two cents :)</p>
<p>Tsjerk</p>
<p><blockquote type="cite">On Dec 29, 2010 7:32 PM, "Justin A. Lemkul" <<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>> wrote:<br><br><p><font color="#500050">
MyLinkka wrote:
>
> Hi Justin,
>
> Thanks for your answer!
>
> As far as I understood, the virial ...</font></p>
I suppose, in theory, it's possible to use an EM algorithm to optimize a box. But for most systems, like a box filled with solvent, during EM the positions move very little and I'd be willing to bet the box wouldn't even change. I guess it depends on what you're looking to do, but for a system that hasn't even been equilibrated under the desired thermodynamic conditions (with a proper velocity distribution applied to the atoms), I don't see any real value in trying to obtain what you might consider an "optimized" box for a system that is not at equilibrium. The box stabilizes during NPT equilibration.<br>
<font color="#888888">
<br>
-Justin</font><p><font color="#500050">
> Ting Li
>
>
> On 12/25/2010 3:41 PM, Justin A. Lemkul wrote:
>>
>>
>>
>> MyLinkka wrote:
>>>
>>>...</font></p><p><font color="#500050">========================================
Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MI...</font></p><p><font color="#500050">gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a>
...</font></p></blockquote></p>