4.2. Self-consistent QED-HF solver