Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phonon calculations”

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 OSTI2024

Phonon Screening of Excitons in Atomically Thin Semiconductors

Atomically thin semiconductors, encompassing both 2D materials and quantum wells, exhibit a pronounced enhancement of excitonic effects due to geometric confinement. Consequently, these materials have become foundational platforms for the exploration and utilization of excitons. Recent ab initio studies have demonstrated that phonons can substantially screen electron-hole interactions in bulk semiconductors and strongly modify the properties of excitons. While excitonic properties of atomically thin semiconductors have been the subject of extensive theoretical investigations, the role of phonon screening on excitons in atomically thin structures remains unexplored. In this Letter, we demonstrate via ab initio GW-Bethe-Salpeter equation calculations that phonon screening can have a significant impact on optical excitations in atomically thin semiconductors. We further show that the degree of phonon screening can be tuned by structural engineering. We focus on atomically thin GaN quantum wells embedded in AlN and identify specific phonons in the surrounding material, AlN, that dramatically alter the lowest-lying exciton in monolayer GaN via screening. Our studies provide new intuition beyond standard models into the interplay among structural properties, phonon characteristics, and exciton properties in atomically thin semiconductors, and have implications for future experiments.

2-dimensional systems↗
DOE OSTI2023

Electron–phonon physics from first principles using the EPW code

EPW is an open-source software for ab initio calculations of electron–phonon interactions and related materials properties. The code combines density functional perturbation theory and maximally localized Wannier functions to efficiently compute electron–phonon coupling matrix elements, and to perform predictive calculations of temperature-dependent properties and phonon-assisted quantum processes in bulk solids and low-dimensional materials. Here, we report on significant developments in the code since 2016, namely: a transport module for the calculation of charge carrier mobility under electric and magnetic fields using the Boltzmann transport equation; a superconductivity module for calculations of phonon-mediated superconductors using the anisotropic multi-band Eliashberg theory; an optics module for calculations of phonon-assisted indirect transitions; a module for the calculation of small and large polarons without supercells; and a module for calculating band structure renormalization and temperature-dependent optical spectra using the special displacement method. For each capability, we outline the methodology and implementation and provide example calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multi-channel direct detection of light dark matter: theoretical framework

We present a unified theoretical framework for computing spin-independent direct detection rates via various channels relevant for sub-GeV dark matter — nuclear re- coils, electron transitions and single phonon excitations. Despite the very different physics involved, in each case the rate factorizes into the particle-level matrix element squared, and an integral over a target material- and channel-specific dynamic structure factor. We show how the dynamic structure factor can be derived in all three cases following the same procedure, and extend previous results in the literature in several aspects. For electron transitions, we incorporate directional dependence and point out anisotropic target materials with strong daily modulation in the scattering rate. For single phonon excitations, we present a new derivation of the rate formula from first principles for generic spin-independent couplings, and include the first calculation of phonon excitation through electron couplings. We also discuss the interplay between single phonon excitations and nuclear recoils, and clarify the role of Umklapp processes, which can dominate the single phonon production rate for dark matter heavier than an MeV. Our results highlight the complementarity between various search channels in probing different kinematic regimes of dark matter scattering, and provide a common reference to connect dark matter theories with ongoing and future direct detection experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Engineering Anomalously Large Electron Transport in Topological Semimetals

Abstract Anomalous transport of topological semimetals has generated significant interest for applications in optoelectronics, nanoscale devices, and interconnects. Understanding the origin of novel transport is crucial to engineering the desired material properties, yet their orders of magnitude higher transport than single‐particle mobilities remain unexplained. This work demonstrates the dramatic mobility enhancements result from phonons primarily returning momentum to electrons due to phonon‐electron dominating over phonon–phonon scattering. Proving this idea, proposed by Peierls in 1932, requires tuning electron and phonon dispersions without changing symmetry, topology, or disorder. This is achieved by combining de Haas ‐ van Alphen (dHvA), electron transport, Raman scattering, and first‐principles calculations in the topological semimetals MX 2 (M = Nb, Ta and X = Ge, Si). Replacing Ge with Si brings the transport mobilities from an order magnitude larger than single particle ones to nearly balanced. This occurs without changing the crystal structure or topology and with small differences in disorder or Fermi surface. Simultaneously, Raman scattering and first‐principles calculations establish phonon–electron dominated scattering only in the MGe 2 compounds. Thus, this study proves that phonon‐drag is crucial to the transport properties of topological semimetals and provides insight to engineer these materials further.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ab Initio Evaluation of Plutonium Dioxide $S(α,β)$ and Thermal Neutron Cross Sections

Plutonium Dioxide (PuO 2 ) is an advanced fuel-cycle nuclear fuel material for thermal and fast reactors, with a melting point at 3261.8 K and thermal conductivity above 3 Wm -1 K -1 for temperatures below 1500 K. Plutonium from spent UO 2 is reprocessed to form chemically stable mixed oxide (MOX) fuel. In contrast to metallic fuels, the low PuO 2 conductivity results in high thermal gradients between fuel center and surface, resulting in structural annealing and efficient fission product migration over burnup. Crystal binding affects the doppler broadening of epithermal resonances for nuclear fuels and scattering behavior of thermal neutrons. Currently, there is no ENDF/B Thermal Scattering Law, i.e., $S(α,β)$, evaluation for PuO 2 , only free atom cross sections for Plutonium and Oxygen, respectively, which neglect lattice contributions to total cross section. Herein, ab initio lattice dynamics (AILD) techniques are employed to calculate the phonon density of states (DOS) using spin-orbit-coupling density functional theory (DFT) to predict the PuO2 paramagnetic ground-state structure. These results highlight the evaluation of $S(α,β)$ for Plutonium and Oxygen, respectively in PuO 2 and consequential generation of thermal neutron scattering cross sections for high fidelity criticality safety analysis and reactor calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

First principles calculations in support of Pu aging: calculating the effects of lattice imperfections on thermodynamics

Plutonium (Pu) has, in theory, well defined crystal structures: its atoms are arranged in regular spatial patterns. But Pu is radioactive, and as its nuclei decay those regular spatial patterns are interrupted. The interruptions are lattice imperfections, which are known to affect how materials respond to their environment. To adequately model Pu, we need to know which lattice imperfections are present, how they interact with each other, and how they affect the material’s response to its environment. The work presented here aims to use density functional theory (DFT) calculations to begin to answer the latter, in particular, how individual lattice imperfections affect measurable effects including thermal expansion (the change in volume in response to a change in temperature), heat capacity (the amount of thermal energy needed to change a material’s temperature), and elastic moduli (a material’s resistance to applied stresses). Pu poses many computational challenges. The $\textit{f}$ electrons require special attention. Of all the elements, Pu has the largest number of electrons that must be included in the calculations. Thermal effects require calculating the phonons (the lattice vibrations), which for systems containing lattice imperfections demand large, complex computational cells - but computational resources limit the size and complexity. With careful restructuring of how the calculations are performed, all these challenges have been met to enable calculations that provide insight into how lattice imperfections affect Pu’s response to its environment. Reported here are the computational challenges and the advances developed to meet them, along with the first results showing the strong effect that one prototype lattice imperfection (an interstitial Pu atom in a delta-phase Pu lattice) has on Pu’s response to its environment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ab Initio Evaluation of Uranium Carbide S(α,β) and Thermal Neutron Cross Sections

Uranium Carbide (UC) is a nuclear fuel material which offers better neutron economy and lower fuel-cycle costs compared to conventional mixed-oxide fuels. UC’s lattice binding and dynamical properties impact thermal neutron scattering and low temperature epithermal resonance absorption. The Thermal Scattering Law (TSL) describes the scattering system available energy and momentum transfer states. There is no TSL evaluation for UC in the ENDF/B-VIII.0 database; herein, ab-initio lattice dynamics (AILD) techniques are invoked to calculate the phonon spectrum for UC using spin-orbit-coupling density functional theory (DFT). The TSLs, inelastic and elastic thermal scattering cross sections for Uranium and Carbon in UC, respectively, are calculated in FLASSH for use in higher fidelity reactor design calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗
DOE OSTI2021

First principles study of the stability and thermal conductivity of novel Li-Be hybrid ceramics

Hybrid Li-Be ceramics, combining both tritium (T) breeder (Li) and neutron multiplier (Be) for use as part of the fuel cycle of future nuclear fusion reactors are proposed. The development of such hybrid materials may reduce thermal gradients through better matching of thermal properties, mitigating the detrimental effects that may accompany traditional breeder systems while maintaining acceptable neutron multiplication and T breeding. Here, first-principles methods are used to investigate stability and thermal transport properties of a set of compounds containing both Li and Be. It is demonstrated that BeLi 2 O 4 Ge and BeLi 2 O 4 Si are mechanically stable and have formation energies comparable with leading candidates for solid state breeder materials, Li 2 TiO 3 , Li 2 ZrO 3 . It is also demonstrated that similar to the leading candidates, these compounds are insulators with thermal transport defined by phonons. The calculated thermal conductivity of BeLi 2 O 4 Ge is slightly higher compared to Li 2 TiO 3 or Li 2 ZrO 3 while in the BeLi 2 O 4 Si it is almost three times higher compared to the rest of compounds due to a higher phonon group velocities and increased phonon lifetimes. These results indicate that hybrid Li-Be ceramics offer a potential route towards better matching of thermal properties with minimal functional property degradation, thereby offering better overall fuel cycle performance.

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

RbB 3 Si 3 : An Alkali Metal Borosilicide that is Metastable and Superconducting at 1 atm

The stability, electronic structure and potential superconductivity in AB 3 Si 3 (A = Na, K, Rb, and Cs) compounds that assume the clathrate-based sodalite structure whose frameworks consists of covalent B-Si bonds is investigated via first-principles calculations. This structure type has recently been predicted in a number of high temperature superconducting hydrides, but these are only stable under megabar pressures. Herein, we predict a novel superconducting phase, RbB 3 Si 3 , that could be synthesized under pressures that are a factor of ten smaller, ~10 GPa, and quenched to atmospheric conditions. Electron-phonon coupling calculations predict that RbB 3 Si 3 possesses a superconducting critical temperature, T c , of 14 K at 1 atm. Here, the dynamic stability of RbB 3 Si 3 and CsB 3 Si 3 at ambient pressure can be explained by considering the chemical pressure exerted on the B-Si framework that is caused by the size effect of the alkali metal atom.

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

Phonon dynamics in the chiral Kitaev spin liquid

Here, we investigate the effect of a magnetic field on the extended Kitaev spin-liquid state through phonon dynamics. Using a constrained fermionic self-consistent mean field method, we analyze the quantum spin liquid (QSL) ground state for the extended Kitaev model with both the Zeeman term and the perturbative three-spin interaction term 𝜅. Our results demonstrate the dependence of the stability of the Kitaev QSL state on the field direction, consistent with findings in the literature. Additionally, we calculate the phonon dynamics for acoustic phonons coupled to the Majorana fermion excitations of the Kitaev spin-liquid state, discussing the temperature and field evolution of these quantities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Topology driven and soft phonon mode enabled Na-ion diffusion in quaternary chalcogenides, Na 3 ZnGaX 4 (X = S, and Se)

The compounds Na 3 ZnGaX 4 (X = S, Se) are potential solid electrolyte materials in sodium-based batteries, which have certain advantages over oxide materials and have shown significant ionic conductivity at ambient temperature. In this paper, we bring out atomic-level features of the diffusion process in these new materials using the microscopic techniques of inelastic neutron scattering (INS), quasielastic neutron scattering (QENS), and ab initio molecular dynamics (AIMD) simulations. The insights obtained from these techniques are unique and not available from other macroscopic experiments. Neutron scattering experiments have been performed at temperatures from 100 to 700 K. The simulations have been carried out up to 900 K. We have calculated the phonon spectra and the space–time correlation functions and found good agreement with the results of the neutron scattering experiments. The simulations enable detailed analysis of the atomic-site dependent dynamical information. We observe low-energy phonon modes of ~6 meV involving the vibrations of certain Na atoms in the lattice. This reveals that the Na at 32g Wyckoff sites (Na2) has sufficiently shallow potential among the two available crystallographic sites. This shallow potential facilitates diffusion. Furthermore, the specific structural topology of the network of interconnected zig-zag chains of the Na2 atomic sites provides the low-barrier energy pathways for diffusion. A small fraction of vacancy defects appears essential for diffusion. We further observe that the Na2 atoms undergo jump-like diffusion to the vacant next or the 2nd next neighbour sites at ~4 Å. While the QENS experiments reveal the jump-like diffusion and its time scale, detailed analysis of the AIMD simulations shows that the jumps appear mostly along zig-zag chains of the Na2 sites in the tetragonal ab-plane, as well as between the chains along the c-axis. In conclusion, the net diffusion is essentially 3-dimensional, with little anisotropy despite the anisotropy of the tetragonal crystal structure.

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

Anharmonic multi-phonon nonradiative transition: An ab initio calculation approach

Nonradiative carrier recombinations at deep centers in semiconductors are of great importance for both fundamental physics and device engineering. Here, we provide a revised analysis of Huang’s original nonradiative multi-phonon (NMP) theory with ab initio calculations. First, we confirmed at the first-principles level that Huang’s concise formula gives the same results as the matrix-based formula, and that Huang’s high-temperature formula provides an analytical expression for the coupling constant in Marcus theory. Secondly, we correct for anharmonic effects by taking into account local phonon-mode variations for different charge states of a defect. The corrected capture rates for defects in GaN and SiC agree well with experiments.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Thickness and Spin Dependence of Raman Modes in Magnetic Layered Fe 3 GeTe 2

2D layered Fe 3 GeTe 2 has attracted increasing attention due to its high magnetic ordering temperature and novel physical properties. Lattice dynamics is a fundamental property of Fe 3 GeTe 2 , and its relationships with the number of layers and interlayer spin ordering have not yet been explored in depth. Here, by first-principles density functional theory calculations, the phonon vibrations and Raman intensities of Fe 3 GeTe 2 are systematically studied from the bulk to monolayer structures. Furthermore, the spin-phonon coupling effect is investigated by considering different interlayer magnetic orderings: ferromagnetic and antiferromagnetic. It is found that the frequencies of Raman modes in Fe 3 GeTe 2 exhibit considerable dependence on the layer number and spin ordering. In this work, the results not only reveal the notable spin-phonon interactions in Fe 3 GeTe 2 , but also demonstrate that Raman modes can be utilized for characterizing the sample thickness and interlayer spin ordering in this 2D magnet.

2D magnets↗
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↗
DOE OSTI2020

Computationally Directed Discovery of MoBi 2

Incorporating bismuth, the heaviest element stable to radioactive decay, into new materials enables the creation of emergent properties such as permanent magnetism, superconductivity, and nontrivial topology. Understanding the factors that drive Bi reactivity is critical for the realization of these properties. Using pressure as a tunable synthetic vector, we can access unexplored regions of phase space to foster reactivity between elements that do not react under ambient conditions. Furthermore, combining computational and experimental methods for materials discovery at high-pressures provides broader insight into the thermodynamic landscape than can be achieved through experiment alone, informing our understanding of the dominant chemical factors governing structure formation. Herein, we report our combined computational and experimental exploration of the Mo–Bi system, for which no binary intermetallic structures were previously known. Using the ab initio random structure searching (AIRSS) approach, we identified multiple synthetic targets between 0–50 GPa. Highpressure in situ powder X-ray diffraction experiments performed in diamond anvil cells confirmed that Mo–Bi mixtures exhibit rich chemistry upon the application of pressure, including experimental realization of the computationally predicted CuAl 2 -type MoBi 2 structure at 35.8(5) GPa. Electronic structure and phonon dispersion calculations on MoBi 2 revealed a correlation between valence electron count and bonding in high-pressure transition metal–Bi structures as well as identified two dynamically stable ambient pressure polymorphs. Our study demonstrates the power of the combined computational–experimental approach in capturing highpressure reactivity for efficient materials discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗