Re: AMBER: SCF convergence issues

From: Steven Winfield <>
Date: Fri, 02 Mar 2007 13:43:58 +0000

Dear Ross

Thanks for your previous help.

>> I built the system from a few TIP3PBOXs, manually adjusted the box
>> dimensions to give the correct density at 300K, performed a
>> few thousand
>> minimisation steps (about half steepest descent), then
>> performed MD from
>> 0K -> 300K over 100ps, followed by another 100ps at 300K, both NVT.
> Was this with QM/MM MD? If so then is definately a problem with
> incompatibility ofr dumpfrc and QMMM.

No, this was done classically. I've bitten the bullet and am now
equilibrating a larger (45A cube) box of water, again classically:

10000 steps minimisation (2000 steepest descent)
200ps NVT 0K -> 300K
1ns NPT 300K

Should I do some minimisation after this with QM/MM turned on before I
do QM/MM MD?

> The safer way would just be to put a write(20,*) f(1:3*natom) at the end of
> the force routine followed by a call to mexit but this obviously involved
> adding a couple of lines to force.f. I don't know if you would be okay to do
> this or not.

Thanks for the suggestion.


The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Sun Mar 04 2007 - 06:07:53 PST
Custom Search