Today I have a question about my solvation. After immersing in a droplet of water, I try to equilbrium, but there is a abnormal termination when heating. This is the code:
DYNA DYN: Step Time TOTEner TOTKe ENERgy TEMPerature DYNA PROP: GRMS HFCTote HFCKe EHFCor VIRKe DYNA INTERN: BONDs ANGLes UREY-b DIHEdrals IMPRopers DYNA EXTERN: VDWaals ELEC HBONds ASP USER DYNA MISCEL: SBOUnd TSM PRMS PANGle other DYNA PRESS: VIRE VIRI PRESSE PRESSI VOLUme ---------- --------- --------- --------- --------- --------- DYNA> 0 0.00000-299487.95558 0.00000-299487.95558 0.00000 DYNA PROP> 28.72222-299327.78729 480.50490 160.16830 -13858.10193 DYNA INTERN> 0.00000 922.19829 207.86836 1105.91786 46.22287 DYNA EXTERN> 10501.32156-312122.45716 0.00000 0.00000 0.00000 DYNA MISCEL> -149.02736 0.00000 0.00000 0.00000 0.00000 DYNA PRESS> 210133.59958-200894.86496 0.00000 0.00000 0.00000 ---------- --------- --------- --------- --------- --------- UPDECI: Nonbond update at step 1 UPDECI: Nonbond update at step 6 UPDECI: Nonbond update at step 10 UPDECI: Nonbond update at step 13 UPDECI: Nonbond update at step 14 WRIDYN: RESTart file was written at step 15
Writing RESTART FILE with previous and current coordinates, which may be read by: READ COOR DYNR .... (see io.doc). NOTE!! THIS FILE C A N N O T BE USED TO RESTART A RUN!!!
***** LEVEL -2 WARNING FROM ***** ***** ENERGY CHANGE TOLERANCE EXCEEDED ****************************************** BOMLEV ( -1) IS REACHED - TERMINATING. WRNLEV IS 0
ABNORMAL TERMINATION MOST SEVERE WARNING WAS AT LEVEL -2
$$$$$ JOB ACCOUNTING INFORMATION $$$$$ ELAPSED TIME: 20.79 MINUTES CPU TIME: 20.70 MINUTES
So how can I deal with it? help would be greatly appreciated wenhao
One possibility is that the structure has not been adequately minimized after solvation. There could also be some issues with your dynamics input script. One thing to note is that the last error or warning message may be a symptom of an earlier problem.
Thank you so much. Now I still try my best to deal with this problem. But I dont know what is your meaning about "an earlier problem"? To be honest,after my first minimization, the first residue and the final resdue were abnormal(bonding of unrelated atoms that needn‘t bond). But I used an order:define heavy sele all .and. .not. ( bynu 1:28 .or. bynu 12547:12569 .or. resname TIP3 .or. hydrogen) end.Then the minimization was fine and the mini file were used as input file for solvation.
Thank you very much ,really. I should tell you about the step of minimization : the out PSF file of minimization before solvation is normal(displaying no bond), but when I open the out PDB-minimization using VMD,the image is abnormal for the first residue and last residue. So my teacher told me to use the order. I will cannot deal with the abnormal bonds without the order, sorry, I am a new learner for Charmm. I really thank you for your reply and your time.There are log files in the attachment. Thank you again and wish you happy every day.
VMD will not display bonds correctly unless you read a PSF before coordinates; for CHARMM, using the COOR file (VMD .cor) is recommended, using PDB is not
Input files are also quite informative, and a lot easier to read than lengthy log files. You should strive to post a very brief input file (and the corresponding log file) that shows the alleged bad behavior. Do not use wrnlev, prnlev or bomlev commands. Your logfiles are so cluttered with mostly uninteresting information that they are unreadable. A "define" command does not in itself do anything. The PSF is not changed by coordinate manipulations such as minimization. And please always use the most recent version of CHARMM.
Lennart Nilsson Karolinska Institutet Stockholm, Sweden
Yes，teacher. I definitely used the COOR file for Charmm. I really cannot deal with this abnormal termination, even if I change the inp file for equilbrium:shake bonh tolerance 1.0e-6 parameter (change the number). How can I do next?
Thank you for your reply and suggestions. I still cannot deal with this complicated termination in heating 300K. There are input files and corresponding log files in the attachment. Thank you again and wish you happy every day.
Please stop using PRNLEV 10 (or 15) for these output logs if you want us to look through them; use the default (PRNLEV 5), or PRNLEV 6 at most.
The CHARMM version is several years out of date.
You also appear to be using out of date, modified parameters, which could be one of the many problems in this set of files. The extraneous character warning indicates an issue with how the file edited at some point. You should resolve the errors that occur during parameter reading.
Code
CHARMM> ! Parameter files
CHARMM> open read unit 2 card name /home/dengwenhao/5ymr/par_all_isethionate.prm
VOPEN> Attempting to open::/home/dengwenhao/5ymr/par_all_isethionate.prm::
OPNLGU> Unit 2 opened for READONLY access to /home/dengwenhao/5ymr/par_all_isethionate.prm
CHARMM> read para card unit 2
PARAMETER FILE BEING READ FROM UNIT 2
TITLE> *>>>>>> COMBINED CHARMM ALL-HYDROGEN PARAMETER FILE FOR <<<<<<<<<
TITLE> *>>>>>>> CHARMM22 PROTEINS AND CHARMM27 NUCLEIC ACIDS <<<<<<<<<<
TITLE> *FROM
TITLE> *>>>> CHARMM22 ALL-HYDROGEN PARAMETER FILE FOR PROTEINS <<<<<<<<
TITLE> *>>>>>>>>>>>>>>>>>>>>>>> AUGUST 1999 <<<<<<<<<<<<<<<<<<<<<<<<<<<
TITLE> *>>>>>>> DIRECT COMMENTS TO ALEXANDER D. MACKERELL JR. <<<<<<<<<
TITLE> *>>>>>> 410-706-7442 OR EMAIL: ALEX,MMIRIS.AB.UMD.EDU <<<<<<<<<
TITLE> *AND
TITLE> * \\\\ CHARMM27 ALL-HYDROGEN NUCLEIC ACID PARAMETER FILE ////
TITLE> * \\\\\\\\\\\\\\\\\\ DEVELOPMENTAL ////////////////////////
TITLE> * ALEXANDER D. MACKERELL JR. AND NICOLAS FOLOPPE
TITLE> * AUGUST 1999
TITLE> * ALL COMMENTS TO ADM JR. EMAIL: ALEX,MMIRIS.AB.UMD.EDU
TITLE> * TELEPHONE: 410-706-7442
TITLE> *
**** Warning **** The following extraneous characters
were found while command processing in PARRDR
£°
PARRDR> WARNING: ATOM FOR NBOND OSL DOESNT EXIST
PARRDR> ALL ANGLES HAVE POSITIVE MINIMA
PARRDR> Error: Repeated BOND parameter ( 19992): SF FE
PARRDR> Error: Repeated BOND parameter ( 14901): ON2 CT2
PARRDR> Error: Repeated BOND parameter ( 11579): NN2 HN2
PARRDR> Error: Repeated BOND parameter ( 4269): FE SS
PARRDR> Error: Repeated BOND parameter ( 4246): FE NPH
PARRDR> Error: Repeated BOND parameter ( 4238): FE NR2
PARRDR> Error: Repeated BOND parameter ( 4223): FE CM
PARRDR> Error: Repeated BOND parameter ( 4222): FE CPM
The initial minimization energies indicate possible issues with the model building; very high VDW energy. You should also be using SHAKE BONH PARAM at this point, and with the default 1.e-9 tolerance; the same applies to the post-solvation minimization.
Why are you using CONS FIX for all the heavy atoms? This is highly inadvisable, and means you aren't really doing much minimization.
This is a good example of the wrong way to condition an initial model to prepare it for dynamics.
