Abstract:Reduced density matrix functional theory (RDMFT) calculations are usually implemented in a decoupled manner, where the orbital and occupation optimizations are repeated alternately. Typically, orbital updates are performed using the unitary optimization method, while occupations are optimized through the explicit-by-implicit (EBI) method. The EBI method addresses explicit constraints by incorporating implicit functions, effectively transforming constrained optimization scenarios into unconstrained minimizations. Although the unitary and EBI methods individually achieve robust performance in optimizing orbitals and occupations, respectively, the decoupled optimization methods often suffer from slow convergence and require dozens of alternations between the orbital and occupation optimizations. To address this issue, this work proposes a coupled optimization method that combines unitary and EBI optimizations to update orbitals and occupations simultaneously at each step. To achieve favorable convergence in coupled optimization using a simple first-order algorithm, an effective and efficient preconditioner and line search are further introduced. The superiority of the new method is demonstrated through numerous tests on different molecules, random initial guesses, different basis sets and different functionals. It outperforms all decoupled optimization methods in terms of convergence speed, convergence results and convergence stability. Even a large system like $\mathrm{C_{60}}$ can converge to $10^{-8}$ au in 154 iterations, which shows that the coupled optimization method can make RDMFT more practical and facilitate its wider application and further development.

Density Matrix Minimization Meets a Combinatorial Optimization Problem

Optimal Reduction of Multivariate Dirac Mixture Densities

Method for Triangular Fuzzy Multiple Attribute Decision Making Based on Two-Dimensional Density Operator Method

The appeals of quadratic majorization–minimization

Global optimization for a class of fractional programming problems

Binary Optimization via Mathematical Programming with Equilibrium Constraints

A Minimization Approach for Minimax Optimization with Coupled Constraints

Enhancing Reduced Density Matrix Functional Theory Calculations by Coupling Orbital and Occupation Optimizations

Convex Matroid Optimization

Combinatorial Bayesian Optimization with Random Mapping Functions to Convex Polytopes

A Hybrid Method of Combinatorial Search and Coordinate Descent for Discrete Optimization

Quadratic Decomposable Submodular Function Minimization: Theory and Practice (Computation and Analysis of PageRank over Hypergraphs)

Diagonalization without Diagonalization: A Direct Optimization Approach for Solid-State Density Functional Theory

A hybrid discrete state transition algorithm for combinatorial optimization problems

High-dimensional Optimal Density Control with Wasserstein Metric Matching

Distance Majorization and Its Applications

A Generalized DCT Compression Based Density Method for Topology Optimization of 2D and 3D Continua

Minimization of Quadratic Binary Functional with Additive Connection Matrix

D-optimal designs via a cocktail algorithm

Binary matrix factorization via collaborative neurodynamic optimization

A Simple Effective Heuristic for Embedded Mixed-Integer Quadratic Programming