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.

We study charged-fluid toroidal structures surrounding a non-rotating charged black hole immersed in a large-scale, asymptotically uniform magnetic field. In continuation of our former study on electrically charged matter in approximation of zero conductivity, we demonstrate existence of orbiting structures in permanent rigid rotation in the equatorial plane, and charged clouds hovering near the symmetry axis. We constrain the range of parameters that allow stable configurations and derive the geometrical shape of equi-pressure surfaces. Our simplified analytical study suggests that these regions of stability may be relevant for trapping electrically charged particles and dust grains in some areas of the black hole magnetosphere, being thus important in some astrophysical situations.

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.

The possible occurrence of equilibrium off-equatorial tori in the gravitational and electromagnetic fields of astrophysical compact objects has been recently proved based on non-ideal MHD theory. These stationary structures can represent plausible candidates for the modelling of coronal plasmas expected to arise in association with accretion discs. However, accretion disc coronae are formed by a highly diluted environment, and so the fluid description may be inappropriate. The question is posed of whether similar off-equatorial solutions can be determined also in the case of collisionless plasmas for which treatment based on kinetic theory, rather than fluid one, is demanded. In this paper the issue is addressed in the framework of the Vlasov-Maxwell description for non-relativistic multi-species axisymmetric plasmas subject to an external dominant spherical gravitational and dipolar magnetic field. Equilibrium configurations are investigated and explicit solutions for the species kinetic distribution function are constructed, which are expressed in terms of generalized Maxwellian functions characterized by isotropic temperature and non-uniform fluid fields. The conditions for the existence of off-equatorial tori are investigated. It is proved that these levitating systems are admitted under general conditions when both gravitational and magnetic fields contribute to shaping the spatial profiles of equilibrium plasma fluid fields. Then, specifically-kinetic effects carried by the equilibrium solution are explicitly provided and identified here with diamagnetic, energy-correction and electrostatic contributions. It is shown that these kinetic terms characterize the plasma equation of state by introducing non-vanishing deviations from the assumption of thermal pressure.

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.

Astrophysical fluids may acquire non-zero electrical charge because of strong irradiation or charge separation in a magnetic field. In this case, electromagnetic and gravitational forces may act together and produce new equilibrium configurations, which are different from the uncharged ones. Following our previous studies of charged test particles and uncharged perfect fluid tori encircling compact objects, we introduce here a simple test model of a charged perfect fluid torus in strong gravitational and electromagnetic fields. In contrast to ideal magnetohydrodynamic models, we consider here the opposite limit of negligible conductivity, where the charges are tied completely to the moving matter. This is an extreme limiting case which can provide a useful reference against which to compare subsequent more complicated astrophysically-motivated calculations. To clearly demonstrate the features of our model, we construct three-dimensional axisymmetric charged toroidal configurations around Reissner-Nordstr\"om black holes and compare them with equivalent configurations of electrically neutral tori.

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.