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

    X射線雙星回旋吸收線研究進(jìn)展

    2019-07-17 02:27:52葉長青王德華張承民吳慶東
    天文學(xué)進(jìn)展 2019年2期
    關(guān)鍵詞:中子星基頻能譜

    葉長青,王德華,張承民,吳慶東

    (1.貴州師范大學(xué) 物理與電子科學(xué)學(xué)院,貴陽550025;2.貴州師范大學(xué) 貴州省射電天文數(shù)據(jù)處理重點(diǎn)實(shí)驗(yàn)室,貴陽550025;3.中國科學(xué)院 國家天文臺(tái),北京100101;4.中國科學(xué)院大學(xué),北京101400)

    1 引 言

    中子星是大質(zhì)量恒星演化到晚期,經(jīng)超新星爆發(fā)而形成的產(chǎn)物。中子星具有較強(qiáng)的磁場、引力場和超高的物質(zhì)密度,是研究極端環(huán)境下物理現(xiàn)象的天然實(shí)驗(yàn)室[1–4]。磁場是中子星重要的物理參數(shù)之一,對(duì)其測量的方法主要包括以下四種:(1)根據(jù)磁通守恒原理,由前身星磁場大小估計(jì)中子星磁場[5];(2)對(duì)于轉(zhuǎn)動(dòng)供能脈沖星,根據(jù)磁偶極輻射假設(shè)估計(jì)中子星磁場[6,7];(3)對(duì)于吸積供能中子星,根據(jù)磁球半徑與X 射線光度關(guān)系估計(jì)中子星磁場[8];(4)對(duì)于X 射線雙星,根據(jù)回旋共振散射特征(cyclotron resonance scattering feature,CRSF),即回旋吸收線特征,測量中子星磁場[9]。回旋吸收線方法是唯一能夠直接測量中子星磁場的方法[10]。

    X 射線雙星中,中子星通過吸積伴星物質(zhì)而發(fā)射X 射線(0.1~100 keV)。系統(tǒng)在寧靜態(tài)的光度范圍是1026~1028J·s?1,爆發(fā)態(tài)的光度能達(dá)到1031J·s?1[11]。在該系統(tǒng)的中子星表面附近,電子能量量子化至朗道能級(jí),并經(jīng)與光子的共振散射作用,在X 射線能譜中呈現(xiàn)吸收線特征:

    其中,Ecyc是回旋吸收線能量(常用單位為千電子伏),是普朗克常數(shù),e是電子電荷量,me是電子質(zhì)量,c是光速,B是磁場強(qiáng)度,B8是以108T 為單位的磁場強(qiáng)度。實(shí)際觀測到的回旋吸收線能量還受到中子星的引力場紅移的影響:

    其中,En為n階能量(n= 1,2,3,···),En,obs為觀測的回旋吸收線能量,z是引力紅移。根據(jù)回旋吸收線的發(fā)射位置,z為0.2~0.3[12]。

    1976年,Tr¨umper 等人[13]利用氣球上搭載的X 射線探測器,首次在Her X-1 的能譜中發(fā)現(xiàn)了回旋吸收線現(xiàn)象。1987年,日本發(fā)射了Ginga X射線衛(wèi)星,并在15 顆源中觀測到了回旋吸收線現(xiàn)象[14]。除此之外,RXTE[15],BeppoSAX[16],INTEGRAL[17],Swift[18],Suzaku[19]等衛(wèi)星也為回旋吸收線的觀測研究作出了巨大貢獻(xiàn)[10]。2012年,美國國家航空航天局(NASA)發(fā)射了高能譜分辨率的X 射線衛(wèi)星NuSTAR[20],并觀測到了一些新的回旋吸收線源。近年來,中國發(fā)射了具有較寬的能段范圍(1~250 keV)的硬X 射線調(diào)制衛(wèi)星——慧眼,為觀測更多高能量(≥80 keV)的回旋吸收線源提供了可能性[21,22]。截至目前,X 射線雙星中探測到的基態(tài)或基頻回旋吸收線能量范圍為10~80 keV,對(duì)應(yīng)的磁場強(qiáng)度范圍為107~109T[11]。

    探測技術(shù)的進(jìn)步使得人們?cè)诨匦站€的觀測方面取得了很多新進(jìn)展,例如,近年來人們發(fā)現(xiàn)一些源的二階諧頻與基頻回旋吸收線能量之比不是嚴(yán)格等于2;回旋吸收線形態(tài)具有復(fù)雜性,且形態(tài)參數(shù)之間存在相關(guān)性;回旋吸收線能量與源的X 射線光度間存在多重依賴關(guān)系;回旋吸收線能量隨著脈沖相位而發(fā)生變化;回旋吸收線存在長時(shí)標(biāo)演化等[11]。理論上,人們通常認(rèn)為,X 射線雙星中,吸積物質(zhì)沿磁力線到達(dá)中子星兩極并形成吸積柱,而回旋吸收線的發(fā)射因受到吸積柱高度、發(fā)射區(qū)磁場強(qiáng)度和吸積物質(zhì)的物理性質(zhì)等多種因素影響,其觀測參數(shù)之間可能呈現(xiàn)復(fù)雜的相關(guān)性[23,24]。例如,理論上人們認(rèn)為,吸積柱幾何和輻射形狀依賴于吸積率或X 射線光度[25]。對(duì)于高X 射線光度源,下落物質(zhì)由輻射主導(dǎo)的激波制動(dòng),因此,光子沿著垂直于吸積柱方向逃逸,并呈扇形的X 射線輻射形狀[26,27]。此時(shí),隨著光度和吸積柱高度的增加,回旋吸收線形成區(qū)的磁場強(qiáng)度下降,觀測到的回旋吸收線能量減小[28,29]。對(duì)于低X 射線光度源,自由下落的物質(zhì)在近中子星表面附近由庫侖相互作用制動(dòng),因此,光子沿著吸積柱方向逃逸,并呈鉛筆形的X 射線輻射形狀[26,27]。此時(shí),隨著光度和吸積柱高度降低,回旋吸收線形成區(qū)的磁場強(qiáng)度增加,觀測到的回旋吸收線能量增加。

    第2 章簡述回旋吸收線的形成原理;第3 章介紹回旋吸收線的最新觀測進(jìn)展及其理論解釋,包括諧頻與基頻回旋吸收線能量的比、回旋吸收線形狀的復(fù)雜性及形態(tài)參數(shù)間的相關(guān)性、回旋吸收線能量與光度的關(guān)系、回旋吸收線的脈沖相位解析譜及回旋吸收線的長時(shí)標(biāo)演化;第4 章進(jìn)行總結(jié)與展望。

    2 回旋吸收線的形成

    在中子星的強(qiáng)磁場環(huán)境下,電子繞磁力線做回旋運(yùn)動(dòng)的拉莫半徑rL和拉莫頻率ωL分別是:

    其中,v⊥是電子速度垂直于磁場的分量[12]。一般認(rèn)為,當(dāng)磁場達(dá)到臨界值Bcrit=m2ec3/(e)=44×108T,即電子的回旋能量與其靜止能量相當(dāng)(ωL≈mec2)時(shí),電子動(dòng)量在垂直于磁場方向的分量P⊥需要進(jìn)行相對(duì)論量子力學(xué)修正,得到量子化的朗道動(dòng)量:

    這時(shí),電子的散射截面在不同朗道能級(jí)差處發(fā)生共振并加強(qiáng)(見圖1),同時(shí),由于朗道能級(jí)的熱展寬作用,光子的能量在接近兩個(gè)電子朗道能級(jí)的差值時(shí),會(huì)因電子的能級(jí)躍遷作用而被吸收,因此,X 射線能譜具有吸收線特征:

    其中,θ是光子入射方向與磁力線方向的夾角,是引力紅移量,G是引力常數(shù),M?和R?分別是中子星的質(zhì)量和半徑,n是量子數(shù),n=1 代表基頻回旋吸收線能量,n=2,3,4,···分別代表2,3,4,···階諧頻回旋吸收線能量。

    圖1 電子躍遷時(shí)散射截面隨靜止能量變化圖[30]

    回旋吸收線在X 射線能譜中表現(xiàn)為疊加的n階吸收線特征。對(duì)能譜擬合時(shí),人們常用冪律函數(shù)加指數(shù)截?cái)嗪瘮?shù)來描述X 射線雙星的連續(xù)譜,例如:

    其中,F(xiàn)是光子數(shù)流量;E是能量;α是譜光子指數(shù)(冪律指數(shù)),其值通常為0.5~2.0;Ec是截?cái)嗄芰?,其值通常?~30 keV;A是系數(shù);Ef是折疊能。以下是幾種常用的能譜擬合模型:(1)冪律函數(shù)加費(fèi)米-狄拉克形式的指數(shù)截?cái)嗪瘮?shù)模型(Fermi-Dirac cutoff,FDCO)[31];(2)負(fù)(和正)冪律函數(shù)加指數(shù)截?cái)嗪瘮?shù)模型(negative and positive power laws with a common exponential cutoff factor,NEPX)[14];(3)改進(jìn)的冪律函數(shù)加高能截?cái)嗪瘮?shù)模型(modified power law cutoff model,MPLCUT)[32]。在此基礎(chǔ)之上,X 射線雙星的能譜還需疊加Fe 的發(fā)射線(能量范圍為6~8 keV)和回旋吸收線(能量范圍為10~80 keV)結(jié)構(gòu),其中回旋吸收線輪廓常用高斯函數(shù)或洛倫茲函數(shù)來擬合。作為例子,圖2展示了平均脈沖相位下4U 0115+634 的X 射線能譜中基頻和五階諧頻回旋吸收線[15,33,34]。

    圖2 4U 0115+634 的X 射線能譜[15]

    3 觀測概論

    觀測上,人們常用以下函數(shù)擬合回旋吸收線:

    其中,Ea,D和W分別是回旋吸收線的中心能量(簡稱能量)、深度和寬度[14]。此外,也可以用高斯函數(shù)擬合回旋吸收線:

    截至目前,人們已在約40 顆X 射線雙星中觀測到了回旋吸收線現(xiàn)象,其中31 顆源的回旋吸收線已被多顆衛(wèi)星的觀測數(shù)據(jù)確認(rèn)。表1列出了這些源的名稱、回旋吸收線能量、中子星自旋周期、雙星軌道周期和參考文獻(xiàn),其中回旋吸收線能量范圍為10~80 keV,中子星自旋周期和雙星軌道周期范圍分別是2~1 100 s 和1~170 d。

    表1 回旋吸收線源及能量等物理參數(shù)

    (續(xù)表)

    從表1可以看出,大多數(shù)源中僅觀測到基頻回旋吸收線,一些源中觀測到了多階諧頻回旋吸收線,如4U 0115+634[135]。近些年,隨著空間X 射線探測技術(shù)的發(fā)展,回旋吸收線觀測方面的研究取得了豐碩的成果。同時(shí),觀測方面的研究進(jìn)展也促進(jìn)了X 射線雙星吸積柱等理論的發(fā)展[16,73,136]。

    3.1 基頻與諧頻回旋吸收線能量的關(guān)系

    人們通過觀測發(fā)現(xiàn),大部分源中僅探測到基頻回旋吸收線,而一些源中則探測到了二階,甚至二階以上諧頻回旋吸收線,且基頻回旋吸收線的輪廓通常比諧頻回旋吸收線的輪廓淺[11]。觀測到二階諧頻回旋吸收線的源包括4U 1907+09[137],4U 1538-52[44,138],Vela X-1[139],Cep X-4[57],Her X-1[140–142]和A 0535+262[143,144];觀測到二階以上回旋吸收線的源包括4U 0115+634 (五階),V 0332+52 (三階)和MAXI J1409-619 (三階),其中4U 0115+634 是唯一一個(gè)觀測到五階諧頻回旋吸收線的源(見圖2)。根據(jù)式(2),二階諧頻回旋吸收線能量與基頻回旋吸收線能量的比值理論上應(yīng)為2。實(shí)際觀測發(fā)現(xiàn),對(duì)于一些源,該比值大于2;而對(duì)于另一些源,該比值小于2。比值大于2 被解釋為多個(gè)回旋吸收線的疊加作用,例如Vela X-1 和A 0535+26[139,143];比值小于2 被解釋為相對(duì)論效應(yīng)修正,例如4U 0115+634[16]。最近人們通過對(duì)Cep X-4 的觀測發(fā)現(xiàn),其二階諧頻回旋吸收線能量與基頻回旋吸收線能量的比值約為1.7。這一觀測結(jié)果被解釋為二階諧頻與基頻回旋吸收線可能形成于不同高度,或者該系統(tǒng)具有較大的觀測角[57]。

    3.2 形態(tài)復(fù)雜性

    人們?cè)鶕?jù)一些模型預(yù)言,回旋吸收線可能存在復(fù)雜的形狀,如基頻回旋吸收線可能存在翼型(wing)形狀[29,30,145]。2010年,RXTE 衛(wèi)星通過對(duì)V 0332+53 的觀測,發(fā)現(xiàn)其回旋吸收線的形狀偏離了高斯函數(shù)或洛倫茲函數(shù)形狀[52]。Doroshenko 等人[146]在2017年用NuSTAR 衛(wèi)星對(duì)該源進(jìn)行重新觀測,卻并沒有觀測到這一現(xiàn)象。他們還指出,回旋吸收線形狀的復(fù)雜性可能源于X 射線能譜擬合模型的不確定性。除了V 0332+53 外,NuSTAR 還觀測到Cep X-4 具有復(fù)雜的回旋吸收線形狀[147]。2017年,Maitra 等人[148]用Suzaku 觀測發(fā)現(xiàn),低光度源Xper 能譜中存在非高斯形狀的回旋吸收線。一般來講,復(fù)雜的回旋吸收線形狀往往在低光度源中被觀測到[11]。Mukherjee 和Bhattacharya[149]認(rèn)為,低光度源中觀測到的復(fù)雜回旋吸收線形狀,可能源于吸積柱局部不穩(wěn)定性導(dǎo)致的偶極磁場扭曲;而Mushtukov等人[150]則認(rèn)為,這種現(xiàn)象可能源于光度變化引起的開普勒增益作用。

    3.3 形態(tài)參數(shù)的相關(guān)性

    回旋吸收線的形態(tài)往往反映了其形成區(qū)的物理環(huán)境。人們研究并發(fā)現(xiàn),回旋吸收線形態(tài)參數(shù)之間存在以下相關(guān)性:(1)回旋吸收線能量與X 射線能譜截?cái)嗄芰康年P(guān)系;(2)回旋吸收線寬度與其能量的關(guān)系;(3)回旋吸收線分?jǐn)?shù)寬度(即回旋吸收線寬度與其能量之比)與其深度的關(guān)系。表2列出了近幾年測量得到的已知源的連續(xù)譜參數(shù)和回旋吸收線形態(tài)參數(shù),包括連續(xù)譜冪律指數(shù)、截?cái)嗄芰浚约盎匦站€的中心能量、寬度和深度。

    表2 回旋吸收線形態(tài)參數(shù)

    早在Ginga 衛(wèi)星時(shí)代,人們通過觀測發(fā)現(xiàn),回旋吸收線能量與X 射線能譜的截?cái)嗄芰看嬖谡嚓P(guān)關(guān)系[14]。這一關(guān)系在RXTE 時(shí)代又被進(jìn)一步證實(shí)[32]。圖3顯示了這一關(guān)系,其中,回旋吸收線能量在35 keV 以下時(shí)相關(guān)性比較明顯,而在高能段相關(guān)性并不很明顯。4U 1626-67 的結(jié)果偏離了該相關(guān)性,這可能是由于其回旋吸收線能量在不同脈沖相位中變化較大,導(dǎo)致其脈沖相位平均譜的回旋吸收線能量的測量不準(zhǔn)確[15]。GX 301-2 和MXB 0656-072的結(jié)果也偏離該相關(guān)性。有人認(rèn)為,這兩顆源中觀測到的回旋吸收線(Ecyc≈50 keV)可能是二階諧頻回旋吸收線,而其基頻回旋吸收線能量約為25 keV。這樣,該基頻回旋吸收線能量就符合與X 射線能譜截?cái)嗄芰康南嚓P(guān)性。也有人認(rèn)為,這兩顆源可能與4U 1626-67 相似,其回旋吸收線能量在不同脈沖相位中變化較大,干擾了其脈沖相位平均譜的回旋吸收線能量的測量。值得注意的是,該相關(guān)性反映了X 射線能譜中特征能量之間的關(guān)系,表明連續(xù)譜的截?cái)嗄芰亢苡锌赡苁且环N磁效應(yīng),而不僅僅由堆積的吸積物質(zhì)的溫度所決定[15]。

    圖3 回旋吸收線能量與連續(xù)譜截?cái)嗄芰康年P(guān)系

    根據(jù)回旋吸收線角度依賴性模型,若回旋吸收線寬度Γcyc由電子溫度kTe決定,則有:

    其中?為觀測方向與電子運(yùn)動(dòng)方向的夾角[154]。式(10)預(yù)言了回旋吸收線能量與寬度之間存在相關(guān)性,且該理論與觀測結(jié)果基本一致(見圖4)。實(shí)際觀測還發(fā)現(xiàn),由不同源推測的電子溫度和觀測角沒有顯著差別。更有趣的是,由于很多回旋吸收線源是高磁傾角系統(tǒng),觀測角的約束間接說明了這些源的磁軸與自轉(zhuǎn)軸幾乎對(duì)齊。對(duì)此現(xiàn)象的一種解釋是,這些源誕生時(shí)就是這樣;另一種解釋是,吸積過程可能使得這些源的磁軸與自旋軸趨于一致。

    有限溫度下的磁場散射截面模型和蒙特卡羅模型等理論預(yù)言,越窄的回旋吸收線應(yīng)該越深[155]。但實(shí)際觀測發(fā)現(xiàn),回旋吸收線分?jǐn)?shù)寬度與其深度存在正相關(guān),即越寬的回旋吸收線往往深度越深(見圖5)。

    圖4 回旋吸收線能量與回旋吸收線寬度的關(guān)系

    圖5 回旋吸收線分?jǐn)?shù)寬度與其深度的關(guān)系

    3.4 回旋吸收線能量的光度依賴性

    早在1976年,Basko 和Sunyaev[26]以愛丁頓光度LEdd為臨界光度,把X 射線雙星分為高光度源和低光度源,并預(yù)言在不同吸積率下,中子星X 射線雙星可能存在兩種不同的吸積模式:對(duì)于高光度源,吸積柱中的下落物質(zhì)受輻射主導(dǎo)的激波作用而制動(dòng),這時(shí)光子將沿垂直于吸積柱的方向逃逸,并形成扇形輻射圖案;對(duì)于低光度源,吸積柱中的物質(zhì)自由下落到中子星表面,并受流體動(dòng)力學(xué)激波作用而制動(dòng),這時(shí)光子將沿平行于吸積柱的方向逃逸,并形成鉛筆形輻射圖案。隨后,以下多種方法被用來模擬回旋吸收線的形成:(1)蒙特卡羅方法[156];(2)Feautrier 方法[154];(3)數(shù)值分析方法[157]。

    根據(jù)該理論預(yù)言,對(duì)于高光度源,隨著吸積率的增加,物質(zhì)的堆積會(huì)使得吸積柱高度增加,從而使回旋吸收線輻射區(qū)逐漸進(jìn)入低磁場區(qū)。這說明吸積率與輻射區(qū)磁場強(qiáng)度可能存在反相關(guān)。通過對(duì)V 0332+53[151],4U 0115+634[158],SMC X-2[55]等回旋吸收線源的幾天到幾個(gè)月的觀測,發(fā)現(xiàn)其X 射線光度與回旋吸收線能量間存在明顯的反相關(guān),從而驗(yàn)證了該理論的預(yù)言(見圖6)。對(duì)于低光度源,通過RXTE 在1996—2005年對(duì)Her X-1 的觀測,發(fā)現(xiàn)其X 射線光度與吸收線能量呈正相關(guān)[140]。類似的源還有GX 304-1[159]和Swift J1626.6-5156[18]。該現(xiàn)象被解釋為,吸積率的增加使得回旋吸收線形成區(qū)的高度降低,磁場強(qiáng)度增大,從而導(dǎo)致X 射線光度與回旋吸收線能量呈正相關(guān)。這一結(jié)果也驗(yàn)證了理論對(duì)不同吸積模式的預(yù)言。

    2012年,Becker 等人[27]在Basko 和Sunyaey 的理論基礎(chǔ)上,考慮光子散射截面和吸積流幾何結(jié)構(gòu)等影響,對(duì)臨界光度進(jìn)行了相應(yīng)修正:

    圖6 回旋吸收線能量與X 射線光度的關(guān)系[9]

    其中,Lcrit是臨界光度,B?為中子星的表面磁場強(qiáng)度。ω與平行于磁場方向傳播的光子的平均能量有關(guān),=ωkTeff,Teff為沖擊后輻射區(qū)域的有效溫度。Λ為參數(shù),Λ<1 代表盤吸積,Λ= 1 代表星風(fēng)吸積[160–162]。對(duì)于低光度情況,即L

    圖7 吸積柱幾何結(jié)構(gòu)和輻射形狀示意圖

    其中,hs是超臨界發(fā)射區(qū)的特征高度,τ?為湯姆孫光深(Thomson optical depth),LX是X射線光度。對(duì)于高光度情況,即下落物質(zhì)受輻射激波作用而制動(dòng),形成扇形輻射圖案,回旋吸收線形成區(qū)高度hs與光度之間呈正相關(guān),如圖7c)和d)所示。

    其中ξ為物質(zhì)下落的有效速度與其流速的比值。

    人們通過近幾年的觀測發(fā)現(xiàn),一些源的回旋吸收線能量與X 射線光度呈不相關(guān)(如A 0535+26[17,74,163]),但在這些源的某些相位中,回旋吸收線能量與X 射線光度呈正相關(guān)。有人認(rèn)為,該現(xiàn)象說明回旋吸收線形成區(qū)的高度可能與吸積率無關(guān),但也有人認(rèn)為,回旋吸收線能量在不同相位中變化較大,這可能導(dǎo)致其能量與X 射線光度的相關(guān)性不明顯[144]。

    3.5 作為幾何追蹤器的脈沖相位解析譜

    由于中子星的自轉(zhuǎn),觀測者在一個(gè)脈沖周期內(nèi)可觀測中子星表面或吸積柱的不同區(qū)域,而回旋吸收線的形成依賴于輻射區(qū)的幾何、磁場結(jié)構(gòu),以及不同發(fā)射區(qū)等離子體的溫度、光深等[11],因此回旋吸收線的觀測性質(zhì)在中子星處于不同脈沖相位時(shí)會(huì)發(fā)生變化。對(duì)該問題的研究需要借助于回旋吸收線的脈沖相位解析譜。目前,人們已對(duì)Her X-1[136,164,165],4U 1538-52[138],4U 0115+634[137],Vela X-1[49,50,139,166],Cen X-3[167],GX 301-2[19],1A 1118-61[168,169],A 0535+26,4U 1907+09,XTE J1946+274[166],4U 1626-67[170]和GX 304-1[55]等源進(jìn)行相關(guān)研究。

    Maitra 和Paul[171]利用Suzaku 數(shù)據(jù)對(duì)一些較亮的回旋吸收線源進(jìn)行觀測,并分析了回旋吸收線參數(shù)隨脈沖相位的變化。首先,他們發(fā)現(xiàn)不同源的回旋吸收線能量會(huì)隨脈沖相位而變化,變化率為10%~30% (見表3);其次,他們分析發(fā)現(xiàn),Vela X-1 的二階諧頻回旋吸收線能量與基頻回旋吸收線能量的比值也會(huì)隨著脈沖相位而變化(見圖8),這說明基頻和諧頻回旋吸收線的形成區(qū)會(huì)隨著觀測角或吸積柱的變化而變化[166];最后,他們發(fā)現(xiàn)4U 1907+09在兩個(gè)不同光度下,回旋吸收線參數(shù)隨著脈沖相位有著相似的變化,這說明在該光度范圍,回旋吸收線形成區(qū)具有相似的發(fā)射高度和幾何結(jié)構(gòu)。理論上,人們一般認(rèn)為回旋吸收線參數(shù)與脈沖相位的依賴性可能與觀測角、發(fā)射區(qū)高度、極冠區(qū)大小、磁場結(jié)構(gòu),以及輻射束形狀(扇形或鉛筆形)等有關(guān)。

    表3 回旋吸收線能量隨脈沖相位的變化[172]

    圖8 a)Vela X-1 回旋吸收線的脈沖相位解析譜;b)不同相位處二階諧頻與基頻回旋吸收線能量之比[166]

    3.6 回旋吸收線能量附近的脈沖輪廓變化

    人們通過觀測發(fā)現(xiàn),一些源的回旋吸收線能量附近的脈沖輪廓發(fā)生了顯著變化(見圖9),這可能是因?yàn)殡娮由⑸浣孛嬖谠撃芰刻幟黠@增強(qiáng),導(dǎo)致脈沖輻射束圖案發(fā)生變化[39,151,173]。對(duì)4U 0115+634 的觀測結(jié)果表明,基頻和諧頻回旋吸收線能量附近的脈沖輪廓均發(fā)生了顯著的變化[39]。對(duì)V 0332+53 的觀測結(jié)果表明,回旋吸收線能量附近脈沖輪廓呈現(xiàn)不對(duì)稱的單峰結(jié)構(gòu),而在其他能段的脈沖輪廓呈現(xiàn)雙峰結(jié)構(gòu)[151]。

    圖9 1A 1118-61,GX 301-2 和Her X-1 的平均脈沖輪廓[166]

    Tsygankov 等人[151]利用INTEGRAL 的觀測數(shù)據(jù),研究了多個(gè)源的脈沖輪廓對(duì)X 射線能量的依賴性,發(fā)現(xiàn)這些源的脈沖輪廓均在回旋吸收線能量附近明顯增大。Maitra和Paul[171]對(duì)GX 301-2,1A 1118-61,Her X-1,XTE J1946+274,A 0535+26,Cen X-3,4U 1907+09 Vela,X-1 和4U 1626-67 等一系列源做了相似的研究,發(fā)現(xiàn)除了源4U 1907+09,Vela X-1 和4U 1626-67 外,其他源在基頻回旋吸收線能量附近的脈沖輪廓均發(fā)生顯著變化[11]。理論上,Sch¨onherr 等人[174]利用數(shù)值分析方法,研究了回旋吸收線能量附近脈沖輪廓的顯著變化,認(rèn)為該現(xiàn)象可能源于電子的回旋共振散射導(dǎo)致的X 射線光子角度的重新分布,且該現(xiàn)象可能受到回旋吸收線形成區(qū)吸積柱幾何的影響。

    3.7 長時(shí)標(biāo)演化

    人們通過觀測發(fā)現(xiàn),一些源的回旋吸收線能量存在長時(shí)標(biāo)演化[11]。Staubert 等人[141]和Soong 等人[164]分析了1996—2012年對(duì)Her X-1 的觀測數(shù)據(jù),發(fā)現(xiàn)其平均脈沖相位譜的回旋吸收線能量呈下降趨勢,其中在2006年之前表現(xiàn)為緩慢衰減,在2006年之后表現(xiàn)為突然衰減(見圖10)。這一結(jié)果被2015年對(duì)NuSTAR 的觀測所證實(shí)[142],并被認(rèn)為該現(xiàn)象可能源于回旋吸收線形成區(qū)的幾何位置變化,或源于持續(xù)吸積過程導(dǎo)致的極冠處磁場結(jié)構(gòu)的變化。Hemphill 等人[175]分析了RXTE 衛(wèi)星在1996—2004年,以及Suzaku 衛(wèi)星在2012年對(duì)源4U 1538-52 的觀測數(shù)據(jù),發(fā)現(xiàn)其平均脈沖相位譜的回旋吸收線能量增加了5%。Mukherjee 和Bhattacharya[149]認(rèn)為,這一現(xiàn)象可能源于吸積柱磁場結(jié)構(gòu)的變化。Cusumano 等人[176]通過觀測發(fā)現(xiàn),在2015年源V 0332+53 爆發(fā)期間,其回旋吸收線能量突然下降,而在爆發(fā)之后,該源的回旋吸收線能量又恢復(fù)到爆發(fā)之前的大小。他們認(rèn)為,這一現(xiàn)象可能是由發(fā)射區(qū)幾何結(jié)構(gòu)變化引起的,而非中子星磁場的衰減導(dǎo)致的[146]。

    圖10 Her X-1 回旋吸收線能量的長時(shí)標(biāo)演化[177]

    4 總結(jié)與展望

    回旋吸收線被發(fā)現(xiàn)已將近40年。隨著空間X 射線探測技術(shù)的提高,回旋吸收線的觀測現(xiàn)象越來越豐富[11]。同時(shí),對(duì)回旋吸收線觀測的進(jìn)步也促進(jìn)了中子星磁場結(jié)構(gòu)、吸積環(huán)境、回旋吸收線發(fā)射機(jī)制等理論的進(jìn)一步完善,其中,通過對(duì)諧頻與基頻回旋吸收線能量的比值的觀測,揭示了發(fā)射區(qū)的局部物理環(huán)境;對(duì)復(fù)雜的回旋吸收線形態(tài)的觀測為探究局部磁場畸變和不穩(wěn)定性提供了線索;對(duì)回旋吸收線形態(tài)參數(shù)間相關(guān)性的觀測,為研究回旋吸收線形成區(qū)等離子體溫度、光學(xué)深度、輻射束形狀和吸積幾何提供了幫助;對(duì)回旋吸收線能量與光度的關(guān)系的觀測,為研究吸積柱幾何隨吸積率變化開辟了一個(gè)新的方向;對(duì)回旋吸收線的脈沖相位解析譜的觀測,有助于分析不同觀測角下回旋吸收線形成區(qū)高度、幾何結(jié)構(gòu)、磁場分布等環(huán)境因素;對(duì)回旋吸收線能量的長期演化的觀測,為探究吸積中子星磁場的長期演化提供了直接證據(jù)。

    未來,回旋吸收線仍是一類熱門的研究課題。首先,NuSTAR 衛(wèi)星等高能譜分辨率X射線衛(wèi)星的研制發(fā)射,有助于進(jìn)一步驗(yàn)證回旋吸收線參數(shù)間的依賴關(guān)系,同時(shí)也有助于深入探測中子星的磁場結(jié)構(gòu)和吸積環(huán)境等。其次,近幾年在GRO J1008-57 中觀測到,其最大基頻回旋吸收線的能量約為76 keV,所對(duì)應(yīng)的磁場強(qiáng)度約為6×108T,接近于磁星磁場的下限[178]。該源的發(fā)現(xiàn),直接證明了強(qiáng)磁場中子星的存在,同時(shí)也為研究吸積磁星與極亮X 射線源(ultra luminous X-ray sources)之間的關(guān)系提供了觀測證據(jù)。此外,對(duì)GRO J1008-57的觀測結(jié)果,還激勵(lì)觀測者進(jìn)一步搜尋高能量回旋吸收線。例如,近年來中國發(fā)射的慧眼衛(wèi)星,因其具有較寬的X 射線能量觀測范圍(1~250 keV),為發(fā)現(xiàn)更多高能量(高于80 keV)的回旋吸收線源和強(qiáng)磁場中子星(不小于109T)提供了觀測上的可能性[21,22]。最后,回旋吸收線方法是直接測量中子星磁場的唯一方法,對(duì)于射電脈沖星,人們常利用偶極輻射假設(shè)估計(jì)中子星磁場。對(duì)比這兩種方法得到的磁場結(jié)果,可以間接佐證射電脈沖星的偶極輻射模型。由于回旋吸收線源中子星的自旋周期大多較長,因此,需要對(duì)長周期射電脈沖星進(jìn)行精細(xì)觀測。目前,中國正在運(yùn)行的FAST 射電望遠(yuǎn)鏡具有極高的射電靈敏度,能夠觀測長周期脈沖星,從而幫助驗(yàn)證射電脈沖星的偶極輻射模型,分析中子星磁場結(jié)構(gòu)[179]。

    猜你喜歡
    中子星基頻能譜
    語音同一認(rèn)定中音段長度對(duì)基頻分析的影響
    能譜CT在術(shù)前預(yù)測胰腺癌淋巴結(jié)轉(zhuǎn)移的價(jià)值
    基于時(shí)域的基頻感知語音分離方法?
    來一勺中子星
    橋面鋪裝層對(duì)中小跨徑橋梁基頻影響分析
    雙中子星合并,又見引力波
    軍事文摘(2017年24期)2018-01-19 03:36:10
    “金箍捧”與“中子星”
    45000kHz基頻晶體濾波器
    電子制作(2017年20期)2017-04-26 06:57:35
    M87的多波段輻射過程及其能譜擬合
    電子材料分析中的能譜干擾峰
    久久精品久久精品一区二区三区| 成人午夜精彩视频在线观看| 色综合欧美亚洲国产小说| 中文字幕人妻熟女乱码| 欧美日韩成人在线一区二区| 亚洲激情五月婷婷啪啪| 美女大奶头黄色视频| 亚洲av日韩精品久久久久久密 | 人人妻人人澡人人爽人人夜夜| 国产成人欧美| 91字幕亚洲| 亚洲成人免费电影在线观看 | 欧美黄色片欧美黄色片| 丝袜美足系列| 男女边吃奶边做爰视频| 捣出白浆h1v1| 啦啦啦在线免费观看视频4| 国产麻豆69| 久久性视频一级片| 午夜福利在线免费观看网站| 欧美老熟妇乱子伦牲交| 考比视频在线观看| 亚洲国产成人一精品久久久| 9热在线视频观看99| 午夜福利视频精品| 午夜久久久在线观看| 国产熟女午夜一区二区三区| 国产精品久久久久久人妻精品电影 | 国产主播在线观看一区二区 | 性色av乱码一区二区三区2| 久久久精品免费免费高清| 777久久人妻少妇嫩草av网站| 校园人妻丝袜中文字幕| 日本wwww免费看| 国产黄色免费在线视频| 老汉色av国产亚洲站长工具| av在线老鸭窝| 国产视频一区二区在线看| 中文字幕高清在线视频| 久久久精品94久久精品| 国产一区二区在线观看av| 男的添女的下面高潮视频| 亚洲,欧美精品.| 美女视频免费永久观看网站| 中文字幕av电影在线播放| 欧美另类一区| 99热国产这里只有精品6| 大陆偷拍与自拍| 欧美老熟妇乱子伦牲交| av在线老鸭窝| 人人妻,人人澡人人爽秒播 | 久久精品aⅴ一区二区三区四区| 男人舔女人的私密视频| 亚洲五月色婷婷综合| 91九色精品人成在线观看| 亚洲av电影在线观看一区二区三区| 亚洲精品国产一区二区精华液| 日韩av在线免费看完整版不卡| 人人妻人人澡人人爽人人夜夜| 日韩欧美一区视频在线观看| 成人国产一区最新在线观看 | 男女下面插进去视频免费观看| 国产精品.久久久| 两个人免费观看高清视频| av在线老鸭窝| 老汉色∧v一级毛片| 亚洲精品成人av观看孕妇| 一本综合久久免费| 久久热在线av| 熟女少妇亚洲综合色aaa.| 菩萨蛮人人尽说江南好唐韦庄| 国产极品粉嫩免费观看在线| 青草久久国产| 99国产精品一区二区蜜桃av | 亚洲第一青青草原| 午夜福利视频在线观看免费| 久久精品国产a三级三级三级| 婷婷成人精品国产| 亚洲一区中文字幕在线| 亚洲欧美清纯卡通| 国产1区2区3区精品| 久久人人爽av亚洲精品天堂| 日本av手机在线免费观看| 成人影院久久| 亚洲国产中文字幕在线视频| 精品久久久久久电影网| 日韩伦理黄色片| 亚洲国产欧美一区二区综合| 大话2 男鬼变身卡| 国产欧美日韩精品亚洲av| 亚洲第一av免费看| av在线app专区| 欧美日韩精品网址| 国产精品熟女久久久久浪| 亚洲国产最新在线播放| 日日摸夜夜添夜夜爱| 久久国产精品人妻蜜桃| 777米奇影视久久| 午夜久久久在线观看| a级毛片在线看网站| 天天添夜夜摸| av有码第一页| 欧美老熟妇乱子伦牲交| 精品国产一区二区三区四区第35| 丝袜喷水一区| 嫩草影视91久久| 国产一级毛片在线| 亚洲一区中文字幕在线| 99热全是精品| 久久久久国产精品人妻一区二区| 在现免费观看毛片| 亚洲一区中文字幕在线| 中文字幕精品免费在线观看视频| netflix在线观看网站| 国产免费视频播放在线视频| 国产成人a∨麻豆精品| 国产成人a∨麻豆精品| 久久精品久久精品一区二区三区| 黑人欧美特级aaaaaa片| 久久热在线av| 国产日韩一区二区三区精品不卡| 午夜免费鲁丝| 中文字幕最新亚洲高清| 青青草视频在线视频观看| 亚洲国产av影院在线观看| 国产精品 国内视频| 久热这里只有精品99| 欧美变态另类bdsm刘玥| 欧美在线黄色| 一区二区日韩欧美中文字幕| av视频免费观看在线观看| 巨乳人妻的诱惑在线观看| 青春草视频在线免费观看| 欧美变态另类bdsm刘玥| 99国产精品免费福利视频| 国产一区二区激情短视频 | 国产淫语在线视频| 天天影视国产精品| 亚洲欧美色中文字幕在线| 国产精品香港三级国产av潘金莲 | 日本vs欧美在线观看视频| 国产成人免费观看mmmm| 国产精品亚洲av一区麻豆| 桃花免费在线播放| 最新在线观看一区二区三区 | 丝袜美足系列| 午夜福利影视在线免费观看| 老司机亚洲免费影院| 老司机在亚洲福利影院| 水蜜桃什么品种好| 亚洲三区欧美一区| 丰满饥渴人妻一区二区三| 每晚都被弄得嗷嗷叫到高潮| 国产深夜福利视频在线观看| 99久久人妻综合| 无遮挡黄片免费观看| 一个人免费看片子| 热99久久久久精品小说推荐| 亚洲国产精品999| 国产精品二区激情视频| 美女高潮到喷水免费观看| 日韩熟女老妇一区二区性免费视频| 欧美黑人欧美精品刺激| 亚洲成色77777| 国产精品久久久久久人妻精品电影 | 成人国产av品久久久| 国产亚洲精品久久久久5区| 久久久亚洲精品成人影院| 超色免费av| 亚洲国产av影院在线观看| 一本—道久久a久久精品蜜桃钙片| 亚洲国产中文字幕在线视频| 两个人看的免费小视频| 久久热在线av| 高潮久久久久久久久久久不卡| 亚洲色图综合在线观看| 中文字幕高清在线视频| 婷婷色综合www| 国产精品久久久久久精品电影小说| 777久久人妻少妇嫩草av网站| 午夜两性在线视频| 久久人人97超碰香蕉20202| av天堂在线播放| 新久久久久国产一级毛片| 欧美在线黄色| av视频免费观看在线观看| 亚洲成av片中文字幕在线观看| 国产在线免费精品| 人体艺术视频欧美日本| 波野结衣二区三区在线| 成人亚洲欧美一区二区av| 交换朋友夫妻互换小说| 亚洲精品日韩在线中文字幕| av又黄又爽大尺度在线免费看| 免费在线观看日本一区| 免费一级毛片在线播放高清视频 | 久久国产亚洲av麻豆专区| 美国免费a级毛片| 亚洲一区中文字幕在线| 久久久久久免费高清国产稀缺| 老汉色∧v一级毛片| 亚洲 国产 在线| 欧美在线黄色| 欧美激情极品国产一区二区三区| 天堂中文最新版在线下载| 男女边摸边吃奶| 亚洲色图 男人天堂 中文字幕| 极品人妻少妇av视频| 成人三级做爰电影| 人人妻人人添人人爽欧美一区卜| 热re99久久国产66热| 亚洲欧美一区二区三区国产| 亚洲精品中文字幕在线视频| 国产精品麻豆人妻色哟哟久久| 久热这里只有精品99| 午夜福利,免费看| 日本a在线网址| 免费高清在线观看视频在线观看| 国产成人欧美在线观看 | 欧美少妇被猛烈插入视频| 看免费av毛片| 午夜福利影视在线免费观看| 亚洲欧美激情在线| 亚洲欧美色中文字幕在线| 日本黄色日本黄色录像| 亚洲欧美一区二区三区黑人| 久久精品亚洲av国产电影网| 日本vs欧美在线观看视频| 久久影院123| 咕卡用的链子| 一级黄色大片毛片| 久久人人97超碰香蕉20202| 一区二区日韩欧美中文字幕| 国产一区二区激情短视频 | 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产欧美一区二区综合| 18禁黄网站禁片午夜丰满| 免费av中文字幕在线| 咕卡用的链子| 在线观看www视频免费| 丝袜喷水一区| 一本一本久久a久久精品综合妖精| 岛国毛片在线播放| 天堂中文最新版在线下载| 69精品国产乱码久久久| 在线av久久热| 国产片内射在线| 精品卡一卡二卡四卡免费| 国产视频一区二区在线看| 日韩制服丝袜自拍偷拍| 80岁老熟妇乱子伦牲交| 国产精品一区二区在线观看99| 亚洲av欧美aⅴ国产| 极品人妻少妇av视频| 久久久精品94久久精品| 9191精品国产免费久久| 国产黄色视频一区二区在线观看| 欧美日韩亚洲高清精品| 久久毛片免费看一区二区三区| 老司机影院毛片| 满18在线观看网站| 一区二区三区精品91| 中文字幕色久视频| xxxhd国产人妻xxx| 99国产精品一区二区三区| 亚洲伊人色综图| 精品人妻在线不人妻| 国产一区二区三区综合在线观看| 亚洲精品一卡2卡三卡4卡5卡 | 黄片播放在线免费| 五月天丁香电影| 狠狠精品人妻久久久久久综合| 久久毛片免费看一区二区三区| 成人亚洲精品一区在线观看| 婷婷色麻豆天堂久久| 99久久精品国产亚洲精品| 高潮久久久久久久久久久不卡| 精品少妇一区二区三区视频日本电影| 国产在线一区二区三区精| 宅男免费午夜| 99国产综合亚洲精品| 欧美成狂野欧美在线观看| 欧美日韩亚洲国产一区二区在线观看 | 欧美+亚洲+日韩+国产| 国产不卡av网站在线观看| 99精国产麻豆久久婷婷| 国产一区二区三区av在线| 一级片'在线观看视频| 国产在视频线精品| 晚上一个人看的免费电影| 久久久精品免费免费高清| 国产精品av久久久久免费| 大香蕉久久成人网| av又黄又爽大尺度在线免费看| 9191精品国产免费久久| 午夜精品国产一区二区电影| 大陆偷拍与自拍| 久9热在线精品视频| 中国国产av一级| 在现免费观看毛片| 国产深夜福利视频在线观看| 欧美黄色片欧美黄色片| 男的添女的下面高潮视频| 免费不卡黄色视频| 精品第一国产精品| 黄色视频在线播放观看不卡| 极品少妇高潮喷水抽搐| 亚洲国产日韩一区二区| 麻豆国产av国片精品| www.熟女人妻精品国产| 久久精品亚洲熟妇少妇任你| 韩国精品一区二区三区| 五月开心婷婷网| 亚洲色图 男人天堂 中文字幕| 一级毛片电影观看| 亚洲自偷自拍图片 自拍| 成年av动漫网址| 中文字幕人妻丝袜一区二区| 欧美精品人与动牲交sv欧美| 七月丁香在线播放| 热re99久久国产66热| 999久久久国产精品视频| 久久精品国产亚洲av涩爱| 一级,二级,三级黄色视频| 精品久久蜜臀av无| 精品卡一卡二卡四卡免费| 亚洲精品国产av蜜桃| 国产一区二区三区综合在线观看| 欧美成狂野欧美在线观看| 一级片'在线观看视频| 精品卡一卡二卡四卡免费| 国产欧美日韩一区二区三区在线| avwww免费| 亚洲欧美一区二区三区黑人| 国产精品亚洲av一区麻豆| 精品久久久精品久久久| av在线app专区| 免费女性裸体啪啪无遮挡网站| 日韩一卡2卡3卡4卡2021年| 热re99久久国产66热| 老汉色∧v一级毛片| 国产av一区二区精品久久| 精品福利永久在线观看| 爱豆传媒免费全集在线观看| 午夜91福利影院| 丝瓜视频免费看黄片| 成人免费观看视频高清| 99热国产这里只有精品6| 久久精品久久精品一区二区三区| 国产精品一二三区在线看| 尾随美女入室| 日韩人妻精品一区2区三区| 十分钟在线观看高清视频www| 91精品伊人久久大香线蕉| 亚洲专区中文字幕在线| 久久精品国产亚洲av涩爱| 亚洲五月婷婷丁香| √禁漫天堂资源中文www| 午夜免费鲁丝| 午夜两性在线视频| 久久久欧美国产精品| 国产亚洲欧美在线一区二区| 少妇裸体淫交视频免费看高清 | 精品国产一区二区三区四区第35| 香蕉丝袜av| 国产极品粉嫩免费观看在线| 久久青草综合色| 亚洲欧美精品综合一区二区三区| 日韩视频在线欧美| 国产免费福利视频在线观看| 久久久久久亚洲精品国产蜜桃av| 欧美 亚洲 国产 日韩一| 精品人妻一区二区三区麻豆| 女人久久www免费人成看片| 在线观看免费午夜福利视频| 久久精品久久精品一区二区三区| 蜜桃国产av成人99| 亚洲成av片中文字幕在线观看| 国产黄频视频在线观看| 国产人伦9x9x在线观看| 亚洲精品第二区| 国产91精品成人一区二区三区 | 亚洲中文字幕日韩| 国产91精品成人一区二区三区 | 亚洲一区中文字幕在线| 日韩电影二区| 亚洲av欧美aⅴ国产| 国产片内射在线| 亚洲三区欧美一区| 欧美黄色片欧美黄色片| 久久99精品国语久久久| 久久久精品94久久精品| 777久久人妻少妇嫩草av网站| 日韩精品免费视频一区二区三区| 亚洲一码二码三码区别大吗| 国产激情久久老熟女| 亚洲精品一卡2卡三卡4卡5卡 | 丝袜脚勾引网站| 久久毛片免费看一区二区三区| 国产成人精品无人区| 国产成人啪精品午夜网站| 亚洲精品一卡2卡三卡4卡5卡 | 黄色视频在线播放观看不卡| a级片在线免费高清观看视频| 十八禁人妻一区二区| 国产成人av激情在线播放| 美女视频免费永久观看网站| 女性生殖器流出的白浆| e午夜精品久久久久久久| www日本在线高清视频| 午夜影院在线不卡| 美女大奶头黄色视频| 免费在线观看日本一区| 美女中出高潮动态图| 精品熟女少妇八av免费久了| 精品卡一卡二卡四卡免费| 国产精品久久久av美女十八| 老鸭窝网址在线观看| 亚洲欧洲精品一区二区精品久久久| 韩国高清视频一区二区三区| 日本黄色日本黄色录像| 亚洲精品国产av蜜桃| 制服诱惑二区| 亚洲精品国产av成人精品| 80岁老熟妇乱子伦牲交| 亚洲国产av新网站| 少妇 在线观看| 女人久久www免费人成看片| 午夜久久久在线观看| www.自偷自拍.com| 国产在线视频一区二区| 最新在线观看一区二区三区 | 乱人伦中国视频| 少妇精品久久久久久久| a级毛片在线看网站| 午夜免费鲁丝| 国产日韩欧美亚洲二区| 日韩av不卡免费在线播放| 人人澡人人妻人| 精品人妻1区二区| 亚洲七黄色美女视频| 国产高清国产精品国产三级| 纯流量卡能插随身wifi吗| 国产免费福利视频在线观看| 国产精品一区二区在线不卡| 久久久久久久大尺度免费视频| 国产av精品麻豆| 成人影院久久| av天堂在线播放| 免费观看a级毛片全部| 国产精品99久久99久久久不卡| 成人亚洲欧美一区二区av| 国产99久久九九免费精品| 丁香六月天网| 亚洲精品久久久久久婷婷小说| 真人做人爱边吃奶动态| 国产精品香港三级国产av潘金莲 | 99九九在线精品视频| 国产av一区二区精品久久| 婷婷色综合大香蕉| 欧美精品啪啪一区二区三区 | 欧美日韩黄片免| 亚洲av电影在线进入| 亚洲免费av在线视频| 蜜桃国产av成人99| 一级片免费观看大全| 女人爽到高潮嗷嗷叫在线视频| 国产成人系列免费观看| 久久久亚洲精品成人影院| 脱女人内裤的视频| 久久国产精品男人的天堂亚洲| 又大又爽又粗| 精品国产国语对白av| 午夜福利一区二区在线看| 一区二区av电影网| 侵犯人妻中文字幕一二三四区| 久久精品人人爽人人爽视色| 夫妻午夜视频| 黑人巨大精品欧美一区二区蜜桃| 日本午夜av视频| 水蜜桃什么品种好| 亚洲人成77777在线视频| 在线观看国产h片| 精品人妻一区二区三区麻豆| 亚洲av国产av综合av卡| 亚洲精品第二区| 男女下面插进去视频免费观看| av不卡在线播放| 丝袜美腿诱惑在线| 精品卡一卡二卡四卡免费| 色94色欧美一区二区| 国产免费视频播放在线视频| 国产成人精品在线电影| 国产成人精品久久二区二区免费| 久久国产精品大桥未久av| 午夜日韩欧美国产| 色网站视频免费| 王馨瑶露胸无遮挡在线观看| 久久精品人人爽人人爽视色| 久久综合国产亚洲精品| 赤兔流量卡办理| www.精华液| 久久亚洲精品不卡| 婷婷丁香在线五月| 一个人免费看片子| 性少妇av在线| 韩国精品一区二区三区| 欧美人与性动交α欧美精品济南到| 天堂俺去俺来也www色官网| 日韩av在线免费看完整版不卡| 国产精品熟女久久久久浪| 你懂的网址亚洲精品在线观看| 一本综合久久免费| 国产精品一区二区在线观看99| 中文字幕av电影在线播放| 亚洲国产日韩一区二区| 欧美另类一区| 国产亚洲av片在线观看秒播厂| 老鸭窝网址在线观看| 精品一区在线观看国产| 飞空精品影院首页| 国产伦人伦偷精品视频| 精品一区二区三区av网在线观看 | 一本一本久久a久久精品综合妖精| av片东京热男人的天堂| 天堂俺去俺来也www色官网| 在线看a的网站| 中文字幕高清在线视频| 欧美黄色片欧美黄色片| 国产精品人妻久久久影院| 91成人精品电影| 两人在一起打扑克的视频| 国产日韩欧美在线精品| av网站在线播放免费| 久久狼人影院| 黑人猛操日本美女一级片| 国产精品久久久人人做人人爽| 国产一区亚洲一区在线观看| 日本五十路高清| 亚洲免费av在线视频| 极品少妇高潮喷水抽搐| 在线av久久热| 一边摸一边做爽爽视频免费| 麻豆av在线久日| 一边摸一边做爽爽视频免费| 十八禁人妻一区二区| 午夜两性在线视频| 久久九九热精品免费| 亚洲精品中文字幕在线视频| 在线观看免费午夜福利视频| 亚洲精品国产av蜜桃| 视频区图区小说| 亚洲 国产 在线| 国产成人一区二区三区免费视频网站 | 亚洲欧美清纯卡通| 欧美老熟妇乱子伦牲交| 色视频在线一区二区三区| 19禁男女啪啪无遮挡网站| 天堂8中文在线网| 高清视频免费观看一区二区| 国产av国产精品国产| www.av在线官网国产| 男男h啪啪无遮挡| 性高湖久久久久久久久免费观看| netflix在线观看网站| 色网站视频免费| 91九色精品人成在线观看| 叶爱在线成人免费视频播放| 精品国产国语对白av| 亚洲人成77777在线视频| 国产免费视频播放在线视频| 美女高潮到喷水免费观看| av在线老鸭窝| 久久久久国产一级毛片高清牌| 一个人免费看片子| 中文精品一卡2卡3卡4更新| 欧美变态另类bdsm刘玥| 丰满迷人的少妇在线观看| 日韩,欧美,国产一区二区三区| 在线观看免费高清a一片| 99国产综合亚洲精品| 亚洲自偷自拍图片 自拍| 国产亚洲精品久久久久5区| av不卡在线播放| 国产精品久久久久久人妻精品电影 | 亚洲图色成人| 国产成人91sexporn| 午夜免费成人在线视频| 一本—道久久a久久精品蜜桃钙片| 精品亚洲成a人片在线观看| av片东京热男人的天堂| 国产野战对白在线观看| 成年女人毛片免费观看观看9 | 久久久久久人人人人人| 国产精品秋霞免费鲁丝片| 午夜福利视频在线观看免费| 少妇人妻久久综合中文| 国产深夜福利视频在线观看| 日本vs欧美在线观看视频| 免费少妇av软件| 精品欧美一区二区三区在线| 波多野结衣av一区二区av| 一边亲一边摸免费视频| 热99久久久久精品小说推荐| 国产精品 欧美亚洲| 午夜福利乱码中文字幕| 两个人看的免费小视频| 免费在线观看完整版高清| 亚洲国产精品一区三区| 18禁国产床啪视频网站| 日本av手机在线免费观看| 91精品三级在线观看| 日本av手机在线免费观看| 国产在线免费精品| 久久精品国产综合久久久|