Electronic structure calculations, in particular the computation of the ground state energy, lead to challenging problems in optimization. These problems are of enormous importance in quantum chemistry for calculations of properties of solids and molecules. Minimization methods for computing the ground state energy can be developed by employing a variational approach, where the second-order reduced density matrix defines the variable. This concept leads to large-scale semidefinite programming problems that provide a lower bound for the ground state energy. Upper bounds of the ground state energy can be calculated for example with the Hartree–Fock method or numerically more exact for a given basis set by full CI. However, Nakata et al. ( J. Chem. Phys.200111482828292) observed that due to numerical errors the semidefinite solver produced erroneous results with a lower bound significantly larger than the full CI energy. For the LiH, CH–, NH–, OH, OH–, and HF molecules violations within one mhartree were observed. We applied the software VSDP which takes all numerical errors due to floating-point arithmetic operations into consideration. For two test libraries VSDP provides tight rigorous error bounds lower than full CI energies reported with an accuracy of 0.1 to 0.01 mhartree. Only little computation work must be spent in order to compute close rigorous error bounds for the ground state energy.