Construction of second-order six-dimensional Hamiltonian-conserving scheme

Abstract

It is shown analytically that the energy-conserving implicit nonsymplectic scheme of Bacchini, Ripperda, Chen and Sironi provides a first-order accuracy to numerical solutions of a six-dimensional conservative Hamiltonian system. Because of this, a new second-order energy-conserving implicit scheme is proposed. Numerical simulations of Galactic model hosting a BL Lacertae object and magnetized rotating black hole background support these analytical results. The new method with appropriate time steps is used to explore the effects of varying the parameters on the presence of chaos in the two physical models. Chaos easily occurs in the Galactic model as the mass of the nucleus, the internal perturbation parameter, and the anisotropy of the potential of the elliptical galaxy increase. The dynamics of charged particles around the magnetized Kerr spacetime is easily chaotic for larger energies of the particles, smaller initial angular momenta of the particles, and stronger magnetic fields. The chaotic properties are not necessarily weakened when the black hole spin increases. The new method can be used for any six-dimensional Hamiltonian problems, including globally hyperbolic spacetimes with readily available (3+1) split coordinates.

26 pages,13 figures

References in corpus (18)

• Lyapunov indices with two nearby trajectories in a curved spacetime
• Transition from Regular to Chaotic Circulation in Magnetized Coronae near Compact Objects
• Symplectic structure of post-Newtonian Hamiltonian for spinning compact binaries
• Construction of Explicit Symplectic Integrators in General Relativity. I. Schwarzschild Black Holes
• On post-Newtonian orbits and the Galactic-center stars
• Construction of explicit symplectic integrators in general relativity. II. Reissner-Nordstrom black holes
• Construction of explicit symplectic integrators in general relativity. IV. Kerr black holes
• Inducing chaos by breaking axial symmetry in a black hole magnetosphere
• Construction of explicit symplectic integrators in general relativity. III. Reissner-Nordstrom-(anti)-de Sitter black holes
• Symplectic Integration of Post-Newtonian Equations of Motion with Spin
• A Symmetric Integrator for non-integrable Hamiltonian Relativistic Systems
• Applying explicit symplectic integrator to study chaos of charged particles around magnetized Kerr black hole
• A new dynamical model for the study of galactic structure
• Extended phase-space symplectic-like integrators for coherent post-Newtonian Euler-Lagrange equations
• A novel energy-conserving scheme for eight-dimensional systems
• Energy-conserving integrator for conservative Hamiltonian systems with ten-dimensional phase space
• FANTASY: User-Friendly Symplectic Geodesic Integrator for Arbitrary Metrics with Automatic Differentiation
• Determining the nature of orbits in a three-dimensional galaxy model hosting a BL Lacertae object