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

    電四極躍遷對電子束離子阱等離子體中離子能級布居的影響*

    2022-10-16 09:23:22孟舉何貞岑顏君吳澤清姚科李冀光吳勇3王建國
    物理學(xué)報 2022年19期
    關(guān)鍵詞:電子密度譜線能級

    孟舉 何貞岑 顏君 吳澤清 姚科? 李冀光? 吳勇3) 王建國

    1) (北京應(yīng)用物理與計算數(shù)學(xué)研究所 計算物理重點(diǎn)實(shí)驗室,北京 100088)

    2) (復(fù)旦大學(xué)核物理與離子束應(yīng)用教育部重點(diǎn)實(shí)驗室,上海 200433)

    3) (北京大學(xué)應(yīng)用物理與技術(shù)研究中心,北京 100871)

    在穩(wěn)態(tài)近似下,通過構(gòu)建三能級體系碰撞輻射模型,解析研究了電四極(E2)躍遷對等離子體中離子能級布居的影響.研究發(fā)現(xiàn),隨著原子序數(shù)增加,E2 躍遷速率逐漸增強(qiáng),E2 躍遷在低電子密度條件下對離子能級布居的影響愈發(fā)顯著.進(jìn)一步地,以電子束離子阱中類鐵鉬(Z=42)和鈾(Z=92)等離子體為例,數(shù)值求解了包含不同退激通道的離子能級布居.在此基礎(chǔ)上,分析了考慮E2 退激通道導(dǎo)致的能級布居變化對基組態(tài)磁偶極(M1)躍遷線強(qiáng)比的影響,并指出在利用高離化態(tài)離子M1 躍遷線強(qiáng)比進(jìn)行等離子體電子密度診斷時E2 退激發(fā)通道的重要性.

    1 引言

    等離子體電子密度是表征等離子體狀態(tài)的一個基本物理量,對該物理量的準(zhǔn)確確定比較困難[1,2].另一方面,等離子體電子密度直接影響等離子體內(nèi)多種原子動力學(xué)過程反應(yīng)速率,各種原子動力學(xué)過程相互競爭會導(dǎo)致特定躍遷譜線強(qiáng)度對等離子體電子密度十分敏感,因此利用等離子體光譜對其診斷是常用的方法之一.黃文忠等[3]根據(jù)類銅銀離子的4d2D5/2-4f2F7/2和4p2P3/2-4d2D5/2兩條譜線的線強(qiáng)比診斷得到了激光加熱銀靶產(chǎn)生等離子體的電子密度.Feldman 等[4]研究了類鈦離子(Xe32+,Gd42+,W52+,Bi61+,U70+)基組態(tài)中處于可見光和近紫外(UV)波段的磁偶極(M1)躍遷線隨等離子體電子密度的變化規(guī)律,并指出來自于J=2 →J=3的M1 躍遷可以作為高溫托卡馬克裝置的密度診斷工具.作為該工作的延續(xù),Doron 等[5]于同年報道了鎢、鉛、鈾等若干個高離化度離子的N 殼層(4pk和4dk)基組態(tài)中處于UV 波段的M1 躍遷線,預(yù)言了若干條能夠在電子束離子阱(EBIT)和托卡馬克裝置中觀測到且較亮的躍遷線,以及部分對電子密度敏感的譜線強(qiáng)度比.Ralchenko 等[6]對類鎳鎢(W46+)的基組態(tài)3d10和第一激發(fā)組態(tài)3d94s 之間的電四極(E2)與磁八極(M3)躍遷之間的線強(qiáng)比進(jìn)行了細(xì)致研究,發(fā)現(xiàn)該線強(qiáng)比呈現(xiàn)出較強(qiáng)的電子密度依賴性,可以用于托卡馬克等離子體的密度診斷.之后,他們又對W48+-W53+的基組態(tài)能級間的M1 躍遷線強(qiáng)比進(jìn)行了研究,報道了若干組對電子密度敏感的線強(qiáng)比[7].近來,丁曉彬等[8]報道了類鎘鎢基組態(tài)間的可用于等離子體電子密度診斷的躍遷線及線強(qiáng)比.Silwal等[1]對高離化態(tài)釔(Y)離子的光譜進(jìn)行了細(xì)致的實(shí)驗和理論研究,并分析了Y26+-Y33+離子之間的M1 躍遷與電偶極躍遷(E1)線強(qiáng)比對等離子體電子密度的依賴關(guān)系,報道了若干組在聚變等離子體密度診斷方面具有潛在應(yīng)用價值的線強(qiáng)比.最近,復(fù)旦大學(xué)現(xiàn)代物理研究所利用新建的低能EBIT裝置(CUBIT)對類鐵和類錳鉬離子(Mo16+,Mo17+)的基組態(tài)M1 躍遷進(jìn)行了觀測,首次測量得到10 條M1 躍遷線的波長和強(qiáng)度,并預(yù)言了若干組對電子密度敏感的線強(qiáng)比,建議這些譜線可以作為潛在的等離子體電子密度診斷工具[9].

    EBIT 裝置被廣泛用于實(shí)驗室產(chǎn)生等離子體,是獲取高精度原子參數(shù)、研究等離子體中各種原子過程的有力實(shí)驗平臺.在EBIT 裝置中,對等離子體的約束時間(電子作用時間)維持在秒量級,而在電子作用下,一般過程的特征時間在毫秒量級或者更短,所以在約束時間內(nèi)EBIT 等離子體很容易達(dá)到穩(wěn)態(tài).因此,為了研究EBIT 等離子體中譜線強(qiáng)度與等離子體電子密度的依賴關(guān)系,理論上需要構(gòu)建碰撞輻射模型(CRM)并求解穩(wěn)態(tài)速率方程得到離子能級布居.基于EBIT 等離子體的特點(diǎn),通常只需要考慮同價離子不同能級間的輻射退激發(fā)及電子碰撞激發(fā)和退激發(fā)過程[5,10-13].對于輻射退激發(fā),理想的情況是包含所有類型的輻射躍遷,但是這會大大增加計算規(guī)模,因此需要根據(jù)具體的研究對象選取主要的輻射退激通道.一般地,由于E2 或者更高階的禁戒躍遷速率較小,它對退激過程的貢獻(xiàn)可以忽略,因此輻射退激過程只需考慮E1 和M1 躍遷.但是E2 躍遷與原子序數(shù)Z和有效電荷Q的標(biāo)度關(guān)系可以達(dá)到AE2∝Z5Q11[14],而M1 和E1 躍遷率對于原子序數(shù)的標(biāo)度關(guān)系分別為AM1∝Z3Q9和AE1∝Z3Q7,對于Z較高的離子,E2躍遷速率可以與M1、甚至E1的躍遷速率相比,這些通道可能引起能級布居變化,此時在CRM 中是否考慮E2 躍遷是一個值得細(xì)致研究的問題.

    本文將圍繞E2 躍遷通道對等離子體中離子能級布居和譜線強(qiáng)度的重要性開展研究.首先通過三能級原子模型解析分析了E2 躍遷通道在不同電子密度下對能級布居的影響;進(jìn)一步以EBIT 中原子序數(shù)較小的Mo(Z=42)和原子序數(shù)較大的U(Z=92)等離子體為例,研究不同電子密度條件下E2躍遷通道對類鐵Mo16+和U66+離子布居、譜線強(qiáng)度和線強(qiáng)比的影響,并討論其對基組態(tài)M1 躍遷強(qiáng)度比進(jìn)行等離子體密度診斷產(chǎn)生的影響.

    2 理論方法和三能級原子模型

    等離子體中處于特定離化態(tài)的離子能級布居可以用碰撞輻射速率方程表示為

    其中Ni表示第i能級的離子布居,Aji表示第j能級到第i能級的輻射躍遷速率,ne表示電子數(shù)密度,cji表示第j能級到第i能級的碰撞激發(fā)速率系數(shù),dji表示對應(yīng)的碰撞退激發(fā)速率系數(shù).碰撞速率系數(shù)可以通過對電子速度分布函數(shù)和碰撞截面的卷積得到.由于EBIT 中電子能量的準(zhǔn)單能性,電子速度分布采用高斯分布,電子的能量在相應(yīng)電離閾附近取值.根據(jù)穩(wěn)態(tài)近似和歸一化條件,通過求解速率方程可以得到EBIT 等離子體的離子能級布居.

    為了能夠得到較為普遍的規(guī)律,首先以三能級原子系統(tǒng)為研究對象,構(gòu)建并解析求解速率方程.在速率方程中,將離子態(tài)用φ 表示,并假定按激發(fā)能升序分別記為φ1,φ2,φ3;同時φ3到φ1和φ2的退激發(fā)為M1 躍遷,φ2到φ1退激發(fā)為E2 躍遷,相應(yīng)的速率系數(shù)標(biāo)記與方程(1)一致;與每個態(tài)相關(guān)的速率和布居分別用下標(biāo)1,2,3 表示,比如3 個態(tài)的布居分別為N1,N2,N3.構(gòu)建兩個碰撞輻射模型,模型I 不包含E2 躍遷,此時的輻射退激通道只包括φ3→φ2,φ3→φ1,即在不考慮E2 躍遷的情況下φ2為亞穩(wěn)態(tài),只能通過碰撞過程退激發(fā);模型II 包含E2 躍遷,此時需要考慮φ2→φ1,φ3→φ2,φ3→φ1輻射退激,即3 個態(tài)之間的躍遷均可發(fā)生,如圖1 所示.除此之外,兩個模型包含的碰撞激發(fā)和退激發(fā)通道是相同的,即3 個態(tài)之間均可通過碰撞過程相互到達(dá).分別求解模型I 和模型II的速率方程,可以得到兩個模型的解析解,分別如下:

    圖1 三能級原子體系能級圖 (a) 模型I 不包含E2 躍遷;(b) 模型II 包含E2 躍遷示意圖Fig.1.Energy levels of the three-level atomic system: (a) Model I without the E2 transition;(b) Model II with inclusion of the E2 transition.

    式中,U和U'分別是模型I 和模型II 中3 個態(tài)布居表達(dá)式中的所有分子項之和,其表達(dá)式分別為

    從(2)—(4)式可以看出,離子態(tài)布居表達(dá)式中的分子和分母項均可以寫成關(guān)于電子數(shù)密度的二次項表達(dá)式,其系數(shù)是關(guān)于躍遷速率或者碰撞速率系數(shù)之間乘積的和,并且能量越高的激發(fā)態(tài),其分子項的個數(shù)越少.

    下面我們分析解析解在低電子密度和高電子密度極限情況下的表達(dá)形式.當(dāng)電子數(shù)密度趨于0 時,模型I 和模型II的布居結(jié)果為

    通過(7)—(9)式可以看出,在模型I 中,φ2能級會分布有一定的布居數(shù),該布居數(shù)值與電子數(shù)密度無關(guān),只與各個躍遷通道和碰撞的速率系數(shù)有關(guān);而對于模型II,所有能級布居都會集中在基態(tài),激發(fā)態(tài)的布居為0.由此表明,E2 躍遷在低電子密度條件下對離子布居會產(chǎn)生較大影響.在等離子體中,離子的能級布居主要是輻射和碰撞過程的競爭所決定.在低電子密度條件下,當(dāng)?shù)入x子體中電子碰撞速率小于輻射躍遷速率時,輻射退激發(fā)過程占主導(dǎo)作用,此時新的輻射躍遷通道對于離子布居就會產(chǎn)生較大的影響,特別是對于只能通過新打開的躍遷通道退激發(fā)的能級.

    當(dāng)電子密度趨向于無窮大時,模型I 和模型II的結(jié)果可以寫成下面的形式:

    從(10)—(12)式可以看出,模型I 和模型II 計算得到的離子布居數(shù)相等,即當(dāng)電子密度很大時,是否包括E2 躍遷通道對最終的離子布居沒有影響,只與碰撞速率系數(shù)相關(guān).這是因為在高電子密度條件下,等離子體內(nèi)部碰撞速率遠(yuǎn)大于輻射躍遷速率,此時等離子體的離子能級布居主要由碰撞過程決定,輻射過程對于布居的影響可以忽略,因此打開新的輻射躍遷通道對布居幾乎沒有影響.

    3 應(yīng)用與討論

    3.1 E2 躍遷對類鐵鉬和鈾等離子體能級布居的影響

    第2 節(jié)中以三能級原子體系為對象,針對電子密度極低和極高的情況,解析研究了E2 退激發(fā)通道對等離子體中離子能級布居的影響.實(shí)際上,原子能級遠(yuǎn)多于3 個,等離子體中離子能級布居需要對速率方程進(jìn)行數(shù)值求解得到.另外,等離子體中的電子密度往往不處于極限情況,各種原子動力學(xué)過程之間的競爭也隨密度改變,對能級布居的影響也會不同.為了分析E2 退激通道對實(shí)際等離子體中離子能級布居的影響,以EBIT 中類鐵鉬和鈾等離子體為例進(jìn)行研究.圖2 分別給出了類鐵鉬離子和鈾離子3d8基組態(tài)能級和躍遷示意圖,圖中實(shí)線表示M1 躍遷,虛線表示E2 躍遷.需要注意的是,由于高Z體系中相對強(qiáng)的相對論效應(yīng)導(dǎo)致類鐵鈾離子基組態(tài)中較大的精細(xì)結(jié)構(gòu)分裂,甚至部分能級順序不同于等電子系列中的低Z體系,因此對于兩個體系序號相同的能級并不一一對應(yīng).從圖2可以看出,根據(jù)躍遷選擇定則類鐵鉬離子的4 號能級類鐵鈾離子的1 號能級和2 號能級均不能通過M1 躍遷通道到達(dá)基態(tài),只有當(dāng)考慮E2 躍遷之后,這些激發(fā)態(tài)才能通過輻射退激發(fā)到達(dá)基態(tài).

    圖2 類鐵鉬(Mo16+)離子(a)和鈾(U66+)離子(b)基組態(tài)內(nèi)的禁戒躍遷,其中實(shí)線為M1 躍遷,虛線為E2 躍遷,圖中的數(shù)字(0,1,2,···)表示能級序號Fig.2.Energy-level diagram of the ground state configuration of Fe-like Mo16+ (a) and Fe-like U66+ (b) ions.Solid lines represent the M1 transitions and dashed lines represent the E2 transitions.The numbers (0,1,2,···) correspond to the energy levels labels.

    利用FAC 程序[15]的CRM 模型[16-19],數(shù)值求解了兩種等離子體的離子能級布居.精確的能級布居依賴于CRM 模型中包括的原子態(tài)是否完備及相關(guān)原子參數(shù)的精度.在目前的CRM 模型中,包含的能級涉及3d8,3d74l(l≤3),3d75l(l≤4),3p53d9,3p53d84l(l≤3),3s13d9,3d64s2,3p43d10組態(tài).相關(guān)的原子參數(shù)采用相對論多體微擾方法(RMBPT)計算得到.該方法的主要思想是將組態(tài)空間分成M 和N 兩個子空間,M 空間包含與當(dāng)前研究相關(guān)的重要組態(tài),N 空間是由M 空間的組態(tài)經(jīng)過單雙激發(fā)產(chǎn)生的組態(tài)構(gòu)成.通過引入投影算符,可以將體系的哈密頓量本征值問題轉(zhuǎn)化成有效哈密頓量的本征值問題,最后對角化有效哈密頓量可以得到本征態(tài)能量及波函數(shù)[20,21].在當(dāng)前計算中,將CRM中的組態(tài)全部包含在M 空間中,對于N 空間,將單電子激發(fā)和雙電子激發(fā)中的第一個電子可以激發(fā)的最大主量子數(shù)(n1)設(shè)置為170,將雙激發(fā)中的第二個電子能夠激發(fā)的最大主量子數(shù)(n2)設(shè)置為65.表1 列出了計算得到的Mo16+和U66+基組態(tài)能級的激發(fā)能和Mo16+離子的NIST 結(jié)果[22,23].從表1可以看出,RMBPT 結(jié)果和NIST的結(jié)果符合得很好,兩組數(shù)據(jù)的最大誤差和平均誤差分別為0.35%和0.09%,由此可以得出當(dāng)前的理論計算結(jié)果是可靠的.對于CRM的完備性驗證,構(gòu)建了新的CR模型,新模型在原有模型的基礎(chǔ)上增加了新的組態(tài),比如3d76l(l≤5),3p53d85l(l≤4),3s13d84l(l≤3).為了保證計算可行,新增加的組態(tài)采用相對論組態(tài)相互作用方法并且忽略與其他組態(tài)間的相互作用得到.計算結(jié)果表明新增加的組態(tài)對離子布居影響很小,可以認(rèn)為CRM 中包含的原子態(tài)是完備的.另外,實(shí)際計算中電子速度分布采用高斯型分布,對于鉬離子和鈾離子,其電子的能量分別為640 eV和8145 eV,半高寬基于EBIT 實(shí)驗參數(shù)采用電子能量的5%.與分析三能級體系類似,構(gòu)建兩個模型分別記為模型I 和模型II,模型I 只考慮E1 和M1 輻射退激通道,模型II 在模型I的基礎(chǔ)上增加了E2 躍遷通道.除此之外,兩個模型采用的其他原子參數(shù)均相同.

    表1 類鐵鉬離子和鈾離子的基組態(tài)精細(xì)能級激發(fā)能Table 1.Excitation energies of the lowest excited levels of Fe-like Mo16+,U66+ ions.

    圖3 分別給出了類鐵鉬離子和鈾離子的3d8基組態(tài)精細(xì)能級布居隨電子密度的變化關(guān)系,圖中實(shí)線是模型I的計算結(jié)果,用E1M1 表示,虛線是模型II的計算結(jié)果,用E1M1E2 表示.圖3(a)展示了類鐵鉬離子的計算結(jié)果,從圖3 可以看出,當(dāng)電子密度較大時(大于1015cm—3),模型I 和模型II 計算得到的各能級布居差異很小,這說明E2 躍遷通道在電子密度較高時對鉬離子布居幾乎沒有影響,該結(jié)論也與三能級原子模型一致.而在電子密度小于106cm—3的區(qū)間,由于E2 躍遷的影響,兩個模型計算的部分能級布居存在一定差異,特別是對于4 號能級模型I 和模型II 計算的布居明顯不同.模型I的計算表明該能級仍會存在一定數(shù)量的布居,占比0.2%,但是根據(jù)模型II的計算,該激發(fā)態(tài)能級的布居已經(jīng)趨近于0.這是因為模型I 中沒有包含E2 躍遷通道,4 號能級只能通過碰撞退激過程退激發(fā),在低電子密度環(huán)境下其碰撞速率很小,因此它們具有一定的布居;而模型II 包含了E2 躍遷,使得該能級可以借助于E2 通道到達(dá)基態(tài),故其布居在低電子密度時趨近于0.同時可以看出,對于類鐵鉬離子,由于E2 躍遷速率較小(最大速率僅47 s—1),只在較低的電子密度下E2 躍遷速率才會超過碰撞退激發(fā)速率發(fā)揮作用,因此在很寬的電子密度區(qū)間E2 躍遷對能級布居的影響都很小.

    圖3 類鐵鉬離子(a)和鈾離子(b)處于基組態(tài)3d8 能級的離子布居對密度的依賴關(guān)系,圖中的數(shù)字(0,1,2,···)對應(yīng)于圖1 中類鐵鉬離子和鈾離子的能級序號Fig.3.The electron-density dependence of the population distribution of energy levels belonging to 3d8 configuration of Fe-like Mo16+(a) and U66+(b) ions.The numbers (0,1,2,···) correspond to the energy levels of Mo16+ and U66+ of Fig.1.

    圖3(b)展示了類鐵鈾離子的計算結(jié)果.從圖3(b)可以看出,在電子密度大于1019cm—3時,兩個模型計算得到的基組態(tài)能級布居差別很小.隨著電子密度的降低,除了1 號能級和2 號能級外,模型I 和模型II 計算得到激發(fā)態(tài)能級布居均很快趨近于0.根據(jù)模型I的計算,當(dāng)電子密度小于107cm—3時,激發(fā)態(tài)仍會占據(jù)一定數(shù)量的布居,分別占比22%和4%,但是模型II的計算結(jié)果表明其布居在該電子密度區(qū)間趨近于0.如前所述,此二能級和不能通過M1 輻射退激到基態(tài),所以當(dāng)CRM 中不包含E2 躍遷通道時,它們只能通過碰撞過程退激發(fā),在低電子密度環(huán)境下其碰撞退激發(fā)速率很小,因此在模型I 中它們具有一定的布居,這也與三能級模型得到的結(jié)論一致.值得注意的是,類鐵鈾離子的E2 躍遷率達(dá)到105s—1,比類鐵鉬離子大了幾個量級,因此對能級布居產(chǎn)生影響時的電子密度也較高.

    通過圖3 兩種離子能級布居結(jié)果的比較可以看出,隨著電子密度的降低,E2 躍遷對類鐵鈾離子布居的影響要大于其對類鐵鉬離子的影響,而這主要是由于E2 躍遷速率增大造成的.對于Mo 離子,原子序數(shù)相對較小,最大的E2 躍遷速率僅為最大的M1 躍遷速率的1%,因此在電子密度從高到低的變化過程中,E2 躍遷對能級布居的影響并不顯著;而對于高Z的鈾離子,E2 躍遷速率增大了近4 個量級,對能級布居產(chǎn)生的影響也相對更大.此外,由于鈾離子中各種輻射退激發(fā)速率比鉬離子的大,與碰撞速率達(dá)到可比擬時所需要的電子密度相應(yīng)增大,因此鈾離子的布居對電子密度敏感的位置相較于鉬離子整體向高密度區(qū)間移動.

    3.2 E2 躍遷對類鐵鉬和鈾離子譜線強(qiáng)度的影響

    對等離子體中電子密度敏感的躍遷線是因為等離子體內(nèi)部輻射、碰撞等原子動力學(xué)過程對產(chǎn)生譜線的上能級布居的貢獻(xiàn)相當(dāng),且相互競爭的結(jié)果.這樣如果其中一個物理量發(fā)生變化就會導(dǎo)致布居出現(xiàn)明顯變化,從而使譜線強(qiáng)度產(chǎn)生改變.因此,E2 躍遷對布居的影響會直接影響譜線強(qiáng)度,同時譜線強(qiáng)度對電子密度的依賴關(guān)系也會發(fā)生改變.基于此原因,在本節(jié)中就E2 躍遷對M1 譜線強(qiáng)度的影響作進(jìn)一步詳細(xì)的研究.

    表2 列出了所研究的類鐵鉬和鈾離子的M1躍遷,其中包括躍遷涉及的上下能級、躍遷能和波長,躍遷速率和振子強(qiáng)度.圖4 呈現(xiàn)了這些M1 躍遷的相對譜線強(qiáng)度隨電子密度的變化關(guān)系.實(shí)線是速率方程求解離子布居中只考慮E1 和M1 躍遷的計算結(jié)果,虛線是包含了E2 躍遷的結(jié)果.對于圖中的譜線線強(qiáng)進(jìn)行歸一化處理.從圖4 中可以看出,E2 躍遷在高電子密度情況下對鉬離子和鈾離子的譜線強(qiáng)度影響很小;但隨著密度降低這種影響在類鐵鈾離子中愈發(fā)顯著,某些譜線強(qiáng)度隨密度的改變有數(shù)倍的變化.由于譜線強(qiáng)度與對應(yīng)躍遷上能級布居和躍遷速率有關(guān),所以E2 躍遷對譜線強(qiáng)度的影響來源于其對布居的影響.也就是說,對于等離子中原子序數(shù)較高的離子,E2 等高階禁戒躍遷對布居和譜線強(qiáng)度的影響是不能被忽略的.

    圖4 類鐵鉬離子(a)和鈾離子(b)基組態(tài)內(nèi)M1 躍遷譜線強(qiáng)度對密度的依賴關(guān)系,圖中的(0 —1,···)是相應(yīng)躍遷,數(shù)字表示該躍遷的下能級和上能級序號Fig.4.The electron-density dependence of the intensity ratios for the M1 transitions from the ground configuration of Mo16+(a) and U66+(b).The numbers (0 —1,···) correspond to the lower and upper energy levels of the lines,respectively.

    表2 類鐵鉬離子和鈾離子的基組態(tài)M1 躍遷的躍遷能ΔE、波長λ、躍遷速率A 和振子強(qiáng)度gfTable 2.Transition energies ΔE,wavelengths λ,transition rates A and oscillator strength gf for the M1 transitions in the ground configuration of Fe-like Mo16+ and U66+.

    3.3 E2 躍遷對類鐵鉬和鈾離子線強(qiáng)比的影響

    如前所述,高離化度離子的M1 躍遷譜線在天體和實(shí)驗室等離子體診斷方面有著重要作用[6,7,24-26].在之前的工作中,針對低Z的16 價和17 價鉬離子基組態(tài)中M1 躍遷譜線開展了實(shí)驗和理論研究,并發(fā)現(xiàn)了若干組對等離子體電子密度敏感的M1 線強(qiáng)比[9].本工作進(jìn)一步研究了E2 躍遷對M1 線強(qiáng)比與等離子體電子密度依賴關(guān)系的影響.圖5 展示了E2 躍遷對類鐵鉬離子和類鐵鈾離子線強(qiáng)比的密度依賴性的影響,圖中的實(shí)線和虛線分別是不包含和包含E2 躍遷的計算結(jié)果.圖中選擇了較為敏感的線強(qiáng)比進(jìn)行說明,對于鉬離子選擇的線強(qiáng)比分別為9/1,9/4,9/7,2/1,2/4,對于鈾離子選擇的線強(qiáng)比為7/2,7/4,7/8,3/2,3/4,這里的羅馬數(shù)字對應(yīng)于表1 中的躍遷序號.從圖5 可以看出,類鐵鉬離子的線強(qiáng)比在109—1013cm—3電子密度區(qū)間最為敏感,該密度區(qū)間與EBIT 裝置的典型密度區(qū)間重合,可用于EBIT 等離子體的密度診斷,鈾離子的密度敏感區(qū)間位于1016—1019cm—3范圍,該敏感區(qū)間趨向于更高的電子密度,其原因是鈾離子中更大的輻射退激發(fā)速率,碰撞過程和輻射過程達(dá)到可比擬時需要的電子密度更大.此外,從圖5 中可以看出,在高電子密度區(qū)間,E2 躍遷對鉬離子和鈾離子的線強(qiáng)比影響都很小;隨著電子密度降低,E2 躍遷對譜線線強(qiáng)比的影響增強(qiáng),并且E2 躍遷對類鐵鈾離子的影響大于類鐵鉬離子的影響.因此,E2 躍遷對能級布居的影響可能會直接導(dǎo)致電子密度診斷的偏差,所以對于原子序數(shù)較大的原子體系,E2 躍遷通道必須予以考慮.

    圖5 類鐵鉬離子(a)和鈾離子(b)的密度敏感線,圖中的數(shù)字(9/7,···)對應(yīng)于躍遷標(biāo)號,與表1 一致Fig.5.Density-sensitive line radios for Mo16+(a) and U66+(b).The numbers (9/7,···) correspond to the transitions in Table 1.

    4 結(jié)論

    本文從解析模型和大規(guī)模數(shù)值模擬兩方面研究了電子束離子阱等離子體中E2 躍遷對離子能級布居的影響.首先,通過求解三能級原子體系穩(wěn)態(tài)近似下的速率方程得到了離子布居的解析表達(dá)式,并對其在高、低電子密度極限的表達(dá)形式進(jìn)行分析,得到如下結(jié)論: E2 退激通道在高電子密度條件下對離子布居影響很小,在低電子密度條件下對離子布居影響較大.進(jìn)一步,本文選擇類鐵鉬和鈾等離子體,構(gòu)建并數(shù)值求解穩(wěn)態(tài)近似下的碰撞輻射模型,討論了E2 躍遷對類鐵鉬和鈾離子基組態(tài)能級布居和相應(yīng)M1 譜線強(qiáng)度對電子密度依賴關(guān)系的影響.研究結(jié)果表明E2 躍遷通道在高電子密度情況下對兩個離子能級布居和譜線強(qiáng)度的影響都較小,隨著電子密度降低,對鈾離子的影響強(qiáng)于鉬離子.此外本文還分析了E2 躍遷對類鐵鉬離子和鈾離子中若干對電子密度敏感的譜線強(qiáng)度比的影響,并指出對于高Z原子體系,E2 躍遷對于M1 譜線線強(qiáng)比有著不可忽略的貢獻(xiàn).

    猜你喜歡
    電子密度譜線能級
    “拼、搶、快、優(yōu)”,展現(xiàn)錢塘“高能級”擔(dān)當(dāng)
    杭州(2023年3期)2023-04-03 07:22:04
    基于HITRAN光譜數(shù)據(jù)庫的合并譜線測溫仿真研究
    顧及地磁影響的GNSS電離層層析不等像素間距算法*
    提升醫(yī)學(xué)教育能級 培養(yǎng)拔尖創(chuàng)新人才
    糾纏Ξ-型三能級原子與糾纏腔場相互作用熵的糾纏演化
    不同GPS掩星電離層剖面產(chǎn)品相關(guān)性分析
    測繪通報(2019年11期)2019-12-03 01:47:34
    等離子體電子密度分布信息提取方法研究
    鐵合金光譜譜線分離實(shí)驗研究
    電子測試(2018年11期)2018-06-26 05:56:00
    一種適用于電離層電子密度重構(gòu)的AMART算法
    鍶原子光鐘鐘躍遷譜線探測中的程序控制
    精品卡一卡二卡四卡免费| 大陆偷拍与自拍| 午夜av观看不卡| 多毛熟女@视频| 99九九在线精品视频| 丝袜在线中文字幕| 欧美人与善性xxx| 麻豆av在线久日| 久久久国产一区二区| 午夜日韩欧美国产| 男人舔女人的私密视频| videosex国产| 少妇人妻 视频| 亚洲国产精品一区三区| 国产精品一二三区在线看| 久久久精品免费免费高清| 51午夜福利影视在线观看| 国产极品粉嫩免费观看在线| 亚洲欧美日韩高清在线视频 | 最新在线观看一区二区三区 | 欧美黄色片欧美黄色片| 久久性视频一级片| 啦啦啦中文免费视频观看日本| 黄网站色视频无遮挡免费观看| 午夜免费观看性视频| 亚洲一码二码三码区别大吗| 多毛熟女@视频| 国产在线视频一区二区| 99香蕉大伊视频| 丰满迷人的少妇在线观看| 97精品久久久久久久久久精品| 天天躁狠狠躁夜夜躁狠狠躁| 大话2 男鬼变身卡| 国产精品免费大片| 精品少妇一区二区三区视频日本电影| 亚洲第一青青草原| 久9热在线精品视频| 久久人人爽av亚洲精品天堂| 晚上一个人看的免费电影| 国产精品二区激情视频| 亚洲人成电影观看| 欧美大码av| 老司机在亚洲福利影院| 国产精品久久久人人做人人爽| 交换朋友夫妻互换小说| 亚洲国产中文字幕在线视频| 国产片特级美女逼逼视频| 欧美97在线视频| 亚洲欧美一区二区三区黑人| 又粗又硬又长又爽又黄的视频| 久久精品国产综合久久久| 黑人巨大精品欧美一区二区蜜桃| 久久国产精品大桥未久av| 精品久久久久久电影网| av不卡在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 日韩av不卡免费在线播放| 熟女少妇亚洲综合色aaa.| 免费日韩欧美在线观看| av视频免费观看在线观看| av网站免费在线观看视频| 国产精品免费大片| 精品福利观看| 久久av网站| 欧美日韩成人在线一区二区| 国产精品一区二区免费欧美 | 人人妻人人爽人人添夜夜欢视频| 久久性视频一级片| 亚洲精品av麻豆狂野| 精品少妇一区二区三区视频日本电影| 久久亚洲国产成人精品v| 国产精品国产三级专区第一集| 伦理电影免费视频| 国产无遮挡羞羞视频在线观看| 亚洲第一青青草原| 亚洲美女黄色视频免费看| 50天的宝宝边吃奶边哭怎么回事| 国产爽快片一区二区三区| 国产一区二区激情短视频 | 热re99久久国产66热| 成年av动漫网址| 亚洲人成77777在线视频| 性少妇av在线| 国产精品熟女久久久久浪| 日韩 欧美 亚洲 中文字幕| 欧美精品一区二区大全| 久久ye,这里只有精品| 秋霞在线观看毛片| 啦啦啦 在线观看视频| 波多野结衣一区麻豆| 看免费av毛片| 国产日韩欧美亚洲二区| 交换朋友夫妻互换小说| 亚洲精品久久午夜乱码| 久久久国产欧美日韩av| 亚洲黑人精品在线| 久久午夜综合久久蜜桃| 另类精品久久| 视频区欧美日本亚洲| 菩萨蛮人人尽说江南好唐韦庄| 国语对白做爰xxxⅹ性视频网站| 国产深夜福利视频在线观看| 国产成人精品久久久久久| 777米奇影视久久| 国产精品熟女久久久久浪| 国产xxxxx性猛交| 最近最新中文字幕大全免费视频 | 久久久久精品人妻al黑| 国产精品亚洲av一区麻豆| 婷婷色综合www| 国产精品久久久av美女十八| 久久性视频一级片| 一区二区三区乱码不卡18| 国产免费福利视频在线观看| 欧美另类一区| 亚洲av综合色区一区| 岛国毛片在线播放| 色婷婷av一区二区三区视频| 97人妻天天添夜夜摸| 1024香蕉在线观看| 人人妻人人澡人人看| 欧美乱码精品一区二区三区| 一区二区日韩欧美中文字幕| 免费在线观看黄色视频的| 大片免费播放器 马上看| 亚洲成人免费av在线播放| 在线看a的网站| 免费av中文字幕在线| 久久久精品区二区三区| 夫妻性生交免费视频一级片| 久久久久视频综合| 国产亚洲一区二区精品| 国产精品一区二区在线观看99| 一本—道久久a久久精品蜜桃钙片| 久久精品亚洲熟妇少妇任你| 国产91精品成人一区二区三区 | 日韩人妻精品一区2区三区| 精品一区在线观看国产| 丝瓜视频免费看黄片| 性色av乱码一区二区三区2| 美女扒开内裤让男人捅视频| 天堂俺去俺来也www色官网| www.av在线官网国产| 一级毛片 在线播放| 老司机在亚洲福利影院| 亚洲成av片中文字幕在线观看| 欧美日韩福利视频一区二区| 美女国产高潮福利片在线看| 在线天堂中文资源库| 日韩一卡2卡3卡4卡2021年| 亚洲一区中文字幕在线| 国产伦理片在线播放av一区| 18禁黄网站禁片午夜丰满| 欧美人与性动交α欧美软件| 精品少妇黑人巨大在线播放| 人人妻人人澡人人看| 久久精品亚洲熟妇少妇任你| 在线天堂中文资源库| 亚洲精品乱久久久久久| 美女视频免费永久观看网站| 久久精品国产亚洲av高清一级| 国产亚洲午夜精品一区二区久久| 欧美日韩亚洲高清精品| 永久免费av网站大全| 黄色视频不卡| 日本猛色少妇xxxxx猛交久久| 黄网站色视频无遮挡免费观看| 国产亚洲欧美在线一区二区| 亚洲伊人久久精品综合| 777久久人妻少妇嫩草av网站| 巨乳人妻的诱惑在线观看| 国产精品 国内视频| 亚洲国产精品一区三区| 国产成人系列免费观看| 别揉我奶头~嗯~啊~动态视频 | 精品亚洲成国产av| 久久性视频一级片| 国产片特级美女逼逼视频| 中文字幕人妻丝袜一区二区| 极品人妻少妇av视频| 国产精品久久久久久精品电影小说| 欧美人与善性xxx| 国产97色在线日韩免费| 国产日韩欧美视频二区| 黄频高清免费视频| www.精华液| 在线亚洲精品国产二区图片欧美| 电影成人av| 亚洲精品国产一区二区精华液| 我要看黄色一级片免费的| 久久亚洲精品不卡| 亚洲国产欧美一区二区综合| 青春草亚洲视频在线观看| 亚洲国产精品国产精品| 久久性视频一级片| 午夜福利影视在线免费观看| 亚洲精品一二三| 精品久久久久久电影网| 蜜桃国产av成人99| 国产有黄有色有爽视频| 两个人看的免费小视频| 日韩视频在线欧美| 亚洲国产av影院在线观看| 亚洲欧美精品综合一区二区三区| 叶爱在线成人免费视频播放| 日本av手机在线免费观看| 欧美av亚洲av综合av国产av| 十八禁人妻一区二区| 亚洲av片天天在线观看| av欧美777| 国产国语露脸激情在线看| 水蜜桃什么品种好| 一级毛片黄色毛片免费观看视频| 亚洲 欧美一区二区三区| 亚洲第一av免费看| 精品一区二区三区四区五区乱码 | 精品免费久久久久久久清纯 | 久久人妻熟女aⅴ| 9色porny在线观看| 成年动漫av网址| 婷婷色综合大香蕉| 久久人人爽人人片av| 成人亚洲欧美一区二区av| 精品一区在线观看国产| 欧美97在线视频| 国产欧美日韩综合在线一区二区| 精品少妇一区二区三区视频日本电影| 久久精品国产亚洲av高清一级| 啦啦啦中文免费视频观看日本| 成年动漫av网址| 少妇人妻 视频| 亚洲精品成人av观看孕妇| 黑丝袜美女国产一区| 中文字幕最新亚洲高清| 午夜福利视频在线观看免费| 婷婷色综合大香蕉| 丰满饥渴人妻一区二区三| 狠狠婷婷综合久久久久久88av| 两个人看的免费小视频| 午夜福利一区二区在线看| 精品久久久久久电影网| 国产欧美日韩一区二区三区在线| 日韩av在线免费看完整版不卡| 亚洲国产最新在线播放| 五月开心婷婷网| 精品一区在线观看国产| 精品视频人人做人人爽| 成年女人毛片免费观看观看9 | 国产精品偷伦视频观看了| 99九九在线精品视频| 欧美黑人欧美精品刺激| 一级毛片电影观看| 国产又爽黄色视频| 国产成人91sexporn| 亚洲精品国产区一区二| 嫁个100分男人电影在线观看 | 欧美精品一区二区免费开放| 久久久久国产精品人妻一区二区| 久久99一区二区三区| 99久久综合免费| 18禁国产床啪视频网站| 国产精品人妻久久久影院| 中文字幕制服av| 久久天躁狠狠躁夜夜2o2o | 婷婷色麻豆天堂久久| 夜夜骑夜夜射夜夜干| 久久鲁丝午夜福利片| 男女边吃奶边做爰视频| 天天影视国产精品| 最新在线观看一区二区三区 | 欧美日韩成人在线一区二区| 亚洲一区中文字幕在线| 久久人妻熟女aⅴ| 亚洲 欧美一区二区三区| 国产午夜精品一二区理论片| av又黄又爽大尺度在线免费看| 自线自在国产av| 大香蕉久久成人网| 男人爽女人下面视频在线观看| 看十八女毛片水多多多| 高清不卡的av网站| 中文字幕人妻丝袜一区二区| 久久午夜综合久久蜜桃| 午夜免费男女啪啪视频观看| 国产精品免费大片| 欧美在线黄色| 国产男人的电影天堂91| 丝袜喷水一区| 国产一区二区在线观看av| 国产一级毛片在线| 波多野结衣一区麻豆| 一级黄色大片毛片| 母亲3免费完整高清在线观看| 日韩免费高清中文字幕av| 国产欧美日韩精品亚洲av| 国产高清videossex| 国产一区有黄有色的免费视频| 中文字幕人妻熟女乱码| 丝袜在线中文字幕| 日本午夜av视频| 国产精品免费视频内射| 男女边摸边吃奶| 99国产精品一区二区三区| 久久女婷五月综合色啪小说| 亚洲 欧美一区二区三区| 亚洲视频免费观看视频| 男女国产视频网站| avwww免费| 亚洲av成人精品一二三区| 精品国产超薄肉色丝袜足j| 亚洲一码二码三码区别大吗| 超碰成人久久| 免费观看人在逋| 精品国产国语对白av| 国产一区二区激情短视频 | 国产成人欧美| 成人手机av| 欧美日韩精品网址| 天天躁夜夜躁狠狠久久av| 欧美精品啪啪一区二区三区 | 亚洲av国产av综合av卡| 久久国产精品影院| 丝袜人妻中文字幕| 亚洲中文字幕日韩| 一二三四社区在线视频社区8| 国产精品欧美亚洲77777| 国产真人三级小视频在线观看| 久久久久久亚洲精品国产蜜桃av| 久久久国产欧美日韩av| 久久久久视频综合| 久久国产精品男人的天堂亚洲| 亚洲伊人久久精品综合| 国产日韩欧美亚洲二区| 国产精品亚洲av一区麻豆| 精品亚洲乱码少妇综合久久| 操美女的视频在线观看| 无遮挡黄片免费观看| 国产成人91sexporn| 黑人巨大精品欧美一区二区蜜桃| 菩萨蛮人人尽说江南好唐韦庄| 亚洲专区中文字幕在线| 99热国产这里只有精品6| 国产无遮挡羞羞视频在线观看| 十八禁人妻一区二区| 一级黄色大片毛片| 大香蕉久久成人网| 国产精品久久久久久精品古装| 永久免费av网站大全| 欧美人与性动交α欧美软件| 在线看a的网站| 免费日韩欧美在线观看| 亚洲七黄色美女视频| 亚洲精品国产av成人精品| 久久中文字幕一级| 高清欧美精品videossex| 99精国产麻豆久久婷婷| 老司机深夜福利视频在线观看 | 狂野欧美激情性xxxx| 免费黄频网站在线观看国产| 91麻豆av在线| 精品国产一区二区三区四区第35| 美女扒开内裤让男人捅视频| 亚洲精品久久成人aⅴ小说| 免费久久久久久久精品成人欧美视频| 午夜福利一区二区在线看| 亚洲欧洲日产国产| 亚洲成人国产一区在线观看 | 黑丝袜美女国产一区| 亚洲成人手机| 国产av精品麻豆| 中文字幕高清在线视频| 成在线人永久免费视频| 久久精品国产综合久久久| 国产精品一区二区在线观看99| 亚洲中文日韩欧美视频| 黑丝袜美女国产一区| 欧美日韩国产mv在线观看视频| 精品久久久久久久毛片微露脸 | 午夜激情av网站| 性色av乱码一区二区三区2| 国产av国产精品国产| 午夜免费观看性视频| 亚洲色图 男人天堂 中文字幕| 人人妻人人澡人人看| 一本色道久久久久久精品综合| 好男人电影高清在线观看| 国产亚洲精品久久久久5区| 99久久综合免费| 欧美激情 高清一区二区三区| e午夜精品久久久久久久| 午夜福利影视在线免费观看| 久久人人爽人人片av| 午夜福利视频在线观看免费| 久久这里只有精品19| 免费看av在线观看网站| 国产黄色免费在线视频| 亚洲国产精品成人久久小说| 在线看a的网站| 不卡av一区二区三区| 久久精品人人爽人人爽视色| 黑人欧美特级aaaaaa片| av视频免费观看在线观看| av片东京热男人的天堂| 日韩一区二区三区影片| 国产在线免费精品| 色视频在线一区二区三区| 美女中出高潮动态图| 日本wwww免费看| 性色av一级| 交换朋友夫妻互换小说| 欧美日韩国产mv在线观看视频| 一二三四社区在线视频社区8| 夫妻午夜视频| 欧美精品高潮呻吟av久久| 99精品久久久久人妻精品| 亚洲精品成人av观看孕妇| 狠狠婷婷综合久久久久久88av| 国产1区2区3区精品| 国产精品国产av在线观看| 各种免费的搞黄视频| www.自偷自拍.com| 波多野结衣一区麻豆| 99热全是精品| 交换朋友夫妻互换小说| 一本大道久久a久久精品| 久久久国产精品麻豆| 黄片播放在线免费| 欧美黑人精品巨大| 少妇的丰满在线观看| 91字幕亚洲| 中文欧美无线码| 国产日韩欧美视频二区| 交换朋友夫妻互换小说| 亚洲精品国产av成人精品| 国产高清不卡午夜福利| 日日爽夜夜爽网站| 午夜福利一区二区在线看| 各种免费的搞黄视频| 嫩草影视91久久| 麻豆国产av国片精品| 国产老妇伦熟女老妇高清| 日韩熟女老妇一区二区性免费视频| 国产成人一区二区三区免费视频网站 | 黄频高清免费视频| 久久国产精品人妻蜜桃| 国产精品香港三级国产av潘金莲 | 中文字幕亚洲精品专区| 母亲3免费完整高清在线观看| av国产久精品久网站免费入址| 欧美黑人欧美精品刺激| 欧美精品一区二区大全| 亚洲一区中文字幕在线| 久久九九热精品免费| 三上悠亚av全集在线观看| 日韩伦理黄色片| 亚洲国产成人一精品久久久| 视频在线观看一区二区三区| 欧美日韩av久久| 黄频高清免费视频| 免费一级毛片在线播放高清视频 | 欧美日韩福利视频一区二区| 又紧又爽又黄一区二区| 伊人亚洲综合成人网| av电影中文网址| 女人久久www免费人成看片| 久久国产亚洲av麻豆专区| 国产伦人伦偷精品视频| 国产精品国产三级国产专区5o| 欧美人与性动交α欧美精品济南到| 成人国产一区最新在线观看 | 亚洲视频免费观看视频| av国产精品久久久久影院| 人人妻人人澡人人爽人人夜夜| 国产在视频线精品| 男女边摸边吃奶| 国产成人欧美| 成年女人毛片免费观看观看9 | 2018国产大陆天天弄谢| www.999成人在线观看| 亚洲国产最新在线播放| 麻豆av在线久日| 亚洲国产成人一精品久久久| 国产精品久久久久久精品古装| 亚洲人成电影免费在线| 巨乳人妻的诱惑在线观看| 丰满迷人的少妇在线观看| 人人妻人人澡人人爽人人夜夜| 美女国产高潮福利片在线看| 日本色播在线视频| 午夜影院在线不卡| 欧美日韩视频高清一区二区三区二| 成人亚洲欧美一区二区av| 国产精品99久久99久久久不卡| 成年人黄色毛片网站| 男女午夜视频在线观看| 欧美大码av| 人人妻人人澡人人爽人人夜夜| 曰老女人黄片| 欧美日韩福利视频一区二区| 美女中出高潮动态图| 亚洲欧美日韩高清在线视频 | 亚洲中文av在线| 精品高清国产在线一区| 夫妻性生交免费视频一级片| 首页视频小说图片口味搜索 | 精品少妇黑人巨大在线播放| 丰满迷人的少妇在线观看| 久久精品aⅴ一区二区三区四区| 午夜激情久久久久久久| 精品国产一区二区三区四区第35| 99久久精品国产亚洲精品| 中文字幕亚洲精品专区| 又大又爽又粗| 久久久久国产精品人妻一区二区| 视频在线观看一区二区三区| 曰老女人黄片| 成年美女黄网站色视频大全免费| 亚洲欧美清纯卡通| 久久中文字幕一级| 我要看黄色一级片免费的| 少妇裸体淫交视频免费看高清 | 观看av在线不卡| 国产日韩欧美视频二区| 中国美女看黄片| 亚洲五月色婷婷综合| 国产黄色免费在线视频| 久久久精品区二区三区| 波多野结衣一区麻豆| 国产爽快片一区二区三区| 日韩一卡2卡3卡4卡2021年| 91成人精品电影| 亚洲精品乱久久久久久| 精品视频人人做人人爽| 欧美日韩亚洲国产一区二区在线观看 | 国产免费一区二区三区四区乱码| 精品久久久久久久毛片微露脸 | 欧美少妇被猛烈插入视频| 久久久欧美国产精品| 国产成人av激情在线播放| 在线天堂中文资源库| 热re99久久国产66热| 亚洲九九香蕉| 五月天丁香电影| 亚洲精品自拍成人| 亚洲中文字幕日韩| 中文欧美无线码| 国产欧美日韩一区二区三区在线| 精品国产一区二区久久| 十八禁高潮呻吟视频| 国产精品人妻久久久影院| 在线观看免费视频网站a站| 啦啦啦在线免费观看视频4| 中文字幕另类日韩欧美亚洲嫩草| 在线观看www视频免费| 国产一区亚洲一区在线观看| 成年人免费黄色播放视频| 美女国产高潮福利片在线看| 亚洲欧洲国产日韩| 国产一卡二卡三卡精品| 国产精品久久久久成人av| 欧美人与性动交α欧美精品济南到| 欧美97在线视频| 少妇猛男粗大的猛烈进出视频| 亚洲国产欧美一区二区综合| 中文字幕亚洲精品专区| svipshipincom国产片| 天天添夜夜摸| 97人妻天天添夜夜摸| 老司机午夜十八禁免费视频| 欧美黄色片欧美黄色片| 五月天丁香电影| 操出白浆在线播放| 午夜福利一区二区在线看| 久久人人爽av亚洲精品天堂| 欧美日本中文国产一区发布| 男人爽女人下面视频在线观看| 亚洲欧洲国产日韩| 久久亚洲国产成人精品v| 不卡av一区二区三区| 精品熟女少妇八av免费久了| 国产精品国产av在线观看| 真人做人爱边吃奶动态| 免费看十八禁软件| 久久久久久久大尺度免费视频| 最新在线观看一区二区三区 | 熟女少妇亚洲综合色aaa.| 午夜福利免费观看在线| 日韩一区二区三区影片| 婷婷丁香在线五月| 久久久国产精品麻豆| 天堂俺去俺来也www色官网| 啦啦啦在线观看免费高清www| 国产一区二区激情短视频 | 高清视频免费观看一区二区| 成年人午夜在线观看视频| 国产麻豆69| 亚洲成av片中文字幕在线观看| 女性被躁到高潮视频| 女警被强在线播放| 亚洲,欧美精品.| 美女国产高潮福利片在线看| 老司机午夜十八禁免费视频| 99国产精品免费福利视频| 国产精品一区二区精品视频观看| 国产97色在线日韩免费| 一二三四社区在线视频社区8| 曰老女人黄片| 国产麻豆69| 久久精品国产亚洲av涩爱| 一二三四在线观看免费中文在| 首页视频小说图片口味搜索 | 日韩免费高清中文字幕av| 一区在线观看完整版| 欧美黑人精品巨大| 一级毛片我不卡|