Search NASA⌕ Search

SEARCH · Search NASA

Results for “BOLTZMANN DISTRIBUTION”

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

Fate of multiparticle resonances: From Q-balls to 3 He droplets

We study a system of N nonrelativistic particles which form a near-threshold resonance. Assuming no subset of these particles can form a bound state, the resonance can only decay through an “explosion” into N particles. We find that the decay width of the resonance scales as E Δ –5/2 in the limit when the energy E of the resonance goes to zero, where Δ is the ground-state energy of a system of N particles in a spherical harmonic trap with unit frequency. Here, the formula remains valid when some pairs of final particles have zero-energy s-wave resonance, but the Efimov effect is not present. In the limit of large N, we show that the final particles follow a Maxwell-Boltzmann distribution if they are bosons and a semicirclelike law if they are fermions. We expect our general result to be applicable to various systems that exist in nature. In particular, we argue that metastable 3 He droplets exist with the lifetime varying over many orders of magnitude ranging from a fraction of a nanosecond to values greatly exceeding the age of the Universe.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

High-temperature ion-thermal behavior from average-atom calculations

Atom-in-jellium calculations of the Einstein frequency were used to calculate the mean displacement of an ion over a wide range of compression and temperature. Expressed as a fraction of the Wigner-Seitz radius, the displacement is a measure of the asymptotic freedom of the ion at high temperature, and thus of the change in heat capacity from six to three quadratic degrees of freedom per atom. A functional form for free energy was proposed based on the Maxwell-Boltzmann distribution as a correction to the Debye free energy, with a single free parameter representing the effective density of potential modes to be saturated. Finally, this parameter was investigated using molecular dynamics simulations, and found to be ~0.2 per atom. In this way, the ion-thermal contribution can be calculated for a wide-range equation of state (EOS) without requiring a large number of molecular dynamics simulations. Example calculations were performed for carbon, including the sensitivity of key EOS loci to ionic freedom.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Urca Cooling in Neutron Star Crusts and Oceans: Effects of Nuclear Excitations

We report the excited-state structure of atomic nuclei can modify nuclear processes in stellar environments. In this Letter, we study the influence of nuclear excitations on Urca cooling (repeated back-and-forth β decay and electron capture in a pair of nuclear isotopes) in the crust and ocean of neutron stars. We provide for the first time an expression for Urca process neutrino luminosity which accounts for a thermal Boltzmann distribution of excited states in both members of an Urca pair. We use our new formula with state-of-the-art nuclear structure inputs to compute neutrino luminosities of candidate Urca cooling pairs. Our nuclear inputs consist of the latest experimental data supplemented with calculations using the projected shell model. We show that, in contrast to previous results that only consider the ground states of both nuclei in the pair, our calculated neutrino luminosities for different Urca pairs vary sensitively with the environment temperature and can be radically different from those obtained in the one-transition approximation. We find that nuclear excitations can lead to an enhancement in total Urca neutrino luminosities in the accreted neutron star crust by about 5 times as compared with the previous Urca results, which is expected to cause significant observational effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Simultaneous single-shot rotation–vibration non-equilibrium thermometry using pure rotational fs/ps CARS coherence beating

Here, we report the development of a simple and sensitive two-beam hybrid femtosecond/picosecond pure rotational coherent anti-Stokes Raman scattering (fs/ps CARS) method to simultaneously measure the rotational and vibrational temperatures of diatomic molecules. Rotation–vibration non-equilibrium plays a key role in the chemistry and thermalization in low-temperature plasmas as well as thermal loading of hypersonic vehicles. This approach uses time-domain interferences between ground state and vibrationally excited N 2 molecules to intentionally induce coherence beating that leads to apparent non-Boltzmann distributions in the pure rotational spectra. These distortions enable simultaneous inference of both the rotational and vibrational temperatures. Coherence beating effects were observed in single-shot fs/ps CARS measurements of a 75 Torr N 2 DC glow discharge and were successfully modeled for rotational and vibrational temperature extraction. We show that this method can be more sensitive than a pure rotational fs/ps CARS approach using a spectrally narrow probe pulse. Lastly, we experimentally measured the beat frequencies via Fourier transform of the time-domain response and obtained excellent agreement with the model.

42 ENGINEERING↗
DOE OSTI2024

One-dimensional temperature measurement of supersonic jet flow by resonantly ionized photoemission thermometry of molecular nitrogen

As the field of fluid dynamics progresses, the demand for sophisticated diagnostic methods to accurately assess flow conditions rises. In this work, resonantly ionized photoemission thermometry (RIPT) has been used to directly target and ionize diatomic nitrogen (N 2 ) to measure one-dimensional (1D) temperature profiles in a supersonic jet flow. This technique can be considered non-intrusive as the premise uses resonantly enhanced multiphoton ionization (REMPI) to target molecular nitrogen. This resonance excites N 2 into absorption bands of the P, Q, and R rotational branches of N 2 ( b 1 Π u ). The ideal (3 + 1) REMPI scheme excites from the ground state and ionizes N 2 ( b 1 Π u ← X 1 Σ g + ) where de-excitation results in photoemission from the first negative band of ionized N 2 + ( B 2 Σ u + → X 2 Σ g + ) as nitrogen returns to the ground state. The resulting emission can be observed using an intensified camera, thus permitting inference of the rotational temperature of ground-state molecular nitrogen. A linearly regressive Boltzmann distribution is applied based on previous calibration data for this technique to quantify the temperature along the ionized line. This work applies this technique to a pure N 2 supersonic jet in cross-flow and counter-flow orientations to demonstrate N 2 RIPT’s applications in a supersonic flow. Temperature variations are observed at different locations downstream of the exit in cross-flow, and axisymmetric in counter-flow, to generate profiles characterizing the flow dynamics. Due to the collisional effects resulting from the number density of N 2 at higher pressures, a (3 + 2) REMPI scheme is observed throughout this text.

Clark, Aleksander (ORCID:0009000855600042)↗
DOE OSTI2020

CI-MBPT and Intensity-Based Lifetime Calculations for Th II

Lifetime calculations of Th II J = 1.5 and 2.5 odd states are performed with configuration–interaction many-body perturbation theory (CI-MBPT). For many J = 2.5 states, lifetimes are quite accurate, but two pairs of J = 2.5 odd states and many groups of J = 1.5 states are strongly mixed, making theoretical predictions unreliable. To solve this problem, a method based on intensities is used. To relate experimental intensities to lifetimes, two parameters, one an overall coefficient of proportionality for transition rates and one temperature of the Boltzmann distribution of populations, are introduced and fitted to minimize the deviation between theoretical and intensity-derived lifetimes. For strongly mixed groups of states, the averaged lifetimes obtained from averaged transition rates were used instead of individual lifetimes in the fit. Close agreement is obtained. Then intensity branching ratios are used to extract individual lifetimes for the strongly mixed states. The resulting lifetimes are compared to available directly measured lifetimes and reasonable agreement is found, considering limited accuracy of intensity measurements. The method of intensity-based lifetime calculations with fit to theoretical lifetimes is quite general and can be applied to many complex atoms where strong mixing between multiple states exists.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Thermal scattering law for structure-dependent-Doppler broadening in FLASSH

Traditional Doppler broadening methods assume that target nuclei velocities follow a Maxwell-Boltzmann distribution which inherently assumes that the atoms are a free gas or that their velocities are independent of each other. This approximation is reasonable at high target temperatures and high neutron energies. However, nuclei are not independent: the lattice temperature and chemical binding will define the atomic motion. The thermal scattering law (i.e., TSL or S(α,β)) is a material property which describes the chemical binding and temperature response in terms of available momentum and energy states. In the thermal energy range, this TSL will define the thermal scattering cross sections. In the resonance region, the TSL offers a first-principles description of the probability distribution function for the velocity of the target nuclei. Using the TSL in Doppler broadening captures the structure of the material consistently from the thermal range into the resonance cross sections. In this work, both traditional free-gas and precise TSL Doppler broadening methods for resonance cross section evaluation have been implemented within the Full Law Analysis Scattering System Hub (FLASSH). This framework provides a generalized formulation for evaluating low-lying resonance data and streamlines nuclear data processing methods. TSL data generated using ab initio lattice dynamics (AILD) can be used to broaden cross section libraries in ENDF format for high-fidelity input into reactor physics calculations. These capabilities are demonstrated for the lowest absorption resonances of {sup 238}U in UO{sub 2}, UC and UN. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

An extended methodology for automated calculations of non-Boltzmann kinetic sequences: H + C 2 H 2 + X and combustion impact

It is generally assumed in phenomenological kinetic models that bimolecular reactions only occur between species whose rovibrational energy follows a Boltzmann (thermal) distribution. That is, any complexes initially formed in non-Boltzmann distributions are assumed to be thermalized by energy-transferring collisions prior to bimolecular reactions. Given the high mole fractions of reactive species, X, in combustion environments, reactive collisions of the complexes with X often occur on the same timescale as energy-transferring collisions – yielding sequences proceeding through non-Boltzmann intermediates across multiple potential energy surfaces. Recent studies have shown that such non-Boltzmann kinetic sequences can have substantial impact on the global reactivity in combustion systems. Simulations of these non-Boltzmann reaction sequences, which can be described in phenomenological kinetic models via chemically termolecular reactions, require that rovibrational excitation from one potential energy surface be carried over to the next. Here this paper presents an extended theoretical and computational methodology that couples multiple master equations and derives rate constants for phenomenological reactions describing the conversion of thermal reactants to thermal products for use in phenomenological kinetic schemes. The methodology is then implemented using in-house scripts for non-Boltzmann sequences involving C2H$^*_3$ + X (with X = O 2 , H, and OH) where C 2 H$^*_3$ is formed via H + C 2 H 2 association – which were identified as having strong potential for influencing combustion predictions in a recent study. The results reveal that non-Boltzmann reaction sequences for X = O 2 (the primary focus of this paper) significantly alters the total conversion rate from H + C 2 H 2 to products and product branching fractions from those of thermal sequential pathways. Furthermore, the present results demonstrate that non-Boltzmann reaction sequences have significant impact – as high as an order of magnitude – on predicted ignition delay times. Similarly, they yield significantly different dependence of ignition delay times with temperature and O 2 mole fraction – yielding signatures that are likely observable experimentally.

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

Adaptive hyperparameter updating for training restricted Boltzmann machines on quantum annealers

Restricted Boltzmann Machines (RBMs) have been proposed for developing neural networks for a variety of unsupervised machine learning applications such as image recognition, drug discovery, and materials design. The Boltzmann probability distribution is used as a model to identify network parameters by optimizing the likelihood of predicting an output given hidden states trained on available data. Training such networks often requires sampling over a large probability space that must be approximated during gradient based optimization. Quantum annealing has been proposed as a means to search this space more efficiently which has been experimentally investigated on D-Wave hardware. D-Wave implementation requires selection of an effective inverse temperature or hyperparameter (β) within the Boltzmann distribution which can strongly influence optimization. Here, we show how this parameter can be estimated as a hyperparameter applied to D-Wave hardware during neural network training by maximizing the likelihood or minimizing the Shannon entropy. We find both methods improve training RBMs based upon D-Wave hardware experimental validation on an image recognition problem. Neural network image reconstruction errors are evaluated using Bayesian uncertainty analysis which illustrate more than an order magnitude lower image reconstruction error using the maximum likelihood over manually optimizing the hyperparameter. The maximum likelihood method is also shown to out-perform minimizing the Shannon entropy for image reconstruction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Particle charge distributions in the effluent of a flow-through atmospheric pressure low temperature plasma

Atmospheric pressure low-temperature plasmas are often utilized to perform particle synthesis, treatment, and removal. It is well-known that dust particles are highly negatively charged in these plasmas; however, little is known about dust particle charging behavior as particles leave the plasma volume and pass through the spatial afterglow region. In this work, monodisperse particles of various sizes and work functions were introduced into an atmospheric pressure radiofrequency capacitively coupled flow-through plasma. Dust particle electrical mobility distributions downstream of the flow-through plasma were measured utilizing a differential mobility analyzer in conjunction with a condensation particle counter at various gas flow velocities. Charge distributions were determined from the measured electrical mobility distributions. Experiments confirm that particles become less negatively charged, and even net-positively charged after leaving the plasma volume, with a distribution that follows a shifted Boltzmann charge distribution. Additionally, particle charge in the effluent of the flow-through plasma is negligibly dependent on work function but highly size and flow velocity dependent. Larger particles were shown to have a higher magnitude of charge under all studied conditions; however, particle polarity was switchable by varying gas flow velocity. The charging dynamics were simulated utilizing a constant number Monte Carlo model that accounts for electron temperature decay and the transition from ambipolar to free diffusion of electrons and ions in the spatial afterglow. Simulation results also suggest that, at the same flow velocity, larger particles obtain a greater magnitude of charge, negative or positive. Furthermore, the decrease in electron mobility and the difference between ion and electron convective loss rates create an ion-rich region in the plasma effluent that promotes ion–particle collisions and drives particle charge removal and even reversal of polarity. Larger particles more favorably collide with energetic species in these environments, which results in higher charge states.

atmospheric pressure low temperature plasma↗
DOE OSTI2023

Double-sided electron energy analyzer for measurement of non-Maxwellian electron energy distributions

Here, a double-sided electron energy analyzer is developed for studies of magnetic reconnection. It can measure electron energy distribution functions along two directions opposite to each other at the same time. Each side is composed of a floating reference grid, an energy selector grid, and a collector plate. The voltage of the selector grid is swept from –40 to 0 V with respect to the reference grid with a frequency of 1 MHz. This fast sweeping is required to resolve sub-Alfvénic changes in plasma quantities of the Magnetic Reconnection Experiment, where the typical Alfvénic time is a few microseconds. The reliability of the energy analyzer is checked in Maxwellian plasmas away from the reconnection region. In this case, the electron temperature computed from the electron energy distribution function agrees with measurements of a reference triple Langmuir probe. When it is located near the reconnection region, the temperatures of the tail electron population from both sides, facing into the electron flow and facing away from it, exceed the bulk electron temperature measured by the Langmuir probe by a factor of about 2.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Analytic insights into nonlocal energy transport: Steady state Fokker Planck theory in arbitrary Z plasmas

The generation of energetic electrons in laser fusion in an important issue. The electrons may either arise from a laser plasma instability, or from the uncoupled high temperature tail of a Maxwellian distribution. To study these in a laser fusion context, it is important to find a method accurate enough to be useful, and simple enough to be incorporated into a radiation hydrodynamics numerical simulation, the main workhorse for studying the laser fusion target. That is why analytic insights become important, they allow one to simplify the Fokker Planck theory so that a solution of it can be incorporated into a radiation hydrodynamic simulation. This work develops and analyzes a steady state Fokker Planck theory for plasmas of arbitrary Z. It develops a method of solving the simplified Fokker Planck method with a technique called sparse eigenfunction analysis. As a result, this method appears to work reasonably well when compared to the experimental results from the Rochester/NIF on plastic spherical targets with and without a silicon layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Establishing criteria for the transition from kinetic to fluid modeling in hollow cathode analysis

Hollow cathodes for plasma switch applications are investigated via 2D3V particle-in-cell simulations of the channel and plume region. The kinetic nature of the plasma within the channel is dependent on the thermalization rate of electrons, emitted from the insert. When Coulomb collisions occur at a much greater rate than ionization or excitation collisions, the electron energy distribution function rapidly relaxes to a Maxwellian and the plasma within the channel can be described accurately via a fluid model. In contrast, if inelastic processes are much faster than Coulomb collisions, then the electron energy distribution function in the channel exhibits a notable high-energy tail, and a kinetic treatment is required. This criterion is applied to hollow cathodes from the literature, revealing that a fluid approach is suitable for most electric propulsion applications, whereas a kinetic treatment can be more critical to accurate modeling of plasma switches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

How to Build an Empirical Speed Distribution for Dark Matter in the Solar Neighborhood

The dark matter flux in a direct detection experiment depends on its local speed distribution. This distribution has been inferred from simulations of Milky Way–like galaxies, but such models serve only as proxies, given that no simulation directly captures the detailed evolution of our own Galaxy. This motivates alternative approaches that obtain this distribution directly from observations. In this work, we utilize 98 Milky Way analogues from the TNG50 simulation to develop and validate a procedure for inferring the dark matter speed distribution using the kinematics of nearby stars. We find that the dark matter that originated from old mergers, plus that from recent nonluminous accretions, is well described by a Maxwell–Boltzmann speed distribution centered at the local standard-of-rest velocity. Meanwhile, recently accreted dark matter from massive mergers has speeds that can be traced from the associated stellar debris of these events. The stellar populations systematically underestimate the velocity dispersion of their dark matter counterparts, but a simple kinematic boost brings the two into good alignment. Using the TNG50 host galaxies, we demonstrate that combining these two contributions provides an accurate reconstruction of the local dark matter speeds. As an application of the procedure to our own Galaxy, we utilize stellar kinematic data from Gaia to quantify how the dark matter remnants from the Milky Way’s last major merger impact its speed distribution in the solar neighborhood.

Shpigel, Tal [Princeton Univ., NJ (United States)]↗
DOE OSTI2022

Hybrid-gyrokinetic simulations of low- n toroidal Alfvén eigenmodes using gKPSP

Here, we report a benchmark study of toroidal Alfvén eigenmode (TAE) simulation using the hybrid-gyrokinetic code GyroKinetic Plasma Simulation Program (gKPSP). A simulation capability for energetic particles based on the gyrokinetic δf method has been newly implemented in the gKPSP code. Benchmark simulations have been performed in both circular and realistic tokamak geometries. Good agreement has been found with previously reported results, demonstrating the new capability of the gKPSP code. We have investigated the effects of the distribution function on TAE stability by examining both isotropic and anisotropic slowing-down distributions of energetic particles. The slowing-down distribution produces a higher linear growth rate than a Maxwellian distribution, while a growth rate scan with the anisotropy shows an opposite trend. This can be attributed to competition between Landau damping and the linear drive, which are correlated with the fraction of resonant passing particles and their distribution in phase space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Novel data interpretation method for DIII-D divertor retarding field energy analyzer with 3-D particle-in-cell simulations

A novel data interpretation process that utilizes comprehensive particle-in-cell (PIC) simulations is developed for the new retarding field energy analyzer (RFEA) currently being constructed at DIII-D for the lower divertor using the Divertor Material Evaluation System. Furthermore, this probe is expected to survive a heat load of up to 100 MW/m 2 for up to 5 s and reliably measure the main ion temperature (T i ) on the divertor target ranging from 10 to 200 eV. These extreme conditions posed significant engineering limitations on the probe geometry, thus extensive validation work has been performed. The conventional fitting method for the RFEA I–V characteristics is based on a simplified 1-D model without considering the ion space charge inside the probe cavity and may not be sufficient for probes designed for the DIII-D divertor environment. In this article, a more realistic description of the particle propagation process within the RFEA cavity is achieved by including both 3-D geometric effects and ion space charge in the PIC simulations, and the capability to reconstruct the ion energy distribution functions is demonstrated with reasonable consistency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗