--- title: 'Na2Co2TeO6: Layered Cobalt Honeycomb Oxide' url: https://www.emergentmind.com/topics/na2co2teo6 type: topic --- # Na2Co2TeO6: Layered Cobalt Honeycomb Oxide Na\(_2\)Co\(_2\)TeO\(_6\) is a layered cobalt honeycomb oxide that has become a central 3\(d\)-electron candidate for bond-directional Kitaev magnetism. In the hexagonal material most commonly studied, Co\(^{2+}\) ions in edge-sharing CoO\(_6\) octahedra form quasi-two-dimensional honeycomb layers, while Na disorder in the interlayer galleries weakens interplane coupling and complicates the low-temperature magnetic state. Across neutron scattering, ESR, NMR, \(\mu\)SR, thermal transport, torque magnetometry, and Raman studies, Na\(_2\)Co\(_2\)TeO\(_6\) is consistently described as a pseudospin-\(\tfrac12\) honeycomb magnet with strong spin-orbit entanglement and substantial frustration; the main open issues are the precise zero-field order, the correct microscopic Hamiltonian, and the nature of the field-induced intermediate phase [2012.00006]. ## 1. Crystal chemistry, layering, and polymorphism Most studies describe Na\(_2\)Co\(_2\)TeO\(_6\) as a hexagonal layered oxide in which edge-sharing CoO\(_6\) octahedra form two-dimensional honeycomb networks in the \(ab\) plane, separated by Na and TeO\(_6\) layers. A commonly reported crystallographic description is space group \(P6_3 22\) (No. 182), with room-temperature lattice parameters roughly \(a = 5.28\) Å and \(c = 11.14\) Å, or closely related values \(a = b = 5.2744\) Å, \(c = 11.2388\) Å, \(V = 270.77\) Å\(^3\) [2012.00940] [2407.14663]. Earlier neutron refinements in the same structural family gave \(a = 5.2770(2)\) Å, \(c = 11.2231(1)\) Å and emphasized that the three Na sites are only partially occupied, producing substantial interlayer disorder [1608.07389]. The local structural motif is unusually consequential for the magnetism. Co\(^{2+}\) ions occupy slightly distorted octahedra and the Co–O–Co geometry is close to \(90^\circ\), a configuration that is repeatedly invoked in the literature as favorable for bond-dependent exchange. Several reports stress that the interlayer Na distribution is disordered: one refinement gave Na occupancies of about \(67\%\), \(25\%\), and \(8\%\) on three partially occupied sites, while a later NMR-oriented refinement described \(67\%\) occupancy on \(6g\) and \(33\%\) on \(2a\) [1608.07389] [2407.14663]. These data are consistent with the repeated observation that magnetic correlations are much longer ranged within the honeycomb planes than perpendicular to them. The structural literature is not fully uniform. Different studies have described the same hexagonal material in \(P6_3\), \(P6_3/mmc\), or \(P6_3/mcm\) settings, always retaining the essential picture of quasi-two-dimensional Co honeycomb layers separated by Na/Te-based spacers [2203.00282] [2307.16451] [2210.13798]. In addition, a distinct monoclinic polymorph in \(C2/m\) has been reported, with room-temperature lattice parameters \(a = 5.33225(6)\) Å, \(b = 9.20808(8)\) Å, \(c = 5.80718(8)\) Å, \(\beta = 108.90837(88)^\circ\), and \(V = 269.745\) Å\(^3\) [2305.10484]. | Structural variant | Reported setting | Key reported feature | |---|---|---| | Hexagonal Na\(_2\)Co\(_2\)TeO\(_6\) | \(P6_3 22\) | Strong Na disorder; quasi-2D honeycomb layers | | Alternative hexagonal descriptions | \(P6_3\), \(P6_3/mmc\), \(P6_3/mcm\) | Same layered Co-honeycomb motif | | Monoclinic polymorph | \(C2/m\) | Lower \(T_N\), single AFM transition, reduced Na/vacancy disorder | The monoclinic polymorph is not a small perturbation of the hexagonal phase. It shows only one AFM transition at \(9.6(6)\) K and no lower-temperature spin reorientation transitions, in contrast to the hexagonal material where ordering near \(27\) K is often accompanied by additional anomalies near \(15\)–\(17\) K and \(4\)–\(7\) K [2305.10484]. This suggests that interlayer Na disorder and stacking details are not secondary crystallographic complications but integral control parameters for the magnetic phenomenology. ## 2. Local electronic structure and effective Hamiltonians In the hexagonal compound, Co\(^{2+}\) is consistently described as high-spin \(3d^7\) in an approximately octahedral crystal field, with the \(t_{2g}^5 e_g^2\) configuration supporting a low-energy Kramers doublet that behaves as an effective pseudospin-\(\tfrac12\) or \(j_{\rm eff}=\tfrac12\) degree of freedom [2012.00006]. One neutron study placed the first excited \(j_{\rm eff}=\tfrac32\) manifold roughly \(21\)–\(27\) meV above the ground doublet, and observed a spin-orbit excitation near \(21.6\) meV; a single-crystal INS study reported two crystal-field bands around \(22\) meV and \(28\)–\(30\) meV; by contrast, an out-of-plane field study described the lowest Kramers doublet as separated from the excited quartet by about \(6\) meV [2203.00282] [2405.13935]. The spread in quoted splittings reflects that different analyses emphasize different effective descriptions and energy windows. The low-energy spin dynamics are usually modeled by a bond-dependent extended Kitaev Hamiltonian on the honeycomb lattice. A representative form is \[ H=\sum_{\langle ij\rangle_\gamma}\Bigl[K\,S_i^\gamma S_j^\gamma + J_1\,\mathbf S_i\!\cdot\!\mathbf S_j + \Gamma\,(S_i^\alpha S_j^\beta+S_i^\beta S_j^\alpha)+\Gamma'\,(S_i^\alpha S_j^\gamma+S_i^\gamma S_j^\alpha+S_i^\beta S_j^\gamma+S_i^\gamma S_j^\beta)\Bigr] +J_2\sum_{\langle\!\langle ij\rangle\!\rangle}\mathbf S_i\!\cdot\!\mathbf S_j +J_3\sum_{\langle\!\langle\!\langle ij\rangle\!\rangle\!\rangle}\mathbf S_i\!\cdot\!\mathbf S_j, \] where \(\langle ij\rangle_\gamma\) labels nearest-neighbor \(x,y,z\) bonds, \(K\) is the Kitaev exchange, \(J_1,J_2,J_3\) are Heisenberg terms, and \(\Gamma,\Gamma'\) are symmetry-allowed off-diagonal exchanges [2012.00006]. Simpler versions retain only \(J_1\), \(K\), and \(\Gamma\), or omit \(J_2\) to reduce over-parameterization [2012.00940] [2112.12254]. Parameter extraction remains one of the most contested aspects of the material: | Source | Model | Reported parameter set for Na\(_2\)Co\(_2\)TeO\(_6\) | |---|---|---| | [2012.00006] | \(K\)–\(J_1\)–\(J_2\)–\(J_3\)–\(\Gamma\)–\(\Gamma'\) | \(K=-9.0(5)\), \(J_1=-0.1(5)\), \(J_2=0.3(3)\), \(J_3=0.9(3)\), \(\Gamma=1.8(5)\), \(\Gamma'=0.3(3)\) meV | | [2012.00940] | \(J_1\)–\(K\)–\(\Gamma\)–\(J_2\)–\(J_3\) | \(J_1=-2.18(5)\), \(J_2\approx0\), \(J_3=2.50(10)\), \(K=+0.125(25)\), \(\Gamma=+0.125(25)\) meV | | [2112.12254] | \(K\)–\(J_1\)–\(\Gamma\)–\(\Gamma'\)–\(J_3\) | FM-\(K\): \(J_1=-0.20\), \(K=-7.00\), \(\Gamma=0.50\), \(\Gamma'=0.15\), \(J_3=1.40\) meV; dual AF-\(K\): \(J_1=-3.50\), \(K=3.20\), \(\Gamma=-3.00\), \(\Gamma'=2.00\), \(J_3=1.60\) meV | | [2409.07959] | \(K\)–\(J_1\)–\(\Gamma\)–\(\Gamma'\)–\(J_3\), VMC | \(K=1.408\pm0.05\), \(J_1=-1.54\pm0.03\), \(\Gamma=-1.32\pm0.04\), \(\Gamma'=+0.88\pm0.02\), \(J_3=+1.32\pm0.04\) meV | Despite the disagreement in fitted couplings, several recurring conclusions appear. First, a purely isotropic Heisenberg model is repeatedly found inadequate for the observed spectra [2112.12254]. Second, many analyses favor a dominant ferromagnetic Kitaev term with subleading \(J_3\), \(\Gamma\), and \(\Gamma'\) [2012.00006] [2112.12254]. Third, other fits obtain very different balances of \(J_1\), \(J_3\), and anisotropies from similar classes of powder data [2012.00940]. This suggests that the microscopic Hamiltonian is underdetermined by powder-averaged spectra alone and is highly sensitive to assumptions about the magnetic unit cell and the relevant field regime. ## 3. Zero-field magnetic order and correlation anisotropy The most frequently reported zero-field ordered state of the hexagonal compound is zigzag antiferromagnetism with propagation vector \(k=(\tfrac12,0,0)\), or equivalently M-point ordering in the honeycomb Brillouin zone. In this picture, ferromagnetic chains run along one in-plane direction and are coupled antiferromagnetically from chain to chain; early neutron powder refinements placed the ordered moments essentially along the crystallographic \(b\) axis, with negligible \(c\)-axis component within uncertainty [1608.07593] [1608.07389]. Representative refined moments are \(2.95(3)\,\mu_B\) and \(2.5(3)\,\mu_B\) on the two Co sites in one study, and \(2.77(3)\,\mu_B\) and \(2.45(2)\,\mu_B\) in another, both below the spin-only \(3\,\mu_B\) expectation for \(S=\tfrac32\) and therefore suggestive of persistent fluctuations [1608.07593] [1608.07389]. A robust feature of the neutron data is the strong anisotropy of magnetic correlations. The first magnetic reflection is resolution-limited, implying very long in-plane correlation lengths, whereas peaks involving nonzero \(L\) are Lorentzian broadened, consistent with much shorter correlations perpendicular to the layers. One analysis quoted \(\xi_{ab}\gg 100\) Å and \(\xi_c\sim10\)–\(20\) Å, while another obtained \(\xi_c \simeq 29\) Å even at \(1.8\) K [1608.07593] [1608.07389]. Above the ordering temperature, diffuse scattering and reverse Monte Carlo analysis showed that short-range correlations are largely confined to the honeycomb planes, with \(\xi_{ab}\approx12\) Å at \(25\) K and a sign structure of positive nearest-neighbor and negative second- and third-neighbor correlations, consistent with the eventual zigzag motif [1608.07389]. The thermal evolution is more complex than a single transition. Several bulk and diffraction studies place the main ordering temperature near \(T_N \approx 24.8\)–\(27\) K and report additional anomalies around \(16\)–\(17\) K and \(4\)–\(7\) K [1608.07593] [2012.00940]. A single-crystal study further resolved a two-step sequence: a transition at \(T_1 = 31.0\) K into a two-dimensional long-range ordered state characterized by Bragg rods at the six M-points, followed by a \(3D\) ordering transition at \(T_2 = 26.7\) K where sharp Bragg peaks at integer \(L\) appear [2012.08781]. In contrast, the linear-spin-wave analysis of one powder INS study discussed zigzag order below \(T_N \approx 17\) K [2012.00006]. The literature therefore supports the existence of multiple thermodynamic and magnetic scales rather than a universally agreed single \(T_N\). The detailed microscopic structure is also debated. The conventional assignment is a collinear zigzag state, but single-crystal INS and diffraction arguments have been used to support a triple-\(q\) description in which the three zigzag ordering vectors \(q_1\), \(q_2\), and \(q_3\) superpose, giving a larger magnetic Brillouin zone and preserving \(C_3\) in a way not expected for single-domain zigzag order [2012.08781]. A later single-crystal neutron diffraction study described the zero-field order as a triple-\(Q\)-type zigzag structure with \(k=(\tfrac12,0,0)\) [2312.04834]. By contrast, an out-of-plane field study described the zero-field state as a canonical ferrimagnetic state with canted Co moments along the \(c\)-axis [2405.13935]. Taken together, these results show that the existence of long-range order is not controversial, but its detailed symmetry content remains unsettled. ## 4. Excitation spectrum from meV to crystal-field scales Powder INS first established that the ordered state supports gapped low-energy modes emerging from \(Q \approx 0.75\) Å\(^{-1}\), a flat band near \(3\) meV, additional weaker flat bands near \(4.5\) and \(7\) meV, and residual magnetic intensity up to about \(10\) meV [2012.00006]. Another powder study in the \(0.4\)–\(7\) meV window found dispersive magnons with a minimum gap of about \(0.4\) meV at the M point, while high-field ESR observed three AFMR modes \(A\), \(B\), and \(C\) with a zero-field gap \(\Delta \approx 0.4\) meV and an emergent mode \(D\) above \(6\) T with slope \(\simeq 0.15\) meV/T and \(g_{\rm eff}\approx2.6\) [2012.00940]. A later neutron/ESR study instead described a zero-field gap of about \(1\) meV, consistent with \(\Delta/k_B \approx \Delta/h \approx 240\) GHz [2312.04834]. High-resolution single-crystal INS substantially complicated this picture. At \(T=5\) K, at least six non-overlapping spin-wave branches were resolved between \(1\) and \(\sim12\) meV, and the lowest mode was reported to have minima of \(\simeq1\) meV at both M and \(\Gamma\), which the authors argued is incompatible with simple nearest-neighbor Kitaev-zigzag scenarios [2203.00282]. In an effective fit restricted to the lowest branch, they obtained \(J_3 = 1.896(9)\) meV and \(\Delta = 0.170(6)\) meV, promoting the third-neighbor distance as an emergent “soft link” [2203.00282]. This was used to argue that all then-available powder-based HK\(\Gamma\) parameterizations fail to reproduce the full branch multiplicity and \(Q\)-dependence. The low-energy magnons coexist with higher-energy orbital and crystal-field excitations. Neutron spectroscopy identified a spin-orbit transition from \(j_{\rm eff}=\tfrac12\) to \(j_{\rm eff}=\tfrac32\) near \(21.6\) meV in one study and two bands around \(22\) meV and \(28\)–\(30\) meV in another [2012.00006] [2203.00282]. Raman spectroscopy extended this hierarchy to much higher energies, reporting broad features at \(130\) meV and \(136\) meV that sharpen and shift below the \(3D\) ordering transition, together with a weaker feature at \(119\) meV; these were assigned to spin-orbit excitons and multi-phonon processes aided by orbital fluctuations [2012.08781]. A later Raman study further reported low-temperature crystal-field excitations at \(182.9\) cm\(^{-1}\), \(188.3\) cm\(^{-1}\), \(967\) cm\(^{-1}\), \(1110\) cm\(^{-1}\), and \(1223\) cm\(^{-1}\), including a \(\sim3\) cm\(^{-1}\) splitting of a Kramers doublet [2502.03970]. Above \(T_N\), the excitation spectrum remains far from featureless. Single-crystal INS found that finite-energy correlations persist up to \(\sim5T_N\), with dynamic intensity remaining concentrated around the M points and equal-time structure factors describable by a single Co\(_6\) hexagon in zigzag-type or noncollinear “tornado-like” arrangements [2203.00282]. This suggests that local hexagonal-cluster correlations survive well into the paramagnetic regime, constraining theories that attempt to identify the material solely through its ordered state. ## 5. Field-induced phase structure and candidate disordered regimes The field response is strongly anisotropic. For fields perpendicular to the \(b\) axis, one thermodynamic/INS/ESR study reported a coexistence regime for \(6\ {\rm T}