Group Theoretical Approach to Pseudo-Hermitian Quantum Mechanics with Lorentz Covariance and $c \rightarrow \infty $ Limit

We present in the article the formulation of a version of Lorentz covariant quantum mechanics based on a group theoretical construction from a Heisenberg-Weyl symmetry with position and momentum operators transforming as Minkowski four-vectors under the Lorentz symmetry. The basic representation is identified as a coherent state representation, essentially an irreducible component of the regular representation, with the matching representation of an extension of the group $C^*$-algebra giving the algebra of observables. The key feature of the formulation is that it is not unitary but pseudo-unitary, exactly in the same sense as the Minkowski spacetime representation. The language of pseudo-Hermitian quantum mechanics is adopted for a clear illustration of the aspect, with a metric operator obtained as really the manifestation of the Minkowski metric on the space of the state vectors. Explicit wavefunction description is given without any restriction of the variable domains, yet with a finite integral inner product. The associated covariant harmonic oscillator Fock state basis has all the standard properties in exact analog to those of a harmonic oscillator with Euclidean position and momentum operators of any `dimension'. Galilean limit of the Lorentz symmetry and the classical limit are retrieved rigorously through appropriate symmetry contractions of the algebra and its representation, including the dynamics described through the symmetry of the phase space.


I. INTRODUCTION
Our group had implemented a, quantum relativity symmetry, group theoretical formulation of the full dynamical theory of the familiar quantum mechanics with rigorous classical limit given as the Newtonian theory, obtained through a contraction of the relativity symmetry applied to the specific representation [1]. The latter is taken as essentially an irreducible component of the regular representation of H(3), the Heisenberg-Weyl group. The full quantum relativity symmetry, denotedG (3), can naturally be seen as a U(1) central extension [2] of the Galilean symmetry. H R (3) is (or is isomorphic to) its subgroup, left after the 'time-translation' is taken out. A H(3) representation is a spin zero, time independent representation ofG (3). The representation is really the one of the canonical coherent states [3][4][5]. The matching representation of the group C * -algebra [6,7], further extended to a proper class of distributions, gives the observable algebra as functions, and distributions, of position and momentum operators,X i = x i ⋆ andP i = p i ⋆, as given by the Weyl-Wigner-Groenewold-Moyal(WWGM) formulation [8][9][10][11]. The operators α(p i ⋆, x i ⋆) = α(p i , x i )⋆ act as differential operators on the wavefunctions on coherent state basis φ(p i , x i ) by the Moyal star-product α ⋆ φ; α ⋆ β⋆ = (α ⋆ β)⋆.X i andP i can be seen as operator coordinates of the quantum phase space [12,13], which has been argued to serve as a proper quantum model for the physical space [1,14]. We naturally seek a Lorentz covariant version of that with a c → ∞ contraction of the symmetry taking the Lorentz boosts to that of the Galilean ones [15]. Such a contraction is the mathematically rigorous way to look at the full approximation of a theory under a certain limit, from the symmetry theoretical perspective.
The relativity symmetry for the quantum theory is one of H R (1, 3), which fits well into the contraction chain, at least at the symmetry and coset space level [16,17]. It has been well known that from a group theoretical perspective, a general overcomplete coherent state basis can naturally be identified with points of the appropriate coset. The latter in our cases corresponds to something like the classical phase space. The formulation of a fully Lorentz covariant version of quantum mechanics, with position and momentum operatorsX µ andP µ transforming as Minkowski four-vectors, has been around since the early days of quantum mechanics. A naive thinking would be to represent those operators as x µ and −i ∂ x µ , respectively, acting on the wavefunctions ψ(x µ ) with the simple inner product giving the squared integral norm, and to take a unitary Schrödinger evolution under the Einstein proper time τ , which gives the Klein-Gordon equation as the τ -independent equation of motion.
Explicit group theoretical picture of that has been available since the sixties [18,19]. The truth is, in any theory of quantum mechanics with wavefunctions on Minkowski four-vector variable(s), the real symmetry behind the system is the H R (1, 3) group instead of only its Poincaré subgroup. There are, however, difficulties with the unitary theory, especially well illustrated in the covariant harmonic oscillator problem [18,20], which we show explicitly in the appendix.
Other than being of interest on its own, the harmonic oscillator problem is of great theoretical importance. For our usual quantum mechanics, we have the well appreciated close connection between the Fock states, as the eigenstates of the harmonic oscillator Hamiltonian, and the canonical coherent states. The set of Fock states is one of the most useful orthonormal basis for the Hilbert space and the latter, as the space of rapidly decreasing functions spanned by their wavefunctions, giving the states on which the position and momentum operators are truly Hermitian in a completely consistent formulation [8]. Upon a more careful inspection, the Fock states are simultaneous eigenstates of the number opera-torsN i , or equivalently ofX 2 i +P 2 i . The subspace spanned by the Fock states of a fixed eigenvalue n of the total number operator N i corresponds exactly to the space of symmetric n-tensors of the three-dimensional Euclidean space. In particular, the three n = 1 states transform exactly as components of a three-vector in a complexified Newtonian space.
A perfectly nice embedding of all that into the space spanned by Fock states of the Lorentz covariant harmonic oscillator problem should be expected to have a fully parallel structure of symmetric n-tensors in the (1 + 3)-dimensional Minkowski spacetime [20]. Unlike for the SO(3) symmetry, however, the noncompact nature of Lorentz SO(1, 3) symmetry means that the corresponding spaces for the symmetric n-tensors, as its irreducible representations, cannot be unitary.
Replacing the full unitarity of the irreducible representation of H R (1, 3) by a pseudounitarity exactly in line with the Minkowski spacetime may be a good direction to formulate a theory of covariant quantum mechanics [21]. The representation as one for the SO (1,3) subgroup would reduce to a sum of finite dimensional irreducible components each labeled by two integers, the n and a nonzero positive integer characterizing the spin independent Casimir invariant. The latter corresponds to one plus the rank of the symmetric n-tensors [22].
We have presented in Ref. [22] the complete set of Fock states wavefunctions of such a pseudo-unitary representation of H R (1, 3) symmetry on a space of rapidly decreasing functions, hence completely free from divergence in themselves as well as in the Lorentz invariant indefinite inner product. Formulation given there corresponds to writing the time coordinate as ict. Here, the representation is rather given in the form of the |p µ , x µ coherent states, with p µ and x µ being real Minkowski four-vectors, from our group theoretical grand framework [1,16,17].
We want to emphasize that quantum dynamics is a symplectic dynamics and the physical Hamiltonian is just one among the many general Hamiltonians with the generated Hamiltonian flows as symmetries of the phase space. It is the symplectic structure of the latter as fixed by the inner product, or the metric for the vector space or its projective space, that is really the key. The actual symmetries of a physical system of course correspond to Hamiltonian flows the generators of which commute with the physical Hamiltonian, with the generators giving the conserved physical quantities.
As a preparation, we first sketch the notion of pseudo-Hermiticity and pseudo-unitarity clearly in the next section. In Sec.III and IV below, we start with an explicit presentation of the regular representation, and its irreducible components, of the H(1, 3) group. A major part of that is also needed to formulate the c → ∞ contraction. Each such component is shown to give essentially the same physical theory of covariant quantum mechanics we present in detail on the coherent state basis, in the abstract form and in wavefunctions, with the Lorentz invariant indefinite inner product. The part that involves the inner product and the pseudo-Hermitian/pseudo-unitary nature of the theory is put in Sec.IV after the presentation of the space of state vectors as the Fock space for covariant harmonic oscillator.
Sec.V deals with the Lorentz to Galilean, c → ∞, contraction of the representation, i.e.
the retrieval of the 'nonrelativistic' limit, the part for the dynamics of which is left to the last subsection of the Sec.VI. The latter is first devoted to the WWGM framework or the observable algebra, focusing on the symmetry transformations and the dynamics as a specific case of such a symmetry flow, with the real parameter characterizing transformation corresponding to an evolution parameter, which is taken as the proper time in the case. In Sec.VII we give a brief description of contraction to the classical theory and conclude in the last section.

II. PSEUDO-HERMITICITY AND PSEUDO-UNITARITY
Pseudo-unitarity is about an inner product that is not positive definite, like the Minkowski metric. Lorentz covariant quantum theory with an indefinite inner product vector space of states was first introduced by Dirac and Pauli [23,24]. However, an explicit detailed formulation of quantum mechanics with a careful attention paid to the covariant and contravariant indices seems not to be available. More interest has been focused on quantum field theories, such as quantum electrodynamics (see Ref. [25] for a review). It has been a common strategy, especially in gauge theories since Gupta-Bleuler [26,27], to formulate a theory on such a Krein space [28] and then project it onto the 'physical' Hilbert space as the positive normed subspace (see also Ref. [20] for the harmonic oscillator case), retrieving a standard probability interpretation. Interest in the related subject matter for quantum mechanics has been brought back to popularity from works on the so-called pseudo-Hermitian quantum mechanics [29][30][31], which we, in a way, rediscovered in our work of Ref. [22].
Let us sketch pseudo-Hermitian quantum mechanics here. A naive direct picture starts with a Hamiltonian operatorÂ H that is not Hermitian with respect to the given inner product of the Hilbert space. If a Hermitian operatorη can be found such that the operatorÂ H is called pseudo-Hermitian andη the (pseudo-)metric operator [24,31]. We think metric operator is the more appropriate name than pseudo-metric operator, especially because in our case it is essentially the exact manifestation of the Minkowski metric η µν .
The interesting thing is that a new inner product η ·|· can be introduced with respect to which the operatorÂ H is really Hermitian, namelyÂ † η H =Â H for the Hermitian conjugation To be more specific, one can call it a η-Hermiticity. 1 We introduce the somewhat unusual notation for a reason. Since we are talking about a second inner product on the same vector space, we want the vectors, kets, to be independent of the inner products, while the sets of bras as functionals can be defined differently [22], giving the different Dirac brackets as the different inner products.
The new inner product is not required to be positive definite. More importantly, the evolution generated byÂ H is 'unitary' [31] in the sense that it preserves the inner product between any two states. In the case of an indefinite inner product, the transformations preserving it are truly represented by the pseudo-unitary, rather than unitary, matrices.
Adopting from the terminology of special relativity, we have states with norms that can be spacelike (+ve), timelike (-ve), or lightlike (0). For the nondegenerate case, explicitly, one can find a countable orthonormal basis, like the Fock basis in our case, with L vectors of the Note that the notion of pseudo-Hermiticity is a relative one.Â H is not Hermitian and is pseudo-Hermitian only with respect to the original inner product ·|· , for which the Hermitian conjugateÂ † H is defined as the operator satisfying Looking at the theory as a dynamical one with the physical Hamiltonian operator, the η ·|· inner product is the only one relevant. The inner product certainly gives a metric to the vector space and its projective space, which also fixes the symplectic structure. That is the meaning of the choice of the (nontrivial) metric operator. The bottom line is, two different inner products on the same vector space really make two different inner product spaces and we generally do not have any necessity to consider two different inner products for a single theory of quantum dynamics. Often time, as in Ref. [22], it is just that the simplest or the most familiar kind of inner product is the 'wrong' one, based on which one can construct 2 Ref. [31] restricted the term 'inner product' to positive-definite products, which is not within its mathematical definition. That is the source of many 'pseudo-' terminology as in 'pseudo-inner product' and 'pseudo-metric', which we see as unnecessary. Defining an absolute pseudo-Hermiticity for the otherwise Hermitian operators which generate pseudo-unitary transformations preserving the indefinite inner product could be quite sensible though.
the 'right' one more easily. In this case, η ·|· = ·|η|· , or equivalently η ·| = ·|η. Obviously, that is the same as ·| = η ·|η −1 , so the two sets of bras are really on the equal footing. The naive perspective that the inner product ·|· is more basic is only a consequence of the presentation. Furthermore, the reality of an eigenvalue for an η-Hermitian operator follows in the same way as for a usual Hermitian one so long as the norm, i.e. η-norm here, of a corresponding eigenstate is nonzero. The latter of course always holds on a Hilbert space or, equivalently, for a positive definite inner product.

III. THE IRREDUCIBLE REPRESENTATIONS OF H R (1, 3)
We give the Lie algebra for H R (1, 3) as where η µν = diag{−1, 1, 1, 1}. The choice of notation with Y µ corresponding essentially to spacetime position observables and E µ to energy-momentum observables is somewhat unusual. The reason for it should be clear from the analysis below. Notice that the generators are all taken to have no physical dimension, and the factor 2 corresponds to in the chosen units, which is at least convenient for the coherent state formulation [1]. In terms of the group element g(p µ , x µ , θ, Λ µ ν ), we have the group product (with the indices suppressed) The story is an extension of what has been done in Ref. [1,14] for The basis of that pseudo-unitarity is the indefinite Minkowski norm associated with the metric η µν extending the Euclidean δ ij [21,22]. In the case of H R (3), the representation is naturally an irreducible component of the regular representation of H(3), which all can be seen actually as physically equivalent. It comes naturally as wavefunctions in the coherent state basis, on which the observables are represented as differential operators, essentially those obtained from the WWGM framework. Details of all that for the case of H(3) group has been presented in Ref. [1].
We present first the results from a harmonic analysis of Heisenberg-Weyl groups adapted to our case of H(1, 3) [32]. We write the left regular representation in = 2 units as are the left-invariant vector fields. In an irreducible representation, the central generator I has to be represented by a multiple of identity. We write the one parameter series V ς where we have with θ (ς) = ςθ, hence in a form formally independent of ς.X L (ς) andP L (ς) are still SO(1, 3) vectors, and so are p (ς) and x (ς) . The (ς) index becomes completely dummy and analysis based on the new operators and parameters is independent of ς so long as we are looking only at a particular irreducible representation. One can even simply drop it. From a physics perspective, we have absorbed the value of ς by a choice of physical unit for measuring the observables corresponding to Y and E, here all in the unit of √ ς. For ς being negative 3 , we transformations, hence each is an independent orbit. That is to say SÔ = SO(1, 3). The fact is of paramount importance for unambiguously identifying the nature of the coherent states below. In view of the discussion above, we can see that for any of the V L ς (p (ς) , x (ς) , θ (ς) ) representation, we can simply write it in the simple notation V L (p, x, θ), like taking the ς = 1 case as a representative. That is essentially what has been done in Ref. [1] for the H(3) or H R (3) case. However, for the reason to be clear below, we keep the explicit ς-notation for the most part of the manuscript.
The standard approach is to introduce the abstract canonical coherent states as for Each group element can be identified with a point in the H R (1, 3)/SO(1, 3) coset space [14,17].X (ς) andP (ς) are operators on the abstract representation space H ς spanned by the vectors, and |0, 0 = |0 (ς) is a fiducial normalized cyclic vector corresponding to the points (0, 0, θ (ς) ) in the coset space, each of which is fixed under SO(1, 3) transformations.

METRY GENERATORS FROM THE FOCK STATES
The kind of operator representation of the four-vector observables given in Eq.(9) would be naively seen as Hermitian, hence the full representation of the H R (1, 3) group as unitary.
To be more careful, the unitarity of a representation is really to be defined with respect to the inner product assumed for the representation space. The operator representation sure is Hermitian with respect to the usual squared-integral inner product, (with bar denoting the complex conjugation), for the wavefunctions that vanish at infinity. In the equivalent formulation in terms of standard Schrödinger wavefunctions ψ(x µ ), that is exactly the unitary representation given explicitly first in 1966 [18,19]. The short-comings of the formulation are best seen in the covariant harmonic oscillator problem [18,20]. We illustrate them explicitly in the appendix. To illuminate the pseudo-Hermitian nature of our representation, we present in the following the pseudo-unitary Fock space and complete the coherent states representation, together with the appropriate Lorentz invariant integral inner product. For convenience, in this section we drop the (ς) and ς subscripts.
We start with summarizing the more transparent abstract algebraic results for the Fock space [20,22] in a better logic and notation. For the HamiltonianX µX µ +P µP µ , we consider Here, we introduce the † η notation for a yet unspecifiedη, not excluding that being trivial, requiring however the η-Hermiticity ofX µ andP µ . Note that the feasible inner product is still to be determined. The Fock states are simultaneous eigenstates of theN (µ) , and hence also theN = N (µ) , operators. The last one is of course Lorentz invariant.
withâ 0 |n 0 ; n 1 , n 2 , n 3 = 2 √ n 0 |n 0 − 1; n 1 , n 2 , n 3 , and the exact corresponding results forâ i |n with |0, 0 = |0 . Moreover, those are exactly the states defined earlier, obtained by an The right inner product to complete the familiar algebra of the problem is, however, nontrivial. While the operatorsâ µ ,â † η µ andN (µ) with the commutation relation of Eq. (14) give a convenient generalization of theâ i ≡â i ,â † i andN (i) system, insensitive to the metric signature yet having the right Lorentz transformation properties of the Fock state solutions built in,â † η 0 is not a naive Hermitian conjugate ofâ 0 . In fact, we can obtain from Eq. (16) that â µ · |· = ·|â † η µ · when the usual orthonormality m|n = δ mn is assumed. We need a new inner product defined as i.e.η = (−1) n 0 |n n|, with the corresponding η-Hermitian conjugation, Note that specifying the inner product for a complete basis uniquely defines the inner product over the whole space. One can easily see that (â µ ) † =â † η µ implies Hermiticity ofX i andP i operators, while from equating (â 0 ) † = X 0 + iP 0 † witĥ a † η 0 , we see thatX 0 andP 0 are anti-Hermitian with respect to the usual inner product, i.e. the one withη being the identity. However, all of the H R (1, 3) generators are represented by pseudo-Hermitian, or η-Hermitian, operators. The operators all have real spectra, as we show explicitly below. Explicitly,X µ ,P µ andĴ µν =X µPν −X νPµ (andÎ) all satisfy Eq. (1) in the place ofÂ H , and we have a pseudo-unitary representation with invariant inner product η ·|· . In particular, the coherent states are normalized as η 0|0 = η p µ , x µ |p µ , x µ = 1, hiding the inner product indefiniteness. We can see the latter either through the explicit use of Eq.(17), or directly from the fact that V(p µ , x µ ) in Eq. (10) is an η-unitary operator. From the definition ofη in Fock basis and Eq.(17) we obtain showing explicitly the metric operatorη is a direct manifestation of the Minkowski metric in the Krein space of our quantum theory, exactly as we are looking for [21].
As the state |0, 0 = |0 has zero expectation values for theX µ andP µ operators, we get The generic wavefunctions can be introduced as exactly in the form of Eq. (9) and We see that the abstract formulation from the set of canonical coherent states based on the H(1, 3) manifold and the one from the irreducible component of the regular representation are really the same one. Wavefunction of a coherent state labeled by A is given by where (x − x A ) 2 and (p − p A ) 2 are the Minkowski vector magnitude squares, and can be seen as a special case of Eq.(23), namely In particular, we have which is the Lorentz invariant symmetric Gaussian.
To obtain the inner product on the space of wavefunctions, one simply has to look for the proper resolution of the identity operator on the Krein space. We havê Therefore, the functional η ψ| is represented on the space of φ(p µ , x µ ) as with the very nontrivial integration measure. The inner product η ψ|φ is then given by Each of the basis functions φ n (p µ , x µ ), and hence any general φ(p µ , x µ ) in the spanned space, is formally divergent at timelike infinity of the four-vector variables. On the other hand, all e (x 0 ) 2 +(p 0 ) 2 , and hence all φ(p µ ,x µ ) e (x 0 ) 2 +(p 0 ) 2 , are rapidly decreasing functions. The factor e −(x 0 ) , which characterizes the class of functions. The integral is finite for all wavefunctions as finite linear combinations of the Fock state basis φ n . Using φ(p µ ,x µ ) e (x 0 ) 2 +(p 0 ) 2 as the wavefunctions cannot be correct, though. That would, for example, make the wavefunction for |0 not Lorentz invariant and mess up the right transformation properties of all those for the Fock states, described in Ref. [22].
Thinking further about ψ * (p i , x i , −p 0 , −x 0 ) as ψ * (p µ , x µ ), one can see in hindsight that the inner product expression is indeed exactly what it should be. Of course we have that here rigorously established.
One can now easily show that the Fock states wavefunctions φ n (p µ , x µ ) have the proper norm ±1, and therefore are non-divergent without restricting the domain. The analytical feature is much better than that of the unitary representation (see the appendix). Note that, other than having a different inner product with a nontrivial integration measure, our formulation in terms of the wavefunctions and differential operator representation of theX µ andP µ really are the same as the usual unitary one. That illustrates clearly that the basic observablesX µ andP µ , as well as other observables in the form of their real functions, all have the same eigenvalues and eigenfunctions. In particular, the spectra are real.

V. LORENTZ TO GALILEAN CONTRACTION
A contraction [34,35] of the Lorentz symmetry SO (1, 3), sitting inside the H R (1, 3), to the Galilean ISO(3) has been discussed in Ref. [17], together with the corresponding coset spaces of interest. The full (quantum) relativity symmetry group obtained by contraction is named H GH (3), with commutators among generators essentially given by Note that the full result for the other commutators beyond the J ij and K i set, originated from SO(1, 3), is essentially fixed by the requirement of having the Galilean K i -H and the Heisenberg X-P commutators. However, for the purpose here, the explicit contraction is to be implemented a bit differently. It is taken as the c → ∞ limit of In the contraction, K i as generators for the Galilean boosts are the basic starting point and we would like to be able to trace physics, including the relative physical dimensions of quantities, by considering the speed of light c as having a physical dimension. Introducing X i = 1 c Y i is to keep the same physical dimensions for X i and P i . However, the essence of the contraction scheme as a formulation to retrieve an approximate physical theory from a more exact one is really to implement the contraction at a representation level.
To implement the contraction on V L ς , or the matching V ς as a representation of the original H(1, 3), it is important to note that the original central charge generator I represented by ςÎ in V ς would give the representation of the contracted I ′ , which remains central, as ς c 2Î . For a sensible result, one needs to consider ς = c 2 χ with χ staying finite at the contraction limit, hence I ′ represented by χÎ (recall:Î is the identity operator). Therefore, V ς contracts into V χ . In another words, the V ς representation of the original H(1, 3), and the full H R (1, 3), , as well as of the full group.
For the c → ∞ limit of V L χ (p µ (ς) , x µ (ς) ), we have to consider first and take that to obtain (with ς = c 2 χ). The above are the basic set of operators acting on the functional space of φ(p (ς) , x (ς) ), with the variables properly rescaled to a new set of variables to match with the operators. There is also the exactly corresponding set of operators,X (χ)i ,P (χ)i ,T (χ) , andĤ (χ) , and V χ on the abstract Hilbert space which are helpful for tracing the proper description.
The proper labels for the states p µ (ς) , x µ (ς) at the contraction limit should be p and hence giving naively We have then, at least formally, The crucial quantities controlling the nature of the representation are the overlaps

given in Eq.(24), we have it as
to be taken at the c → ∞ limit. It holds e should be taken as giving a diverging energy observableĤ L (χ) for any finiteP L (ς)0 . Furthermore, for an Einstein particle of the rest mass m, i.e the particle in Einstein's theory of special relativity, where the neglected terms involve negative powers of c 2 . At the c → ∞ limit, it is indeed diverging. Even p 0 is diverging. That is the result of the rest mass as an energy. Hence, it sure suggests that we should take our variable e (χ) as infinite, and the 'non-relativistic' energy we are interested in is the kinetic energy p i p i 2m given by the limit of e − mc 2 . Taking that feature into our consideration, the Hilbert space of interest under the contraction is really only the space spanned by the H(3) coherent states p for a fixed time t (χ) and a formally infinite e (χ) . To be exact, we should be implementing that logic from an Einstein particle to our quantum observablesĤ (χ) ,P 0 (χ) , andP (χ)i or their expectation values, but the conclusion is the same. The coherent state wavefunction A , hence at the contraction limit there is no more dependence on t (χ) and e (χ) , reducing it essentially to just ) . The operatorT L (χ) acts on the space of wavefunctions only as a multiplication by t (χ) and is just like classical, whileĤ L (χ) is not physically relevant. Note that the full contracted representation is then simply unitary. The part of the inner product η ·|· independent of p 0 (ς) and x 0 (ς) , hence t (χ) and e (χ) , is exactly the usual one, i.e.η essentially reduces to identity under the contraction. The space of wavefunctions spanned by φ A (p i (χ) , x i (χ) ) is a Hilbert space.

VI. GROUP THEORETICALLY BASED WWGM FRAMEWORK WITH WAVE-FUNCTIONS IN COHERENT STATE BASIS
The above analysis gives a successful picture of the phase space of the H R (1, 3) theory, giving in the Galilean limit the phase space of the H R (3) theory at each fixed 'time' value.
The infinite dimensional manifolds give, at the proper relativity symmetry contraction limit, the familiar finite dimensional classical models as approximation. The explicit results of the classical limit for the present case is presented in the section below. The merit of our group theoretical approach is that it gives a full dynamical theory associated with the corresponding spacetime/phase space model for each relativity symmetry, mutually connected through the contraction/deformation pattern. The dynamical theory is naturally a Hamiltonian theory from the symmetry of the phase space as symplectic geometry. The dynamics is better described on the algebra of observables as essentially the matching representation of the group C * -algebra [1,14,21]. Moreover, all those fit in well with the idea of the position and momentum operators as noncommutative coordinates of the phase space [12,13,21].

A. The Algebra of Observables, Symmetries, and Dynamics
The algebra of observables is depicted essentially as the one from a WWGM formalism, as functions and distributions of the position and momentum operatorsX µ andP µ . The basic dynamical variables of our representation on the space of wavefunctions φ(p µ (ς) , x µ (ς) ) arê X L = x + i∂ p = x⋆ andP L = p − i∂ x = p⋆, where we have dropped the µ indices and the subscript (ς) . We may also write a general function of (p µ (ς) , x µ (ς) ) as simply α(p, x), and the ⋆ is as in the Moyal star product with α(p, x)⋆ = α(p⋆, x⋆). Under such notation, the story looks quite the same as the case for H R (3) with onlyX L i , andP L i as x i ⋆ and p i ⋆, given in details in Ref. [1]. Hence, we present here only a summary of the results, leaving the readers to consult the latter paper and references therein.
Let us take a little detour first to clarify our theoretical perspective. What we have is rather like the WWGM put up-side-down [1]. We start with the quantum theory as an irreducible representation of a (quantum) relativity symmetry, including the Heisenberg-Weyl ⋆ is a unitary operator on K, generated by the Hermitian operator in the form of a real function G s (p⋆, x⋆), andŪ ⋆ ⋆ is its inverse obtained by the complex conjugation andŨ ⋆ ∈ O(K R ). We refer to the U ⋆ ⋆ as star-unitary, in particular whenever necessary to highlight it being a function of the p⋆ and x⋆ operators.
The above, illustrated for the case of H R (3) formulation of standard quantum mechanics in Ref. [1], can be applied to our H R (1, 3) case with a slight modification. We need to use the invariant inner product withρ φ = |φ η φ| for normalized |φ , and replace the Hermitian and unitary requirements by η-Hermitian and η-unitary ones. Our relevant symmetry transformations are to be given by η-unitary operator V ⋆(s) ⋆ generated by η-Hermitian G s (p⋆, x⋆), which are real functions of the basic η -Hermitian operators (p⋆, x⋆), i.e. G s (P L µ ,X L µ ) = G s (P L µ ,X L µ ), and we use theᾱ to denote the 'complex conjugate' of α as a function which correspond toᾱ⋆ as the η-Hermitian conjugate of α⋆ as an operators as an element of the observable algebra. The conjugation is the involution of the latter as a * -algebra. V ⋆(s) ⋆ of a η-unitary V ⋆(s) ⋆ is to be interpreted in the same manner. The feature of V ⋆(s) ⋆ to be the inverse of V ⋆(s) ⋆ is exactly η-unitarity. Again, η-Hermiticity is the Hermiticity so long as the algebraic analysis is concerned. Though η-unitarity here is really pseudo-unitarity, only the inner product preserving nature of it is relevant here and it is as where ρ(p, x) ∈K and {·, ·} ⋆ is the Moyal bracket. Hence, with ρ(s) =Ṽ ⋆(s) ρ(s = 0), The equation is the Liouville equation of motion for a mixed state ρ inD, the self-dual cone of K. The class of operators onK representing symmetry generators are important, especially for tracing the symmetries to the classical limit where all G s ⋆ reduce essentially to the commutative G s , as multiplicative operators on the functional space of classical observables.
We can writeG s =Ĝ L s −Ĝ R s , whereĜ L s ≡ G s (p, x)⋆ = G s (P L ,X L ) is a left action andĜ R s is the corresponding right action defined byĜ R s α ≡ α ⋆G s (p, x) = G s (P R ,X R )α. Analogously tô X L andP L coming from the left-invariant vector fields of the Heisenberg-Weyl group, there are those from the right-invariant ones given bŷ From Eq.(23) we see that In the above, for the wavefunctions, we show only the involved pair of variables in each case, and there is always no summation over indices. The other variables are simply not affected by the transformations. In terms of the parameters x µ and p µ , we have with the explicit action (no summation over the indices) where G ω µν are the infinitesimal SO(1, 3) transformation operators corresponding to the coset space action to be obtained from Eq.(5).
All the G −x µ , G p µ and G ω µν (and G θ = 1) make the full set of operators for the generatorŝ 3) group representing the symmetry on K, and constitute a Lie algebra within the algebra of physical observables.Ĝ R s set does the same as a right action, andĜ L s always commute withĜ R s ′ since, in general, [α L ,γ R ] = 0. These fourteen G s as multiplicative operators, of course, all commute among themselves. The commutators forG s are same as those forĜ L s , with however the vanishingG θ giving a vanishing [G p µ ,G −x ν ]. For any function α(p µ , x µ ), there are four associated operators onK. Those are α,α L ,α R andα, but only two of them are linearly independent. For our relativity symmetry operators, the independent set {G −x µ , G p µ , G ω µν ,G −x µ ,G p µ ,G ω µν } has the only non-vanishing commutators among them given by (we also have G θ = 1, the identity, andG θ = 0) Quantum dynamics is completely symplectic, whether described in the Schrödinger picture in terms of real/complex coordinates of the (projective) Hilbert space or the Heisenberg picture as a description in terms of the noncommutative coordinates [12]. The explicit dynamical equation of motion is to be seen as the transformations generated by a physical Hamiltonian characterized by an evolution parameter. In the H R (3) case of the usual ('nonrelativistic') quantum mechanics, it is G t = p i p i 2m . For our H R (1, 3) case, we consider G τ = pµp µ 2m with the parameter τ being the Einstein proper time, which is expected to give the standard covariant description of Einstein particle dynamics, as we see explicitly below.
For some s-dependent operator α(p µ (s), x µ (s))⋆ and a general Hamiltonian G s , Heisenberg equation of motion is given by The right-hand side of the equation is simply the Poisson bracket of α(p⋆, x⋆) and G s (p⋆, x⋆), functions of the noncommutative canonical variables p µ ⋆ and x µ ⋆. The equation can simply be written as and is exactly the differential version of the automorphism flow given in Eq.
which are exactly the standard form of Hamilton's equations of motion for the canonical η-Hermitian operator coordinate pairsX L µ -P L µ . As usual in a Hamiltonian formulation, the constant, or τindependent, momentumP L µ is obtained from the equations of motion as velocity multiplied by the particle mass m. Here, −m 2 is just the constant value of p ν p ν as 2mG τ .
For the Schrödinger picture, as η-unitary flows on K, we have the equation which for G τ ⋆ gives the τ -independent solution for φ in the exact form of the Klein-Gordon equation, provided that the G τ ⋆ eigenvalue is taken to be − m 2 . Explicitly, in terms of the basic variables p µ and x µ , we have giving the wavefunctions φ(p, x) = e i(2kµ−pµ)x µ for eigenvalues 2k µ kµ m . Eigenvalues of the momentum operators p µ ⋆ are 2k µ , satisfying (2k µ )(2k µ ) = −m 2 . The factor of 2 really corresponds to , as in the standard textbook expression. Finally, the τ -dependence is then given by d dτ φ = − m 4i φ, as expected.

B. Lorentz to Galilean Contraction
Contraction to Galilean limit has been presented in Sec. V at the kinematical level. In this section, we present the corresponding contraction in the observable algebra given in the WWGM formalism. Recall that the original Krein space under the contraction becomes reducible into a sum of essentially identical irreducible components, each being spanned by ) for a particular value of 'time' t (χ) . A general operator α(X L µ ,P L µ ) should then be seen as α(X L i ,P L i ,T L (χ) ,Ĥ L (χ) ) withX L i ≡X L (χ)i andP L i ≡P L (χ)i , from results of Eq. (32). Hence, on φ(p i , x i ) we have effectively Hermitian actions of operatorŝ , with the last two reduced to a simple multiplication by the 'variables' t (χ) and (formally infinite) e (χ) , respectively. All α(p µ ⋆, x µ ⋆) operators on φ(p i , x i ) reduce to α(p i ⋆, x i ⋆, t (χ) , e (χ) ), or rather simply to α(p i ⋆, We keep the ∂ e (χ) since theG β i may act on the mixed states. We have the newly relevant nonzero commutators involving a G β i , G t , or G −e , and aG s as well as those involving aG β i , G t , orG −e and a G s , all from the generators of the Lie algebra, as Since on the Hilbert space of the contracted theory we have only φ(p i , x i ) and the corresponding observable algebra as α(p i ⋆, x i ⋆), the loss of p 0 ⋆ and x 0 ⋆, the quantum observables of energy and time, means that the Heisenberg equation of motion, in the form of a differential equation in τ , effectively corresponds to the part of G τ ⋆ involving only p i ⋆. We where G t = p i p i 2m , giving the right time evolution in the 'non-relativistic', or H R (3), quantum theory, as expected. At the c → ∞ limit, the proper time is just the Newtonian time.
One can also see that the quantum Poisson bracket 1 2i [· · · , · · · ] does suggest that the now multiplicative operators t (χ) and e (χ) , from the original p 0 ⋆ and x 0 ⋆, are to be dropped from the canonical coordinates of the noncommutative symplectic geometry, in line with the Hilbert space picture. A G t of the form p i p i 2m + v(x i ), i.e. with a nontrivial interaction potential of course cannot be retrieved from a G τ which does not allow that, so long as the Einstein theory is concerned. If we allow a nontrivial v(x µ ) in G τ , however, everything works fine. For the latter G τ to be taken as a 'relativistic' Hamiltonian, one would have to allow violation of the Einstein relation of p µ p µ = −m 2 .

VII. CONTRACTION TO CLASSICAL THEORY IN BRIEF
In this section we look at the corresponding classical theory at the Lorentz covariant level through the contraction along the line of the one performed in the 'non-relativistic', H R (3), case presented in Ref. [1]. Only a sketch will be presented where the mathematics is essentially the same with the latter. The contraction trivializing the commutators between the position and momentum operators is obtained by rescaling the generators as and taking the limit k x , k p → ∞. The only important difference between k x and k p parameters is their physical dimensions, giving the X c µ and P c µ observables with their different classical units. For the corresponding operators we havê Our study is a part of our fully quantum relativity group-theoretically based program.
The constructed quantum mechanics is the 'relativistic' version of the so-called 'non-relativistic' theory based on the H R (3) group, or on theG (3) (3) they are both subgroups of the c → ∞ approximation of the H R (1, 3), obtained as a symmetry, or Lie algebra, contraction. The study here successfully completes the full picture from 'relativistic' quantum mechanics down to the 'relativistic' classical and the 'non-relativistic' quantum and classical theories as successive contractions/approximations. Some comments on the relation of the work to the noncommutative geometric perspective of quantum physics and quantum spacetime may be in order. Our idea of having the pseudo-unitary metric on the space of states comes mostly from the intuition on the need to take the pseudo-unitary Minkowski metric seriously as a quantum notion, for the noncommutative position and momentum operators as coordinates for the space [21]. The quantum phase space, exactly the projective Hilbert space for the 'non-relativistic' theory has been shown to serve as the quantum model of the physical space [1,14]. Well known as an infinite dimensional symplectic manifold, a noncommutative geometric picture of it has been presented [12] with the position and momentum operators as coordinates. A new conceptual notion of noncommutative values for the quantum observables [37] has been introduced to achieve a consistent interpretation of the values of the six coordinates being able to specify a point in the phase space otherwise described by an infinite number of real number coordinates. The noncommutative value of an observable for a state carries the full information the mathematical formulation of the theory actually contains for that. The notion gives an intuitive, but noncommutative, picture of quantum physics [13]; perhaps also a noncommutative notion of 'local realism'. The current Lorentz covariant theory is what would allow an analogous picture withX µ andP µ as coordinates bearing the Minkowski nature. Our background group theoretical framework has a stable symmetry with X-X and P -P type noncommutativity to which the H R (1, 3) symmetry is a contraction limit [16,17]. Other forms of covariant X-X noncommutative geometry with or without consideration of gravitation have been available in the literature(see for examples Refs. [38][39][40]). We share the belief that the proper theory of quantum gravity has to be a geometrodynamics of quantum, noncommutative, spacetime. Our framework has the unique feature that one has to take the 'phase space' as the model for the spacetime, as one irreducible representation.
The work focuses only on the formulation aspects, establishing such a theory that has all the nice properties mentioned and can successfully address the various concerns raised on such covariant theories. It may be considered mathematically involved, but unfortunately quite necessary. The only practical physical problem we have addressed is the free particle case, and arguably the covariant harmonic oscillator problem, with the latter being of great theoretical importance. However, practical application of 'relativistic' quantum mechanics is generally tricky [41]. Even with dynamics for the electron, quantum field theory is usually preferred. Applying the usual Dirac equation at the particle dynamics level has to confront the tricky issue of the negative energy solutions and the related zitterbewegung. There  [42,43]). On the other hand, our group theoretical framework naturally gives the dynamical theories in a covariant symplectic formalism with dynamical evolution to be described through an invariant parameter, like τ . The equation of motions, spin zero and probably also for the higher spin cases, are really in the form initiated by Stückelberg early the 1940's (see Refs. [44,45] and references therein). The studies stick with unitary representations of the Poincaré symmetry, though having to live with a non-Hermitian γ µP µ for a Dirac spinor. Yet there are some nice theoretical features from the theory, including good position operatorsX µ . Modifying the spin 1 2 theory to our pseudo-Hermitian setting has a good chance of improving its physics picture better. We plan on taking the theory more seriously along the line of our H R (1, 3) representation framework in terms the formulation before going into studies of practical systems. Hope to report on the results in the near future.
Appendix : Illustration of problems in unitary formulation of covariant harmonic oscillator In this appendix, we summarize the standard approach to covariant harmonic oscillator problem, which attempts to construct a unitary Fock space, assuming the position and momentum operators being Hermitian. A special attention will be given to its problems [18,20], here especially as seen in the φ(p µ , x µ ) wavefunction picture, which are all avoided in our pseudo-unitary representation. The difference in the two representations is in the inner product, which is simply for the unitary case. The first sign of the problem arises already in the abstract vector space.
The ladder operators are given essentially in the same way,â µ = η µν (X ν + iP ν ),â † µ = X µ − iP µ , where we drop the corresponding trivialη. As illustrated in the main text, the abstract algebraic analysis is not sensitive to the nature of the metricη. The same conclusion of m|n = (−1) n 0 δ mn cannot be avoided. [20]. So, the Fock space is still the same Krein space, which then cannot be the Hilbert space of the unitary representation of the H R (1, 3) symmetry. The only way to avoid that is to take a |0 state that is not Lorentz invariant [20], meaning that the Lorentz symmetry is spontaneously broken in the system, which hardly sounds like the quantum version of the classical covariant harmonic oscillator system or anything we may have a good reason to be interested in. The key thing is that the noncompact nature of SO (1,3) gives no finite-dimensional unitary representation.
Since the Hamiltonian for the problem, or the operatorN, is Lorentz invariant, the nlevel subspaces are likewise invariant and hence can only be infinite-dimensional, so long as unitary representations are concerned. The states on a fixed n-level do not transform as symmetric Minkowski n-tensors. That is the key issue behind the incompatibility of the latter and the kind of nice physics picture one would like to have for the system [21], which our pseudo-unitary formulation successfully retrieved. x µ x µ + p µ p µ − ∂ 2 ∂p µ ∂p µ − ∂ 2 ∂x µ ∂x µ + 2ix µ ∂ p µ − 2ip µ ∂ x µ − 2 operator. One can easily check that φ n (p µ , x µ ) = e − xµx µ +pµp µ 2 3 µ=0 (x µ − ip µ ) nµ are solutions for the eigenvalue n 0 + n 1 + n 2 + n 3 , with φ o corresponding to |0 state. To stick to the probability interpretation with the trivial measure in the integral inner product, one has to restrict the domain of the wavefunctions to spacelike region of p µ and x µ .
In order to normalize the wavefunction, we need to calculate the integral d 4 xd 4 p π 4φ n (p µ , x µ ) φ n (p µ , x µ ) = to deal with. The integral can be evaluated with coordinates in a polar form as in Ref. [18], by defining r 2 = x µ x µ , ρ = x 0 / √ x i x i , and the spatial angular coordinate of which we skip the details. We obtain for the n = 0 case, with Ω denoting the spatial solid angle. The ρ-integral is still divergent. The integrand for the specific case is ρ-independent. The divergence is simply due to the infinite range of the boost parameter (ρ being its hyperbolic tangent). Hence, it has been suggested to define the integral with 'the infinite volume factor' absorbed [18,20].
However, there is really no sensible way to do that so long as the Fock state wavefunctions are concerned. The corresponding ρ-integral of φ n for a nonzero n 0 has an extra factor of ρ 2n 0 (1−ρ 2 ) n from the integrand, giving a higher order divergence for each larger n value. From the structure of the full integral inner product , it is clear that such contributing terms of the higher order divergence stay. That is to say, none of Fock state wavefunctions are really normalizable under the unitary formulation. This is not an artefact of the coherent state framework. The usual Schrödinger wavefunctions ψ(x µ ) for the Fock states have the same problem.
There is an alternative approach of taking a timelike, instead of spacelike, parameter restriction for the domain of the wavefunctions, really corresponding to defining |0 as satisfying a † |0 = 0. Similar problems persist.