SEARCH · Engineering Papers

Results for “Boltzmann transport”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2024

High-throughput screening of 2D materials identifies p-type monolayer WS2 as potential ultra-high mobility semiconductor

Abstract 2D semiconductors offer a promising pathway to replace silicon in next-generation electronics. Among their many advantages, 2D materials possess atomically-sharp surfaces and enable scaling the channel thickness down to the monolayer limit. However, these materials exhibit comparatively lower charge carrier mobility and higher contact resistance than 3D semiconductors, making it challenging to realize high-performance devices at scale. In this work, we search for high-mobility 2D materials by combining a high-throughput screening strategy with state-of-the-art calculations based on the ab initio Boltzmann transport equation. Our analysis singles out a known transition metal dichalcogenide, monolayer WS 2 , as the most promising 2D semiconductor, with the potential to reach ultra-high room-temperature hole mobilities in excess of 1300 cm 2 /Vs should Ohmic contacts and low defect densities be achieved. Our work also highlights the importance of performing full-blown ab initio transport calculations to achieve predictive accuracy, including spin–orbital couplings, quasiparticle corrections, dipole and quadrupole long-range electron–phonon interactions, as well as scattering by point defects and extended defects.

Chemistry↗

DOE OSTI2021

Oxidation-enhanced thermoelectric efficiency in a two-dimensional phosphorene oxide

We performed density functional theory calculations to investigate the thermoelectric properties of phosphorene oxide (PO) expected to form by spontaneous oxidation of phosphorene. Since thermoelectric features by nature arise from the consequences of the electron-phonon interaction, we computed the phonon-mediated electron relaxation time, which was fed into the semiclassical Boltzmann transport equation to be solved for various thermoelectric-related quantities. It was found that PO exhibits superior thermoelectric performance compared with its pristine counterpart, which has been proposed to be a candidate for the use of future thermoelectric applications. We revealed that spontaneous oxidation of phosphorene leads to a significant enhancement in the thermoelectric properties of n-doped phosphorene oxide, which is attributed to the considerable reduction of lattice thermal conductivity albeit a small decrease in electrical conductivity. Our results suggest that controlling oxidation may be utilized to improve thermoelectric performance in nanostructures, and PO can be a promising candidate for low-dimensional thermoelectric devices.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Relaxation time approximations in PAOFLOW 2.0

Abstract Regardless of its success, the constant relaxation time approximation has limited validity. Temperature and energy dependent effects are important to match experimental trends even in simple situations. We present the implementation of relaxation time approximation models in the calculation of Boltzmann transport in 2.0 and apply those to model band-structures. In addition, using a self-consistent fitting of the model parameters to experimental conductivity data, we provide a flexible tool to extract scattering rates with high accuracy. We illustrate the approximations using simple models and then apply the method to GaAs, Si, $$\hbox {Mg}_3\hbox {Sb}_2$$ Mg 3 Sb 2 , and $$\hbox {CoSb}_3$$ CoSb 3 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Ga 4 C-family crystals, a new generation of star thermoelectric materials, achieved by band degeneracies, valley anisotropy, and strong phonon scattering among others

Using density functional theory combined with the Boltzmann transport equation, we uncover a new class of high-performance thermoelectric materials, i.e., supertetrahedral Ga 4 C-family materials characterized by an ultrahigh thermoelectric figure of merit (ZT), which is produced both by a large power factor and ultralow lattice thermal conductivity. The former is contributed by multiple electronic band degeneracies, flat bands and valley anisotropy, while the latter is contributed both by strong phonon scattering and low phonon group velocity. Furthermore, we find that tensile strain engineering may further suppress lattice thermal conductivity, while well retaining the multi-band degeneracies, valley anisotropy and flat bands. In addition to the band degeneracy, the scattering mechanism of the system can also dominate the power factor in the process of strain implementation. Consequently, a novel thermoelectric phenomenon, i.e., a high thermoelectric figure of merit, is observed in Ca 4 C with a ZT value larger than 4.4 under a small strain ratio (0.75%). In conclusion, our theoretical studies not only uncover a new class of high-performance thermoelectric materials integrating multiple effective mechanisms, but also establish a realistic material plateau to investigate the competitive effect of multiple factors on enhancing the thermoelectric figure of merit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

High electron mobility of Al x Ga1− x N evaluated by unfolding the DFT band structure

We calculate the alloy-disorder-limited electron mobility of AlxGa1−xN from first principles. AlxGa1−xN is a technologically important ultra-wide-bandgap alloy with promise in light emitting diodes and high-power transistors. Alloying introduces statistical disorder, which causes electrons to scatter between different crystal-momentum states, leading to a reduction in mobility for intermediate alloy compositions. The corresponding lifetime, which appears as an energy broadening in the band structure, can be evaluated by unfolding the band structure from the supercell basis to the primitive-cell basis. We fit the first-principles band broadening with a model scattering potential and evaluate the low-field electron mobility using the semiclassical Boltzmann transport equation in the relaxation-time approximation. Our calculated mobility is in agreement with experimental values. We also find the lowest alloy-scattering electron mobility (total electron mobility) across the entire composition range to be 186 cm2/V s (136 cm2/V s), which is comparable to the highest electron mobility predicted in the competitor system, β-(AlxGa1−x)2O3. Our results elucidate the intrinsic limits imposed by alloy disorder on electron transport in AlxGa1−xN.

Pant, Nick (ORCID:0000000254506045)↗

DOE OSTI2021

Thermal conductivity of α -U with point defects

Here, we develop a theoretical model for thermal conductivity of α-U that combines density functional theory calculations and the coupled electron–phonon Boltzmann transport equation. The model incorporates both electron and phonon contributions to thermal conductivity and achieves good agreement with experimental data over a wide temperature range. The dominant scattering mechanism governing thermal transport in α-U at different temperatures is examined. By including phonon–defect and electron–defect scatterings in the model, we study the effect of point defects including U-vacancy, U-interstitial, and Zr-substitution on the thermal conductivity of α-U. The degradation of anisotropic thermal conductivity due to point defects as a function of defect concentration, defect type, and temperature is reported. This model provides insights into the impact of defects on both phonon and electron thermal transport. It will promote the fundamental understanding of thermal transport in α-U and provide a ground for investigation of coupled electron–phonon transport in metallic materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Phonon Olympics: Phonon property and lattice thermal conductivity benchmarking from open-source packages

Three widely used open-source packages for determining phonon properties and lattice thermal conductivities (ALAMODE, phono3py, and ShengBTE) are benchmarked by teams of expert users and the package developers. The phonons for Ge, RbBr, monolayer MoSe 2 , and AlN are modeled at zero temperature, and they scatter through three-phonon and phonon-isotope processes, with thermal conductivities obtained from the linearized Peierls–Boltzmann transport equation with input from density functional theory calculations. Over a wide range of temperatures, the thermal conductivities calculated by the teams fall within at most ±15% of their mean values for each of the four materials. The phonon frequencies, obtained from the harmonic force constants, do not show large differences between the calculations, indicating that the modal heat capacities and group velocities are not responsible for the thermal conductivity variations. It is the lifetimes associated with three-phonon scattering, obtained from the cubic force constants, that drive the variations. The many decisions required to calculate the cubic force constants (e.g., supercell size, atomic displacement, neighbor cutoff, and application of symmetries) make identification of the precise origin of the thermal conductivity variations challenging. The calculated thermal conductivities do not generally show agreement with experimental measurements, which is attributed to the limitations of the density functional theory calculations. Guidance for the development of best practices is provided, which will help to standardize protocols needed for building thermal conductivity databases. The results provide a baseline for future benchmarking of other packages and more advanced calculations.

McGaughey, Alan J. H. [Carnegie Mellon Univ., Pitt↗

DOE OSTI2026

A Proof of the Asymptotic Variance of Path Length Estimators for Single-Collision Monte Carlo Source Iteration in the Thick Diffusion Limit

Here, we prove a theorem relating the variance of path length estimators for single-collision Monte Carlo source iteration to a parameter that becomes infinitesimally small in an important physical regime arising in radiative transfer. In our usage, “single-collision Monte Carlo source iteration” refers to Monte Carlo Boltzmann transport methods in which each Monte Carlo particle history includes no more than a single collision, and the physics of multiple scattering is modeled by lagging the scattering source term and iterating until this term converges. Our theorem can be used to construct variance reduction techniques which improve the order of the estimator variance. This enables calculations that would otherwise require impractically large sample sizes to achieve practical estimator uncertainties. We believe this is the first postulation of a theorem relating estimator variance to a limiting case parameter for single-collision Monte Carlo source iteration, and the first proof of such a theorem. We illustrate the theorem’s value with an example in which the authors of a transport method used the theorem to design a variance reduction technique that improved the uncertainty of their solution by a factor of about 500 for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material.

Mathematics and Computing↗

DOE OSTI2021

Assessment of empirical interatomic potential to predict thermal conductivity in ThO 2 and UO 2

Computing vibrational properties of crystals in the presence of complex defects often necessitates the use of (semi-)empirical potentials, which are typically not well characterized for perfect crystals. In this study we explore the efficacy of a commonly used embedded-atomempirical interatomic potential for the U x Th 1- x O 2 system, to compute phonon dispersion, lifetime, and branch specific thermal conductivity. Our approach for ThO 2 involves using lattice dynamics and the linearized Boltzmann transport equation to calculate phonon transport properties based on second and third order force constants derived from the empirical potential and from first-principles calculations. For UO 2 , to circumvent the accuracy issues associated with first-principles treatments of strong electronic correlations, we compare results derived from the empirical interatomic potential to previous experimental results. It is found that the empirical potential can reasonably capture the dispersion of acoustic branches, but exhibits significant discrepancies for the optical branches, leading to overestimation of phonon lifetime and thermal conductivity. The branch specific conductivity also differs significantly with either first-principles based results (ThO 2 ) or experimental measurements (UO 2 ). These findings suggest that the empirical potential needs to be further optimized for robust prediction of thermal conductivity both in perfect crystals and in the presence of complex defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Refining the two-band model for highly compensated semimetals using thermoelectric coefficients

In studying compensated semimetals, the two-band model has proven extremely useful in capturing electrical conductivity under magnetic field, as a function of density and mobility of electronlike and holelike carriers. However, it rarely offers practical insight into magnetothermoelectric properties. Here, in this work, we report the field dependence of thermoelectric (TE) coefficients in a highly compensated semimetal NbSb 2 , where we find the Seebeck (𝑆 𝑥⁢𝑥 ) and Nernst (𝑆 𝑥⁢𝑦 ) coefficients increase quadratically and linearly with applied magnetic field, respectively. Such field dependence was predicted in previous work that studied a system of two parabolic bands, within semiclassical Boltzmann transport theory when the following two conditions are simultaneously met: 𝜔 𝑐 ⁢𝜏 ≫ 1 and tan ⁡𝜃 𝐻 ≪ 1, where 𝜔 𝑐 and 𝜏 refer to the cyclotron frequency and relaxation time, respectively, and 𝜃 𝐻 is the Hall angle. Under these conditions, we find the field dependence of the TE coefficients directly provides a relation between the electronlike (𝑛 𝑒 ) and holelike (𝑛 ℎ ) carrier densities, which in turn can be used to refine two-band model fitting. With this, we find the compensation factor ($\frac{|Δ⁢𝑛|}{𝑛_{𝑒}}$, where Δ⁢𝑛 = 𝑛 𝑒 −𝑛 ℎ ) of NbSb 2 is two orders of magnitude smaller than what was found in unrestricted fitting, resulting in a larger saturation-field scale for magnetoresistance. Within the same framework of the semiclassical theory, we also deduce that the thermoelectric Hall angle tan⁡ 𝜃 𝛾 = $\frac{𝑆_{𝑥⁢𝑦}}{𝑆_{𝑥⁢𝑥}}$ can be expressed as ($\frac{Δ⁢𝑛}{𝑛_{𝑒}}$ ×𝜔 𝑐 ⁢𝜏) −1 , which serves as a parameter to predict the degree of compensation. Our findings offer crucial insights into identifying empirical conditions for field-induced enhancement of TE performance and into engineering-efficient thermoelectric devices based on semimetallic materials.

electrical conductivity↗

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↗

DOE OSTI2021

Magnetotransport in semiconductors and two-dimensional materials from first principles

Here we demonstrate a first-principles method to study magnetotransport in materials by solving the Boltzmann transport equation (BTE) in the presence of an external magnetic field. Our approach employs ab initio electron-phonon interactions and takes spin-orbit coupling into account. We apply our method to various semiconductors (Si and GaAs) and two-dimensional (2D) materials (graphene) as representative case studies. The magnetoresistance, Hall mobility, and Hall factor in Si and GaAs are in very good agreement with experiments. In graphene, our method predicts a large magnetoresistance, consistent with experiments. Analysis of the steady-state electron occupations in graphene shows the dominant role of optical phonon scattering and the breaking of the relaxation time approximation. Our paper provides a detailed understanding of the microscopic mechanisms governing magnetotransport coefficients, establishing the BTE in a magnetic field as a broadly applicable first-principles tool to investigate transport in semiconductors and 2D materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Ab initio calculation of carrier mobility in semiconductors including ionized-impurity scattering

The past decade has seen the emergence of ab initio computational methods for calculating phonon-limited carrier mobilities in semiconductors with predictive accuracy. More realistic calculations ought to take into account additional scattering mechanisms such as, for example, impurity and grain-boundary scattering. Here, in this paper, we investigate the effect of ionized-impurity scattering on the carrier mobility. We model the analytical impurity potential parameterized from first principles by a collection of randomly distributed Coulomb scattering centers, and we include this relaxation channel into the ab initio Boltzmann transport equation, as implemented in the EPW code. We demonstrate this methodology by considering silicon, silicon carbide, and gallium phosphide, for which detailed experimental data are available. Our calculations agree well with experiments over a broad range of temperatures and impurity concentrations. For each compound investigated here, we compare the relative importance of electron-phonon scattering and ionized-impurity scattering, and we critically assess the reliability of Matthiessen's rule. We also show that an accurate description of dielectric screening and carrier effective masses can improve quantitative agreement with experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Multiscale evolution of charmed particles in a nuclear medium

Parton energy-momentum exchange with the quark gluon plasma (QGP) is a multiscale problem. In this work, we calculate the interaction of charm quarks with the QGP within the higher twist formalism at high virtuality and high energy using the Modular All Twist Transverse-scattering Elastic-drag and Radiation (MATTER) model, while the low-virtuality and high-energy portion is treated via a linearized Boltzmann transport formalism. Coherence effect that reduces the medium-induced emission rate in the MATTER model is also taken into account through a virtuality-dependent qˆ, leaving the simultaneous dependence of qˆ on heavy quark mass and virtuality for future studies. The interplay between these two formalisms is studied phenomenologically and used to produce a first description of the D-meson and charged hadron nuclear modification factor R AA across multiple centralities. As a result, all calculations were carried out utilizing the Jet Energy-loss Tomography with a Statistically and Computationally Advanced Program Envelope framework.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Asymmetric jet shapes with two-dimensional jet tomography

Two-dimensional (2D) jet tomography is a promising tool to study jet medium modification in high-energy heavy-ion collisions. It combines gradient (transverse) and longitudinal jet tomography for selection of events with localized initial jet production positions. It exploits the transverse asymmetry and energy loss that depend, respectively, on the transverse gradient and jet path length inside the quark-gluon plasma (QGP). In this study, we employ 2D jet tomography to study medium modification of the jet shape of γ -triggered jets within the linear Boltzmann transport (LBT) model for jet propagation in heavy-ion collisions. Our results show that jets with small transverse asymmetry ( A N n ⃗ ) or small γ -jet asymmetry ( x J γ = p T jet / p T γ ) exhibit a broader jet shape than those with larger A N n ⃗ or x J γ , since the former are produced at the center and go through longer path lengths while the later are off center and close to the surface of the QGP fireball. In events with finite values of A N n ⃗ , jet shapes are asymmetric with respect to the event plane. Hard partons at the core of the jet are deflected away from the denser region while soft partons from the medium response at large angles flow toward the denser part of QGP. Future experimental measurements of these asymmetric features of the jet shape can be used to study the transport properties of jets and medium responses. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

New metric improving Bayesian calibration of a multistage approach studying hadron and inclusive jet suppression

We study parton energy-momentum exchange with the quark gluon plasma (QGP) within a multistage approach composed of in-medium Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution at high virtuality, and (linearized) Boltzmann transport formalism at lower virtuality. This multistage simulation is then calibrated in comparison with high- p T charged hadrons, D mesons, and the inclusive jet nuclear modification factors, using Bayesian model-to-data comparison, to extract the virtuality-dependent transverse momentum broadening transport coefficient q ̂ . To facilitate this undertaking, we develop a quantitative metric for validating the Bayesian workflow, which is used to analyze the sensitivity of various model parameters to individual observables. The usefulness of this new metric in improving Bayesian model emulation is shown to be highly beneficial for future such analyses. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Hard jet substructure in a multistage approach

We present predictions and postdictions for a wide variety of hard jet-substructure observables using a multistage model within the framework. The details of the multistage model and the various parameter choices are described in []. A novel feature of this model is the presence of two stages of jet modification: a high-virtuality phase [modeled using the modular all twist transverse-scattering elastic-drag and radiation model ()], where modified coherence effects diminish medium-induced radiation, and a lower virtuality phase [modeled using the linear Boltzmann transport model ()], where parton splits are fully resolved by the medium as they endure multiple scattering induced energy loss. Energy-loss calculations are carried out on event-by-event viscous fluid dynamic backgrounds constrained by experimental data. The uniform and consistent descriptions of multiple experimental observables demonstrate the essential role of modified coherence effects and the multistage modeling of jet evolution. Using the best choice of parameters from [], and with no further tuning, we present calculations for the medium modified jet fragmentation function, the groomed jet momentum fraction z g and angular separation r g distributions, as well as the nuclear modification factor of groomed jets. These calculations provide accurate descriptions of published data from experiments at the Large Hadron Collider. Furthermore, we provide predictions from the multistage model for future measurements at the BNL Relativistic Heavy Ion Collider. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Quarkonium polarization in medium from open quantum systems and chromomagnetic correlators

Here, by considering the Markovian condition and applying the Wigner transformation upon the diagonal spin components of the quarkonium density matrix with the semiclassical expansion, we systematically derive the Boltzmann transport equation for quarkonia with polarization dependence in the quantum optical limit. Unlike the spin-independent collision terms governed by certain chromoelectric field correlators, new gauge invariant correlators of chromomagnetic fields determine the recombination and dissociation terms with polarization dependence at the order we are working. We also derive a Lindblad equation describing the in-medium transitions between spin-singlet and spin-triplet heavy quark-antiquark pairs in the quantum Brownian motion limit. The Lindblad equation is governed by new transport coefficients defined in terms of the chromomagnetic field correlators. Our formalism is generic and valid for both weakly coupled and strongly coupled quark gluon plasmas. It can be further applied to study spin alignment of vector quarkonia in heavy ion collisions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗