SEARCH · Engineering Papers

Results for “spectral methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Two-tooth bosonic quantum comb for temporal-correlation sensing

We introduce a two-tooth bosonic quantum comb that captures the sequential interactions between a thermal absorber and a long-lived coherent probe. The comb provides a causal, multi-time description of coherence transport, tracking how the probe records both instantaneous fluctuations and their temporal correlations. Using a process-tensor formulation, we derive closed form expressions showing that interference between the two interaction windows generates a non-monotonic memory response that reflects a fundamental competition between the absorbers thermal population and its dynamical correlations. By sweeping the temporal separation between the interaction windows, the probe directly samples the absorbers population correlator, enabling bosonic noise spectroscopy that discriminates Markovian temperature noise from slow or spectrally structured fluctuations. The approach is readily compatible with circuit-QED platforms and offers a general method for probing fluctuating bosonic environments.

Zhu, Shaojiang [Fermilab]↗

DOE OSTI2024

CNN-Based Phase Fault Classification in Real and Simulated Power Systems Data

This study proposes a convolutional neural network (CNN)–based two-step phase fault detection and identification method to classify anomalies in the power grid signal. Specifically, the first step checks the fault’s existence and determines the need for the second step. Subsequently, in the case of anomalies in the power grid signal, the second step identifies the type of fault, including line-to-line, single-line-to-ground, double-line-to-ground, and triple-line. Accordingly, the CNN architecture is both designed for the classification layers and trained with simulated data. To provide maximum prediction accuracy with minimum processing time, this study investigates the combinations of various feature extraction (FE) techniques, such as fast Fourier transform (FFT), amplitude and phase (AP), auto-correlation function, power spectral density, and wavelet transform (WT). Consequently, simulated and real-world results demonstrate that the proposed two-step method outperforms conventional one-step techniques, with the best performance obtained by using the combination of AP-AP, AP-WT, FFT-AP, and FFT-WT–based FE methods.

Alaca, Ozgur↗

DOE OSTI2025

Sequential spectral line analysis for accurate density and temperature diagnosis of laboratory opacity measurements

The accuracy of iron opacity calculated in stellar interiors has been questioned since the discovery of the “solar problem” and the discrepancies between the measured and modeled iron opacity reported in 2015. Experimental opacity benchmarks require accurate temperature and density measurements, which were inferred by analyzing tracer magnesium spectra in those experiments. Could the observed discrepancy be explained by insufficient accuracy in the inferred temperature, density, and their uncertainties? Previous analyses may have yielded biased results due to three limitations: (1) simultaneous multi-line fitting, (2) approximations in line-shape models, and (3) exclusion of certain spectral lines due to insufficient background characterization. Notably, the first issue is a common concern for many inversion methods, including Bayesian inferences. We present a refined analysis method that overcomes these limitations, applied to three categories of iron opacity experiments (Anchor 1, 2, and 3). In particular, the sequential fitting method yields unbiased results with more realistic uncertainties by accounting for line inconsistencies in the parameter uncertainties. The average electron temperature and density values are 162 ± 6 eV and (7.0 ± 1.9) × 10 21 cm −3 for six Anchor 1 experiments, 189 ± 7 eV and (3.4 ± 0.3) × 10 22 cm −3 for 21 Anchor 2 experiments, and 201 ± 6 eV and (4.8 ± 1.1) × 10 22 cm −3 for nine Anchor 3 experiments. These results show ∼4% temperature and ∼20% density reproducibility over a decade, which also aligns with the inferred parameter uncertainties. In conclusion, the resulting temperature and density uncertainties lead to a quasi-continuum iron opacity variation of ±4%–7% for wavelengths below 9.5 Å, which is insufficient to explain the significant model-data discrepancies reported in 2015.

Absorption spectroscopy↗

DOE OSTI2026

Robust Automatic EXAFS First-Shell Fits

Extended X-ray absorption fine structure (EXAFS) is a widely used technique for atomic structure determination. Fourier transformation connects EXAFS in k space and R space. However, determining the appropriate k-range for the transformation can be challenging, but critical for the first-shell fit. In this study, we present an automatic method to determine the k-range using the Larch package and a Python program. The first step is to estimate spectral noise across a series of k-ranges with a fixed minimum value and identify the optimal maximum value in the k-range (k max ). The k max is determined by an empirical noise threshold that marks the point where the noise level in the Fourier transformed spectrum changes dramatically. Using the obtained k max value, the first shell is modeled to determine the minimum k value (k min ) by optimizing the background function through alignment of the spectrum with theory. The optimal k min corresponds to the point of the minimum R-factor, which quantifies the difference between the experimental and fitted spectrum. Our method was tested on various typical datasets and yielded suitable k-ranges for Fourier transformation and accurate first-shell fits. This approach helps avoid unreliable, irreproducible data analysis, especially for noisy data from diluted samples, and enables robust automatic first-shell EXAFS fitting.

EXAFS analysis↗

DOE OSTI2022

Effective Field Theory of Random Quantum Circuits

Quantum circuits have been widely used as a platform to simulate generic quantum many-body systems. In particular, random quantum circuits provide a means to probe universal features of many-body quantum chaos and ergodicity. Some such features have already been experimentally demonstrated in noisy intermediate-scale quantum (NISQ) devices. On the theory side, properties of random quantum circuits have been studied on a case-by-case basis and for certain specific systems, and a hallmark of quantum chaos—universal Wigner–Dyson level statistics—has been derived. This work develops an effective field theory for a large class of random quantum circuits. The theory has the form of a replica sigma model and is similar to the low-energy approach to diffusion in disordered systems. The method is used to explicitly derive the universal random matrix behavior of a large family of random circuits. In particular, we rederive the Wigner–Dyson spectral statistics of the brickwork circuit model by Chan, De Luca, and Chalker [Phys. Rev. X 8, 041019 (2018)] and show within the same calculation that its various permutations and higher-dimensional generalizations preserve the universal level statistics. Finally, we use the replica sigma model framework to rederive the Weingarten calculus, which is a method of evaluating integrals of polynomials of matrix elements with respect to the Haar measure over compact groups and has many applications in the study of quantum circuits. The effective field theory derived here provides both a method to quantitatively characterize the quantum dynamics of random Floquet systems (e.g., calculating operator and entanglement spreading) and a path to understanding the general fundamental mechanism behind quantum chaos and thermalization in these systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Benzyl trichloroacetimidates as derivatizing agents for phosphonic acids related to nerve agents by EI-GC-MS during OPCW proficiency test scenarios

Abstract The use of benzyl trichloroacetimidates for the benzylation of phosphonic acid nerve agent markers under neutral, basic, and slightly acidic conditions is presented. The benzyl-derived phosphonic acids were detected and analyzed by Electron Ionization Gas Chromatography–Mass Spectrometry (EI-GC–MS). The phosphonic acids used in this work included ethyl-, cyclohexyl- and pinacolyl methylphosphonic acid, first pass hydrolysis products from the nerve agents ethyl N -2-diisopropylaminoethyl methylphosphonothiolate (VX), cyclosarin (GF) and soman (GD) respectively. Optimization of reaction parameters for the benzylation included reaction time and solvent, temperature and the effect of the absence or presence of catalytic acid. The optimized conditions for the derivatization of the phosphonic acids specifically for their benzylation, included neutral as well as catalytic acid (< 5 mol%) and benzyl 2,2,2-trichloroacetimidate in excess coupled to heating the mixture to 60 °C in acetonitrile for 4 h. While the neutral conditions for the method proved to be efficient for the preparation of the p -methoxybenzyl esters of the phosphonic acids, the acid-catalyzed process appeared to provide much lower yields of the products relative to its benzyl counterpart. The method’s efficiency was tested in the successful derivatization and identification of pinacolyl methylphosphonic acid (PMPA) as its benzyl ester when present at a concentration of ~ 5 μg/g in a soil matrix featured in the Organisation for the Prohibition of Chemical Weapons (OPCW) 44th proficiency test (PT). Additionally, the protocol was used in the detection and identification of PMPA when spiked at ~ 10 μg/mL concentration in a fatty acid-rich liquid matrix featured during the 38 th OPCW-PT. The benzyl derivative of PMPA was partially corroborated with the instrument's internal NIST spectral library and the OPCW central analytical database (OCAD v.21_2019) but unambiguously identified through comparison with a synthesized authentic standard. The method’s MDL (LOD) values for the benzyl and the p -methoxybenzyl pinacolyl methylphosphonic acids were determined to be 35 and 63 ng/mL respectively, while the method’s Limit of Quantitation (LOQ) was determined to be 104 and 189 ng/mL respectively in the OPCW-PT soil matrix evaluated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Nuclear-Electronic Orbital Ehrenfest Dynamics

The recently developed real-time nuclear-electronic orbital (RT-NEO) approach provides an elegant framework for treating electrons and selected nuclei, typically pro- tons, quantum mechanically in nonequilibrium dynamical processes. However, the RT-NEO approach neglects the motion of the other nuclei, preventing a complete de- scription of the coupled nuclear-electronic dynamics and spectroscopy. In this work, we describe the dynamical interaction between the other nuclei and the electron-proton subsystem with the mixed quantum-classical Ehrenfest dynamics method. The NEO- Ehrenfest approach propagates the electrons and quantum protons in a time-dependent variational framework, while the remaining nuclei move classically on the correspond- ing average electron-proton vibronic surface. This approach includes the non-Born- Oppenheimer eects between the electrons and quantum protons with RT-NEO and between the classical nuclei and the electron-proton subsystem with Ehrenfest dynam- ics. Spectral features for vibrational modes involving both quantum and classical nuclei are resolved from the time-dependent dipole moments. Moreover, nuclear quantum ef- fects are included directly in nonadiabatic dynamics simulations. The excited state intramolecular proton transfer in the o-hydroxybenzaldehyde molecule is shown to be faster and to exhibit a larger kinetic isotope eect with NEO-Ehrenfest compared to classical Ehrenfest dynamics. This work shows that the NEO-Ehrenfest method is a powerful tool to study dynamical processes with coupled electronic and nuclear degrees of freedom.

Zhao, Luning↗

DOE OSTI2022

Extending water vapor measurement capability of photon-limited differential absorption lidars through simultaneous denoising and inversion

Abstract. The micropulse differential absorption lidar (MPD) was developed at Montana State University (MSU) and the National Center for Atmospheric Research (NCAR) to perform range-resolved water vapor (WV) measurements using low-power lasers and photon-counting detectors. The MPD has proven to produce accurate WV measurements up to 6 km altitude. However, the MPD's ability to produce accurate higher-altitude WV measurements is impeded by the current standard differential absorption lidar (DIAL) retrieval methods. These methods are built upon a fundamental methodology that algebraically solves for the WV using the MPD forward models and noisy observations, which exacerbates any random noise in the lidar observations. The work in this paper introduces the adapted Poisson total variation (PTV) specifically for the MPD instrument. PTV was originally developed for a ground-based high spectral resolution lidar, and this paper reports on the adaptations that were required in order to apply PTV on MPD WV observations. The adapted PTV method, coined PTV-MPD, extends the maximum altitude of the MPD from 6 to 8 km and substantially increases the accuracy of the WV retrievals starting above 2 km. PTV-MPD achieves the improvement by simultaneously denoising the MPD noisy observations and inferring the WV by separating the random noise from the non-random WV. An analysis with 130 radiosonde (RS) comparisons shows that the relative root-mean-square difference (RRMSE) of WV measurements between RS and PTV-MPD exceeds 100 % between 6 and 8 km, whereas the RRMSE between RS and the standard method exceeds 100 % near 3 km. In addition, we show that by employing PTV-MPD, the MPD is able to extend its useful range of WV estimates beyond that of the ARM Southern Great Plains Raman lidar (RRMSE exceeding 100 % between 3 and 4 km); the Raman lidar has a power-aperture product 500 times greater than that of the MPD.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Entanglement Enabled Intensity Interferometry of different wavelengths of light

Highlights: • We use quantum entanglement to erase the color of photons for a new imaging method. • Can achieve higher resolution imaging by accessing hidden phase information. • We detail practical experimental designs to implement our theoretical protocols. We propose methods to perform intensity interferometry of photons having two different wavelengths. Distinguishable particles typically cannot interfere with each other, but we overcome that obstacle by processing the particles via entanglement and projection so that they lead to the same final state at the detection apparatus. Specifically, we discuss how quasi-phase-matched nonlinear crystals can be used to convert a quantum superposition of light of different wavelengths onto a common wavelength, while preserving the phase information essential for their meaningful interference. We thereby gain access to a host of new observables, which can probe subtle frequency correlations and entanglement. Further, we generalize the van Cittert–Zernike formula for the intensity interferometry of extended sources, demonstrate how our proposal supports enhanced resolution of sources with different spectral character, and suggest potential applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Spectral characterization of flash and high flux x-ray radiographic sources with a magnetic Compton spectrometer

In this article, we present a new analysis method applied to revitalize permanent magnet Compton spectrometers used to measure photon energy spectra in the MeV range. The inversion of the measured electron distribution to determine the original photon distribution is achieved via a method of consistent coupled radiation transport and magnetic field mapping of the input photon spectra to the measured electron distribution. The method of linear least squares was used to perform the unfolding of the electron distribution to the initial photon spectra, without any assumptions made regarding the electron distribution. We present an application of this method to data from a nominal 19.4 MeV flash radiographic source (the first axis of the Dual Axis Radiographic Hydro-Test Facility) capable of generating 500 R @ 1 m in ~60 ns and a medical therapy source (a Scanditronix M22, Microtron) capable of variable energies with nominal endpoints of 6, 10, 15, and 20 MeV and an output of ~1000–2000 R/min @ 1 m. The results provide agreement between the modeled and unfolded experimentally measured photon spectra as quantified by statistical tests, from 1.5 to 20 MeV. Experimental results are presented as well as a discussion of the novel MCNP6-based simulations and methods for reconstruction of the spectra.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Spectral signatures of transverse optical modes in semiconductor nanowires

In order to build nanophotonic devices, it is important to understand and ultimately control the optical mode structure within potential components such as nanoscale waveguides. However, experimental characterization of such modes in the optical regime is difficult due to the nanoscale dimensions of such components and the perturbations that would be induced by a near-field probe. Here, we demonstrate a probe-free, far-field method to characterize the optical modes within GaN nanowires (NWs) based on a novel off-axis scanning confocal microscope system. Using this microscope, we observe spectral signatures resulting from lateral leakage of waveguide modes when they exceed their respective cutoff limits. We identify these modes within hyperspectral images using an analytical model coupled with finite element simulations. The model can also be used to predict the spectral signatures for given geometrical and dielectric parameters, which enabled us to deduce the transverse dimension of the NW from hyperspectral images with an accuracy of ∼ <#comment/> 30 nm .

Wang, Maoji↗

DOE OSTI2025

Detection of multi-modal Doppler spectra – Part 1: Establishing characteristic signals in radar moment data

Vertically pointing millimeter-wavelength radars provide a wealth of information about cloud and precipitation particle properties. Doppler spectral data can inform on how particles of varying vertical velocities contribute to the total backscattered power observed. It is more computationally cost effective to process moment data instead of spectra data, but doing so leaves valuable information on the cutting room floor. To confidently identify a multi-modal spectra event, in which two or more modes are present within a layer, Doppler spectral data are essential. This means long-term identification of layers featuring multi-modal spectra can be cost prohibitive. To address this, we explore three multi-modal spectra cases from winter precipitation events to determine characteristic signatures of these layers in the moment data averaged over short time periods (∼ 145 s) and explore how these layers differ from the rest of the vertical profiles. We find that the mean spectrum width and the standard deviation of mean Doppler velocity can be used to determine whether or not a layer is multi-modal. In particular, multi-modal layers in mixed-phase and ice clouds feature larger mean spectrum width (exceeding 0.17 m s −1 ) and smaller standard deviation of the mean Doppler velocity (below 0.1 m s −1 ). In Part 1 of this study, the identification criteria and methods are described. In Part 2 (Wugofski and Kumjian, 2025), we perform a verification of the method for three years of vertically pointing radar data, and explore the meteorological conditions associated with identified multi-modal spectral events.

Wugofski, Sarah [Pennsylvania State Univ., Univers↗

DOE OSTI2022

Comparing the shock sensitivity of insensitive energetic materials

We present a tabletop method to study the shock sensitivity of plastic-bonded explosives that are considered shock insensitive using high dynamic range optical emission spectroscopy with laser-launched km/s flyer plate impacts (2, 3, and 4 km/s), which measures the spectral radiance (the emission spectrum vs a calibrated intensity standard) every 0.8 ns in the nanosecond and microsecond regimes. The explosives were TATB (1,3,5 trinitro, 2,4,6 triamino benzene), FOX-7 (1,1-diamino-2,2-dinitroethylene), LLM-105 (2,6-diamino-3,5-dinitropyrazine-1-oxide), and NTO (nitrotriazolone), all with 20% Sylgard 182 binder. Time resolving emission from shocked explosives produces a unique fingerprint at each shock pressure, and the emission temporal profile can be used to understand each material's shock response and the underlying mechanisms within the framework of the well-known shock ignition and growth models. Ignition was characterized by the emission intensity during nanosecond hot spot ignition. The growth process, which occurs after the shock has unloaded, was characterized by the time derivative of microsecond emission intensity. Singular-value decomposition was used to determine a characteristic spectral radiance over the entire time range, and this radiance was fit to a graybody model to determine a characteristic temperature Tc. Even though these insensitive explosives have similar shock sensitivities, our method can clearly distinguish their shock response. FOX-7 and NTO were the most sensitive, with FOX-7 having the greatest shock ignition sensitivity, growth rate, and Tc. LLM-105 was much less sensitive than FOX-7 and NTO, and TATB consistently showed the lowest sensitivity by any measure.

Akhtar, Meysam (ORCID:0000000243897837)↗

DOE OSTI2022

Rethinking the ill-posedness of the spectral function reconstruction — Why is it fundamentally hard and how Artificial Neural Networks can help

Reconstructing hadron spectral functions through Euclidean correlation functions are of the important missions in lattice QCD calculations. However, in a Källen–Lehmann (KL) spectral representation, the reconstruction is observed to be ill-posed in practice. It is usually ascribed to the fewer observation points compared to the number of points in the spectral function. In this paper, by solving the eigenvalue problem of continuous KL convolution, we show analytically that the ill-posedness of the inversion is fundamental and it exists even for continuous correlation functions. We discussed how to introduce regulators to alleviate the predicament, in which include the Artificial Neural Networks (ANNs) representations recently proposed by the Authors in another study. The uniqueness of solutions using ANNs representations is manifested analytically and validated numerically. Reconstructed spectral functions using different regularization schemes are also demonstrated, together with their eigen-mode decomposition. We observe that components with large eigenvalues can be reliably reconstructed by all methods, whereas those with low eigenvalues need to be constrained by regulators.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Applying Two-Dimensional Correlation Spectroscopy and Principal Component Analysis to Understand How Temperature Affects the Neptunium(V) Absorption Spectrum

The visible-near infrared (Vis-NIR) electronic absorption spectrum of neptunium(V) (NpO 2 + ) comprises numerous f-f electronic transitions with mostly undocumented temperature dependencies. The effect of temperature on the absorption spectrum of the pentavalent neptunyl dioxocation (NpO 2 + ) is an important factor to consider with spectrophotometric applications but has often been overlooked. Optical Vis-NIR absorption spectra (400–1700 nm) of Np(V) (0.017–0.89 M) in 1 M nitric acid were evaluated with varying temperatures (T = 10–80 °C). The intensity, position, and overall shape of the bands were sensitive to interactions with the solvent and coordination environment. Numerous temperature-induced isosbestic points were identified resulting from dynamic, overlapping peak shifts. Spectral variations were characterized using principal component analysis (PCA) and 2D correlation spectroscopy (COS). 2D-COS revealed that the absorption band near 1095 nm likely consists of two bands centered near 1087 and 1096 nm, which cannot be explained by current computational methods. 2D-COS analysis also provided an unambiguous assignment of unresolved peaks in the visible region for comparison with computational predictions. PCA was used to identify nonlinearity in the spectral response at elevated Np(V) concentrations ≥ 0.5 M. This unique experimental data and interpretation will foster a deeper understanding of the absorption spectra for complex actinyl ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Fourier-based methods for removing mesh anomalies from angle resolved photoemission spectra

Recent improvements to spatial resolution in angle-resolved photo-emission spectroscopy (ARPES) have made it common to perform measurements with a very brief dwell time, for the purpose of mapping the spectral function over large surface regions. However, rapid measurement modalities can suffer a grid-like intensity modulation due to a wire mesh that is typically placed in front of the ARPES detector to block stray electrons. Here, we explore Fourier-based methods that can effectively remove this artifact, and improve the quality of ARPES images obtained in rapid scanning modes. In conclusion, an open source software package is provided containing implementations of demonstrated algorithms.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Automatic Échelle Spectrograph Wavelength Calibration

Time-domain astronomy and the increasing number of exoplanet candidates call for reliable, robust, and automatic wavelength calibration. We present an algorithm for wavelength-calibrating échelle spectrographs that uses order-by-order extracted spectra and a list of laboratory wavelengths. Our approach is fully automatic and does not need the pixel locations of certain spectral features with which to anchor the wavelength solution, nor the true order number of each diffraction order. We use spectral features that are duplicated in adjacent orders to establish the scale-invariant component of the wavelength solution. We then match the central wavelengths of spectral features to laboratory wavelengths to establish the scale and higher-order components of the wavelength solution. We demonstrate our method on the four spectrographs of Las Cumbres Observatory’s Network of Robotic Échelle Spectrographs (NRES), on the High Accuracy Radial Velocity Planet Searcher (HARPS) spectrograph, and on synthetic data. We obtain a velocity-equivalent precision of ∼10 m s{sup −1} on NRES. We achieve ∼1 m s{sup −1} on HARPS, which agrees with the precision reported by the HARPS team. On synthetic data, we achieve the velocity precision set by Gaussian centroiding errors. Our algorithm likely holds for a wide range of spectrographs beyond the five presented here. We provide an open-source Python package, xwavecal, which outputs wavelength-calibrated spectra as well as the wavelengths of spectral features.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗