No Arabic abstract
In this paper, an implicit nonsymplectic exact energy-preserving integrator is specifically designed for a ten-dimensional phase-space conservative Hamiltonian system with five degrees of freedom. It is based on a suitable discretization-averaging of the Hamiltonian gradient, with a second-order accuracy to numerical solutions. A one-dimensional disordered discrete nonlinear Schr{o}dinger equation and a post-Newtonian Hamiltonian system of spinning compact binaries are taken as our two examples. We demonstrate numerically that the proposed algorithm exhibits good long-term performance in the preservation of energy, if roundoff errors are neglected. This result is independent of time steps, initial orbital eccentricities, and regular and chaotic orbital dynamical behavior. In particular, the application of appropriately large time steps to the new algorithm is helpful in reducing time-consuming and roundoff errors. This new method, combined with fast Lyapunov indicators, is well suited related to chaos in the two example problems. It is found that chaos in the former system is mainly responsible for one of the parameters. In the latter problem, a combination of small initial separations and high initial eccentricities can easily induce chaos.
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.
Recent advances show that neural networks embedded with physics-informed priors significantly outperform vanilla neural networks in learning and predicting the long term dynamics of complex physical systems from noisy data. Despite this success, there has only been a limited study on how to optimally combine physics priors to improve predictive performance. To tackle this problem we unpack and generalize recent innovations into individual inductive bias segments. As such, we are able to systematically investigate all possible combinations of inductive biases of which existing methods are a natural subset. Using this framework we introduce Variational Integrator Graph Networks - a novel method that unifies the strengths of existing approaches by combining an energy constraint, high-order symplectic variational integrators, and graph neural networks. We demonstrate, across an extensive ablation, that the proposed unifying framework outperforms existing methods, for data-efficient learning and in predictive accuracy, across both single and many-body problems studied in recent literature. We empirically show that the improvements arise because high order variational integrators combined with a potential energy constraint induce coupled learning of generalized position and momentum updates which can be formalized via the Partitioned Runge-Kutta method.
First passage time experiments were used to explore the effects of low amplitude noise as a source of accelerated phase space diffusion in two-dimensional Hamiltonian systems, and these effects were then compared with the effects of periodic driving. The objective was to quantify and understand the manner in which ``sticky chaotic orbits that, in the absence of perturbations, are confined near regular islands for very long times, can become ``unstuck much more quickly when subjected to even very weak perturbations. For both noise and periodic driving, the typical escape time scales logarithmically with the amplitude of the perturbation. For white noise, the details seem unimportant: Additive and multiplicative noise typically have very similar effects, and the presence or absence of a friction related to the noise by a Fluctuation-Dissipation Theorem is also largely irrelevant. Allowing for colored noise can significantly decrease the efficacy of the perturbation, but only when the autocorrelation time becomes so large that there is little power at frequencies comparable to the natural frequencies of the unperturbed orbit. Similarly, periodic driving is relatively inefficient when the driving frequency is not comparable to these natural frequencies. This suggests strongly that noise-induced extrinsic diffusion, like modulational diffusion associated with periodic driving, is a resonance phenomenon. The logarithmic dependence of the escape time on amplitude reflects the fact that the time required for perturbed and unperturbed orbits to diverge a given distance scales logarithmically in the amplitude of the perturbation.
The dynamics of spinning particles in curved space-time is discussed, emphasizing the hamiltonian formulation. Different choices of hamiltonians allow for the description of different gravitating systems. We give full results for the simplest case with minimal hamiltonian, constructing constants of motion including spin. The analysis is illustrated by the example of motion in Schwarzschild space-time. We also discuss a non-minimal extension of the hamiltonian giving rise to a gravitational equivalent of the Stern-Gerlach force. We show that this extension respects a large class of known constants of motion for the minimal case.
As recent work continues to demonstrate, the study of relativistic scattering processes leads to valuable insights and computational tools applicable to the relativistic bound-orbit two-body problem. This is particularly relevant in the post-Minkowskian approach to the gravitational two-body problem, where the field has only recently reached a full description of certain physical observables for scattering orbits, including radiative effects, at the third post-Minkowskian (3PM) order. As an historically instructive simpler example, we consider here the analogous problem in electromagnetism in flat spacetime. We compute for the first time the changes in linear momentum of each particle and the total radiated linear momentum, in the relativistic classical scattering of two point-charges, at sixth order in the charges (analogous to 3PM order in gravity). We accomplish this here via direct iteration of the classical equations of motion, while making comparisons where possible to results from quantum scattering amplitudes, with the aim of contributing to the elucidation of conceptual issues and scalability on both sides. We also discuss further extensions to radiative quantities of recently established relations which analytically continue certain observables from the scattering regime to the regime of bound orbits, applicable for both the electromagnetic and gravitational cases.