Next: Preliminary results on high
Up: Practical implementation of the
Previous: Following the ionic dynamics
Contents
Reducing the number of parameters
After each ionic move it is important to optimize the trial function to be as close as possible to the variational ground state. The point is that there are some parts of the wave function that vary a bit between different ionic configurations, as for example the two-body Jastrow factor or the core orbitals for large ions. So it is not important to move all parameters at each step, but it is fundamental to recognize which ones have to be optimized because they are far from the minimum. Moreover since we are performing an optimization in the presence of noise it is possible to know the exact minimum only within a given statistical error due to the finite sampling.
So we fix a tolerance
and move only such parameters whose generalized forces satisfy:
|
(5.29) |
we have chosen
that amounts to change on average the
of the variational parameters. This procedure has allowed us to follow the ion dynamic with a stable and fast optimization on the variational parameters.
Moreover to stabilize the optimization procedure other two cautions were used:
1. It is very difficult to parametrize an highly correlated wave-function with many parameters. In fact it happens often that in this case the stochastic matrix becomes singular
because there are too many parameters redundant to describe the wave-function and some of them have to be eliminated for a stable optimization scheme, see section 2.1.2.
2. It can happen that the Hessian Matrix is not positive definite, this can be due to different reasons: a statistical fluctuation due the finite sampling or because the trial-function is far from the minimum, in this case the correction, proposed in Ref. (50) and described at the end of section 2.3, is used. In our simulations we have used
as a threshold value for the
parameters that rules the stability of the SRH optimization (see section 2.3). This has achieved a fast and stable convergence of the wave-function.
Next: Preliminary results on high
Up: Practical implementation of the
Previous: Following the ionic dynamics
Contents
Claudio Attaccalite
2005-11-07