The main characteristics in the radio continuum spectra of Supernova Remnants (SNRs) result from simple synchrotron emission. In addition, electron acceleration mechanisms can shape the spectra in specific ways, especially at high radio frequencies. These features are connected to the age and the peculiar conditions of the local interstellar medium interacting with the SNR. Whereas the bulk radio emission is expected at up to $20-50$ GHz, sensitive high-resolution images of SNRs above 10 GHz are lacking and are not easily achievable, especially in the confused regions of the Galactic Plane. In the framework of the early science observations with the Sardinia Radio Telescope in February-March 2016, we obtained high-resolution images of SNRs Tycho, W44 and IC443 that provided accurate integrated flux density measurements at 21.4 GHz: 8.8 $\pm$ 0.9 Jy for Tycho, 25 $\pm$ 3 Jy for W44 and 66 $\pm$ 7 Jy for IC443. We coupled the SRT measurements with radio data available in the literature in order to characterise the integrated and spatially-resolved spectra of these SNRs, and to find significant frequency- and region-dependent spectral slope variations. For the first time, we provide direct evidence of a spectral break in the radio spectral energy distribution of W44 at an exponential cutoff frequency of 15 $\pm$ 2 GHz. This result constrains the maximum energy of the accelerated electrons in the range $6-13$ GeV, in agreement with predictions indirectly derived from AGILE and \textit{Fermi}-LAT gamma-ray observations. With regard to IC443, our results confirm the noticeable presence of a bump in the integrated spectrum around $20-70$ GHz that could result from a spinning dust emission mechanism.

We present our broadband study of GW170817 from radio to hard X-rays, including NuSTAR and Chandra observations up to 165 days after the merger, and a multi-messenger analysis including LIGO constraints. The data are compared with predictions from a wide range of models, providing the first detailed comparison between non-trivial cocoon and jet models. Homogeneous and power-law shaped jets, as well as simple cocoon models are ruled out by the data, while both a Gaussian shaped jet and a cocoon with energy injection can describe the current dataset for a reasonable range of physical parameters, consistent with the typical values derived from short GRB afterglows. We propose that these models can be unambiguously discriminated by future observations measuring the post-peak behaviour, with slope -1.0 for the cocoon and -2.5 for the jet model.

On the wake of the results obtained so far at the SPARC\_LAB test-facility at the Laboratori Nazionali di Frascati (Italy), we are currently investigating the possibility to design and build a new multi-disciplinary user-facility, equipped with a soft X-ray Free Electron Laser (FEL) driven by a $\sim$1 GeV high brightness linac based on plasma accelerator modules. This design study is performed in synergy with the EuPRAXIA design study. In this paper we report about the recent progresses in the on going design study of the new facility.

A long-standing paradigm in astrophysics is that collisions- or mergers- of two neutron stars (NSs) form highly relativistic and collimated outflows (jets) powering gamma-ray bursts (GRBs) of short (< 2 s) duration. However, the observational support for this model is only indirect. A hitherto outstanding prediction is that gravitational wave (GW) events from such mergers should be associated with GRBs, and that a majority of these GRBs should be off-axis, that is, they should point away from the Earth. Here we report the discovery of the X-ray counterpart associated with the GW event GW170817. While the electromagnetic counterpart at optical and infrared frequencies is dominated by the radioactive glow from freshly synthesized r-process material in the merger ejecta, known as kilonova, observations at X-ray and, later, radio frequencies exhibit the behavior of a short GRB viewed off-axis. Our detection of X-ray emission at a location coincident with the kilonova transient provides the missing observational link between short GRBs and GWs from NS mergers, and gives independent confirmation of the collimated nature of the GRB emission.

Observations of supernova remnants (SNRs) are a powerful tool for investigating the later stages of stellar evolution, the properties of the ambient interstellar medium, and the physics of particle acceleration and shocks. For a fraction of SNRs, multi-wavelength coverage from radio to ultra high-energies has been provided, constraining their contributions to the production of Galactic cosmic rays. Although radio emission is the most common identifier of SNRs and a prime probe for refining models, high-resolution images at frequencies above 5 GHz are surprisingly lacking, even for bright and well-known SNRs such as IC443 and W44. In the frameworks of the Astronomical Validation and Early Science Program with the 64-m single-dish Sardinia Radio Telescope, we provided, for the first time, single-dish deep imaging at 7 GHz of the IC443 and W44 complexes coupled with spatially-resolved spectra in the 1.5-7 GHz frequency range. Our images were obtained through on-the-fly mapping techniques, providing antenna beam oversampling and resulting in accurate continuum flux density measurements. The integrated flux densities associated with IC443 are S_1.5GHz = 134 +/- 4 Jy and S_7GHz = 67 +/- 3 Jy. For W44, we measured total flux densities of S_1.5GHz = 214 +/- 6 Jy and S_7GHz = 94 +/- 4 Jy. Spectral index maps provide evidence of a wide physical parameter scatter among different SNR regions: a flat spectrum is observed from the brightest SNR regions at the shock, while steeper spectral indices (up to 0.7) are observed in fainter cooling regions, disentangling in this way different populations and spectra of radio/gamma-ray-emitting electrons in these SNRs.

[Abridged] The Sardinia Radio Telescope (SRT) is the new 64-m dish operated by INAF (Italy). Its active surface will allow us to observe at frequencies of up to 116 GHz. At the moment, three receivers, one per focal position, have been installed and tested. The SRT was officially opened in October 2013, upon completion of its technical commissioning phase. In this paper, we provide an overview of the main science drivers for the SRT, describe the main outcomes from the scientific commissioning of the telescope, and discuss a set of observations demonstrating the SRT's scientific capabilities. One of the main objectives of scientific commissioning was the identification of deficiencies in the instrumentation and/or in the telescope sub-systems for further optimization. As a result, the overall telescope performance has been significantly improved. As part of the scientific commissioning activities, different observing modes were tested and validated, and first astronomical observations were carried out to demonstrate the science capabilities of the SRT. In addition, we developed astronomer-oriented software tools, to support future observers on-site. The astronomical validation activities were prioritized based on technical readiness and scientific impact. The highest priority was to make the SRT available for joint observations as part of European networks. As a result, the SRT started to participate (in shared-risk mode) in EVN (European VLBI Network) and LEAP (Large European Array for Pulsars) observing sessions in early 2014. The validation of single-dish operations for the suite of SRT first light receivers and backends continued in the following years, and was concluded with the first call for shared-risk/early-science observations issued at the end of 2015.

We present high sensitivity polarimetric observations in 6 bands covering the 5.5-38 GHz range of a complete sample of 53 compact extragalactic radio sources brighter than 200 mJy at 20 GHz. The observations, carried out with the Australia Telescope Compact Array (ATCA), achieved a 91% detection rate (at 5 sigma). Within this frequency range the spectra of about 95% of sources are well fitted by double power laws, both in total intensity and in polarisation, but the spectral shapes are generally different in the two cases. Most sources were classified as either steep- or peaked-spectrum but less than 50% have the same classification in total and in polarised intensity. No significant trends of the polarisation degree with flux density or with frequency were found. The mean variability index in total intensity of steep-spectrum sources increases with frequency for a 4-5 year lag, while no significant trend shows up for the other sources and for the 8 year lag. In polarisation, the variability index, that could be computed only for the 8 year lag, is substantially higher than in total intensity and has no significant frequency dependence.

In the frame of the Astronomical Validation activities for the 64m Sardinia Radio Telescope, we performed 5-22 GHz imaging observations of the complex-morphology supernova remnants (SNRs) W44 and IC443. We adopted innovative observing and mapping techniques providing unprecedented accuracy for single-dish imaging of SNRs at these frequencies, revealing morphological details typically available only at lower frequencies through interferometry observations. High-frequency studies of SNRs in the radio range are useful to better characterize the spatially-resolved spectra and the physical parameters of different regions of the SNRs interacting with the ISM. Furthermore, synchrotron-emitting electrons in the high-frequency radio band are also responsible for the observed high-energy phenomenology as -e.g.- Inverse Compton and bremsstrahlung emission components observed in gamma-rays, to be disentangled from hadron emission contribution (providing constraints on the origin of cosmic rays).

We report the results of our observing campaign on GRB140903A, a nearby (z=0.351) short duration (T90~0.3 s) gamma-ray burst discovered by Swift. We monitored the X-ray afterglow with Chandra up to 21 days after the burst, and detected a steeper decay of the X-ray flux after approximately 1 day. Continued monitoring at optical and radio wavelengths showed a similar decay in flux at nearly the same time, and we interpret it as evidence of a narrowly collimated jet. By using the standard fireball model to describe the afterglow evolution, we derive a jet opening angle of 5 deg and a collimation-corrected total energy release of 2E50 erg. We further discuss the nature of the GRB progenitor system. Three main lines disfavor a massive star progenitor: the properties of the prompt gamma-ray emission, the age and low star-formation rate of the host galaxy, and the lack of a bright supernova. We conclude that this event was likely originated by a compact binary merger.

The machine described in this document is an advanced Source of up to 20 MeV Gamma Rays based on Compton back-scattering, i.e. collision of an intense high power laser beam and a high brightness electron beam with maximum kinetic energy of about 720 MeV. Fully equipped with collimation and characterization systems, in order to generate, form and fully measure the physical characteristics of the produced Gamma Ray beam. The quality, i.e. phase space density, of the two colliding beams will be such that the emitted Gamma ray beam is characterized by energy tunability, spectral density, bandwidth, polarization, divergence and brilliance compatible with the requested performances of the ELI-NP user facility, to be built in Romania as the Nuclear Physics oriented Pillar of the European Extreme Light Infrastructure. This document illustrates the Technical Design finally produced by the EuroGammaS Collaboration, after a thorough investigation of the machine expected performances within the constraints imposed by the ELI-NP tender for the Gamma Beam System (ELI-NP-GBS), in terms of available budget, deadlines for machine completion and performance achievement, compatibility with lay-out and characteristics of the planned civil engineering.

The present Report concerns the current status of the Italian Tau/Charm accelerator project and in particular discusses the issues related to the lattice design, to the accelerators systems and to the associated conventional facilities. The project aims at realizing a variable energy Flavor Factory between 1 and 4.6 GeV in the center of mass, and succeeds to the SuperB project from which it inherits most of the solutions proposed in this document. The work comes from a cooperation involving the INFN Frascati National Laboratories accelerator experts, the young newcomers, mostly engineers, of the Cabibbo Lab consortium and key collaborators from external laboratories.

We present polarisation data for 180 extragalactic sources extracted from the Australia Telescope 20 GHz (AT20G) survey catalog, and observed with the Australia Telescope Compact Array during a dedicated, high sensitivity run. For the sake of completeness we extracted the polarisation information for 7 extended sources from the 9-yr WMAP coadded maps at 23 GHz. The full sample of 187 sources constitutes a 99% complete sample of extragalactic sources brighter than S(20 GHz)=500mJy at the selection epoch with declination <-30 deg. The sample has a 91.4% detection rate in polarisation at 20GHz (94% if considering the sub-sample of point like sources). We have measurements also at 4.8 and 8.6 GHz within 1 month of the 20GHz observations for 172 sources to reconstruct the spectral properties of the sample in total intensity and in polarisation: 143 of them have a polarisation detection at all three frequencies. We find that there is no statistically significant evidence of a relationship either between the fraction of polarisation and frequency or between the fraction of polarisation and the total intensity flux density. This indicates that Faraday depolarisation is not very important above 4.8 GHz and that the magnetic field is not substantially more ordered in the regions dominating the emission at higher frequencies (up to 20 GHz). We estimate the distribution of the polarisation fraction and the polarised flux density source counts at 20GHz.

We present follow-up observations at 5, 8 and 30 GHz of the K-band Northern Wide Survey (KNoWS) 20 GHz Bright Sample, performed with the 32-m Medicina Radio Telescope and the 32-m Torun Radio Telescope. The KNoWS sources were selected in the Northern Polar Cap (Delta > 72 deg) and have a flux density limit S(20GHz) = 115 mJy. We include NVSS 1.4 GHz measurements to derive the source radio spectra between 1.4 and 30 GHz. Based on optical identifications, 68 per cent of the sources are QSOs, and 27 per cent are radio galaxies. A redshift measurement is available for 58 per cent of the sources. The radio spectral properties of the different source populations are found to be in agreement with those of other high-frequency selected samples.

This report describes the scientific aims and potentials as well as the preliminary technical design of IRIDE, an innovative tool for multi-disciplinary investigations in a wide field of scientific, technological and industrial applications. IRIDE will be a high intensity 'particle factory', based on a combination of a high duty cycle radio-frequency superconducting electron linac and of high energy lasers. Conceived to provide unique research possibilities for particle physics, for condensed matter physics, chemistry and material science, for structural biology and industrial applications, IRIDE will open completely new research possibilities and advance our knowledge in many branches of science and technology. IRIDE will contribute to open new avenues of discoveries and to address most important riddles: What does matter consist of? What is the structure of proteins that have a fundamental role in life processes? What can we learn from protein structure to improve the treatment of diseases and to design more efficient drugs? But also how does an electronic chip behave under the effect of radiations? How can the heat flow in a large heat exchanger be optimized? The scientific potential of IRIDE is far reaching and justifies the construction of such a large facility in Italy in synergy with the national research institutes and companies and in the framework of the European and international research. It will impact also on R&D work for ILC, FEL, and will be complementarity to other large scale accelerator projects. IRIDE is also intended to be realized in subsequent stages of development depending on the assigned priorities.

We present the high angular resolution catalogue for the Australia Telescope 20 GHz (AT20G) survey, using the high angular resolution 6-km antenna data at the baselines of ~ 4500 m of the Australia Telescope Compact Array (ATCA). We have used the data to produce the visibility catalogue that separates the compact Active Galactic Nuclei (AGNs) from the extended radio sources at the 0.15 arcsec angular scale, corresponding to the linear size scale of 1 kpc at redshifts higher than 0.7. We find the radio population at 20 GHz to be dominated by compact AGNs constituting 77% of the total sources in the AT20G. We introduce the visibility-spectra diagnostic plot, produced using the AT20G cross-matches with lower frequency radio surveys at 1 GHz (the NRAO-VLA Sky Survey (NVSS) and the Sydney University Molonglo Sky Survey (SUMSS)), that separates the 20 GHz population into distinct sub-populations of the compact AGNs, the compact steep-spectrum sources, the extended AGN-powered sources and extended flat-spectrum sources. The extended flat-spectrum sources include a local thermal emitting population of high latitude planetary nebulae and also gravitational lens and binary black hole candidates among the AGNs. We find a smooth transition in properties between the compact CSS sources and the AGN populations. The visibility catalogue, together with the main AT20G survey, provides an estimate of angular size scales for sources in the AT20G and an estimate of the flux arising from central cores of extended radio sources. The identification of the compact AGNs in the AT20G survey provides high quality calibrators for high frequency radio telescope arrays and VLBI observations.

During 2010-2011, the Medicina 32-m dish hosted the 7-feed 18-26.5 GHz receiver built for the Sardinia Radio Telescope, with the goal to perform its commissioning. This opportunity was exploited to carry out a pilot survey at 20 GHz over the area for {\delta} > + 72.3{\deg}. This paper describes all the phases of the observations, as they were performed using new hardware and software facilities. The map-making and source extraction procedures are illustrated. A customised data reduction tool was used during the follow-up phase, which produced a list of 73 confirmed sources down to a flux density of 115 mJy. The resulting catalogue, here presented, is complete above 200 mJy. Source counts are in agreement with those provided by the AT20G survey. This pilot activity paves the way to a larger project, the K-band Northern Wide Survey (KNoWS), whose final aim is to survey the whole Northern Hemisphere down to a flux limit of 50 mJy (5{\sigma}).

We use high angular resolution data, measured from visibility of sources at the longest baseline of 4500 m of the Australia Telescope Compact Array (ATCA), for the Australia Telescope 20 GHz (AT20G) survey to obtain angular size information for > 94% of AT20G sources. We confirm the previous AT20G result that due to the high survey frequency of 20 GHz, the source population is strongly dominated by compact sources (79%). At 0.15 arcseconds angular resolution limit, we show a very strong correlation between the compact and extended sources with flat and steep-spectrum sources respectively. Thus, we provide a firm physical basis for the traditional spectral classification into flat and steep-spectrum sources to separate compact and extended sources. We find the cut-off of -0.46 to be optimum for spectral indices between 1 and 5 GHz and, hence, recommend the continued use of -0.5 for future studies. We study the effect of spectral curvature on redshift cut-off of compact AGNs using recently published redshift data. Using spectral indices at different frequencies, we correct for the redshift effect and produce restframe frequency spectra for compact sources for redshift up to ~5. We show that the flat spectra of most compact sources start to steepen at ~30 GHz. At higher frequencies, the spectra of both populations are steep so the use of spectral index does not separate the compact and extended source populations as well as in lower frequencies. We find that due to the spectral steepening, surveys of compact sources at higher frequencies (>5 GHz) will have redshift cut-off due to spectral curvature but at lower frequencies, the surveys are not significantly affected by spectral curvature, thus, the evidence for a strong redshift cut-off in AGNs found in lower frequency surveys is a real cut-off and not a result of K-correction.

We use AdS_5 x S^5 superstring sigma model perturbation theory to compute the leading one-loop corrections to the dispersion relations of the excitations near a long spinning string in AdS. This investigation is partially motivated by the OPE-based approach to the computation of the expectation value of null polygonal Wilson loops suggested in arXiv:1006.2788. Our results are in partial agreement with the recent asymptotic Bethe ansatz computation in arXiv:1010.5237. In particular, we find that the heaviest AdS mode (absent in the ABA approach) is stable and has a corrected one-loop dispersion relation similar to the other massive modes. Its stability might hold also at the next-to-leading order as we suggest using a unitarity-based argument.

We continue the investigation of two-loop string corrections to the energy of a folded string with a spin S in AdS_5 and an angular momentum J in S^5, in the scaling limit of large J and S with ell=pi J/(lambda^(1/2) ln S)=fixed. We compute the generalized scaling function at two-loop order f_2(ell) both for small and large values of ell matching the predictions based on the asymptotic Bethe ansatz. In particular, in the small ell expansion, we derive an exact integral form for the ell-dependent coefficient of the Catalan's constant term in f_2(ell). Also, by resumming a certain subclass of multi-loop Feynman diagrams we obtain an exact expression for the leading (ln ell) part of f(lambda^(1/2), ell) which is valid to any order in the alpha'~1/lambda^(1/2) expansion. At large ell the string energy has a BMN-like expansion and the first few leading coefficients are expected to be the same at weak and at strong coupling. We provide a new example of this non-renormalization for the term which is generated at two loops in string theory and at one-loop in gauge theory (sub-sub-leading in 1/J). We also derive a simple algebraic formula for the term of maximal transcendentality in f_2(ell) expanded at large ell. In the second part of the paper we initiate the study of 2-loop finite size corrections to the string energy by formally compactifying the spatial world-sheet direction in the string action expanded near long fast-spinning string. We observe that the leading finite-size corrections are of "Casimir" type coming from terms containing at least one massless propagator. We consider in detail the one-loop order (reproducing the leading Landau-Lifshitz model prediction) and then focus on the two-loop contributions to the (1/ln S) term (for J=0). We find that in a certain regularization scheme used to discard power divergences the two-loop coefficient of the (1/ln S) term appears to vanish.

We revisit the computation of the 2-loop correction to the energy of a folded spinning string in AdS_5 with an angular momentum J in S^5 in the scaling limit log S, J >>1 with J / log S fixed. This correction gives the third term in the strong-coupling expansion of the generalized scaling function. The computation, using the AdS light-cone gauge approach developed in our previous paper, is done by expanding the AdS_5 x S^5 superstring partition function near the generalized null cusp world surface associated to the spinning string solution. The result corrects and extends the previous conformal gauge result of arXiv:0712.2479 and is found to be in complete agreement with the corresponding terms in the generalized scaling function as obtained from the asymptotic Bethe ansatz in arXiv:0805.4615 (and also partially from the quantum O(6) model and the Bethe ansatz data in arXiv:0809.4952). This provides a highly nontrivial strong coupling comparison of the Bethe ansatz proposal with the quantum AdS_5 x S^5 superstring theory, which goes beyond the leading semiclassical term effectively controlled by the underlying algebraic curve. The 2-loop computation we perform involves all the structures in the AdS light-cone gauge superstring action of hep-th/0009171 and thus tests its ultraviolet finiteness and, through the agreement with the Bethe ansatz, its quantum integrability. We do most of the computations for a generalized spinning string solution or the corresponding null cusp surface that involves both the orbital momentum and the winding in a large circle of S^5.

We consider the AdS_5 x S^5 superstring in the light-cone gauge adapted to a massless geodesic in AdS5 in the Poincare patch. The resulting action has a relatively simple structure which makes it a natural starting point for various perturbative quantum computations. We illustrate the utility of this AdS light-cone gauge action by computing the 1-loop and 2-loop corrections to the null cusp anomalous dimension reproducing in a much simpler and efficient way earlier results obtained in conformal gauge. This leads to a further insight into the structure of the superstring partition function in non-trivial background.

[Abridged] Physical and evolutionary properties of the sub-mJy radio population are not entirely known. The radio/optical analysis of the ATESP 5 GHz sample has revealed a significant class of compact flat/inverted radio-spectrum sources associated to early-type galaxies up to redshift 2. Such sources are most plausibly triggered by an AGN, but their observational properties are not entirely consistent with those of standard radio galaxy populations. In the present work we aim at a better understanding of the radio spectra of such sources and ultimately of the nature of AGNs at sub-mJy flux levels. We used the ATCA to get multi-frequency (4.8, 8.6 and 19 GHz) quasi-simultaneous observations for a representative sub-sample of ATESP radio sources associated with early-type galaxies (26 objects with S>0.6 mJy). This can give us insight into the accretion/radiative mechanism that is at work, since different regimes display different spectral signatures in the radio domain. From the analysis of the radio spectra, we find that our sources are most probably jet-dominated systems. ADAF models are ruled out by the high frequency data, while ADAF+jet scenarios are still consistent with flat/moderately inverted-spectrum sources, but are not required to explain the data. We compared our sample with high (>20 GHz) frequency selected surveys, finding spectral properties very similar to the ones of much brighter (S>500 mJy) radio galaxies extracted from the Massardi et al. (2008) sample. Linear sizes of ATESP 5 GHz sources associated with early type galaxies are also often consistent with the ones of brighter B2 and 3C radio galaxies, with possibly a very compact component that could be associated at least in part to (obscured) radio-quiet quasar-like objects and/or low power BL Lacs.

We present wide-field imaging and polarimetry at 20GHz of seven of the most extended, bright (Stot >= 0.50 Jy), high-frequency selected radio sources in the southern sky with declinations < -30 deg. Accompanying the data are brief reviews of the literature for each source, The results presented here aid in the statistical completeness of the Australia Telescope 20GHz Survey's bright source sample. The data are of crucial interest for future cosmic microwave background missions as a collection of information about candidate calibrator sources. We are able to obtain data for seven of the nine sources identified by our selection criteria. We report that Pictor A is thus far the best extragalactic calibrator candidate for the Low Frequency Instrument of the Planck European Space Agency mission due to its high level of integrated polarized flux density (0.50+/-0.06 Jy) on a scale of 10 arcmin. Six of the seven sources have a clearly detected compact radio core, with either a null or less than two percent detection of polarized emission from the nucleus. Most sources with detected jets have magnetic field alignments running in a longitudinal configuration, however PKS1333-33 exhibits transverse fields and an orthogonal change in field geometry from nucleus to jets.

We present results from simulations of the extragalactic polarized sky at 1.4 GHz. As the basis for our polarization models, we use a semi-empirical simulation of the extragalactic total intensity (Stokes I) continuum sky developed at the University of Oxford (http://scubed.physics.ox.ac.uk) under the European SKA Design Study (SKADS) initiative, and polarization distributions derived from analysis of polarization observations. By considering a luminosity dependence for the polarization of AGN, we are able to fit the 1.4 GHz polarized source counts derived from the NVSS and the DRAO ELAIS N1 deep field survey down to approximately 1 mJy. This trend is confirmed by analysis of the polarization of a complete sample of bright AGN. We are unable to fit the additional flattening of the polarized source counts from the deepest observations of the ELAIS N1 survey, which go down to ~0.5 mJy. Below 1 mJy in Stokes I at 1.4 GHz, starforming galaxies become an increasingly important fraction of all radio sources. We use a spiral galaxy integrated polarization model to make realistic predictions of the number of polarized sources at microJy levels in polarized flux density and hence, realistic predictions of what the next generation radio telescopes such as ASKAP, other SKA pathfinders and the SKA itself will see.

We discuss our ATCA 18.5 and 22 GHz flux density measurements of the Southern extragalactic sources in the complete 5-GHz sample of Kuehr et al. (1981). The high-frequency (5-18.5 GHz) spectral indices of steep-spectrum sources for which we have 18.5-GHz data (66% of the complete sample) are systematically steeper than the low-frequency (2.7-5 GHz) ones, and there is evidence of an anti-correlation of high-frequency spectral index with luminosity. The completeness of the 18.5-GHz data is much higher (89%) for flat-spectrum sources (mostly quasars), which also exhibit a spectral steepening. Taking advantage of the almost complete redshift information on flat-spectrum quasars, we have estimated their 5-GHz luminosity function in several redshift bins. The results confirm that their radio luminosity density peaks at z_peak \simeq 2.5 but do not provide evidence for deviations from pure luminosity evolution as hinted at by other data sets. A comparison of our 22-GHz flux density with WMAP K-band data for flat-spectrum sources suggests that WMAP flux densities may be low by a median factor of \simeq 1.2. The extrapolations of 5-GHz counts and luminosity functions of flat-spectrum radio quasars using the observed distribution of the 5-18.5 GHz spectral indices match those derived directly from WMAP data, indicating that the high-frequency WMAP survey does not detect any large population of FSRQs with anomalous spectra.