Despite the fact that jets from black holes were first understood to exist over 40 years ago, we are still in ignorance about many primary aspects of these systems – including the radiation mechanism at high energies, the particle makeup of the jets, and how particles are accelerated, possibly to energies as high as 100 TeV and hundreds of kpc from the central engine. We focus in particular on the discovery (and mystery) of strong X-ray emission from radio jets on kpc-scales, enabled by the unequaled high resolution of the Chandra X-ray observatory. We review the main evidence for and against the viable models to explain this X-ray emission over the last 20 years. Finally, we present results of a recent study on the X-ray variability of kpc-scale jets, where we find evidence that between 30-100% of the X-ray jet population is variable at the tens-of-percent level. The short (∼years) variability timescale is incompatible with the IC/CMB model for the X-rays and implies extremely small structures embedded within the kpc-scale jet, and thus requires a reconsideration of many assumptions about jet structure and dynamics.

We present an observational multiwavelength campaign during 2018–19 for PBC J2333.9–2343, a giant radio galaxy with a bright central core associated to a blazar nucleus, whose structure could be due to a significant jet reorientation. We report flux increases by a factor of two or more on timescales shorter than a month, resembling flaring events. The cross correlation between the NIR and optical bands shows quasi-simultaneous variations arising from the jet. The optical variability properties of PBC J2333.9–2343 are more comparable to a sample of blazars than to non-blazar AGN. The SED of the nucleus shows two peaks, with a derived jet angle of 3 degrees, also typical of a blazar. Therefore, we confirm the presence of a blazar-like core in the center of this galaxy.

Thanks to the data of the WISE all-sky survey we discovered that the non-thermal infrared emission of blazars, the largest population of associated γ-ray sources, has peculiar spectral properties. Here we provide a summary of all results achieved on the infrared–γ-ray connection. We also show results on the latest statistical analysis of the tight correlation between the mid-infrared colors and the γ-ray spectral index for Fermi blazars, a connection that links both emitted powers and spectral shapes of particles accelerated in blazar jets over ten decades in frequency of the electromagnetic spectrum. Finally, we outline all developments performed in the last decade achieved using the infrared– γ-ray connection to discover hundreds of new blazars within the sample of unidentified γ-ray sources thanks to optical spectroscopic observations.

We present the results of the gamma-ray flux distribution analysis on 145 gamma-ray bright blazars observed by Fermi-LAT. For the gamma-ray flux distribution, we applied a log-normal distribution to discuss the nature of the high-energy emission processes of blazars and a power-law distribution convolved with a Poisson distribution to investigate the implications of gamma-ray bright blazars for neutrino emission. Both distributions can represent the observed flux distributions as well. The leptonic models, which give the physical relationship between neutrinos and gamma rays, indicate that the flaring contribution to the neutrino emission can be dominant for the power-law index less than ∼2.5. From the power-law distribution analysis, we found that the power-law index < 2.5 accounts for the 82 % blazars. This result suggests that the flaring contribution of blazars is dominant for high-energy neutrino emission.

Studying blazar radio variability on timescales ranging from months to years provides information on the sub-parsec-scale structures of the jets, and the physics of the central active galactic nuclei. In this study, we focus on the radio variability of 1158 blazars observed at 15 GHz through the Owens Valley Radio Observatory Blazar Monitoring Program, where these sources have been observed about twice a week for over a decade. We investigate the dependence of the variability amplitudes and timescales, derived using a simple model fit to the structure function, on the milliarcsecond radio core sizes measured by Very Long Baseline Interferometry. The most compact sources exhibit larger variability amplitudes and shorter variability timescales than the more extended sources; this could be explained by light travel-time effects.

In this work, we investigate the formation and early evolution phase of X-shaped radio galaxies using the Back-flow model. We show how the X-like winged morphology evolves over time in a tri-axial ambient medium, naturally. At this early stage of formation, we demonstrated that both the pair of jet lobes are actively pushing the ambient material out of their path of propagation, forming (X-ray) cavities that are surrounded by a shocked shell (X-ray bright rims) of swept materials. We also noticed how turbulent the wing is in comparison to the active lobe, generating sites of random shocks, indicating that the wings are not passively evolving structures. This study demonstrated that the ambiguous morphology observed in jets is also imprinted over the ambient medium, providing an alternative perspective in understanding the underlying physical process causing such ambiguities. Finally, we indicate that shearing instabilities cause mixing of ambient material at the shearing interface.

We analyse VLBI and optical images of AGNs and their host galaxies and look for statistical correlations between the shape and orientation of the galaxy and the direction of the jet. We utilise the Astrogeo catalogue, which has over 9000 VLBI sources, many of those with a clear core-jet like structure that allows for the jet position angle to be reliably determined. We then use the VLBI source positions to search for optical counterparts within various optical surveys. In order to parameterise the orientation and shape of the host galaxy, we fitted a Gaussian elliptical model to the optical image, taking the PSF into account. We check our own shape parameters from this fit against the ones provided by the optical surveys. As of yet, no clear correlation between the galaxy morphology and the jet direction is seen.

Minkowski’s Object and ‘Death Star galaxy’ are two of the famous cases of rare instances when a radio jet has been observed to directly hit a neighbouring galaxy. RAD12, the RAD@home citizen science discovery with GMRT being presented here, is not only a new system being added to nearly half a dozen rare cases known so far but also the first case where the neighbouring galaxy is not a minor/dwarf companion but a galaxy bigger than the host of the radio jet. Additionally, the jet appears to be one-sided and the jet after interaction completely stops and forms a bubble inflating laterally which is unlike previous cases of minor deviation or loss of collimation. Since the nature of radio jet-ISM coupling is poorly understood so far, more discovery of objects like RAD12 will be important to the understanding of galaxy evolution through merger and AGN feedback.

The RAD@home RGB-maker Tool is a python-based web application that enables citizen science research through collaboration using open-source technology. The tool fetches FITS image data from NASA SkyView and generates false colour images in Red-Green-Blue channels with contour. This helps in the basic multi-wavelength understanding and characterization of extragalactic objects, and further analysis along with reporting of potential new discoveries in a uniform format. Students from Universities or science-graduate citizens gain skills in RGB-C image analysis and identify unique features in objects via either one-day online-weekend events or week-long in-person training. Trained citizen scientists in turn are part of rare discoveries such as the jet-galaxy interaction system RAD12. The tool has been successfully used by thousands of citizen scientists in India since its launch on 26th January 2021 and has demonstrated the potential during IAU symposium 375 to be part of citizen science efforts with international participants.

Unlike most web-based citizen science programmes, RAD@home citizen scientists are trained in an interactive manner, both online and in person. This Collaboratory model can alleviate various geo-political and socio-economic constraints on the growth of people in underdeveloped regions which is in synergy with mandates of the IAU Office of Astronomy for Development. Owing to its Inter-University collaboratory design, supported by over 30 research institutions, it can attempt to convert the ‘Big Data Problem’ into a ‘Big Data Prospect’ by direct training and involving citizens in discoveries from multi-telescope data. Using the RAD-RGB-maker web-tool, e-astronomers learn the basics of multi-wavelength (UV-optical-IR-radio) extragalactic astronomy and create preliminary discovery reports of unusual radio sources thus becoming Co-Investigator in follow-up observations and publications (e.g. Discovery of RAD12). International expansion of RAD@home has been initiated during the recent IAU symposium 375 through its Educational and Outreach programme.

Particle-in-Cell simulations can provide a possible answer to an important key issue for astrophysical plasma jets; namely on how a toroidal magnetic field affects the evolution of pair and electron-ion jets. We show that Weibel, mushroom, and kinetic Kelvin-Helmhotz instabilities excited at the linear stage, generate a quasi-steady x-component of the electric field which accelerates and decelerates electrons. We observe significant differences in the structure of the strong electromagnetic fields that are driven by the kinetic instabilities with the pair jet. We find that the two different jet compositions (e± and e– - i+) generate different instability modes respectively. Moreover, the magnetic field in the non-linear stage generated by different instabilities is dissipated and reorganized into new topologies. A 3D magnetic field topology depiction indicates possible reconnection sites in the non-linear stage where the particles are significantly accelerated by the dissipation of the magnetic field associated to a possible reconnection manifestation.

Recent theoretical considerations and observational evidence evince the spine-sheath morphology of relativistic jets emitted from active galactic nuclei (AGNs) or gamma-ray bursts (GRBs). The resulting shear boundary layers (SBLs) are likely to be an avenue for particle acceleration in relativistic jets. The effect of radiation drag on radiating particles has yet to be addressed in most studies of particle acceleration at shear boundary layers, even though radiative cooling may considerably affect particle dynamics. By using particle-in-cell simulations, we study the effects of inverse Compton cooling on particle dynamics and emerging particle spectra.

We present a new algorithm for the identification and physical characterization of current sheets and reconnection sites as well as the update of post-reconnection particles spectra in 2D and 3D large scale relativistic magnetohydrodynamic simulations. Lagrangian particles, which follow the fluid, are used to sample plasma parameters before entering the reconnection sites that form during the evolution of the different configurations considered. With the sampled parameters and a subgrid model based on results of Particle-in-Cell simulations we introduced in the PLUTO code an algorithm able to describe the post-reconnection spectra associated to the non-thermal component.

With high-sensitivity kiloparsec-scale radio polarimetry, we can examine the jet-medium interactions and get a better understanding of the blazar divide in radio-loud (RL) AGN. We are analyzing the radio polarimetric observations with the EVLA and GMRT of 24 quasars and BL Lacs belonging to the Palomar-Green (PG) sample. The RL quasars show extensive polarisation structures in their cores, jets, lobes, and hotspots, whereas preliminary results suggest that BL Lacs exhibit polarisation primarily in their cores and inner jet regions. These findings imply that both intrinsic (central engine-related) and extrinsic (environment-related) variables are important in the formation of the blazar subclasses. The Fanaroff-Riley (FR) dichotomy can also be studied assuming RL unification and looking through the lens of blazars. Due to the radio-unbiased nature of the optically/UV-selected PG sample, we find a large fraction of the PG quasars are restarted, distorted (S- or X-shaped), or have a hybrid FR morphology.

We performed photo-polarimetry observations for the blazar BL Lacertae (BL Lac) in 2020 to 2021, in which BL Lac showed historical outburst. As a result, we obtained microvariability with a timescale about five minutes and wavelength dependence of polarization degree and angle. These results indicate multiple emission regions and turbulent magnetic field structure.

This report presents the results of optical polarimetric observations carried out with 6-m and 1-m telescopes at SAO RAS. The study of the blazar S5 0716+714 radiation showed the presence of a period of the variability of brightness and polarization vector variations on scales of ∼1.5 hours, constant on a long time scale; multi-colour monitoring of BL Lac polarization before, during and after the flare demonstrates the difference in the patterns of polarization vector variability depending on the wavelength. Several geometrical models and physical descriptions are discussed.

The Fermi γ-ray telescope has detected 6658 sources, with 1845 of them remaining unidentified. We show that polarimetry of γ-ray fields is a powerful asset in the hunt of active galactic nuclei (AGN) as potential optical counterparts for γ-ray sources. We have studied an unidentified Fermi field (3FGL J0221.2+2518) and found a previously-unknown highly-polarized extragalactic object as a potential optical counterpart within the 1-sigma error ellipse of the corresponding γ-ray source. Based on a collection of data, we find that it most probably is a composite object: a star-forming galaxy accompanied by AGN. PASIPHAE is a large polarimetric experiment which will measure the polarisation of sources away from the galactic plane. This will provide an excellent opportunity to study hundreds of unidentified γ-ray sources and unveil potential optical counterparts, using polarimetry.

Relativistic jets from supermassive black holes or stellar mass black holes are among the most powerful astrophysical phenomena. Magnetic field plays an important role in the jet launching and propagation, as well as particle acceleration and radiation. Polarimetry is the only way to observe the magnetic field evolution. The recent launch of the Imaging X-ray Polarimetry Explorer (IXPE) has opened up the X-ray polarization window, which has revealed very interesting phenomena for relativistic jets. However, the field of MeV gamma-ray polarimetry remains largely unexplored. This paper aims to summarize key scientific potentials for MeV polarimetry for blazars and gamma-ray bursts (GRBs) from recent theoretical modeling. These predictions, which are closely related to the cosmic ray acceleration, neutrino production, radiation mechanism, and the jet evolution, can be examined by future MeV polarimeters, such as the Compton Spectrometer and Imager (COSI), the LargE Area burst Polarimeter (LEAP), and the All-sky Medium-Energy Gamma-ray Observatory eXplorer (AMEGO-X).

We present and briefly discuss results of several studies of the source J2102+6015 with tentatively defined redshift z = 4.575 which demonstrates unusual properties in imaging and astrometric VLBI observations. Its properties might be considered as indications on the supermassive black hole binary which can be considered as a so far rare example of a high-redshift source of known electromagnetic and, possibly, predictable gravitational wave emissions.

In recent years, evidence has started piling up that some high-energy cosmic neutrinos can be associated with blazars in flaring states. On 2022 February 26, a new blazar-neutrino coincidence was reported: the track-like neutrino event IC220225A detected by IceCube is spatially coincident with the flat-spectrum radio quasar PKS 0215+015. Like previous associations, this source was found to be in a high optical and γ-ray state. Moreover, the source showed a bright radio outburst, which substantially increases the probability of a true physical association. We have performed six observations with the VLBA shortly after the neutrino event with a monthly cadence and are monitoring the source with the Effelsberg 100m-Telescope, and with the Australia Compact Telescope Array. Here, we present first results on the contemporary parsec-scale jet structure of PKS 0215+015 in total intensity and polarization to constrain possible physical processes leading to neutrino emission in blazars.