Search NASA⌕ Search

SEARCH · Search NASA

Results for “FERMI-DIRAC STATISTICS”

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.

16 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

thornado-transport: IMEX schemes for two-moment neutrino transport respecting Fermi-Dirac statistics

We develop implicit-explicit (IMEX) schemes for neutrino transport in a background material in the context of a two-moment model that evolves the angular moments of a neutrino phase-space distribution function. Considering the upper and lower bounds that are introduced by Pauli's exclusion principle on the moments, an algebraic moment closure based on Fermi-Dirac statistics and a convex-invariant time integrator both are demanded. A finite-volume/first-order discontinuous Galerkin(DG) method is used to illustrate how an algebraic moment closure based on Fermi-Dirac statistics is needed to satisfy the bounds. Several algebraic closures are compared with these bounds in mind, and the Cernohorsky and Bludman closure, which satisfies the bounds, is chosen for our IMEX schemes. For the convex-invariant time integrator, two IMEX schemes named PD-ARS have been proposed. PD-ARS denotes a convex-invariant IMEX Runge-Kutta scheme that is high-order accurate in the streaming limit, and works well in the diffusion limit. Our two PD-ARS schemes use second-and third-order, explicit, strong-stability-preserving Runge-Kutta methods as their explicit part, respectively, and therefore are second-and third-order accurate in the streaming limit, respectively. The accuracy and convex-invariance of our PD-ARS schemes are demonstrated in the numerical tests with a third-order DG method for spatial discretization and a simple Lax-Friedrichs flux. The method has been implemented in our high-order neutrino-radiation hydrodynamics (thornado) toolkit. We show preliminary results employing tabulated neutrino opacities.

Chu, Ran↗
DOE OSTI2019

Realizability-preserving DG-IMEX method for the two-moment model of fermion transport

Building on the framework of Zhang & Shu [1], [2], we develop a realizability-preserving method to simulate the transport of particles (fermions) through a background material using a two-moment model that evolves the angular moments of a phase space distribution function f . The two-moment model is closed using algebraic moment closures; e.g., as proposed by Cernohorsky & Bludman [3] and Banach & Larecki [4]. Variations of this model have recently been used to simulate neutrino transport in nuclear astrophysics applications, including core-collapse supernovae and compact binary mergers. We employ the discontinuous Galerkin (DG) method for spatial discretization (in part to capture the asymptotic diffusion limit of the model) combined with implicit-explicit (IMEX) time integration to stably bypass short timescales induced by frequent interactions between particles and the background. Appropriate care is taken to ensure the method preserves strict algebraic bounds on the evolved moments (particle density and flux) as dictated by Pauli's exclusion principle, which demands a bounded distribution function (i.e., f ϵ [0,1]). This realizability-preserving scheme combines a suitable CFL condition, a realizability-enforcing limiter, a closure procedure based on Fermi-Dirac statistics, and an IMEX scheme whose stages can be written as a convex combination of forward Euler steps combined with a backward Euler step. The IMEX scheme is formally only first-order accurate, but works well in the diffusion limit, and — without interactions with the background — reduces to the optimal second-order strong stability-preserving explicit Runge-Kutta scheme of Shu & Osher [5]. Numerical results demonstrate the realizability-preserving properties of the scheme. As a result, we also demonstrate that the use of algebraic moment closures not based on Fermi-Dirac statistics can lead to unphysical moments in the context of fermion transport.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Modeling the Photoelectrochemical Evolution of Lead-Based, Mixed-Halide Perovskites Due to Photosegregation

Lead-based, mixed-halide perovskites such as methylammonium lead iodide-bromide [MAPb(I 1-x Br x ) 3 ] undergo anion photosegregation under illumination. This is observed as low band gap photoluminescence from photogenerated iodine-rich domains due to favorable band offsets that induce carrier funneling into them. Unfortunately, theoretical rationalizations of mixed-halide photosegregation are complicated by biases inherent to photoluminescence-based observations. Recent compositionally-weighted X-ray diffraction (XRD) measurements now reveal broad distributions of photosegregated stoichiometries not captured by existing photosegregation models. To better bridge experiment and theory, we perform kinetic Monte Carlo (KMC) simulations of photosegregation within the context of a band gap-based thermodynamic model, which has previously accounted for numerous experimental observations. Our KMC simulations are modified to consider high carrier density Fermi-Dirac statistics that result from carrier funneling and accumulation within photosegregated I-rich domains. Obtained KMC results reproduce broad x terminal distributions seen experimentally and illustrate how they are characterized by a central, heavily I-enriched stoichiometry. I-rich domain “drifting” during photosegregation rationalizes the long photosegregation timescales seen experimentally with drifting simultaneously producing a wake of variable stoichiometry I-rich inclusions that form the lion’s share of stoichiometric heterogeneities seen in compositionally-weighted XRD measurements. These simulations and accompanying rationalizations further reveal a general criterion for realizing favorable free energies to induce demixing. Central to the criterion is the statistical occupation of low gap inclusions in the parent alloy by excitations. Furthermore, the resulting model thus provides a general framework for conceptualizing mixed-halide perovskite light and temperature sensitivities, mediated by photocarriers.

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

Electrochemomechanics with flexoelectricity and modelling of electrochemical strain microscopy in mixed ionic-electronic conductors

Recently, a new scanning probe microscopy approach, referred to as electrochemical strain microscopy (ESM), for probing local ionic flows and electrochemical reactions in solids based on the bias-strain coupling was proposed by Morozovska et al. Then, a series of theoretical papers for analyzing the image formation and spectroscopic mechanism of ESM were published within the framework of Fermi-Dirac statistics, the Vegard law, the direct flexoelectric coupling effect, the electrostriction effect, and so on. However, most of the models in these papers are limited to the partial coupling or particular process, and numerically solved by using decoupling approximation. In this paper, to model the ESM measurement with the coupling electrical-chemical-mechanical process, the chemical Gibbs function variational principle for the thermal electrical chemical mechanical fully coupling problem is proposed. The fully coupling governing equations are derived from the variational principle. When the tip concentrates the electric field within a small volume of the material, the inhomogeneous electric field is induced. So, both direct and inverse flexoelectric effects should be taken into account. Here, the bulk defect electrochemical reactions are also taken into account, which are usually omitted in the existing works. This theory can be used to deal with coupling problems in solids, including conductors, semiconductors, and piezoelectric and non-piezoelectric dielectrics. As an application of this work, a developed initial-boundary value problem is solved numerically in a mixed ion-electronic conductor. Numerical results show that it is meaningful and necessary to consider the bulk defect chemical reaction. Besides, the chemical reaction and the flexoelectric effect have an interactive influence on each other. This work can provide theoretical basis for the ESM as well as investigating the bulk chemical reaction process in solids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Efficient numerical method for evaluating normal and anomalous time-domain equilibrium Green's functions in inhomogeneous systems

In this work we develop EPOCH (equilibrium propagator by orthogonal polynomial chain), a computationally efficient method to calculate the time-dependent equilibrium Green's functions, including the anomalous Green's functions of superconductors, to capture the time evolution in large inhomogeneous systems. The EPOCH method generalizes the Chebyshev wave-packet propagation method from quantum chemistry and efficiently incorporates the Fermi-Dirac statistics that is needed for equilibrium quantum condensed matter systems. The computational cost of EPOCH scales only linearly in the system degrees of freedom, generating an extremely efficient algorithm also for very large systems. We demonstrate the power of the EPOCH method by calculating the time evolution of an excitation near a superconductor–normal metal interface in two and three dimensions, capturing transmission as well as normal and Andreev reflections.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Comparison of heavy-ion transport simulations: Collision integral in a box

Simulations by transport codes are indispensable to extract valuable physical information from heavy-ion collisions. In order to understand the origins of discrepancies among different widely used transport codes, we compare 15 such codes under controlled conditions of a system confined to a box with periodic boundary, initialized with Fermi-Dirac distributions at saturation density and temperatures of either 0 or 5 MeV. In such calculations, one is able to check separately the different ingredients of a transport code. In this second publication of the code evaluation project, we only consider the two-body collision term; i.e., we perform cascade calculations. When the Pauli blocking is artificially suppressed, the collision rates are found to be consistent for most codes (to within 1 % or better) with analytical results, or completely controlled results of a basic cascade code. In orderto reach that goal, it was necessary to eliminate correlations within the same pair of colliding particles that can be present depending on the adopted collision prescription. In calculations with active Pauli blocking, the blocking probability was found to deviate from the expected reference values. The reason is found in substantial phase-space fluctuations and smearing tied to numerical algorithms and model assumptions in the representation of phase space. This results in the reduction of the blocking probability in most transport codes, so that the simulated system gradually evolves away from the Fermi-Dirac toward a Boltzmann distribution. Since the numerical fluctuations are weaker in the Boltzmann-Uehling-Uhlenbeck codes, the Fermi-Dirac statistics is maintained there for a longer time than in the quantum molecular dynamics codes. As a result of this investigation, we are able to make judgements about the most effective strategies in transport simulations for determining the collision probabilities and the Pauli blocking. Finally, investigation in a similar vein of other ingredients in transport calculations, like the mean-field propagation or the production of nucleon resonances and mesons, will be discussed in the future publications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Entanglement suppression for $ΩΩ$ scattering

We study entanglement suppression in $s$-wave $ΩΩ$ scattering, where each baryon has spin $3/2$. By treating the $S$-matrix as a quantum operator acting on the spin states, we quantify its ability to generate entanglement and identify the conditions on the phase shifts of the spin channels that minimize entanglement generation in the system. In $ΩΩ$ scattering, only antisymmetric spin channels are allowed due to Fermi-Dirac statistics. Applying the entanglement-suppression framework to $ΩΩ$ scattering, we find two solutions for the phase shifts: one leading to a spin SU(4) symmetry and the other to a nonrelativistic conformal symmetry. We show that the solution associated with the nonrelativistic conformal symmetry originates from the specific structure of the Clebsch-Gordan coefficients in the $3/2 \otimes 3/2$ system.

Sone, Katsuyoshi [Tokyo Metropolitan U.] (ORCID:00↗
DOE OSTI2018

Metamorphic narrow-gap InSb/InAsSb superlattices with ultra-thin layers

Metamorphic strain compensated InSb/InAsSb 0.52 superlattices (SLs) with ultrathin layers and different periods grown on GaSb substrate were designed, fabricated, and characterized. It was shown that a period increase from 3 to 6.2 nm reduced the effective bandgap energy from 70 to 0 meV. A further increase in the period leads to inversion of the valence and conduction bands. Magneto-optical experiments demonstrated that Dirac-like carrier dispersion is characteristic of almost gapless InSb/InAsSb 0.52 SLs. Indication of hole transport enhancement over that found in InAsSb/InAsSb SL structures is presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Nonlinear coupling of electromagnetic and electron acoustic waves in multi-species degenerate astrophysical plasma

Nonlinear wave-coupling is studied in a multispecies degenerate astrophysical plasma consisting of two electron species (at different temperatures): a highly degenerate main component plus a smaller classical relativistic flow immersed in a static neutralizing ion background. It is shown that the high frequency electromagnetic waves through their strong nonlinear interactions with the electron-acoustic waves [sustained by a multielectron component (degenerate) plasma surrounding a compact astrophysical object] can scatter to lower frequencies so that the radiation observed faraway will be spectrally shifted downward. It is also shown that, under definite conditions, the electromagnetic waves could settle into stationary solitonic states. It is expected that the effects of such structures may persist as detectable signatures in forms of modulated micropulses in the radiation observed far away from the accreting compact object. Both these effects will advance our abilities to interpret the radiation coming out of the compact objects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

On the theory of nonhomogeneous nonequilibrium superconductivity in 2D systems with massless fermions

Here we analyze static and nonequilibrium superconducting properties of a 2D relativistic-like model system with local electron-electron interaction, Rashba spin-orbit interaction αR in presence of time-dependent in-plane magnetic field H(t). It is shown that similar to the 2D case with ordinary massive quasiparticle dispersion ε(k)~|k| 2 at large fields, such a system demonstrates a nonhomogeneous superconducting stripe phase with the order parameter Δ(r)=Δ(0)cos(2[μ B B×r]n/$\hbarυ_F$) (B is the magnetic induction, υF is the Fermi velocity, n is the normal to the plane, μ B is the Bohr magneton, and αR$\ll$υF) where the stripes are oriented along the B direction. In the considered system, the inter-stripe period L and the magnitude of the magnetic field B are related by a universal relation BL=$\hbarυ_F$/μ B ≃0.714∙10 -4 Tm. Contrary to the case of massive quasiparticles, where the condition α R ~υ F can be, in principle, satisfied by increasing α R or by charge doping (Fermi velocity decreasing), in a relativistic-like system, where υF is doping-independent and one-two orders of magnitude larger than typical Fermi velocity in the “standard” 2D systems, the stripe phase can be the ground state at a rather low doping level. We also analyzed the nonequilibrium properties of the system with a focus on the melting of the stripe order (when the magnetic field is quenched to a lower value) and stripe dynamics (when the field is rotated by 90° degrees) and found several notable results. In particular, it was shown that the stripe domains melt according to law R~1$\sqrt{t}$ at initial times, while at longer times they shrink exponentially. In the case of the flipped magnetic field, the stripe orientation gradually turns from x- to y-direction, and the intermediate “crossed-stripe” phase takes place during times of order of picoseconds. Such a crossed phase is built of periodic superconducting bubbles that potentially may have applications in modern ultrafast superconducting technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Transient work function gating: A new photoemission regime

We present the theoretical basis for a new photoemission regime, transient work function gating (TWFG), which temporally and energetically gates photoemission and produces near-threshold photoelectrons with thermally limited emittance, percent-level quantum efficiency, and control over temporal coherence. The technique consists of actively gating the work function of a generalized photocathode using a non-ionizing long-wavelength optical field to produce an adiabatic modulation of the carrier density at their surface. We examine TWFG as a means to circumvent the long-standing trade-off between low emittance and high quantum efficiency, untethered to particle source or photocathode specifics. TWFG promises new opportunities in photoemission physics for next generation electron and accelerator-based x-ray photon sources.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Piecewise interaction picture density matrix quantum Monte Carlo

The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact N-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC (IP-DMQMC) method that overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation, thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.

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

Thermal mean-field theories

Several closely related ab initio thermal mean-field theories for fermions, both well-established and new ones, are compared with one another at the formalism level and numerically. The theories considered are Fermi–Dirac theory; thermal Hartree–Fock (HF) theory; two modifications of the thermal single-determinant and the first-order finite-temperature many-body perturbation theory based on a zero-temperature or thermal HF reference. Furthermore, thermal full-configuration-interaction theory is used as the benchmark.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Damped Dirac magnon in the metallic kagome antiferromagnet FeSn

The kagome lattice is a fertile platform to explore topological excitations with both Fermi-Dirac and Bose-Einstein statistics. While relativistic Dirac fermions and flat bands have been discovered in the electronic structure of kagome metals, the spin excitations have received less attention. Here, we report inelastic neutron scattering studies of the prototypical kagome magnetic metal FeSn. The spectra display well-defined spin waves extending to 120 meV. Above this energy, the spin waves become progressively broadened, reflecting interactions with the Stoner continuum. Using linear spin-wave theory, we determine an effective spin Hamiltonian that reproduces the measured dispersion. This analysis indicates that the Dirac magnon at the K point remarkably occurs on the brink of a region where well-defined spin waves become unobservable. Furthermore, our results emphasize the influential role of itinerant carriers on the topological spin excitations of metallic kagome magnets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Tempering stochastic density functional theory

Here, we introduce a tempering approach with stochastic density functional theory (sDFT), labeled t-sDFT, which reduces the statistical errors in the estimates of observable expectation values. This is achieved by rewriting the electronic density as a sum of a "warm" component complemented by "colder" correction(s). Since the warm component is larger in magnitude but faster to evaluate, we use many more stochastic orbitals for its evaluation than for the smaller-sized colder correction(s). This results in a significant reduction in the statistical fluctuations and systematic deviation compared to sDFT for the same computational effort. We demonstrate the method's performance on large hydrogen-passivated silicon nanocrystals, finding a reduction in the systematic deviation in the energy by more than an order of magnitude, while the systematic deviation in the forces is also quenched. Similarly, the statistical fluctuations are reduced by factors of ≈4-5 for the total energy and ≈1.5-2 for the forces on the atoms. Since the embedding in t-sDFT is fully stochastic, it is possible to combine t-sDFT with other variants of sDFT such as energy-window sDFT and embedded-fragmented sDFT.

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

Energy window stochastic density functional theory

We report that linear scaling density functional theory is important for understanding electronic structure properties of nanometer scale systems. Recently developed stochastic density functional theory can achieve linear or even sublinear scaling for various electronic properties without relying on the sparsity of the density matrix. The basic idea relies on projecting stochastic orbitals onto the occupied space by expanding the Fermi-Dirac operator and repeating this for N χ stochastic orbitals. Often, a large number of stochastic orbitals are required to reduce the statistical fluctuations (which scale as N$-1/2\atop{χ}$) below a tolerable threshold. In this work, we introduce a new stochastic density functional theory that can efficiently reduce the statistical fluctuations for certain observable and can also be integrated with an embedded fragmentation scheme. The approach is based on dividing the occupied space into energy windows and projecting the stochastic orbitals with a single expansion onto all windows simultaneously. This allows for a significant reduction of the noise as illustrated for bulk silicon with a large supercell. Finally, we also provide theoretical analysis to rationalize why the noise can be reduced only for a certain class of ground state properties, such as the forces and electron density.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗