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

    光學微腔中倍頻光場演化和光譜特性*

    2020-02-18 03:17:48徐昕金雪瑩胡曉鴻黃新寧
    物理學報 2020年2期
    關(guān)鍵詞:微腔光場倍頻

    徐昕 金雪瑩? 胡曉鴻 黃新寧

    1) (合肥工業(yè)大學,儀器科學與光電工程學院,合肥 230009)

    2) (中國科學院西安光學精密機械研究所,瞬態(tài)光學與光子技術(shù)國家重點實驗室,西安 710119)

    在考慮光學微腔中二階和三階非線性效應的情況下,引入了可同時描述腔內(nèi)基頻和倍頻光場的演化過程的Lugiato-Lefeve方程,分析了SiN微腔中二次諧波的產(chǎn)生,并討論了各參數(shù)對腔內(nèi)基頻和倍頻光場的影響.理論分析結(jié)果表明,失諧參量為0時,穩(wěn)定后的基頻光場為平頂脈沖的形式,而倍頻光場呈正弦分布;失諧參量增加,將導致腔內(nèi)基頻和倍頻光功率在演化過程中出現(xiàn)振蕩,且最終穩(wěn)定的光功率變?nèi)?穩(wěn)定后的光場分布為周期性變化;失諧參量的值過大,會使得微腔光場處于混沌狀態(tài).抽運光強較弱時,腔內(nèi)可形成穩(wěn)定的光場分布;抽運光強較強時,會導致腔內(nèi)色散以及非線性效應過強,最終穩(wěn)定的光場仍然呈周期性變化,且抽運光功率越強,光功率的演化曲線振蕩越強.此外,選取特定的微腔尺寸,微腔可工作于“圖靈環(huán)”狀態(tài).理論分析結(jié)果對研究光學微腔中二次諧波的產(chǎn)生有重要意義.

    1 引 言

    光頻梳技術(shù)已經(jīng)廣泛地應用于絕對距離測量[1]、光譜測量[2]、光學頻率測量[3]等領(lǐng)域.目前,光頻梳主要是基于鎖模激光器實現(xiàn)的,由于受到鎖模激光器腔長的限制,它們的重復頻率一般為幾十到幾百MHz[4,5].但是,隨著光頻梳應用領(lǐng)域的不斷拓展,需要重復頻率更高的光頻梳,例如:在高分辨率天文光譜儀的波長定標中,光頻梳的重復頻率要達到10—30 GHz[6].傳統(tǒng)的基于鎖模激光器的光頻梳已經(jīng)無法滿足需求,因此,一種新型的基于光學微腔的光頻梳在近些年成為了研究熱點[7?9].

    基于微腔的光頻梳是利用級聯(lián)的四波混頻效應,在光學微腔中產(chǎn)生一系列等間距的頻率成分,從而將單一的抽運頻率擴展到一個非常超寬光譜范圍.微腔光頻梳具有體積小、結(jié)構(gòu)簡單、易實現(xiàn)光學集成等優(yōu)點.通過理論研究證實,微腔可工作于多種模式,例如:光孤子、圖靈環(huán)、暗孤子、混沌等模式[10],并且這些工作模式也在實驗中被觀察到[11?15].目前,在光學微腔的理論研究中,大多都是以Lugiato-Lefeve 方程 (Lugiato-Lefeve equation,LLE)為理論模型,通過求解LLE,得到穩(wěn)態(tài)時腔內(nèi)的光場分布[16].關(guān)于光學微腔中二次諧波的產(chǎn)生,已經(jīng)在實驗中被觀察到.Xue等[17]通過四波混頻的方法在正色散的光學微腔中同時觀察到了Kerr光頻梳和二次諧波的產(chǎn)生,此外還建立了光學微腔中二次諧波產(chǎn)生的理論模型,理論分析和實驗結(jié)果相一致.Lin等[18]也在片上集成的光學微腔中實現(xiàn)了高效率的二次諧波產(chǎn)生.

    然而,目前對微腔光頻梳的理論研究采用的LLE,只包含了腔內(nèi)的損耗、色散、抽運項,并沒有考慮微腔中的非線性效應[19].但是由于光學微腔本身體積很小,腔內(nèi)光功率密度較高,因而腔內(nèi)的非線性效應是無法被忽略的.當考慮微腔內(nèi)的二階和三階非線性效應時,腔內(nèi)會由于非線性效應產(chǎn)生二次諧波[20].因此,本文將在二階、三階非線性效應的基礎上,分析光學微腔內(nèi)二次諧波的產(chǎn)生情況,并研究各個參數(shù)對腔內(nèi)基頻光和倍頻光光場產(chǎn)生的影響.

    2 理論模型

    本文討論以SiN為材料的光學微腔.SiN材料的光學微腔一般為以Si為基底的片上結(jié)構(gòu),結(jié)構(gòu)如圖1所示,抽運光通過微腔一側(cè)的波導結(jié)構(gòu)耦合進入微腔,光在微腔傳播過程中耦合進入另一側(cè)的波導結(jié)構(gòu)輸出.在這一過程中,由于二階、三階非線性效應,腔內(nèi)會產(chǎn)生二次諧波.為討論微腔中的二次諧波產(chǎn)生情況,引入包含二次諧波項的LLE來描述腔內(nèi)基頻和倍頻光場在時間和空間上的演化,分別用E1表示基頻光的光場,E2表示倍頻光的光場,基頻和倍頻光場隨時間和空間的變化關(guān)系滿足以下方程[17]:

    其中,z為光沿微腔傳播的距離;τ是光傳播的時間,即快時間;α1和 α2分別為基頻光和倍頻光在腔內(nèi)傳播單位長度的損耗;參量 δ1=(ω0– ωp)tR/L為抽運光頻率ωp和微腔諧振頻率ω0之間的失諧參量,L 為微腔的周長,tR=L/c 為基頻光在腔內(nèi)循環(huán)一周的時間;和分別表示基頻光和倍頻光的二階色散,它們的表達式為

    Δk和Δ k′分別為相位失配和群速度失配,它們通過以下公式定義:

    γ1和 γ2分別為自相位調(diào)制的非線性系數(shù);γ12和γ21是交叉相位調(diào)制的非線性系數(shù);κ為二階耦合系數(shù);η1為抽運光和腔內(nèi)光場的耦合系數(shù);抽運光為連續(xù)光,強度用 Ein表示.此外,對于光學微腔還有一個重要參數(shù),即微腔的自由光譜范圍(free spectral range,FSR),FSR=1/tR,它也表示了微腔光頻梳中梳齒之間的頻率間隔.(1)式和(2)式中各色散、非線性、耦合等參數(shù)是由光學微腔的材料和結(jié)構(gòu)決定的.

    通過求解(1)式和(2)式,可以得到微腔內(nèi)基頻和倍頻光場的分布情況.為了能夠清楚腔內(nèi)光場分布隨時間的演化過程,采用分布傅里葉法求解(1)式和(2)式[19].

    圖1 SiN 微腔結(jié)構(gòu)示意圖Fig.1.Structure of the SiN microresonator.

    3 理論分析

    3.1 光學微腔中二次諧波的產(chǎn)生

    以SiN為材料的光學微腔為正色散.SiN材料的各項參數(shù)如下:

    假設微腔直徑為 100 μm,對應的 FSR=231.3 GHz.微腔采用波長為1550 nm的連續(xù)光抽運.微腔內(nèi)的初始光場是由調(diào)制不穩(wěn)定性產(chǎn)生的,根據(jù)文獻[21]可以假設腔內(nèi)初始的基頻光場E10為高斯脈沖,如由于初始光場源自于腔內(nèi)的調(diào)制不穩(wěn)定性,因此初始光場的強度很弱,因此假設初始光場的振幅 E0=0.1 W1/2,脈沖寬度τ0=0.4 ps.

    首先,忽略抽運頻率和微腔本征頻率之間的頻率失諧,即失諧參量 δ1=0,連續(xù)抽運光的強度Ein=1 W1/2.在非線性效應的作用下,腔內(nèi)有二次諧波產(chǎn)生,二次諧波隨著光在腔內(nèi)循環(huán)次數(shù)產(chǎn)生的光場分布如圖2(a)所示,可見初始產(chǎn)生的倍頻光場較弱,隨著循環(huán)次數(shù)的增加,腔內(nèi)的倍頻光場不斷變強.圖2(b)和(c)表示腔內(nèi)基頻光和倍頻光功率隨循環(huán)次數(shù)不斷增加,并最終穩(wěn)定于某一值,此時腔內(nèi)光場相對穩(wěn)定,穩(wěn)定后基頻光和倍頻光光場分布如圖2(d)和(e)所示,基頻光在腔內(nèi)形成了一個較寬的脈沖,且脈沖頂端較為平坦而倍頻光場則類似正弦信號,其功率和基頻光相比相對較弱.圖2(f)和(g)是對應的光譜,基頻光和倍頻光都在微腔內(nèi)形成了梳狀的光譜,但倍頻光譜寬要小于基頻光的譜寬.

    以上是失諧參數(shù)δ1=0時腔內(nèi)光場的演化情況.下面討論失諧參數(shù) δ1不為 0 的情況.假設 δ1=0.1,抽運光強 Ein=1 W1/2.通過求解 (1)式和(2)式,得到微腔內(nèi)基頻光和倍頻光功率變化如圖3(a)和(b)所示.對比圖2(b)和(c)可以發(fā)現(xiàn),當微腔的抽運頻率和本征頻率發(fā)生失諧時,腔內(nèi)的基頻光和倍頻光功率在光循環(huán)過程中發(fā)生振蕩,且經(jīng)過多次振蕩后腔內(nèi)最終穩(wěn)定的光功率和圖2(b)和(c)中的光功率相比,強度明顯減弱,說明頻率失諧會對腔內(nèi)的光功率造成明顯的損耗.當腔內(nèi)光功率穩(wěn)定后,腔內(nèi)基頻光和倍頻光光場分布隨循環(huán)次數(shù)的演化如圖3(c)和(d)所示,其中任意時刻光場的分布如圖3(e)和(f)所示.可見,基頻光在腔內(nèi)最終形成了穩(wěn)定的暗孤子,而倍頻光在腔內(nèi)的光場仍然近似正弦信號,但正弦信號的相位隨著循環(huán)次數(shù)而變化.

    圖2 SiN 光學微腔中二次諧波的產(chǎn)生 (δ1=0,Ein=1 W1/2) (a) 倍頻光場在腔內(nèi)的演化過程;(b) 基頻光功率隨光在腔內(nèi)循環(huán)次數(shù)的變化曲線;(c) 倍頻光功率隨光在腔內(nèi)循環(huán)次數(shù)的變化曲線;(d)基頻場的穩(wěn)定分布;(e)倍頻光場的穩(wěn)定分布;(f)基頻光的光譜;(g) 基頻光的光譜Fig.2.Second harmonic waves generates inside the SiN microresonator (δ1=0,Ein=1 W1/2):(a) Evolution of the second harmonic waves;(b) curves of the fundamental frequency waves power and (c) the second frequency waves power vary with the round trip number;(d) stationary distribution of the fundamental frequency waves light field and (e) the second frequency waves light field;(f) spectra of the fundamental frequency waves and (g) the second frequency waves.

    圖3 發(fā)生頻率失諧后,微腔內(nèi)光場變化情況 (δ1=0.1,Ein=1 W1/2) (a) 基頻光功率和 (b)倍頻光功率隨光在腔內(nèi)循環(huán)次數(shù)的變化曲線;光場穩(wěn)定后,(c)基頻光場和(d)倍頻光場隨時間的演化;(e) 基頻光場和(f)倍頻廣場的穩(wěn)定分布Fig.3.Light field evolution in the microresonator with the frequency detuning (δ1=0.1 and Ein=1 W1/2):(a) Curves of the fundamental frequency waves power and (b) the second frequency waves power vary with the round trip number;evolution of (c) the fundamental frequency waves (d) the second frequency waves after the light fields are stable;stationary distribution of (e) the fundamental frequency waves light field and (f) the second frequency waves light field.

    通過以上分析可知,微腔的頻率失諧參數(shù)δ1會引起腔內(nèi)基頻光和倍頻光功率的損耗,并導致腔內(nèi)光場分布發(fā)生變化.因此,接下來將分析失諧參數(shù)δ1對微腔內(nèi)二次諧波產(chǎn)生的影響.

    3.2 失諧參量對微腔內(nèi)二次諧波的影響

    假設抽運光場 Ein=1 W1/2不變,依次增加失諧參數(shù)δ1,得到微腔內(nèi)基頻光和倍頻光功率變化曲線如圖4 所示.圖4(a)和 (b)是 δ1取不同值時,基頻光在腔內(nèi)循環(huán)過程中的功率變化曲線,隨著失諧參數(shù)δ1的不斷增加,基頻光的功率在趨于穩(wěn)定的過程中振蕩越來越明顯,當δ1大于0.2后(圖4(b)),功率的振蕩變得非常劇烈.因此,抽運光頻率和微腔本征頻率之間的失諧會導致腔內(nèi)光場分布的不穩(wěn)定.圖4(c)和(d)是失諧參數(shù)δ1對腔內(nèi)倍頻光的影響.和基頻光的變化趨勢一致,倍頻光功率在腔內(nèi)經(jīng)歷先增大再減小的變化過程,隨著δ1的增加,曲線下降沿的抖動也越來越劇烈,腔內(nèi)最終穩(wěn)定存在的倍頻光功率總體呈下降趨勢.可以看到,δ1=1 時,微腔內(nèi)基頻光和倍頻光的功率變化曲線都出現(xiàn)了周期性的抖動,腔內(nèi)光場趨于穩(wěn)定后,光功率仍然保持周期性振蕩.圖5(a)和(b)給出了腔內(nèi)基頻和倍頻光場分布的演化,光場分布隨腔內(nèi)光場循環(huán)的次數(shù)進行周期性的變化.在光場變化周期中,每隔100次循環(huán)選取一個時間點,畫出這些時間點腔內(nèi)的光場分布的形狀,結(jié)果如圖5(c)和(d)所示.由圖5(c)可見,在某些時刻基頻光在腔內(nèi)可形成單脈沖形式的分布(圖5(c)②③),但也仍有脈沖分裂的情況存在(圖5(c)④⑤).圖5(d)為倍頻光在各個時刻的光場分布情況,可以看出,當失諧參量δ1=1時倍頻光在腔內(nèi)仍然近似正弦分布,正弦波的相位隨時間不斷變化.

    圖4 失諧參量 δ1 取不同值時,微腔內(nèi)基頻光和倍頻光功率變化曲線 (a) 0.02 ≤ δ1 ≤ 0.08 時,基頻光功率變化曲線;(b) 0.2 ≤δ1 ≤ 0.8 時,基頻光功率變化曲線;(c) 0.02 ≤ δ1 ≤ 0.08 時,倍頻光功率變化曲線;(d) 0.2 ≤ δ1 ≤ 0.8 時,倍頻光功率變化曲線Fig.4.Influence of the frequency detuning δ1 on the power change curves:(a) Curves of the power variation for the fundamental frequency waves,0.02 ≤ δ1 ≤ 0.08;(b) curves of the power variation for the fundamental frequency waves,0.2 ≤ δ1 ≤ 0.8;(c) curves of the power variation for the second harmonic waves,0.02 ≤ δ1 ≤ 0.08;(d) curves of the power variation for the second harmonic waves,0.2 ≤ δ1 ≤ 0.8.

    研究過程中,我們進一步增大了失諧參量δ1的取值,結(jié)果發(fā)現(xiàn)當 δ1的值過大時,腔內(nèi)的基頻光和倍頻光的光場均呈現(xiàn)無規(guī)則的變化,如圖6(a)和(b)所示.經(jīng)過光場在腔內(nèi)的長時間演化,微腔中的光場在時間和空間上都沒有形成穩(wěn)定分布,此時微腔內(nèi)處于混沌狀態(tài)[22].任取某一時刻,分別繪制基頻光和倍頻光的光場分布,如圖6(c)和(d)所示,圖6(e)和 (f)分別為其對應的光譜.可見,混沌模式的光場表現(xiàn)為光功率隨時間無規(guī)律地變化,其相對應的光譜也是無規(guī)則的,每一個梳齒對應的功率也是時刻發(fā)生變化的.產(chǎn)生混沌的原因為失諧參量δ1的值過大,導致腔內(nèi)光場的分布不穩(wěn)定,在色散以及高階非線性效應的作用下,腔內(nèi)開始隨機地出現(xiàn)毛刺以及功率較高的尖峰[23].因此,在利用光學微腔產(chǎn)生二次諧波光梳的過程中,要盡量避免失諧參量δ1過大,即抽運頻率要盡量接近微腔的本征諧振頻率,以免在腔內(nèi)產(chǎn)生混沌狀態(tài).

    失諧參量δ1是抽運光頻率和微腔諧振頻率之間的失諧,δ1的值越大,表示抽運與微腔諧振之間的頻率失諧越嚴重.這種失諧會導致腔內(nèi)光場不同的模式之間發(fā)生相互耦合,因而使腔內(nèi)基頻和倍頻光場的不穩(wěn)定.特別是當δ1值過大時,腔內(nèi)光場的模式發(fā)生混亂,因而產(chǎn)生不規(guī)則的光場分布.根據(jù)以上分析,δ1的值越小,不僅可以有效避免混沌的現(xiàn)象產(chǎn)生,同時還可獲得較強的倍頻光功率以及腔內(nèi)穩(wěn)定的光場分布.

    圖5 失諧參量 δ1=1 時,腔內(nèi)光場穩(wěn)定后基頻光和倍頻光的光場演化 (a) 基頻光場的演化;(b) 倍頻光場的演化;(c) 基頻光光場變化周期內(nèi),光在腔內(nèi)每循環(huán)100次,繪制其波形;(d) 倍頻光光場變化周期內(nèi),光在腔內(nèi)每循環(huán)100次,繪制其波形Fig.5.Stable evolution of the dual light fields when δ1=1:(a) Evolution of the fundamental frequency waves;(b) evolution of the second harmonic waves;(c) intensity profiles of the fundamental frequency waves at six different moments within a period,the waveforms are plotted every hundred times;(d) intensity profiles of the second harmonic waves at six different moments within a period,the waveforms are plotted every hundred times.

    圖6 失諧參量 δ1 取值過大時,腔內(nèi)光場的混沌狀態(tài) (a) 基頻光光場的混沌狀態(tài);(b) 倍頻光光場的混沌狀態(tài);(c) 某一時刻基頻光場的分布;(d) 某一時刻倍頻光場的分布;(e) 與圖(c)中光場對應的基頻光光譜;(f) 與圖(d)中光場對應的倍頻光光譜Fig.6.Chaos inside the microresonator,when the value of detuning parameter is too large:(a) Chaos of the fundamental frequency waves;(b) chaos of the second harmonic waves;(c) intensity profile of the fundamental frequency waves in a moment;(d) intensity profile of the second harmonic waves in a moment;(e) spectrum of the fundamental frequency waves;(f) spectrum of the second harmonic waves.

    3.3 抽運功率對微腔內(nèi)二次諧波的影響

    微腔的抽運功率同樣也會對腔內(nèi)基頻光和倍頻光的功率以及光場分布產(chǎn)生影響.我們保持失諧參量 δ1=0.1 不變,在此基礎上,不斷增加抽運強度Ein,得到微腔內(nèi)基頻光和倍頻光功率變化曲線分別如圖7(a)和(b)所示.可明顯看出,隨著抽運功率的不斷增大,腔內(nèi)最終穩(wěn)定后光功率也會越大,且抽運功率的大小幾乎不會影響腔內(nèi)形成穩(wěn)定光場所消耗的時間.因此,為了獲得更高功率的倍頻光,可以增強抽運光的功率.

    然而抽運功率的增大也會導致腔內(nèi)光場的不穩(wěn)定.將抽運光場強度 Ein增大到 100 W1/2,腔內(nèi)基頻光和倍頻光功率變化曲線如圖8(a)和(b)所示,可見,腔內(nèi)的光功率經(jīng)過大幅振蕩衰減后趨于穩(wěn)定,穩(wěn)定時期的光功率仍然保持周期性的小幅振蕩.當光功率穩(wěn)定振蕩后,基頻光和倍頻光的光場演化如圖8(c)和(d)所示,光場分布同樣呈周期性變化,并且可以認為在任意時刻,基頻光和倍頻光的光場在腔內(nèi)是均勻分布的.

    在此基礎上進一步增大抽運光強度,抽運光場 Ein分別為 800,1000和1200 W1/2時,基頻光和倍頻光功率變化曲線如圖9所示.當 Ein=800 W1/2時,基頻光的功率在腔內(nèi)經(jīng)過功率振蕩后,最終大約穩(wěn)定于20 kW,而倍頻光功率在循環(huán)過程中增加,達到峰值后迅速衰減,然后趨向于穩(wěn)定,穩(wěn)定部分的光功率也同樣保持周期性振蕩.這是由于抽運光功率過高,導致腔內(nèi)二階非線性效應明顯,較多的基頻光轉(zhuǎn)化為倍頻光,因而腔內(nèi)累積的倍頻光功率越來越大,因此光在腔內(nèi)循環(huán)過程中,倍頻光的功率有明顯的增長.當倍頻光功率達到一定值時,由頻率失諧、色散、損耗引起的倍頻光功率降低,超過了由非線性效應產(chǎn)生的倍頻光強度,因而倍頻光在功率增大到一定程度后會出現(xiàn)功率迅速衰減的現(xiàn)象,所以在圖9(b)中倍頻光的功率會出現(xiàn)一個脈沖形式的變化曲線.由于倍頻光是由基頻光通過非線性效應產(chǎn)生的,因而在圖9(a)中對應了一個向下的脈沖形狀.最終腔內(nèi)形成的穩(wěn)定光場,基頻光和倍頻光功率都會在一定范圍內(nèi)等幅振蕩 (如圖9(b)中放大部分的曲線所示).當Ein=1000 W1/2和Ein=1200 W1/2時,腔內(nèi)功率變化曲線如圖9(c)—(f)所示.抽運功率增加,在圖9(d)和(f)中,倍頻光功率變化曲線中出現(xiàn)了兩個脈沖,并且抽運光功率越強,曲線中第二個脈沖的強度更大,因此對應的基頻光功率振蕩也更加明顯.從LLE的(1)式和(2)式可以看出,非線性效應的強弱和腔內(nèi)光功率的大小是相關(guān)的,腔內(nèi)光功率越大,由自相位調(diào)制和交叉相位調(diào)制引起的倍頻效率也越高.因此,當抽運功率較弱時,腔內(nèi)基頻和倍頻光強也處于相對較弱的水平,增益、損耗、色散、非線性之間相互作用可以達到相對平衡的狀態(tài),最終基頻和倍頻光場均會趨于穩(wěn)定.而當抽運功率較強時,腔內(nèi)增益和非線性均表現(xiàn)為較強的作用,此時腔內(nèi)光場難以達到穩(wěn)定狀態(tài).因此,光學微腔抽運功率過大,會導致腔內(nèi)光場處于不穩(wěn)定的狀態(tài).為保證腔內(nèi)光場和光功率的穩(wěn)定,需要避免抽運功率過強的情況.

    圖7 抽運功率對腔內(nèi)光功率的影響 (a)腔內(nèi)基頻光功率的變化情況;(b)腔內(nèi)倍頻光功率的變化情況Fig.7.Influence of the pump power on the power change curves:(a) Power variation for the fundamental frequency waves;(b) power variation for the second harmonic waves.

    圖8 Ein=100 W1/2 時,微腔內(nèi)光場的演化 (a) 基頻光功率的變化曲線;(b) 倍頻光功率的變化曲線;(c) 光場穩(wěn)定后,腔內(nèi)基頻光光場的周期性演化;(d) 光場穩(wěn)定后,腔內(nèi)倍頻光光場的周期性演化Fig.8.Evolution of the light field in the microresonator at Ein=100 W1/2:(a) Curve of power variation for the fundamental frequency waves;(b) curve of power variation for the second harmonic waves;(c) periodic evolution of the fundamental frequency waves;(d) periodic evolution of the second harmonic waves.

    圖9 基頻光和倍頻光的功率變化曲線 Ein=800 W1/2 時 (a)基頻光和 (b) 倍頻光的功率變化;Ein=1000 W1/2時 (c)基頻光和(d) 倍頻光的功率變化;Ein=1200 W1/2時 (e) 基頻光和 (f)倍頻光的功率變化Fig.9.Curves of the power variation for the fundamental frequency waves and the second harmonic waves:Power variation for(a) the fundamental frequency waves and (b) the second harmonic waves at Ein=800 W1/2;power variation for (c) the fundamental frequency waves and (d) the second harmonic waves at Ein=1000 W1/2;power variation for (e) the fundamental frequency waves and (f) the second harmonic waves at Ein=1200 W1/2.

    3.4 微腔大小對腔內(nèi)光場的影響

    以上都是在微腔 FSR=231.3 GHz (對應的微腔半徑為100 μm)的條件下進行的討論.但是研究中發(fā)現(xiàn),微腔半徑的大小同樣會對腔內(nèi)光場產(chǎn)生影響.設微腔 FSR=300 GHz,δ1=0.1,Ein=100 W1/2,得到腔內(nèi)基頻光和倍頻光功率的隨時間變換,最終都收斂于某一值,如圖10(a),(b)所示,即在圍墻內(nèi)可以得到穩(wěn)定的基頻、倍頻的光場分布.當光功率穩(wěn)定后,光場演化情況如圖10(c),(d)所示,由于微腔半徑增大,在腔內(nèi)基頻光與倍頻光之間非線性作用的距離變長,在較強的非線性作用下,基頻光和倍頻光最終在腔內(nèi)演化出16個脈沖,在對應的光譜圖中 (圖10(e),(f)),每隔 16 FSR,就會有一個功率較大的梳齒,此時微腔工作于“圖靈環(huán)”模式 (turing patterns)[17].

    研究過程中發(fā)現(xiàn),可工作于“圖靈環(huán)”模式的微腔FSR的取值范圍較小.將微腔的FSR減小為298 GHz,腔內(nèi)最終穩(wěn)定的基頻光場分布如圖11(a)所示,基頻光的光場在腔內(nèi)仍然為多脈沖的形式,但脈沖的強度受到了調(diào)制,這種調(diào)制在倍頻光的光場分布中更為明顯 (如圖11(b)所示),多脈沖的光場受到正弦信號的調(diào)制.由于光場中包含了多個頻率成分,因此其光譜圖中噪聲較多,在主梳齒附近會出現(xiàn)很多的毛刺.進一步減小微腔的FSR,當減小為 295 GHz 時,“圖靈環(huán)”工作模式消失,基頻光和倍頻光的光場均演化為正弦分布,且正弦信號的相位隨時間變化,光譜也轉(zhuǎn)化為很窄的梳狀光譜(如圖12所示).

    圖10 微腔 FSR=300 GHz 時,微腔內(nèi)出現(xiàn)“圖靈環(huán)”(Ein=100 W1/2,δ1=0.1) (a)基頻光功率變化曲線;(b)倍頻光功率變化曲線;(c) 腔內(nèi)光場穩(wěn)定后,基頻光光場隨時間的演化;(d) 腔內(nèi)光場穩(wěn)定后,倍頻光光場隨時間的演化;(e) 腔內(nèi)光場穩(wěn)定后,基頻光的光譜;(f) 腔內(nèi)光場穩(wěn)定后,倍頻光的光譜Fig.10.Turning patterns in the microresonator,when FSR=300 GHz (Ein=100 W1/2,δ1=0.1):(a) Curves of the power variation for the fundamental frequency waves;(b) curves of the power variation for the second harmonic waves;(c) evolution of the fundamental frequency waves;(d) evolution of the second harmonic waves;(e)spectra of the fundamental harmonic waves;(f) spectra of the second harmonic waves.

    圖11 微腔 FSR=298 GHz 時,穩(wěn)定后的光場分布及光譜 (a)基頻光的光場分布;(b) 倍頻光的光場分布;(c) 基頻光光譜;(d) 倍頻光光譜Fig.11.Stable intensity profile and spectra of the fundamental frequency waves for FSR=298 GHz:(a) Intensity profile of the fundamental frequency waves;(c) spectrum of the fundamental frequency waves;(b) intensity profile of second harmonic waves;(d) spectrum of second harmonic waves.

    圖12 微腔 FSR=295 GHz 時,穩(wěn)定后的基頻光和倍頻光光場分布及光譜 (a) 基頻光光場隨時間的演化;(b) 倍頻光光場隨時間的演化;(c) 基頻光場的瞬時分布;(d) 倍頻光場的瞬時分布;(e) 基頻光場的光譜;(f) 倍頻光場的光譜Fig.12.Stable intensity profile and spectra of the fundamental frequency waves for FSR=295 GH:(a) Evolution of the fundamental frequency waves;(b) evolution of the second frequency waves;(c) intensity profile of the fundamental frequency waves;(d) spectrum of the fundamental frequency waves;(e) intensity profile of the second harmonic waves;(f) spectrum of the second harmonic waves.

    4 結(jié) 論

    在考慮了光學微腔中二階和三階非線性效應的基礎上,引入了包含二次諧波項的LLE,分析了SiN微腔中二次諧波的產(chǎn)生,并討論了各參數(shù)對腔內(nèi)基頻和倍頻光場的影響.

    研究結(jié)果表明,當抽運頻率和微腔諧振頻率之間沒有發(fā)生頻率失諧,即失諧參量 δ1=0.1 時,腔內(nèi)的基頻光和倍頻光功率不斷增強,并最終達到穩(wěn)定,穩(wěn)定后的基頻光場為平頂脈沖的形式,而倍頻光場呈正弦分布.當抽運頻率和諧振頻率之間發(fā)生頻率失諧后,腔內(nèi)基頻和倍頻光功率在演化過程中出現(xiàn)振蕩,但最終穩(wěn)定于某一值.穩(wěn)定后的光場也發(fā)生變化,基頻光場變?yōu)榘倒伦拥男问?倍頻光場仍然為正弦分布.失諧參量δ1增加,將導致腔內(nèi)基頻和倍頻光功率在演化過程中振蕩變強,且最終穩(wěn)定的光功率變?nèi)?穩(wěn)定后的光場分布為周期性變化.δ1進一步增加會使得微腔光場處于混沌狀態(tài).因此,為保證微腔中二次諧波的產(chǎn)生以及腔內(nèi)光場的穩(wěn)定,實驗中應盡量減小抽運頻率與微腔振蕩頻率之間的頻率失諧.

    當微腔抽運功率較小時,增大抽運功率可有效增強腔內(nèi)基頻和倍頻光場的強度;而當抽運光強較強時,會導致腔內(nèi)色散以及非線性效應過強,最終穩(wěn)定的光場仍然呈周期性變化,且抽運光功率越強,光功率的演化曲線振蕩越強.因此,為保證腔內(nèi)光場和光功率的穩(wěn)定,需要避免抽運功率過強的情況.

    光學微腔的尺寸同樣會對腔內(nèi)光場分布產(chǎn)生影響.當選取微腔 FSR=300 GHz,失諧參量 δ1=0.1,抽運強度 Ein=100 W1/2,微腔工作于“圖靈環(huán)”狀態(tài),腔內(nèi)基頻和倍頻光場均演化出穩(wěn)定的16 個脈沖,相對應光譜每隔 16 FSR,就會有一個功率較大的梳齒.當FSR減小后,腔內(nèi)多脈沖形式的光場受到了正弦信號的調(diào)制,光譜中也出現(xiàn)較多的毛刺.因此,選擇合適的微腔尺寸,可以有效選擇微腔的工作模式,得到需要的光譜.該理論分析結(jié)果對于實驗中微腔參數(shù)和抽運參數(shù)的選取有重要意義.此外,本文中采用的微腔光場演化模型是二維模型,忽略了光場在垂直方向上(y方向)的分布.為了更精確地分析腔光場的分布和演化情況,需要考慮光場在垂直方向分布情況時,因而要引入三維的光場演化模型,此時抽運光場是關(guān)于y的表達式,腔內(nèi)基頻和倍頻光場也和y有關(guān).關(guān)于微腔內(nèi)基頻和倍頻光場三維模型的建立和分析,將在后續(xù)的工作中進一步展開.

    猜你喜歡
    微腔光場倍頻
    基于時域有限差分法的光子晶體波導微腔耦合研究
    利用新型光場顯微鏡高速記錄神經(jīng)元活動和血流動態(tài)變化
    科學(2020年5期)2020-01-05 07:03:12
    微腔結(jié)構(gòu)介質(zhì)阻擋放電的仿真與實驗研究
    壓縮混沌光場的量子統(tǒng)計性質(zhì)研究
    基于LBO晶體三倍頻的激光實驗系統(tǒng)的研究
    脈沖單頻Nd∶YVO4激光器及其倍頻輸出特性研究
    中國光學(2015年5期)2015-12-09 09:00:42
    集成光場三維顯示亮度均勻性校正方法
    瘋攝影
    攝影之友(2014年3期)2014-04-21 20:21:38
    Q開關(guān)倍頻Nd:YAG激光治療激素依賴性皮炎療效觀察
    粒子在二維開放型四分之一圓形微腔中的逃逸研究
    99热只有精品国产| 少妇的丰满在线观看| 亚洲av熟女| 精品免费久久久久久久清纯| 舔av片在线| 91麻豆精品激情在线观看国产| 久久精品国产亚洲av香蕉五月| 国产伦精品一区二区三区视频9 | 搡女人真爽免费视频火全软件 | 不卡一级毛片| 欧美一区二区精品小视频在线| 两个人的视频大全免费| 伊人久久大香线蕉亚洲五| 99国产极品粉嫩在线观看| 男女下面进入的视频免费午夜| 黄色女人牲交| 国产精品久久久久久精品电影| 一a级毛片在线观看| 一进一出好大好爽视频| 欧美午夜高清在线| 亚洲精品色激情综合| 成年女人永久免费观看视频| 亚洲 国产 在线| 两个人的视频大全免费| 97碰自拍视频| 少妇高潮的动态图| 久久久久久久久中文| 午夜福利成人在线免费观看| 内地一区二区视频在线| av黄色大香蕉| 国产精品野战在线观看| 亚洲国产日韩欧美精品在线观看 | 国产一区二区激情短视频| 亚洲性夜色夜夜综合| 91在线精品国自产拍蜜月 | 国产激情欧美一区二区| 小蜜桃在线观看免费完整版高清| 一级毛片高清免费大全| 亚洲最大成人中文| 国产精品1区2区在线观看.| 夜夜看夜夜爽夜夜摸| 91久久精品电影网| 女人高潮潮喷娇喘18禁视频| av欧美777| 麻豆国产97在线/欧美| 中国美女看黄片| 别揉我奶头~嗯~啊~动态视频| 中文字幕精品亚洲无线码一区| 在线观看免费午夜福利视频| 午夜激情福利司机影院| 99久久综合精品五月天人人| 精品久久久久久久久久久久久| 国产激情欧美一区二区| 国产探花极品一区二区| 麻豆成人av在线观看| 白带黄色成豆腐渣| 1000部很黄的大片| 村上凉子中文字幕在线| 国产精品,欧美在线| 亚洲av成人av| 舔av片在线| 国产精华一区二区三区| 午夜激情福利司机影院| 免费人成视频x8x8入口观看| 精品久久久久久久末码| 久久亚洲真实| 真人一进一出gif抽搐免费| 桃色一区二区三区在线观看| 中文字幕人妻熟人妻熟丝袜美 | 亚洲美女黄片视频| 九九久久精品国产亚洲av麻豆| 欧美黑人欧美精品刺激| 国产黄色小视频在线观看| av视频在线观看入口| 亚洲成av人片免费观看| 99热6这里只有精品| 精品日产1卡2卡| 99精品在免费线老司机午夜| 国产精品99久久久久久久久| 日本成人三级电影网站| 精品一区二区三区人妻视频| 亚洲真实伦在线观看| 3wmmmm亚洲av在线观看| 男女做爰动态图高潮gif福利片| 亚洲,欧美精品.| 深爱激情五月婷婷| 18美女黄网站色大片免费观看| 亚洲av免费高清在线观看| 国产成人av教育| 法律面前人人平等表现在哪些方面| 午夜精品一区二区三区免费看| 色视频www国产| 又爽又黄无遮挡网站| 亚洲人成网站高清观看| 亚洲国产欧美网| 老熟妇乱子伦视频在线观看| 免费在线观看影片大全网站| 亚洲精品成人久久久久久| 中文字幕人妻丝袜一区二区| 久久久精品欧美日韩精品| 成人午夜高清在线视频| 国产亚洲精品久久久com| 欧美丝袜亚洲另类 | 国产高清视频在线观看网站| 香蕉丝袜av| www.www免费av| xxxwww97欧美| 中文字幕人妻熟人妻熟丝袜美 | 亚洲内射少妇av| 男女视频在线观看网站免费| 在线观看免费午夜福利视频| 欧美另类亚洲清纯唯美| 可以在线观看毛片的网站| 少妇高潮的动态图| 精品人妻偷拍中文字幕| 一a级毛片在线观看| 欧美一区二区国产精品久久精品| 18禁美女被吸乳视频| 亚洲黑人精品在线| 国产真实乱freesex| 亚洲成a人片在线一区二区| 亚洲精品乱码久久久v下载方式 | 男女做爰动态图高潮gif福利片| 丰满人妻一区二区三区视频av | 俺也久久电影网| 日本撒尿小便嘘嘘汇集6| 精品一区二区三区av网在线观看| 高潮久久久久久久久久久不卡| 在线播放无遮挡| 国产精品99久久99久久久不卡| 国内精品美女久久久久久| 精品国产三级普通话版| 欧美性猛交黑人性爽| 国产精品久久视频播放| 91麻豆av在线| 99久久精品一区二区三区| 夜夜看夜夜爽夜夜摸| 精品国产三级普通话版| 亚洲性夜色夜夜综合| 国产蜜桃级精品一区二区三区| 欧美绝顶高潮抽搐喷水| 成人特级黄色片久久久久久久| 亚洲精品在线观看二区| 久久久成人免费电影| 久久久国产精品麻豆| 一区二区三区免费毛片| 亚洲成a人片在线一区二区| 亚洲在线观看片| 在线观看66精品国产| 亚洲美女黄片视频| 淫秽高清视频在线观看| 欧美xxxx黑人xx丫x性爽| 国产美女午夜福利| 日韩欧美国产在线观看| 一本综合久久免费| 一区二区三区激情视频| 99国产极品粉嫩在线观看| 国产91精品成人一区二区三区| 综合色av麻豆| 桃红色精品国产亚洲av| 欧美极品一区二区三区四区| 真人做人爱边吃奶动态| 国产老妇女一区| 日韩成人在线观看一区二区三区| АⅤ资源中文在线天堂| 欧美zozozo另类| 热99re8久久精品国产| 国产在视频线在精品| 国产精品99久久久久久久久| 亚洲 国产 在线| 日韩欧美三级三区| 亚洲国产精品成人综合色| 最新在线观看一区二区三区| 啦啦啦观看免费观看视频高清| 日韩精品中文字幕看吧| 国产男靠女视频免费网站| 九色成人免费人妻av| 色精品久久人妻99蜜桃| 蜜桃亚洲精品一区二区三区| 午夜日韩欧美国产| 亚洲一区二区三区色噜噜| 亚洲精品成人久久久久久| 九九热线精品视视频播放| 一本久久中文字幕| 熟女少妇亚洲综合色aaa.| 国产免费男女视频| 大型黄色视频在线免费观看| 88av欧美| 日韩av在线大香蕉| 国产真实乱freesex| 精品久久久久久久久久免费视频| 在线十欧美十亚洲十日本专区| 97碰自拍视频| 国产精品三级大全| 性色avwww在线观看| 精品日产1卡2卡| 伊人久久精品亚洲午夜| 成人国产一区最新在线观看| 亚洲精品色激情综合| 天天躁日日操中文字幕| 欧美一级毛片孕妇| netflix在线观看网站| 久久久久亚洲av毛片大全| 欧美黄色片欧美黄色片| 国模一区二区三区四区视频| 亚洲成a人片在线一区二区| netflix在线观看网站| 好男人电影高清在线观看| 国产成人福利小说| 青草久久国产| www国产在线视频色| 在线国产一区二区在线| 亚洲成人久久性| 宅男免费午夜| 日本 欧美在线| 亚洲男人的天堂狠狠| 欧美+日韩+精品| 伊人久久大香线蕉亚洲五| 欧美黄色片欧美黄色片| 18禁美女被吸乳视频| 免费在线观看成人毛片| 成人精品一区二区免费| 麻豆久久精品国产亚洲av| 哪里可以看免费的av片| 毛片女人毛片| 一本精品99久久精品77| 又紧又爽又黄一区二区| 国产真实伦视频高清在线观看 | 特大巨黑吊av在线直播| 精品久久久久久久人妻蜜臀av| 免费在线观看影片大全网站| 嫩草影视91久久| 婷婷六月久久综合丁香| 久久精品影院6| 久久久精品欧美日韩精品| 无人区码免费观看不卡| 国产成人av激情在线播放| 国产午夜福利久久久久久| 久久这里只有精品中国| 久久6这里有精品| 色综合欧美亚洲国产小说| 日韩国内少妇激情av| 村上凉子中文字幕在线| 国产一区二区激情短视频| 免费看十八禁软件| 成熟少妇高潮喷水视频| netflix在线观看网站| 国产精品影院久久| 久久久精品欧美日韩精品| 狂野欧美白嫩少妇大欣赏| 麻豆一二三区av精品| 1000部很黄的大片| 天堂动漫精品| 国产真实乱freesex| 国产伦精品一区二区三区四那| 精品一区二区三区视频在线观看免费| 国产成人系列免费观看| 男女之事视频高清在线观看| 两人在一起打扑克的视频| 久久久久免费精品人妻一区二区| av中文乱码字幕在线| 国产精品爽爽va在线观看网站| 国产中年淑女户外野战色| 亚洲欧美日韩高清在线视频| 在线天堂最新版资源| 久久精品国产自在天天线| 俄罗斯特黄特色一大片| 国产v大片淫在线免费观看| 又黄又爽又免费观看的视频| 久久久成人免费电影| 亚洲欧美日韩高清专用| 中文字幕久久专区| 亚洲电影在线观看av| 十八禁网站免费在线| 国产精品美女特级片免费视频播放器| 一区二区三区国产精品乱码| 久久精品人妻少妇| 中文字幕熟女人妻在线| 国产成人欧美在线观看| 少妇高潮的动态图| 国产精品,欧美在线| 久9热在线精品视频| 美女被艹到高潮喷水动态| av中文乱码字幕在线| 天美传媒精品一区二区| 亚洲av不卡在线观看| 中文字幕av在线有码专区| 国产精品久久久久久亚洲av鲁大| 超碰av人人做人人爽久久 | 一本综合久久免费| 3wmmmm亚洲av在线观看| 国产精品久久视频播放| 人妻丰满熟妇av一区二区三区| 国产免费一级a男人的天堂| 精品久久久久久久毛片微露脸| 成年女人毛片免费观看观看9| 国产毛片a区久久久久| 国产私拍福利视频在线观看| 亚洲成av人片免费观看| 脱女人内裤的视频| av天堂在线播放| 亚洲七黄色美女视频| 国产精品久久电影中文字幕| 亚洲人成网站在线播放欧美日韩| 久久精品影院6| 亚洲国产色片| 91九色精品人成在线观看| 少妇人妻一区二区三区视频| 亚洲五月天丁香| 国内精品久久久久精免费| 精品人妻偷拍中文字幕| 制服丝袜大香蕉在线| 亚洲不卡免费看| 日韩欧美 国产精品| 综合色av麻豆| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品合色在线| 国产精品香港三级国产av潘金莲| 中亚洲国语对白在线视频| av专区在线播放| 日韩成人在线观看一区二区三区| 丝袜美腿在线中文| 淫妇啪啪啪对白视频| av中文乱码字幕在线| 女人十人毛片免费观看3o分钟| 国产 一区 欧美 日韩| 成人av一区二区三区在线看| 天天添夜夜摸| 亚洲国产中文字幕在线视频| 亚洲中文字幕一区二区三区有码在线看| 中文字幕熟女人妻在线| 99热这里只有精品一区| 99久久久亚洲精品蜜臀av| 天美传媒精品一区二区| 日韩欧美国产在线观看| 国产精品综合久久久久久久免费| 免费大片18禁| 日本在线视频免费播放| 午夜福利在线观看免费完整高清在 | av在线天堂中文字幕| 成人亚洲精品av一区二区| 身体一侧抽搐| 听说在线观看完整版免费高清| 欧美另类亚洲清纯唯美| 动漫黄色视频在线观看| 嫩草影院精品99| 精品99又大又爽又粗少妇毛片 | www国产在线视频色| 99视频精品全部免费 在线| 噜噜噜噜噜久久久久久91| 欧美另类亚洲清纯唯美| 成年女人看的毛片在线观看| 欧美性猛交黑人性爽| 欧美中文综合在线视频| 1000部很黄的大片| 日本a在线网址| 中文在线观看免费www的网站| 亚洲精品在线观看二区| 日韩欧美免费精品| 午夜久久久久精精品| 91字幕亚洲| 久久天躁狠狠躁夜夜2o2o| 午夜免费男女啪啪视频观看 | or卡值多少钱| 老司机午夜福利在线观看视频| 全区人妻精品视频| 国产精品女同一区二区软件 | 国产三级黄色录像| 欧美极品一区二区三区四区| 亚洲国产高清在线一区二区三| 欧美一区二区国产精品久久精品| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成人久久性| 国产精品久久久久久人妻精品电影| 免费观看精品视频网站| 国产精品自产拍在线观看55亚洲| 桃色一区二区三区在线观看| 日本黄色视频三级网站网址| 亚洲性夜色夜夜综合| 制服人妻中文乱码| 亚洲乱码一区二区免费版| 高清毛片免费观看视频网站| 国产激情欧美一区二区| 国产伦精品一区二区三区视频9 | 免费在线观看日本一区| 在线视频色国产色| 桃红色精品国产亚洲av| 韩国av一区二区三区四区| 亚洲av免费在线观看| 日本撒尿小便嘘嘘汇集6| 欧美区成人在线视频| 亚洲av电影在线进入| 日韩欧美国产在线观看| 欧美日本亚洲视频在线播放| 国产亚洲精品久久久com| 中文资源天堂在线| 精品乱码久久久久久99久播| 免费看日本二区| 悠悠久久av| 91久久精品国产一区二区成人 | 久久久久久人人人人人| 女人十人毛片免费观看3o分钟| 中文字幕久久专区| av欧美777| 热99re8久久精品国产| 欧美最新免费一区二区三区 | 一a级毛片在线观看| 国产成人av教育| 全区人妻精品视频| 我要搜黄色片| 男女那种视频在线观看| 国产99白浆流出| 最近最新中文字幕大全免费视频| 夜夜躁狠狠躁天天躁| 国产午夜精品久久久久久一区二区三区 | 国产乱人视频| 日韩av在线大香蕉| 免费一级毛片在线播放高清视频| 村上凉子中文字幕在线| 欧美成狂野欧美在线观看| 99国产精品一区二区三区| 19禁男女啪啪无遮挡网站| svipshipincom国产片| 99国产极品粉嫩在线观看| 亚洲av免费高清在线观看| 亚洲国产精品合色在线| 亚洲av日韩精品久久久久久密| 高清在线国产一区| 日本与韩国留学比较| 一本一本综合久久| 网址你懂的国产日韩在线| a级毛片a级免费在线| 一本精品99久久精品77| 亚洲国产精品sss在线观看| 久久九九热精品免费| 99国产精品一区二区三区| a级毛片a级免费在线| 亚洲久久久久久中文字幕| 久久精品人妻少妇| 桃红色精品国产亚洲av| 日本熟妇午夜| 国产精品爽爽va在线观看网站| 国产色婷婷99| 亚洲片人在线观看| 中文字幕精品亚洲无线码一区| 国产精品1区2区在线观看.| 午夜免费观看网址| 欧美av亚洲av综合av国产av| 午夜精品在线福利| 亚洲一区二区三区不卡视频| 九色成人免费人妻av| 国产 一区 欧美 日韩| 国产精品av视频在线免费观看| 99热6这里只有精品| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品久久国产高清桃花| 熟妇人妻久久中文字幕3abv| 国产精品三级大全| 国产高清三级在线| 色综合亚洲欧美另类图片| 欧美一级毛片孕妇| 亚洲乱码一区二区免费版| 狂野欧美白嫩少妇大欣赏| 一区二区三区激情视频| 日韩欧美精品免费久久 | 欧美激情久久久久久爽电影| 听说在线观看完整版免费高清| 国产激情偷乱视频一区二区| 国产精品免费一区二区三区在线| 欧美中文日本在线观看视频| 日韩免费av在线播放| 日韩亚洲欧美综合| 精品国产亚洲在线| 日韩av在线大香蕉| 亚洲最大成人中文| 舔av片在线| 波多野结衣巨乳人妻| 久久亚洲精品不卡| 别揉我奶头~嗯~啊~动态视频| 在线播放无遮挡| 最好的美女福利视频网| 欧美性感艳星| 老司机午夜福利在线观看视频| 久久精品国产综合久久久| 国产乱人视频| a级毛片a级免费在线| 国产成人aa在线观看| 亚洲av成人精品一区久久| 在线天堂最新版资源| 国产精品日韩av在线免费观看| 国内久久婷婷六月综合欲色啪| 国产精品一区二区免费欧美| 伊人久久大香线蕉亚洲五| 免费观看人在逋| 久久伊人香网站| 国内精品美女久久久久久| 叶爱在线成人免费视频播放| 久久精品国产综合久久久| 日本 av在线| 欧美日韩一级在线毛片| 欧美日韩国产亚洲二区| 热99在线观看视频| 日韩亚洲欧美综合| 无限看片的www在线观看| 一进一出好大好爽视频| 亚洲最大成人手机在线| 亚洲午夜理论影院| 啦啦啦免费观看视频1| 日韩欧美精品免费久久 | 青草久久国产| 国产亚洲精品av在线| 91av网一区二区| 亚洲精品久久国产高清桃花| 国产精品野战在线观看| 国产成人系列免费观看| 桃红色精品国产亚洲av| 日韩欧美在线乱码| 精品国产三级普通话版| 国产亚洲欧美98| 亚洲 国产 在线| 亚洲av免费高清在线观看| 国产老妇女一区| 免费看十八禁软件| 亚洲精品一卡2卡三卡4卡5卡| 欧美在线一区亚洲| 欧美日韩瑟瑟在线播放| 久久久久久久精品吃奶| 99久久九九国产精品国产免费| 啪啪无遮挡十八禁网站| 国产精品国产高清国产av| 老鸭窝网址在线观看| 久久国产精品影院| 午夜老司机福利剧场| 欧美日韩综合久久久久久 | 9191精品国产免费久久| 琪琪午夜伦伦电影理论片6080| av专区在线播放| 亚洲国产欧美人成| 色综合婷婷激情| 亚洲av一区综合| 亚洲国产精品合色在线| 免费搜索国产男女视频| 9191精品国产免费久久| 欧美一区二区亚洲| 国内精品久久久久精免费| 久久久久久国产a免费观看| 色综合婷婷激情| 免费看美女性在线毛片视频| 日本熟妇午夜| 一个人观看的视频www高清免费观看| 国产探花极品一区二区| 乱人视频在线观看| 脱女人内裤的视频| 久久精品影院6| 中文字幕高清在线视频| 九九热线精品视视频播放| 午夜日韩欧美国产| 18禁在线播放成人免费| 久久精品国产亚洲av涩爱 | 久久国产乱子伦精品免费另类| 免费在线观看亚洲国产| 人人妻人人澡欧美一区二区| 国产成人欧美在线观看| av国产免费在线观看| 国产精品久久电影中文字幕| 国产精品久久久人人做人人爽| 精品久久久久久久久久免费视频| 亚洲熟妇中文字幕五十中出| 看片在线看免费视频| 偷拍熟女少妇极品色| 成人国产综合亚洲| 亚洲五月天丁香| 久久精品91蜜桃| 精品熟女少妇八av免费久了| 国产极品精品免费视频能看的| 一级作爱视频免费观看| 又粗又爽又猛毛片免费看| 99久久精品国产亚洲精品| 国产成人系列免费观看| 国产成人啪精品午夜网站| 高清日韩中文字幕在线| 国产探花在线观看一区二区| 亚洲内射少妇av| 国产熟女xx| 两个人视频免费观看高清| 免费av不卡在线播放| 一本一本综合久久| 嫩草影视91久久| 久久久国产成人精品二区| 精华霜和精华液先用哪个| 97超级碰碰碰精品色视频在线观看| 在线天堂最新版资源| tocl精华| 午夜福利高清视频| xxxwww97欧美| 亚洲成人中文字幕在线播放| svipshipincom国产片| 日韩 欧美 亚洲 中文字幕| 蜜桃久久精品国产亚洲av| 国产精品一及| 免费观看精品视频网站| 国产精品爽爽va在线观看网站| 久久精品亚洲精品国产色婷小说| 亚洲国产精品成人综合色| 亚洲 欧美 日韩 在线 免费| 欧美色视频一区免费| 久久久久亚洲av毛片大全| 亚洲第一欧美日韩一区二区三区| 热99在线观看视频| 精品久久久久久久久久免费视频| 国产一区二区亚洲精品在线观看| 国产成+人综合+亚洲专区| 亚洲精品影视一区二区三区av| 精品国产超薄肉色丝袜足j| bbb黄色大片|