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

    質(zhì)子碰撞硼原子非輻射的電荷轉(zhuǎn)移過程*

    2023-09-06 02:15:20朱宇豪袁翔吳勇王建國
    物理學(xué)報 2023年16期
    關(guān)鍵詞:質(zhì)子組態(tài)原子

    朱宇豪 袁翔 吳勇 王建國

    1) (西安建筑科技大學(xué)理學(xué)院,西安 710055)

    2) (北京應(yīng)用物理與計算數(shù)學(xué)研究所,計算物理重點實驗室,北京 100084)

    重粒子碰撞中的電子轉(zhuǎn)移涉及復(fù)雜的電子關(guān)聯(lián)機(jī)制,極大地影響等離子體中電荷態(tài)平衡,也是X 射線的輻射的重要來源之一.電子轉(zhuǎn)移截面與速率系數(shù)是國防工業(yè)發(fā)展核聚變等離子體所需要的重要原子參數(shù).基于全量子的非輻射分子軌道密耦合方法,系統(tǒng)研究了質(zhì)子碰撞硼原子在 10-3—103 eV/u 能量區(qū)間內(nèi)的硼原子電子轉(zhuǎn)移過程.計算采用多參考組態(tài)方法得到總共15 個電子轉(zhuǎn)移、激發(fā)以及彈性散射的通道,每個通道對應(yīng)的分子態(tài)能量得到了與實驗符合較好的結(jié)果.分子態(tài)的絕熱勢能曲線間的避免交叉現(xiàn)象明顯,構(gòu)成了電子轉(zhuǎn)移的主要途徑.計算發(fā)現(xiàn),質(zhì)子碰撞硼原子過程中2s 軌道的電子轉(zhuǎn)移是占主導(dǎo)地位,2p 軌道的電子轉(zhuǎn)移貢獻(xiàn)較小.在低能區(qū),電子轉(zhuǎn)移截面出現(xiàn)了明顯的量子共振現(xiàn)象,這些共振主要來源于不同能量通道的耦合.此外,還計算了不同溫度下的質(zhì)子碰撞硼原子的電子轉(zhuǎn)移速率,該速率可為復(fù)雜等離子體環(huán)境的模擬診斷提供重要的原子參數(shù)支持.

    1 引言

    重粒子碰撞中的電子轉(zhuǎn)移涉及多電荷動力學(xué)以及多體相互作用等復(fù)雜物理過程,因為其在天體和等離子體物理中有非常重要的應(yīng)用,吸引了廣泛的理論和實驗研究[1].在磁性約束核聚變中,質(zhì)子對電子的俘獲極大影響了整個系統(tǒng)的電荷態(tài)平衡以及等離子體冷卻過程.電荷態(tài)平衡是診斷復(fù)雜等離子體系統(tǒng)基礎(chǔ)的指標(biāo),從中可以獲得元素豐度、電荷密度等物理信息[2].在天體環(huán)境中,X 射線譜是識別天體可觀測氣體元素重要的手段之一,并且各類元素的X 射線輻射都與電子轉(zhuǎn)移有密切的關(guān)系.特別對于內(nèi)殼層電子的轉(zhuǎn)移過程,是原子系統(tǒng)X 射線輻射的重要來源[3].所以,為了模擬上述X 射線譜和等離子體中電荷態(tài)平衡等過程,獲得重粒子碰撞中的電子轉(zhuǎn)移態(tài)選擇截面和總截面是很有意義的.那么,重粒子碰撞的電子轉(zhuǎn)移截面等相關(guān)信息已經(jīng)成為不可或缺的原子參數(shù),為國防工業(yè)發(fā)展核聚變技術(shù),以及為等離子體物理等學(xué)科中演化模擬復(fù)雜等離子體系統(tǒng)提供強(qiáng)有力的原子參數(shù)支持[4,5].

    近年來,對重粒子碰撞的理論研究也迅猛發(fā)展,在不同的入射能區(qū),理論計算可以做到與實驗測量符合較為精確的結(jié)果.如在中高能區(qū)的原子軌道密耦合方法[6?9],能夠精確計算重粒子碰撞的電子激發(fā)和轉(zhuǎn)移的態(tài)選擇截面以及總截面.因為其只對電子做全量子化的處理,而原子核依然采用經(jīng)典力學(xué)的運動方式,所以該方法為半經(jīng)典的理論方法.在低能入射區(qū),量子效應(yīng)更加明顯,需要采用全量子的碰撞方法才能考慮足夠的物理信息得到精確的電子轉(zhuǎn)移截面.目前,較為常見的全量子的理論方法是非輻射分子軌道密耦合方法,Liu 等[6,10?14]基于此方法成功計算出了He,Li 以及Be 等原子系統(tǒng)的電子轉(zhuǎn)移截面,并與相關(guān)精密實驗測量結(jié)果符合得較好.全量子的分子軌道方法,因為分子軌道模型可以準(zhǔn)確描述較低能區(qū)的散射態(tài)波函數(shù),并能夠有效處理各個分子軌道之間的耦合,其在計算低能區(qū)對電子轉(zhuǎn)移的截面計算屬于國內(nèi)外公認(rèn)的最高精度方法.當(dāng)然,分子軌道方法的精確度主要取決于分子態(tài)能量,即分子電子態(tài)結(jié)構(gòu)的精度.這要求在計算分子結(jié)構(gòu)時,需考慮足夠大的組態(tài)空間和足夠多的電子組態(tài),充分考慮電子關(guān)聯(lián)對分子結(jié)構(gòu)的影響[14].

    在實驗研究方面,重粒子碰撞是闡明量子多體關(guān)聯(lián),揭示靶粒子內(nèi)部結(jié)構(gòu)信息的重要科學(xué)手段.國內(nèi)外已經(jīng)開展了數(shù)十年的重粒子碰撞實驗,諸多質(zhì)子或氫原子與重元素原子不同價態(tài)系統(tǒng)的碰撞實驗都相繼展開.比如,Bruhns 等[15,16]使用合并束技術(shù)研究了氫原子與Si3+離子在地能區(qū)碰撞的電子轉(zhuǎn)移過程,并精密測量了電子轉(zhuǎn)移截面.Gieler等[17,18]研究了質(zhì)子碰撞K 原子的低能區(qū)電子轉(zhuǎn)移總截面,其中Gieler 等[17]的碰撞能區(qū)在 1 —10 keV/u,Ebel 等[18]的研究能區(qū)在 0 .1—1.0 keV/u .雖然諸多原子系統(tǒng)與質(zhì)子碰撞中電子轉(zhuǎn)移過程被實驗陸續(xù)研究[19-22],但是據(jù)目前調(diào)研所知,還從未開展對B原子的實驗研究,然而B 原子在天體中也具有較高豐度,因此迫切需要B 原子參數(shù).此外,B 原子與質(zhì)子的碰撞過程也是研究多體相互作用系統(tǒng)的良好平臺,其中涉及多體相互作用的電子轉(zhuǎn)移機(jī)制還沒有被發(fā)掘.因此,質(zhì)子碰撞B 原子的理論研究,不僅可以更深一步了解其電子轉(zhuǎn)移機(jī)制,還能認(rèn)識碰撞過程中多電子關(guān)聯(lián)效應(yīng).

    本文基于全量子非輻射的分子軌道密耦合方法,系統(tǒng)研究大范圍碰撞能區(qū)內(nèi)質(zhì)子碰撞B 原子電子轉(zhuǎn)移過程.考慮足夠多的電子態(tài),以確保整個計算中通道的收斂性.第2 節(jié)詳細(xì)闡述計算的理論方法;第3 節(jié)展示計算結(jié)果并討論相關(guān)電子轉(zhuǎn)移的物理機(jī)制;第4 節(jié)對研究工作做簡短的總結(jié).

    2 理論方法

    2.1 分子電子態(tài)結(jié)構(gòu)計算

    分子電子態(tài)結(jié)構(gòu)的計算,采用Buenker 等[23?25]提出的Table-direct CI 結(jié)合多參考單雙激發(fā)組態(tài)相互作用方法(MRDCI).該方法在自洽場計算(HF-SCF)之后,能基于單雙激發(fā)的參考組態(tài),可以最優(yōu)化選取占比重要的組態(tài),縮小了組態(tài)空間,且可以保證精度.則說明MRDCI 可以快速優(yōu)化組態(tài)構(gòu)成,得到精確的電子態(tài)波函數(shù),其具體計算公式為

    其中,ψα(PJM) 是電子態(tài)波函數(shù),?r(PJM) 為參考組態(tài)考慮單雙激發(fā)后的優(yōu)化的組態(tài)波函數(shù),cr(r)是混合系數(shù).實際計算中,通常選取對內(nèi)、外殼層和極化作用考慮較好的芯-價關(guān)聯(lián)極化的多Zeta 基組,加入彌散基來彌補(bǔ)由于高激發(fā)態(tài)離核較遠(yuǎn)導(dǎo)致電子云分布較大的影響.做完自洽場計算后,構(gòu)建組態(tài)空間,對于選定的參考組態(tài)進(jìn)行單雙電子激發(fā),但是這樣經(jīng)常會導(dǎo)致組態(tài)空間偏大,因此可以通過設(shè)定閾值、結(jié)合“有效勢場”控制活動電子數(shù)量來縮小組態(tài)空間,選取貢獻(xiàn)較大的組態(tài).

    2.2 重粒子碰撞動力學(xué)計算

    非輻射分子軌道密耦合方法(QMOCC)對于重離子散射系統(tǒng),采用質(zhì)心坐標(biāo)和原子單位,兩原子體系碰撞系統(tǒng)的總哈密頓量表示為

    式中,R為入射離子相對靶粒子核間距,V為所有電子以及原子核的庫侖相互作用之和,質(zhì)量項M=M1+M2,μ=M1M2/M,m=M/(M+1),注意(2)式忽略了極化相互作用項.

    令?i為電子總波函數(shù),在R→∞時候,電子相對于原子核運動,以核為中心,電子總波函數(shù)能在絕熱表象下有如下形式:

    式中,e xp[iγ(R,r)] 稱電子轉(zhuǎn)移因子.在微擾的穩(wěn)態(tài)近似下,可以忽略電子轉(zhuǎn)移因子,把系統(tǒng)的總波函數(shù)近似用絕熱電子波函數(shù)代替,利用Born-Oppenheimer 近似,則總波函數(shù)可以用電子絕熱波函數(shù)展開:

    其中對i的求和限于完備的絕熱矢量集合{ψi}中的子集.而且,每一個ψi在漸進(jìn)區(qū)對于這不同的原子態(tài),這些原子態(tài)也稱為通道.

    將(4)式代入能量本征方程,可得到微擾的穩(wěn)態(tài)方程:

    式中,I為單位矩陣,u為對角的絕熱勢矩陣,uij=εjδij,εj是絕熱態(tài)的本征能量,M(R) 和>P(R)是耦合矩陣分別滿足:

    在球坐標(biāo)系里,(5)式可以寫為

    對于求解(12)式,定義一個新的函數(shù)gl,滿足fl=Cgl,其中C是正交變換矩陣.若f是(12)式的一個解,那么g同樣是另一個解,因此可以這樣構(gòu)建C:

    式中u是對角的絕熱勢矩陣,Uγ,γ′是非絕熱勢矩陣,P是轉(zhuǎn)動耦合矩陣.(14)式滿足以下邊界條件:

    其中jJ和ηJ分別為通道的規(guī)則和非規(guī)則的庫侖函數(shù),Kl為實對稱的K矩陣.從K矩陣可以得到S矩陣:

    其中Sl表示第l個分波對應(yīng)的S矩陣.

    對于速率系數(shù)的計算,采用Maxwellian-Boltzmann 分布對散射截面求平均,那么反應(yīng)速率為

    其中,k為玻爾茲曼常數(shù),T為溫度,μ是約化質(zhì)量.

    3 計算結(jié)果和討論

    對于BH+分子的電子態(tài),總共計算了9 個2Σ+態(tài)和6 個2Π 態(tài),這15 個電子態(tài)的絕熱勢能曲線如圖1 所示,橫坐標(biāo)為原子核間距,縱坐標(biāo)是電子態(tài)能量.從圖1(a)可看到較低的5 個分子態(tài)之間有較大的能量差,較高的激發(fā)態(tài)能級比較密集,處于近簡并的狀態(tài).電子態(tài)之間的避免交叉近區(qū)相護(hù)作用很強(qiáng)烈,這些避免交叉現(xiàn)象在耦合矩陣元中也有明顯的表現(xiàn).圖1(b)中,在R=3.8 ? 處,32Π 和42Π態(tài)有一個顯著的避免交叉,而分別在R=2.3,3.2,4.6 ? 處,42Π 和 52Π 也有比較明顯的避免交叉現(xiàn)象.對于QMOCC 的計算,這些避免交叉能量相近的點起到關(guān)鍵的作用,因為絕熱分子態(tài)的耦合構(gòu)成了非輻射電子轉(zhuǎn)移的主要途徑.

    圖1 (a) BH+分子 Σ 態(tài)勢能曲線;(b) BH+分子 Π 態(tài)勢能曲線Fig.1.(a) Energy curve of BH+ molecule with Σ states;(b) energy curve of BH+ molecule with Π states.

    為了保證后續(xù)動力學(xué)計算的準(zhǔn)確性,表1 列出了這15 個電子態(tài)的理論結(jié)果和實驗結(jié)果,以及各個不同分子態(tài)解離極限所對應(yīng)的原子態(tài).在QMOCC 的理論框架下,散射末態(tài)對應(yīng)的原子能級與分子態(tài)解離極限對應(yīng)的原子態(tài)一一對應(yīng).上述計算構(gòu)建的散射矩陣,可以清晰說明不同散射通道之間躍遷的物理信息.從表1 還可以看到,分子態(tài)的解離極限涉及的原子態(tài)能級最大到2p 電子軌道,電子從B 原子的轉(zhuǎn)移主要來源于2s 和2p 電子軌道.除了非彈性的散射通道之外,計算還包括一些彈性散射通道,如分子態(tài) 32Σ+和 22Π .這些彈性散射通道與非彈性散射通道的耦合,對最終散射態(tài)之間的躍遷有不可忽略的影響.

    表1 不同分子態(tài)對應(yīng)的解離極限下的原子態(tài)以及態(tài)能量的理論結(jié)果和實驗結(jié)果Table 1.Atomic states corresponding to different molecular states under dissociation limit,and the theoretical and experimental results of state energy.

    表1 中所有的散射通道都被考慮在QMOCC計算內(nèi),計算能區(qū)從 10-3— 103eV/u .圖2 和圖3分別為B 原子2s 和2p 軌道的電子轉(zhuǎn)移截面隨入射能量變化的函數(shù)圖像.在圖2 中,對分子的態(tài)的精細(xì)結(jié)構(gòu)也進(jìn)行區(qū)分,分別為3P 和1P 態(tài).因為1P態(tài)的散射通道能量較高,所以需要較大的入射能量,該通道才被打開.在圖3 中,將兩個2p 軌道的電子轉(zhuǎn)移截面加和,得到2p 軌道的電子轉(zhuǎn)移總截面.然而,從圖3 可以得到,2s 和2p 軌道的電子轉(zhuǎn)移截面在計算的能區(qū)范圍內(nèi),有至少兩個數(shù)量級的差異.所以內(nèi)殼層的2s 軌道的電子轉(zhuǎn)移在整個電子轉(zhuǎn)移過程中占主導(dǎo)地位.因為QMOCC 考慮各種電子態(tài)的散射通道之間的耦合,從圖3 可以看到,在低于1 eV/u 能量,電子轉(zhuǎn)移截面具有較多的 共振 現(xiàn)象.尤其在靠近 10-3eV/u 能 量附 近,2s 和2p 軌道的電子轉(zhuǎn)移截面都有強(qiáng)烈的共振,這是低能區(qū)碰撞過程量子效應(yīng)比較強(qiáng)的一種表現(xiàn).重粒子碰撞過程中,往往內(nèi)殼層的電子轉(zhuǎn)移都會伴隨退激發(fā)光輻射現(xiàn)象.

    圖2 質(zhì)子碰撞B 原子過程中2p 軌道電子轉(zhuǎn)移截面Fig.2.Electron transfer cross section of 2p orbitals in the process of protons colliding with B atom.

    圖3 質(zhì)子碰撞B 原子過程中2s 和2p 軌道電子轉(zhuǎn)移截面Fig.3.Electron transfer cross section of 2s and 2p orbitals in the process of protons colliding with B atom.

    對于質(zhì)子碰撞B 原子,如表1 所列,電子轉(zhuǎn)移通道包括 12Σ+,22Σ+,12Π+,42Σ+,32Π+,42Π+,52Σ+,62Σ+,82Σ+和 52Π+這10 個絕熱分子態(tài),其中涉及內(nèi)殼層2s 電子轉(zhuǎn)移的通道有6 個,導(dǎo)致外殼層2p 電子轉(zhuǎn)移的通道有4 個.彈性通道有32Σ+,22Π 兩個絕熱分子態(tài),其開通道能量為5.3 eV左右,所有其在低能區(qū)并未與電子轉(zhuǎn)移通道存在競爭關(guān)系.此外,還有3 個電子從2p 軌道躍遷到3s軌道的激發(fā)通道,分別是 72Σ+,92Σ+和 62Π+絕熱分子態(tài),激發(fā)通道都需要較大的能量才能打開,主要貢獻(xiàn)都在于與其他電子轉(zhuǎn)移通道耦合.圖4 和圖5為部分耦合較強(qiáng)的絕熱分子態(tài)之間的徑向耦合矩陣元和轉(zhuǎn)動耦合矩陣元隨核間距的變化.從圖4 可以看到,徑向耦合中,分子態(tài) 52Σ+-62Σ+之間的耦合最強(qiáng),這兩個通道都對2p 軌道的電子轉(zhuǎn)移有貢獻(xiàn).而且,除了分子態(tài) 22Σ+-32Σ+的耦合為2s軌道的電子轉(zhuǎn)移有影響,剩下幾個耦合比較強(qiáng)的徑向耦合矩陣元,比如 12Σ+-22Σ+和42Σ+-52Σ+也都對2p 軌道的電子轉(zhuǎn)移有影響.然而,2s 軌道的電子轉(zhuǎn)移涉及較多的絕熱分子態(tài)通道,所以其總體截面大于2p 軌道電子轉(zhuǎn)移.雖然,影響2p 軌道電子轉(zhuǎn)移的徑向耦合強(qiáng)度較大,但是圖3 的計算說明,這些耦合更多的造成2p 軌道電子轉(zhuǎn)移截面的減小.

    圖4 徑向耦合矩陣元隨核間距的變化Fig.4.Radial coupling matrix element changes with the nuclear distance.

    圖5 轉(zhuǎn)動耦合矩陣元隨核間距的變化Fig.5.Rotation coupling matrix element varies with the nuclear distance.

    同樣圖5 中最大的兩個轉(zhuǎn)動耦合矩陣元分別是分子態(tài) 42Π—52Π 和 52Π—62Π .其中42Π—52Π之間的轉(zhuǎn)動耦合也是影響2p 軌道電子轉(zhuǎn)移,但52Π—62Π 之間的耦合卻對 2s 軌道的電子轉(zhuǎn)移有貢獻(xiàn).42Π,52Π,62Π 三個態(tài)是較高的激發(fā)態(tài),因為要打開最低能量的態(tài) 42Π 也需要一定入射能量,所以轉(zhuǎn)動耦合雖然大,但是對電子轉(zhuǎn)移過程在低能區(qū)的影響極其有限.涉及2s 軌道的電子轉(zhuǎn)移通道,雖然耦合貢獻(xiàn)不及2p 軌道電子,但是這些通道對應(yīng)的分子態(tài)大多處于較低能級,且數(shù)量較多,所以最終計算的2s 軌道的電子轉(zhuǎn)移截面遠(yuǎn)大于2p 軌道電子.

    在表2 中,考慮整個系統(tǒng)服從麥克斯韋-玻爾茲曼能量分布.計算的電子總的轉(zhuǎn)移速率系數(shù)α(T),計算的溫度點量級從101—107,這個溫度跨度范圍比較廣泛,可以反映宏觀系統(tǒng)電子轉(zhuǎn)移的物理性質(zhì),注意計算的電子轉(zhuǎn)移的速率將包含2s 和2p 軌道.從表2 還可以發(fā)現(xiàn),因為相對高能態(tài),低能態(tài)的布局對速率系數(shù)的影響較大,隨著溫度的升高,電子轉(zhuǎn)移總速率系數(shù)增長幅度越來越平緩.當(dāng)溫度高于105K,系統(tǒng)的電子轉(zhuǎn)移總速率變化幅度明顯小于低溫區(qū).雖然我們考慮總的電子轉(zhuǎn)移速率,但是2s 軌道電子轉(zhuǎn)移占絕對主導(dǎo)地位,所以總的電子轉(zhuǎn)移速率,也可以在一定程度上反映內(nèi)殼層電子轉(zhuǎn)移后導(dǎo)致的X 射線輻射的強(qiáng)度.

    表2 電子轉(zhuǎn)移總的速率系數(shù),其中包含2s 和2p 軌道的電子轉(zhuǎn)移Table 2.Total rate coefficient of electron transfer,including electron transfer of 2s and 2p orbitals.

    4 結(jié)論

    采用QMCOCC 全量子的理論方法,研究了質(zhì)子碰撞B 原子電荷轉(zhuǎn)移過程.通過MRDCI 的第一性原理計算,得到了耦合分子9 個2Σ+態(tài)和6 個2Π 態(tài),取得了與實驗測量較好的結(jié)果.得到的分子態(tài)勢能曲線可避免交叉,這些能量點的避免交叉對絕熱分子態(tài)耦合具有重要的意義.本文所計算的15 個散射通道中,不僅包含彈性散射通道,還有電子激發(fā)和轉(zhuǎn)移的非彈性散射通道.其中,電子轉(zhuǎn)移通道都集中在中低能區(qū),而電子激發(fā)處于較高能區(qū).經(jīng)過計算,得到了從 10-3— 103eV/u 的質(zhì)子入射能區(qū)的電子轉(zhuǎn)移截面,并發(fā)現(xiàn)在低能區(qū)出現(xiàn)較強(qiáng)的共振現(xiàn)象,這主要是因為低能區(qū)通道耦合造成的散射量子效應(yīng).電子轉(zhuǎn)移截面分為2s 和2p 軌道的電子轉(zhuǎn)移截面,其中2s 軌道的電子轉(zhuǎn)移截面占據(jù)絕對主導(dǎo)地位.這些內(nèi)殼層的電子轉(zhuǎn)移過程,往往伴隨這高能量的輻射退激發(fā)發(fā)光想象.從15 個通道分析,發(fā)現(xiàn)包含2s 軌道的電子轉(zhuǎn)移的絕熱分子態(tài)數(shù)量最多,且這些通道的所需要的開通道能量也比較小.此外,還計算了大能量跨度的電子轉(zhuǎn)移速率系數(shù),為等離子體模擬提供重要的原子參數(shù).對于宏觀大尺度系統(tǒng)的電子轉(zhuǎn)移過程,低能區(qū)的散射截面占比較大.通過計算也驗證了這點,研究發(fā)現(xiàn)電子轉(zhuǎn)移速率隨著系統(tǒng)溫度的升高,增長的幅度越來越小.眾所周知,B 元素在天體中的豐度也較高,所以第一性原理計算不僅可以提供精確的原子參數(shù),也可以幫助人們深入認(rèn)識質(zhì)子碰撞B 原子的物理機(jī)制.

    感謝西安建筑科技大學(xué)理學(xué)院和北京應(yīng)用物理與計算數(shù)學(xué)研究所對本工作的支持.

    猜你喜歡
    質(zhì)子組態(tài)原子
    原子究竟有多小?
    原子可以結(jié)合嗎?
    帶你認(rèn)識原子
    基于PLC及組態(tài)技術(shù)的恒溫控制系統(tǒng)開發(fā)探討
    質(zhì)子束放療在腫瘤中的研究新進(jìn)展
    基于PLC和組態(tài)的智能電動擰緊系統(tǒng)
    淺談質(zhì)子守恒
    PLC組態(tài)控制在水箱控制系統(tǒng)的應(yīng)用
    鑄造行業(yè)三維組態(tài)軟件的應(yīng)用
    “質(zhì)子”號一箭發(fā)雙星
    太空探索(2014年6期)2014-07-10 13:06:11
    99国产综合亚洲精品| 99久久99久久久精品蜜桃| 精品不卡国产一区二区三区| 久久精品aⅴ一区二区三区四区| 国产精品爽爽va在线观看网站| 国产精品九九99| 性色av乱码一区二区三区2| cao死你这个sao货| 操出白浆在线播放| 亚洲 欧美一区二区三区| 欧美不卡视频在线免费观看| 男女下面进入的视频免费午夜| 偷拍熟女少妇极品色| 最近在线观看免费完整版| 中文字幕最新亚洲高清| 国产亚洲精品久久久久久毛片| 亚洲国产看品久久| 亚洲av第一区精品v没综合| x7x7x7水蜜桃| 久久久久免费精品人妻一区二区| 欧美又色又爽又黄视频| 成年女人永久免费观看视频| 亚洲人成电影免费在线| 麻豆一二三区av精品| 999精品在线视频| 午夜免费观看网址| 午夜免费激情av| 亚洲国产精品久久男人天堂| 九九在线视频观看精品| 88av欧美| 国产真实乱freesex| 国产精品一区二区免费欧美| 999精品在线视频| 在线免费观看的www视频| 一个人免费在线观看电影 | 一个人看的www免费观看视频| 久久国产乱子伦精品免费另类| 淫秽高清视频在线观看| 欧美色欧美亚洲另类二区| 美女被艹到高潮喷水动态| 999精品在线视频| 亚洲18禁久久av| 俺也久久电影网| 久久伊人香网站| 国内精品一区二区在线观看| 国产高清激情床上av| 国产成人精品久久二区二区91| 精品午夜福利视频在线观看一区| 超碰成人久久| 国产不卡一卡二| 亚洲精品乱码久久久v下载方式 | 日韩 欧美 亚洲 中文字幕| 俄罗斯特黄特色一大片| 中出人妻视频一区二区| 老熟妇乱子伦视频在线观看| 亚洲国产看品久久| 国模一区二区三区四区视频 | 一二三四社区在线视频社区8| 一本一本综合久久| 亚洲 欧美一区二区三区| 国产高清视频在线播放一区| 国产精品爽爽va在线观看网站| 99热6这里只有精品| 国产伦精品一区二区三区四那| 每晚都被弄得嗷嗷叫到高潮| 综合色av麻豆| 成年女人毛片免费观看观看9| 无限看片的www在线观看| 国产野战对白在线观看| 国产成人精品无人区| 精品无人区乱码1区二区| 黑人操中国人逼视频| 欧美又色又爽又黄视频| 免费看美女性在线毛片视频| 夜夜夜夜夜久久久久| 欧美日本视频| 国产精品久久久av美女十八| 亚洲五月天丁香| 成人18禁在线播放| 午夜视频精品福利| 中文字幕熟女人妻在线| 欧美黄色淫秽网站| 亚洲欧美日韩东京热| 精华霜和精华液先用哪个| 亚洲七黄色美女视频| 特大巨黑吊av在线直播| 国产视频一区二区在线看| 香蕉av资源在线| 99riav亚洲国产免费| 最新美女视频免费是黄的| 熟女电影av网| 国产高清视频在线观看网站| 亚洲人成伊人成综合网2020| 亚洲国产日韩欧美精品在线观看 | 制服丝袜大香蕉在线| 国产单亲对白刺激| 又紧又爽又黄一区二区| 99在线人妻在线中文字幕| 欧美日韩一级在线毛片| 搞女人的毛片| 午夜免费观看网址| 女警被强在线播放| 在线观看舔阴道视频| 亚洲美女黄片视频| 国产高清有码在线观看视频| 国内久久婷婷六月综合欲色啪| 在线观看66精品国产| 国产91精品成人一区二区三区| a在线观看视频网站| 国产高清三级在线| 神马国产精品三级电影在线观看| 亚洲九九香蕉| 黄色丝袜av网址大全| 亚洲成人免费电影在线观看| 久久久久久久久久黄片| 日本黄色片子视频| 国产亚洲精品久久久久久毛片| 一二三四社区在线视频社区8| 国产成人av教育| 小蜜桃在线观看免费完整版高清| 国产又色又爽无遮挡免费看| 午夜精品在线福利| 搡老熟女国产l中国老女人| 国产成人系列免费观看| 桃红色精品国产亚洲av| 国产亚洲精品综合一区在线观看| 成人三级黄色视频| 国产精华一区二区三区| 中文资源天堂在线| 在线永久观看黄色视频| 国产欧美日韩一区二区精品| 国产成+人综合+亚洲专区| 欧美在线黄色| 亚洲色图 男人天堂 中文字幕| 国产伦精品一区二区三区四那| 村上凉子中文字幕在线| 日韩高清综合在线| 99国产精品99久久久久| 男女下面进入的视频免费午夜| 午夜福利免费观看在线| 无人区码免费观看不卡| 精品久久久久久久人妻蜜臀av| 欧美日韩精品网址| 久久精品影院6| 非洲黑人性xxxx精品又粗又长| 国产伦在线观看视频一区| 国产主播在线观看一区二区| 国产成人一区二区三区免费视频网站| 老司机在亚洲福利影院| 久久久国产精品麻豆| 欧美一级毛片孕妇| 九九在线视频观看精品| 欧美日韩乱码在线| 色尼玛亚洲综合影院| 免费在线观看影片大全网站| 搡老妇女老女人老熟妇| 国内精品久久久久精免费| 免费看a级黄色片| 两个人视频免费观看高清| 观看免费一级毛片| 色在线成人网| 国产亚洲av嫩草精品影院| 99国产精品99久久久久| 国产私拍福利视频在线观看| 两人在一起打扑克的视频| 岛国在线观看网站| 亚洲,欧美精品.| 可以在线观看的亚洲视频| 国产午夜精品论理片| 国产蜜桃级精品一区二区三区| 久久久色成人| 搡老熟女国产l中国老女人| 国产伦精品一区二区三区视频9 | 日本三级黄在线观看| 国产极品精品免费视频能看的| 女生性感内裤真人,穿戴方法视频| 最新在线观看一区二区三区| 亚洲电影在线观看av| 久久久国产精品麻豆| 成熟少妇高潮喷水视频| 亚洲精品在线美女| 两个人视频免费观看高清| 亚洲在线观看片| 99国产极品粉嫩在线观看| 欧美激情在线99| 怎么达到女性高潮| 九九久久精品国产亚洲av麻豆 | 又紧又爽又黄一区二区| 午夜亚洲福利在线播放| 国产探花在线观看一区二区| 少妇丰满av| 五月伊人婷婷丁香| 久久久色成人| 婷婷精品国产亚洲av| 变态另类丝袜制服| 色尼玛亚洲综合影院| 亚洲 国产 在线| 国产成人精品无人区| 精品国内亚洲2022精品成人| 国产精品久久久久久精品电影| 欧美日韩一级在线毛片| 久久香蕉精品热| 亚洲性夜色夜夜综合| 日本在线视频免费播放| 午夜影院日韩av| www.www免费av| 亚洲国产日韩欧美精品在线观看 | 18禁裸乳无遮挡免费网站照片| 麻豆一二三区av精品| 亚洲片人在线观看| aaaaa片日本免费| 国产亚洲av高清不卡| 听说在线观看完整版免费高清| 亚洲精品乱码久久久v下载方式 | 亚洲欧美一区二区三区黑人| 久久久久久久久中文| 亚洲成人中文字幕在线播放| 国产成人av教育| 一级黄色大片毛片| 一个人免费在线观看的高清视频| 日本一本二区三区精品| 日韩 欧美 亚洲 中文字幕| 99精品久久久久人妻精品| 老熟妇乱子伦视频在线观看| 91av网站免费观看| 国产激情偷乱视频一区二区| 丰满的人妻完整版| 国产97色在线日韩免费| 人妻夜夜爽99麻豆av| 亚洲av第一区精品v没综合| 国产探花在线观看一区二区| 香蕉av资源在线| 日本黄色片子视频| 999久久久精品免费观看国产| 51午夜福利影视在线观看| 精品国内亚洲2022精品成人| 国产精品,欧美在线| www国产在线视频色| 亚洲中文字幕一区二区三区有码在线看 | 在线观看免费午夜福利视频| 变态另类成人亚洲欧美熟女| 午夜亚洲福利在线播放| 欧美高清成人免费视频www| 免费av毛片视频| 国产97色在线日韩免费| 日韩欧美国产一区二区入口| 男人和女人高潮做爰伦理| 在线十欧美十亚洲十日本专区| www.999成人在线观看| 99久久综合精品五月天人人| 国产三级在线视频| 国产成人一区二区三区免费视频网站| 12—13女人毛片做爰片一| 日韩av在线大香蕉| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品美女久久久久99蜜臀| 99久国产av精品| 国产高潮美女av| 最近最新免费中文字幕在线| 免费看光身美女| 床上黄色一级片| 国产久久久一区二区三区| 久久久久亚洲av毛片大全| 两性夫妻黄色片| 欧美日韩瑟瑟在线播放| 国产伦人伦偷精品视频| 亚洲av片天天在线观看| 九色成人免费人妻av| 国产成人啪精品午夜网站| 国产精品久久久久久亚洲av鲁大| 久99久视频精品免费| 在线观看66精品国产| 国产黄色小视频在线观看| 精品国内亚洲2022精品成人| 亚洲18禁久久av| 久久伊人香网站| 亚洲国产精品999在线| 国产真实乱freesex| 小说图片视频综合网站| 欧美成人性av电影在线观看| 国产午夜精品论理片| 19禁男女啪啪无遮挡网站| 国内久久婷婷六月综合欲色啪| 91av网一区二区| 操出白浆在线播放| 老熟妇仑乱视频hdxx| 嫁个100分男人电影在线观看| 免费在线观看视频国产中文字幕亚洲| 琪琪午夜伦伦电影理论片6080| 国产av一区在线观看免费| 巨乳人妻的诱惑在线观看| 黄色片一级片一级黄色片| 全区人妻精品视频| 午夜免费激情av| 网址你懂的国产日韩在线| 制服丝袜大香蕉在线| 少妇的逼水好多| 深夜精品福利| 91老司机精品| 黄片小视频在线播放| 欧美成人免费av一区二区三区| 成人国产一区最新在线观看| 国产麻豆成人av免费视频| 一区二区三区激情视频| 丝袜人妻中文字幕| av福利片在线观看| 男女视频在线观看网站免费| 亚洲自偷自拍图片 自拍| 亚洲国产欧美人成| 精品久久久久久久人妻蜜臀av| 亚洲精品美女久久久久99蜜臀| av黄色大香蕉| 久久久成人免费电影| 精品久久蜜臀av无| 无遮挡黄片免费观看| 国产麻豆成人av免费视频| 我的老师免费观看完整版| 在线观看午夜福利视频| 九九在线视频观看精品| 一本久久中文字幕| 熟女电影av网| 给我免费播放毛片高清在线观看| 久久香蕉精品热| 国产精品久久视频播放| 在线观看免费视频日本深夜| 久久精品91无色码中文字幕| 国产精品亚洲一级av第二区| 黄片大片在线免费观看| 久久香蕉精品热| 午夜福利在线观看吧| 少妇熟女aⅴ在线视频| 极品教师在线免费播放| 国产精品99久久久久久久久| 国产黄色小视频在线观看| 好看av亚洲va欧美ⅴa在| 国产精品爽爽va在线观看网站| 99久久精品一区二区三区| 亚洲熟妇熟女久久| 免费一级毛片在线播放高清视频| 亚洲精品国产精品久久久不卡| 久久久久久久久免费视频了| 1024香蕉在线观看| 亚洲天堂国产精品一区在线| 欧美黄色片欧美黄色片| 两个人看的免费小视频| 97超视频在线观看视频| 99久久成人亚洲精品观看| 国产精品香港三级国产av潘金莲| 欧美极品一区二区三区四区| 国产一区二区在线观看日韩 | 久久久久久久久免费视频了| 97超视频在线观看视频| 一级毛片女人18水好多| 亚洲中文av在线| 脱女人内裤的视频| 久久久水蜜桃国产精品网| 日本一本二区三区精品| 国产av不卡久久| 免费观看精品视频网站| 欧美一区二区国产精品久久精品| 99久久99久久久精品蜜桃| 一个人免费在线观看的高清视频| 好看av亚洲va欧美ⅴa在| 久久中文字幕人妻熟女| 亚洲专区字幕在线| 亚洲欧美一区二区三区黑人| 亚洲色图av天堂| 免费搜索国产男女视频| 好男人电影高清在线观看| 99精品久久久久人妻精品| 一区二区三区高清视频在线| tocl精华| 国产三级中文精品| 日本黄色片子视频| 午夜福利欧美成人| av黄色大香蕉| 一二三四社区在线视频社区8| 人妻丰满熟妇av一区二区三区| 美女高潮喷水抽搐中文字幕| 夜夜躁狠狠躁天天躁| 日韩欧美三级三区| 一级毛片精品| av在线天堂中文字幕| 国产成+人综合+亚洲专区| 国产伦人伦偷精品视频| 网址你懂的国产日韩在线| 亚洲五月天丁香| 黑人巨大精品欧美一区二区mp4| 久久精品91无色码中文字幕| 国产亚洲av高清不卡| 最近最新免费中文字幕在线| 欧美日韩综合久久久久久 | www日本在线高清视频| 99国产精品一区二区蜜桃av| 国产精品98久久久久久宅男小说| 午夜福利视频1000在线观看| 亚洲成人免费电影在线观看| 午夜精品久久久久久毛片777| 免费在线观看影片大全网站| 国内精品久久久久久久电影| 搡老岳熟女国产| 99久国产av精品| 日本 av在线| av黄色大香蕉| 亚洲熟妇熟女久久| 韩国av一区二区三区四区| 亚洲性夜色夜夜综合| 特级一级黄色大片| 中亚洲国语对白在线视频| 日韩中文字幕欧美一区二区| 操出白浆在线播放| 亚洲自拍偷在线| www国产在线视频色| 国产又色又爽无遮挡免费看| 婷婷丁香在线五月| 黄色女人牲交| 99热这里只有是精品50| 久久精品影院6| 中文字幕av在线有码专区| 叶爱在线成人免费视频播放| 免费无遮挡裸体视频| 18禁裸乳无遮挡免费网站照片| 国内精品美女久久久久久| 亚洲国产色片| 午夜免费成人在线视频| 一卡2卡三卡四卡精品乱码亚洲| 色综合亚洲欧美另类图片| 日本一本二区三区精品| 综合色av麻豆| 日韩精品中文字幕看吧| 这个男人来自地球电影免费观看| 久久久久久人人人人人| 日本撒尿小便嘘嘘汇集6| 精品久久久久久,| 成人三级做爰电影| 在线免费观看的www视频| 成人一区二区视频在线观看| 成在线人永久免费视频| 岛国在线观看网站| 久久九九热精品免费| 成人av一区二区三区在线看| 日本黄色片子视频| 亚洲中文日韩欧美视频| 久久中文字幕一级| 国产精品一区二区精品视频观看| 亚洲五月婷婷丁香| 此物有八面人人有两片| www国产在线视频色| 久久精品国产亚洲av香蕉五月| 99热这里只有是精品50| 成人特级av手机在线观看| 久久久久免费精品人妻一区二区| 亚洲国产欧美一区二区综合| 午夜福利在线观看免费完整高清在 | 香蕉av资源在线| 国产探花在线观看一区二区| 免费看美女性在线毛片视频| 久久国产精品影院| 国产成人aa在线观看| 精品久久久久久,| 久久伊人香网站| 99热只有精品国产| 中文字幕高清在线视频| 91久久精品国产一区二区成人 | 亚洲av成人不卡在线观看播放网| 免费看日本二区| 久久精品aⅴ一区二区三区四区| 国产在线精品亚洲第一网站| 午夜福利免费观看在线| 亚洲av五月六月丁香网| 久9热在线精品视频| 最好的美女福利视频网| 99久久成人亚洲精品观看| 国产午夜精品久久久久久| 韩国av一区二区三区四区| 老鸭窝网址在线观看| 三级国产精品欧美在线观看 | 久久精品aⅴ一区二区三区四区| 亚洲在线自拍视频| 免费看美女性在线毛片视频| 精品不卡国产一区二区三区| 国产成人一区二区三区免费视频网站| 一a级毛片在线观看| 亚洲在线自拍视频| 日韩欧美免费精品| 又紧又爽又黄一区二区| 美女 人体艺术 gogo| 少妇丰满av| 久久久精品大字幕| 久久中文字幕人妻熟女| 国产精品亚洲一级av第二区| 精品一区二区三区四区五区乱码| 亚洲美女视频黄频| 国产激情偷乱视频一区二区| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品在线美女| 欧美丝袜亚洲另类 | 成人欧美大片| 欧美成人性av电影在线观看| 色哟哟哟哟哟哟| 成人鲁丝片一二三区免费| 久久精品国产清高在天天线| 99久久99久久久精品蜜桃| 久久亚洲真实| 日韩欧美 国产精品| 中出人妻视频一区二区| 天堂动漫精品| 日本熟妇午夜| 美女cb高潮喷水在线观看 | 国产精品永久免费网站| 床上黄色一级片| 欧美日韩瑟瑟在线播放| 香蕉国产在线看| 亚洲成人精品中文字幕电影| 伦理电影免费视频| 欧美日韩中文字幕国产精品一区二区三区| 成人国产综合亚洲| 成人av在线播放网站| 日韩 欧美 亚洲 中文字幕| 日韩大尺度精品在线看网址| 日本免费a在线| 精品国产亚洲在线| 又爽又黄无遮挡网站| 日韩欧美三级三区| 欧美乱色亚洲激情| 中国美女看黄片| 亚洲av成人不卡在线观看播放网| 国内久久婷婷六月综合欲色啪| 欧美一区二区国产精品久久精品| 国产精品99久久99久久久不卡| 国产精品美女特级片免费视频播放器 | 午夜精品久久久久久毛片777| 精华霜和精华液先用哪个| 观看免费一级毛片| 天天添夜夜摸| 日韩中文字幕欧美一区二区| 国产成人欧美在线观看| 国产精品一区二区三区四区久久| 91在线观看av| 黄片大片在线免费观看| 亚洲欧美日韩无卡精品| 日韩大尺度精品在线看网址| 99国产精品一区二区蜜桃av| 国产乱人视频| 欧洲精品卡2卡3卡4卡5卡区| 一个人免费在线观看的高清视频| 黄色成人免费大全| 在线免费观看不下载黄p国产 | 99热精品在线国产| 国产精品久久久久久人妻精品电影| 一区二区三区高清视频在线| 中出人妻视频一区二区| 亚洲国产欧洲综合997久久,| www日本在线高清视频| 成人av在线播放网站| 日本 av在线| 在线国产一区二区在线| 欧美xxxx黑人xx丫x性爽| 好男人电影高清在线观看| 国内久久婷婷六月综合欲色啪| 久久久精品欧美日韩精品| 中文亚洲av片在线观看爽| 女人被狂操c到高潮| 免费高清视频大片| 成年版毛片免费区| x7x7x7水蜜桃| 日韩欧美三级三区| 亚洲国产精品合色在线| 别揉我奶头~嗯~啊~动态视频| 一区福利在线观看| 一级毛片高清免费大全| 日韩高清综合在线| 91av网一区二区| 国产伦精品一区二区三区四那| 色精品久久人妻99蜜桃| 2021天堂中文幕一二区在线观| 精品久久久久久久久久免费视频| 最近最新中文字幕大全免费视频| 国产三级中文精品| 一本精品99久久精品77| 免费看美女性在线毛片视频| 1000部很黄的大片| 色尼玛亚洲综合影院| 免费大片18禁| 国内精品美女久久久久久| 国产精品久久电影中文字幕| 噜噜噜噜噜久久久久久91| 久久精品91蜜桃| 亚洲av电影不卡..在线观看| 亚洲av电影在线进入| 亚洲人成网站在线播放欧美日韩| 久久草成人影院| 日韩三级视频一区二区三区| 网址你懂的国产日韩在线| 桃色一区二区三区在线观看| 国产欧美日韩精品亚洲av| 曰老女人黄片| 国产视频内射| 男女床上黄色一级片免费看| 成年版毛片免费区| 桃色一区二区三区在线观看| 亚洲精品美女久久av网站| 手机成人av网站| 桃色一区二区三区在线观看| АⅤ资源中文在线天堂| 变态另类成人亚洲欧美熟女| 少妇的丰满在线观看| АⅤ资源中文在线天堂| 三级男女做爰猛烈吃奶摸视频| 桃色一区二区三区在线观看| 99在线人妻在线中文字幕| 日本成人三级电影网站| 99久久综合精品五月天人人| 国产精品爽爽va在线观看网站| 欧美日韩黄片免|