Refine
Has Fulltext
- no (182)
Year of publication
Document Type
- Article (182) (remove)
Language
- English (182) (remove)
Is part of the Bibliography
- yes (182)
Keywords
- Kuramoto model (6)
- synchronization (4)
- Coupled oscillators (3)
- Lyapunov exponent (3)
- Hilbert transform (2)
- Network reconstruction (2)
- Nonlinear coupling (2)
- Oscillator ensembles (2)
- Synchronization (2)
- chaos (2)
Institute
We study the scaling properties of energy spreading in disordered strongly nonlinear Hamiltonian lattices. Such lattices consist of nonlinearly coupled local linear or nonlinear oscillators, and demonstrate a rather slow, subdiffusive spreading of initially localized wave packets. We use a fractional nonlinear diffusion equation as a heuristic model of this process, and confirm that the scaling predictions resulting from a self-similar solution of this equation are indeed applicable to all studied cases. We show that the spreading in nonlinearly coupled linear oscillators slows down compared to a pure power law, while for nonlinear local oscillators a power law is valid in the whole studied range of parameters.
We study populations of globally coupled noisy rotators (oscillators with inertia) allowing a nonequilibrium transition from a desynchronized state to a synchronous one (with the nonvanishing order parameter). The newly developed analytical approaches resulted in solutions describing the synchronous state with constant order parameter for weakly inertial rotators, including the case of zero inertia, when the model is reduced to the Kuramoto model of coupled noise oscillators. These approaches provide also analytical criteria distinguishing supercritical and subcritical transitions to the desynchronized state and indicate the universality of such transitions in rotator ensembles. All the obtained analytical results are confirmed by the numerical ones, both by direct simulations of the large ensembles and by solution of the associated Fokker-Planck equation. We also propose generalizations of the developed approaches for setups where different rotators parameters (natural frequencies, masses, noise intensities, strengths and phase shifts in coupling) are dispersed.
We study populations of globally coupled noisy rotators (oscillators with inertia) allowing a nonequilibrium transition from a desynchronized state to a synchronous one (with the nonvanishing order parameter). The newly developed analytical approaches resulted in solutions describing the synchronous state with constant order parameter for weakly inertial rotators, including the case of zero inertia, when the model is reduced to the Kuramoto model of coupled noise oscillators. These approaches provide also analytical criteria distinguishing supercritical and subcritical transitions to the desynchronized state and indicate the universality of such transitions in rotator ensembles. All the obtained analytical results are confirmed by the numerical ones, both by direct simulations of the large ensembles and by solution of the associated Fokker-Planck equation. We also propose generalizations of the developed approaches for setups where different rotators parameters (natural frequencies, masses, noise intensities, strengths and phase shifts in coupling) are dispersed.
We study an array of activity rhythms generated by a half-center oscillator (HCO), represented by a pair of reciprocally coupled neurons with post-inhibitory rebounds (PIR). Such coupling induced bursting possesses two time scales, one for fast spiking and another for slow quiescent periods, is shown to exhibit an array of synchronization properties. We discuss several HCO configurations constituted by two endogenous bursters, by tonic-spiking and quiescent neurons, as well as mixed-mode configurations composed of neurons of different type. We demonstrate that burst synchronization can be accompanied by complex, often chaotic, interactions of fast spikes within synchronized bursts. (C) 2015 Elsevier B.V. All rights reserved.
We demonstrate the existence of phase synchronization of two chaotic rotators. Contrary to phase synchronization of chaotic oscillators, here the Lyapunov exponents corresponding to both phases remain positive even in the synchronous regime. Such frequency locked dynamics with different ratios of frequencies are studied for driven continuous-time rotators and for discrete circle maps. We show that this transition to phase synchronization occurs via a crisis transition to a band-structured attractor.
Attractor-repeller collision and eyelet intermittency at the transition to phase synchronization
(1997)
The chaotically driven circle map is considered as the simplest model ofphase synchronization of a chaotic continuous-time oscillator by external periodic force. The phase dynamics is analyzed via phase-locking regions of the periodic cycles embedded in the strange attractor. It is shown that full synchronization, where all the periodic cycles are phase locked, disappears via the attractor-repeller collision. Beyond the transition an intermittent regime with exponentially rare phase slips, resulting from the trajectory's hits on an eyelet, is observed.
We demonstrate that a multiple delayed feedback is a powerful tool to control coherence properties of autonomous self-sustained oscillators. We derive the equation for the phase dynamics in presence of noise and delay, and analyze it analytically. In Gaussian approximation a closed set of equations for the frequency and the diffusion constant is obtained. Solutions of these equations are in good agreement with direct numerical simulations.
Supercritical Kuramoto oscillators with distributed frequencies can be separated into two disjoint groups: an ordered one locked to the mean field, and a disordered one consisting of effectively decoupled oscillators-at least so in the thermodynamic limit. In finite ensembles, in contrast, such clear separation fails: The mean field fluctuates due to finite-size effects and thereby induces order in the disordered group. This publication demonstrates this effect, similar to noise-induced synchronization, in a purely deterministic system. We start by modeling the situation as a stationary mean field with additional white noise acting on a pair of unlocked Kuramoto oscillators. An analytical expression shows that the cross-correlation between the two increases with decreasing ratio of natural frequency difference and noise intensity. In a deterministic finite Kuramoto model, the strength of the mean-field fluctuations is inextricably linked to the typical natural frequency difference. Therefore, we let a fluctuating mean field, generated by a finite ensemble of active oscillators, act on pairs of passive oscillators with a microscopic natural frequency difference between which we then measure the cross-correlation, at both super- and subcritical coupling.
We present an alternative approach to finite-size effects around the synchronization transition in the standard Kuramoto model. Our main focus lies on the conditions under which a collective oscillatory mode is well defined. For this purpose, the minimal value of the amplitude of the complex Kuramoto order parameter appears as a proper indicator. The dependence of this minimum on coupling strength varies due to sampling variations and correlates with the sample kurtosis of the natural frequency distribution. The skewness of the frequency sample determines the frequency of the resulting collective mode. The effects of kurtosis and skewness hold in the thermodynamic limit of infinite ensembles. We prove this by integrating a self-consistency equation for the complex Kuramoto order parameter for two families of distributions with controlled kurtosis and skewness, respectively.
We study synchronization properties of three nonlinearly coupled chaotic maps. Coupling is introduced in such a way, that it cannot be reduced to pairwise terms, but includes combined action of all interacting units. For two models of nonlinear coupling we characterize the transition to complete synchrony, as well as partially synchronized states. Relation to hypernetworks of chaotic units is also discussed.
Randomly coupled neural fields demonstrate irregular variation of firing rates, if the coupling is strong enough, as has been shown by Sompolinsky et al. [Phys. Rev. Lett. 61, 259 (1988)]. We present a method for reconstruction of the coupling matrix from a time series of irregular firing rates. The approach is based on the particular property of the nonlinearity in the coupling, as the latter is determined by a sigmoidal gain function. We demonstrate that for a large enough data set and a small measurement noise, the method gives an accurate estimation of the coupling matrix and of other parameters of the system, including the gain function.
I study deterministic dynamics of chiral active particles in two dimensions. Particles are considered as discs interacting with elastic repulsive forces. An ensemble of particles, started from random initial conditions, demonstrates chaotic collisions resulting in their normal diffusion. This chaos is transient, as rather abruptly a synchronous collisionless state establishes. The life time of chaos grows exponentially with the number of particles. External forcing (periodic or chaotic) is shown to facilitate the synchronization transition.
We consider a social-type network of coupled phase oscillators. Such a network consists of an active core of mutually interacting elements, and of a flock of passive units, which follow the driving from the active elements, but otherwise are not interacting. We consider a ring geometry with a long-range coupling, where active oscillators form a fluctuating chimera pattern. We show that the passive elements are strongly correlated. This is explained by negative transversal Lyapunov exponents.
We study dynamics of a serial array of spin-torque oscillators with a parallel inductor-capacitor-resistor (LC R) load. In a large range of parameters the fully synchronous regime, where all the oscillators have the same state and the output field is maximal, is shown to be stable. However, not always such a robust complete synchronization develops from a random initial state; in many cases nontrivial clustering is observed, with a partial synchronization resulting in a quasiperiodic or chaotic mean-field dynamics.
We study numerically secondary modes on top of a chaotic state in disordered nonlinear lattices. Two basic models are considered, with or without a local on-site potential. By performing periodic spatial modulation of displacement and kinetic energy, and following the temporal evolution of the corresponding spatial profiles, we reveal different modes which can be interpreted as first and second sound.
We study how coherence of noisy oscillations can be optimally enhanced by external locking. Based on the condition of minimizing the phase diffusion constant, we find the optimal forcing explicitly in the limits of small and large noise, in dependence of the phase sensitivity of the oscillator. We show analytically that the form of the optimal force bifurcates with the noise intensity; this is confirmed by the analysis of an optimal locking forcing for an experimentally obtained phase sensitivity of a neural cell. In the limit of small noise, the results are compared with purely deterministic conditions of optimal locking.
We consider networks of coupled phase oscillators of different complexity: Kuramoto–Daido-type networks, generalized Winfree networks, and hypernetworks with triple interactions. For these setups an inverse problem of reconstruction of the network connections and of the coupling function from the observations of the phase dynamics is addressed. We show how a reconstruction based on the minimization of the squared error can be implemented in all these cases. Examples include random networks with full disorder both in the connections and in the coupling functions, as well as networks where the coupling functions are taken from experimental data of electrochemical oscillators. The method can be directly applied to asynchronous dynamics of units, while in the case of synchrony, additional phase resettings are necessary for reconstruction.
Topic and aim. Synchronization in populations of coupled oscillators can be characterized with order parameters that describe collective order in ensembles. A dependence of the order parameter on the coupling constants is well-known for coupled periodic oscillators. The goal of the study is to extend this analysis to ensembles of oscillators with chaotic phases, moreover with phases possessing hyperbolic chaos. Models and methods. Two models are studied in the paper. One is an abstract discrete-time map, composed with a hyperbolic Bernoulli transformation and with Kuramoto dynamics. Another model is a system of coupled continuous-time chaotic oscillators, where each individual oscillator has a hyperbolic attractor of Smale-Williams type. Results. The discrete-time model is studied with the Ott-Antonsen ansatz, which is shown to be invariant under the application of the Bernoulli map. The analysis of the resulting map for the order parameter shows, that the asynchronouis state is always stable, but the synchronous one becomes stable above a certain coupling strength. Numerical analysis of the continuous-time model reveals a complex sequence of transitions from an asynchronous state to a completely synchronous hyperbolic chaos, with intermediate stages that include regimes with periodic in time mean field, as well as with weakly and strongly irregular mean field variations. Discussion. Results demonstrate that synchronization of systems with hyperbolic chaos of phases is possible, although a rather strong coupling is required. The approach can be applied to other systems of interacting units with hyperbolic chaotic dynamics.
We consider the dynamics of the Kuramoto ensemble oscillators not included in a common synchronized cluster, where the mean field is subject to fluctuations. The fluctuations can be either related to the finite size of the ensemble or superimposed on the mean field in the form of common noise due to the constructive features of the system. It is shown that the states of such oscillators with close natural frequencies appear correlated with each other, since the mean-field fluctuations act as common noise. We quantify the effect with the synchronization index of two oscillators, which is calculated numerically and analytically as a function of the frequency difference and noise intensity. The results are rigorous for large ensembles with additional noise superimposed on the mean field and are qualitatively true for the systems where the mean-field fluctuations are due to the finite size of the ensemble. In the latter case, the effect is found to be independent of the number of oscillators in the ensemble.
We study the discrete nonlinear Schrodinger equation with a random potential in one dimension. It is characterized by the length, the strength of the random potential, and the field density that determines the effect of nonlinearity. Following the time evolution of the field and calculating the largest Lyapunov exponent, the probability of the system to be regular is established numerically and found to be a scaling function of the parameters. This property is used to calculate the asymptotic properties of the system in regimes beyond our computational power.
We explore ensemble inequivalence in long-range interacting systems by studying an XY model of classical spinswith ferromagnetic and nematic coupling. We demonstrate the inequivalence bymapping themicrocanonical phase diagram onto the canonical one, and also by doing the inverse mapping. We show that the equilibrium phase diagrams within the two ensembles strongly disagree within the regions of first-order transitions, exhibiting interesting features like temperature jumps. In particular, we discuss the coexistence and forbidden regions of different macroscopic states in both the phase diagrams.
We study the dynamics of the excitable Fitz Hugh-Nagumo system under external noisy driving. Noise activates the system producing a sequence of pulses. The coherence of these noise-induced oscillations is shown to be maximal for a certain noise amplitude. This new effect of coherence resonance is explained by different noise dependencies of the activation and the excursion times. A simple one-dimensional model based on the Langevin dynamics is proposed for the quantitative description of this phenomenon.
Comment on "Simple approach to the creation of a strange nonchaotic attractor in any chaotic system"
(2001)
We address the problem of existence of strange nonchaotic attractors (SNAs) in quasiperiodically forced dynamical systems. Recently, Shuai and Wong [Phys. Rev. E 59, 5338 (1999)] suggested a universal method for constructing a SNA in an arbitrary system possessing chaos. We demonstrate here that, in general, this method fails. For arbitrary systems, it gives a SNA only in a vicinity of transition to chaos. We discuss also a special example, where the method by Shuai and Wong indeed produces a SNA.
We study the dynamics of Lyapunov vectors in various models of one-dimensional distributed systems with spacetime chaos. We demonstrate that the vector corresponding to the maximum exponent is always localized and the localization region wanders irregularly. This localization is explained by interpreting the logarithm of the Lyapunov vector as a roughening interface. We show that for many systems, the `interface' belongs to the Kardar-Parisi- Zhang universality class. Accordingly, we discuss the scaling behaviour of finite-size effects and self-averaging properties of the Lyapunov exponents.
Clustering in ensembles of globally coupled identical chaotic oscillators is reconsidered using a twofold approach. Stability of clusters towards "emanation" of the elements is described with the evaporation Lyapunov exponents. It appears that direct numerical simulations of ensembles often lead to spurious clusters that have positive evaporation exponents, due to a numerical trap. We propose a numerical method that surmounts the spurious clustering. We also demonstrate that clustering can be very sensitive to the number of elements in the ensemble.
Phase compactons
(2006)
We study the phase dynamics of a chain of autonomous, self-sustained, dispersively coupled oscillators. In the quasicontinuum limit the basic discrete model reduces to a Korteveg-de Vries-like equation, but with a nonlinear dispersion. The system supports compactons - solitary waves with a compact support - and kovatons - compact formations of glued together kink-antikink pairs that propagate with a unique speed, but may assume an arbitrary width. We demonstrate that lattice solitary waves, though not exactly compact, have tails which decay at a superexponential rate. They are robust and collide nearly elastically and together with wave sources are the building blocks of the dynamics that emerges from typical initial conditions. In finite lattices, after a long time, the dynamics becomes chaotic. Numerical studies of the complex Ginzburg-Landau lattice show that the non-dispersive coupling causes a damping and deceleration, or growth and acceleration, of compactons. A simple perturbation method is applied to study these effects. (c) 2006 Elsevier B.V. All rights reserved
Lai et al. [Phys. Rev. E 62, R29 (2000)] claim that the angular velocity of the phase point moving along the chaotic trajectory in a properly chosen projection (the instantaneous frequency) is intermittent. Using the same examples, namely the Rössler and the Lorenz systems, we show the absence of intermittency in the dynamics of the instantaneous frequency.This is confirmed by demonstrating that the phase dynamics exhibits normal diffusion. We argue that the nonintermittent behavior is generic.
We propose a method for experimental detection of directionality of weak coupling between two self-sustained oscillators from bivariate data. The technique is applicable to both noisy and chaotic systems that can be nonidentical or even structurally different. We introduce an index that quantifies the asymmetry in coupling.
Self-organized partially synchronous dynamics in populations of nonlinearly coupled oscillators
(2009)
We analyze a minimal model of a population of identical oscillators with a nonlinear coupling-a generalization of the popular Kuramoto model. In addition to well-known for the Kuramoto model regimes of full synchrony, full asynchrony, and integrable neutral quasiperiodic states, ensembles of nonlinearly coupled oscillators demonstrate two novel nontrivial types of partially synchronized dynamics: self-organized bunch states and self-organized quasiperiodic dynamics. The analysis based on the Watanabe-Strogatz ansatz allows us to describe the self-organized bunch states in any finite ensemble as a set of equilibria, and the self-organized quasiperiodicity as a two-frequency quasiperiodic regime. An analytic solution in the thermodynamic limit of infinitely many oscillators is also discussed.
We consider general heterogeneous ensembles of phase oscillators, sine coupled to arbitrary external fields. Starting with the infinitely large ensembles, we extend the Watanabe-Strogatz theory, valid for identical oscillators, to cover the case of an arbitrary parameter distribution. The obtained equations yield the description of the ensemble dynamics in terms of collective variables and constants of motion. As a particular case of the general setup we consider hierarchically organized ensembles, consisting of a finite number of subpopulations, whereas the number of elements in a subpopulation can be both finite or infinite. Next, we link the Watanabe-Strogatz and Ott-Antonsen theories and demonstrate that the latter one corresponds to a particular choice of constants of motion. The approach is applied to the standard Kuramoto-Sakaguchi model, to its extension for the case of nonlinear coupling, and to the description of two interacting subpopulations, exhibiting a chimera state. With these examples we illustrate that, although the asymptotic dynamics can be found within the framework of the Ott-Antonsen theory, the transients depend on the constants of motion. The most dramatic effect is the dependence of the basins of attraction of different synchronous regimes on the initial configuration of phases.
In this paper, we discuss recent progress in research of ensembles of mean field coupled oscillators. Without an ambition to present a comprehensive review, we outline most interesting from our viewpoint results and surprises, as well as interrelations between different approaches. (c) 2015 AIP Publishing LLC.
We consider an ensemble of coupled nonlinear noisy oscillators demonstrating in the thermodynamic limit an Ising-type transition. In the ordered phase and for finite ensembles stochastic flips of the mean field are observed with the rate depending on the ensemble size. When a small periodic force acts on the ensemble, the linear response of the system has a maximum at a certain system size, similar to the stochastic resonance phenomenon. We demonstrate this effect of system size resonance for different types of noisy oscillators and for different ensemblesùlattices with nearest neighbors coupling and globally coupled populations. The Ising model is also shown to demonstrate the system size resonance.
Interplay of coupling and common noise at the transition to synchrony in oscillator populations
(2016)
There are two ways to synchronize oscillators: by coupling and by common forcing, which can be pure noise. By virtue of the Ott-Antonsen ansatz for sine-coupled phase oscillators, we obtain analytically tractable equations for the case where both coupling and common noise are present. While noise always tends to synchronize the phase oscillators, the repulsive coupling can act against synchrony, and we focus on this nontrivial situation. For identical oscillators, the fully synchronous state remains stable for small repulsive coupling; moreover it is an absorbing state which always wins over the asynchronous regime. For oscillators with a distribution of natural frequencies, we report on a counter-intuitive effect of dispersion (instead of usual convergence) of the oscillators frequencies at synchrony; the latter effect disappears if noise vanishes.
The connection between the macroscopic description of collective chaos and the underlying microscopic dynamics is thoroughly analysed in mean-field models of one-dimensional oscillators. We investigate to what extent infinitesimal perturbations of the microscopic configurations can provide information also on the stability of the corresponding macroscopic phase. In ensembles of identical one-dimensional dynamical units, it is possible to represent the microscopic configurations so as to make transparent their connection with the macroscopic world. As a result, we find evidence of an intermediate, mesoscopic, range of distances, over which the instability is neither controlled by the microscopic equations nor by the macroscopic ones. We examine a whole series of indicators, ranging from the usual microscopic Lyapunov exponents, to the collective ones, including finite-amplitude exponents. A system of pulse-coupled oscillators is also briefly reviewed as an example of non-identical phase oscillators where collective chaos spontaneously emerges.
Internal signals like one's heartbeats are centrally processed via specific pathways and both their neural representations as well as their conscious perception (interoception) provide key information for many cognitive processes. Recent empirical findings propose that neural processes in the insular cortex, which are related to bodily signals, might constitute a neurophysiological mechanism for the encoding of duration. Nevertheless, the exact nature of such a proposed relationship remains unclear. We aimed to address this question by searching for the effects of cardiac rhythm on time perception by the use of a duration reproduction paradigm. Time intervals used were of 0.5, 2, 3, 7, 10, 14, 25, and 40s length. In a framework of synchronization hypothesis, measures of phase locking between the cardiac cycle and start/stop signals of the reproduction task were calculated to quantify this relationship. The main result is that marginally significant synchronization indices (Sls) between the heart cycle and the time reproduction responses for the time intervals of 2, 3, 10, 14, and 25s length were obtained, while results were not significant for durations of 0.5, 7, and 40s length. On the single participant level, several subjects exhibited some synchrony between the heart cycle and the time reproduction responses, most pronounced for the time interval of 25s (8 out of 23 participants for 20% quantile). Better time reproduction accuracy was not related with larger degree of phase locking, but with greater vagal control of the heart. A higher interoceptive sensitivity (IS) was associated with a higher synchronization index (SI) for the 2s time interval only. We conclude that information obtained from the cardiac cycle is relevant for the encoding and reproduction of time in the time span of 2-25s. Sympathovagal tone as well as interoceptive processes mediate the accuracy of time estimation.
High-frequency (HF) deep brain stimulation (DBS) is the gold standard for the treatment of medically refractory movement disorders like Parkinson’s disease, essential tremor, and dystonia, with a significant potential for application to other neurological diseases. The standard setup of HF DBS utilizes an open-loop stimulation protocol, where a permanent HF electrical pulse train is administered to the brain target areas irrespectively of the ongoing neuronal dynamics. Recent experimental and clinical studies demonstrate that a closed-loop, adaptive DBS might be superior to the open-loop setup. We here combine the notion of the adaptive high-frequency stimulation approach, that aims at delivering stimulation adapted to the extent of appropriately detected biomarkers, with specifically desynchronizing stimulation protocols. To this end, we extend the delayed feedback stimulation methods, which are intrinsically closed-loop techniques and specifically designed to desynchronize abnormal neuronal synchronization, to pulsatile electrical brain stimulation. We show that permanent pulsatile high-frequency stimulation subjected to an amplitude modulation by linear or nonlinear delayed feedback methods can effectively and robustly desynchronize a STN-GPe network of model neurons and suggest this approach for desynchronizing closed-loop DBS.
The transition from fully synchronized behavior to two-cluster dynamics is investigated for a system of N globally coupled chaotic oscillators by means of a model of two coupled logistic maps. An uneven distribution of oscillators between the two clusters causes an asymmetry to arise in the coupling of the model system. While the transverse period-doubling bifurcation remains essentially unaffected by this asymmetry, the transverse pitchfork bifurcation is turned into a saddle-node bifurcation followed by a transcritical riddling bifurcation in which a periodic orbit embedded in the synchronized chaotic state loses its transverse stability. We show that the transcritical riddling transition is always hard. For this, we study the sequence of bifurcations that the asynchronous point cycles produced in the saddle-node bifurcation undergo, and show how the manifolds of these cycles control the magnitude of asynchronous bursts. In the case where the system involves two subpopulations of oscillators with a small mismatch of the parameters, the transcritical riddling will be replaced by two subsequent saddle-node bifurcations, or the saddle cycle involved in the transverse destabilization of the synchronized chaotic state may smoothly shift away from the synchronization manifold. In this way, the transcritical riddling bifurcation is substituted by a symmetry-breaking bifurcation, which is accompanied by the destruction of a thin invariant region around the symmetrical chaotic state.
In the present paper, we study phase waves of self-sustained oscillators with a nearest-neighbor dispersive coupling on an infinite lattice. To analyze the underlying dynamics, we approximate the lattice with a quasi-continuum (QC). The resulting partial differential model is then further reduced to the Gardner equation, which predicts many properties of the underlying solitary structures. Using an iterative procedure on the original lattice equations, we determine the shapes of solitary waves, kinks, and the flat-like solitons that we refer to as flatons. Direct numerical experiments reveal that the interaction of solitons and flatons on the lattice is notably clean. All in all, we find that both the QC and the Gardner equation predict remarkably well the discrete patterns and their dynamics.
We study the phase dynamics of a chain of autonomous oscillators with a dispersive coupling. In the quasicontinuum limit the basic discrete model reduces to a Korteveg-de Vries-like equation, but with a nonlinear dispersion. The system supports compactons: solitary waves with a compact support and kovatons which are compact formations of glued together kink-antikink pairs that may assume an arbitrary width. These robust objects seem to collide elastically and, together with wave trains, are the building blocks of the dynamics for typical initial conditions. Numerical studies of the complex Ginzburg-Landau and Van der Pol lattices show that the presence of a nondispersive coupling does not affect kovatons, but causes a damping and deceleration or growth and acceleration of compactons
We present and study a family of finite amplitude breathers on a genuinely anharmonic Klein-Gordon lattice embedded in a nonlinear site potential. The direct numerical simulations are supported by a quasilinear Schrodinger equation (QLS) derived by averaging out the fast oscillations assuming small, albeit finite, amplitude vibrations. The genuinely anharmonic interlattice forces induce breathers which are strongly localized with tails evanescing at a doubly exponential rate and are either close to a continuum, with discrete effects being suppressed, or close to an anticontinuum state, with discrete effects being enhanced. Whereas the D-QLS breathers appear to be always stable, in general there is a stability threshold which improves with spareness of the lattice.
We introduce and study a family of lattice equations which may be viewed either as a strongly nonlinear discrete extension of the Gardner equation, or a non-convex variant of the Lotka-Volterra chain. Their deceptively simple form supports a very rich family of complex solitary patterns. Some of these patterns are also found in the quasi-continuum rendition, but the more intriguing ones, like interlaced pairs of solitary waves, or waves which may reverse their direction either spontaneously or due a collision, are an intrinsic feature of the discrete realm.