Thermoelectric Na0.8_{0.8}CoO2_2 suppresses its lattice thermal conductivity by a factor of about six through anharmonic rattling vibrations of sodium ions inside ordered multi-vacancy clusters, as established by inelastic X-ray/neutron scattering and density-functional phonon calculations (Voneshen et al., 2013). What remains open is the quantitative microscopic theory behind that observation: (i) a first-principles calculation of phonon linewidths/lifetimes that decomposes the measured widths into anharmonic phonon-phonon and electron-phonon contributions and predicts κL(T)\kappa^L(T); (ii) identification of the mechanism behind the ~50% linewidth broadening and slight dispersion softening that appear when long-range vacancy order is lost at elevated temperature; and (iii) a systematic first-principles and experimental map of how vacancy-cluster arrangement and sodium concentration control rattler energies, anharmonicity, and κL\kappa^L suppression, yielding predictive design rules.

Background

A thermoelectric material converts heat into electricity, and its efficiency is gauged by the dimensionless figure of merit ZT=S2σT/κZT = S^2\sigma T/\kappa, where SS is the Seebeck coefficient, σ\sigma the electrical conductivity, and κ\kappa the thermal conductivity (the sum of an electronic part and a lattice part κL\kappa^L carried by phonons, the quantized lattice vibrations). Good thermoelectrics therefore need high σ\sigma but low κL\kappa^L, a design conflict that motivates the phonon-glass electron-crystal concept: crystals in which charge carriers move freely while loosely bound guest atoms rattle inside oversized structural cages and scatter heat-carrying phonons. A rattling mode is a weakly dispersive, Einstein-like vibration (i.e., behaving like an isolated oscillator at a nearly fixed energy) of such a loosely bound atom.

Layered sodium cobaltate Nax_xCoO2_2 consists of conducting triangular CoO2_2 sheets separated by layers of sodium ions. At thermoelectric compositions the Na layer contains ordered arrays of vacancies: in Na0.8_{0.8}CoO2_2 below about 290 K the vacancies form a square lattice of tri-vacancy clusters (clusters of three neighboring missing Na ions) with a large-period superstructure, and Na ions bordering these clusters are loosely bound. Voneshen et al. (Nature Materials 12, 1028, 2013) combined inelastic X-ray scattering (IXS) and inelastic neutron scattering (INS) - momentum-resolved spectroscopies that measure phonon dispersions E(Q)E(\mathbf{Q}) and, through the energy linewidth Γ\Gamma of a phonon peak, its lifetime τ=ℏ/Γ\tau = \hbar/\Gamma - with density-functional-theory (DFT) phonon calculations on a 114-atom supercell of the tri-vacancy-ordered structure. They directly observed an Einstein-like Na rattling mode near 13 meV with predominantly cubic anharmonicity, measured phonon lifetimes as a function of temperature by INS (resolution-limited below about 80 K), and, inserting the measured lifetimes into the kinetic expression κL=∑jQvjx2 τj(Q) cj(Q)\kappa^L = \sum_{j\mathbf{Q}} v_{jx}^2\,\tau_j(\mathbf{Q})\, c_j(\mathbf{Q}) (group velocity, lifetime, and heat capacity per mode), estimated the lattice thermal conductivity: rattling modes suppress κ\kappa by a factor of six relative to vacancy-free NaCoO2_2.

Three issues were left open by that study. First, the measured linewidths were used empirically but not explained: a first-principles calculation of the phonon linewidths themselves - the contributions from anharmonicity (in practice, three-phonon scattering and, where relevant, four-phonon scattering, obtained from third- and fourth-order force constants) and from electron-phonon scattering (phonons damping by creating electron-hole pairs in this metallic oxide) - remained to be done. Second, IXS at 320 K (disordered stripe phase) and 400 K (fully disordered phase) showed the dispersion softening slightly and linewidths broadening by about 50% once long-range vacancy order is lost; whether this is caused mainly by increased configurational disorder of the vacancies, by thermal lattice expansion, or by increased phonon-phonon scattering is undetermined, and the absence of long-range periodicity in those phases blocks the straightforward supercell DFT treatment used at low temperature. Third, because multi-vacancy clusters exist over a wide range of Na concentrations and most likely host rattling modes throughout the thermoelectric composition range, translating these insights into materials design requires a systematic first-principles and experimental mapping of how vacancy-cluster arrangements and sodium concentrations set the rattler energies, their anharmonicity, and the resulting κL\kappa^L suppression.

Problem Statement

For thermoelectric Nax_xCoO2_2, whose Na0.8_{0.8} member suppresses its lattice thermal conductivity by roughly a factor of six through anharmonic rattling of sodium ions in ordered multi-vacancy clusters, establish the quantitative microscopic theory of phonon damping and its consequences for thermal transport, specifically:

  1. Compute, from first principles, the phonon linewidths (lifetimes) of the vacancy-ordered Na0.8_{0.8}CoO2_2 superstructure - including anharmonic phonon-phonon scattering (three-phonon and, where relevant, four-phonon processes) and electron-phonon scattering - and decompose the measured INS/IXS linewidths into these contributions, quantitatively accounting for their observed temperature dependence and enabling a first-principles prediction of κL(T)\kappa^L(T) to replace the semi-empirical kinetic estimate that inserted measured lifetimes into κL=∑jQvjx2 τj cj\kappa^L = \sum_{j\mathbf{Q}} v_{jx}^2\,\tau_j\, c_j.

  2. Determine the dominant mechanism of the changes observed on heating above the loss of long-range vacancy order - the renormalization of the dispersion to slightly lower energies and the broadening of linewidths by about 50% - by distinguishing quantitatively among increased configurational disorder of the vacancies, thermal lattice expansion, and increased phonon-phonon scattering.

  3. Map, systematically and jointly by first-principles calculation and experiment, how different vacancy-cluster arrangements and sodium concentrations across the thermoelectric composition range 0.5≤x≤0.850.5 \le x \le 0.85 set the rattler mode energies, their anharmonicity, and the resulting suppression of κL\kappa^L, and derive from this map predictive design rules for lowering κL\kappa^L (at preserved electronic transport) in Nax_xCoO2_2 and analogous layered oxides. The map is to name explicitly the vacancy-cluster arrangements it covers - including the low-xx row (zigzag), stripe, and square tri-vacancy-cluster arrangements and any others that occur in this range - and the design rules are to relate the controllable variables (Na content xx and vacancy-cluster arrangement) to the resulting rattler energies, anharmonicity, and κL\kappa^L suppression.

These three components jointly serve one objective: replacing the empirical observation of rattler-induced thermal-conductivity suppression with a predictive, first-principles account linking vacancy superstructure, rattler anharmonicity, phonon lifetimes, and κL\kappa^L, valid at thermoelectric operating temperatures and across compositions. Each item must be resolved for Nax_xCoO2_2 itself: item 1 must deliver the linewidth decomposition together with the resulting first-principles κL(T)\kappa^L(T) prediction, item 2 the quantitative attribution among the three candidate mechanisms, and item 3 a composition- and arrangement-dependent account across the thermoelectric composition range; the extension to analogous layered oxides is optional.

References

  1. D. J. Voneshen, K. Refson, E. Borissenko, M. Krisch, A. Bosak, A. Piovano, E. Cemal, M. J. Enderle, M. J. Gutmann, M. Hoesch, M. Roger, L. Gannon, A. T. Boothroyd, S. Uthayakumar, D. G. Porter, J. P. Goff, Suppression of thermal conductivity by rattling modes in thermoelectric sodium cobaltate, Nature Materials 12, 1028-1032 (2013). DOI: 10.1038/nmat3739, https://doi.org/10.1038/nmat3739
  2. Susumu Fujii, Masato Yoshiya, Craig A. J. Fisher, Quantifying Anharmonic Vibrations in Thermoelectric Layered Cobaltites and Their Role in Suppressing Thermal Conductivity, Scientific Reports 8, 11152 (2018). DOI: 10.1038/s41598-018-29259-z, https://doi.org/10.1038/s41598-018-29259-z
  3. Nina Mattila, Antti J. Karttunen, Lattice thermal conductivity of NaCoO2 and LiCoO2 intercalation materials studied by hybrid density functional theory, Materials Research Express 7, 075502 (2020). DOI: 10.1088/2053-1591/aba3e5, https://doi.org/10.1088/2053-1591/aba3e5
  4. Nina Mattila, Antti J. Karttunen, Electronic Properties and Lattice Dynamics of LixCoO2 and NaxCoO2 (x = 0, 0.5, 1) Studied by Hybrid Density Functional Theory, Physica Status Solidi B 259, 2100665 (2022). DOI: 10.1002/pssb.202100665, https://doi.org/10.1002/pssb.202100665
  5. Ruoming Tian, Gordon J. Kearley, Dehong Yu, Chris D. Ling, Anh Pham, Jan Peter Embs, Elvis Shoko, Sean Li, Phononic Structure Engineering: the Realization of Einstein Rattling in Calcium Cobaltate for the Suppression of Thermal Conductivity, Scientific Reports 6, 30530 (2016). DOI: 10.1038/srep30530, https://doi.org/10.1038/srep30530
  6. M. Medarde, M. Mena, J. L. Gavilano, E. Pomjakushina, J. Sugiyama, K. Kamazawa, V. Yu. Pomjakushin, D. Sheptyakov, B. Batlogg, H. R. Ott, M. Månsson, F. Juranyi, 1D to 2D Na+ Ion Diffusion Inherently Linked to Structural Transitions in Na0.7CoO2, Physical Review Letters 110, 266401 (2013). DOI: 10.1103/PhysRevLett.110.266401, https://doi.org/10.1103/PhysRevLett.110.266401
  7. Fanni Juranyi, Martin Månsson, J. L. Gavilano, Mattia Mena, E. Pomjakushina, M. Medarde, Jun Sugiyama, Kazuya Kamazawa, B. Batlogg, Hans R. Ott, Tilo Seydel, Dynamics across the structural transitions at elevated temperatures in Na0.7CoO2, EPJ Web of Conferences 83, 02008 (2015). DOI: 10.1051/epjconf/20158302008, https://doi.org/10.1051/epjconf/20158302008
  8. Shangfei Wu, Hengxin Tan, Dong Wu, Mingshu Tan, Xinyu Zhou, Tianchen Hu, Tao Dong, Feng Jin, Qingming Zhang, Nanlin Wang, Temperature-driven sodium-ion dynamical-to-static crossover in the zigzag ordered phase of Na0.5CoO2, Physical Review B 114, 105102 (2026). DOI: 10.1103/4w7b-2l3l, https://doi.org/10.1103/4w7b-2l3l