Inhomogeneous SSH models and the doubling of orthogonal polynomials

This paper demonstrates that the doubling method for orthogonal polynomials provides an analytical framework for solving Su-Schrieffer-Heeger models, linking the standard case to Chebyshev polynomials and extending exact solvability to inhomogeneous variants via Krawtchouk and $q$-Racah polynomials.

The Gist

Imagine you are building a long, winding bridge made of stepping stones. In a standard physics model called the SSH model (named after Su, Schrieffer, and Heeger), this bridge has a very specific rhythm: a strong step, a weak step, a strong step, a weak step, all the way down the line.

Physicists love this bridge because it's a perfect playground for studying "topology"—a fancy way of describing how things are connected in a way that can't be easily untangled. But usually, this bridge is perfectly uniform. Every strong step is exactly the same strength, and every weak step is exactly the same.

The Big Idea: Making the Bridge Irregular (But Solvable)

This paper asks a bold question: What happens if we make the bridge messy? What if the strong steps get stronger as you go, and the weak steps get weaker? Or what if they change in a complex, wavy pattern?

Usually, making a system "inhomogeneous" (messy or uneven) makes it impossible to solve mathematically. It's like trying to predict the path of a ball rolling down a bumpy, random hill; the math gets too messy to write down a clean answer.

The authors of this paper found a magical "cheat code" to solve these messy bridges. They used a mathematical trick called Doubling, which relies on a family of special shapes called Orthogonal Polynomials.

The Magic Trick: The "Doubling" Recipe

Think of Orthogonal Polynomials (like Chebyshev, Krawtchouk, or q-Racah) as a set of perfectly tuned musical scales. Each scale has a specific pattern of notes.

The authors discovered that if you take two of these scales and "double" them together—interweaving them in a very specific way—you can create a new, complex pattern. This new pattern acts like a blueprint for building a messy bridge (an inhomogeneous SSH model) that, surprisingly, still has a perfect, clean solution.

Here is how the analogy works:

1. The Standard Bridge (Chebyshev):
   The original, uniform bridge is built using a specific musical scale called Chebyshev polynomials. The authors showed that the "notes" (energy levels) of this bridge are just the roots of these polynomials. It's like knowing that a perfect piano scale always has specific frequencies.

2. The Messy Bridges (Krawtchouk & q-Racah):
   The authors then asked, "What if we use other musical scales?"

   • They tried the Krawtchouk scale. This created a bridge where the step strengths change in a way that looks like a bell curve (getting stronger in the middle, weaker at the ends).
   • They tried the q-Racah scale. This created a bridge with a more exotic, quantum-mechanical pattern of steps.

   Even though these bridges look chaotic and uneven, the "Doubling" method allowed the authors to write down the exact energy levels and the exact path a particle would take on them. They didn't have to guess; they had the recipe.

Why Does This Matter?

You might wonder, "Who cares about messy bridges?"

In the real world, we can't always build perfect, uniform bridges.

• Cold Atoms: Scientists use lasers to trap atoms and make them behave like electrons on a bridge. They can control the "strength" of the links between atoms with incredible precision. They can intentionally make the links get stronger or weaker in a specific pattern.
• Photonic Lattices: Light can be guided through glass fibers arranged in a grid. Engineers can tweak the glass to make the light hop unevenly.

This paper provides the instruction manual for these experiments. If an engineer wants to build a bridge where the steps follow a Krawtchouk pattern, this paper tells them exactly what the energy levels will be and where the "special" particles (called Zero Modes) will hide.

The "Ghost" Particle (Zero Mode)

One of the coolest features of these bridges is a "ghost" particle. In the standard bridge, if you have an odd number of steps, there is a special state where a particle sits perfectly still at one end of the bridge, unaffected by the rest. This is a Zero Mode.

The authors showed that in their messy, inhomogeneous bridges, this ghost particle doesn't just sit at the end. It sits exactly where the "strong" steps and "weak" steps meet in the middle. It's like a hiker who stops exactly where the terrain changes from a steep climb to a gentle slope. The math predicts exactly where this hiker will stand, no matter how complex the terrain is.

Summary

• The Problem: Physics models usually break when you make them messy (inhomogeneous).
• The Solution: The authors used a mathematical "doubling" trick involving special polynomial shapes (like musical scales) to build messy models that are still solvable.
• The Result: They created a new family of "messy" bridges (SSH models) based on Krawtchouk and q-Racah polynomials.
• The Payoff: This gives experimentalists a blueprint to build and understand complex quantum systems in the lab, knowing exactly how they will behave without needing a supercomputer to guess.

In short, they took a simple, rhythmic dance (the SSH model), taught it how to dance to complex, irregular jazz (the inhomogeneous models), and proved that the dancers still know the steps perfectly.

Technical Summary

1. Problem Statement

The Su–Schrieffer–Heeger (SSH) model is a fundamental paradigm in condensed matter physics for studying topological phases, symmetry-protected edge states, and topological phase transitions. While the standard homogeneous SSH model (with constant alternating bond strengths) is exactly solvable using free-fermion techniques, extending this model to include spatial inhomogeneities (where bond strengths vary along the chain) generally renders the system analytically intractable.

The authors address the challenge of constructing exactly solvable inhomogeneous SSH models. Specifically, they aim to determine the coupling constants $t_n^+$ and $t_n^-$ for a chain of $2N+1$ sites such that the Hamiltonian remains diagonalizable with closed-form expressions for the energy spectrum and eigenstates.

2. Methodology: The Doubling Procedure

The core methodological innovation is the application of the doubling procedure for sequences of orthogonal polynomials. This mathematical technique combines two families of polynomials to generate a new sequence that satisfies a specific recurrence relation corresponding to a tridiagonal matrix (the Hamiltonian).

• Orthogonal Polynomials: The authors utilize finite families of orthogonal polynomials from the (q-)Askey scheme (e.g., Chebyshev, Krawtchouk, q-Racah). These polynomials satisfy a three-term recurrence relation:
  $$x P_n(x) = A_n P_{n+1}(x) - (A_n + C_n) P_n(x) + C_n P_{n-1}(x)$$
• The Doubling Construction:
  1. Renormalization: The polynomials $P_n$ are renormalized into a sequence $R_n$ to symmetrize the recurrence coefficients.
  2. Ansatz for Eigenvectors: A new sequence of polynomials $Q_n(x)$ is constructed by "doubling" the index of the original sequence. The ansatz is defined as:
     • $Q_{2n}(x) = t_n^+ R_n(\pi x) + t_{n-1}^- R_{n-1}(\pi x)$
     • $Q_{2n+1}(x) = x R_n(\pi x)$
       Here, $\pi x = \tau_2 x^2 + \tau_0$ is an affine transformation of the spectral variable.
  3. Constraint Matching: By substituting the ansatz into the eigenvalue equation $H \mathbf{Q} = x \mathbf{Q}$, the authors derive a set of algebraic constraints linking the coupling constants ($t_n^\pm$) of the SSH model to the recurrence coefficients ($A_n, C_n$) of the chosen orthogonal polynomials.
  4. Spectral Determination: The eigenvalues are determined by solving the condition $\pi x = \lambda_k$, where $\lambda_k$ are the roots (grid points) of the underlying orthogonal polynomials.

3. Key Contributions and Results

A. Re-derivation of the Standard SSH Model

The authors demonstrate that the standard homogeneous SSH model corresponds to the doubling of Chebyshev polynomials of the second kind ($U_n$).

• Result: This provides a new mathematical framework for the standard model, showing that its spectrum and eigenvectors are naturally derived from the properties of $U_n$.
• Extension: They generalize this to include a non-vanishing chemical potential ($\mu_\pm$) that alternates with site parity, linking it to a modified argument in the Chebyshev polynomials.

B. Construction of Inhomogeneous Models

The primary contribution is the generalization of the doubling procedure to arbitrary finite orthogonal polynomial sequences, yielding new exactly solvable inhomogeneous SSH models. The paper details two specific cases:

1. Krawtchouk-Type SSH Model:

   • Polynomials: Krawtchouk polynomials $K_n(x; p, N)$.
   • Couplings: The bond strengths $t_n^\pm$ vary along the chain according to square roots of linear functions of the site index $n$.
   • Spectrum: The energy spectrum is given by $x_k^\pm = \pm \sqrt{k+1}$ (plus a zero mode).
   • Zero Mode: The zero-energy mode is localized not at the boundary, but at the site $n \approx pN$ where the coupling strengths $t_n^+ \approx t_n^-$. This represents a topological phase transition point within the chain.

2. q-Racah-Type SSH Models:

   • Polynomials: q-Racah polynomials, the most general family in the Askey scheme.
   • Two Distinct Models:
     • Model 1 (Odd sites): Uses a specific contiguity relation to construct a model on $2N+1$ sites. The spectrum involves q-deformed terms: $x_k^\pm = \pm \sqrt{(1-q^{k-N})(\delta - q^{-k})}$.
     • Model 2 (Even sites): By utilizing a different contiguity relation that does not shift the degree $N$, the authors construct a solvable model on $2N$ sites (an even number), which is typically difficult for SSH models due to the absence of a zero mode in the standard formulation.
   • Significance: These models exhibit highly non-trivial, spatially modulated hopping profiles that are analytically tractable.

C. Zero Modes and Topology

The paper rigorously analyzes the zero mode (zero-energy eigenstate).

• In the inhomogeneous case, the zero mode is localized where the local "topological mass" changes sign (i.e., where $t_n^+ \approx t_n^-$).
• This confirms that the inhomogeneous models support topological edge-like states that can be tuned to appear anywhere in the bulk of the chain, depending on the profile of the couplings.

4. Significance and Implications

• Analytical Tractability: The work provides a systematic "recipe" to generate infinite families of exactly solvable inhomogeneous quantum many-body systems. This is rare, as inhomogeneity usually breaks integrability.
• Experimental Relevance: While the specific coupling profiles are mathematical constructs, they are highly relevant for synthetic quantum systems (e.g., cold atoms in optical lattices, photonic crystals, and acoustic lattices) where hopping amplitudes can be engineered with high precision. These models serve as testbeds for studying topological phases in disordered or engineered environments.
• Topological Insights: The ability to analytically solve these models allows for the precise calculation of entanglement measures, correlation functions, and the behavior of zero modes under spatial modulation, which are crucial for understanding topological quantum computing and phase transitions.
• Future Directions: The authors suggest this framework can be extended to other models (e.g., the Kitaev chain) and higher-dimensional lattices using bivariate orthogonal polynomials.

Conclusion

Crampé et al. successfully bridge the gap between the theory of orthogonal polynomials and topological quantum matter. By leveraging the doubling procedure, they transform the problem of finding inhomogeneous couplings into a problem of selecting appropriate polynomial families. This yields a rich class of exactly solvable SSH models with tunable topological properties, offering powerful tools for theoretical analysis and experimental design in topological quantum systems.