Vacancy diffusion and the hydrodynamics of crystals

The hydrodynamics of crystals with vacancies is developed on the basis of local-equilibrium thermodynamics, where the chemical potential of vacancies plays a key role together with a constraint relating the concentration of vacancies to the density of mass and the strain tensor. The microscopic foundations are established, leading to Green-Kubo and Einstein-Helfand formulas for the transport coefficients, including the vacancy conductivities and the coefficients of vacancy thermodiffusion. As a consequence of having introduced the chemical potential of vacancies, a relationship is obtained between the conductivities and the Fickian diffusion coefficients for the vacancies. The macroscopic equations are linearized around equilibrium to deduce the dispersion relations of the eight hydrodynamic modes. The theoretical predictions are confirmed by numerical simulations of the hard-sphere crystal with vacancies. The study explicitly shows that the eighth hydrodynamic mode of nonperfect monatomic crystals is indeed a mode of vacancy diffusion.