We report observations of a transient source \fermi\ from radio to \grs. \fermi\ was discovered by the {\it Fermi-LAT} in May 2017. Follow-up {\it Swift-XRT} observations revealed three flaring episodes through March 2018, and the peak X-ray flux is about $10^3$ higher than the {\it ROSAT all-sky survey (RASS)} flux upper limit. Optical spectral measurements taken by the {\it Magellan 6.5-m telescope} and the {\it Lick-Shane telescope} both show a largely featureless spectrum, strengthening the BL Lac interpretation first proposed by \citet{Bruni18}. The optical and mid-infrared (MIR) emission goes to a higher state in 2018, when the flux in high energies goes down to a lower level. Our {\it RATAN-600m} measurements at 4.8~GHz and 8.2~GHz do not indicate any significant radio flux variation over the monitoring seasons in 2017 and 2018, nor deviate from the archival {\it NVSS} flux level. During GeV flaring times, the spectrum is very hard ($\Gamma_\gamma\sim$1.7) in the GeV band and at times also very hard (($\Gamma_{\rm X}\lesssim2$) in the X-rays, similar to a high-synchrotron-peak (or even an extreme) BL Lac object, making \fermi\ a good target for ground-based {\it Cherenkov telescopes}.

We report the results of decade-long (2008-2018) $\gamma$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $\gamma$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $\gamma$-ray-optical flux-flux relation changes with activity state, ranging from a linear to a more complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz VLBA images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain $\gamma$-ray variability on very short time scales. The Mg II emission line flux in the `blue' and `red' wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands we find progressive delays of the most prominent light curve maxima with decreasing frequency, as expected from the frequency dependence of the $\tau=1$ surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at $\sim5$ GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet.

Gigahertz-Peaked spectrum (GPS) sources are compact active galactic nuclei, presumably young precursors of bright radio sources. The study of GPS radio properties provides information about the features of synchrotron radiation in extragalactic sources. Also in applied research GPS sources are useful as compact stationary radio sources in the sky for astrometric purposes. This paper presents the results of a multifrequency GPS study based on quasi-simultaneous measurements with the RATAN-600 radio telescope during the 2006-2017. The catalog containing spectral flux densities measured at six frequencies (1.1, 2.3, 4.8, 7.7/8.2, 11.2, and 21.7 GHz) have been obtained. In addition, for the analysis of radio spectra, data from the following low-frequency surveys have been used: GLEAM (GaLactic and Extragalactic All-sky Murchison widefield array survey) and TGSS (Tata institute for fundamental research GMRT Sky Survey) and high-frequency measurements from Planck survey. A total number of 164 GPS and candidates to GPS have been identified (17 of them are new discoveries), which makes up a small fraction of GPS in the initial sample of bright AGNs, about 2%. The physical properties and formation conditions of synchrotron radiation is found to be quite different in GPS of different AGN types. The deficit of distant GPS ($z > 2$) with low maximum frequencies (less than 1 GHz) is confirmed. The existing 'size - peak frequency' anticorrelation is continuous. The continuum radio spectra are found to become statistically steeper with increasing redshift.

We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and Fermi ($\gamma$ rays) satellites to study flux and spectral variability and the correlation between flux changes at different wavelengths. Unprecedented $\gamma$-ray flaring activity was observed during 2016 November-2017 February, with four major outbursts. A peak flux of (2158 $\pm$ 63)$\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, corresponding to a luminosity of (2.2 $\pm$ 0.1)$\times$10$^{50}$ erg s$^{-1}$, was reached on 2016 December 28. These four $\gamma$-ray outbursts have corresponding events in the near-infrared, optical, and UV bands, with the peaks observed at the same time. A general agreement between X-ray and $\gamma$-ray activity is found. The $\gamma$-ray flux variations show a general, strong correlation with the optical ones with no time lag between the two bands and a comparable variability amplitude. This $\gamma$-ray/optical relationship is in agreement with the geometrical model that has successfully explained the low-energy flux and spectral behaviour, suggesting that the long-term flux variations are mainly due to changes in the Doppler factor produced by variations of the viewing angle of the emitting regions. The difference in behaviour between radio and higher energy emission would be ascribed to different viewing angles of the jet regions producing their emission.

We present the results of nine years of the blazar observing programme at the RATAN-600 radio telescope (2005-2014). The data were obtained at six frequency bands (1.1, 2.3, 4.8, 7.7, 11.2, 21.7 GHz) for 290 blazars, mostly BL Lacs. In addition, we used data at 37 GHz obtained quasi-simultaneously with the Metsahovi radio observatory for some sources. The sample includes blazars of three types: high-synchrotron peaked (HSP), low-synchrotron peaked (LSP), and intermediate-synchrotron peaked (ISP). We present several epochs of flux density measurements, simultaneous radio spectra, spectral indices and properties of their variability. The analysis of the radio properties of different classes of blazars showed that LSP and HSP BL Lac blazars are quite different objects on average. LSPs have higher flux densities, flatter spectra and their variability increases as higher frequencies are considered. On the other hand, HSPs are very faint in radio domain, tend to have steep low frequency spectra, and they are less variable than LSPs at all frequencies. Another result is spectral flattening above 7.7 GHz detected in HSPs, while an average LSP spectrum typically remains flat at both the low and high frequency ranges we considered.

Observations at 37 GHz, performed at Mets\"ahovi Radio Observatory, are presented for a sample of 78 radio-loud and radio-quiet narrow-line Seyfert 1 (NLS1) galaxies, together with additional lower and higher frequency radio data from RATAN-600, Owens Valley Radio Observatory, and the Planck satellite. Most of the data have been gathered between February 2012 and April 2015 but for some sources even longer lightcurves exist. The detection rate at 37 GHz is around 19%, comparable to other populations of active galactic nuclei presumed to be faint at radio frequencies, such as BL Lac objects. Variability and spectral indices are determined for sources with enough detections. Based on the radio data, many NLS1 galaxies show a blazar-like radio spectra exhibiting significant variability. The spectra at a given time are often inverted or convex. The source of the high-frequency radio emission in NLS1 galaxies, detected at 37 GHz, is most probably a relativistic jet rather than star formation. Jets in NLS1 galaxies are therefore expected to be a much more common phenomenon than earlier assumed.

We have made an estimation of the synchrotron peak frequency ($\nu_{peak}^{s}$) for six very low synchrotron peaked (VLSP) blazars. These objects were selected as VLSP candidates (with the $\nu_{peak}^{s} \leq 10^{13}$ Hz) from the archival data. We have build spectral energy distribution using quasi-simultaneous observations at the Zeiss-1000 and RATAN-600 telescopes of the Special astrophysical observatory of RAS and made an estimation of the $\nu_{peak}^{s}$. We confirmed classification as a VLSP for two sources ([HB89] 1308+326 and 3C 345), for four other blazars we have calculated $\nu_{peak}^{s}>10^{13}$ Hz.

We study the synchrotron component of the spectral energy distribution (SED) on the sample of 877 blazars using the ASDC SED Builder Tool with available broadband data from the literature. Our sample includes 423 flat-spectrum radio sources (FSRQs), 361 BL Lac objects and candidates, and 93 blazars of uncertain type. We have made an estimation of the synchrotron peak frequency ($\nu_{peak}^{s}$) for the 875 objects and further classified them as high, intermediate and low synchrotron peaked sources (HSPs/ISPs/LSPs). There are 42 HSPs with $\nu_{peak}^{s} > 10^{16.5}$ Hz, 222 ISPs with $10^{14.5} < \nu_{peak}^{s} < 10^{16.5}$ Hz, and 611 LSPs with $\nu_{peak}^{s} < 10^{14.5}$ Hz in our sample. We have calculated an average value of $\nu_{peak}^{s}$ to be $10^{13.4 \pm 1.0}$ Hz for FSRQs and $10^{14.6 \pm 1.4}$ Hz for BL Lacs. We found out that $\nu_{peak}^{s}$ and the flux density at 4.8 GHz have a different distribution (as indicated by Kolmogorov-Smirnov test at significance level 0.05) for the FSRQ and BL Lac blazars, and for the RBL and XBL types of BL Lacs. Distribution of $\nu_{peak}^{s}$ values is broader for BL Lacs, than for FSRQs. There are no ultra-high energy peaked objects (with $\nu_{peak}^{s} > 10^{19}$ Hz) in our BL Lac sample according to our estimations. The significant part of FSRQs (41%) and small part of BL Lacs (9%) in our sample could be considered as candidates to the very-low synchrotron peaked blazars (with $\nu_{peak}^{s} < 10^{13}$ Hz). Our foundations confirm results of the previous studies made on samples with significantly smaller number of objects.

We study the correlation between gamma-ray and radio band radiation for 123 blazars, using the Fermi-LAT first source catalog (1FGL) and the RATAN-600 data obtained at the same period of time (within a few months). We found an apparent positive correlation for BL Lac and flat-spectrum radio quasar (FSRQ) sources from our sample through testing the value of the Pearson product-moment correlation coefficient. The BL Lac objects show higher values of the correlation coefficient than FSRQs at all frequencies, except 21.7 GHz, and at all bands, except $10-100$ GeV, typically at high confidence level (> 99%). At higher gamma-ray energies the correlation weakens and even becomes negative for BL Lacs and FSRQs. For BL Lac blazars, the correlation of the fluxes appeared to be more sensitive to the considered gamma-ray energy band, than to the frequency, while for FSRQ sources the correlation changed notably both with the considered radio frequency and gamma-ray energy band. We used a data randomization method to quantify the significance of the computed correlation coefficients. We find that the statistical significance of the correlations we obtained between the flux densities at all frequencies and the photon flux in all gamma-ray bands below 3 GeV is high for BL Lacs (chance probability $\sim 10^{-3} - 10^{-7}$). The correlation coefficient is high and significant for the $0.1-0.3$ GeV band and low and insignificant for the $10-100$ GeV band for both types of blazars for all considered frequencies.

We present a new catalogue of the RATAN-600 multi-frequency measurements for BL Lac objects. The purpose of this catalogue is to compile the BL Lac multi-frequency data that is acquired with the RATAN-600 simultaneously at several frequencies. The BL Lac objects emit a strongly variable and polarized non-thermal radiation across the entire electromagnetic spectrum from radio to gamma-rays and represent about one percent of known AGNs. They belong to the blazar population and differ from other blazars featureless optical spectrum, which sometimes have absorption lines, or have weak and narrow emission lines. One of the most effective ways of studying the physics of BL Lacs is the use of simultaneous multi-frequency data. The multi-frequency broadband radio spectrum was obtained simultaneously with an accuracy of up to 1-2 minutes for four to six frequencies: 1.1, 2.3, 4.8, 7.7, 11.2, and 21.7 GHz. The catalogue is based on the RATAN-600 observations and on the data from: equatorial coordinate and redshift, R-band magnitude, synchrotron peak frequency, SED classes and object type, literature. The present version of the catalogue contains RATAN-600 flux densities measurements over nine years (2006-2014), radio spectra at different epochs, and their parameters for more than 300 BL Lacs objects and candidates. The BL Lacs list is constantly updated with new observational data of RATAN-600.

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and gamma-ray bands. Our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi-LAT, whereas ~40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-mm spectral slope of blazars is quite flat up to ~70GHz, above which it steepens to <\alpha>~-0.65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (\nupS) in the SED of FSRQs is the same in all the blazar samples with <\nupS>=10^13.1 Hz, while the mean inverse-Compton peak frequency, <\nupIC>, ranges from 10^21 to 10^22 Hz. The distributions of \nupS and of \nupIC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs and strongly depend on the selection method. The Compton dominance of blazars ranges from ~0.2 to ~100, with only FSRQs reaching values >3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at ~7 or more, and radio-selected blazars at values ~1, thus implying that the assumption that the blazar power is dominated by high-energy emission is a selection effect. Simple SSC models cannot explain the SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and \nupS predicted by the blazar sequence.

Planck's all sky surveys at 30-857 GHz provide an unprecedented opportunity to follow the radio spectra of a large sample of extragalactic sources to frequencies 2-20 times higher than allowed by past, large area, ground-based surveys. We combine the results of the Planck Early Release Compact Source Catalog (ERCSC) with quasi-simultaneous ground-based observations, as well as archival data, at frequencies below or overlapping Planck frequency bands, to validate the astrometry and photometry of the ERCSC radio sources and study the spectral features shown in this new frequency window opened by Planck. The ERCSC source positions and flux density scales are found to be consistent with the ground-based observations. We present and discuss the spectral energy distributions (SEDs) of a sample of "extreme" radio sources to illustrate the richness of the ERCSC for the study of extragalactic radio sources. Variability is found to play a role in the unusual spectral features of some of these sources.

Spectral energy distributions (SEDs) and radio continuum spectra are presented for a northern sample of 104 extragalactic radio sources, based on the Planck Early Release Compact Source Catalogue (ERCSC) and simultaneous multifrequency data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous observations ranging from radio to gamma-rays. This is the first extensive frequency coverage in the radio and millimetre domains for an essentially complete sample of extragalactic radio sources, and it shows how the individual shocks, each in their own phase of development, shape the radio spectra as they move in the relativistic jet. The SEDs presented in this paper were fitted with second and third degree polynomials to estimate the frequencies of the synchrotron and inverse Compton (IC) peaks, and the spectral indices of low and high frequency radio data, including the Planck ERCSC data, were calculated. SED modelling methods are discussed, with an emphasis on proper, physical modelling of the synchrotron bump using multiplecomponents. Planck ERCSC data also suggest that the original accelerated electron energy spectrum could be much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shock. Furthermore in many cases the Planck data indicate that gamma-ray emission must originate in the same shocks that produce the radio emission.