scispace - formally typeset
Search or ask a question

Showing papers in "arXiv: Chaotic Dynamics in 2000"


Posted Content
TL;DR: In this paper, the authors proposed a method to extract messages masked by a chaotic signal of a time-delay system with very high dimension and many positive Lyapunov exponents using a special embedding coordinate.
Abstract: We show how to extract messages masked by a chaotic signal of a time-delay system with very high dimension and many positive Lyapunov exponents. Using a special embedding coordinate, the infinite dimensional phase space of the time-delay system is projected to a special three-dimensional space, which enables us to identify the time-delay of the system from the transmitted signal, and reconstruct the chaotic dynamics to unmask the hidden message successfully. The message extraction procedure is illustrated by simulations with Mackey-Glass time-delay system for two type of masking schemes and different kinds of messages.

116 citations


Journal ArticleDOI
TL;DR: In this paper, the Fokker-Planck equation can be generalized into a fractional fokker planck equation, i.e. an equation which includes fractional space differentiations, in order to encompass the wide class of anomalous diffusions due to a Levy stable stochastic forcing.
Abstract: We demonstrate that the Fokker-Planck equation can be generalized into a 'Fractional Fokker-Planck' equation, i.e. an equation which includes fractional space differentiations, in order to encompass the wide class of anomalous diffusions due to a Levy stable stochastic forcing. A precise determination of this equation is obtained by substituting a Levy stable source to the classical gaussian one in the Langevin equation. This yields not only the anomalous diffusion coefficient, but a non trivial fractional operator which corresponds to the possible asymmetry of the Levy stable source. Both of them cannot be obtained by scaling arguments. The (mono-) scaling behaviors of the Fractional Fokker-Planck equation and of its solutions are analysed and a generalization of the Einstein relation for the anomalous diffusion coefficient is obtained. This generalization yields a straightforward physical interpretation of the parameters of Levy stable distributions. Furthermore, with the help of important examples, we show the applicability of the Fractional Fokker-Planck equation in physics.

102 citations


Posted Content
TL;DR: In this paper, the authors investigated Lagrangian relative dispersion in direct numerical simulation of two-dimensional inverse cascade turbulence, using both standard fixed time statistics and an exit time approach.
Abstract: We investigate Lagrangian relative dispersion in direct numerical simulation of two-dimensional inverse cascade turbulence. The analysis is performed by using both standard fixed time statistics and an exit time approach. Our results are in good agreement with the original Richardson's description in terms of diffusion equation. The deviations are only of the quantitative nature. These deviations, and the observed strong sensitivity to the finite size effects, are associated with the long-range-correlated nature of the particles' relative motion. The correlation, or persistence, parameter is measured by means of a Lagrangian ``turning point'' statistics.

69 citations


Posted Content
TL;DR: In this article, a two dimensional time-dependent Duffing oscillator model of macroscopic neocortex exhibits chaos for some ranges of parameters, and they embed this model in moderate noise, typical of the context presented in real neocortex, using PATHINT, a non-Monte Carlo pathintegral algorithm that is particularly adept in handling nonlinear Fokker-Planck systems.
Abstract: A two dimensional time-dependent Duffing oscillator model of macroscopic neocortex exhibits chaos for some ranges of parameters. We embed this model in moderate noise, typical of the context presented in real neocortex, using PATHINT, a non-Monte-Carlo path-integral algorithm that is particularly adept in handling nonlinear Fokker-Planck systems. This approach shows promise to investigate whether chaos in neocortex, as predicted by such models, can survive in noisy contexts.

63 citations


Journal ArticleDOI
TL;DR: In this paper, a diagonal summation over correlated periodic orbit pair contributions to the form factor for low-dimensional cases is performed for arbitrary, directed graphs with unitary transfer matrices and the conditions for periodic orbit pairs to be correlated (including correlations due to the unitarity of the transfer matrix) are given explicitly.
Abstract: Quantum graphs have recently been introduced as model systems to study the spectral statistics of linear wave problems with chaotic classical limits. It is proposed here to generalise this approach by considering arbitrary, directed graphs with unitary transfer matrices. An exponentially increasing contribution to the form factor is identified when performing a diagonal summation over periodic orbit degeneracy classes. A special class of graphs, so-called binary graphs, is studied in more detail. For these, the conditions for periodic orbit pairs to be correlated (including correlations due to the unitarity of the transfer matrix) can be given explicitly. Using combinatorial techniques it is possible to perform the summation over correlated periodic orbit pair contributions to the form factor for some low--dimensional cases. Gradual convergence towards random matrix results is observed when increasing the number of vertices of the binary graphs.

56 citations


Posted Content
TL;DR: In this article, an infinite dimensional family of exact solutions of the incompressible three-dimensional Euler equations is presented. These solutions, proposed by Gibbon and Ohkitani, have infinite kinetic energy and blow up in finite time.
Abstract: We present an infinite dimensional family of of exact solutions of the incompressible three-dimensional Euler equations. These solutions, proposed by Gibbon and Ohkitani, have infinite kinetic energy and blow up in finite time.

50 citations


Posted Content
TL;DR: In this article, a new and alternative algorithm for noise reduction using the methods of discrete wavelet transform and numerical differentiation of the data is presented, which is seen to work well for a wide range of noise strengths, even as large as 10% of the signal level.
Abstract: We have presented a new and alternative algorithm for noise reduction using the methods of discrete wavelet transform and numerical differentiation of the data. In our method the threshold for reducing noise comes out automatically. The algorithm has been applied to three model flow systems - Lorenz, Autocatalator, and Rossler systems - all evolving chaotically. The method is seen to work well for a wide range of noise strengths, even as large as 10% of the signal level. We have also applied the method successfully to noisy time series data obtained from the measurement of pressure fluctuations in a fluidized bed, and also to that obtained by conductivity measurement in a liquid surfactant experiment. In all the illustrations we have been able to observe that there is a clean separation in the frequencies covered by the differentiated signal and white noise.

49 citations


Journal ArticleDOI
TL;DR: In this paper, the authors analyze statistical properties of complex eigenvalues of generic sub-unitary random matrices with broken time-reversal invariance and show that the correlation functions acquire the universal form found earlier for weakly non-Hermitian non-random matrices.
Abstract: Random contractions (sub-unitary random matrices) appear naturally when considering quantized chaotic maps within a general theory of open linear stationary systems with discrete time. We analyze statistical properties of complex eigenvalues of generic $N\times N$ random matrices $\hat{A}$ of such a type, corresponding to systems with broken time-reversal invariance. Deviations from unitarity are characterized by rank $M\le N$ and a set of eigenvalues $0 >M,n$ the correlation functions acquire the universal form found earlier for weakly non-Hermitian random matrices.

48 citations


Journal ArticleDOI
TL;DR: In this paper, a spectral analysis of the propagator in a suitable Hilbert space (and Lie algebra) of quantum observables in Heisenberg picture is discussed as an alternative approach to characterize infinite temperature dynamics of non-linear quantum many-body systems or quantum fields.
Abstract: Spectral analysis of the {\em adjoint} propagator in a suitable Hilbert space (and Lie algebra) of quantum observables in Heisenberg picture is discussed as an alternative approach to characterize infinite temperature dynamics of non-linear quantum many-body systems or quantum fields, and to provide a bridge between ergodic properties of such systems and the results of classical ergodic theory. We begin by reviewing some recent analytic and numerical results along this lines. In some cases the Heisenberg dynamics inside the subalgebra of the relevant quantum observables can be mapped explicitly into the (conceptually much simpler) Schr\" odinger dynamics of a single one-(or few)-dimensional quantum particle. The main body of the paper is concerned with an application of the proposed method in order to work out explicitly the general spectral measures and the time correlation functions in {\em a quantum Ising spin 1/2 chain in a periodically kicking transversal magnetic field}, including the results for the simpler autonomous case of a static magnetic field in the appropriate limit. The main result, being a consequence of a purely continuous non-trivial part of the spectrum, is that the general time-correlation functions decay to their saturation values as $t^{-3/2}$.

47 citations


Journal ArticleDOI
TL;DR: In this paper, a statistically stationary turbulence with a mean shear gradient is realized in a flow driven by suitable body forces, where the flow domain is periodic in downstream and spanwise directions and bounded by stress free surfaces in the normal direction.
Abstract: A statistically stationary turbulence with a mean shear gradient is realized in a flow driven by suitable body forces. The flow domain is periodic in downstream and spanwise directions and bounded by stress free surfaces in the normal direction. Except for small layers near the surfaces the flow is homogeneous. The fluctuations in turbulent energy are less violent than in the simulations using remeshing, but the anisotropy on small scales as measured by the skewness of derivatives is similar and decays weakly with increasing Reynolds number.

43 citations


Posted Content
TL;DR: In this article, the connections between microscopic chaos, defined on a dynamical level and arising from collisions between molecules, and diffusion, characterized by a mean square displacement proportional to the time, were investigated.
Abstract: We investigate the connections between microscopic chaos, defined on a dynamical level and arising from collisions between molecules, and diffusion, characterized by a mean square displacement proportional to the time. We use a number of models involving a single particle moving in two dimensions and colliding with fixed scatterers. We find that a number of microscopically nonchaotic models exhibit diffusion, and that the standard methods of chaotic time series analysis are ill suited to the problem of distinguishing between chaotic and nonchaotic microscopic dynamics. However, we show that periodic orbits play an important role in our models, in that their different properties in chaotic and nonchaotic systems can be used to distinguish such systems at the level of time series analysis, and in systems with absorbing boundaries. Our findings are relevant to experiments aimed at verifying the existence of chaoticity and related dynamical properties on a microscopic level in diffusive systems.

Journal ArticleDOI
TL;DR: In this paper, a hierarchy of many-parameter family of chaotic maps of the interval is generated with an invariant measure, by composition of chaotic map of reference, and the measure is used to calculate Kolmogorov-Sinai entropy of these maps analytically.
Abstract: We generate new hierarchy of many-parameter family of maps of the interval [0,1] with an invariant measure, by composition of the chaotic maps of reference [1]. Using the measure, we calculate Kolmogorov-Sinai entropy, or equivalently Lyapunov characteristic exponent, of these maps analytically, where the results thus obtained have been approved with numerical simulation. In contrary to the usual one-dimensional maps and similar to the maps of reference [1], these maps do not possess period doubling or period-n-tupling cascade bifurcation to chaos, but they have single fixed point attractor at certain region of parameters values, where they bifurcate directly to chaos without having period-n-tupling scenario exactly at these values of parameter whose Lyapunov characteristic exponent begins to be positive.

Journal ArticleDOI
TL;DR: In this paper, the authors propose a plasma model for spectral statistics displaying level repulsion without long-range spectral rigidity, i.e. statistics intermediate between random matrix and Poisson statistics similar to the ones found numerically at the critical point of the Anderson metal-insulator transition in disordered systems and in certain dynamical systems.
Abstract: We propose a plasma model for spectral statistics displaying level repulsion without long-range spectral rigidity, i.e. statistics intermediate between random matrix and Poisson statistics similar to the ones found numerically at the critical point of the Anderson metal-insulator transition in disordered systems and in certain dynamical systems. The model emerges from Dysons one-dimensional gas corresponding to the eigenvalue distribution of the classical random matrix ensembles by restricting the logarithmic pair interaction to a finite number $k$ of nearest neighbors. We calculate analytically the spacing distributions and the two-level statistics. In particular we show that the number variance has the asymptotic form $\Sigma^2(L)\sim\chi L$ for large $L$ and the nearest-neighbor distribution decreases exponentially when $s\to \infty$, $P(s)\sim\exp (-\Lambda s)$ with $\Lambda=1/\chi=k\beta+1$, where $\beta$ is the inverse temperature of the gas ($\beta=$1, 2 and 4 for the orthogonal, unitary and symplectic symmetry class respectively). In the simplest case of $k=\beta=1$, the model leads to the so-called Semi-Poisson statistics characterized by particular simple correlation functions e.g. $P(s)=4s\exp(-2s)$. Furthermore we investigate the spectral statistics of several pseudointegrable quantum billiards numerically and compare them to the Semi-Poisson statistics.

Posted Content
TL;DR: In this paper, the authors dealt with feedback control of fractional-order Chua's system with total order less than three which exhibit chaos as well as other nonlinear behavior and theory for control of chaotic systems using sampled data.
Abstract: This paper deals with feedback control of fractional-order Chua's system. The fractional-order Chua's system with total order less than three which exhibit chaos as well as other nonlinear behavior and theory for control of chaotic systems using sampled data are presented. Numerical experimental example is shown to verify the theoretical results.

Posted Content
TL;DR: In this article, the authors consider a simple reaction-advection-diffusion equation with ignition-type nonlinearity and discuss the following question: What kinds of velocity profiles are capable of quenching any given flame, provided the velocity's amplitude is adequately large?
Abstract: We consider a simple scalar reaction-advection-diffusion equation with ignition-type nonlinearity and discuss the following question: What kinds of velocity profiles are capable of quenching any given flame, provided the velocity's amplitude is adequately large? Even for shear flows, the answer turns out to be surprisingly subtle. If the velocity profile changes in space so that it is nowhere identically constant, (or if it is identically constant only in a region of small measure) then the flow can quench any initial data. But if the velocity profile is identically constant in a sizable region, then the ensuing flow is incapable of quenching large enough flames, no matter how much larger is the amplitude of this velocity. The constancy region must be wider across than a couple of laminar propagating front-widths. The proof uses a linear PDE associated to the nonlinear problem and quenching follows when the PDE is hypoelliptic. The techniques used allow the derivation of new, nearly optimal bounds on the speed of traveling wave solutions.

Journal ArticleDOI
TL;DR: In this article, the authors study stochastically forced semilinear parabolic PDEs of the Ginzburg-Landau type and show that the solutions are with high probability analytic in a strip around the real axis and give estimates on the width of that strip.
Abstract: We study stochastically forced semilinear parabolic PDE's of the Ginzburg-Landau type. The class of forcings considered are white noises in time and colored smooth noises in space. Existence of the dynamics in $L^\infty$, as well as existence of an invariant measure are proven. We also show that the solutions are with high probability analytic in a strip around the real axis and give estimates on the width of that strip.

Journal ArticleDOI
TL;DR: In this article, the authors extend the semiclassical theory of scarring of quantum eigenfunctions by classical periodic orbits to include situations where these orbits undergo generic bifurcations and show that the scars are semiclassically wider than those associated with isolated and unstable periodic orbits; moreover, their amplitude is at least as large, and in most cases larger.
Abstract: We extend the semiclassical theory of scarring of quantum eigenfunctions psi_{n}(q) by classical periodic orbits to include situations where these orbits undergo generic bifurcations. It is shown that |psi_{n}(q)|^{2}, averaged locally with respect to position q and the energy spectrum E_{n}, has structure around bifurcating periodic orbits with an amplitude and length-scale whose hbar-dependence is determined by the bifurcation in question. Specifically, the amplitude scales as hbar^{alpha} and the length-scale as hbar^{w}, and values of the scar exponents, alpha and w, are computed for a variety of generic bifurcations. In each case, the scars are semiclassically wider than those associated with isolated and unstable periodic orbits; moreover, their amplitude is at least as large, and in most cases larger. In this sense, bifurcations may be said to give rise to superscars. The competition between the contributions from different bifurcations to determine the moments of the averaged eigenfunction amplitude is analysed. We argue that there is a resulting universal hbar-scaling in the semiclassical asymptotics of these moments for irregular states in systems with a mixed phase-space dynamics. Finally, a number of these predictions are illustrated by numerical computations for a family of perturbed cat maps.

Posted Content
TL;DR: In this article, a signal-driven scheme of the chaotic neural network with the coupling constants corresponding to certain information was set up, and the stochastic resonance-like effects under its deterministic dynamics were investigated.
Abstract: We set up a signal-driven scheme of the chaotic neural network with the coupling constants corresponding to certain information, and investigate the stochastic resonance-like effects under its deterministic dynamics, comparing with the conventional case of Hopfield network with stochastic noise. It is shown that the chaotic neural network can enhance weak subthreshold signals and have higher coherence abilities between stimulus and response than those attained by the conventional stochastic model.

Journal ArticleDOI
TL;DR: In this paper, the authors present a comprehensive investigation of the effect of the exit time statistics on the entropy of a 3-dimensional velocity field in dynamical systems, stochastic processes and turbulence.
Abstract: We present a comprehensive investigation of $\epsilon$-entropy, $h(\epsilon)$, in dynamical systems, stochastic processes and turbulence. Particular emphasis is devoted on a recently proposed approach to the calculation of the $\epsilon$-entropy based on the exit-time statistics. The advantages of this method are demonstrated in examples of deterministic diffusive maps, intermittent maps, stochastic self-affine and multi-affine signals and experimental turbulent data. Concerning turbulence, the multifractal formalism applied to the exit time statistics allows us to predict that $h(\epsilon)\sim \epsilon^{-3}$ for velocity time measurement. This power law is independent of the presence of intermittency and has been confirmed by the experimental data analysis. Moreover, we show that the $\epsilon$-entropy density of a 3-dimensional velocity field is affected by the correlations induced by the sweeping of large scales.

Journal ArticleDOI
TL;DR: In this article, a hierarchy of one-parameter family F(a,x) of maps of the interval [0,1] with an invariant measure was given, where the Kolmogorov-Sinai entropy was calculated analytically.
Abstract: We give hierarchy of one-parameter family F(a,x) of maps of the interval [0,1] with an invariant measure. Using the measure, we calculate Kolmogorov-Sinai entropy, or equivalently Lyapunov characteristic exponent, of these maps analytically, where the results thus obtained have been approved with numerical simulation. In contrary to the usual one-parameter family of maps such as logistic and tent maps, these maps do not possess period doubling or period-n-tupling cascade bifurcation to chaos, but they have single fixed point attractor at certain parameter values, where they bifurcate directly to chaos without having period-n-tupling scenario exactly at these values of parameter whose Lyapunov characteristic exponent begins to be positive.

Journal ArticleDOI
TL;DR: A discrete-time model of reacting evolving fields, transported by a bidimensional chaotic fluid flow, is studied, based on the use of a Lagrangian scheme where fluid particles are advected by a two-dimensional symplectic map possibly yielding Lagrangia chaos.
Abstract: A discrete-time model of reacting evolving fields, transported by a bidimensional chaotic fluid flow, is studied. Our approach is based on the use of a Lagrangian scheme where {\it fluid particles} are advected by a $2d$ symplectic map possibly yielding Lagrangian chaos. Each {\it fluid particle} carries concentrations of active substances which evolve according to its own reaction dynamics. This evolution is also modeled in terms of maps. Motivated by the question, of relevance in marine ecology, of how a localized distribution of nutrients or preys affects the spatial structure of predators transported by a fluid flow, we study a specific model in which the population dynamics is given by a logistic map with space-dependent coefficient, and advection is given by the standard map. Fractal and random patterns in the Eulerian spatial concentration of predators are obtained under different conditions. Exploiting the analogies of this coupled-map (advection plus reaction) system with a random map, some features of these patterns are discussed.

Journal ArticleDOI
TL;DR: In this paper, the results of a numerical and semiclassical analysis of high lying states in a Hamiltonian system were presented, whose classical mechanics is of a generic, mixed type, where the energy surface is split into regions of regular and chaotic motion.
Abstract: In this work we present the results of a numerical and semiclassical analysis of high lying states in a Hamiltonian system, whose classical mechanics is of a generic, mixed type, where the energy surface is split into regions of regular and chaotic motion. As predicted by the principle of uniform semiclassical condensation (PUSC), when the effective $\hbar$ tends to 0, each state can be classified as regular or irregular. We were able to semiclassically reproduce individual regular states by the EBK torus quantization, for which we devise a new approach, while for the irregular ones we found the semiclassical prediction of their autocorrelation function, in a good agreement with numerics. We also looked at the low lying states to better understand the onset of semiclassical behaviour.

Journal ArticleDOI
TL;DR: In this article, the predictability problem in the inverse energy cascade of two-dimensional turbulence is addressed by means of direct numerical simulations, and the growth rate as a function of the error level is determined by a finite size extension of the Lyapunov exponent.
Abstract: The predictability problem in the inverse energy cascade of two-dimensional turbulence is addressed by means of direct numerical simulations. The growth rate as a function of the error level is determined by means of a finite size extension of the Lyapunov exponent. For error within the inertial range, the linear growth of the error energy, predicted by dimensional argument, is verified with great accuracy. Our numerical findings are in close agreement with the result of TFM closure approximation.

Journal ArticleDOI
TL;DR: In this paper, the starting semiclassical approximation to the quantum theory of a classical field can be described by means of an infinite hierarchy of its unstable spatio-temporally periodic solutions.
Abstract: Spatio-temporally chaotic dynamics of a classical field can be described by means of an infinite hierarchy of its unstable spatio-temporally periodic solutions. The periodic orbit theory yields the global averages characterizing the chaotic dynamics, as well as the starting semiclassical approximation to the quantum theory. New methods for computing corrections to the semiclassical approximation are developed; in particular, a nonlinear field transformation yields the perturbative corrections in a form more compact than the Feynman diagram expansions.

Journal ArticleDOI
TL;DR: In this article, an extension of the Poincare-Melnikov method was used to prove that for weak anisotropy, chaos shows up on the zero-energy manifold, and the existence of a large non-chaotic set of uniformly bounded and collisionless solutions.
Abstract: The anisotropic Manev problem, which lies at the intersection of classical, quantum, and relativity physics, describes the motion of two point masses in an anisotropic space under the influence of a Newtonian force-law with a relativistic correction term. Using an extension of the Poincare'-Melnikov method, we first prove that for weak anisotropy, chaos shows up on the zero-energy manifold. Then we put into the evidence a class of isolated periodic orbits and show that the system is nonintegrable. Finally, using the geodesic deviation approach, we prove the existence of a large non-chaotic set of uniformly bounded and collisionless solutions.

Journal ArticleDOI
Ayumu Sugita1
TL;DR: In this paper, the Maslov index of the periodic orbit is reduced to a classification of connections on a vector bundle over a structure group, and a Bohr-Sommerfeld condition including quadratic fluctuation around the orbit is derived using this formula.
Abstract: Semiclassical trace formulas are examined using phase space path integrals. Our main concern in this paper is the Maslov index of the periodic orbit, which seems not fully understood in previous works. We show that the calculation of the Maslov index is reduced to a classification of connections on a vector bundle over $S^{1}$ with structure group $Sp(2n,R)$. We derive a formula for the index of the n-repetition, and show that a Bohr-Sommerfeld type quantization condition including quadratic fluctuation around the orbit is derived using this formula.

Posted Content
TL;DR: In this article, a variational approach and dual maps techniques were used to derive exact results about several models of transport flows, for which they also obtained a complete description of their limit (in time) behavior in the space of configurations.
Abstract: By means of a novel variational approach and using dual maps techniques and general ideas of dynamical system theory we derive exact results about several models of transport flows, for which we also obtain a complete description of their limit (in time) behavior in the space of configurations. Using these results we study the motion of a speedy passive particle (tracer) moving along/against the flow of slow particles and demonstrate that the latter case might be more efficient.

Journal ArticleDOI
TL;DR: In this article, the high lying eigenstates of a mixed type Hamiltonian system were analyzed using the Einstein-Brillouin-Keller quantization, and a prediction for their wave function autocorrelation function was obtained.
Abstract: In this work we semiclassically analyzed the high lying eigenstates of a mixed type Hamiltonian system. For the regular states we employ the Einstein-Brillouin-Keller quantization, while for the chaotic states, following the principle of uniform semiclassical condensation, we obtain a prediction for their wavefunction autocorrelation function.

Book ChapterDOI
TL;DR: The kinetic theory of gases provides methods for calculating Lyapunov exponents and other quantities, such as Kolmogorov-Sinai entropies, that characterize the chaotic behavior of hard ball gases as discussed by the authors.
Abstract: The kinetic theory of gases provides methods for calculating Lyapunov exponents and other quantities, such as Kolmogorov-Sinai entropies, that characterize the chaotic behavior of hard ball gases. Here we illustrate the use of these methods for calculating the Kolmogorov-Sinai entropy, and the largest positive Lyapunov exponent, for dilute hard ball gases in equilibrium. The calculation of the largest Lyapunov exponent makes interesting connections with the theory of propagation of hydrodynamic fronts. Calculations are also presented for the Lyapunov spectrum of dilute, random Lorentz gases in two and three dimensions, which are considerably simpler than the corresponding calculations for hard ball gases. The article concludes with a brief discussion of some interesting open problems.

Journal ArticleDOI
TL;DR: A graphical analysis of the mechanisms underlying the occurrences of bubbling sequences and bistability regions in the bifurcation scenario of a special class of one dimensional two parameter maps concludes that whether it is bubbling or bistable is decided by the sign of the third derivative at the inflection point of the map function.
Abstract: We present a graphical analysis of the mechanisms underlying the occurrences of bubbling sequences and bistability regions in the bifurcation scenario of a special class of one dimensional two parameter maps. The main result of the analysis is that whether it is bubbling or bistability is decided by the sign of the third derivative at the inflection point of the map function.