Arxiv:2009.09441V1

Home , Translational symmetry

arXiv:2009.09441v1 [cond-mat.quant-gas] 20 Sep 2020 r ieaBaasltie h o-nryeetv field effective low-energy the lattice, Bravais a symme- lat- like space-translational optical discrete try or the materials other have crystalline which the of tices study On the Poincaré group in representations. the hand, be irreducible that can ele- fact different helicity of the has or from existence spin deduced The different directly with Poincaré group. particles a the mentary form - which group spacetime transformations Lie arbitrary Lorentz under and invariant translations keep theories field nrcniigqatmmcaiswt oet symme- Lorentz efforts with of try mechanics result quantum the reconciling indeed in is which theory field quantum Lorentz the has that model symmetry. whether lattice to quantum as a arises exists naturally there question condensed A in ignored physics. matter been always has symmetry Lorentz intr prtr nteHletspace Hilbert the an- in or operators unitary tiunitary the by transformations symmetry physical o nesadn h oooia insulators key topological the is the which understanding degeneracy Kramers’ for the topological to leads Floquet years metry recent the in matter classify of to the phases employed and respectively, was quasi-energy, replaced latter and be state energy and Floquet eigenstate discrete by the The sys- that requires periodically-driven tems particles. in of symmetry spectrum time-translational the constrain tries ls h Bose-Hubbard the ples, h otnoseeg fafe atceboe noa into broken particle free bands a Bloch of causes of series symmetry energy the- space-translational continuous lattice discrete a the of The construction all the ory. symmetries in ro- reversal roles translation, time important play spacetime and theories inversion The lattice space tation, quantum properties. new exotic with explore continuously to rcl tm notcllattices optical in atoms tracold h oet ymtyi setal motn in important essentially is symmetry Lorentz The h elzto fitiun atc oes o exam- for models, lattice intriguing of realization The 7 inrstermpvdtewyfrrepresenting for way the paved theorem Wigner’s . hoy h pctm onsweeteGensfnto sn is si a function on Green’s localized stay the particles symmetry. where the propagation, points calcu in the We spacetime Hamiltonians During increases. typical The sites Some between theory. distance operators. the unit field as the the fluctuates both of expressing terms by in a particles the which indistinguishable build then representation, of We of quasi-energy. and existence quasi-momentum the transformationsbetween for Lorentz unitary conditions by discrete represented the The are d which the Poincaré with group, lattice. discrete coexists Bravais that a symmetry recentl Lorentz on the unique the review We is which lattice. Bravais one-dimensional the nti ae,w eeo h unu hoyo atce tha particles of theory quantum the develop we paper, this In .INTRODUCTION I. 4 1 h oainadprt symme- parity and rotation The . eateto hsc,Zein omlUiest,Jnu 3 Jinhua University, Normal Zhejiang Physics, of Department 1 n Haldane and ona´ rsa nteoedmninllattice one-dimensional the on Poincaré crystal 5 3 h iervra sym- reversal time The . norgsresearchers encourages 2 8 oes ihul- with models, h quantum The . 6 Dtd etme 2 2020) 22, September (Dated: u the But . e Wang Pei oe u w ie toiia distance any original neighbors” at sites ”nearest two si- not But of these are relativity of more. one, they the frame first that to so due second the multaneity phases to the different have relative In sites moving two is phase. that same reference are the that in sites nearest-neighbor oscillating two between distance tial rpryo h igepril ibr pc.Orap- Our space. Hilbert single-particle finite-dimensional the the of utilize property and Poincaré group crete vrg seRf 0frarcn eiw.Frstrictly- symmetry For Lorentz review). recently discrete found recent the was a to symmetries, for approach 10 preserved set Ref. on invariance causal (see Lorentz average the preserves in which gravity, introduced quantum was space- process the lattice Poisson a with a following lattice sprinkled random points a time employing by is lattice lattice. be- survive Bravais not invariance, a does Lorentz on symmetry the Lorentz continuous to the conform cause not does theory eidcly noefaeo reference, of shaken frame is lattice one the if In possible is this periodically. view, clas- of From point contraction. sical length of spite in transformations osrc h unu hoista aetediscrete the length lattice have of Bravais one-dimensional that a theories on Poincaré quantum symmetry the construct Poincaré symmetry. say discrete the we has Poincaré group, theory coordinate discrete the arbitrary that the under a in form is same transformation which the physical a Poincaré keeps group If theory Poincaré discrete lattice group. continuous the the spacetime of - the subgroup group on language a translations together the form discrete transformations In Lorentz the discrete lattice. with the spacetime group, a lattice of Bravais into the extended that requires be direction reference time second the the along in becomes distance their frame while now, phase pc-rnltoa ymty rvi atc fcon- of lattice Bravais stant A symmetry. space-translational 1, n a frsretn h oet ymtyo a on symmetry Lorentz the resurrecting of way One nti ae,w eeo eea prahto approach general a develop we paper, this In ∗ geo e ie opeev h Lorentz the preserve to sites few a or ngle r prtr n lqe Hamiltonians Floquet and operators ary a icvrddsrt oet symmetry, Lorentz discrete discovered y prtr naqatmter.W find We theory. quantum a in operators eepesda h ogunerelation congruence the as expressed re aeteGensfntoso h lattice the of functions Green’s the late ep nain ne ylcgopo Lorentz of group cyclic a under invariant keeps n pctm rnltosfr the form translations spacetime and N srt pc rnltoa symmetry translational space iscrete ld h ogrnehpigwhich hopping long-range the clude oet-nain aybd theory many-body Lorentz-invariant a iceePicresmer on Poincaré symmetry discrete has t neodslyaltiestructure. lattice a display onzero esatfo h ersnain fdis- of representations the from start We . a u ocnrcin h periodicity The contraction. to due 11,12 10,China 21004, hc oxsswt h discrete the with coexists which , na aetesame the have a stespa- the is 9 Such . 2 proach is exhaustive. The Lorentz invariance is trans- The discrete Lorentz group is a subgroup of the con- lated into a self-consistent equation of quasi-energy and tinuous Lorentz group. Its elements are quasi-momentum, which is further changed into a congruence relation modulo N between integer variables. This γ(j) sign(j) γ(j)2 1/c L(j)= − , congruence equation has solutions for some but not all sign(j) γ(j)2 1 cp γ(j) − N. As examples, we enumerate the N (N < 100) for p (2) which the solutions exist, and study the corresponding where dispersion relations of a single particle. We then turn 1 j 1 j into the many-body theory of indistinguishable bosons γ(j)= γ(1) + γ(1)2 1 + γ(1) γ(1)2 1 or fermions, obtaining the expression of discrete Lorentz 2 − 2 − − p p (3) transformations and spacetime translations in terms of the field operators. It is appropriate to call such a many- and j = 0, 1, 2, is an integer. γ(1) > 1 is an in- body model a Poincarécrystal. For the convenience teger or half-integer,± ± ··· dubbed the generator. A discrete of designing an experiment, we construct the Floquet Lorentz group is uniquely determined by γ(1). Interest- Hamiltonians of the Poincarécrystals, which might pos- ingly, γ(j) at arbitrary j can be obtained by using the sibly be realized with ultracold atoms in optical lattices. iterative relation γ(j +1)=2γ(1)γ(j) γ(j 1), given The properties of a Poincarécrystal are studied. We cal- γ(0) = 1 and γ(1). In this paper, we consider− − an integer culate the Green’s functions, which exactly reflect the generator (γ(1) = 2, 3, ), hence γ(j) is an integer. We underlying Lorentz symmetry. ··· define ζ(j) = sign(j) γ(j)2 1/ γ(1)2 1, which sat- The paper is organized as follows. In Sec. II, we re- isfies a similar relation ζ(j +1)− = 2γ(1)ζ−(j) ζ(j 1) view the definitions of discrete Lorentz transformations p p with ζ(0) = 0 and ζ(1) = 1. One can prove− that −ζ(j) and discrete Poincarégroup. Sec. III introduces the rep- is an integer. Some useful formulas about γ(j) and ζ(j) resentations of discrete Poincarégroup which determine are listed in Appendix A. Now we reexpress the discrete how a single-particle quantum state transforms under Lorentz transformation as the unitary operators in the group. In Sec. IV, we construct the orthonormal basis of the single-particle Hilbert γ(j) ζ(j) γ(1)2 1/c L(j)= − . (4) space. In this process, the conditions for quasi-energy ζ(j) γ(1)2 1 cp γ(j) and quasi-momentum are found. The many-body the- − p ory of Poincarécrystal is developed in Sec. V in which An important property of L(j) is that it is equal to L(1) the field operators are defined, the symmetry transfor- to the power j. As a consequence, we have L(j + j′) = mations are expressed in terms of them, and the Flo- L(j)L(j′). The discrete Lorentz group is a cyclic group. quet Hamiltonians are presented. Sec. VI discusses the A general coordinate transformation is a combina- Green’s functions of a Poincarécrystal. Finally, Sec. VII tion of Lorentz transformation and spacetime translation, summarizes our results. which is expressed as

t′ t ∆t II. DISCRETE POINCARE´ GROUP = L + , (5) x′ x ∆x

The discrete Poincarégroup is a group of coordinate where (∆t, ∆x)T is the translation vector. We use P transformations introduced in Ref. [11] and [12]. In this to denote such a transformation. P maps (t, x)T into section, we review its definition and properties. T (t , x ) . This fact can be used to derive the rule of mul- We consider the 1+1-dimensional spacetime, and use ′ ′ tiplication for P . All the P s form a continuous group xµ with µ = 0, 1 to denote the time and space coordi- - the so-called Poincarégroup, which is well known in nates, respectively. Sometimes, we use t and x instead special relativity and quantum field theory. of x0 and x1. The spacetime metric η is diagonal with µν The discrete Poincarégroup, denoted by , is a sub- elements η00 = 1 and η11 = 1. A Lorentz transforma- P − group of the continuous Poincarégroup. includes only tion is a rotation in the spacetime. It is represented by a P T discrete Lorentz transformations and discrete spacetime 2-by-2 matrix that acts on the vector (t, x) , reading translations. We use a to denote the period of space γ sign(v) γ2 1/c translation, and define T = γ(1)2 1 a/c to be the L = − , (1) − sign(v) γ2 1 cp γ period of time translation. Thep primitive vectors e = − T f T p (T, 0) and = (0,a) generate a 1+1-dimensional rect- where γ = 1/ 1 v2/c2 1 is dimensionless and v angular lattice, dubbed (see fig. 1(a)). The elements of − ≥ Y is the velocity ofp one reference frame with respect to the are written as P (j,m,n) with j,m,n = 0, 1, 2, other. c is the speed of light in special theory of relativity. beingP integers. P (j,m,n) denotes a Lorentz transforma-± ± ··· In this paper, we treat c as a speed of constant used to tion L(j) followed by a translation of vector me+nf. For define the Lorentz transformation, but the value of c can examples, P (j, 0, 0) is indeed L(j), and P (0, m, 0) and be chosen arbitrarily. The Lorentz transformations with P (0, 0,n) are pure temporal and spatial translations, re- c

t (a) keep some Lorentz symmetry on a lattice model, its symmetry group must be a discrete Poincar´egroup.

t (c) III. REPRESENTATION OF DISCRETE POINCARE´ GROUP O x

In quantum mechanics, a physical state is a vector in the Hilbert space . A symmetry transformation is a unitary operator mappingH to itself. To construct N-2 5 H N-1 4 a quantum theory with discrete Poincarésymmetry is 0 2 3 x E 1 equivalent to find a unitary representation of the group (b) . We set ~ = 1 throughout the paper. P If P (j,m,n) is the transformation from the co- ∈ P ordinate system of observer K to that of observer K′, we use Uˆ(j,m,n) to denote the unitary transformation mapping the quantum states seen by K to those seen O k by K′. Uˆ(j,m,n) obeys the multiplication laws as same as Eq. (6). To save time, we do not write them down again. Let us consider the classification of single-particle states according to their transformation under Uˆ(j,m,n). Since Uˆ(0,m,n) of different (m,n) all commute with each other, they have common eigenstate, denoted as p . Ac- | i cording to Uˆ(0,m,n)= Uˆ m(0, 1, 0)Uˆ n(0, 0, 1), the eigen- FIG. 1. The schematic diagram of (a) a spacetime Bravais ˆ lattice Y, (b) the corresponding reciprocal lattice R and (c) a value of U(0,m,n) for arbitrary (m,n) is determined by lattice of finite size with periodic boundary conditions along the eigenvalues of Uˆ(0, 1, 0) and Uˆ(0, 0, 1). The eigen- the spatial direction. value of a unitary operator must be a complex number of unit magnitude. Without loss of generality, we assume the eigenvalues of Uˆ(0, 1, 0) and Uˆ(0, 0, 1) to be µ ν µ ν under L(j). This guarantees that P (j,m,n) is a group. iηµν e p iηµν f p µ µ { } e− and e− , respectively, where e and f By using the definition of P (j,m,n) and the properties are the components of the translation vectors e and f, of L(j) and by studying the transformation of me + nf respectively, and the summation convention has been ap- under L(j), we obtain the group’s multiplication laws, plied. This assumption makes it clear that p should be which are listed below: a 1+1-dimensional vector with p0 and p1 being its com- P (j,m,n)= P (0,m,n)P (j, 0, 0), ponents. Sometimes, we also use E = p0 and k = p1 to denote the components of p. Now we have P (j, 0, 0) = P j (1, 0, 0), m n µ ν iηµν (me+nf) p P (0,m,n)= P (0, 1, 0)P (0, 0, 1), Uˆ(0,m,n) p =e− p | i | i (7) [P (0, 1, 0), P (0, 0, 1)]= 0, (6) inak+imTE =e− p . 1 | i P (j, 0, 0)P (0,m,n)P − (j, 0, 0) = P (0,m′,n′) In Eq. (7), (me + nf) is recognized as the translation m′ γ(j) ζ(j) m with = 2 , vector of P (0,m,n). Taking a comparison with quantum n′ ζ(j) γ(1) 1 γ(j) n − field theory, one may naively identify E as the ”energy” 1 where P − denotes the inverse of P and [P, P ′] denotes and k as the ”momentum”. We emphasize that the en- the commutator of P and P ′. It is easy to verify that ergy or momentum may not be well-defined if there does m′ and n′ are integers given an integer pair (m,n), since not exist a continuous spacetime-translational symmetry. ζ(j) and γ(j) are both integers. Eq. (6) can be viewed In the presence of discrete translational symmetry, it is as the defining properties of . more appropriate to call E the ”quasi-energy” and k the A few comments are necessary.P (i) P (1, 0, 0), P (0, 1, 0) ”quasi-momentum”. and P (0, 0, 1) are the three generators of . The other To decide how p transforms under Uˆ(j, 0, 0), we uti- | i elements of can be expressed as the productP and power lize the last multiplication law in Eq. (6) by rewriting it of them. (ii)P is uniquely determined by γ(1), c and a. as Two discreteP Poincarégroups with the same γ(1) are iso- ˆ ˆ ˆ ˆ morphic to each other. (iii) The discrete Poincarégroups U(j, 0, 0)U(0,m,n)= U(0,m′,n′)U(j, 0, 0). (8) are the only groups that contain at least one non-identity Lorentz transformation at the same time whose spatial Acting both sides on p , we find that Uˆ(j, 0, 0) p | i | i translations are discrete. In other words, if we want to is an eigenvector of Uˆ(0,m′,n′) with the eigenvalue 4

µ ν iηµν (me+nf) p e− . Without considering the internal de- IV. SINGLE-PARTICLE HILBERT SPACE grees of freedom (spin or helicity), we can simply assume A. How to construct the Hilbert space Uˆ(j, 0, 0) p = p′ , (9) | i | i The eigenvectors of a unitary operator with differ- where p′ satisfies ent eigenvalues are orthogonal to each other. p is the µ ν ′ ′ µ ′ν | i iηµν (me+nf) p iηµν (m e+n f) p eigenvector of Uˆ(0,m,n). For p,p′ B , we then have e− = e− . (10) ∈ f p p′ = δ ′ . Suppose the dimension of the single- h | i p,p Before giving the solution of Eq. (10), we define the particle Hilbert space is finite, denoted as N. We then reciprocal lattice of , denoted by whose two prim- choose N different p with p Bf to form an orthonor- Y R µ ν | i ∈ itive vectors (denoted by g and h) satisfy ηµν e g = mal basis of the single-particle Hilbert space. η f µhν = 2π and η eµhν = η f µgν =− 0. g and h There exist some conditions for the set p being a µν µν µν {| i} are found to be g = (2π/T, 0)T and h = (0, 2π/a)T, basis. To obtain these conditions, we study the single- respectively. Obviously, is a rectangular lattice (see particle basis in real spacetime. Since the single-particle fig. 1(b)). And for arbitraryR lattice vectors r and Hilbert space has N dimensions, our model must con- µ ν ∈ R y , ηµν r y must be integer times of 2π. is a lattice tain N spatially-different sites on which the particle can located∈ Y on the t-x plane (real spacetime), whileY is a be located. Without loss of generality, we number these lattice located on the E-k plane. R sites from 0 to N 1 whose positions are correspond- − p = E, k transforms into p′ = E′, k′ under ingly x =0, a, 2a, , (N 1)a. To maintain the space- | i | i | i | i ··· − the unitary transformation Uˆ(j, 0, 0). The solution of translational symmetry, we adopt the periodic boundary Eq. (10) can be expressed as (see Appendix B for the condition along the x-axis, i.e., acting P (0, 0, 1) on the derivation) site at x = (N 1)a results in the site at x = 0. In other words, we− treat two spacetime sites (mT,na) and E′ E (mT,n′a) in fig. 1(a) as being equivalent to each other if =R L(j) , (11) k′ k n n′ (mod N). The quotient set of with respect to this≡ equivalence relation is topologicallyY equivalent to a where =R denotes an equivalence relation. Two vectors cylinder (see fig. 1(c)), which contains the sites (mT,na) on the E-k plane are equivalent to each other if and only with

B. The case γ(1) = 2

A few examples are helpful for explaining the above ideas. Let us consider the discrete Poincar´egroup with γ(1) = 2 (the smallest integer γ(1)). P The Lorentz matrix depends on γ(j) and ζ(j). Table I lists the ﬁrst few elements of the integer sequences γ(j) and ζ(j) which increase exponentially with j. Indeed, the ratios of γ(j +1) to γ(j) or ζ(j +1) to ζ(j) both converge to γ(1)+ γ(1)2 1 3.732 in the limit j . While − ≈ → ∞ the ratio ofpγ(j) to ζ(j) converges to γ(1)2 1 1.732. − ≈ On the other hand, the periods of spacep and time translations are a and T = √3a/c, respectively. The Lorentz matrix L(j) transforms the spacetime point (mT,na) into (m′T,n′a) with

j m 2 1 m ′ = . (18) n′ 3 2 n We consider a particle on a 1D lattice of N sites under periodic boundary conditions. The quasi-momentum of the particle can only take N possible values which are k = 2πk/˜ (aN) with k˜ =0, 1, ,N 1. The quasi-energy is ˜ ···˜ − Ek =2πEk˜/(TN), where Ek˜ is an integer between 0 and N 1. The unitary operator Uˆ(1, 0, 0) transforms the − FIG. 2. The six Lorentz-invariant dispersion relations at N = single-particle state E , k into E′ ′ , k′ with | k i | k i 18.

E˜′ 2E˜ +3k˜ (mod N) , ≡ (19) k˜′ E˜ +2k˜ (mod N) . ≡ (a) (b) ˜ The dispersion relation Ek˜ must keep invariant under the transformation (19). Let us ﬁrst consider the smallest nontrivial dimension N = 2. There exist two dispersion relations which are ˜ ˜ Lorentz-invariant. We can choose either Ek˜ =0 at k =0 ˜ ˜ ˜ ˜ and Ek˜ =1 at k = 1, or E0 = 1 and E1 = 0. The other two relations (i.e. E˜0 = E˜1 = 1 and E˜0 = E˜1 = 0) do not (c) (d) keep invariant under the transformation (19), therefore, they cannot be the dispersion relations of a theory with discrete Lorentz symmetry. For N = 3, the unique Lorentz-invariant relation is E˜0 = E˜1 = E˜2 = 0, i.e., the quasi-energy is zero at arbitrary quasi-momentum (a ﬂat band). But for N =4 or N = 5, there exists no Lorentz-invariant Ek, hence, no quantum theory can be built. Table II lists all the (e) possible N (N 100) for which the Lorentz-invariant (f) ≤ Ek exists.

j 0 1 2 3 4 5 6 · · ·

γ(j) 1 2 7 26 97 362 1351 · · ·

· · · ζ(j) 0 1 4 15 56 209 780 FIG. 3. The sawtooth Lorentz-invariant dispersion relations at (a) N = 13, (b) N = 23, (c) N = 39, (d) N = 47, (e) TABLE I. The ﬁrst few elements of the sequences γ(j) and N = 59 and (f) N = 73. ζ(j) as γ(1) = 2. 7

For most of N in table II, the Lorentz-invariant dis- results are independent of whether the commutative or persion relations are not unique. Fig. 2 plots the six anticommutative relations stand. Therefore, we will not different Ek for N = 18. They are symmetric with re- distinguish whether the particle is a fermion or boson any spect to k = π/a. By using Eq. (19) one can prove that, more. ˆ if Ek is Lorentz-invariant, then both Ek′ = E2π/a k and The unitary transformation U(j,m,n) does not rotate − Ek′′ =2π/T Ek are Lorentz-invariant (see Appendix C the vacuum state, i.e. Uˆ(j,m,n) vac = vac . This is − | i | i for the proof). Ek′ and Ek′′ are the mirror images of Ek guaranteed if the right-most field operator of Uˆ(j,m,n) with respect to k = π/a and E = π/T , respectively. In is an annihilation operator. Acting Uˆ(0,m,n) on both fig. 2, the six functions are indeed the mirror images of sides of Eq. (20) and then using Eq. (7), we obtain each other with respect to E = π/T . Fig. 3 plots the dispersion relations that have a saw- ˆ ˆ 1 inak+imTEk U(0,m,n)âk† U − (0,m,n)= e− aˆk† . (21) tooth shape for different N. Such a sawtooth Ek is found for many N (e.g. N = 11, 13, 23, 39, 47, 59, Eq. (21) tells us how the field operator transforms under 73, 83, ). The sawtooth Ek is not symmetric with a coordinate translation. For those who are familiar with ··· respect to k = π/a or E = π/T , but it changes into quantum field theory, it is easy to find itself after two successive reflections to k = π/a and ina P kaˆ† aˆ +imT P E aˆ† aˆ E = π/T . Therefore, they always appear in pairs (Ek Uˆ(0,m,n)= e− k k k k k k k . (22) and Ek′ = E2π/a k = 2π/T Ek are both Lorentz- − − invariant). Furthermore, the sawtooth Ek all starts from Taking m = 0 and n = 1, we find the space-translational ia P kaˆ† aˆ E0 = 0, as far as we know. operator to be Uˆ(0, 0, 1) = e− k k k . It is clear that

k kaˆk† aˆk must be the quasi-momentum operator of the quantumP field. V. MANY-BODY THEORY Next we show how to express the unitary operator Uˆ(j, 0, 0) that corresponds to a Lorentz boost. Since A. The field operators and unitary transformations Uˆ(j, 0, 0) = Uˆ j (1, 0, 0), we focus on the expression of Uˆ(1, 0, 0). By using the field operators, Eq. (9) can be The generalization from single-particle Hilbert space reexpressed as to many-body Hilbert space is straightforward. In the ˆ ˆ 1 many-body space, it is convenient to express the unitary U(1, 0, 0)âk† U − (1, 0, 0)=a ˆk† ′ , (23) transformations Uˆ(j,m,n) in terms of the field operators. where the relation between k′ and k has been given We usea ˆ† to denote the field operator that creates one k in Eq. (11) with j = 1. The transformation from particle of quasi-momentum k. The single-particle state ˜ ˜ can then be expressed as k = 2πk/(aN) to k′ = 2πk′/(aN) is also equivalently defined in Eq. (16). k k′ is indeed a permutation of 2π 1 2π 2 → 2π N 1 the set 0, , , , − . As is well known, a Ek, k =â† vac (20) a N a N ··· a N | i k | i permutation of N elements can be represented by a N- with vac denoting the vacuum. The vacuum is a state by-N matrix (so-called permutation matrix), denoted as | i σ, whose entry σk′,k is 1 if k maps to k′ but 0 otherwise. without any particle, which is defined bya ˆk vac = 0 for | i σ is a unitary matrix, therefore, it can always be writ- arbitrary k. Here the particle can be either a fermion or a iK ten as σ = e− with K being a hermitian matrix. The boson. If it is a fermion,a ˆ† satisfies the anticommutative k Lorentz boost can then be expressed as relation, i.e. aˆk, aˆk† ′ =a ˆkaˆk† ′ +âk† ′ aˆk = δk,k′ . If the + † h i i P ′ aˆ ′ Kk′,kaˆk Uˆ(1, 0, 0) = e− k,k k . (24) particle is a boson,a ˆk† satisfies the commutative relation, i.e. aˆ , aˆ† ′ =a ˆ aˆ† ′ aˆ† ′ aˆ = δ ′ . The following k k k k − k k k,k h i B. The Floquet Hamiltonians

2 3 6 11 13 18 22 Let us consider the time evolution operators. Taking m = 1 and n = 0 in Eq. (22), we ﬁnd Uˆ(0, 1, 0) = iT P− E aˆ† aˆ − 23 26 30 33 37 39 46 e− k k k k . Uˆ(0, 1, 0) is the time-translational transformation, mapping− the state seen by an observer 47 54 59 61 66 69 71 K to that seen by K′ (notice that the Heisenberg picture is automatically chosen as we talk about the transforma- 73 74 78 83 90 94 97 tion between observers). When the clock of K reads t =0 (the reference time in the Heisenberg picture), the clock TABLE II. The possible length (N) of a one-dimensional Bra- of K′ reads t′ = t T = T . Therefore, if we turn to the vais lattice that can host a Lorentz-invariant quantum theory Schr¨odinger picture,− Uˆ(0−, 1, 0) must be the evolution with γ(1) = 2. operator over one period T−. In the periodically-driven

2 ¢£ 0 ¡

-1

-2 -3

0.5

-0.5 -3 -1

-2.5 -0.5

FIG. 4. The six coupling functions Jn of the corresponding FIG. 5. The absolute coupling functions of a sawtooth disper- Lorentz-invariant theories at N = 18. sion relation at N = 13, 23, 39, 47, 59 and 73, respectively. problems, people often formally write the evolution oper- at site n. We then define the creation operator at site n iT HˆF ator as e− where HˆF is the so-called Floquet Hamil- according to tonian. Similarly, we introduce the Floquet Hamiltonian ˆ ˆ iT HF ˆ 0,n = b† vac (26) in our theory. By comparing e− with U(0, 1, 0), | i n | i we find − ˆ where n = 0, 1, ,N 1. And bn† is either a fermionic ˆ † ··· − HF = Ekaˆkaˆk. (25) ˆ ˆ ′ operator satisfying bn, bn† ′ = δn,n or a bosonic op- Xk h i+ ˆ ˆ ˆ erator satisfying bn, bn† ′ = δn,n′ . By using Eq. (14) Notice that HF is not defined in a unique way. One does h i not change the evolution operator Uˆ(0, 1, 0) by adding and (20), we connect ˆb† to the field operators in quasi- − n integer times of 2π/T to Ek. In this paper, we choose momentum space: 0 Ek < 2π/T to be in the Brillouin zone when defining ≤ 1 inak HˆF . ˆb† = e− aˆ† . (27) n √ k We emphasize that time is continuous. Only the time- N Xk translational symmetry is discrete. We usually choose 1 inak t = 0 to be the reference time, hence, we can obtain The inverse transformation isa ˆ† = e ˆb . We k √N n n† the quantum states in the Schrödinger picture at t = reexpress the Floquet Hamiltonian (25)P as 0, T, 2T, by acting Uˆ(0, m, 0) on the state at t = 0. ± ± ··· The quantum states not at these times are less interesting ˆ ˆ ˆ ˆ ˆ HF = J0 bn† bn + Jn n′ bn† bn′ + h.c. , (28) to us. They are not unambiguously defined in the theory. − Xn n>nX′ For example, one cannot simply claim that the state at iHˆF T/3 t = T/3 is e− acting on the state at t = 0. Indeed, where the first and second terms are the onsite poten- the trajectory of the state from t = 0 to t = T can be tial and hopping between different sites, respectively. Jn freely chosen, which is not constrained by the discrete with n being an integer is the coupling function, defined 1 iakn Poincarésymmetry. as Jn = N k Eke . From the definition we immedi- Sometimes, it is more convenient for designing an ex- ately obtainPJn = JN∗ n where J ∗ denotes the complex − periment setup if we express the Floquet Hamiltonian conjugate of J. Since the quasi-energy Ek can be in- in terms of the field operators in real space. To define creased or decreased by arbitrary integer times of 2π/T , the real-space field operators, we notice that m,n with the coupling function J is also not uniquely defined. | i n m = 0 represents the quantum state of a particle located Next we study Jn as Ek is in the Brillouin zone Bf . 9

The coupling function is determined by the dispersion Eq. (13), we immediately derive relation, and then depends on N. As a special example, ˆ ˆ ˆ ˆ ˆ ˆ 1 we obtain HF = 0 as N = 3 (Ek = 0). HF = 0 means bn† +n′ (m + m′)= U(0,m′,n′) bn† (m) U − (0,m′,n′). that the quantum state repeats itself after one period of (29) evolution. Fig. 4 plots J at n = 0 for the six different E (see In Eq. (29), if n + n is not in the range [0,N), we need n 6 k ′ fig. 2) as N = 18, in which case Ek = E2π/a k di- to replace it by an integer congruent to n + n′ modulo − rectly indicates that Jn is a real number. Interestingly, N, because we choose the periodic boundary conditions Jn vanishes for most of n. For example, in the top left along the x-axis. Recalling Uˆ(0,N, 0) = 1, we can also panel, only J = J = 4π/(NT ) and J = π/(NT ) choose the periodic boundary conditions along the t-axis. 3 15 − 9 are nonzero. This means that the particle can only hop We then replace m + m′ by an integer congruent to it between two sites at a distance of 3 or 9 from each modulo N if m + m′ is not in the range [0,N). In next, other. Notice that the distance between two sites be- ˆ only the field operators bn† (m) with 0 n,m < N are ing 15 is equivalent to it being 3 due to the periodic considered. ≤ boundary condition. In the bottom right panel, we see There are a few ways to derive the transformation J3 = J15 =4π/(NT ) but J9 = π/(NT ). Indeed, Ek in ˆ − of bn† (m) under a Lorentz boost. One can start from the bottom right panel is the mirror image of that in the Eq. (23) and notice that ˆb is the Fourier transformation top left panel with respect to E =2π/T (see fig. 2). It is n† ofa ˆ† . By using the relation (10), we obtain easy to prove that, if E′ =2π/T E , then J ′ = J . k k − k n − n Fig. 5 plots the absolute value of the coupling functions 1 Uˆ(j, 0, 0)ˆb (m)Uˆ (j, 0, 0) = ˆb† ′ (m ), (30) that correspond to the sawtooth dispersion relations in n† − n ′ fig. 3. Within the different N in fig. 5, N = 39 is spe- where (m,n) (m′,n′) has been defined in Eq. (17). cial, for which Jn is zero at two thirds of n. But Jn is → ˆ nonzero for the other N (N = 13, 23, 47, 59, 73), for The set of operators bn† ′ (m′) transforms under n o which the dispersion relation satisfies Ek = 0 at k = 0 Uˆ(j,m,n) in the same way as the spacetime lattice and E = 2π/T E . From this property6 of6 E k 2π/a k k (m′T,n′a) equipped with a congruence relation (mod- we immediately derive− Re−J = π/(NT ) where ReJ { } n − n ulus N) transforms under P (j,m,n). denotes the real part of Jn. The absolute value Jn displays a peak structure. Especially for N = 47| and| N = 73, we clearly see a few evenly-spaced sharp peaks. B. The Green’s functions The coupling functions we find here are significantly different from those in the other quantum many-body models. Not only the long-range hopping exists, but also the The Green’s functions are the correlations between the coupling strength fluctuates with the distance. field operators at two spacetime points. Let us consider the retarded Green’s function

r ˆ ˆ G (m1,n1,m2,n2)= iθ(m1 m2) bn1 (m1), bn† 2 (m2) , − − hh i i VI. GREEN’S FUNCTIONS (31)∓ where the upper (lower) sign corresponds to bosons In above, we have introduced the general approach (fermions), respectively, the Heaviside function θ(∆m) to construct a quantum many-body theory with discrete is 1 if ∆m 0 but 0 otherwise, and the brackets mean ≥ r Poincarésymmetry on a 1D lattice of length N. The Flo- the expectation with respect to some quantum state. G quet Hamiltonians for some specific N have been shown is the correlation between the field at position x1 = n1a as examples. Next, we discuss the properties of such a at time t1 = m1T and that at position x2 = n2a at ˆ theory. Especially, we focus on the Green’s functions in time t2 = m2T . In this paper, HF has a quadratic form, ˆ ˆ real spacetime. therefore, bn1 (m1), bn† 2 (m2) is a number instead of an h i operator, and Gr is then independent∓ of the quantum state. A. Field operators in real spacetime The discrete Poincarésymmetry manifests itself in the retarded Green’s function. The (anti)commutator transforms in the same way as the field operators, while the It is useful to show how the real-space field operators latter transforms in the same way as the spacetime coor- transform under a symmetry transformation Uˆ(j,m,n). dinates. We then have ˆ In the Heisenberg picture, the operator bn† (t) at t = ˆ ˆ ˆ ˆ ˆ ˆ 1 mT (bn† (m) for short) is naturally defined as bn† (m) = U(j,m,n) bn1 (m1), bn† 2 (m2) U − (j,m,n) ˆ ˆ imT HF ˆ imT HF h i∓ e bn† e− . And its relation with the single- (32) ˆ ′ ˆ ˆ = bn1 (m1′ ), b† ′ (m2′ ) , particle state m,n is bn† (m) vac = m,n . From n2 | i | i | i h i∓ 10

FIG. 7. (Color online) The left panel plots the dispersion

r relation of a Lorentz-invariant theory on a lattice of length FIG. 6. (Color online) The imaginary part of G as a function N = 33. The right panel plots the corresponding Green’s of temporal diﬀerence T ∆m and spatial diﬀerence a∆n for a Gr m, n Gr r function (∆ ∆ ). is nonzero at some dispersed spots Lorentz-invariant theory on a lattice of length N = 73. G is which are in black, purple, red or yellow. nonzero only at a lattice of spacetime points (the black spots).

in the same way under the Lorentz boosts. Therefore, where (m ,n ) (m′ ,n′ ) or (m ,n ) (m′ ,n′ ) is 1 1 → 1 1 2 2 → 2 2 ˜ ˜ the coordinate transformation P (j,m,n). Since Gr is the set k, Ek˜ is Lorentz-invariant. The retarded independent of the quantum state, we obtain Green’sn functiono tells us how a single particle propagate in the spacetime. According to fig. 6, the particle starts r r G (m1,n1,m2,n2)= G (m1′ ,n1′ ,m2′ ,n2′ ) (33) from x =0at t = 0, and jumps from one site to the other after one entire period of evolution. The wave package for m1 m2 and m1′ m2′ . Eq. (33) puts a set of keeps localized as t equals integer times of T . constraints≥ on Gr as we≥ consider different P (j,m,n). r Fig. 7 the right panel gives another example of the The translational symmetry requires that G should de- Green’s function. We choose N = 33 and the correspond- pend only upon the difference of coordinates between r ing dispersion relation is plotted in the left panel. Again, two spacetime points. For simplicity, G (m1,n1,m2,n2) r r Ek =2π/T E2π/a k is satisfied, and then G is purely is rewritten as G (∆m, ∆n) with ∆m = m1 m2 and imaginary. −Gr equals− i, 2i/3, i/3 and i/3 at the ∆n = n n . According to the above discussions,− we − − − 1 − 2 spots colored black, purple, red and yellow, respectively, constrain the arguments ∆m and ∆n to be in the do- but zero everywhere else. According to the Lorentz sym- main [0,N 1). More important, the Lorentz symmetry − metry, the set of black spots must be Lorentz-invariant, requires so are the sets of other colored spots. The difference be- r r r tween fig. 6 and 7 is that G of the latter is nonzero at G (∆m, ∆n)= G (∆m′, ∆n′), (34) multiple sites for a given ∆m. This means that the wave package spreads over a few sites after the evolution over where the map (∆m, ∆n) (∆m′, ∆n′) is defined by Eq. (17). → one period. The black spots reappear at ∆m = 11 or 22. The retarded Green’s function evaluates The wave packages recombine into a single peak located at ∆n = 0 after every 11 periods. N 1 − r 1 i2π(k˜∆n E˜˜ ∆m)/N G (∆m, ∆n)= i e − k . (35) − N X˜ k=0 VII. SUMMARY The Green’s function is uniquely defined. Choosing the quasi-energy-momentum in different Brillouin zone may In this paper, we review the discrete Poincarégroup ˆ r change the definition of HF , but it leaves G unchanged. and systematically develop the quantum theories with Fig. 6 plots the imaginary part of Gr(∆m, ∆n) for N = discreteP Poincarésymmetry on the one-dimensional lat- 73 and the dispersion relation fig. 3(f). This dispersion tice of constant a and length N. contains discrete P relation satisfies Ek = 2π/T E2π/a k, which directly Lorentz transformations and discrete spacetime transla- r − − r indicates G being purely imaginary. We see that G is tions. The latter forms a Bravais lattice on the t-x plane. zero almost everywhere except for N specific (∆m, ∆n) The elements of are represented by unitary operators. where Gr = i (the black spots). The distribution of The common eigenstateP of translations is defined to be − the black spots has the Lorentz symmetry, as required a single-particle state with definite quasi-momentum k by Eq. (34). Their positions are indeed (∆m, ∆n) = and quasi-energy E. The discrete nature of translations ˜ ˜ confines (E, k) in a single Brillouin zone - a primitive k, Ek˜ . This is not surprising, because we have shown cell of the reciprocal lattice. The multiplication law of ˜ ˜ below Eq. (17) that k, Ek˜ and (∆m, ∆n) transform group determines how E and k transform under a dis- 11 crete Lorentz boost. It is similar to the way that energy may be much smaller than the speed of light. As we and momentum transform under a continuous Lorentz have mentioned above, the value of c has no influence transformation. But (E, k) must be mapped back into on our theory and then can be chosen freely. To realize the Brillouin zone after the transformation. The dimen- the Lorentz invariance, the hopping between lattice sites sion of the single-particle Hilbert space is N. A particle needs to be carefully tuned. The hopping as a function of can be localized on one of the N sites in real space. The distance is highly irregular in the models we have studied, N localized states are orthogonal to each other, at the so that it is not very probable in experiments. But we same time, the single-particle Hilbert space keeps invari- have only explored a small fraction of Lorentz-invariant ant under a Lorentz boost. These two conditions require lattice models by limiting ourselves on a short chain that k and E be integer times of 2π/(Na) and 2π/(NT ), (N < 100) with a fixed generator (γ(1) = 2). Whether respectively, where T is the period of time translation, there exists a realizable model on a longer chain or with and the dispersion relation E be Lorentz-invariant. For larger generators (γ(1) = 3, 4, ) is an interesting open k ··· a given N, each Lorentz-invariant Ek defines a quantum problem. Finally, it is possible to add an onsite particle- theory. But for some N, no Lorentz-invariant Ek can be particle interaction to HˆF without breaking the Lorentz found, thereafter, no Lorentz-invariant quantum lattice symmetry. It has been known that a δ-interaction does theory exists. We enumerate the 28 possible N between not break the continuous Lorentz symmetry, because two 2 and 100 for γ(1) = 2 where γ(1) is the generator of particles located at the same spacetime point in one ref- the discrete Lorentz transformations. The corresponding erence frame are also at the same point in the other ref- Ek displays reflection symmetry with respect to k = π/a erence frames. Developing an interacting lattice theory or E = π/T . For some N, there exist a family of saw- with discrete Poincarésymmetry is a challenging work in tooth dispersion relations in which the pairs (Ek, k) form future. a lattice in the Brillouin zone. After we obtain the single-particle dispersion relation, it is straightforward to write down a many-body theory ACKNOWLEDGEMENT of indistinguishable fermions or bosons. The field operators are defined both in real space and in momentum Pei Wang is supported by NSFC under Grant space. The unitary operators of , i.e. the translations Nos. 11774315 and 11835011, and by the Junior Asso- and Lorentz boosts, are expressedP in terms of the field ciates program of the Abdus Salam International Center operators. Especially, the discrete time-translational op- for Theoretical Physics. erator is also the evolution operator over one period, from which we derive the Floquet Hamiltonian HˆF . HˆF depends on Ek, being quadratic before we consider the Appendix A: Formulas about γ(j) and ζ(j) particle-particle interaction. In real space, HˆF includes the onsite potential and hopping term. To preserve the We list some useful formulas about γ(j) and ζ(j) with- Lorentz symmetry, there usually exists long-range hop- out proof. ζ(j) can be redefined as ping and the coupling strength fluctuates with the distance between sites. Some specific models are studied. 1 j ζ(j) γ(1)2 1= γ(1) + γ(1)2 1 An interesting case happens at N = 18, in which the hop- 2 p − p − (A1) ping is nonzero only between two sites at a distance of 3 1 j or 9. We calculate the retarded Green’s function Gr (the γ(1) γ(1)2 1 . − 2 − − propagator), which is the correlation between two field p operators at different positions and times and depends γ(j) (ζ(j)) is an even (odd) function of j, respectively. only upon the difference ∆x and ∆t. Under the Lorentz In other words, γ( j) = γ(j) but ζ( j) = ζ(j). The boost, a field operator in real spacetime transforms in the joint iterative relations− are − − same way as a spacetime point. As a result, the discrete Lorentz symmetry manifests itself in Gr which takes the ζ(j +1)= γ(1)ζ(j)+ γ(j), (A2) same value at a class of (∆t, ∆x) that map into each ζ(j 1) = γ(1)ζ(j) γ(j). other under the Lorentz transformations. In the model − − r with a sawtooth dispersion relation, G is zero almost ζ(j) satisfies a more complicated relation: everywhere except for a lattice of (∆t, ∆x) which has the same shape as the lattice of (Ek, k) in the Brillouin zone. ζ(i + j +1)= ζ(i + 1)ζ(j + 1) ζ(i)ζ(j), (A3) The particle propagates like a soliton even in the absence − of any nonlinear mechanism. or in a matrix form: The ultracold atomic gas in the periodically-driven optical lattice is a possible candidate for observing the dis- ζ(i + 1) ζ(i) ζ(j + 1) ζ(j) ζ(i) ζ(i 1) ζ(j) ζ(j 1) crete Poincarésymmetry. The lattice constant a and − − − − − − (A4) driving period T can be chosen freely. Consequently, the ζ(i + j + 1) ζ(i + j) = . speed of constant becomes c = γ(1)2 1 a/T , which ζ(i + j) ζ(i + j 1) − − − − p 12

Appendix B: Lorentz transformation of the We will prove that Ek′ = E2π/a k and Ek′′ =2π/T Ek − − quasi-energy-momentum are both Lorentz-invariant if Ek is Lorentz-invariant. The relations Ek′ = E2π/a k and Ek′′ = 2π/T Ek can − − Solving Eq. (10) becomes easier if we use the tensor be translated into E˜ = E˜ and E˜ = N E˜ , re- k˜′ N k˜ k˜′′ k˜ notation. From Eq. (10), we immediately obtain spectively. Therefore, we only− need to prove− that the ˜ ˜ ˜ ˜ µ µ sets EN k˜, k and N Ek˜, k are both invari- (me + nf)µ p + l 2π = (m′e + n′f)µ p′ , (B1) n − o n − o ant under the map (C1). where the lower and upper indices denote the compo- For convenience in next, we use to denote the nents of a tensor, and l is an arbitrary integer. The L Lorentz transformation in the tensor form is written as map (C1). Suppose E,˜ k˜ is an element of the set µ µ L ν . L ν is a collection of 4 numbers, which constitutes ˜ ˜ ˜ ˜ ˜ ˜ { } Ek˜, k and maps E, k into E′, k′ . Then the Lorentz matrix (2). In Eq. (6), the map from (m,n) n o L ˜ ˜ ˜ ˜ to (m′,n′) is indeed deﬁned by the coordinate transfor- E′, k′ must be an element of Ek˜, k . The con- mation of the Lorentz boost, reading n o verse is true. If E˜′, k˜′ is an element of the set, then µ µ ν (m′e + n′f) = L (me + nf) . (B2) ν E,˜ k˜ is also an element. This is because is invert- L ′ ′ µ ν ν ible. And the inverse map can be obtained by solving Using L ν Lµ = δν , we obtain Eq. (C1), which reads ν (me + nf)µ = L µ (m′e + n′f)ν . (B3)

Substituting Eq. (B3) into Eq. (B1), we have 2 E˜ γ(1)E˜′ γ(1) 1 k˜′ (mod N) , ≡ − − (C2) µ µ ν (p′ L p ) (m′e + n′f) = l 2π. (B4) k˜ E˜′ + γ(1)k˜′ (mod N) . − ν µ ≡−

(m′e + n′f)µ is the component of a lattice vector of . According to the deﬁnition of the reciprocal lattice Y, µ R From Eq. (C2), we immediately derive if the components of a vector r satisfy r (m′e + n′f)µ being an integer times of 2π, r must be a lattice vector of , i.e. R 2 (p µ Lµ pν ) . (B5) E˜ γ(1)E˜′ + γ(1) 1 N k˜′ (mod N) , ′ ν ≡ − − − ∈ R Rewriting Eq. (B5) in terms of components, we obtain N k˜ E˜′ + γ(1) N k˜′ (mod N) . − ≡ − Eq. (11). (C3)

Appendix C: A theorem on the Lorentz-invariant dispersion relation If we compare Eq. (C1) with Eq. (C3), the latter indeed

tells us that maps E˜′,N k˜′ into E,N˜ k˜ . And Ek is a Lorentz-invariant dispersion relation, if Ek L − − 1 ˜ ˜ ˜ ˜ and k can be expressed as E = 2πE˜ /(TN) and k = − then maps E,N k into E′,N k′ . There- k k˜ L − − ˜ ˜ ˜ 2πk/(aN), respectively, where Ek˜ and k are integers, and 1 ˜ ˜ fore, or − maps each element of the set EN k˜, k ˜ ˜ L L n − o the set Ek˜, k keeps invariant under the map into another element in the same set. In other words, n o E˜ ˜, k˜ keeps invariant under . Similarly, we can 2 N k E˜′ γ(1)E˜ + γ(1) 1 k˜ (mod N) , n − o L ≡ − (C1) ˜ ˜ prove that N Ek˜, k also keeps invariant under . k˜′ E˜ + γ(1)k˜ (mod N) . − L ≡ n o

∗ [email protected] Sen(De), and U. Sen, Adv. Phys. 56, 243 (2007). 1 M. Greiner, O. Mandel, T. Esslinger, T. W. H¨ansch, and 4 S. M. Girvin and K. Yang, Modern Condensed Matter I. Bloch, Nature 415, 39 (2002). Physics (Cambridge University Press, 2019). 2 G. Jotzu, M. Messer, R. Desbuquois, M. Lebrat, T. 5 T. Kitagawa, E. Berg, M. Rudner, and E. Demler, Phys. Uehlinger, D. Greif, and T. Esslinger, Nature 515, 237 Rev. B 82, 235114 (2010). (2014). 6 M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 3 M. Lewenstein, A. Sanpera, V. Ahuﬁnger, B. Damski, A. (2010). 13

7 S. Weinberg, The Quantum Theory of Fields (Cambridge York, 1959). University Press and Beijing World Publishing Corpora- 9 L. Bombelli, J. Lee, D. Meyer, and R. D. Sorkin, Phys. tion, 2013). Rev. Lett. 59, 521 (1987). 8 E. P. Wigner, Gruppentheorie und ihre Anwendung auf die 10 S. Surya, Living Rev. Relativ. 22, 5 (2019). Quantenmechanik der Atomspektren (Braunschweig, 1931): 11 P. Wang, New J. Phys. 20, 023042 (2018). pp. 251-3 (English translation, Academic Press, Inc, New 12 X. Li, J. Chai, H. Zhu, and P. Wang, J. Phys.:Condens. Matter 32, 145402 (2020).