Divided-difference equation and three-term recurrence relations of some systems of bivariate q-orthogonal polynomials

Partial divided-difference equations and three-term recurrence relations satisfied by the bivariate Askey–Wilson and the bivariate q-Racah polynomials are computed in this work. By using limiting processes, partial divided (or q)-difference equations and three-term recurrence relations are also provided for each of the following families of orthogonal polynomials: the bivariate continuous dual q-Hahn, the bivariate Al-Salam-Chihara, the bivariate continuous q-Hahn, the bivariate q-Hahn, the bivariate dual q-Hahn, the bivariate q-Krawtchouk, the bivariate q-Meixner, and the bivariate q-Charlier polynomials.     

Non-synchronous oscillations in four-dimensional nonlinear semelparous Leslie matrix models

In a certain class of semelparous Leslie matrix models, either a positive equilibrium is stable and an invariant set on the boundary of the nonnegative cone is unstable or vice versa generically if the model dimension is two or three. This dynamic dichotomy is expected to be failed in the four-dimensional case. Our study focuses on a semelparous Leslie matrix model with specific nonlinearities and rigorously proves that the dynamic dichotomy does not hold in the four-dimensional case. This result is derived by showing that the four-dimensional semelparous Leslie matrix model can be uniformly persistent with respect to the boundary of the nonnegative cone even if there exists an unstable positive equilibrium. In such a situation, there are no missing age-classes but population oscillation occurs.     

Complex dynamics of a discrete predator–prey model with the prey subject to the Allee effect

In this paper, complex dynamics of the discrete predator–prey model with the prey subject to the Allee effect are investigated in detail. Firstly, when the prey intrinsic growth rate is not large, the basins of attraction of the equilibrium points of the single population model are given. Secondly, rigorous results on the existence and stability of the equilibrium points of the model are derived, especially, by analyzing the higher order terms, we obtain that the non-hyperbolic extinction equilibrium point is locally asymptotically stable. The existences and bifurcation directions for the flip bifurcation, the Neimark–Sacker bifurcation and codimension-two bifurcations with 1:2 resonance are derived by using the center manifold theorem and the bifurcation theory. We derive that the model only exhibits a supercritical flip bifurcation and it is possible for the model to exhibit a supercritical or subcritical Neimark–Sacker bifurcation at the larger positive equilibrium point. Chaos in the sense of Marotto is proved by analytical methods. Finally, numerical simulations including bifurcation diagrams, phase portraits, sensitivity dependence on the initial values, Lyapunov exponents display new and rich dynamical behaviour. The analytic results and numerical simulations demonstrate that the Allee effect plays a very important role for dynamical behaviour.     

Dominant and recessive solutions at infinity and genera of conjoined bases for discrete symplectic systems

In this paper we introduce the theory of dominant solutions at infinity for nonoscillatory discrete symplectic systems without any controllability assumption. Such solutions represent an opposite concept to recessive solutions at infinity, which were recently developed for such systems by the authors. Our main results include: (i) the existence of dominant solutions at infinity for all ranks in a given range depending on the order of abnormality of the system, (ii) construction of dominant solutions at infinity with eventually the same image, (iii) classification of dominant and recessive solutions at infinity with eventually the same image, (iv) limit characterization of recessive solutions at infinity in terms of dominant solutions at infinity and vice versa, and (v) Reid’s construction of the minimal recessive solution at infinity. These results are based on a new theory of genera of conjoined bases for symplectic systems developed for this purpose in this paper.     

Efficient Computation of the Joint Sample Frequency Spectra for Multiple Populations

A wide range of studies in population genetics have employed the sample frequency spectrum (SFS), a summary statistic which describes the distribution of mutant alleles at a polymorphic site in a sample of DNA sequences and provides a highly efficient dimensional reduction of large-scale population genomic variation data. Recently, there has been much interest in analyzing the joint SFS data from multiple populations to infer parameters of complex demographic histories, including variable population sizes, population split times, migration rates, admixture proportions, and so on. SFS-based inference methods require accurate computation of the expected SFS under a given demographic model. Although much methodological progress has been made, existing methods suffer from numerical instability and high computational complexity when multiple populations are involved and the sample size is large. In this article, we present new analytic formulas and algorithms that enable accurate, efficient computation of the expected joint SFS for thousands of individuals sampled from hundreds of populations related by a complex demographic model with arbitrary population size histories (including piecewise-exponential growth). Our results are implemented in a new software package called momi (MOran Models for Inference). Through an empirical study, we demonstrate our improvements to numerical stability and computational complexity.     

N-dimensional Schrödinger equation at finite temperature using the Nikiforov–Uvarov method

The N-radial Schrödinger equation is analytically solved. The Cornell potential is extended to finite temperature. The energy eigenvalues and the wave functions are calculated in the N-dimensional form using the Nikiforov–Uvarov (NV) method. At zero temperature, the energy eigenvalues and the wave functions are obtained in good agreement with other works. The present results are applied on the charmonium and bottomonium masses at finite temperature. The effect of dimensionality number is investigated on the quarkonium masses. A comparison is discussed with other works, which use the QCD sum rules and lattice QCD. The present approach successfully generalizes the energy eigenvalues and corresponding wave functions at finite temperature in the N-dimensional representation. In addition, the present approach can successfully be applied to the quarkonium systems at finite temperature.     

Solitary wave solution of the generalized Hirota–Satsuma coupled KdV system

In this research, we find the exact traveling wave solutions involving parameters of the generalized Hirota–Satsuma couple KdV system according to the modified simple equation method with the aid of Maple 16. When these parameters are taken special values, the solitary wave solutions are derived from the exact traveling wave solutions. It is shown that the modified simple equation method provides an effective and a more powerful mathematical tool for solving nonlinear evolution equations in mathematical physics. Comparison between our results and the well-known results will be presented.