-
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
1015 G, whereas in the inner core of magnetars, the magnitude of the magnetic fields could be of the order of1018 G [1]. Second, it is considered that during the electroweak phase transition in the early universe, extremely strong magnetic fields exceeding1020 G might have been produced [2]. In the course of relativistic non-central heavy-ion collisions, the strength of magnetic fields can reach up to1018 G at the Relativistic Heavy Ion Collider (RHIC) and1020 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 [5–7]. However, there are relatively few studies in regard to the modification of baryon properties in the presence of strong magnetic fields [8–10].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. [12–18] 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 theSU(2) valued field in the Skyrme model to rotate in three-dimensional isospin space and spatial space only in the x – y 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 theB=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. -
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π2∫d5xεijklmTr[RiRjRkRlRm]+116π2∫d4xεμναβTr{eAμQE(RνRαRβ+LνLαLβ)+ie22[AμQERνAαQERβ−AμQELνAαQELβ]−ie2[AμQEU†AνQEURαRβ−AμQEUAνQEU†LαLβ−∂μAνQEU†AαQEURβ−∂μAνQEUAαQEU†Lβ−(AμQE∂νAαQE+∂μAνQEAαQE)(Rβ+Lβ)]+e3AμQEAνQEAαQE(Rβ+Lβ)+e3(AμQEU†AνQEUAαQERβ+AμQEUAνQEU†AαQELβ+AμQEU†AνQE∂αAβQEU+AμQEU†∂νAαQEAβQEU+∂μAνQEAαQEU†AβQEU+AμQE∂νAαQEU†AβQEU)−ie4(AμQEAνQEAαQEU†AβQEU+AμQEU†AνQEAαQEAβQEU+12AμQEU†AνQEUAαQEU†AβQEU)},
(1) where the first term is a five-dimensional integral whose boundary is four-dimensional space-time,
Ri=U†∂iU ,Li=∂iUU† , e is the unit electric charge, andQE=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μU†DμU)+m2πF2π8Tr(U−1)+132g2Tr([U†DμU,U†DνU]2)−9π25m2BμBμ,
(2) with the covariant derivative
DμU defined asDμU=∂μU−ieAμ[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 valuemπ=138 MeV. In this study, the homogeneous magnetic field B along the z direction is considered, and the symmetric gauge of the electromagnetic fieldAμ is chosen asAμ=(0,−yB/2, xB/2, 0) . As anSU(2) valued field, U is expressed asU=exp(iF(r)→τ⋅→n(θ,φ)),
(4) where
F(r) is a real profile function satisfying the boundary conditionsF(0)=π andF(∞)=0 ,|→n|=1 . For baryon numberB>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 asz=tan(θ/2)exp(iφ) . The form of R isp/q , where p and q are polynomials in z with the larger of their degrees being B. For theB=2 skyrmion,R=z2 is chosen. The last term in Eq. (2) is sextic in derivatives, in whichBμ reads asBμ=124π2εμναβTr[(U†DνU)(U†DαU)(U†Dβ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ρ andcz are positive dimensionless parameters, andθ∈[0,π] andφ∈[0,2π] are polar angles. Then, the volume element isdV=dxdydz=c2ρczr2sin(θ)drdθdφ . The parameterscρ andcz have two effects: they deform the shape and scale the volume. Because the scale effect ofcρ andcz can be absorbed by performing the scale transform of r, we follow the choice in Ref. [8] in whichcρ≡1/√cz to seek results with a conserved volume. It should be mentioned that in Ref. [8], which studiesB=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), andcz≠1 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, satisfyingA−1˙A=i2ωaτa,(˜S−1˙˜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 asMd,st=⟨ψd∣M0+M1(eB)+M2(eB)2∣ψd⟩,
(10) in which
M0=(πβg2)∫dr12c2zr2{3cz[16+c3z(−8+3π)]+48[4+c3z(−2+π)−π]r2+96c2zm2πr4−96c2zm2π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)}F′2+4(8+3π)κc2zsin4(F)F′2},
(11) M1=(R3S3πβg2)∫2dr3c2z{sin2(F){−6c2z(−4+π)r2−2c3z(28−9π)sin2(F)−2(−32+9π)sin2(F)−[28−9π+4c3z(−10+3π)]r2F′2}+4κczsin4(F)F′2},
(12) and
M2=(πβg2)∫r2dr480c3z{sin2(F){−40c2z{2(1+S23)[8−3π+(−28+9π)R23]+c3z(−1+S23)[32−9π+(−88+27π)R23]}r2+20sin2(F){8S23(−28+9π+124R23−39πR23)+4c6z(−1+S23)[28−9π+2(−47+15π)R23]+c3z[128−39π+(−448+141π)R23+S23[80−27π+(−400+129π)R23]]}−4{−8S23[48−142R23+15π(−1+3R23)]+4c6z(−1+S23)[48−142R23+15π(−1+3R23)]+c3z[112−45π+3(−136+45π)R23+S23[432−1288R23+135π(−1+3R23)]]}r2F′2}+10κczsin4(F)F′2{2S23[−8+3π+3(24−7π)R23]+c3z(1−S23)[16−3π+3(−24+7π)R23]}}.
(13) Here,
β=2√2/(gFπ) is used as the unit of length,κ=9g2/(5m2β2π2) ,R3 is defined asR3=Tr(τ3A†τ3A)/2 , andS3 is the˜S33 component of the rotation matrix˜S .The replacement of U with
ˆU can also give the kinematical part in the formT=12aiω2i+12biΩ2i−ciωiΩi+D(eB)S3(ω3−2Ω3)+E(eB)(S1Ω1+S2Ω2),
(14) where
i=1,2,3 ,S1 andS2 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 , anda3=b3/4=c3/2 , wherea1=(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, asKa=∂T∂ωa,La=∂T∂Ωa.
(20) Note that in Ref. [8], only the spatial rotation in the x – y plane is considered, that is, only the third component of the spatial angular momentum is defined and
k≡3 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 fieldeB 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 andk3 ) ofJ3 andK3 can simultaneously be half-integers for theB=1 nucleon. However, the integerj3=1 for theB=2 deuteron allowsk3 to be a half-integer, according to Eq. (21), which is prohibited in the quantization of theB=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 andL3 become mutually independent. Consequently, the quantum numbers (l3 andk3 ) ofL3 andK3 are allowed to be zero for theB=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 andi3 are the quantum numbers of the space-fixed isospin momentum (I2 andI3 , physically more meaningful).Now, using the relations
J3=−∑SiLi and∑S2i≡1 , wherei=1,2,3 , we can write the kinematical contribution to the deuteron mass asMd,kin=⟨ψd∣∑a=1,2,3(ωaKa+ΩaLa)−T∣ψd⟩,=⟨ψd∣I22a1+J2+(1−S23)E2(eB)22b1+[(K3−S3D(eB))22a3−(1a1+4b1)K232]+(J3−2S3K3)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 andMd,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 minimizingMd,st , with the help of a certain value ofcz obtained by minimizingMd .The magnetic moment of the deuteron is defined by
μ≡−∂Md/∂(eB) . Becausek3 is equal to zero in the wave function of the deuteron, and the expectation value ofR3 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 asg= 4.84 andFπ= 108 MeV for the SKM, andg= 19.48 andFπ= 129.11 MeV for the gSKM with an additional parameter m chosen asm= 420 MeV [19]. It is known that the wave function of the deuteron, as the ground state3S1 of theB=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. -
In Fig. 1, we illustrate the
|eB| dependence ofcz for the deuteron. We also give the calculations for the first excited state1S0 , 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, withcz=1.03 in the SKM andcz=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 ofcz with|eB| than SKM. Overall, we see that the deviation of the calculated results fromcz=1 in the presence of the magnetic fields is≲10 %.Figure 1. (color online)
|eB| dependence ofcz 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′π+⟨ψd∣r(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 chargeB=2 of the deuteron, and the second termρB(r) is induced by magnetic fields. ForB=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 thatk3=0 and the expectation value ofR3 in Eq. (23) is zero [20] indicates thatρB≡0 in the deuteron. Becausei3=0 for the deuteron, we can evaluate the root mean square (RMS) electric charge radii, defined asR=(∫X2ρ0/2dV)1/2 , in whichX=√x2+y2+z2 , andRx=(∫x2ρ0/2dV)1/2 andRz=(∫z2ρ0/2dV)1/2 for directions of x and z, respectively. The results are illustrated in Fig. 2. The phenomenon thatRx is clearly larger thanRz 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 , andRz 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 andRz 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| narrowsF(r) also indicates a gathering of the electric charges when the magnetic fields increase.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 state3S1 , and 1783 and 1798 MeV for the SKM and gSKM, respectively, for the excited state1S0 (all of which are small compared with the experimental values), and 1876 and 1880 MeV for the3S1 and1S0 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 increasingMd,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 contributionMd,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 thatRx 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 contributionMd,st to the deuteron mass, defined in Eq. (10). The thick curves are for the deuteron massMd , which is the sum of the static and kinematical contributions. -
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 x – y 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 theB=1 skyrmion. These issues are interesting and may be explored in detail in future studies. -
The equations needed in Eqs. (17)–(19) are listed below.
a10=∫drsin2(F)6c2z{2[−8+3π+c3z(16−3π)]sin2(F)+[32−9π+4c3z(−8+3π)]r2F′2+6c2zπr2}+κsin4(F)F′23cz[−8+3π+3c3z(4−π)], a11=∫drr2R3S3{(4−π)sin4(F)+A6c2z[−28+9π+4c3z(10−3π)]sin4(F)F′2},
a12=∫drr4sin4(F)24cz{c3z(1−S23)[30−58R23+9π(−1+2R23)]−2(1+S23)[9−19R23+π(−3+6R23)]}+κsin4(F)F′2240c3z{c3z{5[−26+9π+S23(−86+27π)]−2[−139+45π+3S23(−143+45π)]R23}−4S23[−95+189R23+π(30−60R23)]−2c6z(−1+S23)[95−189R23+30π(−1+2R23)]},
a30=∫drsin2(F)3c2z{6c2zr2(4−π)+2[32−9π+c3z(−28+9π)]sin2(F)+[−28+9π+4c3z(10−3π)]r2F′2}+4κ3czsin4(F)F′2,
a32=∫drr4sin4(F)24cz(1−R23)[2(1+S23)(10−3π)−c3z(1−S23)(28−9π)]−κsin4(F)F′2240c3z(1−R23){8S23(47−15π)+4c6z(−1+S23)(−47+15π)+c3z[−148+45π+S23(−428+135π)]},
b10=∫dr215c3z{−5sin2(F){2[−40+c3z(56−18π)+12π+c6z(−10+3π)]r2+cz[−64+18π−6c3z(−28+9π)+c6z(−104+33π)]sin2(F)}+r2F′2[2(−1+c3z)2r2+5cz[−16+c3z(66−21π)+6π+2c6z(−28+9π)]sin2(F)]}+κsin4(F)F′23cz[8+c3z(−8+3π)],
b11=∫drr2R3S3{415c3zsin2(F){−5[−96+c3z(160−51π)+30π+3c6z(−22+7π)]sin2(F)+2(−1+c3z)2(−47+15π)r2F′2}+κ3c2z[32−9π+c3z(−28+9π)]sin4(F)F′2},
b12=∫drr4sin2(F)3360c4z{−14{−32c3z[330−1178R23+15π(−7+25R23)]+4(1+S23)[800−2976R23+15π(−17+63R23)]+2c9z(−1+S23)[1840−5704R23−15π(39−121R23)]+c6z[5[2176−693π+(−7360+2343π)R23]+S23[45π(43−153R23)+64(−95+338R23)]]}sin2(F)+4(−1+c3z)2[6c3z(−1+S23)[328−105π+(−988+315π)R23]+(1+S23)[1656−4952R23+525π(−1+3R23)]]r2F′2}+κsin4(F)F′2480c3z{4c9z(−1+S23)(236−75π)(−1+R23)−8c6z{4(−2+S23)(−47+15π)+[−706+S23(518−165π)+225π]R23}−4[328−105π+(−808+255π)R23+3S23(−16+5π)(1+9R23)]−c3z{−3536+8336R23+45π(25−59R23)+S23[1552−495π+(−3152+1005π)R23]}}.
Deuteron under strong magnetic fields: Skyrme model study
- Received Date: 2022-09-27
- Available Online: 2023-03-15
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.