ABSTRACT

In this paper, an overview of the theory of Mössbauer effect is covered, and the main hyperfine interactions parameters which affect the shape of the resultant Mössbauer spectrum are explained and illustrated as well. In principle, Mössbauer effect applies to any and all nuclides, but in practice, certain ideal properties are desirable; that is, the conditions for recoil-free emission and absorption of gamma rays must be optimized. Therefore, briefly discussed in this review, one of the most commonly used for practical and fundamental studies the ¹⁵¹Eu Mössbauer isotope. Also, the intermediate valence phenomena and their theoretical treatments are briefly discussed.

Keywords:

Mössbauer Spectroscopy, Resonant Absorption, Recoil-Free Fraction, ¹⁵¹Eu Mössbauer Isotope, Hyperfine Interactions, Valence Fluctuations

1. Introduction

Mössbauer spectroscopy is a method to measure small shifts and splittings with high accuracy in nuclear energy levels caused by interactions between the electrons and the nucleus. This method is based on the Mössbauer effect [1] - [7] which was discovered in 1958 by the German Physicist Rudolph L. Mössbauer and awarded the Nobel Prize in Physics for his discovery in 1961. The effect involves the emission and the absorption of γ-rays without recoil in solids. The issue of this technique is to present a means of measuring some of the comparatively weak interactions between the nucleus and the surrounding electrons using the Mössbauer spectroscopy. It can be used to provide information in many areas of science such as Physics, Chemistry, Biology and Metallurgy. It can give very precise information about the electronic and magnetic properties of the studied material. In the following sections, the Mössbauer effect is described using one of the most used Mössbauer isotopes among rare earth elements ¹⁵¹Eu and a short explanation is given on the hyperfine interactions which determine the structure in a Mössbauer spectrum.

1.1. Recoilless Emission and Absorption of γ-Rays

Mössbauer showed that nuclear radiation can be emitted and absorbed by another nucleus of the same kind as the emitter without recoil if the nucleus placed in a solid [1] . When a free nucleus emits a γ-ray (E_γ) during transition from an excited state level with energy (E_e) to the ground state level with energy (E_g) it will recoil. This recoil energy and velocity are due to conservation of energy and momentum, resulting in the emitted γ-ray being of lower energy than the nuclear transition energy (E₀ = E_e − E_g). The same is observed where the absorbing nucleus recoils, meaning the energy of the resonantly absorbed photon is less than the transition energy. To understand this effect [7],assume a free nucleus of mass M with an excited state level at an energy E₀ and moving with a velocity V along the direction of the emission of E_γ energy (Figure 1).

The energy above the ground state at rest is $E_0+\left(1/2\right)M{V}^2$.When a γ-ray of energy E_γ is emitted,the nucleus recoils and has a velocity (V+υ) and a total energy of $E_0+\left(1/2\right)M{\left(V+\upsilon \right)}^2$. By considering conservation of energy before and after the emission of the a γ-ray:

$E_0+\frac{1}{2}M{V}^2=E_{\gamma }+\frac{1}{2}M{\left(V+\upsilon \right)}^2$ (1)

The actual energy of the photon emitted is given by:

$\begin{array}{c}{E}_{\gamma }=E_0-\frac{1}{2}M{\upsilon }^2-M\upsilon V\\ =E_0-E_R-E_D\end{array}$ (2)

The energy E_γ of the emitted photon in the above equation is the excitation energy E₀ reduced by the recoil kinetic energy E_R of the nucleus and the thermal or Doppler energy (E_D = MυV) which is proportional to the initial velocity V from the thermal motion of the nucleus. The recoil nucleus kinetic energy E_R expressed by:

$E_R=\frac{1}{2}M{\upsilon }^2$ (3)

Equating the components of momentum before and after the emission process of γ-ray gives:

$MV=M\left(V+\upsilon \right)+\frac{E_{\gamma }}{c}$ (4)

where (E_γ/c) is the momentum of the photon and c is the speed of light. The recoil momentum of the γ-ray is:

$M\upsilon =-\frac{E_{\gamma }}{c}$ (5)

Inserting Equation (5) into Equation (3), the expression for the recoil energy E_R of the free nucleus after the emission of the γ-ray is given by:

$E_R=\frac{E_{\gamma }^2}{2M{c}^2}$ (6)

Since the atoms of the emitter as well as the absorber have thermal energy, they move in certain velocities, and for each velocity the emitted wave’s frequency is changed according to the Doppler shift. The Doppler energy E_D is dependent on the thermal motion of the nucleus, and will therefore have a distribution of values which is temperature dependent. A mean value, Ē_D, can be expressed [7] which is related to the mean kinetic energy Ē_k per translational degree of freedom of a free atom.

${\bar{E}}_D=E_{\gamma }\sqrt{\frac{2{\bar{E}}_k}{M{c}^2}};{\bar{E}}_k\cong \frac{1}{2}{k}_{B}T$ (7)

where k_B is Boltzmann’s constant and T is the absolute temperature. As a result, the statistical distribution in energy of the emitted γ-rays is displaced from the true excited state energy by -E_R and broadened by E_D into a Gaussian distribution of width 2Ē_D. The distribution for absorption has the same shape but is displaced by +E_R. Figure 2 illustrated schematically the statistical distributions for the emitting and absorbing atoms. These energy distributions must strongly overlap in order for resonant absorption to occur.

The recoil energy is much larger than the natural linewidth (Г₀ ≈ 10⁻⁸ eV) of the nuclear energy involved which is determined by the Heisenberg’s uncertainty principle as Г₀ = ħ/t, where ħ is Planck’s constant and t is the mean life time of

the excited state, e.g. for E_γ = 21.64 keV of ¹⁵¹Eu nucleus, E_R ≈ 10⁻³ eV. This results in a reduction of the E_γ energy, (beyond the Г₀ limit) therefore, nuclear resonant absorption of γ-ray emitted by the free nuclei can not occur.

Thus, the resonance emission and absorption can take place under the condition that 2E_R ≤ Г₀. This resonance condition is fulfilled when the recoil momentum or energy upon emission or absorption of a γ-ray is taken up by the solid as a whole rather than by a free nucleus. As the equations for recoil energy and Doppler broadening contain the reciprocal mass M, E_R and Ē_D become very small and much lower than Г₀ and thereby fulfill the resonance condition.

If the nucleus is bound in a solid, in which the movements of atoms are superpositions of quantized vibrational modes, then the recoil energy E_R must be almost as large as the energy of the vibration quanta phonon, in order to excite a phonon. In the Einstein model of lattice vibrations, there is only one phonon frequency ω. The lattice vibrations are quantized, and the recoil energy can only be transferred to the lattice if it takes place in integral multiples ħω. During the emission of a gamma quantum, part of the energy E₀ may be lost to the lattice by exciting a vibrational state. In this case the energy of the gamma quantum is reduced by the amount $\left(0,\pm \hslash \omega ,\pm 2\hslash \omega ,\cdots \right)$ and the energy of the gamma quantum will then be too small to be resonantly absorbed by another nucleus. If the recoil energy E_R < ħω, then either zero-phonon event or one-phonon event (ħω) of vibrational energy may be transferred to the lattice. It has been shown [8] through a theoretical consideration that when an average is taken over many emission processes, the energy transferred per event is exactly equal E_R. Using this fact we may write an expression which leads directly to the fraction, f, of emission events which takes place without lattice excitation (zero-phonon transition), and a fraction (1 − f) transfers one phonon energy ħω:

$f=\left(1-f\right)\hslash \omega \text{or}f=1-\frac{E_R}{\hslash \omega }$ (8)

Only these events (i.e. recoilless fraction f) give rise to resonance emission or absorption, i.e. to the Mössbauer effect. f is called Mössbauer-Lamb factor.

1.2. The Lamb-Mössbauer Factor

The probability of recoil-free fraction of γ-rays emitted or absorbed without loss of energy to the lattice vibrations is related to the Mössbauer-Lamb factor f which can be written as [3] :

$f=\exp \left(-k_{\gamma }^2\langle x^2\rangle \right);k_{\gamma }=\frac{E_{\gamma }}{\hslash c}=\frac{2\text{π}}{\lambda }$ (9)

where k_γ is the wave number of γ-ray with the wavelength of the gamma ray (λ = 0.57 Å) in case of ¹⁵¹Eu.

$\langle x^2\rangle$ is the mean square amplitude for thermal vibrations of the emitting or absorbing nucleus in the direction of γ-ray. The Mössbauer Lamb factor f value when it increases, leads to an increase the Mössbauer effect at relatively low γ-ray energy and at low temperature, $\langle x^2\rangle$ decreases with lowering temperature. Therefore, most Mössbauer spectroscopy experiments are performed at low temperature such as liquid nitrogen or liquid helium. For the Debye model of lattice vibrations Equation (9) has the form [3] :

$f=\exp \left(-\frac{3E_R}{2k_B{\Theta }_D}\left(1+4{\left(\frac{T}{{\Theta }_D}\right)}^2{\displaystyle \underset{0}{\overset{\frac{{\Theta }_D}{T}}{\int }}\frac{x}{{\text{e}}^x-1}\text{d}x}\right)\right)$ (10)

where Θ_D is the Debye temperature of the crystal in which the nucleus is bound, T is the absolute temperature of the crystal at which the experiment is performed and k_B is the Boltzmann constant. Since the fraction f is strongly temperature dependent and relative to the Debye temperature, decreases as the temperature increases and conversely, an approximations can be made for very low and high temperature conditions:

$f=\exp \left(\frac{-3E_R}{2k_B{\Theta }_D}\right)\text{for}T\ll {\Theta }_D;$ (11)

$f=\exp \left(\frac{-6E_{R}T}{k_B{\Theta }_D^2}\right)\text{for}T>{\Theta }_D;$ (12)

1.3. Principle of Mössbauer Spectroscopy

In the transmission geometry of Mössbauer experiment, a sample is exposed to a γ-ray source and a detector measures the intensity of the beam transmitted through the sample. This technique is achieved by moving a source relative to an absorber at a variable velocity. In this way, the energy of the emitted γ-rays is shifted due to the Doppler effect in order to have a resonance between a source and an absorber, the energy increases if a source moves toward an absorber and decreases when it moves away. Shifting and splitting of the nuclear energy levels between a source and an absorber can be measured at velocities in the order of mm/s. The resulting Mössbauer spectra is a plot of γ-ray transmission intensity as a function of the source velocity which can be converted to energy shift. A basic set-up for a Mössbauer experiment in transmission geometry and a typical result are shown Figure 3(a) and Figure 3(b) respectively. The Doppler energy shift (E_D) as a function of velocity (υ) is equal to [3] :

$E_D=E_{\gamma }\frac{\upsilon }{c}$ (13)

1.4. Gamma Ray Source

Gamma ray is the most energetic radiation in the electromagnetic spectrum; it has the lowest wavelength and the highest frequency in the spectrum. γ-rays are sent out spontaneously by some radioactive nuclei and having energies in a range 10 - 100 keV. One of the most used Mössbauer isotopes among rare earth elements is ¹⁵¹Eu. The most popular source for Mössbauer spectroscopy with ¹⁵¹Eu is ¹⁵¹SmF₃. Such source has a high recoil free fraction even at room temperature and a single emission line with no appreciable line broadening. The nuclear

decay scheme for ¹⁵¹Sm, 21.64 keV level and the value of relevant parameters are shown in Figure 4. The half-life time of ¹⁵¹Sm is 87 years. It decays first to a metastable state of ¹⁵¹Eu by electron capture to the excited state I_e = 7/2 (21.64 keV) or to the ground state I_g = 5/2 by the γ-ray emission. The 21.64 keV state decays in turn to the ground state by γ-ray emission or internal conversion. The half-life of the 21.64 keV excited state is 9.7(3) ns, giving rise to energy uncertainty. A Mössbauer γ-ray then has a natural linewidth (2Г₀) of a resonant emission and absorption of 0.652 mm/s. Other nuclear parameters of the ¹⁵¹Eu isotope are shown in the Table 1.

2. Hyperfine Interactions

The primary application of Mössbauer is to use the nucleus as a local probe when it is embedded in a solid state lattice. The interaction between a nucleus and its surrounding environment is known as a hyperfine interaction. These interactions are very small compared to the energy levels of the nucleus itself but the extreme energy resolution of the Mössbauer effect enables these interactions to be observed. The very high energy resolution of the Mössbauer spectroscopy permits small energy shifts (~ neV) to be measured. The physical reason for such an energy shift is the hyperfine interactions of a nuclear property with that of an

electronic or atomic property of the atom. These interactions can cause degenerate nuclear excited and ground states to split or they can lead to small shifts in energy levels. Often, the linewidth of the γ-ray 10⁻⁸ eV is less than the hyperfine shifts or splitting of the levels (10⁻⁴ - 10⁻⁷ eV), which then will be reflected in the Mössbauer spectrum, partly as a splitting of the lines and partly as shifts. Three main interactions give rise to the structure of a Mössbauer spectrum. These interactions give rise to the isomer shift, electric quadrupole splitting and the magnetic hyperfine splitting. In the following sections these interactions and associated hyperfine parameters and their influence on the nuclear energy levels are discussed.

2.1. The Electric Monopole Interaction (Isomer Shift)

The isomer shift arises due to the non-zero volume of the nucleus and it is directly proportional to s-electron density at the same nucleus. This leads to an electric monopole or Coulombic interaction, which is the interaction between electrons of the atom and the nucleus. The difference in s-electron density of the emitting and absorbing nuclei and the difference between the radius sizes of the nucleus in the ground and excited states affect this interaction, and lead to only a shift of the nuclear energy levels Figure 5(a). In the resulting Mössbauer spectrum, such as Figure 5(b), the isomer shift is the displacement of the absorption line away from zero velocity and is observed where the Doppler velocity applied to the source correctly modifies the source’s gamma-rays energy to account for the difference between the source and the absorber. This relation between charge electron density and nuclear radius is so called isomer shift (S) and is given by the equation [7] :

$S=\frac{\text{2π}}{3}Z{e}^2\left({\left|{\Psi }_a\left(0\right)\right|}^2-{\left|{\Psi }_s\left(0\right)\right|}^2\right)\left(\langle r_e^2\rangle -\langle r_g^2\rangle \right)$ (14)

where ${\left|{\Psi }_s\left(0\right)\right|}^2$ and ${\left|{\Psi }_a\left(0\right)\right|}^2$ are the electron charge densities at the emitting and absorbing nuclei respectively $\left({\left|{\Psi }_s\left(0\right)\right|}^2\ne {\left|{\Psi }_a\left(0\right)\right|}^2\right)$. Z is the atomic number and e is the elementary electronic charge. $\langle r_e^2\rangle$ and $\langle r_e^2\rangle$ are the mean square nuclear radii of the excited and ground states respectively (r_e ≠ r_g). The isomer shift measurements can provide useful information on the electronic structure surrounding the nucleus and is good for probing the valence state of Mössbauer atom. As the wave function of the s-electrons penetrate into outer shells changes in these shells will directly alter the s-electron charge density at the nucleus. For example, Eu²⁺ and Eu³⁺ have electron configurations of (4f⁷) and (4f⁶) respectively. The Eu²⁺ have less s-electron density at the nucleus due to the greater screening effect of p, d, and f-electrons. This produces a positive isomer shift greater in Eu³⁺ than Eu²⁺. The isomer shift of Eu nuclei in compounds containing stable divalent and trivalent Eu ions falls in two well known separated regions and varies in the range between ≈−8 and ≈−14 mm/s for Eu²⁺ and between ≈ 0 and ≈ +4 mm/s for Eu³⁺ [12] . The ¹⁵¹Eu Mössbauer isomer shift can be used to investigate the valence state in Eu-based compounds, due to this larger difference (∆ ≈ 14 mm/s) between Eu²⁺ and Eu³⁺ compared to the linewidth of the resonance line (3 mm/s).

2.2. The Electric Quadrupole Interaction (Quadrupole Splitting)

If the nucleus has a non-spherical charge distribution (nuclear spin quantum number I > 1/2), then it has electric quadrupole moment [7] . Therefore, the quadrupole moment can interact with an electric field gradient (EFG) at the nucleus due to the surrounding electrons, producing a splitting in the nuclear energy levels. The electric field gradient could have two components: one due to the asymmetry of the electron charge and another to the asymmetry of the lattice charge distribution. Such interaction depends on the squared nuclear spin quantum number and lift only partially the degeneracy of the nuclear levels. The coordinate system of the electric field gradient tensor is chosen such that $\left|V_{xx}\right|\le \left|V_{yy}\right|\le \left|V_{zz}\right|$. The asymmetry parameter η is generally defined as:

$\eta =\frac{V_{zz}-V_{yy}}{V_{zz}};0\le \left|\eta \right|\le 1$ (15)

The electric field gradient can be specified by two parameters η and V_zz. The parameter η describes deviations of the electrical field gradient from rotational symmetry around the principal component (z-axis). The z-axis is often chosen to be the symmetry axis in describing a nucleus which is subjected to this hyperfine interaction, and is the axis along which the EFG is a maximum. There are two contributions to V_zz [13] :

$V_{zz}=\left(1-{\gamma }_{\infty }\right)\cdot V_{zz}^{lat}+\left(1-R\right)\cdot V_{zz}^{val}$ (16)

where γ_∞ and R are Sternheimer antishielding and shielding factors respectively. $V_{zz}^{lat}$ is the contribution of the crystal lattice, i.e. it results from the ionic charges of the atoms which surround the Mössbauer isotope. $V_{zz}^{val}$ describes the contribution of the localized electron of Mössbauer atoms, due to incompletely filled electron shells (valence electrons). In 4f systems this contribution is mostly dominated by a non-spherical charge distribution of the 4f electrons, however in some cases, other valence electrons (e.g. 6p or 5d) have a significant influence. The expression for the electric quardrupole interaction gives rise to the following energy levels scheme:

$E_Q=\frac{eQ{V}_{zz}}{4I\left(2I-1\right)}\left(3M^2-I\left(I+1\right)\right){\left(1+\frac{{\eta }^2}{3}\right)}^{\frac{1}{2}}$ (17)

with eQ being the nuclear quadrupole moment and I the nuclear spin of the nuclear state. V_zz is the max value of the electrical field gradient along a symmetry axis at the nucleus. E_Q dependent on the absolute value of the magnetic quantum number M and it is degenerate in respect of the sign of M. $M=-I,-I+1,\cdots ,I$, I_g = 5/2, I_e = 7/2. It means these two states configurations have the same energy. The magnetic quantum numbers determine the number of energy levels. Therefore, they determine the number of Mössbauer lines. The excited state level spin I_e = 7/2 is split into four M sublevels (±7/2, ±5/2, ±3/2, ±1/2), while the ground state level with spin I_g = 5/2 is split into three M sublevels (±5/2, ±3/2, ±1/2), as shown in Figure 6. The resultant spectrum is from transitions between ground and excited state sublevels subject to the selection rule ∆M = 0, ±1.

2.3. Magnetic Hyperfine Field (Nuclear Zeeman Effect)

The magnetic hyperfine splitting is due to the interaction of magnetic moments μ of the ground and excited states of the nucleus with a magnetic field with the condition that the spin quantum number I > 0 [7] . Such interactions lift completely the degeneracy of the nuclear levels which correspond to different orientations of the spin I. The Hamiltonian for magnetic hyperfine dipole interaction is given as:

$H_m=-\mu \cdot B_{eff}=-g{\mu }_{N}I\cdot B_{eff}$ (18)

where B_eff is referred to as the size of the magnetic hyperfine field at the position of the nucleus. (μ_N = 5.050 × 10⁻²⁷ JT⁻¹) is the nuclear magneton, and g is referred

to as the g-factor (g = μ/Iμ_N) of the nuclear state. The spin I and g are usually different for the excited and ground states. Due to this interaction, the magnetic field splits the nuclear energy levels of spin quantum number I into 2I + 1 equally spaced non-degenerate levels with an energy given by:

$E_m=-\frac{\mu B}{I}M=-g{\mu }_N{B}_{eff}M$ (19)

M is the magnetic quantum numbers and can take the values $M=-I,-I+1,\cdots ,I-1,I$. Transition can take place between sublevels subject to the selection rule that ∆M should change by 0, ±1 in the absorption process. Figure 7 shows the splitting of the excited state (I = 7/2) and ground state (I = 5/2) nuclear levels which gives rise to very complex ¹⁵¹Eu Mössbauer spectra with 18 allowed possible transitions but experimentally an eight-line pattern is seen

in the case of good counting statistics or an intense spectrum. Figure 8 shows the corresponding Mössbauer lines.

The total effective magnetic hyperfine internal field B_eff acting on the atomic nucleus originates primarily from the atom’s own electrons can be described in terms of different contributions [14] . This model has been used to analyze the different contributions to B_eff in several Eu intermetallic compounds [15] [16] .

$B_{eff}=B_c+B_{ce}+B_{thf}+B_{dip}$ (20)

B_c is the Fermi contact field which is caused by the core polarization of the s-electron density at the nucleus. It is originating from the interaction of the S = 7/2 spin of the Eu²⁺ ion with its own (nonmagnetic) electron shells [18] . It has the largest contribution to the internal field B_eff with accepted value of (−34 T) [19] . In the nonrelativistic approximation, this effective field is regarded as being due to the difference of spin up and spin down densities m(0) at the nucleus.

$B_F=\frac{\text{8π}}{3}{\mu }_{B}m\left(0\right)$ (21)

B_F is produced by the polarization of electrons whose wavefunctions overlap the nucleus, i.e. s-electrons. This polarization is due to unpaired electrons in the f orbital and gives an imbalance in spin density at the nucleus from the difference in interaction between the unpaired electron with s-electrons of parallel or antiparallel spin to its own. This can be expressed formally as [20] :

$B_F=-\frac{\text{8π}}{3}{\mu }_B{\displaystyle \sum \left\{{\left|{\Psi }_{\uparrow }^s\left(0\right)\right|}^2-{\left|{\Psi }_{\downarrow }^s\left(0\right)\right|}^2\right\}}$ (22)

where μ_B is the Bohr magneton and ${\left|{\Psi }_{\uparrow }^s\left(0\right)\right|}^2$ and ${\left|{\Psi }_{\downarrow }^s\left(0\right)\right|}^2$ are the spin up and spin down densities for s-electrons at the nucleus respectively. The summation is taken over all s-electrons. The large negative part of this field is owing to core polarization. This results from the different exchange interactions of the 1 s, 2 s, 3 s and 4 s electrons of parallel and antiparallel spin with the f shall [21] . The exchange interaction between the spin up polarized f shell and the spin up s-electron is attractive, while that between the f shell and a spin down s-electron is repulsive. As a result the radial parts of the two s-electron wave functions will be different, one being pushed toward the nucleus, the other pulled outward. Therefore, the spin densities at the nucleus no longer cancel and a Fermi contact interaction is observed [5] .

B_ce is the contribution from the conduction electron spin polarization produced by the 4f electron of the magnetic ion itself.

B_thf is the transferred hyperfine field produced at an Eu nucleus by the spin polarization of the conduction electrons by the 4f magnetic moment of the neighboring atoms. It has been found the transferred hyperfine fields are relatively small in Eu-intermetallic compounds [15] [16] [17] . However, |B_thf | is very sensitive to external pressure and determine the pressure dependent of |B_eff | [14] .

B_dip is contribution from the dipole interaction with the moment of the electron spin. The total B_dip dipolar field is the sum of the Lorentz field, the demagnetization field, and the contribution from neighboring magnetic dipoles; B_dip is usually very small [22] compared to the other contributions.

3. Valence Fluctuations in Eu Intermetallic Compounds

The phenomena of valence fluctuations has been observed in certain 4f systems, such as Ce, Pr, Sm, Eu and Yb compounds. This is due to the fact that the respective rare earth ions can in principle exist in two different valence states with different numbers of 4f electrons (see Figure 9). Temperature and pressure variations, or alloying with proper elements, may lead to a quasidegeneration of the two valence states, giving rise to a nonintegral average occupation number of the 4f shell (i.e. intermediate valence state). The transition from stable (normal) into the intermediate valence state is accompanied by anomalous changes of the physical properties of the 4f-system. An overview on the intermediate valence phenomena and their theoretical treatments have been given by [23] [24] [25] .

In this investigation, we focus on the description of the phenomena of magnetically ordered Eu intermetallic compounds systems. Europium ions are usually found to be either divalent Eu²⁺ state or nonmagnetic trivalent Eu³⁺ state

in metallic compounds. Figure 9 shows the unit cell volume of rare earth elements [26] . The magnetic properties in the two valence states are quite different. The ground state configuration of the Eu²⁺ is 4f⁷ with ionic ground state ⁸S_7/2, with localized magnetic moment of 7μ_B. In the Eu³⁺ the configuration is 4f⁶ with ionic ground state ⁷F₀. In the stable compounds, the Eu²⁺ state is generally located well below the Fermi energy E_F. The intermediate valence (IV) case is encountered in compounds where the excitation energy E_exc (the distance between the Fermi level E_F and localized 4f level) is relatively small E_exc ≤ 300 K, Figure 10. In such compounds the Eu ions exhibit an IV state and undergo a valence change if one changes either temperature or pressure as it has been observed, e.g. in EuPd₂Si₂ [27] [28] [29] and EuCu₂Si₂ [30] [31] . Other example of different rare earths are YbCu₂Si₂ [32] and CePd₃ [33] .

There are varieties of techniques to measure the mean valence in an intermediate valence system. Measurements of macroscopic properties such as magnetic susceptibility and electrical transport are indirectly affected by the existence of the valence fluctuations of the system [34] [35] [36] . On the other hand, measurements of microscopic properties by, e.g. Mössbauer effect, photoelectric emission, nuclear resonant scattering and x-ray absorption techniques provide direct information about the valence fluctuations. The difference between these techniques is set by the timescale to which they are sensitive. Techniques with a timescale that is short with respect to that of the valence fluctuation (≈ 10⁻¹³ sec) [37] , are sensitive to the momentary values. Consequently, the measured spectrum will be a superposition of the spectra corresponding to both discrete valence states. This is the case in, e.g., L_III x-ray absorption edge studies. The techniques with a timescale that is much longer than of the valence fluctuations, such as Mössbauer spectroscopy the measured spectrum gives an average value for the valence.

Ionic Model of Valence Fluctuations

The isomer shift (S) deduced from ¹⁵¹Eu Mössbauer spectroscopy plays important role in the field of intermediate valence Eu compounds due to the large difference in its values between the 4f⁷ (Eu²⁺) and 4f⁶ (Eu³⁺) configurations. In the frame of the ionic the interconfigurational fluctuation (IFC) model [34] it is assumed that the valence of the Eu ions fluctuates rapidly between, the localized 4f level and the conduction band, i.e. fluctuates between the integral values 2 < ν < 3 on a time scale in the order of 10⁻¹³ sec to 10⁻¹⁵ sec [38] [39] . The two configurations are separated by an excitation energy (E_exc) which is the distance between the Fermi level E_F and localized 4f level. In this frame, the measured isomer shift S(T, V) of the Eu ion is converted to the Eu mean valence ν through the formula [31] [39] .

$S\left(T,V\right)=S_2\left(T,V\right)+S\left(S_3-S_2\right)\cdot \Delta \nu \left(T,V\right)$ (23)

The mean valence, ν is given by:

$\nu =2+\Delta \nu \left(T,V\right)$ (24)

where S(T, V) is the measured value of the isomer shift at 4.2 K, S₂(T, V) is the stable isomer shift corresponding to the divalent state. S₂ and S₃ are the integral values of the isomer shift of Eu²⁺ and Eu³⁺ respectively. ∆ν(T, V) is the occupation probability of the Eu³⁺ state, given the mean valence ν and be calculated from the measured values of S at different concentrations and temperatures with suitable choice of the values for S₂ and S₃. Using the following relation [40] , the value of ∆ν can be used to evaluate the corresponding values of E_exc at 300 K.

$\Delta \nu \left(T,V\right)=\frac{8\exp \left[-E_{exc}\left(T,V\right)/K_{B}T{}^{\ast }\right]}{1+3\exp \left[-480/T{}^{\ast }\right]+5\exp \left[-1330/T{}^{\ast }\right]+\cdots }$ (25)

where $T^{\ast }=\sqrt{T^2+T_f^2}$ is the effective temperature, and T_f is being the valence fluctuation temperature, introducing phenomenologically the width of E_exc, i.e. T_f is related to the hybridization energy of the 4f state with the conduction band state, and should tend to zero for E_exc → ±∞ [30] . The multiplicity of the Eu²⁺ ground state (J = 7/2) is 8, and the ground state Eu³⁺ (J = 0) is 1. The first and second excited states of Eu³⁺ (J = 1 and J = 2, respectively) are located at energies of 480 K and 1330 K, respectively [41] .

4. Conclusion

The goal of this review is to create an accessible summary and familiarize the reader who is not an active practitioner of Mössbauer spectroscopy with the basic concepts of the method. First, the theory behind the Mössbauer effect is explained and the key to the success of the technique is the discovery of recoilless gamma ray emission and absorption. It is intended that this information should be sufficient for the reader to understand the application of Mössbauer spectroscopy to analysis the studied system. This paper also presents some ideas based on the nuclear physics hyperfine interactions parameters that can by derived from the Mossbauer effect. These interactions are very small compared to the energy levels of the nucleus itself but the extreme energy resolution of the Mössbauer effect enables these interactions to be observed. The very high energy resolution of the Mössbauer spectroscopy permits small energy shifts (~ neV) to be measured. Therefore, for further studies using Mössbauer spectroscopy technique with external parameters such as temperature and pressure can very precise information about valence state of Mössbauer atom, structural electronic and magnetic properties of the studied material. One could gain information at a microscopic level about valence fluctuations in Eu systems via the ¹⁵¹Eu Mössbauer isomer shift S. Also, could gain a microscope insight into the change of the magnetic moment of the Eu ions via the effective magnetic hyperfine field B_eff at the Eu nucleus. Finally, references of books and articles are given for further and more detailed information.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

Cite this paper

Ahmida, M.A. (2019) Theoretical Review of Mössbauer Effect, Hyperfine Interactions Parameters and the Valence Fluctuations in Eu Systems. Journal of Applied Mathematics and Physics, 7, 254-270. https://doi.org/10.4236/jamp.2019.71021

References