From: Chris Harrison (charris5_at_gmail.com)
Date: Mon Aug 17 2009 - 10:43:24 CDT
Yoaz,
As Eddi said, this is indicative of high nuclear clashes. A suggestion -->
check to see if any atoms are overlapping (literally: have the same
coordinates) in the initial pdb. For various reasons involving topology
definitions and pdb names, this can occur and often one finds these
overlapped atoms at coordinates: 0.0, 0.0, 0.0, so this may be something for
which to specifically watch.
Chris
-- Chris Harrison, Ph.D. Theoretical and Computational Biophysics Group NIH Resource for Macromolecular Modeling and Bioinformatics Beckman Institute for Advanced Science and Technology University of Illinois, 405 N. Mathews Ave., Urbana, IL 61801 char_at_ks.uiuc.edu Voice: 217-244-1733 http://www.ks.uiuc.edu/~char Fax: 217-244-6078 On Mon, Aug 17, 2009 at 10:33 AM, Schreiner Eduard <eschrein_at_ks.uiuc.edu>wrote: > your VDW energy if very high. so, apparently, > you have clashes in your system (atoms too close to each other). > if it does go like that endlessly, it seems to me that these clashes > cannot be resolved in the minimization. > btw, is there a special reason to minimize at constant pressure? > I usually do minimizations at constant volume. > > eddi > > > > On Aug 17, 2009, at 10:20 AM, Yoav Beck wrote: > > Dear all, >> >> Thanks for the help on my previous questions. I am trying to minimize a >> pope membrane (no protein yet, but lipid tails had already been melted) with >> water on top and bottom. I have done it previously without problem. This >> time when I start the minimization, I get this odd result : >> >> TCL: Minimizing for 2500 steps >> PRESSURE: 0 nan nan nan nan nan nan nan nan nan >> GPRESSURE: 0 nan nan nan nan nan nan nan nan nan >> ETITLE: TS BOND ANGLE DIHED >> IMPRP ELECT VDW BOUNDARY MISC >> KINETIC TOTAL TEMP POTENTIAL TOTAL3 >> TEMPAVG PRESSURE GPRESSURE VOLUME >> PRESSAVG GPRESSAVG >> >> ENERGY: 0 18740.4117 2370123.0355 44852.5867 >> 202077.6674 -389262.6409 9999999999.9999 0.0000 >> 0.0000 0.0000 9999999999.9999 0.0000 9999999999.9999 >> 9999999999.9999 0.0000 9999999999.9999 9999999999.9999 >> 1403600.0000 9999999999.9999 9999999999.9999 >> >> OPENING EXTENDED SYSTEM TRAJECTORY FILE >> MINIMIZER SLOWLY MOVING ATOMS WITH BAD CONTACTS DOWNHILL >> PRESSURE: 1 nan nan nan nan nan nan nan nan nan >> GPRESSURE: 1 nan nan nan nan nan nan nan nan nan >> ENERGY: 1 19929.3605 2370494.7192 44860.0041 >> 202085.4401 -389159.7118 9999999999.9999 0.0000 >> 0.0000 0.0000 9999999999.9999 0.0000 9999999999.9999 >> 9999999999.9999 0.0000 9999999999.9999 9999999999.9999 >> 1403600.0000 9999999999.9999 9999999999.9999 >> .. >> LINE MINIMIZER REDUCING GRADIENT FROM nan TO nan >> LINE MINIMIZER BRACKET: DX 0 0 DU 0 0 DUDX nan nan nan >> MINIMIZER RESTARTING CONJUGATE GRADIENT ALGORITHM >> LINE MINIMIZER REDUCING GRADIENT FROM nan TO nan >> LINE MINIMIZER BRACKET: DX 0 0 DU 0 0 DUDX nan nan nan >> MINIMIZER RESTARTING CONJUGATE GRADIENT ALGORITHM >> LINE MINIMIZER REDUCING GRADIENT FROM nan TO nan >> and it goes on endlessly. >> >> Any help will be very appreciate. >> >> Cheers, >> >> Yoav. >> >> > > ===================================================== > Eduard Schreiner > Theoretical and Computational Biophysics Group > NIH Resource for Macromolecular Modeling and Bioinformatics > Beckman Institute for Advanced Science and Technology > University of Illinois, 405 N. Mathews Ave., Urbana, IL 61801 > Phone: 217-244-4361 > Fax: 217-244-6078 > http://www.ks.uiuc.edu/~eschrein/ <http://www.ks.uiuc.edu/%7Eeschrein/> > ============================================= > >
This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:53:11 CST