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

    巡天項(xiàng)目中的白矮星搜尋

    2019-10-24 06:03:52羅阿理
    天文學(xué)進(jìn)展 2019年3期
    關(guān)鍵詞:白矮星巡天恒星

    孔 嘯,羅阿理

    (1.中國科學(xué)院 國家天文臺(tái) 光學(xué)天文重點(diǎn)實(shí)驗(yàn)室,北京 100101;2.中國科學(xué)院大學(xué),北京 100049)

    1 引 言

    當(dāng)恒星內(nèi)部的核反應(yīng)因燃料不足而停止后,主序星一般會(huì)坍縮成白矮星。這些奇特的天體比黯淡的紅色矮星更暗,但它們的顏色更藍(lán),所以才被稱為白矮星。白矮星的初始質(zhì)量為8.5M⊙~10.6M⊙[1],現(xiàn)階段觀測到的樣本質(zhì)量大多為0.5M⊙~0.8M⊙,而半徑卻與地球半徑的數(shù)量級(jí)相同。理論上來講,白矮星的最小質(zhì)量約為0.30M⊙~0.45M⊙[2]。此外,如果一顆恒星能夠演化成質(zhì)量更小的白矮星,那么,它在主序階段演化所需的時(shí)間就會(huì)比當(dāng)前宇宙的年齡還要長。最初,理論學(xué)家稱白矮星為簡并星(degenerate star)或簡并矮星(degenerate dwarf)?,F(xiàn)在,白矮星才是最常用的稱呼。

    根據(jù)觀測數(shù)據(jù)和理論計(jì)算,人們相信在銀河系中超過97% 的恒星最終會(huì)演化成白矮星[3]。由于沒有核聚變,白矮星僅僅依靠殘留的引力收縮輻射來提供能量,所以,白矮星的冷卻時(shí)間一般都非常長。

    人類觀測到的白矮星保留了關(guān)于銀河系形成早期的大量信息,如光度函數(shù)或年老成員的溫度等。因?yàn)榘装堑那吧硇窃跐u近巨星支(AGB)頂部失去了C,N,O,所以它們是銀河系化學(xué)演化的重要貢獻(xiàn)者,甚至有可能是維持生命的化學(xué)物質(zhì)的重要來源。盡管人們對(duì)白矮星有著非常濃厚的興趣,但由于白矮星通常都極其黯淡,以至于上百年來,天文學(xué)家們對(duì)白矮星的觀測和研究都進(jìn)展緩慢。直到20 世紀(jì)末21 世紀(jì)初,隨著斯隆數(shù)字巡天(Sloan Digital Sky Survey,SDSS)和郭守敬望遠(yuǎn)鏡(Large Sky Area Multi-Object Fiber Spectroscopic Telescope,LAMOST)的光譜巡天等大天區(qū)尺度項(xiàng)目的開展,我們才得到較多白矮星樣本,可用來研究分析。

    本文簡述了白矮星的光譜分類、顏色判據(jù)、巡天觀測和物理參數(shù)等信息。

    2 光譜型和顏色特征

    2.1 白矮星的光譜型分類

    如果一顆白矮星的表面大氣主要由H 構(gòu)成,那么它就是DA 型白矮星。DA 型白矮星的數(shù)量占全部白矮星的約80%,另外約20% 主要為大氣成分由He 主導(dǎo)的白矮星,它們包括DB 型和DO 型白矮星。全部白矮星的分類如下。

    (1)DA 型,其表面大氣主要由H 構(gòu)成,光譜表現(xiàn)為巴爾末線系(Balmer lines);

    (2)DB 型,其表面大氣主要由中性He (即HeⅠ)組成;

    (3)DO 型,其表面大氣主要由一次電離的He (即HeⅡ)組成;

    (4)DQ 型,其光譜表現(xiàn)為C2斯旺譜帶(C2Swan band)或中性C (即CⅠ)線;

    (5)DZ 型,其光譜中金屬線(如CaⅡH 和CaⅡK)占主導(dǎo)位置;

    (6)DC 型,其光譜中無明顯譜線,表現(xiàn)為較平滑的連續(xù)譜。

    在這些類型的白矮星中,DO 型白矮星的溫度相對(duì)較高,有效溫度大約為45 000 K;DB 型白矮星的溫度較低,大部分在30 000 K 以下;當(dāng)有效溫度降到10 000 K 以下時(shí),He在其光譜中幾乎不可見,如DC,DQ,DZ 型白矮星[4]。

    由于H 包層的DA 型白矮星數(shù)量眾多,因此,對(duì)DA 型白矮星的研究成果也最多。直到近期的大型巡天望遠(yuǎn)鏡上線后,人們針對(duì)非DA 型白矮星(如DB 型)的專門研究才豐富起來。由于DB 型白矮星的大氣層由近乎純凈的中性He 構(gòu)成,因此,它可以作為宇宙中H 缺失恒星的最典型樣本。許多H 主導(dǎo)的DA 型白矮星會(huì)逐漸轉(zhuǎn)變成DB 型白矮星。隨著白矮星的逐漸冷卻,DA 型和非DA 型白矮星的數(shù)量會(huì)發(fā)生變化,其比例為有效溫度的函數(shù)[3]。如今,DB 型白矮星樣本的擴(kuò)充,給人們提供了研究白矮星演化的可能。

    Cukanovaite 等人[5]為He 主導(dǎo)的DB 型白矮星建立了第一個(gè)3D 模擬網(wǎng)格。這個(gè)模擬是使用CO5BOLD 輻射流體力學(xué)軟件[6]來進(jìn)行計(jì)算的,覆蓋的有效溫度范圍為12 000~34 000 K,表面重力范圍為7.5~9.0 dex。他們對(duì)比了從三維模擬計(jì)算得到的合成光譜與差分方法下一維模型光譜之間的區(qū)別,發(fā)現(xiàn)這種區(qū)別源于不同的對(duì)流模型。

    關(guān)于DB 型白矮星,最讓人熟知的就是所謂的DB 缺失(DB-gap):在有效溫度Teff為30 000~45 000 K 時(shí),幾乎只有DA 型白矮星。同時(shí),He 白矮星的溫度或高于這個(gè)溫度區(qū)間,如DO 型白矮星;或低于這個(gè)溫度區(qū)間,如DB 型白矮星。其實(shí),DB 缺失并不是絕對(duì)的,近些年的巡天項(xiàng)目都曾捕捉到有效溫度在DB 缺失范圍內(nèi)的DB 型白矮星,只是數(shù)量稍微少一些。因此,有人認(rèn)為,把DB 缺失改為DB 缺少(DB-deficiency)更恰當(dāng)一些[7]。

    在已知的白矮星光譜樣本中,約有3% 的光譜存在金屬線(DZ 型),這很可能是由于星體周邊星際塵埃的增加引起的[8–11]??傮w上來說,DZ 型白矮星光譜中Ca 和Mg 譜線是最強(qiáng)的。

    以上的分類以及對(duì)恒星的參數(shù)測量都是基于恒星表面大氣的物理性狀。通過研究恒星自身存在的震蕩現(xiàn)象,我們可以研究它們的內(nèi)部結(jié)構(gòu),完善恒星結(jié)構(gòu)和演化模型。如果一顆白矮星存在星震,人們通常會(huì)在它的分類后加一個(gè)V 的標(biāo)識(shí)[12],比如Landolt[13],Winget等人[14]和McGraw 等人[15]分別第一次發(fā)現(xiàn)了DAV,DBV 和DOV 型變星。

    由于重力擴(kuò)散會(huì)導(dǎo)致最輕的元素逐漸上浮到星體表面,因此,正常的白矮星大氣都是被輕元素覆蓋著。以上對(duì)白矮星的主流光譜學(xué)分類正是根據(jù)表面大氣組成來進(jìn)行的。但也有個(gè)別白矮星,其表面大氣組成以重元素為主,如SDSS J124043.01+671034.68[16]的大氣層幾乎完全由O 組成,其數(shù)量為其他元素的25 倍以上。它的大氣中還有少量的Ne 和Mg,但人們沒有觀測到H 或He 的譜線。根據(jù)目前有關(guān)白矮星的理論,這種情況比較特殊。O,Ne和Mg 都是大質(zhì)量白矮星的前身星形成階段后期C 燃燒的產(chǎn)物。這類白矮星也被稱為Dox,它對(duì)研究白矮星的形成有非常大的幫助。

    2.2 類型之間的轉(zhuǎn)換

    Fontaine 和Wesemael[17]曾經(jīng)提出一種模型,在DB 缺失的高溫端是H 缺失的熱白矮前身星(PG 1159 或者DO 型白矮星)。在恒星大氣層中只有很少量的H 元素。當(dāng)溫度下降后,H 會(huì)逐漸對(duì)流到恒星大氣層表層,形成一種富H 的大氣層。正因?yàn)楫?dāng)白矮星的溫度達(dá)到Teff≈45 000 K 時(shí),該白矮星會(huì)轉(zhuǎn)變成DA 型,所以在熱的前身星中H 的總量有上限MH≈10?16M⊙。如果白矮星大氣中H 層較薄,那么,該H 層會(huì)在He 氣包層頂部漂浮,以達(dá)到擴(kuò)散平衡。此時(shí),星體將呈現(xiàn)為DAO 型白矮星,它具有化學(xué)分層大氣的特征[18]。在DB 缺失的下邊緣(Teff30 000 K),由于少量的H 與厚重的He 包層的對(duì)流,大量白矮星表現(xiàn)為He 白矮星(DB 型白矮星)。

    DB 型白矮星的有效溫度不在30 000~45 000 K 范圍內(nèi),說明白矮星由富H (DA 型)到富He (DB 型)的轉(zhuǎn)換只在有效溫度30 000 K 以下才會(huì)發(fā)生。導(dǎo)致這種轉(zhuǎn)換過程發(fā)生的唯一可能的物理機(jī)制是:在此溫度下,DA 型白矮星的大氣層中較稀薄的H 層與深處濃厚的He 包層發(fā)生輻射對(duì)流。然而,Bergeron 等人[19]曾指出,DB 型白矮星數(shù)量顯著增加的情況只會(huì)發(fā)生在Teff≈20 000 K 以下,而不是所謂的30 000 K。

    嚴(yán)格來說,一半以上已知DB 型白矮星都是具有較弱H 線的DB 型白矮星,即DBA型白矮星[7,20]。如果只考慮有效溫度在20 000 K 以下的樣本,這個(gè)比例還會(huì)更高。其實(shí)在這個(gè)溫度下,在光學(xué)波段仍然會(huì)有很多純He (H/He 10?6)白矮星。

    對(duì)于有效溫度在30 000~45 000 K 的DB 型白矮星,顯然,在大氣包層中并沒有足夠的H 能使它們成為DA 型白矮星,并且它們很可能在其整個(gè)生命周期中演變?yōu)楦籋e 白矮星,Teff25 000 K 的樣本中最熱的DB 型白矮星可能是這樣的情況。另外,如果DA 型白矮星大氣中的H 較薄,含量較少(MH≈10?15M⊙),在溫度達(dá)到Teff≈24 000 K 時(shí),DA型白矮星會(huì)由于大氣層對(duì)流而變成DBA 型白矮星。當(dāng)恒星中存在的少量H 在較低的有效溫度下被He 對(duì)流區(qū)進(jìn)一步稀釋,DBA 型白矮星可能會(huì)很快地轉(zhuǎn)變成純凈的DB 型白矮星。正如Bergeron 等人[19]所指出的,因?yàn)槟壳耙阎顭岬腄Q 型白矮星的溫度接近18 000 K,有一部分必定在較低的有效溫度下以某種方式轉(zhuǎn)變成DB 型白矮星。這一過程的機(jī)制目前尚不清楚。但是,由于已知的熱DQ 型白矮星數(shù)量很少,因此,從整體的角度來看,這個(gè)過程可忽略不計(jì)。

    在SDSS 的白矮星光譜樣本中,較冷的He 白矮星(Teff12 000 K)存在H 的痕跡是很普遍的現(xiàn)象。有些白矮星(如DZA 型白矮星)的光譜中帶有金屬線,這可能是由于彗星、宇宙塵?;蛐⌒行堑任镔|(zhì)的污染引起的。這些冷星的數(shù)量特點(diǎn)也可能是由于選擇效應(yīng)導(dǎo)致的:溫度較低時(shí)探測Hα 譜線會(huì)比較困難。

    當(dāng)有效溫度降到約12 000 K 以下時(shí),白矮星中H 的豐度會(huì)下降到可探測范圍的下限,所以這種情況下DBA 型白矮星不會(huì)轉(zhuǎn)變成冷的富He 的DA 型白矮星,而是會(huì)直接變成沒有譜線信息的DC 型白矮星[7]。相反,薄H 層與更深的He 對(duì)流區(qū)的對(duì)流混合模型最有可能解釋冷的富He 的DA 或DZA 型白矮星中H 的存在。混合后,這些恒星將以幾乎恒定的H豐度演化;當(dāng)Hα 的譜線強(qiáng)度低于探測閾值時(shí),它們最終變成DC 或DZ 型白矮星。在溫度低于約10 000 K 的白矮星中,如果DA 型白矮星大約占一半的比例,那么大約有20% 的白矮星是DB 類型[19],約40% 的冷DA 型白矮星最終會(huì)發(fā)生對(duì)流混合。

    在大型巡天項(xiàng)目誕生之前,已知的白矮星寥寥無幾,因此,對(duì)各類白矮星的研究很大程度上建立在理論模型的基礎(chǔ)之上。究竟DA 型和DB 型白矮星的占比各多少?DBA 型和DAB 型等過渡形態(tài)白矮星的物理性質(zhì)有無特別之處?各種類型的白矮星在銀河系內(nèi)的空間分布和運(yùn)動(dòng)狀態(tài)怎樣?······或許隨著巡天項(xiàng)目的不斷開展,這些問題的答案將會(huì)逐漸清晰。

    2.3 白矮星的顏色特征

    在巡天望遠(yuǎn)鏡拍攝的光譜中,溫度約8 000 K 以上的白矮星與大量星系和主序星的顏色是有區(qū)別的。少量例外的情況是,擁有晚型星伴星的白矮星雙星系統(tǒng),它們的伴星仍然會(huì)在顏色上保留主序星的特點(diǎn)。

    圖1展示的是SDSS 巡天項(xiàng)目數(shù)據(jù)中富H 白矮星(DA 型)和富He 白矮星(DB/DO 型)的u′?g′與g′?r′顏色分布,分別用紅色三角形和藍(lán)色十字標(biāo)識(shí);小黑點(diǎn)代表普通主序星分布在右上方。白矮星的顏色分布與主序星有一定的區(qū)別,可用一些判據(jù)對(duì)其進(jìn)行初步篩選[21]。

    對(duì)于較冷的白矮星(Teff<8 000 K),巴爾末跳變和H 線變得很弱,以至于He 白矮星難以表現(xiàn)出激發(fā)態(tài)的HeⅠ線。這時(shí)白矮星與主序星的顏色可能會(huì)混雜在一起,使挑選特征較弱的白矮星變得更加困難。然而,對(duì)于很多較冷的非DA 型白矮星,如果信噪比較高,也會(huì)顯示出一些較弱的特征,它們可能是C 分子帶(通常是C2)或原子譜線(DQ 型),也可能是重元素的譜線(Ca II H 和Ca II K 等)。之所以會(huì)有C 元素,是由于當(dāng)He 包層下的對(duì)流區(qū)域變得足夠深時(shí),星體內(nèi)一部分富C 的核被翻轉(zhuǎn)到表面。而當(dāng)白矮星偶然間穿過一片云或高密度氣體的區(qū)域時(shí),星體表面的重元素含量可能會(huì)升高;經(jīng)過一段時(shí)間,這些重元素也可能因重力對(duì)流而不可見。

    當(dāng)DQ 或DZ 型白矮星的特征變得明顯,它們的顏色也可能會(huì)偏離正常范圍。在圖1中,g′的C2吸收帶可能會(huì)把較明顯的DQ 型白矮星移到左上角;一條在u′波段的流量較小的典型DZ 光譜,其顏色可能會(huì)朝著其他方向移動(dòng)。類似這樣的白矮星非常稀少,其原因是,它們溫度太低,以至于無法展示出大氣中He 的成分。但是,它們是冷白矮星樣本不可或缺的組成成分。

    然而,當(dāng)有效溫度降低至5 000 K 以下后,幾乎所有的明顯特征似乎都消失了,這樣的星體幾乎都變成了DC 型白矮星。它們的顏色與星族Ⅱ主序亞矮星非常相似,但是它們?cè)谘蔡祉?xiàng)目的海量數(shù)據(jù)中更加難以搜尋到。

    3 觀測及獲取

    3.1 早期觀測

    在20 世紀(jì)80年代以前,白矮星的觀測數(shù)據(jù)主要通過以下幾種途徑來獲?。鹤孕行潜碇械陌敌牵绕涫撬{(lán)顏色的目標(biāo)星體;總體上在高銀緯的區(qū)域?qū)λ{(lán)色暗星的巡天觀測數(shù)據(jù);有自行或者視差的亮星,其伴星更暗。

    由于早期沒有大口徑天文望遠(yuǎn)鏡和高精度觀測儀器,人們所知的白矮星屈指可數(shù)。Luyten[22,23]在20 世紀(jì)70年代通過早期的自行巡天(如Luyten Palomar)觀測,發(fā)布了疑似白矮星的星表。隨后,Green 等人[24]對(duì)偏藍(lán)色天體的巡天觀測結(jié)果,使得被證認(rèn)的簡并星數(shù)量翻倍。McCook 和Sion[25]在1977年發(fā)布了一個(gè)光譜證認(rèn)的白矮星星表。

    限于當(dāng)時(shí)的設(shè)備和技術(shù),已知的白矮星樣本不可避免地偏向于熱星、運(yùn)動(dòng)速度較快的星、有較大半徑和較小質(zhì)量的星[26]等,而且其數(shù)量屈指可數(shù)。

    3.2 近期觀測

    近十幾年,涌現(xiàn)出很多個(gè)高精度大尺度的巡天項(xiàng)目,每個(gè)巡天項(xiàng)目都會(huì)持續(xù)拍攝很多天體,其中就包含數(shù)量巨大的白矮星。于是,隨著已知的白矮星數(shù)量近乎成指數(shù)級(jí)增長,相關(guān)的研究工作也不斷地深入。SDSS,GALEX (Galaxy Evolution Explorer),LAMOST,GAIA(Global Astrometric Interferometer for Astrophysics)這幾個(gè)巡天項(xiàng)目在白矮星的觀測和研究上作出了較大的貢獻(xiàn)。

    3.2.1 SDSS

    正在進(jìn)行的SDSS 始于2000年,到目前為止,它提供的白矮星光譜數(shù)量最多。該項(xiàng)目旨在通過獲取海量測光和光譜數(shù)據(jù),來研究宇宙的大尺度結(jié)構(gòu)、星系的形成和演化等天體物理學(xué)領(lǐng)域的重大前沿課題。SDSS 所使用的望遠(yuǎn)鏡口徑為2.5 m,位于美國新墨西哥州阿帕奇天文臺(tái)(Apache Point Observatory)。對(duì)SDSS 的每個(gè)觀測天區(qū)約7 deg2的面積,人們采用一個(gè)打了640 個(gè)孔(對(duì)應(yīng)觀測目標(biāo)及定標(biāo)星)的焦面金屬板,因此,一次曝光可以拍攝640 條光譜,光譜的波長覆蓋范圍大約是3 800~9 200 ?A,分辨率R ≈1 900。SDSS 光譜巡天的主要目標(biāo)是獲取超過106個(gè)亮星系[27]、105個(gè)亮類星體和數(shù)量可觀的恒星光譜。另外一小部分光纖還會(huì)被用于其他研究項(xiàng)目,例如激變變星或白矮星等。

    從2000年至今,SDSS 先后完成了4 期巡天[28]:SDSSⅠ,Ⅱ[29,30],Ⅲ[31,32],Ⅳ[33]。每期巡天都包含了數(shù)個(gè)巡天計(jì)劃。

    Kleinman 等人[34]在SDSS DR1 (data release 1)中搜尋到了2 551 條白矮星光譜,這些數(shù)據(jù)覆蓋了1 360 deg2的范圍;從DR4[35]的發(fā)布數(shù)據(jù)中,Eisenstein 等人[36]將光譜證認(rèn)的白矮星數(shù)量增加到9 316 顆。借助于Koester[37]計(jì)算得到的大氣模型,他們對(duì)全部的DA 和DB 型白矮星進(jìn)行全譜擬合,得到了大氣參數(shù),其中,重力加速度(lgg)可達(dá)9.0。

    DR7[38]的光譜中包含14 120 條DA 型白矮星和1 011 條DB 型白矮星光譜[21];DR10[32]新增了6 687 個(gè)DA 型白矮星、450 個(gè)DB 型白矮星和若干其他類型的白矮星光譜數(shù)據(jù)[39];到了DR12[40],光譜數(shù)據(jù)總量超過4.35×106條,其中包含了6 576 顆新證認(rèn)的白矮星[41]。Koester 和Kepler[42]選取了其中信噪比大于10 的He 主導(dǎo)的1 107 顆DB 型白矮星,測量并分析了它們的大氣參數(shù)、質(zhì)量和空間分布等信息。

    截至目前,通過巡天望遠(yuǎn)鏡16年的觀測數(shù)據(jù)和科研人員對(duì)數(shù)據(jù)的分析處理,SDSS 發(fā)布了4 851 200 條光譜數(shù)據(jù)(SDSS DR14[43])。在最新發(fā)布的DR14 數(shù)據(jù)集中,DB 型白矮星的數(shù)量增加到了1 999 顆[20]。

    圖2展示的是在銀道坐標(biāo)系下SDSS 的所有光譜數(shù)據(jù),可以看出,白矮星的數(shù)量不在少數(shù)。

    圖2 銀道坐標(biāo)系下SDSS 光譜數(shù)據(jù)的全天腳印圖(綠色點(diǎn))

    3.2.2 GALEX

    GALEX[44–46]項(xiàng)目的主要目的是利用軌道空間望遠(yuǎn)鏡來觀測紫外波段的天體。它于2003年4月28日美國東部時(shí)間上午8:00 被火箭發(fā)射到軌道上。

    通過GALEX 得到的數(shù)據(jù)可用來分析星系是如何發(fā)展和變化的,亦可用于研究恒星形成的原因。在加利福尼亞理工學(xué)院研究團(tuán)隊(duì)的帶領(lǐng)下,人們利用GALEX 進(jìn)行了數(shù)項(xiàng)天空測量,包括銀河系以外的紫外全天測量。在任務(wù)期間,人們將制作出星系宇宙的第一張綜合地圖,以便進(jìn)一步了解像我們銀河系這樣的星系是如何形成的。由于白矮星的特征主要集中在藍(lán)端,所以GALEX 的紫外數(shù)據(jù)對(duì)于研究白矮星有著舉足輕重的作用。

    到目前為止,GALEX 已發(fā)布了7 次數(shù)據(jù)[47],分別是GR1 (GALEX Release 1)到GR7[48]。2013年2月27日,GALEX 發(fā)布了最新的GR6/GR7 數(shù)據(jù),包括34 285 個(gè)AIS(All Sky Imaging Survey)天區(qū)、720 個(gè)DIS (Deep Imaging Survey)天區(qū)、6 964 個(gè)MIS(Medium Imaging Survey)天區(qū)、716 個(gè)NGS (Nearby Galaxy Survey)天區(qū)、2 112 個(gè)GII(Guest Investigator)天區(qū)、87 個(gè)CAI (Calibration Imaging Survey)天區(qū)和311 個(gè)光譜(SPECTRA)天區(qū)。

    在SDSS 探測到的40 000 個(gè)白矮星中,約32 000 個(gè)有GALEX 紫外線觀測值,6 000個(gè)有SST (Spitzer Space Telescope)或WISE (Wide-field Infrared Survey Explorer)紅外觀測值[49]。

    3.2.3 LAMOST

    大天區(qū)面積多目標(biāo)光纖光譜天文望遠(yuǎn)鏡[50](LAMOST)又稱為郭守敬望遠(yuǎn)鏡,是一架中星儀式反射施密特望遠(yuǎn)鏡,坐落于河北省承德市興隆縣的興隆觀測基地[51]。由于同時(shí)擁有大視場(5?)和大口徑(有效口徑達(dá)3.6~4.9 m),LAMOST 可以在1 次的曝光中拍攝4 000個(gè)天體,并通過焦面上的4 000 根光纖在CCD 上成像。作為后起之秀,LAMOST 望遠(yuǎn)鏡是目前世界上光譜獲取率最高的天文望遠(yuǎn)鏡。

    LAMOST 擁有16 臺(tái)光譜儀,每臺(tái)都連接250 根光纖。當(dāng)觀測夜降臨后,每臺(tái)光譜儀通過2 個(gè)CCD 相機(jī),分別觀測紅端(3 700~5 900 ?A)和藍(lán)端(5 700~9 000 ?A)的數(shù)據(jù)流量。直接拍攝到的是星光經(jīng)過色散后的圖像,通過LAMOST 2D Pipeline 軟件[52]的抽譜[53]、定標(biāo)、減天光[54]、紅藍(lán)端數(shù)據(jù)合并等一系列的操作后,才能生成一維光譜。經(jīng)過LAMOST 1D Pipeline 軟件對(duì)光譜數(shù)據(jù)的分析和制作后,光譜和星表等就會(huì)被發(fā)布在網(wǎng)站上[55]。

    從2011年到如今,LAMOST 總共發(fā)布了6 個(gè)完整的數(shù)據(jù),最新的DR6 包含14 種光譜型的9 919 106 條光譜數(shù)據(jù)。DR1 到DR6 的發(fā)布數(shù)據(jù)情況見表1和圖3。

    表1 LAMOST 數(shù)據(jù)發(fā)布情況

    圖3 銀道坐標(biāo)系下LAMOST DR5 發(fā)布數(shù)據(jù)的腳印圖(綠色點(diǎn))

    圖3展示的是LAMOST 的巡天足跡。與圖2類似,藍(lán)色的點(diǎn)為SDSS DR14 和LAMOST DR5 都發(fā)布過的同源光譜數(shù)據(jù)。雖然LAMOST DR5 數(shù)據(jù)總量是SDSS DR14 的2 倍左右,但是LAMOST DR5 中的白矮星(紅點(diǎn))數(shù)量卻少很多,其主要原因是兩個(gè)巡天項(xiàng)目的選源策略和極限星等不同。

    LAMOST 在1 次曝光的過程中可以捕捉4 000 條光譜,但它捕獲暗星的能力卻不如SDSS。圖4展示的是SDSS DR14 和LAMOST DR5 的發(fā)布數(shù)據(jù)中白矮星g 波段的星等分布,其中,兩個(gè)分布直方圖均作了歸一化處理,SDSS DR14 (藍(lán)色)中白矮星的真實(shí)數(shù)量需要乘以因子4 217,LAMOST DR5 (紅色)的因子是536。由于白矮星的表面亮度較低,因此,SDSS 可以比LAMOST 獲取更多白矮星光譜。

    圖4 SDSS DR14 和LAMOST DR5 數(shù)據(jù)中白矮星的g 波段星等分布

    巡天項(xiàng)目的選源策略決定其獲取何種類型的光譜。采用何種策略,不僅要根據(jù)儀器自身的能力,也要根據(jù)巡天目標(biāo)的性質(zhì)。SDSS 在選源上的白矮星數(shù)量就比LAMOST 多很多(如表2所示),由兩個(gè)巡天項(xiàng)目發(fā)布的星表中選源類型和目標(biāo)類型的數(shù)量可以看出,SDSS DR12 的光譜數(shù)據(jù)中,白矮星或矮星的選源數(shù)量有近3.9×104個(gè),而LAMOST DR5 中只有幾百個(gè)。

    表2 SDSS DR12 和LAMOST DR5 的選源

    雖然白矮星的比例較低,但在LAMOST 海量的光譜數(shù)據(jù)中,仍然會(huì)有數(shù)量可觀的白矮星光譜供搜尋和研究。LAMOST 的先導(dǎo)巡天(Pilot)發(fā)布數(shù)據(jù)包含有28 顆白矮主序雙星[56]和230 顆DA 型白矮星[57,58];DR1[59]中的白矮主序雙星數(shù)量上升至121 顆[60];DR2發(fā)布的近1.59×106條光譜數(shù)據(jù)中包含1 056 顆DA 型白矮星、34 顆DB 型白矮星和276顆白矮主序雙星[61];最新的DR5 發(fā)布數(shù)據(jù)中,白矮主序雙星的數(shù)量為876 顆[62],DB 型白矮星的數(shù)量則為287 顆[63]。

    3.2.4 GAIA

    GAIA[64,65]依賴于歐洲航天局的依巴谷衛(wèi)星來幫助解決現(xiàn)代天文學(xué)中最困難,也是最根本的挑戰(zhàn)之一:在整個(gè)銀河系鄰域內(nèi)構(gòu)造一個(gè)精確的恒星三維地圖。這次大規(guī)模的恒星普查將會(huì)提供基本的觀測數(shù)據(jù),以解決與銀河系的起源、結(jié)構(gòu)和進(jìn)化歷史相關(guān)的大量重要問題。GAIA 衛(wèi)星于2013年12月19日升空,在L2 拉格朗日點(diǎn)觀測,并提供了不管是從數(shù)量上,還是在準(zhǔn)確性和完整性上都前所未有的數(shù)據(jù)檔案。

    截至目前,GAIA 已進(jìn)行了2 次數(shù)據(jù)發(fā)布,其中GAIA DR2[66,67]的一部分?jǐn)?shù)據(jù)覆蓋了G 波段(3 300~11 000),另外的數(shù)據(jù)分別覆蓋GBP(3 300~6 700)和GRP(6 100~10 500 ?A)2 個(gè)波段[68]。它所發(fā)布的星表中包含源的位置坐標(biāo)[69]、星等、視差[70]、自行和視向速度等參數(shù),具體如下。

    (1)有超過109個(gè)源擁有赤經(jīng)(RA)、赤緯(DEC)、視差、赤經(jīng)自行和赤緯自行這5 個(gè)參數(shù),這些源的極限星等為Gmag=21 mag。對(duì)于Gmag<15 mag 的亮源,視差的誤差只有0.04 mas;對(duì)于Gmag=17 mag 的源,誤差大約是0.1 mas;對(duì)于Gmag=20 mag 的源,誤差約為0.7 mas。對(duì)應(yīng)的自行誤差分別是0.06 mas·a?1(Gmag<15 mag),0.2 mas·a?1(Gmag=17 mag)和1.2 mas·a?1(Gmag=20 mag)。

    (2)超過7.2×106顆恒星擁有平均視向速度[71,72],它們的G 波段星等在4~13 mag范圍,有效溫度Teff為3 550~6 900 K。包括上述5 個(gè)參數(shù),這些星就擁有了全部的6 個(gè)參數(shù)。對(duì)于亮星,視向速度的精度在200~300 m·s?1量級(jí);對(duì)于較暗的星,當(dāng)有效溫度在4 750 K 和6 500 K 時(shí),視向速度的精度可以分別達(dá)到1.2 km·s?1和2.5 km·s?1。

    (3)超過3×108個(gè)源擁有2 個(gè)參數(shù):天球上的赤經(jīng)和赤緯坐標(biāo)。這些源在Gmag=20 mag 時(shí)的坐標(biāo)誤差大約為2 mas。

    (4)超過1.69×109個(gè)源擁有G 波段的星等,亮星(Gmag<13 mag)的星等精度可達(dá)0.001 mag;暗星(Gmag=20 mag)的星等精度大約為0.02 mag。

    (5)超過1.38×109個(gè)源擁有GBP和GRP波段的星等GBPmag和GRPmag,亮星和暗星的精度分別是10?3mag 量級(jí)和大約0.2 mag。

    GAIA DR1 和GAIA DR2 的發(fā)布數(shù)量如表3所示。

    根據(jù)GAIA 的視差和測光,Hollands 等人[74]對(duì)GAIA DR2 的白矮星進(jìn)行了大氣參數(shù)測量。對(duì)于提供了精確參數(shù)的源,他們重新定義了半徑為20 pc 范圍內(nèi)所有的白矮星雙星。

    表3 GAIA 星表概覽

    3.3 巡天項(xiàng)目中的白矮星搜尋

    在巡天項(xiàng)目釋放的海量(光譜)數(shù)據(jù)中,官方提供的星表已包括了大量白矮星數(shù)據(jù),如SDSS DR14 和LAMOST DR5 中共有上千萬光譜數(shù)據(jù),其中包含了數(shù)萬條白矮星光譜(見表4)。這些分類結(jié)果都是使用Pipeline 軟件[59,75]對(duì)光譜進(jìn)行模板匹配得到的,這種分類法可以在統(tǒng)計(jì)意義上保證分類的準(zhǔn)確性。但是對(duì)于信噪比較低的光譜數(shù)據(jù),模板匹配可能會(huì)出現(xiàn)匹配錯(cuò)誤的情況。所以有很多科研工作者對(duì)白矮星的分類結(jié)果進(jìn)行了完善,并積累了豐富的研究成果。

    表4 SDSS DR14 和LAMOST DR5 發(fā)布星表中的白矮星數(shù)量

    首先,如同第2.3 節(jié)所述,人們對(duì)白矮星的初步篩選方法是采用顏色判據(jù)。針對(duì)SDSS和LAMOST 的光譜數(shù)據(jù),Kleinman 等人[21]、Kepler 等人[39,41]和Guo 等人[61]篩選白矮星時(shí),都使用了先根據(jù)顏色的分布[36]來篩選出白矮星候選體,然后利用白矮星大氣模型[37]來擬合,或者用其他方法來證認(rèn)白矮星。

    除了傳統(tǒng)的顏色判據(jù),也有一些篩選工作是通過機(jī)器學(xué)習(xí)的方法來自動(dòng)挖掘白矮星。Kong 等人[20]利用LASSO 方法[76]提取出了DB 型白矮星光譜區(qū)別于其他類型光譜的典型特征,然后在特征空間中,使用支持向量基[77]對(duì)數(shù)千萬已發(fā)布數(shù)據(jù)進(jìn)行分類,從而把混雜在各種光譜類型中的DB 型白矮星光譜挑選出來。這種方法充分利用了機(jī)器學(xué)習(xí)的優(yōu)勢,尤其對(duì)于信噪比較低,肉眼難辨的數(shù)據(jù),用該方法可以將顯示出HeⅠ譜線的光譜數(shù)據(jù)快速挑出。這種方法也可用于其他稀有類型光譜數(shù)據(jù)的搜尋中。

    4 物理參數(shù)

    4.1 大氣參數(shù)

    Kepler 等人[39]在SDSS DR12 光譜數(shù)據(jù)中搜尋白矮星的同時(shí),也利用修正的三維對(duì)流模型[78]測量了高信噪比DA,DQ,DZ 型白矮星的大氣參數(shù)Teff和lgg。這些白矮星的Teff范圍為5 000~100 000 K。隨后,Koester[37]利用其提出的大氣模型,對(duì)SDSS 數(shù)據(jù)中DB 型白矮星的大氣參數(shù)也進(jìn)行了測量和分析[20,42]。借助于GAIA 的紫外數(shù)據(jù),Kong 等人[20]驗(yàn)證了DB 型白矮星的有效溫度與紫外波段顏色的線性關(guān)系,確定其斜率約為10?4。

    4.2 質(zhì)量

    白矮星的前身星質(zhì)量最大約為9M⊙,其空間分布包含著恒星形成歷史及其在銀河系中演化的信息。由于任何白矮星都沒有足夠的時(shí)間演化,以致冷卻到不可探測的程度,所以人們可從它們的觀測數(shù)據(jù)獲取其各自獨(dú)立的關(guān)于銀河系年齡的信息。通過始末質(zhì)量關(guān)系(initial-final mass relation,IFMR),我們也可以研究星體演化過程中的質(zhì)量損失情況。現(xiàn)今的大規(guī)模光譜巡天項(xiàng)目給我們提供了一個(gè)能夠充分研究白矮星質(zhì)量分布的全新視角。此外,了解白矮星的質(zhì)量分布,也對(duì)研究恒星演化中的質(zhì)量損失過程、IFMR 或緊致雙星演化有重要幫助。

    最常見的DA 型白矮星的質(zhì)量分布已得到了廣泛的研究,特別是對(duì)于Teff10 000 K的DA 型白矮星,其質(zhì)量可以通過光譜中巴爾末線的擬合得出[79,80]。通過對(duì)SDSS DR4 的白矮星樣本數(shù)據(jù)的質(zhì)量分布進(jìn)行分析,Kepler 等人[80]得出,表面溫度較低的DA 型白矮星(Teff12 000 K)和DB 型白矮星(Teff16 000 K)的質(zhì)量較大。這個(gè)結(jié)論依賴于現(xiàn)有的樣本,到目前為止仍然無法推翻,除非現(xiàn)有的模型是不準(zhǔn)確的。另外,DA 型白矮星和DB 型白矮星的質(zhì)量分布是不相同的,DB 型白矮星質(zhì)量的峰值更偏向于大質(zhì)量星。這個(gè)結(jié)論與以往的普遍認(rèn)知相悖,人們?cè)?jīng)認(rèn)為這兩種類型的白矮星質(zhì)量接近。最有可能的原因還是由于樣本的不足導(dǎo)致結(jié)果有偏差。以往可供研究的DB 型白矮星樣本不足50 個(gè),到SDSS DR4時(shí)也只有150 個(gè),而目前的數(shù)量已達(dá)上千個(gè)[20]。未來可能會(huì)有更準(zhǔn)確的結(jié)論。

    白矮星的質(zhì)量-半徑關(guān)系經(jīng)常受限于距離的不確定性。大多數(shù)研究都以巴爾末線為主,因?yàn)檫@些波段的光譜可以通過地基望遠(yuǎn)鏡的觀測得到,而萊曼線系光譜則只能通過空間紫外望遠(yuǎn)鏡的觀測來獲取。Joyce 等人[81]利用來自GAIA DR2 的自行和哈勃空間望遠(yuǎn)鏡的數(shù)據(jù)進(jìn)行研究,他們的樣本支持理論的質(zhì)量-半徑關(guān)系,雖然仍有至少一顆星的情況與理論并不一致。巴爾末和萊曼線系光譜的研究結(jié)果表明理論與觀測是一致的。

    Kepler[85]通過質(zhì)量-半徑關(guān)系[82–84],使用擬合得到的Teff和lgg,估算了SDSS 中所有高信噪比(S/N15)白矮星的質(zhì)量,并通過蒙特卡羅模擬推算出了理論上的質(zhì)量分布。DA 和DB 型白矮星的質(zhì)量分布形狀不同,DA 型白矮星的質(zhì)量分布在大質(zhì)量處有一個(gè)較長的拖尾,然而DB 型白矮星的質(zhì)量分布則在小質(zhì)量處有一個(gè)較長的拖尾。這種現(xiàn)象可能反映了前身星經(jīng)歷熱力學(xué)脈沖進(jìn)而演化成DB 型白矮星的一些限制。通過計(jì)算得出,所有誕生于1.2×1010a 之前的白矮星的有效溫度Teff<10 000 K,且大多數(shù)的溫度都在4 500 K 左右。對(duì)于Teff13 000 K 的白矮星,平均質(zhì)量若在0.64M⊙左右,則它的前身星在2×109a 前一般會(huì)爆發(fā),進(jìn)而形成白矮星,因?yàn)槠淝吧硇琴|(zhì)量只有約2.5M⊙,并且壽命較短。這些結(jié)果表明,大多數(shù)白矮星比我們現(xiàn)在觀測到的樣本更冷,我們必須想辦法將它們與紅矮星區(qū)分開。這種He 白矮星在光譜中沒有顯示出譜線,而DA 型僅顯示出弱的Hα 線。因此,自行和視差是關(guān)鍵,因?yàn)樵诮频臏囟认?,冷白矮星比紅矮星更暗。因此,GAIA 的數(shù)據(jù)對(duì)于它們的識(shí)別至關(guān)重要。

    4.3 運(yùn)動(dòng)學(xué)

    白矮星的運(yùn)動(dòng)學(xué)研究可以為銀河系暗物質(zhì)研究提供非常大的幫助[86,87]。通過對(duì)年輕的DA 型白矮星進(jìn)行分析,Wegg 和Phinney[88]發(fā)現(xiàn)小質(zhì)量的DA 型白矮星(小于0.75M⊙)擁有較高的速度彌散(約46 km·s?1)和非對(duì)稱星流,大質(zhì)量的DA 型白矮星(大于0.75M⊙)速度彌散僅為約19 km·s?1。這應(yīng)該是由于較大質(zhì)量的前身星,其壽命較短導(dǎo)致的。通過對(duì)比SDSS DR14 與LAMOST DR5 中DB 型白矮星的質(zhì)量和空間速度,Kong 等人[63]也發(fā)現(xiàn)了其與DA 型白矮星相似的規(guī)律。

    另外,根據(jù)樣本計(jì)算白矮星的形成率時(shí),白矮星距離銀盤的高度及其隨年齡和質(zhì)量的變化是至關(guān)重要的。

    5 形成和演化

    白矮星是小質(zhì)量(不大于8M⊙)恒星的核心殘留物,由其前身星經(jīng)歷后主序演化之后形成[89]。根據(jù)理論模型[90]計(jì)算,最小質(zhì)量白矮星是通過一顆質(zhì)量為0.50M⊙~0.52M⊙的恒星演化而來,因?yàn)橘|(zhì)量更小的恒星,其壽命要比宇宙年齡更長。根據(jù)模型預(yù)言,這種白矮星應(yīng)該為C/O 白矮星。然而,通過對(duì)SDSS 發(fā)布數(shù)據(jù)中的白矮-主序雙星的研究,Rebassa-Mansergas 等人[91]發(fā)現(xiàn)雙星演化可以改變這種結(jié)果,產(chǎn)生不同類型的白矮星。

    當(dāng)核心的H 和He 在演化早期燃燒殆盡后,質(zhì)量介于0.4M⊙~1.05M⊙的白矮星會(huì)擁有C/O 的核球[92]。質(zhì)量大于1.05M⊙[93]的白矮星的前身星則具有足夠高的溫度來點(diǎn)燃C 核,進(jìn)而形成O/Ne 或Ne/O/Mg 核。具體燃燒到何種程度取決于C 燃燒是否在星體的中央。當(dāng)C 燃燒在偏離星核中心的地方開始后,“火焰”雖然會(huì)向中心傳播,但最終并不一定會(huì)到達(dá)中心。如果這種“C 火”沒有達(dá)到星體核心,那么,C/O 和O/Ne 混合的白矮星就誕生了[94,95],這種情況下核質(zhì)量的下限是1.04M⊙;如果抵達(dá)了中心,則會(huì)產(chǎn)生一顆Ne/O/Mg白矮星,其最小的核質(zhì)量則是1.154M⊙[96]。這些結(jié)論都是依賴于演化模型得出的[93]。不論何種核的白矮星,質(zhì)量的確定對(duì)于我們理解白矮星的形成和演化都是至關(guān)重要的前提條件。

    已有人提出了中小質(zhì)量C/O 核白矮星的演化模型[97,98],也有很多人對(duì)極端小質(zhì)量白矮星的演化模型進(jìn)行了研究[99–101]。但是,目前仍然缺少完善的大質(zhì)量白矮星演化模型。

    隨著巡天項(xiàng)目的開展,光譜證認(rèn)的白矮星數(shù)量已超過4×104顆?;谶@些數(shù)據(jù),Kleinman 等人[34]已對(duì)它們的質(zhì)量分布作了較為細(xì)致的研究和分析。隨著觀測樣本的不斷完備,現(xiàn)有的模型可能會(huì)得到進(jìn)一步的修正和完善,所缺少的模型也會(huì)逐漸被建立起來。

    6 總結(jié)與展望

    白矮星作為銀河系演化的示蹤體,一直以來都備受天文學(xué)家的青睞。它是超過97% 的恒星演化的最終階段。由于星體中沒有了核聚變,僅僅依靠引力收縮輻射來提供能量,所以其亮度非常暗,以致其較難被觀測到。

    根據(jù)大氣層中元素的含量,白矮星可分為DA,DB,DO 等多種光譜型。由于星體內(nèi)部物質(zhì)與表層的對(duì)流,星體表面的元素含量會(huì)發(fā)生變化,從而使得白矮星的類型也會(huì)相互轉(zhuǎn)化。DB 型白矮星的Teff在30 000~45 000 K 范圍內(nèi)的DB 缺失并不是絕對(duì)的,而只是其數(shù)量較少。

    現(xiàn)有H 和He 白矮星樣本的質(zhì)量分布并不完全相同,它們的平均質(zhì)量都在0.5M⊙左右,但有一部分DA 型白矮星的質(zhì)量較大,而DB 型白矮星卻在較小質(zhì)量處有拖尾。另外,溫度較低的白矮星通常具有較大的質(zhì)量,而且小質(zhì)量白矮星的空間運(yùn)動(dòng)速度彌散比大質(zhì)量的大。

    得益于近十幾年來的大型巡天望遠(yuǎn)鏡觀測項(xiàng)目,獲取大量白矮星光譜數(shù)據(jù)成為可能。比較有代表性的是SDSS 巡天項(xiàng)目:從2000年至今的14 次發(fā)布數(shù)據(jù)中,已經(jīng)有超過4×104顆白矮星的光譜。Kepler,Kleinman,Koester 等人利用這些數(shù)據(jù)對(duì)白矮星進(jìn)行了大量系統(tǒng)性的搜尋和研究工作,獲得了關(guān)于其質(zhì)量分布、光度函數(shù)、元素豐度等一系列研究成果。他們?cè)诤A抗庾V數(shù)據(jù)中搜尋白矮星的通用做法就是限制u′?g′和g′?r′的顏色范圍,因?yàn)榘装堑念伾采w范圍與主序星有著較明顯的區(qū)別。另外,Kong 等人使用機(jī)器學(xué)習(xí)的方法在SDSS DR14 和LAMOST DR5 的海量光譜中把DB 型白矮星的數(shù)量增加到了2 000 多顆。

    LAMOST DR5 和SDSS DR14 的數(shù)據(jù)足跡基本上涵蓋了北天球的絕大多數(shù)區(qū)域,兩者同源的數(shù)據(jù)也有很多。所不同的是,SDSS 數(shù)據(jù)中白矮星的數(shù)量比LAMOST 的要多一些,這主要是望遠(yuǎn)鏡的極限星等和選源策略不同導(dǎo)致的。GALEX 巡天也發(fā)布了很多白矮星的紫外數(shù)據(jù),這對(duì)白矮星的研究工作也有很大的幫助。

    隨著各巡天項(xiàng)目的不斷開展,人們可以獲得海量的白矮星數(shù)據(jù)。充分利用好機(jī)器學(xué)習(xí)的優(yōu)勢,可以把搜尋DB 型白矮星的方法應(yīng)用到搜尋其他類型的白矮星上。

    隨著觀測設(shè)備的升級(jí)和觀測精度的提高,人類有望觀測到越來越多的白矮星。研究白矮星對(duì)于理解銀河系中的恒星形成和演化有著非常重要的意義。

    猜你喜歡
    白矮星巡天恒星
    中國的巡天利器有多強(qiáng)
    崔向群:巡天望遠(yuǎn)追星人
    (18)刺殺恒星
    罕見的白矮星和身旁行星
    飛碟探索(2020年2期)2020-10-09 10:49:50
    白矮星在說話
    奧秘(2019年5期)2019-05-14 11:40:24
    DA白矮星視向速度測量*
    恒星的演化
    恒星不恒
    奧秘(2018年10期)2018-10-25 05:38:56
    “巡天·觀地·測月”夢想實(shí)現(xiàn)——“嫦娥三號(hào)”最新科研成果掃描
    嫦娥三號(hào)進(jìn)入第33月夜“測月、巡天、觀地”成果豐碩
    太空探索(2016年9期)2016-07-12 10:00:00
    欧美亚洲日本最大视频资源| 韩国精品一区二区三区| 中国美女看黄片| 日本av免费视频播放| 欧美精品一区二区免费开放| 欧美日韩视频精品一区| 操出白浆在线播放| 国产欧美日韩综合在线一区二区| 十分钟在线观看高清视频www| 亚洲欧洲日产国产| 日韩 亚洲 欧美在线| 搡老岳熟女国产| 俄罗斯特黄特色一大片| 一区二区av电影网| av天堂久久9| 色老头精品视频在线观看| 国产亚洲精品第一综合不卡| 男女之事视频高清在线观看| 亚洲国产欧美在线一区| 超碰成人久久| 午夜福利乱码中文字幕| 狠狠狠狠99中文字幕| 亚洲免费av在线视频| 真人做人爱边吃奶动态| 亚洲伊人色综图| 欧美日韩中文字幕国产精品一区二区三区 | 精品国产国语对白av| 老司机靠b影院| 蜜桃在线观看..| 少妇裸体淫交视频免费看高清 | 国产精品久久久久久精品电影小说| 无遮挡黄片免费观看| 久久国产精品大桥未久av| 国产熟女午夜一区二区三区| 午夜福利视频精品| 丝瓜视频免费看黄片| av有码第一页| e午夜精品久久久久久久| 亚洲成人手机| 欧美av亚洲av综合av国产av| av网站免费在线观看视频| 高清av免费在线| 老司机亚洲免费影院| 欧美精品高潮呻吟av久久| videosex国产| 国产精品亚洲av一区麻豆| 女人爽到高潮嗷嗷叫在线视频| 国产精品影院久久| 天堂中文最新版在线下载| 国产成人一区二区三区免费视频网站| 91九色精品人成在线观看| 99国产极品粉嫩在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 午夜福利在线观看吧| 国产成人啪精品午夜网站| 色婷婷久久久亚洲欧美| 一区二区日韩欧美中文字幕| 免费看十八禁软件| 搡老乐熟女国产| 秋霞在线观看毛片| 国产在线视频一区二区| 亚洲国产av影院在线观看| 男男h啪啪无遮挡| 自线自在国产av| 午夜福利在线免费观看网站| 老司机午夜福利在线观看视频 | 国产精品一二三区在线看| kizo精华| 18禁黄网站禁片午夜丰满| 久久国产亚洲av麻豆专区| 三上悠亚av全集在线观看| 欧美激情极品国产一区二区三区| 久久久久久亚洲精品国产蜜桃av| 在线 av 中文字幕| 精品乱码久久久久久99久播| 国产极品粉嫩免费观看在线| 丰满少妇做爰视频| 久热这里只有精品99| av超薄肉色丝袜交足视频| 一级毛片电影观看| 午夜精品国产一区二区电影| 性色av乱码一区二区三区2| 亚洲精品乱久久久久久| 亚洲专区国产一区二区| 久久久久精品人妻al黑| 精品一区二区三区四区五区乱码| www日本在线高清视频| 久久精品熟女亚洲av麻豆精品| 久久久久久久久久久久大奶| av国产精品久久久久影院| 三上悠亚av全集在线观看| 90打野战视频偷拍视频| 亚洲精品国产一区二区精华液| av免费在线观看网站| 久久久欧美国产精品| a在线观看视频网站| 精品第一国产精品| 伊人久久大香线蕉亚洲五| 精品久久久久久久毛片微露脸 | 成人av一区二区三区在线看 | 欧美国产精品va在线观看不卡| 久久狼人影院| 欧美精品一区二区免费开放| 天天躁日日躁夜夜躁夜夜| 国产精品麻豆人妻色哟哟久久| 久久99一区二区三区| 黑人猛操日本美女一级片| 亚洲欧美激情在线| 日韩视频一区二区在线观看| 久久久国产一区二区| 黄色视频在线播放观看不卡| 精品人妻熟女毛片av久久网站| 国产极品粉嫩免费观看在线| 一边摸一边抽搐一进一出视频| 人妻 亚洲 视频| 超色免费av| 脱女人内裤的视频| 天天影视国产精品| 久久九九热精品免费| 亚洲精品av麻豆狂野| 亚洲av成人不卡在线观看播放网 | 国产亚洲av高清不卡| 久久天堂一区二区三区四区| 后天国语完整版免费观看| 中文字幕高清在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 汤姆久久久久久久影院中文字幕| 脱女人内裤的视频| 精品国产一区二区三区四区第35| 亚洲激情五月婷婷啪啪| 热re99久久国产66热| 国产精品麻豆人妻色哟哟久久| av超薄肉色丝袜交足视频| 国产精品秋霞免费鲁丝片| 亚洲 欧美一区二区三区| 亚洲欧美日韩另类电影网站| 男女无遮挡免费网站观看| 欧美亚洲 丝袜 人妻 在线| 日日夜夜操网爽| 欧美+亚洲+日韩+国产| 91av网站免费观看| 我的亚洲天堂| 天天影视国产精品| 久久香蕉激情| 亚洲午夜精品一区,二区,三区| kizo精华| 国产精品亚洲av一区麻豆| 欧美日韩精品网址| 性高湖久久久久久久久免费观看| 欧美中文综合在线视频| 日韩制服丝袜自拍偷拍| 亚洲激情五月婷婷啪啪| 久久精品国产综合久久久| 国产xxxxx性猛交| 亚洲人成电影观看| 王馨瑶露胸无遮挡在线观看| 欧美老熟妇乱子伦牲交| 亚洲国产av影院在线观看| 岛国在线观看网站| 国产三级黄色录像| 精品一区在线观看国产| 老熟妇仑乱视频hdxx| 99精国产麻豆久久婷婷| 亚洲国产精品一区三区| 在线观看人妻少妇| 亚洲精品久久久久久婷婷小说| 国产成人免费无遮挡视频| 国产精品麻豆人妻色哟哟久久| 又大又爽又粗| 男男h啪啪无遮挡| 国产淫语在线视频| 国产国语露脸激情在线看| 国产精品久久久久成人av| 欧美午夜高清在线| 日韩大片免费观看网站| 婷婷成人精品国产| 欧美精品人与动牲交sv欧美| 脱女人内裤的视频| 午夜福利一区二区在线看| 国产欧美日韩一区二区精品| 男人爽女人下面视频在线观看| 深夜精品福利| 色综合欧美亚洲国产小说| 亚洲中文字幕日韩| 久久精品国产亚洲av高清一级| 久久午夜综合久久蜜桃| 日本欧美视频一区| 国产国语露脸激情在线看| a 毛片基地| 久久久久精品人妻al黑| 欧美少妇被猛烈插入视频| 18禁观看日本| 十八禁人妻一区二区| 国产成人精品久久二区二区免费| 久久av网站| 女人精品久久久久毛片| 老汉色∧v一级毛片| 狠狠婷婷综合久久久久久88av| 搡老乐熟女国产| 99国产极品粉嫩在线观看| 亚洲精品自拍成人| 老司机深夜福利视频在线观看 | 欧美乱码精品一区二区三区| 99re6热这里在线精品视频| 夜夜夜夜夜久久久久| 国产亚洲av高清不卡| 久久国产精品大桥未久av| 亚洲欧洲精品一区二区精品久久久| 午夜精品国产一区二区电影| 一级片免费观看大全| 成人影院久久| 久热爱精品视频在线9| 美女视频免费永久观看网站| √禁漫天堂资源中文www| 丰满少妇做爰视频| 12—13女人毛片做爰片一| 日韩熟女老妇一区二区性免费视频| 一区二区av电影网| 久久精品国产综合久久久| 国产成人系列免费观看| 久久人妻福利社区极品人妻图片| 久久久久国产一级毛片高清牌| 午夜福利影视在线免费观看| 在线亚洲精品国产二区图片欧美| 国产在线视频一区二区| 久久精品久久久久久噜噜老黄| 99久久国产精品久久久| 国产在视频线精品| 久久狼人影院| 一区在线观看完整版| 一区二区三区乱码不卡18| 国产亚洲av片在线观看秒播厂| 欧美激情高清一区二区三区| 亚洲avbb在线观看| 老司机影院成人| 伊人久久大香线蕉亚洲五| 国产伦理片在线播放av一区| 精品国产乱码久久久久久男人| 国产熟女午夜一区二区三区| xxxhd国产人妻xxx| 国产欧美日韩一区二区三区在线| 19禁男女啪啪无遮挡网站| 国产极品粉嫩免费观看在线| 精品卡一卡二卡四卡免费| 日本五十路高清| 天堂中文最新版在线下载| 一级片'在线观看视频| a级毛片黄视频| 久久久欧美国产精品| 欧美成狂野欧美在线观看| 日日摸夜夜添夜夜添小说| 日韩欧美国产一区二区入口| 久久中文字幕一级| 在线av久久热| 亚洲欧美清纯卡通| 亚洲少妇的诱惑av| 亚洲精品美女久久av网站| 成年人黄色毛片网站| 中国国产av一级| 欧美日韩av久久| 80岁老熟妇乱子伦牲交| 亚洲激情五月婷婷啪啪| 久久热在线av| 亚洲国产欧美日韩在线播放| 19禁男女啪啪无遮挡网站| 午夜老司机福利片| 一本综合久久免费| 日韩 欧美 亚洲 中文字幕| 亚洲精品自拍成人| 99国产精品一区二区三区| tube8黄色片| av超薄肉色丝袜交足视频| 久久热在线av| av线在线观看网站| 人成视频在线观看免费观看| 一级毛片电影观看| 久久人人爽人人片av| www.999成人在线观看| 中文精品一卡2卡3卡4更新| www.自偷自拍.com| 在线观看免费午夜福利视频| 欧美亚洲日本最大视频资源| 久久免费观看电影| 免费在线观看视频国产中文字幕亚洲 | 少妇裸体淫交视频免费看高清 | 亚洲成av片中文字幕在线观看| 亚洲精品国产一区二区精华液| 精品人妻一区二区三区麻豆| 亚洲国产av影院在线观看| 美女视频免费永久观看网站| 香蕉国产在线看| 欧美精品亚洲一区二区| 美女大奶头黄色视频| 丝袜美足系列| 精品久久久久久久毛片微露脸 | 国产欧美日韩一区二区三 | 亚洲伊人色综图| 欧美日韩视频精品一区| 欧美日韩黄片免| a级毛片在线看网站| av不卡在线播放| 90打野战视频偷拍视频| 国产成人精品久久二区二区免费| 亚洲精品国产精品久久久不卡| 久久国产精品人妻蜜桃| 国产不卡av网站在线观看| 成人三级做爰电影| 男女无遮挡免费网站观看| 国产97色在线日韩免费| 涩涩av久久男人的天堂| 少妇人妻久久综合中文| 精品亚洲乱码少妇综合久久| 精品久久久久久电影网| 午夜免费成人在线视频| 色婷婷av一区二区三区视频| 97在线人人人人妻| 最新在线观看一区二区三区| 新久久久久国产一级毛片| 久久九九热精品免费| 黑人巨大精品欧美一区二区mp4| 最近最新中文字幕大全免费视频| 国产老妇伦熟女老妇高清| 欧美日韩亚洲综合一区二区三区_| 在线观看免费高清a一片| 亚洲成av片中文字幕在线观看| 搡老乐熟女国产| 午夜免费观看性视频| 老司机午夜福利在线观看视频 | 亚洲视频免费观看视频| 国产高清videossex| 国产主播在线观看一区二区| 日本猛色少妇xxxxx猛交久久| 欧美精品高潮呻吟av久久| 桃花免费在线播放| 免费黄频网站在线观看国产| 久久久久视频综合| 国产真人三级小视频在线观看| 一本—道久久a久久精品蜜桃钙片| 99久久人妻综合| bbb黄色大片| 丝袜在线中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 三上悠亚av全集在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 12—13女人毛片做爰片一| 国产一区二区三区av在线| 一边摸一边做爽爽视频免费| 成人手机av| 天堂俺去俺来也www色官网| 久久精品成人免费网站| 亚洲欧美激情在线| 日本av免费视频播放| 动漫黄色视频在线观看| 欧美人与性动交α欧美精品济南到| 超碰97精品在线观看| 黄片大片在线免费观看| 亚洲精品中文字幕在线视频| 99香蕉大伊视频| 黄色 视频免费看| 国产精品成人在线| 少妇粗大呻吟视频| 国产成人av教育| 三上悠亚av全集在线观看| 中文字幕人妻熟女乱码| 好男人电影高清在线观看| 免费在线观看黄色视频的| 欧美日韩亚洲高清精品| 午夜精品久久久久久毛片777| 老司机靠b影院| 久久久国产成人免费| 12—13女人毛片做爰片一| 成人黄色视频免费在线看| 欧美日韩黄片免| 亚洲熟女精品中文字幕| 男人舔女人的私密视频| 国产一区二区在线观看av| 精品亚洲乱码少妇综合久久| 成年女人毛片免费观看观看9 | 欧美日韩av久久| 电影成人av| 国产精品一区二区精品视频观看| 国产成人影院久久av| 国产视频一区二区在线看| 午夜福利在线观看吧| 亚洲五月色婷婷综合| 这个男人来自地球电影免费观看| 不卡一级毛片| 一进一出抽搐动态| 精品亚洲成a人片在线观看| 成人国产av品久久久| 亚洲成av片中文字幕在线观看| 人人妻人人添人人爽欧美一区卜| 亚洲精品粉嫩美女一区| 欧美日韩国产mv在线观看视频| 最近中文字幕2019免费版| 夜夜骑夜夜射夜夜干| 狠狠婷婷综合久久久久久88av| 日韩免费高清中文字幕av| 精品免费久久久久久久清纯 | 国产在视频线精品| 亚洲黑人精品在线| 午夜两性在线视频| 别揉我奶头~嗯~啊~动态视频 | 国产亚洲欧美精品永久| 考比视频在线观看| 黄色怎么调成土黄色| 亚洲精品国产精品久久久不卡| 国产成人a∨麻豆精品| 大陆偷拍与自拍| 无限看片的www在线观看| 男女之事视频高清在线观看| 国产精品一区二区精品视频观看| 成人手机av| 午夜视频精品福利| 亚洲av男天堂| 精品视频人人做人人爽| 9色porny在线观看| 人妻 亚洲 视频| 可以免费在线观看a视频的电影网站| 久久久久久免费高清国产稀缺| 国产精品.久久久| 熟女少妇亚洲综合色aaa.| 久久国产精品影院| 免费女性裸体啪啪无遮挡网站| 日日夜夜操网爽| 欧美午夜高清在线| 午夜激情久久久久久久| 日本五十路高清| 美女主播在线视频| 亚洲欧美激情在线| 母亲3免费完整高清在线观看| 美国免费a级毛片| www.av在线官网国产| 中国国产av一级| 午夜福利乱码中文字幕| 精品人妻在线不人妻| 亚洲综合色网址| 国产精品国产av在线观看| 成人影院久久| 少妇被粗大的猛进出69影院| 免费久久久久久久精品成人欧美视频| 亚洲视频免费观看视频| 久久热在线av| 亚洲精品一二三| 国产成人啪精品午夜网站| 蜜桃在线观看..| 啦啦啦中文免费视频观看日本| 国产一区二区三区在线臀色熟女 | 1024视频免费在线观看| 欧美 日韩 精品 国产| 性少妇av在线| 一区二区日韩欧美中文字幕| 亚洲欧美日韩高清在线视频 | av片东京热男人的天堂| 欧美黑人精品巨大| 天堂8中文在线网| 成年人免费黄色播放视频| 丝袜美足系列| 亚洲av欧美aⅴ国产| 精品国产超薄肉色丝袜足j| 亚洲欧美成人综合另类久久久| 亚洲欧美色中文字幕在线| av一本久久久久| 日韩三级视频一区二区三区| 老司机深夜福利视频在线观看 | 老熟妇仑乱视频hdxx| 十分钟在线观看高清视频www| 久久亚洲精品不卡| 高潮久久久久久久久久久不卡| 五月开心婷婷网| 两性午夜刺激爽爽歪歪视频在线观看 | 91九色精品人成在线观看| 大型av网站在线播放| 黑人巨大精品欧美一区二区mp4| 99香蕉大伊视频| 午夜福利视频在线观看免费| 成人影院久久| 一个人免费在线观看的高清视频 | 国产一区二区激情短视频 | 久久精品亚洲av国产电影网| 日本五十路高清| 欧美中文综合在线视频| 啦啦啦 在线观看视频| 中文字幕精品免费在线观看视频| 免费高清在线观看视频在线观看| av线在线观看网站| √禁漫天堂资源中文www| 国产欧美日韩精品亚洲av| 人妻人人澡人人爽人人| 1024视频免费在线观看| 中文字幕制服av| 在线av久久热| 精品乱码久久久久久99久播| 在线观看舔阴道视频| 日韩人妻精品一区2区三区| 欧美日韩亚洲高清精品| 中文精品一卡2卡3卡4更新| 亚洲av美国av| 人妻人人澡人人爽人人| 欧美日韩成人在线一区二区| 精品人妻一区二区三区麻豆| 亚洲精品国产色婷婷电影| 岛国毛片在线播放| 99国产精品一区二区三区| 久久精品亚洲av国产电影网| 亚洲精品一区蜜桃| 国产精品 欧美亚洲| 在线精品无人区一区二区三| 国产免费一区二区三区四区乱码| 午夜福利在线观看吧| 久久人人97超碰香蕉20202| 欧美日韩亚洲国产一区二区在线观看 | 精品第一国产精品| 日韩熟女老妇一区二区性免费视频| 在线观看舔阴道视频| 免费看十八禁软件| 精品人妻1区二区| 久久久久精品国产欧美久久久 | 精品少妇久久久久久888优播| 国内毛片毛片毛片毛片毛片| 一区二区三区激情视频| 日韩大片免费观看网站| 人人妻人人爽人人添夜夜欢视频| 在线观看舔阴道视频| 这个男人来自地球电影免费观看| 亚洲伊人色综图| 久久久久精品国产欧美久久久 | 高清在线国产一区| 五月开心婷婷网| 黄色视频在线播放观看不卡| 99久久国产精品久久久| 欧美变态另类bdsm刘玥| 在线观看免费高清a一片| 亚洲一码二码三码区别大吗| 黄色 视频免费看| 自拍欧美九色日韩亚洲蝌蚪91| 在线精品无人区一区二区三| 在线十欧美十亚洲十日本专区| 欧美成狂野欧美在线观看| 精品亚洲成国产av| 亚洲情色 制服丝袜| 黄色视频,在线免费观看| 黑人巨大精品欧美一区二区mp4| 三级毛片av免费| 午夜福利一区二区在线看| 亚洲色图综合在线观看| 大型av网站在线播放| 国产成人欧美| 亚洲国产中文字幕在线视频| 午夜福利免费观看在线| 国产精品一区二区精品视频观看| 日韩制服骚丝袜av| 国产日韩欧美在线精品| 一区二区三区四区激情视频| 欧美在线一区亚洲| 亚洲国产欧美日韩在线播放| 免费久久久久久久精品成人欧美视频| 欧美激情极品国产一区二区三区| 国产精品99久久99久久久不卡| 久久精品久久久久久噜噜老黄| 久久精品国产亚洲av香蕉五月 | 中国国产av一级| 俄罗斯特黄特色一大片| 另类精品久久| 国产区一区二久久| 欧美黄色片欧美黄色片| 99国产精品一区二区三区| 亚洲五月婷婷丁香| 啦啦啦免费观看视频1| 精品国内亚洲2022精品成人 | 亚洲精品久久成人aⅴ小说| av国产精品久久久久影院| 成年美女黄网站色视频大全免费| 999精品在线视频| 欧美一级毛片孕妇| 在线天堂中文资源库| 老熟妇乱子伦视频在线观看 | 性色av一级| 国产成人欧美在线观看 | 欧美午夜高清在线| 精品久久久精品久久久| 一区二区三区精品91| 欧美在线黄色| 国产99久久九九免费精品| 精品福利永久在线观看| 12—13女人毛片做爰片一| 欧美国产精品一级二级三级| 精品少妇黑人巨大在线播放| 9热在线视频观看99| 成人手机av| 成年动漫av网址| 久久青草综合色| 一级片免费观看大全| 久久精品国产a三级三级三级| 国产片内射在线| 国产精品.久久久| 老司机深夜福利视频在线观看 | 菩萨蛮人人尽说江南好唐韦庄| 国产精品熟女久久久久浪| 午夜福利免费观看在线| 亚洲精品久久午夜乱码| 免费高清在线观看日韩| 无限看片的www在线观看| 色精品久久人妻99蜜桃| 亚洲欧美激情在线| 日本一区二区免费在线视频| 在线观看舔阴道视频| 欧美黄色片欧美黄色片| 午夜福利影视在线免费观看| 色精品久久人妻99蜜桃| 91字幕亚洲| 高清av免费在线| 99精品欧美一区二区三区四区| 多毛熟女@视频| 欧美黄色片欧美黄色片|