楊 鵬,李 祥,羅志全,羅楊平

(西華師范大學(xué) 物理與天文學(xué)院,四川 南充 637009)

銀河系是我們的家園星系,理解其結(jié)構(gòu)和演化特點,可以有效地幫助我們了解宇宙中其他星系的結(jié)構(gòu)和演化。銀河系是一個盤狀星系,盤上包含了絕大部分重子物質(zhì),星系中的諸多演化活動,如恒星形成、旋臂、星系棒以及各種形式的長期演化,都發(fā)生在星系盤中[1]。恒星的位置和速度信息反過來又揭示了一些銀盤的內(nèi)部相互作用和外部擾動信息。星系震動學(xué)描述了銀河系在非對稱勢中非平衡和非平穩(wěn)的特征[2-4],具體表現(xiàn)為相空間中豐富的非對稱結(jié)構(gòu)分布[5-15]。隨著對銀河系的認(rèn)識不斷更新,它的一些細(xì)節(jié)結(jié)構(gòu)及演化特點逐漸被發(fā)現(xiàn)。

屋脊結(jié)構(gòu)描述的是恒星的徑向速度(VR)在銀心距(R)與方位角速度(Vφ)的平面分布中(或其他分布中),清晰地呈現(xiàn)出很多斜線結(jié)構(gòu)特征[11]。屋脊由Antoja等[16]在2018年發(fā)現(xiàn)或證實,他們利用Gaia DR2數(shù)據(jù)不僅揭示了太陽鄰域的屋脊結(jié)構(gòu)、拱形和蝸牛殼,還通過研究發(fā)現(xiàn)拱形可能是屋脊在VR-Vφ平面的投影,并且指出對于銀河系的建模需要考慮非對稱的引力勢。最近,Gaia DR3的發(fā)布使我們可以從更加精細(xì)的層面揭露銀河系的相關(guān)細(xì)節(jié)。Gaia DR3相較于Gaia DR1和Gaia DR2,在數(shù)據(jù)量和測量精度上都有很大的提升,總共提供了15億顆恒星的高精度位置、視差和自行的信息,其中3 300萬顆極限星等G≈14 mag的恒星帶有視向速度信息[17-19]。在Gaia DR3幫助下,Recio等[20]在繪制銀河系的化學(xué)分布圖像時,發(fā)現(xiàn)屋脊是一種富金屬的結(jié)構(gòu),且總共發(fā)現(xiàn)了7條屋脊,其中4條較為明顯,其他3條噪聲很大,不夠清晰。Khanna等[21]利用GALAH巡天數(shù)據(jù)在垂直距離(Z)、垂直速度(VZ)、徑向速度(VR)、金屬豐度(Fe/H)和元素豐度(α/Fe)分布中揭示了屋脊結(jié)構(gòu)。Wang等[11]利用LAMOST巡天數(shù)據(jù)中的主序拐點星和OB恒星,在化學(xué)動力學(xué)空間中研究了屋脊的分布,共發(fā)現(xiàn)了3條屋脊,其中2條相對穩(wěn)定,另1條隨著時間演化,但沒有進(jìn)一步在作用量空間中探索屋脊結(jié)構(gòu)。Antoja等[22]利用Gaia EDR3數(shù)據(jù)在銀河系外盤發(fā)現(xiàn)了1條屋脊,但是沒能在徑向速度分布中揭示屋脊的分布。Khoperskov等[23]發(fā)現(xiàn)角動量過密度區(qū)域中的恒星可以追蹤許多千秒差距上的密度屋脊以及太陽鄰域VR-Vφ平面上的過密度區(qū)域,而且當(dāng)它們穿過太陽鄰域時,這些屋脊也可以匹配主要的VR-Vφ特征。

為探究屋脊結(jié)構(gòu)可能的形成機制,人們采用了不同的模型。利用“toy”模型,Antoja等[16]發(fā)現(xiàn)屋脊結(jié)構(gòu)可以由水平方向的相位混合產(chǎn)生,但是沒能排除銀河系中心棒的影響。Hunt等[24]發(fā)現(xiàn),觀測到的屋脊結(jié)構(gòu)可以很好地在構(gòu)建的纏繞螺旋結(jié)構(gòu)的測試粒子模型中再現(xiàn)。Fragkoudi等[25]利用無碰撞N體模擬并結(jié)合軌道積分發(fā)現(xiàn)屋脊結(jié)構(gòu)是由銀河系中心棒的外部Lindblad共振引起的。Barros等[26]發(fā)現(xiàn)螺旋共振捕獲的恒星軌道也能產(chǎn)生屋脊特征?？偟膩碚f,屋脊的產(chǎn)生機制可以大致分為兩類:一類是沒有外部擾動的內(nèi)部動力學(xué)作用,如旋臂或/和銀河系中心棒的共振作用[26-28];另一類則是外部作用,如銀河系與人馬座星系之間的相互作用[16,29-30]。但是,最近Khanna等[21]暗示內(nèi)部和外部共同作用也同樣可以產(chǎn)生屋脊結(jié)構(gòu)。

雖然之前的工作從不同的角度對屋脊結(jié)構(gòu)進(jìn)行了探索,但是對于銀河系外盤的屋脊分布研究較少。本文利用LAMOST巡天數(shù)據(jù)得到的紅團(tuán)簇星對銀河系外盤的屋脊結(jié)構(gòu)進(jìn)行探索,并在作用量和角動量空間中繪制銀河系盤的屋脊結(jié)構(gòu)分布,從而促進(jìn)對銀河系盤動力學(xué)的理解。

1 數(shù) 據(jù)

1.1 數(shù)據(jù)介紹

本文使用的是Ting等[31]利用LAMOST巡天數(shù)據(jù)得到的紅團(tuán)簇星,并且這些紅團(tuán)簇星的距離不確定性在10%以內(nèi)。此星表中包含了所使用紅團(tuán)簇星的自行、視差、距離和視向速度等信息。采用來自Ting等[32]使用機器學(xué)習(xí)方法獲得的年齡信息,精度為25%。此外,還使用了Xiang等[33]利用LAMOST DR5低分辨率光譜結(jié)合Data-Driven Payne方法得到的8 162 566顆恒星的星表,該星表給出了本文中所使用紅團(tuán)簇星的金屬豐度和元素豐度信息。由于紅團(tuán)簇星處于核心氦燃燒階段,其光度基本保持不變,是銀河系中為數(shù)不多的標(biāo)準(zhǔn)燭光,這為后面得到的距離信息提供了可靠的基礎(chǔ)。

1.2 數(shù)據(jù)篩選

以上述數(shù)據(jù)為基礎(chǔ),通過交叉匹配的方式獲取了每顆紅團(tuán)簇星的位置、視差、自行、有效溫度、表面重力、金屬豐度和元素豐度信息。為了進(jìn)一步提高數(shù)據(jù)準(zhǔn)確性以及過濾信息不完整的數(shù)據(jù),采用以下條件進(jìn)行篩選:

1)8 kpc

2)50 km·s-1

3)Parallax >0;

4)0 Gyr

5)[Fe/H]>-1.5 dex;

6)-1.15<[α/Fe]<1.58。

經(jīng)過以上條件的篩選,一共獲得了118 535顆同時具有空間位置、速度和化學(xué)信息的紅團(tuán)簇星。圖1展示了最后樣本的分布情況。圖1(a)展示了樣本在Teff-logg平面的分布;圖1(b)展示了樣本在銀心距R=[8,15] kpc,垂直距離Z=[-1.5,1.5] kpc平面中的分布;圖1(c)則是樣本在天球坐標(biāo)系中的分布。通過觀察紅團(tuán)簇星在R-Z平面的分布,可以看到銀盤北面的恒星數(shù)量比銀盤南面多。此外,還發(fā)現(xiàn)年齡較小的紅團(tuán)簇星位于低銀緯附近,而年齡較大的紅團(tuán)簇星位于高銀緯處。

圖2展示了樣本恒星的三維速度分布圖(頂部)和Bootstrap方法給出的三維速度誤差的分布圖(底部)。對于徑向速度VR在R=[10,13] kpc范圍內(nèi)的分布趨勢,方位角速度Vφ圖中的非對稱漂移以及垂直速度VZ在R=[12,14] kpc范圍內(nèi)的分布趨勢都與Wang等[7,9]和Katz等[18]觀察到的結(jié)果一致。

2 結(jié) 果

2.1 屋脊結(jié)構(gòu)在R-Vφ,Rg-Vφ平面中的分布

圖3展示了紅團(tuán)簇星樣本的數(shù)密度分布。在繪制的密度圖中,可以清楚地看到3條屋脊,分別對應(yīng)恒定角動量曲線的值是2 030、2 200和2 400 km·s-1·kpc。圖3右圖則是紅團(tuán)簇星樣本在Rg-Vφ平面的密度分布,不同顏色的實線代表了3條屋脊被轉(zhuǎn)換后對應(yīng)的引導(dǎo)半徑中位值位置(Rg=8.56、9.28、10.12 kpc)。從圖3中可以看到,R-Vφ平面的屋脊結(jié)構(gòu)在Rg-Vφ平面中被轉(zhuǎn)換為了幾乎垂直的屋脊。

圖4展示了樣本恒星的徑向速度分布。圖4左圖中的屋脊結(jié)構(gòu)分布與Yang等[36]利用紅巨星支發(fā)現(xiàn)的屋脊分布一致。不同的是,在紅團(tuán)簇星樣本繪制的徑向速度分布圖中還存在1條疑似屋脊的信號,其遵守的恒定角動量為2 480 km·s-1·kpc。圖4右圖展示了第三維度為徑向速度時,屋脊結(jié)構(gòu)在Rg-Vφ平面中的分布。與圖3類似,發(fā)現(xiàn)的屋脊在Rg-Vφ平面中都轉(zhuǎn)換為了幾乎垂直的過密度區(qū)域,對應(yīng)的引導(dǎo)半徑中位值分別為9.11、9.79、10.46和11.18 kpc。由于對引導(dǎo)半徑范圍的選擇,轉(zhuǎn)換后角動量為1 570和1 750 km·s-1·kpc的屋脊并沒有展示在圖4右圖中。

圖5與圖4類似,展示了紅團(tuán)簇星的垂直速度分布。在垂直速度圖中同樣也發(fā)現(xiàn)了明顯的屋脊信號,其角動量分別為1 480、1 620和3 000 km·s-1·kpc。角動量為3 000 km·s-1·kpc的屋脊與Yang等[36]在垂直速度圖中發(fā)現(xiàn)的屋脊(LZ=2 920 km·s-1·kpc)相似。在圖5右圖中,角動量為3 000 km·s-1·kpc的屋脊被限制在轉(zhuǎn)換后引導(dǎo)半徑的最大值和最小值之間,對應(yīng)的引導(dǎo)半徑中位值為12.66 kpc。因此,R-Vφ平面中具有斜線特征的屋脊會在Rg-Vφ平面被轉(zhuǎn)換為幾乎垂直的屋脊。

2.2 屋脊結(jié)構(gòu)在LZ-JR,LZ-φ平面中的分布

Recio等[20]利用五百六十多萬顆恒星的數(shù)據(jù)樣本在LZ-JR平面中發(fā)現(xiàn)了7條屋脊。圖6展示了紅團(tuán)簇星在LZ-JR平面的分布圖,顏色分別由恒星數(shù)量、垂直距離、金屬豐度和元素豐度著色。在恒星數(shù)量圖中,存在3條疑似屋脊的信號,分別對應(yīng)于圖中的虛線位置。但在垂直距離、金屬豐度和元素豐度分布中沒能探測到疑似屋脊的信號。在后續(xù)工作中,將會利用數(shù)量更大、恒星種類更多的樣本再次探索屋脊在LZ-JR平面的分布。圖7展示了紅團(tuán)簇星樣本的徑向速度和垂直速度在LZ-φ平面的分布。R-Vφ平面中觀察到的屋脊結(jié)構(gòu),在圖7中全部轉(zhuǎn)化為了水平方向的屋脊。

2.3 引導(dǎo)空間中一些關(guān)于銀河系的印記

在前面兩節(jié)中,通過利用計算所得的引導(dǎo)半徑和角動量,分別在Rg-Vφ、LZ-JR和LZ-φ平面探索了屋脊結(jié)構(gòu)的分布,并發(fā)現(xiàn)了一些有趣的特征。圖8左側(cè)展示了紅團(tuán)簇星樣本在引導(dǎo)空間中的徑向速度分布,右側(cè)是垂直速度分布。在徑向速度分布中,總共可以看到4個明顯的結(jié)構(gòu),其位置和Khoperskov等[23]利用Gaia EDR3和DR2發(fā)現(xiàn)的幾個明顯結(jié)構(gòu)一致。通過對比Khoperskov等[37]的研究結(jié)果,角動量空間中的4個過密度區(qū)域可能是由旋臂和銀河系中心棒共振以及矮星系擾動的恒星構(gòu)成的,但需要在后續(xù)工作中運用模擬的方式來證實這個結(jié)論。此外,在垂直速度分布中還揭示了1個關(guān)于銀河系的印記。

3 結(jié) 論

本文使用118 535顆紅團(tuán)簇星,研究了作用量空間中的屋脊結(jié)構(gòu)分布,并探索了銀河系在引導(dǎo)空間中的一些信號。在密度分布、徑向速度分布和垂直速度分布中都探測到數(shù)量不等的屋脊結(jié)構(gòu)。R-Vφ平面的屋脊在轉(zhuǎn)換到Rg-Vφ平面時,全部轉(zhuǎn)換為了幾乎垂直的屋脊。此外,當(dāng)在LZ-φ平面描繪屋脊結(jié)構(gòu)時,徑向速度和垂直速度圖中的屋脊全部轉(zhuǎn)換為了水平的屋脊。結(jié)合角動量和作用量信息,在LZ-JR平面發(fā)現(xiàn)了3條疑似屋脊結(jié)構(gòu)的信號。最后,在引導(dǎo)空間中揭示了5個關(guān)于銀河系的印記,分別是徑向速度分布中的4個過密度區(qū)域和垂直速度分布中的1個過密度區(qū)域。這些特征信號可能與旋臂和銀河系中心棒的共振等有關(guān)。

現(xiàn)階段的工作重心集中在對于屋脊結(jié)構(gòu)的發(fā)現(xiàn)層面。在后續(xù)的工作中,將利用更加精準(zhǔn)的數(shù)據(jù),結(jié)合模擬的方法探索銀河系屋脊結(jié)構(gòu)的起源,并探索其他星系中是否也存在屋脊這種子結(jié)構(gòu)。