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

    γ射線暴X射線耀發(fā)的研究進展

    2020-05-16 08:21:32劉傳璽毛基榮
    天文學(xué)進展 2020年1期
    關(guān)鍵詞:余輝噴流射線

    劉傳璽,毛基榮

    (1.中國科學(xué)院 云南天文臺,昆明 650011; 2.中國科學(xué)院大學(xué),北京 100049; 3.中國科學(xué)院 天體結(jié)構(gòu)與演化重點實驗室,昆明 650011)

    1 引言

    γ射線暴(gamma-ray burst,GRB)是宇宙中恒星尺度的最劇烈爆發(fā)現(xiàn)象,它在短時間內(nèi)會產(chǎn)生強烈的γ射線輻射。γ射線暴的光變曲線很復(fù)雜,可由單個脈沖構(gòu)成,也可由多個脈沖構(gòu)成。γ射線暴的時標從幾毫秒到幾千秒[1],其中,長暴的持續(xù)時間大于2 s,短暴的持續(xù)時間小于2 s[2]。Kumar和Zhang[3]曾提到,短暴的典型時間為0.3 s,長暴的典型時間是30 s。γ射線暴的光變曲線輪廓可以由 Norris函數(shù)F(t)=Fpλe?τ1/(t?ts)?(t?ts)/τ2描述,其中,F(xiàn)p是X射線耀發(fā)的峰值流量;λ=e2(τ1/τ2)1/2是歸一化因子,其作用是使峰值時間處X射線耀發(fā)流量歸一化為Fp;τ1和τ2是擬合參量,并且(τ1τ2)1/2+ts是X射線耀發(fā)的峰值時間tp;ts是X射線耀發(fā)的觸發(fā)時間;t是從γ射線暴觸發(fā)開始計時的時間[4]。觀測到的γ射線暴光譜可以由Band等人[5]提出的光譜經(jīng)驗函數(shù)描述,該函數(shù)由低能譜指數(shù)、峰值能量和高能譜指數(shù)描述[5]。Preece等人[6]統(tǒng)計了γ射線暴光譜經(jīng)驗函數(shù)中的低能譜指數(shù)、峰值能量和高能譜指數(shù)的分布,發(fā)現(xiàn)峰值能量分布在10~7000 keV,典型值是250 keV;低能譜指數(shù)分布在?2~0,典型值是?1;高能譜指數(shù)分布在?4~?1,典型值是?2.2。Norris等人[7]發(fā)現(xiàn),γ射線暴的峰值流量有光譜延時(spectral lag),低能段的峰值流量滯后于高能段的峰值流量,但脈沖輪廓在低能段和高能段中保持自相似性。我們還注意到,短暴的光譜比長暴的光譜硬[2]。Fishman等人[8]發(fā)現(xiàn),在銀河系的坐標中,γ射線暴分布是各向同性的,因此,γ射線暴應(yīng)該是在宇宙學(xué)距離上發(fā)生的,其平均光度是5×1044J·s?1[9–10]。

    γ射線暴從瞬時輻射轉(zhuǎn)入余輝輻射后,早期X射線余輝輻射的光變曲線依次出現(xiàn)快速衰減(initial steep decay)、緩慢衰減(shallow decay)、標準余輝衰減(normal decay)和晚期快速衰減(late steep decay)[11–12]。X射線耀發(fā)(X-ray flare)是X射線余輝中的閃耀,如圖1所示。

    迄今為止,γ射線暴的物理起源和中心引擎問題依然沒有答案。γ射線暴瞬時輻射與X射線耀發(fā)有相似的光變結(jié)構(gòu)。在γ射線暴瞬時輻射和X射線耀發(fā)脈沖結(jié)構(gòu)中,光譜演化均存在兩種模式:一種是光譜由硬變軟[13–14],另一種是光譜硬度與流量成正相關(guān)[15–16]。因此,通過X射線耀發(fā),我們可以了解中心引擎隨時間的演化過程以及γ射線暴的產(chǎn)生環(huán)境。部分γ射線暴中的X射線耀發(fā)亮度接近于γ射線暴瞬時輻射的亮度[17]。此外,研究高紅移γ射線暴的X射線耀發(fā)對研究早期宇宙也有重要意義。

    2 γ射線暴X射線耀發(fā)的觀測特征

    Swift是用于觀測γ射線暴的多波段望遠鏡,于2004年11月20日發(fā)射。BAT(Burst Alert Telescope)是搭載在Swift上的探測器,它能夠探測并快速定位γ射線暴,其精度可達幾角分。BAT的探測視場達到2.2 sr,探測能段為15~350 keV,定位精度達到3′。XRT(X-ray Telescope)的視場達到23.6′,能段為0.3~10 keV,定位精度達到5′′。XRT可以持續(xù)觀測γ射線暴X射線余輝,觀測的持續(xù)時間可以從BAT觸發(fā)γ射線暴后大約幾十秒到幾周,因此,通過XRT可以觀測早期和晚期的X射線耀發(fā)。XRT可以根據(jù)X射線余輝流量自動切換觀測模式:在WT(windowed timing)模式中,可以觀測計數(shù)率較高的輻射,時間分辨率為1.8×10?3s;在PC(photon counting)模式中,可以觀測計數(shù)率很低的輻射,時間分辨率為2.5 s。需要注意的是,當(dāng)X射線余輝流量增加時,PC模式和WT模式會出現(xiàn)數(shù)據(jù)重疊[18]。利用XRT的模式轉(zhuǎn)化可以有效地觀測X射線耀發(fā),防止X射線耀發(fā)流量過大而產(chǎn)生飽和。

    X射線耀發(fā)的流量是其最直觀的觀測量。我們把X射線耀發(fā)的流量與γ射線暴瞬時輻射的流量作對比,以研究X射線耀發(fā)與γ射線暴瞬時輻射的相似性。GRB 050502B中出現(xiàn)了迄今為止最亮的X射線耀發(fā),該X射線耀發(fā)出現(xiàn)在瞬時輻射后的12 min左右。GRB 050502B中的X射線耀發(fā)流量是X射線余輝流量的500倍,而且X射線耀發(fā)的光譜要比X射線余輝的光譜要硬[17]。O’Brien等人[20]發(fā)現(xiàn),明亮X射線耀發(fā)的流量可以達到瞬時輻射流量的量級。Margutti等人[21]從X射線耀發(fā)的平均光度方面研究了X射線耀發(fā)的性質(zhì)。他們選取的X射線耀發(fā)信息來自于Swift-XRT的觀測數(shù)據(jù)(2005―2010年,Swift-XRT共觀測到44個γ射線暴X射線耀發(fā),其紅移范圍為0.3~6.3)。他們發(fā)現(xiàn),平均光度?L?隨時間的演化關(guān)系如下:?L?=1054.5±0.1t?2.7±0.1,其中,30

    圖2 Swift-XRT觀測到的γ射線暴GRB 051117A余輝中X射線耀發(fā)的光變曲線

    X射線耀發(fā)可以在余輝輻射的任何階段被觀測到[22],其峰值時標和持續(xù)時標是檢驗γ射線暴光變模型的重要觀測量。此外,X射線耀發(fā)的光譜演化性質(zhì)可以用來檢驗X射線耀發(fā)和γ射線暴是否經(jīng)歷了相似的物理過程。GRB 050607 X射線余輝中有2個X射線耀發(fā),其出現(xiàn)時間分別為γ射線暴后135 s和310 s,持續(xù)時間分別是159 s和255 s。它們的上升階段與下降階段具有不對稱性,表現(xiàn)為快速上升和緩慢下降的光譜結(jié)構(gòu)。這2個X射線耀發(fā)都表現(xiàn)出明顯的光譜演化:出現(xiàn)X射線耀發(fā)時,光譜變硬;隨著X射線耀發(fā)流量衰減,光譜逐漸變軟。先出現(xiàn)的X射線耀發(fā)相對于后出現(xiàn)的X射線耀發(fā),其光譜較硬[23]。人們通過Swift-XRT對GRB 050822的觀測發(fā)現(xiàn),在瞬時輻射后有3個X射線耀發(fā),分別出現(xiàn)在γ射線暴后130 s,235 s,420 s。對于每一個X射線耀發(fā),它的硬度與流量正相關(guān),如圖3所示。Butler和Kocevski[15]也發(fā)現(xiàn),在X射線耀發(fā)期間,光譜的硬度隨流量的增加而增加。X射線耀發(fā)出現(xiàn)時光譜變硬,說明X射線耀發(fā)的光譜有時間演化效應(yīng)[24]。Bernardini等人[25]研究了晚期(t>1000 s)X射線耀發(fā)性質(zhì),并對比了早期(t<1000 s)X射線耀發(fā)的性質(zhì)。他們發(fā)現(xiàn),晚期X射線耀發(fā)與早期X射線耀發(fā)具有相似的w-tp關(guān)系和tdecay-trise關(guān)系(tdecay和trise分別為X射線耀發(fā)的下降時標和上升時標)。然而晚期X射線耀發(fā)與早期X射線耀發(fā)仍存在差異:晚期X射線耀發(fā)釋放的能量比早期X射線耀發(fā)釋放的能量低1個量級,且86%的晚期X射線耀發(fā)是由星際介質(zhì)的密度變化引起的[25]。

    γ射線暴中長暴和短暴的物理起源及其周圍的星際介質(zhì)都不相同。區(qū)分長暴和短暴中的X射線耀發(fā)對于了解γ射線暴的前身星和γ射線暴周圍的星際介質(zhì)具有重要意義。Mu等人[26]研究了短暴中明亮的X射線耀發(fā),發(fā)現(xiàn)Fp>3Fb(Fb是γ射線暴X射線余輝在X射線耀發(fā)峰值時間處的流量)。Margutti等人[27]比較了長暴和短暴中X射線耀發(fā)的觀測性質(zhì),發(fā)現(xiàn)短暴中X射線耀發(fā)的光度比長暴中X射線耀發(fā)的光度小1個量級;然而,在短暴和長暴中,X射線耀發(fā)的光譜隨時間的演化是一致的,X射線耀發(fā)在上升段光譜變硬,而在下降段光譜變軟。在短暴和長暴中,X射線耀發(fā)的w-tp關(guān)系也是一致的。

    圖3 GRB 050822余輝中X射線耀發(fā)的光變曲線

    γ射線暴余輝在X射線中出現(xiàn)增亮?xí)r,還可能在光學(xué)波段出現(xiàn)增亮[28]。通過多波段觀測余輝耀發(fā)的峰值時間和持續(xù)時間,可以研究X射線耀發(fā)與余輝在其他波段耀發(fā)的聯(lián)系。不同的輻射機制下,其余輝可能在不同能段發(fā)生增亮現(xiàn)象,因此,輻射機制的不同成分(同步輻射、熱輻射或逆康普頓散射)可能對余輝增亮有不同影響。Perri等人[29]發(fā)現(xiàn),在GRB 050730中同時出現(xiàn)X射線耀發(fā)和光學(xué)耀發(fā)。Yi等人[30]統(tǒng)計了γ射線暴光學(xué)耀發(fā)的時間分布(樣本來自Swift-UVOT 2005年4月―2010年12月的觀測數(shù)據(jù)[31–32],其中共有119個光學(xué)耀發(fā)),發(fā)現(xiàn)光學(xué)耀發(fā)與X射線耀發(fā)的峰值時間都分布在102~106s內(nèi),并有相似的w-tp關(guān)系。相似的時間統(tǒng)計性質(zhì)表明,γ射線暴的光學(xué)耀發(fā)與X射線耀發(fā)有聯(lián)系。然而,Becerra等人[33]用COATLI望遠鏡和XRT望遠鏡分別觀測了GRB 180205的光學(xué)余輝和X射線余輝,發(fā)現(xiàn)在X射線耀發(fā)出現(xiàn)時,光學(xué)余輝未出現(xiàn)增亮現(xiàn)象。此外,Wang和Dai[34]以及He等人[35]都發(fā)現(xiàn),GRB 100728A在X射線耀發(fā)期間還有吉電子伏能段的輻射,這說明在GRB 100728A的X射線耀發(fā)階段可能存在逆康普頓散射。但是,Troja等人[36]的統(tǒng)計研究表明,在X射線耀發(fā)期間兆電子伏至吉電子伏能段和光學(xué)能段沒有出現(xiàn)耀發(fā)現(xiàn)象。

    GRB 170817A在晚期出現(xiàn)不明亮的X射線耀發(fā)事件(Fp<3Fb),其發(fā)生的時間在暴后156 d。該X射線耀發(fā)的持續(xù)時間w是24 d,其峰值光度Lp為2×1032J·s?1[37]。GRB 170817A在爆發(fā)后156 d出現(xiàn)的X射線耀發(fā)的Lp與tp的關(guān)系,與Bernardini等人[25]統(tǒng)計的X射線耀發(fā)的Lp-tp關(guān)系一致,并且該X射線耀發(fā)的w/tp6 0.15。L等人[38]還給出了該X射線耀發(fā)的各向同性能量EX,iso和各向同性光度LX,iso,它們分別是3.13×1038J和1.54×1032J·s?1。L等人[38]對比GRB 170817A晚期X射線耀發(fā)和Swift觀測的X射線耀發(fā)的觀測特征,發(fā)現(xiàn)GRB 170817A晚期X射線耀發(fā)與其他γ射線暴中X射線耀發(fā)的Lp-EX,iso,Lp-LX,iso,Lp-tp,z,w-tp關(guān)系一致,其中,tp,z是紅移改正的峰值時間。該晚期X射線耀發(fā)可能由磁化的噴流引起的。

    2.1 X射線耀發(fā)的統(tǒng)計樣本和統(tǒng)計性質(zhì)

    Chincarini等人和Yi等人都建立了各自的X射線耀發(fā)的統(tǒng)計樣本,并以各自樣本為基礎(chǔ),統(tǒng)計了X射線耀發(fā)的觀測性質(zhì)。

    Chincarini等人[39]選取的樣本來自于Swift-XRT的觀測數(shù)據(jù),時間跨度為2005年4月―2008年3月。Chincarini等人選取X射線耀發(fā)基于下述原則:(1)X射線耀發(fā)要有完整的上升段、峰值和下降段;(2)X 射線耀發(fā)可以用 Norris函數(shù)F(t)=Fpλe?τ1/(t?ts)?(t?ts)/τ2擬合;(3)如果一個γ射線暴有多個X射線耀發(fā),每一個X射線耀發(fā)都可被分辨出來;(4)X 射線耀發(fā)必須是明亮的 (Fp>5×10?16J·m?2·s?1);(5)X 射線耀發(fā)的峰值時間必須小于1000 s。他們最終選取了113個X射線耀發(fā)樣本,對應(yīng)了56個γ射線暴,其中,GRB 051210和GRB 070724是短暴,21個γ射線暴標有紅移信息。

    Chincarini等人和Yi等人統(tǒng)計的X射線耀發(fā)樣本中,X射線耀發(fā)的性質(zhì)具有以下共性:(1)X射線耀發(fā)的持續(xù)時間隨X射線耀發(fā)能量的增加而減少;(2)X射線耀發(fā)的持續(xù)時間隨峰值時間的增加而線性增加;(3)X射線耀發(fā)的光譜比瞬時輻射的光譜軟;(4)X射線耀發(fā)的輻射強度隨峰值時間衰減;(5)X射線耀發(fā)的典型能量是1044J;(6)一個γ射線暴有多個X射線耀發(fā)時,后出現(xiàn)的X射線耀發(fā)的能譜要比先出現(xiàn)的X射線耀發(fā)的能譜軟。

    Yi等人篩選X射線耀發(fā)的標準與Chincarini等人的標準并不完全一致,差別在于:(1)Yi等人沒有對X射線耀發(fā)的峰值時間進行限制;(2)Yi等人擬合的X射線耀發(fā)函數(shù)與Chincarini等人的函數(shù)不同。Yi等人統(tǒng)計的X射線耀發(fā)樣本與Chincarini等人統(tǒng)計的X射線耀發(fā)樣本的不同之處是:(1)Chincarini等人統(tǒng)計的X射線耀發(fā)的峰值時間和持續(xù)時間都小于1000 s,Yi等人統(tǒng)計的X射線耀發(fā)的峰值時間為102~106s,并集中分布在100~1000 s,且X射線耀發(fā)的持續(xù)時間也集中分布在100~1000 s;(2)Chincarini等人統(tǒng)計的X射線耀發(fā)輪廓是非對稱的,trise/tdecay=0.49,Yi等人統(tǒng)計的X射線耀發(fā)輪廓是近似對稱的,trise/tdecay=1.08;(3)Chincarini等人提出,X射線耀發(fā)與γ射線暴瞬時輻射的觀測性質(zhì)相似,Yi等人[40]由統(tǒng)計得到的γ射線暴X射線耀發(fā)觀測性質(zhì)推斷,X射線耀發(fā)是由磁化噴流物質(zhì)釋放能量所引起的。

    Mu等人[41]從Berger[42]的67個短暴中選出了31個有Swift-XRT快速響應(yīng)(t<100 s)的短暴。此外,他們還從2013年1月―2017年8月的觀測數(shù)據(jù)中選出了有Swift-XRT快速響應(yīng)的 18 個短暴。他們采用函數(shù)F(t)=Fpλe?τ1/(t?ts)?(t?ts)/τ2,擬合這 49 個短暴中的 X 射線耀發(fā),發(fā)現(xiàn)GRB 050724,GRB 131004,GRB 161004中有明亮的X射線耀發(fā)(Fp>3Fb)。他們還發(fā)現(xiàn),GRB 131004和GRB 161004中的明亮X射線耀發(fā)服從w2+β(β是譜指數(shù)),這也說明這兩個短暴中的明亮X射線耀發(fā)可能與中心天體的活動有關(guān)。

    在不同天體中(如Sgr A*,M87,γ射線暴)人們都發(fā)現(xiàn)了X射線耀發(fā)[44]。Wang等人[44]統(tǒng)計了不同天體中出現(xiàn)的X射線耀發(fā)的特征,并與太陽產(chǎn)生的X射線耀發(fā)相比較,以尋找其中的聯(lián)系。他們發(fā)現(xiàn),這些天體的X射線耀發(fā)的能量和持續(xù)時間的累計分布都服從冪律函數(shù),且冪律指數(shù)差異不大,這表明噴流是由磁場主導(dǎo)的。

    3 γ射線暴X射線耀發(fā)的理論模型

    GRB 130925A中有X射線耀發(fā),而且在X射線耀發(fā)之后的X射線余輝中,出現(xiàn)了熱輻射和非熱輻射成分。當(dāng)噴流被熱的包層包裹時,X射線余輝中的熱輻射成分來自于包層的熱輻射[45]。

    γ射線暴的理論模型之一是火球模型( fireball model),包括內(nèi)激波機制產(chǎn)生γ射線暴[46–47]和外激波機制產(chǎn)生γ射線暴余輝[48]。Giannios[49]建立了耗散光球模型(photosphere),以解釋γ射線暴瞬時輻射的光變曲線和能譜。致密天體通過吸積過程噴出速度不等的殼層。如果噴流具有不同的磁化區(qū)域,在噴流和殼層運動過程中,噴流在光球半徑處產(chǎn)生熱輻射,磁場和激波的耗散產(chǎn)生非熱輻射[50]。X射線耀發(fā)的理論模型與γ射線暴的理論模型密切相關(guān)。X射線耀發(fā)可以用火球模型和光球模型解釋。磁場主導(dǎo)的噴流中,磁場的不穩(wěn)定性和磁重聯(lián)過程會加速噴流,并為X射線耀發(fā)提供能量。另外,噴流中能量的各向異性分布使得余輝出現(xiàn)增亮現(xiàn)象(增亮的持續(xù)時間與各向異性輻射區(qū)域的尺度有關(guān)),并可能導(dǎo)致X射線耀發(fā)。在致密天體的間歇性吸積過程中,磁場可以調(diào)節(jié)吸積率,并改變噴流能量,因此,我們可能觀測到一個γ射線暴的多次X射線耀發(fā)。

    3.1 火球模型

    3.1.1 內(nèi)激波機制

    3.1.2 外激波機制

    Fan和Wei[53]根據(jù)外激波的輻射,給出了X射線耀發(fā)的流量隨時間的衰減函數(shù)(tobs/trc)?(2+β),其中,trc是反向激波穿過殼層的時間。

    3.2 光球模型

    Beniamini和Kumar[58]考慮用光球模型解釋X射線耀發(fā)。他們根據(jù)該模型得出:(1)中心引擎活動只局限在很短的時間里;(2)瞬時輻射和X射線耀發(fā)是在它們各自光球半徑處發(fā)生的;(3)產(chǎn)生瞬時輻射和X射線耀發(fā)的噴出物質(zhì)是同時從中心引擎里噴出的;(4)兩類物質(zhì)有不同的速度。由此可以得出,瞬時輻射和X射線耀發(fā)的時標范圍和光球半徑的范圍,都與洛倫茲因子的分布有關(guān)。

    Ruffini等人[59]發(fā)現(xiàn),在X射線耀發(fā)的光譜中,有顯著的熱輻射的貢獻。Ruffini等人[60]提出,中子星在超新星噴射介質(zhì)中發(fā)生超吸積,最終坍縮成黑洞,并產(chǎn)生e+-e?等離子體。該等離子體同超新星噴射物質(zhì)一起演化到光學(xué)薄,并在光球半徑處產(chǎn)生了X射線耀發(fā)。熱輻射成分的出現(xiàn)表明,輻射物質(zhì)經(jīng)歷了從光學(xué)厚到光學(xué)薄的演化過程。

    3.3 磁耗散機制

    如果γ射線暴噴流是磁化的,那么,我們可以采用磁場耗散機制解釋X射線耀發(fā)現(xiàn)象。例如,Giannios[61]提出了一個產(chǎn)生短時標X射線耀發(fā)的模型,即當(dāng)噴流與外部介質(zhì)相互作用時,噴流內(nèi)磁場發(fā)生耗散。X射線耀發(fā)的時標與磁耗散的尺度緊密相關(guān),磁耗散的尺度越小,X射線耀發(fā)的時標越短。

    3.4 各向異性輻射機制

    γ射線暴噴流各向異性的輻射會引起X射線耀發(fā)的光變曲線的變化。Geng等人[62]提出,X射線耀發(fā)的快速衰減(α>2+β)是由于各向異性同步輻射引起的。由于磁場垂直于噴流運動方向,電子沿磁場作螺旋運動。如果電子的數(shù)量分布由電子運動方向與磁場方向的夾角決定,那么,電子的輻射就會出現(xiàn)各向異性。如果電子以小傾角沿著磁場螺旋運動,那么,我們就會看到觀測區(qū)域的邊緣要比中心明亮,因此,由于曲率效應(yīng),X射線耀發(fā)的峰值位置要晚出現(xiàn)。他們認為輻射面上的輻射主要集中在以視線方向為軸的圓環(huán)上,環(huán)的寬度會隨時間變寬,而該寬度的時標為T⊙,圓環(huán)劃過能量波動區(qū)域的動力學(xué)時標為TD,因此,各項異性輻射的相對持續(xù)時標為T⊙>TD[63]。噴流中能量波動尺度一旦滿足TD

    3.5 其他物理過程

    還有一些物理過程,主要包括吸積驅(qū)動噴流[64–65]、吸積盤與黑洞磁耦合[66]、磁星旋轉(zhuǎn)驅(qū)動噴流[67]和磁場間歇地驅(qū)動吸積盤[68],這些都可以用來解釋γ射線暴X射線耀發(fā)現(xiàn)象。

    Lazzati等人[64]模擬了γ射線暴噴流傳播的不穩(wěn)定性所產(chǎn)生的X射線耀發(fā)。噴流能量的變化會影響X射線耀發(fā)的光度及其隨時間的演化。在模擬過程中,他們設(shè)置了10°的噴流張角,并設(shè)置X射線耀發(fā)的洛倫茲因子下限為5。噴流光度按照函數(shù)L0x?5/3衰減,其中,L0是初始噴流光度,x是用于調(diào)節(jié)噴流光度的無量綱量。該過程模仿的是吸積過程中的物質(zhì)回落。他們推測,早期X射線耀發(fā)與噴流張角的增加有關(guān),而晚期X射線耀發(fā)與噴流張角的減小有關(guān)。X射線耀發(fā)的持續(xù)時間小于X射線耀發(fā)的峰值時間。

    Luo等人[66]提出,中心引擎為黑洞的低吸積率(0.001M⊙~0.1M⊙s?1)過程可能產(chǎn)生短時標的X射線耀發(fā)。該過程涉及吸積內(nèi)盤與中心黑洞之間的磁耦合機制(magnetic coupling,MC)。磁耦合過程被極偏角限制,并且其磁場強度隨盤半徑呈冪律變化,然而該模型不能解釋長時標的X射線耀發(fā)[66]。

    雙中子星并合后形成磁星。磁星附近扭曲的磁力線所貯存的磁能會轉(zhuǎn)變?yōu)檎撾娮拥入x子體,并以星風(fēng)的形式噴射出磁星表面[69–70]。星風(fēng)驅(qū)動雙中子星并合后的噴流物質(zhì),噴流物質(zhì)又作用在星際介質(zhì)上,最終形成不同物質(zhì)組成的層狀結(jié)構(gòu)。激波化星際介質(zhì)中的相對論電子產(chǎn)生同步輻射,而雙中子星并合后,拋出物質(zhì)的輻射是熱輻射[71]。對于X波段的流量,熱輻射先達到峰值,隨后是同步輻射達到峰值。

    坍縮恒星殘留的碎片或外部介質(zhì)都是吸積驅(qū)動的物質(zhì)來源[72]。吸積體(例如黑洞)周圍的磁場可以改變吸積的速率、吸積的時間以及傳輸?shù)絿娏髦械哪芰俊.?dāng)吸積流體的引力與磁力平衡時,吸積過程基本停止;隨后,吸積流體前端物質(zhì)被不斷累加,引力變大,磁場也不斷地被擠壓變形,最終吸積被重新建立,并產(chǎn)生噴流[68]。X射線耀發(fā)就是由間斷性磁場調(diào)制的吸積過程產(chǎn)生的。

    產(chǎn)生X射線耀發(fā)的中心天體可能是雙中子星并合后的較差旋轉(zhuǎn)的毫秒脈沖星。磁星的較差旋轉(zhuǎn)使極向磁場扭曲并形成環(huán)形磁場。隨著環(huán)形磁場的增強,磁場從恒星表面被拋出,磁場在噴流中的耗散產(chǎn)生X射線耀發(fā)[73]。磁星可以吸積周圍介質(zhì),回落的介質(zhì)會與磁星偶極磁場產(chǎn)生強烈作用。若介質(zhì)與磁星的距離為rm,回落物質(zhì)與磁星的共轉(zhuǎn)半徑為rc,那么,在介質(zhì)回落過程中,磁星吸積介質(zhì)后能否坍縮成黑洞,取決于rc與rm的相對位置:當(dāng)rc>rm時,磁星會坍縮成黑洞;當(dāng)rc

    4 觀測結(jié)果對模型限定的進一步討論

    除了采用理論模型解釋γ射線暴的X射線耀發(fā),我們還可以利用X射線耀發(fā)的觀測結(jié)果限制模型參數(shù)。Jin等人[76]提出,可以采用兩種方法估算X射線耀發(fā)的洛倫茲因子:(1)假設(shè)噴流是重子主導(dǎo)的,并且噴流被熱力學(xué)壓強加速,那么,最終的洛倫茲因子由γ射線暴噴流的光度和X射線耀發(fā)產(chǎn)生處的火球半徑限制。用該方法可以得出洛倫茲因子的上限。(2)假設(shè)X射線耀發(fā)的快速下降是由曲率效應(yīng)導(dǎo)致的,那么,X射線耀發(fā)的流量衰減可以用函數(shù)(?T/trc)?(2+β)描述(?T是噴流從兩個殼層先后噴出的間隔時間),且trc∝RX/(2Γ2c)。結(jié)合光變曲線可以得出,X射線耀發(fā)的洛倫茲因子分布范圍是10~1000。Yi等人[77]也根據(jù)X射線耀發(fā)的持續(xù)時間和噴流的動力學(xué)性質(zhì),給出了X射線耀發(fā)的洛倫茲因子的下限。Mu等人[78]利用曲率效應(yīng)研究了輻射區(qū)域與洛倫茲因子的關(guān)系。在計算流量時,涉及到多普勒效應(yīng)νobs=Γ(1+βccosθshell)νshell,其中,νobs是觀測到的光子頻率,βc=(1?1/Γ2)1/2,θshell是共動系的運動方向與視線方向的夾角,νshell是共動系中光子的頻率。多普勒效應(yīng)中包含了洛倫茲因子和緯度角(緯度角可以由衰減時間和輻射半徑替代),因此,Mu等人[78]建立了洛倫茲因子與光度的關(guān)系。他們利用觀測到的X射線耀發(fā)光度,估算出X射線耀發(fā)的洛倫茲因子分布在17~87,X射線耀發(fā)的輻射區(qū)域范圍為1015~1016cm。此外,我們注意到,GRB 060714的5個X射線耀發(fā)很難用外激波和內(nèi)激波機制解釋[79]。GRB 170817A晚期(156 d)出現(xiàn)的X射線耀發(fā)的觀測特征可以由部分理論模型給出合理的解釋。Lin等人[65]提出,用致密天體吸積薄盤模型可以解釋GRB 170817A晚期X射線耀發(fā)的峰值光度;Piro等人[37]認為,該X射線耀發(fā)可以由環(huán)形磁場的耗散產(chǎn)生;L等人[38]還提出,持續(xù)活動的超大質(zhì)量中子星可以使GRB 170817A在晚期產(chǎn)生X射線耀發(fā)。

    5 總結(jié)與展望

    本文主要介紹γ射線暴X射線耀發(fā)的研究進展,并著重介紹X射線耀發(fā)的觀測性質(zhì)以及相關(guān)理論模型。

    X射線耀發(fā)具有以下主要觀測性質(zhì):(1)快速上升和指數(shù)下降的光變結(jié)構(gòu);(2)整個X射線耀發(fā)的光變輪廓具有不對稱性;(3)X射線耀發(fā)的持續(xù)時間與其峰值時間有線性關(guān)系;(4)X射線耀發(fā)在上升時標內(nèi)的光譜較硬,而在下降時標內(nèi)的光譜較軟;(5)無論是長暴還是短暴,都可能出現(xiàn)X射線耀發(fā),但是短暴中X射線耀發(fā)輻射的能量要低于長暴中X射線耀發(fā)輻射的能量。

    用于研究X射線耀發(fā)的理論模型主要包括火球模型、光球模型、磁耗散機制模型、幾何結(jié)構(gòu)的非各向同性輻射模型和其他物理模型。

    今后對于γ射線暴X射線耀發(fā)的研究中,應(yīng)當(dāng)特別考慮觀測短暴中的X射線耀發(fā)。另外,應(yīng)當(dāng)根據(jù)γ射線暴的能譜分類建立不同的X射線耀發(fā)樣本,并進一步比較γ射線暴的瞬時輻射與X射線耀發(fā)的區(qū)別和聯(lián)系。

    猜你喜歡
    余輝噴流射線
    “直線、射線、線段”檢測題
    “慧眼”發(fā)現(xiàn)迄今距離黑洞最近的高速噴流
    『直線、射線、線段』檢測題
    Zn空位缺陷長余輝發(fā)光材料Zn1-δAl2O4-δ的研究
    光照條件對長余輝材料顯現(xiàn)手印效果影響的研究
    蓄能清潔人造石產(chǎn)品的研制
    赤石脂X-射線衍射指紋圖譜
    中成藥(2017年3期)2017-05-17 06:09:16
    噴流干擾氣動熱數(shù)值模擬的若干影響因素
    耀變體噴流高能電子譜的形成機制
    發(fā)生在活動區(qū)11931附近的重復(fù)噴流?
    美女扒开内裤让男人捅视频| 黄色视频在线播放观看不卡| 国产精品偷伦视频观看了| 岛国毛片在线播放| 这个男人来自地球电影免费观看| 久久天躁狠狠躁夜夜2o2o| 一区在线观看完整版| 欧美精品一区二区免费开放| 一本—道久久a久久精品蜜桃钙片| 伦理电影免费视频| netflix在线观看网站| 国产老妇伦熟女老妇高清| 亚洲精品在线美女| 亚洲精品在线观看二区| 黑丝袜美女国产一区| 亚洲黑人精品在线| 交换朋友夫妻互换小说| 亚洲伊人色综图| 少妇的丰满在线观看| 国产不卡av网站在线观看| 午夜福利视频精品| 国产伦理片在线播放av一区| 男女午夜视频在线观看| 亚洲伊人色综图| 欧美成人午夜精品| 757午夜福利合集在线观看| 色综合婷婷激情| 91大片在线观看| 亚洲精品国产色婷婷电影| 精品久久蜜臀av无| 亚洲色图综合在线观看| 天天躁日日躁夜夜躁夜夜| 女同久久另类99精品国产91| 在线观看一区二区三区激情| 91成年电影在线观看| 欧美日韩成人在线一区二区| √禁漫天堂资源中文www| 久久国产亚洲av麻豆专区| 免费在线观看视频国产中文字幕亚洲| 国产一区二区 视频在线| 人妻久久中文字幕网| 久久中文看片网| 高潮久久久久久久久久久不卡| 国产99久久九九免费精品| 最新美女视频免费是黄的| 麻豆国产av国片精品| 欧美成狂野欧美在线观看| 香蕉丝袜av| 男女免费视频国产| 日韩欧美一区二区三区在线观看 | 国产成人系列免费观看| 飞空精品影院首页| 老司机靠b影院| 婷婷成人精品国产| 欧美日韩黄片免| 欧美成人午夜精品| 一级黄色大片毛片| 曰老女人黄片| 超碰97精品在线观看| 自线自在国产av| 国产午夜精品久久久久久| 如日韩欧美国产精品一区二区三区| 性少妇av在线| 美女福利国产在线| 老汉色∧v一级毛片| 三上悠亚av全集在线观看| 欧美在线黄色| 国产有黄有色有爽视频| 夜夜夜夜夜久久久久| 人人妻人人添人人爽欧美一区卜| 黑人巨大精品欧美一区二区mp4| 在线观看人妻少妇| 99国产综合亚洲精品| 久久性视频一级片| 成年女人毛片免费观看观看9 | 国产亚洲欧美在线一区二区| 欧美精品一区二区大全| 午夜激情久久久久久久| 亚洲av片天天在线观看| 成人特级黄色片久久久久久久 | 久久久久久久国产电影| 亚洲精品久久成人aⅴ小说| 90打野战视频偷拍视频| 日本一区二区免费在线视频| 最新美女视频免费是黄的| 在线十欧美十亚洲十日本专区| 久久久久久人人人人人| 热99久久久久精品小说推荐| 欧美精品一区二区大全| 一进一出抽搐动态| 国产精品 欧美亚洲| 精品高清国产在线一区| 亚洲熟女毛片儿| 国产精品久久久av美女十八| 一个人免费在线观看的高清视频| 亚洲人成电影观看| 午夜激情久久久久久久| 色婷婷久久久亚洲欧美| 欧美+亚洲+日韩+国产| 日韩 欧美 亚洲 中文字幕| 国产精品麻豆人妻色哟哟久久| 国产欧美亚洲国产| 中文字幕人妻熟女乱码| av国产精品久久久久影院| 亚洲五月色婷婷综合| 国产一区有黄有色的免费视频| 国产一区二区激情短视频| 久久久久国内视频| 久久久久精品国产欧美久久久| 在线十欧美十亚洲十日本专区| 乱人伦中国视频| 国产一区有黄有色的免费视频| 色94色欧美一区二区| 如日韩欧美国产精品一区二区三区| 国产精品欧美亚洲77777| 亚洲情色 制服丝袜| 色播在线永久视频| 精品国产乱子伦一区二区三区| 午夜福利乱码中文字幕| 黑人巨大精品欧美一区二区蜜桃| 18禁观看日本| 亚洲伊人色综图| 国产日韩欧美在线精品| 精品一品国产午夜福利视频| 国产主播在线观看一区二区| 午夜福利一区二区在线看| 国产伦理片在线播放av一区| 国产欧美日韩精品亚洲av| 久久精品人人爽人人爽视色| 少妇裸体淫交视频免费看高清 | 黄色怎么调成土黄色| 搡老熟女国产l中国老女人| 在线观看www视频免费| 久久狼人影院| 日韩视频一区二区在线观看| bbb黄色大片| 高潮久久久久久久久久久不卡| 亚洲欧美精品综合一区二区三区| 日韩精品免费视频一区二区三区| 久久午夜亚洲精品久久| 少妇猛男粗大的猛烈进出视频| 丁香六月天网| 香蕉国产在线看| 中文字幕高清在线视频| 在线观看免费视频网站a站| 嫁个100分男人电影在线观看| 久久99热这里只频精品6学生| 嫩草影视91久久| 肉色欧美久久久久久久蜜桃| 欧美大码av| 少妇粗大呻吟视频| 啪啪无遮挡十八禁网站| 天堂8中文在线网| 久久精品国产a三级三级三级| 亚洲七黄色美女视频| 久久久精品94久久精品| 国产精品99久久99久久久不卡| 丰满人妻熟妇乱又伦精品不卡| 国产男女超爽视频在线观看| 99riav亚洲国产免费| 国产97色在线日韩免费| 三上悠亚av全集在线观看| 日韩成人在线观看一区二区三区| av有码第一页| 高清毛片免费观看视频网站 | 免费在线观看黄色视频的| √禁漫天堂资源中文www| 亚洲专区字幕在线| 999久久久国产精品视频| 亚洲视频免费观看视频| 天堂中文最新版在线下载| 黄色片一级片一级黄色片| 成人18禁高潮啪啪吃奶动态图| 夜夜骑夜夜射夜夜干| 成人国产一区最新在线观看| 中文字幕人妻丝袜一区二区| 免费看十八禁软件| 一区二区三区乱码不卡18| 成年动漫av网址| 成人黄色视频免费在线看| 亚洲av成人一区二区三| 国产精品免费一区二区三区在线 | 少妇的丰满在线观看| 亚洲欧美日韩另类电影网站| 大陆偷拍与自拍| 亚洲欧洲日产国产| 国产在线一区二区三区精| 亚洲av欧美aⅴ国产| 精品国产乱子伦一区二区三区| videos熟女内射| 麻豆av在线久日| 热re99久久国产66热| 99re在线观看精品视频| 午夜老司机福利片| 精品一区二区三区视频在线观看免费 | 一本色道久久久久久精品综合| 在线观看免费视频日本深夜| 国产一区二区三区在线臀色熟女 | 一本—道久久a久久精品蜜桃钙片| 在线观看免费视频日本深夜| 日本a在线网址| 久久热在线av| 新久久久久国产一级毛片| 欧美精品一区二区大全| 久久香蕉激情| 大码成人一级视频| 免费一级毛片在线播放高清视频 | 菩萨蛮人人尽说江南好唐韦庄| 最黄视频免费看| 777米奇影视久久| 怎么达到女性高潮| 曰老女人黄片| 韩国精品一区二区三区| 久久av网站| 久久精品91无色码中文字幕| 午夜成年电影在线免费观看| 欧美 亚洲 国产 日韩一| 97在线人人人人妻| 国产一区二区激情短视频| 久久久国产成人免费| 黄片播放在线免费| 黑人巨大精品欧美一区二区mp4| 80岁老熟妇乱子伦牲交| 亚洲国产欧美一区二区综合| 亚洲国产中文字幕在线视频| 一二三四在线观看免费中文在| 日韩欧美一区视频在线观看| 国产精品 国内视频| 18禁黄网站禁片午夜丰满| 老司机影院毛片| 叶爱在线成人免费视频播放| 麻豆国产av国片精品| 91老司机精品| 日韩一区二区三区影片| 他把我摸到了高潮在线观看 | 久久国产亚洲av麻豆专区| 纯流量卡能插随身wifi吗| 久久久水蜜桃国产精品网| 后天国语完整版免费观看| 黄片播放在线免费| 少妇粗大呻吟视频| 欧美日韩中文字幕国产精品一区二区三区 | 18禁国产床啪视频网站| 可以免费在线观看a视频的电影网站| 不卡一级毛片| 十八禁人妻一区二区| 亚洲成人手机| 欧美av亚洲av综合av国产av| 精品免费久久久久久久清纯 | 欧美黄色淫秽网站| 国产成+人综合+亚洲专区| 免费黄频网站在线观看国产| 黄频高清免费视频| 免费久久久久久久精品成人欧美视频| 精品亚洲成a人片在线观看| 国产黄频视频在线观看| av免费在线观看网站| 免费高清在线观看日韩| 久久99热这里只频精品6学生| 久久亚洲精品不卡| 性少妇av在线| 成人国产一区最新在线观看| 成人黄色视频免费在线看| 国产精品久久久久久精品电影小说| 午夜福利免费观看在线| 精品少妇久久久久久888优播| 精品视频人人做人人爽| 精品欧美一区二区三区在线| 精品亚洲成a人片在线观看| 香蕉国产在线看| 曰老女人黄片| 在线看a的网站| 亚洲av片天天在线观看| 女性生殖器流出的白浆| 怎么达到女性高潮| 午夜福利一区二区在线看| 首页视频小说图片口味搜索| 亚洲伊人色综图| 后天国语完整版免费观看| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩亚洲综合一区二区三区_| 亚洲精品久久午夜乱码| 久久av网站| 国产成人啪精品午夜网站| 国产男靠女视频免费网站| 一区在线观看完整版| 久热爱精品视频在线9| 香蕉丝袜av| 91精品国产国语对白视频| 91麻豆精品激情在线观看国产 | 国产色视频综合| 亚洲色图av天堂| 国产精品av久久久久免费| 亚洲情色 制服丝袜| 国产高清激情床上av| 日韩大片免费观看网站| 国产野战对白在线观看| 日本vs欧美在线观看视频| 少妇粗大呻吟视频| 久9热在线精品视频| 一级黄色大片毛片| 久久天躁狠狠躁夜夜2o2o| 不卡av一区二区三区| 国产一区二区三区在线臀色熟女 | 免费av中文字幕在线| 日韩成人在线观看一区二区三区| 十分钟在线观看高清视频www| 亚洲专区中文字幕在线| 国产熟女午夜一区二区三区| 老司机靠b影院| 久久久久国内视频| 久热爱精品视频在线9| 日韩欧美免费精品| 91九色精品人成在线观看| 国产欧美亚洲国产| 国精品久久久久久国模美| 亚洲视频免费观看视频| 国产片内射在线| 日韩 欧美 亚洲 中文字幕| 操美女的视频在线观看| 久久人人97超碰香蕉20202| 男人操女人黄网站| 12—13女人毛片做爰片一| 欧美日韩成人在线一区二区| 建设人人有责人人尽责人人享有的| 亚洲av欧美aⅴ国产| 欧美国产精品va在线观看不卡| 俄罗斯特黄特色一大片| 日韩免费高清中文字幕av| 色94色欧美一区二区| 国产av国产精品国产| 人人妻人人爽人人添夜夜欢视频| 久久久久久久久免费视频了| 99国产精品免费福利视频| 搡老乐熟女国产| 国产免费av片在线观看野外av| 久久久国产一区二区| 丁香欧美五月| 一级片'在线观看视频| 一个人免费看片子| 欧美一级毛片孕妇| 国产成人系列免费观看| 欧美成人免费av一区二区三区 | av视频免费观看在线观看| 亚洲avbb在线观看| 两个人免费观看高清视频| av片东京热男人的天堂| 国产精品国产高清国产av | 久久久水蜜桃国产精品网| 搡老乐熟女国产| av视频免费观看在线观看| 国产精品久久电影中文字幕 | 久久久欧美国产精品| 夫妻午夜视频| 国产不卡av网站在线观看| svipshipincom国产片| avwww免费| 一级片'在线观看视频| 国产又爽黄色视频| 午夜激情久久久久久久| 欧美黑人精品巨大| 国产片内射在线| 国产精品久久久人人做人人爽| 亚洲欧洲日产国产| 国产一卡二卡三卡精品| 9色porny在线观看| 中文字幕人妻丝袜一区二区| 99精品在免费线老司机午夜| 国产精品九九99| 丁香六月欧美| 少妇被粗大的猛进出69影院| 国产精品美女特级片免费视频播放器 | 90打野战视频偷拍视频| 99久久精品国产亚洲精品| 麻豆乱淫一区二区| 亚洲性夜色夜夜综合| 精品第一国产精品| 免费少妇av软件| 男人操女人黄网站| 国产精品久久久久久精品电影小说| 岛国毛片在线播放| www.熟女人妻精品国产| 国产在线视频一区二区| 夜夜夜夜夜久久久久| 999久久久精品免费观看国产| 男女免费视频国产| 少妇猛男粗大的猛烈进出视频| 考比视频在线观看| 久久久久久亚洲精品国产蜜桃av| 国产无遮挡羞羞视频在线观看| 高清黄色对白视频在线免费看| 国精品久久久久久国模美| 80岁老熟妇乱子伦牲交| 夜夜夜夜夜久久久久| 悠悠久久av| 大型av网站在线播放| 黑人欧美特级aaaaaa片| 巨乳人妻的诱惑在线观看| 欧美成狂野欧美在线观看| 老司机深夜福利视频在线观看| 欧美黄色片欧美黄色片| 桃花免费在线播放| 变态另类成人亚洲欧美熟女 | 亚洲熟女精品中文字幕| 色老头精品视频在线观看| 90打野战视频偷拍视频| 国产精品一区二区免费欧美| 久久国产精品男人的天堂亚洲| 在线观看免费午夜福利视频| 亚洲专区国产一区二区| 日韩视频一区二区在线观看| 人人澡人人妻人| 欧美在线黄色| 亚洲国产av新网站| 国产无遮挡羞羞视频在线观看| a级毛片在线看网站| 不卡av一区二区三区| 午夜福利在线免费观看网站| 国产av国产精品国产| 免费在线观看完整版高清| 老司机午夜福利在线观看视频 | 国产高清国产精品国产三级| 欧美 日韩 精品 国产| a级毛片在线看网站| 一本大道久久a久久精品| 99re6热这里在线精品视频| 女警被强在线播放| 午夜激情av网站| 免费高清在线观看日韩| 人人妻人人添人人爽欧美一区卜| 久热这里只有精品99| 久久精品国产99精品国产亚洲性色 | 国产人伦9x9x在线观看| 在线观看一区二区三区激情| 国产成人一区二区三区免费视频网站| 别揉我奶头~嗯~啊~动态视频| 悠悠久久av| 亚洲国产av新网站| 国产精品久久久av美女十八| 黄色丝袜av网址大全| 亚洲久久久国产精品| 少妇猛男粗大的猛烈进出视频| 叶爱在线成人免费视频播放| 欧美精品亚洲一区二区| 亚洲色图 男人天堂 中文字幕| 飞空精品影院首页| 成人国语在线视频| 日韩精品免费视频一区二区三区| 青青草视频在线视频观看| videosex国产| 757午夜福利合集在线观看| 视频在线观看一区二区三区| 亚洲国产成人一精品久久久| 久久久久视频综合| 亚洲成人国产一区在线观看| 久久 成人 亚洲| 久久中文字幕一级| 变态另类成人亚洲欧美熟女 | 99国产极品粉嫩在线观看| 国产成人欧美在线观看 | 亚洲精品自拍成人| 老汉色∧v一级毛片| 国产一区二区在线观看av| 日本欧美视频一区| 制服人妻中文乱码| 一个人免费在线观看的高清视频| 色婷婷久久久亚洲欧美| 国产色视频综合| 韩国精品一区二区三区| 欧美另类亚洲清纯唯美| 亚洲av美国av| 咕卡用的链子| 亚洲精品自拍成人| 亚洲 国产 在线| 国产单亲对白刺激| 久久久欧美国产精品| 日日夜夜操网爽| 中文字幕人妻丝袜制服| 亚洲欧美激情在线| 老司机午夜十八禁免费视频| 日韩欧美一区视频在线观看| 国产av国产精品国产| 狠狠婷婷综合久久久久久88av| 成年版毛片免费区| 国产熟女午夜一区二区三区| 亚洲欧美色中文字幕在线| 欧美大码av| 成人亚洲精品一区在线观看| 精品福利观看| 操出白浆在线播放| 亚洲第一av免费看| 我的亚洲天堂| 狠狠婷婷综合久久久久久88av| 亚洲精品成人av观看孕妇| 日本一区二区免费在线视频| 国产淫语在线视频| 国产av国产精品国产| 亚洲av第一区精品v没综合| 日韩一区二区三区影片| 性少妇av在线| 国产欧美日韩一区二区三| 国产av国产精品国产| 国产免费福利视频在线观看| 国产精品久久久久久精品古装| 国产97色在线日韩免费| a级片在线免费高清观看视频| 久久 成人 亚洲| 久久久久久久精品吃奶| 欧美 日韩 精品 国产| 精品久久蜜臀av无| 男女床上黄色一级片免费看| 欧美人与性动交α欧美软件| 欧美乱妇无乱码| 国产野战对白在线观看| 国产aⅴ精品一区二区三区波| 国产精品.久久久| 激情在线观看视频在线高清 | 黄色a级毛片大全视频| 啦啦啦 在线观看视频| 久久国产亚洲av麻豆专区| 久久国产精品大桥未久av| 久热爱精品视频在线9| 99久久国产精品久久久| 午夜福利一区二区在线看| 午夜福利视频精品| 久久精品91无色码中文字幕| 久久精品国产a三级三级三级| 久久久精品94久久精品| 人人妻人人澡人人爽人人夜夜| 亚洲精品一卡2卡三卡4卡5卡| 视频区图区小说| 80岁老熟妇乱子伦牲交| 精品卡一卡二卡四卡免费| 视频在线观看一区二区三区| 日本五十路高清| 久久午夜综合久久蜜桃| 欧美激情高清一区二区三区| 999久久久精品免费观看国产| 久久精品国产99精品国产亚洲性色 | 欧美乱妇无乱码| 亚洲国产毛片av蜜桃av| 精品国产乱子伦一区二区三区| 国产免费福利视频在线观看| 国产高清视频在线播放一区| 欧美中文综合在线视频| 女人久久www免费人成看片| 精品福利观看| 精品久久久久久电影网| 国产不卡av网站在线观看| 搡老熟女国产l中国老女人| 最近最新中文字幕大全免费视频| 亚洲色图av天堂| 99riav亚洲国产免费| 另类精品久久| 波多野结衣一区麻豆| 国产无遮挡羞羞视频在线观看| 亚洲av国产av综合av卡| 久久久久久久大尺度免费视频| 精品久久蜜臀av无| 黄色视频不卡| 下体分泌物呈黄色| 精品一区二区三区四区五区乱码| 18禁国产床啪视频网站| 亚洲自偷自拍图片 自拍| 亚洲色图 男人天堂 中文字幕| 国产真人三级小视频在线观看| 日日爽夜夜爽网站| 亚洲精品自拍成人| 日本黄色日本黄色录像| 美女国产高潮福利片在线看| 国产欧美日韩综合在线一区二区| 在线永久观看黄色视频| 99久久99久久久精品蜜桃| 亚洲一区二区三区欧美精品| 亚洲美女黄片视频| 精品少妇久久久久久888优播| 欧美精品一区二区大全| 亚洲精品在线观看二区| 国产精品一区二区免费欧美| 满18在线观看网站| 韩国精品一区二区三区| 丰满迷人的少妇在线观看| 亚洲七黄色美女视频| 大型黄色视频在线免费观看| 久久久精品94久久精品| 欧美激情高清一区二区三区| 12—13女人毛片做爰片一| 黑人猛操日本美女一级片| 精品国产一区二区三区久久久樱花| 亚洲精品国产区一区二| 国产精品久久电影中文字幕 | 纵有疾风起免费观看全集完整版| 日韩制服丝袜自拍偷拍| 亚洲情色 制服丝袜| cao死你这个sao货| 精品人妻1区二区| 亚洲欧美日韩高清在线视频 | 丝袜在线中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 国产欧美亚洲国产| 黑人巨大精品欧美一区二区mp4| 国产成人欧美在线观看 | 性少妇av在线| 亚洲精品美女久久久久99蜜臀| 水蜜桃什么品种好| 十分钟在线观看高清视频www| 国产欧美亚洲国产| 国产视频一区二区在线看| 亚洲国产欧美在线一区| 日韩熟女老妇一区二区性免费视频| 中文亚洲av片在线观看爽 | av免费在线观看网站| 18禁裸乳无遮挡动漫免费视频| 国产一区二区激情短视频| 建设人人有责人人尽责人人享有的|