







C++ Codes to Solve Semidefinite Programs (SDPs)
SDPA 6.00 by (M. Yamashita, K. Fujisawa, and M. Kojima)
SDPARA 0.90 by (K. Fujisawa, M. Kojima, K. Nakata, and M. Yamashita)
SDP Benchmark Problems from Electronical Structure CalculationsThe following set of SDPs are the actual ones solved in [1]. The problem consists of determining the ground state energy of an N-electron system in a given external potential through the so-called Reduced Density Matrix Method. A detailed description on the formulation be found in: [1] Z. Zhao, B. J. Braams, M. Fukuda, M. L. Overton, and J. K. Percus, "The reduced density matrix method for electronic structure calculations and the role of three-index representability", October, 2003. The sizes of SDPs depend on the size of the "spin-orbital basis", r, which governs the discretization of the continuous wave function, but do not depend on the number of electrons N of the system. (Accurate electronic structure calculations will have r >= 2N). The Benchmark SDPs are in SDPA format and are divided according to the basis size r. See sizes of corresponding SDPs. In our paper [1], we present several SDP calculations for each molecule, but here we only provide the data for the most accurate SDP calculations: those with the P, Q, G, T1 and T2 N-representability conditions. Our results are given in Tables 1 and 2 in [1]. To obtain the energies displayed in Tables 1 and 2 under column $E_{PQGT1T2}$ [1], it is necessary to add the nuclear repulsion energies to the optimal values of each corresponding SDP. The important information we want to extract from solving the SDPs are the optimal value, and the dual vector variable (which provides the One-Reduced Density Matrix and the Two-Reduced Density Matrix). The applications typically require 4 digits of accuracy after the period for the optimal values (which means 6 or 7 digits of accuracy in total). Electronical Structure Calculation Problems (47 problems)
15 SDPs with r=12 (4.1MB)
6 SDPs with r=14 (3.6MB)
2 SDPs with r=16 (2.4MB)
24 SDPs with r=20 (86.2MB)
| r | # equality constraints | Size of block matrices | Size of the diagonal matrix |
|---|---|---|---|
| 12 | 948 | 6, 6, 6, 6, 15, 15, 36, 15, 15, 36, 72, 36, 36, 20, 90, 90, 20, 306, 306, 90, 90 | 94 |
| 14 | 1743 | 7, 7, 7, 7, 21, 21, 49, 21, 21, 49, 98, 49, 49, 35, 147, 147, 35, 490, 490, 147, 147 | 122 |
| 16 | 2964 | 8, 8, 8, 8, 28, 28, 64, 28, 28, 64, 128, 64, 64, 56, 224, 224, 56, 736, 736, 224, 224 | 154 |
| 20 | 7230 | 10, 10, 10, 10, 45, 45, 100, 45, 45, 100, 200, 100, 100, 120, 450, 450, 120, 1450, 1450, 450, 450 | 230 |
| System | N | r | Repulsion Energy |
|---|---|---|---|
| LiH | 4 | 12 | 0.9948810848 |
| BeH | 5 | 12 | 1.5772794307 |
| BH+ | 5 | 12 | 2.1782219849 |
| BH | 6 | 12 | 2.1469378814 |
| CH+ | 6 | 12 | 2.8075545972 |
| CH | 7 | 12 | 2.8351312564 |
| CH- | 8 | 12 | 2.9398736056 |
| NH+ | 7 | 12 | 3.4619072364 |
| NH | 8 | 12 | 3.5748318307 |
| NH- | 9 | 12 | 3.5379567746 |
| OH+ | 8 | 12 | 4.1145086908 |
| OH | 9 | 12 | 4.3656986614 |
| OH- | 10 | 12 | 4.3643484454 |
| HF+ | 9 | 12 | 4.7578374036 |
| HF | 10 | 12 | 5.1948028370 |
| BH2 | 7 | 14 | 4.7269615539 |
| CH2(1A1) | 8 | 14 | 6.0430274001 |
| CH2(3B1) | 8 | 14 | 6.1757874267 |
| NH2 | 9 | 14 | 7.5640948299 |
| H2O+ | 9 | 14 | 8.7976559250 |
| H2O | 10 | 14 | 9.1949655162 |
| NH3 | 10 | 16 | 11.9103794143 |
| H3O+ | 10 | 16 | 14.0016864732 |
| Li2 | 6 | 20 | 1.7818082386 |
| B2 | 10 | 20 | 8.3203969969 |
| C2+ | 11 | 20 | 13.5783185773 |
| C2 | 12 | 20 | 15.3322985626 |
| C2- | 13 | 20 | 15.0239597508 |
| LiF | 12 | 20 | 9.1359970094 |
| BeO | 12 | 20 | 12.7234743166 |
| NaH | 12 | 20 | 3.0841102782 |
| BeF | 13 | 20 | 13.9973408993 |
| BO | 13 | 20 | 17.5733416023 |
| N2+ | 13 | 20 | 23.2261601585 |
| N2 | 14 | 20 | 23.6218321955 |
| CO+ | 13 | 20 | 22.7786816895 |
| CO | 14 | 20 | 22.5121935230 |
| BF | 14 | 20 | 18.8617633307 |
| AlH | 14 | 20 | 4.1743350953 |
| CF | 15 | 20 | 22.4686046910 |
| O2+ | 15 | 20 | 30.3362091867 |
| O2 | 16 | 20 | 28.0474898021 |
| SiH | 15 | 20 | 4.8749631414 |
| SiH- | 16 | 20 | 5.0261068426 |
| NO- | 16 | 20 | 23.3154413407 |
| NF | 16 | 20 | 25.3156402817 |
| HS+ | 16 | 20 | 6.1603870664 |