Steady state representations for the harmonic process

This paper establishes the previously unavailable matrix product solution for the steady state of the harmonic process by demonstrating its derivation from and clarifying its relationship to existing closed-form and nested integral representations found in recent literature.

The Gist

Imagine a long, narrow hallway with $N$ rooms. In each room, people can gather. In some versions of this story, a room can hold only 0 or 1 person. But in the story told in this paper, a room can hold any number of people—10, 100, or even 1,000. There is no limit.

People move between rooms randomly, and people also enter or leave the hallway from the two ends (the left door and the right door). Over time, the crowd settles into a specific pattern where the number of people in each room stops changing on average. This is called the steady state.

The author of this paper, Rouven Frassek, is a detective trying to solve a mystery: How do we mathematically describe this final crowd pattern?

The Three Ways to Describe the Crowd

Before this paper, scientists knew two ways to describe the final crowd, but both were tricky:

1. The "Closed-Form" Recipe: This is like a complex, multi-step cooking recipe. It gives you the exact answer if you follow the instructions perfectly, but the instructions are long and involve huge numbers (factorials and gamma functions) that are hard to visualize.
2. The "Nested Integral" Recipe: This is like a Russian nesting doll. To find the answer, you have to solve a problem, which gives you a new problem, which gives you another, and so on. You have to solve them one by one, from the inside out. It's elegant but computationally heavy.

The Missing Piece:
Scientists had a third, very powerful tool called the Matrix Product Ansatz (MPA). Think of this as a "Lego instruction manual." Instead of calculating a giant number for the whole hallway, you build the answer by snapping together small, simple blocks (matrices) for each room. If you have the right blocks, you can instantly see how the whole system behaves.

The Problem: No one knew what the "Lego blocks" looked like for this specific "unlimited capacity" hallway. The math was too messy.

The Author's Solution: The Magic Translator

Frassek's paper is the breakthrough that finally provides those Lego blocks. Here is how he did it, using a simple analogy:

1. The "Magic Translator" (Similarity Transformation)

The original problem (the hallway with people moving randomly) is chaotic and hard to solve directly. Frassek realized that if you apply a "magic translator" (a mathematical rotation), you can turn the chaotic hallway into a simpler, quieter version.

In this simpler version:

• The rules for people moving between rooms stay the same.
• But the rules for people entering and leaving the doors become much, much simpler.

It's like taking a noisy, crowded party and putting on noise-canceling headphones that only let you hear the music, not the shouting. The core structure is the same, but the math becomes manageable.

2. Building the Lego Blocks

Once the problem was simplified, Frassek could finally see the pattern. He realized the "Lego blocks" (the matrices) could be built using oscillators.

• Analogy: Imagine a set of musical instruments. Some instruments add a note (create a person), and some remove a note (remove a person).
• Frassek showed that the "blocks" for each room are just specific combinations of these musical notes.
• He proved that if you snap these blocks together in a line, they perfectly recreate the "nested doll" solution and the "complex recipe" solution.

3. The "Hidden Helper" (The Auxiliary Operator)

To make the Lego blocks snap together correctly, Frassek had to invent a "hidden helper" tool (an auxiliary operator). This tool doesn't appear in the final picture of the crowd, but it's essential for the blocks to fit together without falling apart. It's like a scaffolding used to build a bridge; you don't see the scaffolding in the finished bridge, but the bridge wouldn't stand without it.

Why Does This Matter?

1. It Unifies the World: The paper shows that the "Complex Recipe," the "Nesting Dolls," and the "Lego Blocks" are all just different languages describing the exact same reality. They are all connected.
2. It Opens New Doors: Now that we have the Lego blocks, we can easily calculate things that were previously impossible, like how the crowd fluctuates or how long it takes to reach the steady state.
3. It Solves a Long-Standing Puzzle: For years, physicists knew this model existed but couldn't find the "Lego" version for it. This paper fills that gap, showing that even for systems with infinite capacity, there is a neat, structured way to describe them.

The Bottom Line

Rouven Frassek took a messy, complicated problem about people moving in an infinite hallway, simplified it with a mathematical "magic trick," and discovered the hidden "Lego instructions" that describe the system. This allows scientists to build and understand these complex systems much more easily than before.

Technical Summary

1. Problem Statement

The paper addresses the challenge of finding a Matrix Product Ansatz (MPA) solution for the steady state of the harmonic process, a stochastic particle system defined on a one-dimensional lattice with open boundaries.

• Context: The harmonic process is the symmetric (rational) limit of integrable particle processes with unbounded occupation numbers per site. It arises from the nearest-neighbor Hamiltonian of an open Heisenberg XXX spin chain with non-compact $sl(2)$ representations.
• The Gap: While the steady state of this system was previously known in two forms:
  1. A closed-form expression derived via the Quantum Inverse Scattering Method (QISM) using non-local charges [1].
  2. A nested integral representation interpreted as a mixed measure [2, 3].
     ...a Matrix Product Solution (MPS) had remained elusive. The difficulty stems from the unbounded configuration space, which implies the underlying algebra requires infinitely many generators, unlike bounded models like the SSEP (Symmetric Simple Exclusion Process) or ASEP.

2. Methodology

The author employs a multi-step strategy to bridge the gap between known solutions and the desired matrix product form:

• Similarity Transformation: The core methodological innovation is the use of a local similarity transformation to map the stochastic Hamiltonian $H$ (which has complex boundary terms) to a transformed Hamiltonian $\tilde{H}$.
  • The transformation is defined as $\tilde{H} = e^{-\rho_R S_{tot}^+} e^{S_{tot}^-} H e^{-S_{tot}^-} e^{\rho_R S_{tot}^+}$.
  • This transformation simplifies the boundary terms significantly (making the right boundary trivial) while leaving the bulk dynamics unchanged, allowing the steady state $|\nu\rangle$ of $\tilde{H}$ to be found more easily.
• Derivation from Closed-Form: The author constructs the matrix product operators directly from the known closed-form expression of the steady state components $\nu(\vec{m})$. By introducing a pair of oscillators ($a, \bar{a}$) and a Fock space, the telescopic product structure of the closed-form solution is mapped to matrix elements of operators acting on this auxiliary space.
• Derivation from Integral Representation: The paper demonstrates that the nested integral representation can also be converted into a matrix product form. This is achieved by defining integral operators acting on monomials, which are then shown to coincide with the matrix elements derived from the oscillator approach.
• Algebraic Verification: The author explicitly constructs the auxiliary operator $\bar{Y}$ and verifies that the proposed operators satisfy the required bulk and boundary algebraic relations (generalized DEHP algebra) necessary for the MPA to be a valid solution.

3. Key Contributions

• First Matrix Product Solution: The paper provides the first explicit Matrix Product Ansatz for the steady state of the harmonic process with open boundaries.
• Unification of Representations: It rigorously establishes the equivalence between three distinct representations of the steady state:
  1. The closed-form expression (QISM based).
  2. The nested integral representation.
  3. The newly derived matrix product form.
• Explicit Algebraic Structure: The paper derives the specific infinite-dimensional matrix product algebra satisfied by the generators $X(m)$ (and their auxiliary counterparts $\bar{X}$). This includes explicit bulk commutation relations and complex boundary conditions involving the reservoir parameters $\beta_{L,R}$.
• Construction of Generators: The bulk generators $X(m)$ are explicitly constructed in terms of creation/annihilation operators and Gamma functions, involving a similarity rotation of the simpler operators $Y(m)$ found in the transformed frame.

4. Results

• Steady State Form: The steady state $|\mu\rangle$ is expressed as:
  $$|\mu\rangle = Z_N^{-1} \sum_{m_1, \dots, m_N} \langle\langle V | X(m_1) \cdots X(m_N) | W \rangle\rangle |m_1, \dots, m_N\rangle$$
  where $X(m)$ are the bulk generators, and $|W\rangle\rangle, \langle\langle V|$ are boundary states in the auxiliary space.
• Operator Definitions:
  • The auxiliary operators $Y$ (related to the transformed Hamiltonian) are constructed using oscillator algebra: $Y(m) = \kappa(m) \bar{a}^{2s} \frac{\Gamma(N+1)}{\Gamma(N+2s+1)} \bar{a}^m$.
  • The physical operators $X$ are obtained via a rotation: $X = e^{-S_-} e^{\rho_R S_+} Y$.
• Algebraic Relations: The paper proves that these operators satisfy the relation:
  $$H(Y \otimes Y) = Y \otimes \bar{Y} - \bar{Y} \otimes Y$$
  along with specific boundary conditions involving the transformed boundary operators $\tilde{B}_{L,R}$.
• Normalization: The normalization constant $Z_N$ is explicitly calculated, matching the prefactors found in the closed-form and integral representations.

5. Significance

• Theoretical Advancement: This work extends the applicability of the Matrix Product Ansatz to systems with unbounded occupation numbers and non-compact spin representations, a regime where finding finite-dimensional representations is impossible and infinite-dimensional ones are notoriously difficult.
• Methodological Insight: It demonstrates that similarity transformations can be used to "triangularize" stochastic Hamiltonians, simplifying the search for MPA solutions in complex integrable systems.
• Future Directions: The paper suggests that the integral representation of the R-matrix (known for non-compact representations) could be used to derive these matrix product states directly, potentially unifying the integral and matrix approaches. It also opens the door to studying the probabilistic structure of these states and extending the results to fractional spin values ($2s > 0$).
• Broader Impact: By clarifying the relationship between QISM-derived closed forms, integral measures, and matrix product states, this work provides a template for solving other integrable stochastic processes with infinite-dimensional configuration spaces.