[gmx-users] [Spam-FG] [Spam-FG]__Re:_[Spam-FG]_Re:__Negative_pressure_in_NVT_simulation

Neda Rafiee nerafiee at ipm.ir
Mon Jul 16 14:55:50 CEST 2018 

Dear Mark,
I removed different residues from my energy group list, but still the problem of very large negative pressure exist. Any other solution?
Thanks in advance.
Neda


----- Original Message -----
From: Neda Rafiee <nerafiee at ipm.ir>
To: gmx-users at gromacs.org
Sent: Mon, 16 Jul 2018 12:03:34 +0430 (IRDT)
Subject: Re: [gmx-users]	[Spam-FG]	[Spam-FG]__Re:_[Spam-FG]_Re:__Negative_pressure_in_NVT_simulation

Dear Mark,
Thanks a lot for your useful suggestions
Neda


----- Original Message -----
From: Mark Abraham <mark.j.abraham at gmail.com>
To: gmx-users at gromacs.org
Sent: Mon, 16 Jul 2018 11:49:38 +0430 (IRDT)
Subject: Re: [gmx-users] [Spam-FG]	[Spam-FG]__Re:_[Spam-FG]_Re:__Negative_pressure_in_NVT_simulation

Hi,

I presume that's because you still have interactions between groups that
are frozen.

You can make your life easier by using only the groups you need for your
simulation, e.g. those needing a table and those that are frozen. If you
are trying to get per-residue decomposition, then (first) make sure your
force field was parametrized such that this is meaningful (most aren't),
and prefer to do that analysis using mdrun -rerun. mdrun goes slower for
each energy group that has interactions to compute, so minimize that.

I would also have done a simulation with just my slab and found this
pressure issue back then. It is good practice to add complexity to models
incrementally, after validating your expectations on the simulation at the
previous stage.

Mark

On Mon, Jul 16, 2018 at 9:07 AM Neda Rafiee <nerafiee at ipm.ir> wrote:

> Dear Justin,
> I added the following energy exclusion pairs to my mdp file. Now potential
> energy decreased significantly but still I have a very large negative
> pressure (minus 10 to power 28) in my system. It seems that pressure is the
> same as before.Should I do anything else?
>
> energygrp-excl= SIB SIB  SIB OXB  SIB OMB  SIB SiSil  SIB OSil  SIB SiGem
> SIB OGem  OXB OXB  OXB OMB  OXB SiSil  OXB OSil  OXB SiGem  OXB OGem  OMB
> OMB  OMB SiSil  OMB OSil  OMB SiGem  OMB OGem  SiSil SiSil  SiSil OSil
> SiSil SiGem  SiSil OGem  OSil OSil  OSil SiGem  OSil OGem  SiGem SiGem
> SiGem OGem  OGem OGem
>
> Thanks
> Neda
>
> ----- Original Message -----
> From: Neda Rafiee <nerafiee at ipm.ir>
> To: gmx-users at gromacs.org
> Sent: Mon, 16 Jul 2018 00:01:35 +0430 (IRDT)
> Subject: [Spam-FG][gmx-users]
> [Spam-FG]__Re:_[Spam-FG]_Re:__Negative_pressure_in_NVT_simulation
>
> <DIV><a style="color: blue; background-color: ffffdd" href="
> http://nospammer.net/SubmitSpam/submitspam.php?id=I3RWR3EnBQVzelBffkN1
>
> VA__&sig=e3gwAc6BBLhoIzFgaiVlc3l4NCx2MTOLAFdFHFAKSh8WC0RNG1xWBkZNVBhFdhATFBF
>
> AfwZQXlsbVlMBA1EWARgVGlweQQwdDB54XzhyJB08IHZEPFVrMhBFVwckFGh7IR9rc1BCWhNVRgI
> <http://nospammer.net/SubmitSpam/submitspam.php?id=I3RWR3EnBQVzelBffkN1VA__&sig=e3gwAc6BBLhoIzFgaiVlc3l4NCx2MTOLAFdFHFAKSh8WC0RNG1xWBkZNVBhFdhATFBFAfwZQXlsbVlMBA1EWARgVGlweQQwdDB54XzhyJB08IHZEPFVrMhBFVwckFGh7IR9rc1BCWhNVRgI>
>
> LWEQ2RCwUXkgZeEkyBUJYOgZGVh1NNndURHxnOl9iWiYOZjQVNj12HEN1WDUrZgA6LUEiFhM1BlQ
>
> He2QiEEZpC2cvDSIBDH5kEXRRJgNuVycBbD05ESluekcIaiwDZmUNKWgdExQXWHYudwYRLHFAJAF
>
> zLTM*MXlkEXh1IgB2YCooczcpFRBSehFzexctcjxmJWdVBisxU0csZmI5LnF3KTBkMxU6MHVYP2J
>
> EdGgjM2xqJWVxeXg7RhddVhkDVl0sHFFLFB0NNCx2MzN0aCMxbGowAhwBVUFfE0NANA9RXgELRhZ
>
> fFgpTLHYxMXRoIzFsaiWaGQ0MRBZZHl0bG1NQAQdAF18XHUADJURRGQFXYhwLSEoCDBpZRQJCQxU
>
> FDUEEGhoMFUQxB34hcGEbH0ZpNSthDhkfHl9iRzw5Iyl8bkoZTAJMHEtTHVxFeUVZG0EDI18jGhg
>
> RYkAVAl9AJmBJXidxKAA4PlBqPndyPztLRjUsTiM*PhANeTRWVzs9Ygk2KVIKNTEUUCF8YGIMPmZ
> pBC9iM0A4IXUcGX1rHAdtYwBecT1BHRp8ejRkex8MYHQuDUEtIw06cUczVHQMAHtmLSdgMxoKIn
>
> B.IGh3DABGaToBazAaKHU.ajRaRDwpW3AiPXQgEy8sehwXcnYMCmReADB9ZXF5ejQsdjEzdGh9WR
>
> geVV9eVhVVQgNQX1oPUV4BC0YWXxYKUwMSXlABBUZfGAtRDB4XV1dZBENWGhwMRB8PV0gWDBFQSV
>
> lcVwRFTEEYA0oLAlcQQEEaEl0bBg5UHR9MCRgbCl1ZGxxeEGgjMW5qJWVweXg0Nx5FRwRSDB4bHV
>
> JLBREdWEkbWkYYBEJTQglKCHF5eDYsdjEzdGgjCgQeURVLVldDWwEfVAYHTlAPGQsKAx5XZ1kGQV
>
> wGHAx8DQNJDB8eJ3hFBUVAWy9uaUE*VgADCid4RQVFM3RoITFsaiVlcXlUXFgCQQlbR1RGG0RCFx
>
> 4UGVdfWF5BE0dwRBwaShcFVjVVRRpYXRM3b1gfHlZlcXl6NCx2MTN0aGFZGB5VFktWV1lNH11fHR
>
> tXHx8TVksaDRAaXxMeXhUBT1wNBAoJGAoMXUIQXhwTGkxcDQlWSx4LH2tLG0keARtGQx9qJWVzeX
>
> g0LHYxM2kPTklBH1YAAwpVRkkHRFYHHGNWHgVIBBIKVlteETEzdGwjMWxqJWVxYXg0LHcxM3v0yu
> ILmFi9OA8Sb.IVubYVMyMxbG4lZXF5eDQsbjEzdGojMW0HyGY2m9p1eztByEYRoCDQ3g__"
> target="_blank">
> If this email is not spam, click here to submit the signatures to
> FortiGuard - AntiSpam Service.
> </a></DIV>
> <DIV><a style="color: blue; background-color: ffffdd" href="
> http://nospammer.net/SubmitSpam/submitspam.php?id=I3RWAlEXBhJdXFxffkN1
>
> Xg__&sig=e3g2aS5UljdoIzFmaiVlc3l4NCx2MTN7AkJdCQdOEB05DkACE1VGdGgjM2xqJWVxeXg
>
> hSxtJHgEbRkMfKkIXHhQZV19YXkETaCMxbmolZXF5eDRqHkVHBFIMHgIFVhUQFBVRXlhfVgBHcEQ
> <http://nospammer.net/SubmitSpam/submitspam.php?id=I3RWAlEXBhJdXFxffkN1Xg__&sig=e3g2aS5UljdoIzFmaiVlc3l4NCx2MTN7AkJdCQdOEB05DkACE1VGdGgjM2xqJWVxeXghSxtJHgEbRkMfKkIXHhQZV19YXkETaCMxbmolZXF5eDRqHkVHBFIMHgIFVhUQFBVRXlhfVgBHcEQ>
>
> OB0wRIgkZWQMFRFEZAVdCHAtISwERCAtFEgx6Rzp0cD4FUgApIDlwRx5XVR8mEjFsaidlcXl4NCx
>
> 2b1sAHFMLQ0VIBB8MGVgCEUNcGQlAQkIFVwJeHRdXWRtUXQAJV1gDBAoGBAsKUUICHkYHDVEcCx9
>
> MARRWFVBcW15DAAFMXx9ETREcFVtaQxgcVgUdSl0FCFcMBBRVWUh2MTN2aCMxbWolZWoRDEBcTB4
>
> cAx9UHxgCQAkUFBNBQBpQUVoLTFxsaiVncXl4NCx2MQgcHFdBVkUKEgYOVlNeGVxSFxsNXh4NCjY
>
> ECQhbXgIefhUBT1gCDXopGAoMRwMxfGtZPVBUHhl6KRgKDDQsdjMzdGgjMWxqCQ0FDQgOA1lGRAN
>
> GREMDB0QGAlcXRktZYkYEGExDGEVoBBgVEVpLKX1aBxxQMWxqJ2VxeXg0LHZzWwAcU0JWRQoIEBA
>
> UWEUFRR0HEVAfBx5NSwIcV1lNH11eFQYMXQUZUQwfHxcbSwReXhULUB8DGEI6FhQAGVkFVEEHaCM
>
> xbmolZXF5eDQxEVxLWR1QVB4ZCBcUCA1RXwJxVAYHTlAPGQsKAx54NCxyMTN0aCMxbHIlZXF7eDQ
> nAbOD460TYaxz5rvPHU6aKBw_" target="_blank">
> 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
> To: gmx-users at 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
> <http://nospammer.net/SubmitSpam/submitspam.php?id=I3RWARoweXYADkhffkN1WA__&sig=e3g1*tJ118poIzFkaiVlc3l4NCx2MTMqAFdFHFAKShwYFkFNGh9UBgdOUA8ZCwoDHldQQxVEXhEGV1AYA0oLXhoNRl4TX0dbHVBUHkdCEBgdHRtBEkEeGxhXWAMEVksZDRVYDxheXVkNUkQ>
>
> 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 at 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 at gromacs.org.
>
>
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can'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 at gromacs.org.
>
-- 
Gromacs Users mailing list

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

* Can'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 at gromacs.org.

-- 
Gromacs Users mailing list

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

* Can'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 at gromacs.org.

More information about the gromacs.org_gmx-users mailing list