• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于量子計算的高能核物理研究*

    2023-11-16 10:43:20李天胤邢宏喜4張旦波
    物理學報 2023年20期
    關鍵詞:核物理夸克強子

    李天胤 邢宏喜4)? 張旦波

    1)(華南師范大學量子物質研究院,原子亞原子結構與量子調控教育部重點實驗室,廣州 510006)

    2)(華南師范大學量子物質研究院,廣東省核物質科學與技術重點實驗室,廣州 510006)

    3)(華南師范大學南方核科學計算中心,粵港量子物質聯(lián)合實驗室,廣州 510006)

    4)(華南師范大學物理學院,原子亞原子結構與量子調控教育部重點實驗室,廣東省高等學校物質結構與相互作用基礎研究卓越中心,廣州 510006)

    高能核物理旨在探索和理解物質在夸克與膠子層次的組成及演化的基本規(guī)律.然而,從量子色動力學第一性原理出發(fā)來求解高能核物理,經典計算卻存在本質困難.近幾年來,量子計算因給模擬高能核物理提供了潛在的根本性解決方案而受到了較大的關注.本文簡要回顧了高能核物理量子模擬的現狀,介紹了態(tài)制備及光錐關聯(lián)函數測量等典型量子算法,并通過強子散射振幅和有限溫有限密物質相結構的研究,分別展示了量子計算在解決高能核物理中含時問題和符號問題上的優(yōu)勢.

    1 引言

    高能核物理研究夸克與膠子及其強相互作用形成的物質及其演化,對理解物質深層次結構乃至宇宙起源均具有重要意義.高能核物理的底層模型量子色動力學(QCD)具有色禁閉的關鍵特征,即在低能標時QCD 處于強耦合區(qū)域,夸克和膠子等被束縛在色中性的強子中,不能單獨存在.由于夸克禁閉的存在,高能核物理中的許多重要物理量,如強子散射振幅,具有非微擾特性,導致由第一性原理的QCD 理論出發(fā)來理解高能核物理存在基本困難.另一方面,理解核物質在不同溫度及密度下的物態(tài)特性,即核物質相圖,是高能核物理的另一個重大挑戰(zhàn).這是由于在有限密度情形,格點QCD 的蒙特卡羅模擬存在符號問題,從而無法保證精度.總而言之,經典計算機在模擬高能核物理中的含時演化問題[1]及有限密度核物質上[2],均存在本質困難.

    高能核物理經典模擬存在困難的本質在于它屬于強關聯(lián)量子多體系統(tǒng).一個根本性的解決方案是通過可控量子系統(tǒng)來模擬QCD 及相關的高能核物理,此即量子計算模擬方案.量子計算機最先由費曼在1982 年提出[3],Jordan,Lee 和Preskill(JLP)[4]最早在2012 年證明了量子計算機以多項式復雜度模擬標量場中的散射問題.近年來,隨著量子計算機硬件及量子算法的快速發(fā)展,量子計算與高能核物理的交叉研究逐漸成為研究熱點,其中包括核子結構中部分子分布函數的計算[5-10]、部分子簇射問題[11,12]、與環(huán)境有相互作用的開放系統(tǒng)動力學問題[13-15]、非平衡態(tài)熱力學問題[16,17]、高能核物理中散射振幅的計算[18,19]等含時問題[20,21].另一方面,量子計算也用于模擬高能核物理中有限溫有限密系統(tǒng),如有限溫有限密系統(tǒng)的手征相變的量子模擬[22-24]、手征磁效應的量子模擬[25]、以及如何在量子計算機上通過變分量子算法來模擬高溫時夸克的禁閉-退禁閉行為[26].

    本文主要對我們近期基于量子計算模擬高能核物理的工作進行綜述[9,10,19,26].第2 節(jié)介紹了在高能核物理中常用的量子算法;第3 節(jié)介紹了如何基于第2 節(jié)所提的量子算法計算強子散射振幅;第4 節(jié)則介紹如何通過量子計算機來研究有限溫有限密核物質;最后第5 節(jié)是對這一系列工作的總結.

    2 高能核物理相關量子算法

    在量子計算機上模擬高能核物理一般包含兩部分:制備描述強子或核物質的量子態(tài),以及在該態(tài)上進行物理量的測量.本節(jié)基于量子線路模型的數字量子計算機,介紹相關的典型量子算法.

    2.1 強子態(tài)制備

    如何在量子計算機上模擬一個強子?最核心的步驟是制備描述強子的量子態(tài).強子態(tài)是給定量子數下的激發(fā)態(tài).從多體波函數的角度出發(fā),強子態(tài)一般極為復雜,包含夸克及膠子各種構型的疊加,具有較復雜的量子關聯(lián).對于NISQ(noisy intermediate-scale quantum era)量子計算機,變分量子本征值求解器[27-30]是制備強子態(tài)的適用方法,通過采用合適的波函數擬設,可以滿足量子數的要求并精確表達強子態(tài).該算法的框架如圖1所示.首先,采用量子數保持的擬設,即試探波函數由保持系統(tǒng)對稱性的參數化幺正算符U(θ)作用到給定量子數的參考態(tài)|ψlk〉ref上,

    圖1 變分法制備強子態(tài)的量子線路圖.其中量子計算機通過擬設U(θ)來生成試探波函數以及測量試探波函數下El(θ)的值.經典計算機負責更新和優(yōu)化參數θFig.1.Variational quantum algorithm for preparation of the hadronic state.The generation of trial state and the measurement of El(θ) are performed on a quantum computer,while the optimization of parameters θ is done with classical computing.

    其中l(wèi)代表系統(tǒng)除能量外的所有量子數,k代表給定量子數l下的第k個激發(fā)態(tài).參考態(tài)|ψlk〉ref為類Dicke態(tài),其一般的制備方法可在文獻[31]中找到.

    保持系統(tǒng)對稱性的幺正算符U(θ)可以通過交替算符演化擬設[32,33]來構造.為此,先把哈密頓量分成n(≥2)份,HH1+···+Hn,并滿足以下要求:1)每一個Hi都與哈密頓量H具有一樣的對稱性;2) [Hi,Hi+1]0 對所有i成立;3)選取盡量大的n.U(θ)可以通過多層n個哈密頓量H1,···,Hn的交替演化來構造:

    其中p是演化的層數.

    制備強子態(tài)的第二步是最優(yōu)化目標函數El(θ)[28],

    其中wli滿足wl1>wl2>···>wlk >0.在最優(yōu)參數θ*得到之后,量子數為l的第k激發(fā)態(tài)強子可通過|h〉U(θ*)|ψlk〉ref來制備.值得強調的是,El(θ)通過量子-經典混合計算的方式來優(yōu)化(如圖1 所示),其中量子計算機制備試探態(tài)并進行測量,經典計算機收集測量數據并更新參數θ,迭代直至收斂.

    2.2 有限溫度系統(tǒng)的平衡態(tài)制備

    在QCD 相圖中,核物質在不同溫度、不同密度下呈現不同的物態(tài),特別地,在溫度非常高時,夸克出現解禁閉現象,與膠子構成了夸克-膠子等離子體,一直以來是高能核物理的研究重點.然而,傳統(tǒng)的格點色動力學計算方法在模擬有限密度核物質時存在符號問題.在量子計算上,可以直接制備有限溫度系統(tǒng)的平衡態(tài)[34-38],亦即熱態(tài)或吉布斯態(tài),該方法不基于概率來采樣,對零密度及有限密度均適用,從而避免了符號難題.與強子態(tài)制備類似,變分量子算法適用于在NISQ 量子計算機上模擬有限溫系統(tǒng)[36,39,40].然而,制備熱態(tài)需最小化自由能,過程中熵的測量對量子計算機而言是較大的挑戰(zhàn)[36,39].下面介紹一種變分量子算法,通過特定的擬設避免了熵的測量,同時對一般系統(tǒng)的熱態(tài)又有很好的表達力.

    該算法的關鍵在于試探態(tài)可由參數化的幺正算符U(?) 作用到初始混態(tài)生成

    其中ρ0(θ)為直積態(tài),ω(θ,?).由于U(?) 是幺正的,它不會改變熵.又由于初態(tài)具備的直積形式,系統(tǒng)的熵S(ω) 可表達為單比特混態(tài)ρi(θi)上的熵之和.由此,系統(tǒng)的變分自由能可寫為

    其中β1/T.在參數空間中最優(yōu)化自由能,即可得到最優(yōu)參數ω*,則吉布斯態(tài)為ρ(ω*).與強子態(tài)制備的變分量子算法類似,自由能的優(yōu)化通過量子-經典混合計算的方式進行,如圖2 所示.

    圖2 可同時制備吉布斯態(tài)和計算自由能的變分算法量子線路圖.在量子計算機中執(zhí)行的是初始混態(tài)的制備和參數化幺正演化 U(?),而參數(θ,?) 的優(yōu)化通過經典計算機來執(zhí)行Fig.2.A pictorial representation of the variational quantum algorithm,which can prepare thermal states and compute the corresponding free energy.The preparation of initial mixed state and the evolution U(?) are done on a quantum computer,while the optimization of(θ,?) should be done with a classical computer.

    2.3 關聯(lián)函數的測量

    在態(tài)制備之后,量子模擬的另一個重要環(huán)節(jié)是對態(tài)進行測量,得到描述該系統(tǒng)的物理量.物理量一般通過關聯(lián)函數來揭示,它包含靜態(tài)關聯(lián)函數及動力學關聯(lián)函數.在量子計算機上,靜態(tài)關聯(lián)函數只需對量子比特進行關聯(lián)測量即可.動力學關聯(lián)函數的測量更為復雜,具體算法可參考文獻[41].在核子結構的研究中,部分子分布函數、碎裂函數等非微擾量的計算依賴于一種特殊的兩點動態(tài)關聯(lián)函數,即光錐關聯(lián)函數,涉及了光錐上兩個算符的關聯(lián).此外,散射振幅也需要通過計算含時的多點關聯(lián)函數得到.動力學關聯(lián)函數涉及了哈密頓量演化,屬于典型的含時演化問題.相對經典計算機,量子計算機在求解動力學關聯(lián)函數上具備天然的優(yōu)勢.

    先考慮一般形式的兩點關聯(lián)函數

    圖3 計算動力學兩點關聯(lián)函數量子線路.在輔助比特上測量 σ1和σ2 即可得到關聯(lián)函數的實部和虛部Fig.3.Quantum circuit for calculation of the two point correlation function.The real and imaginary part of the correla tion function can be obtained by performing measurements of σ1 and σ2 on the auxiliary qubit,respectively.

    上述即為計算兩點動力學關聯(lián)函數的量子算法,對于多點動力學關聯(lián)函數,只需令輔助量子比特分別控制多點關聯(lián)函數對應的算符并在相鄰的兩算符間插入對應的時間演化算符即可.值得注意的是,哈密頓量演化 e-iHt需通過Trotter 分解的方法表示為量子門的組合[42],從而適用于在數字量子計算機上運行.

    3 強子散射振幅的量子計算

    在高能核物理中,散射截面是聯(lián)系理論和實驗的橋梁,散射振幅是計算散射截面不可或缺的一部分,因此具備重要的地位.強子散射過程是指初態(tài)有強子參與或者末態(tài)有強子產生的散射過程.根據末態(tài)探測粒子數目的不同,可對強子散射過程進行分類.1)e+p →e′+X,在末態(tài)只測量一個電子的信息而忽略其他X粒子,該類過程被稱為單舉過程;e+p →e′+h+X,在末態(tài)測量電子和一個強子h的信息,該類過程被稱為半單舉散射過程;e+γ →e′+h,把末態(tài)粒子信息全部測量完的被稱為遍舉過程.在QCD 因子化定理適用范圍下,初態(tài)有強子參與的單舉過程可被因子化成微擾可計算的硬部分卷積非微擾的部分子分布函數(parton distribution functions,PDFs[43,44]);末態(tài)只生成一個強子的遍舉過程散射振幅則可被因子化為硬部分卷積的非微擾量,即光錐分布振幅(light cone distribution amplitudes,LCDAs[45,46]).從物理含義上看,PDFs 反映了強子內部部分子的縱向動量分布,介子LCDAs 則描述了一對正反夸克在遍舉過程中如何強子化成介子.PDFs 與LCDAs 均為計算強子散射振幅必要的輸入.格點QCD 發(fā)展了許多方法來嘗試計算PDFs[47-62]和LCDAs[63-66].然而,由于PDFs 和LCDAs 均為含時演化問題,依賴于光錐關聯(lián)函數的求解,對格點QCD 而言存在經典計算的本質困難.此外,在碰撞能量較低時,QCD 散射過程可看作是強子間的彈性散射.由于QCD 耦合常數在低能標時較大,微擾論失效,散射振幅可通過格點QCD 間接求解[67-71],但缺乏一般性的直接求解方法.

    總的來說,運用格點QCD 的方法解決上述關于強子散射振幅的非微擾問題存在困難.本節(jié)將介紹如何基于量子計算的方法來求解高能散射緊密關聯(lián)的三個問題,包括單舉過程對應的PDFs、遍舉過程對應的LCDAs,以及低能情形且初末態(tài)強子數目較少過程的散射振幅.

    3.1 Nambu-Jona-Lasinio 模型

    由于近期的量子計算機計算能力有限,規(guī)范場的模擬仍較難在量子計算機上實現.為了從數值上驗證量子算法的有效性,我們選取1+1 維Nambu-Jona-Lasinio(NJL)模型[72-74]作為例子.NJL 模型是一個四費米子自相互作用模型,可視為QCD的低能有效理論,且具有手征對稱性自發(fā)破缺與漸近自由等和QCD 一樣的特點.其拉格朗日量為

    其中Ψ是費米子場,m和g為別為裸夸克質量以及NJL 模型耦合常數.

    在量子計算機上模擬NJL 模型的第一步是對其進行離散化.采用交錯費米子,離散化之后的費米子場可由單分量場χn表示.通過Jordan-Wigner[75]變換,可將費米子場映射成量子比特形式:

    該變換可保持費米子場的反對易關系.在該離散化機制下,兩個量子比特對應著一個空間格點,近期約100 量子比特的量子計算機可模擬約50 個空間格點.NJL 模型的量子模擬對量子比特數資源需求較小,有利于在近期量子計算機上開展研究.在周期邊界條件下,NJL 模型的哈密頓量可寫成HH1+H2+H3+H4,其中

    3.2 部分子分布函數的量子計算

    PDFs 是計算強子初態(tài)單舉過程散射截面的重要輸入量.PDFs 一般被記為fi/h(x,μ),其物理意義是在能標μ下,非極化部分子i占據非極化強子h的縱向動量份額為x的概率密度.夸克PDFs算符定義為

    Pexp 表示路徑順序指數,τa是蓋爾曼矩陣,是SU(3)規(guī)范場.從算符定義可以看出,PDFs 的量子計算包含了3 個關鍵步驟:1)制備強子態(tài)|h(P)〉;2)構造威爾遜鏈W(0) ;3)在量子態(tài)|h(P)〉上測量動力學兩點關聯(lián)函數.

    由于NJL 模型沒有規(guī)范場,我們用NJL 模型來具體討論如何通過步驟1)和步驟3)來實現PDFs 的計算,而對步驟2)的部分僅做一般性的討論.在強子參考系下,光錐關聯(lián)函數D(z)可寫為

    其中|h〉表示的是零動量強子態(tài)|h(p0)〉.PDFs可通過D(z)的傅里葉變換得到

    其中mh為強子的靜質量,x表示部分子的動量份額.

    D(z)的量子計算如圖4 所示,包含強子態(tài)制備和光錐關聯(lián)函數的測量,分別對應圖中虛線的左右部分的量子線路.具體來說,NJL 模型的強子態(tài)|h〉是能量最低的一對正反夸克束縛態(tài),它是與真空具有一樣量子數的第一激發(fā)態(tài).由此,制備|h〉態(tài)需要兩個參考態(tài):

    圖4 部分子分布計算的量子線路.其中虛線左邊部分為制備強子態(tài)的量子線路,右邊部分為測量動力學兩點關聯(lián)函數線路Fig.4.Quantum circuit for calculation of PDFs.The left side of the dashed line of the circuit is for hadronic state preparation,while the right side is for the correlation function.

    在給定參考態(tài)后,試探波函數由(2)式中U(θ)作用到參考態(tài)上生成,為此,只要將U(θ)中的Hi具體化為NJL 模型哈密頓量(10)且選取pN/2.在最優(yōu)化目標函數(3)后即可制備出真空態(tài)|Ω〉U(θ*)|ψΩ,1〉ref和強子態(tài)|h〉U(θ*)|ψΩ,2〉ref.成功制備強子態(tài)之后,只需把強子態(tài)作為圖4 虛線右側量子線路的輸入即可計算出D(z).

    圖4 所示量子線路的復雜度為O(N3)[33,42].量子算法復雜度隨著量子比特數以多項式標度增長說明了在量子計算機上計算PDFs 具有優(yōu)勢.更進一步的分析表明,量子優(yōu)勢在更高維度,如3+1 維且含規(guī)范場的情況下仍能保持[9].

    下面展示1+1 維NJL 模型夸克PDFs 的結果.我們使用Quspin[76]以及ProjectQ[77]在經典計算機模擬圖4 中的量子線路.為了控制有限體積效應和有限格距效應的影響,要求無量綱的強子質量滿足條件 2π/N <mha <π.同時,考慮到量子計算經典模擬的指數復雜度,選取格點數9(對應量子比特數N18)和無量綱質量ma0.8.圖5 分別給出了D(z)的實部和虛部在g0.6 和g1.0情況下的結果.從圖5 可以看出,D(z)的實部與虛部相比基本為零,這意味著fq(x) 滿足奇函數條件fq(x)-fq(-x).PDFs 在減掉真空圖后一定會有fqˉ(x)-fq(-x)[78],結合fq(x) 是奇函數,可得到fq(x)fqˉ(x),這與NJL 模型中強子態(tài)|h〉由一對正反夸克組成相自洽.

    圖5 坐標空間夸克場兩點關聯(lián)函數的實部(實線)和虛部(虛線),其中離散的點是格點計算給出的數值,線由格點數據結果插值得到Fig.5.Real(dashed lines) and imaginary parts(solid lines)of the quark correlation function in position space.The discrete points are the lattice data and the lines are obtained by interpolations.

    部分子分布函數的結果在圖6 中給出.由于格點有限體積效應的存在,fq(x) 在非物理區(qū)域x >1有小的非零值,這一現象同樣存在于格點QCD 的計算當中[79].可預期在增加格點數后,fq(x) 在區(qū)間x >1 處的非零值會進一步被壓低.我們在結果中觀察到了fq(x)在x0.5 處具有一個峰值,這一定性行為與二維格點QCD 的計算結果[80]以及JAM合作組從實驗數據提取出的介子PDF[81]相符.

    圖6 量子計算經典模擬獲得的PDF(空心點)和格點NJL模型精確對角化獲得的PDF(實線),其中不同插值方法帶來的誤差由誤差棒標記出Fig.6.The quark PDF from quantum computing(open markers) and ED(solid lines).The error bars/bands arise from the estimated uncertainties due to different interpolation methods.

    3.3 基于量子計算方法對LCDAs 的研究

    介子LCDAs 是計算遍舉過程散射振幅的重要非微擾輸入量,其物理意義可看作是一對分別攜帶強子縱向動量份額為x和1-x的正反夸克對結合成大動量強子的概率幅.介子LCDAs 的算符定義為

    在強子參考系下,1+1 維單味道NJL 模型的介子LCDA 可寫成

    圖7 模擬左矢為真空,右矢為強子的光錐關聯(lián)函數〈Ω|O|h〉的量子線路圖Fig.7.Quantum circuit for the light-cone correlator〈Ω|O|h〉.

    以下展示數值模擬的結果.如圖8 所示,實線(ED)代表的是對格點NJL 模型進行數值精確對角化的結果,空心點代表的是量子計算經典模擬的結果,二者的符合體現了量子算法的準確性.此外,LCDA在x1/2 處的峰隨著g →0 變得越來越尖銳且趨近于LCDA 的漸近形式[45,46,82].此外,由于有限體積效應,LCDA在x >1 以及x <0 非物 理區(qū)間有振蕩衰減的非零值.這一現象也存在于格點QCD 對LCDA 的計算當中[83].圖9 還給出了LCDA 對強子質量mh的依賴關系.圖中分別選擇mh1.3a-1,1.5a-1,1.7a-1和固定g0.1,N14,可見隨著強子質量的增大,LCDA在x1/2 處的峰變得越來越尖銳,這一定性行為與格點QCD 的結果相符[83].

    圖8 1+1 維NJL 模型LCDA 對耦合常數g 的依賴,其中固定 N 14,mh 1.5a-1Fig.8.LCDA for the 1+1 dimensional NJL model with N 14,mh 1.5a-1.

    圖9 1+1 維NJL模型LCDA 對強子質量mh 的依賴,其中固定 N 14,g 0.1Fig.9.Dependence of the LCDA on the hadron mass mh with fixed bare coupling g 0.1 in the 1+1 dimensional NJL model.

    3.4 LSZ 約化公式及散射振幅的量子計算

    Jordan,Lee 和Preskill(JLP)[4,84,85]的系列論文證明了量子計算機可通過多項式復雜度計算標量場高能過程散射截面.然而,JLP 方案需要制備初態(tài)波包,對計算資源要求較高.另外,JLP 方案里未涉及束縛態(tài)粒子,而QCD 的散射過程的初態(tài)往往會有強子態(tài)出現,因此難以直接使用.為此,基 于Lehmann-Symanzik-Zimmermann(LSZ)約化公式[86],我們提出量子計算束縛態(tài)粒子的散射振幅的一種新途徑.LSZ 公式把散射過程的散射振幅和真空多點關聯(lián)函數聯(lián)系在一起.對于質量為m,可被標量場?湮滅的粒子h參與的散射過程h(k1)+···+→h(p1)+···+,其散射截面可表達成

    其中NNin+Nout.G({pi},{kj}) 是動量空間的連通n點關聯(lián)函數:

    其中T表示的是時序,下標“con”表示關聯(lián)函數的連通部分,|Ω〉是真空態(tài).K(p) 是動量空間的兩點關聯(lián)函數,

    其中R是場強歸一化常數,它的定義為

    其中|h〉表示零動量的單個h粒子態(tài),它可以是強子態(tài).依據(19)式,散射振幅的計算被分為3 個部分:1)計算連通N點關聯(lián)函數G({pi},{kj}) ;2)計算傳播子K(p) ;3)計算場歸一化常數R.1)、2)兩個步驟都是在計算多點關聯(lián)函數,他們可統(tǒng)一通過2.3 節(jié)提到的量子算法來計算.步驟3)中R的計算類似于LCDAs,只不過算符替換成?(0),因此R可用LCDAs 的量子算法來計算.

    下面來討論這一量子算法的復雜度,假設待模擬的格點QFT 系統(tǒng)共有M個空間格點以及M個時間格點,表達一個空間格點上的場自由度需要nq個量子比特,那么在該系統(tǒng)計算N點關聯(lián)函數的復雜度為,它指數依賴于N.因此,該算法適用于初末態(tài)粒子數量較少的散射過程,如2→2散射過程.

    作為一個例子,我們依然以1+1 維NJL 模型這一簡單模型為例來具體給出量子算法的運行結果,并通過這些結果來演示量子算法的合理性.在下面的模擬中選取量子比特數N14,無量綱裸夸克質量為ma0.84,耦合常數為g0.4.

    在動量趨于h的質殼附近時,1+1 維單味道NJL 模型兩點函數KΨ有如下發(fā)散行為:

    依據(23)式可知,當固定p10時,KΨ(p) 作為p0的函數在p0±m(xù)h處會有奇點結構.為驗證這一奇點結構,可通過量子計算的經典模擬來實現3.2 節(jié)所提量子算法,先計算坐標空間兩點關聯(lián)函數,再通過傅里葉變換得到KΨ(p).圖10 展示了以p0為自變量的兩點關聯(lián)函數 TrKΨ(p) 的實部.圖中在p0a±1.18 附近處的峰對應著夸克質量峰.峰位置的準確性已經由精確對角化格點NJL 模型所給出的強子質量譜所驗證.在連續(xù)極限下,這個峰的高度將變?yōu)闊o限大,對應著奇點結構;對于有限格距情況下,這個峰只能有有限的高度.

    圖10 1+1維NJL 模型 TrKΨ(p) 的實部.圖中以 p0 為變量,固定p1 0Fig.10.Real part of TrKΨ(p) for the 1+1-dimensional NJL model as a function of p0a with p1 0.

    NJL 模型的這一簡單例子表明量子算法成功實現了傳播子和n點連通關聯(lián)函數所期望的奇點結構,這對LSZ 約化公式在量子計算機上的實現至關重要.

    4 有限溫格點施溫格模型的量子模擬

    核物質相圖的確定是高能核物理研究中的一大挑戰(zhàn),它依賴于對有限溫度有限密度QCD 系統(tǒng)平衡態(tài)的精確模擬,而格點QCD 方法在有限密度情形存在符號問題的困難.量子計算機在模擬有限溫有限密度平衡態(tài)系統(tǒng)上具有天然的優(yōu)勢,目前已經用于物態(tài)結構的探索[22-26].由于近期量子計算機資源的限制,目前有限溫有限密度QCD 的量子模擬還有待發(fā)展.本節(jié)將通過另一個較簡單的模型來闡述有限溫核物質的模擬,特別是在高溫情形時的夸克解禁閉行為.該模型為施溫格模型,它描述了1+1 維量子電動力學,包含了費米子與U(1)規(guī)范場的耦合.與QCD 類似,施溫格模型存在夸克禁閉[87,88]、手征對稱性自發(fā)破缺[89,90]等特征.因此,施溫格模型被廣泛用于模擬高能核物理問題.

    在選定規(guī)范A00 情況下,施溫格模型的哈密頓量可表達為[91]

    其中Ψ(x) 是費米子場,是電磁矢勢的縱向分量,m是裸費米子質量.

    由于局域規(guī)范不變性,電場和電荷密度算符之間存在約束?1E(x)gΨ(x)?Ψ(x),即滿足高斯定理.在1+1 維時,電場可以通過費米子來表示.這是施溫格模型的一個重要特點,即可表示為長程相互作用的費米子模型,規(guī)范場不出現,從而非常方便用量子比特來模擬.采用交錯費米子[92],并通過Jordan-Wigner 變換[75],施溫格模型可寫為

    其中ε是背景電場,?1/2a.更進一步,若要研究有限密度勢情形,則需在(25)式加一個額外的項 -[93],其中μ為化學勢.

    在零溫時,施溫格模型中正反夸克間具有與長程QCD 相似的線性勢[94],此即夸克禁閉.可看作夸克之間由一根弦連接,弦張力可衡量禁閉強度[95],它會隨著熱漲落的增強而減小.弦張力定義為系統(tǒng)自由能在零和非零背景電場下的差別[96]:

    其中β1/T是溫度的倒數,Fε(β) 是對應哈密頓量Hε的自由能.從弦張力的定義可知,計算弦張力只需計算自由能的差.計算自由能的量子算法已在2.2 節(jié)中給出.我們用開源程序包 QuTip[97]來進行量子計算的經典模擬.為了驗證量子算法的正確性,在結果中把量子模擬的結果與精確對角化的結果進行比對.首先模擬了施溫格模型弦張力對溫度的依賴關系,如圖12(a)所示,弦張力在小β區(qū)間會迅速衰減,這表明弦張力會因熱漲落效應變強而減小.特別地,弦張力在β0,即溫度無限高時趨于零,表明施溫格模型的禁閉-退禁閉轉變發(fā)生在無窮溫度.圖12(b)表明弦張力的大小隨著溫度的增加指數下降,且ε越小下降得越快.弦張力的指數衰減行為與理論預言相一致[96].

    圖12 不同背景電場 ε 下的弦張力.其中給定 m1, g 1, ? 1,N=6(a) 弦張力對 β 的依賴;(b) 弦張力的對數對溫度 T 的依賴Fig.12.String tensions under different ε in the case m1, g 1, ? 1,N=6:(a) The string tension as a function of the inverse temperature β ;(b) logarithm of the string tension as a function of the temperature T.

    進一步對有限溫有限密度情形下的弦張力進行了計算.圖13(a)給出的是弦張力對不同化學勢的依賴關系,可發(fā)現在相同溫度下,化學勢越大對應的弦張力越小.圖13(b)給出了弦張力對溫度和化學勢兩個變量的依賴關系.圖中可見隨著溫度和化學勢的升高弦張力都會減小,這與QCD 相圖的行為類似[98-101].然而,與QCD 不同的是,施溫格模型弦張力只有在T →∞時才會下降到零,這表明在有限溫情形下施溫格模型相圖不存在禁閉-退禁閉相變.盡管如此,施溫格模型相圖仍能從定性上看出禁閉-退禁閉相變的一些性質,并且成功展示出量子計算模擬有限密度區(qū)域相結構的可行性.

    圖13 在給定參數 ε0.5, m1, g 1, ? 1, N 6 下的弦張力(a) 不同化學勢 μ 下弦張力對溫度 T 的依賴;(b) 弦張力對溫度 T 和化學勢μ 的依賴Fig.13.The string tension in the case ε0.5, m1, g 1, ? 1,N 6 :(a) At different μ,the string tension as a function of the temperature T ;(b) the string tension as a function of the temperature T and the chemical potential μ.

    5 總結

    本文對我們近期基于量子計算機來模擬高能核物理的系列工作進行了綜述.針對強子散射問題,系統(tǒng)提出了以多項式復雜度計算PDFs和LCDAs 的量子算法,其中包括量子數分辨的變分量子本征求解器以及光錐關聯(lián)函數的測量兩大部分,并通過NJL 模型的數值模擬,初步展示了量子計算在定性上已經能給出合理的結果.此外,基于LSZ 約化公式,提出了強子散射振幅的量子計算新方案,避免了初態(tài)波包制備的困難,并適用于束縛態(tài)粒子的散射問題.還對有限溫有限密度核物質的量子模擬進行了探索.通過施溫格模型,展示了量子計算機可以模擬夸克在高溫情形出現的禁閉-解禁閉行為.本文的研究初步展示了量子計算在解決高能核物理中含時問題和符號問題的可行性及優(yōu)勢,有助于進一步探索物質深層次結構及核物質相圖.本文提出的強子態(tài)制備及結構探測是普適的量子算法框架,有望在高能核物理領域中廣泛運用,如模擬部分子碎裂函數(parton fragmentation functions,FFs),模擬核效應的產生機制等.同時,針對1+1 維施溫格模型發(fā)展的有限溫量子熱態(tài)的量子算法,有望拓展到2+1和3+1 維下非阿貝爾規(guī)范場,從而模擬有限溫有限密度情形的QCD 相圖,有望揭示QCD 禁閉相到退禁閉相的轉變.總而言之,量子計算為未來研究高能核物理中的含時非微擾問題以及有限溫有限密相結構問題提供了一個潛在的根本性解決方案.

    猜你喜歡
    核物理夸克強子
    大型強子對撞機發(fā)現3 種新奇異粒子
    與塞西爾?德維特-莫雷塔一席談她與中國兩彈元勛彭桓武的戀情
    宇宙中的“煉金術”
    科學技術創(chuàng)新(2020年36期)2020-12-15 08:36:48
    罰你沒商量
    故事會(2020年12期)2020-06-19 08:52:36
    《核物理介紹》教學中構建學生自主學習激勵機制的探索性研究
    我坐在一塊石頭上(外五首)
    都市(2018年12期)2018-09-10 10:13:17
    兩大團隊同時預測四夸克重粒子存在
    核物理研究在我國的主要發(fā)展動向
    利用BS方程研究由兩個粲夸克構成的雙夸克的第一激發(fā)態(tài)
    国产免费视频播放在线视频| 亚洲国产欧美一区二区综合| 午夜激情久久久久久久| 老司机影院成人| 久久人人97超碰香蕉20202| 久久毛片免费看一区二区三区| 天堂中文最新版在线下载| 亚洲欧美成人综合另类久久久| 丰满人妻熟妇乱又伦精品不卡| 大片免费播放器 马上看| 老司机靠b影院| 午夜福利一区二区在线看| 性色av乱码一区二区三区2| bbb黄色大片| 一级毛片女人18水好多| 男人舔女人的私密视频| 99re6热这里在线精品视频| 满18在线观看网站| 十分钟在线观看高清视频www| 国产1区2区3区精品| 精品少妇一区二区三区视频日本电影| 80岁老熟妇乱子伦牲交| 国产黄色免费在线视频| 国产日韩一区二区三区精品不卡| av免费在线观看网站| 欧美av亚洲av综合av国产av| 免费女性裸体啪啪无遮挡网站| 亚洲国产av影院在线观看| 窝窝影院91人妻| 国产一级毛片在线| 欧美日韩中文字幕国产精品一区二区三区 | 国产亚洲精品第一综合不卡| 精品国产乱码久久久久久男人| 午夜老司机福利片| 亚洲专区字幕在线| 国产精品麻豆人妻色哟哟久久| 国产精品自产拍在线观看55亚洲 | 午夜福利一区二区在线看| a级片在线免费高清观看视频| 99国产精品一区二区蜜桃av | 18禁黄网站禁片午夜丰满| 久久综合国产亚洲精品| 欧美日韩亚洲高清精品| 一二三四社区在线视频社区8| 满18在线观看网站| 亚洲欧洲日产国产| 女性被躁到高潮视频| av线在线观看网站| 天天添夜夜摸| 午夜老司机福利片| 中文字幕高清在线视频| 国产精品一二三区在线看| avwww免费| 午夜福利在线免费观看网站| 精品国产一区二区久久| 日日摸夜夜添夜夜添小说| 丰满饥渴人妻一区二区三| 久久亚洲精品不卡| 黄片大片在线免费观看| 久久99热这里只频精品6学生| 亚洲一区中文字幕在线| 亚洲成人免费av在线播放| 亚洲av美国av| 男女边摸边吃奶| 婷婷丁香在线五月| 国产1区2区3区精品| 狂野欧美激情性bbbbbb| 国产免费av片在线观看野外av| 狠狠狠狠99中文字幕| 日韩大片免费观看网站| 欧美日韩av久久| 欧美97在线视频| 香蕉丝袜av| 日韩欧美一区二区三区在线观看 | 一边摸一边做爽爽视频免费| 人妻一区二区av| 国产精品成人在线| 9191精品国产免费久久| 久久国产精品男人的天堂亚洲| 亚洲欧美精品自产自拍| 制服人妻中文乱码| 国产欧美日韩一区二区三区在线| 国产伦理片在线播放av一区| 视频区图区小说| 大片免费播放器 马上看| 新久久久久国产一级毛片| 欧美激情极品国产一区二区三区| 一级片免费观看大全| 国产激情久久老熟女| 一级毛片电影观看| 亚洲成人免费电影在线观看| 十八禁网站网址无遮挡| 老鸭窝网址在线观看| 欧美变态另类bdsm刘玥| 成人av一区二区三区在线看 | 老鸭窝网址在线观看| 国产熟女午夜一区二区三区| 成人三级做爰电影| 国产精品 欧美亚洲| 久久 成人 亚洲| 欧美日韩精品网址| 亚洲专区中文字幕在线| avwww免费| 女人久久www免费人成看片| 伊人亚洲综合成人网| 日韩一卡2卡3卡4卡2021年| 欧美少妇被猛烈插入视频| 99国产精品99久久久久| 久久久精品区二区三区| 久久精品熟女亚洲av麻豆精品| 久久久国产精品麻豆| 中亚洲国语对白在线视频| 视频区图区小说| 法律面前人人平等表现在哪些方面 | 日日爽夜夜爽网站| 岛国在线观看网站| 国产高清视频在线播放一区 | 国产成人啪精品午夜网站| 日日夜夜操网爽| 老熟妇乱子伦视频在线观看 | 肉色欧美久久久久久久蜜桃| 久久性视频一级片| 一本—道久久a久久精品蜜桃钙片| 亚洲av电影在线观看一区二区三区| 亚洲自偷自拍图片 自拍| 国产欧美日韩一区二区三区在线| 久久亚洲国产成人精品v| 久久免费观看电影| 丝袜美足系列| 另类亚洲欧美激情| 免费在线观看黄色视频的| 精品久久久精品久久久| 久久狼人影院| 两性夫妻黄色片| 久久ye,这里只有精品| 国产精品久久久久成人av| 韩国精品一区二区三区| 日韩视频一区二区在线观看| 黑人猛操日本美女一级片| 精品国产超薄肉色丝袜足j| 美女高潮到喷水免费观看| 国产极品粉嫩免费观看在线| 国产日韩一区二区三区精品不卡| 欧美人与性动交α欧美软件| 成人18禁高潮啪啪吃奶动态图| 啦啦啦中文免费视频观看日本| 中文字幕制服av| 黄色视频,在线免费观看| 桃花免费在线播放| 午夜日韩欧美国产| 午夜两性在线视频| 国产野战对白在线观看| 另类亚洲欧美激情| 亚洲av电影在线进入| 精品少妇黑人巨大在线播放| 久久人妻福利社区极品人妻图片| 菩萨蛮人人尽说江南好唐韦庄| 免费在线观看日本一区| 动漫黄色视频在线观看| 日韩欧美一区视频在线观看| 日本a在线网址| 夜夜夜夜夜久久久久| 精品少妇内射三级| 韩国高清视频一区二区三区| 久久亚洲国产成人精品v| 99国产精品一区二区蜜桃av | 男人舔女人的私密视频| 久久亚洲国产成人精品v| 亚洲黑人精品在线| 亚洲免费av在线视频| 美女福利国产在线| 妹子高潮喷水视频| 久久久精品免费免费高清| 精品免费久久久久久久清纯 | 国产又色又爽无遮挡免| 国产福利在线免费观看视频| 狂野欧美激情性xxxx| 99久久99久久久精品蜜桃| 久久人人97超碰香蕉20202| 亚洲国产精品一区二区三区在线| 亚洲精品久久午夜乱码| 国产一区有黄有色的免费视频| 亚洲国产欧美一区二区综合| 精品少妇内射三级| 欧美黄色片欧美黄色片| 老司机福利观看| 天堂俺去俺来也www色官网| 国产精品一区二区在线观看99| 国产男女内射视频| 性色av一级| 黑丝袜美女国产一区| www.精华液| 飞空精品影院首页| 电影成人av| 真人做人爱边吃奶动态| 久久精品国产亚洲av香蕉五月 | 亚洲av欧美aⅴ国产| 久久久久久久精品精品| 欧美变态另类bdsm刘玥| 中文字幕色久视频| 狠狠精品人妻久久久久久综合| 国产欧美亚洲国产| 亚洲精品美女久久av网站| 国产欧美日韩一区二区三区在线| 午夜福利乱码中文字幕| 男女之事视频高清在线观看| 午夜老司机福利片| 中文字幕av电影在线播放| 性色av一级| 国产一卡二卡三卡精品| 久久女婷五月综合色啪小说| 欧美亚洲 丝袜 人妻 在线| 777米奇影视久久| 超碰成人久久| 精品亚洲成国产av| 纵有疾风起免费观看全集完整版| 性少妇av在线| 丝袜美足系列| 亚洲精品国产区一区二| 国产一卡二卡三卡精品| 俄罗斯特黄特色一大片| 一区二区av电影网| 国产成人免费无遮挡视频| 国产成人欧美| 99精品欧美一区二区三区四区| 两性夫妻黄色片| 在线观看免费高清a一片| av天堂在线播放| 一级,二级,三级黄色视频| netflix在线观看网站| 久久天堂一区二区三区四区| 久久久久国产一级毛片高清牌| 蜜桃国产av成人99| 亚洲精品自拍成人| 国产老妇伦熟女老妇高清| 99精品欧美一区二区三区四区| 久久这里只有精品19| 国产精品1区2区在线观看. | 亚洲熟女毛片儿| 人人妻人人澡人人爽人人夜夜| 男女高潮啪啪啪动态图| 日韩视频在线欧美| 十分钟在线观看高清视频www| 婷婷色av中文字幕| 中国美女看黄片| 一本综合久久免费| 19禁男女啪啪无遮挡网站| 夫妻午夜视频| 久久国产精品男人的天堂亚洲| 国产熟女午夜一区二区三区| 一区二区三区激情视频| 电影成人av| 丰满少妇做爰视频| 欧美日本中文国产一区发布| 亚洲精品国产色婷婷电影| 99国产极品粉嫩在线观看| 免费人妻精品一区二区三区视频| 免费在线观看视频国产中文字幕亚洲 | 中文欧美无线码| 久久久久网色| 黑丝袜美女国产一区| 人人妻人人添人人爽欧美一区卜| www.999成人在线观看| 天堂中文最新版在线下载| 国精品久久久久久国模美| 热99久久久久精品小说推荐| 国产亚洲午夜精品一区二区久久| 欧美精品高潮呻吟av久久| 中文字幕人妻熟女乱码| 日韩一区二区三区影片| 性色av乱码一区二区三区2| 国产1区2区3区精品| 考比视频在线观看| 免费久久久久久久精品成人欧美视频| 国产一区有黄有色的免费视频| 人人妻人人澡人人看| 建设人人有责人人尽责人人享有的| 丝袜美腿诱惑在线| 久久精品国产综合久久久| 99精品欧美一区二区三区四区| 亚洲国产中文字幕在线视频| 在线观看一区二区三区激情| 99热网站在线观看| 久久精品国产a三级三级三级| 久久精品亚洲av国产电影网| 精品少妇一区二区三区视频日本电影| 电影成人av| 又黄又粗又硬又大视频| 欧美日韩亚洲综合一区二区三区_| 窝窝影院91人妻| 亚洲视频免费观看视频| 男女国产视频网站| 两个人免费观看高清视频| 老司机午夜福利在线观看视频 | 精品一区二区三卡| 久久久久国产一级毛片高清牌| 韩国精品一区二区三区| 女警被强在线播放| 亚洲专区中文字幕在线| 国产欧美日韩综合在线一区二区| 老汉色∧v一级毛片| 亚洲精品日韩在线中文字幕| 曰老女人黄片| 永久免费av网站大全| 一本久久精品| 99精国产麻豆久久婷婷| 91九色精品人成在线观看| 亚洲精品粉嫩美女一区| 日日夜夜操网爽| 人人妻人人爽人人添夜夜欢视频| 黑人欧美特级aaaaaa片| 99精品欧美一区二区三区四区| 黄色 视频免费看| 国产精品 国内视频| 欧美 亚洲 国产 日韩一| 国产av精品麻豆| 国产成人免费无遮挡视频| 91精品三级在线观看| 久久午夜综合久久蜜桃| 18禁国产床啪视频网站| 国产日韩一区二区三区精品不卡| 亚洲成人国产一区在线观看| 欧美另类一区| 午夜福利一区二区在线看| 午夜激情久久久久久久| 日韩精品免费视频一区二区三区| 午夜激情久久久久久久| 国产一区有黄有色的免费视频| 免费不卡黄色视频| 在线观看人妻少妇| 久久久久视频综合| 在线观看一区二区三区激情| 后天国语完整版免费观看| 午夜激情av网站| 国产成人精品久久二区二区91| 我要看黄色一级片免费的| 亚洲人成77777在线视频| 久久99热这里只频精品6学生| 少妇被粗大的猛进出69影院| 中文字幕最新亚洲高清| 亚洲精品一二三| 免费av中文字幕在线| 日韩精品免费视频一区二区三区| 国产亚洲午夜精品一区二区久久| 国产精品久久久人人做人人爽| 国产激情久久老熟女| 青春草视频在线免费观看| 久久久久国产一级毛片高清牌| 丝袜喷水一区| 777米奇影视久久| 精品亚洲成a人片在线观看| 国产麻豆69| 欧美+亚洲+日韩+国产| 麻豆国产av国片精品| 天堂俺去俺来也www色官网| 一级片免费观看大全| 手机成人av网站| cao死你这个sao货| 一级a爱视频在线免费观看| 日本a在线网址| 免费一级毛片在线播放高清视频 | 精品一品国产午夜福利视频| 久久人妻福利社区极品人妻图片| 午夜久久久在线观看| 91字幕亚洲| 超色免费av| 热re99久久国产66热| 国产精品国产av在线观看| 女警被强在线播放| 国产精品久久久久久精品电影小说| 亚洲精品国产精品久久久不卡| 操美女的视频在线观看| 9热在线视频观看99| 91成年电影在线观看| 国产男女内射视频| 国产一级毛片在线| 黄色视频,在线免费观看| 99热网站在线观看| 亚洲熟女毛片儿| 啦啦啦 在线观看视频| 首页视频小说图片口味搜索| 欧美 日韩 精品 国产| 国内毛片毛片毛片毛片毛片| 免费在线观看日本一区| 搡老岳熟女国产| 69精品国产乱码久久久| 欧美乱码精品一区二区三区| 1024香蕉在线观看| 久久久国产成人免费| 中文字幕精品免费在线观看视频| 女人高潮潮喷娇喘18禁视频| 伊人亚洲综合成人网| 看免费av毛片| 精品一品国产午夜福利视频| 精品国产一区二区久久| 黄色视频在线播放观看不卡| 日韩三级视频一区二区三区| 精品欧美一区二区三区在线| 亚洲 国产 在线| 亚洲av国产av综合av卡| 一区二区三区四区激情视频| 久久狼人影院| 99精品久久久久人妻精品| 精品少妇内射三级| 午夜福利乱码中文字幕| 精品久久久久久电影网| 丰满饥渴人妻一区二区三| 亚洲av日韩在线播放| 真人做人爱边吃奶动态| 亚洲九九香蕉| www日本在线高清视频| 啪啪无遮挡十八禁网站| 丰满饥渴人妻一区二区三| 免费观看a级毛片全部| 久久久久久人人人人人| 国产麻豆69| tube8黄色片| 91麻豆av在线| 国产91精品成人一区二区三区 | 美女福利国产在线| 黑丝袜美女国产一区| 超碰成人久久| 免费看十八禁软件| 亚洲国产欧美日韩在线播放| 日韩制服骚丝袜av| 少妇精品久久久久久久| 欧美中文综合在线视频| 宅男免费午夜| 久久 成人 亚洲| 久久女婷五月综合色啪小说| 免费日韩欧美在线观看| a在线观看视频网站| 丝袜在线中文字幕| 国产精品秋霞免费鲁丝片| 午夜免费成人在线视频| 777米奇影视久久| 国产免费av片在线观看野外av| 国产精品一区二区免费欧美 | 免费在线观看视频国产中文字幕亚洲 | 黄片播放在线免费| 人人妻人人澡人人爽人人夜夜| 亚洲精品久久久久久婷婷小说| 一二三四在线观看免费中文在| 亚洲人成电影观看| 亚洲中文字幕日韩| 一进一出抽搐动态| 国产精品麻豆人妻色哟哟久久| 夫妻午夜视频| 国产一区二区三区综合在线观看| 欧美乱码精品一区二区三区| 电影成人av| 热re99久久精品国产66热6| 久久精品久久久久久噜噜老黄| 十八禁高潮呻吟视频| 成人三级做爰电影| 一本一本久久a久久精品综合妖精| 黄色毛片三级朝国网站| 亚洲欧美一区二区三区黑人| 国产精品一二三区在线看| 操出白浆在线播放| 午夜影院在线不卡| 久久午夜综合久久蜜桃| 女性生殖器流出的白浆| 成年人午夜在线观看视频| 国产黄色免费在线视频| 91字幕亚洲| 亚洲第一欧美日韩一区二区三区 | 视频区欧美日本亚洲| 午夜福利影视在线免费观看| 午夜福利在线观看吧| 国产无遮挡羞羞视频在线观看| 久久久精品国产亚洲av高清涩受| 精品视频人人做人人爽| 精品少妇久久久久久888优播| 老司机影院毛片| 成年人午夜在线观看视频| 国产黄频视频在线观看| av线在线观看网站| 夫妻午夜视频| 大香蕉久久网| 在线十欧美十亚洲十日本专区| 男女边摸边吃奶| 欧美日韩亚洲高清精品| 婷婷丁香在线五月| 99精品欧美一区二区三区四区| netflix在线观看网站| 欧美精品人与动牲交sv欧美| 午夜视频精品福利| 自拍欧美九色日韩亚洲蝌蚪91| 国产伦理片在线播放av一区| 午夜91福利影院| av超薄肉色丝袜交足视频| 老汉色∧v一级毛片| 最新在线观看一区二区三区| 国产一卡二卡三卡精品| 精品高清国产在线一区| 亚洲国产欧美一区二区综合| 亚洲全国av大片| 中文精品一卡2卡3卡4更新| 亚洲国产精品999| 亚洲国产欧美网| 精品福利永久在线观看| 三上悠亚av全集在线观看| 最新在线观看一区二区三区| 久久中文看片网| 亚洲av成人一区二区三| 纵有疾风起免费观看全集完整版| 亚洲国产欧美日韩在线播放| 亚洲国产欧美在线一区| 在线精品无人区一区二区三| 99国产精品一区二区三区| 久久99一区二区三区| 国产亚洲一区二区精品| 亚洲伊人久久精品综合| 热99久久久久精品小说推荐| 黑人操中国人逼视频| 亚洲黑人精品在线| 999久久久国产精品视频| 91麻豆av在线| 国产欧美日韩精品亚洲av| 我要看黄色一级片免费的| 国产精品偷伦视频观看了| 精品亚洲成a人片在线观看| 搡老乐熟女国产| 国产亚洲欧美在线一区二区| a 毛片基地| 在线观看免费午夜福利视频| 久久精品熟女亚洲av麻豆精品| av在线播放精品| 男女下面插进去视频免费观看| 黄色 视频免费看| 中文精品一卡2卡3卡4更新| 久久久久精品人妻al黑| 亚洲一区中文字幕在线| 操出白浆在线播放| 久热这里只有精品99| 亚洲中文日韩欧美视频| 亚洲精品自拍成人| 97在线人人人人妻| 久久热在线av| 男女之事视频高清在线观看| 美女视频免费永久观看网站| 亚洲综合色网址| 免费在线观看完整版高清| 咕卡用的链子| 黄片大片在线免费观看| 国产野战对白在线观看| 黄片播放在线免费| 亚洲av日韩在线播放| 午夜视频精品福利| 国产精品 欧美亚洲| 午夜福利在线免费观看网站| 亚洲欧美激情在线| 亚洲成国产人片在线观看| www.999成人在线观看| 久久久久网色| 三级毛片av免费| 亚洲视频免费观看视频| 亚洲精品国产av成人精品| 成人手机av| 欧美在线黄色| 久久国产亚洲av麻豆专区| 中国美女看黄片| 国产野战对白在线观看| 高潮久久久久久久久久久不卡| 51午夜福利影视在线观看| 悠悠久久av| 法律面前人人平等表现在哪些方面 | 伦理电影免费视频| 久久中文字幕一级| 丝袜在线中文字幕| 午夜成年电影在线免费观看| 老司机影院毛片| 在线观看www视频免费| 精品欧美一区二区三区在线| 美女高潮到喷水免费观看| 侵犯人妻中文字幕一二三四区| 国产97色在线日韩免费| 亚洲精品美女久久av网站| 九色亚洲精品在线播放| 日韩制服骚丝袜av| 亚洲男人天堂网一区| 国产老妇伦熟女老妇高清| 亚洲欧美成人综合另类久久久| 蜜桃在线观看..| 不卡av一区二区三区| 亚洲人成电影免费在线| 热99久久久久精品小说推荐| 中文字幕精品免费在线观看视频| 99精品久久久久人妻精品| 国产又色又爽无遮挡免| 欧美 亚洲 国产 日韩一| 免费在线观看黄色视频的| 亚洲视频免费观看视频| 久久精品亚洲熟妇少妇任你| 99热国产这里只有精品6| av视频免费观看在线观看| 国产在视频线精品| 欧美日韩一级在线毛片| av欧美777| 搡老熟女国产l中国老女人| tocl精华| 少妇被粗大的猛进出69影院| 中文字幕精品免费在线观看视频| 免费在线观看黄色视频的| 亚洲精品一二三| 在线 av 中文字幕| 女警被强在线播放| 久久av网站| 黄色视频在线播放观看不卡| 亚洲免费av在线视频| 欧美激情极品国产一区二区三区| 99热全是精品| 国产亚洲欧美精品永久| 99精国产麻豆久久婷婷|