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  -  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 151Eu 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  . When a free nucleus emits a γ-ray (Eγ) during transition from an excited state level with energy (Ee) to the ground state level with energy (Eg) 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 (E0 = Ee − Eg). 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 ,assume a free nucleus of mass M with an excited state level at an energy E0 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 .When a γ-ray of energy Eγ is emitted,the nucleus recoils and has a velocity (V+υ) and a total energy of . By considering conservation of energy before and after the emission of the a γ-ray:
The actual energy of the photon emitted is given by:
The energy Eγ of the emitted photon in the above equation is the excitation energy E0 reduced by the recoil kinetic energy ER of the nucleus and the thermal or Doppler energy (ED = MυV) which is proportional to the initial velocity V from the thermal motion of the nucleus. The recoil nucleus kinetic energy ER expressed by:
Equating the components of momentum before and after the emission process of γ-ray gives:
Figure 1. Recoil in a free nuclei in emission or absorption of a γ-ray.
where (Eγ/c) is the momentum of the photon and c is the speed of light. The recoil momentum of the γ-ray is:
Inserting Equation (5) into Equation (3), the expression for the recoil energy ER of the free nucleus after the emission of the γ-ray is given by:
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 ED 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  which is related to the mean kinetic energy Ēk per translational degree of freedom of a free atom.
where kB 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 -ER and broadened by ED into a Gaussian distribution of width 2ĒD. The distribution for absorption has the same shape but is displaced by +ER. 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 (Г0 ≈ 10−8 eV) of the nuclear energy involved which is determined by the Heisenberg’s uncertainty principle as Г0 = ħ/t, where ħ is Planck’s constant and t is the mean life time of
Figure 2. The statistical distribution in the γ-rays energy for the resonant overlap for successive emission and absorption.
the excited state, e.g. for Eγ = 21.64 keV of 151Eu nucleus, ER ≈ 10−3 eV. This results in a reduction of the Eγ energy, (beyond the Г0 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 2ER ≤ Г0. 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, ER and ĒD become very small and much lower than Г0 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 ER 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 E0 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 and the energy of the gamma quantum will then be too small to be resonantly absorbed by another nucleus. If the recoil energy ER < ħω, then either zero-phonon event or one-phonon event (ħω) of vibrational energy may be transferred to the lattice. It has been shown  through a theoretical consideration that when an average is taken over many emission processes, the energy transferred per event is exactly equal ER. 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 ħω:
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  :
where kγ is the wave number of γ-ray with the wavelength of the gamma ray (λ = 0.57 Å) in case of 151Eu.
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, 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  :
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 kB 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:
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 (ED) as a function of velocity (υ) is equal to  :
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 151Eu. The most popular source for Mössbauer spectroscopy with 151Eu is 151SmF3. Such source has a high recoil free fraction even at room temperature and a single emission line with no appreciable line broadening. The nuclear
Figure 3. (a) Set-up for a Mössbauer spectroscopy experiment in transmission geometry and (b) Typical Mössbauer spectrum in the absence of splitting or shift of the nuclear levels by hyperfine interactions.
Figure 4. The nuclear decay schemes for the Mössbauer isotope 151Eu  .
decay scheme for 151Sm, 21.64 keV level and the value of relevant parameters are shown in Figure 4. The half-life time of 151Sm is 87 years. It decays first to a metastable state of 151Eu by electron capture to the excited state Ie = 7/2 (21.64 keV) or to the ground state Ig = 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Г0) of a resonant emission and absorption of 0.652 mm/s. Other nuclear parameters of the 151Eu 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
Table 1. The nuclear parameters of the Mössbauer isotope 151Eu   .
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−8 eV is less than the hyperfine shifts or splitting of the levels (10−4 - 10−7 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  :
where and are the electron charge densities at the emitting and absorbing nuclei respectively . Z is the atomic number and e is the elementary electronic charge. and are the mean square nuclear radii of the excited and ground states respectively (re ≠ rg). 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, Eu2+ and Eu3+ have electron configurations of (4f7) and (4f6) respectively. The Eu2+ 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 Eu3+ than Eu2+. 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 Eu2+ and between ≈ 0 and ≈ +4 mm/s for Eu3+  . The 151Eu 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 Eu2+ and Eu3+ 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  . 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 . The asymmetry parameter η is generally defined as:
Figure 5. (a) Nuclear energy levels shifts due to electric monopole interaction and (b) Resulting Mössbauer spectrum of the isomer shift.
The electric field gradient can be specified by two parameters η and Vzz. 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 Vzz  :
where γ∞ and R are Sternheimer antishielding and shielding factors respectively. is the contribution of the crystal lattice, i.e. it results from the ionic charges of the atoms which surround the Mössbauer isotope. 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:
with eQ being the nuclear quadrupole moment and I the nuclear spin of the nuclear state. Vzz is the max value of the electrical field gradient along a symmetry axis at the nucleus. EQ dependent on the absolute value of the magnetic quantum number M and it is degenerate in respect of the sign of M. , Ig = 5/2, Ie = 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 Ie = 7/2 is split into four M sublevels (±7/2, ±5/2, ±3/2, ±1/2), while the ground state level with spin Ig = 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  . 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:
where Beff is referred to as the size of the magnetic hyperfine field at the position of the nucleus. (μN = 5.050 × 10−27 JT−1) 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:
M is the magnetic quantum numbers and can take the values . 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 151Eu Mössbauer spectra with 18 allowed possible transitions but experimentally an eight-line pattern is seen
Figure 6. Nuclear energy levels of 151Eu in the presence of quardrupole interaction.
Figure 7. A schematic representation of magnetic splitting of nuclear energy levels of an 151Eu nucleus.
Figure 8. Resulting Mössbauer spectrum.
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 Beff acting on the atomic nucleus originates primarily from the atom’s own electrons can be described in terms of different contributions  . This model has been used to analyze the different contributions to Beff in several Eu intermetallic compounds   .
Bc 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 Eu2+ ion with its own (nonmagnetic) electron shells  . It has the largest contribution to the internal field Beff with accepted value of (−34 T)  . 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.
BF 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  :
where μB is the Bohr magneton and and 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  . 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  .
Bce is the contribution from the conduction electron spin polarization produced by the 4f electron of the magnetic ion itself.
Bthf 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    . However, |Bthf | is very sensitive to external pressure and determine the pressure dependent of |Beff |  .
Bdip is contribution from the dipole interaction with the moment of the electron spin. The total Bdip dipolar field is the sum of the Lorentz field, the demagnetization field, and the contribution from neighboring magnetic dipoles; Bdip is usually very small  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    .
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 Eu2+ state or nonmagnetic trivalent Eu3+ state
Figure 9. Volume vs. atomic number in rare earth elements. Reconstructed from reference  .
in metallic compounds. Figure 9 shows the unit cell volume of rare earth elements  . The magnetic properties in the two valence states are quite different. The ground state configuration of the Eu2+ is 4f7 with ionic ground state 8S7/2, with localized magnetic moment of 7μB. In the Eu3+ the configuration is 4f6 with ionic ground state 7F0. In the stable compounds, the Eu2+ state is generally located well below the Fermi energy EF. The intermediate valence (IV) case is encountered in compounds where the excitation energy Eexc (the distance between the Fermi level EF and localized 4f level) is relatively small Eexc ≤ 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 EuPd2Si2    and EuCu2Si2   . Other example of different rare earths are YbCu2Si2  and CePd3  .
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    . 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−13 sec)  , 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., LIII 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.
Figure 10. Schematic representation of the position of the 4f7 localized magnetic level relative to the Fermi level EF, in stable and intermediate valent systems.
Ionic Model of Valence Fluctuations
The isomer shift (S) deduced from 151Eu Mössbauer spectroscopy plays important role in the field of intermediate valence Eu compounds due to the large difference in its values between the 4f7 (Eu2+) and 4f6 (Eu3+) configurations. In the frame of the ionic the interconfigurational fluctuation (IFC) model  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−13 sec to 10−15 sec   . The two configurations are separated by an excitation energy (Eexc) which is the distance between the Fermi level EF 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   .
The mean valence, ν is given by:
where S(T, V) is the measured value of the isomer shift at 4.2 K, S2(T, V) is the stable isomer shift corresponding to the divalent state. S2 and S3 are the integral values of the isomer shift of Eu2+ and Eu3+ respectively. ∆ν(T, V) is the occupation probability of the Eu3+ 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 S2 and S3. Using the following relation  , the value of ∆ν can be used to evaluate the corresponding values of Eexc at 300 K.
where is the effective temperature, and Tf is being the valence fluctuation temperature, introducing phenomenologically the width of Eexc, i.e. Tf is related to the hybridization energy of the 4f state with the conduction band state, and should tend to zero for Eexc → ±∞  . The multiplicity of the Eu2+ ground state (J = 7/2) is 8, and the ground state Eu3+ (J = 0) is 1. The first and second excited states of Eu3+ (J = 1 and J = 2, respectively) are located at energies of 480 K and 1330 K, respectively  .
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 151Eu 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 Beff at the Eu nucleus. Finally, references of books and articles are given for further and more detailed information.
 Abd-Elmeguid, M.M., Micklitz, H. and Kaindl, G. (1981) Volume Dependence of Magnetic Hyperfine Fields in Eu-Intermetallic Compounds. Physical Review B, 23, 75.
 Kropp, K., Zipf, W., Dormann, E. and Buschow, K.J.H. (1979) Indirect Exchange in Intermetallic Europium Compounds. Journal of Magnetism and Magnetic Materials, 13, 224-230.
 Sampathkumaran, E.V., Gupta, L.C., Vijayaraghavan, R., Pillay, R.G., Gopalakrishnan, K.V. and Devarem H.G. (1981) A New and Unique Eu-Based Mixed Valence System: EuPd2Si2. Journal of Physics C: Solid State Physics, 14, L237.
 Croft, M., Hodges, J.A., Kemly, E., Krishnan, A., Murgai, V. and Gupta, L.C. (1982) Cooperative Configuration Change in EuPd2Si2. Physical Review Letters, 48, 826-829.
 Stober, D., Saner, C., Zinn, W. and Abd-Elmeguid, M.M. (1987) Magnetic Order versus Intermediate Valence in Eu(Pd1-xAux)2Si2 under High Pressure. Journal of Magnetism and Magnetic Materials, 69, 144-148.
 Abd-Elmeguid, M.M., Sauer, C. and Zinn, W. (1985) Pressure-Induced Valence Change of Eu in Eu(Pd0.8Au0.2)2Si2: Collapse of Magnetic Order. Physical Review Letters, 55, 2467-2470.
 Sales, B.C. and Viswanathan, R. (1976) Demagnetization due to Interconfiguration Fluctuations in the RE-Cu2Si2 Compounds. Journal of Low Temperature Phystcs, 23, 449-467.
 Croft, M., Neifeld, R., Segre, C.U., Raaen, S. and Parkss, R.D. (1984) Ce Valence Variation in Intermetallic Alloys: LIII Absorption Spectroscopy Results. Physical Review B, 30, 4164-4169.
 Hossain, Z., Geibel, C., Senthilkumaran, N., Deppe, M., Baenitz, M., Schiller, F. and Molodtsov, S.L. (2004) Antiferromagnetism, Valence Fluctuation, and Heavy-Fermion Behavior in EuCu2(Ge1-xSix)2. Physical Review B, 69, Acticle ID: 014422.
 Fukuda, S., Nakanuma, Y., Sakurai, J., Mitsuda, A., Isikawa, Y., Ishikawa, F., Goto, T. and Yamamoto, T. (2003) Application of Doniach Diagram on Valence Transition in EuCu2(SixGe1-x)2. Journal of the Physical Society of Japan, 72, 3189-3196.
 Nagarajan, R., Sampathkumaran, E.V., Gupta, L.C., Vijayaraghavan, R. and Padalia, B.D. (1981) X-Ray Absorption Spectroscopic Study of a Mixed Valence System, EuPd2Si2. Physics Letters A, 81, 397-398.
 Bauminger, E.R., Froindlich, D., Nowik, I., Ofer, S., Felner, I. and Mayer, I. (1973) Charge Fluctuations in Europium in Metallic EuCu2Si2. Physical Review Letters, 30, 1053-1056.
 Bauminger, E.R., Felner, I., Froindlich, D., Levron, D., Nowik, I., Ofer, S. and Yanovsky, R. (1974) Mössbauer Effect Studies of Interconfiguration Fluctuations in Metallic Rare Earth Compounds. Journal de Physique Colloques, 35, C6.
 Ksenofontov, V., Kandpal, H.C., Ensling, J., Waldeck, M., Johrendt, D., Mewis, A., Gütlich, P. and Felser, C. (2006) Verwey-Type Transition in EuNiP. Europhysics Letters, 74, 672-678.