Thermal and mechanical properties of the clathrate-II Na24Si136

Abstract

Thermal expansion, lattice dynamics, heat capacity, compressibility, and pressure stability of the intermetallic clathrate Na 24 Si 136 have been investigated by a combination of first-principles calculations and experimentation. Direct comparison of the properties of Na 24 Si 136 with those of the low-density elemental modification Si 136 provide insight into the effects of filling the silicon clathrate framework cages with Na on these properties. Calculations of the phonon dispersion only yield sensible results if the Na atoms in the large cages of the structure are displaced from the cage centers, but the exact nature of off-centering is difficult to elucidate conclusively. Pronounced peaks in the calculated phonon density of states for Na 24 Si 136 , absent for Si 136 , reflect the presence of low-energy vibrational modes associated with the guest atoms, in agreement with prior inelastic neutronscattering experiments and reflected in marked temperature dependence of the guest atom atomic displacement parameters determined by single-crystal x-ray diffraction. The bulk modulus is only weakly influenced by filling the Si framework cages with Na, whereas the phase stability under pressure is significantly enhanced. The roomtemperature linear coefficient of thermal expansion (CTE) is nearly a factor of 3 greater for Na 24 Si 136 compared to Si 136 . Negative thermal expansion (NTE), observed in Si 136 below 100 K, is noticeably absent in Na 24 Si 136 . In contrast to Si 136 , the thermal expansion behavior in Na 24 Si 136 is relatively well described by the conventional Grüneisen-Debye model in the temperature range of 10-700 K. First-principles calculations in the quasiharmonic approximation correctly predict an increase in high-temperature CTE with Na loading, although the increase is less than observed in experiment. The calculations also fail to capture the absence of NTE in Na 24 Si 136 , perhaps due to anharmonic effects and/or inadequateness of the ordered structural model.

Exact factorization of the many-body Green's function theory of electrons and nuclei

“…This renders χ and thus Ψ ≈ χΦ R useless since the nuclear wave function χ is equal to zero. The method nevertheless gives phonon spectra rather closely resembling those obtained experimentally [31] and we can compute various nuclei related physical quantities without using the nuclear wave function in position representation [32][33][34][35]. Problems will appear in this approach if we need for instance the nuclear density, as in the coupled set of equations of motion for the electronic and nuclear Green's functions beyond the BO approximation [19,22].…”