In recent publications, the construction of explicit symplectic integrators for Schwarzschild- and Kerr-type spacetimes is based on splitting and composition methods for numerical integrations of Hamiltonians or time-transformed Hamiltonians associated with these spacetimes. Such splittings are not unique but have various options. A Hamiltonian describing the motion of charged particles around the Schwarzschild black hole with an external magnetic field can be separated into three, four, and five explicitly integrable parts. It is shown through numerical tests of regular and chaotic orbits that the three-part splitting method is the best of the three Hamiltonian splitting methods in accuracy. In the three-part splitting, optimized fourth-order partitioned Runge–Kutta and Runge–Kutta–Nyström explicit symplectic integrators exhibit the best accuracies. In fact, they are several orders of magnitude better than the fourth-order Yoshida algorithms for appropriate time steps. The first two algorithms have a small additional computational cost compared with the latter ones. Optimized sixth-order partitioned Runge–Kutta and Runge–Kutta–Nyström explicit symplectic integrators have no dramatic advantages over the optimized fourth-order ones in accuracy during long-term integrations due to roundoff errors. The idea of finding the integrators with the best performance is also suitable for Hamiltonians or time-transformed Hamiltonians of other curved spacetimes including Kerr-type spacetimes. When the numbers of explicitly integrable splitting sub-Hamiltonians are as small as possible, such splitting Hamiltonian methods would bring better accuracies. In this case, the optimized fourth-order partitioned Runge–Kutta and Runge–Kutta–Nyström methods are worth recommending.

Based on the scalar–tensor–vector modified gravitational theory, a modified gravity Schwarzschild black hole solution has been given in the existing literature. Such a black hole spacetime is obtained through the inclusion of a modified gravity coupling parameter, which corresponds to the modified gravitational constant and the black hole charge. In this sense, the modified gravity parameter acts as not only an enhanced gravitational effect but also a gravitational repulsive force contribution to a test particle moving around the black hole. Because the modified Schwarzschild spacetime is static spherical symmetric, it is integrable. However, the spherical symmetry and the integrability are destroyed when the black hole is immersed in an external asymptotic uniform magnetic field and the particle is charged. Although the magnetized modified Schwarzschild spacetime is nonintegrable and inseparable, it allows for the application of explicit symplectic integrators when its Hamiltonian is split into five explicitly integrable parts. Taking one of the proposed explicit symplectic integrators and the techniques of Poincaré sections and fast Lyapunov indicators as numerical tools, we show that the charged particle can have chaotic motions under some circumstances. Chaos is strengthened with an increase of the modified gravity parameter from the global phase space structures. There are similar results when the magnetic field parameter and the particle energy increase. However, an increase of the particle angular momentum weakens the strength of chaos.

No abstract

A small deformation to the Schwarzschild metric controlled by four free parameters could be referred to as a nonspinning black hole solution in alternative theories of gravity. Since such a non-Schwarzschild metric can be changed into a Kerr-like black hole metric via a complex coordinate transformation, the recently proposed time-transformed, explicit symplectic integrators for the Kerr-type spacetimes are suitable for a Hamiltonian system describing the motion of charged particles around the non-Schwarzschild black hole surrounded with an external magnetic field. The obtained explicit symplectic methods are based on a time-transformed Hamiltonian split into seven parts, whose analytical solutions are explicit functions of new coordinate time. Numerical tests show that such explicit symplectic integrators for intermediate time steps perform well long-term when stabilizing Hamiltonian errors, regardless of regular or chaotic orbits. One of the explicit symplectic integrators with the techniques of Poincaré sections and fast Lyapunov indicators is applied to investigate the effects of the parameters, including the four free deformation parameters, on the orbital dynamical behavior. From the global phase-space structure, chaotic properties are typically strengthened under some circumstances, as the magnitude of the magnetic parameter or any one of the negative deformation parameters increases. However, they are weakened when the angular momentum or any one of the positive deformation parameters increases.

No abstract