单玻色子交换势模型下对Y(4390)的研究

单玻色子交换势模型下对 Y (4390)的研究 Study of Y (4390)in the One-Boson-Exchange Potential Model first-author 黎栩，男，本科生，现从事奇特强子态方面的研究。E-mail: 1170674471@qq.com 黎栩，男，本科生，现从事奇特强子态方面的研究。E-mail: 1170674471@qq.com 1170674471@qq.com corresp 李新强, E-mail: xqli@mail.ccnu.edu.cn 李新强, E-mail: xqli@mail.ccnu.edu.cn xqli@mail.ccnu.edu.cn 华中师范大学 物理科学与技术学院，湖北 武汉 430079 华中师范大学 物理科学与技术学院，湖北 武汉 430079 College of Physics Science and Technology, Central China Normal University, Wuhan 430079, Hubei, China College of Physics Science and Technology, Central China Normal University, Wuhan 430079, Hubei, China 针对2016年BESⅢ合作组发现的新粒子 Y (4390)，采用单玻色子交换势(one-boson-exchange potential，OBEP)模型分析其组成；综合考虑了 π 、 σ 、 ρ 、 ω 、 η 介子的交换效应，轨道波函数 S - D 波的混合效应，以及交换势函数中δ项的贡献，给出 Y (4390)的完整单玻色子交换势。将 Y (4390)理解为量子数为 J PC =1 -- 的 D 1 D * / D 1 D * 系统，利用薛定谔方程求解了系统波函数的数值解，对系统的质量、方均根半径、组分等性质进行了预言，认为 Y (4390)可以看作是 J PC =1 -- 的 D 1 D * / D 1 D * 强子分子态。 In this paper, we study the composition of the new particle Y (4390) observed by the BESⅢ collaboration in 2016, within the framework of the one-boson-exchange potential (OBEP) model. Considering the exchange effects due to π , σ , ρ , ω , η mesons, the contributions of the orbital S - D wave function mixing, as well as the δ term in the exchange potential, the complete oneboson-exchange potential of Y (4390) is given. Interpreting the newly observed Y (4390) as a D 1 D * / D 1 D * system with the quantum number J PC =1 -- , then by solving numerically the Schrodinger equation, we predict the properties of the system, such as the mass, the root-mean-square radius and the proportions. Finally, we conclude that Y (4390) could be regarded as a D 1 D * / D 1 D * hadron molecular state with the quantum number J PC =1 -- . 强子分子态 单玻色子交换势 束缚态 截断参数 hadron molecular state one-boson-exchange potential bound state truncation parameter 引言 1 2016年10月，BESⅢ合作组 [ 1 1 ] 在 e + e - → ( h c + π + π - )过程中发现了2个新粒子 Y (4390)和 Y (4220)，并认为它们的量子数很可能均为1 -- ，其中， Y (4390)的质量和宽度分别被确认为 M =(4 392±6) MeV， Γ =(140±16) MeV。通过分析 Y (4390)的衰变产物，预言它可能是一个包含 c c 味道的重味矢量强子态。 Y (4390)属于近年来发现的大量奇特强子态之一，这些强子态就是被称为 X 、 Y 、 Z 的类粲偶素或类底偶素态 [ 2 2 ~ 4 4 2~4 ] ，其中很多不能用传统的夸克模型进行解释，内部的结构目前也尚无定论。对于这类强子态，常用的研究方法有: QCD求和规则 [ 5 5 ] 、势模型 [ 6 6 ] 、格点QCD [ 7 7 ] 等。这些方法将奇特强子态解释为强子分子态或多夸克态等。2017年，Wang [ 8 8 ] 提出使用QCD求和规则将 Y (4390)解释为由一个反介子 D 1 和一个 D 介子结合形成的 D 1 D 的分子态，但是通过该方法计算得到的 Y (4390)的质量却比分子态阈值(4 290 MeV)高很多，因此该解释不够合理。2018年，Wang [ 9 9 ] 根据数值计算结果提出 Y (4390)还可能为矢量四夸克态 c c q q 。另外，文献[ 10 10 ]中在单玻色子交换模型基础上用BetheSalpeter方程近似的方法进行分析，提出 Y (4390)是 D 1 D * 的分子态，但是此方法得到的分子态质量为4 384 MeV，仍然比实验上 Y (4390)的质量值低一些。目前关于 Y (4390)的结构仍无统一且令人满意的解释。 由于 Y (4390)与 D 1 D * 的阈值很接近，本文尝试从非相对论量子力学的角度出发，将 Y (4390)解释为量子数为1 -- 、结构为 D 1 D * / D 1 D * 的强子分子态。 D 介子是由一个粲夸克和一个轻夸克组成的，可以认为其内在是非相对论的。在非相对论的基础上分析 D 1 D * / D 1 D * 系统，即求解径向薛定谔方程，并使用单玻色子交换势(one-boson-exchange potential，OBEP)模型 [ 11 11 ] 分析方程中的相互作用势，求解 D 1 D * / D 1 D * 系统的波函数和本征能量，对 Y (4390)的性质进行解释和预言。 本文将首先描述OBEP模型的机制，给出标量介子、赝标介子、矢量介子交换势，接着给出模型中的数值计算参数，最后将模型应用到 Y (4390)中，对数值结果进行分析。 OBEP模型及其在 Y (4390)上的应用 1 本节将简要介绍OBEP模型及其相互作用图像并根据 Y (4390)在 D 1 D * / D 1 D * 系统下可能交换的玻色子，求出相关玻色子的通用交换势。 重味粒子的OBEP模型机制 2 OBEP模型 [ 11 11 ] 由Yukawa提出的通过单 π 交换势分析核力的理论 [ 12 12 ] 基础上发展而来。Yukawa认为质子和中子通过交换 π 介子结合在一起；Tornqvist [ 13 13 , 14 14 ] 也曾用单 π 交换势解释了氘核束缚态。同样地，在OBEP模型框架下，强子分子态被看作通过交换多种玻色子而形成的束缚态，利用如下方式可得到分子态中介子之间复杂的相互作用:通过考虑 π 介子的交换可得到长程力，通过考虑 σ 、 ρ 、 ω 、 η 介子的交换将得到中短程力 [ 11 11 , 13 13 , 14 14 ] 。但是强子分子态包含4个夸克，单纯求解4个夸克间玻色子的交换作用会使问题变得非常复杂，所以有必要在OBEP模型基础上进行假设和简化。考虑到重味分子态的系统由2个重夸克和2个轻夸克组成，假设各种玻色子的交换主要在轻夸克之间，而重夸克因过重而忽略玻色子交换的影响。因此，自旋算符和同位旋算符等只作用在轻夸克间。下文中为简单起见，用 A ( Q 1 q 2 )、 B ( Q 3 q 4 )分别表示2个重味介子。其中 q 2 、 q 4 分别表示2个轻夸克， Q 1 、 Q 3 分别表示2个重夸克， Q 1 和 Q 3 是 Q 1 、 Q 3 的反粒子。重味介子 D 1 和 D * 间的玻色子交换如 图 1 图 1 。 OBEP的推导 2 分子系统的量子数限制了能够发生交换的介子种类，下面将对 Y (4390)进行分析并推导各介子交换势。 Y (4390)的量子数是 J PC =1 -- ，对于量子数为1 -- 的 D 1 D * / D 1 D * 系统，它可能的轨道角动量和宇称是0 - 和2 - ，即该分子态的轨道可能处于 S 波或者 D 波； D 1 和 D * 的宇称分别为+1和-1，则系统的宇称自动满足量子数为-1的要求。在OBEP模型框架下，根据费曼图顶点处的宇称和角动量守恒，可以知道所有可能的交换玻色子 [ 6 6 ] ，相应的费曼图如 图 2 图 2 。 图 2 图 2 给出了 D 1 和 D * 散射(及其反粒子散射)之间交换玻色子所有可能的费曼图，并考虑了可能交换的中间玻色子为 σ 、 ρ 、 ω 、 π 、 η 。利用费曼规则 [ 14 14 ] 可分别处理赝标介子( π 、 η )、标量介子( σ )、矢量介子( ρ 、 ω )的交换贡献，从而得到它们在动量空间的 T 矩阵元。由于在玻恩近似下 T 矩阵等价于相互作用的交换势，因此可得到赝标介子、标量介子、矢量介子的交换势 [ 15 15 ] 。 下面以标量介子的交换势为例简要说明推导过程。标量介子的顶点耦合常数是 g sqq [ 16 16 ] ( g sqq 代表标量介子-夸克耦合常数，其中 q 表示夸克， s 表示标量介子)，用 p '和 p ″分别表示夸克的入射和出射四动量。交换的标量介子质量为 m ，四动量为 k ，对应的费曼图如 图 3 图 3 所示。根据 图 3 图 3 ，由费曼规则写出的动量空间的 T 矩阵元为: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465291&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465291&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465291&type=middle 其中， λ ' 2 、 λ ″ 2 分别表示入射和出射轻夸克 q 2 的螺旋度， λ ' 4 、 λ ″ 4 分别表示入射和出射轻夸克 q 4 的螺旋度， u λ ( p )、 u λ ( p )分别表示夸克在定螺旋度λ下对应的狄拉克旋量以及反旋量。 考虑非相对论性约化，将狄拉克旋量近似为: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465295&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465295&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465295&type=middle 其中， mq 表示夸克质量， σ 是泡利自旋矩阵， p '是夸克动量矢量，| i 〉表示夸克 i 的波函数。 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465299&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465299&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465299&type=middle 对于 q 4 也有同样的表达式，再利用介子的动量矢量 k = p ″- p '和( σ · a )( σ · b )= a · b +i σ ·( a × b )，对(2)式加以约化，最后可得动量空间下的有效势: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465303&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465303&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465303&type=middle 其中， http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465338&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465338&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465338&type=middle k 0 2 ， k 0 是介子的静止质量。计算(4)式时，仅保留到 p 的一次项和 k 的平方项。 为得到位置空间的有效势需要对(4)式进行傅里叶变换。但当存在交换势如 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465341&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465341&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465341&type=middle 时，傅里叶变换将会出现 δ 函数项: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465307&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465307&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465307&type=middle 其中， r 表示球坐标系中的半径。 δ 函数项出现的原因是直接把强子当作点粒子处理，而没有考虑强子的内部结构。因此需要在费曼图的顶点函数处引入反映强子内部结构的形状因子，其形式为 [ 16 16 ~ 18 18 16~18 ] : http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465311&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465311&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465311&type=middle 其中， X 2 =Λ 2 - m 2 ，Λ是截断参数，注意 k 和 k 分别表示介子的四动量和动量矢量。形状因子考虑了费曼图的高阶贡献，其中涉及到圈图动量的积分，积分动量在趋于无穷大时往往造成积分的发散，所以一般的做法是在某个大动量Λ处对积分进行截断，以压低大动量转移的贡献。大动量意味着短距离，Λ的值和强子的大小(近似为1/Λ)直接相关，在核子-核子相互作用中通常取在0.8~1.5 GeV范围内。由于本文讨论的是重味介子，其半径比核子更小，Λ的取值应该更大一些。同时当截断参数取值变大时，意味着考虑更多的动量交换，分子态间的重味介子吸引力更强，距离更小。当中间态足以形成分子态时，由3.2节可知，分子态的束缚能的确随着截断参数的增加而缓慢增加，均方根半径也随之变小。 本文用修正的顶点函数等效地描述复杂的胶子-夸克过程，在赝标玻色子和矢量玻色子中都引入形状因子的贡献，树图的顶点函数分别为i g pqq γ 5 和 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465343&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465343&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465343&type=middle [ 16 16 ] ，其中 g pqq 和 g vqq 是赝标与夸克的耦合参数， f vqq 是矢量玻色子与夸克的耦合参数， http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465347&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465347&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465347&type=middle 和 γ 5 均是狄拉克 γ 矩阵。类似(1)式~(5)式计算步骤，省略具体的计算过程，位置空间中各类型的玻色子交换对应的有效势可由 表 1 表 1 所示数据计算得到:将各类型粒子所在行的各项乘上表头的因子(如1、 σ 2 · σ 4 等)，然后相加，再乘上总因子。如赝标玻色子在 μ 2 ＞0下的交换势为 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465349&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465349&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465349&type=middle 。 表 1 表 1 后具体给出了 H 0 、 H 1 、 H 2 、 H 3 、 G 1 和 G 3 的表达式。 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465314&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465314&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465314&type=middle http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465316&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465316&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465316&type=middle 其中， http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465351&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465351&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465351&type=middle 是张量算符， http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465352&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465352&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465352&type=middle 是单位矢量， L 是轨道角动量算符。另外，对于同位旋 I =1的矢量介子，还需要对每一个顶点函数乘上一个同位旋算符 τ [ 15 15 ] 。因此，对于介子 π 和 ρ 的交换势， 表 1 表 1 中的每一项都需乘以因子 τ 2 · τ 4 。 OBEP的参数选取 1 本文采用非相对论薛定谔方程对强子分子态内部的重味介子相互作用进行近似。求解薛定谔方程的关键为势函数的选取，本文将1.2节得到的位置空间OBEP作为势函数的输入，其中OBEP的基本参数输入包括夸克和介子的质量和各耦合常数(如 表 1 表 1 )，粒子的质量由实验得到，可采用粒子数据表(particle data group，PDG)所给出的数据 [ 19 19 ] ，而 表 1 表 1 中的 f vqq 和 g pqq ， g sqq ， g vqq 是各个类型的玻色子-夸克耦合常数，目前没有确切的数值，但人们已通过大量实验和唯象模型分析，得到玻色子-核子耦合常数的数值，通过Goldberger-Treiman关系 [ 18 18 , 20 20 ] 可以推导出这两类耦合常数之间的关系: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465319&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465319&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465319&type=middle 其中， m N 是核子的质量， f ωNN 、 f ρNN 、 g πNN 、 g ηNN 、 g ρNN 、 g ωNN 、 g σNN 均是玻色子-核子耦合常数， f ωqq 、 f ρqq 、 g πqq 、 g ηqq 、 g ρqq 、 g ωqq 、 g σqq 是玻色子-夸克耦合参数。等号右边的耦合常数均是已知的。 Y (4390)在 D 1 D * / D 1 D * 系统下的势函数计算 1 本节将考虑 Y (4390)的结构，写出其系统波函数；并将1.2节得到的一般OBEP势应用在给定的基矢上，计算势函数的矩阵元，将该势函数代入薛定谔方程中，数值求解该方程。 OBEP模型机制下对 Y (4390)的结构分析 2 对于 J PC =1 -- 的 Y (4390)粒子，若作为 D 1 D * / D 1 D * 的强子分子态来处理，则它可能处于 S 波或 D 波轨道态，同时考虑到组分不带电的 D 1 0 D *0 和组分带正负电的 D 1 + D *- 的阈值不同，存在着同位旋破坏，并考虑到电荷共轭变换守恒，可以将 D 1 D * / D 1 D * 系统的波函数写为如下形式 [ 21 21 , 22 22 ] : http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465320&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465320&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465320&type=middle 其中，下标 S 和 D 分别表示 S 波和 D 波， u 1 ( r )、 u 2 ( r )、 u 3 ( r )、 u 4 ( r )表示对应的空间波函数，相应的基矢为 [ 20 20 ] : http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465323&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465323&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465323&type=middle 因为目前实验还没有确定 Y (4390)的同位旋量子数，所以在构造(9)式时还考虑到同位旋对称性破缺的可能，如果轨道波函数 u 1 和 u 3 ， u 2 和 u 4 的符号都相同，则系统同位旋 I =0主导；如果符号都相反，则同位旋 I =1主导。 总的相互作用势 V 应该是 表 1 表 1 中交换粒子势的总和，同时考虑到介子的G宇称 [ 20 20 ] ，有: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465324&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465324&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465324&type=middle 由此求出相互作用势的矩阵元 V Y (4390) [ 20 20 ] ，其中 V C 表示势能的常数项， V S 和 V SI 分别表示与同位旋无关和相关的由自旋算符 σ 2 · σ 4 引起的自旋-自旋相互作用项， V I 是同位旋作用项， V T 和 V TI 分别是与同位旋无关和相关的张量力项， V LS 和 V LSI 分别表示与同位旋无关和相关的自旋-轨道相互作用。 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465327&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465327&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465327&type=middle 在(11)式，位置空间各个有效势 V 分别是: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465329&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465329&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465329&type=middle http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465331&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465331&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465331&type=middle 值得说明的是， μ 2 = m 2 - k 0 2 中的 m 、 k 0 作为输入参数，其中 k 0 2 对应着初态和末态的夸克的静止能差，比如在(11)式的4×4方阵左上角的2×2方块阵中， k 0 2 =( m D * 0 - m D 1 0 )2，而其余的方块阵由于同位旋的破缺， D *+ D 1 - 会与 D * D 1 存在质量差， k 0 2 的取值会有所不同。 数值结果和讨论 2 本文将所得到的总的交换玻色子的有效势((11)式)，利用FORTRAN程序包FESSDE 2.2，计算四耦合道的径向薛定谔方程 [ 23 23 , 24 24 ] 的数值解。FESSDE 2.2利用有限元法可以精确求解二阶偏微分方程的本征值问题，将每一个Λ作为程序的参数输入，通过数值求解将得到对应的本征矢和本征能量，即系统的总能量 M 。二阶偏微分方程的本征解有很多，但由于分子态的势能是负的，本文在 D 1 D * 的阈值(4 430 MeV)之下寻找方程的本征值；同时考虑到截断参数应大于最大的中间玻色子质量( m ω =782 MeV)，以及文献[ 8 8 ]中氘核的Λ取800 MeV，对于重味强子态的Λ应该更大一些，所以本文从850 MeV开始取值，以大约4 MeV的步长，绘制束缚态质量 M 随Λ变化的趋势( 图 4 图 4 )，得到的主要数值解如 表 2 表 2 所示。 A/MeV M /MeV r rms /fm http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465354&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465354&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465354&type=middle 900 4 417.87 1.02 51.0:1.1:47.2:0.7 920 4 413.80 0.88 49.1:1.3:49.1:0.5 950 4 405.74 0.77 49.2:1.1:49.2:0.5 980 4 396.65 0.70 49.7:0.9:48.8:0.6 990 4 392.46 1.41 36.4:18.6:26.8:18.2 992 4 390.45 1.44 32.7:21.0:26.1:20.2 图 4 图 4 绘制了Λ取值为850~1 020 MeV时，在阈值4 430 MeV之下方程的所有根。首先，可以看到， M 随Λ的增加基本呈现平缓下降的变化趋势，符合1.2节的预期；其次，在Λ=850 MeV时，方程只有一个根，而在920 MeV左右方程有两个根。 图 4 图 4 左下的那些根对应最低能级，右上的那些根对应第一能级，而两个能级都可能存在束缚态。 对于 图 4 图 4 左下角的最低能级区域，在Λ=920 MeV附近寻找，即 M ~4 346.7 MeV，可以对应于BaBar合作组 [ 25 25 ] 发现的质量为4 346 MeV的粒子 Y (4360)，而实验上已经确认量子数为 J PC =1 -- ，符合本文在3.1节对基矢的构造要求。文献[ 26 26 ]也曾用单 π 交换势模型，将 Y (4360)解释为 D 1 D * 系统的分子态。 为了找到质量符合 Y (4390)的根，本文在 表 2 表 2 中给出了 图 4 图 4 右上角的第一能级 M 与Λ，以及对应的4个本征轨道波函数的相对占比，并由波函数直接计算出均方根半径(与波函数的平均半径有关)。总体而言，由 表 2 表 2 可知，随Λ的增大，均方根半径变小，符合1.2节的预期。然而Λ取值为990 MeV与980 MeV时，均可以得到质量接近 Y (4390)的两个束缚态，但是均方根半径出现很大的不同，Λ为980 MeV时，对应的均方根半径 r rms (0.70 fm)比Λ为990 MeV时对应的数值(1.41 fm)更小。Λ为980 MeV时，不同的组分比例中， S 波占主要成分，组分|1〉的 S 波和组分|3〉的 S 波占据了波函数的98.5%，轨道波函数处于基态，分子态的结合将更稳定更紧密；而Λ为990 MeV的态的均方根半径更大，分子态结合得较为松弛，且有 S - D 波的大成分混合，即组分|1〉的 S 波和组分|2〉的 D 波比例为36.4%:18.6%，组分|3〉的 S 波和组分|4〉的 D 波比例为26.8%:18.2%，即它处于激发态。为了更详细地进行观察， 图 5 图 5 给出了当截断参数取980 MeV和990 MeV时，系统各组分的空间波函数的分布。 图 5 图 5 中 u ( r )表示与4个基矢相关联的轨道波函数， r 表示轨道半径； 3 S 1 、 3 D 1 表示 S 波和 D 波函数在总角动量 J =1下的轨道态。 由 图 5 图 5 可知，当Λ为980 MeV时系统的态主要处于基态， 3 S 1 的轨道波函数占主导，且重味介子结合紧密，表现得更稳定，更可能形成强子分子态，符合本文最初的预期；同时结合(8)式可知， 3 S 1 的轨道波函数 u 1 和 u 3 ，以及 3 D 1 的轨道波函数 u 2 和 u 4 具有相反的符号，且 u 1 和 u 3 的占比相近，说明系统处于 I =1的同位旋多重态中。而Λ为990 MeV时系统的态不是一个基态，推测是 3 S 1 和 3 D 1 混合构成的激发态，激发态的粒子在非相对论中往往容易衰变，导致态不够稳定；结合 表 2 表 2 可知，Λ为990 MeV时，均方根半径较大，该束缚态束缚较松；同时 u 1 和 u 3 、 u 2 和 u 4 的符号总体相同，对应系统的 I =0部分占主导。综上考虑，本文认为对应截断参数为980 MeV， D 1 D * / D 1 D * 系统形成的质量约4 396.65 MeV的稳定分子束缚态，将更适合作为 Y (4390)的解释。 结论 1 本文采用单玻色子交换势模型从非相对论的薛定谔方程近似的角度研究了 J PC =1 -- 的 D 1 D * / D 1 D * 系统，并数值求解薛定谔方程，得到系统的束缚态质量和波函数，发现在截断参数Λ为980 MeV时，系统的态质量为4 396.65 MeV，处于轨道 3 S 1 态占主导的基态，均方根半径约为0.70 fm，重味介子结合紧密，构成了一个稳定的束缚态，本文希望将这个态用于解释新粒子 Y (4390)，比文献[ 13 13 ]使用Bethe-Salpeter方程近似得到的质量数值(4 384 MeV)更接近实验峰值。此外，本文还预言该粒子的系统中各组分的轨道波函数所占比等信息，并认为 Y (4390)很可能属于同位旋 I =1的多重态，是 Y (4360)在更高能级下的表现。但是目前没有测量到 Y (4390)的同位旋，这些预言有待实验数据的进一步确认。最后，本文得出结论: Y (4390)很可能是量子数为1 -- 的 D 1 D * / D 1 D * 的强子分子态。 参考文献 ABLIKIM M, ACHASOV M N, AHMED S N A, et al . Evidence of two resonant structures in e + e - → π + π - h c [J]. Physical Review Letters , 2017, 118: 092002-1-092002-8. DOI:10.1103/PhysRevLett.118.092002. MIZUK R, CHISTOV R, ADACHI I, et al . Observation of two resonance-like structures in the π + χ c 1 mass distribution in exclusive B 0 → K - π + χ c 1 decays[J]. Physical Review D , 2008, 78: 072004-1-072004-14. DOI: 10.1103/PhysRevD.78.072004. BONDAR A, GARMASH A, MIVUK R, et al . Observation of two charged bottomonium-like resonances in Y (5 S ) decays[J]. Physical Review Letters , 2012, 108: 122001-1-122001-6. DOI: 10.1103/PhysRevLett.108. 122001. GARMASH A, BONDAR A, KUZMI A, et al . Amplitude analysis of e + e - → ϒ( nS ) π + π - at http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465358&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465358&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4465358&type=middle =10.865 GeV[J]. Physical Review D , 2015, 91: 072003-1-072003-16. DOI: 10.1103/PhysRevD.91.072003. IOFFE B L. Calculation of baryon masses in quantum chromo-dynamics[J]. Nuclear Physics B , 1981, 188: 317-341. DOI:10.1016/0550-3213 (81) 90259-5. DING G J. Are Y (4260)and Z 2 + (4250) D 1 D or D 0 D * hadronic molecules?[J]. Physical Review D , 2009, 79: 014001-1-014001-17. DOI: 10.1103/PhysRevD.79.014001. JAFFE R L. Color non-singlet spectroscopy[J]. Physical Review D , 2005, 72:074508-1-074508-12. DOI:10. 1103/PhysRevD.72.074508. WANG Z G. Analysis of the Y (4220)and Y (4390) as molecular states with QCD sum rules[J]. Chinese Physics C , 2017, 41(8):083103-1-083103-8. DOI:10.1088/1674-1137/41/8/083103. WANG Z G. Vector tetraquark state candidates: Y (4260/4220), Y (4360/4320), Y (4390)and Y (4660/4630)[J]. The European Physical Journal C , 2018, 78 (6): 518. DOI: 10.1140/epjc/s10052-018-5996-5. HE J, CHEN D Y. Interpretation of Y (4390)as an isoscalar partner of Z (4430)from D * (2010) D 1 (2420)interaction[J]. The European Physical Journal C , 2017, 77 (6): 398. DOI: 10.1140/epjc/s10052-017-4973-8. THOMAS C E, CLOSE F E. Is X(3872)a molecule? [J]. Physical Review D , 2008, 78: 034007-1-034007-12. DOI:10.1103/PhysRevD.78.034007. YUKAWA H. On the interaction of elementary particles Ⅰ[J]. Supplement of the Progress of Theoretical Physics , 1935, 17: 48-57. DOI: 10.1143/PTPS.1.1. TORNQVIST N A. Possible large deuteron-like meson meson states bound by pions[J]. Physical Review Letters , 1991, 67(5): 556-559. DOI:10.1103/PhysRevLett. 67.556. TORNQVIST N A. From the deuteron to deusons, an analysis of deuteron-like meson-meson bound states[J]. Zeitschrift für Physik C Particles and Fields , 1994, 61 (3): 525-537. DOI: 10.1007/BF01413192. SCADRON M D. Advanced Quantum Theory and Its Applications Through Feynmann Diagrams [M]. New York: Springer-Verlag, 1979: 228. MACHLEIDT R, HOLINDE K, ELSTER C. The bonn meson exchange model for the nucleon nucleon interaction[J]. Physics Reports , 1987, 149(1): 1-89. DOI: 10. 1016/S0370-1573(87)80002-9. MACHLEIDT R. The meson theory of nuclear forces and nuclear structure[J]. Advances in Nuclear Physics , 1989, 19: 189-376. 宁平治.原子核物理基础:核子与核[M].北京:高等教育出版社, 2003: 224. NING P Z. Fundamental Nuclear Physics : Nucleons and Nuclei [M]. Beijing: Higher Education Press, 2003: 224 (Ch). PATRIGNANI C, AGASHE K, AIELLI G, et al . Review of particle physics[J]. Chinese Physics C , 2016, 40 (10):100001-1-100001-5. DOI:10.1088/1674-1137/40/10/100001. DING G J, LIU J F, YAN M L. Dynamics of hadronic molecule in one-boson exchange approach and possible heavy flavor[J]. Physical Review D , 2009, 79:054005-1-054005-18. DOI: 10.1103/PhysRevD.79.054005. LI M T, WANG W L, DONG Y B. A study of one S -and one P -wave heavy meson interaction in a chiral quark model [J]. Communications in Theoretical Physics , 2015, 63 (1): 63-70. DOI:10.1088/0253-6102/63/1/11. LI M T, WANG W L, DONG Y B, et al . Z b (10650) and Z b (10610)states in a chiral quark model[J]. Journal of Physics G : Nuclear and Particle Physics , 2013, 40(1): 1-10. DOI:10.1088/0954-3899/40/1/015003. ABRASHKEVICH A G, ABRASHKEVICH D G, KASCHIEV M S, et al . Finite-element solution of the coupled-channel Schrodinger equation using high-order accuracy approximations[J]. Computer Physics Communications , 1995, 85: 40-64. ABRASHKEVICH A G, ABRASHKEVICH D G, KASCHIEV M S, et al . FESSDE, a program for the finite-element solution of the coupled-channel Schrodinger equation using high-order accuracy approximations[J]. Computer Physics Communications , 1995, 85 (1): 65-81. AUBERT ET B, BARATE R, BONA M, et al . Evidence of a broad structureat an invariant mass of 4.32 GeV/ c 2 in the reaction e + e - → π + π - ψ (2 S ) measured at BaBar[J]. Physical Review Letter , 2007, 98: 212001-1-212001-7. DOI: 10.1103/PhysRevLet.98.212001. CLOSE F, DOWNUM C, THOMAS C E. Novel charmonium and bottomonium spectroscopies due to deeply bound hadronic molecules from single pion exchange[J]. Physical Review D , 2010, 81:074033-1-074033-12. DOI: 10.1103/PhysRevD.81.074033. 10.14188/j.1671-8836.2019.05.011 O571.21 whdxxblxb-65-5-495 中央高校基本科研业务费项目资助(CCNU18TS029) 2018-04-08 2019-10-24 2020-12-10 版权所有© 2019 《武汉大学学报(理学版)》编辑部 Copyright © 2019 Journal of Wuhan University(Natural Science Edition). All rights reserved. 2019 武汉大学学报(理学版) 5 65