Structural Investigation into Magnetic Spin Orders of a Manganese Phosphatic Oxyhydroxide, Mn 5 [(PO 4 ) 2 (PO 3 (OH)) 2 ](HOH) 4

: The ferri- and antiferromagnetic structures of a hureaulite-type synthetic compound, Mn 2+5 (PO 4 ) 2 (PO 3 (OH)) 2 (HOH) 4 , were elucidated by high-resolution neutron powder diffraction in combination with magnetic susceptibility and heat capacity measurements. At 6.17 K, the paramagnetic phase (space group: C 2/c) transforms to inherit a ferrimagnetic order (magnetic space group: C 2 (cid:48) /c (cid:48) ), followed at 1.86 K by an incommensurately modulated antiferromagnetic order (magnetic superspace group: P 2 1 /c.1 (cid:48) ( α 0 γ ) 00s with the propagation vector k (0.523(2), 0, 0.055(1)). In the ferrimagnetic state, antiferromagnetic interactions are dominant for both intra and inter pentamers of Mn 2+ (O, HOH) 6 octahedra. Differently aligned spin-canting sublattices seen in the ferrimagnetic models at 3.4, 4.5, and 6.1 K explain a weak ferromagnetism in the title compound. The observation of magnetic moments vigorously changing in a small temperature range of 6.1–1.5 K adumbrates a high complexity of interplaying structural and magnetic orders in this manganese phosphatic oxyhydroxide.


Introduction
Transition metal phosphatic oxyhydroxides comprise an ample class of dense and porous framework structures occurring in nature [1,2]. Their variety and complexity are directly associated with the presence of (OH) − and (HOH) 0 groups (non-zeolite water) responsible for various configurations of hydrogen bonds (HBs). These (OH) − and (HOH) 0 ligands along with O 2− form tetrahedra and octahedra, where P 5+ and transition metal (TM) cations occupy those respective framework polyhedral centers. With respect to structure-related properties of this material class, fundamental studies of charge transport and proton tunneling over HB networks are application-relevant research topics attracting our interest [3][4][5][6][7]. At the same time, it is supposable to find new (multi)ferroic and magnetic structures depending on the type and arrangement of 3d TM atoms in phosphatic oxyhydroxides. However, up to date, their magnetic structures are either completely unknown or have been poorly investigated. As a systematic approach to unveiling magnetic spin orders in this interesting material system, new antiferromagnetic (AFM) and ferrimagnetic (FM) structures of eosphorite-and rockbridgeite-type oxyhydroxides, respectively, could be determined in our recent studies using neutron diffraction techniques [7,8]. The current Previous Mössbauer spectroscopic studies on magnetic behavior of the hureaulite-type series M 2+ 5 (PO 4 ) 2 ((PO 3 (OH)) 2 (HOH) 4 with M = Mn, Fe, and Co indicated antiferromagnetism in Fe-and Co-rich solid solutions but a tendency to ferromagnetism in Mn-rich members [11][12][13]. In the current study, results from neutron powder diffraction studies of Mn-hureaulite explicate strong temperature-dependent FM spin arrangements below the Curie temperature (T C = 6.17 K). In addition, we present its low-temperature magnetic phase, showing an incommensurately modulated AFM order stable below 1.86 K.

Sample Preparation and Phase Identification
Details of hydrothermal synthesis conditions for a high yield of pure and highly crystalline Mn-hureaulite (protonated, i.e., H-form) can be found in our previous report [3]. The typical deuterated (D-form) sample preparation for neutron powder diffraction is as follows: the mixture of 0.1475 g of SrCO 3 (Merck Technipur, Kenilworth, NJ, USA, ≥99.9%), 0.2516 g of MnCl 2 (Sigma Aldrich, St. Louis, MO, USA, ≥99%), and 0.22517 g of H 3 PO 4 (Merck, 85 wt.%) was dissolved in 9 mL of D 2 O (VWR, Radnor, PA, USA, ≥99.90%) and aged while stirring for one full day. This colloidal solution was charged in a 23 mL Teflon-lined stainless-steel autoclave (Parr Instruments) and allowed to react at 453 K in a convection oven (Nabertherm, Lilienthal, Germany, N30/45HA) for a cyclic heating of 2 × 72 h with a short cooling period in between. The synthesis product was washed with heavy water and dried in a vacuum flask at room temperature. From a 10 mL starting solution, a maximum of about 0.1 g product could be yielded. To obtain several grams of samples for high-resolution neutron powder diffraction (HRNPD), the aforementioned synthesis task was repeated in series using several autoclaves. Every synthesis product from each autoclave was separately prepared for phase analyses using X-ray powder diffraction (XPD) on a Bragg-Brentano diffractometer (XRD3003, GE) with MoK α1 radiation (Ge(111)-monochromator). The phase purity and quality could be verified by comparison of XPD data collected with synthesis products from each autoclave to those of pure Mn-hureaulite samples reported in [3].
The molar heat capacity (C p ) was recorded with 1-2 mg of powdered sample using the semiadiabatic relaxation method on a Quantum Design Physical Properties Measurement System (PPMS), equipped with an 3 He refrigerator, in the temperature range of 0.35-300 K. Samples were prepared as pressed powder pellets and fixed to the sample holder with Apiezon N grease. D.C. magnetic susceptibility was acquired at external fields (100 Oe) on another Quantum Design PPMS available for the temperature range of 2-300 K.

High-Resolution Neutron Powder Diffraction (HRNPD)
HRNPD was conducted with two batch samples of the D-form of Mn-hureaulite (Mn 5 (PO 4 ) 2 (PO 3 (OD)) 2 (DOD) 4 ). The data collection using λ = 2.5360(2) Å was subsequently performed on a batch sample of 2.5 g at 6.5 and 3.4 K on the instrument SPODI at the neutron facility FRM II [14]. Further 6.5, 6.1, 4.5, 1.7, and 1.5 K HRNPD datasets were acquired using λ = 2.4500(1) Å on another batch sample of 1.5 g on the instrument HRPT at the neutron facility SINQ at PSI [15]. Rietveld refinements were carried out using the program package Jana2006 [16] with starting refinement parameters obtained from our recent studies [3]. Atomic and magnetic structure models refined were visualized using the program VESTA [17].

Magnetic Phase Transitions
D.C. magnetic susceptibility (χ DC ) of Mn-hureaulite revealed a nonlinear enhancement distinctively from about 6.2 K to near 2.2 K and then a fast drop down to its lowest measuring temperature 2 K (Figure 2a). This indicates two subsequent magnetic phase transitions; near 6.2 K (Curie temperature: T C ), Mn-hureaulite undergoes a magnetic phase transition from a paramagnetic (PM) to an FM and then to an AFM state near 2.0 K (denoted as T 1 ). Its reciprocal values (1/χ DC ) linearly decreasing from 50 K to 20 K confirm the PM state according to the Curie-Weiss law (Figure 2b). A negative Curie-Weiss temperature (T θ ) about −33.5 K, which was estimated by linear regression with 1/χ DC values between 20 K and 50 K, indicates dominant antiferromagnetic interactions in Mnhureaulite. The inverse slope of 1/χ DC against T corresponds to a Curie constant (C) value of 23.85 emu·K·mol −1 , allowing an estimation of the effective magnetic moment (µ eff ) value in Bohr magneton (µ B ). The resulting µ eff value of 13.81 µ B is comparable with the value 13.53 µ B reported in [12] but much larger than the magnetic moment 5.92 µ B of a single isolated Mn 2+ ion (3d 5 , high-spin). Hence, a Mn pentamer may be regarded as the magnetic unit in Mn-hureaulite [12].
A 2% temperature rise needed for a good accuracy of the heat capacity in the applied relaxation technique gave rise to two peaks clearly resolved at 6.17 K and 1.86 K without compromising the peak height and shape. The molar heat capacity C p (J·K −1 ·mol −1 ) was normed by the gas constant R = 8.31446 J·K −1 ·mol −1 . In the resulting C p /R(T) curve (Figure 3a), these two peaks at 6.17 K and at 1.86 K could be unambiguously assigned to the respective temperatures T C and T 1 in accordance with χ DC data. The relative error of C p showed a noticeable increase around T 1 (Figure 3b) when analyzing the error, which included the deviation from an exponential curve fit. On the other hand, there was a small C p anomaly at about 0.4 K, showing a clear and reproducible change of slope at this temperature. However, the incomplete picture of the anomaly at 0.4 K due to experimental limitation below 0.35 K prevented a categorical assignment. Additional C p measurements through the main transitions upon heating and cooling showed no thermal hysteresis within experiment accuracy (Figure 3c,d).  (a) C p /R(T) shows two distinctive discontinuous points T C at 6.17 and T 1 at 1.86 K. An additional small peak at about 0.4 K could denote a further magnetic phase transition at lower temperatures. An extremely broad bump centered around 10 K indicates a short-range magnetic order. (b) A noticeable increase around T 1 of the percentage of the relative error C p = (100 × error(C p )/C p ) is indicative of a first-order transition at T 1 . Additional C p measurements through T 1 (c) and T C (d) upon heating (circles in red) and cooling (squares in black) show no thermal hysteresis within experiment accuracy.
The magnetic contribution C m to C p was estimated (Figure 4a) after subtracting the lattice contribution C L , without considering electron and nuclear contributions insignificant in the measuring temperature range. The lattice contribution term C L /R = 0.01·T 1.75 was estimated from the linear temperature-dependent region of C p in log-log scale between 30 and 70 K. From integrating C m (T) over T, the magnetic heat Q m (T) was obtained for T < 20 K (Figure 4b). The integrated area of Q m (T) around T 1 (1.5-2.2 K) and T C (5.0-6.5 K) gave rise to the respectively associated latent heat Q L . The resulting Q L (T 1 ) = 10 J·mol −1 and Q L (T C ) = 57 J·mol −1 reflected the respective quantity of heat necessary for the AFM → FM and the FM → PM transition of Mn-hureaulite. The total magnetic entropy S m of an Mn pentamer can be calculated by 5 × R·ln(2S + 1) = 5R·ln 6 = 74.45 J·K −1 ·mol −1 with the spin angular moment S = 5/2 for Mn 2+ . The entropy changes at T 1 and T C were obtained as Q L (T 1 )/T 1 = 5.4 J·K −1 ·mol −1 and Q L (T C )/T C = 9.2 J·K −1 ·mol −1 . A small change in magnetic entropy (∆S m = 3.4 J·K −1 ·mol −1 ) from the AFM to the FM state indicates that the AFM structure of Mn-hureaulite is only slightly higher ordered than the FM one. On the other hand, the total magnetic entropy S m (T) evaluated for a Mn pentamer was much less than 5R·ln(6) = 8.96R expected for complete ordering of localized S = 5/2 spins, as shown by the small values 1.1R and 3.67R for S m (T 1 ) and S m (T C ), respectively ( Figure 4c). These point energetically stable magnetic spin arrangements, where a Mn pentamer acts as magnetic unit in this compound. Consequently, a high degree of long-range magnetic spin orders for both inter and intra pentamers was predicted below T C by heat capacity measurements. Interestingly, C p /R(T) between T C and 20 K in Figure 4a exhibits an extremely broad bump centered around 10 K. This observation resembles a broad magnetic heat capacity peak of ludlamite (Fe 3 (PO 4 ) 2 (HOH) 4 ) at about 15 K, which was explained by a linear-chain Ising model considering a short-range order in the AFM coupling between triads of Fe atoms [18].
Hence, there may be spin orientation disorder between pentamers in Mn-hureaulite before entering into the long-range FM order.
The temperature-induced evolution of magnetic spin orders in Mn-hureaulite was obvious when comparing its raw HRNPD datasets subsequently collected at 6.5, 6.1, 4.5, and 1.5 K with a constant wavelength of 2.45 Å ( Figure 5). A weak intensity for the reflection at 2θ = 29.96 • was exclusively from the atomic order at 6.5 K. The FM order contribution to this Bragg peak increased with decreasing temperatures between T C and T 1 , as highlighted by the red arrow in Figure 5. This noticeable variation in a small temperature range from 6.1 K down to 4.5 K adumbrated vivid spin reorientations in the FM state. Below T 1 (e.g., at 1.5 K), a high number of additional reflections of the AFM phase appeared. Figure 5. Comparison of raw HRNPD data of Mn-hureaulite collected in series at 6.5, 6.1, 4.5, and 1.5 K on the instrument HRPT at the neutron facility SINQ at PSI. The FM order contribution to the reflection 310 (red arrow) is strongly enhanced from 6.1 K to 4.5 K, indicating vivid spin reorientations between T C and T 1 . A high number of additional AFM reflections are seen in HRNPD at 1.5 K, particularly the one near 4.5 • (2θ) highlighted by a blue arrow.

PM and FM Structure Models Refined with HRNPD Data (λ = 2.536 Å)
Rietveld analyses of 6.5 K and 3.4 K HRNPD datasets were performed with constraints on atomic displacement parameters (ADPs); isotropic ADPs for five independent deuterium (D) sites could be individually refined, whereas isotropic ADPs for Mn and O sites were only refined group-specifically. ADPs for P sites were fixed to a small value of 0.001 Å 2 as their values as free parameters became negative. In spite of the application of those constrained ADPs (Supplementary Table S1a), the calculated HRNPD pattern without magnetic contribution at 6.5 K in the space group C2/c was well matched to the observed one, as shown in graphical representation of the final cycle of refinements in Figure 6. This confirms the paramagnetic state at 6.5 K under zero external field. All experimental and refined cell parameters are given along with agreement factors in Table 1.
The refined atomic and profile parameters at 6.5 K were used as starting values in the following refinement with 3.4 K HRNPD data. In comparison to the 6.5 K HRNPD pattern, no extra Bragg peaks were detected in 3.4 K HRNPD data, along with a clear change in relative intensity distribution (Figure 7a). This difference was the most distinguishing feature in HRNPD patterns recoded above and below the Curie temperature T C . In our case, as mentioned in the previous section, the reflection 310 at 2θ = 29.96 • (marked as red arrows in Figures 5 and 6) was mostly sensitive to the occurrence of the FM phase of Mn-hureaulite, as its intensity dramatically increased at lowered temperatures below T C .   (2)   As given in Table 2, keeping the (a × b × c) lattice for the magnetic and atomic unit cells (i.e., the propagation vector k(0, 0, 0)), there were four Shubnikov groups as direct magnetic subgroups of the space group C2/c, i.e., C2 /c , C2/c, C2/c , and C2 /c [19]. Among them, C2/c and C2 /c were out of consideration because they gave rise to AFM orders with zero global magnetization (Table 2). Hence, with 3.4 K, HRNPD data magnetic moments at three unique Mn sites (Table 3) were refined for two FM space groups C2 /c and C2/c. The former allowed for two global magnetization components M x and M z at all three Mn sites whereas the latter gave rise to one nonzero global magnetization component, M y only. Based on agreement factors of Rietveld refinements, C2 /c was identified as the correct magnetic space group for the FM phase at 3.4 K. This can be clearly demonstrated by the difference profile flattened between observed and calculated HRNPD patterns of the model with C2 /c (Figure 7a). In contrary, Rietveld calculations for another FM model with C2/c resulted in large residuals in the difference pattern, particularly for both 310 and 021 reflections highlighted by arrows in Figure 7b. For an explicit comparison of C2 /c to C2/c, three magnetic sublattices resulted from Rietveld refinements with both FM space groups are illustrated in Figure 8. Atomic structure parameters refined with HRNPD data at 3.4 K are given in Table 3 (all atomic parameters obtained from HRNPD data analyses in this study are provided in Supplementary Table S1a-e). Table 2. Magnetic subgroups (Shubnikov groups [19]) of the space group C2/c are given with the magnetic moment M(M x , M y , M z ) allowed at the three Mn sites. Their global magnetization terms (given below the backslash) result in four different FM and AFM types. Mn1 is a special position while Mn2 and Mn3 are general positions (see Table 3).   In conclusion, at 3.4 K, Mn-hureaulite exhibits an FM order in C2 /c , and the easy axis lies on the (a-c) plane (Table 4). More interestingly, its magnetic moments vividly change between a and c by a small temperature range between T C and T 1 , as described below.
where i runs over all points of the profile, y i (obs) and y i (calc) are observed and calculated intensity, respectively, n is the number of profile points, p is the number of refined parameters, and w i = 1/σ 2 (y i (obs)).  Results from Rietveld calculations for both C2 /c and C2/c models with HRNPD datasets collected at 4.5 K and 6.1 K could confirm that the FM phase retains the magnetic space group C2 /c from T C down to T 1 (graphical representations of results from Rietveld calculations are shown in Supplementary Figure S2). However, a strong temperaturedependent magnetic spin-canting reorientation was observed in the FM state. At 6.1 K, all three Mn sites exhibit large values for magnetic moments in the a direction, but at 3.4 K in the c direction (Table 4). Exemplarily, for a small temperature difference of ∆T = 1.6 K, an explicit change of magnetic spin orientations is demonstrated in Figure 9a. The FM arrangement in the (a-b) plane is more obvious at 6.1 K than that at 4.5 K (Figure 9b). At lowered temperatures, the FM character parallel to the a axis is abated due to the reorientation of magnetic moments toward the c axis (especially at Mn1), while the AFM arrangement in the (b-c) plane becomes more obvious (Figure 9b). This strongly temperature-dependent reorientation of magnetic moments is graphically conflated in Figure 10.  Overall, the magnetic system of Mn-hureaulite exhibits a strong anisotropy with three differently aligned spin-canting Mn sublattices. Based on the geometry of inter and intra pentamers in the FM phase (Table 5, Figure 11), superexchange paths via Mn-O-Mn >90 • might favor AFM couplings as dominant magnetic interactions for intra and inter pentamers [20]. Hence, spin-canting alignments varying across the Mn1-Mn2-Mn3 sublattices can be regarded as the main factor for a weak ferromagnetism detected by magnetic susceptibility measurements ( Figure 2). Thus, our FM models between T C and T 1 could validate the tendency to ferromagnetism in Mn-rich hureaulite-type solid solutions, predicted by Mössbauer spectroscopic studies [11][12][13]. The refined magnetic moments at three unique Mn sites are given in Table 4.

The AFM Phase with Incommensurately Modulated Magnetic Spins Below T 1
HRNPD data collected with λ = 2.45 Å at 1.5 and 1.7 K exhibit additional Bragg reflections compared to those collected above T 1 (Figures 5 and 12). This features the presence of an AFM phase by the second phase transition. Rietveld analyses were done with both 1.5 K and 1.7 K HRNPD data, showing the same AFM structure, except for a distinctive decrease in the lattice metric from 1.7 K to 1.5 K (Table 1; Supplementary Figure S3). Hence, results from Rietveld analyses with 1.5 K HRNPD data are described below.
With the incommensurate modulation vector k(0.523(2), 0, 0.055(1)) determined at 1.5 K for an AFM order, its lattice might lose the C-centering to be primitive. Accordingly, there are four possible AFM orders in P2 1 /c.1 , P2 1 /n.1 , P2/c.1 , and P2/n.1 . Within the experimental resolution of our neutron diffraction data, we could exclude the presence of main reflections h0lm (m = 0) with l = odd, as well as 0k0m with k = odd. Thus, only P2 1 /c.1 was taken into account to select possible magnetic superspace groups (MSSG), i.e., P2 1 /c.1 (α0γ)00s and P2 1 /c.1 (α0γ)0ss. Note that s in these one-line magnetic superspace group symbols [21] stands for a translation of 1 2 in the fourth symmetry operator x 4 by time inversion 1 involved in the operation (Table 6) [22]. Results from Rietveld calculations suggested the magnetic structure model in P2 1 /c.1 (α0γ)00s with the weighted residual value of magnetic satellites wR p (mag) = 0.033 ( Figure 12) rather than P2 1 /c.1 (α0γ)0ss with wR p (mag) = 0.043. However, the current model in P2 1 /c.1 (α0γ)00s ( Figure 13) was based on incommensurate modulations of magnetic moments refined without consideration of possible modulations in atomic positions and displacements below T 1 . The difference pattern with the refined AFM model was relatively rough in comparison to those for the FM phase. Hence, it is necessary to improve the current AFM model.  (1)). Short bars in black and in green indicate Bragg reflection positions expected from the atomic arrangement and the magnetic spin order, respectively. There are a large number of additional reflections observed below T 1 , comparing to HRNPD data of the FM phase (Figure 7a), e.g., strong satellite peaks below 5 • (2θ) highlighted with a thick arrow. Table 6. Representative operations of the MSSG P2 1 /c.1 (α0γ)00s. The generalized Seitz notation and symmetry operators used in JANA2006 [16] are listed with the label m (no time inversion) and −m (time inversion) involved in the symmetry operation. On the other hand, with 1.5 K HRNPD data collected with a short wavelength, λ = 1.49 Å, it was impossible to refine the atomic structure with the magnetic modulations simultaneously because of a serious overlap of main and satellite reflections. This led to the so-called continuum problem. Nonetheless, this study could deliver the incommensurable modulation vector k(0.523(2), 0, 0.055(1)) for the AMF phase at 1.5 K, which is the first finding for hureaulite-type solid solution compounds. The current model in the incommensurate MSSG P2 1 /c.1 (α0γ)00s ( Figure 13) agrees with observations in neutron powder diffraction more than with a commensurate AFM model with a propagation vector, k(0.5, 0, 0) (Supplementary Figure S4). Furthermore, a faint peak at about 0.4 K in the heat capacity of Mn-hureaulite may be a sign for such a commensurate AFM phase stable at lower temperatures, i.e., below 0.4 K.
For the lost C-centering, five crystallographically independent Mn sites were present in the atomic structure of the AFM phase (Supplementary Table S1e), where magnetic moments at five independent Mn sites were oriented mainly in the (a-c) plane, profoundly in c direction (Table 7). Mn1 remained at its solo unique site below T 1 , but Mn2_1 and Mn2_2 became independent of each other with respect to the site Mn2 in the FM phase, as well as Mn3_1 and Mn3_2 with respect to Mn 3. The maximal magnitude of the magnetic moments at Mn1, Mn2_2, and Mn3_1 was much larger than that at Mn2_1 and Mn3_2. In other words, for different magnitudes of magnetic moments between Mn2_1 and Mn2_2, as well as between Mn3_1 and Mn3_2, the magnetic spin order could be regarded as the main driving force for the low-temperature phase transition at T 1 . Strongly modulated directions and magnitudes of magnetic moments at five independent Mn sites ( Figure 13) are accompanied by a symmetry lowering of atomic arrangements from C2/c to P2 1 /c. The structural part most sensitive to this phase transition is the HB (O-H···O) geometry ( Figure 14). As given in Table 8, the configuration of each HB pair symmetrically equivalent above T 1 strongly differs below T 1 . Furthermore, the onset of proton tunneling is expected for the presence of short HB distance   For the coexistence of strong variation in magnetic modulations at five Mn sublattices with high geometrical distortion in the HB network, the magnetic spin system seems to strongly interplay with structural dissymmetrization in Mn-hureaulite. These observations presume a high complexity of interplaying structural and magnetic order. In this sense, it is desired to use single-crystal diffraction data to be collected best with polarized neutron or synchrotron X-ray sources. This will improve the current new AFM model, but also reveal subtle spin-lattice correlation, eventually driving the phase transition(s) at or below T 1 .
There was no clear indication of a first-order transition for the FM transition at T C = 6.17 K and the incommensurate AFM transition at T 1 =1.86 K. However, when considering the shapes of the C p peaks, the transition at T 1 was more symmetric, but that at T C was more 'step-like', indicating a first-order and a second-order transition, respectively. A clear step was consistently observed in the entropy change for T 1 , while a soft kink appeared at T C (Figure 4c). On the other hand, when the sample was warmed and cooled through a first-order transition, this affected the temperature response C p curve, which no longer had an exponential temperature dependence, and the fitting procedure gave a larger error which could be monitored through the normalized χ 2 and/or the estimated heat capacity error. An increased error obtained for the data around the T 1 transition (Figure 3b) could be an indication of the first-order type. In contrary, a good exponential fit was obtained in the heating and cooling curves C p (T) around T C , which is typical for a second-order transition (i.e., the sample of Mn-hureaulite was heated while continuously following an exponential law through the transition). Hence, the critical point T C may be associated with the second-order transition, accompanying no abrupt jump in the magnetic entropy S m (T) curve (Figure 4c). On the other hand, the existence of a relatively large latent heat at T C , Q L (T C ) = 57 J·mol −1 necessary for the PM → FM transition implies that the transition at T C features a weak first order in nature [23]. In conclusion, the current study cannot give a decisive statement of phase transition types. To accurately determine the order of phase transitions in Mn-hureaulite, further criterial parameters need to be investigated, such as the order of field-dependence of magnetic entropy change and adiabatic temperature change [24].

Conclusions
This study is the first to successfully provide new magnetic structure models of Mnhureaulite. Its ferrimagnetic phase (MSG: C2 /c ) is stable in a small temperature range of 1.86-6.17 K, showing a continuous change of spin-canting orientations. AFM interactions are dominant for both intra and inter pentamers of Mn 2+ (O, HOH) 6 octahedra in a spincanted ferrimagnetism of the title compound. In the incommensurate AFM phase (MSSG: P2 1 /c.1 (α0γ)00s) present below 1.86 K, overall magnetic moments vigorously vary along with distinctive distortion of the HB network. The propagation vector k(0.523(2), 0, 0.055(1)) at 1.5 K could be determined by Rietveld analyses with HRNPD data. A faint peak at about 0.4 K in the heat capacity of Mn-hureaulite adumbrates the presence of the basic commensurate AFM phase possibly with k(0.5, 0, 0) below 0.4 K.
Mn-hureaulite exhibits an interesting spin-canting configuration similar to crossmagnetic sublattices in LiNiPO 4 , for which field-induced phases show a weak ferromagnetism with strong magnetoelectric effects [25]. Considering the results from the present study, it is worth further investigating the hureaulite-type solid solution compounds under various external fields as possible candidates for application-relevant correlated systems.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/sym13091688/s1. We provide all atomic parameters from HRNPD data analyses carried out in this study (Table S1a-e), graphical representations of Rietveld refinements with HRNPD data sets at 4.5 K ( Figure S2a) and 6.1 K (Figure S2b), temperature-dependent lattice metrics and density ( Figure S3), and graphical representation of Rietveld refinements with HRNPD data 1.5 K ( Figure S4a) for a commensurate AFM model ( Figure S4b). These are available online.