Populations of carbonic acid isomers at 210 K from a fast two-electron reduced-density matrix theory

Parametrization of the 2-electron reduced density matrix (2-RDM) rather than the many-electron wave function yields a new family of electronic-structure methods that are faster and more accurate than traditional coupled electron-pair methods including coupled cluster with single and double excitations. Deriving the parametrization from N-representability conditions generates a 2-RDM that captures significant correlation from triple and higher-order excitations at the cost of double excitations. We apply the parametric 2-RDM method to confirm recent experiments determining the relative thermodynamic populations of the cis-cis and cis-trans isomers of carbonic acid. In 2010 Bernard et al. showed by infrared spectroscopy that the populations of cis-cis and cis-trans isomers have a 10:1 ratio at 210 K. By use of the parametric 2-RDM method, we predict a 8:1 ratio at 210 K. Comparable ab initio methods overestimate the stability of the cis-cis isomer with 24:1 and 21:1 ratios. These 2-RDM-based methods promise to have significant applications throughout chemistry.