CombiningMany-Body Perturbation and Quantum Electrodynamics

It has been a long-sought problem to be able to combine many-body perturbation theory and quantum electrodynamics into a unified, covariant model. Such a model has recently been developed at our laboratory and is outlined in the present paper. The model has potential applications in many areas and opens up the possibility of studying the interplay between various interactions in different system. The model has so far been applied to highly ionized helium-like ions, and some numerical results are given. It is expected that the combined effect—that has never been calculated before—could have a significant effect on certain experimental data. The radiative effects are being regularized using the dimensional regularization in Coulomb gauge, and the first numerical results have been obtained.


Introduction
The many-body perturbation theory (MBPT) is now highly developed, and particularly all-order methods like the coupled-cluster approach (CCA) are very efficient and frequently used in atomic and molecular calculations [1,2].There is nowadays a great research interest also in highly charged ions that can be effectively produced in modern high-energy accelerators.The purpose in this type of research is largely to study quantum electrodynamics (QED) in the presence of strong fields [3].The experimental accuracy can here be quite high-in some cases far beyond the reach of any theoretical models presently available [4,5].The shortcoming of these models is expected to be mainly due to the omission of the combined effect of QED and electron correlation.
The standard methods for MBPT can treat the electron correlation to essentially all orders of perturbation theory, but since they are noncovariant they cannot incorporate QED effects in any systematic way.Several covariant methods have been developed particularly for QED calculations on bound atomic systems, and most frequently used is the Smatrix formulation [6,.These methods, however, have the disadvantage that they cannot be combined with MBPT, and they are for practical reasons limited to twophoton exchange.Particularly for lighter systems, this yields an incomplete treatment of the electron correlation.
It has been a long-sought problem to be able to calculate combined MBPT and QED effects in a systematic way, which could have important implications in various areas.The main obstacle, however, has been that the MBPT and QED procedures have quite different structures.(In 1988 a research programme at the Institute of Theoretical Physics in Santa Barbara, "Relativistic, Quantum electrodynamic, and Weak Interaction Effects in Atoms" was largely devoted to this problem [7].)One of the methods for QED calculations; however, namely, the Covariant-Evolution-Operator (CEO) method, recently developed by the Gothenburg group [8,9], has a structure that is similar to that of MBPT, and it can therefore serve as a basis for a unified procedure.Such a procedure has been developed by the group during the last few years and is described in a number of publications [10,11] as well as in a recent book [6].The procedure is now being implemented [12][13][14][15], and at a later stage interesting comparison can be made with accurate experimental data that presently cannot be explained.The developments made so far are summarized in the present article.

Time-Independent Many-Body Perturbation Theory
As an introduction we shall consider the basic many-body perturbation theory (MBPT), as described, for instance, in the book by Lindgren and Morrison, Atomic Many-Body Theory [1].
2.1.Bloch Equation.Many MBPT calculations in physics and quantum chemistry are based upon the generalized Bloch equation.We have then a set of target states satisfying the Schrödinger equation and for each target state there exist model states, |Ψ α 0 (α = 1 . . .d), which are assumed to be linearly independent and form a model space.The projection operator for the model space is P, which together with the projection operator for the complementary space Q, forms the identity operator ( A single wave operator, Ω, transforms all model states to the corresponding target states, ( An effective Hamiltonian, operating in the model space, generates the exact energy when operating on the model state This leads to which according to the Schrödinger equation ( 1) is equal to HΩ|Ψ α 0 .These relations hold for the entire model space, leading to the generalized Bloch equation In intermediate normalization (IN) the effective Hamiltonian (4) becomes The Hamiltonian is normally partitioned into where H 0 a model Hamiltonian, containing the sum of singleelectron Hamiltonians, and V is a perturbation.This leads to the commonly used form of the generalized Bloch equation [1,16] [Ω, where is the effective interaction.
In the graphical representation the wave operator is under general conditions (complete model space, see [1]) represented by linked diagrams, which is the linked-diagram theorem, first demonstrated for closed-shell systems by Brueckner [17] and Goldstone [18] and later generalized by Brandow [19] and Lindgren [16] to general open-shell systems.This can be expressed a form of the generalized Bloch equation frequently used as the starting point for many-body perturbation calculations in physics and quantum chemistry.It is valid also for a manydimensional model space that does not have to be degenerate.The last term is referred to as the folded part, since it is conventionally represented graphically by folded diagrams [19] (see Figure 1).

All-Order
Procedures.The Bloch equation can be iterated, leading to all-order procedures, as illustrated by the pair function (see Figure 2) Here are the resolvent and reduced resolvent, respectively, and E 0 the energy of the initial state.The all-order pair function satisfies a Dyson-type equation where V is the Coulomb interaction.This is represented graphically in Figure 3.The last term represents the folded contribution with a double resolvent.P E , P E , . . .are projection operators for parts of the model space with energies E , E , . .., respectively.The pair equation (13) represents the important pair-correlation to all orders and is frequently used in calculations on atomic and molecular systems.An even more powerful technique is the coupled-cluster approach (CCA), where the wave operator is expressed in exponential form This has the effect, for instance, that including pair correlation into the cluster operator, yields the most important quadruple effects on the wave operator.This technique has [t!] been extensively used during the last decades, particularly in quantum chemistry.For a review the reader is in particular referred to a recent review [2].For open-shell systems it is convenient to apply the normal-ordered exponential Ansatz [20] where the curly brackets represent normal ordering.In this way spurious contractions between open-shell operators are avoided.
In the CCA all diagrams are "connected" (for a complete model space), and the expansion can be expressed in analogy with (12) as (For the exact definitions of "linked" and "connected" diagrams, see, for instance, [1].) where the three terms represent the sum of single-particle Dirac Hamiltonians, Coulomb, and instantaneous Breit in-teractions between the electrons, respectively.The projection operators, Λ, are inserted in order to avoid negative-energy states that can lead to singularities (Breit-Ravenhall disease [22]).This approximation includes all effect up to order α 2 Rydbergs (or α 4 mc 2 ), where α is the fine-structure constant.
Effects beyond this approximation-of order α 3 Rydbergs and higher-are conventionally referred to as QED effects (see further below).

Covariant-Evolution Operator
Mainly three methods for QED calculations on bound states have been developed (see [9, Section 5], [6,Chapter 6]).The standard procedure is the well-known S-matrix formulation, and more recently two other methods have been developed, the two-times Green's function technique, developed by the St Petersburg group [23], and the covariant-evolutionoperator (CEO) method, developed by the Gothenburg group [9].The CEO method is particularly suitable for combining with MBPT and will be briefly summarized below.

Definition.
The standard time-evolution operator, U(t, t 0 ), describes in the interaction picture the evolution of the nonrelativistic time-dependent state vector where N is a normalization constant.It should be noted that the evolution operator generally does not preserve the normalization, and in addition it can contain singularities.Therefore this equation should be handled with care.The operator can be expanded (Relativistic units are used: where an adiabatic damping is inserted in order to handle the singularities.The damping factor, γ, eventually goes to zero.T is the Wick-time-ordering operator and H (x) is the perturbation density corresponding to the time-dependent perturbation We note that this represents the emission/absorption of a single photon and, hence, operates in the extended Fock space with a variable number of photons.Two such interactions are needed to form the exchange of a virtual photon between the electrons (see Figure 4).The field-theoretical Green's function can in the closedshell case be defined where ψ H , ψ † H are the electron-field operators in the Heisenberg representation.The state |0 H is the "vacuum in the Heisenberg representation," that is, the state in the Heisenberg representation with no particles or holes.The numerator is where 0| is the vacuum state in the interaction picture and S = U(∞, −∞) is the S-matrix.This leads to the expansion The standard evolution operator (20) operates only in the forward direction and is graphically represented by diagrams with outgoing particle lines (of positive energy).The operator for single-photon exchange between the electrons is illustrated in the first diagram in Figure 5.This operator is not covariant and consequently cannot be used in a relativistic treatment.Green's function (24), on the other hand, with electron propagators on the free ends is covariant (second diagram).By attaching free-electron lines to the graphical representation of Green's function, we will represent an evolution operator that is covariant, the covariant-evolution operator (CEO; third diagram).
The covariant-evolution operator is the evolution operator for the relativistic state vector in analogy with ( 20) which also generally lies in the extended space (in the following we shall drop the subscripts "Cov" and "Rel").
The single-particle covariant-evolution operator can generally be expressed involving the same vacuum expectation as in the definition of the single-particle Green's function (24), apart from the fact that we shall here allow for uncontracted photons.The covariant-evolution operator for single-photon exchange, operating on the unperturbed state, |ab , of energy E 0 becomes (leaving out the damping factor) where Γ(E 0 ) is the resolvent (14).When only positive-energy states are involved (no virtual pairs), the potential V sp is given by where f (κ) is a gauge-dependent function of the magnitude of the photon momentum, κ = |k|.In contrast to the situation in the S-matrix formulation, the initial and final states can in the covariant-evolutionoperator formalism have different energies.When the initial and final states do have the same energy, the potential (30) above becomes In a ladder of single-photon-exchange, involving only positive-energy states and interacting to the far right on the unperturbed state of energy E 0 , the general potential becomes (see Figure 6) If there are no model-space states involved, the covariantevolution operator for the complete ladder becomes We note that the energy parameter of all potentials and resolvents as well as the time factor is the energy of the unperturbed state to the far right.

Green's
Operator.The covariant-evolution operator (CEO) becomes singular (or quasisingular) when a state degenerate (or quasidegenerate) with the initial state is involved.In the definition of Green's function in the closed-shell case (24) the corresponding singularities are eliminated by dividing by the expectation value of the S-matrix.In the CEO case a different procedure is needed.We define the Green's operator in the general Fock space by the relation where the operator to the left of the heavy dot does not operate beyond the dot.This implies that the energy parameter of Green's operator above is that of the intermediate modelspace state.The definition leads to the expansion (with As before, P E , P E , . . .are projection operators for the parts of the model space with the energies E , E , . .., and the relevant energy parameters of the operators are indicated.The negative "counterterms" cancel the singularities, making the Green's operator regular.We note that with the perturbation (22) the term P E U (1) (0, E )P E with only a single uncontracted photon vanishes.
It follows from the definition (27) that where is the "parent state."For the time t = 0 we have The target state |Ψ α is intermediately normalized (7) if where |Ψ α 0 is the corresponding model state.This leads to which is a generalization of the Gell-Mann and Low theorem [24], valid also for a multidimensional model space [9].Here, both the numerator and denominator are normally singular, but the ratio is regular in the limit when γ → 0.
Using the definition of Green's operator (34), we have from (36) since is the model state.This implies that Green's operator acts as a wave operator for the relativistic state vector at all times.In particular, the covariant analogue of the MBPT wave operator (3) becomes which is also a Fock-space operator.This gives the connection between the CEO formalism and standard MBPT.

The Many-Body Hamiltonian.
(See [6, Sections 6.4, 5].) The state vector (40), |Ψ α , satisfies a relativistic "Schrödinger-like" eigenvalue equation When the Coulomb gauge is used (as we shall normally assume to be the case here), the perturbation, V F , is given by the Coulomb interaction, V C , and the transverse part of the perturbation (23).The total many-body Hamiltonian then becomes in second quantization where v ext (x) is the external (usually nuclear) potential.This is Fock-space operator, operating in the extended space with unpaired photons.Since the number of photons is not conserved, also the radiation-field Hamiltonian, H Rad , is included.
Using Fock-space operators, the effective interaction (11) becomes The effective interaction still operates in the model space, which is a part of the restricted space with no unpaired photons.
Elimination of the singularity of the second-order evolution operator due to an intermediate model-space state (49) (c.f. Figure 1).

Model-Space Contribution.
(See [10,11], [6, Section 6.7].)The second-order Green's-operator expression (35) is for t = 0 using the fact that G (0) (0) = 1.We assume here that the perturbation is the complete photon exchange (32).The secondorder CEO then becomes This will be (quasi)singular, if the intermediate state lies in the model space.The singularity is cancelled by the counterterm, leaving a finite contribution, where δG (1) (0, E ) = G (1) (0, E ) − G (1) (0, E ) and δE = E − E .This is illustrated in Figure 7.In the case of complete degeneracy the difference ratio turns into a partial derivative.The finite contribution is referred to as the model-space contribution (MSC).If the interaction is energy independent, this becomes identical to the corresponding folded contribution in standard MBPT (Figure 1).Generalizing the procedure above to arbitrary orders, leads to [10, Equation (100)], [11,Equation (64)], [6, Equation (6.96)] where V eff is the effective interaction (11).The first term above represents Green's operator without any folds Consequently the last term of (50) represents the effects of all folds.
Operating with (50) on the model function in the case of a one-dimensional model space (single reference) of energy E 0 , gives where ΔE is the difference between the exact energy E and the unperturbed energy E 0 and V eff |Ψ 0 = ΔE|Ψ 0 .This is a Taylor expansion, (In the limit of complete degeneracy we have (see further [10,Appendix E]) and the result can be expressed implying that the MSC have the effect of shifting the energy parameters from the unperturbed to the perturbed energy.This relation is analogous to the relation between the Brillouin-Wigner and the Rayleigh-Schrödinger expansions in perturbation theory (see, for instance, [1,Chapter 9]).
It can be shown that analogous results hold also for a general model space and for arbitrary times.This implies that also the time dependence of the Green's operator is shifted in a similar way.In view of (33) the time dependence of the GO, operating on a model state Ψ α 0 then becomes e −it(E α −H0) (in interaction picture).Since according to (41) Green's operator has the same time dependence as the relativistic wave function, this result is in accordance with the standard quantum-mechanical picture.It then follows that Projecting onto the model space, yields using the definition (11).

Connection to the Bethe-Salpeter Equation.
(See [10], [6, Section 6.9].)We shall now show that the procedure presented here is compatible with the Bethe-Salpeter equation [24,25].We consider a single target state |Ψ (with energy E) and the model state |Ψ 0 (energy E 0 ).From the relations (43) and (53) we have which is a Brillouin-Wigner expansion of the wave operator.
We then have In analogy with the relation (53) we have for the effective interaction (11) where is the effective interaction without folds.Hence, But according to the definition ( 11) since PΩP = P in intermediate normalization and P commutes with H 0 .Combining this with (57) leads to the relation The treatment can be generalized by replacing the singlephoton exchange by the set of all irreducible multiphoton exchange, V (see Figure 8), leading to the "Schrödinger-like" equation in the restricted Hilbert space [10,Equation (113)], [6,Equation (6.126)] which is equivalent to the Bethe-Salpeter equation [24,25].
In the restricted space we have an equivalent Hamiltonian that is energy dependent.This we can compare with the Gell-Mann-Low Hamiltonian (44) that is energy independent and operates in the extended Fock space.

Combining MBPT with QED
4.1.General.(See [11], [6,Section 8.3].)We shall now describe how the covariant-evolution operator (CEO), developed in the previous section, can be used to combine MBPT and QED into a unified theory.
What is conventionally referred to as QED effects are effects that lie beyond the standard relativistic MBPT treatment, based upon the relativistic Dirac-Coulomb-Breit Hamiltonian (19).These effects are of three kinds, (i) retardation, (ii) virtual electron-positron pairs, and (iii) radiative effects (electron self-energy, vacuum polarization, and vertex correction).
The retardation effect is considered in the photon exchange above (30).Virtual pairs can be treated together with single-photon exchange by generalizing the potential to [6,Equation (8.11)] where t ± , and so forth, represent projection operators for particle/hole states, respectively.The upper or lower sign should be used consistently in each term, inclusive the sign in the front.The 16 combinations of upper and lower signs in the four term lead to the time-ordered diagrams, shown in Figure 9.
Vacuum-polarization effects on electron propagators can be represented by a potential, as discussed in several publications [26].
The single-photon exchange is in the covariant-evolution operator procedure represented by two single-particle interactions (23).The first-order electron self-energy can then be evaluated by closing the second interaction on the same electron (c.f. Figure 4).In both cases-before closing the photon-one or several instantaneous interactions (Coulomb or Breit) could be inserted, leading to the "QED potential", illustrated in Figure 10.It is assumed here that the selfenergy and vertex diagrams are properly renormalized (see further below).This potential contains at most one retarded photon, which is the maximum that for the being can  [t!] Graphical representation of the single-photon Bloch equation (66) (c.f. Figure 3).The last diagram represents the "folded" term, that is, the last term of the equation.This equation can be compared with the Bethe-Salpeter equation, valid only in the single-reference case, where there is no folded contribution.The order-by-order expansion of this equation is illustrated in Figure 11 [6, Figure 8.12].
be handled numerically.It should be noted that by including one or several instantaneous Breit interactions, most of the higher-order effects are actually included.This approximation is therefore quite accurate, as will be demonstrated below.
The QED potential can be iterated, as illustrated in Figure 11.The sum can be expressed in a Dyson-type of equation [6,Equation (8.68)] which is equivalent to a Bloch equation (10).Here, the last term represents the folded contribution.The asterisk indicates that the partial derivation should only include the last interaction.This equation can be generalized by replacing the single-photon potential by the irreducible multiphoton potential (Figure 8).The equation, illustrated in Figure 12, then represents a generalization of the Bethe-Salpeter equation, valid also for a multidimensional model space.Due to its resemblance with the standard Bloch equation of MBPT (10), it is referred to as the Bethe-Salpeter-Bloch equation.
The Bethe-Salpeter equation as it stands, however, is not suitable for numerical work, since already in second order (two-photon exchange) it is essentially beyond reach for computers of today.This problem can be circumvented by using the covariant-evolution operator with (to start with) one retarded photon and a number of instantaneous ones, which can easily be generated.We may then start with an iterated pair function, Ω I (13), and add a transverse photon, with or without crossing Coulomb interactions, and finish with new Coulomb iterations.This leads to the modified Bethe-Salpeter equation [11], [6,Equation (8.70)], illustrated in Figure 13.This procedure leads to much faster convergence and is already with a single retarded photon in most cases more accurate than the second-order BSE.

4.2.
Coupled-Cluster-QED Procedure.The procedure described above can also be combined with the coupled-cluster approach, mentioned earlier (Section 2.2).This leads to a CC-QED procedure, described in a recent publication [27].It is also briefly described in the book [6,Section 8.4].The pair functions of the classical procedure are then replaced by pair functions including QED effects, and the single-particle interactions will include electron self-energy and other single-particle effects.Since the classical and the QED procedures are completely compatible, the QED effects need to be included only in the components, where the effect is expected to significant.

Renormalization.
(See [6,Chapter 12].)The radiative QED effects involve divergences that have to be eliminated by regularization and renormalization processes.In order to take full advantage of the developments in MBPT, it is necessary to perform the calculations using the Coulomb gauge.The renormalization process is more complicated in the non-covariant Coulomb gauge than in the covariant gauges normally used in QED calculations.However, formulas have been derived by Adkins for the dimensional regularization in Coulomb gauge for the lowest-order free-electron self-energy and vertex correction [28,29].This procedure is now being implemented in our program structure [14,15].It has been tested on the self-energy of highly charged hydrogenlike ions, and the results agree very well with those obtained earlier using Feynman gauge.In addition, the Coulombgauge results exhibit higher numerical stability.To our knowledge this represents the first numerical dimensional regularization in the Coulomb gauge.

Implementation Procedure and Numerical Results
We illustrate the results of the numerical implementation of the procedures briefly described in the previous sections by the calculations performed by Hedendahl in his doctorial thesis [12] and illustrated in Figure 14.The implementation procedure and some results are also summarized in [6, Chapter 10] (more details will be given in a forthcoming publication [13]).
The picture shows the effect of electron correlation beyond two-photon exchange for the Coulomb-Breit interaction on helium-like ions.The results are normalized to the nonrelativistic ionization energy, and the vertical scale is logarithmic.The top line represents the total correlation effect in the no-pair approximation.This effect is largely present in relativistic MBPT (19).The second and third lines show the effect of correlation upon the retardation, without and with Coulomb crossings (without virtual pairs), and the last heavy line represents the correlation effect on the virtual pairs.These effects lie beyond all previously existing procedures and, therefore, have never been calculated.
As a comparison, the thin line at the bottom indicates the estimated effect of two retarded interactions without correlation.This is included in the second-order Bethe-Salpeter equation while the other effects are not.This comparison demonstrates clearly that-for light and medium-heavy elements-even a single retarded photon together with instantaneous Coulomb interactions contains substantial effects not present in second-order BSE.The accuracy of the combined MBPT-QED procedure can be improved even further by including the instantaneous Breit interaction.

Summary and Outlook
The present article summarizes the research activity of the Gothenburg Atomic-Theory Group over the last decade.For the first time a unified theory of MBPT and QED has been developed and implemented on highly charged helium-like ions.The procedure is compatible with the Bethe-Salpeter equation (BSE), which is exact but intractable for practical work and leads already in second order to almost insurmountable problems.The unified MBPT-QED procedure is a way around this problem.Here, a limited number of retarded interactions are combined with numerous instantaneous interactions (Coulomb and Breit).It has been demonstrated that already a single retarded interaction together with instantaneous interactions leads for light and medium-heavy elements to much higher accuracy than second-order BSE (Figure 14).This implies that the combined MBPT-QED procedure is superior to a perturbation expansion of the original BSE.Therefore, much more accurate calculations can be performed than have so far been possible.This opens up possibilities of performing calculations that was previously inaccessible.Then comparison can be made with very accurate measurements of energy separations [4,5], results that so far have not been explained.Obviously, the model proposed has implications also beyond the current project.

Figure 1 :Figure 2 :
Figure 1: (a) Folded diagrams appear in the standard (Goldstone) graphical formulation of MBPT and represent the finite remainder after eliminating singularities due to intermediate model-space state, |cd .In a relativistic treatment, Feynman diagrams are normally used, and then the corresponding diagrams will be drawn straight (right).(b) The second-order diagram is separated into a regular part with the intermediate state in the model space and a remainder, due to the intermediate model-space state, represented by the folded diagram.

Figure 3 :
Figure 3: Graphical representation of the all-order pair equation (15).The last diagram represents the "folded" term.The double line represents the double denominator (double resolvent).

Figure 4 :
Figure4: The single-photon exchange between the electrons is in the covariant-evolution operator represented by two perturbations(23).

Figure 5 :
Figure 5: Comparison between the standard evolution operator, Green's function, and the covariant-evolution operator for single-photon exchange in the equal-time approximation [6, Figure 8.1].

Figure 6 :
Figure 6: The evolution-operator diagram for single-photon exchange.

Figure 11 :
Figure 11: Iteration of the QED potential in Figure 10.This can be expressed by the self-consistent Dyson-type equation in Figure 12 [6, Figure 8.13].