射電脈沖星PSRB0329+54模式持續(xù)時標分布的觀測研究

陳建玲，王洪光，王 娜，蘇江波

(1. 運城學(xué)院 物理與電子工程系，山西 運城 044000；2. 中國科學(xué)院 新疆天文臺，新疆 烏魯木齊 830011)(3. 廣州大學(xué) 天體物理中心，廣東 廣州 510006)

射電脈沖星PSRB0329+54模式持續(xù)時標分布的觀測研究

陳建玲1,2，王洪光3，王 娜2，蘇江波1

(1. 運城學(xué)院 物理與電子工程系，山西 運城 044000；2. 中國科學(xué)院 新疆天文臺，新疆 烏魯木齊 830011)(3. 廣州大學(xué) 天體物理中心，廣東 廣州 510006)

利用烏魯木齊25米望遠鏡在2004年3月份將近20天的數(shù)據(jù)，我們研究了PSRB0329+54在1.5GHz的模式變化。通過測量得出每個子積分輪廓前導(dǎo)成分和后隨成分的相對強度的比值R的分布。結(jié)果發(fā)現(xiàn)，R的分布很好的符合高斯分布。我們統(tǒng)計了在這段觀測中正常模式和反常模式的內(nèi)稟持續(xù)時標分布，結(jié)果發(fā)現(xiàn)，正常模式和反常模式的時標分布均服從伽馬分布，最佳時標參數(shù)分別為154.2分鐘和31.5分鐘。

脈沖星；模式變化；持續(xù)時標；伽馬分布

0. 引言

脈沖星是能夠發(fā)出周期性信號的致密天體，大多數(shù)射電脈沖星的累積輪廓的形狀非常穩(wěn)定，但少數(shù)一些脈沖星的輪廓可以觀測到兩種或者兩種以上的形狀，這就是脈沖星的模式變化現(xiàn)象。PSRB0329+54是具有模式變化的典型的一顆源，這顆源流量很強，且模式變化頻繁。Lyne(1971)[1]第一次在408MHz探測到這顆源的反常模式，在每次幾個小時的連續(xù)觀測中，反常模式可以持續(xù)長達幾十分鐘。而在2.7GHz和14.8GHz頻率上的觀測說明，較高頻率反常模式的輪廓與其在408MHz的輪廓并不相同(Hesse1973，Bartel1978)[2，3]。1977年5月14.8GHz的觀測中，正常模式和反常模式分別觀測到的總時間為2.6hr和0.7hr(Bartel1978)[3]。

Bartel(1982)綜合了在410M-14.8GHz很寬頻率范圍內(nèi)的觀測結(jié)果，總結(jié)出脈沖星模式變化發(fā)生時，由于不同的脈沖成分相位發(fā)生變化而導(dǎo)致總的脈沖寬度變窄，反常模式不同成分的譜會變平或者變陡[4]。所有的偏振性質(zhì)，尤其是線偏振位置角都要受到模式變化的影響，并且發(fā)現(xiàn)模式變化在1.4和9.0GHz兩個頻率上同時發(fā)生。在對PSRB0329+54在1.4GHz的脈沖輪廓進行了長期監(jiān)測后(約2.3天)，他們得出結(jié)論，反常模式的持續(xù)時間一般從幾個脈沖周期到幾個小時，而模式變化的切換時間甚至可以短至一個周期。該項工作首次把1.4GHz的反常模式的輪廓分成三種。Suleimanova&Pugachev(2002)將觀測延伸到較低的頻率111.4MHz,并且發(fā)現(xiàn)，在模式變化發(fā)生時，只有一半的情況下外面成分相對于中心成分的相對強度也在同時變化[5]。

Xilourisetal.(1995)[6]首次對這顆源在10.5GHz進行了偏振觀測，發(fā)現(xiàn)在這個頻率上，反常模式出現(xiàn)的概率大約也為15%，這與1.4GHz的觀測結(jié)果是一致的，但這個頻率上的偏振性質(zhì)與1.7GHz的有很大差別(Barteletal.1978)。在32GHz頻率上，也可以很清楚地看到PSRB0329+54的模式變化(Krameretal. 1996)[7]。對這顆源模式變化的觀測最高頻率可以達到43GHz(Krameretal.1997)[8]。

利用烏魯木齊南山25米望遠鏡，Alietal.(2003)對PSRB0329+54在327MHz和610MHz兩個頻率的模式變化進行了觀測，分別得到了正常模式和反常模式所對應(yīng)的前導(dǎo)成分和后隨成分的相對譜指數(shù)[9]，之后，Liuetal. (2006)用同樣的系統(tǒng)對這顆源進行了單脈沖觀測，并且同時觀測到了正常模式和反常模式[10]。我們總結(jié)了以前關(guān)于這顆源的模式變化的工作，見表1。從中可以看出，以前的觀測工作主要限制在幾個小時到幾天的時間范圍內(nèi)，而觀測我們的工作研究了PSRB0329+54在1.5GHz2004年3月兩段八天的準連續(xù)數(shù)據(jù),這是目前同類數(shù)據(jù)中最長的,為我們深入分析模式時標特征提供了必備條件。

表1 PSR B0329+54的模式變化的觀測總結(jié)

1．觀測數(shù)據(jù)以及處理方法

我們的觀測數(shù)據(jù)利用的是南山25米望遠鏡得到的，其中心頻率為1540MHz，總帶寬為320MHz，一個周期采樣點 256 個bins，最短積分時間 59.9秒，更多關(guān)于觀測系統(tǒng)的詳細介紹見Wangetal.2001[11]。觀測時間為2004年3月12日到3月31日,由兩段準連續(xù)觀測組成， 每段持續(xù)8天(中間間斷 2.3 天)，前八天的持續(xù)時間為180.2小時，后八天的持續(xù)時間為178.2小時，總時間360 小時。觀測窗口有很多小的間隙，間隙的時間大約是幾分鐘，這可能是由于觀測程序原因造成的。在1.5GHz，PSRB0329+54輪廓不同模式的主要特征是成分Ⅰ和Ⅲ的相對強度的變化，如圖1所示，給出了這顆源的正常模式和反常模式的平均輪廓(該圖是對所有正常模式和反常模式持續(xù)時間分別長于25分鐘和15分鐘的輪廓序列進行疊加后得到的)。

圖1 PSR B0329+54正常模式和反常模式的平均輪廓圖

圖1中PSRB0329+54正常模式(實線)和反常模式(虛線)的平均輪廓，分別是對2000個以上的脈沖周期疊加得到的結(jié)果。圖中的點表示觀測的256個采樣點在脈沖輻射窗口和中間以及附近的部分。點線為劃分輪廓成分Ⅰ，Ⅱ，和Ⅲ時的邊界。

2. 結(jié)果

2.1 比例R隨時間的變化

圖2中給出了為2004年3月兩段數(shù)據(jù)前八天的連續(xù)觀測所對應(yīng)的R隨時間的變化。橫坐標表示觀測時間，縱坐標起伏較大的線表示比值R，起伏較小的線表示比值R′，圖中圓圈的部分表示觀測間隙。

從圖中可以很明顯地看出，比值R在(0.35-0.8)之間和(0.7-1.8)兩個區(qū)域之間變化，分別對應(yīng)著輪廓的正常模式和反常模式，比值R′和R也有相同的變化趨勢。根據(jù)比例R的明顯區(qū)別，我們可以很容易的區(qū)分出PSR B0329+54的脈沖輪廓的正常模式和反常模式。在兩段長時間的準連續(xù)觀測中，我們可以分別識別出244次正常模式和133次反常模式序列。

圖2

2.2 比例R的分布

圖3

圖3給出了兩段準連續(xù)數(shù)據(jù)的正常模式和反常模式比例R的分布。左圖為正常模式，右圖為反常模式。

我們用了2004年3月12日到3月31日的準連續(xù)數(shù)據(jù)來研究比例R的分布，這段時期內(nèi)總的積分時間為360.2小時，占了所有總觀測時間的69%。正常模式和反常模式比例R的直方圖分別見圖3，該圖表明R的分布很好地符合高斯分布，所以可以得出脈沖輪廓存在兩種主要模式的結(jié)論。反常模式的比例的分布要比正常模式的比例分布離散很多，這表明反常模式?jīng)]有正常模式穩(wěn)定，這與以前的工作是一致的，即反常模式比正常模式有更長的穩(wěn)定時標(Helfand1975)[13]。

4.3 模式持續(xù)時標分布

兩段觀測數(shù)據(jù)的時間窗口由148個觀測子窗口和47個時間間隙組成。圖4給出了準連續(xù)數(shù)據(jù)觀測到的正常模式和反常模式的持續(xù)時標分布。左圖為正常模式，右圖為反常模式，橫坐標為時間，單位為分鐘，縱坐標為出現(xiàn)的次數(shù)。由于觀測窗口時間為3個小時，所以正常模式的持續(xù)時間在175分鐘附近出現(xiàn)一個峰值。從中可以看出，反常模式比正常模式典型的時間尺度要短得多，且無論對于正常模式還是反常模式，在較長的時間尺度出現(xiàn)的次數(shù)較少。

圖4

由于觀測窗口存在間隙，所以我們要從觀測到的脈沖星模式持續(xù)時標的分布得到其內(nèi)稟的持續(xù)時標分布，需要用貝葉斯推斷的方法。Chenetal.2011年在工作中[14]，利用這兩段持續(xù)八天的準連續(xù)數(shù)據(jù)，并結(jié)合JodrellBank天文臺15米射電望遠鏡610MHz的每天觀測數(shù)據(jù)，利用貝葉斯推斷的方法對反常模式內(nèi)稟的持續(xù)時標做了限定，結(jié)果發(fā)現(xiàn)，相比正態(tài)分布，比對數(shù)正態(tài)分布以及帕累托分布，形狀參數(shù)略小于1的伽馬分布能夠更好地描述模式內(nèi)稟的持續(xù)時標分布。對于正常模式和反常模式，伽馬分布的最佳時標參數(shù)分別為154.2分鐘和31.5分鐘。而形狀參數(shù)也有類似的值，對于正常模式為0.75(-0.17～+0.22)，而對于反常模式則為0.84(-0.22～+0.28)，這說明兩種模式的物理機制可能是相同的。這項工作首次成功地限定了這顆脈沖星內(nèi)稟的持續(xù)時標分布，這為我們從物理上理解脈沖星模式變化現(xiàn)象提供了非常有用的信息。

3. 總結(jié)和討論

利用烏魯木齊南山25米射電望遠鏡在2004年3月將近20天的準連續(xù)觀測數(shù)據(jù)，我們研究了PSRB0329+54的模式變化現(xiàn)象。通過對比例R的分布的研究，結(jié)果發(fā)現(xiàn)PSRB0329+54有兩種主要模式，分別為正常模式和反常模式。兩種模式的主要區(qū)別在于前導(dǎo)成分和后隨成分相對強度的變化以及模式的持續(xù)時標。通過研究我們發(fā)現(xiàn)，正常模式和反常模式分別占了觀測時間的84.9%和15.1%，這與以前的工作是一致的。正常模式和反常模式觀測到的最長的持續(xù)時間分別達4.4小時和2.47小時。平均脈沖輪廓的正常模式和反常模式的前導(dǎo)和后隨成分的相對強度的比值都很好地符合高斯分布，其中正常模式和反常模式對應(yīng)的平均值分別為0.56和1.17。反常模式比正常模式對應(yīng)著更寬的高斯分布的寬度，這說明反常模式?jīng)]有正常模式穩(wěn)定。

[1]LyneA.G,Modechanginginpulsarradiation[J],MNRAS, 1971 ,153:27-32.

[2]HesseK.H.,PulseProfileVariationsinPSR0329 + 54at2695MHz[J],AstronomyandAstrophysics, 1973，27:373-377.

[3]BartelN.,SieberW.,WielebinskiR.,Observationsofpulsarsat14.8and22.7GHz[J],Astron.Astrophys, 1978,68:361-365.

[4]BartelN,MorrisD,SieberWetal,Themode-switchingphenomenoninpulsars[J],ApJ,1982 , 258:776-789.

[5]SuleimanovaS.A,PugachevV.D,BehaviorofPulsarB0329+54PulseCharacteristicsintheImmediateVicinityofModeSwitchTimesat111.4MHz[J],AstronomyReports, 2002, 46(l): 34-48.

[6]Xilourisetal.Pulsarpolarimetricobservationsat10.55GHz[J],AstronomyandAstrophysics, 1995, 293: 153-165.

[7]KramerM.,XilourisK.M.,JessnerA.etal.,Aturn-upinpulsarspectraatmm-wavelengths[J] ,AstronomyandAstrophysics, 1996, 306:867-876.

[8]KramerM.,JessnerA.,DoroshenkoO.etal.,ObservationsofPulsarsat7Millimeters[J],ApJ, 1997,488:364-367.

[9]Esamdin,Ali,Wu,X.J.,Zhang,X.Z.,Two-FrequencySimultaneousObservationsofMode-ChangingPhenomenaofPSR0329+54[J],ACTA,2003,44:223-226 .

[10]Liu,Z.Y.,Wang,N.,Urama,J.O.etal.,MonitoringofPulseIntensityandModeChangingforPSRB0329+54[J],ChJAA, 2006,Vol.6,Suppl.2, 64-67.

[11]WangN.,Manchester,R.N.,Zhang,J.etal.,PulsartimingatUrumqiAstronomicalObservatory:observingsystemandresults[J],MNRAS, 2001, 328:855-866.

[12]Rankin,J.M.,Towardanempiricaltheoryofpulsaremission.VI-Thegeometryoftheconalemissionregion:AppendixandtablesApJS[J], 1993, 85: 145-161.

[13]HelfandD.J,Manchester,R.N,TaylorJ.H,Observationsofpulsarradioemission.III-Stabilityofintegratedprofiles[J],ApJ,1975, 198:661-670.

[14]ChenJ.L.,Wang,H.G.,Wang,N.;etal.,LongtermmonitoringofmodeswitchingforPSRB0329+54[J],ApJ, 2011, 741, 48.

【責(zé)任編輯 馬太來】

The Observational Study of Timescale Distributions for PSR B0329+54 Mode

Chen Jian-Ling1，2，Wang Hong-Guang3, Wang-Na2, Su Jiang-Bo1

(1.DepartmentofPhysicsandElectronicEngineering，YunchengUniversity,Yuncheng044000,China) (2.XinjiangAstronomicalObservatory,Urumqi,Xinjiang,830011,China) (3.CenterforAstrophysics,GuangzhouUniversity,Guangzhou51006,China)

he mode-switching phenomenon of PSR B0329+54 was studied by using the 20 days’ data of March in 2004 by Urumqi 25 m radio telescope at 1540 MHz. The distributions of relative intensity between the leading and trailing component were measured for the individual profiles. It is found that the ratio follows Gauss distribution well. The intrinsic mode timescale distributions for normal and abnormal mode were also calculated with our data. The results show that both timescale distributions of normal and abnormal mode are according with Gamma distribution, where the best scales time parameters is 154.2 minutes for normal mode and 31.5 minutes for abnormal mode separately.

Pulsars；Mode change; Timescale distribution

2015-08-11

西部之光(XBBS201422)；國家自然科學(xué)基金(U1431125)；運城學(xué)院博士啟動基金(YQ-2012011,YQ-2014013)

陳建玲(1981-)，女，山西新絳人，運城學(xué)院物理與電子工程系講師，博士，研究方向為射電脈沖量的觀測與理論。

O

A

1008-8008(2015)06-0038-04