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.

arXiv comment17 pages, 2 figures, html version not recommended

References in corpus (19)

• Perspective on density functional theory
• Recent developments in the PySCF program package
• A natural orbital functional for the many-electron problem
• The Pauli principle revisited
• Global Method for Electron Correlation
• Interpretations of ground-state symmetry breaking and strong correlation in wavefunction and density functional theories
• Benchmark calculations for reduced density-matrix functional theory
• An overview of self-consistent field calculations within finite basis sets
• 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
• Performance of one-body reduced density matrix functionals for the homogeneous electron gas
• 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
• Reduced density-matrix functionals applied to the Hubbard dimer
• 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