ABSTRACT
Cavity quantum electrodynamics (QED) generalizations of time-dependent (TD) density functional theory (DFT) and equation-of-motion (EOM) coupled-cluster (CC) theory are used to model small molecules strongly coupled to optical cavity modes. We consider two types of calculations. In the first approach (termed "relaxed"), we use a coherent-state-transformed Hamiltonian within the ground- and excited-state portions of the calculations, and cavity-induced orbital relaxation effects are included at the mean-field level. This procedure guarantees that the energy is origin-invariant in post-self-consistent-field calculations. In the second approach (termed "unrelaxed"), we ignore the coherent-state transformation and the associated orbital relaxation effects. In this case, ground-state unrelaxed QED-CC calculations pick up a modest origin dependence but otherwise reproduce relaxed QED-CC results within the coherent-state basis. On the other hand, a severe origin dependence manifests in ground-state unrelaxed QED mean-field energies. For excitation energies computed at experimentally realizable coupling strengths, relaxed and unrelaxed QED-EOM-CC results are similar, while significant differences emerge for unrelaxed and relaxed QED-TDDFT. First, QED-EOM-CC and relaxed QED-TDDFT both predict that electronic states that are not resonant with the cavity mode are nonetheless perturbed by the cavity. Unrelaxed QED-TDDFT, on the other hand, fails to capture this effect. Second, in the limit of large coupling strengths, relaxed QED-TDDFT tends to overestimate Rabi splittings, while unrelaxed QED-TDDFT underestimates them, given splittings from relaxed QED-EOM-CC as a reference, and relaxed QED-TDDFT generally does the better job of reproducing the QED-EOM-CC results.
ABSTRACT
The electron attachment variant of equation-of-motion coupled-cluster theory (EOM-EA-CC) is generalized to the case of strong light-matter coupling within the framework of cavity quantum electrodynamics (QED). The resulting EOM-EA-QED-CC formalism provides an ab initio, correlated, and non-perturbative description of cavity-induced effects in many-electron systems that complements other recently proposed cavity-QED-based extensions of CC theory. Importantly, this work demonstrates that QED generalizations of EOM-CC theory are useful frameworks for exploring particle-non-conserving sectors of Fock space, thereby establishing a path forward for the simultaneous description of both strong electron-electron and electron-photon correlation effects.
ABSTRACT
The direct variational optimization of the two-electron reduced density matrix (2RDM) can provide a reference-independent description of the electronic structure of many-electron systems that naturally capture strong or nondynamic correlation effects. Such variational 2RDM approaches can often provide a highly accurate description of strong electron correlation, provided that the 2RDMs satisfy at least partial three-particle N-representability conditions (e.g., the T2 condition). However, recent benchmark calculations on hydrogen clusters [N. H. Stair and F. A. Evangelista, J. Chem. Phys. 153, 104108 (2020)] suggest that even the T2 condition leads to unacceptably inaccurate results in the case of two- and three-dimensional clusters. We demonstrate that these failures persist under the application of full three-particle N-representability conditions (3POS). A variety of correlation metrics are explored in order to identify regimes under which 3POS calculations become unreliable, and we find that the relative squared magnitudes of the cumulant three- and two-particle reduced density matrices correlate reasonably well with the energy error in these systems. However, calculations on other molecular systems reveal that this metric is not a universal indicator for the reliability of the reduced-density-matrix theory with 3POS conditions.
