Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boltzmann Transport Equation”

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

Optical phonon dominated heat transport: A first-principles thermal conductivity study of BaSn 2

Acoustic phonons with long mean free paths have long been believed to control the lattice thermal conductivity κ L in solids dominantly. In this study, however, we demonstrate an optical phonon dominated κ L in BaSnS 2 . By solving the Peierls-Boltzmann transport equation, we predict a low diagonal lattice thermal conductivity κ L (D) of 0.34 W m -1 K -1 at 850 K, which is less than half the κ L (D) of SnS at the same temperature. Further calculations following the Allen-Feldman model suggest the additional off-diagonal lattice thermal conductivity κ L (OD) contributed by wavelike tunneling phonons. The κ L (OD) becomes pronounced at the high temperature ( 0.17 W m -1 K -1 at 850 K) and leads to a deviation of the temperature dependence of κ L from T -1 to T -0.76 , suggesting the potential lattice anharmonicity in BaSnS 2 . Further analyses indicate BaSnS 2 has over 68% of κ L contributed by optical phonons. We show this uncommon optical phonon dominated κ L is due to the relatively high group velocities of optical phonons in BaSnS 2 . The phonon mode visualization suggests these relatively high-velocity optical phonons correspond to the antiphase vibrations in BaSnS 2 monolayers, which is originated from the unique permutation of SnS 3 tetrahedra. Finally, by investigating the mode-resolved group velocity, relaxation time, and Grüneisen parameter, we attribute the intrinsic low κ L of BaSnS 2 to the soft lattice and the relatively high lattice anharmonicity induced by the Ba-S weak bonding and Sn(II) lone-pair electrons. Our study explicitly analyzes the microscopic mechanism of optical phonon dominated heat transport in BaSnS 2 and suggests it worthy of further experimental studies as an intrinsic low-κ L material.

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 OSTI2024

Dynamics of nonequilibrium magnons in gapped Heisenberg antiferromagnets

Nonequilibrium dynamics in spin systems is a topic currently under intense investigation as it provides fundamental insights into thermalization, universality, and exotic transport phenomena. While most of the studies have been focused on ideal closed quantum many-body systems such as ultracold atomic quantum gases and one-dimensional spin chains, driven-dissipative Bose gases in steady states away from equilibrium in classical systems also lead to intriguing nonequilibrium physics. Here, in this work, we theoretically investigate out-of-equilibrium dynamics of magnons in a gapped Heisenberg quantum antiferromagnet based on Boltzmann transport theory. We show that, by treating scattering terms beyond the relaxation-time approximation in the Boltzmann transport equation, energy and particle number conservation mandate that nonequilibrium magnons cannot relax to equilibrium, but decay to other nonequilibrium stationary states. The only decay channel for these stationary states back to equilibrium is through the nonconserving interactions (i.e., changing particle number and/or energy within the magnon system) such as boundary or magnon-phonon scattering. At low temperatures, these nonconserving interactions are much slower processes than intrinsic magnon-magnon interaction in a gapped spin system. Using magnon-phonon interaction as a quintessential type of nonconserving interaction, we then propose that nonequilibrium steady states of magnons can be maintained and tailored using periodic driving at frequencies faster than relaxation due to phonon interactions. These findings reveal a class of classical material systems that are suitable platforms to study nonequilibrium statistical physics and macroscopic phenomena such as classical Bose-Einstein condensation of quasiparticles and magnon supercurrents that are relevant for spintronic applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗
DOE OSTI2020

Particlelike Phonon Propagation Dominates Ultralow Lattice Thermal Conductivity in Crystalline Tl 3 VSe 4

We explore the microscopic mechanisms of ultralow lattice thermal conductivity ($\kappa_{l}$) in Tl$_{3}$VSe$_{4}$~by combining a first-principles density-functional theory (DFT) based framework of anharmonic lattice dynamics with the Peierls-Boltzmann transport equation (PBTE) for phonons. We include contributions of the three- and four-phonon scattering processes to the phonon lifetimes as well as the temperature-dependent anharmonic renormalization of phonon energies arising from an unusually strong quartic anharmonicity in Tl$_{3}$VSe$_{4}$. In contrast to a recent report by Mukhopadhyay~\etal [\textcolor{blue}{Science 360, 1455 (2018)}] which suggested that a significant contribution to $\kappa_{l}$ arises from random walks among uncorrelated oscillators, we show that particle-like propagation of phonon excitations can successfully explain the experimentally observed ultralow $\kappa_{l}$. Our results are also supported by explicit calculations of the off-diagonal terms of the heat-current operator, which are found to be small and indicate that wave-like tunneling of heat-carrying vibrations is of minor importance. Our results (i) resolve the discrepancy between the theoretical and experimental $\kappa_{l}$, (ii) offer new insights into the minimum $\kappa_{l}$ achievable in \TlVSe, and (iii) highlight the importance of high-order anharmonicity in low-$\kappa_{l}$ systems. The methodology demonstrated here may be used to resolve the discrepancies between the experimentally measured and the theoretically calculated $\kappa_{l}$ in skutterides and perovskites, as well as to understand the glasslike $\kappa_{l}$ in complex crystals with strong anharmonicity, leading towards the goal of rational design of new materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Thermal conductivity of the quasi-one-dimensional materials TaSe 3 and ZrTe 3

The high breakdown current densities and resilience to scaling of the metallic transition metal trichalcogenides TaSe 3 and ZrTe 3 make them of interest for possible interconnect applications, and it motivates this study of their thermal conductivities and phonon properties. These crystals consist of planes of strongly bonded one-dimensional chains more weakly bonded to neighboring chains. Phonon dispersions and the thermal conductivity tensors are calculated using density functional theory combined with an iterative solution of the phonon Boltzmann transport equation. The phonon velocities and the thermal conductivities of TaSe 3 are considerably more anisotropic than those of ZrTe 3 . The maximum LA velocity in ZrTe 3 occurs in the cross-chain direction, and this is consistent with the strong cross-chain bonding that gives rise to large Fermi velocities in that direction. The thermal conductivities are similar to those of other metallic two-dimensional transition metal dichalcogenides. At room temperature, a significant portion of the heat is carried by the optical modes. In the low frequency range, the phonon lifetimes and mean free paths in TaSe 3 are considerably shorter than those in ZrTe 3 . Here, we find that the shorter lifetimes in TaSe 3 are consistent with the presence of lower frequency optical branches and zone-folding features in the acoustic branches that arise due to the doubling of the TaSe 3 unit cell within the plane.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mesoscopic finite-size effects of unconventional electron transport in PdCoO 2

A wide range of unconventional transport phenomena has recently been observed in single-crystal delafossite metals. Here, we present a theoretical framework to elucidate electron transport using a combination of first-principles calculations and numerical modeling of the anisotropic Boltzmann transport equation. Using PdCoO 2 as a model system, we study different microscopic electron and phonon scattering mechanisms and establish the mean free path hierarchy of quasiparticles at different temperatures. We treat the anisotropic Fermi surface explicitly to numerically obtain experimentally-accessible transport observables, which bridge between the “diffusive,” “ballistic,” and “hydrodynamic” transport regime limits. Further, we illustrate that the distinction between the “quasiballistic” and “quasihydrodynamic” regimes is challenging and often needs to be quantitative in nature. From first-principles calculations, we populate the resulting transport regime plots and demonstrate how the Fermi surface orientation adds complexity to the observed transport signatures in micrometer-scale devices. Our work provides key insights into microscopic interaction mechanisms on open hexagonal Fermi surfaces and establishes their connection to the macroscopic electron transport in finite-size channels.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-principles determination of the phonon-point defect scattering and thermal transport due to fission products in ThO 2

This work presents the first-principles calculations of the lattice thermal conductivity degradation due to point defects in thorium dioxide using an iterative solution of the Pierels-Boltzmann transport equation. We have used the non perturbative Green’s function methodology to compute the phonon- point defect scattering rates that consider the local distortion around the point defect, including the mass difference changes, interatomic force constants, and structural relaxation near the point defects. The point defects considered in this work include the vacancy of thorium ($V$ $TH$ ) and oxygen ($V$ $O$ ), substitutions of helium ($He$ $Th$ ), krypton ($Kr$ $Th$ ), zirconium ($Zr$ $Th$ ), iodine ($I$ $Th$ ) and xenon ($Xe$ $Th$ ) in the thorium site, and the three different configurations of the Schottky defects. The results of the phonon-defect scattering rate reveal that, among all the considered intrinsic defects, the thorium vacancy and helium substitution in the thorium site scatter the phonon most due to the substantial changes in the force constant and structural distortions. The scattering of phonons due to the substitutional defects unveils that the zirconium atom scatters phonons the least, followed by xenon, iodine, krypton, and helium. This is contrary to the intuition that the scattering strength follows $He$ $Th$ > $Kr$ $Th$ > $Zr$ $Th$ > $I$ $Th$ > $Xe$ $Th$ based on the mass difference. This striking difference in the zirconium phonon scattering is due to the local chemical environment changes. Zirconium is an electropositive element with valency similar to thorium and, therefore, can bond with the oxygen atoms, thus creating less force constant variance compared to iodine, an electronegative element noble gas helium, xenon, and krypton. In conclusion, these results can serve as the benchmark for the analytical models and help the engineering-scale modeling effort for nuclear design.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Intrinsic Limits of Charge Carrier Mobilities in Layered Halide Perovskites

Layered halide perovskites have emerged as potential alternatives to three-dimensional (3D) halide perovskites due to their improved stability and larger material phase space, allowing fine tuning of structural, electronic, and optical properties. However, their charge carrier mobilities are significantly smaller than those of 3D halide perovskites, which has a considerable impact on their application in optoelectronic devices. Here, we employ state-of-the-art approaches to unveil the electron-phonon mechanisms responsible for the diminished transport properties of layered halide perovskites. Starting from a prototypical A M X 3 halide perovskite, we model the case of n = 1 and n = 2 layered structures and compare their electronic and transport properties to the 3D reference. The electronic and phononic properties are investigated within density functional theory (DFT) and density functional perturbation theory (DFPT), while transport properties are obtained via the Boltzmann transport equation. The vibrational modes contributing to charge carrier scattering are investigated and associated with polar-phonon scattering mechanisms arising from the long-range Fröhlich coupling and deformation-potential scattering processes. Our investigation reveals that the lower mobilities in layered systems primarily originate from the increased electronic density of states at the vicinity of the band edges, while the electron-phonon coupling strength remains similar. Such an increase is caused by the dimensionality reduction and the break in octahedra connectivity along the stacking direction. Our findings provide a fundamental understanding of the electron-phonon coupling mechanisms in layered perovskites and highlight the intrinsic limitations of the charge carrier transport in these materials. Published by the American Physical Society 2024

Cucco, Bruno (ORCID:0000000331564143)↗
DOE OSTI2023

Effect of Strain on the Resistivity and Thermal Conductivity of High Purity Niobium

High purity niobium (Nb) is a technologically important material for large-scale accelerator and nano-scale quantum computing applications in microwave frequency range. The high thermal conductivity and low resistivity of Nb are critical to the high performance at temperature range of 0.01–4.0 K. The presence of interstitials such as O, N, H, and C act as scattering centers and alter the mean free path, reducing resistivity and thermal conductivity, and contributing significantly to Nb's thermal performance for temperatures of 2.0 K and above. The residual resistivity ratio (RRR), defined as the ratio of the normal state resistivity at 300 K to that at 4.2 K (Nb, T c = 9.2 K), is an accepted direct estimate of the impurity content of fully recrystallized Nb. Complete re-crystallization of Nb is challenging unless very high temperatures are employed, which is often impractical, hence, in practice, dislocation and dislocation structures impact the thermal performance of Nb due to strong phonon scattering contributions. Furthermore, this paper reports on the degradation of thermal conductivity and RRR of high purity Nb large grain, single crystal with fixed impurity, varying strain, and dislocation content levels. Experimental thermal conductivity data fits the Boltzmann transport equation incorporating dislocation density.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thermal transport in Yb-based 1-2-20 materials

Abstract Lattice thermal-transport properties of heavy-fermion YbT 2 Zn 20 -based (T = Co, Rh, and Ir) 1-2-20 compounds are calculated with the single-mode relaxation-time approximation and the full solution of the linearized phonon Boltzmann transport equation from first-principles anharmonic phonon calculations. We predict low lattice thermal conductivity κ L with a maximum value of ∼ 24.16 · K at 85.14 K, 137.29 · K at 9.67 K, and 23.55 · K at 67.50 K and ∼ 12.02, 40.04, and 10.30 · K at room temperature for YbCo 2 Zn 20 , YbRh 2 Zn 20 , and YbIr 2 Zn 20 , respectively. Based on the analysis of the cumulative κ L as a function of the phonon mean free path and the frequency-dependent joint density of states at various temperatures, we attribute the low κ L to a rattling mode with an average characteristic rattling frequency that enhanced phonon scattering processes, which reduced the phonon mean free paths, suppressed the phonon lifetime, and enhanced the probability of three-phonon scattering events. The predicted low κ L , especially in YbCo 2 Zn 20 and YbIr 2 Zn 20 , makes them promising candidate materials for thermoelectric applications and thermal management.

Physics↗
DOE OSTI2020

(U) Transport Corrections Implemented in SENSMG

Truncating the spherical harmonics expansion of the neutron scattering source in the Boltzmann transport equation leads to truncation error. Methods for approximately correcting for this truncation error are called transport corrections. The PARTISN multigroup discrete ordinates neutron transport code has three transport correction options. These options can be specified in the SENSMG multigroup neutron sensitivity code and passed to PARTISN for the neutron transport.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

MPACT Software Management Plan (V.4.3)

The MPACT code solves a discretized form of the Boltzmann transport equation on a wide variety of geometries and is distributed with a multigroup neutron cross section library. MPACT provides an advanced geometrically resolved neutral-particle transport capability to solve the flux distribution throughout the entire problem geometry, and it can model the isotopic depletion, decay, and activation of materials. The flux solution in MPACT is provided using a 2D/1D synthesis method within the framework of the 3D coarse mesh finite difference (CMFD) method for which axial and radial correction factors are obtained from 2D method of characteristics (MOC) and 1D nodal expansion method (NEM), PN, or SN. Other key characteristics of the MPACT code include the subgroup method and the embedded self-shielding method (ESSM) for resonance treatment, depletion capability based on the ORIGEN exponential matrix method, and a simplified thermal-hydraulics method for temperature/fluid feedback. The sole purpose of the simplified feedback model is to provide a mechanism for testing during code development and to provide a limited capability for educational applications. Work performed at the code level supports the VERA-QA-001, quality assurance program plan (QAPP) and VERA-QA-002, VERA Software Quality Assurance Plan.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Transport @ LANL: An Overview [Slides]

Deterministic neutron transport requires discretization of the Boltzmann transport equation in space, angle, energy, and time. Discretization in energy is typically performed with the multigroup approximation, which tries to preserve continuous energy cross section information by using a weighted average over discrete energy groups. There are two important requirements for accurate multigroup cross sections: appropriate energy group boundaries and an accurate weight function. The former can be selected with appropriate optimizers; the latter either requires flux solutions from comparable simulations or expert knowledge.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Semiconductor Nanostructures by Scientific Design (Final Technical Report)

We developed and used quantum simulation techniques to predict structural, electronic, optical and thermal properties affecting the scientific design of semiconductor nanostructures. We focused on the study of surfaces and interfaces at the nanoscale, in particular embedding and solvation of semiconductor nanoparticles, and provided a microscopic description of nanostructured materials in realistic environments, directly comparable with experimental conditions. We addressed three major problems: (i) the design of realistic structural models of integrated nanostructures, e.g. quantum dots embedded in solid matrices and nanostructures immersed in fluids; (ii) the development of efficient and accurate calculations to describe the excited state properties of nanostructures, beyond standard Density Functional Theory, so as to understand and predict electronic properties involved in renewable energy applications, for example solar conversion devices. In particular we developed and used methods based on a many body perturbation theory. (iii) the development of calculations of heat transport properties of nanostructured materials, for thermoelectric applications, including predictive calculations of thermal and electrical conductivity and of Seebeck coefficients. Specifically, we developed and used methods based on the Boltzmann transport equation and molecular dynamics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

MPACT 4.4 Theory Manual

MPACT is a three-dimensional (3D) full-core neutron transport code capable of calculating subpin power distributions. Calculations are based on the Boltzmann transport equation for neutron fluxes for problems in which the detailed geometrical configuration of fuel components such as the pellet and cladding are explicitly retained. The cross-section data needed for the neutron transport calculation are obtained directly from a multigroup cross section library, which has traditionally been used by lattice physics codes to generate few-group homogenized cross sections for nodal core simulators. Hence, MPACT assumes neither a priori homogenization nor group condensation for the full core spatial solution. The 3D MPACT transport solution can be obtained using the method of characteristics (MOC), which employs discrete ray tracing within each fuel pin. However, for practical reactor applications, the direct application of MOC to 3D core configurations requires an excessive amount of memory and computing time due to the very large number of rays. For practical 3D full-core calculations, MPACT commonly uses an approximate “2D/1D” method that treats the radial (x and y) variables differently from the axial (z) variable. In particular, the radial dependence of the solution is calculated using transport theory, and the axial dependence is calculated using diffusion or P 3 theory. The 2D/1D method requires the core to be divided into a vertical stack of axial slices with a thickness of Δ z ≈ 5–10 cm. Each axial slice is divided radially into coarse spatial cells with boundaries that usually constitute the pin cell boundaries, for which Δ x = Δ y ≈ 1.5 cm. Then, each coarse radial cell (pin cell) is divided into 50–100 fine radial cells, which resolve the angular flux in the fuel, cladding, and moderator regions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

SCEPTRE 2.7 User's Guide

Sandia’s Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) is a computer code that solves the linear Boltzmann transport equation, particularly targeting coupled photon-electron problems. It uses unstructured finite element meshes in space, multigroup in energy, and discrete ordinates (Sn) or other methods in angle. SCEPTRE uses an xml-based input file to specify the problem. This report documents the options and syntax of that input file.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

SCEPTRE 2.6 User's Guide

Sandia’s Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) is a computer code that solves the linear Boltzmann transport equation, particularly targeting coupled photon-electron problems. It uses unstructured finite element meshes in space, multigroup in energy, and discrete ordinates (S n ) or other methods in angle. SCEPTRE uses an xml-based input file to specify the problem. This report documents the options and syntax of that input file.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗