The regularity of the multiple higher-order poles solitons of the NLS equation

Based on the inverse scattering method, the formulae of one higher-order pole solitons and multiple higher-order poles solitons of the nonlinear Schrödinger equation (NLS) equation are obtained. Their denominators are expressed as , where is a matrix frequently constructed for solving the Riemann-Hilbert problem, and the asterisk denotes complex conjugate. We take two methods for proving is invertible. The first one shows matrix is equivalent to a self-adjoint Hankel matrix , proving . The second one considers the block-matrix form of , proving . In addition, we prove that the dimension of is equivalent to the sum of the orders of pole points of the transmission coefficient and its diagonal entries compose a set of basis.     

Finite‐size corrections at the hard edge for the Laguerre β ensemble

A fundamental question in random matrix theory is to quantify the optimal rate of convergence to universal laws. We take up this problem for the Laguerre β ensemble, characterized by the Dyson parameter β, and the Laguerre weight , in the hard edge limit. The latter relates to the eigenvalues in the vicinity of the origin in the scaled variable . Previous work has established the corresponding functional form of various statistical quantities—for example, the distribution of the smallest eigenvalue, provided that . We show, using the theory of multidimensional hypergeometric functions based on Jack polynomials, that with the modified hard edge scaling , the rate of convergence to the limiting distribution is , which is optimal. In the case , general the explicit functional form of the distribution of the smallest eigenvalue at this order can be computed, as it can for and general . An iterative scheme is presented to numerically approximate the functional form for general .     

A difference equation approach to Plancherel-Rotach asymptotics for -orthogonal polynomials

In this paper, we employ a difference equation approach to study the Plancherel-Rotach asymptotics of -orthogonal polynomials about their largest zeros. Our method for -difference equations is an analogue to the turning point problem for Hermite differential equations. It works well in the toy problems of Stieltjes-Wigert polynomials and -Hermite polynomials.     

The Riemann–Hilbert analysis to the Pollaczek–Jacobi type orthogonal polynomials

In this paper, we study polynomials orthogonal with respect to a Pollaczek–Jacobi type weight The uniform asymptotic expansions for the monic orthogonal polynomials on the interval (0,1) and outside this interval are obtained. Moreover, near , the uniform asymptotic expansion involves Airy function as , and Bessel function of order α as in the neighborhood of , the uniform asymptotic expansion is associated with Bessel function of order β as . The recurrence coefficients and leading coefficient of the orthogonal polynomials are expressed in terms of a particular Painlevé III transcendent. We also obtain the limit of the kernel in the bulk of the spectrum. The double scaled logarithmic derivative of the Hankel determinant satisfies a σ‐form Painlevé III equation. The asymptotic analysis is based on the Deift and Zhou's steepest descent method.     

Uniform asymptotics and zeros of the associated Pollaczek polynomials

Using a difference-equation method in a previous paper, we study the associated Pollaczek polynomials defined by a three-term recurrence relation. Two asymptotic approximations are derived for these polynomials; one holds for with and , and the other holds for with t in a neighborhood of . An asymptotic formula is also provided for their largest zeros.     

Stable blow-up dynamics in the -critical and -supercritical generalized Hartree equation

We study stable blow-up dynamics in the generalized Hartree equation with radial symmetry, which is a Schrödinger-type equation with a nonlocal, convolution-type nonlinearity: First, we consider the -critical case in dimensions and obtain that a generic blow-up has a self-similar structure and exhibits not only the square root blowup rate , but also the log-log correction (via asymptotic analysis and functional fitting), thus, behaving similarly to the stable blow-up regime in the -critical nonlinear Schrödinger equation. In this setting, we also study blow-up profiles and show that generic blow-up solutions converge to the rescaled , a ground state solution of the elliptic equation . We also consider the -supercritical case in dimensions . We derive the profile equation for the self-similar blow-up and establish the existence and local uniqueness of its solutions. As in the NLS -supercritical regime, the profile equation exhibits branches of nonoscillating, polynomially decaying (multi-bump) solutions. A numerical scheme of putting constraints into solving the corresponding ordinary differential equation is applied during the process of finding the multi-bump solutions. Direct numerical simulation of solutions to the generalized Hartree equation by the dynamic rescaling method indicates that the is the profile for the stable blow-up. In this supercritical case, we obtain the blow-up rate without any correction. This blow-up happens at the focusing level , and thus, numerically observable (unlike the -critical case). In summary, we find that the results are similar to the behavior of stable self-similar blowup solutions in the corresponding settings for the nonlinear Schrödinger equation. Consequently, one may expect that the form of the nonlinearity in the Schrödinger-type equations is not essential in the stable formation of singularities.     

Development and analysis of a malaria transmission mathematical model with seasonal mosquito life-history traits

In this paper, we develop and analyze a malaria model with seasonality of mosquito life-history traits: periodic-mosquitoes per capita birth rate, -mosquitoes death rate, -probability of mosquito to human disease transmission, -probability of human to mosquito disease transmission, and -mosquitoes biting rate. All these parameters are assumed to be time dependent leading to a nonautonomous differential equation system. We provide a global analysis of the model depending on two threshold parameters and (with ). When , then the disease-free stationary state is locally asymptotically stable. In the presence of the human disease-induced mortality, the global stability of the disease-free stationary state is guarantied when . On the contrary, if , the disease persists in the host population in the long term and the model admits at least one positive periodic solution. Moreover, by a numerical simulation, we show that a sub-critical (backward) bifurcation is possible at . Finally, the simulation results are in accordance with the seasonal variation of the reported cases of a malaria-epidemic region in Mpumalanga province in South Africa.     

Inverse scattering transform for the defocusing Ablowitz–Ladik system with arbitrarily large nonzero background

In this work we develop the inverse scattering transform (IST) for the defocusing Ablowitz–Ladik (AL) equation with arbitrarily a large nonzero background at space infinity. The IST was developed in previous works under the assumption that the amplitude of the background satisfies a “small norm” condition . On the other hand, Ohta and Yang recently showed that the defocusing AL system, which is modulationally stable for , becomes unstable if , and exhibits discrete rogue wave solutions, some of which are regular for all times. Here, we construct the IST for the defocusing AL with , analyze the spectrum, and characterize the soliton and rational solutions from a spectral point of view. We formulate the direct and inverse problems by using a suitable uniformization variable, and pose the inverse problem as an RHP across a simple contour in the complex plane of the uniform variable. As a by‐product of the IST, we also obtain explicit soliton solutions, which are the discrete analog of the celebrated Kuznetsov–Ma, Akhmediev, Peregrine solutions, and which mimic the corresponding solutions for the focusing AL equation. Soliton solutions that are the analog of the dark soliton solutions of the defocusing AL equation in the case are also presented.     

Using Symmetries à Rebours

We consider nonclassical symmetries of partial differential equations (PDEs) in dimensions. Given a th‐order ordinary differential equation in the unknown we are able to find the most general scalar PDE of a given order which can be reduced via a nonclassical symmetry to .     

Commutativity of quaternion‐matrix–valued functions and quaternion matrix dynamic equations on time scales

In this paper, we obtain some basic results of quaternion algorithms and quaternion calculus on time scales. Based on this, a Liouville formula and some related properties are derived for quaternion dynamic equations on time scales through conjugate transposed matrix algorithms. Moreover, we introduce the quaternion matrix exponential function by homogeneous quaternion matrix dynamic equations. Also a corresponding existence and uniqueness theorem is proved. In addition, the commutativity of quaternion‐matrix–valued functions is investigated and some sufficient and necessary conditions of commutativity and noncommutativity are established on time scales. Also the fundamental solution matrices of some basic quaternion matrix dynamic equations are obtained. Examples are provided to illustrate the results, which are completely new on hybrid domains particularly when the time scales are the quantum case and the discrete case ; , both of which are significant for the study of quaternion q‐dynamic equations and quaternion difference dynamic equations. Finally, we present several applications including multidimensional rotations and transformations of the submarine, the gyroscope, and the planet whose dynamical behaviors are depicted by quaternion dynamics on time scales and the corresponding iteration numerical solution for homogeneous quaternion dynamic equations are provided on various time scales.     

Connection formulae for asymptotics of the fifth Painlevé transcendent on the imaginary axis: I

Leading terms of asymptotic expansions for the general complex solutions of the fifth Painlevé equation as are found. These asymptotics are parameterized by monodromy data of the associated linear ODE, The parameterization allows one to derive connection formulae for the asymptotics. We provide numerical verification of the results. Important special cases of the connection formulae are also considered.     

Invariant discrete flows

In this paper, we investigate the evolution of joint invariants under invariant geometric flows using the theory of equivariant moving frames and the induced invariant discrete variational complex. For certain arc length preserving planar curve flows invariant under the special Euclidean group , the special linear group , and the semidirect group , we find that the induced evolution of the discrete curvature satisfies the differential‐difference mKdV, KdV, and Burgers' equations, respectively. These three equations are completely integrable, and we show that a recursion operator can be constructed by precomposing the characteristic operator of the curvature by a certain invariant difference operator. Finally, we derive the constraint for the integrability of the discrete curvature evolution to lift to the evolution of the discrete curve itself.     

Well‐posedness,positivity, and time asymptotics properties for a reaction–diffusion model of plankton communities,involving a rational nonlinearity with singularity

In this work, we consider a reaction–diffusion system, modeling the interaction between nutrients, phytoplankton, and zooplankton. Using a semigroup approach in , we prove global existence, uniqueness, and positivity of the solutions. The nonlinearity is handled by providing estimates in , allowing to deal with most of the functional responses that describe predator/prey interactions (Holling I, II, III, Ivlev) in ecology. The paper finally exhibits some time asymptotic properties of the solutions.     

Solutions with concentration for conservation laws with discontinuous flux and its applications to numerical schemes for hyperbolic systems

Measure-valued weak solutions for conservation laws with discontinuous flux are proposed and explicit formulae have been derived. We propose convergent discontinuous flux-based numerical schemes for the class of hyperbolic systems that admit nonclassical -shocks, by extending the theory of discontinuous flux for nonlinear conservation laws to scalar transport equation with a discontinuous coefficient. The article also discusses the concentration phenomenon of solutions along the line of discontinuity, for scalar transport equations with a discontinuous coefficient. The existence of the solutions for transport equation is shown using the vanishing viscosity approach and the asymptotic behavior of the solutions is also established. The performance of the numerical schemes for both scalar conservation laws and systems to capture the -shocks effectively is displayed through various numerical experiments.     

Epidemic dynamics and spatial segregation driven by cognitive diffusion and nonlinear incidence

We propose susceptible-infected-susceptible epidemic reaction–diffusion models with cognitive movement and nonlinear incidence $S^q{I}^p$ $(p,q>0)$ in a spatially heterogeneous environment. The cognitive dispersal term takes either random diffusion or symmetric diffusion. Building upon the $L^{\infty }$-estimates of positive solutions under $p,q>0$, we state the asymptotic dynamics for , $q>0$. The numerical results reveal spatial segregation of susceptible and infected populations: (a) the heterogeneous random diffusion can segregate the population and reduce the infection fraction significantly; (b) the segregation phenomenon disappears as the ratio $p/q$ approaches one from below; (c) the disease-free region strengthens the segregation induced by heterogeneous random diffusion; (d) the segregation governed by random diffusion is more sensitive to the incidence mechanism; (e) the distribution of steady states driven by symmetric diffusion is always similar to that by homogeneous diffusion.     

Gaussian unitary ensembles with two jump discontinuities,PDEs, and the coupled Painlevé II and IV systems

We consider the Hankel determinant generated by the Gaussian weight with two jump discontinuities. Utilizing the results of Min and Chen [Math. Methods Appl Sci. 2019;42:301‐321] where a second‐order partial differential equation (PDE) was deduced for the log derivative of the Hankel determinant by using the ladder operators adapted to orthogonal polynomials, we derive the coupled Painlevé IV system which was established in Wu and Xu [arXiv: 2002.11240v2] by a study of the Riemann‐Hilbert problem for orthogonal polynomials. Under double scaling, we show that, as , the log derivative of the Hankel determinant in the scaled variables tends to the Hamiltonian of a coupled Painlevé II system and it satisfies a second‐order PDE. In addition, we obtain the asymptotics for the recurrence coefficients of orthogonal polynomials, which are connected with the solutions of the coupled Painlevé II system.     

A Szegő limit theorem for translation‐invariant operators on polygons

We prove Szeg?‐type trace asymptotics for translation‐invariant operators on polygons. More precisely, consider a Fourier multiplier on with a sufficiently decaying, smooth symbol . Let be the interior of a polygon and, for , define its scaled version . Then we study the spectral asymptotics for the operator , the spatial restriction of A onto : for entire functions h with we provide a complete asymptotic expansion of as . These trace asymptotics consist of three terms that reflect the geometry of the polygon. If P is replaced by a domain with smooth boundary, a complete asymptotic expansion of the trace has been known for more than 30 years. However, for polygons the formula for the constant order term in the asymptotics is new. In particular, we show that each corner of the polygon produces an extra contribution; as a consequence, the constant order term exhibits an anomaly similar to the heat trace asymptotics for the Dirichlet Laplacian.     

General soliton and (semi-)rational solutions to the nonlocal Mel'nikov equation on the periodic background

In this paper, the Hirota's bilinear method and Kadomtsev-Petviashvili hierarchy reduction method are applied to construct soliton, line breather and (semi-)rational solutions to the nonlocal Mel'nikov equation with nonzero boundary conditions. These solutions are expressed as Gram-type determinants. When N is even, soliton, line breather and (semi-)rational solutions on the constant background are derived while these solutions are located on the periodic background for odd N. Regularity of these solutions and their connections with the local Mel'nikov equation are analyzed for proper choices of parameters that appear in the solutions. The dynamics of the solutions are discussed in detail. All possible configurations of soliton and lump solutions are found for . Several interesting dynamical behaviors of semi-rational solutions are observed. It is shown that certain lumps may exhibit fusion and fission phenomena during their interactions with solitons while some lump may change its direction of movement after it collides with solitons.