Equilibrium-torus bifurcation in nonsmooth systems

Considering a set of two coupled nonautonomous differential equations with discontinuous right-hand sides describing the behavior of a DC/DC power converter, we discuss a border-collision bifurcation that can lead to the birth of a two-dimensional invariant torus from a stable node equilibrium point. We obtain the chart of dynamic modes and show that there is a region of parameter space in which the system has a single stable node equilibrium point. Under variation of the parameters, this equilibrium may disappear as it collides with a discontinuity boundary between two smooth regions in the phase space. The disappearance of the equilibrium point is accompanied by the soft appearance of an unstable focus period-1 orbit surrounded by a resonant or ergodic torus.Detailed numerical calculations are supported by a theoretical investigation of the normal form map that represents the piecewise linear approximation to our system in the neighbourhood of the border. We determine the functional relationships between the parameters of the normal form map and the actual system and illustrate how the normal form theory can predict the bifurcation behaviour along the border-collision equilibrium-torus bifurcation curve.     

Center manifolds,normal forms and bifurcations of vector fields with application to coupling between periodic and steady motions

We study the instabilities known to aeronautical engineers as flutter and divergence. Mathematically, these states correspond to bifurcations to limit cycles and multiple equilibrium points in a differential equation. Making use of the center manifold and normal form theorems, we concentrate on the situation in which flutter and divergence become coupled, and show that there are essentially two ways in which this is likely to occur. In the first case the system can be reduced to an essential model which takes the form of a single degree of freedom nonlinear oscillator. This system, which may be analyzed by conventional phase-plane techniques, captures all the qualitative features of the full system. We discuss the reduction and show how the nonlinear terms may be simplified and put into normal form. Invariant manifold theory and the normal form theorem play a major role in this work and this paper serves as an introduction to their application in mechanics.Repeating the approach in the second case, we show that the essential model is now three dimensional and that far more complex behavior is possible, including nonperiodic and ‘chaotic’ motions. Throughout, we take a two degree of freedom system as an example, but the general methods are applicable to multi- and even infinite degree of freedom problems.     

A new method of analysis of the effect of weak colored noise in nonlinear dynamical systems

A systematic method for obtaining the asymptotic behavior of a dynamical system forced by colored noise in the limit of small intensity is developed. It is based on the search of WKB solutions to the Fokker-Planck equation for the joint probability density of the system and noise, in which the perturbation expansion is continued to the first correction beyond the Hamilton-Jacobi limit. The method can be applied to noise with correlation time of order unity. It is illustrated on the normal form of a pitchfork bifurcation, where it is pointed out that additive noise can induce a shift of the most probable value. This prediction is confirmed by numerical simulation of the stochastic differential equations.     

Synchronization of globally coupled phase oscillators: singularities and scaling for general couplings

The onset of collective behavior in a population of globally coupled oscillators with randomly distributed frequencies is studied for phase dynamical models with arbitrary coupling; the effect of a stochastic temporal variation in the frequencies is also included. The Fokker-Planck equation for the coupled Langevin system is reduced to a kinetic equation for the oscillator distribution function. Instabilities of the phase-incoherent state are studied by center manifold reduction to the amplitude dynamics of the unstable modes. Depending on the coupling, the coefficients in the normal form can be singular in the limit of weak instability when the diffusive effect of the noise is neglected. A detailed analysis of these singularities to all orders in the normal form expansion is presented. Physically, the singularities are interpreted as predicting an altered scaling of the entrained component near the onset of synchronization. These predictions are verified by numerically solving the kinetic equation for various couplings and frequency distributions.     

Asymptotic chaos

We provide in table I a list of normal forms of ordinary differential equations describing the dynamics of physical systems in conditions near to the simultaneous onset of up three instabilities. The first (quadratic) terms in the Taylor series for the nonlinear terms in these amplitude equations (as they are called in fluid dynamics) are given in each case. We focus on a particular example involving three instabilities and derive an asymptotic version of the corresponding normal form as a limit of small dissipation is approached. The numerical investigation of this asymptotic normal form strongly suggests that chaotic behavior occurs as close as one wants to the onset of the triple instability. This chaotic behavior is also exhibited by a return map constructed by direct numerical experiments on the amplitude equation. We also derive by analytic methods a model return map that qualitatively reproduces much of the dynamics observed numerically in the solutions of the asymptotic normal form in nearly homoclinic conditions. In the limit of strong contraction, this model map of the plane reduces to a unidimensional map that is valuable in understanding the dynamics of the original system.     

A normal form for excitable media

We present a normal form for traveling waves in one-dimensional excitable media in the form of a differential delay equation. The normal form is built around the well-known saddle-node bifurcation generically present in excitable media. Finite wavelength effects are captured by a delay. The normal form describes the behavior of single pulses in a periodic domain and also the richer behavior of wave trains. The normal form exhibits a symmetry preserving Hopf bifurcation which may coalesce with the saddle node in a Bogdanov-Takens point, and a symmetry-breaking spatially inhomogeneous pitchfork bifurcation. We verify the existence of these bifurcations in numerical simulations. The parameters of the normal form are determined and its predictions are tested against numerical simulations of partial differential equation models of excitable media with good agreement.     

A study of the asymmetric Malkus waterwheel: the biased Lorenz equations

In this work, the asymmetric case of the Malkus waterwheel is studied, where the water inflow to the system is biasing the system toward stable motion in one direction, like a Pelton wheel. The governing equations of this system, when expressed in Fourier space and decoupled to form a closed set, can be mapped into a four-dimensional space where they form a quasi-Lorenz system. This set of equations is analyzed in light of analogues of the Rayleigh Bernard convection and conclusions are drawn. The properties and behavior of the equations are studied and correlated to the physical model. Phase space behavior and linear stability analysis are used for this. Spectral analysis is used as a qualitative measure of chaos. Chaotic behavior is quantified through the calculation of the Lyapunov exponents and these are further correlated to the bifurcation diagrams for a conclusive analysis of the dynamical behavior of the system.     

Control of degenerate Hopf bifurcations in three-dimensional maps

A feedback control method is proposed to create a degenerate Hopf bifurcation in three-dimensional maps at a desired parameter point. The particularity of this bifurcation is that the system admits a stable fixed point inside a stable Hopf circle, between which an unstable Hopf circle resides. The interest of this solution structure is that the asymptotic behavior of the system can be switched between stationary and quasi-periodic motions by only tuning the initial state conditions. A set of critical and stability conditions for the degenerate Hopf bifurcation are discussed. The washout-filter-based controller with a polynomial control law is utilized. The control gains are derived from the theory of Chenciner's degenerate Hopf bifurcation with the aid of the center manifold reduction and the normal form evolution.     

Hypergeometric Equation in Modeling Relativistic Isotropic Sphere

We study the Einstein system of equations in static spherically symmetric spacetimes. We obtained classes of exact solutions to the Einstein system by transforming the condition for pressure isotropy to a hypergeometric equation choosing a rational form for one of the gravitational potentials. The solutions are given in simple form that is a desirable requisite to study the behavior of relativistic compact objects in detail. A physical analysis indicate that our models satisfy all the fundamental requirements of realistic star and match smoothly with the exterior Schwarzschild metric. The derived masses and densities are consistent with the previously reported experimental and theoretical studies describing strange stars. The models satisfy the standard energy conditions required by normal matter.     

Uphill anomalous transport in a deterministic system with speed-dependent friction coefficient

We investigate the transport of a deterministic Brownian particle theoretically, which moves in simple onedimensional, symmetric periodic potentials under the influence of both a time periodic and a static biasing force. The physical system employed contains a friction coefficient that is speed-dependent. Within the tailored parameter regime, the absolute negative mobility, in which a particle can travel in the direction opposite to a constant applied force, is observed.This behavior is robust and can be maximized at two regimes upon variation of the characteristic factor of friction coefficient. Further analysis reveals that this uphill motion is subdiffusion in terms of localization(diffusion coefficient with the form D(t) ~t~(-1) at long times). We also have observed the non-trivially anomalous subdiffusion which is significantly deviated from the localization; whereas most of the downhill motion evolves chaotically, with the normal diffusion.     

A perturbation method for computing the simplest normal forms of dynamical systems

A previously developed perturbation method is generalized for computing the simplest normal form (at each level of computation, the minimum number of terms are retained) of general n-dimensional differential equations. This “direct” approach, combining the normal form theory with center manifold theory in one unified procedure, can be used to systematically compute the simplest (or unique) normal form. Two particular singularities of the Jacobian of the system are considered in this paper: the first one is associated with one pair of purely imaginary eigenvalues (Hopf-type singularity), and the other corresponds to a simple zero and a pair of purely imaginary eigenvalues (Hopf-zero-type singularity). The approach can be easily formulated and implemented using a computer algebra system. Maple programs have been developed in this paper which can be “automatically” executed by a user without the knowledge of computer algebra. A physical oscillator model is studied in detail to show the computational efficiency of the “direct” method, and the advantage of using the simplest normal form, which greatly simplifies the analysis on dynamical systems, in particular, for bifurcations and stability.     

Aharonov-Bohm shift as a shift in normal coordinates

The Aharonov-Bohm shift in a closed system is considered. The solenoid is a charged, rotating cylinder which is electrically neutral. This model of Henneberger and Opatrny has a Hamiltonian which is a quadratic form. This quadratic form is transformed to normal coordinates, so that the stationary states become self-evident. It is shown that, in the original system, it is the kinetic angular momentum which is quantized. Solutions of the problem for an electron inside the solenoid are discussed. It is shown that the rotating cylinder exhibits different behavior if the electron is in the magnetic field or if it is in the external region. An external field approximation which replaces the cylinder by a constant magnetic field therefore cannot yield a correct solution of the Schrödinger equation which is continuous at the surface of the solenoid.     

The critical behavior of Kac-type models for semi-infinite systems

The critical behavior of the layer magnetizations and local susceptibilities of theD-vector lattice models with Kac-type ferromagnetic interactions for a semi-infinite system is studied. These local quantities behave less singularly than the bulk ones, showing that this is not peculiar to the two-dimensional Ising model. Moreover, the limiting form (at the critical point) of the magnetization profile can be obtained, which, when properly scaled, satisfies the minimum condition in the Landau theory for a semi-infinite continuous system. Landau-type critical behavior is thus recovered.     

The numerical method for three-dimensional impact with friction of multi-rigid-body system

The differential equations for planar impacts reduce to an algebraic form, and can be easily solved. For three dimensional impacts with friction, there is no closed-form solution, and numerical integration is required due to the swerve behavior of tangential impulse during collisions. The dynamic governing equations in the impact process are built up in impulse space based on the Lagrangian equation in this paper. The coefficient of restitution defined by Poisson is used as the condition of impact termination. A valid numerical method for solving three-dimensional frictional impact of multi-rigid body system is established. The singular cases of tangential movement in sticking point are especially noticed and analyzed. Several examples are present to reveal the different kinds of tangential movement modes varied with the normal impulse during collision.     

Circuit QED and sudden phase switching in a superconducting qubit array

Superconducting qubits connected in an array can form quantum many-body systems such as the quantum Ising model. By coupling the qubits to a superconducting resonator, the combined system forms a circuit QED system. Here, we study the nonlinear behavior in the many-body state of the qubit array using a semiclassical approach. We show that sudden switchings as well as a bistable regime between the ferromagnetic phase and the paramagnetic phase can be observed in the qubit array. A superconducting circuit to implement this system is presented with realistic parameters.     

HIGHER ORDER NORMAL FORM AND PERIOD AVERAGING

Calculation of the higher order averaging equations of a non-linear oscillator is very tedious using the classical averaging method. This is also true for higher order normal forms. This paper presents an alternative method, which is a combination of the method of normal form and the classical averaging method. A simple and efficient program is given to calculate the higher order averaging equations by using the symbolic computer algebra system Mathematica. Furthermore, the program can be used to calculate the higher order coefficients of normal form. Four examples are given and compare well with the existing results.     

Model for the extreme wetting hysteresis of liquid helium on cesium

Although the liquid- 4He-cesium system is a nearly ideal one for studying wetting phenomena, it can show extreme hysteresis which is profoundly nonideal in behavior. It is suggested that this is due to the roughness of these Cs surfaces. We show that stable micropuddles of liquid 4He can form in shallow cavities on a Cs surface. It is the potential to form micropuddles, as the liquid tries to recede, which pins the contact line due to the large energy needed to create the surface of a micropuddle. This model also accounts for the memory that these surfaces have of being in contact with liquid 4He.