Here is how the geometry changes during the optimization
procedure:
|
O-H bond
length (Å) |
H-O-H angle
(º) |
Initial model |
0.9600 |
109.471 |
Step 1 |
0.9888 |
95.817 |
Step 2 |
0.9932 |
100.450 |
Step 3 |
0.9893 |
100.035 |
How does the program know when to stop?
i.e. how does it know when it has
reached a minimum?
Convergence
criteria are used to test if the minimization has converged
For example:
· monitoring change
in energy from one step to the next, stop the minimization when the change
falls below a threshold
· monitor change in
coordinates, stop when structure is changing only a little
· stop when the
root mean square of the gradient is small,
, i.e. the square
root of the gradient (with respect to the coordinates) squared divided by the
number of coordinates.
· The biggest
element of g is also monitored