荔 寧,高探彪,吳亞杰
(西安工業(yè)大學(xué) 理學(xué)院,西安 710021)
自然界中存在四種相互作用,強(qiáng)相互作用、電磁相互作用、弱相互作用以及引力相互作用.在這四種相互作用中強(qiáng)相互作用最強(qiáng),正是由于有強(qiáng)相互作用的存在,才能使夸克、膠子形成介子、重子等,從而形成世界萬(wàn)物.而研究強(qiáng)相互作用最重要的理論就是量子色動(dòng)力學(xué).因此,強(qiáng)子譜的研究在量子色動(dòng)力學(xué)中是最基本也是最重要的部分.在強(qiáng)子譜的研究過(guò)程中會(huì)碰到非微擾的問(wèn)題,這促使人們尋找一種非微擾的手段來(lái)研究該問(wèn)題.而格點(diǎn)量子色動(dòng)力學(xué)是從量子色動(dòng)力學(xué)的第一性原理出發(fā)來(lái)進(jìn)行數(shù)值計(jì)算的一種方法.通過(guò)在格點(diǎn)上計(jì)算各種關(guān)聯(lián)函數(shù),就可以得到目標(biāo)物理量.由于格點(diǎn)上的計(jì)算是在有限的格子體積上,而真實(shí)世界的物理量是連續(xù)空間的,因此,對(duì)于格點(diǎn)上得到的物理量需要進(jìn)行連續(xù)極限的外推.
近年來(lái),BaBar,BELLE,BES是含粲夸克態(tài)[1].例如,2002年,BaBar研究組發(fā)現(xiàn)了一些宇稱(chēng)為正的標(biāo)量介子DsJ(2317).同時(shí),CLEO觀察到了1+的態(tài)DsJ(2460).這兩個(gè)共振態(tài)的衰變寬度都很窄,并且在DK以及D*K的閾值以下[2-5].
格點(diǎn)上要得到某個(gè)物理量的值,必須先計(jì)算與該物理量相關(guān)聯(lián)的關(guān)聯(lián)函數(shù),由于有格距的存在且采用周期邊界條件,使得關(guān)聯(lián)函數(shù)所對(duì)應(yīng)的動(dòng)量的最小間隔是2π/L,要減小這個(gè)動(dòng)量間隔,從而獲得更多動(dòng)量的關(guān)聯(lián)函數(shù),需要引入扭轉(zhuǎn)的邊界條件,文中討論了格點(diǎn)上的關(guān)聯(lián)函數(shù)以及扭轉(zhuǎn)邊界條件,為在格點(diǎn)上計(jì)算粒子的質(zhì)量和能量提供了一定的理論基礎(chǔ)[7].
文中主要計(jì)算介子的質(zhì)量以及能量,因此,需要在格點(diǎn)上計(jì)算兩點(diǎn)關(guān)聯(lián)函數(shù),通過(guò)計(jì)算兩點(diǎn)關(guān)聯(lián)函數(shù),從而提取出該粒子的質(zhì)量以及能量.在計(jì)算關(guān)聯(lián)函數(shù)之前需要在格點(diǎn)上構(gòu)造能在真空中產(chǎn)生與該粒子JPC相同的內(nèi)插場(chǎng)算符.文中采用局域的夸克與反夸克介子算符,其算符形式為Ω.這時(shí),歐氏空間的關(guān)聯(lián)函數(shù)為
(1)
(2)
其中n為能級(jí)級(jí)數(shù)。
(3)
(4)
(5)
(6)
(7)
(8)
通過(guò)調(diào)節(jié)扭轉(zhuǎn)的角度,可得到任意的動(dòng)量.
文中的數(shù)值模擬細(xì)節(jié)分為兩部分:組態(tài)及算符構(gòu)造和質(zhì)量譜及色散關(guān)系.在組態(tài)及算符構(gòu)造中主要論述組態(tài)的具體信息以及如何構(gòu)造具有確定量子數(shù)的算符,并將物理基下的算符通過(guò)一定的關(guān)系
式轉(zhuǎn)到使用于Twisted Mass組態(tài)中的扭轉(zhuǎn)基下面.而在質(zhì)量譜及色散關(guān)系中展示B1和C1兩套組態(tài)下的質(zhì)量平臺(tái)圖以及不同動(dòng)量模式下的能量圖,通過(guò)擬合不同的動(dòng)量所對(duì)應(yīng)的能量得到色散關(guān)系中的Z因子.
數(shù)值模擬基于Nf=2的Twisted Mass組態(tài),該組態(tài)由ETMC小組生成,在海夸克中包含了u夸克和d夸克是非淬火的組態(tài),且這個(gè)組態(tài)對(duì)應(yīng)的扭轉(zhuǎn)角度是π/2,即最大的扭轉(zhuǎn)角度.這個(gè)組態(tài)最大的優(yōu)勢(shì)在于其是格距次方改進(jìn)的,也就是說(shuō)在該組態(tài)上算出的物理量的格距誤差是格距的平方[6].具體的組態(tài)參數(shù)見(jiàn)表1[9].
表1 組態(tài)參數(shù)
為了在格點(diǎn)上進(jìn)行兩點(diǎn)關(guān)聯(lián)函數(shù)的計(jì)算,引入了三個(gè)扭轉(zhuǎn)doublets,(u,d),(s,s′),(c,c′)作為價(jià)夸克,質(zhì)量分別是μl,μs,μc:在格子243×48上的具體數(shù)值分別是0.004,0.018,0.215;在格子323×64上的具體數(shù)值分別是0.003,0.015,0.185.在每個(gè)doublets中,威爾遜參數(shù)具有相反的符號(hào)(r=-r′=1).進(jìn)行變換后可以將物理基下的夸克場(chǎng)變成扭轉(zhuǎn)基下的夸克場(chǎng)[10-14]:
(9)
(10)
2) 由于計(jì)算得到的兩點(diǎn)關(guān)聯(lián)函數(shù)在時(shí)間方向上是反周期邊界條件,且式(10)關(guān)于T/2是對(duì)稱(chēng)的.因此在處理數(shù)據(jù)的時(shí)候可以將兩點(diǎn)關(guān)聯(lián)函數(shù)關(guān)于T/2作對(duì)稱(chēng)化處理,這樣相當(dāng)于間接地增加了兩點(diǎn)關(guān)聯(lián)函數(shù)的統(tǒng)計(jì)量.
圖1 組態(tài)B1上介子的質(zhì)量平臺(tái)圖
圖2 組態(tài)C1上介子的質(zhì)量平臺(tái)圖
圖3 組態(tài)B1上介子r(t)的平臺(tái)圖
圖4 組態(tài)C1上介子r(t)的平臺(tái)圖
表2 組態(tài)B1,C1以及格距外推結(jié)果和PDG上介子質(zhì)量的對(duì)比
利用格點(diǎn)上所得到的能量來(lái)擬合格點(diǎn)上的色散關(guān)系[10],即
(11)
(12)
表3 組態(tài)B1上不同動(dòng)量所對(duì)應(yīng)的介子的能量
表4 格點(diǎn)上擬合的Z因子的值以及理論值
表4中列出了利用表3中格點(diǎn)計(jì)算出來(lái)的能量通過(guò)式(11)在兩套組態(tài)上擬合得到的Z因子的值,擬合后得到的Z因子與理論值1非常接近.