🔬 physics

Vibrational spectroscopy identifies the bond asymmetry of hexagonal diamond

By combining first-principles lattice dynamics with Raman spectroscopy, this study resolves the structural controversy of bulk hexagonal diamond by identifying a small positive bond asymmetry where interlayer bonds are longer than intralayer bonds, a finding that contradicts previous refinements and suggests the 1,529 cm⁻¹ Raman feature does not originate from ideal 2H diamond.

Original authors: Li Zhu

Published 2026-08-04

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Li Zhu

Original paper licensed under CC BY 4.0 (https://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

The Diamond Detective Story: Why Shape Matters More Than You Think

Imagine a world built entirely of tiny, super-strong Lego bricks. In this world, the most famous structure is the "cubic diamond," a perfect, rigid cube where every brick is connected to its neighbors in a staggered, zig-zag pattern. This is the hardest natural material on Earth, used in everything from jewelry to drill bits. But scientists have been hunting for a different, hexagonal version of this structure for decades, often called "lonsdaleite" or "hexagonal diamond." Think of it as the hexagonal cousin of the standard diamond. While the standard diamond is staggered like a staircase, this hexagonal cousin has layers that line up directly on top of each other, like a stack of pancakes.

Why does this matter? Because the way these layers stack changes how the "glue" (the chemical bonds) holds them together. In the standard diamond, all the bonds are the same length and strength. But in the hexagonal version, the bonds connecting the layers might be different from the bonds inside the layers. This tiny difference—measured in thousandths of an angstrom (a unit so small it's hard to imagine)—could make the material even harder or stronger than regular diamond. For years, scientists have argued over whether the "pancake glue" is shorter and tighter, or longer and looser than the rest. It's a bit like trying to figure out if a specific link in a chain is slightly stretched or slightly squished, just by looking at the whole chain. Getting this right is crucial because if we can control this structure, we might be able to create the ultimate super-material.

The Great Bond Length Debate

For a long time, the scientific community was stuck in a "he said, she said" situation regarding hexagonal diamond. Two recent studies, both claiming to have synthesized large, pure samples of this material, looked at the same thing but saw completely opposite things. One team said the bond connecting the layers was significantly shorter than the bonds inside the layers. The other team said it was much longer. They were off by a massive margin in the world of atoms—about 238 millionths of a nanometer (238 mÅ). It was as if one group said a bridge was 10 feet long and the other said it was 10 feet shorter than that, and both were sure they were right.

Enter Li Zhu, a researcher at Rutgers University, who decided to settle the argument not by building more diamonds, but by using a super-powerful computer simulation to act as a detective. Zhu didn't just guess; they used "first-principles lattice dynamics," which is a fancy way of saying they calculated how the atoms should vibrate based on the laws of physics, without needing to tweak the numbers to fit a specific result.

The Vibration Test: Listening to the Atoms

Here is the clever trick Zhu used. Atoms in a solid aren't static; they are constantly jiggling and vibrating. These vibrations create a unique "fingerprint" called a vibrational spectrum, which can be read using a technique called Raman spectroscopy. Think of it like tuning a guitar: if you tighten a string (shorten a bond), the pitch goes up. If you loosen it (lengthen a bond), the pitch goes down.

Zhu's simulations revealed a very specific rule: the "pitch" of a particular vibration (called the A1g mode) is directly tied to the length of the bond connecting the layers. The computer showed that this pitch changes by about -2,100 cm⁻¹ for every angstrom the bond length changes. This means the vibration acts like a super-precise ruler.

When Zhu applied this "ruler" to the data from the two conflicting studies, the results were clear. Neither of the recently proposed structures matched the actual sounds (spectra) the materials were making.

The structure claiming the bond was shorter predicted vibrations that were way too high-pitched compared to what was measured.
The structure claiming the bond was longer predicted vibrations that were way too low-pitched.

In fact, the only structure that matched the "song" of the atoms was one where the interlayer bond was slightly longer than the intralayer bonds. Specifically, the computer simulations showed the bond connecting the layers is about 24 mÅ longer than the bonds inside the layers. This result aligns with a much older study from 2003, which had suggested the same thing but was largely ignored in the recent debate.

Why the "Short Bond" Idea Fails

The paper doesn't just say the "short bond" idea is wrong; it explains why it's physically impossible for the material to be in that state under normal conditions. The reason the layers want to be slightly stretched apart is due to a phenomenon called "eclipsed conformation." Imagine two people standing on top of each other, both with their arms raised. If they stand directly on top of each other (eclipsed), their arms bump into each other, creating a little bit of repulsion. This pushes them slightly apart, making the bond longer. In the hexagonal diamond, the atoms are stacked in this "eclipsed" way, so the bond naturally wants to be longer.

To force that bond to be shorter (as the other study claimed), you would need to squeeze the material with tens of gigapascals of pressure—like the weight of a mountain on a postage stamp. But the study's own measurements of the material's size (diffraction data) prove it wasn't under that kind of pressure. Therefore, the "short bond" structure is a mathematical possibility but a physical impossibility for the sample they had.

The Mystery of the Extra Note

There was one final puzzle. The study that claimed the bond was short also reported a strange, high-pitched "note" in their data at 1,529 cm⁻¹ that didn't fit the hexagonal diamond pattern at all. Zhu's investigation suggests this isn't a feature of the diamond itself. Instead, it's likely caused by tiny bits of "graphite" (a softer, layered form of carbon) that got stretched or disordered during the creation process. Graphite has a known vibration that shifts to exactly 1,529 cm⁻¹ when it is under tension. Since this high-pitched note appeared in some samples but not in the purest ones, and since it doesn't fit the rules of a perfect hexagonal diamond, the paper concludes it's an "imposter" signal from a different material, not a new type of diamond bond.

The Verdict

By combining computer simulations with real-world sound measurements, this paper effectively solves the mystery. The hexagonal diamond does exist, and its structure is confirmed to have a slightly longer bond between its layers than within them (by about 24 mÅ). The two recent studies that claimed the opposite were likely misled by difficult-to-measure experimental data. The paper also highlights that the "eclipsed" stacking of the layers is the key reason for this stretch, a rule that seems to apply to other similar carbon structures too. While the debate was fierce, the "song" of the atoms has finally been heard clearly, pointing back to the older, correct understanding of how this super-hard material is built.

Problem Statement
Hexagonal diamond (HD, or lonsdaleite) has been synthesized in bulk form through independent routes, yet its fundamental crystal structure remains contested. Specifically, the sign and magnitude of the bond asymmetry between the three intralayer bonds ( $O_A$ ) and the single interlayer bond ( $O_B$ ) are unresolved. Two recent structure refinements on bulk samples provide mutually exclusive answers: one suggests a significantly shorter interlayer bond ( $O_B - O_A \approx -55$ mÅ), while the other suggests a significantly longer one ( $O_B - O_A \approx +183$ mÅ). Both recent determinations also deviate from an earlier 2003 refinement. Furthermore, the literature contains conflicting interpretations of Raman spectra, including a high-frequency feature at 1,529 cm $^{-1}$ that has been difficult to assign to an ideal 2H diamond structure.

Methodology
The study employs first-principles lattice dynamics based on Density Functional Theory (DFT) to adjudicate the structural controversy. The author utilizes two functionals: a dispersion-corrected meta-GGA (r2SCAN+rVV10) and the PBE functional.

Structural Relaxation: The hexagonal diamond structure (space group $P6_3/mmc$ ) was fully relaxed to determine the equilibrium bond lengths and internal coordinates.
Polytype and Strain Analysis: The study constructs various polytypes (2H, 4H, 6H, and twinned supercells) and applies uniaxial and biaxial strains to map the relationship between stacking sequence, stress, and bond asymmetry.
Lattice Dynamics and Raman Simulation: Zone-center phonon frequencies and Raman activities were calculated along a path of varying internal coordinates. The study establishes a "sensitivity map" linking the frequency of the $A_{1g}$ mode (interlayer bond stretch) to the interlayer bond length ( $O_B$ ).
Spectral Inversion: The measured Raman spectra from two distinct bulk HD samples (one phase-pure, one twinned) and an impact diamond sample were inverted using the DFT-calibrated sensitivity map to infer the bond asymmetry.
Cross-Validation: The inferred structures were tested against the measured X-ray diffraction (XRD) lattice constants and compared against the published Rietveld refinements.

Key Contributions and Results

Equilibrium Structure: DFT relaxation reveals that bulk hexagonal diamond possesses a small, robust positive bond asymmetry where the interlayer bond is longer than the intralayer bonds by approximately 24 mÅ ( $O_B - O_A \approx +24$ mÅ). This result aligns with the sign and magnitude of the 2003 refinement but contradicts the two recent refinements.
Mechanism: The asymmetry is identified as a local effect of the "eclipsed" conformation of the interlayer bond in hexagonal stacking, as opposed to the "staggered" conformation in cubic diamond. This effect scales with polytype hexagonality and is erased at cubic stacking faults.
Vibrational Gauge: The bright zone-center $A_{1g}$ mode is established as a direct gauge of the interlayer bond, with a frequency sensitivity of approximately $-2,100$ cm $^{-1}$ Å $^{-1}$ .
Spectral Inversion: Inverting the Raman spectrum of the phase-pure sample yields a bond asymmetry of $+24 \pm 3$ mÅ (95% interval). Independent determinations on impact diamonds and diffraction data from 2003 yield consistent positive values ( $+33 \pm 8$ mÅ and $+60 \pm 45$ mÅ, respectively).
Refutation of Recent Models:
- The two recent Rietveld refinements (Yang et al. and Lai et al.) are shown to be inconsistent with the measured Raman spectra of their own samples. Neither refined coordinate reproduces the full pattern of measured modes within the ideal-crystal harmonic model.
- The structure proposed by Yang et al. to explain the 1,529 cm $^{-1}$ band (requiring a short interlayer bond of $\approx 1.44$ Å) is mathematically possible but physically unviable. A structure matching all three of Yang's bands requires tens of gigapascals of confining stress and lattice constants that are excluded by the sample's own ambient XRD data.
Origin of the 1,529 cm $^{-1}$ Feature: The paper concludes that the 1,529 cm $^{-1}$ feature cannot be assigned to a homogeneous ideal 2H diamond. Instead, it is likely an extrinsic feature arising from minority sp $^2$ carbon (strained or disordered) at grain boundaries or twins, or potentially locally compressed regions that are mechanically inconsistent with the bulk sample's diffraction data.

Significance
The paper claims to resolve the structural controversy of hexagonal diamond by demonstrating that vibrational spectroscopy, when calibrated by first-principles lattice dynamics, provides a sharper discrimination of structure than the available diffraction data. The study asserts that:

The intrinsic structure of hexagonal diamond has a long interlayer bond ( $O_B > O_A$ ), consistent with the 2003 refinement and theoretical predictions of eclipsed-bond repulsion.
The two recent diffraction refinements are inconsistent with each other and with the vibrational spectra of the samples they describe.
The $A_{1g}$ mode serves as a reliable gauge for the interlayer bond, allowing for the determination of bond asymmetry directly from Raman spectra.
The high-frequency Raman feature (1,529 cm $^{-1}$ ) observed in some samples is not a fundamental mode of ideal hexagonal diamond but likely stems from extrinsic factors such as strained sp $^2$ carbon.

The work effectively removes the bond asymmetry question from the domain where theory must arbitrate between conflicting experiments, showing instead that the calculation agrees with the oldest measurement and all measured spectra, while disagreeing with the two most recent refined coordinates.