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

    用于飛秒激光納米加工的TiO2粒子?陣列誘導多種基底表面近場增強

    2017-08-08 06:57:50焦悅陶海巖季博宇宋曉偉林景全
    物理學報 2017年14期
    關鍵詞:近場飛秒折射率

    焦悅 陶海巖季博宇 宋曉偉林景全

    (長春理工大學理學院,長春 130022)

    用于飛秒激光納米加工的TiO2粒子?陣列誘導多種基底表面近場增強

    焦悅 陶海巖?季博宇 宋曉偉?林景全

    (長春理工大學理學院,長春 130022)

    (2017年1月20日收到;2017年3月11日收到修改稿)

    利用納米粒子輔助對飛秒激光能量進行空間局域化,使其在基底表面誘導產(chǎn)生納米尺度的近場增強,這對超衍射極限微結構加工具有重要意義.目前對于粒子陣列誘導飛秒激光納米孔加工的研究僅限于金屬Au粒子及低折射率聚苯乙烯介電粒子等,本文提出并開展了應用高折射率TiO2介電粒子陣列作為輔助誘導激光近場增強從而進行飛秒激光超衍射納米孔加工的研究.對TiO2介電粒子陣列在Si,Pt及SiO2表面的近場強度分布進行了數(shù)值模擬,研究其基底表面近場增強的規(guī)律及物理過程.研究結果發(fā)現(xiàn),使用硅基底時,陣列與單一TiO2球形粒子相比其近場增強僅下降約30%;相對于入射激光強度而言,在直徑約為100 nm的空間范圍內(nèi)獲得140倍的近場增強,這一現(xiàn)象可用于百納米孔的激光加工.同時在其他典型基底的理論計算結果中也表明,幾乎在所有金屬及介電材料表面均可以實現(xiàn)良好的百納米空間范圍內(nèi)的近場增強,并且具有近場隨著基底折射率變大而增強的規(guī)律.這些現(xiàn)象的產(chǎn)生歸因于TiO2粒子中磁四極振蕩產(chǎn)生的激光前向場增強及粒子與基底的耦合作用.進一步引入鏡像電荷模型對基底光學參數(shù)對其表面近場增強的影響規(guī)律進行了分析和解釋.本文的模擬結果對飛秒激光近場超衍射極限納米加工的應用有著重要的意義.

    飛秒激光,超衍射極限加工,近場光學,TiO2粒子

    1 引 言

    器件微小型化一直是醫(yī)療、生物、信息等領域的發(fā)展需求,也是先進制造發(fā)展的重要方向.短脈沖激光同時具有超強和超快的獨特優(yōu)勢,特別是飛秒激光,非熱的激光燒蝕機制使其具有加工精度高及低損傷等諸多優(yōu)點.因此飛秒激光微納制備技術在多個領域都有著廣泛的應用[1-3].但由于光學衍射極限的制約(≈λ/2,λ為入射激光波長),使超衍射極限尺度下納米結構的激光加工應用受到了限制.近年來,研究人員提出了一種基于納米結構輔助的超衍射極限激光加工方案:利用納米結構可對激光近場進行空間局域化的特點,在待加工基底表面產(chǎn)生超衍射極限空間尺度的激光近場增強,從而借助此區(qū)域內(nèi)增強激光場的燒蝕效應形成納米孔.例如:利用激光照射顯微系統(tǒng)的金屬探針尖端(如原子力顯微鏡[4,5]等)與基底的間隙,形成局域近場增強,這一方式可實現(xiàn)百納米特征尺寸的納米加工.但這種方法速度非常緩慢,實現(xiàn)大面積陣列加工困難.在待加工的硅基底表面放置Au納米粒子也能取得相似的燒蝕效果[6].研究發(fā)現(xiàn),這種通過放置納米粒子的方式更有利于大面積、高效率的加工,因此納米粒子輔助激光近場納米加工得到了廣泛的關注.到目前為止,國外已經(jīng)開展了對單一粒子近場增強及應用于納米燒蝕的研究[7],實現(xiàn)了不同形狀Au納米結構(Au納米棒[8]和Au納米球等[9])對硅表面的輔助加工.另外,還研究了Au納米球?qū)Σ煌纵o助加工效果的影響[10].實驗結果表明,Au粒子可有效應用于半導體(Si)[10]和金屬(Pt)等[11]基底表面的納米孔加工.這是由于基底表面可提供自由電子參與等離激元共振,使基底與粒子間局域化的近場強度得到大幅度提升[12].但是由于等離激元共振的產(chǎn)生對粒子所處介電環(huán)境要求很高(對基底材料選擇局限)[13,10],因此,Au納米粒子在應用于其他低折射率基底上(如SiO2,聚甲基丙烯酸甲酯(PMMA)等)時,表面等離激元幾乎不能被激發(fā),其近場強度要較Si基底表面的情況低幾個量級,實驗上也證實無法得到清晰的納米孔加工效果[14].

    對于單一粒子研究的主要目的是對其機理的研究,而粒子輔助激光納米加工的提出,是基于其陣列應用時的快速制備優(yōu)勢,因此陣列情況下的研究十分必要.然而研究發(fā)現(xiàn),粒子陣列緊密排列時與單一放置情況下相比,其在基底表面近場增強效果并不相同,而且相差甚遠.例如Au粒子,即使Au粒子在單一放置情況下,在近場增強效果最為明顯的Pt基底上[11]緊密排布時,基底表面的近場增強效果不但不能保持而且還會變得微弱,這是由于粒子間強的場耦合作用,使得基底表面加工難以實現(xiàn).目前報道的有效陣列納米孔制備方式是聚苯乙烯(PS)粒子陣列排布在PMMA基底表面, PS粒子的直徑接近入射激光波長時可以產(chǎn)生粒子透鏡效應,對光進行“聚焦”[15,16],但是這種方式對基底介電環(huán)境要求苛刻,PS粒子與其他基底組合后聚焦位置會脫離靶材,進而失去加工能力.因此,如硅等具有廣泛應用價值的基底如何實現(xiàn)粒子陣列輔助激光納米加工成為亟待解決的難題.最近人們發(fā)現(xiàn),單一的小尺寸高折射率TiO2介電粒子在低折射率基底SiO2上通過飛秒激光照射可實現(xiàn)直徑為百納米孔的加工[17].然而到目前為止,尚未開展對TiO2粒子在不同基底表面近場增強的研究.同時,由于金屬粒子與高折射率介電粒子對激光吸收及近場增強機理存在差異[18],TiO2粒子陣列誘導納米孔加工效果有望區(qū)別于Au粒子陣列的情況,有必要開展其加工能力的研究.

    本文利用時域有限差分法(FDTD)數(shù)值計算方法對Si,SiO2及Pt基底表面TiO2納米粒子的近場分布進行數(shù)值模擬,研究了當TiO2球形粒子陣列排布時,在基底表面實現(xiàn)超衍射極限近場納米加工的能力.模擬結果發(fā)現(xiàn),在硅基底表面,納米球形粒子陣列與單一情況下獲得的近場增強相比,并沒有出現(xiàn)Au粒子陣列那樣因粒子間相互作用使基底表面處近場發(fā)生大幅度下降的情況.所獲得的增強因子使加工區(qū)域與其周圍的環(huán)境具有很好的強度差,這對于納米加工是十分理想的.TiO2納米粒子在SiO2及Pt等典型基底表面均可以實現(xiàn)納米加工的近場增強.進一步對基底光學參數(shù)對其表面近場增強的影響規(guī)律進行了分析和解釋.研究結果表明,原理上TiO2粒子可以應用于各種金屬及介電基底上超衍射極限納米孔的飛秒激光加工.這一工作對納米結構輔助的激光超衍射極限納米加工的發(fā)展具有重要意義.

    2 數(shù)值模擬方法

    本文通過基于FDTD的數(shù)值計算軟件FDTD solutions(Lumerical,Ltd)對基底表面單一及陣列TiO2納米球形粒子周圍的近場分布進行了模擬計算[19].因在實際應用中,通過一定的技術工藝(通常為旋涂法)可獲得基底表面單層球緊湊排布的陣列[15],因此在本文研究過程中使用的粒子陣列均為這種排列方式,如圖1(a)所示.在模擬研究中,將單一或陣列排列的TiO2納米粒子分別放置在Si,SiO2及Pt三種不同材料的基底表面.粒子直徑為D,計算中使用參考文獻[20,21]中給出的Si, SiO2及Pt的介電常數(shù)值,TiO2介電常數(shù)的實部和虛部分別取為2.66和0.024[22].入射激光場波長選擇了400 nm,這樣選擇的原因主要是考慮了這一波長在TiO2磁四極共振譜線內(nèi),同時也是近紅外飛秒激光輸出經(jīng)倍頻后的工作波長.電場強度設置為1 V/m,為了確?;妆砻娴慕鼒龇植紴閷ΨQ分布,所有模擬中均使用圓偏振光.為了更好地從近場強度推測加工能力,在這里引用增強因子的概念,將基底表面有無粒子情況下同一位置光強的比值定義為增強因子,即

    其中E為有粒子情況下基底表面在入射光照射后獲得的電場值,Esub為同一位置無粒子情況下的基底表面在入射光照射后獲得的電場值(區(qū)別于入射光場電場值E0).通過增強因子值,可以從理論上獲得近場增強點和基底之間的激光場強度差,從而在加工過程中更好地控制激光能量密度,使其在確保不損傷基底無粒子區(qū)域表面的情況下,在有粒子區(qū)域進行納米加工.

    圖1 (網(wǎng)刊彩色)納米粒子陣列(黑色圓)六角排列在基底表面(藍色區(qū)域)的示意圖 模擬計算時選擇周期性邊界條件(紅色線)和完美吸收邊界條件(紫色線),入射光為波長400 nm的圓偏振光,逆著z方向垂直入射,電場強度為1 V/mFig.1.(color online)Arrayed particles(black circles)hexagonally close-packed on the substrate(blue area)is shown in the(a)and(b).The boundaries with periodic conditions and perfectly matched layers(PMLs)are shown with red and purple lines,respectively.The total field/scattered field technique was employed to simulate the incident field. The incident electromagnetic wave is circularly polarized in x-y plane,and it propagates along the-z direction. The wavelength is 400 nm.The incident electric field intensity is 1 V/m.

    如圖1所示,在模擬計算中,入射光逆著z方向入射.將基底表面所在的平面設為x-y平面,基底表面上方為z軸正方向,粒子與基底接觸點在x-z平面內(nèi).x和y方向的邊界條件為周期性邊界(period boundary),z方向為完美吸收邊(PML boundary),監(jiān)視器位于粒子中心和基底接觸點所在的x-z平面和基底表面x-y平面.

    3 結果與分析

    圖2(a)為單一TiO2納米粒子在Si基底表面近場增強因子隨粒子直徑D的變化曲線,其中D的取值范圍是100—800 nm.從圖2(a)中可以看出,隨粒子直徑D的增大,增強因子Q先快速增大,后又快速減弱,之后趨于穩(wěn)定.在D=300 nm時,Q值最大.圖2(b)為基底表面的俯視近場分布圖,可以發(fā)現(xiàn)基底表面近場增強主要區(qū)域直徑小于100 nm,其場聚集的空間尺度已遠突破入射激光的衍射極限.

    為了對比粒子陣列排布與單一放置時對基底表面近場分布的影響,進一步對TiO2納米粒子陣列的增強因子及近場分布進行了模擬研究,如圖3所示.圖3(a)為Si基底表面TiO2納米粒子陣列的近場增強因子Q隨粒子直徑D的變化曲線,其中D值與單一粒子的情況取值相同,可以看出,列陣情況與單一粒子的情況增強因子的變化趨勢相似,也是粒子直徑為300 nm時,增強因子值最大.圖3(b)為粒子直徑為300 nm時基底表面的近場分布圖,從圖中可以觀察到,當納米粒子陣列緊密排列在基底表面時,可以獲得規(guī)則排列的近場分布,且基底表面光場增強區(qū)域由直徑均小于100 nm的陣列組成,與單一粒子時近場增強情況相似.但是粒子陣列相對單一粒子情況,基底表面近場增強強度下降約30%.根據(jù)米氏散射理論可知,當激光照射納米球形粒子時,納米球形粒子會對其吸收截面放大,增加能量的吸收[14].當球形粒子緊密排列時,陣列中的單個球體的散射截面會相對減小[11],因此對入射激光的吸收也會減小,球形粒子整體近場總強度也會隨之減小,這樣球與基底間電場增強幅度也會隨之減小.

    圖2 (網(wǎng)刊彩色)(a)Si基底表面單一TiO2納米粒子的光強增強因子隨粒子直徑D的變化曲線;(b)D=300 nm(圖(a)中增強因子最大值對應的粒子尺寸)時基底表面的近場分布圖,其中虛線代表粒子的外圍輪廓,入射光為波長400 nm的圓偏振光,逆著z方向垂直入射,入射電場強度為1 V/mFig.2.(color online)A single TiO2particle with di ff erent diameter deposited on Si substrate surface:(a)Simulated optical intensity enhancement factor as a function of diameter of TiO2particle;(b)electric field intensity distribution on the substrate surface when a single particle with diameter D=300 nm is placed(large peak value of enhancement factor is obtained at diameter of 300 nm),the white dotted line represents the particle delineation.The incident electromagnetic wave is circularly polarized in x-y plane,and it propagates along the-z direction;the wavelength is 400 nm,the incident electric field intensity is 1 V/m.

    圖3 (網(wǎng)刊彩色)(a)TiO2納米粒子陣列排列在Si基底表面的光強增強因子隨粒子直徑D的變化曲線;(b)D=300 nm (圖(a)中增強因子最大值對應的粒子尺寸)時基底表面的近場分布圖,其中虛線代表球的外圍輪廓,入射光為波長400 nm的圓偏振光,逆著z方向垂直入射,入射電場強度為1 V/mFig.3.(color online)Arrayed TiO2particles with di ff erent diameter deposited on Si substrate surface:(a)Simulated optical intensity enhancement factor as a function of diameter of TiO2particle;(b)electric field intensity distribution on the substrate surface when arrayed particles with diameter D=300 nm are placed(large peak value of enhancement factor is obtained at diameter of 300 nm),the white dotted line represents the particle delineation. The incident electromagnetic wave is circularly polarized in x-y plane,and it propagates along the-z direction,the wavelength is 400 nm,the incident electric field intensity is 1 V/m.

    基底表面近場強度的增強幅度與球形粒子整體的電荷分布有關,為探知TiO2納米粒子內(nèi)近場的分布,分別對單一及陣列TiO2納米粒子x-z面的近場分布進行了數(shù)值模擬.圖4(a)和圖4(b)給出了TiO2單一納米粒子以及陣列排布時在硅表面及兩相鄰粒子接觸點附近的近場增強模擬圖.由圖4(c)和圖4(d)可知,球形粒子陣列排列后,球與球接觸點位置出現(xiàn)了近場增強,同時球形粒子底部與基底接觸處的近場有所減小.高折射率的TiO2是一種半導體材料,在外界光電場的作用下,粒子表面會誘導產(chǎn)生少量表面電荷.當球形粒子緊密排列時,球的接觸點處會發(fā)生類似金屬尖端的避雷針效應(lightening rod e ff ect)[23],誘導近場的增強.這一近場的產(chǎn)生與增強會對粒子內(nèi)的電荷分布產(chǎn)生影響,進而對基底表面的近場增強產(chǎn)生作用,使整個球體的電場分布發(fā)生變化,削弱了基底表面的近場增強.由于兩粒子接觸點的近場增強相對基底并不高,因此對基底近場的影響有限.對比于Au球形粒子陣列排布的情況,金屬材料富含自由電子,在激光輻照Au粒子后,由于發(fā)生電偶極振蕩,大量電荷會分布在球體表面,由于粒子間隙的存在會產(chǎn)生十分顯著的避雷針效應,使大量的電荷發(fā)生匯聚,從而產(chǎn)生極強的增強區(qū)域,此時基底表面的光強相較于單一粒子情況會被削弱很多.這會嚴重影響甚至喪失對基底進行燒蝕加工的能力,這一現(xiàn)象已被實驗證實[11].最后對比于空氣中單一粒子,如圖4(e)所示,通過在入射光照射粒子后x-z面的近場分布結果發(fā)現(xiàn),粒子內(nèi)有較為明顯的場增強現(xiàn)象,但在球底部并未形成近場增強區(qū)域.當基底材料引入后,與球底接觸處產(chǎn)生了局域近場增強效果(如圖4(a)—圖4(d),這說明基底材料的引入對納米球的近場分布產(chǎn)生了顯著的影響.

    圖4 (網(wǎng)刊彩色)Si基底表面(a)單一TiO2納米粒子(D=300 nm)和(b)TiO2納米粒子陣列的x-z截面(通過粒子中心點與基底接觸點的面)近場分布圖,(c)沿著圖(a)中基底表面z=0 nm(白色虛線)的光強分布,(d)為圖(b)中基底表面z=0 nm和粒子之間接觸點z=125 nm的光強分布;(e)無基底時單一TiO2納米粒子近場分布圖;入射光為波長400 nm的圓偏振光,逆著z方向垂直入射,入射電場強度為1 V/mFig.4.(color online)Electric fi led intensity distribution on an x-z plane for simulation system with(a)a single TiO2particle,(b)arrayed TiO2particles with diameter of 300 nm deposited on Si substrate surface.The intensity distribution is shown on the plane through the center of the particle and the contact point between the particle and the substrate consisting of two touching gold nanoparticle placed on silicon substrate.The intensity distribution is shown in the(c)along the white dashed line in Fig.(a)(on the substrate surface,z=0 nm)and(d)along the white dashed line in Fig.(b)(through the center of the particles and the contact point between the particles and the substrate,z=125 nm and z=0 nm).(e)shows the electric field intensity distribution with a single TiO2particle in vacuum.The incident electromagnetic wave is circularly polarized in x-y plane,and it propagates along the-z direction.The wavelength is 400 nm.The incident electric field intensity is 1 V/m.

    在納米粒子輔助激光近場納米加工中,能流密度分布是決定基底燒蝕質(zhì)量的關鍵因素之一[18],為此,對在球內(nèi)及Si基底的能流分布進行了計算,結果如圖5所示.圖中箭頭方向為能流矢量的方向,箭頭的寬度和長度表示能流密度大小.從圖5(a)中可以觀察到,單一球形粒子內(nèi)的能流是從球上半部中心繞行后在下半部分進行匯聚,并流向基底表面,激光場的能流密度通過粒子后強度有明顯增強.這是由于TiO2納米粒子存在磁四極共振模式,在入射光的激發(fā)下會有前向延伸激發(fā)模式的產(chǎn)生[17],使電磁場通過粒子后產(chǎn)生強度增強.當粒子以陣列排布時,如圖5(b)所示,由于粒子與粒子接觸處場強的增加(如圖4(b)),使能流路徑相對于單一粒子更集中于球的中間區(qū)域,能流方向未發(fā)生改變,仍集中向基底流動.因此,陣列排列的粒子內(nèi)的能流的方向及大小沒有受到很大的影響,從而保證了基底表面獲得較強的近場增強.進一步發(fā)現(xiàn),陣列情況下的能流更深入到基底的內(nèi)部,這樣的能流空間分布對表面加工更加有利.而對于Au粒子的情況,其外界光場照射誘導的電偶極振蕩使電荷分布于球表面,對球內(nèi)形成了電磁屏蔽,使光電磁場(光能流)不能進入球的內(nèi)部[18],其基底表面的近場增強主要依靠等離激元振蕩.可見,TiO2和Au粒子不同的振蕩機制導致了其場增強產(chǎn)生了顯著差異.

    以上工作中基底材料為半導體Si,接下來進一步對TiO2納米粒子陣列在SiO2和Pt典型基底表面(金屬和透明介質(zhì)材料)的近場分布進行數(shù)值模擬研究.模擬結果如圖6(a)和圖6(b)所示,可以發(fā)現(xiàn)Pt和SiO2基底表面增強因子及近場分布與Si基底類似.同時對比于單一粒子近場分布情況(如圖6(c)和圖6(d)所示)也發(fā)現(xiàn),陣列粒子排布的基底表面近場強度小于單一粒子情況,這也與Si基底情況的規(guī)律相一致.基于計算的增強因子有望在這兩種材料表面實現(xiàn)超衍射納米孔加工[17].綜上,在三種不同類型的典型基底表面均可實現(xiàn)較好的近場增強效果.同時也發(fā)現(xiàn),基底材料不同時,基底對增強近場的耦合效果也有所不同.

    圖5 (網(wǎng)刊彩色)(a)單一TiO2納米粒子和(b)納米粒子陣列在Si基底表面的能流密度分布,D=300 nmFig.5.(color online)Intensity fl ow distributions on Si substrate using(a)a single TiO2particle and(b)arrayed TiO2particles,both with 300 nm diameter.

    為探索基底材料參數(shù)對TiO2納米粒子陣列近場增強的影響機制,在以上材料參數(shù)的基礎上,將基底材料折射率(N=n+ik)中的n及k分別取不同值,并對其增強因子進行計算.在金屬材料中, 400 nm對應的折射率實部n變化較小,因此在這里保持Pt基底材料的n值不變,改變其虛部k值(k取值參考常用材料的取值范圍[20]).模擬結果如圖7所示,增強因子與基底折射率虛部k成線性關系——隨著折射率虛部k數(shù)值的增加,增強因子整體呈上升趨勢.而在一般介電材料中,折射率虛部k很小或不存在,因此,保持SiO2基底材料的k值不變(k=0),改變折射率實部n值來選定介電材料(n取值范圍參考常用材料的范圍[20]).模擬結果如圖8所示,隨著基底材料折射率實部n數(shù)值的增加,近場增強因子也逐漸增加.根據(jù)圖7和圖8的模擬結果可知,原理上TiO2納米粒子陣列在不同材料基底表面均可以誘導產(chǎn)生較好的近場增強效果.

    激光輻照球形粒子時,其粒子本身的散射效率及近場分布在有無基底的情況下具有很大的差別[18].因此,在研究粒子輔助基底表面獲得增強近場分布時,必須要考慮基底的影響.基底表面近場增強的獲得主要源于粒子輔助獲得的增強電磁場與基底表面的耦合.因此,耦合過程可引入鏡像電荷模型對光強增強因子與基底折射率關系進行分析.表面形成鏡像電荷越多表示耦合得越好.因此基底材料不同時,通過粒子內(nèi)增強電場產(chǎn)生的鏡像電荷效應也不同,也就是電場耦合的程度不同.對于金屬材料,表面存在大量的自由電子,很容易被外加電場所極化,因此可近似認為吸收的激光可以使所有的表面自由電子形成鏡像電荷.而對于入射光的吸收(在趨膚深度內(nèi))取決于其折射率的虛部k,所以k也稱為吸收系數(shù)(the index of absorption),k越大對入射光的吸收越強,而n僅貢獻于相位效應(phase e ff ects),對光吸收沒有貢獻[20].

    圖6 (網(wǎng)刊彩色)TiO2納米粒子陣列在(a)Pt基底和(b)SiO2基底表面光強增強因子隨粒子直徑D的變化曲線,圖中插圖分別為Pt和SiO2基底表面粒子直徑D=300 nm時的近場分布圖;(c),(d)分別為單一情況下基底表面近場分布,其中虛線代表球的外圍輪廓;入射光為波長400 nm的圓偏振光,逆著z方向垂直入射,入射電場強度為1 V/mFig.6.(color online)Simulated optical intensity enhancement factors with di ff erent diameter arrayed particles deposited on(a)Pt substrate and(b)SiO2substrate surfaces as a function of diameter of TiO2particle,the inset shows the electric field intensity distribution on the substrate surface with diameter of 300 nm of particle;(c), (d)show the electric field intensity distribution on the Pt and SiO2substrate surface in the case of a single TiO2particle,separately,the white dotted line represents the particle delineation.The incident electromagnetic wave is circularly polarized in x-y plane,and it propagates along the-z direction.The wavelength is 400 nm.The incident electric field intensity is 1 V/m.

    通過之前的研究結果可知,基底表面獲得的近場源于激光的前向場增強,更多增強的激光近場吸收會使基底表面極化電荷數(shù)量相應增加,從而其表面近場也會隨之增強(如圖7).而對于介電材料基底,基底表面受外場作用產(chǎn)生感應電荷,此時基底內(nèi)鏡像電荷表達式為[18]

    其中,q′是鏡像電荷,q為原始電荷;εs和εa分別為基底和空氣的介電常數(shù),并且εs=n2s,ns為基底折射率.由(1)式可知,對于介電材料基底而言,當基底介電常數(shù)εs越大,也就是ns越大基底表面鏡像電荷數(shù)量越多,此時基底與近場耦合會越強.因此,增強因子Q會隨著折射率的增加而增強(如圖8).綜上,基底材料折射率的實部n值和虛部k值都會對基底表面近場增強產(chǎn)生影響,這可以歸結于粒子增強的近場與基底耦合過程——基底折射率越大(折射率實部或虛部增加)越有利于其表面耦合TiO2球形粒子輔助獲得的球內(nèi)增強近場.目前飛秒激光技術發(fā)展迅速,可以實現(xiàn)充足的激光強度供給.因此只要通過粒子陣列排布獲得相當?shù)募す鈴姸葘Ρ榷?增強因子Q)就可將其應用于激光納米超衍射加工.圖7和圖8所取基底光學參數(shù)(折射率)范圍幾乎包含了所有常見的典型固體材料,通過參考已有實驗研究結果[17],本文計算所得的近場增強模擬結果(圖3及圖6—圖8)在實驗中都有望達到很好的加工效果.因此,本文的模擬結果說明TiO2粒子陣列輔助飛秒激光近場納米加工可以適用于多種基底.

    圖7 TiO2納米粒子陣列在金屬材料基底表面的光強增強因子與基底折射率虛部(k)的函數(shù)關系(n=1.7)Fig.7.Simulated optical intensity enhancement factors with 300 nm diameter arrayed TiO2particles deposited on mental substrate surface as a function of extinction coefficient(k)of substrate with refractive index n=1.7.

    圖8 TiO2納米粒子陣列在非金屬介電材料基底表面的光強增強因子與基底折射率實部(n)的函數(shù)關系(k=0)Fig.8.Simulated optical intensity enhancement factors with 300 nm diameter arrayed TiO2particles deposited on mental substrate surface as a function of refractive index(n)of substrate with extinction coeffi cient k=0.

    4 結 論

    本文通過FDTD數(shù)值計算方法對TiO2粒子陣列輔助提高Si,Pt及SiO2基底表面局域近場進行了模擬研究.結果發(fā)現(xiàn),在Si基底時,陣列與單一TiO2相比近場增強下降小于30%,但相對于入射激光場強具有140倍的增強,可實現(xiàn)飛秒激光的超衍射加工所需的場強度差.同時,Pt和SiO2(常用于納米孔加工的基底表面)基底表面TiO2粒子陣列近場分布的計算結果表明均獲得了很好的增強效果.通過能流密度分布計算發(fā)現(xiàn),基于粒子磁四極振蕩吸收的激光電磁場,能流透射粒子過程中電場會逐漸增強,產(chǎn)生激光前向場增強效應,在陣列排布時其增強及能流方向也不受影響,這為各種基底的近場增強提供了保障.這也是TiO2粒子陣列對各種基底材料均可以產(chǎn)生近場增強的原因.進一步對基底材料的折射率取不同值,發(fā)現(xiàn)隨著折射率的增加,近場增強效果也隨之增強.這一現(xiàn)象的發(fā)生歸因于粒子內(nèi)增強的近場和基底相互作用的結果,鏡像電荷模型分析表明該現(xiàn)象與基底表面等效電荷形成有關.這一現(xiàn)象的發(fā)現(xiàn)也驗證了在其他金屬和介電材料基底表面,原理上都可以獲得良好的增強效果.本研究結果對TiO2粒子陣列輔助飛秒激光近場超衍射納米加工的應用有著重要的意義.

    [1]Tao H,Song X,Hao Z,Lin J 2015 Chin.Opt.lett.13 061402

    [2]Tao H,Lin J,Hao Z,Gao X,Song X,Sun C,Tan X 2012 Appl.Phys.Lett.100 201111

    [3]Li G Q,Li X H,Yang H D,Qiu R,Huang W H 2011 Chin.Opt.4 72(in Chinese)[李國強,李曉紅,楊宏道,邱榮,黃文浩2011中國光學4 72]

    [4]Zenhausern F,Martin Y,Wickramasinghe H K 1995 Science 269 1083

    [5]Merlein J,Kahl M,Zuschlag A,Sell A,Halm A, Boneberg J,Leiderer P,Leitenstorfer A,Bratschitsch R 2008 Nature Photon.2 230

    [6]Wang Z B,Luk’yanchuk B S,Li L,Crouse P L,Liu Z, Dearden G,Watkins K G 2007 Appl.Phys.A 89 363

    [7]Plech A,Kotaidis V,Lorenc M,Boneberg J 2006 Nature Phys.2 44

    [8]Robitaille A,Boulais é,Meunier M 2013 Opt.Express 21 9703

    [9]Nedyalkov N,Miyanishi T,Obara M 2007 Appl.Surf. Sci.253 6558

    [10]Atanasov P A,Nedyalkov N N,Sakai T,Obara M 2007 Appl.Surf.Sci.254 794

    [11]Nedyalkov N,Sakai T,Miyanishi T,Obara M 2007 Appl. Phys.Lett.90 123106

    [12]Quan S,Kosei U,Han Y,Atsushi K,Yasutaka M,Hiroaki M 2013 Light Sci.Appl.2 e118

    [13]Pe?arodríguez O,Pal U,Campoyquiles M,Rodríguezfernández L,Garriga M,Alonso M I 2011 J.Phys.Chem. C 115 6410

    [14]Terakawa M,Takeda S,Tanaka Y,Obara G,Miyanishi T,Sakai T,Sumiyoshi T,Sekita H,Hasegawa M,Viktorovitch P,Obara M 2012 Prog.Quantum Electron 36 194

    [15]Afanasiev A,Bredikhin V,Pikulin A,Ilyakov I,Shishkin B,Akhmedzhanov R,Bityurin N 2015 Appl.Phys.Lett. 106 183102

    [16]Pikulin A,Afanasiev A,Agareva N,Alexandrov A P, Bredikhin V,Bityurin N 2012 Opt.Express 20 9052

    [17]Tanaka Y,Obara G,Zenidaka A,Terakawa M,Obara M 2010 Appl.Phys.Lett.96 261103

    [18]Tanaka Y,Obara M 2009 Jpn.J.Appl.Phys.48 122002

    [19]Ta fl ove A,Hagness S C 2000 Computational Electrodynamics:The Finite-Di ff erence Time-Domain Method (Boston:Artech House)

    [20]Palik E D 1998 Handbook of Optical Constants of Solids (Vol.1)(San Diego,CA:Academic)p333

    [21]Johnson P B,Christy R W 1972 Phys.Rev.B 6 4370

    [22]Messinger B J,Raben K U,Chang R K,Barber P W 1981 Phys.Rev.B 24 649

    [23]Maier S A 2007 Plasmonics:Fundamentals and Applications(New York:Springer)

    PACS:42.62.—b,79.20.Eb,11.10.—z DOI:10.7498/aps.66.144203

    Near field enhancement of TiO2nanoparticle array on di ff erent substrates for femtosecond laser processing?

    Jiao Yue Tao Hai-Yan?Ji Bo-Yu Song Xiao-Wei?Lin Jing-Quan
    (School of Science,Changchun University of Science and Technology,Changchun 130022,China)

    20 January 2017;revised manuscript

    11 March 2017)

    Optical near field enhancement on substrate can be achieved by localizing femtosecond laser energy with nanoparticles.The enhanced field is located in the region between nanoparticles and the substrate.The localized femtosecond optical field is of great signi fi cance for fabricating the micro/nano structure with characteristic size beyond the di ff raction limit.Up to now,femtosecond processing nanohole assisted by particle array is only possible for metal particle(Au)and low-refractive-index dielectric polystyrene particle.However,previous research results show that it cannot be realized for metal particle arrays(Au)to form periodic nanohole arrays,and it is limited for polystyrene particle to choose the corresponding substrate.In this paper,a novel method is proposed,in which high refractive index TiO2arrayed particles are placed on the substrate to achieve laser induced near field enhancement.This makes feasible the nanoscale processing beyond the di ff raction limit.In this paper,near field distributions of TiO2particle array on Si,Pt and SiO2substrates are simulated by the fi nite-di ff erence time-domain(FDTD)method.The results show that TiO2particles concentrate the laser energy to a region with a diameter of 100 nm around the particle and the near field enhancement is 140 times higher than the incident laser intensity,which is bene fi cial to fabricating the nanostructure of super di ff raction limit,such as sub-hundred nanometer nanohole ablation by femtosecond laser.For Si substrate,the near field enhancement is only about 30%lower for TiO2particle array than that for single TiO2particle.In order to explore the in fl uence mechanism of the substrate material parameters on the near field enhancement of TiO2nanoparticle array,we further simulate the enhancement factor for the substrates of di ff erent refractive indices.It is found that the near field is enhanced with the increase of substrate refractive index,and this is attributed to an increased interaction of the particle with the near field of substrate and the scattering e ff ect in which the TiO2particle supports forward near field intensity pattern.Moreover, the image charge model is introduced to analyze the e ff ect of substrate optical parameters on local field enhancement. Results in this paper can be applied to most metals as well as dielectric substrate surfaces,and they open a new way for femtosecond laser near field nano-processing with characteristic size beyond the di ff raction limit.

    femtosecond laser,beyond di ff raction limit processing,near field optics,TiO2nanoparticles

    :42.62.—b,79.20.Eb,11.10.—z

    10.7498/aps.66.144203

    ?國家自然科學基金(批準號:61605017)和長春理工大學青年科學基金(批準號:XQNJJ-2015-01)資助的課題.

    ?通信作者.E-mail:hytao@cust.edu.cn

    ?通信作者.E-mail:songxiaowei@cust.edu.cn

    ?2017中國物理學會Chinese Physical Society

    http://wulixb.iphy.ac.cn

    *Project supported by the National Natural Science Foundation of China(Grant No.61605017)and the Foundation for Young Scientists of Chuangchun University of Science and Technology,China(Grant No.XQNJJ-2015-01).

    ?Corresponding author.E-mail:hytao@cust.edu.cn

    ?Corresponding author.E-mail:songxiaowei@cust.edu.cn

    猜你喜歡
    近場飛秒折射率
    超大規(guī)模智能反射面輔助的近場移動通信研究
    全飛秒與半飛秒的區(qū)別
    人人健康(2021年16期)2021-12-01 07:08:33
    基于反射型超表面的近場聚焦研究
    基于飛秒激光的固體危化品切割機床設計與開發(fā)
    溴丙烯在800nm和400nm飛秒激光強場下的解離電離
    一種基于PDV的近場沖擊波高壓測量技術
    中國測試(2018年10期)2018-11-17 01:58:50
    基于飛秒脈沖相關法的高精度時間同步測量
    近場RCS測量不確定度分析
    制導與引信(2016年3期)2016-03-20 16:02:01
    單軸晶體雙折射率的測定
    用Z-掃描技術研究量子點的非線性折射率
    物理實驗(2015年9期)2015-02-28 17:36:49
    禁无遮挡网站| 嫩草影院入口| 亚洲精品乱码久久久v下载方式| 亚洲欧美日韩卡通动漫| 亚洲国产色片| 国产高清三级在线| 色综合色国产| 亚洲自拍偷在线| 丰满人妻一区二区三区视频av| 九九热线精品视视频播放| av在线播放精品| 国产人妻一区二区三区在| 成年av动漫网址| 赤兔流量卡办理| 欧美最新免费一区二区三区| 波多野结衣高清无吗| 爱豆传媒免费全集在线观看| 亚洲五月天丁香| 韩国av在线不卡| 人妻夜夜爽99麻豆av| 国产一级毛片七仙女欲春2| 国产精品无大码| 中文精品一卡2卡3卡4更新| 偷拍熟女少妇极品色| 久久国内精品自在自线图片| 久久精品国产亚洲av天美| 免费av观看视频| 小说图片视频综合网站| 26uuu在线亚洲综合色| 91久久精品国产一区二区三区| 丝袜美腿在线中文| 性欧美人与动物交配| 国产成人a区在线观看| 伊人久久精品亚洲午夜| 国产一区二区三区av在线 | 美女 人体艺术 gogo| 亚洲欧美精品综合久久99| 精品久久久久久久末码| 久久热精品热| 久久人妻av系列| 久久6这里有精品| 99riav亚洲国产免费| 成人二区视频| 国产伦一二天堂av在线观看| 99国产极品粉嫩在线观看| 国产成年人精品一区二区| 亚洲欧美日韩卡通动漫| 老司机福利观看| a级一级毛片免费在线观看| 免费无遮挡裸体视频| 男人舔奶头视频| 亚洲自偷自拍三级| 久久久久九九精品影院| 99久久成人亚洲精品观看| 国产伦一二天堂av在线观看| 日韩精品青青久久久久久| av免费在线看不卡| 一级毛片我不卡| 亚洲欧美精品综合久久99| 看片在线看免费视频| 小说图片视频综合网站| 在线观看一区二区三区| 国产高清有码在线观看视频| 男人狂女人下面高潮的视频| 中国美白少妇内射xxxbb| 不卡一级毛片| 亚洲熟妇中文字幕五十中出| 插阴视频在线观看视频| 欧美极品一区二区三区四区| 欧美激情国产日韩精品一区| 亚洲乱码一区二区免费版| 亚洲av成人av| 国产老妇伦熟女老妇高清| 真实男女啪啪啪动态图| 亚洲电影在线观看av| 国产私拍福利视频在线观看| 久久国内精品自在自线图片| 国产色婷婷99| 成人毛片a级毛片在线播放| 国产高清视频在线观看网站| 可以在线观看的亚洲视频| 联通29元200g的流量卡| 亚洲精品日韩在线中文字幕 | av在线亚洲专区| 国产极品天堂在线| 天美传媒精品一区二区| 欧洲精品卡2卡3卡4卡5卡区| 国产精品一区二区在线观看99 | 特级一级黄色大片| 免费观看a级毛片全部| 嫩草影院新地址| 蜜桃亚洲精品一区二区三区| 少妇熟女aⅴ在线视频| 深夜a级毛片| 校园人妻丝袜中文字幕| 国产乱人偷精品视频| 天天一区二区日本电影三级| 欧美三级亚洲精品| 国产精品嫩草影院av在线观看| 日韩精品有码人妻一区| 一区二区三区四区激情视频 | 国产高潮美女av| 精品午夜福利在线看| 禁无遮挡网站| 久久亚洲国产成人精品v| 久久99精品国语久久久| 国产黄片视频在线免费观看| 老司机福利观看| 少妇高潮的动态图| 国产精品免费一区二区三区在线| 国产淫片久久久久久久久| 联通29元200g的流量卡| 不卡一级毛片| 久久久午夜欧美精品| 少妇猛男粗大的猛烈进出视频 | 成年女人看的毛片在线观看| 精品国内亚洲2022精品成人| 欧美高清成人免费视频www| 国产精品99久久久久久久久| 亚洲性久久影院| 啦啦啦观看免费观看视频高清| 国产伦理片在线播放av一区 | 日韩高清综合在线| 久久久久久久亚洲中文字幕| 日韩人妻高清精品专区| 国产视频首页在线观看| 久久精品久久久久久久性| 国产成人午夜福利电影在线观看| 一级黄色大片毛片| 91久久精品国产一区二区成人| 亚洲国产精品成人综合色| 人人妻人人看人人澡| 日韩视频在线欧美| 在线播放无遮挡| 国产精品人妻久久久久久| 老女人水多毛片| 精品欧美国产一区二区三| 男女视频在线观看网站免费| 久久久久九九精品影院| 欧美成人一区二区免费高清观看| 熟女人妻精品中文字幕| 久久草成人影院| 中文字幕制服av| 欧美+日韩+精品| 国产 一区精品| 91在线精品国自产拍蜜月| 欧美成人精品欧美一级黄| 内射极品少妇av片p| www日本黄色视频网| 久久韩国三级中文字幕| 婷婷色av中文字幕| 小说图片视频综合网站| 国产男人的电影天堂91| 国产片特级美女逼逼视频| 麻豆一二三区av精品| 亚洲真实伦在线观看| 97超视频在线观看视频| 美女国产视频在线观看| 日本熟妇午夜| 国产伦一二天堂av在线观看| 久久鲁丝午夜福利片| 三级经典国产精品| 日韩欧美 国产精品| 精华霜和精华液先用哪个| 亚洲内射少妇av| 欧美成人一区二区免费高清观看| 欧美3d第一页| 99久久九九国产精品国产免费| 国产蜜桃级精品一区二区三区| 午夜免费激情av| 九九爱精品视频在线观看| 大又大粗又爽又黄少妇毛片口| 亚州av有码| 免费观看人在逋| 狂野欧美白嫩少妇大欣赏| 日韩视频在线欧美| 久久精品91蜜桃| 爱豆传媒免费全集在线观看| 在线观看午夜福利视频| 欧美成人免费av一区二区三区| 成人永久免费在线观看视频| ponron亚洲| 久久精品综合一区二区三区| 最近2019中文字幕mv第一页| av在线老鸭窝| 夫妻性生交免费视频一级片| 久久精品国产自在天天线| 超碰av人人做人人爽久久| 国模一区二区三区四区视频| 日韩 亚洲 欧美在线| 久久草成人影院| 亚洲人成网站高清观看| 成人高潮视频无遮挡免费网站| 亚洲中文字幕一区二区三区有码在线看| 色综合亚洲欧美另类图片| 欧美性感艳星| 免费看日本二区| 亚洲人与动物交配视频| 99热全是精品| 亚洲在线观看片| 人人妻人人澡人人爽人人夜夜 | 亚洲精品日韩在线中文字幕 | 欧美高清性xxxxhd video| 久久久久久久久中文| 波野结衣二区三区在线| 久久午夜福利片| 亚洲av成人av| 久久精品夜夜夜夜夜久久蜜豆| 精品久久久久久久久av| 男人的好看免费观看在线视频| 成人av在线播放网站| 亚洲av成人精品一区久久| 精品国内亚洲2022精品成人| 黄色配什么色好看| 有码 亚洲区| 蜜桃亚洲精品一区二区三区| 三级国产精品欧美在线观看| 亚洲精品影视一区二区三区av| 国产乱人偷精品视频| 午夜福利在线观看免费完整高清在 | 欧美日韩在线观看h| 一边摸一边抽搐一进一小说| 亚洲性久久影院| 亚洲一级一片aⅴ在线观看| 国产人妻一区二区三区在| 久久久久久久久久成人| 少妇丰满av| 好男人在线观看高清免费视频| 欧美高清成人免费视频www| 国产高清激情床上av| 国产亚洲av嫩草精品影院| 国产伦理片在线播放av一区 | 中文字幕免费在线视频6| 麻豆久久精品国产亚洲av| 91在线精品国自产拍蜜月| 中文字幕久久专区| 久久久久久久久大av| 亚洲成a人片在线一区二区| 熟女电影av网| 深爱激情五月婷婷| 天堂网av新在线| 日韩欧美精品v在线| 能在线免费看毛片的网站| 国产成人午夜福利电影在线观看| 亚洲国产欧美在线一区| 亚洲美女搞黄在线观看| 边亲边吃奶的免费视频| 在线观看午夜福利视频| 国产午夜福利久久久久久| 亚洲七黄色美女视频| 精品久久久久久久人妻蜜臀av| 男女做爰动态图高潮gif福利片| 久久午夜亚洲精品久久| 日本熟妇午夜| 丰满人妻一区二区三区视频av| 男人和女人高潮做爰伦理| 亚洲av免费在线观看| 免费看光身美女| 啦啦啦观看免费观看视频高清| www.av在线官网国产| av免费观看日本| 免费不卡的大黄色大毛片视频在线观看 | 国内精品宾馆在线| 国产精品久久久久久久久免| 99久久精品一区二区三区| 国产精品.久久久| 久久精品国产自在天天线| 亚洲av一区综合| 亚洲av不卡在线观看| 国语自产精品视频在线第100页| 日韩欧美一区二区三区在线观看| 禁无遮挡网站| 国产午夜精品一二区理论片| 亚洲精品乱码久久久久久按摩| 国产黄a三级三级三级人| 三级毛片av免费| www.av在线官网国产| 久久精品综合一区二区三区| 日本一本二区三区精品| 久久精品国产自在天天线| 好男人在线观看高清免费视频| 日本熟妇午夜| 波多野结衣高清作品| 亚洲av不卡在线观看| 国内少妇人妻偷人精品xxx网站| av视频在线观看入口| 99热精品在线国产| 免费电影在线观看免费观看| 国产一级毛片在线| 国产三级中文精品| 一本一本综合久久| .国产精品久久| 淫秽高清视频在线观看| 国产高清三级在线| 日日撸夜夜添| 成人亚洲精品av一区二区| 男人舔奶头视频| 亚洲四区av| 欧美激情久久久久久爽电影| 欧美一区二区亚洲| 国产午夜精品一二区理论片| 亚州av有码| 亚洲av不卡在线观看| 成年女人看的毛片在线观看| 亚洲精品国产成人久久av| 蜜桃久久精品国产亚洲av| 国内精品美女久久久久久| 男女边吃奶边做爰视频| 午夜久久久久精精品| 亚洲乱码一区二区免费版| 国产精品美女特级片免费视频播放器| 成人毛片60女人毛片免费| av卡一久久| 中国美白少妇内射xxxbb| 99在线人妻在线中文字幕| 色尼玛亚洲综合影院| 亚洲欧美精品专区久久| 亚洲成a人片在线一区二区| 一区二区三区高清视频在线| 1024手机看黄色片| 两个人的视频大全免费| 午夜福利视频1000在线观看| 三级男女做爰猛烈吃奶摸视频| 偷拍熟女少妇极品色| 国产精品野战在线观看| 欧美精品国产亚洲| 国产乱人偷精品视频| 3wmmmm亚洲av在线观看| 国产单亲对白刺激| 91aial.com中文字幕在线观看| 熟女人妻精品中文字幕| ponron亚洲| 毛片女人毛片| 午夜福利在线观看免费完整高清在 | 99热只有精品国产| 边亲边吃奶的免费视频| 久久人人精品亚洲av| 听说在线观看完整版免费高清| 亚洲自偷自拍三级| 亚洲国产精品成人综合色| 日韩欧美国产在线观看| 午夜福利在线在线| 亚洲成人久久爱视频| 波多野结衣巨乳人妻| 欧美日本亚洲视频在线播放| 女人十人毛片免费观看3o分钟| 日韩视频在线欧美| eeuss影院久久| 一边摸一边抽搐一进一小说| 99热这里只有是精品50| 国产高清激情床上av| 国产淫片久久久久久久久| 国产伦理片在线播放av一区 | 岛国毛片在线播放| 免费观看a级毛片全部| 少妇猛男粗大的猛烈进出视频 | 国产精品av视频在线免费观看| 3wmmmm亚洲av在线观看| 1000部很黄的大片| 国产老妇女一区| 欧美在线一区亚洲| 亚洲国产欧美人成| 国产黄片视频在线免费观看| 国产亚洲av嫩草精品影院| 午夜激情福利司机影院| 尤物成人国产欧美一区二区三区| 中文字幕久久专区| 亚洲av中文av极速乱| 卡戴珊不雅视频在线播放| 成人美女网站在线观看视频| av天堂在线播放| 哪里可以看免费的av片| 国产高清激情床上av| 91狼人影院| 日韩在线高清观看一区二区三区| 精品少妇黑人巨大在线播放 | 22中文网久久字幕| 国产av在哪里看| 日日撸夜夜添| 观看美女的网站| 亚洲精品影视一区二区三区av| 少妇熟女欧美另类| 丝袜喷水一区| 亚洲人成网站在线播放欧美日韩| 国产精品福利在线免费观看| 99热只有精品国产| 天堂√8在线中文| a级毛片a级免费在线| 国产成人a∨麻豆精品| 亚洲成a人片在线一区二区| 天堂影院成人在线观看| 亚洲欧美成人精品一区二区| 亚洲久久久久久中文字幕| 啦啦啦啦在线视频资源| 亚洲精品乱码久久久v下载方式| 在线播放国产精品三级| 高清毛片免费看| 中文字幕人妻熟人妻熟丝袜美| av黄色大香蕉| 欧美潮喷喷水| 成人午夜精彩视频在线观看| 国产欧美日韩精品一区二区| 国产午夜精品论理片| 欧美不卡视频在线免费观看| 久久久久久久久久久免费av| 午夜福利高清视频| av免费观看日本| 长腿黑丝高跟| 卡戴珊不雅视频在线播放| 国产精品久久久久久亚洲av鲁大| 成人三级黄色视频| 国产精品,欧美在线| 狂野欧美激情性xxxx在线观看| 蜜桃亚洲精品一区二区三区| 久久精品国产亚洲网站| 观看免费一级毛片| 日韩,欧美,国产一区二区三区 | 偷拍熟女少妇极品色| 美女 人体艺术 gogo| 春色校园在线视频观看| 国产av不卡久久| 夜夜爽天天搞| 男人狂女人下面高潮的视频| av专区在线播放| 国模一区二区三区四区视频| 欧美一区二区国产精品久久精品| 成人高潮视频无遮挡免费网站| 国产成年人精品一区二区| 日日干狠狠操夜夜爽| 少妇熟女欧美另类| 久久久久免费精品人妻一区二区| av福利片在线观看| 黄片wwwwww| 成人漫画全彩无遮挡| 国产亚洲精品久久久com| 亚洲图色成人| 日本熟妇午夜| 男人舔女人下体高潮全视频| 女人十人毛片免费观看3o分钟| 性欧美人与动物交配| 国产亚洲精品av在线| 波多野结衣高清无吗| 淫秽高清视频在线观看| 乱系列少妇在线播放| 99热精品在线国产| 精品人妻偷拍中文字幕| 欧美在线一区亚洲| 亚洲欧美精品综合久久99| 神马国产精品三级电影在线观看| 蜜桃久久精品国产亚洲av| 免费看a级黄色片| 啦啦啦啦在线视频资源| 在线天堂最新版资源| 国产成人影院久久av| 中文字幕av在线有码专区| 色吧在线观看| 99久久久亚洲精品蜜臀av| 日韩av在线大香蕉| 日本熟妇午夜| 国产视频首页在线观看| 国产成人freesex在线| 综合色av麻豆| 日韩一区二区视频免费看| 丰满的人妻完整版| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美成人综合另类久久久 | 国产私拍福利视频在线观看| av在线老鸭窝| 国产成人91sexporn| 成人二区视频| 在线天堂最新版资源| 久久九九热精品免费| 三级男女做爰猛烈吃奶摸视频| 亚洲七黄色美女视频| 国内精品久久久久精免费| 久久精品久久久久久噜噜老黄 | 久久久久性生活片| 午夜福利在线观看免费完整高清在 | 国产精品精品国产色婷婷| 在线播放无遮挡| 亚洲国产精品成人综合色| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品粉嫩美女一区| 亚洲性久久影院| 综合色av麻豆| 天堂av国产一区二区熟女人妻| 久久久精品欧美日韩精品| 99热精品在线国产| 国产成人精品婷婷| 国产精品蜜桃在线观看 | 亚洲精品国产av成人精品| 亚洲欧美日韩东京热| 狂野欧美激情性xxxx在线观看| 看免费成人av毛片| 免费大片18禁| 亚洲精华国产精华液的使用体验 | 国产视频内射| 成人永久免费在线观看视频| 久久综合国产亚洲精品| 欧美高清性xxxxhd video| 美女黄网站色视频| 又爽又黄a免费视频| 亚洲综合色惰| 国产伦精品一区二区三区四那| 又粗又硬又长又爽又黄的视频 | 成人美女网站在线观看视频| 舔av片在线| 日韩欧美在线乱码| 色哟哟·www| 亚洲在久久综合| 久久综合国产亚洲精品| 亚洲国产高清在线一区二区三| 日韩高清综合在线| 久久精品91蜜桃| 国产日本99.免费观看| 此物有八面人人有两片| 一本一本综合久久| 亚洲av中文av极速乱| av女优亚洲男人天堂| 12—13女人毛片做爰片一| 性插视频无遮挡在线免费观看| 我要看日韩黄色一级片| 亚洲av电影不卡..在线观看| 天堂网av新在线| 久久精品影院6| 99热精品在线国产| 哪个播放器可以免费观看大片| 最近最新中文字幕大全电影3| 国产一级毛片七仙女欲春2| 午夜福利视频1000在线观看| 直男gayav资源| 欧美精品一区二区大全| 久久午夜亚洲精品久久| 亚洲精品日韩在线中文字幕 | 中文字幕人妻熟人妻熟丝袜美| 欧美性猛交黑人性爽| 国产亚洲av片在线观看秒播厂 | 国产久久久一区二区三区| 精品国产三级普通话版| 欧洲精品卡2卡3卡4卡5卡区| 草草在线视频免费看| 免费一级毛片在线播放高清视频| 麻豆国产97在线/欧美| 国产成人精品一,二区 | 成人鲁丝片一二三区免费| 亚洲三级黄色毛片| 国产精品综合久久久久久久免费| 国产黄片视频在线免费观看| a级一级毛片免费在线观看| 日韩欧美一区二区三区在线观看| 亚洲久久久久久中文字幕| 九九久久精品国产亚洲av麻豆| 晚上一个人看的免费电影| www日本黄色视频网| 啦啦啦韩国在线观看视频| 九色成人免费人妻av| 亚洲人与动物交配视频| 亚洲精品久久久久久婷婷小说 | 91久久精品电影网| 欧美性猛交╳xxx乱大交人| 欧美日本视频| 听说在线观看完整版免费高清| 国产伦在线观看视频一区| 18禁在线播放成人免费| 桃色一区二区三区在线观看| 五月伊人婷婷丁香| 蜜桃亚洲精品一区二区三区| 国产一区二区在线观看日韩| 国产成人aa在线观看| 日韩人妻高清精品专区| 三级毛片av免费| 高清午夜精品一区二区三区 | 亚洲第一区二区三区不卡| 国产精品嫩草影院av在线观看| 久久国产乱子免费精品| 中文字幕精品亚洲无线码一区| 久久久成人免费电影| 禁无遮挡网站| 免费人成在线观看视频色| 亚洲一级一片aⅴ在线观看| 国产又黄又爽又无遮挡在线| 免费大片18禁| 亚洲欧美日韩东京热| 亚洲人与动物交配视频| 国产爱豆传媒在线观看| 桃色一区二区三区在线观看| 日韩欧美精品v在线| 欧美另类亚洲清纯唯美| 日韩制服骚丝袜av| 十八禁国产超污无遮挡网站| 22中文网久久字幕| 啦啦啦啦在线视频资源| 日产精品乱码卡一卡2卡三| 国产高清有码在线观看视频| 久久亚洲国产成人精品v| 欧美一区二区国产精品久久精品| 亚洲性久久影院| www日本黄色视频网| 好男人在线观看高清免费视频| ponron亚洲| 蜜桃久久精品国产亚洲av| 久久精品国产亚洲av香蕉五月| a级一级毛片免费在线观看| 中文字幕熟女人妻在线| 你懂的网址亚洲精品在线观看 | 五月伊人婷婷丁香| 天堂av国产一区二区熟女人妻| 舔av片在线| 成人毛片60女人毛片免费| 亚洲一级一片aⅴ在线观看| 婷婷亚洲欧美| 国产蜜桃级精品一区二区三区| 久久久精品大字幕| av.在线天堂| 亚洲av电影不卡..在线观看|