COSMOS场中高红移大质量星系的形态分类研究∗

方官文马仲阳 孔 旭

(1大理大学天文与科技史研究所大理671003)

(2中国科学院星系与宇宙学重点实验室合肥230026)

(3中国科学技术大学天文学系合肥230026)

COSMOS场中高红移大质量星系的形态分类研究∗

方官文1†马仲阳2,3孔 旭2,3

(1大理大学天文与科技史研究所大理671003)

(2中国科学院星系与宇宙学重点实验室合肥230026)

(3中国科学技术大学天文学系合肥230026)

基于COSMOS(Cosmic Evolution Survey)天区的多波段测光数据和HST (Hubble Space Telescope)近红外高分辨率观测图像,利用质量限(恒星质量M∗≥1010.5M⊙)选取了362个红移分布在1≤z≤3的星系样本,并对这些大质量星系的形态特征进行了分类研究.来自UVJ(U−V和V−J)双色图分类系统、目视分类系统、非模型化分类系统(基尼系数G和矩指数M20)和模型化分类系统(S´ersic index,n)的分类结果彼此相一致.相比较于恒星形成星系(SFGs),通过UVJ双色图定义的宁静星系(QGs)表现出致密的椭圆结构,而且G和n值偏大,但M20和星系有效半径(re)偏小.不同星系分类系统(双色图分类系统、非模型化分类系统和模型化分类系统)定义的SFGs和QGs样本,都明显存在星系的大小随红移的演化关系,这种演化趋势QGs比SFGs更剧烈,而且不依赖于星系分类方法的选择.

星系:演化,星系:基本参数,星系:结构,星系:高红移

1 引言

高红移大质量(恒星质量M∗≥1010M⊙)星系的观测和研究是当代天文学的一个热点，它对于我们了解星系的形成和演化，约束不同的星系形成模型等都至关重要.通过研究星系的结构随红移的演化,我们可以了解宇宙早期的原初物质是怎样通过引力作用逐渐塌缩形成原初星系,以及高红移星系是通过怎样的过程演变为现在能够观测到的各类星系,从而了解和约束星系的形成和演化过程.表征星系结构的最直接的观测特征为星系的形态，即星系中发光物质(组成星系的恒星和尘埃)在观测者视线方向的投影.因此,研究不同红移处星系的形态特征,并与模型结果相比较,可以了解星系、乃至宇宙的形成和演化[1−2].星系形态是反映星系物理性质的一个重要特征量,它与星系的其他物理性质有一定的关系,如早型星系年老、颜色红、消光弱、多处于高密度区、恒星形成率低,而晚型星系有着相反的特性[3].为了研究不同类型星系的物理特征,我们首先要做的是有效区分星系的不同形态,这样一个合理的分类方法就显得非常重要.

不同的星系分类方法依据的分类标准不尽相同,这里主要介绍一下目视分类系统、模型化分类系统、非模型化分类系统和双色图分类系统.(1)利用目视的分类方法,很多学者提出了不同的星系分类方案,其中最著名的为哈勃在1926年提出的“音叉图”星系形态分类系统[4].该分类系统将星系形态分类成:椭圆星系、透镜星系、旋涡星系和不规则星系.(2)模型化分类系统是利用星系的面亮度轮廓对星系形态进行分类,不同形态星系的面亮度(即星系每单位面积发射的视光强度)分布是不同的[5].椭圆星系的面亮度轮廓可以较好地利用r1/4(r是星系面亮度分布半径)律进行拟合,而盘主导的星系面亮度满足指数律分布.(3)非模型化分类系统是基于表征星系形态的结构参数来对星系形态进行分类,这些参数包括基尼系数(G)、矩指数(M20)、聚集度指数(C)以及非对称指数(A)等[6].早型椭圆状星系有较大的G和C,但M20和A较小,而对于晚型旋涡星系(或不规则星系)则有着相反的性质.(4)双色图分类系统是利用不同类型星系能谱分布的差异,将星系分类成恒星形成星系(SFGs)和宁静星系(QGs).例如在静止坐标系中的UVJ(U−V和V−J)双色图上,Muzzin等[7]定义QGs(对于红移分布在11.3,V−J<1.5且U−V>0.88(V−J)+0.59,而SFGs则分布在双色图中相反的区域.

以上我们介绍了星系的目视分类系统、模型化分类系统、非模型化分类系统和双色图分类系统,更多关于星系分类方法的描述和分析可参考文献Sandage[8]和汪敏等[9].相对于目视分类方法只适合小样本和具有很强的主观性,模型化分类方法需要设定星系光度分布函数形式.在星系的具体面亮度分布未知的情况下,通常利用S´ersic提出的面亮度轮廓函数(r1/n,n称为S´ersic指数)来对星系面亮度分布进行拟合.对于在光学波段近邻早型星系的面亮度分布一般满足n>2,而晚型星系对应n<2[10−11].非模型化分类方法的优势是基于计算机自动测量,不需要考虑星系光度分布形式,所以对大样本星系形态的研究很适用,而且这些形态结构参数有着丰富的物理意义[2].例如:M20和A与星系相互作用相关,G和C反映星系辐射流量分布性质.至于双色图分类系统定义的SFGs和QGs,污染源主要来自颜色分界线附近的星系.

来自星系分类系统定义的恒星形成星系(SFGs)和宁静星系(QGs),它们的形态结构差异反映出星系不同的形成和演化过程[2].基于高分辨率的光学和近红外波段观测图像,最近很多研究工作发现在相同恒星质量条件下,在红移13时,大质量星系的大部分恒星已经形成,并在随后的演化过程中是被动演化,没有并合事件发生,星系的大小随红移没有明显的演化关系.

本文我们基于COSMOS(Cosmic Evolution Survey)/UltraVISTA(Ultra-deep Visible and Infrared Survey Telescope for Astronomy)场中多波段测光数据和HST WFC3 (Hubble Space Telescope Wide Field Camera 3)近红外高分辨率的观测图像,利用目视分类系统、模型化分类系统、非模型化分类系统和双色图分类系统对不同类型的高红移大质量星系进行分类研究,并比较不同分类方法之间的差异.另外,我们也会讨论不同星系分类系统定义的SFGs和QGs样本,是否会导致不同类型星系的大小随红移的演化关系发生明显的变化.

2 数据和样本选择

COSMOS/UltraVISTA巡天场中既有窄带光谱巡天数据[29],又有多波段宽带观测图像(波长范围从X射线到射电波段)[30],特别是来自CANDELS[31−32](Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey)巡天的HST WFC3近红外(F125W和F160W)高分辨率(0.06′′pixel−1)的观测图像.我们的研究工作是基于Muzzin等[33]提供的星表,该星表不仅包括多波段测光数据,而且提供了星系各种物理参数测量.例如,测光红移,恒星质量,静止坐标系中的颜色(U−V和V−J)等.关于COSMOS巡天以及多波段观测数据的描述也可参考Fang等[34−35]文献.

为了使质量限选取样本的完备性高于95%,我们利用M∗≥1010.5M⊙的选择判据从CANDELS-COSMOS场(天区有效面积∼210 arcmin2)中选取了362个红移分布在1≤z≤3星系样本,这些源都有HST WFC3近红外(F125W和F160W)观测数据.考虑到星系形态结构性质对观测波长的依赖,对于红移处在1≤z<1.76范围内的星系,我们利用F125W波段观测图像研究它们的形态结构特征.而对于红移分布在1.76

3 星系的形态分类

基于Muzzin等[7]定义的UVJ(U−V和V−J)双色图分类判据,我们将样本中的星系分类成恒星形成星系(SFGs)和宁静星系(QGs).如图1所示,其中152个QGs分布在UVJ双色图上的左上角,210个SFGs则分布在双色图中相反的区域,它们占整个样本的58%.图中的实线对应的颜色判据为:U−V=1.3,V−J=1.5和U−V=0.88(V−J)+0.59.为了在UVJ双色图上研究星系的S´ersic指数分布情况,我们将样本中的星系依据红移大小分成4个红移区间,每个红移区间大小是0.5,如图2所示.从该图中我们发现被UVJ划分为QGs的源S´ersic指数偏大,主要分布在n>2区域,而SFGs的S´ersic指数明显偏小,大部分星系的n值小于2.对于1.5

图1 星系在UVJ双色图上的分布(1≤z≤3).图中的实线是来自Muzzin等[7]提供的颜色判据(U−V=1.3, V−J=1.5和U−V=0.88(V−J)+0.59).Fig.1 The distribution of galaxies in the(U−V)restvs.(V−J)restdiagram(1≤z≤3).The solid lines correspond to the color criteria from Muzzin et al.[7](U−V=1.3,V−J=1.5,and U−V=0.88(V−J)+0.59).

利用哈勃空间望远镜WFC3近红外(F125W和F160W)高分辨率观测图像,我们对样本中的星系进行了目视形态分类.分类的依据类似于哈勃“音叉图”星系形态分类系统,其形态类型包括椭圆星系、盘星系和不规则星系.(1)椭圆星系:从图像观察星系表现出椭球或球状特征,而且没有延展子结构,星系辐射的光度聚集在星系中心区域.(2)盘星系:该类型星系从外表上看有单一的亮核并有延展的盘状结构.(3)不规则星系:结构上缺少对称性或没有明显旋臂特征,而且有一些团块状子结构,有些源明显带有相互作用特征(如潮汐臂、气体桥、双不对称体等).图3提供了上述星系形态目视分类的典型例子,每个源对应的G和M20也被标注在每个小图像上.在图4中,我们进一步比较了星系目视分类和双色图分类的结果.在星系被UVJ判据划分为SFGs的样本中,椭圆星系、盘星系和不规则星系所占比例依次为:27%,31%,42%,而在QGs样本中,相应的比例是:80%,14%,6%.这表明恒星形成星系的形态主要是盘状或不规则结构(73%),但对于宁静星系来说,主导的结构特征是椭圆状(80%).另外,从图4中我们也可以看到,盘星系主要分布在颜色分界线附近,但在更高红移区间内(1.76

图2 不同红移区间内星系S´ersic指数在UVJ双色图上的分布.不同大小的S´ersic指数用不同颜色表示.每个红移区间内星系的数目,以及QGs和SFGs所占的比例也标记在该图中.Fig.2 The distribution of S´ersic indexes of galaxies for di ff erent redshift bins in the(U−V)restvs. (V−J)restdiagram.Di ff erent values of S´ersic indexes are plotted in di ff erent colors.The number of galaxies,and the fraction of QGs and SFGs in each redshift bin are respectively shown in this figure.

图3 星系形态目视分类的典型例子.上图:椭圆星系,中图:盘星系,下图:不规则星系.每个源对应的G和M20也被标注在每个小图像上.每个小图像的尺寸是4′′×4′′.Fig.3 Examples of three general morphological types from visual inspection.Top:spheroid,middle:disk, bottom:irregular.Values of G and M20are shown in each stamp image.The size of each image is 4′′×4′′.

图4 星系目视分类结果在UVJ双色图上的分布.红色菱形、绿色方块和蓝色三角形分别表示椭圆星系、盘星系和不规则星系.Fig.4 The distribution of visual morphological types of galaxies for di ff erent redshift bins in the (U−V)restvs.(V−J)restdiagram.The red diamond,green square,and blue triangle represent spheroid, disk,and irregular galaxies by visual inspection,respectively.

为了定量描述样本中星系的形态性质,我们利用HST WFC3近红外观测数据测量了样本中星系的形态参数(G和M20).图5显示了所有源的形态参数分布,其中红色和蓝色的点代表UVJ分类结果,而形状不同的数据点是由目视分类给出.图中黑色实线表示的是非模型化分类系统定义的形态划分标准,星系的形态参数(G和M20)值满足G>0.125M20+0.8判据则被归类到早型星系,反之则属于晚型星系.从该图中我们发现,目视分类中的椭圆星系、盘星系和不规则星系,对应的G值逐渐减小,而M20则增大,盘星系的形态参数分布位于椭圆星系和不规则星系之间.这说明椭圆星系光的聚集度最大,不对称性最小,而不规则星系却恰恰相反.另外,相比较于双色图定义的SFGs,样本中QGs的G偏大而M20偏小.以上结论表明:来自双色图分类系统、目视分类系统和非模型化分类系统的星系分类结果彼此相一致.

4 星系的结构演化

当前关于星系形成和演化主要有等级成团模型[40]和纯光度演化模型[28]两种模型.等级成团模型认为星系起源于暗物质晕的塌缩和并合,重子物质在暗物质晕中冷却并塌缩形成星系盘.盘星系通过并合破坏盘状结构形成椭圆星系.在这种模型中先有小质量的星系形成,而后通过不断地并合形成近邻的大质量星系.检验这两种模型一个最直接的有效方法是研究星系的大小re随红移演化关系.基于来自CANDELS-COSMOS场中所选的大质量星系样本,我们利用UVJ定义的颜色判据将星系分类成SFGs和QGs.在每一类型样本当中,又根据红移大小将其划分成不同的区间.如图6所示,这里给出了两个红移区间内星系有效半径与恒星质量的关系,黑色实线是来自Shen等[11]对近邻星系的研究结果,其中LTGs和ETGs分别表示晚型星系和早型星系.从该图我们可以得到这样一些信息:(1)高红移星系有着类似于近邻星系的re∝Mα∗关系,指数α的差异取决于星系的形态类型;(2)相比较于同类型的近邻星系,高红移星系有效半径整体偏小,而且质量比较大的星系re偏大,这点对于QGs星系更加明显;(3)在高红移处不论SFGs还是QGs,都存在极其致密的源(re<1 kpc),这些星系的形成很可能来自更高红移处富气体星系间的并合.

图5 样本中的星系在G−M20图上的分布.红色(QGs)和蓝色(SFGs)的点代表UVJ分类结果.形状不同的数据点是由目视分类给出,菱形、方块和三角形分别代表椭圆星系、盘星系和不规则星系.图中黑色实线表示的是非模型化分类系统定义的形态划分标准:G=0.125M20+0.8.Fig.5Distribution of galaxies in the M20vs.G plane.The SFGs are plotted in blue color,while QGs are in red.Spheroid,disk,and irregular galaxies which are classified by visual inspection,are plotted in diamond,square,and triangle symbols,respectively.The solid lines represent the defined criterion of nonparametric morphology(G=0.125M20+0.8).

图7进一步显示了星系有效半径随红移的演化关系,红色点和蓝色点分别代表UVJ定义的QGs和SFGs,其中带有误差棒的方形点表示不同红移区间内(红移区间大小是0.5)re的中值,红色线(QGs)和蓝色线(SFGs)是对这些中值的最佳拟合.对于QGs,由于在红移2.5≤z<3范围内源的数目太少,因此拟合时并没有采用这个红移区间内的数据.与此同时,我们在图7中也给出了其他分类系统对应的不同类型星系大小随红移的演化关系,正如上一节所描述的那样,模型化分类系统采用的是n=2判据,而非模型化分类系统则利用G=0.125M20+0.8标准,为了不影响图的显示效果,这里我们就没有画出这两种分类方法产生的数据点,只给出最佳拟合曲线.从图7我们可以得出这些结论:(1)高红移大质量星系(QGs和SFGs)的大小随红移有明显的演化关系,这种演化趋势QGs比SFGs更剧烈,这表明QGs大小的增长可能是由于星系间的微并合导致的,这种并合在QGs演化过程中可能发生多次;(2)不同星系分类方法得到的QGs和SFGs样本都存在(1)的结论,这说明星系大小随红移演化的关系不依赖于星系分类方法的选择;(3)我们证实了以前的研究结果,在相同恒星质量条件下,高红移星系(1

图6 不同红移区间内星系有效半径与恒星质量的关系.黑色实线是来自Shen等[11]对近邻星系的研究结果,其中LTGs和ETGs分别表示晚型星系和早型星系.Fig.6 The e ff ective radius vs.stellar mass diagram for our sample.SFGs are plotted in(a)and(b) panels,while QGs are plotted in(c)and(d)panels.The solid lines represent the fitting curves for late type galaxies(LTGs)and early type galaxies(ETGs)at z∼0 from Shen et al.[11].

5 总结

基于CANDELS-COSMOS场中多波段测光数据,利用质量限(M∗≥1010.5M⊙),我们构建了一个高红移(1≤z≤3)大质量星系样本,样本包括362个星系.这些源都有HST WFC3近红外(F125W和F160W)高分辨率的观测数据,我们的形态结构分析是依据星系红移大小的差异而选择不同的观测图像进行研究,对于红移处在1≤z<1.76范围内的星系,我们利用F125W波段观测数据研究它们的形态结构特征.而对于红移分布在1.76

图7 大质量(M∗>1010.5M⊙)星系的有效半径(re)随红移的演化关系.带有误差棒的方形点表示不同红移区间内re的中值,红色线(QGs)和蓝色线(SFGs)是对这些中值的最佳拟合.绿线和青色线分别表示模型化分类系统和非模型化分类系统对应的不同类型星系大小随红移的演化关系.Fig.7 Evolution of the e ff ective radius(re)with redshift for the massive galaxies with estimated stellar masses M∗>1010.5M⊙.The median values of e ff ective radius of QGs and SFGs in each redshift bin are represented by red and blue squares,respectively.The red and blue lines are the best-fit curves respectively for QGs and SFGs in our sample.The green and cyan lines represent the size evolution with redshift for di ff erent divisions of QG and SFG samples,respectively.

在本工作中,我们采用了目视分类系统、模型化分类系统、非模型化分类系统和双色图分类系统对不同类型的高红移大质量星系进行了分类研究,并比较了不同分类方法之间的差异,结果表明不同星系分类方法给出的结果彼此相一致.相比较于恒星形成星系(SFGs),通过UVJ双色图定义的宁静星系(QGs)表现出致密的椭圆结构,而且G和n值偏大,但M20和星系有效半径(re)偏小.与此同时,我们也发现不同星系分类方法得到的QGs和SFGs样本,对于同一类型星系都有相似的星系大小随红移的演化关系,而且这种演化趋势QGs比SFGs更剧烈,这说明星系大小随红移演化的关系不依赖于星系分类方法的选择.

[1]Ravindranath S,Giavalisco M,Ferguson H C,et al.ApJ,2006,652:963

[2]Conselice C J.ARA&A,2014,52:291

[3]Ball N M,Loveday J,Brunner R J.MNRAS,2008,383:907

[4]Hubble E.ApJ,1926,64:321

[5]Peng C Y,Ho L C,Impey C D,et al.AJ,2002,124:266

[6]Lotz J M,Primack J,Madau P.AJ,2004,128:163

[7]Muzzin A,Marchesini D,Stefanon M,et al.ApJ,2013,777:18

[8]Sandage A.ARA&A,2005,43:581

[9]汪敏,孔旭.天文学进展,2007,25:215

[10]Blanton M R,Hogg D W,Bahcall N A,et al.ApJ,2003,592:819

[11]Shen S,Mo H J,White S D M,et al.MNRAS,2003,343:978

[12]van der Wel A,Franx M,van Dokkum P G,et al.ApJ,2014,788:28

[13]Trujillo I,Conselice C J,Bundy K,et al.MNRAS,2007,382:109

[14]Cimatti A,Cassata P,Pozzetti L,et al.A&A,2008,482:21

[15]van Dokkum P G,Kriek M,Franx M.ApJL,2008,677:L5

[16]Akiyama M,Minowa Y,Kobayashi N,et al.ApJS,2008,175:1

[17]Franx M,van Dokkum P G,F¨orster Schreiber N M,et al.ApJL,2008,689:L101

[18]Szomoru D,Franx M,van Dokkum P G,et al.ApJL,2010,714:L244

[19]Szomoru D,Franx M,van Dokkum P G,et al.ApJ,2012,749:121

[20]van der Wel A,Rix H W,Wuyts S,et al.ApJ,2011,730:38

[21]Barro G,Faber S M,P´erez-Gonz´alez P G,et al.ApJ,2013,765:104

[22]Patel S G,van Dokkum P G,Franx M,et al.ApJ,2013,766:15

[23]Morishita T,Ichikawa T,Kajisawa M.ApJ,2014,785:18

[24]Naab T,Johansson P H,Ostriker J P.ApJL,2009,699:L178

[25]Hopkins P F,Bundy K,Hernquist L,et al.MNRAS,2010,401:1099

[26]Carollo C M,Bschorr T J,Renzini A,et al.ApJ,2013,773:112

[27]Poggianti B M,Moretti A,Calvi R,et al.ApJ,2013,777:125

[28]Larson R.MNRAS,1975,173:671

[29]Lilly S J,Le F`evre O,Renzini A,et al.ApJS,2007,172:70

[30]Scoville N,Aussel H,Brusa M,et al.ApJS,2007,172:1

[31]Grogin N A,Kocevski D D,Faber S M,et al.ApJS,2011,197:35

[32]Koekemoer A M,Faber S M,Ferguson H C,et al.ApJS,2011,197:36

[33]Muzzin A,Marchesini D,Stefanon M,et al.ApJS,2013,206:8

[34]方官文,马仲阳,陈洋,等.天文学报,2015,56:7

[35]Fang G W,Ma Z Y,Chen Y,et al.ChA&A,2015,39:307

[36]方官文,马仲阳,陈洋,等.天文学报,2014,55:458

[37]Fang G W,Ma Z Y,Chen Y,et al.ChA&A,2015,39:156

[38]方官文.天文学报,2014,55:444

[39]van der Wel A,Bell E F,H¨aussler B,et al.ApJS,2012,203:24

[40]White S D M,Frenk C S.ApJ,1991,379:52

Morphological Classification of High-redshift Massive Galaxies in the COSMOS/UltraVISTA Field

FANG Guan-wen1MA Zhong-yang2,3KONG Xu2,3
(1 Institute for Astronomy and History of Science and Technology,Dali University,Dali 671003)
(2 Key Laboratory for Research in Galaxies and Cosmology,Chinese Academy of Sciences,Hefei 230026)
(3 Department of Astronomy,University of Science and Technology of China,Hefei 230026)

Utilizing the multi-band photometry catalog of the COSMOS(Cosmic Evolution Survey)/UltraVISTA(Ultra-deep Visible and Infrared Survey Telescope for Astronomy)field and the high-resolution HST WFC3(Hubble Space Telescope Wide Field Camera 3)near-infrared imaging from the CANDELS(Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey)field,we present a quantitative study of the morphological classification of galaxy for a large mass-selected sample.Our sample includes 362 galaxies within photometric redshift 1≤z≤3 and stellar massM∗≥1010.5M⊙. The results from the rest-frame(U−V)vs.(V−J)(UVJ)colors classification,visualinspection,nonparametric morphology analysis,and structural parameters study are in good agreement with each other.Quiescent galaxies(QGs)classified by UVJ colors generally have larger S´ersic index(n)and Gini coefficient(G),smaller size(re)and moment(M20),and they are visually compact.While star-forming galaxies(SFGs) are reversed.In the meantime,we explore the size evolution with redshift for various divisions of QG and SFG samples,and confirm that both of size will enlarge with time, but QGs are rapider than SFGs.Moreover,we find that the choice of division between QGs and SFGs(i.e.colour,shape,morphology)is not particularly critical.

galaxies:evolution,galaxies:fundamental parameters,galaxies:structure, galaxies:high-redshift

P157;

A

10.15940/j.cnki.0001-5245.2015.05.002

2015-03-27收到原稿,2015-04-21收到修改稿

∗国家自然科学基金项目(11303002,11225315,11320101002)、中国科学院战略性先导科技专项(XDB09000000)和云南省应用基础研究计划项目(2014FB155)资助

†wen@mail.ustc.edu.cn