We performed a long-term optical (B, V, R bands), infra-red (J and K bands) and radio band (15, 22, 37 GHz band) study on the flat spectrum radio quasar, 3C 454.3, using the data collected over a period of more than 8 years (MJD 54500--57500). The temporal variability, spectral properties and inter-waveband correlations were studied by dividing the available data into smaller segments with more regular sampling. This helped us constrain the size and the relative locations of the emission regions for different wavebands. Spectral analysis of the source revealed the interplay between the accretion disk and jet emission. The source predominantly showed a redder-when-brighter trend, though we observed a bluer-when-brighter trend at high flux levels which could be signatures of particle acceleration and radiative cooling. Significant correlations with near-zero lag were seen between various optical/infra-red bands, indicating that these emission regions are co-spatial. Correlations with a time lag of about 10--100 days are seen between optical/infra-red and radio bands indicating these emissions arise from different regions. We also observe the DCF peak lag change from year to year. We try to explain these differences using a curved jet model where the different emission regions have different viewing angles resulting in a frequency dependent Doppler factor. This variable Doppler factor model explains the variability timescales and the variation in DCF peak lag between the radio and optical emissions in different segments. Lags of 6-180 days are seen between emissions in various radio bands, indicating a core-shift effect

We present a study of a sample of 223 radio-loud quasars (up to redshift $<$0.3) in order to investigate their spectral properties. Twenty-six of these radio-loud quasars are identified as Flat Spectrum Radio Quasars (FSRQs) and fifty-four are identified as Steep Spectrum Radio Quasars (SSRQs) based on their radio spectral index. We study the [O III] line properties of these quasars to investigate the origin and properties of blue wings (shift of the profile towards lower wavelengths) and blue outliers (shift of the whole spectroscopic feature). Most of the quasars show blue wings with velocities up to 420 km $s^{-1}$. We find that around 17\% of the quasars show outliers with velocities spanning 419 to -315 km $s^{-1}$. Finally, we revisit the $\it M_{\rm BH} - \sigma$ relation of our sample using [S II]$\lambda$6716, 6731 and [O III] linewidths as surrogates for stellar velocity dispersions, $\sigma$, to investigate their location on the $\it M_{\rm BH} - \sigma$ relation for quiescent galaxies. Due to the strong blending of [S II] with $\rm H_{\alpha}$, we could estimate $\sigma_{[\rm SII]}$ of only 123 quasars. We find that the radio-loud quasars do not show a relationship between $\it M_{\rm BH}$ and $\sigma_{\rm [SII]/[OIII]}$ up to a redshift of 0.3, although they cluster around the local relation. We find an overall offset of 0.12$\pm$0.05 dex of our sample of radio-loud quasars from the $\it M_{\rm BH} - \sigma$ relation of quiescent galaxies. Quasars in our highest redshift bin (z=0.25-0.3) show a deviation of $\sim$0.33 $\pm$ 0.06 dex with respect to the local relation.

We present X-ray flux and spectral analyses of the three pointed Suzaku observations of the TeV high synchrotron peak blazar Mrk 421 taken throughout its complete operational duration. The observation taken on 5 May 2008 is, at 364.6 kiloseconds (i.e., 101.3 hours), the longest and most evenly sampled continuous observation of this source, or any blazar, in the X-ray energy 0.8 - 60 keV until now. We found large amplitude intra-day variability in all soft and hard bands in all the light curves. The discrete correction function analysis of the light curves in soft and hard bands peaks on zero lag, showing that the emission in hard and soft bands are cospatial and emitted from the same population of leptons. The hardness ratio plots imply that the source is more variable in the harder bands compared to the softer bands. The source is harder-when-brighter, following the general behavior of high synchrotron peak blazars. Power spectral densities of all three light curves are red noise dominated, with a range of power spectra slopes. If one assumes that the emission originates very close to the central super massive black hole, a crude estimate for its mass, of ~ 4 * 10^{8} M_{\odot}, can be made; but if the variability is due to perturbations arising there that are advected into the jet and are thus Doppler boosted, substantially higher masses are consistent with the quickest seen variations. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux and spectral variability.

We present the results of photometric observations of three TeV blazars, 3C 66A, S5 0954+658 and BL Lacertae, during the period 2013--2017. Our extensive observations were performed in a total of 360 nights which produced $\sim$6820 image frames in BVRI bands. We study flux and spectral variability of these blazars on these lengthy timescales. We also examine the optical Spectral Energy Distributions of these blazars, which are crucial in understanding the emission mechanism of long-term variability in blazars. All three TeV blazars exhibited strong flux variability during our observations. The colour variations are mildly chromatic on long timescales for two of them. The nature of the long-term variability of 3C 66A and S5 0954+658 is consistent with a model of a non-thermal variable component that has a continuous injection of relativistic electrons with power law distributions around 4.3 and 4.6, respectively. However, the long-term flux and colour variability of BL Lac suggests that these can arise from modest changes in velocities or viewing angle toward the emission region, leading to variations in the Doppler boosting of the radiation by a factor ~1.2 over the period of these observations.

We report on a recent multi-band optical photometric and polarimetric observational campaign of the blazar OJ 287 which was carried out during September 2016 -- December 2017. We employed nine telescopes in Bulgaria, China, Georgia, Japan, Serbia, Spain and the United States. We collected over 1800 photometric image frames in BVRI bands and over 100 polarimetric measurements over ~175 nights. In 11 nights with many quasi-simultaneous multi-band (V, R, I) observations, we did not detect any genuine intraday variability in flux or color. On longer timescales, multiple flaring events were seen. Large changes in color with respect to time and in a color--magnitude diagram were seen, and while only a weak systematic variability trend was noticed in color with respect to time, the color--magnitude diagram shows a bluer-when-brighter trend. Large changes in the degree of polarization, and substantial swings in the polarization angle were detected. The fractional Stokes parameters of the polarization showed a systematic trend with time in the beginning of these observations, followed by chaotic changes and then an apparently systematic variation at the end. These polarization changes coincide with the detection and duration of the source at very high energies as seen by VERITAS. The spectral index shows a systematic variation with time and V-band magnitude. We briefly discuss possible physical mechanisms that could explain the observed flux, color, polarization, and spectral variability.

We present a multi-wavelength spectral and temporal investigation of OJ 287 emission during its strong optical-to-X-ray activity between July 2016 - July 2017. The daily $\gamma$-ray fluxes from \emph{Fermi}-LAT are consistent with no variability. The strong optical-to-X-ray variability is accompanied by a change in power-law spectral index of the X-ray spectrum from $< 2$ to $>2$, with variations often associated with changes in optical polarization properties. Cross-correlations between optical-to-X-ray emission during four continuous segments show simultaneous optical-ultraviolet (UV) variations while the X-ray and UV/optical are simultaneous only during the middle two segments. In the first segment, the results suggest X-rays lag the optical/UV, while in the last segment X-rays lead by $\sim$ 5-6 days. The last segment also shows a systematic trend with variations appearing first at higher energies followed by lower energy ones. The LAT spectrum before the VHE activity is similar to preceding quiescent state spectrum while it hardens during VHE activity period and is consistent with the extrapolated VHE spectrum during the latter. Overall, the broadband spectral energy distributions (SEDs) during high activity periods are a combination of a typical OJ 287 SED and an HBL SED, and can be explained in a two-zone leptonic model, with the second zone located at parsec scales, beyond the broad line region, being responsible for the HBL-like spectrum. The change of polarization properties from systematic to chaotic and back to systematic, before, during and after the VHE activity, suggest dynamic roles for magnetic fields and turbulence.

We monitored BL Lacertae frequently during 2014 - 2016 when it was generally in a high state. We searched for intra-day variability for 43 nights using quasi-simultaneous measurements in the B, V, R, and I bands (totaling 143 light curves); the typical sampling interval was about eight minutes. On hour-like timescales, BL Lac exhibited significant variations during 13 nights in various optical bands. Significant spectral variations are seen during most of these nights such that the optical spectrum becomes bluer when brighter. The amplitude of variability is usually greater for longer observations but is lower when BL Lac is brighter. No evidence for periodicities or characteristic variability time-scales in the light curves was found. The color variations are mildly chromatic on long timescales.

The synchrotron peak of the X-ray bright High Energy Peaked Blazar (HBL) PKS 2155$-$304 occurs in the UV-EUV region and hence its X-ray emission (0.6--10 keV) lies mostly in the falling part of the synchrotron hump. We aim to study the X-ray emission of PKS 2155$-$304 during different intensity states in 2009$-$2014 using XMM$-$Newton satellite. We studied the spectral curvature of all of the observations to provide crucial information on the energy distribution of the non-thermal particles. Most of the observations show curvature or deviation from a single power-law and can be well modeled by a log parabola model. In some of the observations, we find spectral flattening after 6 keV. In order to find the possible origin of the X-ray excess, we built the Multi-band Spectral Energy distribution (SED). We find that the X-ray excess in PKS 2155--304 is difficult to fit in the one zone model but, could be easily reconciled in the spine/layer jet structure. The hard X-ray excess can be explained by the inverse Comptonization of the synchrotron photons (from the layer) by the spine electrons.

Blazars are classified into high, intermediate and low energy peaked sources based on the location of their synchrotron peak. This lies in infra-red/optical to ultra-violet bands for low and intermediate peaked blazars. The transition from synchrotron to inverse Compton emission falls in the X-ray bands for such sources. We present the spectral and timing analysis of 14 low and intermediate energy peaked blazars ob- served with XMMNewton spanning 31 epochs. Parametric fits to X-ray spectra helps constrain the possible location of transition from the high energy end of the syn- chrotron to the low energy end of the inverse Compton emission. In seven sources in our sample, we infer such a transition and constrain the break energy in the range 0.6 10 keV. The Lomb-Scargle periodogram is used to estimate the power spectral density (PSD) shape. It is well described by a power law in a majority of light curves, the index being flatter compared to general expectation from AGN, ranging here between 0.01 and 1.12, possibly due to short observation durations resulting in an absence of long term trends. A toy model involving synchrotron self-Compton (SSC) and exter- nal Compton (EC; disk, broad line region, torus) mechanisms are used to estimate magnetic field strength 6 0.03 - 0.88 G in sources displaying the energy break and infer a prominent EC contribution. The timescale for variability being shorter than synchrotron cooling implies steeper PSD slopes which are inferred in these sources.

We present a multi-wavelength spectral and temporal analysis of the blazar OJ 287 during its recent activity between December 2015 -- May 2016, showing strong variability in the near-infrared (NIR) to X-ray energies with detection at $\gamma$-ray energies as well. Most of the optical flux variations exhibit strong changes in polarization angle and degree. All the inter-band time lags are consistent with simultaneous emissions. Interestingly, on days with excellent data coverage in the NIR--UV bands, the spectral energy distributions (SEDs) show signatures of bumps in the visible--UV bands, never seen before in this source. The optical bump can be explained as accretion-disk emission associated with the primary black hole of mass $\sim \rm 1.8 \times10^{10} M_{\odot}$ while the little bump feature in the optical-UV appears consistent with line emission. Further, the broadband SEDs extracted during the first flare and during a quiescent period during this span show very different $\gamma$-ray spectra compared to previously observed flare or quiescent spectra. The probable thermal bump in the visible seems to have been clearly present since May 2013, as found by examining all available NIR-optical observations, and favors the binary super-massive black hole model. The simultaneous multi-wavelength variability and relatively weak $\gamma$-ray emission that shows a shift in the SED peak is consistent with $\gamma$-ray emission originating from inverse Compton scattering of photons from the line emission that apparently contributes to the little blue bump.

The blazar 3C454.3 exhibited a strong flare seen in gamma-rays, X-rays, and optical/NIR bands during 3--12 December 2009. Emission in the V and J bands rose more gradually than did the gamma-rays and soft X-rays, though all peaked at nearly the same time. Optical polarization measurements showed dramatic changes during the flare, with a strong anti-correlation between optical flux and degree of polarization (which rose from ~ 3% to ~ 20%) during the declining phase of the flare. The flare was accompanied by large rapid swings in polarization angle of ~ 170 degree. This combination of behaviors appear to be unique. We have cm-band radio data during the same period but they show no correlation with variations at higher frequencies. Such peculiar behavior may be explained using jet models incorporating fully relativistic effects with a dominant source region moving along a helical path or by a shock-in-jet model incorporating three-dimensional radiation transfer if there is a dominant helical magnetic field. We find that spectral energy distributions at different times during the flare can be fit using modified one-zone models where only the magnetic field strength and particle break frequencies and normalizations need change. An optical spectrum taken at nearly the same time provides an estimate for the central black hole mass of ~ 2.3 * 10^9 M_sun. We also consider two weaker flares seen during the $\sim 200$ d span over which multi-band data are available. In one of them, the V and J bands appear to lead the $\gamma$-ray and X-ray bands by a few days; in the other, all variations are simultaneous.

We present recent optical photometric observations of the blazar OJ 287 taken during September 2015 -- May 2016. Our intense observations of the blazar started in November 2015 and continued until May 2016 and included detection of the large optical outburst in December 2016 that was predicted using the binary black hole model for OJ 287. For our observing campaign, we used a total of 9 ground based optical telescopes of which one is in Japan, one is in India, three are in Bulgaria, one is in Serbia, one is in Georgia, and two are in the USA. These observations were carried out in 102 nights with a total of ~ 1000 image frames in BVRI bands, though the majority were in the R band. We detected a second comparably strong flare in March 2016. In addition, we investigated multi-band flux variations, colour variations, and spectral changes in the blazar on diverse timescales as they are useful in understanding the emission mechanisms. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux, colour and spectral variability.

Blazars can be divided into two sub-classes namely high energy and low energy peaked blazars. In spectral energy distribution, the first synchrotron hump of the former class peaks in UV/X-rays and in IR/optical bands for the latter class. The peak of the spectral energy distribution seems to be responsible for variability properties of these classes of blazars in X-ray and optical bands. Since, in low energy peaked blazars, the X-ray bands lies well below the synchrotron hump, one expects that the highest energy electrons available for the synchrotron emission would have slower effect of variability on X-ray intra-day timescale. In this paper, by taking the advantage of a sample of 12 low energy peaked blazars with total 50 observations from XMM$-$Newton since its launch, we confirm that this class is less variable in X-ray bands. We found that out of 50 observational light curves, genuine intra-day variability is present in only two of light curves i.e 4%. Similar results we obtained from our earlier optical intra-day variability studies of high energy peaked blazars where out of 144 light curves, only genuine intra-day variability was detected in 6 light curves i.e ~ 4%. Since, X-ray bands lie below the peak of the spectral energy distribution of LSPs where inverse Compton mechanism is dominating rather than synchrotron radiation at the peak of the optical band, leads to slower variability in the X-ray bands. Hence, reducing their intra-day variability in X-ray bands as compared to the variability in optical bands.

We present our optical photometric observations of three TeV blazars, PKS 1510-089, PG 1553+113 and Mrk 501 taken using two telescopes in India, one in Bulgaria, one in Greece and one in Serbia during 2012 - 2014. These observations covered a total of 95 nights with a total of 202 B filter frames, 247 images in V band, 817 in R band while 229 images were taken in the I filter. This work is focused on multi-band flux and colour variability studies of these blazars on diverse timescales which are useful in understanding the emission mechanisms. We studied the variability characteristics of above three blazars and found all to be active over our entire observational campaigns. We also searched for any correlation between the brightness of the sources and their colour indices. During the times of variability, no significant evidence for the sources to display spectral changes correlated with magnitude was found on timescales of a few months. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux variability.

We present the results of extensive multi-band intra-night optical monitoring of BL Lacertae during 2010--2012. BL Lacertae was very active in this period and showed intense variability in almost all wavelengths. We extensively observed it for a total for 38 nights; on 26 of them observations were done quasi-simultaneously in B, V, R and I bands (totaling 113 light curves), with an average sampling interval of around 8 minutes. BL Lacertae showed significant variations on hour-like timescales in a total of 19 nights in different optical bands. We did not find any evidence for periodicities or characteristic variability time-scales in the light curves. The intranight variability amplitude is generally greater at higher frequencies and decreases as the source flux increases. We found spectral variations in BL Lacertae in the sense that the optical spectrum becomes flatter as the flux increases but in several flaring states deviates from the linear trend suggesting different jet components contributing to the emission at different times.

We observed the prototype blazar, BL Lacertae, extensively in optical and radio bands during an active phase in the period 2010--2013 when the source showed several prominent outbursts. We searched for possible correlations and time lags between the optical and radio band flux variations using multifrequency data to learn about the mechanisms producing variability. During an active phase of BL Lacertae, we searched for possible correlations and time lags between multifrequency light curves of several optical and radio bands. We tried to estimate any possible variability timescales and inter-band lags in these bands. We performed optical observations in B, V, R and I bands from seven telescopes in Bulgaria, Georgia, Greece and India and obtained radio data at 36.8, 22.2, 14.5, 8 and 4.8 GHz frequencies from three telescopes in Ukraine, Finland and USA. Significant cross-correlations between optical and radio bands are found in our observations with a delay of cm-fluxes with respect to optical ones of ~250 days. The optical and radio light curves do not show any significant timescales of variability. BL Lacertae showed many optical 'mini-flares' on short time-scales. Variations on longer term timescales are mildly chromatic with superposition of many strong optical outbursts. In radio bands, the amplitude of variability is frequency dependent. Flux variations at higher radio frequencies lead the lower frequencies by days or weeks. The optical variations are consistent with being dominated by a geometric scenario where a region of emitting plasma moves along a helical path in a relativistic jet. The frequency dependence of the variability amplitude supports an origin of the observed variations intrinsic to the source.

To search for optical variability on a wide range of timescales, we have carried out photometric monitoring of two flat spectrum radio quasars, 3C 454.3 and 3C 279, plus one BL Lac, S5 0716+714, all of which have been exhibiting remarkably high activity and pronounced variability at all wavelengths. CCD magnitudes in B, V, R and I pass-bands were determined for $\sim$ 7000 new optical observations from 114 nights made during 2011 - 2014, with an average length of $\sim$ 4 h each, at seven optical telescopes: four in Bulgaria, one in Greece, and two in India. We measured multiband optical flux and colour variations on diverse timescales. Discrete correlation functions were computed among B, V, R, and I observations, to search for any time delays. We found weak correlations in some cases with no significant time lags. The structure function method was used to estimate any characteristic time-scales of variability. We also investigated the spectral energy distribution of the three blazars using B, V, R, I, J and K pass-band data. We found that the sources almost always follows a bluer-when-brighter trend. We discuss possible physical causes of the observed spectral variability.

We present the results of photometric (V band) and polarimetric observations of the blazar BL Lac during 2008--2010 using TRISPEC attached to the KANATA 1.5-m telescope in Japan. The data reveal a great deal of variability ranging from days to months with detection of strong variations in fractional polarization. The V band flux strongly anti-correlates with the degree of polarization during the first of two observing seasons but not during the second. The direction of the electric vector, however, remained roughly constant during all our observations. These results are consistent with a model with at least two emission regions being present, with the more variable component having a polarization direction nearly perpendicular to that of the relatively quiescent region so that a rising flux can produce a decline in degree of polarization. We also computed models involving helical jet structures and single transverse shocks in jets and show that they might also be able to agree with the anti-correlations between flux and fractional polarization.

We report the results of optical monitoring for a sample of 11 blazars including 10 BL Lacs and 1 Flat Spectrum Radio Quasar (FSRQ). We have measured the multiband optical flux and colour variations in these blazars on intra-day and short-term timescales of months and have limited data for 2 more blazars. These photometric observations were made during 2009 to 2011, using six optical telescopes, four in Bulgaria, one in Greece and one in India. On short-term timescales we found significant flux variations in 9 of the sources and colour variations in 3 of them. Intra-day variability was detected on 6 nights for 2 sources out of the 18 nights and 4 sources for which we collected such data. These new optical observations of these blazars plus data from our previous published papers (for 3 more blazars) were used to analyze their spectral flux distributions in the optical frequency range. Our full sample for this purpose includes 6 high-synchrotron-frequency-peaked BL Lacs (HSPs), 3 intermediate-synchrotron-frequency-peaked BL Lacs (ISPs) and 6 low-synchrotron-frequency-peaked BL Lacs (LSPs; including both BL Lacs and FSRQs). We also investigated the spectral slope variability and found that the average spectral slopes of LSPs show a good accordance with the Synchrotron Self-Compton (SSC) loss dominated model. Our analysis supports previous studies that found that the spectra of the HSPs and FSRQs have significant additional emission components. The spectra of all these HSPs and LSPs get flatter when they become brighter, while for FSRQs the opposite appears to hold. This supports the hypothesis that there is a significant thermal contribution to the optical spectrum for FSRQs.

In this paper we present results of a short-term optical monitoring of 13 blazars. The objects were monitored mostly in the R-band for a total of ~ 160 hours between 2006 and 2011. We study the nature of the short-term variations and show that most of them could be described as slow, smooth, and (almost) linear changes of up to ~ 0.1 mag/hour, but many objects show no short-term variations at all. In fact, we found only ~ 2 per cent chance to observe variability of more than 0.1 mag/hour for the sample we observed. Hints for quasi-periodic oscillations at very low amplitude levels are also found for some objects. We briefly discuss some of the possible mechanisms to generate the intra-night variability and the quasi-periodic oscillations.

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

We report the results of photometric observations of the blazars Mrk 421 and 3C 454.3 designed to search for intraday variability (IDV) and short-term variability (STV). Optical photometric observations were spread over eighteen nights for Mrk 421 and seven nights for 3C 454.3 during our observing run in 2009-2010 at the 1.04 m telescope at ARIES, India. Genuine IDV is found for the source 3C 454.3 but not for Mrk 421. Genuine STV is found for both sources. Mrk 421 was revealed by the MAXI X-ray detector on the International Space Station to be in an exceptionally high flux state in 2010 January - February. We performed a correlation between the X-ray and optical bands to search for time delays and found a weak correlation with higher frequencies leading the lower frequencies by about ten days. The blazar 3C 454.3 was found to be in high flux state in November-December 2009. We performed correlations in optical observations made at three telescopes, along with X-ray data from the MAXI satellite and public release gamma-ray data from the Fermi space telescope. We found strong correlations between the gamma-ray and optical bands at a time lag of about four days but the X-ray flux is not correlated with either. We briefly discuss the possible reasons for the time delays between these bands within the framework of existing models for X-ray and gamma-ray emission mechanisms.

We report the results of quasi-simultaneous two filter optical monitoring of two high-energy peaked blazars, 1ES 1959+650 and 1ES 2344+514, to search for microvariability and short-term variability (STV). We carried out optical photometric monitoring of these sources in an alternating sequence of B and R pass-bands, and have 24 and 19 nights of new data for these two sources, respectively. No genuine microvariability (intra-night variability) was detected in either of these sources. This non-detection of intra-night variations is in agreement with the conclusions of previous studies that high-energy peaked BL Lacs are intrinsically less variable than low-energy peaked BL Lacs in the optical bands. We also report the results of STV studies for these two sources between July 2009 and August 2010. Genuine STV is found for the source 1ES 1959+650 but not for 1ES 2344+514. We briefly discuss possible reasons for the difference between the intra-night variability behaviour of high- and low-energy peaked blazars.

We selected a sample of 24 XMM-Newton light curves (LCs) of four high energy peaked blazars, PKS 0548-322, ON 231, 1ES 1426+428 and PKS 2155-304. These data comprise continuous light curves of 7.67h to 18.97h in length. We searched for possible quasi-periodic oscillations (QPO) and intra-day variability (IDV) timescales in the LCs of these blazars. We found a likely QPO in one LC of PKS 2155-304 which was reported elsewhere (Lachowicz et al. 2009). In the remaining 23 LCs we found hints of possible weak QPOs in one LC of each of ON 231 and PKS 2155-304, but neither is statistically significant. We found IDV timescales that ranged from 15.7 ks to 46.8 ks in 8 LCs. In 13 LCs any variability timescales were longer than the length of the data. Assuming the possible weak QPO periods in the blazars PKS 2155-304 and ON 231 are real and are associated with the innermost portions of their accretion disk, we can estimate that their central black hole masses exceed 1.2 $\times$ 10$^{7}$ M$_{\odot}$. Emission models for radio-loud active galactic nuclei (AGN) that could explain our results are briefly discussed.

We report a possible detection of an ~4.6-hour quasi-periodic oscillation (QPO) in the 0.3-10 keV emission of the high-energy peaked blazar PKS 2155-304 from a 64 ks observation by the XMM-Newton EPIC/pn detector. We identify a total modulation of ~5% in the light curve and confirm that nominal period by periodogram, structure function and wavelet analyses. The limited light curve duration allows the capture of only 3.8 cycles of this oscillation and thus precludes a very strong claim for this QPO, despite a nominally high (>3 sigma) statistical significance. We briefly discuss models capable of producing an X-ray QPO of such a period in a blazar.