The interaction between molecular excitations and nanoconfined photons can produce the requisite strong interactions for polaritonic chemistry . Motivated by a desire to provide a realistic picture of the molecular structure under the influence of strong photonic interaction, there has been a recent surge in activity focused on merging ab initio molecular electronic structure theory with cavity quantum electrodynamics (ab initio CQED) to provide an accurate and predictive model of polaritonic chemistry . Such efforts include combining CQED with density functional theory (DFT) or it's time-dependent extension (TDDFT) , reduced density matrix mechanics , and wavefunction theory using the coupled-cluster ansatz . These approaches can provide access to potential energy surfaces, couplings, and other properties of interest for simulating the structure and reactivity of polaritonic chemical systems.