where the usual HF mean-field procedure is used to iteratively modify the electronic HF molecular orbitals and photon coefficients, eventually reaching a variational solution to the ground state energy of the molecule-cavity hybrid system. In practice, one can first obtain the bare molecular HF energy, E HF , outside the influence of the cavity and variationally optimize the photonic expansion coefficients of the partially evaluated PF Hamiltonian as follows 100