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 .     

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.     

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.     

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.     

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.     

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 .     

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.     

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.     

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.     

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.     

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.     

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.     

A complete solution to existence of H designs

An H design is a triple , where is a set of points, a partition of into disjoint sets of size , and a set of ‐element transverses of , such that each ‐element transverse of is contained in exactly one of them. In 1990, Mills determined the existence of an H design with . In this paper, an efficient construction shows that an H exists for any integer with . Consequently, the necessary and sufficient conditions for the existence of an H design are , , and , with a definite exception .     

The maximum number of columns in supersaturated designs with

Cheng and Tang [Biometrika, 88 (2001), pp. 1169–1174] derived an upper bound on the maximum number of columns that can be accommodated in a two‐symbol supersaturated design (SSD) for a given number of rows () and a maximum in absolute value correlation between any two columns (). In particular, they proved that for (mod ) and . However, the only known SSD satisfying this upper bound is when . By utilizing a computer search, we prove that for , and . These results are obtained by proving the nonexistence of certain resolvable incomplete blocks designs. The combinatorial properties of the RIBDs are used to reduce the search space. Our results improve the lower bound for SSDs with rows and columns, for , and . Finally, we show that a skew‐type Hadamard matrix of order can be used to construct an SSD with rows and columns that proves . Hence, we establish for and for all (mod ) such that . Our result also implies that when is a prime power and (mod ). We conjecture that for all and (mod ), where is the maximum number of equiangular lines in with pairwise angle .     

On Turán problems for Cartesian products of graphs

Let be disjoint sets of sizes and . Let be a family of quadruples, having elements from and from , such that any subset with and contains one of the quadruples. We prove that the smallest size of is as . We also solve asymptotically a more general two‐partite Turán problem for quadruples.     

The short pulse equation: Bäcklund transformations and applications

A Bäcklund transformation (BT), which involves both independent and dependent variables, is established and studied for the short pulse (SP) equation. Based it, the nonlinear superposition formulae for 2-, 3-, and 4-BT are presented. The general result for the composition of -BTs is achieved and given in terms of determinants. As applications, various solutions including loop solitons, breather solutions, and their interaction solutions are worked out.     

Phase transitions in graphs on orientable surfaces

Let be the orientable surface of genus and denote by the class of all graphs on vertex set with edges embeddable on . We prove that the component structure of a graph chosen uniformly at random from features two phase transitions. The first phase transition mirrors the classical phase transition in the Erd?s‐Rényi random graph chosen uniformly at random from all graphs with vertex set and edges. It takes place at , when the giant component emerges. The second phase transition occurs at , when the giant component covers almost all vertices of the graph. This kind of phenomenon is strikingly different from and has only been observed for graphs on surfaces.     

Painlevé III Asymptotics of Hankel Determinants for a Perturbed Jacobi Weight

We study the Hankel determinants associated with the weight where , , , is analytic in a domain containing [ ? 1, 1] and for . In this paper, based on the Deift–Zhou nonlinear steepest descent analysis, we study the double scaling limit of the Hankel determinants as and . We obtain the asymptotic approximations of the Hankel determinants, evaluated in terms of the Jimbo–Miwa–Okamoto σ‐function for the Painlevé III equation. The asymptotics of the leading coefficients and the recurrence coefficients for the perturbed Jacobi polynomials are also obtained.