Search

The temperature-dependent phonons are a generalization of interatomic force constants varying in T, which as found widespread use in computing the thermal transport of materials. A formal justification for using this combination to access thermal conductivity in anharmonic crystals, beyond the harmonic approximation and perturbation theory, is still lacking. In this work, we derive a theory of heat transport of anharmonic crystals, using the mode-coupling theory of anharmonic lattice dynamics. Starting from the Green-Kubo formula, we develop the thermal conductivity tensor based on the system's dynamical susceptibility, or spectral function. Our results account for both the diagonal and off-diagonal contributions of the heat current, with and without collective effects. We implement our theory in the TDEP package, and have notably introduced a Monte Carlo scheme to compute phonon scattering due to third- and fourth-order interactions, achieving a substantial reduction in computational cost which enables full convergence of such calculations for the first time. We apply our methodology to systems with varying regimes of anharmonicity and thermal conductivity to demonstrate its universality. These applications highlight the importance of the phonon renormalizations and their interactions beyond the harmonic order. Overall, our work advances the understanding of thermal conductivity in anharmonic crystals and provides a theoretically robust framework for predicting heat transport in complex materials.

The anharmonic lattice dynamics of oxide and halide perovskites play a crucial role in their mechanical and optical properties. Raman spectroscopy is one of the key methods used to study these structural dynamics. However, despite decades of research, existing interpretations cannot explain the temperature dependence of the observed Raman spectra. We demonstrate the non-monotonic evolution with temperature of the scattering intensity and present a model for 2nd-order Raman scattering that accounts for this unique trend. By invoking a low-frequency anharmonic feature, we are able to reproduce the Raman spectral line-shapes and integrated intensity temperature dependence. Numerical simulations support our interpretation of this low-frequency mode as a transition between two minima of a double-well potential surface. The model can be applied to other dynamically disordered crystal phases, providing a better understanding of the structural dynamics, leading to favorable electronic, optical, and mechanical properties in functional materials.

A hybrid ab initio theoretical approach for examining thermal properties in magnetic systems of unknown entropy is presented. Commonly used theoretical approaches interrogate thermal properties from Gibbs/Helmholtz free energies, which require an accurate model of magnetic interactions. The present approach avoids this requirement by instead calculating system pressure from thermally disordered microstates that properly incorporate vibrational and spin subsystems at each temperature as well as the coupling between these subsystems. In place of a specific model for magnetic interactions, the approach integrates measurements of temperature dependent magnetization of the studied material. We apply the approach to calculate phonon modes and to investigate the anomalously low thermal expansion of the classical Invar alloy, Fe_0.65Ni_0.35. The calculated phonon dispersions for Invar are in excellent agreement with measured data. The Invar thermal expansion is shown to remain small between 50 K and room temperature, consistent with the experimentally observed low thermal expansion value in this same temperature range. This anomalously small thermal expansion is directly connected to a small positive contribution from lattice thermal disorder that is nearly canceled by a large negative magnetic disorder contribution. By contrast, calculations for bcc Fe show a much larger thermal expansion, consistent with experiment, which is dominated by a large contribution from lattice thermal disorder that is reduced only slightly by a small negative contribution from that of magnetism. These findings give insights into the unusual nature of magnetism and spin-lattice coupling in Invar and Fe, and they support the presented new methodology as a complementary way to investigate thermal properties of magnetic materials.

We have found that the polarization dependence of Raman scattering in organic crystals at finite temperatures can only be described by a fourth-rank formalism. This generalization of the second-rank Raman tensor $\mathcal{R}$ stems from the effect of off-diagonal components in the crystal self-energy on the light scattering mechanism. We thus establish a novel manifestation of phonon-phonon interaction in inelastic light scattering, markedly separate from the better-known phonon lifetime.

Phonon chirality has attracted intensive attention since it breaks the traditional cognition that phonons are linear propagating bosons. This new quasiparticle property has been extensively studied theoretically and experimentally. However, characterization of the phonon chirality throughout the full Brillouin zone is still not possible due to the lack of available experimental tools. In this work, phonon dispersion and chirality of tungsten carbide were investigated by millielectronvolt energy-resolution inelastic X-ray scattering. The atomistic calculation indicates that in-plane longitudinal and transverse acoustic phonons near K and K$^\prime$ points are circularly polarized due to the broken inversion symmetry. Anomalous inelastic X-ray scattering by these circularly polarized phonons was observed and attributed to their chirality. Our results show that inelastic X-ray scattering can be utilized to characterize phonon chirality in materials and suggest that a revision to the phonon scattering function is necessary.

We combine ab initio simulations and Raman scattering measurements to demonstrate explicit anharmonic effects in the temperature dependent dielectric response of a NaCl single crystal. We measure the temperature evolution of its Raman spectrum and compare it to both a quasi-harmonic and anharmonic model. Results demonstrate the necessity of including anharmonic lattice dynamics to explain the dielectric response of NaCl, as it is manifested in Raman scattering. Our model fully captures the linear dielectric response of a crystal at finite temperatures and may therefore be used to calculate the temperature dependence of other material properties governed by it., Comment: Manuscript: 8 pages, 4 figures. Supplememtary material: 16 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.

In this study, we use the transient thermal grating optical technique \textemdash a non-contact, laser-based thermal metrology technique with intrinsically high accuracy \textemdash to investigate room-temperature phonon-mediated thermal transport in two nanoporous holey silicon membranes with limiting dimensions of 100 nm and 250 nm respectively. We compare the experimental results to ab initio calculations of phonon-mediated thermal transport according to the phonon Boltzmann transport equation (BTE) using two different computational techniques. We find that the calculations conducted within the Casimir framework, i.e. based on the BTE with the bulk phonon dispersion and diffuse scattering from surfaces, are in quantitative agreement with the experimental data, and thus conclude that this framework is adequate for describing phonon-mediated thermal transport through holey silicon membranes with feature sizes on the order of 100 nm.

The lead-free halide double perovskite class of materials offers a promising venue for resolving issues related to toxicity of Pb and long-term stability of the lead-containing halide perovskites. We present a first-principles study of the lattice vibrations in Cs$_2$AgBiBr$_6$ , the prototypical compound in this class, and show that the lattice dynamics of Cs$_2$AgBiBr$_6$ is highly anharmonic, largely in regards to tilting of AgBr$_6$ and BiBr$_6$ octahedra. Using an energy and temperature dependent phonon spectral function, we then show how the experimentally observed cubic-to-tetragonal phase transformation is caused by the collapse of a soft phonon branch. We finally reveal that the softness and anharmonicity of Cs$_2$AgBiBr$_6$ yield an ultra-low thermal conductivity, unexpected of high symmetry cubic structures.

Theoretical calculations of sound-wave velocities of materials at extreme conditions are of great importance to various fields, in particular geophysics. For example, the seismic data on sound-wave propagation through the solid iron-rich Earth's inner core have been the main source for elucidating its properties and building models. As the laboratory experiments at very high temperatures and pressures are non-trivial, ab initio predictions are invaluable. The latter, however, tend to disagree with experiment. We notice that many attempts to calculate sound-wave velocities of matter at extreme conditions in the framework of quantum-mechanics based methods have not been taking into account the effect of anharmonic atomic vibrations. We show how anharmonic effects can be incorporated into ab initio calculations and demonstrate that in particular they might be non-negligible for iron in Earth's core. Therefore, we open an avenue to reconcile experiment and ab initio theory.

The thermoelectric properties of 54 different group 4 half-Heusler (HH) alloys have been studied from first principles. Electronic transport was studied with density functional theory using hybrid functionals facilitated by the $\mathbf{k} \cdot \mathbf{p}$ method, while the temperature dependent effective potential method was used for the phonon contributions to the figure of merit $ZT$. The phonon thermal conductivity was calculated including anharmonic phonon-phonon, isotope, alloy and grain-boundary scattering. HH alloys have an ${\it XYZ}$ composition and those studied here are in the group 4-9-15 (Ti,Zr,Hf)(Co,Rh,Ir)(As,Sb,Bi) and group 4-10-14 (Ti,Zr,Hf)(Ni,Pd,Pt)(Ge,Sn,Pb). The electronic part of the thermal conductivity was found to significantly impact $ZT$ and thus the optimal doping level. Furthermore, the choice of functional was found to significantly affect thermoelectric properties, particularly for structures exhibiting band alignment features. The intrinsic thermal conductivity was significantly reduced when alloy and grain boundary scattering were accounted for, which also reduced the spread in thermal conductivity. It was found that sub-lattice disorder on the ${\it Z}$-site, i.e. the site occupied by group 14 or 15 elements, was more effective than ${\it X}$-site substitution, occupied by group 4 elements. The calculations confirmed that ZrNiSn, ZrCoSb and ZrCoBi based alloys display promising thermoelectric properties. A few other n-type and p-type compounds were also predicted to be potentially excellent thermoelectric materials, given that sufficiently high charge carrier concentrations can be achieved. This study provides insight into the thermoelectric potential of HH alloys and casts light on strategies to optimize thermoelectric performance of multicomponent alloys.

With the examples of the C $K$-edge in graphite and the B $K$-edge in hexagonal BN, we demonstrate the impact of vibrational coupling and lattice distortions on the X-ray absorption near-edge structure (XANES) in 2D layered materials. Theoretical XANES spectra are obtained by solving the Bethe-Salpeter equation of many-body perturbation theory, including excitonic effects through the correlated motion of core-hole and excited electron. We show that accounting for zero-point motion is important for the interpretation and understanding of the measured X-ray absorption fine structure in both materials, in particular for describing the $\sigma^*$-peak structure., Comment: & pages, 4 figures

The prototypical antiferroelectric PbZrO$_3$ has several unsettled questions, such as the nature of the antiferroelectric transition, possible intermediate phase and the microscopic origin of the Pbam ground state. Using first principles, we show that no phonon becomes truly soft at the cubic-to-Pbam transition temperature, and the order-disorder character of this transition is clearly demonstrated based on molecular dynamics simulations and potential energy surfaces. The out-of-phase octahedral tilting is an important degree of freedom, which can collaborate with other phonon distortions and form a complex energy landscape with multiple minima. Candidates of the possible intermediate phase are suggested based on the calculated kinetic barriers between energy minima, and the development of a first-principles-based effective Hamiltonian. The use of this latter scheme further reveals that specific bi-linear interactions between local dipoles and octahedral tiltings play a major role in the formation of the Pbam ground state, which contrasts with most of the previous explanations.

Structural phase transitions and soft phonon modes pose a longstanding challenge to computing electron-phonon (e-ph) interactions in strongly anharmonic crystals. Here we develop a first-principles approach to compute e-ph scattering and charge transport in materials with anharmonic lattice dynamics. Our approach employs renormalized phonons to compute the temperature-dependent e-ph coupling for all phonon modes, including the soft modes associated with ferroelectricity and phase transitions. We show that the electron mobility in cubic SrTiO$_{3}$ is controlled by scattering with longitudinal optical phonons at room temperature and with ferroelectric soft phonons below 200~K. Our calculations can accurately predict the temperature dependence of the electron mobility in SrTiO$_{3}$ between 150$-$300~K, and reveal the microscopic origin of its roughly $T^{-3}$ trend. Our approach enables first-principles calculations of e-ph interactions and charge transport in broad classes of crystals with phase transitions and strongly anharmonic phonons., Comment: 5 pages, 3 figures. Updated to published version

Lead chalcogenides such as PbS, PbSe, and PbTe are of interest for their exceptional thermoelectric properties and strongly anharmonic lattice dynamics. Although PbTe has received the most attention, PbSe has a lower thermal conductivity despite being stiffer, a trend that prior first-principles calculations have not reproduced. Here, we use ab-initio calculations that explicitly account for strong anharmonicity to identify the origin of this low thermal conductivity as an anomalously large anharmonic interaction, exceeding in strength that in PbTe, between the transverse optic and longitudinal acoustic branches. The strong anharmonicity is reflected in the striking observation of an intrinsic localized mode that forms in the acoustic frequencies. Our work shows the deep insights into thermal phonons that can be obtained from ab-initio calculations that are not confined to the weak limit of anharmonicity., Comment: 18 pages, 5 figures

We revisit the color-diffusion algorithm [P. C. Aeberhard et al., Phys. Rev. Lett. 108, 095901 (2012)] in nonequilibrium ab initio molecular dynamics (NE-AIMD), and propose a simple efficient approach for the estimation of monovacancy jump rates in crystalline solids at temperatures well below melting. Color-diffusion applied to monovacancy migration entails that one lattice atom (colored-atom) is accelerated toward the neighboring defect-site by an external constant force F. Considering bcc molybdenum between 1000 and 2800 K as a model system, NE-AIMD results show that the colored-atom jump rate k_{NE} increases exponentially with the force intensity F, up to F values far beyond the linear-fitting regime employed previously. Using a simple model, we derive an analytical expression which reproduces the observed k_{NE}(F) dependence on F. Equilibrium rates extrapolated by NE-AIMD results are in excellent agreement with those of unconstrained dynamics. The gain in computational efficiency achieved with our approach increases rapidly with decreasing temperatures, and reaches a factor of four orders of magnitude at the lowest temperature considered in the present study.

We have developed a method to accurately and efficiently determine the vibrational free energy as a function of temperature and volume for substitutional alloys from first principles. Taking Ti$_{1-x}$Al$_x$N alloy as a model system, we calculate the isostructural phase diagram by finding the global minimum of the free energy, corresponding to the true equilibrium state of the system. We demonstrate that the anharmonic contribution and temperature dependence of the mixing enthalpy have a decisive impact on the calculated phase diagram of a Ti$_{1-x}$Al$_x$N alloy, lowering the maximum temperature for the miscibility gap from 6560 K to 2860 K. Our local chemical composition measurements on thermally aged Ti$_{0.5}$Al$_{0.5}$N alloys agree with the calculated phase diagram., Comment: 4 pages, 5 figures, supplementary material

We present a theoretical first-principles method to calculate the free energy of a magnetic system in its high-temperature paramagnetic phase, including vibrational, electronic, and magnetic contributions. The method for calculating free energies is based on ab-initio molecular dynamics and combines a treatment of disordered magnetism using disordered local moments molecular dynamics (DLM-MD) with the temperature dependent effective potential (TDEP) method to obtain the vibrational contribution to the free energy. We illustrate the applicability of the method by obtaining the anharmonic free energy for the paramagnetic cubic and the antiferromagnetic orthorhombic phases of chromium nitride. The influence of lattice dynamics on the transition between the two phases is demonstrated by constructing the temperature-pressure phase diagram., Comment: 10 pages, 9 figures

We investigate the harmonic and anharmonic contributions to the phonon spectrum of lead telluride, and perform a complete characterization of how the anharmonic effects dominate the phonons in PbTe as temperature increases. This effect is the strongest factor in the favorable thermoelectric properties of PbTe: an optical-acoustic phonon band crossing reduces the speed of sound and the intrinsic thermal conductivity. We present the detailed temperature dependence of the dispersion relation and compare our calculated neutron scattering cross section with recent experimental measurements. We analyze the thermal resistivity's variation with temperature and clarify misconceptions about existing experimental literature. This quantitative prediction opens the way to phonon phase space engineering, to tailor the lifetimes of crucial heat carrying phonons.

The temperature dependent effective potential (TDEP) method is generalized beyond pair interactions. The second and third order force constants are determined consistently from ab initio molecular dynamics simulations at finite temperature. The reliability of the approach is demonstrated by calculations of the Mode Gr\"uneisen parameters for Si. We show that the extension of TDEP to higher order allows for an efficient calculation of the phonon life time, in Si as well as in $\epsilon$-FeSi, a system that exhibits anomalous softening with temperature.

We have developed a thorough and accurate method of determining anharmonic free energies, the temperature dependent effective potential technique (TDEP). It is based on \emph{ab initio} molecular dynamics followed by a mapping onto a model Hamiltonian that describes the lattice dynamics. The formalism and the numerical aspects of the technique are described in details. A number of practical examples are given, and results are presented, which confirm the usefulness of TDEP within \emph{ab initio} and classical molecular dynamics frameworks. In particular, we examine from first-principles the behavior of force constants upon the dynamical stabilization of body centered phase of Zr, and show that they become more localized. We also calculate phase diagram for $^4$He modeled with the Aziz \emph{et al.} potential and obtain results which are in favorable agreement both with respect to experiment and established techniques.

An accurate and easily extendable method to deal with lattice dynamics of solids is offered. It is based on first-principles molecular dynamics simulations and provides a consistent way to extract the best possible harmonic - or higher order - potential energy surface at finite temperatures. It is designed to work even for strongly anharmonic systems where the traditional quasiharmonic approximation fails. The accuracy and convergence of the method are controlled in a straightforward way. Excellent agreement of the calculated phonon dispersion relations at finite temperature with experimental results for bcc Li and bcc Zr is demonstrated.