arxiv: v1 [cond-mat.mtrl-sci] 14 Oct 2016

Transcription

1 Magnetic resonances of multiferroic TbFe BO arxiv:11.55v1 [cond-mat.mtrl-sci] 1 Oct 1 Dávid Szaller, 1, Vilmos Kocsis, 1, Sándor Bordács, 1, Titusz Fehér, 1,, Toomas Rõõm, Urmas Nagel, Hans Engelkamp, 5 Kenya Ohgushi, and István Kézsmárki 1, 1 Department of Physics, Budapest University of Technology and Economics, 1111 Budapest, Hungary MTA-BME Lendület Magneto-optical Spectroscopy Research Group, 1111 Budapest, Hungary Condensed Matter Research Group of the Hungarian Academy of Sciences, 1111 Budapest, Hungary National Institute of Chemical Physics and Biophysics, Akadeemia tee, 118 Tallinn, Estonia 5 High Field Magnet Laboratory HFML-EMFL, Radboud University, Toernooiveld 7, 55 ED Nijmegen, The Netherlands Institute for Solid State Physics ISSP, University of Tokyo, Kashiwa-no-ha, Kashiwa, Chiba , Japan Dated: September 7, 18 Low-energy magnetic excitations of the easy-axis antiferromagnet TbFe BO are investigated by farinfrared absorption and reflection spectroscopy in high magnetic fields up to T. The observed field dependence of the resonance frequencies and the magnetization are reproduced by a mean-field spin model for magnetic fields applied both along and perpendicular to the easy axis. Based on this model we determined the full set of magnetic interactions, including Fe-Fe and Fe-Tb exchange interactions, single-ion anisotropy for Tb ions and g-factors, which describe the ground-state spin texture and the low-energy spin excitations of TbFe BO. Compared to earlier studies we allow a small canting of the nearly Ising-like Tb moments to achieve a quantitative agreement with the magnetic susceptibility measurements. The additional high energy magnetic resonance lines observed, besides the two resonances expected for a two-sublattice antiferromagnet, suggest a more complex six-sublattice magnetic ground state for TbFe BO. I. INTRODUCTION Magnetoelectric multiferroics, i.e. materials hosting both ferroelectric and ferromagnetic orders, attracted enormous interest due to their potential in information technology applications. 1 The magnetoelectric effect emerges not only in the static limit but also in the optical regime, as it is manifested in the difference between the refractive indices of counter-propagating light beams. 7 9 Indeed, strong directional dichroism, i.e. different absorption coefficient for light beams travelling in opposite directions, has been reported for spin excitations in multiferroics and was proposed as a new principle of directional light switch operating in the GHz-THz range. 7,1 1 Recently a new family of magnetoelectric multiferroic crystals, RFe BO rare-earth ferroborates, attracted much attention from the scientific community. Their uniqe crystal structure possessing magnetic iron and rare-earth sites R in a chiral arrangement allows the investigation of a wide variety of exotic magnetic and magnetoelectric phenomena. 17 The Fe + ions are surrounded by edge-sharing O octahedra and form quasi one-dimensional helical chains along the trigonal c axis of the crystal., These helices are expected to be only weakly connected by the Fe-O-O-Fe superexchange paths. 5 Thus the magnetic interaction between the iron chains and the rare-earth ions located between the iron helices plays an important role in tuning the effective dimensionality of the magnetic system. The single-ion anisotropy of the rare earth spins is transmitted to the otherwise nearly-isotropic Fe spins via J f d exchange interactions. The dominant magnetic interaction is the J dd exchange coupling between Fe spins, which leads to an antiferromagnetic ordering of the Fe subsystem. The rare earth spins largely separated from each other remain paramagnetic and are only polarized by the ordered Fe moments.,7 The strong spin-orbit coupling of the rare-earth ions plays a key role in the multiferroicity of these materials, i.e. their magneto-electric response is dominated by the rare-earth sites. 18,19,8 Correspondingly, the magnetoelectric properties of rare-earth ferroborates can be efficiently tuned by the selection of different rare earth elements characterized by different magnetic anisotropies. TbFe BO is a particularly interesting member of the rare-earth ferroborate family. Due to the effect of the strong crystal field, the ground-state doublet of Tb + ion is separated from the excited states by a considerable energy gap of 5 mev THz. 9, Thus at low temperatures the S Tb = spin of the Tb + ion behaves like an Ising moment pointing along the trigonal c axis of the crystal. As in the sister compounds R=Pr Ref. [1] and R=Dy Ref. [], in TbFe BO the easy-axis anisotropy of the rare-earth ion is transmitted to the antiferromagnetic S Fe = 5 iron system. Consequently, TbFe BO shows a collinear antiferromagnetic order below T N = K with all spins lying along the c axis.,9, Magnetic field along the c axis induces spin flop transition and all the Tb moments become parallel to the field while the sublattice magnetization of the antiferromagnetic Fe subsystem rotates to the ab plane and a weak canting develops along the c axis. The spin-flop transition field is B S F =.5 T at T = K and increases with increasing temperature. The magnetic order of TbFe BO was widely investigated by both magnetization and elastic neutron scattering experiments.,9 The temperature and field dependent behavior of the static magnetization was reproduced by former mean-field calculations 9,. However, the collective magnetic excitations of the ground state were only studied in zero magnetic field via optical spectroscopy. Here we investigate the low-energy magnetic excitations of TbFe BO using far-infrared optical spectroscopy up to high magnetic fields applied along and perpendicular to the trigonal c axis. The observed field dependence of the reso- Typeset by REVTEX

2 Magnetic Field [T] a B c = 1 cm -1 1B 1C 1A = 1 cm -1 5 b c = 1 cm -1 B c = 1 cm -1 R = = 1 cm Magnetic Field [T] Magnetic Field [T] FIG. 1. Color online Magnetic field dependence of the low-frequency magnetic resonances in TbFe BO at T = K. In panel a and c the magnetic field points along the trigonal c axis, while in panel b the field is perpendicular to the c axis. Light propagation was always parallel or antiparallel to the applied field Faraday geometry. Magnetic absorption and reflection spectra are vertically shifted in proportion to the magnetic field. In panel a, red and blue absorption curves measured in HFML correspond to light propagation parallel and antiparallel to the magnetic field, respectively. Black absorption curves on all of the panels and green reflectivity spectra on panel c were measured on the TeslaFIR setup. Insets of panels a and b show the weaker resonances of the corresponding spectra on a 1 times magnified absorption scale. Grey shaded lanes are guides for the eye. nance frequencies and the magnetization are reproduced by a mean-field spin model. In contrast to earlier studies in our model we allow a small canting of the quasi-ising Tb moments and a Tb-Fe exchange interaction is introduced for ions located in the same ab plane to achieve a quantitative agreement with the static and dynamic magnetic properties. Fourier transform spectroscopy was used to study the optical absorption and reflection of TbFe BO in theν =. THz frequency range with 8 GHz resolution. The magnetic field dependence of the spectra in the B= 17 T magnetic field range was investigated using the TeslaFIR setup of the National Institute of Chemical Physics and Biophysics in Tallinn. 7 Optical absorption experiments up to T were carried out in the High Field Magnet Laboratory in Nijmegen HFML. The spectra were measured in the Faraday configuration, i.e. in magnetic fields parallel to the direction of light propagation, using oriented single crystal samples with a typical thicknesses of 1 mm. The crystals were grown by a flux method.,8,5 II. EXPERIMENTAL DETAILS III. RESULTS AND DISCUSSION A. Experimental results In Figs. 1a-c, the optical absorption spectrum measured at T = K and in B = T shows a clear resonance at ν 1 =. THz. Though this resonance has already been observed and assigned as an antiferromagnetic resonance of the

3 Fe system, its field dependence has not been investigated so far. In magnetic fields B< B S F =.5 T parallel to the c axis the ν 1 resonance shows a V-shape splitting toν 1A andν 1B modes, as shown in Figs. 1a and 1c. Above B S F these resonances are replaced by a single mode,ν 1C, which hardens linearly with increasing field. These resonances are also visible in the field dependence of the reflectivity spectra whose inverse line shape dip in the reflectivity as opposed to dielectric resonances supports their magnetic nature. Two further resonances with field independent ν =.9 THz andν = 1.17 THz frequencies, best visible in Fig. 1a, also appear in the spin-flop phase. Despite their constant frequency, their oscillator strength grows with increasing field, indicating the magnetic origin of these modes. Besides these pronounced resonances some weaker ones can also be observed in the low-field phase, as shown in the inset of Fig. 1a. These resonance lines start atν =.98 THz and ν 5 = 1. THz frequencies in B= T and also show a V- shape splitting with increasing magnetic field. The splitting has the same slope as for theν 1 =. THz mode. The weakν andν 5 resonances, and the field independentν and ν modes are not visible in the reflectivity spectra of Fig. 1c. The absorption coefficient of magnetic excitations is the same for light propagation parallel and antiparallel to the magnetic field within the accuracy of the measurement. It means that despite the chiral crystal structure, the homochiral sample does not show considerable magneto-chiral dichroism MChD in the studied frequency window. 7,8,11 Static studies revealed that the magnetoelectric effect is mainly associated with Tb sites in TbFe BO. 8 Since static and dynamic magnetoelectric effects are closely related by a sum rule and the latters are responsible for the directional dichroism in the THz spectral range, the absence of MChD suggests that the observed excitations are of purely magnetic origin, and belong to the Fe subsystem. This is in accordance with the strong Ising character of the Tb moments, which act on the Fe spins as a static internal magnetic field but otherwise do not contribute to the spin dynamics. In contrast, a recent study 1 found strong directional dichroism in the easy plane antiferromagnet Sm.5 La.5 Fe BO, which originates from the coupled dynamics of the rare-earth sites with magnetoelectric character and the Fe spins. As shown in Fig. 1b, theν 1 resonance exhibits a quadratic shift towards higher frequencies with increasing magnetic field applied perpendicular to the c axis, as expected for an easy axis antiferromagnet. The other excitation discernible in this configuration is theν mode, whose frequency is nearly frequency independent in the range B = 1 T although it gains oscillator strength with increasing magnetic field. The temperature dependence of theν 1 resonance frequency in zero field shows the behavior expected for an anisotropic antiferromagnet, namely it gets softer and broader as the temperature approaches T N = K. Our results measured in transmission and reflection geometries are in good agreement with the previous study of A. M. Kuzmenko et al., as shown in Fig. c. B. Classical mean-field model Here we propose a classical mean-field model, which captures the magnetic field dependence of the ground state and that of the strong magnetic modes. We model the Fe and Tb moments, two sublattices for each, as classical vectors with different lengths, 5 for Fe S FeA and S FeB and for Tb moments S TbA and S TbB. The energy of a single magnetic unit cell can be written as: H=Λ S c + TbA S c TbB J dd S FeA S FeB J f d1 S FeA S TbB + S FeB S TbA J f d S FeA S TbA + S FeB S TbB B S FeA + S FeB g Tb µ B B S TbA + S TbB. 1 We include a strong negative uniaxial single-ion anisotropy termλfor the Tb sites to model their Ising-like nature. In case of the terms describing the exchange energy contributions, numeric prefactors correspond to the coordination numbers. An isotropic negative J dd exchange term connecting the Fe sublattices is responsible for the antiferromagnetic order. It is reasonable to consider the magnetic anisotropy of the rare-earth system only since in non-magnetic YFe BO and magnetically isotropic GdFe BO members of the crystal family the Fe spin system is nearly isotropic Heisenberg-like. 7 The Tb and Fe moments are coupled by the J f d1 and J f d exchange terms, where the dominant J f d1 connects moments on adjacent ab layers and the weaker J f d links ions on the same ab layers. This kind of coupling was neglected so far, 17,9, since the possible superexchange path corresponding to J f d is quite long. The last two terms of Eq. 1 describe the Zeeman energy, where B is the external magnetic field. 1. Estimation of model parameters based on magnetization data Approximate values for the parameters in Eq. 1 can be deducted from the static magnetization data and can be further tuned to fit the lowest excitation frequencies of the system. For the first estimation of the exchange and anisotropy constants we use the Landé values g Fe = and g Tb = 1.5 given for the free ions, which could be modified by crystal-field effects. The slope of the magnetization in B > B S F magnetic field parallel to the c axis is governed by the susceptibility of the Fe subsystem, because in the spin flop phase all the

4 Magnetization [ B /f.u.] 1 8 a c Magnetic field [T] 1 Temperature [K] b 5B A B 1A 1B 5A 1 1C Magnetic field [T] FIG.. Color online Comparison of experiments and model calculations. a The red and blue curves are experimental magnetization data reproduced from [] and [9], respectively, while black curves resulted from our classical Monte Carlo simulation. When the magnetic field B is perpendicular to the trigonal c axis the three lines coincide within the linewidth. b The magnetic field dependence of the low frequency Fe resonances; symbols are used for the experimental data and lines for the simulation results. Black color is used when the magnetic field is parallel and red when it is perpendicular to the trigonal axis. The nearly field-independent dashed mode was not observed in our Faraday geometry measurements light propagation parallel to the magnetic field, but is allowed in the Voigt geometry light propagation perpendicular to the magnetic field. Smaller symbols indicate the weaker resonances. c The temperature dependence of the zero field resonance, as seen in our optical experiment black squares, in backward-wave tube experiments gray circles, and as calculated solid black line from the temperature dependence of the magnetic moment lengths using Eq.. Tb moments point along the magnetic field and do not contribute to the susceptibility. The susceptibility per formula unitχ c =.1±.1 µ B T is determined by the J dd exchange, hence J dd g Fe χ =.9±. mev.,9 This value of c J dd is in good agreement with most of the previous studies J dd.1 mev,. mev, mev.,9, At the spin-flop transition the S TbB moments occuping one half of the Tb sites flip from S c TbB = to S c TbB =, giving rise to a g Tb jump in the magnetization per formula unit. In addition, the Fe moments show some canting in the spin-flop phase which also gives a minor contribution to the jump of the magnetization. The spin-flop transition takes place when the energy of the collinear and spin flop phases are equal, resulting in an approximate expression for the Tb-Fe coupling: J f d1 J f d / g Tbµ B B S F g Tbµ B B S F χ c, 9 where we neglected higher order terms inχ c. Due to the different exchange paths J f d / is expected to be much smaller than J f d1. Correspondingly all of the previous studies neglected the contribution of J f d and attributed the whole Tb-Fe coupling to the J f d1 exchange. 9, On the other hand, static magnetization data are not sufficient to unambiguously determine both J f d1 and J f d, thus, in this section we assume J f d = following previous works. The effective field along the c axis acting on the antiferromagnetic Fe system can be approximated by the sum of the external field and the effective field of the Tb moments. Since in the spin-flop phase both of the Tb sublattices are in the S c Tb =+ state, their effective field on the Fe site is J f d1 + J f d / J f d1 = g Tb 5g Fe B S F =.1 T. Thus, the M c = 9.1±.1µ B magnetization jump,9 at the spinflop transition is M c B S F =g Tb +χ c B S F+ J f d1 + J f d / a g Tb +χ c B S F 1+ g Tb, b 5g Fe giving rise to a refined value of g Tb 1.5±.5 which is significantly lower than the g Tb = 1.5 Landé value. Neutron scattering studies reported g Tb S c Tb = 8.5µ B ordered Tb moment at T= K, corresponding to g Tb = 1., which is in good accordance with our analysis.

5 5 TABLE I. Parameters of Eq. 1 obtained by using the static magnetization data,9 first row and refined parameters using both magnetization and magnetic resonance data second row. In the last column upper and lower error bounds are indicated by + and, respectively. g Fe g Tb J dd [mev] J f d1 [mev] J f d [mev] Λ[meV] 1.5±.5.9±..± ±.5.7±.15.5±..± Using Eq. the strength of the Tb-Fe exchange can be determined, J f d1.±.1 mev. This is in good agreement with the. mev and.9 mev values of previous magnetization studies. 17,9 The same coupling constant was determined from the splitting of the ground quasi-doublet of the Tb ions, which was observed as a splitting of the infrared transitions, corresponding to.5 mev. For fields perpendicular to the c axis, the susceptibility of the system isχ ab =.1±. µ B T, which is about % larger thanχ c.,9 As the Fe system is expected to be isotropic, the anisotropy of the susceptibility indicates the small canting of the Tb moments and thus can be used to estimate the anisotropy of the Tb sites: Λ g Tb =.1 χ ab χ c +1.. mev upper and lower error bounds are indicated by+and, respectively. The uncertainty of Λ comes from the variation of the experimental values forχ ab andχ c. However, due to the length of the Tb moments,λ S c Tb gives the dominant energy scale of the system in the studied magnetic field range. This justifies the approximation that Tb moments behave almost like Ising spins. The values obtained for Λ correspond to the range of the lowest excited crystal field energy levels calculated for the Tb + ion. The model parameter set obtained above is presented in the first row of Table I, which reproduces the static magnetization data. However, the static magnetization data only supports a rough estimation of the model parameters. Moreover, in the former expressions only one combination of the two types of Tb-Fe coupling appears, namely J f d1 J f d /, and therefore in studies based on the magnetization data the minor J f d was simply neglected. 9, In contrast, the magnetic field dependence of the dominant low-frequency magnetic excitations allows us to separate J f d1 and J f d and refine the values of all parameters in the Hamiltonian in Eq. 1.. Determination of model parameters based on magnetic resonances Assuming Ising-like Tb moments Λ, the zero temperature resonance frequencies of the Fe system can be calculated 8 using the S Tb = and S Fe = 5 values: ν 1 B== ν 1A/B B= < B S F = ν 1C B= > B S F = ν 1 B= B ab = J f d1 J f d J f d1 J f d J f d1 + J f d J f d1 J f d S Tb J dd S Fe J f d1 J f d S Tb S Tb J dd S Fe J f d1 J f d S Tb ± B 5 S Tb + B S Tb J dd S Fe J f d1 J f d S Tb + B. 7 Using Eq., the J dd Fe-Fe exchange can be determined with higher accuracy than our previous estimation from theχ c magnetic susceptibility. Based on the experimentalν 1 B==.±.5 THz frequency value we get J dd =.7±.15 mev. In the spin flop phase at = B S F the effective magnetic field acting on the antiferromagnetic Fe system is B S F + J f d1 + J f d /, giving information about the sum of J f d1 and J f d, thus can be used to unambiguously determine J f d. The experimental value of theν 1C resonance frequency in the flop phase can be extrapolated toν 1C = B S F =195 GHz, corresponding to an effective field of 7 T. This results in J f d1 =.5 mev and J f d =. mev, and refines the Tb g-factor to g Tb = 1.5 according to Eq. a. The J f d exchange is indeed weaker than J f d1 but does not have a ferromagnetic character, in contrast to former expectations based on the crystal and magnetic structure. Thus, in the zero-field ground state the bond corresponding to J f d is frustrated. According to Eqs. 5 and, the slope of theν 1A/B andν 1C modes yields the g-factor of the Fe system, which within the error of the measurement is equal to the spin-only g= value. Using Eq. 7 to fit the resonance frequencies measured in the B c case we get the same g-factor, thus the spin-only Fe g-factor is isotropic, as expected.

6 When considering the finite temperature excitations of the system, in the zero field case Eq. remains valid, only the temperature dependence of the lengths of the S Tb and S Fe ordered moments needs to be to be taken into account. However, a mean-field model like Eq. 1 is not able to properly describe the temperature dependence of the magnetic properties due to neglected thermal fluctuations, thus for a quantitative description additional experimental input is needed. In the elastic neutron scattering studies the lengths of the S Tb and S Fe ordered moments were reported in the whole temperature range of the antiferromagnetic phase. By substituting these temperature dependent ordered moments into Eq. the temperature dependence of the zero field resonance can be well reproduced with g Fe =, J dd =.7 mev, J f d1 =.5 mev and J f d =. mev, as shown in Fig. c. Earlier backwardwave oscillator spectroscopy studies reported the same temperature dependence. For finite values of the Tb single-ion anisotropy Λ, the analytical solution corresponding to Eqs. -7 is too complicated. Thus we calculated the field dependence of the zero temperature resonances numerically. We used a classical Monte Carlo approach to find the minimal energy configuration of the fourspin system, and determined the resonances by calculating the response to small perturbations. The Tb single-ion anisotropy was set to Λ g Tb χ ab + g Fe J dd = 8.1 mev, 8 and the other exchange parameters used in the simulation are listed in the second row of Table I. The calculated field dependence of the magnetization and antiferromagnetic resonance frequencies reproduce the experimental curves, as shown in Fig. a and b. With finiteλthe Tb moments are not static any more but oscillate with a zero-field resonance frequency ofν Tb =ΛS Tb = 5 mev, which agrees well with the frequency range of the lowest excited crystal field energy levels calculated for the Tb + ion. The field independentν andν resonances and the weak ν andν 5 modes cannot be explained by this simple classical four-sublattice mean-field spin model. Their presence shows that the Fe sites are crystallographically not equivalent, as is expected for the low temperature P 1 1 space group 5,9 of TbFe BO. Thus the proper description of the magnetic resonances is possible only with six magnetic Fe sublattices which are connected by the various, non-equivalent Fe-O-Fe intrachain and Fe-O-B-O-Fe interchain superexchange paths. Distinction between intrachain and interchain coupling would allow the tuning of the dimensionality of the system, thus the Monte Carlo approach could probably also reproduce the magnetic properties and resonance frequencies at finite temperatures. Nevertheless, due to the weak structural distortion from the room-temperature R structure to the low temperature P 1 1, the magnetic properties can be approximated by assuming crystallographically equivalent Fe sites. IV. SUMMARY In this study we have investigated the low frequency magnetic excitations of the multiferroic TbFe BO using farinfrared spectroscopy. We developed a classical mean-field spin model which quantitatively describes the main features in field dependence of the magnetization data,9 and that of the resonance frequencies with a minimal set of magnetic interactions including exchange couplings and single-ion anisotropy. Our far-infrared experiments also pointed out that the magnetic structure of TbFe BO is more complicated than previously expected. There are six inequivalent magnetic Fe sublattices, thus a more detailed neutron diffraction study is necessary to clarify the real magnetic ground state. V. ACKNOWLEDGMENT The authors are grateful for enlightening discussions with K. Penc and L. Mihály. This project was funded by Hungarian Research Funds OTKA K18918, PD 11175, K178 and Bolyai 55/1/11, by the institutional research funding IUT- of the Estonian Ministry of Education and Research and the European Regional Development Fund project TK1 TR and UN. We acknowledge the support of the HFML- RU/FOM, member of the European Magnetic Field Laboratory. D. Sz. was supported by the ÚNKP-1-/III. New National Excellence Program of the Ministry of Human Capacities. 1 A. J. Freeman and H. Schmid eds. Magnetoelectric interaction phenomena in crystals Gordon and Breach, London, M. Fiebig, J. Phys. D: Appl. Phys. 8, R1 5. W. Eerenstein, N. D. Mathur and J. F. Scott, Nature, 759. R. Ramesh and N. A. Spaldin, Nat. Mater., L. W. Martin, Y.-H. Chuc and R. Ramesh, Materials Science and Engineering R 8, S. M. Wu, Shane A. Cybart, D. Yi, James M. Parker, R. Ramesh and R. C. Dynes, Phys. Rev. Lett. 11, I. Kézsmárki, D. Szaller, S. Bordács, V. Kocsis, Y. Tokunaga, Y. Taguchi, H. Murakawa, Y. Tokura, H. Engelkamp, T. Rõõm and U. Nagel, Nat. Commun. 5, 1. 8 G. L. J. A. Rikken and E. Raupach, Nature 9, G. L. J. A. Rikken, C. Strohm and P. Wyder, Phys. Rev. Lett. 89, I. Kézsmárki, N. Kida, H. Murakawa, S. Bordács, Y. Onose, and Y. Tokura, Phys. Rev. Lett. 1, S. Bordács, I. Kézsmárki, D. Szaller, L. Demkó, N. Kida, H. Murakawa, Y. Onose, R. Shimano, T. Rõõm, U. Nagel, S. Miyahara, N. Furukawa and Y. Tokura, Nat. Phys. 8, D. Szaller, S. Bordács and I. Kézsmárki, Phys. Rev. B 87, Y. Takahashi, R. Shimano, Y. Kaneko, H. Murakawa and Y. Tokura, Nat. Phys. 8, 11 1.

7 7 1 Y. Takahashi, Y. Yamasaki, and Y. Tokura, Phys. Rev. Lett. 111, Y. Okamura, F. Kagawa, M. Mochizuki, M. Kubota, S. Seki, S. Ishiwata, M. Kawasaki, Y. Onose and Y. Tokura, Nat. Commun., A. M. Kuzmenko, V. Dziom, A. Shuvaev, Anna Pimenov, M. Schiebl, A. A. Mukhin, V. Yu. Ivanov, I. A. Gudim, L. N. Bezmaternykh and A. Pimenov, Phys. Rev. B 9, A. K. Zvezdin, A. M. Kadomtseva, Yu. F. Popov, G. P. Vorobev, A. P. Pyatakov, V. Yu. Ivanov, A. M. Kuzmenko, A. A. Mukhin, L. N. Bezmaternykh and I. A. Gudim, JETP 19, A. K. Zvezdin, S. S. Krotov, A. M. Kadomtseva, G. P. Vorobev, Yu. F. Popov, A. P. Pyatakov, L. N. Bezmaternykh and E. A. Popova, JETP Lett. 81, A. K. Zvezdin, G. P. Vorobev, A. M. Kadomtseva, Yu. F. Popov, A. P. Pyatakov, L. N. Bezmaternykh, A. V. Kuvardin and E. A. Popova, JETP Lett. 8, 59. A. M. Kadomtseva, A. K. Zvezdin, A. P. Pyatakov, A. V. Kuvardin, G. P. Vorobev, Yu. F. Popov and L. N. Bezmaternykh, JETP 15, A. M. Kadomtseva, Yu. F. Popov, G. P. Vorobev, A. P. Pyatakov, S. S. Krotov and K. I. Kamilov, Low Temp. Phys., U. Adem, L. Wang, D. Fausti, W. Schottenhamel, P. H. M. van Loosdrecht, A. Vasiliev, L. N. Bezmaternykh, B. Buchner, C. Hess and R. Klingeler, Phys. Rev. B 8, 1. J. C. Joubert, W. B. White and R. Roy, J. Appl. Cryst. 1, J. A. Campá, C. Cascales, E. Gutiérrez-Puebla, M. A. Monge, I. Rasines and C. Ru z-valero, Chem. Mater. 9, S. A. Klimin, D. Fausti, A. Meetsma, L. N. Bezmaternykh, P. H. M. van Loosdrecht and T. T. M. Palstra, Acta Cryst. B 1, C. Ritter, A. Balaev, A. Vorotynov, G. Petrakovskii, D. Velikanov, V. Temerov, and I. Gudim, J. Phys.: Condens. Matter 19, A. I. Pankrats, G. A. Petrakovskii, L. N. Bezmaternykh and O. A. Bayukov, JETP T. Kurumaji, K. Ohgushi, and Y. Tokura, Phys. Rev. B 89, E. A. Popova, D. V. Volkov, A. N. Vasiliev, A. A. Demidov, N. P. Kolmakova, I. A. Gudim, L. N. Bezmaternykh, N. Tristan, Yu. Skourski, B. Büchner, C. Hess, and R. Klingeler, Phys. Rev. B 75, 1 7. M. N. Popova, T. N. Stanislavchuk, B. Z. Malkin and L. N. Bezmaternykh, J. Phys.: Condens. Matter C. Ritter, A. Vorotynov, A. I. Pankrats, G. A. Petrakovskii, V. Temerov, I. Gudim and R. Szymczak, J. Phys.: Condens. Matter 1. C. Ritter, A. I. Pankrats, I. Gudim and A Vorotynov, J. Phys.: Conf. Ser Y. Hinatsu, Y. Doi, K. Ito, M. Wakeshima and A. Alemi, J. Solid State Chem A. M. Kuzmenko, A. A. Mukhin, V. Yu. Ivanov, A. M. Kadomtseva, S. P. Lebedeva, and L. N. Bezmaternykh, JETP 11, S. Hayashida, M. Soda, S. Itoh, T. Yokoo, K. Ohgushi, D. Kawana, H. M. Rønnow and T. Masuda Phys. Rev. B 9, D. Szaller, S. Bordács, V. Kocsis, T. Rõõm, U. Nagel and I. Kézsmárki, Phys. Rev. B 89, A. I. Pankrats, G. A. PetrakovskiI, L. N. Bezmaternykh and V. L. Temerov, Phys. Solid State E. A. Turov, Physical properties of magnetically ordered crystals, Moscow: Izdat. Acad. Sci. SSSR, D. Fausti, A. A. Nugroho, P. H. M. van Loosdrecht, S. A. Klimin, M. N. Popova, and L. N. Bezmaternykh, Phys. Rev. B 7,.