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.

Balancing error budget for fermionic k-RDM estimation

Error-mitigated fermionic classical shadows on noisy quantum devices

The Cost of Entanglement Renormalization on a Fault-Tolerant Quantum Computer

Geometric Constraints on Two-electron Reduced Density Matrices

Fermionic Partial Tomography via Classical Shadows

Efficient Strategies for Reducing Sampling Error in Quantum Krylov Subspace Diagonalization

Fermionic reduced density low-rank matrix completion, noise filtering, and measurement reduction in quantum simulations

Error-Robust Quantum Signal Processing using Rydberg Atoms

Quantum-classical hybrid algorithm using an error-mitigating $N$-representability condition to compute the Mott metal-insulator transition

Demonstrating Bayesian Quantum Phase Estimation with Quantum Error Detection

Classical shadows of fermions with particle number symmetry

Post-processing noisy quantum computations utilizing N-representability constraints

Error-mitigated Quantum Gates Exceeding Physical Fidelities in a Trapped-Ion System

Measurement-efficient quantum Krylov subspace diagonalisation

Sampling Error Analysis in Quantum Krylov Subspace Diagonalization

Qubit frugal entanglement determination with the deep multi-scale entanglement renormalization ansatz

An efficient adaptive variational quantum solver of the Schrödinger equation based on reduced density matrices

Robust design under uncertainty in quantum error mitigation

Error-mitigated quantum simulation of interacting fermions with trapped ions

Relativistic VQE calculations of molecular electric dipole moments on trapped ion quantum hardware

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