“…We consider a molecular system composed of N e electrons and N n nuclei coupled to 2 N c quantized transverse field modes of an infrared Fabry–Pérot cavity. This light–matter hybrid system is described by the molecular Pauli–Fierz Hamiltonian in minimal-coupling form ,, Ĥ PF = prefix∑ i N e ( true p̂ ̲ i + e  true̲ false( r̲ i false) ) 2 2 m e + prefix∑ a N n ( true P̂ ̲ a − Q a  true̲ false( R̲ a false) ) 2 2 M a + V ( r̲ , R̲ ) + prefix∑ λ , k 2 N normalc ℏ ω k ( â λ k † â λ k + 1 2 ) In the first line, we have kinetic energy contributions with electronic and nuclear momentum operators, p̂ i = − i ℏ ∇ i , and, P̂ a = − i ℏ ∇ a , electronic and nuclear masses, m e and M a , elementary charge, e , as well as nuclear charges, Q a = Z a e , with charge number, Z a . The first term in the second line is the molecular Coulomb interaction potential V ( r̲ , R̲ …”