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

    H2的中紅外波段研究進展

    2019-04-10 02:00:24周葉艷
    天文學進展 2019年1期
    關鍵詞:質量

    周葉艷,郝 蕾

    (1.中國科學院 上海天文臺,星系與宇宙學重點實驗室,上海200030; 2.中國科學院大學,北京100049)

    1 引 言

    星際介質是星系中重要的組成成分,它在恒星形成過程中起著至關重要的作用,恒星形成就發(fā)生在星際介質中[1]。一個星系通常由恒星、星際氣體和塵埃組成。氣體和塵埃組成了彌散的星際介質(其中氣體約占99%)。星際介質中的元素,70% 由H 組成,剩余的大多為He。在星際介質中,致密的物質形成氣體云[2],在這些致密的氣體云中,H 大多以分子的形態(tài)存在,所以這些形成恒星的氣體云又稱為分子云[3]。這些分子云就是恒星誕生的搖籃[3]。根據位力定理[4],當氣體分子云的引力勢能等于內部熱能的兩倍時,氣體壓力產生的熱能與引力勢能達到平衡,此時,該氣體分子云可以保持流體靜力學平衡。反之,當氣體分子云質量足夠大,乃至于氣體壓強無法抵抗引力勢能,這時,分子云開始坍縮,最終經過一系列復雜的過程,形成恒星。而恒星又生存于星際介質中,一些物質可能以星風的形式被噴射出去,而當恒星消亡之后,部分元素和能量又回到星際介質中。因此,星際介質為恒星形成、演化和結束提供場所。對星際介質的研究對于我們研究星系演化具有重要意義。

    星際介質包含從離子到分子等各種狀態(tài)的粒子和塵埃,H2是宇宙中含量最豐富的分子。從地球大氣到類星體的宿主星系,H2存在于不同種類的環(huán)境之中。同時,它也是星系中分子云的最主要成分[5]。一方面,恒星形成于冷的分子云的收縮,H2可以被視作恒星形成主要的原初燃料。Kennicutt[6]和Schmidt[7]發(fā)現恒星形成率與中性氣體(包括原子H 氣體和H2氣體的總和)之間存在緊密的關系(簡稱KS 關系),但后續(xù)人們逐漸認識到恒星形成與分子氣體更直接相關。在KS 關系中,如果使用分子氣體表面密度,則可得斜率為1,即分子云在不同星系中的恒星形成率基本相同。另一方面,在恒星形成過程中,H2作為冷卻劑也推動了恒星的形成[5]。有研究表明,在宇宙早期,H2在原初氣體冷卻的過程中起主導作用,加快了氣體的冷卻[8]。為了研究恒星形成過程以及H2在整個宇宙恒星形成歷史中如何變化和參與,準確測定H2的質量和分布是非常必要的。

    雖然H2的質量在富氣體星系中能達到1010M⊙[9],但是在冷星際介質中,H2不能有效地進行輻射,因此難以被直接觀測。主要原因有:(1)H2是同核雙原子分子,不具有永久性的電偶極矩和磁偶極矩,因此電子基態(tài)的電、磁偶極矩的振動-旋轉躍遷是禁戒的,只發(fā)生極微弱的電四級矩躍遷;(2)H2的振動轉動躍遷激發(fā)溫度較高,振動基態(tài)的最低能級轉動躍遷的等效溫度為510 K[8,10],冷氣體分子云(T <100 K)很難將分子激發(fā)到H2的轉動激發(fā)態(tài)[11]。基于以上原因,人們通過觀測CO 分子的旋轉發(fā)射線等間接的方式來研究分子氣體的質量[12–15]。由于CO 分子是星際氣體分子中除了H2以外含量最豐富的分子,并且CO 分子譜線較容易被觀測到(CO 分子具有永久性的偶極矩),CO 分子的豐度會隨著金屬豐度的降低而極速下降[16,17],所以在低金屬豐度的環(huán)境下,利用CO 分子間接獲得H2質量的方法將失去可靠性;另外,這種間接的追蹤行為都依賴于一定的物理環(huán)境,也就是說,這些間接得到H2質量的方法并不適用于所有環(huán)境。在一些特定的星際環(huán)境中測得的結果會存在較大偏差。因此,找到不受環(huán)境因素影響并直接測量H2質量的方法極為必要。

    近年來的研究表明,在溫度較高的介質中,人們可觀測到H2輻射。這些輻射主要集中在紅外波段(這些輻射都來自于電四級矩躍遷)。Gautier 等人在1976年首次在獵戶星云中觀測到H2的振動-旋轉譜線[18],這一譜線(2.12 μm)被證明是在溫度范圍為2 000~3 000 K 的條件下產生的。隨著大型紅外望遠鏡的投入使用,人們能夠在中紅外波段探測到越來越多H2振動-旋轉譜線。H2在中紅外波段的旋轉譜線波長分布在5~30 μm 范圍內,溫度范圍在100~1 000 K 之間(如表1 所示)。表1 列出常見的中紅外波段H2旋轉譜線的波長,對應的轉換躍遷在單位溫度下的能級上限,以及它們的輻射率系數,對應的H2譜線的名稱縮寫為S(0)—S(7)。現在我們常把溫度范圍在20~50 K 內,很難直接探測而只能通過CO 分子來間接研究的H2稱為“冷”的H2,而對于溫度在50~100 K 以上,到幾百開爾文溫度的,能在紅外波段探測到的H2稱為“溫”的H2。這些“溫”的H2輻射能否跟蹤總分子氣體,以及與恒星形成是否相關,已逐漸成為目前H2研究領域的熱點。

    已有學者針對這些H2輻射是否直接與恒星形成區(qū)的冷分子氣體相關這一問題給出不同的結果。例如,Roussel 等人[22]于2007年證實H2輻射與已被公認能夠表征恒星形成的多環(huán)芳烴(polycyclic aromatic hydrocarbon, PAH)輻射相關;Zakamska 等人[23]則于2010年發(fā)現,中紅外波段探測的H2輻射很可能來自于塵埃區(qū)之外,因此與恒星形成活動不直接相關,而更可能由激波引起;2016年,Togi 和Smith[19]又發(fā)現,通過對H2輻射反映出來的H2溫度建立模型,可以外推整體的H2質量。本文通過比較這些工作,討論其中是否存在矛盾之處。在第2 章,我們利用低光度星系樣本和極亮紅外星系樣本,著重討論中紅外H2發(fā)射強度與恒星形成的關系;第3 章介紹H2輻射與總分子氣體質量之間的關系,還根據H2在中紅外波段的譜線強度推算H2總質量的方法,并進行討論。

    表1 觀測得到的H2譜線[19]

    2 中紅外波段H2發(fā)射強度與恒星形成的關系

    星系中的恒星形成與星系氣體含量,特別是與致密分子氣體含量之間存在緊密的關系[6,7](如KS 關系)。那么,中紅外波段探測到的“溫”或“熱”H2譜線與星系的恒星形成有相關性嗎?

    最初人們認為,觀測到的H2振動-旋轉輻射可能來源于年輕恒星周圍的光致電離區(qū)域[24,25]。在這種情況下,主要由于新形成的大質量恒星中的紫外光子激發(fā)H2,促使H2直接轉換到激發(fā)電離態(tài),或是紫外光子先加熱塵埃,而后塵埃加熱氣體。Roussel 等人[22]在2007年研究低光度恒星形成星系時也證實了這一觀點。然而,Appleton 等人[26]用多波段SPITZER IRS 光譜發(fā)現,在無塵埃的河外天體史蒂芬五重星系(Stephan’s Quintet, SQ)中存在罕見的高光度旋轉H2輻射,而在該天體中并未發(fā)現有很強的恒星活動。爾后,相似的觀測結果被多人證實[27–29],在這些天體中,少數甚至根本沒有恒星形成區(qū)域。這一實驗事實說明H2譜線的激發(fā)可能存在很多種不同的物理機制。事實上,H2的發(fā)射線在各種不同的環(huán)境中都可探測到,如年輕恒星的外流[18,30]、光解離區(qū)域[31,32]、行星狀星云[33,34]、超新星遺跡[35,36]、星系中心的大區(qū)域[37,38],以及來源于星系碰撞時產生的激波[2,39,40]。因此,所處的環(huán)境不同,H2譜線的激發(fā)機制也可能不相同。目前,已經提出的H2激發(fā)的主要物理機制包括:(1)與年輕恒星周圍的恒星形成區(qū)域有關[24,25],(2)與大尺度的激波有關。大尺度激波來源于相對論性的噴射[41]或者星系并合[26,28],又或者是星系團間氣體掉落至中心星系[27];(3)在星系際介質的纖維狀結構里,H2的激發(fā)輻射可能來自于宇宙線堆積凝結成分子氣體[42,43];(4)星系團中心星系中的H2可能來自于整個星系團熱氣體云的電子加熱[44]。

    在不同的物理機制下,中紅外波段發(fā)射的H2譜線是否可以用來追蹤恒星形成,人們所持觀點也各不相同。目前主要幾個觀點為:(1)Roussel 等人[22]于2007年通過研究發(fā)現,對于低光度(相比于極亮紅外星系而言)的星系,在恒星形成區(qū)域占主導的這些星系中,“溫”H2譜線是恒星形成活動的產物;而在低電離星系核(Low-ionization nuclear emission-line region, LINER)和賽弗特(Seyfert)天體中,“溫”H2譜線是激波的示蹤者。(2)對于極亮紅外星系(Ultra-luminous infrared galaxy, ULIRG,其紅外光度通常更高),Higdon 等人[5]認為中紅外波段觀測到的H2旋轉譜線不受塵埃影響,Zakamska 等人[23]證實了這一結果,并深入研究后發(fā)現,中紅外波段的H2譜線與恒星形成不直接相關而與激波有關[44]。

    2.1 低光度星系的H2研究

    在H2研究初期,紅外光譜技術還無法提供靈敏的較低光度星系的觀測數據,大部分觀測都只關注極亮的星系,并不具有普適意義。因此,Roussel 等人[22]在2007年利用SINGS樣本[45](Spitzer Infrared Nearby Galaxies Survey)中部分正常星系的數據(包括一組較寬區(qū)間范圍紅外光度的星系,且包括不同種類的星系),并加入其他相對較暗的天體樣本,期望在擴大研究樣本的條件下得到一條適合更多普通星系的H2譜線規(guī)律。

    我們選取了SINGS 項目中的57 個星系作為樣本。樣本中的每個星系都至少含有4 條H2S(0)—S(3)的最低能量旋轉譜線觀測數據,且其中有極少數星系可觀測到能級S(7)。該樣本包括大質量星系、矮星系、低電離星系核、賽弗特星系,以及恒星形成星系和漩渦星系等多種類型。因此,得到的結果更具有普遍適用的意義。樣本中,17 μm S(1)線在恒星質量超過109.5M⊙的星系中均有探測到。譜線測量的覆蓋區(qū)域平均達到0.9 kpc,覆蓋了大部分的恒星形成區(qū)。

    文章推導出的單個星系中的“溫”H2(占據1%~30% 的總H2質量),在千秒差距尺度上的柱密度平均值與單個恒星形成區(qū)的光解離區(qū)域的值一致。而作為恒星形成追蹤者(被星際塵埃所遮蔽)的PAH 發(fā)射線,被認為主要來源于光解離區(qū)域。Roussel 等人的主要結論如圖1 a)所示。在恒星形成星系樣本中,H2譜線強度與PAH 發(fā)射線強度有緊密的相關性,S(0)—S(2)的總譜線強度積分與PAH 發(fā)射線強度的比值不隨H2總譜線強度的增加而增加,而是在一個平均值上下浮動,這表明兩種譜線的來源相同。由于PAH 發(fā)射線已經被證實來源于恒星形成區(qū),所以H2在恒星形成星系樣本中的激發(fā)來源也可能是由恒星形成引起的。而樣本中的25 個低電離星系核/賽弗特星系與恒星形成星系樣本則有著不同的結果,與樣本中的其他星系相比,它們有更多更“溫”的H2旋轉譜線強度。我們假設LH2/LPAH的比值存在一個閾值,這個值可以作為識別低光度活動星系核(active galactic nucleus, AGN)的標志。對于這部分特殊的星系,激波加熱可以用來解釋它們的特殊結果。

    根據Roussel 等人的結論可知,大部分的恒星在星系中形成,H2的激發(fā)機制由恒星形成區(qū)引起。而在活躍的活動星系核中,H2的激發(fā)由激波的加熱導致。

    圖1 恒星形成星系和活動星系核中H2與PAH 發(fā)射線的關系[22]

    2.2 ULIRG的H2研究

    ULIRG 是鄰近宇宙中整體光度最大的天體,這些天體被認為是由富氣體塵埃的盤星系合并而成[23,46]。Higdon 等人在2006年利用Spitzer 紅外望遠鏡的數據研究ULIRG 時發(fā)現,由9.7μm 硅酸鹽吸收線光深得到的高消光在觀測H2的視線中并未探測到,從而推測中紅外波段觀測到的H2旋轉譜線不受塵埃的影響[5]。與Roussel 等人[22]的研究類似,Zakamska等人[23]同樣用PAH 譜線作為恒星形成的示蹤者,于2010年也驗證了這一結論。如果PAH與H2的輻射都由恒星形成激發(fā),那么這些譜線都應該受到塵埃的消光影響。由于消光強度在很大程度上取決于波長,不同譜線受消光的影響程度也會不同(如圖2)。而落在9.7μm 硅酸鹽吸收線的不透明度峰值之內的譜線應該會受消光影響最嚴重?;谶@個事實,人們利用吸收譜的峰之內與之外的譜線流量的比值,就能夠測試出譜線激發(fā)區(qū)域相對于不透明度源的位置。通過判斷H2的輻射與不透明度的輻射區(qū)域是否一致,可以判斷H2所在的區(qū)域。即使ULIRG 在中紅外波段也是光學厚的[23],這就意味著這種不透明度的影響效應會非常強,所以對于ULIRG 而言,這種判斷H2所在區(qū)域的方法是可行的。

    Zakamska 等人[23]研究發(fā)現,PAH 譜線的比值與吸收線的強度相關(見圖3),這也證實了PAH 受塵埃的影響,即PAH 的發(fā)射區(qū)域是在硅酸鹽與水冰的位置之間;與之相反,H2的流量和比值與硅酸鹽不透明度之間沒有顯示出任何的相關關系。即使是離不透明度峰值最近的S(3),本應該被嚴重影響,也沒有顯示出H2的所在區(qū)域,即H2不受吸收線的影響。因此,他們得到結論:恒星形成區(qū)域隱藏在光學厚的氣體塵埃云內,而H2的分子旋轉譜線強度是基于恒星形成率而得到值的3 倍多。因此,H2與恒星形成區(qū)不在同一空間位置上,而是位于塵埃區(qū)以外。故塵埃消光影響了恒星形成示蹤者PAH 但不影響H2旋轉譜線,這一結果與傳統上認為H2譜線與恒星形成活動直接相關的結論相矛盾。

    在2014年,Hill 和Zakamska 等人[44]又進一步研究了極亮紅外星系中旋轉H2譜線的起源。在這些星系中,他們發(fā)現溫的H2都與[Fe II]譜線顯著相關(見圖4)。而在之前的研究中,強[Fe II]發(fā)射線通常會與激波有很強相關性。在激波破壞塵埃粒子的過程中,Fe 元素會變?yōu)橹行越橘|,在碰撞激發(fā)下進一步變成[Fe II][47]。同樣的觀點還被H2與激波的光學譜線比值之間的關系所證實。由于[O I]λ6 300/Hα 已經被證明與激波有很強的相關性,而如圖5 所示,隨著這一光學譜線比值的增加(激波的貢獻比重也對應增大),LH2/LPAH也隨之增大。這一結果也同樣證實H2與激波的相關性。恒星形成星系和活動星系核主導星系里都有類似H2與[Fe II]的關系。這一超強的H2發(fā)射現象被解釋為來自超新星和活動星系核外流進入ULIRGs 的中性星際介質激發(fā)。

    這些天體中的H2譜線可以用來追蹤周圍環(huán)境中的那些依次被星系與鄰近星系相互作用而激發(fā)的激波,那些大規(guī)模的強激波通過H2輻射冷卻的現象可能比之前更加普遍。在早期的宇宙中,這一過程中H2輻射的增加可以加速物質的冷卻,最后它收縮形成第一代恒星和星系,因此超強的H2輻射有助于第一代恒星和星系的形成。

    圖2 ULIRG 光譜和代表性的不透明度曲線中的特征譜線波長[22]

    圖3 PAH 特征譜線與塵埃消光的相關性[23]

    圖4 FS(3) 和F[FeII]λ 5.34 μm 的光度關系[44]

    圖5 LH2/LPAH 與光學譜線強度比值之間的關系[44]

    2.3 討論

    根據上述針對低光度星系和極亮紅外天體中H2在中紅外波段的輻射是否由恒星形成激發(fā)的研究,人們得到不盡相同的結論。Roussel 等人[22]認為,對于低光度的正常星系,在電離氫區(qū)占主導的這些星系中,溫H2譜線是恒星形成過程的產物;而在低電離核發(fā)射星系LINERs 和賽弗特星系中,溫H2譜線是激波的示蹤者;而Zakamska 等人[23]認為,在ULIRG 中,中紅外波段的H2譜線與激波有關,而與恒星形成不直接相關。上述兩個工作選取的樣本不同,使用的方法也不完全相同。為了理解他們之間的關系,檢驗這兩個工作是否一致,我們嘗試利用Zakamska 的樣本去重復Roussel 等人的研究。

    在Roussel 等人的工作中,他們利用PAH 發(fā)射線作為恒星形成的示蹤者,且用S(0)—S(2)譜線強度的和代表H2。由于Zakamska 等人的樣本中沒有S(0)的數據,因此我們利用數據較多的S(3)以及S(1)—S(3)總光度的和來代替Roussel 等人工作中的H2譜線強度和,試圖尋找Zakamska 等人的樣本星系中H2輻射與恒星形成的關系。結果如圖6 所示,在Zakamska 等人的樣本中,H2譜線與PAH 發(fā)射線之間并沒有強烈的相關性(隨著H2譜線強度的增加,H2譜線光度與PAH 發(fā)射線光度的比值并不是穩(wěn)定不變的)。

    為進一步研究H2譜線與PAH 發(fā)射線之間的關系,我們給出S(3)以及S(1)—S(3)的總光度和與PAH 發(fā)射線光度之間的關系(見圖7)。與圖6 類似,H2譜線與PAH 發(fā)射線光度之間沒有明顯的相關性,此結果與Roussel 等人的賽弗特星系/LINERs 樣本結果類似,也與Zakamska 等人得出的H2譜線與恒星形成不直接相關而與激波有關的結論一致。

    因為ULIRG 作為一個類別,可能包括活動星系核,也可能包括恒星形成星系。因此,我們進一步探索Zakamska 等人的樣本中每一個星系的類型。從圖7 可以看出,在Zakamska 等人的樣本中,星系類型大多為AGN 和復合星系,極少數為恒星形成的星系。因此Zakamska等人的樣本更接近Roussel 等人的AGN 樣本。雖然這兩個工作針對不同光度的樣本,但基本結論是一致的。

    綜合以上研究結果,我們可以發(fā)現:在恒星形成的星系中,中紅外波段的H2譜線主要由恒星形成激發(fā),因此這些H2譜線應該可以成為恒星形成的示蹤者;而在非恒星形成主導的星系中,H2譜線激發(fā)極有可能來源于激波。

    圖6 Zakamska 樣本中H2譜線與PAH 發(fā)射線強度之間的關系

    圖7 H2譜線與PAH 發(fā)射線強度之間的關系

    3 中紅外波段H2譜線與總分子氣體質量的研究

    長期以來,人們認為由中紅外波段的H2輻射反映的溫熱分子氣體只占總分子氣體的一小部分,而與恒星形成直接相關的冷分子氣體占總分子氣體的大多數[5]。因此,中紅外波段觀測的H2輻射并不能很好地反映總分子氣體質量。同時,因為H2在冷分子氣體環(huán)境中不能進行有效的輻射,在很長的時間里,我們無法通過H2來獲得總分子氣體的質量。人們通常先使用CO 的分子旋轉譜線得到CO 的分子質量,然后再推算出星際介質中分子氣體的總質量[12,13]。然而,近年來,隨著越來越多的星際氣體被觀測到,這一觀點正在逐漸受到挑戰(zhàn)。我們知道CO 存在于最冷的分子氣體里,但現在很多研究發(fā)現,有大部分H2(可能超過一半)存在于溫度更亮的氣體(50~100 K)中。在這種溫度下,CO 將被光致解離為C 或C+,因此這些氣體是“CO 暗”的,即這部分H2的質量無法通過CO 的方法獲得[8,16,47,50–53]。在這樣的情況下,利用中紅外波段觀測的H2輻射(約100 K)通過一定的插值外推可得到分子氣體的總質量。這些進展有助于人們直接利用中紅外波段探測到的H2輻射來推測總分子氣體質量,特別是對那些利用傳統CO 來示蹤總分子氣體可能會具有很大偏差的星系,如低金屬豐度星系。

    3.1 模型和測量方法

    Togi 和Smith 于2016年通過中紅外探測的多條H2輻射建立了H2的溫度分布模型,并利用這一模型外推到低溫(約50 K)來獲取總H2質量[19]。他們認為H2的溫度分布可以用一條平滑的冪律譜來擬合(見圖8)。假設dN=mT ?ndT,用一條連續(xù)的冪律溫度分布擬合H2的旋轉譜線,這里dN是激發(fā)溫度在T和T+ dT之間的H2柱密度,n是冪律指數,與氣體的加熱機制有關(如氣體加熱來源可能為激波、光電效應和紫外光子等)。m是一個常數。在已知冪律指數n的情況下,對H2的溫度分布模型做積分。設定上限溫度Tu和外推的下限溫度Tl,即可得到總的分子氣體質量。根據Togi 和Smith 的研究結果,我們設定上限溫度為Tu= 2 000 K,而另外兩個參數Tl和n可以根據Levenberg-Marquardt 優(yōu)化訓練[52]來匹配,整個訓練樣本的平均外推溫度為=(49±9)K。

    圖8 NGC 5033 和3c 293 的激發(fā)圖[19]

    這個模型擬合的結果需要一個可靠的樣本(已知分子氣體質量)進行修正。修正以后,它將與任何一種間接的分子氣體示蹤者(如CO)無關。下限溫度設定為T?= 49 K,把擬合得到的樣本中所有星系分子氣體總質量與利用可信CO 示蹤得到的分子氣體總質量相比較,結果如圖9 所示??梢钥闯?,兩種方法得到的總分子氣體質量值比較接近,從而很好地驗證了Togi 和Smith 研究方法的有效性。圖9 中列出的天體均為利用CO 示蹤可以得到總分子氣體質量的星系。而如前面提到的,CO 示蹤方法在低金屬豐度或其他一些極端環(huán)境下將失效。Togi 和Smith 對于此類極端環(huán)境的星系也做了驗證,結果證明他們的方法也同樣有效。

    圖9 模型擬合得到的分子氣體總質量與CO 示蹤方法得到的分子氣體總質量對比[19]

    Togi 和Smith 認為該方法的適用性與星系類型無關。他們的樣本包括了恒星形成星系、活動星系核、低電離星系核、矮星系、射電星系、極亮紅外星系和亮紅外星系等多種電離機制可能完全不同的星系,他們發(fā)現該方法都可以適用。推導出的總分子氣體質量與CO 推導出的分子氣體質量都有很好的相關性。

    3.2 結果與分析

    Togi 和Smith 于2016年提出了利用中紅外波段的H2譜線強度推導總分子氣體質量的方法,該方法適用于任何類型的星系,也說明星系(無論類型)中總分子氣體質量與H2譜線強度相關。如前所述,中紅外H2輻射可能由激波激發(fā),與恒星形成的關系不大。人們通常認為恒星形成的強弱與分子氣體質量面密度緊密相關[6],而對這些星系,觀測得到的中紅外H2譜線強度卻與總分子氣體質量面密度關系不大。這是否與Togi 等人的結果相矛盾呢?

    為探討這一問題,我們首先研究Togi 和Smith 的模型,驗證作為模型輸入量的H2發(fā)射線輻射強度與作為模型輸出量的總分子氣體質量是否具有相關性。我們試圖研究樣本中H2(S(0)—S(3))譜線強度與分子氣體總質量M(H2)之間的關系,見圖10 a)。由于S(0)—S(3)都有測量值的星系數量較少,而同時有S(1)和S(2)的星系數量較多,我們也試圖用S(1)+S(2)的譜線強度和來代表H2的譜線輻射尋找H2輻射與分子氣體總質量之間的關系(見圖10 b))。整體上無論圖10 a)還是圖10 b),均沒有顯示出H2輻射與分子氣體總質量有明顯的關系。但是單獨分不同星系類型來看,如單獨看紫色方框(恒星形成星系),可以看到整體呈現出一定的正相關趨勢,對此,還需要進一步研究。同時,由ULIRG, LIRG 等得出的總分子氣體質量要比恒星形成星系和AGN 高。Togi 和Smith 的模型雖然利用H2輻射作為輸入量,但對一個包括多種不同類型星系的樣本而言,得出的總分子氣體質量已經與輸入量沒有明顯的相關性。那么模型中的哪些參量可導致模型的輸入量與輸出量之間弱相關?

    圖10 分子氣體的總質量與H2輻射強度和的關系圖

    H2溫度分布的冪律指數n是Togi 和Smith 的模型中主要的可變參量。在圖11 中,我們用顏色指示n的大小,再次展示圖10 a)的結果,發(fā)現冪律指數n并未在任何方向上表現出趨勢性,所以模型中輸入量與輸出量之間的弱相關不是由星系中H2溫度的分布導致的。

    圖11 分子氣體的總質量與H2輻射(S(0)—S(3))強度和的關系

    在圖10 b)中,我們發(fā)現ULIRG 中分子氣體的總質量明顯比恒星形成星系和活動星系核偏高,因此我們必須考慮模型中的另一個重要參數:星系的距離d。在研究輻射的總光度時,該參數以d2的形式引入到模型中。因此,在圖12 中,我們把總分子氣體質量除以d2,研究與觀測得到的H2輻射的關系。如果不考慮不同星系角大小的差別,該物理量基本等同于分子氣體的質量面密度,我們就稱其為分子氣體的質量面密度。我們發(fā)現圖12 顯示較強的相關性(整個樣本的相關性為0.83)。

    圖12 分子氣體的表面密度與H2輻射強度和的關系圖

    如果圖12 展示的相關性是真實的,而且Togi 和Smith 的方法是可行的,那么我們可以推斷中紅外波段觀測得到的H2輻射強度與總分子氣體的表面密度是相關的,而且這種相關性與星系類型無關。假定總分子氣體的表面密度是與恒星形成緊密相關,那么就會與AGN(包括高光度AGN 主導的ULIRG)中H2輻射主要由激波激發(fā)而與恒星形成關系不大相矛盾??梢詮囊韵聨讉€方面解釋這種矛盾:(1)傳統觀點認為分子氣體面密度與恒星形成緊密相關源于對致密分子云的探索。這些分子氣體溫度都比較冷,與恒星從致密冷的分子云中產生的傳統觀點相一致。近期人們發(fā)現越來越多的分子氣體可能大部分在溫度更高(約50 K)的區(qū)域內[8,16,47,50–53]??偟姆肿託怏w面密度與恒星形成的相關性如何,特別是對AGN 這樣的特殊星系,需重新進行研究,也許不能籠統地認為它們相關。(2)在第2 章的討論中,AGN 中的恒星形成強弱多用PAH 來表征。雖然通常人們認為PAH 是恒星形成強度的表征,但在AGN 中,PAH 會不會因為受到AGN 輻射的影響而不能準確地表征恒星形成的強弱?該問題并沒有被系統地探索和驗證過。(3)圖12 中的相關性在H2面密度方向還有一定的彌散性,最大可達一個數量級。第2 章中樣本都比較小,為更深入地探索這一問題,需要一個更大的樣本,該樣本應具有統一測量的多條H2輻射、PAH 輻射、硅酸鹽吸收強度等,并最好同時具有用其他手段測量的H2質量和面密度,以及其他表征恒星形成的物理量。這樣的樣本收集和分析工作目前還在進行中,我們將在今后的工作中進一步探討這一問題。

    4 總結與展望

    本文介紹了近20年來在中紅外波段對H2的研究進展。重點討論了兩個關鍵問題:(1)H2輻射與恒星形成的相關性;(2)針對H2輻射能否跟蹤總分子氣體質量。同時,我們還介紹了能夠通過中紅外H2譜線直接得到總分子氣體質量的最新方法。具體總結如下。

    (1)截至目前,H2的中紅外旋轉譜線已經在年輕恒星的外流、光解離區(qū)域、行星狀星云、超新星遺跡,以及星系中心的大區(qū)域等多種環(huán)境中被探測到。不同環(huán)境中探測到的H2譜線的激發(fā)機制也不同,因此關于中紅外波段發(fā)射的H2譜線是否可以用來追蹤恒星形成的觀點也不同。一般認為在恒星形成星系中,中紅外波段的H2譜線來源于恒星形成,因此可以成為恒星形成的示蹤者;而在非恒星形成主導的星系中,H2譜線激發(fā)大多來源于激波。

    (2)之前人們認為絕大部分H2存在于極低(10~20 K)的溫度下,無法進行有效的輻射。因此一般使用CO 分子的旋轉譜線來間接追蹤分子氣體的總質量。然而,在很多極端星際環(huán)境下,這種間接的方法將會失效。第3 章介紹從Togi 和Smith 提出的模型推得分子氣體總質量的方法,結果說明中紅外波段觀測得到的H2輻射可以追蹤總分子氣體質量,這一擬合結果獨立于其他的間接示蹤方法,并且不依賴于星系的類型。

    (3)雖然第3 章列出利用H2輻射強度推得總分子氣體質量的方法,但H2輻射與總分子氣體質量并不直接相關。在尋找H2的表面密度與H2的輻射強度的關系時,可以看到各類星系均表現出相同的關系(H2的表面密度與H2的輻射強度正相關)。

    未來的研究中,我們期望在以下幾個方面取得進展。

    (1)現有的數據都是以小樣本為基礎的,我們急需一個更大的能夠包含大部分普遍及特殊星系的,且具有一定的精度的樣本,去做大數據的中紅外H2與恒星形成之間的關系研究。在恒星形成示蹤者的選擇上可以采用多樣化的示蹤者(適用于不同類型的星系),以期望得到能夠普遍適用的方法。

    (2)Togi 和Smith 的模型還有待進一步的改正,可以利用這一模型得到大樣本(包括各種各樣的星系類型)星系中的分子氣體總質量與利用其他方法間接得到的分子氣體總質量對比,驗證模型的有效性。

    猜你喜歡
    質量
    聚焦質量守恒定律
    “質量”知識鞏固
    “質量”知識鞏固
    質量守恒定律考什么
    做夢導致睡眠質量差嗎
    焊接質量的控制
    關于質量的快速Q&A
    初中『質量』點擊
    質量投訴超六成
    汽車觀察(2016年3期)2016-02-28 13:16:26
    你睡得香嗎?
    民生周刊(2014年7期)2014-03-28 01:30:54
    国产乱来视频区| 亚洲精品aⅴ在线观看| 午夜视频国产福利| 国产精品久久久久久精品古装| 最近手机中文字幕大全| 精品一区二区三区四区五区乱码 | 人人澡人人妻人| 午夜av观看不卡| 亚洲成人av在线免费| av在线老鸭窝| 日韩 亚洲 欧美在线| 综合色丁香网| 视频区图区小说| 国产成人精品久久久久久| 久久国产精品男人的天堂亚洲 | 人体艺术视频欧美日本| 秋霞伦理黄片| 亚洲精品456在线播放app| 国产亚洲精品久久久com| 免费日韩欧美在线观看| 国产高清国产精品国产三级| 日本av免费视频播放| 人成视频在线观看免费观看| 亚洲经典国产精华液单| 亚洲精品久久成人aⅴ小说| 美女视频免费永久观看网站| 夜夜骑夜夜射夜夜干| 亚洲国产精品一区三区| 国产国拍精品亚洲av在线观看| 成人黄色视频免费在线看| 国产有黄有色有爽视频| 精品少妇黑人巨大在线播放| av在线老鸭窝| 在线亚洲精品国产二区图片欧美| 国产精品熟女久久久久浪| 国产成人av激情在线播放| 免费av中文字幕在线| 日产精品乱码卡一卡2卡三| 伊人亚洲综合成人网| 在线天堂最新版资源| 看免费成人av毛片| 高清视频免费观看一区二区| 日韩在线高清观看一区二区三区| 男人添女人高潮全过程视频| 又黄又爽又刺激的免费视频.| 三上悠亚av全集在线观看| 久久久久国产精品人妻一区二区| 老司机亚洲免费影院| 性色avwww在线观看| 亚洲av电影在线观看一区二区三区| 18禁裸乳无遮挡动漫免费视频| 亚洲美女黄色视频免费看| 久久人妻熟女aⅴ| 国产毛片在线视频| 在线观看一区二区三区激情| 高清毛片免费看| 欧美成人午夜精品| 99久国产av精品国产电影| 青青草视频在线视频观看| 精品午夜福利在线看| 国产黄色免费在线视频| 99国产综合亚洲精品| 亚洲成色77777| 亚洲性久久影院| 亚洲精品久久成人aⅴ小说| 久久久a久久爽久久v久久| 婷婷成人精品国产| 少妇人妻精品综合一区二区| 热re99久久国产66热| 国产av国产精品国产| 人人妻人人添人人爽欧美一区卜| 午夜久久久在线观看| 2018国产大陆天天弄谢| 日韩大片免费观看网站| 久久精品久久久久久噜噜老黄| 免费播放大片免费观看视频在线观看| 国产极品粉嫩免费观看在线| 免费在线观看黄色视频的| 国产成人免费无遮挡视频| 各种免费的搞黄视频| 亚洲欧洲国产日韩| 国产国语露脸激情在线看| 日本午夜av视频| 国产永久视频网站| 男女免费视频国产| 国产高清不卡午夜福利| 日韩熟女老妇一区二区性免费视频| 国产在线视频一区二区| 国产午夜精品一二区理论片| 秋霞伦理黄片| 久久精品熟女亚洲av麻豆精品| 久久精品国产a三级三级三级| 秋霞在线观看毛片| 午夜91福利影院| www.熟女人妻精品国产 | 卡戴珊不雅视频在线播放| 制服人妻中文乱码| 各种免费的搞黄视频| 超色免费av| 国产女主播在线喷水免费视频网站| 亚洲精品美女久久久久99蜜臀 | 久久这里只有精品19| 国产极品天堂在线| 午夜福利在线观看免费完整高清在| 尾随美女入室| 美女中出高潮动态图| 午夜影院在线不卡| 久久久久久久大尺度免费视频| 亚洲精品成人av观看孕妇| 亚洲国产精品一区三区| 亚洲精品色激情综合| 超色免费av| 色网站视频免费| 少妇熟女欧美另类| 久久久国产精品麻豆| 欧美日韩成人在线一区二区| 欧美xxⅹ黑人| 国产精品熟女久久久久浪| 亚洲伊人久久精品综合| 一区二区三区精品91| 如日韩欧美国产精品一区二区三区| 狠狠婷婷综合久久久久久88av| 99国产综合亚洲精品| 人体艺术视频欧美日本| 久久精品人人爽人人爽视色| 女人久久www免费人成看片| 亚洲精华国产精华液的使用体验| 久久狼人影院| 在线观看一区二区三区激情| 精品国产一区二区三区四区第35| 国产色婷婷99| 欧美97在线视频| 久久国产精品男人的天堂亚洲 | 最近最新中文字幕大全免费视频 | a级毛片黄视频| 欧美人与性动交α欧美软件 | 99re6热这里在线精品视频| 国产成人午夜福利电影在线观看| 18+在线观看网站| 90打野战视频偷拍视频| 精品熟女少妇av免费看| 欧美精品一区二区大全| 亚洲精品视频女| 制服诱惑二区| 少妇人妻 视频| 女人精品久久久久毛片| xxx大片免费视频| 欧美另类一区| 亚洲欧美日韩卡通动漫| av国产久精品久网站免费入址| 一级毛片我不卡| 亚洲国产最新在线播放| 五月天丁香电影| 亚洲精品中文字幕在线视频| 国产精品一二三区在线看| 中文字幕制服av| 在现免费观看毛片| 久久久久精品性色| 侵犯人妻中文字幕一二三四区| 下体分泌物呈黄色| 欧美精品亚洲一区二区| 又黄又粗又硬又大视频| 国产片特级美女逼逼视频| 男女高潮啪啪啪动态图| 又粗又硬又长又爽又黄的视频| 777米奇影视久久| 高清av免费在线| 久久影院123| 在线观看美女被高潮喷水网站| 亚洲色图综合在线观看| 国产精品熟女久久久久浪| 一级黄片播放器| 天堂俺去俺来也www色官网| 亚洲精品成人av观看孕妇| 免费黄色在线免费观看| 国产乱来视频区| 香蕉精品网在线| 乱码一卡2卡4卡精品| 午夜免费观看性视频| 国产毛片在线视频| 国产白丝娇喘喷水9色精品| 九九在线视频观看精品| 大话2 男鬼变身卡| 亚洲精品视频女| 久久久久精品久久久久真实原创| 日韩大片免费观看网站| 激情五月婷婷亚洲| 久久热在线av| 嫩草影院入口| 秋霞在线观看毛片| 国产爽快片一区二区三区| 亚洲国产毛片av蜜桃av| 国产又色又爽无遮挡免| 十八禁网站网址无遮挡| 永久免费av网站大全| 久久久久精品久久久久真实原创| 啦啦啦在线观看免费高清www| 一区二区三区四区激情视频| 一本久久精品| 人妻系列 视频| 高清毛片免费看| 日韩电影二区| 大话2 男鬼变身卡| 这个男人来自地球电影免费观看 | 国产精品99久久99久久久不卡 | 日本色播在线视频| 中国三级夫妇交换| 亚洲 欧美一区二区三区| 大香蕉97超碰在线| 美女视频免费永久观看网站| 日本欧美视频一区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 97精品久久久久久久久久精品| 国产精品蜜桃在线观看| 亚洲人成网站在线观看播放| 久久久久久久久久人人人人人人| 亚洲精品,欧美精品| 免费不卡的大黄色大毛片视频在线观看| 精品第一国产精品| 国产精品嫩草影院av在线观看| 欧美成人精品欧美一级黄| 成年动漫av网址| 少妇高潮的动态图| 日本黄大片高清| 99香蕉大伊视频| 男女高潮啪啪啪动态图| 亚洲高清免费不卡视频| 侵犯人妻中文字幕一二三四区| 男女午夜视频在线观看 | 中文欧美无线码| 久久精品国产自在天天线| 各种免费的搞黄视频| 日本av手机在线免费观看| 人妻少妇偷人精品九色| 欧美成人午夜免费资源| 久久婷婷青草| 波野结衣二区三区在线| 日产精品乱码卡一卡2卡三| av免费在线看不卡| 国精品久久久久久国模美| 51国产日韩欧美| 丝袜喷水一区| 国产精品人妻久久久久久| 一级毛片黄色毛片免费观看视频| 成人黄色视频免费在线看| 十分钟在线观看高清视频www| 丰满乱子伦码专区| 亚洲精品国产av成人精品| 久久久久久人人人人人| 亚洲第一av免费看| 亚洲在久久综合| 中文字幕人妻丝袜制服| 大香蕉久久成人网| 亚洲av福利一区| 91精品伊人久久大香线蕉| 国产精品国产三级国产专区5o| 久久久a久久爽久久v久久| 不卡视频在线观看欧美| 精品福利永久在线观看| 亚洲伊人久久精品综合| 最近中文字幕高清免费大全6| 水蜜桃什么品种好| av在线app专区| 久久久欧美国产精品| 国产熟女欧美一区二区| 亚洲天堂av无毛| 国产精品欧美亚洲77777| 精品国产国语对白av| 三级国产精品片| 两个人免费观看高清视频| 亚洲av综合色区一区| 久久久久久人妻| 久久毛片免费看一区二区三区| 大陆偷拍与自拍| 国产 一区精品| 国产激情久久老熟女| 一本大道久久a久久精品| 久久精品国产鲁丝片午夜精品| 在线精品无人区一区二区三| 午夜福利网站1000一区二区三区| 亚洲av成人精品一二三区| 国产亚洲最大av| 一本大道久久a久久精品| 妹子高潮喷水视频| 精品一区二区三卡| 成人免费观看视频高清| 亚洲美女黄色视频免费看| 国产亚洲午夜精品一区二区久久| 啦啦啦中文免费视频观看日本| 少妇 在线观看| 久久狼人影院| 精品熟女少妇av免费看| 草草在线视频免费看| 视频中文字幕在线观看| 国产乱来视频区| 一区二区日韩欧美中文字幕 | 亚洲精品久久久久久婷婷小说| 一边亲一边摸免费视频| 美女福利国产在线| 99热这里只有是精品在线观看| 男女无遮挡免费网站观看| 狠狠婷婷综合久久久久久88av| 少妇被粗大猛烈的视频| 美女xxoo啪啪120秒动态图| 亚洲欧美日韩另类电影网站| 少妇精品久久久久久久| 亚洲经典国产精华液单| av福利片在线| 校园人妻丝袜中文字幕| 啦啦啦啦在线视频资源| 国产午夜精品一二区理论片| 午夜福利网站1000一区二区三区| 日韩三级伦理在线观看| 亚洲婷婷狠狠爱综合网| 一级爰片在线观看| 欧美激情国产日韩精品一区| 天天操日日干夜夜撸| 又大又黄又爽视频免费| 日日摸夜夜添夜夜爱| 亚洲欧美一区二区三区国产| 亚洲欧美色中文字幕在线| 少妇高潮的动态图| 日本与韩国留学比较| 狂野欧美激情性xxxx在线观看| 欧美精品亚洲一区二区| 美女中出高潮动态图| 高清黄色对白视频在线免费看| 777米奇影视久久| 亚洲av免费高清在线观看| 久久久a久久爽久久v久久| 欧美日韩一区二区视频在线观看视频在线| www.色视频.com| 国产免费一区二区三区四区乱码| 老熟女久久久| 国产av精品麻豆| 日韩大片免费观看网站| 欧美 亚洲 国产 日韩一| 91成人精品电影| 国产精品久久久久久久电影| 亚洲成国产人片在线观看| 性色avwww在线观看| 国产精品麻豆人妻色哟哟久久| 久久人妻熟女aⅴ| 国产精品国产三级国产av玫瑰| 2021少妇久久久久久久久久久| 亚洲国产毛片av蜜桃av| 在线天堂最新版资源| 亚洲精品色激情综合| 午夜老司机福利剧场| 国产在线一区二区三区精| 波野结衣二区三区在线| 国产成人免费观看mmmm| 2021少妇久久久久久久久久久| 日韩精品免费视频一区二区三区 | 亚洲国产av影院在线观看| 午夜福利视频精品| 国产深夜福利视频在线观看| 欧美日韩综合久久久久久| 欧美bdsm另类| 国产免费一区二区三区四区乱码| 久久久久久久亚洲中文字幕| 纵有疾风起免费观看全集完整版| 啦啦啦在线观看免费高清www| a 毛片基地| 亚洲丝袜综合中文字幕| 男女免费视频国产| 免费av中文字幕在线| 建设人人有责人人尽责人人享有的| 久久精品aⅴ一区二区三区四区 | 在线观看人妻少妇| 性高湖久久久久久久久免费观看| 女人精品久久久久毛片| 久久久久网色| 国产又色又爽无遮挡免| 波多野结衣一区麻豆| 国产免费一区二区三区四区乱码| 久久国产精品男人的天堂亚洲 | 美女福利国产在线| 97在线人人人人妻| 国产福利在线免费观看视频| 久久ye,这里只有精品| 久久精品久久久久久久性| 岛国毛片在线播放| 亚洲国产av影院在线观看| 成人无遮挡网站| 一级爰片在线观看| 这个男人来自地球电影免费观看 | 免费少妇av软件| 亚洲精品日韩在线中文字幕| 97在线视频观看| 美女视频免费永久观看网站| 免费av不卡在线播放| av片东京热男人的天堂| 国内精品宾馆在线| 国产av码专区亚洲av| 国产麻豆69| 在线 av 中文字幕| 国产毛片在线视频| 搡女人真爽免费视频火全软件| 一区二区三区乱码不卡18| 秋霞伦理黄片| av播播在线观看一区| 国产免费一区二区三区四区乱码| 国产免费又黄又爽又色| a级毛片在线看网站| av在线播放精品| av黄色大香蕉| 永久网站在线| 伦理电影免费视频| 飞空精品影院首页| 中文字幕人妻丝袜制服| 午夜久久久在线观看| 一边摸一边做爽爽视频免费| 亚洲在久久综合| 亚洲av男天堂| 韩国av在线不卡| 亚洲,欧美,日韩| 色婷婷av一区二区三区视频| 咕卡用的链子| 久久久欧美国产精品| 欧美日韩精品成人综合77777| 咕卡用的链子| 国产午夜精品一二区理论片| 69精品国产乱码久久久| 亚洲婷婷狠狠爱综合网| 久久精品aⅴ一区二区三区四区 | 亚洲欧美成人精品一区二区| 99香蕉大伊视频| 在线 av 中文字幕| 国产毛片在线视频| 亚洲伊人久久精品综合| 水蜜桃什么品种好| 99re6热这里在线精品视频| 99热网站在线观看| 欧美精品亚洲一区二区| 岛国毛片在线播放| 亚洲精品色激情综合| 国产乱来视频区| 国产男女超爽视频在线观看| 精品久久国产蜜桃| 啦啦啦啦在线视频资源| 天堂中文最新版在线下载| 两个人看的免费小视频| 久久精品国产亚洲av天美| 精品人妻偷拍中文字幕| 久久精品久久精品一区二区三区| 高清毛片免费看| 久久久精品区二区三区| 国产精品.久久久| 777米奇影视久久| 另类精品久久| 亚洲第一区二区三区不卡| 久久人人爽av亚洲精品天堂| 又黄又爽又刺激的免费视频.| av在线观看视频网站免费| 99九九在线精品视频| 久久毛片免费看一区二区三区| 成人国产麻豆网| 侵犯人妻中文字幕一二三四区| 大香蕉久久网| 美女国产视频在线观看| 免费黄色在线免费观看| 国产av精品麻豆| 人人澡人人妻人| 99久久人妻综合| 亚洲av中文av极速乱| 日韩欧美精品免费久久| 国产伦理片在线播放av一区| 成人国语在线视频| 99热这里只有是精品在线观看| 99国产综合亚洲精品| 色视频在线一区二区三区| 亚洲国产看品久久| 中文字幕人妻熟女乱码| 水蜜桃什么品种好| 90打野战视频偷拍视频| 91久久精品国产一区二区三区| 免费黄频网站在线观看国产| 国产精品成人在线| a级毛片黄视频| 波野结衣二区三区在线| 欧美少妇被猛烈插入视频| av在线app专区| 国产精品一国产av| 亚洲欧美中文字幕日韩二区| 夫妻性生交免费视频一级片| 天堂8中文在线网| av在线老鸭窝| 国产一区二区在线观看日韩| 97精品久久久久久久久久精品| 插逼视频在线观看| 亚洲色图综合在线观看| 国产av一区二区精品久久| 国产成人91sexporn| 超碰97精品在线观看| 中文天堂在线官网| 亚洲精品av麻豆狂野| 国产精品人妻久久久久久| 中国三级夫妇交换| 18禁观看日本| 国产黄频视频在线观看| 五月天丁香电影| 国产精品久久久久久久电影| 观看美女的网站| 国产精品国产av在线观看| 在线观看三级黄色| 乱人伦中国视频| 日本午夜av视频| 观看av在线不卡| 天堂俺去俺来也www色官网| 日产精品乱码卡一卡2卡三| 精品久久久精品久久久| 亚洲久久久国产精品| 欧美另类一区| 精品一区二区三卡| 亚洲伊人久久精品综合| 免费看av在线观看网站| 一级片免费观看大全| 欧美日韩国产mv在线观看视频| kizo精华| 日韩av免费高清视频| 极品少妇高潮喷水抽搐| av福利片在线| 天堂俺去俺来也www色官网| 成年美女黄网站色视频大全免费| 亚洲成人手机| 日韩,欧美,国产一区二区三区| 精品一区在线观看国产| 丁香六月天网| 精品国产乱码久久久久久小说| 人妻系列 视频| 在线看a的网站| 夜夜骑夜夜射夜夜干| www.av在线官网国产| 久久午夜综合久久蜜桃| 成年女人在线观看亚洲视频| 亚洲国产欧美在线一区| 91成人精品电影| 桃花免费在线播放| 伊人久久国产一区二区| 最新中文字幕久久久久| 黄网站色视频无遮挡免费观看| 永久免费av网站大全| 日韩免费高清中文字幕av| 精品酒店卫生间| 日韩成人av中文字幕在线观看| 国产午夜精品一二区理论片| 又粗又硬又长又爽又黄的视频| 国产深夜福利视频在线观看| 99香蕉大伊视频| 啦啦啦啦在线视频资源| 日产精品乱码卡一卡2卡三| 狂野欧美激情性bbbbbb| 亚洲精品aⅴ在线观看| 国产精品嫩草影院av在线观看| 色网站视频免费| 美女内射精品一级片tv| 大香蕉久久网| 色视频在线一区二区三区| 熟女av电影| 在线观看国产h片| 日韩欧美一区视频在线观看| 777米奇影视久久| 国产色婷婷99| 中文精品一卡2卡3卡4更新| 国产精品麻豆人妻色哟哟久久| 99热国产这里只有精品6| 国产成人欧美| 一区在线观看完整版| 啦啦啦啦在线视频资源| 一区二区三区精品91| 一边摸一边做爽爽视频免费| av在线观看视频网站免费| 久久人人爽人人片av| 欧美亚洲 丝袜 人妻 在线| a 毛片基地| 日韩人妻精品一区2区三区| 人妻系列 视频| 免费高清在线观看视频在线观看| 少妇熟女欧美另类| 亚洲精品自拍成人| 亚洲精品成人av观看孕妇| 欧美日韩成人在线一区二区| 99热这里只有是精品在线观看| 大香蕉久久成人网| 在线观看www视频免费| 91aial.com中文字幕在线观看| 天堂俺去俺来也www色官网| 久久99精品国语久久久| 丰满饥渴人妻一区二区三| av播播在线观看一区| 超碰97精品在线观看| 哪个播放器可以免费观看大片| 亚洲国产欧美日韩在线播放| 最新中文字幕久久久久| 亚洲国产成人一精品久久久| 免费播放大片免费观看视频在线观看| 午夜精品国产一区二区电影| 99香蕉大伊视频| 在线亚洲精品国产二区图片欧美| 草草在线视频免费看| 精品亚洲成a人片在线观看| 各种免费的搞黄视频| 国产色婷婷99| 欧美另类一区| 免费在线观看完整版高清| 国产成人av激情在线播放| 十八禁网站网址无遮挡| 在线观看免费高清a一片| 久久女婷五月综合色啪小说| av电影中文网址| 国产福利在线免费观看视频| 下体分泌物呈黄色| 婷婷色麻豆天堂久久| 亚洲精品中文字幕在线视频|