Using the NOvA neutrino detectors, a broad search has been performed for any signal coincident with 28 gravitational wave events detected by the LIGO/Virgo Collaboration between September 2015 and July 2019. For all of these events, NOvA is sensitive to possible arrival of neutrinos and cosmic rays of GeV and higher energies. For five (seven) events in the NOvA Far (Near) Detector, timely public alerts from the LIGO/Virgo Collaboration allowed recording of MeV-scale events. No signal candidates were found.

Using two years of data from the NOvA Near Detector at Fermilab, we report a seasonal variation of cosmic ray induced multiple-muon event rates which has an opposite phase to the seasonal variation in the atmospheric temperature. The strength of the seasonal multipl$ increase as a function of the muon multiplicity. However, no significant dependence of the strength of the seasonal variation of the multiple-muon variation is seen as a function of the muon zenith angle, or the spatial or angular separation between the correlated muons.

The {\gamma}-ray sky can be decomposed into individually detected sources, diffuse emission attributed to the interactions of Galactic cosmic rays with gas and radiation fields, and a residual all-sky emission component commonly called the isotropic diffuse {\gamma}-ray background (IGRB). The IGRB comprises all extragalactic emissions too faint or too diffuse to be resolved in a given survey, as well as any residual Galactic foregrounds that are approximately isotropic. The first IGRB measurement with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) used 10 months of sky-survey data and considered an energy range between 200 MeV and 100 GeV. Improvements in event selection and characterization of cosmic-ray backgrounds, better understanding of the diffuse Galactic emission, and a longer data accumulation of 50 months, allow for a refinement and extension of the IGRB measurement with the LAT, now covering the energy range from 100 MeV to 820 GeV. The IGRB spectrum shows a significant high-energy cutoff feature, and can be well described over nearly four decades in energy by a power law with exponential cutoff having a spectral index of $2.32\pm0.02$ and a break energy of $(279\pm52)$ GeV using our baseline diffuse Galactic emission model. The total intensity attributed to the IGRB is $(7.2\pm0.6) \times 10^{-6}$ cm$^{-2}$ s$^{-1}$ sr$^{-1}$ above 100 MeV, with an additional $+15$%/$-30$% systematic uncertainty due to the Galactic diffuse foregrounds.

Current theories predict relativistic hadronic particle populations in clusters of galaxies in addition to the already observed relativistic leptons. In these scenarios hadronic interactions give rise to neutral pions which decay into $\gamma$ rays, that are potentially observable with the Large Area Telescope (LAT) on board the Fermi space telescope. We present a joint likelihood analysis searching for spatially extended $\gamma$-ray emission at the locations of 50 galaxy clusters in 4 years of Fermi-LAT data under the assumption of the universal cosmic-ray model proposed by Pinzke & Pfrommer (2010). We find an excess at a significance of $2.7\sigma$ which upon closer inspection is however correlated to individual excess emission towards three galaxy clusters: Abell 400, Abell 1367 and Abell 3112. We discuss these cases in detail and conservatively attribute the emission to unmodeled background (for example, radio galaxies within the clusters). Through the combined analysis of 50 clusters we exclude hadronic injection efficiencies in simple hadronic models above 21% and establish limits on the cosmic-ray to thermal pressure ratio within the virial radius, $R_{200}$, to be below 1.2-1.4% depending on the morphological classification. In addition we derive new limits on the $\gamma$-ray flux from individual clusters in our sample.

Recent accurate measurements of cosmic-ray (CR) species by ATIC-2, CREAM, and PAMELA reveal an unexpected hardening in the proton and He spectra above a few hundred GeV, a gradual softening of the spectra just below a few hundred GeV, and a harder spectrum of He compared to that of protons. These newly-discovered features may offer a clue to the origin of high-energy CRs. We use the ${\it Fermi}$ Large Area Telescope observations of the $\gamma$-ray emission from the Earth's limb for an indirect measurement of the local spectrum of CR protons in the energy range $\sim 90~$GeV-$6~$TeV (derived from a photon energy range $15~$GeV-$1~$TeV). Our analysis shows that single power law and broken power law spectra fit the data equally well and yield a proton spectrum with index $2.68 \pm 0.04$ and $2.61 \pm 0.08$ above $\sim 200~$GeV, respectively.

Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial pulse up to 2.5 s is possibly the brightest well-isolated pulse observed to date. A fine time resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly correlated power-law behavior is observed between the luminosity and the spectral peak energy that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the existing models to account for all of the observed spectral and temporal behaviors simultaneously.

Cosmic rays are particles (mostly protons) accelerated to relativistic speeds. Despite wide agreement that supernova remnants (SNRs) are the sources of galactic cosmic rays, unequivocal evidence for the acceleration of protons in these objects is still lacking. When accelerated protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma rays. This offers a compelling way to detect the acceleration sites of protons. The identification of pion-decay gamma rays has been difficult because high-energy electrons also produce gamma rays via bremsstrahlung and inverse Compton scattering. We detected the characteristic pion-decay feature in the gamma-ray spectra of two SNRs, IC 443 and W44, with the Fermi Large Area Telescope. This detection provides direct evidence that cosmic-ray protons are accelerated in SNRs.

We have performed an analysis of the diffuse gamma-ray emission with the Fermi Large Area Telescope in the Milky Way Halo region searching for a signal from dark matter annihilation or decay. In the absence of a robust dark matter signal, constraints are presented. We consider both gamma rays produced directly in the dark matter annihilation/decay and produced by inverse Compton scattering of the e+e- produced in the annihilation/decay. Conservative limits are derived requiring that the dark matter signal does not exceed the observed diffuse gamma-ray emission. A second set of more stringent limits is derived based on modeling the foreground astrophysical diffuse emission using the GALPROP code. Uncertainties in the height of the diffusive cosmic-ray halo, the distribution of the cosmic-ray sources in the Galaxy, the index of the injection cosmic-ray electron spectrum and the column density of the interstellar gas are taken into account using a profile likelihood formalism, while the parameters governing the cosmic-ray propagation have been derived from fits to local cosmic-ray data. The resulting limits impact the range of particle masses over which dark matter thermal production in the early Universe is possible, and challenge the interpretation of the PAMELA/Fermi-LAT cosmic ray anomalies as annihilation of dark matter.

Millisecond pulsars, old neutron stars spun-up by accreting matter from a companion star, can reach high rotation rates of hundreds of revolutions per second. Until now, all such "recycled" rotation-powered pulsars have been detected by their spin-modulated radio emission. In a computing-intensive blind search of gamma-ray data from the Fermi Large Area Telescope (with partial constraints from optical data), we detected a 2.5-millisecond pulsar, PSR J1311-3430. This unambiguously explains a formerly unidentified gamma-ray source that had been a decade-long enigma, confirming previous conjectures. The pulsar is in a circular orbit with an orbital period of only 93 minutes, the shortest of any spin-powered pulsar binary ever found.

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

We report on the gamma-ray observations of giant molecular clouds Orion A and B with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. The gamma-ray emission in the energy band between \sim100 MeV and \sim100 GeV is predicted to trace the gas mass distribution in the clouds through nuclear interactions between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production cross-section for the nuclear interaction is known to \sim10% precision which makes the LAT a powerful tool to measure the gas mass column density distribution of molecular clouds for a known CR intensity. We present here such distributions for Orion A and B, and correlate them with those of the velocity integrated CO intensity (WCO) at a 1{\deg} \times1{\deg} pixel level. The correlation is found to be linear over a WCO range of ~10 fold when divided in 3 regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The Wco-to-mass conversion factor, Xco, is found to be \sim2.3\times10^20 cm-2(K km s-1)-1 for the high-longitude part of Orion A (l > 212{\deg}), \sim1.7 times higher than \sim1.3 \times 10^20 found for the rest of Orion A and B. We interpret the apparent high Xco in the high-longitude region of Orion A in the light of recent works proposing a non-linear relation between H2 and CO densities in the diffuse molecular gas. Wco decreases faster than the H2 column density in the region making the gas "darker" to Wco.

Recent detections of the starburst galaxies M82 and NGC 253 by gamma-ray telescopes suggest that galaxies rapidly forming massive stars are more luminous at gamma-ray energies compared to their quiescent relatives. Building upon those results, we examine a sample of 69 dwarf, spiral, and luminous and ultraluminous infrared galaxies at photon energies 0.1-100 GeV using 3 years of data collected by the Large Area Telescope (LAT) on the \textit{Fermi Gamma-ray Space Telescope} (\textit{Fermi}). Measured fluxes from significantly detected sources and flux upper limits for the remaining galaxies are used to explore the physics of cosmic rays in galaxies. We find further evidence for quasi-linear scaling relations between gamma-ray luminosity and both radio continuum luminosity and total infrared luminosity which apply both to quiescent galaxies of the Local Group and low-redshift starburst galaxies (conservative $P$-values $\lesssim0.05$ accounting for statistical and systematic uncertainties). The normalizations of these scaling relations correspond to luminosity ratios of $\log(L_{0.1-100 \rm{GeV}}/L_{1.4 \rm{GHz}}) = 1.7 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ and $\log(L_{0.1-100 \rm{GeV}}/L_{8-1000 \mu\rm{m}}) = -4.3 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ for a galaxy with a star formation rate of 1 $M_{\odot}$ yr$^{-1}$, assuming a Chabrier initial mass function. Using the relationship between infrared luminosity and gamma-ray luminosity, the collective intensity of unresolved star-forming galaxies at redshifts $0<z<2.5$ above 0.1 GeV is estimated to be 0.4-2.4 $\times 10^{-6}$ ph cm$^{-2}$ s$^{-1}$ sr$^{-1}$ (4-23% of the intensity of the isotropic diffuse component measured with the LAT). We anticipate that $\sim10$ galaxies could be detected by their cosmic-ray induced gamma-ray emission during a 10-year \textit{Fermi} mission.

Dark matter particle annihilation or decay can produce monochromatic gamma-ray lines and contribute to the diffuse gamma-ray background. Flux upper limits are presented for gamma-ray spectral lines from 7 to 200 GeV and for the diffuse gamma-ray background from 4.8 GeV to 264 GeV obtained from two years of Fermi Large Area Telescope data integrated over most of the sky. We give cross section upper limits and decay lifetime lower limits for dark matter models that produce gamma-ray lines or contribute to the diffuse spectrum, including models proposed as explanations of the PAMELA and Fermi cosmic-ray data.

We measured separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope. Because the instrument does not have an onboard magnet, we distinguish the two species by exploiting the Earth's shadow, which is offset in opposite directions for opposite charges due to the Earth's magnetic field. We estimate and subtract the cosmic-ray proton background using two different methods that produce consistent results. We report the electron-only spectrum, the positron-only spectrum, and the positron fraction between 20 GeV and 200 GeV. We confirm that the fraction rises with energy in the 20-100 GeV range. The three new spectral points between 100 and 200 GeV are consistent with a fraction that is continuing to rise with energy.

The contribution of unresolved sources to the diffuse gamma-ray background could induce anisotropies in this emission on small angular scales. We analyze the angular power spectrum of the diffuse emission measured by the Fermi LAT at Galactic latitudes |b| > 30 deg in four energy bins spanning 1 to 50 GeV. At multipoles \ell \ge 155, corresponding to angular scales \lesssim 2 deg, angular power above the photon noise level is detected at >99.99% CL in the 1-2 GeV, 2-5 GeV, and 5-10 GeV energy bins, and at >99% CL at 10-50 GeV. Within each energy bin the measured angular power takes approximately the same value at all multipoles \ell \ge 155, suggesting that it originates from the contribution of one or more unclustered source populations. The amplitude of the angular power normalized to the mean intensity in each energy bin is consistent with a constant value at all energies, C_P/<I>^2 = 9.05 +/- 0.84 x 10^{-6} sr, while the energy dependence of C_P is consistent with the anisotropy arising from one or more source populations with power-law photon spectra with spectral index \Gamma_s = 2.40 +/- 0.07. We discuss the implications of the measured angular power for gamma-ray source populations that may provide a contribution to the diffuse gamma-ray background.

Numerical simulations based on the Lambda-CDM model of cosmology predict a large number of as yet unobserved Galactic dark matter satellites. We report the results of a Large Area Telescope (LAT) search for these satellites via the gamma-ray emission expected from the annihilation of weakly interacting massive particle (WIMP) dark matter. Some dark matter satellites are expected to have hard gamma-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We sought to identify LAT sources with these characteristics, focusing on gamma-ray spectra consistent with WIMP annihilation through the $b \bar b$ channel. We found no viable dark matter satellite candidates using one year of data, and we present a framework for interpreting this result in the context of numerical simulations to constrain the velocity-averaged annihilation cross section for a conventional 100 GeV WIMP annihilating through the $b \bar b$ channel.

Satellite galaxies of the Milky Way are among the most promising targets for dark matter searches in gamma rays. We present a search for dark matter consisting of weakly interacting massive particles, applying a joint likelihood analysis to 10 satellite galaxies with 24 months of data of the Fermi Large Area Telescope. No dark matter signal is detected. Including the uncertainty in the dark matter distribution, robust upper limits are placed on dark matter annihilation cross sections. The 95% confidence level upper limits range from about 1e-26 cm^3 s^-1 at 5 GeV to about 5e-23 cm^3 s^-1 at 1 TeV, depending on the dark matter annihilation final state. For the first time, using gamma rays, we are able to rule out models with the most generic cross section (~3e-26 cm^3 s^-1 for a purely s-wave cross section), without assuming additional boost factors.

The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009.

We investigate the Fermi LAT gamma-ray and 15 GHz VLBA radio properties of a joint gamma-ray- and radio-selected sample of AGNs obtained during the first 11 months of the Fermi mission (2008 Aug 4 - 2009 Jul 5). Our sample contains the brightest 173 AGNs in these bands above declination -30 deg. during this period, and thus probes the full range of gamma-ray loudness (gamma-ray to radio band luminosity ratio) in the bright blazar population. The latter quantity spans at least four orders of magnitude, reflecting a wide range of spectral energy distribution (SED) parameters in the bright blazar population. The BL Lac objects, however, display a linear correlation of increasing gamma-ray loudness with synchrotron SED peak frequency, suggesting a universal SED shape for objects of this class. The synchrotron self-Compton model is favored for the gamma-ray emission in these BL Lacs over external seed photon models, since the latter predict a dependence of Compton dominance on Doppler factor that would destroy any observed synchrotron SED peak - gamma-ray loudness correlation. The high-synchrotron peaked (HSP) BL Lac objects are distinguished by lower than average radio core brightness temperatures, and none display large radio modulation indices or high linear core polarization levels. No equivalent trends are seen for the flat-spectrum radio quasars (FSRQ) in our sample. Given the association of such properties with relativistic beaming, we suggest that the HSP BL Lacs have generally lower Doppler factors than the lower-synchrotron peaked BL Lacs or FSRQs in our sample.

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.

The high-frequency-peaked BL Lacertae object RGB J0710+591 was observed in the very high-energy (VHE; E > 100 GeV) wave band by the VERITAS array of atmospheric Cherenkov telescopes. The observations, taken between 2008 December and 2009 March and totaling 22.1 hr, yield the discovery of VHE gamma rays from the source. RGB J0710+591 is detected at a statistical significance of 5.5 standard deviations (5.5{\sigma}) above the background, corresponding to an integral flux of (3.9 +/- 0.8) x 10-12 cm-2 s-1 (3% of the Crab Nebula's flux) above 300 GeV. The observed spectrum can be fit by a power law from 0.31 to 4.6 TeV with a photon spectral index of 2.69 +/- 0.26stat +/- 0.20sys. These data are complemented by contemporaneous multiwavelength data from the Fermi Large Area Telescope, the Swift X-ray Telescope, the Swift Ultra-Violet and Optical Telescope, and the Michigan-Dartmouth-MIT observatory. Modeling the broadband spectral energy distribution (SED) with an equilibrium synchrotron self-Compton model yields a good statistical fit to the data. The addition of an external-Compton component to the model does not improve the fit nor brings the system closer to equipartition. The combined Fermi and VERITAS data constrain the properties of the high-energy emission component of the source over 4 orders of magnitude and give measurements of the rising and falling sections of the SED.

The first published Fermi large area telescope (Fermi-LAT) measurement of the isotropic diffuse gamma-ray emission is in good agreement with a single power law, and is not showing any signature of a dominant contribution from dark matter sources in the energy range from 20 to 100 GeV. We use the absolute size and spectral shape of this measured flux to derive cross section limits on three types of generic dark matter candidates: annihilating into quarks, charged leptons and monochromatic photons. Predicted gamma-ray fluxes from annihilating dark matter are strongly affected by the underlying distribution of dark matter, and by using different available results of matter structure formation we assess these uncertainties. We also quantify how the dark matter constraints depend on the assumed conventional backgrounds and on the Universe's transparency to high-energy gamma-rays. In reasonable background and dark matter structure scenarios (but not in all scenarios we consider) it is possible to exclude models proposed to explain the excess of electrons and positrons measured by the Fermi-LAT and PAMELA experiments. Derived limits also start to probe cross sections expected from thermally produced relics (e.g. in minimal supersymmetry models) annihilating predominantly into quarks. For the monochromatic gamma-ray signature, the current measurement constrains only dark matter scenarios with very strong signals.

Dark matter (DM) particle annihilation or decay can produce monochromatic $\gamma$-rays readily distinguishable from astrophysical sources. $\gamma$-ray line limits from 30 GeV to 200 GeV obtained from 11 months of Fermi Large Area Space Telescope data from 20-300 GeV are presented using a selection based on requirements for a $\gamma$-ray line analysis, and integrated over most of the sky. We obtain $\gamma$-ray line flux upper limits in the range $0.6-4.5\times 10^{-9}\mathrm{cm}^{-2}\mathrm{s}^{-1}$, and give corresponding DM annihilation cross-section and decay lifetime limits. Theoretical implications are briefly discussed.

We report on the observations of 14 dwarf spheroidal galaxies with the Fermi Gamma-Ray Space Telescope taken during the first 11 months of survey mode operations. The Fermi telescope provides a new opportunity to test particle dark matter models through the expected gamma-ray emission produced by pair annihilation of weakly interacting massive particles (WIMPs). Local Group dwarf spheroidal galaxies, the largest galactic substructures predicted by the cold dark matter scenario, are attractive targets for such indirect searches for dark matter because they are nearby and among the most extreme dark matter dominated environments. No significant gamma-ray emission was detected above 100 MeV from the candidate dwarf galaxies. We determine upper limits to the gamma-ray flux assuming both power-law spectra and representative spectra from WIMP annihilation. The resulting integral flux above 100 MeV is constrained to be at a level below around 10^-9 photons cm^-2 s^-1. Using recent stellar kinematic data, the gamma-ray flux limits are combined with improved determinations of the dark matter density profile in 8 of the 14 candidate dwarfs to place limits on the pair annihilation cross-section of WIMPs in several widely studied extensions of the standard model. With the present data, we are able to rule out large parts of the parameter space where the thermal relic density is below the observed cosmological dark matter density and WIMPs (neutralinos here) are dominantly produced non-thermally, e.g. in models where supersymmetry breaking occurs via anomaly mediation. The gamma-ray limits presented here also constrain some WIMP models proposed to explain the Fermi and PAMELA e^+e^- data, including low-mass wino-like neutralinos and models with TeV masses pair-annihilating into muon-antimuon pairs. (Abridged)

(Abridged) We have conducted a detailed investigation of the broad-band spectral properties of the \gamma-ray selected blazars of the Fermi LAT Bright AGN Sample (LBAS). By combining our accurately estimated Fermi gamma-ray spectra with Swift, radio, infra-red, optical and other hard X-ray/gamma-ray data, collected within three months of the LBAS data taking period, we were able to assemble high-quality and quasi-simultaneous Spectral Energy Distributions (SED) for 48 LBAS blazars.The SED of these gamma-ray sources is similar to that of blazars discovered at other wavelengths, clearly showing, in the usual Log $\nu $ - Log $\nu$ F$_\nu$ representation, the typical broad-band spectral signatures normally attributed to a combination of low-energy synchrotron radiation followed by inverse Compton emission of one or more components. We have used these SEDs to characterize the peak intensity of both the low and the high-energy components. The results have been used to derive empirical relationships that estimate the position of the two peaks from the broad-band colors (i.e. the radio to optical and optical to X-ray spectral slopes) and from the gamma-ray spectral index. Our data show that the synchrotron peak frequency $\nu_p^S$ is positioned between 10$^{12.5}$ and 10$^{14.5}$ Hz in broad-lined FSRQs and between $10^{13}$ and $10^{17}$ Hz in featureless BL Lacertae objects.We find that the gamma-ray spectral slope is strongly correlated with the synchrotron peak energy and with the X-ray spectral index, as expected at first order in synchrotron - inverse Compton scenarios. However, simple homogeneous, one-zone, Synchrotron Self Compton (SSC) models cannot explain most of our SEDs, especially in the case of FSRQs and low energy peaked (LBL) BL Lacs. (...)