Symmetry of the Relativistic Two-Body Bound State

We show that in a relativistically covariant formulation of the two-body bound state problem, the bound state spectrum is in agreement, up to relativistic corrections, with the non-relativistic bound state spectrum. The solution is achieved by solving the problem with support of the wave function in an O(2,1) invariant submanifold of the Minkowski spacetime. The O(3,1) invariance of the differential equation requires, however, that the solutions provide a representation of O(3,1). Such solutions are obtained by means of the method of induced representations, providing a basic insight into the subject of the symmetries of relativistic dynamics.


Introduction
In the nonrelativistic Newtonian-Galilean view, two particles may be thought of as interacting through a potential function V (x 1 (t), x 2 (t)); for Galiliean invariance, V must be a scalar function of the difference, i.e.,V (x 1 (t) − x 2 (t)). In such a potential model, x 1 and x 2 are taken to be at equal time, corresponding to a correlation between the two particles consistent with the Newtonian-Galilean picture.
For the relativistic theory, two world lines with action at a distance interaction between two points x µ 1 and x µ 2 cannot be correlated by the variable t in every frame. The Stueckelberg (SHP) theory [1] provides an effective and systematic way of dealing with the N body problem, and has been applied in describing relativistic fluid mechanics [2], the Gibbs ensembles in statistical mechanics and the Boltzmann equation [3], systems of many identical particles [4], and other applications.
The basic idea of the SHP theory is the parametrization of the world lines of particles with a universal parameter τ [5](see also [6] [7]). Stueckelberg [8] described classical pair annihilation with a world line that proceeds, in τ , in the positive direction of the time t (the observable time of Einstein [9]) and then passes to a motion in the negative direction of time for τ proceeding in its monotonic development, precisely as postulated by Newton [10] [11]. The transition is caused by interaction, such as emission of a photon. Although this process was considered to be classical, it occurs in a diagram in Feynman's perturbative expansion of the S-matrix [12].
He assumed the equations of motioṅ It then follows from (1) that the proper time ds 2 = −dx µ dx µ satisfies The theory implies that the particle mass m is a dynamical variable, reflecting the fact that the Einstein time t is an observable, and therefore that E = ± p 2 + m 2 , conjugate to t, must be an observable as well [5]. For m 2 = M 2 , (3) implies that the square of the proper time interval is equal to (dτ ) 2 , but in general, this relation cannot be maintained for non-trivial interaction. The Poisson bracket structure then follows from (2). The τ derivative of a function of x, p is given by With this,we see that [x µ , p ν ] P B = δ µ ν .
* For c → ∞, ct may remain finite for t → 0 and can be taken to be an arbitrary constant.
Following Dirac [13], it is assumed that the operator commutation relations, following the group action of translation implied by the Poisson bracket, as the basis for the construction Pof the quantum theory [5].
The correspondig Stueckelberg-Schrödinger equation is then taken to be, derived from the unitary evolution of the wave function ψ τ (x), with the operators p µ in K represented as −ih ∂ ∂x µ , self-adjoint in the scalar product (ψ, χ) = d 4 xψ τ * (x)χ τ (x). Eq. (7) corresponds to the quantum one particle problem. We now proceed to discuss the two body problem.

The two body bounPPPd state.
We review here the relativistic two body problem with invariant action at a distance potentials, for bound states.
As a candidate for an invariant action at a distance potential for the two body relativistic bound state we take for the potential V the function V (ρ), for where x µ 1 and x µ 2 are taken at equal τ , acting as a correlation parameter as well as the global generating parameter of evolution. This "relative coordinate" (squared) reduces to (x 1 − x 2 ) 2 ≡ x 2 at equal time for the two particles in the nonrelativistic limit, so that ρ becomes r in this limit (for simultaneous t 1 and t 2 ). Clearly, the solutions of a problem with this potential must then reduce to the solutions of the corresponding nonrelativistic problem in that limit.
The two body Stueckelberg Hamiltonian, is and the relative motion momentum The pairs P µ , X µ and p µ , x µ satisfy separately the canonical Poisson bracket (classically) and commutation relations (quantum mechanically), and commute with each other. Then where , and x = x 1 − x 2 . Both K CM and K rel are constants of the motion; the total and relative momenta for the quantum case may be represented by partial derivatives with respect to the corresponding coordinates. This problem was solved explicitly for the classical case by Horwitz and Piron [5], where it was shown that there is no precession of the typeP predicted by Sommerfeld [14], who used the nonrelativistic form 1/r for the potential (and obtained a period for the precession of Mercury that does not fit the data). The corresponding quantum problem was solved by Cook [15], with support for the wave functions in the full spacelike region; however, he obtained a spectrum of the form 1/(n + 1 2 ) 2 , with n an integer, that does not agree with the Balmer spectrum for hydrogen. Zmuidzinas [16], brought to our attention by P. Winternitz [17]), however, proved that there is no complete orthogonal set of functions in the full spacelike region, and separated the spacelike region into two submanifolds, in each of which there could be complete orthogonal sets. The region for which x 2 > t 2 , in particular, permits the solution of the differential equations corresponding to the problem posed by (2.2) by separation of variables and provides spectra that coincide, up to relativistic corrections, with the corresponding nonrelativistic probPPPlems with potentials depending on r alone. We shall call this sector the RMS (reduced Minkowski space) [18] [19].
We may see, moreover, that the RMS carries an important physical interpretation for the nature of the solutions of the differential equations by examining the appropriate variables describing the full spacelike and RMS regions. The full spacelike region is spanned by over all ρ from 0 to ∞, β in (−∞, ∞), φ in (0, 2π) and θ in (0, π). Separation of variables in this choice, however, leaves the variable β for last; the quantum number (separation constant) obtained in this way has no obvious physical interpretation. Moreover, as found by Cook [15], the resulting spectrum for the Coulomb type potential (proportional to 1/ρ) does not agree with the Balmer series.
On the other hand, the set of variables describing the RMS, running over the same range of parameters [16], cover the entire space within the RMS (for x 2 1 + x 2 2 > t 2 ). In this coordinatization, the separation constant for θ (at the last stage), which enters the radial equation and determines the corresponding spectrum, has the interpretation of the angular momentum quantum number ℓ(ℓ + 1).
As for (14), for β → 0, these coordinates become the standard spherical representation of the three dimensional space (at the "simultaneity" point t = 0, where ρ becomes r). Independently of the form of the potential V (ρ), one obtains the same radial equation (in ρ) as for the nonrelativistic Schrödinger equation (in r), and therefore the same spectra (the two-body mass squared) for the reduced Hamiltonian. We shall discuss the relation of these results to the energy spectrum after writing the solutions. We summarize in the following the basic mathematical steps.
Assuming the total wavefunction (for P → P ′ , a point on the continuum of the spectrum of the conserved operator P ) the evolution equation for each value of the total energy momentum of the system is then For the case of discrete eigenvalues K a of K rel . We then have the eigenvalue equation (cancelling the center of mass wave function factor and K CM on both sides) Using the O(3, 1) Casimir operator, in a way quite analogous to the the use of the square of the total angular momentum operator, the Casimir operator of the rotation group O(3) in the nonrelativistic case, we may separate the angular and hyperbolic angular degrees of freedom from the ρ dependence. There are two Casimir operators defining the representations of O(3, 1) [20][21] [22]. The first Casimir operator is the second Casimir operator 1 2 ǫ µνλσ M µν M λσ is identically zero for two particles without spin. Recalling that our separation into center of mass and relative motion is canonical, and that using the canonical commutation relations, one finds that Since so that Eq. (18) can then be written as Choosing the RMS variables as we have defined them in (15), and with corresponding to the definition of the nonrelativistic angular momentum L, and corresponding to the boost generator A, We then find that where is the Casimir operator of the O(2, 1) subgroup of O(3, 1) leaving the z axis (and the RMS submanifold) invariant [18]. In terms of the RMS variables that we have defined above, We now proceed to separate variables and find the eigenfunctions. The solution of the general eigenvalue problem (24) can be written with invariant measure in the L 2 (R 4 ) of the RMS dµ = ρ 3 sin 2 θ cosh βdρdφdβdθ.
To satisfy the φ derivatives in (30), it is necessary to take where we have indexed the solutions by the separation constant m. For the case m an integer, this is a double valued function. To be compatible with the conditions on the other factors, this is the necessary choice; one must use, in fact,Φ m (φ) for m ≥ 0 and Φ * m (φ) for m < 0.
It has been suggested by M. Bacry [23] that the occurrence of the half-integer in the phase is associated with the fact that the RMS is a connected, but not simply connected manifold. One can see this by considering the projective form of the restrictions assuring that the events are relatively spacelike, and assuring, in addition, that the relative coordinates lie in the RMS. Dividing (34) and (35) by t 2 , and calling the corresponding projective variables X, Y, Z, we have from (34) the exterior of the unit sphere in the projective space, and from (35), the exterior of the unit cylinder along the z-axis. Identifying the points at infinity of the cylinder, we see that this corresponds to a torus with the unit sphere imbedded in the torus at the origin. Such a topological structure is associated with half integer phase (e.g. [24]). We now continue with our discussion of the structure of the solutions. The operator Λ contains the O(2, 1) Casimir N 2 ; with our solution (2.23), we then have where n 2 is the separation constant for the variable β. The term (m+ 1 2 ) 2 must be replaced by (m − 1 2 ) 2 = (|m| + 1 2 ) 2 for m < 0. We study only the case m ≥ 0 in what follows. The remaining equation for Λ is then For the treatment of Eq (38), it is convenient to make the substitution so that −1 ≤ ζ ≤ 1. One then finds that for The solutions are the associated Legendre functions of the first and second kind (Gel'fand [21]; see also Merzbacher [25]), P n m (ζ) and Q n m (ζ). The normalization condition on these solutions, with the measure (42) is The second kind Legendre functions do not satisfy this condition. For the condition on the P n m (ζ), it is simplest to write the known result [26] 1 −1 The normalized solutions (it is sufficient to consider n ≥ 0) may be written aŝ where m ≥ n. The case n = 0 must be treated with special care; it requires a regularization. For n = 0, the associated Legendre functions become the Legendre polynomials P m (ζ). In terms of the integration on β, the factor cosh β = (1 − ζ 2 ) −1/2 in the measure is cancelled by the square of the factor (1 − ζ 2 ) 1/4 in the norm, so that the integration appears as The Legendre polynomials do not vanish at ζ = ±1, so ifB m and P m are related by a finite coefficient, the integral would diverge. When n goes to zero, associated with the ground state, the wave function spreads along the hyperbola labelled by ρ, going asymptotically to the light plane; the probability density with respect to intervals of β becomes constant for large β. The (regularized) expectation values reproduce the distribution of the Schrödinger bound states, although the spacetime wave function approaches that of a generalized eigenfunction.
To carry out the regularization, we take the limit as n goes continuously to zero after computation of scalar products. Thus, we assume the form with ǫ → 0 after computation of scalar products. This formula is essentially a residue of the Rodrigues formula for n → 0. The operator for the differential equation (24) for the eigenvalue of the reduced motion is invariant under the action of the Lorentz group. It follows from acting on the equation with the unitary representation of the Lorentz group that the eigenfunctions must be representations of that group [24] for each value of the eigenvalue. However, as one can easily see, the solutions that we found are, in fact, irreducible representations of O(2, 1), not, a priori, representations of the Lorentz group O(3, 1). We discuss below how to construct such a representation.
We have required that the wave functions be eigenfunctions of the Casimir operator (29) of the O(2, 1) subgroup. For the generators of O(2, 1), we note that It then follows that H ± are raising and lowering operators for m on the functions where it is convenient to replace m by n + k. With the relation one can show [19] that and that The complex conjugate of χ −n n+k transforms in a similar way, resulting in a second (inequivalent) representation of O(2, 1) with the same value of the O(2, 1) Casimir operator (these states correspond to replacement of m + 1 2 by m − 1 2 for m < 0, and are the result of charge conjugation. Since the operators A 1 , A 2 and L 3 are Hermitian, complex conjugation is equivalent to the transpose. Replacing these operators by their negative transpose (defined by C), leaves the commutation relations invariant. Thus the action on the complex conjugate states involves However, attempting to apply this method to the representations of O(2, 1) fails because the application of the Lie algebra to this set connects the lowest state of the tower with the ground state which, as we have shown, requires regularization. The action of the algebra does not provide such a regularization, and therefore the method is inapplicable.
We therefore turn to the method of induced representations [27]. We may apply this method to contructing the representations of O(3, 1) based on an induced representation with the O(2, 1) "little group", based on a spacelike vector corresponding to the choice of the z axis. We shall discuss his method in detail below.
We first record the solutions of the equation (18).
and the functionsΘ Eq. (39) becomes where we have defined The solutions are proportional to the associated Legendre functions of the first or second kind, P n ℓ (ξ) or Q n ℓ (ξ). For n = 0, the second kind functions are not normalizable. We therefore reject these.
The normalizable irreducible representations of O(2, 1) are single or double valued, and hence m must be integer or half integer. As we have seen, k is integer valued, and therefore n must be integer or half integer also. Normalizability conditions on the associated Legendre functions then require that ℓ be respectively, positive half-integer or integer. The lowest mass state, as we shall see from the spectral results, corresponds to ℓ = 0, and hence we shall consider only integer values of ℓ.Therefore, n and m must be integer.
We now turn to the solution of the radial equations, containing the spectral content of the theory. With the evaluation of Λ in (57), we may write the radial equation as If we put Eq. (58) becomes precisely the nonrelativistic Schrödinger equation forR (a) in the variable ρ, with potential V (ρ) (the measure for these functions is, from (32), just ρ 2 dρ, as for the nonrelativistic theory)

The spectrum
The lowest eigenvalue K a , as for the energy in the nonrelativistic Schrödinger equation, corresponds to the ℓ = 0 state of the sequence ℓ = 0, 1, 2, 3, ..., and therefore the quantum number ℓ plays a role analogous to the orbital angular momentum. This energy is of a lower value than achievable with wave functions with support in the full spacelike region [15] and the relaxation of the system to wave functions with support in the RMS may be thought of, in this sense, as a spontaneous symmetry breaking (we thank A. Ashtekar for his remark on this point [28]).
The value of the full generator K is then determined by these eigenvalues and the value of the center of mass total mass squared operator, i.e., The first term corresponds to the total effective rest mass of the system. In particular, the invariant mass squared of the system is given by (sometimes called the Mandelstam variable s [29]) s a ≡ −P 2 a = 2M (K a − K).
This total center of mass momentum is observed in the laboratory in scattering and decay processes, where it is defined as the sum of the outgoing momenta squared. In the case of two particles, it would be given by −(p µ 1 + p µ 2 )(p 1µ + p 2µ ), as we have defined it in (62). This quantity is given in terms of total energy and momentum by and in the center of momentum frame, for P = 0, is just E 2 T . In order to extract information about the energy spectrum,we must therefore make some assumption on the value of the conserved quantity K. In the case of a potential that vanishes for large ρ, we may consider the two particles to be asymptotically free, so the effective Hamiltonian in this asymptotic region Further, assuming that the two particles at very large distances, in accordance with our experience, undergo a relaxation to their mass shells, so that p 2 i ∼ = −M 2 i . In this case, K would be assigned the value The two particles in this asympotic state would, for the bound state problem, be at the ionization point. If these assumptions are approximately valid, we find for the total energy, which we now label where we have restored the factors c.
In the case of excitations small compared to the total mass of the system, we may factor out M c and represent the result in a power series expansion so that the energy spectrum is just the set {K a } up to relativistic corrections. Thus, the spectrum for the 1/ρ potential is just that of the nonrelativistic hydrogen problem up to relativistic corrections, of order 1/c 2 . If the spectral set {K a } includes large negative values, the result (66) could become imaginary, indicating the possible onset of instability. However, the asymptotic condition imposed on the evaluation of K must be re-examined in this case. If the potential grows very rapidly as ρ → 0, then at large spacelike distances, where the hyperbolic surfaces ρ = const approach the lightcone, the Euclidean measure d 4 x (thought of, in this context, as small but finite) on the R 4 of spacetime starts to cover very singular values and the expectation values of the Hamiltonian at large spacelike distances may not permit the contribution of the potential to become negligible; it may have an effectively very long range. This effect can occur in the transverse direction to the z axis along the tangent to the light cone; the hyperbolas cannot reach the light cone in the z direction, which may play an important role in the modelling the behavior of the transverse scattering amplitudes in high energy scattering studied, for example, by Hagedorn [30].

Some examples
In this section we give the examples of the Coulomb potential and the oscillator. For the analog of the Coulomb potential, we take As we have remarked above,for c → ∞, this potential reduces to 1/r, the usual Coulomb, and therefore the spectrum must reduce to the usual Balmer series in this limit. In this case the spectrum, according to the solutions above, is given by where n a = 0, 1, 2, 3..... The wave functionsR(ρ) a are the usual hydrogen functionŝ R n a ℓ (ρ) = Zn a ! (n a + ℓ + 1) 2 (n a + 2ℓ + 1) where L 2ℓ+1 n a are the Laguerre polynomials, and the variable x is defined by and a 0 =h 2 /me 2 ). The size of the bound state, which is related to the atomic form factor, is measured according to the variable ρ [31]. For the lowest level (using the regularized functions) n a = ℓ = 0, The total mass spectrum, given by (62), is then For the case that the nonrelativistic spectrum has value small compared to the sum of the particle rest masses, we may use the approximate relation (66) to obtain The lowest order relativistic correction to the rest energy of the two body system with Coulomb like potential is then For positronium, ∆(E − M c 2 ) ∼ 2 × 10 −5 eV it is about one part in 10 5 , about 2% of the positronium hyperfine splitting of 8.4 × 10 −4 eV [32]. We see quantitatively that the relativistic theory gives results that are consistent with the known data on these experimentally well studied bound state systems. For the four dimensional oscillator, with V (ρ) = 1 2 mω 2 ρ 2 , Eq.(60) takes the form With the transformationR we obtain the equation Normalizable solutions, the Laguerre polynomials L ℓ+1/2 n a (x), exist [18] when the coefficient of w (a) (x) is a negative integer, so that the eignvalues are where n a = 0, 1, 2, 3, . . . The total mass spectrum is given by (62) as Note that the "zero point" term is 3 2 , indicating that in the RMS, in the covariant equations there are effectively three intrinsic degrees of freedom, as for the nonrelativistic oscillator.
The choice of K is arbitrary here, since there is no ionization point for the oscillator, and no a priori way of assigning it a value; setting K = − M c 2 2 as for the Coulomb problem (a choice that may be justified by setting the spring constant equal to zero and adiabatically increasing it to its final value), one obtains, for small excitations relative to the particle masses, Feynman, Kislinger and Ravndal [33], Kim and Noz [34] and Leutwyler and Stern [35] have studied the relativistic oscillator and obtained a positive spectrum by imposing a subsidiary condition suppressing timelike excitations, which lead, in the formalism of annihilationcreation operators to generate the spectrum, to negative norm states ("ghosts"). There are no ghost states in the covariant treatment we discuss here, and no extra constraints invoked in finding the spectrum. The solutions are given in terms of Laguerre poynomials, but unlike the case of the standard treatment of the 4D oscillator, in which x µ ± ip µ are considered annihilation-creation operators, the spectrum generating algebra (for example, Dothan [36]) for the covariant SHP oscillator has been elusive [37].

The induced representation
We have remarked that the solutions of the invariant two body problem results in solutions that are irreducible representations of O(2, 1), in fact, the complex representations of its covering group SU (1, 1), and pointed out that the ladder representations generated by the action of the Lorentz group on these states cannot be used to obtain representations of the full Lorentz group O(3, 1) or its covering SL(2, C). Since the differential equations defining the physical states are covariant under the action of O(3, 1), the solutions must be representations of O(3, 1). To solve this problem, one observes [1] that the O(2, 1) solutions are constructed in the RMS which is referred to the spacelike z axis. Under a Lorentz boost, the entire RMS turns, leaving the light cone invariant. After this transformation the new RMS is constructed on the basis of a new spacelike direction which we call here m µ . However, the differential equations remain identically the same since the operator form of these equations is invariant undr Lorentz transformations.. The change of coordinates to RMS variables has the same form as well, and therefore the set of solutions of these equations have the same structure. These functions are now related to the new z axis. Under the action of the full Lorentz group the wave functions undergo a transformation involving a linear combination of the set of eigenfunctions found in the previous section; this action does not change the value of the SU (1, 1) (or O(2, 1)) Casimir operator; together with the change in direction of the vector m µ , they provide an induced representation of SL(2, C) (or O(3, 1) with little group SU (1, 1) in the same way that relativistic spin is a representation of SL(2, C) with SU (2) little group [27].
Let us define the coordinates {y µ }, isomorphic to the set {x µ }, defined in an accompanying frame for the RMS(m µ )), with y 3 along the axis m µ . Along with infinitesimal operators of the O(2, 1) generating changes within the RMS(m µ )), there are generators on O(3, 1) which change the direction of m µ ; as for the induced representations for systems with spin [27], the Lorentz group contains these two actions, and therefore both Casimir operators are essential to defining the representations,i.e., both and which is not identically zero, and commutes with c 1 .
In the following, we construct functions on the orbit of the SU (1, 1) little group representing the full Lorentz group; along with the designation of the point on the orbit, labelled by m µ , these functions constitute a description of the physical state of the system.
It is a quite general result that the induced representation of a noncompact group contains all of the irreducible representations. We decompose the functions along the orbit into basis sets corresponding to eigenfunctions for the O(3) subgroup Casimir operator L(m) 2 → L(L + 1) and L 1 → q that take on values that persist along the orbit; these solutions correspond to the principal series of Gel'fand [21]. These quantum numbers for the induced representation do not correspond directly to the observed angular momenta of the system. The values that correspond to spectra and wavefunctions with nonrelativistic limit coinciding with those of the nonrelativistic problem problem, are those with L half-integer for the lowest Gel'fand L level. The Gel'fand classification, the two Casimir operators take on the values c 1 = L 2 0 + L 2 1 − 1, c 2 = −iL 0 L 1 , where L 1 is pure imaginary and , in general, L 0 is integer or half-integer. In the nonrelativistic limit, the action of the group on the relative coordinates becomes deformed in such a way that the O(3, 1) goes into the nonrelativistic O(3), and the O(2, 1) into the O(2) subgroup in the initial configuration of the RMS based on the z axis.
The representations that we shall obtain, in the principal series of Gel'fand [21], are unitary in a Hilbert space with scalar product that is defined by an integration invariant under the full SL(2, C), including an integration over the the measure space of SU (1, 1), carried out in the scalar product in L 2 (R 4 ⊆ RMS(m µ )), for each m µ (corresponding to the orientation of the new z axis, and an integration over the measure of the coset space SL(2, C)/SU (1, 1); the complete measure is d 4 yd 4 mδ(m 2 − 1), i.e., a probability measure on R 7 , where y µ ∈ RMS(m µ ). The coordinate description of the quantum state therefore corresponds to an ensemble of (relatively defined) events lying in a set of RMS(m µ )'s over all possible spacelike {m µ }.
A coordinate system oriented with its z axis along the direction m µ , as referred to above, can be constructed by means of a coordinate transformation of Lorentz type (here m represents the spacelike orientation of the transformed RMS, not to be confused with a magnetic quantum number), For example, if we take a vector x µ parallel to m µ , with x µ = λm µ , then the corresponding y µ is λm 0 µ , with m 0 µ in the direction of the initial orientation of the orbit, say, the z axis. This definition may be replaced by another by right multiplication of an element of the stability group of mµ and left multiplication by an element of the stability group of m 0 µ , constituting an isomorphism in the RMS. The variables y µ may be parametrized by the same trigonometric and hyperbolic functions as in (15) since they span the RMS, and provide a complete characterization of the configuration space in the RMS(m µ ) that is universal in the sense that it is the same in every Lorentz frame. It is convenient to define the functions We can then define the map of the Hilbert spaces associate with each m µ in the foliation H m → H Λm such that the state vectors are related by the norm preserving transformation In the new Lorentz frame (with y = L(Λm)x), If φ m (x) is scalar under Lorentz transformation, so that (we assume no additional phase) it follows from (88) that The wave function φ Λ Λm (x) describes a system in a Lorentz frame in motion with respect to the frame in which the state is described by φ m (x), and for which the support is in the RMS((Λm) µ ). The value of this function at x in the new frame is determined by its value at Λ −1 x in the original frame; moreover, the subensemble associated with values of m µ over the orbit in the new frame is determined by the subensemble associated with the values of (Λ −1 m) µ in the old frame. We define the description of the state of the system in the new frame in terms of the set (over {m µ }) of transformed wave functions where we have used (88) (the transformed function has support oriented with m µ ) and defined the (pseudo) orthogonal matrix (we define a "matrix" A as {A µ ν }) D(Λ, m) = L(m)ΛL T (Λ −1 m).
The transformation D −1 (Λ, m) stabilizes m 0 µ , and is therefore in the O(2, 1) subgroup that leaves the RMS of the original system invariant. Eq. (91) defines an induced representation of SL(2, C), the double covering of O(3, 1).
Classification of the orbits of the induced representation are determined by the Casimir operators of SL(2, C), defined as differential operators on the functions ψ m (y) of (86), i.e., the operators defined in (83) and (84). To define these variables as differential operators on the space {y}, we study the infinitesimal Lorentz transformations for which ψ 1+λ m (y) = ψ m−λm (D −1 (1 + λ, n)y), and λ is an infinitesimal Lorentz transformation (antisymmetric). To first order,the little group transformation is where d m is a derivative with respect to m µ holding y µ fixed, For the transformation of ψ m we then obtain ψ 1+λ m (y) ∼ = ψ m (y) − d m (λ + g m (λ))ψ m (y), where g m (λ) = G m (λ) µ ν y ν ∂ ∂y µ .
Eq. (99) displays explicitly the effect of the transformation along the orbit and the transformation within the little group. The algebra of these generators of the Lorentz group are investigated in [1]; the closure of this algebra follows from the remarkable property of compensation for the derivatives of the little group generators along the orbit (behaving in a way similar to a covariant derivative in differential geometry). The general structure we have exhibited here is a type of fiber bundle, sometimes called a Hilbert bundle, consisting of a set of Hilbert spaces on the base space of the orbit; in this case, the fibers, corresponding to these Hilbert spaces, transform under the little group O(2, 1).
There are functions on the orbit with definite values of the two Casimir operators, as well as L(m) 2 and L 1 (m); one finds the Gel'fand Naimark canonical represention with decomposition over the SU (2) subgroup of SL(2, C), enabling an identification of the angular momentum content of the representations [17]. With a consistency relation between the Casimir operators (for the solution of the finite set of equations involving functions on the hyperbolic parameters of the spacelike four vector m µ ), we find that we are dealing with the principal series of Gel'fand [20] [21].