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.     

A Cartesian treecode for screened coulomb interactions

A treecode algorithm is presented for evaluating electrostatic potentials in a charged particle system undergoing screened Coulomb interactions in 3D. The method uses a far-field Taylor expansion in Cartesian coordinates to compute particle–cluster interactions. The Taylor coefficients are evaluated using new recurrence relations which permit efficient computation of high order approximations. Two types of clusters are considered, uniform cubes and adapted rectangular boxes. The treecode error, CPU time and memory usage are reported and compared with direct summation for randomly distributed particles inside a cube, on the surface of a sphere and on an 8-sphere configuration. For a given order of Taylor approximation, the treecode CPU time scales as O(NlogN)$O(N\log N)$ and the memory usage scales as O(N)$O(N)$, where N is the number of particles. Results show that the treecode is well suited for non-homogeneous particle distributions as in the sphere and 8-sphere test cases.     

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.     

A Variational Iteration Solving Method for a Class of Generalized Boussinesq Equations

We study a generalized nonlinear Boussinesq equation by introducing a proper functional and constructing the variational iteration sequence with suitable initial approximation. The approximate solution is obtained for the solitary wave of the Boussinesq equation with the variational iteration method.     

A high-order energy-conserving integration scheme for Hamiltonian systems

A new energy-conserving numerical integration method for Hamiltonian systems is presented. The method is constructed by a parallel connection of n multi-stage schemes of order 2 and its order of accuracy is 2n.     

Method of self-similar factor approximants

The method of self-similar factor approximants is completed by defining the approximants of odd orders, constructed from the power series with the largest term of an odd power. It is shown that the method provides good approximations for transcendental functions. In some cases, just a few terms in a power series make it possible to reconstruct a transcendental function exactly. Numerical convergence of the factor approximants is checked for several examples. A special attention is paid to the possibility of extrapolating the behavior of functions, with arguments tending to infinity, from the related asymptotic series at small arguments. Applications of the method are thoroughly illustrated by the examples of several functions, nonlinear differential equations, and anharmonic models.     

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.     

Linearization of Systems of Nonlinear Diffusion Equations

We investigate the linearization of systems of n-component nonlinear diffusion equations; such systems have physical applications in soil science, mathematical biology and invariant curve flows. Equivalence transformations of their auxiliary systems are used to identify the systems that can be linearized. We also provide several examples of systems with two-component equations, and show how to linearize them by nonlocal mappings.     

A New Method for Constructing Travelling Wave Solutions to the modified Benjamin--Bona--Mahoney Equation

We present a new method to find the exact travelling wave solutions of nonlinear evolution equations, with the aid of the symbolic computation. Based on this method, we successfully solve the modified BenjaminBona-Mahoney equation, and obtain some new solutions which can be expressed by trigonometric functions and hyperbolic functions, It is shown that the proposed method is direct, effective and can be used for many other nonlinear evolution equations in mathematical physics.     

A Time-Dependent Wavepacket Method for Photodissociation Dynamics of Triatomic Molecule

We report a time-dependent quantum wavepacket theory employed to interpret the photoabsorption spectrum of the N20 molecule in terms of the nuclear motion on the upper 21A＇ and 11A＂ potential energy surfaces. The N2-O bond breaks upon excitation leading to dissociation. The total angular momentum is treated correctly taking into account the vector property of the electric field of the exciting radiation.     

Numerical solution of 3D Green's function for the dynamic system of anisotropic elasticity

In this Letter, we used homotopy perturbation method to obtain numerical solution of the 3D Green's function for the dynamic system of anisotropic elasticity. Application of homotopy perturbation method to this problem shows the rapid convergence of the sequence constructed by this method to the exact solution. The numerical results obtained from convolution of Green's function and data of the Cauchy problem are compared with the exact solution for cubic media. The results reveal that the proposed method is very effective and simple.     

Numerical computation of the asymptotic size of the rotation domain for the Arnold family

We consider the Arnold Tongue of the Arnold family of circle maps associated to a fixed Diophantine rotation number θ. The corresponding maps of the family are analytically conjugate to a rigid rotation. This conjugation is defined on a (maximal) complex strip of the circle and, after a suitable scaling, the size of this strip is given by an analytic function of the perturbative parameter.The main purpose of this paper is to perform a numerical accurate computation of this function and of its Taylor expansion. This allows us to verify previous theoretical results. The rotation numbers we select are quadratic irrationals, mainly the Golden Mean.By introducing a nonstandard extrapolation process, especially suited for the problem, we compute all the quantities required (rotation numbers, Arnold Tongues, Fourier and Taylor coefficients) with high precision.     

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.     

Iterative approximation of limit cycles for a class of Abel equations

This article considers the analytical approximation of limit cycles that may appear in Abel equations written in the normal form. The procedure uses an iterative approach that takes advantage of the contraction mapping theorem. Thus, the obtained sequence exhibits uniform convergence to the target periodic solution. The effectiveness of the technique is illustrated through the approximation of an unstable limit cycle that appears in an Abel equation arising in a tracking control problem that affects an elementary, second-order bilinear power converter.     

Two Kinds of Square-Conservative Integrators for Nonlinear Evolution Equations

Based on the Lie-group and Gauss-Legendre methods, two kinds of square-conservative integrators for square- conservative nonlinear evolution equations are presented. Lie-group based square-conservative integrators are linearly implicit, while Gauss-Legendre based square-conservative integrators are nonlinearly implicit and iterative schemes are needed to solve the corresponding integrators. These two kinds of integrators provide natural candidates for simulating square-conservative nonlinear evolution equations in the sense that these integrators not only preserve the square-conservative properties of the continuous equations but also are nonlinearly stable. Numerical experiments are performed to test the presented integrators.     

A New Solution to Einstein's Field Equations

We construct a new exact solution to the vacuum Einstein field equations. This solution possesses a naked physical singularity. The norm of the Riemann curvature tensor of the solution takes infinity at some points and the solution does not have any event horizon around the singularity. A detailed analysis of this new singularity is also presented.     

Numerical Solution of the Three-Dimensional Helmholtz Equation

We apply the homotopy analysis method to solve the nonhomogeneous multidimensional partial differential equation model problem. The analytic solutions are calculated in terms of convergent series with easily computable components. The nonhomogeneous problem is quickly solved by observing the self-canceling ＂noise＂ terms whose sum vanishes in the limit. Numerical results clearly reveal the complete reliability and efficiency of the proposed algorithm.     

Two-component description of dynamical systems that can be approximated by solitons: The case of the ion acoustic wave equations of plasma physics

A new approach to the perturbative analysis of dynamical systems, which can be described approximately by soliton solutions of integrable non-linear wave equations, is employed in the case of small-amplitude solutions of the ion acoustic wave equations of plasma physics. Instead of pursuing the traditional derivation of a perturbed KdV equation, the ion velocity is written as a sum of two components: elastic and inelastic. In the single-soliton case, the elastic component is the full solution. In the multiple-soliton case, it is complemented by the inelastic component. The original system is transformed into two evolution equations: An asymptotically integrable Normal Form for ordinary KdV solitons, and an equation for the inelastic component. The zero-order term of the elastic component is a single-soliton or multiple-soliton solution of the Normal Form. The inelastic component asymptotes into a linear combination of single-soliton solutions of the Normal Form, with amplitudes determined by soliton interactions, plus a second-order decaying dispersive wave. Satisfaction of a conservation law by the inelastic component and of mass conservation by the disturbance to the ion density is determined solely by the initial data and/or boundary conditions imposed on the inelastic component. The electrostatic potential is a first-order quantity. It is affected by the inelastic component only in second order. The charge density displays a triple-layer structure. The analysis is carried out through the third order.