Coupling between ferroelasticity and magnetization in two-dimensional organic- inorganic perovskites (C6H5C2H4NH3)2MCl4(M = Mn, Cu, Fe)

ISSN 2050-7526 Materials for optical, magnetic and electronic devices Journal of Materials Chemistry C rsc.li/materials-c PAPER Katsuya Inoue et al . Coupling between ferroelasticity and magnetization in two-dimensional organic–inorganic perovskites (C 6 H 5 C 2 H 4 NH 3 ) 2 M Cl 4 ( M = Mn, Cu, Fe) Volume 13 Number 6 14 February 2025 Pages 2539–3074 This journal is © The Royal Society of Chemistry 2025 J. Mater. Chem. C, 2025, 13, 2661–2672 | 2661 Cite this: J. Mater. Chem. C, 2025, 13, 2661 Coupling between ferroelasticity and magnetization in two-dimensional organic– inorganic perovskites (C 6 H 5 C 2 H 4 NH 3 ) 2 MCl 4 (M= Mn, Cu, Fe)† Naoto Tsuchiya, a Saya Aoki, a Yuki Nakayama, a Goulven Cosquer, bc Sadafumi Nishihara, abd Miguel Pardo-Sainz, ef Jose ´Alberto Rodrı ´guez-Velamaza ´n, g Javier Campo ce and Katsuya Inoue * abc Materials with coexistence of two or more ferroic orders are known as multiferroics. Magneto-elastic multiferroics, where ferromagnetism and ferroelasticity coexist, have been rarely reported previously. We studied the magneto-elastic multiferroic properties of two-dimensional organic–inorganic perovskites having the formulas (PEA) 2 MnCl 4 , (PEA) 2 CuCl 4 and (PEA) 2 FeCl 4 (PEA = C 6 H 5 C 2 H 4 NH 3 ). All three exhibited ferroelasticity but the manganese and iron compounds showed canted antiferromagnetism and the copper one showed ferromagnetism. Also, only (PEA) 2 FeCl 4 displayed a shift of magnetization when the sample was cooled in a magnetic field from above the magnetic ordering temperature. We propose that the magnetization shift originates from the coupling between ferroelasticity and magnetization via spin– orbit coupling (SOC). This work would shed light on understanding the coupling mechanism between ferroelasticity and magnetization towards the interesting role of SOC in ferroelastic materials. Introduction Organic–inorganic perovskites have gained worldwide attention for their multiferroic, photovoltaic and semiconducting properties. 1–4 They are composed of organic ammonium cations and metal halide octahedra. While the organic group offers various properties such as elasticity and efficient luminescence, the inorganic parts provide important properties such as thermal stability and magnetic and dielectric properties. 5–7 Recent studies in organic–inorganic perovskites have focused not only on their respective characteristics but also on the interactions between them. Among the organic–inorganic perovskites, two-dimensional organic–inorganic perovskites can be a promising platform for functional material design due to their ability to introduce various sizes of organic cations. 8,9 Recently, these perovskites have been recognised as candidates for multiferroic materials. 10–15 Multiferroic materials exhibit a direct correlation between two or more ferroic orders such as ferromagnetism, ferroelectricity, ferroelasticity and ferrotoroidicity. 1,2,16,17 For example, materials coupling ferromagnetism and ferroelectricity display electric fieldinduced magnetization and magnetic field-induced electric polarization. Such properties are known as the magnetoelectric effect (ME effect), which has been reported in inorganic perovskites and organic–inorganic hybrid materials. 18–23 The ME effect enables potential applications in high-efficiency memories due to controllable magnetic properties by the electric field without energy dissipation. 24 The correlation between magnetic and elastic orders is known as the magnetoelastic effect (MA effect). However, observations of the MA effect are extremely challenging and have been sparsely investigated. 25–27 The MA effect gives rise to the possibility of applications in strain-assisted logic memory and magnetic position sensors. 28,29 Recently, the coupling of crystal a Chemistry Program, Graduate School of Advanced Science and Engineering, Hiroshima University, 1-3-1 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-8526, Japan. E-mail: [email protected] b Chirality Research Center (CResCent), Hiroshima University, 1-3-1 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-8526, Japan c International Institute for Sustainability with Knotted Chiral Meta Matter (WPISKCM 2 ), Hiroshima University, 1-3-1 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-8526, Japan d Precursory Research for Embryonic Science and Technology (PRESTO), Japan Science and Technology Agency (JST), 4-1-8, Honcho, Kawaguchi, Saitama 332-0012, Japan e Arago ´n Nanoscience and Materials Institute (CSIC – University of Zaragoza) and Physics Condensed Matter Dept., C/Pedro Cerbuna 12, 50009 Zaragoza, Spain f Graduate School of Science, Osaka Metropolitan University, 1-1, Gakuen-chou, Naka-ku, Sakai, Osaka 599-8531, Japan g Institut Laue-Langevin, 71 avenue des Martyrs, CS 20156, 38042 Grenoble Cedex 9, France †Electronic supplementary information (ESI) available: Thermogravimetricdifferential thermal analysis, differential scanning calorimetry measurement, stress test, magnetic properties, neutron diffraction, and crystallographic data. CCDC 2411272–2411280. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4tc04445b Received 17th October 2024, Accepted 17th January 2025 DOI: 10.1039/d4tc04445b rsc.li/materials-c Journal of Materials Chemistry C PAPER Open Access Article. Published on 29 January 2025. Downloaded on 5/19/2025 9:58:11 AM. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence. View Article Online View Journal | View Issue 2662 | J. Mater. Chem. C, 2025, 13, 2661–2672 This journal is © The Royal Society of Chemistry 2025 symmetry and spin structures has been observed in chiral magnets. 30 A crystal without an inversion center generates an asymmetry of the orbital angular momentum due to the nonsymmetry of the electric field in the solid, which twists each neighboring spin via spin–orbital coupling (SOC). 31 This is called the Dzyaloshinskii–Moriya interaction (DMI), 32,33 which is the key to chiral magnetism. Therefore, by extension, the MA effect should be observed in multiferroic materials, given the ability of coherent coupling of crystal symmetry and spin structures in chiral magnets. We previously reported that the two-dimensional organic– inorganic perovskite (C 6 H 5 C 2 H 4 NH 3 ) 2 FeCl 4 (PEA-Fe), composed of a 2-phenylethylammonium cation (PEA) and an FeCl 42 anion, is a material with ferroelasticity and canted antiferromagnetism (CAF). 14 PEA-Fe underwent a successive structural phase transitions with rotation of the PEA and tilting of the FeCl 6 octahedron (Fig. 1). The reported crystal structures for each phase of PEA-Fe are I4/mmm above 433 K, Bbcm between 433 K and 323 K, and Pbca from 323 K to at least 90 K, with a ferroelastic phase transition from I4/mmm to Bbcm, clearly identifiable by the observation of ferroelastic domains. The magnetic properties of PEA-Fe exhibited CAF below 98 K (T N ). Although ferroelastic and CAF in PEA-Fe were established in previous work, the crystal structure below T N and the magnetic structure have not been extensively discussed. To investigate the crystal and magnetic structures of two-dimensional organic–inorganic perovskites at liquid helium temperature, neutron diffraction experiments are a powerful tool. 34 In this paper, we evaluated the structural and magnetic properties of two-dimensional organic–inorganic perovskites PEA-M(M=Mn, Cu, Fe) by using magnetometry and powder neutron diffraction (PND) techniques to obtain insight into the MA effect. PEA-Mn and PEA-Cu were synthesized to further discuss the MA effect of PEA-Fe. The three compounds are isostructural and show ferroelasticity and ferromagnetism (M= Cu) or CAF (M=Mn,Fe).In PEA-Fe, we discovered a magnetization shift after cooling in a magnetic field from above T N , and the crystal and magnetic structures down to 2 K were determined by PND measurements. Several models to explain the origin of the magnetization shift will be discussed. The most plausible mechanism is likely associated with the coupling between ferroelasticity and magnetic order, which is similar to the mechanism of chiralityassisted direct coupling between the lattice and magnetic degrees of freedom in chiral magnets. Experimental Crystallization All reagents and solvents were used as purchased. (C 6 H 5 C 2 H 4 NH 3 ) 2 FeCl 4 crystals were prepared as reported previously. 14 The yield was 22%. Elemental analysis: calc. (%) for C 16 H 24 N 2 FeCl 4 : C, 43.47; N, 6.34; H, 5.47. Found: C, 43.66; N, 6.34; H, 5.35. (C 6 H 5 C 2 H 4 NH 3 ) 2 MnCl 4 crystals were prepared using a slow evaporation method. PEACl and MnCl 2 4H 2 O were dissolved in methanol according to the molar ratio. After being kept open to air for several days, pale-pink transparent plate-like single crystals were obtained, and the yield was 72%. Elemental analysis: calc. (%) for C 16 H 24 N 2 MnCl 4 : C, 43.56; N, 6.35; H, 5.48. Found: C, 43.68; N, 6.31; H, 5.49. (C 6 H 5 C 2 H 4 NH 3 ) 2 CuCl 4 crystals were prepared by a similar method to that of (C 6 H 5 C 2 H 4 NH 3 ) 2 MnCl 4 . PEACl and CuCl 2  2H 2 O were dissolved in distilled water according to the molar ratio. After several days of slow evaporation, yellow transparent plate-like single crystals were obtained, and the yield was 82%. Elemental analysis: calc. (%) for C 16 H 24 N 2 CuCl 4 : C, 42.73; N, 6.23; H, 5.38. Found: C, 42.56; N, 6.12; H, 5.35. Characterization CHN spectroscopy was carried out using a PerkinElmer series II CHNS/O Analyzer 2400 or an Exeter Analytical series CE440 Analyzer. Thermogravimetry and differential thermal analysis (TG-DTA) were performed on a Seiko Instruments SII Exstar TG/ DTA 6200 with a temperature range of 293–673 K and a heating rate of 5 K min 1 . Differential scanning calorimetry (DSC) measurements were carried out on a Rigaku Thermo plus DSC8230 with a temperature range of 273–453 K and a scanning rate of 5 K min 1 . The phase transition temperature was determined from the beginning of the thermal anomaly. Crystal images under polarized light were collected using a Meiji Techno EMZ-5HPOL2 polarized microscope to evaluate the ferroelastic behavior. X-ray crystallography Unit-cell determinations were performed on a Bruker D8QUEST equipped with a CMOS area detector or a Rigaku Fig. 1 Comparison of the (a)–(c) side and (d)–(f) top views of the crystal structures of PEA-M(M= Mn, Cu, Fe) in (a) and (d) Pbca, (b) and (d) Bbcm, and (c) and (f) I4/mmm space groups. MCl 6 units show polyhedra. H atoms are omitted for clarify. Color code: orange, M; gray, C; blue, N; green, Cl. Paper Journal of Materials Chemistry C Open Access Article. Published on 29 January 2025. Downloaded on 5/19/2025 9:58:11 AM. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence. View Article Online This journal is © The Royal Society of Chemistry 2025 J. Mater. Chem. C, 2025, 13, 2661–2672 | 2663 XtaLAB Synergy-DW equipped with a HyPix diffractometer. Both employed graphite-monochromated Mo K a radiation (l= 0.71073 Å). Using Olex2, the structures were solved with the SHELXS or the SHELXT structure solution programs and refined with the SHELXL refinement package. 35–38 The nonhydrogen atoms were refined with anisotropic thermal parameters, and hydrogen atoms were added and refined using a riding model. Drawings of crystal structures were performed using the VESTA program. 39 Stress tests The stress tests were performed as reported previously. 14 A single crystal was sandwiched between two stainless plates and pressurized by loading the metal in the tube. The stress application device was heated to higher temperature in an oven (AVO-310V, ETTAS). Magnetic measurements Dc magnetic susceptibilities were collected using a MPMS-5S, MPMS-7, MPMS-XL7 or MPMS3 superconducting quantum interference device magnetometer (Quantum Design) with temperature and dc field ranges of 2–300 K and 50 to +50 kOe, respectively, for single crystals and powder samples. Powder samples were fixed on gelatin capsules after mixing with a small amount of n-eicosane, while single crystals were fixed on gelatin capsules with quartz cotton. The n-eicosane was added to prevent crystallite torquing during sweeping the magnetic field. The quartz cotton was used to avoid stress due to thermal expansion and contraction of adhesive during heating and cooling processes. Single crystals were indexed using X-ray diffraction prior to or after the measurements. Diamagnetic contributions of the sample holder and sample were corrected by the measurement of the sample holder and calculation of Pascal’s constants, respectively. 40 Temperature dependences of the magnetization were measured in the zero-field-cooled warming (ZFCW) and field-cooled cooling (FCC) processes. In the ZFCW process, the sample was first cooled from 150 K to 2 K in the absence of an external magnetic field, and then measurements were made with increasing temperature under the fixed magnetic field. In the FCC process, the magnetization was recorded in the presence of a fixed applied magnetic field with decreasing temperature to 2 K. Fieldsweep measurements were done after the ZFC and FC processes. For the ZFC process, the sample was cooled from 150 to the target temperature under a zero magnetic field, while for the FC process the sample was cooled under a desired magnetic field. For example, the FC process under the 50 kOe magnetic field is referred to as FC 50kOe . Powder neutron diffraction (PND) experiments PND measurements on PEA-Fe on non-deuterated samples were performedatD20andD2B,highfluxandhigh-resolutiondiffractometers, respectively, at the Institut Laue-Langevin (ILL) in France. Powder samples were introduced in a cylindrical vanadium holder of 8 mm in diameter. PND patterns at D20 were collected with l= 2.41 Å at fixed temperatures of 2, 50, 120, 300, 400 and 460 K using a cryo-furnace. Also, thermo-diffractograms were collected, in each temperature ramp between each of these temperatures (thermodiffractograms). High resolution diffractograms were collected at D2B with l= 1.594 Å for different fixed temperatures of 10, 300, 400 and 450 K. Crystalline and magnetic structural parameters at each temperature were determined by Rietveld refinement using the FullProf suite. 41 The raw diffraction data at each temperature are provided in Fig. S1 (ESI†). Results and discussion Structural phase transitions From TG of PEA-Mn and PEA-Cu, weight loss starting from 534 K and 497 K was observed, respectively (Fig. S2, ESI†). DSC curvesofPEA-MnandPEA-Cushowedtwoendothermiceventsat 364 K and 415 K, and 340 K and 408 K, respectively, for the heating process. Two exothermic events for PEA-Mn and PEA-Cu were observed at 368 K and 418 K, and 334 K and 409 K, respectively, for the cooling process. The enthalpy changes DH and the corresponding of the entropy changes DSfor PEA-M were obtained as determined from the area under the heat flow vs.the temperature curve and the equation of DS=DH/T,respectively (Table S1, ESI†). According to the Boltzmann equation of DS= Rln(N), where Ris the gas constant and Nis the variation of the number of disorder positions during the transition, the Nvalues for the low-temperature and high-temperature sides of PEA-Mn and PEA-Cu were calculated to be 2.5 and 1.8, and 1.5 and 2.0, respectively. DSC measurements indicated that two structural phase transitions occurred in PEA-Mn and PEA-Cu. To further investigate the structural phase transitions, single crystal X-ray structural analyses were performed at 293 K, 393 K, and 433 K for PEA-Mn and PEA-Cu (Tables S2 and S3, ESI†). At 293 K, both compounds crystallized in the Pbca space group where 2 PEA units and the corner-shared inorganic layers of MCl 42 (M= Cu, Mn) were alternately stacked along the c-axis. The bond lengths between Mn and the bridging Cl ions in MnCl 6 octahedra were 2.5747(8) and 2.5763(8) Å and that between Mn and the nonbridging Cl ions was 2.4827(9) Å. The bond distances between Cu and the bridging Cl ions in CuCl 6 octahedra were 2.2853(10) and 2.9004(10) Å, and the bond distance between Cu and the nonbridging Cl ions was 2.2950(13) Å. Such M–Cl bond lengths are consistent with those reported values in PEA-Mn and PEA-Cu. 42,43 Upon heating, both compounds were kept in an orthorhombic system but had a Bbcm space group. The C 2 H 4 NH 3 group of the PEA was disordered in two crystallographic equivalent positions. At the highest temperature, the crystal structure changed to a tetragonal system of an I4/mmm space group. In this state, the PEA was fully disordered giving a four-fold symmetry along the c-axis. PEA-Mn and PEA-Cu underwent the same structural phase transformations as PEA-Fe with the increasing temperature (Table S4, ESI†). According to the DSC data, Nwas close to two for both structural phase transitions. It indicates that the transition is close to a simple 2-fold order–disorder model. On the basis of our crystallographic data, the energetically equivalent positions of the PEA groups are changed from one to two then to four Journal of Materials Chemistry C Paper Open Access Article. Published on 29 January 2025. Downloaded on 5/19/2025 9:58:11 AM. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence. View Article Online 2664 | J. Mater. Chem. C, 2025, 13, 2661–2672 This journal is © The Royal Society of Chemistry 2025 during consecutive structural phase transitions, allowing plausible Nvalues. Observations of ferroelastic domains The structural phase transition from I4/mmm to Bbcm is a ferroelastic phase transition with an Aizu notation of 4/ mmmFmmm. 44 The symmetry element is divided by two through the ferroelastic phase transition, from 16 symmetry elements for I4/mmm to 8 for Bbcm, indicating the appearance of two possible orientation states in the ferroelastic phase. From a microscopic point of view, the formation of ferroelastic domains is one of the effective methods to confirm the ferroelastic behavior. Thus, observations of ferroelastic domains for PEA-Mn and PEA-Cu crystals were carried out by polarizing microscopy along the c-axis (Fig. 2). No ferroelastic domain structures were observed at room temperature in fresh PEA-Mn and PEA-Cu crystals at room temperature. However, when the crystal structure was changed to the I4/mmm space group by heating to 430 K and then cooled to room temperature, the linear ferroelastic domain structures appeared. The ferroelastic domain walls of PEA-Mn and PEA-Cu were parallel and orthogonal to the (110) directions. The appearance of domain structures for both PEA-Mn and PEA-Cu means that these compounds exhibit ferroelasticity, similarly to PEA-Fe. For the ferroelastic state with the species 4/mmmFmmm, there are two orientation states of S 1 and S 2 . Then the spontaneous strain tensor ein S 1 and S 2 states is, respectively, 45 eS1 ðÞ¼ e00 0e0 000 2 6 6 6 4 3 7 7 7 5 ;eS2 ðÞ¼ e00 0e0 000 2 6 6 6 4 3 7 7 7 5 where eis the element of the e. The magnitude of eis calculated using the lattice parameters of aand bin the ferroelastic phase. e¼ab aþb  With respect to the lattice parameters in the Bbcm phase for PEA-M(M= Mn, Cu, Fe), the magnitudes of eare 1.5 10 3 , 3.0 10 3 and 2.6 10 3 , respectively. The spontaneous strain e s is defined as the following equation. ðesÞ2¼X 3 i¼1X 3 j¼1 eij  2 As a result, the e s values for PEA-M(M= Mn, Cu, Fe) are 2.1  10 3 , 4.3 10 3 and 3.6 10 3 , respectively. To further investigate the ferroelastic properties of PEA-Mn and PEA-Cu, stress tests were performed on the multi-domain ferroelastic crystal (Fig. S3, ESI†). By application of stress at room temperature, the ferroelastic domains of PEA-Mn were moved, while for PEA-Cu, the crystals broke before observation of domain wall motion, similar to PEA-Fe. 14 After thermal treatment up to 373 K under stress, the ferroelastic domains of PEA-Cu were moved. Therefore, the stress tests for both PEAMn and PEA-Cu confirmed the ferroelastic behavior of these compounds. In contrast to PEA-Cu and PEA-Fe, the ferroelastic domains of PEA-Mn are mobile at room temperature, which can be attributed to the larger spontaneous strain e s for PEA-Cu and PEA-Fe compared to PEA-Mn. Magnetic properties of PEA-M PEA-Mn showed a magnetic transition at 44.3 K with a rapid increase of the magnetization when the temperature decreased and saturation at low temperatures (Fig. S4, ESI†). The magnetic transition temperature was consistent with that in previous study. 42 Our experimental results included the magnetic properties of PEA-Mn along all crystallographic axes. The ZFCW and FCC processes exhibited similar magnetization, with values at 2 K for the FCC process of 4.2 emu Oe mol 1 and 6.2 emu Oe mol 1 , along the band c-axes, respectively. However, along the a-axis, the FCC magnetization reached a value of 31.70 emu Oe mol 1 at 2 K whereas a ZFCW magnetization of 0.48 emu Oe mol 1 at 2 K and a cusp reaching 12 emu Oe mol 1 at 44.0 K were observed. A magnetic hysteresis with a remanent magnetization (M rem )of510 3 m B and a coercive field of 260 Oe was observed when the magnetic field was applied along the a-axis. In the high field region, a linear increase of magnetization was observed, reaching a value of 0.18 m B at 50 kOe. This value is much smaller than the theoretical saturated magnetization value of M sat =gSm B =5m B (Mn 2+ ,g=2,S= 5/2). Application of a magnetic field along the b-axis exhibited a linear increase of magnetization, reaching a value of 0.17 m B at 50 kOe. Along the c-axis, spin-flop transition was observed around 35 kOe, which was in agreement with the reported study. 42 During this spin-flop transition, the easy Fig. 2 Polarized optical microscopy images of different ferroelastic domains in single crystals of (a) and (b) PEA-Mn and (c) and (d) PEA-Cu at room temperature before (a) and (c) and after (b) and (d) heating. Scale bar: 0.5 mm. Paper Journal of Materials Chemistry C Open Access Article. Published on 29 January 2025. Downloaded on 5/19/2025 9:58:11 AM. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence. View Article Online This journal is © The Royal Society of Chemistry 2025 J. Mater. Chem. C, 2025, 13, 2661–2672 | 2665 axis of the antiferromagnetic (AFM) arrangement was altered from the c-axis to the direction perpendicular to the applied field. For PEA-Cu, observation for the ZFCW and FCC processes was similar (Fig. S4, ESI†). Along all axes, abrupt changes in the magnetization were observed around 9.5 K demonstrating a long-range magnetic order, in agreement with the previous studies. 43,46,47 The magnetic field dependences along all axes were described as soft magnetic behavior with a saturated magnetization of around 1.1 m B at 50 kOe, and hysteresis could not be observed in spite of the measurements at a field step of 5 Oe. The experimental saturation value is in good agreement with a value of M sat =1m B (Cu 2+ ,g=2,S= 1/2). For PEA-Fe, a magnetic transition at 98 K with a rapid increase of the magnetization along all axes in both ZFCW and FCC processes was observed (Fig. S4, ESI†). This transition temperature was in good agreement with the reported one. 14 A magnetic hysteresis at 5 K with a M rem of 0.02 m B and a coercive field of 5 kOe was observed when the magnetic field was applied along the a-axis. At 50 kOe, the magnetization reached a value of 0.10 m B , which was much smaller than the saturation value of gSm B = 6.84m B (Fe 2+ ,g= 3.42, 48 S= 2). Application of the magnetic field along band c-axes exhibited a linear increase of magnetization, reaching values of 0.01 m B and 0.08 m B at 50 kOe, respectively. Results from temperature and magnetic field dependency measurements suggested that PEA-Mn and PEA-Fe display CAF with a spin-canted angle yestimated at 0.051and 0.171, respectively, using M rem =M sat siny. 49 In comparison to PEAMn and PEA-Fe, PEA-Cu displays ferromagnetic (FM) ordering below 9.5 K. Magnetization shift of PEA-M The field-sweep measurements after the ZFC process exhibited a symmetric magnetization loop with respect to the origin for PEA-Mn and PEA-Cu. For PEA-Fe, the symmetric magnetization loop was observed along band c-axes while not along the aaxis. The shift of the magnetization curve after ZFC reflects on the weak cooling field remaining in the system. In fact, this shift is dependent on the cooling magnetic field (see ‘‘Relationship between SOC and the value of magnetization shift’’ section for details). Therefore, the application of exact ZFC would introduce the absence of the magnetization shift along the a-axis. After the FC process, a similar behavior to that of the ZFC process was observed for PEA-Mn and PEA-Cu (Fig. 3 and Fig. S5, ESI†). In the case of PEA-Fe, a displacement of the magnetization curves was observed, along all axes, after the FC process. To quantify the displacements of the magnetization curve, a shift value (M Shift ) can be defined as (M + +M  )/2, where M + and M  are the magnetization at +50 and 50 kOe magnetic fields, respectively. Additionally, a shifted field (H Shift ) is calculated using the equation H Shift =(H + +H  )/2, where H + and H  correspond to the upper and lower magnetic fields of the magnetization zero point. The values of H Shift during the FC 50kOe process were 810 kOe along aand b-axes and only 80.5 kOe along the c-axis. The M Shift values after FC 50kOe along the a-, band c-axes were 2.04 10 2 m B ,2.7 10 3 m B Fig. 3 Field dependence of the magnetization of PEA-M, with the magnetic field applied along (top) a-, (middle) band (bottom) c-axes at 5 K after ZFC, FC +50kOe and FC 50kOe . (a) M= Mn, (b) M= Cu, and (c) M= Fe. Journal of Materials Chemistry C Paper Open Access Article. Published on 29 January 2025. Downloaded on 5/19/2025 9:58:11 AM. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence. View Article Online 2666 | J. Mater. Chem. C, 2025, 13, 2661–2672 This journal is © The Royal Society of Chemistry 2025 and 0.8 10 3 m B , respectively (Fig. 4). For powder samples, the magnetization shift was also observed for PEA-Fe after FC 50kOe while not for PEA-Mn and PEA-Cu (Fig. 4 and Fig. S6, ESI†). The field-sweep measurements were also measured at 105 K where PEA-Fe showed paramagnetic, confirming the lack of the magnetization shift (Fig. S7, ESI†). The M Shift ,H Shift ,H + and H  for PEA-Fe along the a-axis were plotted as a function of the temperature (Fig. S8, ESI†). The H + and H  became presence below T N while the M Shift and H Shift were observable at least blow 30 K. The increase in M Shift was accompanied by an improvement in the H Shift . PND study of PEA-Fe The thermo-diffractograms collected at D20 for PEA-Fe as the system was heated from 2 K to 460 K allow two structural phase transitions to be observed (Fig. 5). The first one appeared around 340(10) K, where some peaks merged as the temperature was increased and a second one was visible at around 432(5) K. Although our data do not allow a high-quality refinement due to the presence of the incoherent scattering of the H atom, and a high degree of disorder at high temperatures, the space group in each phase has been identified, which corresponds to the Pbca for 2 K oTo343 K, Bbcm for Tin the interval of 343 K oTo433 K and I4/mmm for 433 K oT,in agreement with the anomalies observed in DSC measurements. The analysis of the thermo-diffractograms allows the thermal evolution of the lattice parameters to be determined as a function of the temperature (Fig. S9, ESI†). The change from tetragonal to orthorhombic was clearly visible at 433 K whereas the change from Bbcm to Pbca was noticed at 343 K as a drastic Fig. 4 (a) The values of the magnetization shift for single crystals and powder samples of PEA-M(Mn, Cu and Fe) at 5 K after ZFC, FC +50kOe and FC 50kOe . (b) The enlarged figure of (a). The gray bidirectional arrows present in both graphs correspond to the same range. Fig. 5 (a) Thermo-diffractograms measured warming up the PEA-Fe from 2 K to 460 K. (b) and (c) Neutron diffraction patterns of the regions (b) 2.2 o Qo2.6 Å 1 and (c) 3.3 oQo3.7 Å 1 at 120, 300, 400 and 460 K, corresponding to green, yellow, orange and red double-headed arrows in panel (a), respectively. In Table S5 (ESI†), the refinement parameters for different temperatures in the Pbca space group phase below 343 K are given. Paper Journal of Materials Chemistry C Open Access Article. Published on 29 January 2025. Downloaded on 5/19/2025 9:58:11 AM. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence. View Article Online This journal is © The Royal Society of Chemistry 2025 J. Mater. Chem. C, 2025, 13, 2661–2672 | 2667 change in the aand b-axes. The unit cell volume is decreasing progressively, but not linearly, as temperature decreases and shows two small anomalies at the critical temperatures. The high statistics data collected at fixed temperatures on D20 confirmed that there was no structural phase transition in the range 2 K to 343 K. However, some very small additional peaks were observed at the diffractograms at 2 K and 50 K around Q= 0.9, 1.8 and 2 Å 1 that were not present at the diffractograms measured for T4100Kandthatwereindexed with a propagation vector - k= (0,0,0) (Fig. 6a). This suggested the onset of a long-range magnetic order, which was estimated to appear around T= 97(7) K from the thermo-diffractograms data. PND measurements of PEA-Fe revealed the absence of the structural phase transition from the ferroelastic Pbca to a higher symmetry space group, suggesting the persistence of ferroelasticity below T N . The remaining of the ferroelastic phase below the magnetic phase transition temperatures in PEA-Mn and PEA-Cu can be assumed, due to their structural similarity with PEA-Fe. Magnetic structure analysis of PEA-Fe To describe the symmetry of the magnetic phase, the irreducible representation theory is employed. For the Fe atom in Wyckoff position 4b of the Pbca space group and propagation vector - k= (0,0,0), the magnetic representation was decomposed as the direct sum of irreducible representations (Irreps) as follows: G M =3mG + 1 (1) + 3mG + 2 (1) + 3mG + 3 (1) + 3mG + 4 (1) (1) which indicated that each 1-dimensional Irreps appears 3 times. In Table 1, the Fourier coefficient for each Irreps was shown. In all the above Irreps, the magnetic moments were allowed to have components (u,v,w) along the a-, b-, and c-axes. From the magnetization data, we know that the magnetic structure has to be a CAF along the a-axis and the only Irreps that allows a FM component along the a-axis and AFM along the others is G + 4 . The fit of the high statistics diffractogram at 2 K was done with each one of the previous Irreps, and, as was Fig. 6 (a) and (b) Refinements of the diffractogram for PEA-Fe at (a) 120 K and (b) 2 K. The red points are the experimental points, the solid black line is the fit, the blue line is the difference between the experimental data and the fit, and the green lines correspond to the hkl positions for the nuclear and magnetic phases. In the inset of panel (a), the subtraction between the high statistics (more than 4 hours) diffractograms collected at 2 K and 120 K in D20 is shown in order to better visualize those small magnetic peaks. The inset in panel (b) represents the fitting curve by using the model Pb0c0a. (c) Magnetic structure of PEA-Fe. Red and blue lines show the (blue) ferromagnetic J FM and (red) antiferromagnetic J AFM interactions between neigboring magnetic ions. Yellow arrows indicate the orientation of the magnetic dipoles. Journal of Materials Chemistry C Paper Open Access Article. Published on 29 January 2025. Downloaded on 5/19/2025 9:58:11 AM. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence. View Article Online 2668 | J. Mater. Chem. C, 2025, 13, 2661–2672 This journal is © The Royal Society of Chemistry 2025 expected, only the G + 4 was able to fit the peaks shown in the inset of Fig. 6a giving a value of the magnetic moment of 2.6(2) m B along the b-axis and zero along the other directions. In fact, the other possible models allowed by the symmetry, described by G + 1 ,G + 2 and G + 3 , failed to reproduce the intensity found experimentally at the magnetic Bragg peak (1 0 1). At this point, it is important to remark that the neutron diffraction experiments are done under zero magnetic field and that the measured magnetic signal, in spite of the high acquisition time, was very small. Therefore, we can conclude that the magnetic structure is the one given by the G + 4 irreducible representation, where the Fe magnetic moments form a CAF structure along the a-axis. In fact, the structure obtained from the fit is a pure AFM along the b-axis, with no component along the a-axis. However, this is in agreement with our magnetization data, since it is mentioned that the component parallel to the a-axis is around 0.02 m B , with a canting angle of around 0.171. In this case, such a low magnetic component will be undetectable with PND, even with high acquisition times. The magnetic subgroup which corresponds with the fitted magnetic structure is Pb0c0a(No. 61.4.500), with a transformation matrix given by: (b,c,a; 0, 0, 0). The fits of the diffractograms at 120 K in the paramagnetic phase and 2 K in the ordered phase are shown in Fig. 6, respectively, with continuous black lines. The experimental data are the red points whereas the blue line represents the difference between the experimental and the fitted model at each temperature. The most intense magnetic peak is only correctly fitted by considering the model Pb0c0a(the inset in Fig. 6b). Relationship between SOC and the value of the magnetization shift Comparing the field cooling effect for PEA-M, only PEA-Fe shows the shift in both the powder samples and single crystals. Such magnetic behavior has been observed in the iron-based 2D perovskite compound (C 2 H 5 NH 3 ) 2 FeCl 4 . 12 The authors tentatively associated the magnetization shift with the interaction of magnetic, electric, and elastic domains. In this study, we considered three possible mechanisms to explain the magnetization shift: a minor loop effect, an exchange bias effect and the interaction between magnetic and elastic orders. The first possible mechanism to explain the magnetization shift involves minor loops. Minor loops are typically defined as a hysteresis loop without saturation within the applied magnetic field. After cooling under the magnetic field, minor loops appear as shifted hysteresis loops when the maximal external field is too small to completely reverse the magnetization. Saturation of the magnetization of PEA-Fe within an external field of 50 kOe was not achieved, making it difficult to distinguish between magnetization shifts from the minor loop phenomena. To do so, the field dependence magnetization was performed at 5 K after various field cooling (Fig. S10a, ESI†). The value of the magnetization shift was increased and reached in saturation by increasing the cooling field (Fig. S10b, ESI†). The magnetization shift is observed even when the external field is over 1000 times larger than the cooling field, and magnetic behaviors under the field sweep cannot be minor loops. Therefore, the scenario based on the minor loop phenomenon is not suitable. The second plausible mechanism to explain the magnetization shift involves the coexisting (anti)ferromagnetic and spin glass orders. Such a property is known as the exchange bias. 50 The exchange bias has been observed in various systems including FM or ferrimagnetic nanoparticles embedded in an AFM matrix, FM/AFM thin film heterostructures and materials with coexistence of (anti)ferromagnetic and spin-glass orders. 50–52 The theory for the exchange bias in their systems has been described asvariousmodelssuchastheuncompensatedspins,andthe pinning of domain walls at the interface. 53,54 Many materials with the exchange bias exhibit the training effect, i.e., a reduction of the displacement value of the magnetization curve upon repeated measurements without new field cooling. The origin of the training effect is explained as a rearrangement of the AFM domain structures with repetition of field cycles. 55,56 For the powder sample of PEA-Fe, we repeated the magnetic field-sweep measurements six times at 5 K after field cooling under the 50 kOe magnetic field (Fig. S11, ESI†). The magnetization shift was invariant of the number of field cycles, suggesting that the second mechanism also cannot describe the magnetization shift. The third plausible mechanism to explain the magnetization shift involves the coupling of different ferroic orders. In contrast to the second mechanism, where the behavior is attributed to the interaction between different spin orders, here the behavior is attributed to the interaction between different ferroic orders. ME materials show an electrically controllable hysteresis loops shift, which is attributed to the coupling between magnetic and electric orders. 57,58 Materials with the ME effect tend to have an incommensurate magnetic structure; however, PND experiments on PEA-Fe revealed a commensurate magnetic structure with the propagation vector - k= (0,0,0). In addition, the crystal structure of PEA-Fe below T N is nonpolar (Pbca). From this finding, we assume the absence of the Table 1 Irreps contained in the magnetic representation G M for the Fe atom in the Wyckoff position 4b of the Pbca space group and the propagation vector k -= (0,0,0) Fe1 Fe2 Fe3 Fe4 x,y,zx+ 1/2, y,z1/2 x,y+ 1/2, z+ 1/2 x+ 1/2, y+ 1/2, z+1 G + 1 (u,v,w)(u,v,w)(u,v,w)(u,v,w) G + 2 (u,v,w)(u,v,w)(u,v,w)(u,v,w) G + 3 (u,v,w)(u,v,w)(u,v,w)(u,v,w) G + 4 (u,v,w)(u,v,w)(u,v,w)(u,v,w) Paper Journal of Materials Chemistry C Open Access Article. Published on 29 January 2025. Downloaded on 5/19/2025 9:58:11 AM. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence. View Article Online