[gmx-users] Why are the bond and angle energy at the first step of energy minimization not equal 0?
trinhlanhoa at gmail.com
Sat Mar 17 20:09:09 CET 2018
Great explanation! Thank you very much.
2018-03-17 13:02 GMT-05:00 Justin Lemkul <jalemkul at vt.edu>:
> On 3/17/18 1:56 PM, Hoa Trinh wrote:
>> Hi Gromacs user,
>> I am following a Gromacs tutorial on the website of Dr. Justin:
>> After performing energy minimization, I open the log file and see that at
>> time step = 0 (time = 0.000), bond energy and angle energy are not equal
>> Given that: E_bond = 0.5*k*(r-r0)^2
>> At t = 0: r = r0, therefore E_bond should be approximate 0.
>> Please correct me.
> Simply, r =/= r0 at t = 0 (and in energy minimization it's not time, it's
> step number).
> Unless a structure was built *exactly* according to the force field
> parameters in terms of equilibrium bond lengths and valence angles, you
> will never have a zero value for these terms. Experimental structures
> certainly never conform exactly to these parameters, nor will they
> throughout the simulation.
> 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
> 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