The optical Su-Schrieffer-Heeger model on a triangular lattice

Using determinant quantum Monte Carlo simulations, this study reveals that the triangular lattice optical Su-Schrieffer-Heeger model exhibits distinct doping-dependent behaviors, including a symmetry-breaking bond-order-wave transition at half-filling and a potential $s$-wave superconducting phase at three-quarters filling for large phonon energies, while showing no enhanced magnetic correlations.

The Gist

Imagine a crowded dance floor where the dancers are electrons and the floor itself is made of springs (the atoms in a crystal lattice). Usually, we think of the dancers moving around a rigid floor. But in this study, the floor is alive: when a dancer steps on a spring, the spring stretches or compresses, changing the shape of the dance floor for everyone else.

This is the world of the Su-Schrieffer-Heeger (SSH) model. The researchers in this paper decided to put this dance floor on a triangular grid (like a honeycomb or a triangle of triangles) instead of the usual square grid, and they asked: How does the floor's movement change the way the dancers behave?

Here is the breakdown of their findings, translated into everyday language:

1. The Setup: A Triangular Dance Floor

In physics, the shape of the grid matters a lot.

• Square Grid: Like a checkerboard. It's orderly.
• Triangular Grid: Like a honeycomb. It's "frustrated." Imagine trying to sit at a round table with three friends; if you all want to face each other, it's impossible to please everyone at once. This "frustration" makes the physics much more interesting and chaotic.

The researchers used a super-powerful computer simulation (called Quantum Monte Carlo) to watch millions of these "dance moves" happen, looking for patterns.

2. The Two Special Crowd Densities

They tested two different crowd sizes on the dance floor:

Scenario A: The Half-Full Floor (1/4 Filling)

Imagine the dance floor is only half full.

• What happened: The dancers started organizing themselves into a rigid pattern. They stopped moving freely and started locking arms in a specific way, creating a "Bond-Order Wave" (BOW).
• The Analogy: Think of a crowd of people in a hallway. Suddenly, everyone decides to stop walking and form a perfect, alternating line of "left-right-left-right" hands. They freeze into a solid, insulating state.
• The Twist: On a square floor, this pattern is usually simple. But on this triangular floor, the pattern broke a specific symmetry (called $C_6$). It's like the dancers suddenly decided to all face only one specific corner of the room, breaking the perfect hexagonal balance of the floor.

Scenario B: The Almost-Full Floor (3/4 Filling)

Now, imagine the dance floor is almost packed, with only a few empty spots.

• What happened: This is where things got exciting. The outcome depended on how "heavy" or "fast" the springs (phonons) were.
  • Slow Springs (Adiabatic limit): If the floor moves slowly, the dancers again formed that rigid, frozen pattern (BOW), similar to the half-full case.
  • Fast Springs (Anti-adiabatic limit): If the floor vibrates very quickly, the dancers did something magical: They paired up.
• The Analogy: Imagine the floor is shaking so fast that the dancers can't walk in a straight line. Instead, they grab hands and start dancing in pairs, spinning together. In physics, this is Superconductivity (specifically s-wave). The pairs can move through the crowd without any friction or resistance.

3. The "Sign Flip" Secret

Why did the fast springs cause pairing?
The researchers found a clever trick. When the floor vibrates hard enough, the springs stretch so far that they actually flip the sign of the connection between dancers.

• Analogy: Imagine a rule that says "If you step here, you move forward." But if the spring stretches too far, the rule flips to "If you step here, you move backward." This sudden flip in the rules creates a perfect environment for the dancers to pair up and glide effortlessly.

4. What Didn't Happen (The Surprise)

In similar studies on square dance floors, scientists often see the dancers get angry and start fighting (forming magnetism or anti-ferromagnetism).

• The Result: On this triangular floor, the researchers found no fighting. The magnetic correlations were weak or non-existent.
• Why it matters: It shows that the triangular shape and the specific way the floor moves (optical phonons) completely suppress the "anger" (magnetism) that usually ruins superconductivity in other models.

Summary: The Big Picture

This paper is like a recipe book for a new kind of superconductivity.

• The Ingredients: Electrons, a triangular lattice, and a floor that vibrates.
• The Secret Sauce: If you fill the floor just right (3/4 full) and make the floor vibrate fast enough, the electrons will naturally pair up and become superconductors.
• The Takeaway: The shape of the world (triangular vs. square) and the speed of the vibrations are the keys to unlocking whether the electrons freeze into a solid block or flow like a frictionless super-current.

In short: By shaking a triangular dance floor just right, the researchers found a way to make electrons dance in perfect pairs, potentially paving the way for new types of superconducting materials.

Technical Summary

1. Problem Statement

The paper investigates the interplay between electrons and phonons in a triangular lattice using the optical Su-Schrieffer-Heeger (SSH) model. While electron-phonon (e-ph) interactions are known to stabilize various states of matter (e.g., superconductivity, bond-order waves, polarons), the specific effects of geometric frustration inherent to the triangular lattice remain less explored compared to bipartite lattices (like the square lattice).

Key questions addressed include:

• How does the triangular geometry, which introduces kinetic energy frustration, compete with e-ph coupling to determine the ground state?
• What phases emerge at different carrier concentrations (filling fractions), specifically at one-quarter ($\langle n \rangle = 0.5$) and three-quarters ($\langle n \rangle = 1.5$) filling?
• Does the model exhibit enhanced magnetic correlations (antiferromagnetism) similar to square lattice SSH models, or does frustration suppress them?
• How do the adiabatic (low phonon energy) and anti-adiabatic (high phonon energy) limits affect the competition between Bond-Order-Wave (BOW) and superconducting phases?

2. Methodology

The authors employed Determinant Quantum Monte Carlo (DQMC) simulations to solve the model non-perturbatively.

• Model Hamiltonian: The optical SSH model on a triangular lattice where atomic displacements modulate the nearest-neighbor hopping integrals ($t$). The Hamiltonian includes:
  • Kinetic energy with displacement-dependent hopping: $t(R_i - R_j) \approx t_0 + \nabla t \cdot (Q_i - Q_j)$.
  • Chemical potential ($\mu$).
  • Einstein phonon oscillators with frequency $\Omega$ and mass $M$.
  • Dimensionless coupling constant: $\lambda = \alpha^2 / (M\Omega^2 t)$.
• Simulation Details:
  • Lattice: $L \times L$ triangular lattices with periodic boundary conditions.
  • Algorithm: Hybrid Monte Carlo (HMC) updates implemented via the SmoQyDQMC.jl package.
  • Sign Problem: The model is free of the sign problem because the two spin species couple symmetrically to the phonon coordinates, resulting in equal real-valued determinants.
  • Parameters: Systematic variation of carrier concentration ($\langle n \rangle$), e-ph coupling ($\lambda$), and phonon energy ($\Omega$).
• Observables:
  • Superconductivity: $s$-wave pair-field susceptibility ($\chi_{sc}$) and structure factor.
  • Bond Order: Bond-order-wave (BOW) susceptibility ($\chi_b$) and structure factors for broken $C_2, C_3, C_6$ rotational symmetries.
  • Magnetism: Spin-spin correlation functions and structure factors.
  • Phase Transitions: Critical temperatures and couplings estimated using the correlation ratio method ($R_\Theta$) and finite-size scaling.

3. Key Contributions & Results

A. Phase Diagram at One-Quarter Filling ($\langle n \rangle = 0.5$)

• Fermi Surface: Corresponds to a circular non-interacting Fermi surface.
• Ground State: The system undergoes a transition from a metal to an insulating Bond-Order-Wave (BOW) phase as the e-ph coupling ($\lambda$) increases.
• Symmetry Breaking: The BOW phase breaks the local $C_6$ rotational symmetry of the lattice.
• Coupling Dependence: The critical coupling $\lambda_c$ required to induce the BOW phase increases linearly with phonon energy $\Omega$.
• Magnetism: No evidence of enhanced magnetic correlations; the system remains non-magnetic.

B. Phase Diagram at Three-Quarters Filling ($\langle n \rangle = 1.5$)

• Fermi Surface: Corresponds to a hexagonal non-interacting Fermi surface.
• Competing Phases:
  • Adiabatic Limit (Small $\Omega$): The system favors a BOW phase (though weaker than at $\langle n \rangle = 0.5$).
  • Anti-adiabatic Limit (Large $\Omega$): The system transitions to a robust $s$-wave superconducting phase.
• Mechanism: The tendency toward pairing in the anti-adiabatic limit is associated with the possibility of a sign change in the effective intersite hopping. This occurs when lattice displacements are large enough to invert the sign of the hopping integral, a regime where the linear approximation of the SSH model begins to break down.
• Competition: In the intermediate regime, BOW and superconducting correlations compete, separated by a metallic region.

C. Magnetic Correlations

• Contrary to findings in square lattice SSH models (where weak antiferromagnetism is often observed), the triangular lattice optical SSH model shows no enhanced magnetic correlations.
• Introducing e-ph coupling ($\lambda$) actually suppresses spin correlations across all momenta and regimes. This suggests that magnetic order is not a dominant feature in this geometry, likely due to the interplay of frustration and longer-range exchange couplings mediated by phonons.

D. Validity of Linear Approximation

• The study monitors the breakdown of the linear approximation (Eq. 1) by tracking sign flips in the effective hopping integrals.
• The superconducting phase at $\langle n \rangle = 1.5$ and the upper portion of the BOW phase occur in regions where the linear approximation begins to fail (1–5% sign flips), suggesting that nonlinear e-ph interactions may be crucial for stabilizing these phases.

4. Significance

• Role of Frustration: The work demonstrates that geometric frustration on a triangular lattice fundamentally alters the competition between ordered phases compared to bipartite lattices. It suppresses charge order and magnetism while allowing for the stabilization of $s$-wave superconductivity at high phonon frequencies.
• Mechanism for Superconductivity: The identification of a superconducting phase driven by large lattice displacements (sign-flipping hopping) provides a new perspective on phonon-mediated pairing, distinct from standard Holstein or small-polaron mechanisms.
• Methodological Benchmark: The successful application of sign-problem-free DQMC to the optical SSH model on a frustrated lattice provides a robust framework for studying complex e-ph systems.
• Future Directions: The results suggest that nonlinear electron-phonon interactions are likely essential for a complete description of the phase boundaries, particularly for the superconducting state, and that doping away from the identified special fillings could reveal further rich physics.

In summary, the paper establishes that the triangular optical SSH model hosts a rich phase diagram where frustration suppresses magnetism, BOW phases dominate at low phonon energies, and $s$-wave superconductivity emerges at high phonon energies via a mechanism involving significant lattice distortions.