
A theoretical model for tellurite-sulfates Na2Cu5(TeO3)(SO4)3(OH)4and
K2Cu5(TeO3)(SO4)3(OH)4
I. L. Bartolom´e,1L. Errico,1, 2 V. Fernandez,1M. Matera,1, ∗A.V. Gil Rebaza,1and C.A. Lamas1, †
1IFLP - CONICET. Departamento de F´ısica, Facultad de Ciencias Exactas.
Universidad Nacional de La Plata,C.C. 67, 1900 La Plata, Argentina.
2Universidad Nacional del Noroeste de la Provincia de Buenos Aires (UNNOBA),
Monteagudo 2772, CP 2700 Pergamino, Buenos Aires, Argentina
A theoretical model for two new tellurite-sulfates, namely Na2Cu5(TeO3)(SO4)3(OH)4and
K2Cu5(TeO3)(SO4)3(OH)4is determined to be compatible with ab-initio calculations. The results
obtained in this work show that some previous speculations in the literature about the couplings
are correct, obtaining a model with a mixture of ferromagnetic and antiferromagnetic couplings.
We use a combination of numerical techniques to study the magnetic properties of the model. Our
numerical calculations based on the density-matrix renormalization group method reveal that the
system presents Ising-like magnetization plateaux at rational values of the saturation magnetization.
PACS numbers: 05.30.Rt,03.65.Aa,03.67.Ac
I. INTRODUCTION
Recently, Yingying Tang et al.1synthesized
by hydrothermal reaction, two new tellurite-
sulfates (TS) with a distorted Kagom´e strip
structure: Na2Cu5(TeO3)(SO4)3(OH)4and
K2Cu5(TeO3)(SO4)3(OH)4(Na-TS and K-TS in the
following). In both compounds, the magnetic behavior
is associated with the single unpaired electron associated
with each Cu+2 ions, localized over a 1D kagom´e strip
sub-lattice. This particular geometry corresponds to
the one dimensional version of the paradigmatic two
dimensional Kagom´e lattice for which some experimen-
tal realizations for S= 1/2 as the Herbertsmithite
ZnCu3(OH)6Cl22, the α-vesignieite BaCu3V2O8(OH)23,
and [NH4]2[C7H14N][V7O6F18]54were estudied.
The crystal structure of the compounds is schematized
in Fig. 1 and the simplified magnetic geometry we con-
sider is shown in Fig. 2. We show that several magnetic
properties like magnetic plateaux are determined by the
geometry of the plaquette.
The synthesis of these compounds has aroused great in-
terest in the study of the magnetic phase diagram of mod-
els with this Kagom´e strip geometry5–14. In this sense,
the presence of magnetization plateaux15, a Haldane-like
phase9and localized magnon crystal phases have been
detected6,16. The studies carried out so far describe gen-
eral phase diagrams in a parameter space that, a priori, is
not related to the couplings that describe these materials.
Improving the theoretical description then requires esti-
mating the coupling constants of the effective magnetic
model. As proposed by Noodelman17, a way to determine
these coupling constants is by comparing the spectrum of
the reduced model to those obtained by setting the corre-
sponding magnetic configurations in Density Functional
Theory (DFT) based calculations. The original method
was successfully applied in the literature to compute the
magnetic coupling constants of many compounds. How-
ever, as the number of coupling constants and atoms in
FIG. 1: Crystal structure corresponding to
Na2Cu5(TeO3)(SO4)3(OH)4(Na-TS). The compound
K2Cu5(TeO3)(SO4)3(OH)4is isostructural with Na-TS.
the supercell grows, the direct application of the method
becomes challenging: since the number of possible mag-
netic configurations grows exponentially with the num-
ber of magnetic atoms, and the evaluation of the energy
of each configuration is computationally expensive, to ex-
haust the full set of magnetic configurations becomes im-
practical even for a small number of magnetic atoms. On
the other hand, choosing a small set of magnetic configu-
rations could introduce a large bias in the determination
of the coupling constants. To overcome these issues, a
novel strategy based on Noodelman’s breaking symmetry
method was proposed18. In this work, that methodology
is used to determine the couplings in the magnetic model
describing the tellurite-sulfates. A discussion about this
coupling determination is presented and we show that
the S= 1/2 Heisenberg model with these couplings de-
scribes the magnetic properties of the system and allows
to determine qualitatively the behavior of the magnetic
transitions.
Inspired by the experimental determination of the
atomic distance we propose a model with five different
arXiv:2210.15416v1 [cond-mat.str-el] 27 Oct 2022