中国科学院大学学报  2024, Vol. 41 Issue (1): 35-41   PDF    
光生反应γp→π-Δ++中的核子激发态
朱一鸣, 杨富中     
中国科学院大学物理科学学院, 北京 100049
摘要: 运用有效拉氏量方法分析LEPS合作组首次发布的关于γp → π-Δ++反应的高精度微分截面和极化观测量Beam asymmetry(Σ)的实验数据。把t道的π和ρ交换,u道的Δ交换,s道的核子N交换以及使光生反应振幅满足规范不变性的推广接触项作为背景项,尝试在s道进一步引入核子激发态交换的贡献以描述此实验数据。分析表明,如果在模型计算中不引入s道核子激发态交换的贡献,背景项的理论计算结果虽然能够大体上描述γp→π-Δ++反应的微分截面的实验数据,但是无法描述极化观测量Σ的实验数据;在s道中引入核子激发态N(1860)5/2+,模型计算结果就能够很好地同时描述LEPS实验组的微分截面和极化观测量Σ的实验数据。
关键词: 光生反应    核子激发态    有效拉氏量方法    
Nucleon resonances in γp→π-Δ++ photoproduction
ZHU Yiming, YANG Fuzhong     
School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
Abstract: The first measured high-precision data on the differential cross sections and polarization observable beam asymmetry Σ for the γp→π-Δ++ reaction from the LEPS collaboration are analyzed within a tree-level effective lagrangian approach. In addition to the t-channel π and ρ exchanges, the u-channel Δ exchange, the s-channel N exchange, and the interaction current that is required by the gauge invariance, the s-channel nucleon resonance exchanges are further taken into account to reproduce the data of LEPS. Numerical results show that, the data on Σ can not be well reproduced when no contributions from the s-channel resonance exchanges are introduced. Further analyses show that, by including the contribution from the resonance N(1860)5/2+, both the data on the differential cross sections and beam asymmetry Σ can be well reproduced.
Keywords: photoproduction reaction    nucleon resonances    effective lagrangian approach    

研究强子的性质,例如研究强子的质量,宽度,形状因子等,可以帮助进一步理解强子的内部结构以及量子色动力学(QCD)的非微扰特性。当前,实验上观测到的核子激发态数目要远远少于夸克模型计算所预言的激发态数目,这就是强子谱研究中的“丢失的激发态”问题。早期关于核子激发态的信息主要来自πN散射实验以及单π产生的光生反应实验。那些丢失的激发态可能是因为和πN的耦合比较弱,因此没有在传统的πN散射实验以及单π产生的光生反应实验中被观察到。因此,不少实验和理论工作已经尝试在不包含πN末态的强子产生反应过程中寻找那些丢失的激发态,例如对矢量介子,η′介子和KY*(Y = Σ,Λ)光生反应的研究[1-15]

本工作感兴趣的反应过程是πN光生反应,即γp→π-Δ++反应过程。2018年,LEPS合作组发布此反应的高精度实验数据,其中包括散射角度0.7<cos θ<1,光子能量Eγ=1.5~2.95 GeV的微分截面和极化观测量beam asymmetry(Σ)的实验数据[14]。目前还未有已发表的理论工作对此实验数据进行分析。

文献[14]中,作者给出γp→π-Δ++反应的最新实验数据与早期基于树图近似的Regge模型的理论预言[15]的对比。对比发现,理论预言[15]大体上能描述LEPS实验组最新的微分截面数据,但是对关于极化观测量Σ的数据描述较差。文献[15]对γp→π-Δ++的模型计算只考虑了包含s道核子N交换,t道介子π和ρ交换,u道重子Δ交换以及接触项的背景项的贡献,并没有考虑来自s道核子激发态的贡献,微分截面的实验数据比模型理论预言要大,意味着可能需要在模型中引入s道核子激发态交换的贡献才能够描述实验数据。

本工作利用有效拉式量的方法分析LEPS合作组最新发布的γp→π-Δ++反应的高精度微分截面和极化观测量Σ的实验数据[14]。除文献[15]中所考虑的背景项以外,尝试进一步在s道引入核子激发态交换的贡献以同时描述关于微分截面和极化观测量Σ的实验数据。本工作的主要目的是结合此实验数据,研究该反应过程中核子激发态的贡献,提取相关的激发态的信息以及分析激发态在该反应中的作用机制。

1 理论模型

依据文献[16-17]中所描述的Full Field Theory,反应γp→π-Δ++的总振幅表达为

Mμν=Mμνs+Mμνt+Mμνu+Mμνint. (1)

前3项MsμνMtμνMuμν分别代表的是反应s道、t道、u道振幅的贡献,最后1项Mintμν为接触项。总振幅表达式中的4项在图 1中以费曼图的形式表现出来。如图所示,将对γp→π-Δ++反应过程考虑如下的交换:(ⅰ) s道的核子N与核子激发态N*交换,(ⅱ) t道的介子π和ρ交换,(ⅲ) u道的重子Δ交换。

Download:
图 1 光生反应γp→π-Δ++的费曼图 Fig. 1 Generic structures of the photoproduction amplitudes for γp→π-Δ++

依据文献[16-19]对Mintμν项的处理方法,把Mintμν写成如下形式

Mνμint=ΓνπNΔ(k2)Cμ+MνμKRft, (2)

其中ΓπNΔν是由拉氏量LπNΔ得到的πNΔ的耦合顶角,

ΓνπNΔ(k2)=gπNΔmπk2ν. (3)

Cμ为非奇异的辅助流函数,它能够保证总振幅的规范不变,具体形式表达为

Cμ=QπftˆFtk22(2k2k1)μQΔfuˆFup22×(2p2k1)μQNfsˆFsp21(2p1+k1)μ, (4)
ˆF=1ˆh(1fs)(1fu)(1ft). (5)

其中p1p2k1k2分别为入射核子N,出射重子Δ,入射光子γ,出射介子π的4动量; QNQπQΔ分别为N,π-和Δ++的电荷;fsfuft分别是s道N交换, u道Δ交换和t道π交换的形状因子;ˆh是任意一个在高能极限下归一的函数,为了计算简便,把它定为1。MKRνμ是反应γp→π-Δ++的Kroll-Ruderman项:

MμνKR=gπNΔmπgμνQπ. (6)

本节的剩下部分给出本工作使用的有效拉氏量,核子激发态的传播子和形状因子。

1.1 有效拉氏量

本工作使用的拉氏量如下,为了符号上的简便记

Γ(+)=γ5 和 Γ()=1 (7)

计算非核子激发态部分费曼图的电磁相互作用的拉氏量为

LγNN=eˉN[(ˆeγμˆκN2mNσμvv)Aμ]N, (8)
Lγπρ=egγπρmπεαμλv(αAμ)(λπ)ρv, (9)
Lγππ=ie[π+(μπ)π(μπ+)]Aμ, (10)
LγΔΔ=QΔi2ˉΔμ{σμν,γλ}AλΔv (11)

这里ˆe为电荷算符; ˆκN=κp(1+τ3)/2+κn(1τ3)/2,反常磁矩κp=1.793和κn=-1.913;εαμλν是全反称Levi-Civita张量,约定ε0123=1。介子-介子电磁耦合常数gγπρ=0.11取自文献[19]。

介子-重子作用的有效拉氏量为

LρNΔ=igρNΔ2mNˉΔμγvγ5TρμνN+H.c., (12)
LπNΔ=gπNΔmπˉΔμT(μπ)N+H.c.. (13)

其中耦合常数gπNΔ=2.13,gρNΔ=-39.1取自文献[19]。激发态-核子-光子相互作用的有效拉氏量为

L1/2±RNγ=eg(1)RNγ2mNˉRΓ()σμν(vAμ)N+H.c., (14)
L3/2±RNγ=ieg(1)RNγ2mNˉRμγvΓ(±)FμνN+eg(2)RNγ(2mN)2ˉRμΓ(±)FμνvN+H.c., (15)
L5/2±RNγ=eg(1)RNγ(2mN)2ˉRμαγvΓ()(αFμν)N± ie g(2)RNγ(2mN)3ˉRμαΓ()(αFμν)vN+H.c., (16)
L7/2±RNγ=ieg(1)RNγ(2mN)3ˉRμαβγvΓ(±)(αβFμν)Neg(2)RNγ(2mN)4ˉRμαβΓ(±)(αβFμv)vN+H.c.. (17)

这里R表示核子激发态的场,拉氏量LRNγ的上标表示对应激发态的自旋与宇称。拉氏量中gRNγ(i)(i=1,2)为激发态-电磁耦合常数。激发态-强子相互作用的有效拉氏量为

L1/2±πRΔ=g(1)πRΔmπˉΔμΓ()(μπ)R+H.c., (18)
L3/2±πRΔ=g(1)πRΔmπˉΔμγαΓ(±)(απ)Rμ+ig(2)πRΔm2πˉΔμΓ(±)(μαπ)Rμ+H.c., (19)
L5/2±πRΔ=ig(1)πRΔm2πˉΔαγμΓ()(μβπ)Rαβg(2)πRΔm3πˉΔμΓ()(μαβπ)Rαβ+H.c., (20)
L7/2±πRΔ=g(1)πRΔm3πˉΔαγμΓ(±)(μβλπ)Rαβλig(2)πRΔm4πˉΔμΓ(±)(μαβλπ)Rαβλ+H.c.. (21)

其中gπRΔ(i)(i=1,2)为激发态-强子耦合常数。拟合实验数据过程中把gπRΔ(1)与激发态-电磁耦合参数gRNγ(i)(i=1,2)的乘积作为自由参数,因为这些参数在振幅的表达式中以乘积的形式存在。为了计算上的简便并减少自由参数的数目,本工作略去包含gπRΔ(2)的项。

1.2 核子激发态的传播子

自旋1/2的激发态的传播子为

S1/2(p)=iMR+iΓR2 (22)

依据文献[20-22],本工作中自旋3/2,5/2与7/2的核子激发态的传播子取下式的形式

S3/2(p)=iMR+iΓR2(˜gμν+13˜γμ˜γv), (23)
S5/2(p)=iMR+iΓR2[12(˜gμα˜gvβ+˜gμβ˜gvα)15˜gμ˜gαβ+110(˜gμα˜γv˜γβ+˜gμβ˜γv˜γα+˜gvα˜γμ˜γβ+˜gvβ˜γμ˜γα)], (24)
S7/2(p)=iMR+iΓR2136PμPv(˜gμ1v1˜gμ2v2˜gμ3v37˜gμ1μ2˜gv1v2˜gμ3v3+37˜γμ1˜γv1˜gμ2v2˜gμ3v3335˜γμ1˜γv1˜gμ2μ3˜gv2v3) (25)

其中˜gμν˜γμ分别是

˜gμν=gμν+pμpvm2R, (26)
˜γμ=γv˜gνμ=γμ+pμpm2R. (27)

式(25) 中对Pμ(Pν)的求和取遍μ1μ2μ3(ν1ν2ν3)所有可能的排列(共有3!= 6项)。

1.3 费曼振幅的形状因子

由拉氏量得到的强相互作用顶角需要添加一个唯象的形状因子使强子的结构参数化。依照文献[15],本工作使用如下的形状因子

f(p2)=Λ4Λ4+(p2m2)2. (28)

其中pm分别代表中间态交换粒子的4动量和质量。s道核子N,t道ρ和u道Δ交换的截断质量在本文中定为450 MeV,t道π和s道激发态N*的截断质量作为自由参数拟合实验数据。注意到计算过程中最终的总振幅的规范不变性并不依赖于形状因子的具体表达形式。

2 计算结果与讨论

2018年LEPS实验组发布了关于光生反应γp→π-Δ++的前角部分(0.7<cos θ<1,这里θ为质心系下的散射角度)光子能量Eγ=1.5~2.95 GeV的微分截面与极化观测量Σ的高精度实验数据[14]。其中,微分截面的实验数据点较多,有145个实验数据点,而极化观测量Σ的实验数据点较少,有56个实验数据点。在此实验结果发布之后迄今还未有已发表的理论工作对此实验数据进行分析。文献[14]中展示了2011年理论工作[15]对反应γp→π-Δ++所做的理论预言与此实验数值的比较,其中微分截面部分与实验数值大体符合,散射角度0.966<cos θ<1处,实验值比理论值高,极化观测量Σ在散射角度0.7<cos θ<0.9处理论对实验的描述较差,可能是由于文献[15]的理论模型中没有考虑来自核子激发态的贡献。

本工作基于有效拉氏量的方法,分析LEPS实验组发布的极化观测量Σ和微分截面的实验数据,构建同时描述这2个观测量实验数据的理论模型。特别地,将在s道加入可能的核子激发态以使模型对实验的描述更好。使用χ2/N表征理论值与实验值的偏差,值越小表示计算结果与实验值越接近,同时考虑总的包括微分截面与极化观测量Σ的χ2/N和单独极化观测量Σ的χ2/N。首先用不加入任何核子激发态的背景项的理论模型来拟合LEPS实验组的实验数据。图 2(a)2(b)分别展示微分截面与极化观测量Σ的拟合结果。对于微分截面,散射角度0.7<cos θ<0.966处理论值与实验值符合得很好;散射角度0.966<cos θ<1,光子能量Eγ=1.5~1.8 GeV处,理论值比实验值低。对于极化观测量Σ,散射角度0.7<cos θ<0.9处,理论描述较差。此时总的χ2/N=13.2,极化观测量Σ的χ2/N为14.1。

Download:
图 2 微分截面与极化观测量Σ Fig. 2 Differential cross sections and beam asymmetry Σ

进一步尝试加入核子激发态以期使模型同时描述微分截面与极化观测量Σ的实验数据。对粒子数据表上的所有核子激发态一一进行尝试,加入它们之后拟合并且计算拟合之后的结果。在此过程中,发现加入核子激发态的模型与只考虑背景项的模型相比,对微分截面实验数据的描述改善较小,散射角度0.966<cos θ<1,光子能量Eγ=1.5~1.8 GeV处,只考虑背景项拟合的结果与实验的偏差在加入激发态拟合之后略微减少;对背景项加入核子激发态N(1860)5/2+能够很大程度上改善极化观测量Σ的结果。加入激发态N(1860)5/2+后拟合得到的总χ2/N=7.6且极化观测量Σ的χ2/N由只考虑背景项时的14.1降低到5.4,而加入其他激发态后拟合得到的这2个χ2/N都会比加入N(1860)5/2+后得到的大。加入N(1860)5/2+拟合的参数见表 1。表中激发态的名字下面表示PDG对此激发态评测的星级[23](粒子数据表上的激发态被PDG评测为从1星到4星4个星级,星数越多代表对此激发态的研究越充分),MRΓR分别代表此激发态的质量和宽度,ΛR和Λπ分别代表s道激发态和t道π交换的截断质量,s道N交换,t道ρ交换和u道Δ交换的截断质量在本工作中定为450 MeV。gRNγgπRΔ为出现在拉氏量LRNγLπRΔ中的耦合常数,本工作中拟合的参数为这些耦合常数的乘积,计算中出现的其他参数在上一节已给出。下面对模型拟合的计算结果进行讨论。

表 1 反应γp→π-Δ++模型拟合的参数 Table 1 Fitted parameters used in the model for the reaction γp→π-Δ++

对于微分截面,如上所述,加入激发态N(1860)5/2+的拟合结果与单纯背景项的拟合结果差别不大。图 2(c)实线展示的是对背景项加入N(1860)5/2+后模型在固定散射角度下微分截面随光子能量Eγ的变化,图中画出主要项各自的贡献。其中实线代表总的包括背景项与N(1860)5/2+的贡献,虚线、点线与虚点线分别表示来自激发态、t道π交换与接触项的贡献,实验值取自文献[14]。对于微分截面,t道π交换与接触项占主要贡献,s道N(1860)5/2+交换的贡献也很大,越前角的部分,N(1860)5/2+交换的贡献越小。散射角度0.966<cos θ<1,光子能量Eγ=1.5~1.8 GeV处,N(1860)5/2+对于微分截面有一定的提升作用,使得理论计算结果与实验值更加接近,但与实验值还是有一些偏差。

对于极化观测量Σ,图 2(d)中实线展示的是加入激发态N(1860)5/2+后的模型在固定散射角度下极化观测量Σ随光子能量Eγ的变化。图中点线、虚点线、虚线分别代表在此拟合中去掉t道π交换、接触项和N(1860)5/2+时的计算结果。从图中看出,s道加入激发态N(1860)5/2+会使模型对极化观测量Σ实验数据的理论描述有很大的改进。如上述对微分截面部分的讨论一致,对于极化观测量Σ在本工作中所考虑的散射角度0.7<cos θ<1内t道π交换和接触项占主要贡献,N(1860)5/2+在其中的贡献也不可忽略(点线、虚线、虚点线与实线的偏差很大)。越靠近前角(cos θ~1)的部分,激发态的贡献越小(虚线与实线越在前角的部分越接近)。

3 总结

本工作基于有效拉氏量方法,在树图层次上分析2018年LEPS实验组发布的关于γp→π-Δ++反应的高精度微分截面和极化观测量Σ的实验数据[14]。s道核子N交换,t道介子π和ρ交换,u道重子Δ交换以及接触项作为背景项,尝试在s道进一步引入核子激发态N*交换的贡献以同时描述关于微分截面和极化观测量Σ的实验数据。结果表明,只考虑背景项贡献的模型计算虽然大体上能够描述微分截面的实验数据,但是无法描述极化观测量Σ的实验数据。通过在s道引入激发态N(1860)5/2+交换,理论模型就能很好地同时描述微分截面和极化观测量Σ的实验数据。进一步数值计算表明,t道π交换与接触项在微分截面和极化观测量Σ中占主要贡献,N(1860)5/2+的贡献也很大,越靠近前角的部分,贡献越小。与只考虑背景项的结果相比,加入激发态N(1860)5/2+之后,散射角度0.966<cos θ<1,光子能量Eγ=1.5~1.8 GeV处微分截面有一定的提升,但与实验值还是有一些偏差。文献[14]中的实验数据虽然精度很高,但只包含散射角度为前角部分(散射角度0.7<cos θ<1)的数据,为获得对反应γp→π-Δ++中激发态N(1860)5/2+在其中的作用机制更全面的分析和了解,需要在分析中把更宽广的散射角度和能量范围的实验数据考虑进来做进一步研究。

感谢黄飞教授与韦能昌博士在研究过程中给予的指导和建议。
参考文献
[1]
Anisovich A V, Burkert V, Dugger M, et al. Proton-η' interactions at threshold[J]. Physics Letters B, 2018, 785: 626-630. Doi:10.1016/j.physletb.2018.06.034
[2]
Kim S H, Hosaka A, Kim H C. Effects of N(2000) 5/2+, N (2060) 5/2-, N (2120) 3/2-, and N (2190) 7/2- on K*Λ photoproduction[EB/OL]. (2014-02-07)[2021-04-10]. http://arxiv.org/abs/1402.1753.
[3]
Moriya K, Schumacher R A, Adhikari K P, et al. Differential photoproduction cross sections of the Σ0(1385), Λ(1405), and Λ(1520)[J]. Physical Review C, 2013, 88(4): 045201. Doi:10.1103/PhysRevC.88.045201
[4]
Tiator L, Gorchtein M, Kashevarov V L, et al. Eta and etaprime photoproduction on the nucleon with the isobar model EtaMAID2018[J]. The European Physical Journal A, 2018, 54(12): 210. Doi:10.1140/epja/i2018-12643-x
[5]
Wang A C, Wang W L, Huang F, et al. Nucleon resonances in γp → K*+Λ[EB/OL]. 2017: arXiv: 1704.04562. (2017-04-15)[2021-04-08]. http://arxiv.org/abs/1704.04562.
[6]
Wang A C, Wang W L, Huang F. Effects of Δ(1905)5/2+ on K*Σ photoproduction[EB/OL]. 2018: arXiv: 1810.07012. (2018-10-14)[2021-03-15]. http://arxiv.org/abs/1810.07012/.
[7]
Wang A C, Huang F, Wang W L, et al. Photoproduction γp → K*+Λ in a Reggeized model[J]. Physical Review C, 2020, 102: 015203. Doi:10.1103/physrevc.102.015203
[8]
Wang A C, Wang W L, Huang F. Nucleon and Δ resonances in γp → K+Σ0(1385) photoproduction[EB/OL]. (2020-02-11)[2021-03-10]. https://doi.org/10.48550/arXiv.2002.04213.
[9]
Wang E, Xie J J, Liang W H, et al. Role of a triangle singularity in the γp → K+Λ(1405) reaction[EB/OL]. 2016. arXiv: 1610.07117. (2016-10-23)[20210-03-20]. https://arxiv.org/abx/1610.07117.
[10]
Wei N C, Huang F, Nakayama K, et al. Nucleon resonances in γp → ωp reaction[J]. Physical Review D, 2019, 100(11): 114026. Doi:10.1103/physrevd.100.114026
[11]
Wei N C, Wang A C, Huang F, et al. Analysis of the data on spin density matrix elements for γp → K*+Λ[J]. Physical Review C, 2020, 101: 014003. Doi:10.1103/physrevc.101.014003
[12]
Wei N C, Zhang Y, Huang F, et al. Photoproduction γp→K+Λ(1520) in an effective Lagrangian approach[J]. Physical Review D, 2021, 103(3): 034007. Doi:10.1103/physRevD.103.034007
[13]
Zhang Y, Huang F. Effects of N(2000)5/2+ on γp → K+Λ(1405)[EB/OL]. 2021: arXiv: 2102.10404. (2021-02-20)[2021-04-10]. https://arxiv.org/abs/2102.10404.
[14]
Kohri H, Shiu S H, Chang W C, et al. Differential Cross Section and Photon-Beam Asymmetry for the γp→π-Δ++(1232) Reaction at Forward π- Angles for Eγ=1.5~2.95 GeV[J]. Physical Review Letters, 2018, 120(20): 202004. Doi:10.1103/PhysRevlett.120.202004
[15]
Nam S I, Yu B G. Charged pion photoproduction with the Δ(1232) baryon beyond the resonance region[J]. Physical Review C, 2011, 84(2): 025203. Doi:10.1103/physrevc.84.025203
[16]
Haberzettl H, Nakayama K, Krewald S. Gauge-invariant approach to meson photoproduction including the final-state interaction[J]. Physical Review C, 2006, 74(4): 045202. Doi:10.1103/physrevc.74.045202
[17]
Haberzettl H. Gauge-invariant theory of pion photoproduction with dressed hadrons[J]. Physical Review C, 1997, 56(4): 2041-2058. Doi:10.1103/physrevc.87.054004
[18]
Huang F, Haberzettl H, Nakayama K. Combined analysis of η' production reactions: γN→η'N, NN→NNη', and πN → η'N[J]. Physical Review C, 2013, 87(5): 054004. Doi:10.1103/physrevc.87.054004
[19]
Huang F, Döring M, Haberzettl H, et al. Pion photoproduction in a dynamical coupled-channels model[J]. Physical Review C, 2012, 85(5): 054003. Doi:10.1103/physrevc.85.054003
[20]
Behrends R E, Fronsdal C. Fermi decay of higher spin particles[J]. Physical Review, 1957, 106(2): 345. Doi:10.1103/physrev.106.345
[21]
Fronsdal C. On the theory of higher spin fields[J]. Il Nuovo Cimento (1955-1965), 1958, 9(2): 416-443.
[22]
Zhu J J, Yan M L. Covariant amplitudes for mesons[EB/OL]. arXiv: 9903349. (1999-03-15)[2020-10-12]. https://arxiv.org/abs/hep-ph/9903349.
[23]
Tanabashi M, Hagiwara K, Hikasa K, et al. Review of Particle Physics[J]. Physical Review D, 2018, 98(3): 030001. Doi:10.1103/physrevd.98.030001
光生反应γp→π-Δ++中的核子激发态
朱一鸣, 杨富中