We provide a detailed description of the Herschel-SPIRE Fourier Transform Spectrometer (FTS) Spectral Feature Finder (FF). The FF is an automated process designed to extract significant spectral features from SPIRE FTS data products. Optimising the number of features found in SPIRE-FTS spectra is challenging. The wide SPIRE-FTS frequency range (447-1568 GHz) leads to many molecular species and atomic fine structure lines falling within the observed bands. As the best spectral resolution of the SPIRE-FTS is ~1.2 GHz, there can be significant line blending, depending on the source type. In order to find, both efficiently and reliably, features in spectra associated with a wide range of sources, the FF iteratively searches for peaks over a number of signal-to-noise ratio (SNR) thresholds. For each threshold, newly identified features are rigorously checked before being added to the fitting model. At the end of each iteration, the FF simultaneously fits the continuum and features found, with the resulting residual spectrum used in the next iteration. The final FF products report the frequency of the features found and the associated SNRs. Line flux determination is not included as part of the FF products, as extracting reliable line flux from SPIRE-FTS data is a complex process that requires careful evaluation and analysis of the spectra on a case-by-case basis. The FF results are 100% complete for features with SNR greater than 10 and 50-70% complete at SNR of 5. The FF code and all FF products are publicly available via the Herschel Science Archive.

The advent of the Fermi-GST with its unprecedented capability to monitor the entire 4 pi sky within less than 2-3 hours, introduced new standard in time domain gamma-ray astronomy. To explore this new avenue of extragalactic physics the F-GAMMA programme undertook the task of conducting nearly monthly, broadband radio monitoring of selected blazars from January 2007 to January 2015. In this work we release all the light curves at 2.64, 4.85, 8.35, 10.45, 14.6, 23.05, 32, and 43 GHz and present first order derivative data products after all necessary post-measurement corrections and quality checks; that is flux density moments and spectral indices. The release includes 155 sources. The effective cadence after the quality flagging is around one radio SED every 1.3 months. The coherence of each radio SED is around 40 minutes. The released dataset includes more than $4\times10^4$ measurements. The median fractional error at the lowest frequencies (2.64-10.45 GHz) is below 2%. At the highest frequencies (14.6-43 GHz) with limiting factor of the atmospheric conditions, the errors range from 3% to 9%, respectively.

Context. In recent years, an in-depth gamma-ray analysis of the Orion region has been carried out by the AGILE and Fermi-LAT (Large Area Telescope) teams with the aim of estimating the H2-CO conversion factor, XCO. The comparison of the data from both satellites with models of diffuse gamma-ray Galactic emission unveiled an excess at (l,b)=[213.9, -19.5], in a region at a short angular distance from the OB star k-Ori. Possible explanations of this excess are scattering of the so-called "dark gas", non-linearity in the H2-CO relation, or Cosmic-Ray (CR) energization at the k-Ori wind shock. Aims. Concerning this last hypothesis, we want to verify whether cosmic-ray acceleration or re-acceleration could be triggered at the k-Ori forward shock, which we suppose to be interacting with a star-forming shell detected in several wavebands and probably triggered by high energy particles. Methods. Starting from the AGILE spectrum of the detected gamma-ray excess, showed here for the first time, we developed a valid physical model for cosmic-ray energization, taking into account re-acceleration, acceleration, energy losses, and secondary electron contribution. Results. Despite the characteristic low velocity of an OB star forward shock during its "snowplow" expansion phase, we find that the Orion gamma-ray excess could be explained by re-acceleration of pre-existing cosmic rays in the interaction between the forward shock of k-Ori and the CO-detected, star-forming shell swept-up by the star expansion. According to our calculations, a possible contribution from freshly accelerated particles is sub-dominant with respect the re-acceleration contribution. However, a simple adiabatic compression of the shell could also explain the detected gamma-ray emission. Futher GeV and TeV observations of this region are highly recommended in order to correctly identify the real physical scenario.

Narrow line Seyfert 1 galaxies (NLSy1s) constitute the AGN subclass associated with systematically smaller black hole masses. A few radio loud ones have been detected in MeV -- GeV energy bands by Fermi and evidence for the presence of blazar-like jets has been accumulated. In this study we wish to quantify the temporal behaviour of the optical polarisation, fraction and angle, for a selected sample of radio loud NLSy1s. We also search for rotations of the polarisation plane similar to those commonly observed in blazars. We have conducted R-band optical polarisation monitoring of a sample of 10 RL NLSy1s 5 of which have been previously detected by Fermi. The dataset includes observations with the RoboPol, KANATA, Perkins and Steward polarimeters. In the cases where evidences for long rotations of the polarisation plane are found, we carry out numerical simulations to assess the probability that they are caused by intrinsically evolving EVPAs instead of observational noise. Even our moderately sampled sources show indications of variability, both in polarisation fraction and angle. For the four best sampled objects in our sample we find multiple periods of significant polarisation angle variability. In the two best sampled cases, namely J1505+0326 and J0324+3410, we find indications for three long rotations. We show that although noise can induce the observed behaviour, it is much more likely that the apparent rotation is caused by intrinsic evolution of the EVPA. To our knowledge this is the very first detection of such events in this class of sources. In the case of the largest dataset (J0324+3410) we find that the EVPA concentrates around a direction which is at 49.3\degr to the 15-GHz radio jet implying a projected magnetic field at an angle of 40.7\degr to that axis.

We present a catalogue of ALMA flux density measurements of 754 calibrators observed between August 2012 and September 2017, for a total of 16,263 observations in different bands and epochs. The flux densities were measured reprocessing the ALMA images generated in the framework of the ALMACAL project, with a new code developed by the Italian node of the European ALMA Regional Centre. A search in the online databases yielded redshift measurements for 589 sources ($\sim$78 per cent of the total). Almost all sources are flat-spectrum, based on their low-frequency spectral index, and have properties consistent with being blazars of different types. To illustrate the properties of the sample we show the redshift and flux density distributions as well as the distributions of the number of observations of individual sources and of time spans in the source frame for sources observed in bands 3 (84$-$116 GHz) and 6 (211$-$275 GHz). As examples of the scientific investigations allowed by the catalogue we briefly discuss the variability properties of our sources in ALMA bands 3 and 6 and the frequency spectra between the effective frequencies of these bands. We find that the median variability index steadily increases with the source-frame time lag increasing from 100 to 800 days, and that the frequency spectra of BL Lacs are significantly flatter than those of flat-spectrum radio quasars. We also show the global spectral energy distributions of our sources over 17 orders of magnitude in frequency.

Diffuse galactic gamma-ray emission is produced by the interaction of cosmic rays (CRs) with the interstellar environment. The study of gamma-ray emission is therefore a powerful tool to investigate the origin of CRs and the processes through which they are accelerated. We aim to gain deeper insights of the nature of gamma-ray emission in the region of Orion, which is one of the best studied sites of on-going star formation, by analysing data from the AGILE satellite. The diffuse gamma-ray emission expected from the Orion region is relatively high. Its separation from the galactic plane also ensures a very low contribution from foreground or background emission, which makes it an ideal site for studying the processes of particle acceleration in star forming environments. The AGILE data are modelled through a template that quantifies the gamma-ray diffuse emission expected from atomic and molecular hydrogen. Other sources of emission are modelled as an isotropic contribution. Gamma-ray emission exceeding the amount expected by the diffuse emission model is detected with high level of significance. The main excess is in the high-longitude part of Orion A. A thorough analysis of this feature suggests a connection between the observed gamma-ray emission and the B0.5 Ia star k Orionis. The location of the gamma-ray excess is compatible with the site where stellar wind collides with the ISM. Both scattering on dark gas and cosmic-ray acceleration at the shock between the two environments are discussed as possible explanations, with the latter hypothesis being supported by the hardness of the energy spectrum of the emission. If confirmed, this would be the first direct detection of gamma-ray emission from the interaction between ISM and a single star's stellar wind.

We describe an update to the Herschel-SPIRE Fourier-Transform Spectrometer (FTS) calibration for extended sources, which incorporates a correction for the frequency-dependent far-field feedhorn efficiency, $\eta_\mathrm{FF}$. This significant correction affects all FTS extended-source calibrated spectra in sparse or mapping mode, regardless of the spectral resolution. Line fluxes and continuum levels are underestimated by factors of 1.3-2 in the Spectrometer Long-Wavelength band (SLW, 447-1018 GHz; 671-294 $\mu$m) and 1.4-1.5 in the Spectrometer Short-Wavelength band (SSW, 944-1568 GHz; 318-191 $\mu$m). The correction was implemented in the FTS pipeline version 14.1 and has also been described in the SPIRE Handbook since Feb 2017. Studies based on extended-source calibrated spectra produced prior to this pipeline version should be critically reconsidered using the current products available in the Herschel Science Archive. Once the extended-source calibrated spectra are corrected for $\eta_\mathrm{FF}$, the synthetic photometry and the broadband intensities from SPIRE photometer maps agree within 2-4% -- similar levels to the comparison of point-source calibrated spectra and photometry from point-source calibrated maps. The two calibration schemes for the FTS are now self-consistent: the conversion between the corrected extended-source and point-source calibrated spectra can be achieved with the beam solid angle and a gain correction that accounts for the diffraction loss.

Black holes, anywhere in the stellar-mass to supermassive range, are often associated with relativistic jets. Models suggest that jet production may be a universal process common in all black hole systems regardless of their mass. Although in many cases observations support such hypotheses for microquasars and Seyfert galaxies, little is known on whether boosted blazar jets also comply with such universal scaling laws. We use uniquely rich multiwavelength radio light curves from the F-GAMMA program and the most accurate Doppler factors available to date to probe blazar jets in their emission rest frame with unprecedented accuracy. We identify for the first time a strong correlation between the blazar intrinsic broad-band radio luminosity and black hole mass, which extends over $\sim$ 9 orders of magnitude down to microquasars scales. Our results reveal the presence of a universal scaling law that bridges the observing and emission rest frames in beamed sources and allows us to effectively constrain jet models. They consequently provide an independent method for estimating the Doppler factor, and for predicting expected radio luminosities of boosted jets operating in systems of intermediate or tens-of-solar mass black holes, immediately applicable to cases as those recently observed by LIGO.

We present a pipeline that allows recovering reliable information for all four Stokes parameters with high accuracy. Its novelty relies on the treatment of the instrumental effects already prior to the computation of the Stokes parameters contrary to conventional methods, such as the M\"uller matrix one. The instrumental linear polarization is corrected across the whole telescope beam and significant Stokes $Q$ and $U$ can be recovered even when the recorded signals are severely corrupted. The accuracy we reach in terms of polarization degree is of the order of 0.1-0.2 %. The polarization angles are determined with an accuracy of almost 1$^{\circ}$. The presented methodology was applied to recover the linear and circular polarization of around 150 Active Galactic Nuclei. The sources were monitored from July 2010 to April 2016 with the Effelsberg 100-m telescope at 4.85 GHz and 8.35 GHz with a cadence of around 1.2 months. The polarized emission of the Moon was used to calibrate the polarization angle. Our analysis showed a small system-induced rotation of about 1$^{\circ}$ at both observing frequencies. Finally, we identify five sources with significant and stable linear polarization; three sources remain constantly linearly unpolarized over the period we examined; a total of 11 sources have stable circular polarization degree $m_\mathrm{c}$ and four of them with non-zero $m_\mathrm{c}$. We also identify eight sources that maintain a stable polarization angle over the examined period. All this is provided to the community for polarization observations reference. We finally show that our analysis method is conceptually different from the traditionally used ones and performs better than the M\"uller matrix method. Although it was developed for a system equipped with circularly polarized feeds it can easily be modified for systems with linearly polarized feeds as well.

Hi-GAL is a large-scale survey of the Galactic plane, performed with Herschel in five infrared continuum bands between 70 and 500 $\mu$m. We present a band-merged catalogue of spatially matched sources and their properties derived from fits to the spectral energy distributions (SEDs) and heliocentric distances, based on the photometric catalogs presented in Molinari et al. (2016a), covering the portion of Galactic plane $-71.0^{\circ}< \ell < 67.0^{\circ}$. The band-merged catalogue contains 100922 sources with a regular SED, 24584 of which show a 70 $\mu$m counterpart and are thus considered proto-stellar, while the remainder are considered starless. Thanks to this huge number of sources, we are able to carry out a preliminary analysis of early stages of star formation, identifying the conditions that characterise different evolutionary phases on a statistically significant basis. We calculate surface densities to investigate the gravitational stability of clumps and their potential to form massive stars. We also explore evolutionary status metrics such as the dust temperature, luminosity and bolometric temperature, finding that these are higher in proto-stellar sources compared to pre-stellar ones. The surface density of sources follows an increasing trend as they evolve from pre-stellar to proto-stellar, but then it is found to decrease again in the majority of the most evolved clumps. Finally, we study the physical parameters of sources with respect to Galactic longitude and the association with spiral arms, finding only minor or no differences between the average evolutionary status of sources in the fourth and first Galactic quadrants, or between "on-arm" and "inter-arm" positions.

Recent population studies have shown that the variability Doppler factors can adequately describe blazars as a population. We use the flux density variations found within the extensive radio multi-wavelength datasets of the F-GAMMA program, a total of 10 frequencies from 2.64 up to 142.33 GHz, in order to estimate the variability Doppler factors for 58 $\gamma$-ray bright sources, for 20 of which no variability Doppler factor has been estimated before. We employ specifically designed algorithms in order to obtain a model for each flare at each frequency. We then identify each event and track its evolution through all the available frequencies for each source. This approach allows us to distinguish significant events producing flares from stochastic variability in blazar jets. It also allows us to effectively constrain the variability brightness temperature and hence the variability Doppler factor as well as provide error estimates. Our method can produce the most accurate (16\% error on average) estimates in the literature to date.

The SPIRE Fourier Transform Spectrometer on board the Herschel Space Observatory had two standard spectral resolution modes for science observations: high resolution (HR) and low resolution (LR), which could also be performed in sequence (H+LR). A comparison of the HR and LR resolution spectra taken in this sequential mode, revealed a systematic discrepancy in the continuum level. Analysing the data at different stages during standard pipeline processing, demonstrates the telescope and instrument emission affect HR and H+LR observations in a systematically different way. The origin of this difference is found to lie in the variation of both the telescope and instrument response functions, while it is triggered by fast variation of the instrument temperatures. As it is not possible to trace the evolution of the response functions through auxiliary housekeeping parameters, the calibration cannot be corrected analytically. Therefore an empirical correction for LR spectra has been developed, which removes the systematic noise introduced by the variation of the response functions.

To fully exploit the scientific potential of the Fermi mission, we initiated the F-GAMMA program. Between 2007 and 2015 it was the prime provider of complementary multi-frequency monitoring in the radio regime. We quantify the radio variability of gamma-ray blazars. We investigate its dependence on source class and examine whether the radio variability is related to the gamma-ray loudness. Finally, we assess the validity of a putative correlation between the two bands. The F-GAMMA monitored monthly a sample of about 60 sources at up to twelve radio frequencies between 2.64 and 228.39 GHz. We perform a time series analysis on the first 2.5-year dataset to obtain variability parameters. A maximum likelihood analysis is used to assess the significance of a correlation between radio and gamma-ray fluxes. We present light curves and spectra (coherent within ten days) obtained with the Effelsberg 100-m and IRAM 30-m telescopes. All sources are variable across all frequency bands with amplitudes increasing with frequency up to rest frame frequencies of around 60 - 80 GHz as expected by shock-in-jet models. Compared to FSRQs, BL Lacs show systematically lower variability amplitudes, brightness temperatures and Doppler factors at lower frequencies, while the difference vanishes towards higher ones. The time scales appear similar for the two classes. The distribution of spectral indices appears flatter or more inverted at higher frequencies for BL Lacs. Evolving synchrotron self-absorbed components can naturally account for the observed spectral variability. We find that the Fermi-detected sources show larger variability amplitudes as well as brightness temperatures and Doppler factors, than non-detected ones. Flux densities at 86.2 and 142.3 GHz correlate with 1 GeV fluxes at a significance level better than 3sigma, implying that gamma rays are produced very close to the mm-band emission region.

The IDV observations of S4 0917+624 were carried out monthly, from August 2005 to January 2010, with the Urumqi 25m radio telescope at 4.8 GHz. The quasar S4 0917+624 exhibits only very weak or no IDV during our 4.5 year observing interval. Prior to the year 2000, the source S4 0917+624 was one of the most prominent IDV sources. Our new data indicate that the previous strong IDV has ceased. We analyzed the long-term VLBI structural variability using Gaussian model-fitting. From this we obtained the flux densities and the deconvolved sizes of core and inner-jet components of the source. We studied the properties such as core fraction, angular size, spectral index, and brightness temperature of VLBI core for S4 0917+624, as well as the time delay between 5 and 15 GHz variations, and compared them with the IDV properties of S4 0917+624. The source shows ejection of several jet components that are suspected to have partially reduced the IDV amplitude of S4 0917+624. However, during 2005-2006, the VLBI core size was comparable to the size before the year 2000, but no strong IDV was detected in the period, suggesting that the quenching effect due to source size changes may not be responsible for the lack of strong IDV after the year 2000. The refractive scattering properties for the strong IDV phase of S4 0917+624 before the year 2000 are discussed. The disappearance of strong IDV in S4 0917+624 after the year 2000 is a mystery and cannot be explained via the quenching effect by changes in the observable VLBI structure. However, it may be caused by changes in the interstellar medium, i.e. by interstellar weather, which induces changes in the scintillation pattern on timescales of several years. Further coordinated multi-frequency observations will be required to distinguish between the effect of source-intrinsic variability and changing properties of the interstellar medium.

A systematic programme of calibration observations was carried out to monitor the performance of the SPIRE FTS instrument on board the Herschel Space Observatory. Observations of planets (including the prime point-source calibrator, Uranus), asteroids, line sources, dark sky, and cross-calibration sources were made in order to monitor repeatability and sensitivity, and to improve FTS calibration. We present a complete analysis of the full set of calibration observations and use them to assess the performance of the FTS. Particular care is taken to understand and separate out the effect of pointing uncertainties, including the position of the internal beam steering mirror for sparse observations in the early part of the mission. The repeatability of spectral line centre positions is <5km/s, for lines with signal-to-noise ratios >40, corresponding to <0.5-2.0% of a resolution element. For spectral line flux, the repeatability is better than 6%, which improves to 1-2% for spectra corrected for pointing offsets. The continuum repeatability is 4.4% for the SLW band and 13.6% for the SSW band, which reduces to ~1% once the data have been corrected for pointing offsets. Observations of dark sky were used to assess the sensitivity and the systematic offset in the continuum, both of which were found to be consistent across the FTS detector arrays. The average point-source calibrated sensitivity for the centre detectors is 0.20 and 0.21 Jy [1 sigma; 1 hour], for SLW and SSW. The average continuum offset is 0.40 Jy for the SLW band and 0.28 Jy for the SSW band.

We studied the radio emission from four radio-loud and gamma-ray-loud narrow-line Seyfert 1 galaxies. The goal was to investigate whether a relativistic jet is operating at the source, and quantify its characteristics. We relied on the most systematic monitoring of such system in the cm and mm radio bands which is conducted with the Effelsberg 100 m and IRAM 30 m telescopes and covers the longest time-baselines and the most radio frequencies to date. We extract variability parameters and compute variability brightness temperatures and Doppler factors. The jet powers were computed from the light curves to estimate the energy output. The dynamics of radio spectral energy distributions were examined to understand the mechanism causing the variability. All the sources display intensive variability that occurs at a pace faster than what is commonly seen in blazars. The flaring events show intensive spectral evolution indicative of shock evolution. The brightness temperatures and Doppler factors are moderate, implying a mildly relativistic jet. The computed jet powers show very energetic flows. The radio polarisation in one case clearly implies a quiescent jet underlying the recursive flaring activity. Despite the generally lower flux densities, the sources appear to show all typical characteristics seen in blazars that are powered by relativistic jets.

The Herschel SPIRE instrument consists of an imaging photometric camera and an imaging Fourier Transform Spectrometer (FTS), both operating over a frequency range of 450-1550 GHz. In this paper, we briefly review the FTS design, operation, and data reduction, and describe in detail the approach taken to relative calibration (removal of instrument signatures) and absolute calibration against standard astronomical sources. The calibration scheme assumes a spatially extended source and uses the Herschel telescope as primary calibrator. Conversion from extended to point-source calibration is carried out using observations of the planet Uranus. The model of the telescope emission is shown to be accurate to within 6% and repeatable to better than 0.06% and, by comparison with models of Mars and Neptune, the Uranus model is shown to be accurate to within 3%. Multiple observations of a number of point-like sources show that the repeatability of the calibration is better than 1%, if the effects of the satellite absolute pointing error (APE) are corrected. The satellite APE leads to a decrement in the derived flux, which can be up to ~10% (1 sigma) at the high-frequency end of the SPIRE range in the first part of the mission, and ~4% after Herschel operational day 1011. The lower frequency range of the SPIRE band is unaffected by this pointing error due to the larger beam size. Overall, for well-pointed, point-like sources, the absolute flux calibration is better than 6%, and for extended sources where mapping is required it is better than 7%.

The Fourier Transform Spectrometer (FTS) of the Spectral and Photometric Imaging REceiver (SPIRE) on board the ESA Herschel Space Observatory has two detector setting modes: (a) a nominal mode, which is optimized for observing moderately bright to faint astronomical targets, and (b) a bright-source mode recommended for sources significantly brighter than 500 Jy, within the SPIRE FTS bandwidth of 446.7-1544 GHz (or 194-671 microns in wavelength), which employs a reduced detector responsivity and out-of-phase analog signal amplifier/demodulator. We address in detail the calibration issues unique to the bright-source mode, describe the integration of the bright-mode data processing into the existing pipeline for the nominal mode, and show that the flux calibration accuracy of the bright-source mode is generally within 2% of that of the nominal mode, and that the bright-source mode is 3 to 4 times less sensitive than the nominal mode.

We present a method to derive the relative pointing offsets for SPIRE Fourier-Transform Spectrometer (FTS) solar system object (SSO) calibration targets, which were observed regularly throughout the Herschel mission. We construct ratios of the spectra for all observations of a given source with respect to a reference. The reference observation is selected iteratively to be the one with the highest observed continuum. Assuming that any pointing offset leads to an overall shift of the continuum level, then these ratios represent the relative flux loss due to mispointing. The mispointing effects are more pronounced for a smaller beam, so we consider only the FTS short wavelength array (SSW, 958-1546 GHz) to derive a pointing correction. We obtain the relative pointing offset by comparing the ratio to a grid of expected losses for a model source at different distances from the centre of the beam, under the assumption that the SSW FTS beam can be well approximated by a Gaussian. In order to avoid dependency on the point source flux conversion, which uses a particular observation of Uranus, we use extended source flux calibrated spectra to construct the ratios for the SSOs. In order to account for continuum variability, due to the changing distance from the Herschel telescope, the SSO ratios are normalised by the expected model ratios for the corresponding observing epoch. We confirm the accuracy of the derived pointing offset by comparing the results with a number of control observations, where the actual pointing of Herschel is known with good precision. Using the method we derived pointing offsets for repeated observations of Uranus (including observations centred on off-axis detectors), Neptune, Ceres and NGC7027. The results are used to validate and improve the point-source flux calibration of the FTS.

Emission from the Herschel telescope is the dominant source of radiation for the majority of SPIRE Fourier transform spectrometer (FTS) observations, despite the exceptionally low emissivity of the primary and secondary mirrors. Accurate modelling and removal of the telescope contribution is, therefore, an important and challenging aspect of FTS calibration and data reduction pipeline. A dust-contaminated telescope model with time invariant mirror emissivity was adopted before the Herschel launch. However, measured FTS spectra show a clear evolution of the telescope contribution over the mission and strong need for a correction to the standard telescope model in order to reduce residual background (of up to 7 Jy) in the final data products. Systematic changes in observations of dark sky, taken over the course of the mission, provide a measure of the evolution between observed telescope emission and the telescope model. These dark sky observations have been used to derive a time dependent correction to the telescope emissivity that reduces the systematic error in the continuum of the final FTS spectra to ~0.35 Jy.

The Herschel SPIRE Fourier transform spectrometer (FTS) performs spectral imaging in the 447-1546 GHz band. It can observe in three spatial sampling modes: sparse mode, with a single pointing on sky, or intermediate or full modes with 1 and 1/2 beam spacing, respectively. In this paper, we investigate the uncertainty and repeatability for fully sampled FTS mapping observations. The repeatability is characterised using nine observations of the Orion Bar. Metrics are derived based on the ratio of the measured intensity in each observation compared to that in the combined spectral cube from all observations. The mean relative deviation is determined to be within 2%, and the pixel-by-pixel scatter is ~7%. The scatter increases towards the edges of the maps. The uncertainty in the frequency scale is also studied, and the spread in the line centre velocity across the maps is found to be ~15 km/s. Other causes of uncertainty are also discussed including the effect of pointing and the additive uncertainty in the continuum.

Herschel/SPIRE Fourier transform spectrometer (FTS) observations contain emission from both the Herschel Telescope and the SPIRE Instrument itself, both of which are typically orders of magnitude greater than the emission from the astronomical source, and must be removed in order to recover the source spectrum. The effects of the Herschel Telescope and the SPIRE Instrument are removed during data reduction using relative spectral response calibration curves and emission models. We present the evolution of the methods used to derive the relative spectral response calibration curves for the SPIRE FTS. The relationship between the calibration curves and the ultimate sensitivity of calibrated SPIRE FTS data is discussed and the results from the derivation methods are compared. These comparisons show that the latest derivation methods result in calibration curves that impart a factor of between 2 and 100 less noise to the overall error budget, which results in calibrated spectra for individual observations whose noise is reduced by a factor of 2-3, with a gain in the overall spectral sensitivity of 23% and 21% for the two detector bands, respectively.

The quasar 1156+295 (4C +29.45) is one of the targets in the Urumqi monitoring program which aimed to search for evidence of annual modulation in the timescales of intra-day variable (IDV) sources. The IDV observations of 1156+295 were carried out nearly monthly from October 2007 to October 2009, with the Urumqi 25m radio telescope at 4.8 GHz. The source has shown prominent IDV of total flux density in most observing sessions with variability timescales of <= 1 day at 4.8 GHz. The estimated IDV timescales seem to follow an annual cycle that can be fitted with an anisotropic interstellar scintillation (ISS) model suggesting that a significant part of the flux density variations is due to ISS. The source underwent a dramatic flare in 2008. We studied the possible consequences of the flare on the IDV of 1156+295 by comparing the changes in its IDV characteristics with the evolution of the 43 GHz Very Long Baseline Array (VLBA) core size of the source. Conclusions: The quasar 1156+295 shows evidence for an annual modulation of its IDV timescales at 4.8 GHz, the ISS-induced IDV timescales and variability strength might be affected by the overall activity state of the source core. More frequent IDV and VLBI measurements are required to confirm the relation between the IDV appearance and the core-size evolution of the source.

The detection of gamma rays from a small number of Narrow Line Seyfert 1 galaxies by the LAT instrument onboard Fermi seriously challenged our understanding of AGN physics. Among the most important findings associated with their discovery has been the realisation that smaller-mass black holes seem to be hosted by these systems. Immediately after their discovery a radio multi- frequency monitoring campaign was initiated to understand their jet radio emission. Here the first results of the campaign are presented. The light curves and some first variability analyses are discussed, showing that the brightness temperatures and Doppler factors are moderate. The phenomenologies are typically blazar-like. The frequency domain on the other hand indicates intense spectral evolution and the variability patterns indicate mechanisms similar to those acting in the jets of BL Lacs and FSRQs. Finally, the linear polarisation also reveals the presence of a quiescent, optically thin jet in certain cases.

We present the results of a series of radio, optical, X-ray and gamma-ray observations of the BL Lac object S50716+714 carried out between April 2007 and January 2011. The multi-frequency observations were obtained using several ground and space based facilities. The intense optical monitoring of the source reveals faster repetitive variations superimposed on a long-term variability trend at a time scale of ~350 days. Episodes of fast variability recur on time scales of ~ 60-70 days. The intense and simultaneous activity at optical and gamma-ray frequencies favors the SSC mechanism for the production of the high-energy emission. Two major low-peaking radio flares were observed during this high optical/gamma-ray activity period. The radio flares are characterized by a rising and a decaying stage and are in agreement with the formation of a shock and its evolution. We found that the evolution of the radio flares requires a geometrical variation in addition to intrinsic variations of the source. Different estimates yield a robust and self-consistent lower limits of \delta > 20 and equipartition magnetic field B_eq > 0.36 G. Causality arguments constrain the size of emission region \theta < 0.004 mas. We found a significant correlation between flux variations at radio frequencies with those at optical and gamma-rays. The optical/GeV flux variations lead the radio variability by ~65 days. The longer time delays between low-peaking radio outbursts and optical flares imply that optical flares are the precursors of radio ones. An orphan X-ray flare challenges the simple, one-zone emission models, rendering them too simple. Here we also describe the spectral energy distribution modeling of the source from simultaneous data taken through different activity periods.

We present a comprehensive data description for Ks-band measurements of Sgr A*. We characterize the statistical properties of the variability of Sgr A* in the near-infrared, which we find to be consistent with a single-state process forming a power-law distribution of the flux density. We discover a linear rms-flux relation for the flux-density range up to 12 mJy on a timescale of 24 minutes. This and the power-law flux density distribution implies a phenomenological, formally non-linear statistical variability model with which we can simulate the observed variability and extrapolate its behavior to higher flux levels and longer timescales. We present reasons why data with our cadence cannot be used to decide on the question whether the power spectral density of the underlying random process shows more structure at timescales between 25 min and 100 min compared to what is expected from a red noise random process.

We report results from a 1 week multi-wavelength campaign to monitor the BL Lac object S5 0716+714 (on December 9-16, 2009). In the radio bands the source shows rapid (~ (0.5-1.5) day) intra-day variability with peak amplitudes of up to ~ 10 %. The variability at 2.8 cm leads by about 1 day the variability at 6 cm and 11 cm. This time lag and more rapid variations suggests an intrinsic contribution to the source's intraday variability at 2.8 cm, while at 6 cm and 11 cm interstellar scintillation (ISS) seems to predominate. Large and quasi-sinusoidal variations of ~ 0.8 mag were detected in the V, R and I-bands. The X-ray data (0.2-10 keV) do not reveal significant variability on a 4 day time scale, favoring reprocessed inverse-Compton over synchrotron radiation in this band. The characteristic variability time scales in radio and optical bands are similar. A quasi-periodic variation (QPO) of 0.9 - 1.1 days in the optical data may be present, but if so it is marginal and limited to 2.2 cycles. Cross-correlations between radio and optical are discussed. The lack of a strong radio-optical correlation indicates different physical causes of variability (ISS at long radio wavelengths, source intrinsic origin in the optical), and is consistent with a high jet opacity and a compact synchrotron component peaking at ~= 100 GHz in an ongoing very prominent flux density outburst. For the campaign period, we construct a quasi-simultaneous spectral energy distribution (SED), including gamma-ray data from the FERMI satellite. We obtain lower limits for the relativistic Doppler-boosting of delta >= 12-26, which for a BL\,Lac type object, is remarkably high.

The BLLac object S4 0954+65 is one of the main targets of the Urumqi monitoring program targeting IntraDay Variable (IDV) sources. Between August 2005 and December 2009, the source was included in 41 observing sessions, carried out at a frequency of 4.8 GHz. The time analysis of the collected light curves, performed by applying both a structure function analysis and a specifically developed wavelet-based algorithm, discovered an annual cycle in the variability timescales, suggesting that there is a fundamental contribution by interstellar scintillation to the IDV pattern of the source. The combined use of the two analysis methods also revealed that there was a dramatic change in the variability characteristics of the source between February and March 2008, at the starting time of a strong outburst phase. The analysis' results suggest that the flaring state of the source coincides with the appearance of multiple timescales in its light curves, indicating that changes in the structure of the relativistically moving emitting region may strongly influence the variability observed on IDV timescales.

We aim to search for evidence of annual modulation in the time scales of the BL Lac object S5 0716+714. The intra-day variability (IDV) observations were carried out monthly from 2005 to 2009, with the Urumqi 25m radio telescope at 4.8 GHz. The source has shown prominent IDV as well as long-term flux variations. The IDV time scale does show evidence in favor of an annual modulation, suggesting that the IDV of 0716+714 is dominated by interstellar scintillation. The source underwent a strong outburst phase between mid-2008 and mid-2009; a second intense flare was observed in late 2009, but no correlation between the total flux density and the IDV time scale is found, implying that the flaring state of the source does not have serious implications for the general characteristics of its intra-day variability. However, we find that the inner-jet position angle is changing throughout the years, which could result in an annual modulation noise in the anisotropic ISS model fit. There is also an indication that the lowest IDV amplitudes (rms flux density) correspond to the slowest time scales of IDV, which would be consistent with an ISS origin of the IDV of 0716+714.

The F-GAMMA program is among the most comprehensive programs that aim at understanding the physics in active galactic nuclei through the multi-frequency monitoring of Fermi blazars. Here we discuss monthly sampled broad-band radio spectra (2.6 - 142 GHz). Two different studies are presented. (a) We discuss that the variability patterns traced can be classified into two classes: (1) to those showing intense spectral-evolution and (2) those showing a self-similar quasi-achromatic behaviour. We show that a simple two-component model can very well reproduce the observed phenomenologies. (b) We present the cm-to-mm behaviour of three gamma-ray bright Narrow Line Seyfert 1 galaxies over time spans varying between ~1.5 and 3 years and compare their variability characteristics with typical blazars.

The F-GAMMA program is a coordinated effort to investigate the physics of Active Galactic Nuclei (AGNs) via multi-frequency monitoring of Fermi blazars. In the current study we show and discuss the evolution of broad-band radio spectra, which are measured at ten frequencies between 2.64 and 142 GHz using the Effelsberg 100-m and the IRAM 30-m telescopes. It is shown that any of the 78 sources studied can be classified in terms of their variability characteristics in merely 5 types of variability. It is argued that these can be attributed to only two classes of variability mechanisms. The first four types are dominated by spectral evolution and can be described by a simple two-component system composed of: (a) a steep quiescent spectral component from a large scale jet and (b) a time evolving flare component following the "Shock-in-Jet" evolutionary path. The fifth type is characterised by an achromatic change of the broad band spectrum, which could be attributed to a different mechanism, likely involving differential Doppler boosting caused by geometrical effects. Here we present the classification, the assumed physical scenario and the results of calculations that have been performed for the spectral evolution of flares.

The F-GAMMA program is a coordinated effort to investigate the physics of Active Galactic Nuclei (AGNs) via multi-frequency monitoring of {\em Fermi} blazars. The current study is concerned with the broad-band radio spectra composed of measurement at ten frequencies between 2.64 and 142 GHz. It is shown that any of the 78 sources studied can be classified in terms of their variability characteristics in merely 5 types of variability. The first four types are dominated by spectral evolution and can be reproduced by a simple two-component system made of the quiescent spectrum of a large scale jet populated with a flaring event evolving according to Marscher & Gear (1985). The last type is characterized by an achromatic change of the broad-band spectrum which must be attributed to a completely different mechanism. Here are presented, the classification, the assumed physical system and the results of simulations that have been conducted.

Using the Urumqi 25m radio telescope, sources from the first three-month Fermi-LAT detected AGN catalog with declination >0 were observed in 2009 at 4.8 GHz. The radio flux density appears to correlate with the gamma-ray intensity. Intra-day variability (IDV) observations were performed in March, April and May in 2009 for selected 42 gamma-ray bright blazars, and 60% of them show evident flux variability at 4.8 GHz during the IDV observations, the IDV detection rate is higher than that in previous flat-spectrum AGN samples. The IDV appears more often in the VLBI-core dominant blazars, and the non-IDV blazars show relatively `steeper' spectral indices than the IDV blazars. Pronounced inter-month variability has been found in two BL Lac objects: J0112+2244 and J0238+1636.

Starting from December 2004, a program for the monitoring of intraday variable sources at a frequency of 5 GHz was performed at the Urumqi Observatory. The analysis of the variability characteristics of the flat-spectrum radio source AO 0235+164 revealed the existence of an annual cycle in the variability amplitude. This appears to correlate with the solar elongation of the source. A thorough analysis of the results of the MASIV IDV survey --- which provides the variability characteristics of a large sample of compact radio sources --- confirms that there is a small but detectable component of the observed fractional modulation which increases with decreasing solar elongation. We discuss the hypothesis that the phenomenon is related to interplanetary scintillation.

Blazars from the first-three-months Fermi-AGN list were observed with the Urumqi 25m radio telescope at 5GHz in IDV (Intra-Day Variability) mode and inter-month observation mode. A significant correlation between the flux density at 5GHz and the gamma-ray intensity for the Fermi-LAT detected blazars is seen. There is a higher IDV detection rate in Fermi detected blazars than those reported for other samples. Stronger variability appears at lower Galactic latitudes; IDV appears to be stronger in weaker sources, indicating that the variability is affected by interstellar scintillation.

The F -GAMMA program is a coordinated effort of several observing facilities to understand the AGN/blazar phenomenon via a multi-frequency monitoring approach, especially in the era of Fermi-GST. Some 60 prominent sources are monitored monthly with the Effelsberg 100-m telescope, the IRAM 30-m telescope and more frequently but in a less uniform fashion at the APEX 12-m telescope, covering from 2.64 to 345 GHz. The program has been running since January 2007 and here some of its findings are summarized. (a) There are two major variability patterns that the spectra of sources follow, one spectral-evolution-dominated and one achromatic. (b) The FSRQs show higher brightness temperatures indicative of larger Doppler factors at play and (c) a statistically significant radio-gamma-ray correlation has been found with a method recently suggested by Pavlidou et al. (in prep.).

Following the recent discovery of gamma rays from the radio-loud narrow-line Seyfert 1 galaxy PMN J0948+0022 (z=0.5846), we started a multiwavelength campaign from radio to gamma rays, which was carried out between the end of March and the beginning of July 2009. The source displayed activity at all the observed wavelengths: a general decreasing trend from optical to gamma-ray frequencies was followed by an increase of radio emission after less than two months from the peak of the gamma-ray emission. The largest flux change, about a factor of about 4, occurred in the X-ray band. The smallest was at ultraviolet and near-infrared frequencies, where the rate of the detected photons dropped by a factor 1.6-1.9. At optical wavelengths, where the sampling rate was the highest, it was possible to observe day-scale variability, with flux variations up to a factor of about 3. The behavior of PMN J0948+0022 observed in this campaign and the calculated power carried out by its jet in the form of protons, electrons, radiation and magnetic field are quite similar to that of blazars, specifically of flat-spectrum radio quasars. These results confirm the idea that radio-loud narrow-line Seyfert 1 galaxies host relativistic jets with power similar to that of average blazars.

Short time-scale flux density variations of flat spectrum radio sources are often explained by the scattering of radio waves in the turbulent, ionized Interstellar Matter of the Milky Way. One of the most convincing observational arguments in favor of this is the annual modulation of the variability time-scale caused by the Earth orbital motion around the Sun. J1128+592 is an IDV source with a possible annual modulation in its variability time-scale. We observed the source in 6 epochs with the VLBA at 5, 8 and 15 GHz in total intensity and polarization. The VLBA observations revealed an east-west oriented core-jet structure. Its position angle agrees with the angle of anisotropy derived from the annual modulation model. No significant long-term structural changes were observed with VLBI on mas-scales, however, compared to archival data, the VLBI core size is expanded. This expansion offers a possible explanation to the observed decrease of the strength of IDV. VLBI polarimetry revealed significant changes in the electric vector position angle and Rotation Measure of the core and jet. Part of the observed RM variability could be attributed to a scattering screen (37 pc distance), which covers the source (core and jet) and which may be responsible for the IDV. Superposition of polarized sub-components below the angular resolution limit may affect the observed RM as well.

The F-GAMMA-project is the coordinated effort of several observatories to understand the AGN phenomenon and specifically blazars via multi-frequency monitoring in collaboration with the {\sl Fermi}-GST satellite since January 2007. The core observatories are: the Effelsberg 100-m, the IRAM 30-m and the OVRO 40-m telescope covering the range between 2.6 and 270 GHz. Effelsberg and IRAM stations do a monthly monitoring of the cm to sub-mm radio spectra of 60 selected blazars whereas the OVRO telescope is observing roughly 1200 objects at 15 GHz with a dense sampling of 2 points per week. The calibration uncertainty even at high frequencies, is of a few percent. 47% of the Effelsberg/Pico Veleta sample is included in the LBAS list. An update of the monitored sample is currently underway.

We report on a set of 5 GHz Urumqi observations of the galaxy M82, made between August 2005 and May 2009. From the resulting flux densities, we detect a strong flare, starting in March or April and peaking in June 2008. We identify this flare with supernova SN 2008iz. The time sampling of the radio light curve allows us to obtain information on the precursor mass-loss rate, the strength of the magnetic field in the radiating region, the explosion date, and the deceleration of the expanding shock. We also check the possible contribution of Synchrotron Self Absorption (SSA) to the radio light curve and compare our model with other observations of the supernova at 22 GHz.

The analysis of the spectral energy distribution variability at frequencies from radio to TeV is a powerful tool in the investigation of the dynamics, the physics and the structure evolution occurring in the most exotic flavour of active galaxies, the blazars. In particular, the presence of {\sl Fermi-GST} is providing a unique opportunity for such studies delivering $\gamma$-ray data of unprecedented quality. Here we introduce a monitoring program that has been running at the Effelsberg 100m telescope since January 2007, underpinning a broad multi-frequency collaboration of facilities that cover the band from radio to infrared. Sixty one selected blazars are observed monthly between 2.64 GHz and 43 GHz. The calibration accuracy is better than a few percent as it is demonstrated with some preliminary examples.

A sample of classical IntraDay Variable (IDV) and IDV candidate sources has been monitored with the Urumqi 25m telescope and the Effelsberg 100m telescope. Aim of the project is to investigate the origin of IntraDay Variability, a phenomenon which has been observed in about 30% of flat spectrum radio quasars. Simultaneous Effelsberg-Urumqi observations demonstrated that the Urumqi antenna, although relatively small in diameter, is well suitable for IDV experiments. A few Urumqi datasets, however, turned out to be affected by a spurious $\sim 24$ hours modulation, an effect which has been removed by means of a new procedure for data reduction. In about 14 months, 12 epochs of observation have been collected, for a total observing time of more than 45 days. The epochs are regularly distributed over the whole year, in order to check for the presence of systematic annual changes in the variability time scales - a crucial test for verifying the consistency of source-extrinsic models of the variability. Preliminary time-analysis of the monitored sources revealed some hint for a slowing down of the characteristic time scales of S5~0716+714, a result that, if confirmed, would be compatible with a source-extrinsic origin of the variability, in contrast to previous IDV studies. No significant modulation of the time scales has been detected for S4~0954+658.

Short time-scale radio variations of compact extragalactic radio quasars and blazars known as IntraDay Variability (IDV) can be explained in at least some sources as a propagation effect; the variations are interpreted as scintillation of radio waves in the turbulent interstellar medium of the Milky Way. One of the most convincing observational arguments in favor of a propagation-induced variability scenario is the observed annual modulation in the characteristic time scale of the variation due to the Earth's orbital motion. So far there are only two sources known with a well-constrained seasonal cycle. Annual modulation has been proposed for a few other less well-documented objects. However, for some other IDV sources source-intrinsic structural variations which cause drastic changes in the variability time scale were also suggested. J1128+592 is a recently discovered, highly variable IDV source. Previous, densely time-sampled flux-density measurements with the Effelsberg 100-m radio telescope (Germany) and the Urumqi 25-m radio telescope (China), strongly indicate an annual modulation of the time scale. The most recent 4 observations in 2006/7, however, do not fit well to the annual modulation model proposed before. In this paper, we investigate a possible explanation of this discrepancy.

Short time-scale radio variations of compact extragalactic radio sources, known as IntraDay Variability, can be explained in at least some sources by a source-extrinsic effect, in which the variations are interpreted as scintillation of radio waves caused by the turbulent ISM of the Milky Way. One of the most convincing observational arguments in favour of propagation-induced variability is the so called annual modulation of the characteristic variability time-scale, which is due to the orbital motion of the Earth. Data for the recently discovered and highly variable IDV source J1128+5925 are presented. We study the frequency and time dependence of the IDV in this compact quasar. We measure the characteristic variability time-scale of the IDV throughout the year, and analyze whether the observed changes in the variability time-scale are consistent with annual modulation. We monitored the flux density variability of J1128+5925 with dense time sampling between 2.7 and 10.45GHz with the 100m Effelsberg radio telescope of the MPIfR and with the 25m Urumqi radio telescope. From ten observing sessions, we determine the variability characteristics and time-scales. The observed pronounced changes of the variability time-scale of J1128+5925 are modelled with an anisotropic annual modulation model. The observed frequency dependence of the variation is in good agreement with the prediction from interstellar scintillation. Adopting a simple model for the annual modulation model and using also the frequency dependence of the IDV, we derive a lower limit to the distance of the scattering screen and an upper limit to the scintillating source size. The latter is found to be consistent with the measured core size from VLBI.

We present the results of a series of Swift and quasi simultaneous ground-based infra-red observations of the blazar 3C 454.3 carried out in April-May 2005 when the source was 10 to 30 times brighter than previously observed. We found 3C 454.3 to be very bright and variable at all frequencies covered by our instrumentation. The broad-band Spectral Energy Distribution (SED) shows the usual two-bump shape (in nu-nu f(nu) space) with the Infra-red, optical and UV data sampling the declining part of the synchrotron emission that, even during this extremely large outburst, had its maximum in the far-infrared. The X-ray spectral data from the XRT and BAT instruments are flat and due to inverse Compton emission. The remarkable SED observed implies that at the time of the Swift pointings 3C 454.3 was one of the brightest objects in the extragalactic sky with a gamma-ray emission similar or brighter than that of 3C 279 when observed in a high state by EGRET. Time variability in the optical-UV flux is very different from that in the X-ray data: while the first component varied by about a factor two within a single exposure, but remained approximately constant between different observations, the inverse Compton component did not vary on short time-scales but changed by more than a factor of 3 between observations separated by a few days. This different dynamical behaviour illustrates the need to collect simultaneous multi-frequency data over a wide range of time-scales to fully constrain physical parameters in blazars.

The flat-spectrum radio quasar 3C 454.3 is well known to be a highly active and variable source with outbursts occurring across the whole electromagnetic spectrum over the last decades. In spring 2005, 3C 454.3 has been reported to exhibit a strong optical outburst which subsequently triggered multi-frequency observations of the source covering the radio up to gamma-ray bands. Here, we present first results of our near-IR/optical (V, R, I, H band) photometry performed between May 11 and August 5, 2005 with the Rapid Eye Mount (REM) at La Silla in Chile and the Automatic Imaging Telescope (AIT) of the Perugia University Observatory. 3C 454.3 was observed during an exceptional and historical high state with a subsequent decrease in brightness over our 86 days observing period. The continuum spectral behaviour during the flaring and declining phase suggests a synchrotron peak below the near-IR band as well as a geometrical origin of the variations e.g. due to changes in the direction of forward beaming.

We present preliminary results of a statistical analysis obtained with a sample of blazars observed at the Perugia Astronomical Observatory since 1992. We briefly show how these statistical results can be useful to discriminate faint variable sources against the background noise. This technique, together with the more traditional ones, may be used to discover and identify high-energy point sources.