Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phonon density of states”

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 OSTI2025

Generalized approach for rapid entropy calculation of liquids and solids

We build a comprehensive methodology for the fast computation of entropy across both solid and liquid phases. The proposed method utilizes a single trajectory of molecular dynamics (MD) to facilitate the calculation of entropy, which is composed of three components. The electronic entropy is determined through the temporal average acquired from density functional theory MD simulations. The vibrational entropy, typically the predominant contributor to the total entropy, even within the liquid state, is evaluated by computing the phonon density of states via the velocity autocorrelation function. The most arduous component to quantify, the configurational entropy, is assessed by probability analysis of the local structural arrangement and atomic distribution. We illustrate, through a variety of examples, that this method is both a versatile and valid technique for characterizing the thermodynamic states of both solids and liquids. Furthermore, this method is employed to expedite the calculation of melting temperatures, demonstrating its practical utility in computational thermodynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermodynamic origin of the pressure-induced Invar effect: General criterion and experimental study of Fe 68 ⁢Pd 32

Synchrotron X-ray diffraction measurements were performed on 57 Fe 68 Pd 32 at multiple pressures and two temperatures in a diamond-anvil cell. Between 4 and 11GPa, the thermal expansion was zero or slightly negative. This pressure-induced Invar effect was studied further with 57 Fe nuclear forward scattering and nuclear resonant inelastic X-ray scattering to obtain information on the pressure-induced changes in both the magnetization and the phonon density of states. Magnetic entropy and phonon entropy were obtained from these results, the latter with additional measurements from inelastic neutron scattering to account for the contributions from Pd atoms. The dependencies of these entropies on pressure gave the magnetic and phonon contributions to thermal expansion. These canceled in the region of pressure-induced Invar behavior, even though they individually increased by more than a factor of 2 below the Curie pressure. The behavior of phonons gives evidence for spin-phonon interactions in Fe 68 Pd 32 . A general explanation of the pressure-induced Invar effect is presented, showing that Invar behavior is typically expected at a pressure P* below a magnetic transition at pressure PC. The difference in pressure P C − P* scales with the fractional reduction in magnetic exchange interaction divided by the average Gr¨uneisen parameter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A noncollinear density functional theory ansatz for the phononic and thermodynamic properties of α -Pu

Plutonium's phase diagram is host to complex structures and interactions that make the description of its ground state properties elusive. Using all-electron density functional theory calculations, we study the thermodynamic properties of α-Pu. To do this, we build on recent work in the literature by introducing a novel noncollinear magnetic ansatz for α-Pu's ground state. We show that the noncollinear ansatz improves the description of the experimental phonon density of states, heat capacity, and thermal expansion over previous results obtained with a collinear ansatz. Here, despite lacking a bond-equivalent picture of the bonding, these new results for α-Pu, along with recent results on δ-Pu, demonstrate the efficacy of noncollinear ansatzes for the description of plutonium.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermal Neutron Scattering Law for Beryllium Hydride and Critical Mass Calculations

The thermal neutron scattering law (TSL) for crystalline beryllium hydride (BeH 2 ) is developed from first-principles ab initio lattice dynamics calculations and the impact of neutron thermalization in this material on critical mass is estimated. BeH 2 has a body-centered orthorhombic crystal structure with 12 molecules per unit cell and a theoretical density of 0.755 g/cm 3 . The vibrational (phonon) densities of states for H and Be bound in BeH 2 are determined using VASP density functional theory and PHONON lattice dynamics calculations. The TSLs for H bound in BeH 2 , H(BeH 2 ), and Be bound in BeH 2 , Be(BeH 2 ), are then evaluated in the incoherent approximation from the calculated H and Be partial phonon density of states using FLASSH. Finally, critical mass as a function of 235 U loading density for bare and reflected BeH 2 moderated spheres is predicted from MC21 Monte Carlo neutron transport calculations using ENDF/B-VIII.0 cross sections and the H(BeH 2 ) and Be(BeH 2 ) TSL evaluations. Comparisons are made to water (H 2 O), polyethylene (CH 2 ), and beryllium oxide (BeO) as moderators. These critical mass predictions are a refinement upon the prior work by Rao and Srinivasan that neglected thermal neutron scattering effects. The minimum critical mass of a BeH 2 moderated assembly is estimated to be 0.207 kg 235 U for a 0.20 m thick BeO reflected sphere and 0.178 kg 235 U for a 0.40 m thick BeO reflected sphere.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

The influence of phonon softening on the superconducting critical temperature of Sn nanostructures

Abstract The increase in superconducting transition temperature (T C ) of Sn nanostructures in comparison to bulk, was studied. Changes in the phonon density of states (PDOS) of the weakly coupled superconductor Sn were analyzed and correlated with the increase in T C measured by magnetometry. The PDOS of all nanostructured samples shows a slightly increased number of low-energy phonon modes and a strong decrease in the number of high-energy phonon modes in comparison to the bulk Sn PDOS. The phonon densities of states, which were determined previously using nuclear resonant inelastic X-ray scattering, were used to calculate the superconducting transition temperature using the Allen-Dynes-McMillan (ADMM) formalism. Both the calculated as well as the experimentally determined values of T C show an increase compared to the bulk superconducting transition temperature. The good agreement between these values indicates that phonon softening has a major influence on the superconducting transition temperature of Sn nanostructures. The influence of electron confinement effects appears to be minor in these systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Machine learned potential for high-throughput phonon calculations of metal—organic frameworks

Metal–organic frameworks (MOFs) are highly porous and versatile materials studied extensively for applications such as carbon capture and water harvesting. However, computing phonon-mediated properties in MOFs, like thermal expansion and mechanical stability, remains challenging due to the large number of atoms per unit cell, making traditional Density Functional Theory (DFT) methods impractical for high-throughput screening. Recent advances in machine learning potentials have led to foundation atomistic models, such as MACE-MP-0, that accurately predict equilibrium structures but struggle with phonon properties of MOFs. In this work, we developed a workflow for computing phonons in MOFs within the quasi-harmonic approximation with a fine-tuned MACE model, MACE-MP-MOF0. The model was trained on a curated dataset of 127 representative and diverse MOFs. The fine-tuned MACE-MP-MOF0 improves the accuracy of phonon density of states and corrects the imaginary phonon modes of MACE-MP-0, enabling high-throughput phonon calculations with state-of-the-art precision. The model successfully predicts thermal expansion and bulk moduli in agreement with DFT and experimental data for several well-known MOFs. These results highlight the potential of MACE-MP-MOF0 in guiding MOF design for applications in energy storage and thermoelectrics.

Elena, Alin Marin↗
DOE OSTI2021

Vibrational properties and thermal transport in quaternary chalcogenides: The case of Te-based compositions

Vibrational thermal properties of CuZn 2 InTe 4 , AgZn 2 InTe 4 , and Cu 2 CdSnTe 4 , derived from binary II-VI zinc-blendes, are reported based on first-principles calculations. While the chalcogenide atoms in these materials have the same lattice positions, the cation atom arrangements vary, resulting in different crystal symmetries and subsequent properties. The compositional differences have important effects on the vibrational thermal characteristics of the studied materials, which demonstrate that low-frequency optical phonons hybridize with acoustic phonons and lead to enhanced phonon-phonon scattering and low lattice thermal conductivities. The phonon density of states, mode Grüneisen parameters, and phonon scattering rates are also calculated, enabling deeper insight into the microscopic thermal conduction processes in these materials. Compositional variations drive differences among the three materials considered here; nonetheless, their structural similarities and generally low thermal conductivities (0.5–4 W/m K at room temperature) suggest that other similar II-VI zinc-blende derived materials will also exhibit similarly low values, as also corroborated by experimental data. Finally, this, combined with the versatility in designing a variety of motifs on the overall structure, makes quaternary chalcogenides interesting for thermal management and energy conversion applications that require low thermal conductivity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phonons in Copper Diphosphide (CuP 2 ): Raman Spectroscopy and Lattice Dynamics Calculations

Copper diphosphide (CuP 2 ) is an emerging binary semiconductor with promising properties for energy conversion and storage applications. While functionality and possible applications of CuP 2 have been studied, there is a curious gap in the investigation of its vibrational properties. In this work, we provide a reference Raman spectrum of CuP 2 , with a complete analysis of all Raman active modes from both experimental and theoretical perspectives. Raman measurements have been performed on polycrystalline CuP 2 thin films with close to stoichiometric composition. Detailed deconvolution of the Raman spectrum with Lorentzian curves has allowed identification of all theoretically predicted Raman active modes (9A g and 9B g ), including their positions and symmetry assignment. Furthermore, calculations of the phonon density of states (PDOS), as well as the phonon dispersions, provide a microscopic understanding of the experimentally observed phonon lines, in addition to the assignment to the specific lattice eigenmodes. We further provide the theoretically predicted positions of the infrared (IR) active modes, along with the simulated IR spectrum from density functional theory (DFT). Overall good agreement is found between the experimental and DFT-calculated Raman spectra of CuP 2 , providing a reference platform for future investigations on this material.

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

Interplay between phonon downconversion efficiency, density of states at Fermi energy, and intrinsic energy resolution for microwave kinectic inductance detectors

Microwave Kinetic Inductance detectors (MKIDs) have been recognized as a powerful new tool for single photon detection. These highly multiplexed superconducting devices give timing and energy measurement for every detected photon. However, the full potential of MKID single photon spectroscopy has not been reached , the achieved energy resolution is lower than expected from first principles. Here, we study the efficiency in the phonon downconversion process following the absorption of energetic photons in MKIDs. Assuming previously published material properties, we measure an average downconversion efficiency for three TiN resonators is $\eta$=0.14. We discuss how this efficiency can impact the intrinsic energy resolution of MKID, and how any uncertainty in the unknown density of electron states at the Fermi energy directly affects the efficiency estimations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Methodology for Generating Covariance Data of Thermal Neutron Scattering Cross Sections

This paper details and implements a framework for evaluating thermal neutron scattering cross sections that provide S(α,β) data and covariance data for hydrogen in light water. This methodology involves perturbing model parameters of molecular dynamics potentials and fitting the simulation results to experimental data. The framework is general and can be applied to any material or simulation method. The fit is made using the Unified Monte Carlo method to experimentally measure double-differential scattering cross sections of light water at the Spallation Neutron Source at Oak Ridge National Laboratory. Mean values and covariance data were generated for model parameters, phonon density of states, double-differential cross sections, and total scattering cross sections. These posterior parameter values were very similar to their prior values with a maximum relative error of 0.54%. This falls within in the Unified Monte Carlo–calculated uncertainties on the order of 2.7%. Additionally, posterior double-differential cross sections agree favorably with ENDF/B-VIII.0 cross sections. The new thermal scattering law was tested by comparing it against benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook, which showed a slight improvement over the ENDF/B-VIII.0 library. Additionally, the covariance matrix of the phonon density of states was validated to confirm that the spread of k eff from the density of states used to generate the covariance matrix was similar to the spread of k eff from the density of states of the sampled covariance matrix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2020

Desorption Kinetics of O and CO from Graphitic Carbon Surfaces

The desorption of O/CO from graphitic carbon surfaces is investigated using a one-dimensional model describing the adsorbate interactions with the surface phonon bath. The kinetics of desorption are described through the solution of a master equation for the time-dependent population of the adsorbate in an oscillator state, which is modified through thermal fluctuations at the surface. The interaction of the adsorbate with the surface phonons is explicitly captured by using the computed phonon density Of states (PDOS) of the surface. The coupling of the adsorbate with the phonon bath results in the transition of the adsorbate up and down a vibrational ladder. The adsorbate-surface interaction is represented in the model using a Morse potential, which allows for the desorption process to be directly modeled as a transition from bound to free (continuum) state. The PDOS is a property of the material and the lattice; and is highly sensitive to the presence of defects. The effect of etch pits along with random surface defects on the PDOS is considered in the present work. The presence of defects causes a redshift and broadening of the PDOS, which in turn changes the phonon frequency modes available for adsorbate coupling at the surface. Using the realistic PDOS distributions, the phonon-induced desorption (PID) model was used to compute the transition and desorption rates for both pristine and defective systems. Mathissen’s rule is used to compute the phonon relaxation time for pristine and defective systems based on the phonon scattering times for each of the different scattering processes. First, the desorption rates of the pristine system is fitted against the experimental values to obtain the Morse potential parameters for each of the observed adatoms. These Morse potential parameters are used along with the defective PDOS and phonon relaxation time to compute the desorption rates for the defective system. The defective system rates (both transition and desorption) were consistently lower in comparison with the pristine system. The difference between the transition rates is more significant at lower initial states due to higher energy spacing between the levels. In the case of the desorption rates, the difference between the defective and pristine system is more significant at higher temperatures. The desorption rates for each of the system shows an order of magnitude decrease with the strongly bound systems exhibiting the greatest reduction in the desorption rates.

Swaminathan-Gopalan, Krishnan↗
DOE OSTI2023

Vibronic Relaxation Pathways in Molecular Spin Qubit Na 9 [Ho(W 5 O 18 ) 2 ]·35H 2 O under Pressure

In order to explore how spectral sparsity and vibronic decoherence pathways can be controlled in a model qubit system with atomic clock transitions, we combined diamond anvil cell techniques with synchrotron-based far infrared spectroscopy and first-principles calculations to reveal the vibrational response of Na 9 [Ho(W 5 O 18 ) 2 ]·35H 2 O under compression. Because the hole in the phonon density of states acts to reduce the overlap between the phonons and f manifold excitations in this system, we postulated that pressure might move the HoO 4 rocking, bending, and asymmetric stretching modes that couple with the M J = ±5, ±2, and ±7 levels out of resonance, reducing their interactions and minimizing decoherence processes, while a potentially beneficial strategy for some molecular qubits, pressure slightly hardens the phonons in Na 9 [Ho(W 5 O 18 ) 2 ]·35H 2 O and systematically fills in the transparency window in the phonon response. The net result is that the vibrational spectrum becomes less sparse and the overlap with the various M J levels of the Ho 3+ ion actually increases. These findings suggest that negative pressure, achieved using chemical means or elongational strain, could further open the transparency window in this rare earth-containing spin qubit system, thus paving the way for the use of device surfaces and interface elongational/compressive strains to better manage decoherence pathways.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Influence of hydrogenation on the vibrational density of states of magnetocaloric LaFe 11.4 Si 1.6 H 1.6

We report on the impact of magnetoelastic coupling on the magnetocaloric properties of LaFe 11.4 Si 1.6 H 1.6 in terms of the vibrational (phonon) density of states (VDOS), which we determined with 57 Fe nuclear resonant inelastic x-ray scattering (NRIXS) measurements and with density functional theory (DFT) based first-principles calculations in the ferromagnetic (FM) low-temperature and paramagnetic (PM) high-temperature phase. In experiments and calculations, we observe pronounced differences in the shape of the Fe-partial VDOS between nonhydrogenated and hydrogenated samples. Furthermore, this shows that hydrogen not only shifts the temperature of the first-order phase transition, but also affects the elastic response of the Fe subsystem significantly. In turn, the anomalous redshift of the Fe VDOS, observed by going to the low-volume PM phase, survives hydrogenation. As a consequence, the change in the Fe-specific vibrational entropy ΔS lat across the phase transition has the same sign as the magnetic and electronic contribution. DFT calculations show that the same mechanism, which is a consequence of the itinerant electron metamagnetism associated with the Fe subsystem, is effective in both the hydrogenated and the hydrogen-free compounds. Although reduced by 50% as compared to the hydrogen-free system, the measured change ΔS lat of 3.2 ± 1.9)J/ kg K across the FM-to-PM transition contributes with ~35 % significantly and cooperatively to the total isothermal entropy change ΔS iso . Hydrogenation is observed to induce an overall blueshift of the Fe VDOS with respect to the H-free compound; this effect, together with the enhanced Debye temperature observed, is a fingerprint of the hardening of the Fe sublattice by hydrogen incorporation. In addition, the mean Debye velocity of sound of LaFe 11.4 Si 1.6 H 1.6 was determined from the NRIXS and the DFT data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Probing charge density wave phases and the Mott transition in 1 T -TaS 2 by inelastic light scattering

We present a polarization-resolved, high-resolution Raman scattering study of the three consecutive charge density wave (CDW) regimes in 1T - TaS 2 single crystals, supported by ab initio calculations. Our analysis of the spectra within the low-temperature commensurate (C-CDW) regime shows P$\bar3$ symmetry of the system, thus excluding the previously proposed triclinic stacking of the “star-of-David” structure, and promoting trigonal or hexagonal stacking instead. The spectra of the high-temperature incommensurate (IC-CDW) phase directly project the phonon density of states due to the breaking of the translational invariance, supplemented by sizable electron-phonon coupling. Between 200 and 352 K, our Raman spectra show contributions from both the IC-CDW and the C-CDW phases, indicating their coexistence in the so-called nearly commensurate (NC-CDW) phase. The temperature dependence of the symmetry-resolved Raman conductivity indicates the stepwise reduction of the density of states in the CDW phases, followed by a Mott transition within the C-CDW phase. We determine the size of the Mott gap to be Ω gap ≈ 170 – 190 meV, and track its temperature dependence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Controlling interactions between high-frequency phonons and single quantum systems using phononic crystals

The ability to control phonons in solids is key in many fields of quantum science, ranging from quantum information processing to sensing. Phonons often act as a source of noise and decoherence when solid-state quantum systems interact with the phonon bath of their host matrix. In this study, we demonstrate the ability to control the phononic local density of states of the host matrix using phononic crystals and measure its positive impact on single quantum systems. We design and fabricate diamond phononic crystals with features down to around 20 nm, resulting in a high-frequency complete phononic bandgap from 50 to 70 GHz. The engineered local density of states is probed using single silicon-vacancy colour centres embedded in the phononic crystals. We observe an 18-fold reduction in the phonon-induced orbital relaxation rate of the emitters compared to bulk, thereby demonstrating that the phononic crystal suppresses spontaneous single-phonon processes. Furthermore, we show that our approach can efficiently suppress single-phonon-emitter interactions up to 20 K, allowing the investigation of multi-phonon processes in the emitters. Our results represent an important step towards the realization of efficient phonon-emitter interfaces that can be used for quantum acoustodynamics and quantum phononic networks.

quantum information↗
DOE OSTI2021

Phonon confinement and interface lattice dynamics of ultrathin high- k rare earth sesquioxide films: the case of Eu 2 O 3 on YSZ(001)

The spatial confinement of atoms at surfaces and interfaces significantly alters the lattice dynamics of thin films, heterostructures and multilayers. Ultrathin films with high dielectric constants (high-k) are of paramount interest for applications as gate layers in current and future integrated circuits. Here we report a lattice dynamics study of high-k Eu 2 O 3 films with thicknesses of 21.3, 2.2, 1.3, and 0.8 nm deposited on YSZ(001). The Eu-partial phonon density of states (PDOS), obtained from nuclear inelastic scattering, exhibits broadening of the phonon peaks accompanied by up to a four-fold enhancement of the number of low-energy states compared to the ab initio calculated PDOS of a perfect Eu 2 O 3 crystal. Our analysis demonstrates that while the former effect reflects the reduced phonon lifetimes observed in thin films due to scattering from lattice defects, the latter phenomenon arises from an ultrathin EuO layer formed between the thin Eu 2 O 3 film and the YSZ(001) substrate. Thus, our work uncovers another potential source of vibrational anomalies in thin films and multilayers, which has to be cautiously considered.

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

Strong Proton‐Phonon Coupling Drives Fast Ion Transport in Perovskites

Conduction of protons in solids is a cooperative process propelled by phonons, with molecular details obscured by the irregular movements in the thermal bath. It is shown that substitution with Y forms an imaginary phonon mode, instrumental for the function as proton conductor and effectively lowering the activation barrier for proton transport. To untangle the interplay in the exemplary proton conductor BaSn 0.9 Y 0.1 O 3 , its crystallographic structure is determined with high resolution neutron diffractometry and its phonon density of states with density functional theory calculations, experimentally validated by element specific nuclear resonant vibration spectroscopy. Based on phonon analysis, a quantitative transport model is present, which predicts the activation energy and performance by the ratio of ionic radii. Rather than individual vibrational modes, it is the oxygen sub-lattice which exerts its momentum on the protons. The extent of this momentum transfer is governed by the ratio of ionic radii. This model extends the transition state theory by the phonon-phonon interaction and complements the previously proposed idea that lattice dynamics is decisive for proton transport and specifies which properties of the material exactly define the vibration properties.

activation barriers↗
DOE OSTI2021

DFTTK: Density Functional Theory ToolKit for high-throughput lattice dynamics calculations

In this work, we present a software package in Python for high-throughput first-principles calculations of thermodynamic properties at finite temperatures, which we refer to as DFTTK (Density Functional Theory ToolKit). DFTTK is based on the atomate package and integrates our experiences in the last decades on the development of theoretical methods and computational softwares. It includes task submissions on all major operating systems and task executions on high-performance computing environments. Furthermore, the distribution of the DFTTK package comes with examples of calculations of phonon density of states, heat capacity, entropy, enthalpy, and free energy under the quasi-harmonic phonon scheme for the stoichiometric phases of Al, Ni, Al 3 Ni, AlNi, AlNi 3 , Al 3 Ni 4 , and Al 3 Ni 5 , and the fcc solution phases treated using the special quasirandom structures at the compositions of Al 3 Ni, AlNi, and AlNi 3 .

97 MATHEMATICS AND COMPUTING↗