Carbon diffusion in cementite was explored by using molecular dynamics simulations. The assumption that carbon atoms interacted with each other only indirectly via neighbouring iron atoms was used. An interstitial mechanism of carbon diffusion in cementite was revealed. Carbon diffusion occurred via interstitial sites, which formed four positions per Fe3C unit cell (0.0, 0.0, 0.0), (0.5, 0.5, 0.0), (0.0, 0.0, 0.5) and (0.5, 0.5, 0.5) in units of the lattice parameters: a, b, c. The principal tracer diffusion coefficients and activation parameters of carbon diffusion in cementite were calculated at 1273 to 1373K, and compared with available experimental data. It was argued that carbon diffusion was predominantly a consecutive chain of jumps:

original carbon site → interstitial carbon site → original carbon site → …

The principal tracer diffusion coefficients of carbon atoms for this mechanism were obtained in a form predicted by random-walk theory. The formation energy (about 0.3eV/atom) of defects (carbon atom on an interstitial position and vacant site on original carbon position) as well as the migration energy (about 1.3eV/atom) for an elementary carbon atom jump in cementite were estimated at 1273 to 1373K from the molecular dynamics data.

Molecular Dynamics Simulation and Theoretical Analysis of Carbon Diffusion in Cementite. E.V.Levchenko, A.V.Evteev, I.V.Belova, G.E.Murch: Acta Materialia, 2009, 57[3], 846-53