The BetheāSalpeter equation (named after Hans Bethe and Edwin Salpeter) [1] describes the bound states of a two-body (particles) quantum field theoretical system in a relativistically covariant formalism. The equation was first published in 1950 at the end of a paper by Yoichiro Nambu, but without derivation. [2]
Due to its generality and its application in many branches of theoretical physics, the BetheāSalpeter equation appears in many different forms. One form, that is quite often used in high energy physics is
where Ī is the BetheāSalpeter amplitude, K the interaction and S the propagators of the two participating particles.
In quantum theory, bound states are objects with lifetimes that are much longer than the time-scale of the interaction ruling their structure (otherwise they are called resonances). Thus the constituents interact essentially infinitely many times. By summing up, infinitely many times, all possible interactions that can occur between the two constituents, the BetheāSalpeter equation is a tool to calculate properties of bound states. Its solution, the BetheāSalpeter amplitude, is a description of the bound state under consideration.
As it can be derived via identifying bound-states with poles in the S-matrix, it can be connected to the quantum theoretical description of scattering processes and Green's functions.
The BetheāSalpeter equation is a general quantum field theoretical tool, thus applications for it can be found in any quantum field theory. Some examples are positronium (bound state of an electronā positron pair), excitons (bound states of an electronāhole pairs [3]), and mesons (as quark-antiquark bound states). [4]
Even for simple systems such as the positronium, the equation cannot be solved exactly, although in principle it can be formulated exactly. A classification of the states can be achieved without the need for an exact solution. If one of the particles is significantly more massive than the other, the problem is considerably simplified as one solves the Dirac equation for the lighter particle under the external potential of the heavier particle.
The starting point for the derivation of the BetheāSalpeter equation is the two-particle (or four point) Dyson equation
in momentum space, where "G" is the two-particle Green function , "S" are the free propagators and "K" is an interaction kernel, which contains all possible interactions between the two particles. The crucial step is now, to assume that bound states appear as poles in the Green function. One assumes, that two particles come together and form a bound state with mass "M", this bound state propagates freely, and then the bound state splits in its two constituents again. Therefore, one introduces the BetheāSalpeter wave function , which is a transition amplitude of two constituents into a bound state , and then makes an Ansatz for the Green function in the vicinity of the pole as
where P is the total momentum of the system. One sees, that if for this momentum the equation holds, which is exactly the Einstein energy-momentum relation (with the Four-momentum and ), the four-point Green function contains a pole. If one plugs that Ansatz into the Dyson equation above, and sets the total momentum "P" such that the energy-momentum relation holds, on both sides of the term a pole appears.
Comparing the residues yields
This is already the BetheāSalpeter equation, written in terms of the BetheāSalpeter wave functions. To obtain the above form one introduces the BetheāSalpeter amplitudes "Ī"
and gets finally
which is written down above, with the explicit momentum dependence.
In principle the interaction kernel K contains all possible two-particle-irreducible interactions that can occur between the two constituents. In order to carry out practical calculations one has to model it by choosing a subset of the interactions. As in quantum field theories, interaction is described via the exchange of particles (e.g. photons in quantum electrodynamics, or gluons in quantum chromodynamics), other than contact interactions the most simple interaction is modeled by the exchange of only one of these force-carrying particles with a known propagator.
As the BetheāSalpeter equation sums up the interaction infinitely many times from a perturbative view point, the resulting Feynman graph resembles the form of a ladder (or rainbow), hence the name of this approximation.
While in quantum electrodynamics (QED) the ladder approximation caused problems with crossing symmetry and gauge invariance, indicating the inclusion of crossed-ladder terms. In quantum chromodynamics (QCD) this approximation is frequently used phenomenologically to calculate hadron mass and its structure in terms of Bethe--Salpeter amplitudes and Faddeev amplitudes, a well-known Ansatz of which is proposed by Maris and Tandy [4]. Such an Ansatz for the dressed quark-gluon vertex within the rainbow-ladder truncation respects chiral symmetry and its dynamical breaking, which therefore is an important modeling of the strong nuclear interaction. As an example the structure of pions can be solved applying the Maris--Tandy Ansatz from the Bethe--Salpeter equation in Euclidean space [5].
As for solutions of any homogeneous equation, that of the BetheāSalpeter equation is determined up to a numerical factor. This factor has to be specified by a certain normalization condition. For the BetheāSalpeter amplitudes this is usually done by demanding probability conservation (similar to the normalization of the quantum mechanical wave function), which corresponds to the equation [6]
Normalizations to the charge and energy-momentum tensor of the bound state lead to the same equation. In the rainbow-ladder approximation this Interaction kernel does not depend on the total momentum of the BetheāSalpeter amplitude, in which case the second term of the normalization condition vanishes.
The Bethe--Salpeter equation applies to all kinematic region of the Bethe--Salpeter amplitude. Consequently it determines the amplitudes where the functions are not continuous. Such singularities are usually located when the constituent momentum is timelike, which are not directly accessible from Euclidean-space solutions of this equation. Instead one develop methods to solve these type of integral equations directly in the timelike region [7]. In the case of scalar bound states through a scalar-particle exchange in the rainbow-ladder truncation, the Bethe--Salpeter equation in the Minkowski space can be solved with the assistance of Nakanishi integral representation [8].
Many modern quantum field theory textbooks and a few articles provide pedagogical accounts for the BetheāSalpeter equation's context and uses. See:
Still a good introduction is given by the review article of Nakanishi
For historical aspects, see
For a more comprehensive list of first principles codes see here: List_of_quantum_chemistry_and_solid-state_physics_software
The BetheāSalpeter equation (named after Hans Bethe and Edwin Salpeter) [1] describes the bound states of a two-body (particles) quantum field theoretical system in a relativistically covariant formalism. The equation was first published in 1950 at the end of a paper by Yoichiro Nambu, but without derivation. [2]
Due to its generality and its application in many branches of theoretical physics, the BetheāSalpeter equation appears in many different forms. One form, that is quite often used in high energy physics is
where Ī is the BetheāSalpeter amplitude, K the interaction and S the propagators of the two participating particles.
In quantum theory, bound states are objects with lifetimes that are much longer than the time-scale of the interaction ruling their structure (otherwise they are called resonances). Thus the constituents interact essentially infinitely many times. By summing up, infinitely many times, all possible interactions that can occur between the two constituents, the BetheāSalpeter equation is a tool to calculate properties of bound states. Its solution, the BetheāSalpeter amplitude, is a description of the bound state under consideration.
As it can be derived via identifying bound-states with poles in the S-matrix, it can be connected to the quantum theoretical description of scattering processes and Green's functions.
The BetheāSalpeter equation is a general quantum field theoretical tool, thus applications for it can be found in any quantum field theory. Some examples are positronium (bound state of an electronā positron pair), excitons (bound states of an electronāhole pairs [3]), and mesons (as quark-antiquark bound states). [4]
Even for simple systems such as the positronium, the equation cannot be solved exactly, although in principle it can be formulated exactly. A classification of the states can be achieved without the need for an exact solution. If one of the particles is significantly more massive than the other, the problem is considerably simplified as one solves the Dirac equation for the lighter particle under the external potential of the heavier particle.
The starting point for the derivation of the BetheāSalpeter equation is the two-particle (or four point) Dyson equation
in momentum space, where "G" is the two-particle Green function , "S" are the free propagators and "K" is an interaction kernel, which contains all possible interactions between the two particles. The crucial step is now, to assume that bound states appear as poles in the Green function. One assumes, that two particles come together and form a bound state with mass "M", this bound state propagates freely, and then the bound state splits in its two constituents again. Therefore, one introduces the BetheāSalpeter wave function , which is a transition amplitude of two constituents into a bound state , and then makes an Ansatz for the Green function in the vicinity of the pole as
where P is the total momentum of the system. One sees, that if for this momentum the equation holds, which is exactly the Einstein energy-momentum relation (with the Four-momentum and ), the four-point Green function contains a pole. If one plugs that Ansatz into the Dyson equation above, and sets the total momentum "P" such that the energy-momentum relation holds, on both sides of the term a pole appears.
Comparing the residues yields
This is already the BetheāSalpeter equation, written in terms of the BetheāSalpeter wave functions. To obtain the above form one introduces the BetheāSalpeter amplitudes "Ī"
and gets finally
which is written down above, with the explicit momentum dependence.
In principle the interaction kernel K contains all possible two-particle-irreducible interactions that can occur between the two constituents. In order to carry out practical calculations one has to model it by choosing a subset of the interactions. As in quantum field theories, interaction is described via the exchange of particles (e.g. photons in quantum electrodynamics, or gluons in quantum chromodynamics), other than contact interactions the most simple interaction is modeled by the exchange of only one of these force-carrying particles with a known propagator.
As the BetheāSalpeter equation sums up the interaction infinitely many times from a perturbative view point, the resulting Feynman graph resembles the form of a ladder (or rainbow), hence the name of this approximation.
While in quantum electrodynamics (QED) the ladder approximation caused problems with crossing symmetry and gauge invariance, indicating the inclusion of crossed-ladder terms. In quantum chromodynamics (QCD) this approximation is frequently used phenomenologically to calculate hadron mass and its structure in terms of Bethe--Salpeter amplitudes and Faddeev amplitudes, a well-known Ansatz of which is proposed by Maris and Tandy [4]. Such an Ansatz for the dressed quark-gluon vertex within the rainbow-ladder truncation respects chiral symmetry and its dynamical breaking, which therefore is an important modeling of the strong nuclear interaction. As an example the structure of pions can be solved applying the Maris--Tandy Ansatz from the Bethe--Salpeter equation in Euclidean space [5].
As for solutions of any homogeneous equation, that of the BetheāSalpeter equation is determined up to a numerical factor. This factor has to be specified by a certain normalization condition. For the BetheāSalpeter amplitudes this is usually done by demanding probability conservation (similar to the normalization of the quantum mechanical wave function), which corresponds to the equation [6]
Normalizations to the charge and energy-momentum tensor of the bound state lead to the same equation. In the rainbow-ladder approximation this Interaction kernel does not depend on the total momentum of the BetheāSalpeter amplitude, in which case the second term of the normalization condition vanishes.
The Bethe--Salpeter equation applies to all kinematic region of the Bethe--Salpeter amplitude. Consequently it determines the amplitudes where the functions are not continuous. Such singularities are usually located when the constituent momentum is timelike, which are not directly accessible from Euclidean-space solutions of this equation. Instead one develop methods to solve these type of integral equations directly in the timelike region [7]. In the case of scalar bound states through a scalar-particle exchange in the rainbow-ladder truncation, the Bethe--Salpeter equation in the Minkowski space can be solved with the assistance of Nakanishi integral representation [8].
Many modern quantum field theory textbooks and a few articles provide pedagogical accounts for the BetheāSalpeter equation's context and uses. See:
Still a good introduction is given by the review article of Nakanishi
For historical aspects, see
For a more comprehensive list of first principles codes see here: List_of_quantum_chemistry_and_solid-state_physics_software