8582 Inorganic Chemistry, Vol. 49, No. 18, 2010
Ben Yahia et al.
quantum number of each spin site (i.e., S = 5/2 in the present
case). Thus, for NaMnVO4, θ can be approximated by
35ð2J1 þ 4J2 þ 4J3 þ 2J4Þ
θ ꢁ
ð3Þ
12kB
The θ value is estimated to be -56.6 and -44.8 by using the
spin exchange parameters from the GGAþU calculations
with U = 4 and 5 eV, respectively. This is in good agreement
with experiment (i.e., -62 K), given that GGAþU electronic
structure calculations generally overestimate the magnitude
of spin exchange interactions by a factor approximately up to
four.25
In NaMnVO4, the intrachain exchange J1 is by far the
strongest exchange, which forms uniform antiferromagnetic
chains along the b-direction. This explains the occurrence of a
broad maximum in the magnetic susceptibility at 24 K. The
spin exchanges J1-J4 are all antiferromagnetic. Therefore,
the exchanges in the (J1, J2, J3) and (J2, J2, J4) triangles are
spin-frustrated (Figure 5a,b), namely, the interchain spin
exchanges are spin-frustrated. This explains why TN is con-
siderably smaller than |θ| in NaMnVO4 (11.8 vs 62 K). It is of
interest to consider the probable magnetic structure for the
ordered antiferromagnetic state below TN. Among the inter-
chain spin exchanges, J4 is stronger than J2 and J3 by a factor
of ∼2. Thus, the magnetic structure shown in Figure 5c, in
which the spins are antiferromagnetically coupled in all J1
and J4 magnetic bonds, is energetically most favorable and
hence is expected to represent the ordered magnetic structure
of NaMnVO4 below TN.
Figure 6. Five ordered spin states of NaMnVO4 employed to extract the
spin exchanges J1-J4 by using a (a, b, 2c) supercell. The unshaded and
shaded circles represent the Mn2þ ions with up-spin and down-spin,
respectively. The two numbers in the parentheses (from left to right) for
each state are the relative energies (in meV per 4 FUs) obtained from the
GGAþU calculations (with U = 4 and 5 eV, respectively).
To determine the values ofJ1-J4, we determine the relative
energies of the five ordered spin states on the basis of density
functional calculations. Our calculations employed the fro-
zen-core projector augmented wave method20 encoded in the
Vienna ab initio simulation packages21 and the generalized-
gradient approximation (GGA)22 with the plane-wave-cutoff
energy of 500 eV and a set of 36 k points for the irreducible
Brillouin zone. To properly describe the effect of electron
correlation in the Mn 3d states, the GGA plus on-site
repulsion method (GGAþU)23 was used with the effective
U values of 4 and 5 eV. The relative energies (per 4 FUs) of
the five ordered spin states obtained from GGAþU calcula-
tions are summarized in Figure 6. Thus, by mapping the
relative energies of the five spin ordered states onto the
corresponding energies expected from the total spin exchange
energies listed in eq 1, we obtain the values of J1-J4
summarized in Table 3.
6. Concluding remarks
NaMnVO4 crystallizes in the maricite-type structure as
does NaMnPO4 and undergoes a three-dimensional antifer-
romagnetic ordering at TN = 11.8 K, with a broad maximum
in the magnetic susceptibility at Tmax = 24 K. The broad
maximum reflects the fact that the strongest spin exchange J1
forms uniform antiferromagnetic chains, and TN is consider-
ably smaller inmagnitude than the Curie-Weiss temperature
of -62 K because the interchain spin exchanges are spin-
frustrated. Our study predicts that, in the ordered magnetic
structure below TN, the spins are antiferromagnetically
coupled in all J1 and J4 magnetic bonds. It would be
interesting to verify this prediction by neutron diffraction
measurements.
To see how reasonable the extracted spin exchanges are, we
calculate the Curie-Weiss temperature θ, which in the mean
field theory24 is related to spin exchanges as
X
SðS þ 1Þ
θ ¼
ziJi
ð2Þ
3kB
i
Acknowledgment. The authors acknowledge the work
that occurred at NCSU funded by the Office of Basic
Energy Sciences, Division of Materials Sciences,
U.S. Department of Energy, under Grant DE-FG02-
86ER45259, and also the computing resources of the
NERSC center and the HPC center of NCSU.
where the summation runs over all nearest neighbors of a
given spin site, zi is the number of nearest neighbors con-
nected by the spin exchange parameter Ji, and S is the spin
(19) (a) Dai, D.; Whangbo, M.-H. J. Chem. Phys. 2001, 114, 2887. (b) Dai,
D.; Whangbo, M.-H. J. Chem. Phys. 2003, 118, 29.
€
(20) Blochl, P. E. Phys. Rev. B 1994, 50, 17953.
(21) (a) Kresse, G.; Hafner, J. Phys. Rev. B 1993, 47, 558. (b) Kresse, G.;
Furthm€uller, J. Comput. Mater. Sci. 1996, 6, 15. (c) Kresse, G.; Furthmu€ller, J.
Phys. Rev. B 1996, 54, 11169.
(22) Perdew, J. P.; Burke, K.; Ernzerhof, M. Phys. Rev. Lett. 1996, 77,
3865.
Supporting Information Available: X-ray diffraction details
(CIF) as well as Tables S1 and S2 and Figures S1 and S2. This
material is available free of charge via the Internet at http://
pubs.acs.org.
(23) Dudarev, S. L.; Botton, G. A.; Savrasov, S. Y.; Humphreys, C. J.;
Sutton, A. P. Phys. Rev. B 1998, 57, 1505.
(24) Smart, J. S. Effective Field Theory of Magnetism; Saunders: Philadelphia,
1966.
(25) (a) Xiang, H. J.; Lee, C.; Whangbo, M.-H. Phys. Rev. B: Rapid
Commun 2007, 76, 220411(R). (b) Koo, H.-J.; Whangbo, M.-H. Inorg. Chem.
2008, 47, 128. (c) Koo, H.-J.; Whangbo, M.-H. Inorg. Chem. 2008, 47, 4779.