Impact of Parametrizations of the One-Body Reduced Density Matrix on the Energy Landscape

Abstract

Many electronic structure methods rely on the minimization of the energy of the system with respect to the one-body reduced density matrix (1-RDM). To formulate a minimization algorithm, the 1-RDM is often expressed in terms of its eigenvectors via an orthonormal transformation and its eigenvalues. This transformation drastically alters the energy landscape. Especially in 1-RDM functional theory this means that the convexity of the energy functional is lost. We show that degeneracies in the occupation numbers can lead to additional critical points which are classified as saddle points. Using a Cayley or Householder parametrization for the orthonormal transformation, no extra critical points arise. In case of Given's rotations or the exponential, additional critical points can arise, which are of no concern in practical minimization. These findings provide an explanation for the success of recent minimization procedures using second-order information.

17 pages, 2 figures, html version not recommended

References in corpus (12)

• The Pauli principle revisited
• Global Method for Electron Correlation
• Interpretations of ground-state symmetry breaking and strong correlation in wavefunction and density functional theories
• Global Natural Orbital Functional: Towards the Complete Description of the Electron Correlation
• Kohn-Sham calculations with the exact functional
• Relating the pure and ensemble density matrix functional
• Generalized Pauli constraints in reduced density matrix functional theory
• Open shells in reduced-density-matrix-functional theory
• DoNOF: an open-source implementation of natural-orbital-functional-based methods for quantum chemistry
• Exploiting the hessian for a better convergence of the SCF RDMFT procedure
• Co-iterative augmented Hessian method for orbital optimization
• Sylvester's Minorant Criterion, Lagrange-Beltrami Identity, and Nonnegative Definiteness