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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Exploring the Vibrational Side of Spin-Phonon Coupling in Single-Molecule Magnets via 161 Dy Nuclear Resonance Vibrational Spectroscopy

Synchrotron-based nuclear resonance vibrational spectroscopy (NRVS) using the Mossbauer isotope 161 Dy has been employed for the first time to study the vibrational properties of a single-molecule magnet (SMM) incorporating Dy III , namely [Dy(Cy3PO)(2)(H2O)(5)]Br-3•2 (Cy3PO)•2 H2O•2 EtOH. The experimental partial phonon density of states (pDOS), which includes all vibrational modes involving a displacement of the Dy III ion, was reproduced by means of simulations using density functional theory (DFT), enabling the assignment of all intramolecular vibrational modes. This study proves that 161 Dy NRVS is a powerful experimental tool with significant potential to help to clarify the role of phonons in SMMs.

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

Lattice Dynamics and Structural Phase Transitions in Eu 2 O 3

Using the density functional theory, we study the structural and lattice dynamical properties of europium sesquioxide (Eu 2 O 3 ) in the cubic, trigonal, and monoclinic phases. The obtained lattice parameters and energies of the Raman modes show a good agreement with the available experimental data. The Eu-partial phonon density of states calculated for the cubic structure is compared with the nuclear inelastic scattering data obtained from a 20 nm thick Eu 2 O 3 film deposited on a YSZ substrate. A small shift of the experimental spectrum to higher energies results from a compressive strain induced by the substrate. On the basis of lattice and phonon properties, we analyze the mechanisms of structural transitions between different phases of Eu 2 O 3 .

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

Investigating the non-classical M-H 2 bonding in OsClH 3 (PPh 3 ) 3

Dihydrogen ligation to metal centers results in a wide range of metal-hydrogen interactions. Compounds containing non-classical dihydrogen ligands in which the H-H contacts are significantly longer than the 0.74 A distance found in solid H 2 are classified as Kubas-like compounds, stretched dihydrogen ligands, or compressed dihydrides, depending on the H-H contact length. Compounds containing stretched dihydrogen and compressed dihydride ligands are somewhat rare. Developing a more complete understanding of the metal-hydrogen and hydrogen-hydrogen interactions in these species may provide insights into the mechanisms for the oxidative addition reaction and hydrogen storage. Here, we use diffraction and inelastic neutron scattering (INS) measurements, paired with both molecular and solid-state density functional theory (DFT), to fully characterize the nature of the metal-hydrogen and hydrogen-hydrogen interactions in a previously reported compound, OsClH 3 (PPh 3 ) 3 . In that report, the H-H distance was measured as 1.48(2) A by single-crystal neutron diffraction and was described as a "stretched" dihydrogen. The INS-generated phonon density of states is well described by the DFT calculations, including those normal modes dominated by the vibrational motions of the three Os-ligated hydrogen atoms. Additionally, the resulting calculated electron charge density landscape indicates that there is no remnant bonding character between the Os-ligated hydrogen atoms in the compound. Based on these findings, we update the understanding of the metal-hydrogen bonding in this complex and reclassify it as a compressed dihydride. As such, OsClH 3 (PPh 3 ) 3 represents a rare example of a non-classical compressed dihydride species.

08 HYDROGEN↗
DOE OSTI2022

Thermal neutron scattering cross sections for amorphous carbon

Carbon materials are commonly found in both nuclear reactors and experimental systems. Various carbon structures occur in nuclear applications ranging from crystalline and nuclear graphite to the amorphous carbon seen in next-generation advanced reactor designs. Amorphous carbon is based on a randomized graphite-like structure and offers the unique ability to disperse impurities throughout the bulk composition. A graphite-like amorphous carbon system was modeled using the classical molecular dynamics (MD) code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). An improved version of the temperature-dependent Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) potential was used to model the carbon-carbon atomic interactions for the temperature at 300 K along with densities 1.60, 1.70, 1.85, and 2.23 g/cm{sup 3}. From the normalized velocity autocorrelation function (VACF), the phonon density of state (DOS) was then calculated as the Fourier transform of the normalized VACF. This DOS was then used as the primary input for the evaluation of the thermal scattering law (TSL, i.e. S(α,β)) and associated neutron thermal scattering cross sections. The TSL was analyzed using the Full Law Analysis Scattering System Hub (FLASSH). The amorphous structure results in shifts of the phonon DOS to lower energy modes than typically displayed for ideal crystalline graphite. This impact on the DOS is directly reflected in the TSL. Furthermore, the typical optical peak at 0.25 eV for the ideal graphite disappears for amorphous carbon, in good agreement with the expected structure. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Thermal Neutron Scattering Cross Sections for Graphitic Amorphous Carbon

Carbon materials are commonly found in both nuclear reactors and experimental systems. Various carbon structures occur in nuclear applications ranging from crystalline and nuclear graphite to the amorphous carbon seen in next-generation advanced reactor designs. Amorphous carbon is based on a randomized graphite-like structure and offers the unique ability to disperse impurities throughout the bulk composition. A graphite-like amorphous carbon system was modeled using the classical molecular dynamics (MD) code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). An improved version of the temperature-dependent Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) potential was used to model the carbon-carbon atomic interactions for the temperature at 300 K along with densities 1.60, 1.70, 1.85, and 2.23 g/cm 3 . From the normalized velocity autocorrelation function (VACF), the phonon density of state (DOS) was then calculated as the Fourier transform of the normalized VACF. This DOS was then used as the primary input for the evaluation of the thermal scattering law (TSL, i.e. S(α,β)) and associated neutron thermal scattering cross sections. The TSL was analyzed using the Full Law Analysis Scattering System Hub (FLASSH). The amorphous structure results in shifts of the phonon DOS to lower energy modes than typically displayed for ideal crystalline graphite. This impact on the DOS is directly reflected in the TSL. Furthermore, the typical features and the optical graphitic peak at 0.25 eV for the ideal graphite DOS disappear for graphite-like amorphous carbon, which shows good agreement with the expected structure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

INSPIRED: Inelastic neutron scattering prediction for instantaneous results and experimental design

Inelastic neutron scattering (INS) has unique advantages in probing how atoms vibrate and how the vibrations propagate and interact. Such dynamic information is crucial in understanding various material properties, from heat capacity, thermal conductivity, phase transitions, and chemical reactions to more exotic quantum behavior. The analysis and interpretation of the INS spectra often start from a model structure of the sample, followed by a series of calculations to obtain the simulated spectra to compare with experiments. The conventional way to perform such calculations usually requires significant time, computing resources, and specialized expertise. Here, we present a new program named INSPIRED (Inelastic Neutron Scattering Prediction for Instantaneous Results and Experimental Design), which enables users to perform rapid INS simulations in several different ways on their personal computers in just a few clicks, with the crystal structure as the only input file. Specifically, the users can choose a pre-trained symmetry-aware neural network (coupled with an autoencoder) to predict the phonon density of states (DOS), 1D S(E) and 2D S(|Q|,E) spectra for any given structure. One can also choose an existing density functional theory (DFT) calculation from a database (containing over 12,000 crystals), and quickly obtain the simulated INS spectra for single crystals and powders. It is also possible to use pre-trained universal machine learning force fields to relax a given crystal structure, calculate the phonon dispersion and DOS, and, subsequently, the INS spectra. All these functions are implemented with a PyQt graphic user interface. Finally, we expect these new tools will benefit broad user communities and significantly improve the efficiency of experiment design, execution, and data analysis for INS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

How Proton Incorporation Reshapes Lattice Dynamics In BaSnO 3 ‐Type Proton Conductors

Proton conduction in acceptor-doped perovskites is fundamentally a vibronic process: mobile H + and D + do not move independently, but dynamically co-vibrate with the surrounding oxygen–metal framework. Direct experimental evidence for this behavior is presented using in situ 119⁢ 𝑆⁢𝑛 nuclear resonance vibrational spectroscopy (NRVS) on hydrated, deuterated, and dry 𝐵⁢𝑎⁢𝑆⁢𝑛 1−𝑥⁢ 𝑌 𝑥 ⁢𝑂 3−𝛿 . Hydration induces systematic redistributions in the Sn-projected phonon density of states (PDOS), including an upshift of the first spectral moment by about 0.4 meV, indicating a stiffening of the extended Sn–O network. H/D isotopic substitution leaves the Sn-projected PDOS largely unchanged, with only subtle isotope-dependent spectral reweighting, demonstrating that protonic degrees of freedom are not localized oscillators but are embedded in collective lattice modes. These results are rationalized using a classical coupled proton–phonon oscillator model that links the observed PDOS variations to changes in effective force constants and vibrational mass terms. The model captures how H + and D + participate in cooperative lattice dynamics rather than forming isolated OH/OD entities. Overall, NRVS probes proton–lattice coupling in ceramic proton conductors and quantitatively describes how protonic defects modulate host lattice dynamics to enable phonon-assisted long-range proton transport.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An analytic and complete equation of state for condensed phase materials

Analytic equations of state (EOS) are intended to reproduce theoretical and experimental data in a single phase portion of the thermodynamic space. We devise a complete and thermodynamically consistent model with four distinct features: (1) a reference isotherm that remains thermodynamically stable, (2) a flexible specific heat model based on a fourth-order rational polynomial, (3) a Grüneisen parameter that depends on specific volume and temperature, and (4) pressure and internal energy functions that can be inverted analytically in temperature. The model aims to improve the accuracy of existing equations of state while remaining computationally efficient. To demonstrate its features, we include calibrations for single-crystal pentaerythritol tetranitrate (PETN), liquid nitromethane (NM), and hexagonal close-packed beryllium (Be) metal. The parameter optimization uses the specific heat capacity, Grüneisen parameter, and static compression curves obtained from density functional theory for the crystalline solids and molecular dynamics simulations for liquid NM. We also present a velocity autocorrelation function that yields accurate phonon densities of states for the EOS calibration from the molecular dynamics trajectories. Each of the three calibrations is constrained to enforce the ambient state from experimental measurements and validated against experimental Hugoniot data from multiple sources. We also include one-dimensional hydrodynamic simulations of the isentropic compression experiments for beryllium conducted at the Z facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Lattice dynamics of 119 Sn impurity in a bcc-Cr crystal

The chromium crystal doped with 119 Sn isotope was studied using the nuclear resonance inelastic x-ray scattering and first principles calculations. The Sn partial phonon density of states (PDOS) was obtained for three temperatures that correspond to different magnetic states of Cr. At all temperatures, the energy spectrum consists of a broad band around 18 meV and a narrow peak at 43 meV. The additional peak around 39 meV is observed only in the magnetically ordered phases, indicating the influence of magnetic order in chromium on lattice dynamics. The partial PDOS calculated with the antiferromagnetic order on Cr atoms show a very good agreement with the experimental data. It is revealed that the high-energy peak is lying above the phonon spectra of the pure bcc-Cr crystal. These are the local modes with the increased energies due to a strongly reduced distance between Sn and the nearest-neighbor Cr atoms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Structure and lattice excitations of the copper substituted lead oxyapatite Pb 9.06⁢(7)⁢ Cu 0.94⁢(6)⁢ (PO 3.92⁢(4) ) 6 ⁢O 0.96⁢(3)

The copper substituted lead oxyapatite, Pb 10-x⁢ Cu x (PO 3.92⁢(4) ) 6 ⁢O 0.96⁢(3) [x = 0.94(6)] was studied using neutron and x-ray diffraction and neutron spectroscopy techniques. The crystal structure of the main phase of our sample, which has come to be colloquially known as LK-99, is verified to possess a hexagonal structure with space group P6 3 /m, alongside the presence of impurity phases Cu and Cu 2 ⁢S. We determine the primary substitution location of the Cu as the Pb1 (6⁢ℎ) site, with a small substitution at the Pb2 (4⁢f) site. Consequently, no clear Cu-doping-induced structural distortion was observed in the investigated temperature region between 10 K and 300 K. Specially, we did not observe a reduction of coordinate number at the Pb2 site or a clear tilting of PO 4 tetrahedron. Magnetic characterization reveals a diamagnetic signal in the specimen, accompanied by a very weak ferromagnetic component at 2 K. No long-range magnetic order down to 10 K was detected by the neutron diffraction. Inelastic neutron scattering measurements did not show magnetic excitations for energies up to 350 meV. There is no sign of a superconducting resonance in the excitation spectrum of this material. The measured phonon density of states compares well with density functional theory calculations performed for the main LK-99 phase and its impurity phases. Our study may shed some insight into the role of the favored substitution site of copper in the absence of structural distortion and superconductivity in LK-99.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Determination of single-crystal elastic moduli of Li RE F 4 ( RE =Y, Gd, and Tb) by resonant ultrasound spectroscopy

The tetragonal fluoro-scheelite Li RE F 4 compounds (RE = rare earth) have been shown to exhibit a variety of useful optical and magnetic properties. While LiYF 4 has been widely studied, many of the fundamental thermodynamic properties of other members of this family remain unknown. Here, we report the complete elastic tensors (C ij ) of single-crystalline LiYF 4 , LiTbF 4 , and LiGdF 4 using resonant ultrasound spectroscopy and density functional theory (DFT). Here we compare the results for LiYF 4 with prior experimental results using time-of-flight ultrasound methods. This is the first report, however, of the experimental elastic tensors of LiTbF 4 and LiGdF 4 . The present results point to a softening of the elastic moduli of the Li RE F 4 system when Y is replaced by the larger ionic radius of Tb or Gd. Furthermore, we find that just 0.3% doping with Nd on the Y site also leads to a slight softening of the moduli. The variation of the elastic moduli as a function of temperature up to 216°C was also measured. A nearly linear softening of all seven independent elastic moduli was observed with increasing temperature. Phonon dispersions and phonon density of states obtained by DFT support the experimental finding of a significantly higher sound velocity due to lighter Y atoms in LiYF 4 , as compared to heavier LiTbF 4 and LiGdF 4 .

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

Thermal Neutron Scattering Law Evaluation for Zirconium Carbide and Critical Mass Calculations

Zirconium carbide (ZrC) is a candidate material for use in advanced high temperature reactors, including space nuclear thermal propulsion applications. Thermal neutron scattering laws (TSLs) are generated for carbon bound in ZrC, C(ZrC), and zirconium bound in ZrC, Zr(ZrC), using ab initio lattice dynamics methods. These evaluations are to be submitted for inclusion in ENDF/B-VIII.1 and use the incoherent approximation for inelastic scattering as well as the new mixed elastic scattering treatment. The application of disordered alloy theory is introduced to appropriately capture the isotopic composition of Zr and C in the elastic scattering cross section. Localized higher energy vibrations in the C(ZrC) phonon density of states that are separated from lower energy modes result in quantized oscillations in the inelastic contributions to the TSL with a significant likelihood of large energy down-scattering and up-scattering interactions, where the latter increases in probability with temperature. The quanta of energy transfer during neutron thermalization is substantially greater than classically expected within the thermal neutron energy range. MC21 critical mass calculations of ZrC mixtures with high-enriched uranium demonstrate an impact of the TSLs when compared to free-gas treatment for 235 U concentrations less than 0.2 g/cm 3 . Additional MC21 critical mass calculations with homogenous mixtures of ZrC and reactor-grade graphite also demonstrate sensitivity to the ZrC TSL for thermal spectrum driven fission systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Lattice dynamics of high-pressure hydrides studied by inelastic neutron scattering

Due to the small mass and anomalously large neutron scattering cross-section of proton (about 80 barns compared to a few barns for other nuclei), inelastic neutron scattering is considered as one of the most effective tools in studying optical vibrations of hydrogen atoms in metal hydrides. The current review is focused on the binary hydrides of 3d- and 4d-metals of groups VI–VIII, which were produced at high hydrogen pressures of several gigapascals in relatively large quantities of hundreds of mg, quenched to low temperature and studied by INS ex situ at ambient pressure with high statistical accuracy. One of the unusual effects revealed by INS is a strong increase in the strength of the metal-hydrogen interactions with decreasing atomic number of the d-metal accompanied by an increase in the Me-H distance. Based on the available experimental results, the spectra g(E) of the phonon density of states and temperature dependencies CV(T) of the heat capacity at constant volume at T up to 1000 K have been derived in this paper and presented both in the figures and in digital form. This provides the reference data for the theoretical investigations of the crystal structures and compositions of new practically important hydrides giving the opportunity to validate calculation methods by comparing the calculated g(E) and CV(T) with the accurate experimental dependencies for the binary hydrides. Recent INS studies showed [R.A. Klein et al., J. Alloy. Compd. 894 (2022) 162381] that the fingerprints of anomalously short H-H separations of 1.6 Å violating the “2 Å rule” can be easily and unambiguously identified in the complex INS spectra of quaternary hydrides (La,Ce)NiInH 1+x . This makes neutron spectroscopy an attractive means for obtaining valuable data in the search for novel hydrides with a record high hydrogen capacity.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Reversible hydrogen storage in multilayer graphane: Lattice dynamics, compressibility, and heat capacity studies

Multilayer graphane (hydride of graphite) is a crystalline hydrocarbon of composition CH, which can be synthesized from graphite and molecular hydrogen at pressures above 2GPa [V.E. Antonov et al. Carbon 100 (2016) 465]. Using X-ray diffraction, this compound was tentatively identified as the “graphane II” phase of 3D-graphane predicted by ab initio calculations [X.-D. Wen et al. PNAS 108 (2011) 6833] and consisting of layers of 2D-graphane in the “chair” conformation. When heated in a vacuum, the compound does not form any intermediate hydrocarbons and reversibly decomposes back into graphite and hydrogen at 770–920 K. In the present work, almost single-phase samples of graphite hydride and deuteride were synthesized at 7.4 GPa and 870 K. Their investigation by inelastic neutron scattering supplemented by ab initio calculations gave spectra g(E) of the phonon density of states with a gap of about 15 meV at approx. 100 meV, which is a unique identifier for the chair form of graphane. The equation of state V(P) of the hydride was studied at room temperature and hydrogen pressures up to 53 GPa by synchrotron X-ray diffraction in a diamond anvil cell. Further, the graphane II phase did not react with the surrounding hydrogen and did not undergo any phase transformations upon the compression and after heating to 1500 K at 53GPa. The high thermal and pressure stability of this exotic phase makes it an important part of the C–H system. The obtained g(E) spectra of graphite hydride and deuteride were used to calculate temperature dependences of their heat capacity. Measurements of the heat capacity at temperatures 120–673 K confirmed the good accuracy of these calculations.

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

Operando NRIXS and XAFS Investigation of Segregation Phenomena in Fe-Cu and Fe-Ag Nanoparticle Catalysts during CO 2 Electroreduction

Operando nuclear resonant inelastic X-ray scattering (NRIXS) and X-ray absorption fine-structure spectroscopy (XAFS) measurements were used to gain insight into the structure and surface composition of FeCu and FeAg nanoparticles (NPs) during the electrochemical CO 2 reduction (CO 2 RR) and to extract correlations with their catalytic activity and selectivity. The formation of a core–shell structure during CO 2 RR for FeAg NPs was inferred from the analysis of the operando NRIXS data (phonon density of states, PDOS) and XAFS measurements. Electrochemical analysis of the FeAg NPs revealed a faradaic selectivity of 36 % for CO in 0.1 M KHCO 3 at -1.1 V vs. RHE, similar to that of pure Ag NPs. In contrast, a predominant selectivity towards H 2 evolution is obtained in the case of the FeCu NPs, analogous to the results obtained for pure Fe NPs, although small Cu NPs have also been shown to favor H 2 production.

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

Broadband phonon production from axion absorption

We show that axion dark matter in the range meV ≲ m a ≲ 100 meV can incoherently excite phonons in crystal targets with unpolarised nuclear spins. This can occur through its coupling to nuclear spins and/or through its induced time-dependent electric dipole moment in nuclei. Due to the random orientation of the nuclear spins, translation symmetry is broken in the phonon effective theory, allowing axion absorption to create phonons with unrestricted momentum. The absorption rate is therefore proportional to the phonon density of states, which generically has support across a wide range of energies, allowing for a broadband detection scheme. We calculate the absorption rate for solid H 2 , D 2 , Al 2 O 3 , GaAs, H 2 O, D 2 O, Be and Li 2 O, and find that materials containing light, non-zero spin nuclei are the most promising. The predicted rates for the QCD axion are of the order of a few events / 10 kg-year exposure, setting an ambitious target for the required exposure and background suppression.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Evolutionary multi-objective optimization and Pareto-frontal uncertainty quantification of interatomic forcefields for thermal conductivity simulations

Predictive Molecular Dynamics simulations of thermal transport require forcefields that can simultaneously reproduce several structural, thermodynamic and vibrational properties of materials like lattice constants, phonon density of states, and specific heat. This requires a multi-objective optimization approach for forcefield parameterization. Existing methodologies for forcefield parameterization use ad-hoc and empirical weighting schemes to convert this into a single-objective optimization problem. Here, we provide and describe software to perform multi-objective optimization of Stillinger–Weber forcefields (SWFF) for two-dimensional layered materials using the recently developed 3rd generation non-dominated sorting genetic algorithm (NSGA-III). NSGA-III converges to the set of optimal forcefields lying on the Pareto front in the multi-dimensional objective space. This set of forcefields is used for uncertainty quantification of computed thermal conductivity due to variability in the forcefield parameters. We demonstrate this new optimization scheme by constructing a SWFF for a representative two-dimensional material, 2H-MoSe 2 and quantifying the uncertainty in their computed thermal conductivity.

97 MATHEMATICS AND COMPUTING↗