arxiv: v1 [hep-th] 10 Nov 2017

Transcription

1 Velocity operator approach to a Fermion system arxiv: v1 [hep-th] 10 ov 017 Seiya ISHIYAMA and João da PROVIDÊCIA Centro de Física, Departamento de Física, Universidade de Coimbra, P Coimbra, Portugal ovember 15, 017 Abstract In this paper, we formulate a velocity operator approach to a three-dimensional (3D Fermion system. Following Sunaawa, introducing density and velocity operators, we treat 3D quantum fluid dynamics in the system. We get a collective Hamiltonian in terms of collective variables. The lowest order collective Hamiltonian is diagonalized. This diagonalization leads us to a Bogoliubov transformation for Boson-lie operators. Keywords: Collective motion of a three-dimensional Fermion system; velocity operator; vortex motion; Grassmann numbers PACS umber (s: n, 1.60.Ev 1 Introduction To approach elementary excitations in a Fermi system, Tomonaga and Emery gave basic ideas in their collective motion theories [1, ]. On the other hand, Sunaawa s discrete integral equation method for a Fermi system [3] may be expected to also wor well for a collective motion problem. In the preceding papers [4, 5], introducing density operators ρ and associated variables π and defining exact momenta Π (collective variables, we could get an exact canonically conjugate momenta approach to one- and three-dimensional (1D and 3D Fermion systems. In the present paper, we formulate a velocity operator approach to a 3D Fermion system. Following Sunaawa, after introducing momentum density operators g, we define velocity operators v which denote classical fluid velocities. We derive a collective Hamiltonian in terms of the collective variables v and ρ for the elementary excitations including vortex motions. The lowest order collective Hamiltonian is diagonalized. This diagonalization leads to a Bogoliubov transformation for Boson-lie operators [6]. In firstweintroducecollectivevariablesρ andassociatedmomentumdensityvariables g and give the commutation relations between them. The velocity operator v is defined by a discrete integral equation and commutation relations between the velocity operators are also given. In 3 the dependence of the original Hamiltonian on v and ρ is determined. This section is also devoted to a calculation of a constant term in the collective Hamiltonian. Its lowest order is diagonalized and leads to a Bogoliubov transformation for Boson-lie operators. Finally in 4 some discussions and further outloo are given. 1 Corresponding author. address: seioceu@he.biglobe.ne.jp; nisiyama@teor.fis.uc.pt address: providencia@teor.fis.uc.pt 1

2 Collective variable and velocity operator Let us define the Fourier component of the density operator (ρ(x=ψ (xψ(x as ρ 1 p a p+ a p, ρ 0 = 1 p a p a p=, (: total number of Fermion. (.1 We here consider a spin-less Fermion and have used the anti commutation relations (CRs among a s and a s. The Hamiltonian H in a cubic box of volume Ω (=L3 is given as H=T+V = m a a + 0 Ω ν(ρ ρ 0 Ω ν(+( 1 Ω ν(0. (. TheFouriercomponent oftheinteractionv isgivenbytherelationv(x 1 Ω ν(ei x. Following Sunaawa [7], we introduce a momentum density operator given by g p pa p a p+, g 0 = p pa p a p=0. (.3 CRs among ρ and vector g (=(g (1,g(,g(3 with vector (=( 1,, 3 are obtained as [g,ρ ]= ρ, [g(i,g(j ]= ( g (i + j g (j + i. (.4 Following Sunaawa [7], we define the modified momentum density operator v by v g 1 p ρ p v p ( 0, v 0 =0. (.5 AlongthesamewayastheoneinIandIIandtheSunaawa s[3], weprovetwoimportantcrs [v,ρ ]= δ,, (.6 [v (i,v(j ] pall ρ p (p i v p (j p jv p (i + 1 p,q ρ p ρ q [v(i p,v(j q ]. (.7 The symbol means that the term p= and q= at the same time should be omitted. The CRs (.6 and (.7 are quite the same as the Sunaawa s [7]. As pointed out by him, Fourier transforms of the operators ρ and v and the CRs among them are identical with those found by Landau [8] for the fluid dynamical density operator and the velocity operator. Then, it turns out that the quantum mechanical operatorv(x, which satisfies the famous CR [v (i (x,v (j (x ]= i m δ(x x ρ(x 1 (rotv(x ( (i, j, cyclic, (.8 corresponds to the fluid dynamical velocity. We also have [v(x,ρ(x ]= i m xδ(x x. 3 v - and ρ -dependence of the Hamiltonian We derive here a collective Hamiltonian in terms of the v and ρ. Following Sunaawa, we expand the inetic operator T in a power series of the velocity operator v as follows: T=T 0 (ρ+ p 0 T 1(ρ;p vp+ p 0,q 0 T (ρ;p,qvp vq+, T (ρ;p,q=t (ρ;q,p, (3.1 in whicht n (n 0are unnown expansion coefficients. In order to determine their explicit expressions, we tae the CRs between T and ρ as follows: [T,ρ ]= T 1 (ρ; + p 0 T (ρ;p,vp +, [[T,ρ ],ρ ]= T (ρ;, + (3.. From (.3 and (.5, we can calculate the CRs between T and ρ as follows: [T,ρ ]= m g = m v + m p ρ p+ v p, [[T,ρ ],ρ ]= m δ,, (3.3 m ρ +,, [[[T,ρ ],ρ ],ρ ]=0,.

3 Comparing the above results with the CRs (3., we can determine the coefficients T n (n 0. Then we can express the inetic part T in terms of the ρ and v as follows: T =T 0 (ρ+ 1 m 0 v v + 1 m p+q 0 ρ p+qvp vq, (v 0 =0. (3.4 Up to the present stage, all the expressions have been derived without any approximation. To determine T 0 (ρ, we expand it inapower series of ρ and taethe CRs [T 0 (ρ,v (i ] as follows: T 0 (ρ=c 0 + p 0 C 1(pρp + p 0,q 0 C (p,qρpρq +, C (p,q=c (q,p, [v (i,t 0(ρ]= i C 1 (+ i p 0 C (p;ρp +, (3.5 [v (j,[v (i,t 0(ρ]]= i j C ( ;+6 i j p 0 C 3(p; ;ρp +. From (.5, we have a discrete integral equation [v (i,t 0(ρ]=[g (i,t 0(ρ] 1 ( p ρ p [v(i p,t 0(ρ], Denote [g (i,t 0(ρ] as f (i (ρ;.(3.6 With the aid of (3.4 and using two CRs of (.4, f (i (ρ; is calculated approximately as f (i (ρ; [g (i,t] 1 m p+q 0 [g(i,ρ p+q]vp vq 1 p 0 m [g(i,v p v p] =[g (i {,T] p,qall m ρ p+q ( v p } v (i q + ivp vq m { } (3.7 p 0 ( v pv (i +p +v(i +p ( v p p i v p v +p p iv +p v p =[g (i,t] Λ1(i Λ(i Λ3(i, Λ1(i p,qall m ρ ( p+q v p v (i q, where { ( [g (i 3,T]= m pall p i p+ ( p } a p a p+ = 3 m pall p i(p a p a p+, (3.8 Λ (i p,qall m ρ p+q ivq vp, (3.9 Λ 3(i m { p 0 ( v pv (i +p +v(i +p ( v p p i v p v +p p iv +p v p Λ 1(i = 3 }. (3.10 From now on, using ρ 0 =, we mae approximations for v,ρ,a 0 and a 0 as v = (θa a θ, ρ =θa +a θ, a 0= θ, a 0= θ, (3.11 where θ and θ are the Grassmann numbers [9]. Using ρ p+q = δ q, p, Λ1(i in (3.7 is computed as p all p( i p i (θa p (θap a a ( p θ p θ 4m 4m p all p ( i p i ρ pρ p 3 +θθ m p all {(p = 3 = 3 4m ( p i + i pall p ( i p i ρ pρ p +θθ 3 m 3 + (p + ( i 4 i ( p i i ρ, } a p a p+ (3.1

4 where the last term in the last line is obtained by extracting the terms withp= orp= in the θθ term of the third line. As is clear from its structure, it evidently vanishes. The remaining terms in the θθ term are small and can be neglected. The termsλ (i andλ3(i in (3.9 and (3.10, which have never been seen in the previous paper [5], are also calculated similarly as the above, respectively and are given as Λ (i = 3 i 4m p all (θa ( p p p (θap a a θ ( p p θ = 3 i 4m pall p (p ρ (3.13 i pρ p θθ3 3 16m ρ, Λ 3(i 3 = 4m p 0 {( i+p i p p i (+p p} (θa +p (θap a a (+p θ p θ (3.14 = 3 4m p 0 {( i+p i p p i (+p p}ρpρ. (p+ Here in (3.13 we use the relation θθ=1. Substituting [g (i,t] = 3 i 4m ρ andλ(i s, wegetanapproximateformulaforf(i (ρ;as f (i (ρ;= 7 3 i 16m ρ 3 4m = 7 3 i 4 4m ρ 3 4m p {p ( i p i i p (p ( i +p i p+p i (+p p}ρ pρ p (3.15 p {(p ( i p i p ( i p i }ρ pρ p. In the last line, the numerical factor 7 4 and the term p ( i p i, which have also never been seen in the previous paper [5], arise due to the consideration of the termsλ (i andλ3(i. Further substituting (3.15 into (3.6 and picing up the next leading term, we have [v (i,t 0(ρ] 7 3 i 16m ρ 3 4m p {(p p ( i p i + 74 } p p i ρ pρ p. (3.16 To eep a favored form p (p ( i p i ρ pρ p occurred in [5], in the curly bracets, we had better to treat i p i and 7 4 p i equally and to lead p (p+. So, we tae p i = 8 15 i and then i p i = 7 15 i. Thus, we have a final expression for [v (i,t 0(ρ] up to the order of 1 as [v (i,t 0(ρ]= 7 3 i 4 4m ρ 14 3 i 158m p p (p+ρ pρ p. (3.17 From (3.16 and (.6, we get the following CRs between v (i and T 0(ρ: 4 i j 4m ( + + ρ,,t 0(ρ]]]= 14 5 i j ( m ( + +, [v (j,[v (i,t 0(ρ]]= 7 4 i j 4 4m δ, [v (,[v (j,[v (i [v (l,[v (,[v (j,[v (i,t 0(ρ]]]]=0. 4 (3.18

5 By the procedure similar to the previous one, we can determine the coefficients C and C 3 in (3.5 as C ( ;= 7 4 8m δ, and C 3( ; ;= 7 ( 151m + + δ, which are the same forms as the Sunaawa s C and C 3 [7] but having different numerical factors. This is also due to the consideration of the termsλ (i andλ3(i. Then get an approximate form of T 0 (ρ in terms of variables ρ as T 0 (ρ=c m 0 ρ ρ m p 0,q 0,p+q 0 (p +p q+q ρpρqρ p q.(3.19 Substituting (3.19 into (3.4 and using the underlying identity, p 0,q 0,p+q 0 (p +p q+q ρpρqρ p q= 3 p 0,q 0,p+q 0 p qρ pρqρ p q, (3.0 the constant term C 0 is computed as C 0 = T 7 48m 0 ρ ρ 1 m v v 1 m p+q 0 ρ p+qvp vq 14 ( m p 0,q 0,p+q 0 p qρ pρqρ p q. Using (3.11, we can calculate the third term in (3.1 and then reach to a result such as 1 m v v = 8m (θa a (θa θ a θ (3. = 8m ρ ρ +θθ m θθ m a a. As for the forth and last terms in (3.1, they become vanishing due to θθ=0 and θθ=0. Substituting these into (3.1 and replacing ρ ρ by <ρ ρ > Ave, we have C 0 = ( 1 θθ T+θθ m 3 3m ρ ρ = m 3 3m <ρ ρ > Ave, (3.3 where the quantity <ρ ρ > Ave stands for the average value of ρ ρ and we have used the relation θθ=1. It is surprising to see that the C 0 coincides with the constant term in the ground-state energy given by Tomonaga [1] except the last term. Using (3.4, (3.19 and (3.3, we can express the Hamiltonian H(= H I +H II in terms of ρ and v as follows: [ H I = ( 1 7 Ω ν( m 7 Ω ν(+ 1 ( 7 m v v + ] 3m + Ω ν( ρ ρ ( ( 3 7 3m <ρ ρ > Ave m Ω ν(, (3.4 0 H II = 1 m p 0,q 0,p+q 0 ρ p+qvp vq m p 0,q 0,p+q 0 p qρ pρ p qρq, which is similar to the Sunaawa s Hamiltonian [3] except the term <ρ ρ > Ave in H I. ow, let us introduce the Boson annihilation and creation operators defined as α me ρ + 1 m E v, α me ρ + 1 m E v, ( 0. (3.5 5

6 Using(3.5, λ and(3.11, thecollectivevariablesρ me andv ( 0areexpressedas ρ = ( λ α =θa +α +a θ, v = (α λ α = (θa a θ. (3.6 The H I in (3.4 is considered as the lowest order Hamiltonian within the scope of the present approximation. With the use of an inverse transformation of (3.5, the lowest order Hamiltonian H I denoted as H 0 is diagonalized as H 0 =E0 G + E 0 α α, E G 0 E ν( ( 1 Ω 7 4 ε + m Ω ν(, ε m, ( 7 7 E 4m 7Ω ν( 3 16m <ρ ρ > Ave where E is the quasi-particle energy but with the modified single-particle energy, (3.7 7 ε in the lowest approximation and we have used [v,ρ ]=. In this sense, the zero point energy of the collective mode is included in the above diagonalization. The quantity E0 G in (3.7 corresponds to the lowest order ground-state energy [10, 6] but with the average value <ρ ρ > Ave. The part of the Hamiltonian, H II, is also represented in terms of the Boson operators α and α though we omit here its explicit expression. Further we have a Bogoliubov transformation for Boson-lie operators θa and a θ as α = ( E +ε θa + ( E ε a θ ε E, α = ( E ε θa + ( E +ε a θ ε E, (3.8 which already appeared in the preceding papers [4, 5] and is of the same form as the Bogoliubov transformation for the usual Bosons [6] except the introduction of the Grassmann numbers θ and θ. The above ind of the diagonalization (3.7 has also been given by Sunaawa in the Boson system [11]. 4 Discussions and further outloo In this paper, we have proposed a velocity operator approach to a 3D Fermion system. After introducing collective variables, the velocity operator approach to the 3D Fermion system could be provided. Particularly, an interesting problem of describing elementary excitations occurring in a strongly correlated system may be possibly treated as an elementary exercise. For such problems, for example, see textboos [1, 13]. By applying the velocity operator approach to such a problem, an excellent description of the elementary excitations in the 3D Fermion system will be expected to reproduce various possible behaviors including excited energies. Because the present theory is constructed to tae into account important many-body correlations which have not been investigated sufficiently for a long time in the historical ways for such a problem. Here the velocity operator v (i is defined through the discrete integral equation which is necessarily accompanied by the inhomogeneous term f (i (ρ;. It is a very important tas how to evaluate this term which brings an essential effect to the construction of the theory for the velocity operator approach. More precisely, in the last line 6

7 of (3.7, we have taen into account the contributions from all of the second term Λ 1(i, the third one Λ (i and the forth one Λ3(i. Such an evaluation is quite different from the manner adopted in the preceding papers [4, 5] in which only the term Λ 1(i had been taen into account. As a result, the present evaluation leads to another approximate formula for f (i (ρ; and another final form of [v (i,t 0(ρ] which are of course different from those given in those papers. Under the use of the CRs (3.18 between v (i and T 0(ρ, thus, the correct expression for T 0 (ρ could be determined in terms of ρ and the constant C 0 appearing in the expansion of a modified T 0 (ρ turns out to involve the term <ρ ρ > Ave which is not determined yet. However, now, maing use of the inverse transformation of (3.5, we determine a notable term as <ρ ρ > Ave= me <(α +α (α +α > Ave= me = ε E. Using the last two energies E and ε in (3.7, here we approximate E and the term 3 16m <ρ ρ > Ave in 7 7 E0 G of (3.7 as ε and 3 8 ε, respectively. Then the lowest order ground-state energy ( E0 G is rewritten ase0= G ( 1 Ω ν(0+1 E ε 7 Ω ν(. Further we derive thebogoliubovtransformationforboson-lie operatorsθa anda θ whichalreadyappeared in the preceding papers [4, 5] and is the same form as the famous Bogoliubov transformation 7 for the usual Bosons [7] but which brings the modified single-particle energy ε in the quasi-particle excitation energy E in the lowest order approximation. Thus, we could provide the modified and then the correct theory for the velocity operator approach. Based on such approach, it is expected that a new field of exploration of elementary excitation in a 3D Fermi system may open up with aid of the velocity operator approach whose new development may appear elsewhere. The connection of the present theory with the fluid dynamics was been mentioned briefly by Sunaawa. Restrict the Hilbert space to a subspace in which the vortex operator satisfies rotv(x >=0, ( v >=0 leading effectively to [v (i,v(j ]=0. He transformed the quantum-fluid Hamiltonian (3.4 to one in configuration space and obtained the classical-fluid Hamiltonian for the case of irrotational flow [7, 11]. From this view point, very recently, we have developed a quantum-fluid approach to vortex motion in nuclei [14]. In this wor we have introduced Clebsch parameterization maing it possible to describe canonically quantum-fluid dynamics and have done Clebsch-parameterized gauge potential possessing Chern-Simons number (quantized helicity [15, 16, 17]. A study of Chern-Simons form shed new light on the aspect of the behavior of the quantum fluid. Acnowledgements One of the authors (S.. would lie to express his sincere thans to Professor Constança Providência for ind and warm hospitality extended to him at the Centro de Física, Universidade de Coimbra. This wor was supported by FCT (Portugal under the project CER/FP/83505/008. The authors than the Yuawa Institute for Theoretical Physics at Kyoto University. Discussions during the worshop YITP-W on Strings and Fields 017 are useful to complete this wor. 7

8 References [1] S. Tomonaga, Prog. Theor. Phys. 5, 544 (1950. [] V.J. Emery, Theory of the One-Dimensional Electron Gas in Highly Conducting One- Dimensional Solids, Editors, J.T. Devreese, R.P. Evrard and V.E. van Dore, Physics of Solids and Liquids, Springer US, 1979, pp [3] S. Sunaawa, Y. Yoo-o and H. aatani, Prog. Theor. Phys. 7, 589 (196. [4] S. ishiyama and J. da Providência, Int. J. Mod. Phys. E4, (015. [5] S. ishiyama and J. da Providência, Int. J. Mod. Phys. E5, (016. [6].. Bogoliubov, J. Phys. 11, 3 (1947. [7] S. Sunaawa, Y. Yoo-o and H. aatani, Prog. Theor. Phys. 7, 600 (196. [8] L. Landau, J. Phys. 5, 71 (1941. [9] F.A. Berezin, The Method of Second Quantization, Academic Press, ew Yor and London, [10].. Bogoliubov and D.. Zubarev, Soviet Phys. JETP 1, 83 (1955. [11] S. Sunaawa, Y. Yoo-o and H. aatani, Prog. Theor. Phys. 8, 17 (196. [1] P. Ring and P. Schuc, The nuclear many body problem, Texts and monographs in Physics (Springer-Verlag, Berlin, Heidelberg and ew Yor [13] J.W. egele and H. Orland, Quantum Many-Particle Systems, Frontiers in Physics, Volume 68, Addison-Wesley Publishing Company, [14] S. ishiyama and J. da Providência, Int. J. Mod. Phys. E6, (017. [15] A. Jaciw and S.-Y. Pi, Phys. Rev. D61, (000. [16] A. Jaciw, V.P. air and S.-Y. Pi, Phys. Rev. D6, (000. [17] A. Jaciw, CRM Series in Mathematical Physics, Lectures on fluid dynamics, A Particle Theorist s View of Supersymmetric, on-abelian, oncommutative Fluid Mechanics and d-branes, 00 Springer-Verlag, ew Yor, Inc. 8