Search NASA⌕ Search

SEARCH · Search NASA

Results for “electron root”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Kinetics of Magnetic Skyrmion Crystal Formation from the Conical Phase

The particle-like magnetic skyrmion or skyrmion lattice (SkX) formation has promoted strong application and fundamental science interests. Despite extensive research, the kinetic of the SkX development is much less understood because of the ultrafast spin rotation and high sensitivity to external perturbations. In this work, using in situ Lorentz transmission electron microscopy, we successfully measured the dynamics of SkX formation from the conical phase with precise control of both the temperature and the magnetic field. We discovered that the Avrami equation can accurately describe the transition process with an initial Avrami constant around 1, suggesting that the rate-limiting step for the quasiparticle lattice formation is one-dimensional heterogeneous nucleation of individual skyrmions. A modified Arrhenius rate law is established, with an energy barrier that has a square-root dependence on temperature and a quadratic dependence on the magnetic field. This study paves the way toward precise and predictable manipulation of topological spin structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Block Krylov Subspace Methods for Functions of Matrices II: Modified Block FOM

We analyze an expansion of the generalized block Krylov subspace framework of [Electron. Trans. Numer. Anal., 47 (2017), pp. 100--126]. This expansion allows the use of low-rank modifications of the matrix projected onto the block Krylov subspace and contains, as special cases, the block GMRES method and the new block Radau--Arnoldi method. Within this general setting, we present results that extend the interpolation property from the nonblock case to a matrix polynomial interpolation property for the block case, and we relate the eigenvalues of the projected matrix to the latent roots of these matrix polynomials. Some error bounds for these modified block FOM methods for solving linear systems are presented. We then show how cospatial residuals can be preserved in the case of families of shifted linear block systems. This result is used to derive computationally practical restarted algorithms for block Krylov approximations that compute the action of a matrix function on a set of several vectors simultaneously. Finally, we prove some error bounds and present numerical results showing that two modifications of FOM, the block harmonic and the block Radau--Arnoldi methods for matrix functions, can significantly improve the convergence behavior.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The N‐terminal intrinsically disordered region of Ncb5or docks with the cytochrome b 5 core to form a helical motif that is of ancient origin

Abstract NADH cytochromeb 5 oxidoreductase (Ncb5or) is a cytosolic ferric reductase implicated in diabetes and neurological conditions. Ncb5or comprises cytochromeb 5 (b 5 ) and cytochromeb 5 reductase (b 5 R) domains separated by a CHORD‐Sgt1 (CS) linker domain. Ncb5or redox activity depends on proper inter‐domain interactions to mediate electron transfer from NADH or NADPH via FAD to heme. While full‐length human Ncb5or has proven resistant to crystallization, we have succeeded in obtaining high‐resolution atomic structures of theb 5 domain and a construct containing the CS andb 5 R domains (CS/b 5 R). Ncb5or also contains an N‐terminal intrinsically disordered region of 50 residues that has no homologs in other protein families in animals but features a distinctive, conserved L 34 MDWIRL 40 motif also present in reduced lateral root formation (RLF) protein in rice and increased recombination center 21 in baker's yeast, all attaching to ab 5 domain. After unsuccessful attempts at crystallizing a human Ncb5or construct comprising the N‐terminal region naturally fused to theb 5 domain, we were able to obtain a high‐resolution atomic structure of a recombinant rice RLF construct corresponding to residues 25–129 of human Ncb5or (52% sequence identity; 74% similarity). The structure reveals Trp 120 (corresponding to invariant Trp 37 in Ncb5or) to be part of an 11‐residue α‐helix (S 116 QMDWLKLTRT 126 ) packing against two of the four helices in theb 5 domain that surround heme (α2 and α5). The Trp 120 side chain forms a network of interactions with the side chains of four highly conserved residues corresponding to Tyr 85 and Tyr 88 (α2), Cys 124 (α5), and Leu 47 in Ncb5or. Circular dichroism measurements of human Ncb5or fragments further support a key role of Trp 37 in nucleating the formation of the N‐terminal helix, whose location in the N/b 5 module suggests a role in regulating the function of this multi‐domain redox enzyme. This study revealed for the first time an ancient origin of a helical motif in the N/b 5 module as reflected by its existence in a class of cytochromeb 5 proteins from three kingdoms among eukaryotes.

Biochemistry & Molecular Biology↗
DOE OSTI2026

Microstructure and Composition of Passivating Interfaces in Silicon Heterojunction Solar Modules Weathered in Different Climates

Sanyo/Panasonic patented silicon heterojunction with intrinsic thin layer (HIT) solar cells in the 1990's, which demonstrated world record photovoltaic (PV) efficiency around 2014 and inspired the newer generations of silicon heterojunction technology (SHJ) as well as the current world record back-contact PV cell designs. The high-quality passivation strategy, utilizing ultra-thin layers of hydrogenated amorphous silicon (a-Si:H), leads to high voltages and long carrier lifetimes, but these qualities may degrade over time as the modules operate outdoors. Here, we investigate local microstructure and composition at the interface layers of cells from HIT modules weathered in a hot, humid climate (Florida, USA) and temperate climate (Colorado, USA) for 10 years. We employ a comprehensive set of high-resolution electron microscopy imaging and spectroscopy to directly resolve structure and composition in these cells down to the nanoscale. We show features such as alignment of the In2O3 transparent conducting oxide and a-Si:H layers on the textured c-Si facets, interfacial oxidation at the a-Si:H/c-Si and In2O3/a-Si:H interfaces, and twinned c-Si resulting from epitaxial growth into the a-Si:H layer. Our results raise potential degradation mechanisms in these outdoor-weathered modules, but the root cause of electrical loss remains uncertain due to the very small changes in aged samples. This study shows that many atomic scale features at the In2O3/a-Si:H/c-Si interfaces are surprisingly robust, and the electrical losses with aging may either be attributed to other pathways or be very sensitive to the subtle chemical and microstructural features observed here.

14 SOLAR ENERGY↗
DOE OSTI2021

Atom Equivalent Energies for the Rapid Estimation of the Heat of Formation of Explosive Molecules from Density Functional Tight Binding Theory

Atom equivalent energies have been derived from which the gas-phase heat of formation of explosive molecules can be estimated from fast, semiempirical density functional tight binding total energy calculations. The root-mean-square deviation and maximum deviation of the heats of formation from the experimental values for the set of 45 energetic molecules compiled by Byrd and Rice [ J. Phys. Chem. A, 2006, 110, 1005–1013] are 10.4 and 25.5 kcal/mol, respectively, using 4 atom equivalent energies and 7.4 and 15.0 kcal/mol, respectively, using 7 atom equivalent energies. These errors are around a factor of 2–3 larger than those obtained from density functional theory calculations but are smaller than those obtained from other semiempirical electronic structure methods. Heats of formation calculated with density functional tight binding theory using the 4 and 7 atom equivalent energies, the Byrd and Rice scheme, and the atom pair contribution method for a new set of 531 energetic molecules that contain only carbon, hydrogen, nitrogen, and oxygen are provided.

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

Densification and microstructure features of lithium hydride fabrication

The manufacturing of lithium hydride (LiH) utilizing uniaxial pressing, which offers fabrication with tailorable properties via microstructure control, can lead to the expansion in application of LiH while bypassing the challenges presented by historical casting manufacturing techniques. Through control of consolidation conditions such as pressure, temperature, dwell time and powder load, the presented work highlights the densification of LiH, with an emphasis on quantifying oxygen content, for applications requiring a specific density range necessary for optimized material performance. Karl Fischer Titration and X-ray Diffraction proved useful in determining oxygen and phase content while Computed Tomography and Scanning Electron Microscopy provided structural analysis. The temperature dependent densification of LiH fit with an Arrhenius term resulted in an activation energy of 21.2 kJ/mol. Images of fractured surfaces of LiH pressed at 500 °C revealed drastic grain coarsening, aided by the presence of oxygen impurities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Using respiration quotients to track changing sources of soil respiration seasonally and with experimental warming

Developing a more mechanistic understanding of soil respiration is hampered by the difficulty in determining the contribution of different organic substrates to respiration and in disentangling autotrophic-versus-heterotrophic and aerobic-versus-anaerobic processes. Here, we use a relatively novel tool for better understanding soil respiration: the apparent respiration quotient (ARQ). The ARQ is the amount of CO 2 produced in the soil divided by the amount of O 2 consumed, and it changes according to which organic substrates are being consumed and whether oxygen is being used as an electron acceptor. We investigated how the ARQ of soil gas varied seasonally, by soil depth, and by in situ experimental warming ( +4 ° C) in a coniferous-forest whole-soil-profile warming experiment over 2 years. We then compared the patterns in ARQ to those of soil δ 13 CO 2 . Our measurements showed strong seasonal variations in ARQ, from ≈0.9 during the late spring and summer to ≈0.7 during the winter. This pattern likely reflected a shift from respiration being fueled by oxidized substrates like sugars and organic acids derived from root and root respiration during the growing season to more reduced substrates such as lipids and proteins derived from microbial necromass during the winter. This interpretation was supported by δ 13 CO 2 values, which were lower, like lipids, in the winter and higher, like sugars, in the summer. Furthermore, experimental warming significantly changed how both ARQ and δ 13 CO 2 responded to soil temperature. Wintertime ARQ and δ 13 CO 2 values were higher in heated than in control plots, probably due to the warming-driven increase in microbial activity that may have utilized oxidized carbon substrates, while growing-season values were lower in heated plots. Experimental warming and phenology change the sources of soil respiration throughout the soil profile. The sensitivity of ARQ to these changes demonstrates its potential as a tool for disentangling the biological sources contributing to soil respiration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

CIGS photovoltaics: reviewing an evolving paradigm

Copper indium selenide chalcopyrite-structure alloys with gallium (CIGS) are unique among the highest performing photovoltaic (PV) semiconductor technologies. They are structurally disordered, nonstoichiometric materials that have been engineered to achieve remarkably low bulk nonradiative recombination levels. Nevertheless, their performance can be further improved. This review adopts a fundamental thermodynamic perspective to comparatively assess the root causes of present limitations on CIGS PV performance. The topics of selectivity and passivation of contacts to CIGS and its multinary alloys are covered, highlighting pathways to maximizing the electrochemical potential between those contacts under illumination. Further, an overview of absorber growth methods and resulting properties is also provided. We recommend that CIGS researchers consider strategies that have been successfully implemented in the more mature wafer-based GaAs and Si PV device technologies, based on the paradigm of an idealized PV device design using an isotropic absorber with minimal nonradiative recombination, maximal light trapping, and both electron-selective and hole-selective passivated contacts. We foresee that CIGS technology will reach the 25% efficiency level within the next few years through enhanced collection and reduced recombination. To significantly impact power-generation applications, cost-effective, manufacturable solutions are also essential.

14 SOLAR ENERGY↗
DOE OSTI2021

Characterization of injection and confinement improvement through impurity induced profile modifications on the Wendelstein 7-X stellarator

Pulsed injections of boron carbide granules into Wendelstein 7-X stellarator (W7-X) plasmas transiently increase the plasma stored energy and core ion temperatures above the reference W7-X experimental programs by up to 30%. In a series of 4 MW electron cyclotron resonance heating experiments, the PPPL Probe Mounted Powder Injector provided 50 ms bursts of 100 μm granules every 350 ms at estimated quantities ranging from approximately 1 mg/pulse to over 30 mg/pulse. For each injection, the stored energy was observed to initially drop and the radiated power transiently increased, while the radial electron density profile rose at the edge as material was assimilated. Once the injected boron carbide was fully absorbed, the density rise transitioned to the core while the stored energy increased above the previous baseline level by an amount linearly correlated with the injection quantity. During the injection, the ion temperature gradient steepened with peak core ion temperatures observed to increase from a nominal 1.7 keV to over 2.6 keV for the largest injection amounts. Enhanced performance is accompanied by a reversal of the radial electric field at ρ < 0.3, indicating that the core transport has switched to the ion root. These observations are suggestive of a change in transport and provide further evidence that externally induced profile modifications provide a possible path to enhanced W7-X performance metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Explainable machine learning for hydrogen diffusion in metals and random binary alloys

Hydrogen diffusion in metals and alloys plays an important role in the discovery of new materials for fuel cell and energy storage technology. While analytic models use hand-selected features that have clear physical ties to hydrogen diffusion, they often lack accuracy when making quantitative predictions. Machine learning models are capable of making accurate predictions, but their inner workings are obscured, rendering it unclear which physical features are truly important. To develop interpretable machine learning models to predict the activation energies of hydrogen diffusion in metals and random binary alloys, we create a database for physical and chemical properties of the species and use it to fit six machine learning models. Our models achieve root-mean-squared errors between 98–119 meV on the testing data and accurately predict that elemental Ru has a large activation energy, while elemental Cr and Fe have small activation energies. By analyzing the feature importances of these fitted models, we identify relevant physical properties for predicting hydrogen diffusivity. While metrics for measuring the individual feature importances for machine learning models exist, correlations between the features lead to disagreement between models and limit the conclusions that can be drawn. Instead grouped feature importance, formed by combining the features via their correlations, agree across the six models and reveal that the two groups containing the packing factor and electronic specific heat are particularly significant for predicting hydrogen diffusion in metals and random binary alloys. In conclusion, this framework allows us to interpret machine learning models and enables rapid screening of new materials with the desired rates of hydrogen diffusion.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Limited thermal acclimation of photosynthesis in tropical montane tree species

Abstract The temperature sensitivity of physiological processes and growth of tropical trees remains a key uncertainty in predicting how tropical forests will adjust to future climates. In particular, our knowledge regarding warming responses of photosynthesis, and its underlying biochemical mechanisms, is very limited. We grew seedlings of two tropical montane rainforest tree species, the early‐successional species Harungana montana and the late‐successional species Syzygium guineense , at three different sites along an elevation gradient, differing by 6.8℃ in daytime ambient air temperature. Their physiological and growth performance was investigated at each site. The optimum temperature of net photosynthesis ( T optA ) did not significantly increase in warm‐grown trees in either species. Similarly, the thermal optima ( T optV and T optJ ) and activation energies ( E aV and E aJ ) of maximum Rubisco carboxylation capacity ( V cmax ) and maximum electron transport rate ( J max ) were largely unaffected by warming. However, V cmax , J max and foliar dark respiration ( R d ) at 25℃ were significantly reduced by warming in both species, and this decline was partly associated with concomitant reduction in total leaf nitrogen content. The ratio of J max / V cmax decreased with increasing leaf temperature for both species, but the ratio at 25℃ was constant across sites. Furthermore, in H . montana , stomatal conductance at 25℃ remained constant across the different temperature treatments, while in S . guineense it increased with warming. Total dry biomass increased with warming in H . montana but remained constant in S . guineense . The biomass allocated to roots, stem and leaves was not affected by warming in H . montana , whereas the biomass allocated to roots significantly increased in S . guineense . Overall, our findings show that in these two tropical montane rainforest tree species, the capacity to acclimate the thermal optimum of photosynthesis is limited while warming‐induced reductions in respiration and photosynthetic capacity rates are tightly coupled and linked to responses of leaf nitrogen.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Jordan–Wigner Transformation for the Description of Strong Correlation in Fermionic Systems

Seniority is a useful way of organizing Hilbert space for strongly correlated systems. The exact zero-seniority wave function, doubly occupied configuration interaction (DOCI), provides accurate results (given the right orbitals) for many strongly correlated electronic systems but has a combinatorial computational cost. In many cases, pair coupled cluster doubles provide a polynomial-cost approximation that closely reproduces the energies of DOCI, but it breaks down in some cases and, as shown herein, it does not provide particularly good density matrices. In this article, we demonstrate that by using the Jordan–Wigner transformation to turn the seniority zero problem back into a Fermionic one, we can provide mean-field variational results of DOCI quality for the Hubbard model and a few small molecular dissociation examples, with polynomial cost, both for the energies and for density matrices, all while being protected from collapse. This success is rooted in the proof we provide, showing that the Hartree–Fock wave function on the Jordan–Wigner-transformed Hamiltonian transforms back to variational coupled cluster doubles in the seniority zero representation, but restricted to have determinant rather than permanent amplitude coefficients, without compromising its overall accuracy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Beyond Volume Variation: Anisotropic and Protrusive Lithiation in Bismuth Nanowire

Materials storing energy via an alloying reaction are promising anode candidates in rechargeable lithium-ion batteries (LIBs) due to their much higher energy density than the current graphite anode. Until now, the volumetric expansion of such electrode particles during lithiation has been considered as solely responsible for cycling-induced structural failure. Here, we report different structural failure mechanisms using single-crystalline bismuth nanowires as the alloying-based anode. The Li–Bi alloying process exhibits a two-step transition, that is, Bi–Li 1 Bi and Li 1 Bi–Li 3 Bi. Interestingly, the Bi–Li 1 Bi phase transition occurs not only in the bulk Bi nanowire but also on the particle surface showing its characteristic behavior. The bulk alloying kinetics favors a Bi-(012)-facilitated anisotropic lithiation, whose mechanism and energetics are further studied using the density functional theory calculations. More importantly, the protrusion of Li 1 Bi nanograins as a result of anisotropic Li–Bi alloying is found to dominate the surface morphology of Bi particles. The growth kinetics of Li 1 Bi protrusions is understood atomically with the identification of two different controlling mechanisms, that is, the dislocation-assisted strain relaxation at the Bi/Li 1 Bi interface and the short-range migration of Bi supporting the off-Bi growth of Li 1 Bi. As loosely rooted to the bulk substrate and easily peeled off and detached into the electrolyte, these nanoscale protrusions developed during battery cycling are believed to be an important factor responsible for the capacity decay of such alloying-based anodes at the electrode level.

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

Antihydrogen chemistry

A survey of antimatter reactions is presented, including the formation of the antihydrogen atom and anionic, cationic, and molecular species by collisional and radiative processes. Our approach is rooted in the detailed knowledge available for many matter counterpart (hydrogenic) reactions, due to their importance in controlling early Universe chemistry. We point out that the availability of trapped antihydrogen at densities similar to those pertaining to the epoch of hydrogen chemistry will soon be available. In addition, using modern atomic physics techniques, it should be feasible to control antimatter in the laboratory to facilitate antihydrogen chemistry. Our purpose is to summarize what is known from hydrogen chemistry that is of relevance for antimatter and to indicate, based on possible reaction rates, which processes may be fruitful to pursue to create new antimatter entities as probes of fundamental symmetries. We include antihydrogen, positrons, and antiprotons in our discussion and additionally the electron due to its propensity to form positronium and perhaps to participate in certain reactions. We attempt to indicate whether further theoretical/computational work is necessary to add to the assessment of reaction rates, and we discount processes where the projected rates are too low to be of interest, given foreseeable experimental capabilities.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Evaluation and Demonstration of Intrusion Detection for Spent Fuel Storage Facilities

The IAEA recommends dual containment and surveillance (C/S) systems for difficult to access or difficult to measure spent fuel storage areas. However, many storage areas have limitations that prevent traditional secondary C/S systems, such as radiation levels too high for inspectors to apply individual Tamper Indicating Devices (TID) or have physical impedances to the use of TIDs. Recent developments in C/S technology include Laser Curtain for Containment (LCCT) technology implemented at Atucha-1 power plant in Argentina. The LCCT is a new approach to perimeter-wide C/S and the IAEA has granted provisional approval for its use as one of two required C/S systems. This proposal focuses on investigating and testing new and evolving completely passive (or ultra-low power) tamper indicating systems that will provide backup C/S to the LCCT even during power outages. An investigation of Commercial Off The Shelf (COTS) technology will be conducted to identify readily available options. Optically Stimulated Luminescence (OSL) fibers will also be investigated. OSL fibers are a passive technology that can measure radiation dose by knocking an electron into a metastable state where it remains until the fiber is interrogated with light. The goal of this project would be to investigate COTS technologies that can be used in low or no power scenarios for radiation detection as well as test and evaluate OSL technology, which can operate passively in no power scenarios. Limitations for each technology will be determined and documented in an effort to determine the best technology for spent fuel storage monitoring.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Accurate calculation of many-body energies in water clusters using a classical geometry-dependent induction model

We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The expansion is carried out up to hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA interaction potentials. The classical induction model with no adjustable parameters reproduces the ab-initio 3- and 4-body terms contained in these teo Databases with a Root-Mean-Square-Error (RMSE) of 0.104/0.058 and a Mean-Absolute-Error (MAE) of 0.054/0.026 kcal/mol, respectively, results that are on a par with those obtained 1 by fitting the same data using tens of thousands of Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy was included in the 3-body energy and was found to be small but not quite negligible. The model represents a practical, efficient and transferable approach for obtaining accurate non-additive interactions for multi-component systems without the need of performing tens of thousands of high level electronic structure calculations and fitting them with tens of thousands of PIPs.

molecular analysis water, aqueous↗
DOE OSTI2025

Skipper-CCD Quantum Efficiency Analysis

Scientific skipper-CCDs with single-electron resolution present dozens of possibilities for detecting dark matter candidates. Fermilab's Cosmic Physics Center contributes to the DarkNESS mission, which aims to place a skipper multi-chip module in a 6U CubeSat designed for low earth orbit. The DarkNESS nanosat will have the capability to search for 1-10keV band X-rays that may originate from DM decays. One of the challenges of detection in space is the large amount of cosmic radiation contributing to sensor noise. To mitigate this, an aluminum shield is proposed to be placed on the sensor. This project aims to test and characterize the energy resolution with a shield of various thicknesses (0-100nm) using a single CCD and an iron-55 x-ray source. ROOT analysis was used to parse data from several runs into sections based on shield thickness, create strategic data cuts, and characterize Fano plus signal shot noise in the sensor.

Wells, Megan E. [U. Illinois, Chicago]↗
DOE OSTI2026

Equivalent equilibrium ion distributions in ion sources and beams

Ion sources and ion beams often contain ions of more than one charge-to-mass ratio, due to the presence of multiple ionization states and/or multiple atomic and molecular ion species in the source. We show that, under certain conditions, analysis of these systems can be simplified using the notion of “equivalent equilibrium ion distributions” (EEIDs). These distributions have identical ion trajectories, space-charge distributions, and electric-field distributions, but can be composed of ion species with arbitrary charge-to-mass ratios. We use this result to obtain multi-species generalizations for the perveance of an ion diode, which relates the diode voltage to the extracted beam current, and for the perveance of a drifting ion beam, which determines the axial length-scale for beam expansion due to the self-electric field. Applied to a quasi-neutral plasma ion source containing cold ions and hot electrons, we find that the location of the beam-plasma interface (meniscus), where ions are extracted from the plasma, depends on the electron density and ion energy per unit charge at the interface, independent of the ion composition. For certain types of numerical simulations, EEIDs allow one to calculate the beam optics of an intense multi-species ion source by doing a simulation for a single ion species with an artificial charge-to-mass ratio chosen to minimize computation time. We illustrate these results using the 3D particle-in-cell, finite-element code Aleph. The simulations show how ions in a plasma source self-consistently accelerate to the supersonic speed needed for steady, radially expanding flow. Multi-species equilibria are subject to streaming instabilities because, while the different species have the same kinetic energy per unit charge, their velocities scale with the square root of the charge-to-mass ratio. The simulations show that a strong ion streaming instability can occur in a multi-species plasma ion source.

74 ATOMIC AND MOLECULAR PHYSICS↗