An Application of the Extended Global SO ( 3 ) × SO ( 3 ) × U ( 1 ) Symmetry of the Hubbard Model on a Square Lattice : The Spinon , η-Spinon , and c Fermion Description

In this paper we review recent results on the preliminary applications of the new-found extended global SO(3) × SO(3) × U(1) symmetry of the Hubbard model on a bipartite lattice. Our results refer to the particular case of the bipartite square lattice. Specifically, we review a general description for such a model with nearest-neighbor transfer integral t and on-site repulsion U on a square lattice with N a 1 sites consistent with its extended global symmetry. It refers to three types of elementary objects whose occupancy configurations generate the state representations of the model extended global symmetry. Such objects emerge from a suitable electron-rotated-electron unitary transformation. An application to the spin spectrum of the parent compound La2CuO4 is shortly reviewed.


Introduction
The Hubbard model on a bipartite lattice is the simplest realistic toy model for description of the electronic correlation effects in general many-electron problems with short-range interaction on such a lattice.The model involves two effective parameters: the in-plane nearest-neighbor transfer integral t and the effective on-site repulsion U .Despite that it is among the mostly studied models in condensed matter and ultra-cold atom physics, except for the one-dimensional (1D) bipartite lattice [1][2][3] there is no exact solution and few controlled approximations exist for finite U/4t values.
An exact result for the Hubbard model on any bipartite lattice is that in addition to the spin SU (2) symmetry it has a second global SU (2) symmetry [4], called by some authors and in this paper η-spin symmetry [5].A trivial result is that at vanishing onsite interaction U = 0 the global symmetry of the Hubbard model on a bipartite lattice at vanishing chemical potential and magnetic field is O(4) = SO(4) × Z 2 .Here the factor Z 2 refers to the particle-hole transformation on a single spin under which the model Hamiltonian is not invariant for U = 0 and SO(4) = [SU (2) × SU (2)]/Z 2 contains the two SU (2) symmetries.Yang and Zhang considered the most natural possibility that the SO(4) symmetry inherited from the U = 0 Hamiltonian O(4) = SO(4) × Z 2 symmetry was the model global symmetry for U = 0 [5].
However, a recent study of the problem by the authors and collaborator reported in Reference [6] reveals an exact extra hidden global U (1) symmetry emerging for U = 0 in addition to SO (4).It is related to the U = 0 local SU (2) × SU (2) × U (1) gauge symmetry of the Hubbard model on a bipartite lattice with vanishing transfer integral [7].Such a local SU (2)×SU ( 2)×U (1) gauge symmetry becomes for finite U and t a group of permissible unitary transformations.It is such that the corresponding local U (1) canonical transformation is not the ordinary U (1) gauge subgroup of electromagnetism.Instead it is a "nonlinear" transformation [7].
Consistently, for U = 0 the related new found global symmetry of the model on any bipartite lattice is larger than SO(4) and given by [SU (2)×SU (2)×U (1)]/Z 2 2 = [SO(4)×U (1)]/Z 2 = SO(3)×SO(3)× U (1).The factor 1/Z 2 (and 1/Z 2  2 ) in SO(4) = [SU (2)×SU (2)]/Z 2 (and in [SU (2)×SU (2)×U (1)]/Z 2 2 ) imposes that [S η + S s ] is an integer number (and both [S s + S c ] and [S η + S c ] are integer numbers).Here S η , S s , and S c are the η-spin, the spin, and one-half the eigenvalue 2S c of the generator of the new global U (1) symmetry, respectively.The latter is found in Reference [6] to be the number of rotated-electron singly occupied sites.This refers to any of the infinite electron-rotated-electron unitary transformations of Reference [8], such that rotated-electron single and double occupancy are good quantum numbers for U/4t = 0. 2S c is then the number of rotated-electron singly occupied sites.Within the present notation, S An important point is that although addition of chemical-potential and magnetic-field operator terms to the Hubbard model on a bipartite lattice Hamiltonian lowers its symmetry, such terms commute with it.Therefore, the global symmetry being SO(3) × SO(3) × U (1) implies that the set of independent rotated-electron occupancy configurations that generate the model energy and momentum eigenstates generate state representations of that global symmetry for all values of the electronic density n and spin density m.It then follows that the total number of such independent representations must equal the Hilbert-space dimension, 4 N D a .The results of Reference [6] confirm that for the model on a bipartite lattice in its full Hilbert space the number of independent representations of the group SO(3) × SO(3) × U (1) is indeed 4 N D a .In contrast, the number of independent representations of the group SO(4) is found to be smaller than the Hilbert-space dimension 4 N D a .In this paper we review an application of the extended global SO(3) × SO(3) × U (1) symmetry of the Hubbard model on a square lattice.It refers to a rotated-electron related operator description in terms of spin-1/2 spinons, η-spin-1/2 η-spinons, and spin-less and η-spin-less c fermions [9,10].Such an operator description involves an electron-rotated-electron unitary transformation of the type considered in Reference [8].A property specific to the specific transformation associated with the operator description reviewed in this paper is that the U > 0 energy eigenstates can be generated from suitable chosen U/4t → ∞ energy eigenstates upon application onto the latter states of the corresponding electron-rotated-electron unitary operator.The related set of state representations of the group SO(3) × SO(3) × U (1) that emerge from our description are energy eigenstates yet are generated by exactly the same electron-rotated-electron unitary transformation from corresponding U/4t → ∞ states.
Their rotated electron occupancy configurations are simpler to describe in terms of the above-mentioned three types of elementary objects directly related to the rotated electrons whose numbers and designations are: a −2S c ] sites doubly occupied and unoccupied by rotated electrons.Specifically, the η-spinons of η-spin projection −1/2 and +1/2 refer to the sites doubly occupied and unoccupied, respectively, by rotated electrons.The remaining degrees of freedom of such rotated-electron occupancy configurations are described by the N h c = [N D a − 2S c ] c fermion holes.The expression of the rotated-electron operators in terms of these three elementary objects is for U > 0 identical to that of the electron operators for large U values in terms of the three objects obtained from the exact transformation for separation of spin-1/2 fermions without constraints considered in Reference [11].
Interestingly, η-spinon (and spinons) that are not invariant under the electron-rotated-electron unitary transformation considered in this paper have η-spin 1/2 (and spin 1/2) but are anti-bound within η-spin-neutral (and bound within spin-neutral) 2ν-η-spinon (and 2ν-spinon) composite ην fermions (and sν fermions).Here ν = 1, 2, . . . is the number of anti-bound η-spinon (and bound spinon) pairs.In this paper we follow the notation of Reference [9] and call such spinons and η-spinons, confined spinons and confined η-spinons, respectively.We emphasize though that by "confinement" is meant here that the spinons and η-spinons are bound and anti-bound within such composite sν and ην fermions, respectively, alike for instance protons and neutrons are bound within the nucleus.Whether the potential associated with such a behavior is in some limit confining remains an interesting open question.The reason for the notation of Reference [9] is that the term confined can be used both for spinons and η-spinons.A possible alternative notation would be bound spinons and anti-bound η-spinons, which would be different for spinons and η-spinons.We thus stress that within our notation the term confined is not equivalent to the high-energy physics concept of confinement.
Moreover, again we follow the notation of Reference [9] and call deconfined spinons and deconfined η-spinons those that are invariant under the electron-rotated-electron unitary transformation.Here deconfined means just that such objects are not bounded and anti-bounded within sν fermions and ην fermions, respectively.In case that future research reveals that the potential behind the energetic relations of sν fermions and ην fermions discussed below in Section 4.3 is at least in some limit confining, any possible relation to the concept of deconfined quantum criticality of Reference [12] should be clarified.
The paper is organized as follows.A uniquely defined rotated-electron description for the Hubbrad model on the square lattice and its relation to the global SO(3)×SO( 3)×U (1) symmetry are the subjects of Section II.In Section III the c fermion, η-spin-1/2 η-spinon, and spin-1/2 spinon and corresponding c, η-spin, and spin effective lattices are reviewed.The vacua of the theory, the transformation laws under the electron-rotated-electron unitary transformation of such objects, and the subspaces they refer to are issues also addressed in that section.The composite αν bond particles and αν fermions, corresponding αν effective lattices, ground-state occupancy configurations, and a complete set of momentum eigenstates are the problems studied in Section IV.The c and s1 fermion square-lattice quantum liquid and corresponding one-and two-electron subspace are the subjects discussed in Section V.Moreover, in that section a preliminary application of the present operator description to the inelastic neutron scattering of the Mott-Hubbard insulator parent compound La 2 CuO 4 (LCO) [13] is presented.Our description leads to simple expressions that agree both with the LCO inelastic neutron scattering and the results of Reference [14], which involved the summation of an infinite number of diagrams.Finally, the concluding remarks are presented in Section VI.The Hubbard model on a square (or 1D) lattice with a very large number N D a 1 of sites reads, Periodic boundary conditions and torus periodic boundary conditions are considered for the 1D lattice for which D = 1 and the square lattice for which D = 2, respectively.Moreover, in Equation (1) t is the nearest-neighbor transfer integral, T is the kinetic-energy operator in units of t, and Q is the operator that counts the number of electron singly occupied sites.Hence the operator D = [ N − Q]/2 counts the number of electron doubly occupied sites.Moreover, n r j ,σ = c † r j ,σ c r j ,σ where −σ =↑ (and −σ =↓) for σ =↓ (and σ =↑), N = σ Nσ , and Nσ = N D a j=1 n r j ,σ .We focus our attention onto ground states with hole concentration x = [N D a − N ]/N D a ≥ 0 and spin density m = [N ↑ − N ↓ ]/N D a = 0 and their excited states.Such an hole-concentration range corresponds to electronic densities n = N/N D a ∈ (0, 1).We are particularly interested in the LWS subspace spanned by the lowest-weight states (LWSs) of both the η-spin and spin algebras.Such energy eigenstates refer to values of S α and S x 3 α such that S α = −S x 3 α for α = η, s.
The kinetic-energy operator T can be expressed in terms of the operators, as T = T0 + T+1 + T−1 .These three kinetic operators play an important role in the physics.The operator T0 does not change electron double occupancy whereas the operators T+1 and T−1 change it by +1 and −1, respectively.The studies of Reference [6] consider unitary operators V = V (U/4t) and corresponding rotated-electron operators, Those are such that rotated-electron single and double occupancy are good quantum numbers for U/4t > 0. The global U (1) symmetry generator 2 Sc of eigenvalue 2S c reads [6], It follows that 2 Sc = V † 2 Ŝc V where 2 Ŝc = Q and the operator Q is given in Equation (1).As mentioned in the previous section, 2S c is the number of rotated-electron singly occupied sites.Most choices of the unitary operators V correspond to choices of U/4t → ∞ sets {|Ψ ∞ } of 4 N D a energy eigenstates such that the states |Ψ U/4t = V † |Ψ ∞ are not energy and momentum eigenstates for finite U/4t values yet belong to a subspace with fixed and well-defined values of S c , S η , and S s .
Let {|Ψ U/4t } be a complete set of 4 N D a energy, momentum, η-spin, η-spin projection, spin, and spin-projection eigenstates for U/4t > 0. In the limit U/4t → ∞ such states correspond to one of the many choices of sets {|Ψ ∞ } of 4 N D a U/4t-infinite energy eigenstates.Both the sets of states {|Ψ U/4t } and {|Ψ ∞ }, respectively, are complete and the model Hilbert space is the same for all U/4t > 0 values considered here.Hence it follows from basic quantum mechanics Hilbert-space and operator properties that for this choice there exists exactly one unitary operator Here we consider such a unitary operator and corresponding generator Sc = V † Ŝc V given in Equation ( 4) and rotated-electron operators provided in Equation (3).The states one for each value of U/4t > 0) that are generated from the same initial state |Ψ ∞ belong to the same V tower.
. Such states are generated by occupancy configurations of the quantum objects of the general description considered in this paper.The unitary operator V † appearing in the general expression of such states is what also appears in the general expression |Ψ U/4t = V † |Ψ ∞ of the above energy and momentum eigenstates.However the states |Φ U/4t are not in general energy eigenstates of the Hubbard model on the square lattice.The interest of the states |Φ U/4t = V † |Φ ∞ is that those contained in the one-and two-electron subspace defined below are both momentum and energy eigenstates of that model.Fortunately, the present general description is physically most useful and important for the Hubbard model on the square lattice in that subspace, which refers to the square-lattice quantum liquid further investigated in Reference [10].
Alike in Reference [6], we associate with any operator Ô an operator Õ = V † Ô V that has the same expression in terms of rotated-electron creation and annihilation operators as Ô in terms of electron creation and annihilation operators, respectively.Our convention is that marks placed over letters being a caret or a tilde denote operators.(An exception are the electron operators of Equation ( 3), which we denote by c † r j ,σ and c r j ,σ rather than by ĉ † r j ,σ and ĉ r j ,σ , respectively.)Any operator Ô can then be written in terms of rotated-electron creation and annihilation operators as, The operator Ŝ appearing in this equation is related to the unitary operator as Although the general expression of Ŝ remains unknown, an exact result is that it involves only the kinetic operators T0 , T+1 , and T−1 of Equation ( 2) and numerical U/4t dependent coefficients [8,9].For U/4t > 0 that expression can be expanded in a series of t/U .The corresponding first-order term has the universal form given in Equation (5).To arrive to the expression of Ô in terms of the operator S also given in Equation ( 5), the property that the operator V commutes with itself is used.It implies that V = e − Ŝ = Ṽ = e − S and Ŝ = S. Hence both the operators V and Ŝ have the same expression in terms of electron and rotated-electron creation and annihilation operators.This justifies why the expansion given in that equation for the operator S has the same form as that of Ŝ.The higher-order terms of the operator S expression can be written as a product of operator factors.Their expressions involve the rotated kinetic operators T0 , T+1 , and T−1 .The full expression of the operator Ŝ = S can for U/4t > 0 be written as Here Ŝ(∞) = S(∞) corresponds to the operator S(l) at l = ∞ defined in Equation (61) of Reference [8].Moreover, δ Ŝ = δ S has the general form provided in Equation (64) of that paper.The unitary operator V = e − Ŝ = e − S considered here corresponds to exactly one choice of the coefficients D (k) (m) of that equation.Here the index k = 1, 2, . . .refers to the number of rotated-electron doubly occupied sites.The problem of finding the explicit form of the operators V = Ṽ and Ŝ = S is equivalent to finding all coefficients D (k) (m) associated with the electron-rotated-electron unitary transformation as defined above.
For finite U/4t values the Hamiltonian Ĥ of Equation ( 1) does not commute with the unitary operator V = e − Ŝ .Hence when expressed in terms of the rotated-electron creation and annihilation operators of Equation ( 3) it has an infinite number of terms.According to Equation (5) it reads, The commutator [ H, S ] does not vanish except for U/4t → ∞ so that Ĥ = H for finite values of U/4t.Fortunately, there is strong evidence that for approximately U/4t > 1, out of the infinite terms on the right-hand-side of Equation ( 6) only the first few Hamiltonian terms play an active role in the physics of the Hubbard model on the square lattice in the one-and two-electron subspace [10].This follows in part from in that subspace the rotated-electron configurations expressed in Reference [9] in terms of related c and s1 fermion operators referring to energy eigenstates.For instance, the spin spectrum used below in Section 5.4 was derived in Reference [10] under the assumption that for approximately U/4t > 1 only the first few Hamiltonian terms on the right-hand side of Equation ( 6) play an active role in the model physics.For U/4t = 1.525 the validity of that assumption is consistent with the very good agreement with both results obtained by the standard formalism of many-body physics for the model spin spectrum [14] and the inelastic neutron scattering of LCO reported in that Section.Note, however, that this property is expected to hold provided that the unitary transformation used to generate the rotated electrons from the electrons is that considered above.Alike the operator S, the Hamiltonian expression in terms of rotated-electron operators Equation ( 6) can be expanded in a series of t/U .Its terms generated up to fourth order in t/U are within a unitary transformation the equivalent to the t−J model with ring exchange and various correlated hoppings [15].Furthermore, at half filling the terms of the Hamiltonian Equation ( 6) expansion in t/U with odd powers in t vanish due to the particle-hole symmetry and the resulting invariance of the spectrum under t → − t.In turn, for finite hole concentration x > 0 the expansion in powers of t/U of the Hamiltonian Equation ( 6) involves terms with odd powers in t, absent at x = 0.This is consistent with the effects of increasing U/4t being often different at x = 0 and for x > 0 [16].
It is shown in Reference [8] that all infinite electron-rotated-electron unitary transformations are well defined except at U = 0.This applies to the specific choice of electron-rotated-electron unitary transformation associated with the operational description reviewed in this paper.Given the unitarity character of that transformation, it may be used in two different physical problems, which however are mathematically and technically fully equivalent.These two physical problems are considered in Reference [8] and References [9,10], respectively.Specifically, in Reference [8] it is considered that the rotated creation and annihilation operators of Equation (3) refer to electrons.Within that choice, except for U/4t → ∞ the Hamiltonian on the right-hand side of Equation ( 6) is not the Hubbard Hamiltonian.Instead it is an Hamiltonian for which electron double occupancy and single occupancy are good quantum numbers.In turn, within the alternative physical problem studied here and in References [9,10], the rotated creation and annihilation operators of Equation (3) refer to rotated electrons, which except for U/4t → ∞ are objects different from electrons.In this case, provided that both the on-site interaction U = 0 is finite and one accounts for all higher-order terms on the right-hand-side of Equation ( 6), the corresponding expression refers to the Hubbard model in terms of rotated-electron creation and annihilation operators.Hence within the operational description reviewed in this paper the general expression given in Equation ( 6) is not a mapping of the Hubbard model onto some generalized t − J model.Indeed the electron-rotated-electron unitary transformation considered here has been constructed to inherently rotated-electron double occupancy and singly occupancy being good quantum numbers for the Hubbard model with U = 0. Hence unlike in the derivation of t − J models of References [17,18], the range of applicability of the Hubbard model expression given in Equation ( 6) refers to all finite U = 0 values.
In Reference [6] it is found that, in contrast to the Hamiltonian, the three components of the momentum operator ˆ P , three generators of the spin SU (2) symmetry, and three generators of the η-spin SU (2) symmetry commute with the electron-rotated-electron unitary operator V = Ṽ .This also holds for the specific choice of that operator associated with the rotated-electron description considered in this paper.Hence the above operators have the same expression in terms of electron and rotated-electron creation and annihilation operators, so that the momentum operator reads, Furthermore, the above-mentioned six generators are given by, However, except for U/4t → ∞ the local rotated-electron operators, [ñ r j ,↑ − ñ r j ,↓ ] ; s+ r j ,s = c † r j ,↓ c r j ,↑ ; s− r j ,s = c † r j ,↑ c r j ,↓ , j = 1, 2, ..., N D a (9) and corresponding local electron operators, [n r j ,↑ − n r j ,↓ ] ; ŝ+ r j ,s = ĉ † r j ,↓ ĉ r j ,↑ ; ŝ− r j ,s = ĉ † r j ,↑ ĉ r j ,↓ , j = 1, 2, ..., N D a (10) appearing in the expressions given in Equation ( 8) of such six generators are different operators.We emphasize though that in the U/4t → ∞ limit such two sets of six local operators become the same operators and in addition are the six out of seven generators of the local SU (2) × SU (2) × U (1) gauge symmetry reported in Reference [7].The vector π appearing in Equations ( 9) and (10) has Cartesian components π = [π, π] (and component π) for the model on the square (and 1D) lattice.
In contrast, the generator 2 Sc of the charge independent U (1) symmetry given in Equation (4) does not commute with the unitary operator V .This is behind the hidden character of such a symmetry.On the contrary of the Hamiltonian, that generator has a complicated expression in terms of electron creation and annihilation operators and a simple expression given in that equation in terms of rotated-electron creation and annihilation operators.The operator of Equation (4) plus the six operators provided in Equation ( 8) are the seven generators of the group [SO(4) × U (1)]/Z 2 = SO(3) × SO(3) × U (1) associated with the global symmetry of the Hamiltonian Equation (1).

Three Elementary Quantum Objects and Corresponding c, η-Spin, and Spin Effective Lattices
Within the operational description reviewed this paper, three rotated-electron related elementary objects naturally emerge from the local operators that in the U/4t → ∞ limit generate the local SU (2)×SU (2)×U (1) gauge symmetry of Reference [7].Summation over the N a sites of such operators gives for U/4t > 0 the generators of the corresponding model extended global symmetry.Such elementary objects are mapped from the rotated electrons by an exact local transformation that does not introduce constraints.Given their direct relation to the generators of the model extended symmetry, their occupancy configurations naturally generate representations of the global symmetry algebra.
In the U/4t → ∞ limit and except for unimportant phase factors, the c fermion creation operators and the spinon and η-spinon operators introduced below become the quasicharge annihilation operators and the local spin and pseudospin operators, respectively, of Reference [11].The representation associated with the latter operators refers to all dimensions and lattices.For instance, in Reference [11] it is used in the study of the Hubbard model on the square lattice.The unitarity and remaining properties of the electron-rotated-electron transformation considered in this paper assures that the corresponding c fermion, spinon, and η-spinon operator representation introduced in the following for U/4t > 0 applies as well to the Hubbard model on the square lattice.

Elementary Quantum Objects and Their Operators
The local c fermion annihilation operator f r j ,c and corresponding creation operator f † r j ,c = (f r j ,c ) † are constructed in terms of the rotated-electron operators of Equation (3).Within the LWS representation the latter operator reads, In Equation (11) we have introduced as well the corresponding c fermion momentum-dependent operators and e i π• r j is ±1 depending on which sub-lattice site r j is on.(For the 1D lattice that phase factor can be written as (−1) j .)The c momentum band is studied in Reference [10] and has the same shape and momentum area as the electronic first-Brillouin zone.
The three spin local rotated operators sl r j ,s and the three η-spin local rotated operators sl r j ,η such that l = ±, x 3 , given in Equation ( 9) are associated with the spinons and η-spinons, respectively, as defined in this paper.Consistently with the f † r j ,c expression given in Equation (11), the local rotated operator s r j ,c and such spin and η-spin local rotated operators may be rewritten as, Here q± r j = qx 1 r j ± i qx 2 r j where x 1 , x 2 , x 3 denotes the Cartesian coordinates of the rotated quasi-spin operators qx i r j with i = 1, 2, 3.The former rotated quasi-spin operators ql r j = sl r j ,s + sl r j ,η such that l = ±, x 3 have for U/4t > 0 the following expression in terms of rotated-electron operators, Since the electron-rotated-electron transformation generated by the operator V is unitary, the operators c † r j ,σ and c r j ,σ of Equation ( 3) have the same anticommutation relations as c † r j ,σ and c r j ,σ , respectively.Straightforward manipulations based on Equations ( 11)-( 13) then lead to the following algebra for the c fermion operators, and c fermion operators and rotated quasi-spin operators, In turn, the rotated quasi-spin operators q x 1 r j , q x 2 r j , and q x 3 r j and the related operators q ± r j = q x 1 r j ± i q x 2 r j obey the following algebra, [q + r j , q − r j ] = δ j,j q x 3 r j ; [q ± r j , q ± r j ] = 0 Hence the operators q ± r j anticommute on the same site and commute on different sites.The same applies to three spinon operators s l r j and three η-spinon operators p l r j , respectively, whose expressions are given in Equation (12).
The relations provided in Equations ( 14)- (18) confirm that the c fermions associated with the hidden global U (1) symmetry are η-spinless and spinless fermionic objects.They are consistent as well with the spinons and η-spinons being spin-1/2 and η-spin-1/2 objects, respectively, whose local operators obey the usual corresponding SU (2) algebras.
On inverting the relations given in Equations ( 11) and ( 13), the rotated-electron creation and/or annihilation operators of Equation (3) are written in terms of the operators of the c fermions and rotated quasi-spin operators.For the LWS subspace considered here this leads to, The relations given in Equations ( 11)-( 13) and (19) between the c fermion, spinon, and η-spinon operators and the rotated-electron operators differ from those of References [9,10] by unimportant phase factors.For U/4t → ∞ the rotated electrons become electrons.In that limit and as mentioned above, except for unimportant phase factors the c fermion creation operators become the quasicharge annihilation operators of Reference [11] and the spinon and η-spinon operators become the local spin and pseudospin operators, respectively, of that reference.Consistently, in that limit and again unimportant phase factors Equations ( 11)-( 19) are equivalent to Equations ( 1)-(3) of Reference [11] with the rotated-electron operators replaced by the corresponding electron operators and the c fermion creation operator f † r j ,c replaced by the quasicharge annihilation operator ĉr .Alike the transformation considered in Reference [11], that considered here does not introduce Hilbert-space constraints.Hence suitable occupancy configurations of c fermions, spinons, and η-spinons generate a complete set of finite-U/4t states of the form |Φ U/4t = V † |Φ ∞ .Those are both state representations of the model global SO(3) × SO(3) × U (1) symmetry and momentum eigenstates.In general the energy and momentum eigenstates |Ψ U/4t = V † |Ψ ∞ of the Hubbard model on the square lattice are a superposition of a sub-set of states |Φ U/4t with the same momentum eigenvalue.

Interplay of the Global Symmetry with the Transformation Laws Under the Operator V : Three Basic Effective Lattices and the Theory Vacua
It follows from the results of Reference [6] that the η-spin SU (2) and spin SU (2) state representations correspond to the η-spin and spin degrees of freedom of independent rotated-electron occupancy configurations of [N D a − 2S c ] sites and 2S c sites, respectively, of the original lattice.This applies to any choice of electron-rotated-electron unitary transformation and thus applies as well to that associated with the present description.In turn, the state representations associated with the new-found hidden U (1) 2 symmetry refer to the relative positions in the original lattice of the sites involved in each of these two types of configurations.In the present N D a 1 limit it is useful to consider the following numbers, such that N D a = N D aη + N D as .Here the integer numbers N D aη and N D as are the number of sites of the original lattice singly occupied by rotated electrons and unoccupied plus doubly occupied by rotated electrons, respectively.Below such numbers are found to play the role of number of sites of a η-spin and spin effective lattice, respectively.For the D = 2 square lattice the number N 2 a = N a × N a is chosen so that the number N a of sites in a row or column is an integer.However, the designations N D aη and N D as do not imply that the corresponding numbers N aη and N as are integers.In general they are not integers.For finite values of x and (1 − x) and N D a → ∞ we use for the numbers N aη and N as of sites in a row and column of the η-spin and spin effective square lattices, respectively, the corresponding closest integer numbers.
Importantly, the numbers given in Equation ( 20) are fully controlled by the eigenvalue 2S c of the generator Equation (4) of the hidden global U (1) symmetry.This confirms the important role played by such a symmetry in the present general description.The degrees of freedom of the rotated-electron occupancy configurations of each of the sets of N D aη = [N D a − 2S c ] and N D as = 2S c sites of the original lattice that generate the S η , S x 3 η , S s , S x 3 s , S c , and momentum eigenstates studied in Reference [9], which refer to state representations of the model global SO(3) × SO(3) × U (1) symmetry, naturally separate as follows: (i) The occupancy configurations of the c fermions associated with the operators f † r j ,c = (f r j ,c ) † of Equation (11) where j = 1, ..., N D a correspond to the state representations of the hidden global U (1) symmetry found in Reference [6].Such c fermions live on the c effective lattice.It is identical to the original lattice.Its occupancies are related to those of the rotated electrons: The number of c fermion occupied and unoccupied sites is given by N respectively.Indeed, the c fermions occupy the sites singly occupied by the rotated electrons.In turn, the rotated-electron doubly-occupied and unoccupied sites are those unoccupied by the c fermions.Hence the c fermion occupancy configurations describe the relative positions in the original lattice of the The former configurations refer to the operators p l r j of Equation ( 12), which act only onto the ] sites of the η-spin effective lattice.The latter configurations correspond to the operators s l r j given in the same equation, which act onto the N D as = 2S c sites of the spin effective lattice.This is assured by the operators (1 − n r j ,c ) and n r j ,c in their expressions provided in that equation, which play the role of projectors onto the η-spin and spin effective lattice, respectively.
For U/4t > 0 the degrees of freedom of each rotated-electron singly occupied site then separate into a spin-less c fermion carrying the electronic charge and a spin-down or spin-up spinon.Furthermore, the degrees of freedom of each rotated-electron doubly-occupied or unoccupied site separate into a η-spin-less c fermion hole and a η-spin-down or η-spin-up η-spinon, respectively.The η-spin-down or η-spin-up η-spinon refers to the η-spin degrees of freedom of a rotated-electron doubly-occupied or unoccupied site, respectively, of the original lattice.
Hence a key feature of the operational description reviewed in this paper is that for U/4t > 0 its quantum objects correspond to rotated-electron configurations whose numbers of spin-down and spin-up single occupied sites, double occupied sites, and unoccupied sites are good quantum numbers.This is in contrast to descriptions in terms of electron configurations, whose numbers of spin-down and spin-up single occupied sites, double occupied sites, and unoccupied sites are good quantum numbers only for U/4t 1 [19].The transformation laws under the electron-rotated-electron unitary transformation of the η-spinons (and spinons) play a major role in the description of the η-spin (and spin) SU (2) state representations in terms of the occupancy configurations of the ] sites of the η-spin effective lattice (and set of N D as = 2S c sites of the spin effective lattice).Indeed, a well-defined number of η-spinons (and spinons) remain invariant under that unitary transformation.Those are called deconfined ±1/2 η-spinons (and deconfined ±1/2 spinons).As further discussed below, they play the role of unoccupied sites of the η-spin (and spin) effective lattice.The values of the numbers M de η,±1/2 of deconfined ±1/2 η-spinons and M de s,±1/2 of deconfined ±1/2 spinons are fully controlled by the η-spin S η and η-spin projection S x 3 η = − x 2 N D a and spin S s and spin projection S x 3 s = − m 2 N D a , respectively, of the state under consideration as follows, The invariance of such deconfined η-spinons (and deconfined spinons) stems from the off diagonal generators of the η-spin (and spin) algebra, which flip their η-spin (and spin), commuting with the unitary operator V .This justifies why such generators have for U/4t > 0 the same expressions in terms of electron and rotated-electron operators, as given in Equation ( 8).It follows that the number of sites of the η-spin (α = η) and spin (α = s) effective lattice can be written as, As justified below, here M co α is the confined η-spinon (α = η) or confined spinon (α = s) even number.The η-spin or spin degrees of freedom of the occupancy configurations in a set of 2S α = M de α sites out of the whole set of N D aα sites have for U/4t > 0 the same form in terms of electrons and rotated electrons.In turn, for finite values of U/4t and spin density −(1 − x) < m < (1 − x) the corresponding c fermion occupancy configurations that store the information on the relative positions in the original lattice of the N D as sites of the spin effective lattice and sites of the η-spin effective lattice are not invariant under the electron-rotated-electron unitary transformation.
Such an invariance of the η-spin degrees of freedom of the above occupancy configurations implies that in each state representation there are exactly sites are doubly occupied and M de η,+1/2 sites are unoccupied both by electrons and rotated electrons.Furthermore, the invariance of the spin degrees of freedom of the sites singly occupied by rotated electrons implies that there are exactly 2S s = [M de s,−1/2 + M de s,+1/2 ] = M de s sites of the original lattice such that M de s,−1/2 sites are singly occupied both for spin-down electrons and spin-down rotated electrons and M de s,+1/2 sites are singly occupied both for spin-up electrons and spin-up rotated electrons.The state representation is in general a superposition of such occupancy configurations whose positions of the [2S η + 2S s ] = [M de η + M de s ] sites are different.In turn, out of the set of [M co η + M co s ] sites of the original lattice left over, a set of M co η /2 sites are unoccupied by rotated electrons, a set of M co η /2 sites are doubly occupied by rotated electrons, a set of M co s /2 sites are singly occupied by spin-up rotated electrons, and a set of M co s /2 sites are singly occupied by spin-down rotated electrons.However, in terms of electrons that set of [M co η + M co s ] sites of the original lattice has for finite U/4t values very involved occupancies.Indeed for electrons and except for U/4t → ∞ singly and doubly occupancy are not good quantum numbers and thus are not conserved.The important point brought about by the present description is that due to the transformation laws under the electron-rotated-electron unitary transformation this refers only to a sub-set of [M co η + M co s ] sites out of the N D a sites of the original lattice.Indeed for the remaining sites singly and doubly occupancy are good quantum numbers both for electrons and rotated electrons.
The site numbers M co η ≥ 0 and M co s ≥ 0 are good quantum numbers given by, Hence their values are fully determined by those of the eigenvalue 2S c of the hidden global U (1) symmetry generator and η-spin S η or spin S s , respectively.This reveals that M co η and M co s are not independent quantum numbers.
It follows from the above discussions that the total numbers of η-spinons (α = η) and spinons (α = s) read, The η-spinon and spinon operator algebra refers to well-defined subspaces spanned by states whose number of each of these basic objects is conserved and given in Equations ( 21), (23), and (24).Hence in such subspaces the number 2S c of rotated-electron singly occupied sites and the numbers N D aη and N D as of sites of the η-spin and spin effective lattices, respectively, are fixed.For hole concentrations 0 ≤ x < 1 and maximum spin density m = (1 − x) (or m = −(1 − x)) reached at some critical magnetic field H c aligned parallel to the square-lattice plane for D = 2 and pointing along the chain for D = 1 the c fermion operators are invariant under the electron-rotated-electron unitary transformation.There is a corresponding fully polarized vacuum |0 ηs that remains invariant under such a transformation.It reads, Here the η-spin SU (2) vacuum |0 η ; N D aη associated with N D aη deconfined +1/2 η-spinons, the spin SU (2) vacuum |0 s ; N D as with N D as deconfined +1/2 spinons, and the c U (1) vacuum |GS c ; 2S c with N c = 2S c c fermions remain invariant under the electron-rotated-electron unitary transformation.The latter vacuum may be expressed as q f † q,c |GS c ; 0 .Here the vacuum |GS c ; 0 corresponds to the electron and rotated-electron vacuum.Its form is that given in Equation ( 25 deconfined +1/2 η-spinons refer to N D aη = M de η,+1/2 sites of the original lattice unoccupied by rotated electrons (and the M s = M de s,+1/2 deconfined +1/2 spinons to the spins of N D as = M de s,+1/2 spin-up rotated electrons that singly occupy sites of such a lattice).At maximum spin density m = (1−x) the c fermions are the non-interacting spinless fermions that describe the charge degrees of freedom of the electrons of the fully polarized ground state.At that spin density there are no electron doubly occupied sites and the quantum problem is non-interacting for U/4t > 0.
Below we confirm that within the present description, out of the sites of the η-spin (α = η) and spin (α = s) effective lattice, the 2S α = M de α sites whose occupancies η-spin (α = η) and spin (α = s) degrees of freedom are invariant under electron-rotatedelectron unitary transformation play the role of unoccupied sites.In turn, the remaining M co α sites play the role of occupied sites.This is a natural consequence of the η-spin SU (2) vacuum |0 η ; N D aη (and spin SU (2) vacuum |0 s ; N D as ) being for all U/4t and m values invariant under the electron-rotated-electron unitary transformation.Indeed that vacuum is such that N D aη = 2S η = M de η (and N D as = 2S s = M de s ), so that M co η = 0 (and M co s = 0.)

Important Energy Scales
For processes associated with one-and two-electron excitations the large-U/4t energy scale U is for finite U/4t values and x = 0 replaced by twice the absolute value of chemical potential, 2|µ|.This is because for x ∈ (0, 1) (and x ∈ (−1, 0)) and finite U/4t values the minimum energy for creation onto the ground state of a rotated-electron doubly occupied site (and unoccupied site) is given by 2|µ|.In turn, at x = 0 the chemical potential µ belongs to the range µ ∈ (−µ 0 , µ 0 ) whose energy width 2µ 0 ≡ lim x→0 2µ equals the Mott-Hubbard gap.We use the convention that for x = 0 the chemical potential µ sign is that of the hole concentration x ∈ (−1, 1).For spin densities m = 0 we consider that the corresponding magnetic field points in directions such that H > 0 for spin density m > 0 and H < 0 for m < 0. The minimum magnitude of the energy ∆ Drot for creation of a number D rot = M η,−1/2 = M de η,−1/2 + M co η /2 of rotated-electron doubly occupied sites onto a m = 0 and x = 0 ground state and a m = 0 and x > 0 ground state is, respectively.Similar expressions apply to the minimum magnitude of the energy ∆ D h rot for creation of a number D h rot = M η,+1/2 = M de η,+1/2 + M co η /2 of rotated-electron unoccupied sites onto a m = 0 and x = 0 ground state and a m = 0 and x < 0 ground state, respectively, provided that µ is replaced by −µ and M de η,−1/2 by M de η,+1/2 .Hence for µ ∈ (−µ 0 , µ 0 ) the minimum energy for creation onto the m = 0 and x = 0 ground state of one rotated-electron doubly occupied site and rotated-electron unoccupied site is (µ 0 + µ) and (µ 0 − µ), respectively.For the chemical potential in the middle of the Mott-Hubbard gap so that µ = 0 at x = 0 the minimum energy of either process is given by µ 0 .
For the model on the 1D (and square) lattice at m = 0 the important x = 0 energy scale 2µ 0 appearing in Equation (26), which is the Mott-Hubbard gap, has the following exact (and approximate) limiting behaviors, These half-filling results are consistent with the properties of the x = 0 and m = 0 absolute ground state [10].In turn, for 0 < x < 1 and m = 0 the chemical potential reads µ ≈ U/2 for U/4t → ∞ and for finite U/4t values is an increasing function of the hole concentration x such that, Here µ 1 ≡ lim x→1 µ.The related energy scale 2µ 1 reads, The expression given in Equation ( 29) is exact both for the Hubbard model on a 1D and square lattice.It can be explicitly derived for both lattices.It refers to the non-interacting limit of vanishing electronic density.For the model on the square lattice the limiting behaviors reported in Equation ( 27) read 2µ 0 ≈ 64 t e −2π √ t/U and 2µ 0 ≈ U for U/4t 1 and U/4t 1, respectively.They are those of the related zero-temperature gap of Equation ( 13) of Reference [20], which although showing up in the spin degrees of freedom equals half the charge Mott-Hubbard gap [10].
sites refers to rotated-electron unoccupied or doubly occupied sites.Obviously, in 2D there are many other shapes for compact domains of The square shape will be justified below.In this configuration and considering processes of energy lower than 2|µ|, either the Pauli principle or the onsite repulsion prevents the rotated electrons to move (except for those on the lines referring to the square edges separating the two 2D domains).Although in this configuration there is no direct interaction between the rotated electrons, through an intermediate state with energy 2|µ| neighboring rotated electrons can see each other's spins, and rotated electrons with different spins can change position.Alike for 1D, the situation with the rotated-electron unoccupied and doubly occupied sites is similar for a original-lattice compact square domain of Unfortunately and unlike for 1D, for a finite system the corresponding spinon and η-spinon distributions may change if the square-shape compact domain of rotated-electron singly occupied sites is "diluted" by rotated-electron unoccupied and doubly occupied sites.Indeed, for a 2D system there is no order equivalent to the 1D uniquely defined chain order.
However, within the N 2 a → ∞ limit that the present description refers to the dominant c fermion occupancy configurations refer to an average uniform distribution of the c effective lattice occupied sites.In contrast to a finite system, configurations where for instance the c fermion occupied and unoccupied sites refer to two compact domains have vanishing weight as N 2 a → ∞.Specifically, the c fermion positions of such a uniform configuration correspond to the average positions of the c fermions in an energy eigenstate.Indeed, the c fermion momentum occupancy configurations of such a state are a superposition of all compatible real-space c effective lattice occupancy configurations.Fortunately, for the average configuration in which the rotated-electron singly occupied sites are uniformly "diluted" by rotated-electron unoccupied and doubly occupied sites the spinon and η-spinon distributions of the above original-lattice square-shape compact domain of N 2 as = 2S c = [2S s + M co s ] = N c sites do not change.They may be described by corresponding occupancy configurations of an effective square lattice with N 2 as = 2S c = [2S s + M co s ] = N c sites and edge length L whose average spacing is for x = 0 larger than that of the original lattice.The same arguments apply to the η-spinon distributions.This reveals that for the Hubbard model on the square lattice with N 2 a → ∞ sites the concepts of a spin effective square lattice and η-spin square effective lattice apply provided that n = (1 − x) and x are finite, respectively.
That the above considered initial reference original-lattice compact domain of ) sites has a square shape follows from only that compact domain shape transforming into a corresponding uniform "diluted" domain fully contained in the original square lattice.Any other initial reference compact-domain shape would transform into a uniform "diluted" domain that is not fully contained in the original square lattice and thus is neither physically nor mathematically acceptable.Here we are considering that the 2D crystal has a square shape whose edge is L = a N a .If the square-lattice crystal has any other shape the compact domain should have a similar shape with an area reduced by a factor of 2S c /N 2 a for the spin effective lattice (and of [N 2 a − 2S c ]/N 2 a for the η-spin effective lattice).Only then it transforms into a corresponding uniform "diluted" domain fully contained in the crystal.The final result is though the same, the obtained effective lattice being a square lattice whose spacing is given by, where N aα = (N D aα ) 1/D .This spinon and η-spinon distributions "average invariance" under "site dilution" emerging for the square lattice in the N D a → ∞ limit is behind the description of such distributions in terms of the occupancy configurations of independent η-spin and spin effective square lattices, respectively, and the direct propagation of c fermions in terms of occupancy configurations in an independent c effective lattice.While for 0 < x < 1 and 0 < m < n the former two lattices have both for 1D and 2D a number of sites smaller than N D a , the c effective lattice is identical to the original lattice.Within the N D a 1 limit of the present description the concept of a spin (and η-spin) effective lattice is well defined for finite values of the electronic density n = (1 − x) (and hole concentration x).The reasoning for the validity of the use of the corresponding effective lattices occupancy configurations may be summarized by the two following statements: (1) The representation associated with the present description contains full information about the relative positions of the sites of the η-spin and spin effective lattices in the original lattice.For each energy-eigenstate rotated-electron real-space occupancy configuration, that information is stored in the corresponding occupancy configurations of the c fermions in their c effective lattice.The latter lattice is identical to the original lattice.Such configurations correspond to the state representations of the U (1) symmetry in the subspaces spanned by states with fixed values of S c , S η , and S s .Indeed, the sites of the η-spin (and spin) effective lattice have in the original lattice the same real-space coordinates as the sites of the c effective lattice unoccupied (and occupied) by c fermions.(2) Within the N D a 1 limit that the rotated-electron related operational description refers to, provided that the electronic density n = (1 − x) (and hole concentration x) is finite, the dominant c effective lattice occupancy configurations of an energy eigenstate of the Hubbard model on the square lattice refer to a nearly uniform distribution of the c fermions occupied sites (and unoccupied sites).Hence due to the spinon and η spinon distribution "average order" emerging for the model on the square lattice in the N D a → ∞ limit, the spin and η-spin effective lattices may be represented by square lattices.Moreover, the chain order invariance occurring for the 1D model both for the finite system and in that limit justifies why such effective lattices are 1D lattices.For both models the corresponding spin effective lattice spacing a s and η-spin effective lattice spacing a η refers to the average spacing between the c effective lattice occupied sites and between such a lattice unoccupied sites, respectively, given in Equation (30).Such spin and η-spin effective lattices obey the physical requirement condition that in the x → 0 and x → ±1 limit, respectively, equal the original lattice.Note that in the x → 0 (and x → ±1) limit one has that N D as = N D a and the η-spin effective lattice does not exist (and N D aη = N D a and the spin effective lattice does not exist.) The validity for N D a 1 of the concept of a spin effective square lattice as constructed in this paper is confirmed by the behavior of the expectation value δd = Ψ|δ d|Ψ of any energy eigenstate |Ψ .Here δ d is the operator associated with the distance in real space of any of the N D as sites of the spin effective lattice from the rotated-electron singly occupied site of the original lattice closest to it.The state |Ψ belongs to a subspace with fixed number of rotated-electron singly occupied sites.The point is that δd = Ψ|δ d|Ψ vanishes in the N D a → ∞ limit.We recall that the N D as = 2S c sites of the spin effective lattice occupancies describe the spin degrees of freedom of the N c = 2S c rotated-electron singly occupied sites of the original lattice.The same applies to the η-spin effective lattice.
Consistently with the expression a α = L/N aα = [N a /N aα ] a of Equation ( 30) where α = η, s, the η-spin (and spin) effective lattice has both for 1D and 2D the same length and edge length L, respectively, as the original lattice.Furthermore, the requirement that for the 2D case when going through the whole crystal of square shape along the ox 1 or ox 2 directions a η-spinon (and spinon) passes an overall distance L is met by an effective η-spin (and spin) square lattice.Since the number of sites sum-rule [N D aη + N D as ] = N D a holds, the η-spin and spin effective lattices have in general a number of sites N D aη and N D as , respectively, smaller than that of the original lattice, N D a .It follows that their lattice spacings Equation (30) are larger than that of the original lattice.

The Composite αν Fermions
The global SO(3) × SO(3) × U (1) symmetry of the model on any bipartite lattice implies that some features of the corresponding state representations are common to all such lattices.Our extension of the 1D model -site η-spin-neutral rotated-electron occupancy configurations in terms of those of composite sν fermions and ην fermions, respectively, to the model on the square lattice accounts for the basic differences between the physics of the two models.This holds as well for the related momentum eigenstates |Φ U/4t = V † |Φ ∞ considered in Reference [9].For the model on the square lattice, the latter states refer in general to a partition different from that of the energy eigenstates |Ψ U/4t = V † |Ψ ∞ of the η-spin-neutral (and spin-neutral) configurations of such M co η (and M co s ) η-spinons (and spinons), in terms of smaller configurations.Specifically, for a given momentum eigenstate there is for each branch involving ν = 1, 2, ... pairs of η-spinons (and spinons) a well-defined number N ην of η-spin-neutral 2ν-η-spinon composite ην fermions (and N sν of spin-neutral 2ν-spinon composite sν fermions).The set of such composite objects describes the η-spinon (and spinon) occupancy configurations of exactly M co η (and M co s ) sites of the original lattice.Hence the following two sum rules hold for all momentum eigenstates, Such sum rules refer to subspaces spanned by states with fixed values of S c , S η , and S s .One η-spin neutral 2ν-η-spinon composite ην fermion describes the η-spin degrees of freedom of a η-spin-singlet occupancy configuration involving ν ≤ M co η /2 sites of the original lattice unoccupied by rotated electrons and an equal number of sites doubly occupied by rotated electrons.The remaining degrees of freedom of such a rotated-electron occupancy configuration are described by 2ν unoccupied sites of the c effective lattice whose spatial coordinates are those of the corresponding 2ν sites of the original lattice.
Similarly, one spin neutral 2ν-spinon composite sν fermion describes the spin degrees of freedom of a spin-singlet occupancy configuration involving ν ≤ M co s /2 sites of the original lattice singly occupied by spin-up rotated electrons and an equal number of sites singly occupied by spin-down rotated electrons.The remaining degrees of freedom of that rotated-electron occupancy configuration are described by 2ν occupied sites of the c effective lattice whose spatial coordinates are those of the corresponding 2ν sites of the original lattice.
For each ην fermion branch (and sν fermion branch), one may consider a ην effective lattice (and sν effective lattice).It refers to occupancy configurations of η-spin-neutral (and spin-neutral) bonds of 2ν confined η-spinons (and spinons).Hence each "occupied site" of such an effective lattice corresponds to 2ν sites of the η-spin (and spin) effective lattice.In turn, the M de η = 2S η (and M de s = 2S s ) η-spin (and spin) effective lattice sites referring to the deconfined η-spinons (and deconfined spinons) and some of such a lattice sites referring to 2ν -η-spinon composite ην fermions (and 2ν -spinon composite sν fermions) of ν > ν branches are found below to play the role of "unoccupied sites" of such a ην effective lattice (and sν effective lattice).The conjugate of the αν effective lattice site space variables are the αν band discrete momentum values.Their number equals that of the αν effective lattice sites.(For 1D, such a momentum values are good quantum numbers.) Within chromodynamics the quarks have color but all quark-composite physical particles are color-neutral [21].Here the η-spinon (and spinons) that are not invariant under the electron-rotatedelectron unitary transformation have η-spin 1/2 (and spin 1/2) but are confined within η-spin-neutral (and spin-neutral) 2ν-η-spinon (and 2ν-spinon) composite ην fermions (and sν fermions).The exact and detailed internal 2ν-spinon configuration and 2ν-η-spinon configuration of a composite sν fermion and ην fermion, respectively, is in general an involved unsolved problem.(In 1D the exact BA solution takes implicitly into account such internal configurations.)Fortunately, however, the problem simplifies for the model in the one-and two-electron subspace reviewed below for which the only composite object that plays an active role is the two-spinon s1 fermion.Its internal structure is an issue studied in Reference [9].Such a two-spinon object is related to the resonating-valence-bond pictures for spin-singlet occupancy configurations of ground states [22,23].
In the following we consider all multi-η-spinon and multi-spinon composite objects.Full information on the 2ν-η-spinon (α = η) or 2ν-spinon (α = s) configurations associated with the internal degrees of freedom of the composite αν fermions is not needed for the goals of this paper.Indeed, within the present N D a 1 limit the problem of the internal degrees of freedom of the composite αν fermions and αν bond particles separates from that of their positions in the corresponding effective lattices.The partial information on the internal degrees of freedom of the composite αν fermions needed for the studies reported in this paper is accessed in the following by suitable use of their transformation laws under the electron-rotated-electron unitary transformation.In turn, the deconfined η-spinons and deconfined spinons are invariant under that transformation.Thus they are non-interacting deconfined objects that are not part of composite ην fermions and composite sν fermions, respectively.

The αν Fermion Operators
For the αν effective lattice, one local αν fermion "occupied site" refers to 2ν sites of the η-spin (α = η) or spin (α = s) effective lattice.Alike for the corresponding latter lattice, the 2S α sites occupied by deconfined η-spinons (α = η) or deconfined spinons (α = s) are among those playing the role of the αν effective lattice "unoccupied sites".Unlike for the former lattice, it is found below that a number 2(ν − ν) of sites of each αν fermion with a number ν > ν of confined η-spinons (α = η) or confined spinons (α = s) play as well the role of αν effective lattice "unoccupied sites".
The conjugate variables of the αν effective lattice real-space coordinates are the discrete momentum values of the αν band.As shown in Reference [9], for the 1D model such discrete momentum values are the quantum numbers of the exact BA solution.For the Hubbard model on the square lattice the s1 band momentum discrete values of state representations belonging to the one-and two-electron subspace reviewed below are good quantum numbers as well.
As reported in Reference [10] for the s1 fermion operators, the αν fermion operators can be generated from the operators of corresponding hard-core αν bond-particle operators as follows, f † r j ,αν = e iφ j,αν g † r j ,αν ; φ j,αν = j =j f † r j ,αν f r j ,αν φ j ,j,αν ; φ j ,j,αν = arctan Such expressions are valid provided that (1 − x) > 0 for the αν = s1 branch and S α /N D a > 0 for the remaining αν = s1 branches.Explicit expressions for the operators f † r j ,s1 and g † r j ,s1 are given in Reference [9].The quantity φ j,αν appearing in Equation ( 32) is the Jordan-Wigner phase [10,24] operator, the indices j and j refer to sites of the αν effective lattice, and f † q j ,αν are the corresponding momentum-dependent αν fermion operators.The number N D aαν of discrete momentum values of the αν momentum band equals that of sites of the αν effective lattice.Its expression is derived below.
The η-spin-neutral 2ν-η-spinon composite ην bond-particle operators and spin-neutral 2ν-spinon composite sν bond-particle operators denoted in Equation ( 32) by g † r j ,αν where α = η, s are constructed to inherently upon acting onto their αν effective lattice anticommuting on the same site and commuting on different sites.Hence they are hard-core like and can be transformed into fermionic operators, as given in that equation.For ν > 1 the algebra behind their construction in terms of the elementary η-spinon or spinon operators of Equations ( 12) and ( 13) is much more cumbersome than that of the two-spinon s1 bond particles studied in Reference [9].Fortunately, the only property needed for the goals of this paper is that upon acting onto their αν effective lattice they are hard-core like.
The expressions given in Equation ( 32) apply to αν = s1 branches provided that S α /N D a > 0. One can also handle the problem when S α = 0 and N αν /N D a 1 for a given αν = s1 branch.Then provided that N αν = 0 for all remaining αν branches with a number of η-spinon or spinon pairs ν > ν one finds below that N D aαν = N αν .The momentum of the operators f † q j ,αν is in that limiting case given by q j ≈ 0. For the states that span the corresponding S α = 0 and N αν /N D a 1 subspace all sites of the αν effective lattice are occupied and the αν momentum band is full.If N αν is finite one has N D aαν = N αν discrete momentum values q j ≈ 0 compactly distributed around zero momentum.Their Cartesian components momentum spacing is 2π/L.A case of interest is when S α = 0 and N αν = 1.Then the αν effective lattice has a single site and the corresponding αν band a single discrete momentum value, q = 0.In that case φ αν = φ j,αν = 0. Hence f † r,αν = g † r,αν and f † q,αν = f † r,αν where q = 0.The operators f † q j ,αν act onto subspaces with fixed values for the set of numbers S α , N αν , and {N αν } for ν > ν branches.(Below it is shown that this is equivalent to fixed values for the set of numbers S c and {N αν } for all ν = 1, 2, . . .including ν.)Such subspaces are spanned by mutually neutral states.Those are states with fixed values for the numbers of αν fermions and αν fermion holes.Hence such states can be transformed into each other by αν band particle-hole processes.Creation of one αν fermion is a process that involves the transition between two states belonging to different such subspaces.It is a well-defined process whose generator is the product of two operators.The first operator may add sites to or remove sites from the αν effective lattice.Alternatively, it may introduce corresponding changes in the αν momentum band.The second operator is the creation operator f † r,αν or f † q,αν appropriate to the excited-state subspace.
Provided that (1−x) is finite for s1 fermions [10] and S α /N D a > 0 as N D a → ∞ for αν = s1 fermions, the phases φ j,αν given in Equation ( 32) are associated with an effective vector potential [24,25], For the model on the square lattice the vector e x 3 appearing here is the unit vector perpendicular to the plane.(Often we use units such that the fictitious magnetic flux quantum is given by Φ 0 = 1.) The components of the microscopic momenta of the αν fermions are eigenvalues of the two (and one for 1D) αν translation generators qαν x 1 and qαν x 2 in the presence of the fictitious magnetic field B αν ( r j ).That seems to imply that for the model on the square lattice the components q x1 and q x2 of the microscopic momenta q = [q x1 , q x2 ] refer to operators that do not commute.However, in the subspaces where the operators f † q j ,αν act onto such are commuting operators.Indeed, those subspaces are spanned by neutral states [25].Since [q αν x 1 , qαν x 2 ] = 0 in such subspaces, for the model on the square lattice the αν fermions carry a microscopic momentum q = [q x1 , q x2 ] where the components q x1 and q x2 are well-defined simultaneously.

Ranges of the c and αν Fermion Energies, Their Transformation Laws, and the Ground-State Occupancies
Since the microscopic momenta of the c fermions are good quantum numbers both for the Hubbard model on the square and 1D lattices, one may define an energy dispersion c ( q) [10].For the 1D model also the αν energy dispersions αν (q) are well defined.In turn, for the model on the square lattice an energy dispersion αν ( q) is for αν = s1 branches well defined for the momentum values for which such objects are invariant under the electron-rotated-electron unitary transformation.For general αν fermion momentum values there is not in general a dispersion αν = αν ( q) defining a one-to-one correspondence between the energy αν and momentum q.However, the range of the αν fermion energy αν remains well defined.A s1 fermion energy dispersion s1 = s1 ( q) is well defined for the square-lattice quantum liquid [10].Indeed for the Hubbard model in the one-and two-electron subspace the s1 band momentum values become good quantum numbers [9].
The quantum-object occupancy configurations of the ground state are found below.They are consistent with the energies for creation onto such states of our description quantum objects.For instance, the elementary energies s,−1/2 = 2µ B H and η,−1/2 = 2µ correspond to creation onto a m ≥ 0 and x > 0 ground state of a deconfined −1/2 spinon and a deconfined −1/2 η-spinon, respectively.Here µ B is the Bohr magneton and as above H and µ are the magnetic field and the chemical potential, respectively.The energy s,−1/2 = 2µ B H (and η,−1/2 = 2µ) refers to an elementary spin-flip (and η-spin-flip) process.It transforms a deconfined +1/2 spinon (and deconfined +1/2 η-spinon) into a deconfined −1/2 spinon (and deconfined −1/2 η-spinon).Such elementary energies control the range of several physically important energy scales.Within the LWS representation, a deconfined +1/2 spinon (and deconfined +1/2 η-spinon) has vanishing energy so that s,+1/2 = 0 (and η,+1/2 = 0).It follows that the energy of a pair of deconfined spinons (and deconfined η-spinons) with opposite projections is 2µ B H (and 2µ).Indeed due to the invariance of such objects under the electron-rotated-electron unitary transformation, they are not energy entangled and the total energy is the sum of their individual energies.
In the following we confirm that ground states have no ην fermions and no sν fermions with ν > 1 spinon pairs.The corresponding energies ην and sν , respectively, considered below refer to creation onto the ground state of one of such objects.We start by providing a set of useful properties.We emphasize that some of these properties are not valid for descriptions generated by rotated-electron operators associated with the general unitary operators V considered in Reference [6].The following properties rather refer to the specific operator description associated with the rotated-electron operators c † r j ,σ = V † c † r j ,σ V of Equation (3) as defined above.
Some of the following results are obtained from extension to the model on the square lattice of exact results extracted from the 1D model BA solution.However such an extension accounts for the different physics of such models.Both for the model on the square and 1D lattice, the range of the energy αν for addition onto the ground state of one αν fermion derived below is that consistent with the interplay of the transformation laws of the αν fermions under the electron-rotated-electron unitary transformation with the model global SO(3) × SO(3) × U (1) symmetry.

The ην Fermion Energy Range
Alike for 1D, for the model on the square lattice one ην fermion is a η-spin-neutral anti-bounding configuration of a number ν = 1, 2, ... of confined −1/2 η-spinons and an equal number of confined +1/2 η-spinons.Symmetry implies that for U/4t > 0 there is no energy overlap between the ην ranges corresponding to different ν = 1, 2, ... branches.Fermions belonging to neighboring ην and ην + 1 branches differ in the number of η-spinon pairs by one.The requirement for the above lack of energy overlap is then that the energy bandwidth of the ην range is smaller than or equal to 2|µ|.For all x values, the energy scale 2|µ| = [ η,−1/2 + η,+1/2 ] where µ = µ 0 at x = 0 refers to the energy of a pair of deconfined η-spinons of opposite η-spin projection.Such properties imply the following range for the energy ην , 2ν|µ| where µ = µ 0 at x = 0. Deconfined η-spinons are invariant under the electron-rotated-electron unitary transformation.Consistently and as justified below, they are non-interacting and their energies are additive.For all x values 2ν|µ| is the energy of ν deconfined −1/2 η-spinons and ν deconfined +1/2 η-spinons.For instance, for x > 0 the energy 2νµ is as well that for creation of a number ν = 1, 2, . . . of deconfined −1/2 η-spinons onto a S η = ν ground state with 2ν deconfined +1/2 η-spinons.Such a creation refers to ν η-spin-flip processes (transformation of ν rotated-electron unoccupied sites into ν rotated-electron doubly occupied sites.)For m = 0 and x > 0 the number i ην decreases continuously for increasing values of U/4t.It has the limiting behaviors i ην → 1 for U/4t → 0 and i ην → 0 for U/4t → ∞.Hence the ην range vanishes for U/4t → ∞.The latter behavior is associated with the full degeneracy of the η-spin configurations reached as U/4t → ∞.In turn, at m = 0 and x = 0 the number i ην vanishes and ην = 2νµ 0 for the whole finite interaction range U/4t > 0. As discussed below, this behavior follows from the invariance under the electron-rotated-electron unitary transformation of a ην fermion created onto a x = 0 and m = 0 ground state.
For the 1D Hubbard model one can define momentum dependent ην fermion energy bands [26].The η1 and η2 fermion bands are plotted in Figures 1 and 2, respectively, in units of t for electronic densities n = (1 − x) = 1/2, 5/6, spin density m = 0, and different U/t values.In the figures the lower energy level of the bands was shifted to zero.According to the above analysis, for U/t → 0 and U/t → ∞ the energy bandwidth of all ην bands should read 2|µ| and 0, respectively.For electronic densities 1/2 and 5/6 and U/t → 0 one has that twice the chemical-potential absolute value reads 2|µ| = 2 √ 2 t ≈ 2.828 t and 2|µ| = 4 t cos(5π/12) ≈ 1.035 t, respectively.Consistently, notice that in the U/t → 0 limit the energy bandwidths of the η1 and η2 fermion bands plotted in the figures are 2 √ 2 t ≈ 2.828 t and 4 t cos(5π/12) ≈ 1.035 t for electronic densities 1/2 and 5/6, respectively, Furthermore, in the limit U/t → ∞ the energy bandwidths of the η1 and η2 fermion bands indeed vanish.The parameter i ην appearing in Equation (34), which controls such energy bandwidths, has the same limiting behaviors i ην → 1 for U/t → 0 and i ην → 0 for U/t → ∞ for all remaining ην bands corresponding to a number ν = 1, 2, 3, . . . of confined η-spinon pairs.Figure 1.The η1 fermion band for the 1D Hubbard model in units of t plotted for electronic densities n = (1 − x) = 1/2, 5/6, spin density m = 0, and a set of U/t values.The band U/t → ∞ limit corresponds to the horizontal line chosen as zero-energy level.[26] -0.5 0 0.5 q/π 0 0.5 The sν Fermion Energy Range A sν fermion is a spin-neutral bounding configuration of a number ν = 1, 2, . . . of confined −1/2 spinons and an equal number of confined +1/2 spinons.Again, symmetry implies that for U/4t > 0 there is no energy overlap between the sν ranges of different ν = 1, 2, . . .branches.For sν branches with a number of spinon pairs ν > 1 such an energy range bandwidth is for the present bounding configurations and for the same reasoning as for the ην fermion smaller than or equal to 2µ B |H|.For all m values, 2µ B |H| = [ s,−1/2 + s,+1/2 ] equals the energy of a pair of deconfined spinons of opposite spin projection.Hence the range of the energy sν for addition onto the ground state of one sν fermion with ν > 1 spinon pairs is, Deconfined spinons are invariant under the electron-rotated-electron unitary transformation.
Since for a non fully polarized state 2νµ B H → 0 as U/4t → ∞, the energy bandwidth of the sν range vanishes for U/4t → ∞.Such a behavior is associated with the full degeneracy of the spin configurations reached for U/4t → ∞. 1 and 2 for the η1 and η2 fermion bands, respectively, confirm such a i sν behavior for the sν bands associated with ν > 1 confined spinon pairs, but with the U/4t → 0 energy bandwidth 2|µ| replaced by 2µ B |H| and the electronic densities n = (1 − x) < 1 replaced by spin densities m > 0. Note that at m = 0 and thus H = 0 one has that sν = 0 for the whole finite interaction range U/4t > 0. This behavior follows from the invariance under the electron-rotated-electron unitary transformation of a sν fermion with ν > 1 spinon pairs created onto a m = 0 ground state.

Figures similar to those of Figures
In turn, it is found below that for a m = 0 and x ≥ 0 ground state all sites of the s1 effective lattice are occupied.Hence the corresponding s1 momentum band is full.The range of the energy − s1 for removal from that state of one s1 fermion then is, Here |∆| vanishes for 1D.The s1 fermion energy band is plotted in units of t in Figure 3 for the 1D Hubbard model at electronic densities n = (1 − x) = 1/2 and n = 5/6, spin density m = 0, and a set of U/t values.For the Hubbard model on the square lattice, |∆| denotes the s1 fermion pairing energy per spinon considered in Reference [10].There is strong evidence of the occurrence in the half-filling Hubbard model on the square lattice of strong short-range antiferromagnetic correlations for T > 0 and below a crossover temperature called T x in Reference [20].This is consistent with then the system being driven into a phase with short-range spin order.Within the representation of Reference [9], the occurrence of long-range antiferromagnetic order requires that T = 0, N 2 aη = 0, and N 2 as = N 2 a .This is then consistent with the short-range spin order occurring for m = 0, 0 < x 1, and 0 ≤ T < T x being similar to that occurring for m = 0, x = 0, and 0 < T < T x .The latter order was studied previously in Reference [20] for 0 < T T x .For small finite hole concentrations 0 < x 1 the s1 fermion pairing energy 2|∆| is identified with 2∆ 0 ≈ 2k B T x .In turn, for U/4t > 0 and x = 0 its magnitude 2|∆| = µ 0 is larger than for x → 0 and the s1 fermion spinon pairing refers to an antiferromagnetic long-range order [10].
For 1D there is no spin short-range spin order, so that 2|∆| = 0.The temperature T x of Reference [20] is plotted in Figure 3 of that reference.Its U/4t dependence is qualitatively correct.T x vanishes both in the limits U/4t → 0 and U/4t → ∞.It goes through a maximum magnitude at an intermediate value 5/4 < U/4t < 3/2.Nevertheless, the interpolation function used to produce it, provided in Ref. 74 of such a paper, is poor for intermediate values of U/4t.That used in the studies of Reference [10] is such that 2∆ 0 ≈ 2k B T x vanishes both in the U/4t → 0 and U/4t → ∞ limits and goes through a maximum value max {2∆ 0 } ≈ 2t/π at U/4t ≈ 1.3.At fixed U/4t values, 2|∆| decreases for increasing x as 2|∆| ≈ 2∆ 0 (1 − x/x * ) and vanishes for x > x * .Here x * ∈ (0.23, 0.28) for U/4t ∈ (1.3, 1.6) is a critical hole concentration below which the s1 fermion pairing refers to a spin short-range order [10].
On the other hand, at m = 0 the energy scale W s1 is the s1 fermion energy nodal bandwidth defined in Reference [10].Its maximum magnitude is reached at U/4t = 0.For U/4t > 0 it decreases monotonously for increasing values of U/4t, vanishing for U/4t → ∞.That W s1 → 0 for U/4t → ∞ is associated with the full degeneracy of the spin configurations reached in that limit.In it the spectrum of the two-spinon composite s1 fermions becomes dispersionless for the square-lattice quantum liquid of Reference [10].

The c Fermion Energy Range
The energy c for addition onto the ground state of one c fermion of a given momentum and the energy − c for removal from that state of such a c fermion have the following ranges, respectively.Here W h c = [4Dt − W p c ] ∈ (0, 4Dt) increases monotonously for increasing values of hole concentration x ∈ (0, 1).The energy bandwidth W p c depends little on U/4t.For U/4t > 0 it has the following limiting behaviors, The behaviors reported here for c are justified in Reference [10].The c fermion energy band is plotted in units of t in Figure 4 for the 1D Hubbard model at electronic densities n = (1 − x) = 1/2, 5/6, spin density m = 0, and a set of U/t values.Note that for 1D (and the square lattice) its energy bandwidth is 4t (and 8t) for all density values and the whole U/4t range.The minimum magnitude of the energy ∆ Drot for creation of a number D rot = M η,−1/2 of rotated-electron doubly occupied sites onto a m = 0 and x ≥ 0 ground state given in Equation ( 26) may be expressed in terms of both the numbers of deconfined η-spinons and ην fermions as follows, energy sν = 0 (and ην = 2νµ 0 ).Furthermore, its momentum is q = q s = 0 (and q = q η = 0).Hence it obeys indeed to the condition given in Equation (42).In turn, sν fermions with ν > 1 spinon pairs (and ην fermions) created onto initial S s > 0 (and S η > 0) ground states may have energy sν smaller than 2νµ B |H| (and energy ην larger than 2ν|µ|), so that they are not necessarily invariant under the electron-rotated-electron unitary transformation.
For the 1D model the αν fermion energy αν (q) depends on the αν fermion momentum q, which is a good quantum number.For U/4t > 0 and the sν branches with ν > 1 spinon pairs and all ην branches such momenta belong to the range q ∈ (−mπ, +mπ) and q ∈ (−xπ, +xπ), respectively.Only at the limiting momenta q = q η = ±xπ (and q = q s = ±mπ) is the invariance condition given in Equation (42) met by the ην fermion energy ην (q) (and sν fermion energy sν (q) for ν > 1 branches).The ην fermion energy dispersions ην (q) are plotted for the η1 (and η2) branches in Figures 8(a) and 9(a) (and Figures 8(b) and 9(b)) of Reference [26] as a function of q for several U/4t values and electronic density n = 1/2 and n = 5/6, respectively.(In that reference the ην fermions are called c, ν pseudoparticles.) The zero-energy level of Figures 8(a For the Hubbard model on the square lattice the αν fermions whose energy obeys the invariance condition given in that equation have a well-defined momentum q α , which can point in different directions.At x = 0 (and m = 0) the momentum q η (and q s ) vanishes.This is alike for the above corresponding momentum q η = ±xπ (and q s = ±mπ) of the 1D model, which vanishes at x = 0 (and m = 0).While at 1D the momentum q η = ±xπ (and q s = ±mπ) can for x > 0 (and m > 0) have two values, for the square-lattice model the momentum q η (and q s ) can for x > 0 (and m > 0) point to several directions.For instance, for x < x * and m = 0 where x * is the hole concentration considered in Reference [10] below which the ground-state s1 fermion spinon pairing is associated with a short-range spin order, a good approximation for the momentum q η centered at − π = [−π, −π] is, Here q h F c is the c fermion hole Fermi momentum considered in Reference [10] and e φc is a unit vector centered at − π of Cartesian components [cos φ c , sin φ c ].As required, q η → 0 as x → 0.

Ground State Occupancies
Consistent with the energy values given in Equations ( 40) and (41), one finds that both for the model on the 1D and square lattice in the subspace spanned by the x > 0 and m > 0 LWS ground states and their excited energy eigenstates of energy ω < min {2µ, 2µ B H} the η-spinon and spinon numbers are given by Hence for such energy eigenstates the numbers of ην fermions and sν fermions with ν > 1 spinon pairs vanish and the number of c fermions is N c = 2S c = N and that of s1 fermions is conserved and reads N s1 = N ↓ .
Except that the −1/2 η-spinons and −1/2 spinons play the role of the +1/2 η-spinons and +1/2 spinons, respectively, and vice versa, similar results are reached for highest-weight states (HWSs) of the η-spin and spin algebras.Comparison of the occupancies of the spin LWS ground states (m > 0) and spin HWS ground states (m < 0) provides useful information.From it and again consistently with the energy values of Equations ( 40) and (41), one finds that a m = 0 ground state for which N is even and x ≥ 0 has M η,±1/2 and N c values as given above whereas M s,±1/2 = M co s /2 = N/2 so that M de s,±1/2 = 0 and N s1 = N/2.

The Site Numbers and Spacing of the αν Effective Lattices
The number of sites of the αν effective lattice, which equals that of αν band discrete momentum values, has the general form, where N h αν is the number of unoccupied sites whose expression is derived in Reference [9] on use of symmetry representation counting and reads, where M de α = 2S α .Importantly, this expression is also that obtained from the BA exact solution for the 1D Hubbard model.However, it is fully determined by state-representation dimension requirements of the global SO(3) × SO(3) × U (1) symmetry that apply to the Hubbard model on the square lattice as well.(Note that the equivalence of the two expressions given in Equation ( 45) confirms that the N c = 2S c and {N αν } values with ν = 1, 2, . . .remaining fixed is equivalent to the η-spin S η (α = η) or spin S s (α = s) and values of the set of numbers {N αν } with ν = ν, ν + 1, . . .remaining fixed as well.In both cases that implies that the N D aαν value remains fixed.)Provided that N D aαν /N D a remains finite as N D a → ∞, the related αν effective lattices can for the 1D and square-lattice models be represented by 1D and square lattices, respectively, of spacing, Here ν = 1, 2, . . .and α = η, s.The arguments behind the lattice geometry and the average distance a αν between the sites of the αν effective lattice playing the role of lattice spacing are similar to those used for the lattice geometry and spacing of the η-spin and spin effective lattices.In turn, the corresponding αν bands whose number of discrete momentum values is also given by N D aαν are well defined even when N D aαν is given by a finite small number, N D aαν = 1, 2, . . .In this section we review a suitable one-and two-electron subspace.(For hole concentrations x > 0 our studies refer to excitation energies below 2µ for which M η,−1/2 = 0 for that subspace.)For x ≥ 0 the picture that emerges is that of a two-component quantum liquid of charge c fermions and spin neutral two-spinon s1 fermions.It refers to the square-lattice quantum liquid further investigated in Reference [10].Moreover, a preliminary application of the description reviewed here is presented.

The One-and Two-Electron Subspace
General N -Electron Subspaces We consider a x ≥ 0 and m = 0 ground state |Ψ GS .Application onto it of a N -electron operator ÔN generates a state, contained in a N -electron subspace.This is a subspace spanned by the set of energy eigenstates It is thus associated with a given N -electron operator ÔN .Generally, such operators can be written as a product of one-electron creation and annihilation operators.Here N lη,ls is the number of electron creation and annihilation operators in the operator ÔN expression for l η = −1 and l η = +1, respectively, and with spin down and spin up for l s = −1 and l s = +1, respectively.A general local N -electron operator ÔN,j refers to a product of N local electron creation and annihilation operators.For N > 1 such an operator has a well defined local structure.It involves N −1,ls electron creation operators of spin projection l s /2 and N l +1,ls electron annihilation operators of spin projection l s /2 whose real-space coordinates refer in general to a compact domain of neighboring lattice sites.Such a local N -electron operator ÔN,j may be labelled by the real-space coordinate r j of a corresponding central site.A second type of N -electron operator is denoted by ÔN ( k) and carries momentum k.It is related to a local operator ÔN,j by a Fourier transform.
The general N -electron operators ÔN considered here belong to one of these two types and are such that the ratio N /N D a vanishes in the thermodynamic limit.The operators ÔN ( k) of physical interest correspond in general to operators ÔN,j whose N elementary electronic operators create or annihilate electrons in a compact domain of lattice sites.The more usual cases for the description of experimental studies correspond to the N = 1 one-electron and N = 2 two-electron operators.Therefore, in this section we are mostly interested in the corresponding one-and two-electron subspace.
Application onto a x ≥ 0 and m ≥ 0 ground state of a general N -electron operator ÔN leads to electron number deviations δN = δN ↑ +δN ↓ and δN ↑ −δN ↓ .As a result of the expressions and relations given above in this paper, such deviations may be expressed in terms of corresponding deviations in the number of c fermions, αν fermions, deconfined η-spinons, and deconfined spinons as follows, and respectively.(Note that M de η,−1/2 = M co η = N ην = M de s,−1/2 = N sν | ν>1 = 0 for the initial ground state so that δM de η,−1/2 = M de η,−1/2 , δM co η = M co η , δN ην = N ην , δM de s,−1/2 = M de s,−1/2 , and δN sν | ν>1 = N sν | ν>1 .)Hence only transitions to excited states associated with deviations obeying the sum rules Equations ( 49) and (50) are permitted.The electron number deviations Equations ( 49) and (50) are associated with sum rules obeyed by the numbers N lη,ls of Equation (48) specific to the N -electron operator ÔN under consideration, which read, Furthermore, it is straightforward to show that useful exact selection rules hold for excitations of well-defined initial ground states.For instance, the values of the numbers M de η = [M de η,−1/2 + M de η,+1/2 ] = 2S η of deconfined η-spinons and M de η,±1/2 of ±1/2 deconfined η-spinons generated by application onto a S η = S s = 0 ground state of a N -electron operator ÔN are restricted to the following ranges, respectively.Here S η denotes the excited-state η-spin and the numbers N lη,ls are those of Equation ( 48) specific to the N -electron operator.
In the case of an initial S s = 0 ground state with hole concentration x > N /N D a or x = 0 one finds that the numbers M de s = [M de s,−1/2 + M de s,+1/2 ] = 2S s and M de s,±1/2 of deconfined spinons generated by application onto that state of a N -electron operator are restricted to the ranges, The range restrictions of Equations ( 52) and ( 53) are exact for both the model on the square and 1D lattice, as well as for any other bipartite lattice.For x > 0 we limit our considerations to the vanishing rotated-electron double occupancy subspace.Consistently with the ∆ Drot energy spectrum of Equations ( 26) and (39), this is accomplished merely by limiting the excitation energy to values below 2µ, so that the M η,−1/2 = 0 constraint is automatically fulfilled: It follows from the form of such a spectrum that excited states with ∆ Drot < 2µ have vanishing rotated-electron double occupancy.Indeed creation of one rotated-electron doubly occupied site onto an initial x > 0 and m = 0 ground state is a process of minimum energy 2µ.In turn, at x = 0 we consider both states with vanishing and finite rotated-electron double occupancy.
Creation onto the S η = 0, µ = 0, and S s = 0 ground state of one ην fermion is a vanishing momentum process whose finite energy is exactly given by ην = 2νµ 0 .That object then obeys the criterion of Equation (42) for invariance under the electron-rotated-electron unitary transformation.It follows that the η-spin degrees of freedom of such a ην fermion exactly involve ν electron doubly occupied sites.Furthermore, creation onto an initial x ≥ 0 and S s = 0 ground state of one sν fermion with a number ν > 1 of spinon pairs is a vanishing energy and momentum process.Since vanishing spin S s = 0 refers to vanishing magnetic field H = 0, such an object obeys the criterion sν = 2νµ B |H| = 0 of Equation (42).Thus it is invariant under the electron-rotated-electron unitary transformation.It follows that for U/4t > 0 creation of such an object involves occupancy configurations whose spin degrees of freedom are similar in terms of both rotated-electron and electron occupancy configurations.That reveals that such a sν fermion describes the spin degrees of freedom of a number 2ν of electrons.
It then follows from the invariance under the electron-rotated-electron unitary transformation of the η-spin and spin degrees of freedom of the above ην fermion and sν fermion, respectively, that for x > N /N D a where N /N D a → 0 as N D a → ∞ and excitation energy ω < 2µ and any excitation energy at x = 0 nearly the whole spectral weight generated by application onto the above ground states of N -electron operators refers to a subspace spanned by excited states with numbers in the following range, The selection rules given in Equations ( 52) and (53) are exact.In turn, for x > 0 and m = 0 initial ground states and excitation energy ω < 2µ (and an initial x = 0 and m = 0 ground state and any excitation energy) nearly the whole N -electron spectral weight is generated by excited states whose numbers obey the approximate selection rules given in Equation (54).Indeed excited states with numbers [M co − 2B s ] > N for x > 0 (and M η > N and [M s − 2B s > ls=±1 N +1,ls ] at x = 0) generate a very small amount yet non-vanishing N -electron spectral weight.
Why, in spite of the invariance under the electron-rotated-electron unitary transformation of the sν fermions with a number ν > 1 of spinon pairs created onto an initial x ≥ 0 and m = 0 ground state (and that of the ην fermions created onto an initial x = 0 and m = 0 ground state), are the selection rules provided in Equation ( 54) not exact?The reason is that while the spin degrees of freedom of the 2ν-electron occupancy configurations involved in a sν fermion are exactly described by that object, their hidden U (1) symmetry degrees of freedom are not invariant under the electron-rotated-electron unitary transformation.The same applies to the η-spin degrees of freedom of the 2ν-electron occupancy configurations exactly described by a ην fermion created onto a x = 0 ground state.Their corresponding hidden U (1) symmetry degrees of freedom are not in general invariant under that transformation.The spin degrees of freedom of the 2ν-electron occupancy configurations involved in a sν fermion are for U/4t > 0 exactly the same as those of the corresponding 2ν-rotated-electron occupancy configurations.In turn, the occupancy configurations of the 2ν c fermions that describe the hidden U (1) symmetry degrees of freedom of the 2ν rotated electrons under consideration are slightly different from those of the corresponding 2ν electrons.The same applies to the η-spin degrees of freedom and hidden U (1) symmetry degrees of freedom of the 2ν electrons and corresponding 2ν rotated electrons involved in a ην fermion created onto a x = 0 ground state.The former and the latter are and are not invariant under that transformation.Hence that the selection rules of Equation (54) are not exact but a very good approximation stems from the lack of invariance under the electron-rotated-electron unitary transformation of the degrees of freedom associated with the hidden U (1) symmetry of the Hubbard model.This applies both to the model on the 1D and square lattice.
The N -electron spectral weight generated by excited states of initial x > 0 and m = 0 ground states of excitation energy ω < 2µ and numbers [M s − 2B s ] > N is extremely small.The same applies to excited states of initial x = 0 and m = 0 ground states of numbers M η > N and [M s − 2B s ] > ls=±1 N +1,ls .Therefore, in this paper we define the N -electron subspace as that spanned by an initial x ≥ 0 and m = 0 ground state plus the set of excited states whose numbers obey the approximate selection rules given in Equation (54).Note that the latter set of excited states depends on the specific N -electron operator under consideration.(For hole concentrations x > 0 this definition refers to excitation energy ω < 2µ.) The One-and Two-Electron Subspace The concept of a N -electron subspace as defined above refers to a specific operator.In contrast, rather than referring to a specific N -electron operator, the one-and two subspace is the set of N = 1 and N = 2 subspaces associated with the one-electron operator and all simple two-electron operators, respectively.Besides the N = 1 one-electron operator Ô1 ( k) = c k,σ (measured in the angle-resolved photoelectron spectroscopy), this includes a set of N = 2 operators ÔN ( k) such as the spin-projection σ density operator Ôσsd , and the charge density operator (measured in density-density electron energy loss spectroscopy and inelastic X-ray scattering).The latter operator is written in terms of the above spin-up and spin-down density operators.Moreover, the set of N = 2 operators includes several superconductivity operators whose pairing symmetries are in general different at 1D and for the square lattice.The local operators ÔN,j corresponding to the operators ÔN ( k) whose explicit expression is provided above read Ô1,j = c r j ,σ , Ôσsd 2,j = c † r j ,σ c r j ,σ , and Ôsdw 2,j = c † r j ,↑ c r j ,↓ , respectively.In the case of excitations of x > 0 and m = 0 ground states, in the remaining of this paper we are mostly interested in the subspace obtained from the overlap of the one-and two-electron subspace with the vanishing rotated-electron double occupancy subspace.Such a subspace is the one-and two-electron subspace for excitation energy ω < 2µ.For finite hole concentrations this is the subspace of interest for the one-and two-electron physics.In turn, concerning the excitations of a x = 0 and m = 0 ground state we consider the whole one-and two-electron subspace, which refers both to the spin lower-Hubbard band physics and one-electron and charge upper-Hubbard band physics.
As discussed above concerning the general N -electron spectral weight, for the model on the 1D and square lattices there is for x > 0 an extremely small amount of one-and two-electron spectral weight that for excitation energy ω < 2µ is generated by states that do not obey the approximate selection rules of Equation (54) for N = 1, 2. Nearly all such very small amount of spectral weight refers to N s3 = 1 excited states.(States with N s4 = 1 or N s2 = 2 generate nearly no spectral weight.)That very small weight is neglected within the present definition of the one-and two-electron subspace, which refers to excitation energies below 2µ.In turn, concerning the one-and two-electron excitations of a x = 0 and Reference [10], it is convenient to express the one-and two-electron excitation spectrum relative to initial x ≥ 0 and m = 0 ground states in terms of the deviations in the numbers of c effective lattice unoccupied sites and s1 effective lattice unoccupied sites.Those are given explicitly in Equation (60) and equal the corresponding deviations in the numbers of c band fermion holes and s1 band fermion holes, respectively.Note that for x > 0 and ω < 2µ states the s1 fermion related numbers provided in Equation ( 57 As discussed above, for N s2 = 1 spin-singlet excited energy eigenstates the single s2 fermion has vanishing energy and momentum.Consistent with Equation (42), for vanishing magnetic field H = 0 it is invariant under the electron-rotated-electron unitary transformation.The same applies to the single η1 fermion of N η1 = 1 η-spin-singlet excited states of the x = 0, µ = 0, and m = 0 ground state.Therefore, the only effect of creation and annihilation of such two objects is in the numbers of sites and occupied sites of the s1 effective lattice.Their creation can then be merely accounted for by small changes in the occupancies of the discrete momentum values of the s1 band.Hence the only composite object whose internal occupancy configurations are important for the physics of the Hamiltonian Equation ( 1) in the one-and two-electron subspace is the spin-neutral two-spinon s1 fermion and related spin-singlet two-spinon s1 bond particle [10].
It turns out that for the Hubbard model in the one-and two-electron subspace and alike for the s2 fermion and/or the η1 fermion, the presence of deconfined spinons is felt through the numbers of occupied and unoccupied sites of the s1 effective lattice.For excited states of x ≥ 0 and m = 0 ground states the number of deconfined η-spinons equals that of the unoccupied sites of the c effective lattice.For excited states of the x = 0, µ = 0, and m = 0 ground state the presence of deconfined η-spinons is felt in addition through the above numbers of sites and occupied sites of the s1 effective lattice.Those may be rewritten as respectively, where M de η = 2S η = 0, 1, 2. Therefore, when acting onto the one-and two-electron subspace, the Hubbard model refers to a two-component quantum liquid that can be described only in terms of c fermions and s1 fermions.For excited states of x > 0 and m = 0 ground states this analysis applies to excitation energies ω < 2µ.For those of the x = 0, µ = 0, and m = 0 ground state the η-spin degrees of freedom are in addition behind the finite energy ∆ Drot given in Equation (39) and the related energy ∆ h Drot .Those are associated with rotated-electron doubly occupied sites and rotated-electron unoccupied sites, respectively, of excited states with finite occupancy M de η = 1, 2 of deconfined η-spinons or N η1 = 1 of a single η1 fermion.However, the magnitude of that energy is fixed for each branch of excitations of that ground state.Hence for it and its excited states the square-lattice quantum liquid may again be described solely in terms of c fermions and s1 fermions.
For excited states of x > 0 and m = 0 ground states belonging to the one-and two-electron subspace the spin effective lattice has a number of sites given by N D as = (1 − x) N D a .For those of the x = 0, µ = 0, and m = 0 ground state it reads . For x > 0 its value is smaller than that of the original lattice.Within the N D a 1 limit one may neglect corrections of the order 1/N D a so that for both types of excited states the lattice spacing a s is that provided in Equation ( 30) for α = s.For the model on the square lattice it reads a s ≈ a/ √ 1 − x.Both it and its general expression given in Equation ( 30) are such that the area L 2 = [a s × N as ] 2 = [a × N a ] 2 of the system is preserved.The concept of a spin effective lattice is valid only within the N D a 1 limit that the description reviewed in this paper refers to.For the model on the square lattice, the corresponding s1 fermion spinon occupancy configurations are expected to be a good approximation provided that the ratio N 2 as /N 2 a and thus the electronic density n = (1 − x) remain finite as N 2 a → ∞.This is met for the hole concentration range x ∈ (0, x * ) considered in the studies of Reference [10].
Within the present N 2 a 1 limit there is for the one-and two-electron subspace of the model on the square lattice commensurability between the real-space distributions of the N 2 a s1 ≈ N s1 sites of the s1 effective lattice and the N 2 as ≈ 2N s1 sites of the spin effective lattice.For (1 − x) ≥ 1/N 2 a and N 2 a 1 the spin effective lattice has N 2 as ≈ (1 − x) N 2 a sites and from the use of the expression given in Equation (57) for the number of s1 effective lattice sites N 2 a s1 and Equation (46) for the corresponding spacing a s1 we find, Our general description refers to a very large number of sites N 2 a 1.Although very large, we assume that N D a is finite and only in the end of any calculation take the N D a → ∞ limit.For N D a 1 very large but finite the m = 0 ground state spin effective lattice is full and both at x = 0 and for x > 0 such a state is a spin-singlet state.(For m = 0 and x = 0 this agrees with the exact theorem of Reference [4].)For N h s1 = 0 states such as x ≥ 0 and m = 0 ground states and their charge excited states the spin effective square lattice has two well-defined sub-lattices, which we call sub-lattice 1 and 2, respectively.(For the N h s1 = 1, 2 states of the present subspace the spin effective lattice has two bipartite lattices as well, with the one or two extra sites accounted for by suitable boundary conditions.)The two spin effective sub-lattices have spacing a s1 ≈ √ 2 a s .The fundamental translation vectors of the sub-lattices 1 and 2 read, respectively.Here e x 1 and e x 2 are unit vectors pointing in the direction associated with the Cartesian coordinates x 1 and x 2 , respectively.The vectors given in this equation are the fundamental translation vectors of the s1 effective square lattice.
In the case of x ≥ 0, m = 0, and N h s1 = 0 ground states whose s1 momentum band is full and all N 2 a s1 = N a s1 × N a s1 sites of the s1 effective square lattice are occupied, we consider that the square root N as of the number N 2 as = N as × N as of sites of the corresponding spin effective square lattice is an integer.Although the square root N a s1 of the number N 2 a s1 = N a s1 × N a s1 of sites of the s1 effective lattice is not in general an integer number, within the present N 2 a 1 limit we consider that it is the closest integer to it.

The Square-Lattice Quantum Liquid of c and s1 Fermions
It follows from the above results that when acting onto the one-and two-electron subspace, the Hubbard model on a 1D or square lattice refers to a two-component quantum liquid described in terms of two types of objects on the corresponding effective lattices and momentum bands: The charge c fermions and spin-neutral two-spinon s1 fermions.The one-and two-electron subspace can be divided into smaller subspaces that conserve S c and S s .When expressed in terms of c and s1 fermion operators, the Hubbard model on a square lattice in the one-and two-electron subspace is the square-lattice quantum liquid further studied in Reference [10].
The quantum-liquid c fermions are η-spinless and spinless objects without internal degrees of freedom and structure whose effective lattice is identical to the original lattice.For the complete set of U/4t > 0 energy eigenstates that span the full Hilbert space, the occupied sites (and unoccupied sites) of the c effective lattice correspond to those singly occupied (and doubly occupied plus unoccupied) by the rotated electrons.The corresponding c band has the same shape and momentum area as the first Brillouin zone [10].
In contrast, the quantum-liquid composite spin-neutral two-spinon s1 fermions have internal structure.Thus the spinon occupancy configurations that describe such objects are a more complex problem for the one-and two-electron subspace.It is simplified by the property that the number of unoccupied sites of the s1 effective lattice is in that subspace limited to the values N h s1 = 0, 1, 2. That the square-lattice quantum liquid is constructed to inherently the c and s1 fermion discrete momentum values being good quantum numbers is behind the suitability of the present description in terms of occupancy configurations of the c and s1 effective lattices and corresponding c and s1 band discrete momentum values.The latter c and s1 values are the conjugate of the real-space coordinates of the c and s1 effective lattice, respectively.Are the approximations used in the construction of the s1 effective lattice inconsistent with the s1 band discrete momentum values being good quantum numbers?The answer is no.Indeed, such approximations concern the relative positions of the j = 1, . .., N 2 a s1 sites of the s1 effective lattice [10].Those control the shape of the s1 momentum band boundary.They do not affect the s1 band discrete momentum values being good quantum numbers.At x = 0 the spin effective lattice is identical to the original square lattice and the s1 effective lattice is one of its two sub-lattices.Consistently, at x = 0 and m = 0 the boundary of the s1 momentum band is accurately known.Then the s1 band coincides with an antiferromagnetic reduced Brillouin zone of momentum area 2π 2 such that |q In turn, it is known that for x > 0 and m = 0 the s1 band boundary encloses a smaller momentum area (1 − x)2π 2 yet its precise shape remains an open issue.
The general problem of expressing the 1D Hubbard model in terms of c and αν fermion operators is beyond the goals of this paper.Fortunately, the problem of physical interest is simpler.It refers to the expression of the Hamiltonian Ĥ of Equation ( 6) in the one-and two-electron subspace alone.This leads to an expression in terms of only c fermion and s1 fermion operators.However, this remains a problem of some complexity.It is discussed in Section 4.1 of Reference [10].The corresponding general energy spectrum involves c and s1 momentum bands and corresponding energy dispersions.Such dispersions and associated velocities are studied in that reference.

A Preliminary Application: The Inelastic Neutron Scattering of LCO
It is desirable that the results of the square-lattice quantum liquid studied in this paper are compared with those of the standard formalism of many-body physics.Unfortunately, such a quantum liquid is non-perturbative in terms of electron operators so that, in contrast to a 3D isotropic Fermi liquid [27], rewriting the theory in terms of it is an extremely complex problem.
In spite of the lack of an exact solution for the model on a square lattice and the non-perturbative character of the quantum problem in terms of electrons, in this subsection results achieved by the square-lattice quantum liquid description are compared with those obtained by methods relying on the standard formalism of many-body physics.Unfortunately, there are not many controlled results for the Hubbard model on the square lattice from approximations relying on that formalism.Here we consider the interesting problem of the spin-excitation spectrum of the half-filling Hubbard model on the square lattice.
Within the present description and as discussed below, since the c and s1 fermion momentum values are for the model in the one-and two-electron subspace good quantum numbers, that problem refers to an effectively non-interacting limit whereas in terms of electrons it is an involved many-body problem.Fortunately, there are reliable results on that particular problem obtained by controlled approximations of the standard formalism of many-electron physics: within such approximations its solution requires summing up an infinite set of ladder diagrams, to find the spin-wave dispersion of the half-filled Hubbard model on the square lattice in a spin-density-wave-broken symmetry ground state [14].In turn, within the present description the spin spectrum involves the creation of two holes in the s1 band.Agreement between the two methods is both a further checking of the validity of the description and a confirmation that the c and s1 fermion interactions are indeed residual and their momentum values good quantum numbers for the model on the square lattice in the one-and two-electron subspace.
The spin-triplet excitations relative to the x = 0 and m = 0 absolute ground state involve creation of two holes in the s1 band along with a shift π/N 2 a of all discrete momentum values of the full c band so that the their spectrum reads, where π = ±[π, ±π] and for approximately U/4t ≥ 1.3, excitation energy ω below µ 0 , and temperatures T below ω/k B , the s1 fermion energy dispersion s1 ( q) reads [10], Here the maximum gap magnitude µ 0 /2 of the s1 fermion energy dispersion equals the maximum s1 fermion pairing energy per spinon for x = 0 and m = 0.It corresponds to the s1 band momentum values belonging to the s1 boundary line and pointing in the anti-nodal directions.The x = 0 and m = 0 s1 energy dispersion Equation ( 64) is derived in Reference [10].The main point is that in the absence of any spin order, the s1 fermion spinon-pairing energy would vanish and the energy dispersion s1 ( q) of the s1 fermions would be given by the auxiliary energy dispersion 0 s1 ( q), also provided in that equation.An important issue for the Reference [10] derivation of the s1 fermion energy dispersion s1 ( q) in the presence of such an order as N 2 a → ∞, is that for excitation energy below the energy scale µ 0 one can ignore the amplitude fluctuations of the order parameter and the problem can be handled by a suitable mean-field theory where the occurrence of that order is described for the x = 0 and m = 0 problem by a s1 energy dispersion of the general form − | 0 s1 ( q)| 2 + |∆ s1 ( q)| 2 , where |∆ s1 ( q)| is the s1 fermion pairing energy per spinon given in Equation (64).
We emphasize that for small U/4t values the energy scale µ 0 /2 becomes small and given by µ 0 /2 ≈ 16 t e −π √ 4t/U , so that the amplitude fluctuations of the order parameter cannot be ignored and thus the energy dispersion provided in Equation ( 64) is not expected to be a good approximation.Indeed and as discussed in Reference [10], within the c and s1 fermion description the small-U/4t physics corresponds to a non-trivial problem.This applies to the derivation of the s1 fermion dispersion s1 ( q) for small values of U/4t, which remains an unsolved issue.Fortunately, the energy dispersion Equation ( 64) is expected to be a good approximation for U/4t ≥ 1.3.Furthermore, as confirmed in the following the relation of the present quantum problem to the unusual physics of the Mott-Hubbard insulator parent compounds such as LCO refers to a value U/4t ≈ 1.525 > 1.3, so that the description introduced in References [9,10] is of interest for the study of such materials.For x = 0 and m = 0 both the c and s1 bands are full for the initial ground state and since the c band remains full for the excited states one can ignore the s1-s1 and s1-c fermion interactions.Indeed, then the residual fermion interactions have little effect on the occupancy configurations of the two holes created in the s1 band upon the two-electron spin-triplet excitations.This is consistent with the lack of a c Fermi line for the initial ground state and the lack of s1 band holes other than the two holes created upon the spin-triplet excitation so that in spite of the s1-s1 fermion long-range interactions associated with the effective vector potential of Equation (33) the exclusion principle, phase-space restrictions, and momentum and energy conservation drastically limit the number of available momentum occupancy configurations of the final excited states.
The excitation spectrum Equation (63) refers both to coherent and incoherent spin spectral weight.In contrast to the 1D case where a suitable c and s1 fermion dynamical theory is available [28,29], for the square-lattice quantum liquid there are within the present status of the theory no suitable tools to calculate matrix elements between the ground state and one-and two-electron excited states.Hence, one cannot calculate explicitly spin-spin correlation functions.Within the 1D c and s1 fermion dynamical theory, the sharp features of the spin two-electron spectral weight distributions result from processes where one of the two s1 fermion holes is created at the s1 boundary line.
The coherent spin spectral weight is here associated with a Goldstone-mode-like gapless spin-wave spectrum.It consists of sharp δ-peaks having as background the incoherent spectral-weight distribution.From comparison with the results of Reference [14], we have confirmed that such spectral weight is generated by processes corresponding to well-defined values of the momenta q and q of the general spectrum Equation (63) such that one hole is created at a momentum pointing in the nodal directions of the s1 band and the other hole at a momentum belonging the s1 band boundary line, as expected from analogy with the 1D spectral-weight distributions.The incoherent part corresponds to the remaining values of q and q of the excitation spectrum Equation (63).The occurrence of the Goldstone-mode-like gapless spin-wave spectrum follows from the long-range antiferromagnetic order of the initial x = 0 and Figure 5.The theoretical spin spectrum Equation (65) (solid lines) plotted for the high symmetry directions in the second Brillouin zone for µ 0 = 565.6 meV and W 0 s1 = 49.6 meV and the experimental data of Reference [13] (circles) in meV.Such theoretical magnitudes correspond to t ≈ 0.295 eV and U ≈ 1.800 eV, so that U/4t ≈ 1.525.The momentum is given in units of 2π.The corresponding theoretical lines plotted in Figure 5 of Reference [14] are very similar to those plotted here yet are obtained within the standard formalism of many-body physics by summing up an infinite number of ladder diagrams.The curves plotted in Figure 5 of Reference [14] refer to the high symmetry directions in the Brillouin zone.The use of the above general expression given in Equation (65) leads for U/4t = 1.525 to an excellent agreement with such curves.These directions correspond also to those measured by highresolution inelastic neutron scattering in LCO, as plotted in Figure 3(A) of Reference [13].We denote such symmetry directions by M O, ΓO, XM , ΓX, and XO.They connect the momentum-space points M = [π, π], O = [π/2, π/2], Γ = [0, 0], and X = [π, 0] of the general spin-wave spectrum provided in Equation (65).
Such a spin-wave excitation spectrum is plotted for the above symmetry directions in Figure 1 (solid line) for µ 0 = 565.6 meV and W 0 s1 = 49.6 meV together with the experimental results (circles) for T = 10 K.This gives a Mott-Hubbard gap 2µ 0 = 1131.2meV.The spin-spectrum expression given in Equation (65) refers to the first Brillouin zone.In Figure 1 we plot it in the second Brillouin zone, alike in Figure 3(A) of Reference [13] and Figure 5 of Reference [14].An excellent quantitative agreement is reached for these magnitudes of the involved energy scales, which according to the complementary results of Reference [14] correspond to U/4t ≈ 1.525 and t ≈ 295 meV.
The above results confirm that the predictions of the square-lattice quantum liquid theory concerning the spin spectrum at half filling agree both with experiments on the parent compound LCO and results obtained by the standard formalism of many-body physics.

Concluding Remarks
In this paper a general operator description valid for both the Hubbard model on the 1D and square lattice in terms of three types of elementary quantum objects introduced in References [9,10] was reviewed.Such a rotated-electron related operational description is an application of the Hubbard model on a bipartite lattice extended global SO(3) × SO(3) × U (1) symmetry recently found in Reference [6].The occupancy configurations of such objects correspond to state representations of such an extended global symmetry.The physical interest of the present rotated-electron related operational description refers mostly to the model in the one-and two-electron subspace.However, the definition of that subspace requires the use of general properties of the full-Hilbert-space description summarized in this paper.
Although the operator description reviewed in this paper is compatible with and in part inspired in the exact solution of the 1D model, it accounts for the basic differences between the physics of the Hubbard model on the 1D and square lattice, respectively.For instance, in 1D the occurrence of an infinite set of conservation laws associated with the model integrability implies that the c-s1 fermion residual interactions refer only to zero-momentum forward-scattering.They merely give rise to phase shifts whose expressions may be extracted from the BA solution.This allows the introduction of a pseudofermion dynamical theory, which provides finite-energy spectral and correlation function expressions involving phase shifts [28,29].Hence in 1D such interactions do not involve interchange of energy and momentum.In contrast, they do for the Hubbard model on the square lattice, yet they are much simpler than the corresponding electronic correlations.Indeed the quantum problem is non-perturbative in terms of electron operators.It follows that in contrast to a 3D isotropic Fermi liquid [27], rewriting the square-lattice quantum-liquid theory in terms of the standard formalism of many-electron physics is in general an extremely complex problem.Fortunately, such a quantum liquid dramatically simplifies when expressed in terms of the c fermion and s1 fermion operators [10].
The problem is simplest at x = 0 for spin excitations for which the c band remains full and the effects of the c-s1 fermion interactions are frozen.The preliminary investigations reported in this paper on the physical consequences of the model on the square lattice new found global symmetry in actual materials in terms of the c and s1 fermion description refer to x = 0.Such results confirm that the present description is useful for the further understanding of the role played by the electronic correlations in the spin spectrum of the parent compound LCO [13].Indeed, it is quantitatively described by the corresponding spin spectrum of the square-lattice quantum liquid at U/4t ≈ 1.525 and t ≈ 295 meV.
A similar good quantitative agreement with the LCO spin-wave spectrum is reached by use of the Hubbard model on the square lattice for U/4t values in the range U/4t ∈ (1.5, 2.0) and suitable different magnitudes of t.For instance, the studies of References [10,14] find U/4t = 1.525 alike here whereas those of References [30] and [31] use U/4t = 1.625 and U/4t = 2.000, respectively.Consistent with the results of References [17,18], the Hubbard model in the subspace of vanishing rotated-electron double occupancy is related to the t-J model only for larger U/4t values.Importantly, for arbitrary U/4t > 0 values the usual large-U energy scale 4t 2 /U is replaced by the s1 fermion energy nodal bandwidth W s1 of Equation (36).Its maximum magnitude W s1 = 4t is reached for U/4t → 0. For U/4t > 0 it decreases monotonously for increasing values of U/4t, vanishing for U/4t → ∞ as W s1 ∝ 4t 2 /U [10].Hence the 4t 2 /U magnitude found here and in References [10,14,30,31] for the description of the LCO spin-wave spectrum within the Hubbard model should not be compared with that used by t-J model descriptions of the same spectrum.Indeed for the above intermediate U/4t values there is no direct relation between the two models.
Further investigations on more complex x > 0 and m = 0 2D problems for which the effects of doping are accounted for in terms of c-s1 fermion residual inelastic interactions will be also fulfilled elsewhere.
x 3 η = −[N D a − N ]/2 and S x 3 s = −[N ↑ − N ↓ ]/2 are the η-spin projection and spin projection, respectively, and N D a 1 denotes the number of lattice sites.For the bipartite 1D and square lattices considered in this paper the labeling index D in N D a ≡ [N a ] D reads D = 1 and D = 2, accounting for the N a 1 and N 2 a = N a × N a 1 lattice sites, respectively.The square and 1D lattices have spacing a and length edge and chain length L = N a a, respectively.
spin-1/2 η-spinons, and N c = 2S c spin-less and η-spin-less charge c fermions.The latter live on a lattice with N D a = [N c + N h c ] sites identical to the original lattice.Here N h c = [N D a − 2S c ] gives the number of c fermion holes.The relation of such objects to the rotated electrons is as follows.The M s = 2S c spin-1/2 spinons describe the spin degrees of freedom of the 2S c rotated electrons that singly occupy sites.The charge degrees of freedom of such rotated electrons are described by the N c = 2S c c fermions.The M η = [N D a − 2S c ] η-spin-1/2 η-spinons describe the η-spin degrees of freedom of the [N D

2 .
The Model, a Suitable Rotated-Electron Description, and Relation to the Global SO(3) × SO(3) × U (1) Symmetry sites of the η-spin effective lattice and N D as = 2S c sites of the spin effective lattice.(ii) The remaining degrees of freedom of rotated-electron occupancies of the sets of N D aη = [N D a − 2S c ] and N D as = 2S c original-lattice sites correspond to the occupancy configurations associated with the η-spin SU (2) symmetry and spin SU (2) symmetry state representations, respectively.The occupancy configurations of the set of N D aη = [N D a − 2S c ] sites of the η-spin effective lattice and set of N D as = 2S c sites of the spin effective lattice are independent.
) with N D aη = N D a and N D as = 2S c = 0.For U/4t > 0 only for a m = (1 − x) fully polarized state are the occupancy configurations of the state |GS c ; 2S c and the corresponding N c = 2S c c fermions invariant under the electron-rotated-electron unitary transformation.For the corresponding vacuum |0 η ; N D aη (and |0 s ; N D as ), the M η = M de η,+1/2

3. 4 .
Spacing and Occupied and Unoccupied Sites of the η-Spin and Spin Effective Lattices One may consider a configuration in which a original-lattice compact square domain of N D as = 2S c = [2S s + M co s ] = N c sites whose edge contains N as sites is singly occupied by rotated electrons, and the complementary two

4. 1 .
The M co η -η-Spinon and M co s -Spinon Configuration PartitionsWithin the rotated-electron occupancy configurations that generate the exact energy eigenstates |Ψ U/4t = V † |Ψ ∞ considered above, there are [M co η +M co s ] sites out of the N D a sites of the original lattice whose rotated-electron occupancy configurations are not invariant under the electron-rotated-electron unitary transformation.The η-spin (and spin) degrees of freedom of M co η (and M co s ) of such sites refer to η-spin-neutral (and spin-neutral) configurations involving M co η /2 +1/2 η-spinons (and M co s /2 + 1/2 spinons) and an equal number M co η /2 of −1/2 η-spinons (and M co s /2 of −1/2 spinons).

Figure 4 .Transformation
Figure 4.The c fermion band for the 1D Hubbard model in units of t plotted for electronic densities n = (1 − x) = 1/2, 5/6, spin density m = 0, and a set of U/t values.The ground-state energy level was chosen to correspond to zero energy.It is marked by a horizontal line which overlaps the c band at the c Fermi points q = ±q F c = ±2k F = π n = π (1 − x).[26] ) and 9(a) (and Figures 8(b) and 9(b)) refers to the energy 2|µ| (and 4|µ|) of the invariance condition given in Equation (42) for ν = 1 (and ν = 2).
Figure 5.The theoretical spin spectrum Equation (65) (solid lines) plotted for the high symmetry directions in the second Brillouin zone for µ 0 = 565.6 meV and W 0 s1 = 49.6 meV and the experimental data of Reference[13] (circles) in meV.Such theoretical magnitudes correspond to t ≈ 0.295 eV and U ≈ 1.800 eV, so that U/4t ≈ 1.525.The momentum is given in units of 2π.The corresponding theoretical lines plotted in Figure5of Reference[14] are very similar to those plotted here yet are obtained within the standard formalism of many-body physics by summing up an infinite number of ladder diagrams.[10] Thus their energies are additive.It follows that for all m values 2νµ B |H| is the energy of ν deconfined −1/2 spinons and ν deconfined +1/2 spinons.For example, for m > 0 the energy 2νµ B H is that for creation of a number ν = 1, 2, . . . of deconfined −1/2 spinons onto a S s = ν ground state with 2ν deconfined +1/2 spinons.Such a creation refers to ν spin-flip processes.The number i sν decreases continuously for increasing values of U/4t.For any fixed m value it has the limiting behaviors i sν → 1 for U/4t → 0 and i sν