星系團(tuán)數(shù)據(jù)分析方法及Python語(yǔ)言實(shí)現(xiàn)

茍?zhí)熳?/p>

（貴州大學(xué) 理學(xué)院，貴州 貴陽(yáng) 550025）

目前，天文學(xué)研究領(lǐng)域內(nèi)已經(jīng)展開了許多星系巡天的觀測(cè)研 究， 如Sloan數(shù)字巡天計(jì)劃 （SDSS）、兩微米全天巡天（2MASS）、Arecibo 中性氫巡天（ALFA）等。不同的巡天項(xiàng)目，有不同的數(shù)據(jù)庫(kù)，這給多波段整合研究帶來(lái)了困難。

采用多個(gè)波段觀測(cè)數(shù)據(jù)的對(duì)比研究方法是研究星系特征的一個(gè)有效手段。在星系團(tuán)的特征研究中，能夠在星系團(tuán)內(nèi)證認(rèn)出星系成員是非常關(guān)鍵的一步。 最廣泛的一種證認(rèn)方法是Huchra 和Geller 的找朋友算法（Friends of Friends）。

本文主要描述了統(tǒng)一不同數(shù)據(jù)庫(kù)坐標(biāo)格式的方法、尋找中性氫觀測(cè)的光學(xué)對(duì)應(yīng)體方法和證認(rèn)星系團(tuán)成員的方法。 第一章介紹Sloan 數(shù)字巡天計(jì)劃（SDSS）和Arecibo 中性氫巡天（ALFA）項(xiàng)目，第二章介紹星系團(tuán)數(shù)據(jù)分析方法和Python 語(yǔ)言實(shí)現(xiàn)，第三章進(jìn)行簡(jiǎn)單的總結(jié)。

1 光學(xué)巡天與中性氫巡天

由于不同的巡天項(xiàng)目有著不同的觀測(cè)原理，從而形成了對(duì)同一位置的觀測(cè)目標(biāo)不一致現(xiàn)象，并且觀測(cè)結(jié)果表示形式不一，因此在采用多波段觀測(cè)研究星系特征的方法 （如采用光學(xué)巡天和中性氫氣巡天多波段研究星系特征）時(shí)，就需要找到中性氫巡天的光學(xué)對(duì)應(yīng)體并且將單位表示統(tǒng)一化。

1.1 Sloan 數(shù)字巡天計(jì)劃（SDSS）

Sloan 數(shù)字巡天第九批釋放的數(shù)據(jù) （Data Release 9,DR9）[3]實(shí)現(xiàn)了14 555 平方度的天空覆蓋，含有恒星光譜668054 個(gè)、星系光譜1457 002 個(gè)。 DR9 數(shù)據(jù)庫(kù)內(nèi)包含有許多的表格，如PhotoObjAll、Photoz、Region、StellarMassPassivePor 和 TwoMass等。 在DR9 數(shù)據(jù)庫(kù)內(nèi)查詢表格內(nèi)的信息需要使用SQL 查詢語(yǔ)言。 查詢時(shí)需要整理出需要的觀測(cè)數(shù)據(jù)， 然后在casjobs 內(nèi)用SQL 語(yǔ)言實(shí)現(xiàn)表格的交叉查詢。 星系的位置坐標(biāo)，查詢結(jié)果顯示的是以角度為單位，即赤經(jīng)和赤緯單位都是角度。

1.2 Arecibo 中性氫巡天（ALFA）

ALFALFA 觀測(cè)數(shù)據(jù)已經(jīng)釋放了40%的數(shù)據(jù)庫(kù)（a.40）[4]內(nèi)包含的觀測(cè)參數(shù)包括有星系的位置（赤經(jīng)和赤緯）、中性氫譜線流量密度、觀測(cè)速度達(dá)到峰值的50%的觀測(cè)值W50 和中性氫的觀測(cè)類型等。 其中，星系的位置坐標(biāo)表示形式為，赤經(jīng)(hhmmss.s)+赤緯(ddmmss)的形式。

2 星系團(tuán)數(shù)據(jù)分析方法及Python 語(yǔ)言實(shí)現(xiàn)

采用整合多波段觀測(cè)數(shù)據(jù)的方法來(lái)對(duì)比研究星系團(tuán)的特征， 需要統(tǒng)一多波段的觀測(cè)數(shù)據(jù)格式和尋找中性氫觀測(cè)的光學(xué)對(duì)應(yīng)體，之后對(duì)星系團(tuán)成員進(jìn)行證認(rèn)。本文主要運(yùn)用Python 語(yǔ)言來(lái)實(shí)現(xiàn)統(tǒng)一不同數(shù)據(jù)庫(kù)的數(shù)據(jù)格式、 尋找光學(xué)對(duì)應(yīng)體和證認(rèn)星系團(tuán)成員的算法（F0F）。 Python[2]是一種非常流行的面向?qū)ο蟮母呒?jí)解釋型語(yǔ)言，其語(yǔ)法非常簡(jiǎn)潔并且便于人理解，擁有便捷的文本處理方法和很多強(qiáng)大的數(shù)學(xué)庫(kù)、繪圖庫(kù)，天然支持高精度計(jì)算，所以越來(lái)越多的應(yīng)用在了科學(xué)計(jì)算中。

2.1 尋找中性氫觀測(cè)的光學(xué)對(duì)應(yīng)體

光學(xué)巡天數(shù)據(jù)庫(kù)（DR9）內(nèi)星系數(shù)據(jù)查詢結(jié)果中，星系的位置坐標(biāo)（赤經(jīng)和赤緯）是以角度單位。中性氫巡天數(shù)據(jù)庫(kù)（a.40）中星系位置坐標(biāo)表示形式是， 赤經(jīng)+赤緯， 單位表示為， 時(shí)分秒(hhmmss.s)+度分秒(ddmmss)的形式。將星系位置坐標(biāo)格式統(tǒng)一為以角度為單位的方法，我們是利用python 語(yǔ)言來(lái)實(shí)現(xiàn)。

ra 和dec 分別表示星系轉(zhuǎn)換以后的赤經(jīng)和赤緯。 將a.40 的位置坐標(biāo)轉(zhuǎn)換為角度制的方法如下：

ra=（hh+mm÷60+ss÷3 600）×15

dec=（dd×1）+（mm÷60）+（ss÷3 600）

統(tǒng)一位置坐標(biāo)單位以后，將中性氫巡天數(shù)據(jù)庫(kù)（a.40）內(nèi)的星系位置坐標(biāo)與光學(xué)巡天數(shù)據(jù)庫(kù) （DR9）內(nèi)的星系位置坐標(biāo)相減，滿足小于1 個(gè)角分，即小于約0.016 度的兩個(gè)星系視為可能為同一個(gè)星系。

這個(gè)星系在中性氫巡天數(shù)據(jù)庫(kù) （a.40）和光學(xué)巡天數(shù)據(jù)庫(kù)（DR9）內(nèi)，紅移的比較滿足以下關(guān)系的確定為同一個(gè)星系。

|ZHI-Zoc|

式中的ZHI、Zoc和zErr 分別表示， 中性氫巡天數(shù)據(jù)庫(kù)（a.40）中的星系紅移、光學(xué)巡天數(shù)據(jù)庫(kù)（DR9）中的光學(xué)紅移和光學(xué)巡天數(shù)據(jù)庫(kù)（DR9）中的光學(xué)紅移誤差。

2.2 找朋友算法（Friends Of Friends）

圖1 FOF 算法流程圖

本文引用的是Huchra 和Geller 的找朋友算法（FOF）[1]來(lái)實(shí)現(xiàn)證認(rèn)星系團(tuán)成員，這是最常用的一種驗(yàn)證星系成團(tuán)性的算法。具體的算法流程圖如下：

找朋友算法主要是計(jì)算星系的3 個(gè)物理參量 （星系的赤經(jīng)、赤緯、星系的測(cè)光紅移）來(lái)尋找星系的成員。

其中，θij表示星系i 和星系j 之間的角距離，H0表示哈勃常量， 這對(duì)星系的平均速度為，vi和vj分別表示星系i 和星系j 的速度。

星系i 和星系j 必須滿足以下三個(gè)條件才能視為屬于一個(gè)星系群或星系團(tuán)內(nèi)：

參數(shù)DL和VL主要是根據(jù)使用的星系數(shù)據(jù)庫(kù)來(lái)具體取值。一般將研究中將群的最小成員數(shù)Nmin設(shè)置為3 個(gè)。

找朋友算法主要是計(jì)算星系的3 個(gè)物理參量 （星系的赤經(jīng)、赤緯、星系的測(cè)光紅移）來(lái)尋找星系的成員。

其中，θij表示星系i 和星系j 之間的角距離，H0表示哈勃常量， 這對(duì)星系的平均速度為和vj分別表示星系i 和星系j 的速度。

星系i 和星系j 必須滿足以下三個(gè)條件才能視為屬于一個(gè)星系群或星系團(tuán)內(nèi)：

參數(shù)DL和VL主要是根據(jù)使用的星系數(shù)據(jù)庫(kù)來(lái)具體取值。一般將研究中將群的最小成員數(shù)Nmin設(shè)置為3 個(gè)。

用以上方法就可以輸出在一個(gè)天區(qū)內(nèi)的星系團(tuán)成員情況表3FOF 算法的Python 實(shí)現(xiàn)通過(guò)上述過(guò)程，執(zhí)行Python 程序，本文從光學(xué)巡天數(shù)據(jù)(DR9)和中性氫巡天數(shù)據(jù)(a.40)中獲得了星系團(tuán)的特征數(shù)據(jù)。 通過(guò)該數(shù)據(jù)，研究者可以做出一系列的星系團(tuán)特征圖，來(lái)分析星系團(tuán)的結(jié)構(gòu)與環(huán)境特征。 本文通過(guò)這個(gè)特征數(shù)據(jù)的分析，得到了與黃珊等人[5]做出的a.40 星系的特性研究結(jié)果一致。

3 總結(jié)

由于不同的天文觀測(cè)項(xiàng)目有不同的觀測(cè)原理，所以在不同的數(shù)據(jù)庫(kù)內(nèi)查詢的星系數(shù)據(jù)，通常需要經(jīng)過(guò)轉(zhuǎn)換坐標(biāo)單位格式、辨認(rèn)同一星系、搜集同一星系的多波段數(shù)據(jù)、整理數(shù)據(jù)計(jì)算物理參量來(lái)研究星系特征的過(guò)程。 研究星系團(tuán)特征前，一般需要證認(rèn)其星系團(tuán)成員后才能研究其特征。 本文介紹的證認(rèn)方法是找朋友算法（FOF）。 通過(guò)python 語(yǔ)言來(lái)實(shí)現(xiàn)a.40 內(nèi)的光學(xué)對(duì)應(yīng)體的找尋過(guò)程，進(jìn)而整合出星系的光學(xué)數(shù)據(jù)和中性氫的觀測(cè)數(shù)據(jù)，以及實(shí)現(xiàn)了找朋友算法。 數(shù)據(jù)的分析結(jié)果表明了我們的研究方法是正確可靠的，為下一步的科學(xué)研究打下了良好的數(shù)據(jù)基礎(chǔ)。

[1]C.S. Botzler, J. Snigula, R.Bender, et al. Finding structures in photometric redshift galaxy survey: An extended friends-offriends algorithm.Mon.Not.Roy.Astron.Soc.349(2004)425.

[2]Python:http://www.python.org/.

[3]dr9：http://skyserver.sdss3.org/CasJobs/.

[4]a.40：http://egg.astro.cornell.edu/alfalfa/data/index.php.

[5]Shan Huang，Martha P.Haynes，Riccardo Giovanelli，et al.THE ARECIBO LEGACY FAST ALFA SURVEY:THE GALAXY POPULATION DETECTED BY ALFALFA.The Astrophysical Journal,Volume 756,Issue 2,arti cle id.113,29 pp.(2012).