Search NASA⌕ Search

SEARCH · Search NASA

Results for “Anharmonic lattice dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

124 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Matryoshka phonon twinning in α-GaN

Understanding lattice dynamics is crucial for effective thermal management in electronic devices because phonons dominate thermal transport in most semiconductors. α -GaN has become a focus of interest as one of the most important third-generation power semiconductors, however, the knowledge on its phonon dynamics remains limited. Here we show a Matryoshka phonon dispersion of α -GaN with the complementary inelastic X-ray and neutron scattering techniques and the first-principles calculations. Such Matryoshka twinning throughout the basal plane of the reciprocal space is demonstrated to amplify the anharmonicity of the related phonons through creating abundant three-phonon scattering channels and cutting the lifetime of affected modes by more than 50%. Such phonon topology contributes to reducing the in-plane thermal transport, thus the anisotropic thermal conductivity of α -GaN. The results not only have implications for engineering the thermal performance of α -GaN, but also offer valuable insights on the role of anomalous phonon topology in thermal transport of other technically semiconductors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

High-throughput computational framework for high-order anharmonic thermal transport in cubic and tetragonal crystals

Accurate first-principles prediction of lattice thermal conductivity (κL) remains challenging in identifying materials with extreme thermal behavior. While the harmonic approximation with three-phonon scattering (HA + 3ph) is now routine, reliable κL prediction often requires higher-order anharmonic effects, including self-consistent phonon renormalization, three- and four-phonon scattering, and off-diagonal heat flux (SCPH + 3, 4ph + OD). We present a state-of-the-art high-throughput workflow that unifies these effects and apply it to 773 cubic and tetragonal crystals spanning diverse chemistries and structures. From 562 dynamically stable compounds, we assess the hierarchical impacts of higher-order anharmonicity. For around 60% of materials, HA + 3ph predictions closely match those from SCPH + 3, 4ph + OD. SCPH generally increases κ L , by over 8 times in extreme cases, whereas four-phonon scattering universally suppresses κ L , sometimes to 15% of the HA + 3ph value. Off-diagonal contributions are negligible in high-κ L systems but can rival diagonal terms in highly anharmonic low-κ L compounds. We highlight four case studies, Rb 2 TlAlH 6 , Cu 3 VSe 4 , CuBr, and KTlCl 4 , that exhibit distinct extreme behaviors. This work delivers not only a robust workflow for high-fidelity κL dataset but also a quantitative framework to determine when higher-order effects are essential. The hierarchy of κ L results, from the HA + 3ph to SCPH + 3, 4ph + OD level, offers a scalable, interpretable route to discovering next-generation extreme thermal materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Robust cycloid crossover driven by anisotropy in the skyrmion host GaV 4 S 8

We report on the anomalous magnetization dynamics of the cycloidally modulated spin textures under the influence of uniaxial anisotropy in multiferroic GaV 4 S 8 . The temperature and field dependence of the linear ac susceptibility [$χ'_{1ω}$ (T, H)], AC magnetic loss [$χ''_{1ω}$ (T, H)], and nonlinear AC magnetic response [$M_{3ω}$ (T, H)] are examined across the magnetic phase diagram in the frequency range f = 10 - 10 000 Hz. According to recent theory, skyrmion vortices under axial crystal symmetry are confined along specific orientations, resulting in enhanced robustness against oblique magnetic fields and altered spin dynamics. We characterize the magnetic response of each spin texture and find that the dynamic rigidity of the Néel skyrmion lattice appears enhanced compared to Bloch-type skyrmions in cubic systems, even in the multidomain state. Anomalous M 3 ω and strong dissipation emerge over the same phase regime where strong variations in the cycloid pitch were observed on lowering temperature in recent small-angle neutron-scattering experiments [White et al., Phys. Rev. B 97, 020401(R) (2018)]. Here, we show that strong anisotropy also drives an extended crossover of the zero-field cycloid texture in GaV 4 S 8 . The frequency dependence of these dynamic signatures is consistent with that of a robust anharmonic spin texture exhibiting a correlated domain arrangement. The results underpin the essential role of magnetic anisotropy in enhancing the rigidity of topological spin textures for diverse applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Quantum stabilization and flat hydrogen-based bands of nitrogen-doped lutetium hydride

We explore the electronic and structural properties of Lu-H-N structures with specific N,H ordering due to their remarkably narrow hydrogen-based bands at the Fermi level. Although exhibits a classical instability persisting up to 17 GPa, it is anharmonically stable near ambient pressure when accounting for quantum nuclear effects. The presence of flat bands near is understood to arise from destructive quantum interference between and surrounding orbitals, with certain types of defects leaving the flat bands unaffected. The results suggest that if this system is superconducting, there is an optimal pressure near ambient where the superconducting is maximized by anharmonically stabilized low-frequency and nonadiabatically coupled high-frequency hydrogen modes. Despite the metastability of this structure, its dynamical stability when modeled beyond a classical harmonic approach suggests the narrow bands near do not cause structural distortions in certain Lu-H-N stoichiometries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Temperature-dependent lattice dynamics in iridium

The characterization of simple elemental systems is key to benchmarking first-principles modeling of electronic and vibrational behaviors of materials. A large body of literature has been built for most elemental systems in the periodic table. However, surprisingly little neutron work has been performed to probe the vibrational properties of iridium, likely due to its large neutron absorption cross section. Nonetheless, iridium is of significant scientific and technological interest due to large relativistic electron effects and electron-phonon coupling, particularly in strongly correlated iridate compounds. In this paper, we report temperature-dependent inelastic neutron scattering measurements of the vibrational properties of iridium, from which we extract key thermodynamic properties. To overcome the challenge of the large neutron absorption of iridium, we developed a simple post-processing correction procedure. The measured densities of phonon states compare well with quasiharmonic density functional theory calculations, although the obtained experimental phonon Grüneisen parameters are much larger than expected, reaching as high as γ=4.5, indicating substantial anharmonicity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Origin of Intrinsically Low Thermal Conductivity in a Garnet-Type Solid Electrolyte: Linking Lattice and Ionic Dynamics with Thermal Transport

Understanding thermal transport in solid electrolytes is essential for improving the performance, reliability, and safety of all-solid-state batteries. Garnet-type lithium-ion conductors are promising candidates for solid electrolytes, yet their thermal-transport mechanisms remain poorly understood. Here, we connect the lattice and ion dynamics of single-crystal garnet-type Li 6.5 La 3 Zr 1.5 Ta 0.5 O 12 to its intrinsically low thermal conductivity. Our study reveals that the single crystals grown by the floating-zone method exhibit remarkably low glasslike thermal conductivity. Using first-principles calculations and inelastic-neutron-scattering measurements, we identify both the acoustic and numerous optical phonon modes, which stem from the complex crystal structure of the material. Notably, a low-energy optical branch exhibits an avoided crossing with acoustic phonons near 7 meV. These optical modes can enhance the scattering of heat-carrying acoustic phonons and reduce thermal conductivity. Furthermore, the calculated Grüneisen parameters are large, especially for the vibrational modes around 6 meV, indicating strong anharmonicity, with a noticeable contribution from lithium-ion vibrations. A two-channel thermal-transport model is employed to describe the weak temperature dependence of the thermal conductivity, which can be attributed to the substantial contribution of diffuson transport facilitated by the abundance of optical phonons and intrinsic anharmonicity. These results offer valuable insights into the thermal transport in a broad class of ionic conductors of interest for energy conversion and storage applications.

Ab initio calculations↗
DOE OSTI2020

Response of the mode Grüneisen parameters with anisotropic compression: A pressure and temperature dependent Raman study of β-Sn

The lattice dynamic response of body-centered tetragonal beta-Sn (I4 1 /amd) under high-pressure and temperature conditions is determined using experimental optical vibration modes. Here, Raman scattering is used to map the phase stability region of beta-Sn to perform mode Gruneisen analysis, and we demonstrate the necessity of an optical intensity calibration for Raman thermometry. The Gruneisen tensor is evaluated along a set of isotherms to address shortcomings of single-mode Gruneisen parameters with respect to anisotropic deformations of this tetragonal structured soft metal. The changes observed here in the Gruneisen tensor as a function of temperature are related to anharmonicity and denote potential criteria for the onset of premelting.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Anomalous thermal transport under high pressure in boron arsenide

High pressure represents extreme environments and provides opportunities for materials discovery. Thermal transport under high hydrostatic pressure has been investigated for more than 100 years and all measurements of crystals so far have indicated a monotonically increasing lattice thermal conductivity. Here we report in situ thermal transport measurements in the newly discovered semiconductor crystal boron arsenide, and observe an anomalous pressure dependence of the thermal conductivity. We use ultrafast optics, Raman spectroscopy and inelastic X-ray scattering measurements to examine the phonon bandstructure evolution of the optical and acoustic branches, as well as thermal conductivity under varied temperatures and pressures up to 32 gigapascals. Using atomistic theory, we attribute the anomalous high-pressure behaviour to competitive heat conduction channels from interactive high-order anharmonicity physics inherent to the unique phonon bandstructure. Our study verifies ab initio theory calculations and we show that the phonon dynamics—resulting from competing three-phonon and four-phonon scattering processes—are beyond those expected from classical models and seen in common materials. In conclusion, this work uses high-pressure spectroscopy combined with atomistic theory as a powerful approach to probe complex phonon physics and provide fundamental insights for understanding microscopic energy transport in materials of extreme properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Local structure, bonding, and asymmetry of $\mathrm{((NH_2)_2CH)PbBr_3}$, $\mathrm{CsPbBr_3}$, and $\mathrm{(CH_3MH_3)PbBr_3}$

Here we report local structure measurements for CsPbBr 3 and ((NH 2 ) 2 CH)PbBr 3 (FAPbBr 3 ) and compare them with recent results for (CH 3 NH 3 )PbBr 3 (MAPbBr 3 ). The Pb-Br bonding is similar for all three systems; the effective spring constants, κ, are comparable (ranging from 1.20 to 1.95 eV/Å 2 ), but small in magnitude indicating very soft materials. However, there are also important differences between the three systems. Static disorder is very small for CsPbBr3 but increases somewhat with the size of the organic molecular ions MA + and FA + . At room temperature, dynamic disorder dominates in all compounds. The thermal disorder of the Pb-Br pair distribution function (PDF), i.e., the Debye-Waller factor σ 2 follows a correlated Debye or Einstein model up to 300 K in CsPbBr 3 (orthorhombic phase), but for FAPbBr 3 and MAPbBr 3 , there is a break in the σ 2 (T) curve at the orthorhombic-tetragonal transition (o-t) near 150 K, indicating a small change in the spring constant κ. κ increases for MAPbBr 3 but decreases in FAPbBr 3 at this transition. These changes are attributed to changes in the H-bonding between Br - and MA + or FA + at this transition, as a result of librations or rotations of these molecular cations. In addition, the Pb-Br PDF becomes asymmetric at a relatively low temperature for FAPbBr 3 and MAPbBr 3 , while this effect is significantly smaller for CsPbBr 3 . Finally, we address the question of a model to explain the asymmetric PDF. Two main models are discussed in the literature, an anharmonic pair potential and a split-pair distribution, possibly driven by the presence of a lone pair on the Pb ion. We show that the fourth cumulant C 4 can differentiate between these two models and other possible models. Experimentally C 4 is positive at 250 K and above, for all three systems and that is inconsistent with a split-peak model, for which C 4 is negative for splittings larger than 0.12 Å.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Phonon-induced rotation of the electronic nematic director in superconducting Bi 2 Se 3

The doped topological insulator A x Bi 2 Se 3 , with A = {Cu,Sr,Nb}, becomes a nematic superconductor below T c ≈ 3–4 K. The associated electronic nematic director is described by an angle α and is experimentally manifested in the elliptical shape of the in-plane critical magnetic field H c2 . Because of the threefold rotational symmetry of the lattice, α is expected to align with one of three high-symmetry directions corresponding to the in-plane nearest-neighbor bonds, consistent with a Z 3 -Potts nematic transition. Here, we show that the nematic coupling to the acoustic phonons, which makes the nematic correlation length tend to diverge along certain directions only, can fundamentally alter this phenomenology in trigonal lattices. Compared to hexagonal lattices, the former possesses a sixth independent elastic constant c14 due to the fact that the in-plane shear strain doublet (ε xx – ε yy , – 2ε xy ) and the out-of-plane shear strain doublet (2ε yz , – 2ε zx ) transform as the same irreducible representation. We find that, when c 14 overcomes a threshold value, which is expected to be the case in doped Bi 2 Se 3 , the nematic director α unlocks from the high-symmetry directions due to the competition between the quadratic phonon-mediated interaction and the cubic nematic anharmonicity. This implies the breaking of the residual in-plane twofold rotational symmetry (C 2x ), resulting in a triclinic phase. Furthermore, we discuss the implications of these findings for the structure of nematic domains, for the shape of the in-plane H c2 in A x Bi 2 Se 3 , and for the presence of nodes inside the superconducting state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Thermal conductivity and its relation to atomic structure for symmetrical tilt grain boundaries in silicon

We perform a systematic study of thermal resistance and conductance of tilt grain boundaries (GBs) in Si using classical molecular dynamics. The GBs studied are naturally divided into three groups according to the structural units forming the GB core. In this work we find that, within each group, the GB thermal conductivity strongly correlates with the excess GB energy. All three groups predict nearly the same GB conductivity extrapolated to the high-energy limit. This limiting value is close to the thermal conductivity of amorphous Si, suggesting similar heat transport mechanisms. While the lattice thermal conductivity decreases with temperature, the GB conductivity slightly increases. However, at high temperatures it turns over and starts decreasing if the GB structure undergoes a premelting transformation. Analysis of vibrational spectra of GBs resolved along different directions sheds light on the mechanisms of their thermal resistance. The existence of alternating tensile and compressive atomic environments in the GB core gives rise to localized vibrational modes, frequency gaps creating acoustic mismatch with lattice phonons, and anharmonic vibrations of loosely bound atoms residing in open atomic environments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Anharmonic theory of superconductivity in the high-pressure materials

Electron-phonon superconductors at high pressures have displayed the highest values of critical superconducting temperature Tc on record, now rapidly approaching room temperature. Despite the importance of high-P superconductivity in the quest for room-temperature superconductors, a mechanistic understanding of the effect of pressure and its complex interplay with phonon anharmonicity and superconductivity is missing, as numerical simulations can bring only system-specific details, clouding out key players controlling the physics. In this work, we develop a minimal model of electron-phonon superconductivity under an applied pressure which takes into account the anharmonic decoherence of the optical phonons. We find that T c behaves nonmonotonically as a function of the ratio Γ/ω 0 , where Γ is the optical phonon damping and ω 0 is the optical phonon energy at zero pressure and momentum. Optimal pairing occurs for a critical ratio Γ/ω 0 when the phonons are on the verge of decoherence (“diffusonlike” limit). Our framework gives insights into recent experimental observations of T c as a function of pressure in the complex BCS material TlInTe 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Efficient perturbation-tracking method for directly probing the spectral phonon properties from molecular dynamics simulations

Existing methods for directly extracting the spectral phonon properties from molecular dynamics (MD) simulations, like the normal mode analysis (NMA) and spectral energy density analysis, all require a very long simulation time to produce reliable results with good convergence. So far, these methods are mainly applied in studies using small systems and with empirical potentials, as the heavy computational load has greatly hindered their further applications. Here we propose a perturbation-tracking (PT) method for directly probing the mode-wise phonon anharmonic frequencies and lifetimes. We show that results obtained from our method are in excellent agreement with those from the conventional NMA approach, using Si as the model material system. Comparing with the NMA approach, the PT method offers a greater accuracy and significant improvement of efficiency. It takes an average of two orders of magnitude and up to three orders of magnitude less simulation time to obtain the same lifetime result of a phonon mode with intermediate to high accuracy. Meanwhile, our method preserves all the dynamics of probed phonon mode from a particular state, which means it is capable of studying the transient thermal transport processes in a nonequilibrium system. Besides the exceptional efficiency, our method also comes with freedom to choose to probe only those modes of interest. This makes it ideal for use with large systems and in computationally demanding applications, such as ab initio MD simulations. Moreover, the PT method we propose here is very straightforward and easy to implement.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Intrinsic anharmonicity and thermal properties of ultralow thermal conductivity Ba 6 Sn 6 Se 13

Ultralow thermal conductivity materials continue to be of great interest for technologically important applications such as thermal insulators and thermoelectrics. Often, such materials possess constituents with extensive dynamic disorder, or “rattlers,” or structural disorder such as metallic glasses. Nevertheless, targeted crystalline bonding order and anharmonicity can provide effective means for the scattering of the phonons in certain materials. Our combined experimental and theoretical investigation of Ba 6 Sn 6 Se 13 reveals such a material, with ultralow thermal conductivity measured over a large temperature range. Optic phonon modes hybridize with acoustic modes at relatively low energies resulting in strong acoustic-optic scattering and limited phonon lifetimes. Moreover, our investigations reveal various avoided crossings at low energies that contribute to phonon scattering. These are derived from chiral phase symmetries and may result in nontrivial topological behaviors for phonon band crossings in this chiral material. This work contributes to the ongoing research on low thermal conductivity materials and the underlying mechanisms that affect their thermal behaviors. Manipulation of these effects may provide pathways for enhancement of material properties for targeted thermal applications using similar materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electronic structure and lattice dynamics of 1⁢𝑇−VSe 2 : Origin of the three-dimensional charge density wave

To characterize in detail the charge density wave (CDW) transition of 1⁢𝑇−VSe 2 , its electronic structure and lattice dynamics are comprehensively studied by means of x-ray diffraction, muon spectroscopy, angle resolved photoemission (ARPES), diffuse and inelastic x-ray scattering, and state-of-the-art first-principles density functional theory calculations. Resonant elastic x-ray scattering does not show any resonant enhancement at either V or Se, indicating that the CDW peak at the 𝐾 edges describes a purely structural modulation of the electronic ordering. ARPES experiments identify (i) a pseudogap at 𝑇 > 𝑇⁢ CDW , which leads to a depletion of the density of states in the ML-M'L' plane at 𝑇 < 𝑇⁢ CDW , and (ii) anomalies in the electronic dispersion reflecting a sizable impact of phonons on it. A diffuse scattering precursor, characteristic of soft phonons, is observed at room temperature (RT) and leads to the full collapse of the low-energy phonon (𝜔 1 ) with propagation vector (0.25 0 −0.3) r.l.u. Here, we show that the frequency and linewidth of this mode are anisotropic in momentum space, reflecting the momentum dependence of the electron-phonon interaction (EPI), hence demonstrating that the origin of the CDW is, to a much larger extent, due to the momentum dependent EPI with a small contribution from nesting. The pressure dependence of the 𝜔 1 soft mode remains nearly constant up to 13 GPa at RT, with only a modest softening before the transition to the high-pressure monoclinic C2/m phase. The wide set of experimental data is well captured by our state-of-the art first-principles anharmonic calculations with the inclusion of van der Waals corrections in the exchange-correlation functional. The comprehensive description of the electronic and dynamical properties of VSe 2 reported here adds important pieces of information to the understanding of the electronic modulations in the family of transition-metal dichalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Generalized quasiharmonic approximation via space group irreducible derivatives

The quasiharmonic approximation (QHA) is the simplest nontrivial approximation for interacting phonons under constant pressure, bringing the effects of anharmonicity into temperature-dependent observables. Nonetheless, the QHA is often implemented with additional approximations due to the complexity of computing phonons under arbitrary strains, and the generalized QHA, which employs constant stress boundary conditions, has not been completely developed. In this work we formulate the generalized QHA, providing a practical algorithm for computing the strain state and other observables as a function of temperature and true stress. We circumvent the complexity of computing phonons under arbitrary strains by employing irreducible second-order displacement derivatives of the Born-Oppenheimer potential and their strain dependence, which are efficiently and precisely computed using the lone irreducible derivative approach. We formulate two complementary strain parametrizations: a discretized strain grid interpolation and a Taylor series expansion in symmetrized strain. We illustrate our approach by evaluating the temperature and pressure dependence of select elastic constants and the thermal expansion in thoria (ThO 2 ) using density functional theory with three exchange-correlation functionals. The QHA results are compared to our measurements of the elastic constant tensor using time-domain Brillouin scattering and inelastic neutron scattering. Our irreducible derivative approach simplifies the implementation of the generalized QHA, which will facilitate reproducible, data-driven applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗