999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

兩個獨立全光纖多通道光子糾纏源的Hong-Ou-Mandel干涉?

2017-08-07 07:59:24李銀海許昭懷王雙許立新周志遠史保森
物理學報 2017年12期
關鍵詞:測量

李銀海 許昭懷 王雙 許立新?周志遠? 史保森

1)(中國科學技術大學光學與光學工程系,合肥 230026)

2)(中國科學技術大學,中國科學院量子信息重點實驗室,合肥 230026)

3)(中國科學技術大學,量子信息和量子物理協作創新中心,合肥 230026)

兩個獨立全光纖多通道光子糾纏源的Hong-Ou-Mandel干涉?

李銀海1)許昭懷1)王雙2)3)許立新1)?周志遠2)3)?史保森2)3)

1)(中國科學技術大學光學與光學工程系,合肥 230026)

2)(中國科學技術大學,中國科學院量子信息重點實驗室,合肥 230026)

3)(中國科學技術大學,量子信息和量子物理協作創新中心,合肥 230026)

(2017年2月21日收到;2017年3月29日收到修改稿)

獨立光子源的干涉是實現復雜量子體系應用(比如多光子糾纏態產生和量子隱形傳態等)的核心技術.利用100 GHz密集波分復用技術,實現了1.55μm全光纖多通道獨立糾纏光子源的Hong-Ou-Mandel干涉,在不去除暗符合(隨機符合計數)的情況下,可見度為53.2%±8.4%,去除暗符合可見度可達到82.9%±5.3%.給出了關于色散位移光纖中基于自發四波混頻過程產生的單光子光譜純度嚴格的理論描述,模擬了抽運脈沖寬度和濾波器帶寬對單光子光譜純度的影響,并給出了理論上的最佳條件(最佳的抽運脈沖寬度為8 ps,高斯濾波器帶寬為40 GHz及以下).在測量Hong-Ou-Mandel干涉之前,先測量了液氮冷卻狀態下的色散位移光纖關聯光子源的符合和隨機符合比率,在抽運功率為23μW的情況下,最大比率可以達到131.Hong-Ou-Mandel干涉在高精度光學測量、測量裝置無關的量子密鑰分配等應用中扮演著極為重要的角色.

密集波分復用,Hong-Ou-Mandel干涉,色散位移光纖

1 引 言

獨立光子源之間的Hong-Ou-Mandel干涉(HOMI)是實現高精度光學測量[1]、復雜的多光子GHz態和團簇態產生、量子密鑰分配和光子NOON態測量[2?7]的核心技術,而這些技術正是量子科學的基石,因此獨立源之間的HOMI已被廣泛研究.最早在可見光波段的研究,是利用非線性晶體中的自發參量下轉換過程[8?13]或者是光子晶體光纖中的自發四波混頻過程展開[14?16];后來對通信波段的研究是利用非線性晶體中的自發參量下轉換過程展開[17].從這些工作中,我們了解到提高干涉可見度的關鍵是提高光子源的單光子光譜純度(簡記為純度),可采用的手段為管理色散、控制相位匹配條件和盡量減小濾波器帶寬.

當今基于波導的平臺(比如光纖和硅基波導)相較于傳統的空間光學系統表現出了巨大的潛力:低成本、高集成和光纖通信系統無縫連接等.許多科研工作者將精力集中到基于色散位移光纖(dispersion shifted fi ber,DSF)[18?23]和硅波導[24?26]的糾纏光子源產生和量子態管理中.為了在DSF中產生更復雜的量子態,獨立源的HOMI勢在必行.

對于DSF中的HOMI,Takesue[27]利用100 MHz重復頻率,脈沖寬度小于1 ps的鎖模激光器(波長1551.1 nm)和25 GHz的濾波器實現了53%的干涉可見度.在本文的工作中,首先詳細給出了DSF中關聯光子光譜純度的理論描述,并給出了理論上的最佳實驗條件;然后利用基于密集波分復用(DWDM)技術在DSF中產生的多通道關聯光子源[28],理論上實現了多通道的HOMI.

2 基于DSF中自發四波混頻過程的關聯光子光譜純度理論

2.1 理論模型

首先我們給出了基于DSF中自發四波混頻產生光子對的理論描述.產生光子對的量子態可以表述為[9,12,29]

每個模式都有一個由施密特數gn[29]決定的權重.施密特分解告訴我們,如果一個光子位于態則預示了對應的其余光子的態為對應的純度定義為而且純度由聯合譜振幅(joint spectral amplitude,JSA)f(ωs,ωi)[30]的分解因式決定. 雖然施密特分解不能解析地實現,但是我們可以數值模擬并計算奇異值.純度由奇異值的平方求和給出,與施密特量[31]相等.因此,一個必要條件是在單個光譜模式下,我們將f(ωs,ωi)分解為f(ωs,ωi)=fs(ωs)fi(ωi).

在DSF中產生的光子對的JSA決定于抽運包絡ε(ωs+ ωi)和光纖相位匹配條件Γ(ωs,ωi),因此f(ωs,ωi)= ε(ωs+ ωi)Γ(ωs,ωi). 假設抽運激光是高斯模式,帶寬為σp,則|ε(ωs+ ωi)|2=exp{?2[(ωs+ ωi? ωp)/σp]2}, 相較于非線性晶體中的自發參量下轉換過程,有一個因數2的不同,這是由于自發四波混頻產生信號光和閑頻光的過程不同[32].DSF中自發四波混頻的相位匹配條件可以寫成|Γ(ωs,ωi)|2=[sinc(?kL/2)]2,其中?k=2kp?ks?ki+2γP是相位失配,kp,ks,ki分別為抽運光、信號光和閑頻光的波數,γ=2.6 km?1·W?1為DSF的三階非線性系數,P為抽運功率.我們用高斯函數來逼近這個位相匹配條件|Γ(ωs,ωi)|2∝ exp(?2α2?k2L2), 其中α=0.220.在抽運光、信號光和閑頻光的中心頻率對相位失配展開,聯合譜強度(joint spectral intensity,JSI)可以寫成

2.2 模擬結果

結合實驗條件,我們數字模擬了抽運光的寬度和濾波器帶寬對JSI和光子源純度的影響,結果參見圖1.圖1(a)—(d)為相位匹配函數、抽運包絡函數、JSI不加40 GHz濾波器和加40 GHz濾波器的結果,圖中橫坐標和縱坐標為信號光和閑頻光分別相對其中心頻率的頻率偏移量,圖形表示了強度分布和頻率失諧量的關系.模擬的參數為其中Td=25 ps為抽運激光的脈沖寬度;L=300 m為光纖長度,理論上,不加濾波器時單光子純度為0.5693,加上100 GHz DWDM(等效高斯濾波帶寬約為40 GHz)濾波之后為0.8765.圖1(e)為不加濾波器的情況下抽運脈沖寬度對純度的影響,我們得到最佳的抽運脈沖寬度為8 ps.圖1(f)為高斯濾波對純度的影響,對應濾波器帶寬的降低,光子純度增加(抽運脈沖寬度固定為25 ps,和我們實際使用的激光器的脈沖寬度一致).

圖1 (網刊彩色)(a)相位匹配函數數學模擬;(b)抽運函數;(c)不加濾波器的JSI;(d)加40 GHz高斯濾波器的JSI;(e)光子純度和抽運脈沖寬度的關系;(f)光子純度和濾波器帶寬的關系Fig.1.(color online)(a)Numerical simulations of phase matching function;(b)pump functions;(c)JSI without spectral fi ltering;(d)JSI fi ltered with 40 GHz Gaussian fi lters;(e)single photon spectral purity against pump pulse width;(f)single photon spectral purity against fi lter bandwidth of the photon pairs.

3 實驗系統及測量結果

實驗裝置如圖2所示.我們用一個自制的脈寬25 ps、重復頻率27.9 MHz的鎖模激光器作為抽運光,然后用可調衰減器(TA)來調節光子源的抽運功率,利用級聯的100 GHz的DWDM濾波器來過濾抽運光.然后用50:50分束器將濾波后的抽運光分成兩路分別抽運兩段液氮冷卻的300 m長DSF光纖來產生兩個獨立的光子源.經過DSF之后,使用兩組級聯的200 GHz的DWDM濾波器以及光纖偏振旋轉器和偏振片過濾掉抽運光及其正交偏振方向的拉曼散射.然后兩個32通道的100 GHz DWDM器件被用來將對應通道的關聯光子對分開(選用信道C31和C37,對應中心波長1552.52 nm和1547.72 nm).兩個獨立源之間的延遲用可調光纖延遲線來調整.從這兩個32通道的DWDM器件相同信道出射的光子被引入同一個50:50耦合器進行干涉,另外兩路光子被用作觸發信號,用兩個連續觸發InGaAs單光子探測器(APD1,APD4,ID Quanta,ID220,20%探測效率,3μs死時間,暗計數率4k cps),它們輸出的電信號用來觸發兩個門觸發模式的單光子探測器(APD2,APD3,Qasky,合肥,中國,100 MHz,20%探測效率,每個門暗計數4×10?5)來做光子符合測量.APD2和APD3探測的光子信號輸出到我們的符合計數設備(Pico quanta,TimeHarp 260,1.6 ns符合窗口)來進行四光子符合測量.

圖2 兩個基于300 m DSF中自發四波混頻的獨立光子源的HOMI實驗裝置(DSF,色散位移光纖;APD1-4,雪崩光電二極管單光子探測器;PFR&P,光纖偏振旋轉器和檢偏器;PC1(2),偏振控制器)Fig.2.Experiment setup for HOMI between two independent photon sources generated by spontaneous four-wave mixing in two 300 m DSFs(DSF,dispersion shifted fi ber;APD1-4,avalanched photon detector;FPR and P, fi ber polarization rotator and polarizer;PC1(2),polarization controller).

4 討 論

圖3 (a)CAR和抽運功率的關系;(b)符合時間為1000 s,四光子符合和兩個單光子源的相對延遲的關系Fig.3.(a)CAR as a function of pump power;(b)four-fold coincidences in 1000 s as a function of the relative delay between the two single photon sources.

在測量兩個獨立單光子源的HOMI之前,我們首先測量了單光子源的符合和隨機符合比率(coincidence to accidental coincidence ratio,CAR)關于抽運功率的關系,結果見圖3(a).對于相對較低的抽運功率,CAR受限于由單光子探測器暗計數導致的暗符合,隨著抽運功率的上升,CAR隨之升高;在抽運功率為23μW時達到最大,CAR為131.抽運功率繼續升高,多光子效應開始出現并限制CAR的提高,所以CAR開始下降.然后我們進行了HOMI測量過程,濾波后的抽運功率為0.25 mW,每段DSF的抽運功率是0.12 mW.采用相對較高的功率可以縮短測量時間,但是相對較大的多光子效應會導致較低的干涉可見度.不去除暗符合,可見度為53.2%±8.4%,去除暗符合之后,為82.9%±5.3%,去除暗符合之后可見度的值非常接近我們理論計算的結果0.8765,分束器的分束比誤差和兩段DSF的不完全一致導致了和理論值的微小差異.通過斷開耦合器的一臂進行測量,我們可以得到暗符合的值.單光子探測器較高的暗計數和多光子效應導致的暗符合嚴重降低了HOMI的可見度.最大可見度可以表示為

5 結 論

我們給出了兩個獨立的多通道全光纖光子源的HOMI的理論研究和實驗結果,詳細討論了在DSF中抽運脈沖寬度和濾波器帶寬對單光子純度的影響,并給出了理論上最佳的實驗條件.需要進一步說明的是:1)我們的光子源抽運脈沖寬度并沒有處在理論最佳的實驗位置;2)試驗中使用的單光子探測器探測效率相對較低,暗計數相對較高,這是影響實驗結果的重要因素,如果能提高探測效率,降低暗計數,四光子符合實驗結果將能大大提高;3)在實驗中我們用了相對較高的抽運功率(0.12 mW)來大幅縮短測量時間(符合測量時間為1000 s),而在抽運功率為23μW時,DSF單光子源的CAR最佳可以達到131;4)我們在實驗中采用了32通道DWDM的一對信道(C31-C37),根據之前實現的多通道全光纖糾纏光子源,我們有理由相信可以采用其他信道(比如C30-C38,C32-C36)來拓展在量子通信任務中的信道容量.

總之,兩個獨立光子源的HOMI是光子對在多光子態操作、量子隱形傳態和測量裝置無關量子密鑰分配等應用中的關鍵技術,我們實現了兩個獨立的多通道全光纖糾纏光子源的HOMI.基于光纖的糾纏光子系統在未來的量子信息技術中扮演著極為重要的角色.

感謝中國科學技術大學的李海鷗博士提供了液氮、陳巍博士提供了單光子探測器以及武漢工程大學的金銳博博士提供給予的建議和幫助.

[1]Nagata T,Okamoto R,O’Brien J L,Sasaki K,Takeuchi S 2007 Science 316 726

[2]Pan J W,Chen Z B,Lu C Y,Weinfurter H,Zeilinger A,Zukowski M 2012 Rev.Mod.Phys.84 777

[3]Pirandola S,Eisert J,Weedbrook C,Furusawa A,Braunstein S L 2015 Nat.Photon.9 641

[4]Wang C,Song X T,Yin Z Q,Wang S,Chen W,Zhang C M,Guo G C,Han Z F 2015 Phys.Rev.Lett.115 160502

[5]Afek I,Ambar O,Silberberg Y 2010 Science 328 879

[6]Tan Y G,Liu Q 2016 Chin.Phys.Lett.33 090303

[7]Ma H X,Bao W S,Li H W,Chou C 2016 Chin.Phys.B 25 080309

[8]KaltenbaekR,BlauensteinerB,ZukowskiM,Aspelmeyer M,Zeilinger A 2006 Phys.Rev.Lett.96 240502

[9]Mosley P J,Lundeen J S,Smith B J,Wasylczyk P,U’Ren A B,Silberhorn C,Walmsley I A 2008 Phys.Rev.Lett.100 133601

[10]Mosley P J,Lundeen J S,Smith B J,Walmsley I A 2008 New J.Phys.10 093011

[11]Tanida M,Okamoto R,Takeuchi S 2012 Opt.Express 20 15275

[12]Kim Y H,Grice W P 2005 Opt.Lett.30 908

[13]Wasilewski W,Wasylczyk P,Kolenderski P,Banaszek K,Radzewicz C 2006 Opt.Lett.31 1130

[14]Fulconis J,Alibart O,O’Brien J L,Wadsworth W J,Rarity J G 2007 Phys.Rev.Lett.99 120501

[15]Halder M,Fulconis J,Cemlyn B,Clark A,Xiong C,Wadsworth W J,Rarity J G 2009 Opt.Express 17 4670

[16]Soller C,Cohen O,Smith B J,Walmsley I A,Silberhorn C 2011 Phys.Rev.A 83 031806

[17]Jin R B,Wakui K,Shimizu R,Benichi H,Miki S,Yamashita T,Terai H,Wang Z,Fujiwara M,Sasaki M 2013 Phys.Rev.A 87 063801

[18]McMillan A R,Labonte L,Clark A S,Bell B,Alibart O,Martin A,Wadsworth W J,Tanzilli S,Rarity J G 2013 Sci.Rep.3 2032

[19]Li X,Voss P L,Sharping J E,Kumar P 2005 Phys.Rev.Lett.94 053601

[20]Takesue H,Inoue K 2004 Phys.Rev.A 70 031802

[21]Takesue H,Inoue K 2005 Phys.Rev.A 72 041804

[22]Wang S X,Kanter G S 2009 IEEE J.Selected Topics in Quantum Electronics 15 1733

[23]Yang L,Li X Y,Wang B S 2008 Acta Phys.Sin.57 4933(in Chinese)[楊磊,李小英,王寶善 2008物理學報57 4933]

[24]Silverstone J W,Bonneau D,Ohira K,Suzuki N,Yoshida H,Iizuka N,Ezaki M,Natarajan C M,Tanner M G,Had fi eld R H,Zwiller V,Marshall G D,Rarity J G,O’Brien J L,Thompson M G 2014 Nat.Photon.8 104

[25]Reimer C,Kues M,Caspani L,Wetzel B,Roztocki P,Clerici M,Jestin Y,Ferrera M,Peccianti M,Pasquazi A,Little B E,Chu S T,Moss D J,Morandotti R 2015 Nat.Comunn.6 8236

[26]Reimer C,Kues M,Roztocki P,Wetzel B,Grazioso F,Little B E,Chu S T,Johnston T,Bromberg Y,Caspani L,Moss D J,Morandotti R 2016 Science 351 1176

[27]Takesue H 2007 Appl.Phys.Lett.90 204101

[28]Li Y H,Zhou Z Y,Xu Z H,Xu L X,Shi B S,Guo G C 2016 Phys.Rev.A 94 043810

[29]U’Ren A B,Silberhorn C,Banaszek K,Walmsley I A,Erdmann R,Grice W P,Raymer M G 2005 Laser Phys.15 146

[30]Grice W P,U’Ren A B,Walmsley I A 2001 Phys.Rev.A 64 063815

[31]Law C K,Walmsley I A,Eberly J H 2000 Phys.Rev.Lett.84 5304

[32]Jin R B,Shimizu R,Wakui K,Benichi H,Sasaki M 2013 Opt.Express 21 10659

PACS:03.67.Hk,42.79.Sz,42.81.UvDOI:10.7498/aps.66.120302

Hong-Ou-Mandel interference between two independent all- fi ber multiplexed photon sources?

Li Yin-Hai1)Xu Zhao-Huai1)Wang Shuang2)3)Xu Li-Xin1)?Zhou Zhi-Yuan2)3)?Shi Bao-Sen2)3)

1)(Department of Optics and Optical Engineering,University of Science and Technology of China,Hefei 230026,China)
2)(Key Laboratory of Quantum Information,University of Science and Technology of China,Hefei 230026,China)
3)(Synergetic Innovation Center of Quantum Information and Quantum Physics,University of Science and Technology of China,Hefei 230026,China)

21 February 2017;revised manuscript

29 March 2017)

Interference between independent photon sources is the key technique to realize complex quantum systems for more sophisticated applications such as multi-photon entanglement generation and quantum teleportation.Here,we report Hong-Ou-Mandel interference(HOMI)between two independent 1.55μm all- fi ber photon pair sources over two 100 GHz dense wave division multiplexing(DWDM)channels,whose visibility reaches 53.2%±8.4%(82.9%±5.3%)without(with)back ground counts subtracted.In addition,we theoretically describe in detail the single photon spectral purity of the photon source generated in dispersion shifted fi ber(DSF),simulate the in fl uences of the pulse width and fi lter bandwidth on the purity,and obtain the optimized condition.The optimized pump pulse width is 8 ps and fi lter bandwidth is about 40 GHz or less.A home-made 1550.1 nm mode-locked fi ber laser source,whose pulse width and repetition rate are 25 ps and 27.9 MHz respectively,acts as a pump of photon source.A tunable attenuator is used to adjust the pump power of the photon source,and the broad band background fl uorescence photons are fi ltered out by cascade 100 GHz DWDM fi lters.The clean pump beam is divided into two equal parts by the 50:50 optical coupler to pump two 300 m DSFs(cooled by liquid nitrogen)to generate independent photon sources.Then the strong pump beam and noise photon from Raman scattering in orthogonal polarization are removed by 2 groups of 200 GHz DWDM fi lters and fi ber polarization rotator and polarizer.Then two 100 GHz DWDMs are used for separating photons at correlated channel pairs.The relative delay between the two independent photons is adjusted by tunable fi ber delay line.Photons from the same channels are combined in a second beam splitter for interference,and the other two photons are used as trigger signals.The two triggered photons are detected by two free running InGaAs avalanched single photon detectors(APD1,APD4,ID Quanta,ID220,20%detection efficiency,3μs dead time,dark count rate 4k cps),and the outputs of detectors APD1 and APD4 are used to trigger two single-photon detectors running in the gated mode(APD2,APD3,Qasky,Hefei,China,100 MHz,free gating single photon detectors,20%detection efficiency,dark count probability 4×10?5per gate)for twophoton coincidence measurement.Detection output signals from APD2 and APD3 are sent to our coincidence count device(Pico quanta,TimeHarp 260,1.6 ns coincidence window)for four-photon coincidence measurement.Before measuring the HOMI,we obtain a maximum-coincidence-to-accidental-coincidence ratio(CAR)of 131 by cooling the fi ber in liquid nitrogen when the pump power is 23μW.There are a few remarks we want to point out.Firstly,the photon sources are not operated at the optimized pump pulse width for pure single photon generation,but narrow band 100 GHz fi lters are used in the experiments to increase the purity of the sources.Secondly,single photon detectors used in our experiment have lower detection efficiency and much higher dark counts than nano-wire single photon detectors,if we have high-performance nano-wire single photon detector,experimental results will be greatly improved due to the four-fold coincidences and dark coincidences scaling quadruplicate with the detection efficiency and dark count probability of a single detector.Thirdly,we use relatively high pump power for each DSF(0.12 mW)to reduce measurement time for photon coincidence,which will lead to a very poor raw visibility certainly.Finally,though only a 100 GHz channel pair is used in our experiment,we can use other channels for multiplexing such interference processes to improve the channel capacity in future quantum communication tasks theoretically.Our study shows greatly promising integrated optical elements for future scalable quantum information processing.

dense wave division multiplexing,Hong-Ou-Mandel interference,dispersion shifted fi ber

10.7498/aps.66.120302

?國家自然科學基金(批準號:11174271,61275115,61435011,61525504)和中央高校基本科研業務費專項資金(批準號:WK2030380009)資助的課題.

?通信作者.E-mail:xulixin@ustc.edu.cn

?通信作者.E-mail:zyzhouphy@ustc.edu.cn

?2017中國物理學會Chinese Physical Society

http://wulixb.iphy.ac.cn

*Project supported by the National Natural Science Foundation of China(Grant Nos.11174271,61275115,61435011,61525504)and the Fundamental Research Funds for the Central Universities,China(Grant No.WK2030380009).

?Corresponding author.E-mail:xulixin@ustc.edu.cn

?Corresponding author.E-mail:zyzhouphy@ustc.edu.cn

猜你喜歡
測量
測量重量,測量長度……
把握四個“三” 測量變簡單
滑動摩擦力的測量和計算
滑動摩擦力的測量與計算
測量的樂趣
二十四節氣簡易測量
日出日落的觀察與測量
滑動摩擦力的測量與計算
測量
測量水的多少……
主站蜘蛛池模板: 97国产在线播放| 一本大道香蕉中文日本不卡高清二区 | aa级毛片毛片免费观看久| 少妇高潮惨叫久久久久久| 亚洲精品无码久久久久苍井空| 一边摸一边做爽的视频17国产| 伊人久久大香线蕉aⅴ色| 18禁高潮出水呻吟娇喘蜜芽| 韩国自拍偷自拍亚洲精品| 国产在线观看精品| 欧美视频在线第一页| 好紧好深好大乳无码中文字幕| 九色综合伊人久久富二代| A级毛片无码久久精品免费| 欧美一级高清视频在线播放| …亚洲 欧洲 另类 春色| 最新国产高清在线| 国产高颜值露脸在线观看| 亚欧美国产综合| 97视频在线精品国自产拍| 欧美成人午夜视频免看| 五月婷婷亚洲综合| 国产精品网址你懂的| 伊在人亞洲香蕉精品區| 美女毛片在线| 亚洲日本中文综合在线| 久久久精品国产SM调教网站| 国产十八禁在线观看免费| 少妇高潮惨叫久久久久久| 亚洲精品大秀视频| 制服丝袜一区二区三区在线| 亚洲天堂网在线视频| 亚洲精品va| 91精品免费高清在线| 亚洲国产精品一区二区第一页免| 欧美色丁香| 亚洲欧美另类色图| 中国美女**毛片录像在线| 一级香蕉人体视频| 欧美狠狠干| 亚洲最新网址| 中日韩欧亚无码视频| 国产成人永久免费视频| 亚洲综合片| 麻豆精选在线| 波多野结衣久久高清免费| 日韩中文字幕免费在线观看| 国产精品欧美亚洲韩国日本不卡| 在线观看无码a∨| 亚洲AV无码乱码在线观看代蜜桃| 日本在线国产| 2020久久国产综合精品swag| 内射人妻无套中出无码| 好紧好深好大乳无码中文字幕| 国产精品无码一区二区桃花视频| 成人无码区免费视频网站蜜臀| 欧美日韩国产综合视频在线观看| 亚洲另类国产欧美一区二区| 玖玖精品在线| 亚洲黄色激情网站| 国产成人盗摄精品| 无码高潮喷水在线观看| 成人国产一区二区三区| 亚洲成人高清在线观看| 亚国产欧美在线人成| 美女视频黄又黄又免费高清| 免费观看国产小粉嫩喷水| 国产日韩精品欧美一区灰| 国产精品美女网站| 国产成人精品第一区二区| 中字无码精油按摩中出视频| 日韩不卡高清视频| 激情网址在线观看| 中字无码精油按摩中出视频| 精品国产中文一级毛片在线看| 欧美日韩精品一区二区视频| 无码精油按摩潮喷在线播放| 欧美成人看片一区二区三区| a级高清毛片| 亚洲视频免| 制服丝袜国产精品| 亚洲人成影院午夜网站|