Search

We tackle the regularisation of a differential system related to generalised Krawtchouk polynomials. We show a straightforward connection between certain auxiliary quantities involving the recurrence coefficients of these polynomials and Painlev\'e V equation via iterative regularisation. Furthermore, we explore how iterative regularisation yields polynomial systems and enables us to find decompositions of certain birational transformations.

General Geronimus transformations, defined by regular matrix polynomials that are neither required to be monic nor restricted by the rank of their leading coefficients, are applied through both right and left multiplication to a rectangular matrix of measures associated with mixed multiple orthogonal polynomials. These transformations produce Christoffel-type formulas that establish relationships between the perturbed and original polynomials. Moreover, it is proven that the existence of Geronimus-perturbed orthogonality is equivalent to the non-cancellation of certain $\tau$-determinants. The effect of these transformations on the Markov-Stieltjes matrix functions is also determined. As a case study, we examine the Jacobi-Pi\~neiro orthogonal polynomials with three weights., Comment: 40 pages

Mixed orthogonal Laurent polynomials on the unit circle of CMV type are constructed utilizing a matrix of moments and its Gauss--Borel factorization and employing a multiple extension of the CMV ordering. A systematic analysis of the associated multiple orthogonality and biorthogonality relations, and an examination of the degrees of the Laurent polynomials is given. Recurrence relations, expressed in terms of banded matrices, are found. These recurrence relations lay the groundwork for corresponding Christoffel-Darboux kernels and relations, as well as for elucidating the ABC theorem. The paper also develops the theory of diagonal Christoffel and Geronimus perturbations of the matrix of measures. Christoffel formulas are found for both perturbations., Comment: 48 pages

This paper demonstrates how to explicitly construct a bidiagonal factorization of the banded recurrence matrix that appears in mixed multiple orthogonality on the step-line in terms of the coeffcients of the mixed multiple orthogonal polynomials. The construction is based on the \(LU\) factorization of the moment matrix and Christoffel transformations applied to the matrix of measures and the associated mixed multiple orthogonal polynomials., Comment: 14 pages

Understanding open quantum systems using information encoded in its complex eigenvalues has been a subject of growing interest. In this paper, we study higher-order gap ratios of the singular values of generic open quantum systems. We show that $k$-th order gap ratio of the singular values of an open quantum system can be connected to the nearest-neighbor spacing ratio of positions of classical particles of a harmonically confined log-gas with inverse temperature $\beta'(k)$ where $\beta'(k)$ is an analytical function that depends on $k$ and the Dyson's index $\beta=1,2,$ and $4$ that characterizes the properties of the associated Hermitized matrix. Our findings are crucial not only for understanding long-range correlations between the eigenvalues but also provide an excellent way of distinguishing different symmetry classes in an open quantum system. To highlight the universality of our findings, we demonstrate the higher-order gap ratios using different platforms such as non-Hermitian random matrices, random dissipative Liouvillians, Hamiltonians coupled to a Markovian bath, and Hamiltonians with in-built non-Hermiticity., Comment: 12 pages, 12 figures

This work explores classical discrete multiple orthogonal polynomials, including Hahn, Meixner of the first and second kinds, Kravchuk, and Charlier polynomials, with an arbitrary number of weights. Explicit expressions for the recursion coefficients of Hahn multiple orthogonal polynomials are derived. By leveraging the multiple Askey scheme and the recently discovered explicit hypergeometric representation of type I multiple Hahn polynomials, corresponding explicit hypergeometric representations are provided for the type I polynomials and recursion coefficients of all the aforementioned descendants within the Askey scheme. Additionally, integral representations for these families within the Hahn class in the Askey scheme are presented. The multiple Askey scheme is further completed by providing the corresponding limits for the weights, polynomials, and recurrence coefficients., Comment: 37 pages, 5 figures

This paper addresses two primary objectives in the realm of classical multiple orthogonal polynomials with an arbitrary number of weights. Firstly, it establishes new and explicit hypergeometric expressions for type I Hahn multiple orthogonal polynomials. Secondly, applying the residue theorem and the Mellin transform, the paper derives contour integral representations for several families of orthogonal polynomials. Specifically, it presents contour integral formulas for both type I and type II multiple orthogonal polynomials in the Laguerre of the first kind, Jacobi-Pi\~neiro, and Hahn families. The evaluation of these integrals leads to explicit hypergeometric representations., Comment: 21 pages

In this paper, we extend our investigation into semiclassical multiple discrete orthogonal polynomials by considering an arbitrary number of weights. We derive multiple versions of the Toda equations and the Laguerre-Freud equations for the multiple generalized Charlier and multiple generalized Meixner II families., Comment: 23 pages

Performing both right and left multiplication operations using general regular matrix polynomials, which need not be monic and may possess leading coefficients of arbitrary rank, on a rectangular matrix of measures associated with mixed multiple orthogonal polynomials, reveals corresponding Christoffel formulas. These formulas express the perturbed mixed multiple orthogonal polynomials in relation to the original ones. Utilizing the divisibility theorem for matrix polynomials, we establish a criterion for the existence of perturbed orthogonality, expressed through the non-cancellation of certain $\tau$ determinants., Comment: 32 pages. Some minor typos corrected in the second version

This paper serves as an introduction to banded totally positive matrices, exploring various characterizations and associated properties. A significant result within is the demonstration that the collection of such matrices forms a semigroup, notably including a subset permitting positive bidiagonal factorization. Moreover, the paper applies this concept to investigate step line normality concerning the degrees of associated recursion polynomials. It presents a spectral Favard theorem, ensuring the existence of measures, thereby guaranteeing that these recursion polynomials represent mixed multiple orthogonal polynomials that maintain normality on the step line indices.

This paper delves into classical multiple orthogonal polynomials with an arbitrary number of weights, including Jacobi-Pi\~neiro, Laguerre of both first and second kinds, as well as multiple orthogonal Hermite polynomials. Novel explicit expressions for nearest-neighbor recurrence coefficients, as well as the step line case, are provided for all these polynomial families. Furthermore, new explicit expressions for type I multiple orthogonal polynomials are derived for Laguerre of the second kind and also for multiple Hermite polynomials.

We investigate the full counting statistics (FCS) of a harmonically confined 1d short-range Riesz gas consisting of $N$ particles in equilibrium at finite temperature. The particles interact with each other through a repulsive power-law interaction with an exponent $k>1$ which includes the Calogero-Moser model for $k=2$. We examine the probability distribution of the number of particles in a finite domain $[-W, W]$ called number distribution, denoted by $\mathcal{N}(W, N)$. We analyze the probability distribution of $\mathcal{N}(W, N)$ and show that it exhibits a large deviation form for large $N$ characterised by a speed $N^{\frac{3k+2}{k+2}}$ and by a large deviation function of the fraction $c = \mathcal{N}(W, N)/N$ of the particles inside the domain and $W$. We show that the density profiles that create the large deviations display interesting shape transitions as one varies $c$ and $W$. This is manifested by a third-order phase transition exhibited by the large deviation function that has discontinuous third derivatives. Monte-Carlo (MC) simulations show good agreement with our analytical expressions for the corresponding density profiles. We find that the typical fluctuations of $\mathcal{N}(W, N)$, obtained from our field theoretic calculations are Gaussian distributed with a variance that scales as $N^{\nu_k}$, with $\nu_k = (2-k)/(2+k)$. We also present some numerical findings on the mean and the variance. Furthermore, we adapt our formalism to study the index distribution (where the domain is semi-infinite $(-\infty, W])$, linear statistics (the variance), thermodynamic pressure and bulk modulus., Comment: 36 pages, 7 figures

In our contribution we study stochastic models in one space dimension with two conservation laws. One model is the coupled continuum stochastic Burgers equation, for which each current is a sum of quadratic non-linearities, linear diffusion, and spacetime white noise. The second model is a two-lane stochastic lattice gas. As distinct from previous studies, the two conserved densities are tuned such that the flux Jacobian, a $2 \times 2$ matrix, has coinciding eigenvalues. In the steady state, investigated are spacetime correlations of the conserved fields and the time-integrated currents at the origin. For a particular choice of couplings the dynamical exponent 3/2 is confirmed. Furthermore, at these couplings, continuum stochastic Burgers equation and lattice gas are demonstrated to be in the same universality class., Comment: 34 pages, 9 figures

This paper concerns a boundary integral formulation for the two-dimensional massive Dirac equation. The mass term is assumed to jump across a one-dimensional interface, which models a transition between two insulating materials. This jump induces surface waves that propagate outward along the interface but decay exponentially in the transverse direction. After providing a derivation of our integral equation, we prove that it has a unique solution for almost all choices of parameters using holomorphic perturbation theory. We then extend these results to a Dirac equation with two interfaces. Finally, we implement a fast numerical method for solving our boundary integral equations and present several numerical examples of solutions and scattering effects.

This paper explores a factorization using bidiagonal matrices of the recurrence matrix of Hahn multiple orthogonal polynomials. The factorization is expressed in terms of ratios involving the generalized hypergeometric function ${}_3F_2$ and is proven using recently discovered contiguous relations. Moreover, employing the multiple Askey scheme, a bidiagonal factorization is derived for the Hahn descendants, including Jacobi-Pi\~neiro, multiple Meixner (kinds I and II), multiple Laguerre (kinds I and II), multiple Kravchuk, and multiple Charlier, all represented in terms of hypergeometric functions. For the cases of multiple Hahn, Jacobi-Pi\~neiro, Meixner of kind II, and Laguerre of kind I, where there exists a region where the recurrence matrix is nonnegative, subregions are identified where the bidiagonal factorization becomes a positive bidiagonal factorization., Comment: 14 pages, 2 figures

This paper investigates stochastic finite matrices and the corresponding finite Markov chains constructed using recurrence matrices for general families of orthogonal polynomials and multiple orthogonal polynomials. The paper explores the spectral theory of transition matrices, utilizing both orthogonal and multiple orthogonal polynomials. Several properties are derived, including classes, periodicity, recurrence, stationary states, ergodicity, expected recurrence times, time-reversed chains, and reversibility. Furthermore, the paper uncovers factorization in terms of pure birth and pure death processes. The case study focuses on hypergeometric orthogonal polynomials, where all the computations can be carried out effectively. Particularly within the Askey scheme, all descendants under Hahn (excluding Bessel), such as Hahn, Jacobi, Meixner, Kravchuk, Laguerre, Charlier, and Hermite, present interesting examples of recurrent reversible birth and death finite Markov chains. Additionally, the paper considers multiple orthogonal polynomials, including multiple Hahn, Jacobi-Pi\~neiro, Laguerre of the first kind, and Meixner of the second kind, along with their hypergeometric representations and derives the corresponding recurrent finite Markov chains and time-reversed chains., Comment: 40 pages, 3 figures

In this paper, the authors investigate the case of discrete multiple orthogonal polynomials with two weights on the step line, which satisfy Pearson equations. The discrete multiple orthogonal polynomials in question are expressed in terms of tau functions, which are double Wronskians of generalized hypergeometric series. The shifts in the spectral parameter for type II and type I multiple orthogonal polynomials are described using banded matrices. It is demonstrated that these polynomials offer solutions to multicomponent integrable extensions of the nonlinear Toda equations. Additionally, the paper characterizes extensions of the Nijhoff-Capel totally discrete Toda equations. The hypergeometric $\tau$-functions are shown to provide solutions to these integrable nonlinear equations. Furthermore, the authors explore Laguerre-Freud equations, nonlinear equations for the recursion coefficients, with a particular focus on the multiple Charlier, generalized multiple Charlier, multiple Meixner II, and generalized multiple Meixner II cases., Comment: 31 pages

Anomalous transport in low dimensional spin chains is an intriguing topic that can offer key insights into the interplay of integrability and symmetry in many-body dynamics. Recent studies have shown that spin-spin correlations in spin chains, where integrability is either perfectly preserved or broken by symmetry-preserving interactions, fall in the Kardar-Parisi-Zhang (KPZ) universality class. Similarly, energy transport can show ballistic or diffusive-like behaviour. Although such behaviour has been studied under equilibrium conditions, no results on nonequilibrium spin transport in classical spin chains has been reported so far. In this work, we investigate both spin and energy transport in classical spin chains (integrable and non-integrable) when coupled to two reservoirs at two different temperatures/magnetization. In both the integrable case and broken-integrability (but spin-symmetry preserving), we report anomalous scaling of spin current with system size ($\mathbb{J}^s \propto L^{-\mu}$) with an exponent, $\mu \approx 2/3$, falling under the KPZ universality class. On the other hand, it is noteworthy that energy current remains ballistic ($\mathbb{J}^e \propto L^{-\eta}$ with $\eta \approx 0$) in the purely integrable case and there is departure from ballistic behaviour ($\eta > 0$) when integrability is broken regardless of spin-symmetry. Under nonequilibrium conditions, we have thoroughly investigated spatial profiles of local magnetization and energy. We find interesting nonlinear spatial profiles which are hallmarks of anomalous transport. We also unravel subtle striking differences between the equilibrium and nonequilibrium steady state through the lens of spatial spin-spin correlations., Comment: 13 pages, 10 figures (including supplementary material)

The out-of-time-ordered correlator (OTOC) has emerged as an interesting object in both classical and quantum systems for probing the spatial spread and temporal growth of initially local perturbations in spatially extended chaotic systems. Here, we study the (classical) OTOC and its ``light-cone'' in the nonlinear Kuramoto-Sivashinsky (KS) equation, using extensive numerical simulations. We also show that the linearized KS equation exhibits a qualitatively similar OTOC and light-cone, which can be understood via a saddle-point analysis of the linearly unstable modes. Given the deep connection between the KS (deterministic) and the Kardar-Parisi-Zhang (KPZ which is stochastic) equations, we also explore the OTOC in the KPZ equation. While our numerical results in the KS case are expected to hold in the continuum limit, for the KPZ case it is valid in a discretized version of the KPZ equation. More broadly, our work unravels the intrinsic interplay between noise/instability, nonlinearity and dissipation in partial differential equations (deterministic or stochastic) through the lens of OTOC., Comment: 8 pages (including supplemental material), 5 figures, updated

We investigate the possibility of a Markovian quantum master equation (QME) that consistently describes a finite-dimensional system, a part of which is weakly coupled to a thermal bath. In order to preserve complete positivity and trace, such a QME must be of Lindblad form. For physical consistency, it should additionally preserve local conservation laws and be able to show thermalization. We search of Lindblad equations satisfying these additional criteria. First, we show that the microscopically derived Bloch-Redfield equation (RE) violates complete positivity unless in extremely special cases. We then prove that imposing complete positivity and demanding preservation of local conservation laws enforces the Lindblad operators and the lamb-shift Hamiltonian to be `local', i.e, to be supported only on the part of the system directly coupled to the bath. We then cast the problem of finding `local' Lindblad QME which can show thermalization into a semidefinite program (SDP). We call this the thermalization optimization problem (TOP). For given system parameters and temperature, the solution of the TOP conclusively shows whether the desired type of QME is possible up to a given precision. Whenever possible, it also outputs a form for such a QME. For a XXZ chain of few qubits, fixing a reasonably high precision, we find that such a QME is impossible over a considerably wide parameter regime when only the first qubit is coupled to the bath. Remarkably, we find that when the first two qubits are attached to the bath, such a QME becomes possible over much of the same paramater regime, including a wide range of temperatures.

We consider the analysis of singular waveguides separating insulating phases in two-space dimensions. The insulating domains are modeled by a massive Schr\"odinger equation and the singular waveguide by appropriate jump conditions along the one-dimensional interface separating the insulators. We present an integral formulation of the problem and analyze its mathematical properties. We also implement a fast multipole and sweeping-accelerated iterative algorithm for solving the integral equations, and demonstrate numerically the fast convergence of this method. Several numerical examples of solutions and scattering effects illustrate our theory.

Explicit expressions for the Hahn multiple polynomials of type I, in terms of Kamp\'e de F\'eriet hypergeometric series, are given. Orthogonal and biorthogonal relations are proven. Then, part of the Askey scheme for multiple orthogonal polynomials type I is completed. In particular, explicit expressions in terms of generalized hypergeometric series and Kamp\'e de F\'eriet hypergeometric series, are given for the multiple orthogonal polynomials of type I for the Jacobi-Pi\~neiro, Meixner I, Meixner II, Kravchuk, Laguerre I, Laguerre II and Charlier families., Comment: 27 pages. Minor typos corrected

Spectral and factorization properties of oscillatory matrices leads to a spectral Favard theorem for bounded banded matrices, that admit a positive bidiagonal factorization, in terms of sequences of mixed multiple orthogonal polynomials with respect to a set positive Lebesgue-Stieltjes measures. A mixed multiple Gauss quadrature formula with corresponding degrees of precision is given., Comment: 33 pages

The Cholesky factorization of the moment matrix is considered for discrete orthogonal polynomials of hypergeometrical type. We derive the Laguerre-Freud equations when the first moments of the weights are given by the ${}_1F_2$, ${}_2F_2$ and ${}_3F_2$ generalized hypergeometrical functions., Comment: 19 pages. arXiv admin note: substantial text overlap with arXiv:2107.02177, arXiv:2107.01747

We study the average density of $N$ spinless noninteracting fermions in a $2d$ harmonic trap rotating with a constant frequency $\Omega$ and in the presence of an additional repulsive central potential $\gamma/r^2$. The average density at zero temperature was recently studied in Phys. Rev. A $\textbf{103}$, 033321 (2021) and an interesting multi-layered "wedding cake" structure with a "hole" at the center was found for the density in the large $N$ limit. In this paper, we study the average density at finite temperature. We demonstrate how this "wedding-cake" structure is modified at finite temperature. These large $N$ results warrant going much beyond the standard Local Density Approximation. We also generalize our results to a wide variety of trapping potentials and demonstrate the universality of the associated scaling functions both in the bulk and at the edges of the "wedding-cake"., Comment: 18 pages, 7 figures, 1 table

Recent investigations have observed superdiffusion in integrable classical and quantum spin chains. An intriguing connection between these spin chains and Kardar-Parisi-Zhang (KPZ) universality class has emerged. Theoretical developments (e.g. generalized hydrodynamics) have highlighted the role of integrability as well as spin-symmetry in KPZ behaviour. However understanding their precise role on superdiffusive transport still remains a challenging task. The widely used quantum spin chain platform comes with severe numerical limitations. To circumvent this barrier, we focus on a classical integrable spin chain which was shown to have deep analogy with the quantum spin-$\frac{1}{2}$ Heisenberg chain. Remarkably, we find that KPZ behaviour prevails even when one considers integrability-breaking but spin-symmetry preserving terms, strongly indicating that spin-symmetry plays a central role even in the non-perturbative regime. On the other hand, in the non-perturbative regime, we find that energy correlations exhibit clear diffusive behaviour. We also study the classical analog of out-of-time-ordered correlator (OTOC) and Lyapunov exponents. We find significant presence of chaos for the integrability-broken cases even though KPZ behaviour remains robust. The robustness of KPZ behaviour is demonstrated for a wide class of spin-symmetry preserving integrability-breaking terms., Comment: 10 pages, 9 figures (including supplementary material)

A spectral Favard theorem for bounded banded lower Hessenberg matrices that admit a positive bidiagonal factorization is found. The large knowledge on the spectral and factorization properties of oscillatory matrices leads to this spectral Favard theorem in terms of sequences of multiple orthogonal polynomials of types I and II with respect to a set of positive Lebesgue-Stieltjes~measures. Also a multiple Gauss quadrature is proven and corresponding degrees of precision are found. This spectral Favard theorem is applied to Markov chains with $(p+2)$-diagonal transition matrices, i.e. beyond birth and death, that admit a positive stochastic bidiagonal factorization. In the finite case, the Karlin-McGregor spectral representation is given. It is shown that the Markov chains are recurrent and explicit expressions in terms of the orthogonal polynomials for the stationary distributions are given. Similar results are obtained for the countable infinite Markov chain. Now the Markov chain is not necessarily recurrent, and it is characterized in terms of the first measure. Ergodicity of the Markov chain is discussed in terms of the existence of a mass at $1$, which is an eigenvalue corresponding to the right and left eigenvectors., Comment: 38 pages. This a very much improved version of the initial one. We have divided the original paper in three parts. Two more parts to follow: arXiv:2210.10728, Positive bidiagonal factorization of tetradiagonal Hessenberg matrices arXiv:2210.10727, Bidiagonal factorization of tetradiagonal matrices and Darboux transformations

We study the distribution of the position of the rightmost particle $x_{\max}$ in a $N$-particle Riesz gas in one dimension confined in a harmonic trap. The particles interact via long-range repulsive potential, of the form $r^{-k}$ with $-2-2$. We also find that these large deviation functions describe a pulled to pushed type phase transition as observed in Dyson's log-gas ($k\to 0$) and $1d$ one component plasma ($k=-1$). Remarkably, we find that the phase transition remains $3^{\rm rd}$ order for the entire regime. Our results demonstrate the striking universality of the $3^{\rm rd}$ order transition even in models that fall outside the paradigm of Coulomb systems and the random matrix theory. We numerically verify our analytical expressions of the large deviation functions via Monte Carlo simulation using an importance sampling algorithm., Comment: 42 pages, 13 figures

In this paper we analyze the dynamics of single-excitation states, which model the scattering of a single photon from multiple two level atoms. For short times and weak atom-field couplings we show that the atomic amplitudes are given by a sum of decaying exponentials, where the decay rates and Lamb shifts are given by the poles of a certain analytic function. This result is a refinement of the "pole approximation" appearing in the standard Wigner-Weisskopf analysis of spontaneous emission. On the other hand, at large times, the atomic field decays like $O(1/t^3)$ with a known constant expressed in terms of the coupling parameter and the resonant frequency of the atoms. Moreover, we show that for stronger coupling, the solutions also feature a collection of oscillatory exponentials which dominate the behavior at long times. Finally, we extend the analysis to the continuum limit in which atoms are distributed according to a given density., Comment: 22 pages, 4 figures

We consider the $N$ particle classical Riesz gas confined in a one-dimensional external harmonic potential with power law interaction of the form $1/r^k$ where $r$ is the separation between particles. As special limits it contains several systems such as Dyson's log-gas ($k\to 0^+$), Calogero-Moser model ($k=2$), 1d one component plasma ($k=-1$) and the hard-rod gas ($k\to \infty$). Despite its growing importance, only large-$N$ field theory and average density profile are known for general $k$. In this Letter, we study the fluctuations in the system by looking at the statistics of the gap between successive particles. This quantity is analogous to the well-known level spacing statistics which is ubiquitous in several branches of physics. We show that the variance goes as $N^{-b_k}$ and we find the $k$ dependence of $b_k$ via direct Monte Carlo simulations. We provide supporting arguments based on microscopic Hessian calculation and a quadratic field theory approach. We compute the gap distribution and study its system size scaling. Except in the range $-1-2$ with both Gaussian and non-Gaussian scaling forms., Comment: 13 pages, 11 figures

Contiguous hypergeometric relations for semiclassical discrete orthogonal polynomials are described as Christoffel and Geronimus transformations. Using the Christoffel-Geronimus-Uvarov formulas quasi-determinatal expressions for the shifted semiclassical discrete orthogonal polynomials are obtained.

The Cholesky factorization of the moment matrix is considered for the generalized Charlier, generalized Meixner, and Gauss hypergeometric discrete orthogonal polynomials. For the generalized Charlier, we present an alternative derivation of the Laguerre-Freud relations found by Smet and Van Assche. Third-order and second-order nonlinear ordinary differential equations are found for the recursion coefficient $\gamma_n$, that happen to be forms of the Painlev\'e $\text{deg-P}_{\text V}$ in disguise. Laguerre-Freud relations are also found for the generalized Meixner case, which are compared with those of Smet and Van Assche. Finally, the Gauss hypergeometric discrete orthogonal polynomials, also known as generalized Hahn of type I, are also studied. Laguerre-Freud equations are found, and the differences with the equations found by Dominici and by Filipuk and Van Assche are provided., Comment: Some minor typos corrected In the new version the bibliography has been updated and the text improved. Typo in Theorem 8 corrected. arXiv admin note: substantial text overlap with arXiv:2107.01747

The Cholesky factorization of the moment matrix is applied to discrete orthogonal polynomials on the homogeneous lattice. In particular, semiclassical discrete orthogonal polynomials, which are built in terms of a discrete Pearson equation, are studied. The Laguerre-Freud structure semi-infinite matrix that models the shifts by $\pm 1$ in the independent variable of the set of orthogonal polynomials is introduced. In the semiclassical case it is proven that this Laguerre-Freud matrix is banded. From the well known fact that moments of the semiclassical weights are logarithmic derivatives of generalized hypergeometric functions, it is shown how the contiguous relations for these hypergeometric functions translate as symmetries for the corresponding moment matrix. It is found that the 3D Nijhoff-Capel discrete Toda lattice describes the corresponding contiguous shifts for the squared norms of the orthogonal polynomials. The continuous Toda for these semiclassical discrete orthogonal polynomials is discussed and the compatibility equations are derived. It also shown that the Kadomtesev-Petvishvilii equation is connected to an adequate deformed semiclassical discrete weight, but in this case the deformation do not satisfy a Pearson equation.

In this paper, we compute exactly the average density of a harmonically confined Riesz gas of $N$ particles for large $N$ in the presence of a hard wall. In this Riesz gas, the particles repel each other via a pairwise interaction that behaves as $|x_i - x_j|^{-k}$ for $k>-2$, with $x_i$ denoting the position of the $i^{\rm th}$ particle. This density can be classified into three different regimes of $k$. For $k \geq 1$, where the interactions are effectively short-ranged, the appropriately scaled density has a finite support over $[-l_k(w),w]$ where $w$ is the scaled position of the wall. While the density vanishes at the left edge of the support, it approaches a nonzero constant at the right edge $w$. For $-1

The recently found hypergeometric multiple orthogonal polynomials on the step-line by Lima and Loureiro are shown to be random walk polynomials. It is proven that the corresponding Jacobi matrix and its transpose, which are nonnegative matrices and describe higher recurrence relations, can be normalized to two stochastic matrices, dual to each other. Using the Christoffel-Darboux formula on the step-line and the Poincar\'e theory for non-homogeneous recurrence relations it is proven that both stochastic matrices are related by transposition in the large $n$ limit. These random walks are beyond birth and death, as they describe a chain in where transitions to the two previous states are allowed, or in the dual to the two next states.The corresponding Karlin-McGregor representation formula is given for these new Markov chains. The regions of hypergeometric parameters where the Markov chains are recurrent or transient are given. Stochastic factorizations, in terms of pure birth and of pure death factors, for the corresponding Markov matrices of types I and II, are provided.Twelve uniform Jacobi matrices and the corresponding random walks, related to a Jacobi matrix of Toeplitz type, and theirs stochastic or semi-stochastic matrices (with sinks and sources), that describe Markov chains beyond birth and death, are found and studied. One of these uniform stochastic cases, which is a recurrent random walk, is the only hypergeometric multiple random walk having a uniform stochastic factorization. The corresponding weights, Jacobi and Markov transition matrices and sequences of type II multiple orthogonal polynomials are provided. Chain of Christoffel transformations connecting the stochastic uniform tuples between them, and the semi-stochastic uniform tuples, between them, are presented., Comment: In the revision we added explicit expressions for the type I multiple orthogonal polynomials associated with the uniform tuples

Multiple orthogonal polynomials with respect to two weights on the step-line are considered. A connection between different dual spectral matrices, one banded (recursion matrix) and one Hessenberg, respectively, and the Gauss-Borel factorization of the moment matrix is given. It is shown a hidden freedom exhibited by the spectral system related to the multiple orthogonal polynomials. Pearson equations are discussed, a Laguerre-Freud matrix is considered, and differential equations for type I and II multiple orthogonal polynomials, as well as for the corresponding linear forms are given. The Jacobi-Pi\~neiro multiple orthogonal polynomials of type I and type II are used as an illustrating case and the corresponding differential relations are presented. A permuting Christoffel transformation is discussed, finding the connection between the different families of multiple orthogonal polynomials. The Jacobi-Pi\~neiro case provides a convenient illustration of these symmetries, giving linear relations between different polynomials with shifted and permuted parameters. We also present the general theory for the perturbation of each weight by a different polynomial or rational function aka called Christoffel and Geronimus transformations. The connections formulas between the type II multiple orthogonal polynomials, the type I linear forms, as well as the vector Stieltjes-Markov vector functions is also presented. We illustrate these findings by analyzing the special case of modification by an even polynomial., Comment: Revised version, completely new section on general Christoffel and Geronimus for multiple orthogonal polynomials on the stepline

Given a non-negative Jacobi matrix describing higher order recurrence relations for multiple orthogonal polynomials of type~II and corresponding linear forms of type I, a general strategy for constructing a pair of stochastic matrices, dual to each other, is provided. The corresponding Markov chains (or 1D random walks) allow, in one transition, to reach for the N-th previous states, to remain in the state or reach for the immediately next state. The dual Markov chains allow, in one transition, to reach for the N-th next states, to remain in the state or reach for immediately previous state. The connection between both dual Markov chains is discussed at the light of the Poincar\'e's theorem on ratio asymptotics for homogeneous linear recurrence relations and the Christoffel-Darboux formula within the sequence of multiple orthogonal polynomials and linear forms of type I. The Karlin-McGregor representation formula is extended to both dual random walks, and applied to the discussion of the corresponding generating functions and first-passage distributions. Recurrent or transient character of the Markov chain is discussed. Steady state and some conjectures on its existence and the relation with mass points are also given. The Jacobi-Pi\~neiro multiple orthogonal polynomials are taken as a case study of the described results. For the first time in the literature, an explicit formula for the type~I Jacobi--Pi\~neiro polynomials is determined. The region of parameters where the Markov chains are recurrent or transient is given, and it is conjectured that when recurrent, the Markov chains are null recurrent and, consequently, the expected return times are infinity. Examples of recurrent and transient Jacobi--Pi\~neiro random walks are constructed explicitly.

We present a general scheme for constructing robust excitations (soliton-like) in non-integrable multicomponent systems. By robust, we mean localised excitations that propagate with almost constant velocity and which interact cleanly with little to no radiation. We achieve this via a reduction of these complex systems to more familiar effective chiral field-theories using perturbation techniques and the Fredholm alternative. As a specific platform, we consider the generalised multicomponent Nonlinear Schr\"{o}dinger Equations (MNLS) with arbitrary interaction coefficients. This non-integrable system reduces to uncoupled Korteweg-de Vries (KdV) equations, one for each sound speed of the system. This reduction then enables us to exploit the multi-soliton solutions of the KdV equation which in turn leads to the construction of soliton-like profiles for the original non-integrable system. We demonstrate that this powerful technique leads to the coherent evolution of excitations with minimal radiative loss in arbitrary non-integrable systems. These constructed coherent objects for non-integrable systems bear remarkably close resemblance to true solitons of integrable models. Although we use the ubiquitous MNLS system as a platform, our findings are a major step forward towards constructing excitations in generic continuum non-integrable systems., Comment: 21 pages, 3 figures

We present exact results for the classical version of the Out-of-Time-Order Commutator (OTOC) for a family of power-law models consisting of $N$ particles in one dimension and confined by an external harmonic potential. These particles are interacting via power-law interaction of the form $\propto \sum_{\substack{i, j=1 (i\neq j)}}^N|x_i-x_j|^{-k}$ $\forall$ $k>1$ where $x_i$ is the position of the $i^\text{th}$ particle. We present numerical results for the OTOC for finite $N$ at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many body ground state. In the large-$N$ limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the "light cone". We find remarkably distinct features in OTOC below and above $k=3$ in terms of going from non-Airy behaviour ($13$). We present certain additional rich features for the case $k=2$ that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatio-temporal spread of perturbations in long-range interacting systems., Comment: 14 pages, 4 figures (including supplementary material)

We compute exactly the average spatial density for $N$ spinless noninteracting fermions in a $2d$ harmonic trap rotating with a constant frequency $\Omega$ in the presence of an additional repulsive central potential $\gamma/r^2$. We find that, in the large $N$ limit, the bulk density has a rich and nontrivial profile -- with a hole at the center of the trap and surrounded by a multi-layered "wedding cake" structure. The number of layers depends on $N$ and on the two parameters $\Omega$ and $\gamma$ leading to a rich phase diagram. Zooming in on the edge of the $k^{\rm th}$ layer, we find that the edge density profile exhibits $k$ kinks located at the zeroes of the $k^{\rm th}$ Hermite polynomial. Interestingly, in the large $k$ limit, we show that the edge density profile approaches a limiting form, which resembles the shape of a propagating front, found in the unitary evolution of certain quantum spin chains. We also study how a newly formed droplet grows in size on top of the last layer as one changes the parameters., Comment: 17 pages, 10 figures. Published version, with typos corrected

Metrics and pseudometrics are defined on the group of unitary operators in a Hilbert space. Several explicit formulas are derived. A special feature of the work is investigation of pseudometrics in unitary groups. The rich classes of pseudometrics have many interesting applications. Three such applications, distinguishibility of unitary operators, quantum coding and quantum search problems are discussed., Comment: 18 pages, 1 figure

We study the dynamics of multi-component Bose gas described by the Vector Nonlinear Schr\"{o}dinger Equation (VNLS), aka the Vector Gross--Pitaevskii Equation (VGPE) . Through a Madelung transformation, the VNLS can be reduced to coupled hydrodynamic equations in terms of multiple density and velocity fields. Using a multi-scaling and a perturbation method along with the Fredholm alternative, we reduce the problem to a Korteweg de-Vries (KdV) system. This is of great importance to study more transparently, the obscure features hidden in VNLS. This ensures that hydrodynamic effects such as dispersion and nonlinearity are captured at an equal footing. Importantly, before studying the KdV connection, we provide a rigorous analysis of the linear problem. We write down a set of theorems along with proofs and associated corollaries that shine light on the conditions of existence and nature of eigenvalues and eigenvectors of the linear problem. This rigorous analysis is paramount for understanding the nonlinear problem and the KdV connection. We provide strong evidence of agreement between VNLS systems and KdV equations by using soliton solutions as a platform for comparison. Our results are expected to be relevant not only for cold atomic gases, but also for nonlinear optics and other branches where VNLS equations play a defining role., Comment: 21 pages, 2 Figures

We study an interacting system of $N$ classical particles on a line at thermal equilibrium. The particles are confined by a harmonic trap and repelling each other via pairwise interaction potential that behaves as a power law $\propto \sum_{\substack{i\neq j}}^N|x_i-x_j|^{-k}$ (with $k>-2$) of their mutual distance. This is a generalization of the well known cases of the one component plasma ($k=-1$), Dyson's log-gas ($k\to 0^+$), and the Calogero-Moser model ($k=2$). Due to the competition between harmonic confinement and pairwise repulsion, the particles spread over a finite region of space for all $k>-2$. We compute exactly the average density profile for large $N$ for all $k>-2$ and show that while it is independent of temperature for sufficiently low temperature, it has a rich and nontrivial dependence on $k$ with distinct behavior for $-21$ and $k=1$., Comment: Main text: 6 pages + 1 Fig., Supp. Mat.: 11 pages + 3 Figs. Accepted for publication in Physical Review Letters

The Gauss-Borel or $LU$ factorization of Gram matrices of bilinear forms is the pivotal element in the discussion of the theory of biorthogonal polynomials. The construction of biorthogonal families of polynomials and its second kind functions, of the spectral matrices modeling the multiplication by the independent variable $x$, the Christoffel-Darboux kernel and its projection properties, are discussed from this point of view. Then, the Hankel case is presented and different properties, specific of this case, as the three terms relations, Heine formulas, Gauss quadrature and the Christoffel-Darboux formula are given. The classical orthogonal polynomial of Hermite, Laguerre and Jacobi type are discussed and characterized within this scheme. Finally, it is shown who this approach is instrumental in the derivation of Christoffel formulas for general Christoffel and Geronimus perturbations of the bilinear forms., Comment: 23 pages, 2 figures

In this paper the Riemann-Hilbert problem, with jump supported on a appropriate curve on the complex plane with a finite endpoint at the origin, is used for the study of corresponding matrix biorthogonal polynomials associated with Laguerre type matrices of weights ---which are constructed in terms of a given matrix Pearson equation. First and second order differential systems for the fundamental matrix, solution of the mentioned Riemann-Hilbert problem are derived. An explicit and general example is presented to illustrate the theoretical results of the work. Related matrix eigenvalue problems for second order matrix differential operators and non-Abelian extensions of a family of discrete Painlev\'e IV equations are discussed., Comment: arXiv admin note: text overlap with arXiv:1807.07119

Models that remain integrable even in confining potentials are extremely rare and almost non-existent. Here, we consider a one-dimensional hyperbolic interaction model, which we call as the Hyperbolic Calogero (HC) model. This is classically integrable even in confining potentials (which have box-like shapes). We present a first-order formulation of the HC model in an external confining potential. Using the rich property of duality, we find multi-soliton solutions of this confined integrable model. Absence of solitons correspond to the equilibrium solution of the model. We demonstrate the dynamics of multi-soliton solutions via brute-force numerical simulations. We studied the physics of soliton collisions and quenches using numerical simulations. We have examined the motion of dual complex variables and found an analytic expression for the time period in a certain limit. We give the field theory description of this model and find the background solution (absence of solitons) analytically in the large-N limit. Analytical expressions of soliton solutions have been obtained in the absence of external confining potential. Our work is of importance to understand the general features of trapped interacting particles that remain classically integrable and can be of relevance to the collective behaviour of trapped cold atomic gases as well., Comment: 46 pages, 14 figures

In this work we discuss connections between a one-dimensional system of $N$ particles interacting with a repulsive inverse square potential and confined in a harmonic potential (Calogero-Moser model) and the log-gas model which appears in random matrix theory. Both models have the same minimum energy configuration, with the particle positions given by the zeros of the Hermite polynomial. Moreover, the Hessian describing small oscillations around equilibrium are also related for the two models. The Hessian matrix of the Calogero-Moser model is the square of that of the log-gas. We explore this connection further by studying finite temperature equilibrium properties of the two models through Monte-Carlo simulations. In particular, we study the single particle distribution and the marginal distribution of the boundary particle which, for the log-gas, are respectively given by the Wigner semi-circle and the Tracy-Widom distribution. For particles in the bulk, where typical fluctuations are Gaussian, we find that numerical results obtained from small oscillation theory are in very good agreement with the Monte-Carlo simulation results for both the models. For the log-gas, our findings agree with rigorous results from random matrix theory., Comment: 23 pages, 11 figures

In this paper we use the Riemann-Hilbert problem, with jumps supported on appropriate curves in the complex plane, for matrix biorthogonal polynomials and apply it to find Sylvester systems of differential equations for the orthogonal polynomials and its second kind functions as well. For this aim, Sylvester type differential Pearson equations for the matrix of weights are shown to be instrumental. Several applications are given, in order of increasing complexity. First, a general discussion of non-Abelian Hermite biorthogonal polynomials in the real line, understood as those whose matrix of weights is a solution of a Sylvester type Pearson equation with coefficients first order matrix polynomials, is given. All these is applied to the discussion of possible scenarios leading to eigenvalue problems for second order linear differential operators with matrix eigenvalues. Nonlinear matrix difference equations are discussed next. Firstly, for the general Hermite situation a general non linear relation (non trivial because the non commutativity features of the setting) for the recursion coefficients is gotten. In the next case of higher difficulty, degree two polynomials are allowed in the Pearson equation, but the discussion is simplified by considering only a left Pearson equation. In the case, the support of the measure is on an appropriate branch of an hyperbola. The recursion coefficients are shown to fulfill a non-Abelian extension of the alternate discrete Painlev\'e I equation. Finally, a discussion is given for the case of degree three polynomials as coefficients in the left Pearson equation characterizing the matrix of weights. However, for simplicity only odd polynomials are allowed. In this case, a new and more general matrix extension of the discrete Painlev\'e I equation is found., Comment: 38 pages, 1 figure

In Stokes flow, the stream function associated with the velocity of the fluid satisfies the biharmonic equation. The detailed behavior of solutions to the biharmonic equation on regions with corners has been historically difficult to characterize. The problem was first examined by Lord Rayleigh in 1920; in 1973, the existence of infinite oscillations in the domain Green's function was proven in the case of the right angle by S.~Osher. In this paper, we observe that, when the biharmonic equation is formulated as a boundary integral equation, the solutions are representable by rapidly convergent series of the form $\sum_{j} ( c_{j} t^{\mu_{j}} \sin{(\beta_{j} \log{(t)})} + d_{j} t^{\mu_{j}} \cos{(\beta_{j} \log{(t)})} )$, where $t$ is the distance from the corner and the parameters $\mu_{j},\beta_{j}$ are real, and are determined via an explicit formula depending on the angle at the corner. In addition to being analytically perspicuous, these representations lend themselves to the construction of highly accurate and efficient numerical discretizations, significantly reducing the number of degrees of freedom required for the solution of the corresponding integral equations. The results are illustrated by several numerical examples.

We present a first-order formulation of the Calogero model in external potentials in terms of a generating function, which simplifies the derivation of its dual form. Solitons naturally appear in this formulation as particles of negative mass. Using this method, we obtain the dual form of Calogero particles in external quartic, trigonometric and hyperbolic potentials, which were known to be integrable but had no known dual formulation. We derive the corresponding soliton solutions, generalizing earlier results for the harmonic Calogero system, and present numerical results that demonstrate the integrable nature of the soliton motion. We also give the collective fluid mechanical formulation of these models and derive the corresponding fluid soliton solutions in terms of meromorphic fields, commenting on issues of stability and integrability., Comment: 29 pages, 3 tables, 4 figures