An Exactly Solvable Asymmetric Simple Inclusion Process

Basic Information

• Paper ID: 2510.09191
• Title: An exactly solvable asymmetric simple inclusion process
• Authors: Arvind Ayyer (Indian Institute of Science), Samarth Misra (Indian Institute of Science)
• Classification: cond-mat.stat-mech math-ph math.CO math.MP math.PR
• Publication Date: October 13, 2025
• Paper Link: https://arxiv.org/abs/2510.09191

Abstract

This paper investigates a generalization of the asymmetric simple inclusion process (ASIP) on a periodic one-dimensional lattice, where integers in particle transition rates are deformed into their t-analogues. The authors term this process the (q, t, θ) ASIP, where q is the asymmetric hopping parameter and θ is the diffusion parameter. The study demonstrates that this process is a misanthrope process, and consequently the stationary distribution is independent of q. The paper computes the stationary distribution, single-point correlation functions, and stationary currents. Notably, it proves that at t=1, the single-point occupation probability follows a beta-binomial distribution. The authors also compute two-dimensional phase diagrams for different parameter intervals (t, θ) and verify them through simulations, discovering that the corrected form of stationary weights at t≠1 satisfies interesting palindromic and anti-palindromic symmetries.

Research Background and Motivation

Problem Background

1. Limitations of ASEP: The asymmetric simple exclusion process (ASEP), as a classical model of one-dimensional interacting particle systems, allows at most one particle per lattice site, exhibiting "fermionic" characteristics with mutual particle exclusion.
2. Introduction of Inclusion Processes: To study the "bosonic" counterpart, researchers introduced the symmetric inclusion process (SIP) and asymmetric inclusion process (ASIP), which allow multiple particles per lattice site, promoting particle clustering.
3. Limitations of Existing Models: Existing ASIP research primarily focuses on condensation phenomena under specific boundary conditions, lacking a more generalized theoretical framework.

Research Motivation

This paper aims to generalize the ASIP model by introducing t-analogues, constructing a more general (q, t, θ) ASIP framework to:

• Unify different existing ASIP variants
• Provide an exactly solvable mathematical framework
• Reveal new mathematical structures and physical phenomena

Core Contributions

1. Proposed the (q, t, θ) ASIP model: Generalized the classical ASIP model by replacing integers in transition rates with t-analogues
2. Proved misanthrope process property: Rigorously demonstrated that this model belongs to the misanthrope process class, with stationary distribution independent of asymmetric parameter q
3. Obtained exact stationary solutions: Computed explicit expressions for the stationary distribution, particularly closed-form solutions at t=1
4. Discovered beta-binomial distribution: Proved that single-point occupation probability follows a beta-binomial distribution at t=1
5. Constructed complete phase diagram: Analyzed the two-dimensional phase diagram in parameter space (t, θ), identifying different physical phases
6. Revealed palindromic symmetries: Discovered that stationary weights possess palindromic and anti-palindromic symmetries under specific parameterization
7. Constructed enriched processes: Defined enriched processes projecting onto the original process when t=1 and θ is an integer

Detailed Methodology

Task Definition

Study particle dynamics on a periodic one-dimensional lattice where:

• Input: L lattice sites, n indistinguishable particles, parameters (q, t, θ)
• Output: Stationary distribution, correlation functions, currents, and other observables
• Constraints: Particle number conservation, periodic boundary conditions

Model Architecture

Configuration Space

The configuration space is defined as: $\Omega_{L,n} = \left\{\eta = (\eta_1, \eta_2, \ldots, \eta_L) \in \{0, 1, \ldots, n\}^L \mid \sum_{i=1}^L \eta_i = n\right\}$

where $\eta_i$ denotes the number of particles at site i, with total configurations numbering $\binom{n+L-1}{n}$.

Transition Rates

For adjacent sites $(i, i+1)$ with occupation numbers $(\alpha, \beta)$:

• Forward transition: $(\alpha, \beta) \to (\alpha-1, \beta+1)$, rate $[\alpha]_t(\theta + [\beta]_t)$
• Backward transition: $(\alpha, \beta) \to (\alpha+1, \beta-1)$, rate $q[\beta]_t(\theta + [\alpha]_t)$

where $[k]_t = \frac{1-t^k}{1-t}$ is the t-analogue.

Misanthrope Process Verification

By verifying the conditions: $\frac{u(\beta, \alpha)}{u(\alpha+1, \beta-1)} = \frac{u(1, \alpha)u(\beta, 0)}{u(\alpha+1, 0)u(1, \beta-1)}$$u(\beta, \alpha) - u(\alpha, \beta) = u(\beta, 0) - u(\alpha, 0)$

the process is proven to belong to the misanthrope process class.

Technical Innovations

1. t-analogue generalization: Extended classical integer rates to t-analogues while maintaining exact solvability
2. Unified framework: Recovered known ASIP variants in different parameter limits
3. Product-form stationary distribution: Utilized misanthrope process theory to obtain product-form stationary distributions
4. Special parameterization: Introduced $\theta = \frac{1-a}{a(1-t)}$ parameterization to reveal symmetries

Experimental Setup

Numerical Simulation Parameters

• System size: Primarily L=20 sites, n=20 particles
• Parameter range: t and θ varied over multiple orders of magnitude
• Simulation method: Gillespie algorithm for stochastic simulation
• Statistics: 20 independent runs per parameter combination, collecting 10,000 stationary configurations per run

Evaluation Metrics

1. Stationary probability distribution: Comparison with theoretical predictions
2. Single-point occupation variance: Comparison of theoretical and simulation results
3. Phase diagram analysis: Dominant configuration types in different parameter regions
4. Convergence: Time scales for reaching stationarity

Experimental Results

Main Results

Stationary Distribution Verification

For small system L=3, n=4, theoretical predictions perfectly match numerical results:

• Configuration (0,0,4): theoretical probability 0.33, simulation results match
• Probability distributions of all other configurations agree with theory

Variance Verification

Theoretical formula for single-point occupation variance: $\text{Var} = \frac{n(L-1)(L\theta + n)}{L^2(L\theta + 1)}$

Numerical verification shows relative error between theory and simulation less than 2.5%.

Phase Diagram Analysis

The following main phase regions are identified:

1. Strong condensation phase (θ≪1): All particles aggregate at a single site
2. Weak condensation phase (θ≫1, uniform distribution): Effective aggregation due to combinatorial factors
3. Flat phase (θt≫1): Particles distributed as uniformly as possible
4. Uniform phase (θt=1): All configurations equally probable

Special Findings

Beta-Binomial Distribution

At t=1, single-point occupation probability strictly follows a beta-binomial distribution: $\langle\eta_i = \alpha\rangle = \binom{n}{\alpha} \frac{B(\alpha + \theta, n - \alpha + (L-1)\theta)}{B(\theta, (L-1)\theta)}$

Palindromic Symmetry

The stationary weight $w^t_{a,t}(\eta)$ as a polynomial in a and t is palindromic when n is even and anti-palindromic when n is odd.

Related Work

Historical Development

1. ASEP foundations: Classical exclusion process theory established by Spitzer et al.
2. Inclusion processes: Systematic study of SIP by Giardinà-Redig-Vafayi
3. ASIP development: ASIP proposed by Grosskinsky-Redig-Vafayi with condensation phenomena studies
4. Misanthrope processes: General theoretical framework established by Cocozza-Thivent

Innovations of This Paper

Compared to existing work, this paper achieves through t-analogue generalization:

• More general mathematical framework
• New exactly solvable models
• Rich symmetry structures
• Unified theoretical description

Conclusions and Discussion

Main Conclusions

1. The (q, t, θ) ASIP is exactly solvable: Complete analytical solutions obtained through misanthrope process theory
2. Stationary distribution independent of asymmetry: Stationary distribution does not depend on parameter q
3. Rich phase structure: Parameter space exhibits multiple physical phases and phase transitions
4. Deep mathematical structure: Discovery of palindromic symmetries and enriched processes

Limitations

1. One-dimensional restriction: Model applicable only to one-dimensional systems
2. Periodic boundaries: Only periodic boundary conditions considered
3. Numerical challenges: Simulation faces numerical stability issues at large t values
4. Phase transition understanding: Physical mechanisms in certain phase regions require deeper investigation

Future Directions

1. Higher-dimensional generalization: Extension to two-dimensional or higher-dimensional lattices
2. Open boundaries: Study behavior under open boundary conditions
3. Dynamical properties: Investigation of relaxation times and dynamical phase transitions
4. Experimental realization: Search for possible physical implementations

In-Depth Evaluation

Strengths

1. Mathematical rigor: All results supported by rigorous mathematical proofs
2. Theoretical completeness: Complete theoretical system from basic definitions to complex properties
3. Computational feasibility: Provides computable explicit formulas
4. Physical insights: Reveals rich physical phenomena and phase structures
5. Numerical verification: Comprehensive numerical simulations support theoretical predictions

Weaknesses

1. Limited applicability: Primarily theoretical research with unclear practical applications
2. Complexity: Introduction of t-analogues increases model complexity
3. Numerical difficulties: Simulation challenges in certain parameter regions
4. Physical interpretation: Physical meaning of some mathematical results requires further clarification

Impact

1. Theoretical contribution: Provides new examples and methods for integrable systems theory
2. Mathematical value: Discoveries such as palindromic symmetries have independent mathematical significance
3. Methodological: t-analogue generalization method applicable to other models
4. Educational value: Excellent example of exactly solvable models

Applicable Scenarios

1. Statistical physics: Non-equilibrium statistical mechanics theory research
2. Mathematical physics: Integrable systems and exact solution methods
3. Probability theory: Markov processes and stochastic process theory
4. Combinatorics: q-analogues and special functions research

References

The paper cites 31 important references covering:

• Classical ASEP theory (Spitzer, Liggett, et al.)
• Inclusion process research (Giardinà, Redig, Vafayi, et al.)
• Misanthrope process theory (Cocozza-Thivent, Evans, Waclaw, et al.)
• Special functions and q-analogues (Gasper, Rahman, et al.)

This paper demonstrates excellence in theoretical depth and mathematical rigor, providing a new exactly solvable model for non-equilibrium statistical mechanics with significant theoretical value and methodological importance.