Recently the IceCube collaboration detected very high energy (VHE) neutrinos and associated them with the blazar \txs{}, raising a possible association of VHE neutrinos with this and other individual blazars. Very Long Baseline Interferometry (VLBI) is so far the only technique enabling the imaging of the innermost jet at milli-arcsec resolution (parsec scale), where the high energy emission possibly originates from. Here, we report on the radio properties of the parsec scale jet in \txs{} derived from the analysis of multi-epoch multi-frequency archive VLBI data. The half opening angle of the jet beam is about 3.8\degr, and the jet inclination angle is about 20\degr. The overall jet structure shows a helical trajectory with a precessing period of 5--6 years, likely originating from instabilities operating at parsec scales. The calculated beaming parameters (Doppler boosting factor, bulk Lorentz factor) suggest a moderately relativistic jet. The pc-scale magnetic field strength is estimated in the contexts of core-shift and variability, and is in general agreement in the range of 0.2 - 0.7 G. And it is found to decrease from a relatively larger value during the quiescent period before the ongoing flare. This suggests a conversion of magnetic field energy density to particle energy density that help accelerate injected particles at the jet base and result in variable shocked emission. The neutrino event could be associated with the onset of energetic particle injection into the jet. This scenario then supports the lepto-hadronic origin of the VHE neutrinos and $\gamma$-ray emission owing to a co-spatial origin.

Highly accreting quasars are quite luminous in the X-ray and optical regimes. While, they tend to become radio quiet and have optically thin radio spectra. Among the known quasars, IRAS F11119+3257 is a supercritical accretion source because it has a bolometric luminosity above the Eddington limit and extremely powerful X-ray outflows. To probe its radio structure, we investigated its radio spectrum between 0.15 and 96.15 GHz and performed very-long-baseline interferometric (VLBI) observations with the European VLBI Network (EVN) at 1.66 and 4.93 GHz. The deep EVN image at 1.66 GHz shows a two-sided jet with a projected separation about two hundred parsec and a very high flux density ratio of about 290. Together with the best-fit value of the integrated spectral index of -1.31+/-0.02 in the optically thin part, we infer that the approaching jet has an intrinsic speed at least 0.57 times of the light speed. This is a new record among the known all kinds of super-Eddington accreting sources and unlikely accelerated by the radiation pressure. We propose a scenario in which IRAS F11119+3257 is an unusual compact symmetric object with a small jet viewing angle and a radio spectrum peaking at 0.53+/-0.06 GHz mainly due to the synchrotron self-absorption.

We present the observational results from the 43-GHz Very Long Baseline Array (VLBA) observations of 124 compact radio-loud active galactic nuclei (AGNs) that were conducted between 2014 November and 2016 May. The typical dimensions of the restoring beam in each image are about 0.5 mas $\times$ 0.2 mas. The highest resolution of 0.2 mas corresponds to a physical size of 0.02 pc for the lowest redshift source in the sample. The 43-GHz very long baseline interferometry (VLBI) images of 97 AGNs are presented for the first time. We study the source compactness on milli-arcsec (mas) and sub-mas scales, and suggest that 95 sources in our sample are suitable for future space VLBI observations. By analyzing our data supplemented with other VLBA AGN surveys from literature, we find that the core brightness temperature increases with increasing frequency below a break frequency ~ 7 GHz, and decreases between ~7--240~GHz but increases again above~240 GHz in the rest frame of the sources. This indicates that the synchrotron opacity changes from optically thick to thin. We also find a strong statistical correlation between radio and $\gamma$-ray flux densities. Our correlation is tighter than those in literature derived from lower-frequency VLBI data, suggesting that the $\gamma$-ray emission is produced more co-spatially with the 43-GHz VLBA core emission. This correlation can also be extrapolated to the un-beamed AGN population, implying that a universal $\gamma$-ray production mechanism might be at work for all types of AGNs.

The fast-rising blue optical transient AT2018cow indicated unusual early phase characteristics unlike relatively better studied explosive transients. Its afterglow may be produced by either a relativistically beamed (jetted) or intrinsically luminous (non-jetted) ejecta and carries observational signatures of the progenitor and environment. High resolution monitoring can distinguish between these scenarios and clarify the progenitor nature. We present very long baseline interferometry (VLBI) observations of AT2018cow at 5 GHz involving 21 radio telescopes from the European VLBI Network with five sessions spanning ~ 1 year. With an astrometric precision up to 25 micro-arcseconds per epoch, the rapidly fading compact milli arcsecond scale source is found to be non-jetted with a proper motion of <= 0.15 mas/yr (0.14 c). This and a dense (number density ~ 10^4 - 10^5/cm^3) magnetized environment (magnetic field strength >= 0.84 G) are characteristic of a newly formed magnetar driven central engine, originating in the successful explosion of a low-mass star.

Blazars are a sub-class of quasars with Doppler boosted jets oriented close to the line of sight, and thus efficient probes of supermassive black hole growth and their environment, especially at high redshifts. Here we report on Very Long Baseline Interferometry observations of a blazar J0906+6930 at z = 5.47, which enabled the detection of polarised emission and measurement of jet proper motion at parsec scales. The observations suggest a less powerful jet compared with the general blazar population, including lower proper motion and bulk Lorentz factor. This coupled with a previously inferred high accretion rate indicate a transition from an accretion radiative power to a jet mechanical power based transfer of energy and momentum to the surrounding gas.While alternative scenarios could not be fully ruled out, our results indicate a possibly nascent jet embedded in and interacting with a dense medium resulting in a jet bending.

Space very long baseline interferometry (VLBI) has unique applications in high-resolution imaging of fine structure of astronomical objects and high-precision astrometry due to the key long space-Earth or space-space baselines beyond the Earth's diameter. China has been actively involved in the development of space VLBI in recent years. This review briefly summarizes China's research progress in space VLBI and the future development plan.

Apparently periodic optical variations of the luminous high-redshift (z=2.06) quasar PSO J334.2028+1.4075 led Liu et al. (2015) to interpret the variability as the orbital period of a binary supermassive black hole (SMBH) residing in a single circumbinary accretion disk. The proposed orbital separation was around 0.006 pc, and the possible inspiral time about 7 yr in the rest frame of the quasar. Such objects would be of high interest as the difficult-to-find end products of binary SMBH evolution, and potential sources of low-frequency gravitational waves. However, extending the time baseline of the variability study, Liu et al. (2016) later found that the periodicity of PSO J334.2028+1.4075 does not remain persistent. Foord et al. (2017) did not find evidence for the binary active galactic nucleus scenario based on Chandra X-ray observations. The object has also been studied in detail in the radio (Mooley et al. 2018) with the Karl G. Jansky Very Large Array (VLA) and the Very Long Baseline Array (VLBA), revealing a lobe-dominated quasar at kpc scales, and possibly a precessing jet, which might retain PSO J334.2028+1.4075 as a binary SMBH candidate. Here we report on our 1.7-GHz observation with the European VLBI Network (EVN) which complements the VLBA data taken at higher frequencies, and discuss the current knowledge about the nature of this interesting object.

A space-based very long baseline interferometry (VLBI) programme, named as the Cosmic Microscope, is proposed to involve dual VLBI telescopes in the space working together with giant ground-based telescopes (e.g., Square Kilometre Array, FAST, Arecibo) to image the low radio frequency Universe with the purpose of unraveling the compact structure of cosmic constituents including supermassive black holes and binaries, pulsars, astronomical masers and the underlying source, and exoplanets amongst others. The operational frequency bands are 30, 74, 330 and 1670 MHz, supporting broad science areas. The mission plans to launch two 30-m-diameter radio telescopes into 2,000 km x 90,000 km elliptical orbits. The two telescopes can work in flexibly diverse modes: (i) space-ground VLBI. The maximum space-ground baseline length is about 100,000 km; it provides a high-dynamic-range imaging capacity with unprecedented high resolutions at low frequencies (0.4 mas at 1.67 GHz and 20 mas at 30 MHz) enabling studies of exoplanets and supermassive black hole binaries (which emit nanoHz gravitational waves); (ii) space-space single-baseline VLBI. This unique baseline enables the detection of flaring hydroxyl masers, and more precise position measurement of pulsars and radio transients at milli-arcsecond level; (iii) single dish mode, where each telescope can be used to monitor transient bursts and rapidly trigger follow-up VLBI observations. The large space telescope will also contribute in measuring and constraining the total angular power spectrum from the Epoch of Reionization. In short, the Cosmic Microscope offers astronomers the opportunity to conduct novel, frontier science.

The search and study of close pairs of supermassive black holes (SMBHs) is important in the study of galaxy mergers which can possibly trigger active galactic nucleus (AGN) activity, and in the context of their evolution into the gravitational wave emitting regime. The quasar SDSS~J1048+0055 was identified as a SMBH binary candidate based on the observed double-peaked \OIII$\lambda\lambda$4959,5007 emission lines and two distinct radio components separated by $\sim 20$~pc \citep{2004ApJ...604L..33Z}. To ascertain the binary nature of this source, we analyzed multi-frequency, multi-epoch very long baseline interferometry (VLBI) data to investigate its pc-scale radio properties. The source shows double components with the western feature being brighter than the eastern one. This brighter component has a brightness temperature of $\sim 10^{10}$~K, spectral index of $\alpha = -0.09 \pm 0.09$ (flat) and is indicative of mildly relativistic beaming. In contrast, the faint component has a lower brightness temperature of $\sim 10^{8-9}$~K and steep spectrum. These clues are consistent with a core--jet structure, moreover, the apparent separation speed between the two components is much higher than the expected orbital motion in a binary SMBH. Thus the present study excludes the association of the two VLBI components with the cores of a SMBH binary, although the SMBH binary possibility (e.g., a pair of radio-loud and radio-quiet AGNs) is not fully ruled out. In the single active galactic nucleus (AGN) scenario, the double-peaked optical emission lines can originate from the jet interacting with the narrow-line region as indicated by a change in the jet direction at $\sim$ 140 pc.

Active galactic nuclei (AGN) accreting at rates close to the Eddington limit can host radiatively driven mildly relativistic outflows. Some of these X-ray absorbing but powerful outflows may produce strong shocks resulting in a significant non-thermal emission. This outflow-driven radio emission may be detectable in the radio-quiet quasar PDS 456 since it has a bolometric luminosity reaching the Eddington limit and a relativistic wide-aperture X-ray outflow with a kinetic power high enough to quench the star formation in its host galaxy. To investigate this possibility, we performed very-long-baseline interferometric (VLBI) observations of the quasar with the European VLBI Network (EVN) at 5 GHz. The EVN image with the full resolution reveals two faint and diffuse radio components with a projected separation of about 20 pc and an average brightness temperature of around two million Kelvin. In relation to the optical sub-mas-accuracy position measured by the Gaia mission, the two components are very likely on opposite sides of an undetected radio core. The VLBI structure at the deca-pc scale can thus be either a young jet or a bidirectional radio-emitting outflow, launched in the vicinity of a strongly accreting central engine. Two diffuse components at the hecto-pc scale, likely the relic radio emission from the past AGN activity, are tentatively detected on each side in the low-resolution EVN image.

Dual active galactic nuclei (DAGN) and supermassive black hole binaries (SMBHBs) at kpc and pc-scale separations, respectively, are expected during stages of galaxy merger and evolution. Their observational identification can address a range of areas of current astrophysics frontiers including the final parsec problem and their contribution towards the emission of low-frequency gravitational waves. This has however been difficult to achieve with current spectroscopy and time domain strategies. Very long baseline interferometry (VLBI) as a method of directly imaging radio structures with milli-arcsecond (mas) and sub-mas resolutions is introduced as a possible means of detecting DAGN and SMBHBs. We motivate its usage with expected observational signatures and cite some studies from literature to illustrate its current status, and present an updated list of candidates imaged with high-resolution radio observations. We then recall some shortcomings of the method with possible solutions and discuss future directions, relevant to large surveys with the upcoming Square Kilometer Array and future space VLBI missions.

The $\gamma$-ray bright blazar CTA 102 is studied using imaging (new 15 GHz and archival 43 GHz Very Long Baseline Array, VLBA data) and time variable optical flux density, polarization degree and electric vector position angle (EVPA) spanning between 2015 June 1 and 2016 October 1, covering a prominent $\gamma$-ray flare during 2016 January. The pc-scale jet indicates expansion with oscillatory features upto 17 mas. Component proper motions are in the range 0.04 - 0.33 mas/yr with acceleration upto 1.2 mas followed by a slowing down beyond 1.5 mas. A jet bulk Lorentz factor $\geq$ 17.5, position angle of 128.3 degrees, inclination angle $\leq$ 6.6 degrees and intrinsic half opening angle $\leq$ 1.8 degrees are derived from the VLBA data. These inferences are employed in a helical jet model to infer long term variability in flux density, polarization degree, EVPA and a rotation of the Stokes Q and U parameters. A core distance of $r_{\rm core,43 \ GHz}$ = 22.9 pc, and a magnetic field strength at 1 pc and the core location of 1.57 G and 0.07 G respectively are inferred using the core shift method. The study is useful in the context of estimating jet parameters and in offering clues to distinguish mechanisms responsible for variability over different timescales.

Radio-loud active galactic nuclei (AGNs), hosting powerful relativistic jet outflows, provide an excellent laboratory for studying jet physics. Very long baseline interferometry (VLBI) enables high-resolution imaging on milli-arcsecond (mas) and sub-mas scales, making it a powerful tool to explore the inner jet structure, shedding light on the formation, acceleration and collimation of AGN jets. In this paper, we present Very Long Baseline Array (VLBA) observations of ten radio-loud AGNs at 43 and 86~GHz, which were selected from the {\it Planck} catalogue of compact sources and are among the brightest in published VLBI images at and below 15 GHz. The image noise levels in our observations are typically 0.3 mJy beam$^{-1}$ and 1.5 mJy beam$^{-1}$ at 43 and 86 GHz, respectively. Compared with the VLBI data observed at lower frequencies from the literature, our observations with higher resolution (the highest resolution up to 0.07 mas at 86 GHz and 0.18 mas at 43 GHz) and at higher frequencies detected new jet components at sub-parsec scales, offering valuable data for studies of the physical properties of innermost jets. These include compactness factor of the radio structure (the ratio of core flux density to total flux density), and core brightness temperature ($T_{\rm b}$). In all these sources, the compact core accounts for a significant fraction ($> 60\%$) of the total flux density. Their correlated flux density at the longest baselines is higher than 0.16 Jy. The compactness of these sources make them good phase calibrators of mm-wavelength ground-based and space VLBI.

Radio astronomy observational facilities are under constant upgradation and development to achieve better capabilities including increasing the time and frequency resolutions of the recorded data, and increasing the receiving and recording bandwidth. As only a limited spectrum resource has been allocated to radio astronomy by the International Telecommunication Union, this results in the radio observational instrumentation being inevitably exposed to undesirable radio frequency interference (RFI) signals which originate mainly from terrestrial human activity and are becoming stronger with time. RFIs degrade the quality of astronomical data and even lead to data loss. The impact of RFIs on scientific outcome is becoming progressively difficult to manage. In this article, we motivate the requirement for RFI mitigation, and review the RFI characteristics, mitigation techniques and strategies. Mitigation strategies adopted at some representative observatories, telescopes and arrays are also introduced. We also discuss and present advantages and shortcomings of the four classes of RFI mitigation strategies, applicable at the connected causal stages: preventive, pre-detection, pre-correlation and post-correlation. The proper identification and flagging of RFI is key to the reduction of data loss and improvement in data quality, and is also the ultimate goal of developing RFI mitigation techniques. This can be achieved through a strategy involving a combination of the discussed techniques in stages. Recent advances in high speed digital signal processing and high performance computing allow for performing RFI excision of large data volumes generated from large telescopes or arrays in both real time and offline modes, aiding the proposed strategy.

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.

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 studied the pc-scale core shift effect using radio light curves for three blazars, S5 0716+714, 3C 279 and BL Lacertae, which were monitored at five frequencies ($\nu$) between 4.8 GHz and 36.8 GHz using the University of Michigan Radio Astronomical Observatory (UMRAO), the Crimean Astrophysical Observatory (CrAO), and Metsahovi Radio Observatory for over 40 years. Flares were Gaussian fitted to derive time delays between observed frequencies for each flare ($\Delta t$), peak amplitude ($A$), and their half width. Using $A \propto \nu^{\alpha}$ we infer $\alpha$ in the range $-$16.67 to 2.41 and using $\Delta t \propto \nu^{1/k_r}$, we infer $k_r \sim 1$, employed in the context of equipartition between magnetic and kinetic energy density for parameter estimation. From the estimated core position offset ($\Omega_{r \nu}$) and the core radius ($r_{\rm core}$), we infer that opacity model may not be valid in all cases. The mean magnetic field strength at 1 pc ($B_1$) and at the core ($B_{\rm core}$), are in agreement with previous estimates. We apply the magnetically arrested disk model to estimate black hole spins in the range $0.15-0.9$ for these blazars, indicating that the model is consistent with expected accretion mode in such sources. The power law shaped power spectral density has slopes $-$1.3 to $-$2.3 and is interpreted in terms of multiple shocks or magnetic instabilities.

Radio-loud high-redshift quasars (HRQs), although only a few of them are known to date, are crucial for the studies of the growth of supermassive black holes (SMBHs) and the evolution of active galactic nuclei (AGN) at early cosmological epochs. Radio jets offer direct evidence of SMBHs, and their radio structures can be studied with the highest angular resolution using Very Long Baseline Interferometry (VLBI). Here we report on the observations of three HRQs (J0131-0321, J0906+6930, J1026+2542) at z>5 using the Korean VLBI Network and VLBI Exploration of Radio Astrometry Arrays (together known as KaVA) with the purpose of studying their pc-scale jet properties. The observations were carried out at 22 and 43 GHz in 2016 January among the first-batch open-use experiments of KaVA. The quasar J0906+6930 was detected at 22 GHz but not at 43 GHz. The other two sources were not detected and upper limits to their compact radio emission are given. Archival VLBI imaging data and single-dish 15-GHz monitoring light curve of J0906+6930 were also acquired as complementary information. J0906+6930 shows a moderate-level variability at 15 GHz. The radio image is characterized by a core-jet structure with a total detectable size of ~5 pc in projection. The brightness temperature, 1.9x10^{11} K, indicates relativistic beaming of the jet. The radio properties of J0906+6930 are consistent with a blazar. Follow-up VLBI observations will be helpful for determining its structural variation.

The quasar 3C~286 is one of two compact steep spectrum sources detected by the {\it Fermi}/LAT. Here, we investigate the radio properties of the parsec(pc)-scale jet and its (possible) association with the $\gamma$-ray emission in 3C~286. The Very Long Baseline Interferometry (VLBI) images at various frequencies reveal a one-sided core--jet structure extending to the southwest at a projected distance of $\sim$1 kpc. The component at the jet base showing an inverted spectrum is identified as the core, with a mean brightness temperature of $2.8\times 10^{9}$~K. The jet bends at about 600 pc (in projection) away from the core, from a position angle of $-135^\circ$ to $-115^\circ$. Based on the available VLBI data, we inferred the proper motion speed of the inner jet as $0.013 \pm 0.011$ mas yr$^{-1}$ ($\beta_{\rm app} = 0.6 \pm 0.5$), corresponding to a jet speed of about $0.5\,c$ at an inclination angle of $48^\circ$ between the jet and the line of sight of the observer. The brightness temperature, jet speed and Lorentz factor are much lower than those of $\gamma$-ray-emitting blazars, implying that the pc-scale jet in 3C~286 is mildly relativistic. Unlike blazars in which $\gamma$-ray emission is in general thought to originate from the beamed innermost jet, the location and mechanism of $\gamma$-ray emission in 3C~286 may be different as indicated by the current radio data. Multi-band spectrum fitting may offer a complementary diagnostic clue of the $\gamma$-ray production mechanism in this source.

The quasar PG 1302$-$102 is believed to harbour a supermassive binary black hole (SMBBH) system. Using the available 15 GHz and $2-8$ GHz, multi-epoch Very Long Baseline Array data, we constrain the pc-scale jet properties based on the inferred mean proper motion, including a bulk Lorentz factor $\geq 5.1 \pm 0.8$, jet inclination angle $\leq (11.4 \pm 1.7)$ degrees, projected position angle $= 31.8$ degrees, intrinsic half opening angle $\leq (0.9 \pm 0.1)$ degrees and a mean $2-8$ GHz spectral index of 0.31. A general relativistic helical jet model is presented and applied to predict quasi-periodic oscillations of $\sim$ 10 days, power law power spectrum shape and a contribution of up to $\sim$ 53 percent to the observed variable core flux density. The model is used to make a case for high resolution, moderately sampled, long duration radio interferometric observations to reveal signatures due to helical knots and distinguish them from those due to SMBBH orbital activity including a phase difference $\sim \pi$ and an amplitude ratio (helical light curve amplitude/SMBBH light curve amplitude) of $0.2-3.3$. The prescription can be used to identify helical kinematic signatures from quasars, providing possible candidates for further studies with polarization measurements. It can also be used to infer promising SMBBH candidates for the study of gravitational waves if there are systematic deviations from helical signatures.

PKS 0202$+$149 is a low-power radio source with blazar-like $\gamma$-ray AGN characteristics and its properties and classification are investigated in relation to its $\gamma$-ray characteristics. It shows a hint of low frequency turnover at about 200 MHz. Radio imaging data of 0202$+$149 at different frequencies show differing morphologies on both kilo-parsec (kpc) and parsec (pc) scales. The overall source shows a triple structure of core and double lobes with a total projected size of $\sim$1.3 kpc. The compact source structure of 0202+149 is reminiscent of a compact steep spectrum (CSS) source. At pc scales a core-jet structure extends $\sim$25 pc (in projection) at a position angle perpendicular to the kpc-scale structure. The curved pc-scale structure with jet and inner lobe suggests that the CSS nuclear activity has recently re-started although its power has been decreasing, while the kpc-scale lobes are relics of earlier activity. A maximum apparent superluminal motion of $\sim16\,c$ is detected in the jet components, indicating a highly relativistic jet flow. The brightness temperature of the core is lower than the average value found for highly-beamed $\gamma$-ray AGN, indicating a lower radio power and a relatively lower Doppler boosting factor. The CSS radio classification indicates that blazar-like $\gamma$-ray properties can also be manifested in low-power CSS radio sources with appropriate jet and beaming properties.

Increasingly there is a need to develop astronomical visualisation and manipulations tools which allow viewers to interact with displayed data directly, in real time and across a range of platforms. In addition, increases in dynamic range available for astronomical images with next generation telescopes have led to a desire to develop enhanced visualisations capable of presenting information across a wide range of intensities. This paper describes three new tools for astronomical visualisation and image manipulation that are the result of a collaboration between software engineers and radio astronomers. The first tool, FITS3D, is a fast, interactive 3D data cube viewer designed to allow real-time interactive comparisons of multiple spectral line data cubes simultaneously. It features region specific selection manipulation including smoothing. The second tool, FITS2D, aids the visualisation and manipulation of 2D fits images. The tool supports the interactive creation of free-form masks which allow the user to extract any (potentially non-contiguous) subset of a fits image. It also supports annotations which can be placed without affecting the underlying data. The final tool is an R package for applying high dynamic range compression techniques to 2D fits images. This allows the full range of pixel brightness to be imaged in a single image, simultaneously showing the detail in bright sources while preserving the distinction of faint sources. Here we will present these three tools and demonstrate their capability using images from a range of astronomical images.

We study the long term Kepler light curve of the blazar W2R 1926+42 ($\sim$ 1.6 years) which indicates a variety of variability properties during different intervals of observation. The normalized excess variance, $F_{\rm var}$ ranges from 1.8 % in the quiescent phase and 43.3 % in the outburst phase. We find no significant deviation from linearity in the $F_{\rm var}$-flux relation. Time series analysis is conducted using the Fourier power spectrum and the wavelet analysis methods to study the power spectral density (PSD) shape, infer characteristic timescales and statistically significant quasi-periodic oscillations (QPOs). A bending power law with an associated timescale of $T_B = 6.2^{+6.4}_{-3.1}$ hours is inferred in the PSD analysis. We obtain a black hole mass of $M_\bullet = (1.5 - 5.9) \times 10^7 M_\odot$ for the first time using $F_{\rm var}$ and the bend timescale for this source. From a mean outburst lifetime of days, we infer a distance from the jet base $r \leq 1.75$ pc indicating that the outburst originates due to a shock. A possible QPO peaked at 9.1 days and lasting 3.4 cycles is inferred from the wavelet analysis. Assuming that the QPO is a true feature, $r = (152 - 378)~ G M_\bullet/c^2$ and supported by the other timing analysis products such as a weighted mean PSD slope of $-1.5 \pm 0.2$ from the PSD analysis, we argue that the observed variability and the weak and short duration QPO could be due to jet based processes including orbital features in a relativistic helical jet and others such as shocks and turbulence.

We study the core shift effect in the parsec scale jet of the blazar 3C 454.3 using the 4.8 GHz - 36.8 GHz radio light curves obtained from three decades of continuous monitoring. From a piecewise Gaussian fit to each flare, time lags $\Delta t$ between the observation frequencies $\nu$ and spectral indices $\alpha$ based on peak amplitudes $A$ are determined. From the fit $\Delta t \propto \nu^{1/k_r}$, $k_r = 1.10 \pm 0.18$ indicating equipartition between the magnetic field energy density and the particle energy density. From the fit $A \propto \nu^\alpha$, $\alpha$ is in the range $-0.24$ to $1.52$. A mean magnetic field strength at 1 pc, $B_1 = 0.5 \pm 0.2$ G, and at the core, $B_{\rm core} = 46 \pm 16$ mG, are inferred, consistent with previous estimates. The measure of core position offset is $\Omega_{r\nu} = 6.4 \pm 2.8$ pc GHz$^{1/k_r}$ when averaged over all frequency pairs. Based on the statistical trend shown by the measured core radius $r_{\rm core}$ as a function of $\nu$, we infer that the synchrotron opacity model may not be valid for all cases. A Fourier periodogram analysis yields power law slopes in the range $-1.6$ to $-3.5$ describing the power spectral density shape and gives bend timescales in the range $0.52 - 0.66~$yr. This result, and both positive and negative $\alpha$, indicate that the flares originate from multiple shocks in a small region. Important objectives met in our study include: the demonstration of the computational efficiency and statistical basis of the piecewise Gaussian fit; consistency with previously reported results; evidence for the core shift dependence on observation frequency and its utility in jet diagnostics in the region close to the resolving limit of very long baseline interferometry observations.

We present a general relativistic (GR) model of jet variability in active galactic nuclei due to orbiting blobs in helical motion along a funnel or cone shaped magnetic surface anchored to the accretion disk near the black hole. Considering a radiation pressure driven flow in the inner region, we find that it stabilizes the flow, yielding Lorentz factors ranging between 1.1 and 7 at small radii for reasonable initial conditions. Assuming these as inputs, simulated light curves (LCs) for the funnel model include Doppler and gravitational shifts, aberration, light bending, and time delay. These LCs are studied for quasi-periodic oscillations (QPOs) and the power spectral density (PSD) shape and yield an increased amplitude ($\sim$ 12 %); a beamed portion and a systematic phase shift with respect to that from a previous special relativistic model. The results strongly justify implementing a realistic magnetic surface geometry in Schwarzschild geometry to describe effects on emission from orbital features in the jet close to the horizon radius. A power law shaped PSD with a typical slope of $-2$ and QPOs with timescales in the range of $(1.37 - 130.7)$ days consistent with optical variability in Blazars, emerges from the simulations for black hole masses $M_{\bullet} = (0.5 - 5) \times 10^8 M_{\odot}$ and initial Lorentz factors $\gamma_{jet,i} = 2 - 10$. The models presented here can be applied to explain radio, optical, and X-ray variability from a range of jetted sources including active galactic nuclei, X-ray binaries and neutron stars.

We present theoretical models of X-ray variability attributable to orbital signatures from an accretion disk including emission region size, quasi-periodic oscillations (QPOs) and its quality factor $Q$, and the emergence of a break frequency in the power spectral density shape. We find a fractional variability amplitude of $F_{var}\propto M^{-0.4}_{\bullet}$. We conduct a time series analysis on X-ray light curves ($0.3-10$ keV) of a sample of active galactic nuclei (AGNs). A statistically significant bend frequency is inferred in 9 of 58 light curves (16 %) from 3 AGNs for which the break timescale is consistent with the reported BH spin but not with the reported BH mass. Upper limits of $2.85 \times 10^7 M_\odot$ in NGC 4051, $8.02 \times 10^7 M_\odot$ in MRK 766 and $4.68 \times 10^7 M_\odot$ in MCG-6-30-15 are inferred for maximally spinning BHs. For REJ 1034+396, where a QPO at 3733 s was reported, we obtain an emission region size of $(6 - 6.5) M$ and a BH spin $a\lesssim$ 0.08. The relativistic inner region of a thin disk, dominated by radiation pressure and electron scattering is likely to host the orbital features as the simulated $Q$ ranges from $6.3 \times 10^{-2}$ to $4.25 \times 10^6$, containing the observed $Q$. The derived value of $Q \sim$ 32 for REJ 1034+396 therefore suggests that the AGN hosts a thin disk.

Variability in active galactic nuclei is observed in ultraviolet to X-ray emission based light curves. This could be attributed to orbital signatures of the plasma that constitutes the accretion flow on the putative disk or in the developing jet close to the inner region of the central black hole. We discuss some theoretical models which build on this view. These models include general relativistic effects such as light bending, aberration effects, gravitational and Doppler redshifts. The novel aspects relate to the treatment of helical flow in cylindrical and conical geometries in the vicinity of a Schwarzschild black hole that leads to amplitude and frequency modulations of simulated light curves as well as the inclusion of beaming effects in these idealized geometries. We then present a suite of time series analysis techniques applicable to data with varied properties which can extract detailed information from them for their use in theoretical models.

The maximum likelihood estimator is used to determine fit parameters for various parametric models of the Fourier periodogram followed by the selection of the best fit model amongst competing models using the Akaike information criteria. This analysis, when applied to light curves of active galactic nuclei can be used to infer the presence of quasi-periodicity and break or knee frequencies. The extracted information can be used to place constraints on the mass, spin and other properties of the putative central black hole and the region surrounding it through theoretical models involving disk and jet physics.

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.