Hi,everyone,

I have a question about conjugate gradient minimization procedure: the forces on

atoms, espicially free surface atoms, are not close to zero. For example, the

forces on interior atoms are 10^(-4)ev/A(metal units), but the force normal to

the free surface can reach 0.1ev/A. The maximum iterations and force

evaluations are set to be very large and are not excedeed in simulations. Since

the force is closely associated with the gradient of the energy, do you think

the conjugate gradient method works well for my case?

Besides, there are other two criteria to stop the minimization

1. the change in energy between outer iterations is less than the tolerance

2. the 3N dimensional force vector goes to zero

I decrease the tolerance but cannot improve the result.

What kind of modification I should make if I want the second criterion (the

force vector goes to zero) to be the only one can stop the minimization?

Thanks,