Phonon Thermal Transport in UO_{2} via Self-Consistent Perturbation Theory.

Computing thermal transport from first-principles in UO_{2} is complicated due to the challenges associated with Mott physics. Here, we use irreducible derivative approaches to compute the cubic and quartic phonon interactions in UO_{2} from first principles, and we perform enhanced thermal transport computations by evaluating the phonon Green's function via self-consistent diagrammatic perturbation theory. Our predicted phonon lifetimes at T=600 K agree well with our inelastic neutron scattering measurements across the entire Brillouin zone, and our thermal conductivity predictions agree well with previous measurements. Both the changes due to thermal expansion and self-consistent contributions are nontrivial at high temperatures, though the effects tend to cancel, and interband transitions yield a substantial contribution.

Thermal Energy Transport in Oxide Nuclear Fuel.

To efficiently capture the energy of the nuclear bond, advanced nuclear reactor concepts seek solid fuels that must withstand unprecedented temperature and radiation extremes. In these advanced fuels, thermal energy transport under irradiation is directly related to reactor performance as well as reactor safety. The science of thermal transport in nuclear fuel is a grand challenge as a result of both computational and experimental complexities. Here we provide a comprehensive review of thermal transport research on two actinide oxides: one currently in use in commercial nuclear reactors, uranium dioxide (UO2), and one advanced fuel candidate material, thorium dioxide (ThO2). In both materials, heat is carried by lattice waves or phonons. Crystalline defects caused by fission events effectively scatter phonons and lead to a degradation in fuel performance over time. Bolstered by new computational and experimental tools, researchers are now developing the foundational work necessary to accurately model and ultimately control thermal transport in advanced nuclear fuels. We begin by reviewing research aimed at understanding thermal transport in perfect single crystals. The absence of defects enables studies that focus on the fundamental aspects of phonon transport. Next, we review research that targets defect generation and evolution. Here the focus is on ion irradiation studies used as surrogates for damage caused by fission products. We end this review with a discussion of modeling and experimental efforts directed at predicting and validating mesoscale thermal transport in the presence of irradiation defects. While efforts in these research areas have been robust, challenging work remains in developing holistic tools to capture and predict thermal energy transport across widely varying environmental conditions.

Vibrations and Phase Stability in Mixed Valence Antimony Oxide.

α-Sb2O4 (cervantite) and ß-Sb2O4 (clinocervantite) are mixed valence compounds with equal proportions of SbIII and SbV as represented in the formula SbIIISbVO4. Their structure and properties can be difficult to calculate owing to the SbIII lone-pair electrons. Here, we present a study of the lattice dynamics and vibrational properties using a combination of inelastic neutron scattering, Mössbauer spectroscopy, nuclear inelastic scattering, and density functional theory (DFT) calculations. DFT calculations that account for lone-pair electrons match the experimental densities of phonon states. Mössbauer spectroscopy reveals the ß phase to be significantly harder than the α phase. Calculations with O vacancies reveal the possibility for nonstoichiometric proportions of SbIII and SbV in both phases. An open question is what drives the stability of the α phase over the ß phase, as the latter shows pronounced kinetic stability and lower symmetry despite being in the high-temperature phase. Since the vibrational entropy difference is small, it is unlikely to stabilize the α phase. Our results suggest that the α phase is more stable only because the material is not fully stoichiometric.

Phonon Spectroscopy in Antimony and Tellurium Oxides.

α-Sb2O3 (senarmontite), ß-Sb2O3 (valentinite), and α-TeO2 (paratellurite) are compounds with pronounced stereochemically active Sb and Te lone pairs. The vibrational and lattice properties of each have been previously studied but often lead to incomplete or unreliable results due to modes being inactive in infrared or Raman spectroscopy. Here, we present a study of the relationship between bonding and lattice dynamics of these compounds. Mössbauer spectroscopy is used to study the structure of Sb in α-Sb2O3 and ß-Sb2O3, whereas the vibrational modes of Sb and Te for each oxide are investigated using nuclear inelastic scattering, and further information on O vibrational modes is obtained using inelastic neutron scattering. Additionally, vibrational frequencies obtained by density functional theory (DFT) calculations are compared with experimental results in order to assess the validity of the utilized functional. Good agreement was found between DFT-calculated and experimental density of phonon states with a 7% scaling factor. The Sb-O-Sb wagging mode of α-Sb2O3 whose frequency was not clear in most previous studies is experimentally observed for the first time at ∼340 cm-1. Softer lattice vibrational modes occur in orthorhombic ß-Sb2O3 compared to cubic α-Sb2O3, indicating that the antimony bonds are weakened upon transforming from the molecular α phase to the layer-chained ß structure. The resulting vibrational entropy increase of 0.45 ± 0.1 kB/Sb2O3 at 880 K accounts for about half of the α-ß transition entropy. The comparison of experimental and theoretical approaches presented here provides a detailed picture of the lattice dynamics in these oxides beyond the zone center and shows that the accuracy of DFT is sufficient for future calculations of similar material structures.

Metallization of vanadium dioxide driven by large phonon entropy.

Phase competition underlies many remarkable and technologically important phenomena in transition metal oxides. Vanadium dioxide (VO2) exhibits a first-order metal-insulator transition (MIT) near room temperature, where conductivity is suppressed and the lattice changes from tetragonal to monoclinic on cooling. Ongoing attempts to explain this coupled structural and electronic transition begin with two alternative starting points: a Peierls MIT driven by instabilities in electron-lattice dynamics and a Mott MIT where strong electron-electron correlations drive charge localization. A key missing piece of the VO2 puzzle is the role of lattice vibrations. Moreover, a comprehensive thermodynamic treatment must integrate both entropic and energetic aspects of the transition. Here we report that the entropy driving the MIT in VO2 is dominated by strongly anharmonic phonons rather than electronic contributions, and provide a direct determination of phonon dispersions. Our ab initio calculations identify softer bonding in the tetragonal phase, relative to the monoclinic phase, as the origin of the large vibrational entropy stabilizing the metallic rutile phase. They further reveal how a balance between higher entropy in the metal and orbital-driven lower energy in the insulator fully describes the thermodynamic forces controlling the MIT. Our study illustrates the critical role of anharmonic lattice dynamics in metal oxide phase competition, and provides guidance for the predictive design of new materials.

Hybrid magnon-phonon localization enhances function near ferroic glassy states.

Ferroic materials on the verge of forming ferroic glasses exhibit heightened functionality that is often attributed to competing long- and short-range correlations. However, the physics underlying these enhancements is not well understood. The Ni45Co5Mn36.6In13.4 Heusler alloy is on the edge of forming both spin and strain glasses and exhibits magnetic field-induced shape memory and large magnetocaloric effects, making it a candidate for multicaloric cooling applications. We show using neutron scattering that localized magnon-phonon hybrid modes, which are inherently spread across reciprocal space, act as a bridge between phonons and magnons and result in substantial magnetic field-induced shifts in the phonons, triple the caloric response, and alter phase stability. We attribute these modes to the localization of phonons and magnons by antiphase boundaries coupled to magnetic domains. Because the interplay between short- and long-range correlations is common near ferroic glassy states, our work provides general insights on how glassiness enhances function.

Giant Modulation of Refractive Index from Picoscale Atomic Displacements.

It is shown that structural disorder-in the form of anisotropic, picoscale atomic displacements-modulates the refractive index tensor and results in the giant optical anisotropy observed in BaTiS3, a quasi-1D hexagonal chalcogenide. Single-crystal X-ray diffraction studies reveal the presence of antipolar displacements of Ti atoms within adjacent TiS6 chains along the c-axis, and threefold degenerate Ti displacements in the a-b plane. 47/49Ti solid-state NMR provides additional evidence for those Ti displacements in the form of a three-horned NMR lineshape resulting from a low symmetry local environment around Ti atoms. Scanning transmission electron microscopy is used to directly observe the globally disordered Ti a-b plane displacements and find them to be ordered locally over a few unit cells. First-principles calculations show that the Ti a-b plane displacements selectively reduce the refractive index along the ab-plane, while having minimal impact on the refractive index along the chain direction, thus resulting in a giant enhancement in the optical anisotropy. By showing a strong connection between structural disorder with picoscale displacements and the optical response in BaTiS3, this study opens a pathway for designing optical materials with high refractive index and functionalities such as large optical anisotropy and nonlinearity.

Implications of phonon anisotropy on thermal conductivity of fluorite oxides.

Fluorite oxides are attractive ionic compounds for a range of applications with critical thermal management requirements. In view of recent reports alluding to anisotropic thermal conductivity in this face-centered cubic crystalline systems, we perform a detailed analysis of the impact of direction-dependent phonon group velocities and lifetimes on the thermal transport of fluorite oxides. We demonstrate that the bulk thermal conductivity of this class of materials remains isotropic despite notable anisotropy in phonon lifetime and group velocity. However, breaking the symmetry of the phonon lifetime under external stimuli including boundary scattering present in nonequilibrium molecular dynamics simulations of finite size simulation cell gives rise to apparent thermal conductivity anisotropy. We observe that for accurate determination of thermal conductivity, it is important to consider phonon properties not only along high symmetry directions commonly measured in inelastic neutron or x-ray scattering experiments but also of those along lower symmetry. Our results suggests that certain low symmetry directions have a larger contribution to thermal conductivity compared to high symmetry ones.

(Mg,Mn,Fe,Co,Ni)O: A rocksalt high-entropy oxide containing divalent Mn and Fe.

High-entropy oxides (HEOs) have aroused growing interest due to fundamental questions relating to their structure formation, phase stability, and the interplay between configurational disorder and physical and chemical properties. Introducing Fe(II) and Mn(II) into a rocksalt HEO is considered challenging, as theoretical analysis suggests that they are unstable in this structure under ambient conditions. Here, we develop a bottom-up method for synthesizing Mn- and Fe-containing rocksalt HEO (FeO-HEO). We present a comprehensive investigation of its crystal structure and the random cation-site occupancy. We show the improved structural robustness of this FeO-HEO and verify the viability of an oxygen sublattice as a buffer layer. Compositional analysis reveals the valence and spin state of the iron species. We further report the antiferromagnetic order of this FeO-HEO below the transition temperature ~218 K and predict the conditions of phase stability of Mn- and Fe-containing HEOs. Our results provide fresh insights into the design and property tailoring of emerging classes of HEOs.

Assessment of empirical interatomic potential to predict thermal conductivity in ThO2and UO2.

Computing vibrational properties of crystals in the presence of complex defects often necessitates the use of (semi-)empirical potentials, which are typically not well characterized for perfect crystals. Here we explore the efficacy of a commonly used embedded-atomempirical interatomic potential for the UxTh1-xO2system, to compute phonon dispersion, lifetime, and branch specific thermal conductivity. Our approach for ThO2involves using lattice dynamics and the linearized Boltzmann transport equation to calculate phonon transport properties based on second and third order force constants derived from the empirical potential and from first-principles calculations. For UO2, to circumvent the accuracy issues associated with first-principles treatments of strong electronic correlations, we compare results derived from the empirical interatomic potential to previous experimental results. It is found that the empirical potential can reasonably capture the dispersion of acoustic branches, but exhibits significant discrepancies for the optical branches, leading to overestimation of phonon lifetime and thermal conductivity. The branch specific conductivity also differs significantly with either first-principles based results (ThO2) or experimental measurements (UO2). These findings suggest that the empirical potential needs to be further optimized for robust prediction of thermal conductivity both in perfect crystals and in the presence of complex defects.

High frequency atomic tunneling yields ultralow and glass-like thermal conductivity in chalcogenide single crystals.

Crystalline solids exhibiting glass-like thermal conductivity have attracted substantial attention both for fundamental interest and applications such as thermoelectrics. In most crystals, the competition of phonon scattering by anharmonic interactions and crystalline imperfections leads to a non-monotonic trend of thermal conductivity with temperature. Defect-free crystals that exhibit the glassy trend of low thermal conductivity with a monotonic increase with temperature are desirable because they are intrinsically thermally insulating while retaining useful properties of perfect crystals. However, this behavior is rare, and its microscopic origin remains unclear. Here, we report the observation of ultralow and glass-like thermal conductivity in a hexagonal perovskite chalcogenide single crystal, BaTiS3, despite its highly symmetric and simple primitive cell. Elastic and inelastic scattering measurements reveal the quantum mechanical origin of this unusual trend. A two-level atomic tunneling system exists in a shallow double-well potential of the Ti atom and is of sufficiently high frequency to scatter heat-carrying phonons up to room temperature. While atomic tunneling has been invoked to explain the low-temperature thermal conductivity of solids for decades, our study establishes the presence of sub-THz frequency tunneling systems even in high-quality, electrically insulating single crystals, leading to anomalous transport properties well above cryogenic temperatures.

Giant electromechanical coupling of relaxor ferroelectrics controlled by polar nanoregion vibrations.

Relaxor-based ferroelectrics are prized for their giant electromechanical coupling and have revolutionized sensor and ultrasound applications. A long-standing challenge for piezoelectric materials has been to understand how these ultrahigh electromechanical responses occur when the polar atomic displacements underlying the response are partially broken into polar nanoregions (PNRs) in relaxor-based ferroelectrics. Given the complex inhomogeneous nanostructure of these materials, it has generally been assumed that this enhanced response must involve complicated interactions. By using neutron scattering measurements of lattice dynamics and local structure, we show that the vibrational modes of the PNRs enable giant coupling by softening the underlying macrodomain polarization rotations in relaxor-based ferroelectric PMN-xPT {(1 - x)[Pb(Mg1/3Nb2/3)O3] - xPbTiO3} (x = 30%). The mechanism involves the collective motion of the PNRs with transverse acoustic phonons and results in two hybrid modes, one softer and one stiffer than the bare acoustic phonon. The softer mode is the origin of macroscopic shear softening. Furthermore, a PNR mode and a component of the local structure align in an electric field; this further enhances shear softening, revealing a way to tune the ultrahigh piezoelectric response by engineering elastic shear softening.