基于EMD-AR谱的3C 273光变周期分析

光变是Blazar最显著的观测特征之一,有的光变呈准周期性. Blazar周期光变可能为其中心黑洞和吸积盘等内部物理结构提供重要信息.因此,周期分析对于讨论Blazar的光变有重要意义.结合经验模态分解(Empirical Mode Decomposition, EMD)和自回归(Auto-regressive, AR)模型谱估计方法的优点,提出了一种基于EMD-AR谱的组合光变周期分析新方法.先对观测数据进行EMD分解,得到各阶模态分量,并计算其与原始光变曲线的相关系数,再取相关度高的分量求和后估算AR谱,最后将此功率谱的周期分析方法用于天体的光变数据分析.讨论了EMD-AR谱方法的原理和应用步骤,并用此方法对类星体3C 273在1887—2016年的观测数据进行了分析,得到3C 273的长时标周期分别为:21.23、13.51、11.02、5.51、4.69、3.79、2.76 yr,与文献中常用方法的分析结果有很好的一致性;短时标周期分别为:(30±1)、(15±0.3)、(7.5±0.2)、(10±0.1)、(5±0.6)、(6±0.4)和(3±0.5) min,此结果未见相关...     

Time Lags between the 22 and 37 GHz Bursts of 48 Radio-loud AGNs

Based on the light curves at 22 and 37 GHz from the Metsahovi monitoring pro-gram, we investigate the time lags between the two radio bands for 48 radio-loud AGNs. DCF and ZDCF analyses are applied to the data. Our results show that there is a strong correlation between the two radio frequencies for all the sources, with the variations in the light curvesat 37 GHz leading the ones at 22 GHz in general. There is no obvious differences between different sub-class AGNs as regards the time lag. In two sources, it was found that the bursts at the lower frequency lead the ones at the higher frequency. One possible explanation is that electron acceleration dominates the light curve until the radiation reaches the maximum. Some sources, such as 3C 273, 3C 279, 3C 345 and 3C 454.3, have good enough data, so we can calculate their lags burst-by-burst. Our calculations show that different outbursts have dif- ferent lags. Some bursts have positive lags, most of bursts have no clear lags, and a few have negative lags. This result means that different bursts are triggered by different mechanisms, and the interpretation for the result involves both an intrinsic and a geometric mechanism. The positive lags are well consistent with the shock model, and we use these lags to calculate the typical magnetic field strength of the radiating region.     

Soft X-ray variability of the bright quasar 3C273

We present the results from ROSAT observations of 3C273 in the soft X-ray band. The light variation of 3C273 was investigated for three different energy bands of soft, medium, and hard. The maximum variability with a factor of 2 for 551 days was confirmed at all three different bands. This appears to be a periodic variation within the period of roughly 6 months. However, the short-term or micro variation was not so distinct and the light variation of each band did not show any correlation between them. The hardness ratio for hard and soft bands shows irregular variation but there was no correlation between them. There is no distinct variation of the photon index in the case of simple power law model fitting. For power law + free absorption model fitting, the average photon index (Γ) is 2.08.     

Constraints on a soft X-ray excess in the quasar 3C 279

We present the results of a three-week daily monitoring campaign of the quasar 3C 279 by the X-ray satellites RXTE ROSAT . A cross-correlation provides no evidence for any time lag between the very similar soft and hard X-ray light curves, and the source shows no significant spectral variability over the observing period. There is no evidence to support the presence of a soft excess, with a 99 per cent upper limit on any such component of 25 per cent of the total observed luminosity in the 0.1–2 keV band (<3 × 10³⁸ W). This fraction (but not the luminosity) is significantly lower than that of the soft excess observed in 3C 273.     

Simultaneous X-ray and IR variability in the quasar 3C 273

From a combination of high-quality X-ray observations from the NASA Rossi X-ray Timing Explorer ( RXTE ) and IR observations from the UK Infrared Telescope (UKIRT) we show that the medium-energy X-ray (3–20 keV) and near-IR fluxes in the quasar 3C 273 are highly correlated. It is widely believed that the X-ray emission in quasars like 3C 273 arises from Compton scattering of low-energy seed photons, and our observations provide the first reliable detection of correlated variations in 3C 273 between the X-ray band and any lower energy band. For a realistic electron distribution we demonstrate that it is probable that each decade of the seed-photon distribution from the mm to IR wavebands contributes roughly equally to the medium-energy X-ray flux. However, the expected mm variations are too small to be detected above the noise, probably explaining the lack of success of previous searches for a correlation between X-ray and mm variations. In addition, we show that the IR leads the X-rays by  0.75±0.25 d  . These observations rule out the 'External Compton' emission process for the production of the X-rays.     

The infrared variability of QSO 3C 273

The near-infrared data of the quasar 3C273 are collected from available literature and are used to analyze its variability properties. The largest amplitudes of variations are ΔJ = 1.0, ΔH = 0.9, and ΔK = 0.86. The color indices are J − H = 0.82, J − K = 1.96, and H − K = 1.13. Analysis with the discrete correlation function (DCF) method indicates that the variability between any two infrared bands is correlated without any time lag. The relation between the color index and brightness is also investigated and the result indicates that the spectrum steepens when the source first dims, while it flattens after the source has dimmed to a some extent, suggesting that the emission of 3C 273 at near-IR consists of at least two components. Some discussions are presented.     

3C 345和3C 273的光变分析

用结构函数法,对两个经典的blazar 3C 273和3C 345的22 GHz、37 GHz和光学波段的数据进行了周期分析.结果显示,3C 273的周期性比较明显,22 GHz经分析存在7.0年和14.7年的周期,37 GHz存在7.2年和14.5年的周期,同时它的光学波段可能存在0.7年和8.7年的周期.3C 345的37 GHz存在比较明显的8.8年的周期,22 GHz和光学波段分别存在不太明显的9.3年和10.2年的周期.     

K-band polarimetry of seven high-redshift radio galaxies

We present the results of K -band imaging polarimetry of seven 3CR radio galaxies with 0.7 <  z  < 1.2. We find strong evidence for polarization in three sources: 3C 22, 3C 41 and 3C 114. Of these, 3C 41 shows strong evidence of having a quasar core the infrared light of which is scattered by dust. We also find some evidence for polarization in 3C 54 and in 3C 356. The two point-like sources (3C 22 and 3C 41) and the barely-elongated 3C 54 appear to have of order 10 per cent of their K -band flux contributed by scattered light from the active nucleus. We conclude that scattered nuclear light can form a significant component of the near-infrared light emitted by high-redshift radio galaxies, and discuss models in which the scattering particles are electrons and dust grains.     

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.     

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的周期则是不稳定的.     

类星体3C 446射电光变非线性特性分析

类星体有剧烈、大幅度的光变现象, 光变研究有助于建立与观测相符的理论模型. 这篇文章从密歇根大学射电天文台数据库收集了类星体3C 446射电4.8、8.0和14.5GHz波段的长期观测数据. 传统的线性方法难以分析复杂的光变现象, 文章采用了集合经验模态分解(Ensemble Empirical Mode Decomposition, EEMD)方法和非线性分析方法相结合, 从混沌动力学特性、分形特性和周期性多角度对类星体光变随时间演化的规律进行了较全面的分析, 并重点对比分析了除去周期成分或混沌成分前后, 光变的周期性和非线性特性是否存在明显区别. 分析结果表明, 类星体3C 446射电波段光变资料由周期成分、趋势成分和混沌成分组成, 光变具有周期性、混沌性和分形特性. 除去混沌成分和趋势成分后的光变周期与原始光变资料的周期完全相同, 而两者的混沌和分形特性有明显不同. 从饱和关联维数来看, 重构动力学系统时, 除去周期成分和趋势成分后的光变资料比原始光变资料需要更多的独立参量, Kolmogorov熵值表明前者信息的损失率比后者大, 系统的混沌程度更高, 系统也更复杂, Hurst值表明后者自相似性和长程相关性比前者略强.     

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.     

Radio and X-ray observations of an exceptional radio flare in the extreme z= 4.72 blazar GB B1428+4217

We report on the extreme behaviour of the high-redshift blazar GB B1428+4217 at   z = 4.72  . A continued programme of radio measurements has revealed an exceptional flare in the light curve, with the 15.2-GHz flux density rising by a factor of ∼3 from ∼140 to ∼430  mJy in a rest-frame time-scale of only ∼4 months – much larger than any previous flares observed in this source. In addition to new measurements of the 1.4–43  GHz radio spectrum, we also present the analysis and results of a target-of-opportunity X-ray observation using XMM–Newton , made close to the peak in radio flux. Although the X-ray data do not show a flare in the high-energy light curve, we are able to confirm the X-ray spectral variability hinted at in previous observations. GB B1428+4217 is one of several high-redshift radio-loud quasars that display a low-energy break in the X-ray spectrum, probably due to the presence of excess absorption in the source. X-ray spectral analysis of the latest XMM–Newton data is shown to be consistent with the warm-absorption scenario which we have hypothesized previously. Warm absorption is also consistent with the observed X-ray spectral variability of the source, in which the spectral changes can be successfully accounted-for with a fixed column density of material in which the ionization state is correlated with hardness of the underlying power-law emission.     

High resolution imaging of parsec-scale jets: the impact of the VSOP observations

Resolution improvements achieved with the ongoing space VLBI mission VSOP provide an excellent opportunity for studying compact, parsec-scale jets in active galactic nuclei. Visibilities on space baselines contain structural information that cannot be recovered by super-resolving ground-array data at the same frequency. We illustrate this by discussing the results of recent VSOP observations of parsec-scale jets in 0836+710 and 3C273.     

Orbital structure in N-body models of barred-spiral galaxies

We study the orbital structure in a series of self-consistent N -body configurations simulating rotating barred galaxies with spiral and ring structures. We perform frequency analysis in order to measure the angular and the radial frequencies of the orbits at two different time snapshots during the evolution of each N -body system. The analysis is done separately for the regular and the chaotic orbits. We thereby identify the various types of orbits, determine the shape and percentages of the orbits supporting the bar and the ring/spiral structures, and study how the latter quantities change during the secular evolution of each system. Although the frequency maps of the chaotic orbits are scattered, we can still identify concentrations around resonances. We give the distributions of frequencies of the most important populations of orbits. We explore the phase-space structure of each system using projections of the 4D surfaces of section. These are obtained via the numerical integration not only of the orbits of test particles, but also of the real N -body particles. We thus identify which domains of the phase space are preferred and which are avoided by the real particles. The chaotic orbits are found to play a major role in supporting the shape of the outer envelope of the bar as well as the rings and the spiral arms formed outside corotation.     

Predicting the second caustic crossing in binary microlensing events

We fit binary lens models to the data covering the initial part of real microlensing events in an attempt to predict the time of the second caustic crossing. We use approximations during the initial search through the parameter space for light curves that roughly match the observed ones. Exact methods for calculating the lens magnification of an extended source are used when we refine our best initial models. Our calculations show that the reliable prediction of the second crossing can only be made very late, when the light curve has risen appreciably after the minimum between the two caustic crossings. The best observational strategy is therefore to sample as frequently as possible once the light curve starts to rise after the minimum.     

Stellar disruption by a supermassive black hole: is the light curve really proportional to t−5/3?

In this paper, we revisit the arguments for the basis of the time evolution of the flares expected to arise when a star is disrupted by a supermassive black hole. We present a simple analytic model relating the light curve to the internal density structure of the star. We thus show that the standard light curve proportional to   t ^−5/3  only holds at late times. Close to the peak luminosity the light curve is shallower, deviating more strongly from   t ^−5/3  for more centrally concentrated (e.g. solar type) stars. We test our model numerically by simulating the tidal disruption of several stellar models, described by simple polytropic spheres with index γ. The simulations agree with the analytical model given two considerations. First, the stars are somewhat inflated on reaching pericentre because of the effective reduction of gravity in the tidal field of the black hole. This is well described by a homologous expansion by a factor which becomes smaller as the polytropic index becomes larger. Secondly, for large polytropic indices wings appear in the tails of the energy distribution, indicating that some material is pushed further away from parabolic orbits by shocks in the tidal tails. In all our simulations, the   t ^−5/3  light curve is achieved only at late stages. In particular, we predict that for solar-type stars, this happens only after the luminosity has dropped by at least 2 mag from the peak. We discuss our results in the light of recent observations of flares in otherwise quiescent galaxies and note the dependence of these results on further parameters, such as the star/hole mass ratio and the stellar orbit.     

Long-term optical and X-ray variability of the Seyfert galaxy Markarian 79

We present the results of concurrent X-ray and optical monitoring of the Seyfert 1 galaxy Mrk 79 over a period of more than 5 yr. We find that on short to medium time-scales (days to a few tens of days) the 2–10 keV X-ray and optical u - and V -band fluxes are significantly correlated, with a delay between the bands consistent with 0 d. We show that most of these variations may be well reproduced by a model where the short-term optical variations originate from reprocessing of X-rays by an optically thick accretion disc. The optical light curves, however, also display long time-scale variations over thousands of days, which are not present in the X-ray light curve. These optical variations must originate from an independent variability mechanism and we show that they can be produced by variations in the (geometrically) thin disc accretion rate as well as by varying reprocessed fractions through changes in the location of the X-ray corona.     

The long-term multiband optical observations and colour index for the quasar 3C 273

The long baseline simultaneous multiband ( BVRI ) observations of blazar 3C 273 are presented. We have made 758 optical multiband observations with the Yunnan Astronomical Observatory (YAO) 103 cm and Shanghai Astronomical Observatory (SHAO) 156 cm from 2003 January to 2005 April. In our observational campaign, the average colour indices are   B − V = 0.21, V − R = 0.26  and   R − I = 0.31  . The mean magnitudes in B, V, R and I are 12.895, 12.698, 12.441 and 12.139, respectively. The source is in a steady state over observational campaign. The correlations between colour index and brightness are discussed. We find that the correlations of   B − V   versus B ,   V − R   versus V and   R − I   versus R have significant negative correlation, while   B − V   versus V ,   V − R   versus R and   R − I   versus I have positive correlation. These strong correlations imply that the spectrum becomes bluer (flatter) when the source becomes brighter, and redden (softer) when the source fades both in intraday and long-term variability. The spectral evolution trends of 3C 273 are consistent with those of BL Lac objects.     

Submillimetre detection of a high-redshift type 2 QSO

We report on the first SCUBA detection of a type 2 QSO at   z = 3.660  in the Chandra Deep Field South. This source is X-ray-absorbed, shows only narrow emission lines in the optical spectrum and is detected in the submillimetre: it is the ideal candidate in an evolution scheme for active galactic nuclei (AGN) (e.g. Fabian 1999 ; Page et al. 2004 ) of an early phase corresponding to the main growth of the host galaxy and formation of the central black hole. The overall photometry (from the radio to the X-ray energy band) of this source is well reproduced by the spectral energy distribution (SED) of NGC 6240, while it is incompatible with the spectrum of a type 1 QSO (3C 273) or a starburst galaxy (Arp 220). Its submillimetre (850 μm) to X-ray (2 keV) spectral slope  (α_SX)  is close to the predicted value for a Compton-thick AGN in which only 1 per cent of the nuclear emission emerges through scattering. Using the observed flux at 850 μm we have derived a star formation rate of  550–680 M_⊙ yr⁻¹  and an estimate of the dust mass   M _dust= 4.2 × 10⁸ M_⊙  .