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

    耀變體中紅外日內(nèi)光變特征研究

    2023-11-17 06:57:28普治元易庭豐陳雨潼
    天文研究與技術(shù) 2023年6期

    普治元,易庭豐*,張 順,王 亮,陳雨潼

    (1. 云南師范大學(xué)物理與電子信息學(xué)院,云南 昆明 650500;2. 廣西相對(duì)論天體物理重點(diǎn)實(shí)驗(yàn)室,廣西 南寧 530004)

    耀變體是活動(dòng)星系核(Active Galactic Nucleus, AGN)中具有特殊觀測(cè)性質(zhì)的一個(gè)子類,已知它擁有強(qiáng)大的相對(duì)論性帶電粒子噴流,其方向很接近地球觀測(cè)者的視線[1]。耀變體的輻射主要來(lái)源于自身噴流的非熱輻射,少數(shù)來(lái)源于自身吸積盤或塵埃環(huán)的熱輻射,其范圍從射電波段一直延伸至伽馬波段,相應(yīng)的多波段光譜能量分布呈現(xiàn)一個(gè)典型的雙峰結(jié)構(gòu)。低能峰位于光學(xué)到紫外線或者軟X射線波段之間,來(lái)源于噴流內(nèi)部極端相對(duì)論性電子的同步輻射,所以又稱為 “同步峰”;高能峰位于硬X射線到Tev波段之間,它的產(chǎn)生目前有兩種解釋模型:輕子模型和強(qiáng)子模型[2-3]。

    光變是耀變體最基本的觀測(cè)特征之一,是指耀變體在不同波段、不同時(shí)間尺度中觀測(cè)星等或輻射流量發(fā)生變化。一般認(rèn)為,耀變體在所有可觀測(cè)波段中具有強(qiáng)烈、快速的大幅度光變,且光變時(shí)標(biāo)從幾分鐘到幾年不等。若耀變體的光變發(fā)生在幾分鐘到一天的時(shí)間范圍內(nèi),稱為日內(nèi)光變[8];發(fā)生在幾周到幾個(gè)月的時(shí)間范圍內(nèi),稱為短時(shí)標(biāo)光變(Short-Term Variability, STV);發(fā)生在幾個(gè)月到幾年的時(shí)間范圍內(nèi),稱為長(zhǎng)時(shí)標(biāo)光變(Long-Term Variability, LTV)[9]。3種光變中,耀變體日內(nèi)光變的產(chǎn)生仍然是一個(gè)存在爭(zhēng)議的問(wèn)題,有幾種模型認(rèn)為可能與自身的相對(duì)論性噴流或吸積盤存在緊密聯(lián)系,也可能與剪切層中等離子體的不穩(wěn)定有關(guān)[10]。文[11]先后對(duì)高峰頻耀變體在光學(xué)波段的日內(nèi)光變和低峰頻耀變體在X射線波段的日內(nèi)光變進(jìn)行了搜索,在高峰頻耀變體的144條光變曲線和低峰頻耀變體的50條光變曲線中,兩者都僅有4%的光變曲線存在真正的日內(nèi)光變。該結(jié)果似乎暗示了光譜能量分布的同步峰在耀變體日內(nèi)光變中扮演了一個(gè)重要的角色。

    關(guān)于耀變體紅外輻射的產(chǎn)生,目前有兩種主流觀點(diǎn),一種認(rèn)為它來(lái)源于吸積盤或塵埃環(huán)的熱輻射,另一種則認(rèn)為它來(lái)源于噴流內(nèi)部的非熱輻射[12]。對(duì)耀變體在紅外波段的光變特征進(jìn)行系統(tǒng)性的研究,可以幫助我們進(jìn)一步了解其內(nèi)部結(jié)構(gòu)、輻射區(qū)域和輻射物理機(jī)制。文[13]研究了費(fèi)米探測(cè)到的耀變體(以下簡(jiǎn)稱為費(fèi)米耀變體)和非費(fèi)米探測(cè)到的耀變體(以下簡(jiǎn)稱為非費(fèi)米耀變體)在中紅外波段的長(zhǎng)時(shí)標(biāo)光變。文[14]基于費(fèi)米大視場(chǎng)望遠(yuǎn)鏡探測(cè)到的第3個(gè)活動(dòng)星系核目錄(The Third Catalog of Active Galactic Nuclei Detected by the Fermi Large Area Telescope,3LAC)研究了伽馬射線耀變體在中紅外波段的日內(nèi)時(shí)標(biāo)光變和長(zhǎng)時(shí)標(biāo)光變,并根據(jù)他們自己的樣本探究了平譜射電類星體和蝎虎座BL型天體之間中紅外光變特征的相似性和差異性,比較了高峰頻耀變體、中峰頻耀變體和低峰頻耀變體的中紅外光變[15]。在本文的研究工作中,我們將費(fèi)米大視場(chǎng)望遠(yuǎn)鏡探測(cè)到的第4個(gè)活動(dòng)星系核目錄的增量版本(The Incremental Version of the Fourth Catalog of AGN Detected by the Fermi-LAT,4LAC-DR2)與第5期Roma-BZCAT多波段目錄(The 5th Edition of the Roma-BZCAT Multifrequency Catalogue of Blazars)[16]相結(jié)合,獲取了一個(gè)只包含耀變體的樣本目錄,主要目的是:(1)在擴(kuò)大樣本的基礎(chǔ)上,比較平譜射電類星體和蝎虎座BL型天體在中紅外波段的日內(nèi)亮度變化和顏色變化,分析其中的差異性和相似性;(2)按照高峰頻耀變體、中峰頻耀變體和低峰頻耀變體的分類,探究耀變體在中紅外波段的日內(nèi)光變是否與光譜能量分布的同步峰存在聯(lián)系;(3)研究費(fèi)米耀變體和非費(fèi)米耀變體在中紅外波段的日內(nèi)光變,并在不同的時(shí)間尺度上與文[13]的研究結(jié)果形成對(duì)比,旨在為兩者的不同提供更多的信息。

    1 數(shù)據(jù)樣本

    費(fèi)米伽馬射線太空望遠(yuǎn)鏡發(fā)射于2008年6月,其搭載的大視場(chǎng)望遠(yuǎn)鏡(Large Area Telescope, LAT)在GeV能量范圍內(nèi)的靈敏度比EGRET(Energetic Gamma-Ray Experiment Telescope)提高了20倍。為了及時(shí)探測(cè)到伽馬射線源,大視場(chǎng)望遠(yuǎn)鏡在每3 h內(nèi)不間斷地對(duì)整個(gè)天空進(jìn)行掃描,目前為止,F(xiàn)ermi官方已經(jīng)發(fā)布了多個(gè)伽馬射線源目錄,其中的4LAC-DR2目錄與4LAC[17]目錄(The Fourth Catalog of Active Galactic Nuclei Detected by the Fermi Large Area Telescope)相比,新增了285個(gè)活動(dòng)星系核。在4LAC-DR2目錄中,耀變體的占比最多,共3 437個(gè),其中包括1 308個(gè)蝎虎座BL型天體,744個(gè)平譜射電類星體,1 385個(gè)不確定類型的耀變體候選者(Blazar Candidates of Uncertain type, BCUs)。

    文[18]基于多波段測(cè)量和廣泛的文獻(xiàn)回顧,發(fā)布了目前為止最全面的耀變體目錄,即Roma-BZCAT目錄。第5期Roma-BZCAT目錄包含3 561個(gè)源的坐標(biāo)和多波段數(shù)據(jù),目錄中所有的源都可在射電波段探測(cè)到[16]。我們將第5期Roma-BZCAT目錄中的耀變體與4LAC-DR2目錄中的耀變體進(jìn)行了詳細(xì)的交叉匹配,發(fā)現(xiàn)有1 910個(gè)源不與4LAC-DR2目錄中的任何源存在聯(lián)系,我們將這些源作為非費(fèi)米耀變體,4LAC-DR2目錄中的3 437個(gè)源作為費(fèi)米耀變體,這5 347個(gè)源就是本文的樣本目錄。

    美國(guó)宇航局于2009年12月發(fā)射的廣域紅外巡天探測(cè)衛(wèi)星(Wide-Field Infrared Survey Explorer, WISE)搭載了一臺(tái)40 cm的紅外望遠(yuǎn)鏡,以6.1″,6.4″,6.5″和12.0″的空間分辨率在3.4 μm,4.6 μm, 12 μm和22 μm 4個(gè)中紅外波段(分別簡(jiǎn)稱為W1,W2,W3和W4波段)進(jìn)行了巡天觀測(cè)[19]。WISE在大約1天的時(shí)間內(nèi)繞軌道運(yùn)行15次,可以得到目標(biāo)源的天量級(jí)光變曲線,其中的星等數(shù)據(jù)基于織女星系統(tǒng)進(jìn)行校正[20]。ALLWISE項(xiàng)目結(jié)合了來(lái)自WISE和NEOWISE探測(cè)階段的所有數(shù)據(jù),基于4LAC-DR2目錄和第5期Roma-BACAT目錄提供的坐標(biāo),我們通過(guò)NASA/IPAC紅外科學(xué)檔案館的GATOR查詢服務(wù)從 “ALLWISE Multiepoch Photometry Table” 中(https://irsa.ipac.caltech.edu/),以2″為搜索半徑,對(duì)樣本目錄中耀變體的中紅外波段測(cè)光數(shù)據(jù)進(jìn)行了查找。

    在獲取的數(shù)據(jù)中,由于W3和W4波段的測(cè)光精度較差,數(shù)據(jù)多為 “null” 值,所以本文重點(diǎn)圍繞W1和W2波段的測(cè)光數(shù)據(jù)展開。為了盡量減少低質(zhì)量測(cè)光點(diǎn)對(duì)結(jié)果的影響,我們參考ALLWISE數(shù)據(jù)的官方解釋文檔和文[21]的做法,使用以下標(biāo)準(zhǔn)構(gòu)建光變曲線:(1)對(duì)于含有多個(gè)不同坐標(biāo)的中紅外波段測(cè)光數(shù)據(jù),我們僅保留與搜索時(shí)所用坐標(biāo)相差較小的同時(shí)測(cè)光點(diǎn);(2)剔除W1星等(W1星等誤差)或W2星等(W2星等誤差)等于 “null” 的同時(shí)測(cè)光點(diǎn);(3)剔除不滿足w1rchi2_ep(w2rchi2_ep)<5, nb<3, qi_fact=1, cc_flags=“0000”,na=0的同時(shí)測(cè)光點(diǎn)。

    根據(jù)上述標(biāo)準(zhǔn)初步構(gòu)建的光變曲線,為了消除宇宙射線對(duì)測(cè)光數(shù)據(jù)的影響,我們應(yīng)用3σ(σ為光變曲線的標(biāo)準(zhǔn)偏差)原則剔除每條光變曲線中的異常值。為了后續(xù)所研究的日內(nèi)亮度變化和顏色變化具有討論意義,對(duì)于只包含5個(gè)以下測(cè)光點(diǎn)或非日內(nèi)時(shí)間范圍的光變曲線,我們沒(méi)有統(tǒng)計(jì)。

    在最終的樣本數(shù)據(jù)中,共計(jì)6 424條中紅外波段天量級(jí)光變曲線,對(duì)應(yīng)3 114個(gè)耀變體,每條光變曲線中的每?jī)蓚€(gè)連續(xù)測(cè)光點(diǎn)的時(shí)間間隔均不超過(guò)1.2天。我們從4LAC-DR2目錄、第5期Roma-BZCAT目錄以及文[22]的研究工作中,獲取了部分源的常規(guī)類別和同步峰峰值頻率。將3 114個(gè)耀變體基于常規(guī)類別進(jìn)行分類,共有1 054個(gè)蝎虎座BL型天體、1 305個(gè)平譜射電類星體和755個(gè)不確定類型的耀變體候選者;基于同步峰峰值頻率進(jìn)行分類,共有346個(gè)高峰頻耀變體、356個(gè)中峰頻耀變體和1 054個(gè)低峰頻耀變體;基于是否被費(fèi)米衛(wèi)星探測(cè)進(jìn)行分類,共有1 815個(gè)費(fèi)米耀變體和1 299個(gè)非費(fèi)米耀變體。在3 114個(gè)耀變體中,具有紅移(z)的源共有2 092個(gè),范圍為1.26×10-3-6.80。

    表1總結(jié)了樣本數(shù)據(jù)中耀變體的相關(guān)信息。圖1按照常規(guī)類別展示了樣本數(shù)據(jù)中耀變體的紅移分布。圖2展示了樣本數(shù)據(jù)中所有光變曲線的測(cè)光點(diǎn)在W1和W2波段的觀測(cè)質(zhì)量。W1和W2波段共計(jì)150 214個(gè)測(cè)光點(diǎn),W1波段星等誤差分布的第25,50,75和95百分位分別為0.05 mag,0.07 mag,0.11 mag和0.22 mag;W2波段星等誤差分布的第25,50,75和95百分位分別為0.06 mag,0.10 mag,0.17 mag和0.33 mag。圖3展示了樣本數(shù)據(jù)中的3條光變曲線,從左到右分別是4FGL J0124.8-0625(BL Lac),5BZQ J1107+5219(FSRQ)和4FGL J0110.9+4344(BCU)。

    圖1 本次研究所使用的耀變體的紅移分布Fig.1 The redshift distribution of the Blazars used in this study

    圖3 中紅外波段天量級(jí)光變曲線的3個(gè)例子。(a)4FGL J0124.8-0625(BL Lac);(b)5BZQ J1107+5219(FSRQ);(c)4FGL J0110.9+4344(BCU)Fig.3 Three examples of Intra-Day light curves in the mid-infrared band. (a)4FGL J0124.8-0625 (BL Lac);(b)5BZQ J1107+5219 (FSRQ);(c)4FGL J0110.9+4344 (BCU)

    表1 本次研究所使用的耀變體的詳細(xì)信息Table 1 Details of the Blazars used in this study

    2 方 法

    對(duì)于樣本數(shù)據(jù)中的中紅外波段天量級(jí)光變曲線,W1星等的變化曲線可以用于研究亮度的變化。因?yàn)閃1和W2波段的觀測(cè)是同時(shí)的,所以某一測(cè)光點(diǎn)對(duì)應(yīng)的W1與W2的星等差值即為該測(cè)光點(diǎn)的顏色,W1-W2的變化曲線可以用于研究顏色的變化。本文利用3種方法探討樣本數(shù)據(jù)中耀變體的日內(nèi)亮度變化和顏色變化。

    2.1 光變幅度

    為了描述一個(gè)源的亮度變化,我們引入了文[23]定義的光變幅度。對(duì)于擁有N個(gè)測(cè)光點(diǎn)的光變曲線,可以先計(jì)算方差,然后通過(guò)移除方差中的觀測(cè)不確定性來(lái)計(jì)算其光變幅度σm,若σm>0,則表明亮度的變化幅度大于0。σm的計(jì)算公式為

    (1)

    其中,∑是光變曲線的方差,其計(jì)算公式為

    (2)

    其中,mi為觀測(cè)星等,〈m〉為觀測(cè)星等的加權(quán)平均值;的計(jì)算公式為

    (3)

    2.2 標(biāo)準(zhǔn)額外方差

    (4)

    其中,Ci為顏色;〈C〉為顏色的加權(quán)平均值;i(C)為顏色的誤差值。的誤差值通過(guò)

    (5)

    (6)

    2.3 參數(shù)V

    文[27]定義的參數(shù)V可以作為衡量亮度變化和顏色變化的指標(biāo),

    (7)

    其中,P為下不完全伽馬函數(shù);χ2為光變曲線的卡方值,

    (8)

    該式中yi為觀測(cè)星等(顏色),εi為觀測(cè)星等(顏色)的誤差值,〈y〉為觀測(cè)星等(顏色)的加權(quán)平均值。通過(guò)(7)式計(jì)算得到的V值越大,亮度(顏色)變化的概率越高,歷史文獻(xiàn)中常采用的標(biāo)準(zhǔn)為V>1.3,即變化概率高于95%[26]。

    3 結(jié) 果

    根據(jù)上述3種方法,我們對(duì)樣本數(shù)據(jù)中所有中紅外波段天量級(jí)光變曲線進(jìn)行了分析,若W1星等的變化曲線滿足V>1.3和σm>0,認(rèn)為表現(xiàn)出統(tǒng)計(jì)意義上顯著的亮度變化;若W1-W2的變化曲線滿足V>1.3和Δ>0,認(rèn)為表現(xiàn)出統(tǒng)計(jì)意義上顯著的顏色變化[28]。

    本文中,我們重點(diǎn)根據(jù)蝎虎座BL型天體和平譜射電類星體;高峰頻耀變體、中峰頻耀變體和低峰頻耀變體;費(fèi)米耀變體和非費(fèi)米耀變體的分類對(duì)光變曲線的分析結(jié)果進(jìn)行討論(見(jiàn)表2)。

    表2 W1波段中,表現(xiàn)出日內(nèi)亮度變化的耀變體的σm統(tǒng)計(jì)結(jié)果Table 2 Statistics of σm for Blazars exhibiting Intra-Day variability in the W1 band

    3.1 平譜射電類星體和蝎虎座BL型天體的日內(nèi)亮度變化和顏色變化

    我們對(duì)平譜射電類星體和蝎虎座BL型天體的σm進(jìn)行了K-S檢驗(yàn),其中D統(tǒng)計(jì)值和P值分別為0.31和6.59×10-28。圖4展示了兩者的σm直方分布和累積分布。

    圖4 W1波段中,平譜射電類星體和蝎虎座BL型天體在日內(nèi)時(shí)間尺度上σm的直方分布和累積分布Fig.4 Histogram and cumulative distributions of σm for FSRQs and BL Lacs on the Intra-Day time-scale in the W1 band

    為了進(jìn)一步了解平譜射電類星體和蝎虎座BL型天體的日內(nèi)顏色變化,我們對(duì)表現(xiàn)出亮度變化和顏色變化的光變曲線建立了顏色-星等圖,其中W1星等作為橫坐標(biāo),顏色(W1-W2)作為縱坐標(biāo)??紤]到顏色和星等的誤差值,我們采用了文[29]提出的算法擬合顏色-星等散點(diǎn)圖:(W1-W2)=AW1+B,A為斜率,B為截距。此外,我們還計(jì)算了W1星等和顏色(W1-W2)之間的斯皮爾曼等級(jí)相關(guān)系數(shù)rs和偶然概率p,采用以下標(biāo)準(zhǔn)表征顏色的變化:(1)如果A>0,rs>0.5,p< 0.05,則源表現(xiàn)為變亮變藍(lán)的趨勢(shì);(2)如果A<0,rs<-0.5,p<0.05,則源表現(xiàn)為變亮變紅的趨勢(shì)。按照此標(biāo)準(zhǔn),在蝎虎座BL型天體表現(xiàn)出亮度變化和顏色變化的光變曲線中,共有52條光變曲線傾向于變亮變藍(lán)趨勢(shì);在平譜射電類星體表現(xiàn)出亮度變化和顏色變化的光變曲線中,共有90條光變曲線傾向于變亮變藍(lán)趨勢(shì);在兩者表現(xiàn)出亮度變化和顏色變化的光變曲線中,均未發(fā)現(xiàn)有光變曲線傾向于變亮變紅趨勢(shì)。

    圖5展示了在平譜射電類星體和蝎虎座BL型天體的顏色-星等圖中,最佳擬合直線斜率A的直方分布。圖6展示了平譜射電類星體和蝎虎座BL型天體表現(xiàn)出變亮變藍(lán)趨勢(shì)的兩個(gè)例子,從左到右分別是4FGL J2054.8+0015(BL Lac)和4FGL J0634.9-2335(FSRQ)。

    圖5 在日內(nèi)時(shí)間尺度上,對(duì)蝎虎座BL型天體和平譜射電類星體的顏色-星等圖進(jìn)行最佳線性擬合的斜率A的直方分布Fig.5 Histogram distribution of the slope A of the best linear fit to the color magnitude plot of on the Intra-Day time-scale

    圖6 在日內(nèi)時(shí)間尺度上,蝎虎座BL型天體和平譜射電類星體表現(xiàn)出變亮變藍(lán)趨勢(shì)的兩個(gè)例子。(a)4FGL J2054.8+0015(BL Lac);(b)4FGL J0634.9-2335(FSRQ)Fig.6 Two examples of BL Lacs and FSRQs showing a brightening and bluing trend on the Intra-Day time-scale. (a)4FGL J2054.8+0015 (BL Lac);(b)4FGL J0634.9-2335 (FSRQ)

    3.2 高峰頻耀變體、中峰頻耀變體和低峰頻耀變體,費(fèi)米耀變體和非費(fèi)米耀變體的日內(nèi)亮度變化

    圖7 W1波段中,高峰頻耀變體、中峰頻耀變體和低峰頻耀變體在日內(nèi)時(shí)間尺度上σm的直方分布和累積分布Fig.7 Histogram and cumulative distributions of σm for HSPs,ISPs,and LSPs on the Intra-Day time-scale in the W1 band

    我們也對(duì)費(fèi)米耀變體、平譜射電類星體和蝎虎座BL型天體與非費(fèi)米耀變體、平譜射電類星體和蝎虎座BL型天體的σm進(jìn)行了K-S檢驗(yàn),其中D統(tǒng)計(jì)值、P值分別為0.20,0.19,0.34和1.50×10-16,1.51×10-6,1.06×10-12,圖8展示了兩者σm的直方分布和累計(jì)分布。

    圖8 W1波段中,費(fèi)米耀變體、平譜射電類星體和蝎虎座BL型天體與非費(fèi)米耀變體、平譜射電類星體和蝎虎座BL型天體在日內(nèi)時(shí)間尺度上σm的直方分布和累積分布Fig.8 Histogram and cumulative distributions of σm for Fermi-detected Blazars/FSRQs/BL Lacs and Non-Fermi-detected Blazars/FSRQs/BL Lacs on the Intra-Day time-scale in the W1 band

    4 總結(jié)與討論

    在本文的研究中,我們將4LAC-DR2目錄與第5期Roma-BZCAT目錄相結(jié)合,利用ALLWISE項(xiàng)目的公開數(shù)據(jù),構(gòu)建了一個(gè)包含6 432條中紅外波段天量級(jí)光變曲線的樣本數(shù)據(jù)。我們對(duì)該樣本數(shù)據(jù)中耀變體的日內(nèi)亮度變化和顏色變化進(jìn)行了分析,并按照平譜射電類星體、蝎虎座BL型天體,高峰頻耀變體、中峰頻耀變體和低峰頻耀變體,費(fèi)米耀變體、非費(fèi)米耀變體的分類,比較了其對(duì)應(yīng)的光變幅度,以此來(lái)研究不同類別耀變體的中紅外日內(nèi)光變的特征。

    文[15]研究了伽馬射線耀變體在中紅外波段的日內(nèi)光變和長(zhǎng)時(shí)標(biāo)光變,結(jié)果顯示,在中紅外波段,耀變體的日內(nèi)光變幅度整體上比長(zhǎng)時(shí)標(biāo)光變幅度小,在與文[13]研究的耀變體的中紅外長(zhǎng)時(shí)標(biāo)光變比較后,我們也得到了相同的結(jié)果。

    表現(xiàn)出日內(nèi)亮度變化的耀變體具體由平譜射電類星體、蝎虎座BL型天體和不確定類型的耀變體候選者組成,通過(guò)分析對(duì)應(yīng)的天量級(jí)光變曲線,我們發(fā)現(xiàn)3種類型光變曲線的連續(xù)數(shù)據(jù)點(diǎn)之間的最大時(shí)間間隔分別為0.66天、0.79天和0.66天。因?yàn)樵?個(gè)時(shí)間間隔中均檢測(cè)到了星等的變化,所以各類型光變曲線的光變時(shí)標(biāo)應(yīng)該遠(yuǎn)遠(yuǎn)小于這3個(gè)時(shí)間尺度。我們使用0.66天、0.79天和0.66天作為相應(yīng)類型光變曲線觀察到的光變時(shí)標(biāo)的上限,通過(guò)

    R

    (9)

    估算平譜射電類星體、蝎虎座BL型天體和不確定類型的耀變體候選者的中紅外波段輻射區(qū)域,其中,Δtint=Δtobs/(1+Z);c為光速。表現(xiàn)出日內(nèi)亮度變化的蝎虎座BL型天體、平譜射電類星體和不確定類型的耀變體候選者的紅移范圍分別為2.31×10-3~2.83,0.12~5.47和5.52×10-3~1.86。根據(jù)(9)式,蝎虎座BL型天體的中紅外波段輻射可能來(lái)自于1.45×10-4~5.53×10-4pc的空間尺度范圍;平譜射電類星體的中紅外波段輻射可能來(lái)自于1.03×10-4~5.92×10-4pc的空間尺度范圍;不確定類型的耀變體候選者的中紅外波段輻射可能來(lái)自于1.94×10-4~5.51×10-4pc的空間尺度范圍。該結(jié)果說(shuō)明樣本數(shù)據(jù)中耀變體的中紅外輻射區(qū)域遠(yuǎn)小于吸積盤的經(jīng)典半徑,故日內(nèi)光變并不是來(lái)源于吸積盤或塵埃環(huán),這或許暗示了觀測(cè)到的中紅外波段輻射是由源的非熱輻射主導(dǎo),考慮到部分耀變體并不是伽馬射線源,所以這一結(jié)論還需要更多的觀測(cè)數(shù)據(jù)檢驗(yàn)。目前為止,伽馬射線耀變體的中紅外日內(nèi)光變可以用基于噴流的模型解釋[30-31]。

    通過(guò)對(duì)樣本數(shù)據(jù)中耀變體的σm進(jìn)行統(tǒng)計(jì),我們發(fā)現(xiàn)在中紅外波段,表現(xiàn)出日內(nèi)亮度變化的低峰頻耀變體相比于中峰頻耀變體和高峰頻耀變體展現(xiàn)出更大的光變幅度,而中峰頻耀變體和高峰頻耀變體的光變幅度總體上相差不大,兩者的累計(jì)分布圖也很難區(qū)分,出現(xiàn)這一現(xiàn)象的原因可能與三者的同步峰有關(guān):低峰頻耀變體的同步峰峰值頻率在中紅外波段附近,因此其在中紅外波段具有較高的流量密度,相對(duì)能量輸出較大,而中峰頻耀變體和高峰頻耀變體的中紅外波段對(duì)應(yīng)于同步輻射成分的上升部分,因此兩者在中紅外波段具有較低的流量密度,相對(duì)能量輸出較小[7,32]。盡管樣本數(shù)據(jù)中的部分平譜射電類星體和蝎虎座BL型天體缺少具體的同步峰峰值頻率,但在表現(xiàn)出日內(nèi)亮度變化的平譜射電類星體中,大部分是低峰頻耀變體,而蝎虎座BL型天體則是中峰頻耀變體和高峰頻耀變體居多,因此平譜射電類星體相比蝎虎座BL型天體展現(xiàn)出更大的光變幅度。一些相關(guān)研究報(bào)道了平譜射電類星體具有更顯著的噴流和更大的多普勒因子[33-34],這些因素也可能對(duì)平譜射電類星體和蝎虎座BL型天體的中紅外日內(nèi)光變產(chǎn)生影響。

    對(duì)耀變體顏色變化的相關(guān)研究主要集中于光學(xué)和近紅外波段[35-36],本文基于中紅外波段的數(shù)據(jù),研究了平譜射電類星體和蝎虎座BL型天體的顏色變化趨勢(shì):在日內(nèi)時(shí)間尺度上,表現(xiàn)出亮度變化和顏色變化的平譜射電類星體和蝎虎座BL型天體主要傾向于變亮變藍(lán)趨勢(shì)。耀變體吸積盤中的局部溫度隨著吸積率的改變而發(fā)生變化,可能導(dǎo)致源表現(xiàn)出變亮變藍(lán)趨勢(shì)[37],由于觀察到的中紅外日內(nèi)光變?cè)谏鲜龇治鲋姓J(rèn)為源于耀變體自身的噴流,因此,我們發(fā)現(xiàn)的顏色變化趨勢(shì)可能與平譜射電類星體和蝎虎座BL型天體自身噴流的內(nèi)部物理過(guò)程有關(guān)。

    一些相關(guān)的工作研究了費(fèi)米耀變體和非費(fèi)米耀變體之間的區(qū)別:(1)在費(fèi)米蝎虎座BL型天體和非費(fèi)米蝎虎座BL型天體中,兩者的顯著區(qū)別就是前者擁有更長(zhǎng)的噴流;(2)費(fèi)米平譜射電類星體的噴流張角、核心區(qū)域偏振度通常比非費(fèi)米平譜射電類星體更大,并且擁有更高的核心區(qū)域亮溫度[38];(3)費(fèi)米耀變體似乎具有比非費(fèi)米耀變體更極端的物理參數(shù),例如更快的視噴流速度[39]、更大的多普勒因子[40]、更高的核心區(qū)域流量密度和亮溫度[41]等。文[13]研究了費(fèi)米耀變體和非費(fèi)米耀變體的中紅外長(zhǎng)時(shí)標(biāo)光變的差異性,結(jié)果顯示,在長(zhǎng)時(shí)標(biāo)時(shí)間尺度上,費(fèi)米耀變體相比于非費(fèi)米耀變體、平譜射電類星體和蝎虎座BL型天體表現(xiàn)出更大的光變幅度[13],而本文的研究結(jié)果顯示,表現(xiàn)出日內(nèi)亮度變化的費(fèi)米耀變體相比于非費(fèi)米耀變體、平譜射電類星體和蝎虎座BL型天體表現(xiàn)出更小的光變幅度。這兩個(gè)在不同時(shí)間尺度上的相反結(jié)果似乎暗示了耀變體的中紅外日內(nèi)光變和長(zhǎng)時(shí)標(biāo)光變具有不同的產(chǎn)生機(jī)制,雖然有研究表明耀變體存在光學(xué)、近紅外和伽馬射線輻射不相關(guān)的情況[42-44],但考慮到數(shù)據(jù)的選擇存在差異,所以本文得到的結(jié)果還需要對(duì)更多的中紅外波段數(shù)據(jù)進(jìn)行研究。

    亚洲精品日韩av片在线观看 | 成年版毛片免费区| 美女 人体艺术 gogo| 午夜精品在线福利| 欧美日韩黄片免| 亚洲av中文字字幕乱码综合| 精品国内亚洲2022精品成人| 欧美不卡视频在线免费观看| e午夜精品久久久久久久| 国产精品免费一区二区三区在线| 亚洲成人中文字幕在线播放| 麻豆国产97在线/欧美| 制服丝袜大香蕉在线| 给我免费播放毛片高清在线观看| 成年女人毛片免费观看观看9| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美zozozo另类| 午夜精品久久久久久毛片777| 日本在线视频免费播放| 国产伦在线观看视频一区| 色哟哟哟哟哟哟| 久久欧美精品欧美久久欧美| 18禁黄网站禁片免费观看直播| 国产三级在线视频| 亚洲专区国产一区二区| 亚洲天堂国产精品一区在线| 尤物成人国产欧美一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲中文字幕日韩| 天天一区二区日本电影三级| 亚洲欧美日韩高清在线视频| av片东京热男人的天堂| 中文字幕久久专区| 小说图片视频综合网站| 91久久精品电影网| 美女免费视频网站| 日韩欧美三级三区| 男女床上黄色一级片免费看| 国产一区二区在线观看日韩 | 91久久精品电影网| 亚洲自拍偷在线| 色吧在线观看| 一个人免费在线观看电影| 桃红色精品国产亚洲av| 亚洲国产精品久久男人天堂| 一本综合久久免费| 国产精品影院久久| 国产蜜桃级精品一区二区三区| 色吧在线观看| 3wmmmm亚洲av在线观看| 久久精品91蜜桃| 成人特级黄色片久久久久久久| 啦啦啦观看免费观看视频高清| 日本a在线网址| 九色成人免费人妻av| 久久久久亚洲av毛片大全| 亚洲熟妇中文字幕五十中出| 丝袜美腿在线中文| 亚洲欧美激情综合另类| 国内精品久久久久精免费| 成人永久免费在线观看视频| 国产午夜福利久久久久久| 国产精品久久久人人做人人爽| 亚洲av电影在线进入| 久久6这里有精品| 欧美黄色淫秽网站| 啦啦啦韩国在线观看视频| 精品国产亚洲在线| 久久精品夜夜夜夜夜久久蜜豆| 久久久久久久久久黄片| 亚洲国产欧洲综合997久久,| 国产一区二区亚洲精品在线观看| 丰满人妻一区二区三区视频av | 国产高清有码在线观看视频| ponron亚洲| 午夜福利高清视频| 国产午夜精品久久久久久一区二区三区 | 他把我摸到了高潮在线观看| 欧美日韩黄片免| 午夜福利成人在线免费观看| 欧美高清成人免费视频www| 内射极品少妇av片p| 精品电影一区二区在线| 99精品欧美一区二区三区四区| 国产亚洲精品综合一区在线观看| 成人一区二区视频在线观看| 99国产精品一区二区三区| 免费av观看视频| 国内久久婷婷六月综合欲色啪| 夜夜躁狠狠躁天天躁| 国产亚洲精品一区二区www| 国产探花极品一区二区| 亚洲真实伦在线观看| 欧美性猛交黑人性爽| 午夜精品久久久久久毛片777| 亚洲欧美日韩高清在线视频| 在线天堂最新版资源| 五月伊人婷婷丁香| av在线天堂中文字幕| 99精品久久久久人妻精品| 亚洲欧美日韩高清在线视频| 亚洲中文字幕日韩| 丰满的人妻完整版| 免费搜索国产男女视频| 成年免费大片在线观看| 成人欧美大片| 在线天堂最新版资源| 一级黄色大片毛片| 黄色日韩在线| 两个人看的免费小视频| 久久久久久久久久黄片| 欧美日韩乱码在线| 99在线视频只有这里精品首页| 老司机午夜十八禁免费视频| 一卡2卡三卡四卡精品乱码亚洲| 无人区码免费观看不卡| 亚洲精品456在线播放app | 日韩欧美国产一区二区入口| 天堂动漫精品| 性色av乱码一区二区三区2| 国产精品一及| 18禁裸乳无遮挡免费网站照片| 免费人成视频x8x8入口观看| 又爽又黄无遮挡网站| 国产精品野战在线观看| 亚洲专区中文字幕在线| 精品乱码久久久久久99久播| 免费观看的影片在线观看| 欧美3d第一页| 中文字幕高清在线视频| 熟妇人妻久久中文字幕3abv| 国产高潮美女av| 97碰自拍视频| 变态另类丝袜制服| 一级作爱视频免费观看| 99精品在免费线老司机午夜| 国产亚洲精品久久久久久毛片| 欧美又色又爽又黄视频| 婷婷六月久久综合丁香| 美女大奶头视频| 热99在线观看视频| 中国美女看黄片| 国产成人av教育| 亚洲国产精品久久男人天堂| 亚洲最大成人中文| ponron亚洲| 黄片小视频在线播放| 天美传媒精品一区二区| 真人一进一出gif抽搐免费| 国产 一区 欧美 日韩| 波多野结衣巨乳人妻| 久久久久久久久中文| 国产成人啪精品午夜网站| 最近视频中文字幕2019在线8| 欧美性猛交╳xxx乱大交人| 51国产日韩欧美| 级片在线观看| 日本一本二区三区精品| 精品福利观看| 色综合亚洲欧美另类图片| 午夜免费激情av| 色综合欧美亚洲国产小说| 午夜视频国产福利| 搞女人的毛片| 国内少妇人妻偷人精品xxx网站| www日本黄色视频网| 天天添夜夜摸| 国产成人欧美在线观看| 人人妻人人看人人澡| 久久精品91蜜桃| 午夜福利在线观看免费完整高清在 | 国产精品久久久久久久久免 | 成年女人永久免费观看视频| 亚洲18禁久久av| 久久久久久久久久黄片| 亚洲五月婷婷丁香| 色尼玛亚洲综合影院| 欧美日韩精品网址| 好男人电影高清在线观看| 亚洲熟妇中文字幕五十中出| 欧美色视频一区免费| 国内精品久久久久精免费| 在线天堂最新版资源| 亚洲激情在线av| 偷拍熟女少妇极品色| 夜夜看夜夜爽夜夜摸| 91九色精品人成在线观看| 国产午夜精品久久久久久一区二区三区 | 国内精品久久久久久久电影| 国产探花在线观看一区二区| 黄色日韩在线| 激情在线观看视频在线高清| 丰满乱子伦码专区| 午夜福利免费观看在线| 桃红色精品国产亚洲av| 首页视频小说图片口味搜索| 亚洲国产精品久久男人天堂| 国产亚洲精品一区二区www| 欧美绝顶高潮抽搐喷水| 亚洲国产精品合色在线| 黄色片一级片一级黄色片| 精品熟女少妇八av免费久了| 亚洲av不卡在线观看| 国产精品98久久久久久宅男小说| 免费av毛片视频| 最近最新中文字幕大全免费视频| x7x7x7水蜜桃| 国产色婷婷99| 日本成人三级电影网站| 亚洲18禁久久av| 成人性生交大片免费视频hd| 精品久久久久久久人妻蜜臀av| 国产午夜福利久久久久久| 一区福利在线观看| 国产精品免费一区二区三区在线| 色综合站精品国产| 禁无遮挡网站| 欧美极品一区二区三区四区| 又紧又爽又黄一区二区| 亚洲精品美女久久久久99蜜臀| 国产亚洲av嫩草精品影院| 最近最新中文字幕大全电影3| 在线天堂最新版资源| 男人舔女人下体高潮全视频| 免费大片18禁| 亚洲精品影视一区二区三区av| 国产 一区 欧美 日韩| 国产欧美日韩精品一区二区| 一进一出抽搐gif免费好疼| 少妇丰满av| 亚洲五月天丁香| 最近最新免费中文字幕在线| 99在线视频只有这里精品首页| 一级毛片高清免费大全| 日本黄色视频三级网站网址| 在线观看av片永久免费下载| 国产一区在线观看成人免费| 国产精品女同一区二区软件 | 欧美另类亚洲清纯唯美| 淫妇啪啪啪对白视频| 国产乱人伦免费视频| 久久国产精品影院| 亚洲专区中文字幕在线| 宅男免费午夜| 久久天躁狠狠躁夜夜2o2o| 国产久久久一区二区三区| 欧美日韩瑟瑟在线播放| 一a级毛片在线观看| 欧美成人a在线观看| 国产av一区在线观看免费| 一本久久中文字幕| 床上黄色一级片| 18禁裸乳无遮挡免费网站照片| 欧美大码av| 久久香蕉精品热| 国产蜜桃级精品一区二区三区| 日本a在线网址| 亚洲av不卡在线观看| 国产乱人伦免费视频| 欧美黄色片欧美黄色片| 欧美乱妇无乱码| 九色国产91popny在线| 国产黄片美女视频| 国产麻豆成人av免费视频| 欧美在线一区亚洲| 又爽又黄无遮挡网站| 精品福利观看| 黄色女人牲交| tocl精华| 午夜精品一区二区三区免费看| 精品国产超薄肉色丝袜足j| 热99在线观看视频| 日韩高清综合在线| 高清在线国产一区| 欧美日韩瑟瑟在线播放| 日韩成人在线观看一区二区三区| 日韩精品青青久久久久久| 又黄又粗又硬又大视频| 无遮挡黄片免费观看| 欧美黑人巨大hd| av黄色大香蕉| 九九在线视频观看精品| 天堂动漫精品| 两个人视频免费观看高清| 色综合欧美亚洲国产小说| 精品不卡国产一区二区三区| 亚洲激情在线av| 在线观看日韩欧美| 亚洲av成人精品一区久久| 成人特级黄色片久久久久久久| 亚洲av成人av| 久久欧美精品欧美久久欧美| 精品电影一区二区在线| 国产一区二区激情短视频| 国产麻豆成人av免费视频| 午夜视频国产福利| 欧美av亚洲av综合av国产av| 天天躁日日操中文字幕| 国产色婷婷99| 欧美一区二区精品小视频在线| 亚洲av一区综合| 少妇的丰满在线观看| 91麻豆av在线| 一级作爱视频免费观看| 国产麻豆成人av免费视频| 久久人妻av系列| 1000部很黄的大片| 国产黄a三级三级三级人| 午夜精品久久久久久毛片777| 中文资源天堂在线| 国产主播在线观看一区二区| 色精品久久人妻99蜜桃| 国产中年淑女户外野战色| 国产成人系列免费观看| 日本精品一区二区三区蜜桃| 青草久久国产| 国产精品野战在线观看| 午夜a级毛片| 国产精品嫩草影院av在线观看 | 人妻夜夜爽99麻豆av| 午夜亚洲福利在线播放| 天堂网av新在线| 3wmmmm亚洲av在线观看| 亚洲av免费高清在线观看| 国产精品乱码一区二三区的特点| 天堂动漫精品| 精品99又大又爽又粗少妇毛片 | 久久久久久久亚洲中文字幕 | 国产精品一及| 性欧美人与动物交配| 亚洲成人久久性| 久久天躁狠狠躁夜夜2o2o| 9191精品国产免费久久| 国产免费av片在线观看野外av| 亚洲欧美日韩高清专用| 国产毛片a区久久久久| 九九热线精品视视频播放| 日本一本二区三区精品| 日韩欧美精品v在线| 久久性视频一级片| 国产精品99久久久久久久久| 久久久久久久久久黄片| 色视频www国产| 不卡一级毛片| 亚洲aⅴ乱码一区二区在线播放| 国产精品野战在线观看| 国产欧美日韩一区二区三| 99久久精品国产亚洲精品| 国产一级毛片七仙女欲春2| 岛国在线免费视频观看| 看黄色毛片网站| 97超视频在线观看视频| 国产精品亚洲av一区麻豆| 亚洲 欧美 日韩 在线 免费| 亚洲七黄色美女视频| 91在线观看av| 99热只有精品国产| 特大巨黑吊av在线直播| 国产av在哪里看| 18禁黄网站禁片午夜丰满| 日本成人三级电影网站| 在线观看美女被高潮喷水网站 | 最近在线观看免费完整版| 一个人看视频在线观看www免费 | 免费人成在线观看视频色| 一a级毛片在线观看| 国产中年淑女户外野战色| 亚洲欧美日韩高清在线视频| 在线观看av片永久免费下载| 亚洲五月天丁香| 91在线观看av| 五月玫瑰六月丁香| aaaaa片日本免费| 九九热线精品视视频播放| 久久午夜亚洲精品久久| 少妇熟女aⅴ在线视频| 深夜精品福利| 久久精品国产99精品国产亚洲性色| 成人国产综合亚洲| 国产一区二区激情短视频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲最大成人中文| 国产三级黄色录像| 欧美成人免费av一区二区三区| 我的老师免费观看完整版| 99国产精品一区二区蜜桃av| 久久人妻av系列| 97碰自拍视频| 午夜免费激情av| 香蕉久久夜色| 国产成人系列免费观看| 精品人妻偷拍中文字幕| 国产伦人伦偷精品视频| 亚洲av中文字字幕乱码综合| 亚洲成人中文字幕在线播放| 热99re8久久精品国产| 国产午夜福利久久久久久| 18美女黄网站色大片免费观看| 欧美+日韩+精品| 亚洲无线在线观看| 高潮久久久久久久久久久不卡| 婷婷六月久久综合丁香| 少妇的丰满在线观看| 2021天堂中文幕一二区在线观| 亚洲国产色片| 一级黄色大片毛片| 天美传媒精品一区二区| 淫妇啪啪啪对白视频| 99久久久亚洲精品蜜臀av| 日本a在线网址| 亚洲中文日韩欧美视频| av国产免费在线观看| 在线观看免费午夜福利视频| 99视频精品全部免费 在线| 69av精品久久久久久| 亚洲成av人片在线播放无| 丰满乱子伦码专区| 国产成人a区在线观看| 国产免费av片在线观看野外av| 99久久综合精品五月天人人| 一区二区三区高清视频在线| 免费观看人在逋| www.熟女人妻精品国产| 啦啦啦免费观看视频1| 99在线视频只有这里精品首页| 欧美日本视频| 日韩欧美三级三区| 99久国产av精品| 怎么达到女性高潮| 久久久久免费精品人妻一区二区| 亚洲精品一区av在线观看| 性色avwww在线观看| 日韩国内少妇激情av| 长腿黑丝高跟| 99热这里只有精品一区| 国产99白浆流出| 岛国在线免费视频观看| 国产熟女xx| 久久久久久久久久黄片| 亚洲久久久久久中文字幕| 成人三级黄色视频| av片东京热男人的天堂| 两个人的视频大全免费| 国产精品爽爽va在线观看网站| 蜜桃久久精品国产亚洲av| 日本五十路高清| 午夜精品在线福利| 老司机深夜福利视频在线观看| 噜噜噜噜噜久久久久久91| 精品国内亚洲2022精品成人| 亚洲精品色激情综合| 久久久久精品国产欧美久久久| 久久婷婷人人爽人人干人人爱| 免费高清视频大片| 国产精品一区二区免费欧美| 色综合亚洲欧美另类图片| 日本三级黄在线观看| 国产精品香港三级国产av潘金莲| 1000部很黄的大片| 校园春色视频在线观看| 成人av一区二区三区在线看| 亚洲av一区综合| 国产探花在线观看一区二区| 国产成人系列免费观看| 亚洲五月婷婷丁香| www日本在线高清视频| 有码 亚洲区| 精品日产1卡2卡| 好男人电影高清在线观看| 午夜免费激情av| www日本在线高清视频| 免费观看人在逋| 无限看片的www在线观看| 国产精品精品国产色婷婷| 在线播放无遮挡| 琪琪午夜伦伦电影理论片6080| 国产精品一区二区免费欧美| 精品国产美女av久久久久小说| 少妇的丰满在线观看| 岛国在线观看网站| 欧美精品啪啪一区二区三区| 99久久精品热视频| 国产久久久一区二区三区| 女人高潮潮喷娇喘18禁视频| 国产美女午夜福利| 国产高潮美女av| 日本一二三区视频观看| 欧美日韩精品网址| 日韩高清综合在线| 国产一区二区三区视频了| 亚洲人成电影免费在线| 亚洲专区中文字幕在线| 一级黄片播放器| 男女午夜视频在线观看| 亚洲第一欧美日韩一区二区三区| 日韩 欧美 亚洲 中文字幕| 亚洲av不卡在线观看| 在线观看一区二区三区| 欧美成人免费av一区二区三区| 国产男靠女视频免费网站| 国产视频内射| 亚洲精品日韩av片在线观看 | 日韩欧美国产在线观看| 搡女人真爽免费视频火全软件 | 精品久久久久久成人av| 97人妻精品一区二区三区麻豆| 51国产日韩欧美| 日韩免费av在线播放| 欧美zozozo另类| 欧美精品啪啪一区二区三区| 亚洲av第一区精品v没综合| 欧美激情久久久久久爽电影| 日本a在线网址| 网址你懂的国产日韩在线| 亚洲精华国产精华精| 网址你懂的国产日韩在线| 男女视频在线观看网站免费| 老熟妇乱子伦视频在线观看| 欧美一区二区亚洲| 日本熟妇午夜| 99国产精品一区二区蜜桃av| 国语自产精品视频在线第100页| 啪啪无遮挡十八禁网站| 午夜免费激情av| 白带黄色成豆腐渣| 亚洲精华国产精华精| 99视频精品全部免费 在线| 午夜福利视频1000在线观看| 精品国产亚洲在线| 精品一区二区三区视频在线 | 操出白浆在线播放| 一进一出抽搐gif免费好疼| 午夜激情欧美在线| 人人妻人人澡欧美一区二区| 国产成人欧美在线观看| 午夜福利欧美成人| 久久伊人香网站| 99热只有精品国产| 中文字幕久久专区| 国产一区二区激情短视频| 亚洲成av人片在线播放无| 亚洲av一区综合| 婷婷亚洲欧美| 国模一区二区三区四区视频| 亚洲成人免费电影在线观看| 久久九九热精品免费| 在线观看66精品国产| 一本久久中文字幕| 国产成人a区在线观看| 国产午夜精品论理片| 成人特级av手机在线观看| 亚洲午夜理论影院| 亚洲av一区综合| 欧美乱色亚洲激情| 狠狠狠狠99中文字幕| 亚洲av五月六月丁香网| 亚洲五月天丁香| av福利片在线观看| 国模一区二区三区四区视频| 久久香蕉精品热| 国产淫片久久久久久久久 | 国产高清三级在线| 亚洲18禁久久av| 看片在线看免费视频| 中文字幕熟女人妻在线| 国产亚洲精品一区二区www| 成人午夜高清在线视频| 日本成人三级电影网站| 给我免费播放毛片高清在线观看| 国产精品98久久久久久宅男小说| 男插女下体视频免费在线播放| 国产伦一二天堂av在线观看| 久久久久久久精品吃奶| 国产乱人伦免费视频| 桃红色精品国产亚洲av| 每晚都被弄得嗷嗷叫到高潮| 黄色片一级片一级黄色片| 女人高潮潮喷娇喘18禁视频| 美女高潮的动态| 男人舔女人下体高潮全视频| 尤物成人国产欧美一区二区三区| 深夜精品福利| a级一级毛片免费在线观看| 国产一区二区三区在线臀色熟女| 日韩国内少妇激情av| 国模一区二区三区四区视频| 麻豆久久精品国产亚洲av| 亚洲欧美激情综合另类| 欧美bdsm另类| 99热这里只有是精品50| 露出奶头的视频| 欧美日韩瑟瑟在线播放| 黄片大片在线免费观看| 大型黄色视频在线免费观看| 亚洲第一欧美日韩一区二区三区| 色综合亚洲欧美另类图片| 久久精品国产综合久久久| 欧美成人性av电影在线观看| 在线天堂最新版资源| 在线免费观看不下载黄p国产 | 国产aⅴ精品一区二区三区波| 99在线人妻在线中文字幕| 欧美激情久久久久久爽电影| 99热6这里只有精品| 男女床上黄色一级片免费看| 老熟妇乱子伦视频在线观看| 99久久精品国产亚洲精品| 久久精品影院6| 在线a可以看的网站| 女人被狂操c到高潮| 久久久久久国产a免费观看| 夜夜看夜夜爽夜夜摸| 精品99又大又爽又粗少妇毛片 | 久久香蕉国产精品| 久久香蕉精品热|