Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coulomb”

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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Verification of a Monte Carlo binary collision model for simulating elastic and inelastic collisions in particle-in-cell simulations

We present the development and verification of a Monte Carlo binary collision model for simulating elastic and inelastic collisions in particle-in-cell simulations. We apply the corrected binary collision model originally developed for charged-particles collisions to all considered scattering channels, including Coulomb collisions, elastic neutral–neutral and charged–neutral collisions, ionization, excitation, and fusion. The model's implementation is described and verified through a series of simulations, including charged- and neutral-particle thermal equilibration, slowing of electrons in warm solid-density aluminum, collisional damping of a Langmuir wave, helium gas breakdown in an applied electric field, and thermonuclear and beam–target fusion. Then, we demonstrate the model within simulations of hydrogen plasma formation in the Princeton Field-Reversed Configuration as well as of the burning of aneutronic fusion fuel p-11B. The latter includes measurement of the fusion power density in a low-density plasma and fusion production due to the stopping of a proton ignitor beam in a compressed boron target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Ultrafast laser triggered electron emission from ultrananocrystalline diamond pyramid tip cathode

Nitrogen-incorporated ultrananocrystalline diamond [(N)UNCD] pyramid tip cathode has been considered as a next-generation high peak current electron source for dielectric laser accelerators as well as other high peak current particle accelerator applications. In this work, we study non-linear photoemission from an (N)UNCD pyramid tip cathode using an ultrafast laser with the pulse length of 150 fs with the central wavelength of 800 nm in the peak intensity range of 10 9 –10 10 W/cm 2 . We demonstrated that as the incident laser intensity increases, the current emitted from the nano-tip first increases as a power function with an exponent of about 5 and then starts to roll over to an exponent of 3. This roll over is attributed to the Coulomb interaction between electrons emitted from the tip also known as the space charge. We also measured the photoemission electron energy spectra that show electrons with energies as high as ~10 eV. Based on the shape of the electron energy spectra, we conclude that the high-energy electrons are thermally emitted electrons due to ultrafast laser heating at the tip of the (N)UNCD pyramid tip cathode.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Acceleration mechanisms of energetic ion debris in laser-driven tin plasma EUV sources

Laser-driven tin plasmas are driving new-generation nanolithography as sources of extreme ultraviolet (EUV) radiation centered at 13.5 nm. A major challenge facing industrial EUV source development is predicting energetic ion debris produced during the plasma expansion that may damage the sensitive EUV channeling multilayer optics. Gaining a detailed understanding of the plasma dynamics and ion acceleration mechanisms in these sources could provide critical insights for designing debris mitigation strategies in future high-power EUV sources. Here, we develop a fully kinetic model of tin-EUV sources using one-dimensional particle-in-cell simulations to study ion debris acceleration, which will be valuable for cross-validation of radiation-hydrodynamic simulations. An inverse-bremsstrahlung heating operator is used to model the interaction of a tin target with an Nd:YAG laser, and thermal conduction is included through a Monte Carlo Coulomb collision operator. While the large-scale evolution is in reasonable agreement with analogous hydrodynamic simulations, the significant timescale for collisional equilibration between electrons and ions allows for the development of prominent two-temperature features. A collimated flow of energetic ions is produced with a spectrum that is significantly enhanced at high energies compared to fluid simulations. The dominant acceleration mechanism is found to be a large-scale electric field supported mainly by the electron pressure gradient, which is enhanced in the kinetic simulations due to the increased electron temperature. We discuss the implications of these results for future modeling of tin-EUV sources and the development of debris mitigation schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Mean force emission theory for classical bremsstrahlung in strongly coupled plasmas

This work presents mean force emission theory, which extends the classical theory of bremsstrahlung emission to strongly coupled plasmas. In the high-frequency limit, the theory reduces to solving for the electron trajectory during a binary collision, but where the electron–ion interactions occur through the potential of mean force. In the low-frequency limit, it uses an autocorrelation formalism that captures effects of multiple collisions and strongly correlated motion. The predictions are benchmarked by comparison with first-principles classical molecular dynamics simulations of a fully ionized hydrogen plasma in which all interactions are repulsive. The comparison shows good agreement up to Coulomb coupling strengths of Γ∼30. The theory improves upon traditional models by including strong coupling effects and systematically including the effect of multiple collisions. Furthermore, mean force emission theory provides evidence that the Drude correction factor commonly used in quantum calculations of optical quantities may not be adequate at strong coupling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Verification of the kinetic electron role in the microinstabilities in a negative triangularity model equilibrium

Effect of kinetic electrons on negative triangularity plasmas has been investigated and compared against the corresponding positive triangularity plasmas, using the global gyrokinetic code X-point Gyrokinetic Code with scale-separated delta-f option without Coulomb collisions. Our model magnetic equilibria have strong positive and negative triangularities and weak magnetic shear. However, unusually large ρ i > a and low density plasmas are chosen to maximize the nonlocal effect to investigate the finite ρ i effect and to be clearly away from kinetic ballooning modes. Similar conclusions to previous flux tube and global simulations have been obtained in this highly nonlocal model plasma: it is essential to include kinetic electrons in the micro-instability study of negative triangularity plasmas. Most physics findings agree with existing reports, with some disagreement. We offer a new “effective trapping fraction” concept that can add to the explanation of the growth rate difference between NT and PT plasmas, pointing to the significant variation in trapped particle fractions that have turning points in the mode growth regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Alternative CNDOL Fockians for fast and accurate description of molecular exciton properties

CNDOL is an a priori, approximate Fockian for molecular wave functions. In this study, we employ several modes of singly excited configuration interaction (CIS) to model molecular excitation properties by using four combinations of the one electron operator terms. Those options are compared to the experimental and theoretical data for a carefully selected set of molecules. The resulting excitons are represented by CIS wave functions that encompass all valence electrons in the system for each excited state energy. The Coulomb–exchange term associated to the calculated excitation energies is rationalized to evaluate theoretical exciton binding energies. This property is shown to be useful for discriminating the charge donation ability of molecular and supermolecular systems. Multielectronic 3D maps of exciton formal charges are showcased, demonstrating the applicability of these approximate wave functions for modeling properties of large molecules and clusters at nanoscales. This modeling proves useful in designing molecular photovoltaic devices. Our methodology holds potential applications in systematic evaluations of such systems and the development of fundamental artificial intelligence databases for predicting related properties.

Chemistry↗
DOE OSTI2024

GW with hybrid functionals for large molecular systems

A low-cost approach for stochastically sampling static exchange during time-dependent Hartree–Fock-type propagation is presented. This enables the use of an excellent hybrid density functional theory (DFT) starting point for stochastic GW quasiparticle energy calculations. Generalized Kohn–Sham molecular orbitals and energies, rather than those of a local-DFT calculation, are used for building the Green function and effective Coulomb interaction. The use of an optimally tuned hybrid diminishes the starting point dependency in one-shot stochastic GW, effectively avoiding the need for self-consistent GW iterations.

Chemistry↗
DOE OSTI2024

A molecular ground electronic state with an occupied 5g spinor—The superheavy (E125)F molecule

Fully relativistic calculations, primarily at the 4-component coupled-cluster singles and doubles with perturbative triples [CCSD(T)] level of theory with the Dirac–Coulomb (DC) Hamiltonian, have been carried out for the superheavy (E125)F molecule using large Gaussian basis sets. The electronic ground state is determined to have an [Og]8s 2 5g 1 6f 3 configuration on E125 with an Ω = 6 ground state and an 8p electron largely donated to F. A Mulliken population analysis indicates that the ground state is mainly ionic with a partial charge of +0.79 on E125 and a single sigma bond involving the F 2p and E125 8p spinors. The occupied g spinor is not involved in the bonding. With the largest basis set used in this work, the (0 K) dissociation energy was calculated at the DC-CCSD(T) level of theory to be 7.02 eV. Analogous calculations were also carried out for the E125 atom, both the neutral and its cation. The lowest energy electron configuration of E125 + , [Og]$8s$$^{2}_{1/2}$$5g$$^{1}_{7/2}$$6f$$^{3}_{5/2}$ with a J = 6 ground state, was found to be similar to that in (E125)F, while the neutral E125 atom has an [Og]$8s$$^{2}_{1/2}$$5g$$^{1}_{7/2}$$6f$$^{2}_{5/2}$$7d$$^{1}_{3/2}$$8p$$^{1}_{1/2}$ ground state electron configuration with a J = 17/2 ground state. The ionization energy (IE) of E125 is reported for the first time and is calculated to be 4.70 eV at the DC-CCSD(T) level of theory. Non-relativistic calculations were also carried out on the E125 atom and the (E125)F molecule. Here, the non-relativistic ground state of the E125 atom was calculated to have a 5g 5 ground state with an IE of just 3.4 eV. The net effect of relativity on (E125)F is to stabilize its bonding.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

The XPS of pyridine: A combined theoretical and experimental analysis

A detailed analysis of the N(1s) and C(1s) X-Ray Photoelectron Spectroscopy (XPS) is made, where the measured XPS is compared with theoretical Sudden Approximation (SA) intensities and theoretical XPS Binding Energies (BEs). There is remarkably good agreement between the theoretical predictions and the measured XPS; in particular, the different full width at half maximum values for the C(1s) and N(1s) BEs are explained in terms of unresolved C(1s) BEs for the different C atoms in pyridine. Here, this work demonstrates that the combination of theory and XPS measurements can extract analysis of the XPS relevant to the molecular electronic structure. The theory used is based on fully relativistic self-consistent field solutions of the Dirac–Coulomb Hamiltonian, and the SA is used to determine relative XPS intensities.

Binding Energies (BEs)↗
DOE OSTI2025

Revealing the potential of nickel zinc ferrite: Facile synthesis and cost-effective anode for lithium-ion batteries

The increasing demand for lithium-ion batteries (LIBs) underscores the need for cost-effective alternative anode materials to ensure efficient Li-ion storage, given their pivotal role in various industries. This research focuses on the facile synthesis of nickel zinc ferrite (NZFO: Ni 0.65 Zn 0.35 Fe 2 O 4 ) and conducts comprehensive electrochemical analyses to evaluate its potential as a high-capacity alternative anode material for LIBs. The NZFO-CMR [sodium carboxymethyl cellulose (2%) and styrene-butadiene rubber (1%)] exhibited an initial delithiated capacity of ∼1232 mA h g −1 and maintained a stable capacity of around 358 mA h g −1 , along with an average Coulombic efficiency of 99.6% over 200 cycles. Cyclic voltammetry analysis revealed that Li-ion insertion was predominantly governed by ion diffusion, and the consistent correlation observed in electrochemical impedance spectroscopy spectra indicated stable electrochemical behavior throughout cycling. The facile synthesis approach and reasonable electrochemical performance of NZFO suggest its potential as an alternative anode material for advancing LIB's technology.

25 ENERGY STORAGE↗
DOE OSTI2025

Impact of quantum well thickness on efficiency loss in InGaN/GaN LEDs: Challenges for thin-well designs

We investigate the impact of quantum well (QW) thickness on efficiency loss in c-plane InGaN/GaN LEDs using a small-signal electroluminescence technique. Multiple mechanisms related to efficiency loss are independently examined, including injection efficiency, carrier density vs current density relationship, phase space filling, quantum-confined Stark effect, and Coulomb enhancement. An optimal QW thickness of around 2.7 nm in these InGaN/GaN LEDs was determined for QWs having constant In composition. Despite improved control of deep-level defects and lower carrier density at a given current density, LEDs with thin QWs still suffer from an imbalance in enhancement effects on the radiative and intrinsic Auger–Meitner recombination coefficients. The imbalance in enhancement effects results in a decline in internal quantum efficiency and radiative efficiency with decreasing QW thickness at low current density in LEDs with QW thicknesses below 2.7 nm. Here, we also investigate how LED modulation bandwidth varies with QW thickness, identifying the key trends and their implications for device performance.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Similarity for downscaled kinetic simulations of electrostatic plasmas: Reconciling the large system size with small Debye length

A simple similarity has been proposed for kinetic (e.g., particle-in-cell) simulations of plasma transport that can effectively address the long-standing challenge of reconciling the tiny Debye length with the vast system size. This applies to both transport in unmagnetized plasma and parallel transport in magnetized plasmas, where the characteristics length scales are given by the Debye length, collisional mean free paths, and the system or gradient lengths. The controlled scaled variables are the configuration space, x/L, and an artificial Coulomb Logarithm, L ln Λ, for collisions, while the scaled time, t/L, and electric field, LE, are automatic outcomes. The similarity properties are examined, demonstrating that the macroscopic transport physics is preserved through a similarity transformation while keeping the microscopic physics at its original scale of Debye length. To showcase the utility of this approach, two examples of 1D plasma transport problems were simulated using the VPIC code: the plasma thermal quench in tokamaks [Li et al., Nuclear Fusion 63, 066030 (2023)] and the plasma sheath in the high-recycling regime [Li et al., Physics of Plasmas 30, 063505 (2023)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Proton radiation effects in indium oxide using cascade molecular dynamics simulations

Metal oxide (MO) semiconductors, characterized by their wide band gaps and notable charge transport properties, are promising candidates for electronic applications in extreme environments, including near-Earth space. However, atomistic simulations of radiation–matter interactions in MOs remain challenging due to the limitations of existing interatomic potentials, which often fail to capture both the short-range repulsive forces essential for radiation damage modeling and the long-range electrostatic effects governing defect evolution. In this work, we develop a customized interatomic potential tailored for radiation damage simulations in indium oxide (In 2 O 3 ) as a model system, a representative MO material. Our potential integrates the Ziegler-Biersack-Littmark potential to accurately describe short-range interactions with Buckingham and Coulombic potentials to account for long-range forces. We perform molecular dynamics simulations of low-energy proton irradiation using this custom potential. We employ the primary knock-on atom (PKA) cascade method to study atomic displacements and primary defect formation. Simulations were conducted for 1 keV proton irradiation in four randomly chosen directions, and PKA-driven defect analyses at 5, 10, and 15 keV to examine the effects of direction and energy level on damage generation. Our results provide insight into the impact of irradiation direction and energy level on the cascade evolution and defect formation mechanisms.

Atomistic simulations↗
DOE OSTI2025

Imaging transient molecular configurations in UV-excited diiodomethane

Femtosecond structural dynamics of diiodomethane (CH 2 I 2 ) triggered by ultraviolet (UV) photoabsorption at 290 and 330 nm are studied using time-resolved coincident Coulomb explosion imaging driven by a near-infrared probe pulse. Here, we map the dominant single-photon process, the cleavage of the carbon–iodine bond producing rotationally excited CH 2 I radical, identify the contributions of the three-body (CH 2 + I + I) dissociation and molecular iodine formation channels, which are primarily driven by the absorption of more than one UV photon, and demonstrate the existence of a weak reaction pathway involving the formation of short-lived transient species resembling iso-CH 2 I 2 -like geometries with a slightly shorter I–I separation compared to the ground-state CH 2 I 2 . These transient molecular configurations, which can be separated from the other channels by applying a set of conditions on the correlated momenta of three ionic fragments, are formed within ∼100 fs after the initial photoexcitation and decay within the next 100 fs.

Atomic and molecular processes↗
DOE OSTI2025

Simulations of muon imaging with the LANL GMT detector for spent nuclear fuel cask content verification

Atmospheric muons are typically high energy, highly penetrating charged particles. They interact with matter primarily through multiple Coulomb scatterings. Muon scattering intensities can be used to characterize the density and atomic number of the matter that they pass through. Previously, the Los Alamos National Laboratory (LANL) muon tomography team performed muon imaging of the partially filled MC-10 spent nuclear fuel (SNF) cask at Idaho National Laboratory (INL). This experiment demonstrated the feasibility of muon imaging for the verification of spent fuel container contents. That original effort used the mini muon tracker array, consisting of two arrays of drift tubes on either side of the SNF cask. The reconstructed image quality was limited by statistics, largely due to low muon flux at high zenith angles. A LANL led team will perform new measurements with a larger array, the Giant Muon Tracker (GMT), to improve data collection rates and statistics. In this work, simulations were performed with the GMT near the partially filled INL MC-10 cask. For more general fuel diversion detection, a full MC-10 cask and casks with a singular missing fuel bundle were also simulated. To understand minimum measurement times needed for missing bundle identification, 100 000 to millions of tracked muons (corresponding to 1.4 days to several weeks measurement time) were analyzed. Simulated images were then analyzed visually and numerically to explore techniques designed to minimize the collection time needed to identify the diversion of fuel in each scenario.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Scattering theory in noncanonical phase space: A Drift-Kinetic collision operator for weakly collisional plasmas

After developing a scattering theory for grazing collisions in general noncanonical phase spaces, we introduce a guiding center collision operator in five-dimensional phase space designed for plasma regimes characterized by long wavelengths (relative to the Larmor radius), low frequencies (relative to the cyclotron frequency), and weak collisionality (where repeated Coulomb collisions induce cumulatively small changes in particle magnetic moment). The collision operator is fully determined by the noncanonical Hamiltonian structure of guiding center dynamics and exhibits a metriplectic structure, ensuring the conservation of particle number, momentum, energy, and interior Casimir invariants. It also satisfies an H-theorem, allowing for deviations from an equilibrium Maxwellian distribution due to the nontrivial kernel of the noncanonical guiding center Poisson tensor, spanned by the magnetic moment. We propose that this collision operator and its underlying mathematical structure may offer valuable insight into the study of turbulence, transport, and self-organizing phenomena in both laboratory and astrophysical plasmas.

Hamiltonian mechanics↗
DOE OSTI2025

Shock wave bending around a dusty plasma void

We report on experimental observations of the bending of a dust acoustic shock wave around a dust void region. This phenomenon occurs as a planar shock wavefront encounters a compressible obstacle in the form of a void whose size is larger than the wavelength of the wave. As they collide, the central portion of the wavefront, that is the first to touch the void, is blocked while the rest of the front continues to propagate, resulting in an inward bending of the shock wave. The bent shock wave eventually collapses, leading to the transient trapping of dust particles in the void. Subsequently, a Coulomb explosion of the trapped particles generates a bow shock. The experiments have been carried out in a Direct Current (DC) glow discharge plasma, where the shock wave and the void are simultaneously created as self-excited modes of a three-dimensional dust cloud. The salient features of this phenomenon are reproduced in molecular dynamics simulations, which provide valuable insights into the underlying dynamics of this interaction.

Dust acoustic waves↗
DOE OSTI2026

Impact of atomic substitution on core-hole relaxation dynamics: A study of Br 2 and IBr

Understanding inner-shell decay processes in heavy-element molecules is essential for unraveling x-ray-induced photodynamics and advancing molecular imaging techniques. Here, in this study, we investigate the influence of atomic substitution on core-hole relaxation dynamics and molecular fragmentation in Br 2 and IBr, initiated by x-ray ionization absorption at the Br K-edge. Using a combination of x-ray/ion coincidence measurements and Monte Carlo/molecular dynamics simulations, we track the charge distribution and the kinetic energy release (KER) of fragment ions with a total charge from 2+ to 8+. For both molecules, the simulated KER values show good agreement with experiment across different fragmentation channels. Our comparison reveals that substituting Br with the heavier I atom in IBr has a minimal impact on the inner-shell electronic decay process but significantly influences nuclear motion, leading to slower dissociation and thus a KER close to the Coulomb limit—an effect attributed to the atomic mass. These findings highlight the interplay between electronic and nuclear effects in molecular fragmentation, particularly in heavy-element species, and provide new insights into medical therapies, structural biology, and astrophysics.

Bhat, Nivedita [Argonne National Laboratory (ANL),↗