Thermal conductivity of an imperfect anharmonic crystal

The thermal conductivity of an anharmonic crystal containing randomly distributed substitutional defects due to impurity-phonon scattering is theoretically investigated with the use of the method of double-time thermal Green's functions and the Kubo formalism considering all the terms, i.e., diagonal, nondiagonal, cubic anharmonic, and imperfection terms in the energy-flux operator as propounded by Hardy. The study uses cubic, quartic anharmonic, and defect terms in the Hamiltonian. Mass changes as well as force-constant changes between impurity and host-lattice atoms are taken into account explicitly. It is shown that the total conductivity can be written as a sum of contributions, namely diagonal, nondiagonal, anharmonic, and imperfection contributions. For phonons of small halfwidth, the diagonal contribution has precisely the same form which is obtained from Boltzmann's transport equation for impurity scattering in the relaxation-time approximation. The present study shows that there is a finite contribution of the nondiagonal term, cubic anharmonic term, and the term due to lattice imperfections in the energy-flux operator to the thermal conductivity although the contribution is small compared with that from the diagonal part. We have also discussed the feasibility of numerical evaluation of the various contributions to the thermal conductivity.