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

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗
DOE OSTI2022

Multilayer Chitin–Chitosan–Cellulose Barrier Coatings on Poly(ethylene terephthalate)

Currently, effective barrier properties over a wide range of humidity are achieved by layering multiple plastics, resulting in structures that are nearly impossible to recycle through current systems. There is strong interest in developing barriers that are recyclable, biodegradable, or compostable. Chitin and cellulose nanomaterials are promising alternatives because of their high crystallinity, renewable sourcing, strong electrostatic interactions, and potential for biodegradation or composting. In this paper, we explore the potential to create multilayer materials with effective oxygen and water vapor barrier properties at elevated humidity by combining recyclable poly(ethylene terephthalate) (PET) films with single-layer and bilayer coatings of chitin nanowhiskers (ChNWs), chitosan (CS), and cellulose nanocrystals (CNCs). The addition of CS to the ChNW suspension improved the gas barrier properties of ChNW coatings, likely by filling voids, and short ChNWs (SChNW, 114 nm) performed better than long ChNWs (LChNW, 230 nm). At the optimal mass ratio of 1:1 SChNW/CS, the oxygen permeability (OP) of SChNW/CS-coated PET films was 5.1 cm 3 μm m –2 day –1 kPa –1 versus 15.9 cm 3 μm m –2 day –1 kPa –1 for neat PET. To take advantage of electrostatic attraction and hydrogen bonding, CNC was first coated on the PET, followed by coating with ChNW/CS, resulting in an OP as low as 3.6 cm 3 μm m –2 day –1 kPa –1 at 50% relative humidity (RH), and 6.4 cm 3 μm m –2 day –1 kPa –1 at 80% RH. Furthermore, after thermal treatment (120 °C, 2.5 h), the OP dropped to 2.2 cm 3 μm m –2 day –1 kPa –1 (50% RH) and 4.8 cm 3 μm m –2 day –1 kPa –1 (80% RH). The ChNW/CS and CNC coating can be efficiently removed with alkali treatment enabling possible recycling of the coated PET films.

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

Toward accurate measurement of electromagnetic field by retrieving and refining the center position of non-uniform diffraction disks in Lorentz 4D-STEM

Recent advancement in scanning transmission electron microscopy (STEM) allows the use of 4D-STEM, a technique that captures an electron diffraction pattern at each scan point in STEM, to measure electrostatic and magnetic potential and field in materials. However, accurate measurement, separation of the magnetic and electric signals, and removal of artifacts remain challenging, especially in the presence of complex non-uniform diffraction contrast within the disks. In this work, based on dynamic simulations of 4D-STEM patterns built upon superstructures consisting of millions of atoms to account for different sample thickness and edge geometries, we show how the shape and intensity distribution of the central disk are affected by multiple scattering. We propose a robust refinement procedure through iteration of the spin-sensitive peak position of the disk-center in the circular Hough transform filtered images from experimental Lorentz 4D-STEM dataset after minimizing the possible artifacts, such as those due to the change of thickness, dynamic scattering, and scanning process. We verify that caution must be taken as in practice the rigid-disk-shift model used to reconstruct induction maps can easily break down due to disk-protrusion when there exists a nonconstant phase gradient or thickness within the width of the probe. Through quantitative analysis and comparing experiment with calculation the effect of the non-spin-related intensity distribution inside the disk as well as that causes the disk shift due to the intensity-protrusion can be removed, and high-quality magnetic field mapping is possible.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Calculation of self-diffusion coefficients in supercritical carbon dioxide using mean force kinetic theory

This paper presents an application of mean force kinetic theory (MFT) to the calculation of the self-diffusivity of CO 2 in the supercritical fluid regime. Two modifications to the typical application of MFT are employed to allow its application to a system of molecular species. Therefore, the first is the assumption that the inter-particle potential of mean force can be obtained from the molecule center-of-mass pair correlation function, which in the case of CO 2 is the C–C pair correlation function. The second is a new definition of the Enskog factor that describes the effect of correlations at the surface of the collision volume. The new definition retains the physical picture that this quantity represents a local density increase, resulting from particle correlations, relative to that in the zero density homogeneous fluid limit. These calculations are facilitated by the calculation of pair correlation functions from molecular dynamics (MD) simulations using the FEPM2 molecular CO2 model. The self-diffusivity calculated from theory is in good agreement with that from MD simulations up to and slightly beyond the density at the location of the Frenkel line. The calculation is compared with and is found to perform similarly well to other commonly used models but has a greater potential for application to systems of mixed species and to systems of particles with long range interatomic potentials due to electrostatic interactions.

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

Non-linear Plasma Wake Growth of Electron Holes

An object’s wake in a plasma with small Debye length that drifts across the magnetic field is subject to electrostatic electron instabilities. Such situations include, for example, the moon in the solar wind wake and probes in magnetized laboratory plasmas. The instability drive mechanism can equivalently be considered drift down the potential-energy gradient or drift up the density gradient. The gradients arise because the plasma wake has a region of depressed density and electrostatic potential into which ions are attracted along the field. The non-linear consequences of the instability are analysed in this paper. At

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Anisotropic surface potentials induced by competitive ion adsorption enable the synthesis of branched cubic Pt mesocrystals

Creation of complex nanostructured materials through oriented attachment (OA) requires the manipulation of interparticle forces, including electrostatic repulsion, which depends strongly on surface potentials and can be modified through the effect of solution environment on interfacial chemistry. Here we show that time-dependent anisotropies in surface potential driven by competitive ion adsorption can alter facet-selectivity during OA. This phenomenon enables the synthesis of branched cubic Pt mesocrystals. Initially, Pt nanoparticles attach preferentially at their {100} facets to form a well-defined cubic core. Over time, changes in ion adsorption shift the attachment preference to the {111} facets, promoting branch formation. In both stages, anisotropic surface potentials generate electrostatic torques that align the particles prior to attachment. These findings demonstrate a generalizable strategy for directing the architecture of nanomaterials through time-resolved control of interfacial chemistry during OA, offering new pathways for the design of complex mesoscale structures.

Bae, Yuna [Pacific Northwest National Laboratory (↗
DOE OSTI2025

Quantum confining excitons with an electrostatic moiré superlattice

Quantum confining excitons has been a persistent challenge in the pursuit of strong exciton interactions and quantum light generation. Unlike electrons, which can be readily controlled via electric fields, imposing strong nanoscale potentials on excitons to enable quantum confinement has proven challenging. In this study, we utilize piezoelectric force microscopy to image the domain structures of twisted hexagonal boron nitride (h – BN), revealing evidence of strong in-plane electric fields at the domain boundaries. By placing a monolayer MoSe 2 only one to two nanometers away from the twisted h – BN interface, we observe energy splitting of neutral excitons and Fermi polarons by several millielectronvolts at the moiré domain boundaries. By directly correlating local structural and optical properties, we attribute such observations to excitons confined in a nanoscale one-dimensional electrostatic potential created by the strong in-plane electric fields at the moiré domain boundaries. Intriguingly, this 1D quantum confinement results in pronounced polarization anisotropy in the excitons’ reflection and emission, persistent to temperatures as high as ~80 Kelvins. Furthermore, these findings open new avenues for exploring and controlling strongly interacting excitons for classical and quantum optoelectronics.

2-dimensional systems↗
DOE OSTI2020

Optimized utilization of COMB3 reactive potentials in LAMMPS

An investigation to optimize the application of the third-generation charge optimized many-body (COMB3) interatomic potential and associated input parameters was carried out through the study of solid–liquid interactions in classical molecular dynamics simulations. The rates of these molecular interactions are understood through the wetting rates of water nano-droplets on a bare copper (111) surface. Implementing the Langevin thermostat, the influence of simulation time step, the number of atoms in the system, the frequency at which charge equilibration is performed, and the temperature relaxation rate are all examined. The results indicate that time steps of 0.4 fs are possible when using longer relaxation times for the system temperature, which is almost double the typical time step used for reactive potentials. Additionally, the use of the charge equilibration allows for a fewer atomic layers to be used in the Cu slab. In addition, charge equilibrium schemes do not need to be performed every time step to ensure accurate charge transfer. Interestingly, the rate of wetting for the nanodroplets is dominantly dependent on the temperature relaxation time, which is predicted to significantly change the viscosity of the water droplets. This work provides a pathway for optimizing simulations using the COMB3 reactive interatomic potential.

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

Shadow molecular dynamics for flexible multipole models

Shadow molecular dynamics provide an efficient and stable atomistic simulation framework for flexible charge models with long-range electrostatic interactions. Shadow molecular dynamics simulations are driven by approximate “shadow” Born–Oppenheimer potentials for which the exact charges and forces are directly accessible without relying on costly (and approximate) iterative solvers. While previous implementations have been limited to atomic monopole charge distributions, we extend this approach to flexible multipole models. We derive detailed expressions for the shadow energy functions, potentials, and force terms, explicitly incorporating monopole–monopole, dipole–monopole, and dipole–dipole interactions. In our formulation, both atomic monopoles and atomic dipoles are treated as extended dynamical variables alongside the propagation of the nuclear degrees of freedom. We demonstrate that introducing the additional dipole degrees of freedom preserves the stability and accuracy previously seen in monopole-only shadow molecular dynamics simulations. In addition, we present a shadow molecular dynamics scheme where the monopole charges are held fixed while the dipoles remain flexible. Our extended shadow dynamics provide a framework for stable, computationally efficient, and versatile molecular dynamics simulations involving long-range interactions between flexible multipoles. This is of particular current interest in combination with machine-learned interatomic potentials, including long-range electrostatic interactions.

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

The emergence of three-dimensional chiral domain walls in polar vortices

Abstract Chirality or handedness of a material can be used as an order parameter to uncover the emergent electronic properties for quantum information science. Conventionally, chirality is found in naturally occurring biomolecules and magnetic materials. Chirality can be engineered in a topological polar vortex ferroelectric/dielectric system via atomic-scale symmetry-breaking operations. We use four-dimensional scanning transmission electron microscopy (4D-STEM) to map out the topology-driven three-dimensional domain walls, where the handedness of two neighbor topological domains change or remain the same. The nature of the domain walls is governed by the interplay of the local perpendicular (lateral) and parallel (axial) polarization with respect to the tubular vortex structures. Unique symmetry-breaking operations and the finite nature of domain walls result in a triple point formation at the junction of chiral and achiral domain walls. The unconventional nature of the domain walls with triple point pairs may result in unique electrostatic and magnetic properties potentially useful for quantum sensing applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tokamak ITG-KBM transition benchmarking with the mixed variables/pullback transformation electromagnetic gyrokinetic scheme

Electromagnetic gyrokinetic simulation of high temperature plasma is required to predict confinement in magnetic fusion devices and has posed challenges for existing codes. In this paper, we demonstrate successful global gyrokinetic simulation of the ion temperature gradient-driven mode-kinetic ballooning mode transition in a toroidal fusion plasma test case using the mixed variables/pullback transformation (MV/PT) scheme with the particle-in-cell codes XGC and ORB5, and compare to results from a conventional continuum code from the literature. Furthermore, the MV/PT scheme combines explicit time integration with mitigation of the well-known electromagnetic gyrokinetic “cancelation problem.” We calculate eigenmodes in the electrostatic and parallel vector potentials, and find good agreement in growth rate, real frequency, and the normalized plasma pressure of mode transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Partial charge transfer and absence of induced magnetization in EuS(111)/Bi 2 Se 3 heterostructures

We report that heterostructures made from topological and magnetic insulators promise to form excellent platforms for new electronic and spintronic functionalities mediated by interfacial effects. We report the results of a first-principles density functional theory study of the geometric, electronic structure, and magnetic properties of the EuS( 111 )/Bi 2 Se 3 interface, including van der Waals and relativistic spin-orbit effects. In contrast to previous theoretical studies, we find no appreciable magnetic anisotropy in such a heterostructure. We also do not see additional induced magnetization at the interface or the magnetic proximity effect on the topological states. This is due to the localized nature of Eu moments and because of a partial charge transfer of ~ 0.5 electron from Eu to Se. The formation of the surface dipole shifts the Dirac cone about 0.4 eV below the chemical potential, and the associated electrostatic screening moves the topological state from the first to the second quintuple layer of Bi 2 Se 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Surface Functionalized Barium Titanate Nanoparticles: A Combined Experimental and Computational Study

Barium titanate (BTO) nanoparticles show great potential for use in electrostatic capacitors with high energy density. This includes both polymer composite and sintered capacitors. However, questions about the nanoparticles' size distribution, amount of agglomeration, and surface ligand effect on performance properties remain. Reducing particle agglomeration is a crucial step to understanding the properties of nanoscale particles, as agglomeration has significant effects on the composite dielectric constant. BTO surface functionalization using phosphonic acids is known reduce BTO nanoparticle agglomeration. We explore solution synthesized 10 nm BTO particles with tert-butylphosphonic acid ligands. Recent methods to quantifying agglomeration using an epoxy matrix before imaging shows that tert-butylphosphonic acid ligands reduce BTO agglomeration by 33%. Thermometric, spectroscopic, and computational methods provide confirmation of ligand binding and provide evidence of multiple ligand binding modes on the BTO particle surface.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-Fidelity Arc-Discharge Model for Hydrogen-Plasma-Smelting-Reduction of Iron Ore

Electrification and use of renewable hydrogen is currently a necessity for decarbonizing the iron-and-steel industry. In this regard, hydrogen plasma smelting reduction (HPSR) is a novel pathway that is being explored for reduction of iron ore. HPSR provides several decarbonization merits compared to conventional blast furnaces. Firstly, the use of renewable hydrogen drastically reduces the CO2 emissions compared to the use of coke. Secondly, renewable electricity in the form of a thermal plasma for making reactive hydrogen species (radicals, ions) are more efficient at reducing iron ore compared to neutral H2. Thirdly, a molten product compatible with downstream processes is obtained from the intense heat transfer from the plasma. However, the scale-up of this technology requires fundamental exploration of hydrogen plasma dynamics and its interaction with complex solid material that include phase changing iron-ore and slag. In this work, we present a first principles continuum scale model for thermal plasmas in Ar/H2 gas mixtures typically used for HPSR. The thermal plasma governing equations for mass, momentum and energy with Lorentz force and Joule heating source terms are solved along with electromagnetic equations for electrostatic and magnetic vector potential. Our solver will be based on Pele, a suite of reacting flow solvers designed for advanced scientific computing architectures (Henry De Frahan et al., Proceedings of SIAM Parallel Processing, 13-25, 2024), and will utilize adaptive mesh generation for enhanced resolutions at locations of intense physicochemical interactions. This study will present the impact of Ar to H2 ratios on excited/dissociated hydrogen species concentrations, plasma temperature and conductivity along with the impact of outgassed species (water, metal vapor, O, OH radicals) from ore surface on gas phase chemistry. Furthermore, the heat and species flux to the surface will be quantified as a function of applied voltages in a transferred arc configuration.

hydrogen plasma↗
DOE OSTI2023

Probing the self-ionization of liquid water with ab initio deep potential molecular dynamics

The chemical equilibrium between self-ionized and molecular water dictates the acid–base chemistry in aqueous solutions, yet understanding the microscopic mechanisms of water self-ionization remains experimentally and computationally challenging. Herein, Density Functional Theory (DFT)–based deep neural network (DNN) potentials are combined with enhanced sampling techniques and a global acid–base collective variable to perform extensive atomistic simulations of water self-ionization for model systems of increasing size. The explicit inclusion of long-range electrostatic interactions in the DNN potential is found to be crucial to accurately reproduce the DFT free energy profile of solvated water ion pairs in small (64 and 128 H 2 O) cells. The reversible work to separate the hydroxide and hydronium to a distance S is found to converge for simulation cells containing more than 500 H 2 O, and a distance of ∼ 8 Å is the threshold beyond which the work to further separate the two ions becomes approximately zero. The slow convergence of the potential of mean force with system size is related to a restructuring of water and an increase of the local order around the water ions. Calculation of the dissociation equilibrium constant illustrates the key role of long-range electrostatics and entropic effects in the water autoionization process.

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

Potential-induced wetting and dewetting in hydrophobic nanochannels for mass transport control

Wetting and dewetting transitions play a central role in controlling the hydrophobicity of the lining of biological channels in order to regulate aqueous solution permeation. Understanding of the operational characteristics of biological nanochannels led to experimental efforts to mimic their behavior and to achieve potential-induced, repeatedly-switchable wettability transitions in synthetic nanochannels in the early 2010s. Since then, research has identified conditions needed to produce reversible wettability transitions using a number of different environmental stimuli - such as light, pH, and electrostatic charge - in addition to potential. Furthermore, nascent understanding of the underlying phenomena in synthetic nanochannels was rapidly followed by practical applications, including oil-water separations, drug release, and electroactive flow control based on switchable wettability. Furthermore, more practical applications are being developed continuously, as the physical and chemical principles that govern hydrophobic gating at the nanoscale are further elucidated, making it possible to exploit wettability as a design element in nanofluidic systems.

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

Extending Badger's rule. I. The relationship between energy and structure in hydrogen bonds

Here, we derive a new expression for the strength of a hydrogen bond (V HB ) in terms of the elongation of the covalent bond of the donor fragment participating in the hydrogen bond (Δr HB ) and the intermolecular coordinates R (separation between the heavy atoms) and θ (deviation of the hydrogen bond from linearity). The expression includes components describing the covalent D–H bond of the hydrogen bond donor via a Morse potential, the Pauli repulsion, and electrostatic interactions between the constituent fragments using a linear expansion of their dipole moment and a quadratic expansion of their polarizability tensor. We fitted the parameters of the model using ab initio electronic structure results for six hydrogen bonded dimers, namely, NH 3 –NH 3 , H 2 O–H 2 O, HF–HF, H 2 O–NH 3 , HF–H 2 O, and HF–NH 3 , and validated its performance for extended parts of their potential energy surfaces, resulting in a mean absolute error ranging from 0.07 to 0.31 kcal/mol. The derived expression describes the energy–structure relationship in terms of a single structural parameter, namely, the elongation of the donor’s covalent bond (Δr HB ), and suggests the novel relationship of 8.0 kcal/mol pm −1 (or 0.8 kcal/mol per 0.001 Å elongation). This structural parameter is easily obtained from theory and can serve as the single descriptor of the strength of individual hydrogen bonds.

Santis, Garrett D. [Univ. of Washington, Seattle, ↗
DOE OSTI2024

X-ray photoelectron spectroscopy of epitaxial films and heterostructures

X-ray photoelectron spectroscopy is a powerful experimental technique that yields invaluable information on a range of phenomena that occur in solids, liquids, and gasses. The binding energy and shape of a photoemission peak is sensitive not only to the atomic number, valence and orbital from which the electron is ejected, but also to complex many-body effects that accompany photoemission. Provided the influences of these different drivers of spectral line shapes can be disentangled, a great deal can be learned about the electronic structures of specific atoms in the material of interest. In addition to these largely local effects, the long-range electrostatic environment and resulting electric potential at the emitting atom have a direct effect on the measured binding energies. Furthermore, this fact opens the door to extracting information about the dependence of the valence and conduction band minima on depth below the surface, which in turn allows both vertical and lateral electrical transport data to be better understood. One purpose of this Report is to describe how the different physical forces described above impact spectral properties of complex oxide epitaxial films. This class of materials typically incorporates transition metal cations in different valences and of all the elements, these exhibit the most complex core-level spectra. A second purpose is to show how a comprehensive understanding of local physical effects in x-ray photoemission allows one to extract detailed information on internal electric fields and band edge discontinuities in heterostructures involving complex oxides from core-level line shapes.

47 OTHER INSTRUMENTATION↗