pith. sign in

arxiv: 2604.02193 · v1 · submitted 2026-04-02 · ❄️ cond-mat.mtrl-sci

Terahertz optical activity near crystal field transitions of Tm3+ ions in magnetoelectric alumoborates

Pith reviewed 2026-05-13 20:36 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords terahertz spectroscopyoptical activitycrystal field transitionsTm3+ ionsmagnetoelectric materialsalumoboratesdynamic magnetoelectric susceptibilitypolarization rotation
0
0 comments X

The pith

Strong natural optical activity near crystal field transitions of Tm3+ ions produces polarization plane rotations up to 25 degrees in terahertz spectra of magnetoelectric alumoborates.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines crystal field excitations within the ground multiplet of Tm3+ ions through terahertz transmission spectra measured on TmAl3(BO3)4 and a lightly doped variant. These excitations are identified as mainly magnetic dipole transitions from the ground singlet A1 to the excited doublet E split by D3 symmetry. Strong natural optical activity appears near the transitions and rotates the polarization plane by as much as 25 degrees. The activity is accounted for by the combined magnetic and electric dipole contributions to the dynamic magnetoelectric susceptibility, with the observed fine structure tied to differing local distortions of the D3 crystal field in pure versus doped crystals.

Core claim

The central claim is that strong natural optical activity near the crystal field transitions of Tm3+ ions arises from magnetic and electric dipole transitions that contribute to the dynamic magnetoelectric susceptibility, with the fine structure of the modes reflecting different local distortions of the D3 symmetry that are resolved at low temperatures and differ between pure and doped samples.

What carries the argument

Dynamic magnetoelectric susceptibility formed by magnetic and electric dipole transitions between crystal-field levels of Tm3+ ions split under D3 symmetry.

If this is right

  • Polarization rotation reaches 25 degrees near the identified transitions.
  • The mode fine structure differs between pure Tm borates and lightly doped samples because of varying local crystal-field distortions.
  • The transitions are classified as primarily magnetic dipole from the A1 singlet to the E doublet.
  • Quantitative modeling of the spectra matches experiment once the local distortions and both dipole types are included in the susceptibility.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar polarization effects could appear in other rare-earth magnetoelectric crystals at terahertz frequencies when comparable crystal-field transitions are present.
  • The rotation magnitude suggests possible use for compact terahertz polarization elements if the material response can be engineered or scaled.
  • Temperature-dependent studies would test how thermal population of levels alters the strength of the magnetoelectric contribution.

Load-bearing premise

The fine structure of the modes is produced by different local distortions of the D3 crystal field, and the transitions remain predominantly magnetic dipole from the A1 ground singlet to the E doublet.

What would settle it

A terahertz measurement on a non-magnetoelectric analog crystal or on the same samples at temperatures high enough to thermally quench the transitions and eliminate the polarization rotation would show whether the claimed link to these specific crystal-field transitions holds.

Figures

Figures reproduced from arXiv: 2604.02193 by A.A. Mukhin, A.M. Kuzmenko, A. Pimenov, A. Shuvaev, I. A. Gudim, K.N. Boldyrev, S.V. Garnov, V.Yu. Ivanov.

Figure 1
Figure 1. Figure 1: Trigonal crystal structure (R32) of rare-earth aluminum borate projected onto the ab-plane (a) and b*c-plane (b) (b* is perpendicular to the a and c axes). (c) The orientation of the crystallographic axes with respect to natural crystal faces in the source crystal used to prepare a- and c-cut samples. The present work is devoted to terahertz study of a natural optical activity at low-energy Tm3+ crystal fi… view at source ↗
Figure 2
Figure 2. Figure 2: Scheme of the optical activity measurements us [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Transmission spectra (left) and polariza [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Transmission spectra (left) and polarization rotation angle (right) of the c-cut Tm0.05Yb0.1Y0.85Al3(BO3)4 sample (d=3.547 mm) for ac magnetic field hc-axis. Points – experiment, solid lines – theory as described in the text. For the present geometry with k||c-axis (c-cut) the relevant components of magnetic, dielectric and magnetoelectric permittivities, µxx,yy, εxx,yy, xx,yy can be represented as super… view at source ↗
Figure 5
Figure 5. Figure 5: Local arrangement of six possible positions of strongly distorted Tm3+ ions around a single Bi￾impurity for nearest neighbors (site #3) and next nearest neighbors (site #2). All remaining Tm³⁺ ions occupy sites #1, with larger Bi–Tm separations than at sites #2 and #3 (not shown) [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Splitting of the crystal field states of Tm3+ ions due to distortions in the vicinity of a Bi-impurity for different sites in pure TmAl3(BO3)4 (right) and simulated continuous distribution of the E-doublet splitting in diluted Tm0.05Yb0.1Y0.85Al3(BO3)4 (middle). The violation of the local symmetry due to distortions near Bi3+ impurities leads to the splitting of the E-doublet and to changes in the wave fun… view at source ↗
read the original abstract

Crystal field (CF) excitations in the ground multiplet $^3H_6$ of Tm$^{3+}$ ions were investigated using terahertz transmission spectra of magnetoelectric TmAl$_3$(BO$_3$)$_4$ and Tm$_{0.05}$Yb$_{0.1}$Y$_{0.85}$Al$_3$(BO$_3$)$_4$. These excitations were identified as mainly magnetic dipole transitions from the ground singlet A$_1$ to the next excited doublet E, split by the crystal field of the D$_3$ symmetry. The fine structure of the modes was resolved at low temperatures. It manifested differently in lightly doped and in pure Tm borates, consistent with different distortions of the local crystal field with the D$_3$ symmetry. Strong natural optical activity was observed near the CF transitions resulting in a polarization plane rotation up to 25 degrees. The optical activity is quantitatively described by contributions of magnetic and electric dipole transitions to dynamic magnetoelectric susceptibility and taking into account the classification of local distortions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript reports terahertz transmission spectroscopy on pure TmAl3(BO3)4 and lightly doped Tm0.05Yb0.1Y0.85Al3(BO3)4, identifying crystal-field excitations within the 3H6 multiplet of Tm3+ ions as predominantly magnetic-dipole transitions from the ground A1 singlet to the excited E doublet under D3 local symmetry. Fine structure is resolved at low temperature and attributed to differing local distortions of the D3 environment in the pure versus doped compounds. Strong natural optical activity is observed near these transitions, producing polarization-plane rotation up to 25°, and is modeled quantitatively via contributions of magnetic and electric dipole transitions to the dynamic magnetoelectric susceptibility.

Significance. If the local-distortion assignment is confirmed, the work supplies a concrete experimental demonstration of large THz optical activity in a magnetoelectric rare-earth borate together with a symmetry-based susceptibility framework that reproduces the rotation magnitude. This strengthens the connection between crystal-field fine structure and magnetoelectric response and may guide design of materials for THz polarization manipulation. The raw observation of 25° rotation and its temperature dependence constitute the clearest strength.

major comments (2)
  1. [Discussion of fine structure and susceptibility modeling] The quantitative reproduction of the 25° rotation rests on the assignment of the observed fine structure to distinct D3-symmetry local distortions that differ between the pure and doped crystals. The manuscript provides no explicit calculation of the expected splitting magnitudes, no comparison of the resulting magnetoelectric susceptibility tensors, and no test against alternative origins such as inter-ion coupling or defect-induced symmetry lowering. Without this, the claim that the susceptibility model is predictive rather than post-hoc remains unverified.
  2. [Quantitative modeling section] The abstract and modeling section state that the optical activity is 'quantitatively described' by magnetic- and electric-dipole contributions, yet no fitting parameters, covariance matrix, or direct overlay of measured versus calculated rotation spectra are supplied. Inclusion of these would allow assessment of whether the model is over- or under-constrained.
minor comments (2)
  1. [Experimental results] Error bars are absent from the reported rotation angles and transmission spectra; their addition would clarify the statistical significance of the 25° maximum and the fine-structure splittings.
  2. [Symmetry analysis] The classification of local distortions is invoked repeatedly but never tabulated or illustrated; a short supplementary table listing the distinct D3 subgroups and their expected transition intensities would improve clarity.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the positive assessment of our work and the constructive comments on the fine-structure assignment and modeling details. We address each major comment below and have revised the manuscript to improve clarity and transparency where possible.

read point-by-point responses
  1. Referee: The quantitative reproduction of the 25° rotation rests on the assignment of the observed fine structure to distinct D3-symmetry local distortions that differ between the pure and doped crystals. The manuscript provides no explicit calculation of the expected splitting magnitudes, no comparison of the resulting magnetoelectric susceptibility tensors, and no test against alternative origins such as inter-ion coupling or defect-induced symmetry lowering. Without this, the claim that the susceptibility model is predictive rather than post-hoc remains unverified.

    Authors: The assignment relies on the observed differences in fine structure between the pure and doped samples together with their temperature dependence, which match the expected behavior for D3 local environments under varying strain. We have added a dedicated paragraph discussing why inter-ion coupling is improbable at the low Tm concentration and why defect-induced symmetry lowering does not reproduce the observed splitting pattern. A qualitative comparison of how different D3 distortion parameters modify the magnetoelectric susceptibility tensor has also been included. However, quantitative ab-initio evaluation of the splitting energies lies outside the present experimental scope and is noted as a limitation for future work. revision: partial

  2. Referee: The abstract and modeling section state that the optical activity is 'quantitatively described' by magnetic- and electric-dipole contributions, yet no fitting parameters, covariance matrix, or direct overlay of measured versus calculated rotation spectra are supplied. Inclusion of these would allow assessment of whether the model is over- or under-constrained.

    Authors: We agree that greater transparency is needed. The revised manuscript now lists the numerical values of the dynamic magnetoelectric susceptibility components employed in the calculation. A new supplementary figure overlays the measured and calculated polarization-rotation spectra for both compounds, and estimated uncertainties on the key parameters are provided. These additions permit direct evaluation of the model’s constraints and fit quality. revision: yes

standing simulated objections not resolved
  • Explicit numerical calculation of the expected crystal-field splitting magnitudes arising from specific local D3 distortions.

Circularity Check

0 steps flagged

No load-bearing circularity; observations and model rest on independent dipole theory and direct spectra

full rationale

The derivation chain begins with measured THz transmission and polarization rotation data (up to 25°). These are interpreted via standard magnetic/electric dipole contributions to dynamic magnetoelectric susceptibility under D3 symmetry classification of local distortions. No equations reduce a fitted parameter to a renamed prediction, no self-citation supplies a uniqueness theorem that forces the result, and the quantitative description does not collapse to the input spectra by construction. Minor self-citations to prior dipole formalism are present but non-load-bearing. The skeptic concern about distortion assignment affects interpretive correctness, not circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The description rests on standard crystal-field symmetry for D3 point group and the decomposition of magnetoelectric susceptibility into electric and magnetic dipole channels; no new free parameters or invented entities are introduced beyond conventional fitting of transition strengths.

axioms (2)
  • domain assumption Crystal field of D3 symmetry splits the Tm3+ ground multiplet into A1 singlet and E doublet.
    Invoked to assign the observed modes as magnetic dipole transitions from ground A1 to excited E.
  • standard math Dynamic magnetoelectric susceptibility receives additive contributions from magnetic and electric dipole transitions.
    Standard framework used to quantitatively describe the observed optical activity.

pith-pipeline@v0.9.0 · 5529 in / 1397 out tokens · 38752 ms · 2026-05-13T20:36:09.110725+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

36 extracted references · 36 canonical work pages

  1. [1]

    R. P. Chaudhury, F. Yen, B. Lorenz, Y. Y. Sun, L. N. Bezmaternykh, V. L. Temerov, and C. W. Chu, Magnetoelectric Effect and Spont aneous Polarization in HoFe 3(BO3)4 and Ho0.5Nd0.5Fe3(BO3)4, Phys. Rev. B 80, 104424 (2009)

  2. [2]

    Liang, R

    K.-C. Liang, R. P. Chaudhury, B. Lorenz, Y. Y. Sun, L. N. Bezmaternykh, V. L. Temerov, and C.W. Chu, Giant magnetoelectric effect in HoAl3(BO3)4, Phys. Rev. B 83, 180417 (2011)

  3. [3]

    R. P. Chaudhury, B. Lorenz, Y. Y. Sun, L. N. Bezmaternykh, V. L. Temerov, and C. W. Chu, Magnetoelectricity and magnetostriction due to the rare-earth moment in TmAl 3(BO3)4, Phys. Rev. B 81, 220402(R) (2010)

  4. [4]

    A. M. Kadomtseva, Yu. F. Popov, G. P. Vorob’ev, N. V. Kostyuchenko, A.I. Popov, A. A. Mukhin, V. Yu. Ivanov, L. N. Bezmaternykh, I. A. Gudim, V. L. Temerov, A. P. Pyatakov and A. K. Zvezdin, High-temperature magnetoelectricity of terbium aluminum borate: The role of excited states of the rare- earth ion, Phys. Rev. B 89, 014418 (2014)

  5. [5]

    Ivanov, A.M.Kuzmenko, & A.A

    V.Yu. Ivanov, A.M.Kuzmenko, & A.A. Mukhin, Magne toelectric effect in yt terbium aluminum borate YbAl3(BO3)4 JETP Lett. 105, 435 (2017)

  6. [6]

    Ballman, A new series of synthetic borat es isostructural with th e carbonate mineral huntite American Mineralogist 47, 1380 (1962)

    A.A. Ballman, A new series of synthetic borat es isostructural with th e carbonate mineral huntite American Mineralogist 47, 1380 (1962)

  7. [7]

    A. I. Popov, D. I. Plokhov, and A. K. Zvezdin, Quantum theory of magnetoel ectricity in rare-earth multiferroics: Nd, Sm, and Eu ferroborates, Phys. Rev. B 87, 024413 (2013)

  8. [8]

    N. V. Kostyuchenko, V. Yu. Ivanov, A. A. Mukhin, A. I. Popov, and A. K. Zvezdin, Features of Magnetic and Magnetoelectric Properties in Rare-Earth Aluminum Borates RAl3(BO3)4 , IEEE Trans. Magn. 50, 2504504 (2014)

  9. [9]

    Phys.: Condens

    T Arima, Magneto-electric optics in non-centrosymmetric ferromagnets, J. Phys.: Condens. Matter 20, 434211 (2008)

  10. [10]

    Shimada, H

    Y. Shimada, H. Kiyama, and Y. Tokura, Magnetoelectric emission in rare-earth doped ferroelectric crystals La2Ti2O7 :R3+(R =Er, Eu, and Nd), Phys. Rev. B 75, 245125 (2007)

  11. [11]

    Yokosuk, Heung-Sik Kim, Kendall D

    Michael O. Yokosuk, Heung-Sik Kim, Kendall D. Hughey, Jaewook Kim, Andreas V. Stier, Kenneth R. O’Neal, Junjie Yang,Scott A. Crooker, Kristjan Haule , Sang-Wook Cheong, David Vanderbilt and Janice L. Musfeldt, Nonreciprocal directional dichroism of a chiral magnet in the visible range, npj Quantum Materials 5:20 (2020)

  12. [12]

    Malakhovskii, V.V

    A.V. Malakhovskii, V.V. Sokolov, I.A. Gudim, Natural optical activity of f-f transitions in ErAl3(BO3)4 single crystal, Chemical Physics 493, 102 (2017)

  13. [13]

    Bordács, I

    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, Chirality of matter shows up via spin excitations, Nature Physics, 8, 734 (2012)

  14. [14]

    Kocsis , S

    V. Kocsis , S. Bordács, Y. Tokunaga, J. Viirok, L. Peedu , T. Rõõm , U. Nage, Y. Taguchi, Y. Tokura, and I. Kézsmárki, Magnetoelectric spectroscopy of spin excitations in LiCoPO4, Phys. Rev. B. 100, 155124 (2019)

  15. [15]

    Takahashi and Y

    Y. Takahashi and Y. Tokura, Electromagnons with Giant Nonreciprocal Responses in Spin-Origin Multiferroics, J. Phys. Soc. Jpn. 94, 052001 (2025)

  16. [16]

    Okamura, F

    Y. Okamura, F. Kagawa, M. Mochizuki, M. Kubota, S. Seki, S. Ishiwata, M. Kawasaki, Y. Onose, and Y. Tokura, Microwave magnetoelectric effect via skyrmion resonance modes in a helimagnetic multiferroic, Nat. Commun. 4, 2391 (2013)

  17. [17]

    Okamura, F

    Y. Okamura, F. Kagawa, S. Seki, M. Kubota, M. Kawasaki, and Y.Tokura, Phys. Rev. Lett. 114, 197202 (2015)

  18. [18]

    Yoichi Nii, Ryo Sasaki, Yusuke Iguchi, and Yoshinori Onose, Microwave Magnetochiral Effect in the Non-centrosymmetric Magnet CuB2O4, Journal of the Physical Society of Japan 86, 024707 (2017)

  19. [19]

    A. M. Kuzmenko, A. Shuvaev, V. Dziom, Anna Pimenov, M.Schiebl, A. A. Mukhin, V. Yu. Ivanov, L. N. Bezmaternykh, and A. Pimenov, Giant gigahertz optical activity in multiferroic ferroborate, Phys. Rev. B 89, 174407 (2014). 10

  20. [20]

    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, Large directional optical anisotropy in multiferroic ferroborate, Phys. Rev. B 92, 184409 (2015)

  21. [21]

    Mukhin, A.M

    A.A. Mukhin, A.M. Kuzmenko, V.Yu. Ivanov, A.G. Pimenov, A.M. Shuvaev, V.E. Dziom, Dynamic magnetoelectric phenomena within electromagnons in rare-earth borate multiferroics, Physics-Uspehi 58, 993 (2015)

  22. [22]

    Pimenov, A

    A. Pimenov, A. A. Mukhin, V. Yu. Ivanov, V. D. Travkin, A. M. Balbashov, and A. Loidl, Possible evidence for electromagnons in multiferroic manganites, Nat. Phys. 2, 97 (2006)

  23. [23]

    A. M. Kuzmenko, V. Dziom, A. Shuvaev, Anna Pimenov, D. Szaller, A. A. Mukhin, V. Yu. Ivanov, and A. Pimenov, Sign change of polarization rotation under time or space inversion in magnetoelectric YbAl3(BO3)4, Phys. Rev. B 99, 224417 (2019)

  24. [24]

    Kebaili, M

    I. Kebaili, M. Dammak, E. Cavalli, M. Bettinelli, Energy levels and crystal-field analysis of Tm 3+ in YAl3(BO3)4 crystals, Journal of Luminescence 131, 2010–2015 (2011)

  25. [25]

    A. M. Kuzmenko, A. A. Mukhin, V. Yu. Ivanov, G. A. Komandin, A. Shuvaev, A. Pimenov, V. Dziom, L. N. Bezmaternykh, and I. A. Gudim, Terahertz spectroscopy of crystal-field transitions in magnetoelectric TmAl3(BO3)4, Phys. Rev. B 94, 174419 (2016)

  26. [26]

    V. L. Temerov, A. É. Sokolov, A. L. Sukhachev, A. F. Bovina, I. S. Édel’man, A. V. Malakhovskiĭ, Optical properties of trigonal single crystals (Yb,Tm)Al3(BO3)4 grown from fluxes based on the bismuth and lithium molybdates, Crystallog. Rep. 53, 1157 (2008)

  27. [27]

    M. N. Popova, K. N. Boldyrev, P. O. Petit, B. Viana, and L. N. Bezmaternykh, High-resolution spectroscopy of YbAl₃(BO₃)₄ stoichiometric nonlinear laser crystals, J. Phys.: Condens. Matter 20, 455210 (2008)

  28. [28]

    K. N. Boldyrev, M. N. Popova, L. N. Bezmaternykh, and M. Bettinelli, Nonequivalent Yb³⁺ centers in Y₁₋ₓYbₓAl₃(BO₃)₄ laser crystals, Quantum Electron. 41, 120–124 (2011)

  29. [29]

    K. N. Boldyrev, M. N. Popova, M. Bettinelli, V. L. Temerov, I. A. Gudim, L. N. Bezmaternykh, P. Loiseau, G. Aka, and N. I. Leonyuk, Quality of the rare-earth aluminium borate crystals for laser applications probed by high-resolution spectroscopy of the Yb³⁺ ion, Opt. Mater. 34, 1885–1889 (2012)

  30. [30]

    Oreshonkov, Evgenii M

    Aleksandr S. Oreshonkov, Evgenii M. Roginskii, Ni kolai P. Shestakov, Irina A. Gudim, Vladislav L. Temerov, Ivan V. Nemtsev, Maxim S. Molokeev, Sergey V. Adi htchev, Alexey M. Pugachev and Yuriy G. Denisenko, Structural, Elect ronic and Vibrationa l Properties of YAl 3(BO3)4 , Materials 13, 545 (2020)

  31. [31]

    Cohere nt Source Submillimeter Wave Spectroscopy

    G. V. Kozlov and A. A. Volkov, “Cohere nt Source Submillimeter Wave Spectroscopy”, in “Millimeter and Submillimeter Wave Spectroscopy of Solids”, Springer, Berlin (1998), p.51

  32. [32]

    Klimin, D

    S.A. Klimin, D. S. Pytalev, M. N. Popova, B. Z. Malkin, M.n. Vanyuin, S.L. Korableva , High- resolution optical spectroscopy of Tm3+ ions in LiYF4: Crystal-field energies, hyperfine and deformation splittings, and the isotopic structure Phys. Rev. B 81, 045113 (2010)

  33. [33]

    B. Z. Malkin, D. S. Pytalev, M. N. Popova, E. I. Baibekov, M. L. Falin, K. I. Gerasimov, and N. M. Khaidukov, Random lattice deformations in rare-earth-doped cubic hexafluoroelpasolites: High- resolution optical spectroscopy and theoretical studies, Phys. Rev. B 86, 134110 (2012)

  34. [34]

    A Ivanov, V

    M. A Ivanov, V. Y. Mitrofanov, L. D. Falkovskaya, and A. Y. Fishman, Random fields in disordered magnetics with Jahn-Teller ions, J. Magn. Magn. Mater. 36, 26 (1983)

  35. [35]

    Arteaga, Optical activity of oriented molecular systems in terms of the magnetoelectric tensor of gyrotropy, J

    O. Arteaga, Optical activity of oriented molecular systems in terms of the magnetoelectric tensor of gyrotropy, J. Opt. 16 , 125707 (2014)

  36. [36]

    M. N. Popova, B. Z. Malkin, K. N. Boldyrev, T. N. Stanislavchuk, D. A. Erofeev, V. L. Temerov, and I. A. Gudim, Evidence for a collinear easy-plane magnetic structure of multiferroic EuFe₃(BO₃)₄: Spectroscopic and theoretical studies, Phys. Rev. B 94, 184418 (2016)