Source author record

I. Liodakis

I. Liodakis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2026arXiv

Polarization echoes from past nuclear activity in the quasi-periodic eruption source GSN 069

Context. X-ray quasi-periodic eruptions (QPEs) are repeating, high-amplitude, soft X-ray bursts observed from the nuclei of a dozen nearby low-mass galaxies. Their origin remains a puzzle in the physics of accretion variability. Observational data indicate that X-ray and/or optical tidal disruption events (TDEs) may precede QPE detections. Although both kinds of outburst are driven by supermassive black holes, they are more frequently detected in faded active galactic nuclei (AGNs), when the TDE is not happening in a dormant galaxy. In the case of the QPE discovery source, GSN 069, it remains debated whether its past activity arose from a previous AGN phase or from an enhanced TDE rate. Aims. We investigated the origin of the past nuclear activity in GSN 069. Methods. Past AGN activity imprints detectable polarization in optical light, due to the expected delay between direct and scattered light. On 6 September 2019, we targeted GSN 069 with VLT/FORS2 in both imaging polarimetry and spectropolarimetry modes so that its optical polarization could be investigated while the first detected QPE phase was still active. Results. We measured a rising polarization, from ~0% to ~1.5%, as moving away from the nucleus of GSN 069. This rise is probed to be intrinsic to the central engine, confirming the already detected extended emission line region (EELR) by integral field unit data. Conclusions. The increasing radial polarization demonstrates a switched-off nucleus. The polarization angle traces an axis aligned with elongated [OIII], [NII], and Hα gas distributions, revealing an EELR that may be consistent with relic polarization cones, therefore suggesting the presence of a torus-like structure in the past. Thus, optical polarization echoes geometrically favor a faded AGN as the origin of the EELR rather than a past elevated TDE rate, although the latter cannot be excluded.

preprint2023arXiv

Optical/$γ$-ray blazar flare correlations: understanding the high-energy emission process using ASAS-SN and Fermi light curves

Using blazar light curves from the optical All-Sky Automated Survey for Supernovae (ASAS-SN) and the $γ$-ray \textit{Fermi}-LAT telescope, we performed the most extensive statistical correlation study between both bands, using a sample of 1,180 blazars. This is almost an order of magnitude larger than other recent studies. Blazars represent more than 98\% of the AGNs detected by \textit{Fermi}-LAT and are the brightest $γ$-ray sources in the extragalactic sky. They are essential for studying the physical properties of astrophysical jets from central black holes. However, their $γ$-ray flare mechanism is not fully understood. Multi-wavelength correlations help constrain the dominant mechanisms of blazar variability. We search for temporal relationships between optical and $γ$-ray bands. Using a Bayesian Block Decomposition, we detect 1414 optical and 510 $γ$-ray flares, we find a strong correlation between both bands. Among all the flares, we find 321 correlated flares from 133 blazars, and derive an average rest-frame time delay of only 1.1$_{-8.5}^{+7.1}$ days, with no difference between the flat-spectrum radio quasars, BL Lacertae-like objects or low, intermediate, and high-synchrotron peaked blazar classes. Our time-delay limit rules out the hadronic proton-synchrotron model as the driver for non-orphan flares and suggests a leptonic single-zone model. Limiting our search to well-defined light curves and removing 976 potential but unclear ``orphan'' flares, we find 191 (13\%) and 115 (22\%) clear ``orphan'' optical and $γ$-ray flares. The presence of ``orphan'' flares in both bands challenges the standard one-zone blazar flare leptonic model and suggests multi-zone synchrotron sites or a hadronic model for some blazars.

preprint2022arXiv

Accretion geometry of the neutron star low mass X-ray binary Cyg X-2 from X-ray polarization measurements

We report spectro-polarimetric results of an observational campaign of the bright neutron star low-mass X-ray binary Cyg X-2 simultaneously observed by IXPE, NICER and INTEGRAL. Consistently with previous results, the broad-band spectrum is characterized by a lower-energy component, attributed to the accretion disc with $kT_{\rm in} \approx$ 1 keV, plus unsaturated Comptonization in thermal plasma with temperature $kT_{\rm e} = 3$ keV and optical depth $τ\approx 4$, assuming a slab geometry. We measure the polarization degree in the 2-8 keV band $P=1.8 \pm 0.3$ per cent and polarization angle $ϕ= 140^{\circ} \pm 4^{\circ}$, consistent with the previous X-ray polarimetric measurements by OSO-8 as well as with the direction of the radio jet which was earlier observed from the source. While polarization of the disc spectral component is poorly constrained with the IXPE data, the Comptonized emission has a polarization degree $P =4.0 \pm 0.7$ per cent and a polarization angle aligned with the radio jet. Our results strongly favour a spreading layer at the neutron star surface as the main source of the polarization signal. However, we cannot exclude a significant contribution from reflection off the accretion disc, as indicated by the presence of the iron fluorescence line.

preprint2022arXiv

Detecting intermediate mass black holes in midiquasars with current and future surveys

The lack of detected intermediate mass black holes poses a gap in our understanding of the growth and evolution of the most exotic of astrophysical objects. Here we investigate the possibility of low-luminosity relativistic jets launched by intermediate mass black holes in the centers of dwarf galaxies. We built population models that allow us to make predictions for their radio emission and quantify their detectability by current and future surveys. We find that the upcoming instruments in optical and radio like the SKA, ngVLA and the Vera C. Rubin Observatory will likely be able to detect a significant fraction ($>38\%$) of such sources population if they exist. In addition, our results suggest that it is not unlikely a small number of midiquasars possibly masquerading as low-luminosity active galactic nuclei may have already been detected by existing surveys.

preprint2022arXiv

Investigating the blazar TXS 0506+056 through sharp multi-wavelength eyes during 2017-2019

The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $σ$ level to a $γ$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the IceCube-170922A event. To better characterize its broad-band emission, we organized a multi-wavelength campaign lasting 16 months (November 2017 to February 2019), covering the radio-band (Metsähovi, OVRO), the optical/UV (ASAS-SN, KVA, REM, Swift/UVOT), the X-rays (Swift/XRT, NuSTAR), the high-energy $γ$ rays (Fermi/LAT) and the very-high-energy (VHE) $γ$ rays (MAGIC). In $γ$ rays, the behaviour of the source was significantly different from the 2017 one: MAGIC observations show the presence of flaring activity during December 2018, while the source only shows an excess at the 4$σ$ level during the rest of the campaign (74 hours of accumulated exposure); Fermi/LAT observations show several short (days-to-week timescale) flares, different from the long-term brightening of 2017. No significant flares are detected at lower energies. The radio light curve shows an increasing flux trend, not seen in other wavelengths. We model the multi-wavelength spectral energy distributions in a lepto-hadronic scenario, in which the hadronic emission emerges as Bethe-Heitler and pion-decay cascade in the X-rays and VHE $γ$ rays. According to the model presented here, the December 2018 $γ$-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments.

preprint2022arXiv

New Tests of Millilensing in the Blazar PKS 1413+135

Symmetric Achromatic Variability (SAV) is a rare form of radio variability in blazars that has been attributed to gravitational millilensing by a ~$10^2 - 10^5$ $M_\odot$ mass condensate. Four SAVs have been identified between 1980 and 2020 in the long-term radio monitoring data of the blazar PKS 1413+135. We show that all four can be fitted with the same, unchanging, gravitational lens model. If SAV is due to gravitational millilensing, PKS 1413+135 provides a unique system for studying active galactic nuclei with unprecedented microarcsecond resolution, as well as for studying the nature of the millilens itself. We discuss two possible candidates for the putative millilens: a giant molecular cloud hosted in the intervening edge-on spiral galaxy, and an undetected dwarf galaxy with a massive black hole. We find a significant dependence of SAV crossing time on frequency, which could indicate a fast shock moving in a slower underlying flow. We also find tentative evidence for a 989-day periodicity in the SAVs, which, if real, makes possible the prediction of future SAVs: the next three windows for possible SAVs begin in August 2022, May 2025, and February 2028.

preprint2022arXiv

The Unanticipated Phenomenology of the Blazar PKS~2131$-$021: A Unique Super-Massive Black Hole Binary Candidate

Most large galaxies host supermassive black holes in their nuclei and are subject to mergers, which can produce a supermassive black hole binary (SMBHB), and hence periodic signatures due to orbital motion. We report unique periodic radio flux density variations in the blazar PKS~2131$-$021, which strongly suggest an SMBHB with an orbital separation of $\sim 0.001-0.01$ pc. Our 45.1-year radio light curve shows two epochs of strong sinusoidal variation with the same period and phase to within $<2\%$ and $\sim 10\%$, respectively, straddling a 20-year period when this variation was absent. Our simulated light curves accurately reproduce the ``red noise'' of this object, and Lomb-Scargle, weighted wavelet Z-transform, and least-squares sine wave analyses demonstrate conclusively, at the $4.6σ$ significance level, that the periodicity in this object is not due to random fluctuations in flux density. The observed period translates to $2.082\pm 0.003$ years in the rest frame at the $z=1.285$ redshift of PKS~2131$-$021. The periodic variation in PKS~2131$-$021 is remarkably sinusoidal. We present a model in which orbital motion, combined with the strong Doppler boosting of the approaching relativistic jet, produces a sine-wave modulation in the flux density which easily fits the observations. Given the rapidly-developing field of gravitational wave experiments with pulsar timing arrays, closer counterparts to PKS~2131$-$021 and searches using the techniques we have developed are strongly motivated. These results constitute a compelling demonstration that the phenomenology, not the theory, must provide the lead in this field.

preprint2020arXiv

Fermi and Swift Observations of GRB 190114C: Tracing the Evolution of High-Energy Emission from Prompt to Afterglow

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

preprint2020arXiv

Proton Synchrotron $γ$-rays and the Energy Crisis in Blazars

The origin of high-energy emission in blazars jets (i.e., leptonic versus hadronic) has been a long-standing matter of debate. Here, we focus on one variant of hadronic models where proton synchrotron radiation accounts for the observed steady $γ$-ray blazar emission. Using analytical methods, we derive the minimum jet power ($P_{j,\min}$) for the largest blazar sample analyzed to date (145 sources), taking into account uncertainties of observables and jet's physical parameters. We compare $P_{j,\min}$ against three characteristic energy estimators for accreting systems, i.e., the Eddington luminosity, the accretion disk luminosity, and the power of the Blandford-Znajek process, and find that $P_{j,\min}$ is about 2 orders of magnitude higher than all energetic estimators for the majority of our sample. The derived magnetic field strengths in the emission region require either large amplification of the jet's magnetic field (factor of 30) or place the $γ$-ray production site at sub-pc scales. The expected neutrino emission peaks at $\sim 0.1-10$ EeV, with typical peak neutrino fluxes $\sim 10^{-4}$ times lower than the peak $γ$-ray fluxes. We conclude that if relativistic hadrons are present in blazar jets, they can only produce a radiatively subdominant component of the overall spectral energy distribution of the blazar's steady emission.

preprint2020arXiv

The Relativistic Jet Orientation and Host Galaxy of the Peculiar Blazar PKS 1413+135

PKS 1413+135 is one of the most peculiar blazars known. Its strange properties led to the hypothesis almost four decades ago that it is gravitationally lensed by a mass concentration associated with an intervening galaxy. It exhibits symmetric achromatic variability, a rare form of variability that has been attributed to gravitational milli-lensing. It has been classified as a BL Lac object, and is one of the rare objects in this class with a visible counterjet. BL Lac objects have jet axes aligned close to the line of sight. It has also been classified as a compact symmetric object, which have jet axes not aligned close to the line of sight. Intensive efforts to understand this blazar have hitherto failed to resolve even the questions of the orientation of the relativistic jet, and the host galaxy. Answering these two questions is important as they challenge our understanding of jets in active galactic nuclei and the classification schemes we use to describe them. We show that the jet axis is aligned close to the line of sight and PKS 1413+135 is almost certainly not located in the apparent host galaxy, but is a background object in the redshift range $0.247 < z < 0.5$. The intervening spiral galaxy at $z = 0.247$ provides a natural host for the putative lens responsible for symmetric achromatic variability and is shown to be a Seyfert 2 galaxy. We also show that, as for the radio emission, a "multizone" model is needed to account for the high-energy emission.

preprint2016arXiv

Detecting the Elusive Blazar Counter-Jets

Detection of blazar pc scale counter-jets is difficult, but it can provide invaluable insight into the relativistic effects, radiative processes and the complex mechanisms of jet production, collimation and accelation in blazars. We build on recent populations models (optimized using the MOJAVE apparent velocity and redshift distributions) in order to derive the distribution of jet-to-counter-jet ratios and the flux densities of the counter-jet at different frequencies, in an effort to set minimum sensitivity limits required for existing and future telescope arrays in order to detect these elusive counter-jets. We find that: for the BL Lacs $5\%$ of their counter-jets have a flux-density higher than 100mJy, $15\%$ are higher than 10 mJy, and $32\%$ have higher flux-density than 1 mJy, whereas for the FSRQs $8\%$ have a flux-density higher than 10mJy, $17\%$ are higher than 1 mJy, and $32\%$ are higher than 0.1 mJy (at 15 GHz). Future telescopes like the SKA and newly operating like e-MERLIN and JVLA may detect up to $99\%$ of the BL Lac and $77\%$ of the FSRQ counter-jets. Sources with both low apparent velocity and a low Doppler factor make prime candidates for counter-jet detection. Combining our findings with literature values we have identified five such counter-jet detection candidates. Finally, we discuss possible effects beyond relativistic deboosting that may complicate the detection of counter-jets and that need to be accounted for in the interpretation of detections.

preprint2016arXiv

Estimating the distribution of rest-frame timescales for blazar jets: a statistical approach

In any flux-density limited sample of blazars, the distribution of the timescale modulation factor $Δt'/Δt$, which quantifies the change in observed timescales compared to the rest-frame ones due to redshift and relativistic compression follows an exponential distribution with a mean depending on the flux limit of the sample. In this work we produce the mathematical formalism that allows us to use this information in order to uncover the underlining rest-frame probability density function of measurable timescales of blazar jets. We extensively test our proposed methodology using a simulated FSRQ population with a 1.5 Jy flux-density limit in the simple case (where all blazars share the same intrinsic timescale), in order to identify limits of applicability and potential biases due to observational systematics and sample selection. We find that for monitoring with time intervals between observations longer than $\sim$30\% of the intrinsic timescale under investigation the method loses its ability to produce robust results. For time intervals of $\sim$3\% of the intrinsic timescale the error of the method is as low as 1\% in recovering the intrinsic rest-frame timescale. We applied our method to rotations of the optical polarization angle of blazars observed by RoboPol. We found that the intrinsic timescales of the longest-duration rotation event in each blazar follows a narrow distribution, well-described by a normal distribution with mean 87 days and standard deviation 5 days. We discuss possible interpretations of this result.

preprint2016arXiv

Optical polarization map of the Polaris Flare with RoboPol

The stages before the formation of stars in molecular clouds are poorly understood. Insights can be gained by studying the properties of quiescent clouds, such as their magnetic field structure. The plane-of-the-sky orientation of the field can be traced by polarized starlight. We present the first extended, wide-field ($\sim$10 $\rm deg^2$) map of the Polaris Flare cloud in dust-absorption induced optical polarization of background stars, using the RoboPol polarimeter at the Skinakas Observatory. This is the first application of the wide-field imaging capabilities of RoboPol. The data were taken in the R-band and analysed with the automated reduction pipeline of the instrument. We present in detail optimizations in the reduction pipeline specific to wide-field observations. Our analysis resulted in reliable measurements of 641 stars with median fractional linear polarization 1.3%. The projected magnetic field shows a large scale ordered pattern. At high longitudes it appears to align with faint striations seen in the Herschel-SPIRE map of dust emission (250 $μm$), while in the central 4-5 deg$^2$ it shows an eddy-like feature. The overall polarization pattern we obtain is in good agreement with large scale measurements by Planck of the dust emission polarization in the same area of the sky.

preprint2016arXiv

Optical polarization of high-energy BL Lac objects

We investigate the optical polarization properties of high-energy BL Lac objects using data from the RoboPol blazar monitoring program and the Nordic Optical Telescope. We wish to understand if there are differences in the BL Lac objects that are detected with the current-generation TeV instruments compared to those that have not yet been detected. The mean polarization fraction of the TeV-detected BL Lacs is 5% while the non-TeV sources show a higher mean polarization fraction of 7%. This difference in polarization fraction disappears when the dilution by the unpolarized light of the host galaxy is accounted for. The TeV sources show somewhat lower fractional polarization variability amplitudes than the non-TeV sources. Also the fraction of sources with a smaller spread in the Q/I - U/I -plane and a clumped distribution of points away from the origin, possibly indicating a preferred polarization angle, is larger in the TeV than in the non-TeV sources. These differences between TeV and non-TeV samples seems to arise from differences between intermediate and high spectral peaking sources instead of the TeV detection. When the EVPA variations are studied, the rate of EVPA change is similar in both samples. We detect significant EVPA rotations in both TeV and non-TeV sources, showing that rotations can occur in high spectral peaking BL Lac objects when the monitoring cadence is dense enough. Our simulations show that we cannot exclude a random walk origin for these rotations. These results indicate that there are no intrinsic differences in the polarization properties of the TeV-detected and non-TeV-detected high-energy BL Lac objects. This suggests that the polarization properties are not directly related to the TeV-detection, but instead the TeV loudness is connected to the general flaring activity, redshift, and the synchrotron peak location. (Abridged)

preprint2016arXiv

RoboPol: Do optical polarization rotations occur in all blazars?

We present a new set of optical polarization plane rotations in blazars, observed during the third year of operation of RoboPol. The entire set of rotation events discovered during three years of observations is analysed with the aim of determining whether these events are inherent in all blazars. It is found that the frequency of the polarization plane rotations varies widely among blazars. This variation cannot be explained either by a difference in the relativistic boosting or by selection effects caused by a difference in the average fractional polarization. We conclude that the rotations are characteristic of a subset of blazars and that they occur as a consequence of their intrinsic properties.

preprint2016arXiv

RoboPol: First season rotations of optical polarization plane in blazars

We present first results on polarization swings in optical emission of blazars obtained by RoboPol, a monitoring program of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events. A possible connection of polarization swing events with periods of high activity in gamma rays is investigated using the dataset obtained during the first season of operation. It was found that the brightest gamma-ray flares tend to be located closer in time to rotation events, which may be an indication of two separate mechanisms responsible for the rotations. Blazars with detected rotations have significantly larger amplitude and faster variations of polarization angle in optical than blazars without rotations. Our simulations show that the full set of observed rotations is not a likely outcome (probability $\le 1.5 \times 10^{-2}$) of a random walk of the polarization vector simulated by a multicell model. Furthermore, it is highly unlikely ($\sim 5 \times 10^{-5}$) that none of our rotations is physically connected with an increase in gamma-ray activity.

preprint2016arXiv

RoboPol: optical polarization-plane rotations and flaring activity in blazars

We present measurements of rotations of the optical polarization of blazars during the second year of operation of RoboPol, a monitoring programme of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events, and we analyse the large set of rotation events discovered in two years of observation. We investigate patterns of variability in the polarization parameters and total flux density during the rotation events and compare them to the behaviour in a non-rotating state. We have searched for possible correlations between average parameters of the polarization-plane rotations and average parameters of polarization, with the following results: (1) there is no statistical association of the rotations with contemporaneous optical flares; (2) the average fractional polarization during the rotations tends to be lower than that in a non-rotating state; (3) the average fractional polarization during rotations is correlated with the rotation rate of the polarization plane in the jet rest frame; (4) it is likely that distributions of amplitudes and durations of the rotations have physical upper bounds, so arbitrarily long rotations are not realised in nature.

preprint2016arXiv

RoboPol: The optical polarization of gamma-ray--loud and gamma-ray--quiet blazars

We present average R-band optopolarimetric data, as well as variability parameters, from the first and second RoboPol observing season. We investigate whether gamma- ray--loud and gamma-ray--quiet blazars exhibit systematic differences in their optical polarization properties. We find that gamma-ray--loud blazars have a systematically higher polarization fraction (0.092) than gamma-ray--quiet blazars (0.031), with the hypothesis of the two samples being drawn from the same distribution of polarization fractions being rejected at the 3σ level. We have not found any evidence that this discrepancy is related to differences in the redshift distribution, rest-frame R-band lu- minosity density, or the source classification. The median polarization fraction versus synchrotron-peak-frequency plot shows an envelope implying that high synchrotron- peaked sources have a smaller range of median polarization fractions concentrated around lower values. Our gamma-ray--quiet sources show similar median polarization fractions although they are all low synchrotron-peaked. We also find that the random- ness of the polarization angle depends on the synchrotron peak frequency. For high synchrotron-peaked sources it tends to concentrate around preferred directions while for low synchrotron-peaked sources it is more variable and less likely to have a pre- ferred direction. We propose a scenario which mediates efficient particle acceleration in shocks and increases the helical B-field component immediately downstream of the shock.

preprint2015arXiv

Population statistics of beamed sources. I: A new model for blazars

Observations of blazar jets are shrouded in relativistic effects, thus hindering our understanding of their intrinsic properties and dominant physical processes responsible for their generation, evolution,radiation and particle emission. In this work we focus on extracting information about timescales in the jet rest frame using a population-modeling approach. We employ Monte Carlo simulations to derive a simple population model for the intrinsic unbeamed luminosities and the Lorentz $Γ$ of blazar jets that adequately describe the observed redshift and apparent superluminal speed distributions for flux-limited blazar samples. We derive separate models for BL Lacs and Flat Spectrum Radio Quasars. We then use these models to compute the predicted distribution of Doppler factors in each blazar class, and address the following questions: (a) What is the relativistically induced spread in observed timescales (e.g., event duration, time lags) in a flux-limited sample of relativistic jets? (b) Could differences between BL Lacs and FSRQs observed in the time domain be attributed to differences in beaming between the two populations? (c) Is there a statistically preferred amount of beaming in a flux-limited sample? How large are statistical deviations from that preferred value? We use our findings to propose promising approaches in phenomenological studies of timescales in blazar jets.

preprint2015arXiv

Population statistics of beamed sources. II: Evaluation of Doppler factor estimates

In a companion paper we presented a statistical model for the blazar population, consisting of distributions for the unbeamed radio luminosity function and the Lorentz factor distribution of each of the BL Lac and Flat Spectrum Radio Quasar (FSRQ) classes. Our model has been optimized so that it reproduces the MOJAVE distributions of apparent speeds and redshifts when the appropriate flux limit is applied and a uniform distribution of jet viewing angles is assumed for the population. Here we use this model to predict the Doppler factor distribution for various flux-limited samples for which Doppler factors have been estimated in a variety of ways (equipartition, variability + equipartition, inverse Compton dominance) on a blazar-by-blazar basis. By comparing the simulated and data-estimated Doppler factor distributions in each case, we evaluate the different methods of estimating blazar Doppler factors. We find that the variability Doppler factors assuming equipartition are the ones in the best agreement with our statistical model, whereas the inverse Compton Doppler factor method is only suitable for FSRQs. In the case of variability Doppler factors, we find that while random errors are relatively low ($\sim 30\%$), uncertainties are dominated by systematic effects. In the case of inverse Compton Doppler factors, random errors appear to dominate, but are significantly larger ($\sim 60\%$).

preprint2014arXiv

Jet bending and velocity profiles

In this study, we consider the problem of geometry of curved jets. We work with a simplified two-dimensional kinematical model assuming that the jet velocity is varying. We prescribe the velocity variation and study two cases concerning parsec and kilo-parsec scale jets by applying some simplified numerical simulations. For parsec scale jets, we use some general expressions and derive the shape of jets for certain numerical cases. For kilo-parsec scale jets, we study the cases of 4C53.37 and 1159+583 of Abell 2220 and Abell 1446 galaxy clusters, respectively, in an effort to test our model and compare our numerical results with observational data.