Gamma-ray Luminosity Research Articles

Radio pulsars with expected gamma radiation and gamma-ray pulsars as pulsating radio emitters

Pulsars play a crucial astrophysical role as highly energetic compact radio, X-ray and gamma-ray sources. Our previous works show that radio pulsars identified as pulsing gamma-ray sources by the Large Area Telescope (LAT) on board the Fermi Gamma-Ray Space Telescope have high values of magnetic field near the light cylinder, two-three orders of magnitude stronger compared with the magnetic fields of radio pulsars: log Blc (G) are 3.60–3.95 and 1.75 correspondingly. Moreover, their losses of rotational energy are also three orders higher than the corresponding values for the main group of radio pulsars on average: log Ė (erg s−1) = 35.37–35.53 and 32.64. The correlation between gamma-ray luminosities and radio luminosities is found. It allows us to select those objects from all sets of known radio pulsars that can be detected as gamma-ray pulsars with high probability. We provide a list of such radio pulsars and propose to search for gamma emission from these objects. On the other hand, the known catalog of gamma-ray pulsars contains some sources which are not currently identified as radio pulsars. Some of them have large values of gamma-ray luminosities and according to the obtained correlation, we can expect marked radio emission from these objects. We give the list of such pulsars and expected flux densities to search for radiation at frequencies 1400 and 111 MHz.

Millisecond Pulsar Origin of the Galactic Center Excess and Extended Gamma-Ray Emission from Andromeda: A Closer Look

Abstract A new measurement of a spatially extended gamma-ray signal from the center of the Andromeda galaxy (M31) has recently been published by the Fermi-LAT collaboration, reporting that the emission broadly resembles the so-called Galactic center excess (GCE) of the Milky Way (MW). The weight of the evidence is steadily accumulating on a millisecond pulsar (MSPs) origin for the GCE. These elements prompt us to compare these observations with what is, perhaps, the simplest model for an MSP population, which is solely obtained by rescaling of the MSP luminosity function that is determined in the local MW disk via the respective stellar mass of the systems. Remarkably, we find that without free fitting parameters, this model can account for both the energetics and the morphology of the GCE within uncertainties. For M31, the estimated luminosity due to primordial MSPs is expected to only contribute about a quarter of the detected emission, although a stronger contribution cannot be excluded given the large uncertainties. If correct, the model predicts that the M31 disk emission due to MSPs is not far below the present upper bound. We also discuss additional refinements of this simple model. Using the correlation between globular cluster gamma-ray luminosity and stellar encounter rate, we gauge the dynamical MSP formation in the bulge. This component is expected to contribute to the GCE only at a level of ≲5%, it could affect the signal’s morphology. We also comment on the limitations of our model and on future perspectives for improved diagnostics.

The broad-band properties of the intermediate synchrotron peaked BL Lac S2 0109+22 from radio to VHE gamma-rays

The MAGIC telescopes observed S2 0109+22 in 2015 July during its flaring activity in high energy gamma rays observed by Fermi-LAT. We analyse the MAGIC data to characterise the very high energy (VHE) gamma-ray emission of S2 0109+22, which belongs to the subclass of intermediate synchrotron peak (ISP) BL Lac objects. We study the multi-frequency emission in order to investigate the source classification. Finally, we compare the source long-term behaviour to other VHE gamma-ray emitting (TeV) blazars. We performed a temporal and spectral analysis of the data centred around the MAGIC interval of observation (MJD 57225-57231). Long-term radio and optical data have also been investigated using the discrete correlation function. The redshift of the source is estimated through optical host-galaxy imaging and also using the amount of VHE gamma-ray absorption. The quasi-simultaneous multi-frequency spectral energy distribution (SED) is modelled with the conventional one-zone synchrotron self-Compton (SSC) model. MAGIC observations resulted in the detection of the source at a significance level of $5.3\,\sigma$. The VHE gamma-ray emission of S2 0109+22 is variable on a daily time scale. VHE gamma-ray luminosity of the source is lower than the average of TeV BL Lacs. The optical polarization, and long-term optical/radio behaviour of the source are different from the general population of TeV blazars. All these findings agree with the classification of the source as an ISP BL Lac object. We estimate the source redshift as $z = 0.36 \pm 0.07$. The SSC parameters describing the SED are rather typical for blazars.

On the existence of a luminosity threshold of GRB jets in massive stars

ABSTRACT Motivated by the many associations of gamma-ray bursts (GRBs) with energetic supernova (SN) explosions, we study the propagation of relativistic jets within the progenitor star in which an SN shock wave may be launched briefly before the jets start to propagate. Based on analytic considerations and verified with an extensive set of 2D axisymmetric relativistic hydrodynamic simulations, we have estimated a threshold intrinsic jet luminosity, $L_{\rm j}^{\rm thr}$, for successfully launching a jet. This threshold depends on the structure of the progenitor and, thus, it is sensitive to its mass and to its metallicity. For a prototype host of cosmological long GRBs, a low-metallicity star of 35 M⊙, it is $L_{\rm j}^{\rm thr}\simeq 1.35\times 10^{49}$ erg s−1. The observed equivalent isotropic gamma-ray luminosity, $L_{\rm \gamma ,iso,BO} \simeq 4 \epsilon _\gamma L_{\rm j} \theta _{\rm BO}^{-2}$, crucially depends on the jet opening angle after breakout, θBO, and on the efficiency for converting the intrinsic jet luminosity into γ-radiation, εγ. Highly energetic jets can produce low-luminosity events if either their opening angle after the breakout is large, which is found in our models, or if the conversion efficiency of kinetic and internal energy into radiation is low enough. Beyond this theoretical analysis, we show how the presence of an SN shock wave may reduce this luminosity threshold by means of numerical simulations. We foresee that the high-energy transients released by jets produced near the luminosity threshold will be more similar to llGRBs or X-ray flashes than to GRBs.

Evolution of jets driven by relativistic radiation hydrodynamics as long and low-luminosity GRBs

We present numerical simulations of jets modelled with relativistic radiation hydrodynamics (RRH), which evolve across two environments: (i) a stratified surrounding medium and (ii) a 16TI progenitor model. We consider opacities consistent with various processes of interaction between the fluid and radiation, specifically free–free, bound–free, bound–bound, and electron scattering. We explore various initial conditions, with different radiation energy densities of the beam in hydrodynamical and radiation-pressure-dominated scenarios, considering only highly relativistic jets. In order to investigate the impact of the radiation field on the evolution of the jets, we compare our results with purely hydrodynamical jets. Comparing among jets driven by an RRH, we find that radiation-pressure-dominated jets propagate slightly faster than gas pressure dominated ones. Finally, we construct the luminosity light curves (LCs) associated with the two cases. The construction of LCs uses the fluxes of the radiation field that is fully coupled to the hydrodynamics equations during the evolution. The main properties of the jets propagating on the stratified surrounding medium are that the LCs show the same order of magnitude as the gamma-ray luminosity of typical Long gamma-ray bursts 1050–1054 erg s−1, and the difference between the radiation and gas temperatures is of nearly one order of magnitude. The properties of jets breaking out from the progenitor star model are that the LCs are of the order of magnitude of low-luminosity GRBs 1046−1049 erg s−1, and in this scenario, the difference between the gas and radiation temperature is of four orders of magnitude, which is a case far from thermal equilibrium.

Multiwavelength photometric and spectropolarimetric analysis of the FSRQ 3C 279

In this paper, we present light curves for 3C 279 over a time period of six years; from 2008 to 2014. Our multiwavelength data comprise 1 mm to gamma-rays, with additional optical polarimetry. Based on the behaviour of the gamma-ray light curve with respect to other bands, we identified three different activity periods. One of the activity periods shows anomalous behaviour with no gamma-ray counterpart associated with optical and NIR flares. Another anomalous activity period shows a flare in gamma-rays, 1 mm and polarization degree, however, it does not have counterparts in the UV continuum, optical and NIR bands. We find a significant overall correlation of the UV continuum emission, the optical and NIR bands. This correlation suggests that the NIR to UV continuum is co-spatial. We also find a correlation between the UV continuum and the 1 mm data, which implies that the dominant process in producing the UV continuum is synchrotron emission. The gamma-ray spectral index shows statistically significant variability and an anti-correlation with the gamma-ray luminosity. We demonstrate that the dominant gamma-ray emission mechanism in 3C 279 changes over time. Alternatively, the location of the gamma-ray emission zone itself may change depending on the activity state of the central engine.

An MAD explanation for the correlation between bulk Lorentz factor and minimum variability time-scale

We offer an explanation for the anti-correlation between the minimum variability timescale ($MTS$) in the prompt emission light curve of gamma-ray bursts (GRBs) and the estimated bulk Lorentz factor of these GRBs, in the context of a magnetically arrested disk (MAD) model. In particular, we show that previously derived limits on the maximum available energy per baryon in a Blandford-Znajek jet leads to a relationship between the characteristic MAD timescale, $t_{MAD}$, in GRBs and the maximum bulk Lorentz factor: $t_{MAD} \propto \Gamma^{-6}$, somewhat steeper than (although within the error bars of) the fitted relationship found in the GRB data. Similarly, the MAD model also naturally accounts for the observed anti-correlation between $MTS$ and gamma-ray luminosity $L$ in the GRB data, and we estimate the accretion rates of the GRB disks (given these luminosities) in the context of this model. Both of these correlations ($MTS-\Gamma$ and $MTS-L$) are also observed in the AGN data, and we discuss the implications of our results in the context of both GRB and blazar systems.

High-energy gamma-ray and neutrino production in star-forming galaxies across cosmic time: Difficulties in explaining the IceCube data

Abstract We present new theoretical modeling to predict the luminosity and spectrum of gamma-ray and neutrino emission of a star-forming galaxy, from the star formation rate (ψ), gas mass (Mgas), stellar mass, and disk size, taking into account production, propagation, and interactions of cosmic rays. The model reproduces the observed gamma-ray luminosities of nearby galaxies detected by Fermi better than the simple power-law models as a function of ψ or ψMgas. This model is then used to predict the cosmic background flux of gamma-rays and neutrinos from star-forming galaxies, by using a semi-analytical model of cosmological galaxy formation that reproduces many observed quantities of local and high-redshift galaxies. Calibration of the model using gamma-ray luminosities of nearby galaxies allows us to make a more reliable prediction than previous studies. In our baseline model, star-forming galaxies produce about 20% of the isotropic gamma-ray background unresolved by Fermi, and only 0.5% of IceCube neutrinos. Even with an extreme model assuming a hard injection cosmic-ray spectral index of 2.0 for all galaxies, at most 22% of IceCube neutrinos can be accounted for. These results indicate that it is difficult to explain most of the IceCube neutrinos by star-forming galaxies, without violating the gamma-ray constraints from nearby galaxies.

The population of TeV pulsar wind nebulae in the H.E.S.S. Galactic Plane Survey

The nine-year H.E.S.S. Galactic Plane Survey (HGPS) has yielded the most uniform observation scan of the inner Milky Way in the TeV gamma-ray band to date. The sky maps and source catalogue of the HGPS allow for a systematic study of the population of TeV pulsar wind nebulae found throughout the last decade. To investigate the nature and evolution of pulsar wind nebulae, for the first time we also present several upper limits for regions around pulsars without a detected TeV wind nebula. Our data exhibit a correlation of TeV surface brightness with pulsar spin-down power Ė. This seems to be caused both by an increase of extension with decreasing Ė, and hence with time, compatible with a power law RPWN(Ė) ~Ė−0.65±0.20, and by a mild decrease of TeV gamma-ray luminosity with decreasing Ė, compatible with L1−10 TeV ~Ė0.59±0.21. We also find that the offsets of pulsars with respect to the wind nebula centre with ages around 10 kyr are frequently larger than can be plausibly explained by pulsar proper motion and could be due to an asymmetric environment. In the present data, it seems that a large pulsar offset is correlated with a high apparent TeV efficiency L1−10 TeV∕Ė. In addition to 14 HGPS sources considered firmly identified pulsar wind nebulae and 5 additional pulsar wind nebulae taken from literature, we find 10 HGPS sources that are likely TeV pulsar wind nebula candidates. Using a model that subsumes the present common understanding of the very high-energy radiative evolution of pulsar wind nebulae, we find that the trends and variations of the TeV observables and limits can be reproduced to a good level, drawing a consistent picture of present-day TeV data and theory.

Polarimetric Monitoring of Jets with Kanata Telescope

The polarization of relativistic jets is of interest for the understanding of their origin, confinement, and propagation. However, even though numerous measurements have been performed, the mechanisms behind jet variability, creation, and composition are still debated. We have performed simultaneous gamma-ray and optical photopolarimetry observations of 45 blazars with the Kanata telescope since July 2008 to investigate the mechanisms of variability and search for a basic relation between the several subclasses of relativistic jets. Consequently, we found that a correlation between the gamma-ray and optical flux might be related to gamma-ray luminosity, and the maximum polarization degree might be related to gamma-ray luminosity or the ratio of gamma-ray to optical flux. These results imply that low gamma-ray luminosity blazars emit from multiple regions.

A Detailed Observational Analysis of V1324 Sco, the Most Gamma-Ray-luminous Classical Nova to Date

Abstract It has recently been discovered that some, if not all, classical novae emit GeV gamma-rays during outburst, but the mechanisms involved in the production ofgamma-rays are still not well understood. We present here a comprehensive multiwavelength data set—from radio to X-rays—for the most gamma-ray-luminous classical nova to date, V1324 Sco. Using this data set, we show that V1324 Sco is a canonical dusty Fe ii-type nova, with a maximum ejecta velocity of 2600 km s−1 and an ejecta mass of a few . There is also evidence for complex shock interactions, including a double-peaked radio light curve which shows high brightness temperatures at early times. To explore why V1324 Sco was so gamma-ray luminous, we present a model of the nova ejecta featuring strong internal shocks and find that higher gamma-ray luminosities result from higher ejecta velocities and/or mass-loss rates. Comparison of V1324 Sco with other gamma-ray-detected novae does not show clear signatures of either, and we conclude that a larger sample of similarly well-observed novae is needed to understand the origin and variation of gamma-rays in novae.

Discovery of GeV gamma-ray emission from the LMC B0443-6657 with the Fermi Large Area Telescope

We report the discovery of gamma-ray detection from the Large Magellanic Cloud (LMC) B0443-6657 using the Large Area Telescope (LAT) on board the \textit{Fermi Gamma-ray Space Telescope}. LMC B0443-6657 is a flat spectrum radio source, possibly associated with a supernova remnant in the Large Magellanic Cloud (LMC N4). Employing the LAT data of 8 years, our results show a significant excess ($>9.4\sigma$) of gamma-ray in the range of 0.2--100\,GeV above the gamma-ray background. A power-law function is found to be adequate to describe the $0.2-100$ GeV $\gamma$-ray spectrum, which yields a photon flux of $3.27\pm0.53\ \mathrm{photon\ \,cm}^2\ \mathrm{s}^{-1}$ with a photon index of $2.35\pm0.11$, corresponding to an isotropic gamma-ray luminosity of $5.3\times10^{40}\ \mathrm{erg\ \,s}^{-1}$. The hadronic model predicts a low X-ray and TeV flux while the leptonic model predicts an observable flux in these two energy bands. The follow-up observations of the LMC B0443-6657 in X-ray or TeV band would distinguish the radiation models of gamma-rays from this region.

High-energy Emission from Nonrelativistic Radiative Shocks: Application to Gamma-Ray Novae

Abstract The observation of GeV gamma-rays from novae by Fermi/LAT demonstrates that the nonrelativistic radiative shocks in these systems can accelerate particles to energies of at least ∼10 GeV. The low-energy extension of the same nonthermal particle distribution inevitably gives rise to emission in the hard X-ray band. Above , this radiation can escape the system without significant absorption/attenuation, and can potentially be detected by NuSTAR. We present theoretical models for hard X-ray and gamma-ray emission from radiative shocks in both leptonic and hadronic scenarios, accounting for the rapid evolution of the downstream properties due to the fast cooling of thermal plasma. We find that due to strong Coulomb losses, only a fraction of of the gamma-ray luminosity is radiated in the NuSTAR band; nevertheless, this emission could be detectable simultaneously with the LAT emission in bright gamma-ray novae with a ∼50 ks exposure. The spectral slope in hard X-rays is for typical nova parameters, thus serving as a testable prediction of the model. Our work demonstrates how combined hard X-ray and gamma-ray observations can be used to constrain properties of the nova outflow (velocity, density, and mass outflow rate) and particle acceleration at the shock. A very low X-ray to gamma-ray luminosity ratio ( ) would disfavor leptonic models for the gamma-ray emission. Our model can also be applied to other astrophysical environments with radiative shocks, including SNe IIn and colliding winds in massive star binaries.

Probing Intrinsic Properties of Short Gamma-Ray Bursts with Gravitational Waves.

Progenitors of short gamma-ray bursts are thought to be neutron stars coalescing with their companion black hole or neutron star, which are one of the main gravitational wave sources. We have devised a Bayesian framework for combining gamma-ray burst and gravitational wave information that allows us to probe short gamma-ray burst luminosities. We show that combined short gamma-ray burst and gravitational wave observations not only improve progenitor distance and inclination angle estimates, they also allow the isotropic luminosities of short gamma-ray bursts to be determined without the need for host galaxy or light-curve information. We characterize our approach by simulating 1000 joint short gamma-ray burst and gravitational wave detections by Advanced LIGO and Advanced Virgo. We show that ∼90% of the simulations have uncertainties on short gamma-ray burst isotropic luminosity estimates that are within a factor of two of the ideal scenario, where the distance is known exactly. Therefore, isotropic luminosities can be confidently determined for short gamma-ray bursts observed jointly with gravitational waves detected by Advanced LIGO and Advanced Virgo. Planned enhancements to Advanced LIGO will extend its range and likely produce several joint detections of short gamma-ray bursts and gravitational waves. Third-generation gravitational wave detectors will allow for isotropic luminosity estimates for the majority of the short gamma-ray burst population within a redshift of z∼1.

A New Look at Distances and Velocities of Neutron Stars

We take a fresh look at the determination of distances and velocities of neutron stars. The conversion of a parallax measurement into a distance, or distance probability distribution, has led to a debate quite similar to the one involving Cepheids, centering on the question whether priors can be used when discussing a single system. With the example of PSRJ0218+4232 we show that a prior is necessary to determine the probability distribution for the distance. The distance of this pulsar implies a gamma-ray luminosity larger than 10% of its spindown luminosity. For velocities the debate is whether a single Maxwellian describes the distribution for young pulsars. By limiting our discussion to accurate (VLBI) measurements we argue that a description with two Maxwellians, with distribution parameters sigma1=77 and sigma2=320 km/s, is significantly better. Corrections for galactic rotation, to derive velocities with respect to the local standards of rest, are insignificant.

Gamma rays from clumpy wind-jet interactions in high-mass microquasars

The stellar winds of the massive stars in high-mass microquasars are thought to be inhomogeneous. The interaction of these inhomogeneities, or clumps, with the jets of these objects may be a major factor in gamma-ray production. Our goal is to characterize a typical scenario of clump-jet interaction, and calculate the contribution of these interactions to the gamma-ray emission from these systems. We use axisymmetric, RHD simulations to model the emitting flow in a typical clump-jet interaction. Using the simulation results we perform a numerical calculation of the high-energy emission from one of these interactions. The radiative calculations are done for relativistic electrons locally accelerated at the jet shock, and the synchrotron and IC radiation spectra are computed for different stages of the shocked clump evolution. We also explore different parameter values, such as viewing angle and magnetic field strength. We generalize phenomenologically the results derived from one clump-jet interaction to multiple interactions under different wind models, estimating the clump-jet interaction rates, and the resulting luminosities in the GeV range. If particles are efficiently accelerated in clump-jet interactions, the apparent gamma-ray luminosity through IC scattering with the stellar photons can be significant even for rather strong magnetic fields and thus efficient synchrotron cooling. Moreover, despite the standing nature or slow motion of the jet shocks for most of the interaction stage, Doppler boosting in the postshock flow is relevant even for mildly relativistic jets. For clump-to-average wind density contrasts $\gtrsim 10$, clump-jet interactions could be bright enough to match the observed GeV luminosity in Cyg~X-1 and Cyg~X-3 when a jet is present in these sources, with required non-thermal-to-total available power fractions $\sim 0.01$ and $0.1$, respectively.

Enhanced gamma radiation towards the rotation axis from the immediate vicinity of extremely rotating black holes

Abstract We investigate the acceleration of electrons and positrons by magnetic-field-aligned electric fields in the polar funnel of an accreting black hole (BH). Applying the pulsar outer-gap theory to BH magnetospheres, we find that such a lepton accelerator arises in the immediate vicinity of the event horizon due to frame-dragging, and that their gamma-ray luminosity increases with decreasing accretion rate. Furthermore, we demonstrate that the gamma-ray flux is enhanced along the rotation axis by more than an order of magnitude if the BH spin increases from a = 0.90M to a = 0.9999M. As a result, if a ten-solar-mass, almost-maximally rotating BH is located within 3 kpc, when its accretion rate is between 0.005 and 0.01 per cent of the Eddington rate, its high-energy flare becomes detectable with the Fermi/Large Area Telescope, provided that the flare lasts longer than 1.2 months and that we view the source nearly along the rotation axis. In addition, its very high energy flux is marginally detectable with the Cherenkov Telescope Array, provided that the flare lasts longer than a night and that our viewing angle is about 45 deg with respect to the rotation axis.

Cherenkov telescope array extragalactic survey discovery potential and the impact of axion-like particles and secondary gamma rays

The Cherenkov Telescope Array (CTA) is about to enter construction phase and one of its main key science projects is to perform an unbiased survey in search of extragalactic sources. We make use of both the latest blazar gamma–ray luminosity function and spectral energy distribution to derive the expected number of detectable sources for both the planned Northern and Southern arrays of the CTA observatory. We find that a shallow, wide survey of about 0.5 hour per field of view would lead to the highest number of blazar detections. Furthermore, we investigate the effect of axion-like particles and secondary gamma rays from propagating cosmic rays on the source count distribution, since these processes predict different spectral shape from standard extragalactic background light attenuation. We can generally expect more distant objects in the secondary gamma-ray scenario, while axion-like particles do not significantly alter the expected distribution. Yet, we find that, these results strongly depend on the assumed magnetic field strength during the propagation. We also provide source count predictions for the High Altitude Water Cherenkov observatory (HAWC), the Large High Altitude Air Shower Observatory (LHAASO) and a novel proposal of a hybrid detector.

MOJAVE – XIV. Shapes and opening angles of AGN jets

We present 15 GHz stacked VLBA images of 373 jets associated with active galactic nuclei (AGN) having at least five observing epochs within a 20 yr time interval 1994-2015 from the MOJAVE programme and/or its precursor, the 2 cm VLBA Survey. These data are supplemented by 1.4 GHz single-epoch VLBA observations of 135 MOJAVE AGNs to probe larger scale jet structures. The typical jet geometry is found to be close to conical on scales from hundreds to thousands of parsecs, while a number of galaxies show quasi-parabolic streamlines on smaller scales. A true jet geometry in a considerable fraction of AGNs appears only after stacking epochs over several years. The jets with significant radial accelerated motion undergo more active collimation. We have analysed total intensity jet profiles transverse to the local jet ridgeline and derived both apparent and intrinsic opening angles of the flows, with medians of $21.5\deg$ and $1.3\deg$, respectively. The Fermi LAT-detected gamma-ray AGNs in our sample have, on average, wider apparent and narrower intrinsic opening angle, and smaller viewing angle than non LAT-detected AGNs. We have established a highly significant correlation between the apparent opening angle and gamma-ray luminosity, driven by Doppler beaming and projection effects.

The Fermi blazar sequence

We revisit the blazar sequence exploiting the complete, flux limited sample of blazars with known redshift detected by the Fermi satellite after 4 years of operations (the 3LAC sample). We divide the sources into gamma-ray luminosity bins, collect all the archival data for all blazars, and construct their spectral energy distribution (SED). We describe the average SED of blazars in the same luminosity bin through a very simple, completely phenomenological function consisting of two broken power laws connecting with a power law of fixed slope describing the radio emission. We do that separately for BL Lacs and for flat spectrum radio quasars (FSRQs) and also for all blazars together. The main results are: i) FSRQs display approximately the same SED as the luminosity increases, except for the fact that the relative importance of the high energy peak increases; ii) as a consequence, X-ray spectra of FSRQs become harder for larger luminosities; iii) BL Lacs form indeed a sequence: they become redder (i.e. the peak frequencies becomes smaller) for increasing luminosities, with a steeper gamma-ray slope and a larger dominance of the high energy peak; iv) for all blazars (BL Lacs+FSRQs) these properties becomes more prominent, as the highest luminosity bin is populated mostly by FSRQs and the lowest luminosity bin mostly by BL Lacs. This agrees with the original blazar sequence, although BL Lacs never have an average gamma-ray slope as hard as found in the original sequence. v) At high luminosities, a large fraction of FSRQs shows signs of thermal emission from the accretion disc, contributing in the optical-UV.