An application of symplectic integration for general relativistic planetary orbitography subject to non-gravitational forces

Spacecraft propagation tools describe the motion of near-Earth objects and interplanetary probes using Newton’s theory of gravity supplemented with the approximate general relativistic n-body Einstein–Infeld–Hoffmann equations of motion. With respect to the general theory of relativity and the long-standing recommendations of the International Astronomical Union for astrometry, celestial mechanics and metrology, we believe modern orbitography software is now reaching its limits in terms of complexity. In this paper, we present the first results of a prototype software titled General Relativistic Accelerometer-based Propagation Environment (GRAPE). We describe the motion of interplanetary probes and spacecraft using extended general relativistic equations of motion which account for non-gravitational forces using end-user supplied accelerometer data or approximate dynamical models. We exploit the unique general relativistic quadratic invariant associated with the orthogonality between four-velocity and acceleration and simulate the perturbed orbits for Molniya, Parker Solar Probe and Mercury Planetary Orbiter-like test particles subject to a radiation-like four-force. The accuracy of the numerical procedure is maintained using a 5-stage, $$10^\mathrm{th}$$ 10 th -order structure-preserving Gauss collocation symplectic integration scheme. GRAPE preserves the norm of the tangent vector to the test particle worldline at the order of $$10^{-32}$$ 10 - 32 .

Energy Mechanisms of Free Vibrations and Resonance in Elastic Bodies

The scientific novelty of this work is determined by the rationale for the participation in transformations, along with the kinetic energy of particles, of four types of elastic energy, identified by the peculiarities of their phase changes in the oscillation process. Two types are converted into kinetic energy, while the other two types change the deformed state of particles in accordance with the equations of motion due to internal sources. The result is obtained based on the use of the superposition principle in the space of Lagrange variables with the imposition of forced and free oscillations, as well as a new model of mechanics based on the concepts of space, time, and energy with a new scale of average stresses that takes into account the energy of particles in the initial state. In such a model of mechanics, a generalized measure of the elastic energy of particles is a quadratic invariant of asymmetric tensor whose components are partial derivatives of Euler variables with respect to Lagrange variables. The concept of kinematic energy parameters is introduced, which differ from the corresponding volumetric energy densities by a multiplier equal to the modulus of elasticity, which is directly proportional to the density and heat capacity of the material, and inversely proportional to the volumetric compression coefficient. Comparison of the values of kinematic parameters shows that most of the energy required for oscillations is associated with the deformation of particles and comes from internal sources. The mechanisms of transformation of forced vibrations into their own for transverse, torsional, and longitudinal vibrations are considered, as well as the occurrence of resonance when free and forced vibrations are superimposed with the same or a similar frequency. The formation of a new free wave after each cycle of external influences with an increase in amplitude, which occurs mainly due to internal, and not external, energy sources is justified.

Dynamical Invariant Applied on General Time-Dependent Three Coupled Nano-Optomechanical Oscillators

A quadratic invariant operator for general time-dependent three coupled nano-optomechanical oscillators is investigated. We show that the invariant operator that we have established satisfies the Liouville-von Neumann equation and coincides with its classical counterpart. To diagonalize the invariant, we carry out a unitary transformation of it at first. From such a transformation, the quantal invariant operator reduces to an equal, but a simple one which corresponds to three coupled oscillators with time-dependent frequencies and unit masses. Finally, we diagonalize the matrix representation of the transformed invariant by using a unitary matrix. The diagonalized invariant is just the same as the Hamiltonian of three simple oscillators. Thanks to such a diagonalization, we can analyze various dynamical properties of the nano-optomechanical system. Quantum characteristics of the system are investigated as an example, by utilizing the diagonalized invariant. We derive not only the eigenfunctions of the invariant operator, but also the wave functions in the Fock state.

Linear System of Differential Equations with a Quadratic Invariant as the Schrödinger Equation

Linear systems of differential equations with an invariant in the form of a positive definite quadratic form in a real Hilbert space are considered. It is assumed that the system has a simple spectrum and the eigenvectors form a complete orthonormal system. Under these assumptions, the linear system can be represented in the form of the Schrödinger equation by introducing a suitable complex structure. As an example, we present such a representation for the Maxwell equations without currents. In view of these observations, the dynamics defined by some linear partial differential equations can be treated in terms of the basic principles and methods of quantum mechanics.

On the realization of explicit Runge-Kutta schemes preserving quadratic invariants of dynamical systems

We implement several explicit Runge-Kutta schemes that preserve quadratic invariants of autonomous dynamical systems in Sage. In this paper, we want to present our package ex.sage and the results of our numerical experiments. In the package, the functions rrk_solve, idt_solve and project_1 are constructed for the case when only one given quadratic invariant will be exactly preserved. The function phi_solve_1 allows us to preserve two specified quadratic invariants simultaneously. To solve the equations with respect to parameters determined by the conservation law we use the elimination technique based on Grbner basis implemented in Sage. An elliptic oscillator is used as a test example of the presented package. This dynamical system has two quadratic invariants. Numerical results of the comparing of standard explicit Runge-Kutta method RK(4,4) with rrk_solve are presented. In addition, for the functions rrk_solve and idt_solve, that preserve only one given invariant, we investigated the change of the second quadratic invariant of the elliptic oscillator. In conclusion, the drawbacks of using these schemes are discussed.

Construction of rogue waves and conservation laws of the complex coupled Kadomtsev–Petviashvili equation

In this work, the solutions of the real coupled Kadomtsev–Petviashvili equations (KPEs) are obtained and they are used to find the ones of the complex system. To this end, the extension theorem is used. The formation of rogue wave (RW) is shown to occur via inelastic collisions of lines-lump solitons and periodic waves and also of two elliptic waves. The quadratic invariant is used to show the phase portrait and bifurcation through the contour plot. Further, the conservation laws (CLs) are constructed by means of new conservation theorem by finding the governing system of equation. The results are given, work consolidated from a previous work, by the first two authors, for suggesting a mechanism for the formation of RW when the complex Korteweg–de Vries equation was studied. Based on the exact solutions found here, numerical computations are carried out and they are represented graphically. Finally, the unified method (UM) can be applied to solve a variety of partial differential equations in science and engineering.

Decreasing height along continued fractions

The fact that the euclidean algorithm eventually terminates is pervasive in mathematics. In the language of continued fractions, it can be stated by saying that the orbits of rational points under the Gauss map$x\mapsto x^{-1}-\lfloor x^{-1}\rfloor$eventually reach zero. Analogues of this fact for Gauss maps defined over quadratic number fields have relevance in the theory of flows on translation surfaces, and have been established via powerful machinery, ultimately relying on the Veech dichotomy. In this paper, for each commensurability class of non-cocompact triangle groups of quadratic invariant trace field, we construct a Gauss map whose defining matrices generate a group in the class; we then provide a direct and self-contained proof of termination. As a byproduct, we provide a new proof of the fact that non-cocompact triangle groups of quadratic invariant trace field have the projective line over that field as the set of cross-ratios of cusps. Our proof is based on an analysis of the action of non-negative matrices with quadratic integer entries on the Weil height of points. As a consequence of the analysis, we show that long symbolic sequences in the alphabet of our maps can be effectively split into blocks of predetermined shape having the property that the height of points which obey the sequence and belong to the base field decreases strictly at each block end. Since the height cannot decrease infinitely, the termination property follows.