BL Lac天體ON 231的光變周期研究*

徐云冰，張 雄，李孝攀，董富通

(1.云南師范大學(xué)物理與電子信息學(xué)院，昆明 650092；2.中國科學(xué)院國家天文臺，北京 100012)

ON 231是活動激烈的BL Lac天體，天體的周期性光變現(xiàn)象可以給出天體內(nèi)部的一些參數(shù)，進而可以研究其內(nèi)部物理機制[1]。因此，人們對ON 231周期性的研究具有濃厚興趣，很多人已經(jīng)對其光變周期進行了研究并取得了很多成果。1987年，Xie等人發(fā)現(xiàn)了約6天的周期[2]，這一結(jié)果與OJ 287不存在幾天短時標光變的結(jié)果不一致[3]；1988年，Barbieri等人用Deeming(1975)使用的功率譜分析方法沒有發(fā)現(xiàn)ON 231有明顯的周期[4-5]；1995年，Liu F K等人找到了ON 231有(0.981±0.005)年、(13.6±1.3)年的周期[3]；1998年，Zhang等人找到了(13.6±1.3)年、(26.6±0.8)年的周期[1]。但是，由于ON 231內(nèi)部結(jié)構(gòu)復(fù)雜且活動劇烈，其光變曲線也很復(fù)雜，另外，由于觀測手段和天氣等因素的限制，目前還沒有很完備的長期樣本數(shù)據(jù)庫，因此，對ON 231周期性的明確認證總體來說還比較困難，其長周期光變有可能大于現(xiàn)有的觀測歷史，有待于進一步監(jiān)測認證[1]。

本文主要利用Jurkevich方法[6]和小波分析法[7-8]研究了ON 231的歷史光變資料中可能存在的周期。

1 ON 231的光變周期性分析

1.1 ON 231的光變曲線

ON 231的測光資料來源于文獻[4，9-29]，共收集了1374個B波段觀測數(shù)據(jù)，獲得了其天平均歷史光變曲線如圖1。從光變曲線可以看出天體活動十分劇烈，在1998年最亮?xí)r達mB=13.1mag，1984年最暗時達mB=17.80mag，在近110多年的觀測中出現(xiàn)過多次大的爆發(fā)，從天平均的歷史光變曲線很難看出明顯的周期變化。

1.2 ON 231周期性分析

ON 231的觀測資料表明它內(nèi)部確實存在劇烈變化。它的長期光變周期一直沒有最后確認，觀測已經(jīng)發(fā)現(xiàn)了很多次大的爆發(fā)結(jié)果。從這些結(jié)果可以看出，它有可能存在光變的周期性，有許多人用多種方法研究周期性,但由于計算方法的缺陷和觀測數(shù)據(jù)的連續(xù)性差,確定的周期均未被觀測驗證,選用小波分析法研究其周期性,用Jurkevich方法進行比較,獲得了較好的結(jié)果。當然，它的實際周期也有可能大于現(xiàn)有的觀測資料[1]。

圖1 ON 231 1898～2005年天平均的歷史光變曲線

1.2.1 Jurkevich方法

如果f≥0.5，則表明樣本數(shù)據(jù)有很強的周期性；如果f<0.25，則表明樣本數(shù)據(jù)可能有弱的周期性；如果f=0，則表明樣本數(shù)據(jù)沒有周期性。

圖2 用圖1中的數(shù)據(jù)由Jurkevich方法尋找ON 231周期結(jié)果

1.2.2 小波分析法

小波分析法[7-8]是為了克服Fourier分析的不足而發(fā)展起來的，F(xiàn)ourier分析反映的是信號的整體特征，但不能刻畫信號的局部特性。小波分析法是一種時頻分析方法[31]，既能在時域很好地刻畫信號的局部性，同時也能在頻域反映信號的局部性。它的基本思想是通過小波變換對信號進行多尺度分析，從而突出信號的局部特征。該方法已廣泛應(yīng)用于地震勘探、氣象學(xué)等方面的周期性研究。

選用Morlet小波：ψ(t)=e-t2/2eiωt，容許性條件要求ω≥5。根據(jù)Morlet小波尺度因子a與周期T之間的關(guān)系[32]：

取ω=6.2，則周期T≈a。利用離散小波變換求出小波系數(shù)，然后給出小波系數(shù)的實部等值線圖，等值線的變化即可反映出所分析數(shù)據(jù)在小波域中的波動情況。對時間域上的尺度因子a積分得到小波方差，再畫出小波方差圖。通過該圖可以對一些復(fù)雜的非線性過程進行解釋，從而確定其中的主要周期。

ON 231 1962年以前的數(shù)據(jù)點很少，只取了1962年以后的數(shù)據(jù)進行分析。為了平滑觀測誤差，對該時間段中的數(shù)據(jù)進行30天平均，再利用線性插值的方法對不等間隔的數(shù)據(jù)進行填充[8]，以消除邊界效應(yīng)。圖3給出了ON 231 1962～2005年B波段數(shù)據(jù)小波變換系數(shù)實部等值線圖。圖中顯示星等值在小波域中高低起伏的變換情況，即反映了ON 231在時間域上的總體光變過程。

圖3 ON 231 B波段數(shù)據(jù)小波變換系數(shù)實部等值線圖

從圖中等值線的變化可以看出在周期尺度T=4000～6000天之間和T=7500～10500天之間存在明顯的周期性。為了確定其中的主要周期,作出其變換系數(shù)的小波方差圖如圖4。

圖4 變換系數(shù)的小波方差圖

從圖4可以看出，在T=4940天和T=9520天處達到極大值，說明這兩個值就是分析信號中含有的周期。T=4940天(13.5年),T=9520天(26.1年)與Liu F K等人[3]、Zhang等人[1]26.6±0.8年的結(jié)果基本一致。

1.3 中心黑洞質(zhì)量估計

2 討 論

作者認為，時標為幾年的光變特性可能是由于通常在活動星系核中心區(qū)黑洞周圍的吸積盤不穩(wěn)定所致，特別是薄盤的熱不穩(wěn)定性可能導(dǎo)致爆發(fā)式振蕩[36]；雙黑洞模型不是解釋BL Lac天體長周期光變最理想的模型,長周期爆發(fā)不可能完全起源于雙黑洞[1]。

3 結(jié) 論

綜上所述，ON 231具有0.978年、4940天(13.5年)和9520天(26.1年)的周期。其中T=0.978年與Liu F K等人[3]的0.981±0.005年一致；T=13.5年和T=26.1年，與Liu F K等人[3]、Zhang等人[1]的13.6±1.3年和26.6±0.8年的結(jié)果一致。

[1] 張雄,謝光中,白金明.用Jurkevich方法計算BL Lac天體的光變周期[J].天體物理學(xué)報,1998,18(3):260-268.

Zhang Xiong,Xie Guangzhong,Bai Jinming.Computing the Period of Light Variability in BL Lac Objects Using the Jurkevich Method[J].Acta Astrophysica Sinica,1998,18(3):260-268.

[2] Xie G Z,Li K H,Bao M X,et al.The optical variability of seven BL Lacertae objects[J].A &AS,1987,67:17-24.

[3] Liu F K,Xie G Z,Bai J M.An historical light curve of ON 231 and its periodic analysis[J].A&A,1995,295:1-10.

[4] Barbieri C,Cappellaro E,Romano G,et al.Optical variability of 16 quasi stellar objects[J].A&AS,1988,76:477-487.

[5] Deeming T J.Fourier analysis with unequally-spaced data[J].APSS,1975,36(1):137-158.

[6] Jurkevich I N.A method of computing periods of cyclic phenomena[J].APSS,1971,13(1):154-167.

[7] Christopher Torrence,Gilbert P Compo.A Practical Guide to Wavelet Analysis[J].Bulletin of the American Meteorological Society,1998,79(1):61-78.

[8] 董富通,張雄.類星體3C345的光變周期特性[J].物理學(xué)報,2009,58(11):8116-8122.

Dong Fu-Tong,Zhang Xiong.Analysis of the variability periodicities of quasar 3C345[J].Acta physica sinica,2009,58(11):8116-8122.

[9] Craine E R,Warner J W.Optical variations of OJ 287,ON 231,and OQ 208[J].ApJ,1973,179:L53-L56.

[10] Wing R F.UBV sequences in the fields of radio sources[J].AJ,1973,(78):684-686.

[11] Pollock J T,Hall D J,Ambruster C,et al. A Sustained Pulse from ON 231[J].A&A,1974,30:41-43.

[12] McGimsey B Q,Smith A G,Scott R L,et al.Optical behavior of 20 violently variable extragalactic radio sources[J].AJ,1975,80:895-922.

[13] V’eron P,V’eron M P.Photographic Photometry of Five BL Lacertae-type Objects[J].A&A,1975,39:281-288.

[14] Tapia S,Craine E R,Jahnson K.Four-point optical energy distributions for faint BL lacertae objects[J].ApJ,1976,203:291-296.

[15] Pollock J T,Pica A J,Smith A G,et al. Long-term optical variations of 20 violently variable extragalactic radio sources[J].AJ,1979,84:1658-1676.

[16] Barbieri C,Romano G.Survey of the optical variability of compact extragalactic objects.IV.Objects from 12#to 16#[J].A&AS,1981,44:159-162.

[17] Schaefer B E.A search for rapid variability in eight quasars and BL Lacertae objects[J].PASP,1980,92:255-258.

[18] Zekl H,Klare G,Appenzeller I.Optical Brightness Variations of BL Lac Objects[J].A&A,1981,103:342-348.

[19] Gruze-Gonzalez I,Huchra J P.Continuum distributions of an x-ray observed sample of BL Lac objects[J].AJ,1984,89:441-465.

[20] Webb J R,Smith A G,Leacock R J,et al.Optical observations of 22 violently variable extragalactic sources:1968-1986[J].AJ,1988,95:374-397.

[21] Xie G Z,Li K H, Bao M X,et al.The optical variability of seven BL Lacertae objects[J].A&AS,1987,67:17-24.

[22] Xie G Z,Li K H,Zhou Y,et al. The optical variability of two x-ray-selected and nine radio-selected BL Lac objects[J].AJ,1988,96:24-29.

[23] Xie G Z,Li K H,Cheng F Z，et al.Search for short variability time-scales of BL Lacertae objects[J].A&A,1990,229(2):329-339.

[24] Xie G Z,Li K H,Cheng F Z，et al.CCD photometry of 10 BL Lacertae objects[J].A&AS,1991,87(3):461-469.

[25] Xie G Z,Li K H,Liu F K，et al.CCD photometry of 14 BL Lacertae objects and theoretical model[J].ApJ S,1992,80(2):683-699.

[26] Tosti G,Fiorucci M,Luciani M,et al.BVRcIcmonitoring of ON 231 during the great outburst in 1994-1997[J].A&AS,1998,130:109-115.

[27] Fan J H,Lin R G.Optical variability and periodicity analysis for blazars I Light curves for radio-selected BL Lacertae objects[J].ApJ,2000,537:101-122.

[28] Tosti G,Massaro E,Nesci R, et al. The optical behaviour of ON 231(W Comae)during and after the great outburst of spring 1998[J].A&A,2002,395(1):11-15.

[29] Babadzhanyants M K,Belokon’ E T.Optical monitoring of the BL Lac object ON 231 in 1972-1990[J].ARep,2002,46(8):609-625.

[30] Kidger M,Takalo L,Sillanpaa A.A new analysis of the 11-year period in OJ287:Confirmation of its existence[J].A&A,1992,264(1):32-36.

[31] 彭玉華.小波變換與工程應(yīng)用[M].北京：科學(xué)出版社, 2002.

[32] Meyers S D,Kelly B G,O’Brien J J.An Introduction to Wavelet Analysis in Oceanography and Meteorology:With Application to the Dispersion of Yanai Waves[J].Mon Wen Rev,1993,121(10):2858-2866.

[33] Wallinder F H,Abramowicz M A,Kato S.Variability of the central region in active galactic nuclei[J].A&AR,1992,4(2):79-122.

[34] Honma F,Kato S,Matsunioto R.Nonlinear oscillations of thermally unstable slim accretion disks around a neutron star or a black hole[J].PASJ,1991,43(1):147-168.

[35] Horiuch T,Kato S.A model of hydromagnetic turbulent viscosity in radiation-pressure-dominated disks[J].PASJ,1990,42(5):661-674.

[36] 樊軍輝,G E Romero,林瑞光.3C273的光學(xué)光變周期[J].天文學(xué)報,2001,42(1):9-15.

Fan Junhui,G E Romero,Lin Ruiguang.The Optical Variability Periodicity Analysis of 3C273[J].ACTA Astronomica Sinica,2001,42(1):9-15.