Bound-State Variational Wave Equation For Fermion Systems In QED

Andrei G. Terekidi; Jurij W. Darewych; Marko Horbatsch

REVIEW

Bound-State Variational Wave Equation For Fermion Systems In QED

Andrei G. Terekidi, Jurij W. Darewych, Marko Horbatsch

High Energy Physics - Theory · Submitted 10 Apr 2006 · 19 references

Pith's verdict: A genuine variational route to the M-matrix kernel for two-fermion bound states, but the all-orders claim rests on an unproven adiabatic-stationarity step; the one-loop results are solid and the paper deserves a serious referee.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 28, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

T0 AI review An AI read the paper and answered a fixed list of questions.
T1 Audit Automatic checks that the references and files the paper relies on exist and match.
T2 Compute The paper's code or calculations re-run with the same results.
T3 Formal statement The claims written out in a formal math language a computer can read.
T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →

Challenge this review Comment (0) Read on arXiv ↗

Cite ▾

arXiv IDhep-th/0604078 (version 1)

Pith IDUWAWIJKD, Pith's permanent record of this paper and its review

CitationAndrei G. Terekidi, Jurij W. Darewych, Marko Horbatsch. "Bound-State Variational Wave Equation For Fermion Systems In QED." arXiv:hep-th/0604078 (2006). https://pith.science/paper/UWAWIJKD

Cite this reviewBibTeX for Pith's review of this paper

Claim paper

Why Pith says this

Highlight any text to comment on it.

Editor's note

read the letter

Dear X,

The paper deserves a serious look, but referee it with the main formal claim in mind: the all-orders M-matrix kernel is asserted, not proven. What is actually new is that Terekidi et al. obtain the two-fermion wave equation with the generalized invariant M-matrix kernel from the variational principle plus adiabatic hypothesis, rather than inserting it by hand as Zhang and Koniuk did. The radial equations for arbitrary mass ratio and the full state classification are worked out in detail, and the nonrelativistic limit correctly reduces to the Schrödinger equation. The one-loop radiative corrections for nS singlet states agree with Fulton and Martin; that is real evidence the machinery gets the right physics at O(alpha^5).

The soft spot is precisely the step the abstract calls rigorous. Equation (39) is derived from (32) under the assumption that the interaction-picture state is stationary, Eq. (35). Section 4 then sets t1=-infinity, t2=+infinity and inserts a trial state that is not an eigenstate of the evolution operator and contains no photon Fock components. The adiabatic switching (40) is invoked, but the epsilon, eta to zero limit is never taken, and no argument shows that the state at t=+infinity is related to the one at t=-infinity by a phase. Without that, the cancellation of disconnected diagrams and the identification of the kernel with the connected M-matrix do not go through. This is a genuine gap, not a nitpick. The one-loop application does not rescue it: that calculation uses first-order perturbation theory with the renormalized one-loop kernels, where the stationarity issue is much less relevant.

There are smaller issues. The substitution (124)-(126) for the renormalized vertex, vacuum polarization, and self-energy is imported from Ref. [16] without derivation, and the algebra from (121)-(123) to (127)-(128) is labeled straightforward but not shown. A referee will want those filled in. The quasi-spin label for the mixed states is a mild flourish, but harmless.

Who is this for? People working on relativistic two-body bound states in QED and on variational methods in QFT. It is a useful reference for how the M-matrix structure can emerge from a Hamiltonian approach, but it should be read with skepticism about the all-orders claim.

Verdict: send to peer review, yes. The flaws are addressable, the physics is right where it is tested, and the derivation, once repaired, would be a real contribution. In the current form the abstract overstates the rigor, and the adiabatic step needs either a properly taken limit or an explicit statement that it is a heuristic assumption. I would not cite it as a rigorous derivation until that is fixed, but I would want to see the revision.

Referee report

Significance

If the derivation of the all-orders wave equation were valid, the paper would provide a unified variational framework for relativistic two-fermion bound states, with the prospect of systematically including higher-order QED corrections beyond one loop. The paper has several genuine strengths: the reformulation of QED in terms of mediating-field Green's functions is explicit; the nonrelativistic reduction to the Schrödinger equation in Eq. (72) is demonstrated; the state classification via symmetry operators is systematic; and the O(alpha^5) energy shifts in Eqs. (129)-(130) agree with the existing literature. However, the central all-orders claim is not supported by the derivation as written, because the key stationarity assumption in Eq. (35) is not satisfied by the trial state used in Section 4. The significance of the paper is therefore conditional on a corrected derivation of the master variational equation.

Major objections (3)

[Section 3, Eqs. (35)-(39)] The passage from Eq. (32) to Eq. (39) relies on the assumption in Eq. (35) that U(t2,t1)|psi_i(t1)> = e^{i theta}|psi_i(t1)>, which is then used in Eq. (37) to replace U^{-1}(t2,t1) by a c-number. The trial state (41) used in Section 4 is a superposition of free two-fermion Fock states; it is not an eigenstate of the full interacting Hamiltonian, and it is not shown to be an eigenstate of U(+infinity,-infinity). The adiabatic switching ansatz (40) only makes the state at t=-infinity an eigenstate of the weak 'bait' Hamiltonian in the idealized epsilon,eta -> 0 limit, and this limit is never taken. Consequently, the denominator cancellation leading to Eq. (39) and the identification of the kernel with the all-orders M-matrix in Eq. (47) are not established. The authors should either prove the phase relation for the trial state in a well-defined adiabatic limit or rederive Eq. (39) using a standard linked-cluster expansion that does not require the trial state to be an eigenstate of U.
[Section 4 and Section 7, Eqs. (42), (47), (119)-(120)] The O(alpha^5) application does not test the all-orders equation (47). The calculation in Section 7 uses the nonrelativistic hydrogenic radial functions fJ(p) as zeroth-order input and employs first-order perturbation theory with only the one-loop kernels (53)-(56). It never uses the all-orders structure of Eq. (47) beyond these one-loop terms. Thus the agreement with Ref. [18] in Eqs. (129)-(130) is a perturbative consistency check and cannot validate the central claim that the kernel is the all-orders invariant M-matrix. This limitation should be stated explicitly in the paper.
[Section 7, Eqs. (121)-(128)] The derivation of the nonvanishing kernels (127)-(128) is omitted with the statement 'the calculations are straightforward', and the replacements (124)-(126) for the vertex, vacuum polarization, and self-energy functions are asserted with a reference to Ref. [16]. Since these kernels are the quantitative output of the paper, an outline of the trace reductions and the justification for retaining only the anomalous-magnetic-moment part of the vertex and dropping the self-energy contributions should be provided, for instance in an appendix.

Minor objections (5)

[Section 3, Eq. (40)] The order of the limits epsilon -> 0 and eta -> 0 and the precise sense in which the state at t=-infinity is an eigenstate of the baited Hamiltonian should be specified; the text simply says 'we put t1 = -infinity' without taking these limits.
[Section 4, Eq. (45)] The notation T(U(+infinity,-infinity) H_Ii(t) delta(t)) is confusing; it should be explained that the delta function fixes the time argument of the variational expectation to t=0, reducing the d^4x integral to a spatial integral.
[Appendix B, after Eq. (149)] The cross-reference 'expression (160)' appears to be a typo; it should refer to the displayed equation immediately above, likely Eq. (150).
[Section 5, around Eq. (84)] The quasi-spin quantum number t with t3 = +/- 1/2 is introduced but is not used anywhere in the subsequent derivations; it should either be used or removed.
[Throughout] There are several typographical errors, including 'Shrödinger' in the abstract, 'absorbtion' in Section 4, 'trial a tate' in Section 4, and 'McGrow-Hill' in Ref. [13]; a careful proofreading pass is needed.

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

The two-fermion bound-state kernel is shown to be the complete invariant M-matrix, summing all Feynman diagrams to every order.

The reading

This paper claims that the Hamiltonian variational method, combined with the adiabatic hypothesis, can turn QED's two-body bound-state problem into a single relativistic wave equation whose interaction kernel is the complete invariant M-matrix — the sum of all Feynman diagrams, reducible and irreducible. The derivation is carried out for two fermions of arbitrary mass ratio, and the equation is shown to reduce to the Schrödinger equation with a Coulomb kernel in the non-relativistic limit. If correct, the result provides a first-principles two-body framework that can include vacuum polarization, vertex corrections, self-energy, and multi-photon exchange at any order of the coupling constant, with a natural extension to three or more fermions. As a concrete demonstration, the authors compute one-loop radiative corrections to the binding energy for singlet l=J=0 states and recover the known O($alpha^{5}$) shifts.

What carries the argument

The load-bearing object is the generalized invariant M-matrix, M-tilde_{s1 s2 sigma1 sigma2}(p,q) = M^(1) + M^(2) + ... (Eq. 46), which assembles the one-photon tree diagram and all higher-order Feynman diagrams into a single kernel. The argument is carried by the identity in equations (38)-(39): the matrix element of T(U(t2,t1) H_I(t)) between interaction-picture states factorizes into connected and disconnected parts, and the disconnected part is exactly the denominator left over from the phase, so the variational equation keeps only connected diagrams. The adiabatic switching function $\alpha$(t) = $\alpha$((1-eta)$e^{{-epsilon |t-t0|}}$ + eta) (Eq. 40) prepares the trial state as an eigenstate of a weak 'bait' Hamiltonian at t = -infinity and transports it to the full interaction at t = 0, which is what licenses writing the state at early and late times as the same stationary state up to a phase.

Extended reading notes

Core claim

On the paper's own terms, the central discovery is equation (47): a variational equation for the two-fermion amplitude F_{s1 s2}(p) in which the free two-body energy (omega_p + Omega_p - E) is the kinetic term and the kernel is (-i) times the generalized invariant M-matrix M-tilde_{s1 s2 sigma1 sigma2}(p,q), defined by expanding the interaction-picture time-evolution operator and containing all Feynman diagrams, reducible and irreducible. The authors show that this equation follows from the variational principle after rewriting the expectation value in the interaction picture and using the factorization of matrix elements into connected and disconnected parts; the disconnected part cancels a denominator, leaving a connected-diagram-only variational equation. The second-order term M^(2) decomposes into vacuum polarization, vertex, self-energy, and two-photon box and crossed-box contributions, and the non-relativistic limit of the full equation is the Coulomb Schrödinger equation. The paper presents this as the first derivation, rather than postulation, of a two-body equation with an M-matrix kernel.

Load-bearing premise

The key load-bearing premise is that the trial states at two different times are stationary and related by a pure phase under the evolution operator — equation (35) — with the adiabatic switching (40) as justification; the paper never sends epsilon and eta to zero, so this phase relation is assumed rather than demonstrated.

Editorial extensions

If this is right

The two-fermion equation is in principle non-perturbative: the kernel contains all orders of the coupling constant, so the method is not limited to weakly coupled systems.
Higher-order radiative corrections to bound-state energies, including vacuum polarization, vertex corrections, fermion self-energy, and two-photon exchange, follow automatically from the variational principle rather than being inserted by hand.
The equation has the correct non-relativistic limit, reducing to the Schrödinger equation with Coulomb kernel, so standard hydrogen-like spectroscopy sits at the leading order of the formalism.
Positronium-like and muonium-like trial states yield radial equations for singlet, triplet, mixed-spin, and l-mixing states, and the one-loop O(alpha^5) singlet J=0 shifts reproduce earlier results.
The method is designed to generalize to three- and multi-fermion systems, which the authors state as a forthcoming application.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

I infer that the same variational machinery, with an N-fermion trial state, would produce N-body relativistic wave equations with an M-matrix-type kernel; the paper announces but does not demonstrate this generalization.
The paper leaves the adiabatic phase relation untested; a direct numerical check of (35) for a simple two-fermion trial state would either close the gap or point to a needed modification of the switching prescription.
The formalism suggests a comparison with the standard covariant bound-state treatment of two-body QFT: the M-matrix kernel can be viewed as a mass-shell projected, single-time reduction of the full covariant kernel, a connection the paper does not make.

Original abstract

We present a formulation of the Hamiltonian variational method for QED which enables the derivation of relativistic few-fermion wave equation that can account, at least in principle, for interactions to any order of the coupling constant. We derive a relativistic two-fermion wave equation using this approach. The interaction kernel of the equation is shown to be the generalized invariant M-matrix including all orders of Feynman diagrams. The result is obtained rigorously from the underlying QFT for arbitrary mass ratio of the two fermions. Our approach is based on three key points: a reformulation of QED, the variational method, and adiabatic hypothesis. As an application we calculate the one-loop contribution of radiative corrections to the two-fermion binding energy for singlet states with arbitrary principal quantum number $n$, and $l =J=0$. Our calculations are carried out in the explicitly covariant Feynman gauge.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 5,555 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.04 · states limitation 0.52 · human participants 0.01 · controlled comparison 0.01 · randomized assignment 0.01 · states quantitative result 0.10 · compares with prior results 0.06 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; show every role.

We present a formulation of the Hamiltonian variational method for QED which enables the derivation of relativistic few-fermion wave equation that can account, at least in principle, for interactions to any order of the coupling constant. We derive a relativistic two-fermion wave equation using this approach. The interaction kernel of the equation is shown to be the generalized invariant M-matrix including all orders of Feynman diagrams. The result is obtained rigorously from the underlying QFT for arbitrary mass ratio of the two fermions. Our approach is based on three key points: a reformulation of QED, the variational method, and adiabatic hypothesis. As an application we calculate the one-loop contribution of radiative corrections to the two-fermion binding energy for singlet states with arbitrary principal quantum number $n$, and $l =J=0$. Our calculations are carried out in the explicitly covariant Feynman gauge.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#	sentence	role	the authors' own work	states the main conclusion	strength
s0	We present a formulation of the Hamiltonian variational method for QED which enables the derivation of relativistic few-fermion wave equation that can account, at least in principle, for interactions to any order of the coupling constant.	contribution 0.94	0.94	0.45	asserted 0.66
s1	We derive a relativistic two-fermion wave equation using this approach.	contribution 0.77	0.97	0.16	asserted 1.00
s2	The interaction kernel of the equation is shown to be the generalized invariant M-matrix including all orders of Feynman diagrams.	result 0.95	0.90	0.28	asserted 1.00
s3	The result is obtained rigorously from the underlying QFT for arbitrary mass ratio of the two fermions.	result 0.47	0.94	0.13	asserted 1.00
s4	Our approach is based on three key points: a reformulation of QED, the variational method, and adiabatic hypothesis.	method 0.99	0.78	0.06	asserted 1.00
s5	As an application we calculate the one-loop contribution of radiative corrections to the two-fermion binding energy for singlet states with arbitrary principal quantum number $n$, and $l =J=0$.	contribution 0.61	0.96	0.10	asserted 1.00
s6	Our calculations are carried out in the explicitly covariant Feynman gauge.	method 1.00	0.93	0.03	asserted 1.00

What kind of paper the classifier reports this is: method or tool 0.52 · evidence from proof or derivation 1.00.

What was examined, and what was not

sources read

arXiv listing abstract v1

not examined

full text
figures and tables
references and citing sentences
relations to other papers

levels run

P abstract

not run

S sections
T body sentences
C citations
R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

References

Works cited (19)

[16]

T. Zhang and R. Koniuk, Can. J. Phys. 70, 683 (1992)

View on Pith 1992
[18]

T. Fulton, P. Martin Phys. Rev. 95, 811 (1954)

View on Pith 1954
[1]

A. G. Terekidi, J. W. Darewych, Journal of Mathematical Physic s 45, 1474 (2004)

View on Pith 2004
[2]

A. G. Terekidi, J. W. Darewych, Journal of Mathematical Physic s 46, 032302 (2005)

View on Pith 2005
[3]

J. W. Darewych, Annales Fond. L. de Broglie (Paris) 23, 15 (1998 )

View on Pith 1998
[4]

J. W. Darewych, in Causality and Locality in Modern Physics , G Hunter et al. (eds.), p. 333, (Kluwer, 1998)

View on Pith 1998
[5]

J. W. Darewych, Can. J. Phys. 76, 523 (1998)

View on Pith 1998
[6]

M. Barham and J. W. Darewych, J. Phys. A 31, 3481 (1998)

View on Pith 1998
[7]

B. Ding and J. Darewych, J. Phys. G 26, 907 (2000)

View on Pith 2000
[8]

J. D. Jackson, Classical Electrodynamics (John Wiley, New York, 1975)

View on Pith 1975
[9]

A. O. Barut, Electrodynamics and Classical Theory of Fields and Particl es (Dover, New York, 1980)

View on Pith 1980
[10]

W. T. Grandy, Relativistic Quantum Mechanics of Leptons and Fields (Kluwer, 1991)

View on Pith 1991
[11]

A. O. Barut, in Geometrical and Algebraic Aspects of Nonlinear Field Theor y, edited by S. De Filippo, M. Marinaro, G. Marmo and G. Vilasi, (Elsevier New York , 1989), p. 37

View on Pith 1989
[12]

M. Gell-Mann, F. Low, Phys. Rev. 84, 350 (1951)

View on Pith 1951
[13]

A. Fetter, J. Walecka, Quantum Theory of Many-Particle Systems (McGrow-Hill, 1971)

View on Pith 1971
[14]

J. W. Darewych and M. Horbatsch, J. Phys. B: At. Mol. Opt. 22, 973 (1989); 23, 337 (1990)

View on Pith 1989
[15]

W. Dykshoorn and R. Koniuk, Phys. Rev. A 41, 64 (1990)

View on Pith 1990
[17]

H. A. Bethe and E. E. Salpeter, Quantum Mechanics of One- and Two-Electron Atoms (Springer, 1957)

View on Pith 1957
[19]

T. Zhang and R. Koniuk, Can. J. Phys. 70, 670 (1992). 26

View on Pith 1992

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.

[1] [16]

T. Zhang and R. Koniuk, Can. J. Phys. 70, 683 (1992)

View on Pith 1992

[2] [18]

T. Fulton, P. Martin Phys. Rev. 95, 811 (1954)

View on Pith 1954

[3] [1]

A. G. Terekidi, J. W. Darewych, Journal of Mathematical Physic s 45, 1474 (2004)

View on Pith 2004

[4] [2]

A. G. Terekidi, J. W. Darewych, Journal of Mathematical Physic s 46, 032302 (2005)

View on Pith 2005

[5] [3]

J. W. Darewych, Annales Fond. L. de Broglie (Paris) 23, 15 (1998 )

View on Pith 1998

[6] [4]

J. W. Darewych, in Causality and Locality in Modern Physics , G Hunter et al. (eds.), p. 333, (Kluwer, 1998)

View on Pith 1998

[7] [5]

J. W. Darewych, Can. J. Phys. 76, 523 (1998)

View on Pith 1998

[8] [6]

M. Barham and J. W. Darewych, J. Phys. A 31, 3481 (1998)

View on Pith 1998

[9] [7]

B. Ding and J. Darewych, J. Phys. G 26, 907 (2000)

View on Pith 2000

[10] [8]

J. D. Jackson, Classical Electrodynamics (John Wiley, New York, 1975)

View on Pith 1975

[11] [9]

A. O. Barut, Electrodynamics and Classical Theory of Fields and Particl es (Dover, New York, 1980)

View on Pith 1980

[12] [10]

W. T. Grandy, Relativistic Quantum Mechanics of Leptons and Fields (Kluwer, 1991)

View on Pith 1991

[13] [11]

A. O. Barut, in Geometrical and Algebraic Aspects of Nonlinear Field Theor y, edited by S. De Filippo, M. Marinaro, G. Marmo and G. Vilasi, (Elsevier New York , 1989), p. 37

View on Pith 1989

[14] [12]

M. Gell-Mann, F. Low, Phys. Rev. 84, 350 (1951)

View on Pith 1951

[15] [13]

A. Fetter, J. Walecka, Quantum Theory of Many-Particle Systems (McGrow-Hill, 1971)

View on Pith 1971

[16] [14]

J. W. Darewych and M. Horbatsch, J. Phys. B: At. Mol. Opt. 22, 973 (1989); 23, 337 (1990)

View on Pith 1989

[17] [15]

W. Dykshoorn and R. Koniuk, Phys. Rev. A 41, 64 (1990)

View on Pith 1990

[18] [17]

H. A. Bethe and E. E. Salpeter, Quantum Mechanics of One- and Two-Electron Atoms (Springer, 1957)

View on Pith 1957

[19] [19]

T. Zhang and R. Koniuk, Can. J. Phys. 70, 670 (1992). 26

View on Pith 1992