Computability of entropy and information in classical Hamiltonian systems

We consider the computability of entropy and information in classical Hamiltonian systems. We define the information part and total information capacity part of entropy in classical Hamiltonian systems using relative information under a computable discrete partition. Using a recursively enumerable nonrecursive set it is shown that even though the initial probability distribution, entropy, Hamiltonian and its partial derivatives are computable under a computable partition, the time evolution of its information capacity under the original partition can grow faster than any recursive function. This implies that even though the probability measure and information are conserved in classical Hamiltonian time evolution we might not actually compute the information with respect to the original computable partition.     

Two new solutions to the third-order symplectic integration method

Two new solutions are obtained for the symplecticity conditions of explicit third-order partitioned Runge-Kutta time integration method. One of them has larger stability limit and better dispersion property than the Ruth's method.     

Computation of derivatives of the rotation number for parametric families of circle diffeomorphisms

In this paper we present a numerical method to compute derivatives of the rotation number for parametric families of circle diffeomorphisms with high accuracy. Our methodology is an extension of a recently developed approach to compute rotation numbers based on suitable averages of iterates of the map and Richardson extrapolation. We focus on analytic circle diffeomorphisms, but the method also works if the maps are differentiable enough. In order to justify the method, we also require the family of maps to be differentiable with respect to the parameters and the rotation number to be Diophantine. In particular, the method turns out to be very efficient for computing Taylor expansions of Arnold Tongues of families of circle maps. Finally, we adapt these ideas to study invariant curves for parametric families of planar twist maps.     

Numerical computation of rotation numbers of quasi-periodic planar curves

Recently, a new numerical method has been proposed to compute rotation numbers of analytic circle diffeomorphisms, as well as derivatives with respect to parameters, that takes advantage of the existence of an analytic conjugation to a rigid rotation. This method can be directly applied to the study of invariant curves of planar twist maps by simply projecting the iterates of the curve onto a circle. In this work we extend the methodology to deal with general planar maps. Our approach consists in computing suitable averages of the iterates of the map that allow us to obtain a new curve for which the direct projection onto a circle is well posed. Furthermore, since our construction does not use the invariance of the quasi-periodic curve under the map, it can be applied to more general contexts. We illustrate the method with several examples.     

A Maple Package to Compute Lie Symmetry Groups and Symmetry Reductions of （1＋1）-Dimensional Nonlinear Systems

Armed with the computer algebra system Maple, using a direct algebraic substitution method, we obtain Lie point symmetries, Lie symmetry groups and the corresponding symmetry reductions of one component nonlinear integrable and nonintegrable equations only by clicking the ‘Enter＇ key. Abundant （1＋1）-dimensional nonlinear mathematical physical systems are analysed effectively by using a Maple package LieSYMGRP proposed by us.     

Unwrapping eigenfunctions to discover the geometry of almost-invariant sets in hyperbolic maps

The numerical approximation of Perron-Frobenius operators allows efficient determination of the physical invariant measure of chaotic dynamical systems as a fixed point of the operator. Eigenfunctions of the Perron-Frobenius operator corresponding to large subunit eigenvalues have been shown to describe “almost-invariant” dynamics in one-dimensional expanding maps. We extend these ideas to hyperbolic maps in higher dimensions. While the eigendistributions of the operator are relatively uninformative, applying a new procedure called “unwrapping” to regularised versions of the eigendistributions clearly reveals the geometric structures associated with almost-invariant dynamics. This unwrapping procedure is applied to a uniformly hyperbolic map of the unit square to discover this map’s dominant underlying dynamical structure, and to the standard map to pinpoint clusters of period 6 orbits.     

A Field Method for Integrating Equations of Motion of Nonlinear Mechanico-Electrical Coupling Dynamical Systems

We deal with the generalization of the field method to weakly non-linear mechanico-electrical coupling systems. The field co-ordinates and field momenta approaches are combined with the method of multiple time scales in order to obtain the amplitudes and phase of oscillations in the first approximation. An example in mechanico-electrical coupling systems is given to illustrate this method.     

Detecting resonances in conservative maps using evolutionary algorithms

A numerical method is proposed for detecting resonances of conservative maps which reduces this task to an optimization problem. We then solve this problem using evolutionary algorithms, which are methods for global optimization inspired by biological evolution. The proposed methodology is simple and can be easily applied to maps of arbitrary dimensions. In this Letter we apply it to several examples of 2- and 4-dimensional conservative maps, with quite promising results concerning integrability, the location of resonances and the presence of chaotic regions surrounding the island chains that correspond to these resonances.     

Analytical approximation schemes for solving exact renormalization group equations in the local potential approximation

The relation between the Wilson–Polchinski and the Litim optimized ERGEs in the local potential approximation is studied with high accuracy using two different analytical approaches based on a field expansion: a recently proposed genuine analytical approximation scheme to two-point boundary value problems of ordinary differential equations, and a new one based on approximating the solution by generalized hypergeometric functions. A comparison with the numerical results obtained with the shooting method is made. A similar accuracy is reached in each case. Both two methods appear to be more efficient than the usual field expansions frequently used in the current studies of ERGEs (in particular for the Wilson–Polchinski case in the study of which they fail).     

Analytical approximation schemes for solving exact renormalization group equations. II Conformal mappings

We present a new efficient analytical approximation scheme to two-point boundary value problems of ordinary differential equations (ODEs) adapted to the study of the derivative expansion of the exact renormalization group equations. It is based on a compactification of the complex plane of the independent variable using a mapping of an angular sector onto a unit disc. We explicitly treat, for the scalar field, the local potential approximations of the Wegner–Houghton equation in the dimension d=3$d=3$ and of the Wilson–Polchinski equation for some values of d∈]2,3]$d\in ]2,3]$. We then consider, for d=3$d=3$, the coupled ODEs obtained by Morris at the second order of the derivative expansion. In both cases the fixed points and the eigenvalues attached to them are estimated. Comparisons of the results obtained are made with the shooting method and with the other analytical methods available. The best accuracy is reached with our new method which presents also the advantage of being very fast. Thus, it is well adapted to the study of more complicated systems of equations.     

Solution of singular two-point boundary value problems using differential transformation method

In this Letter, we implemented relatively new, exact series method of solution known as the differential transform method for solving singular two-point boundary value problems. Several illustrative examples are given to demonstrate the effectiveness of the present method.     

Finite difference approximations for the fractional advection-diffusion equation

Fractional order diffusion equations are viewed as generalizations of classical diffusion equations, treating super-diffusive flow processes. In this Letter, in order to solve the two-sided fractional advection-diffusion equation, the fractional Crank-Nicholson method (FCN) is given, which is based on shifted Grünwald-Letnikov formula. It is shown that this method is unconditionally stable, consistent and convergent. The accuracy with respect to the time step is of order ²(Δt). A numerical example is presented to confirm the conclusions.     

A simple approach to nonlinear oscillators

A very simple and effective approach to nonlinear oscillators is suggested. Anyone with basic knowledge of advanced calculus can apply the method to finding approximately the amplitude-frequency relationship of a nonlinear oscillator. Some examples are given to illustrate its extremely simple solution procedure and an acceptable accuracy of the obtained solutions.     

Solution of the MHD Falkner-Skan flow by Hankel-Padé method

In this Letter, an analysis is performed to find the unknown skin friction coefficient of the boundary layer Falkner-Skan equation for wedge. The governing equation takes into account the effect of a non-uniform magnetic field. The results are obtained by using Hankel-Padé method.     

Exact and numerical solutions of generalized Drinfeld-Sokolov equations

In this Letter, we consider a system of generalized Drinfeld-Sokolov (gDS) equations which models one-dimensional nonlinear wave processes in two-component media. We find some exact solutions of gDS by using tanh function method and we also obtain a numerical solution by using the Adomian's Decomposition Method (ADM).     

Differential transform method for solving linear and non-linear systems of partial differential equations

In this Letter, we propose a reliable algorithm to develop exact and approximate solutions for the linear and non-linear systems of partial differential equations. The approach rest mainly on two-dimensional differential transform method which is one of the approximate methods. The method can easily be applied to many linear and non-linear problems and is capable of reducing the size of computational work. Exact solutions can also be achieved by the known forms of the series solutions. Several illustrative examples are given to demonstrate the effectiveness of the present method.     

Maintaining the stability of nonlinear differential equations by the enhancement of HPM

Homotopy perturbation method is an effective method to find a solution of a nonlinear differential equation. In this method, a nonlinear complex differential equation is transformed to a series of linear and nonlinear parts, almost simpler differential equations. These sets of equations are then solved iteratively. Finally, a linear series of the solutions completes the answer if the convergence is maintained. In this Letter, the need for stability verification is shown through some examples. Consequently, HPM is enhanced by a preliminary assumption. The idea is to keep the inherent stability of nonlinear dynamic, even the selected linear part is not.     

Self-similar factor approximants for evolution equations and boundary-value problems

The method of self-similar factor approximants is shown to be very convenient for solving different evolution equations and boundary-value problems typical of physical applications. The method is general and simple, being a straightforward two-step procedure. First, the solution to an equation is represented as an asymptotic series in powers of a variable. Second, the series are summed by means of the self-similar factor approximants. The obtained expressions provide highly accurate approximate solutions to the considered equations. In some cases, it is even possible to reconstruct exact solutions for the whole region of variables, starting from asymptotic series for small variables. This can become possible even when the solution is a transcendental function. The method is shown to be more simple and accurate than different variants of perturbation theory with respect to small parameters, being applicable even when these parameters are large. The generality and accuracy of the method are illustrated by a number of evolution equations as well as boundary value problems.     

New universal bifurcation scenario in one-dimensional trimodal maps

By means of star products and high precision numerical calculation, an abnormal phenomenon is found in period-p-tupling bifurcation processes in one-dimensional trimodal maps. A route of transition to chaos, presented by a right-associative non-normal star product, breaks the Feigenbaum's metric universality, namely, the conventional Feigenbaum's successive rates exhibit a strong divergence. To overcome the divergence, an approximate scheme of accelerating convergence is proposed; and the Feigenbaum scenario is included as a special case in the new bifurcation scenario. It will provide access to understanding non-normal star products and their corresponding renormalization.     

A butterfly-shaped localization set for the Lorenz attractor

In this Letter, we construct a new localization set for the Lorenz system by using an extended Lyapunov method. This set reflects the butterfly-shape of the Lorenz attractor.