[gmx-users] Problem in PMF simulation
Billy Williams-Noonan
billy.williams-noonan at monash.edu
Fri Dec 2 02:55:54 CET 2016
Thanks :)
On 1 December 2016 at 23:40, Justin Lemkul <jalemkul at vt.edu> wrote:
>
>
> On 12/1/16 4:53 AM, Billy Williams-Noonan wrote:
>
>> Also, you may want to check that your rcoulomb and rvdw values are
>> correct... Didn't look at your mdp files until just now. Why are they
>> set
>> to zero?
>>
>>
> Vacuum simulations are done with infinite cutoffs (rlist=rcoulomb=rvdw=0)
> and no PBC.
>
> -Justin
>
> On 29 November 2016 at 16:56, Billy Williams-Noonan <
>> billy.williams-noonan at monash.edu> wrote:
>>
>> My PMF curves look similar to that. What you're seeing towards the end of
>>> that curve is either an artifact of the algorithms used to model the
>>> system, a very minor sampling problem, or a mix of both. Your Gaussians
>>> may be split because they are bi-modal and the most energetically
>>> favourable position around your constraint is in two places. Try having
>>> a
>>> look at your individual umbrella windows to figure out what's going on.
>>> :)
>>>
>>> Cheers,
>>>
>>> Billy
>>>
>>> On 29 November 2016 at 15:20, Shi Li <sli259 at g.uky.edu> wrote:
>>>
>>> Dear Gromacs users,
>>>>
>>>> I posted this question a few days ago, but still couldn’t solve it. This
>>>> problem has bothered my for a long time, so I am reposting this for
>>>> help.
>>>>
>>>> I am currently doing some PMF simulations for two small molecular types
>>>> in vacuum. I applied an energy minimization, a NVT equilibrium and a
>>>> production run on my system. Following are the mdp files for each. I
>>>> used a
>>>> commend in my run script to change the Window in pull_coord1_init in the
>>>> pull code for each sample window.
>>>>
>>>> Energy minimization
>>>> ---
>>>> integrator = steep
>>>> nsteps = 20000
>>>>
>>>> nstenergy = 500
>>>> nstlog = 500
>>>> nstxout-compressed = 1000
>>>>
>>>> cutoff-scheme = group
>>>>
>>>> coulombtype = Cut-off
>>>> rcoulomb = 0
>>>>
>>>> vdwtype = Cut-off
>>>> rvdw = 0
>>>> rlist = 0
>>>> pbc = no
>>>>
>>>> ; Pull code
>>>> pull = umbrella
>>>> pull_geometry = distance
>>>> pull_start = no
>>>> pull-ncoords = 1
>>>> pull_group1_name = ADT1
>>>> pull_group2_name = ADT2
>>>> pull-coord1-groups = 1 2
>>>> pull_coord1_init = WINDOW
>>>> pull_coord1_rate = 0.0
>>>> pull_coord1_k = 5000
>>>>
>>>> —
>>>> NVT
>>>> ----
>>>> integrator = md
>>>> dt = 0.002 ; 2 fs
>>>> nsteps = 500000 ; 1.0 ns
>>>>
>>>> nstenergy = 200
>>>> nstlog = 2000
>>>> nstxout-compressed = 10000
>>>>
>>>> continuation = yes
>>>> constraint-algorithm = lincs
>>>> constraints = h-bonds
>>>>
>>>> cutoff-scheme = group
>>>> rlist = 0
>>>>
>>>> rcoulomb = 0
>>>>
>>>> rvdw = 0
>>>> pbc = no
>>>> tcoupl = V-rescale
>>>> tc-grps = System
>>>> tau-t = 2.0
>>>> ref-t = 298.15
>>>> nhchainlength = 1
>>>>
>>>> ; Pull code
>>>> pull = umbrella
>>>> pull_geometry = distance
>>>> ;pull_dim = N N Y
>>>> pull_start = no
>>>> pull-ncoords = 1
>>>> pull_ngroups = 2
>>>> pull_group1_name = ADT1
>>>> pull_group2_name = ADT2
>>>> pull-coord1-groups = 1 2
>>>> pull_coord1_init = WINDOW
>>>> pull_coord1_rate = 0.0
>>>> pull_coord1_k = 5000
>>>>
>>>> ——
>>>> Production run:
>>>> ----
>>>> integrator = md
>>>> dt = 0.002 ; 2 fs
>>>> nsteps = 10000000 ; 20.0 ns
>>>>
>>>> nstenergy = 5000
>>>> nstlog = 5000
>>>> nstxout-compressed = 5000
>>>>
>>>> continuation = yes
>>>> constraint-algorithm = lincs
>>>> constraints = h-bonds
>>>>
>>>> cutoff-scheme = Group
>>>>
>>>> rcoulomb = 0
>>>> rlist = 0
>>>>
>>>> rvdw = 0
>>>> pbc = no
>>>>
>>>> tcoupl = V-rescale; Nose-Hoover
>>>> tc-grps = System
>>>> tau-t = 2.0
>>>> ref-t = 298.15
>>>> nhchainlength = 1
>>>>
>>>> ; Pull code
>>>> pull = umbrella
>>>> pull_geometry = distance
>>>> pull_start = no
>>>> pull-ncoords = 1
>>>> pull_ngroups = 2
>>>> pull_group1_name = ADT1
>>>> pull_group2_name = ADT2
>>>> pull-coord1-groups = 1 2
>>>> pull_coord1_init = WINDOW
>>>> pull_coord1_rate = 0.0
>>>> pull_coord1_k = 5000
>>>>
>>>> ———
>>>> I have a very strange profile and histogram. I uploaded to these links.
>>>>
>>>> PMF: https://www.dropbox.com/s/i10hn1p30w2j71r/F-anti-adt.eps?dl=0
>>>> HISTO: https://www.dropbox.com/s/1xi9i4cj60g97n2/F-anti-histo.eps?dl=0
>>>>
>>>> The molecule is about 1.5 nm long, so the smooth profile before 1.5 nm
>>>> looks fine, but I don’t understand why the profile had this fluctuation
>>>> after 1.5 nm and the histogram seems to have some split peaks after this
>>>> length. Can anyone tell me why this happened? Is there something wrong
>>>> in
>>>> my mdp file or something wrong in my topology? I have used this similar
>>>> mdp
>>>> files for my other types of molecular systems and the results were
>>>> fine. So
>>>> it is really confused me now.
>>>>
>>>> I appreciate for any help! Thanks!
>>>>
>>>> Shi
>>>>
>>>>
>>>>
>>>> --
>>>> 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.
>>>>
>>>
>>>
>>>
>>>
>>> --
>>> Billy Noonan* | *PhD Student *|* Bsci ( *Adv* ), IA Hon
>>>
>>> *LinkedIn Profile
>>> <http://www.linkedin.com/profile/preview?locale=en_US&trk=
>>> prof-0-sb-preview-primary-button>
>>> **|* +61420 382 557 <+61%20420%20382%20557>
>>>
>>> Monash Institute for Pharmaceutical Sciences ( *MIPS* )
>>> Royal Parade, Parkville, 3052
>>>
>>>
>>>
>>
>>
> --
> ==================================================
>
> Justin A. Lemkul, Ph.D.
> Ruth L. Kirschstein NRSA Postdoctoral Fellow
>
> Department of Pharmaceutical Sciences
> School of Pharmacy
> Health Sciences Facility II, Room 629
> University of Maryland, Baltimore
> 20 Penn St.
> Baltimore, MD 21201
>
> jalemkul at outerbanks.umaryland.edu | (410) 706-7441
> http://mackerell.umaryland.edu/~jalemkul
>
> ==================================================
>
> --
> 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.
>
--
Billy Noonan* | *PhD Student *|* Bsci ( *Adv* ), IA Hon
*LinkedIn Profile
<http://www.linkedin.com/profile/preview?locale=en_US&trk=prof-0-sb-preview-primary-button>
**|* +61420 382 557
Monash Institute for Pharmaceutical Sciences ( *MIPS* )
Royal Parade, Parkville, 3052
More information about the gromacs.org_gmx-users
mailing list