Synchronization transitions induced by partial time delay in a excitatory–inhibitory coupled neuronal network

Information transmission delays are an inherent factor of neuronal systems as a consequence of the finite propagation speeds and time lapses occurring by both dendritic and synaptic processes. In real neuronal systems, some delay between two neurons is too small and can be ignored, which results in partial time delay. In this paper, we focus on investigating influences of partial time delay on synchronization transitions in a excitatory–inhibitory (E–I) coupled neuronal networks. Here, we suppose time delay between two neurons equals to \(\tau \) with probability \(p_{\mathrm{delay}}\) and investigate effect of partial time delay on synchronization transitions of the neuronal networks by controlling \(\tau \) and \(p_{\mathrm{delay}}\) under three cases. In these three cases, excitatory synapses are always considered to delayed with probability \(p_{\mathrm{delay}}\), while inhibitory synapses are considered to be without delays (case I), delayed with probability \(p_{\mathrm{delay}}\) (case II), and always delayed (case III), respectively. It is revealed that, in the first two cases, partial time delay has little influences on synchronization of the neuronal network for small \(p_{\mathrm{delay}}\), while it could induce synchronization transitions at \(\tau \) around integer multiples of the period of individual neuron T when \(p_{\mathrm{delay}}\) is large enough, while in the case III, partial time delay could induce synchronization transitions at \(\tau \) being around odd integer multiples of T / 2 for small \(p_{\mathrm{delay}}\) and at \(\tau \) being around integer multiples of T for large \(p_{\mathrm{delay}}\). Most interesting observation is that partial time delay could induce frequent synchronization transitions at \(\tau \) being around integer multiples of T / 2 for intermediate \(p_{\mathrm{delay}}\). Moreover, effect of rewiring probability on synchronization transitions induced by partial time delay has been discussed. It is found that synchronization transitions induced by partial time delay are robust to rewiring probability for large \(p_{\mathrm{delay}}\) under the three cases.     

A practical synchronization approach for fractional-order chaotic systems

A practical synchronization approach is proposed for a class of fractional-order chaotic systems to realize perfect \(\delta \)-synchronization, and the nonlinear functions in the fractional-order chaotic systems are all polynomials. The \(\delta \)-synchronization scheme in this paper means that the origin in synchronization error system is stable. The reliability of \(\delta \)-synchronization has been confirmed on a class of fractional-order chaotic systems with detailed theoretical proof and discussion. Furthermore, the \(\delta \)-synchronization scheme for the fractional-order Lorenz chaotic system and the fractional-order Chua circuit is presented to demonstrate the effectiveness of the proposed method.     

Exact solution and invariant for fractional Cattaneo anomalous diffusion of cells in two-dimensional comb framework

This paper presents an investigation on anomalous diffusion of cells in a two-dimensional comb framework with effects of fractional Cattaneo flux. Formulated governing equation is an evolution equation with the coexisting characteristics of parabolic (diffusion) and hyperbolic (wave) for \(\alpha \) in (0, 1). Exact solution is obtained by the special fractional integral transformations, and a novel invariant is established, i.e., \(\left\langle {x^{2}\left( t \right) } \right\rangle \cdot \left\langle P \right\rangle = 0.5\) (the mean square displacement multiplied by the total number of cells along the x-axis = 0.5). Moreover, the characteristics of cells distribution, the total number and the mean square displacement of cells along the x-axis with different involved parameters, especially with the fractional parameter evolution, are shown graphically and analyzed in detail. For the cells distribution versus x, it turns from parabolic and hyperbolic with the decrease in t or the increase in \(\alpha \) or \(\xi \). It is monotonically decreasing for the cells distribution versus \(\alpha \) with different x, t and \(\xi \). For the distribution versus t with different \(\alpha \) and \(\xi \) or versus \(\alpha \) with different t, it is monotonically decreasing for the distribution of total number while monotonically increasing for the distribution of mean square displacement. It is remarkable that the anomalous subdiffusion happens along the x-axis for arbitrary parameters which is different from the classical Cattaneo diffusion.     

On the Justification of Viscoelastic Elliptic Membrane Shell Equations

We consider a family of linearly viscoelastic shells with thickness \(2\varepsilon\), clamped along their entire lateral face, all having the same middle surface \(S=\boldsymbol{\theta}(\bar{\omega})\subset \mathbb{R}^{3}\), where \(\omega\subset\mathbb{R}^{2}\) is a bounded and connected open set with a Lipschitz-continuous boundary \(\gamma\). We make an essential geometrical assumption on the middle surface \(S\), which is satisfied if \(\gamma\) and \(\boldsymbol{\theta}\) are smooth enough and \(S\) is uniformly elliptic. We show that, if the applied body force density is \(O(1)\) with respect to \(\varepsilon\) and surface tractions density is \(O(\varepsilon)\), the solution of the scaled variational problem in curvilinear coordinates, \(\boldsymbol{u}( \varepsilon)\), defined over the fixed domain \(\varOmega=\omega\times (-1,1)\) for each \(t\in[0,T]\), converges to a limit \(\boldsymbol{u}\) with \(u_{\alpha}(\varepsilon)\rightarrow u_{\alpha}\) in \(W^{1,2}(0,T,H ^{1}(\varOmega))\) and \(u_{3}(\varepsilon)\rightarrow u_{3}\) in \(W^{1,2}(0,T,L^{2}(\varOmega))\) as \(\varepsilon\to0\). Moreover, we prove that this limit is independent of the transverse variable. Furthermore, the average \(\bar{\boldsymbol{u}}= \frac{1}{2}\int_{-1}^{1} \boldsymbol{u}dx_{3}\), which belongs to the space \(W^{1,2}(0,T, V_{M}( \omega))\), where

$$V_{M}(\omega)=H^{1}_{0}(\omega)\times H^{1}_{0}(\omega)\times L ^{2}(\omega), $$

satisfies what we have identified as (scaled) two-dimensional equations of a viscoelastic membrane elliptic shell, which includes a long-term memory that takes into account previous deformations. We finally provide convergence results which justify those equations.

     

Further Dense Properties of the Space of Circle Diffeomorphisms with a Liouville Rotation Number

In continuation of Matsumoto’s paper (Nonlinearity 25:1495–1511, 2012) we show that various subspaces are \(C^{\infty }\)-dense in the space of orientation-preserving \(C^{\infty }\)-diffeomorphisms of the circle with rotation number \(\alpha \), where \(\alpha \in {\mathbb {S}}^1\) is any prescribed Liouville number. In particular, for every odometer \({\mathcal {O}}\) of product type we prove the denseness of the subspace of diffeomorphisms which are orbit-equivalent to \({\mathcal {O}}\).     

Oriented Shadowing Property and $$\Omega $$-Stability for Vector Fields

We call that a vector field has the oriented shadowing property if for any \(\varepsilon >0\) there is \(d>0\) such that each \(d\)-pseudo orbit is \(\varepsilon \)-oriented shadowed by some real orbit. In this paper, we show that the \(C^1\)-interior of the set of vector fields with the oriented shadowing property is contained in the set of vector fields with the \(\Omega \)-stability.     

Continuity and Invariance of the Sacker–Sell Spectrum

The Sacker–Sell (also called dichotomy or dynamical) spectrum \(\varSigma \) is a fundamental concept in the geometric, as well as for a developing bifurcation theory of nonautonomous dynamical systems. In general, it behaves merely upper-semicontinuously and a perturbation theory is therefore delicate. This paper explores an operator-theoretical approach to obtain invariance and continuity conditions for both \(\varSigma \) and its dynamically relevant subsets. Our criteria allow to avoid nonautonomous bifurcations due to collapsing spectral intervals and justify numerical approximation schemes for \(\varSigma \).     

Initial and Boundary Blow-U Problem for $$$$-Lalacian Parabolic Equation with General Absortion

In this article, we investigate the initial and boundary blow-up problem for the \(p\)-Laplacian parabolic equation \(u_t-\Delta _p u=-b(x,t)f(u)\) over a smooth bounded domain \(\Omega \) of \(\mathbb {R}^N\) with \(N\ge 2\), where \(\Delta _pu=\mathrm{div}(|\nabla u|^{p-2}\nabla u)\) with \(p>1\), and \(f(u)\) is a function of regular variation at infinity. We study the existence and uniqueness of positive solutions, and their asymptotic behaviors near the parabolic boundary.     

A novel approach for generating multi-direction multi-double-scroll attractors

Attention is focused in this work on quasiperiodic motion of nonlinear systems whose spectrum contains uniformly spaced sideband frequencies with a distance \(\omega _{d}\) apart, around a frequency \(\omega \) with \(\omega \gg \omega _{d}\) and its integer multiples, which are referred to as carrier frequencies. The ratio of the two frequencies \(\omega \) and \(\omega _{d}\) is an irrational number. A new method based on the traditional incremental harmonic balance (IHB) method with multiple timescales, referred to as Lau method, where two timescales, \(\tau _{1}=\omega t\) (a fast timescale) and \(\tau _{2}=\omega _{d}t\) (a slow timescale), are introduced, is presented to analyze quasiperiodic motion of nonlinear systems. An amplitude increment algorithm is adapted to deal with cases where the two frequencies \(\omega \) and \(\omega _{d}\) are    unknown a priori, in order to automatically trace frequency response of quasiperiodic motion of nonlinear systems and accurately calculate all frequency components and their corresponding amplitudes. Results of application of the present IHB method to quasiperiodic free vibration of a hinged–clamped beam with internal resonance between two transverse modes are shown and compared with previously published results with Lau method and those from numerical integration. While differences are noted between results predicted by the present IHB method and Lau method, excellent agreement is achieved between results from the present IHB method and numerical integration even in cases of strongly nonlinear vibration. The present IHB method is also used to analyze quasiperiodic free vibration of high-dimensional models of the hinged–clamped beam.     

Spreading in a Shifting Environment Modeled by the Diffusive Logistic Equation with a Free Boundary

We investigate the influence of a shifting environment on the spreading of an invasive species through a model given by the diffusive logistic equation with a free boundary. When the environment is homogeneous and favourable, this model was first studied in Du and Lin (SIAM J Math Anal 42:377–405, 2010), where a spreading–vanishing dichotomy was established for the long-time dynamics of the species, and when spreading happens, it was shown that the species invades the new territory at some uniquely determined asymptotic speed \(c_0>0\). Here we consider the situation that part of such an environment becomes unfavourable, and the unfavourable range of the environment moves into the favourable part with speed \(c>0\). We prove that when \(c\ge c_0\), the species always dies out in the long-run, but when \(0<c<c_0\), the long-time behavior of the species is determined by a trichotomy described by (a) vanishing, (b) borderline spreading, or (c) spreading. If the initial population is written in the form \(u_0(x)=\sigma \phi (x)\) with \(\phi \) fixed and \(\sigma >0\) a parameter, then there exists \(\sigma _0>0\) such that vanishing happens when \(\sigma \in (0,\sigma _0)\), borderline spreading happens when \(\sigma =\sigma _0\), and spreading happens when \(\sigma >\sigma _0\).     

Threshold Behavior and Non-quasiconvergent Solutions with Localized Initial Data for Bistable Reaction–Diffusion Equations

We consider bounded solutions of the semilinear heat equation \(u_t=u_{xx}+f(u)\) on \(R\), where \(f\) is of the unbalanced bistable type. We examine the \(\omega \)-limit sets of bounded solutions with respect to the locally uniform convergence. Our goal is to show that even for solutions whose initial data vanish at \(x=\pm \infty \), the \(\omega \)-limit sets may contain functions which are not steady states. Previously, such examples were known for balanced bistable nonlinearities. The novelty of the present result is that it applies to a robust class of nonlinearities. Our proof is based on an analysis of threshold solutions for ordered families of initial data whose limits at infinity are not necessarily zeros of \(f\).     

A Note on the Convergence of Singularly Perturbed Second Order Potential-Type Equations

In this paper we study the limit as \(\varepsilon \rightarrow 0\) of the singularly perturbed second order equation \(\varepsilon ^2 \ddot{u}_\varepsilon + \nabla _{\!x} V(t,u_\varepsilon (t))=0\), where V(t, x) is a potential. We assume that \(u_0(t)\) is one of its equilibrium points such that \(\nabla _{\!x}V(t,u_0(t))=0\) and \(\nabla _{\!x}^2V(t,u_0(t))>0\). We find that, under suitable initial data, the solutions \(u_\varepsilon \) converge uniformly to \(u_0\), by imposing mild hypotheses on V. A counterexample shows that they cannot be weakened.     

Models of Elastic Shells in Contact with a Rigid Foundation: An Asymptotic Approach

We consider a family of linearly elastic shells with thickness \(2\varepsilon\) (where \(\varepsilon\) is a small parameter). The shells are clamped along a portion of their lateral face, all having the same middle surface \(S\), and may enter in contact with a rigid foundation along the bottom face.We are interested in studying the limit behavior of both the three-dimensional problems, given in curvilinear coordinates, and their solutions (displacements \(\boldsymbol{u}^{\varepsilon}\) of covariant components \(u_{i}^{\varepsilon}\)) when \(\varepsilon\) tends to zero. To do that, we use asymptotic analysis methods. On one hand, we find that if the applied body force density is \(O(1)\) with respect to \(\varepsilon\) and surface tractions density is \(O(\varepsilon)\), a suitable approximation of the variational formulation of the contact problem is a two-dimensional variational inequality which can be identified as the variational formulation of the obstacle problem for an elastic membrane. On the other hand, if the applied body force density is \(O(\varepsilon^{2})\) and surface tractions density is \(O(\varepsilon^{3})\), the corresponding approximation is a different two-dimensional inequality which can be identified as the variational formulation of the obstacle problem for an elastic flexural shell. We finally discuss the existence and uniqueness of solution for the limit two-dimensional variational problems found.     

L.A.S. and Negative Schwarzian Derivative Do Not Imply G.A.S. in Clark’s Equation

It has been conjectured that a locally attracting equilibrium of Clark’s equation \(x_{n+1}=\alpha x_n + (1-\alpha )h(x_{n-k})\) is also globally attracting whenever h is a unimodal or decreasing map with negative Schwarzian derivative. The main aim of this paper is to show that the conjecture is false when \(k\ge 3\). This is done by studying the (Neimark–Sacker) bifurcation at the parameter \(\alpha \) where the locally attracting equilibrium of the equation becomes unstable. Our results, on the other hand, reinforce the validity of the conjecture in the cases \(k=1,2\).     

The Mixed Convection Boundary Layer on a Vertical Melting Front in a Non-Darcian Porous Medium

The effect of surface melting on the dual solutions that can arise in the problem of the mixed convection boundary-layer flow past a vertical surface embedded in a non-Darcian porous medium is considered. The problem is described by M, melting parameter, \(\lambda \), mixed convection parameter, and \(\gamma \), the flow inertia coefficient, numerical results being obtained in terms of these three parameters. It is seen that the melting phenomenon reduces the heat transfer rate and enhances the boundary-layer separation at the solid–liquid interface. Asymptotic solutions for the forced convection, \(\lambda =0\), and free convection, large \(\lambda \), limits are derived.     

A semi-locally scaled eddy viscosity formulation for LES wall models and flows at high speeds

We show that the mean wall-shear stresses in wall-modeled large-eddy simulations (WMLES) of high-speed flows can be off by up to \(\approx 100\%\) with respect to a DNS benchmark when using the van-Driest-based damping function, i.e., the conventional damping function. Errors in the WMLES-predicted wall-shear stresses are often attributed to the so-called log-layer mismatch, which, albeit also an error in wall-shear stresses \(\tau _\mathrm{w}\), is an error of about \(15\%\). The larger error identified here cannot be removed using the previously developed remedies for the log-layer mismatch. This error may be removed by using the semi-local scaling, i.e., \(l_\nu =\mu /\sqrt{\rho \tau _\mathrm{w}}\), in the damping function, where \(\mu \) and \(\rho \) are the local mean dynamic viscosity and density, respectively.     

Convergence of Radially Symmetric Solutions for (<Emphasis Type="Italic">p</Emphasis>, <Emphasis Type="Italic">q</Emphasis>)-Laplacian Elliptic Equations with a Damping Term

This study considers the quasilinear elliptic equation with a damping term,

$$\begin{aligned} \text {div}(D(u)\nabla u) + \frac{k(|{\mathbf {x}}|)}{|{\mathbf {x}}|}\,{\mathbf {x}}\cdot (D(u)\nabla u) + \omega ^2\big (|u|^{p-2}u + |u|^{q-2}u\big ) = 0, \end{aligned}$$

where \({\mathbf {x}}\) is an N-dimensional vector in \(\big \{{\mathbf {x}} \in \mathbb {R}^N: |{\mathbf {x}}| \ge \alpha \big \}\) for some \(\alpha > 0\) and \(N \in {\mathbb {N}}\setminus \{1\}\); \(D(u) = |\nabla u|^{p-2} + |\nabla u|^{q-2}\) with \(1 < q \le p\); k is a nonnegative and locally integrable function on \([\alpha ,\infty )\); and \(\omega \) is a positive constant. A necessary and sufficient condition is given for all radially symmetric solutions to converge to zero as \(|{\mathbf {x}}|\rightarrow \infty \). Our necessary and sufficient condition is expressed by an improper integral related to the damping coefficient k. The case that k is a power function is explained in detail.

     

Long-time Dynamics of Resonant Weakly Nonlinear CGL Equations

Consider a weakly nonlinear CGL equation on the torus \(\mathbb {T}^d\):

$$\begin{aligned} u_t+i\Delta u=\epsilon [\mu (-1)^{m-1}\Delta ^{m} u+b|u|^{2p}u+ ic|u|^{2q}u]. \end{aligned}$$

(*)

Here \(u=u(t,x)\), \(x\in \mathbb {T}^d\), \(0<\epsilon <<1\), \(\mu \geqslant 0\), \(b,c\in \mathbb {R}\) and \(m,p,q\in \mathbb {N}\). Define \(I(u)=(I_{\mathbf {k}},\mathbf {k}\in \mathbb {Z}^d)\), where \(I_{\mathbf {k}}=v_{\mathbf {k}}\bar{v}_{\mathbf {k}}/2\) and \(v_{\mathbf {k}}\), \(\mathbf {k}\in \mathbb {Z}^d\), are the Fourier coefficients of the function \(u\) we give. Assume that the equation \((*)\) is well posed on time intervals of order \(\epsilon ^{-1}\) and its solutions have there a-priori bounds, independent of the small parameter. Let \(u(t,x)\) solve the equation \((*)\). If \(\epsilon \) is small enough, then for \(t\lesssim {\epsilon ^{-1}}\), the quantity \(I(u(t,x))\) can be well described by solutions of an effective equation:

$$\begin{aligned} u_t=\epsilon [\mu (-1)^{m-1}\Delta ^m u+ F(u)], \end{aligned}$$

where the term \(F(u)\) can be constructed through a kind of resonant averaging of the nonlinearity \(b|u|^{2p}+ ic|u|^{2q}u\).

     

A Shilnikov Phenomenon Due to State-Dependent Delay,by Means of the Fixed Point Index

The first part of this paper is a general approach towards chaotic dynamics for a continuous map \(f:X\supset M\rightarrow X\) which employs the fixed point index and continuation. The second part deals with the differential equation

$$\begin{aligned} x'(t)=-\alpha \,x(t-d_{{\varDelta }}(x_t)). \end{aligned}$$

with state-dependent delay. For a suitable parameter \(\alpha \) close to \(5\pi /2\) we construct a delay functional \(d_{{\varDelta }}\), constant near the origin, so that the previous equation has a homoclinic solution, \(h(t)\rightarrow 0\) as \(t\rightarrow \pm \infty \), with certain regularity properties of the linearization of the semiflow along the flowline \(t\mapsto h_t\). The third part applies the method from the beginning to a return map which describes solution behaviour close to the homoclinic loop, and yields the existence of chaotic motion.

     

New Derived from Anosov Diffeomorphisms with Pathological Center Foliation

In this paper we focused our study on derived from Anosov diffeomorphisms (DA diffeomorphisms ) of the torus \(\mathbb {T}^3,\) it is, an absolute partially hyperbolic diffeomorphism on \(\mathbb {T}^3\) homotopic to a linear Anosov automorphism of the \(\mathbb {T}^3.\) We can prove that if \(f: \mathbb {T}^3 \rightarrow \mathbb {T}^3 \) is a volume preserving DA diffeomorphism homotopic to a linear Anosov A,  such that the center Lyapunov exponent satisfies \(\lambda ^c_f(x) > \lambda ^c_A > 0,\) with x belongs to a positive volume set, then the center foliation of f is non absolutely continuous. We construct a new open class U of non Anosov and volume preserving DA diffeomorphisms, satisfying the property \(\lambda ^c_f(x) > \lambda ^c_A > 0\) for \(m-\)almost everywhere \(x \in \mathbb {T}^3.\) Particularly for every \(f \in U,\) the center foliation of f is non absolutely continuous.