• <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算法
    鍶原子光鐘鐘躍遷譜線探測中的程序控制
    校园春色视频在线观看| 麻豆成人午夜福利视频| www.999成人在线观看| 午夜福利欧美成人| 在线视频色国产色| 精品久久久久久,| 嫩草影视91久久| 全区人妻精品视频| 男女下面进入的视频免费午夜| 天堂√8在线中文| 久久久久国内视频| 美女扒开内裤让男人捅视频| 亚洲午夜精品一区,二区,三区| 人人妻,人人澡人人爽秒播| 亚洲熟妇中文字幕五十中出| 午夜精品在线福利| 国内精品久久久久精免费| 在线a可以看的网站| 97人妻精品一区二区三区麻豆| 视频区欧美日本亚洲| 久久性视频一级片| 欧美中文综合在线视频| 国产精品一区二区精品视频观看| 两个人免费观看高清视频| 丝袜美腿诱惑在线| 国内精品久久久久久久电影| 99久久99久久久精品蜜桃| 每晚都被弄得嗷嗷叫到高潮| 日本精品一区二区三区蜜桃| 最近最新中文字幕大全免费视频| 少妇的丰满在线观看| 在线观看舔阴道视频| 18禁观看日本| 九九热线精品视视频播放| 好男人在线观看高清免费视频| 级片在线观看| 在线观看免费视频日本深夜| 美女高潮喷水抽搐中文字幕| 1024手机看黄色片| 免费观看精品视频网站| 18禁观看日本| 叶爱在线成人免费视频播放| 国产成人欧美在线观看| 岛国在线观看网站| 99re在线观看精品视频| 一个人观看的视频www高清免费观看 | 久久香蕉精品热| 两个人免费观看高清视频| 亚洲国产精品成人综合色| 欧美3d第一页| 欧美精品啪啪一区二区三区| 18禁黄网站禁片免费观看直播| av在线天堂中文字幕| 久久久精品国产亚洲av高清涩受| 9191精品国产免费久久| 欧美av亚洲av综合av国产av| 国产黄色小视频在线观看| 99久久国产精品久久久| 亚洲免费av在线视频| 亚洲av中文字字幕乱码综合| 手机成人av网站| 99国产极品粉嫩在线观看| 色综合站精品国产| 久久久久久亚洲精品国产蜜桃av| 婷婷亚洲欧美| 亚洲激情在线av| a级毛片在线看网站| 亚洲avbb在线观看| 久久国产精品影院| 日韩欧美三级三区| 操出白浆在线播放| 女人爽到高潮嗷嗷叫在线视频| 亚洲18禁久久av| 在线观看美女被高潮喷水网站 | 亚洲成人免费电影在线观看| 久久久久久久久久黄片| av福利片在线| 国产99久久九九免费精品| 欧美日韩黄片免| 99精品欧美一区二区三区四区| a在线观看视频网站| 亚洲欧洲精品一区二区精品久久久| 午夜激情av网站| 日本 av在线| 一级片免费观看大全| 99热这里只有精品一区 | 老司机靠b影院| 国产成年人精品一区二区| 国产精华一区二区三区| 亚洲人成网站高清观看| 久久精品人妻少妇| 性色av乱码一区二区三区2| 欧美中文日本在线观看视频| 欧美精品啪啪一区二区三区| 欧美色视频一区免费| 午夜福利成人在线免费观看| 亚洲人成网站在线播放欧美日韩| 亚洲精品粉嫩美女一区| 国产成人aa在线观看| 99国产精品99久久久久| 国产成人系列免费观看| 亚洲精品久久国产高清桃花| 脱女人内裤的视频| 亚洲真实伦在线观看| 一个人免费在线观看的高清视频| 97超级碰碰碰精品色视频在线观看| 91九色精品人成在线观看| 亚洲精品一区av在线观看| 国产精品自产拍在线观看55亚洲| 成熟少妇高潮喷水视频| 成人特级黄色片久久久久久久| 日韩欧美三级三区| 国产成人aa在线观看| 午夜日韩欧美国产| 国产又色又爽无遮挡免费看| 91国产中文字幕| 精品一区二区三区av网在线观看| 久久久久亚洲av毛片大全| 嫩草影院精品99| 国产视频内射| 黄色成人免费大全| 久久天堂一区二区三区四区| 国产亚洲精品久久久久5区| 岛国在线免费视频观看| 免费av毛片视频| 国产高清视频在线观看网站| 国产精品久久久久久人妻精品电影| 亚洲欧美激情综合另类| 99热这里只有是精品50| 国产熟女xx| 亚洲全国av大片| 欧美黑人欧美精品刺激| 母亲3免费完整高清在线观看| 久久久久九九精品影院| 日韩欧美国产在线观看| 国语自产精品视频在线第100页| 午夜福利高清视频| 一a级毛片在线观看| 欧美日韩亚洲综合一区二区三区_| 亚洲成人久久爱视频| 国产97色在线日韩免费| 精品久久久久久久久久久久久| 国产1区2区3区精品| 白带黄色成豆腐渣| 老汉色∧v一级毛片| 少妇人妻一区二区三区视频| av在线天堂中文字幕| 亚洲人成伊人成综合网2020| 一进一出抽搐动态| 国产v大片淫在线免费观看| 99热这里只有是精品50| 中文字幕av在线有码专区| 777久久人妻少妇嫩草av网站| 亚洲国产精品999在线| 日日干狠狠操夜夜爽| 欧美中文日本在线观看视频| www日本在线高清视频| 欧美在线黄色| 日韩大码丰满熟妇| 亚洲av中文字字幕乱码综合| 两个人免费观看高清视频| 国产亚洲av嫩草精品影院| 天堂av国产一区二区熟女人妻 | 极品教师在线免费播放| 国产又黄又爽又无遮挡在线| 欧美日本亚洲视频在线播放| 男人的好看免费观看在线视频 | 51午夜福利影视在线观看| 国产人伦9x9x在线观看| 国产男靠女视频免费网站| 国产成人精品无人区| 九色国产91popny在线| 99国产极品粉嫩在线观看| 母亲3免费完整高清在线观看| 一级毛片女人18水好多| 婷婷六月久久综合丁香| 美女 人体艺术 gogo| 亚洲色图av天堂| 亚洲精品美女久久av网站| 两人在一起打扑克的视频| 婷婷精品国产亚洲av| 99在线人妻在线中文字幕| 老司机午夜福利在线观看视频| 亚洲熟女毛片儿| or卡值多少钱| 最新美女视频免费是黄的| 男人舔女人的私密视频| 一本久久中文字幕| 天天一区二区日本电影三级| 色综合站精品国产| 欧美久久黑人一区二区| 欧美在线黄色| 老司机午夜十八禁免费视频| 国产午夜福利久久久久久| 99国产综合亚洲精品| 国产男靠女视频免费网站| 两个人视频免费观看高清| 两个人看的免费小视频| 亚洲真实伦在线观看| 在线观看午夜福利视频| 亚洲最大成人中文| 国产精品一区二区三区四区久久| 成人亚洲精品av一区二区| 国产精品免费视频内射| 一个人免费在线观看的高清视频| 99热这里只有精品一区 | 成年女人毛片免费观看观看9| 欧美高清成人免费视频www| 日本熟妇午夜| 中文字幕人妻丝袜一区二区| 国产成人欧美在线观看| 亚洲午夜精品一区,二区,三区| 亚洲自拍偷在线| 国产精华一区二区三区| 日日爽夜夜爽网站| 人人妻人人澡欧美一区二区| 九九热线精品视视频播放| 国内久久婷婷六月综合欲色啪| 又紧又爽又黄一区二区| 美女高潮喷水抽搐中文字幕| 视频区欧美日本亚洲| 亚洲人成电影免费在线| 一级毛片高清免费大全| 老熟妇仑乱视频hdxx| 99国产精品一区二区蜜桃av| 久久久久精品国产欧美久久久| 小说图片视频综合网站| 精品乱码久久久久久99久播| 老司机福利观看| 国产精品国产高清国产av| 亚洲美女视频黄频| 禁无遮挡网站| 老熟妇乱子伦视频在线观看| 1024香蕉在线观看| 欧美日韩一级在线毛片| 在线国产一区二区在线| 一二三四社区在线视频社区8| 麻豆国产av国片精品| 一级毛片女人18水好多| 亚洲精品一区av在线观看| 麻豆国产97在线/欧美 | 国产亚洲精品av在线| 又粗又爽又猛毛片免费看| 欧美日韩中文字幕国产精品一区二区三区| 在线看三级毛片| 看片在线看免费视频| 国产主播在线观看一区二区| 91国产中文字幕| 亚洲欧美日韩高清在线视频| 不卡一级毛片| 国产蜜桃级精品一区二区三区| 国产高清视频在线观看网站| 天天添夜夜摸| 蜜桃久久精品国产亚洲av| 校园春色视频在线观看| 精品少妇一区二区三区视频日本电影| 最近视频中文字幕2019在线8| 99久久精品国产亚洲精品| 中文资源天堂在线| 亚洲国产欧美人成| 欧美中文日本在线观看视频| 人人妻,人人澡人人爽秒播| 黄色视频不卡| 狂野欧美激情性xxxx| 香蕉久久夜色| or卡值多少钱| 午夜福利成人在线免费观看| 99国产极品粉嫩在线观看| 亚洲av第一区精品v没综合| 真人做人爱边吃奶动态| 亚洲成人中文字幕在线播放| xxxwww97欧美| 在线观看免费视频日本深夜| 日本 欧美在线| 美女扒开内裤让男人捅视频| 免费在线观看影片大全网站| 性色av乱码一区二区三区2| 精品日产1卡2卡| 日本成人三级电影网站| 欧美在线一区亚洲| 久久久精品大字幕| 少妇被粗大的猛进出69影院| 美女 人体艺术 gogo| 国产一级毛片七仙女欲春2| 国产一区二区在线av高清观看| 毛片女人毛片| av福利片在线观看| 女生性感内裤真人,穿戴方法视频| 国产精品 欧美亚洲| 身体一侧抽搐| www.www免费av| 日日摸夜夜添夜夜添小说| 91字幕亚洲| 亚洲五月婷婷丁香| 国产亚洲精品av在线| 97人妻精品一区二区三区麻豆| 美女 人体艺术 gogo| 久久久久免费精品人妻一区二区| 男女下面进入的视频免费午夜| 一边摸一边抽搐一进一小说| 亚洲专区字幕在线| 日韩精品免费视频一区二区三区| 亚洲欧美日韩高清专用| 午夜久久久久精精品| 免费看a级黄色片| 国产成人精品久久二区二区91| 蜜桃久久精品国产亚洲av| 一二三四在线观看免费中文在| 1024香蕉在线观看| av福利片在线| 搞女人的毛片| 亚洲精品色激情综合| 欧美极品一区二区三区四区| 动漫黄色视频在线观看| 精品久久久久久久毛片微露脸| 国产av又大| 欧美成人午夜精品| 最近视频中文字幕2019在线8| 国产单亲对白刺激| 国产亚洲精品一区二区www| 一二三四社区在线视频社区8| 很黄的视频免费| 最近视频中文字幕2019在线8| 激情在线观看视频在线高清| 日韩欧美 国产精品| 两个人免费观看高清视频| 国产亚洲av嫩草精品影院| 一本精品99久久精品77| av福利片在线观看| 亚洲精品中文字幕在线视频| 国产成年人精品一区二区| 国内揄拍国产精品人妻在线| 在线观看日韩欧美| xxx96com| 国产精品久久久久久久电影 | 成人三级黄色视频| 久久久久国内视频| 国产精品一区二区免费欧美| 国产一区在线观看成人免费| 99精品久久久久人妻精品| 在线a可以看的网站| 国产欧美日韩精品亚洲av| 午夜精品久久久久久毛片777| 三级男女做爰猛烈吃奶摸视频| 国产成年人精品一区二区| 丰满人妻熟妇乱又伦精品不卡| 久久久久久免费高清国产稀缺| 欧美日韩中文字幕国产精品一区二区三区| 亚洲人成网站高清观看| 日韩成人在线观看一区二区三区| 久久久久九九精品影院| 亚洲人成网站在线播放欧美日韩| 90打野战视频偷拍视频| 丰满人妻熟妇乱又伦精品不卡| 久久精品91蜜桃| 桃色一区二区三区在线观看| 午夜激情av网站| 19禁男女啪啪无遮挡网站| 精品熟女少妇八av免费久了| 国产一级毛片七仙女欲春2| 90打野战视频偷拍视频| 国产又黄又爽又无遮挡在线| 亚洲片人在线观看| 老司机午夜十八禁免费视频| 亚洲国产日韩欧美精品在线观看 | 日韩中文字幕欧美一区二区| 天天添夜夜摸| 很黄的视频免费| 91字幕亚洲| 老鸭窝网址在线观看| 成人一区二区视频在线观看| 中文字幕人成人乱码亚洲影| 成熟少妇高潮喷水视频| 国产精品综合久久久久久久免费| 99久久无色码亚洲精品果冻| 日日摸夜夜添夜夜添小说| 特大巨黑吊av在线直播| 欧美日韩黄片免| 久久国产乱子伦精品免费另类| 久久中文看片网| 91成年电影在线观看| 国产高清激情床上av| 丁香欧美五月| 久久精品国产清高在天天线| 国产一区二区激情短视频| av超薄肉色丝袜交足视频| 99久久无色码亚洲精品果冻| 午夜激情福利司机影院| 正在播放国产对白刺激| 特级一级黄色大片| 动漫黄色视频在线观看| 亚洲真实伦在线观看| 男人舔女人下体高潮全视频| 天天躁狠狠躁夜夜躁狠狠躁| 99re在线观看精品视频| 亚洲精品美女久久久久99蜜臀| av国产免费在线观看| 国产精品,欧美在线| 日韩av在线大香蕉| 可以在线观看毛片的网站| 精品久久久久久,| 成人手机av| 日本在线视频免费播放| 国产免费av片在线观看野外av| 国内精品久久久久久久电影| av福利片在线| 激情在线观看视频在线高清| 色精品久久人妻99蜜桃| 亚洲中文日韩欧美视频| 国产精品一及| 久久人妻av系列| 亚洲国产精品sss在线观看| 国内精品久久久久久久电影| 国内少妇人妻偷人精品xxx网站 | 嫁个100分男人电影在线观看| 久久久久久久久中文| 欧美黑人欧美精品刺激| 精品欧美一区二区三区在线| 国产精品香港三级国产av潘金莲| 三级毛片av免费| 亚洲全国av大片| 色精品久久人妻99蜜桃| 三级国产精品欧美在线观看 | 日韩av在线大香蕉| 在线观看66精品国产| 麻豆国产97在线/欧美 | 亚洲av成人不卡在线观看播放网| bbb黄色大片| 精品少妇一区二区三区视频日本电影| 久久久久久国产a免费观看| 九色国产91popny在线| 午夜影院日韩av| 可以免费在线观看a视频的电影网站| 桃色一区二区三区在线观看| 亚洲精品中文字幕在线视频| 亚洲人成电影免费在线| 在线观看免费视频日本深夜| 91成年电影在线观看| 啦啦啦观看免费观看视频高清| 九色成人免费人妻av| 天堂影院成人在线观看| 亚洲真实伦在线观看| 精品日产1卡2卡| 久久久久久九九精品二区国产 | 舔av片在线| 亚洲国产欧洲综合997久久,| 亚洲精品国产一区二区精华液| 变态另类成人亚洲欧美熟女| www.www免费av| 亚洲国产日韩欧美精品在线观看 | 国产精品久久久久久久电影 | 国产成人av激情在线播放| 91麻豆精品激情在线观看国产| 一本一本综合久久| 欧美最黄视频在线播放免费| 欧美黑人精品巨大| 黄片大片在线免费观看| 看免费av毛片| 免费电影在线观看免费观看| 日本撒尿小便嘘嘘汇集6| 无遮挡黄片免费观看| 久久久国产成人免费| 亚洲专区国产一区二区| 亚洲免费av在线视频| 午夜福利在线在线| 国产av又大| 欧美激情久久久久久爽电影| 天天添夜夜摸| www日本在线高清视频| 99精品在免费线老司机午夜| 亚洲美女视频黄频| 性欧美人与动物交配| 美女免费视频网站| 非洲黑人性xxxx精品又粗又长| 熟女电影av网| 亚洲七黄色美女视频| 手机成人av网站| 久久中文字幕一级| 免费在线观看日本一区| 嫩草影院精品99| 50天的宝宝边吃奶边哭怎么回事| 国产一区二区在线观看日韩 | 欧美乱码精品一区二区三区| 午夜a级毛片| 制服人妻中文乱码| 久久这里只有精品中国| 亚洲欧美激情综合另类| 一区福利在线观看| 亚洲人与动物交配视频| 久久婷婷成人综合色麻豆| 精品无人区乱码1区二区| 亚洲电影在线观看av| 人成视频在线观看免费观看| 欧美乱色亚洲激情| 夜夜爽天天搞| 国产伦在线观看视频一区| 国产精品,欧美在线| 少妇的丰满在线观看| 国产私拍福利视频在线观看| 色av中文字幕| 在线观看66精品国产| 亚洲色图 男人天堂 中文字幕| 在线a可以看的网站| av视频在线观看入口| 亚洲午夜理论影院| 国产av不卡久久| 又黄又粗又硬又大视频| 91成年电影在线观看| 男女那种视频在线观看| 俄罗斯特黄特色一大片| 亚洲av片天天在线观看| 色综合站精品国产| 成年版毛片免费区| 国产高清有码在线观看视频 | 曰老女人黄片| 19禁男女啪啪无遮挡网站| 99久久综合精品五月天人人| 欧美+亚洲+日韩+国产| 制服诱惑二区| 亚洲成人免费电影在线观看| 欧美色视频一区免费| 一区二区三区高清视频在线| 久久久久国产一级毛片高清牌| 人妻夜夜爽99麻豆av| 免费看日本二区| 亚洲男人天堂网一区| 国产人伦9x9x在线观看| 中文字幕熟女人妻在线| 波多野结衣巨乳人妻| av视频在线观看入口| www日本在线高清视频| 精品欧美国产一区二区三| 夜夜看夜夜爽夜夜摸| a级毛片在线看网站| 99久久久亚洲精品蜜臀av| 久久性视频一级片| 国产av在哪里看| 国产成人一区二区三区免费视频网站| 亚洲精品久久成人aⅴ小说| 法律面前人人平等表现在哪些方面| 日韩三级视频一区二区三区| 国产单亲对白刺激| 精品国产亚洲在线| 亚洲一卡2卡3卡4卡5卡精品中文| 无遮挡黄片免费观看| 性欧美人与动物交配| 亚洲精品色激情综合| 黄色a级毛片大全视频| 精品免费久久久久久久清纯| 欧美绝顶高潮抽搐喷水| 欧美性猛交黑人性爽| 一二三四在线观看免费中文在| 欧美成人一区二区免费高清观看 | 三级男女做爰猛烈吃奶摸视频| 1024手机看黄色片| 97人妻精品一区二区三区麻豆| 欧美另类亚洲清纯唯美| 国产亚洲欧美在线一区二区| 亚洲欧美精品综合一区二区三区| 又大又爽又粗| 村上凉子中文字幕在线| 在线看三级毛片| 精品乱码久久久久久99久播| 国产熟女xx| 两个人视频免费观看高清| 亚洲中文字幕一区二区三区有码在线看 | 午夜福利在线在线| 国产69精品久久久久777片 | 国产精品免费一区二区三区在线| 成人精品一区二区免费| 久久久久国内视频| 国产精品一区二区精品视频观看| x7x7x7水蜜桃| 又粗又爽又猛毛片免费看| 国内久久婷婷六月综合欲色啪| 草草在线视频免费看| 黑人操中国人逼视频| 天天添夜夜摸| 日韩精品中文字幕看吧| 亚洲国产欧美一区二区综合| 久久香蕉精品热| 少妇的丰满在线观看| 色噜噜av男人的天堂激情| 精品高清国产在线一区| 久久人妻福利社区极品人妻图片| 色哟哟哟哟哟哟| bbb黄色大片| 成年版毛片免费区| 亚洲性夜色夜夜综合| 国产av一区在线观看免费| 亚洲无线在线观看| 一级毛片精品| 90打野战视频偷拍视频| 日本黄色视频三级网站网址| 麻豆成人午夜福利视频| 国产主播在线观看一区二区| 久久 成人 亚洲| 日韩欧美免费精品| 欧美日韩黄片免| 国产熟女午夜一区二区三区| 婷婷六月久久综合丁香| 国产精品永久免费网站| 我要搜黄色片| 一二三四在线观看免费中文在| 久久久久九九精品影院| 啦啦啦免费观看视频1| 国产激情久久老熟女| 亚洲性夜色夜夜综合| 欧美日本视频| 最近在线观看免费完整版| 亚洲中文字幕日韩| 欧美成人免费av一区二区三区| 五月伊人婷婷丁香| 国产激情偷乱视频一区二区| 99国产精品一区二区蜜桃av|