# The Effect of the Anharmonicity in a Perturbation Estimate of the Self-Diffusion Coefficient

Authors: Takaki Indoh; Masaru Aniya
Journal: Physics Research International
 (2011)

---

## Abstract

A perturbative expression for the self diffusion coefficient is derived by considering that the particles are executing anharmonic vibrations around the sites with local potential minima and jump to neighbor sites after a certain waiting time. The derived expression is a natural extension of the expression obtained by Zwanzig that assumes harmonic oscillations. The relation between the low frequency oscillation and the deviation from the harmonic behavior is also shown.


---

## Body

## 1. Introduction

The dynamics of atoms in solids and liquids have been studied for many years. It is well known that the atomic transport properties in materials are described in terms of the velocity autocorrelation  function (VACF). The VACF carries the information on the many body interactions of the system. Therefore, its calculation is not straightforward. Usually, the VACF is evaluated numerically with the help of computers. However, for the understanding of the physical essence of atomic transport properties, it is very useful if we have at hand an analytical expression for the VACF or related quantity such as the self diffusion coefficient. In the past, some works have been performed along this spirit [1–3]. Among these, the model proposed by Zwanzig is well known [1]. It has been widely used in the literature [4–6] and extended [7]. To describe the dynamics of ions within the materials, the anharmonic effect is also an important factor. For instance, the mobile ion in superionic conductors execute its motion in an anharmonic potential [8]. Therefore, the understanding of the effect of anharmonicity on the atomic transport is directly connected with the understanding of superionic behavior. Visco and Sen discussed the relation between the anharmonic oscillation effect and the normalized VACF [9]. The effect of anharmonicity on the diffusion in glass forming binary Lennard-Jones systems was discussed by Bordat et al. using molecular dynamics simulations [10]. In this study, we modified the assumption used by Zwanzig by incorporating the effect of the anharmonicity and derive a new perturbative expression for the self diffusion coefficient.


## 2. Theory

To describe the effect of anharmonicity, a modified version of the Zwanzig model for the self diffusion coefficient is used. The model is based on the following assumptions.(1)
The particles that form the liquid are executing anharmonic vibrations around the positions where the potential is minima. A fixed frequencyω characterizes adequately the anharmonic oscillations.(2)
When the particle diffuses, it jumps in a very short time to a new position which is almost equivalent to the old position. The waiting time distribution for the diffusion is given byexp(-t/τ), where τ is the lifetime that reflects the interruption of the normal mode [1].(3)
After the jump, the particle loses its memory.In our calculation, the effect of the anharmonicity is incorporated through the magnitude of the anharmonicityf. It causes the deviation from the harmonic oscillation and gives a new expression for the lifetime τ.To discuss the transport property, it is necessary at first to derive the VACF of a particle in an anharmonic potential. The Hamiltonian for a particle with massm that executes anharmonic oscillations is given by(1)H=p22m+12mω02x2-14mfx4.
Here, p is the particle momentum, ω0 is the harmonic oscillator frequency as in the Zwanzig model, f is a positive constant that gives the magnitude of the anharmonicity, and x is the displacement of the particle. The equation of motion obtained from the Hamiltonian is(2)d2xdt2+ω02x-fx3=0.
By solving the equation, we obtain the following expression [9, 11, 12](3)x(t)=a∑n=0∞bnsin⁡(2n+1)ωt,(4)ω=ω0-3fa28ω0+O(f2a4),
where a and bn are constants. The outline of the derivation is given in the appendix. At this point, it is interesting to note that the expression given by (4) may have connection with the low energy excitation and the high anharmonic behavior that some materials exhibit [13–15]. By using (3) and (4), the VACF of a particle in an anharmonic potential is written as
(5)〈pi(t)⋅pi(0)〉〈pi2(0)〉=(∑n=0∞bn2(2n+1)2)-1×∑n=0∞bn2(2n+1)2cos⁡⁡[(2n+1)×(ω0-3fa2/8ω0)+O(f2a4)⋅t].Consider a system that consists ofN atoms with mass m. Under the assumptions that the oscillatory component of the VACF is described by (5) and that the waiting time distribution for the particle jump is accounted for by introducing the factor exp(-t/τ), the following expression for the VACF of a diffusing particle is obtained. This expression is a new expression, which reduces to the expression derived by Zwanzig in the case where f=0. (6)C(t)=1N∑i=1N〈pi(t)⋅pi(0)〉〈pi2(0)〉=cos⁡⁡[(ω0-3fa28ω0)t]exp⁡⁡(-tτ).
In the derivation, only the leading term of the expansion in (3) was retained. The accuracy of the expression can be improved by considering the higher terms.By using (6), we obtain the following expression for the self diffusion coefficient D through the Green-Kubo formula(7)D=kBTm∫0∞C(t)dt=kTmτ1+(ω0-(3fa2)/(8ω0))2τ2.
In particular, if f≪1, we have
(8)D=kBTmτ1+ω02τ2(1+34a2τ2f1+ω02τ2).
The first term of (8) is the usual expression used in many works [1, 4, 5] and arises from the harmonic contributions. The second term of (8) is the correction that arises from the contribution of the anharmonicity.From the Hamiltonian of (1), the maximum value of the potential is mω04/4f. This potential maxima corresponds to the activation energy. On the other hand, τ-1 gives the jumping frequency of a hopping particle. The height of potential and τ-1 is usually related as(9)1τ=Aexp⁡⁡(-mω044fkBT),
where A is a constant. By inserting this expression into (8), we obtain(10)D=kBTm×{A-1exp⁡⁡(mω044fkBT)⋅(1+ω02A-2exp⁡(mω042fkBT))-1+34a2fA-3exp⁡(3mω044fkBT)⋅(1+ω02A-2exp⁡⁡(mω042fkBT))-2}.
At low temperature, (10) becomes(11)D=kBTmAexp⁡⁡(-mω044fkBT)(1ω02+3a2f4ω04).
We can note that this expression generalizes the harmonic expression derived by Zwanzig. The anharmonic contribution is expressed by the second term that contains the magnitude of the anharmonicity f.


## 3. Conclusion

A perturbative expression for the self diffusion coefficient that incorporates the anharmonicity of the atomic vibration has been derived. It is shown that the self diffusion coefficient increases with the magnitude of the anharminicity. In other words, the increase of the anharmonicity results in the decrease of the activation energy. This finding could have implications in diverse topics related to atomic transport, such as in cases where the interatomic potential is soft or in the cases where the material is disordered. In particular, it should be noted that our model catches the essence of ion transport in superionic conductors. For instance, the derived expression shows clearly that the low frequency oscillation is related to the deviations from the harmonic behavior. The low frequency oscillation is one of the phenomena observed in many superionic conductors.



---

*Source: 645264-2011-03-20.xml*