Thermal Drude weight in an integrable chiral clock model

✨

This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you have a long line of people holding hands, passing a secret message down the line. In the world of quantum physics, this "message" is energy, and the "people" are tiny particles. Usually, when you try to pass a message down a crowded line, it gets jumbled, delayed, or lost. This is like heat conduction in a normal material: energy diffuses slowly, like a drop of ink spreading in water.

However, in some special, perfectly ordered lines (called integrable systems), the message can zip down the line without getting lost. This is called ballistic transport. The speed and efficiency of this "perfect delivery" are measured by something physicists call the Drude weight. Think of the Drude weight as a score: a high score means the energy moves like a bullet (ballistic), while a zero score means it moves like a slow leak (diffusive).

This paper investigates a specific, exotic quantum line called the $Z_3$ Chiral Clock Model. Here's a breakdown of what the researchers did and found, using simple analogies:

1. The Setup: A Twisted Clock

Most quantum models are like a standard clock where the hands move in a predictable circle. This model is a "chiral" clock, meaning it has a built-in twist or "handedness." Imagine a clock where the numbers are arranged in a spiral, and moving clockwise is slightly different energetically than moving counter-clockwise.

The researchers focused on a specific "sweet spot" in the model's settings where the system is integrable. In our analogy, this is like tuning the clock so perfectly that the gears never jam, and the hands move in a mathematically perfect, predictable rhythm.

2. The Mystery: Is the Energy Current "Conserved"?

In a famous model called the XXZ chain (a standard benchmark in physics), the energy current is a "conserved charge." Think of this like a VIP pass that never expires; once you have it, you keep it forever, and the energy flows perfectly.

The researchers asked: Does this VIP pass exist in our twisted Chiral Clock model?

The Answer: No, not exactly. The energy current isn't a perfect VIP pass. It's more like a VIP pass that gets stamped with a "partial validity" every time it moves.
The Discovery: Even though the current isn't perfectly conserved, it has a strong "friendship" (overlap) with a specific, hidden rule of the system called $Q(2)$ . This rule acts like a shadow of the VIP pass. Because of this connection, the energy still flows very efficiently, resulting in a finite Drude weight (a high score) even at warm temperatures.

3. The Method: The "Ancilla Disentangler" Trick

To study this, the researchers used a powerful computer simulation technique called tDMRG. Imagine trying to simulate a crowded room where everyone is holding hands. As time goes on, the "entanglement" (the complexity of who is holding whose hand) grows so fast that your computer runs out of memory.

To fix this, they used a trick called the Ancilla Disentangler.

The Analogy: Imagine the simulation is a messy party. The "ancilla" is a group of invisible helpers. The "disentangler" is a bouncer who gently tells the helpers, "You don't need to hold hands with everyone; just stand near the person you're talking to."
The Result: This keeps the party organized. The researchers found this trick works amazingly well in the "perfect" (integrable) clock model, keeping the simulation running for a long time. However, if they broke the perfect rhythm (made the model non-integrable), the bouncer couldn't keep up, and the simulation got messy again.

4. The Findings: What Happens at Different Temperatures?

The team calculated how well energy flows at different temperatures:

High Temperatures: As the system gets hotter, the Drude weight drops, but it follows a predictable mathematical curve. It's like the VIP pass gets slightly less effective as the crowd gets rowdier, but the system still remembers the rules.
Low Temperatures:
- In the gapless (critical) phase, the flow drops linearly as it gets colder.
- In the gapped (ordered) phase, the flow drops exponentially (very fast) as it gets colder, almost stopping.
The "Regular" Part: While the main flow is ballistic (fast), there is also a tiny bit of "leakage" (diffusive transport). This is because the energy current isn't perfectly aligned with the conserved rule $Q(2)$ . It's like a runner who mostly sprints but occasionally stumbles.

5. The Big Picture

The paper confirms that integrability (perfect order) leads to anomalous transport (super-fast energy flow), even in this twisted, chiral model.

They proved that the "Drude weight" is real and finite.
They showed that a specific mathematical rule ( $Q(2)$ ) explains almost all of this efficiency.
They demonstrated that their computer simulation tricks (the disentangler) are great for studying these perfect systems but struggle when the system gets chaotic.

In summary: The researchers took a complex, twisted quantum clock, found a way to make it run perfectly, and discovered that even though the energy current isn't perfectly "conserved," it still zips through the system with incredible efficiency because it's tightly linked to a hidden, conserved rule of the universe. They also showed that their new computer tricks are the best way to watch this happen, as long as the clock stays perfectly tuned.

1. Problem Statement

The paper investigates the anomalous transport properties of one-dimensional quantum integrable systems, specifically focusing on the thermal conductivity and the Drude weight (the zero-frequency delta-function peak in conductivity).

Context: While the spin-1/2 XXZ chain is a well-studied model where thermal current is a conserved quantity (leading to ballistic transport), less is known about other integrable models, particularly those with chiral symmetry and time-reversal invariance.
Specific Gap: The authors aim to characterize the thermal Drude weight in the $Z_3$ chiral clock model along its integrable line. A key question is whether the thermal current operator in this model is a conserved charge (as in the XXZ model) or if it has a finite overlap with conserved charges derived from the transfer matrix, leading to a finite Drude weight.
Computational Challenge: Calculating finite-temperature dynamical properties in large systems requires handling entanglement growth during time evolution. The paper also seeks to evaluate the effectiveness of ancilla disentangler techniques in reducing computational costs for both integrable and non-integrable regimes.

2. Methodology

The authors employ a combination of analytical theory and advanced numerical simulations:

Model: The study focuses on the 1D $Z_3$ chiral clock model with Hamiltonian $H = -\sum [\bar{\alpha}_1\tau_j + \bar{\alpha}_2\tau^\dagger_j + \alpha_1\sigma_j\sigma^\dagger_{j+1} + \alpha_2\sigma^\dagger_j\sigma_{j+1}]$ . They restrict the study to the time-reversal invariant regime ( $\phi=0$ ) along the integrable line defined by $f \cos(3\phi) = J \cos(3\theta)$ .
Numerical Technique:
- Finite-Temperature tDMRG: They use Time-Dependent Density Matrix Renormalization Group (tDMRG) based on Matrix Product States (MPS).
- Ancilla Purification: The system is purified by entangling physical sites with ancilla sites to represent the thermal density matrix $\rho = e^{-\beta H}/Z$ .
- Ancilla Disentangler: To mitigate the exponential growth of bond dimension during real-time evolution, they apply a specific unitary disentangler on the ancilla subsystem. They empirically determine the optimal disentangler parameters to localize entanglement growth around the quench site.
Analytical Framework:
- Mazur Bound: They utilize Mazur's inequality to establish a lower bound for the Drude weight based on the overlap between the thermal current operator $I$ and local conserved charges $Q^{(k)}$ .
- Conserved Charges: They derive local conserved charges $Q^{(k)}$ from the logarithmic expansion of the transfer matrix of the classical chiral clock model. Specifically, they construct $Q^{(2)}$ and $Q^{(3)}$ .
- Sum Rules: They validate numerical results using the thermal conductivity sum rule: $\int_0^\infty \text{Re}\,\kappa(\Omega) d\Omega = \frac{\pi}{2\hbar T L} \langle \Theta \rangle$ .

3. Key Contributions

Identification of Conserved Charge Overlap: The authors demonstrate that unlike the XXZ model (where thermal current is conserved), the thermal current in the $Z_3$ chiral clock model is not a conserved quantity. However, it has a finite overlap specifically with the local conserved charge $Q^{(2)}$ derived from the transfer matrix, while having zero overlap with $Q^{(3)}$ .
Saturation of the Mazur Bound: They show that the Mazur bound calculated using only $Q^{(2)}$ saturates the numerically calculated Drude weight. This proves that $Q^{(2)}$ is the sole conserved charge responsible for the ballistic component of thermal transport in this model.
Asymptotic Expression: Leveraging the saturation of the Mazur bound, they derive an asymptotically exact expression for the Drude weight at high temperatures ( $T \to \infty$ ).
Efficiency of Disentangler: They provide a detailed analysis of the ancilla disentangler, showing it significantly improves computational efficiency in the integrable regime at high temperatures but offers diminishing returns in non-integrable regimes and at low temperatures.

4. Key Results

Finite Drude Weight: The thermal Drude weight $D_{th}$ is finite at all non-zero temperatures for the integrable points, confirming anomalous (ballistic) transport.
Temperature Dependence:
- High Temperature: $D_{th}$ decays as a power law $\sim 1/T^2$ . The derived asymptotic formula is $D_{th}(T=\infty) = \frac{1}{2T^2} \frac{24f^2}{5-f^2}$ .
- Low Temperature:
  - In the gapped phase ( $f \neq 1$ ), $D_{th}$ decays exponentially as $e^{-\Delta/T}$ , where $\Delta$ is the spectral gap.
  - At the critical point (gapless, $f=1$ ), $D_{th}$ scales linearly with $T$ ( $D_{th} \propto T$ ).
Regular Part of Conductivity: The regular part of the thermal conductivity, $\kappa_{reg}(\Omega)$ , is non-zero. It exhibits a broad peak at finite frequencies (related to domain wall energies) and a small finite value at $\Omega \to 0$ . This indicates a coexistence of ballistic and diffusive components in energy transport, unlike the XXZ model where $\kappa_{reg}=0$ .
Numerical Validation: The numerical results for the current-current correlator match the Mazur bound estimates perfectly. The sum rule for thermal conductivity is satisfied, validating the accuracy of the tDMRG calculations.
Disentangler Performance:
- Integrable/High-T: The disentangler reduces the maximum bond dimension growth by up to 2.5x, allowing access to longer time scales.
- Non-Integrable/Low-T: The improvement is marginal. In non-integrable regimes, the transient decay times required to observe diffusive behavior are too long for the disentangler to bridge the gap.

5. Significance

Universality of Integrability: The results reinforce the connection between integrability and anomalous transport across different model classes (XXZ vs. Chiral Clock). Even when the current is not strictly conserved, the existence of a single overlapping conserved charge ( $Q^{(2)}$ ) is sufficient to sustain a finite Drude weight.
Transport Mechanism: The paper clarifies that in chiral clock models, energy transport is a mix of ballistic (due to $Q^{(2)}$ ) and diffusive (due to the non-conserved component of the current), providing a more nuanced view than the purely ballistic transport seen in the XXZ thermal channel.
Methodological Insight: The study offers critical benchmarks for tDMRG techniques, specifically highlighting the limitations and optimal usage of ancilla disentanglers. It suggests that while disentanglers are powerful for integrable high-temperature dynamics, they may not be sufficient to solve the long-time dynamics of non-integrable systems without further algorithmic advancements.
Future Directions: The work opens avenues for studying time-reversal breaking regimes ( $\phi \neq 0$ ) and superintegrable points, as well as the effects of integrability-breaking perturbations on transport.