Relativistic effects of compact objects onto electromagnetic fields in their vicinity are investigated using the test-field approximation. In particular, we study the possible emergence of magnetic null points which are astrophysically relevant for the processes of magnetic reconnection. While the magnetic reconnection occurs in the presence of plasma and may lead to violent mass ejection, we show here that strong gravitation of the supermassive black hole may actively support the process by suitably entangling the field lines even in the electro-vacuum description. In this contribution we further discuss the case of a dipole-type magnetic field of the neutron star on the plunging trajectory to the supermassive black hole. While we have previously shown that given model in principle admits the formation of magnetic null points, here we explore whether and where the null points appear for the astrophysically relevant values of the parameters.

We explore the structure of a dipole-type vacuum field of a slowly rotating magnetic star near the horizon of a supermassive black hole, where the structure of field lines becomes highly distorted by effects of strong gravity. Such a situation may occur near a neutron star in the final stages of a plunging trajectory into a galactic center. We solve Maxwell's equations in the Rindler approximation for the rotating conducting source of dipolar magnetic field arbitrarily inclined with respect to the axis of rotation. For the fixed inclination angle we calculate the field including the radiative terms while in the general case we discuss the electromagnetic field considering the near-field terms only. In the latter case we investigate the emergence of magnetic null points within the vacuum magnetosphere. Null points become highly relevant in the presence of astrophysical plasma where they are connected with processes of magnetic reconnection and mass ejection.

An interplay of magnetic fields and gravitation drives accretion and outflows near black holes. However, a specific mechanism is still a matter of debate; it is very likely that different processes dominate under various conditions. In particular, for the acceleration of particles and their collimation in jets, an ordered component of the magnetic field seems to be essential. Here we discuss the role of large-scale magnetic fields in transporting the charged particles and dust grains from the bound orbits in the equatorial plane of a rotating (Kerr) black hole and the resulting acceleration along trajectories escaping the system in a direction parallel to the symmetry axis (perpendicular to the accretion disk). We consider a specific scenario of destabilization of circular geodesics of initially neutral matter by charging (e.g., due to photoionization). Some particles may be set on escaping trajectories and attain relativistic velocity. The case of charged particles differs from charged dust grains by their charge-to-mass ratio, but the acceleration mechanism operates in a similar manner. It appears that the chaotic dynamics controls the outflow and supports the formation of near-horizon escape zones. We employ the technique of recurrence plots to characterize the onset of chaos in the outflowing medium. We investigate the system numerically and construct the basin-boundary plots, which show the location and the extent of the escape zones. The effects of black hole spin and magnetic field strength on the formation and location of escape zones are discussed, and the maximal escape velocity is computed.

We present a detailed comparison of several integration schemes applied to the dynamic system consisting of a charged particle on the Kerr background endowed with the axisymmetric electromagnetic test field. In particular, we compare the performance of the symplectic integrator with several non-symplectic routines and discuss under which circumstances we should choose the symplectic one and when we should switch to some other scheme. We are basically concerned with two crucial, yet opposing aspects - accuracy of the integration and CPU time consumption. The latter is generally less critical in our application while the highest possible accuracy is strongly demanded.

Dynamics of charged matter in the oblique black hole magnetosphere is investigated. In particular, we adopt a model consisting of a rotating black hole embedded in the external large-scale magnetic field that is inclined arbitrarily with respect to the rotation axis. Breaking the axial symmetry appears to have profound consequences regarding the dynamics of particles and it also poses some methodological difficulties. In this contribution we discuss the applicability of the method of effective potential for the non-axisymmetric model and show that it may only be applied in the appropriate reference frame.

Dynamics of charged particles in the vicinity of a rotating black hole embedded in the external large-scale magnetic field is numerically investigated. In particular, we consider a non-axisymmetric model in which the asymptotically uniform magnetic field is inclined with respect to the axis of rotation. We study the effect of inclination onto the prevailing dynamic regime of particle motion, i.e. we ask whether the inclined field allows regular trajectories or if instead, the deterministic chaos dominates the motion. In this contribution we further discuss the role of initial condition, particularly, the initial azimuthal angle. To characterize the measure of chaoticness we compute maximal Lyapunov exponents and employ the method of Recurrence Quantification Analysis.

Circular motion of particles, dust grains and fluids in the vicinity of compact objects has been investigated as a model for accretion of gaseous and dusty environment. Here we further discuss, within the framework of general relativity, figures of equilibrium of matter under the influence of combined gravitational and large-scale magnetic fields, assuming that the accreted material acquires a small electric charge due to interplay of plasma processes and photoionization. In particular, we employ an exact solution describing the massive magnetic dipole and we identify the regions of stable motion. We also investigate situations when the particle dynamics exhibits the onset of chaos. In order to characterize the measure of chaoticness we employ techniques of Poincar\'e surfaces of section and of recurrence plots.

While the motion of particles near a rotating, electrically-neutral (Kerr), and charged (Kerr--Newman) black hole is always strictly regular, a perturbation in the gravitational or the electromagnetic field generally leads to chaos. The transition from regular to chaotic dynamics is relatively gradual if the system preserves axial symmetry, whereas non-axisymmetry induces chaos more efficiently. Here we study the development of chaos in an oblique (electro-vacuum) magnetosphere of a magnetized black hole. Besides the strong gravity of the massive source represented by the Kerr metric we consider the presence of a weak, ordered, large-scale magnetic field. An axially-symmetric model consisting of a rotating black hole embedded in an aligned magnetic field is generalized by allowing an oblique direction of the field having a general inclination, with respect to the rotation axis of the system. The inclination of the field acts as an additional perturbation to the motion of charged particles as it breaks the axial symmetry of the system and cancels the related integral of motion. The axial component of angular momentum is no longer conserved and the resulting system thus has three degrees of freedom. Our primary concern within this contribution is to find out how sensitive the system of bound particles is to the inclination of the field. We employ the method of the maximal Lyapunov exponent to distinguish between regular and chaotic orbits and to quantify their chaoticity. We find that even a small misalignment induces chaotic motion.

Within the framework of Bonnor's exact solution describing a massive magnetic dipole, we study the motion of neutral and electrically charged test particles. In dependence on the Bonnor spacetime parameters, we determine regions enabling the existence of stable circular orbits confined to the equatorial plane and of those levitating above the equatorial plane. Constructing Poincar\'e surfaces of section and recurrence plots, we also investigate the dynamics of particles moving along general off-equatorial trajectories bound in effective potential wells forming around the stable circular orbits. We demonstrate that the motion in the Bonnor spacetime is not integrable. This extends previous investigations of generalized St\"ormer's problem into the realm of exact solutions of Einstein-Maxwell equations, where the gravitational and electromagnetic effects play a comparable role on the particle motion.

Cosmic black holes can act as agents of particle acceleration. We study properties of a system consisting of a rotating black hole immersed in a large-scale organized magnetic field. Electrically charged particles in the immediate neighborhood of the horizon are influenced by strong gravity acting together with magnetic and induced electric components. We relax several constraints which were often imposed in previous works: the magnetic field does not have to share a common symmetry axis with the spin of the black hole but they can be inclined with respect to each other, thus violating the axial symmetry. Also, the black hole does not have to remain at rest but it can instead perform fast translational motion together with rotation. We demonstrate that the generalization brings new effects. Starting from uniform electro-vacuum fields in the curved spacetime, we find separatrices and identify magnetic neutral points forming in certain circumstances. We suggest that these structures can represent signatures of magnetic reconnection triggered by frame-dragging effects in the ergosphere. We further investigate the motion of charged particles in these black hole magnetospheres. We concentrate on the transition from the regular motion to chaos, and in this context we explore the characteristics of chaos in relativity. For the first time, we apply recurrence plots as a suitable technique to quantify the degree of chaoticness near a black hole.

We study the motion of charged test particles around a Kerr black hole immersed in the asymptotically uniform magnetic field, concluding that off-equatorial stable orbits are allowed in this system. Being interested in dynamical properties of these astrophysically relevant orbits we employ rather novel approach based on the analysis of recurrences of the system to the vicinity of its previous states. We use recurrence plots (RPs) as a tool to visualize recurrences of the trajectory in the phase space. Construction of RPs is simple and straightforward regardless of the dimension of the phase space, which is a major advantage of this approach when compared to the "traditional" methods of the numerical analysis of dynamical systems (for instance the visual survey of Poincar\'{e} surfaces of section, evaluation of the Lyapunov spectra etc.). We show that RPs and their quantitative measures (obtained from recurrence quantification analysis -- RQA) are powerful tools to detect dynamical regime of motion (regular vs. chaotic) and precisely locate the transitions between these regimes.

Accretion onto black holes and compact stars brings material in a zone of strong gravitational and electromagnetic fields. We study dynamical properties of motion of electrically charged particles forming a highly diluted medium (a corona) in the regime of strong gravity and large-scale (ordered) magnetic field. We start our work from a system that allows regular motion, then we focus on the onset of chaos. To this end, we investigate the case of a rotating black hole immersed in a weak, asymptotically uniform magnetic field. We also consider a magnetic star, approximated by the Schwarzschild metric and a test magnetic field of a rotating dipole. These are two model examples of systems permitting energetically bound, off-equatorial motion of matter confined to the halo lobes that encircle the central body. Our approach allows us to address the question of whether the spin parameter of the black hole plays any major role in determining the degree of the chaoticness. To characterize the motion, we construct the Recurrence Plots (RP) and we compare them with Poincar\'e surfaces of section. We describe the Recurrence Plots in terms of the Recurrence Quantification Analysis (RQA), which allows us to identify the transition between different dynamical regimes. We demonstrate that this new technique is able to detect the chaos onset very efficiently, and to provide its quantitative measure. The chaos typically occurs when the conserved energy is raised to a sufficiently high level that allows the particles to traverse the equatorial plane. We find that the role of the black-hole spin in setting the chaos is more complicated than initially thought.