Anomalous Thermal Transport Reveals Weak First-Order Melting of Charge Density Waves in 2H-TaSe2

This study reveals that the weak first-order melting of the charge density wave state in 2H-TaSe2 is driven by dislocations and fluctuations, evidenced by anomalous V-shaped thermal conductivity and persistent short-range lattice distortions up to 300 K.

The Gist

The Big Picture: A Mystery in a Crystal

Imagine a crystal made of layers of atoms, like a stack of pancakes. In this specific crystal (called 2H-TaSe₂), the atoms like to organize themselves into a special pattern called a Charge Density Wave (CDW). Think of this pattern like a group of people in a stadium doing "The Wave." When they are organized, the wave moves smoothly across the crowd.

Scientists have known for decades that when you heat this crystal up, the "Wave" eventually breaks down, and the atoms go back to being messy and disorganized. This is called "melting."

Usually, scientists think melting happens in one of two ways:

1. The "Snap": A sudden, sharp break (like ice melting into water).
2. The "Slow Fade": A gradual loosening where the wave just gets weaker and weaker until it's gone.

But this new paper says: "Wait a minute. It's doing something weird."

The Strange Clue: The "V" Shape

The researchers didn't just look at the atoms; they measured how well the crystal conducts heat (thermal conductivity). Imagine heat as a crowd of people trying to walk through a hallway.

• Normal Hallway: As the hallway gets hotter, people get more jittery and bump into each other, making it harder to walk. Heat flow goes down.
• This Crystal's Hallway: The researchers found something bizarre. As they heated the crystal, the heat flow went down, hit a bottom, and then started going back up at high temperatures.

On a graph, this looks like a V-shape. It's like walking down a hill, hitting a valley, and then suddenly starting to walk uphill again without any stairs. This "uphill" part is impossible to explain with normal physics. It meant there was a hidden force helping the heat move.

The Detective Work: Finding the "Ghost" Pattern

To solve the mystery, the team used two powerful tools:

1. Electron Microscopy (SAED): Like taking a super-magnified photo of the atoms.
2. X-ray Diffraction (XRD): Like shining a light through the crystal to see how the atoms are arranged.

What they found:
Even when the crystal was hot enough that the "Wave" should have completely disappeared (the "normal" state), the photos showed that small patches of the Wave were still there.

Imagine a stadium where the big "Wave" has stopped, but if you look closely, you can still see small groups of 5 or 10 people doing the wave locally. These are short-range correlations. They are "ghosts" of the pattern that refuse to die.

The Solution: A "Weak" Melting

The researchers realized the crystal isn't melting like ice. It's undergoing a "Weak First-Order Melting."

Here is the analogy:

• Normal Melting (Ice): One second it's solid, the next it's liquid. Total chaos.
• This Crystal's Melting: It's like a crowd of people at a concert. The main organized dance stops, but instead of everyone running away, they break into small, chaotic dance circles that keep moving around.
  • The "V" Shape Explained:
    • Going Down: As the crystal heats up, the big organized wave breaks into these smaller, chaotic circles. These circles act like speed bumps, slowing down the heat flow.
    • The Bottom: At a certain temperature (around 210 K), there are so many of these chaotic circles that heat flow is at its slowest.
    • Going Up: As it gets even hotter, these small circles start to dissolve completely. Once they are gone, the "speed bumps" disappear, and the heat can flow freely again.

Why Does This Matter?

1. It's a New Kind of Physics: This proves that in some materials, order doesn't just vanish; it lingers in a "zombie" state for a long time, influencing how the material behaves even when it looks disordered.
2. Heat is a Better Detective: The paper shows that measuring heat flow is a super-sensitive way to find these hidden "ghost" patterns. Other tools (like looking at electricity) missed them because the heat carriers (phonons) bump into these hidden patterns, but electricity doesn't.
3. Connection to Superconductors: The authors suggest this might be similar to what happens in high-temperature superconductors (materials that conduct electricity with zero resistance). If we can understand how these "ghost waves" behave in this crystal, we might learn how to make better superconductors.

Summary

The scientists found that a special crystal doesn't just "melt" when heated. Instead, its organized pattern breaks into tiny, lingering fragments that act like speed bumps for heat. By measuring how heat flows, they discovered a hidden "V-shaped" behavior that revealed a new, messy, and fascinating way that matter changes from order to chaos. It's a reminder that sometimes, the things you can't see (hidden fluctuations) are the most important things to measure.

Technical Summary

1. Problem Statement

The nature of phase transitions in low-dimensional quantum materials, particularly Charge Density Wave (CDW) transitions, remains a subject of debate. While textbook models distinguish between continuous second-order transitions and discontinuous first-order transitions, low-dimensional systems often exhibit complex behaviors involving topological defects and strong fluctuations (e.g., the Kosterlitz-Thouless-Halperin-Nelson-Young or KTHNY framework).

In the layered transition metal dichalcogenide 2H-TaSe₂, the CDW transition has been studied for decades, yet its fundamental mechanism is controversial. Key unresolved issues include:

• Hidden Fluctuations: The relevant excitations (charge-neutral lattice distortions) are difficult to probe with conventional electronic or optical methods.
• Contradictory Observations: The material exhibits a large CDW gap (~0.1 eV) relative to its transition temperature ($T_{CDW} \approx 122$ K), short coherence lengths, and metallic optical conductivity below $T_{CDW}$, which challenge standard Fermi surface nesting models.
• Melting Mechanism: It is unclear whether the CDW melting is a continuous topological unbinding (KTHNY), a strict second-order transition (Ising-like), or a weak first-order transition involving latent heat and defect proliferation.

2. Methodology

The authors employed a multi-modal approach combining thermal transport measurements with structural characterization to probe charge-neutral lattice dynamics:

• Thermal Transport (Transient Thermal Grating - TTG):
  • Used a non-contact optical pump-probe technique to measure in-plane thermal diffusivity from 77 K to 350 K.
  • Combined diffusivity with specific heat (Debye model) and density to calculate thermal conductivity ($\kappa$).
  • This method is sensitive to phonon scattering mechanisms, including those caused by charge-neutral fluctuations.
• Structural Characterization:
  • Selected Area Electron Diffraction (SAED): Performed on exfoliated flakes to detect short-range CDW order and superlattice peaks at high temperatures (up to 300 K).
  • X-ray Diffraction (XRD): Conducted at CHESS with cooling and warming cycles to detect thermal hysteresis in the CDW wavevector ($q_{CDW}$).
  • Inelastic X-ray Scattering (IXS): Performed at the Advanced Photon Source to probe lattice dynamics and phonon softening near the CDW wavevector.
• Theoretical Modeling:
  • First-Principles Calculations: Used Temperature Dependent Effective Potentials (TDEP) with VASP to model finite-temperature lattice dynamics.
  • Phenomenological Model: Developed a model linking phonon scattering rates to the spatial gradients of the CDW order parameter ($\psi$), accounting for amplitude fluctuations ($A$) and wavefront distortions ($u$).

3. Key Results

A. Anomalous Thermal Conductivity

The thermal conductivity ($\kappa$) of 2H-TaSe₂ exhibits a striking V-shaped temperature dependence that defies conventional phonon-phonon scattering models:

1. Peak at 122 K: A sharp peak coincides with the incommensurate-to-normal CDW transition, driven by a $\lambda$-shaped specific heat peak (critical fluctuations).
2. Dip at ~210 K: $\kappa$ decreases upon warming, reaching a minimum.
3. High-Temperature Upturn: Above 210 K, $\kappa$ anomalously increases with temperature. In typical crystals, $\kappa$ decreases monotonically due to enhanced Umklapp scattering. This recovery suggests a reduction in a specific scattering channel at high temperatures.
   • Note: The electronic contribution (Wiedemann-Franz law) is negligible (~5 W/m·K), confirming the anomaly is phononic.

B. Persistence of Short-Range Order

• SAED Findings: Distinct CDW superlattice satellite peaks persist up to 300 K, far above the long-range order transition (122 K).
• Coherence Length: The width of these peaks increases linearly with temperature, indicating a decreasing coherence length that follows a scaling behavior ($1/|T-T_{CDW}|$) extending over 150 K above $T_{CDW}$.
• Absence of Hexatic Phase: The peaks fade as lobes rather than forming a diffuse ring, suggesting the CDW wavevector remains locked to crystallographic axes (orientational pinning) even as translational order is lost. This rules out a canonical KTHNY hexatic melting scenario.

C. Thermodynamic Hysteresis

• XRD Findings: Clear thermal hysteresis was observed in the CDW wavevector ($q_{CDW}$) during cooling and warming cycles near the transition.
• Implication: Hysteresis indicates latent energy, a definitive signature of a first-order transition.

D. Phenomenological Model Success

The authors modeled the phonon scattering rate ($\tau_{CDW}^{-1}$) as proportional to the variance of the order parameter gradients:
$$\tau_{CDW}^{-1} \propto \langle |\nabla A|^2 \rangle + q^2 \langle A^2 (\nabla u)^2 \rangle$$

• Mechanism: At intermediate temperatures (122–210 K), growing spatial fluctuations and disorder enhance phonon-CDW scattering, reducing $\kappa$. At higher temperatures (>210 K), the CDW amplitude ($A$) collapses, and the domains become dilute, reducing scattering events and allowing $\kappa$ to recover.
• Fit: The model successfully reproduces the experimental V-shaped curve, confirming that the anomaly arises from the competition between fluctuation-driven scattering and the dilution of local order.

4. Key Contributions

1. Discovery of Weak First-Order Melting: The study provides the first experimental evidence of a weak first-order melting scenario for CDWs in a layered TMDC. This regime is characterized by a subtle discontinuity (latent heat/hysteresis) coexisting with strong, long-range fluctuations and defect proliferation.
2. Thermal Transport as a Probe: Demonstrates that longitudinal thermal conductivity is a highly sensitive, non-destructive probe for charge-neutral lattice fluctuations and hidden order, capable of detecting phenomena invisible to electronic probes.
3. Resolution of the "Melting" Debate: The results argue against a pure KTHNY (continuous) or strict Ising (second-order) model. Instead, they propose a unified picture where dislocation-driven melting and fluctuation-driven persistence of local order coexist.
4. Implications for Correlated Materials: The finding that short-range CDW correlations persist well into the "normal" phase suggests that conventional phase diagrams may underestimate the temperature range where CDW physics influences material properties. This has direct relevance to high-$T_c$ superconductors, where fluctuating charge orders persist in the pseudogap regime.

5. Significance

This work fundamentally alters the understanding of phase transitions in low-dimensional quantum materials. By identifying a weak first-order melting mechanism driven by the interplay of defects and fluctuations, it bridges the gap between theoretical frameworks of topological melting and traditional first-order transitions.

Furthermore, the study establishes thermal transport as a critical diagnostic tool for "hidden" degrees of freedom in quantum materials. The observation that local CDW order survives up to 300 K implies that electron-phonon coupling and electronic structure modifications persist far beyond the nominal transition temperature, offering new avenues for understanding and engineering correlated electron systems like cuprate superconductors.