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

Mark Abraham mark.j.abraham at gmail.com
Mon Jul 16 17:13:51 CEST 2018


Hi,

No idea sorry. I never attempt anything like this, in part because frozen
groups are anyway unphysical.

Mark

On Mon, Jul 16, 2018 at 2:56 PM Neda Rafiee <nerafiee at ipm.ir> wrote:

> 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 thin
> <https://maps.google.com/?q=tend+to+do+a+simulation+with+so+many+thin&entry=gmail&source=g>
> gs
> > 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.
>
> --
> 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