🔬 materials science

Anisotropic Phonon Heat Flow and Thermoelectric Response in Tetragonal GeS $_2$ and GeSe $_2$

This study employs first-principles calculations to demonstrate that tetragonal GeS $_2$ and GeSe $_2$ exhibit pronounced anisotropic thermoelectric properties, characterized by significantly suppressed cross-plane lattice thermal conductivity and moderate figure-of-merit ( $zT$ ) values, particularly for n-type GeS $_2$ at high temperatures.

Original authors: Neeraj Kulhari, Krishna Swaroop Sharma, K. C. Bhamu

Published 2026-08-05

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Neeraj Kulhari, Krishna Swaroop Sharma, K. C. Bhamu

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ 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 the world is constantly trying to turn waste heat into electricity, like a magical engine that runs on the warmth of a summer day or a car's exhaust. This is the dream of thermoelectric materials. To make this work, scientists look for special crystals that are good at letting electricity flow through them but terrible at letting heat flow. Think of it like a bouncer at a club: you want the electricity (the cool kids) to get in, but you want to stop the heat (the rowdy crowd) from following them. The measure of how good a material is at this job is called zT. A high zT means you get a lot of electricity for the heat you put in.

But here's the tricky part: usually, if a material lets electricity flow easily, it also lets heat flow easily. It's like a wide-open door that lets everyone in. To fix this, scientists look for materials where the rules are different depending on which way you look. Some materials are like a maze: if you try to walk straight through them, you get stuck (heat stops), but if you walk sideways, you can zip right through (electricity flows). This paper dives into two specific materials, Germanium Sulfide (GeS2) and Germanium Selenide (GeSe2), to see if they are the perfect "mazes" for turning heat into power.

The Crystal Maze: A Tale of Two Geometries

In this study, the researchers acted like digital architects, building and testing two imaginary crystal cities: one made of Germanium and Sulfur (GeS2), and the other made of Germanium and Selenium (GeSe2). They didn't just build them; they simulated how atoms vibrate and how electrons zoom around inside them using powerful computer models.

These crystals have a very specific shape called "tetragonal." Imagine a stack of flat, square pancakes. The atoms are arranged in layers, like sheets of paper in a notebook. The researchers wanted to know: if you try to push heat through the layers (sideways) versus pushing it through the stack (up and down), does it behave differently?

The Heat Traffic Jam

The results were surprisingly dramatic. The computer simulations showed that these materials are incredibly picky about which direction heat can travel.

For the Sulfur-based material (GeS2), the heat flows like a highway when moving sideways across the layers. At room temperature (300 K), the heat conductivity is a whopping 26.86 W m⁻¹ K⁻¹. But if you try to push that same heat up through the stack (the "cross-plane" direction), it hits a massive traffic jam. The heat conductivity drops to a tiny 1.19 W m⁻¹ K⁻¹. That's a difference of more than 20 times!

The Selenium-based material (GeSe2) acts similarly, though the "highway" is a bit slower. Sideways, it conducts 18.74 W m⁻¹ K⁻¹, but up and down, it drops to 1.52 W m⁻¹ K⁻¹.

Why does this happen? The researchers found that the atoms in these crystals are arranged in a way that makes them very "soft" when you try to squish them between the layers. It's like the layers are held together by weak springs. When heat tries to travel up and down, the atoms wobble and scatter the heat energy, stopping it in its tracks. However, moving sideways is like skating on ice; the atoms are tightly locked together, letting heat zip through easily.

The Electron Dance

While the heat was getting stuck, the researchers checked how electricity behaved. They found that these materials are semiconductors, meaning they can conduct electricity if you add a few extra electrons (n-type) or remove some (p-type).

However, the electricity didn't get the same "traffic jam" treatment as the heat. The electrons could still move, but the material wasn't a super-conductor. The researchers calculated that for the Sulfur version (GeS2), if you tune the number of electrons just right and heat it up to 800 K, you can get a performance score (zT) of 0.257 when looking at the cross-plane direction.

For the Selenium version (GeSe2), the score was lower, around 0.066 for p-type transport at the same temperature.

What This Means (and What It Doesn't)

The big takeaway is that these materials are "anisotropic," which is a fancy word for "directional." They are excellent at blocking heat in one direction while letting electricity do its thing. This is exactly what you want for a thermoelectric device.

However, the paper is careful to point out that while this is a great start, these materials aren't quite the "holy grail" of energy yet. A zT of 0.257 is considered "moderate." To be truly useful in real-world engines or power plants, scientists usually look for numbers closer to 1.0 or higher. The researchers also noted that their results are based on computer simulations. They haven't built a physical machine with these crystals yet to prove it works in the real world.

The study also ruled out the idea that the materials are weak or unstable. Even though the layers are soft and stop heat, the atoms are actually held together by strong chemical bonds (a mix of covalent and ionic, like a strong handshake with a little bit of a magnet). This means the material is sturdy enough to handle the heat without falling apart.

The Bottom Line

In simple terms, this paper discovered that tetragonal GeS2 and GeSe2 are like one-way streets for heat. They are fantastic at stopping heat from flowing up and down their layers, which is a huge advantage for turning waste heat into electricity. While they aren't the perfect solution just yet, they show a very promising path forward. If scientists can figure out how to boost the electricity flow even more while keeping that heat-blocking superpower, these materials could become key players in our future clean-energy toolkit.

Technical Summary: Anisotropic Phonon Heat Flow and Thermoelectric Response in Tetragonal GeS2 and GeSe2

Problem Statement
Thermoelectric materials, which convert heat directly into electricity, require a high dimensionless figure of merit ($zT$), defined by the ratio of the power factor ( $\sigma S^2$ ) to the total thermal conductivity ( $\kappa_e + \kappa_l$ ). A central challenge in thermoelectric design is the mutual coupling of these parameters: increasing electrical conductivity ( $\sigma$ ) often increases electronic thermal conductivity ( $\kappa_e$ ). While reducing lattice thermal conductivity ( $\kappa_l$ ) through nanostructuring or alloying is a common strategy, there is a need to explore non-toxic, earth-abundant compounds with intrinsically low $\kappa_l$ and favorable electronic properties. Specifically, the electronic, optical, and phonon-transport properties of bulk tetragonal germanium dichalcogenides (GeS $_2$ and GeSe $_2$ ) remain underexplored from a thermoelectric perspective, despite their structural anisotropy and potential for decoupling electronic and phonon transport.

Methodology
The authors employed a comprehensive first-principles framework to investigate the structural, electronic, and thermal properties of tetragonal GeS $_2$ and GeSe $_2$ (space group $P42/nmc$).

Electronic Structure: Calculations were performed using Density Functional Theory (DFT) with the PBE functional and Grimme D2 dispersion corrections for structural relaxation and elastic properties. To obtain accurate band gaps, screened HSE03 hybrid functional calculations were conducted. Electronic transport coefficients were calculated using AMSET (Ab initio Boltzmann Transport), incorporating state-dependent scattering mechanisms including acoustic deformation potentials, polar optical phonons, and ionized impurities.
Lattice Dynamics and Thermal Transport: Harmonic and anharmonic lattice dynamics were investigated using Density Functional Perturbation Theory (DFPT) and finite-displacement methods. Lattice thermal conductivity ( $\kappa_l$ ) was computed using the ShengBTE code, solving the phonon Boltzmann transport equation within the Relaxation Time Approximation (RTA).
Bonding Analysis: The microscopic origin of transport behavior was analyzed using LOBSTER to evaluate Integrated Crystal Orbital Hamilton Population (ICOHP) and Crystal Orbital Overlap Population (ICOOP) descriptors, characterizing the covalent-ionic nature of Ge–X bonds.
Thermoelectric Evaluation: The figure of merit ($zT$) was calculated by combining the AMSET electronic transport tensors with the ShengBTE RTA lattice thermal conductivity tensors, ensuring component-wise consistency (pairing in-plane electronic coefficients with in-plane $\kappa_l$ and cross-plane coefficients with cross-plane $\kappa_l$ ).

Key Contributions and Results

Structural and Mechanical Stability: Both tetragonal GeS $_2$ and GeSe $_2$ were confirmed to be mechanically and dynamically stable. The structures exhibit strong anisotropy, characterized by large in-plane elastic constants ( $C_{11}$ ) and significantly smaller out-of-plane constants ( $C_{33}$ ). Notably, both materials display soft shear deformation channels (low $C_{66}$ ), which is favorable for reducing lattice thermal conductivity.
Electronic Properties: HSE03/Wannier interpolation yielded semiconductor band gaps of 2.48 eV for GeS $_2$ and 1.23 eV for GeSe $_2$ . LOBSTER analysis revealed mixed covalent–ionic bonding, with Ge–S bonds showing a larger stabilizing ICOHP magnitude (−5.27 eV/bond) compared to Ge–Se bonds (−4.74 eV/bond).
Anisotropic Lattice Thermal Conductivity: The study revealed pronounced anisotropy in heat transport. At 300 K, the in-plane lattice thermal conductivity ( $\kappa_{ab}$ ) is 26.86 W m $^{-1}$ K $^{-1}$ for GeS $_2$ and 18.74 W m $^{-1}$ K $^{-1}$ for GeSe $_2$ . In contrast, the cross-plane conductivity ( $\kappa_c$ ) is drastically lower, at 1.19 W m $^{-1}$ K $^{-1}$ for GeS $_2$ and 1.52 W m $^{-1}$ K $^{-1}$ for GeSe $_2$ . This suppression of cross-plane transport is attributed to the layered tetrahedral framework, reduced out-of-plane phonon group velocities, and moderate anharmonicity. Substitution of S with Se lowers the upper phonon frequency from ~13.6 THz to ~10.3 THz.
Thermoelectric Performance: Combining the electronic and lattice transport data, the calculated $zT$ values are moderate. The highest performance is observed for cross-plane n-type GeS $_2$ at 800 K, yielding a $zT_c$ of 0.257 at a carrier concentration of $1 \times 10^{19}$ cm $^{-3}$ . For GeSe $_2$ , the PBE-AMSET estimate gives a $zT_c$ of 0.066 for p-type cross-plane transport at 800 K. The authors note that the GeSe $_2$ results should be treated as qualitative due to the near-zero band gap in the PBE description, which may require gap correction for quantitative accuracy.

Significance and Claims
The paper claims that tetragonal GeX $_2$ compounds represent a class of strongly anisotropic thermoelectric materials where the cross-plane response benefits significantly from suppressed lattice heat transport. The study demonstrates that while low $\kappa_l$ is a necessary condition, it is not sufficient for high $zT$; the performance is ultimately controlled by the simultaneous balance of the power factor, mobility degradation at high carrier concentrations, and electronic heat transport. The authors conclude that pristine tetragonal GeS $_2$ and GeSe $_2$ exhibit moderate calculated $zT$ values, with the most promising transport channel along the $c$ -axis. They suggest that future improvements could be achieved through gap-corrected scattering calculations, defect engineering, alloying, and strain control to further leverage the favorable cross-plane thermal response.