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

    B2+和B+離子的靜態(tài)偶極極化率和超極化率的理論研究*

    2023-07-27 10:59:12陳池婷吳磊王霞王婷劉延君蔣軍董晨鐘
    物理學(xué)報 2023年14期
    關(guān)鍵詞:電偶基態(tài)振子

    陳池婷 吳磊 王霞 王婷 劉延君 蔣軍? 董晨鐘

    1) (西北師范大學(xué)物理與電子工程學(xué)院,蘭州 730070)

    2) (蘭州理工大學(xué)理學(xué)院,蘭州 730050)

    利用相對論組態(tài)相互作用模型勢方法計算了B2+和B+離子的波函數(shù)、能級和振子強度,進一步得到B2+離子2s1/2,2p1/2,2p3/2,3s1/2 態(tài)的電偶極極化率和基態(tài)2s1/2 的超極化率,以及B+離子2s21S0 與2s2p 3P0 態(tài)的電偶極極化率.B2+離子2p1/2 和2p3/2 的偶極極化率為負(fù)值,基態(tài)2s1/2 的超極化率的貢獻主要來自于與極化率相關(guān)的 α 1β0 項.對于B+離子,鐘躍遷2s2p 3P0 → 2s21S0 在室溫下的黑體輻射頻移是0.01605 Hz,該黑體輻射頻移比堿土金屬原子的鐘躍遷黑體輻射頻移小1—2 個數(shù)量級.

    1 引言

    近年來,隨著激光冷卻和囚禁技術(shù)的快速發(fā)展,光學(xué)原子鐘的精度和穩(wěn)定性得到極大的提高[1-7].高精度的原子光鐘可用于精確測量基本物理常數(shù)[8],測試局部洛倫茲不變性[9,10],探索精細(xì)結(jié)構(gòu)常數(shù)α隨時間的變化[11,12],探測暗物質(zhì)和暗能量[13,14],探測引力波[15]以及檢測超越粒子物理標(biāo)準(zhǔn)模型的新作用力[14,15].

    在外電場中,原子的能級會發(fā)生斯塔克(Stark)效應(yīng),靜態(tài)Stark 效應(yīng)導(dǎo)致的能級移動可以寫為[16]

    其中α1為靜態(tài)電偶極極化率,FZ是靜電場沿Z軸的分量,γ0為靜態(tài)電偶極超極化率.電偶極極化率和超極化率都是描述外電場中原子的電子云偏離正常分布程度的物理量.極化率是電荷分布對外電場的最低階響應(yīng),極化率在原子鐘[8,17]的研制方面有著十分重要的應(yīng)用.例如,利用原子的極化率可以估算原子鐘鐘躍遷的黑體輻射(BBR)頻移,BBR 頻移和Stark 效應(yīng)是影響原子鐘精度的主要來源[1,3,18-20].極化率可用來構(gòu)建超冷原子之間的長程相互作用勢[21-23]等.原子的超極化率可以評估高階Stark 頻移,例如,超極化率對Sr 和Yb 原子鐘可以產(chǎn)生0.1—1.0 mHz 的Stark 頻移[24-26],對Al+離子量子邏輯鐘的影響可以達到10—19量級[1].超極化率已經(jīng)成為影響新一代原子鐘精度的重要因素.然而實驗上精確測量原子極化率和超極化率非常困難.因此需要理論方面提供精確的原子極化率以及超極化率.

    B2+離子是三電子體系,基態(tài)為1s22s2S1/2,B2+離子是精密測量物理非常重要的研究體系[27,28].B+離子的基態(tài)為2s21S0,第一激發(fā)態(tài)2s2p3P0是一個長壽命的亞穩(wěn)態(tài),2s21S0→2s2p3P0是開發(fā)高精度光鐘的潛在候選體系之一[29].因此,關(guān)于B2+離子和B+離子的原子結(jié)構(gòu)參數(shù)(例如能級、振子強度、電偶極極化率以及超極化率)的精確計算對于精密測量物理具有十分重要的意義.

    本文采用相對論組態(tài)相互作用模型勢方法(記為RCICP),計算了B2+和B+離子的能級和振子強度,進一步利用求和規(guī)則得到了B2+離子2s1/2,2p1/2,2p3/2,3s1/2態(tài)的電偶極極化率、基態(tài)2s1/2的超極化率,以及B+離子2s21S0和2s2p3P0態(tài)的電偶極極化率,得到鐘躍遷2s2p3P0→ 2s21S0的黑體輻射頻移,本工作結(jié)果與已有的理論結(jié)果進行詳細(xì)的比較.

    2 理論方法

    RCICP 方法的核心思想是將原子體系簡化為原子實部分和價電子部分.對于B2+離子,原子實是1s2,價電子是2s.對于B+離子,價電子是2s2.原子實的軌道波函數(shù)通過求解Dirac-Fock 方程獲得.然后,用半經(jīng)驗極化勢方法計算了價電子波函數(shù).兩個價電子的有效哈密頓量可以寫成[27]

    其中Z是原子序數(shù).Vdir(r)和Vexc(r) 分別表示價電子與原子實電子之間的直接相互作用和交換相互作用.單價電子極化勢Vp1(r) 可以寫為

    3 結(jié)果與討論

    3.1 B2+離子

    3.1.1 能級

    表2 列出了利用RCICP 方法計算的B2+離子基態(tài)和部分低激發(fā)態(tài)相對于原子實的能級和NIST 推薦的實驗值[33].在RCICP 方法中,通過調(diào)節(jié)表1 中的截斷參數(shù)ρl,j,使2s,2pj和3dj能級都非常接近于NIST 值.從表2 也可以看出,ns(n≥3),npj(n≥3),ndj(n≥4)和nfj(n≥4) 等 更高激發(fā) 態(tài)能級值與實驗值也都符合得非常好,僅在小數(shù)點后第5 位有差別.例如,4 p1/2和4 p3/2的結(jié)果分別為—0.2874707 a.u.和—0.2874514 a.u.,這與實驗結(jié)果—0.2875098 a.u.和—0.2874920 a.u.符合得很好.

    表1 B2+離子的截斷參數(shù) ρ l,j (單位: a.u.)Table 1. Cut-off parameters ρ l,j of B2+ ions (in a.u.).

    表2 B2+離子的基態(tài)和部分低激發(fā)態(tài)相對于原子實的能級,實驗值(Expt.) [33]是來自于NIST 的數(shù)據(jù)(單位: a.u.),“Diff.”表示用RCICP 方法計算的結(jié)果與NIST 結(jié)果之差的百分比Table 2. Energy levels of the ground state and some low-lying states of B2+ ions relative to atomic core.Experimental values (Expt.) [33] are from the NIST data (in a.u.).“Diff.” denotes the difference in percentage from calculated by RCICP method and NIST results.

    3.1.2 振子強度

    從態(tài)i到態(tài)m電偶極躍遷的吸收振子強度可以表示為[34-37]

    式中n,l,j分別為主量子數(shù)、軌道角動量量子數(shù)和總角動量量子數(shù),ΔEm→i=Em-Ei表示躍遷能.D是電偶極躍遷算符,可以表示為[38-40]

    表3 列出了利用RCICP 方法計算的B2+離子基態(tài)和部分激發(fā)態(tài)之間的振子強度并與相應(yīng)NIST 的值[41]、相對論多體微擾(RMBPT)方法[42]以及Hylleraas 方法(HR)[43]的計算結(jié)果進行了比較.從表3 可以看出,對于2s1/2→3p1/2,2p1/2→3d3/2,2p3/2→3s1/2,2p3/2→3d3/2,5/2,3s1/2→3p1/2,3/2的躍遷,RCICP 計算的結(jié)果與NIST 值[41]以及RMBPT方法[42]的結(jié)果符合得非常好,差別小于0.1%.對于2s1/2→2p1/2,3/2,2s1/2→3p3/2,2p1/2→3s1/2,2p1/2→4s1/2,2p1/2→4d3/2,2p3/2→4s1/2,2p3/2→4d3/2,5/2,3s1/2→4p1/2,3/2的躍遷,本工作計算的振子強度與NIST 值[41]、RMBPT 方法[42]以及HR 方法[43]的計算結(jié)果差別在0.6%以內(nèi).

    表3 B2+離子基態(tài)和部分低激發(fā)態(tài)之間躍遷的振子強度,“Diff.”表示用RCICP 方法計算的結(jié)果與NIST 結(jié)果[41]之差的百分比Table 3. Oscillator strengths of transitions between the ground state and some low-lying states of B2+ ions.“Diff.” represents the difference in percentage form calculated by RCICP method and NIST results.

    3.1.3 極化率

    原子態(tài)i的靜態(tài)電偶極極化率可以表示為

    其中角動量ji=1/2 的原子態(tài)的電偶極極化率與磁投影Mi無關(guān),而對于ji >1/2 的原子態(tài),它取決于Mi,即有標(biāo)量極化率(αS)和張量極化率(αT).標(biāo)量極化率和張量極化率通常定義為所有可能中間態(tài)的總和,而且包括連續(xù)態(tài),如下式所示:

    表4 列出了采用RCICP 方法計算的B2+離子基態(tài)2s1/2,以及部分低激發(fā)態(tài)2p1/2,2p3/2和3s1/2的靜態(tài)電偶極極化率和占主導(dǎo)躍遷的貢獻,并與其他理論方法計算結(jié)果進行比較.表中的“Remains”表示高激發(fā)態(tài)以及連續(xù)態(tài)的貢獻,“Core”表示原子實(1s2)電子的貢獻,“Total”表示總的極化率.對于基態(tài)2s1/2的靜態(tài)電偶極極化率,2s1/2→2p1/2,3/2躍遷占主導(dǎo),RCICP 計算的結(jié)果與full-core plus correlation (FCPC)[44],非相對論組態(tài)相互作用加核勢(CICP)[45],關(guān)聯(lián)組態(tài)(SCC)[46]以及完全關(guān)聯(lián)高斯(FCG)[47]方法計算的結(jié)果符合得非常好,差別在0.1%以內(nèi).對于2p1/2態(tài),貢獻主要來自于2p1/2→2s1/2和2p1/2→3d3/2的躍遷,但是由于2p1/2→2s1/2吸收能為負(fù)值,這一項對2p1/2的貢獻為負(fù),因此2p1/2的電偶極極化率為負(fù)值,RCICP 計算的結(jié)果為—0.5737 a.u.對于2p3/2態(tài)的標(biāo)量極化率,主要貢獻來自于2p3/2→ 2s1/2和2p3/2→ 3d5/2的躍遷,RCICP 計算的結(jié)果為—0.5713 a.u.,和2p1/2態(tài)類似,由于2p3/2→2s1/2躍遷對2p3/2態(tài)極化率的貢獻為負(fù),所以2p3/2態(tài)極化率小于零.RCICP 計算的結(jié)果與CICP 方法[45]計算的結(jié)果符合得非常好,其差別在0.8%以內(nèi),這里需要指出的是CICP[45]是非相對論結(jié)果,其值僅僅是2p 態(tài)的,并不是精細(xì)能級的極化率.對于2p3/2態(tài)的張量極化率,其主要貢獻來自于2p3/2→2s1/2和2p3/2→3d5/2的躍遷,分別為2.4963 a.u.與—0.2537 a.u.,RCICP 計算的結(jié)果為2.1683 a.u.,目前還沒有可比較的相關(guān)數(shù)據(jù).對于3s1/2態(tài),主要貢獻來自于3s1/2→3p1/2,3/2的躍遷,本文計算結(jié)果為182.90 a.u.,與CICP[45]方法計算的結(jié)果182.94 a.u.符合得非常好,差別在0.02%以內(nèi).

    表4 B2+離子基態(tài)與部分低激發(fā)態(tài)的靜態(tài)電偶極標(biāo)量極化率與張量極化率以及主要躍遷的貢獻(單位: a.u.)Table 4. Static electric-dipole scalar and tensor polarizability of the ground state and some low-lying state of B2+ ions and breakdowns of the contributions of individual transitions (in a.u.).

    3.1.4 超極化率

    根據(jù)四階微擾理論,對于B2+離子基態(tài)2s1/2的超極化率可以簡化為[48]

    其中α1是基態(tài)的靜態(tài)電偶極極化率,β0是電偶極極化率的一階非絕熱修正[35].T(s1/2,pj′,s1/2,pj′′)和T(s1/2,pj′,dj,pj′′) 可以表示為

    表5 列出了B2+的基態(tài)2s1/2的超極化率和中間態(tài)的貢獻,并與Roy 和Bhattacharya[49]使用耦合Hatree-Fock (CHF)及 Drake 和Cohen[50]使用未耦合Hatree-Fock (UCHF) 近似方法得到的經(jīng)驗公式計算的結(jié)果進行了比較.根據(jù)(11)式,B2+離子基態(tài)超極化率的計算由3 部分組成,分別為T(s,pj′,s,pj′′),T(s,pj′,dj,pj′′)和α1β0.為了評 估該計算結(jié)果的誤差,將對超極化率有重要貢獻的2s→2pj,2pj→3dj躍遷矩陣元替換為 NIST 推薦的躍遷矩陣元,將3s,3pj,4s,4pj,4dj,5s,5pj和5dj的能級值也替換成NIST 推薦值,并重新計算了超極化率(標(biāo)記為γ0C(2s) ),研究發(fā)現(xiàn)γ0(2s) 與γ0C(2s)的結(jié)果僅僅相差 6.645 a.u.,這個差別僅占目前計算的超極化率的 0.6%.因此,采用 RCICP計算的超極化率的精度在 1%以內(nèi).在所有貢獻中,與極化率相關(guān)的α1β0項貢獻最大,結(jié)果 為134.364(586) a.u..RCICP 方法計算的B2+基態(tài)的超極化率結(jié)果為—1063.346(6.645) a.u.,與CHF[49]以及UCHF[50]近似方法計算的結(jié)果符合得比較好.

    表5 B2+離子基態(tài)的超極化率及其中間態(tài)對超極化率的貢獻(單位: a.u.)Table 5. Hyperpolarizability of the ground state of B2+ ion and the contributions to the hyperpolarizability (in a.u.).

    3.2 B+離子

    3.2.1 能 級

    利用計算得到的B2+的單電子軌道,耦合得到B+的組態(tài)波函數(shù),進一步全對角化了B+離子兩電子組態(tài)空間的哈密頓矩陣,得到了B+離子原子態(tài)的波函數(shù).表6 列出了基態(tài)以及部分低激發(fā)態(tài)相對于原子實 1 s21S0的能級,并與NIST 推薦的實驗值[51]進行比較.從表6 可以發(fā)現(xiàn),RCICP 方法計算的結(jié)果與NIST 推薦的實驗值符合得非常好,差異不超過0.05%.

    表6 B+基態(tài)和部分低激發(fā)態(tài)相對于原子實的能級值,實驗值(Expt.) [51]是來自于NIST 的數(shù)據(jù)(單位: a.u.),“Diff.”表示用RCICP 方法計算的結(jié)果與NIST 結(jié)果之差的百分比Table 6. Energy levels of the ground state and some low-lying states of B+ ions relative to atomic core.Experimental values (Expt.) are from the NIST data (in a.u.).“Diff.” denotes the difference in percentage from calculated by RCICP method and NIST results.

    3.2.2 振子強度

    基于以上的能級和波函數(shù),利用(6)式,得到了B+離子基態(tài)和部分低激發(fā)態(tài)之間的電偶極躍遷振子強度,表7 列出了本文計算結(jié)果,并與相應(yīng)的NIST[41]推薦值、CICP[45]方法、B-樣條基組態(tài)相互作用加核勢(BCICP)方法[52]、多組態(tài)Hartree-Fock 與Breit-Pauli (MCHF-BP)方法[53]、多組態(tài)Hartree-Fock (MCHF)方法[54]計算的結(jié)果進行比較.對于2s21S0→2s2p1P1共振躍遷,RCICP 的結(jié)果與NIST 推薦值[41]和其他理論結(jié)果[45,52-54]符合得非常好,差別在0.4%以內(nèi).從表7 還可以看出,除2s21S0→2s4p1P1和2s21S0→2s5p1P1躍遷外,RCICP 方法計算的其余躍遷的振子強度與NIST 推薦值[41]之間的差別在2%以內(nèi),2s21S0→2s4p1P1和2s21S0→2s5p1P1躍遷的振子強度與NIST 推薦值[41]的差別分別約為4%和7%.造成這種差異的原因可能有3 個,其中一個是,單電子軌道是利用有限個基矢的線性組合得到的,基態(tài)的收斂性非常好,單電子高激發(fā)態(tài)的收斂性要比基態(tài)差.另外一個原因是,在本文方法中根據(jù)原子態(tài)的對稱性(宇稱和總角動量J)將計算分為不同的組,計算2s2p,2s3p,2s4p,2s5p 態(tài)所使用的組態(tài)空間相同,2s4p1P1和2s5p1P1的能級相對于2s2p1P1態(tài)收斂性變差,相應(yīng)的2s4p1P1和2s5p1P1態(tài)到基態(tài)的躍遷振子強度與NIST 推薦值之間的相對差異變大.最后一個原因是,2s4p1P1和2s5p1P1態(tài)到基態(tài)的躍遷振子強度的值非常小.例如,對于2s21S0→2s5p1P1躍遷的振子強度,本文計算的結(jié)果是0.0224,而NIST 推薦值為0.0241.RCICP 的結(jié)果與其他理論方法[45,52-54]計算的結(jié)果也符合得比較好,除2s2p3P0→2p23P1的躍遷外,其余躍遷振子強度的差別在2%以內(nèi),2s2p3P0→2p23P1躍遷的振子強度與BCICP 結(jié)果[52]的差別小于7%.

    表7 B+離子基態(tài)和部分低激發(fā)態(tài)之間電偶極躍遷的振子強度(單位: a.u.)Table 7. Oscillator strengths of electric-dipole transitions between the ground state and some low-lying states of B+ ions(in a.u.).

    3.2.3 極化率

    表8 所示為B+離子與鐘躍遷相關(guān)的2s21S0和2s2p3P0態(tài)的靜態(tài)電偶極極化率,并與組態(tài)相互作用(CI)方法[55]、組態(tài)相互作用和多體微擾理論(CI+MBPT)方法[55]、組態(tài)相互作用加全階(CI+all-orders)方法[55]、有限場耦合簇(CCD+ST)方法[56]、CICP 方法[45]、微擾相對論耦合簇(PRCC)方法[29]、三重激發(fā)微擾的耦合簇(CCSDpT)方法[57]以及瑞利-里茲變分(RRV)方法[58]計算的結(jié)果進行比較.從表8 可以看出,對于2s21S0態(tài)的極化率,2s21S0→2s2p1P1的躍遷占主導(dǎo),貢獻約為93%,RCICP 的結(jié)果與CI+all-orders 方法[55]、CI[55]、CI+MBPT[55]、CCD+ST[56]、RRV[58]以及CICP[45]方法的結(jié)果差別不超過0.6%,與PRCC[29]結(jié)果差異在2%左右,與CCSDpT[57]結(jié)果的差異在7%之間.對于2s2p3P0態(tài)的極化率,2s2p3P0→2p23P1的躍遷和2s2p3P0→2s3d3D1的躍遷占主導(dǎo),貢獻分別約為56%和24%.目前RCICP 結(jié)果與CI[55],CI+MBPT[55],CI+all-order[55]以及CICP 方法[45]計算的結(jié)果差異不超過0.3%.

    表8 B+離子2s21S0 和2s2p 3P0 的電偶極極化率Table 8. Electric-dipole polarizability of 2s21S0 and 2s2p 3P0 states of B+ ions.

    與溫度T有關(guān)的BBR 頻移可以寫為

    式中,Δα1是躍遷初末態(tài)的靜態(tài)電偶極極化率之差,系數(shù)η是一個可以忽略的微小的動態(tài)修正[20].T是環(huán)境溫度(K),當(dāng)T為室溫300 K 時,鐘躍遷2s2p3P0→2s21S0的BBR 頻移為0.01605 Hz,與CI+MBPT[55]方法計算的結(jié)果0.0159(16) Hz 符合得非常好,該BBR 值比Ca+,Sr+離子鐘躍遷2D5/2→2S1/2的BBR 頻移(Ca+: 0.38(1) Hz[59],Sr+:0.250(9) Hz[60])小1 個數(shù)量級,比Ca,Sr,Yb 原子鐘的躍遷3P0→1S0的BBR 頻移 (Ca: 1.171 Hz,Sr: 2.354 Hz,Yb: 1.25 Hz)小2 個數(shù)量級[20].這一特性對B+離子鐘的研究具有重要意義.

    4 結(jié)論

    利用相對論組態(tài)相互作用模型勢方法計算了B2+和B+離子的波函數(shù)、能級和振子強度,進一步得到了B2+離子2s1/2,2p1/2,2p3/2,3s1/2態(tài)的電偶極極化率和基態(tài)2s1/2的超極化率,以及B+離子2s21S0和2s2p3P0態(tài)的電偶極極化率,并與其他理論結(jié)果進行比較,本文結(jié)果與已有結(jié)果符合得非常好.對于B2+離子2p1/2和2p3/2態(tài)的極化率為負(fù)值,是因為2p1/2,3/2→2s1/2躍遷的吸收能為負(fù)值,這項對2p1/2和2p3/2態(tài)的極化率的貢獻都為負(fù).對于2p3/2態(tài)的張量極化率,其主要貢獻來自于2p3/2→2s1/2和2p3/2→3d5/2的躍遷,分別為2.4963 a.u.和—0.2537 a.u.,RCICP 計算的結(jié)果為2.1683 a.u..基態(tài)2s1/2的超極化率的貢獻主要來自于與極化率相關(guān)的α1β0項.對于B+離子基態(tài)2s21S0的偶極極化率是9.6220 a.u.,2s2p3P0態(tài)的偶極極化率是7.7594 a.u.,計算的鐘躍遷2s2p3P0→ 2s21S0的黑體輻射頻移是0.01605 Hz,這個黑體輻射頻移比堿土金屬原子的鐘躍遷黑體輻射頻移小1—2 個數(shù)量級.

    猜你喜歡
    電偶基態(tài)振子
    更正說明
    腐蝕與防護(2022年6期)2023-01-10 11:25:26
    一類非線性Choquard方程基態(tài)解的存在性
    擬相對論薛定諤方程基態(tài)解的存在性與爆破行為
    一類反應(yīng)擴散方程的Nehari-Pankov型基態(tài)解
    非線性臨界Kirchhoff型問題的正基態(tài)解
    彈簧振子問題的分析與求解
    非線性Duffing擾動振子共振機制的研究
    Q245R/316L與Q345R/316L電偶對在3.5%NaCl溶液中的電偶腐蝕
    鋼質(zhì)接地材料電偶加速腐蝕評價方法
    功能材料(2016年8期)2016-09-12 06:22:54
    基于近似熵和混沌振子的電力諧波檢測與估計
    精品国产一区二区三区久久久樱花| 日本av手机在线免费观看| 一区在线观看完整版| 国产精品成人在线| 色精品久久人妻99蜜桃| 91av网站免费观看| 精品国产乱码久久久久久男人| 精品久久久久久久毛片微露脸 | 麻豆av在线久日| 老鸭窝网址在线观看| 成人亚洲精品一区在线观看| 国产一区二区激情短视频 | 午夜久久久在线观看| 婷婷丁香在线五月| 狂野欧美激情性xxxx| 高清av免费在线| 日韩制服骚丝袜av| 国产成人av教育| 在线观看免费午夜福利视频| av超薄肉色丝袜交足视频| 欧美在线一区亚洲| 男女边摸边吃奶| 电影成人av| tocl精华| 黄色片一级片一级黄色片| 亚洲天堂av无毛| 免费高清在线观看视频在线观看| 黄片大片在线免费观看| a在线观看视频网站| 成人手机av| 国产av又大| 国产av又大| 国产1区2区3区精品| 多毛熟女@视频| 欧美黄色淫秽网站| 成人国产av品久久久| 99re6热这里在线精品视频| 日韩欧美一区视频在线观看| 久久精品aⅴ一区二区三区四区| 他把我摸到了高潮在线观看 | www.自偷自拍.com| 久久人人97超碰香蕉20202| 男人舔女人的私密视频| 狠狠精品人妻久久久久久综合| 男人舔女人的私密视频| 国产在线一区二区三区精| 久久国产精品男人的天堂亚洲| 国产又爽黄色视频| 中文字幕人妻丝袜一区二区| 1024香蕉在线观看| 国产伦理片在线播放av一区| 亚洲中文av在线| 天堂中文最新版在线下载| 亚洲国产欧美一区二区综合| 波多野结衣一区麻豆| 亚洲熟女精品中文字幕| 亚洲伊人久久精品综合| 亚洲国产欧美在线一区| 亚洲综合色网址| 法律面前人人平等表现在哪些方面 | 51午夜福利影视在线观看| 精品国产一区二区三区久久久樱花| 交换朋友夫妻互换小说| 欧美日韩成人在线一区二区| 中文字幕人妻熟女乱码| 亚洲免费av在线视频| 日本猛色少妇xxxxx猛交久久| 黑人猛操日本美女一级片| 我的亚洲天堂| 亚洲中文日韩欧美视频| 亚洲午夜精品一区,二区,三区| 一级a爱视频在线免费观看| 丝袜在线中文字幕| 女人高潮潮喷娇喘18禁视频| 精品视频人人做人人爽| 99精品欧美一区二区三区四区| 成人亚洲精品一区在线观看| 国产精品一区二区在线观看99| 十八禁网站网址无遮挡| 91成年电影在线观看| cao死你这个sao货| 欧美日韩视频精品一区| 两人在一起打扑克的视频| 亚洲国产中文字幕在线视频| 人人澡人人妻人| 咕卡用的链子| 欧美国产精品va在线观看不卡| 午夜福利视频精品| 亚洲午夜精品一区,二区,三区| 性色av一级| 人人妻,人人澡人人爽秒播| 国产xxxxx性猛交| 国产又爽黄色视频| 另类亚洲欧美激情| 国产一区二区三区综合在线观看| 肉色欧美久久久久久久蜜桃| 一区二区三区精品91| 国产国语露脸激情在线看| 欧美日韩亚洲国产一区二区在线观看 | 男女国产视频网站| 母亲3免费完整高清在线观看| 大香蕉久久成人网| 久久香蕉激情| 热99re8久久精品国产| 中文精品一卡2卡3卡4更新| 午夜免费鲁丝| 国产亚洲精品一区二区www | 悠悠久久av| av在线app专区| 久久狼人影院| 80岁老熟妇乱子伦牲交| 宅男免费午夜| 女人久久www免费人成看片| 丝袜美腿诱惑在线| 亚洲精华国产精华精| 久久国产精品大桥未久av| avwww免费| 啦啦啦 在线观看视频| 狠狠狠狠99中文字幕| 亚洲av国产av综合av卡| 日本一区二区免费在线视频| www.精华液| 脱女人内裤的视频| 高清欧美精品videossex| 欧美成人午夜精品| 亚洲精品国产av蜜桃| 中文字幕另类日韩欧美亚洲嫩草| 9191精品国产免费久久| av免费在线观看网站| 涩涩av久久男人的天堂| 国产av一区二区精品久久| 久久久国产精品麻豆| 国产黄色免费在线视频| 国产成人免费观看mmmm| 午夜成年电影在线免费观看| 精品少妇内射三级| 亚洲五月色婷婷综合| 18禁观看日本| 国产淫语在线视频| 亚洲av日韩精品久久久久久密| 纵有疾风起免费观看全集完整版| 免费高清在线观看视频在线观看| 大码成人一级视频| 亚洲av电影在线观看一区二区三区| 黑人操中国人逼视频| 97人妻天天添夜夜摸| 亚洲中文字幕日韩| 91字幕亚洲| 少妇精品久久久久久久| 美女高潮喷水抽搐中文字幕| 精品福利永久在线观看| 精品熟女少妇八av免费久了| 一级片'在线观看视频| 亚洲av成人一区二区三| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品久久久人人做人人爽| 亚洲国产日韩一区二区| 永久免费av网站大全| 亚洲天堂av无毛| 国产男女内射视频| 久久国产亚洲av麻豆专区| 亚洲精品国产av成人精品| 一区二区三区激情视频| 女人高潮潮喷娇喘18禁视频| 99re6热这里在线精品视频| 电影成人av| 国产精品秋霞免费鲁丝片| 久久影院123| 超碰成人久久| 久久久久久久大尺度免费视频| 考比视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲熟女精品中文字幕| 午夜久久久在线观看| 夜夜骑夜夜射夜夜干| 久久久精品区二区三区| 精品一品国产午夜福利视频| 老司机影院毛片| 国产精品偷伦视频观看了| 国产亚洲欧美精品永久| 国产成人精品久久二区二区免费| 99热全是精品| 欧美精品高潮呻吟av久久| 国产亚洲欧美在线一区二区| 大香蕉久久成人网| 国产男女内射视频| 午夜91福利影院| 性色av一级| 国产主播在线观看一区二区| 最近中文字幕2019免费版| 欧美精品一区二区免费开放| 美女扒开内裤让男人捅视频| 亚洲成人手机| 少妇粗大呻吟视频| 精品人妻在线不人妻| 大型av网站在线播放| 国产国语露脸激情在线看| 制服人妻中文乱码| 亚洲一区中文字幕在线| 极品少妇高潮喷水抽搐| 国产欧美日韩精品亚洲av| 久久久久久久久久久久大奶| 99久久精品国产亚洲精品| 国产日韩一区二区三区精品不卡| a 毛片基地| 97在线人人人人妻| 国产在线视频一区二区| av网站免费在线观看视频| 精品久久蜜臀av无| 又黄又粗又硬又大视频| 日韩欧美国产一区二区入口| 午夜福利,免费看| 最近最新免费中文字幕在线| 亚洲免费av在线视频| 国产亚洲精品一区二区www | 久久精品国产亚洲av高清一级| av福利片在线| 9色porny在线观看| 老司机亚洲免费影院| 亚洲成人免费av在线播放| 成人亚洲精品一区在线观看| 久9热在线精品视频| 亚洲精品国产色婷婷电影| 黑人操中国人逼视频| 国产精品一二三区在线看| 80岁老熟妇乱子伦牲交| 日韩三级视频一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 免费高清在线观看视频在线观看| 免费在线观看黄色视频的| 久久青草综合色| 激情视频va一区二区三区| 他把我摸到了高潮在线观看 | 亚洲欧美一区二区三区黑人| 色婷婷av一区二区三区视频| 精品国内亚洲2022精品成人 | 18禁国产床啪视频网站| 精品福利永久在线观看| 国产男人的电影天堂91| 五月天丁香电影| 天天影视国产精品| 亚洲七黄色美女视频| 欧美另类亚洲清纯唯美| 国产成+人综合+亚洲专区| 免费少妇av软件| cao死你这个sao货| 女人精品久久久久毛片| 亚洲欧美一区二区三区久久| 国产1区2区3区精品| 亚洲伊人色综图| 美女高潮喷水抽搐中文字幕| 超色免费av| 最新在线观看一区二区三区| 国产成+人综合+亚洲专区| 精品卡一卡二卡四卡免费| 18禁观看日本| 老汉色av国产亚洲站长工具| av福利片在线| 男男h啪啪无遮挡| 韩国精品一区二区三区| 捣出白浆h1v1| 女警被强在线播放| 脱女人内裤的视频| 男人舔女人的私密视频| 女性生殖器流出的白浆| 亚洲一卡2卡3卡4卡5卡精品中文| 啦啦啦视频在线资源免费观看| 日韩中文字幕欧美一区二区| 国产91精品成人一区二区三区 | 两性夫妻黄色片| 国产一区二区 视频在线| 天堂俺去俺来也www色官网| 热99久久久久精品小说推荐| 亚洲激情五月婷婷啪啪| 欧美大码av| 亚洲成人免费电影在线观看| 中文字幕人妻丝袜一区二区| 精品国产乱码久久久久久小说| 亚洲av成人一区二区三| 亚洲伊人色综图| av有码第一页| 亚洲一卡2卡3卡4卡5卡精品中文| xxxhd国产人妻xxx| 超碰97精品在线观看| 高清视频免费观看一区二区| 18禁黄网站禁片午夜丰满| 青青草视频在线视频观看| 不卡一级毛片| 亚洲中文av在线| 女人精品久久久久毛片| 五月天丁香电影| 久久久精品免费免费高清| 国产精品成人在线| 亚洲精品av麻豆狂野| 日韩制服丝袜自拍偷拍| 99精国产麻豆久久婷婷| 夜夜夜夜夜久久久久| 日韩免费高清中文字幕av| 啦啦啦免费观看视频1| 男女下面插进去视频免费观看| 高清av免费在线| 国产精品免费视频内射| 欧美午夜高清在线| 在线观看免费日韩欧美大片| 日日摸夜夜添夜夜添小说| 亚洲av成人不卡在线观看播放网 | 久久国产精品人妻蜜桃| 天堂俺去俺来也www色官网| 黄片小视频在线播放| 亚洲美女黄色视频免费看| 亚洲国产欧美网| 在线观看人妻少妇| 久久精品国产亚洲av高清一级| 国产激情久久老熟女| 久久久精品区二区三区| 日韩三级视频一区二区三区| 99久久精品国产亚洲精品| 黄色视频在线播放观看不卡| 1024香蕉在线观看| 免费人妻精品一区二区三区视频| 亚洲av电影在线观看一区二区三区| 午夜两性在线视频| 亚洲久久久国产精品| 极品人妻少妇av视频| 丝袜脚勾引网站| 一区福利在线观看| 桃红色精品国产亚洲av| 久久人人爽av亚洲精品天堂| 18禁观看日本| 亚洲七黄色美女视频| 国产高清国产精品国产三级| 老司机靠b影院| 欧美日韩亚洲综合一区二区三区_| av网站在线播放免费| 精品少妇内射三级| 国产亚洲av片在线观看秒播厂| 日本a在线网址| 精品第一国产精品| 免费在线观看影片大全网站| 俄罗斯特黄特色一大片| 久久久精品94久久精品| 亚洲精品中文字幕在线视频| 久久国产精品大桥未久av| 日本猛色少妇xxxxx猛交久久| 亚洲熟女毛片儿| 一级毛片女人18水好多| 97人妻天天添夜夜摸| 精品乱码久久久久久99久播| 国产成人一区二区三区免费视频网站| 亚洲欧美清纯卡通| 国产av一区二区精品久久| 后天国语完整版免费观看| 免费一级毛片在线播放高清视频 | 十分钟在线观看高清视频www| 中文字幕最新亚洲高清| 国产一级毛片在线| av有码第一页| 欧美黄色淫秽网站| 九色亚洲精品在线播放| 亚洲男人天堂网一区| 最黄视频免费看| 青青草视频在线视频观看| 好男人电影高清在线观看| 美女脱内裤让男人舔精品视频| 日韩中文字幕视频在线看片| av免费在线观看网站| 捣出白浆h1v1| 色综合欧美亚洲国产小说| 一区二区三区激情视频| 国产精品av久久久久免费| av一本久久久久| 国产成人欧美| 国产福利在线免费观看视频| 黑人欧美特级aaaaaa片| 精品亚洲成a人片在线观看| 中国美女看黄片| 午夜福利视频精品| 午夜免费观看性视频| 欧美精品一区二区免费开放| 人人妻人人添人人爽欧美一区卜| 乱人伦中国视频| 亚洲av日韩在线播放| 9色porny在线观看| 少妇精品久久久久久久| 最近最新中文字幕大全免费视频| 我的亚洲天堂| 国产野战对白在线观看| 免费观看人在逋| 男女国产视频网站| 男人舔女人的私密视频| 亚洲黑人精品在线| 在线观看免费高清a一片| 十八禁高潮呻吟视频| 亚洲精品日韩在线中文字幕| 中文字幕av电影在线播放| 在线观看一区二区三区激情| 十八禁人妻一区二区| 亚洲av片天天在线观看| 考比视频在线观看| 国产欧美日韩综合在线一区二区| 咕卡用的链子| 国产91精品成人一区二区三区 | 精品国产一区二区久久| 亚洲精品一卡2卡三卡4卡5卡 | 女警被强在线播放| 曰老女人黄片| 精品人妻熟女毛片av久久网站| 中文字幕色久视频| 悠悠久久av| 电影成人av| 久久 成人 亚洲| 在线永久观看黄色视频| 一区在线观看完整版| 中文字幕精品免费在线观看视频| 人人妻人人爽人人添夜夜欢视频| 日韩一卡2卡3卡4卡2021年| 亚洲欧美精品综合一区二区三区| 视频区欧美日本亚洲| 搡老乐熟女国产| 欧美精品高潮呻吟av久久| 大香蕉久久成人网| 无限看片的www在线观看| 99国产精品免费福利视频| 女警被强在线播放| 黄网站色视频无遮挡免费观看| 免费观看av网站的网址| 国精品久久久久久国模美| 久热这里只有精品99| 日韩大片免费观看网站| 又紧又爽又黄一区二区| 肉色欧美久久久久久久蜜桃| 人人妻人人添人人爽欧美一区卜| 国产精品影院久久| 国产色视频综合| 国产精品熟女久久久久浪| 国产亚洲欧美在线一区二区| 亚洲国产精品一区二区三区在线| 黑人猛操日本美女一级片| 日韩欧美国产一区二区入口| 91国产中文字幕| 一边摸一边抽搐一进一出视频| 国产欧美亚洲国产| √禁漫天堂资源中文www| 黑人巨大精品欧美一区二区蜜桃| 午夜福利在线观看吧| 丰满少妇做爰视频| 久久久久精品人妻al黑| 午夜成年电影在线免费观看| 丝袜喷水一区| 精品国产一区二区久久| 建设人人有责人人尽责人人享有的| 最近最新中文字幕大全免费视频| 午夜日韩欧美国产| 99精品欧美一区二区三区四区| 人妻久久中文字幕网| 99国产极品粉嫩在线观看| 成人手机av| 精品福利观看| 免费女性裸体啪啪无遮挡网站| av在线播放精品| 亚洲全国av大片| 国产精品国产av在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 黑人巨大精品欧美一区二区蜜桃| 大香蕉久久成人网| 美女国产高潮福利片在线看| 天天添夜夜摸| 精品国内亚洲2022精品成人 | 人妻久久中文字幕网| 亚洲一区中文字幕在线| 亚洲久久久国产精品| 色综合欧美亚洲国产小说| 女人高潮潮喷娇喘18禁视频| 国产色视频综合| 亚洲五月色婷婷综合| 天堂中文最新版在线下载| 久久国产亚洲av麻豆专区| 亚洲第一青青草原| 夫妻午夜视频| 搡老岳熟女国产| 国产成人av教育| 国产精品一二三区在线看| 麻豆av在线久日| 国产精品99久久99久久久不卡| 国内毛片毛片毛片毛片毛片| 在线观看免费高清a一片| 亚洲视频免费观看视频| 亚洲国产中文字幕在线视频| 亚洲欧美色中文字幕在线| 视频区图区小说| 老司机靠b影院| 国产色视频综合| 亚洲欧美精品综合一区二区三区| 精品免费久久久久久久清纯 | 久久九九热精品免费| 日韩 亚洲 欧美在线| 丝袜美腿诱惑在线| 国产成人精品久久二区二区免费| 亚洲精品美女久久av网站| 制服人妻中文乱码| av视频免费观看在线观看| 高清av免费在线| 男女无遮挡免费网站观看| 欧美一级毛片孕妇| 精品亚洲乱码少妇综合久久| 国产国语露脸激情在线看| 亚洲精品中文字幕在线视频| 老鸭窝网址在线观看| 成年人午夜在线观看视频| 99国产综合亚洲精品| 亚洲国产欧美网| 亚洲中文av在线| av国产精品久久久久影院| 精品人妻熟女毛片av久久网站| 国产精品一区二区在线不卡| 欧美黑人精品巨大| 午夜激情久久久久久久| 美女午夜性视频免费| 日韩大片免费观看网站| av天堂在线播放| 两人在一起打扑克的视频| 欧美+亚洲+日韩+国产| 精品国产一区二区久久| 自线自在国产av| 中文字幕人妻熟女乱码| 欧美亚洲日本最大视频资源| 亚洲色图 男人天堂 中文字幕| 日韩中文字幕欧美一区二区| 中亚洲国语对白在线视频| a 毛片基地| 亚洲av欧美aⅴ国产| 亚洲国产精品成人久久小说| 欧美在线一区亚洲| 丝袜脚勾引网站| 51午夜福利影视在线观看| 欧美 亚洲 国产 日韩一| 亚洲欧美精品自产自拍| 亚洲自偷自拍图片 自拍| 免费观看av网站的网址| 久久精品亚洲熟妇少妇任你| 香蕉丝袜av| 国产精品熟女久久久久浪| 亚洲欧美精品自产自拍| 免费日韩欧美在线观看| 久久精品国产亚洲av高清一级| 亚洲一区中文字幕在线| 欧美精品一区二区大全| 成年人午夜在线观看视频| 美女大奶头黄色视频| 国产亚洲av高清不卡| 超色免费av| 亚洲精品国产区一区二| 99香蕉大伊视频| 老熟妇仑乱视频hdxx| 国产成人av教育| 99久久人妻综合| 免费不卡黄色视频| 新久久久久国产一级毛片| 可以免费在线观看a视频的电影网站| 国产真人三级小视频在线观看| 成人国产av品久久久| 国产精品1区2区在线观看. | 一区二区三区四区激情视频| av有码第一页| 美国免费a级毛片| 久久人妻福利社区极品人妻图片| 一本综合久久免费| 满18在线观看网站| 一二三四在线观看免费中文在| 中国美女看黄片| 美女视频免费永久观看网站| 热99re8久久精品国产| 极品少妇高潮喷水抽搐| 少妇的丰满在线观看| 久久精品熟女亚洲av麻豆精品| 久久久精品94久久精品| 亚洲精品一二三| 国产成人av激情在线播放| 国产欧美日韩精品亚洲av| 精品高清国产在线一区| 国产精品一区二区在线观看99| 99九九在线精品视频| 日本a在线网址| 国产精品免费视频内射| 国产成人影院久久av| 欧美 日韩 精品 国产| 国产成人一区二区三区免费视频网站| 女人高潮潮喷娇喘18禁视频| 亚洲国产欧美在线一区| h视频一区二区三区| 久久ye,这里只有精品| 在线精品无人区一区二区三| 在线观看免费高清a一片| 性高湖久久久久久久久免费观看| 精品少妇内射三级| 国产免费福利视频在线观看| 黄色视频在线播放观看不卡| 成人手机av| 日韩欧美国产一区二区入口| 亚洲成国产人片在线观看| 久久久精品免费免费高清| 久久国产精品人妻蜜桃| 日本黄色日本黄色录像| 免费一级毛片在线播放高清视频 | 韩国精品一区二区三区| 国产成人av教育| 日韩视频在线欧美| 成人18禁高潮啪啪吃奶动态图| 亚洲精品美女久久av网站| 热re99久久国产66热| 中文字幕另类日韩欧美亚洲嫩草| 一边摸一边抽搐一进一出视频| 亚洲人成电影观看| 久久久久久久大尺度免费视频| 首页视频小说图片口味搜索|