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.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Anharmonic Eigenvectors and Acoustic Phonon Disappearance in Quantum Paraelectric SrTiO 3

Pronounced anomalies in the SrTiO 3 dynamical structure factor, S(Q,E), including the disappearance of acoustic phonon branches at low temperatures, were uncovered with inelastic neutron scattering (INS) and simulations. The striking effect reflects anharmonic couplings between acoustic and optic phonons and the incipient ferroelectric instability near the quantum critical point. It is rationalized using a first-principles renormalized anharmonic phonon approach, pointing to nonlinear Ti-O hybridization causing unusual changes in real-space phonon eigenvectors, frequencies, group velocities, and scattering phase space. Lastly, our method is general and establishes how T dependences beyond the harmonic regime, assessed by INS mapping of large reciprocal-space volumes, provide real-space insights into anharmonic atomic dynamics near phase transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Kohn anomalies and phonon anharmonicity in iridium

Elemental iridium presents surprising challenges for both inelastic neutron scattering (INS) and theoretical thermal transport calculations due to its high neutron absorption cross-section and strong electron-phonon interactions, respectively. Here, in this study, we overcome these challenges to measure temperature-dependent phonon dispersion curves, compare these with calculations based on density functional theory (DFT), and ultimately examine the electron-phonon limited transport behaviors of this material. Our DFT calculations demonstrate Kohn anomalies, near the 𝐾 point of the iridium Brillouin zone, indicating coupling between electrons and phonons. Strong electron-phonon coupling can compete with anharmonic effects to determine electrical and thermal transport behaviors and make the Kohn anomalies challenging to observe. Nonetheless, our INS measurements map these anomalies and other dispersion features over the Brillouin zone from 100 to 700 K. These measurements also uncover unexpectedly large mode specific Grüneisen parameters obtained from the temperature-dependent phonon energies, highlighting strong anharmonicity in iridium. DFT-based Boltzmann transport calculations demonstrate how anharmonicity and electron-phonon couplings determine electronic and lattice transport behaviors. Furthermore, we correlate the Kohn anomalies with calculated electron-phonon nesting functions, Fermi surfaces, and DFT-derived coupling strengths. This study provides detailed insights into the temperature-dependent mode-resolved lattice dynamics and anharmonicity, transport behavior, and electron-phonon interactions.

DFT↗
DOE OSTI2022

Stability and distortion of fcc LaH 10 with path-integral molecular dynamics

The synthesis of the high-temperature superconductor LaH 10 requires pressures in excess of 100 GPa, wherein it adopts a face-centered cubic structure. Upon decompression, this structure undergoes a distortion, which still supports superconductivity, but with a lower critical temperature. Previous calculations have shown that quantum and anharmonic effects are necessary to stabilize the cubic structure, but have not resolved the low pressure distortion. Using large scale path-integral molecular dynamics enabled by a machine learned potential, we show that a rhombohedral distortion appears at sufficiently low pressures. Here, we also highlight the importance of quantum zero-point motion in stabilizing the cubic structure.

08 HYDROGEN↗
DOE OSTI2021

Onset of anharmonicity and thermal conductivity in SnSe

The anharmonicity in SnSe is investigated through the analysis of moments of 119Sn nuclear resonant inelastic x-ray scattering and ab initio molecular dynamics calculations. Experimental evidences show that the anharmonic behavior started around 300 K, substantially lower than the usually suggested structural transition at 800 K. Both theory and experiments reveal substantial lifetime broadening and frequency renormalization of the optical phonons. Additionally, thermal conductivities calculated from the temporal energy moment using the Einstein diffusion equation are in good agreement with previous experiments. The abrupt increase of the thermal power near 800 K is driven by an electronic factor and not by the enhanced anharmonicity due to structural change.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Ultrafast Suppression of the Ferroelectric Instability in KTaO 3

We use an x-ray free-electron laser to study the lattice dynamics following photoexcitation with ultrafast near-UV light (wavelength 266 nm, 50 fs pulse duration) of the incipient ferroelectric potassium tantalate, KT a O 3 . By probing the lattice dynamics corresponding to multiple Brillouin zones through the x-ray diffuse scattering with pulses from the Linac Coherent Light Source (LCLS) (wavelength 1.3 Å and < 10 fs pulse duration), we observe changes in the diffuse intensity associated with a hardening of the transverse acoustic phonon branches along Γ to X and Γ to M. Using force constants from density functional theory, we fit the quasiequilibrium intensity and obtain the instantaneous lattice temperature and density of photoexcited charge carriers. Additionally, the density functional theory calculations demonstrate that photoexcitation transfers charge from oxygen 2p derived π-bonding orbitals to Ta 5d derived antibonding orbitals, further suppressing the ferroelectric instability and increasing the stability of the cubic, paraelectric structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Anharmonicity and Phase Diagram of Magnesium Oxide in the Megabar Regime

With density functional molecular dynamics simulations, we computed the phase diagram of MgO from 50 to 2000 GPa up to 20000 K. Via thermodynamic integration (TDI), we derive the Gibbs free energies of the B1, B2, and liquid phases and determine their phase boundaries. With TDI and a pseudo-quasi-harmonic approach, we show that anharmonic effects are important and stabilize the B1 phase in particular. As a result, the B1-B2 transition boundary in the pressure-temperature plane exhibits a steep slope. We predict the B1-B2-liquid triple point to occur at approximately T= 10000 K and P = 370 GPa, which is higher in pressure than was inferred with quasi-harmonic methods alone. We predict the principal shock Hugoniot curve to enter the B2 phase stability domain but only over a very small range of parameters. Here, this may render it difficult to observe this phase with shock experiments because of kinetic effects.

3-dimensional systems↗
DOE OSTI2025

Strong anharmonicity dictates ultralow thermal conductivities of type-I clathrates

Type-I clathrate solids have attracted significant interest due to their ultralow thermal conductivities and sub- sequent promise for thermoelectric applications, yet the mechanisms underlying these properties are not well understood. Here, we extend the framework of vibrational dynamical mean-field theory (VDMFT) to calculate temperature-dependent thermal transport properties of solids using a many-body Green’s function approach. When applied to a coarse-grained description of 𝑋 8 Ga 16 Ge 30 , where 𝑋= Ba, Sr, we find that nonresonant scattering between cage acoustic modes and rattling modes leads to a reduction of acoustic phonon lifetimes and thus thermal conductivities. Moreover, we find that the moderate temperature dependence of conductivities above 300 K, which is consistent with experimental measurements, cannot be reproduced by textbook perturbation theory calculations, which predict a 𝑇 −1 dependence. Therefore, we suggest that nonperturbative anharmonic effects, including four- and higher-phonon scattering processes, are responsible for the ultralow thermal conductivities of type-I clathrates.

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

Anharmonic phonon behavior via irreducible derivatives: Self-consistent perturbation theory and molecular dynamics

Cubic phonon interactions are now regularly computed from first principles, and the quartic interactions have begun to receive more attention. Given this realistic anharmonic vibrational Hamiltonian, the classical phonon Green's function can be precisely measured using molecular dynamics, which can then be used to rigorously assess the range of validity for self-consistent diagrammatic approaches in the classical limit. Here we use the bundled irreducible derivative approach to efficiently and precisely compute the cubic and quartic phonon interactions of CaF 2 , systematically obtaining the vibrational Hamiltonian purely in terms of irreducible derivatives. We demonstrate that the 4-phonon sunset diagram has an important contribution to the optical phonon linewidths beyond T = 500 K. Reasonable results are obtained even at T = 900 K when performing self-consistency using the 4-phonon loop diagram and evaluating the 3-phonon bubble and 4-phonon sunset diagrams post-self-consistency. Further improvements are obtained by performing quasiparticle perturbation theory, where both the 4-phonon loop and the real part of the 3-phonon bubble are employed during self-consistency. Finally, our irreducible derivative approach to self-consistent perturbation theory is a robust tool for studying anharmonic phonons in both the quantum and classical regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Thermal conductivity of CaSiO 3 perovskite at lower mantle conditions

Thermal conductivity (κ) of mantle minerals is key to understanding dynamics in the deep Earth. It controls the style of mantle convection and the timescale of cooling both the mantle and the core. Cubic CaSiO 3 perovskite (CaPv) is the third most abundant mineral in the lower mantle (7 vol %). Despite its importance, no theoretical or experimental estimate of CaPv's κ is available. Theoretical investigations of its properties are challenging because of its strong anharmonicity. Experimental measurements at relevant pressures and temperatures are equally challenging. Here we present ab initio results for CaPv's κ obtained using the phonon quasiparticle approach to address its strong anharmonicity. We also offer experimental measurements of κ up to 67 GPa and 1950 K. Predictions and measurements are in good agreement and reveal a surprisingly large κ for cubic CaPv that can be explained on the basis of its high crystal structure symmetry. Despite its relatively low abundance, CaPv's κ increases the lower mantle κ by ~10 %, if accounted for. κ of mantle regions enriched in subducted crustal materials will be more strongly impacted.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Phonon second harmonic generation in NaBr studied by inelastic neutron scattering and computer simulation

The phenomenon of second harmonic generation (SHG) was found for phonons in anharmonic NaBr by inelastic neutron scattering. The temperature dependence of this phonon SHG was measured from 300 K to 650 K. At 300 K the second harmonic (SH) is seen as a high-energy branch around 33 meV, nearly independent of $\overrightarrow{Q}$. The temperature effective potential (TDEP) method and classical molecular dynamics (MD) simulation with machine learning interatomic potential were able to reproduce the SH, and showed that SHG occurs with the flat transverse optical (TO) phonon branch. A classical model of a nonlinear medium explains the intensity and lifetime of the SH, compared to those of the TO modes. Also successful was a quantum model based on the Heisenberg-Langevin equation for interacting phonons coupled to a thermal bath, which also predicts a spectral distribution of the SH. In conclusion, the measured temperature dependence of the intensity of the second harmonic showed that it follows the Planck distribution of a one-phonon quasiparticle, and not two TO phonons.

36 MATERIALS SCIENCE↗
DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Three-Dimensional Lead Iodide Perovskitoid Hybrids with High X-ray Photoresponse

Large organic A cations cannot stabilize the 3D perovskite AMX 3 structure because they cannot be accommodated in the cubo-octhedral cage (do not follow the Goldschmidt tolerance factor rule), and they generally template low-dimensional structures. In this paper, we report that the large dication aminomethylpyridinium (AMPY) can template novel 3D structures which resemble conventional perovskites. They have the formula (xAMPY)M 2 I 6 (x = 3 or 4, M = Sn 2+ or Pb 2+ ) which is double of the AMX3 formula. However, because of the steric requirement of the Goldschmidt tolerance factor rule, it is impossible for (xAMPY)M 2 I 6 to form proper perovskite structures. Instead, a combination of corner-sharing and edge-sharing connectivity is adopted in these compounds leading to the new 3D structures. DFT calculations reveal that the compounds are indirect band gap semiconductors with direct band gaps presenting at slightly higher energies and dispersive electronic bands. The indirect band gaps of the Sn and Pb compounds are ~1.7 and 2.0 eV, respectively, which is slightly higher than the corresponding AMI 3 3D perovskites. The Raman spectra for the compounds are diffuse, with a broad rising central peak at very low frequencies around 0 cm –1 , a feature that is characteristic of dynamical lattices, high anharmonicity, and dissipative vibrations very similar to the 3D AMX 3 perovskites. Devices of (3AMPY)Pb 2 I 6 crystals exhibit clear photoresponse under ambient light without applied bias, reflecting a high carrier mobility (μ) and long carrier lifetime (τ). The devices also exhibit sizable X-ray generated photocurrent with a high μτ product of ~1.2×10 -4 cm 2 /V and a X-ray sensitivity of 207 μC∙Gy -1 ∙cm -2 .

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

A unified understanding of minimum lattice thermal conductivity

Here, we propose a first-principles model of minimum lattice thermal conductivity ($κ^{min}_L$) based on a unified theoretical treatment of thermal transport in crystals and glasses. We apply this model to thousands of inorganic compounds and find a universal behavior of $κ^{min}_L$ in crystals in the high-temperature limit: The isotropically averaged $κ^{min}_L$ is independent of structural complexity and bounded within a range from ~0.1 to ~2.6 W/(m K), in striking contrast to the conventional phonon gas model which predicts no lower bound. We unveil the underlying physics by showing that for a given parent compound, $κ^{min}_L$ is bounded from below by a value that is approximately insensitive to disorder, but the relative importance of different heat transport channels (phonon gas versus diffuson) depends strongly on the degree of disorder. Moreover, we propose that the diffuson-dominated $κ^{min}_L$ in complex and disordered compounds might be effectively approximated by the phonon gas model for an ordered compound by averaging out disorder and applying phonon unfolding. With these insights, we further bridge the knowledge gap between our model and the well-known Cahill–Watson–Pohl (CWP) model, rationalizing the successes and limitations of the CWP model in the absence of heat transfer mediated by diffusons. Finally, we construct graph network and random forest machine learning models to extend our predictions to all compounds within the Inorganic Crystal Structure Database (ICSD), which were validated against thermoelectric materials possessing experimentally measured ultralow κ L . Our work offers a unified understanding of $κ^{min}_L$, which can guide the rational engineering of materials to achieve .

42 ENGINEERING↗
DOE OSTI2020

Anharmonic phonon damping enhances the T c of BCS-type superconductors

In this study, a theory of superconductivity is presented where the effect of anharmonicity, as entailed in the acoustic, or optical, phonon damping, is explicitly considered in the pairing mechanism. The gap equation is solved including diffusive Akhiezer damping for longitudinal acoustic phonons or Klemens damping for optical phonons, with a damping coefficient which, in either case, can be directly related to the Grüneisen parameter and hence to the anharmonic coefficients in the interatomic potential. The results show that the increase of anharmonicity has a strikingly nonmonotonic effect on the critical temperature T c . The optimal damping coefficient yielding maximum T c is set by the velocity of the bosonic mediator. This theory may open up unprecedented opportunities for material design where T c may be tuned via the anharmonicity of the interatomic potential, and presents implications for the superconductivity in the recently discovered hydrides, where anharmonicity is very strong and for which the anharmonic damping is especially relevant.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ab initio anharmonic thermodynamic properties of cubic Ca Si O 3 perovskite

We present an ab initio study of the thermodynamic properties of cubic CaSiO 3 perovskite (CaPv) over the pressure and temperature range of the Earth's lower mantle. We compute the anharmonic phonon dispersions throughout the Brillouin zone by utilizing the phonon quasiparticle approach, which characterizes the intrinsic temperature dependence of phonon frequencies and, in principle, captures full anharmonicity. Such temperature-dependent phonon dispersions are used to calculate ab initio free energy in the thermodynamic limit (N → ∞) within the framework of the phonon gas model. Accurate free energy calculations enable us to investigate cubic CaPv's thermodynamic properties, e.g., thermal expansivity, Grüneisen parameter, bulk modulus, heat capacity and thermal equation of state, where anharmonic effects are demonstrated. Furthermore, the present methodology provides an important theoretical approach for exploring phase boundaries, thermodynamic, and thermoelastic properties of strongly anharmonic materials at high pressures and temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Relaxation timescales and electron-phonon coupling in optically pumped $\mathrm{YBa_2Cu_3O_{6+\mathcal{x}}}$ revealed by time-resolved Raman scattering

Time-resolved measurements provide a new way to disentangle complex interactions in quantum materials due to their different timescales. We used pump-probe Raman scattering to investigate the apical oxygen vibration in $\mathrm{YBa_2Cu_3O_{6+\mathcal{x}}}$ under nonequilibrium conditions. Time dependence of the phonon population demonstrated strong electron-phonon coupling. Most importantly, the phonon shifts to a higher energy due to transient smearing of the Fermi surface in a remarkable agreement with diagrammatic theory. Here we also discuss insights into photoinduced superconductivity reported at lower doping that follow from these results.

36 MATERIALS SCIENCE↗