We will discuss two complementary variational approaches to ai-QED that seek to find accurate eigenstates of the Pauli-Fierz (PF) Hamiltonian represented in the length gauge and within the dipole and Born-Oppenheimer approximations. Here we write down the Pauli-Fierz Hamiltonian for a molecular system coupled to a single photonic mode in atomic units as