CC BY
42ISSN 2072-5981 doi: 10.26907/mrsej
aänetic Resonance in Solids
Electronic Journal
Volume 21 Special Issue 3 Paper No 19303
1-6 pages
2019
doi: 10.26907/mrsej-19303
http: //mrsej. kpfu. ru http: //mrsej. ksu. ru
Established and published by Kazan University Endorsed by International Society of Magnetic Resonance (ISMAR) Registered by Russian Federation Committee on Press (#015140),
August 2, 1996 First Issue appeared on July 25, 1997
© Kazan Federal University (KFU)*
"Magnetic Resonance in Solids. Electronic Journal" (MRSey) is a
peer-reviewed, all electronic journal, publishing articles which meet the highest standards of scientific quality in the field of basic research of a magnetic resonance in solids and related phenomena.
Indexed and abstracted by Web of Science (ESCI, Clarivate Analytics, from 2015), Scopus (Elsevier, from 2012), RusIndexSC (eLibrary, from 2006), Google Scholar, DOAJ, ROAD, CyberLeninka (from 2006), SCImago Journal & Country Rank, etc.
Editor-in-Chief Boris Kochelaev (KFU, Kazan)
Honorary Editors
Jean Jeener (Universite Libre de Bruxelles, Brussels) Raymond Orbach (University of California, Riverside)
Executive Editor
Yurii Proshin (KFU, Kazan) mrsej@kpfu. ru
This work is licensed under a Creative Commons Attribution-Share Alike 4.0
International License.
This is an open access journal which means that all content is freely available without charge to the user or his/her institution. This is in accordance with the BOAI definition of open access.
Editors
Vadim Atsarkin (Institute of Radio Engineering and Electronics, Moscow) Yurij Bunkov (CNRS, Grenoble) Mikhail Eremin (KFU, Kazan) David Fushman (University of Maryland, College Park) Hugo Keller (University of Zürich,
Zürich)
Yoshio Kitaoka (Osaka University,
Osaka)
Boris Malkin (KFU, Kazan) Alexander Shengelaya (Tbilisi State University, Tbilisi) Jörg Sichelschmidt (Max Planck Institute for Chemical Physics of Solids, Dresden) Haruhiko Suzuki (Kanazawa University, Kanazava) Murat Tagirov (KFU, Kazan) Dmitrii Tayurskii (KFU, Kazan) Valentine Zhikharev (KNRTU,
Kazan)
Technical Editors of Issue
Maxim Avdeev (KFU) Alexander Kutuzov (KFU)
* In Kazan University the Electron Paramagnetic Resonance (EPR) was discovered by Zavoisky E.K. in 1944.
Short cite this: Magn. Reson. Solids 21, 19303 (2019)
doi: 10.26907/mrsej-19303
Magnetoelectric coupling in noncollinear ferrimagnets
M.V. Eremin
Kazan Federal University, Kremlevskaya 18, 420008 Kazan, Russia
E-mail: [email protected]
(Received March 30, 2019; accepted March 31, 2019; published April 19, 2019)
The analytical expressions for the energy coupling of exchange coupled spins with electric field, induced by odd crystal field has been derived in the third order perturbation theory, combining the action of electric field, spin-orbit and exchange interactions. The magnetoelectric coupling due to the partial replacement of the Cu2+ positions by Li+ in LiCuVO4 has been discussed.
PACS: 77.80.-e, 76.30.Fc.
Keywords: multiferroics, LiCuVO4, incommensurate magnetic structure.
In honor of Professor Boris Ivanovich Kochelaev
1. Introduction
The materials, which simultaneously exhibit a long-range magnetic order and a spontaneous electric polarization, have attracted a great attention due to their potential applications in promising electronic devices as well as the source of fundamental knowledge of the nature of spins and electric field coupling [1-3]. Recent reviews of such kind materials and possible microscopic mechanisms are presented in Refs. [4-7]. Among most relevant mechanisms are magnetostriction and the so called inverse Dzyaloshinskii-Moriya (DM) mechanism [8,9]. Currently the suggested expressions of spin-electric field coupling have a phenomenological status since the estimated parameters of suggested spin-Hamiltonian parameters are still do not to correspond those which determined experimentally [10]. Moreover, in particular for LiCuVO4 and LiCu2O2 the inverse Dzyaloshinskii-Moriya (DM) mechanism [11,12] fails to reproduce the direction of spontaneous electric polarization induced by spin order at T < TN. Therefore, the phenomenological expressions for spin-electric field coupling should be revised as well.
In the present communication we study the effect of an electric field, induced by an odd crystal field, on a system of exchange-coupled spins. It is known that odd crystal field plays a key role in the theory of electric field effect in electron paramagnetic resonance [13] and in the theory of induced electric dipole transitions in optical spectroscopy [14]. Recently it was found that its role is important for explanation of the peculiarity of magnetoelectric coupling in FeCr2O4 ferrimagnet [15]. The ground state of Fe2+ in FeCr2O4 possesses the orbital degree of freedom. In present communication we study the exchange coupled systems like LiCuVO4 and LiCu2O2 in which the ground states of magnetic ions are nondegenerate and assume that the energy intervals between ground and excited states are much larger with respect to exchange interactions.
2. Matrix elements of electric field enhanced by odd crystal field
The energy operator of interaction of an electron with applied electric field is written as follows:
V = E (nl C<?> 10 , (1)
where Cyfc) = -\/4n/(2k + 1) Ykq are the spherical tensor components, Dq,1) are quantities related to the applied electric field components; = |e| rezEz, = T |e| r (e
T ieyEy) /x/2,
ex, ey, ez are dielectric permittivity parameters. For the case of the isotropic medium ex = ey =
ez = (e + 2) /3.
The energy operator corresponding to the odd crystal field acting on the magnetic ions can be written in the same form:
Hcr = e Bfc)< (nl Cf | e) ac. (2)
The odd crystal field leads to the admixture of the opposite-parity excited 3dn-14p states to the ground states of the 3dn electron configuration. As a result, it induces the effect of the applied electric field on the orbital degrees of freedom of 3d-electrons.
The effective operator describing the coupling of orbital states for 3d-electrons with enhanced electric field is derived in the second order of perturbation theory. Using second quantization technique one has
He = -e BkD1 )a+ (n| Cf) | £') ac, a+ (el C^ | n') a^ + h.c.
Apd
= eEM1)a+ [(nl Cf |e) (el cq,1) |n') + (nl C1 le) (el Cf |n')] a* (3)
= e a+ (nl he | n') an'.
Here for short we have introduced the quantities = —B^k)/ l Apdl, where Apd is the energy interval between excited 3dn-14p and ground 3dn states. We have calculated the expressions for the matrix elements in the basis of wave functions |3dmi) which are given below:
(2l He l2) = (—2l He l—2) = || (Ei!)D-!> + E-> D11)) — | (e«1)D-1,) + E-'d^ ,
(2 l He l1) = — yV/|E«3)D-'1) — 33EM" - ^D«')E<3)
+ ^ (E<1)D<1) + D^E™) ,
(2 l He|0) = 1 EM" + ^E<3)D«1) + E<3)D-1' — ^«> (2 |He|—1) = — 3 /1 Ef)D<1) — 7 /| EfD««, (2 IHeI—2) = 6 /2E<3)D<1),
(1 | He |1) = | (2E<1)D<1) — E<1)D-» + E-M") — ||E<3)D<1)
|| 1 - 1 - 1 1 (1 I He I 0) = |/1 (E'1)D<1) + E(1)D01)) + fE<3)D<1) + E03)D11) + 7 /| EfD-1',
(1 IHeI—1) = — | E11) D11) — 4§ E13)D11) — 7 /|E<3) D01),
(1 IHeI—2) = 3 /3E<3)D<1) + 3 /1 E<3)D<1),
(4)
(|)
<0| He |0) = 15 (4E03)D01) - E^D-1} - E-M^ +
+ 32E01)D01) + (E(3)D-11 + E-31Ds1)) ,
<0| He |-1) = -i^i (E01)D(1) + D01)E(1)) - ^E^
5
—E03)D(1) - 1 ^/2e23)D(11, 35 0 1 7 V 5 2 -1'
<0| He |-2) = -2 J2E<')D<') + 1 E<3)D<') + E«3)+ Ef D/
7 V 5
(6)
,(3) n(1)
7 V 5
<-1| He |-1) = <1| HE |1),
<-1| He |-2) =
2) = -f (e01)D(1) + E(1)D01)) +3352E03)D(1)
1 ^ e23)d(11 + ^ e03)d(1),
7 V 5 2 -1 + 35 0 1 ,
(7)
+ -\ -j
+ 7 V 5
<-2| He |-2) = <2| HE |2).
Combining these matrix elements with those which are allowed for spin-orbit and exchange interaction between magnetic ions one can calculate the coupling energy of the spins with electric field, using third order perturbation theory.
3. Effective Hamiltonian for interaction of exchange coupled spins with electric field
Let us consider two magnetic ions (a and b) with one electron (hole) per each site. The wave function of the ground state in second quantization representation is written as
|n) = a+b+ |0),
(8)
where Greek symbols denote occupied orbital states. Quantum numbers of spins are dropped for short. The states, when the electron at site a (b) occupies the excited state t are
in
= a+b+ |0),
In
= a+b+ |0).
(9)
Below, in contrast to Ref. [15], we consider the cases when the energy intervals between excited and ground state are larger with respect to spin-orbit and exchange interactions. Those energy intervals are formed by even crystal fields.
Two types of terms can be selected in third order perturbation theory. The first one is that when matrix elements of operators HE(a) and Hso(a) are nondiagonal, whereas the exchange coupling operators H6X(t^) = Jr^(sasb) and Hex(p^) = JP^(sasb) are diagonal, i.e. they are acting within the equal energy states. One can speculate about the effect of the exchange coupling of magnetic ion with another one in excited state. The second type term of virtual excitation process is that when the matrix elements <p| HE(a) |p) and <^| HE(b) are diagonal. These matrix elements can be interpreted as a result of electric polarisation of the ions in their excited states.
Therefore we can discuss two possible mechanisms of the spins and electric field coupling. The effective spin-Hamiltonian, which corresponds to the first case, is written as
H (1) _ • Ao Jp<
Heff _ ' 2Apr
|t ) <t | He (a) |p) + h.cA [sa x s6]
+ i
. AbJ
<p
— i
2A2
Ao Jt<
2A2
2AP<
4a) |p) <p| He(b) |p) + h.cA[s6 x s0]
i
. AbJ
PM
2A2
l0a) |t) <T| He |p)- h.c.} [So x Sb]c
1(a) |p) <p| He(b) |p) - h.c. ^ [Sb x So]
Here the index a = ; Aa and Ab are spin-orbit coupling parameters, ApT and are the energies of excitations at the sites a and b, respectively.
The spin-Hamiltonian, containing electric polarization of the ground and excited states is
Hff _ iAo <P|2HE(a)|P) {<P| l0a) |T) Jt<,p< + h.c.} [So x Sb]a
+ i
;Ab <p| He(b) |p)
2A2
i
l(a) |p) JMP,<P + h.c^ [sb x So]c
iAo <T^E(a)|T) {<p| ioa) |T) Jt<,p< - h.c.} [So x Sb]£
(11)
TP
i
.Ab <p| He(b) |p)
2A2
l(a) JMP,<P - h.c^ [Sb x So]a.
Note that deriving the expression (11) we used the exchange coupling operators as follows:
Hex (pt) _ Jp<,T<a+aTb+b< + h.c., Hex (^>p) _ J^p.T^b+b^a+ap/ + h.c.
(12)
As one can see from (10) and (11), both expressions vanish if the direction of spins at site a and b states are parallel. The vector product of spins reminds those which enter in the inverse Dzyaloshinskii-Moriya mechanism [8,9]. However, quantities which stay in front of components [Sb x So]a have different nature and values. Moreover, the expressions (10), (11) do not vanish for number of noncollinear spin systems in which the Dzyaloshinskii-Moriya interaction is absent.
4. The origin of spin-electric field coupling in LiCuVO4
The nature of spin-electric field coupling in LiCuVO4 has been discussed in number of papers [16-20]. In this compounds copper spins (S _ 1/2) build up the chains along the b axis. Below Neel temperature (TN _ 2.4K), spiral structures in Cu-O chins are set due to the competition between the nearest neighbor ferromagnetic and next-nearest neighbor antiferromagnetic exchange interactions along the chain. Simultaneously, an electric polarization appears along the a axis [16,17].
The symmetry analysis of crystal structure rules out the presence of Dzyaloshinskii-Moriya interaction [18,19]. Authors of Ref. [20] suggested to take into account the non-stoichiometry of LiCuVO4. Indeed, the ionic radius of Li+ and Cu2+ ions are very close to each other. As a result even the best samples of LiCuVO4 show up the non-stoichiometry and structural defects [21,22]. The model of the structural defect considered in Ref. [20] is that when the
a
a
a
Figure 1. Fragment of Cu-O chain with Li+. Arrows illustrate the direction of odd crystal fields on copper ions.
Li+ position is occupied by Cu2+ ion (out-of-chain Cu2+ centers). Here we focus on the other relevant model of the structural defect in LiCuVO4, which is schematically shown in Fig. 1; i.e. Li+ sets instead Cu2+ in Cu-chain.
As one can see there is strong odd crystal field action on 3d-hole at Cu2+ site neighboring to Li1+ ion. In the coordinate system with the y axis directed along the Cu-O bonds the ground state of Cu2+ is the hole state |xy) = |Z). The nearest excited state is |x2 - y2) = |e). The matrix element <Z| lz k) = 2i is the largest one among others. On the other hand there is matrix element <e| He (a) |Z) when electric field is oriented along the a axis. This follows from the symmetry selection rule. Indeed, in this case the operator He (a) contains the term which transforms as y3x. Therefore one has the following contribution to the energy of spin-electric field coupling per one Li+ defect position:
Hff = -2i
A
<Z | lZa) k) <e| He (a) |Z) JeC [s« x sb]z = P*Ex[sa x s^
(13)
Derived phenomenological form, of the expression (13) well corresponds to experimentally observed situation in LiCuVO4 [16-18]. At T < TN the helical spin order sets in xy plane and simultaneously electric polarization appears along the a axis.
In conclusion, the analytical expressions for the energy coupling of exchange coupled spins with electric field, induced by odd crystal field has been derived in the third order perturbation theory, combining the action of electric field, spin-orbit and exchange interactions. As an example, I have analyzed the magnetoelectric coupling, caused by in-plane defects for LiCuVO4 and found that the resulting expression for the binding energy of spins with an electric field corresponds, at least qualitatively, to the results of experimental studies.
Acknowledgments
I am grateful to Kirill Vasin for his help. This work is partially supported by the Ministry of Science and Higher Education of the Russian Federation, project no. 3.6722.2017/8.9 for the Kazan Federal University.
References
1. Tokura Y. Science 312, 1481 (2006)
2. Cheong S.-W., Mostovoy M. Nat. Mater. 6, 13 (2007)
3. Khomskii D. Physics 2, 20 (2009)
4. Tokura Y., Seki S., Nagaosa N. Rep. Prog. Phys. 77, 076501 (2014)
5. Spaldin N.A. Fiebig M. Science 309, 391 (2005)
6. Fiebig M., Lottermoser T., Meier D., Trassin M. Nat. Rev. Mater. 1, 16046 (2016)
7. Spaldin N.A. MRS Bulletin 42, 385 (2017)
8. Katsura H., Nagaosa N., Balatsky A.V. Phys. Rev. Lett. 95, 057205 (2005)
9. Sergienko I.A., Dagotto E. Phys. Rev. B 73, 094434 (2006)
10. Solovyev I.V. Phys. Rev. B 95, 214406 (2017)
11. Moskvin A.S., Drechsler S.-L. Phys. Rev. B 78, 024102 (2008)
12. Moskvin A.S., Drechsler S.L. Eur. Phys. J. B 71, 331 (2009)
13. Judd B.R. Phys Rev. B 127, 750 (1961)
14. Mims W.B. The Linear Electric Field Effect in Paramagnetic Resonance. Clarendon Press, Oxford (1976)
15. Eremin M.V. JETP Letters 109, 249 (2019)
16. Naito Y., Sato K., Yasui Y., Kobayashi Yu., Kobayashi Yo., Sato M. J. Phys. Soc. Jpn. 76, 023708 (2007)
17. Yasui Y., Naito Y., Sato K., Moyoshi T., Sato M., Kakurai K. J. Phys. Soc. Jpn. 77, 023712 (2008)
18. Schrettle F., Krohns S., Lunkenheimer P., Hemberger J., Buttgen N., Krug von Nidda H.-A., Prokofiev A.V., Loidl A. Phys. Rev. B 77, 144101 (2008)
19. Mourigal M., Enderle M., Kremer R.K., Law J.M., Fak B. Phys. Rev. B 83, 100409(R) (2011)
20. Moskvin A.S., Drechsler S.-L. EPL 81, 57004 (2008)
21. Prokofiev A.V., Vasilyeva I.G., Ikorskii V.N., Malakhov V.V., Asanov I.P., Assmus W. J. Solid State Chem. 177, 3131 (2004)
22. Grams Ch.P., Kopatz S., Bruning D., Biesenkamp S., Becker P., Bohaty L., Lorenz Th., Hemberger J. Scientific Reports 9, 4391 (2019)
