AMBER Archive (2006)

Subject: RE: AMBER: rgbmax--something wrong with my input file?

From: Ross Walker (
Date: Tue Jun 27 2006 - 17:52:06 CDT

Hi Kenley,
You have something very wrong with your starting structure:
      1 6.8283E+12 4.3714E+11 3.8864E+13 CD 6756
 BOND = 55589.4701 ANGLE = 12184.7419 DIHED =
 VDWAALS = ************* EEL = -179753.6450 EGB =
 1-4 VDW = 2621189.7868 1-4 EEL = 113699.2076 RESTRAINT =

Note the *****'s in the VDWAALS term, this means the number is too large to
fit in the output space, plus your BOND energy is huge as is your 1-4 VDW
energy. You will need to work out what is wrong with this and fix this. You
probably have 1 or more atoms sitting on top of each other. I would start by
looking at atom 6756 (CD) as this is reported to have the biggest force on
it which is 3.886 * 10^13!!!! KCal/mol/A. This is HUGE... Anything else past
this point in the output file is pretty much meaningless as there is no way
the minimiser can cope with such a strained system.
As for the second problem with RGBMAX I cannot see anything obviously wrong
with your input file but you could try specifying the values that are floats
as floats instead of integers. I.e. rgbmax=30.0 rather than just 30. Same
with cut, and restraint_wt. You might also have some strange line endingings
in your file. Try creating a brand new file and typing in the options again.
On another note a restraint_wt of 500.0 is probably way too much. You can
get away with this for minimisation but it would be way too high for MD. Try
a value of 10.0 or so, you may get better results. You still need to fix
your input file first though.
All the best

|\oss Walker

| HPC Consultant and Staff Scientist |
| San Diego Supercomputer Center |
| Tel: +1 858 822 0854 | EMail:- |
| <> | PGP Key
available on request |

Note: Electronic Mail is not secure, has no guarantee of delivery, may not
be read every day, and should not be used for urgent or sensitive issues.



From: [] On Behalf Of
Kenley Barrett
Sent: Tuesday, June 27, 2006 14:10
Subject: AMBER: rgbmax--something wrong with my input file?

Dear AMBER community,

I am having problems running a simple implicit solvent minimization. I have
tried to run the job on two different computers. On the first computer, the
job completed but produced a very strange .out file in which, after a point,
the energy began rotating between -1, -9, and -3. I have attached this .out
file. On the second computer the job crashed with the error: " The NAMELIST
READ statement cannot be completed because item rgbmax is not a member of
the NAMELIST group cntrl. The program will recover by
discontinuing further processing of the READ statement..."

I have examined my .in file and I can't figure out how I used "rgbmax"
incorrectly, or what else I might have done wrong. If anyone could offer any
advice on what mistake I have made, I would be very grateful. My .in file is
pasted below.

Thank you in advance.


  restraintmask=:* & !@H= & !:WAT

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to