Posts Tagged time lag

Today's Postings

Joint spectral-timing modelling of the hard lags in GX 339-4: constraints on reflection models

The X-ray variations of hard state black hole X-ray binaries above 2 keV show ‘hard lags’, in that the variations at harder energies follow variations at softer energies, with a time-lag \tau depending on frequency \nu approximately as \tau \propto \nu^{-0.7}. Several models have so far been proposed to explain this time delay, including fluctuations propagating through an accretion flow, spectral variations during coronal flares, Comptonisation in the extended hot corona or a jet, or time-delays due to large-scale reflection from the accretion disc. In principle these models can be used to predict the shape of the energy spectrum as well as the frequency-dependence of the time-lags, through the construction of energy-dependent response functions which map the emission as a function of time-delay in the system. Here we use this approach to test a simple reflection model for the frequency-dependent lags seen in the hard state of GX 339-4, by simultaneously fitting the model to the frequency-dependent lags and energy spectrum measured by XMM-Newton in 2004 and 2009. Our model cannot simultaneously fit both the lag and spectral data, since the relatively large lags require an extremely flared disc which subtends a large solid angle to the continuum at large radii, in disagreement with the observed Fe K\alpha emission. Therefore, we consider it more likely that the lags > 2 keV are caused by propagation effects in the accretion flow, possibly related to the accretion disc fluctuations which have been observed previously.

Evidence for Widespread Cooling in an Active Region Observed with the SDO Atmospheric Imaging Assembly

A well known behavior of EUV light curves of discrete coronal loops is that the peak intensities of cooler channels or spectral lines are reached at progressively later times than hotter channels. This time lag is understood to be the result of hot coronal loop plasma cooling through these lower respective temperatures. However, loops typically comprise only a minority of the total emission in active regions. Is this cooling pattern a common property of active region coronal plasma, or does it only occur in unique circumstances, locations, and times? The new SDO/AIA data provide a wonderful opportunity to answer this question systematically for an entire active region. We measure the time lag between pairs of SDO/AIA EUV channels using 24 hours of images of AR 11082 observed on 19 June 2010. We find that there is a time-lag signal consistent with cooling plasma, just as is usually found for loops, throughout the active region including the diffuse emission between loops for the entire 24 hour duration. The pattern persists consistently for all channel pairs and choice of window length within the 24 hour time period, giving us confidence that the plasma is cooling from temperatures of greater than 3 MK, and sometimes exceeding 7 MK, down to temperatures lower than ~ 0.8 MK. This suggests that the bulk of the emitting coronal plasma in this active region is not steady; rather, it is dynamic and constantly evolving. These measurements provide crucial constraints on any model which seeks to describe coronal heating.

Continuum reverberation mapping in a z = 1.41 radio-loud quasar

Q0957+561 was the first discovered gravitationally lensed quasar. The mirage shows two images of a radio-loud quasar at redshift z = 1.41. The time lag between these two images is well established around one year. We detected a very prominent variation in the optical brightness of Q0957+561A at the beginning of 2009, which allowed us to predict the presence of significant intrinsic variations in multi-wavelength light curves of Q0957+561B over the first semester of 2010. To study the predicted brightness fluctuations of Q0957+561B, we conducted an X-ray, NUV, optical and NIR monitoring campaign using both ground-based and space-based facilities. The continuum NUV-optical light curves revealed evidence of a centrally irradiated, standard accretion disk. In this paper, we focus on the radial structure of the standard accretion disk and the nature of the central irradiating source in the distant radio-loud active galactic nucleus (AGN).

NGC454: unveiling a new "changing look" AGN

We present a detailed analysis of the X-ray spectrum of the Seyfert 2 galaxy NGC454E, belonging to the interacting system NGC454. Observations performed with Suzaku, XMM-Newton and Swift allowed us to detect a dramatic change in the curvature of the 2-10 keV spectrum, revealing a significant variation of the absorbing column density along the line of sight (from ~ 1 x10^{24}cm^{-2} to ~ 1×10^{23}cm^{-2}). Consequently, we propose this source as a new member of the class of “changing look” AGN, i.e. AGN that have been observed both in Compton-thin (NH =10^{23 cm^{-2}) and reflection dominated states (Compton-thick, NH >10^{24} cm^{-2}). Due to the quite long time lag (6 months) between the Suzaku and XMM-Newton observations we cannot infer the possible location of the obscuring material causing the observed variability. In the 6-7 keV range the XMM-Newton observation also shows a clear signature of the presence of an ionized absorber. Since this feature is not detected during the Suzaku observation (despite its detectability), the simplest interpretation is that the ionized absorber is also variable; its location is estimated to be within ~10^{-3} pc from the central black hole, probably much closer in than the rather neutral absorber.

Ensemble Variability of Near-Infrared-Selected Active Galactic Nuclei

We present the properties of the ensemble variability $V$ for nearly 5000 near-infrared (NIR) AGNs selected from the catalog of Quasars and Active Galactic Nuclei (13th Ed.) and the SDSS-DR7 quasar catalog. From 2MASS, DENIS, and UKIDSS/LAS point source catalogs, we extract 2MASS-DENIS and 2MASS-UKIDSS counterparts for cataloged AGNs by catalog cross-identification. We further select variable AGNs based on an optimal criterion for selecting the variable sources. The sample objects are divided into subsets according to whether NIR light originates by optical or NIR emission in the rest frame; and we examine the correlations of the ensemble variability with the rest-frame wavelength, redshift, luminosity, and rest-frame time lag. In addition, we also examine the correlations of variability amplitude with optical variability, radio intensity, and radio-to-optical flux ratio. The rest-frame optical variability of our samples shows known negative correlations with luminosity and positive correlations with rest-frame time lag (i.e., the structure function, SF). However, no well-known negative correlation exists between the rest wavelength and optical variability. This inconsistency might be due to a biased sampling of high-redshift AGNs. NIR variability in the rest frame is anticorrelated with the rest wavelength, which is consistent with previous suggestions. However, correlations of NIR variability with luminosity and rest-frame time lag are the opposite of these correlations of the optical variability; that is, the NIR variability is positively correlated with luminosity but negatively correlated with the rest-frame time lag. Because these trends are qualitatively consistent with the properties of radio-loud quasars reported by some previous studies, most of our sample objects are probably radio-loud quasars. Finally, we also discuss the negative correlations seen in the NIR SFs.

Discovery of the black hole mass and soft lag relation in active galactic nuclei

We carried out a systematic analysis of time lags between X-ray energy bands in a large sample (32 sources) of unabsorbed, radio quiet active galactic nuclei (AGN), observed by XMM-Newton. The analysis of X-ray lags is performed in the Fourier-frequency domain, between energy bands where the soft excess (soft band) and the primary power law (hard band) dominate the emission. We detected a soft/negative lag in a total of 15 out of 32 sources. Considering that 7 of these have not been previously reported in the literature, this work more than doubles the number of known sources with a soft/negative lag. The characteristic time-scales (i.e. frequency and amplitude) of the soft/negative lag do show a highly significant (i.e. $\gsim 4\sigma$) correlation with the black hole mass. The measured correlations indicate that soft lags are systematically shifted to lower frequencies and higher absolute amplitudes as the mass of the source increases. To first approximation, all the sources in the sample are consistent with having similar mass-scaled lag properties. These results clearly demonstrate the existence of a mass-scaling law for the soft/negative lag, that holds for AGN spanning a large range of masses (about 2.5 orders of magnitude), thus supporting the idea that soft lags originate in the innermost regions of AGN and are powerful tools for testing their physics and geometry.

Multiwavelength Variability of the Blazars Mrk 421 and 3C 454.3 in the High-State [Replacement]

We report the results of photometric observations of the blazars Mrk 421 and 3C 454.3 designed to search for intraday variability (IDV) and short-term variability (STV). Optical photometric observations were spread over eighteen nights for Mrk 421 and seven nights for 3C 454.3 during our observing run in 2009-2010 at the 1.04 m telescope at ARIES, India. Genuine IDV is found for the source 3C 454.3 but not for Mrk 421. Genuine STV is found for both sources. Mrk 421 was revealed by the MAXI X-ray detector on the International Space Station to be in an exceptionally high flux state in 2010 January – February. We performed a correlation between the X-ray and optical bands to search for time delays and found a weak correlation with higher frequencies leading the lower frequencies by about ten days. The blazar 3C 454.3 was found to be in high flux state in November-December 2009. We performed correlations in optical observations made at three telescopes, along with X-ray data from the MAXI satellite and public release gamma-ray data from the Fermi space telescope. We found strong correlations between the gamma-ray and optical bands at a time lag of about four days but the X-ray flux is not correlated with either. We briefly discuss the possible reasons for the time delays between these bands within the framework of existing models for X-ray and gamma-ray emission mechanisms.

Multiwavelength Variability of the Blazars Mrk 421 and 3C 454.3 in the High-State

We report the results of photometric observations of the blazars Mrk 421 and 3C 454.3 designed to search for intraday variability (IDV) and short-term variability (STV). Optical photometric observations were spread over eighteen nights for Mrk 421 and seven nights for 3C 454.3 during our observing run in 2009-2010 at the 1.04 m telescope at ARIES, India. Genuine IDV is found for the source 3C 454.3 but not for Mrk 421. Genuine STV is found for both sources. Mrk 421 was revealed by the MAXI X-ray detector on the International Space Station to be in an exceptionally high flux state in 2010 January – February. We performed a correlation between the X-ray and optical bands to search for time delays and found a weak correlation with higher frequencies leading the lower frequencies by about ten days. The blazar 3C 454.3 was found to be in high flux state in November-December 2009. We performed correlations in optical observations made at three telescopes, along with X-ray data from the MAXI satellite and public release gamma-ray data from the Fermi space telescope. We found strong correlations between the gamma-ray and optical bands at a time lag of about four days but the X-ray flux is not correlated with either. We briefly discuss the possible reasons for the time delays between these bands within the framework of existing models for X-ray and gamma-ray emission mechanisms.

A Description of Quasar Variability Measured Using Repeated SDSS and POSS Imaging

We provide a quantitative description and statistical interpretation of the optical continuum variability of quasars. The Sloan Digital Sky Survey (SDSS) has obtained repeated imaging in five UV-to-IR photometric bands for 34,727 spectroscopically confirmed quasars. About 10,000 quasars have an average of 60 observations in each band obtained over a decade along stripe 82 (S82), whereas the remaining ~25,000 have 2-3 observations due to scan overlaps. The observed time lags span the range from a day to almost 10 years, and constrain quasar variability at rest-frame time lags of up to 4 years, and at rest-frame wavelengths from 1000A to 6000A. We publicly release a user-friendly catalog of quasars from the SDSS Data Release 7 that have been observed at least twice in SDSS or once in both SDSS and the Palomar Observatory Sky Survey, and we use it to analyze the ensemble properties of quasar variability. Based on a damped random walk (DRW) model defined by a characteristic time scale and an asymptotic variability amplitude that scale with the luminosity, black hole mass, and rest wavelength for individual quasars calibrated in S82, we can fully explain the ensemble variability statistics of the non-S82 quasars. The distribution of magnitude changes is exponential at large values and constrains the intrinsic scatter in DRW model parameters among quasars with fixed luminosity, black hole mass, and rest wavelength. All available data are consistent with the DRW model as a viable description of the optical continuum variability of quasars on time scales of ~1-2000 days in the rest frame. We use these models to predict the incidence of quasar contamination in transient surveys such as those from PTF and LSST.

Ensemble X-ray variability of Active Galactic Nuclei from serendipitous source catalogues [Replacement]

The X-ray variability of the Active Galactic Nuclei (AGN) has been most often investigated with studies of individual, nearby, sources, and only a few ensemble analyses have been applied to large samples in wide ranges of luminosity and redshift. We want to determine the ensemble variability properties of two serendipitously selected AGN samples extracted from the catalogues of XMM-Newton and Swift, with redshift between ~0.2 and ~4.5, and X-ray luminosities, in the 0.5-4.5 keV band, between ~10^43 erg/s and ~10^46 erg/s. We use the structure function (SF), which operates in the time domain, and allows for an ensemble analysis even when only a few observations are available for individual sources and the power spectral density (PSD) cannot be derived. SF is also more appropriate than fractional variability and excess variance, because such parameters are biased by the duration of the monitoring time interval in the rest-frame, and thus by cosmological time dilation. We find statistically consistent results for the two samples, with the SF described by a power law of the time lag, approximately as SF \propto tau^0.1. We do not find evidence of the break in the SF, at variance with the case of lower luminosity AGNs. We confirm a strong anti-correlation of the variability with X-ray luminosity, accompanied by a change of the slope of the SF. We find evidence in support of a weak, intrinsic, average increase of X-ray variability with redshift. The change of amplitude and slope of the SF with X-ray luminosity provides new constraints on both single oscillator models and multiple subunits models of variability.

Ensemble X-ray variability of Active Galactic Nuclei from serendipitous source catalogs

The X-ray variability of the Active Galactic Nuclei (AGN) has been most often investigated with studies of individual, nearby, sources, and only a few ensemble analyses have been applied to large samples in wide ranges of luminosity and redshift. We want to determine the ensemble variability properties of two serendipitously selected AGN samples extracted from the catalogs of XMM-Newton and Swift, with redshift between ~0.2 and ~4.5, and X-ray luminosities, in the 0.5-4.5 keV band, between ~10^43 erg/s and ~10^46 erg/s. We use the structure function (SF), which operates in the time domain, and allows for an ensemble analysis even when only a few observations are available for individual sources and the power spectral density (PSD) cannot be derived. SF is also more appropriate than fractional variability and excess variance, because such parameters are biased by the duration of the monitoring time interval in the rest-frame, and thus by cosmological time dilation. We find statistically consistent results for the two samples, with the SF described by a power law of the time lag, approximately as SF tau^0.1. We do not find evidence of the break in the SF, at variance with the case of lower luminosity AGNs. We confirm a strong anti-correlation of the variability with X-ray luminosity, accompanied by a change of the slope of the SF. We find evidence in support of a weak, intrinsic, average increase of X-ray variability with redshift. The change of amplitude and slope of the SF with X-ray luminosity provides new constraints on both single oscillator models and multiple subunits models of variability.

Oscillations in the Solar Faculae. III. The Phase Relations between Chromospheric and Photospheric LOS Velocities

An analysis of line-of-sight velocity oscillation in nine solar faculae was undertaken with the aim of studying of phase relations between chromosphere (He i 10830 A line) and photosphere (Si i 10827 A line) five-minute oscillations. We found that time lag of the chromospheric signal relative to photospheric one varies from -12 to 100 seconds and is about 50 seconds on average.We assume that the small observed lag can have three possible explanations: i) convergence of formation levels of He i 10830 A and Si i 10827 A in faculae; ii) significant increase of five-minute oscillation propagation velocity above faculae; iii) simultaneous presence of standing and travelling waves.

Constraining properties of dusty environments by infrared variability

We present model simulations of time-variable infrared (IR) emission from dust as a consequence of variability of the incident radiation. For that we introduce a generalized treatment for temperature variations in a dusty environment, which is not limited to any specific astronomical source. The treatment has been incorporated into a simplified clumpy torus model, with the radial brightness distribution as the main parameter, to study the IR emission of type 1 active galactic nuclei (AGN). We show that any variability signal in the optical is smoothened stronger if the brightness distribution is very extended, and this smoothing strongly depends on wavelength. This also affects time lags between the optical and near-/mid-IR emission, which can be up to 10s of sublimation radii for long wavelengths and extended brightness distributions. The dependence of time lag on wavelength and distribution can be used to quantify the brightness distribution in an AGN torus, either by comparing optical light curves to near-IR and mid-IR light curves, or by directly comparing near-IR to mid-IR light curves. Moreover, our model has been applied to near-IR data of the nearby Seyfert 1 galaxy NGC 4151. We show that the simple model can reproduce the overall observed variability signal and found that about 40% of the energy in the variability signal in the V-band has been converted into K-band variability. This low value may be explained by a “snowball” model of gradually-sublimating clouds at the inner edge of the torus. We also note that our modeling does not support a change of time lag/sublimation radius over the observed light curve epoch in spite of a significant change in V-band emission.

Photometric AGN reverberation mapping - an efficient tool for BLR sizes, black hole masses and host-subtracted AGN luminosities

Photometric reverberation mapping employs a wide bandpass to measure the AGN continuum variations and a suitable band, usually a narrow band (NB), to trace the echo of an emission line in the broad line region (BLR). The narrow band catches both the emission line and the underlying continuum, and one needs to extract the pure emission line light curve. We performed a test on two local AGNs, PG0003+199 (=Mrk335) and Ark120, observing well-sampled broad- (B, V) and narrow-band light curves with the robotic 15cm telescope VYSOS-6 on Cerro Armazones, Chile. In PG0003+199, H_alpha dominates the flux in the NB by 85%, allowing us to measure the time lag of H_alpha against B without the need to correct for the continuum contribution. In Ark120, H_beta contributes only 50% to the flux in the NB. The cross correlation of the B and NB light curves shows two distinct peaks of similar strength, one at lag zero from the autocorrelated continuum and one from the emission line at tau_cent = 47.5 +/- 3.4 days. We constructed a synthetic H_beta light curve, by subtracting a scaled V light curve, which traces the continuum, from the NB light curve. The cross correlation of this synthetic H_beta light curve with the B light curve shows only one major peak at tau_cent = 48.0 +/- 3.3 days, while the peak from the autocorrelated continuum at lag zero is absent. We conclude that, as long as the emission line contributes at least 50% to the bandpass, the pure emission line light curve can be reconstructed from photometric monitoring data so that the time lag can be measured. For both objects the lags we find are consistent with spectroscopic reverberation results. While the dense sampling (median 2 days) enables us to determine tau_cent with small (10%) formal errors, we caution that gaps in the light curves may lead to much larger systematic uncertainties. (Abstract shortened, see the manuscript.)

Reverberation Mapping of the Intermediate Mass Nuclear Black Hole in SDSS J114008.71+030711.4

We present the results of a reverberation mapping (RM) campaign on the black hole (BH) associated with the active galactic nucleus (AGN) in SDSS J114008.71+030711.4 (hereafter GH08). This object is selected from a sample of 19 candidate intermediate mass BHs (M_{BH} < 10^{6} Msun) found by Greene & Ho 2004 in the Sloan Digital Sky Survey (SDSS). We used the Hobby-Eberly Telescope to obtain 30 spectra over a period of 178 days in an attempt to resolve the reverberation time lag (tau) between the continuum source and the broad line region (BLR) in order to determine the radius of the BLR (R_{BLR}) in GH08. We measure tau to be 2 days with an upper limit of 6 days. We estimate the AGN luminosity at 5100 Angstroms to be approximately 1.1 x 10^{43} erg s^{-1} after deconvolution from the host galaxy. The most well calibrated R_{BLR}-L relation predicts a time lag which is 4 times larger than what we measure. Using the measured H\beta\ full-width-at-half-maximum of 703 (+/-) 110 km s^{-1} and an upper limit for R_{BLR} = 6 light days, we find M_{BH} < 5.8 x 10^{5} Msun as an upper limit to the BH virial mass in GH08, which implies super-Eddington accretion. Based on our measured M_{BH} we propose that GH08 may be another candidate to add to the very short list of AGNs with M_{BH} < 10^{6} Msun determined using RM.

Constraining broad-line regions from time lags of broad emission lines relative to radio emission

In this paper, a new method is proposed to estimate the broad-line region sizes of UV lines $R^{\rm{uv}}_{\rm{BLR}}$. It is applied to 3C 273. First, we derive the time lags of radio emission relative to broad emission lines Ly$\alpha$ and C IV by the ZDCF method. The broad lines lag the 5, 8, 15, 22 and 37 GHz emission. The measured lags $\tau^{\rm{uv}}_{\rm{ob}}$ are of the order of years. For a given line, $\tau^{\rm{uv}}_{\rm{ob}}$ decreases as the radio frequency increases. This trend results from the radiative cooling of relativistic electrons. Both UV lines have a lag of $\tau^{\rm{uv}}_{\rm{ob}}=-2.74^{+0.06}_{-0.25}$ yr relative to the 37 GHz emission. These results are consistent with those derived from the Balmer lines in Paper I. Second, we derive the time lags of the lines Ly$\alpha$, CIV, H$\gamma$, H$\beta$ and H$\alpha$ relative to the 37 GHz emission by the FR/RSS Monte Carlo method. The measured lags are $\tau_{\rm{ob}}=-3.40^{+0.31}_{-0.05}$, $-3.40^{+0.41}_{-0.14}$, $-2.06^{+0.36}_{-0.92}$, $-3.40^{+1.15}_{-0.20}$ and $-3.56^{+0.35}_{-0.18}$ yr for the lines Ly$\alpha$, CIV, H$\gamma$, H$\beta$ and H$\alpha$, respectively. These estimated lags are consistent with those derived by the ZDCF method within the uncertainties. Based on the new method, we derive $R^{\rm{uv}}_{\rm{BLR}}=2.54^{+0.71}_{-0.35}$–$4.01^{+0.90}_{-1.16}$ and $2.54^{+0.80}_{-0.43}$–$4.01^{+0.98}_{-1.24}$ light-years for the Ly$\alpha$ and CIV lines, respectively. Considering the uncertainties, these estimated sizes are consistent with those obtained in the classical reverberation mapping for the UV lines and the Balmer lines. This indicates that their emitting regions are not separated so large as in the classical mapping of the UV and optical lines. These results seem to depart from the stratified ionization structures obtained in the classical mapping.

Spectro-timing analysis of Cygnus X-1 during a fast state transition

We present the analysis of two long, quasi-uninterrupted RXTE observations of Cygnus X-1 that span several days within a 10 d interval. The spectral characteristics during this observation cover the region where previous observations have shown the source to be most dynamic. Despite that the source behavior on time scales of hours and days is remarkably similar to that on year time scales. This includes a variety of spectral/temporal correlations that previously had only been observed over Cyg X-1’s long-term evolution. Furthermore, we observe a full transition from a hard to a soft spectral state that occurs within less than 2.5 hours – shorter than previously reported for any other similar Cyg X-1 transition. We describe the spectra with a phenomenological model dominated by a broken power law, and we fit the X-ray variability power spectra with a combination of a cutoff power law and Lorentzian components. The spectral and timing properties are correlated: the power spectrum Lorentzian components have an energy-dependent amplitude, and their peak frequencies increase with photon spectral index. Averaged over 3.2-10 Hz, the time lag between the variability in the 4.5-5.7 keV and 9.5-15 keV bands increases with decreasing hardness when the variability is dominated by the Lorentzian components during the hard state. The lag is small when there is a large power law noise contribution, shortly after the transition to the soft state. Interestingly, the soft state not only shows the shortest lags, but also the longest lags when the spectrum is at its softest and faintest. We discuss our results in terms of emission models for black hole binaries.

Soft time lags in the X-ray emission of Mrk 1040

Temporal analysis of X-ray binaries and Active Galactic Nuclei have shown that hard X-rays react to variation of soft ones after a time delay. The opposite trend, or soft lag, has only been seen in a few rare Quasi-periodic Oscillations in X-ray binaries and recently for the AGN, 1H 0707-495, on short timescales of ~ 10^3 secs. Here, we report analysis of a XMM-Newton observation of Mrk~1040, which reveals that on the dominant variability timescale of ~ 10^4 secs, the source seems to exhibit soft lags. If the lags are frequency independent, they could be due to reverberation effects of a relativistically blurred reflection component responding to a varying continuum. Alternatively, they could be due to Comptonization delays in the case when high energy photons impinge back on the soft photon source. Both models can be verified and their parameters tightly constrained, because they will need to predict the photon spectrum, the r.m.s variability and time lag as a function of energy. A successful application of either model will provide unprecedented information on the radiative process, geometry and more importantly the size of the system, which in turn may provide stringent test of strong general relativistic effects.

Negative X-ray reverberation time delays from MCG-6-30-15 and Mrk 766

We present an X-ray time lag analysis, as a function of Fourier frequency, for MCG-6-30-15 and Mrk 766 using long term XMM-Newton light curves in the 0.5-1.5 keV and the 2-4 keV energy bands, together with some physical modelling of the corresponding time lag spectra. Both the time lag spectra of MCG-6-30-15 and Mrk 766 show negative values (i.e. soft band variations lag behind the corresponding hard band variations) at high frequencies, around 0.001 Hz, similar to those previously observed from 1H 0707-495. The remarkable morphological resemblance between the time lag spectra of MCG-6-30-15 and Mrk 766 indicate that the physical processes responsible for the observed soft time delays is very similar in the two sources, favouring a reflection scenario from material situated very nearby to the central black hole.

An algorithm for Monte-Carlo time-dependent radiation transfer

A new Monte-Carlo algorithm for calculating time-dependent radiative-transfer under the assumption of LTE is presented. Unlike flux-limited diffusion the method is polychromatic, includes scattering, and is able to treat the optically thick and free-streaming regimes simultaneously. The algorithm is tested on a variety of 1-d and 2-d problems, and good agreement with benchmark solutions is found. The method is used to calculate the time-varying spectral energy distribution from a circumstellar disc illuminated by a protostar whose accretion luminosity is varying. It is shown that the time lag between the optical variability and the infrared variability results from a combination of the photon travel time and the thermal response in the disc, and that the lag is an approximately linear function of wavelength.

The updated luminosity correlations of gamma-ray bursts and cosmological implications [Replacement]

Several interesting luminosity correlations among gamma-ray burst (GRB) variables have been recently discussed extensively. In this paper, we derive the six luminosity correlations ($\tlag-L$, $V-L$, $\epkk-L$, $\epkk-E_\gamma$, $\trt-L$, $\epkk-E_{\gamma, \mathrm{iso}}$) from the light curves and spectra of the latest 116 long GRBs, including the time lag ($\tlag$) between low and high photon energy light curves, the variability ($V$) of the light curve, the peak energy of the spectrum ($\epkk$), and the minimum rise time ($\trt$) of the peaks. We find that the intrinsic scatter of the $V-L$ correlation is too large and there seems no inherent correlation between the two parameters using the latest GRB data. The other five correlations indeed exist when the sample is enlarged. The $\epkk-E_\gamma$ correlation has a significantly lower intrinsic scatter compared to the other correlations. We divide the full data into four redshift bins when testing possible evolution of the correlations with redshift. We find no statistically significant evidence for the redshift evolution of the luminosity correlations. To avoid the circularity problem when constraining the cosmological parameters, we simultaneously minimize $\chi^2$ with respect to both correlation parameters $a$, $b$ and the cosmological parameters using the maximum likelihood method. For the flat $\Lambda$CDM, the best fit is $\omm=0.31^{+0.13}_{-0.10}$. We also constrain the possible evolution of the equation of state (EOS) of the dark energy using the GRBs together with the Union2 compilation of SNe Ia and the $H(z)$ data. The result is consistent with the cosmological constant at $2 \sigma$ confidence level and mainly due to the GRB data, the dark energy EOS shows slight deviation from -1 at $z \geq0.5$ as was persistently presented with many previous data sets.

The updated luminosity correlations of gamma-ray bursts and cosmological implications

Several interesting luminosity correlations among gamma-ray burst (GRB) variables have been recently discussed extensively. In this paper, we derive the six luminosity correlations ($\tlag-L$, $V-L$, $\epkk-L$, $\epkk-E_\gamma$, $\trt-L$, $\epkk-E_{\gamma, \mathrm{iso}}$) from the light curves and spectra of the latest 116 long GRBs, including the time lag ($\tlag$) between low and high photon energy light curves, the variability ($V$) of the light curve, the peak energy of the spectrum ($\epkk$), and the minimum rise time ($\trt$) of the peaks. We find that the intrinsic scatter of the $V-L$ correlation is too large and there seems no inherent correlation between the two parameters using the latest GRB data. The other five correlations indeed exist when the sample is enlarged. The $\epkk-E_\gamma$ correlation has a significantly lower intrinsic scatter compared to the other correlations. We divide the full data into four redshift bins when testing possible evolution of the correlations with redshift. We find no statistically significant evidence for the redshift evolution of the luminosity correlations. To avoid the circularity problem when constraining the cosmological parameters, we simultaneously minimize $\chi^2$ with respect to both correlation parameters $a$, $b$ and the cosmological parameters using the maximum likelihood method. For the flat $\Lambda$CDM, the best fit is $\omm=0.31^{+0.13}_{-0.10}$. We also constrain the possible evolution of the equation of state (EOS) of the dark energy using the GRBs together with the Union2 compilation of SNe Ia and the $H(z)$ data. The result is consistent with the cosmological constant at $2 \sigma$ confidence level and mainly due to the GRB data, the dark energy EOS shows slight deviation from -1 at $z \geq0.5$ as was persistently presented with many previous data sets.

The causal connection between disc and power-law variability in hard state black hole X-ray binaries

We use the XMM-Newton EPIC-pn instrument in timing mode to extend spectral time-lag studies of hard state black hole X-ray binaries into the soft X-ray band. We show that variations of the disc blackbody emission substantially lead variations in the power-law emission, by tenths of a second on variability time-scales of seconds or longer. The large lags cannot be explained by Compton scattering but are consistent with time-delays due to viscous propagation of mass accretion fluctuations in the disc. However, on time-scales less than a second the disc lags the power-law variations by a few ms, consistent with the disc variations being dominated by X-ray heating by the power-law, with the short lag corresponding to the light-travel time between the power-law emitting region and the disc. Our results indicate that instabilities in the accretion disc are responsible for continuum variability on time-scales of seconds or longer and probably also on shorter time-scales.

Water vapour at high redshift: Arecibo monitoring of the megamaser in MG J0414+0534

The study of water masers at cosmological distances would allow us to investigate the parsec-scale environment around powerful radio sources, to probe the physical conditions of the molecular gas in the inner parsecs of quasars, and to estimate their nuclear engine masses in the early universe. To derive this information, the nature of the maser source, jet or disk-maser, needs to be assessed through a detailed investigation of the observational characteristics of the line emission. We monitored the maser line in the lensed quasar MGJ0414+0534 at z = 2.64 with the 300-m Arecibo telescope for ~15 months to detect possible additional maser components and to measure a potential velocity drift of the lines. In addition, we follow the maser and continuum emissions to reveal significant variations in their flux density and to determine correlation or time-lag, if any, between them. The main maser line profile is complex and can be resolved into a number of broad features with line widths of 30-160 km/s. A new maser component was tentatively detected in October 2008 that is redshifted by 470 km/s w.r.t the systemic velocity of the quasar. The line width of the main maser feature increased by a factor of two between the Effelsberg and EVLA observations reported by Impellizzeri et al. (2008) and the first epoch of the Arecibo monitoring campaign. After correcting for the lens magnification, we find that the total H2O isotropic luminosity of the maser in MGJ0414+0534 is now ~30,000 Lsun, making this source the most luminous ever discovered.[Abridged]

Weighing Super-Massive Black Holes with Narrow Fe K$\alpha$ Line [Replacement]

It has been suggested that the narrow cores of the Fe K$\alpha$ emission lines in Active Galactic Nuclei (AGNs) are likely produced in the torus, the inner radius of which can be measured by observing the lag time between the $V$ and $K$ band flux variations. In this paper we compare the virial products of the infrared time lags and the narrow Fe K$\alpha$ widths for 10 type 1 AGNs with the black hole masses from other techniques. We find the narrow Fe K$\alpha$ line width is in average 2.6$^{+0.9}_{-0.4}$ times broader than expected assuming an isotropic velocity distribution of the torus at the distance measured by the infrared lags. We propose the thick disk model of the torus could explain the observed larger line width. Another possibility is the contamination by emission from the broad line region or the outer accretion disk. Alternatively, the narrow iron line might originate from the inner most part of the obscuring torus within the sublimation radius, while the infrared emission from outer cooler part. We note the correlation between the black hole masses based on this new technique and those based on other known techniques is statistically insignificant. We argue that this could be attributed to the small sample size and the very large uncertainties in the measurements of iron K line widths. The next generation of X-ray observatories could help verify the origin of the narrow iron K$\alpha$ line and the reliability of this new technique.

Weigh the Super-Massive Black Holes with Narrow Fe K$\alpha$ Line?

It has been suggested that the narrow cores of the Fe K$\alpha$ emission lines in Active Galactic Nuclei (AGNs) are likely produced in the torus, the inner radius of which can be measured by observing the lag time between the $V$ and $K$ band flux variations. In this paper we compare the virial products of the infrared time lags and the narrow Fe K$\alpha$ widths for 10 type 1 AGNs with the black hole masses from other techniques. We find the narrow Fe K$\alpha$ line width is in average 2.6$^{+0.9}_{-0.4}$ times broader than expected assuming an isotropic velocity distribution of the torus at the distance measured by the infrared lags. We propose the thick disk model of the torus could explain the observed larger line width. Another possibility is the contamination by emission from the broad line region or the outer accretion disk. Alternatively, the narrow iron line might originate from the inner most part of the obscuring torus within the sublimation radius, while the infrared emission from outer cooler part. We note the correlation between the black hole masses based on this new technique and those based on other known techniques is statistically insignificant. We argue that this could be attributed to the small sample size and the very large uncertainties in the measurements of iron K line widths. The next generation of X-ray observatories could help verify the origin of the narrow iron K$\alpha$ line and the reliability of this new technique.

Simultaneous monitoring of the photometric and polarimetric activity of the young star PV Cep in the optical/near-infrared bands

We present the results of a simultaneous monitoring, lasting more than 2 years, of the optical and near-infrared photometric and polarimetric activity of the variable protostar PV Cep. During the monitoring period, an outburst has occurred in all the photometric bands, whose declining phase ($\Delta$J $\approx$ 3 mag) lasted about 120 days. A time lag of $\sim$ 30 days between optical and infrared light curves has been measured and interpreted in the framework of an accretion event. This latter is directly recognizable in the significant variations of the near-infrared colors, that appear bluer in the outburst phase, when the star dominates the emission, and redder in declining phase, when the disk emission prevails. All the observational data have been combined to derive a coherent picture of the complex morphology of the whole PV Cep system, that, in addition to the star and the accretion disk, is composed also by a variable biconical nebula. In particular, the mutual interaction between all these components is the cause of the high value of the polarization ($\approx$ 20%) and of its fluctuations. The observational data concur to indicate that PV Cep is not a genuine EXor star, but rather a more complex object; moreover the case of PV Cep leads to argue about the classification of other recently discovered young sources in outburst, that have been considered, maybe over-simplifying, as EXor.

The 2009 outburst of IGR J17511-3057 as observed by SWIFT and RXTE [Replacement]

The accretion-powered millisecond pulsar IGR J17511-3057 was discovered in September 2009. This class of sources is growing steadily and examining the recently obtained data is important to reveal and understand their observational properties and underlying physics. Our spectral analysis of the source indicated only slight spectral shape evolution during the entire outburst. Pulse profile analysis revealed smooth pulsations with moderate decrease in amplitude and absence of dramatic shape evolution, however we have found that the time lag between soft and hard energy pulses increases by factor of two during the outburst. We modelled the lag effect and concluded that the only likely cause of this evolution is change of the Comptonized radiation emissivity pattern. It is also important to trace dramatic changes in the object characteristics, as such facts are necessary to understand the physics of accretion in these objects. We have found that, similarly to its progenitor SAX J1808.4-3658, this source demonstrates change of outburst stages from so-called “slow decay” into “rapid drop”, as we can see from the dramatic flux drop and pulse profile evolution at the very end of the outburst.

The 2009 outburst of IGR J17511-3057 as observed by SWIFT and RXTE

The accretion-powered millisecond pulsar IGR J17511-3057 was discovered in September 2009. This class of sources is growing steadily and examining the recently obtained data is important to reveal and understand their observational properties and underlying physics. Our spectral analysis of the source indicated only slight spectral shape evolution during the entire outburst. Pulse profile analysis revealed smooth pulsations with moderate decrease in amplitude and absence of dramatic shape evolution, however we have found that the time lag between soft and hard energy pulses increases by factor of two during the outburst. We modelled the lag effect and concluded that the only likely cause of this evolution is change of the Comptonized radiation emissivity pattern. It is also important to trace dramatic changes in the object characteristics, as such facts are necessary to understand the physics of accretion in these objects. We have found that, similarly to its progenitor SAX J1808.4-3658, this source demonstrates change of outburst stages from so-called “slow decay” into “rapid drop”, as we can see from the dramatic flux drop and pulse profile evolution at the very end of the outburst.

Near-infrared polarimetry as a tool for testing properties of accreting super-massive black holes

Several massive black holes exhibit flux variability on time scales that correspond to source sizes of the order of few Schwarzschild radii. We survey the potential of near-infrared and X-ray polarimetry to constrain physical properties of such black hole systems, namely their spin and inclination. We have focused on a model where an orbiting hot spot is embedded in an accretion disk. A new method of searching for the time-lag between orthogonal polarization channels is developed and applied to an ensemble of hot spot models that samples a wide range of parameter space. We found that the hot spot model predicts signatures in polarized light which are in the range to be measured directly in the near future. However, our estimations are predicted upon the assumption of a Keplerian velocity distribution inside the flow where the dominant part of the magnetic field is toroidal. We also found that if the right model of the accretion flow can be chosen for each source (e.g. on the basis of magnetohydrodynamics simulations) then the black hole spin and inclination can be constrained to a small two-dimensional area in the spin-inclination space. The results of the application of the method to the available near-infrared polarimetric data of Sagittarius A* (Sgr A*) is presented. It is shown that even with the currently available data the spin and inclination of Sgr A* can be constrained. Next generations of near-infrared and X-ray polarimeters should be able to exploit this tool.

The phase relation between sunspot numbers and soft X-ray flares [Replacement]

To better understand long-term flare activity, we present a statistical study on soft X-ray flares from May 1976 to May 2008. It is found that the smoothed monthly peak fluxes of C-class, M-class, and X-class flares have a very noticeable time lag of 13, 8, and 8 months in cycle 21 respectively with respect to the smoothed monthly sunspot numbers. There is no time lag between the sunspot numbers and M-class flares in cycle 22. However, there is a one-month time lag for C-class flares and a one-month time lead for X-class flares with regard to sunspot numbers in cycle 22. For cycle 23, the smoothed monthly peak fluxes of C-class, M-class, and X-class flares have a very noticeable time lag of one month, 5 months, and 21 months respectively with respect to sunspot numbers. If we take the three types of flares together, the smoothed monthly peak fluxes of soft X-ray flares have a time lag of 9 months in cycle 21, no time lag in cycle 22 and a characteristic time lag of 5 months in cycle 23 with respect to the smoothed monthly sunspot numbers. Furthermore, the correlation coefficients of the smoothed monthly peak fluxes of M-class and X-class flares and the smoothed monthly sunspot numbers are higher in cycle 22 than those in cycles 21 and 23. The correlation coefficients between the three kinds of soft X-ray flares in cycle 22 are higher than those in cycles 21 and 23. These findings may be instructive in predicting C-class, M-class, and X-class flares regarding sunspot numbers in the next cycle and the physical processes of energy storage and dissipation in the corona.

The phase relation between sunspot numbers and soft X-ray flares

To better understand long-term flare activity, we present a statistical study on soft X-ray flares from May 1976 to May 2008. It is found that the smoothed monthly peak fluxes of C-class, M-class, and X-class flares have a very noticeable time lag of 13, 8, and 8 months in cycle 21 respectively with respect to the smoothed monthly sunspot numbers. There is no time lag between the sunspot numbers and M-class flares in cycle 22. However, there is a one-month time lag for C-class flares and a one-month time lead for X-class flares with regard to sunspot numbers in cycle 22. For cycle 23, the smoothed monthly peak fluxes of C-class, M-class, and X-class flares have a very noticeable time lag of one month, 5 months, and 21 months respectively with respect to sunspot numbers. If we take the three types of flares together, the smoothed monthly peak fluxes of soft X-ray flares have a time lag of 9 months in cycle 21, no time lag in cycle 22 and a characteristic time lag of 5 months in cycle 23 with respect to the smoothed monthly sunspot numbers. Furthermore, the correlation coefficients of the smoothed monthly peak fluxes of M-class and X-class flares and the smoothed monthly sunspot numbers are higher in cycle 22 than those in cycles 21 and 23. The correlation coefficients between the three kinds of soft X-ray flares in cycle 22 are higher than those in cycles 21 and 23. These findings may be instructive in predicting C-class, M-class, and X-class flares regarding sunspot numbers in the next cycle and the physical processes of energy storage and dissipation in the corona.

An unified polar cap/striped wind model for pulsed radio and gamma-ray emission in pulsars

(abridged) Thanks to the recent discovery by Fermi of about fifty new gamma-ray pulsars, it becomes possible to look for statistical properties of their pulsed high-energy emission, especially their light-curves and phase-resolved spectra. These pulsars emit by definition mostly gamma-ray photons but some of them are also detected in the radio band. For those seen in these two extreme energies, the relation between time lag of radio/gamma-ray pulses and gamma-ray peak separation, in case both high-energy pulses are seen, helps to put some constrain on the magnetospheric emission mechanisms and location. This idea is analyzed in detail in this paper, assuming a polar cap model for the radio pulses and the striped wind geometry for the pulsed high-energy counterpart. Combining the time-dependent emissivity in the wind, supposed to be inverse Compton radiation, with a simple polar cap emission model along and around the magnetic axis, we compute the radio and gamma-ray light-curves, summarizing the results in several phase plots. The phase lag as well as the gamma-ray peak separation dependence on the pulsar inclination angle and on the viewing angle are studied. Using the gamma-ray pulsar catalog compiled from the Fermi data, we are able to predict the radio lag/peak separation relation and compare it with available observations taken from this catalog.

The Jet in the Galactic Center: An Ideal Laboratory for Magnetohydrodynamics and General Relativity

In this paper we review and discuss some of the intriguing properties of the Galactic Center supermassive black hole candidate Sgr A*. Of all possible black hole sources, the event horizon of Sgr A*, subtends the largest angular scale on the sky. It is therefore a prime candidate to study and image plasma processes in strong gravity and it even allows imaging of the shadow cast by the event horizon. Recent mm-wave VLBI and radio timing observations as well as numerical GRMHD simulations now have provided several breakthroughs that put Sgr A* back into the focus. Firstly, VLBI observations have now measured the intrinsic size of Sgr A* at multiple frequencies, where the highest frequency measurements have approached the scale of the black hole shadow. Moreover, measurements of the radio variability show a clear time lag between 22 GHz and 43 GHz. The combination of size and timing measurements, allows one to actually measure the flow speed and direction of magnetized plasma at some tens of Schwarzschild radii. This data strongly support a moderately relativistic outflow, consistent with an accelerating jet model. This is compared to recent GRMHD simulation that show the presence of a moderately relativistic outflow coupled to an accretion flow Sgr A*. Further VLBI and timing observations coupled to simulations have the potential to map out the velocity profile from 5-40 Schwarzschild radii and to provide a first glimpse at the appearance of a jet-disk system near the event horizon. Future submm-VLBI experiments would even be able to directly image those processes in strong gravity and directly confirm the presence of an event horizon.

Jets at lowest mass accretion rates

We present results of recent observations and theoretical modeling of data from black holes accreting at very low luminosities (L/L_Edd ~ 10^{-8}). We discuss our newly developed time-dependent model for episodic ejection of relativistic plasma within a jet framework, and a successful application of this model to describe the origin of radio flares seen in Sgr A*, the Galactic center black hole. Both the observed time lags and size-frequency relationships are reproduced well by the model. We also discuss results from new Spitzer data of the stellar black hole X-ray binary system A0620-00. Complemented by long term SMARTS monitoring, these observations indicate that once the contribution from the accretion disk and the donor star are properly included, the residual mid-IR spectral energy distribution of A0620-00 is quite flat and consistent with a non-thermal origin. The results above suggest that a significant fraction of the observed spectral energy distribution originating near black holes accreting at low luminosities could result from a mildly relativistic outflow. The fact that these outflows are seen in both stellar-mass black holes as well as in supermassive black holes at the heart of AGNs strengthens our expectation that accretion and jet physics scales with mass.

The transition from quasar radio-loud to radio-quiet state in the framework of the black hole scalability hypothesis [Replacement]

There are several lines of evidence that active galactic nuclei (AGN) can be regarded as scaled-up X-ray binaries (XRB). The timescales of the evolutionary phenomena in these two classes are proportional to the black hole (BH) masses. Consequently, unlike in the case of XRBs, the evolution of AGNs is too slow to be followed directly. What could be done, however, is to assign particular types of active galaxies to different evolutionary stages observable in XRBs. We studied such an assignment for three quasars with clear signatures of a recent transition from the radio-loud to the radio-quiet state. The quasars we investigated have large-scale radio lobes that are clearly asymmetric — one lobe is of Fanaroff-Riley II type, while the other one is a diffuse relic devoid of a hotspot. We suggest that the prime cause of the asymmetry of these radio sources is that the nuclei of their host galaxies currently produce no jets. To prove that, we observed them with milliarcsecond resolution to check if they are similar to those in radio-quiet quasars. The observations carried out with the EVN revealed that the nuclei of the quasars under investigation are not of a core-jet type that is characteristic for radio-loud, lobe-dominated quasars. It follows that the lobes are no longer fuelled and that the apparent asymmetry results from the orientation, which causes a time lag of the order of 10^6 years between their images: the lobe perceived as a relic is nearer than the lobe with a hotspot and so it is observed in a later stage of the decay. The three AGNs under investigation were radio-loud earlier, but now they have switched to the radio-quiet state. In the framework of the XRB/AGN unification, the above means that they have left the very high state and have moved now to the high/soft state. (abridged)

The transition from quasar radio-loud to radio-quiet state in the framework of the black hole scalability hypothesis

There are several lines of evidence that active galactic nuclei (AGN) can be regarded as scaled-up X-ray binaries (XRB). The timescales of the evolutionary phenomena in these two classes are proportional to the black hole (BH) masses. Consequently, unlike in the case of XRBs, the evolution of AGNs is too slow to be followed directly. What could be done, however, is to assign particular types of active galaxies to different evolutionary stages observable in XRBs. We studied such an assignment for three quasars with clear signatures of a recent transition from the radio-loud to the radio-quiet state. The quasars we investigated have large-scale radio lobes that are clearly asymmetric — one lobe is of Fanaroff-Riley II type, while the other one is a diffuse relic devoid of a hotspot. We suggest that the prime cause of the asymmetry of these radio sources is that the nuclei of their host galaxies currently produce no jets. To prove that, we observed them with milliarcsecond resolution to check if they are similar to those in radio-quiet quasars. The observations carried out with the EVN revealed that the nuclei of the quasars under investigation are not of a core-jet type that is characteristic for radio-loud, lobe-dominated quasars. It follows that the lobes are no longer fuelled and that the apparent asymmetry results from the orientation, which causes a time lag of the order of 10^6 years between their images: the lobe perceived as a relic is nearer than the lobe with a hotspot and so it is observed in a later stage of the decay.The three AGNs under investigation were radio-loud earlier, but now they have switched to the radio-quiet state. In the framework of the XRB/AGN unification, the above means that they have left the very high state and have moved now to the high/soft state. (abridged)

Observational evidence for matter propagation in accretion flows

We study simultaneous X-ray and optical observations of three intermediate polars EX Hya, V1223 Sgr and TV Col with the aim to understand the propagation of matter in their accretion flows. We show that in all cases the power spectra of flux variability of binary systems in X-rays and in optical band are similar to each other and the majority of X-ray and optical fluxes are correlated with time lag <1 sec. These findings support the idea that optical emission of accretion disks, in these binary systems,largely originates as reprocessing of X-ray luminosity of their white dwarfs. In the best obtained dataset of EX Hya we see that the optical lightcurve unambiguously contains some component, which leads the X-ray emission by ~7 sec. We interpret this in the framework of the model of propagating fluctuations and thus deduce the time of travel of the matter from the innermost part of the truncated accretion disk to the white dwarf surface. This value agrees very well with the time expected for matter threaded onto the magnetosphere of the white dwarf to fall to its surface. The datasets of V1223 Sgr and TV Col in general confirm these findings,but have poorer quality.

The bulk of the black hole growth since z~1 occurs in a secular universe: No major merger-AGN connection [Replacement]

What is the relevance of major mergers and interactions as triggering mechanisms for active galactic nuclei (AGN) activity? To answer this longstanding question, we analyze 140 XMM-selected AGN host galaxies and a matched control sample of 1264 inactive galaxies over z~0.3-1.0 and log(M_*/M_sun)85%) of the AGN host galaxies do not show strong distortions, and (2) there is no significant difference in the distortion fractions between active and inactive galaxies. Our findings provide the best direct evidence that, since z~1, the bulk of black hole accretion has not been triggered by major galaxy mergers, therefore arguing that the alternative mechanisms, i.e., secular processes and minor interactions, are the leading triggers for the episodes of major black hole growth. We also exclude an alternative interpretation of our results: a significant time lag between merging and the observability of the AGN phase could wash out the most significant merging signatures, explaining the lack of enhancement of strong distortions on the AGN hosts. We show that this alternative scenario is unlikely due to: (1) recent major mergers being ruled out for the majority of sources due to the high fraction of disk-hosted AGN, (2) the lack of a significant X-ray signal in merging inactive galaxies as a signature of a potential buried AGN, and (3) the low levels of soft X-ray obscuration for AGN hosted by interacting galaxies, in contrast to model predictions.

The bulk of the black hole growth since z~1 occurs in a secular universe: No merger-AGN connection

What is the relevance of major mergers and interactions as triggering mechanisms for active galactic nuclei (AGN) activity? To answer this longstanding question, we analyze 140 XMM-selected AGN host galaxies and a matched control sample of 1264 inactive galaxies over z~0.3-1.0 and log(M_*/M_sun)<11.7 with high-resolution HST/ACS imaging from the COSMOS field. The visual analysis of their morphologies by 10 independent human classifiers yields a measure of the fraction of distorted morphologies in the AGN and control samples, i.e. quantifying the signature of recent mergers which might potentially be responsible for fueling/triggering the AGN. We find that (i) the vast majority (>85%) of the AGN host galaxies do not show strong distortions, and (ii) there is no significant difference in the distortion fractions between active and inactive galaxies. Our findings provide the best direct evidence that, since z~1, the bulk of black hole accretion has not been triggered by major galaxy mergers, therefore arguing that the alternative mechanisms, i.e., secular processes and minor interactions, are the leading triggers for the episodes of major black hole growth. We also exclude an alternative interpretation of our results: a significant time lag between merging and the observability of the AGN phase could wash out the most significant merging signatures, explaining the lack of enhancement of strong distortions on the AGN hosts. We show that this alternative scenario is unlikely due to: (i) recent major mergers being ruled out for the majority of sources due to the high fraction of disk-hosted AGN, (ii) the lack of a hidden X-ray signal in merging inactive galaxies as a signature of a potential buried AGN, and (iii) the low levels of soft X-ray obscuration for AGN hosted by interacting galaxies, in contrast to model predictions.

X-ray and radio variability in the low luminosity Active Galactic Nucleus NGC 7213

We present the results of a ~ 3 year campaign to monitor the low luminosity active galactic nucleus (LLAGN) NGC 7213 in the radio (4.8 and 8.4 GHz) and X-ray bands (2-10 keV). With a reported X-ray Eddington ratio of 7 x 10^-4 L_Edd, NGC 7213 can be considered to be comparable to a hard state black hole X-ray binary. We show that a weak correlation exists between the X-ray and radio light curves. We use the cross-correlation function to calculate a global time lag between events in the X-ray and radio bands to be 24 +/- 12 days lag (8.4 GHz radio lagging X-ray), and 40 +/- 13 days lag (4.8 GHz radio lagging X-ray). The radio-radio light curves are extremely well correlated with a lag of 20.5 +/- 12.9 days (4.8 GHz lagging 8.4 GHz). We explore the previously established scaling relationship between core radio and X-ray luminosities and black hole mass L_{R} \propto M^{0.6-0.8} L_{X}^{0.6}, known as the `fundamental plane of black hole activity’, and show that NGC 7213 lies very close to the best-fit `global’ correlation for the plane as one of the most luminous LLAGN. With a large number of quasi-simultaneous radio and X-ray observations, we explore for the first time the variations of a single AGN with respect to the fundamental plane. Although the average radio and X-ray luminosities for NGC 7213 are in good agreement with the plane, we show that there is intrinsic scatter with respect to the plane for the individual data points.

An Alternative Approach To Measuring Reverberation Lags in Active Galactic Nuclei [Replacement]

Motivated by recent progress in the statistical modeling of quasar variability, we develop a new approach to measuring emission-line reverberation lags to estimate the size of broad-line regions (BLRs) in active galactic nuclei. Assuming that all emission-line light curves are scaled, smoothed, and displaced versions of the continuum, this alternative approach fits the light curves directly using a damped random walk model and aligns them to recover the time lag and its statistical confidence limits. We introduce the mathematical formalism of this approach and demonstrate its ability to cope with some of the problems for traditional methods, such as irregular sampling, correlated errors, and seasonal gaps. We redetermine the lags for 87 emission lines in 31 quasars and reassess the BLR size–luminosity relationship using 60 H-beta lags. We confirm the general results from the traditional cross-correlation methods, with a few exceptions. Our method, however, also supports a broad range of extensions. In particular, it can simultaneously fit multiple lines and continuum light curves which improves the lag estimate for the lines and provides estimates of the error correlations between them. Determining these correlations is of particular importance for interpreting emission-line velocity–delay maps. We can also include parameters for luminosity-dependent lags or line responses. We use this to detect the scaling of the BLR size with continuum luminosity in NGC 5548.

A New Approach To Reverberation Mapping

Motivated by recent progress in the statistical modeling of quasar variability, we develop a new approach to reverberation mapping for estimating the size of broad-line regions (BLRs) in AGNs. Assuming that all light curves are scaled, smoothed, and displaced versions of the continuum, the new approach fits the light curves directly using a damped random walk model and aligns them to recover the time lag and its statistical confidence limits. We introduce the mathematical formalism of the new approach and demonstrate its ability to cope with some of the problems for traditional methods, such as irregular sampling, correlated errors, and seasonal gaps. We redetermine the lags for 86 emission lines in 31 quasars and estimate a new BLR size-luminosity relationship using 59 H\beta lags. On a positive note, we confirm the general results from the traditional methods, with few exceptions. Our method, however, also supports a broad range of extensions. In particular, it can simultaneously fit multiple lines and continuum light curves which improves the lag estimate for the lines and provides estimates of the error correlations between them. Determining these correlation is of particular importance for interpreting line velocity-delay maps. We can also include parameters for luminosity-dependent lags or line responses. We use this to detect the scaling of the BLR size with continuum luminosity in NGC 5548.

Lorentz Symmetry breaking studies with photons from astrophysical observations [Replacement]

Lorentz Invariance Violation (LIV) may be a good observational window on Quantum Gravity physics. Within last few years, all major Gamma-ray experiments have published results from the search for LIV with variable astrophysical sources: gamma-ray bursts with detectors on-board satellites and Active Galactic Nuclei with ground-based experiments. In this paper, the recent time-of-flight studies with unpolarized photons published from the space and ground based observations are reviewed. Various methods used in the time delay searches are described, and their performance discussed. Since no significant time-lag value was found within experimental precision of the measurements, the present results consist of 95% confidence cevel limits on the Quantum Gravity scale on the linear and quadratic terms in the standard photon dispersion relations.

Lorentz Symmetry Breaking studies with photons from Astrophysical Observations

Lorentz Invariance Violation (LIV) may be a good observational window on Quantum Gravity physics. Within last few years, all major Gamma-ray experiments have published results from the search for LIV with variable astrophysical sources: Gamma-ray Bursts with detectors on-board satellites and Active Galactic Nuclei with ground-based experiments. In this paper, the recent time-of-flight studies with unpolarized photons published from the space and ground based observations are reviewed. Various methods used in the time delay searches are described, and their performance discussed. Since no significant time-lag value was found within experimental precision of the measurements, the present results consist of 95% Confidence Level limits on the Quantum Gravity scale on the linear and quadratic terms in the standard photon dispersion relations.

Magnetospheric accretion and spin-down of the prototypical classical T Tauri star AATau

From observations collected with the ESPaDOnS & NARVAL spectropolarimeters at CFHT and TBL, we report the detection of Zeeman signatures on the prototypical classical TTauri star AATau, both in photospheric lines and accretion-powered emission lines. Using time series of unpolarized and circularly polarized spectra, we reconstruct at two epochs maps of the magnetic field, surface brightness and accretion-powered emission of AATau. We find that AATau hosts a 2-3kG magnetic dipole tilted at ~20deg to the rotation axis, and of presumably dynamo origin. We also show that the magnetic poles of AATau host large cool spots at photospheric level and accretion regions at chromospheric level. The logarithmic accretion rate at the surface of AATau at the time of our observations is strongly variable, ranging from -9.6 to -8.5 and equal to -9.2 in average (in Msun/yr); this is an order of magnitude smaller than the disc accretion rate at which the magnetic truncation radius (below which the disc is disrupted by the stellar magnetic field) matches the corotation radius (where the Keplerian period equals the stellar rotation period) – a necessary condition for accretion to occur. It suggests that AATau is largely in the propeller regime, with most of the accreting material in the inner disc regions being expelled outwards and only a small fraction accreted towards the surface of the star. The strong variability in the observed surface mass-accretion rate and the systematic time-lag of optical occultations (by the warped accretion disc) with respect to magnetic and accretion-powered emission maxima also support this conclusion. Our results imply that AATau is being actively spun-down by the star-disc magnetic coupling and appears as an ideal laboratory for studying angular momentum losses of forming Suns in the propeller regime.

The peculiar optical spectrum of 4C+22.25: Imprint of a massive black hole binary?

We report the discovery of peculiar features in the optical spectrum of 4C+22.25, a flat spectrum radio quasar at z=0.4183 observed in the SDSS and in a dedicated spectroscopic follow-up from the Nordic Optical Telescope. The Hbeta and Halpha lines show broad profiles (FWHM~12,000 km/s), faint fluxes and extreme offsets (Delta v=8,700+/-1,300 km/s) with respect to the narrow emission lines. These features show no significant variation in a time lag of ~3.1 yr (rest frame). We rule out possible interpretations based on the superposition of two sources or on recoiling black holes, and we discuss the virtues and limitations of a massive black hole binary scenario.

X-ray reverberation in 1H0707-495 revisited

The narrow-line Seyfert 1 galaxy 1H0707-495 has previously been identified as showing time lags between flux variations in the soft- (0.3-1 keV) and medium-energy (1-4 keV) X-ray bands that oscillate between positive and negative values as a function of the frequency of the mode of variation. Here we measure and analyse the lags also between a harder X-ray band (4-7.5 keV) and the soft and medium bands, using existing XMM-Newton data, and demonstrate that the entire spectrum of lags, considering both the full energy range, 0.3-7.5 keV, and the full frequency range, 10^-5 < nu < 10^-2 Hz, are inconsistent with previous claims of arising as reverberation associated with the inner accretion disk. Instead we demonstrate that a simple reverberation model, in which scattering or reflection is present in all X-ray bands, explains the full set of lags without requiring any ad hoc explanation for the time lag sign changes. The range of time delays required to explain the observed lags extends up to about 1800 s in the hard band. The results are consistent with reverberation caused by scattering of X-rays passing through an absorbing medium whose opacity decreases with increasing energy and that partially-covers the source. A high covering factor of absorbing and scattering circumnuclear material is inferred.

Analytical error estimate for the cross-correlation, phase and time lag between two light curves

Temporal analysis of radiation from Astrophysical sources like Active Galactic Nuclei, X-ray Binaries and Gamma-ray bursts provide information on the geometry and sizes of the emitting regions. Robustly establishing that two light-curves in different energy bands are correlated and measuring the phase and time-lag between them is an important and frequently used temporal diagnostic. Analytical expressions to estimate the errors on the cross-correlation, phase and time-lag between two light-curves are presented. Earlier estimates depended upon numerically expensive simulations or on dividing the light-curves in large number of segments to find the variance. Thus, the analytical estimates presented here allow for analysis of light-curves with relatively small (~ 1000) number of points, as well as to obtain information on the longest time-scales available. The error estimation is verified using simulations of light-curves derived from both white and 1/f stochastic processes with measurement errors. As a demonstration, we apply this technique to the XMM-Newton light-curves of the Active Galactic Nucleus, Akn 564.

Thermal Radiation from GRB Jets [Replacement]

In this study, the light curves and spectrum of the photospheric thermal radiation from ultrarelativistic gamma-ray burst (GRB) jets are calculated using 2D relativistic hydrodynamic simulations of jets from a collapsar. As the jet advances, the density around the head of the jet decreases, and its Lorentz factor reaches as high as 200 at the photosphere and 400 inside the photosphere. For an on-axis observer, the photosphere appears concave shaped due to the low density and high beaming factor of the jet. The luminosity varies because of the abrupt change in the position of the photosphere due to the internal structure of the jet. Comparing our results with GRB090902B, the flux level of the thermal-like component is similar to our model, although the peak energy looks a little bit higher (but still within a factor of 2). From the comparison, we estimate that the bulk Lorentz factor of GRB090902B is $\Gamma \sim 2.4 \times 10^2 (r/10^{12} \rm cm$) where $r$ is the radius of the photosphere. The spectrum for an on-axis observer is harder than that for an off-axis observer. There is a time lag of a few seconds for high energy bands in the light curve. This may be the reason for the delayed onset of GeV emission seen in GRB080916C. The spectrum below the peak energy is a power law and the index is $2.3 \sim 2.6$ which is softer than that of single temperature plank distribution but still harder than that of typical value of observed one.

Thermal Radiation from GRB Jets [Replacement]

In this study, the light curves and spectrum of the photospheric thermal radiation from ultrarelativistic gamma-ray burst (GRB) jets are calculated using 2D relativistic hydrodynamic simulations of jets from a collapsar. As the jet advances, the density around the head of the jet decreases, and its Lorentz factor reaches as high as 200 at the photosphere and 400 inside the photosphere. For an on-axis observer, the photosphere appears concave shaped due to the low density and high beaming factor of the jet. The light curve for the first phases is flared for all viewing angles. For an on-axis observer, luminosity varies because of the abrupt change in the position of the photosphere due to the internal structure of the jet. Comparing our results with GRB090902B, the flux level of the thermal(-like) component is similar to our model, although the peak energy looks a little bit higher (but still within a factor of $2$). From the comparison, we estimate that the bulk Lorentz factor of GRB090902B is $\Gamma \sim 2.4 \times 10^2 (r/10^{12} \rm cm$) where $r$ is the radius of the photosphere. The spectrum for an on-axis observer is harder than that for an off-axis observer. There is a time lag of a few seconds for high energy bands in the light curve. This may be the reason for the delayed onset seen in GRB080916C. The spectrum below the peak energy is power law and the index is $2 \sim 2.6$ which is softer than that of single temperature plank distribution but still harder than that of typical value of observed one.

Thermal Radiation from GRB Jets [Replacement]

In this study, the light curves and spectrum of the photospheric thermal radiation from ultrarelativistic gamma-ray burst (GRB) jets are calculated using 2D relativistic hydrodynamic simulations of jets from a collapsar. As the jet advances, the density around its head of the jet decreases, and its Lorentz factor reaches as high as 200 at the photosphere and 400 inside the photosphere. For an on-axis observer, the photosphere appears concave shaped due to the low density and high beaming factor of the jet. The light curve for the first phases is flared for all viewing angles. For an on-axis observer, luminosity varies because of the abrupt change in the position of the photosphere due to the internal structure of the jet. Comparing our results with those of GRB090902B, we conclude that GRB090902B had a more collimated jet with a higher temperature and/or beaming factor than that of our simulation. The spectrum for an on-axis observer is harder than that for an off-axis observer. There is a time lag of a few seconds for high energy bands in the light curve, which may be the reason for the delayed onset seen in GRB080916C.

Thermal Radiation from GRB Jets [Replacement]

In this study, the light curves and spectrum of the photospheric thermal radiation from ultrarelativistic gamma-ray burst (GRB) jets are calculated using 2D relativistic hydrodynamic simulations of jets from a collapsar. As the jet advances, the density around its head of the jet decreases, and its Lorentz factor reaches as high as 200 at the photosphere and 400 inside the photosphere. For an on-axis observer, the photosphere appears concave shaped due to the low density and high beaming factor of the jet. The light curve for the first phases is flared for all viewing angles. For an on-axis observer, luminosity varies because of the abrupt change in the position of the photosphere due to the internal structure of the jet. Comparing our results with those of GRB090902B, we conclude that GRB090902B had a more collimated jet with a higher temperature and/or beaming factor than that of our simulation. The spectrum for an on-axis observer is harder than that for an off-axis observer. There is a time lag of a few seconds for high energy bands in the light curve, which may be the reason for the delayed onset seen in GRB080916C.

Thermal Radiation from GRB Jets

The light curves and spectrum of the photospheric thermal radiation from ultrarelativistic gamma-ray burst (GRB) jets are calculated, using the results of the 2D relativistic hydrodynamic simulations of the jets from a collapsar. The density around the head of the jet decreases, as the jet proceeds, and the Lorentz factor of the jet reaches as high as 200 at the photospheer and as high as 400 inside the photosphere. The shape of the photosphere for on-axis observer gets concave due to low density and high beaming factor of the jet. The light curve is flare like for the first phases for all viewing angles. Then the luminosity for on-axis observer varies caused by the quick transition of the position of the photosphere due to the internal structure in the jet. We compare our results with GRB090902B, and conclude that GRB090902B had a more collimated jet with higher temperature and/or beaming factor than our simulation. The spectrum for on-axis observer is harder than that for off-axis observer. There is a few seconds time lag for high energy bands in the light curve, which may be related with the delayed onset seen in GRB080916C.

The variability of optical \feii emission in PG QSO 1700+518

It is found that \feii emission contributes significantly to the optical and ultraviolet spectra of most active galactic nuclei. The origin of the optical/UV \feii emission is still a question open to debate. The variability of \feii would give clues to this origin. Using 7.5 yr spectroscopic monitoring data of one Palomer-Green (PG) quasi-stellar object (QSO), PG 1700+518, with strong optical \feii emission, we obtain the light curves of the continuum \lv, \feii, the broad component of \hb, and the narrow component of \hb by the spectral decomposition. Through the interpolation cross-correlation method, we calculate the time lags for light curves of \feii, the total \hb, the broad component of \hb, and the narrow component of \hb with respect to the continuum light curve. We find that the \feii time lag in PG1700+518 is $209^{+100}_{-147}$ days, and the \hb time lag cannot be determined. Assuming that \feii and \hb emission regions follow the virial relation between the time lag and the FWHM for the \hb and \feii emission lines, we can derive that the \hb time lag is $148^{+72}_{-104}$ days. The \hb time lag calculated from the empirical luminosity–size relation is 222 days, which is consistent with our measured \feii time lag. Considering the optical \feii contribution, PG 1700+518 shares the same characteristic on the spectral slope variability as other 15 PG QSOs in our previous work, i.e., harder spectrum during brighter phase.

The Lick AGN Monitoring Project: Reverberation Mapping of Optical Hydrogen and Helium Recombination Lines

We have recently completed a 64-night spectroscopic monitoring campaign at the Lick Observatory 3-m Shane telescope with the aim of measuring the masses of the black holes in 12 nearby (z < 0.05) Seyfert 1 galaxies with expected masses in the range ~10^6-10^7M_sun and also the well-studied nearby active galactic nucleus (AGN) NGC 5548. Nine of the objects in the sample (including NGC 5548) showed optical variability of sufficient strength during the monitoring campaign to allow for a time lag to be measured between the continuum fluctuations and the response to these fluctuations in the broad Hbeta emission, which we have previously reported. We present here the light curves for the Halpha, Hgamma, HeII 4686, and HeI 5876 emission lines and the time lags for the emission-line responses relative to changes in the continuum flux. Combining each emission-line time lag with the measured width of the line in the variable part of the spectrum, we determine a virial mass of the central supermassive black hole from several independent emission lines. We find that the masses are generally consistent within the uncertainties. The time-lag response as a function of velocity across the Balmer line profiles is examined for six of the AGNs. Finally we compare several trends seen in the dataset against the predictions from photoionization calculations as presented by Korista & Goad. We confirm several of their predictions, including an increase in responsivity and a decrease in the mean time lag as the excitation and ionization level for the species increases. Further confirmation of photoionization predictions for broad-line gas behavior will require additional monitoring programs for these AGNs while they are in different luminosity states. [abridged]

Discovery of a high confidence soft lag from an X-ray flare of Markarian 421

We present the X-ray variability properties of the X-ray and TeV bright blazar Mrk 421 with a ~60 ks long XMM-Newton observation performed on November 9-10, 2005. The source experienced a pronounced flare, of which the inter-band time lags were determined with a very high confidence level. The soft (0.6-0.8 keV) X-ray variations lagged the hard (4-10 keV) ones by 1.09 ks, and the soft lag increases with increasing difference in the photon energy. The energy-dependent soft lags can be well fitted with the difference of the energy-dependent cooling timescales of the relativistic electron distribution responsible for the observed X-ray emission, which constrains the magnetic field strength and Doppler factor of the emitting region to be B\delta^(1/3)~1.78 Gauss.

The Correlated Multi-color Optical Variations of BL Lac Object S5 0716+714

S5 0716+714 is a well-studied BL Lac object in the sky. Verifying the existence of correlations among the flux variations in different bands serves as an important tool to investigate the emission processes. To examine the possible existence of a lag between variations in different optical bands on this source, we employ a discrete correlation function (DCF) analysis on the light curves. In order to obtain statistically meaningful values for the cross-correlation time lags and their related uncertainties, we perform Monte Carlo simulations called "flux redistribution/random subset selection" (FR/RSS). Our analysis confirms that the variations in different optical light curves are strongly correlated. The time lags show a hint of the variations in high frequency band leading those in low frequency band of the order of a few minutes.

Solar dynamo and geomagnetic activity [Replacement]

The correlation between geomagnetic activity and the sunspot number in the 11-year solar cycle exhibits long-term variations due to the varying time lag between the sunspot-related and non-sunspot related geomagnetic activity, and the varying relative amplitude of the respective geomagnetic activity peaks. As the sunspot-related and non-sunspot related geomagnetic activity are caused by different solar agents, related to the solar toroidal and poloidal fields, respectively, we use their variations to derive the parameters of the solar dynamo transforming the poloidal field into toroidal field and back. We find that in the last 12 cycles the solar surface meridional circulation varied between 5 and 20 m/s (averaged over latitude and over the sunspot cycle), the deep circulation varied between 2.5 and 5.5 m/s, and the diffusivity in the whole of the convection zone was ~10**8 m2/s. In the last 12 cycles solar dynamo has been operating in moderately diffusion dominated regime in the bulk of the convection zone. This means that a part of the poloidal field generated at the surface is advected by the meridional circulation all the way to the poles, down to the tachocline and equatorward to sunspot latitudes, while another part is diffused directly to the tachocline at midlatitudes, "short-circuiting" the meridional circulation. The sunspot maximum is the superposition of the two surges of toroidal field generated by these two parts of the poloidal field, which is the explanation of the double peaks and the Gnevyshev gap in sunspot maximum. Near the tachocline, dynamo has been operating in diffusion dominated regime in which diffusion is more important than advection, so with increasing speed of the deep circulation the time for diffusive decay of the poloidal field decreases, and more toroidal field is generated leading to a higher sunspot maximum. During the Maunder minimum the dynamo was operating in advection dominated regime near the tachocline, with the transition from diffusion dominated to advection dominated regime caused by a sharp drop in the surface meridional circulation which is in general the most important factor modulating the amplitude of the sunspot cycle.

Solar dynamo and geomagnetic activity

The correlation between geomagnetic activity and the sunspot number in the 11-year solar cycle exhibits long-term variations due to the varying time lag between the sunspot-related and non-sunspot related geomagnetic activity, and the varying relative amplitude of the respective geomagnetic activity peaks. As the sunspot-related and non-sunspot related geomagnetic activity are caused by different solar agents, related to the solar toroidal and poloidal fields, respectively, we use their variations to derive the parameters of the solar dynamo transforming the poloidal field into toroidal field and back. We find that in the last 12 cycles the solar surface meridional circulation varied between 5 and 20 m/s (averaged over latitude and over the sunspot cycle), the deep circulation varied between 2.5 and 5.5 m/s, and the diffusivity in the whole of the convection zone was ~10**12 m2/s. In the last 12 cycles solar dynamo has been operating in moderately diffusion dominated regime in the bulk of the convection zone. This means that a part of the poloidal field generated at the surface is advected by the meridional circulation all the way to the poles, down to the tachocline and equatorward to sunspot latitudes, while another part is diffused directly to the tachocline at midlatitudes, "short-circuiting" the meridional circulation. The sunspot maximum is the superposition of the two surges of toroidal field generated by these two parts of the poloidal field, which is the explanation of the double peaks and the Gnevyshev gap in sunspot maximum. Near the tachocline, dynamo has been operating in diffusion dominated regime in which diffusion is more important than advection, so with increasing speed of the deep circulation the time for diffusive decay of the poloidal field decreases, and more toroidal field is generated leading to a higher sunspot maximum. During the Maunder minimum the dynamo was operating in advection dominated regime near the tachocline, with the transition from diffusion dominated to advection dominated regime caused by a sharp drop in the surface meridional circulation which is in general the most important factor modulating the amplitude of the sunspot cycle.

Spectral optical monitoring of 3C390.3 in 1995-2007: I. Light curves and flux variation of the continuum and broad lines [Replacement]

Here we present the results of the long-term (1995-2007) spectral monitoring of the broad line radio galaxy \object{3C~390.3}, a well known AGN with the double peaked broad emission lines, usually assumed to be emitted from an accretion disk. To explore dimensions and structure of the BLR, we analyze the light curves of the broad H$\alpha$ and H$\beta$ line fluxes and the continuum flux. In order to find changes in the BLR, we analyze the H$\alpha$ and H$\beta$ line profiles, as well as the change in the line profiles during the monitoring period. First we try to find a periodicity in the continuum and H$\beta$ light curves, finding that there is a good chance for quasi-periodical oscillations. Using the line shapes and their characteristics (as e.g. peaks separation and their intensity ratio, or FWHM) of broad H$\beta$ and H$\alpha$ lines, we discuss the structure of the BLR. Also, we cross-correlate the continuum flux with H$\beta$ and H$\alpha$ lines to find dimensions of the BLR. We found that during the monitoring period the broad emission component of the H$\alpha$ and H$\beta$ lines, and the continuum flux varied by a factor of $\approx $4-5. Also, we detected different structure in the line profiles of H$\alpha$ and H$\beta$. It seems that an additional central component is present and superposed to the disk emission. In the period of high activity (after 2002), H$\beta$ became broader than H$\alpha$ and red wing of H$\beta$ was higher than the one of H$\alpha$. We found time lags of $\sim$95 days between the continuum and H$\beta$ flux, and about 120 days between the continuum and H$\alpha$ flux. Variation in the line profiles, as well as correlation between the line and continuum flux during the monitoring period is in the favor of the disk origin of the broad lines with the possible contribution of some additional region and/or some kind of perturbation in the disk.

Spectral optical monitoring of 3C390.3 in 1995-2007: I. Light curves and flux variation of the continuum and broad lines

Here we present the results of the long-term (1995-2007) spectral monitoring of the broad line radio galaxy \object{3C~390.3}, a well known AGN with the double peaked broad emission lines, usually assumed to be emitted from an accretion disk. To explore dimensions and structure of the BLR, we analyze the light curves of the broad H$\alpha$ and H$\beta$ line fluxes and the continuum flux. In order to find changes in the BLR, we analyze the H$\alpha$ and H$\beta$ line profiles, as well as the change in the line profiles during the monitoring period. First we try to find a periodicity in the continuum and H$\beta$ light curves, finding that there is a good chance for quasi-periodical oscillations. Using the line shapes and their characteristics (as e.g. peaks separation and their intensity ratio, or FWHM) of broad H$\beta$ and H$\alpha$ lines, we discuss the structure of the BLR. Also, we cross-correlate the continuum flux with H$\beta$ and H$\alpha$ lines to find dimensions of the BLR. We found that during the monitoring period the broad emission component of the H$\alpha$ and H$\beta$ lines, and the continuum flux varied by a factor of $\approx $4-5. Also, we detected different structure in the line profiles of H$\alpha$ and H$\beta$. It seems that an additional central component is present and superposed to the disk emission. In the period of high activity (after 2002), H$\beta$ became broader than H$\alpha$ and red wing of H$\beta$ was higher than the one of H$\alpha$. We found time lags of $\sim$95 days between the continuum and H$\beta$ flux, and about 120 days between the continuum and H$\alpha$ flux. Variation in the line profiles, as well as correlation between the line and continuum flux during the monitoring period is in the favor of the disk origin of the broad lines with the possible contribution of some additional region and/or some kind of perturbation in the disk.

The hard state of black hole candidates: XTE J1752-223

We present a two-day long RXTE observation and simultaneous Swift data of the bright X-ray transient XTE J1752-223. Spectral and timing properties were stable during the observation. The energy spectrum is well described by a broken power-law with a high energy cut-off. A cold disc (~ 0.3 keV) is observed when Swift/XRT data are considered. The fractional rms amplitude of the aperiodic variability (0.002-128 Hz) is 48.2 +- 0.1% and it is not energy dependent. The continuum of the power density spectrum can be fitted by using four broad-band Lorentzians. A high frequency (~ 21 Hz) component and two weak QPO-like features are also present. Time-lags between soft and hard X-rays roughly follow the relation time-lag~vnu^ {-0.7}, with delays dropping from ~ 0.5 (0.003 Hz) to ~ 0.0015 (>10 Hz) seconds. Our results are consistent with XTE J1752-223 being a black-hole candidate, with all timing and spectral components very similar to those of Cyg X-1 during its canonical hard state.

Accretion and ejection in Sgr A*

We review our current understanding to the accretion and ejection processes in Sgr A*. Roughly speaking, they correspond to the quiescent and flare states of the source respectively. The high-resolution {\it Chandra} observations to the gas at the Bondi radius combined with the Bondi accretion theory, the spectral energy distribution from radio to X-ray, and the radio polarization provide us strict constraints and abundant information to the theory of accretion. We review these observational results and describe how the advection-dominated accretion flow model explains these observations. Recently more attentions have been paid to flares in Sgr A*. Many simultaneous multi-wavelength campaigns have been conducted, aiming at uncovering the nature of flares. The main observational properties of flares are briefly reviewed. Especially, the time lag between the peaks of flare at two radio frequencies strongly indicates that the flare is associated with ejection of radio-emitting blobs from the underlying accretion flow. Such kind of episodic jets is distinctive from the continuous jets and are quite common in black hole systems. We introduce the magnetohydrodynamical model for the formation of episodic jets recently proposed based on the analogy with the theory of coronal mass ejection in the Sun. We point out that the various observational appearances of flares should be explained in the framework of this model, since ejection and flare originate from the same physical process.

Tidal effects on brown dwarfs: Application to the eclipsing binary 2MASSJ05352184-0546085 - The anomalous temperature reversal in the context of tidal heating

2MASSJ05352184-0546085 (2M0535-05) is the only known eclipsing brown dwarf (BD) binary, and so may serve as an important benchmark for models of BD formation and evolution. However, theoretical predictions of the system’s properties seem inconsistent with observations: i. The more massive (primary) component is observed to be cooler than the less massive (secondary) one. ii. The secondary is more luminous (by roughly 10^{24} W) than expected. We study the impact of tidal heating to the energy budget of both components. We also compare various plausible tidal models to determine a range of predicted properties. We apply two versions of two different, well-known models for tidal interaction, respectively, (i.) the ‘constant-phase-lag’ model and (ii.) the ‘constant-time-lag’ model, and incorporate the predicted tidal heating into a model of BD structure. We find that the contribution of heat from tides in 2M0535-05 alone may only be large enough to account for the discrepancies between observation and theory in an unlikely region of the parameter space. The tidal quality factor of BDs, Q_{BD}, would have to be 10^{3.5} and the secondary needs a spin-orbit misalignment greater than 50 degrees. However, tidal synchronization time scales for 2M0535-05 restrict the tidal dissipation function Q_{BD} to values greater than 10^{4.5} and rule out intense tidal heating in 2M0535-05. We provide the first constraint on Q_{BD}. Tidal heating alone is unlikely to be responsible for the surprising temperature reversal within 2M0535-05. But an evolutionary embedment of tidal effects and a coupled treatment with the structural evolution of the BDs is necessary to corroborate or refute this result.

Multiwavelength observations of 3C 454.3. III. Eighteen months of AGILE monitoring of the "Crazy Diamond"

We report on 18 months of multiwavelength observations of the blazar 3C 454.3 (Crazy Diamond) carried out in July 2007-January 2009. We show the results of the AGILE campaigns which took place on May-June 2008, July-August 2008, and October 2008-January 2009. During the May 2008-January 2009 period, the source average flux was highly variable, from an average gamma-ray flux F(E>100MeV) > 200E-8 ph/cm2/s in May-June 2008, to F(E>100MeV)~80E-8 ph/cm2/s in October 2008-January 2009. The average gamma-ray spectrum between 100 MeV and 1 GeV can be fit by a simple power law (Gamma_GRID ~ 2.0 to 2.2). Only 3-sigma upper limits can be derived in the 20-60 keV energy band with Super-AGILE. During July-August 2007 and May-June 2008, RXTE measured a flux of F(3-20 keV)= 8.4E-11 erg/cm2/s, and F(3-20 keV)=4.5E-11 erg/cm2/s, respectively and a constant photon index Gamma_PCA=1.65. Swift/XRT observations were carried out during all AGILE campaigns, obtaining a F(2-10 keV)=(0.9-7.5)E-11 erg/cm2/s and a photon index Gamma_XRT=1.33-2.04. BAT measured an average flux of ~5 mCrab. GASP-WEBT monitored 3C 454.3 during the whole 2007-2008 period from the radio to the optical. A correlation analysis between the optical and the gamma-ray fluxes shows a time lag of tau=-0.4 days. An analysis of 15 GHz and 43 GHz VLBI core radio flux observations shows an increasing trend of the core radio flux, anti- correlated with the higher frequency data. The modeling SEDs, and the behavior of the long-term light curves in different energy bands, allow us to compare the jet properties during different emission states, and to study the geometrical properties of the jet on a time-span longer than one year.

Suzaku and Optical Spectroscopic Observations of SS 433 in the 2006 April Multiwavelength Campaign

We report results of the 2006 April multi-wavelengths campaign of SS 433, focusing on X-ray data observed with Suzaku at two orbital phases (in- and out-of- eclipse) and simultaneous optical spectroscopic observations. By analyzing the Fe25 K_alpha lines originating from the jets, we detect rapid variability of the Doppler shifts, dz/dt ~ 0.019/0.33 day^-1, which is larger than those expected from the precession and/or nodding motion. This phenomenon probably corresponding to "jitter" motions observed for the first time in X-rays, for which significant variability both in the jet angle and intrinsic speed is required. From the time lag of optical Doppler curves from those of X-rays, we estimate the distance of the optical jets from the base to be ~(3-4) \times 10^14 cm. Based on the radiatively cooling jet model, we determine the innermost temperature of the jets to be T_0 = 13 +/- 2 keV and 16 +/- 3 keV (the average of the blue and red jets) for the out-of-eclipse and in-eclipse phase, respectively, from the line intensity ratio of Fe25 K_alpha and Fe26 K_alpha. While the broad band continuum spectra over the 5–40 keV band in eclipse is consistent with a multi-temperature bremsstrahlung emission expected from the jets, and its reflection component from cold matter, the out-of-eclipse spectrum is harder than the jet emission with the base temperature determined above, implying the presence of an additional hard component.

Timing Signatures of the Internal-Shock Model for Blazars

We investigate the spectral and timing signatures of the internal-shock model for blazars. For this purpose, we develop a semi-analytical model for the time-dependent radiative output from internal shocks arising from colliding relativistic shells in a blazar jet. The emission through synchrotron and synchrotron-self Compton (SSC) radiation as well as Comptonization of an isotropic external radiation field are taken into account. We evaluate the discrete correlation function (DCF) of the model light curves in order to evaluate features of photon-energy dependent time lags and the quality of the correlation, represented by the peak value of the DCF. The almost completely analytic nature of our approach allows us to study in detail the influence of various model parameters on the resulting spectral and timing features. This paper focuses on a range of parameters in which the gamma-ray production is dominated by Comptonization of external radiation, most likely appropriate for gamma-ray bright flat-spectrum radio quasars (FSRQs) or low-frequency peaked BL Lac objects (LBLs). In most cases relevant for FSRQs and LBLs, the variability of the optical emission is highly correlated with the X-ray and high-energy (HE: > 100 MeV) gamma-ray emission. Our baseline model predicts a lead of the optical variability with respect to the higher-energy bands by 1 – 2 hours and of the HE gamma-rays before the X-rays by about 1 hour. We show that variations of certain parameters may lead to changing signs of inter-band time lags, potentially explaining the lack of persistent trends of time lags in most blazars.

A useful guide for gravitational wave observers to test modified gravity models [Replacement]

We present an extension of a previously suggested test of all modified theories of gravity that would reproduce MOND at low accelerations. In a class of models, called “dark matter emulators”, gravitational waves and other particles couple to different metrics. This leads to a detectable time lag between their detection at Earth from the same source. We calculate this time lag numerically for any event that occurs in our galaxy up to 400 kpc, and present a graph of this possible time lag. This suggests that, gravitational wave observers might have to consider the possibility of extending their analysis to non-coincident gravitational and electromagnetic signals, and the graph that we present might be a useful guideline for this effort.

A useful guide for gravitational wave observers to test modified gravity models [Replacement]

We present an extension of a previously suggested test of all modified theories of gravity that would reproduce MOND at low accelerations. In a class of models, called "dark matter emulators", gravitational waves and other particles couple to different metrics. This leads to a detectable time lag between their detection at Earth from the same source. We calculate this time lag numerically for any event that occurs in our galaxy up to 400 kpc, and present a graph of this possible time lag. This suggests that, gravitational wave observers might have to consider the possibility of extending their analysis to non-coincident gravitational and electromagnetic signals, and the graph that we present might be a useful guideline for this effort.

A useful guide for gravitational wave observers to test modified gravity models [Cross-Listing]

We present an extension of a previously suggested test of all modified theories of gravity that would reproduce MOND at low accelerations. In a class of models, called "dark matter emulators", gravitational waves and other particles couple to different metrics. This leads to a detectable time lag between their detection at Earth from the same source. We calculate this time lag numerically for any event that occurs in our galaxy up to 400 kpc, and present a graph of this possible time lag. This suggests that, gravitational wave observers might have to consider the possibility of extending their analysis to non-coincident gravitational and electromagnetic signals, and the graph that we present might be a useful guideline for this effort.

Limits on Spacetime Foam [Replacement]

Plausibly spacetime is “foamy” on small distance scales, due to quantum fluctuations. We elaborate on the proposal to detect spacetime foam by looking for seeing disks in the images of distant quasars and AGNs. This is a null test in the sense that the continued presence of unresolved “point” sources at the milli-arc second level in samples of distant compact sources puts severe constraints on theories of quantized spacetime foam at the Planckian level. We discuss the geometry of foamy spacetime, and the appropriate distance measure for calculating the expected angular broadening. We then deal with recent data and the constraints they put on spacetime foam models. While time lags from distant pulsed sources such as GRBs have been posited as a possible test of spacetime foam models, we demonstrate that the time-lag effect is rather smaller than has been calculated, due to the equal probability of positive and negative fluctuations in the speed of light inherent in such models. Thus far, images of high-redshift quasars from the Hubble Ultra-Deep Field (UDF) provide the most stringent test of spacetime foam theories. While random walk models ($\alpha = 1/2$) have already been ruled out, the holographic ($\alpha=2/3$) model remains viable. Here $\alpha \sim 1$ parametrizes the different spacetime foam models according to which the fluctuation of a distance $l$ is given by $\sim l^{1 – \alpha} l_P^{\alpha}$ with $l_P$ being the Planck length. Indeed, we see a slight wavelength-dependent blurring in the UDF images selected for this study. Using existing data in the {\it Hubble Space Telescope (HST)} archive we find it is impossible to rule out the $\alpha=2/3$ model, but exclude all models with $\alpha<0.65$. By comparison, current GRB time lag observations only exclude models with $\alpha<0.3$.

Limits on Spacetime Foam [Replacement]

Plausibly spacetime is “foamy” on small distance scales, due to quantum fluctuations. We elaborate on the proposal to detect spacetime foam by looking for seeing disks in the images of distant quasars and AGNs. This is a null test in the sense that the continued presence of unresolved “point” sources at the milli-arc second level in samples of distant compact sources puts severe constraints on theories of quantized spacetime foam at the Planckian level. We discuss the geometry of foamy spacetime, and the appropriate distance measure for calculating the expected angular broadening. We then deal with recent data and the constraints they put on spacetime foam models. While time lags from distant pulsed sources such as GRBs have been posited as a possible test of spacetime foam models, we demonstrate that the time-lag effect is rather smaller than has been calculated, due to the equal probability of positive and negative fluctuations in the speed of light inherent in such models. Thus far, images of high-redshift quasars from the Hubble Ultra-Deep Field (UDF) provide the most stringent test of spacetime foam theories. While random walk models ($\alpha = 1/2$) have already been ruled out, the holographic ($\alpha=2/3$) model remains viable. Here $\alpha \sim 1$ parametrizes the different spacetime foam models according to which the fluctuation of a distance $l$ is given by $\sim l^{1 – \alpha} l_P^{\alpha}$ with $l_P$ being the Planck length. Indeed, we see a slight wavelength-dependent blurring in the UDF images selected for this study. Using existing data in the {\it Hubble Space Telescope (HST)} archive we find it is impossible to rule out the $\alpha=2/3$ model, but exclude all models with $\alpha<0.65$. By comparison, current GRB time lag observations only exclude models with $\alpha<0.3$.

Absolute timing of the Crab pulsar with the INTEGRAL/SPI telescope

We have investigated the pulse shape evolution of the Crab pulsar emission in the hard X-ray domain of the electromagnetic spectrum. In particular, we have studied the alignment of the Crab pulsar phase profiles measured in the hard X-rays and in other wavebands. To obtain the hard X-ray pulse profiles, we have used six year (2003-2009, with a total exposure of about 4 Ms) of publicly available data of the SPI telescope on-board of the INTEGRAL observatory, folded with the pulsar time solution derived from the Jodrell Bank Crab Pulsar Monthly Ephemeris. We found that the main pulse in the hard X-ray 20-100 keV energy band is leading the radio one by $8.18\pm0.46$ milliperiods in phase, or $275\pm15 \mu s$ in time. Quoted errors represent only statistical uncertainties.Our systematic error is estimated to be $\sim 40 \mu s$ and is mainly caused by the radio measurement uncertainties. In hard X-rays, the average distance between the main pulse and interpulse on the phase plane is $0.3989\pm0.0009$. To compare our findings in hard X-rays with the soft 2-20 keV X-ray band, we have used data of quasi-simultaneous Crab observations with the PCA monitor on-board the Rossi X-Ray Timing Explorer (RXTE) mission. The time lag and the pulses separation values measured in the 3-20 keV band are $0.00933\pm0.00016$ (corresponding to $310\pm6 \mu s$) and $0.40016\pm0.00028$ parts of the cycle, respectively. While the pulse separation values measured in soft X-rays and hard X-rays agree, the time lags are statistically different. Additional analysis show that the delay between the radio and X-ray signals varies with energy in the 2 – 300 keV energy range. We explain such a behaviour as due to the superposition of two independent components responsible for the Crab pulsed emission in this energy band.

A time-dependent jet model for the emission from Sagittarius A*

The source of emission from Sgr A*, the supermassive black hole at the Galactic Center, is still unknown. Flares and data from multiwavelength campaigns provide important clues about the nature of Sgr A* itself. Here we attempt to constrain the physical origin of the broadband emission and the radio flares from Sgr A*. We developed a time-dependent jet model, which for the first time allows one to compare the model predictions with flare data from Sgr A*. Taking into account relevant cooling mechanisms, we calculate the frequency-dependent time lags and photosphere size expected in the jet model. The predicted lags and sizes are then compared with recent observations. Both the observed time lags and size-frequency relationships are reproduced well by the model. The combined timing and structural information strongly constrain the speed of the outflow to be mildly relativistic, and the radio flares are likely to be caused by a transient increase in the matter channelled into the jets. The model also predicts light curves and structural information at other wavelengths which could be tested by observations in the near future. We show that a time-dependent relativistic jet model can successfully reproduce: (1) the quiescent broadband spectral energy distribution of Sgr A*, (2) the observed 22 and 43 GHz light curve morphologies and time lags, and (3) the frequency-size relationship. The results suggest that the observed emission at radio frequencies from Sgr A* is most easily explained by a stratified, optically thick, mildly relativistic jet outflow. Frequency-dependent measurements of time-lags and intrinsic source size provide strong constraints on the bulk motion of the jet plasma.

The optical microvariability and spectral changes of the BL Lacertae object S5 0716+714

We monitored the BL Lac object S5 0716+714 in the optical band during October 2008, December 2008 and February 2009 with a best temporal resolution of about 5 minutes in the BVRI bands. Four fast flares were observed with amplitudes ranging from 0.3 to 0.75 mag. The source remained active during the whole monitoring campaign, showing microvariability in all days except for one. The overall variability amplitudes are delta B ~ 0.89 mag, delta V ~ 0.80 mag, delta R ~ 0.73 mag and delta I ~0.51 mag. Typical timescales of microvariability range from 2 to 8 hours. The overall V – R color index ranges from 0.37 to 0.59. Strong bluer- when-brighter chromatism was found on internight timescales. However, different spectral behavior was found on intranight timescales. A possible time lag of ~ 11 mins between B and I bands was found on one night. The shock-in-jet model and geometric effects can be applied to explain the source’s intranight behavior.

Long-term variability of radio-bright BL Lacertae objects

Radio-bright BL Lacertae objects (BLOs) are typically variable and exhibit prominent flaring. We use a sample of 24 BLOs to get a clear idea of their flaring behavior and to find possible commonalities in their variability patterns. Our goal was to compare the results given by computational time scales and the observed variability parameters determined directly from the flux curves. Also, we wanted to find out if the BLO flares adhere to the generalized shock model. We use long-term monitoring data from 4.8, 8, 14.5, 22, 37, 90 and 230 GHz. The structure function, discrete correlation function and Lomb-Scargle periodogram time scales, calculated in a previous study, are analyzed in more detail. We determine flare durations, rise and decay times, absolute and relative peak fluxes from the monitoring data. We find that BLOs demonstrate a wide range of variability behavior. BLOs include sources with fast and strong variability, such as OJ 287, PKS 1749+096 and BL Lac, but also sources with more rolling fluctuations like PKS 0735+178. The most extreme flares can last for up to 13 years or have peak fluxes of approximately 12 Jy in the observer’s frame. When the Doppler boosting effect is taken into account, the peak flux of a flare does not depend on the duration of the flare. A rough analysis of the time lags and peak flux evolution indicates that BLO flares in the mm – cm wavelengths are high-peaking, i.e., are in the adiabatic stage. Thus, the results concur with the generalized shock model.

The Lick AGN Monitoring Project: Photometric Light Curves and Optical Variability Characteristics

The Lick AGN Monitoring Project targeted 13 nearby Seyfert 1 galaxies with the intent of measuring the masses of their central black holes using reverberation mapping. The sample includes 12 galaxies selected to have black holes with masses roughly in the range 10^6-10^7 solar masses, as well as the well-studied AGN NGC 5548. In conjunction with a spectroscopic monitoring campaign, we obtained broad-band B and V images on most nights from 2008 February through 2008 May. The imaging observations were carried out by four telescopes: the 0.76-m Katzman Automatic Imaging Telescope (KAIT), the 2-m Multicolor Active Galactic Nuclei Monitoring (MAGNUM) telescope, the Palomar 60-in (1.5-m) telescope, and the 0.80-m Tenagra II telescope. Having well-sampled light curves over the course of a few months is useful for obtaining the broad-line reverberation lag and black hole mass, and also allows us to examine the characteristics of the continuum variability. In this paper, we discuss the observational methods and the photometric measurements, and present the AGN continuum light curves. We measure various variability characteristics of each of the light curves. We do not detect any evidence for a time lag between the B- and V-band variations, and we do not find significant color variations for the AGNs in our sample.

Simultaneous MITSuME gRI monitoring of S5 0716+714

We present results of our intra-night optical flux monitoring observations of S5 0716+714 done simultaneously in gRI filters. The observations were done using Multicolor Imaging Telescopes for Survey and Monstrous Explosions (MITSuME) instrument on the 50 cm telescope at the Okayama Astrophysical Observatory over 30 nights between 11 March 2008 and 8 May 2008. Of these 30 nights, 22 nights have continuous (without any break) observations with duration ranging from 1 to 6 hours and hence were considered for intra-night optical variability (INOV). In total we have 4888 datapoints which were simultaneous in gR and I filters. Of the 22 nights considered for INOV, the object showed flux variability on 19 nights with the amplitude of variability in the I-band ranging from ~4% to ~55%. The duty cycle for INOV was thus found to be 83%. No time lag between different bands was noticed on most of the nights, except for 3 nights where the variation in g was found to lead that of the I band by 0.3 to 1.5 hrs. On inter-night timescales, no lag was found between g and I bands. On inter-night timescales as well as intra-night timescales on most of the nights, the amplitude of variability was found to increase toward shorter wavelengths. The flux variations in the different bands were not achromatic, with the blazar tending to become bluer when brighter both on inter-night and intra-night timescales; and this might be attributed to the larger amplitude variation at shorter wavelengths. A clear periodic variation of 3.3 hrs was found on 1 April 2008 and a hint for another possible periodic variability of 4 hrs was found on 31 March 2008. During our 30 days of observations over a 2 month period the source has varied with an amplitude of variability as large as ~80%.

Energy Dependent time lags in the Seyfert 1 galaxy NGC 4593

We investigate the energy-time lag dependence of the source NGC 4593 using XMM-{\it Newton}/EPIC-pn data. We found that the time lag dependency is linear in nature with respect to the logarithm of different energy bands. We also investigate the frequency dependent time lags and identify that at some frequency range (5 $\times$ 10$^{-5}$ Hz — 2 $\times$ 10$^{-4}$ Hz) the X-ray emission is highly coherent, mildly frequency dependent and very strongly energy dependent. These observations can be explained in the frame work of the thermal Comptonization process and they indicate a truncated accretion disk very close to the black hole. We discuss the plausible spectral state to explain the phenomenon and conclude that the observed properties bear a close resemblance to the intermediate state or the steep power-law state, found in galactic black hole sources.

The correlated optical and radio variability of BL Lacertae. WEBT data analysis 1994-2005

Since 1997, BL Lacertae has undergone a phase of high optical activity, with the occurrence of several prominent outbursts. Starting from 1999, the Whole Earth Blazar Telescope (WEBT) consortium has organized various multifrequency campaigns on this blazar, collecting tens of thousands of data points. One of the main issues in the study of this huge dataset has been the search for correlations between the optical and radio flux variations, and for possible periodicities in the light curves. The analysis of the data assembled during the first four campaigns (comprising also archival data to cover the period 1968-2003) revealed a fair optical-radio correlation in 1994-2003, with a delay of the hard radio events of ~100 days. Moreover, various statistical methods suggested the existence of a radio periodicity of ~8 years. In 2004 the WEBT started a new campaign to extend the dataset to the most recent observing seasons, in order to possibly confirm and better understand the previous results. In this campaign we have collected and assembled about 11000 new optical observations from twenty telescopes, plus near-IR and radio data at various frequencies. Here, we perform a correlation analysis on the long-term R-band and radio light curves. In general, we confirm the ~100-day delay of the hard radio events with respect to the optical ones, even if longer (~200-300 days) time lags are also found in particular periods. The radio quasi-periodicity is confirmed too, but the "period" seems to progressively lengthen from 7.4 to 9.3 years in the last three cycles. The optical and radio behaviour in the last forty years suggests a scenario where geometric effects play a major role. In particular, the alternation of enhanced and suppressed optical activity (accompanied by hard and soft radio events, respectively) can

Blazar Optical Variability in the Palomar-QUEST Survey

We study the ensemble optical variability of 276 FSRQs and 86 BL Lacs in the Palomar-QUEST Survey with the goal of searching for common fluctuation properties, examining the range of behavior across the sample, and characterizing the appearance of blazars in such a survey so that future work can more easily identify such objects. The survey, which covers 15,000 square degrees multiple times over 3.5 years, allows for the first ensemble blazar study of this scale. Variability amplitude distributions are shown for the FSRQ and BL Lac samples for numerous time lags, and also studied through structure function analyses. Individual blazars show a wide range of variability amplitudes, timescales, and duty cycles. Of the best sampled objects, 35% are seen to vary by more than 0.4 magnitudes; for these, the fraction of measurements contributing to the high amplitude variability ranges constantly from about 5% to 80%. Blazar variability has some similarities to that of type I quasars but includes larger amplitude fluctuations on all timescales. FSRQ variability amplitudes are particularly similar to those of QSOs on timescales of several months, suggesting significant contributions from the accretion disk to the variable flux at these timescales. Optical variability amplitudes are correlated with the maximum apparent velocities of the radio jet for the subset of FSRQs with MOJAVE VLBA measurements, implying that the optically variable flux’s strength is typically related to that of the radio emission. We also study CRATES radio-selected FSRQ candidates, which show similar variability characteristics to known FSRQs; this suggests a high purity for the CRATES sample.

Correlated optical and gamma emissions from GRB 081126

We present an analysis of time-resolved optical emissions observed from the gamma-ray burst GRB 081126 during the prompt phase. The analysis employed time-resolved photometry using optical data obtained by the TAROT telescope, using BAT data from the Swift spacecraft, and time-resolved spectroscopy at high energies from the GBM instrument onboard the Fermi spacecraft. The optical emission of GRB 081126 is found to be compatible with the second gamma emission pulse shifted by a positive time lag of 8.4 $\pm$ 3.9 s. This is the first well-resolved observation of a time lag between optical and gamma emissions during a gamma-ray burst. Our observations could potentially provide new constraints on the fireball model for gamma-ray burst early emissions. Furthermore, observations of time lags between optical and gamma ray photons provides an exciting opportunity to constrain quantum gravity theories.

Cosmological perturbations from an inhomogeneous phase transition [Cross-Listing]

A mechanism for generating metric perturbations in inflationary models is considered. Long-wavelength inhomogeneities of light scalar fields in a decoupled sector may give rise to superhorizon fluctuations of couplings and masses in the low-energy effective action. Cosmological phase transitions may then occur that are not simultaneous in space, but occur with time lags in different Hubble patches that arise from the long-wavelength inhomogeneities. Here an interesting model in which cosmological perturbations may be created at the electroweak phase transition is considered. The results show that phase transitions may be a generic source of non-Gaussianity.

Probing a Possible Vacuum Refractive Index with Gamma-Ray Telescopes

We have used a stringy model of quantum space-time foam to suggest that the vacuum may exhibit a non-trivial refractive index depending linearly on gamma-ray energy: eta -1 ~ E_gamma/M_QG1, where M_QG1 is some mass scale typical of quantum gravity that may be ~ 10^18 GeV: see Phys. Lett. B 665, 412 (2008) and references therein. The MAGIC, HESS and Fermi gamma-ray telescopes have recently probed the possible existence of such an energy-dependent vacuum refractive index. All find indications of time-lags for higher-energy photons, but cannot exclude the possibility that they are due to intrinsic delays at the sources. However, the MAGIC and HESS observations of time-lags in emissions from AGNs Mkn 501 and PKS 2155-304 are compatible with each other and a refractive index depending linearly on the gamma-ray energy, with M_QG1 ~ 10^18 GeV. We combine their results to estimate the time-lag Delta t to be expected for the highest-energy photon from GRB 080916c measured by the Fermi telescope, which has an energy ~ 13.2 GeV, assuming the redshift z = 4.2 \pm 0.3 measured by GROND. In the case of a refractive index depending linearly on the gamma-ray energy we predict Delta t = 25 \pm 11 s. This is compatible with the time-lag Delta t <= 16.5 s reported by the Fermi Collaboration, whereas the time-lag would be negligible in the case of a refractive index depending quadratically on the gamma-ray energy. We suggest a strategy for future observations that could distinguish between a quantum-gravitational effect and other interpretations of the time-lags observed by the MAGIC, HESS and Fermi gamma-ray telescopes.

Jet-lag in Sgr A*: What size and timing measurements tell us about the central black hole in the Milky Way

The black hole at the Galactic Center, Sgr A*, is the prototype of a galactic nucleus at a very low level of activity. Its radio through submm-wave emission is known to come from a region close to the event horizon, however, the source of the emission is still under debate. A successful theory explaining the emission is based on a relativistic jet model scaled down from powerful quasars. We want to test the predictive power of this established jet model against newly available measurements of wavelength-dependent time lags and the size-wavelength structure in Sgr A*. Using all available closure amplitude VLBI data from different groups, we again derived the intrinsic wavelength-dependent size of Sgr A*. This allowed us to calculate the expected frequency-dependent time lags of radio flares, assuming a range of in- and outflow velocities. Moreover, we calculated the time lags expected in the previously published pressure-driven jet model. The predicted lags are then compared to radio monitoring observations at 22, 43, and 350 GHz. The combination of time lags and size measurements imply a mildly relativistic outflow with bulk outflow speeds of gamma*beta ~ 0.5-2. The newly measured time lags are reproduced well by the jet model without any major fine tuning. The results further strengthen the case for the cm-to-mm wave radio emission in Sgr A* as coming from a mildly relativistic jet-like outflow. The combination of radio time lag and VLBI closure amplitude measurements is a powerful new tool for assessing the flow speed and direction in Sgr A*. Future VLBI and time lag measurements over a range of wavelengths will reveal more information about Sgr A*, such as the existence of a jet nozzle, and measure the detailed velocity structure of a relativistic jet near its launching point for the first time.

Correlated variability in the blazar 3C 454.3

The blazar 3C 454.3 was revealed by the Fermi Gamma-ray Space Telescope to be in an exceptionally high flux state in July 2008. Accordingly, we performed a multi-wavelength monitoring campaign on this blazar using IR and optical observations from the SMARTS telescopes, optical, UV and X-ray data from the Swift satellite, and public-release gamma-ray data from Fermi. We find an excellent correlation between the IR, optical, UV and gamma-ray light curves, with a time lag of less than one day. The amplitude of the infrared variability is comparable to that in gamma-rays, and larger than at optical or UV wavelengths. The X-ray flux is not strongly correlated with either the gamma-rays or longer wavelength data. These variability characteristics find a natural explanation in the external Compton model, in which electrons with Lorentz factor gamma~10^(3-4) radiate synchrotron emission in the infrared-optical and also scatter accretion disk or emission line photons to gamma-ray energies, while much cooler electrons (gamma~10^(1-2)) produce X-rays by scattering synchrotron or other ambient photons.

Quasar optical variability: searching for interband time delays

Aims. The main purpose of this paper is to study time delays between the light variations in different wavebands for a sample of quasars. Measuring a reliable time delay for a large number of quasars may help constraint the models of their central engines. The standard accretion disk irradiation model predicts a delay of the longer wavelengths behind the shorter ones, a delay that depends on the fundamental quasar parameters. Since the black hole masses and the accretion rates are approximately known for the sample we use, one can compare the observed time delays with the expected ones. Methods. We applied the interpolation cross-correlation function (ICCF) method to the Giveon et al. sample of 42 quasars, monitored in two (B and R) colors, to find the time lags represented by the ICCF peaks. Different tests were performed to assess the influence of photometric errors, sampling, etc., on the final result. Results. We found that most of the objects show a delay in the red light curve behind the blue one (a positive lag), which on average for the sample is about +4 days (+3 for the median), although the scatter is significant. These results are broadly consistent with the reprocessing model, especially for the well-sampled objects. The normalized time-lag deviations do not seem to correlate significantly with other quasar properties, including optical, radio, or X-ray measurables. On the other hand, many objects show a clear negative lag, which, if real, may have important consequences for the variability models.

New two-colour light curves of Q0957+561: time delays and the origin of intrinsic variations

We extend the gr-band time coverage of the gravitationally lensed double quasar Q0957+561. New gr light curves permit us to detect significant intrinsic fluctuations, to determine new time delays, and thus to gain perspective on the mechanism of intrinsic variability in Q0957+561. We use new optical frames of Q0957+561 in the g and r passbands from January 2005 to July 2007. These frames are part of an ongoing long-term monitoring with the Liverpool robotic telescope. We also introduce two photometric pipelines that are applied to the new gr frames of Q0957+561. The transformation pipeline incorporates zero-point, colour, and inhomogeneity corrections to the instrumental magnitudes, so final photometry to the 1-2% level is achieved for both quasar components. The two-colour final records are then used to measure time delays. The gr light curves of Q0957+561 show several prominent events and gradients, and some of them (in the g band) lead to a time delay between components of 417 +/- 2 d (1 sigma). We do not find evidence of extrinsic variability in the light curves of Q0957+561. We also explore the possibility of a delay between a large event in the g band and the corresponding event in the r band. The gr cross-correlation reveals a time lag of 4.0 +/- 2.0 d (1 sigma; the g-band event is leading) that confirms a previous claim of the existence of a delay between the g and r band in this lensed quasar. The time delays (between quasar components and between optical bands) from the new records and previous ones in similar bands indicate that most observed variations in Q0957+561 (amplitudes of about 100 mmag and timescales of about 100 d) are very probably due to reverberation within the gas disc around the supermassive black hole.

Study on X-ray Spectra of Obscured AGNs based on Monte Carlo simulation - an interpretation of observed wide-band spectra

Monte Carlo simulation is one of the best tools to study the complex spectra of Compton-thick AGNs and to figure out the relation between their nuclear structures and X-ray spectra. We have simulated X-ray spectra of Compton-thick AGNs obscured by an accretion torus whose structure is characterized by a half-opening angle, an inclination angle of the torus relative to the observer, and a column density along the equatorial plane. We divided the simulated spectra into three components: one direct component, an absorbed reflection component and an unabsorbed reflection component. We then deduced the dependencies of these components on the parameters describing the structure of the torus. Our simulation results were applied to fit the wide-band spectrum of the Seyfert 2 galaxy Mrk 3 obtained by $Suzaku$. The spectral analysis indicates that we observe the nucleus along a line of sight intercepting the torus near its edge, and the column density along the equatorial plane was estimated to be ~10^24 cm^-2. Using this model, we can estimate the luminosities of both the direct emission and the emission irradiating the surrounding matter. This is useful to find the time variability and time lag between the direct and reflected light.

Limits on an Energy Dependence of the Speed of Light from a Flare of the Active Galaxy PKS 2155-304

In the past few decades, several models have predicted an energy-dependence of the speed of light in the context of quantum gravity. For cosmological sources such as active galaxies, this minuscule effect can add up to measurable photon-energy dependent time lags. In this paper a search for such time lags during the H.E.S.S. observations of the exceptional very high energy flare of the active galaxy PKS 2155-304 on 28 July in 2006 is presented. Since no significant time lag is found, lower limits on the energy scale of speed of light modifications are derived.

GOLF: a new proxy for solar magnetism

Solar magnetism is measured with different indexes: for instance the MPSI and the MWSI, number of sunspots, radio flux at 10.7 cm, Ca II K, Mg II K, EUV, He I or L_alpha. Bachmann & White (1994) had compared these indicators of the solar activity showing a hysteresis of the solar cycle variations and a time lag between these indices not related to instrumental effects. Later on, Ozguc & Atac (2001) extended this study of hysteresis phenomenon between Flare index and other solar indices (mean magnetic field, coronal index). In its original working configuration, GOLF/SoHO was able to measure during 26 days the solar mean magnetic field (Garcia et al. 1999). We check here if the velocity data could be used as another solar magnetism proxy with the advantage of having a duty cycle >95% during the last 12 years. We will compare the GOLF data with some of the above-mentioned solar activity indexes.

A Study of the X-Ray Dust Scattering Halo of Cyg X-1 with a Cross-Correlation Method

X-ray photons scattered by the interstellar medium carry information about dust distribution, dust grain model, scattering cross section, and the distance of the source; they also take longer time than unscattered photons to reach the observer. Using a cross-correlation method, we study the light curves of the X-ray dust scattering halo of Cyg X-1, observed with the \textit{Chandra X-ray Observatory}. Significant time lags are found between the light curves of the point source and its halo. This time lag increases with the angular distance from Cyg X-1, implying a dust concentration at a distance along the line of sight (LOS) of 2.0 kpc $\times$ (0.876 $\pm$ 0.002) from the Earth. By fitting the observed light curves of the halo at different radii with simulated light curves, we obtain a width of $\mathit{\Delta L}=33_{-13}^{+18}$ pc of this dust concentration. The origin of this dust concentration is still not clearly known. The advantage of our method is that we need no assumption of scattering cross section, dust grain model, or dust distribution along the LOS. Combining the derived dust distribution from the cross-correlation study with the surface brightness distribution of the halo, we conclude that the two commonly accepted models of dust grain size distribution need to be modified significantly.

 

You need to log in to vote

The blog owner requires users to be logged in to be able to vote for this post.

Alternatively, if you do not have an account yet you can create one here.

Powered by Vote It Up

^ Return to the top of page ^