The two lowest-lying excited states of the NH-CQED-CIS with both photon frequency values are plotted as a function of bondlength for values between 1.3 and 2.7 Angstroms with increments of 0.1 Angstroms (see Ref. 53) in Figure . In addition to computing these polariton surfaces at the NH-CQED-CIS/cc-pVDZ level, we compute these polaritonic surfaces at the same values of the MgH+ bondlength R using a model 3-level Hamiltonian: