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

Pressure-induced superconductivity in the hydrogen-rich pseudobinary CaB - H n compounds

Here, the crystal structures of CaB – Hn compounds with n = 1 – 12 in a pressure range of 50–300 GPa were studied using the genetic algorithm method and first-principles density-function theory calculations. Stable structures with stoichiometry of CaB H 6 and CaB H 7 were predicted in different pressure range. BH 4 , BH 5 , and BH 6 units were found to be the main motifs in these compounds. Moreover, metastable Imm2 CaB H 7 is dynamically stable above 180 GPa, with the formation of tetrahedral BH 4 unit surrounded by Ca atom and H 3 unit. Electron-phonon coupling (EPC) calculations reveal that the superconducting properties are closely related to the strong hydrogen-boron bonding of the BH 4 unit in Imm2, and T c can reach ~200 K at 200 GPa. As the major units, BH 4 , BH 5 , and BH 6 units exist in the Imm2, P 2 , and P2 1 /m phases of CaB H 7 at 300 GPa, respectively; the corresponding EPC parameter λ decreases with the increase of hydrogen content in CaB H 7 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Interplay of coulomb and exciton–phonon coupling controls singlet fission dynamics in two pentacene polymorphs

Pentacene is an important model organic semiconductor in both the singlet exciton fission (SF) and organic electronics communities. We have investigated the effect of changing crystal structure on the SF process, generating multiple triplet excitons from an initial singlet exciton, and subsequent triplet recombination. Unlike for similar organic semiconductors that have strong SF sensitive to polymorphism, we find almost no quantitative difference between the kinetics of triplet pair (TT) formation in the two dominant polymorphs of pentacene. Both pairwise dimer coupling and momentum-space crystal models predict much faster TT formation from the bright singlet excited state of the Bulk vs ThinFilm polymorph, contrasting with the experiment. GW and Bethe–Salpeter equation calculations, including exciton–phonon coupling, reveal that ultrafast phonon-driven transitions in the ThinFilm polymorph compensate the intrinsically slower purely Coulomb-mediated TT formation channel, rationalizing the similarity in observed rates. Taking into account the influence of subtle structural distinctions on both the direct and phonon-mediated SF pathways reveals a predictive capability to these methods, expected to be applicable to a wide variety of molecular crystals.

14 SOLAR ENERGY↗
DOE OSTI2023

Electron-phonon coupling from GW perturbation theory: Practical workflow combining BerkeleyGW, ABINIT, and EPW

Here we present a workflow of practical calculations of electron-phonon (e-ph) coupling with many-electron correlation effects included using the GW perturbation theory (GWPT). This workflow combines BerkeleyGW, ABINIT, and EPW software packages to enable accurate e-ph calculations at the GW self-energy level, going beyond standard calculations based on density functional theory (DFT) and density-functional perturbation theory (DFPT). This workflow begins with DFT and DFPT calculations (ABINIT) as starting point, followed by GW and GWPT calculations (BerkeleyGW) for the quasiparticle band structures and e-ph matrix elements on coarse electron k- and phonon q-grids, which are then interpolated to finer grids through Wannier interpolation (EPW) for computations of various e-ph coupling determined physical quantities such as the electron self-energies or solutions of anisotropic Eliashberg equations, among others. A gauge-recovering symmetry unfolding technique is developed to reduce the computational cost of GWPT (as well as DFPT) while fulfilling the gauge consistency requirement for Wannier interpolation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Monte-Carlo modeling of phonon thermal transport using DFT-based anisotropic dispersion relations over the full Brillouin zone

In this work, we present a Monte Carlo (MC) approach to solve the phonon Boltzmann transport equation (BTE) in which the anisotropic phonon dispersion relations over the full Brillouin zone (BZ) are used. In this approach, the discretization of the BZ used to compute the phonon relaxation time places constraints on the direction of scattered phonons in the real-space simulation domain. The phonon dispersion and phonon relaxation times are calculated using the density functional theory (DFT) approach. The modified MC approach is validated by a close examination of its ability to simulate phonon transport in both the ballistic and diffusive regimes for multiple materials including GaAs, InAs, ThO 2 , and α-U. In doing so, the phonon thermal conductivities from 100 K to 1000 K are calculated and compared with traditional non-transport solution of the phonon BTE. It is found that the phonon thermal conductivities of α-U and ThO 2 obtained from MC simulations using isotropic dispersion are larger than the values obtained using anisotropic phonon dispersion relations over the full BZ. The effect of phonon-defect scattering on the thermal conductivity of ThO 2 is also studied as an application of the current MC approach and found to agree with previously computed values in the literature. The MC solver developed here has been parallelized as a step to demonstrate its potential to solving computationally intensive phonon thermal transport problems at the mesoscale.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High critical field superconductivity at ambient pressure in MoB 2 stabilized in the P6/mmm structure via Nb substitution

Recently it was discovered that, under elevated pressures, MoB 2 exhibits superconductivity at a critical temperature T c as high as 32 K. The superconductivity appears to develop following a pressure-induced structural transition from the ambient pressure $\text{R}\bar{3}$⁢m structure to an MgB 2 -like P6/mmm structure. This suggests that remarkably high T c values among diborides are not restricted to MgB 2 as previously appeared to be the case, and that similarly high T c values may occur in other diborides if they can be coerced into the MgB 2 structure. In this paper, we show that density functional theory calculations indicate that phonon free energy stabilizes the P6/mmm structure over the $\text{R}\bar{3}$⁢m at high temperatures across the Nb 1–x ⁢Mo x B 2 series. X-ray diffraction confirms that the synthesized Nb-substituted MoB 2 adopts the MgB 2 crystal structure. Finally, high magnetic field electrical resistivity measurements and specific heat measurements demonstrate that Nb x ⁢Mo 1–x⁢ B 2 exhibits superconductivity with T c as high as 8 K and critical fields approaching 6 T.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Structure of high-lying levels populated in the 96 Y → 96 Zr β decay

The nature of the high-lying final levels of the 96 Y gs β decay, one of the three most important contributors to the high-energy reactor antineutrino spectrum, has been investigated in high-resolution γ-ray spectroscopy following the β decay as well as in a campaign of inelastic photon scattering experiments. The comprehensive approach establishes 1 - levels associated with the Pygmy Dipole Resonance as high-lying final levels in the β decay. Branching ratios extracted from β decay complement photon scattering and allow the absolute E1 excitation strength to be determined for levels populated in both reactions. The combined data represents a comprehensive approach to the wavefunction of the 1 - levels below the Qβ value, which are investigated in the Quasiparticle Phonon Model. The calculations reveal that the components populated in β decay contribute only with small amplitudes to the complex wavefunction of these 1 - levels. A comparison of the β decay results to data from total absorption γ -ray spectroscopy demonstrates a good agreement between both measurements.

Pygmy dipole resonance↗
DOE OSTI2021

High-Pressure Synthesis of Bulk Cobalt Cementite, Co 3 C

Transition metal carbides find use in a wide range of advanced high-resilience applications including high-strength steels, heat shields, and deep-earth drills. However, carbides of the mid-to-late transition metals remain difficult to isolate and characterize on account of their metastability, which precludes the preparation of high-quality bulk single crystal samples using traditional solid-state methods. Herein, we report a combined computational and experimental survey of the cobalt-carbon binary system under high pressures and demonstrate that pressure offers a route toward the bulk synthesis of the metastable cementite-type cobalt carbide, Co 3 C, which under ambient conditions can only be prepared in low-dimensional thin film or nanoparticle forms. First-principles calculations reveal two competitive low-energy stoichiometric phases under ambient pressures - Pnnm-Co 2 C (Fe 2 C-type) and Pnma-Co 3 C (Fe 3 Ctype) - consistent with the known low-dimensional phases that have been studied for their promising magnetic properties. However, the calculated formation enthalpy of Pnma-Co 3 C decreases steadily with the applied pressure, while that of Pnnm-Co 2 C increases. Here, we pursue these results using high-pressure laser-heated synthesis methods coupled with in situ X-ray diffraction and observe the formation of Pnma-Co 3 C above 4.8 GPa. We determine the experimental bulk modulus of Co 3 C to be K o = 237 GPa (K p = 4.0). First-principles calculations of the phonon modes in Co 3 C reveal dynamical instabilities at ambient pressure that are absent under compression. These results offer a promising new route for the synthesis of rare-earth-free magnets.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Large-gap insulating dimer ground state in monolayer IrTe2

Abstract Monolayers of two-dimensional van der Waals materials exhibit novel electronic phases distinct from their bulk due to the symmetry breaking and reduced screening in the absence of the interlayer coupling. In this work, we combine angle-resolved photoemission spectroscopy and scanning tunneling microscopy/spectroscopy to demonstrate the emergence of a unique insulating 2 × 1 dimer ground state in monolayer 1 T -IrTe 2 that has a large band gap in contrast to the metallic bilayer-to-bulk forms of this material. First-principles calculations reveal that phonon and charge instabilities as well as local bond formation collectively enhance and stabilize a charge-ordered ground state. Our findings provide important insights into the subtle balance of interactions having similar energy scales that occurs in the absence of strong interlayer coupling, which offers new opportunities to engineer the properties of 2D monolayers.

36 MATERIALS SCIENCE↗
DOE OSTI2020

ABINIT: Overview and focus on selected capabilities

ABINIT is probably the first electronic-structure package to have been released under an open-source license about 20 years ago. It implements density functional theory, density-functional perturbation theory (DFPT), many-body perturbation theory (GW approximation and Bethe–Salpeter equation), and more specific or advanced formalisms, such as dynamical mean-field theory (DMFT) and the “temperature-dependent effective potential” approach for anharmonic effects. Relying on planewaves for the representation of wavefunctions, density, and other space-dependent quantities, with pseudopotentials or projector-augmented waves (PAWs), it is well suited for the study of periodic materials, although nanostructures and molecules can be treated with the supercell technique. The present article starts with a brief description of the project, a summary of the theories upon which ABINIT relies, and a list of the associated capabilities. It then focuses on selected capabilities that might not be present in the majority of electronic structure packages either among planewave codes or, in general, treatment of strongly correlated materials using DMFT; materials under finite electric fields; properties at nuclei (electric field gradient, Mössbauer shifts, and orbital magnetization); positron annihilation; Raman intensities and electro-optic effect; and DFPT calculations of response to strain perturbation (elastic constants and piezoelectricity), spatial dispersion (flexoelectricity), electronic mobility, temperature dependence of the gap, and spin-magnetic-field perturbation. The ABINIT DFPT implementation is very general, including systems with van der Waals interaction or with noncollinear magnetism. Community projects are also described: generation of pseudopotential and PAW datasets, high-throughput calculations (databases of phonon band structure, second-harmonic generation, and GW computations of bandgaps), and the library libpaw. ABINIT has strong links with many other software projects that are briefly mentioned.

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

Precisely computing phonons via irreducible derivatives

Computing phonons from first principles is typically considered a solved problem, yet inadequacies in existing techniques continue to yield deficient results in systems with sensitive phonons. Here, in this study, we circumvent this issue using the lone irreducible derivative (LID) and bundled irreducible derivative (BID) approaches to computing phonons via finite displacements, where the former optimizes precision via energy derivatives and the latter provides the most efficient algorithm using force derivatives. A condition number optimized basis for BID is derived which guarantees the minimum amplification of error. Additionally, a hybrid LID-BID approach is formulated, in which select irreducible derivatives computed using LID replace BID results. We illustrate our approach on two prototypical systems with sensitive phonons: the shape memory alloy AuZn and metallic lithium. Comparing our resulting phonons in the aforementioned crystals to calculations in the literature reveals nontrivial inaccuracies. Our approaches can be fully automated, making them well suited for both niche systems of interest and high-throughput approaches.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Lamellar fluctuations melt ferroelectricity.

We consider a standard Ginzburg-Landau model of a ferroelectric whose electrical polarization is coupled to gradients of elastic strain. At the harmonic level, such flexoelectric interaction is known to hybridize acoustic and optic phonon modes and lead to phases with modulated lattice structures that precede the state with spontaneously broken inversion symmetry. Here, we use the self-consistent phonon approximation to calculate the effects of thermal and quantum polarization fluctuations on the bare hybridized modes to show that such long-range modulated order is unstable at all temperatures. We discuss the implications for the nearly ferroelectric SrTiO3 and KTaO3, and we propose that these systems are melted versions of an underlying modulated state that is dominated by nonzero momentum thermal fluctuations except at the very lowest temperatures.

Guzman-Verri, G. G.↗
DOE OSTI2024

Deep-ultraviolet transparent conducting SrSnO 3 via heterostructure design

Exploration and advancements in ultrawide bandgap (UWBG) semiconductors are pivotal for next-generation high-power electronics and deep-ultraviolet (DUV) optoelectronics. Here, we used a thin heterostructure design to facilitate high conductivity due to the low electron mass and relatively weak electron-phonon coupling, while the atomically thin films ensured high transparency. We used a heterostructure comprising SrSnO 3 /La:SrSnO 3 /GdScO 3 (110), and applied electrostatic gating, which allow us to effectively separate charge carriers in SrSnO 3 from dopants and achieve phonon-limited transport behavior in strain-stabilized tetragonal SrSnO 3 . This led to a modulation of carrier density from 10 18 to 10 20 cm −3 , with room temperature mobilities ranging from 40 to 140 cm 2 V −1 s −1 . The phonon-limited mobility, calculated from first principles, closely matched experimental results, suggesting that room temperature mobility could be further increased with higher electron density. In addition, the sample exhibited 85% optical transparency at a 300-nm wavelength. These findings highlight the potential of heterostructure design for transparent UWBG semiconductor applications, especially in DUV regime.

Science & Technology - Other Topics↗
DOE OSTI2025

An experimental and computational investigation of the structure and spectroscopic signatures of α -UO 3

α-UO 3 is a common intermediate compound found in the nuclear fuel cycle, yet the exact crystal structure of this material has long been debated. Inconsistent computational and experimental data in previous works has led to varying conclusions between authors. Furthermore, to ensure the validity of our results in this work, the structural and spectroscopic signatures of pure phase α-UO 3 are investigated using powder X-ray diffraction and optical vibrational spectroscopy (infrared and Raman). Rietveld refinement of powder X-ray diffraction data on pure phase α-UO 3 collected in this work allows us to propose an alteration to the currently accepted C2mm structure (a = 3.9705 Å, b = 6.8553 Å, c = 4.15955 Å, α = β = γ = 90°) for α-UO 3 with no uranyl [UO 2 2+ ] bonds. Raman spectra collected using two excitation wavelengths (two instruments using 532 nm and one 785 nm) are presented, and differences with recently published results are discussed. Infrared spectra from two instruments used here agree well with recently published results, but the spectral range encompassed in our data extends past what has been reported with modern techniques. Additionally, we provide tentative vibrational mode assignments based on density functional perturbation theory calculations and resulting phonon eigenvector visualizations. Unexpected features in the optical vibrational spectra of α-UO 3 are explained by unique features in the structure we present.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRSDate not supplied

Towards Accurate Predictions of Martensitic Transition Temperatures for Shape Memory Alloys from Ab Initio Simulations

Experimentally, NiTi undergoes a single martensitic phase transition around 341 K from the lowtemperature (T) monoclinic B19’ phase (P21/m) to the high-temperature cubic B2 phase (Pm3m). Theoretically, an orthorhombic B33 (Cmcm) has also been proposed as the T=0 ground state structure, although this phase has never been observed. Accurate predictions of martensitic transition temperatures (MTT) have remained elusive in part due to several well-known theoretical complexities of these systems including low temperature instabilities of the B2 phase. Recently, we proposed a rigorous thermodynamic integration approach based on ab initio simulations to resolve many of these difficulties [1,2]. However, an unsatisfying overprediction of the MTT relative to experiment (by ~100 K) means a fully quantitative theory is still lacking. In this work, we report several new developments to our method that bring first principles theory and experiment much closer into agreement. Our calculations indicate that phonon free energies at low temperature stabilizes B19’ over B33, rationalizing B19’ as the ground state down to T=0. We also find that accurate computations of the electronic free energy, i.e. the change in energy and the appearance of electronic configurational entropy due to finite temperature, is crucial to obtain accurate MTT. Incorporating these corrections results in an MTT prediction of 365 K for binary NiTi, which is in very close agreement with experiment. Our theoretical approach is expected to be a broadly applicable and predictive theory for MTT of SMAs.

Zhigang Wu↗
DOE OSTI2023

Computational Study of Mixing Solid Materials for CO2 Capture Technology

CO2 is one of the major combustion products that, once released into the air, can contribute to global climate change. Solid sorbents have been reported as promising candidates for CO2 sorbent applications due to their high CO2 absorption capacities at moderate working temperatures. By combining thermodynamic database mining with first-principles density functional theory and phonon lattice dynamics calculations, the National Energy Technology Laboratory has proposed and validated a theoretical screening methodology to identify the most promising CO2 sorbent candidates from a vast array of possible solid materials. The advantage of this method is that it identifies the thermodynamic properties of the CO2 capture reaction as a function of temperature and gas pressure without any experimental input beyond crystallographic structural information of the solid phases involved. The calculated thermodynamic properties of solid materials versus temperature and pressure changes were further used to evaluate the equilibrium properties for the CO2 adsorption/desorption cycles. The selected CO2 sorbent candidates were further considered for experimental validations. In this presentation, the results of M2O/Al2O3 (M=Li, Na, K) mixtures (MAlO2 and M5AlO4) capturing CO2 will be demonstrated in detail.

Duan, Yuhua↗
DOE OSTI2024

Generation of Enrichment-Dependent Thermal Neutron Scattering Data

This work details the generation of enrichment-dependent thermal neutron scattering cross sections for several crucial uranium fuel compounds. The evaluations of the thermal scattering law (TSL) and associated cross sections for uranium dioxide (UO 2 ), uranium carbide (UC), and uranium nitride (UN) were performed using standard ab initio lattice dynamics (AILD) methods. The data for uranium metal was produced using a novel hybrid approach of molecular dynamics combined with lattice dynamics methods. 235 U enrichments of 5%, 10% (LEU+), 19.75% (HALEU), 93% (HEU), and 100% were considered, in addition to natural uranium. The enrichment-dependent masses and free atom cross sections were used in the generation of elastic and inelastic thermal neutron scattering cross sections, while the calculation of the phonon density of states (DOS) and resulting TSL considered only the natural isotopic composition of uranium. The use of an identical DOS for all enrichments is expected to have minimal impact on the final data, as the small change in uranium mass should not significantly affect lattice vibrations. The cross sections are shown to exhibit significant dependence on 235 U enrichment. The submission of this data to the National Nuclear Data Center (NNDC) for release in the ENDF/B-VIII.1 database should support the design of advanced reactor concepts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

PyDislocDyn: A Python code for calculating dislocation drag and other crystal properties

PyDislocDyn is a suite of python programs designed to perform various calculations for dislocation dynamics in the continuum limit. In particular, one of its main purposes is to calculate dislocation drag from phonon wind. Additional features include the averaging of elastic constants for polycrystals, the calculation of the dislocation field including its limiting velocities, and the calculation of dislocation self-energy and line tension.

36 MATERIALS SCIENCE↗