吳淑梅，張福鵬，陸燁，陸由俊

(中國(guó)科學(xué)院國(guó)家天文臺(tái)，北京 100012)

低紅移SDSS類星體輻射效率的估計(jì)

吳淑梅，張福鵬，陸燁，陸由俊

(中國(guó)科學(xué)院國(guó)家天文臺(tái)，北京 100012)

研究了低紅移SDSS類星體(0.025＜z＜0.5)的輻射效率。首先利用類星體中心黑洞的質(zhì)量估計(jì)，通過(guò)薄盤吸積模型估計(jì)了每一個(gè)類星體的吸積率。其次根據(jù)光學(xué)波段的觀測(cè)光度，利用經(jīng)驗(yàn)的類星體全波段模板譜估計(jì)了每個(gè)類星體的熱光度。最后由估計(jì)的吸積率和熱光度得到了每個(gè)類星體的輻射效率。發(fā)現(xiàn)低紅移SDSS類星體的輻射效率與黑洞質(zhì)量強(qiáng)相關(guān)，并滿足ε∝M?0.63。這一結(jié)果與Davis和Laor(2011)對(duì)PG類星體樣本得到的結(jié)果一致。進(jìn)一步討論由于輻射效率估計(jì)中各種假設(shè)可能引入的偏差及其對(duì)ε與M?關(guān)系的可能影響。

吸積盤；輻射效率；黑洞；吸積率；自旋

CN53－1189/P ISSN1672－7673

一般認(rèn)為大質(zhì)量黑洞的物質(zhì)吸積是活動(dòng)星系核和類星體的能量來(lái)源。類星體的光學(xué)紫外輻射近似地可用幾何薄光學(xué)厚的標(biāo)準(zhǔn)薄吸積盤描述[1]。有關(guān)大質(zhì)量黑洞和類星體的觀測(cè)統(tǒng)計(jì)表明其中心黑洞是通過(guò)類星體階段的物質(zhì)吸積獲得絕大部分質(zhì)量[2－3]。大質(zhì)量黑洞的形成演化和類星體的吸積歷史密切相關(guān)。例如，如果吸積過(guò)程是長(zhǎng)期的薄盤吸積，其中心黑洞的自旋會(huì)快速演化至接近于1的一個(gè)平衡態(tài)值；但若吸積過(guò)程為具有隨機(jī)角動(dòng)量的多個(gè)短時(shí)標(biāo)吸積，吸積階段測(cè)量其黑洞自旋很可能接近于0[4－5]。

在標(biāo)準(zhǔn)薄盤模型中，類星體的光譜發(fā)射主要取決于幾個(gè)基本的參數(shù)，即黑洞的質(zhì)量M?、自旋a、質(zhì)量吸積率acc以及吸積盤相對(duì)于觀測(cè)視線的傾角i等。測(cè)量類星體的這些基本參數(shù)例如M?和a及估計(jì)其統(tǒng)計(jì)分布成為當(dāng)今類星體研究的一個(gè)重要目標(biāo)。這些測(cè)量不僅可以幫助理解單個(gè)類星體的基本特征、限制吸積盤模型，也可以幫助理解大質(zhì)量黑洞和類星體的宇宙學(xué)形成和演化歷史。

類星體中心黑洞的質(zhì)量可以通過(guò)反響映射法進(jìn)行測(cè)量[6－11]，但需要積累高精度的測(cè)光和分光觀測(cè)資料。通常人們采用由反響映射法得到的經(jīng)驗(yàn)關(guān)系即維里質(zhì)量估算方法估計(jì)黑洞的質(zhì)量。維里質(zhì)量估算法可以簡(jiǎn)單通過(guò)類星體光學(xué)波段光度以及寬發(fā)射線(如Hβ、CⅣ、MgⅡ等)的線寬估計(jì)黑洞的質(zhì)量，因而便于應(yīng)用。對(duì)于黑洞自旋的測(cè)量目前還沒有可靠的辦法，盡管對(duì)于一些特別的類星體，如MCG-6-30-15[12]，其中心黑洞的自旋可以通過(guò)Fe Kα線或其它方法測(cè)得，但是絕大多數(shù)的類星體中心黑洞的自旋難以直接測(cè)量。

黑洞吸積周圍氣體的過(guò)程中會(huì)把一部分物質(zhì)轉(zhuǎn)化為能量，這種質(zhì)能轉(zhuǎn)化的效率稱為輻射效率(記為ε)。在薄吸積盤理論中，輻射效率是由黑洞的自旋決定的[13]。如果能得到黑洞的輻射效率，那么原則上可以進(jìn)一步限制黑洞的自旋。類星體的平均輻射效率可以通過(guò)對(duì)比類星體的總能量發(fā)射密度與近鄰黑洞的質(zhì)量密度獲得，大約為ε～0.1[2－3，14－18]。但是對(duì)于單個(gè)類星體，由于受觀測(cè)水平的限制，導(dǎo)致其輻射效率很難被精確地估計(jì)。如果已知黑洞的質(zhì)量M?，可以通過(guò)吸積盤模型擬合光學(xué)波段的光度估計(jì)黑洞的質(zhì)量吸積率acc。若同時(shí)可以獨(dú)立地估計(jì)類星體的熱光度Lbol，則可以得到類星體的輻射效率，即ε＝Lbol/(accc2)。Davis＆Laor[19](2011，簡(jiǎn)寫為DL)采用這種方法估計(jì)了80個(gè)PG類星體的輻射效率。他們發(fā)現(xiàn)類星體中心黑洞質(zhì)量M?與其輻射效率呈現(xiàn)一個(gè)明顯的正相關(guān)的關(guān)系，即ε∝M?1/2。

如果這一關(guān)系是真實(shí)可靠的，那就意味著質(zhì)量較大的黑洞的自旋要比質(zhì)量小的黑洞的自旋系統(tǒng)偏大，即大質(zhì)量和小質(zhì)量黑洞的演化過(guò)程很可能不同。但是DL的PG類星體樣本只有80個(gè)類星體，且分布在有限的光度范圍內(nèi)。因此他們得到的ε和M?之間的關(guān)系可能受有限的樣本數(shù)目、較小的光度范圍等選擇效應(yīng)的影響，未必代表輻射效率與質(zhì)量真的相關(guān)。

本文首先在SDSS中選擇一個(gè)低紅移類星體樣本，其紅移值在0.025至0.5之間，與DL采用的PG類星體的紅移范圍一致。樣本的每個(gè)成員光譜均有半高全寬FHWM＞1 000 km s－1的寬線發(fā)射，樣本總數(shù)為5 786，這遠(yuǎn)遠(yuǎn)大于PG類星體的數(shù)目。針對(duì)此類星體樣本，利用標(biāo)準(zhǔn)薄盤模型估計(jì)其中每一個(gè)類星體的吸積率，并通過(guò)全波段光譜的經(jīng)驗(yàn)?zāi)０遄V估計(jì)其熱光度，從而得到每一個(gè)類星體的輻射效率，并進(jìn)一步檢驗(yàn)DL發(fā)現(xiàn)的輻射效率與黑洞質(zhì)量之間的關(guān)系是否在SDSS類星體樣本中同樣存在。

1 SDSS類星體的質(zhì)量吸積率acc的估計(jì)

從公式ε＝Lbol/(accc2)可知，要估計(jì)類星體的輻射效率，需要在一定程度上準(zhǔn)確估計(jì)質(zhì)量吸積率acc和熱光度Lbol。在標(biāo)準(zhǔn)的吸積盤模型中[20]，類星體各個(gè)波段的光譜是由中心黑洞的質(zhì)量、自旋和質(zhì)量吸積率決定的。類星體的光學(xué)波段的輻射主要由吸積盤外區(qū)的輻射主導(dǎo)，并且主要取決于吸積盤的質(zhì)量吸積率，其受自旋影響較小[19，21]。因此，可以用標(biāo)準(zhǔn)薄吸積盤模型擬合光學(xué)波段觀測(cè)來(lái)估計(jì)類星體的吸積率。

在標(biāo)準(zhǔn)吸積盤模型中，吸積盤是光學(xué)厚的，出射譜在光學(xué)紫外波段近似為多溫黑體輻射。吸積盤上半徑R處的輻射流量為:

式中，d為吸積盤距觀測(cè)者的距離；h為普朗克常量；k為波爾茲曼常量；i為觀測(cè)視線與盤面的交角；T(R)為吸積盤的表面溫度。對(duì)于標(biāo)準(zhǔn)薄盤吸積，T(R)可以近似表達(dá)為Teff＝T0R－3/4，這里:

式中，G為引力常量；σ為斯特藩－波爾茲曼常量；c為光速。將Teff代入(1)式，并做合理的近似，可以反解得到acc的值。如在DL中，采用486.1 nm處的光度，可直接得到質(zhì)量吸積率:

模型主要考慮4個(gè)參數(shù)，即黑洞質(zhì)量M?，黑洞的自旋a，質(zhì)量吸積率acc和吸積盤與觀測(cè)視線方向的夾角i。絕大多數(shù)SDSS類星體中心黑洞質(zhì)量M?已經(jīng)通過(guò)維里方法估計(jì)得到[8－9]。吸積盤的內(nèi)半徑由黑洞的自旋a決定。在高頻段吸積盤的能譜受盤的內(nèi)半徑影響較大，但在光學(xué)波段能譜則幾乎不受自旋影響。因此在本文的計(jì)算中只簡(jiǎn)單地采用一個(gè)固定的自旋值。當(dāng)前研究表明，類星體的平均輻射效率ε～0.1[2－3，14－18]，其對(duì)應(yīng)的黑洞自旋為a＝0.67，因此選定a＝0.67。在SDSS觀測(cè)中，類星體的傾角無(wú)法直接測(cè)量，但是由于SDSS數(shù)據(jù)中的類星體是Ⅰ型類星體，其軸向與觀測(cè)者視線夾角較小，一般認(rèn)為cos i在0.5到1之間均勻分布。在本文的計(jì)算中選擇固定的傾角，其余弦cos i＝0.8，這對(duì)于質(zhì)量吸積率acc的估計(jì)可能會(huì)引入一定的誤差。在給定黑洞質(zhì)量、黑洞自旋和吸積盤傾角的情況下，吸積盤的光學(xué)輻射由吸積率唯一確定。通過(guò)吸積盤模型擬合光學(xué)波段的光度Lobs，510.0nm可得到每一個(gè)類星體的質(zhì)量吸積率。

2 SDSS類星體的熱光度Lbol的估計(jì)

如果給出類星體各個(gè)波段的觀測(cè)，那么可以直接對(duì)光譜從射電到X波段積分得到其熱光度Lbol。但是，對(duì)于大多數(shù)SDSS類星體，觀測(cè)上只給出某些波段的光度，如光學(xué)波段。對(duì)于一些特殊的波段，如遠(yuǎn)紫外，其輻射往往被氣體或星系塵埃遮擋，很難觀測(cè)到。由于類星體在各個(gè)波段的輻射機(jī)制較為復(fù)雜，不可能簡(jiǎn)單地用吸積盤模型直接估計(jì)熱光度。幸運(yùn)的是對(duì)一些類星體和活動(dòng)星系核的多波段光譜研究發(fā)現(xiàn)在各個(gè)不同波長(zhǎng)段其譜分布遵從一定的冪律分布，而不同的類星體在同一波長(zhǎng)段的冪律譜指數(shù)相差很小[26－29]。因此可以構(gòu)建類星體的模板譜，從而通過(guò)一個(gè)波段或有限幾個(gè)波段的觀測(cè)估計(jì)熱光度。文[19]和文[30]中，其對(duì)熱光度的估計(jì)采取的方法本質(zhì)上與此類似。

與文[31]、文[3]等類似，本文中也采用經(jīng)驗(yàn)的方法。通過(guò)各個(gè)波段由觀測(cè)給出的經(jīng)驗(yàn)冪律譜指數(shù)分布估計(jì)不同波段的光度，并對(duì)整個(gè)波段積分得到類星體的熱光度。具體的方法如下:在光學(xué)－紫外波段，即1 μm＜λ＜130.0 nm，假設(shè)類星體的光度呈冪律分布，冪律譜譜指數(shù)呈高斯分布，均值為αO＝－0.44，方差為σαO＝0.125[26]，統(tǒng)一采用單色波段510.0 nm處的光度即L510.0nm做計(jì)算；在紫外波段，即120.0 nm＜λ＜50.0 nm，光譜也由冪律描述，高斯分布的冪律譜譜指數(shù)均值為αuv＝－1.76，方差為σαUV＝0.12[27]；在X射線波段(0.5～500 KeV)，同樣假設(shè)一個(gè)高斯分布的冪律譜譜指數(shù)，均值為Γuv＝－1.8[28－29]，方差為σΓ＝0.18。根據(jù)公式αOX＝－0.384 log[Lν，250.0nm/Lν，2keV]，可以計(jì)算出2 keV波段的光度，其中αOX由250.0 nm波段的光度決定，采用文[32]的結(jié)果αOX＝－0.154[log Lν，250.0nm/(erg s－1Hz－1)]＋3.176。假設(shè)盤的傾角為cos i＝0.8，并且盤金屬豐度為一個(gè)太陽(yáng)豐度。根據(jù)這些參數(shù)用XSPEC中的PEXRAV模型能計(jì)算出X射線波段的反射光譜[33]；50.0 nm到5.0 nm之間的光譜用冪律譜連接起來(lái)。在λ＞1 μm的紅外波段，其主要由塵埃環(huán)對(duì)中心光學(xué)紫外連續(xù)譜的再輻射，本文不將其計(jì)入熱光度估計(jì)中。

SDSS數(shù)據(jù)中提供了每個(gè)類星體光學(xué)段光度L510.0nm，通過(guò)這些譜指數(shù)可以得到各個(gè)波段的光譜，對(duì)整個(gè)光譜做積分就得到每個(gè)類星體的熱光度Lbol。根據(jù)估計(jì)的熱光度可以得到在510.0 nm處的熱改正，這里熱改正定義為Kbol＝Lbol/Lopt(如圖1)。圖1顯示的是低紅移SDSS類星體和PG類星體的熱改正，其中灰色的點(diǎn)表示低紅移SDSS樣本中5 786個(gè)類星體的熱改正，十字表示文[19]中80個(gè)PG類星體的熱改正。實(shí)線表示在文[9]、文[34]結(jié)果中固定的熱改正的值。Lopt均采用510.0 nm處的光度。從圖1可以看出低紅移SDSS類星體的熱改正不依賴于黑洞質(zhì)量，這與DL的結(jié)果非常一致。

圖1 低紅移SDSS類星體和PG類星體的熱改正KbolFig.1 The bolometric corrections Kbolfor the low-redshift SDSS QSOs and PG QSOs

3 SDSS類星體的輻射效率ε

根據(jù)前文描述的方法，通過(guò)相對(duì)論性的黑體譜模型擬合每個(gè)SDSS類星體光學(xué)波段的光度L510.0nm，從而得到其質(zhì)量吸積率acc，再根據(jù)第2節(jié)中的方法估計(jì)每個(gè)SDSS類星體的熱光度Lbol，這樣可以得到輻射效率ε＝Lbol/(accc2cos i)。同時(shí)也得到每個(gè)類星體的愛丁頓比率，即λ＝Lbol/Ledd，這里L(fēng)edd＝1.3×1046(M?/108M☉) erg s－1表示質(zhì)量為M?的黑洞對(duì)應(yīng)的愛丁頓光度。

圖2和圖3分別展示了SDSS樣本的輻射效率ε和愛丁頓比率λ隨質(zhì)量的分布。由圖2所示，顯然類星體的輻射效率與黑洞質(zhì)量強(qiáng)相關(guān)。采用Spearman相關(guān)分析，發(fā)現(xiàn)這一相關(guān)系數(shù)為rε＝0.68，相應(yīng)的此關(guān)系由隨機(jī)樣本產(chǎn)生的概率P＝0。由圖3所示類星體的愛丁頓比率與黑洞質(zhì)量的相關(guān)性也很強(qiáng)，Spearman相關(guān)系數(shù)為rλ＝－0.75，相應(yīng)的此關(guān)系由隨機(jī)樣本產(chǎn)生的概率P＝0。分別在對(duì)數(shù)空間內(nèi)對(duì)上述關(guān)系做線性擬合得到:

和

這里M8＝M?/108M☉。

圖2 SDSS類星體輻射效率ε隨質(zhì)量M?的分布Fig.2 The distribution of the low-redshift SDSS QSOs in the ε－M?plane

圖3 SDSS類星體愛丁頓比率λ隨質(zhì)量M?的分布Fig.3 The distribution of the low-redshift SDSS QSOs in the λ－M?plane

與DL結(jié)果類似，發(fā)現(xiàn)低紅移SDSS類星體的輻射效率看起來(lái)也隨黑洞維里質(zhì)量的增長(zhǎng)而增長(zhǎng)。如果這一相關(guān)是內(nèi)稟的，則可能表明黑洞的自旋也隨其質(zhì)量的增大而增大。因而較大質(zhì)量黑洞的成長(zhǎng)歷史可能與相對(duì)較小質(zhì)量的黑洞的成長(zhǎng)史有差異。但是我們注意到在類星體輻射效率的估計(jì)中涉及多個(gè)假設(shè)，這些假設(shè)均有可能引入輻射效率估計(jì)的系統(tǒng)誤差，從而有可能導(dǎo)致ε與M?間的偽相關(guān)性。下面一一討論這些假設(shè)對(duì)結(jié)果的影響:

(1)在估計(jì)吸積率時(shí)，假設(shè)維里質(zhì)量就是真實(shí)的質(zhì)量。但實(shí)際上，維里質(zhì)量相對(duì)于黑洞的真實(shí)質(zhì)量可能有0.3～0.4 dex的偏差，這一偏差可能導(dǎo)致吸積率估計(jì)的偏差。由(2)式可以看出，如果維里質(zhì)量高于真實(shí)的質(zhì)量，那么質(zhì)量吸積率會(huì)被低估，反之則被高估。在SDSS樣本中，偏高(或偏低)的黑洞維里質(zhì)量往往是被高估(或低估)的，這就會(huì)導(dǎo)致質(zhì)量吸積率比真實(shí)值偏低(或偏高)，從而使大質(zhì)量(或小質(zhì)量)樣本的輻射效率被高估(或低估)，這樣在計(jì)算的結(jié)果中就可能會(huì)出現(xiàn)低質(zhì)量處輻射效率低，高質(zhì)量處輻射效率高的趨勢(shì)。Davis和Laor(2011)也討論了維里質(zhì)量相對(duì)于黑洞質(zhì)量的偏差而造成的對(duì)輻射效率估計(jì)的影響。但是他們認(rèn)為ε與M?間的強(qiáng)相關(guān)不太可能是完全由此造成的。進(jìn)一步的考慮應(yīng)該是建立一個(gè)SDSS模擬樣本，對(duì)模擬樣本中真實(shí)黑洞質(zhì)量加入維里質(zhì)量偏差，對(duì)比加入質(zhì)量偏差前后模擬樣本質(zhì)量和輻射效率的關(guān)系。

(2)熱光度估計(jì)的不確定性。對(duì)于熱光度的估計(jì)，由于無(wú)法得到類星體每一個(gè)波段的光度，只能采用經(jīng)驗(yàn)的方法，假設(shè)SDSS各個(gè)波段的光譜均服從小樣本統(tǒng)計(jì)得來(lái)的經(jīng)驗(yàn)冪律分布，通過(guò)可觀測(cè)波段的光度推測(cè)整個(gè)光譜。這樣得到的熱改正彌散性很大。而且由于受觀測(cè)技術(shù)水平的限制，得到的樣本會(huì)有一定的選擇效應(yīng)，例如，對(duì)于愛丁頓比率較低的樣本，很可能觀測(cè)不到，導(dǎo)致我們不可能用完備的樣本分析輻射效率和質(zhì)量之間的關(guān)系。

(3)黑洞的自旋和傾角。樣本自旋和傾角對(duì)于吸積盤光譜的計(jì)算也會(huì)有影響。對(duì)于SDSS類星體很難得到樣本真實(shí)的自旋值，雖然在計(jì)算中擬合吸積率時(shí)只用光學(xué)波段的光度，且吸積盤模型中光學(xué)波段的光譜對(duì)于自旋的取值不敏感，但是簡(jiǎn)單地采用一個(gè)固定的自旋值擬合也可能會(huì)造成質(zhì)量吸積率的偏差?？紤]到SDSS類星體多數(shù)是Ⅰ型活動(dòng)星系核，所以假設(shè)樣本的傾角范圍在cos i＝[0.5～1]間均勻分布。在計(jì)算過(guò)程中采用一個(gè)固定的值cos i＝0.8估計(jì)樣本的吸積率，也可能給輻射效率的估計(jì)帶來(lái)誤差。不過(guò)總的來(lái)說(shuō)由此引入的輻射效率估計(jì)誤差不會(huì)超過(guò)30%，且這一誤差應(yīng)是隨機(jī)發(fā)散的，因而也不太可能導(dǎo)致ε與M?間的強(qiáng)相關(guān)。在考慮了傾角對(duì)質(zhì)量吸積率估計(jì)的影響后，發(fā)現(xiàn)其造成的影響不大。

(4)吸積盤模型。盡管標(biāo)準(zhǔn)的吸積盤模型在實(shí)際應(yīng)用中獲得很大的成功，但實(shí)際的黑洞吸積過(guò)程較復(fù)雜，通過(guò)簡(jiǎn)單吸積盤模型計(jì)算得到的質(zhì)量吸積率很可能會(huì)有系統(tǒng)偏差。例如，多溫黑體譜模型估計(jì)的吸積率相對(duì)于由TLUSTY模型估計(jì)的吸積率要大1.5倍，這就表明由于吸積盤模型的選取導(dǎo)致輻射效率的估計(jì)可能會(huì)系統(tǒng)偏大或偏小1.5倍左右。但這種系統(tǒng)偏差也應(yīng)該不會(huì)導(dǎo)致ε與M?間的強(qiáng)相關(guān)。

(5)樣本的選擇效應(yīng)。文[30]作者認(rèn)為DL工作中PG類星體的輻射效率和質(zhì)量的相關(guān)可能是由于樣本的選擇效應(yīng)導(dǎo)致的。他們認(rèn)為采用分布在更廣光度范圍的大樣本，可能會(huì)得到非常不同的結(jié)果。除此以外，輻射效率的估計(jì)還受很多其它因素的影響，如宿主星系、塵埃以及熱光度估計(jì)的不確定性等。從本文的研究來(lái)看，單純地加大樣本數(shù)目不會(huì)減弱DL發(fā)現(xiàn)的ε與M?間強(qiáng)相關(guān)。

在后續(xù)研究中，將進(jìn)一步建立一個(gè)SDSS類星體的模擬樣本，從而可以有效地考慮上述各種假設(shè)引入的對(duì)輻射效率估計(jì)的影響。根據(jù)此模擬樣本與真實(shí)SDSS類星體樣本的比對(duì)研究可以區(qū)分以上各種不同的偏差對(duì)本文發(fā)現(xiàn)的ε與M?間的強(qiáng)相關(guān)貢獻(xiàn)的多少，從而回答究竟輻射效率是否會(huì)隨黑洞質(zhì)量增大而增大這個(gè)問(wèn)題。也將輻射效率的估計(jì)擴(kuò)展至所有的SDSS類星體，并考察類星體的輻射效率是否隨紅移演化，從而幫助理解不同質(zhì)量黑洞的吸積歷史以及大質(zhì)量黑洞的宇宙學(xué)演化。

4 結(jié) 論

本文研究了低紅移SDSS類星體(0.025＜z＜0.5)的輻射效率，利用樣本提供的黑洞質(zhì)量，通過(guò)薄吸積盤模型估計(jì)了每一個(gè)類星體的吸積率。通過(guò)光學(xué)波段的觀測(cè)光度，利用經(jīng)驗(yàn)的類星體全波段模板譜估計(jì)了每個(gè)類星體的熱光度。根據(jù)估計(jì)的吸積率和熱光度得到了每個(gè)類星體的輻射效率。發(fā)現(xiàn)低紅移SDSS類星體的輻射效率與黑洞質(zhì)量強(qiáng)相關(guān)，并滿足ε∝M?0.63。這一結(jié)果與文[19]對(duì)PG類星體樣本研究得到的結(jié)果一致。進(jìn)一步討論由輻射效率估計(jì)中各種假設(shè)可能引入的偏差及其對(duì)ε與M?關(guān)系的貢獻(xiàn)。將通過(guò)進(jìn)一步的系統(tǒng)研究揭示ε與M?之間的相關(guān)是否是內(nèi)稟的。

[1] Shakura N I，Sunyaev R A.Black holes in binary systems.Observational appearance[J]. Astronomy＆Astrophysics，1973，24(1):337－355.

[2] Yu Q，Tremaine S.Observational constraints on growth of massive black holes[J].Monthly Notices of the Royal Astronomical Society，2002，335:965－976.

[3] Marconi A，Risaliti G，Gillin R，et al.Local super massive black holes，relics of active galactic nuclei and the X-ray background[J].Monthly Notices of the Royal Astronomical Society，2004，351:169－185.

[4] Gammie C F，Shapiro S L，McKinney J C.Black hole spin evolution[J].The Astrophysical Journal，2004，602:L312－319.

[5] Volonteri M，Madau P，Quataert E，et al.The distribution and cosmic evolution of massive black hole spins[J].The Astrophysical Journal，2005，620(1):69－77.

[6] Kaspi S，Smith P S，Netzer H，et al.Reverberation measurements for 17 quasars and the size-massluminosity relations in active galactic nuclei[J].The Astrophysical Journal，2000，533(2):631－649.

[7] Peterson B M，F(xiàn)errarese L，Gilbert K M，et al.Central masses and broad line region sizes of active galactic nuclei.II.A homogeneous analysis of a large reverberation-mapping database[J]. The Astrophysical Journal，2004，613(2):682－699.

[8] Shen Y，Greene J F，Strauss M A，et al.Biases in virial black hole masses:An SDSS pers pective[J].The Astrophysical Journal，2008，680(1):169－190.

[9] Shen Y，Richards G T，Strauss M A，et al.A catalog of quasar properties from sloan digital sky survey data release 7[J].The Astrophysical Journal Supplement Series，2011，194:45－63.

[10] Vestergaard M.Determining central black hole masses in distant active galaxies[J].The Astrophysical Journal，2002，571(2):733－752.

[11] Vestergaard M，Peterson B M.Determining central black hole masses in distant active galaxies and quasars.II.Improved optical and UV scaling relationships[J].The Astrophysical Journal，2006，641(2):689－709.

[12] Fabian A C，Iwasawa K，Reynolds C S，et al.Broad iron lines in active galactic nuclei[J]. Publications of the Astronomical Society of the Pacific，2000，112:1145－1161.

[13] Lynden-Bell D.Galactic nuclei as collapsed old quasars[J].Nature，1969，223:690－694.

[14] Yu Q，Lu Y.Constraints on QSO models from a relation between the QSO luminosity function and the local black hole mass function[J].The Astrophysical Journal，2004，602(2):603－624.

[15] Yu Q，Lu Y.Toward precise constraints on the growth of massive black holes[J].The Astrophysical Journal，2008，689(2):732－754.

[16] Barger A J，Cowie L L，Mushotzky R F，et al.The cosmic evolution of hard X-Ray-selected active galactic nuclei[J].The Astronomical Journal，2005，129(2):578－609.

[17] Shankar F，Weinberg D H，Miralda-Escudé J.Self-consistent models of the AGN and black hole populations:duty cycles，accretion rates，and the mean radiative efficiency[J].The Astrophysical Journal，2009，690(1):20－41.

[18] Elvis M，Risaliti G，Zamorani G.Most supermassive black holes must be rapidly rotating[J]. The Astrophysical Journal，2002，565(2):L75－L77.

[19] Davis S W，Laor A.The radiative efficiency of accretion flows in individual active galactic nuclei [J].The Astrophysical Journal，2011，728(1):98－117.

[20] Novikov I D，Thorne K S.Black holes(Les Astres Occlus)[M].New York:Gordon and Breach，1973:343－450.

[21] Krolik J H.Active Galactic Nuclei:From the central black hole to the galactic environment [M].Princeton:Princeton University Press，1999.

[22] Gierliński M，Maciolek-Niedzwiecki A，Ebisawa K.Application of a relativistic accretion disc model to X-ray spectra of LMC X-1 and GRO J1655-40[J].Monthly Notices of the Royal Astronomical Society，2001，325:1253－1265.

[23] Page D N，Thorne K S.Disk-Accretion onto a Black Hole.Time-averaged structure of accretion disk[J].The Astrophysical Journal，1974，191(2):499－506.

[24] Hubeny I，Lanz T.Non-LTE line-blanketed model atmospheres of hot stars.1:Hybrid complete linearization/accelerated lambda iteration method[J].The Astrophysical Journal，1995，439(2):875－904.

[25] Hubeny I，Agol E，Blaes O，et al.Non-LTE models and theoretical spectra of accretion disks in active galactic nuclei.III.Integrated spectra for hydrogen-helium disks[J].The Astrophysical Journal，2000，533(2):710－728.

[26] Vanden Berk D E，Richards G T，Bauer A，et al.Composite quasar spectra from the SloanDigital Sky Survey[J].The Astronomical Journal，2001，122(2):549－564.

[27] Telfer R C，Zheng W，Kriss G A，et al.The rest-frame extreme-ultraviolet spectral properties of quasi-stellar objects[J].The Astrophysical Journal，2002，565(2):773－785.

[28] Jin C C，Ward M，Done C，et al.A combined optical and X-ray study of unobscured type I active galactic nuclei-I.Optical spectra and spectral energy distribution modelling[J].Monthly Notices of the Royal Astronomical Society，2012，420(3):1825－1847.

[29] Tozzi P，Gilli R，Mainieri V，et al.X-ray spectral properties of active galactic nuclei in the Chandra Deep Field South[J].Astronomy＆Astrophysics，2006，451(2):457－474.

[30] Raimundo S I，F(xiàn)abian A C，Vasudevan R V，et al.Can we measure the accretion efficiency of active galactic nuclei?[J].Monthly Notices of the Royal Astronomical Society，2012，419(3): 2529－2544.

[31] Hopkins P F，Richards G T，Hernquist L.An observational determination of the bolometric quasar luminosity function[J].The Astrophysical Journal，2007，654(2):731－753.

[32] Lusso E，Comastri A，Vignali C，et al.The X-ray to optical-UV luminosity ratio of X-ray selected type 1 AGN in XMM-COSMOS[J].Astronomy＆Astrophysics，2010，2512(1):34－48.

[33] Magdziarz P，Zdziarski A A.Angle-dependent Compton reflection of X-rays and gamma-rays [J].Monthly Notices of the Royal Astronomical Society，1995，273(2):837－848.

[34] Richards G T，Lacy M，Storrie-Lombardi L J，et al.Spectral energy distributions and multiwavelength selection of type 1 quasars[J].The Astrophysical Journal Supplement Series，2006，166(2):470－497.

Estimates of Radiation Efficiencies of Low-Redshift SDSS QSOs

Wu Shumei，Zhang Fupeng，Lu Ye，Lu Youjun
(National Astronomical Observatories，Chinese Academy of Sciences，Beijing 100012，China，Email:shmwu＠bao.ac.cn)

Radiation efficiencies of QSOs are believed to be determined by the angular momenta of their central Massive Black Holes(MBH).Assembly histories of MBH may have shaped their angular momenta.In this paper we estimate radiation efficiencies of selected low-redshift SDSS QSOs(0.025＜z＜0.5).Using estimated Virial-mass values to approximate the MBH mass(M?)values，we fit predictions of the standard model of thin accretion disks to observed optical luminosities of the QSOs to individually estimate their accretion rates.The bolometric luminosity of each QSO is estimated by extrapolating its observed optical luminosities in different bandpasses using empirical panchromatic templates established through various observations.With the estimated accretion rate and bolometric luminosity we obtain the radiation efficiency of each QSO.We find a strong correlation between the radiation efficiency(ε)and the black-hole mass，ε∝M?0.63，which is similar to the result of Davis＆Laor(2011).The estimated radiation efficiency of an individual MBH has a large error which includes errors propagated from various sources.We further discuss various possible biases involved in the estimates of the radiation efficiencies，as these biases potentially contribute to the result ε-M?relation.These include the uncertainty of approximating the MBH mass with the estimated Virial mass，the oversimplification of the accretion-disk model，the sample selection effects，and the errors of the estimated bolometric luminosities due to the limited database for templates and unknown dust/gas obscuration.

Accretion disk；Radiation efficiency；Black hole；Accretion rate；Spin

P142.9

A

1672－7673(2014)02－0095－07

2013－04－28；修定日期:2013－06－24

吳淑梅，女，碩士.研究方向:類星體與活動(dòng)星系核.Email:shmwu＠bao.ac.cn