Variational calculation of second-order reduced density matrices by strong N-representability conditions and an accurate semidefinite programming solver
- 28 April 2008
- journal article
- research article
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 128 (16), 164113
- https://doi.org/10.1063/1.2911696
Abstract
The reduced density matrix (RDM) method, which is a variational calculation based on the second-order reduced density matrix, is applied to the ground stateenergies and the dipole moments for 57 different states of atoms, molecules, and to the ground stateenergies and the elements of 2-RDM for the Hubbard model. We explore the well-known N -representability conditions ( P , Q , and G ) together with the more recent and much stronger T 1 and T 2 ′ conditions. T 2 ′ condition was recently rederived and it implies T 2 condition. Using these N -representability conditions, we can usually calculate correlationenergies in percentage ranging from 100% to 101%, whose accuracy is similar to CCSD(T) and even better for high spin states or anion systems where CCSD(T) fails. Highly accurate calculations are carried out by handling equality constraints and/or developing multiple precision arithmetic in the semidefinite programming (SDP) solver. Results show that handling equality constraints correctly improves the accuracy from 0.1 to 0.6 mhartree . Additionally, improvements by replacing T 2 condition with T 2 ′ condition are typically of 0.1 – 0.5 mhartree . The newly developed multiple precision arithmetic version of SDP solver calculates extraordinary accurate energies for the one dimensional Hubbard model and Be atom. It gives at least 16 significant digits for energies, where double precision calculations gives only two to eight digits. It also provides physically meaningful results for the Hubbard model in the high correlation limit.Keywords
This publication has 68 references indexed in Scilit:
- Large-scale semidefinite programs in electronic structure calculationMathematical Programming, 2006
- Variational reduced-density-matrix method using three-particle-representability conditions with application to many-electron moleculesPhysical Review A, 2006
- Local Minima and Convergence in Low-Rank Semidefinite ProgrammingMathematical Programming, 2004
- Realization of Quantum Chemistry without Wave Functions through First-Order Semidefinite ProgrammingPhysical Review Letters, 2004
- The reduced density matrix method for electronic structure calculations and the role of three-index representability conditionsThe Journal of Chemical Physics, 2004
- Variational calculations of fermion second-order reduced density matrices by semidefinite programming algorithmThe Journal of Chemical Physics, 2001
- Contracted Schrödinger equation: Determining quantum energies and two-particle density matrices without wave functionsPhysical Review A, 1998
- Approximatingq-order reduced density matrices in terms of the lower-order ones. II. ApplicationsPhysical Review A, 1993
- Approximatingq-order reduced density matrices in terms of the lower-order ones. I. General relationsPhysical Review A, 1993
- Equation for the direct determination of the density matrixPhysical Review A, 1976