Non-relativistic expansion of single-nucleon Dirac equation: Comparison between Foldy-Wouthuysen transformation andsimilarity renormalization group

  • By following the Foldy-Wouthuysen (FW) transformation of the Dirac equation, we derive the exact analytic expression up to the 1/M4 order for general cases in the covariant density functional theory. The results are compared with the corresponding ones derived from another novel non-relativistic expansion method, the similarity renormalization group (SRG). Based on this comparison, the origin of the difference between the results obtained with the FW transformation and the SRG method is explored.
  • The Foldy-Wouthuysen (FW) transformation, also called the Pryce-Tani-Foldy-Wouthuysen transformation [14], was initially formulated around 1950 to describe the non-relativistic limit of relativistic electrons. A great advantage of the FW transformation is the simple form of operators corresponding to classical observables. This transformation for spin-1/2 particles was soon extended to the cases of spin-0 and spin-1 particles [5], and later generalized to the case of arbitrary spins [6]. Recently, the exact form of the exponential FW operator for a particle with an arbitrary spin was given in Ref. [7].

    The FW transformation is also an elegant method for non-relativistic expansion of the Dirac equation [8]. One milestone of its applications in atomic physics was that, using the FW transformation in the presence of external fields, the fine structure [9, 10], describing the splitting of spectral lines of atoms due to the electron spin and the relativistic corrections to the non-relativistic Schrödinger equation, was successfully testified [8].

    Due to its elegant and powerful scheme, the FW transformation has been applied in very diverse areas, such as electrodynamics [11, 12], atomic physics [13, 14], condensed matter physics [1517], optics [1820], acoustics [21], quantum chemistry [2227], nuclear physics [28, 29], quantum field theory [30], gravity [3134], in the theory of the weak interaction [35, 36], and even in the scheme of supersymmetric quantum mechanics [37]. In recent years, the FW transformation has also been discussed in a gravitational wave background of the generalized Dirac Hamiltonian in the non-relativistic limit [3840].

    Nevertheless, in most of the above studies and applications of the FW transformation, only the leading-order terms were derived and treated, such as the relativistic correction to the kinetic energy appearing up to the order of 1/M3 (with M the bare mass of the particle), and the correction due to the spin-orbit coupling and the Darwin term appearing in the order of 1/M2 [8]. Recently, the FW non-relativistic expansions have been carried out up to the order of 1/M7 [41], and even 1/M14 [42]. Nevertheless, such derivations are valid only for a family of potentials with certain properties.

    For systems of atomic nuclei, the covariant density functional theory [4345] is one of the state-of-the-art methodologies. During the past decades, it has achieved great success in describing various nuclear phenomena in both stable and exotic nuclei. In this theoretical framework, the equation of motion for nucleons is described by the Dirac equation which includes not only the strong vector potential V(r) , but also the strong scalar potential S(r), whose typical values are V(r)300 MeV and S(r)350 MeV, while the mass of nucleon is M=939 MeV. As a result, the non-relativistic expansion of the single-nucleon Dirac equation is non-trivial, and its high-order terms have to be derived and verified carefully.

    In 2012, Guo [46] applied the similarity renormalization group (SRG) [47, 48], instead of the FW transformation, for the non-relativistic expansion of the single-nucleon Dirac equation up to the 1/M3 order for investigating the nuclear pseudospin symmetry [49]. Using the SRG method, the Dirac Hamiltonian is gradually transformed into a diagonal form with the evolution controlled by the flow equation. As a result, the eigenequations for the upper and lower components of the Dirac spinors are decoupled at the end of the flow. Also, the non-relativistic reduced Hamiltonian thus obtained is Hermitian, and can also be expanded into the series in terms of 1/Mi. To achieve the accuracy of single-particle energies at the 0.1 MeV level, we demonstrated recently that the 1/M4-order terms are needed, and presented a detailed derivation and results in Ref. [50].

    Both the FW transformation and the SRG method provide a systematic way to carry out a non-relativistic expansion of the Dirac equation up to an arbitrary order. However, the FW transformation employs finite steps of unitary transformation, while the SRG method employs infinite infinitesimal unitary transformations along the flow. Therefore, it is interesting and important to investigate the similarities and differences between these two approaches.

    In this paper, we perform an exact non-relativistic expansion of the single-nucleon Dirac equation up to the 1/M4 order by using the FW transformation. We also investigate the difference between such results and those obtained by the SRG method.

    The paper is organized as follows. In Section 2.1, the SRG method is first recalled and the expansion up to the 1/M4 order is shown. The new expansion up to the same order using the FW transformation in a general case, such as the inclusion of the scalar potential, is derived in Section 2.2. The specific form of the result in spherical symmetry is presented in Section 3.1, and the difference between the SRG and FW results, as well as the origin of the difference, is investigated in Section 3.2. Finally, a summary and perspectives are given in Section 4.

    In the theoretical framework of covariant density functional theory [4345, 5156], the Dirac Hamiltonian for nucleons reads

    H=αp+β(M+S)+V,

    (1)

    where α and β are the Dirac matrices, M is the mass of the nucleon, and S and V are the scalar and vector potentials, respectively. In this section, we recall the technique of SRG [47, 48] and show the results up to the 1/M4 order [46, 50].

    The basic idea of SRG is that the Dirac Hamiltonian in Eq. (1) is transformed by a unitary operator U(l) as

    H(l)=U(l)HU(l),H(0)=H,

    (2)

    with a flow parameter l. While U(l) is not known explicitly, an anti-Hermitian generator η(l) related to U(l) as η(l)=dU(l)dlU(l) is chosen, and the Hamiltonian H(l) evolves with the differential flow equation

    dH(l)dl=[η(l),H(l)].

    (3)

    One of the convenient and appropriate choices of η(l) reads [48]

    η(l)=[βM,H(l)],

    (4)

    which transforms the Dirac Hamiltonian into a diagonal form.

    For solving Eq. (3), the Hamiltonian is first divided into two parts:

    H(l)=ε(l)+o(l),

    (5)

    according to the commutation and anti-commutation relations with respect to β, i.e., [ε,β]=0 and {o,β}=0. As a result,

    dε(l)dl=4Mβo2(l),

    (6a)

    do(l)dl=2Mβ[o(l),ε(l)],

    (6b)

    with the initial conditions

    ε(0)=β(M+S)+V,o(0)=αp.

    (7)

    Equations (6a) and (6b) can then be solved with the expansion into the series of 1/Mi [48], which gives

    ε(λ)M=i=0εi(λ)Mi,

    (8a)

    o(λ)M=j=1oj(λ)Mj,

    (8b)

    while the flow parameter is replaced by a dimensionless parameter λ=lM2. The solutions are obtained as [46]

    εn(λ)=εn(0)+4βλ0n1k=1ok(λ)onk(λ)dλ,

    (9a)

    on(λ)=on(0)e4λ+2βe4λλ0n1k=1[e4λok(λ),εnk(λ)]dλ,

    (9b)

    with the initial conditions,

    ε0(0)=β,ε1(0)=βS+V,εn(0)=0ifn2,o1(0)=αp,on(0)=0ifn2.

    (10)

    Therefore, at the end of the flow λ, all the off-diagonal parts vanish, i.e., on()=0 for all n. The diagonalized Dirac operator up to the 1/M4 order is not difficult to obtain as [50]

    HSRG=ε()=Mε0()+ε1()+ε2()M+ε3()M2+ε4()M3+ε5()M4+=Mε0(0)+ε1(0)+12Mβo21(0)+18M2[[o1(0),ε1(0)],o1(0)]+132M3β(4o41(0)+[[o1(0),ε1(0)],ε1(0)]o1(0)+o1(0)[[o1(0),ε1(0)],ε1(0)]2[o1(0),ε1(0)][o1(0),ε1(0)])+1128M4(9[[o1(0),ε1(0)],o31(0)]+3[o1(0),ε1(0)]2ε1(0)+3ε1(0)[o1(0),ε1(0)]26[o1(0),ε1(0)]ε1(0)[o1(0),ε1(0)]+3[o1(0)[o1(0),ε1(0)]o1(0),o1(0)+[[[[o1(0),ε1(0)],ε1(0)],ε1(0)],o1(0)])+

    (11)

    In this section, we perform the exact non-relativistic expansion of the single-nucleon Dirac equation up to the 1/M4 order by using the FW transformation for a general case, in particular for the scalar potential.

    According to the FW transformation in the presence of external fields [2, 4, 8], the corresponding operators are defined as

    O=αp,ε=βS+V,Λ=i2MβO,

    (12)

    where the operators O and ε satisfy Oβ=βO and εβ=βε, respectively. The Hamiltonian (1) reads

    H=βM+O+ε.

    (13)

    Based on the above definitions, the Dirac Hamiltonian is transformed by a unitary operation into [8]

    H=eiΛHeiΛ=H+i[Λ,H]+i22![Λ,[Λ,H]]++inn![Λ,[Λ,,[Λn,H]]]+

    (14)

    where

    inn![Λ,[Λ,,[Λn,H]]]=(1)n(n1)2βnn!Mn(OnβM+On+1+12n[O,[O,,[On,ε]]]).

    (15)

    Keeping all the terms up to the 1/Mn order, the unitary transformed Hamiltonian reads

    H1/Mn=βM+ε+n+1k=1(1)(k1)(k2)2k1k!Mk1βk1Ok+nk=0(1)k(k1)2βk2kk!Mk[O,[O,,[Ok,ε]]].

    (16)

    For example, the corresponding result up to the 1/M4 order is

    H1/M4=βM+ε+12MβO213M2O318M3βO4+130M4O5+12Mβ[O,ε]18M2[O,[O,ε]]148M3β[O,[O,[O,ε]]]+1384M4[O,[O,[O,[O,ε]]]].

    (17)

    In the unitary transformed Hamiltonian, e.g., Eq. (17), the off-diagonal parts are not zero but raised by one order from O to 12Mβ[O,ε] (plus higher-order terms). In order to make the off-diagonal parts vanish, more precisely speaking to make them higher than a given order, one should repeat the FW transformation until the accuracy is achieved [8]. For that purpose, the operators O and ε can be redefined according to Eq. (17), and one has

    ε=ε+12MβO218M2[O,[O,ε]]18M3βO4+1384M4[O,[O,[O,[O,ε]]]],

    (18a)

    O=12Mβ[O,ε]13M2O3148M3β[O,[O,[O,ε]]]+130M4O5,

    (18b)

    Λ=i2MβO.

    (18c)

    The corresponding FW transformation reads

    H=eiΛHeiΛ.

    (19)

    It can be seen that the off-diagonal parts in H are of the order of 1/M2. Therefore, one should repeat this procedure to H′′′′′ to make the off-diagonal parts of the order of 1/M5. Of course, to get the resultant diagonal parts in H′′′′′ up to the 1/M4 order, a recursion technique makes the calculation less complicated than it looks.

    As a result, the non-relativistic expansion with the FW transformation up to the 1/M4 order reads

    HFW=βM+ε+12MβO218M2[O,[O,ε]]18M3βO418M3β[O,ε][O,ε]+1384M4[O,[O,[O,[O,ε]]]]+112M4[O3,[O,ε]]+132M4[O,ε][O,ε]ε116M4[O,ε]ε[O,ε]+132M4ε[O,ε][O,ε].

    (20)

    For the systems with spherical symmetry, i.e., for spherical nuclei, the corresponding radial single-nucleon Dirac equation reads [4345]

    (Σ(r)+Mddr+κrddr+κrΔ(r)M)(G(r)F(r))=E(G(r)F(r)),

    (21)

    where κ is a good quantum number defined as κ=(j+1/2) for j=l±1/2, and Σ(r)=V(r)+S(r) and Δ(r)=V(r)S(r) are the sum and the difference of the vector and scalar potentials, respectively. The single-particle energy E=ε+M includes the mass of the nucleon. The operators ε and O read

    ε=(Σ(r)00Δ(r)),O=(0ddr+κrddr+κr0).

    (22)

    According to Eq. (20), the Dirac Hamiltonian is transformed by the FW transformation as

    (H(F)FW+MO(1M5)O(1M5)H(D)FWM),

    (23)

    where the superscripts (F) and (D) represent the components of the Dirac Hamiltonian acting on single-particle states in the Fermi sea and Dirac sea, respectively. It is seen that the off-diagonal parts are not strictly zero, which is different from the results of SRG, but are of higher order than the required one. Focusing on single-particle states in the Fermi sea which correspond to their counterparts in the non-relativistic framework, the explicit expansion of H(F)FW up to the 1/M4 order is worked out in detail as

    H(F)0,FW=Σ(r),

    (24a)

    H(F)1,FW=12Mp2,

    (24b)

    H(F)2,FW=18M2(4Sp2+4Sddr2κrΔ+Σ),

    (24c)

    H(F)3,FW=18M3(p4+4S2p28SSddr2SΣ+4SΔκr+ΣΔ),

    (24d)

    H(F)4,FW=1384M4{144Sp4288Sp2ddr+[72Δκr24(4Σ+9S)192S3]p2+[72Δκr272Δκr+24(3S+4Σ)+576S2S]ddr+[12(5Σ36S)κ(κ+1)r3+12(5Σ+12S)κ(κ+1)r272Δκr3+72Δκr236Δκr288S2Δκr+33Σ+144S2Σ+24SΣ(2Σ6Δ)]},

    (24e)

    where

    p2=d2dr2+κ(κ+1)r2,

    (25)

    and

    p4=d4dr42κ(κ+1)r2d2dr2+4κ(κ+1)r3ddr+κ(κ+1)(κ+3)(κ2)r4.

    (26)

    In contrast, the Dirac Hamiltonian transformed by the SRG method reads

    (H(F)SRG+M00H(D)SRGM).

    (27)

    Its off-diagonal parts strictly vanish up to infinite order. According to Eqs. (11), the explicit expansion of H(F)SRG up to the 1/M4 order reads [50]

    H(F)0,SRG=Σ(r),

    (28a)

    H(F)1,SRG=12Mp2,

    (28b)

    H(F)2,SRG=18M2(4Sp2+4Sddr2κrΔ+Σ),

    (28c)

    H(F)3,SRG=132M3(4p4+16S2p232SSddr8SΣ+16SΔκr2Σ2+4ΣΔ),

    (28d)

    H(F)4,SRG=1128M4{48Sp496Sp2ddr+[24Δκr24(Σ+3S)64S3]p2+[24Δκr224Δκr+24(S+Σ)+192S2S]ddr

    +[12(Σ12S)κ(κ+1)r3+12(Σ+4S)κ(κ+1)r224Δκr3+24Δκr212Δκr96S2Δκr+9Σ+48S2Σ+24SΣ(Σ2Δ)]}.

    (28e)

    By comparing Eqs. (11) and (20), it is found that

    H0,SRGH0,FW=0,

    (29a)

    H1,SRGH1,FW=0,

    (29b)

    H2,SRGH2,FW=0,

    (29c)

    H3,SRGH3,FW=β32M3[[O2,ε],ε],

    (29d)

    H4,SRGH4,FW=164M4[[O2,ε],O2]+1128M4[([ε2,O2]2[ε,OεO]),ε].

    (29e)

    This result seems to indicate that the FW transformation and the SRG method agree with each other up to the 1/M2 order, but they lead to different results starting from the 1/M3 order. However, considering the infinite mass limit M, all expressions should be strictly organized order by order, and thus the difference between Eqs. (27) and (23) in the 1/M5 order cannot lead to the differences of the results in the 1/M3 and 1/M4 orders. This puzzle needs to be further investigated.

    It turns out that the difference shown in Eq. (29) comes from an additional unitary transformation, after the Dirac Hamiltonian is decoupled into the upper and lower parts. Let

    Ξ=iβ32M3[O2,ε]i128M4([ε2,O2]2[ε,OεO]).

    (30)

    This is a Hermitian and diagonal operator, i.e., Ξ=Ξ and βΞ=Ξβ. Acting an additional unitary transformation on HFW, it reads

    eiΞHFWeiΞ=HFW+i[Ξ,HFW]+

    (31)

    Keeping all terms up to the 1/M4 order, the result is

    eiΞHFWeiΞ=HFW+β32M3[[O2,ε],ε]+164M4[[O2,ε],O2]+1128M4[([ε2,O2]2[ε,OεO]),ε],

    (32)

    which is nothing else but HSRG.

    In the spherical case, the explicit expressions for operators β32M3[[O2,ε],ε], 164M4[[O2,ε],O2], and 1128M4[([ε2,O2]2[ε,OεO]),ε] , acting on single-particle states in the Fermi sea, respectively read

    116M3Σ2,

    (33a)

    164M4[4Σp24Σddr+4Σκ(κ+1)r34Σκ(κ+1)r2Σ],

    (33b)

    116M4SΣ2.

    (33c)

    These equations explain the difference between Eqs. (24d) and (28d), as well as between Eqs. (24h) and (28h).

    Since eiΞ is an additional unitary transformation acting on the already decoupled Hamiltonian, it does not affect the non-relativistic expansion of the Dirac equation, and the single-particle spectra obtained by the FW transformation and the SRG method are the same.

    The unitary transformation eiΞ is equivalent to a rotation in a linear and closed complex space, and in quantum mechanics, the operation or rotation in such a space should be associated with a corresponding generator. However, in the present context, it is nontrivial to write down a general expression for such a generator, or derive a conserved quantity that is asymptotically related to a generator up to an arbitrary order. This is because the number of unitary transformations from the FW results to the SRG ones would in principle be infinite, since the SRG method employs infinite infinitesimal unitary transformations along the flow. Therefore, the physics behind the introduced unitary transformation is open for future studies.

    In this work, we investigated the non-relativistic expansion of the single-nucleon Dirac equation in the theoretical framework of covariant density functional theory for a general case where the scalar potential is included. We worked out the exact analytical expansion up to the 1/M4 order using the FW transformation. A further investigation of the difference between Eqs. (27) and (23) explained the puzzle that a disagreement seems to appear between the results obtained with the FW transformation and the SRG method, i.e. the non-relativistic expansion of the Dirac equation is affected, but the single-particle spectra obtained by the FW transformation and the SRG method are the same.

    The non-relativistic expansion of the Dirac equation for nuclear systems using the FW transformation has been extended to as high orders as needed. Similarly to the SRG method, we anticipate that this novel non-relativistic expansion method will establish in future studies a potential bridge between the relativistic and non-relativistic density functional theories. In particular, with the present FW transformation, all unitary transformations involved are in an explicit form. This leads in a straightforward way to the corresponding transformations of the other relevant operators, such as the one-body density operators.

    [1] M. H. L. Pryce and S. Chapman, Proc. Royal Soc. A, 195: 62 (1948) doi: 10.1098/rspa.1948.0103
    [2] L. L. Foldy and S. A. Wouthuysen, Phys. Rev., 78: 29 (1950) doi: 10.1103/PhysRev.78.29
    [3] S. Tani, Prog. Theor. Phys., 6: 267 (1951) doi: 10.1143/ptp/6.3.267
    [4] L. L. Foldy, Phys. Rev., 87: 688 (1952) doi: 10.1103/PhysRev.87.688
    [5] K. M. Case, Phys. Rev., 95: 1323 (1954) doi: 10.1103/PhysRev.95.1323
    [6] J. Jayaraman, J. Phys. A: Math. Gen., 8: L1 (1975) doi: 10.1088/0305-4470/8/1/001
    [7] A. J. Silenko, Phys. Rev. A, 94: 032104 (2016) doi: 10.1103/PhysRevA.94.032104
    [8] W. Greiner, Relativistic Quantum Mechanics Wave Equations (Springer, 1990) pp. 277−290
    [9] A. A. Michelson and E. W. Morley, Philos. Mag. (Series 5), 24: 463 (1887) doi: 10.1080/14786448708628131
    [10] A. Sommerfeld, Sci. Nat., 28: 417 (1940) doi: 10.1007/BF01490583
    [11] A. Y. Silenko, Theor. Math. Phys., 105: 1224 (1995) doi: 10.1007/BF02067491
    [12] A. J. Silenko, J. Math. Phys., 44: 2952 (2003) doi: 10.1063/1.1579991
    [13] M. Hiramoto, T. Fujita, and T. Asaga, Prog. Theor. Phys., 106: 1223 (2001) doi: 10.1143/PTP.106.1223
    [14] K. Pachucki, Phys. Rev. A, 71: 012503 (2005) doi: 10.1103/PhysRevA.71.012503
    [15] K. Y. Bliokh, EPL, 72: 7 (2005) doi: 10.1209/epl/i2005-10205-1
    [16] P. Gosselin, A. Bérard, and H. Mohrbach, Eur. Phys. J. B, 58: 137 (2007) doi: 10.1140/epjb/e2007-00212-6
    [17] P. Gosselin, A. Bérard, H. Mohrbach, and S. Ghosh, Eur. Phys. J. C, 59: 883 (2009) doi: 10.1140/epjc/s10052-008-0839-4
    [18] M. Lippert, T. Bruckel, T. Kohler, and J. R. Schneider, EPL, 27: 537 (1994) doi: 10.1209/0295-5075/27/7/008
    [19] S. A. Khan, Optik, 117: 481 (2006) doi: 10.1016/j.ijleo.2005.11.010
    [20] Y. Hinschberger and P.-A. Hervieux, Phys. Lett. A, 376: 813 (2012) doi: 10.1016/j.physleta.2012.01.023
    [21] L. Fishman, J. Math. Phys., 33: 1887 (1992) doi: 10.1063/1.529666
    [22] M. Reiher, Theor. Chem. Acc., 116: 241 (2006) doi: 10.1007/s00214-005-0003-2
    [23] J. Autschbach, Coord. Chem. Rev., 251: 1796 (2007) doi: 10.1016/j.ccr.2007.02.012
    [24] W. Liu, Mol. Phys., 108: 1679 (2010) doi: 10.1080/00268971003781571
    [25] T. Nakajima and K. Hirao, Chem. Rev., 112: 385 (2012) doi: 10.1021/cr200040s
    [26] D. Peng and M. Reiher, J. Chem. Phys., 136: 244108 (2012) doi: 10.1063/1.4729788
    [27] D. Peng and M. Reiher, Theor. Chem. Acc., 131: 1081 (2012) doi: 10.1007/s00214-011-1081-y
    [28] H. Arodź and K. Golec, Phys. Rev. D, 33: 534 (1986) doi: 10.1103/PhysRevD.33.534
    [29] S. Scherer, G. Poulis, and H. Fearing, Nucl. Phys. A, 570: 686 (1994) doi: 10.1016/0375-9474(94)90079-5
    [30] V. P. Neznamov, Phys. Part. Nucl., 37: 86 (2006) doi: 10.1134/S1063779606010023
    [31] A. J. Silenko and O. V. Teryaev, Phys. Rev. D, 71: 064016 (2005) doi: 10.1103/PhysRevD.71.064016
    [32] P. Gosselin and H. Mohrbach, Eur. Phys. J. C, 71: 1739 (2011) doi: 10.1140/epjc/s10052-011-1739-6
    [33] Y. N. Obukhov, A. J. Silenko, and O. V. Teryaev, Phys. Rev. D, 88: 084014 (2013) doi: 10.1103/PhysRevD.88.084014
    [34] Y. N. Obukhov, A. J. Silenko, and O. V. Teryaev, Phys. Rev. D, 90: 124068 (2014) doi: 10.1103/PhysRevD.90.124068
    [35] A. J. Silenko, Nucl. Instrum. Methods Phys. Res. Sect. B, 114: 259 (1996) doi: 10.1016/0168-583X(96)00212-1
    [36] A. J. Silenko, Theor. Math. Phys., 112: 922 (1997) doi: 10.1007/BF02634108
    [37] G. Junker and A. Inomata, J. Math. Phys., 59: 052301 (2018) doi: 10.1063/1.5020545
    [38] B. Gonçalves, Y. N. Obukhov, and I. L. Shapiro, Phys. Rev. D, 75: 124023 (2007) doi: 10.1103/PhysRevD.75.124023
    [39] J. Q. Quach, Phys. Rev. D, 92: 084047 (2015) doi: 10.1103/PhysRevD.92.084047
    [40] J. Q. Quach, Phys. Rev. D, 93: 104048 (2016) doi: 10.1103/PhysRevD.93.104048
    [41] T.-W. Chen and D.-W. Chiou, Phys. Rev. A, 82: 012115 (2010) doi: 10.1103/PhysRevA.82.012115
    [42] T.-W. Chen and D.-W. Chiou, Phys. Rev. A, 90: 012112 (2014) doi: 10.1103/PhysRevA.90.012112
    [43] J. Meng, H. Toki, S. Zhou et al, Prog. Part. Nucl. Phys., 57: 470 (2006) doi: 10.1016/j.ppnp.2005.06.001
    [44] J. Meng and S. G. Zhou, J. Phys. G: Nucl. Part. Phys., 42: 093101 (2015) doi: 10.1088/0954-3899/42/9/093101
    [45] J. Meng, ed., Relativistic Density Functional for Nuclear Structure (World Scientific, 2016)
    [46] J.-Y. Guo, Phys. Rev. C, 85: 021302 (2012) doi: 10.1103/PhysRevC.85.021302
    [47] F. Wegner, Ann. Phys. (Berl.), 506: 77 (1994) doi: 10.1002/andp.19945060203
    [48] A. Bylev and H. Pirner, Phys. Lett. B, 428: 329 (1998) doi: 10.1016/S0370-2693(98)00417-1
    [49] H. Liang, J. Meng, and S.-G. Zhou, Phys. Rep., 570: 1 (2015) doi: 10.1016/j.physrep.2014.12.005
    [50] Y. Guo and H. Liang, Phys. Rev. C, 99: 054324 (2019) doi: 10.1103/PhysRevC.99.054324
    [51] Y. Chen, L. Li, H. Liang et al, Phys. Rev. C, 85: 067301 (2012) doi: 10.1103/PhysRevC.85.067301
    [52] W. Zhang and Y.-F. Niu, Chin. Phys. C, 41: 094102 (2017) doi: 10.1088/1674-1137/41/9/094102
    [53] Z. Xu and Z.-P. Li, Chin. Phys. C, 41: 124107 (2017) doi: 10.1088/1674-1137/41/12/124107
    [54] Y.-Y. Zong and B.-Y. Sun, Chin. Phys. C, 42: 024101 (2018) doi: 10.1088/1674-1137/42/2/024101
    [55] R. Swain, S. K. Patra, and B. B. Sahu, Chin. Phys. C, 42: 084102 (2018) doi: 10.1088/1674-1137/42/8/084102
    [56] Z.-Z. Li, S.-Y. Chang, Q. Zhao et al, Chin. Phys. C, 43: 074107 (2019) doi: 10.1088/1674-1137/43/7/074107
  • [1] M. H. L. Pryce and S. Chapman, Proc. Royal Soc. A, 195: 62 (1948) doi: 10.1098/rspa.1948.0103
    [2] L. L. Foldy and S. A. Wouthuysen, Phys. Rev., 78: 29 (1950) doi: 10.1103/PhysRev.78.29
    [3] S. Tani, Prog. Theor. Phys., 6: 267 (1951) doi: 10.1143/ptp/6.3.267
    [4] L. L. Foldy, Phys. Rev., 87: 688 (1952) doi: 10.1103/PhysRev.87.688
    [5] K. M. Case, Phys. Rev., 95: 1323 (1954) doi: 10.1103/PhysRev.95.1323
    [6] J. Jayaraman, J. Phys. A: Math. Gen., 8: L1 (1975) doi: 10.1088/0305-4470/8/1/001
    [7] A. J. Silenko, Phys. Rev. A, 94: 032104 (2016) doi: 10.1103/PhysRevA.94.032104
    [8] W. Greiner, Relativistic Quantum Mechanics Wave Equations (Springer, 1990) pp. 277−290
    [9] A. A. Michelson and E. W. Morley, Philos. Mag. (Series 5), 24: 463 (1887) doi: 10.1080/14786448708628131
    [10] A. Sommerfeld, Sci. Nat., 28: 417 (1940) doi: 10.1007/BF01490583
    [11] A. Y. Silenko, Theor. Math. Phys., 105: 1224 (1995) doi: 10.1007/BF02067491
    [12] A. J. Silenko, J. Math. Phys., 44: 2952 (2003) doi: 10.1063/1.1579991
    [13] M. Hiramoto, T. Fujita, and T. Asaga, Prog. Theor. Phys., 106: 1223 (2001) doi: 10.1143/PTP.106.1223
    [14] K. Pachucki, Phys. Rev. A, 71: 012503 (2005) doi: 10.1103/PhysRevA.71.012503
    [15] K. Y. Bliokh, EPL, 72: 7 (2005) doi: 10.1209/epl/i2005-10205-1
    [16] P. Gosselin, A. Bérard, and H. Mohrbach, Eur. Phys. J. B, 58: 137 (2007) doi: 10.1140/epjb/e2007-00212-6
    [17] P. Gosselin, A. Bérard, H. Mohrbach, and S. Ghosh, Eur. Phys. J. C, 59: 883 (2009) doi: 10.1140/epjc/s10052-008-0839-4
    [18] M. Lippert, T. Bruckel, T. Kohler, and J. R. Schneider, EPL, 27: 537 (1994) doi: 10.1209/0295-5075/27/7/008
    [19] S. A. Khan, Optik, 117: 481 (2006) doi: 10.1016/j.ijleo.2005.11.010
    [20] Y. Hinschberger and P.-A. Hervieux, Phys. Lett. A, 376: 813 (2012) doi: 10.1016/j.physleta.2012.01.023
    [21] L. Fishman, J. Math. Phys., 33: 1887 (1992) doi: 10.1063/1.529666
    [22] M. Reiher, Theor. Chem. Acc., 116: 241 (2006) doi: 10.1007/s00214-005-0003-2
    [23] J. Autschbach, Coord. Chem. Rev., 251: 1796 (2007) doi: 10.1016/j.ccr.2007.02.012
    [24] W. Liu, Mol. Phys., 108: 1679 (2010) doi: 10.1080/00268971003781571
    [25] T. Nakajima and K. Hirao, Chem. Rev., 112: 385 (2012) doi: 10.1021/cr200040s
    [26] D. Peng and M. Reiher, J. Chem. Phys., 136: 244108 (2012) doi: 10.1063/1.4729788
    [27] D. Peng and M. Reiher, Theor. Chem. Acc., 131: 1081 (2012) doi: 10.1007/s00214-011-1081-y
    [28] H. Arodź and K. Golec, Phys. Rev. D, 33: 534 (1986) doi: 10.1103/PhysRevD.33.534
    [29] S. Scherer, G. Poulis, and H. Fearing, Nucl. Phys. A, 570: 686 (1994) doi: 10.1016/0375-9474(94)90079-5
    [30] V. P. Neznamov, Phys. Part. Nucl., 37: 86 (2006) doi: 10.1134/S1063779606010023
    [31] A. J. Silenko and O. V. Teryaev, Phys. Rev. D, 71: 064016 (2005) doi: 10.1103/PhysRevD.71.064016
    [32] P. Gosselin and H. Mohrbach, Eur. Phys. J. C, 71: 1739 (2011) doi: 10.1140/epjc/s10052-011-1739-6
    [33] Y. N. Obukhov, A. J. Silenko, and O. V. Teryaev, Phys. Rev. D, 88: 084014 (2013) doi: 10.1103/PhysRevD.88.084014
    [34] Y. N. Obukhov, A. J. Silenko, and O. V. Teryaev, Phys. Rev. D, 90: 124068 (2014) doi: 10.1103/PhysRevD.90.124068
    [35] A. J. Silenko, Nucl. Instrum. Methods Phys. Res. Sect. B, 114: 259 (1996) doi: 10.1016/0168-583X(96)00212-1
    [36] A. J. Silenko, Theor. Math. Phys., 112: 922 (1997) doi: 10.1007/BF02634108
    [37] G. Junker and A. Inomata, J. Math. Phys., 59: 052301 (2018) doi: 10.1063/1.5020545
    [38] B. Gonçalves, Y. N. Obukhov, and I. L. Shapiro, Phys. Rev. D, 75: 124023 (2007) doi: 10.1103/PhysRevD.75.124023
    [39] J. Q. Quach, Phys. Rev. D, 92: 084047 (2015) doi: 10.1103/PhysRevD.92.084047
    [40] J. Q. Quach, Phys. Rev. D, 93: 104048 (2016) doi: 10.1103/PhysRevD.93.104048
    [41] T.-W. Chen and D.-W. Chiou, Phys. Rev. A, 82: 012115 (2010) doi: 10.1103/PhysRevA.82.012115
    [42] T.-W. Chen and D.-W. Chiou, Phys. Rev. A, 90: 012112 (2014) doi: 10.1103/PhysRevA.90.012112
    [43] J. Meng, H. Toki, S. Zhou et al, Prog. Part. Nucl. Phys., 57: 470 (2006) doi: 10.1016/j.ppnp.2005.06.001
    [44] J. Meng and S. G. Zhou, J. Phys. G: Nucl. Part. Phys., 42: 093101 (2015) doi: 10.1088/0954-3899/42/9/093101
    [45] J. Meng, ed., Relativistic Density Functional for Nuclear Structure (World Scientific, 2016)
    [46] J.-Y. Guo, Phys. Rev. C, 85: 021302 (2012) doi: 10.1103/PhysRevC.85.021302
    [47] F. Wegner, Ann. Phys. (Berl.), 506: 77 (1994) doi: 10.1002/andp.19945060203
    [48] A. Bylev and H. Pirner, Phys. Lett. B, 428: 329 (1998) doi: 10.1016/S0370-2693(98)00417-1
    [49] H. Liang, J. Meng, and S.-G. Zhou, Phys. Rep., 570: 1 (2015) doi: 10.1016/j.physrep.2014.12.005
    [50] Y. Guo and H. Liang, Phys. Rev. C, 99: 054324 (2019) doi: 10.1103/PhysRevC.99.054324
    [51] Y. Chen, L. Li, H. Liang et al, Phys. Rev. C, 85: 067301 (2012) doi: 10.1103/PhysRevC.85.067301
    [52] W. Zhang and Y.-F. Niu, Chin. Phys. C, 41: 094102 (2017) doi: 10.1088/1674-1137/41/9/094102
    [53] Z. Xu and Z.-P. Li, Chin. Phys. C, 41: 124107 (2017) doi: 10.1088/1674-1137/41/12/124107
    [54] Y.-Y. Zong and B.-Y. Sun, Chin. Phys. C, 42: 024101 (2018) doi: 10.1088/1674-1137/42/2/024101
    [55] R. Swain, S. K. Patra, and B. B. Sahu, Chin. Phys. C, 42: 084102 (2018) doi: 10.1088/1674-1137/42/8/084102
    [56] Z.-Z. Li, S.-Y. Chang, Q. Zhao et al, Chin. Phys. C, 43: 074107 (2019) doi: 10.1088/1674-1137/43/7/074107
  • 加载中

Cited by

1. Shi, X.-X., Huang, B., Cao, X.-N. Exploring relativistic symmetry in shell evolution in the N=32 and N=34 isotones within a relativistic point-coupling model and the similarity renormalization group[J]. Physical Review C, 2025, 111(3): 034304. doi: 10.1103/PhysRevC.111.034304
2. Haouam, I.. Foldy-Wouthuysen transformation of Dirac equation in the context of conformable fractional derivative[J]. Modern Physics Letters A, 2024, 39(33-34): 2450155. doi: 10.1142/S0217732324501554
3. Li, Z., Chen, S., Zhou, M. et al. Covariant density functional theory for nuclear fission based on a two-center harmonic oscillator basis[J]. Physical Review C, 2024, 109(6): 064310. doi: 10.1103/PhysRevC.109.064310
4. Haouam, I.. Foldy-Wouthuysen transformation of Dirac equation in the context of Magueijo-Smolin model of the doubly special relativity[J]. International Journal of Modern Physics A, 2024, 39(07n08): 2450037. doi: 10.1142/S0217751X24500374
5. Zhao, Q., Ren, Z., Zhao, P. et al. Covariant density functional theory with localized exchange terms[J]. Physical Review C, 2022, 106(3): 034315. doi: 10.1103/PhysRevC.106.034315
6. Huang, B., Guo, J.-Y., Chen, S.-W. Investigation of pseudospin and spin symmetries in relativistic mean field theory combined with a similarity renormalization group approach[J]. Physical Review C, 2022, 105(5): 054313. doi: 10.1103/PhysRevC.105.054313
7. Wang, Z., Naito, T., Liang, H. et al. Exploring effects of tensor force and its strength via neutron drops[J]. Chinese Physics C, 2021, 45(6): 064103. doi: 10.1088/1674-1137/abf036
8. Ren, Z.X., Zhao, P.W. Toward a bridge between relativistic and nonrelativistic density functional theories for nuclei[J]. Physical Review C, 2020, 102(2): 021301. doi: 10.1103/PHYSREVC.102.021301
9. Xu, Q., Luo, L.-F. Investigation of relativistic symmetry from different fields with the similarity renormalization group[J]. European Physical Journal A, 2020, 56(5): 129. doi: 10.1140/epja/s10050-020-00159-z
10. Guo, Y., Liang, H. Nonrelativistic expansion of the Dirac equation with spherical scalar and vector potentials by a reconstituted Foldy-Wouthuysen transformation[J]. Physical Review C, 2020, 101(2): 024304. doi: 10.1103/PhysRevC.101.024304

Figures(1)

Get Citation
Yixin Guo and Haozhao Liang. Non-relativistic expansion of single-nucleon Dirac equation: Comparison between Foldy-Wouthuysen transformation andsimilarity renormalization group[J]. Chinese Physics C. doi: 10.1088/1674-1137/43/11/114105
Yixin Guo and Haozhao Liang. Non-relativistic expansion of single-nucleon Dirac equation: Comparison between Foldy-Wouthuysen transformation andsimilarity renormalization group[J]. Chinese Physics C.  doi: 10.1088/1674-1137/43/11/114105 shu
Milestone
Received: 2019-06-23
Article Metric

Article Views(3981)
PDF Downloads(66)
Cited by(10)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Non-relativistic expansion of single-nucleon Dirac equation: Comparison between Foldy-Wouthuysen transformation andsimilarity renormalization group

    Corresponding author: Haozhao Liang, haozhao.liang@riken.jp
  • 1. Department of Physics, Graduate School of Science, The University of Tokyo, Tokyo 113-0033, Japan
  • 2. RIKEN Nishina Center, Wako 351-0198, Japan
  • 3. Department of Modern Physics, University of Science and Technology of China, Hefei 230026, China

Abstract: By following the Foldy-Wouthuysen (FW) transformation of the Dirac equation, we derive the exact analytic expression up to the 1/M4 order for general cases in the covariant density functional theory. The results are compared with the corresponding ones derived from another novel non-relativistic expansion method, the similarity renormalization group (SRG). Based on this comparison, the origin of the difference between the results obtained with the FW transformation and the SRG method is explored.

    HTML

    1.   Introduction
    • The Foldy-Wouthuysen (FW) transformation, also called the Pryce-Tani-Foldy-Wouthuysen transformation [14], was initially formulated around 1950 to describe the non-relativistic limit of relativistic electrons. A great advantage of the FW transformation is the simple form of operators corresponding to classical observables. This transformation for spin-1/2 particles was soon extended to the cases of spin-0 and spin-1 particles [5], and later generalized to the case of arbitrary spins [6]. Recently, the exact form of the exponential FW operator for a particle with an arbitrary spin was given in Ref. [7].

      The FW transformation is also an elegant method for non-relativistic expansion of the Dirac equation [8]. One milestone of its applications in atomic physics was that, using the FW transformation in the presence of external fields, the fine structure [9, 10], describing the splitting of spectral lines of atoms due to the electron spin and the relativistic corrections to the non-relativistic Schrödinger equation, was successfully testified [8].

      Due to its elegant and powerful scheme, the FW transformation has been applied in very diverse areas, such as electrodynamics [11, 12], atomic physics [13, 14], condensed matter physics [1517], optics [1820], acoustics [21], quantum chemistry [2227], nuclear physics [28, 29], quantum field theory [30], gravity [3134], in the theory of the weak interaction [35, 36], and even in the scheme of supersymmetric quantum mechanics [37]. In recent years, the FW transformation has also been discussed in a gravitational wave background of the generalized Dirac Hamiltonian in the non-relativistic limit [3840].

      Nevertheless, in most of the above studies and applications of the FW transformation, only the leading-order terms were derived and treated, such as the relativistic correction to the kinetic energy appearing up to the order of 1/M3 (with M the bare mass of the particle), and the correction due to the spin-orbit coupling and the Darwin term appearing in the order of 1/M2 [8]. Recently, the FW non-relativistic expansions have been carried out up to the order of 1/M7 [41], and even 1/M14 [42]. Nevertheless, such derivations are valid only for a family of potentials with certain properties.

      For systems of atomic nuclei, the covariant density functional theory [4345] is one of the state-of-the-art methodologies. During the past decades, it has achieved great success in describing various nuclear phenomena in both stable and exotic nuclei. In this theoretical framework, the equation of motion for nucleons is described by the Dirac equation which includes not only the strong vector potential V(r) , but also the strong scalar potential S(r), whose typical values are V(r)300 MeV and S(r)350 MeV, while the mass of nucleon is M=939 MeV. As a result, the non-relativistic expansion of the single-nucleon Dirac equation is non-trivial, and its high-order terms have to be derived and verified carefully.

      In 2012, Guo [46] applied the similarity renormalization group (SRG) [47, 48], instead of the FW transformation, for the non-relativistic expansion of the single-nucleon Dirac equation up to the 1/M3 order for investigating the nuclear pseudospin symmetry [49]. Using the SRG method, the Dirac Hamiltonian is gradually transformed into a diagonal form with the evolution controlled by the flow equation. As a result, the eigenequations for the upper and lower components of the Dirac spinors are decoupled at the end of the flow. Also, the non-relativistic reduced Hamiltonian thus obtained is Hermitian, and can also be expanded into the series in terms of 1/Mi. To achieve the accuracy of single-particle energies at the 0.1 MeV level, we demonstrated recently that the 1/M4-order terms are needed, and presented a detailed derivation and results in Ref. [50].

      Both the FW transformation and the SRG method provide a systematic way to carry out a non-relativistic expansion of the Dirac equation up to an arbitrary order. However, the FW transformation employs finite steps of unitary transformation, while the SRG method employs infinite infinitesimal unitary transformations along the flow. Therefore, it is interesting and important to investigate the similarities and differences between these two approaches.

      In this paper, we perform an exact non-relativistic expansion of the single-nucleon Dirac equation up to the 1/M4 order by using the FW transformation. We also investigate the difference between such results and those obtained by the SRG method.

      The paper is organized as follows. In Section 2.1, the SRG method is first recalled and the expansion up to the 1/M4 order is shown. The new expansion up to the same order using the FW transformation in a general case, such as the inclusion of the scalar potential, is derived in Section 2.2. The specific form of the result in spherical symmetry is presented in Section 3.1, and the difference between the SRG and FW results, as well as the origin of the difference, is investigated in Section 3.2. Finally, a summary and perspectives are given in Section 4.

    2.   Theoretical framework

      2.1.   Non-relativistic expansion by SRG

    • In the theoretical framework of covariant density functional theory [4345, 5156], the Dirac Hamiltonian for nucleons reads

      H=αp+β(M+S)+V,

      (1)

      where α and β are the Dirac matrices, M is the mass of the nucleon, and S and V are the scalar and vector potentials, respectively. In this section, we recall the technique of SRG [47, 48] and show the results up to the 1/M4 order [46, 50].

      The basic idea of SRG is that the Dirac Hamiltonian in Eq. (1) is transformed by a unitary operator U(l) as

      H(l)=U(l)HU(l),H(0)=H,

      (2)

      with a flow parameter l. While U(l) is not known explicitly, an anti-Hermitian generator η(l) related to U(l) as η(l)=dU(l)dlU(l) is chosen, and the Hamiltonian H(l) evolves with the differential flow equation

      dH(l)dl=[η(l),H(l)].

      (3)

      One of the convenient and appropriate choices of η(l) reads [48]

      η(l)=[βM,H(l)],

      (4)

      which transforms the Dirac Hamiltonian into a diagonal form.

      For solving Eq. (3), the Hamiltonian is first divided into two parts:

      H(l)=ε(l)+o(l),

      (5)

      according to the commutation and anti-commutation relations with respect to β, i.e., [ε,β]=0 and {o,β}=0. As a result,

      dε(l)dl=4Mβo2(l),

      (6a)

      do(l)dl=2Mβ[o(l),ε(l)],

      (6b)

      with the initial conditions

      ε(0)=β(M+S)+V,o(0)=αp.

      (7)

      Equations (6a) and (6b) can then be solved with the expansion into the series of 1/Mi [48], which gives

      ε(λ)M=i=0εi(λ)Mi,

      (8a)

      o(λ)M=j=1oj(λ)Mj,

      (8b)

      while the flow parameter is replaced by a dimensionless parameter λ=lM2. The solutions are obtained as [46]

      εn(λ)=εn(0)+4βλ0n1k=1ok(λ)onk(λ)dλ,

      (9a)

      on(λ)=on(0)e4λ+2βe4λλ0n1k=1[e4λok(λ),εnk(λ)]dλ,

      (9b)

      with the initial conditions,

      ε0(0)=β,ε1(0)=βS+V,εn(0)=0ifn2,o1(0)=αp,on(0)=0ifn2.

      (10)

      Therefore, at the end of the flow λ, all the off-diagonal parts vanish, i.e., on()=0 for all n. The diagonalized Dirac operator up to the 1/M4 order is not difficult to obtain as [50]

      HSRG=ε()=Mε0()+ε1()+ε2()M+ε3()M2+ε4()M3+ε5()M4+=Mε0(0)+ε1(0)+12Mβo21(0)+18M2[[o1(0),ε1(0)],o1(0)]+132M3β(4o41(0)+[[o1(0),ε1(0)],ε1(0)]o1(0)+o1(0)[[o1(0),ε1(0)],ε1(0)]2[o1(0),ε1(0)][o1(0),ε1(0)])+1128M4(9[[o1(0),ε1(0)],o31(0)]+3[o1(0),ε1(0)]2ε1(0)+3ε1(0)[o1(0),ε1(0)]26[o1(0),ε1(0)]ε1(0)[o1(0),ε1(0)]+3[o1(0)[o1(0),ε1(0)]o1(0),o1(0)+[[[[o1(0),ε1(0)],ε1(0)],ε1(0)],o1(0)])+

      (11)
    • 2.2.   Non-relativistic expansion with the FW transformation

    • In this section, we perform the exact non-relativistic expansion of the single-nucleon Dirac equation up to the 1/M4 order by using the FW transformation for a general case, in particular for the scalar potential.

      According to the FW transformation in the presence of external fields [2, 4, 8], the corresponding operators are defined as

      O=αp,ε=βS+V,Λ=i2MβO,

      (12)

      where the operators O and ε satisfy Oβ=βO and εβ=βε, respectively. The Hamiltonian (1) reads

      H=βM+O+ε.

      (13)

      Based on the above definitions, the Dirac Hamiltonian is transformed by a unitary operation into [8]

      H=eiΛHeiΛ=H+i[Λ,H]+i22![Λ,[Λ,H]]++inn![Λ,[Λ,,[Λn,H]]]+

      (14)

      where

      inn![Λ,[Λ,,[Λn,H]]]=(1)n(n1)2βnn!Mn(OnβM+On+1+12n[O,[O,,[On,ε]]]).

      (15)

      Keeping all the terms up to the 1/Mn order, the unitary transformed Hamiltonian reads

      H1/Mn=βM+ε+n+1k=1(1)(k1)(k2)2k1k!Mk1βk1Ok+nk=0(1)k(k1)2βk2kk!Mk[O,[O,,[Ok,ε]]].

      (16)

      For example, the corresponding result up to the 1/M4 order is

      H1/M4=βM+ε+12MβO213M2O318M3βO4+130M4O5+12Mβ[O,ε]18M2[O,[O,ε]]148M3β[O,[O,[O,ε]]]+1384M4[O,[O,[O,[O,ε]]]].

      (17)

      In the unitary transformed Hamiltonian, e.g., Eq. (17), the off-diagonal parts are not zero but raised by one order from O to 12Mβ[O,ε] (plus higher-order terms). In order to make the off-diagonal parts vanish, more precisely speaking to make them higher than a given order, one should repeat the FW transformation until the accuracy is achieved [8]. For that purpose, the operators O and ε can be redefined according to Eq. (17), and one has

      ε=ε+12MβO218M2[O,[O,ε]]18M3βO4+1384M4[O,[O,[O,[O,ε]]]],

      (18a)

      O=12Mβ[O,ε]13M2O3148M3β[O,[O,[O,ε]]]+130M4O5,

      (18b)

      Λ=i2MβO.

      (18c)

      The corresponding FW transformation reads

      H=eiΛHeiΛ.

      (19)

      It can be seen that the off-diagonal parts in H are of the order of 1/M2. Therefore, one should repeat this procedure to H′′′′′ to make the off-diagonal parts of the order of 1/M5. Of course, to get the resultant diagonal parts in H′′′′′ up to the 1/M4 order, a recursion technique makes the calculation less complicated than it looks.

      As a result, the non-relativistic expansion with the FW transformation up to the 1/M4 order reads

      HFW=βM+ε+12MβO218M2[O,[O,ε]]18M3βO418M3β[O,ε][O,ε]+1384M4[O,[O,[O,[O,ε]]]]+112M4[O3,[O,ε]]+132M4[O,ε][O,ε]ε116M4[O,ε]ε[O,ε]+132M4ε[O,ε][O,ε].

      (20)
    3.   Results and discussion

      3.1.   Results for spherical symmetry

    • For the systems with spherical symmetry, i.e., for spherical nuclei, the corresponding radial single-nucleon Dirac equation reads [4345]

      (Σ(r)+Mddr+κrddr+κrΔ(r)M)(G(r)F(r))=E(G(r)F(r)),

      (21)

      where κ is a good quantum number defined as κ=(j+1/2) for j=l±1/2, and Σ(r)=V(r)+S(r) and Δ(r)=V(r)S(r) are the sum and the difference of the vector and scalar potentials, respectively. The single-particle energy E=ε+M includes the mass of the nucleon. The operators ε and O read

      ε=(Σ(r)00Δ(r)),O=(0ddr+κrddr+κr0).

      (22)

      According to Eq. (20), the Dirac Hamiltonian is transformed by the FW transformation as

      (H(F)FW+MO(1M5)O(1M5)H(D)FWM),

      (23)

      where the superscripts (F) and (D) represent the components of the Dirac Hamiltonian acting on single-particle states in the Fermi sea and Dirac sea, respectively. It is seen that the off-diagonal parts are not strictly zero, which is different from the results of SRG, but are of higher order than the required one. Focusing on single-particle states in the Fermi sea which correspond to their counterparts in the non-relativistic framework, the explicit expansion of H(F)FW up to the 1/M4 order is worked out in detail as

      H(F)0,FW=Σ(r),

      (24a)

      H(F)1,FW=12Mp2,

      (24b)

      H(F)2,FW=18M2(4Sp2+4Sddr2κrΔ+Σ),

      (24c)

      H(F)3,FW=18M3(p4+4S2p28SSddr2SΣ+4SΔκr+ΣΔ),

      (24d)

      H(F)4,FW=1384M4{144Sp4288Sp2ddr+[72Δκr24(4Σ+9S)192S3]p2+[72Δκr272Δκr+24(3S+4Σ)+576S2S]ddr+[12(5Σ36S)κ(κ+1)r3+12(5Σ+12S)κ(κ+1)r272Δκr3+72Δκr236Δκr288S2Δκr+33Σ+144S2Σ+24SΣ(2Σ6Δ)]},

      (24e)

      where

      p2=d2dr2+κ(κ+1)r2,

      (25)

      and

      p4=d4dr42κ(κ+1)r2d2dr2+4κ(κ+1)r3ddr+κ(κ+1)(κ+3)(κ2)r4.

      (26)

      In contrast, the Dirac Hamiltonian transformed by the SRG method reads

      (H(F)SRG+M00H(D)SRGM).

      (27)

      Its off-diagonal parts strictly vanish up to infinite order. According to Eqs. (11), the explicit expansion of H(F)SRG up to the 1/M4 order reads [50]

      H(F)0,SRG=Σ(r),

      (28a)

      H(F)1,SRG=12Mp2,

      (28b)

      H(F)2,SRG=18M2(4Sp2+4Sddr2κrΔ+Σ),

      (28c)

      H(F)3,SRG=132M3(4p4+16S2p232SSddr8SΣ+16SΔκr2Σ2+4ΣΔ),

      (28d)

      H(F)4,SRG=1128M4{48Sp496Sp2ddr+[24Δκr24(Σ+3S)64S3]p2+[24Δκr224Δκr+24(S+Σ)+192S2S]ddr

      +[12(Σ12S)κ(κ+1)r3+12(Σ+4S)κ(κ+1)r224Δκr3+24Δκr212Δκr96S2Δκr+9Σ+48S2Σ+24SΣ(Σ2Δ)]}.

      (28e)
    • 3.2.   Comparison between FW and SRG approaches

    • By comparing Eqs. (11) and (20), it is found that

      H0,SRGH0,FW=0,

      (29a)

      H1,SRGH1,FW=0,

      (29b)

      H2,SRGH2,FW=0,

      (29c)

      H3,SRGH3,FW=β32M3[[O2,ε],ε],

      (29d)

      H4,SRGH4,FW=164M4[[O2,ε],O2]+1128M4[([ε2,O2]2[ε,OεO]),ε].

      (29e)

      This result seems to indicate that the FW transformation and the SRG method agree with each other up to the 1/M2 order, but they lead to different results starting from the 1/M3 order. However, considering the infinite mass limit M, all expressions should be strictly organized order by order, and thus the difference between Eqs. (27) and (23) in the 1/M5 order cannot lead to the differences of the results in the 1/M3 and 1/M4 orders. This puzzle needs to be further investigated.

      It turns out that the difference shown in Eq. (29) comes from an additional unitary transformation, after the Dirac Hamiltonian is decoupled into the upper and lower parts. Let

      Ξ=iβ32M3[O2,ε]i128M4([ε2,O2]2[ε,OεO]).

      (30)

      This is a Hermitian and diagonal operator, i.e., Ξ=Ξ and βΞ=Ξβ. Acting an additional unitary transformation on HFW, it reads

      eiΞHFWeiΞ=HFW+i[Ξ,HFW]+

      (31)

      Keeping all terms up to the 1/M4 order, the result is

      eiΞHFWeiΞ=HFW+β32M3[[O2,ε],ε]+164M4[[O2,ε],O2]+1128M4[([ε2,O2]2[ε,OεO]),ε],

      (32)

      which is nothing else but HSRG.

      In the spherical case, the explicit expressions for operators β32M3[[O2,ε],ε], 164M4[[O2,ε],O2], and 1128M4[([ε2,O2]2[ε,OεO]),ε] , acting on single-particle states in the Fermi sea, respectively read

      116M3Σ2,

      (33a)

      164M4[4Σp24Σddr+4Σκ(κ+1)r34Σκ(κ+1)r2Σ],

      (33b)

      116M4SΣ2.

      (33c)

      These equations explain the difference between Eqs. (24d) and (28d), as well as between Eqs. (24h) and (28h).

      Since eiΞ is an additional unitary transformation acting on the already decoupled Hamiltonian, it does not affect the non-relativistic expansion of the Dirac equation, and the single-particle spectra obtained by the FW transformation and the SRG method are the same.

      The unitary transformation eiΞ is equivalent to a rotation in a linear and closed complex space, and in quantum mechanics, the operation or rotation in such a space should be associated with a corresponding generator. However, in the present context, it is nontrivial to write down a general expression for such a generator, or derive a conserved quantity that is asymptotically related to a generator up to an arbitrary order. This is because the number of unitary transformations from the FW results to the SRG ones would in principle be infinite, since the SRG method employs infinite infinitesimal unitary transformations along the flow. Therefore, the physics behind the introduced unitary transformation is open for future studies.

    4.   Summary and perspectives
    • In this work, we investigated the non-relativistic expansion of the single-nucleon Dirac equation in the theoretical framework of covariant density functional theory for a general case where the scalar potential is included. We worked out the exact analytical expansion up to the 1/M4 order using the FW transformation. A further investigation of the difference between Eqs. (27) and (23) explained the puzzle that a disagreement seems to appear between the results obtained with the FW transformation and the SRG method, i.e. the non-relativistic expansion of the Dirac equation is affected, but the single-particle spectra obtained by the FW transformation and the SRG method are the same.

      The non-relativistic expansion of the Dirac equation for nuclear systems using the FW transformation has been extended to as high orders as needed. Similarly to the SRG method, we anticipate that this novel non-relativistic expansion method will establish in future studies a potential bridge between the relativistic and non-relativistic density functional theories. In particular, with the present FW transformation, all unitary transformations involved are in an explicit form. This leads in a straightforward way to the corresponding transformations of the other relevant operators, such as the one-body density operators.

Reference (56)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return