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

    脈沖星射電輻射束半徑與周期關(guān)系的研究?

    2014-11-29 05:11:29魏丁丁游霄鵬
    天文學(xué)報(bào) 2014年6期
    關(guān)鍵詞:實(shí)線脈沖星偏振

    魏丁丁 游霄鵬

    (西南大學(xué)物理科學(xué)與技術(shù)學(xué)院重慶400715)

    脈沖星射電輻射束半徑與周期關(guān)系的研究?

    魏丁丁 游霄鵬?

    (西南大學(xué)物理科學(xué)與技術(shù)學(xué)院重慶400715)

    脈沖星的周期非常穩(wěn)定.前人通過研究發(fā)現(xiàn)脈沖星射電輻射束半徑大小應(yīng)該正比于周期的?1/3或?1/2次方.通過整理收集了87顆信噪比很好的、偏振位置角明顯為“S”曲線的脈沖星數(shù)據(jù),用脈沖星輻射束幾何模型計(jì)算得到了每顆脈沖星的輻射束半徑及誤差.發(fā)現(xiàn)脈沖星周期小于0.85 s時(shí),輻射束半徑正比于周期的?1/2次方.周期大于0.85 s時(shí),輻射束半徑與周期之間并沒有明顯的冪律關(guān)系,但這可能與選擇效應(yīng)有關(guān).當(dāng)假設(shè)脈沖星的磁傾角為90°時(shí),計(jì)算得到的輻射束半徑隨周期分布的下邊界正比于周期的?1/2次方,與前人的結(jié)論一致.

    脈沖星:普通,方法:統(tǒng)計(jì)

    1 引言

    脈沖星具有很強(qiáng)的磁場(chǎng),是高速自轉(zhuǎn)的磁中子星,其自轉(zhuǎn)周期P很穩(wěn)定.脈沖星的輻射一般用燈塔模型來解釋.當(dāng)脈沖星每自轉(zhuǎn)1周,其輻射束若掃過觀測(cè)者視線,就能觀測(cè)到脈沖信號(hào).單個(gè)脈沖信號(hào)的輪廓變化很大,但多個(gè)脈沖信號(hào)的累積輪廓非常穩(wěn)定.一般呈現(xiàn)出單峰、雙峰、三峰以及多峰等形態(tài)[1],可用經(jīng)驗(yàn)的核錐模型[2]或補(bǔ)丁模型[3]來解釋.最近,還提出了輻射束的風(fēng)扇模型[4].前人研究了核錐模型中核成分與錐成分的寬度關(guān)系[5?7].

    脈沖星的輻射束模型最早由Radhakrishnan等[8]提出,其幾何形狀如圖1所示.圖中3條軸線分別代表脈沖星的自轉(zhuǎn)軸、磁軸和觀測(cè)者的視線方向.自轉(zhuǎn)軸與磁軸之間的夾角是磁傾角α.當(dāng)輻射束繞自轉(zhuǎn)軸旋轉(zhuǎn)時(shí),磁軸與視線距離最近時(shí)所張的角度為撞擊角β.圖中ψ是偏振位置角,而?為輻射輪廓的相位.根據(jù)此模型,ψ隨?的變化應(yīng)該呈“S”曲線.輻射束半徑大小為ρ,圖中ST曲線在赤道圈上的投影對(duì)球心的張角值為輻射輪廓的寬度W.

    脈沖星輻射束半徑ρ可以通過脈沖寬度W以及磁傾角α、撞擊角β計(jì)算得到.前人研究了輻射束半徑與脈沖星自轉(zhuǎn)周期P的關(guān)系.Lyne等[3]取α為90°,得出ρ的下限值正比于周期的?1/3次方,并且通過分析說明輻射束更像是補(bǔ)丁模型.Rankin[2,9]基于核錐模型研究得出了在1 GHz附近,核成分與錐成分的角半徑的下限值都正比于周期的?1/2次方.最近,Maciesiak等[10?12]用具有中間脈沖的脈沖星進(jìn)行分析,同樣得出了核成分與錐成分的角半徑下限值都正比于周期的?1/2次方.

    圖1 脈沖星輻射束的幾何模型.圖形來自Lyne等[3]Fig.1 The geometry of the radiation beam model for pulsars taken from Lyne et al[3]

    前人的工作得出的結(jié)論主要是輻射束半徑的下限值與周期的關(guān)系,依賴于α接近于90°,而且輪廓寬度值的測(cè)量方式也有不同,并且他們的工作沒有給出具體的脈沖星輪廓選擇標(biāo)準(zhǔn).因此本文將用更多的嚴(yán)格標(biāo)準(zhǔn)選擇的脈沖星偏振輪廓來重新分析數(shù)據(jù),并且希望得出不依賴于α取值的輻射束半徑與周期之間的關(guān)系.

    2 數(shù)據(jù)的收集與處理

    本文收集的數(shù)據(jù)主要來自歐洲脈沖星多頻輪廓數(shù)據(jù)庫(EPN),其中包含了1 000多顆脈沖星的多種頻段的偏振數(shù)據(jù).偏振數(shù)據(jù)中包含了脈沖星的相位?,斯托克斯參量I、Q、U、V的信息.我們從中嚴(yán)格挑選出了87顆很好的脈沖星數(shù)據(jù)[13?17],選擇標(biāo)準(zhǔn)如下:(1)具有偏振輪廓;(2)輪廓的信噪比高,若為多頻觀測(cè),則選擇其中信噪比最好的頻率的輪廓圖;(3)偏振位置角曲線呈現(xiàn)比較明顯的“S”型.我們對(duì)選取的每一顆脈沖星作了偏振輪廓圖.圖2展示了其中的6顆.每個(gè)脈沖星輪廓圖的上部分為偏振流量輪廓圖,其中實(shí)線代表總流量I,點(diǎn)線代表線偏振強(qiáng)度虛線代表圓偏振強(qiáng)度V.中間部分為偏振位置角ψ隨相位?的變化,ψ=1/2tan?1(U/Q),偏振位置角的誤差由誤差傳遞公式計(jì)算得出,其中σU和σQ是用U和Q在脈沖輪廓外的噪音值計(jì)算得到.偏振位置角曲線有時(shí)會(huì)出現(xiàn)180°的躍變,或由于正交模式的影響[18]出現(xiàn)90°的躍變,我們對(duì)此躍變作了修正后,由圖2下半部分給出,并計(jì)算出dψ/d?的最大值,用直線畫出.

    圖2 6個(gè)脈沖星的偏振輪廓及偏振位置角曲線圖.其中每個(gè)脈沖星圖的上半部分為偏振輪廓曲線,實(shí)線為總流量,點(diǎn)線為線偏振,虛線為圓偏振.中間部分為偏振位置角曲線圖.下半部分為進(jìn)行180°或90°躍變修正后的偏振位置角曲線圖,直線為dψ/d?的最大值.Fig.2 Polarization and position angle profiles for 6 pulsars.For each pulsar,the upper panel shows the polarization profiles.The solid,dotted,and dashed lines represent the total intensity,the linear and circular polarization,respectively.The middle panel shows the position angle curves.The lower panel shows the position angle curves after the correction for 180° or 90° jumping,and the straight line represents dψ/d?|max.

    我們能從偏振輪廓圖得到輪廓的寬度,而脈沖寬度W的測(cè)量也有多種方式.一般情況下,W10為輪廓峰值取10%后輪廓所對(duì)應(yīng)的寬度,W50是峰值取50%后所對(duì)應(yīng)的寬度.若為雙峰以上型輪廓,則取距離最遠(yuǎn)的兩個(gè)峰之間的寬度為峰峰寬Wpp.對(duì)于脈沖輪廓為單峰的情況,脈沖寬度就只有W10和W50兩個(gè)值.表1列出了篩選出的87顆脈沖星的參數(shù),第1列是脈沖星的名字,第2列為脈沖星的觀測(cè)頻率,第3列為周期,第4列為參考文獻(xiàn),第5列至第7列分別為W10、W50和Wpp,第8列為偏振位置角曲線斜率最大值.

    由圖1所示的輻射束幾何圖可知,輻射束半徑ρ與磁傾角α、撞擊角β、輪廓寬度W之間的幾何關(guān)系為[3]

    而偏振位置角曲線的斜率的最大值為

    由偏振位置角曲線斜率的最大值公式,如果知道α的值,就由(2)式可得β值.再根據(jù)(1)式,測(cè)量出輪廓的寬度W,就能得出輻射束半徑ρ.考慮目前文獻(xiàn)中通過不同方法給出的α值通常有很大的不確定度,且樣本較少.為了使結(jié)果不依賴于α的取值,我們?cè)O(shè)定α的值分別為90°、60°、45°和30°.再根據(jù)脈沖星輪廓的形狀,采用3個(gè)(多峰)或兩個(gè)(單峰)寬度值,這樣每顆脈沖星就可以求得12或8個(gè)脈沖星輻射束半徑的值.我們把這12或8個(gè)值取平均作為脈沖星輻射束半徑的值,均方根作為誤差,就可以檢驗(yàn)脈沖星輻射束半徑與周期之間的關(guān)系.

    3 結(jié)果和分析

    圖3給出了87顆脈沖星輻射束半徑的分布圖.從圖中可以看出脈沖星輻射束半徑的值大部分分布在10°以下,在我們的樣本中占了83%,說明脈沖星的輻射束一般較窄.但也有例外.如PSR J1932+1059的ρ達(dá)到了32°.該脈沖星的dψ/d?|max值很小,說明β值較大,視線比較靠近輻射束的外圍.

    圖4給出了87顆脈沖星的輻射束半徑ρ隨周期P變化的關(guān)系圖.很容易看出在周期小于0.85 s時(shí)(圖4點(diǎn)線以左),ρ明顯有隨P增大而減小的趨勢(shì).但P大于0.85 s時(shí),這一變化不明顯.因此我們首先分析P小于0.85 s的脈沖星,共有49顆,如圖5所示.我們用加權(quán)最小二乘法擬合ρ=cPb,得出c=3.00°±0.47°,b=?0.47±0.17,接近于Rankin[9]得出的ρ=2.45°P?0.5的關(guān)系,如圖5中的實(shí)線所示.如果固定b的值分別為?0.5和?0.33進(jìn)行擬合,得出結(jié)果如圖5中的點(diǎn)劃線和虛線所示.可以看出我們的擬合結(jié)果非常接近于ρ∝P?0.5,而離ρ∝P?0.33的線較遠(yuǎn).

    圖3 脈沖星輻射束半徑分布圖Fig.3 The distribution histogram of the pulsar radiation beam radii

    圖4 87顆脈沖星的輻射束半徑與周期關(guān)系圖.其中點(diǎn)線表示P=0.85 s的分界線.Fig.4 The relation between pulsar beam radii and periods for 87 pulsars.The dotted line represents P=0.85 s.

    對(duì)于P大于0.85 s的脈沖星,共38顆,我們做了類似的擬合,結(jié)果如圖6所示.實(shí)線為擬合出的b=?0.02±0.25,接近于0,說明ρ與P的關(guān)系確實(shí)并不明顯,而且與b=?0.5和b=?0.33的點(diǎn)劃線與虛線都相差很遠(yuǎn).因此當(dāng)P大于0.85 s時(shí),我們的樣本顯示ρ與P不成冪律關(guān)系.但這可能與選擇效應(yīng)有關(guān).如果滿足ρ∝P?0.5的關(guān)系,則P越大,要求ρ越小,則觀測(cè)到脈沖星的概率就較小.按照前面所述的選擇標(biāo)準(zhǔn),我們選擇的較長(zhǎng)周期的脈沖星可能ρ會(huì)相對(duì)偏大,造成不滿足ρ∝P?0.5的關(guān)系.而且為了保證信噪比和偏振位置角的形狀,我們選擇的輪廓的觀測(cè)頻率并不相同,而有些脈沖星的輪廓寬度會(huì)隨著頻率有明顯的變化,這樣會(huì)導(dǎo)致寬度的彌散,這也有可能影響統(tǒng)計(jì)結(jié)果.

    圖5 49顆P<0.85 s脈沖星的輻射束半徑與周期關(guān)系圖.實(shí)線表示ρ=3.00°P?0.47,點(diǎn)劃線為ρ=2.94°P?0.5,虛線為ρ =3.36°P?0.33.Fig.5 The relation between pulsar beam radii and periods for 49 pulsars with P<0.85 s.The solid,dot-dashed,and dashed lines represent ρ =3.00°P?0.47,ρ =2.94°P?0.5,and ρ =3.36°P?0.33,respectively.

    圖6 38顆P>0.85 s脈沖星的輻射束半徑與周期關(guān)系圖.實(shí)線表示ρ=4.45°P?0.02,點(diǎn)劃線為ρ=5.17°P?0.5,虛線為ρ =4.94°P?0.33.Fig.6 The relation between pulsar beam radii and periods for 38 pulsars with P>0.85 s.The solid,dot-dashed,and dashed lines represent ρ =4.45°P?0.02,ρ =5.17°P?0.5,and ρ =4.94°P?0.33,respectively.

    為了與前人工作中所取ρ的下限值的結(jié)果比較,我們計(jì)算了當(dāng)α=90°和輪廓寬度為W10時(shí)的輻射束半徑值ρ90.圖7顯示了ρ90與P的關(guān)系.Rankin[9]提出ρ=2.45°P?0.5/sinα,若α=90°,則公式簡(jiǎn)化為ρ90=2.45°P?0.5. 圖7中的實(shí)線為ρ90=2.45°P?0.5,可以看出我們的結(jié)果與該下限關(guān)系符合得非常好.

    圖7 脈沖星輻射束半徑下限值與周期關(guān)系圖.實(shí)線為ρ90=2.45°P?0.5.Fig.7 The relation between the lower boundary of pulsar radiation beam radii and periods.The solid line represents ρ90=2.45°P?0.5.

    4 總結(jié)

    我們利用在歐洲脈沖星多頻輪廓數(shù)據(jù)庫網(wǎng)站已經(jīng)公開的數(shù)據(jù),篩選了87顆信噪比較好、偏振位置曲線較明顯為“S”型的脈沖星進(jìn)行了研究.對(duì)每顆脈沖星都作了偏振位置角曲線,找出最大斜率值.為了驗(yàn)證不依賴于磁傾角的輻射束半徑值ρ與周期P之間的關(guān)系,我們假設(shè)脈沖星的磁傾角α的值為90°、60°、45°和30°,取W10、W50和Wpp(多峰輪廓)或W10和W50(單峰輪廓)作為寬度值,求出了脈沖星輻射束半徑值及誤差,并用最小二乘法進(jìn)行擬合.發(fā)現(xiàn)在我們的樣本中,周期小于0.85 s的脈沖星符合ρ正比于P的?1/2次方的關(guān)系,大于0.85 s的脈沖星ρ未隨P的增大而明顯減小,但這在一定程度上受到選擇效應(yīng)的影響.若只取α的值為90°,我們求得的輻射束半徑隨周期分布的下邊界符合Rankin[9]提出的ρ90=2.45°P?0.5的下限關(guān)系.

    表1 87顆脈沖星的相關(guān)參數(shù)值Table 1 Parameters for 87 pulsars

    表1 續(xù)1Table 1 Continued 1

    表1 續(xù)2Table 1 Continued 2

    [1]Rankin J M.ApJ,1983,274:333

    [2]Rankin J M.ApJ,1993,405:285

    [3]Lyne A G,Manchester R N.MNRAS,1988,234:477

    [4]Wang H G,Pi F P,Zheng X P,et al.arXiv:1405.6825

    [5]Gould D M.Structure and Polarization of Pulsar Radio Emission Beams.Manchester:Univ.Manchester,1994:120

    [6]蔣雪,游霄鵬.天文學(xué)報(bào),2013,54:93

    [7]Jiang X,You X P.ChA&A,2014,38:32

    [8]Radhakrishnan V,Cooke D J.ApL,1969,3:225

    [9]Rankin J M.ApJ,1990,352:247

    [10]Maciesiak K,Gil J,Ribeiro V A R M.MNRAS,2011,414:1314

    [11]Maciesiak K,Gil J.MNRAS,2011,417:1444

    [12]Maciesiak K,Gil J,Melikidze G.MNRAS,2012,424:1762

    [13]Gould D M,Lyne A G.MNRAS,1998,301:235

    [14]Manchester R N,Han J L,Qiao G J.MNRAS,1998,295:280

    [15]Wu X J,Manchester R N,Lyne A G,et al.MNRAS,1993,261:630

    [16]Qiao G J,Manchester R N,Lyne A G,et al.MNRAS,1995,274:572

    [17]von Hoensbroech A,Xilouris K M.A&AS,1997,324:981

    [18]Cordes J M,Rankin J M,Backer D C.ApJ,1978,223:961

    A Study on the Relation between Radio Radiation Beam Radii and Periods of Pulsars

    WEI Ding-ding YOU Xiao-peng

    (School of Physical Science and Technology,Southwest University,Chongqing 400715)

    The period of pulsar is very stable.Some previous works have found that the pulsar radio radiation beam radius is proportional to the period with a power law index of?1/2 or?1/3.In this paper,87 pulsars with high signal-noise-ratio and obvious S shape curves of position angle are selected.The radius and error of radiation beam are calculated for each pulsar.We found that the distribution of pulsar beam radii is proportional to the periods with a power law index of?1/2 when the periods are smaller than 0.85 s.When periods are greater than 0.85 s,there is no signi ficant relation between pulsar beam radii and periods,but this is probably because of selection effect.When the inclination angle is assumed to be 90°,the lower boundary of beam radii is found to be proportional to the periods with a power law index of?1/2,which is consistent with the conclusions of previous works.

    pulsars:general,methods:statistical

    P145;

    A

    2014-06-03收到原稿,2014-06-16收到修改稿

    ?國家自然科學(xué)天文聯(lián)合基金項(xiàng)目(U1231120)資助

    ?yxp0910@swu.edu.cn

    猜你喜歡
    實(shí)線脈沖星偏振
    關(guān)于調(diào)整上海道路非必要超長(zhǎng)實(shí)線及高速監(jiān)控探頭強(qiáng)光燈建議
    “中國天眼”已發(fā)現(xiàn)740余顆新脈沖星
    軍事文摘(2023年12期)2023-06-12 07:51:00
    發(fā)現(xiàn)脈沖星的女天文學(xué)家——貝爾
    科學(xué)(2022年4期)2022-10-25 02:43:42
    秋天來啦
    戒煙
    詩潮(2019年8期)2019-08-23 05:39:48
    偏振糾纏雙光子態(tài)的糾纏特性分析
    電子制作(2019年12期)2019-07-16 08:45:20
    疊疊看 真神奇
    啟蒙(3-7歲)(2019年3期)2019-04-03 01:39:28
    基于LabVIEW的偏振調(diào)制激光測(cè)距測(cè)量控制系統(tǒng)
    偏振旋轉(zhuǎn)效應(yīng)在全光緩存器中的應(yīng)用與實(shí)現(xiàn)
    基于虛擬觀測(cè)值的X射線單脈沖星星光組合導(dǎo)航
    国产欧美亚洲国产| 天堂俺去俺来也www色官网| 欧美精品一区二区免费开放| 青春草亚洲视频在线观看| 国产成人精品福利久久| 精品卡一卡二卡四卡免费| 又大又爽又粗| 日韩 欧美 亚洲 中文字幕| 99精品久久久久人妻精品| 少妇精品久久久久久久| www日本在线高清视频| 精品久久久久久电影网| 看十八女毛片水多多多| 777久久人妻少妇嫩草av网站| 成年女人毛片免费观看观看9 | 亚洲第一区二区三区不卡| 在线观看免费视频网站a站| 女性被躁到高潮视频| 国产一区二区在线观看av| 亚洲精品日韩在线中文字幕| 一本一本久久a久久精品综合妖精| 蜜桃国产av成人99| 久久青草综合色| 久久99精品国语久久久| 亚洲人成77777在线视频| videos熟女内射| 九色亚洲精品在线播放| 日韩成人av中文字幕在线观看| 日韩大片免费观看网站| 国产免费福利视频在线观看| 十八禁网站网址无遮挡| 国产爽快片一区二区三区| 久久国产亚洲av麻豆专区| 狠狠婷婷综合久久久久久88av| 99久久人妻综合| 一区二区三区四区激情视频| 91aial.com中文字幕在线观看| 深夜精品福利| 黄色 视频免费看| 青春草亚洲视频在线观看| 肉色欧美久久久久久久蜜桃| 两个人免费观看高清视频| 菩萨蛮人人尽说江南好唐韦庄| 欧美激情极品国产一区二区三区| 在线观看www视频免费| 波多野结衣一区麻豆| a级片在线免费高清观看视频| 日韩欧美精品免费久久| 桃花免费在线播放| 日韩人妻精品一区2区三区| 午夜福利乱码中文字幕| 91精品伊人久久大香线蕉| 久久久久精品久久久久真实原创| 国产高清不卡午夜福利| 亚洲美女搞黄在线观看| 好男人视频免费观看在线| 亚洲天堂av无毛| 亚洲av日韩精品久久久久久密 | 亚洲男人天堂网一区| 国产精品国产三级国产专区5o| 久久99热这里只频精品6学生| 如日韩欧美国产精品一区二区三区| 久久久久久人人人人人| 97人妻天天添夜夜摸| 精品午夜福利在线看| 午夜福利网站1000一区二区三区| 国产精品.久久久| 在线观看www视频免费| 国产成人免费观看mmmm| 伊人久久国产一区二区| 韩国高清视频一区二区三区| 人人澡人人妻人| 满18在线观看网站| 日韩欧美精品免费久久| 中文字幕人妻熟女乱码| 高清黄色对白视频在线免费看| 日韩电影二区| 亚洲婷婷狠狠爱综合网| 国产精品偷伦视频观看了| 久久久国产精品麻豆| 考比视频在线观看| 蜜桃国产av成人99| 亚洲国产欧美网| svipshipincom国产片| 香蕉国产在线看| 精品国产乱码久久久久久男人| 青春草视频在线免费观看| 国产乱来视频区| 国产精品久久久久成人av| av一本久久久久| 91aial.com中文字幕在线观看| 亚洲成人av在线免费| 国产精品久久久人人做人人爽| 国产乱人偷精品视频| 久久久精品区二区三区| 啦啦啦中文免费视频观看日本| 国产亚洲午夜精品一区二区久久| 久久这里只有精品19| av网站免费在线观看视频| 黄色视频不卡| 成年人免费黄色播放视频| 免费看不卡的av| 一个人免费看片子| 嫩草影视91久久| 日韩欧美精品免费久久| 中文字幕最新亚洲高清| 亚洲精品乱久久久久久| 婷婷色综合大香蕉| 日本av免费视频播放| 国产av精品麻豆| 青春草视频在线免费观看| 国产乱人偷精品视频| 97精品久久久久久久久久精品| 一区在线观看完整版| 又大又爽又粗| 色94色欧美一区二区| 搡老岳熟女国产| 韩国精品一区二区三区| 人人妻,人人澡人人爽秒播 | 满18在线观看网站| 最近最新中文字幕大全免费视频 | 日韩欧美精品免费久久| 国产精品欧美亚洲77777| 爱豆传媒免费全集在线观看| 国产精品.久久久| 在线观看www视频免费| 一本大道久久a久久精品| 欧美久久黑人一区二区| 久久久精品国产亚洲av高清涩受| 久久久亚洲精品成人影院| 午夜福利在线免费观看网站| 久久精品国产亚洲av涩爱| 日韩欧美一区视频在线观看| 丝袜喷水一区| 亚洲精品中文字幕在线视频| 一边摸一边抽搐一进一出视频| 国产免费一区二区三区四区乱码| 久久97久久精品| 一区二区三区四区激情视频| 一边亲一边摸免费视频| 色吧在线观看| 黄片无遮挡物在线观看| 久久久国产精品麻豆| 久久精品aⅴ一区二区三区四区| 美女视频免费永久观看网站| 亚洲欧美中文字幕日韩二区| 高清欧美精品videossex| 99精品久久久久人妻精品| 精品国产一区二区三区久久久樱花| 九草在线视频观看| 无遮挡黄片免费观看| 国产日韩欧美亚洲二区| 久久久久久免费高清国产稀缺| 色婷婷av一区二区三区视频| 亚洲色图 男人天堂 中文字幕| 51午夜福利影视在线观看| 日韩,欧美,国产一区二区三区| 丁香六月天网| 女人被躁到高潮嗷嗷叫费观| 亚洲五月色婷婷综合| 人人澡人人妻人| 国产99久久九九免费精品| 国产高清国产精品国产三级| 又大又爽又粗| 国产在线免费精品| 丝袜喷水一区| av线在线观看网站| 亚洲五月色婷婷综合| 亚洲欧美日韩另类电影网站| 欧美日韩av久久| 日韩欧美精品免费久久| 久久久久精品国产欧美久久久 | 亚洲欧洲国产日韩| 亚洲一码二码三码区别大吗| 欧美日韩综合久久久久久| 国产精品.久久久| 亚洲第一青青草原| 久久久国产欧美日韩av| 亚洲国产欧美日韩在线播放| 妹子高潮喷水视频| 飞空精品影院首页| 老司机在亚洲福利影院| 欧美人与性动交α欧美软件| 久久99热这里只频精品6学生| 欧美少妇被猛烈插入视频| 成人三级做爰电影| 爱豆传媒免费全集在线观看| 亚洲精品,欧美精品| av网站免费在线观看视频| 丝袜美腿诱惑在线| 久久女婷五月综合色啪小说| 精品酒店卫生间| 黄片小视频在线播放| 国产成人午夜福利电影在线观看| 秋霞伦理黄片| 久久精品熟女亚洲av麻豆精品| 丰满饥渴人妻一区二区三| 国产 精品1| 久久人人爽人人片av| 欧美成人精品欧美一级黄| 国产在线一区二区三区精| 丁香六月天网| 国产精品香港三级国产av潘金莲 | 丰满饥渴人妻一区二区三| 999久久久国产精品视频| videos熟女内射| 十八禁高潮呻吟视频| 亚洲少妇的诱惑av| 麻豆精品久久久久久蜜桃| 下体分泌物呈黄色| 免费少妇av软件| 黑人欧美特级aaaaaa片| 欧美日韩视频高清一区二区三区二| 国产精品一区二区在线不卡| 欧美精品av麻豆av| 高清av免费在线| 少妇的丰满在线观看| 日韩中文字幕视频在线看片| 亚洲av福利一区| 大片电影免费在线观看免费| 亚洲一级一片aⅴ在线观看| 中国三级夫妇交换| 18禁裸乳无遮挡动漫免费视频| av免费观看日本| 日本色播在线视频| 丝袜人妻中文字幕| 超色免费av| 久久久久精品久久久久真实原创| 一本色道久久久久久精品综合| 男女边吃奶边做爰视频| 欧美日韩视频精品一区| 国产又色又爽无遮挡免| 日韩 欧美 亚洲 中文字幕| 久久精品熟女亚洲av麻豆精品| 欧美变态另类bdsm刘玥| 免费少妇av软件| 美女午夜性视频免费| 80岁老熟妇乱子伦牲交| 男女床上黄色一级片免费看| 亚洲精华国产精华液的使用体验| 深夜精品福利| 大香蕉久久网| 免费看不卡的av| 久久久精品94久久精品| 一级毛片电影观看| 青青草视频在线视频观看| 黄色毛片三级朝国网站| 午夜福利影视在线免费观看| av女优亚洲男人天堂| 欧美 亚洲 国产 日韩一| 高清不卡的av网站| 日韩免费高清中文字幕av| 亚洲欧美一区二区三区国产| 精品少妇久久久久久888优播| 亚洲欧洲精品一区二区精品久久久 | 国产片内射在线| 国产视频首页在线观看| 一区二区三区四区激情视频| 国产免费视频播放在线视频| 午夜福利视频精品| 欧美 亚洲 国产 日韩一| 亚洲欧美一区二区三区国产| 久热爱精品视频在线9| 国产精品一二三区在线看| 国产精品一二三区在线看| 自线自在国产av| 精品亚洲成a人片在线观看| 亚洲国产看品久久| 亚洲国产精品999| 亚洲av日韩精品久久久久久密 | av天堂久久9| 久久99热这里只频精品6学生| 国语对白做爰xxxⅹ性视频网站| 午夜激情av网站| 99热全是精品| 久久久久久人妻| 男女无遮挡免费网站观看| 欧美变态另类bdsm刘玥| 亚洲久久久国产精品| 免费女性裸体啪啪无遮挡网站| 一级片免费观看大全| 亚洲天堂av无毛| 亚洲精品乱久久久久久| 中文字幕最新亚洲高清| 美女福利国产在线| www日本在线高清视频| 国产一级毛片在线| 国产免费又黄又爽又色| 国产在线视频一区二区| 午夜福利,免费看| 免费观看av网站的网址| 超色免费av| 久久这里只有精品19| 久久久国产一区二区| 免费观看性生交大片5| 亚洲精品在线美女| 午夜福利一区二区在线看| av视频免费观看在线观看| 国产一区二区三区综合在线观看| 美女高潮到喷水免费观看| 国产成人精品福利久久| 女性生殖器流出的白浆| 欧美xxⅹ黑人| 久久精品国产亚洲av高清一级| 久久久亚洲精品成人影院| 欧美精品一区二区免费开放| 日本午夜av视频| 国产97色在线日韩免费| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品久久午夜乱码| 在线观看免费高清a一片| 在现免费观看毛片| 亚洲欧美激情在线| 如何舔出高潮| 极品少妇高潮喷水抽搐| 国产精品麻豆人妻色哟哟久久| 久久久国产欧美日韩av| 国产一卡二卡三卡精品 | 黄片小视频在线播放| 国产又色又爽无遮挡免| 精品国产露脸久久av麻豆| 一级,二级,三级黄色视频| 在线观看免费午夜福利视频| 中文字幕色久视频| av在线播放精品| 老司机影院成人| 亚洲欧美中文字幕日韩二区| 最近2019中文字幕mv第一页| 欧美日韩av久久| 欧美成人精品欧美一级黄| 亚洲欧美色中文字幕在线| 成人亚洲精品一区在线观看| 亚洲熟女精品中文字幕| svipshipincom国产片| 亚洲精品第二区| 日韩一区二区视频免费看| 亚洲精品aⅴ在线观看| 肉色欧美久久久久久久蜜桃| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩av免费高清视频| 日本爱情动作片www.在线观看| 国产精品一区二区精品视频观看| 又大又黄又爽视频免费| 十分钟在线观看高清视频www| av一本久久久久| 国产免费一区二区三区四区乱码| 老司机靠b影院| 国产一卡二卡三卡精品 | 啦啦啦观看免费观看视频高清 | 久久精品影院6| 免费看十八禁软件| 欧美中文日本在线观看视频| 久久国产亚洲av麻豆专区| 亚洲电影在线观看av| 国产精品久久久久久精品电影 | av中文乱码字幕在线| 亚洲av成人av| 久久精品成人免费网站| 一卡2卡三卡四卡精品乱码亚洲| 亚洲中文日韩欧美视频| 天天躁夜夜躁狠狠躁躁| 欧美最黄视频在线播放免费| 精品不卡国产一区二区三区| 久久精品aⅴ一区二区三区四区| 人成视频在线观看免费观看| www国产在线视频色| 欧美乱码精品一区二区三区| 亚洲精品在线观看二区| 久久人妻福利社区极品人妻图片| 97超级碰碰碰精品色视频在线观看| 精品久久久久久成人av| 亚洲精品av麻豆狂野| 国语自产精品视频在线第100页| 激情视频va一区二区三区| 88av欧美| 亚洲成av人片免费观看| 国产日韩一区二区三区精品不卡| 精品不卡国产一区二区三区| 国产精品久久久人人做人人爽| 91国产中文字幕| 亚洲av成人一区二区三| 男女下面进入的视频免费午夜 | 超碰成人久久| 美女大奶头视频| 电影成人av| 91老司机精品| 亚洲精品国产区一区二| 久久久久久国产a免费观看| 国产精品秋霞免费鲁丝片| netflix在线观看网站| www.自偷自拍.com| 香蕉丝袜av| 午夜福利成人在线免费观看| 亚洲成人免费电影在线观看| 国产一区二区在线av高清观看| 可以在线观看的亚洲视频| 黄色毛片三级朝国网站| 村上凉子中文字幕在线| 男女下面插进去视频免费观看| 在线播放国产精品三级| 午夜视频精品福利| or卡值多少钱| 久久精品国产99精品国产亚洲性色 | 中文字幕色久视频| 女同久久另类99精品国产91| 欧美黄色片欧美黄色片| 一级作爱视频免费观看| 黑人操中国人逼视频| 精品一品国产午夜福利视频| www.www免费av| 日日爽夜夜爽网站| 欧美av亚洲av综合av国产av| 国产精品 国内视频| 亚洲国产精品久久男人天堂| 欧美日韩瑟瑟在线播放| 中文字幕人成人乱码亚洲影| 欧美日韩乱码在线| 成熟少妇高潮喷水视频| 一边摸一边抽搐一进一出视频| 久久中文看片网| 免费在线观看亚洲国产| 国产精品国产高清国产av| 国产精品,欧美在线| 非洲黑人性xxxx精品又粗又长| 在线国产一区二区在线| 国产真人三级小视频在线观看| 母亲3免费完整高清在线观看| 天堂影院成人在线观看| 一区二区三区激情视频| 精品久久久久久久毛片微露脸| 可以免费在线观看a视频的电影网站| 午夜两性在线视频| 中文字幕人成人乱码亚洲影| 电影成人av| 成人欧美大片| 国产xxxxx性猛交| a在线观看视频网站| 国产1区2区3区精品| 午夜a级毛片| 国产精品美女特级片免费视频播放器 | 免费观看人在逋| 午夜免费鲁丝| 黄色 视频免费看| 精品免费久久久久久久清纯| 高清在线国产一区| 国产高清视频在线播放一区| 日韩欧美免费精品| 亚洲欧美日韩高清在线视频| 涩涩av久久男人的天堂| 日本 欧美在线| 高潮久久久久久久久久久不卡| 久久欧美精品欧美久久欧美| 亚洲精品久久国产高清桃花| 搞女人的毛片| 国产一区二区三区综合在线观看| 国产97色在线日韩免费| 亚洲 欧美 日韩 在线 免费| 欧美日韩亚洲综合一区二区三区_| 国产人伦9x9x在线观看| 久久欧美精品欧美久久欧美| 亚洲av第一区精品v没综合| 最新美女视频免费是黄的| ponron亚洲| 久久久久久久精品吃奶| 精品国产乱子伦一区二区三区| 国产av在哪里看| 成人手机av| 咕卡用的链子| 在线观看免费视频网站a站| 嫁个100分男人电影在线观看| 咕卡用的链子| 亚洲激情在线av| 国产成人欧美在线观看| 性色av乱码一区二区三区2| 成人18禁高潮啪啪吃奶动态图| 韩国精品一区二区三区| 久久久久九九精品影院| 午夜福利影视在线免费观看| 国产91精品成人一区二区三区| 两人在一起打扑克的视频| 国产精品亚洲一级av第二区| 国产精品九九99| √禁漫天堂资源中文www| 嫩草影视91久久| 久久久久九九精品影院| 九色亚洲精品在线播放| 91av网站免费观看| 成人18禁高潮啪啪吃奶动态图| 亚洲第一av免费看| e午夜精品久久久久久久| 少妇熟女aⅴ在线视频| 亚洲国产精品999在线| tocl精华| 国产精品一区二区精品视频观看| 淫秽高清视频在线观看| 久久热在线av| 国产在线观看jvid| 三级毛片av免费| 亚洲熟女毛片儿| 国内毛片毛片毛片毛片毛片| 人妻丰满熟妇av一区二区三区| 啦啦啦免费观看视频1| 欧美色视频一区免费| 成年版毛片免费区| 精品久久久久久久毛片微露脸| 日韩欧美国产在线观看| 亚洲国产毛片av蜜桃av| 亚洲成av人片免费观看| 熟妇人妻久久中文字幕3abv| 日韩欧美三级三区| 大码成人一级视频| 999久久久国产精品视频| 十八禁人妻一区二区| 亚洲专区字幕在线| 国产精品九九99| 亚洲情色 制服丝袜| 亚洲片人在线观看| 女人被狂操c到高潮| 999久久久国产精品视频| 大香蕉久久成人网| 欧美精品亚洲一区二区| 制服丝袜大香蕉在线| 中文字幕高清在线视频| 久久精品影院6| 午夜久久久久精精品| 99re在线观看精品视频| 国产亚洲欧美98| 国产精品影院久久| 男人操女人黄网站| 亚洲 欧美 日韩 在线 免费| 高清毛片免费观看视频网站| 亚洲人成网站在线播放欧美日韩| 丝袜人妻中文字幕| 国产亚洲精品久久久久久毛片| 大香蕉久久成人网| 亚洲人成网站在线播放欧美日韩| 国产成人精品久久二区二区91| 久久精品亚洲精品国产色婷小说| videosex国产| 黄片播放在线免费| 亚洲视频免费观看视频| 人妻久久中文字幕网| 国产成人系列免费观看| aaaaa片日本免费| 一区二区三区高清视频在线| 波多野结衣一区麻豆| 亚洲精品在线美女| 国产精华一区二区三区| 久久精品91蜜桃| 午夜福利在线观看吧| 少妇熟女aⅴ在线视频| 久久中文字幕一级| 国产亚洲精品综合一区在线观看 | 久99久视频精品免费| 亚洲中文字幕日韩| av超薄肉色丝袜交足视频| 又黄又粗又硬又大视频| 亚洲狠狠婷婷综合久久图片| 啦啦啦韩国在线观看视频| 亚洲第一青青草原| 国产精品永久免费网站| 免费少妇av软件| 国产单亲对白刺激| 免费看a级黄色片| 国产麻豆69| 国产精品综合久久久久久久免费 | 亚洲精品国产一区二区精华液| 亚洲成人精品中文字幕电影| 免费在线观看黄色视频的| 国产成人欧美| 桃红色精品国产亚洲av| 麻豆成人av在线观看| 免费人成视频x8x8入口观看| 搞女人的毛片| 久久人妻福利社区极品人妻图片| 9191精品国产免费久久| 琪琪午夜伦伦电影理论片6080| 黄片播放在线免费| 国产精品免费一区二区三区在线| 色尼玛亚洲综合影院| 亚洲国产看品久久| 可以免费在线观看a视频的电影网站| 午夜免费成人在线视频| 久久久精品国产亚洲av高清涩受| ponron亚洲| 操出白浆在线播放| 日本免费一区二区三区高清不卡 | 九色亚洲精品在线播放| 99国产精品一区二区蜜桃av| 一边摸一边抽搐一进一出视频| 欧美成人性av电影在线观看| 欧美绝顶高潮抽搐喷水| 日韩欧美免费精品| 久久精品亚洲精品国产色婷小说| 天天添夜夜摸| 女人爽到高潮嗷嗷叫在线视频| 国产一区二区在线av高清观看| 日日夜夜操网爽| 国产亚洲欧美在线一区二区| 亚洲色图综合在线观看| 叶爱在线成人免费视频播放| 男男h啪啪无遮挡| 久久久久久久午夜电影| 亚洲午夜理论影院| 国产成人精品久久二区二区91| 亚洲精品一卡2卡三卡4卡5卡| 久久久久国内视频| 亚洲精品国产区一区二| 国产熟女午夜一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 久久久精品国产亚洲av高清涩受| 人人妻人人澡欧美一区二区 | 老司机在亚洲福利影院| 久久天躁狠狠躁夜夜2o2o| 亚洲中文字幕日韩|