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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Inverse scattering transforms and soliton solutions of nonlocal reverse-space nonlinear Schrödinger hierarchies

The aim of the paper is to construct nonlocal reverse-space nonlinear Schrödinger (NLS) hierarchies through nonlocal group reductions of eigenvalue problems and generate their inverse scattering transforms and soliton solutions. The inverse scattering problems are formulated by Riemann-Hilbert problems which determine generalized matrix Jost eigenfunctions. The Sokhotski-Plemelj formula is used to transform the Riemann-Hilbert problems into Gelfand-Levitan-Marchenko type integral equations. A solution formulation to special Riemann-Hilbert problems with the identity jump matrix, corresponding to the reflectionless transforms, is presented and applied to -soliton solutions of the nonlocal NLS hierarchies.     

Well‐posedness of fractional integro‐differential equations in vector‐valued functional spaces

We study the well‐posedness of the fractional differential equations with infinite delay on Lebesgue–Bochner spaces and Besov spaces , where A and B are closed linear operators on a Banach space X satisfying ,  and . Under suitable assumptions on the kernels a and b, we completely characterize the well‐posedness of in the above vector‐valued function spaces on by using known operator‐valued Fourier multiplier theorems. We also give concrete examples where our abstract results may be applied.     

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 .     

Decomposable super‐simple NRBIBDs with block size 4 and index 6

Necessary conditions for the existence of a super‐simple, decomposable, near‐resolvable ‐balanced incomplete block design (BIBD) whose 2‐component subdesigns are both near‐resolvable ‐BIBDs are (mod ) and . In this paper, we show that these necessary conditions are sufficient. Using these designs, we also establish that the necessary conditions for the existence of a super‐simple near‐resolvable ‐RBIBD, namely (mod ) and , are sufficient. A few new pairwise balanced designs are also given.     

A note on the eigenvalues of p‐fractional Hardy–Sobolev operator with indefinite weight

In this article, we study the eigenvalues of p‐fractional Hardy operator where , , , and Ω is an unbounded domain in with Lipschitz boundary containing 0. The weight function V may change sign and may have singular points. We also show that the least positive eigenvalue is simple and it is uniquely associated to a nonnegative eigenfunction. Moreover, we proved that there exists a sequence of eigenvalues as .     

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 .     

‐Theory for Schrödinger systems

In this article, we study for the ‐realization of the vector‐valued Schrödinger operator . Using a noncommutative version of the Dore–Venni theorem due to Monniaux and Prüss, we prove that the ‐realization of , defined on the intersection of the natural domains of the differential and multiplication operators which form , generates a strongly continuous contraction semigroup on . We also study additional properties of the semigroup such as extension to L¹, positivity, ultracontractivity and prove that the generator has compact resolvent.