Anharmonicity Driven by Vacancy Ordering Unlocks High-performance Thermoelectric Conversion in Defective Chalcopyrites II-III$_2$-VI$_4$

This study reveals that vacancy ordering in II-III$_2$-VI$_4$ defective chalcopyrites induces strong lattice anharmonicity and four-phonon scattering to ultralow thermal conductivity, while anion substitution optimizes electronic transport, collectively enabling high-performance thermoelectric conversion exemplified by CdGa$_2$Te$_4$ with a room-temperature $ZT$ of 0.957.

The Gist

Imagine you are trying to build the ultimate thermoelectric generator. Think of this device as a magical bridge that turns waste heat (like the heat from a car engine or a computer chip) directly into electricity.

To make this bridge work well, you need a material that acts like a one-way street for heat: it should let electricity flow through easily (like a highway for cars) but block heat from flowing (like a brick wall). This is a very difficult balancing act because, in most materials, if electricity flows well, heat usually flows well too.

This paper introduces a new family of materials called Defective Chalcopyrites (specifically a compound called CdGa₂Te₄) that solves this problem using two clever tricks: Vacancy Ordering and Anion Tuning.

Here is the breakdown of how they did it, using simple analogies:

1. The "Ghost" Trick: Vacancy Ordering (Stopping the Heat)

Most crystals are like a perfectly organized dance floor where everyone (atoms) holds hands in a rigid, symmetrical pattern. In this perfect order, heat (which is just atoms vibrating) can travel quickly across the room, like a wave passing through a crowd doing "the wave."

The researchers found that in these special materials, some of the dancers are missing. These missing spots are called vacancies. But here's the magic: the vacancies aren't random; they are ordered. They form a specific, repeating pattern.

• The Analogy: Imagine a dance floor where every third dancer is missing, but they are missing in a perfect, predictable rhythm. This breaks the symmetry of the floor.
• The Result: When heat tries to travel through this "broken" floor, it gets confused. The vibrations (phonons) hit the empty spots and the distorted bonds, causing them to scatter wildly.
• The "Four-Phonon" Effect: Usually, heat particles bounce off each other in pairs (like two billiard balls hitting). But because the lattice is so distorted by these vacancies, the heat particles start colliding in groups of four at once. It's like a chaotic mosh pit where four people crash into each other simultaneously. This "four-phonon scattering" is so effective that it stops heat almost dead in its tracks.
• The Outcome: The material becomes an ultra-insulator for heat. The heat conductivity drops to an incredibly low level (0.19 W·m⁻¹K⁻¹), meaning the heat stays put instead of leaking away.

2. The "Tuning Knob" Trick: Anion Substitution (Boosting the Electricity)

While the vacancies stop the heat, the material still needs to let electricity flow. The researchers realized they could tune the electrical properties by swapping out one specific ingredient: the Anion (the negative part of the molecule, like Sulfur, Selenium, or Tellurium).

• The Analogy: Think of the electrons (electricity) as runners on a track. In some materials, the track is full of hurdles and the runners are heavy, making them slow.
• The Switch: By swapping a smaller, "greedy" atom (like Sulfur) for a larger, "lazy" atom (like Tellurium), the researchers loosened the grip the atoms have on the electrons.
• The Result: This change narrows the "energy gap" (the hurdle height). Now, the electrons can jump the hurdle much easier. The material becomes a highway for electricity, allowing current to flow with very little resistance.

3. The Grand Finale: The Perfect Balance

By combining these two tricks, the researchers created a material that is a thermal brick wall but an electrical superhighway.

• The Star Player: The specific compound CdGa₂Te₄ (Cadmium-Gallium-Tellurium) turned out to be the champion.
• The Score: It achieved a "ZT" score of 0.957 at room temperature. In the world of thermoelectrics, a score near 1.0 is considered excellent and commercially viable. This is much better than many traditional materials.

Why Does This Matter?

Think of this discovery as finding a new type of smart fabric.

• Old fabrics let heat and electricity pass through equally (bad for generators).
• This new fabric is like a thermos flask for heat (keeping the heat where you want it) but a super-conductor for electricity (letting power flow freely).

In summary: The paper shows that by intentionally leaving "ghost" holes in the crystal structure (vacancies) and swapping ingredients to loosen the atomic bonds, we can create materials that turn waste heat into electricity with unprecedented efficiency. This could lead to better batteries, more efficient car engines, and computers that don't overheat.

Technical Summary

1. Problem Statement

Thermoelectric (TE) materials convert heat into electricity, but their efficiency is limited by the strong coupling between electrical and thermal transport. High-performance TE materials require high electrical conductivity ($\sigma$) and Seebeck coefficient ($S$) but low lattice thermal conductivity ($\kappa_L$).

• The Challenge: Traditional strategies to reduce $\kappa_L$ often rely on extrinsic methods like point-defect doping or nanostructuring, which can degrade electrical properties.
• The Opportunity: Defective chalcopyrites (II-III$_2$-VI$_4$) possess intrinsic ordered vacancies that lower crystal symmetry and induce lattice distortion. However, the microscopic mechanisms governing their thermal transport (specifically the role of higher-order phonon scattering) and the optimal electronic structure for carrier transport remain unclear.
• Goal: To systematically investigate the structure-property relationships in II-III$_2$-VI$_4$ compounds to identify design principles for ultralow $\kappa_L$ and high TE figure of merit ($ZT$).

2. Methodology

The authors employed a multi-scale computational approach combining First-Principles Density Functional Theory (DFT) with Machine-Learning Interatomic Potentials (MLIPs):

• Electronic Structure: DFT calculations were performed using the VASP code with the HSE06 hybrid functional to accurately predict band gaps. Transport properties were calculated using the AMSET code, considering acoustic deformation potential (ADP), polar optical phonon (POP), and ionized impurity (IMP) scattering mechanisms.
• Thermal Transport:
  • Force Constants: Third- and fourth-order interatomic force constants (IFCs) were calculated using the finite displacement method. Due to the computational cost of fourth-order IFCs, Moment Tensor Potentials (MTP) and Deep Potential (DP) models were trained on DFT data to generate forces for thousands of displacement configurations.
  • Phonon Scattering: Lattice thermal conductivity ($\kappa_L$) was calculated using the Linearized Wigner Transport Equation (LWTE) framework. This included both three-phonon (3ph) and four-phonon (4ph) scattering processes, as well as coherent contributions ($\kappa_C$).
• Materials Studied: A series of defective chalcopyrites including Zn/Ga/In/Cd combinations with S/Se/Te anions (e.g., CdGa$_2$Te$_4$, ZnGa$_2$S$_4$, ZnIn$_2$Te$_4$).

3. Key Contributions & Findings

A. Mechanism of Ultralow Lattice Thermal Conductivity

The study reveals that ordered intrinsic vacancies are the primary driver of ultralow $\kappa_L$, acting through several synergistic mechanisms:

1. Symmetry Breaking & Anharmonicity: Vacancy ordering lowers the space group symmetry from $I\bar{4}2d$ (chalcopyrite) to $I\bar{4}$, intensifying tetrahedral distortion. This induces strong lattice anharmonicity, characterized by a large anharmonicity index ($\sigma_A$) and strongly negative Gr"uneisen parameters.
2. Metavalent Bonding (MVB): The vacancies promote a metavalent bonding character (intermediate between covalent and metallic). This bonding nature softens low-frequency optical phonons and enhances phonon-phonon coupling.
3. Dominance of Four-Phonon Scattering:
   • The inclusion of 4ph scattering reduces $\kappa_L$ by 50–70% at 300 K compared to calculations considering only 3ph scattering.
   • For CdGa$_2$Te$_4$, 4ph scattering reduces $\kappa_L$ from 0.61 to 0.19 W·m$^{-1}$K$^{-1}$ (a 68% reduction).
   • The weighted phase space for 4ph scattering ($WP_4$) in CdGa$_2$Te$_4$ is two orders of magnitude larger than in ZnGa$_2$S$_4$, indicating a massive number of available scattering channels.
   • Coherent phonon transport ($\kappa_C$) was found to be negligible ($\sim 10^{-6}$ W·m$^{-1}$K$^{-1}$), confirming that incoherent 4ph scattering is the dominant heat-blocking mechanism.

B. Electronic Structure Tuning via Anion Engineering

The study establishes a clear trend for optimizing electrical transport by substituting the VI-site anion (S $\to$ Se $\to$ Te):

• Band Gap Narrowing: Decreasing anion electronegativity (heavier anions) weakens metal-anion hybridization and shifts anion $p$-states upward. This narrows the band gap ($E_g$).
• Enhanced Conductivity: The reduced band gap lowers the excitation energy for carriers, increasing intrinsic carrier concentration and electrical conductivity ($\sigma$).
• Orbital Hybridization: Analysis of Crystal Orbital Hamilton Population (COHP) shows that weaker hybridization in heavier anion systems reduces bonding-antibonding splitting, further facilitating carrier transport.

C. Performance of CdGa$_2$Te$_4$

CdGa$_2$Te$_4$ emerges as the optimal material in this family due to the synergy of the above mechanisms:

• Ultralow $\kappa_L$: 0.19 W·m$^{-1}$K$^{-1}$ at 300 K (driven by heavy atomic mass, strong anharmonicity, and dominant 4ph scattering).
• High Electrical Performance: High carrier mobility and optimized power factor ($PF$) due to the narrowed band gap.
• Figure of Merit ($ZT$): Achieves a room-temperature $ZT$ of 0.957.
  • Without 4ph scattering, $ZT$ would be only 0.416.
  • This value significantly outperforms traditional TE materials like PbSe and AgSbTe$_2$ at room temperature.

4. Significance

• Paradigm Shift in Thermal Transport: The work demonstrates that four-phonon scattering is not a minor correction but a dominant mechanism in defective chalcopyrites, essential for accurately predicting $\kappa_L$.
• Design Principle: It establishes a unified design strategy: Vacancy Ordering (for phonon suppression) combined with Anion Engineering (for carrier optimization).
• Material Discovery: Identifies CdGa$_2$Te$_4$ as a high-performance, room-temperature thermoelectric candidate, validating defective chalcopyrites as a promising platform for next-generation energy conversion technologies.
• Methodological Advance: Successfully integrates MLIPs (MTP/DP) with high-order anharmonic calculations to handle the computational complexity of 4ph scattering in complex crystal structures.

In conclusion, the paper provides a microscopic framework linking vacancy-induced structural distortion to metavalent bonding, higher-order phonon scattering, and electronic band engineering, offering a robust pathway to high-efficiency thermoelectric materials.