Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boltzmann Theory”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The elphbolt ab initio solver for the coupled electron-phonon Boltzmann transport equations

elphbolt is a modern Fortran (2018 standard) code for efficiently solving the coupled electron–phonon Boltzmann transport equations from first principles. Using results from density functional and density functional perturbation theory as inputs, it can calculate the effect of the non-equilibrium phonons on the electronic transport (phonon drag) and non-equilibrium electrons on the phononic transport (electron drag) in a fully self-consistent manner and obeying the constraints mandated by thermodynamics. It can calculate the lattice, charge, and thermoelectric transport coefficients for the temperature gradient and electric fields, and the effect of the mutual electron–phonon drag on these transport properties. The code fully exploits the symmetries of the crystal and the transport-active window to allow the sampling of extremely fine electron and phonon wave vector meshes required for accurately capturing the drag phenomena. The corray feature of modern Fortran, which offers native and convenient support for parallelization, is utilized. The code is compact, readable, well-documented, and extensible by design.

electron drag↗
DOE OSTI2022

First principles thermal transport modeling in GaN and related materials

Gallium nitride is a wide bandgap material utilized in a variety of technologies, including high-power electronics and light-emitting diodes, partly due to its favorable thermal properties. This chapter describes modern first-principles-based modeling of phonons and lattice thermal conductivity (k) of GaN, III-nitrides and related materials. In particular, we describe the theoretical underpinnings of calculating phonon dispersions, intrinsic phonon interactions, and other lattice dynamical properties from quantum perturbation theory and density functional theory (DFT) methods. Description of how these methods are then coupled with the Peierls-Boltzmann transport (PBT) equation to determine phonon distributions and lifetimes relevant for thermal transport is given. These theoretical and numerical methods have demonstrated quantitative accuracy and predictive power for calculating pristine k and defect-limited k from first principles for a variety of materials. We present a review of the literature utilizing DFT-PBT methods to understand novel k behaviors in III-nitrides and related materials.

Lindsay, Lucas↗
DOE OSTI2023

Electron–phonon physics from first principles using the EPW code

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

36 MATERIALS SCIENCE↗
DOE OSTI2025

ForestFlow: predicting the Lyman-α forest clustering from linear to nonlinear scales

On large scales, the Lyman-α forest provides insights into the expansion history of the Universe, while on small scales, it imposes strict constraints on the growth history, the nature of dark matter, and the sum of neutrino masses. This work introduces ForestFlow, a novel framework that bridges the gap between large- and small-scale analyses, which have traditionally relied on distinct modeling approaches. Using conditional normalizing flows, ForestFlow predicts the two Lyman-α linear biases (b δ and b η ) and six parameters describing small-scale deviations of the three-dimensional flux power spectrum (P 3D ) from linear theory as a function of cosmology and intergalactic medium physics. These are then combined with a Boltzmann solver to make consistent predictions, from arbitrarily large scales down to the nonlinear regime, for P 3D and any other statistics derived from it. Trained on a suite of 30 fixed-and-paired cosmological hydrodynamical simulations spanning redshifts from z = 2 to 4.5, ForestFlow achieves 3 and 1.5% precision in describing P 3D and the one-dimensional flux power spectrum (P 1D ) from linear scales to k = 5 Mpc −1 and k ∥ = 4 Mpc −1 , respectively. Thanks to its conditional parameterization, ForestFlow shows similar performance for ionization histories and two ΛCDM model extensions – massive neutrinos and curvature – even though none of these are included in the training set. This framework will enable full-scale cosmological analyses of Lyman-α forest measurements from the DESI survey.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Thermodynamics of the Inner Heliosheath

We derive annual sky maps of the proton temperature in the inner heliosheath (IHS), and track their temporal evolution over the years 2009–2016 of Interstellar Boundary Explorer observations. Other associated thermodynamic parameters also determined are the density, kappa (the parameter that characterizes kappa distributions), temperature rate, polytropic index, and entropy. We exploit the theory of kappa distributions and their connection with polytropes, to (i) express a new polytropic quantity Π that remains invariant along streamlines where temperature and density may vary, (ii) parameterize the proton flux in terms of the Π invariant and kappa, and (iii) derive the temperature and density, respectively, from the slope and intercept of the linear relationship between kappa and logarithm of Π. We find the following thermodynamic characteristics: (1) temperature sky maps and histograms shifted to their lowest values in 2012 and their highest in 2015; (2) temperature negatively correlated with density, reflecting the subisothermal polytropic behavior; (3) temperature positively correlated with kappa, revealing characteristics of the mechanism responsible for generating kappa distributions; (4) processes in IHS are subisothermal tending toward isobaric, consistent with previously published results; (5) linear relationship between kappa and polytropic indices, revealing characteristics of the particle potential energy; and (6) entropy positively correlated with polytropic index, aligned with the underlying theory that entropy increases toward the isothermal state where the kappa distribution reduces to the Maxwell–Boltzmann description.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Linear source method of characteristics in HELIOS2

The Studsvik Scandpower code HELIOS2 is a fast lattice transport code that has gained wide use for reactor fuel cycle analysis in commercial, regulatory, and research environments. It contains two transport solvers - Collision Probabilities (CP) and Method of Characteristics (MoC) - with a general geometry ability and generalized depletion for a variety of fuel types, experimental facilities, and isotope production applications. The HELIOS2 lattice code has been extended with the development of a linear source MoC option. The linear source equations are derived within the scope of the HELIOS2 general geometry and non-cyclic ray tracing. The implementation was tested for fine and coarse spatial mesh models of the 2D C5G7 MOX computational benchmark, which provides a pure and extreme test of a transport solver. The linear source yields at least comparable accuracy with reduced computational resource requirements relative to the flat-source approximation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Microscopic mechanism of unusual lattice thermal transport in TlInTe 2

We investigate the microscopic mechanism of ultralow lattice thermal conductivity (κ l ) of TlInTe 2 and its weak temperature dependence using a unified theory of lattice heat transport, that considers contributions arising from the particle-like propagation as well as wave-like tunneling of phonons. While we use the Peierls–Boltzmann transport equation (PBTE) to calculate the particle-like contributions (κ l (PBTE)), we explicitly calculate the off-diagonal (OD) components of the heat-flux operator within a first-principles density functional theory framework to determine the contributions (κ l (OD)) arising from the wave-like tunneling of phonons. At each temperature, T, we anharmonically renormalize the phonon frequencies using the self-consistent phonon theory including quartic anharmonicity, and utilize them to calculate κ l (PBTE) and κ l (OD). With the combined inclusion of κ l (PBTE), κ l (OD), and additional grain-boundary scatterings, our calculations successfully reproduce the experimental results. Our analysis shows that large quartic anharmonicity of TlInTe 2 (a) strongly hardens the low-energy phonon branches, (b) diminishes the three-phonon scattering processes at finite T, and (c) recovers the weaker than T –1 decay of the measured κ l .

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

Low-temperature breakdown of many-body perturbation theory for thermodynamics

It is shown analytically and numerically that the finite-temperature many-body perturbation theory in the grand canonical ensemble has zero radius of convergence at zero temperature when the energy ordering or degree of degeneracy for the ground state changes with the perturbation strength. When the degeneracy of the reference state is partially or fully lifted at the first-order Hirschfelder-Certain degenerate perturbation theory, the grand potential and internal energy diverge as T → 0. Contrary to earlier suggestions of renormalizability by the chemical potential μ, this nonconvergence, first suspected by W. Kohn and J. M. Luttinger, is caused by the nonanalytic nature of the Boltzmann factor e –E/k e T at T = 0, also plaguing the canonical ensemble, which does not involve chemical potential. The finding reveals a fundamental flaw in perturbation theory, which is deeply rooted in the mathematical limitation of power-series expansions and is unlikely to be removed within its framework.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Nonlocal suppression of Biermann battery magnetic-field generation for arbitrary atomic numbers and magnetization

The Biermann battery term of magnetohydrodynamics (MHD) generates a magnetic field where electron density gradients and electron temperature gradients are perpendicular to one another. Kinetic simulations and experiments have shown that the rate of magnetic-field generation is lower than Biermann when the electron mean free path becomes comparable to or greater than the temperature gradient scale length, known as the nonlocal regime. We investigate the nonlocal suppression of the Biermann term using simplified Fokker–Planck simulations covering a wide range of parameters. We provide the first fit for nonlocal Biermann suppression that has physically accurate behavior for small and large values of a suitable nonlocality parameter, valid for an arbitrary atomic number, and that includes the effect of magnetization on nonlocality. Lastly, the fit is intended to provide an approximate method to account for reduced magnetic-field generation in MHD codes and theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Impact of classical statistics on thermal conductivity predictions of BAs and diamond using machine learning molecular dynamics

Machine learning interatomic potentials (MLIPs) have greatly enhanced molecular dynamics (MD) simulations, achieving near-first-principles accuracy in thermal conductivity studies. In this work, we reveal that this accuracy, observed in BAs and diamond at sub-Debye temperatures, stems from an accidental error cancelation: classical statistics overestimates specific heat while underestimating phonon lifetimes, balancing out in thermal conductivity predictions. However, this balance is disrupted when isotopes are introduced, leading MLIP-based MD to significantly underpredict thermal conductivity compared to experiments and quantum statistics-based Boltzmann transport equation. This discrepancy arises not from classical statistics affecting phonon–isotope scattering rates but from its impact on the interplay between phonon–isotope and phonon–phonon scattering in the normal scattering-dominated BAs and diamond. In conclusion, this work underscores the limitations of MLIP-based MD for thermal conductivity studies at sub-Debye temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A new approach to include electron interaction effects in super transition array opacity theory

A method is presented for the improved calculation of super-shell partition functions which include the repulsive electron–electron interaction energy terms in the Boltzmann factor. Heretofore these interaction terms were approximately treated via the use of Feynman–Jensen inequalities. Further, such investigations are of particular interest for the super-transition-array approach of hot-plasma radiative opacity.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Reexamination of hydrodynamic phonon transport in thin graphite

The recent hydrodynamic phonon transport theory for graphitic materials has been supported by the measurements of the second sound at temperatures up to about 100 K. When boundary scattering becomes comparable to momentum-conserving normal phonon scattering processes that are responsible for phonon hydrodynamics, Poiseuille phonon flow phenomena can emerge to give rise to unique size-dependent thermal conductivity in thin graphite. Here, we examine the thickness range for the Poiseuille phonon flow to become observable in thin graphite with the use of both deviational Monte Carlo simulation of the Peierls-Boltzmann transport equation and four-probe thermal transport measurements. As the basal-plane thermal conductivity calculated by prior first-principles theories saturates to the graphite value when the thickness is increased to five graphene layers, the phonon dispersion of graphite is used in the current calculations of thin graphite of micrometer thickness and a 23-layer thick ultrathin graphite (UTG) sample. The calculations show that diffuse surface scattering by surface defects can lead to Poiseuille phonon flow at 50 K in thin graphite with the thickness close to several micrometers but not in the 65 μm thin graphite and 23-layer UTG, where phonon scattering by the top and bottom surfaces become, respectively, much less and more frequent than the normal processes. In addition, the calculation results with the bulk graphite dispersion and diffuse surface scattering show decreased basal-plane thermal conductivity with decreasing thickness, opposite to recent thermocouple measurements of thin graphite samples. In comparison, the calculation results reveal that partially diffuse surface defect scattering can yield the four-probe measurement results of UTG samples, which are prepared here with an improved process to minimize surface contamination.

36 MATERIALS SCIENCE↗
DOE OSTI2020

(Photo)Electrocatalytic CO 2 Reduction at the Defective Anatase TiO 2 (101) Surface

Excessive carbon dioxide (CO 2 ) emissions by combustion of fossil fuels are linked to global warming and rapid climate change. One promising route to lowering the concentration of CO 2 in the atmosphere is to reduce it to useful small molecules via photoelectrocatalytic hydrogenation, which would enable solar energy storage with a zero-carbon emission cycle and perform a more efficient separation of the photogenerated electron and hole pair than pure photocatalysis. Indeed, photoelectrocatalytic CO 2 reduction has been an intense focus of research. Using the density functional theory (DFT), we studied the CO 2 reduction reaction on the defective anatase TiO 2 (101) surface, at both the solvent/catalyst and the electrolyte/catalyst interfaces. The analysis of the electronic structure of the surface shows a contrast between the solvent/catalyst and the electrolyte/catalyst interfaces, which results in the two corresponding catalytic cycles being distinct. Furthermore, our study explains at the electronic and mechanistic levels why methanol is the main product in the presence of the electrolyte and why the overpotential is not only controlled by the reaction process but also by the diffusion process.

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

QCD axion: Some like it hot

We compare the quantum chromodynamics (QCD) axion phase-space distribution from unitarized next-to-leading order chiral perturbation theory with the one extracted from pion-scattering data. We derive a robust bound by confronting momentum-dependent Boltzmann equations against up-to-date observations of the cosmic microwave background, of the baryonic acoustic oscillations and of primordial abundances. These datasets imply 𝑚 𝑎 ≤ 0.16 eV for the 95% credible interval, i.e., ∼ 30% stronger bound than what previously found. We present forecasts using dedicated likelihoods for future cosmological surveys and the sphaleron rate from unquenched lattice QCD.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A new approach for resonance treatment of doubly heterogeneous fuel using the RSE method

A new resonance calculation method for the doubly-heterogeneous (DH) fuels such as high- temperature gas-cooled reactor fuel is proposed based on the resonance calculation based on Spectral Expansion (RSE) method. The concept of pointwise disadvantage factor for fuel grain is taken into account to treat the DH fuels. The verification calculation is carried out for simplified single fuel cell and fuel compact consisting of five fuel cells and graphite moderator. The calculation results indicate that the present method can appropriately handle the space- dependent self-shielding effect for DH fuels. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Performance enhancement of planar MOC calculation in hexagonal geometries by assembly-wise domain decomposition

The parallelization of planar Method of Characteristics (MOC) calculation with the higher order scattering in the hexagonal version of the nTRACER direct whole core calculation code is enhanced by adopting assembly-wise domain decomposition and other techniques. In order to augment the hexagonal domain decomposition, the three-color scheme, the angular decomposition, and the angular flux storage scheme are examined. The parallel performance assessed for the VVER benchmark problems demonstrates that the hexagonal assembly-wise domain decomposition significantly reduces the computing time in ray tracing by up to 64 %. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Variance reduction techniques for Monte Carlo neutron noise simulations

The small fluctuations of the neutron flux caused by small perturbations of the macroscopic cross-sections take the name of neutron noise. Advanced Monte Carlo methods have been recently proposed in order to solve the neutron noise equations in the frequency domain, which allows establishing reference solutions to validate faster but approximate deterministic solvers. Due to the presence of particles carrying two statistical weights (for the real and imaginary components of the noise field), both of which may be positive or negative, Monte-Carlo simulations of neutron noise pose distinct challenges in terms of variance reduction. In this work we investigate two variance-reduction techniques, namely branchless collisions and weight cancellation, and probe their effectiveness for a benchmark con- figuration concerning the noise field induced by a pin with oscillating cross sections in a fuel assembly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗