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 .     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Triangle-free planar graphs with the smallest independence number

Steinberg and Tovey proved that every -vertex planar triangle-free graph has an independent set of size at least , and described an infinite class of tight examples. We show that all -vertex planar triangle-free graphs except for this one infinite class have independent sets of size at least .     

Counting restricted orientations of random graphs

We count orientations of avoiding certain classes of oriented graphs. In particular, we study , the number of orientations of the binomial random graph in which every copy of is transitive, and , the number of orientations of containing no strongly connected copy of . We give the correct order of growth of and up to polylogarithmic factors; for orientations with no cyclic triangle, this significantly improves a result of Allen, Kohayakawa, Mota, and Parente. We also discuss the problem for a single forbidden oriented graph, and state a number of open problems and conjectures.     

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 .     

Extending edge-colorings of complete hypergraphs into regular colorings

Let be the collection of all -subsets of an -set . Given a coloring (partition) of a set , we are interested in finding conditions under which this coloring is extendible to a coloring of so that the number of times each element of appears in each color class (all sets of the same color) is the same number . The case was studied by Sylvester in the 18th century and remained open until the 1970s. The case is extensively studied in the literature and is closely related to completing partial symmetric Latin squares. For , we settle the cases , and completely. Moreover, we make partial progress toward solving the case where . These results can be seen as extensions of the famous Baranyai’s theorem, and make progress toward settling a 40-year-old problem posed by Cameron.     

Maximizing the mean subtree order

This article focuses on properties and structures of trees with maximum mean subtree order in a given family; such trees are called optimal in the family. Our main goal is to describe the structure of optimal trees in and , the families of all trees and caterpillars, respectively, of order . We begin by establishing a powerful tool called the Gluing Lemma, which is used to prove several of our main results. In particular, we show that if is an optimal tree in or for , then every leaf of is adjacent to a vertex of degree at least . We also use the Gluing Lemma to answer an open question of Jamison and to provide a conceptually simple proof of Jamison's result that the path has minimum mean subtree order among all trees of order . We prove that if is optimal in , then the number of leaves in is and that if is optimal in , then the number of leaves in is . Along the way, we describe the asymptotic structure of optimal trees in several narrower families of trees.     

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.