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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Scale-invariant machine-learning model accelerates the discovery of quaternary chalcogenides with ultralow lattice thermal conductivity

We design an advanced machine-learning (ML) model based on crystal graph convolutional neural network that is insensitive to volumes (i.e., scale) of the input crystal structures to discover novel quaternary chalcogenides, AMM'Q 3 (A/M/M' = alkali, alkaline earth, post-transition metals, lanthanides, and Q = chalcogens). These compounds are shown to possess ultralow lattice thermal conductivity (κ l ), a desired requirement for thermal-barrier coatings and thermoelectrics. Upon screening the thermodynamic stability of ~1 million compounds using the ML model iteratively and performing density-functional theory (DFT) calculations for a small fraction of compounds, we discover 99 compounds that are validated to be stable in DFT. Taking several DFT-stable compounds, we calculate their κ l using Peierls–Boltzmann transport equation, which reveals ultralow κ l (<2 Wm -1 K -1 at room temperature) due to their soft elasticity and strong phonon anharmonicity. Our work demonstrates the high efficiency of scale-invariant ML model in predicting novel compounds and presents experimental-research opportunities with these new compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Primitive to conventional geometry projection for efficient phonon transport calculations

Abstract The primitive Wigner-Seitz cell and corresponding first Brillouin zone (FBZ) are typically used in calculations of lattice vibrational and transport properties as they contain the smallest number of degrees of freedom and thus have the cheapest computational cost. However, in complex materials, the FBZ can take on irregular shapes where lattice symmetries are not apparent. Thus, conventional cells (with more atoms and regular shapes) are often used to describe materials, though dynamical and transport calculations are more expensive. Here we discuss an efficient anharmonic lattice dynamic method that maps conventional cell dynamics to primitive cell dynamics based on translational symmetries. Such symmetries have not been utilized in typical lattice dynamical calculations. This leads to phase interference conditions that act like conserved quantum numbers and a conservation rule for phonon scattering that is hidden in conventional dynamics which significantly reduces the computational cost. We demonstrate this method for phonon transport in a variety of materials with inputs from first-principles calculations and attribute its efficiency to reduced scattering phase space and fewer summations in scattering matrix element calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A high-throughput framework for lattice dynamics

We develop an automated high-throughput workflow for calculating lattice dynamical properties from first principles including those dictated by anharmonicity. The pipeline automatically computes interatomic force constants (IFCs) up to 4th order from perturbed training supercells, and uses the IFCs to calculate lattice thermal conductivity, coefficient of thermal expansion, and vibrational free energy and entropy. It performs phonon renormalization for dynamically unstable compounds to obtain real effective phonon spectra at finite temperatures and calculates the associated free energy corrections. The methods and parameters are chosen to balance computational efficiency and result accuracy, assessed through convergence testing and comparisons with experimental measurements. Deployment of this workflow at a large scale would facilitate materials discovery efforts toward functionalities including thermoelectrics, contact materials, ferroelectrics, aerospace components, as well as general phase diagram construction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Dynamic crystallography reveals spontaneous anisotropy in cubic GeTe

Cubic energy materials such as thermoelectrics or hybrid perovskite materials are often understood to be highly disordered. In GeTe and related IV–VI compounds, this is thought to provide the low thermal conductivities needed for thermoelectric applications. Since conventional crystallography cannot distinguish between static disorder and atomic motions, we develop the energy-resolved variable-shutter pair distribution function technique. This collects structural snapshots with varying exposure times, on timescales relevant for atomic motions. In disagreement with previous interpretations, we find the time-averaged structure of GeTe to be crystalline at all temperatures, but with anisotropic anharmonic dynamics at higher temperatures that resemble static disorder at fast shutter speeds, with correlated ferroelectric fluctuations along the <100> c direction. We show that this anisotropy naturally emerges from a Ginzburg–Landau model that couples polarization fluctuations through long-range elastic interactions. By accessing time-dependent atomic correlations in energy materials, we resolve the long-standing disagreement between local and average structure probes and show that spontaneous anisotropy is ubiquitous in cubic IV–VI materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Nonlinear nanomechanical resonators approaching the quantum ground state

It is an open question whether mechanical resonators can be made nonlinear with vibrations approaching the quantum ground state. This requires the engineering of a mechanical nonlinearity far beyond what has been realized so far. Here we discover a mechanism to boost the Duffing nonlinearity by coupling the vibrations of a nanotube resonator to single-electron tunnelling and by operating the system in the ultrastrong-coupling regime. We find that thermal vibrations become highly nonlinear when lowering the temperature. The average vibration amplitude at the lowest temperature is 13 times the zero-point motion, with approximately 42% of the thermal energy stored in the anharmonic part of the potential. Our work may enable the realization of mechanical Schrödinger cat states, mechanical qubits and quantum simulators emulating the electron–phonon coupling.

36 MATERIALS SCIENCE↗
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 OSTI2022

Molecular dynamics on quantum annealers

Abstract In this work we demonstrate a practical prospect of using quantum annealers for simulation of molecular dynamics. A methodology developed for this goal, dubbed Quantum Differential Equations (QDE), is applied to propagate classical trajectories for the vibration of the hydrogen molecule in several regimes: nearly harmonic, highly anharmonic, and dissociative motion. The results obtained using the D-Wave 2000Q quantum annealer are all consistent and quickly converge to the analytical reference solution. Several alternative strategies for such calculations are explored and it was found that the most accurate results and the best efficiency are obtained by combining the quantum annealer with classical post-processing (greedy algorithm). Importantly, the QDE framework developed here is entirely general and can be applied to solve any system of first-order ordinary nonlinear differential equations using a quantum annealer.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Ultrafast solvent-to-solute proton transfer mediated by intermolecular coherent vibrations

Ultrafast photoinduced excited-state proton transfer (ESPT) plays a crucial role in protecting biomolecules and functional materials from photodamage. However, the influence of solute-solvent interactions on ESPT dynamics remains under active investigation. Here, we present an ultrafast spectroscopic study of ESPT in the photobase 2-(2´-pyridyl)benzimidazole (PBI) in methanol. Ultrafast absorption spectroscopy, supported by quantum chemical calculations, reveals three distinct kinetic steps: (1) a 2.2 ps solvent-to-solute proton transfer, (2) subsequent nonradiative relaxation to the ground state within 31 ps, producing a vibrationally hot ensemble with substantial excess kinetic energy, and (3) equilibration as this energy dissipates into the surrounding solvent bath over 186 ps. Femtosecond-resolved dynamics exhibit oscillatory signals indicative of coherent wavepacket motion on the S 1 potential energy surface. A phase flip in the excited-state absorption maximum confirms this assignment. Fourier analysis resolves two dominant periods (∼117 fs and ∼340 fs), corresponding to in-plane and out-of-plane vibrational modes coupled between PBI and the hydrogen-bonded methanol molecule. The rapid dephasing ( < 300 fs) suggests that the nuclear wavefunction evolves on an anharmonic potential energy surface while traversing the ESPT reaction coordinate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2022

Tracking carrier and exciton dynamics in mixed-cation lead mixed-halide perovskite thin films

Mixed-cation lead mixed-halide perovskites simultaneously possess structural stability and high power conversion efficiency. A thorough study of both carrier and exciton dynamics is needed to understand the photophysical properties that underpin its superior photovoltaic performance. By utilizing a broadband transient absorption spectroscopy, we observe the carrier and exciton dynamics in a FA 0.85 Cs 0.15 Pb(I 0.97 Br 0.03 ) 3 (FCPIB) perovskite by simultaneously resolving the carrier and exciton contribution to the transient change of the absorption spectra, from which the carrier density and exciton oscillator strength can be determined. Our data reveal a quick and significant conversion of the photogenerated carriers to excitons, on top of the usual carrier recombination process. Moreover, the decay of carrier density shows a change of kinetics from a second-order recombination at high pump fluence to a third-order recombination at low pump fluence. Our analysis utilizes band anharmonicity, presents an independent determination of electronic temperature and quasi-Fermi energy, and reveals an interesting interplay among the processes of carrier cooling, exciton formation/decay and carrier recombination, all as a function of time after photoexcitation. Our work demonstrates the use of pump fluence as a knob to tune the relative populations of carriers and excitons in halide perovskite materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mapping the nonequilibrium order parameter of a quasi-two dimensional charge density wave system

The driving force of a charge density wave (CDW) transition in quasi-two dimensional systems is still debated, while being crucial in understanding electronic correlation in such materials. Here we use femtosecond time- and angle-resolved photoemission spectroscopy combined with computational methods to investigate the coherent lattice dynamics of a prototypical CDW system. The photo-induced temporal evolution of the periodic lattice distortion associated with the amplitude mode reveals the dynamics of the free energy functional governing the order parameter. Our approach establishes that optically-induced screening rather than CDW melting at the electronic level leads to a transiently modified potential which explains the anharmonic behaviour of the amplitude mode and discloses the structural origin of the symmetry-breaking phase transition.

electronic properties and materials↗
DOE OSTI2020

Stability and metallization of solid oxygen at high pressure

The phase diagram of oxygen is investigated for pressures from 50 to 130 GPa and temperatures up to 1200 K using first-principles theory. A metallic molecular structure with the P6 3 /mmc symmetry (η' phase) is determined to be thermodynamically stable in this pressure range at elevated temperatures above the ε(O 8 ) phase. Crucial for obtaining this result is the inclusion of anharmonic lattice dynamics effects and accurate calculations of exchange interactions in the presence of thermal disorder. Overall, we present analysis of electronic, structural, and thermodynamic properties of solid oxygen at 0 K and finite temperature with hybrid exchange functionals, including a comparison with available experimental data.

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

Challenges for density functional theory: calculation of CO adsorption on electrocatalytically relevant metals

Density Functional Theory (DFT) is currently the most tractable choice of theoretical model used to understand the mechanistic pathways for electrocatalytic processes such as CO 2 or CO reduction. In this work, we assess the performance of two DFT functionals designed specifically to describe surface interactions, RTPSS and RPBE, as well as two popular meta-GGA functionals, SCAN and B97M-rV, that have not been a priori optimized for better interfacial properties. We assess all four functionals against available experimental data for prediction of bulk and bare surface properties on four electrocatalytically relevant metals, Au, Ag, Cu, and Pt, and for binding CO to surfaces of these metals. To partially mitigate issues such as thermal and anharmonic corrections associated with comparing computations with experiments, molecular benchmarks against high level quantum chemistry are reported for CO complexes with Au, Ag, Cu and Pt atoms, as well as the CO–water complex and the water dimer. Overall, we find that the surface modified RPBE functional performs reliably for many of the benchmarks examined here, and the meta-GGA functionals also show promising results. Specifically B97M-rV predicts the correct site preference for CO binding on Ag and Au (the only functional tested here to do so), while RTPSS performs well for surface relaxations and binding of CO on Pt and Cu.

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

A theoretical investigation of the hydrolysis of uranium hexafluoride: the initiation mechanism and vibrational spectroscopy

Depleted uranium hexafluoride (UF 6 ), a stockpiled byproduct of the nuclear fuel cycle, reacts readily with atmospheric humidity, but the mechanism is poorly understood. Here we compare several potential initiation steps at a consistent level of theory, generating underlying structures and vibrational modes using hybrid density functional theory (DFT) and computing relative energies of stationary points with double-hybrid (DH) DFT. A benchmark comparison is performed to assess the quality of DH-DFT data using reference energy differences obtained using a complete-basis-limit coupled-cluster (CC) composite method. The associated large-basis CC computations were enabled by a new general-purpose pseudopotential capability implemented as part of this work. Dispersion-corrected parameter-free DH-DFT methods, namely PBE0-DH-D3(BJ) and PBE-QIDH-D3(BJ), provided mean unsigned errors within chemical accuracy (1 kcal mol -1 ) for a set of barrier heights corresponding to the most energetically favorable initiation steps. The hydrolysis mechanism is found to proceed via intermolecular hydrogen transfer within van der Waals complexes involving UF 6 , UF 5 OH, and UOF 4 , in agreement with previous studies, followed by the formation of a previously unappreciated dihydroxide intermediate, UF 4 (OH) 2 . The dihydroxide is predicted to form under both kinetic and thermodynamic control, and, unlike the alternate pathway leading to the UO 2 F 2 monomer, its reaction energy is exothermic, in agreement with observation. Finally, harmonic and anharmonic vibrational simulations are performed to reinterpret literature infrared spectroscopy in light of this newly identified species.

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

Ab initio molecular dynamics with enhanced sampling in heterogeneous catalysis

Ab initio molecular dynamics simulations combined with enhanced sampling techniques are becoming widespread methods to investigate chemical phenomena in catalytic systems. These techniques automatically include finite temperature effects, anharmonicity, and collective dynamics in their robust description of enthalpic and entropic contributions, which can have significant impact on reaction free energy landscapes. This contrasts with standard ab initio static approaches that are based on assessing reaction free energies from various coarse-grained descriptions of the reaction potential energy surface. Enhanced sampling ab initio molecular dynamics opens the way to first principles simulations of systems of increasing complexity like solid/liquid catalytic interfaces. Here, we aim at guiding the reader through the basis of these techniques, summarizing their fundamental theoretical and practical aspects, and reviewing the relevant literature in the field. After a brief introduction to the problem, we will illustrate the advantage of using molecular simulations to include finite temperature effects, examine the most common ab initio techniques currently in use, describe their application to solid state heterogeneous catalysts, and finally critically review the most popular enhanced sampling techniques used in computational catalysis.

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

Synthesis, structure, electronic and thermal properties of sphalerite CuZn 2 InS 4

Quaternary chalcogenides continue to be of interest due to the variety of physical properties they possess, as well as their potential for different applications of interest. Investigations on materials with the sphalerite crystal structure have only recently begun. In this study we have synthesized sulfur-based sphalerite quaternary chalcogenides, including off-stoichiometric compositions, and investigated the temperature-dependent electronic, thermal and structural properties of these materials. Insulating to semiconducting transport is observed with stoichiometric variation, and analyses of heat capacity and thermal expansion revealed lattice anharmonicity that contributes to the low thermal conductivity these materials possess. Finally, we include similar analyses for CuZn 2 InSe 4 and compare these sphalerite quaternary chalcogenides to that of zinc blende binaries in order to fully understand the origin of the low thermal conductivity these quaternary chalcogenides possess.

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

Structural and thermal properties of ultralow thermal conductivity Ba 3 Cu 2 Sn 3 Se 10

The thermal properties of Ba 3 Cu 2 Sn 3 Se 10 were investigated by measurement of the thermal conductivity and heat capacity. The chemical bonding in this diamagnetic material was investigated using structural data from Rietveld refinement and calculated electron localization. This quaternary chalcogenide is monoclinic (P2 1 /c), has a large unit cell with 72 atoms in the primitive cell, and a high local coordination environment. Additionally, the Debye temperature (162 K) and average speed of sound (1666 m s -1 ) are relatively low with a very small electronic contribution to the heat capacity. Ultralow thermal conductivity (0.46 W m -1 K -1 at room temperature) is attributed to the relatively weak chemical bonding and intrinsic anharmonicity, in addition to a large unit cell. This work is part of the continuing effort to explore quaternary chalcogenides with intrinsically low thermal conductivity and identify the features that result in a low thermal conductivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗