Na2Cu2TeO6: Quasi-1D Cuprate Tellurate

• Na2Cu2TeO6 is a layered cuprate tellurate defined by a distorted honeycomb lattice that converts into weakly coupled one-dimensional spin chains via Jahn–Teller distortion.
• Advanced studies using DFT+U, neutron scattering, and other techniques reveal alternating ferromagnetic and antiferromagnetic exchanges with a characteristic spin gap around 127 K.
• Experimental and modeling approaches indicate a Mott insulating state with cuprate-like charge transfer behavior and incipient hole pairing tendencies.

Na$_2$Cu$_2$TeO$_6$ is a layered copper tellurate whose Cu$^{2+}$ network is crystallographically a distorted honeycomb lattice within Cu$_2$TeO$_6$ layers, but whose low-energy magnetic sector is more naturally described as a set of weakly coupled, dimerized $S=\tfrac12$ chains running along the crystallographic $b$ axis. Across neutron scattering, diffraction, susceptibility, DFT+$U$, Wannier downfolding, DMRG, Lanczos, specific-heat, and UV–Vis studies, it is described as a spin-gap or more generally gapped low-dimensional antiferromagnet, a cuprate-like charge-transfer Mott insulator, and a system in which Jahn–Teller distortion and anisotropic exchange geometry convert nominal honeycomb connectivity into quasi-one-dimensional magnetism (Sato et al., 2014, Lin et al., 2022, Patil et al., 22 Aug 2025).

1. Crystal-chemical framework

Na$_2$Cu$_2$0TeO$_2$1 is reported as monoclinic, space group $_2$2, with alternating Cu$_2$3TeO$_2$4 and Na layers stacked along $_2$5. Within each Cu$_2$6TeO$_2$7 layer, the Cu ions occupy edge-sharing CuO$_2$8 octahedra, while a later diffraction study also describes Cu as forming CuO$_2$9 plaquettes and, simultaneously, as part of axially elongated CuO$_6$0 octahedra due to Jahn–Teller distortion. These descriptions are complementary rather than contradictory: the same Cu$_6$1 environment is being viewed from the perspectives of local square-planar bonding and elongated octahedral coordination (Lin et al., 2022, Patil et al., 22 Aug 2025).

The neutron-diffraction refinement reported lattice parameters $_6$2, $_6$3, $_6$4, and $_6$5. The refined atomic positions were given as Cu on 4g $_6$6, Te on 2a $_6$7, O1 on 8j $_6$8, O2 on 4i $_6$9, and Na on 4h $^{2+}$0. The same study states that the refinement converged without secondary phases, which it interprets as evidence for phase purity (Patil et al., 22 Aug 2025).

The structural literature repeatedly emphasizes a dual description. In one language, the compound contains a distorted honeycomb Cu$^{2+}$1 lattice with Te$^{2+}$2 ions at the centers of Cu hexagons. In another, magnetically more consequential language, the Cu ions form chains along $^{2+}$3 with alternating couplings along the chain and a weaker interchain coupling. The distinction is attributed to monoclinic distortion, preferential Cu–O bonding, and strong Jahn–Teller activity of Cu$^{2+}$4, which break ideal honeycomb equivalence and generate exchange anisotropy (Sato et al., 2014, Patil et al., 22 Aug 2025).

Selected bond lengths and angles reinforce that distortion. Reported distances include Te–O1 $^{2+}$5 \AA, Te–O2 $^{2+}$6 \AA, Cu–O1 $^{2+}$7 \AA, Cu–O2 $^{2+}$8 \AA, an elongated Cu–O1 $^{2+}$9 \AA, and Cu–Cu distances of $_2$0 \AA\ and $_2$1 \AA. The same work also highlights nearly equal short in-plane Cu–O values of roughly $_2$2 \AA\ $_2$3 and $_2$4 \AA). Angles such as $_2$5 and $_2$6, and $_2$7 and $_2$8, were used to argue for substantial local distortion and anisotropic exchange geometry (Patil et al., 22 Aug 2025).

2. Electronic structure and correlated-insulator description

The Cu valence is treated as Cu$_2$9, corresponding to a nominal $_6$0 configuration with one half-filled hole-like orbital per Cu. DFT identifies the states near the Fermi level as dominated by Cu $_6$1 character strongly hybridized with O $_6$2 states, whereas the other Cu $_6$3 orbitals lie lower in energy and are fully occupied. In the nonmagnetic calculation, the active $_6$4-derived bands extend approximately from $_6$5 to $_6$6 eV, giving a narrow bandwidth $_6$7 eV. Because the local axes are nearly aligned with Cu–O bonds, the $_6$8 orbital spans the CuO$_6$9-like plaquette and hybridizes strongly with O $S=\tfrac12$0, motivating the reduction to an effective single-orbital low-energy model (Lin et al., 2022).

The resulting Wannier orbital is explicitly not a pure atomic Cu $S=\tfrac12$1 state, but an antibonding Cu $S=\tfrac12$2-O $S=\tfrac12$3 molecular orbital. Maximally localized Wannier analysis yields dominant in-plane hoppings $S=\tfrac12$4, $S=\tfrac12$5, and $S=\tfrac12$6, with only very small interlayer hoppings of approximately $S=\tfrac12$7 eV. A central result is that $S=\tfrac12$8, connecting the longer Cu–Cu path, is much larger than $S=\tfrac12$9, even though the $b$0 Cu–Cu distance is much shorter. The paper attributes this to substantial overlap of the effective Cu–O-based Wannier functions along a Cu–O–O–Cu path for $b$1, versus nearly orthogonal overlap along a Cu–O–Cu geometry close to $b$2 for $b$3 (Lin et al., 2022).

Because $b$4 eV is small, the ratio $b$5 is large for physically relevant $b$6, and the interacting system is described as a Mott insulator. In LSDA+$b$7, all magnetic configurations studied are insulating, and the gap increases with $b$8. In the D-AFM state at $b$9 eV, the half-occupied $U$0 sector splits into lower and upper Hubbard-like features, while the oxygen $U$1 DOS remains close to the Fermi level, supporting a charge-transfer-type copper-oxide description. The same study reports a small oxygen moment of about $U$2 on the oxygens along the $U$3 Cu–O–O–Cu path, again reflecting strong hybridization (Lin et al., 2022).

Two distinct gap scales are reported in the literature. The electronic-structure study finds a small nonmagnetic gap of about $U$4 eV, attributed to dimerization of the antibonding $U$5-combination formed by Cu $U$6 and O $U$7 states in the distorted lattice. A later UV–Vis/Tauc analysis estimates a direct optical band gap of approximately $U$8 eV and reports no clear indication of an indirect transition under the stated experimental conditions. The papers present these as different outputs of different analyses rather than as a single unified gap parameter (Lin et al., 2022, Patil et al., 22 Aug 2025).

Quantity	Reported value
Space group	$U$9
Active band window	$_2$0 to $_2$1 eV
$_2$2 bandwidth	$_2$3 eV
Nonmagnetic dimerization gap	$_2$4 eV
$_2$5	$_2$6 eV
$_2$7	$_2$8 eV
$_2$9	$_2$00 eV
Interlayer hopping	$_2$01 eV
Optical gap $_2$02	$_2$03 eV

3. Exchange topology and magnetic Hamiltonians

The magnetic literature uses two related but not identical notational schemes. In the earlier neutron-based description, Na$_2$04Cu$_2$05TeO$_2$06 is modeled as a one-dimensional alternating Heisenberg chain with two relevant nearest exchange couplings,

$_2$07

where $_2$08 acts within a structural dimer and $_2$09 acts between neighboring dimers along the chain. In that convention, the established sign pattern is

$_2$10

namely ferromagnetic within a structural dimer and antiferromagnetic between dimers. The compound is therefore described as a spin-gap system whose magnetism is captured by an alternating $_2$11 Heisenberg chain in a distorted honeycomb structural setting (Sato et al., 2014).

A later DFT+$_2$12 study parameterizes the in-plane exchange network by three couplings $_2$13, $_2$14, and $_2$15, extracted by mapping total energies of ordered states onto a Heisenberg model. In the sign language used for the final comparison to neutron data, $_2$16 and $_2$17 are antiferromagnetic and $_2$18 is ferromagnetic. At $_2$19 eV, identified there as the most realistic value, the extracted exchanges are

$_2$20

to be compared with neutron values

$_2$21

The corresponding ratios are $_2$22 and $_2$23 from DFT, versus $_2$24 and $_2$25 experimentally. In that formulation, the system lies in the regime of weakly coupled alternating AFM–FM chains (Lin et al., 2022).

The microscopic origin of the exchange hierarchy is a recurring theme. The largest AFM coupling does not occur on the shortest Cu–Cu bond; rather, it is associated with a Cu–O–O–Cu super-superexchange path, where overlap of oxygen-tailed Wannier functions is substantial. By contrast, the shorter-bond ferromagnetic interaction is associated with a Cu–O–Cu geometry close to $_2$26, for which Goodenough–Kanamori reasoning favors ferromagnetism because virtual hopping involves nearly orthogonal O $_2$27 orbitals and Hund alignment on oxygen. The weaker $_2$28 is an AFM interchain coupling consistent with the quasi-1D character (Lin et al., 2022).

A later structure–property study recasts the same exchange anisotropy in explicitly geometrical terms. It assigns $_2$29 and $_2$30 as alternating couplings along the chain and $_2$31 between chains, with reported geometries $_2$32: Cu–Cu $_2$33 \AA, Cu–O1–Cu $_2$34; $_2$35: Cu–Cu $_2$36 \AA, Cu–O2–Cu $_2$37; and $_2$38: Cu–Cu $_2$39 \AA, O1–O1 $_2$40 \AA, Cu–O1/O2–Cu $_2$41 and $_2$42. That paper infers a stronger antiferromagnetic $_2$43, a ferromagnetic $_2$44, and a weaker interchain $_2$45, again concluding that the monoclinic distortion converts the honeycomb layer into weakly coupled alternating Cu chains (Patil et al., 22 Aug 2025).

4. Experimental signatures of low-dimensional magnetism

Neutron magnetic scattering provides the sharpest momentum-space evidence for quasi-one-dimensionality. In the reciprocal-space notation $_2$46, with $_2$47 and $_2$48 lying in the honeycomb plane and $_2$49, the low-energy magnetic dispersion and intensity are tracked primarily along $_2$50. At $_2$51, stated to be close to the minimum spin-gap energy, intensity peaks of the triplet excitation were observed for $_2$52 at $_2$53 in the region $_2$54. For $_2$55 and $_2$56, the peaks occurred at the same $_2$57 positions, while no peaks were observed at $_2$58 in scans along $_2$59, $_2$60, and $_2$61. These observations were interpreted as evidence for one-dimensionality and for the spin-correlation pattern associated with ferromagnetic $_2$62 and antiferromagnetic $_2$63. The same neutron analysis emphasizes that the minimum spin-gap energy occurs at the same $_2$64-points where $_2$65 has maxima near $_2$66 meV (Sato et al., 2014).

Thermodynamic and bulk magnetic measurements are consistent with strong short-range antiferromagnetic correlations but no established long-range magnetic order. DC susceptibility measured between $_2$67 and $_2$68 K shows a broad maximum near $_2$69 K, a ZFC–FC bifurcation onset near $_2$70 K, and an FC upturn near $_2$71 K extending to $_2$72 K; the bifurcation is reduced at $_2$73 Oe. Fitting the $_2$74–$_2$75 K inverse susceptibility to the paper’s printed form $_2$76 yields $_2$77 and $_2$78, from which the paper reports $_2$79. The authors interpret the broad maximum as a hallmark of low-dimensional antiferromagnetism and the strongly negative $_2$80 and small reported moment as signatures of strong antiferromagnetic interactions, frustration, and enhanced quantum fluctuations. Earlier work cited in the electronic-structure study reports a spin gap $_2$81 K (Lin et al., 2022, Patil et al., 22 Aug 2025).

Specific heat and diffraction further constrain the ground state. The heat capacity increases smoothly from $_2$82 to $_2$83 K with no sharp anomaly or plateau. A low-temperature fit to

$_2$84

gives $_2$85 and $_2$86, corresponding to a Debye temperature reported as $_2$87 K, while the figure caption gives $_2$88 K. Temperature-dependent neutron diffraction between $_2$89 and $_2$90 K finds no magnetic peaks or magnetic Bragg reflections, and therefore no long-range magnetic order down to $_2$91 K. Over the same interval the unit-cell volume changes from $_2$92 at $_2$93 K to $_2$94 at $_2$95 K, and the Cu–O bond lengths change only slightly, which the authors interpret as weak spin-lattice coupling (Patil et al., 22 Aug 2025).

Isothermal $_2$96 measurements add a more ambiguous low-temperature component. At $_2$97 and $_2$98 K the magnetization is linear in field, at $_2$99 K there is slight nonlinearity at higher fields, and at $_6$00 K the curve is strongly nonlinear, especially from $_6$01 to $_6$02 T, which the authors describe as partial magnetization saturation. The same paper states in the $_6$03 section that an inset hysteresis loop at $_6$04 K shows a coercive field of $_6$05 kOe, whereas its abstract describes the coercivity as negligible. The paper does not resolve that discrepancy (Patil et al., 22 Aug 2025).

5. Effective one-band Hubbard model and hole pairing tendencies

The low-energy carrier sector is modeled by a single-orbital Hubbard Hamiltonian,

$_6$06

where $_6$07 labels the three hoppings and the active filling is $_6$08. The rationale is that the isolated Cu $_6$09-derived low-energy manifold supports a one-band description, although the same work explicitly notes that the FM sign of the shorter-bond exchange is better understood from the underlying Cu–O geometry than from a naive $_6$10 estimate alone (Lin et al., 2022).

The many-body analysis uses DMRG primarily on open chains of length $_6$11, with checks at $_6$12, at least 3000 kept states, and up to 17 finite sweeps; Lanczos exact diagonalization was also performed up to $_6$13. The diagnostic observable is the real-space spin correlation $_6$14. Both methods find rapidly decaying correlations, indicating a gapped short-range dimerized phase rather than long-range magnetic order. The strongest bond correlations are AFM on the $_6$15 dimers, while the weaker interdimer correlations are FM, producing the chain pattern described as AFM–FM and explicitly as an $_6$16-$_6$17-$_6$18-$_6$19 arrangement along the chain. These results are reported as robust for $_6$20 to $_6$21 eV (Lin et al., 2022).

The same effective model was used to study hole doping through the two-hole binding energy

$_6$22

In the paper’s criterion, $_6$23 indicates that two holes lower their energy by forming a bound state. DMRG finds that $_6$24 becomes negative for $_6$25 eV and is most negative around $_6$26 eV, while Lanczos on $_6$27 reproduces the same nonmonotonic dependence. Real-space charge densities show that at small $_6$28 the holes remain spread apart, whereas in the regime with $_6$29 they move closer together and form a tight pair; at $_6$30 eV the pair is most compact. The authors describe these results as incipient or tentative pairing tendencies rather than evidence for a robust superconducting state, emphasizing that $_6$31 is quite small because $_6$32 is tiny and that any superconductivity, if realizable at all, would likely have a very low $_6$33 (Lin et al., 2022).

6. Historiography, notation, and placement within copper tellurates

A notable feature of the Na$_6$34Cu$_6$35TeO$_6$36 literature is that part of the record is organized around priority and notation rather than disagreement over the basic magnetic picture. The 2014 comment on Schmitt et al. states that the sign problem for the alternating-chain exchanges had already been experimentally settled in two earlier papers: a first paper that observed the spin-gap phenomenon but could not yet determine which bond was FM or AFM, while already ruling out $_6$37, and a second paper that unambiguously determined $_6$38 and $_6$39 from the momentum dependence of the dynamical structure factor measured by neutron magnetic scattering. The comment explicitly says that Schmitt et al. confirmed the experimental result theoretically rather than overturning it (Sato et al., 2014).

The later first-principles study uses a different exchange notation, $_6$40, $_6$41, $_6$42, and shows excellent agreement with neutron-derived exchange magnitudes. The coexistence of the $_6$43 and $_6$44 conventions is therefore a feature of the literature rather than a contradiction in the underlying physics. The robust common ground is that Na$_6$45Cu$_6$46TeO$_6$47 is not well described as an isotropic 2D honeycomb magnet, but as a quasi-1D alternating-chain system with one ferromagnetic chain bond, one antiferromagnetic chain bond, and, in the three-coupling description, a much smaller antiferromagnetic interchain coupling (Sato et al., 2014, Lin et al., 2022).

Within the broader copper tellurium oxide landscape, Na$_6$48Cu$_6$49TeO$_6$50 occupies a somewhat unusual position. A selective 2017 review of copper tellurium oxides does not discuss it explicitly at all and therefore provides no direct crystallographic, magnetic, or thermodynamic data for it. The review is nonetheless relevant because it emphasizes a family-level motif: in copper tellurates, Te frequently enters exchange pathways indirectly, so super-superexchange processes such as Cu–O–Te–O–Cu or Cu–O–O–Cu often shape the magnetic energy scale and dimensionality (Norman, 2017).

A plausible family-level implication is reinforced by the hydroflux study of K-based analogs. That work treats K$_6$51Cu$_6$52TeO$_6$53 and K$_6$54Cu$_6$55TeO$_6$56$_6$57H$_6$58O) as alternating-chain antiferromagnetic Heisenberg systems and therefore suggests that layered Cu$_6$59TeO$_6$60 frameworks can remain structurally honeycomb-derived while becoming magnetically chain-dominated. For Na$_6$61Cu$_6$62TeO$_6$63, that comparison does not establish new facts, but it places the compound within a broader $_6$64Cu$_6$65TeO$_6$66 context in which alkali-ion chemistry, interlayer structure, and exchange anisotropy jointly control the crossover between layered crystallography and quasi-one-dimensional magnetism (Iwanicki et al., 2024).