Towards investigation of evolution of dynamical systems with independence of time accuracy: more classes of systems

The recently developed method (Paper 1) enabling one to investigate the evolution of dynamical systems with an accuracy not dependent on time is developed further. The classes of dynamical systems which can be studied by that method are much extended, now including systems that are; (1) non-Hamiltonian, conservative; (2) Hamiltonian with time-dependent perturbation; (3) non-conservative (with dissipation). These systems cover various types of N-body gravitating systems of astrophysical and cosmological interest, such as the orbital evolution of planets, minor planets, artificial satellites due to tidal, non-tidal perturbations and thermal thrust, evolving close binary stellar systems, and the dynamics of accretion disks.


Introduction
The typical procedure used in simulations of N-body gravitating systems [1] is the numerical iterative integration of the equations of motion, but this has a well-known difficulty: the inevitable accumulation of errors as the number of iterations increases and hence beginning from a distinct time further computations become meaningless. The difficulty is serious in the case of typical nonlinear and/or many dimensional systems covering, as now is understood, most of the physical systems. The systems of N gravitating bodies involving various astrophysical problems possess chaotic properties (e.g. [2,3,4,5,6,7]), which make their long term simulations rather complex.
Recently a method of numerical investigation was developed [8], providing a principal possibility for determination of the evolution of dynamical systems with an accuracy not dependent on time: an example of an evolution of each parameter of a system at given accuracy was exhibited in [8]. This method was here shown to be applicable for Hamiltonian systems with small periodic perturbation not dependent on time [9] H(I, ϑ, β) = H 0 (I) + βH 1 (I, ϑ) , where H 0 (I) is the integrable Hamiltonian, I, ϑ are the action-angle variables and β is the small parameter of the non-integrable Hamiltonian H 1 .
Here we represent a further development of that method extending the class of dynamical systems which can be investigated, paying mainly attention to the systems having primary physical and astrophysical interest. We show that systems of rather general classes can be principally studied, including systems that are: 1. non-Hamiltonian, conservative; 2. non-conservative (with dissipation); 3. Hamiltonian, with time-dependent perturbation.
In our analysis we essentially use a result by Chernoff [10] concerning the properties of operators on a complete Riemannian manifold.
Our initial aim is to investigate a rather general type first order differential equationṡ The main idea of the method developed in [8] was the reduction of this equation to the following one where L is a differential operator. The solution of this equation by means of the computation of the resolvent using computer algebraic analytical codes gives the evolution of the function u(t) in time. Then the functions x a (t), being determined by the function u(t) and describing the initial system, can be found (see [8]).

The resolvent
Let M be a d-dimensional complete Riemannian manifold with a metric g. Denote by F (M ) the set of smooth functions determined on M . Let F 0 g(M ) be the subset of F (M ) which includes functions having a compact support. On F 0 (M ) one has the operation of inner product where u and v are elements of F 0 (M ). Denote by L 2 (M ) the completion of F 0 (M ) with respect to that product. Consider the operator in the form where f a (x) are smooth functions. One can show that the following relation for this operator is fulfilled [9] (u, where u, v ∈ F 0 (M ), divf = (det g) −1/2 ∂ a ((det g) 1/2 f a ) = f a ; a .
When divf = 0 , then (u, Lv) = (Lu, v) , i.e. L is a formal self-adjoint operator. In our case the velocity of the propagation c(x) introduced by Chernoff [10] is equal to According to Chernoff (Theorem 2.2 in [10]), if M is a complete Riemannian manifold and where c(r) = sup{c(x); x ∈ S(y, r)}, S(y, r) is a ball of radius r with center y, then the operator L defined on space L 2 (M ) with domain F 0 (M ) is essentially self-adjoint. Consider now Eq.(3), to which, as mentioned above, Eq.(2) is reduced. The solution of that equation can be found by means of the Laplace transform [8], where R λ (β) is the resolvent, Let us first consider the case when f a (x) = ξ a + βζ a (x) (13) and then the resolvent is well approximated by formula [8] where In the case when for the dynamical system defined on M we shall consider another one, having in local coordinates the forṁ Here the manifold M is reduced to M × R with the metric In the particular case that andθ a = ξ a (I), a = 1, . . . , k , the resolvent R λ (β) can be found as in [8].

Conclusions
Let us summarize the dynamical systems covered by the above analysis. The systems represented by equationsẋ a = f a (x) , include: 1. Systems for which the Liouville theorem is fulfilled, This problem is reduced to finding out of the resolvent R λ (0), including the case when: a) the numbers of action-angle variables do not equal to each other. The calculations of R λ (0) should be considered separately at each given case. For example, R λ (0) is easily found and therefore the method is trivially applicable when b) f a is a periodic function; 2. Systems with perturbed Hamiltonian H(I, ϑ, t, β) = H 0 (I) + βH 1 (I, ϑ, t).
When H(I, ϑ, t, β) is a periodic function of time the resolvent is trivially represented via Fourier series (cf. [8]). In other cases the R λ (0) must be calculated separately; 3. Non-conservative systems, when divf = 0 .
This case is reduced to P.1, as shown above. The evolution of all these dynamical systems can be investigated by the method developed in [8] with an accuracy not dependent on the value of time. Evidently, now the classes of physical and astrophysical problems are far extended. Namely, besides traditional astrophysical nonlinear dynamical systems -N-body gravitating systems -now this will include an entire bunch of problems, e.g. the non-Hamiltonian motion of minor planets and space probes at tidal, nontidal perturbation modes, and thermal thrust; these are particularly important, for example, for accurately testing of frame dragging, e.g. [11]. Other examples of non-conservative dissipative systems include e.g. the orbital evolution of close binary stellar systems due to evolving tidal effects, the magnetohydrodynamical evolution of non-stationary accretion disks near black holes, see e.g. [12,13,14].