Core Concepts
Improving convergence properties of the single-reference coupled cluster method through an augmented Lagrangian formalism.
Abstract
The article introduces an augmented Lagrangian method (alm-CC) to enhance the stability and convergence of the coupled-cluster (CC) approach. It addresses issues with conventional CC methods, such as non-convergence and convergence to unphysical states. The proposed method minimizes an energy term alongside solving CC equations, transforming root-finding into a minimization problem. Numerical experiments demonstrate superior convergence performance with comparable computational costs to traditional CC methods. The content is structured around Introduction, Previous Works, Contribution, Optimizing Coupled Cluster, The Cluster Matrices, Numerics, Agreement with Quasi-Newton Method, Global Convergence Behavior sections.
Introduction:
- Describes the long-standing challenge of describing electrons in quantum many-body systems.
Previous Works:
- Discusses numerical methods in computational chemistry and physics community for tackling large problems efficiently.
Contribution:
- Proposes an augmented Lagrangian formulation (alm-CC) to address convergence issues in CC methods by minimizing energy terms alongside solving CC equations.
Optimizing Coupled Cluster:
- Introduces optimization framework as an alternative to conventional CC method using augmented Lagrangian approach.
The Cluster Matrices:
- Details the mathematical structures underlying the CC ansatz focusing on fermionic creation and annihilation matrices.
Numerics:
- Compares alm-CCSD with conventional CCSD for various scenarios like helium atom and dissociation processes showing improved convergence stability.
Agreement with Quasi-Newton Method:
- Demonstrates agreement between alm-CCSD and conventional CCSD for well-known scenarios where CCSD works reliably.
Global Convergence Behavior:
- Investigates multiple roots issue in H4 model system showcasing enhanced convergence stability of alm over Newton's method.
Stats
The number K is related to the number of basis functions used for discretizing the electronic Schrödinger equation.
For more complex systems we may easily obtain matrices of size 280 × 280 which would naively require about 10^22 GB to merely store the matrix.
Quotes
"An augmented Lagrangian formalism can help address these issues without a significant performance overhead."
"We propose to minimize an energy term in addition to solving the CC equations."