Variational calculation of second-order reduced density matrices by strong N -representability conditions and an accurate semidefinite programming solver

Maho Nakata, Bastiaan J. Braams, Katsuki Fujisawa, Mituhiro Fukuda, Jerome K. Percus, Makoto Yamashita, Zhengji Zhao

Research output: Contribution to journalArticle

51 Citations (Scopus)

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 state energies and the dipole moments for 57 different states of atoms, molecules, and to the ground state energies 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 T1 and T 2′ conditions. T 2′ condition was recently rederived and it implies T2 condition. Using these N -representability conditions, we can usually calculate correlation energies 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 T2 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.

Original languageEnglish
Article number164113
JournalJournal of Chemical Physics
Volume128
Issue number16
DOIs
Publication statusPublished - May 8 2008
Externally publishedYes

Fingerprint

Hubbard model
programming
Ground state
Atoms
Dipole moment
digits
Anions
energy
Molecules
ground state
matrix methods
atoms
dipole moments
anions
matrices

All Science Journal Classification (ASJC) codes

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Cite this

Variational calculation of second-order reduced density matrices by strong N -representability conditions and an accurate semidefinite programming solver. / Nakata, Maho; Braams, Bastiaan J.; Fujisawa, Katsuki; Fukuda, Mituhiro; Percus, Jerome K.; Yamashita, Makoto; Zhao, Zhengji.

In: Journal of Chemical Physics, Vol. 128, No. 16, 164113, 08.05.2008.

Research output: Contribution to journalArticle

Nakata, Maho ; Braams, Bastiaan J. ; Fujisawa, Katsuki ; Fukuda, Mituhiro ; Percus, Jerome K. ; Yamashita, Makoto ; Zhao, Zhengji. / Variational calculation of second-order reduced density matrices by strong N -representability conditions and an accurate semidefinite programming solver. In: Journal of Chemical Physics. 2008 ; Vol. 128, No. 16.
@article{2d6187d2534845e7b4f21b1a05601b56,
title = "Variational calculation of second-order reduced density matrices by strong N -representability conditions and an accurate semidefinite programming solver",
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 state energies and the dipole moments for 57 different states of atoms, molecules, and to the ground state energies 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 T1 and T 2′ conditions. T 2′ condition was recently rederived and it implies T2 condition. Using these N -representability conditions, we can usually calculate correlation energies 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 T2 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.",
author = "Maho Nakata and Braams, {Bastiaan J.} and Katsuki Fujisawa and Mituhiro Fukuda and Percus, {Jerome K.} and Makoto Yamashita and Zhengji Zhao",
year = "2008",
month = "5",
day = "8",
doi = "10.1063/1.2911696",
language = "English",
volume = "128",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "16",

}

TY - JOUR

T1 - Variational calculation of second-order reduced density matrices by strong N -representability conditions and an accurate semidefinite programming solver

AU - Nakata, Maho

AU - Braams, Bastiaan J.

AU - Fujisawa, Katsuki

AU - Fukuda, Mituhiro

AU - Percus, Jerome K.

AU - Yamashita, Makoto

AU - Zhao, Zhengji

PY - 2008/5/8

Y1 - 2008/5/8

N2 - The reduced density matrix (RDM) method, which is a variational calculation based on the second-order reduced density matrix, is applied to the ground state energies and the dipole moments for 57 different states of atoms, molecules, and to the ground state energies 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 T1 and T 2′ conditions. T 2′ condition was recently rederived and it implies T2 condition. Using these N -representability conditions, we can usually calculate correlation energies 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 T2 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.

AB - The reduced density matrix (RDM) method, which is a variational calculation based on the second-order reduced density matrix, is applied to the ground state energies and the dipole moments for 57 different states of atoms, molecules, and to the ground state energies 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 T1 and T 2′ conditions. T 2′ condition was recently rederived and it implies T2 condition. Using these N -representability conditions, we can usually calculate correlation energies 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 T2 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.

UR - http://www.scopus.com/inward/record.url?scp=42949093565&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=42949093565&partnerID=8YFLogxK

U2 - 10.1063/1.2911696

DO - 10.1063/1.2911696

M3 - Article

C2 - 18447427

AN - SCOPUS:42949093565

VL - 128

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 16

M1 - 164113

ER -