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.

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.

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.

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.

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.

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.

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.

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.

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.

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).

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 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.

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.

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.

We consider the Schrodinger equations with arbitrary (large) power non-linearity on the three-dimensional torus. We construct non-trivial probability measures supported on Sobolev spaces and show that the equations are globally well-posed on the supports of these measures, respectively. Moreover, these measures are invariant under the flows that are constructed. Therefore, the constructed solutions are recurrent in time. Also, we show slow growth control on the time evolution of the solutions. A generalization to any dimension is given. Our proof relies on a new approach combining the fluctuation-dissipation method and some features of the Gibbs measures theory for Hamiltonian PDEs. The strategy of the paper applies to other contexts.

This paper studies a dissipative heat equation with Delta sources of non-linear strength located on a periodic lattice. The model arises from intracellular waves in continuum excitable media with discrete release sites. Due to the presence of Delta sources, the solution of the model has discontinuous spatial derivatives. We focus on the bistable regime of the model, determined by the decay strength parameter a and the separation distance L between release sites, in which the model admits exactly three L-periodic steady states. We establish the existence of pulsating waves spatially connecting them. For the case of waves connecting two stable L-periodic steady states, the uniqueness and global exponential stability of pulsating waves are shown. Also a new technique is introduced to find the fine structure of the tails of pulsating waves.

This paper focuses on the Dysthe equation which is a higher order approximation of the water waves system in the modulation (Schrodinger) regime and in the infinite depth case. We first review the derivation of the Dysthe and related equations. Then we study the initial-value problem. We prove a small data global well-posedness and scattering result in the critical space L 2 ( R 2 ) . This result is sharp in view of the fact that the flow map cannot be C 3 continuous below L 2 ( R 2 ) . Our analysis relies on linear and bilinear Strichartz estimates in the context of the Fourier restriction norm method. Moreover, since we are at a critical level, we need to work in the framework of the atomic space U S 2 and its dual V S 2 of square bounded variation functions. We also prove that the initial-value problem is locally well-posed in H s ( R 2 ) , s > 0 . Our results extend to the finite depth version of the Dysthe equation.

We consider the initial boundary value problem for the isentropic compressible Navier-Stokes equations with cylindrical symmetry. The existence of boundary layers is well-known when the shear viscosity vanishes. In this paper, we derive explicit Prandtl type boundary layer equations and prove the global in time stability of the boundary layer profile together with the optimal convergence rate of the vanishing shear viscosity limit without any smallness assumption on the initial and boundary data.

In this paper we study the structure of positive vector solutions for nonlinear Schrodinger systems with 3 components when all interspecies interaction forces are positive and large while all intraspecies interaction forces are positive and fixed. We will show that the structure strongly depends on some relation of large interspecies interaction forces.

We consider the following k-coupled nonlinear Schrodinger systems: { − Δ u j + λ j u j = μ j u j 3 + ∑ i = 1 , i ≠ j k β i , j u i 2 u j in R N , u j > 0 in R N , u j ( x ) → 0 as | x | → + ∞ , j = 1 , 2 , ⋯ , k , where N ≤ 3 , k ≥ 3 , λ j , μ j > 0 are constants and β i , j = β j , i ≠ 0 are parameters. There have been intensive studies for the above systems when k = 2 or the systems are purely attractive ( β i , j > 0 , ∀ i ≠ j ) or purely repulsive ( β i , j 0 , ∀ i ≠ j ); however very few results are available for k ≥ 3 when the systems admit mixed couplings and the components are organized into groups, i.e., there exist ( i 1 , j 1 ) and ( i 2 , j 2 ) such that β i 1 , j 1 > 0 and β i 2 , j 2 0 . In this paper we give the first systematic and an (almost) complete study on the existence of ground states when the systems admit mixed couplings and the components are organized into groups. We first divide these systems into repulsive-mixed and total-mixed cases. In the first case we prove nonexistence of ground states. In the second case we give a necessary condition for the existence of ground states and also provide estimates for Morse index. The key idea is the block decomposition of the systems (optimal block decompositions, eventual block decompositions), and the measure of total interaction forces between different blocks. Finally the assumptions on the existence of ground states are shown to be optimal in some special cases.

In this paper, we study large time behavior of complex-valued solutions to nonlinear Klein-Gordon equation with a gauge invariant quadratic nonlinearity in two spatial dimensions. To find a possible asymptotic behavior, we consider the final value problem. It turns out that one possible behavior is a linear solution with a logarithmic phase correction as in the real-valued case. However, the shape of the logarithmic correction term has one more parameter which is also given by the final data. In the real case the parameter is constant so one cannot see its effect. However, in the complex case it varies in general. The one dimensional case is also discussed., Comment: 25 papges, 2 figures

This paper investigates the dynamic behavior of a diffusive plankton model with defensive and offensive effects in two cases. For the single compartment model, we first derive the sufficient conditions for the stability and Hopf bifurcation of coexisting equilibrium, which implies that the changes of defense and offense can cause oscillation of planktonic population. Then the properties of Hopf bifurcation are discussed by center manifold theorem. For the spatially extended model, we obtain the sufficient conditions for Turing instability and Hopf bifurcation. It is observed that spatial patterns put in place, under the interaction of diffusion, defense and offense. Finally, some numerical simulations are carried out to support the analytical results.

This paper is concerned with an operational matrix method based on the shifted Legendre cardinal functions for solving the nonlinear variable-order time fractional Schrodinger equation. The variable-order fractional derivative operator is defined in the Atangana–Baleanu–Caputo sense. Through the way, a new operational matrix of variable-order fractional derivative is derived for the shifted Legendre cardinal functions and used in the established method. More precisely, the unknown solution is separated into the real and imaginary parts, and then these parts are expanded in terms of the shifted Legendre cardinal functions with undetermined coefficients. These expansions are substituted into the main equation and the generated operational matrix is utilized to extract a system of nonlinear algebraic equations. Thereafter, the yielded system is solved to obtain an approximate solution for the problem. The precision of the established approach is examined through various types of test examples. Numerical simulations confirm that the suggested approach is high accurate in providing satisfactory results.

This paper aims to study the critical sectional area (CSA) of surge chamber considering the nonlinearity of head loss of diversion tunnel and steady output of turbine. Firstly, three basic equations for hydropower station with surge chamber are established. Four mathematical models for the derivation of CSA of surge chamber are constructed. Then, the stability of hydropower station with surge chamber is analyzed by Hopf bifurcation. Based on the critical stable state of hydropower station, the formulas for the CSA of surge chamber are derived. Finally, the verification and comparison of different CSAs are conducted. The correctness and rationality of obtained formulas are explained. The results indicate that, under load increase operation condition, both the nonlinearity of head loss of diversion tunnel and nonlinearity of steady output of turbine can reduce the value of CSA of surge chamber and are favorable for the stability of hydropower station. Under load decrease operation condition, the rules are opposite. Under both load increase and load decrease operation conditions, the effect of the nonlinearity of head loss of diversion tunnel on CSA of surge chamber is much more significant than that of the nonlinearity of steady output of turbine. The formula for CSA of surge chamber considering both the nonlinearity of head loss of diversion tunnel and nonlinearity of steady output of turbine can be expressed as an amplification coefficient times of the Thoma formula. That formula has a higher precision than Thoma formula.

In this paper we are established the existence of positive solutions (EPS) and the Hyers-Ulam (HU) stability of a general class of nonlinear Atangana-Baleanu-Caputo (ABC) fractional differential equations (FDEs) with singularity and nonlinear p-Laplacian operator in Banach’s space. To find the solution for the EPS, we use the Guo-Krasnoselskii theorem. The fractional differential equation is converted into an alternative integral structure using the Atangana-Baleanu fractional integral operator. Also, HU-stability is analyzed. We include an example with specific parameters and assumptions to show the results of the proposal.

Bessel processes $(X_{t,k})_{t\ge0}$ in $N$ dimensions are classified via associated root systems and multiplicity constants $k\ge0$. They describe interacting Calogero-Moser-Suther\-land particle systems with $N$ particles and are related to $\beta$-Hermite and $\beta$-Laguerre ensembles. Recently, several central limit theorems were derived for fixed $t>0$, fixed starting points, and $k\to\infty$. In this paper we extend the CLT in the A-case from start in 0 to arbitrary starting distributions by using a limit result for the corresponding Bessel functions. We also determine the eigenvalues and eigenvectors of the covariance matrices of the Gaussian limits and study applications to CLTs for the intermediate particles for $k\to\infty$ and then $N\to\infty$., Comment: 20 pages

In this paper, the third order partial differential equation defined by Caputo fractional derivative with Atangana–Baleanu derivative has been investigated. The stability estimates are proved for the exact solution. Difference schemes for Crank–Nicholson finite difference scheme method is constructed. The stability of difference schemes for this problem is shown by Von Neumann method (Fourier analysis method). Numerical results with respect to the exact solution confirm the accuracy and effectiveness of the technique. The reproducing kernel function for the problem has been found.