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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Fast and accurate machine learning prediction of phonon scattering rates and lattice thermal conductivity

Abstract Lattice thermal conductivity is important for many applications, but experimental measurements or first principles calculations including three-phonon and four-phonon scattering are expensive or even unaffordable. Machine learning approaches that can achieve similar accuracy have been a long-standing open question. Despite recent progress, machine learning models using structural information as descriptors fall short of experimental or first principles accuracy. This study presents a machine learning approach that predicts phonon scattering rates and thermal conductivity with experimental and first principles accuracy. The success of our approach is enabled by mitigating computational challenges associated with the high skewness of phonon scattering rates and their complex contributions to the total thermal resistance. Transfer learning between different orders of phonon scattering can further improve the model performance. Our surrogates offer up to two orders of magnitude acceleration compared to first principles calculations and would enable large-scale thermal transport informatics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ab initio anharmonic thermodynamic properties of cubic Ca Si O 3 perovskite

We present an ab initio study of the thermodynamic properties of cubic CaSiO 3 perovskite (CaPv) over the pressure and temperature range of the Earth's lower mantle. We compute the anharmonic phonon dispersions throughout the Brillouin zone by utilizing the phonon quasiparticle approach, which characterizes the intrinsic temperature dependence of phonon frequencies and, in principle, captures full anharmonicity. Such temperature-dependent phonon dispersions are used to calculate ab initio free energy in the thermodynamic limit (N → ∞) within the framework of the phonon gas model. Accurate free energy calculations enable us to investigate cubic CaPv's thermodynamic properties, e.g., thermal expansivity, Grüneisen parameter, bulk modulus, heat capacity and thermal equation of state, where anharmonic effects are demonstrated. Furthermore, the present methodology provides an important theoretical approach for exploring phase boundaries, thermodynamic, and thermoelastic properties of strongly anharmonic materials at high pressures and temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Niobium in the clean limit: An intrinsic type-I superconductor

Niobium is one of the most researched superconductors, both theoretically and experimentally. It is enormously significant in all branches of superconducting applications, from powerful magnets to quantum computing. It is therefore imperative to understand its fundamental properties in great detail. In this work, we use the results of recent microscopic calculations of anisotropic electronic, phonon, and superconducting properties, and apply a thermodynamic criterion for the type of superconductivity, more accurate and straightforward than a conventional Ginzburg-Landau parameter κ-based delineation, to show that pure niobium is a type-I superconductor in the clean limit. However, disorder (impurities, defects, strain, stress) pushes it to become a type-II superconductor.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Electronic, direct optical, and phonon-assisted optical properties of 4H Si from first principles

The cubic polytype of silicon (Si) is the most commercialized semiconductor material and finds applications in numerous electronic and optoelectronic devices, such as solar cells. However, recent reports on the synthesis of the hexagonal 4H Si polytype have attracted the attention of the scientific community to understand its functional properties. Here, we report the electronic, vibrational, and optical properties of the 4H Si polytype obtained with predictive first-principles calculations, with an emphasis of phonon-assisted absorption in the indirect regime. Compared to the cubic polytype, 4H Si shows a slightly narrower indirect gap by ∼0.05 eV. The calculated phonon-assisted optical spectra show that 4H Si exhibits a stronger absorption coefficient than cubic Si across the visible and IR spectral regions. We further evaluate the short-circuit current density of textured thin-films. We demonstrate that 4H Si can achieve the same short-circuit current density for a five times thinner film compared to the cubic polytype, which mainly resulted from absorption in the indirect gap regime. Our work demonstrates the advantages of 4H Si for thin-film silicon-based solar-cell applications.

Materials Science↗
NASA NTRS1992

Thermal conductivity of thermoelectric Si(0.8)-Ge(0.2) alloys

The thermal conductivity of heavily doped, n-type Si-Ge alloys has been studied from 300 to 1200 K. The scattering rate of several phonon scattering mechanisms has been calculated, including intrinsic scattering, mass defect and distortion scattering, phonon-electron scattering, and scattering by inclusions. These rates were then used to calculate the lattice thermal conductivity. The electronic component of the thermal conductivity was calculated from the calculated Lorenz ratio and measured values of the electrical conductivity. The total thermal conductivity was then compared to measured values for a specimen studied by Vining et al. (1991).

White, D. P.↗
DOE OSTI2021

Onset of anharmonicity and thermal conductivity in SnSe

The anharmonicity in SnSe is investigated through the analysis of moments of 119Sn nuclear resonant inelastic x-ray scattering and ab initio molecular dynamics calculations. Experimental evidences show that the anharmonic behavior started around 300 K, substantially lower than the usually suggested structural transition at 800 K. Both theory and experiments reveal substantial lifetime broadening and frequency renormalization of the optical phonons. Additionally, thermal conductivities calculated from the temporal energy moment using the Einstein diffusion equation are in good agreement with previous experiments. The abrupt increase of the thermal power near 800 K is driven by an electronic factor and not by the enhanced anharmonicity due to structural change.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

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

Temperature-dependent lattice dynamics in iridium

The characterization of simple elemental systems is key to benchmarking first-principles modeling of electronic and vibrational behaviors of materials. A large body of literature has been built for most elemental systems in the periodic table. However, surprisingly little neutron work has been performed to probe the vibrational properties of iridium, likely due to its large neutron absorption cross section. Nonetheless, iridium is of significant scientific and technological interest due to large relativistic electron effects and electron-phonon coupling, particularly in strongly correlated iridate compounds. In this paper, we report temperature-dependent inelastic neutron scattering measurements of the vibrational properties of iridium, from which we extract key thermodynamic properties. To overcome the challenge of the large neutron absorption of iridium, we developed a simple post-processing correction procedure. The measured densities of phonon states compare well with quasiharmonic density functional theory calculations, although the obtained experimental phonon Grüneisen parameters are much larger than expected, reaching as high as γ=4.5, indicating substantial anharmonicity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Electron-phonon drag enhancement of transport properties from a fully coupled ab initio Boltzmann formalism

We present a combined treatment of the nonequilibrium dynamics and transport of electrons and phonons by carrying out ab initio calculations of the fully coupled electron and phonon Boltzmann transport equations. We find that the presence of mutual drag between the two carriers causes the thermopower to be enhanced and dominated by the transport of phonons, rather than electrons as in the traditional semiconductor picture. Drag also strongly boosts the intrinsic electron mobility, thermal conductivity and the Lorenz number. Impurity scattering is seen to suppress the drag enhancement of the thermal and electrical conductivities, while having weak effects on the enhancement of the Lorenz number and thermopower. We demonstrate these effects in n-doped 3C-SiC at room temperature, and explain their origins. This work establishes the roles of microscopic scattering mechanisms in the emergence of strong drag effects in the transport of the interacting electron-phonon gas.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Clathrate BaNi 2 P 4 : An Interplay of Heat and Charge Transport Due to Strong Host–Guest Interactions

Heat and charge transport properties of the metallic unconventional clathrate BaNi 2 P 4 , hosting Ba cations in oversized Ni 8 P 16 cages, are investigated. A novel method of single-crystal growth was developed, yielding 2–3 mm sized crystals of BaNi 2 P 4 . We also developed a setup to accurately measure thermal conductivity and electrical resistivity of the synthesized single crystals in a wide temperature range avoiding crystal remounting. BaNi 2 P 4 has a metallic temperature dependence of its electrical resistivity (decreasing with decreasing temperature) and manifests an unconventional T 2 power law for 50 K < T < 300 K; below 50 K, the power-law exponent increases gradually such that below 10 K the power law is T 5 , a predicted but extremely rarely experimentally observed dependence for peculiar electron–phonon interactions. Furthermore, electronic band structure calculations, consistent with measurements of de Haas–van Alphen oscillations, show large band dispersions with significant contributions of Ba orbitals to states near the Fermi level, which is atypical for clathrates. The thermal properties of BaNi 2 P 4 were probed using a combination of variable-temperature single-crystal X-ray diffraction experiments, heat capacity measurements, first-principles phonon dispersion calculations, and inelastic neutron scattering measurements. BaNi 2 P 4 exhibits significant hybridization of the Ba-guest and Ni–P-framework vibrational modes, which may be enhanced via the detected split of the Ba position, which results in strong Ba-framework interactions.

36 MATERIALS SCIENCE↗
NASA NTRS1983

Prediction of structures and magnetic orientations in solid alpha and beta-O2

A quasi-harmonic-lattice-dynamics method coupled with a pattern-recognition optimization scheme is used to determine the minimum energy structures and magnetic orientations of solid oxygen. It is shown that the magnetic interaction is responsible for the stability of alpha-O2 with respect to beta-O2 at zero temperature and pressure. The calculated alpha-O2 lattice parameters, magnetic orientations, and sublimation energy are in good agreement with experiment. Phonon dispersion curves are calculated at vector k not equal to zero and the acoustic sound velocities are determined. The rms translational and orientational fluctuations from equilibrium are also calculated. The beta-O2 phase is described by constraining the magnetic moments so that the magnetic Hamiltonian preserves the hexagonal symmetry of the crystal. The calculated lattice parameters are in good agreement with the experiments, and a three-sublattice, quasi-helical magnetic orientation is predicted from structural and energetic considerations.

Etters, R. D.↗
DOE OSTI2021

Spectral properties and enhanced superconductivity in renormalized Migdal-Eliashberg theory

Migdal-Eliashberg theory describes the properties of the normal and superconducting states of electron-phonon-mediated superconductors based on a perturbative treatment of the electron-phonon interactions. It is necessary to include both electron and phonon self-energies self-consistently in Migdal-Eliashberg theory in order to match numerically exact results from determinantal quantum Monte Carlo in the adiabatic limit. Here in this work we provide a method to obtain the real-axis solutions of the Migdal-Eliashberg equations with electron and phonon self-energies calculated self-consistently. Our method avoids the typical challenge of computing cumbersome singular integrals on the real axis and is numerically stable and exhibits fast convergence. Analyzing the resulting real-frequency spectra and self-energies of the two-dimensional Holstein model, we find that self-consistently including the lowest-order correction to the phonon self-energy significantly affects the solution of the Migdal-Eliashberg equations. The calculation captures the broadness of the spectral function, renormalization of the phonon dispersion, enhanced effective electron-phonon coupling strength, minimal increase in the electron effective mass, and the enhancement of superconductivity which manifests as a superconducting ground state despite strong competition with charge-density-wave order. We discuss surprising differences in two common definitions of the electron-phonon coupling strength derived from the electron mass and the density of states, quantities which are accessible through experiments such as angle-resolved photoemission spectroscopy and electron tunneling. An approximate upper bound on 2Δ/T c for conventional superconductors mediated by retarded electron-phonon interactions is proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Soft-mode enhanced type-I superconductivity in LiPd 2 Ge

The synthesis, crystal structure, and physical properties (magnetization, resistivity, heat capacity) in combination with theoretical calculations of the electronic structure and phonon properties are reported for intermetallic compounds LiPd 2 X ( X = Si, Ge, and Sn). LeBail refinement of powder x-ray diffraction data confirms that all compounds belong to the Heusler family (space group $\textit{F m-3m}$, No. 225). The lattice parameter increases with atomic size of X , and its value varies from $\textit{a}$ = 5.9059(4) Å for LiPd 2 Si and $\textit{a}$ = 6.0082(3)Å for LiPd 2 Ge, to $\textit{a}$ = 6.2644(1) Å for LiPd 2 Sn. The first compound, LiPd 2 Si, has apparently not been previously reported. All measured quantities demonstrate that LiPd 2 Ge exhibits superconductivity below $T_c$ = 1.96 K and the normal- and superconducting-state data indicate that it is a weak-strength type-I superconductor ($C/γT_c$ = 1.38) with electron-phonon coupling constant $λ_{e–p}$ = (0.53–0.56). LiPd 2 Si and LiPd 2 Sn are not superconducting above 1.68 K. The experimental observations are supported by theoretical calculations which show that LiPd 2 Ge has the highest computed $λ_{e–p}$ and $T_c$ of the group. A strong softening of the acoustic phonon mode is calculated, and in the case of X = Ge and Sn, imaginary phonon frequencies were computed. In this work, the soft mode is most pronounced in the case of LiPd 2 Ge, which suggests its correlation with superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1991

A model for the high-temperature transport properties of heavily doped n-type silicon-germanium alloys

A model is presented for the high-temperature transport properties of large-grain-size, heavily doped n-type silicon-germanium alloys. Electron and phonon transport coefficients are calculated using standard Boltzmann equation expressions in the relaxation time approximation. Good agreement with experiment is found by considering acoustic phonon and ionized impurity scattering for electrons, and phonon-phonon, point defect, and electron-phonon scattering for phonons. The parameters describing electron transport in heavily doped and lightly doped materials are significantly different and suggest that most carriers in heavily doped materials are in a band formed largely from impurity states. The maximum dimensionless thermoelectric figure of merit for single-crystal, n-type Si(0.8)Ge(0.2) at 1300 K is estimated at ZT about 1.13 with an optimum carrier concentration of n about 2.9 x 10 to the 20th/cu cm.

Vining, Cronin B.↗
DOE OSTI2022

express: Extensible, high-level workflows for swifter ab initio materials modeling

In this work, we introduce an open-source Julia project, express, an extensible, lightweight, high-throughput, high-level workflow framework that aims to automate ab initio calculations for the materials science community. express is shipped with well-tested workflow templates, including structure optimization, equation of state (EOS) fitting, phonon spectrum (lattice dynamics) calculation, and thermodynamic property calculation in the framework of the quasi-harmonic approximation (QHA). It is designed to be highly modularized so that its components can be reused across various occasions, and customized workflows can be built on top of that. Users can also track the status of workflows in real-time, and rerun failed jobs thanks to the data lineage feature express provides. Finally, two working examples, i.e., all workflows applied to lime and akimotoite, are also presented in the code and this paper.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Anomalous thermal transport under high pressure in boron arsenide

High pressure represents extreme environments and provides opportunities for materials discovery. Thermal transport under high hydrostatic pressure has been investigated for more than 100 years and all measurements of crystals so far have indicated a monotonically increasing lattice thermal conductivity. Here we report in situ thermal transport measurements in the newly discovered semiconductor crystal boron arsenide, and observe an anomalous pressure dependence of the thermal conductivity. We use ultrafast optics, Raman spectroscopy and inelastic X-ray scattering measurements to examine the phonon bandstructure evolution of the optical and acoustic branches, as well as thermal conductivity under varied temperatures and pressures up to 32 gigapascals. Using atomistic theory, we attribute the anomalous high-pressure behaviour to competitive heat conduction channels from interactive high-order anharmonicity physics inherent to the unique phonon bandstructure. Our study verifies ab initio theory calculations and we show that the phonon dynamics—resulting from competing three-phonon and four-phonon scattering processes—are beyond those expected from classical models and seen in common materials. In conclusion, this work uses high-pressure spectroscopy combined with atomistic theory as a powerful approach to probe complex phonon physics and provide fundamental insights for understanding microscopic energy transport in materials of extreme properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Temperature-dependent phonon lifetimes and thermal conductivity of silicon by inelastic neutron scattering and ab initio calculations

In this paper, inelastic neutron scattering on a single crystal of silicon was performed at temperatures from 100 to 1500 K. These experimental data were reduced to obtain phonon spectral intensity at all wave vectors $\overset{→}{Q}$ and frequencies ω in the first Brillouin zone. Thermal broadenings of the phonon peaks were obtained by fitting and by calculating with an iterative ab initio method that uses thermal atom displacements on an ensemble of superlattices. Agreement between the calculated and experimental broadenings was good, with possible discrepancies at the highest temperatures. Distributions of phonon widths versus phonon energy had similar shapes for computation and experiment. These distributions grew with temperature but maintained similar shapes. Parameters from the ab initio calculations were used to obtain the thermal conductivity from the Boltzmann transport equation, which was in good agreement with experimental data. Despite the high group velocities of longitudinal acoustic phonons, their shorter lifetimes reduced their contribution to the thermal conductivity, which was dominated by transverse acoustic modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗