Anisotropic magnetoelastic coupling in the honeycomb magnet Na$_3$Co$_2$SbO$_6$

This study combines magnetization and dilatometry measurements with *ab initio* calculations to map the field-temperature phase diagram of the honeycomb magnet Na$_3$Co$_2$SbO$_6$, revealing strongly anisotropic magnetoelastic coupling driven by Co--O--Co bond angle variations and demonstrating that its field-induced transitions are first-order without evidence of a quantum spin liquid state.

The Gist

Here is an explanation of the research paper, translated into everyday language with some creative analogies.

The Big Picture: Hunting for a "Quantum Ghost"

Imagine you are a detective trying to find a ghost. In the world of physics, this "ghost" is called a Quantum Spin Liquid (QSL).

Normally, when you cool down a magnet, the tiny atomic magnets (spins) inside it line up in an orderly pattern, like soldiers marching in formation. This is called "magnetic order." But a Quantum Spin Liquid is a weird state of matter where the atoms never line up, even when it's super cold. They keep jiggling and dancing in a chaotic, fluid-like state forever. Scientists think this state might hold the key to building super-fast quantum computers.

For years, scientists have been looking for this ghost in a special type of material called a "honeycomb magnet" (named because the atoms are arranged in a hexagonal pattern, like a beehive). One of the most promising suspects is a material called Na₃Co₂SbO₆ (let's call it NCSO for short).

This paper is the report from a team of detectives who went to NCSO with very sensitive tools to see if they could catch this ghost.

The Investigation: Stretching and Squeezing the Material

To find out what's happening inside NCSO, the scientists didn't just look at it; they played with it. They used two main tricks:

1. The Magnet Test: They applied strong magnetic fields from different directions (like pushing a box from the left or the top) to see how the atoms reacted.
2. The Stretch Test (Dilatometry): This is the star of the show. They measured the material's length with extreme precision as they changed the temperature and magnetic field. Think of it like measuring a rubber band to see if it stretches or shrinks when you pull on it.

The Analogy: Imagine the atoms in NCSO are like a group of people holding hands in a hexagon dance.

• When you apply a magnetic field, it's like a bouncer telling them to change their dance moves.
• The scientists measured if the whole dance floor (the crystal) got bigger or smaller when the dancers changed their moves.

The Findings: A Tale of Two Directions

The most exciting discovery was that the material behaves very differently depending on which way you push it.

• Direction A (The "Stretchy" Side): When they pushed the magnetic field along one direction, the material actually expanded (got longer) as the atoms changed their dance moves.
• Direction B (The "Shrinking" Side): When they pushed from the other direction, the material contracted (got shorter).

This is called anisotropic magnetoelastic coupling. In plain English: The material is "directionally sensitive." It reacts to magnetic forces like a piece of wood that expands differently depending on whether you push it with the grain or against it.

Why does this happen?
The scientists used computer simulations (like a digital twin of the material) to figure out why. They found that the atoms are connected by oxygen bridges (like links in a chain). When the magnetic field changes, these bridges twist and turn, changing the angles between the atoms. This twisting physically stretches or shrinks the whole crystal. It's like a mechanical linkage in a car engine; when one part moves, the whole shape changes.

The Verdict: No Ghost Found (Yet)

So, did they find the Quantum Spin Liquid ghost? No.

Here is what they found instead:

1. Orderly Transitions: As they increased the magnetic field, the material jumped from one orderly state to another. It wasn't a smooth, chaotic fluid; it was a series of sharp, distinct steps.
2. No "Critical Point": In some materials, scientists see a "divergence" (a spike in data) that suggests the material is teetering on the edge of becoming a Quantum Spin Liquid. In NCSO, they saw a small bump, but it didn't get bigger as the temperature dropped. This suggests the "ghost" isn't hiding there.
3. The "Step" Mystery: At very low temperatures, the magnetization (how magnetic the material is) showed little "steps" or jumps. The scientists think this isn't a new quantum state, but rather the material getting "stuck" in temporary, messy configurations (like a door that gets stuck halfway open and then suddenly snaps shut).

The Conclusion

The team concluded that Na₃Co₂SbO₆ is a fascinating material, but it is not the Quantum Spin Liquid they were hoping for.

Instead of a chaotic, fluid quantum state, NCSO is a very structured, highly directional magnet. It shows us that the relationship between magnetism and the physical shape of the crystal is incredibly complex and depends entirely on the direction you look at it.

The Takeaway:
While they didn't find the "holy grail" of quantum computing materials in this specific sample, they learned a lot about how these honeycomb magnets work. They proved that the shape of the crystal is tightly linked to how the spins behave, and that sometimes, what looks like a mysterious quantum effect is actually just the material getting stuck in a temporary, messy state.

It's like searching for a unicorn in a forest and finding a very interesting, color-shifting chameleon instead. You didn't find the unicorn, but the chameleon is still pretty amazing.

Technical Summary

Here is a detailed technical summary of the paper "Anisotropic magnetoelastic coupling in the honeycomb magnet Na3Co2SbO6".

1. Problem Statement and Motivation

The search for Quantum Spin Liquid (QSL) states, particularly those described by the Kitaev honeycomb model, is a central theme in modern condensed matter physics. While 4d/5d transition metal compounds (like $\alpha$-RuCl$_3$) are leading candidates, they often suffer from structural disorder and complex phase diagrams. Recently, 3d$^7$ Co$^{2+}$ honeycomb cobaltates (e.g., Na$_3$Co$_2$SbO$_6$, Na$_2$Co$_2$TeO$_6$) have emerged as promising alternatives due to their ability to be synthesized as high-quality, large single crystals.

Na$_3$Co$_2$SbO$_6$ (NCSO) is a specific focus of interest because:

• It exhibits a honeycomb lattice of CoO$_6$ octahedra with nearly 90° Co–O–Co bond angles, theoretically favoring dominant Kitaev interactions.
• It possesses a complex magnetic ground state (recently identified as a double-q structure rather than a simple zigzag) and strong in-plane magnetic anisotropy.
• Previous studies identified field-induced transitions at critical fields $B_{c1}$ and $B_{c2}$, but the nature of the high-field state (specifically whether it hosts a field-induced QSL or quantum criticality) remains debated.
• The Gap: There is a lack of understanding regarding the magnetoelastic coupling (interaction between spins and the crystal lattice) in NCSO. Understanding this coupling is crucial for determining the nature of phase transitions and the potential existence of exotic quantum states.

2. Methodology

The authors employed a multi-probe experimental approach on high-quality, twin-free single crystals of NCSO to map the field-temperature ($B-T$) phase diagram down to sub-Kelvin temperatures.

• Sample Preparation: Twin-free crystals (1 $\times$ 1 $\times$ 0.1 mm$^3$) were grown via an optimized flux method. Twinning was rigorously excluded using Laue diffraction and angle-dependent susceptibility measurements.
• Thermal Expansion & Magnetostriction: High-resolution capacitance dilatometry was used to measure length changes ($\Delta L/L_0$) along the $c^*$-axis (out-of-plane) under magnetic fields applied along the in-plane $a$ and $b$ directions. This allowed for the calculation of the linear thermal expansion coefficient ($\alpha_{c^*}$) and magnetostriction ($\lambda$).
• Magnetization ($M$): Measurements were performed using a Quantum Design MPMS3 (with He$_3$ attachment for $T < 2$ K) to detect magnetization steps and hysteresis.
• Specific Heat & Magnetic Grüneisen Parameter: Specific heat ($C_p$) was measured to identify phase transitions. The magnetic Grüneisen parameter ($\Gamma_B$), derived from the magnetocaloric effect, was used to probe quantum criticality.
• Theoretical Calculations: Density Functional Theory (DFT) calculations (VASP code with DFT+U+SO) were performed to model the total energy dependence on lattice strain and spin configurations (ferromagnetic, zigzag, and double-q states).

3. Key Contributions

• First comprehensive magnetoelastic study of NCSO: The paper provides the first detailed mapping of the out-of-plane lattice response to in-plane magnetic fields in this material.
• Resolution of the $B_{c2}$ transition nature: The study clarifies that the transition at $B_{c2}$ changes from second-order at higher temperatures to first-order at sub-Kelvin temperatures, evidenced by hysteresis in magnetization and magnetostriction.
• Exclusion of Quantum Criticality/QSL: The authors provide strong thermodynamic evidence against the existence of a field-induced Quantum Spin Liquid state or a genuine quantum critical point near $B_{c2}$.
• Mechanism of Anisotropy: The paper links the observed anisotropic lattice response directly to changes in Co–O–Co bond angles via ab initio calculations, rather than simple single-ion effects.

4. Key Results

A. Anisotropic Magnetoelastic Coupling

• Strong Anisotropy: The lattice response along the $c^*$-axis is highly dependent on the field direction ($a$ vs. $b$).
  • For $B \parallel a$: As the field approaches $B_{c2}$ ($\sim 1.8$ T), the lattice exhibits a broad minimum in length change, followed by lattice expansion upon cooling in high fields.
  • For $B \parallel b$: The lattice continues to shrink monotonically as temperature decreases, even above $B_{c2}$ ($\sim 0.85$ T).
• Origin: DFT calculations reveal that different spin configurations (FM along $a$ vs. $b$, zigzag, double-q) result in different equilibrium $c^*$-axis lengths. This is driven by the anisotropic variation of Co–O–Co bond angles under strain, which modulates the exchange interactions. The bond angles change by $\sim 0.4^\circ$ depending on spin direction, a signature of anisotropic exchange (Kitaev/off-diagonal terms).

B. Nature of Phase Transitions

• $B_{c1}$ Transition: Observed as a first-order transition with clear hysteresis in both magnetization and magnetostriction, persisting up to $\sim 4-5$ K.
• $B_{c2}$ Transition:
  • At $T > 2$ K, the transition appears second-order.
  • At $T < 2$ K (specifically at 0.4 K), a small hysteresis emerges in $M(B)$ and magnetostriction for $B \parallel b$, indicating a crossover to a first-order transition.
  • The magnetic Grüneisen parameter ($\Gamma_B$) shows a divergence-like feature near $B_{c2}$, but the prefactor is small, and the divergence weakens upon cooling (unlike a true quantum critical point where it should strengthen).

C. Low-Temperature Magnetization ($T = 0.4$ K)

• The magnetization curve $M(B)$ exhibits step-like features (sharp jumps in $dM/dB$), particularly for $B \parallel b$.
• These steps suggest abrupt changes in spin configuration rather than gradual canting. They are interpreted as signatures of metastable states and slow domain/defect dynamics, similar to observations in other Co-based magnets like BaCo$_2$(AsO$_4$)$_2$.

D. Phase Diagram Construction

• A detailed $B-T$ phase diagram was constructed.
• Phase I: Double-q magnetic order (zero field).
• Phase II: Field-induced ordered states ($1/3$-AFM for$B \parallel a$, canted zigzag for$B \parallel b$).
• High Field: A field-polarized state. The crossover from the saturated state to the paramagnetic state is marked by a broad maximum/minimum in $\alpha_{c^*}$.

5. Significance and Conclusion

This work fundamentally advances the understanding of honeycomb cobaltates by:

1. Ruling out QSL: It definitively excludes the presence of a field-induced QSL state in Na$_3$Co$_2$SbO$_6$ near $B_{c2}$, as there is no thermodynamic evidence for quantum criticality (the divergence of $\Gamma_B$ is weak and temperature-suppressed).
2. Clarifying Lattice-Spin Coupling: It demonstrates that the magnetoelastic response is dominated by the modulation of Co–O–Co bond angles, providing a microscopic mechanism for the strong in-plane anisotropy observed in these materials.
3. Refining the Phase Diagram: It establishes that the high-field transitions in NCSO are first-order at low temperatures, driven by metastable states and domain dynamics, rather than a continuous quantum phase transition.

The study concludes that while NCSO is a rich platform for studying anisotropic magnetic interactions and Kitaev physics, the intermediate-field regime is governed by classical magnetic order and metastability, not a quantum spin liquid.