[gmx-users] umbrella sampling
Rose
rose.rhmn93 at gmail.com
Sun Feb 4 22:18:41 CET 2018
Sorry nothing mysterious,i'm rose(rose.rhmn93) in these all conversations, but i have some problems in my gmail which i can't open my emails every times(i don't why but it's sth related to country ip&...) because Of that i have to use my roomate account,sorrrry again
Sent from my iPhone
> On Feb 5, 2018, at 0:35, Justin Lemkul <jalemkul at vt.edu> wrote:
>
>
> I don't know whether you're intentionally changing your email address or display name, but please refrain from doing so as it is very confusing who I'm talking to.
>
>> On 2/4/18 4:00 PM, Nick Johans wrote:
>> maybe this is irrelative but i expeienced in pulling that when i decrease
>> dt from 2fs to 1fs i could make AA much closer to Sheet (0.1 nm)in
>> comparison with 2fs(max 0.4 nm). it's normal but these pulling runs were
>> crushed before(when dt=2fs) so i had had to decrease nsteps or decrease
>> pull_rate , but i don't face with these crush anymore! i can pull it too
>> close. i know maybe these beginner confessions seem unacceptable ;)))) and
>> silly.
>> i know this question is not accurate and can have many answers but, as a
>> professional and experienced person, would you please tell me, whenever you
>> face with a crushed simulation, what will you check?doubt? in your
>> system(separate of anything related to cpu&.... computers limitations) .
>> what cause simulations to be crushed usually?
>
> If your simulation is crashing, you're doing something inappropriate, by causing particles to collide with one another. This is an entirely separate issue from the relationship between constraints and time step, dt. Your simulations are not unstable because of excessively fast bond vibrations, they're unstable because the nonbonded interactions in your system are unfavorable when allowed to run with dt = 2 fs.
>
> If the system is unstable at any point, you should verify the parametrization of any non-standard species, which in your case is the sheet. I recall long ago telling you this would be a challenging entity to parametrize, yet you quickly obtained a model and began using it. What verification did you perform to assure yourself that interactions of this ZnS sheet are reasonable in the context of protein moieties? Did you perform any QM interaction energies? Are there experimentally known binding free energies for model compounds? How do you know this model is appropriate?
>
> These are all things you need to consider, and questions that are not within the scope of this mailing list.
>
> -Justin
>
>>> On Mon, Feb 5, 2018 at 12:05 AM, Justin Lemkul <jalemkul at vt.edu> wrote:
>>>
>>>
>>>> On 2/4/18 3:25 PM, Rose wrote:
>>>>
>>>> Jeah,i meant grompp :)))
>>>> and last question, i set constraint = h-Bonds, and dt=0.002 because Of
>>>> that(since usually h-Bonds have Max frequency and.... if i assume them
>>>> constraint so it's not reasonable to set dt=0.001 ) , is this assumption
>>>> true?
>>> You can always decrease dt, that's "reasonable," but it makes no sense and
>>> is inefficient. Without any form of constraints, you'd likely be limited to
>>> dt = 0.0005 in many cases.
>>>
>>>
>>> -Justin
>>>
>>> Sent from my iPhone
>>>>> On Feb 4, 2018, at 23:30, Justin Lemkul <jalemkul at vt.edu> wrote:
>>>>>
>>>>>
>>>>>> On 2/4/18 1:56 PM, Rose wrote:
>>>>>> Thank you so much,
>>>>>> Is it difference between using gmx tpbconv -extend 5ps OR use .gro
>>>>>> file of last simulation as an input For new 5nS
>>>>>> simulation?(continiuation=yes in .mdp file)
>>>>> The only reason to invoke grompp is if you're changing something about
>>>>> the system. For a simple continuation use tpbconv/convert-tpr and mdrun
>>>>> -cpi.
>>>>>
>>>>> -Justin
>>>>>
>>>>> Sent from my iPhone
>>>>>>> On Feb 4, 2018, at 17:00, Justin Lemkul <jalemkul at vt.edu> wrote:
>>>>>>>
>>>>>>>
>>>>>>>> On 2/2/18 10:44 AM, Rose wrote:
>>>>>>>> Sent from my iPhone
>>>>>>>>
>>>>>>>>> On Feb 2, 2018, at 19:01, Justin Lemkul <jalemkul at vt.edu> wrote:
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>> On 2/2/18 8:43 AM, rose rahmani wrote:
>>>>>>>>>> Hi, i'm doing umbrella sampling.(pulling AA toward nanosheet).and
>>>>>>>>>> 4nS
>>>>>>>>>> simulation for each window.
>>>>>>>>>>
>>>>>>>>>> g_energy -f umbrella0.edr -o out.xvg
>>>>>>>>>>> 1 Bond 2 Angle 3 Proper-Dih. 4
>>>>>>>>>> Improper-Dih.
>>>>>>>>>> 5 LJ-14 6 Coulomb-14 7 LJ-(SR) 8
>>>>>>>>>> Coulomb-(SR)
>>>>>>>>>> 9 Coul.-recip. 10 COM-Pull-En. 11 Potential 12
>>>>>>>>>> Kinetic-En.
>>>>>>>>>> 13 Total-Energy 14 Conserved-En. 15 Temperature 16
>>>>>>>>>> Pressure
>>>>>>>>>>
>>>>>>>>>> 17 Constr.-rmsd 18 Vir-XX 19 Vir-XY 20
>>>>>>>>>> Vir-XZ
>>>>>>>>>>
>>>>>>>>>> 21 Vir-YX 22 Vir-YY 23 Vir-YZ 24
>>>>>>>>>> Vir-ZX
>>>>>>>>>>
>>>>>>>>>> 25 Vir-ZY 26 Vir-ZZ 27 Pres-XX 28
>>>>>>>>>> Pres-XY
>>>>>>>>>>
>>>>>>>>>> 29 Pres-XZ 30 Pres-YX 31 Pres-YY 32
>>>>>>>>>> Pres-YZ
>>>>>>>>>>
>>>>>>>>>> 33 Pres-ZX 34 Pres-ZY 35 Pres-ZZ 36
>>>>>>>>>> #Surf*SurfTen
>>>>>>>>>> 37 Mu-X 38 Mu-Y
>>>>>>>>>>
>>>>>>>>>> 39 Mu-Z 40 Coul-SR:SOL-SOL
>>>>>>>>>>
>>>>>>>>>> 41 LJ-SR:SOL-SOL 42 Coul-14:SOL-SOL
>>>>>>>>>>
>>>>>>>>>> 43 LJ-14:SOL-SOL 44 Coul-SR:SOL-WAL
>>>>>>>>>>
>>>>>>>>>> 45 LJ-SR:SOL-WAL 46 Coul-14:SOL-WAL
>>>>>>>>>>
>>>>>>>>>> 47 LJ-14:SOL-WAL 48 Coul-SR:SOL-ZnS
>>>>>>>>>>
>>>>>>>>>> 49 LJ-SR:SOL-ZnS 50 Coul-14:SOL-ZnS
>>>>>>>>>>
>>>>>>>>>> 51 LJ-14:SOL-ZnS 52 Coul-SR:SOL-Protein
>>>>>>>>>>
>>>>>>>>>> 53 LJ-SR:SOL-Protein 54 Coul-14:SOL-Protein
>>>>>>>>>>
>>>>>>>>>> 55 LJ-14:SOL-Protein 56 Coul-SR:SOL-NA
>>>>>>>>>>
>>>>>>>>>> 57 LJ-SR:SOL-NA 58 Coul-14:SOL-NA 59 LJ-14:SOL-NA 60
>>>>>>>>>> Coul-SR:SOL-CL
>>>>>>>>>> 61 LJ-SR:SOL-CL 62 Coul-14:SOL-CL
>>>>>>>>>>
>>>>>>>>>> 63 LJ-14:SOL-CL 64 Coul-SR:SOL-wall0
>>>>>>>>>>
>>>>>>>>>> 65 LJ-SR:SOL-wall0 66 Coul-14:SOL-wall0
>>>>>>>>>>
>>>>>>>>>> 67 LJ-14:SOL-wall0 68 Coul-SR:SOL-wall1
>>>>>>>>>>
>>>>>>>>>> 69 LJ-SR:SOL-wall1 70 Coul-14:SOL-wall1
>>>>>>>>>>
>>>>>>>>>> 71 LJ-14:SOL-wall1 72 Coul-SR:WAL-WAL
>>>>>>>>>>
>>>>>>>>>> 73 LJ-SR:WAL-WAL 74 Coul-14:WAL-WAL
>>>>>>>>>>
>>>>>>>>>> 75 LJ-14:WAL-WAL 76 Coul-SR:WAL-ZnS
>>>>>>>>>>
>>>>>>>>>> 77 LJ-SR:WAL-ZnS 78 Coul-14:WAL-ZnS
>>>>>>>>>>
>>>>>>>>>> 79 LJ-14:WAL-ZnS 80 Coul-SR:WAL-Protein
>>>>>>>>>>
>>>>>>>>>> 81 LJ-SR:WAL-Protein 82 Coul-14:WAL-Protein
>>>>>>>>>>
>>>>>>>>>> 83 LJ-14:WAL-Protein 84 Coul-SR:WAL-NA
>>>>>>>>>>
>>>>>>>>>> 85 LJ-SR:WAL-NA 86 Coul-14:WAL-NA 87 LJ-14:WAL-NA 88
>>>>>>>>>> Coul-SR:WAL-CL
>>>>>>>>>> 89 LJ-SR:WAL-CL 90 Coul-14:WAL-CL
>>>>>>>>>>
>>>>>>>>>> 91 LJ-14:WAL-CL 92 Coul-SR:WAL-wall0
>>>>>>>>>>
>>>>>>>>>> 93 LJ-SR:WAL-wall0 94 Coul-14:WAL-wall0
>>>>>>>>>>
>>>>>>>>>> 95 LJ-14:WAL-wall0 96 Coul-SR:WAL-wall1
>>>>>>>>>>
>>>>>>>>>> 97 LJ-SR:WAL-wall1 98 Coul-14:WAL-wall1
>>>>>>>>>>
>>>>>>>>>> 99 LJ-14:WAL-wall1 100 Coul-SR:ZnS-ZnS
>>>>>>>>>>
>>>>>>>>>> 101 LJ-SR:ZnS-ZnS 102 Coul-14:ZnS-ZnS
>>>>>>>>>>
>>>>>>>>>> 103 LJ-14:ZnS-ZnS 104 Coul-SR:ZnS-Protein
>>>>>>>>>>
>>>>>>>>>> 105 LJ-SR:ZnS-Protein 106 Coul-14:ZnS-Protein
>>>>>>>>>>
>>>>>>>>>> 107 LJ-14:ZnS-Protein 108 Coul-SR:ZnS-NA
>>>>>>>>>>
>>>>>>>>>> 109 LJ-SR:ZnS-NA 110 Coul-14:ZnS-NA 111 LJ-14:ZnS-NA 112
>>>>>>>>>> Coul-SR:ZnS-CL
>>>>>>>>>> 113 LJ-SR:ZnS-CL 114 Coul-14:ZnS-CL
>>>>>>>>>>
>>>>>>>>>> 115 LJ-14:ZnS-CL 116 Coul-SR:ZnS-wall0
>>>>>>>>>>
>>>>>>>>>> 117 LJ-SR:ZnS-wall0 118 Coul-14:ZnS-wall0
>>>>>>>>>>
>>>>>>>>>> 119 LJ-14:ZnS-wall0 120 Coul-SR:ZnS-wall1
>>>>>>>>>>
>>>>>>>>>> 121 LJ-SR:ZnS-wall1 122 Coul-14:ZnS-wall1
>>>>>>>>>>
>>>>>>>>>> 123 LJ-14:ZnS-wall1 124 Coul-SR:Protein-Protein
>>>>>>>>>>
>>>>>>>>>> 125 LJ-SR:Protein-Protein 126 Coul-14:Protein-Protein
>>>>>>>>>>
>>>>>>>>>> 127 LJ-14:Protein-Protein 128 Coul-SR:Protein-NA
>>>>>>>>>>
>>>>>>>>>> 129 LJ-SR:Protein-NA 130 Coul-14:Protein-NA
>>>>>>>>>>
>>>>>>>>>> 131 LJ-14:Protein-NA 132 Coul-SR:Protein-CL
>>>>>>>>>>
>>>>>>>>>> 133 LJ-SR:Protein-CL 134 Coul-14:Protein-CL
>>>>>>>>>>
>>>>>>>>>> 135 LJ-14:Protein-CL 136 Coul-SR:Protein-wall0
>>>>>>>>>>
>>>>>>>>>> 137 LJ-SR:Protein-wall0 138 Coul-14:Protein-wall0
>>>>>>>>>>
>>>>>>>>>> 139 LJ-14:Protein-wall0 140 Coul-SR:Protein-wall1
>>>>>>>>>>
>>>>>>>>>> 141 LJ-SR:Protein-wall1 142 Coul-14:Protein-wall1
>>>>>>>>>>
>>>>>>>>>> 143 LJ-14:Protein-wall1 144 Coul-SR:NA-NA
>>>>>>>>>>
>>>>>>>>>> 145 LJ-SR:NA-NA 146 Coul-14:NA-NA 147 LJ-14:NA-NA 148
>>>>>>>>>> Coul-SR:NA-CL
>>>>>>>>>> 149 LJ-SR:NA-CL 150 Coul-14:NA-CL
>>>>>>>>>>
>>>>>>>>>> 151 LJ-14:NA-CL 152 Coul-SR:NA-wall0
>>>>>>>>>>
>>>>>>>>>> 153 LJ-SR:NA-wall0 154 Coul-14:NA-wall0
>>>>>>>>>>
>>>>>>>>>> 155 LJ-14:NA-wall0 156 Coul-SR:NA-wall1
>>>>>>>>>>
>>>>>>>>>> 157 LJ-SR:NA-wall1 158 Coul-14:NA-wall1
>>>>>>>>>>
>>>>>>>>>> 159 LJ-14:NA-wall1 160 Coul-SR:CL-CL 161 LJ-SR:CL-CL 162
>>>>>>>>>> Coul-14:CL-CL
>>>>>>>>>> 163 LJ-14:CL-CL 164 Coul-SR:CL-wall0
>>>>>>>>>>
>>>>>>>>>> 165 LJ-SR:CL-wall0 166 Coul-14:CL-wall0
>>>>>>>>>>
>>>>>>>>>> 167 LJ-14:CL-wall0 168 Coul-SR:CL-wall1
>>>>>>>>>>
>>>>>>>>>> 169 LJ-SR:CL-wall1 170 Coul-14:CL-wall1
>>>>>>>>>>
>>>>>>>>>> 171 LJ-14:CL-wall1 172 Coul-SR:wall0-wall0
>>>>>>>>>>
>>>>>>>>>> 173 LJ-SR:wall0-wall0 174 Coul-14:wall0-wall0
>>>>>>>>>>
>>>>>>>>>> 175 LJ-14:wall0-wall0 176 Coul-SR:wall0-wall1
>>>>>>>>>>
>>>>>>>>>> 177 LJ-SR:wall0-wall1 178 Coul-14:wall0-wall1
>>>>>>>>>>
>>>>>>>>>> 179 LJ-14:wall0-wall1 180 Coul-SR:wall1-wall1
>>>>>>>>>>
>>>>>>>>>> 181 LJ-SR:wall1-wall1 182 Coul-14:wall1-wall1
>>>>>>>>>>
>>>>>>>>>> 183 LJ-14:wall1-wall1 184 T-System
>>>>>>>>>>
>>>>>>>>>> 185 Lamb-System
>>>>>>>>>>
>>>>>>>>>> how can i understand that this 4nS simulation is enough for all 30
>>>>>>>>>> windows
>>>>>>>>>> which i selected after pulling? as you see i'm doing NVT.i tried to
>>>>>>>>>> calculate temperature(15) and total energy(13) and potential(11) for
>>>>>>>>>> example for first window;
>>>>>>>>>> Energy Average Err.Est. RMSD Tot-Drift
>>>>>>>>>> ------------------------------------------------------------
>>>>>>>>>> -------------------
>>>>>>>>>> Potential -72708 9.2 227.038 -38.8081
>>>>>>>>>> (kJ/mol)
>>>>>>>>>> Total Energy -63763.9 8.7 273.747 -34.9771
>>>>>>>>>> (kJ/mol)
>>>>>>>>>> Temperature 299.977 0.022 5.04346
>>>>>>>>>> 0.128492 (K)
>>>>>>>>>> plots show that temperature is ok and total energy doesn't have
>>>>>>>>>> considerable fluctuation (except after first moment). but potential
>>>>>>>>> None of those values tell you anything about the convergence of your
>>>>>>>>> simulations.
>>>>>>>>>
>>>>>>>>> fluctuates. since i don't see sharp peaks in PMF (but good trends),
>>>>>>>>>> is it
>>>>>>>>>> rational to do long time simulation? what is the main clue for
>>>>>>>>>> knowing
>>>>>>>>>> that? is it potential?
>>>>>>>>> If you want to know if the PMF has converged, calculate it for
>>>>>>>>> consecutive, non-overlapping time periods. If they are statistically
>>>>>>>>> indistinguishable, you are done. If they are still changing over time, you
>>>>>>>>> need longer runs.
>>>>>>>> You mean i should use the output Of 4nS simulation, as an input
>>>>>>>> for,for example next 10nS and at the end i have simulated it for 14nS,yes?
>>>>>>> I don't have any idea about what arbitrary amount of time you should
>>>>>>> target for your simulation, but you need to let convergence testing be your
>>>>>>> guide.
>>>>>>>
>>>>>>> What do you mean "statistically", would you please give me an example?
>>>>>>> Use the bootstrapping method of gmx wham to produce error estimates.
>>>>>>> Then plot the PMF curves and if they are overlapping
>>>>>>> <https://maps.google.com/?q=he+PMF+curves+and+if+they+are+overlapping&entry=gmail&source=g>
>>>>>>> within error, then you have an argument for convergence.
>>>>>>>
>>>>>>> -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.biochem.vt.edu/people/faculty/JustinLemkul.html
>>>>>>>
>>>>>>> ==================================================
>>>>>>>
>>>>>>> --
>>>>>>> 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.
>>>>>> --
>>>>> ==================================================
>>>>>
>>>>> 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.biochem.vt.edu/people/faculty/JustinLemkul.html
>>>>>
>>>>> ==================================================
>>>>>
>>>>> --
>>>>> 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.
>>> --
>>> ==================================================
>>>
>>> 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.biochem.vt.edu/people/faculty/JustinLemkul.html
>>>
>>> ==================================================
>>>
>>> --
>>> 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.
>
> --
> ==================================================
>
> 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.biochem.vt.edu/people/faculty/JustinLemkul.html
>
> ==================================================
>
> --
> 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