Search

Connectivity is the safety network for biodiversity conservation because connected habitats are more effective for saving the species and ecological functions. The nature of coupling for connectivity also plays an important role in the co-existence of species in cyclic-dominance. The rock-paper-scissors game is one of the paradigmatic mathematical model in evolutionary game theory to understand the mechanism of biodiversity in cyclic-dominance. In this paper, the metapopulation model for rock-paper-scissors with mutations is presented in which the total population is divided into patches and the patches form a network of complete graph. The migration among patches is allowed through simple random walk. The replicator-mutator equations are used with the migration term. When migration is allowed then the population of the patches will synchronized and attain stable state through Hopf bifurcation. Apart form this, two phases are observed when the strategies of one of the species mutate to other two species: co-existence of all the species phase and existence of one kind of species phase. The transition from one phase to another phase is taking place due to transcritical bifurcation. The dynamics of the population of species of rock, paper, scissors is studied in the environment of homogeneous and heterogeneous mutation. Numerical simulations have been performed when mutation is allowed in all the patches (homogeneous mutation) and some of the patches (heterogeneous mutation). It has been observed that when the number of patches is increased in the case of heterogeneous mutation then the population of any of the species will not extinct and all the species will co-exist.

Classic frameworks of rock-paper-scissors game have been assumed in a closed community that a density of each group is only affected by internal factors such as competition interplay among groups and reproduction itself. In real systems in ecological and social sciences, however, the survival and a change of a density of a group can be also affected by various external factors. One of common features in real population systems in ecological and social sciences is population flow that is characterized by population inflow and outflow in a group or a society, which has been usually overlooked in previous works on models of rock-paper-scissors game. In this paper, we suggest the rock-paper-scissors system by implementing population flow and investigate its effect on biodiversity. For two scenarios of either balanced or imbalanced population flow, we found that the population flow can strongly affect group diversity by exhibiting rich phenomena. In particular, while the balanced flow can only lead the persistent coexistence of all groups which accompanies a phase transition through supercritical Hopf bifurcation on different carrying simplices, the imbalanced flow strongly facilitates rich dynamics such as alternative stable survival states by exhibiting various group survival states and multistability of sole group survivals by showing not fully covered but spirally entangled basins of initial densities due to local stabilities of associated fixed points. In addition, we found that, the system can exhibit oscillatory dynamics for coexistence by relativistic interplay of population flows which can capture the robustness of the coexistence state. Applying population flow in the rock-paper-scissors system can ultimately change a community paradigm from closed to open one, and our foundation can eventually reveal that population flow can be also a significant factor on a group density which is independent to fundamental interactions among groups.

We complete a gap in the proof that exceptional polynomials are complete orthogonal systems in the associated Hilbert spaces.

In this paper, it is shown that von Neumann continuous geometry may be regarded as the first attempt towards formulating a general quantum spacetime geometry akin to that of Cantorian spacetime E(∞) and noncommutative geometry.

This paper is dedicated to L p bounds on eigenfunctions of a Schrodinger-type operator ( − Δ g ) α / 2 + V on closed Riemannian manifolds for critically singular potentials V. The operator ( − Δ g ) α / 2 is defined spectrally in terms of the eigenfunctions of − Δ g . We obtain also quasimodes and spectral clusters estimates. As an application, we derive Strichartz estimates for the fractional wave equation ( ∂ t 2 + ( − Δ g ) α / 2 + V ) u = 0 . The wave kernel techniques recently developed by Bourgain-Shao-Sogge-Yao [4] and Shao-Yao [27] play a key role in this paper. We construct a new reproducing operator with several local operators and some good error terms. Moreover, we shall prove that these local operators satisfy certain variable coefficient versions of the “uniform Sobolev estimates” by Kenig-Ruiz-Sogge [18] . These enable us to handle the critically singular potentials V and prove the quasimode estimates.

International audience; In this paper, we consider a semi-linear fourth order parabolic equation in a bounded smooth domain Ω with homogeneous Dirichlet and Neumann boundary conditions. The main result of this paper is the null controllability and the exact controllability to the trajectories at any time T > 0 for the associated control system with a control function acting at the interior.; Dans ce papier, on considère uneéquation parabolique semi-linéaire de quatrième ordre dans un domaine borné régulier Ω avec des conditions aux limites de type Dirichlet et Neumann homogènes. Le résultat principal de ce papier concerne la contrôlabilitéà zéro et la contrôlabilité exacte pour tout T > 0 du système de contrôle associé avec un contrôle agissantà l'interieur.

The purpose of this paper is to show well-posedness results in L 2 -based Sobolev spaces for transmission, Dirichlet, Neumann, and mixed boundary value problems for the Brinkman system with L ∞ coefficients in Lipschitz domains on a compact Riemannian manifold of dimension m ≥ 2 . The Dirichlet, transmission, and mixed problems for the nonlinear Darcy-Forchheimer-Brinkman system with L ∞ coefficients are also analyzed. First, we focus on the well-posedness of linear transmission, Dirichlet and mixed boundary value problems for the Brinkman system with L ∞ coefficients in Lipschitz domains on compact Riemannian manifolds by using a variational approach that reduces such a boundary value problem to a mixed variational formulation defined in terms of two bilinear continuous forms, one of them satisfying a coercivity condition and another one the inf-sup condition. Further, we show the equivalence between each boundary value problem for the Brinkman system with L ∞ coefficients and its mixed variational counterpart, and then the well posedness in L 2 -based Sobolev spaces by using the Necas-Babuska-Brezzi technique. The second goal of this paper is the construction of the Newtonian and layer potential operators for the Brinkman system with L ∞ coefficients in Lipschitz domains on compact Riemannian manifolds by using the well-posedness results for the analyzed linear transmission problems. Various mapping properties of these operators are also obtained and used to describe the weak solutions of the Poisson problems with Dirichlet and Neumann conditions for the nonsmooth Brinkman system in terms of such potentials. Finally, we combine the well-posedness results of the Poisson problems of Dirichlet, transmission, and mixed type for the nonsmooth Brinkman system with a fixed point theorem in order to show the existence of a weak solution of the Poisson problem of Dirichlet, transmission, or mixed type for the (nonlinear) Darcy-Forchheimer-Brinkman system with L ∞ coefficients in L 2 -based Sobolev spaces in Lipschitz domains on compact Riemannian manifolds of dimension m ∈ { 2 , 3 } .

Though the boundary layer formation in the chemotactic process has been observed in experiment (cf. [63] ), the mathematical study on the boundary layer solutions of chemotaxis models is just in its infant stage. Apart from the sophisticated theoretical tools involved in the analysis, how to impose/derive physical boundary conditions is a state-of-the-art in studying the boundary layer problem of chemotaxis models. This paper will proceed with a previous work [24] in one dimension to establish the convergence of boundary layer solutions of the Keller–Segel model with singular sensitivity in a two-dimensional space (half-plane) with respect to the chemical diffusion rate denoted by e ≥ 0 . Compared to the one-dimensional boundary layer problem, there are many new issues arising from multi-dimensions such as possible Prandtl type degeneracy, curl-free preservation and well-posedness of large-data solutions. In this paper, we shall derive appropriate physical boundary conditions and gradually overcome these barriers and hence establish the convergence of boundary layer solutions of the singular Keller–Segel system in the half-plane as the chemical diffusion rate vanishes. Specially speaking, we justify that the boundary layer converges to the outer layer (solution with e = 0 ) plus the inner layer as e → 0 , where both outer and inner layer profiles are precisely derived and well understood. By doing this, the structure of boundary layer solutions is clearly characterized. We hope that our results and methods can shed lights on the understanding of underlying mechanisms of the boundary layer patterns observed in the experiment for chemotaxis such as the work by Tuval et al. [63] , and open a new window in the future theoretical study of chemotaxis models.

Diaconis and Griffiths (2014) study the multivariate Krawtchouk polynomials orthogonal on the multinomial distribution. In this paper we derive the reproducing kernel orthogonal polynomials Q n ( x , y ; N , p ) on the multinomial distribution which are sums of products of orthonormal polynomials in x and y of fixed total degree n = 0 , 1 , … , N . The Poisson kernel ∑ n = 0 N ρ n Q n ( x , y ; N , p ) arises naturally from a probabilistic argument. An application to a multinomial goodness of fit test is developed, where the chi-squared test statistic is decomposed into orthogonal components which test the order of fit. A new duplication formula for the reproducing kernel polynomials in terms of the 1-dimensional Krawtchouk polynomials is derived. The duplication formula allows a Lancaster characterization of all reversible Markov chains with a multinomial stationary distribution whose eigenvectors are multivariate Krawtchouk polynomials and where eigenvalues are repeated within the same total degree. The χ 2 cutoff time, and total variation cutoff time is investigated in such chains. Emphasis throughout the paper is on a probabilistic understanding of the polynomials and their applications, particularly to Markov chains.

This paper is concerned with the spatiotemporal behaviors of a Gierer–Meinhardt system in discrete time and space form. Through the linear stability analysis, the parametric conditions are gained to ensure the stability of the homogeneous steady state of the system. Based on the bifurcation theory, as well as center manifold theorem, we derive the critical parameter values of the flip, Neimark–Sacker and Turing bifurcation respectively. Besides, the specific parameter expression to form patterns are also determined. In order to identify chaos among regular behaviors, we calculate the Maximum Lyapunov exponents. The results obtained in this paper are illustrated by numerical simulations. From the simulations, we can see some complex dynamics, such as period doubling cascade, invariant cycles, periodic windows, chaotic behaviors, and some striking Turing patterns, e.g. circle, mosaic, spiral, spatiotemporal chaotic patterns and so on, which can be produced by flip-Turing instability, Neimark–Sacker–Turing instability and chaos.

From spline theory it is well-known that univariate cubic spline interpolation, if carried out in its natural Hilbert space W 2 2 [ a , b ] and on point sets with fill distance h , converges only like O ( h 2 ) in L 2 [ a , b ] if no additional assumptions are made. But superconvergence up to order h 4 occurs if more smoothness is assumed and if certain additional boundary conditions are satisfied. This phenomenon was generalized in 1999 to multivariate interpolation in Reproducing Kernel Hilbert Spaces on domains Ω ⊂ R d for continuous positive definite Fourier-transformable shift-invariant kernels on R d . But the sufficient condition for superconvergence given in 1999 still needs further analysis, because the interplay between smoothness and boundary conditions is not clear at all. Furthermore, if only additional smoothness is assumed, superconvergence is numerically observed in the interior of the domain, but a theoretical foundation still is a challenging open problem. This paper first generalizes the “improved error bounds” of 1999 by an abstract theory that includes the Aubin–Nitsche trick and the known superconvergence results for univariate polynomial splines. Then the paper analyzes what is behind the sufficient conditions for superconvergence. They split into conditions on smoothness and localization, and these are investigated independently. If sufficient smoothness is present, but no additional localization conditions are assumed, it is numerically observed that superconvergence always occurs in the interior of the domain, and some supporting arguments are provided. If smoothness and localization interact in the kernel-based case on R d , weak and strong boundary conditions in terms of pseudodifferential operators occur. A special section on Mercer expansions is added, because Mercer eigenfunctions always satisfy the sufficient conditions for superconvergence. Numerical examples illustrate the theoretical findings.

In this paper, we present and analyze a prey-predator system, in which prey species can be infected with some disease. The model presented in this paper is motivated from D. Mukherjee’s model in which he has considered an SI model for the prey species. There are substantial evidences that infected individuals have the ability to recover from the disease if vaccinated/ treated properly. In this regard, Mukherjee’s model is modified by considering SIS model for prey species. Theoretical and numerical simulations show that the recovery of infected prey species plays a crucial role in eliminating the limit cycle oscillations and thus making the interior equilibrium point stable. The possibility of Hopf bifurcation around non zero equilibrium point using the recovery rate as a bifurcation parameter, is discussed. Further, the model is extended by incorporating the harvesting of predator population. A monitory agency has been introduced which monitors the exploitation of resources by implementing certain taxes for each unit biomass of the predator population. The main purpose of the present research is to explore the effect of recovery rate of prey on the dynamics of the system and to optimize the total economical net profits from harvesting of predator species, taking taxation as control parameter.

It has been shown that natural interval extensions (NIE) can be used to calculate the largest positive Lyapunov exponent (LLE). However, the elaboration of NIE are not always possible for some dynamical systems, such as those modelled by simple equations or by Simulink-type blocks. In this paper, we use rounding mode of floating-point numbers to compute the LLE. We have exhibited how to produce two pseudo-orbits by means of different rounding modes; these pseudo-orbits are used to calculate the Lower Bound Error (LBE). The LLE is the slope of the line gotten from the logarithm of the LBE, which is estimated by means of a recursive least square algorithm (RLS). The main contribution of this paper is to develop a procedure to compute the LLE based on the LBE without using the NIE. Additionally, with the aid of RLS the number of required points has been decreased. Eight numerical examples are given to show the effectiveness of the proposed technique.

In a recent paper [R. Qesmi, Dynamics of an opinion model with threshold-type delay, Chaos, Solitons & Fractals 98 (2020). (https://doi.org/10.1016/j.chaos.2020.110379)], we proposed a mathematical model of threshold-type delay differential equations describing the relationship between two subpopulations with opposite opinions and the opinion spread dynamics. The study there showed the possibility of a transcritical forward and backward bifurcations of positive equilibria. In the present paper, we show that the opinion model undergoes a Hopf bifurcation through which one of the bifurcation branches loses the stability and periodic solutions appear. One of the important consequences of the obtained dynamics is that the consensus of the both opinions could be lost by maintaining the balance between the time taken for an individual to become convinced of the outsider opinion, which need be short, and the number of individuals converted to the local opinion which need be low. Finally, we provide numerical simulations to illustrate and support our theoretical results. results.

Focusing on the unpredictability of person-to-person contacts and the complexity of random variations in nature, this paper will formulate a stochastic SIR epidemic model with nonlinear incidence rate and general stochastic noises. First, we derive a stochastic critical value R 0 S related to the basic reproduction number R 0 . Via our new method in constructing suitable Lyapunov function types, we obtain the existence and uniqueness of an ergodic stationary distribution of the stochastic system if R 0 S > 1 . Next, via solving the corresponding Fokker-Planck equation, it is theoretically proved that the stochastic model has a log-normal probability density function when another critical value R 0 H > 1 . Then the exact expression of the density function is obtained. Moreover, we establish the sufficient condition R 0 C 1 for disease extinction. Finally, several numerical simulations are provided to verify our analytical results. By comparison with other existing results, our developed theories and methods will be highlighted to end this paper.

Semi-Markov processes established a rich framework for many real-world problems. Semi-Markov processes include Markov processes, Markov chains, renewal processes, Markov renewal processes, Poisson processes, birth and death processes, and etc. Since semi-Markov processes describe more complex mathematical models of various objects, in many cases the obtained mathematical models are investigated numerically. But in some cases it is possible to study such mathematical models using analytical methods. In the presented paper, the authors chose to do the second way. In this paper we study the semi-Markov random walk processes with negative drift and positive jumps. The random variable or the moment, at which the process for the time reaches in zero level is introduced. An integral equation for the Laplace transform of conditional distribution of this random variable is obtained. In this paper length of jump is given by the gamma distribution with parameters α and β resulting in a fractional order integral equation. In the class of gamma distributions, the resulting general integral equation of convolution type is reduced to a fractional order differential equation with constant coefficients. And also, the exact solution of the resulting fractional differential equation with constant coefficients has been found. Finally, using form of Laplace transform the expectation and variance of the random variable are found.

The present work is devoted to giving new insights into the extended Malkus-Robbins (EMR) dynamo. Firstly, based on differential geometry method, i.e. Kosambi-Cartan-Chern (KCC) theory, the paper investigates the Jacobi stability of the equilibrium and periodic orbit by the eigenvalues of the deviation curvature tensor. The deviation vector is applied to analyze the trajectory behaviors near the equilibrium and periodic orbit. Secondly, the zero-zero-Hopf bifurcation is investigated. The paper obtains the conditions that two periodic solutions appear at the bifurcation point and discusses their stability. Finally, on the global dynamics, the ultimate bound sets of the system are estimated. Numerical simulations are given to verify and visualize the corresponding theoretical results.

Randomness is a common phenomenon in nonlinear systems. And conditions to reach synchronization are more complex and difficult when chaotic systems have random parameters. So in this paper, an adaptive scheme for synchronization of chaotic system with random parameters by using the fuzzy impulsive method and combining the properties of Wiener process and Ito differential is investigated. The main concepts of this paper are applying fuzzy method to approximate the nonlinear part of system, then using Ito differential to study the Wiener process of random parameters of chaotic system, and realizing synchronization under fuzzy impulsive method. The stability is analyzed by Lyapunov stability theorem. At the end of the paper, numerical simulation is presented to illustrate the effectiveness of the results obtained in this paper.

A periodic predator-prey model has been introduced in [5] to study the effect of water level on persistence or extinction of fish populations living in an artificial lake. By using the continuation theorem of Mawhin’s coincidence degree theory, the authors give sufficient conditions for the existence of at least one positive periodic solution. In this paper we study the problem in the general case. We begin by analyzing the invariance, permanence, non-persistence and the globally asymptotic stability for the system. Most interestingly, under additional conditions, we find that the periodic solution obtained in [5] is unique. Finally, in order to make the model system more realistic, we consider the special case when the periodicity in [5] is replaced by almost periodicity. We obtain conditions for existence, uniqueness and stability of a positive almost periodic solution. The methods used in this paper will be comparison theorems and Lyapunov functions. An example is employed to illustrate our result.

This paper presents a new chaotic system that exhibits a two wing (2W) chaotic attractor in its integer order dynamics, three-wing (3W) and four-wing (4W) chaotic attractors in its fractional-order (FO) dynamics, and an eight-wing (8W) attractor in its nonautonomous fractional dynamics. An interesting feature of the proposed system is that two distinct periodic orbits coexist with a strange attractor that gradually evolves into a 4W attractor. The asymmetry, dissimilarity and topological structure of this proposed system with respect to those available in literature, manifest increased irregularity, which in turn indicate more chaos. Besides, the authors have drawn its comparison with various well-known fractional-order chaotic systems (FOCS)s to prove its enhanced features in terms of higher Lyapunov Exponent, fractional order orbital velocities, bandwidth, density, range of dynamical behaviour, etc. A control scheme is proposed to enable switching hybrid synchronisation between the 8W nonautonomous FOCS and the 4W autonomous FOCS, using the former as master and the latter as slave. This work throws light on the potential practical applicability of the proposed system by designing a circuit using minimum circuit components possible, thus signifying the objectives of the paper are finally achieved.

Motivated by the paper Gerard-Varet and Dormy (2010) [6] [JAMS, 2010] about the linear ill-posedness for the Prandtl equations around a shear flow with exponential decay in normal variable, and the recent study of well-posedness on the Prandtl equations in Sobolev spaces, this paper aims to extend the result in [6] to the case when the shear flow has general decay. The key observation is to construct an approximate solution that captures the initial layer to the linearized problem motivated by the precise formulation of solutions to the inviscid Prandtl equations.

Yilmaz et al. (Fluct. Noise Lett. 17, 1830002, 2018) investigated the stochastic phenomenological bifurcations of a generalized Chua circuit driven by Skew-Gaussian distributed noise. They proved it is possible to decrease the number of scrolls by properly choosing the stochastic excitation, manipulating the skewness and noise intensity parameters. Based on the latter, this paper proposes an extension of skew-gaussian noise based on the family of Scale Mixtures of Skew-normal (SMSN) distributions, which includes the skew- t , the skew-gaussian, and the gaussian noises as particular cases. The Lorenz, Generalized Lorenz, Proto–Lorenzand Rossler attractors driven by skew- t distributed noise are considered. Results show that the chaotic regime’s behavior is influenced by the freedom parameter degrees of skew- t noise, increasing the noise variance. This paper concludes that noise intensity increases by rescaling the skew- t distribution at zero mean, rather than by increasing the asymmetry parameter.

In a published paper of Journal of Chaos, Solitons and Fractals, some miss prints were found. One is about the calculation of the solution and also about the expression of the solution using Green’s function of the fractional differential equation studied in this paper; The others are for the properties for Green’s function, these miss prints affected the deriving of the main results in Khan’s paper, some corrections for Khan’s paper should be needed. In this paper, we make some corrections and give the correct proof proceeding for the results, a new example is given to validate part of the proven results.

Reduction theory has played a major role in the study of Hamiltonian systems. Whilst the Hamilton–Jacobi theory is one of the main tools to integrate the dynamics of certain Hamiltonian problems and a topic of research on its own. Moreover, the construction of several symplectic integrators relies on approximations of a complete solution of the Hamilton–Jacobi equation. The natural question that we address in this paper is how these two topics (reduction and Hamilton–Jacobi theory) fit together. We obtain a reduction and reconstruction procedure for the Hamilton–Jacobi equation with symmetries, even in a generalized sense to be clarified below. Several applications and relations to other reduction of the Hamilton–Jacobi theory are shown in the last section of the paper. It is remarkable that as by-product we obtain a generalization of the Ge–Marsden reduction procedure [18] and the results in [17] . Quite surprisingly, the classical ansatze available in the literature to solve the Hamilton–Jacobi equation (see [2] , [19] ) are also particular instances of our framework.

We study the approximation of compact linear operators defined over certain weighted tensor product Hilbert spaces. The information complexity is defined as the minimal number of arbitrary linear functionals needed to obtain an e -approximation for the d -variate problem which is fully determined in terms of the weights and univariate singular values. Exponential tractability means that the information complexity is bounded by a certain function that depends polynomially on d and logarithmically on e − 1 . The corresponding unweighted problem was studied in Hickernell et al. (2020) with many negative results for exponential tractability. The product weights studied in the present paper change the situation. Depending on the form of polynomial dependence on d and logarithmic dependence on e − 1 , we study exponential strong polynomial, exponential polynomial, exponential quasi-polynomial, and exponential ( s , t ) -weak tractability with max ( s , t ) ≥ 1 . For all these notions of exponential tractability, we establish necessary and sufficient conditions on weights and univariate singular values for which it is indeed possible to achieve the corresponding notion of exponential tractability. The case of exponential ( s , t ) -weak tractability with max ( s , t ) 1 is left for future study. The paper uses some general results obtained in Hickernell et al. (2020) and Kritzer and Woźniakowski (2019).

This review paper aims at answering a basic question on the sign of Lyapunov exponents. A few recent papers reported hyperchaotic system having the sign of Lyapunov exponents as (+, +, −, −). Such (+, +, −, −) sign of Lyapunov exponents is in contradiction with the well known (+, +, 0, −) sign of Lyapunov exponents for a 4-D hyperchaotic system. This paper thus discusses various issues related to Lyapunov exponents and proves that the reported sign of (+, +, −, −) Lyapunov exponents is actually (+, +, 0, −) or (+, 0, −, −). This clarification is very important and essential since Lyapunov exponents are the only quantitative measure for the existence of hyperchaos. Three different algorithms are used for calculating the Laypunov exponents to prove the actual sign of Lyapunov exponents for a hyperchaotic system.

The paper explores an eco-epidemiological model with weak Allee in predator, and the disease in the prey population. We consider a predator-prey model with type II functional response. The curiosity of this paper is to consider different competition coefficients within the prey population, which leads to the emergent carrying capacity. We perform the local and global stability analysis of the equilibrium points and the Hopf bifurcation analysis around the endemic equilibrium point. Further we pay attention to the chaotic dynamics which is produced by disease. Our numerical simulations reveal that the three species eco-epidemiological system without weak-Allee induced chaos from stable focus for increasing the force of infection, whereas in the presence of the weak-Allee effect, it exhibits stable solution. We conclude that chaotic dynamics can be controlled by the Allee parameter as well as the competition coefficients. We apply basic tools of non-linear dynamics such as Poincare section and maximum Lyapunov exponent to identify chaotic behavior of the system.

The paper addresses the mathematical analysis of the fishery model in the context of the fractional derivative operator. We use the Caputo–Fabrizio derivative in the investigations. We first prove the fishery model is biologically well definite by proposing the existence and the uniqueness of its solution. The main objective of this paper is to study the dynamics of the predator and the prey in the fishery model when the fractional-order derivative is used. Notably, we analyze the impact of the fractional-order derivative on the dynamics of the fishery model explicitly. To answer this issue, we introduce a new numerical scheme based on the discretization of the fractional integral associated with the Caputo–Fabrizio derivative. The numerical simulations of the solutions of the fractional model are intended to illustrate the numerical scheme presented in our paper. We finish by analyzing the local and global asymptotic stability of the equilibrium points using the Jacobian matrix and the Lyapunov direct method. The Lyapunov function is constructed using standard construction. We notice here that the solutions of the fractional fishery model with different values of the orders describes a cycle when we depict them in three dimensional spaces in time, and furthermore the marine reserves ensure the sustainability of fractional system.

In this paper, we consider a time fractional-order derivative for a diffusive mussel–algae model. The existence of pattern formation was the subject of interest of many previous research works in the case of the diffusive mussel–algae model. Examples include the Turing instability, Hopf bifurcation, Turing-Hopf bifurcation, and others. The presence of the time–fractional–order derivative never been investigated in this model. Next to it ecological relevant, it can generate some important patterns. One of these patterns is produced by the presence of the Turing-Hopf bifurcation. Therefore, our main interest is to analyze the effect of the time fractional–order derivative on the spatiotemporal behavior of the solution, which never been achieved for the mussel-algae model. Besides, Turing–Hopf was studied exclusively on the classical reaction-diffusion systems, where it was also considered for the diffusive mussel-algae model. Thus, our paper puts the fist steps on proving the existence of this type of codimension bifurcation on the diffusive systems with time fractional–order–derivative systems. Further, a suitable numerical simulations are used for confirming the theoretical obtained results.

This paper considers the problem of approximating the spectral factor of continuous spectral densities with finite Dirichlet energy based on finitely many samples of these spectral densities. Although there exists a closed form expression for the spectral factor, this formula shows a very complicated behavior because of the non-linear dependency of the spectral factor from spectral density and because of a singular integral in this expression. Therefore approximation methods are usually applied to calculate the spectral factor. It is shown that there exists no sampling-based method which depends continuously on the samples and which is able to approximate the spectral factor for all densities in this set. Instead, to any sampling-based approximation method there exists a large set of spectral densities so that the approximation method does not converge to the spectral factor for every spectral density in this set as the number of available sampling points is increased. The paper will also show that the same results hold for sampling-based algorithms for the calculation of the outer function in the theory of Hardy spaces.

In this paper, we investigate the dynamics and the global exponential stability of a new class of Hopfield neural network with time-varying and distributed delays. In fact, the properties of norms and the contraction principle are adjusted to ensure the existence as well as the uniqueness of the pseudo almost periodic solution, which is also its derivative pseudo almost periodic. This results are without resorting to the theory of exponential dichotomy. Furthermore, by employing the suitable Lyapunov function, some delay-independent sufficient conditions are derived for exponential convergence. The main originality lies in the fact that spaces considered in this paper generalize the notion of periodicity and almost periodicity. Lastly, two examples are given to demonstrate the validity of the proposed theoretical results.

This paper studies the approximation of the Hilbert transform f ? = H f of continuous functions f with continuous conjugate f ? based on a finite number of samples. It is known that every sequence { H N f } N ? N which approximates f ? from samples of f diverges (weakly) with respect to the uniform norm. This paper conjectures that all of these approximation sequences even contain no convergent subsequence. A property which is termed strong divergence.The conjecture is supported by two results. First it is proven that the sequence of the sampled conjugate Fejer means diverges strongly. Second, it is shown that for every sample based approximation method { H N } N ? N there are functions f such that ? H N f ? ∞ exceeds any given bound for any given number of consecutive indices N .As an application, the later result is used to investigate a problem associated with a question of Ul'yanov on Fourier series which is related to the possibility to construct adaptive approximation methods to determine the Hilbert transform from sampled data. This paper shows that no such approximation method with a finite search horizon exists.

In this paper we report for the first time on the binary generalized synchronization, when for the certain values of the coupling strength two unidirectionally coupled dynamical systems generating the aperiodic binary sequences are in the generalized synchronization regime. The presence of the binary generalized synchronization has been revealed with the help of both the auxiliary system approach and the largest conditional Lyapunov exponent calculation. The mechanism resulting in the binary generalized synchronization has been explained. The finding discussed in this paper gives a strong potential for new applications under many relevant circumstances.

This paper studies the problem of finding optimal parameters for a Poincare section used for determining the type of behavior of a time series: a deterministic or stochastic one. To reach that goal optimization algorithms are coupled with the Poincare & Higuchi (P&H) method, which calculates the Higuchi dimension using points obtained by performing a Poincare section of a certain attractor. The P&H method generates distinctive patterns that can be used for determining if a given attractor is produced by a deterministic or a stochastic system, but this method is sensitive to the parameters of the Poincare section. Patterns generated by the P&H method can be characterized using numerical measures which in turn can be used for finding such parameters for the Poincare section for which the patterns produced by the P&H method are the most prominent. This paper studies several approaches to parameterization of the Poincare section. Proposed approaches are tested on twelve time series, six produced by deterministic chaotic systems and six generated randomly. The obtained results show, that finding good parameters of the Poincare section is important for determining the type of behavior of a time series. Among the tested methods the evolutionary algorithm was able to find the best Poincare sections for use with the P&H method.

In a recent series of papers, Kavitha et al. [2,3,4] solved three inhomogeneous nonlinear Schrodinger (INLS) integro-differential equation under the influence of a variety of nonlinear inhomogeneities and nonlocal damping by the modified extended tangent hyperbolic function method. They obtained several kinds of exact solitary solutions accompanied by the shape changing property. In this paper, we demonstrate that most of exact solutions derived by them do not satisfy the nonlinear equations and consequently are wrong. Furthermore, we study a generalized Hirota equation with spatially-inhomogenetiy and nonlocal nonlinearity. Its integrability is explored through Painleve analysis and N-soliton solutions are obtained based on the Hirota bilinear method. Effects of linear inhomogeneity on the profiles and dynamics of solitons are also investigated graphically.

In this paper, we obtain an explicit formula for the discrepancy between the limit of the second inner variations of p -Laplace Allen–Cahn energies and the second inner variation of their Γ -limit which is the area functional. Our analysis explains the mysterious discrepancy term found in our previous paper [8] in the case p = 2 . The discrepancy term turns out to be related to the convergence of certain 4-tensors which are absent in the usual Allen–Cahn functional. These (hidden) 4-tensors suggest that, in the complex-valued Ginzburg–Landau setting, we should expect a different discrepancy term which we are able to identify. Along the way, we partially answer a question of Kohn and Sternberg [6] by giving a relation between the limit of second variations of the Allen–Cahn functional and the second inner variation of the area functional at local minimizers. Moreover, our analysis reveals an interesting identity connecting second inner variation and Poincare inequality for area-minimizing surfaces with volume constraint in the work of Sternberg and Zumbrun [16] .

In this paper we investigate quasi-Monte Carlo (QMC) integration using digital nets over Z b in reproducing kernel Hilbert spaces. The tent transformation (previously called baker’s transform) was originally used for lattice rules by Hickernell (2002) to achieve higher order convergence of the integration error for smooth non-periodic integrands, and later, has been successfully applied to digital nets over Z 2 by Cristea et al. (2007) and Goda (2015). The aim of this paper is to generalize the latter two results to digital nets over Z b for an arbitrary prime b . For this purpose, we introduce the b -adic tent transformation for an arbitrary positive integer b greater than 1, which is a generalization of the original (dyadic) tent transformation. Further, again for an arbitrary positive integer b greater than 1, we analyze the mean square worst-case error of QMC rules using digital nets over Z b which are randomly digitally shifted and then folded using the b -adic tent transformation in reproducing kernel Hilbert spaces. Using this result, for a prime b , we prove the existence of good higher order polynomial lattice rules over Z b among a smaller number of candidates as compared to the result by Dick and Pillichshammer (2007), which achieve almost the optimal convergence rate of the mean square worst-case error in unanchored Sobolev spaces of smoothness of arbitrary high order.

Growing biological evidence from various ecosystems (Berec et al., 2007) suggests that the Allee effect generated by two or more mechanisms can act simultaneously on a single population. Surprisingly, prey–predator system incorporating multiple Allee effect is relatively poorly studied in literature. In this paper, we consider a ratio dependent prey–predator system with a double Allee effect in prey population growth. We have taken into account the case of the strong and the weak Allee effect separately and study the stability and complete bifurcation analysis. It is shown that the model exhibits the bi-stability and there exists separatrix curve(s) in the phase plane implying that dynamics of the system is very sensitive to the variation of the initial conditions. Generic normal forms at double zero singularity are derived by a rigorous mathematical analysis to show the different types of bifurcation of the proposed model system including Bogdanov–Takens bifurcation, saddle–node bifurcation curve, a Hopf bifurcation curve and a homoclinic bifurcation curve. The complete analysis of possible topological structures including elliptic, parabolic, or hyperbolic orbits, and any combination of them in a neighborhood of the complicated singular point ( 0 , 0 ) in the interior of the first quadrant is explored by using a blow up transformation. These structures have important implications for the global behavior of the model. Numerical simulation results are given to support our theoretical results. Finally, the paper concludes with a discussion of the ecological implications of our analytic and numerical findings.

In this paper, we introduce and study an implicit iterative method to approximate a common solution of split equilibrium problem and fixed point problem for a nonexpansive semigroup in real Hilbert spaces. Further, we prove that the nets generated by the implicit iterative method converge strongly to the common solution of split equilibrium problem and fixed point problem for a nonexpansive semigroup. This common solution is the unique solution of a variational inequality problem and is the optimality condition for a minimization problem. Furthermore, we justify our main result through a numerical example. The results presented in this paper extend and generalize the corresponding results given by Plubtieng and Punpaeng [S. Plubtieng, R. Punpaeng, Fixed point solutions of variational inequalities for nonexpansive semigroups in Hilbert spaces, Math. Comput. Model. 48 (2008) 279–286] and Cianciaruso et al. [F. Cianciaruso, G. Marino, L. Muglia, Iterative methods for equilibrium and fixed point problems for nonexpansive semigroups in Hilbert space, J. Optim. Theory Appl. 146 (2010) 491–509].

Third part of the paper is devoted to analysis of the hyper, hyper-hyper and spatial–temporal chaos of continuous mechanical systems using the Lyapunov exponents. The constructed algorithms for the Lyapunov exponents’ computation allowed detecting and analysing novel phase transitions from chaos through hyper chaos to hyper-hyper chaos. In addition, a novel characteristic “maximal deflection versus excitation amplitude” has been introduced to study stability properties of the investigated continuous systems. It should be emphasized that the latter characteristic yields results in full agreements with those obtained via the Lyapunov exponents’ spectrum estimation. The introduced methods and tools of analysis allowed detecting the Sharkovskii windows of periodicity in all continuous mechanical systems investigated in this paper. Finally, the approach to study the space-temporal chaos exhibited by shell structural-members is also proposed.

In this work chaotic dynamics of continuous mechanical systems such as flexible plates and shallow shells is studied. Namely, a wide class of the mentioned objects is analyzed including flexible plates and cylinder-like panels of infinite length, rectangular spherical and cylindrical shells, closed cylindrical shells, axially symmetric plates, as well as spherical and conical shells. The considered problems are solved by the Bubnov–Galerkin and higher approximation Ritz methods. Convergence and validation of those methods are studied. The Cauchy problems are solved mainly by the fourth Runge-Kutta method, although all variants of the Runge-Kutta methods are considered. New scenarios of transition from regular to chaotic orbits are detected, analyzed and discussed. First part of the paper is devoted to the validation of results obtained. This is why the same infinite length problem is reduced to that of a finite dimension through the FDM (Finite Difference Method) with the approximation order of O ( c 2 ), BGM (Bubnov–Galerkin Method) or RM (Ritz Method) with higher approximations. We pay attention not only to convergence of the mentioned methods regarding the number of partitions of the interval [0, 1] in the FDM or regarding the number of terms in the series applied either in the BGM or RM methods, but we also compare the results obtained via the mentioned different approaches. Furthermore, a so called practical convergence of different Runge-Kutta type methods are tested starting from the second and ending with the eighth order. Second part of the work is devoted to a study of routes to chaos in the so far mentioned mechanical objects. For this purpose the so-called “dynamical charts” are constructed versus control parameters { q 0 , ω p }, where q 0 denotes the loading amplitude, and ω p is the loading frequency. The charts are constructed through analyses of frequency power spectra and the largest Lyapunov exponent (LE). Analysis of the mentioned charts indicates clearly that different routes to chaos exist and allow us to control the objects being investigated. In some cases we detect the classical Feigenbaum scenario and we compute also the Feigenbaum constant. This scenario accompanied all problems which we studied. In addition, we detect and illustrate novel scenarios of transition from regularity into chaos including the Ruelle–Takens–Newhouse–Feigenbaum scenario, and the so called modified Pomeau–Manneville scenario. Third part of the paper is devoted to analysis of the Lyapunov exponents. Namely, while investigating evolutions of vibration regimes of a shell associated with an increase of excitation amplitude q 0 phase transitions chaos–hyper chaos as well as chaos-hyper chaos–hyper–hyper chaos dynamics are illustrated and studied. Furthermore, for all investigated plates and shells the Sharkovskiy windows of periodicity are detected. In particular, a space-temporal chaos/turbulence is studied.

The purpose of this paper is to study the dynamics of a monopolistic firm in a continuous-time framework. The firm is assumed to be boundedly rational and to experience time delays in obtaining and implementing information on output. The dynamic adjustment process is based on the gradient of the expected profit. The paper is divided into three parts: we examine delay effects on dynamics caused by one-time delay and two-time delays in the first two parts. Global dynamics and analytical results on local dynamics are numerically confirmed in the third part. Four main results are demonstrated. First, the stability switch from stability to instability occurs only once in the case of a single delay. Second, the alternation of stability and instability can continue if two time delays are involved. Third, the occurence of Hopf bifurcation is analytically shown if stability is lost. Finally, in a bifurcation process, there are a period-doubling cascade to chaos and a period-halving cascade to the equilibrium point in the case of two time delays if the difference between the two delays is large.

We say that the polynomial sequence $(Q^{(\lambda)}_n)$ is a semiclassical Sobolev polynomial sequence when it is orthogonal with respect to the inner product $$ _S= +\lambda , $$ where ${\bf u}$ is a semiclassical linear functional, ${\mathscr D}$ is the differential, the difference or the $q$--difference operator, and $\lambda$ is a positive constant. In this paper we get algebraic and differential/difference properties for such polynomials as well as algebraic relations between them and the polynomial sequence orthogonal with respect to the semiclassical functional $\bf u$. The main goal of this article is to give a general approach to the study of the polynomials orthogonal with respect to the above nonstandard inner product regardless of the type of operator ${\mathscr D}$ considered. Finally, we illustrate our results by applying them to some known families of Sobolev orthogonal polynomials as well as to some new ones introduced in this paper for the first time., Comment: 23 pages, special issue dedicated to Professor Guillermo Lopez lagomasino on the occasion of his 60th birthday, accepted in Journal of Approximation Theory

In this paper our aim is to show some mean value inequalities for the modified Bessel functions of the first and second kinds. Our proofs are based on some bounds for the logarithmic derivatives of these functions, which are in fact equivalent to the corresponding Tur\'an type inequalities for these functions. As an application of the results concerning the modified Bessel function of the second kind we prove that the cumulative distribution function of the gamma-gamma distribution is log-concave. At the end of this paper several open problems are posed, which may be of interest for further research., Comment: 14 pages

In this paper, we rigorously prove the existence and stability of single-peaked patterns for the singularly perturbed Gierer–Meinhardt system on a compact two-dimensional Riemannian manifold without boundary which are far from spatial homogeneity. Throughout the paper we assume that the activator diffusivity ϵ 2 is small enough. We show that for the threshold ratio D ∼ 1 ϵ 2 of the activator diffusivity ϵ 2 and the inhibitor diffusivity D, the Gaussian curvature and the Green's function interact. A convex combination of the Gaussian curvature and the Green's function together with their derivatives are linked to the peak locations and the o ( 1 ) eigenvalues. A nonlocal eigenvalue problem (NLEP) determines the O ( 1 ) eigenvalues which all have negative part in this case.

At first, the paper described a brief history of communication in China. In 2008, the third generation (3G) telecommunication technologies began to develop. In Chinese 3G telecommunication market, there are three oligopolistic competitors. In order to study competing process of the triopoly, the paper considered a Bertrand model with bounded rational. A triopoly game is modeled by three nonlinear difference equations. By using the theory of bifurcations of dynamical systems, the existence and stability for the equilibria of this system were obtained. Numerical simulations were used to show bifurcations diagrams, the maximum Lyapunov exponent and sensitive dependence on initial conditions. We observed that increasing the speed of adjustment of bounded rational player might change the stability of Nash equilibrium point and cause bifurcation and chaos to occur. Thus, it will be difficult to make strategy and to forecast for each enterprise. It might cause the telecom markets to be irregular. The analysis and results in this paper are important for Chinese telecommunication markets, the same to mathematics.

Period-doubling bifurcation and its route to chaos have been thoroughly investigated in voltage-mode and current-mode controlled DC motor drives under simple proportional control. In this paper, the phenomena of Hopf bifurcation and chaos from torus breakdown in a voltage-mode controlled DC drive system is reported. It has been shown that Hopf bifurcation may occur when the DC drive system adopts a more practical proportional-integral control. The phenomena of period-adding and phase-locking are also observed after the Hopf bifurcation. Furthermore, it is shown that the stable torus can breakdown and chaos emerges afterwards. The work presented in this paper provides more complete information about the dynamical behaviors of DC drive systems.

This paper illustrates an approach to generate multi-wing attractors in coupled Lorenz systems. In particular, novel four-wing (eight-wing) hyperchaotic attractors are generated by coupling two (three) identical Lorenz systems. The paper shows that the equilibria of the proposed systems have certain symmetries with respect to specific coordinate planes and the eigenvalues of the associated Jacobian matrices exhibit the property of similarity. In analogy with the original Lorenz system, where the two-wings of the butterfly attractor are located around the two equilibria with the unstable pair of complex-conjugate eigenvalues, this paper shows that the four-wings (eight-wings) of these attractors are located around the four (eight) equilibria with two (three) pairs of unstable complex-conjugate eigenvalues.

In this paper, several classes of impulsive Cohen–Grossberg neural networks with continuously distributed delays are considered. Global exponential stability and robust global exponential stability of the equilibrium solution are investigated by using Lyapunov function and integro-differential inequality. Moreover, sufficient conditions are also given to guarantee the existence of ϖ-periodic solution and that all other solutions are convergent to it globally exponentially. Finally, two examples are given to demonstrate the effectiveness of our results in this paper.

In this paper, the anti-synchronization of a general class of chaotic systems is investigated. A simple adaptive feedback scheme is proposed to anti-synchronize many familiar chaotic systems, including autonomous systems and non-autonomous systems. Lyapunov analysis for the error system gives the asymptotic stability conditions based on the invariance principle of differential equations. The schemes are successfully applied to three groups of examples: the van der Pol–Duffing oscillator, the parametrically harmonically excited 4D system, and the additionally harmonically excited Murali–Lakshmanan–Chua circuit. Numerical results are presented to justify the theoretical analysis in this paper.