We study behaviour of ionized region of a Keplerian disk orbiting a Schwarzschild black hole immersed in an asymptotically uniform magnetic field. In dependence on the magnetic parameter B, and inclination angle θ of the disk plane with respect to the magnetic field direction, the charged particles of the ionized disk can enter three regimes: a) regular oscillatory motion, b) destruction due to capture by the magnetized black hole, c) chaotic regime of the motion. In order to study transition between the regular and chaotic type of the charged particle motion, we generate time series of the solution of equations of motion under various conditions, and study them by non-linear (box counting, correlation dimension, Lyapunov exponent, recurrence analysis, machine learning) methods of chaos determination. We demonstrate that the machine learning method appears to be the most efficient in determining the chaotic region of the θ − r space. We show that the chaotic character of the ionized particle motion increases with the inclination angle. For the inclination angles θ ∼ 0 whole the ionized internal part of the Keplerian disk is captured by the black hole.

On 2019/07/30.86853 UT, IceCube detected a high-energy astrophysical neutrino candidate. The Flat Spectrum Radio Quasar PKS 1502+106 is located within the 50 per cent uncertainty region of the event. Our analysis of 15 GHz Very Long Baseline Array (VLBA) and astrometric 8 GHz VLBA data, in a time span prior and after the IceCube event, reveals evidence for a radio ring structure which develops with time. Several arc-structures evolve perpendicular to the jet ridge line. We find evidence for precession of a curved jet based on kinematic modelling and a periodicity analysis. An outflowing broad line region (BLR) based on the C IV line emission (Sloan Digital Sky Survey, SDSS) is found. We attribute the atypical ring to an interaction of the precessing jet with the outflowing material. We discuss our findings in the context of a spine-sheath scenario where the ring reveals the sheath and its interaction with the surroundings (narrow line region, NLR, clouds). We find that the radio emission is correlated with the γ-ray emission, with radio lagging the γ-rays. Based on the γ-ray variability timescale, we constrain the γ-ray emission zone to the BLR (30-200 rg) and within the jet launching region. We discuss that the outflowing BLR provides the external radiation field for γ-ray production via external Compton scattering. The neutrino is most likely produced by proton-proton interaction in the blazar zone (beyond the BLR), enabled by episodic encounters of the jet with dense clouds, i.e. some molecular cloud in the NLR.

The X-ray emission from blazars has been widely investigated using several space telescopes. In this work, we explored statistical properties of the X-ray variability in the blazars S5 0716+714, OJ 287, Mkn 501 and RBS 2070 using the archival observations from the XMM-Newton telescope between the period 2002–2020. Several methods of timing and spectral analyses including fractional variability, minimum variability timescale, power spectral density analyses and countrate distribution were performed. In addition, we fitted various spectral models to the observations as well as estimated hardness ratio. The results show that the sources are moderately variable within the intra-day timescale. Three of the four sources exhibited a clear bi-modal pattern in their countrate distribution revealing possible indication of two distinct countrate states, that is, hard and soft countrate states. The slope indices of the power spectral density were found to be centered around 0.5. Furthermore, the spectra of the sources were fitted with single power-law, broken power-law, log-parabolic and black-body+log-parabolic models (the latter only for OJ 287). We conclude that for most of the observations log-parabolic model was the best fit. The power spectral density analysis revealed the variable nature of PSD slopes in the source light curves. The results of this analysis could indicate the non-stationary nature of the blazar processes on intra-day timescales. The observed features can be explained within the context of current blazar models, in which the non-thermal emission mostly arises from kilo-parsec scale relativistic jets.