Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “anharmonicity”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
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 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 OSTI2025

First-principles thermodynamics of Al 10 ⁢V: An analytical treatment of localized anharmonic modes

Many complex intermetallic structures possess cagelike environments that can host additional guest atoms. In Al 10 ⁢V, these atoms give rise to low-frequency, localized vibrations (Einstein modes) that dominate the thermodynamic response at low temperature. They become imaginary under volume expansion as temperature rises, invalidating the harmonic approximation. Here, we develop a framework to incorporate these strongly anharmonic vibrational modes into first-principles thermodynamic calculations. By explicitly modeling the cage potential and solving the associated Schrödinger equation numerically, we compute the full anharmonic free energy contribution and demonstrate its impact on the thermodynamic behavior of Al 10⁢ V. This allows us to examine structures with different cage fillings and construct the Al-V phase diagram in the relevant composition range. Our results reproduce key experimental signatures, including the anomalous rise in the thermal expansion coefficient and heat capacity at low temperatures, and reveal that the presence and the extent of cage filling by guest atoms is essential to stabilizing the Al 10 ⁢V phase at elevated temperatures.

anharmonic lattice dynamics↗
DOE OSTI2023

Structure and Anharmonicity of α- and β -Sb 2 O 3 at Low Temperature

Antimony oxides are important materials for catalysis and flame-retardant applications. The two most common phases, α-Sb 2 O 3 (senarmontite) and β-Sb 2 O 3 (valentinite), have been studied extensively. Specific focus has been placed recently on their lattice dynamics properties and how they relate to the α-β phase transformation and their potential anharmonicity. However, there has not been any direct investigation of anharmonicity in these systems, and a surprising lack of low-temperature structural information has prevented further study. Here, we report the powder neutron diffraction data of both phases of Sb 2 O 3 , as well as structural information. α-Sb 2 O 3 behaved as expected, but β-Sb 2 O 3 revealed a small region of zero thermal expansion along the c axis. Additionally, while the β phase matched well with reported atomic displacement parameters, the α phase displayed a marked deviation. This data will enable further investigations into these systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Strongly Anharmonic Phonons and Their Role in Superionic Diffusion and Ultralow Thermal Conductivity of Cu 7 PSe 6

We report the quest for advanced superionic materials requires understanding their complex atomic dynamics, but detailed studies of the interplay between lattice vibrations and ionic diffusion remain scarce. Here inelastic and quasielastic neutron scattering measurements in the superionic argyrodite Cu 7 PSe 6 are reported, combined with molecular dynamics (MD) based on ab initio and machine-learned potentials, providing critical insights into the atomistic mechanisms underlying fast ion conduction. The results reveal how long-range Cu diffusion is limited by intercluster hopping, controlled by selective anharmonic phonons of the crystalline framework. Further, the Green–Kubo simulations reproduce the ultralow lattice thermal conductivity and identify contributions from mobile ions, phonons, and their cross-correlations. The mode resolved analysis shows that the thermal conductivity is dominated by low-energy acoustic phonon modes of the overall crystal framework. The analysis of mode-resolved spectral functions further show that vibrational modes with significant Cu contributions are strongly damped, corresponding to the breakdown of associated phonon quasiparticles. These results highlight the importance of strongly anharmonic effects in superionic systems, in which the traditional quasiharmonic phonon picture is insufficient, and pave the way toward combining machine-learning accelerated simulations with neutron scattering experiments to rationalize the complex atomic dynamics underlying ionic and thermal transport.

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

Large Enhancements in Optical and Piezoelectric Properties in Ferroelectric Zn 1- x Mg x O Thin Films through Engineering Electronic and Ionic Anharmonicities

Multifunctionality as a paradigm requires materials exhibiting multiple superior properties. Integrating second-order optical nonlinearity and large bandgap with piezoelectricity can, for example, enable broadband, strain-tunable photonics. Though very different phenomena at distinct frequencies, both second-order optical nonlinearity and piezoelectricity are third-rank polar tensors present only in acentric crystal structures. However, simultaneously enhancing both phenomena is highly challenging since it involves competing effects with tradeoffs. Recently, a large switchable ferroelectric polarization of ≈80 μC cm –2 was reported in Zn 1-x Mg x O films. Here, ferroelectric Zn 1-x Mg x O is demonstrated to be a platform that hosts simultaneously a 30% increase in the electronic bandgap, a 50% enhancement in the second harmonic generation (SHG) coefficients, and a near 200% improvement in the piezoelectric coefficients over pure ZnO. These enhancements are shown to be due to a 400% increase in the electronic anharmonicity and a ≈200% decrease in the ionic anharmonicity with Mg substitution. Precisely controllable periodic ferroelectric domain gratings are demonstrated down to 800 nm domain width, enabling ultraviolet quasi-phase-matched optical harmonic generation as well as domain-engineered piezoelectric devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Active Thermochemical Tables: the thermophysical and thermochemical properties of methyl, CH 3 , and methylene, CH 2 , corrected for nonrigid rotor and anharmonic oscillator effects

The thermophysical properties (isobaric heat capacity, entropy, enthalpy increment) of two prominent radicals, methyl, CH 3 , and methylene, CH 2 , were computed using the Nonrigid Rotor Anharmonic Oscillator (NRRAO) approach and compared to their RRHO counterparts, demonstrating significant differences between the results from the two approaches. Methylene presents a typical case in which the NRRAO thermophysical properties have significantly higher values than their RRHO counterparts at higher temperatures. In the case of methyl, the positive anharmonicity of the umbrella motion causes an opposite effect, and the NRRAO corrected thermophysical properties have lower values than their RRHO counterparts. The NRRAO corrected thermophysical properties, in turn, affect the resulting thermochemical properties. Here, two reactions important in combustion modelling were tested: the recombination of methyl radical with hydrogen atoms to form methane, and the recombination of two methyl radicals to form ethane. The related NRRAO equilibrium constants differ significantly from their RRHO analogs, and the consequences for chemical modelling are discussed. Also reported are the most current ATcT enthalpies of formation for CH n (n = 4-0) species and for C 2 H 6 , together with the tightly related sequential bond dissociation enthalpies along with the CHn series.

74 ATOMIC AND MOLECULAR PHYSICS↗
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 OSTI2022

Hybridization-driven strong anharmonicity in Yb-filled skutterudites

A high-pressure study of the structural and thermal transport properties is carried out on one of the most efficient filled skutterudites, $\mathrm{Yb_{0.3}Co_4Sb_{12}}$, to understand the relatively low thermal conductivity behavior in this family. By combining x-ray diffraction and Raman scattering measurements, we detect a phase transition at around 12.4 GPa. The mode Grüneisen parameters of the observed phonon modes are obtained from the determined bulk modulus and the phonon frequency shifts with pressure. The strong anharmonicity in this material is demonstrated by the obtained large average Grüneisen parameter. We also find the depressed group velocity within the low-frequency range, the flat guest mode avoided crossing with the acoustic-phonon mode, and the significant contribution of optical phonons. The hybridization of the guest atom and host lattice and the related enhanced anharmonicity are suggested to account for the low lattice thermal conductivity in the studied system. So, these findings provide new insight into how phonon-phonon interactions lower lattice thermal conductivity in this important thermoelectric family.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Toward 1% single-photon anharmonicity with periodically poled lithium niobate microring resonators

The absence of the single-photon nonlinearity has been a major roadblock in developing quantum photonic circuits at optical frequencies. In this paper, we demonstrate a periodically poled thin film lithium niobate microring resonator (PPLNMR) that reaches 5,000,000%/W second-harmonic conversion efficiency—almost 20-fold enhancement over the state-of-the-art—by accessing its largest χ <#comment/> ( 2 ) tensor component d 33 via quasi-phase matching. The corresponding single-photon coupling rate g / 2 π <#comment/> is estimated to be 1.2 MHz, which is an important milestone as it approaches the dissipation rate κ <#comment/> / 2 π <#comment/> of best-available lithium niobate microresonators developed in the community. Using a figure of merit defined as g / κ <#comment/> , our device reaches a single-photon nonlinear anharmonicity approaching 1%. We show that, by further scaling of the device, it is possible to improve the single-photon anharmonicity to a regime where photon blockade effect can be manifested.

Lu, Juanjuan (ORCID:0000000319389119)↗
DOE OSTI2024

Decoupling of the onset of anharmonicity between a protein and its surface water around 200 K

The protein dynamical transition at ~200 K, where the biomolecule transforms from a harmonic, non-functional form to an anharmonic, functional state, has been thought to be slaved to the thermal activation of dynamics in its surface hydration water. Here, by selectively probing the dynamics of protein and hydration water using elastic neutron scattering and isotopic labeling, we found that the onset of anharmonicity in the two components around 200 K is decoupled. The one in protein is an intrinsic transition, whose characteristic temperature is independent of the instrumental resolution time, but varies with the biomolecular structure and the amount of hydration, while the one of water is merely a resolution effect.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Calculation of mode lifetimes in weakly anharmonic solids using self-consistent ensemble eigenstates of the Liouvillian

Recent schemes for the efficient calculation of vibrational mode frequencies and lifetimes in weakly anharmonic systems such as solid crystal structures have shown promise, with simple relationships being observed between low order moments of the Liouvillian and the observed lifetime. While only shown empirically, two parameters seemed sufficient to determine the lifetime to high accuracy over a large range of temperatures for a few simple low degree of freedom models and a Lennard–Jones solid. In this paper we analytically reproduce the numerically observed relationship using a few generally applicable assumptions about the density of states. This advance is made by studying the autocorrelation based on an “ensemble eigenstate” of the Liouvillian. We propose the use of this ensemble eigenstate as the basis for autocorrelations to simplify calculations in richer systems. Finally, an analysis is also performed of the relationship between canonical and microcanonical autocorrelation functions and a long standing error in the long time behavior of such functions originally presented by Sen, Sinkovits, and Chakravarti is corrected.

36 MATERIALS SCIENCE↗
DOE OSTI2023

MSTor 2023 : A new version of the computer code for multistructural torsional anharmonicity, now with automatic torsional identification using redundant internal coordinates

The MSTor program is a computer program for calculating partition functions and thermodynamic functions of complex gas-phase molecules with multiple torsions; the multi-structural approximation with torsional anharmonicity (MS-T) can be used based on either a coupled torsional potential or an uncoupled torsional potential. The program can also carry out calculations in the multiple-structure local harmonic (MS-LH) approximation or multi-structural local quasiharmonic (MS-LQ) approximation or by the dual-level MS-T method. Furthermore, the program package includes eight utility codes that can be used as stand-alone programs to help the user to generate the input files for the code or to generate comparison results. The 2023 version of MSTor includes a new capability, namely automatic identification of torsional modes with redundant internal coordinates.

97 MATHEMATICS AND COMPUTING↗
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

Anharmonic phonons, thermal expansion, and nuclear quantum effects in Zn

The anharmonic behavior of phonons and thermal expansion of hexagonal zinc were studied from 15 to 690 K by inelastic neutron scattering (INS) and ab initio simulations. Phonon spectra were measured for Q-points covering the full Brillouin zone, giving the phonon density of states (DOS) and dispersions along high-symmetry directions. The dispersions were sharp at 15 K, but diffuse intensity was observed at energies above them. The dispersions broadened with temperature T and the diffuse intensity became stronger. This diffuse intensity appeared in all INS measurements and simulations, except for classical molecular dynamics at 15 K. The temperature-dependent effective potential (TDEP) method, which included the nuclear quantum effect from zero-point vibrational dynamics, was used to calculate the free energy and thermal expansion. For T < 100 K nuclear quantum effects were important for obtaining the correct negative thermal expansion, and path integral molecular dynamics (PIMD) was particularly effective for obtaining the negative thermal expansion in the basal plane. A Heisenberg-Langevin model for interacting phonons coupled to a thermal bath was able to reproduce the shape and intensity of the diffuse spectral features.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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

Phonon Anharmonicity and Spin–Phonon Coupling in CrI 3

We report on the far-infrared, temperature-dependent optical properties of a CrI 3 transition metal halide single crystal, a van der Waals ferromagnet (FM) with a Curie temperature of 61 K. In addition to the expected phonon modes determined by the crystalline symmetry, the optical reflectance and transmittance spectra of CrI 3 single crystals show many other excitations as a function of temperature as a consequence of the combination of a strong lattice anharmonicity and spin–phonon coupling. This complex vibrational spectrum highlights the presence of entangled interactions among the different degrees of freedom in CrI 3 .

36 MATERIALS SCIENCE↗