Integrable Research Articles

The bifurcation, chaotic behavior and exact solutions of the fractional stochastic Jimbo–Miwa equations

The Jimbo–Miwa equation is the second equation in the KP hierarchy of integrable systems. In this paper, this equation is extended and introduced with the stochastic process and fractional derivatives. Firstly, the phase portrait of the Hamiltonian system generated by it is studied to understand its bifurcation behavior. Additionally, non-periodic and periodic perturbation terms are added to this system. Different values are assigned to the parameters in the perturbation terms to analyze its sensitivity and the resulting chaos is obtained. Finally, through integration techniques, the expression of the solution of this equation is obtained. These solutions are related to rational functions, trigonometric functions, exponential functions and Jacobi elliptic functions. To observe the form of the solutions more intuitively, 3D and 2D numerical simulations are conducted on the solutions and the solution images of the stochastic fractional differential equation are given by Matlab software. Compared with the existing literature, the research on the stochastic fractional equation of this equation is relatively rare and the analysis of the phase portrait is even scarcer. Our solution method is quite different from that in the previous literature. Therefore, this paper is novel. The conclusion of this paper will be of great help for the practical application of this equation.

Exact Solutions of Boussinesq Equations By Hirota Direct Method

Boussinesq Equations (BSQ) are the focus of this article. First, we provide a basic overview of Hirota's D operator, which is used to build multi-soliton solutions for equations involving nonlinear evolution. After that, some details regarding fourth-order BSQ are provided, and we use Hirota's direct method to find a one-solution solution. Hirota's bilinear approach is also used to solve a class of sixth-order Boussinesq-like equations with nonlinear evolution. The outcomes verified that this approach requires complete integrability.

More on bilinearization of supersymmetric coupled dispersionless (SUSY-CD) integrable system

We study the bilinearization of [Formula: see text] Supersymmetric Coupled Dispersionless (SUSY-CD) integrable system by introducing transformations defined in terms of bosonic and fermionic superfields. From the bosonic and fermionic superfield bilinear equations, we obtain the superfield soliton solutions of SUSY-CD integrable system.

Commutative families in DIM algebra, integrable many-body systems and q, t matrix models

We extend our consideration of commutative subalgebras (rays) in different representations of the W1+∞ algebra to the elliptic Hall algebra (or, equivalently, to the Ding-Iohara-Miki (DIM) algebra Uq,tgl̂̂1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {U}_{q,t}\\left({\\hat{\\hat{\\mathfrak{gl}}}}_1\\right) $$\\end{document}). Its advantage is that it possesses the Miki automorphism, which makes all commutative rays equivalent. Integrable systems associated with these rays become finite-difference and, apart from the trigonometric Ruijsenaars system not too much familiar. We concentrate on the simplest many-body and Fock representations, and derive explicit formulas for all generators of the elliptic Hall algebra en,m. In the one-body representation, they differ just by normalization from znqmD̂\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {z}^n{q}^{m\\hat{D}} $$\\end{document} of the W1+∞ Lie algebra, and, in the N -body case, they are non-trivially generalized to monomials of the Cherednik operators with action restricted to symmetric polynomials. In the Fock representation, the resulting operators are expressed through auxiliary polynomials of n variables, which define weights in the residues formulas. We also discuss q, t-deformation of matrix models associated with constructed commutative subalgebras.

Decorating the gauge/YBE correspondence

In this paper, we aim to study the three-dimensional N=2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathcal {N}}=2$$\\end{document} supersymmetric dual gauge theories on Sb3/Zr\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$S_b^3/\\mathbb {Z}_r$$\\end{document} in the context of the gauge/YBE correspondence. We consider hyperbolic integral identities acquired via the equality of supersymmetric lens partition functions as solutions to the decoration transformation and the flipping relation in statistical mechanics. The solutions of those transformations aim at investigating various decorated lattice models possessing the Boltzmann weights of integrable Ising-like models obtained via the gauge/YBE correspondence. We also constructed The Bailey pairs for the decoration transformation and the flipping relation.

Nonlinear evolution of disturbances in higher time-derivative theories

We investigate the evolution of localized initial value profiles when propagated in integrable versions of higher time-derivative theories. In contrast to the standard cases in nonlinear integrable systems, where these profiles evolve into a specific number of N-soliton solutions as dictated by the conservation laws, in the higher time-derivative theories the theoretical prediction is that the initial profiles can settle into either two-soliton solutions or into any number of N-soliton solutions. In the latter case this implies that the solutions exhibit oscillations that spread in time but remain finite. We confirm these analytical predictions by explicitly solving the associated Cauchy problem numerically with multiple initial profiles for various higher time-derivative versions of integrable modified Korteweg-de Vries equations. In the case with the theoretical possibility of a decay into two-soliton solutions, the emergence of underlying singularities may prevent the profiles from fully developing or may be accompanied by oscillatory, chargeless standing waves at the origin.

Bound-state confinement after trap-expansion dynamics in integrable systems

Abstract Integrable systems possess stable families of quasiparticles, which are composite objects (bound states) of elementary excitations. Motivated by recent quantum computer experiments, we investigate bound-state transport in the spin- 1 / 2 anisotropic Heisenberg chain (XXZ chain). Specifically, we consider the sudden vacuum expansion of a finite region A prepared in a non-equilibrium state. In the hydrodynamic regime, if interactions are strong enough, bound states remain confined in the initial region. Bound-state confinement persists until the density of unbound excitations remains finite in the bulk of A. Since region A is finite, at asymptotically long times bound states are ‘liberated’ after the ‘evaporation’ of the unbound excitations. Fingerprints of confinement are visible in the space-time profiles of local spin-projection operators. To be specific, here we focus on the expansion of the p-Néel states, which are obtained by repetition of a unit cell with p up spins followed by p down spins. Upon increasing p, the bound-state content is enhanced. In the limit p → ∞ one obtains the domain-wall initial state. We show that for p < 4, only bound states with n > p are confined at large chain anisotropy. For p ≳ 4 , bound states with n = p are also confined, consistent with the absence of transport in the limit p → ∞ . The scenario of bound-state confinement leads to a hierarchy of timescales at which bound states of different sizes are liberated, which is also reflected in the dynamics of the von Neumann entropy.

[formula omitted]-periodic wave solutions of the [formula omitted] supersymmetric KdV equation

In this paper, the N-periodic wave solutions of the N=2 supersymmetric KdV equation are studied by combining the super Hirota bilinear form with the super Riemann-theta function, which can be used to describe new phenomena on super quasi-periodic waves with the fermionic field. With the aid of the Gauss–Newton method, the three-periodic and four-periodic wave solutions are obtained. In particular, these quasi-periodic waves can produce parallel, crossed and degenerated patterns. The analytical method related to the characteristic lines is used to analyze the dynamic characteristics of the three-periodic and four-periodic waves. In addition, it has been indicated that N-periodic waves can exist in the supersymmetric integrable systems.

Polynomial dynamical systems associated with the KdV hierarchy

In 1974, S.P. Novikov introduced the stationary n-equations of the Korteweg–de Vries hierarchy, namely the n-Novikov equations. These are associated with integrable polynomial dynamical systems, with polynomial 2n integrals, in ℂ3n. In this paper, we construct an infinite-dimensional polynomial dynamical system that is universal for all dynamical systems corresponding to the n-Novikov equations. Thus, we solve the well-known problem of the relationship between the n-Novikov equations for different n.

Flat limit of massless scalar scattering in AdS2

We explore the flat limit of massless scalar scattering in AdS2. We derive the 1→1S-matrix from the CFT 2-point function. We show a key property of the 2→2S-matrix in 2d, where the contact interaction in the flat limit gives momentum conserving delta function. We show the factorization of the n→nS-matrix for integrable models in the flat limit, focusing on contact interactions. We calculate the S-matrix by linking the CFT operator on the AdS boundary to the scattering state in flat-space. We use bulk operator reconstruction to study massless scalar scattering in the flat limit and solve the Klein-Gordon equation in global AdS2 for the massless scalar field. The solution is simple, involving a pure phase in global time and a sinusoidal function in the radial coordinate. This simplicity also extends to the smearing function, allowing us to map the scattering state to the CFT operator while taking AdS corrections into account.

Development and Analysis of novel Integrable Nonlinear Dynamical Systems on Quasi-One-Dimensional Lattices. Parametrically Driven Nonlinear System of Pseudo-Excitations on a Two-Leg Ladder Lattice

Following the main principles of developing the evolutionary nonlinear integrable systems on quasi-one-dimensional lattices, we suggest a novel nonlinear integrable system of parametrically driven pseudo-excitations on a regular two-leg ladder lattice. The initial (prototype) form of the system is derivable in the framework of semi-discrete zero-curvature equation with the spectral and evolution operators specified by the properly organized 3 × 3 square matrices. Although the lowest conserved local densities found via the direct recursive method do not prompt us the algebraic structure of system’s Hamiltonian function, but the heuristically substantiated search for the suitable two-stage transformation of prototype field functions to the physically motivated ones has allowed to disclose the physically meaningful nonlinear integrable system with time-dependent longitudinal and transverse inter-site coupling parameters. The time dependencies of inter-site coupling parameters in the transformed system are consistently defined in terms of the accompanying parametric driver formalized by the set of four homogeneous ordinary linear differential equations with the time-dependent coefficients. The physically meaningful parametrically driven nonlinear system permits its concise Hamiltonian formulation with the two pairs of field functions serving as the two pairs of canonically conjugated field amplitudes. The explicit example of oscillatory parametric drive is described in full mathematical details.

Statistics of Matrix Elements of Local Operators in Integrable Models

We study the statistics of matrix elements of local operators in the basis of energy eigenstates in a paradigmatic, integrable, many-particle quantum theory, the Lieb-Liniger model of bosons with repulsive delta-function interactions. Using methods of quantum integrability, we determine the scaling of matrix elements with system size. As a consequence of the extensive number of conservation laws, the structure of matrix elements is fundamentally different from, and much more intricate than, the predictions of the eigenstate thermalization hypothesis for generic models. We uncover an interesting connection between this structure for local operators in interacting integrable models and the one for local operators that are not local with respect to the elementary excitations in free theories. We find that typical off-diagonal matrix elements ⟨μ|O|λ⟩ in the same macrostate scale as exp(−cOLln(L)−LMμ,λO), where the probability distribution function for Mμ,λO is well described by Fréchet distributions and cO depends only on macrostate information. In contrast, typical off-diagonal matrix elements between two different macrostates scale as exp(−dOL2), where dO depends only on macrostate information. Diagonal matrix elements depend only on macrostate information up to finite-size corrections. Published by the American Physical Society 2024

Quantum next generation reservoir computing: an efficient quantum algorithm for forecasting quantum dynamics

Next Generation Reservoir Computing (NG-RC) is a modern class of model-free machine learning that enables an accurate forecasting of time series data generated by dynamical systems. We demonstrate that NG-RC can accurately predict full many-body quantum dynamics in both integrable and chaotic systems. This is in contrast to the conventional application of reservoir computing that concentrates on the prediction of the dynamics of observables. In addition, we apply a technique which we refer to as skipping ahead to predict far future states accurately without the need to extract information about the intermediate states. However, adopting a classical NG-RC for many-body quantum dynamics prediction is computationally prohibitive due to the large Hilbert space of sample input data. In this work, we propose an end-to-end quantum algorithm for many-body quantum dynamics forecasting with a quantum computational speedup via the block-encoding technique. This proposal presents an efficient model-free quantum scheme to forecast quantum dynamics coherently, bypassing inductive biases incurred in a model-based approach.

Integrable sigma model with generalized [formula omitted] structure, Yang-Baxter sigma model with generalized complex structure and multi-Yang-Baxter sigma model

We construct an integrable sigma model with a generalized F structure, which involves a generalized Nijenhuis structure J satisfying J3=−J. Utilizing the expression of the generalized complex structure on the metric Lie group manifold G in terms of operator relations on its Lie algebra g, we formulate a Yang-Baxter sigma model with a generalized complex structure. Additionally, we present multi-Yang-Baxter sigma models featuring two and three compatible Nijenhuis structures. Examples for each of these models are provided.

Periodic orbits in Fermi-Pasta-Ulam-Tsingou systems.

The Fermi-Pasta-Ulam-Tsingou (FPUT) paradox is the phenomenon whereby a one-dimensional chain of oscillators with nonlinear couplings shows long-lived nonergodic behavior prior to thermalization. The trajectory of the system in phase space, with a long-wavelength initial condition, closely follows that of the Toda model over short times, as both systems seem to relax quickly to a non-thermal, metastable state. Over longer times, resonances in the FPUT spectrum drive the system toward equilibrium, away from the Toda trajectory. Similar resonances are observed in q-breather spectra, suggesting that q-breathers are involved in the route toward thermalization. In this article, we first review previous important results related to the metastable state, solitons, and q-breathers. We then investigate orbit bifurcations of q-breathers and show that they occur due to resonances, where the q-breather frequencies become commensurate as mΩ1=Ωk. The resonances appear as peaks in the breather energy spectrum. Furthermore, they give rise to new "composite periodic orbits," which are nonlinear combinations of multiple q-breathers that exist following orbit bifurcations. We find that such resonances are absent in integrable systems, as a consequence of the (extensive number of) conservation laws associated with integrability.

Chaotic roots of the modular multiplication dynamical system in Shor's algorithm

Shor's factoring algorithm, believed to provide an exponential speedup over classical computation, relies on finding the period of an exactly periodic quantum modular multiplication operator. This exact periodicity is the hallmark of an integrable system, which is paradoxical from the viewpoint of quantum chaos, given that the classical limit of the modular multiplication operator is a highly chaotic system that occupies the “maximally random” Bernoulli level of the classical ergodic hierarchy. In this work, we approach this apparent paradox from a quantum dynamical systems viewpoint, and consider whether signatures of ergodicity and chaos may indeed be encoded in such an “integrable” quantization of a chaotic system. We show that Shor's modular multiplication operator, in specific cases, can be written as a superposition of quantized A-baker's maps exhibiting more typical signatures of quantum chaos and ergodicity. This work suggests that the integrability of Shor's modular multiplication operator may stem from the interference of other “chaotic” quantizations of the same family of maps, and paves the way for deeper studies on the interplay of integrability, ergodicity, and chaos in and via quantum algorithms. Published by the American Physical Society 2024

FateNet: an integration of dynamical systems and deep learning for cell fate prediction.

Understanding cellular decision-making, particularly its timing and impact on the biological system such as tissue health and function, is a fundamental challenge in biology and medicine. Existing methods for inferring fate decisions and cellular state dynamics from single-cell RNA sequencing data lack precision regarding decision points and broader tissue implications. Addressing this gap, we present FateNet, a computational approach integrating dynamical systems theory and deep learning to probe the cell decision-making process using scRNA-seq data. By leveraging information about normal forms and scaling behavior near bifurcations common to many dynamical systems, FateNet predicts cell decision occurrence with higher accuracy than conventional methods and offers qualitative insights into the new state of the biological system. Also, through in-silico perturbation experiments, FateNet identifies key genes and pathways governing the differentiation process in hematopoiesis. Validated using different scRNA-seq data, FateNet emerges as a user-friendly and valuable tool for predicting critical points in biological processes, providing insights into complex trajectories. github.com/ThomasMBury/fatenet.

Tautological relations and integrable systems

Tautological relations and integrable systems

One-dimensional Fermi polaron after a kick: Two-sided singularity of the momentum distribution, Bragg reflection and other exact results

A mobile impurity particle immersed in a quantum fluid forms a polaron - a quasiparticle consisting of the impurity and a local disturbance of the fluid around it. We ask what happens to a one-dimensional polaron after a kick, i.e. an abrupt application of a force that instantly delivers a finite impulse to the impurity. In the framework of an integrable model describing an impurity in a one-dimensional gas of fermions or hard-core bosons, we calculate the distribution of the polaron momentum established when the post-kick relaxation is over. A remarkable feature of this distribution is a two-sided power-law singularity. It emerges due to one of two processes. In the first process, the whole impulse is transferred to the polaron, without creating phonon-like excitations of the fluid. In the second process, the impulse is shared between the polaron and the center-of-mass motion of the fluid, again without creating any fluid excitations. The latter process is, in fact, a Bragg reflection at the edge of the emergent Brillouin zone. We carefully analyze the conditions for each of the two processes. The asymptotic form of the distribution in the vicinity of the singularity is derived.

Multi-pole soliton of discrete integrable equations and modified Riemann-Hilbert approach: discrete Hirota equation

In this study, the Riemann-Hilbert approach was developed and applied to the discrete Hirota equation. We constructed a modified Riemann-Hilbert problem compatible with the discrete Hirota equation and derived a reconstruction formula for its solutions. Because the characteristic function contains a potential, we modify the Riemann-Hilbert approach to make the Riemann-Hilbert matrix have good asymptotic properties. We believe that the modified Riemann-Hilbert approach can also be applied to other discrete integrable models. By using the direct method of Laurent series, we obtained the expression of multi-pole solutions for the discrete Hirota equation and demonstrated the dynamic behavior of some solutions.