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

    外擴(kuò)型電磁場控制筒形陰極內(nèi)等離子體放電輸運(yùn)特性的仿真研究*

    2019-10-22 02:01:28崔歲寒吳忠振肖舒陳磊李體軍劉亮亮傅勁裕田修波朱劍豪譚文長
    物理學(xué)報(bào) 2019年19期
    關(guān)鍵詞:靶面陰極等離子體

    崔歲寒 吳忠振? 肖舒 陳磊 李體軍 劉亮亮 傅勁裕 田修波 朱劍豪 譚文長

    1) (北京大學(xué)深圳研究生院新材料學(xué)院, 深圳 518055)

    2) (香港城市大學(xué)物理與材料科學(xué)系, 香港 999077)

    筒形陰極由于具有向內(nèi)放電的特性, 可改善高功率脈沖磁控濺射技術(shù)放電不穩(wěn)定、濺射材料離化率差異大等缺陷.然而其產(chǎn)生的等離子體僅能依靠濃度差擴(kuò)散的方式向基體運(yùn)動, 沉積速率并沒有明顯改善, 尤其是在遠(yuǎn)離陰極區(qū)域.使用外擴(kuò)型磁場對離子運(yùn)動進(jìn)行引導(dǎo), 可實(shí)現(xiàn)等離子體的聚焦和遠(yuǎn)距離輸運(yùn), 從而減少離子損失, 提高沉積效率.本文從模擬和實(shí)驗(yàn)的角度對磁場的布局與設(shè)置進(jìn)行研究, 并獲得不同磁場條件下的等離子體空間和時(shí)間輸運(yùn)特性及其對薄膜沉積的影響.結(jié)果表明電磁場的引入不僅可以大幅提高筒形陰極內(nèi)等離子體的引出效率, 實(shí)現(xiàn)不同程度的引出或聚焦, 而且對等離子體放電也產(chǎn)生明顯的增強(qiáng)或減弱, 可根據(jù)不同的需求或材料進(jìn)行精確調(diào)控.通過控制磁場, 可獲得較強(qiáng)的HiPIMS放電和較高的沉積速率, 實(shí)驗(yàn)結(jié)果與仿真預(yù)測相符合.該工作完善了HiPIMS沉積技術(shù)在沉積效率上的不足, 拓寬了筒形陰極的濺射工藝窗口和適用范圍, 有助于HiPIMS更進(jìn)一步的推廣與應(yīng)用.

    1 引 言

    高功率脈沖磁控濺射 (high power impulse magnetron sputtering, HiPIMS)技術(shù)[1]是使用高脈沖電壓實(shí)現(xiàn)高強(qiáng)度放電的新型PVD技術(shù), 其較高濺射材料離化率[2,3]可增加沉積粒子的能量可控性, 進(jìn)而實(shí)現(xiàn)薄膜結(jié)構(gòu)的精確調(diào)控[4].但其瞬時(shí)的高強(qiáng)度脈沖放電模式使得其放電穩(wěn)定性和重復(fù)性變差[4], 不同的濺射材料表現(xiàn)出明顯的離化率差異[5],同時(shí)離子在高靶電壓作用下的回吸也使得沉積速率大幅下降[6], 這一系列的不足導(dǎo)致其工業(yè)應(yīng)用發(fā)展較慢[7].

    目前, 有關(guān)HiPIMS的研究主要是圍繞上述問題進(jìn)行的, 為了改善其放電穩(wěn)定性, 包括直流、射頻、中頻在內(nèi)的多種HiPIMS復(fù)合技術(shù)[8?11]被提出, 通過低功率的預(yù)離化基本消除了脈沖電流相對于脈沖電壓的延遲現(xiàn)象, 同時(shí)減緩了“打火”頻率,提高了該技術(shù)的放電穩(wěn)定性.但是低功率的復(fù)合技術(shù)在等離子體束流中引入了太多未離化濺射原子,使其離化率嚴(yán)重下降[12].?apek等[13]通過改變靶背面磁場強(qiáng)度來實(shí)現(xiàn)HiPIMS的穩(wěn)定放電, 但是磁場強(qiáng)度大小和分布卻十分依賴放電材料, 即不同的材料需要的磁場不同, 無法滿足陰極的通用性.多家電源公司或研究單位, 如 Huettinger, Hauzer,Melec和哈工大等, 嘗試通過改進(jìn)電源技術(shù), 減小“打火”反饋時(shí)間或通過電流歸一化來實(shí)現(xiàn)放電穩(wěn)定性和重復(fù)性, 事實(shí)上, 反饋再快, “打火”仍在發(fā)生, 電流的歸一化也并不代表放電的完全一致[14,15].兩年前我們提出了筒形陰極設(shè)計(jì)理念[16,17], 一方面利用類空心陰極效應(yīng)[18]增強(qiáng)等離子體碰撞, 進(jìn)而減緩離化對放電峰值功率的需求, 另一方面通過筒形結(jié)構(gòu)緩解放電“打火”和材料差異對涂層產(chǎn)生的影響, 從而保證了不同材料的放電過程中等離子體束流離化率的一致性, 并實(shí)現(xiàn)了高質(zhì)量薄膜的生長[19].然而, 在使用筒形陰極進(jìn)行HiPIMS放電時(shí),其等離子體向基片上的輸運(yùn)只能依靠濃度差來驅(qū)動, 且溢出筒形陰極后呈發(fā)散狀分布, 離子運(yùn)動方向的不一致性導(dǎo)致?lián)p失較大, 薄膜沉積速率并沒有得到明顯的改善.

    除了通過改變放電增強(qiáng)離子束流以外, 提高沉積速率的主要方法就是通過電場或磁場控制等離子體的輸運(yùn)過程, 減少離子損失, 使產(chǎn)生的離子盡可能多的沉積為薄膜.其中電場方法主要是對基體施加偏壓吸引離子向基體運(yùn)動并沉積[20,21]和增加輔助陽極改變電子運(yùn)動軌跡, 并由電子牽引離子向基體運(yùn)動并沉積[17,22].但由于HiPIMS放電強(qiáng)度大, 等離子體密度高, 電場鞘層較薄, 對離子的作用有限, 且較高的偏壓也會引起反向?yàn)R射, 因此對沉積速率的提高效果并不明顯[23,24].磁場可以有效地約束電子運(yùn)動路徑, 進(jìn)而牽引離子束流, 形成聚焦或擴(kuò)散引出[25,26].比如擴(kuò)散型的非平衡磁控濺射就是通過磁感線約束電子運(yùn)動, 從而牽引離子增加薄膜沉積速率[27,28]; Anders和 Brown[29]以 及Li等[30]也曾通過添加直流電磁線圈在陰極前形成徑向磁場, 顯著增強(qiáng)了離子運(yùn)動的方向性, 使得沉積速率提高.

    鑒于此, 本文針對筒形陰極, 提出用電磁場約束離子束流的引出方案, 由于筒形陰極的放電平面與離子引出方向平行, 外加電磁場不僅可以約束輸出離子的運(yùn)動方向, 而且能夠調(diào)控靶面的磁場, 對放電進(jìn)行調(diào)節(jié).結(jié)果表明合適的線圈結(jié)構(gòu)和磁場強(qiáng)度可以有效地增強(qiáng)HiPIMS放電, 且能明顯提高薄膜沉積效率.

    2 實(shí) 驗(yàn)

    2.1 電磁場的實(shí)現(xiàn)

    圖1 筒形陰極及電磁線圈結(jié)構(gòu)示意圖Fig.1.Schematic diagram of the cylindrical cathode and the electromagnetic coils.

    電磁線圈由直徑為5 mm的銅導(dǎo)線纏繞而成,通過控制線圈扎數(shù)分布實(shí)現(xiàn)對磁場強(qiáng)度分布的調(diào)節(jié).電磁線圈安裝在筒形陰極兩側(cè), 為了增加離子束流的沉積范圍, 將電磁線圈設(shè)計(jì)成外擴(kuò)型, 即不同位置的纏繞直徑不同, 具體尺寸如圖1所示.在距離筒形陰極較近的位置用較厚的電磁線圈和較小的纏繞直徑, 而在遠(yuǎn)端線圈的纏繞直徑較大, 不同位置的磁場強(qiáng)度隨到筒形陰極的距離增加而逐漸減弱, 其磁感線可與靶面背后永磁鐵的磁感線相互作用, 形成對等離子體的連續(xù)作用.電磁線圈形成的軸向(z方向)磁場, 一方面束縛進(jìn)入線圈的電子, 使其沿磁感線方向做螺旋運(yùn)動, 進(jìn)而引導(dǎo)離子在等離子體自洽電勢的作用下同時(shí)沿磁感線向基體運(yùn)動; 另一方面作用于陰極靶面, 當(dāng)電磁場與永磁體磁場方向一致時(shí), 靶面橫向磁場增強(qiáng), 當(dāng)電磁場與永磁體磁場方向相反時(shí), 靶面橫向磁場減弱,且陰極與線圈之間的縱向磁場增強(qiáng), 進(jìn)而對等離子體放電產(chǎn)生不同影響.

    2.2 模擬方法

    針對外擴(kuò)型的電磁線圈的仿真研究可分為兩部分, 即電磁場仿真和等離子體輸運(yùn)過程仿真.根據(jù)筒形陰極內(nèi)外不同的物理過程, 對筒形陰極內(nèi)的輝光和擴(kuò)散區(qū)域和電磁線圈內(nèi)的引出區(qū)域分別進(jìn)行仿真.由于筒形陰極具有旋轉(zhuǎn)對稱性, 可將其簡化為一個(gè)旋轉(zhuǎn)截面, 并設(shè)置柱坐標(biāo)系進(jìn)行計(jì)算.圖1中z軸(軸向)為中央旋轉(zhuǎn)軸,r軸(徑向)表示與z軸垂直的徑向方向, 原點(diǎn)選取在筒形陰極的邊緣.

    2.2.1 磁場模擬

    在進(jìn)行磁場仿真時(shí), 所做簡化如下:選取尺寸為半徑40 cm, 高80 cm的圓柱形空間作為電磁場的仿真區(qū)域, 固定其中心與筒形陰極的中心重合.假設(shè)筒形陰極的永磁鐵和電磁鐵磁場只在仿真區(qū)域內(nèi)起作用, 因此, 仿真區(qū)域的邊界條件為磁絕緣邊界條件.仿真區(qū)域內(nèi)的總磁場B相當(dāng)于永磁鐵的磁場Bmagnet和電磁鐵磁場Bcoil的矢量和, 如(1)式所示, 其中永磁鐵磁場Bmagnet和電磁線圈磁場Bcoil可分別通過方程組(2)和(3)進(jìn)行求解.

    (2)式中Hmagnet(A/m)為永磁鐵的磁場強(qiáng)度,Vmagnet(A)為永磁鐵的磁標(biāo)勢,μ0= 4π×10–7N/A2為真空磁導(dǎo)率,μr為相對磁導(dǎo)率.(3)式中Hcoil(A/m)為電磁線圈的磁場強(qiáng)度,J(A/m2)為電流密度,Ncoil,Acoil(m2)和Icoil(A)分別為電磁線圈的匝數(shù), 截面積和電流,e= 1.6×1019C 為單位電荷.在中央旋轉(zhuǎn)軸邊界上, 磁場在徑向的上的導(dǎo)數(shù)為 0, 即

    在仿真區(qū)域其他邊界上的磁絕緣邊界條件可描述為磁場在邊界法向上的分量為0, 如(5)式所示, 式中n表示邊界上的法向單位向量.

    2.2.2 筒形陰極內(nèi)等離子體運(yùn)動仿真

    利用粒子網(wǎng)格/蒙特卡羅法 (particle in cell/Monte Carlo collision, PIC/MCC)[31,32]對筒形陰極內(nèi)部進(jìn)行模擬, 分別計(jì)算不同電磁線圈電流(15,10, 5, 0, –5, –10 和–15 A)條件下的電磁場狀態(tài)和等離子體在磁場作用下的運(yùn)動和分布情況.模型采用Matlab軟件進(jìn)行編輯, 為簡化計(jì)算量, 模型只考慮簡單的Ar氣放電[33], 包括4種主要組分:電子e, 氬原子Ar, 氬離子Ar+和激發(fā)態(tài)氬原子Arm,各組分之間主要的反應(yīng)如表1所列.筒形陰極兩側(cè)的屏蔽罩為陽極并且接地, 靶面為陰極, 施加靶電壓為–800 V, 放電氣壓為 1 Pa, 初始等離子體密度設(shè)置為 1×1013m–3.仿真開始時(shí), 在模型中分別隨機(jī)取3600個(gè)仿真離子和仿真電子, 每個(gè)粒子代表2×106個(gè)實(shí)際粒子.由于仿真區(qū)域?yàn)橐?guī)則的矩形,初始時(shí)刻將仿真區(qū)域劃分成多個(gè)尺寸為1 mm ×1 mm的正方形網(wǎng)格.隨著仿真的進(jìn)行, 由于等離子體的德拜長度逐漸減小, 網(wǎng)格尺寸也隨之減小并始終小于德拜長度.為保證模型計(jì)算的準(zhǔn)確性, 需要在仿真過程中考慮仿真粒子的三維速度.模型的計(jì)算時(shí)間步長為1×1012s, 每個(gè)步長后都利用泊松方程計(jì)算一次等離子體的自洽電勢.

    表1 Ar氣放電的主要反應(yīng)表Table 1.The main reactions of simple Ar gas discharge.

    2.2.3 筒形陰極外等離子體輸運(yùn)仿真

    本文以 Cr+離子為例, 采用 Comsol Multiphysics軟件對等離子體中Cr+離子束流在磁場中的輸運(yùn)特性進(jìn)行仿真, 計(jì)算在不同電磁線圈電流的條件下, 筒形陰極輝光區(qū)域外Cr+離子的密度分布, 并換算為沉積速率, 研究其隨空間的變化.由于離開離化區(qū)域后等離子體密度會顯著下降, 因此可認(rèn)為此時(shí)等離子體的德拜長度較大, 將仿真區(qū)域劃分成尺寸為 1 mm × 1 mm 的正方形網(wǎng)格.同時(shí)可認(rèn)為仿真區(qū)域內(nèi)等離子體整體呈電中性, 且Cr+離子只與背景氣體Ar(1 Pa)發(fā)生彈性碰撞, 即忽略擴(kuò)散過程中的電離和電荷交換.因此Cr+離子密度n(m–3)可由連續(xù)性方程((6)式)描述.

    其中D(m2/s)為Cr+離子的擴(kuò)散系數(shù), 可分為平行于磁感線的分量D//和垂直于磁感線的分量D⊥,可由(7)式求出[34].

    其中e= 1.6×1019C 為單位電荷,m= 8.68×10–26kg 為 Cr+離子質(zhì)量,n= 2×105Hz 為 Cr+離子與背景 Ar原子的碰撞頻率,Te= 3 V 為電子溫度.wt= 251 為霍爾系數(shù), 因此D//?D⊥, 表示離子傾向沿著磁感線方向運(yùn)動.借鑒之前工作[19],將等離子體離化區(qū)域的邊界設(shè)置為源邊界, 形狀為半圓形, 滿足第一類邊界條件, Cr+離子密度恒定為 1×1020m–3.在離化區(qū)域以外, 認(rèn)為等離子體呈電中性.其他邊界被設(shè)置為自由邊界, 離子以玻姆速度vbohm[35]溢出, 其通量Г(1/m2s)為:

    2.3 放電測試

    將配有電磁線圈的筒形陰極安裝在自主設(shè)計(jì)的矩形真空腔 (600 mm × 600 mm × 500 mm)內(nèi), 使用哈爾濱工業(yè)大學(xué)研制的HiPIMS電源進(jìn)行放電[15], 其輸出電壓為 0—1000 V, 脈寬為 50—400 μs, 頻率為 50—300 Hz.放電實(shí)驗(yàn)使用 99.99%的高純氬氣作為工作氣體, 濺射靶為高純鉻(99.9%).放電時(shí)氣壓設(shè)置為 1 Pa, 靶電壓脈沖頻率為 100 Hz, 脈沖寬度為 300 μs.使用等離子體發(fā)射 光 譜 儀 (optical emission spectroscopy, OES)測試筒形陰極中央的等離子體強(qiáng)度.為了得到沉積速率隨空間的變化, 將 16 cm × 16 cm 的基片放置在距離筒形陰極不同的位置進(jìn)行Cr薄膜沉積,并選取基片上不同的點(diǎn)進(jìn)行厚度測量.

    3 結(jié)果與討論

    圖2所示為不同的電磁線圈電流對筒形陰極內(nèi)外磁場強(qiáng)度分布的仿真結(jié)果, 其中圖2(a)為電磁線圈電流為0 A時(shí)的情況, 等效于不加線圈的筒形陰極磁場分布, 此時(shí)靶面的磁感線形狀呈飽滿的半圓形, 如右側(cè)放大圖所示, 此時(shí)磁感線形狀與HiPIMS放電時(shí)的輝光形狀符合[19].當(dāng)在電磁線圈內(nèi)施加電流后, 電磁線圈內(nèi)部形成外擴(kuò)型磁場分布, 其外擴(kuò)尺寸與電磁線圈設(shè)計(jì)尺寸近似.當(dāng)線圈電流沿順時(shí)針方向流動時(shí), 電磁線圈與筒形陰極相鄰位置處, 兩者的磁極方向相同, 形成的磁感線相互排斥, 此時(shí)靶面的磁感線形狀的軸向(z方向)尺寸逐漸增大, 徑向(r方向)尺寸逐漸減小, 說明此時(shí)靶面的軸向磁場Bz增大, 徑向磁場Br減小.隨著電流的增大, 磁感線幾乎貫穿整個(gè)仿真區(qū)域, 并在其內(nèi)部形成了幾乎平行的分布, 如圖2(b)—圖2(d)所示, 本文定義此時(shí)的電磁線圈電流方向?yàn)檎?反之, 當(dāng)筒形陰極永磁鐵磁極和相鄰電磁鐵磁極方向相反時(shí), 二者磁場相互吸引, 此時(shí)從電磁鐵磁極發(fā)散出的磁感線絕大多數(shù)進(jìn)入永磁鐵, 靶面磁感線形狀的軸向尺寸逐漸減小, 徑向尺寸逐漸增大, 說明靶面的軸向磁場Bz減小, 徑向磁場Br增大.隨電磁線圈電流的增大, 仿真區(qū)域內(nèi)的磁感線被分為兩部分, 如圖2(e)—圖2(g) 所示, 本文定義此時(shí)的電磁線圈電流方向?yàn)樨?fù)向.

    圖2 磁感線分布與電磁線圈電流的關(guān)系, 右側(cè)為筒形陰極內(nèi)的磁感線分布放大圖 (a) 0 A; (b) 5 A; (c) 10 A; (d) 15 A; (e)–5 A; (f) –10 A; (g) –15 A.圖中紫色虛線為筒形陰極中央軸, 藍(lán)色方框中的部分為筒形陰極的放大區(qū)域Fig.2.The distribution of the magnetic induction line vs.the current in electromagnetic coils:(a) 0 A; (b) 5 A; (c) 10 A; (d) 15 A;(e) –5 A; (f) –10 A and (g) –15 A.The purple dotted line is the center axis of the cylindrical cathode and the magnified area of the cylindrical cathode is marked by the blue square.

    筒形陰極靶面的軸向磁場Bz和徑向磁場Br分別對靶面的放電刻蝕和電子溢出的速度有重要影響[36,37], 因此, 施加不同的電磁線圈電流, 其等離子體放電狀態(tài)和演化趨勢也將產(chǎn)生明顯變化.以圖2中的磁場分布作為背景場, 用PIC/MCC模型對筒形陰極內(nèi)Ar等離子體的演變趨勢進(jìn)行仿真 (800 V, 1 Pa), 得到 Ar+離子密度分布如圖3所示, 為減少計(jì)算量, 數(shù)據(jù)取自等離子體演變基本穩(wěn)定時(shí), 即計(jì)算 0.5 μs后.仿真結(jié)果如圖3 所示,仿真區(qū)域中的Ar+離子密度分布與對應(yīng)的靶面附近磁感線的形狀基本符合, 說明磁感線對等離子體分布產(chǎn)生了明顯作用.以等離子體密度高于1×1017m–3的區(qū)域作為輝光區(qū)域, 其邊界即離子密度等于 1×1017m–3的等值線, 如圖3 中藍(lán)線所示.當(dāng)電磁線圈電流為正時(shí), 如圖3(b)—圖3(d) 所示,隨電流的增加, 主要輝光區(qū)域的徑向尺寸減小, 軸向尺寸增大, 說明靶面的放電寬度逐漸增大, 但電子從筒形陰極兩側(cè)邊緣的逃逸更加容易.隨著等離子體放電的主要離化區(qū)域逐漸被磁感線壓縮, 高密度等離子體遠(yuǎn)離中央軸, 導(dǎo)致電子在筒形陰極內(nèi)部的徑向震蕩效果減弱, 即筒形陰極結(jié)構(gòu)引入的空心陰極效應(yīng)的作用效果隨之減弱.反之, 當(dāng)電磁線圈電流為負(fù)時(shí), 隨線圈電流的增加, 如圖3(e)—圖3(g)所示, 主要輝光區(qū)域的徑向尺寸增大, 軸向尺寸減小, 說明靶面放電寬度減小, 但電子從筒形陰極邊緣的逃逸減弱.隨著等離子體離化區(qū)域被徑向拉伸, 并逐漸靠近中央軸, 電子在筒形陰極內(nèi)部往復(fù)運(yùn)動的效果增強(qiáng), 離子在中央軸附近形成聚焦, 密度大幅上升, 表明空心陰極效應(yīng)的作用隨負(fù)向電流的增大而逐漸增強(qiáng).

    對圖3中PIC/MCC模型的計(jì)算結(jié)果進(jìn)行統(tǒng)計(jì), 得到靶面放電寬度和離化區(qū)域內(nèi)的最大等離子密度如表2所列.當(dāng)電磁線圈電流為正時(shí), 隨著電流的增大, 靶面的放電寬度逐漸增大, 靶面刻蝕更加均勻, 靶材利用率會增加.在電磁線圈電流為15 A的條件下, 放電寬度最大達(dá)到4.2 cm, 靶面刻蝕區(qū)占整個(gè)靶面積的比例超過80%.與靶面的放電寬度不同, 等離子體放電強(qiáng)度呈現(xiàn)出隨線圈電流的增加而逐漸減弱的趨勢, 這主要是因?yàn)殡娮釉谕残侮帢O兩側(cè)邊緣處的軸向逃逸嚴(yán)重, 等離子體得不到快速累積, 故離化區(qū)域內(nèi)等離子體密度較低.反之, 當(dāng)電磁線圈電流為負(fù)時(shí), 隨著電流的增大, 靶面的放電寬度減小, 靶面刻蝕區(qū)面積減少, 在電磁線圈電流為–15 A的條件時(shí), 放電寬度減小到只有2.3 cm,靶面刻蝕區(qū)占整個(gè)靶面積的比例不足50%.但電子逃逸被限制, 等離子體密度大大增加, 空心陰極效應(yīng)的作用也增強(qiáng), 在電磁線圈電流為–5 A 時(shí), 離化區(qū)域內(nèi)的等離子體密度達(dá)到最高.但隨著電磁線圈電流繼續(xù)增大, 靶面軸向磁場強(qiáng)度減小, 使電子的徑向逃逸加快[37], 等離子體密度轉(zhuǎn)而減小.

    圖3 不同電磁線圈電流條件下筒形陰極內(nèi) Ar+離子分布 (a) 0 A; (b) 5 A; (c) 10 A; (d) 15 A; (e) –5 A; (f) –10 A; (g)–15 A.圖中紫色虛線為中央旋轉(zhuǎn)軸, 灰色區(qū)域?yàn)榘胁? 藍(lán)線為 1×1017 m–3 的密度等值線Fig.3.The Ar+ ion distribution in the cylindrical cathode vs.the current in electromagnetic coils:(a) 0 A; (b) 5 A; (c) 10 A; (d) 15 A; (e) –5 A; (f) –10 A and (g) –15 A.The purple dotted line is the center axis of the cylindrical cathode, the gray square is the target and the blue line is the contour line of 1×1017 m–3.

    為驗(yàn)證磁場對筒形陰極放電閾值的影響, 測量了1.5 Pa條件下筒形陰極的HiPIMS起輝電壓,其隨電磁線圈電流的變化如圖4所示.在電磁線圈電流由–15 A 逐漸減小到 0 A, 再增加到 15 A 的過程中, 起輝電壓先減小后上升, 說明放電難度先下降后上升, 與表2中等離子體密度變化趨勢相符合.在電磁線圈電流為–5 A 的條件下, 起輝電壓最小, 為 391 V, 此時(shí), 徑向磁場的增加減小了電子的軸向逃逸, 靶面軸向磁場均勻且足夠大, 磁場強(qiáng)度保持在45 mT左右, 可有效阻止電子的徑向逃逸,使得盡可能多的電子在靶面積累[36], 等離子體密度迅速增加.可見, 通過改變電磁線圈電流的大小和方向, 可對等離子體放電產(chǎn)生明顯的增強(qiáng)或減弱作用, 改變等離子體放電區(qū)域的具體位置和幾何形狀, 調(diào)節(jié)空心陰極效應(yīng)的作用效果, 實(shí)現(xiàn)對筒內(nèi)等離子體放電狀態(tài)的調(diào)控.

    表2 靶面放電寬度和最大的等離子密度Table 2.Target discharge width and maximum plasma density.

    圖4 HiPIMS放電的起輝電壓隨電磁線圈電流的變化Fig.4.Threshold voltage of HiPIMS discharge vs.the current in the electromagnetic coil.

    圖5為不同電磁線圈電流條件下整個(gè)模擬空間的等離子體密度擴(kuò)散分布圖, 用來表征外加磁場對等離子體擴(kuò)散路徑的影響.可以看出, 等離子體在模擬區(qū)域的運(yùn)動仍然以擴(kuò)散為主, 隨著擴(kuò)散距離的增加, 等離子體密度逐漸降低, 但是擴(kuò)散過程中受到磁場作用, 其擴(kuò)散路徑和聚集位置發(fā)生了明顯變化.與電磁線圈電流為零時(shí)(即無外加磁場)相比(圖5(a)), 電磁線圈提供的外加磁場減緩了等離子體的擴(kuò)散損耗速度, 使得筒形陰極內(nèi)的等離子體密度大大增加, 且效果隨電磁線圈電流的增大而提高, 這是由于磁場使得電子在向外擴(kuò)散過程中產(chǎn)生螺旋運(yùn)動, 延長了運(yùn)動路徑, 從而產(chǎn)生累積, 使等離子體密度增加[26].在電磁線圈電流為正時(shí), 由于筒形陰極兩端的徑向磁場被大大降低, 而靶面的軸向磁場卻增加, 使得電子更容易從筒形陰極兩端溢出, 而較難繼續(xù)向筒形陰極中央軸運(yùn)動.在電子的帶動下, 等離子體在放電區(qū)產(chǎn)生后即可一邊向外輸運(yùn), 一邊向中央軸擴(kuò)散, 在筒形陰極側(cè)面形成相對均勻的等離子體密度分布, 如圖5(b)—圖5(d)所示.隨線圈電流的增加, 該趨勢越來越明顯, 在筒形陰極內(nèi)徑范圍內(nèi)均勻輸出的距離也不斷增加.當(dāng)電磁線圈電流為負(fù)時(shí), 其等離子體密度分布結(jié)果與線圈電流為正時(shí)相反, 其靶面的軸向磁場減弱, 而陰極兩端的徑向磁場增加, 使得電子從筒形陰極邊緣向外擴(kuò)散的難度增大, 而向中央軸擴(kuò)散更容易.故在電子的帶動下, 等離子體呈現(xiàn)出先向中央軸擴(kuò)散, 然后再向外擴(kuò)散的趨勢, 在筒形陰極外形成一個(gè)高度聚焦的等離子體束流, 中心濃度高, 四周濃度低, 如圖5(e)—圖5(g)所示.

    磁場能改變等離子體的沉積路徑, 因此對等離子體沉積速率及其分布也會產(chǎn)生重要影響.根據(jù)等離子密度數(shù)據(jù)計(jì)算特定位置的離子通量, 再除以計(jì)算出的離子通量最大值進(jìn)行歸一化處理, 即可獲得相應(yīng)位置上的歸一化沉積速率, 如圖6所示.圖6(a)為筒形陰極中央軸上距離陰極側(cè)面出口不同處的沉積速率, 可見, 隨擴(kuò)散距離的增加, 其沉積速率逐漸減小, 兩種線圈電流模式都增加了沉積速率,并隨電磁線圈電流的增大而增大, 說明電磁線圈的引入可以明顯優(yōu)化等離子體的沉積效率.主要原因是磁感線的束縛作用限制了運(yùn)動方向較偏而沉積到真空室壁上的等離子體損失, 使大部分離子束流沿磁感線運(yùn)動并形成了有效沉積.電磁線圈電流為正時(shí), 等離子體向外運(yùn)動的阻力較小, 而向中央軸運(yùn)動的阻力較大, 故即使靠近陰極側(cè)面出口位置處 (0—3 cm 時(shí)), 其沉積速率提升并不明顯, 甚至在線圈電流較小時(shí)還出現(xiàn)少量的降低.但隨擴(kuò)散距離的增加, 中央軸上的離子濃度得到不斷補(bǔ)充, 使得中央軸上的沉積速率隨擴(kuò)散距離的增加衰減逐漸減弱, 因此其沉積速率相對不加線圈電流時(shí)有了明顯提高, 尤其是距出口4 cm以上時(shí)沉積速率的改善效果更加明顯.但在電磁線圈電流為負(fù)時(shí), 由于陰極兩端的徑向磁場大, 而靶面的軸向磁場小,使得電子更傾向于向中央軸聚集, 然后向外擴(kuò)散,于是中央軸上的沉積速率遠(yuǎn)大于無線圈電流或線圈電流為正時(shí), 尤其是距離陰極側(cè)面出口越近, 其沉積速率越高.

    圖5 等離子體密度分布與電磁線圈電流的關(guān)系 (a) 0 A; (b) 5 A; (c) 10 A; (d) 15 A; (e) –5 A; (f) –10 A; (g) –15 A.圖中紫色虛線為中央旋轉(zhuǎn)軸Fig.5.The plasma density distribution vs.the current in electromagnetic coils:(a) 0 A; (b) 5 A; (c) 10 A; (d) 15 A; (e) –5 A; (f)–10 A and (g) –15 A.The purple dotted line is the center axis of the cylindrical cathode.

    圖6(b)所示為在不同電磁線圈電流下, 距離筒形陰極側(cè)面出口8 cm處沉積速率的徑向分布,相對無線圈電流的情況, 有線圈電流時(shí)的沉積速率均有增加.當(dāng)線圈電流為正時(shí), 隨線圈電流增加,沉積速率增加, 其徑向分布曲線逐漸出現(xiàn)一個(gè)平臺, 寬度隨電磁線圈電流的增大而增大, 說明此時(shí)在陰極側(cè)面出口范圍內(nèi)離子沉積速率的均勻性較好, 尤其是線圈電流 15 A 時(shí), 可在–4—4 cm 的范圍內(nèi)保持較高速度的均勻沉積.而當(dāng)線圈電流為負(fù)時(shí), 雖然沉積速率增加更明顯, 可實(shí)現(xiàn)快速沉積,但其沉積速率曲線在輸出范圍內(nèi)成峰狀, 中央軸上的沉積速率遠(yuǎn)大于兩側(cè), 意味著輸出范圍內(nèi)沉積均勻性較差.但正是由于等離子體向中央軸的聚集,使得其因漂移到真空室壁而損失的離子數(shù)量大大減少, 可實(shí)現(xiàn)更高的離子利用率.

    為進(jìn)一步研究距離陰極出口不同位置處的沉積速率的空間分布, 計(jì)算了線圈電流為15 A時(shí)不同位置的沉積速率徑向分布, 如圖6(c)所示.在距離較小 (3 cm)的截面上, 沉積速率分布不均勻, 曲線呈現(xiàn)兩個(gè)峰, 對應(yīng)于離子束流被引出的位置, 說明此時(shí)等離子體尚未在中央軸位置形成最大聚集.隨著擴(kuò)散距離的增大, 中央軸處的等離子體得到補(bǔ)充, 曲線中央的凹陷逐漸消失形成平臺, 在6 cm處形成最均勻的沉積.隨著截面距離的進(jìn)一步增大, 兩側(cè)等離子體的消耗較快, 形成中心沉積速率大于兩側(cè)的現(xiàn)象.

    圖6 沉積速率的空間分布 (a) 不同線圈電流下中央軸沉積速率隨擴(kuò)散距離的變化; (b) 不同線圈電流下距筒形陰極8 cm處截面上的沉積速率分布; (c) 線圈電流15 A時(shí)不同距離的截面上的沉積速率分布Fig.6.Spatial distribution of the deposition rates:(a) Deposition rate on the central axis depends on the diffusion distance and the coil currents; (b) deposition rate at the distance of 8 cm from the cylindrical cathode at different coil currents; (c) deposition rate distribution at different distances from the cylindrical cathode when the coil current is 15 A.

    圖7 不同線圈電流時(shí)的等離子體束流輝光照片 (a) 線圈電流為正; (b) 線圈電流為負(fù)Fig.7.The plasma flow pictures vs.the electromagnetic coil currents:(a) The coil current is positive; (b) the coil current is negative.

    為驗(yàn)證仿真結(jié)果, 在不同的線圈電流下進(jìn)行HiPIMS 放電測試, 放電總功率為 3 kW, 其放電輝光照片如圖7所示.可見, 照片中的等離子體聚集情況基本與仿真結(jié)果(圖5)一致, 即筒形陰極外等離子體輝光強(qiáng)度隨線圈電流增加而增強(qiáng).當(dāng)線圈電流為正時(shí), 筒形陰極出口處的輝光在出口尺寸范圍內(nèi)分布均勻, 隨著線圈電流的增加, 輝光范圍變大,輝光變亮, 說明等離子體束流密度增大, 且沿磁感線均勻輸出, 如圖7(a)所示.當(dāng)線圈電流為負(fù)時(shí),輝光區(qū)域呈現(xiàn)收縮狀態(tài), 集中從中央位置往外擴(kuò)展.隨著線圈電流的增加, 輝光區(qū)域的長度增加,輝光更亮, 意味著離子束流的聚焦程度增大, 引出效率大幅提高, 如圖7(b)所示.

    圖8 不同的線圈電流時(shí)筒形陰極中央軸上的光譜強(qiáng)度Fig.8.The OES intensity at the center axis of cylindrical cathode vs.the electromagnetic coil currents.

    為了進(jìn)一步驗(yàn)證線圈電流對空間的等離子體密度分布的作用, 采用等離子體發(fā)射光譜儀探測了筒形陰極中央軸上的發(fā)射光譜隨線圈電流的變化規(guī)律, 如圖8 所示.結(jié)果表明:線圈電流為負(fù)時(shí), 隨線圈電流的增大, 光譜強(qiáng)度大幅增強(qiáng), 說明筒形陰極中央的等離子體密度上升明顯, 這主要來源于等離子體在中央軸位置的聚集以及空心陰極效應(yīng)的增強(qiáng).在線圈電流為–15 A 時(shí), 光譜強(qiáng)度達(dá)到了無線圈電流時(shí)的4倍以上.而線圈電流為正時(shí), 隨線圈電流的增加, 光譜強(qiáng)度逐漸減弱, 說明此時(shí)等離子體向中央軸的擴(kuò)散受到抑制, 中央軸位置處的等離子體密度較低, 空心陰極效應(yīng)因此減弱.

    根據(jù)我們之前的工作[19], 無外加磁場情況下筒形陰極只能在直徑為1.5 cm的范圍內(nèi)保持均勻沉積, 且在距離陰極出口15 cm處沉積速率不足40 nm/min.配備電磁線圈后, 以純金屬 Cr靶進(jìn)行HiPIMS放電為例, 研究了線圈電流為15和–15 A時(shí)的Cr薄膜沉積速率, 結(jié)果分別如圖9(a)和圖9(b)所示.可見, 在線圈電流為 15 A 時(shí), 不同距離上的薄膜沉積速率較均勻, 且大體上隨著與筒形陰極端口距離的增加呈現(xiàn)下降的趨勢, 與仿真結(jié)果一致.但當(dāng)距離陰極出口較近 (3 cm)時(shí), 實(shí)驗(yàn)與仿真結(jié)果差異較大, 此時(shí)基片上的沉積速率相對仿真結(jié)果較低且均勻性較差.這是由于實(shí)驗(yàn)中距離陰極較近的基片上會形成懸浮電位, 從而影響等離子體的擴(kuò)散, 使離子向基片的運(yùn)動受阻, 從而減少了沉積速率, 而仿真中忽略了基片懸浮電位對擴(kuò)散的影響.隨著基片與陰極出口距離的增大, 整個(gè)基片上薄膜沉積均勻性得到明顯改善, 其演化趨勢與仿真結(jié)果(圖5(c))相符合.與無線圈電流相比, 線圈電流為15 A時(shí)薄膜的沉積速率略有提高, 但薄膜的均勻沉積范圍顯著增大.圖9(b)是線圈電流為–15 A 時(shí), 基片在距離筒形陰極出口 15 cm 處薄膜徑向沉積結(jié)果, 可見, 基片中央部分的沉積速率遠(yuǎn)高于兩側(cè), 與仿真結(jié)果相符合.基片中央位置的沉積速率相對無線圈電流時(shí)提高4倍以上, 而兩側(cè)位置沉積速率與無線圈電流時(shí)類似, 只是略微提高.

    圖9 (a)線圈電流為 15 A 時(shí)沉積速率的分布; (b)線圈電流為–15 A 時(shí)沉積速率的分布Fig.9.Distribution of the deposition rate when the electromagnetic coil current is (a) 15 A and (b) –15 A.

    4 結(jié) 論

    為了改善并控制筒形陰極在HiPIMS放電時(shí)的等離子體狀態(tài)和沉積效率, 本文提出電磁場控制和牽引方案, 通過研究線圈電流大小和方向?qū)ν残侮帢OHiPIMS放電等離子體的作用, 實(shí)現(xiàn)了放電區(qū)域內(nèi)等離子體密度可控性和等離子體沉積效率的大幅度提高.結(jié)果表明:增加電磁場后, 等離子體密度都獲得不同程度的改變, 通過磁感線對電子的束縛作用, 可實(shí)現(xiàn)等離子體在不同位置的聚集和不同的擴(kuò)散路徑, 從而實(shí)現(xiàn)薄膜沉積速率和沉積均勻性的有效調(diào)控.當(dāng)線圈電流為正時(shí), 電子的軸向運(yùn)動阻力小, 徑向運(yùn)動阻力大, 空心陰極效應(yīng)減弱,等離子體束流傾向于均勻輸出.當(dāng)線圈電流為負(fù)時(shí), 電子的軸向運(yùn)動阻力增大, 徑向運(yùn)動阻力減小,電子更傾向于在中央軸位置聚集, 然后向外擴(kuò)散,形成沿中央軸的高速輸出.該電磁場設(shè)計(jì)策略不僅實(shí)現(xiàn)了等離子體的增強(qiáng)和有效控制, 而且完善了HiPIMS沉積技術(shù)在沉積效率上的明顯不足, 結(jié)合筒形陰極的放電穩(wěn)定性和高純度離子輸出特性, 基本為HiPIMS技術(shù)的工業(yè)應(yīng)用奠定了基礎(chǔ).

    猜你喜歡
    靶面陰極等離子體
    Evaluation of Arctic Sea Ice Drift and its Relationship with Near-surface Wind and Ocean Current in Nine CMIP6 Models from China
    連續(xù)磁活動對等離子體層演化的影響
    基于低溫等離子體修飾的PET/PVC浮選分離
    船舶高壓水射流除銹過程有限元仿真分析
    廣東造船(2020年1期)2020-05-03 14:01:49
    傾斜螺紋孔射流傳熱特性的數(shù)值研究
    場發(fā)射ZrO/W肖特基式場發(fā)射陰極研究進(jìn)展
    電子制作(2018年12期)2018-08-01 00:47:46
    實(shí)驗(yàn)研究類螺紋孔旋流沖擊射流的冷卻特性
    等離子體種子處理技術(shù)介紹
    多彈丸時(shí)空散布參數(shù)測試方法仿真研究
    IT-SOFCs陰極材料Sm0.8La0.2Ba1-xSrxFe2O5+δ的制備與表征
    99精品欧美一区二区三区四区| 我的亚洲天堂| 最近最新免费中文字幕在线| 国产亚洲精品久久久久5区| 日韩欧美免费精品| 久久精品成人免费网站| 18禁黄网站禁片午夜丰满| 亚洲av电影不卡..在线观看| 嫁个100分男人电影在线观看| 大型av网站在线播放| 长腿黑丝高跟| 搡老岳熟女国产| 欧美精品亚洲一区二区| 国产1区2区3区精品| 国产精品久久电影中文字幕| 99精品久久久久人妻精品| videosex国产| 精华霜和精华液先用哪个| 午夜免费鲁丝| 操出白浆在线播放| 日韩三级视频一区二区三区| 亚洲国产看品久久| 国产亚洲精品第一综合不卡| 岛国在线观看网站| 精品高清国产在线一区| 伦理电影免费视频| 国产久久久一区二区三区| av片东京热男人的天堂| 香蕉国产在线看| 在线观看午夜福利视频| 亚洲一码二码三码区别大吗| 亚洲狠狠婷婷综合久久图片| 日韩欧美一区视频在线观看| 妹子高潮喷水视频| 99国产精品一区二区三区| 制服人妻中文乱码| 国产精品亚洲av一区麻豆| 亚洲性夜色夜夜综合| 久久国产精品男人的天堂亚洲| 成年版毛片免费区| 亚洲欧美日韩高清在线视频| 夜夜夜夜夜久久久久| 欧美日韩福利视频一区二区| 中亚洲国语对白在线视频| 精品国产乱码久久久久久男人| 久久精品影院6| 久热这里只有精品99| 日韩精品青青久久久久久| 日本黄色视频三级网站网址| 一级黄色大片毛片| 一二三四在线观看免费中文在| 中文字幕精品免费在线观看视频| 欧美最黄视频在线播放免费| 在线播放国产精品三级| 在线观看免费午夜福利视频| 亚洲七黄色美女视频| 免费av毛片视频| 久久人人精品亚洲av| 精品少妇一区二区三区视频日本电影| 日日摸夜夜添夜夜添小说| 中文字幕精品亚洲无线码一区 | 97人妻精品一区二区三区麻豆 | 欧美一区二区精品小视频在线| 岛国在线观看网站| 9191精品国产免费久久| 国产激情欧美一区二区| 中文在线观看免费www的网站 | 欧美一级a爱片免费观看看 | 香蕉国产在线看| 丝袜人妻中文字幕| 国产精品久久久久久人妻精品电影| 91av网站免费观看| 久久久久久久午夜电影| 亚洲av美国av| 男女做爰动态图高潮gif福利片| 大型黄色视频在线免费观看| 国产精品98久久久久久宅男小说| 91在线观看av| 欧美中文综合在线视频| 99热6这里只有精品| 免费人成视频x8x8入口观看| 男女床上黄色一级片免费看| 国产高清视频在线播放一区| 老汉色∧v一级毛片| 欧美性长视频在线观看| 夜夜躁狠狠躁天天躁| 欧洲精品卡2卡3卡4卡5卡区| 久久性视频一级片| 亚洲五月天丁香| 制服丝袜大香蕉在线| a级毛片a级免费在线| 99re在线观看精品视频| 欧美+亚洲+日韩+国产| 国产99白浆流出| 国产视频一区二区在线看| 国产精品1区2区在线观看.| 亚洲va日本ⅴa欧美va伊人久久| 搡老熟女国产l中国老女人| 亚洲午夜精品一区,二区,三区| 在线观看66精品国产| 色哟哟哟哟哟哟| 亚洲五月色婷婷综合| 欧美日韩瑟瑟在线播放| 757午夜福利合集在线观看| 亚洲精品一卡2卡三卡4卡5卡| 熟女电影av网| 成人特级黄色片久久久久久久| 国语自产精品视频在线第100页| 国产精品一区二区精品视频观看| 麻豆成人av在线观看| 曰老女人黄片| av电影中文网址| netflix在线观看网站| 午夜久久久久精精品| 久久精品国产99精品国产亚洲性色| 黄色片一级片一级黄色片| 一级片免费观看大全| 日韩欧美免费精品| 丝袜美腿诱惑在线| 中文字幕人妻丝袜一区二区| 国产精品久久久人人做人人爽| 国产三级在线视频| 久久久国产欧美日韩av| 色在线成人网| 亚洲 国产 在线| 99久久久亚洲精品蜜臀av| 天天躁夜夜躁狠狠躁躁| 男人舔女人的私密视频| 97碰自拍视频| 99国产极品粉嫩在线观看| 自线自在国产av| 久久国产精品男人的天堂亚洲| 久久国产精品男人的天堂亚洲| 国产主播在线观看一区二区| 免费av毛片视频| 日本精品一区二区三区蜜桃| 成人免费观看视频高清| 俄罗斯特黄特色一大片| 国产亚洲精品久久久久5区| 免费搜索国产男女视频| 大型av网站在线播放| 精品久久久久久久末码| 成人免费观看视频高清| 亚洲无线在线观看| 国产91精品成人一区二区三区| 日韩欧美一区二区三区在线观看| 亚洲成人精品中文字幕电影| 国产男靠女视频免费网站| 亚洲av成人不卡在线观看播放网| 男人舔女人的私密视频| 国产爱豆传媒在线观看 | 一本综合久久免费| 757午夜福利合集在线观看| 美女国产高潮福利片在线看| 国语自产精品视频在线第100页| 国产精品久久电影中文字幕| 最新美女视频免费是黄的| 欧美一级a爱片免费观看看 | 久久精品国产综合久久久| 国产乱人伦免费视频| 久久久久久免费高清国产稀缺| 欧美亚洲日本最大视频资源| 亚洲真实伦在线观看| 日韩精品免费视频一区二区三区| 国产精品,欧美在线| 亚洲熟妇熟女久久| 亚洲精品美女久久av网站| 国产一区二区激情短视频| 久久久水蜜桃国产精品网| 日韩三级视频一区二区三区| 久久99热这里只有精品18| 久久婷婷人人爽人人干人人爱| 男女下面进入的视频免费午夜 | 亚洲精品久久成人aⅴ小说| 久热这里只有精品99| 亚洲成av片中文字幕在线观看| av有码第一页| 亚洲av第一区精品v没综合| a在线观看视频网站| 国产极品粉嫩免费观看在线| 免费看日本二区| 窝窝影院91人妻| 亚洲国产日韩欧美精品在线观看 | 一级a爱视频在线免费观看| 色哟哟哟哟哟哟| 久久人人精品亚洲av| 国产亚洲精品第一综合不卡| 草草在线视频免费看| 国产欧美日韩精品亚洲av| 精品久久久久久久毛片微露脸| 成人午夜高清在线视频 | 精品一区二区三区av网在线观看| 一级毛片女人18水好多| 韩国av一区二区三区四区| 一级黄色大片毛片| 久久人人精品亚洲av| 美女高潮喷水抽搐中文字幕| 一个人免费在线观看的高清视频| 欧美日韩乱码在线| 成人午夜高清在线视频 | 欧美大码av| 日韩一卡2卡3卡4卡2021年| 久久香蕉激情| 国产麻豆成人av免费视频| 少妇的丰满在线观看| 久久久久久久久免费视频了| 超碰成人久久| 999精品在线视频| 久久欧美精品欧美久久欧美| 午夜免费成人在线视频| 国产主播在线观看一区二区| av有码第一页| 国产成人精品无人区| 国产熟女午夜一区二区三区| 亚洲欧美精品综合久久99| 欧美一级a爱片免费观看看 | 国内精品久久久久精免费| 亚洲av中文字字幕乱码综合 | av有码第一页| 丝袜人妻中文字幕| 久久精品国产99精品国产亚洲性色| 亚洲欧洲精品一区二区精品久久久| 免费在线观看日本一区| 又黄又爽又免费观看的视频| 欧美+亚洲+日韩+国产| 宅男免费午夜| 精品久久久久久久毛片微露脸| 69av精品久久久久久| 免费看十八禁软件| 欧美色视频一区免费| 两人在一起打扑克的视频| 亚洲欧美精品综合久久99| 精品日产1卡2卡| 观看免费一级毛片| 波多野结衣巨乳人妻| 欧美日韩精品网址| 中文亚洲av片在线观看爽| 天天躁夜夜躁狠狠躁躁| 中文字幕人成人乱码亚洲影| 99热只有精品国产| 国产精品九九99| 欧美性猛交╳xxx乱大交人| 这个男人来自地球电影免费观看| 无人区码免费观看不卡| 国产人伦9x9x在线观看| 亚洲专区字幕在线| 国产精品综合久久久久久久免费| 婷婷六月久久综合丁香| 嫩草影院精品99| 国产色视频综合| 日韩欧美国产在线观看| 亚洲va日本ⅴa欧美va伊人久久| 麻豆久久精品国产亚洲av| 一区二区日韩欧美中文字幕| 国内毛片毛片毛片毛片毛片| 99国产精品一区二区三区| 一级a爱视频在线免费观看| 久久欧美精品欧美久久欧美| 欧美黑人精品巨大| 一个人免费在线观看的高清视频| 最近在线观看免费完整版| 国产成人啪精品午夜网站| 在线观看舔阴道视频| 国产精品综合久久久久久久免费| 国产不卡一卡二| av有码第一页| 日韩大码丰满熟妇| 久9热在线精品视频| www.精华液| svipshipincom国产片| 久久久久久久久中文| 不卡av一区二区三区| 亚洲av成人不卡在线观看播放网| 黄片小视频在线播放| 亚洲精品美女久久久久99蜜臀| 日本a在线网址| 桃红色精品国产亚洲av| 制服诱惑二区| 中文字幕精品亚洲无线码一区 | 91大片在线观看| 美女国产高潮福利片在线看| 欧美日本亚洲视频在线播放| 日韩三级视频一区二区三区| 99久久无色码亚洲精品果冻| 夜夜夜夜夜久久久久| 黄色毛片三级朝国网站| 日韩欧美 国产精品| 亚洲精品久久成人aⅴ小说| 国产精品自产拍在线观看55亚洲| 久久久国产成人免费| 黄色a级毛片大全视频| 久久香蕉精品热| 中文字幕人成人乱码亚洲影| 午夜福利免费观看在线| 亚洲色图 男人天堂 中文字幕| 久久午夜综合久久蜜桃| 国产黄片美女视频| 男人舔奶头视频| 久久精品国产亚洲av香蕉五月| 久久香蕉国产精品| 欧美乱妇无乱码| 国产亚洲欧美精品永久| 黑人欧美特级aaaaaa片| 国产精品国产高清国产av| 国产又爽黄色视频| 成人午夜高清在线视频 | 少妇熟女aⅴ在线视频| 嫩草影视91久久| 一区二区日韩欧美中文字幕| 丝袜美腿诱惑在线| 欧美色视频一区免费| 国产激情久久老熟女| 国产精品亚洲美女久久久| 亚洲自拍偷在线| 亚洲精华国产精华精| 亚洲在线自拍视频| 国产亚洲精品久久久久久毛片| 久久精品aⅴ一区二区三区四区| 视频区欧美日本亚洲| 一边摸一边抽搐一进一小说| 精品久久久久久久久久久久久 | av欧美777| 制服诱惑二区| 老司机深夜福利视频在线观看| 色综合亚洲欧美另类图片| 日韩大码丰满熟妇| 国产精品野战在线观看| 亚洲真实伦在线观看| 国产伦人伦偷精品视频| 国产区一区二久久| 日日摸夜夜添夜夜添小说| 亚洲一区二区三区色噜噜| 亚洲精品国产精品久久久不卡| 成年人黄色毛片网站| 午夜成年电影在线免费观看| 不卡av一区二区三区| 嫩草影视91久久| АⅤ资源中文在线天堂| 国内精品久久久久久久电影| 久久久久免费精品人妻一区二区 | 久热这里只有精品99| 18禁美女被吸乳视频| 午夜老司机福利片| 亚洲欧美日韩高清在线视频| 国产亚洲av高清不卡| 久久精品夜夜夜夜夜久久蜜豆 | 在线观看www视频免费| 亚洲精品中文字幕在线视频| 日韩欧美一区二区三区在线观看| 国产亚洲av嫩草精品影院| 亚洲欧美日韩高清在线视频| 国产一级毛片七仙女欲春2 | xxx96com| 国产精品久久久久久亚洲av鲁大| cao死你这个sao货| 色播在线永久视频| 午夜福利一区二区在线看| 国产极品粉嫩免费观看在线| 麻豆成人av在线观看| 国内毛片毛片毛片毛片毛片| 免费观看精品视频网站| 男女下面进入的视频免费午夜 | 人人妻人人看人人澡| 黄片小视频在线播放| 两人在一起打扑克的视频| 欧美日韩黄片免| 精品免费久久久久久久清纯| 国产精品一区二区免费欧美| 无遮挡黄片免费观看| 99在线视频只有这里精品首页| 国产亚洲精品综合一区在线观看 | 亚洲成人久久性| 亚洲狠狠婷婷综合久久图片| 成人三级做爰电影| a级毛片在线看网站| 波多野结衣巨乳人妻| 草草在线视频免费看| 久久性视频一级片| 9191精品国产免费久久| 一a级毛片在线观看| 久久九九热精品免费| 中文字幕精品免费在线观看视频| 久久亚洲精品不卡| 最近在线观看免费完整版| 天天一区二区日本电影三级| 美女高潮喷水抽搐中文字幕| 精品无人区乱码1区二区| 老熟妇乱子伦视频在线观看| 黄色视频不卡| 熟女少妇亚洲综合色aaa.| 日本成人三级电影网站| 天堂√8在线中文| 国产成人欧美在线观看| 国产精品亚洲一级av第二区| 黄片小视频在线播放| 免费在线观看黄色视频的| netflix在线观看网站| 精品国产一区二区三区四区第35| 一区二区三区精品91| 精品久久久久久久人妻蜜臀av| 人人妻人人澡欧美一区二区| 国产亚洲精品久久久久5区| 中文字幕久久专区| 欧美日本亚洲视频在线播放| 色精品久久人妻99蜜桃| 欧美一级a爱片免费观看看 | 999久久久精品免费观看国产| 90打野战视频偷拍视频| 亚洲成a人片在线一区二区| 日韩精品中文字幕看吧| 男女床上黄色一级片免费看| 国产视频内射| 亚洲最大成人中文| av在线播放免费不卡| 久久九九热精品免费| 亚洲中文日韩欧美视频| 97碰自拍视频| 欧美人与性动交α欧美精品济南到| 精品高清国产在线一区| 99re在线观看精品视频| av超薄肉色丝袜交足视频| 国产成人av激情在线播放| 欧美又色又爽又黄视频| 最近最新免费中文字幕在线| 一级毛片精品| 巨乳人妻的诱惑在线观看| 麻豆国产av国片精品| 岛国在线观看网站| 成人免费观看视频高清| 日本一本二区三区精品| 国产不卡一卡二| 99热6这里只有精品| 精品乱码久久久久久99久播| 嫩草影院精品99| 午夜福利成人在线免费观看| 精品久久久久久,| 在线观看免费午夜福利视频| 精品不卡国产一区二区三区| 日韩欧美免费精品| 亚洲在线自拍视频| 精品国产一区二区三区四区第35| 国产日本99.免费观看| 亚洲精品国产区一区二| 欧美在线一区亚洲| 日日夜夜操网爽| 国产一区在线观看成人免费| 18禁黄网站禁片午夜丰满| 午夜激情av网站| 成人免费观看视频高清| 久久午夜综合久久蜜桃| 亚洲国产看品久久| 一本一本综合久久| 久久精品国产综合久久久| 国产av一区二区精品久久| 精品无人区乱码1区二区| 一本综合久久免费| 久久久久久久久久黄片| 十分钟在线观看高清视频www| 女同久久另类99精品国产91| 黄色视频不卡| 99久久综合精品五月天人人| 99久久久亚洲精品蜜臀av| 大香蕉久久成人网| 黑人巨大精品欧美一区二区mp4| 男女之事视频高清在线观看| 久久精品国产综合久久久| 女生性感内裤真人,穿戴方法视频| 人人澡人人妻人| 国产成人系列免费观看| 两人在一起打扑克的视频| 亚洲中文字幕日韩| 18禁美女被吸乳视频| 我的亚洲天堂| 99re在线观看精品视频| 欧美人与性动交α欧美精品济南到| 老司机深夜福利视频在线观看| 国产成年人精品一区二区| 国产97色在线日韩免费| 黄网站色视频无遮挡免费观看| 99国产综合亚洲精品| 中文在线观看免费www的网站 | 欧美一级a爱片免费观看看 | 久久婷婷成人综合色麻豆| 欧美日韩乱码在线| 精品国产乱子伦一区二区三区| 精品国产一区二区三区四区第35| 后天国语完整版免费观看| 日韩国内少妇激情av| 日韩三级视频一区二区三区| 成熟少妇高潮喷水视频| 大型黄色视频在线免费观看| av超薄肉色丝袜交足视频| 在线国产一区二区在线| 男女做爰动态图高潮gif福利片| 99久久国产精品久久久| 欧美成人午夜精品| 精品久久久久久久毛片微露脸| 国产人伦9x9x在线观看| 亚洲成人精品中文字幕电影| 后天国语完整版免费观看| 美女国产高潮福利片在线看| 国产一区二区激情短视频| 色尼玛亚洲综合影院| 国产野战对白在线观看| 亚洲成av人片免费观看| 女性被躁到高潮视频| 久久久久久大精品| 欧美大码av| 18禁黄网站禁片免费观看直播| 99精品久久久久人妻精品| 欧美精品亚洲一区二区| 久久中文字幕一级| 中文字幕久久专区| 亚洲av成人av| 中出人妻视频一区二区| 嫩草影视91久久| 久久性视频一级片| 手机成人av网站| 久久中文字幕一级| 一级毛片女人18水好多| av在线播放免费不卡| 欧美色欧美亚洲另类二区| 国产一级毛片七仙女欲春2 | 真人做人爱边吃奶动态| 国产精品,欧美在线| 亚洲av电影不卡..在线观看| 亚洲精品在线观看二区| 精品少妇一区二区三区视频日本电影| 久久中文字幕一级| 国产精品野战在线观看| 午夜免费观看网址| 国产单亲对白刺激| 91字幕亚洲| 成年人黄色毛片网站| 一二三四在线观看免费中文在| 免费观看人在逋| 麻豆一二三区av精品| 亚洲成人国产一区在线观看| 久久久久精品国产欧美久久久| 亚洲精品一区av在线观看| 一级作爱视频免费观看| 女同久久另类99精品国产91| 最近最新中文字幕大全免费视频| 久久精品国产综合久久久| 99在线人妻在线中文字幕| 伦理电影免费视频| 日本 av在线| 日本三级黄在线观看| 色在线成人网| 夜夜夜夜夜久久久久| 欧美在线黄色| 亚洲欧美精品综合久久99| 午夜福利在线在线| 最近在线观看免费完整版| 夜夜看夜夜爽夜夜摸| 久久久久国内视频| 日本成人三级电影网站| 日本 欧美在线| 国语自产精品视频在线第100页| 久久精品国产亚洲av高清一级| 91在线观看av| 国产黄色小视频在线观看| 在线十欧美十亚洲十日本专区| 一夜夜www| 在线十欧美十亚洲十日本专区| 午夜亚洲福利在线播放| 欧美日韩精品网址| 国产精品 国内视频| 国产精品一区二区精品视频观看| av视频在线观看入口| 国产精品免费视频内射| 亚洲精品久久国产高清桃花| 色老头精品视频在线观看| av视频在线观看入口| 亚洲五月色婷婷综合| a在线观看视频网站| 国产99白浆流出| 国产麻豆成人av免费视频| 自线自在国产av| 欧美中文日本在线观看视频| 日韩欧美一区视频在线观看| 夜夜躁狠狠躁天天躁| 色综合站精品国产| 国产精品乱码一区二三区的特点| 99久久无色码亚洲精品果冻| 日本黄色视频三级网站网址| 中文字幕人妻熟女乱码| 男女午夜视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 91麻豆精品激情在线观看国产| 欧美最黄视频在线播放免费| 叶爱在线成人免费视频播放| www日本在线高清视频| 中亚洲国语对白在线视频| av免费在线观看网站| 亚洲五月天丁香| 999久久久精品免费观看国产| 久久久精品国产亚洲av高清涩受| 波多野结衣av一区二区av| 亚洲男人的天堂狠狠| 久久久久久久久久黄片| 国产激情偷乱视频一区二区| 51午夜福利影视在线观看| 亚洲欧美精品综合一区二区三区| 亚洲成av人片免费观看| 大型黄色视频在线免费观看| 久久久久国产一级毛片高清牌| 久久99热这里只有精品18| 欧美日韩亚洲国产一区二区在线观看| 欧美日本亚洲视频在线播放| 亚洲精品国产一区二区精华液| 免费av毛片视频| 亚洲欧美精品综合久久99| 日韩成人在线观看一区二区三区| 精品日产1卡2卡|