Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Electrostatic Potential”

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

Triboelectric backgrounds to radio-based polar ultra-high energy neutrino (UHEN) experiments

Here, in the hopes of observing the highest-energy neutrinos (E EeV) populating the Universe, both past (RICE, AURA, ANITA) and current (RNO-G, ARIANNA, ARA and TAROGE-M) polar-sited experiments exploit the impulsive radio emission produced by neutrino interactions. In such experiments, rare single event candidates must be unambiguously identified above backgrounds. Background rejection strategies to date primarily target thermal noise fluctuations and also impulsive radio-frequency signals of anthropogenic origin. In this paper, we consider the possibility that ‘fake’ neutrino signals may also be generated naturally via the ‘triboelectric effect.’ This broadly describes any process in which force applied at a boundary layer results in displacement of surface charge, leading to the production of an electrostatic potential difference V. Wind blowing over granular surfaces such as snow can induce such a potential difference, with subsequent coronal discharge. Discharges over timescales as short as nanoseconds can then lead to radio-frequency emissions at characteristic MHz–GHz frequencies. Using data from various past (RICE, AURA, SATRA, ANITA) and current (RNO-G, ARIANNA and ARA) neutrino experiments, we find evidence for such backgrounds, which are generally characterized by: (a) a threshold wind velocity which likely depends on the experimental trigger criteria and layout; for the experiments considered herein, this value is typically $\mathscr{O}$(10 m/s), (b) frequency spectra generally shifted to the low-end of the frequency regime to which current radio experiments are typically sensitive (100–200 MHz), (c) for the strongest background signals, an apparent preference for discharges from above-surface structures, although the presence of more isotropic, lower amplitude triboelectric discharges cannot be excluded.

Background rejection↗

Quantifying the anion effect of gas solubility within ionic liquids using the solvation affinity index

The solubility of gas molecules within ionic liquids (ILs) is strongly dictated by the anion effect and the free volume effect, and of these two influences, the anion effect can be challenging to predict. We evaluate a wide range of different gas solutes (O 2 , NH 3 , CH 4 , H 2 S, SO 2 , Ar, H 2 , N 2 , CO 2 ) in different ILs, and we find that the solubility is strongly correlated to one of three different IL properties, which can be directly calculated from the anion electrostatic potential (ESP) surfaces. Thus, the solutes can be grouped into three fundamental types, leading to a general descriptor, the solvation affinity index (SAI), which can identify the affinity between a wide range of gas solutes and IL solvents. The SAI of nine common industrial gases with 47 common anions are generated. Furthermore, this can provide a fast approach for pre-screening IL solvents or even for ILs embedded within porous supports.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An immersed interface method for microstructure-scale electrochemical battery models: numerical formulation and performance portable implementation

We present the numerical formulation, verification, and performance portable implementation of an immersed interface method for microstructure scale electrochemical modeling of batteries. The innovation in this approach is the resolution of chemical species and electrostatic potential discontinuities at active interfaces without the use of interface conforming unstructured grids. A unified formulation on Cartesian grids for all domains (electrodes and electrolyte) is used with interfacial flux conditions applied using volume fraction or “color” function gradients. We have developed one dimensional and two dimensional test cases with analytic solutions for electrochemical modeling using which we verified the consistency and accuracy of our scheme. Our solver is also validated against solutions from a macroscale model and an unstructured multi-subdomain solver for a full lithium ion cell. We then demonstrated the utility of our solver on an image-based complex battery electrode microstructure at high charging rate. Our technique also exhibits good scalability on distributed memory architectures using central processing units (CPU), with problem sizes up to 1.8 billion degrees of freedom and with 5400 ranks. Initial performance studies of our open-source performance portable solver showed about 70 times speed up using a graphics processing unit (GPU) compared to single compute core for a problem with 4 million cells.

25 ENERGY STORAGE↗

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An ion specific continuum model on ionic surfactant's binary phase diagram in aqueous solution

Here, a continuum aggregation model is proposed to account for specific ion effects, enabling accurate phase diagrams calculations for the micellar, cylindrical and lamellar aggregates of sodium/potassium carboxylate surfactants in aqueous solution across a range of temperatures. Three groups of concentrations at distinctive temperatures are fitted to build empirical temperature dependence given the limited availability of relevant experimental measurements. The specific ion effects are manifested in the aggregates' surface tension as well as in the distributions of counter-ions' concentrations in the vicinity of the aggregates. The aggregates' geometric sizes are well-reproduced. The differential evolution algorithm is applied to address boundary conditions of the electrostatic potential, aggregate size optimization as well as the equilibrium of monomers transferring between the aggregate and aqueous region.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fabrication and characterization of boron-terminated tetravacancies in monolayer hBN using STEM, EELS and electron ptychography

Tetravacancies in monolayer hexagonal boron nitride (hBN) with consistent edge termination (boron or nitrogen) form triangular nanopores with electrostatic potentials that can be leveraged for applications such as selective ion transport and neuromorphic computing. In order to quantitatively predict the properties of these structures, an atomic-level understanding of their local electronic and chemical environments is required. Moreover, robust methods for their precision manufacture are needed. Here we use electron irradiation in a scanning transmission electron microscope (STEM) at a high dose rate to drive the formation of boron-terminated tetravacancies in monolayer hBN. Characterization of the defects is achieved using aberration-corrected STEM, monochromated electron energy-loss spectroscopy (EELS), and electron ptychography. Z-contrast in STEM and chemical fingerprinting by core-loss EELS enable identification of the edge terminations, while electron ptychography gives insight into structural relaxation of the tetravacancies and provides evidence of enhanced electron density around the defect perimeters indicative of bonding effects.

2D hBN↗

Zonally dominated dynamics and Dimits threshold in curvature-driven ITG turbulence

The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$ -pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a finite-amplitude saturated state dominated by strong zonal-flow and zonal temperature perturbations, and a blow-up state that fails to saturate on a box-independent scale. We argue that this transition is equivalent to the Dimits transition from a low-transport to a high-transport state seen in gyrokinetic numerical simulations (Dimits et al. , Phys. Plasmas , vol. 7, 2000, 969). A quasi-static staircase-like structure of the temperature gradient intertwined with zonal flows, which have patch-wise constant shear, emerges near the Dimits threshold. The turbulent heat flux in the low-collisionality near-marginal state is dominated by turbulent bursts, triggered by coherent long-lived structures closely resembling those found in gyrokinetic simulations with imposed equilibrium flow shear (van Wyk et al. , J. Plasma Phys. , vol. 82, 2016, 905820609). The breakup of the low-transport Dimits regime is linked to a competition between the two different sources of poloidal momentum in the system – the Reynolds stress and the advection of the diamagnetic flow by the $\boldsymbol {E}\times \boldsymbol {B}$ flow. By analysing the linear ion-temperature-gradient modes, we obtain a semi-analytic model for the Dimits threshold at large collisionality.

Physics↗

Bispectral analysis of broadband turbulence and geodesic acoustic modes in the T-10 tokamak

Local fluctuations of electrostatic potential, poloidal electric field, magnetic potential and electron density are simultaneously measured in the T-10 tokamak by a heavy ion beam probe (HIBP) having a five-slit energy analyser, which allows an estimate of the turbulent particle flux and $\boldsymbol {E}\times \boldsymbol {B}$ rotation velocity in the off-minor-axis gradient zone of the toroidal plasma column. The high spatial and temporal resolution of the modern multichannel HIBP makes it an effective tool to study plasma oscillations. Motivated by previous work that has documented time-resolved interactions between measured plasma parameters using correlation analysis (coherence of $E_{\textrm {pol}}$ and density $n_e$ , and cross-phase), a new result from bicorrelation analysis (bicoherence of magnetic potential $A_\zeta$ and density $n_e$ , and biphase) is reported for documenting the evidence of wave–wave coupling and energy transfer associated with the interaction between geodesic acoustic modes (GAM) and broadband, quasi-coherent modes.

Physics↗

Chaotic and integrable magnetic fields in one-dimensional hybrid Vlasov–Maxwell equilibria

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗

Drift-cyclotron loss-cone instability in 3-D simulations of a sloshing-ion simple mirror

The kinetic stability of collisionless, sloshing beam-ion (45° pitch angle) plasma is studied in a three-dimensional (3-D) simple magnetic mirror, mimicking the Wisconsin high-temperature superconductor axisymmetric mirror experiment. The collisional Fokker–Planck code CQL3D-m provides a slowing-down beam-ion distribution to initialize the kinetic-ion/fluid-electron code Hybrid-VPIC, which then simulates free plasma decay without external heating or fuelling. Over 1 – 10 μs, drift-cyclotron loss-cone (DCLC) modes grow and saturate in amplitude. The DCLC scatters ions to a marginally stable distribution with gas-dynamic rather than classical-mirror confinement. Sloshing ions can trap cool (low-energy) ions in an electrostatic potential well to stabilize DCLC, but DCLC itself does not scatter sloshing beam-ions into the said well. Instead, cool ions must come from external sources such as charge-exchange collisions with a low-density neutral population. Manually adding cool ∼1keV ions improves beam-ion confinement several-fold in Hybrid-VPIC simulations, which qualitatively corroborates prior measurements from real mirror devices with sloshing ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

From Cation Order to Disorder: Unlocking Ion Transport Pathways in Li–Zn–Zr–Cl Halospinels

Lithium metal chloride halospinels of the general formula Li 2 MCl 4 are a promising class of earth-abundant ion conductors for all-solid-state batteries. However, poor roomtemperature ionic conductivity has historically limited their use in practical applications. Here, we substitute Zr 4+ into Li 2 ZnCl 4 along the series Li 2−2x/3 Zn 1−x Zr 2x/3 Cl 4 (x = 0, 0.1, 0.3, 0.4, 0.6, 0.9, and 1.0) to understand how cation disorder and vacancy tuning impact ion transport in “normal” halospinels. Aliovalent Zr 4+ substitution increases ionic conductivity by nearly 5 orders of magnitude, from 1.320(3) × 10 −9 S cm −1 in Li 2 ZnCl 4 to 6.74(1)× 10 −5 S cm −1 for x = 0.6. Average and local structure characterization through synchrotron X-ray diffraction (SXRD) and neutron pair distribution function (nPDF) analysis reveal that Zr 4+ redistributes the Zn 2+ and Li + sublattices into previously unoccupied interstitial sites, which form new low-energy hopping pathways that facilitate ion transport. We rationalize the dramatic rearrangement of the cation local structure by considering the coordination preferences of the cations and the potential electrostatic penalties incurred by the higher-valent Zr 4+ cations. This work delivers an atomistic understanding of substitution-induced cation disorder and ion transport properties in a new family of earth-abundant halospinels.

Cardoza, Abby M. [Colorado School of Mines, Golden↗

Impact of Noncovalent Interactions on the Structural Chemistry of Thorium(IV)-Aquo-Chloro Complexes

Five novel tetravalent thorium (Th) compounds that consist of Th(H 2 O) x Cl y structural units were isolated from acidic aqueous solutions using a series of nitrogen-containing heterocyclic hydrogen (H) bond donors. Taken together with three previously reported phases, the compounds provide a series of monomeric Th IV complexes wherein the effects of noncovalent interactions (and H-bond donor identity) on Th structural chemistry can be examined. Seven distinct structural units of the general formulas [Th(H 2 O) x Cl 8–x ] x-4 (x = 2, 4) and [Th(H 2 O) x Cl 9–x] x-5 (x = 5–7) are described. The complexes range from chloride-deficient [Th(H 2 O) 7 Cl 2 ] 2+ to chloride-rich [Th(H 2 O) 2 Cl 6 ] 2– species, and theory was used to understand the relative energies that separate complexes within this series via the stepwise chloride addition to an aquated Th cation. Electronic structure theory predicted the reaction energies of chloride addition and release of water through a series of transformations, generally highlighting an energetic driving force for chloride complexation. To probe the role of the counterion in the stabilization of these complexes, electrostatic potential (ESP) surfaces were calculated. Finally, the ESP surfaces indicated a dependence of the chloride distribution about the Th metal center on the pK a of the countercation, highlighting the directing effects of noncovalent interactions (e.g., Hbonding) on Th speciation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

InteraChem: Virtual Reality Visualizer for Reactive Interactive Molecular Dynamics

Interactive molecular dynamics in virtual reality (IMD-VR) simulations provide a digital molecular playground for students as an alternative or complement to traditional molecular modeling kits or 2D illustrations. Previous IMD-VR studies have used molecular mechanics to enable simulations of macromolecules such as proteins and nanostructures for the classroom setting with considerable success. Here, we present the InteraChem molecular visualizer, intended for reactive IMD-VR simulation using semiempirical and ab initio methods. InteraChem visualizes not only the molecular geometry but also (1) isosurfaces such as molecular orbitals and electrostatic potentials and (2) two-dimensional graphs of time-varying simulation quantities such as kinetic/potential energy, internal coordinates, and user-applied force. Additionally, InteraChem employs speech recognition to facilitate user interaction and introduces a novel “atom happiness” visualization using emojis to indicate the energetic feasibility of a particular bonding arrangement. Furthermore, we include a set of accompanying exercises that we have used to teach chemical reactivity in small molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of Intermolecular Interactions in Deep Eutectic Solvents for CO 2 Capture: Vibrational Spectroscopy and Quantum Chemical Studies

Recent research and reviews on CO 2 capture methods, along with advancements in industry, have highlighted high costs and energy-intensive nature as the primary limitations of conventional direct air capture and storage (DACS) methods. In response to these challenges, deep eutectic solvents (DESs) have emerged as promising absorbents due to their scalability, selectivity, and lower environmental impact compared to other absorbents. However, the molecular origins of their enhanced thermal stability and selectivity for DAC applications have not been explored before. Therefore, the current study focuses on a comprehensive investigation into the molecular interactions within an alkaline DES composed of potassium hydroxide (KOH) and ethylene glycol (EG). Combining Fourier transform infrared (FT-IR) and quantum chemical calculations, the study reports structural changes and intermolecular interactions induced in EG upon addition of KOH and its implications on CO 2 capture. Experimental and computational spectroscopic studies confirm the presence of noncovalent interactions (hydrogen bonds) within both EG and the KOH-EG system and point to the aggregation of ions at higher KOH concentrations. Additionally, molecular electrostatic potential (MESP) surface analysis, natural bond orbital (NBO) analysis, quantum theory of atoms-in-molecules (QTAIM) analysis, and reduced density gradient-noncovalent interaction (RDG-NCI) plot analysis elucidate changes in polarizability, charge distribution, hydrogen bond types, noncovalent interactions, and interaction strengths, respectively. Evaluation of explicit and hybrid models assesses their effectiveness in representing intermolecular interactions. This research enhances our understanding of molecular interactions in the KOH-EG system, which are essential for both the absorption and desorption of CO 2 . The study also aids in predicting and selecting DES components, optimizing their ratios with salts, and fine-tuning the properties of similar solvents and salts for enhanced CO 2 capture efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

STEPs-SOL, a Peptoid Force Field Parameterization to Include Solvent Effects

As peptoids (N-substituted glycines) continue to gain popularity as a class of biomimetic polymers, the importance and demand for accurate force fields in molecular simulations also grow. Building on the vacuum-optimized Systematic and Extensible Force Field for Peptoids (STEPs) force field, here we present STEPs-SOL, a novel peptoid force field parametrization that effectively incorporates solvent effects to enhance the accuracy of peptoid simulations. The development of STEPs-SOL is based on the need for precise electrostatic modeling achieved through solvent-specific partial charge optimization. Here, our systematic approach significantly improves agreement with experimental measurements, reducing the mean absolute error in cis/trans ratio predictions (ΔG c/t ) by an average of 38% across multiple peptoid residues and solvent environments. This improved parametrization addresses computational challenges associated with nonbonded energies while maintaining a workflow that relies on high-level quantum mechanical data rather than depending solely on limited experimental equilibrium properties. By evaluating the effects of conformational bias in restrained electrostatic potential (RESP) charge generation and examining their impact on peptoid conformations in various solvents, we enhance our understanding of peptoid structural dynamics while providing a more accurate modeling framework.

force field↗

Electronic Structure Modification of Rectangular Phosphorene Quantum Dots Via Edge Passivation

Phosphorene is an elemental two-dimensional material that shows great promise as a metal-free photocatalyst owing to its high carrier mobility and direct band gap. Doping and edge functionalization of phosphorene can be exploited to tune its band gap and alter other properties. Herein, the electronic properties of a series of edge-passivated monolayer phosphorene quantum dots are investigated using an array of functional groups (OH, SH, NH 2 , SCH 3 , OCN, CN, and Cl) that are known to have electron-donating or electron-withdrawing capabilities. By employing density functional theory calculations, we demonstrate that edge functionalization can perturbatively shift the HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) energy levels while simultaneously preserving the geometric structure and other favorable properties of pristine phosphorene. We use the average overall electrostatic potential of these systems to successfully predict the relative positioning of the HOMO and LUMO energy levels while presenting qualitative explanations on potential jumps introduced by the dipole layers formed due to edge functionalization. Furthermore, our results suggest that selective functionalization permits considerable control over the band edge positions, which could be utilized for applications in energy and optoelectronic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Ab Initio -Derived Force Field for Amorphous Silica Interfaces for Use in Molecular Dynamics Simulations

Here, we present a classical interatomic force field, silica-DDEC, to describe the interactions of amorphous and crystalline silica surfaces, parametrized using density functional theory-based charges. Charge schemes for silica surfaces were developed using the density-derived electrostatic and chemical (DDEC) method, which reproduces atomic charges of the periodic models as well as the electrostatic potential away from the atom sites. Lennard–Jones parameters were determined by requiring the correct description of (i) the amorphous silica density, coordination defects, and local coordination geometry, relative to experimental measurements, and (ii) water-silica interatomic distances compared with ab initio results. Deprotonated surface silanol sites are also described within the model based on DDEC charges. The result is a general electronic structure-derived model for describing fully flexible amorphous and crystalline silica surfaces and interactions of liquids with silica surfaces of varying structure and protonation state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Voltage-Dependent First-Principles Barriers to Li Transport within Li-Ion Battery Solid Electrolyte Interphases

Charging a Li-ion battery requires Li-ion transport between the cathode and the anode. This Li-ion transport is dependent on (among other factors) the electrostatic environment that the ion encounters within the solid electrolyte interphase (SEI), which separates the anode from the surrounding electrolyte. A previous first-principles work has illuminated the reaction barriers through likely atomistic SEI environments but has had difficulty accurately reflecting the larger electrostatic potential landscape that an ion encounters moving through the SEI. In this work, we apply the recently developed quantum continuum approximation (QCA) technique to provide an equilibrium electronic potentiostat for first-principles interface calculations. Using QCA, we calculate the potential barrier for Li-ion transport through LiF, Li 2 O, and Li 2 CO 3 SEIs along with LiF–LiF and LiF–Li 2 O grain boundaries, all paired with Li metal anodes. Here, we demonstrate that the SEI potential barrier is dependent on the electrochemical potentials of the anode in each system. Finally, we use these techniques to estimate the change in the diffusion barrier for a Li ion moving in a LiF SEI as a function of the anode potential. We find that properly accounting for interface and electronic voltage effects significantly lowers reaction barriers compared with previous literature results.

25 ENERGY STORAGE↗