Crystalline and magnetic structures, magnetization, heat capacity, and anisotropic magnetostriction effect in a yttrium-chromium oxide
Yinghao Zhu,1, ∗ Ying Fu,1, ∗ Bao Tu,1, 2, ∗ Tao Li,3 Jun Miao,4 Qian Zhao,1 Si Wu,1 Junchao Xia,1 Pengfei Zhou,1 Ashfia Huq,5 Wolfgang Schmidt,6 Defang Ouyang,7 Defang Ouyang,7 Zikang Tang,1, † Zhubing He,2, ‡ and Hai-Feng Li1, § 1Joint Key Laboratory of the Ministry of Education, Institute of Applied Physics and Materials Engineering, University of Macau, Avenida da Universidade, Taipa, Macao SAR 999078, China 2Department of Materials Science and Engineering, Shenzhen Key Laboratory of Full Spectral Solar Electricity Generation (FSSEG), Southern University of Science and Technology, No. 1088, Xueyuan Rd., Shenzhen 518055, Guangdong, China 3Neutron Scattering Technical Engineering Research Center, School of Mechanical Engineering, Dongguan University of Technology, Dongguan 523808, China 4School of Materials Science and Engineering, University of Science and Technology Beijing, Beijing 100083, China 5Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 6Forschungszentrum Ju¨lich GmbH, Ju¨lich Centre for Neutron Science at ILL, 71 avenue des Martyrs, Grenoble 38042, France 7State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences (ICMS), University of Macau, Avenida da Universidade, Taipa, Macao SAR 999078, China (Dated: September 22, 2020)
We have studied a nearly stoichiometric insulating Y0.97(2)Cr0.98(2)O3.00(2) single crystal by per- forming measurements of magnetization, heat capacity, and neutron diffraction. Albeit that the YCrO3 compound behaves like a soft ferromagnet with a coersive force of ∼0.05 T, there exist strong antiferromagnetic (AFM) interactions between Cr3+ spins due to a strongly negative paramagnetic Curie-Weiss temperature, i.e., –433.2(6) K. The coexistence of ferromagnetism and antiferromag- netism may indicate a canted AFM structure. The AFM phase transition occurs at TN = 141.5(1) K, which increases to TN(5T) = 144.5(1) K at 5 T. Within the accuracy of the present neutron- diffraction studies, we determined a G-type AFM structure with a propagation vector k = (1 1 0) and Cr3+ spin directions along the crystallographic c axis of the orthorhombic structure with space group Pnma below TN. At 12 K, the refined moment size is 2.45(6) µB, ∼82% of the theoretical 3+ saturation value 3 µB. The Cr spin interactions are probably two-dimensional Ising like within the reciprocal (1 1 0) scattering plane. Below TN, the lattice configuration (a, b, c, and V ) deviates largely downward from the Gr¨uneisenlaw, displaying an anisotropic magnetostriction effect and a magnetoelastic effect. Especially, the sample contraction upon cooling is enhanced below the AFM transition temperature. There is evidence to suggest that the actual crystalline symmetry of YCrO3 compound is probably lower than the currently assumed one. Additionally, we compared the t2g YCrO3 and the eg La7/8Sr1/8MnO3 single crystals for a further understanding of the reason for the possible symmetry lowering.
I. INTRODUCTION coupled families [5, 6]. Therefore, to build the actual structural and magnetic models may shed light on the nature of the ferroelectric phase transition. Magnetic materials with ferroelectricity are very inter- esting because a spontaneous electric polarization exists Most ferroelectric materials are perovskite-based ox- within them to give the substance extraordinary physical ides. In 1954, Looby and Katz replaced lanthanum in and electronic properties and a wide variety of applica- LaCrO3 compound with yttrium during searching for tions such as data storage, catalysts, fuel cells, and sen- new perovskite-type families and synthesized the YCrO3 arXiv:2007.08985v2 [cond-mat.str-el] 19 Sep 2020 sors [1–3]. The magnetic and ferroelectric properties can compound with an impurity of ∼2.5% Cr2O3 using NaCl be coupled with each other, therefore, an applied mag- as the flux under a hydrogen atmosphere [7]. Based on netic/electric field is able to switch the electric polariza- the observation of a very weak extra Bragg peak, they in- tion/magnetization. Above the ferroelectric phase tran- dexed the x-ray powder-diffraction pattern with a mon- sition, the ferroelectric materials are usually centrosym- oclinic cell (a = c = 7.61 A,˚ and b = 7.54 A)˚ by dou- metric structurally and behave as an ordinary dielec- bling the fundamental perovskite unit cell [7]. One year tric; below the phase transition, an electrically-polarized later, the crystal structure (Fig. 1) was determined to phase forms spontaneously with a noncentrosymmetric be orthorhombic (with space group P bnm) with unit- structure in the weakly-coupled systems [4] or can be cell constants a = 5.238, b = 5.518, and c = 7.54 A˚ [8]. induced by a magnetic phase transition in the strongly- This structure becomes more distorted as the pressure in- 2 creases [9]. The lanthanide orthochromites of general for- Y mula RECrO3 (RE = rare earth and Y) can be prepared O1 by four different self-propagating high-temperature syn- theses [10], i.e., the amorphous citrate precursor method [11], the conventional solid-state reaction method [12– Cr 14], the hydrothermal synthesis [15], and the microwave- O2 assisted technique [16]. Chemical substitution effect in the Y1−xMxCrO3 (M = Mg, Ca, Sr, Ba) compound [17] and defect chemistry of the Ca-doped YCrO3 compound [18] were investigated. The YCrO3 compound in forms of bulk and thin film was suggested to be a candidate ma- terial for high-temperature thermistors [19–21]. Among the catalysts of ABO3-type perovskite oxides (A = La, Y, Nd, Gd; B = Fe, Mn, Cr, Co) for the oxidation of 1,2-dichlorobenzene, the YCrO3 compound was found to be the most active catalyst and was the only one that b displayed no loss of its initial activity after several hours on stream [11]. The studies of nanocrystalline (un)doped a YCrO3 materials were reported [22–24]. The Nd-doped c YCrO3 nanoparticles display a semiconducting feature and an enhanced dc conductivity as the Nd content in- FIG. 1. Orthorhombic crystal structure (with space group creases, following the Dyre’s free energy barrier model P nma) with one unit cell (solid lines) and the AFM structure [24]. The magnetic configuration of the YCrO3 com- in one AFM unit cell with the propagation vector at k = (1 pound below TN = 140 K was proposed to be a canted 1 0) below TN = 141.5(1) K of the YCrO3 single crystal. The antiferromagnetic (AFM) structure with antisymmetric arrows on the Cr ions represent the spins of chromium. Both spin superexchanges [25–27]. It was reported that there the unit cells of orthorhombic and AFM structures are (a b existed a spin reorientation of the Cr3+ moments in the c). YCrO3 compound at ∼60 K, corresponding to a rotation of the AFM easy axis [28]. Ferrimagnetism was found KF flux in a platinum crucible [34–36]. in the half-doped YMn0.5Cr0.5O3 compound [29]. The In this paper, we have synthesized a centimeter-sized magnetic interactions in bulk YCrO3 compound were classified as classical three-dimensional isotropic Heisen- YCrO3 single crystal with a laser diode floating-zone berg universality according to the ab initio calculations (FZ) furnace [38] and performed measurements of the and Monte Carlo simulations based on a cubic structure chemical compositions, resistivity, and magnetization as with space group P m3m and lattice constant a = 3.76 A˚ functions of temperature and applied-magnetic field, heat [30]. To understand the ferroelectric anomaly occurring capacity, time-of-flight neutron-powder diffraction based at ∼473 K, the first-principles density functional theory on a spallation neutron source, and single-crystal neu- calculations found that the noncentrosymmetric mono- tron diffraction based on a reactor. The chemical com- positions of the grown YCrO single crystal are nearly clinic structure (with space group P 21) was the stablest 3 stoichiometric, and the YCrO compound is a robust in- one in view of its lowest energy [31]. The YCrO3 com- 3 pound was reported to be a relaxor ferroelectric material sulator. at about 450 K because of the local noncentrosymmetric structure [28]. The high-temperature magnetism (300– 980 K) and crystallographic information (321–1200 K) II. EXPERIMENTAL were studied by a time-of-flight neutron powder diffrac- tion [14], and it is of great interest that the structural The preparation of polycrystalline samples of YCrO3 information such as lattice constants, space group, bond compound has been described previously in detail [14]. angles, bond lengths, and the local distortion parame- The single crystals of YCrO3 compound were grown by ter have no response to the dielectric anomaly observed the FZ method [13] with a laser diode FZ furnace (Model: around 473 K [31]. There has been a long-standing de- LD-FZ-5-200W-VPO-PC-UM) [38, 39]. We employed bate about the decision as to which structural symmetry inductively coupled plasma with optical emission spec- is correct [7–38], which necessitates a growth of the high- troscopy (ICP-OES) analysis to quantitatively determine quality YCrO3 single crystal [34–38] and a careful study chemical compositions of the investigated single crystals. of its crystalline and magnetic properties. Previously, We measured resistivity, magnetization, and heat ca- small yttrium chromite single crystals with millimeter in pacity of the YCrO3 samples with a Quantum Design size were grown from the PbF2-B2O3 or the PbF2-B2O3- physical property measurement system (PPMS Dyna- 3 ] 2 2 ZFC, 0.01 T ZFC, 0.01 T 10 FC, 0.01 T 1 4 FC, 0.01 T -
) CW-fitting, ZFC, 0.01 T 2 -
T N /Cr) 3 B μ r 10 T /C 1 N [T ( [T B 2 μ ( M / M H 1 0 μ
(a) = (b) 1 -
0 0 0 50 100 150 200 250 300 -400 -300 -200 -100 0 100 200 300 T (K) T (K)
FIG. 2. (a) ZFC and FC magnetization (M ) of chromium ions in the single-crystal YCrO3 compound as a function of −1 temperature measured at µ0H = 0.01 T. (b) Corresponding ZFC and FC inverse magnetic susceptibility χ (circles) of chromium ions in the single-crystal YCrO3 compound versus temperature. The dash-dotted line indicates a CW behavior of the ZFC data at elevated temperatures between 200 and 300 K, which was extrapolated to χ−1 = 0 to show the PM Curie temperature θCW. The fit results were listed in Table I. In (a) and (b), TN = 141.5(1) K labels the AFM transition temperature at µ0H = 0.01 T, and the solid lines are guides to the eye.
10 4 (a) TABLE I. Theoretical quantum numbers of YCrO3 com- pound: spin S, orbital L, total angular momentum J, as well 2S+1 2 as the ground-state term LJ . Due to a quenching by the 5 hosted crystal field, the actual orbital angular momentum L ) 2 - 0 = 0 for the 3d ions in most cases, leading to the Land´efactor gJ = 2. We also summarized the theoretical (theo.) and mea- -2 sured (meas.) [Fig. 2(b)] values of effective (eff) chromium
/Cr /Cr 10 0 moment µ , PM Curie temperature θ , theoretical satura- B -4 eff CW μ
( tion (sat) chromium moment µsat theo., and AFM transition N´eeltemperatures (TN) at µ0H = 0.01 and 5 T. The refined
M -0.1 0.0 0.1 -5 chromium moment size (µmeas.) at 12 K with the AFM model as shown in Fig. 1 from our POWGEN study was listed. The Hysteresis loop @ 2 K numbers in parenthesis are the estimated standard deviations -10 of the last significant digit.
-7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 A YCrO3 single crystal 3+ μ0H (T) 3d ion Cr 3dn 3 FIG. 3. ZFC magnetic hysteresis loop of the single-crystal S 3/2 YCrO3 compound measured at 2 K. Inset (a) is the enlarged L 3 image of the narrow loop. J = L − S (Hund’s rule for free Cr3+) 3/2 2S+1 4 LJ F 3 2 Cool instrument). Resistivity measurements were carried gJ (quenched L = 0, J = S) 2 1 out with the standard four-probe method at zero field µeff theo. = gJ [S(S + 1)] 2 (µB) 3.873 from 2 to 300 K. The dc magnetization measurements µsat theo. = gJ S (µB) 3 at an applied-magnetic field of 0.01 T were carried out µeff meas. (µB) 3.95(2) with two modes from 5 to 295 K: One was after cooling θCW (K) –433.2(6) with zero magnetic field (ZFC), and the other was un- TN (at 0.01 T) 141.5(1) der the applied-magnetic field (FC). Magnetic hysteresis TN (at 5 T) 144.5(1) loop from 7 to –7 T and then back to 7 T was measured µ (12 K, POWGEN) (µ ) 2.45(6) at 2 K. Heat capacities were measured at 0 (2–273 K) meas. B and 5 T (2–205 K). We pulverized one grown YCrO3 single crystal (∼4 g) with a Vibratory Micro Mill (FRITSCH PUL- powder diffraction study on the POWGEN diffractome- VERISETTE 0) and performed a time-of-flight neutron- ter (SNS, USA) from 12 to 302 K at 0 T. The d band cov- 4
88 TABLE II. Refined structural parameters of the pulverized 0 T YCrO3 single crystal at 12, 145, and 300 K, including lat- tice constants, unit-cell volume, atomic positions, isotropic 66 5 T thermal parameters (B), bond lengths, bond angles, and the ) distortion parameter ∆ [14]. We listed the Wyckoff site of
.K each ion and the goodness of fit. The numbers in parenthesis are the estimated standard deviations of the last significant T mol 70 (a) N 44 digit. (J/
p 0 T T N(5T)
C 65 A pulverized YCrO3 single crystal (Orthorhombic, space group P nma (No. 62), Z = 4) 22 5 T 60 T (K) 12 145 300 ˚ 130 140 150 a (A) 5.5189(1) 5.5181(1) 5.5198(1) 0 b (A)˚ 7.5205(1) 7.5213(1) 7.5297(1) 0 40 80 120 160 200 240 280 c (A)˚ 5.2323(1) 5.2328(1) 5.2392(1) T (K) α(β, γ)(◦) 90 90 90 V (A˚3) 217.17(1) 217.18(1) 217.75(1) FIG. 4. Heat capacities of the single-crystal YCrO com- 3 Y(4c) x 0.0682(4) 0.0672(1) 0.0665(1) pound measured at 0 T (solid circles) and 5 T (void circles). The solid lines are guides to the eye. Inset (a) is the enlarged Y(4c) y 0.25 0.25 0.25 image around the AFM transition temperatures. The vertical Y(4c) z –0.0172(4) –0.0177(1) –0.0174(2) dashed lines show the detailed transition temperatures at the Y(4c) B (A˚2) 0.2 0.28(2) 0.47(2) fields of 0 and 5 T. Here, TN(0T) = 141.5(1) K at 0 T; by Cr(4b)(x, y, z) (0, 0, 0.5) (0, 0, 0.5) (0, 0, 0.5) comparison, at 5 T, T (5T) = 144.5(1) K. The solid lines are N ˚2 guides to the eye. Cr(4b) B (A ) 0.2 0.29(3) 0.39(3) O1(4c) x 0.4643(5) 0.4646(2) 0.4647(2) O1(4c) y 0.25 0.25 0.25 ers a range of 0.78–7.77 A.˚ Neutron holds a magnetic mo- O1(4c) z 0.1039(5) 0.1052(2) 0.1050(2) ment, thus, neutron scattering is a powerful technique for O1(4c) B (A˚2) 0.2 0.34(2) 0.51(2) solving magnetic structures [40–43]. The higher d band O2(8d) x 0.3020(4) 0.3020(1) 0.3021(1) (1.7–7.77 A)˚ is able to monitor all magnetic Bragg reflec- O2(8d) y 0.0539(2) 0.0538(1) 0.0536(1) tions that are used to refine the low-temperature mag- O2(8d) z –0.3065(4) –0.3067(1) –0.3066(1) netic structure. Data from the lower d band (0.78–3.00 O2(8d) B (A˚2) 0.2 0.34(2) 0.51(2) A)˚ give indication of all possible structural phase tran- Y-O11 (A)˚ 2.237(3) 2.231(1) 2.233(1) sitions. We analyzed all collected time-of-flight neutron- Y-O12 (A)˚ 2.276(4) 2.286(1) 2.290(1) powder diffraction data with the software of FULLPROF ˚ SUITE [44]. We refined scale factor, lattice constants, Y-O21 (A) (×2) 2.277(3) 2.272(1) 2.273(1) ˚ zero epithermal shift, background, peak profile shape, Y-O22 (A) (×2) 2.476(3) 2.479(1) 2.485(1) atomic positions, isotropic thermal parameters, and pre-
By quantitative ICP-OES measurements, we deter- Rp 6.40 5.46 6.03 mined the chemical compositions of the studied single Rwp 9.33 3.67 3.74 crystal as Y0.97(2)Cr0.98(2)O3.00(2). This implies that the Rexp 7.33 2.69 2.83 resultant single crystals of YCrO3 compound by our FZ χ2 1.62 1.86 1.75 method are nearly stoichiometric within the experimental accuracy. Therefore, during analyzing magnetization and time-of-flight neutron-powder diffraction data, we kept the stoichiometry of the synthesized YCrO3 samples be- ing the ideal one (i.e., 1:1:3). 5
300 K 12 8 POWGEN (SNS, USA) 8 4 4 0 0 ) s 145 K 12
count 8 Time-of-flight NPD 4 8 4 4 0 0 tensity (10 tensity
In Space group: Pnma 12 K 16 8 12
(1 1 0) (1 8
4 1 1) (0 4 0 0
1 2 3 2 4 6 8 d (Å)
FIG. 5. Observed (circles) and calculated (solid lines) time-of-flight neutron-powder diffraction (NPD) patterns of a pulverized YCrO3 single crystal, collected from the POWGEN diffractometer (SNS, USA) at 12, 145, and 300 K. The vertical bars mark the positions of nuclear (up, space group P nma) and magnetic (down, space group P -1) Bragg reflections, and the lower curves represent the difference between observed and calculated patterns.
B. Resistivity measurements tion curves measured at 0.01 T show very small values down to temperature ∼141.5 K, e.g., ZFC magnetization −7 −7 3+ We tried to measure possible resistivity in the YCrO3 = 4.860(4) × 10 and 7.310(6) × 10 µB/Cr at 295 single crystals with a multimeter at room temperature. and 142.3 K, respectively. Around 141.5 K, they increase Unfortunately, it was beyond the maximum range (106 sharply by ∼38% in a small thermal range of ∼5 K, fol- ohm) of the ohmmeter. In addition, our attempt to lowed by a smooth increase down to 5 K. This resembles measure it by the standard four-probe method with our the characteristic feature of a reasonable canted antifer- PPMS DynaCool system from 2 to 300 K was fruitless. romagnet and rules out the possibility for a ferrimagnet. −2 3+ Therefore, we conclude that the YCrO3 compound is a At 5 K, ZFC magnetization = 4.530(4) × 10 µB/Cr . robust insulator in our studied temperature range. A We calculated the inverse magnetic susceptibility −1 deeper understanding of the electronic states of conduct- χ = µ0H/M as shown in Fig. 2(b), where the nearly −1 ing VO [45] and insulating YCrO3 compounds necessi- linear increase of χ in the PM state at high tempera- tates more experimental work and theoretical band struc- tures obeys well the molar susceptibility via a CW law ture calculations. Perhaps both samples are the only two 2 C NAµ pure 3d3 compounds. χ(T ) = = eff , (1) T − θCW 3kB(T − θCW)
where C is the Curie constant, θCW is the PM Curie tem- 23 −1 C. Magnetization versus temperature perature, NA = 6.022 × 10 mol is the Avogadro’s p number, µeff = gµB J(J + 1) is the effective PM mo- −23 Fig. 2(a) shows magnetization measurements of a small ment, and kB = 1.38062 × 10 J/K is the Boltzmann piece of randomly-orientated YCrO3 single crystal. We constant. The fit by Eq. (1) was shown as the dash- 3+ transferred the unit of vertical axis into µB per Cr dotted line in Fig. 2(b). The fit parameters were listed ion. There exists no obvious difference between ZFC in Table I. and FC data. Upon cooling, ZFC and FC magnetiza- Figures 2(a) and 2(b) clearly indicate a sharp mag- 6
2.5 E. Heat capacity
Figure 4 shows the heat capacity measurements. At 0 2.0 T N T, with decreasing temperature, the measured heat ca- pacity decreases until TN = 141.5 K, followed by an ap-
/Cr) 1.5 pearance of a λ-shape peak. Below this, heat capacity B μ
( continues to decrease and gets flat below ∼24 K. The
ref 1.0 observation of the λ-shape peak indicates a phase tran- M sition. To reveal the nature of the phase transition, we measured heat capacity under an applied-magnetic field 0.5 of 5 T. As shown in Fig. 4(a), at 5 T, it was noted that the intensity of the λ-shape peak was reduced, accompanied 0.0 by a shift of the peak position from ∼138.8 K (0 T) to an 0 35 70 105 140 175 210 elevated temperature ∼139.9 K (5 T). This is the char- T (K) acteristic feature of a canted antiferromagnet. Thus, the phase transition is magnetic rather than structural. We FIG. 6. Refined chromium-moment size Mref (at 0 T) of determined TN (5 T) = 144.5(1) K, ∼3 K higher than the a pulverized YCrO3 single crystal versus temperature by the TN at 0 T. Quantitative analysis of the relationship be- software of FULLPROF SUITE [44]. The solid line is a guide tween values of T and applied-magnetic-field strengths to the eye. Error bars are standard deviations obtained from N necessitates more measurements. It is pointed out that our FULLPROF refinements in the Pnma symmetry. TN = 141.5(1) K labels the AFM transition temperature at zero with the measurement of magnetization versus temper- applied-magnetic field. ature, it is easy to determine the value of TN. From heat capacity measurements, the phase transition tem- perature is at the temperature point at which a kink netic phase transition. We determined the magnetic exists in the Cp − T curve as marked in Fig. 4(a). phase transition temperature as TN = 141.5(1) K. The Albeit that the magnetization and heat capacity mea- resultant PM CW temperature θCW = –433.2(6) K, in- surements show FM behaviors below TN, the net mag- dicating an existence of strong AFM correlations. We netic interaction strength inside YCrO3 compound is of calculated the frustrating parameter [46–51], i.e., f =| strongly AFM because θCW = –433.2(6) K, indicating a θCW | /TN = 3.061(5), which was consistent with our complex low-temperature magnetic structure. high-temperature magnetization study [14]. This value indicates that the low-temperature magnetic moments of 3+ Cr ions in YCrO3 compound are frustrated by compet- F. Time of flight neutron powder diffraction itive spin exchanges. As listed in Table I, the extracted effective PM moment µeff meas = 3.95(2) µB, a little bit To make the nature of the observed weak ferromag- larger than the calculated theoretical value µeff theo = netism clear and explore possible structural phase tran- 3.873 µB, which was acceptable within the present ex- sitions in the YCrO3 single crystal, we performed a time- perimental accuracy. of-flight neutron-powder diffraction study. The results were shown in Fig. 5. At the three temperatures as la- beled, i.e., below (12 K), around (145 K), and above (300 D. Magnetization versus applied magnetic field K) the magnetic transition temperature (∼141.5 K), all time-of-flight neutron-powder diffraction patterns were Figure 3 shows the measurement of magnetic hystere- well indexed by an orthorhombic structure with the space sis loop at 2 K. Figure 3(a) clearly exhibits the hysteresis group Pnma. There was no detectable peak splitting or loop whose shape is of a parallelogram. The loop locates an appearance of satellite reflections. This indicates that in a magnetic field range of ∼–0.1 to 0.1 T with a coer- no structural phase transition occurs in the YCrO3 single sive force of ∼0.05 T and a residual magnetism of ∼4.54 crystal as a function of temperature in the studied ther- −2 3+ ×10 µB/Cr . These small values demonstrate that mal regime. This is in agreement with our heat capacity the YCrO3 compound is a soft canted antiferromagnet measurements. Based on the observed magnetic Bragg at low temperatures. At 2 K and 7 T, the measured ZFC (0 1 1) and (1 1 0) peaks, as labeled in the right-bottom −2 3+ magnetization M = 9.603×10 µB/Cr . From 0.1 to 7 of Fig. 5, we established an AFM model with the propa- T, the measured magnetization almost increases linearly gation vector at k = (1 1 0) and the moment directions −3 −1 3+ with χ = M/µ0H = 7.24(1) ×10 µBT /Cr . We along the crystallographic c axis. The extracted mag- therefore estimated that reaching a complete magnetic netic structure was schematically drawn in Fig. 1. It is saturation state, an applied-magnetic field µ0H ≥ ∼408 pointed out that the magnetic (1 1 0) reflection is struc- T is required [52]. turally forbidden by the space group Pnma. We tried 7
300 (a) (b)
ts) 250 80 T N
coun 200 2 K
2 60 195 K 150 300 K (10 (1 1 0) 40 100 (1 1 0)
50 20 Intensity 0 0
0.92 0.96 1.00 1.04 1.08 units) (arb. Intensity Integrated 0 50 100 150 200 250 300 h in hk-scan at (1 1 0) (r.l.u.) T (K)
FIG. 7. (a) Representative longitudinal scans of the magnetic Bragg (1 1 0) reflection at three temperatures of 2, 195, and 300 K, from the D23 (ILL, France) study on a YCrO3 single crystal. The solid lines are guides to the eye. (b) Corresponding temperature-dependent integrated intensities of the magnetic Bragg (1 1 0) reflection. TN = 141.5(1) K labels the AFM transition temperature. The solid line was a fit to Eqs. (2) and (4) in the affiliated thermal regime. It was extrapolated to overall temperatures (dash-dotted line). The error bars in (a) and (b) are the standard deviations based on our measurements and fits. all possible canted AFM models, unfortunately, the cor- responding FULLPROF refinements were not successful. 80 nits) (1 1 0) The refined structural parameters were listed in Table II. It is pointed out that for the refinement of the data at 12 K, we constrained the isotropic thermal parameters (B) 60 40 of Y, Cr, O1, and O2 ions being the same as 0.2. Within the present experimental accuracy, we can only 40 tensity (arb. u (arb. tensity 20 determine a G-type AFM structure as shown in Fig. 1, where the nearest-neighbor Cr3+ spins are aligned an- (a) tiparallel. The directions of the AFM submoments are 0 20 132 136 140 144 along with the crystallographic c axis, i.e., the direction with the smallest lattice constant. It is hard to determine the possible canting angle. It was suggested that includ- In Integrated 0 ing a spin-orbital coupling on the quenched Cr3+ ground 0 70 140 210 280 350 state, an antisymmetric exchange interaction would cant T (K) the AFM moments along the crystallographic b axis, i.e., the direction with the largest lattice constant, according FIG. 8. Subtracted integrated intensity from the pure mag- to the Dzialoshinski-Moriya theory [53]. This leads to an netic contribution at the Bragg (1 1 0) peak position (void appearance of the weak ferromagnetism. pentagons), to see detailed analysis in the text. Inset (a) shows the enlarged image around TN from 132–145 K. The Figure 6 shows the refined moment size of Cr3+ ions solid line was a fit to the power-law Eq. (5) in the affiliated in the YCrO3 single crystal, extracted from our time- thermal regime. The error bars are the propagated standard of-flight neutron-powder diffraction study. As listed in deviations based on our calculations. Table I, the refined moment size at 12 K is 2.45(6) µB, ∼82% of the theoretical saturation moment (3 µB), in agreement with our conclusion that there exists a mag- netic frustration in YCrO3 compound and the studies with x-ray magnetic circular dichroism and absorption spectroscopies [30] where the computed values of spin and orbital moments are 2.38 µB and –0.094 µB, respectively, and that the total magnetic moment has little contri- ture upon warming until a radical disappearance around bution from the orbital component. As temperature in- 141.5 K, the onset temperature of the AFM transition creases, the refined moment size remains a plateau up to (Figs. 2 and 6). Above TN (0 T) = 141.5 K, intensities ∼50 K, followed by a gradual diminution with tempera- of the two magnetic Bragg reflections were undetectable. 8
G. Single crystal neutron diffraction 5.24 5.524 (a) c 7.533 An orientated YCrO3 single crystal with the (HK 0) reciprocal-lattice vector in the scattering plane was used (Å) nts 5.522 b to perform a single-crystal neutron diffraction study on 7.524 5.23 the D23 single-crystal diffractometer (ILL, France). Fig- a ure 7(a) shows some longitudinal scans of the magnetic 5.520 Bragg (1 1 0) reflection. At 2 K (below TN), we ob- 7.515 served a very strong peak, indicating a formation of the consta Lattice 5.518 5.22 AFM structure. At 195 and 300 K (above TN), the inten- sity of the magnetic (1 1 0) peak decreased sharply but ) (b) did not disappear. It is interesting that there still ex- 3 217.8 ists detectable intensity of the Bragg (1 1 0) peak above (Å V 217.6 TN. We ruled out the λ/2 contamination. As the fore- T N going remark, this reflection is forbidden by the space group Pnma. Therefore, the existence of the Bragg (1 217.4
1 0) reflection above TN indicates that the actual crys- volume ll
ce 217.2 talline structure of YCrO3 compound may be lower than - the orthorhombic structure with space group Pnma. Our POWGEN (SNS, USA) studies also demonstrate that the scattering ability of a Unit 217.0 single crystal is much higher than that of the correspond- 0 70 140 210 280 ing pulverized powder sample. T (K) At 195 and 300 K, the observed Bragg (1 1 0) re- flection forbidden structurally by the Pnma symmetry FIG. 9. (a) Temperature-dependent lattice constants a, b, was treated to be from a pure nuclear contribution. and c of a pulverized YCrO3 single crystal. (b) Correspond- The temperature variation of this contribution depends ing anomalous unit-cell volume V expansion with tempera- mainly on the thermal dynamic vibrations of related ture. The solid lines in (a) and (b) are theoretical estimates atoms, i.e., Debye-Waller (DW) factors. The falloff of of the variation of structural parameters using the Gr¨uneisen the temperature-weakened intensity at a certain scatter- model with Debye temperature of θD = 580 K that is the same value as the one reported previously [14]. TN = 141.5(1) K ing vector Q almost decays exponentially and can be labels the AFM transition temperature. The error bars in (a) estimated by and (b) are the standard deviations obtained from our FULL- −2W (Q,T ) PROF refinements with the Pnma structural symmetry. I = I0e , (2) where the exponential part is the DW factor, and contribution below TN as shown in Fig. 8. The resultant 2 2 Z ~ Q ~ω Z(ω) magnetic intensity above TN is approximately zero within 2W (Q, T ) = coth dω, (3) accuracy, which in turn supports the above subtraction. 2M 2kBT ω The extracted integrated intensity (I ) of Bragg (1 1 0) where ~ = 1.054589 × 10−34 J.s is the Planck constant reflection from the pure magnetic contribution can be fit divided by 2π, M is the atomic mass, and Z(ω) is the to a power-law equation [60, 61] phonon density of states [54]. At high temperatures, it 2β is given simply by |T − TN| I(T ) = I0 , (5) TN 3 2Q2 2W = ~ T, (4) Mk Θ2 where TN is the value of the N´eeltemperature, and β B W is the critical exponent. Our fit with Eq. (5) to the ex- where ΘW is the effective Debye temperature [55]. There- tracted data in a narrow thermal range from 125 to 140 fore, we treated practically the DW factor with a linear K, shown as the solid line in Fig. 8, produces a N´eel response to temperature, as verified in previous studies temperature TN1 = 140.0(1) K, and the critical exponent [56–59]. We fit the integrated intensities of the Bragg 2β = 0.215(6), indicating a second-order type phase tran- (1 1 0) reflection above TN to Eqs. (2) and (4), shown sition and probably two-dimensional Ising-like spin inter- as the solid line in Fig. 7(b), and extrapolated the fit to actions existing within the reciprocal (1 1 0) scattering the entire temperature range (shown as the dash-dotted plane [41, 60]. The fitting procedure was as follow: First, line). we kept I0 = 80 and TN = 141.5 K and allowed β to vary; Furthermore, we subtracted the corresponding nuclear finally, we fit all parameters together. For comparison, component from the total scattered intensity at the we further fit the data in three temperature ranges of Bragg (1 1 0) peak position to extract the pure magnetic 125–140 K, 130–140 K, and 135–140 K. No clear differ- 9
2.04 La7/8Sr1/8MnO3 YCrO3
) T N Å 2.01 Mn-O22 T JT Cr-O21 Cr-O22 1.98 Mn-O1 ond lengths ( lengths ond B 1.95 Mn-O21 Cr-O1 (a) (b) 1.92 0 50 100 150 200 250 300 0 50 100 150 200 250 300 T (K)
FIG. 10. (a) Three bond lengths of Mn-O1, Mn-O21, and Mn-O22 in the La 7 Sr 1 MnO3 single crystal [40] versus temperature. 8 8 TJT ≈ 180–270 K denotes the regime of the Jahn-Teller effect. (b) Corresponding bond lengths in the YCrO3 single crystal as a function of temperature from this study. TN = 141.5(1) K labels the AFM transition temperature. The error bars in (a) and (b) were from our FULLPROF refinements.
ences exist in the values of the refined TN and β, which O21 validates our choice of the temperature range of 125–140 K for the final fitting. O2 O22 Cr Cr As shown in Fig. 8(a), it is interesting to note that above TN1, there exists weaker critical scattering over a O22 range of temperature up to 145 K. O21
H. Anisotropic magnetostriction effect O1 O22 The refined lattice parameters a, b, and c, as well as the unit-cell volume V, from our time-of-flight neutron- O22 O21 Cr powder diffraction studies were shown in Fig. 9 (void O21 symbols). Upon cooling, the refined (Re) a, b, c, and V almost shrink linearly down to TN at which a cusp appears.
As in the foregoing discussions, the YCrO3 compound is an insulator. We thus neglected the electronic contri- bution to the thermal expansion of its lattice configura- O1 tion (ε). The temperature variation of the nonmagnetic contribution component is then mainly from phonons. FIG. 11. Schematic illustration of the three Cr-O bonds (Cr- This can approximately be estimated according to the O1, Cr-O21, and Cr-O22), as well as the two bond angles Cr- Gr¨uneisenrules at zero pressure with a second-order fash- O-Cr (Cr-O1-Cr and Cr-O2-Cr) in the orthorhombic structure ion [62–64] of a YCrO3 single crystal. In this structural symmetry (with space group Pnma), Cr ions in YCrO3 compound have the same Wyckoff site, 4b (0 0 0.5), as that of the Mn ions in U ε(T ) = ε0 + ε0 , (6) La 7 Sr 1 MnO3 compound [13, 40, 56, 57]. Q − BU 8 8
where ε0 is the lattice parameter at zero Kelvin, and the internal energy U can be calculated with the Debye 10
162 2.395 TJT (a) La7/8Sr1/8MnO3 ) ) La Sr MnO 7/8 1/8 3 Å Å ( ( 2.775 > 161 > O O
2.390 - - Y La Mn-O2-Mn TN 160 T 2.770 JT 2.385 ength ength
) 159 Mn-O1-Mn Bond l Bond ( Bond length length Bond 2.380 YCrO3 2.765 158 (b) YCrO3 0 50 100 150 200 250 300 145.4 T N T (K)
Bond angles Bond Cr-O2-Cr 145.2 FIG. 13. Comparison of the averaged bond lengths of Y-O in YCrO3 (left, from the present study) and La-O in La 7 Sr 1 MnO3 [40] (right) single crystals. TJT ≈ 180–270 K 145.0 8 8 denotes the regime of the Jahn-Teller effect of La 7 Sr 1 MnO3 8 8 144.8 compound. TN = 141.5(1) K labels the AFM transition tem- perature of YCrO3 compound. The error bars are the calcu- Cr-O1-Cr lated standard deviations. The solid lines are guides to the 144.6 eye. It is clear that the bond length of hY-Oi is shorter than that of the hLa-Oi bond beyond statistics. 144.4 0 50 100 150 200 250 300 T (K) pound. Below TN, the magnetically-driven additional de- creases of a, b, and c jointly result in an enhanced sam- 12K 12K FIG. 12. (a) Temperature-dependent bond angles of Mn-O1- VRe −VGr ple contraction upon cooling, e.g., 12K ≈ –7.43 × Mn and Mn-O2-Mn in the La 7 Sr 1 MnO3 single crystal [40]. VGr 8 8 −4 TJT ≈ 180–270 K denotes the regime of the Jahn-Teller effect. 10 , signifying a magnetoelastic effect [65] and a local- (b) Temperature-dependent bond angles of Cr-O1-Cr and Cr- ized nature of the t2g moments, opposite the case in the O2-Cr in the YCrO3 single crystal from the present study. GdSi metallic compound [62]. TN = 141.5(1) K labels the AFM transition temperature. The error bars in (a) and (b) are the standard deviations from refinements. The solid lines in (a) and (b) are guides to the I. Comparison between t2g YCrO3 and eg eye. La7/8Sr1/8MnO3 compounds approximations As shown in Fig. 10, it is of interest to compare the 3+ 4+ 1 3 0 3 bond lengths of the 3d Mn /Mn (egt2g/egt2g) ions 3 ΘD Z T 3 in La 7 Sr 1 MnO3 compound [40, 66] with those of the T x 8 8 U(T ) = 9Nk T dx, (7) 3+ 0 3 B x Cr (e t ) ions in YCrO3 compound. The structural ΘD 0 e − 1 g 2g parameters of La 7 Sr 1 MnO3 compound [40] were from 8 8 where N = 5 is the number of atoms per formula unit, a neutron-powder diffraction study on samples pulver- and ΘD is the Debye temperature. With Eqs. (6) and ized from a single crystal. This is thus comparable to (7), we fit the lattice parameters of YCrO3 compound the results from the present study. Both bond lengths of in the PM state (above TN ≈ 141.5 K) and extrapolated Mn-O21 and Mn-O22 [Fig. 10(a)] and the bond angle of the fits to overall temperatures as shown in Fig. 9 (solid Mn-O1-Mn [Fig. 12(a)] respond readily to the Jahn-Teller lines). For example, the fitting for the unit-cell volume effect that occurs in La 7 Sr 1 MnO3 compound within a 3 −18 8 8 V results in V0 ≈ 217.14 A˚ , Q ≈ 7.57 × 10 J, and temperature range of ∼180–270 K, whereas, those in B ≈ –40.73. The different variations in a, b, and c be- YCrO3 compound keep nearly constants, and both Cr- low TN in contrast to our theoretical estimates by the O bond lengths [Fig. 10(b)] and Cr-O-Cr bond angles 12K 12K aRe −aGr Gr¨uneisen(Gr) law [63, 64] (solid lines), e.g., 12K [Fig. 12(b)] exhibit no response to the AFM transition, aGr 3+ 12K 12K 12K 12K consistent with the fact that Cr ions don’t have the −4 bRe −bGr −4 cRe −cGr ≈ –2.73 × 10 , 12K ≈ –2.19 × 10 , and 12K bGr cGr Jahn-Teller effect. The values of Cr-O1 bond lengths are ≈ –2.67 × 10−4, indicate an anisotropic magnetostriction similar to those of the Mn-O21, and Cr-O21 and Cr-O22 effect and that magnetic anisotropy exists in YCrO3 com- to Mn-O1 (Fig. 10). 11
As shown in Figs. 11 and 12, in contrast to the bond parison, below TN, the lattice configuration (a, b, c, and ◦ angle of Mn-O-Mn, ∠Cr-O-Cr decreases hugely by ∼15 , V ) deviates largely downward from the Gr¨uneisenlaw, which in our opinion corresponds intimately to a possible displaying an anisotropic magnetostriction effect along lowering of the crystalline symmetry in the YCrO3 sin- the crystallographic a, b, and c axes and a magnetoelas- gle crystal. The relatively shorter bond length of hY-Oi tic effect with the unit-cell volume V. Especially, upon 3+ (Fig. 13) introduces an immense mismatch between Y cooling, the sample contraction is enhanced below TN. and Cr3+ sites, leading to a huge chemical pressure and In the whole studied temperature range of 12–302 K, driving the subsequent rotating and tilting of the CrO6 we did not find any crystalline structural phase transition octahedra. with the neutron-powder diffraction study, whereas by our single-crystal neutron diffraction study, we observed clearly the existence of the Bragg (1 1 0) peak above IV. CONCLUSIONS the magnetic phase transition temperature 141.5(1) K. This peak persists up to 300 K and is forbidden by In summary, we have grown a nearly stoichiometric the crystalline orthorhombic structure (with space group Y0.97(2)Cr0.98(2)O3.00(2) single crystal by a laser diode FZ P nma). This implies that the actual crystalline structure furnace. There are three electrons locating on the 3d t2g of YCrO3 compound is probably lower than the present 3+ orbitals of Cr ions, therefore, the YCrO3 compound is one. To figure out the reasons for a possible lowering of a robust insulator. Although the measurements of the the structural symmetry in the YCrO3 single crystal, we applied-magnetic-field dependent heat capacity as well compared the t2g YCrO3 and the eg La7/8Sr1/8MnO3 sin- as the magnetization versus temperature and applied- gle crystals. It is pointed out that with a limited number magnetic field show the character of a very soft ferro- of the observed magnetic Bragg peaks, it is hard to de- magnet with a coersive force of ∼0.05 T, the extracted termine the canting degree of the AFM Cr3+ spins. This PM CW temperature, θCW = –433.2(6) K, by the fit can be addressed by a further time-of-flight single-crystal with a CW law is strongly negative with the frustrating neutron-diffraction study. To determine the actual crys- parameter f =| θCW | /TN = 3.061(5), and the mea- talline structure is not an easy job, but it would shed sured magnetization at 2 K and 7 T is only ∼3.2% of light on the dielectric anomaly of YCrO3 compound. the theoretical saturation moment. These indicate that 3+ the spin moments of Cr ions in YCrO3 compound are magnetically frustrated. The consistency between ACKNOWLEDGEMENTS the effective PM moment, µeff meas = 3.95(2) µB, and the theoretically-calculated value, µ = 3.873 µ , eff theo B T.L. acknowledges the National Natural Science Foun- validates our results concluded in the framework of the dation of China (Grant No. 11604214), the Foundation of CW-law fitting. By magnetization measurements, we de- Department of Education of Guangdong Province (Grant termined the magnetic phase transition temperature as No. 2018KTSCX223), and the Foundation of Depart- T = 141.5(1) K at an applied-magnetic field of 0.01 T. N ment of Science and Technology of Guangdong Province This is in agreement with our neutron-powder and single- (Grant No. 2020A1515010814). J.M. acknowledges crystal diffraction studies. The magnetic transition tem- the National Science Foundation of China (Grant No. perature was pushed upward to T (5T) = 144.5(1) K at N 11974042). D.O. acknowledges financial support from 5 T, increased by ∼3 K. the Science and Technology Development Fund, Macao With our neutron-powder diffraction study, we have SAR (File No. 0029/2018/A1). H.-F.L. acknowledges established an AFM structure with the propagation vec- the University of Macau (Files No. SRG2016–00091–FST tor at k = (1 1 0) and the same unit cell as that of and No. MYRG2020–00278–IAPME), the Science and the crystalline structure (with space group P nma). The Technology Development Fund, Macao SAR (Files No. direction of the Cr3+ spin moments is along the crys- 063/2016/A2, No. 064/2016/A2, No. 028/2017/A1, and tallographic c axis. The refined moment size is 2.45(6) No. 0051/2019/AFJ), and the Guangdong–Hong Kong– µ at 12 K, ∼82% of the theoretical saturation value B Macao Joint Laboratory for Neutron Scattering Science 3 µ . This is consistent with the fact that a magnetic B and Technology. This research used resources at the Spal- frustration exists in YCrO compound. By fitting inte- 3 lation Neutron Source, a Department of Energy Office of grated intensities of the magnetic Bragg (1 1 0) reflec- Science User Facility operated by Oak Ridge National tion extracted from the pure magnetic contribution with Laboratory. a power law, we found that the Cr3+ spin interactions were probably two-dimensional Ising like within the re- ciprocal (1 1 0) scattering plane. Above TN = 141.5(1) K, the refined lattice constants a, b, and c, as well as the unit-cell volume V, agree well with the Gr¨uneisenrules ∗ These three authors contributed equally. at zero pressure with a second-order fashion. By com- † [email protected] 12
‡ [email protected] Chem. 183, 1863 (2010). § [email protected] [29] L. Hao, L. Yang, M.-H. Lee, T.-H. Lin, Z. Zhang, X. Xie, [1] O. Auciello, J. F. Scott, and R. Ramesh, Phys. Today and H. Zhu, J. Alloys Compd. 601, 14 (2014). 51, 22 (1998). [30] S. Sidi Ahmed, G. D. Ngantso, M. Boujnah, A. Beny- [2] B. B. Vanaken, T. T. M. Palstra, A. Filippetti, and N. oussef, and A. El Kenz, Phys. Lett. A 383, 121 (2019). A. Spaldin, Nat. Mater. 3, 164 (2004). [31] C. R. Serrao, A. K. Kundu, S. B. Krupanidhi, U. V. [3] H. J. Zhao, W. Ren, Y. Yang, J. ´I˜niguez,X. M. Chen, Waghmare, and C. N. R. Rao, Physical Review B 72, and L. Bellaiche, Nat. Commun. 5, 4021 (2014). 220101(R) (2005). [4] A. S. Gibbs, K. S. Knight, and P. Lightfoot, Phys. Rev. [32] L. Katz, Acta Cryst. 8, 121 (1955). B 83, 094111 (2011). [33] J. P. Van Der Ziel and L. G. Van Uitert, Phys. Rev. 179, [5] M. Kenzelmann, A. B. Harris, S. Jonas, C. Broholm, 343 (1969). J. Schefer, S. B. Kim, C. L. Zhang, S.-W. Cheong, O. [34] J. P. Remeika, J. Am. Chem. Soc. 78, 4259 (1956). P. Vajk, and J. W. Lynn, Phys. Rev. Lett. 95, 087206 [35] W. H. Grodkiewicz and D. J. Nitti, J. Am. Ceram. Soc. (2005). 49, 576 (1966). [6] M. V. Mostovoy, Phys. Rev. Lett. 96, 067601 (2006). [36] N. D. Todorov, M. V. Abrashev, V. G. Ivanov, G. G. [7] J. T. Looby and L. Katz, J. Am. Chem. Soc. 76, 6029 Tsutsumanova, V. Marinova, Y.-Q. Wang, and M. N. (1954). Iliev, Phys. Rev. B 83, 224303 (2011). [8] S. Geller and E. A. Wood, Acta Cryst. 9, 563 (1956). [37] S. Sugano, K. Aoyagi, and K. Tsushima, J. Phys. Soc. [9] M. Ardit, G. Cruciani, M. Dondi, M. Merlini, and P. Jpn. 31, 706 (1971). Bouvier, Phys. Rev. B 82, 064109 (2010). [38] Y. H. Zhu, S. Wu, Z. K. Tang, and H.-F. Li, A method of [10] M. V. Kuznetsov and I. P. Parkin, Polyhedron 17, 4443 centimeter-sized single crystal growth of chromate com- (1998). pounds and related storage device, China Patent No. [11] K. Poplawski, K. Schnitzlein, J. Lichtenberger, F. J. Keil, 201911281088.8, (2019). and M. D. Amiridis, Catal. Today 62, 329 (2000). [39] S. Wu, Y. H. Zhu, H. S. Gao, Y. Xiao, J. C. Xia, P. F. [12] D. Westphal, A. Laske, S. Jakobs, and U. Guth, Ionics Zhou, D. F. Ouyang, Z. Li, Z. Q. Chen, Z. K. Tang, and 6, 346 (2000). H.-F. Li, ACS Omega 5, 16584 (2020). [13] H. Li, Synthesis of CMR manganites and order- [40] H.-F. Li, Y. Su, Y. G. Xiao, J. Persson, P. Meuffels, and ing phenomena in complex transition metal oxides Th. Br¨uckel, Eur. Phys. J. B 67, 149 (2009). (Forschungszentrum J¨ulich GmbH Press, J¨ulich, 2008). [41] Y. Xiao, Y. Su, H.-F. Li, C. M. N. Kumar, R. Mittal, J. [14] Y. Zhu, S. Wu, B. Tu, S. Jin, A. Huq, J. Persson, H. S. Persson, A. Senyshyn, K. Gross, and Th. Brueckel, Phys. Gao, D. Ouyang, Z. B. He, D.-X. Yao, Z. K. Tang, and Rev. B 82, 094437 (2010). H.-F. Li, Phys. Rev. B 101, 014114 (2020). [42] Q. Zhang, W. Tian, H. Li, J.-W. Kim, J. Yan, R. W. Mc- [15] K. Sardar, M. R. Lees, R. J. Kashtiban, J. Sloan, and R. Callum, T. A. Lograsso, J. L. Zarestky, S. L. Bud0ko, R. I. Walton, Chem. Mater. 23, 48 (2011). J. McQueeney, and D. Vaknin, Phys. Rev. B 88, 174517 [16] J. Prado-Gonjal, R. Schmidt, J.-J. Romero, D. Avila,´ U. (2013). Amador, and E. Mor´an,Inorg. Chem. 52, 313 (2013). [43] W. H. Ji, L. Yin, W. M. Zhu, C. M. N. Kumar, C. Li, H.- [17] W. J. Weber, C. W. Griffin, and J. L. Bates, J. Mater. F. Li, W. T. Jin, S. Nandi, X. Sun, Y. Su, Th. Br¨uckel, Res. 1, 675 (1986). Y. Lee, B. N. Harmon, L. Ke, Z. W. Ouyang, and Y. [18] G. F. Carini II, H. U. Anderson, D. M. Sparlin, and M. Xiao, Phys. Rev. B 100, 134420 (2019). M. Nasrallah, Solid State lonics 49, 233 (1991). [44] J. Rodr´ıguez-Carvajal, Phys. B 192, 55 (1993). [19] M. Kagawa, Y. Kato, and Y. Syono, J. Aeroso Sci. 28, [45] N. N. Greenwood and A. Earnshaw, 475 (1997). ChemistryoftheElements (Elsevier Ltd, Oxford, [20] J. H. Kim, H.-S. Shin, S.-H. Kim, J.-H. Moon, and B.-T. 1997). Lee, Jpn. J. Appl. Phys. 42, 575 (2003). [46] A. P. Ramirez, HandbookofMagneticMaterials, Vol. [21] J. Seo, Y. Ahn, and J. Y. Son, Ceram. Int. 41, 12471 13, edited by K. Buschow (Elsevier, Amsterdam, 2001), (2015). p. 423. [22] V. Bedekar, R. Shukla, and A. K. Tyagi, Nanotechnology [47] A. P. Ramirez, Annu. Rev. Mater. Sci. 24, 453 (1994). 18, 155706 (2007). [48] H. T. Diep, F rustratedSpinSystems (World Scientific, [23] J. Bahadur, D. Sen, S. Mazumder, R. Shukla, and A. K. Singapore, 2004). Tyagi, J. Phys. Condens. Mat. 20, 345201 (2008). [49] M. P. Shores, E. A. Nytko, B. M. Bartlett, and D. G. [24] R. Sinha, S. Basu, and A. K. Meikap, Phys. E Low- Nocera, J. Am. Chem. Soc. 127, 13462 (2005). Dimens. Syst. Nanostruct. 113, 194 (2019). [50] H.-F. Li, C. Zhang, A. Senyshyn, A. Wildes, K. Schmalzl, [25] V. M. J¨udinand A. B. Sherman, Solid State Commun. W. Schmidt, M. Boehm, E. Ressouche, B. Hou, P. Meuf- 4, 661 (1966). fels, G. Roth, and Th. Br¨uckel, Front. Phys. 2, 42 (2014). [26] E. F. Bertaut, G. Bassi, G. Buisson, P. Burlet, J. Chap- [51] H.-F. Li, A. Senyshyn, O. Fabelo, J. Persson, B. Hou, pert, A. Delapalme, J. Mareschal, G. Roult, R. Aleonard, M. Boehm, K. Schmalzl, W. Schmidt, J.-P. Vassalli, P. R. Pauthenet, and J. P. Rebouillat, J. Appl. Phys. 37, Thakuria, X. Sun, L. Wang, G. Khazaradze, B. Schmitz, 1038 (1966). C. Zhang, G. Roth, J. G. Roca, and A. Wildes, J. Mater. [27] T. Morishita and K. Tsushima, Phys. Rev. B 24, 341 Chem. C 3, 7658 (2015). (1981). [52] H.-F. Li, npj Comput. Mater. 2, 16032 (2016). [28] A. Dur´an,A. M. Ar´evalo-L´opez, E. Castillo-Mart´ınez, [53] A. H. Cooke, D. M. Martin, and M. R. Wells, J. Phys. M. Garc´ıa-Guaderrama, E. Moran, M. P. Cruz, F. C: Solid State Phys. 7, 3133 (1974). Fern´andez, and M. A. Alario-Francob, J. Solid State [54] W. L. Stephen, T heoryofneurtronscatteringfromcondensedmatter (Clarendon Press, Oxford, 1984). 13
[55] C. Kittel, QuantumT heoryofSolids (Wiley, New York, [60] M. F. Collins, MagneticCriticalScattering (Oxford Uni- 1987). versity Press, London, 1989). [56] H. F. Li, Y. Su, J. Persson, P. Meuffels, J. M. Walter, R. [61] S. J. Payne, M. E. Hagen, and M. J. Harris, J. Phys.: Skowronek, and Th. Br¨uckel, J. Phys.: Condens. Matter Condens. Matter 8 91 (1996). 19, 016003 (2006). [62] H.-F. Li, Y. Xiao, B. Schmitz, J. Persson, W. Schmidt, [57] H. F. Li, Y. Su, J. Persson, P. Meuffels, J. M. Walter, R. P. Meuffels, G. Roth, and Th. Br¨uckel, Sci. Rep. 2, 750 Skowronek, and Th. Br¨uckel, J. Phys.: Condens. Matter (2012). 19, 176226 (2007). [63] D. C. Wallace, Phys. Rev. B 58, 15433 (1998). [58] A. Tari, T heSpecificHeatofMatteratLowT emperatures [64] L. Vo˘cadlo,K. S. Knight, G. D. Price, and I. G. Wood, (Imperial College Press, London, 2003). Phys. Chem. Miner. 29, 132 (2002). [59] H.-F. Li, A. Wildes, B. Y. Hou, C. Zhang, B. Schmitz, P. [65] H.-F. Li, J. Q. Yan, J. W. Kim, R. W. McCallum, T. A. Meuffels, G. Rothb, and Th. Br¨uckel, RSC Adv. 4, 53602 Lograsso, and D. Vaknin, Phys. Rev. B 84, 220501(R) (2014). (2011). [66] A. Kotani, H. Nakajima, K. Harada, Y. Ishii, and S. Mori, Phys. Rev. B 95, 144403 (2017).