Hello Charles, > what are the values you get? Da=-20 R=0.67 > I would suggest > using the anneal.py script as a template- this script only allows Da and > Rh to vary only during final minimzation. Yes, that what I am using. R does not change. It is 0.667. But Da=-50. Irina.