Modified Rational Six Vertex Model on a Rectangular Lattice: New Formula, Homogeneous and Thermodynamic Limits

Basic Information

• Paper ID: 2509.11797
• Title: Modified rational six vertex model on a rectangular lattice: new formula, homogeneous and thermodynamic limits
• Authors: Matthieu Cornillault, Samuel Belliard (Institut Denis Poisson, CNRS-UMR 7013, Université de Tours)
• Classification: math-ph (Mathematical Physics), cond-mat.stat-mech (Statistical Mechanics), math.MP (Mathematical Physics), quant-ph (Quantum Physics)
• Publication Date: October 15, 2025 (arXiv v2)
• Paper Link: https://arxiv.org/abs/2509.11797v2

Abstract

This paper extends the work of Belliard, Pimenta, and Slavnov (2024) on the modified rational six vertex model. The authors discover an alternative formula for the partition function of the non-homogeneous model, based on expressions mixing modified Izergin determinants and Vandermonde determinants. This expression enables the calculation of the partition function for rectangular lattices in the homogeneous limit, thereby facilitating the study of the thermodynamic limit. The research yields a new result: the leading-order term of the free energy including boundary effects in the thermodynamic limit.

Research Background and Motivation

Problem Background

The six vertex model is a classical integrable model in statistical mechanics, closely related to quantum spin chains. This model is crucial for understanding phase transitions, critical phenomena, and quantum many-body systems. Specifically:

1. Model Development History: Beginning from XXX or XXZ models on square lattices with domain wall boundary conditions (DWBC), whose partition function is represented by the celebrated Izergin determinant
2. Generalization of Boundary Conditions: Introduction of partial domain wall boundary conditions (pDWBC) and the Kostov determinant related to N=4 supersymmetric Yang-Mills theory
3. General Boundary Conditions: The modified rational six vertex (MR6V) model introduces more general boundary conditions, where boundaries no longer contain purely up or down spins, but rather linear combinations of up and down spins

Research Motivation

1. Theoretical Completeness: Existing partition function expressions are primarily designed for square lattices; the rectangular lattice case requires new mathematical tools
2. Thermodynamic Limit: Understanding how boundary effects influence macroscopic properties of the system in the thermodynamic limit
3. Integrability Studies: Establishing connections between modified Bethe vectors constructed via modified algebraic Bethe ansatz (MABA) and the partition function

Limitations of Existing Methods

• Known modified Izergin determinant formulas are primarily designed for square lattices
• Lack of effective methods for handling the thermodynamic limit of rectangular lattices
• The behavior of boundary effects in the thermodynamic limit remains unclear

Core Contributions

1. New Partition Function Formula: Proposes a new determinantal expression applicable to n×m rectangular lattices, mixing modified Izergin determinants and Vandermonde determinants
2. Homogeneous Limit Calculation: Successfully computes the partition function in the homogeneous limit using the new formula
3. Thermodynamic Limit Analysis: Systematically studies two thermodynamic limit scenarios: semi-infinite and infinite lattices
4. Quantitative Description of Boundary Effects: First obtains the free energy expression including boundary effects in the thermodynamic limit
5. Correspondence with pDWBC: Establishes explicit connections between the MR6V model and models with partial domain wall boundary conditions

Methodology Details

Problem Definition

Study the modified rational six vertex model on an n×m rectangular lattice with general boundary conditions, where:

• Input: Lattice size parameters n, m, non-homogeneity parameters ū={u₁,...,uₙ} and v̄={v₁,...,vₘ}, boundary condition parameters B and Ĉ
• Output: Partition function Zₙₘ(ū|v̄|B|Ĉ) and its behavior in various limits
• Constraint: Maintain integrability of the Yang-Baxter equation

Core Mathematical Framework

1. R-Matrix and Statistical Weights

The R-matrix in the rational case is:

R_{ab}(u,v) = (u-v)/c · I + P_{ab}

with statistical weights:

• a(u,v) = (u-v+c)/c
• b(u,v) = (u-v)/c
• c(u,v) = 1

2. General Boundary Conditions

Boundary conditions are defined by four vectors:

|x⟩ = x₁|1⟩ + x₂|2⟩, x ∈ {n,s,e,w}

with twist operators:

B = |e⟩⊗⟨w|, Ĉ = |s⟩⊗⟨n|

3. New Partition Function Formula

Main Result (Proposition 1.3):

Zₙₘ(ū|v̄|B|Ĉ) = tr(B)ⁿtr(Ĉ)ᵐ (1-β)^{min(n,m)} h(ū,v̄)/(g(ū,v̄)Δ(ū)Δ'(v̄))

multiplied by the following determinant:

• When n≤m: det_{1≤i,j≤m}({φ_β(uᵢ-vⱼ) if i≤n; ψ_(-vⱼ) if i≥n+1})
• When n≥m: det_{1≤i,j≤n}({φ_β(uᵢ-vⱼ) if j≤m; ψ_(uᵢ) if j≥m+1})

where β = 1 - tr(B)tr(Ĉ)/tr(BĈ)

Technical Innovations

1. Mixed Determinant Structure

The key innovation of the new formula is decomposing the max(n,m)×max(n,m) determinant into two parts:

• First min(n,m) rows/columns: Modified Izergin determinant type
• Subsequent rows/columns: Vandermonde determinant type

2. Partial Cauchy Matrix Method

Introduces partial Cauchy matrices Cₙₘ(ū|v̄) and their inverses, establishing equivalence proofs between different partition function expressions.

3. Generalization of the Foda-Wheeler Method

Adopts and generalizes the method developed by Foda and Wheeler for pDWBC, deriving results for rectangular lattices from larger square lattices through limiting procedures.

Experimental Setup

Mathematical Verification Framework

This paper primarily conducts theoretical analysis and mathematical derivations, with verification methods including:

1. Formula Equivalence Proofs: Proving the correctness of the new formula through two independent methods
2. Limiting Behavior Analysis: Systematically studying behavior under various parameter limits
3. Special Case Verification: Verifying consistency of the new formula with known special cases

Analysis Methods

1. Inductive Proofs: Used to prove determinantal properties of partial Cauchy matrices
2. Complex Analysis Techniques: Utilizing residue theorem to calculate properties of complex rational functions
3. Asymptotic Analysis: Studying asymptotic behavior in the thermodynamic limit

Main Results

1. Homogeneous Limit Results

When uᵢ→u, vⱼ→v, the partition function simplifies to:

Zₙₘ(x|B|Ĉ) = tr(B)ⁿtr(Ĉ)ᵐ (1-β)^{min(n,m)} 
× det_{1≤i,j≤min(n,m)}((d+i+j-2 choose i-1)((1+x/c)^{d+i+j-1} - β(x/c)^{d+i+j-1}))

where d = |n-m|, x = u-v.

2. Completely Homogeneous Lattice Physical Interpretation

The partition function for a completely homogeneous lattice (x=0) is:

Z⁰ₙₘ(B|Ĉ) = tr(BĈ)^{min(n,m)} × {tr(Ĉ)^d if n≤m; tr(B)^d if n≥m}

The total free energy contains volume and boundary effects:

F^{tot}_{nm} = nmε - k_BT(min(n,m)ln(tr(BĈ)) + d×boundary term)

3. Thermodynamic Limit Analysis

Semi-Infinite Lattice Limit

When max(n,m)→+∞:

Zₙₘ(x|B|Ĉ) ∼ Z⁰ₙₘ(B|Ĉ)(1+x/c)^{nm}

The asymptotic behavior resembles all vertices being of type a (lowest energy state).

Infinite Lattice Limit

When n,m→+∞ with d=|n-m| held constant, the bulk free energy satisfies the Toda equation:

f''(x) = exp(2f(x))

The solution has the form:

exp(f(x)) = α/(sinh(α(x-x₀)))

Depending on the value of d, the expression for α is:

• d>1: α=0
• d=1: α²=3β/c²
• d=0: α²=3β(β-2)/c²

4. Quantitative Description of Boundary Effects

The final expression for bulk free energy is:

F̃(x̃) = {
  ln(sinh(√β̃x̃)/(√β̃x̃) × 1/(1+x̃))  if β̃>0
  ln(1/(1+x̃))                        if β̃=0  
  ln(sin(√-β̃x̃)/(√-β̃x̃) × 1/(1+x̃))   if β̃<0
}

where x̃=x/c, and β̃ depends on the values of d and β.

Related Work

1. Izergin Determinant Theory: Fundamental representation of six vertex model partition functions
2. Partial Domain Wall Boundary Conditions: Pioneering work by Foda and Wheeler
3. Modified Algebraic Bethe Ansatz: Theoretical framework by Belliard et al.
4. Thermodynamic Limit Studies: Classical results by Korepin and Zinn-Justin

Conclusions and Discussion

Main Conclusions

1. Successfully derives a new partition function formula for the MR6V model on rectangular lattices
2. Provides complete analysis of the progression from non-homogeneous to homogeneous to thermodynamic limits
3. Discovers significant contributions of boundary effects to free energy in the thermodynamic limit
4. Establishes explicit correspondence with pDWBC models

Limitations

1. The asymptotic approximation (2.24) may neglect important effects of strong boundary fields
2. Higher-order correction terms are needed for more complete descriptions
3. Currently limited to the rational case; trigonometric and elliptic cases require further investigation

Future Directions

1. Generalization to Other Models: Trigonometric six vertex model (XXZ spin chain) and eight vertex model (XYZ spin chain)
2. Elliptic Case: Seeking similar formulas for the elliptic case based on the existence of elliptic source identities
3. Phase Transition Studies: Investigating free energy behavior in different phases, particularly antiferroelectric phases and arctic curve phenomena
4. Reflection Boundary Conditions: Applying linear system methods to study partition functions with reflection boundary conditions

In-Depth Evaluation

Strengths

1. Mathematical Rigor: Provides two independent proof methods, ensuring reliability of results
2. Methodological Innovation: The partial Cauchy matrix method and mixed determinant structure are original
3. Result Completeness: Complete analytical chain from partition function formulas to thermodynamic properties
4. Physical Insight: Reveals the important role of boundary effects in the thermodynamic limit
5. Clarity of Presentation: Well-structured paper with detailed mathematical derivations

Weaknesses

1. Computational Complexity: The new formula has high computational complexity, potentially limiting practical applications
2. Scope of Applicability: Currently applicable only to the rational case; generalization requires additional work
3. Numerical Verification: Lacks numerical calculations to verify theoretical results
4. Physical Applications: Connections to real physical systems require further clarification

Impact

1. Theoretical Contribution: Provides new mathematical tools for integrable systems theory
2. Method Value: The partial Cauchy matrix method may be applicable to other related problems
3. Foundation for Future Research: Establishes foundation for investigations of trigonometric and elliptic cases
4. Interdisciplinary Value: Connects mathematical physics, statistical mechanics, and quantum information fields

Applicable Scenarios

1. Theoretical Research: Theoretical analysis of integrable systems, quantum spin chains, and statistical mechanics models
2. Numerical Computation: Provides theoretical foundation for large-scale system numerical simulations
3. Phase Transition Studies: Analysis of critical phenomena and phase transition behavior
4. Quantum Information: Investigation of quantum entanglement and quantum phase transitions

References

The paper cites 44 important references covering multiple fields including integrable systems, statistical mechanics, and mathematical physics, particularly:

• Baxter's "Exactly Solved Models in Statistical Mechanics"
• Pioneering work by Izergin et al. on six vertex models
• Recent research by Belliard et al. on modified Bethe ansatz
• Important contributions by Foda and Wheeler on partial domain wall boundary conditions