arxiv:cond-mat/ v1 30 Sep 1996

Transcription

1 Complete Pseudoparticle Operator Representation for Integrable Electronic Luttinger Liquids J. M. P. Carmelo 1, N. M. R. Peres 1,, and D. K. Campbell arxiv:cond-mat/960991v1 30 Sep Department of Physics, University of Évora, Apartado 94, P-7001 Évora Codex, Portugal Department of Physics, UIUC, 1110 West Green Street, Urbana, IL (Received August 1996) We introduce the pseudoparticle operator algebra that generates all Hubbardchain eigenstates from a single reference vacuum. The occurence of only zeromomentum pseudoparticle forward scattering at all energy scales combined with the infinite pseudoparticle-number conservation laws establish this operator basis as the most suitable for dealing with transport properties at finite frequencies. We discuss the explicit form of the finite-frequency conductivity in these models. PACS numbers: Nj, Fk, y, Ca Typeset Using REVTEX 1

2 The unusual spectral properties of Luttinger liquids [1] imply that they can not be described by one-electron, Fermi-liquid-like models. Further, descriptions of these quantum systems in terms of exotic excitations such as holons and spinons, or pseudoparticles, have been limited to the low-energy Hilbert subspace, and both bosonization and conformal-field theory techniques also apply only in that limit. Thus developing a microscopic operator description of these systems for arbitrary energies remains an important theoretical challenge. In this Letter we show that, for the integrable Luttinger liquid associated with the Hubbard chain [,3] in a magnetic field, a pseudoparticle operator algebra provides the desired microscopic description for all states and all energies. This algebra is a generalization of the low-energy pseudoparticle operator representation of Refs. [4 6]. It goes beyond the BA solution [,3] because it i)- generates from a single vacuum all Hamiltonian eigenstates, including the η-spin and spin non-lowest-weight states (non-lws s) and non-highest-weight states (non-hws s) of the SO(4) algebra [7] absent in that solution; ii)- provides relevant information on non-diagonal operators; and iii)- permits the study of the finite-energy transport properties. Using this algebra, we obtain explicit expressions for the Drude peak (for the non-half-filled cases) and for a prominent feature that sets in with a fractional power law at ω µ, where µ is the chemical potential. Our results confirm the recent indications that the integrability of BA solvable many-electron problems is closely related to peculiar spectral and transport properties [8]: we find that at all energy scales there is only zero-momentum pseudoparticle forward scattering in the Hubbard chain. Although this exact property holds only for integrable systems, we believe that our pseudoparticle operator description will provide the most suitable starting point for description of the unusual finite-frequency transport properties of both integrable and non-integrable electronic Luttinger liquids. We consider the Hubbard chain [,3] in a magnetic field H and chemical potential µ [ 6], Ĥ = ĤSO(4) + µˆη z + µ 0 HŜz, where ĤSO(4) = t [ j,σ c jσc j+1σ + c ] j+1σc jσ + U j[c j c j 1 ][c j c j 1]. The generators of the η-spin and spin SO(4) algebra are [7], ˆη z = 1[N a

3 ˆN σ σ ], ˆη = j( 1) j c j c j, ˆη + = j( 1) j c j c j for η-spin and Ŝz = 1 σ σ ˆN σ, Ŝ = j c j c j, Ŝ + = j c j c j for spin. Here c jσ and c jσ are electron operators of spin projection σ at site j, ˆNσ = j c jσc jσ and, t, U, µ, H, and µ 0 are the number of sites, the transfer integral, the onsite Coulomb interaction, the chemical potential, the magnetic field, and the Bohr magneton, respectively. We consider the electronic density n = n + n and the spin density m = n n, with n σ = Nσ. The related holon and spinon study of Ref. [9] only refers and applies to states with eigenvalues η z 0 and S z 0. Below we generate all Hamiltonian eigenstates from the pseudohole vacuum V introduced in Ref. [6], which is the half-filled and zero-magnetic-field ground state. Our general description uses the α, β pseudoholes introduced in Ref. [6], which have creation operators a q,α,β. In addition, the general description involves an infinite number of α, γ pseudoparticle branches, with creation operators b q,α,γ. Here α = c, s and β = ±1 [6] and the quantum numbers γ are the positive integers (γ = 1,,..., ). The pseudoparticle/hole quantum numbers α, β, and γ can in some limits be related to electronic numbers: in the n 1 and η 0 (or m 0 and S 0) limit, α = c (or α = s) becomes charge (or spin) [5], and the c, β (or s, β) pseudoholes are such that β = η z = ± 1 (or β = S z = 1 ) [6]. Moreover, for U, the creation of one c, γ pseudoparticle (γ > 0) is associated with the creation of D = γ doubly occupied sites ([Ĥ, ˆD] = 0 in that limit). The above operators obey an anticommuting algebra such that {a q,α,β, a q,α,β } = δ q,q δ α,α δ β,β and {b q,α,γ, b q,α,γ } = δ q,q δ α,α δ γ,γ and all remaining anticommutators vanish. In contrast to their holes, the usual α pseudoparticles (of number N α [4,5]) are not labeled by the β pseudohole quantum number [6]. The α, + 1 and α, 1 pseudoholes occupy different pseudomomentum values in the same α band [as the right-hand side of Eq. () below confirms]. For convenience we shall call it the α, 0 band. We thus have a set of α, γ pseudomomentum bands which are occupied by the α, β pseudoholes (and α pseudoparticles) for γ = 0 and by the α, γ pseudoparticles for γ > 0. [Below when we refer α, γ pseudoparticles (bands) we consider γ = 1,,... (γ = 0, 1,,...) whereas in the case of the N α,γ numbers, we can consider either 3

4 case, with N α,0 N α ]. A crucial point is that the pseudohole and pseudoparticle numbers are constant for each Hamiltonian eigenstate. The N α,γ numbers are related to the BA numbers of Ref. [3] as N c,0 N M, N c,γ M γ (with γ > 0), and N s,γ M 1+γ (with γ = 0, 1,...). However, to study non-lws s outside the BA solution [7] requires the introduction of the new α, β pseudohole numbers. These are given by N h c,β = Nh c β[ N] and N h s,β = Nh s β[n N ] where N h c = N c,0 and N h s = N c,0 γ=0 N s,γ. These important expressions are generalizations of the pseudohole expressions of Ref. [6], which refer only to the η and S LWS s and HWS s I (which we call simply states I). The operators ˆN h α,β (q) = a q,α,β a q,α,β and ˆN α,γ (q) = b q,α,γ b q,α,γ (with γ = 1,,...), and ˆN α,0 (q) 1 ˆN β α,β h (q) play a central role in the present theory. The number operators read ˆN h α,β = q ˆN h α,β(q) and ˆN α,γ = q ˆN α,γ (q). We introduce the notations η = S c and S = S s for η-spin and spin, respectively. The SO(4) operator generators, Ŝα z and Ŝα ±, and the η-spin and spin numbers, S α, are given by Ŝz α = β ˆN α,β h (q) ; q,β Ŝα ± = q a q,α,± a 1 q,α, 1 ; S α = 1 [Nh α γn α,γ ]. (1) γ=1 We emphasize that the expressions for the off-diagonal SO(4) generators go beyond BA, which provides the expressions only for diagonal operators. The momentum simply reads ˆP = γ=0 q,α qc α,γ ˆNα,γ (q) + π[ ˆN h + ˆN c, 1 γ=1 c,γ ], where C c,γ = 1 for γ > 0 and C α,γ = 1 otherwise. The usual first-quantized BA wave function expression for the Hamiltonian eigenstates has an involved form which includes many permutations [,3]. Although the pseudoparticles/holes obey generalized [10] rather than Fermi or Bose statistics, one can use their formal anticommutation relations (together with suitable topological-momenton operators [6]) to develop an extremely simple second-quantized representation of all 4 Na orthonormal Hamiltonian eigenstates: ψ; {N h α,β}, {N α,γ } = 1 C α [Ŝα +] Nh α, 1 γ=1 a b q,α, 1 q,α,γ q,q V, () 4

5 where the generator Ŝα + is given in Eq. (1), {Nh α,β } and {N α,γ} are abbreviations for the sets {N h, N h, N h, N h } and {N c, 1 c, 1 s, 1 s, 1 c,1,..., N c,, N s,1,..., N s, }, respectively, and C = α[nα!/n h h!]. The number of discrete pseudomomentum values in each α, γ band (α, γ- α, 1 band Fock-space dimension) is given by d Fα,γ = N h α +N α,γ γ =1 [ min{γ, γ }]N α,γ. In the above summations and products, the pseudomomentum takes on values q j = π in contrast to the usual momentum, I α,γ j N α,γ odd and even, respectively. Here N c,0 = α I α,γ j, where are consecutive integers or half-odd integers for N h α γ=1 N c,γ and N α,γ = d Fα,γ otherwise. It follows that for each α, γ band q α,γ ( ) q q(+) α,γ, with the limits of the pseudo-brillouin zones given by q (±) c,0 = ±π[1 1 ] for N c,0 even and q (+) c,0 = π (or q (+) c,0 = π[1 ]) and q ( ) c,0 = π[1 ] (or q ( ) c,0 = π) for N c,0 odd and simply by q (±) α,γ = ± π [d Fα,γ 1] for all the remaining α, γ bands. A S α LWS (or HWS) is such that N h α, 1 = γ=1 γn α,γ (or N h α, 1 = γ=1 γn α,γ ) with N h α, 1 N h α, 1 (or N h α, 1 N h α, 1 ). For the states I, we have that N α,γ = 0 for all γ = 1,,... branches in the former equation and for S α LWS (or HWS s) II [5,6] (which we call states II), N α,γ > 0 for at least one γ = 1,,... branch. On the other hand, the S α non-lws s and non-hws s outside the BA are such that N h α γ=1 γn α,γ > N h α, 1 N h. Moreover, α, 1 the conservation of the pseudohole and pseudoparticle numbers permits dividing the Hilbert space into subspaces spanned by the set of Hamiltonian eigenstates () with the same {Nα,β h } and {N α,γ } numbers, which correspond to the same sub-canonical ensemble. Obviously, a canonical ensemble (with constant electron numbers) is usually realized by several {Nα,β}, h {N α,γ } sub-canonical ensembles. Let us introduce the generalized ground state (GGS) as the Hamiltonian eigenstate(s) () of lowest energy in each Hilbert subpace. We find it is of the form GGS; {N h α,β}, {N α,γ } = 1 C α [Ŝα +] Nh α, 1 γ=1 q ( ) F α,0 q=q ( ) α,0 q (+) α,0 q= q (+) F α,0 q (+) F α,γ q =q ( ) F α,γ a q,α, b 1 q,α,γ V, (3) where q (±) Fα,γ = q (±) Fα,γ ± π and when N α,γ is odd (even) and I α,γ j are integers (half-odd integers) the pseudo-fermi points are symmetric and given by q (±) Fα,γ = ±[q Fα,γ π ] where 5

6 q Fα,γ = πnα,γ. On the other hand, when N α,γ is odd (even) and I α,γ j are half-odd integers (integers) we have that q (+) Fα,γ = q Fα,γ and q ( ) Fα,γ = [q Fα,γ π ], or q (+) Fα,γ = q Fα,γ π and q ( ) Fα,γ = q Fα,γ. The ground state (GS) associated with a given canonical ensemble is always a state I of the form presented in Ref. [6], which is a particular case of the general GGS expression (3). To classify the statistics of our pseudoholes and pseudoparticles according to the generalized Pauli principle of Ref. [10], we introduce the pseudohole and pseudoparticle dimensions d h α 1+d F α,0 N h α and d α,γ 1+d Fα,γ N α,γ, respectively. We find d h α = α g α,α Nh α α γ =1 g α;α,γ N α,γ and d α,γ = α g α,γ;α Nh α α γ =1 g α,γ;α,γ N α,γ where the the statistical-interaction matrix has infinite dimension and reads g c,α = δ c,α, g c;α,γ = 0, g s,α = 1, g s;α,γ = δ s,α, g α,γ;α = δ α,α, g α,γ;α,γ = δ α,α min{γ, γ }. This fully defines the statistics, the pseudoholes and pseudoparticles being neither fermions nor bosons [10]. Let us introduce the α and α, γ topological-momenton operators (for the states I see Refs. [6,11]) which generate topological excitations associated with a shift of all the pseudomomenta of the corresponding band by ± π. Creation or annihilation of odd numbers of pseudoholes or pseudoparticles always requires the occurrence of such topological momentons. Their unitary generators are such that U α ±1a q,α,β = a ±1 q± π,α,βu α and U α,γ ±1b q,α,γ = b ±1 q± π,α,γu α,γ Na Na and read U ±1 α = exp[ q,β a q± π Na,α,βa q,α,β] and U ±1 α,γ = exp[ q b q π Na,α,γb q,α,γ]. A remarkable property is that the transition between a GS and any eigenstate () can be separated into two types of excitations: (a) a topological GS GGS transition which involves the creation or annihilation of pseudoholes and (or) pseudoparticles as well as the occurrence of topological momentons and (b) a Landau-liquid excitation [5,6] associated with pseudoparticle - pseudohole processes relative to the GGS. The topological transitions (a) are basically superpositions of three kinds of elementary transitions: (i) GS - GS transitions involving changes in the σ electron numbers by one, (ii) single β flip processes of α, β pseudoholes which lead to non-lws s and non-hws s outside the BA, and (iii) creation of single α, γ pseudoparticles at constant values of η, η z, S, and S z. While the transitions (i) are gapless, the elementary 6

7 excitations (ii) and (iii) require a finite amount of energy. The generators of the excitations (i) were studied in Ref. [6] and the ones of (ii) involve only topological momentons and β flips (which describe either creation of electron pairs or spin flip processes). Consider as an example of a transition (iii) the creation of one c, γ pseudoparticle which we find below to be the relevant process for the finite-frequency conductivity. (We examine here explicitly topological momentons of occupied bands only, i.e. these which generate momentum.) When γ > 0 is even (or odd) the generator reads U c ±1 G c,γ (or G c,γ U s ±1 ), where the topological momenton operator was defined above and G c,γ = [Ŝc +] γ q( ) F c,0 + π Na (γ 1) q (+) F c,0 q=q ( ) q (+) F c,0 F c,0 π Na (γ 1) a b q,c, 1 q =0,c,γ. Therefore, such transitions involve one c or s topological momenton, the creation of a number γ of c, 1 pseudoholes and γ of c, 1 pseudoholes, and the creation of one c, γ pseudoparticle at q = 0. Given a sub-canonical ensemble, we can write the Hamiltonian expression in normal order relative to the GS of the associated canonical ensemble. We find : Ĥ := Ĥtop + j=1 Ĥ (j) where Ĥtop has eigenvalue ω 0 = E GGS E GS and corresponds to the GS - GGS transition (a) and the remaining terms describe the Landau-liquid excitations (b) relative to the GGS. This Hamiltonian is suitable to study the physics at low positive energy above the gap, i.e., small ω ω 0. When several GGS s have the same energy gap (ω 0 ) the physics involves the Hilbert subspace spanned by all eigenstates () associated with the corresponding different sub-canonical ensembles. For sub-canonical ensembles such that S α z 0 and characterized by small values of N h /N α, lα α h (where l α = Sz α/ Sα z = ±1) we find Ĥtop = µ ˆN h c, lc µ 0 H [ ˆN h + s, ls γ=1 ˆNs,γ ] and thus the gap becomes + ω 0 = µ N h c, lc + µ 0 H [N h s, ls + N s,γ ]. (4) γ=1 For transitions to pure LWS s or HWS s II, this leads to ω 0 = µ γ=1 γn c,γ + µ 0 H γ=1 [1 + γ]n s,γ. Therefore, the single (γ > 0) α, γ heavy-pseudoparticle gap is γ µ and [1 + γ]µ 0 H for c and s, respectively. The relevant terms of the Landauliquid Hamiltonian at low energy ω ω 0 are Ĥ(1) = q,α,γ ǫ 0 α,γ(q) : ˆNα,γ (q) : and Ĥ () = 1 q,α,γ q,α,γ 1f α,γ;α,γ (q, q ) : ˆNα,γ (q) :: ˆNα,γ (q ) :. We can also write the 7

8 expression of Ĥ (j) for higher scattering orders j >. The main feature is that for all sub-canonical ensembles, energy scales, and pseudoparticle scattering orders, the Landauliquid terms of the above Hamiltonian involve only zero-momentum forward-scattering. The pseudoparticle bands ǫ 0 α,γ (q), associated group velocities v α,γ(q), two-pseudoparticle forward-scattering phase shifts Φ α,γ;α,γ (q, q ) (whose expressions are similar to the ones of the corresponding γ = 0 functions [4] and will be studied elsewhere) and light velocities v α,γ v α,γ (q Fα,γ ) control the physics at low energies ω ω 0. For instance, the f functions read f α,γ;α,γ (q, q ) = v α,γ (q)φ α,γ;α,γ (q, q ) + v α,γ (q )Φ α,γ ;α,γ(q, q) + γ j=±1 α =0 v α,γ Φ α,γ ;α,γ(jq Fα,γ, q)φ α,γ ;α,γ (jq Fα,γ, q ). Let us apply our pseudoparticle operator basis to charge transport in this system. We can write the charge-current operator as Ĵ = Ĵd+Ĵoff, where Ĵd = l l Ĵ l l l and Ĵoff = l l l Ĵ l l l and the summations run over all Hamiltonian eigenstates. By applying a vector potential, we can derive the expression for Ĵd. As at low energy [4], the pseudoparticles are the transport carriers at all energy scales and the electronic degrees of freedom couple to the potential through them. We find that the couplings of the α, 0 α and α, γ pseudoparticles to the potential read ec ρ α,γ where e is the electronic charge, the superscript ρ indicates the coupling to charge, and C ρ c,γ = δ γ,0+γ and C ρ s,γ = 0. Like the Hamiltonian, the GS-normal-ordered expression for Ĵd has an infinite number of terms, : Ĵ d := i=1 Ĵ (i) d. In this case only the first-order term is relevant at low positive energy ω ω 0. It reads Ĵ(1) d = e q,α γ=0 j α,γ (q) : ˆNα,γ (q) : where the pseudoparticle elementary charge currents are given by j α,γ (q) = [ γ α =0 Cρ α,γ vα,γ (q)δ α,α δ γ,γ + 1 π j=±1 jf α,γ;α,γ (q, jq Fα,γ )], and involve only the group velocities and f functions. The pseudoparticle current j α,γ j α,γ (q Fα,γ ) and charge transport mass m ρ α,γ q Fα,γ/[C ρ α,γ j α,γ], such that m ρ c,γ = q Fα,γ/[(δ γ,0 + γ)j α,γ ] and m ρ s,γ =, play central roles in charge transport. (The γ = 0 expressions recover the results of [4].) The frequency-dependent conductivity, Reσ(ω), involves transitions induced by Ĵoff and for the Hubbard model (in any dimension) in principle consists of three different pieces: (1) 8

9 the Drude absorption peak (absent for n = 1), () a finite-frequency absorption associated with creation of doubly-occupied sites, and (3) an absorption continuum between the Drude peak and absorption () (absent for n = 1). However, numerical studies show that (in contrast to D) the 1D Hubbard model has only an extremely small type-(3) absorption [13]. Since the conductivity sum rule [1] proves that there is spectral weight away of the Drude peak, it must be mostly concentrated in the absorption (). This picture is confirmed by our exact study. Although we can not yet derive a closed-form expression for Ĵoff, our operator basis does provide important information about the conductivity spectrum: namely, 1) The occurrence of only forward scattering in the Landau-liquid term of the pseudoparticle Hamiltonian is consistent with the extremely small type-(3) absorption found numerically [1,13] and indicates that the topological term and the associated finite-energy transitions (a) [followed by excitations (b)] are relevant for that spectrum at finite energy; ) That the current commutes with η, S, ˆη z, and Ŝz restricts the relevant topological excitations (a) to superpositions of elementary transitions (iii); 3) That the s, γ pseudoparticles do not couple to charge further limits the transitions (iii) to the ones associated with creation of c, γ pseudoparticles; and 4) Finally, based on the fully non-dissipative character of the pseudoparticle interactions at all energy scales combined with the occurrence of the infinite pseudohole- and pseudoparticle-number conservation laws, the important contributions to the conductivity spectrum occur for small values of ω and for values of ω just above the finite frequencies of general form ω 0 = µ γ=1 γn c,γ. For small ω, the charge - charge correlation function leads to the Drude peak [4,1] only, ( e) π[q Fc,0 /m ρ c,0]δ(ω). For low positive energies ω ω 0, we find the following general expression for type-() absorption: Reσ(ω) [ω ω 0 ] ζ, (5) where the exponent ζ is a non-classical combination of pseudoparticle phase shifts and ζ > 1 for all c, γ transitions except for the single c, 1 elementary transitions, for which 1 ζ < 1. (For large U such elementary processes create one doubly-occupied site.) Since the absorption associated with other energies and excitations is extremely small, we thus 9

10 conclude that in the 1D model (in contrast to the D case where the absorption (3) is important [13]) the spectral weight is mostly concentrated at the ω = 0 Drude peak (for n 1) and at an absorption starting at ω µ, where the conductivity shows a Luttingerliquid band edge with infinite slope at ω = µ. Based on the form of the c, 1 excitation spectrum we predict this absorption to have a width W such that 0 < W 8t. W vanishes both for n 0 and n and reaches its maximum value, W = 8t, for n 1. At n = 1 this fully agrees with the numerical studies of Ref. [14]. [Considerably smaller weight exists at the remaining absorptions of higher energy and vanishing slope; the slope also vanishes for the extremely small absorption (3), which for the doped Mott-Hubbard insulator and ω 0, is predicted to have an ω 3 onset [13].] For fixed values of U and m, µ changes between µ = MH for n 1 and a maximum value of µ = MH + 8t for both n 0 and n, where MH is the half-filled Mott-Hubbard gap []. The exponent approaches its minimal value, ζ = 1, when the amount of spectral weight away of the Drude peak is largest (at zero field and finite U, ζ = 1 for half filling) and ζ 1 in the limits U 0, U (at m = 0), and m n when all spectral weight is transferred to the Drude peak [and the absorption () disappears]. We thank J.M.B. Lopes dos Santos and J. Tinka Gammel for illuminating discussions and the US NSF (Grant DMR ) for partial support. N.M.R.P. thanks the Luso- American and Gulbenkian Foundations for financial support. 10

11 REFERENCES 1 F. D. M. Haldane, Phys. Rev. Lett. 45, 1358 (1980); J. Phys. C (1981) Elliott H. Lieb and F. Y. Wu, Phys. Rev. Lett. 0, 1445 (1968). 3 M. Takahashi, Prog. Theor. Phys. 47, 69 (197). 4 J. M. P. Carmelo and P. Horsch, Phys. Rev. Lett. 68, 871 (199); J. M. P. Carmelo, P. Horsch, and A. A. Ovchinnikov, Phys. Rev. B 46, (199). 5 J. M. P. Carmelo, A. H. Castro Neto, and D. K. Campbell, Phys. Rev. Lett. 73, 96 (1994); ibid 74, 3089 (1995). 6 J. M. P. Carmelo and N. M. R. Peres, Nucl. Phys. B 458, 579 (1996). 7 Fabian H. L. Essler, Vladimir E. Korepin, and Kareljan Schoutens, Phys. Rev. Lett. 67, 3848 (1991). 8 H. Castella, X. Zotos, and P. Prelovsek, Phys. Rev. Lett. 74, 97 (1995). 9 Fabian H. L. Essler and Vladimir E. Korepin, Phys. Rev. Lett. 7, 908 (1994). 10 F. D. M. Haldane, Phys. Rev. Lett. 67, 937 (1991). 11 J. M. P. Carmelo and A. H. Castro Neto, to appear in Phys. Rev. B. 1 D. Baeriswyl, J. Carmelo, and A. Luther, Phys. Rev. B 33, 747 (1986); X. Zotos, P. Prelovsek, and I. Sega, ibid 4, 8445 (1990); H. Schulz, Phys. Rev. Lett. 64, 831 (1990). 13 W. Stephan and P. Horsch, Phys. Rev. B 4, 8736 (1990); P. Horsch and W. Stephan, ibid 48, (1993). 14 E. Y. Loh, JR. and D. K. Campbell, Synth. Metals 7, A499 (1988). 11