Perturbed molecular dynamics for calculating thermal conductivity of Zirconia

M. Yoshiya, A. Harada, M. Takeuchi, K. Matsunaga, H. Matsubara

Research output: Contribution to journalArticlepeer-review

24 Citations (Scopus)


A perturbation method is modified to calculate the thermal conductivity of ionic crystals by molecular dynamics simulation. The energy flux and perturbation tensor are formulated so as to contain only phase space variables and a single convergence parameter. The characteristics of the perturbation method are studied using ZrO2 as a model ionic crystal. The energy flux due to the perturbation is found to exhibit long time-scale oscillations due to the inclusion of long-range Coulombic interactions. In the linear response regime, the calculated thermal conductivity is independent of the external force field parameter, which is the only arbitrary parameter in the method. However, it is found that the external force field parameter plays an important role in minimizing thermal noise in the energy flux. An exponential relationship between thermal conductivity and the maximum external force field parameter is found, by which one can select an external force field parameter suitable for performing thermal conductivity calculations on different ionic materials without needing to carry out numerous preliminary simulations.

Original languageEnglish
Pages (from-to)953-961
Number of pages9
JournalMolecular Simulation
Issue number13-15
Publication statusPublished - 2004 Nov 15


  • Molecular dynamics
  • Non-equilibrium
  • Perturbation
  • Thermal conductivity
  • Zirconia

ASJC Scopus subject areas

  • Chemistry(all)
  • Information Systems
  • Modelling and Simulation
  • Chemical Engineering(all)
  • Materials Science(all)
  • Condensed Matter Physics


Dive into the research topics of 'Perturbed molecular dynamics for calculating thermal conductivity of Zirconia'. Together they form a unique fingerprint.

Cite this