Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coulomb scattering”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

𝐴𝑏 initio density-matrix approach to exciton coherence: Phonon scattering, Coulomb interactions, and radiative recombination

Relaxation processes following light excitation in semiconductors are key in materials-based quantum technology applications. These processes are broadly studied in atomically thin transition-metal dichalcogenides, quasi-two-dimensional excitonic semiconductors in which atomistic design allows for tunable excited-state properties, such as relaxation lifetimes and photoinduced coherence. In this work, we present a density-matrix-based approach to compute exciton relaxation within a many-body ab initio perspective. We expand our previously developed Lindblad density-matrix formalism to capture multichannel electron-hole pair relaxation processes, including phonon and Coulomb scattering as well as radiative recombination, and we study their effect on the time-resolved excited-state propagation. Using monolayer MoSe 2 as a prototypical example, we examine many-body effects on the time-dependent dynamics of photoactive excitations, exploring how the electron-hole pair interactions are reflected in variations of the excitation energy, spectral signature, and state coherence. In conclusion, our method supplies a detailed understanding of exciton relaxation mechanisms in realistic materials, offering a previously unexplored pathway to study excited-state dynamics in semiconductors from first principles.

Band structure methods↗
DOE OSTI2025

Multiple Coulomb scattering in acrylic of a 221.3 MeV therapeutic proton beam

Measurements of multiple Coulomb scattering (MCS) distributions for 221.3 MeV therapeutic protons are presented using a novel detector system comprised of a thin scintillator, a pellicle mirror, and a digital camera. The MCS distributions were characterized for three acrylic phantoms of varying lengths and for two biological density-equivalent phantoms simulating bone and muscle. Additionally, beam profiles were measured across an energy range of 80.3–221.3 MeV in 20 MeV increments. The observed energy dependence of the photon yields is consistent with the tabulated stopping power values. Finally, the experimental results are benchmarked against Geant4 simulations, demonstrating consistent agreement and validating the capability of the detector system for radiology measurements.

Digital camera↗
DOE OSTI2022

A Cartesian-diffusion Langevin method for hybrid kinetic-fluid Coulomb scattering in particle-in-cell plasma simulations

A novel, drag-diffusion Langevin method of hybrid, kinetic-fluid Coulomb scattering in plasmas is presented. Unlike previous methods, the frictional drag is always applied in the simulation frame of reference. The velocity-space diffusion is performed in the stationary-fluid frame of reference when anisotropic, and in the laboratory frame when isotropic. While the general method is mass-ratio independent, we focus on interactions of kinetic-ions and fluid-electrons to show first-order modifications to the electron velocity distribution function that are an important correction for the accurate calculation of electric resistivity. Inclusion of sub-cycling and a limit to the maximum collision frequency is shown to allow for arbitrarily large timesteps without numerical failure. Here the Langevin method is compared with a grid-based binary method and found to require a much less restrictive timestep in cases of ion–electron slowing and temperature equilibration; this finding differs from previous work and is dependent on the mass ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Improved Muon Energy Estimation Using a Detailed Model of Multiple Coulomb Scattering in the MicroBooNE LArTPC

We present an improved technique for estimating a muon's energy by measuring the deflections along its path inside the MicroBooNE detector from multiple Coulomb scattering (MCS). This approach implements several innovations that better capture detector non-idealizations compared to previous MCS-based muon energy estimators. As a result, it achieves improved resolution, reduced bias, and better data-model agreement. Using model simulation, for fully contained events the estimated bias is within 1\% and the estimated resolution narrows from 10\% to 4.3\% as muon energy increases from 0.1\,GeV to 2\,GeV. For events with particles exiting the detector volume, at least a meter of reconstructed muon track, and a muon energy below 2\,GeV, the estimated bias is less than 2\% and the estimated resolution varies from 7\% to 17\% over muon energy. These demonstrate significant improvements over the performance of previous work using an MCS-based energy estimator at MicroBooNE~\cite{mcs_2017}, which exhibited approximately twice worse resolution and a bias of 20\% over the same energy region. Data-model goodness-of-fit studies are used to validate the estimator's performance on data, showing good agreement within model uncertainties.

Cooper-Troendle, London [U. Pittsburgh (main); Fer↗
DOE OSTI2026

Dynamics of ballistic photocurrents driven by Coulomb scattering in a two-dimensional material

First-principles real-time time-dependent density-functional theory (rt-TDDFT) calculations reveal the existence of ballistic photocurrents generated by Coulomb scattering, which has not previously been considered as a mechanism for the bulk photovoltaic effect. With monolayer GeS as an example, it is predicted that ballistic currents can be comparable to shift currents under experimentally accessible conditions.

2-dimensional systems↗
DOE OSTI2026

Multiple Coulomb scattering for muon energy reconstruction in SBND Monte Carlo simulations

The Short Baseline Near Detector (SBND) is a liquid argon time projection chamber (LArTPC) located 110 meters downstream of the Booster Neutrino Beam target at Fermilab \cite{SBN}. The LArTPC is a detector technology with unprecedented position and calorimetric resolution for precise neutrino studies. The most dominant neutrino interaction topology in SBND is charged-current quasi-elastic muon neutrino-argon interactions, where the final state consists of a muon and a proton. Here, muons often travel long distances, sometimes exiting the detector before stopping. This loss of information creates a challenge for accurate energy estimation. There have been efforts in the MicroBooNE and ICARUS collaborations to improve energy estimation using a phenomenon known as Multiple Coulomb Scattering (MCS), which describes the jagged path taken by muons through a detector due to continuous Coulomb interactions with argon atoms. Building on these efforts, my project aims to improve an MCS-based approach to estimating the energies of exiting muons in SBND.

Puig, Xilonem Angeles [Mt. Holyoke Coll.; Fermilab↗
DOE OSTI2020

A corrected method for Coulomb scattering in arbitrarily weighted particle-in-cell plasma simulations

A physically and statistically accurate method to perform binary, pairwise Coulomb collisions for arbitrarily-weighted macroparticles in particle-in-cell simulations is described. The method is shown to conserve global energy and momentum on average and to produce physically accurate scattering frequencies. This is justified using the statistically expected values of the method and through a series of numerical tests showing that the method is valid independently of macroparticle weight. This novel method is required due to an error in previous work on binary scattering, initially by Nanbu and Yonemura [J. Comp. Phys. 145, 639 (1998)] and applied to subsequent work, that causes inaccuracies in the scattering frequencies for unequally weighted simulations; a problem that occurs both for scattering within the group and between different groups.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Distribution function of nuclei from $\mathcal{e}^±$ scattering in the presence of a strong primordial magnetic field

The amplitude of the primordial magnetic field (PMF) is constrained from observational limits on primordial nuclear abundances. Within this constraint, it is possible that nuclear motion is regulated by Coulomb scattering with electrons and positrons (e ± s), while e ± s are affected by a PMF rather than collisions. For example, at a temperature of 10 9 K, thermal nuclei typically experience ~ 10 21 scatterings per second that are dominated by very small angle scattering leading to minuscule changes in the nuclear kinetic energy of order $\mathcal{O}$(1) eV. In this paper the upper limit on the effects of a possible discretization of the e ± momenta by the PMF on the nuclear momentum distribution is estimated under the extreme assumptions that the momentum of the e ± is relaxed before and after Coulomb scattering to Landau levels, and that during Coulomb scattering the PMF is neglected. This assumption explicitly breaks the time reversal invariance of Coulomb scattering, and the Maxwell-Boltzmann distribution is not a trivial steady state solution of the Boltzmann equation under these assumptions. Here we numerically evaluate the collision terms in the Boltzmann equation, and show that the introduction of a special direction in the e ± distribution by the PMF generates no directional dependence of the collisional destruction term of nuclei. Large anisotropies in the nuclear distribution function are then constrained from big bang nucleosynthesis. Ultimately, we conclude that a PMF does not significantly affect the isotropy or big bang nucleosynthesis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Estimating Residual Gas Beam Loss in the BNL AGS Booster:Model, Numerical Tools, and Implications for a Polarimeter Chamber

Beamloss from residual gas scattering is calculated for a storage ring of circumference C = 200m, covering a linear magnetic rigidity ramp from B ρ inj = 1.25Tm to B ρ ext = 9.50Tm over tramp = 1.00s. Three loss mechanisms are considered: nuclear inelastic reactions, single Coulomb scattering beyond the acceptance angle of the machine, and charge-state exchange for heavy ions. Calculations are performed for protons and Au 32+ ions, the latter being the charge state delivered by the EBIS pre-injector into the BNL AGS Booster. At a baseline vacuum of P 0 = 5 × 10 −11 mbar the integrated fractional beam loss over the full ramp is negligible for protons (< 0.001%) and 33.1% for Au 32+ . The impact of a 1m polarimeter section at elevated pressure P s = 5 × 10 −9 mbar is evaluated: the proton loss remains negligible, while the Au 32+ loss rises to 45.1%, establishing the vacuum requirement for the polarimeter chamber. In the present compact approximate model, Au 32+ losses are dominated by charge exchange at injection, while proton losses are dominated by single Coulomb scattering. The design, outgassing budget, pumping strategy, and commissioning sequence for a UHV polarimeter chamber housing silicon detectors and ASICs are described in detail.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Hyperdiffusion of dust particles in a turbulent tokamak plasma

In this work, the effect of plasma turbulence on the trajectories of dust particles is investigated for the first time. The dynamics of dust particles is computed using the ad hoc developed Dust Injection Simulator code, using a 3D turbulent plasma background computed with the TOKAM3X code. As a result, the evolution of the particle trajectories is governed by the ion drag force, and the shape of the trajectory is set by the Stokes number St∝a d /n 0 , with a d the dust radius and n 0 the density at the separatrix. The plasma turbulence is observed to scatter the dust particles, exhibiting a hyperdiffusive regime in all cases. The amplitude of the turbulent spread of the trajectories Δr 2 is shown to depend on the ratio Ku/St, with Ku∝u rms the Kubo number and urms the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting Δr 2 ∝(Ku/St) 2 t 3 , or Δr 2 ∝(u rms n 0 /a d ) 2 t 3 . As the dust is heated by the plasma fluxes, thermionic emission sets the dust charge, originally negative, to slightly positive values. This results in a substantial reduction of the ion drag force through the suppression of its Coulomb scattering component. The dust grain inertia is then no longer negligible and drives the transition from a hyperdiffusive regime toward a ballistic one.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Performance evaluation of cosmic ray muon trajectory estimation algorithms

Muons, being elementary particles with minimal interaction with nuclear materials and abundant at sea level, have sparked interest in utilizing them for imaging various applications, such as mining [Borselli et al., Sci. Rep. 12, 22329 (2022)], volcano imaging [Nagamine et al., Nucl. Instrum. Meth. A, 356, 585(1995)], and underground tunnel detection [Guardincerri et al., Pure Appl. Geophys. 174, 2133 (2017)]. Recently, their use in nuclear nonproliferation and safeguard verification has gained attention, particularly in cargo screening for nuclear waste smuggling [Baesso et al., J. Instrum. 9, C10041 (2014)], source localization [L. J. Schultz et al., Nucl. Instrum. Meth. A 519, 687 (2004)], and locating nuclear fuel debris in reactors [Borozdin et al., Phys. Rev. Let. 109, 152501 (2012)]. However, the resolution of muon image reconstruction techniques is limited due to multiple Coulomb scattering (MCS) within the target object. To achieve robust muon tomography, it is crucial to develop efficient and flexible physics-based algorithms that can model the MCS process accurately and estimate the most probable trajectory of muons as they pass through the target object. To address this limitation, in this study, a novel algorithmic approach utilizing the Bayesian probability theory and Gaussian approximation of MCS is chosen. Different energy levels, materials, and target sizes were considered in the evaluations. The results demonstrate that the Generalized Muon Trajectory Estimation (GMTE) algorithm offers significant improvements over currently used algorithms. Across all test scenarios, the GMTE algorithm demonstrated ~50% and 38% increase in precision compared to Straight Line Path (SLP) and Point of Closest Approach (PoCA) algorithms, respectively. Furthermore, it exhibited 10%–35% and 10%–15% increases in muon flux utilization for high and medium Z materials, respectively, compared to the PoCA algorithm. In conclusion, the extensive simulations confirm the enhanced performance and efficiency of the GMTE algorithm, offering improved resolution and reduced measurement time for cosmic ray muon imaging compared to the current SLP and PoCA algorithms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Absorption and Transport Effects Induced in Plasmas by the Interaction of Electrons with Laser Speckles

Here, we show that the ponderomotive force associated with laser speckles can scatter electrons in a laser-produced plasma in a manner similar to Coulomb scattering. Analytic expressions for the effective collision rates are given. The electron-speckle collisions become important at high laser intensity or during filamentation, affecting both long- and short-pulse laser intensity regimes. As an example, we find that the effective collision rate in the laser-overlap region of hohlraums on the National Ignition Facility is expected to exceed the Coulomb collision rate by 1 order of magnitude, leading to a fundamental change to the electron transport properties. At the high intensities characteristic of short-pulse laser-plasma interactions (I ≳ 10 17 W cm –2 ), the scattering is strong enough to cause the direct absorption of laser energy, generating hot electrons with energy scaling as E ≈ 1.44 (I/10 18 W cm –2 ) 1/2 MeV , close to experimentally observed results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

High temperature characteristics of nitric oxide annealed p-channel 4H-SiC metal oxide semiconductor field effect transistors

In this work, p-channel 4H-SiC MOSFETs were fabricated and analyzed at high temperature. It is demonstrated that nitridation of the gate oxide enables enhancement mode operation in these devices. Nitrogen incorporation at the 4H-SiC/SiO2 interface by nitric oxide annealing reduces the interface trap density energetically located in the lower half of the 4H-SiC bandgap, resulting in viable high temperature p-channel devices. In the 27–300 °C temperature range, the threshold voltage decreases with increasing temperature, consistent with the reduction of occupied interface traps at higher temperatures. The hole channel mobility is weakly temperature dependent under strong inversion conditions. Hall measurements support that above the threshold voltage, the mobility is limited by surface roughness scattering. In weaker inversion, the channel conductivity is limited by interface hole trapping and Coulomb scattering. In addition, high temperature bias stress measurements confirm a temperature activated hole trapping under negative gate bias, which requires further investigations.

Das, Suman (ORCID:0000000317146287)↗
DOE OSTI2023

Ab initio calculation of carrier mobility in semiconductors including ionized-impurity scattering

The past decade has seen the emergence of ab initio computational methods for calculating phonon-limited carrier mobilities in semiconductors with predictive accuracy. More realistic calculations ought to take into account additional scattering mechanisms such as, for example, impurity and grain-boundary scattering. Here, in this paper, we investigate the effect of ionized-impurity scattering on the carrier mobility. We model the analytical impurity potential parameterized from first principles by a collection of randomly distributed Coulomb scattering centers, and we include this relaxation channel into the ab initio Boltzmann transport equation, as implemented in the EPW code. We demonstrate this methodology by considering silicon, silicon carbide, and gallium phosphide, for which detailed experimental data are available. Our calculations agree well with experiments over a broad range of temperatures and impurity concentrations. For each compound investigated here, we compare the relative importance of electron-phonon scattering and ionized-impurity scattering, and we critically assess the reliability of Matthiessen's rule. We also show that an accurate description of dielectric screening and carrier effective masses can improve quantitative agreement with experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Muon momentum reconstruction in ICARUS-T600 LArTPC via Multiple Scattering

The SBN program at Fermilab investigates the possible existence of eV-scale sterile neutrinos. This work focuses on reconstructing muon momentum in $\nu$-Ar interactions within the ICARUS-T600 LArTPC detector using Multiple Coulomb Scattering (MCS), the only viable method for non-contained muons. Two algorithms are evaluated: the Gran Sasso algorithm (2D, Collection view) and the MicroBooNE algorithm (3D angles), applied to 2391 simulated and 2255 real stopping muon tracks (0.4-1 GeV/c). Performance is assessed via bias and resolution, using range-based momentum as a reference. Two enhanced versions of the Gran Sasso algorithm are also developed: one combining all three wire plane views, and one using 3D angles. Results show that the MicroBooNE algorithm performs slightly better, with resolution down to 14\% in simulation and 16\% in real data. These improvements significantly boost $\nu_\mu$ CC statistics, enhancing sensitivity to $\nu_\mu$ disappearance and 3+1 oscillation signals.

Chiello, Giovanni [Pisa U.] (ORCID:000900064571071↗
DOE OSTI2020

Thermal transport in compensated semimetals: Effect of electron-electron scattering on Lorenz ratio

It is well known that the electronic thermal conductivity of clean compensated semimetals can be greatly enhanced over the electric conductivity by the availability of an ambipolar mechanism of conduction, whereby electrons and holes flow in the same direction experiencing negligible Coulomb scattering as well as negligible impurity scattering. This enhancement—resulting in a breakdown of the Wiedemann-Franz law with an anomalously large Lorenz ratio—has been recently observed in two-dimensional monolayer and bilayer graphene near the charge neutrality point. In contrast to this, three-dimensional compensated semimetals such as WP 2 and Sb are typically found to show a reduced Lorenz ratio. Here, we investigate the reasons for this difference, focusing on the low-temperature regime where the electron-electron scattering is expected to dominate over other scattering mechanisms. We show that the different regimes of Fermi statistics (nondegenerate electron-hole liquid in graphene versus degenerate electron-hole liquid in compensated semimetals) are not sufficient to explain the reduction of the Lorenz ratio in the latter. We propose that the solution of the puzzle lies in the large separation of electron and hole pockets in momentum space, which allows compensated semimetals to sustain sizable regions of electron-hole accumulation near the contacts. These accumulations suppress the ambipolar conduction mechanism and effectively split the system into two independent electron and hole conductors. We present a quantitative theory of the crossover from ambipolar to unipolar conduction as a function of the size of the electron-hole accumulation regions, and show that it naturally leads to a sample-size-dependent thermal conductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗