Search

In crystalline materials, low lattice thermal conductivity is often associated with strong anharmonicity, which can cause significant deviations from the expected Lorentzian lineshape of phonon spectral functions. These deviations, occurring in an overdamped regime, raise questions about the applicability of the Boltzmann transport equation. Furthermore, strong anharmonicity can trigger structural phase transitions with temperature, which cannot be adequately described by the standard harmonic approximation. To address these challenges, we propose a novel approach for computing the lattice thermal conductivity. Our method combines the Green-Kubo linear response theory with the stochastic self-consistent harmonic approximation. The latter allows us to describe the temperature-dependent evolution of the crystal structure, including first- and second-order phase transitions, as well as the vibrational properties in highly anharmonic materials. The Green-Kubo method considers the entire lineshapes of phonon spectral functions in the calculation of the lattice thermal conductivity, thus eliminating the questionable use of phonon lifetimes in the overdamped regime, as well as naturally including coherent transport effects. Additionally, we extend our theory to model complex dynamical lattice thermal conductivity, enhancing our understanding of time-dependent thermoreflectance experiments. As a practical application, we employ this approach to calculate the lattice thermal conductivity of CsPbBr$_3$, a complex crystal known for its anomalous thermal transport behavior with a complex phase diagram. Our method is able to determine the thermal conductivity across different phases in good agreement with experiments.

The recent report of near-ambient conditions superconductivity in a nitrogen-doped lutetium hydride has inspired a large number of experimental studies with contradictory results. We model from first principles the physical properties of the possible parent structures of the reported superconductor, LuH$_2$ and LuH$_3$. We show that only the phonon band structure of LuH$_3$ can explain the reported Raman spectra due to the presence of hydrogens at the interstitial octahedral sites. However, this structure is stabilized by anharmonicity only above 6 GPa. We find that the intriguing color change with pressure in the reported superconductor is consistent with the optical properties of LuH$_2$, which are determined by the presence of an undamped interband plasmon. The plasmon blue-shifts with pressure and modifies the color of the sample without requiring any structural phase transition. Our findings suggest that the main component in the experiments is LuH$_2$ with some extra hydrogen atoms at octahedral sites. None of LuH$_2$ and LuH$_3$ superconduct at high temperatures.

Phonon anharmonicity plays a crucial role in determining the stability and vibrational properties of high-pressure hydrides. Furthermore, strong anharmonicity can render phonon quasiparticle picture obsolete questioning standard approaches for modeling superconductivity in these material systems. In this work, we show the effects of non-Lorentzian phonon lineshapes on the superconductivity of high-pressure solid hydrogen. We calculate the superconducting critical temperature T$_\mathrm{C}$ \emph{ab initio} considering the full phonon spectral function and show that it overall enhances the T$_\mathrm{C}$ estimate. The anharmonicity-induced phonon softening exhibited in spectral functions increases the estimate of the critical temperature, while the broadening of phonon lines due to phonon-phonon interaction decreases it. Our calculations also reveal that superconductivity emerges in hydrogen in the $Cmca-12$ molecular phase VI at pressures between 450 and 500 GPa and explain the disagreement between the previous theoretical results and experiments.

Superconductivity in the new polymerized C60 fullerite clathrate doped with simple metals was investigated through density functional theory methods. The phonon dispersion curves were systematically computed for the clathrate structures containing the guest dopants (Li, Na, K, Be, Mg, Ca, Al, Ga, Ge) in one of the two distinct cages, fullerene-like C60 and sodalite-like C24, exhibited by this structure. Only four of these structures, Li@C24, Na@C24, Ga@C24 and Be@C60, are dynamically stable in the harmonic approximation. They all show superconducting behavior, although their critical temperatures are predicted to be below 2 K.

Experimental investigations of the phase transition in GeTe provide contradictory conclusions regarding the nature of the phase transition. Considering growing interest in technological applications of GeTe, settling these disputes is of great importance. To that end, we present a molecular dynamics study of the structural phase transition in GeTe using a machine-learned interatomic potential with ab-initio accuracy. First, we calculate the asymmetric shape of the radial distribution function of the nearest-neighbor bonds above the critical temperature, in agreement with previous studies. However, we show that this effect is not necessarily linked with the order-disorder phase transition and can occur as a result of large anharmonicity. Next, we study in detail the static and dynamic properties of the order parameter in the vicinity of the phase transition and find fingerprints of both order-disorder and displacive phase transition.

IV-VI materials are some of the most efficient bulk thermoelectric materials due to their proximity to soft-mode phase transitions, which leads to low lattice thermal conductivity. It has been shown that the lattice thermal conductivity of PbTe can be considerably reduced by bringing PbTe closer to the phase transition e.g. via lattice expansion. However, the effect of soft phonon modes on the electronic thermoelectric properties of such system remains unknown. Using first principles calculations, we show that the soft zone center transverse optical phonons do not deteriorate the electronic thermoelectric properties of PbTe driven closer to the phase transition via lattice expansion due to external stress, and thus enhance the thermoelectric figure of merit. We find that the optical deformation potentials change very weakly as the proximity to the phase transition increases, but the population and scattering phase space of soft phonon modes increase. Nevertheless, scattering between electronic states near the band edge and soft optical phonons remains relatively weak even very near the phase transition., Comment: 8 pages, 7 figures

The proximity to structural phase transitions in IV-VI thermoelectric materials is one of the main reasons for their large phonon anharmonicity and intrinsically low lattice thermal conductivity $\kappa$. However, the $\kappa$ of GeTe increases at the ferroelectric phase transition near $700$ K. Using first-principles calculations with the temperature dependent effective potential method, we show that this rise in $\kappa$ is the consequence of negative thermal expansion in the rhombohedral phase and increase in the phonon lifetimes in the high-symmetry phase. Negative thermal expansion increases phonon group velocities, which counteracts enhanced anharmonicity of phonon modes and boosts $\kappa$ close to the phase transition in the rhombohedral phase. A drastic decrease in the anharmonic force constants in the cubic phase increases the phonon lifetimes and $\kappa$. Strong anharmonicity near the phase transition induces non-Lorentzian shapes of the phonon power spectra. To account for these effects, we implement a novel method of calculating $\kappa$ based on the Green-Kubo approach and find that the Boltzmann transport equation underestimates $\kappa$ near the phase transition. Our findings elucidate the influence of structural phase transitions on $\kappa$ and provide guidance for design of better thermoelectric materials.

The temperature renormalization of the bulk band structure of a topological crystalline insulator, SnTe, is calculated using first principles methods. We explicitly include the effect of thermal-expansion-induced modification of electronic states and their band inversion on electron-phonon interaction. We show that the direct gap decreases with temperature, as both thermal expansion and electron-phonon interaction drive SnTe towards the phase transition to a topologically trivial phase as temperature increases. The band gap renormalization due to electron-phonon interaction exhibits a non-linear dependence on temperature as the material approaches the phase transition, while the lifetimes of the conduction band states near the band edge show a non-monotonic behavior with temperature. These effects should have important implications on bulk electronic and thermoelectric transport in SnTe and other topological insulators., Comment: 10 pages, 8 figures. Accepted for publication in Phys. Rev. B on June 8, 2020

Ferroelectric domain walls are boundaries between regions with different polarization orientations in a ferroelectric material. Using first principles calculations, we characterize all different types of domain walls forming on ($11\bar{1}$), ($111$) and ($1\bar{1}0$) crystallographic planes in thermoelectric GeTe. We find large structural distortions in the vicinity of most of these domain walls, which are driven by polarization variations. We show that such strong strain-order parameter coupling will considerably reduce the lattice thermal conductivity of GeTe samples containing domain walls with respect to single crystal. Our results thus suggest that domain engineering is a promising path for enhancing the thermoelectric figure of merit of GeTe.

GeTe is a well-known ferroelectric and thermoelectric material that undergoes a structural phase transition from a rhombohedral to the rocksalt structure at $\sim 600-700$ K. We present a first principles approach to calculate the thermal expansion of GeTe in the rhombohedral phase up to the Curie temperature. We find the minimum of the Helmholtz free energy with respect to each structural parameter in a manner similar to the traditional Gr{\"u}neisen theory, and account for the temperature dependence of elastic constants. We obtain the temperature variation of the structural parameters of rhombohedral GeTe in very good agreement with experiments. In particular, we correctly reproduce a negative volumetric thermal expansion of GeTe near the phase transition. We show that the negative thermal expansion is induced by the coupling between acoustic and soft transverse optical phonons, which is also responsible for the low lattice thermal conductivity of GeTe.