Search papers, labs, and topics across Lattice.
This paper introduces a novel damped alternating direction method of multipliers (DASSP) for constrained density functional theory (CDFT) that effectively addresses the challenges posed by Stiefel manifold and nonconvex quadratic constraints. By reformulating the optimization problem to decouple constraints and leveraging rotation invariance, DASSP achieves rigorous convergence guarantees and high feasibility accuracy. Numerical experiments show that DASSP not only improves efficiency but also maintains robustness, marking a significant advancement for large-scale CDFT applications.
The first algorithm for constrained density functional theory with rigorous convergence guarantees achieves high accuracy and efficiency without compromising robustness.
Constrained density functional theory (CDFT) provides a powerful framework for describing electronically excited and charge-localized states, which underlie a broad range of physical and chemical phenomena. However, the discretized optimization problems arising from CDFT calculations remain challenging, owing to the presence of both the Stiefel manifold constraint and additional nonconvex quadratic constraints. Existing algorithms either fail to enforce the quadratic constraints with high accuracy or face convergence issues due to double-loop iterative structures. In this paper, we first derive a subspace-splitting reformulation that decouples the two groups of constraints, by exploiting the inherent rotation invariance and introducing a nonlinear subspace alignment constraint. Based on this reformulation, we propose a single-loop damped alternating direction method of multipliers, called DASSP. To the best of our knowledge, DASSP is the first algorithm for CDFT calculations with rigorous convergence guarantees. Each iteration of DASSP comprises a spectral minimization step, a projected gradient step, and a damped dual ascent step, all of which admit efficient implementations. Numerical results on synthetic and realistic CDFT problems demonstrate that DASSP attains high feasibility accuracy and exhibits favorable efficiency without compromising robustness. We expect that this work will pave the way toward reliable and efficient large-scale CDFT applications.
