Spin dynamics in linear magnetoelectric material Mn3Ta2O8
Abstract
We performed inelastic neutron scattering experiments on single crystal samples of a linear magnetoelectric material Mn3Ta2O8, which exhibits a collinear antiferromagnetic order, to reveal the spin dynamics. Numerous modes observed in the neutron spectra were reasonably reproduced by linear spin-wave theory on the basis of the spin Hamiltonian including eight Heisenberg interactions and an easy-plane type single-ion anisotropy. The presence of strong frustration was found in the identified spin Hamiltonian.
-
Apr 2024
Keywords: Inelastic neutron scattering, linear magnetoelectric effect, frustration
1 Introduction
![Refer to caption](x1.png)
The interaction between magnetism and electricity induces a magnetoelectric (ME) effect, in which the electric dipole moment and magnetic moment are simultaneously induced. This phenomenon was first observed in Cr2O3[1, 2], characterized by an electric polarization and a magnetic moment that respond linearly to applied magnetic and electric fields, respectively. To exhibit the linear ME effect, it is necessary to break both spatial inversion and time-reversal symmetries. Since the discovery in Cr2O3, the effect has been studied in various materials[3, 4, 5, 6], becoming one of the most active research topics in the field of condensed matter physics in recent years.
In frustrated magnets exhibiting noncollinear orders such as cycloidal structures[7, 8, 9, 10, 11], spontaneous electric polarization and ordered magnetic moments are simultaneously induced via the spin-current mechanism[12]. Other contributing mechanisms include exchange striction[13, 14] and spin-dependent metal-ligand hybridization[15, 16, 17]. In MnO3 ( = rare earth), numerous inelastic neutron scattering (INS) experiments have been conducted, revealing the presence of the frustration [18, 19, 20].
Mn3Ta2O8 is a tetragonal crystal with space group , where Mn2+ ions carry spin = [21, 22, 23, 24]. Figure 1 illustrates the crystal and magnetic structures, projecting the Mn2+ ions onto the plane, with different Mn sites indicated by different colors. Mn2+ ions at Mn1 sites represented by yellow symbols form spirals along the -axis, with the black arrows indicating the chiralities of the spirals. At = 24 K, a magnetic phase transition to a collinear antiferromagnetic structure occurs. The moment sizes of Mn2+ ions at Mn1, Mn2, and Mn3 sites are 3.91(3), 4.17(8) and 4.04(8) , respectively[24]. Previous studies [22, 24] have revealed that the Weiss temperature() is 120 K, suggesting strong frustration, as is five times larger than . The magnetic structure belongs to the space group, breaking both spatial inversion and time-reversal symmetries, leading to the observation of linear magnetoelectric effect. An antiferroelectric-like phase transition due to the off-center displacement of magnetic Mn2+ ions has been observed [24].
In the present study, INS experiments were performed using single crystal samples to elucidate the spin dynamics in Mn3Ta2O8. The observed spectra were reasonably explained by the spin Hamiltonian composed of eight Heisenberg terms and an easy-plane single-ion anisotropy term. Strong frustration was found between the intrachain antiferromagnetic interaction and interchain ferromagnetic interaction.
2 Experimental details
Single crystals of Mn3Ta2O8 were grown by the floating-zone method[24]. The quality of the samples was evaluated using High Energy Resolution (HER) triple-axis spectrometer installed at the C11 beam port at Japan Research Reactor 3 (JRR-3).
To collect INS spectra in wide momentum () - energy () space, INS experiments were carried out by High Resolution Chopper spectrometer (HRC) [25] installed at BL12 in J-PARC/MLF. A single crystal of Mn3Ta2O8 with a mass of 2.20 g was used. Temperatures were controlled using a Gifford-McMahon (GM) cryostat. Frequency of Fermi chopper was 100 Hz, and incident neutron energy () was 15 meV. The instrumental energy resolution at = 0 is 0.49 meV. The data reduction was performed by HANA software [26].
To detail the spectra in low-energy range, INS experiments were carried out by cold neutron triple-axis spectrometer (CTAX) installed at High Flux Isotope Reactor (HFIR) in Oak Ridge National Laboratory (ORNL). A single crystal with a mass of 1.57 g was used. Cooling was achieved with an ORANGE cryostat utilizing liquid helium, maintaining the measurement temperature at 1.5 K. The collimator setup was configured as guide - open - radial collimator - 120’. Final neutron energy was fixed at meV. BeO filter was employed downstream of the analyzer to cut off the higher harmonics. The instrumental energy resolution ranges from 0.22 meV to 0.40 meV, depending on the energy transfer.
3 Experimental results
![Refer to caption](x5.png)
![Refer to caption](x6.png)
First, we conducted measurements over a wide range of - space using the HRC spectrometer. In Figs 2(a) and 2(b), gapless excitations were observed, arising from in the range of 3 meV. Other excitations were observed in the range of 3 meV 5 meV, and flat excitations were observed near 8 meV. Note that the scale of the intensity is logarithmic. In Figs. 2(c) and 2(d), low energy excitations are observed in the range of meV. Further observations showed modes in the range of 3 meV 8 meV. The solid curves in Figs. 2(a)-2(d) represent the dispersion relations in each direction, calculated using the linear spin-wave theory (LSWT), the details of which will be described later.
Next, we conducted measurements using CTAX, which provides a higher energy resolution, to elucidate the details of low-energy excitations. The constant scans were performed centered at (3, 1, 0), along and as shown in Figs 3(a) and 3(b), respectively. The results are consistent with those measured by HRC spectrometer in Figs. 2(a) and 2(b). Solid curves represent the fitting curves obtained through multiple Gaussian fittings. The full widths at half maximum (FWHM) of most of the Gaussians are broader than the instrumental resolution. Due to the large number of modes and some having weak intensities, it is not feasible to separate all of them for fitting.
4 Analysis
![Refer to caption](x7.png)
The observed INS spectra were analyzed based on LSWT by using the SpinW package[27]. The spin Hamiltonian is as follows;
(1) | |||||
The magnetic structure used for the analysis is depicted in Fig. 1. Even though Mn2+ ions at Mn1, Mn2, and Mn3 sites have different moment sizes, we assumed a uniform moment size of 4 for all Mn2+ ions. Each interaction is defined in Fig. 4. For convenience, we designated the structure composed of Mn2+ ions at Mn1 site as “chain” and the structure composed of Mn2+ ions at Mn2 and Mn3 sites as “square”. Note that the projection onto the plane appears “square” as shown in Fig. 4(a); however, the structure is not a true square as depicted in Fig.4(b). In Figs 4(c) and 4(d) the interactions between chain and square are projected onto the and planes, respectively. In addition to the eight interactions shown here, an easy-plane type single-ion anisotropy was incorporated into the Hamiltonian.
LSWT analysis was performed on the dispersion relation measured along and by the CTAX spectrometer and that along and by the HRC spectrometer. The weighted sum of squared residuals is calculated as follows:
(2) |
Here, is the observed energy of the -th peak in a constant -scan/slice at = , is the calculated excitation energy that is the closest to , represents the FWHM of the Gaussian fit to the excitation at the -th data point, and is the number of data points. Given that multiple modes are in close proximity to one another, the FWHM was used as the weighting factor in the fitting process for parameter optimization in the analysis.
-
(meV) (meV) (meV) (meV) (meV) (meV) (meV) (meV) (meV)
The values were computed for all nine parameters varied within a wide range, and the parameters with the lowest values are summarized in Table 1. The values in parentheses represent the step size of each parameter used during the fitting. The dispersion curves calculated using the values from Table 1 are shown in Figs. 2(a)-2(d), and the INS spectra are shown in Figs. 2(e)-2(h). The calculations reasonably reproduce the experimental spectra, and we identified the spin Hamiltonian.
5 Discussion
In Fig. 1, , indicated by a yellow bond, is antiferromagnetic, and , indicated by a light blue bond, is ferromagnetic. The neighboring spins connected by are antiparallel to each other, which demonstrates the presence of frustration between and . Our calculations suggest that a noncollinear magnetic structure emerges when exceeds 30% of . In the identified spin Hamiltonian, however, is 16% of and it is not strong enough to induce the noncollinear structure.
The curvatures of INS spectra near 4 meV around (3, 1, 0) in Figs. 2(a), 2(b), 2(e), and 2(f) are negative. We confirmed that the curvature is, in fact, influenced by the presence or absence of frustration. In cases of positive –indicative of no frustration in the magnetic structure depicted in Fig. 1– the calculated curvature was positive, contrasting with the experimental observations. The presence of frustration, thus, qualitatively affects the spectra.
The magnitude of the ME tensor, which determines the size of the linear magnetoelectric effect, is related to the magnitude of the magnetic susceptibility[28]. It is posited that the ME tensor increases as the susceptibility increases. Magnetic frustration, which suppresses magnetic transitions, leads to an enhancement of the susceptibility at low temperatures. Therefore, the magnetic frustration identified in Mn3Ta2O8 in the present study is key to understand the linear magnetoelectric effect in the system.
6 Summary
We conducted inelastic neutron scattering experiments on linear ME material Mn3Ta2O8, which exhibits a collinear antiferromagnetic structure. The spin Hamiltonian was successfully identified by LSWT calculation. The presence of strong frustration was demonstrated.
acknowledgements
We are grateful to D. Kawana, T. Asami, and R. Sugiura for supporting us in the neutron scattering experiment at HRC and HER. The neutron experiment using HRC spectrometer at the Materials and Life Science Experimental Facility of the J-PARC was performed under a user program (Proposal No. 2021S01). A portion of this research used resources at the High Flux Isotope Reactor, a DOE Office of Science User Facility operated by the Oak Ridge National Laboratory (ORNL). Travel expenses for the neutron scattering experiments performed using CTAX at ORNL, USA was supported by the U.S.-Japan Cooperative Research Program on Neutron Scattering (proposal no. 2021-4). The neutron experiment using HER at JRR-3 was carried out by the joint research in the Institute for Solid State Physics, the University of Tokyo (Proposal No. 21403). H. Kikuchi was supported by Support for Pioneering Research Initiated by Next Generation (SPRING) of Japan Science and Technology Agency (JST). S. Hasegawa was supported by the Japan Society for the Promotion of Science through the Leading Graduate Schools (MERIT) This project was supported by JSPS KAKENHI Grant Numbers 19KK0069 and 21H04441.
References
References
- [1] ASTROV D 1960 JETP 11 708–709
- [2] Rado G T and Folen V J 1961 Phys. Rev. Lett. 7 310
- [3] Fiebig M 2005 J. Phys. D: Appl. Phys. 38 R123
- [4] Kornev I, Bichurin M, Rivera J, Gentil S, Schmid H, Jansen A and Wyder P 2000 Phys. Rev. B 62 12247
- [5] Mufti N, Blake G R, Mostovoy M, Riyadi S, Nugroho A A and Palstra T T M 2011 Phys. Rev. B 83 104416
- [6] Jiang S, Shan J and Mak K F 2018 Nat. Mater. 17 406
- [7] Kimura T, Goto T, Shintani H, Ishizaka K, Arima T and Tokura Y 2003 Nature 426 55–58
- [8] Taniguchi K, Abe N, Takenobu T, Iwasa Y and Arima T 2006 Phys. Rev. Lett. 97 097203
- [9] Sergienko I A and Dagotto E 2006 Phys. Rev. B 73 094434
- [10] Kenzelmann M, Harris A B, Jonas S, Broholm C, Schefer J, Kim S B, Zhang C L, Cheong S W, Vajk O P and Lynn J W 2005 Phys. Rev. Lett. 95 087206
- [11] Arkenbout A H, Palstra T T M, Siegrist T and Kimura T 2006 Phys. Rev. B 74 184431
- [12] Katsura H, Nagaosa N and Balatsky A V 2005 Phys. Rev. Lett. 95 057205
- [13] Arima T, Tokunaga A, Goto T, Kimura H, Noda Y and Tokura Y 2006 Phys. Rev. Lett. 96 097202
- [14] Okuyama D, Ishiwata S, Takahashi Y, Yamauchi K, Picozzi S, Sugimoto K, Sakai H, Takata M, Shimano R, Taguchi Y, Arima T and Tokura Y 2011 Phys. Rev. B 84 054440
- [15] Arima T 2007 J. Phys. Soc. Jpn. 76 073702
- [16] Murakawa H, Onose Y, Miyahara S, Furukawa N and Tokura Y 2010 Phys. Rev. Lett. 105 137202
- [17] Soda M, Matsumoto M, Månsson M, Ohira-Kawamura S, Nakajima K, Shiina R and Masuda T 2014 Phys. Rev. Lett. 112 127205
- [18] Kajimoto R, Mochizuki H, Yoshizawa H, Shintani H, Kimura T and Tokura Y 2005 J. Phys. Soc. Jpn. 74 2430–2433
- [19] Oh J, Le M D, Nahm H H, Sim H, Jeong J, Perring T G, Woo H, Nakajima K, Ohira-Kawamura S, Yamani Z, Yoshida Y, Eisaki H, Cheong S W, Chernyshev A L and Park J G 2016 Nat. Commun. 7 13146
- [20] Holbein S, Steffens P, Biesenkamp S, Ollivier J, Komarek A C, Baum M and Braden M 2023 Phys. Rev. B 108 104404
- [21] Esmaeilzadeh S, Grins J and Fitch A 1998 J. Mater. Chem. 8 2493–2497
- [22] Grins J, Esmaeilzadeh S, Berastegui P and Eriksson S 1999 J. Mater. Chem. 9 1575
- [23] Rickert K, Pozzi E A, Khanal R, Onoue M, Trimarchi G, Medvedeva J E, Hersam M C, Van Duyne R P and Poeppelmeier K R 2015 Inorg. Chem. 54 6513
- [24] Kimura K, Yagi N, Hasegawa S, Hagihala M, Miyake A, Tokunaga M, Cao H, Masuda T and Kimura T 2021 Inorg. Chem. 60 15078
- [25] Itoh S, Yokoo T, Satoh S, Yano S, Kawana D, Suzuki J and Sato T J 2011 Nucl. Instrum. Meth. A 631 90
- [26] Kawana D, Soda M, Yoshida M, Ikeda Y, Asami T, Sugiura R, Yoshizawa H, Masuda T, Hawai T, Ibuka S, Yokoo T and Itoh S 2018 J. Phys. Conf. Ser. 1021 012014
- [27] Toth S and Lake B 2015 J. Phys.: Condens. Matter 27 166002
- [28] Brown W F, Hornreich R M and Shtrikman S 1968 Phys. Rev. 168 574