Toda Chain Research Articles

Particle scattering and fusion for the Ablowitz–Ladik chain

The Ablowitz–Ladik (AL) chain is an integrable discretized version of the nonlinear Schrödinger equation. We report on a novel underlying Hamiltonian particle system with properties similar to the ones known for the classical Toda chain and Calogero fluid with 1/sinh2 pair interaction. Boundary conditions are imposed such that, both in the distant past and future, particles have a constant velocity. We establish the many-particle scattering for the AL chain and obtain properties known for generic integrable many-body systems. For a specific choice of the chain, real initial data remain real in the course of time. Then, asymptotically, particles move in pairs with a velocity-dependent size and scattering shifts are governed by the fusion rule.

Non-abelian discrete Toda chains and related lattices

We have derived a non-abelian analog for the two-dimensional discrete Toda lattice which possesses solutions in terms of quasideterminants and admits Lax pairs of different forms. Its connection with non-abelian analogs for several well-known (1+1) and one-dimensional lattices is discussed. In particular, we consider a non-commutative analog of the scheme: discrete Toda equations → Somos-N sequences → discrete Painlevé equations.

Weak noise theory of the O'Connell-Yor polymer as an integrable discretization of the nonlinear Schrödinger equation.

We investigate and solve the weak noise theory for the semidiscrete O'Connell-Yor directed polymer. In the large deviation regime, the most probable evolution of the partition function obeys a classical nonlinear system which is a nonstandard discretization of the nonlinear Schrödinger equationwith mixed initial-final conditions. We show that this system is integrable and find its general solution through an inverse scattering method and a non-standard Fredholm determinant framework that we develop. This allows us to obtain the large deviation rate function of the free energy of the polymer model from its conserved quantities and to study its convergence to the large deviations of the Kardar-Parisi-Zhang equation. Our model also degenerates to the classical Toda chain, which further substantiates the applicability of our Fredholm framework.

Dynamical chaos in the integrable Toda chain induced by time discretization.

We use the Toda chain model to demonstrate that numerical simulation of integrable Hamiltonian dynamics using time discretization destroys integrability and induces dynamical chaos. Specifically, we integrate this model with various symplectic integrators parametrized by the time step τ and measure the Lyapunov time TΛ (inverse of the largest Lyapunov exponent Λ). A key observation is that TΛ is finite whenever τ is finite but diverges when τ→0. We compare the Toda chain results with the nonintegrable Fermi-Pasta-Ulam-Tsingou chain dynamics. In addition, we observe a breakdown of the simulations at times TB≫TΛ due to certain positions and momenta becoming extremely large ("Not a Number"). This phenomenon originates from the periodic driving introduced by symplectic integrators and we also identify the concrete mechanism of the breakdown in the case of the Toda chain.

Lagrangian multiforms on coadjoint orbits for finite-dimensional integrable systems

Lagrangian multiforms provide a variational framework to describe integrable hierarchies. The case of Lagrangian 1-forms covers finite-dimensional integrable systems. We use the theory of Lie dialgebras introduced by Semenov-Tian-Shansky to construct a Lagrangian 1-form. Given a Lie dialgebra associated with a Lie algebra g\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathfrak {g}$$\\end{document} and a collection Hk\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H_k$$\\end{document}, k=1,⋯,N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$k=1,\\dots ,N$$\\end{document}, of invariant functions on g∗\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathfrak {g}^*$$\\end{document}, we give a formula for a Lagrangian multiform describing the commuting flows for Hk\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H_k$$\\end{document} on a coadjoint orbit in g∗\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathfrak {g}^*$$\\end{document}. We show that the Euler–Lagrange equations for our multiform produce the set of compatible equations in Lax form associated with the underlying r-matrix of the Lie dialgebra. We establish a structural result which relates the closure relation for our multiform to the Poisson involutivity of the Hamiltonians Hk\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H_k$$\\end{document} and the so-called double zero on the Euler–Lagrange equations. The construction is extended to a general coadjoint orbit by using reduction from the free motion of the cotangent bundle of a Lie group. We illustrate the dialgebra construction of a Lagrangian multiform with the open Toda chain and the rational Gaudin model. The open Toda chain is built using two different Lie dialgebra structures on sl(N+1)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathfrak {sl}(N+1)$$\\end{document}. The first one possesses a non-skew-symmetric r-matrix and falls within the Adler–Kostant–Symes scheme. The second one possesses a skew-symmetric r-matrix. In both cases, the connection with the well-known descriptions of the chain in Flaschka and canonical coordinates is provided.

Dubrovin method and the Toda chain

A hierarchy of Lax pairs with 2 × 2 2\times 2 matrix coefficients is presented. The compatibility conditions for these pairs include the Toda chain equation, and other differential-difference integrable systems. Various kinds of finite gap solutions for such systems are constructed. Examples of simplest one- and two-phase solutions are given, together with the corresponding spectral curves.

Asymptotic expansion of matrix models in the multi-cut regime

Abstract We establish the asymptotic expansion in $\beta $ matrix models with a confining, off-critical potential in the regime where the support of the equilibrium measure is a finite union of segments. We first address the case where the filling fractions of these segments are fixed and show the existence of a $\frac {1}{N}$ expansion. We then study the asymptotics of the sum over the filling fractions to obtain the full asymptotic expansion for the initial problem in the multi-cut regime. In particular, we identify the fluctuations of the linear statistics and show that they are approximated in law by the sum of a Gaussian random variable and an independent Gaussian discrete random variable with oscillating center. Fluctuations of filling fractions are also described by an oscillating discrete Gaussian random variable. We apply our results to study the all-order small dispersion asymptotics of solutions of the Toda chain associated with the one Hermitian matrix model ( $\beta = 2$ ) as well as orthogonal ( $\beta = 1$ ) and skew-orthogonal ( $\beta = 4$ ) polynomials outside the bulk.

Spectral analysis of an open $q$-difference Toda chain with two-sided boundary interactions on the finite integer lattice

A quantum n -particle model consisting of an open q -difference Toda chain with two-sided boundary interactions is placed on a finite integer lattice. The spectrum and eigenbasis are computed by establishing the equivalence with a previously studied q -boson model from which the quantum integrability is inherited. Specifically, the q -boson-Toda correspondence in question yields Bethe Ansatz eigenfunctions in terms of hyperoctahedral Hall–Littlewood polynomials and provides the pertinent solutions of the Bethe Ansatz equations via the global minima of corresponding Yang–Yang-type Morse functions.

Integrable hydrodynamics of Toda chain: case of small systems

AbstractPassing from a microscopic discrete lattice system with many degrees of freedom to a mesoscopic continuum system described by a few coarse-grained equations is challenging. The common folklore is to take the thermodynamic limit so that the physics of the discrete lattice describes the continuum results. The analytical procedure to do so relies on defining a small length scale (typically the lattice spacing) to coarse grain the microscopic evolution equations. Moving from the microscopic scale to the mesoscopic scale then requires careful approximations. In this work, we numerically test the coarsening in a Toda chain, which is an interacting integrable system, i.e., a system with a macroscopic number of conserved charges. Specifically, we study the spreading of fluctuations by computing the spatio-temporal thermal correlations with three different methods: (a) using microscopic molecular dynamics simulation with a large number of particles; (b) solving the generalized hydrodynamics equation; (c) solving the linear Euler scale equations for each conserved quantities. Surprisingly, the results for the small systems (c) match the thermodynamic results in (a) and (b) for macroscopic systems. This reiterates the importance and validity of integrable hydrodynamics in describing experiments in the laboratory, where we typically have microscopic systems.

Hermite–Padé approximation and integrability

We show that solution to the Hermite–Padé type I approximation problem leads in a natural way to a subclass of solutions of the Hirota (discrete Kadomtsev–Petviashvili) system and of its adjoint linear problem. Our result explains the appearance of various ingredients of the integrable systems theory in application to multiple orthogonal polynomials, numerical algorithms, random matrices, and in other branches of mathematical physics and applied mathematics where the Hermite–Padé approximation problem is relevant. We present also the geometric algorithm, based on the notion of Desargues maps, of construction of solutions of the problem in the projective space over the field of rational functions. As a byproduct we obtain the corresponding generalization of the Wynn recurrence. We isolate the boundary data of the Hirota system which provide solutions to Hermite–Padé problem showing that the corresponding reduction lowers dimensionality of the system. In particular, we obtain certain equations which, in addition to the known ones given by Paszkowski, can be considered as direct analogs of the Frobenius identities. We study the place of the reduced system within the integrability theory, which results in finding multidimensional (in the sense of number of variables) extension of the discrete-time Toda chain equations.

Integration of the finite complex Toda lattice with a self-consistent source

In the paper, we derive a finite complex Toda lattice with a self-consistent source. We discuss the complete integrability of the constructed systems that is based on the transformation to the spectral data of an associated finite Jacobi matrix. We show that the finite complex Toda lattice with a self-consistent source is also an important theoretical model as it is a completely integrable system. Namely, we have determined the time evolution of the spectral data for the Jacobi matrix associated with the solution of a finite complex Toda chain with a self-consistent source. Then, using the solution of the inverse spectral problem with respect to the time-dependent spectral data, we recover the time-dependent Jacobi matrix and hence the desired solution of the finite complex Toda chain with a self-consistent source. For the case N=2, explicit formulas for the solution of the problem under consideration are obtained.

Effective Stochastic Model for Chaos in the Fermi–Pasta–Ulam–Tsingou Chain

Understanding the interplay between different wave excitations, such as phonons and localized solitons, is crucial for developing coarse-grained descriptions of many-body, near-integrable systems. We treat the Fermi–Pasta–Ulam–Tsingou (FPUT) non-linear chain and show numerically that at short timescales, relevant to the largest Lyapunov exponent, it can be modeled as a random perturbation of its integrable approximation—the Toda chain. At low energies, the separation between two trajectories that start at close proximity is dictated by the interaction between few soliton modes and an intrinsic, apparent bath representing a background of many radiative modes. It is sufficient to consider only one randomly perturbed Toda soliton-like mode to explain the power-law profiles reported in previous works, describing how the Lyapunov exponent of large FPUT chains decreases with the energy density of the system.

Chaotic string motion in a near pp-wave limit

We revisit classical string motion in a near pp-wave limit of AdS5 × S5. It is known that the Toda lattice models are integrable. But if the exponential potential is truncated at finite order, then the system may become non-integrable. In particular, when the exponential potential in a three-particle periodic Toda chain is truncated at the third order of the dynamical variables, the resulting system becomes a well-known non-integrable system, Henon-Heiles model. The same thing may happen in a near pp-wave limit of AdS5 × S5, on which the classical string motion becomes chaotic.

Toda equations for surface defects in SYM and instanton counting for classical Lie groups

The partition function of super Yang-Mills theories with arbitrary simple gauge group coupled to a self-dual Ω background is shown to be fully determined by studying the renormalization group equations relevant to the surface operators generating its one-form symmetries. The corresponding system of equations results in a non-autonomous Toda chain on the root system of the Langlands dual, the evolution parameter being the RG scale. A systematic algorithm computing the full multi-instanton corrections is derived in terms of recursion relations whose gauge theoretical solution is obtained just by fixing the perturbative part of the IR prepotential as its asymptotic boundary condition for the RGE. We analyze the explicit solutions of the τ-system for all the classical groups at the diverse levels, extend our analysis to affine twisted Lie algebras and provide conjectural bilinear relations for the τ-functions of linear quiver gauge theory.

Regular Representation of the Group GL(N, ℝ): Factorization, Casimir Operators and Toda Chain

Regular Representation of the Group GL(N, ℝ): Factorization, Casimir Operators and Toda Chain

Mirror symmetry for extended affine Weyl groups

We give a uniform, Lie-theoretic mirror symmetry construction for the Frobenius manifolds defined by Dubrovin--Zhang in arXiv:hep-th/9611200 on the orbit spaces of extended affine Weyl groups, including exceptional Dynkin types. The B-model mirror is given by a one-dimensional Landau--Ginzburg superpotential constructed from a suitable degeneration of the family of spectral curves of the affine relativistic Toda chain for the corresponding affine Poisson--Lie group. As applications of our mirror theorem we give closed-form expressions for the flat coordinates of the Saito metric and the Frobenius prepotentials in all Dynkin types, compute the topological degree of the Lyashko--Looijenga mapping for certain higher genus Hurwitz space strata, and construct hydrodynamic bihamiltonian hierarchies (in both Lax--Sato and Hamiltonian form) that are root-theoretic generalisations of the long-wave limit of the extended Toda hierarchy.

From quantum groups to Liouville and dilaton quantum gravity

We investigate the underlying quantum group symmetry of 2d Liouville and dilaton gravity models, both consolidating known results and extending them to the cases with mathcal{N} = 1 supersymmetry. We first calculate the mixed parabolic representation matrix element (or Whittaker function) of Uq( mathfrak{sl} (2, ℝ)) and review its applications to Liouville gravity. We then derive the corresponding matrix element for Uq( mathfrak{osp} (1|2, ℝ)) and apply it to explain structural features of mathcal{N} = 1 Liouville supergravity. We show that this matrix element has the following properties: (1) its q → 1 limit is the classical OSp+(1|2, ℝ) Whittaker function, (2) it yields the Plancherel measure as the density of black hole states in mathcal{N} = 1 Liouville supergravity, and (3) it leads to 3j-symbols that match with the coupling of boundary vertex operators to the gravitational states as appropriate for mathcal{N} = 1 Liouville supergravity. This object should likewise be of interest in the context of integrability of supersymmetric relativistic Toda chains. We furthermore relate Liouville (super)gravity to dilaton (super)gravity with a hyperbolic sine (pre)potential. We do so by showing that the quantization of the target space Poisson structure in the (graded) Poisson sigma model description leads directly to the quantum group Uq( mathfrak{sl} (2, ℝ)) or the quantum supergroup Uq( mathfrak{osp} (1|2, ℝ)).

The anti-Fermi–Pasta–Ulam–Tsingou problem in one-dimensional diatomic lattices

We study the thermalization dynamics of one-dimensional diatomic lattices (which represents the simplest system possessing multi-branch phonons), exemplified by the famous Fermi–Pasta–Ulam–Tsingou (FPUT)-β and the Toda models. Here we focus on how the system relaxes to the equilibrium state when part of highest-frequency optical modes are initially excited, which is called the anti-FPUT problem comparing with the original FPUT problem (low frequency excitations of the monatomic lattice). It is shown numerically that the final thermalization time T eq of the diatomic FPUT-β chain depends on whether its acoustic modes are thermalized, whereas the T eq of the diatomic Toda chain depends on the optical ones; in addition, the metastable state of both models have different energy distributions and lifetimes. Despite these differences, in the near-integrable region, the T eq of both models still follows the same scaling law, i.e. T eq is inversely proportional to the square of the perturbation strength. Finally, comparisons of the thermalization behavior between different models under various initial conditions are briefly summarized.

Large deviations for Gibbs ensembles of the classical Toda chain

We prove large deviation principles for the distribution of the empirical measure of the eigenvalues of Lax matrices following the Generalized Gibbs ensembles of the classical Toda chain introduced in [9]. We deduce the almost sure convergence of this empirical measure towards a limit which we describe in terms of the limiting empirical measure of Beta-ensembles. Our results apply to general smooth potentials.

On the quantum $$\mathfrak{osp}(1|2\ell)$$ Toda chain

The orthosymplectic super Lie algebra $\mathfrak{osp}(1|\,2\ell)$ is the closest analog of standard Lie algebras in the world of super Lie algebras. We demonstrate that the corresponding $\mathfrak{osp}(1|\,2\ell)$-Toda chain turns out to be an instance of a $BC_\ell$-Toda chain. The underlying reason for this relation is discussed.