Mössbauer Spectroscopy of Triphylite (LiFePO 4 ) at Low Temperatures

: Low temperature magnetic ordering in the LiFePO 4 compound is investigated experimentally using Mössbauer spectroscopy and theoretically via ﬁrst principles calculations. The evaluation of experiment carried out on a powder sample is compatible with an antiferromagnetic order of Fe ion magnetic moments. When an external magnetic ﬁeld is applied, Fe magnetic moments start to deviate slightly from the [010] easy magnetization direction. These ﬁndings are conﬁrmed by means of ﬁrst principles calculations, which also suggest the magnitude of single ion magnetic anisotropy and orbital and spin-dipolar contributions to the magnetic hyperﬁne ﬁeld, which is eventually in a good agreement with the experiment. Diffraction and magnetic measurements complement the study.


Introduction
Triphylite or lithium iron phosphate (LiFePO 4 , LFP) compound is frequently used as a cathode material for Li-ion batteries [1,2]. It offers some advantages, but also shows a few disadvantages when compared to the most frequently used lithium cobalt oxide (LiCoO 2 ) based batteries. Triphylite exhibits a relatively low-symmetry orthorhombic structure [1, 3], which is the origin of many interesting and sometimes puzzling properties. For example, LFP is changing its bond character from ionic to more covalent upon delithiation, thus affecting various battery characteristics [4][5][6]. Another example can be Li diffusion capability, which seems to be hampered by grain boundaries [6,7].
Even if magnetic properties of LiFePO 4 do not affect directly the operation of batteries since the working temperature is well above the ordering temperature, it is worth studying such magnetic phenomena because of their involvement in structural and electronic properties [8]. Magnetic properties of LFP were studied by Santoro and Newnham [9] and Rousse et al. [10] using neutron scattering on polycrystalline/powdered samples. Both works found an antiferromagnetic structure with the Néel temperature (T N ) close to 50 K. Later more precise experiments [11][12][13] performed on single crystals have unveiled a more complicated magnetic structure. Magnetic moments (MMs) are not exactly oriented along the crystallographic axis b (Pnma space group considered), but there are both collinear rotations from this axis and canting of MMs. The overall deviation from the axis b is rather small and amounts to only 1.3 • [13]. In this context, it is worth mentioning that the magnetoelectric effect is expected to appear in LiFePO 4 , as also discussed in Reference [13]. The high magnetic field phase diagram was studied as well [14]. 57 Fe Mössbauer spectroscopy (MS) is capable of detecting very small hyperfine splitting due to the interactions of 57 Fe nuclei with their surroundings. Such interactions comprise specific (static) electric and magnetic coupling of nuclei to electronic charge and spin distributions [15]. This makes MS of Fe sites/nuclei in FeLiPO 4 a valuable tool to complement and check results of other techniques. For instance, magnetization measurements or neutron diffraction studies bring an "integral"/collective response from a sample measured. Moreover, it is certainly meaningful to inspect the sample from a "local" perspective obtained with Mössbauer spectroscopy (by accumulating resonant absorption events originating from individual nuclei). The main subject of the present study is thus low temperature magnetic properties as seen by Mössbauer spectroscopy and the way how MS can help in understanding aforementioned characteristics. Previous MS experimental investigations (see, e.g., References [16][17][18]) indicate that above T N there is no magnetic ordering (only quadrupolar splitting is observed), whereas below T N magnetic order is seen (magnetic hyperfine splitting is combined with the quadrupolar one). The question now arises: What happens if an external magnetic field is applied? Then, is the magnetic anisotropy large enough to prevent magnetic moments becoming aligned along the field? Ab initio calculations may substantially help when answering such questions. The electronic structure, magnetic order and hyperfine interactions in LiFePO 4 are examined computationally within the framework of the density functional theory (DFT) [19], assisting thereby the evaluation and interpretation of experimental MS data. Even if not all queries can be handled perfectly quantitatively, DFT-based approaches give usually valid responses with an unprecedented level of details, which is hard to obtain from experiment.
The paper is organized as follows. Section 2 describes experimental methods used, the studied LiFePO 4 sample, and computational/theoretical approaches. In addition to Mössbauer experiments, the crystal structure and magnetic behavior of the sample are also examined. Section 3 then presents and discusses the obtained results, including the correlation between theory and experiment. The paper is concluded in Section 4, where also an outlook is given. Figure 1 shows the orthorhombic structure of LiFePO 4 . The corresponding space group is Pnma (standard setting, No. 62) [3]. This implies that the lattice parameter order is a > b > c, with the structure mirror planes perpendicular to the axis b (or [010] direction). These mirror planes (with y = 1/4, 3/4) go through the Fe nuclei, which means the C 1h (C s ) local point symmetry of Fe cations. The octahedral coordination of Fe cations (nominally 2+ charge) is also exhibited in Figure 1, where Fe atomic positions are numbered Fe1, Fe2, Fe3, and Fe4. We note that these Fe positions are all crystallographically equivalent (Wyckoff position 4c). Concerning magnetic order (occurring below T N ), MMs at Fe1 and Fe3 sites are parallel, as well as those at Fe2 and Fe4. But MMs at the first group of sites are antiparallel to those at the second group of sites, if a small noncollinearity with the crystallographic axis b is neglected (see Reference [13] for details). Effectively, there are alternating Fe ion planes (y = 1/4 and 3/4) with opposite MM directions collinear with b. At the same time, these planes represent mirrors of the LiFePO 4 crystal. A powder sample obtained from Sigma-Aldrich (Product No. 759546) was studied. The particle size was smaller than 5 µm (BET) and the sample purity was better than 97 % (XRF). During the measurements, the powder particles were fixed using instant glue (solid polymer matrix). In order to check the crystal structure of the sample and find out the corresponding lattice parameters, X-ray diffraction (XRD) measurements were carried out. For this purpose, a Panalytical Empyrean diffractometer equipped with a Mo anode (wavelength λ = 0.70932 Å) working in the parafocusing Brag-Brentano geometry was employed. The diffractometer was furnished with a planar multilayer mirror in the primary beam and axial divergence eliminating Soller slit collimators with an acceptance of 0.02 rad in both the primary and diffracted beams. The diffracted intensity was detected with a 2D hybrid pixel single photon counting detector Galipix equipped with CdTe chip. Measured data were processed using the full powder pattern fitting procedure-Rietveld method. The computer program MStruct [20] was used for the fitting. The magnetic behavior of the studied sample at low temperatures represents also an important aspect in the current investigations. Magnetic properties of the powder LiFePO 4 sample were examined using a Quantum Design superconducting quantum interference device MPMS XL 7T in fields up to 7 T and temperatures down to 2 K.

Experimental and Computational Methods
The 57 Fe Mössbauer spectra of the sample were collected in transmission geometry using a constant-acceleration spectrometer with a 57 Co/Rh source. The calibration of the spectrometer and determination of isomer shifts are given with respect to room temperature Mössbauer spectrum of α-Fe foil. The Mössbauer spectra were acquired in a Janis bath cryostat at temperatures T = 60 K and T = 4.2 K. The in-field spectrum was received at 4.2 K in the external magnetic field B ext = 6 T oriented perpendicularly to the direction of the γ-rays. Measured spectra were evaluated using the current version of MossWinn R fitting programme [21]. The spectra at liquid helium temperature were analyzed by a procedure that makes use of the diagonalization of the full (static) hyperfine interaction Hamiltonian to calculate the relevant energy levels and the transition probabilities between these levels for the 57 Fe isotope, assuming a random powder absorber. The following hyperfine parameters entering the Hamiltonian can, in principle, be extracted from measured Mössbauer 57 Fe spectra: hyperfine magnetic field (B hf ), originating mainly from the Fermi contact interaction, and electric field gradient (EFG) characterized by its principal tensor components V xx , V yy , and V zz (with the increasing magnitude), due to aspherical charge distribution around nucleus. If the EFG asymmetry is defined as η = (V xx − V yy )/V zz , the quadrupole splitting of the 57 Fe nucleus translates into the EFG via the relation QS = eQV zz 1 + η 2 /3 /2 (e and Q being the elementary charge and nuclear quadrupole moment, respectively). The quadrupole moment of the excited 57 Fe nucleus was taken from Reference [22] to be Q = 0.16 barn. If an external magnetic field (B ext ) is applied, then it sums up with B hf in the vectorial way, i.e., one gets an effective magnetic field B eff = |B eff | = |B hf + B ext |, which determines the overall magnetic splitting of the 57 Fe nuclear levels. For powder samples, ∆B eff , the distribution width of the effective magnetic field, is another useful parameter to characterize Mössbauer spectra. When such spectra are fitted, it is yet important to obtain the geometrical relationship of B eff or B hf and the EFG tensor. This relationship is defined using the parameter ϑ, the polar angle between the direction of B eff and the principal axis of the main component V zz of the EFG, and parameter φ, the azimuthal angle between B eff and the V xx 's principal axis. Finally, the isomer shift, IS, characterizes the change of the 57 Fe nucleus' chemical environment with respect to a pure α-iron sample. A negligible spectral component with an intensity of ∼1 % corresponding likely to Fe 3+ ions was subtracted from the measured spectra. This component can be due to a minor content of other phases and/or presence of point defects (in the vicinity of Fe ions, misplaced Fe ions, e.g., at Li sites) [23]. The origin of this spectral component will be studied in the future. The subtraction of the component does not affect the current evaluation of the magnetic structure in powdered LiFePO 4 .
In order to calculate hyperfine parameters that are related to the electronic and magnetic structure of the studied material, a DFT-based approach was employed. In particular, the WIEN2k code [24], which is an implementation of the augmented plane wave plus local orbital [25,26] concept for electronic structure calculations for crystalline solids. In the course of calculations, the spin-polarized exchange-correlation functional-within the generalized gradient approximation [27]-after Perdew, Burke and Ernzerhof (PBE) [28] is used. The spin-orbit coupling (SOC) is also taken into account in specific investigations. The WIEN2k implementation is described in Reference [29], though in the present study the p 1/2 local orbitals were not considered in SOC calculations. The calculations of the EFG and magnetic hyperfine field follows References [30,31] and [32,33], respectively. Figure 2 shows the XRD patterns of the LiFePO 4 sample. According to the powder XRD Rietveld analysis, the sample was assigned to be a single phase of an orthorhombic symmetry with the space group Pnma (No. 62). The observed lattice parameters a = 10.3234(3) Å, b = 6.0045(2) Å and c = 4.6915(1) Å is are in a good agreement with the previous data [1, 3,34,35] summarized also in Reference [2]. In this way, the expected structure of the investigated sample was confirmed. Moreover, we can infer that the Li stoichiometry of the sample is likely good since reduced Li content results in a decrease of the parameters a and b and an increase of the parameter c.  For temperatures T > T N , the molar susceptibility, χ m , see Figure 3, was fitted to the Curie-Weiss law χ m = C/(T + Θ); C and Θ are constants related to the measured system (LiFePO 4 ) and observed magnetic phase transition. The Curie constant C = µ 0 N A µ 2 eff /3k B , with µ 0 being the vacuum permeability, N A the Avogadro constant, k B the Boltzmann constant, and µ eff the effective magnetic moment. The fitted values are Θ = −91(2) K and C = 4.6(1) × 10 −5 m 3 K/mol. From the Curie constant C, the value µ eff = 5.5(1) Bohr magnetons (µ B ) is obtained, in good agreement with Θ = −92(1) K and µ eff = 5.58(1) µ B reported in Reference [36] (see also Reference [17]). The effective magnetic moment is slightly higher than µ eff = 4.90 µ B for Fe 2+ in the high-spin state (S = 2) with the orbital angular momentum quenched (L ≈ 0) by the crystal field. The theoretical value of the high spin state of the free Fe 2+ ion (S = 2, L = 2) is µ eff = 6.71 µ B . A higher value of µ eff observed here thus indicates that the orbital angular momentum was not fully quenched by the crystal field. This observation is consistent with a non-zero orbital Fe 2+ magnetic moment deduced from ab initio calculations discussed below. We mention an increase of the magnetization below temperature ∼ 17 K (Figure 3). Such an increase of the magnetization in LiFePO 4 has also been reported in References [14,17,36]. Rhee et al. [17] related this effect to the influence of the spin-orbit coupling when it becomes comparable with the thermal energy at about 20 K, which results in an 'unquenching' of the orbital magnetic moments of Fe 2+ ions, thus increasing their total magnetic moment. The magnetization curves of the LiFePO 4 sample were measured at temperatures 4.2 K and 60 K and are presented in Figure 4. A linear dependence of the magnetic moment on the magnetic field typical for an AF material was observed at the liquid helium temperature. A larger value of the magnetic moment at 60 K is due to the transition from the AF to the paramagnetic state, which is clearly seen in the temperature dependence of the magnetization in a field of 1 T (Figure 3). We can state that the AF character of the magnetic order in the studied LiFePO 4 sample is confirmed by the magnetic measurements discussed above, and we can take it into account in the subsequent analysis of the Mössbauer spectra.   Table 1 and discussed in detail below. The overall observation is that theory-as discussed in more detail the subsequent section-underestimates the V zz and η parameters and overestimates B hf (when only the Fermi contact interaction contributes to B hf ). The 57 Fe Mössbauer spectrum of LiFePO 4 above the Néel temperature at temperature T = 60 K, see Figure 5, was fitted with a symmetrical, Lorentzian-shaped quadrupole doublet D. The doublet D with reasonably narrow line widths Γ = 0.30(1) mm/s, quadrupole splitting QS = 3.02(2) mm/s and isomer shift IS = 1.33(2) mm/s (in agreement with work [16]), see Table 1, was assigned to ferrous Fe 2+ ions in the high spin state S = 2 [15,37].

Structural and Magnetic Measurements
At temperature T = 4.2 K the 57 Fe Mössbauer spectrum was fitted with one octet O (eight absorption lines [38]) with a narrow distribution of B hf in the static Hamiltonian (see Figure 6). The octet component O with B hf = 12.4(2) T, QS = 3.05(2) mm/s, and IS = 1. 35(2) mm/s were ascribed, similarly to the doublet D at T = 60 K, to Fe 2+ ions with the high spin state S = 2. The orientation of B hf on 57 Fe nuclei is antiparallel to the orientation of the magnetic moment of Fe 2+ ions (like in pure iron), which is parallel to the crystallographic direction [010] [10]. The principal axis of the EFG main component V zz = 16.7(1) × 10 21 V/m 2 , with respect to the zero polar angle ϑ ∼ 0 • , see Table 1, is collinear with the direction [010], which is expectable in some respect (because of the Fe site symmetry, one principal EFG axis should be parallel with this direction). At the Fe 2+ sites, due to a non-zero value of the asymmetry parameter η ∼ 0.8, there is no ordinary local symmetry axis, except the normal to the mirror plane (i.e., C 1h point symmetry).  The in-field Mössbauer spectrum of LiFePO 4 at T = 4.2 K and B ext = 6 T is presented in Figure 7. The application of the external field provokes a broadening of the octet O lines due to the vectorial sum of randomly oriented local B hf and the applied external field B ext perpendicular to the γ-ray direction. The average value of the polar angle ϑ ∼ 13 • , see Table 1, and a higher probability above the mean value of B hf , see inset in Figure 7, indicates the canting of the Fe 2+ moments arranged originally antiferromagnetically. In fact, there are four nonequivalent Fe ion sites with respect to hyperfine interactions when an external magnetic field is applied. Since the single Fe ion magnetic anisotropy is strong, which is explained below, the external field does not cause significant deviation of the Fe ion MM from its original direction along [010], and in the fit we can consider just one Fe site where its MM is slightly deviated from [010], while B eff possesses a broad distribution.

Theoretical Exploration of Magnetism and Hyperfine Interactions
The magnetic structure of LFP was deduced to be antiferromagnetic (AF)-as also confirmed here for the sample studied-with the easy magnetization axis along the [010] direction, based on neutron scattering experiments [9,10]. Other experiments [11][12][13] confirmed this finding, except small deviations from this direction due to anisotropic exchange, which is not taken into account in the following considerations. In addition to the easy magnetization direction, another question arises. Namely, what is the arrangement of the magnetic moments in the LiFePO 4 unit cell? There are three possible MM orders which result in an AF order. These are: Fe1 and Fe2 MMs are parallel (↑↑) being antiparallel (↑↓) to ↓↓ Fe3 and Fe4 (AF1), (Fe1 ↑↑ Fe3) ↑↓ (Fe2 ↓↓ Fe4) (AF2), and (Fe1 ↑↑ Fe4) ↑↓ (Fe2 ↓↓ Fe3) (AF3). Our DFT PBE spin-polarized calculations (without SOC) show that the AF2 order exhibits the lowest energy, followed by AF1 and the highest energy occurs for AF3. In all calculations the lattice parameters from Reference [3], i.e., a = 10.332 Å, b = 6.010 Å and c = 4.692 Å, were employed. Table 2 collects results of these calculations (see the second column). The AF2 order is observed in reality (see, e.g., Reference [10]) and was also confirmed in other computations [8,18].  Table 2 also contains results of hyperfine parameter calculations. One can see that B hf is significantly overestimated (see Table 1 for comparison with experiment), which is likely the effect of neglecting other contributions to the hyperfine field, which will be yet discussed below (see also Reference [17]). The V zz value is underestimated by about 15 %, which is acceptable and a common feature when experimental and computational EFG results are related. The values of asymmetry parameter (η) are also underestimated compared to experiment. AF1 and AF2 show very similar hyperfine parameters, whereas AF3 differs slightly from them. Nevertheless, this indicates a weak sensitivity of the EFG and Fermi contact contribution to B hf to various antiferromagnetic arrangements (Reference [18] reports similar observation).
With regard to the evaluations of more complex experimental situations, we note that different Fe ions have generally different orientations of the EFG principal axes with respect to the lattice translation vectors. This fact is documented in Figure 8 where EFG axes are schematically shown for each of Fe1, Fe2, Fe3, and Fe4 ions together with the directions of translation vectors a, b, and c. Thus, on one hand, Fe1 and Fe2 have the same EFG axis orientations, and, on the other hand, Fe3 and Fe4 have also the same orientations, but distinct from that of Fe1 and Fe2. The V zz axis is always along b, and V yy and V xx axes are always in the a-c plane, though they are tilted from the a and c directions. For B ext = 0, all Fe ions are equivalent with respect to hyperfine interactions. When the external magnetic field is applied along b, EFG principal axis orientations are unimportant since ϑ = 0 • [15]. In this case, Fe1 and Fe3 are pairwise equivalent as well as are Fe2 and Fe4, because of the "magnetic" equivalence within the AF2 order, as discussed above. Considering the most general situation, when B ext is deviated from the easy magnetization axis b, results in the outcome that all Fe sites are nonequivalent regarding the hyperfine interactions.  Figure 1). Coloring of EFG V zz , V yy , and V xx axes/vectors is illustrated for the Fe1 case and kept identical for other cations. The length of vectors corresponds to the size of V zz , V yy , and V xx EFG tensor components.
When the SOC is enabled in the DFT calculations, it is possible to find the easy and hard magnetization directions. Table 3 lists results of three calculations with the magnetization kept collinear with a, b, and c directions while retaining the AF2 order. The total energies per unit cell are given relative to the lowest energy in the second column. One can see that the easy magnetization direction is b ([010]), whereas c ([001]) appears to be the hard direction. The a ([100]) magnetization direction is somewhat softer than the previous one. This is in good agreement with experimental observations [13]. The spin magnetic moment of Fe ions calculated in the muffin-tin sphere is 3.49 µ B , which is reasonable considering that no correction to account for correlated 3d-electrons of Fe was applied (cf., e.g., Reference [39]). The orbital contributions (m l ) to the magnetic moment are also given in Table 3 for all three investigated cases (magnitude is shown only). The largest contribution 0.11 µ B occurs for the easy magnetization direction. When the energy difference (per Fe ion) between the hard and soft direction is compared to the energy of the Fe 2+ ion with a magnetic moment of 4 µ B in the magnetic field B ext , one gets that these energies match at the field B ext 40 T. Then, a significant deviation from the easy magnetization direction may happen only when B ext approaches 40 T, demonstrating a large single ion magnetic anisotropy of Fe ions in LiFePO 4 . This is consistent with the in-field MS results discussed above where a relatively small (average) deviation of the Fe ion MMs from the easy axis direction is observed at the field B ext = 6 T. Moreover, the above considerations about distinct hyperfine interactions at four Fe ion sites confirm the approach regarding how the corresponding MS spectrum was fitted. Table 3. Results of easy/hard magnetization direction calculations with the spin-orbit coupling (SOC) taken into account. The first column indicates the direction with which Fe ion MMs are collinear. ∆E eh is the unit cell total energy relative to the easy direction case. B hf also includes orbital and spin dipolar contributions. When coming to the hyperfine parameters, the EFG-characterized by V zz and η-is almost not affected by switching on the SOC (cf. Table 2). On the other hand, B hf is strongly affected and depends on the magnetization direction. B hf values (see Table 3) now include the orbital magnetic moment and dipolar spin contributions [33] (whereas previous numbers originated from the Fermi contact interaction only). The value B hf = 11.7 T calculated for the magnetization direction collinear with the b axis ([010]) is now very close to the corresponding experimental value 12.4 T (see Table 1).
From these considerations, it becomes obvious that for an adequate description of magnetic hyperfine interactions it is necessary to include their anisotropy. A sufficiently general expression for the nuclear Hamiltonian could be written asĤ whereÎ is the operator of the nuclear spin, Ŝ is the quantum-mechanical and thermodynamic mean value of the Fe ion spin operator, and A is the second order tensor (matrix) describing such type of hyperfine interaction. In principle, each Fe ionic site may have different tensor A in the global coordinate system, which is important in the case when an external magnetic field is applied. A possibly anisotropic (dipolar) contribution from other Fe sites to the hyperfine field should also be considered. Currently, the code [21] used to fit the Mössbauer spectra does not allow us to take into account anisotropic magnetic hyperfine interactions. In the case of an isotropic interaction, the tensor A reduces to a constant A. When the nuclear magnetic moment operator isμ = g N µ NÎ (g N and µ N being, respectively, the g-factor of the 57 Fe nucleus and nuclear magneton), B hf = −A Ŝ /(g N µ N ). Since Ŝ is proportional to the ion magnetic moment, so is B hf . An example of how to handle the anisotropic magnetic hyperfine interaction can be found in Reference [40] where iron phosphide (FeP) having an orthorhombic structure with the same space group as LiFePO 4 was investigated. B hf values from Table 3 could be used to estimate the tensor A (supposed for simplicity to be diagonal in the coordinate system with axes parallel to orthorhomic a, b, and c translation vectors) considering that Fe ions are in the high spin state S = 2 regardless of the magnetic moment orientation. When studying magnetic anisotropy, the single ion anisotropy may play an important role in addition to the exchange interaction anisotropy. The general formula to describe the magnetic anisotropy energy originating from the individual atoms of a crystal reads with K (i) lm being the anisotropy constants of i-th atom/ion and Y lm are spherical harmonics expressed in terms of spherical coordinates (θ polar and ϕ azimutal angle [41]). The summation over atoms goes through all atoms contributing to the anisotropy. Since atomic sites have usually some non trivial symmetry, the effective number of anisotropy constants is reduced. The summation over l is typically limited by 2, 4, or 6.
Using the work [42] and the so-called real spherical harmonics, we can deduce the formula describing magnetocrystalline anisotropy (up to l = 2) per one Fe atom reflecting its C 1h (mirror) point symmetry in LiFePO 4 : + K 20 cos 2 θ + K 22c sin 2 θ cos 2ϕ + K 22s sin 2 θ sin 2ϕ , where K lm [c|s] are anisotropy constants and θ, ϕ are defined with respect to a suitable coordinate system. Namely, θ is the deviation from the axis b, whereas ϕ does not need to be fixed at this moment. Equation (3) is valid for one Fe ion. We have, however, same equations for other ions where anisotropy constants are the same, but θ and ϕ angles are related using symmetry operations (rotations) which mutually transform Fe ion surroundings. In this way, there are just 6 anisotropy constants. They could, in principle, be obtained using ab initio calculations by taking into account the spin-orbit coupling. By considering the external magnetic field and single ion magnetic anisotropy, the direction of MMs at individual Fe sites can be determined. This may serve as a starting point to fit more precisely Mössbauer spectra considering also the anisotropic magnetic hyperfine interaction and EFG tensor principal axis geometrical relationships discussed above.

Conclusions and Outlook
The LiFePO 4 powder sample was first examined experimentally. Its crystal structure at room temperature and magnetic order below the Néel temperature were determined, demonstrating a fair agreement with previous investigations. The Mössbauer study of this LiFePO 4 sample was then performed, showing the good potential of this technique to investigate complex magnetic structures. Measurements done above and below the Néel temperature agree with earlier Mössbauer studies. A doublet and an octet are detected, respectively, in the corresponding spectra. Applying an external magnetic field to the sample at liquid helium temperature results in the deviations of Fe magnetic moments from the crystallographic axis b, being the easy magnetization direction. A small average size of such deviations points to a strong magnetic anisotropy. Moreover, a rather flat distribution of the magnetic effective fields ranging from about 6 to 19 T shows that the external field is added to or subtracted from the hyperfine field, indicating again a large anisotropy. Whereas the distribution should be narrow and peaked at about 6.5 T when the anisotropy is weak-magnetic moments follow the external field, thus destroying the antiferromagnetic order, which was not observed.
By studying theoretically/computationally various features of the magnetic order and hyperfine interactions in LiFePO 4 , we may state a good predictability of experimental results and suggest an advanced method of Mössbauer spectra evaluations, preferably from measurements carried out for single crystal specimens. It is observed that in addition to the Fermi contact term, the orbital magnetic moments of Fe ions contribute significantly to the hyperfine field. The spin-dipolar contribution is important too. The single Fe ion anisotropy is predicted to be high in agreement with experiment. An expression describing the single ion magnetic anisotropy is suggested and its parameters could be obtained in the future by means of SOC DFT calculations with varying directions of the magnetization. This would allow us to determine the direction of magnetic moments at individual Fe sites upon application of the magnetic field, which can be used to fit the Mössbauer spectra more adequately.