<DIV><a style="color: blue; background-color: ffffdd"
href="http://nospammer.net/SubmitSpam/submitspam.php?id=I3RWAlEXBhJdXFxffkN1
Xg__&sig=e3g2aS5UljdoIzFmaiVlc3l4NCx2MTN7AkJdCQdOEB05DkACE1VGdGgjM2xqJWVxeXg
hSxtJHgEbRkMfKkIXHhQZV19YXkETaCMxbmolZXF5eDRqHkVHBFIMHgIFVhUQFBVRXlhfVgBHcEQ
OB0wRIgkZWQMFRFEZAVdCHAtISwERCAtFEgx6Rzp0cD4FUgApIDlwRx5XVR8mEjFsaidlcXl4NCx
2b1sAHFMLQ0VIBB8MGVgCEUNcGQlAQkIFVwJeHRdXWRtUXQAJV1gDBAoGBAsKUUICHkYHDVEcCx9
MARRWFVBcW15DAAFMXx9ETREcFVtaQxgcVgUdSl0FCFcMBBRVWUh2MTN2aCMxbWolZWoRDEBcTB4
cAx9UHxgCQAkUFBNBQBpQUVoLTFxsaiVncXl4NCx2MQgcHFdBVkUKEgYOVlNeGVxSFxsNXh4NCjY
ECQhbXgIefhUBT1gCDXopGAoMRwMxfGtZPVBUHhl6KRgKDDQsdjMzdGgjMWxqCQ0FDQgOA1lGRAN
GREMDB0QGAlcXRktZYkYEGExDGEVoBBgVEVpLKX1aBxxQMWxqJ2VxeXg0LHZzWwAcU0JWRQoIEBA
UWEUFRR0HEVAfBx5NSwIcV1lNH11eFQYMXQUZUQwfHxcbSwReXhULUB8DGEI6FhQAGVkFVEEHaCM
xbmolZXF5eDQxEVxLWR1QVB4ZCBcUCA1RXwJxVAYHTlAPGQsKAx54NCxyMTN0aCMxbHIlZXF7eDQ
nAbOD460TYaxz5rvPHU6aKBw_">
If this email is not spam, click here to submit the signatures to 
FortiGuard - AntiSpam Service.
</a></DIV>

<PRE>  
Dear Justin,
Thanks a lot for your response.
Neda

----- Original Message -----
From: Justin Lemkul <jalemkul@vt.edu>
To: gmx-users@gromacs.org
Sent: Sun, 15 Jul 2018 23:09:11 +0430 (IRDT)
Subject: [Spam-FG] Re: [gmx-users] Negative pressure in NVT simulation



<DIV><a style="color: blue; background-color: ffffdd" href="http://nospammer.net/SubmitSpam/submitspam.php?id=I3RWARoweXYADkhffkN1
WA__&sig=e3g1*tJ118poIzFkaiVlc3l4NCx2MTMqAFdFHFAKShwYFkFNGh9UBgdOUA8ZCwoDHld
QQxVEXhEGV1AYA0oLXhoNRl4TX0dbHVBUHkdCEBgdHRtBEkEeGxhXWAMEVksZDRVYDxheXVkNUkQ
FBkwHAxANWQEbVTN0aCExbGolZXF5d15NGlReHx1PcRoeCwAVDHg0LHQxM3RpIzFscU0RBQlCGwM
BRkRaHEtUAA9IDgQVFFVOWFJcGWgjMW5qJWVxeXg0Fx5FRwRSDB4bHVJLFgsXWU0VQh0bGkQePx9
VFR4LDBthF1hfHQZEbiADVhECVj95dFtkQBEaUG4gA1YRcXl4Nix2MTN0aCMdBB5RFUtWV0NbAR9
UBgdOUA8ZCwoDHldnWQZBXAYcDHwNA0kMHx4neEUFRUB0aCMzbGolZXF5eHZEAkVDB1IMHgELTAk
dEAtAAgVIQFoDV1lCGUBKHBgRWEEXXxwYAVBFBQRDCl4eCltBF1JAWgdRVjMNSB1cDAtRXgUxM3R
qIzFsaiVlcWQfWVRbREARGlAcHg9UEBQKDHRLBF5eFQtQHwMYQmVxeXw0LHYxM3RoOzFsaidlcX*
QXtr1BsjkCP*XyL2dAAIh5g__" target="_blank">
If this email is not spam, click here to submit the signatures to 
FortiGuard - AntiSpam Service.
</a></DIV>
<PRE>  
On 7/15/18 2:32 PM, Neda Rafiee wrote:
&gt; ;
&gt; energygrps          = LYS LEU VAL PHE1 PHE2 ALA GLU ASP NaB ClB OXB OSil OGem OMB SIB SiSil SiGem HSil HGem OW HW
&gt; energygrp_table = NaB OXB NaB OSil NaB OGem NaB OMB OMB OW OW SIB OW SiSil OW SiGem OW HSil OW HGem HW OMB HW OXB OW OSil OW OGem HW OSil HW OGem OXB OW
&gt;
&gt; freezegrps          = SIB OXB OMB SiSil OSil SiGem OGem
&gt; freezedim           = Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y

You&#39;re not using energy group exclusions, which the manual recommends 
when freezing:

&quot;To avoid spurious contributions to the virial and pressure due to large 
forces between completely frozen atoms you need to use energy group 
exclusions, this also saves computing time.&quot;

http://manual.gromacs.org/documentation/current/user-guide/mdp-options.html#non-equilibrium-md

Freezing is a very severe artificial perturbation, so you will have to 
use hacks like that if you intend to do a simulation with so many things 
frozen.

-Justin

-- 
==================================================

Justin A. Lemkul, Ph.D.
Assistant Professor
Virginia Tech Department of Biochemistry

303 Engel Hall
340 West Campus Dr.
Blacksburg, VA 24061

jalemkul@vt.edu | (540) 231-3129
http://www.thelemkullab.com

==================================================

-- 
Gromacs Users mailing list

* Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

* Can&#39;t post? Read http://www.gromacs.org/Support/Mailing_Lists

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request@gromacs.org.