SEARCH · Search NASA

Results for “Charge density waves”

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.

414 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Anode-less Solid-State Li–S Batteries Enabled by Fe-Stabilized Polysulfides

Anode-less solid-state lithium-sulfur batteries (SSLSBs) with lithium sulfide (Li 2 S) as the cathode promise a high energy density and ease of manufacturing. However, Li 2 S is plagued by poor conductivity, sluggish activation kinetics, and a poor cycle life. Here, in this study, we report an FeCl 3 -activated Li 2 S (FLS) cathode with solid-state polysulfide intermediates generated through a redox reaction between FeCl 3 and Li 2 S. This strategy is shown to boost the electrical conductivity of Li 2 S by 7 orders of magnitude and lower the activation barrier. During cycling, Fe plays a significant role in stabilizing the highly active polysulfide species, contributing to the exceptional electrochemical performance. The FLS cathode achieves 80% capacity retention over 500 cycles with >99% Li 2 S utilization. Furthermore, a Li-metal-free (anode-less) full cell retained over 80% of its initial capacity after 240 cycles. This work underscores the promise of leveraging Fe-stabilized polysulfides in enabling high-energy, long-lasting, solid-state Li-S batteries.

25 ENERGY STORAGE↗

DOE OSTI2025

Interleaved bond frustration in a triangular lattice antiferromagnet

Frustration of long-range order via lattice geometry amplifies fluctuations and generates ground states that are highly sensitive to perturbations. Traditionally, geometric frustration is used to engineer unconventional magnetic states; however, the charge degree of freedom and bond order can be similarly frustrated. Finding materials that host both frustrated magnetic and bond networks holds promise for engineering structural and magnetic states with the potential of coupling to one another via either magnetic or strain fields. Here we identify an unusual instance of this coexistence in the triangular lattice antiferromagnets LnCd 3 P 3 (Ln = lanthanides). These compounds feature two-dimensional planes of unique trigonal planar CdP 3 units with an underlying bond instability that is frustrated via emergent kagome ice correlations. This bond instability is interleaved in between layers of frustrated magnetic moments. Furthermore, our results establish LnCd 3 P 3 as a rare materials class in which frustrated magnetism is embedded within a dopable semiconductor with a frustrated bond order instability.

Gomez Alvarado, S. J. [University of California, S↗

DOE OSTI2022

Efficient electrocatalytic valorization of chlorinated organic water pollutant to ethylene

Electrochemistry can provide an efficient and sustainable way to treat environmental waters polluted by chlorinated organic compounds. However, the electrochemical valorization of 1,2-dichloroethane (DCA) is currently challenged by the lack of a catalyst that can selectively convert DCA in aqueous solutions into ethylene. Here, in this work, we report a catalyst comprising cobalt phthalocyanine molecules assembled on multiwalled carbon nanotubes that can electrochemically decompose aqueous DCA with high current and energy efficiencies. Ethylene is produced at high rates with unprecedented ~100% Faradaic efficiency across wide electrode potential and reactant concentration ranges. Kinetic studies and density functional theory calculations reveal that the rate-determining step is the first C–Cl bond breaking, which does not involve protons—a key mechanistic feature that enables cobalt phthalocyanine/carbon nanotube to efficiently catalyse DCA dechlorination and suppress the hydrogen evolution reaction. The nanotubular structure of the catalyst enables us to shape it into a flow-through electrified membrane, which we have used to demonstrate >95% DCA removal from simulated water samples with environmentally relevant DCA and electrolyte concentrations.

Choi, Chungseok [Yale Univ., New Haven, CT (United↗

DOE OSTI2024

Coercivity enhancement in hematite/permalloy heterostructures across the Morin transition

Interfacial effects between antiferromagnetic (AFM) and ferromagnetic (FM) materials have long been a center of magnetism studies. Aside from the exchange bias occurring at the AFM/FM interface, controlling the coercivity is another significant topic in magnetic recordings. The coercivity of FM materials is often determined through varying grain size, alloy composition, density of defects, etc., which is set during material growth and offers limited room for modification after growth. Hematite (α - Fe 2 O 3 ) is an AFM material that undergoes a temperature-controlled spin-flip transition, the so-called Morin transition. This transition gives an extra degree of freedom making hematite an intriguing component to study the exchange coupling when interfaced with an FM material. Here, in this work, changes in the magnetic properties of soft magnetic permalloy (Ni 81 Fe 19 , or Py) thin films grown on hematite were studied across the Morin transition. Surprisingly, these samples showed a remarkable change in coercivity during the Morin transition. We attribute this effect to the magnetic domain mixture of hematite during the Morin transition. Our findings present a novel method of controlling the coercivity of plain ferromagnetic thin films.

Antiferromagnet↗

DOE OSTI2026

Flame stabilization in DME spray flames under engine-relevant conditions characterized by OH* chemiluminescence and formaldehyde laser-induced fluorescence

The transient and quasi-steady flame structures of Dimethyl Ether (DME) fuel sprays, produced by a single-hole injector (Spray D), were investigated using Planar Laser-Induced Fluorescence (PLIF) and chemiluminescence imaging in a constant-volume chamber under Engine Combustion Network (ECN) Spray A conditions (900 K ambient temperature, 60 bar ambient pressure, 1500 bar injection pressure, and 22.8 kg/m 3 ambient density). Low-temperature chemical reaction zones were visualized using formaldehyde (CH 2 O) PLIF with 355 nm excitation, while high-temperature flame regions were captured via chemiluminescence imaging of excited-state hydroxyl radicals (OH*). Both transient and quasi-steady flame structures clearly show the transition from CH 2 O to OH*, highlighting the progression from low- to high-temperature combustion, while the position of the flame is displaced for DME compared to reference hydrocarbon n-dodecane. Homogeneous reactor calculations with detailed chemistry and using adiabatic mixing for initial temperature show that CH 2 O peaks are significantly higher for DME at the same equivalence ratio, with a higher heat-release during the cool-flame regime with respect to the fuel heating value. Thus, the cool-flame dynamic as a precursor to high-temperature combustion and flame stabilization exhibit distinct behavior for DME relative to conventional hydrocarbons, and these phenomena are effectively resolved through the soot-free nature of DME and the high-speed, time-resolved diagnostics.

CH2O laser-induced fluorescence↗

DOE OSTI2026

Pinning ångström-size solid ionic channels for rare-earth element separation

High-purity rare-earth elements are essential for modern technologies, yet current solvent extraction processes are energy-intensive and environmentally harmful because of inadequate selectivity and ligand toxicity. Although combining size exclusion and binding affinity can improve lanthanide separation, the role of long-range confinement remains underexplored. Here we report lanthanide separation in aqueous systems using extremely confined manganese oxide solid ionic channels with optimized layer spacing. Different lanthanides induce distinct solid-state phase transformations in manganese oxide, creating a strong driving force for separation. Two lanthanide groups, differing by ~1.4 Å in spacing, were identified and confirmed to be stable by density functional theory. The narrower confinement of heavier Group II lanthanides improves cross-group separation by increasing the dehydration barrier for lighter Group I lanthanides without inducing strong binding. Here, we further developed a strategy to pin the confinement dimensions and enhance same-group separation, increasing enrichment factors for La–Nd and La–Pr pairs from 1.6 ± 0.1 and 1.5 ± 0.1 to 5.4 ± 0.1 and 4.2 ± 0.1, respectively.

Chemical engineering↗

DOE OSTI2026

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗

DOE OSTI2026

SOLPS-ITER modeling of a dedicated divertor for negative triangularity operation on DIII-D

The design of a new dedicated divertor for negative triangularity (NT) operation on DIII-D with neutral baffles and pumping is informed by SOLPS-ITER transport modeling. This dedicated NT divertor is the latest step in a progression of NT shapes with various divertor characteristics explored on DIII-D, including NT shapes at reduced triangularity and a campaign with stronger shaping that included new armored components on the outboard side. SOLPS simulations played a key role in these divertor designs. Interpretive simulations, using cross-field diffusivities constrained by experimental data in the NT Shelf shape were used to inform the design of the 2023 armor campaign components. A similar procedure used armor campaign data to predict conditions for the dedicated NT divertor. The predictive simulations were used to assess the divertor fluxes, detachment threshold, pumped flux, and neutral leakage. For the dedicated NT divertor, SOLPS simulations and two-point-modeling were used to show the relative impact of magnetic topology (mainly longer connection length) and divertor closure on the divertor conditions relative to the armor campaign. It is predicted that the dedicated NT divertor reaches detachment (measured by target ion flux rollover) at a lower upstream density (≈(1.75−−2.4)×1⁢019m−3) as compared to the armor campaign shape. For the preliminary design geometry, divertor closure reduces the neutral leakage by ≈10%. Parametric optimization indicating further ≈20%–60% improvement in the leakage flux and recycled flux crossing the pump entrance is possible for relatively minor changes to the divertor and baffle layout.

Lore, Jeremy [ORNL] (ORCID:000000029192465X)↗

DOE OSTI2026

Insulating ground state and 2−𝑘 magnetic structure of candidate Weyl hydrogen-atom K2⁢Mn3⁢(AsO4)3

The ideal Weyl "hydrogen-atom" semimetal exhibits only a single pair of Weyl nodes and no other trivial states at the Fermi energy. Such a material would be a panacea in the study of Weyl quasi-particles, allowing direct unambiguous observation of their topological properties. The alluaudite-like K2⁢Mn3⁢(AsO4)3 compound was recently proposed as such a material. Here, we use comprehensive experimental work and first-principle calculations to assess this prediction. We find K2⁢Mn3⁢(AsO4)3 crystallizes in the 𝐶⁢2/𝑐 symmetry with a quasi-one-dimensional Mn sublattice, growing as small needle-like crystals. Bulk property measurements reveal magnetic transitions at ≈8 and ≈4 K, which neutron scattering experiments show correspond to two distinct magnetic orders, first a partially ordered ferrimagnetic 𝐤𝟏=(0,0,0) structure at 8 K and a second transition of 𝐤𝟐=(1,0,0) at 4 K to a fully ordered state. Below the second transition, both ordering vectors are necessary to describe the complex magnetic structure with modulated spin magnitudes. Both of the best-fit magnetic structures in this work are found to break the symmetry necessary for the generation of Weyl nodes, though one of the magnetic structures allowed by 𝐤𝟏 does preserve this symmetry. However, the crystals are optically transparent and ellipsometry measurements reveal a large band gap, undermining expectations of semimetallic behavior. Density functional theory calculations predict an insulating antiferromagnetic ground state, in contrast to previous reports, and suggest potential frustration on the magnetic sublattice. Given the wide tunability of the alluaudite structure, we consider ways to push the system closer to a semimetallic state.

Taddei, Keith M [Argonne National Laboratory]↗

DOE OSTI2026

Manipulating Na/TM Ratio‐Driven Structural Heterogeneity of O3‐NaNi 1/3 Fe 1/3 Mn 1/3 O 2 Cathode for High‐Voltage Sodium‐Ion Batteries

The stability of O3-type NaNi 1/3 Fe 1/3 Mn 1/3 O 2 under high-voltage cycling is dictated by how synthesis encodes lattice strain and redox heterogeneity. Here, in this study, the role of Na:TM stoichiometry is systematically resolved by tuning the NaOH:precursor ratio during solid-state synthesis. The stoichiometric condition (Na:TM = 1.00) yields minimized microstrain, enabling uniform O3–P3 phase evolution and homogeneous multi-metal redox with preserved octahedral symmetry. In contrast, Na-excess compositions inherit disordered intermediates and heterogeneous distortion fields that trigger abrupt multiphase transitions and promote localized charge redistribution. In situ XRD captures the divergence in phase-transition pathways, TXM resolves particle-level redox heterogeneity, and XANES corroborates a stronger and more reversible Fe redox contribution at stoichiometry, shifting to diminished Fe participation and spatially inhomogeneous redox at higher Na content. These results establish Na:TM stoichiometry as a critical synthesis parameter controlling both structural coherence and redox stability. Electrochemically, the stoichiometric composition exhibits smooth voltage profiles with minimal polarization growth and retains nearly 80% of its initial capacity after 100 cycles even at an extended 4.2 V cutoff, whereas Na-excess compositions show significantly reduced initial coulombic efficiency and rapid voltage fade. Precise stoichiometric tuning provides a scalable route to defect-suppressed O3 frameworks, enabling structurally resilient, high-voltage sodium-layered cathodes.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Review on Preprocessing Strategies, Deactivation, Thermal Safety, and Future Perspectives in Lithium-Ion Battery Recycling

The rapid growth in the use of lithium-ion batteries (LIBs) in electric vehicles, consumer electronics, and renewable energy storage has made effective end-of-life management essential. Recycling LIBs is critical not only for resource recovery and environmental protection but also for ensuring safety and economic viability. This review focuses on the preprocessing technologies that precede typical recovery processes, including disassembly, sorting, discharging, electrolyte removal, dismantling, thermal treatment, separation, and flotation. These steps play a foundational role in determining the efficiency, safety, and environmental impact of LIB recycling. LIBs pose substantial fire and explosion risks due to residual charge, flammable electrolytes, and reactive materials. The conditions and successive progression of the exothermic reactions which lead to thermal runaway has been discussed. It also explores secure deactivation techniques such as external circuit discharge, saline immersion, and thermomechanical methods, alongside fire prevention strategies including the use of flame retardants, elimination of oxidants, and reduction of heat generation and accumulation. Challenges and future directions are outlined, highlighting the need for standardized designs, automation, and safer, more sustainable recycling infrastructure. Furthermore, this review is distinguished by its focused analysis of preprocessing and deactivation steps, with particular attention to the thermal safety engineering aspects of LIB recycling.

Battery deactivation↗

DOE OSTI2024

Photoelectrochemical Proton-Coupled Electron Transfer of TiO 2 Thin Films on Silicon

TiO 2 thin films are often used as protective layers on semiconductors for applications in photovoltaics, molecule–semiconductor hybrid photoelectrodes, and more. Experiments reported here show that TiO 2 thin films on silicon are electrochemically and photoelectrochemically reduced in buffered acetonitrile at potentials relevant to photoelectrocatalysis of CO 2 reduction, N 2 reduction, and H 2 evolution. On both n-type Si and irradiated p-type Si, TiO 2 reduction is proton-coupled with a 1e – :1H + stoichiometry, as demonstrated by the Nernstian dependence of the Ti 4+/3+ E 1/2 on the buffer pK a . Experiments were conducted with and without illumination, and a photovoltage of ∼0.6 V was observed across 20 orders of magnitude in proton activity. The 4 nm films are almost stoichiometrically reduced under mild conditions. The reduced films catalytically transfer protons and electrons to hydrogen atom acceptors, based on cyclic voltammogram, bulk electrolysis, and other mechanistic evidence. TiO 2 /Si thus has the potential to photoelectrochemically generate high-energy H atom carriers. Characterization of the TiO 2 films after reduction reveals restructuring with the formation of islands, rendering TiO 2 films as a potentially poor choice as protecting films or catalyst supports under reducing and protic conditions. Altogether, this work demonstrates that atomic layer deposition TiO 2 films on silicon photoelectrodes undergo both chemical and morphological changes upon application of potentials only modestly negative of RHE in these media. While the results should serve as a cautionary tale for researchers aiming to immobilize molecular monolayers on “protective” metal oxides, the robust proton-coupled electron transfer reactivity of the films introduces opportunities for the photoelectrochemical generation of reactive charge-carrying mediators.

Electrodes↗

DOE OSTI2026

How Proton Incorporation Reshapes Lattice Dynamics In BaSnO 3 ‐Type Proton Conductors

Proton conduction in acceptor-doped perovskites is fundamentally a vibronic process: mobile H + and D + do not move independently, but dynamically co-vibrate with the surrounding oxygen–metal framework. Direct experimental evidence for this behavior is presented using in situ 119⁢ 𝑆⁢𝑛 nuclear resonance vibrational spectroscopy (NRVS) on hydrated, deuterated, and dry 𝐵⁢𝑎⁢𝑆⁢𝑛 1−𝑥⁢ 𝑌 𝑥 ⁢𝑂 3−𝛿 . Hydration induces systematic redistributions in the Sn-projected phonon density of states (PDOS), including an upshift of the first spectral moment by about 0.4 meV, indicating a stiffening of the extended Sn–O network. H/D isotopic substitution leaves the Sn-projected PDOS largely unchanged, with only subtle isotope-dependent spectral reweighting, demonstrating that protonic degrees of freedom are not localized oscillators but are embedded in collective lattice modes. These results are rationalized using a classical coupled proton–phonon oscillator model that links the observed PDOS variations to changes in effective force constants and vibrational mass terms. The model captures how H + and D + participate in cooperative lattice dynamics rather than forming isolated OH/OD entities. Overall, NRVS probes proton–lattice coupling in ceramic proton conductors and quantitatively describes how protonic defects modulate host lattice dynamics to enable phonon-assisted long-range proton transport.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Giant Exfoliation Induced Magnetic Coercivity in Fe 3 GaTe 2

Permanent magnets with strong anisotropy and high coercivity underpin modern information and energy technologies, yet rare-earth-free alternatives remain limited. Here, we show that thickness engineering via mechanical exfoliation induces hard magnetic behavior in the van der Waals ferromagnet Fe 3 GaTe 2 . Bulk crystals exhibit Curie temperatures above 350 K but negligible room-temperature coercivity. When thinned below ∼100 nm, the coercive field is dramatically enhanced, reaching nearly 1 T at room temperature for in-plane fields—comparable to conventional hard magnets. Micromagnetic analysis reveals a crossover in magnetization reversal from domain-mediated processes in bulk samples to quasi-coherent rotation in thin flakes, driven by increased effective anisotropy and suppressed domain formation. This thickness-dependent transition enables tuning of magnetic hardness without chemical modification. Combined with high saturation magnetization and robust room-temperature performance, Fe 3 GaTe 2 emerges as a promising rare-earth-free material for spintronic applications. Its layered structure further allows integration into van der Waals heterostructures, where large in-plane coercivity can stabilize magnetic states against perturbations and interlayer coupling, offering potential for high-density nonvolatile memory and domain-wall-based devices.

36 MATERIALS SCIENCE↗

DOE OSTI1991

Recommendations for Establishing Local Rural Electrification Programs Using Photovoltaic Systems [Recomendaciones Para la Implantacidn de Programas Locales de Electrificacidn Rural con Sistemas Fotovoltaicos]

Photovoltaic (PV) technology is becoming one the best options for supplying electricity to rural communities far from electrical networks and which have a small and disperse demand for electricity. As the cost of the technology goes down and its development advances, opportunities for applying it continue to grow. As a result, in the near future we will probably see an increase in massive application programs of this technology in rural Mexico. This document presents some of the most relevant problems needing immediate solution to increase the possibility of success in establishing photovoltaic programs for rural electrification in Mexico. These problems are influenced by economic, technological, engineering, infrastructural, social and political issues. The paper presents the requirements for an industrial center based on PV technology, the accessibility of replacement parts and maintenance service, development of financial solutions, training of the users and introducing norms and technical regulations. The document summarizes the issue of rural electrification in Mexico, emphasizing the low population density in these areas. A brief description is given of photovoltaics showing different configurations that can be utilized for supplying electricity to the rural areas. The issues involved and recommendations for introducing PV systems to the rural area are described.

14 SOLAR ENERGY↗

DOE OSTI2026

Metal–bicarbonate ion pairing in alkaline aqueous solutions from multilevel embedded correlated wavefunction theory and molecular dynamics

In this work, we examine ion-pairing mechanisms of bicarbonates in alkaline aqueous solutions with the divalent metal ions most abundantly present in seawater, namely, Ca 2+ and Mg 2+ . We employ a rare-event enhanced sampling approach within first-principles molecular dynamics to explore regions of phase space spanning solvent-shared to contact ion pairs. Second-order Møller–Plesset perturbation theory (MP2) corrections are subsequently applied in an embedding framework (EMB) to refine the electronic structure of stationary states and associated reaction barriers along the free-energy profiles while retaining the extended solvent effects at the density functional theory (DFT) level. Ca 2+ –HCO 3 − was previously hypothesized to exist in a solvent-shared ion pair (SSHIP) by DFT studies with an endergonic contact ion pair (CIP) formation; however, our EMB-MP2 refinement of the DFT ion-pairing pathways reveals that Ca 2+ and HCO 3 − form a virtually barrier-free CIP in alkaline solutions, with even more energetic ease than the widely studied Ca–CO 3 ion pair. We find qualitative agreement between DFT and EMB-MP2 for Mg 2+ —unlike Ca 2+ , Mg 2+ refuses to shed its strong hydration shell, thereby preferring a SSHIP state with a significant activation barrier to crossover to the CIP forms—a trait reminiscent of ion pairing in Mg–CO 3 and closely related to the kinetic limitations underlying the famous subject of the dolomite problem. Our study highlights the importance of improved electronic structure descriptions of liquids, modeled as a condensed phase of matter lacking in long-range crystalline order. It also strongly suggests that Ca 2+ –HCO 3 − CIPs are likely precursors involved in prenucleation of CaCO 3 mineral formation in seawater.

Sharma, Vidushi [Princeton Plasma Physics Laborato↗

DOE OSTI2026

Passive Confirmation of the Presence Of High-Explosive Material Via Neutron Transmission Spectroscopy

The goal of this project was to explore a novel approach for identifying presence and potential types of high explosives (HE) in treaty-controlled items with passive neutron sources by using passive neutron transmission spectroscopy. We predominantly look to measure relative elemental abundance of Carbon, Hydrogen, Oxygen, and Nitrogen (CHON) since most relevant materials are certain mix of these elements, as shown in Table 1. Previously, most studies exploring potential identification of CHON elemental content of targets in the vicinity of passive neutron source were focused solely on gamma spectroscopy using high resolution gamma detectors. Using neutron transmission spectroscopy with pulse-shape discrimination (PSD) capable organic scintillators is not only a novel idea in of itself, but arguably necessary for accurate identification of the CHON elemental content of an interrogated target. While high resolution gamma spectroscopy in principle can determine a presence of hydrogen and nitrogen, it is effectively blind to a quantitatively measuring areal densities of C, O. The initially proposed approach, described in detail in the project proposal, was to leverage previous Monte Carlo study on neutron transmission spectroscopy using slab shaped target interrogated with well-collimated (“pencil”) neutron beam. In this study, a simulated test CO 2 target was interrogated by a pencil neutron beam and transmitted neutron spectra were measured using PSD capable organic scintillator using MLEM based unfolding technique. To establish relative elemental content of carbon and oxygen, the unfolded neutron spectrum was fit with a parametrized combination of their respective elemental spectral templates. The individual spectral template was calculated by convoluting ENDF neutron crosssection with the known resolution of the PSD detector used in the study. This approach in the studied configuration successfully quantitatively established relative elemental composition of carbon and oxygen.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2026

Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code

Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.

Nikic, Dejan [University of New Mexico, Albuquerqu↗