4.3. Variational Transformation based QED-HF solver