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

    LIGO的O1/O2引力波事件對Einstein-aether理論限制能力探討

    2021-03-29 15:08:24劉剛強曹周鍵
    天文學(xué)進(jìn)展 2021年1期
    關(guān)鍵詞:雙星引力波引力

    劉剛強,曹周鍵

    (北京師范大學(xué) 天文學(xué)系,100875)

    1 引 言

    廣義相對論是繼牛頓萬有引力理論后,在強引力場、強動態(tài)時空區(qū)域更精準(zhǔn)的引力理論。廣義相對論預(yù)言的很多物理現(xiàn)象,包括星光偏折、引力紅移、夏皮羅時間延遲、參考系被拖拽的效應(yīng)和黑洞等都一一得到驗證。在廣義相對論發(fā)表100年后,2015年9月15日引力波被美國的激光干涉儀引力波天文臺(the laser interferometer gravitational-wave observatory,LIGO)成功直接探測到。引力波的成功探測完成了廣義相對論最后一個理論預(yù)言的實驗檢驗[1,2],也為人們觀測天體和宇宙提供了一種全新手段[3,4]。

    然而,廣義相對論仍然存在各種各樣的問題。暗物質(zhì)和暗能量始終是廣義相對論這片廣袤天空上的兩朵烏云。對于推動宇宙加速膨脹起到重要作用的暗能量[5],我們只能考慮有效場論,唯象層面上在標(biāo)準(zhǔn)廣義相對論下加入一些修正(通常在標(biāo)準(zhǔn)廣義相對論下額外考慮一個標(biāo)量場),最典型的是k-essence理論[6]。對于暗物質(zhì)同樣如此,1933年,F(xiàn)irtz Zwicy發(fā)現(xiàn)后發(fā)座星系團中各個星系團的速度彌散太高,最開始認(rèn)為只是一些不發(fā)光的物質(zhì)。20世紀(jì)80年代左右,大量的觀測表明,這種不發(fā)光的物質(zhì)并非我們之前所認(rèn)識到的基本粒子,這對基礎(chǔ)物理提出了一個很大的挑戰(zhàn),這種物質(zhì)雖然具有很強的宏觀引力效應(yīng),但是卻在局域上與物質(zhì)相互作用很弱[7]。另外,廣義相對論與量子理論在邏輯體系上不相容[8];而且廣義相對論自身預(yù)言時空奇點的存在,這些都暗示著廣義相對論的理論缺陷。

    對廣義相對論所有理論預(yù)言的成功檢驗既表明了廣義相對論的成功,也給人們帶來機會尋找廣義相對論失效的物理條件。牛頓萬有引力理論適用范圍對應(yīng)的弱場和低速,都是相對于廣義相對論黑洞致密程度對應(yīng)的強場和相對論預(yù)言的極限光速而言。可以想象,廣義相對論失效的物理條件是相對于一種新理論給出來的特征物理量而言,所以尋找廣義相對論失效物理條件的最好方法就是對比不同于廣義相對論的引力理論,用實驗數(shù)據(jù)來進(jìn)行理論選擇[9]。引力波實驗正好對應(yīng)強引力場、強動態(tài)時空區(qū)域的物理條件。引力波實驗正在以前所未有的標(biāo)準(zhǔn)檢驗著廣義相對論,同時眾多的其他引力理論被引力波實驗否定[5,10–14]。毫無疑問,科學(xué)研究進(jìn)展使人們對引力基本物理有了更深刻的理解。

    本文以Einstein-aether理論為例,探討引力波對超越廣義相對論引力理論的實驗檢驗。Einstein-aether理論的引力作用量可寫為:

    其中,Uμ是以太對應(yīng)的四速,gμν是時空度規(guī),R是里奇標(biāo)量是Kronecker delta記號[15]。GEA是Einstein-aether理論的有效引力常數(shù),有別于牛頓的萬有引力常數(shù),c1,c2,c3,c4是Einstein-aether理論的4個參數(shù),如果這4個參數(shù)都等于0,則Einsteinaether理論回到廣義相對論。Einstein-aether理論的參數(shù)組合為[16,17]:

    脈沖雙星的觀測[16]和宇宙學(xué)的觀測[17]將其限制為:

    文章余下部分的安排如下:第2章介紹筆者提出的Einstein-aether理論中圓軌道雙星系統(tǒng)的后愛因斯坦引力波波形形式;第3章介紹本文使用的Fisher Matrix參數(shù)估計方法;第4章給出引力波對Einstein-aether理論中參數(shù)的限制結(jié)果和檢驗Einstein-aether理論的結(jié)果;第5章是小結(jié)和討論。

    2 Einstein-aether理論中圓軌道雙星系統(tǒng)的后愛因斯坦引力波波形

    根據(jù)文獻(xiàn)[18]的計算結(jié)果,結(jié)合圓軌道雙星系統(tǒng)的TaylorF2引力波波形模型[19,20],我們把Einstein-aether理論中圓軌道雙星系統(tǒng)的后愛因斯坦引力波波形寫成如下形式[21]:

    其中,γ是歐拉常數(shù),M=m1+m2是雙星總質(zhì)量,是雙星系統(tǒng)的對稱質(zhì)量比,flso是最后穩(wěn)定軌道的頻率,f是引力波的頻率,p0,p1,p2,p3,p4,p5,p6,p7和l0,l1,l2,l3,l4,l5,l6,l7是引力波相位Ψ1和Ψ2的展開系數(shù)。本文提出的波形模型與文獻(xiàn)[21]略有不同。文獻(xiàn)[21]關(guān)注后牛頓近似,波形公式(包括式(7),(9),(20)和(27))中的質(zhì)量都是啁啾質(zhì)量(chirp mass)。我們關(guān)注波形的準(zhǔn)確性,使用了總質(zhì)量。我們的波形模型在α=β?1=β0=0時,嚴(yán)格回到LIGO數(shù)據(jù)處理使用的TaylorF2波形模型。在唯象層次上,我們的上述模型很容易直接推廣到LIGO數(shù)據(jù)分析軟件庫中最先進(jìn)的IMRPhenom波形模型[22,23]。文獻(xiàn)[24]中的公式(2)(即文獻(xiàn)[25]的公式(2.5))里的ψn和ψ(l)n,以及文獻(xiàn)[2]中的ψ?2和ψ0與本文采用的參數(shù)對應(yīng)關(guān)系如下:

    考慮到GW170817的多信使天文學(xué)對引力波速度的精確測定(精度約為10?15)[5,10–12,26],我們直接把引力波速度設(shè)定為光速。于是將上述波形表達(dá)式與文獻(xiàn)[18]的計算結(jié)果相結(jié)合,我們得到:

    其中,ψ是引力波偏振角,F(xiàn)b,FL,FX,FY是引力波探測器對相應(yīng)引力波偏振模式的響應(yīng)函數(shù),γb,γL,γX1,γX2,γY1,γY2是Einstein-aether理論中相應(yīng)引力波偏振模式的強度系數(shù)。相應(yīng)函數(shù)和偏振模式強度系數(shù)可以表達(dá)為Einstein-aether理論參數(shù)c1,c2,c3,c4的函數(shù),相關(guān)表達(dá)式的信息參考文獻(xiàn)[18]。?x和κ3是Einstein-aether理論給出的參數(shù),也是c1,c2,c3,c4的函數(shù),分別定義在文獻(xiàn)[18]中的公式(4.6)和(6.7)。廣義相對論對應(yīng)?x=0和κ3=1。本文接下來介紹如何利用引力波探測來限制α,β?1和β0。

    高級LIGO探測器自2015年9月開機以來,已經(jīng)歷了三次觀測運行,其探測靈敏度還在不斷提高。高級LIGO探測器的噪聲功率譜可以近似成下面的形式[25]:

    其中,f為引力波頻率,x為約化頻率,S0為總的噪聲水平。

    3 Fisher Matrix的參數(shù)估計方法

    引力波探測器所測得的數(shù)據(jù)d包含信號s和噪聲n,即d(t)=s(t)+n(t)。在理想情況下,探測器的噪聲是高斯穩(wěn)態(tài)噪聲,即d(t)?s(t)是一個隨機過程,而且各種聯(lián)合概率分布都是與時間無關(guān)的高斯分布。

    基于引力波探測器的靈敏度,我們定義兩組時間序列h1(t)和h2(t)的內(nèi)積如下:

    這里h1和h2頭上的“?”代表傅立葉變換,“?”代表取復(fù)共軛,“?”代表取實部,(fmin,fmax)對應(yīng)探測器的敏感頻率范圍。對于高斯穩(wěn)態(tài)噪聲n(t),其內(nèi)積應(yīng)該滿足高斯分布:

    于是我們就有條件概率分布:

    對于引力波探測,我們關(guān)心的是基于探測數(shù)據(jù)d(t),獲得信號h(t)的概率是多大,也就是條件概率p(h|d)。根據(jù)貝葉斯定理,我們可得:

    上式分母實為分子對應(yīng)的概率歸一化因子。在貝葉斯統(tǒng)計學(xué)里,p(h|d)被稱為后驗概率,p(h)被稱為先驗概率,p(d|h)被稱為似然函數(shù)。

    對于某一種類引力波源,如雙星系統(tǒng),引力波波形除了是時間的函數(shù)外,還依賴于波源的參數(shù),如雙黑洞的質(zhì)量、自旋和軌道離心率等。我們把這些參數(shù)記作,則可得h()。概率分布p(h|d)和p(h)是相對于的函數(shù),即是參數(shù)空間上的函數(shù)?!襭(d|h)p(h)是對整個參數(shù)空間積分。

    后驗概率p(h|d)作為參數(shù)空間上的函數(shù),其最大值所對應(yīng)的點被稱為最可幾參數(shù)0。在引力波的匹配濾波數(shù)據(jù)處理中,0對應(yīng)最好匹配的那組參數(shù)[1]。我們再進(jìn)一步假設(shè)信號很強,即信噪比很高,噪聲n基本可忽略不計,則有d≈h(0),?n,n?≈?n,h?≈?n,?λih?≈?n,?λi?λjh?≈0。如果我們簡化認(rèn)為先驗均勻分布p(h)∝1,則根據(jù)式(36)和(37)可得:

    其中,C是概率分布的歸一化因子。我們考慮點處的泰勒展開近似:

    上式右邊的矩陣??λih,?λjh?剛好對應(yīng)統(tǒng)計學(xué)里的Fisher Matrix。也就是說在高信噪比條件下,后驗概率分布在最可幾參數(shù)點可近似為協(xié)方差矩陣Σij,等于Fisher Matrix逆的高斯分布。本文接下來使用Fisher Matrix的方法來估計引力波測量的參數(shù)誤差

    本文使用的Einstein-aether理論中,圓軌道雙星系統(tǒng)的后愛因斯坦引力波波形模型包含參數(shù)(A,M,η,tc,?c,α,β?1,β0),它們分別是引力波幅度(與引力波源的光度距離和方位角等有關(guān))、啁啾質(zhì)量、對稱質(zhì)量比、并合時間與相位,以及Einstein-aether理論對應(yīng)的后愛因斯坦參數(shù)。我們在計算中發(fā)現(xiàn)A與其他參數(shù)的關(guān)聯(lián)度很弱,該參數(shù)是否參與Fisher Matrix,所得其他參數(shù)的誤差基本相同。使用高級LIGO探測器的靈敏度,對典型雙星系統(tǒng)參數(shù)誤差結(jié)果如表1所示。我們的結(jié)果與文獻(xiàn)[21]的一致。由于我們的波形模型與其的細(xì)微差別導(dǎo)致了定量上的細(xì)微差別,從差別可以看出,我們采用更精確的波形模型(等效到3.5階后牛頓)求出的參數(shù)誤差比采用文獻(xiàn)[21]的波形模型(等效到2階后牛頓)算得的略小。

    表1 使用Fisher Matrix方法針對高級LIGO靈敏度對典型雙星系統(tǒng)進(jìn)行參數(shù)估計的結(jié)果

    關(guān)于tc和?c這兩個參數(shù),我們還嘗試了使用F統(tǒng)計方法來處理[19]。對于參數(shù)誤差,使用F統(tǒng)計與否,Einstein-aether理論對應(yīng)的后愛因斯坦參數(shù)限制結(jié)果基本沒有差別。在下文計算中,我們得到的結(jié)果都是沒有使用F統(tǒng)計的計算結(jié)果。

    4 對Einstein-aether理論進(jìn)行限制

    高級LIGO主導(dǎo)的LIGO科學(xué)合作組織第一次觀測運行從2015年9月12日到2016年1月19日,共發(fā)現(xiàn)GW150914,GW151012和GW151226三個雙黑洞并合事件。高級LIGO主導(dǎo)的LIGO科學(xué)合作組織第二次觀測運行從2016年11月30日開始到2017年8月25日結(jié)束。其中高級VIRGO從2017年8月1日開始加入第二次觀測運行。LIGO科學(xué)合作組織第二次觀測運行共發(fā)現(xiàn)GW170104,GW170608,GW170729,GW170809,GW170814,GW170817,GW170818和GW170823八個雙星并合事件。

    我們把LIGO科學(xué)合作組織第一次和第二次觀測運行發(fā)現(xiàn)的11個雙星并合事件對Einstein-aether理論進(jìn)行參數(shù)估計的結(jié)果列在了表2中。因為只有GstLAL數(shù)據(jù)處理流水線發(fā)現(xiàn)了所有的11個事件,所以表2里所列信噪比SNR和雙星系統(tǒng)的質(zhì)量m1和m2都是由GstLAL數(shù)據(jù)處理流水線所得[27]。我們將LIGO-VIRGO第一次和第二次探測到的11個引力波事件的參數(shù)代入到式(5)—(8)中,得出帶有Einstein-aether理論的3個參數(shù)的模擬波形,再使用Fisher Matrix方法估計引力波事件對這3個參數(shù)的限制能力。從該表結(jié)果可以看出,雙星系統(tǒng)的質(zhì)量越小,則對Einstein-aether理論的參數(shù)限制越嚴(yán)格;這是因為質(zhì)量越小,落到高級LIGO和高級VIRGO探測器頻段中的引力波信號越多,對Einstein-aether理論的驗證就越有效。當(dāng)雙星系統(tǒng)的質(zhì)量增大,直到并合頻率落到探測器頻段中,雖然并合信號會較強,但我們所使用的引力波波形模型無法利用上并合短信號,導(dǎo)致所得限制不強。

    表2 11個雙星并合事件對Einstein-aether理論進(jìn)行參數(shù)估計的結(jié)果

    綜合LIGO科學(xué)合作組織第一次和第二次觀測運行發(fā)現(xiàn)的11個雙星并合事件對Einstein-aether理論進(jìn)行參數(shù)限制的結(jié)果,可以得到:

    但我們要注意,以上限制結(jié)果基于假定廣義相對論的引力波模板可以完全匹配探測器所測得的信號。引力波對α參數(shù)限制不好,對β?1和β0限制不錯;這是因為α與引力波的幅度相關(guān),β?1和β0與引力波的相位相關(guān)。而匹配濾波本來就只是敏感依賴于相位。

    為了估計廣義相對論波形模板匹配Einstein-aether理論引力波信號的程度,我們引入引力波匹配因子FF[28,29]。給定兩個引力波波形h1(t)和h2(t),它們間的匹配因子為:

    當(dāng)用一個波形近似另一個波形信號時,信噪比會被降低FF倍。如果降得太多,這樣的波形近似就會失效。

    我們以GW170817的雙星參數(shù)為例,研究廣義相對論與Einstein-aether理論的引力波波形匹配情況。鑒于匹配濾波與引力波幅度不是敏感依賴,我們只討論匹配因子與β?1和β0的關(guān)系,結(jié)果見圖1。從圖1我們可以看出,只要β?1和β0偏離0分別在10?6和10?3的量級,Einstein-aether理論的引力波波形與廣義相對論的波形完全不匹配。所以只要GW170817的引力波信號是廣義相對論描述的波形,Einstein-aether理論的β?1和β0參數(shù)就可以得到很好的限制。

    圖1 類似GW170817的廣義相對論引力波波形和Einstein-aether理論的引力波波形匹配情況

    這里還存在一個問題:是否GW170817實為Einstein-aether理論描述,只是因為β?1和β0與啁啾質(zhì)量的參數(shù)簡并使得廣義相對論模板沒有遺漏這個信號?接下來我們研究這個問題。為此,我們設(shè)定β?1=6×10?5,β0=0和β?1=0,β0=0.002,然后改變啁啾質(zhì)量,研究波形匹配因子,結(jié)果見圖2。從圖2可以看出,的確可以通過改變啁啾質(zhì)量來使得廣義相對論的引力波波形與Einstein-aether理論的引力波波形相匹配,即β?1和β0與啁啾質(zhì)量存在參數(shù)簡并。所以在不確定所探測到的引力波信號可由廣義相對論描述的情況下,引力波的探測不能被認(rèn)已經(jīng)對Einstein-aether理論給出了很好的限制,現(xiàn)在我們只能說廣義相對論能很好地解釋所有探測到的引力波信號[2]。

    圖2 Einstein-aether理論中β?1和β0與啁啾質(zhì)量的參數(shù)簡并情況

    5 結(jié)論和討論

    自2015年LIGO首次直接探測到引力波以來,高級LIGO和高級VIRGO一起進(jìn)行了三次觀測,共探測到幾十例引力波事件,包括雙黑洞、雙中子星和黑洞-中子雙星系統(tǒng)的并合,以引力波為觀測手段的引力波天文學(xué)被很好地建立起來。

    引力波天文學(xué)可以用來研究基本物理學(xué)中的引力相互作用[22,24,30]。GW170817就被成功應(yīng)用到引力理論的檢驗中,若干不同于廣義相對論的引力理論得到了很強的限制。

    在文獻(xiàn)[18]的計算結(jié)果基礎(chǔ)上,結(jié)合TaylorF2引力波波形的模式,我們提出了Einstein-aether理論中雙星系統(tǒng)的引力波波形模型?;谒岢龅牟ㄐ文P?,我們利用Fisher Matrix方法,結(jié)合LIGO科學(xué)合作組織公開的第一、第二次觀測運行結(jié)果數(shù)據(jù)中的11個雙星系統(tǒng),對Einstein-aether理論進(jìn)行了限制、研究。Fisher Matrix方法只對參數(shù)微小偏離,或者說只在Einstein-aether理論微小修正廣義相對論引力波波形的情形下有效。這種情形下,使用Einstein-aether理論的引力波波形模板對引力波實際數(shù)據(jù)的處理會得出與廣義相對論相似的結(jié)果,所以我們只需要使用Fisher Matrix方法對后驗概率分布做微擾研究即可。這樣的條件使得我們可以把LIGO已經(jīng)測出的11個雙星對應(yīng)的參數(shù)輸入到Fisher Matrix方法進(jìn)行分析。反之,一旦參數(shù)微小偏離不成立,F(xiàn)isher Matrix方法將失效。我們采用匹配濾波中匹配因子的方法,探討廣義相對論引力波波形與Einstein-aether理論的引力波波形存在顯著偏離的可能性及其對引力波數(shù)據(jù)處理的影響。

    結(jié)果表明,如果Einstein-aether理論只是廣義相對論的微小偏離(在引力波波形意義上),則我們可以把Einstein-aether理論中的參數(shù)限制為α<0.166,β?1<8.07×10?6,β0<7.79×10?2。但只是基于引力波探測,我們還不能完全排除Einstein-aether理論顯著偏離廣義相對論的可能性,這時Einstein-aether理論的參數(shù)和雙星系統(tǒng)的啁啾質(zhì)量存在簡并。我們計劃在后續(xù)工作中,使用貝葉斯模型選擇的方法來研究這種簡并的可能性有多大。

    猜你喜歡
    雙星引力波引力
    雙星啟示錄
    李雙星 一心為民拔“窮根”
    河北畫報(2020年21期)2020-12-14 03:16:08
    黃浦江邊的“引力波”
    EN菌的引力波探測器
    引力
    初中生(2017年3期)2017-02-21 09:17:40
    發(fā)現(xiàn)引力波
    感受引力
    新春“引力波”一觸即發(fā)
    海峽姐妹(2016年2期)2016-02-27 15:15:18
    A dew drop
    雙星計劃中的雙星
    太空探索(2015年12期)2015-07-18 10:59:19
    中文字幕熟女人妻在线| 久久精品国产综合久久久| 成人欧美大片| 嫩草影院入口| 久久久久久久久久黄片| 少妇人妻一区二区三区视频| 免费人成视频x8x8入口观看| av福利片在线观看| 亚洲 欧美 日韩 在线 免费| 国产高潮美女av| 午夜福利在线观看免费完整高清在 | 亚洲性夜色夜夜综合| 香蕉丝袜av| 嫩草影院入口| av片东京热男人的天堂| 成人18禁在线播放| 日日摸夜夜添夜夜添小说| 亚洲无线观看免费| av福利片在线观看| 婷婷精品国产亚洲av在线| 美女免费视频网站| 亚洲内射少妇av| 深爱激情五月婷婷| 国产精品久久视频播放| 国产黄色小视频在线观看| 国产成人影院久久av| 99热6这里只有精品| 熟女少妇亚洲综合色aaa.| 久久久久久久亚洲中文字幕 | 久久久色成人| 亚洲人成网站在线播| 午夜久久久久精精品| 一区二区三区国产精品乱码| 在线观看66精品国产| 日韩人妻高清精品专区| 久久精品国产亚洲av香蕉五月| 久久久久久久午夜电影| 国产一区二区激情短视频| 18禁在线播放成人免费| 日韩高清综合在线| 精品无人区乱码1区二区| 成熟少妇高潮喷水视频| 成人亚洲精品av一区二区| 黑人欧美特级aaaaaa片| 波野结衣二区三区在线 | 宅男免费午夜| 欧美性感艳星| 男女床上黄色一级片免费看| 欧美成人性av电影在线观看| 青草久久国产| 两人在一起打扑克的视频| 国产成人aa在线观看| 老司机午夜福利在线观看视频| av黄色大香蕉| 亚洲午夜理论影院| 亚洲色图av天堂| 一个人看的www免费观看视频| 在线天堂最新版资源| 国产三级中文精品| 两个人的视频大全免费| 日韩免费av在线播放| 亚洲精品日韩av片在线观看 | 少妇熟女aⅴ在线视频| 欧美区成人在线视频| 日韩精品青青久久久久久| 成人国产一区最新在线观看| 欧美乱码精品一区二区三区| 99久久精品热视频| 国产伦精品一区二区三区视频9 | 一区福利在线观看| 欧美性猛交╳xxx乱大交人| 俺也久久电影网| 成人av在线播放网站| 精品99又大又爽又粗少妇毛片 | 日本 欧美在线| 色av中文字幕| 熟妇人妻久久中文字幕3abv| 此物有八面人人有两片| 99热只有精品国产| 啪啪无遮挡十八禁网站| 国产日本99.免费观看| 日日摸夜夜添夜夜添小说| 1000部很黄的大片| 宅男免费午夜| 国产精品久久久久久人妻精品电影| 老司机午夜十八禁免费视频| 夜夜躁狠狠躁天天躁| 国产精品1区2区在线观看.| 美女免费视频网站| 极品教师在线免费播放| 欧美黑人巨大hd| 国产日本99.免费观看| 亚洲av成人精品一区久久| 一本综合久久免费| 中文亚洲av片在线观看爽| 精品电影一区二区在线| 香蕉av资源在线| 久久精品夜夜夜夜夜久久蜜豆| 色综合亚洲欧美另类图片| 精品人妻1区二区| 无遮挡黄片免费观看| 欧美在线黄色| 十八禁网站免费在线| 两个人视频免费观看高清| tocl精华| 色尼玛亚洲综合影院| 久久精品国产综合久久久| 精品电影一区二区在线| 亚洲国产欧洲综合997久久,| 女人十人毛片免费观看3o分钟| 97超级碰碰碰精品色视频在线观看| 国产成人a区在线观看| 99国产精品一区二区三区| 亚洲精品亚洲一区二区| 变态另类丝袜制服| 国产成人aa在线观看| 日日夜夜操网爽| 99国产精品一区二区三区| 亚洲精品一区av在线观看| 国产成年人精品一区二区| 69av精品久久久久久| 99久久精品国产亚洲精品| 日本撒尿小便嘘嘘汇集6| tocl精华| 观看免费一级毛片| 亚洲av中文字字幕乱码综合| 久久精品91蜜桃| 日韩成人在线观看一区二区三区| 琪琪午夜伦伦电影理论片6080| 九九热线精品视视频播放| 欧美日韩中文字幕国产精品一区二区三区| 欧美性感艳星| 18禁裸乳无遮挡免费网站照片| 久9热在线精品视频| 一区二区三区免费毛片| 国产美女午夜福利| av片东京热男人的天堂| 久久久国产成人精品二区| 国产高清videossex| 国内精品久久久久精免费| netflix在线观看网站| 国产精品乱码一区二三区的特点| 亚洲 国产 在线| 欧美av亚洲av综合av国产av| 长腿黑丝高跟| 一个人免费在线观看的高清视频| 午夜福利免费观看在线| 久久香蕉国产精品| 日韩高清综合在线| 国产精品 欧美亚洲| 午夜视频国产福利| av天堂中文字幕网| 久久久久久人人人人人| 国产主播在线观看一区二区| 丰满乱子伦码专区| 色精品久久人妻99蜜桃| 久久精品国产99精品国产亚洲性色| 国产激情欧美一区二区| 激情在线观看视频在线高清| 欧美又色又爽又黄视频| 99在线人妻在线中文字幕| 久久性视频一级片| 国产精品久久久久久久久免 | 男女午夜视频在线观看| 999久久久精品免费观看国产| 久久精品91蜜桃| 成人特级黄色片久久久久久久| 精品人妻偷拍中文字幕| 亚洲国产精品sss在线观看| www.色视频.com| 丝袜美腿在线中文| 国产69精品久久久久777片| 免费观看精品视频网站| 欧美最新免费一区二区三区 | 亚洲专区国产一区二区| 婷婷六月久久综合丁香| 国产精品精品国产色婷婷| 少妇人妻精品综合一区二区 | 国产精品一及| 他把我摸到了高潮在线观看| 免费大片18禁| 精品一区二区三区视频在线 | 婷婷精品国产亚洲av| 欧美日韩精品网址| 老熟妇仑乱视频hdxx| 久久这里只有精品中国| 18禁黄网站禁片免费观看直播| 又黄又粗又硬又大视频| 国产中年淑女户外野战色| 日本黄大片高清| 午夜激情福利司机影院| h日本视频在线播放| 久久天躁狠狠躁夜夜2o2o| 男人舔奶头视频| 成人亚洲精品av一区二区| 亚洲七黄色美女视频| 深爱激情五月婷婷| 久久欧美精品欧美久久欧美| 两个人看的免费小视频| 变态另类丝袜制服| 欧美日韩国产亚洲二区| 成年人黄色毛片网站| 国产精品一区二区三区四区久久| 麻豆成人av在线观看| 成年免费大片在线观看| 高清日韩中文字幕在线| 欧美日本视频| 俄罗斯特黄特色一大片| 手机成人av网站| 欧美3d第一页| 国产成人福利小说| 国产午夜福利久久久久久| 国产精品久久久久久久久免 | bbb黄色大片| 97超视频在线观看视频| 国产乱人视频| 国产午夜精品论理片| 波野结衣二区三区在线 | 精品乱码久久久久久99久播| 亚洲精品国产精品久久久不卡| 久久99热这里只有精品18| 97人妻精品一区二区三区麻豆| 亚洲欧美日韩无卡精品| 国产精品香港三级国产av潘金莲| 一个人免费在线观看电影| 精品一区二区三区视频在线 | 婷婷精品国产亚洲av在线| 五月伊人婷婷丁香| 乱人视频在线观看| 男人和女人高潮做爰伦理| 亚洲avbb在线观看| 欧美+日韩+精品| 欧美区成人在线视频| 人人妻人人澡欧美一区二区| 亚洲av美国av| 69av精品久久久久久| 精品欧美国产一区二区三| 在线观看日韩欧美| av片东京热男人的天堂| 窝窝影院91人妻| 国产精品98久久久久久宅男小说| 亚洲国产精品合色在线| 亚洲成人免费电影在线观看| 我的老师免费观看完整版| 久久人人精品亚洲av| 久久久久久国产a免费观看| 久久久精品欧美日韩精品| 午夜福利视频1000在线观看| 精品人妻一区二区三区麻豆 | 国产成人aa在线观看| 精品久久久久久成人av| 午夜免费观看性视频| 看免费成人av毛片| 日韩av在线免费看完整版不卡| 女人十人毛片免费观看3o分钟| 天天躁日日操中文字幕| 久久久久精品性色| 成年免费大片在线观看| 国产综合懂色| 80岁老熟妇乱子伦牲交| 久久久久久久久久成人| 亚洲成色77777| 精品一区二区免费观看| 久久精品久久精品一区二区三区| 成年免费大片在线观看| 国产午夜精品论理片| 国产av在哪里看| 国产av国产精品国产| 成人美女网站在线观看视频| 午夜福利网站1000一区二区三区| 日韩欧美三级三区| 99久久精品国产国产毛片| 一级爰片在线观看| 美女脱内裤让男人舔精品视频| 色综合亚洲欧美另类图片| 一区二区三区免费毛片| 国产乱人偷精品视频| 白带黄色成豆腐渣| 搡女人真爽免费视频火全软件| 国产高清有码在线观看视频| 免费黄网站久久成人精品| 日本-黄色视频高清免费观看| 97在线视频观看| 国产精品一二三区在线看| www.av在线官网国产| 亚洲美女视频黄频| 丰满少妇做爰视频| 国产91av在线免费观看| 亚洲欧美成人精品一区二区| 精品国产露脸久久av麻豆 | 亚洲精品色激情综合| 能在线免费看毛片的网站| 久久久久久久国产电影| 天天一区二区日本电影三级| 国产av在哪里看| 国产精品爽爽va在线观看网站| 午夜亚洲福利在线播放| 亚洲成人久久爱视频| 亚洲av中文字字幕乱码综合| 午夜福利在线在线| 亚洲国产欧美在线一区| 午夜福利视频1000在线观看| 精品人妻一区二区三区麻豆| 国产伦精品一区二区三区四那| 十八禁网站网址无遮挡 | 波多野结衣巨乳人妻| 亚洲经典国产精华液单| 七月丁香在线播放| 青青草视频在线视频观看| 高清欧美精品videossex| 精品久久久噜噜| 日韩av在线免费看完整版不卡| 在线观看美女被高潮喷水网站| 欧美xxxx性猛交bbbb| 联通29元200g的流量卡| 国产有黄有色有爽视频| 亚洲av福利一区| 精品欧美国产一区二区三| 日韩电影二区| 自拍偷自拍亚洲精品老妇| 国产黄色小视频在线观看| 免费观看精品视频网站| 日日撸夜夜添| videos熟女内射| 少妇丰满av| 韩国高清视频一区二区三区| av在线蜜桃| 日本一二三区视频观看| 精品国产露脸久久av麻豆 | 亚洲av中文字字幕乱码综合| 又粗又硬又长又爽又黄的视频| 亚洲久久久久久中文字幕| 91aial.com中文字幕在线观看| 午夜激情福利司机影院| 婷婷色麻豆天堂久久| 人妻少妇偷人精品九色| 午夜福利网站1000一区二区三区| 2021少妇久久久久久久久久久| 美女脱内裤让男人舔精品视频| 国产成人a∨麻豆精品| 九草在线视频观看| 超碰av人人做人人爽久久| 国产在线男女| 又粗又硬又长又爽又黄的视频| 男女那种视频在线观看| 99久久中文字幕三级久久日本| 亚洲国产成人一精品久久久| 又大又黄又爽视频免费| 狠狠精品人妻久久久久久综合| 美女大奶头视频| 免费黄色在线免费观看| ponron亚洲| 精品一区二区三区视频在线| 一区二区三区高清视频在线| 麻豆成人午夜福利视频| 国产免费福利视频在线观看| 国产片特级美女逼逼视频| 高清午夜精品一区二区三区| 性色avwww在线观看| 老司机影院毛片| 大片免费播放器 马上看| 熟妇人妻久久中文字幕3abv| 超碰av人人做人人爽久久| 久久亚洲国产成人精品v| 久久久久久国产a免费观看| 成人美女网站在线观看视频| 日韩三级伦理在线观看| 国产综合懂色| 国产男人的电影天堂91| 精品国内亚洲2022精品成人| 2018国产大陆天天弄谢| 波多野结衣巨乳人妻| 成人特级av手机在线观看| 80岁老熟妇乱子伦牲交| 青春草亚洲视频在线观看| 欧美日本视频| 乱人视频在线观看| 久久久久免费精品人妻一区二区| 99九九线精品视频在线观看视频| 午夜免费男女啪啪视频观看| 搡老乐熟女国产| 在线播放无遮挡| 精品久久久久久久久av| 国产在线男女| 噜噜噜噜噜久久久久久91| 成人午夜精彩视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 久久久久久久国产电影| 亚洲怡红院男人天堂| 亚洲精品一区蜜桃| 国产一区亚洲一区在线观看| 十八禁国产超污无遮挡网站| 日韩av免费高清视频| 亚洲av中文字字幕乱码综合| av在线天堂中文字幕| 久久亚洲国产成人精品v| 99九九线精品视频在线观看视频| 伦精品一区二区三区| 亚洲精品久久午夜乱码| 国产精品1区2区在线观看.| 18禁在线无遮挡免费观看视频| 黄色日韩在线| 久久久久九九精品影院| 熟女电影av网| 男人和女人高潮做爰伦理| 中文天堂在线官网| 非洲黑人性xxxx精品又粗又长| 欧美日韩亚洲高清精品| 99热网站在线观看| 国产探花在线观看一区二区| 2018国产大陆天天弄谢| 午夜免费观看性视频| 久久久久久国产a免费观看| 国产又色又爽无遮挡免| 日本免费a在线| 99re6热这里在线精品视频| 两个人视频免费观看高清| 国产精品.久久久| 91精品伊人久久大香线蕉| 欧美不卡视频在线免费观看| 插阴视频在线观看视频| 免费观看av网站的网址| 在线免费十八禁| 欧美精品一区二区大全| 高清欧美精品videossex| 欧美xxⅹ黑人| 能在线免费观看的黄片| 91aial.com中文字幕在线观看| 欧美3d第一页| 亚洲美女视频黄频| 久久久久久久亚洲中文字幕| 色综合站精品国产| 亚洲在久久综合| 男女视频在线观看网站免费| 亚洲成色77777| 春色校园在线视频观看| 91午夜精品亚洲一区二区三区| 国产中年淑女户外野战色| 亚洲自拍偷在线| 色视频www国产| 久久国产乱子免费精品| 国产乱人偷精品视频| xxx大片免费视频| 久久久a久久爽久久v久久| 久久久国产一区二区| 久久久久久久久中文| 欧美一区二区亚洲| 国产成人freesex在线| 免费观看性生交大片5| .国产精品久久| 国语对白做爰xxxⅹ性视频网站| 亚洲av日韩在线播放| 久久久久网色| 亚洲av国产av综合av卡| 亚洲av不卡在线观看| 九九久久精品国产亚洲av麻豆| 六月丁香七月| 亚洲,欧美,日韩| 国产高清有码在线观看视频| 久久久久免费精品人妻一区二区| 久久久久久国产a免费观看| 日韩视频在线欧美| 看非洲黑人一级黄片| 日韩一区二区三区影片| 我的老师免费观看完整版| 少妇猛男粗大的猛烈进出视频 | 蜜桃亚洲精品一区二区三区| 国产免费又黄又爽又色| 尾随美女入室| av网站免费在线观看视频 | 一级黄片播放器| av在线老鸭窝| 欧美xxⅹ黑人| 91av网一区二区| 亚洲av成人av| 99热网站在线观看| 一区二区三区四区激情视频| 校园人妻丝袜中文字幕| 国产在线男女| av一本久久久久| 成人综合一区亚洲| 国产精品爽爽va在线观看网站| 国产成人freesex在线| 能在线免费看毛片的网站| 91久久精品电影网| 亚洲精品456在线播放app| 美女国产视频在线观看| 人人妻人人看人人澡| 成人午夜高清在线视频| 国产成人精品久久久久久| 狂野欧美白嫩少妇大欣赏| 亚洲精品成人久久久久久| 狠狠精品人妻久久久久久综合| 黑人高潮一二区| 久久久精品免费免费高清| 男女边吃奶边做爰视频| 九草在线视频观看| 中文乱码字字幕精品一区二区三区 | 精品一区在线观看国产| 欧美不卡视频在线免费观看| 一本一本综合久久| 日产精品乱码卡一卡2卡三| 亚洲av电影不卡..在线观看| 亚洲av免费在线观看| 人人妻人人看人人澡| 久久97久久精品| 寂寞人妻少妇视频99o| 麻豆乱淫一区二区| 欧美bdsm另类| 免费人成在线观看视频色| 最近中文字幕高清免费大全6| 国产男女超爽视频在线观看| 最近的中文字幕免费完整| 日本午夜av视频| 国产亚洲av片在线观看秒播厂 | 天堂网av新在线| 国产高清国产精品国产三级 | 97超碰精品成人国产| 性插视频无遮挡在线免费观看| 欧美日韩综合久久久久久| 狂野欧美白嫩少妇大欣赏| av国产免费在线观看| 国产色婷婷99| 91在线精品国自产拍蜜月| 日本-黄色视频高清免费观看| 建设人人有责人人尽责人人享有的 | 国产日韩欧美在线精品| 少妇裸体淫交视频免费看高清| 日本免费在线观看一区| 啦啦啦啦在线视频资源| 一级毛片 在线播放| 亚洲精品一二三| 精品久久久精品久久久| 2021天堂中文幕一二区在线观| 午夜福利成人在线免费观看| 高清av免费在线| 成人高潮视频无遮挡免费网站| 好男人在线观看高清免费视频| 蜜臀久久99精品久久宅男| 日日啪夜夜爽| 深爱激情五月婷婷| 成年女人看的毛片在线观看| 91精品国产九色| 七月丁香在线播放| 91狼人影院| 一级毛片我不卡| 在线观看av片永久免费下载| av线在线观看网站| 青春草亚洲视频在线观看| 国产av码专区亚洲av| 成人鲁丝片一二三区免费| 亚州av有码| 亚洲最大成人中文| 简卡轻食公司| 婷婷六月久久综合丁香| 99热这里只有是精品50| 国产黄片视频在线免费观看| 国产av码专区亚洲av| 麻豆av噜噜一区二区三区| 国产精品人妻久久久影院| 免费看光身美女| 国产人妻一区二区三区在| 哪个播放器可以免费观看大片| 日韩电影二区| 色播亚洲综合网| 男人舔奶头视频| 一级av片app| 一级片'在线观看视频| 国产精品.久久久| 亚洲欧美一区二区三区国产| 99热网站在线观看| 特级一级黄色大片| 成年人午夜在线观看视频 | 婷婷色av中文字幕| 国产成年人精品一区二区| 亚洲欧美日韩东京热| 久99久视频精品免费| 日韩精品有码人妻一区| 80岁老熟妇乱子伦牲交| 欧美日韩国产mv在线观看视频 | 日韩欧美 国产精品| 免费黄网站久久成人精品| 激情 狠狠 欧美| 免费观看在线日韩| 中文字幕亚洲精品专区| 一级毛片电影观看| 三级男女做爰猛烈吃奶摸视频| 啦啦啦啦在线视频资源| 国产av不卡久久| 夫妻性生交免费视频一级片| 天堂俺去俺来也www色官网 | 欧美成人午夜免费资源| 久久久精品欧美日韩精品| 亚洲欧美一区二区三区黑人 | 超碰av人人做人人爽久久| 亚洲欧美一区二区三区黑人 | 亚洲欧洲日产国产| 成人亚洲欧美一区二区av| 建设人人有责人人尽责人人享有的 | 国产亚洲精品av在线| 男女视频在线观看网站免费| 国产亚洲5aaaaa淫片| 男女那种视频在线观看| 亚洲国产精品成人综合色| 国产亚洲5aaaaa淫片| 成人美女网站在线观看视频| 毛片一级片免费看久久久久| 亚洲精品自拍成人| av在线老鸭窝| 在线观看免费高清a一片| 国产精品福利在线免费观看| 大又大粗又爽又黄少妇毛片口| 亚洲aⅴ乱码一区二区在线播放| 婷婷色麻豆天堂久久| 日韩视频在线欧美|