CC BY
113ISSN 2072-5981 doi: 10.26907/mrsej
aänetic Resonance in Solids
Electronic Journal
Volume 21 Special Issue 4 Paper No 19405 1-6 pages
2019
doi: 10.26907/mrsej-19405
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
0 © I This work is licensed under a Creative
* m-we» 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.
Special Editor of Issue
Eduard Baibekov (KFU)
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
Nurbulat Abishev (KFU) Maxim Avdeev (KFU) Eduard Baibekov (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, 19405 (2019)
doi: 10.26907/mrsej-19405
Advances of crystal field theory and exchange charge model
M.G.Brik1'2'3'*, N.M. Avram4'5, C.N.Avram4
1College of Sciences, Chongqing University of Posts and Telecommunications, Chongqing 400065, People's Republic of China 2Institute of Physics, University of Tartu, W. Ostwald Str. 1, Tartu 50411, Estonia 3Institute of Physics, Jan Dlugosz University, Armii Krajowej 13/15, Czestochowa PL-42200, Poland 4Department of Physics, West University of Timisoara, Bd. V. Parvan, No. 4,
Timisoara 300223, Romania 5Academy of Romanian Scientist, Independentei 54, Bucharest 050094, Romania
*E-mail: mikhail.brik@ut.ee (Received May 27, 2019; accepted May 28, 2019; published June 6, 2019)
A short overview of the exchange charge model of crystal field is given. Several most important applications of the model to the calculations of the crystal field parameters, energy levels of impurity ions, parameters of the electron-vibrational interaction and transition intensities supported by the recent literature publications are discussed.
PACS: 78.20.Bh, 71.70.Ch, 71.70.Ej.
Keywords: crystal field theory, transition metal and rare earth ions, crystal field splitting. Preface
This short article was written on occasion of the 80th anniversary of the author of exchange charge model (ECM), Professor Boris Zalmanovich Malkin from Kazan University. The authors of this article are honored to know Boris Zalmanovich personally. We would like to congratulate him sincerely on occasion of his jubilee, which marks his bright scientific career and profound contribution to the development of crystal field theory and deeper understanding of properties of the impurity ions in solids. Happy 80th anniversary, and a lot of joy, good health and happiness, dear Boris Zalmanovich!
1. Introduction
Fascinating variety of colors of gemstones, crystals and glasses has always admired people through the whole long history of different civilizations. Since ancient times people could tune the colors of glasses by changing their constituting elements and/or preparation conditions (of course, without having realized what was happening on the microscopic level), which resulted
in a qualitative understanding of the relations between the colors and chemical composition. Still
a very long way was left to be gone through until it became possible to get a quantitative description of optical properties of solids, understand the role played by the impurities - both intentional or unintentional - in the formation of those optical properties and learn how to control amount of those impurities (or defects) to produce crystals of high optical quality for numerous applications.
Development of quantum mechanics, theory of free atoms/ions spectra and group theory created all required conditions and set the stage for emerging crystal field theory (CFT) that allowed to uncover the main physical mechanisms responsible for varying optical properties of impurities and calculate their energy levels splitting.
Historically, the CFT takes its roots from 1929, when H. Bethe published his famous paper [1], in which he linked the group theory and splitting of the spectral terms of impurity ions in crystal fields (CF) of different symmetries. The next developments were made by Van Vleck [2].
One of the CFT keystones is the possibility of decomposing the CF Hamiltonian into irreducible tensor operators [3,4] as
H = EE Bk ok , (1)
p=0,2,4,6 k=-p
where Ok are the linear combinations of the irreducible tensor operators acting on the angular parts of the impurity ion wave functions. The first term Bpk stands for the crystal field parameters (CFPs). They have dimension of energy, and all information about the point symmetry of the crystal lattice site occupied by an impurity center is "hidden" in those CFPs. The term with p = 0 can be dropped since it just produces the overall shift of the whole energy level scheme and does not affect the relative intervals between the calculated energy levels. In the case of the impurity ions with the d-electron shell p = 2,4, and in the case of the impurity ions with the f-electron shell p = 2, 4, 6.
If the calculations of the Ok matrix elements are reasonably easy to be performed, the CFPs are harder to evaluate. The first model of CFT was the so-called the point charge model (PCM), which considered only the Coulomb or point charge interaction between the impurity ion and all surrounding crystal lattice ions. The PCM equations allow to calculate the splitting of the impurity ions spectral terms. However, it has been recognized immediately that these splittings are severely underestimated and there was an urgent need in the development of more advanced models capable of taking into account more interactions between impurity ions and their neighbors, differences in the electronic structure of different ions, nature of ligands, spatial distribution of their electron densities, formation of chemical bonds between the impurities and ligands etc.
The number of non-zero CFPs in Hamiltonian (1) depends on the point symmetry around an impurity ion; the lists on those non-zero CFPs for different point groups can be found in the literature [5,6].
There are two possible ways to get the numerical values of the CFPs. The first one is based on the non-linear fit of the calculated eigenvalues of the Hamiltonian (1) to the experimentally observed energy levels. In this case, the information about the point symmetry of the impurity ion site should be known before the calculations - this will determine which CFPs are not zero. Then the initial values of those non-zero CFPs should be taken as an input for the fitting procedure. After that the root-mean-square (rms) deviation between the calculated and experimental energy levels should be minimized by varying the CFPs values. This approach is rather mathematical; it allows to get a very small rms value (of the order of 10-20 cm-1 in the case of the rare earth ions, for example). There are, of course, certain disadvantages of this method. First of all, the number of fitting parameters is very large, because all CFPs are treated as independent parameters to be varied freely. Secondly, it may be difficult to get really global rms deviation minimum, since the fitting procedure may end up in one of possible local rms deviation minima. In addition, it may happen - especially when the CF has a rather low symmetry - that several CFPs sets result in the same or very close rms values, making them mathematically indistinguishable and hindering their physical analysis.
The second method of getting the CFPs numerical values is to calculate them from the known crystal structure. This approach is more physical than the previously described one because it does not imply any a priori knowledge of the point symmetry of the impurity ion site. Moreover,
M.G. Brik, N.M. Avram, C.N. Avram
it would allow to decrease considerably the number of the fitting parameters and establish clear trends that link together symmetry of the crystal lattice site occupied by an impurity ion and its spectroscopic properties.
One of the CFT models that give a possibility of calculating CFPs values is the so-called superposition model developed by Newman [7,8]. With the positions of the nearest neighbors of the impurity ion, it is possible to calculate the CFPs. The superposition model if applied to the d- or f -ions employs at least four and six fitting parameters, respectively, the so-called intrinsic parameters Bp and power-law exponents tp (p = 2,4 for the d-ions and 2,4, 6 for the f-ions). Detailed description of this model can be found in a review chapter by Y.Y. Yeung [9].
The exchange charge model (ECM) of crystal field was developed by B.Z. Malkin [10]. It is based even on a smaller (sometimes only one!) number of fitting parameters. The model contains explicit expressions that allow to calculate all CFPs from the crystal structure without any a priori assumption about the CF symmetry. The number of ions taken for such calculations is practically unlimited and can include contributions of ions located very far from the impurity ion site, which allows to achieve convergence of the oscillating sign-changing crystal lattice sums.
The ECM has such a wide potential of applications to different systems/phenomena, that it makes it one of the most successful CF models developed so far. In the next paragraph we briefly discuss the model and its usage.
2. A short description of ECM and examples of its applications
The main ECM idea is to represent the CFPs as a sum of two terms:
^fc Jjfc | ofc rn\
p = Bp,q + Bp,S. (2)
The first term Bpq describes the electrostatic interaction between an impurity ion and ions of crystal lattice, which, at this point, are considered as the point charges with the electrical charges equal to the formal valences of the ions. The second term BpS is proportional to the overlap integrals between the wave functions of an impurity ion and ligands. Its proper calculation implies the quantum-mechanical treatment of both central ion and ligands, to distinguish between the overlap integrals corresponding to the electron orbitals with different l and m quantum numbers. This term models the effects of the covalent bond formation and considerably improves the original point charge model of CF. As a result, good agreement between the ECM calculated and experimentally observed energy levels of the impurity ions can be obtained. The analytical expressions for calculating both terms are as follows [10]:
Bpkq = -Kp(1 - ,p)eV> £ qi, (3)
i Ri
Bp'S = Kpe2£ (GsS2(s) + S2(.) + GnS2(n)) ^^^. (4)
ii
The summation in both equations is performed over the crystal lattice ions denoted by index i with charges qi. The spherical coordinates of the i-th ion of crystal lattice in the coordinate system, whose origin coincides with the impurity ion position, are denoted as Ri,di,^i. The averaged values (rp> of p-th power of the impurity ion electrons coordinate can be found in the literature [11]. Alternatively, they can be calculated numerically using the corresponding ions' radial wave functions. The values of the numerical factors Kp and expressions for the polynomials Vp are given in Ref. [10]. The value of l in Eq. (4) is 2 (3) for the d-(f-) electrons, respectively. The shielding constant ap is required only for the f ions, to take into account
the screening effects of the filled 5s and 5p orbitals; it is identically zero for the d-electrons, since the d-shells are the outer ones. The overlap integrals between the d-functions of the central ion and p-, s-functions of the ligands are denoted by Ss, Sa, Sn. In the (1m|1'm!) notation: Ss = (d0|s0), SCT = (d0|p0), Sn = (d1|p1). These integrals can be calculated numerically; the main mathematical difficulty is that the wave functions involved are centered on different ions, and a coordinate transformation should be made to express all coordinates in one common system of reference [12]. The ECM dimensionless adjustable parameters are denoted by Gs, Ga, Gn. Their numerical values are determined by matching the positions of three lowest in energy experimentally observed spectral bands to the corresponding calculated energy levels. Very often they can be taken equal: Gs = Ga = Gn = G. In this case, the number of the ECM fitting parameters is reduced to one only, and only the lowest in energy absorption band is required to determine the value of G by adjusting the corresponding calculated energy level to its experimental counterpart. The numerical factors yp are y2 = 1, y4 = -4/3 for the d-electrons and y2 = 3/2, y4 = 1/3, Y6 = -3/2 for the f-electrons [10].
Further details of the ECM and description of its application can be readily found in the relevant literature. For the sake of brevity, we just mention the main applications of the ECM, supporting those examples by the corresponding references, where the reader can get more information.
i) The first and the most obvious application of the ECM is the calculation of the CFPs with subsequent diagonalization of the CF Hamiltonian and comparison of the calculated energy levels with the experimental absorption spectra of the corresponding systems. The coordinates of the crystal ions can be obtained from the structural data; the overlap integrals between the impurity ions and ligands can be calculated numerically for different interionic distances and extrapolated then by the exponential or linear decreasing functions of distance. Remarkably, even with one ECM fitting parameter, the agreement between the calculated and observed energy levels is very good. Energy levels of both d- and f-ions can be calculated with the ECM. Examples of ECM application to the d-ions can be given as follows: Li2MgSiO4:Cr4+ and Li2ZnSiO4:Cr4+ [13], MgAl2O4:Ni2+ [14], Mn4+ in a number of red-light emitting phosphor materials such as Mg2TiO4 [15], Ba2LaNbO6 [16], LaAlO3 [17] etc. Spectra of many crystals doped with f-ions or containing them as a part of their chemical formula and crystal lattice were analyzed by B.Z.Malkin et al: LiYF4:Tm3+ [18], LiYF4:Ho3+ [19], Pr4+ in BaPrO3 [20], Er3+ in (ErxY1-x)2BaNiO5 [21], Pr3+ in CsCdBr3 [22], Tb3+ in Tb2Ti2Oz [23] etc.
ii) As an additional result of the CF diagonalization, it is possible to get the wave functions of all CF split levels, which allows to make a few more steps when analyzing the spectroscopic properties of impurities in crystals. For example, magnetic susceptibilities of trivalent lanthanides in rare earth titanate pyrochlores in Ref. [24] and in NdFe3(BO3)4 [25] were calculated in the ECM framework. Intensities of the lanthanide 5d-4f transitions were successfully modeled in Refs. [26,27].
iii) Explicit expressions for calculating CFPs allow for finding their derivatives with respect to the Cartesian coordinates (or, finally, in terms of the normal modes), which immediately opens a way to study the phenomena related to the electron-vibrational interaction. Examples of such an analysis can be found in Refs. [28,29].
iv) Equations for the calculations of the CFPs can be used to get the CFPs numerical values for the decreased or increased interatomic distances, which allows to model the pressure effects
M.G. Brik, N.M. Avram, C.N. Avram
on the behavior of the calculated energy levels. In particular, it can be possible to find the dependence of the CF strength on the interatomic distance and estimate the bulk modulus of the host material, parameters of the electron-vibrational interaction, Stokes shift and the Jahn-Teller stabilization energy [30-32].
v) Calculations of the CFPs in the framework of the ECM allow to link the model with other CF models. Thus, the superposition model intrinsic parameters and power exponents for the Co2+ ions were extracted from the distance dependence of the CFPs calculated by application of the ECM equations in Ref. [33].
vi) Finally, the unified method of CFPs and energy levels calculations that is offered by the ECM allows for meaningful comparison of the CFPs sets and calculated energy levels of dopants for isostructural crystals and isoelectronic/isovalent impurity ions. The trends in variation of these properties are very important for smart search for new materials with desired characteristics.
3. Conclusions
Success of the ECM is determined by its firm physical grounds, rigorous derivation of all necessary equations and clear meaning of each term in all mathematical expressions. Wide applicability of the model, its ability to tackle complicated crystal systems and physical phenomena
make the ECM to be at the top of the CF models developed and used so far.
References
1. Bethe H. Ann. Physik 3, 133 (1929)
2. Van Vleck J. Phys. Rev. 41, 208 (1932)
3. Stevens K.W.H. Proc. Phys. Soc. A 65, 209 (1952)
4. Judd B.R. Proc. Royal Soc. A 232, 458 (1955); ibid 241, 414 (1957)
5. Liu G.K., in: Spectroscopic Properties of Rare Earth in Optical Materials, edited by G.K. Liu, B. Jacquier, p. 1, Springer, New York (2005)
6. Su P., Ma C.-G., Brik M.G., Srivastava A.M. Opt. Mater. 79, 129 (2018)
7. Newman D.J. Adv. Phys. 20, 197 (1971)
8. Newman D.J., Ng B. Rep. Prog. Phys. 52, 699 (1989)
9. Yeung Y.Y., in: Optical Properties of 3d-Ions in Crystals: Spectroscopy and Crystal Field, Analysis, edited by N.M. Avram, M.G. Brik, Ch. 3, p. 95, Springer and Tsinghua University Press (2013)
10. Malkin B.Z., Crystal Field and Electron-Phonon Interaction in Rare-Earth Ionic Paramagnets, in: Spectroscopy of Solids Containing Rare-Earth Ions, edited by A.A. Kaplyanskii, B.M. Macfarlane, p. 13, North Holland, Amsterdam (1987)
11. Abragam A.G., Bleaney B. Electron Paramagnetic Resonance of Transition Ions, Oxford University Press, Oxford, Clarendon (1970)
12. Brik M.G., Avram N.M. AIP Conf. Proc. 1131, 86 (2009)
13. Jousseaume C., Vivien D., Kahn-Harari A., Malkin B.Z. Opt. Mater. 24, 143 (2003)
14. Brik M.G., Avram N.M., Avram C.N., Rudowicz C., Yeung Y.Y., Gnutek P. J. Alloys Compds. 432, 61 (2007)
15. Medic M.M., Brik M.G., Drazic G., Antic Z., Lojpur V.M., Dramicanin M.D. J. Phys. Chem. C 119, 724 (2015)
16. Brik M.G., Srivastava A.M. J. Lumin. 132, 579 (2012)
17. Srivastava A.M., Brik M.G. Opt. Mater. 35, 1544 (2013)
18. Klimin S.A., Pytalev D.S., Popova M.N., Malkin B.Z., Vanyunin M.V., Korableva S.L. Phys. Rev. B 81, 045113 (2010)
19. Shakurov G.S., Vanyunin M.V., Malkin B.Z., Barbara B., Abdulsabirov R.Y., Korableva S.L. Appl. Magn. Res. 28, 251 (2005)
20. Popova M.N., Klimin S.A., Malkin B.Z., Kasatkina L.A., Cao G., Crow J. Phys. Lett. A 223, 308 (1996)
21. Popova M.N., Klimin S.A., Chukalina E.P., Malkin B.Z., Leviitin R.Z., Mill B.V., Antic-Fidancev E. Phys. Rev. B 68, 155103 (2003)
22. Popova M.N., Chukalina E.P., Malkin B.Z., Iskhakova A.I., Antic-Fidancev E., Porcher P., Chaminade J.P. Phys. Rev. B 63, 075103 (2001)
23. Klekovkina V.V., Malkin B.Z. Opt. Spectr. 116, 849 (2014)
24. Malkin B.Z., Lummen T.T.A., van Loosdrecht P.H.M., Dhalenne G., Zakirov A.R. J. Phys.: Condens. Matter 22, 276003 (2010)
25. Popova M.N., Chukalina E.P., Stanislavchuk T.N., Malkin B.Z., Zakirov A.R., Antic-Fidancev E., Popova E.A., Bezmaternykh L.N., Temerov V.L. Phys. Rev. B 75, 224435 (2007)
26. Kirm M., Stryganyuk G., Vielhauer S., Zimmerer G., Makhov V.N., Malkin B.Z., Solovyev O.V., Abdulsabirov R.Y., Korableva S.L. Phys. Rev. B 75, 075111 (2007)
27. Makhov V.N., Kirm M., Struganyuk G., Vielhauer S., Zimmerer G., Malkin B.Z., Solovyev O.V., Korableva S.L. J. Lumin. 132, 418 (2012)
28. Malkin B.Z., Solovyev O.V., Malishev A.Y., Saikin S.K. J. Lumin. 125, 175 (2007)
29. Brik M.G., Avram C.N. J. Lumin. 102-103, 283 (2003)
30. Brik M.G., Avram N.M. J. Phys.: Condens. Matter 21, 155502 (2009)
31. Avram N.M., Brik M.G., Kityk I.V. Opt. Mater. 32, 1668 (2010)
32. Brik M.G., Avram C.N. J. Lumin. 131, 2642 (2011)
33. Brik M.G., Yeung Y.Y. J. Phys. Chem. Solids 69, 2401 (2008)
