Boltzmann-Gibbs thermal equilibrium distribution for classical systems and Newton law: a computational discussion

We implement a general numerical calculation that allows for a direct comparison between nonlinear Hamiltonian dynamics and the Boltzmann-Gibbs canonical distribution in Gibbs Γ-space. Using paradigmatic first-neighbor models, namely, the inertial XY ferromagnet and the Fermi-Pasta-Ulam β-model, we show that at intermediate energies the Boltzmann-Gibbs equilibrium distribution is a consequence of Newton second law (F=ma). At higher energies we discuss partial agreement between time and ensemble averages.     

Geometric and dynamic perspectives on phase-coherent and noncoherent chaos

Statistically distinguishing between phase-coherent and noncoherent chaotic dynamics from time series is a contemporary problem in nonlinear sciences. In this work, we propose different measures based on recurrence properties of recorded trajectories, which characterize the underlying systems from both geometric and dynamic viewpoints. The potentials of the individual measures for discriminating phase-coherent and noncoherent chaotic oscillations are discussed. A detailed numerical analysis is performed for the chaotic Ro?ssler system, which displays both types of chaos as one control parameter is varied, and the Mackey-Glass system as an example of a time-delay system with noncoherent chaos. Our results demonstrate that especially geometric measures from recurrence network analysis are well suited for tracing transitions between spiral- and screw-type chaos, a common route from phase-coherent to noncoherent chaos also found in other nonlinear oscillators. A detailed explanation of the observed behavior in terms of attractor geometry is given.     

Canonical methods for Hamiltonian systems: numerical experiments

A Hamiltonian system possesses dynamics (e.g. preservation of volume in phase space and symplectic structure) that call for special numerical integrators, namely canonical methods. Recent research on this aspect have shown that canonical numerical integrators may be needed for Hamiltonian systems. In this paper, we focus on numerical experiments that compare canonical and non-canonical numerical integrators. Test problems are taken from different areas in physical sciences. These experiments help to buttress the claims that canonical numerical integrators give results that mimic the qualitative behavior of the original system and that canonical numerical integrators are suitable for long time integrations. Our experiments indicate that higher-order canonical methods allow for larger timestep than lower-order canonical methods.     

Chaotic diffusion for particles moving in a time dependent potential well

The chaotic diffusion for particles moving in a time dependent potential well is described by using two different procedures: (i) via direct evolution of the mapping describing the dynamics and; (ii) by the solution of the diffusion equation. The dynamic of the diffusing particles is made by the use of a two dimensional, nonlinear area preserving map for the variables energy and time. The phase space of the system is mixed containing both chaos, periodic regions and invariant spanning curves limiting the diffusion of the chaotic particles. The chaotic evolution for an ensemble of particles is treated as random particles motion and hence described by the diffusion equation. The boundary conditions impose that the particles can not cross the invariant spanning curves, serving as upper boundary for the diffusion, nor the lowest energy domain that is the energy the particles escape from the time moving potential well. The diffusion coefficient is determined via the equation of the mapping while the analytical solution of the diffusion equation gives the probability to find a given particle with a certain energy at a specific time. The momenta of the probability describe qualitatively the behavior of the average energy obtained by numerical simulation, which is investigated either as a function of the time as well as some of the control parameters of the problem.     

Molecular dynamics and time reversibility

We present a time-symmetrical integer arithmetic algorithm for numerical (molecular dynamics) simulations of classical fluids. This algorithm is used to illustrate, through concrete examples, that time-asymmetric evolutions are typical for systems of many particles evolving according to reversible microscopic dynamics and to calculate the asymptotic behavior of the velocity autocorrelation function with an improved accuracy. The equivalence between equilibrium time averages and microcanonical ensemble averages is checked via two new sampling methods for computing microcanonical averages of classical systems.     

Chaotic and Arnold stripes in weakly chaotic Hamiltonian systems

The dynamics in weakly chaotic Hamiltonian systems strongly depends on initial conditions (ICs) and little can be affirmed about generic behaviors. Using two distinct Hamiltonian systems, namely one particle in an open rectangular billiard and four particles globally coupled on a discrete lattice, we show that in these models, the transition from integrable motion to weak chaos emerges via chaotic stripes as the nonlinear parameter is increased. The stripes represent intervals of initial conditions which generate chaotic trajectories and increase with the nonlinear parameter of the system. In the billiard case, the initial conditions are the injection angles. For higher-dimensional systems and small nonlinearities, the chaotic stripes are the initial condition inside which Arnold diffusion occurs.     

Quantum and classical spin clusters: disappearance of quantum numbers and Hamiltonian chaos

We present a direct link between manifestations of classical Hamiltonian chaos and quantum nonintegrability effects as they occur in quantum invariants. In integrable classical Hamiltonian systems, analytic invariants (integrals of the motion) can be constructed numerically by means of time averages of dynamical variables over phase-space trajectories, whereas in near-integrable models such time averages yield nonanalytic invariants with qualitatively different properties. Translated into quantum mechanics, the invariants obtained from time averages of dynamical variables in energy eigenstates provide a topographical map of the plane of quantized actions (quantum numbers) with properties which again depend sensitively on whether or not the classical integrability condition is satisfied. The most conspicuous indicator of quantum chaos is the disappearance of quantum numbers, a phenomenon directly related to the breakdown of invariant tori in the classical phase flow. All results are for a system consisting of two exchange-coupled spins with biaxial exchange and single-site anisotropy, a system with a nontrivial integrability condition.     

Three types of chaos in the forced nonlinear Schrödinger equation

Three different types of chaotic behavior and instabilities (homoclinic chaos, hyperbolic resonance, and parabolic resonance) in Hamiltonian perturbations of the nonlinear Schr?dinger (NLS) equation are described. The analysis is performed on a truncated model using a novel framework in which a hierarchy of bifurcations is constructed. It is demonstrated numerically that the forced NLS equation exhibits analogous types of chaotic phenomena. Thus, by adjusting the forcing frequency, the behavior near the plane wave solution may be set to any one of the three different types of chaos for any periodic box length.     

Chaos suppression of uncertain gyros in a given finite time

The gyro is one of the most interesting and everlasting nonlinear dynamical systems,which displays very rich and complex dynamics,such as sub-harmonic and chaotic behaviors.We study the chaos suppression of the chaotic gyros in a given finite time.Considering the effects of model uncertainties,external disturbances,and fully unknown parameters,we design a robust adaptive finite-time controller to suppress the chaotic vibration of the uncertain gyro as quickly as possible.Using the finite-time control technique,we give the exact value of the chaos suppression time.A mathematical theorem is presented to prove the finite-time stability of the proposed scheme.The numerical simulation shows the efficiency and usefulness of the proposed finite-time chaos suppression strategy.     

Using the extended Melnikov method to study the multi-pulse global bifurcations and chaos of a cantilever beam

The aim of this paper is to investigate the multi-pulse global bifurcations and chaotic dynamics for the nonlinear non-planar oscillations of a cantilever beam subjected to a harmonic axial excitation and two transverse excitations at the free end by using an extended Melnikov method in the resonant case. First, the extended Melnikov method for studying the Shilnikov-type multi-pulse homoclinic orbits and chaos in high-dimensional nonlinear systems is briefly introduced in the theoretical frame. Then, this method is utilized to investigate the Shilnikov-type multi-pulse homoclinic bifurcations and chaotic dynamics for the nonlinear non-planar oscillations of the cantilever beam. How to employ this method to analyze the Shilnikov-type multi-pulse homoclinic bifurcations and chaotic dynamics of high-dimensional nonlinear systems in engineering applications is demonstrated through this example. Finally, the results of numerical simulation are given and also show that the Shilnikov-type multi-pulse chaotic motions can occur for the nonlinear non-planar oscillations of the cantilever beam, which verifies the analytical prediction.     

An Improved Calculation Formula of the Extended Entropic Chaos Degree and Its Application to Two-Dimensional Chaotic Maps

The Lyapunov exponent is primarily used to quantify the chaos of a dynamical system. However, it is difficult to compute the Lyapunov exponent of dynamical systems from a time series. The entropic chaos degree is a criterion for quantifying chaos in dynamical systems through information dynamics, which is directly computable for any time series. However, it requires higher values than the Lyapunov exponent for any chaotic map. Therefore, the improved entropic chaos degree for a one-dimensional chaotic map under typical chaotic conditions was introduced to reduce the difference between the Lyapunov exponent and the entropic chaos degree. Moreover, the improved entropic chaos degree was extended for a multidimensional chaotic map. Recently, the author has shown that the extended entropic chaos degree takes the same value as the total sum of the Lyapunov exponents under typical chaotic conditions. However, the author has assumed a value of infinity for some numbers, especially the number of mapping points. Nevertheless, in actual numerical computations, these numbers are treated as finite. This study proposes an improved calculation formula of the extended entropic chaos degree to obtain appropriate numerical computation results for two-dimensional chaotic maps.     

Communicating with optical hyperchaos: information encryption and decryption in delayed nonlinear feedback systems

Recent theoretical studies and experimental demonstrations have shown the possibility of using chaos for the encryption of message signals in communication systems. Chaos is generated by systems with delayed nonlinear feedback, which feature hyperchaotic (i.e., of high dimensionality) dynamics. The different ways for the injection of the information in the emitter and the process of the synchronization of the receiver are considered. The analysis of all the possibilities can be used to choose the correct topology of communication systems and, more generally, to explain the behavior of any chaotic systems ruled by nonlinear difference-differential equations.     

基于非线性控制的超混沌Chen系统混沌同步

研究了基于非线性控制的超混沌Chen系统的混沌同步问题.基于Lyapunov稳定性理论，设计了非线性控制器，改进了Jiang和Huang等设计的同步误差系统的Lyapunov函数形式，理论证明了超混沌Chen系统的自同步和超混沌Chen系统与超混沌R?ssler系统的异结构同步.数值模拟进一步验证了所提出方案的有效性. 关键词： 超混沌Chen系统 自同步 异结构同步 非线性控制器     

Chaotic solitons in the quadratic-cubic nonlinear Schrödinger equation under nonlinearity management

We analyze the response of rational and regular (hyperbolic-secant) soliton solutions of an extended nonlinear Schro?dinger equation (NLSE) which includes an additional self-defocusing quadratic term, to periodic modulations of the coefficient in front of this term. Using the variational approximation (VA) with rational and hyperbolic trial functions, we transform this NLSE into Hamiltonian dynamical systems which give rise to chaotic solutions. The presence of chaos in the variational solutions is corroborated by calculating their power spectra and the correlation dimension of the Poincare? maps. This chaotic behavior (predicted by the VA) is not observed in the direct numerical solutions of the NLSE when rational initial conditions are used. The solitary-wave solutions generated by these initial conditions gradually decay under the action of the nonlinearity management. On the contrary, the solutions of the NLSE with exponentially localized initial conditions are robust solitary-waves with oscillations consistent with a chaotic or a complex quasiperiodic behavior.     

A Novel Measure Inspired by Lyapunov Exponents for the Characterization of Dynamics in State-Transition Networks

The combination of network sciences, nonlinear dynamics and time series analysis provides novel insights and analogies between the different approaches to complex systems. By combining the considerations behind the Lyapunov exponent of dynamical systems and the average entropy of transition probabilities for Markov chains, we introduce a network measure for characterizing the dynamics on state-transition networks with special focus on differentiating between chaotic and cyclic modes. One important property of this Lyapunov measure consists of its non-monotonous dependence on the cylicity of the dynamics. Motivated by providing proper use cases for studying the new measure, we also lay out a method for mapping time series to state transition networks by phase space coarse graining. Using both discrete time and continuous time dynamical systems the Lyapunov measure extracted from the corresponding state-transition networks exhibits similar behavior to that of the Lyapunov exponent. In addition, it demonstrates a strong sensitivity to boundary crisis suggesting applicability in predicting the collapse of chaos.     

Stochastic-dynamical thermostats for constraints and stiff restraints

A broad array of canonical sampling methods are available for molecular simulation based on stochastic-dynamical perturbation of Newtonian dynamics, including Langevin dynamics, Stochastic Velo- city Rescaling, and methods that combine Nosé-Hoover dynamics with stochastic perturbation. In this article we discuss several stochastic-dynamical thermostats in the setting of simulating systems with holonomic constraints. The approaches described are easily implemented and facilitate the recovery of correct canonical averages with minimal disturbance of the underlying dynamics. For the purpose of illustrating our results, we examine the numerical application of these methods to a simple atomic chain, where a Fixman term is required to correct the thermodynamic ensemble.