  Deformed correlated Gaussian basis functions are introduced and their matrix
elements are calculated. These basis functions can be used to solve problems
with nonspherical potentials. One example of such potential is the dipole
self-interaction term in the Pauli-Fierz Hamiltonian. Examples are presented
showing the accuracy and necessity of deformed Gaussian basis functions to
accurately solve light-matter coupled systems in cavity QED.
