董愛軍

(滁州學(xué)院機(jī)電學(xué)院，安徽 滁州 239000)

通過功率譜分析，人們發(fā)現(xiàn)黑洞候選體X射線雙星Cyg X-1存在420 d和150 d的光變周期［1］。雖然功率譜分析法是當(dāng)前處理天體光變周期行之有效的方法之一，但是，由于該方法只能處理等時(shí)間間隔的觀測(cè)數(shù)據(jù)，分析天體的光變周期往往需要進(jìn)行插值或曲線擬合，這必然帶來一些不真實(shí)的信息，增加了光變周期準(zhǔn)確性認(rèn)證的困難。Jurkevich法［2］是一種建立在期望值均方誤差基礎(chǔ)上的頻譜分析方法，該方法最大的優(yōu)點(diǎn)是能夠直接處理非等時(shí)間間隔的觀測(cè)數(shù)據(jù)，真實(shí)反應(yīng)天體光變曲線的周期特征［3－4］，而且利用該方法已經(jīng)成功討論了許多天體的光變周期［5－12］。另外，Kidger等［1，3］引入的表征周期特性可信度的Kidger判據(jù)，擴(kuò)大了Jurkevich法的適用范圍。

為了討論黑洞候選體X射線雙星Cyg X-1的光變周期，本文利用Jurkevich法和Kidger判據(jù)，分析了黑洞候選體X射線雙星Cyg X-1近16年Rossi X-ray Timing Explorer(RXTE)All-Sky Monitor(ASM)觀測(cè)資料。本文第1節(jié)給出數(shù)據(jù)處理方法，第2節(jié)給出分析結(jié)果，第3節(jié)對(duì)結(jié)果進(jìn)行了討論，文章的最后給出了本文的結(jié)論。

1 數(shù)據(jù)處理

1.1 Jurkevich 法

Jurkevich法［2］是一種基于期望值均方誤差的統(tǒng)計(jì)方法，它把觀測(cè)數(shù)據(jù)按試驗(yàn)周期折疊，分為m組，并對(duì)每組求均方誤差，然后取總方差。若試驗(yàn)周期等于實(shí)際周期，則出現(xiàn)最小值［2，4］。

對(duì)于某段光變曲線{xi，i=1，2，3，…，N}，利用Jurkevich法分析其光變周期主要包括以下幾個(gè)步驟:

(1)根據(jù)試驗(yàn)周期p把光變曲線{xi，i=1，2，3，…，N}分成m組

給定觀測(cè)點(diǎn)進(jìn)入的分組索引號(hào)為:

式中，char是取整運(yùn)算;p為試驗(yàn)周期;t0為觀測(cè)數(shù)據(jù)起始時(shí)間。則組索引號(hào)的取值范圍NDGR={0，1，2，…，m－1}。

則m組的總方差:

(3)歸一化

由(1)和(2)步可知若光變曲線{xi，i=1，2，3，…，N}存在某一光變周期T，則在p－圖上，p=T處，應(yīng)出現(xiàn)最小值。

(5)誤差分析

本文取峰值的半最大值半寬度(Half Width at Half Maxium，HWHM)作為對(duì)應(yīng)光變周期的誤差范圍。

1.2 Kidger判據(jù)

1992年Kidger給出了 Jurkevich法判別p－圖中光變周期真?zhèn)涡缘?Kidger判據(jù)［13］:若為歸一化的，則

f≥0.5:光變曲線存在很強(qiáng)的周期特性; f≤0.25:光變曲線存在較弱的周期特性;f=0:光變曲線不存在周期特性。

1.3 分組m及觀測(cè)時(shí)間跨度對(duì)觀測(cè)結(jié)果的影響

在計(jì)算過程中，m的選擇十分重要，m越大靈敏度越高，但會(huì)產(chǎn)生越大的噪音;m太小則找不到相應(yīng)的光變周期［2］。同時(shí)，觀測(cè)時(shí)間跨度對(duì)光變周期的確定影響也十分大，一般情況下，觀測(cè)時(shí)間跨度要大于6倍有效周期［13］，但是如果數(shù)據(jù)分布較好，1.5個(gè)周期的觀測(cè)時(shí)間跨度就可以了［4，14］。

2 結(jié)果與分析

利用Jurkevich法分析處理了黑洞候選體X射線雙星 Cyg X-1 MJD=50087－55841期間的ASM觀測(cè)資料，結(jié)果表明，Cyg X-1存在明顯的光變周期特性。利用Jurkevich法分析光變曲線的周期特性時(shí)，分組數(shù)m的選擇通常會(huì)影響光變周期特性的認(rèn)證。為了準(zhǔn)確測(cè)定Cyg X-1的光變周期，本文取m=(10，15，20)，試驗(yàn)周期p從1 d到1200 d，步長(zhǎng)為1 d。

圖1給出了Cyg X-1 MJD=50087－55841期間軟X射線波段的光變曲線以及硬度比(Hardness Ratio，HR)的時(shí)變特性。從圖中可以看出，Cyg X-1光子計(jì)數(shù)率變化十分劇烈，經(jīng)歷了數(shù)次態(tài)的轉(zhuǎn)換，但光變周期特性并不顯著。傳統(tǒng)上，一般采用功率譜法和相位分析法等手段分析其光變曲線的特征頻率。這些方法處理具有固定步長(zhǎng)的離散觀測(cè)數(shù)據(jù)的光變周期較為方便，但對(duì)于非均勻采樣的觀測(cè)數(shù)據(jù)，往往需要通過插值法或曲線擬合，這必然帶來一些不真實(shí)的信息，增加了光變周期準(zhǔn)確性的認(rèn)證困難。利用Jurkevich法很好地解決了這些問題。圖2給出了不同m時(shí)，試驗(yàn)周期p與總方差的關(guān)系曲線。從圖中可以看出，總是在某些試驗(yàn)周期處出現(xiàn)最小值。

圖1 1996年1月至2011年10月黑洞候選體X射線雙星1.5－12 KeV能段上，Cyg X-1的光變曲線以及對(duì)應(yīng)的硬度比HR。上圖:光變曲線，下圖:硬度比HRFig.1 The Light Curve and Hardness Ratios of the black-hole candidate X-ray binary Cyg X-1 from January 1996 to October 2011 in the 1.5－12KeV energy band

從圖2和表1可以看出，在m=(10，15，20)3種情況下，黑洞候選體X射線雙星Cyg X-1的光變曲線都存在363±11.5 d的光變周期。同時(shí)，在試驗(yàn)周期為728±24 d和1120±42 d時(shí)，也出現(xiàn)極小值，而且周期特性十分明顯(f＞0.5)。由于該周期近似為363±11.5 d的倍數(shù)關(guān)系，因此可能是363±11.5 d周期的觀測(cè)表象。另外，從圖2和表1還可以發(fā)現(xiàn)，Cyg X-1存在51±1.4 d的光變周期，其Kidger判據(jù)0.25≤f≤0.5，同時(shí)還發(fā)現(xiàn)了該光變周期的倍數(shù)周期(102±0.9 d、153±0.9 d以及314±5 d)。

表1 不同m值，試驗(yàn)周期p對(duì)應(yīng)的Kidger判據(jù)fTable 1 The trial period p and the Kidger f for different values of m

圖2 m=(10，15，20)時(shí)，試驗(yàn)周期步長(zhǎng)為1 d的p－曲線Fig.2 The p－curves for m=(10，15，20)and step of tiral period=1 day

3 討論

文［1］作者利用1996～2001年Cyg X-1處于硬態(tài)時(shí)的ASM和PCA資料，計(jì)算分析了該天體光變曲線的功率密度譜，發(fā)現(xiàn)Cyg X-1存在p=150 d和420 d的光變周期(圖3)。通過對(duì)不同波段功率密度譜分析，他們認(rèn)為150 d的光變周期只是光度變化的一種反映，而420 d的光變周期是軟X射線流的內(nèi)在反映，同時(shí)影響著X射線譜的變化。但是功率密度譜通常只能處理固定時(shí)間步長(zhǎng)的觀測(cè)數(shù)據(jù)，在處理不等時(shí)間步長(zhǎng)的觀測(cè)數(shù)據(jù)時(shí)，往往需要進(jìn)行插值或曲線擬合，這必然會(huì)影響信息的真實(shí)性。Jurkevich法是一種建立在期望值均方誤差基礎(chǔ)上，能夠處理非等時(shí)間間隔數(shù)據(jù)的頻譜分析方法，能夠充分反映觀測(cè)數(shù)據(jù)的真實(shí)性。

Kidger［13］指出，利用Jurkevich法計(jì)算光變曲線的周期特性時(shí)，數(shù)據(jù)的最小時(shí)間跨度為6倍有效周期長(zhǎng)，圖2給出的51±1.4 d和363±11.5 d的光變周期滿足這個(gè)條件。

從圖3可以看出，(1)Cyg X-1存在420 d的光變周期，但是誤差范圍較大(HWHM≈120 d);(2)功率譜分析法容易帶來不真實(shí)的信息。利用Jurkevich法分析1997年12月～2001年7月間的ASM資料及圖3［1］中采用的數(shù)據(jù)，MJD=50786－52121)的結(jié)果表明(圖4)，當(dāng)滿足Kidger最小時(shí)間跨度時(shí)，得到了與圖2(b)近似一致的結(jié)果，不滿足時(shí)，最小值明顯發(fā)散。另外，從圖4還可以看出，p≈363和p≈420時(shí)出現(xiàn)了最小值。由于363 d和420 d不能滿足Kidger最小時(shí)間跨度要求，二者都有可能是真實(shí)光變周期處最小值的發(fā)散導(dǎo)致的。從圖2可以看出，Cyg X-1的光變曲線存在363±11.5 d的光變周期及其倍數(shù)周期，但420 d的光變周期并沒有發(fā)現(xiàn)。鑒于以上原因認(rèn)為420 d的光變周期可能是觀測(cè)數(shù)據(jù)分布帶來的，而不能準(zhǔn)確反映Cyg X-1的光變特性，363±11.5 d才是Cyg X-1的光變周期。

圖3 Cyg X-11997年12月至2001年7月間的功率密度譜［1］Fig.3 The Power Spectrum Density for Cyg X-1 during December 1997 to July 2001

圖4 m=15時(shí)，試驗(yàn)周期步長(zhǎng)為1 d的p－曲線Fig.4 The p－curve for m=15 and step of tiral period=1 day

圖5給出了試驗(yàn)周期從1 d到1200 d，步長(zhǎng)為10 d的p－曲線。從圖中可以看出Cyg X-1存在51±1.4 d的光變周期，但是150 d的光變周期并不存在。因此認(rèn)為圖2出現(xiàn)的51±1.4 d的光變周期才是Cyg X-1真實(shí)的光變周期，而102±0.9 d、153±0.9 d和314±5 d的光變周期只是其觀測(cè)表象而已。

圖5 m=15時(shí)，試驗(yàn)周期步長(zhǎng)為10 d的p－曲線Fig.5 The p－curve for m=15 and step of tiral period=10 day

4 結(jié)論

Jurkevich法是一種處理非等時(shí)間間隔觀測(cè)數(shù)據(jù)十分有效的頻譜分析法。本文利用該方法發(fā)現(xiàn)黑洞候選體X射線雙星存在363±11.5 d的光變周期。同時(shí)認(rèn)為，利用功率譜得到的150 d的光變周期，只不過是51±1.4 d光變周期的觀測(cè)表象。

［1］Benlloch S，Pottschmidt K，Wilms J，et al.Different kinds of long-term variability from cygnus X-1 ［J］.AIP Conference Proceedings，2004，714:61－64.

［2］Jurkevich I.A method of computing periods of cyclic phenomena ［J］.Astrophysics and Space Science，1971，13(1):154－167.

［3］Babadzhanyants M K，Belokon E T.In variability blazars［M］.Cambridge:cambridge university press，1991:384.

［4］Fan J H，Xie G Z，Pecontal E，et al.Historic light curve and long-term optical variation of BL lacertae 2200+420 ［J］.The Astrophysical Journal，1998，507(1):173－178.

［5］Jurkevich I，Usher P D，Shen B S P.On the light variations of seyfert galaxy 3C120 ［J］.Astrophysics and Space Science，1971，10(3):402－410.

［6］樊軍輝，蘇成悅.NGC 4151的長(zhǎng)期光變行為 ［J］.天文學(xué)報(bào)，1998，39(3):251－255.Fan Junhui，Su Chengyue.The long-term behavior of optical variability of NGC 4151 ［J］.Acta Astronomica Sinica，1998，39(3):251－255.

［7］張雄，謝光中，白金明.類星體3C345的光變周期特性研究 ［J］.天文學(xué)報(bào)，1998，39(3):256－264.Zhang Xiong，Xie Guangzhong，Bai Jinming.A study of the light curve periodic behavior of quasar 3C345 ［J］.Acta Astronomica Sinica，1998，39(3):256－264.

［8］樊軍輝，Romero G E，林瑞光.3C273的光學(xué)光變周期 ［J］.天文學(xué)報(bào)，2001，42(1):9－15.Fan Junhui，Romero G E，Lin Ruiguang.The optical variability periodicity analysis of 3C273［J］.Acta Astronomica Sinica，2001，42(1):9－15.

［9］陶雋，錢伯辰，梁兆剛，等.判定光變周期的一種方法 ［J］.中國(guó)科學(xué)院上海天文臺(tái)年刊，2004，25:97－103.Tao Jun，Qian Bochen，Liang Zhaogang，et al.A method used for finding the period of the light curve ［J］.Annals of Shanghai Observatory Academia Sinica，2004，25:97－103.

［10］毛李勝.3C120射電光變曲線中的一個(gè)可能周期 ［J］.天文研究與技術(shù)——國(guó)家天文臺(tái)臺(tái)刊，2007，4(4):307－312.Mao Lisheng.A possible periodicity in the radio light curves of 3C120 ［J］.Astronomical Research ＆ Technology——Publications of National Astronomical Observatories of China，2007，4(4):307－312.

［11］劉凱博，楊東末，侯德東，等.Jurkevic方法在3C120天體光變周期分析中的應(yīng)用 ［J］.云南師范大學(xué)學(xué)報(bào)，2009，29(5):1－16.Liu Kaibo，Yang Dongmo，Hou Dedong，et al.Computing the period of light variability in seyfert galaxy 3C120 using the jurkevich method ［J］.Journal of Yunnan Normal University，2009，29(5):1－16.

［12］徐云冰，張雄，李孝攀，等.BL Lac天體ON 231的光變周期研究 ［J］.天文研究與技術(shù)——國(guó)家天文臺(tái)臺(tái)刊，2010，7(4):283－289.Xu Yunbing，Zhang Xiong，Li Xiaopan，et al.A study of the periodicity of the light curve of the BL lac object ON 231 ［J］.Astronomical Research ＆ Technology——Publications of National Astronomical Observatories of China，2010，7(4):283－289.

［13］Kidger M，Takalo L，Sillanpaa A.A new analysis of the 11-year period in OJ287-confirmation of its existence ［J］.Astronomy and Astrophysics，1992，264(1):32－36.

［14］Fan J H，Xie G Z，Wen S L.The relativistic beaming model for active galactic nuclei［J］.Astronomy and Astrophysics Supplement Series，1996，116(3):409－415.