Search NASA⌕ Search

SEARCH · Search NASA

Results for “Phonon dispersion”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

First-principles investigation of uranium mononitride (UN): Effect of magnetic ordering, spin-orbit interactions and exchange correlation functional

Uranium mononitride (UN) is a promising nuclear fuel that combines the advantageous properties of readily used UO 2 and uranium alloys. Various properties of UN have been previously studied using different density functional theory (DFT) methodologies; however, there are still inconsistencies when it comes to the dynamical stability and defect properties of UN. We address these inconsistencies by studying the UN phonons and defect properties using DFT calculations employing two generalized gradient approximation (GGA) exchange-correlation functionals: PBE and AM05, with and without an added on-site Coulomb repulsion term (+U). Furthermore, we investigate the importance of spin-orbit coupling (SOC) when calculating the properties of UN. We use the different methodologies to determine the preference of UN to have antiferromagnetic (AFM) ordering, as seen in experiments, or ferromagnetic (FM) ordering of the uranium spins. We compare the crystallographic properties, density of states, the DFT X-ray photoelectron spectra and phonon dispersions calculated using the different methodologies. We demonstrate that GGA + U reproduces the AFM ordering in UN, but the crystal structure is dynamically unstable. We also show that magnetic ordering is important in finding the lowest energy defective structure, and that SOC has a distinct influence on the energy of the different uranium interstitial defects. Lastly, we discuss the point defect formation energies under U-rich and N-rich conditions, and the stoichiometric formation energies calculated with the different methodologies, providing insight into the observed tendency for forming hypostoichiometric UN.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High pressure melt line of nickel using a generalized embedded atomic method potential

As the second most abundant metal in the Earth's core, nickel plays an important role in determining the structure and temperature of the Earth's core. Yet, the melt line of Ni at pressures corresponding to the Earth's core has not been explored in the literature. Many previous experimental and simulation efforts have reported the melting point of Ni at pressures below 100 GPa, but there exist large discrepancies, most of which have persisted due to various experimental and simulation bottlenecks in handling extreme pressure and temperature conditions. We adopted the generalized embedded atom method, which overcomes the limitations of existing interatomic potentials, to probe phase stability and phase boundaries of Ni at pressures between 50 and 500 GPa. Further, the potential was validated by comparing the cold curves, phonon dispersion curves, and enthalpies of fusion with ab initio density functional theory calculations. Our analysis shows that face centered cubic (FCC) is stable, and the hexagonal close packed (HCP) and body centered cubic (BCC) phases are metastable close to the melt line. Melting temperatures at different pressures were obtained from two-phase co-existence simulations and take the following functional form: $T$ m = $1969.23+19.15P-0.012P$ 2 . In contrast to iron, differences between the melting points of the stable and metastable phases of Ni are less than 250 K at 300 GPa, and the difference in melting points of the metastable BCC and HCP phases changes sign at 500 GPa, which implies that the phase transition mechanisms during solidification can be very complex.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Kohn anomaly and elastic softening in body-centered cubic molybdenum at high pressure

Transition metals in body-centered cubic (bcc) structures under compression can display several novel physical properties because of their complex electronic structures and electron-phonon interactions. Here, we used inelastic x-ray scattering experiments in a diamond-anvil cell up to ∼45 GPa and density-functional theory calculations up to 210 GPa to investigate the phonon dispersions, and electronic and elastic properties of single-crystal molybdenum (Mo). Our results show a pressure-induced Kohn anomaly at 𝑞∼0.5 along the [ξ00] direction in the longitudinal acoustic mode at ∼45 GPa; this anomaly is triggered by the pressure-enhanced Fermi-surface nesting effect. Theoretical calculations show that electron redistributions in the 𝑠-to-𝑑 orbitals of bcc-Mo contribute to the shear modulus anomaly at ∼50 GPa. In contrast, the Young's modulus anomaly in bcc-Mo at ∼210 GPa results from a Lifshitz-type electronic topological transition. In conclusion, our results shed light on the complex electronic behaviors that are associated with macroscopic elastic properties in typical bcc 𝑑-block transition metals under compression.

3-dimensional systems↗
DOE OSTI2022

Theory of thermal properties of magnetic materials with unknown entropy

Theoretical approaches to study thermal properties of magnetic materials typically require accurate models of magnetic interactions in order to define the entropy. Here we introduce a complementary approach for examining thermal properties in magnetic systems where an accepted model for such interactions does not exist. In place of a specific model for magnetic interactions, the approach integrates measurements of temperature dependent magnetization of the studied material into a first principles computational scheme. The approach calculates system pressure from thermally disordered microstates that properly incorporate vibrational and spin subsystems at each temperature as well as the coupling between these subsystems. We apply the approach to calculate phonon modes and to investigate the anomalously low thermal expansion of the classical Invar alloy, Fe 0.65 Ni 0.35 . Here, the calculated phonon dispersions for Invar are in excellent agreement with measured data. The Invar thermal expansion is shown to remain small between 50 K and room temperature, consistent with the experimentally observed low thermal expansion value in this same temperature range. This anomalously small thermal expansion is directly connected to a small positive contribution from lattice thermal disorder that is nearly canceled by a large negative magnetic disorder contribution. By contrast, calculations for bcc Fe show a much larger thermal expansion, consistent with experiment, which is dominated by a large contribution from lattice thermal disorder that is reduced only slightly by a small negative contribution from that of magnetism. These findings give insights into the unusual nature of magnetism and spin-lattice coupling in Invar and Fe. In addition, they give promising preliminary support to the presented new methodology as a complementary way to investigate thermal properties of magnetic materials. The success achieved on Invar and Fe motivates future testing of the approach on other magnetic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electronic and magnetic properties of iridium-based novel Heusler alloys

We report half-metallicity and magnetism including exchange splitting are the most significant physical parameters to predict and design a candidate material for spintronic applications. We report here an ab-initio investigation on chemical formation and dynamical stability along with electronic structure and magnetic properties of Ir 2 Cr (Si, Ge) and IrRhCr (Si, Ge) Heusler alloys. The negative formation and cohesive energies with positive phonon dispersions confirm the stabilities of these alloys. Electronic structure calculations reveal that Ir 2 Cr (Si, Ge) and IrRhCrSi alloys are half-metallic ferromagnets with unprecedented exchange splitting. In addition, IrRhCrGe also shows semi-metallic nature. All of these materials follow Slater Pauling rule with large magnetic moments and 100% spin-polarization. With Cr bearing the majority of the local magnetic moment and exchange splitting, a ferromagnetic state is more stable than a nonmagnetic state. The electronic charge distribution and population analysis confirm mixed ionic and covalent bonding. The magnetocrystalline anisotropy energy, with the easy magnetization along the [1 1 1] direction, is significantly high in Ir 2 CrGe. Elastic constants such as shear (G), bulk (B), Young’s moduli, and Poisson’s ratio indicate that the IrRhCrSi and IrRhCrGe alloys are mechanically stable, and Ir 2 CrSi and Ir 2 CrGe are mechanically unstable. The Pugh’s (B/G) and Poisson’s ratios confirm that the stable alloys are ductile.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine learning interatomic potential for predicting the thermal properties of uranium nitride

We present a combined computational and experimental investigation of the thermal properties of uranium nitride (UN), focusing on the development of a machine learning interatomic potential (MLIP) using the moment tensor potential framework. The MLIP was trained on density functional theory (DFT) data and validated against various quantities including energies, forces, elastic constants, phonon dispersion, and defect formation energies, achieving excellent agreement with DFT calculations, prior experimental results, and our thermal conductivity measurement. The potential was then employed in molecular dynamics simulations to predict key thermal properties such as melting point, thermal expansion, specific heat, and lattice thermal conductivity. To further assess model accuracy, we fabricated a UN sample and performed new thermal conductivity measurements representative of single-crystal properties, which showed strong agreement with the MLIP predictions. This work confirms the reliability and predictive capability of the developed potential for determining the thermal properties of UN.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Lattice-Charge Coupling in a Trilayer Nickelate with Intertwined Density Wave Order

Intertwined charge and spin correlations are ubiquitous in a wide range of transition metal oxides and are often perceived as intimately related to unconventional superconductivity. Theoretically envisioned as driven by strong electronic correlations, the intertwined order is usually found to be strongly coupled to the lattice as signaled by pronounced phonon softening. Recently, both charge and spin density waves (CDW and SDW) and superconductivity have been discovered in several Ruddlesden-Popper (RP) nickelates, in particular trilayer nickelates 𝑅⁢𝐸 4⁢ Ni 3 ⁢O 10 (𝑅⁢𝐸 = Pr, La). The nature of the intertwined order and the role of lattice-charge coupling are at the heart of the debate about these materials. Using inelastic x-ray scattering, we mapped the low-energy phonon dispersions in 𝑅⁢𝐸 4⁢ Ni 3 ⁢O 10 and found no evidence of softening near the CDW wave vector over a wide temperature range, which contrasts with the pronounced anomalies frequently observed in cuprate superconductors. Calculations of the electronic susceptibility revealed a peak at the observed SDW ordering vector but not at the CDW wave vector. Our experimental and theoretical findings highlight the crucial role of the spin degree of freedom and establish a foundation for understanding the interplay between superconductivity and density-wave transitions in RP nickelate superconductors and beyond.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analysis of phase stability and chemical segregation in the Mo-V alloys using a generalized embedded atom method potential

A new interatomic potential for the Mo-V system is introduced to facilitate the study of phase stability and mechanical properties at lower temperatures. This potential is based on a generalization of the embedded atom method and includes contributions from embedding energy, explicit two- and three-body interactions and nonlocal many-body interaction terms. The parameters of the potential are optimized by using data from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, such as elastic constants, equation of states, phonon dispersion curves, point defect properties and melting temperatures for different compositions. Even though our potential is trained on a small dataset, its accuracy is comparable to available machine learning potentials for Mo and V. Furthermore, our results show that an ordered B2 phase is stable at low temperatures in alloys containing 50% V, but the solid solution phase is stable above 800 K. However, such long-range ordering is not observed in V-rich or Mo-rich alloys. In addition, our results show that V segregates to dislocation cores and grain boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗
DOE OSTI2022

Symmetry-correct bonding in density functional theory calculations for delta phase Pu

The long-held conclusion that magnetic order in delta phase Pu makes the structure mechanically unstable in density functional theory calculations is incorrect if the magnetic order is allowed to be non-collinear. The non-collinear 3Q spin structure is intrinsically cubic and makes all structurally equivalent bonds equivalent in their bonding character. Applied to density functional theory calculations on delta phase Pu, the 3Q spin structure results in elastic constants and phonons with the correct symmetry. Finally, the calculated phonon dispersion agrees better with experiment than previous calculations, and the calculated thermal expansion shows the unique behavior seen experimentally: it is negative.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Verification of stability and unraveling the electronic and physical properties of bulk and (001)-surfaces of newly synthesized Ti 2 ZnX (X = C, N) MAX phases

MAX phase family has been extended by the addition of late transition metals at the A-site with the expectation of diverse functional properties. Here, we present our systematic density functional investigation on the thermodynamic and phonon stabilities, elastic properties, including elastic constants, elastic moduli and elastic anisotropy of newly synthesized Ti 2 ZnX (X = C, N) phases in comparison with conventional Ti 2 AlX (X = C, N). Due to the smaller size of N as compared to C, the unit cell dimension is reduced when C atoms are replaced by N atoms at the X-site. Furthermore, the Ti 2 ZnC and Ti 2 ZnN are stable at the equilibrium volume of 110.84 Å 3 and 105.70 Å 3 . The thermodynamic, mechanical and dynamical stabilities are validated by estimating the formation energies, elastic constants and phonon dispersions, respectively. The elastic properties of Ti 2 ZnN are less anisotropic as compared to those of Ti 2 ZnC. To understand the thin-film characteristics in Ti 2 ZnX, the surface properties with (001)-terminated slabs are investigated. Both Ti 2 ZnX bulk and (001)-surfaces exhibit metal-like electronic structures. There is a strong covalent bonding between Ti-X and Ti-Zn atoms confirmed by the charge density map and Mulliken population analysis. Additional states are generated at the Fermi level (EF) due to the unusual d-p states hybridization between Ti and Zn atoms. The anisotropy in chemical bonding is confirmed by the cleavage energy difference between Ti-X and Ti-Zn. Here, Ti(X)-001 and Zn-001 terminations are stable surfaces; however, in terms of chemical potentials, Zn-001 termination is the most favourable in Ti 2 ZnX.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Structural Phase Transitions of NbO 2 : Bulk versus Surface

The metal to insulator transition of NbO 2 has been predicted to be a result of a structural phase transition (SPT) governed by Peierls physics. However, direct observation of the SPT using experimental techniques is still restricted by the extremely high transition temperature (810 °C) and the proclivity for NbO 2 to oxidize into Nb 2 O 5 above 400 °C when exposed to air. In this work, we address these issues and employ temperature-dependent X-ray spectroscopy to describe the SPT of NbO 2 from the bulk to surface. Temperature-dependent extended X-ray absorption fine structure spectroscopy (T-EXAFS) reveals a gradual weakening of the bulk Nb dimers over a large temperature range, which is indicative of a second-order Peierls mechanism. From these measurements, we determine the critical dimer distance to be 2.77 Å. Our T-EXAFS observations are supported by density functional theory of the phonon dispersion and the electronic density of states of NbO 2 , which conclude that the dimerization is responsible for the insulating phase. The dimerization does not extend to the topmost layers, where an oxygen rich surface reconstruction is preferred irrespective of temperature even in extremely reducing environments; changes in the low-energy electron diffraction patterns are attributed to oxygen concentration and are independent of the underlying bulk phase transitions of NbO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Role of Phonons and Oxygen Vacancies in Non-Cubic SrVO 3

Combining neutron diffraction with pair distribution function analysis, we have uncovered hidden reduced symmetry in the correlated metallic d 1 perovskite, SrVO 3 . Specifically, we show that both the local and global structures are better described using a GdFeO 3 distorted (orthorhombic) model as opposed to the ideal cubic ABO 3 perovskite type. Recent reports of imaginary phonon frequencies in the density functional theory (DFT)-calculated phonon dispersion for cubic SrVO 3 suggest a possible origin of this observed non-cubicity. Namely, the imaginary frequencies computed could indicate that the cubic crystal structure is unstable at T = 0 K. However, our DFT calculations provide compelling evidence that point defects in the form of oxygen vacancies, and not an observable symmetry breaking associated with calculated imaginary frequencies, primarily result in the observed non-cubicity of SrVO 3 . These experimental and computational results are broadly impactful because they reach into the thin-film and theoretical communities who have shown that SrVO 3 is a technologically viable transparent conducting oxide material and have used SrVO 3 to develop theoretical methods, respectively.

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

Helical Organic and Inorganic Polymers

Despite being a staple of synthetic plastics and biomolecules, helical polymers are scarcely studied with Gaussian-basis-set ab initio electron-correlated methods on an equal footing with molecules. This article introduces an ab initio second-order many-body Green’s-function [MBGF(2)] method with nondiagonal, frequency-dependent Dyson self-energy for infinite helical polymers using screw-axis-symmetry-adapted Gaussian-spherical-harmonics basis functions. Together with the Gaussian-basis-set density-functional theory for energies, analytical atomic forces, translational-period force, and helical-angle force, it can compute correlated energy, quasiparticle energy bands, structures, and vibrational frequencies of an infinite helical polymer, which smoothly converge at the corresponding oligomer results. These methods can handle incommensurable structures, which have an infinite translational period and are hard to characterize by any other method, just as efficiently as commensurable structures. We apply them to polyethylene (2/1 helix), polyacetylene (Peierls’ system), and polytetrafluoroethylene (13/6 helix) to establish the quantitative accuracy of MBGF(2)/cc-pVDZ in simulating their (angle-resolved) ultraviolet photoelectron spectra, and of B3LYP/cc-pVDZ or 6-31G** in reproducing their structures, infrared and Raman band positions, phonon dispersions, and (coherent and incoherent) inelastic neutron scattering spectra. We then predict the same properties for infinitely catenated chains of nitrogen or oxygen and discuss their possible metastable existence under ambient conditions. They include planar zigzag polyazene (N2)x (Peierls’ system), 11/3-helical isotactic polyazane (NH)x, 9/4-helical isotactic polyfluoroazane (NF)x, and 7/2-helical polyoxane (O)x as potential high-energy-density materials.

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

Optical and Plasmonic Properties of High-Electron-Density Epitaxial and Oxidative Controlled Titanium Nitride Thin Films

The present paper reports on the fabrication, detailed structural characterizations, and theoretical modeling of titanium nitride (TiN) and its isostructural oxide derivative, titanium oxynitride (TiNO) thin films that have excellent plasmonic properties and that also have the potential to overcome the limitation of noble metal and refractory metals. The TiNO films deposited at 700 °C in high vacuum conditions have the highest reflectance (R = ~ 95%), largest negative dielectric constant (ε 1 = –161), and maximal plasmonic figure of merit (FoM = –ε 1 /ε 2 ) of 1.2, followed by the 600 °C samples deposited in a vacuum (R = ~ 85%, ε 1 = –145, FoM = 0.8) and 700 °C–5 mTorr sample (R = ~ 82%, ε 1 = –8, FoM = 0.3). To corroborate our experimental observations, we calculated the phonon dispersions and Raman active modes of TiNO by using the virtual crystal approximation. From the experimental and theoretical studies, a multilayer optical model has been proposed for the TiN/TiNO epitaxial thin films for obtaining individual complex dielectric functions from which many other optical parameters can be calculated. The advantages of oxide derivatives of TiN are the continuation of similar free electron density as in TiN and the acquisition of additional features such as oxygen-dependent semiconductivity with a tunable bandgap.

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

Multi-reward Reinforcement Learning Based Bond-Order Potential to Study Strain-Assisted Phase Transitions in Phosphorene

Here, we introduce a multi-reward reinforcement learning (RL) approach to train a flexible bond-order potential (BOP) for 2D phosphorene based on ab initio training data sets. Our approach is based on a continuous action space Monte Carlo tree search algorithm that is general and scalable and presents an efficient multiobjective optimization scheme for high-dimensional materials design problems. As a proof-of-concept, we deploy this scheme to parametrize multiple structural and dynamical properties of 2D phosphorene polymorphs. Our RL-trained BOP model adequately captures the structure, energetics, transformation barriers, equation of state, elastic constants, and phonon dispersions of various 2D P polymorphs. We use this model to probe the impact of temperature and strain rate on the phase transition from black (α-P) to blue phosphorene (β-P) through molecular dynamics simulations. A decrease in critical strain for this phase transition with increase in temperature is observed, and the underlying atomistic mechanisms are discussed.

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

Signatures of Coherent Phonon Transport in Ultralow Thermal Conductivity Two-Dimensional Ruddlesden–Popper Phase Perovskites

An emerging class of methylammonium lead iodide (MAPbI 3 )-based Ruddlesden-Popper (RP) phase perovskites, BA 2 MA n-1 Pb n I 3n+1 (n = 1-7), exhibit enhanced stability to environmental conditions relative to MAPbI 3 , yet still degrade at elevated temperatures. Here, we experimentally determine the thermal conductivities of these layered RP phases for n = 1-6, where n defines the number of repeated perovskite octahedra per layer. We measure thermal conductivities of 0.37 ± 0.13/0.12, 0.17 ± 0.08/0.07, 0.21 ± 0.05/0.04, and 0.19 ± 0.04/0.03 W/m∙K in thin films of n = 1-4 and 0.08 ± 0.06/0.04, 0.06 ± 0.04/0.03, 0.06 ± 0.03/0.03, and 0.08 ± 0.07/0.04 W/m∙K in single crystals of n = 3-6. With the exception of n = 1, these thermal conductivities are lower than the range of 0.34-0.50 W/m∙K reported for single-crystal MAPbI 3 . Reduced-order lattice dynamics modeling suggests that the initially decreasing trend of thermal conductivity in similarly oriented perovskites with increasing n may result from the transport properties of coherent phonons, emergent from the superstructure, that do not scatter at the interfaces of organic butylammonium chains and perovskite octahedra. Reduced group velocity of coherent phonons in n = 3-6, a consequence of band flattening in the phonon dispersion, is primarily responsible for their ultralow thermal conductivities. Similar effects on thermal conductivity have been experimentally demonstrated in deposited superlattices, but never in naturally defined materials such as RP phases. GIWAXS measurements reveal that higher n RP phase thin films are less orientationally controlled and therefore possess apparently elevated thermal conductivities relative to single crystals of the same n.

2D perovskites↗