Search NASA⌕ Search

SEARCH · Search NASA

Results for “NEUTRAL ELEMENT”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

FLAG Review 2019: Flavour Lattice Averaging Group (FLAG)

Abstract We review lattice results related to pion, kaon, D -meson, B -meson, and nucleon physics with the aim of making them easily accessible to the nuclear and particle physics communities. More specifically, we report on the determination of the light-quark masses, the form factor $$f_+(0)$$ f + ( 0 ) arising in the semileptonic $$K \rightarrow \pi $$ K → π transition at zero momentum transfer, as well as the decay constant ratio $$f_K/f_\pi $$ f K / f π and its consequences for the CKM matrix elements $$V_{us}$$ V us and $$V_{ud}$$ V ud . Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of $$SU(2)_L\times SU(2)_R$$ S U ( 2 ) L × S U ( 2 ) R and $$SU(3)_L\times SU(3)_R$$ S U ( 3 ) L × S U ( 3 ) R Chiral Perturbation Theory. We review the determination of the $$B_K$$ B K parameter of neutral kaon mixing as well as the additional four B parameters that arise in theories of physics beyond the Standard Model. For the heavy-quark sector, we provide results for $$m_c$$ m c and $$m_b$$ m b as well as those for D - and B -meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. We review the status of lattice determinations of the strong coupling constant $$\alpha _s$$ α s . Finally, in this review we have added a new section reviewing results for nucleon matrix elements of the axial, scalar and tensor bilinears, both isovector and flavor diagonal.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Recovering Rare Earth Elements from Coal Mine Drainage Using Industrial Byproducts: Environmental and Economic Consequences

Coal mine drainage (CMD) impairs tens of thousands of kilometers of U.S. waterways each year, in part with the leaching of low concentrations of rare earth elements (REEs). REEs are essential for modern technologies, yet economically viable natural deposits are geospatially limited, thus engendering geopolitical concerns, and their mining is energy intense and environmentally destructive. This work summarizes laboratory-scale experimentalresults of a trap-extract-precipitate (TEP) process and uses the mass and energy balances to estimate the economic costs and environmental impacts of the TEP. The TEP process uses the alkalinity and filtering capacity of stabilized flue gas desulfurization (sFGD) material or water treatment plant (WTP) sludge to remediate CMD waters and extract REEs. Passive treatment systems that use WTP sludge are cheaper than those that use sFGD material ($\$$89,300/year or $\$$86/gT-REE vs. $\$$89,800/year or $\$$278/gT-REE) and have improved environmental performance across all indicators from two different impact assessment methods. These differences are largely attributable to the larger neutralizing capacity of WTP sludge in the treatment application.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

A reactive electrochemomechanical theory for growth and remodeling of polyelectrolyte hydrogels and application to dynamic polymerization of $\text{DNA}$ hydrogels

Here this study develops a framework for growth and remodeling of active polyelectrolyte hydrogels that accounts for effects of compositional changes on the mechanical response. By developing a reactive electrochemomechanical theory, thermodynamical constraints upon reactive and remodeling processes are elucidated within a general framework that allows any number of chemical reactions to evolve the response of the gel and transfer mass and charge between constituents. Fully coupled, nonlinear constitutive relations are adopted for molar fluxes, allowing exploration of effects including cross-diffusion, electrophoresis, and electro-osmosis. A robust finite element implementation is developed in the open source FEBio software (febio.org) by exploiting an equivalence between electrochemomechanics and mixture theory. The implementation is verified against analytical solutions for free swelling, and a proper reduction to a prior chemomechanical theory is demonstrated for neutral gels swollen only by a solvent with no solutes. The theory and implementation are then applied to model the tunable large swelling achieved through dynamic polymerization of DNA crosslinkers seen in our recently developed experimental hydrogel system (Cangialosi et al., 2017). A novel constitutive model for reaction-driven evolution of the locking stretch λ L in a non-Gaussian mechanical free energy was developed, where the increasing concentration of DNA crosslinkers makes further swelling energetically favorable. With a single free parameter, excellent agreement was found between measured and predicted equilibrium swelling ratios. This study demonstrated the ability to extend the electrochemomechanical framework to include chemical reactions and composition-aware constitutive models, and showed that development of reactive models allows simulation of complex dynamic polymerization phenomena not treated before. The theoretical frame here can be further expanded in scope to incorporate additional non-ideal and nonlinear phenomena.

42 ENGINEERING↗
DOE OSTI2022

FLAG Review 2021

Abstract We review lattice results related to pion, kaon, D -meson, B -meson, and nucleon physics with the aim of making them easily accessible to the nuclear and particle physics communities. More specifically, we report on the determination of the light-quark masses, the form factor $$f_+(0)$$ f + ( 0 ) arising in the semileptonic $$K \rightarrow \pi $$ K → π transition at zero momentum transfer, as well as the decay constant ratio $$f_K/f_\pi $$ f K / f π and its consequences for the CKM matrix elements $$V_{us}$$ V us and $$V_{ud}$$ V ud . Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of $$SU(2)_L\times SU(2)_R$$ S U ( 2 ) L × S U ( 2 ) R and $$SU(3)_L\times SU(3)_R$$ S U ( 3 ) L × S U ( 3 ) R Chiral Perturbation Theory. We review the determination of the $$B_K$$ B K parameter of neutral kaon mixing as well as the additional four B parameters that arise in theories of physics beyond the Standard Model. For the heavy-quark sector, we provide results for $$m_c$$ m c and $$m_b$$ m b as well as those for the decay constants, form factors, and mixing parameters of charmed and bottom mesons and baryons. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. We review the status of lattice determinations of the strong coupling constant $$\alpha _s$$ α s . We consider nucleon matrix elements, and review the determinations of the axial, scalar and tensor bilinears, both isovector and flavor diagonal. Finally, in this review we have added a new section reviewing determinations of scale-setting quantities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Software For Advanced Large-scale Analysis Of Magnetic Confinement For Numerical Design, Engineering & Research (salamander)

As magnetic confinement fusion energy gains traction internationally to enable abundant energy production, designing components for fusion systems is a pressing challenge. During the planned lifetime of a fusion device, components evolve in extreme environments and must withstand large, repeated thermal loads and bombardment by 14 MeV neutrons, plasma ions, and neutral particles (deuterium, tritium, and helium), corrosive conditions, etc. All these physical processes take place simultaneously, interact in intricate ways, and impose important constraints that can affect performance. Experimental data is rare and costly to obtain, making design particularly challenging. Predictive computational frameworks must be an integral part of an accelerated and cost-effective design process by modeling fusion system performance in simulated environments. To better understand component degradation and operational impacts on their performance, the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER) is designed as an open-source, fully integrated, multiphysics, multiscale, NQA-1 compliant framework facilitating 3D, high-fidelity fusion system modeling. To that end, SALAMANDER is a MOOSE-based framework, and therefore leverages MOOSE upstream libraries such as PETSc and libMesh to deliver sophisticated finite element, finite volume, and nonlinear solver technology for fusion energy simulations. SALAMANDER couples MOOSE physics module capabilities—such as thermal hydraulics, heat conduction, Navier-Stokes, and thermomechanics—with tritium transport via TMAP8, neutronics via Cardinal, and nascent particle-in-cell capabilities. Direct simulation Monte Carlo methods will be used to address neutral transport near the walls. By coupling all these physics in an integrated application, SALAMANDER will enable high-fidelity modeling of irradiation levels and plasma exposure conditions of plasma facing components and their impact on heat and tritium distributions, as well as the resulting mechanical constraints experienced by the plasma facing components and performance of blanket systems. Furthermore, SALAMANDER will be particularly suited for engineering studies thanks to the stochastic tool module readily available in MOOSE, allowing for extended uncertainty quantification and risk analysis studies. It is also able to use computer-aided design (CAD) meshes to model complex geometries, which is indispensable for fusion systems. SALAMANDER therefore supports design, safety, engineering, and research projects for magnetic confinement fusion systems

Simon, Pierre-Clement [Idaho National Laboratory (↗
DOE OSTI2025

Computational Investigation of the Chemical Bond between An(III) Ions and Soft-Donor Ligands

The chemical bonding of actinide ions with arene and borohydride ligands is explored via quantum chemical methods to understand how the transuranium elements interact with softdonor ligands. Specifically, the [An(C 6 Me 6 )(BH 4 ) 3 ] complexes (An = U, Np, and Pu) and their reduced congeners are studied. Density functional theory (DFT) shows that the metal–ligand interactions in the neutral complexes are governed by electrostatic interactions. Both DFT and complete active space (CASSCF) results show that as one moves from U to Pu, the 5f-orbitals are stabilized leading to a poorer energy match with the ligand orbitals. This contributes to progressively weaker metal-arene and metal-borohydride interactions across the series due to a decrease in energy-driven covalency. A reduction in orbital contributions to bonding is obtained for the transuranium-arene interactions as well. Upon reduction, the arene is reduced, forming a δ-bond. This causes the An–arene distances to contract by 0.1–0.2 Å compared to the neutral complexes. The ground state is assigned as the intermediate-spin state where the arene radical is antiferromagnetically coupled to the metal-centered f-electrons in Np and Pu. On the other hand, the ferromagnetically and antiferromagnetically coupled states are close in energy in the uranium complex, but do not mix when spin– orbit coupling is included using a state-interaction approach (SO-CASPT2). The population of the CASSCF δ*-antibonding natural orbital increases from U to Pu consistent with the increased An–arene distances, weaker interactions, and decreasing covalency across the series. Although the An–B distance increases by ca. 0.06 Å upon reduction, both the neutral and reduced species involve an An(III)–borohydride bond and as such are qualitatively similar. The Np complexes can be assigned to have slightly weaker bonding than the uranium analogs but are overall “uranium-like”. The Pu complexes are predicted to have less covalent contributions to bonding in both the Pu–arene and Pu–borohydride interactions; however, the Pu–arene interaction is predicted to be particularly weak.

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

Design and Engineering of LUPIN: A Test-Bed Radio-Frequency Ion Source for Enhanced Neutral Beam Injection on DIII-D

The Large, Uniform Plasma for Ionizing Neutrals (LUPIN) is a radio-frequency (RF) inductively coupled plasma (ICP) chamber for demonstrating plasma performance of an RF ICP positive ion source upgrade for the DIII-D neutral beam injection (NBI) system. LUPIN will be used to investigate ion source physics, including neutral gas dynamics, plasma density uniformity, interactions with Faraday shields, and power coupling to novel RF antenna designs. LUPIN has an RF generator capable of delivering 20 kW of power at 2 MHz, which is coupled into a cylindrical quartz vessel measuring 20 cm in length and 10 cm in radius. This configuration matches the power density requirements for a full-scale ion source. Target hydrogen and deuterium plasma densities exceeding 10 18 m -3 would relate to extracted ion current densities of 2100 A/m 2 for 10s. Vacuum conductance and gas flow calculations predict a maximum achievable neutral gas flow rate of 1675 Pa ⋅ L/s at 5 Pa of He, which mimics the gas flow of the DIII-D NBI system. Designs have been developed for an internal Faraday shield to mitigate heat flux and ion sputtering on the dielectric vessel. Thermomechanical finite element simulations demonstrated the Faraday shield design to be capable of withstanding anticipated heat loads from worst-case operation scenarios. Finally, results of upcoming experimental investigations on LUPIN will guide the design of a full-scale prototype for DIII-D integration.

Faraday shield↗
DOE OSTI2022

Dispersive optical systems for scalable Raman driving of hyperfine qubits

Hyperfine atomic states are among the most promising candidates for qubit encoding in quantum information processing. In atomic systems, hyperfine transitions are typically driven through a two-photon Raman process by a laser field which is amplitude modulated at the hyperfine qubit frequency. Here we introduce a method for generating amplitude modulation by phase modulating a laser and reflecting it from a highly dispersive optical element known as a chirped Bragg grating. This approach is passively stable, offers high efficiency, and is compatible with high-power laser sources, enabling large Rabi frequencies and improved quantum coherence. We benchmark this approach by globally driving an array of approximately 300 neutral 87 Rb atomic qubits trapped in optical tweezers and obtain Rabi frequencies of 2 MHz with photon-scattering error rates of less than 2×10 -4 per π pulse. This robust approach can be directly integrated with local addressing optics in both neutral atom and trapped ion systems to facilitate high-fidelity single-qubit operations for quantum information processing.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Kinetic Plasma Simulation in the MOOSE Framework: Verification of Electrostatic Particle In Cell Capabilities

In magnetic confinement nuclear fusion reactors, the interaction between the plasma edge and plasma facing components is extremely important. At the plasma edge, a kinetic representation such as particle-in-cell (rather than a fluid representation) is required to accurately capture the plasma behavior. General purpose particle-in-cell plasma simulation capabilities have been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. This new capability is a part of the development of a new MOOSE-based framework for modeling plasma facing components, the Fusion ENergy Integrated multiphys-X (FENIX) framework. In this work, the verification of foundational particle-in-cell capabilities in FENIX is presented. This new plasma simulation capability has three main components: moving particles in discrete steps on the finite element mesh, mapping charge density from the particle's location to the finite element mesh, and solving for the electrostatic potential based on the charge density mapped from particles to the mesh. In this paper, simple verification problems demonstrating each of these new capabilities are presented, and future work includes electromagnetic capabilities and Monte Carlo collisions with neutral gas particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A Comparison of Three Neodymium Atomic Data Sets for Kilonova Modeling

We examine the impact of input neodymium (Nd) atomic data on the light curves and spectra of kilonovae (KNe), probing the sensitivity of kilonova observables to the atomic physics of this important lanthanide element. We use the SuperNu Monte Carlo radiative transfer code, simulating a simple semianalytic 1D kilonova (KN) with a pure Nd atmosphere, fixing the radiative transfer method while using input atomic data generated by three different codes: the LANL suite of atomic physics codes, HULLAC, and Autostructure. We see that the choice of atomic data significantly shapes the resulting light curves and spectra. Peak bolometric luminosities differ by a ratio of nearly 1.5 between HULLAC/Autostructure and LANL data sets. Moreover, we observe significant near- to mid-IR differences in the structure of the spectra. We specifically attribute these differences to the choice of atomic data for neutral Nd I. Many of the results here have been adapted from a presentation at “Radiative Transfer and Atomic Physics of Kilonovae” in Stockholm, 2023. We additionally present a LANL data set with energies calibrated to available values in the NIST Atomic Spectra Database, and demonstrate that this calibration also significantly affects IR spectral structure at late time. The substantial differences in KN observables that arise from tuning the atomic data of just one lanthanide element highlight the special attention that must be paid to atomic physics uncertainties when modeling KNe, from AT2017gfo to beyond.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A combined experimental and modelling study of granite hydrothermal alteration

Geochemical reactions can induce significant changes of rock reservoir porosity and permeability via mineral dissolution and precipitation processes, affecting the long-term fluid behaviour within various geological systems. Here, the understanding and quantification of these reactions rely on field and experimental studies and on the predictions of reactive transport models. The present study was aimed at assessing the extent to which current geochemical models integrating available mineral dissolution/precipitation rate equations can reproduce the experimental data obtained from 4 to 17-day long hydrothermal alteration experiments of a muscovite-biotite granite and, thus, help provide an accurate description of the evolution of geothermal systems within granitic reservoirs. The experiments were conducted at a constant temperature of 180 °C and over an aqueous fluid pH range of 2 to 8.5, using both mixed-flow and static batch reactors. Modelled major element (K, Al, Si, Ca, and Mg) concentrations were generally in satisfactory agreement with the corresponding measured elemental fluxes – the differences between modelled and experimental values were generally within a factor of 5 – and the predicted identity and mass of formed secondary phases were consistent with the microscopic observations of the reacted solids. However, larger differences between measured and modelled element concentrations were observed when significant amounts of secondary phases formed, notably at pH 2 to 3, and for longer-term batch experiments. Much of this concentration difference stems from the underestimation of the amounts of Al-phases formed at acid to near-neutral pH. Although an idealized rock composition was considered, the observed mismatch between model calculations and experimental data can be attributed to inadequate mineral precipitation reaction rates and a poor description of reactive surface areas in existing geochemical modelling codes. More accurate quantification of precipitation kinetics, including nucleation and growth, and improved descriptions of the temporal change of mineral surface area would enhance the predictive capabilities of reactive transport models and benefit, particularly, the efforts aimed at increasing the sustainability of EGS reservoirs.

58 GEOSCIENCES↗
DOE OSTI2023

Heterogenization of Homogeneous Ruthenium(II) Catalysts for Carbon-Neutral Dehydrogenation of Polyalcohols

Liquid organic hydrogen carrier (LOHC) systems are an excellent alternative to pressurized gas and liquid hydrogen storage technologies due to their high volumetric storage capacities and straightforward adaptation to existing infrastructure. Homogeneous catalysts are promising for the selective and reversible release of hydrogen from LOHC. However, separation from product mixtures and recycling inhibit their use, particularly when comprised of costly low-abundance elements, motivating the development of heterogeneous versions that are more easily recovered and reused. Here, we describe two methods for the heterogenization of molecular Ru catalysts that efficiently dehydrogenate the polyalcohol LOHCs ethylene glycol (EG) and 1,2-propanediol (1,2-PDO). The heterogeneous versions of these catalysts maintain catalytic activity for hydrogen production comparable to the homogeneous complexes, with up to 81% conversion and 99% selectivity. Further, DFT modeling indicates mechanistic similarities for the dehydrogenations of EG and 1,2-PDO, with the rate-limiting steps associated with protonation of the Ru–H bond to form H 2 and the alkoxide species coordinated to Ru(II), followed by β-hydride elimination to regenerate the Ru–H bond. Overall, the data suggest these heterogenized molecular catalysts have potential for practical use in polyalcohol-based LOHC systems.

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

The deconstruction of flavor in the privately democratic Higgs sector

The Standard Model (SM) of particle physics fails to explain the observed hierarchy in fermion masses or the origin of fermion-flavor structure. We construct a model to explain these observations in the quark sector. We introduce a spectrum of new particles consisting of six of each — massive singlet vector-like quarks (VLQs), singlet scalars, and SU(2)-doublet scalars. SM quark masses are generated when the neutral components of the SU(2)-doublet scalars acquire non-zero vacuum expectation values (VEVs). We impose global symmetries to ensure that Yukawa couplings stay roughly flavor diagonal and democratic (of the same order), as well as to suppress tree-level flavor-changing neutral currents. Quark-mass hierarchy then follows from a hierarchy in scalar VEVs. The singlet scalars also acquire weak-scale VEVs. Together with the VLQs, they act as messengers between different generations of quarks in the SM. These messenger particles are responsible for generating the elements of the Cabibbo-Kobayashi-Masakawa (CKM) matrix which depend on the ratios of the singlet VEVs and VLQ masses. Constructed this way, the CKM matrix is found to be independent of the SM fermion masses. Using the measured values of the CKM matrix elements and assuming order-one couplings, we derive constraints on the masses of the VLQs and discuss prospects for probing our model in the near future.

CKM Parameters↗
DOE OSTI2025

Experimental Setup to Investigate Hydrogen Isotope Retention on Powder Samples as Slag/Dust Proxies for Advanced Fusion Reactors

Slag/dust formation is one of the critical issues regarding tritium retention in advanced fusion reactors. We upgraded an ultrahigh vacuum device, the sample exposure station (SES), to perform temperature-programmed desorption (TPD), also known as thermal desorption spectroscopy, to investigate deuterium (D) retention in powder samples that are utilized as proxies for slag/dust particles formed in fusion reactors. TPD analysis, calibrated for the desorption rate, was successfully performed on commercially available crystalline/amorphous boron (B) powder after D 2 neutral gas or D + ion exposure. We calibrated the B powder temperature measurement for TPD analysis utilizing a D 2 desorption peak from B-D bonding known for its desorption temperature ~700 K. We introduce a separation procedure of observed D TPD signals, convoluted for corresponding D retention pathways, that was necessary for TPD analysis. This is the first report on the retention of D 2 neutral gas to B, and this observation was enabled by utilizing the B powder sample, which has a larger surface area than films commonly used in laboratory experiments. In conclusion, the upgraded SES system will enable us to investigate hydrogen isotope behaviors for powder/dust samples of various elemental compositions, including metals.

Plasma-facing component debris↗
DOE OSTI2021

Commissioning the H i Observing Mode of the Beam Former for the Cryogenically Cooled Focal L-band Array for the GBT (FLAG)

We present the results of commissioning observations for a new digital beam-forming back end for the Focal plane L-band Array for the Robert C. Byrd Green Bank Telescope (FLAG), a cryogenically cooled Phased Array Feed (PAF) with the lowest measured T {sub sys}/η of any PAF outfitted on a radio telescope to date. We describe the custom software used to apply beam-forming weights to the raw element covariances to create research-quality spectral-line images for the new fine-channel mode, study the stability of the beam weights over time, characterize FLAG’s sensitivity over a frequency range of 150 MHz, and compare the measured noise properties and observed distribution of neutral hydrogen emission from several extragalactic and Galactic sources with data obtained with the current single-pixel L-band receiver. These commissioning runs establish FLAG as the preeminent PAF receiver currently available for spectral-line observations on the world’s major radio telescopes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Design workflow of a symmetric traveling wave antenna for fast ion production on DD tokamaks

Initial computational plasma physics scoping and a finite element method antenna modeling design workflow for a symmetric center-fed high-field side high harmonic fast wave traveling wave array (TWA) antenna are reported here. The TWA is designed to generate a test population of fast deuterium ions in an existing D–D tokamak by heating neutral beam deuterium ions, accelerating them from 80 keV to several hundred keV. The resulting fast particles are tailored to mimic key reactor energetic particle parameters with regards to exciting Alfven eigenmode instabilities, allowing for a D–D tokamak like DIII-D or ASDEX-U to replicate reactor-relevant conditions experimentally. Initial scenario scoping for high single-pass absorption as well as good preferential fast ion damping relative to electron damping was completed using the ray-tracing/Fokker–Planck codes GENRAY and CQL3D. Python RF network analysis packages were used to create a custom TWA optimization tool to inform a COMSOL flat antenna design, and Petra-M was used to study cold plasma effects. The TWA produced by this workflow has several novel features when compared to previous TWA studies, including symmetric center feeding, and passive end straps for image current cancellation for reduced impurity production. We show here that the antenna design workflow can readily produce TWA antennas optimized for reflection coefficient, image current cancellation, and launched power spectrum shape; and that a population of fast ions can be generated in the correct region of parameter space, warranting future more detailed studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell

Development of particle-in-cell (PIC) methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While this development is indeed a significant advance, it should be noted that implicit FEM transient solvers support an underlying null space that corresponds to a gradient of a scalar potential ∇Φ(r) (or t∇Φ(r) in the case of wave equation solvers). While explicit schemes do not suffer from this drawback, they are only conditionally stable with time step sizes that are mesh dependent and very small. Furthermore, the null space produces spurious charge that can corrupt the desired physics of a PIC simulation. The way to overcome this bottleneck, and indeed, satisfy all four Maxwell's equation is to use a quasi-Helmholtz formulation on a tessellation. In the re-formulation presented, we strictly satisfy the equation of continuity and Gauss' laws for both the electric and magnetic flux densities. Results illustrating the efficacy of this scheme will be demonstrated by analyzing non-neutral beams with space-charge and the adiabatic expansion of a neutral plasma with realistic parameters (Debye length and real mass ratios).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗