mersenneforum.org > Math Energy Minimization
 Register FAQ Search Today's Posts Mark Forums Read

 2008-04-12, 12:54 #1 ShiningArcanine     Dec 2005 9210 Posts Energy Minimization I am doing undergraduate research in a lab that does molecular dynamics simulations and it struck me yesterday that the minimization of energy functions could be done with the methods I learned in vector calculus if one was to apply constraints to a system of equations with 3N variables, where N is the number of atoms.I went through a book on molecular modeling yesterday to see if it said anything regarding this and it stated that it is not generally possible to do that "for molecular systems due to the complicated way in which the energy varies with the coordinates." It then states that minimization is done with numerical methods, but it does not elaborate on why the way that the energy is modeled prevents use of the extreme value theorem. Would someone elaborate on the nature of problems where the extreme value theorem cannot be used for minimization of a system of variables? If I do not get a response here, I will ask in the lab on Monday, but people here seem to be really good at applied mathematics, so I was hoping that someone here could elaborate on what the book said.
2008-04-12, 13:14   #2
jasonp
Tribal Bullet

Oct 2004

DD916 Posts

Quote:
 Originally Posted by ShiningArcanine Would someone elaborate on the nature of problems where the extreme value theorem cannot be used for minimization of a system of variables?
My understanding is that the complexity of the minimization job depends on how accurately you describe the forces involved, and whether subatomic forces are included in that.

Your book says that because the common wisdom holds that using gradient-or Hessian-based minimization is impractical when there are tens of thousands of variables that all contribute to the minimum. Just calculating the forces that every atom exerts on all the others is a big job, and there is a large research literature on doing the computation efficiently. Finally, if you're looking for a global minimum somewhere inside a space that is very full of local minima, then you may never find it without massively sampling the entire space for the correct set of initial conditions needed. Certainly nobody expects a non-numercial solution to such large problems, at least without simplifying away most of the detail that makes the problem relevant.

Last fiddled with by jasonp on 2008-04-12 at 13:15

 2008-04-16, 13:47 #3 ShiningArcanine     Dec 2005 22·23 Posts Thank you so much for the information. I did not have much of a chance to go to my lab on Monday. Yesterday, I was there to get a book on Fortran. Today, I will be going there for a lab meeting, so I will ask then. Last fiddled with by ShiningArcanine on 2008-04-16 at 13:47

 Similar Threads Thread Thread Starter Forum Replies Last Post em99010pepe Lounge 31 2011-02-14 01:57 science_man_88 Miscellaneous Math 8 2010-05-29 04:14 biwema Hardware 0 2008-08-24 00:16 mfgoode Science & Technology 3 2006-11-29 07:46 markhl Hardware 5 2004-02-04 13:33

All times are UTC. The time now is 17:26.

Sun May 22 17:26:57 UTC 2022 up 38 days, 15:28, 1 user, load averages: 2.14, 1.72, 1.42