Numerical search of discontinuities in self-consistent potential energy surfaces

Abstract: Potential energy surfaces calculated with self-consistent mean-field methods are a very powerful tool, since their solutions are, in theory, global minima of the non-constrained subspace. However, this minimization leads to an incertitude concerning the saddle points, that can sometimes be no more saddle points in bigger constrained subspaces (fake saddle points), or can be missing on a trajectory (it missing saddle points). These phenomena are the consequences of discontinuities of the self-consistent potential energy surfaces (SPES). These discontinuities may have important consequences, since they can for example hide the real height of an energy barrier, and avoid any use of a SPES for further dynamical calculations, barrier penetrability estimations, or trajectory predictions. Discontinuities are not related to the quality of the production of a SPES, since even a perfectly converged SPES with an ideally fine mesh can be discontinuous. In this paper we explain what are the discontinuities, their consequences, and their origins. We then propose a numerical method to detect and identify discontinuities on a given SPES, and finally we discuss what are the best ways to transform a discontinuous SPES into a continuous one.