Chaos in a Magnetized Brane-World Spacetime Using Explicit Symplectic Integrators

A brane-world metric with an external magnetic field is a modified theory of gravity. It is suitable for the description of compact sources on the brane such as stars and black holes. We design a class of explicit symplectic integrators for this spacetime and use one of the integrators to investigate how variations of the parameters affect the motion of test particles. When the magnetic field does not vanish, the integrability of the system is destroyed. Thus, the onset of chaos can be allowed under some circumstances. Chaos easily occurs when the electromagnetic parameter becomes large enough. Dark matter acts as a gravitational force, so that chaotic motion can become more obvious as dark matter increases. The gravity of the black hole is weakened with an increasing positive cosmological parameter; therefore, the extent of chaos can be also strengthened. The proposed symplectic integrator is applied to a ray-tracing method and the study of such chaotic dynamics will be a possible reference for future studies of brane-world black hole shadows with chaotic patterns of self-similar fractal structures based on the Event Horizon Telescope data for M87* and Sagittarius A*.

Charged Particle Motions near Non-Schwarzschild Black Holes with External Magnetic Fields in Modified Theories of Gravity

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.

Dynamics of Charged Particles Moving around Kerr Black Hole with Inductive Charge and External Magnetic Field

We mainly focus on the effects of small changes of parameters on the dynamics of charged particles around Kerr black holes surrounded by an external magnetic field, which can be considered as a tidal environment. The radial motions of charged particles on the equatorial plane are studied via an effective potential. It is found that the particle energies at the local maxima values of the effective potentials increase with an increase in the black hole spin and the particle angular momenta, but decrease with an increase of one of the inductive charge parameter and magnetic field parameter. The radii of stable circular orbits on the equatorial plane also increase, whereas those of the innermost stable circular orbits decrease. On the other hand, the effects of small variations of the parameters on the orbital regular and chaotic dynamics of charged particles on the non-equatorial plane are traced by means of a time-transformed explicit symplectic integrator, Poincaré sections and fast Lyapunov indicators. It is shown that the dynamics sensitivity depends on small variations in the inductive charge parameter, magnetic field parameter, energy, and angular momentum. Chaos occurs easily as each of the inductive charge parameter, magnetic field parameter, and energy increases but is weakened as the angular momentum increases. When the dragging effects of the spacetime increase, the chaotic properties are not always weakened under some circumstances.

Adjustment of Force–Gradient Operator in Symplectic Methods

Many force–gradient explicit symplectic integration algorithms have been designed for the Hamiltonian H=T(p)+V(q) with kinetic energy T(p)=p2/2 in the existing references. When a force–gradient operator is appropriately adjusted as a new operator, it is still suitable for a class of Hamiltonian problems H=K(p,q)+V(q) with integrable part K(p,q)=∑i=1n∑j=1naijpipj+∑i=1nbipi, where aij=aij(q) and bi=bi(q) are functions of coordinates q. The newly adjusted operator is not a force–gradient operator but is similar to the momentum-version operator associated to the potential V. The newly extended (or adjusted) algorithms are no longer solvers of the original Hamiltonian, but are solvers of slightly modified Hamiltonians. They are explicit symplectic integrators with symmetry or time reversibility. Numerical tests show that the standard symplectic integrators without the new operator are generally poorer than the corresponding extended methods with the new operator in computational accuracies and efficiencies. The optimized methods have better accuracies than the corresponding non-optimized counterparts. Among the tested symplectic methods, the two extended optimized seven-stage fourth-order methods of Omelyan, Mryglod and Folk exhibit the best numerical performance. As a result, one of the two optimized algorithms is used to study the orbital dynamical features of a modified Hénon–Heiles system and a spring pendulum. These extended integrators allow for integrations in Hamiltonian problems, such as the spiral structure in self-consistent models of rotating galaxies and the spiral arms in galaxies.