Reaction‐Diffusion Processes and Their Connection with Integrable Quantum Spin Chains

This is a pedagogical account on reaction-diffusion systems and their relationship with integrable quantum spin chains. Reaction-diffusion systems are paradigmatic examples of non-equilibrium systems. Their long-time behaviour is strongly influenced through fluctuation effects in low dimensions which renders the habitual mean-field cinetic equations inapplicable. Starting from the master equation rewritten as a Schr\"odinger equation with imaginary time, the associated quantum hamiltonian of certain one-dimensional reaction-diffusion models is closely related to integrable magnetic chains. The relationship with the Hecke algebra and its quotients allows to identify integrable reaction-diffusion models and, through the Baxterization procedure, relate them to the solutions of Yang-Baxter equations which can be solved via the Bethe ansatz. Methods such as spectral and partial integrability, free fermions, similarity transformations or diffusion algebras are reviewed, with several concrete examples treated explicitly. An outlook on how the recently-introduced concept of local scale invariance might become useful in the description of non-equilibrium ageing phenomena is presented, with particular emphasis on the kinetic Ising model with Glauber dynamics.

We review the dynamics after quantum quenches in integrable quantum spin chains. We give a pedagogical introduction to relaxation in isolated quantum systems, and discuss the description of the steady state by (generalized) Gibbs ensembles. We then turn to general features in the time evolution of local observables after the quench, using a simple model of free fermions as an example. In the second part we present an overview of recent progress in describing quench dynamics in two key paradigms for quantum integrable models, the transverse field Ising chain and the anisotropic spin-1/2 Heisenberg chain.

We construct integrable Hamiltonian systems on <f>$G/K$</f>, where <f>$G$</f> is a coboundary Poisson–Lie group and <f>$K$</f> is a Lie subgroup arising as the fixed point set of a group automorphism <f>$\\sigma $</f> of <f>$G$</f> satisfying the classical reflection equation. We show that the time evolution of these systems is described by a Lax equation, and under a factorizability assumption, present its solution in terms of a factorization problem in <f>$G$</f>. Our construction is closely related to the semiclassical limit of Sklyanin's integrable quantum spin chains with reflecting boundaries.Communicated by Anton Alekseev

Motivated by a study of the crossing symmetry of the asymmetric twin or ‘gemini’representation of the affine Hecke algebra we give a construction for crossing tensor spacerepresentations of ordinary Hecke algebras. These representations build solutions to theYang–Baxter equation satisfying the crossing condition (that is, integrable quantum spinchains). We show that every crossing representation of the Temperley–Lieb algebra appearsin this construction, and in particular that this construction builds new representations. Weextend these to new representations of the blob algebra, which build new solutions to theboundary Yang–Baxter equation (i.e. open spin chains with integrable boundaryconditions).We prove that the open spin chain Hamiltonian derived from Sklyanin’s commutingtransfer matrix using such a solution can always be expressed as the representation of anelement of the blob algebra, and determine this element. We determine the representationtheory (irreducible content) of the new representations and hence show that all suchHamiltonians have the same spectrum up to multiplicity, for any given value of thealgebraic boundary parameter. (A corollary is that our models have the samespectrum as the open XXZ chain with nondiagonal boundary—despite differingfrom this model in having reference states.) Using these multiplicity data, andother ideas, we investigate the underlying quantum group symmetry of the newHamiltonians. We derive the form of the spectrum and the Bethe ansatz equations.

Integrable quantum spin chains display distinctive physical properties making them a laboratory to test and assess different states of matter. The study of the finite temperature properties is possible by use of the thermodynamic Bethe ansatz, however at the expense of dealing with non-linear integral equations for, in general, infinitely many auxiliary functions. The definition of an alternative finite set of auxiliary functions allowing for the complete description of their thermodynamic properties at finite temperature and fields has been elusive. Indeed, in the context of sln-invariant models satisfactory auxiliary functions have been established only for n≤4. In this paper we take a step further by proposing a systematic approach to generate finite sets of auxiliary functions for sln-invariant models. We refer to this construction as the canonical formulation. The numerical efficiency is illustrated for n=5, for which we present some of the thermodynamic properties of the corresponding spin chain.

A family of A2n(2) integrable open spin chains with Uq(Cn) symmetry was recently identified in arXiv:1702.01482. We identify here in a similar way a family of A2n−1(2) integrable open spin chains with Uq(Dn) symmetry, and two families of Dn+1(2) integrable open spin chains with Uq(Bn) symmetry. We discuss the consequences of these symmetries for the degeneracies and multiplicities of the spectrum. We propose Bethe ansatz solutions for two of these models, whose completeness we check numerically for small values of n and chain length N. We find formulas for the Dynkin labels in terms of the numbers of Bethe roots of each type, which are useful for determining the corresponding degeneracies. In an appendix, we briefly consider Dn+1(2) chains with other integrable boundary conditions, which do not have quantum group symmetry.

The exact solution of an integrable anisotropic Heisenberg spin chain with nearest-neighbor, next-nearest-neighbor and scalar chirality couplings is studied, where the boundary condition is the antiperiodic one. The detailed construction of Hamiltonian and the proof of integrability are given. The antiperiodic boundary condition breaks the U(1)-symmetry of the system and we use the off-diagonal Bethe Ansatz to solve it. The energy spectrum is characterized by the inhomogeneous T−Q relations and the contribution of the inhomogeneous term is studied. The ground state energy and the twisted boundary energy in different regions are obtained. We also find that the Bethe roots at the ground state form the string structure if the coupling constant J=−1 although the Bethe Ansatz equations are the inhomogeneous ones.

Evolution of reaction-diffusion patterns in infinite and bounded domains

The authors have generalized Sklyanin's approach of constructing open integrable quantum spin chains to the case of PT-invariant R matrices. They formulate a fusion procedure for such chains. In particular, they show that the fused transfer matrix can be expressed in terms of products of the original transfer matrix and products of certain quantum determinants which can be explicitly evaluated. Applications of these results include constructing open integrable higher-spin chains, as well as obtaining functional equations for transfer-matrix eigenvalues, which may be solved by an analytical Bethe ansatz.

In this contribution we briefly review recent developments in the theory of long-range integrable spin chains. These spin chains constitute a natural generalisation of the well-studied integrable nearest-neighbour chains and are of particular relevance to the integrability in the AdS/CFT correspondence since the dilatation operator in the asymptotic region is conjectured to be a Hamiltonian of an integrable long-range psu spin chain.

Asymptotic stability results of generalized discrete time reaction diffusion system applied to Lengyel-Epstein and Dagn Harrison models

Tooth patterning and evolution

We generalize our recent discrete cellular automata (CA) model of excitable media [Y. B. Chernyak, A. B. Feldman, and R. J. Cohen, Phys. Rev. E 55, 3215 (1997)] to incorporate the effects of inhibitory processes on the propagation of the excitation wave front. In the common two variable reaction-diffusion (RD) models of excitable media, the inhibitory process is described by the $v$ ``controller'' variable responsible for the restoration of the equilibrium state following excitation. In myocardial tissue, the inhibitory effects are mainly due to the inactivation of the fast sodium current. We represent inhibition using a physical model in which the ``source'' contribution of excited elements to the excitation of their neighbors decreases with time as a simple function with a single adjustable parameter (a rate constant). We sought specific solutions of the CA state transition equations and obtained (both analytically and numerically) the dependence of the wave-front speed $c$ on the four model parameters and the wave-front curvature $\ensuremath{\kappa}$. By requiring that the major characteristics of $c(\ensuremath{\kappa})$ in our CA model coincide with those obtained from solutions of a specific RD model, we find a unique set of CA parameter values for a given excitable medium. The basic structure of our CA solutions is remarkably similar to that found in typical RD systems (similar behavior is observed when the analogous model parameters are varied). Most notably, the ``turn-on'' of the inhibitory process is accompanied by the appearance of a solution branch of slow speed, unstable waves. Additionally, when $\ensuremath{\kappa}$ is small, we obtain a family of ``eikonal'' relations $c(\ensuremath{\kappa})$ that are suitable for the kinematic analysis of traveling waves in the CA medium. We compared the solutions of the CA equations to CA simulations for the case of plane waves and circular (target) waves and found excellent agreement. We then studied a spiral wave using the CA model adjusted to a specific RD system and found good correspondence between the shapes of the RD and CA spiral arms in the region away from the tip where kinematic theory applies. Our analysis suggests that only four physical parameters control the behavior of wave fronts in excitable media.

In the present paper, I propose a simple generic discretized model for the formation of various patterns based on a reaction–diffusion model of the growth patterns of bacterial colonies. Bacterial colonies show several patterns, with a characteristic pattern being a branching dendritic one with various branch densities. In addition to living systems, these branching patterns can also be observed in physical and chemical systems. Several reaction–diffusion models have been proposed to describe their formation. The model proposed here is a simplified discretized one based on previous models and involves three variables; namely, activators, inactivators, and the substrate. By changing only the mobility of the activators and the initial substrate, the proposed model can reproduce different patterns that are based on the branching patterns. I also investigate the characteristics of the shape and growth rate of the pattern formation.

The most general nonuniform reaction-diffusion models on a one-dimensional lattice with boundaries, for which the time evolution equations of corre- lation functions are closed, are considered. A transfer matrix method is used to find the static solution. It is seen that this transfer matrix can be obtained in a closed form, if the reaction rates satisfy certain conditions. We call such models superautonomous. Possible static phase transitions of such models are investigated. At the end, as an example of superau- tonomous models, a nonuniform voter model is introduced, and solved explicitly.