Deuteron under strong magnetic fields: Skyrme model study

Figures(4)

Get Citation
Yan-jun Chen and Jiang Wei. Deuteron under strong magnetic flelds: Skyrme model study[J]. Chinese Physics C. doi: 10.1088/1674-1137/aca8f2
Yan-jun Chen and Jiang Wei. Deuteron under strong magnetic flelds: Skyrme model study[J]. Chinese Physics C.  doi: 10.1088/1674-1137/aca8f2 shu
Milestone
Received: 2022-09-27
Article Metric

Article Views(1314)
PDF Downloads(32)
Cited by(0)
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:

Deuteron under strong magnetic fields: Skyrme model study

  • 1. Department of Physics and Electronic Science, Changsha University of Science and Technology, Changsha 410114, China
  • 2. Hunan Provincial Key Laboratory of Flexible Electronic Materials Genome Engineering, Changsha 410114, China

Abstract: We explore the deuteron under strong magnetic fields in Skyrme models. The effects of the derivative dependent sextic term in the Skyrme Lagrangian are investigated, and the rational map approximation is used to describe the deuteron. The influences of strong magnetic fields on the electric charge distribution and mass of the deuteron are discussed.

    HTML

    I.   INTRODUCTION
    • Recent investigations have shown that a strong magnetic field may exist in at least two high-energy astrophysical objects. First, the magnetic fields on the surface of magnetars are expected to reach 1015G, whereas in the inner core of magnetars, the magnitude of the magnetic fields could be of the order of 1018 G [1]. Second, it is considered that during the electroweak phase transition in the early universe, extremely strong magnetic fields exceeding 1020 G might have been produced [2]. In the course of relativistic non-central heavy-ion collisions, the strength of magnetic fields can reach up to 1018 G at the Relativistic Heavy Ion Collider (RHIC) and 1020 G at the Large Hadron Collider (LHC) [3, 4]. Although such strong magnetic fields in these heavy-ion colliding experiments decay during a very short period, they can still induce significant effects on the phase structure of QCD matter. Many theoretical investigations based on the lattice QCD method or effective QCD theories, such as NJL models, have shown that strong magnetic fields can modify the quark potential and meson properties owing to complex QED-QCD interactions [57]. However, there are relatively few studies in regard to the modification of baryon properties in the presence of strong magnetic fields [810].

      It is difficult to study low-energy strong interaction physics directly from QCD theory owing to its non-perturbative character at the low-energy scale. Alternative approaches, such as lattice QCD or effective field theories (EFTs), must instead be considered. The Skyrme model, a promising EFT, was first proposed by Skyrme in 1961 [11]. In this model, meson fields are the basic degrees of freedom, and baryons emerge as collective excitations or topological solitons, known as skyrmions, from nonlinear interactions between mesons, the topological charge of which is identified by the baryon number. The Lagrangian of the standard Skyrme model consists of three terms. The first is the sigma-model term, quadratic in its first derivatives, which provides the kinetic energy of π mesons, the second is the Skyrme term, which is quartic in its first derivatives, and relates to the stable solutions of models, and the third is the potential term, which gives explicit breaking of chiral symmetry and relates to the π meson mass. Although this standard Skyrme model can successfully describe many nuclear phenomena, there are still several defects. For example, it provides excessively high binding energies of atomic nuclei and thus is generally difficult to apply to the investigation of matter with a large baryon number. Many efforts have been made to improve the standard Skyrme model (see, for example, Refs. [1218] and the references therein). In recent years, so-called generalized Skyrme models have been proposed [19, 20], where a term sextic in derivatives is added and may be related to vector meson exchange based on phenomenological arguments, which can overcome the shortcomings of the original Skyrme models.

      In this study, we modify the method proposed by Ref. [8] to study the deuteron as a B=2 skyrmion. Besides the conventional Skyrme model (labeled as the SKM), the effects of the derivative dependent sextic Lagrangian term are also considered (labeled as the gSKM). This article is organized as follows: In Sec. II, we describe the formulas necessary. Here, we quantize the skyrmion within the semi-classical quantization method [21, 22]. Note that the quantization method in Ref. [8] consists of allowing the SU(2)valued field in the Skyrme model to rotate in three-dimensional isospin space and spatial space only in the xy plane, which is inconsistent with the constraint of Eq. (21) resulting from the cylindrical symmetry implicitly imposed on the solution of the skyrmion in this study when we quantize the B=2 skyrmion. Therefore, in the present study, we modify the quantization method in Ref. [8] and consider three-dimensional rotation in spatial space as well as in isospin space. In Sec. III, the calculated results and discussions are given. Finally, a summary is presented in Sec. IV.

    II.   FORMALISM
    • In this study, the action of the Skyrme model coupled with magnetic fields is ΓWZW(U)+d4x£0. Here, the Wess-Zumino-Witten (WZW) term represents chiral anomaly effects [23, 24], and the gauged WZW action ΓWZW(U)=Γc(U)Γc(1) in the presence of electromagnetic fields is given in Ref. [25], with

      Γc(U)=i80π2d5xεijklmTr[RiRjRkRlRm]+116π2d4xεμναβTr{eAμQE(RνRαRβ+LνLαLβ)+ie22[AμQERνAαQERβAμQELνAαQELβ]ie2[AμQEUAνQEURαRβAμQEUAνQEULαLβμAνQEUAαQEURβμAνQEUAαQEULβ(AμQEνAαQE+μAνQEAαQE)(Rβ+Lβ)]+e3AμQEAνQEAαQE(Rβ+Lβ)+e3(AμQEUAνQEUAαQERβ+AμQEUAνQEUAαQELβ+AμQEUAνQEαAβQEU+AμQEUνAαQEAβQEU+μAνQEAαQEUAβQEU+AμQEνAαQEUAβQEU)ie4(AμQEAνQEAαQEUAβQEU+AμQEUAνQEAαQEAβQEU+12AμQEUAνQEUAαQEUAβQEU)},

      (1)

      where the first term is a five-dimensional integral whose boundary is four-dimensional space-time, Ri=UiU, Li=iUU, e is the unit electric charge, and QE=I/6+τ3/2 is the electric charge matrix, where I is the rank 2 unit matrix, and τ3 is the third Pauli matrix. The generalized Skyrme Lagrangian £0 with local gauge invariance reads as

      £0=F2π16Tr(DμUDμU)+m2πF2π8Tr(U1)+132g2Tr([UDμU,UDνU]2)9π25m2BμBμ,

      (2)

      with the covariant derivative DμU defined as

      DμU=μUieAμ[QE,U].

      (3)

      Here, the Lagrangian has three free parameters (Fπ, g, m), which are used to fit several experimental values. The pion mass is set to its physical value mπ=138 MeV. In this study, the homogeneous magnetic field B along the z direction is considered, and the symmetric gauge of the electromagnetic field Aμ is chosen as Aμ=(0,yB/2, xB/2, 0). As an SU(2) valued field, U is expressed as

      U=exp(iF(r)τn(θ,φ)),

      (4)

      where F(r) is a real profile function satisfying the boundary conditions F(0)=π and F()=0, |n|=1. For baryon number B>1, we utilize the so-called rational map ansatz [26], and n is expressed as

      (n1,n2,n3)=11+|R|2(2Re(R),2Im(R),1|R|2),

      (5)

      where R(z) is the rational map between two Riemann spheres, with z as the Riemann sphere coordinate defined as z=tan(θ/2)exp(iφ). The form of R is p/q, where p and q are polynomials in z with the larger of their degrees being B. For the B=2 skyrmion, R=z2 is chosen. The last term in Eq. (2) is sextic in derivatives, in which Bμ reads as

      Bμ=124π2εμναβTr[(UDνU)(UDαU)(UDβU)],

      (6)

      which degenerates into the topological current in the absence of magnetic fields.

      Considering the axial symmetry imposed by magnetic fields, we follow the ansatz in Ref. [8] to deform the spatial space as

      x=cρrsin(θ)cos(φ),y=cρrsin(θ)sin(φ),z=czrcos(θ),

      (7)

      where cρ and cz are positive dimensionless parameters, and θ[0,π] and φ[0,2π] are polar angles. Then, the volume element is dV=dxdydz=c2ρczr2sin(θ)drdθdφ. The parameters cρ and cz have two effects: they deform the shape and scale the volume. Because the scale effect of cρ and cz can be absorbed by performing the scale transform of r, we follow the choice in Ref. [8] in which cρ1/cz to seek results with a conserved volume. It should be mentioned that in Ref. [8], which studies B=1 baryons, cz characterizes the deviation degree from a sphere. However, in the present study, the deuteron is toroidal and based on the rational map ansatz, which has already deviated from a sphere. Thus, cz in this study can be understood as a new parameter introduced into the conventional rational map approximation (cRMA), and cz1 means a deviation from the cRMA.

      In the semi-classical quantization method, to obtain the spin and isospin of baryons, U is given time dependence by rotating U in both space and isospin as

      ˆU=A(U(˜S))A,

      (8)

      in which A and ˜S are the rotation matrices for the isospin and spatial rotations, respectively, satisfying

      A1˙A=i2ωaτa,(˜S1˙˜S)ij=εijkΩk,

      (9)

      where a,i,j,k=1,2,3. By replacing U with ˆU in £0, the static contribution to the deuteron mass can be expressed as

      Md,st=ψdM0+M1(eB)+M2(eB)2ψd,

      (10)

      in which

      M0=(πβg2)dr12c2zr2{3cz[16+c3z(8+3π)]+48[4+c3z(2+π)π]r2+96c2zm2πr496c2zm2πr4cos(F)4{cz[16+c3z(8+3π)]+12[4+c3z(2+π)π]r2}cos(2F)+cz[16+c3z(8+3π)]cos(4F)+4r2{6cz[2+c3z(4+π)π]+4(1+2c3z)r2+6cz[2+c3z(4+π)π]cos(2F)}F2+4(8+3π)κc2zsin4(F)F2},

      (11)

      M1=(R3S3πβg2)2dr3c2z{sin2(F){6c2z(4+π)r22c3z(289π)sin2(F)2(32+9π)sin2(F)[289π+4c3z(10+3π)]r2F2}+4κczsin4(F)F2},

      (12)

      and

      M2=(πβg2)r2dr480c3z{sin2(F){40c2z{2(1+S23)[83π+(28+9π)R23]+c3z(1+S23)[329π+(88+27π)R23]}r2+20sin2(F){8S23(28+9π+124R2339πR23)+4c6z(1+S23)[289π+2(47+15π)R23]+c3z[12839π+(448+141π)R23+S23[8027π+(400+129π)R23]]}4{8S23[48142R23+15π(1+3R23)]+4c6z(1+S23)[48142R23+15π(1+3R23)]+c3z[11245π+3(136+45π)R23+S23[4321288R23+135π(1+3R23)]]}r2F2}+10κczsin4(F)F2{2S23[8+3π+3(247π)R23]+c3z(1S23)[163π+3(24+7π)R23]}}.

      (13)

      Here, β=22/(gFπ) is used as the unit of length, κ=9g2/(5m2β2π2), R3 is defined as R3=Tr(τ3Aτ3A)/2, and S3 is the ˜S33 component of the rotation matrix ˜S.

      The replacement of U with ˆU can also give the kinematical part in the form

      T=12aiω2i+12biΩ2iciωiΩi+D(eB)S3(ω32Ω3)+E(eB)(S1Ω1+S2Ω2),

      (14)

      where i=1,2,3, S1 and S2 are the ˜S13 and ˜S23 components of the rotation matrix ˜S, respectively,

      D=drr22czπ(4+π)sin2(F)F,

      (15)

      E=drr22czπ[4π+2c3z(3+π)]sin2(F)F,

      (16)

      and a1=a2, b1=b2, c1=c2=0, and a3=b3/4=c3/2, where

      a1=(2πβg2)[a10+a11(eB)+a12(eB)2],

      (17)

      a3=(2πβg2)[a30+a32(eB)2],

      (18)

      b1=(2πβg2)[b10+b11(eB)+b12(eB)2].

      (19)

      The equations needed in Eqs. (17)–(19) are listed in Appendix A. We can then obtain the body-fixed canonical conjugate momenta of the isospin and spin by taking functional derivatives with respect to ωa and Ωa, respectively, as

      Ka=Tωa,La=TΩa.

      (20)

      Note that in Ref. [8], only the spatial rotation in the xy plane is considered, that is, only the third component of the spatial angular momentum is defined and k3 in Eq. (9), which results in the equivalence between the third component of the body-fixed (L) and space-fixed (J, physically more meaningful) spatial angular momenta, that is, J3=L3. When the constraint [The distinction between the baryon number B here and the magnetic field eB should always be clear by context.]

      (BK3+L3)|ψ=0,

      (21)

      which results from the cylindrical symmetry implicitly imposed on the solution of the skyrmion in this study, is cast into the calculation for the B=1 nucleon, as in Ref. [8], there is no trouble. This is because the quantum numbers (j3 and k3) of J3 and K3 can simultaneously be half-integers for the B=1 nucleon. However, the integer j3=1 for the B=2 deuteron allows k3 to be a half-integer, according to Eq. (21), which is prohibited in the quantization of the B=2 deuteron. Therefore, in this study, we modify the quantization method in Ref. [8] and consider three-dimensional rotation in the spatial space as well as in the isospin space. As a result, J3 and L3 become mutually independent. Consequently, the quantum numbers (l3 and k3) of L3 and K3 are allowed to be zero for the B=2 deuteron, and the constraint of Eq. (21) is not violated. Here, the wave function of the deuteron is given by the direct product |ii3k3|jj3l3, where i and i3 are the quantum numbers of the space-fixed isospin momentum (I2 and I3, physically more meaningful).

      Now, using the relations J3=SiLi and S2i1, where i=1,2,3, we can write the kinematical contribution to the deuteron mass as

      Md,kin=ψda=1,2,3(ωaKa+ΩaLa)Tψd,=ψdI22a1+J2+(1S23)E2(eB)22b1+[(K3S3D(eB))22a3(1a1+4b1)K232]+(J32S3K3)E(eB)b1ψd.

      (22)

      The deuteron mass is the sum of the static and kinematical contributions, that is, Md=Md,st+Md,kin. Md,st and Md,kin can be obtained with the help of the wave function, as stated in Ref. [20]. F(r) in Eq. (4) can be obtained at the O(NC) order by minimizing Md,st, with the help of a certain value of cz obtained by minimizing Md.

      The magnetic moment of the deuteron is defined by μMd/(eB). Because k3 is equal to zero in the wave function of the deuteron, and the expectation value of R3 is zero [20], most of the linear terms of magnetic fields in Eqs. (10) and (22) offer no contribution to the mass, except the last term in Eq. (22). Thus the magnetic moment is given by μ=β2j3E/b1 for |eB|=0 (the presence of β2 results from the consideration of the unit of length in this study). This equation for the magnetic moment is equivalent to that deduced in Ref. [27]. We take the model parameters as g= 4.84 and Fπ= 108 MeV for the SKM, and g= 19.48 and Fπ= 129.11 MeV for the gSKM with an additional parameter m chosen as m= 420 MeV [19]. It is known that the wave function of the deuteron, as the ground state 3S1 of the B=2 skyrmion, is |000|110. We then obtain the results of the magnetic moment for |eB|=0, which in the unit of the Bohr magneton μB=1/(2MN) for the SKM and gSKM are 0.69 and 0.77, which are 19% and 10% smaller, respectively, than the experimental value of 0.857.

    III.   RESULTS AND DISCUSSION
    • In Fig. 1, we illustrate the |eB| dependence of cz for the deuteron. We also give the calculations for the first excited state 1S0, the wave function of which is |100|000, for comparison. In the absence of magnetic fields, the results slightly deviate along the direction of the magnetic field from those in cRMA, with cz=1.03 in the SKM and cz=1.05 in the gSKM for the deuteron, which are smaller than those for the first excited state. When the magnetic fields increase, cz decreases continuously. The gSKM exhibits a larger descent velocity of cz with |eB| than SKM. Overall, we see that the deviation of the calculated results from cz=1 in the presence of the magnetic fields is 10 %.

      Figure 1.  (color online) |eB| dependence of cz for the deuteron in the ground (color black) and excited states (color red) for the gSKM (solid curves) and SKM (dashed curves).

      In the deuteron under the magnetic fields, the radial density of baryon charge is given by

      ρ(r)=ρ0(r)+ρB(r),=4sin2(F)Fπ+ψdr(eB)24πcz(4π)R3S3×[sin(2F)+rcos(2F)F]ψd,

      (23)

      where the integral of the first term ρ0(r) over r is 2, which is consistent with the baryon charge B=2 of the deuteron, and the second term ρB(r) is induced by magnetic fields. For B=1, the presence of ρB(r) under magnetic fields does not break baryon number conservation because integrating it over the entire space is equal to zero [8]. The fact that k3=0 and the expectation value of R3 in Eq. (23) is zero [20] indicates that ρB0 in the deuteron. Because i3=0 for the deuteron, we can evaluate the root mean square (RMS) electric charge radii, defined as R=(X2ρ0/2dV)1/2, in which X=x2+y2+z2, and Rx=(x2ρ0/2dV)1/2 and Rz=(z2ρ0/2dV)1/2 for directions of x and z, respectively. The results are illustrated in Fig. 2. The phenomenon that Rx is clearly larger than Rz is shown, which is consistent with the toroidal structure of the B=2 skyrmion. The RMS charge radii for the gSKM are larger than those for the SKM, and the calculated results for the ground state are similar to those for the excited state. Furthermore, Fig. 2 shows that R, Rx, and Rz decrease with increasing magnetic field strength, which indicates a gathering of the electric charges when the magnetic fields increase.

      Figure 2.  (color online) |eB| dependence of the RMS electric charge radii (labeled as R) of the deuteron, Rx and Rz for the x and z directions, in the ground (color black) and excited states (color red) for the gSKM (solid curves) and SKM (dashed curves).

      The results on the RMS electric charge radii in Fig. 2 can also be understood from the profile functions F(r), which are shown in Fig. 3. Figure 3 shows that the profile functions under magnetic fields for the gSKM are wider than those for the SKM, which is consistent with the larger RMS radii for the gSKM in Fig. 2. Moreover, the phenomenon that the increase in |eB| narrows F(r) also indicates a gathering of the electric charges when the magnetic fields increase.

      Figure 3.  (color online) Profile functions F(r) for the deuteron in the ground state.

      In Fig. 4, we illustrate the |eB| dependence of the mass for the deuteron. It is shown that the calculated masses in the absence of magnetic fields are 1750 and 1760 MeV for the SKM and gSKM, respectively, for the ground state 3S1, and 1783 and 1798 MeV for the SKM and gSKM, respectively, for the excited state 1S0 (all of which are small compared with the experimental values), and 1876 and 1880 MeV for the 3S1 and 1S0 states. With increasing magnetic field strength, the calculated masses are shown to increase, and more specifically, they increase faster under the strong magnetic fields in the gSKM than in the SKM. This is mainly a result of the continuously increasing Md,st, the static contribution to the deuteron mass, defined in Eq. (10), the increasing velocity with |eB| of which is larger in the gSKM. Even though the kinematical contribution Md,kin, defined in Eq. (22), plays a minor role, it can still be seen that in the gSKM, Md,kin, that is, the difference between the thick and thin curves in Fig. 4, is smaller under strong magnetic fields than that in the SKM. Physically speaking, Md,kin can be regarded as the rotational energy in the deuteron mass, and the small kinematical energy in the gSKM can be interpreted as the large moment of inertia, which might conform with the phenomenon that Rx decreases with |eB| more slowly in the gSKM than in the SKM, as shown in Fig. 2.

      Figure 4.  (color online) |eB| dependence of the masses of the deuteron in the ground (color black) and excited states (color red) for the gSKM and SKM. The thin curves are only for the static contribution Md,st to the deuteron mass, defined in Eq. (10). The thick curves are for the deuteron mass Md, which is the sum of the static and kinematical contributions.

    IV.   SUMMARY
    • In summary, we modify the method proposed in Ref. [8] to adapt to the investigation for the deuteron under strong magnetic fields. Besides the conventional Skyrme model, the effects of the derivative dependent sextic Lagrangian term are also considered. The main conclusions are listed below.

      (i) The deviation of the calculated results in the presence of the magnetic fields from the cRMA, that is, cz=1 in this study, is 10 %.

      (ii) The phenomenon that the RMS radii decrease with increasing magnetic fields indicates a gathering of the electric charges when the magnetic fields increase.

      (iii) The static contribution of the deuteron mass increases with increasing magnetic field strength, and the increasing velocity is larger for the gSKM. This results in the continuous increasing of the mass with |eB|, which is faster for the gSKM than for the SKM. The kinematical contribution plays a minor role and is smaller under strong magnetic fields in the gSKM than that in the SKM.

      Unlike in Ref. [8], where only the spatial rotation in the xy plane is considered when the skyrmion is quantized, we consider three-dimensional rotation in spatial space as well as in isospin space in this study. This quantization method can extend to the investigation of the B>2 skyrmion under magnetic fields more conveniently than that used in Ref. [8] because both spin and isospin rotations can be considered independently, and the quantization rules for skyrmions, the so-called Finkelstein-Rubinstein constraints [28, 29], which constrain the spins and isospins of quantum states, are satisfied more easily. This method can also be used to study the B=1 skyrmion. These issues are interesting and may be explored in detail in future studies.

    APPENDIX
    • The equations needed in Eqs. (17)–(19) are listed below.

      a10=drsin2(F)6c2z{2[8+3π+c3z(163π)]sin2(F)+[329π+4c3z(8+3π)]r2F2+6c2zπr2}+κsin4(F)F23cz[8+3π+3c3z(4π)],

      a11=drr2R3S3{(4π)sin4(F)+A6c2z[28+9π+4c3z(103π)]sin4(F)F2},

      a12=drr4sin4(F)24cz{c3z(1S23)[3058R23+9π(1+2R23)]2(1+S23)[919R23+π(3+6R23)]}+κsin4(F)F2240c3z{c3z{5[26+9π+S23(86+27π)]2[139+45π+3S23(143+45π)]R23}4S23[95+189R23+π(3060R23)]2c6z(1+S23)[95189R23+30π(1+2R23)]},

      a30=drsin2(F)3c2z{6c2zr2(4π)+2[329π+c3z(28+9π)]sin2(F)+[28+9π+4c3z(103π)]r2F2}+4κ3czsin4(F)F2,

      a32=drr4sin4(F)24cz(1R23)[2(1+S23)(103π)c3z(1S23)(289π)]κsin4(F)F2240c3z(1R23){8S23(4715π)+4c6z(1+S23)(47+15π)+c3z[148+45π+S23(428+135π)]},

      b10=dr215c3z{5sin2(F){2[40+c3z(5618π)+12π+c6z(10+3π)]r2+cz[64+18π6c3z(28+9π)+c6z(104+33π)]sin2(F)}+r2F2[2(1+c3z)2r2+5cz[16+c3z(6621π)+6π+2c6z(28+9π)]sin2(F)]}+κsin4(F)F23cz[8+c3z(8+3π)],

      b11=drr2R3S3{415c3zsin2(F){5[96+c3z(16051π)+30π+3c6z(22+7π)]sin2(F)+2(1+c3z)2(47+15π)r2F2}+κ3c2z[329π+c3z(28+9π)]sin4(F)F2},

      b12=drr4sin2(F)3360c4z{14{32c3z[3301178R23+15π(7+25R23)]+4(1+S23)[8002976R23+15π(17+63R23)]+2c9z(1+S23)[18405704R2315π(39121R23)]+c6z[5[2176693π+(7360+2343π)R23]+S23[45π(43153R23)+64(95+338R23)]]}sin2(F)+4(1+c3z)2[6c3z(1+S23)[328105π+(988+315π)R23]+(1+S23)[16564952R23+525π(1+3R23)]]r2F2}+κsin4(F)F2480c3z{4c9z(1+S23)(23675π)(1+R23)8c6z{4(2+S23)(47+15π)+[706+S23(518165π)+225π]R23}4[328105π+(808+255π)R23+3S23(16+5π)(1+9R23)]c3z{3536+8336R23+45π(2559R23)+S23[1552495π+(3152+1005π)R23]}}.

Reference (29)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return