## Abstract

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 ZrO_{2} 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 language | English |
---|---|

Pages (from-to) | 953-961 |

Number of pages | 9 |

Journal | Molecular Simulation |

Volume | 30 |

Issue number | 13-15 |

DOIs | |

Publication status | Published - 2004 Nov 15 |

## Keywords

- 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