Construction of second-order six-dimensional Hamiltonian-conserving scheme
arXiv:2108.09882 · doi:10.3847/1538-4365/ac1ff3
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.
arXiv comment26 pages,13 figures
References in corpus (29)
- Lyapunov indices with two nearby trajectories in a curved spacetime
- Transition from Regular to Chaotic Circulation in Magnetized Coronae near Compact Objects
- Chaotic motion of Charged Particles in an Electromagnetic Field Surrounding a Rotating Black Hole
- RAPTOR I: Time-dependent radiative transfer in arbitrary spacetimes
- Symplectic structure of post-Newtonian Hamiltonian for spinning compact binaries
- Construction of Explicit Symplectic Integrators in General Relativity. I. Schwarzschild Black Holes
- Construction of explicit symplectic integrators in general relativity. II. Reissner-Nordstrom black holes
- On post-Newtonian orbits and the Galactic-center stars
- 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
- Generalized, energy-conserving numerical simulations of particles in general relativity. I. Time-like and null geodesics
- RAPTOR II: Polarized radiative transfer in curved spacetime
- Generalized, energy-conserving numerical simulations of particles in general relativity. II. Test particles in electromagnetic fields and GRMHD
- The Midpoint Rule as a Variational--Symplectic Integrator. I. Hamiltonian Systems
- A Symmetric Integrator for non-integrable Hamiltonian Relativistic Systems
- Second post-Newtonian Lagrangian dynamics of spinning compact binaries
- 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
- Exploring the nature of orbits in a galactic model with a massive nucleus
- Extended phase-space symplectic-like integrators for coherent post-Newtonian Euler-Lagrange equations
- A novel energy-conserving scheme for eight-dimensional systems
- The use of Kepler solver in numerical integrations of quasi-Keplerian orbits
- Energy-conserving integrator for conservative Hamiltonian systems with ten-dimensional phase space
- Order and Chaos in a Three-dimensional Binary System of Interacting Galaxies
- Revealing the character of orbits in a binary system consisting of a primary galaxy and a satellite companion
- 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