We present near-infrared 1.1-1.3 and 1.3-1.6 $\mu$m spectra of the Ly$\alpha$-luminous quasar PSO J006.1240+39.2219 at $z=6.617$ obtained with the NIRSPEC spectrograph at the Keck-II telescope. The spectra cover the CIV $\lambda$1549, CIII] $\lambda$1909 emission lines and part of the UV continuum of the quasar. From the NIRSPEC observations of PSO J006.1240+39.2219, we constrain the spectral slope of its UV continuum to be $\alpha_{\lambda}=-1.35\pm0.26$ and measure an absolute magnitude of $M_{1450}=-25.60$. Using the scaling relation between black hole mass, width of the CIV line and ultraviolet continuum luminosity, we derive a black hole mass of $(2.19\pm0.30)\times 10^8 M_{sun}$, which is consistent but somewhat smaller than the typical black hole masses of $z\gtrsim6$ quasars of similar luminosities. The inferred accretion rate of $L_{bol}/L_{edd}\gtrsim2$ indicates that PSO J006.1240+39.2219 is in the phase of the rapid growth of its supermassive black hole characterized by the high NV/CIV line ratio, NV/CIV$>1$, and lower level of ionization of its circumnuclear gas than in other high-redshift luminous quasars. The NV/CIV line ratio of PSO J006.1240+39.2219 implies relatively high abundance of nitrogen in its circumnuclear gas. This abundance might be produced by the post-starburst population of stars that provide the fuel for black hole accretion.

We report the first detection of molecular gas beyond redshift 7. The molecular gas is associated with the host galaxy of the quasar candidate PSO J145.5964+19.3565 and its companion PSO J145.5964+19.3565N separated by 20.7 kpc. The molecular gas of both companions is detected in two rotational transition lines of carbon monoxide, CO(6-5) and CO(7-6). We also detected the Lyman-alpha line of PSO J145.5964+19.3565 at z=7.08. The Lyman-alpha line emission is extended and might represent the blended emission of two different sources at a separation of <5 kpc. The detected CO and Lyman-alpha emission likely originate from a system of interacting star-forming galaxies that might host a quasar(s). We also report the detection of a new emission line from the system which is a possible 793.62-GHz water maser line.

Distant luminous quasars provide important information on the growth of the first supermassive black holes, their host galaxies and the epoch of reionization. The identification of quasars is usually performed through detection of their Lyman-$\alpha$ line redshifted to $\sim$ 0.9 microns at z>6.5. Here, we report the discovery of a very Lyman-$\alpha$ luminous quasar, PSO J006.1240+39.2219 at redshift z=6.618, selected based on its red colour and multi-epoch detection of the Lyman-$\alpha$ emission in a single near-infrared band. The Lyman-$\alpha$-line luminosity of PSO J006.1240+39.2219 is unusually high and estimated to be 0.8$\times$10$^{12}$ Solar luminosities (about 3% of the total quasar luminosity). The Lyman-$\alpha$ emission of PSO J006.1240+39.2219 shows fast variability on timescales of days in the quasar rest frame, which has never been detected in any of the known high-redshift quasars. The high luminosity of the Lyman-$\alpha$ line, its narrow width and fast variability resemble properties of local Narrow-Line Seyfert 1 galaxies which suggests that the quasar is likely at the active phase of the black hole growth accreting close or even beyond the Eddington limit.

We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.

After several years of quiescence, the blazar CTA 102 underwent an exceptional outburst in 2012 September-October. The flare was tracked from gamma-ray to near-infrared frequencies, including Fermi and Swift data as well as photometric and polarimetric data from several observatories. An intensive GASP-WEBT collaboration campaign in optical and NIR bands, with an addition of previously unpublished archival data and extension through fall 2015, allows comparison of this outburst with the previous activity period of this blazar in 2004-2005. We find remarkable similarity between the optical and gamma-ray behaviour of CTA 102 during the outburst, with a time lag between the two light curves of ~1 hour, indicative of co-spatiality of the optical and gamma-ray emission regions. The relation between the gamma-ray and optical fluxes is consistent with the SSC mechanism, with a quadratic dependence of the SSC gamma-ray flux on the synchrotron optical flux evident in the post-outburst stage. However, the gamma-ray/optical relationship is linear during the outburst; we attribute this to changes in the Doppler factor. A strong harder-when-brighter spectral dependence is seen both the in gamma-ray and optical non-thermal emission. This hardening can be explained by convexity of the UV-NIR spectrum that moves to higher frequencies owing to an increased Doppler shift as the viewing angle decreases during the outburst stage. The overall pattern of Stokes parameter variations agrees with a model of a radiating blob or shock wave that moves along a helical path down the jet.

We study the multi-band variability and correlations of the TeV blazar Mrk 421 on year time scales, which can bring additional insight on the processes responsible for its broadband emission. We observed Mrk 421 in the very high energy (VHE) gamma-ray range with the Cherenkov telescope MAGIC-I from March 2007 to June 2009 for a total of 96 hours of effective time after quality cuts. The VHE flux variability is quantified with several methods, including the Bayesian Block algorithm, which is applied to data from Cherenkov telescopes for the first time. The 2.3 year long MAGIC light curve is complemented with data from the Swift/BAT and RXTE/ASM satellites and the KVA, GASP-WEBT, OVRO, and Mets\"ahovi telescopes from February 2007 to July 2009, allowing for an excellent characterisation of the multi-band variability and correlations over year time scales. Mrk 421 was found in different gamma-ray emission states during the 2.3 year long observation period. Flares and different levels of variability in the gamma-ray light curve could be identified with the Bayesian Block algorithm. The same behaviour of a quiet and active emission was found in the X-ray light curves measured by Swift/BAT and the RXTE/ASM, with a direct correlation in time. The behaviour of the optical light curve of GASP-WEBT and the radio light curves by OVRO and Mets\"ahovi are different as they show no coincident features with the higher energetic light curves and a less variable emission. The fractional variability is overall increasing with energy. The comparable variability in the X-ray and VHE bands and their direct correlation during both high- and low-activity periods spanning many months show that the electron populations radiating the X-ray and gamma-ray photons are either the same, as expected in the Synchrotron-Self-Compton mechanism, or at least strongly correlated, as expected in electromagnetic cascades.

We present infrared galaxy luminosity functions (LFs) in the AKARI North Ecliptic Pole (NEP) deep field using recently-obtained, wider CFHT optical/near-IR images. AKARI has obtained deep images in the mid-infrared (IR), covering 0.6 deg$^2$ of the NEP deep field. However, our previous work was limited to the central area of 0.25 deg$^2$ due to the lack of optical coverage of the full AKARI NEP survey. To rectify the situation, we recently obtained CFHT optical and near-IR images over the entire AKARI NEP deep field. These new CFHT images are used to derive accurate photometric redshifts, allowing us to fully exploit the whole AKARI NEP deep field. AKARI's deep, continuous filter coverage in the mid-IR wavelengths (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24$\mu$m) exists nowhere else, due to filter gaps of other space telescopes. It allows us to estimate restframe 8$\mu$m and 12$\mu$m luminosities without using a large extrapolation based on spectral energy distribution (SED) fitting, which was the largest uncertainty in previous studies. Total infrared luminosity (TIR) is also obtained more reliably due to the superior filter coverage. The resulting restframe 8$\mu$m, 12$\mu$m, and TIR LFs at $0.15<z<2.2$ are consistent with previous works, but with reduced uncertainties, especially at the high luminosity-end, thanks to the wide field coverage. In terms of cosmic infrared luminosity density ($\Omega_{\mathrm{IR}}$), we found that the $\Omega_{\mathrm{IR}}$ evolves as $\propto (1+z)^{4.2\pm 0.4}$.

A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameters, and evaluate whether the observed broadband SED variability can be associated to variations in the relativistic particle population. Flux variability was remarkable in the X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone SSC model can describe reasonably well the SED of each day for the 13 consecutive days. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission while the other smaller zone, which is spatially separated from the first one, contributes to the daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly to the underlying particle population. This shows that the particle acceleration and cooling mechanism producing the radiating particles could be the main one responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by the variation of the parameters related to the emitting region itself ($\delta$, $B$ and $R$), in addition to the parameters related to the particle population.

We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.

Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.

Aims: We performed a detailed photometric analysis of the lensed system UM673 (Q0142-100) and an analysis of the tentative lens models. Methods: High-resolution adaptive optics images of UM673 taken with the Subaru telescope in the H band were examined. We also analysed the J, H and K-band observational data of UM673 obtained with the 1.3m telescope at the CTIO observatory. Results: We present photometry of quasar components A and B of UM673, the lens, and the nearby bright galaxy using H-band observational data obtained with the Subaru telescope. Based on the CTIO observations of UM673, we also present J- and H-band photometry and estimates of the J, H and K-band flux ratios between the two UM673 components in recent epochs. The near-infrared fluxes of the A and B components of UM673 and their published optical fluxes are analysed to measure extinction properties of the lensing galaxy. We estimate the extinction-corrected flux ratio between components A and B to be about 2.14 mag. We discuss lens models for the UM673 system constrained with the positions of the UM673 components, their flux ratio, and the previously measured time delay

The quasar-type blazar 4C 38.41 (B3 1633+382) experienced a large outburst in 2011, which was detected throughout the entire electromagnetic spectrum. We present the results of low-energy multifrequency monitoring by the GASP project of the WEBT consortium and collaborators, as well as those of spectropolarimetric/spectrophotometric monitoring at the Steward Observatory. We also analyse high-energy observations of the Swift and Fermi satellites. In the optical-UV band, several results indicate that there is a contribution from a QSO-like emission component, in addition to both variable and polarised jet emission. The unpolarised emission component is likely thermal radiation from the accretion disc that dilutes the jet polarisation. We estimate its brightness to be R(QSO) ~ 17.85 - 18 and derive the intrinsic jet polarisation degree. We find no clear correlation between the optical and radio light curves, while the correlation between the optical and \gamma-ray flux apparently fades in time, likely because of an increasing optical to \gamma-ray flux ratio. As suggested for other blazars, the long-term variability of 4C 38.41 can be interpreted in terms of an inhomogeneous bent jet, where different emitting regions can change their alignment with respect to the line of sight, leading to variations in the Doppler factor \delta. Under the hypothesis that in the period 2008-2011 all the \gamma-ray and optical variability on a one-week timescale were due to changes in \delta, this would range between ~ 7 and ~ 21. If the variability were caused by changes in the viewing angle \theta\ only, then \theta\ would go from ~ 2.6 degr to ~ 5 degr. Variations in the viewing angle would also account for the dependence of the polarisation degree on the source brightness in the framework of a shock-in-jet model.

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi-wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope mission. While investigating the previously reported gamma-ray/optical flare accompanied by a change in optical polarization, we found that the optical emission appears delayed with respect to the gamma-ray emission by about 10 days. X-ray observations reveal a pair of `isolated' flares separated by ~90 days, with only weak gamma-ray/optical counterparts. The spectral structure measured by Spitzer reveals a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared band during the gamma-ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected spectral energy distributions are fitted with leptonic models including Comptonization of external radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the polarization swing involving propagation of the emitting region along a curved trajectory, we can explain the evolution of the broad-band spectra during the gamma-ray flaring event by a shift of its location from ~ 1 pc to ~ 4 pc from the central black hole. On the other hand, if the gamma-ray flare is generated instead at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking synchrotron component, and suggest that the corresponding inverse-Compton component explains the steady X-ray emission.

We study brightness variations in the double lensed quasar UM673 (Q0142-100) with the aim of measuring the time delay between its two images. In the paper we combine our previously published observational data of UM673 obtained during the 2003 - 2005 seasons at the Maidanak Observatory with archival and recently observed Maidanak and CTIO UM673 data. We analyze the V, R and I-band light curves of the A and B images of UM673, which cover ten observational seasons from August 2001 to November 2010. We also analyze the time evolution of the difference in magnitudes between images A and B of UM673 over more than ten years. We find that the quasar exhibits both short-term (with amplitude of \sim 0.1 mag in the R band) and high-amplitude (\sim 0.3 mag) long-term variability on timescales of about several months and several years, respectively. These brightness variations are used to constrain the time delay between the images of UM673. From cross-correlation analysis of the A and B quasar light curves and error analysis we measure the mean time delay and its error of 89 \pm11 days. Given the input time delay of 88 days, the most probable value of the delay that can be recovered from light curves with the same statistical properties as the observed R-band light curves of UM673 is 95{+5/-16}{+14/-29} days (68 and 95 % confidence intervals). Analysis of the V - I color variations and V, R and I-band magnitude differences of the quasar images does not show clear evidence of the microlensing variations between 1998 and 2010.

For more than 15 years, since the days of the Energetic Gamma-Ray Experiment Telescope (EGRET) on board the Compton Gamma-Ray Observatory (CGRO; 1991-2000), it has remained an open question why the prominent blazar 3C 345 was not reliably detected at gamma-ray energies <=20 MeV. Recently a bright gamma-ray source (0FGL J1641.4+3939/1FGL J1642.5+3947), potentially associated with 3C 345, was detected by the Large Area Telescope (LAT) on Fermi. Multiwavelength observations from radio bands to X-rays (mainly GASP-WEBT and Swift) of possible counterparts (3C 345, NRAO 512, B3 1640+396) were combined with 20 months of Fermi-LAT monitoring data (August 2008 - April 2010) to associate and identify the dominating gamma-ray emitting counterpart of 1FGL J1642.5+3947. The source 3C 345 is identified as the main contributor for this gamma-ray emitting region. However, after November 2009 (15 months), a significant excess of photons from the nearby quasar NRAO 512 started to contribute and thereafter was detected with increasing gamma-ray activity, possibly adding flux to 1FGL J1642.5+3947. For the same time period and during the summer of 2010, an increase of radio, optical and X-ray activity of NRAO 512 was observed. No gamma-ray emission from B3 1640+396 was detected.

We report on the extreme gamma-ray activity from the FSRQ PKS 1510-089 observed by AGILE in March 2009. In the same period a radio-to-optical monitoring of the source was provided by the GASP-WEBT and REM. Moreover, several Swift ToO observations were triggered, adding important information on the source behaviour from optical/UV to hard X-rays. We paid particular attention to the calibration of the Swift/UVOT data to make it suitable to the blazars spectra. Simultaneous observations from radio to gamma rays allowed us to study in detail the correlation among the emission variability at different frequencies and to investigate the mechanisms at work. In the period 9-30 March 2009, AGILE detected an average gamma-ray flux of (311+/-21)x10^-8 ph cm^-2 s^-1 for E>100 MeV, and a peak level of (702+/-131)x10^-8 ph cm^-2 s^-1 on daily integration. The gamma-ray activity occurred during a period of increasing activity from near-IR to UV, with a flaring episode detected on 26-27 March 2009, suggesting that a single mechanism is responsible for the flux enhancement observed from near-IR to UV. By contrast, Swift/XRT observations seem to show no clear correlation of the X-ray fluxes with the optical and gamma-ray ones. However, the X-ray observations show a harder photon index (1.3-1.6) with respect to most FSRQs and a hint of harder-when-brighter behaviour, indicating the possible presence of a second emission component at soft X-ray energies. Moreover, the broad band spectrum from radio-to-UV confirmed the evidence of thermal features in the optical/UV spectrum of PKS 1510-089 also during high gamma-ray state. On the other hand, during 25-26 March 2009 a flat spectrum in the optical/UV energy band was observed, suggesting an important contribution of the synchrotron emission in this part of the spectrum during the brightest gamma-ray flare, therefore a significant shift of the synchrotron peak.

The GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) monitored BL Lacertae in 2008-2009 at radio, near-IR, and optical frequencies. During this period, high-energy observations were performed by XMM-Newton, Swift, and Fermi. We analyse these data with particular attention to the calibration of Swift UV data, and apply a helical jet model to interpret the source broad-band variability. The GASP-WEBT observations show an optical flare in 2008 February-March, and oscillations of several tenths of mag on a few-day time scale afterwards. The radio flux is only mildly variable. The UV data from both XMM-Newton and Swift seem to confirm a UV excess that is likely caused by thermal emission from the accretion disc. The X-ray data from XMM-Newton indicate a strongly concave spectrum, as well as moderate flux variability on an hour time scale. The Swift X-ray data reveal fast (interday) flux changes, not correlated with those observed at lower energies. We compare the spectral energy distribution (SED) corresponding to the 2008 low-brightness state, which was characterised by a synchrotron dominance, to the 1997 outburst state, where the inverse-Compton emission was prevailing. A fit with an inhomogeneous helical jet model suggests that two synchrotron components are at work with their self inverse-Compton emission. Most likely, they represent the radiation from two distinct emitting regions in the jet. We show that the difference between the source SEDs in 2008 and 1997 can be explained in terms of pure geometrical variations. The outburst state occurred when the jet-emitting regions were better aligned with the line of sight, producing an increase of the Doppler beaming factor. Our analysis demonstrates that the jet geometry can play an extremely important role in the BL Lacertae flux and spectral variability.

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.

We report the detection by the AGILE satellite of a rapid gamma-ray flare from the powerful gamma-ray quasar PKS 1510-089, during a pointing centered on the Galactic Center region from 1 March to 30 March 2008. This source has been continuosly monitored in the radio-to-optical bands by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT). Moreover, the gamma-ray flaring episode triggered three ToO observations by the Swift satellite in three consecutive days, starting from 20 March 2008. In the period 1-16 March 2008, AGILE detected gamma-ray emission from PKS 1510-089 at a significance level of 6.2-sigma with an average flux over the entire period of (84 +/- 17) x 10^{-8} photons cm^{-2} s^{-1} for photon energies above 100 MeV. After a predefined satellite re-pointing, between 17 and 21 March 2008, AGILE detected the source at a significance level of 7.3-sigma, with an average flux (E > 100 MeV) of (134 +/- 29) x 10^{-8} photons cm^{-2} s^{-1} and a peak level of (281 +/- 68) x 10^{-8} photons cm^{-2} s^{-1} with daily integration. During the observing period January-April 2008, the source also showed an intense and variable optical activity, with several flaring episodes and a significant increase of the flux was observed at millimetric frequencies. Moreover, in the X-ray band the Swift/XRT observations seem to show an harder-when-brighter behaviour of the source spectrum. The spectral energy distribution of mid-March 2008 is modelled with a homogeneous one-zone synchrotron self Compton emission plus contributions from inverse Compton scattering of external photons from both the accretion disc and the broad line region. Indeed, some features in the optical-UV spectrum seem to indicate the presence of Seyfert-like components, such as the little blue bump and the big blue bump.

We present and analyse new R-band frames of the gravitationally lensed double quasar FBQ 0951+2635. These images were obtained with the 1.5m AZT-22 Telescope at Maidanak (Uzbekistan) in the 2001-2006 period. Previous results in the R band (1999-2001 period) and the new data allow us to discuss the dominant kind of microlensing variability in FBQ 0951+2635. The time evolution of the flux ratio A/B does not favour the continuous production of short-timescale (months) flares in the faintest quasar component B (crossing the central region of the lensing galaxy). Instead of a rapid variability scenario, the observations are consistent with the existence of a long-timescale fluctuation. The flux ratio shows a bump in the 2003-2004 period and a quasi-flat trend in more recent epochs. Apart from the global behaviour of A/B, we study the intra-year variability over the first semester of 2004, which is reasonably well sampled. Short-timescale microlensing is not detected in that period. Additional data in the i band (from new i-band images taken in 2007 with the 2m Liverpool Robotic Telescope at La Palma, Canary Islands) also indicate the absence of short-timescale events in 2007.

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.

To go into the details about the variability of the double quasar SBS 0909+532, we designed a monitoring programme with the 2 m Liverpool Robotic Telescope in the r Sloan filter, spanning 1.5 years from 2005 January to 2006 June. The r-band light curves of the A and B components, several cross-correlation techniques and a large number of simulations (synthetic light curves) lead to a robust delay of 49 +/- 6 days (1-sigma interval) that agrees with our previous results (the B component is leading). Once the time delay and the magnitude offset are known, the magnitude- and time-shifted light curve of image A is subtracted from the light curve of image B. This difference light curve of SBS 0909+532 is consistent with zero, so any possible extrinsic signal must be very weak, i.e., the observed variability in A and B is basically due to observational noise and intrinsic signal. We then make the combined light curve and analyse its statistical properties (structure functions). The structure function of the intrinsic luminosity is fitted to predictions of simple models of two physical scenarios: accretion disc instabilities and nuclear starbursts. Although no simple model is able to accurately reproduce the observed trend, symmetric triangular flares in an accretion disc seems to be the best option to account for it.

The aim of this contribution is to present the two first phases of the optical monitoring programme of the Gravitational Lenses group at the Universidad de Cantabria (GLUC, http://grupos.unican.es/glendama/). In an initial stage (2003 March-June), the Estacion de Observacion de Calar Alto (EOCA) was used to obtain VR frames of SBS 0909+532 and QSO 0957+561. These observations in 2003 led to accurate fluxes of the two components of both double QSOs, which are being compared and complemented with data from other 1-1.5 m telescopes located in the North Hemisphere: Fred Lawrence Whipple Observatory (USA), Maidanak Observatory (Uzbekistan) and Wise Observatory (Israel). On the other hand, the GLUC started the second phase of its monitoring programme in 2005 January. In this second phase, they are using the 2 m fully robotic Liverpool Telescope (LT). The key idea is the two-band photometric follow-up of four lensed QSOs with different main lensing galaxies: SBS 0909+532 (elliptical), QSO 0957+561 (giant cD), B1600+434 (edge-on spiral) and QSO 2237+0305 (face-on spiral). Thus, the light rays associated with the components of the four gravitational mirages cross different galaxy environments, and the corresponding light curves could unveil the content of these environments. While SBS 0909+532 and QSO 0957+561 are the targets for the two first years with the LT (2005-2006), the rest of targets (B1600+434 and QSO 2237+0305) will be monitored starting from 2007.

The time delays between the components of a lensed quasar are basic tools to analyze the expansion of the Universe and the structure of the main lens galaxy halo. In this paper, we focus on the variability and time delay of the double system SBS 0909+532A,B as well as the time behaviour of the field stars. We use VR optical observations of SBS 0909+532A,B and the field stars in 2003. The frames were taken at Calar Alto, Maidanak and Wise observatories, and the VR light curves of the field stars and quasar components are derived from aperture and point-spread function fitting methods. We measure the R-band time delay of the system from the chi-square and dispersion techniques and 1000 synthetic light curves based on the observed records. One nearby field star (SBS 0909+532c) is found to be variable, and the other two nearby field stars are non-variable sources. With respect to the quasar components, the R-band records seem more reliable and are more densely populated than the V-band ones. The observed R-band fluctuations permit a pre-conditioned measurement of the time delay. From the chi-square minimization, if we assume that the quasar emission is observed first in B and afterwards in A (in agreement with basic observations of the system and the corresponding predictions), we obtain a delay of - 45 (+ 1)/(- 11) days (95% confidence interval). The dispersion technique leads to a similar delay range. A by-product of the analysis is the determination of a totally corrected flux ratio in the R band (corrected by the time delay and the contamination due to the galaxy light). Our 95% measurement of this ratio (0.575 +/- 0.014 mag) is in excellent agreement with previous results from contaminated fluxes at the same time of observation.

The analysis of the high magnification events in the A and C components of the quadruple gravitational lens QSO2237+0305 observed by OGLE and GLITP collaborations in V band was carried out. The significant light amplifications of the components are interpreted as the effect of microlensing with a fold caustic. For the reconstruction of the one-dimensional source profile the technique based on Tikhonov regularization method was used. The estimates of the effective radius of the quasar's emitting region (the radius within which half of the light is emitted) based on reconstructed profile of the source from microlensing of the A and C components are in the range of 31 and 21 days and correspond to the linear sizes 0.62e+15 cm and 0.42e+15 cm. For the A component the positive crossing of the caustic and for the C component the negative crossing of the caustic was confirmed.

We have observed the gravitational lens system Q2237+0305 from the Maidanak Observatory over the period from August 2002 to November 2003. Here we report the results of our observations. We implemented a two-stage technique that has been developed specifically for the purpose of gravitational lens image reconstruction. The technique is based on the Tikhonov regularization approach and allows one to obtain astrometric and photometric characteristics of the gravitational lens system. Light curves with 78 data points for the four quasar components are obtained. Slow brightness variations over the observational period are found in all components. Images A, C, D have a tendency to decrease in brightness. Image B does not vary more than 0.05mag. The observations did not reveal evidence for large variations in brightness of the components due to microlensing effects. To provide an overall picture of the photometry behaviour, our data are combined with the Maidanak observations published for 1995 -- 2000.

This note addresses possible applications of the Tikhonov regularization to image reconstruction of gravitational lens systems. Several modifications of the regularization algorithm are discussed. Our illustrative example is the close quadruple gravitational lens QSO 2237+0305 (Einstein Cross). The restored image of the lens is decomposed into two parts - the quasar components and the background galaxy.