Analysis of the T-point-Hopf bifurcation

In a parameterized three-dimensional system of autonomous differential equations, a T-point is a point of the parameter space where a special kind of codimension-2 heteroclinic cycle occurs. If the parameter space is three-dimensional, such a bifurcation is located generically on a curve. A more degenerate scenario appears when this curve reaches a surface of Hopf bifurcations of one of the equilibria involved in the heteroclinic cycle. We are interested in the analysis of this codimension-3 bifurcation, which we call T-point-Hopf. In this work we propose a model, based on the construction of a Poincaré map, that describes the global behavior close to a T-point-Hopf bifurcation. The existence of certain kinds of homoclinic and heteroclinic connections between equilibria and/or periodic orbits is proved. The predictions deduced from this model strongly agree with the numerical results obtained in a modified van der Pol-Duffing electronic oscillator.     

Complex oscillations and waves of calcium in pancreatic acinar cells

We perform a bifurcation analysis of a model of Ca²⁺ wave propagation in the basal region of pancreatic acinar cells. The model we consider was first presented in Sneyd et al. [J. Sneyd, K. Tsaneva-Atanasova, J.I.E. Bruce, S.V. Straub, D.R. Giovannucci, D.I. Yule, A model of calcium waves in pancreatic and parotid acinar cells, Biophys. J. 85 (2003) 1392–1405], where a partial bifurcation analysis was given of the model in the absence of diffusion. We obtain more complete information about bifurcations of the diffusionless model via numerical studies, then analyse the spatially extended model by numerical investigation of the travelling wave equations and direct numerical solution of the model equations. We find solitary waves in the model equations arising from homoclinic bifurcations in the travelling wave equations. The solitary waves exist and appear to be stable for a significant interval of the primary bifurcation parameter (i.e., the concentration of inositol trisphosphate) but are eventually replaced by irregular spatio-temporal behaviour. The homoclinic bifurcations are related to a number of complicated mathematical structures in the travelling wave equations, including an anomalous homoclinic-Hopf bifurcation, heteroclinic bifurcations between an equilibrium and a periodic orbit, and homoclinic bifurcations of periodic orbits.     

Basic structures of the Shilnikov homoclinic bifurcation scenario

We find numerically small scale basic structures of homoclinic bifurcation curves in the parameter space of the Chua circuit. The distribution of these basic structures in the parameter space and their geometrical properties constitute a complete homoclinic bifurcation scenario of this system. Furthermore, these structures and the scenario are theoretically demonstrated to be generic to a large class of dynamical systems that presents, as the Chua circuit, Shilnikov homoclinic orbits. We classify the complexity of primary and subsidiary homoclinic orbits by their order given by the number of their returning loops. Our results confirm previous predictions of structures of homoclinic bifurcation curves and extend this study to high order primary orbits. Furthermore, we identify accumulations of bifurcation curves of subsidiary homoclinic orbits into bifurcation curves of both primary and subsidiary orbits.     

T-points: A codimension two heteroclinic bifurcation

The local bifurcation structure of a heteroclinic bifurcation which has been observed in the Lorenz equations is analyzed. The existence of a particular heteroclinic loop at one point in a two-dimensional parameter space (a T point) implies the existence of a line of heteroclinic loops and a logarithmic spiral of homoclinic orbits, as well as countably many other topologically more complicatedT points in a small neighborhood in parameter space.     

谐和激励与有界噪声作用下具有同宿和异宿轨道的Duffing振子的混沌运动

研究了具有同宿轨道、异宿轨道的双势阱Duffing振子在谐和激励与有界噪声摄动下的混沌运动.基于同宿分叉和异宿分叉，由Melnikov理论推导了系统出现混沌运动的必要条件及出现分形边界的充分条件.结果表明：当Wiener过程的强度参数大于某一临界值时，噪声增大了诱发混沌运动的有界噪声的临界幅值，相应地缩小了参数空间的混沌域，且产生混沌运动的临界幅值随着噪声强度的增大而增大.同时数值计算了最大Lyapunov指数，由最大Lyapunov指数为零从另一角度得到了系统出现混沌运动的有界噪声的临界幅值，发现在Wi 关键词： 混沌 同宿和异宿分叉 随机Melnikov方法 最大Lyapunov指数     

Primary resonant optimal control for nested homoclinic and heteroclinic bifurcations in single-dof nonlinear oscillators

A unified control theorem is presented in this paper, whose aim is to suppress the transversal intersections of stable and unstable manifolds of homoclinic and heteroclinic orbits in the Poincarè map embedding in system dynamics. Based on the control theorem, a primary resonant optimal control technique (PROCT for short) is applied to a general single-dof nonlinear oscillator. The novelty of this technique is able to obtain the unified analytical expressions of the control gain and the control parameters for suppressing the homoclinic and heteroclinic bifurcations, where the control gain can guarantee that the control region where the homoclinic and heteroclinic bifurcations do not occur can be enlarged as much as possible at least cost. The technique is applied to a nonlinear oscillator with a pair of nested homoclinic and heteroclinic orbits. By the PROCT, the transversal intersections of homoclinic and heteroclinic orbits can be suppressed, respectively. The hopping phenomenon that there coexist two kinds of chaotic attractors of Duffing-type and pendulum-type can be suppressed. On the contrary, if the first amplitude coefficient is greater than the critical heteroclinic bifurcation value, then another degenerate hopping behavior of chaos will take place again. Therefore, the phenomenon of hopping is the dominant type of chaos in this oscillator, whose suppressing or inducing is admissible from the points of practical and theoretical view.     

Homoclinic (Heteroclinic) Orbit of Complex Dynamical System and Spiral Structure

Starting from iterated systems, it is shown that the homoclinic (heteroclinic) orbit is a kind of spiral structure. The emphasis is laid to show that there are homoclinic or heteroclinic orbits in complex discrete and continuous systems, and these homoclinic or heteroclinic orbits are some kind of spiral structure.     

Phase space structure and chaotic scattering in near-integrable systems

We investigate the bifurcation phenomena and the change in phase space structure connected with the transition from regular to chaotic scattering in classical systems with unbounded dynamics. The regular systems discussed in this paper are integrable ones in the sense of Liouville, possessing a degenerated unstable periodic orbit at infinity. By means of a McGehee transformation the degeneracy can be removed and the usual Melnikov method is applied to predict homoclinic crossings of stable and unstable manifolds for the perturbed system. The chosen examples are the perturbed radial Kepler problem and two kinetically coupled Morse oscillators with different potential parameters which model the stretching dynamics in ABC molecules. The calculated subharmonic and homoclinic Melnikov functions can be used to prove the existence of chaotic scattering and of elliptic and hyperbolic periodic orbits, to calculate the width of the main stochastic layer and of the resonances, and to predict the range of initial conditions where singularities in the scattering function are found. In the second example the value of the perturbation parameter at which channel transitions set in is calculated. The theoretical results are supplemented by numerical experiments.     

Slow passage through a transcritical bifurcation for Hamiltonian systems and the change in action due to a nonhyperbolic homoclinic orbit

One-degree of freedom conservative slowly varying Hamiltonian systems are analyzed in the case in which a saddle-center pair undergo a transcritical bifurcation. We analyze the case in which the method of averaging predicts the solution crosses the unperturbed homoclinic orbit at the precise time at which the transcritical bifurcation occurs. For the slow passage through the nonhyperbolic homoclinic orbit associated with a transcritical bifurcation, the solution consists of a large sequence of nonhyperbolic homoclinic orbits surrounded by autonomous nonlinear saddle approaches. The change in action is computed by matching these solutions to those obtained by averaging, valid before and after crossing the nonhyperbolic homoclinic orbit. For initial conditions near the stable manifold of the nonhyperbolic saddle point, one saddle approach has particularly small energy and instead satisfies a nonautonomous nonlinear equation, which provides a transition between nonhyperbolic homoclinic orbits, centers, and saddles. (c) 2000 American Institute of Physics.     

A variational approach to connecting orbits in nonlinear dynamical systems

We propose a variational method for determining homoclinic and heteroclinic orbits including spiral-shaped ones in nonlinear dynamical systems. Starting from a suitable initial curve, a homotopy evolution equation is used to approach a true connecting orbit. The procedure is an extension of a variational method that has been used previously for locating cycles, and avoids the need for linearization in search of simple connecting orbits. Examples of homoclinic and heteroclinic orbits for typical dynamical systems are presented. In particular, several heteroclinic orbits of the steady-state Kuramoto–Sivashinsky equation are found, which display interesting topological structures, closely related to those of the corresponding periodic orbits.     

Bifurcation to heteroclinic cycles and sensitivity in three and four coupled phase oscillators

We study the bifurcation and dynamical behaviour of the system of N globally coupled identical phase oscillators introduced by Hansel, Mato and Meunier, in the cases N=3 and N=4. This model has been found to exhibit robust ‘slow switching’ oscillations that are caused by the presence of robust heteroclinic attractors. This paper presents a bifurcation analysis of the system in an attempt to better understand the creation of such attractors. We consider bifurcations that occur in a system of identical oscillators on varying the parameters in the coupling function. These bifurcations preserve the permutation symmetry of the system. We then investigate the implications of these bifurcations for the sensitivity to detuning (i.e. the size of the smallest perturbations that give rise to loss of frequency locking).For N=3 we find three types of heteroclinic bifurcation that are codimension-one with symmetry. On varying two parameters in the coupling function we find three curves giving (a) an S₃-transcritical homoclinic bifurcation, (b) a saddle-node/heteroclinic bifurcation and (c) a Z₃-heteroclinic bifurcation. We also identify several global bifurcations with symmetry that organize the bifurcation diagram; these are codimension-two with symmetry.For N=4 oscillators we determine many (but not all) codimension-one bifurcations with symmetry, including those that lead to a robust heteroclinic cycle. A robust heteroclinic cycle is stable in an open region of parameter space and unstable in another open region. Furthermore, we verify that there is a subregion where the heteroclinic cycle is the only attractor of the system, while for other parts of the phase plane it can coexist with stable limit cycles. We finish with a discussion of bifurcations that appear for this coupling function and general N, as well as for more general coupling functions.     

Comment on "Sil'nikov chaos of the Liu system" [Chaos 18, 013113 (2008)

In the referenced paper, the authors use the undetermined coefficient method to prove analytically the existence of homoclinic and heteroclinic orbits in a Lorenz-like system. If the proof was correct, the existence of horseshoe chaos would be guaranteed via the Sil'nikov criterion. However, we hereby show that their demonstration is incorrect for two reasons. On the one hand, they wrongly use a symmetry the Lorenz-like system exhibits. On the other hand, they try to find structurally unstable global bifurcations by means of a series that is uniformly convergent in an open set of the parameter space: this would imply that the dynamical object they have found is structurally stable.     

Shilnikov homoclinic orbit bifurcations in the Chua's circuit

We analytically describe the complex scenario of homoclinic bifurcations in the Chua's circuit. We obtain a general scaling law that gives the ratio between bifurcation parameters of different nearby homoclinic orbits. As an application of this theoretical approach, we estimate the number of higher order subsidiary homoclinic orbits that appear between two consecutive lower order subsidiary orbits. Our analytical finds might be valid for a large class of dynamical systems and are numerically confirmed in the parameter space of the Chua's circuit.     

连续时间系统同宿轨的搜索算法及其应用

同宿轨的求解是非线性系统领域的核心问题之一, 特别是对动力系统分岔与混沌的研究有重要意义. 根据同宿轨的几何特点, 采用轨线逼近的方式, 通过定义逼近轨线与鞍点的距离, 将同宿轨的求解转化为求距离最小值的无约束非线性优化问题. 为了提高优化结果的完整性, 还提出了基于区间细分的搜索算法和实现方法, 并找出了Lorenz系统, Shimizu-Morioka系统和超混沌Lorenz系统等的多个同宿轨道和对应参数, 验证了本文方法的有效性. 关键词： 混沌 同宿轨 非线性系统 数值计算     

Bifurcations and chaotic threshold for a nonlinear system with an irrational restoring force

Nonlinear dynamical systems with an irrational restoring force often occur in both science and engineering, and always lead to a barrier for conventional nonlinear techniques. In this paper, we have investigated the global bifurcations and the chaos directly for a nonlinear system with irrational nonlinearity avoiding the conventional Taylor's expansion to retain the natural characteristics of the system. A series of transformations are proposed to convert the homoclinic orbits of the unperturbed system to the heteroclinic orbits in the new coordinate, which can be transformed back to the analytical expressions of the homoclinic orbits. Melnikov's method is employed to obtain the criteria for chaotic motion, which implies that the existence of homoclinic orbits to chaos arose from the breaking of homoclinic orbits under the perturbation of damping and external forcing. The efficiency of the criteria for chaotic motion obtained in this paper is verified via bifurcation diagrams, Lyapunov exponents, and numerical simulations. It is worthwhile noting that our study is an attempt to make a step toward the solution of the problem proposed by Cao Q J et al. (Cao Q J, Wiercigroch M, Pavlovskaia E E, Thompson J M T and Grebogi C 2008 Phil. Trans. R. Soc. A 366 635).     

Scarring by homoclinic and heteroclinic orbits

In addition to the well-known scarring effect of periodic orbits, we show here that homoclinic and heteroclinic orbits, which are cornerstones in the theory of classical chaos, also scar eigenfunctions of classically chaotic systems when associated closed circuits in phase space are properly quantized, thus introducing strong quantum correlations. The corresponding quantization rules are also established. This opens the door for developing computationally tractable methods to calculate eigenstates of chaotic systems.     

Subharmonic and homoclinic bifurcations in a parametrically forced pendulum

Depending on the parameters of a parametrically forced pendulum system the boundaries of subharmonic and homoclinic bifurcations are calculated on the basis of the Melnikov method and of averaging methods. It is shown that, as a parameter is varied, repeated resonances of successively higher periods occur culminating in homoclinic orbits. According to the theorem of Smale homoclinic bifurcation is the source of the unstable chaotic motions observed. For some selected parameter sets the theoretical predictions are tested by numerical calculations. Very good agreement is found between analytical and numerical results.     

Mixed-mode oscillations in a stochastic,piecewise-linear system

We analyse a piecewise-linear FitzHugh–Nagumo model. The system exhibits a canard near which both small amplitude and large amplitude periodic orbits exist. The addition of small noise induces mixed-mode oscillations (MMOs) in the vicinity of the canard point. We determine the effect of each model parameter on the stochastically driven MMOs. In particular we show that any parameter variation (such as a modification of the piecewise-linear function in the model) that leaves the ratio of noise amplitude to time-scale separation unchanged typically has little effect on the width of the interval of the primary bifurcation parameter over which MMOs occur. In that sense, the MMOs are robust. Furthermore, we show that the piecewise-linear model exhibits MMOs more readily than the classical FitzHugh–Nagumo model for which a cubic polynomial is the only nonlinearity. By studying a piecewise-linear model, we are able to explain results using analytical expressions and compare these with numerical investigations.     

Torus destruction via global bifurcations in a piecewise-smooth, continuous map with square-root nonlinearity

It has been shown recently that torus formation in piecewise-smooth maps can occur through a special type of border collision bifurcation in which a pair of complex conjugate Floquet multipliers “jump” from the inside to the outside of the unit circle. It has also been shown that a large class of impacting mechanical systems yield piecewise-smooth maps with square-root singularity. In this Letter we investigate the dynamics of a two-dimensional piecewise-smooth map with square-root type nonlinearity, and describe two new routes to chaos through the destruction of two-frequency torus. In the first scenario, we identify the transition to chaos through the destruction of a loop torus via homoclinic bifurcation. In the other scenario, a change of structure in the torus occurs via heteroclinic saddle connections. Further parameter changes lead to a homoclinic bifurcation resulting in the creation of a chaotic attractor. However, this scenario is much more complex, with the appearance of a sequence of heteroclinic and homoclinic bifurcations.