用周期图谱分析方法计算Blazar天体的光变周期

对适用于等间隔时间序列的周期图谱分析方法进行了研究,用模拟信号进行了检验,证实了它适用于有噪声的时间序列.并把这种方法应用到类星体3C 279、3C 345,BL Lac天体OJ 287、ON 231的光变周期分析中,得到它们的光变周期分别为7.14年、10.00年、11.76年、6.80年.结果表明用周期图谱分析方法得到的结论和其他文献用Jurkevich方法得到的结论一致.获得的周期对进一步研究Blazar天体的物理机制很有帮助.还分析了窗函数的影响,指出它们的优缺点,便于在实际工作中做出正确选择.     

Computing the Periods of Light Variation of Blazar Objects 3C279 and OJ 287 with Autoregressive Spectral Analysis Method

The light variability is one of the main characteristics of blazar objects. Because of the complexity of their light curves, the present periodicity analysis methods are not yet perfect. Based on the modern spectral estimate theory, this paper has described in details the principles of the maximum entropy spectral estimate and autoregressive (AR) spectral estimate, analyzed the effect of the order number selection on the resultant model. Applying these methods to the periodicity analysis of the quasar 3C 279 and BL Lac object OJ 287, their light periods are obtained to be 7.14 and 11.76 yr, respectively. As is verified by experiments, the AR spectral estimate has a high resolution and is a rather good periodicity analysis method. Finally, the items noteworthy for the application of these spectrum estimation methods to the periodicity analysis of the light variations of blazars are mentioned.     

An Analysis of Light Periods of BL Lac Object S5 0716+714 with the MUSIC Algorithm

The multiple signal classi?cation (MUSIC) algorithm is introduced to the estimation of light periods of BL Lac objects. The principle of the MUSIC algorithm is given, together with a testing on its spectral resolution by using a simulative signal. From a lot of literature, we have collected a large number of effective observational data of the BL Lac object S5 0716+714 in the three optical wavebands V, R, and I from 1994 to 2008. The light periods of S5 0716+714 are obtained by means of the MUSIC algorithm and average periodogram algorithm, respectively. It is found that there exist two major periodic components, one is the period of (3.33±0.08) yr, another is the period of (1.24±0.01) yr. The comparison of the performances of periodicity analysis of two algorithms indicates that the MUSIC algorithm has a smaller requirement on the sample length, as well as a good spectral resolution and anti-noise ability, to improve the accuracy of periodicity analysis in the case of short sample length.     

基于自回归方法的Blazar天体3C 279和OJ 287光变周期研究

探索Blazar光变资料中的周期或准周期变化是1个有待深入研究的领域.由于光变资料的复杂性,目前寻找周期的算法还不够完善.文中以现代谱估计为基础,详细论述了自回归(Auto-regressive,AR)模型谱估计方法和最大熵谱估计的基本原理,分析了阶数选择对模型的重要影响,并把这些方法应用到类星体3C 279和BL Lac天体OJ287的光变周期分析中,得到它们的光变周期分别为7.14和11.76yr.通过验证,自回归 模型谱估计方法由于其分辨率高,可以作为1种较好的分析周期的方法.最后指出在应用谱估计分析Blazar天体的光变周期注意事项.     

Spectral methods for analyzing the rotation of solar structures

The application of spectral analysis methods for studying the rotation of solar structures is considered. The time series characterizing the time variation of the solar He I 1083 nm emission in 5° latitude zones have been used. Three types of spectral analysis have been tested: the Welch method, the multitaper method, and the Schuster periodogram method. The first two methods have been chosen for the analysis of observing time intervals 26 and 3 years in length. The Schuster periodogram method is more suitable for the sliding spectral analysis in a 1-year-long temporal window with a shift by half a year. The chosen methods for analyzing the power spectra allow one to obtain the spectral densities, the powers of significant peaks in them and the corresponding periods, the total powers in the specified intervals of periods and to estimate the significance of the peaks found and the intervals in which the true periods corresponding to the peaks can be located.     

利用MUSIC算法分析BL Lac天体S5 0716+714光变周期

将基于多重信号分类的MUSIC谱估计算法引入BL Lac天体光变周期分析中.给出了MUSIC算法的基本原理,利用模拟信号检测了算法的频谱分辨率.从大量文献中收集了BL Lac天体S5 0716+714光学V、R、I 3个波段从1994年到2008年的有效观测数据,用MUSIC算法和平均周期图算法分别计算了它们的光变周期,发现存在两个主要光变周期:一个是(3.33±0.08)yr的周期,另一个是(1.24±0.01)yr的周期.对这两种算法的周期估计性能进行了比较,结果表明,MUSIC谱估计算法对样本长度要求较低,具有良好的分辨特性和抗噪声能力,能提高在样本长度较短情况下光变周期分析的准确性.     

Blazar天体3C 66A光学波段准周期光变分析

研究了Blazar天体3C 66A光学波段的准周期光变行为.收集了3C 66A光学V波段将近18 yr (2003—2021年)的测光数据,观测数据主要来源是:上海天文台(ShAO)、 AAVSO (The American Association of Variable Star Observers)数据库、Steward天文台.使用了Jurkevich和Lomb-Scargle两种方法分析了光变数据.Jurkevich方法得到了(850±90) d (～2.3 yr)和(1150±140) d (～3.2 yr)的光变周期,而Lomb-Scargle方法在充分考虑了“红噪声”效应之后同样得到了(869±70) d和(1111±90) d的光变周期,它们的置信水平分别为>99%和> 95%.通过与之前的研究结果比较,发现～2.3 yr的光变周期在3C 66A的历史光变数据中是一个稳定的周期,而～3.2 yr的周期则是不稳定的.     

Re-Analysis of QPO in 3C 273 Light Curve

We have developed analysis tools to search for quasi periodic oscillations in light curves from active galactic nuclei, using the following time series techniques: Wavelets, periodogram, Lomb–Scargle periodogram, structure function and multi-harmonic analysis of variance. The analysis tools incorporate different noise models with significant levels for all the techniques that is an improvement over the previous work. By looking for consistently high significance, we make the detection of periodicities more robust. We apply this tool to a previously reported QPO (Espaillat et al. 2008) in the X-ray light curve of 3C 273 with a periodicity of ~3300 s and find that the significance is only 74% in the wavelet and fails to show up above 95% significance in the periodogram and multi-harmonic analysis of variance.     

Significance of periodogram peaks and a pulsation mode analysis of the Beta Cephei star V403 Car

We discuss some commonly used methods for determining the significance of peaks in the periodograms of time series. We review methods for constructing the classical significance tests, their corresponding false alarm probability functions and the role played in these by independent random variables and by empirical and theoretical cumulative distribution functions. We discuss the concepts of independent frequencies and oversampling in periodogram analysis. We then compare the results of new Monte Carlo simulations for evenly spaced time series with results obtained previously by other authors, and present the results of Monte Carlo simulations for a specific unevenly spaced time series obtained for V403 Car.     

Long-term Variability Properties of Mkn501

Optical photometric observations at the I and R wavebands were carried out towards Mkn501 using the 1.56 meter optical telescope of Shanghai Astronomical Observatory at Sheshan. Combining our new observations with the published historical data, we have obtained the light curves of Mkn501 at the optical, infrared and radio wavebands with a time coverage of nearly 80 years. The relationship between the light variability and the color index is discussed, it is found that a strong correlation exists between the color indices (B − V) and (B − R), with the correlation coefficient reaching r = 0.73. The correlations of multi-band light variabilities are analyzed by the DCF method, it is found that certain positive correlations of B-band light variability with the 4.8 GHz and infrared light variabilities exist. And the spectral analysis on the B-band light curve with the CLEANest method indicates that the light curve of Mkn501 contains probably two periodical components of (10.06 ± 0.04) yr and (21.60 ± 0.17) yr.     

Comparison of irregular variations in the Earth’s rotation rate and parameters of geophysical and atmospheric processes with events on the sun

A study of the Earth’s rotation in space reveals a complex pattern of variations in its orientation, the excitation mechanisms of these variations, and their manifestations in various natural processes. The Earth’s rotation rate is not constant but exhibits complex fluctuations that account for some fraction of 10⁸ (corresponding to variations of several milliseconds (ms) in the length of the day). These variations span a wide spectrum of time scales, from hours to centuries or longer, reflecting the fact that they are produced by a wide variety of geophysical and astronomical processes. We discuss the results of our statistical comparison of long series of observations to reveal the most coherent variations. The spectral composition of the experimental time series has been determined using modified periodogram and single-channel autoregression methods. A comparative analysis has been performed by a two-channel autoregression spectral estimation method. The results of our comparison of the time series suggest that the fluctuations with periods of about 73 years are highly coherent.     

Search for orbital motion of the pulsar 4U 1626–67: Candidate for a neutron star with a supernova fall-back accretion disk

We report here results from a new search for orbital motion of the accretion powered X-ray pulsar 4U 1626–67 using two different analysis techniques. X-ray light curve obtained with the Proportional Counter Array of the Rossi X-ray Timing Explorer during a long observation carried out in February 1996, was used in this work. The spin period and the local period derivative were first determined from the broad 2–60 keV energy band light curve and these were used for all subsequent timing analysis. In the first technique, the orbital phase dependent pulse arrival times were determined for different trial orbital periods in the range of 500 to 10,000 s. We have determined a 3σ upper limit of 13 lt-ms on the projected semimajor axis of the orbit of the neutron star for most of the orbital period range, while in some narrow orbital period ranges, covering about 10% of the total orbital period range, it is 20lt-ms. In the second method, we have measured the pulse arrival times at intervals of 100 s over the entire duration of the observation. The pulse arrival time data were used to put an upper limit on any periodic arrival time delay using the Lomb-Scargle periodogram. We have obtained a similar upper limit of 10 lt-ms using the second method over the orbital period range of 500–10,000 s. This puts very stringent upper limits for the mass of the compact object except for the unlikely case of a complete face-on orientation of the binary system with respect to our line-of-sight. In the light of this measurement and the earlier reports, we discuss the possibility of this system being a neutron star with a supernovae fall-back accretion disk.     

Infrared variability properties of the blazar 3C 279

The long-term (about 27 years) near-infrared K light curve is constructed from the published literature for the blazar 3C 279. The Jurkevich method is adopted to analyse the periodicity, and a strong 7.1±0.44 yr period is found, suggesting that the next near-infrared outburst will occur in 2002/03. The correlation between colour index (spectral index) and magnitude is discussed, and a significant correlation between ( J − H ) and K is found with a correlation coefficient r =0.72 ( p =2.0×10⁻¹⁰), which is consistent with Brown et al.'s proposal.     

Wavelet Analysis of the New Solar Neutrino Capture Rate Data for the Homestake Experiment

Wavelet analysis offers an alternative to Fourier based time- series analysis and is particularly useful when spectral features are time dependent. We analyze the solar neutrino capture rate of the radiochemical Homestake chlorine experiment with abbreviated Morlet wavelets, using Foster's (AJ, 111,1709(1996)) rescaled wavelet technique. We emphasize the complementarity of wavelet analysis to Fourier analysis. Wavelet analysis confirms the results of previously undertaken Fourier analysis. The Homestake data seem to contain a harmonic content with periodicities of 4.76 yr, 1.89 yr, 0.85 yr, and 0.51 yr. Wavelet analysis reveals that the 4.76 yr and 1.89 yr periods show an almost constant behavior over the 25 yr Homestake data record, while the 0.85 yr and 0.51 yr periods exhibit a transient phenomenon. The analysis does not show strong evidence for a period of the solar 11 yr cycle. This revised version was published online in July 2006 with corrections to the Cover Date.     

Long‐term photometry of three active red giants in close binary systems: V2253 Oph,IT Com and IS Vir

We present and analyze long‐term optical photometric measurements of the three active stars V2253 Oph, IT Com and IS Vir. All three systems are single‐lined spectroscopic binaries with an early K giant as primary component but in different stages of orbital‐rotational synchronization. Our photometry is supplemented by 2MASS and WISE near‐IR and mid‐IR magnitudes and then used to obtain more accurate effective temperatures and extinctions. For V2253 Oph and IT Com, we found their spectral energy distributions consistent with pure photospheric emission. For IS Vir, we detect a marginal mid‐IR excess which hints towards a dust disk. The orbital and rotational planes of IT Com appear tobe coplanar, contrary to previous findings in the literature. We apply a multiple frequency analysis technique to determine photometric periods, and possibly changes of periods, ranging from days to decades. New rotational periods of 21.55±0.03 d, 65.1±0.3 d, and 23.50±0.04 d were determined for V2253 Oph, IT Com, and IS Vir, respectively. Splitting of these periods led to tentative detections of differential surface rotations of δP/P ≈ 0.02 for V2253 Oph and 0.07 for IT Com. Using a time‐frequency technique based on short‐term Fourier transforms we present evidence of cyclic light variations of length ≈ 10 yr for V2253 Oph and 5–6 yr for IS Vir. A single flip‐flop event has been observed for IT Com of duration 2–3 yr. Its exchange of the dominant active longitude had happened close to a time of periastron passage, suggesting some response of the magnetic activity from the orbital dynamics. The 21.55‐d rotational modulation of V2253 Oph showed phase coherence also with the orbital period, which is 15 times longer than the rotational period, thus also indicating a tidal feedback with the stellar magnetic activity. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Detecting non-sinusoidal periodicities in observational data using multiharmonic periodograms

We address the problem of assessing the statistical significance of candidate periodicities found using the so-called 'multiharmonic' periodogram, which is being used for detection of non-sinusoidal signals and is based on the least-squares fitting of truncated Fourier series. The recent author's investigation made for the Lomb–Scargle periodogram is extended to the more general multiharmonic periodogram. As a result, closed and efficient analytic approximations to the false alarm probability, associated with multiharmonic periodogram peaks, are obtained. The resulting analytic approximations are tested under various conditions using Monte Carlo simulations. The simulations showed a nice precision and robustness of these approximations.     

New spectroscopic and polarimetric observations of the A0 supergiant HD 92207

Our recent search for the presence of a magnetic field in the bright early A‐type supergiant HD 92207 using FORS 2 in spectropolarimetric mode revealed the presence of a longitudinal magnetic field of the order of a few hundred Gauss. However, the definite confirmation of the magnetic nature of this object remained pending due to the detection of shortterm spectral variability probably affecting the position of line profiles in left‐ and right‐hand polarized spectra. We present new magnetic field measurements of HD 92207 obtained on three different epochs in 2013 and 2014 using FORS 2 in spectropolarimetric mode. A 3σ detection of the mean longitudinal magnetic field using the entire spectrum, 〈B_z〉_all = 104 ± 34 G, was achieved in observations obtained in 2014 January. At this epoch, the position of the spectral lines appeared stable. Our analysis of spectral line shapes recorded in opposite circularly polarized light, i.e. in light with opposite sense of rotation, reveals that line profiles in the light polarized in a certain direction appear slightly split. The mechanism causing such a behaviour in the circularly polarized light is currently unknown. Trying to settle the issue of short‐term variability, we searched for changes in the spectral line profiles on a time scale of 8–10 min using HARPS polarimetric spectra and on a time scale of 3–4 min using time series obtained with the CORALIE spectrograph. No significant variability was detected on these time scales during the epochs studied. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

大行星轨道运动与太阳黑子数的中长周期变化

本文对不同序列的太阳黑子数资料作了分析研究,计算得到了可能的太阳黑子活动的中长周期变化,并分别与由大行星轨道运动引起的日心轨道角动量变化的周期进行比较,发现二者具有比较一致的谱结构。基于本文的讨论和文[17]的结论,我们进一步认为大行星轨道运动是太阳黑子数周期性变化的可能的外部因素。     

The Periodicity Analysis of Radio Flux Variability of PKS 1510-089 with the Period04 Method

We have collected the 37 GHz radio data of PKS 1510-089 from 1990 to 2005, and obtained its long-term light curve after making data processing. From the light curve we can find that the activity in PKS 1510-089 is very strong. In this paper, the Period04 method is used for the first time to analyze the light periodicity of PKS 1510-089. The result indicates that its radio emission at 37 GHz has the periods of (1.87±0.13) yr and (0.87±0.07) yr. This result is consistent with the results obtained by Xie et al. in 2004, 2005, 2008, and Wu et al. in 2005 with other methods of data analysis.     

A period study and light-curve synthesis for the Algol-type semidetached binary XX Cephei

We obtained CCD photometric observations of the Algol-type semidetached binary XX Cephei (XX Cep) during 15 nights from 2002 September 17 to 2003 February 2, and also on 2005 January 21. Except for those data taken on the last night of the concentrated observing season, the 3881 measurements were obtained over an interval of only 106 nights. From these data, four new times of minimum light were calculated. The  (O− C)  diagram formed from all available timings, and thus the orbital period of the system, can be partly represented as a beat effect between two cyclical variations with different periods (      yr,      yr) and amplitudes  ( K ₁=0.015 d, K ₂=0.103 d)  , respectively. Both physical and non-physical interpretations of these cycles were investigated. The long-term sinusoidal variation is too long for magnetic cycling in solar-type single and close binary stars. In addition, we have studied the effect of a possible secular period variation. By analysing the residuals from our Wilson–Devinney (WD) binary model, we found small light variations with a period of 5.99 d with amplitudes growing toward longer wavelengths. We think that these oscillations may be produced by instabilities at the systemic L ₁ point (also occupied by the point of the cool star) and that these instabilities are, in turn, caused by non-uniform and sporadic convection. There is also a short-period oscillation of about 45 min in the WD light residuals that is attributed to accretion on to the mass-gaining primary component from a feeble gas stream originating on the cool donor star.