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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

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

Efficient and Photostable Organic Solar Cells Achieved by Alloyed Dimer Acceptors with Tailored Linker Structures

Abstract High power conversion efficiency (PCE) and long‐term stability are prerequisites for commercialization of organic solar cells (OSCs). Herein, two dimer acceptors (DYTVT and DYTCVT) are developed with different properties through linker engineering, and study their effects as alloy‐like acceptors on the photovoltaic performance and photostability of OSCs. These ternary OSCs effectively combine the advantages of both dimer acceptors. DYTVT, characterized by its high backbone planarity, ensures elevated electron mobility and high glass‐transition temperature (T g ), leading to efficient charge transport and enhanced photostability of OSCs. Conversely, DYTCVT, with its significant dipole moment and electrostatic potential, enhances compatibility of the alloy acceptors with donors and refines the blend morphology, facilitating efficient charge generation in OSCs. Consequently, D18:DYTVT:DYTCVT OSCs exhibit higher PCE (18.4%) compared to D18:MYT (monomer acceptor, PCE = 16.5%), D18:DYTVT (PCE = 17.4%), and D18:DYTCVT (PCE = 17.0%) OSCs. Furthermore, owing to higherT g of alloy acceptors (133 °C) than MYT (T g = 80 °C) and DYTCVT (T g = 120 °C), D18:DYTVT:DYTCVT OSCs have significantly higher photostability (t 80% lifetime = 4250 h under 1‐sun illumination) compared to D18:MYT (t 80% lifetime = 40 h) and D18:DYTCVT OSCs (t 80% lifetime = 2910 h).

Chemistry↗
DOE OSTI2025

Mechanistic Insights into the Demethylation of Lignin‐Derived Structures Using Protic Ionic Liquids: A Density Functional Theory Study

Lignin valorization is restricted by the stability of its methoxy groups, creating a critical need for efficient demethylation strategies. Here, in this study, density functional theory (DFT) is employed to dissect the mechanistic pathways of demethylation in lignin model compounds, guaiacol and syringol, using protic ionic liquids (PILs) that act as both solvent and catalyst. Conductor‐like Screening Model for Real Solvents (COSMO‐RS) analysis identifies monoethanolammonium acetate ([MEOA][Ace]) as the most promising medium, attributed to its strong hydrogen bonding network and solvation ability. By integrating implicit and explicit solvation models, it is revealed that an acid‐catalyzed hydrolytic mechanism governs demethylation, with PILs stabilizing crucial transition states and intermediates. Complementary electronic structure evaluations, including highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO‐LUMO) gap analysis, charge distribution, and electrostatic potential mapping, demonstrate how PILs lower energetic barriers and enhance reactivity. To simulate realistic environments, this study is extended to lignin dimer complexes with varying water content, uncovering how water‐bridged solvation reverses demethylation preference from guaiacyl to syringyl units. This mechanistic shift aligns with experimental observations showing faster S‐unit reactivity in hydrated systems. Together, these findings provide atomic‐level insight into lignin demethylation dynamics and highlight how tuning acid concentration in PILs can accelerate kinetics, enabling a rational pathway toward next‐generation biomass conversion technologies.

COSMO-RS↗
DOE OSTI2026

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↗
DOE OSTI2024

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↗
DOE OSTI2025

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↗
DOE OSTI2025

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↗
DOE OSTI2026

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↗
DOE OSTI2024

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↗
DOE OSTI2025

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↗
DOE OSTI2024

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↗
DOE OSTI2025

Polarization-Controlled Structural Modulation in the Single Atomic Layer at the PbZr 0.2 Ti 0.8 O 3 /LaNiO 3 Interface

Conductivity modulation via ferroelectric polarization coupling with LaNiO 3 (LNO) is demonstrated at an epitaxial ferroelectric–LNO interface. Conductivity measurements, varying the thickness of the LNO channel, show that this phenomenon is confined to a few atomic layers at the interface. Combining in situ biasing and off-axis holography, we mapped out electrostatic potentials at the PbZr 0.2 Ti 0.8 O 3 (PZT)/LNO/SrTiO 3 (STO) heterostructure upon polarization switching. Using aberration-corrected STEM, the interfacial atomic structures were investigated for the two different PZT polarization states. Polarization in PZT induces a significant change in the in-plane O–Ni–O bond angles, with a 37° modulation in the topmost 1 or 2 LNO unit cells, driven by strain in the oxygen sublattice for the two opposite polarization directions in PZT. Both oxygen and cation sublattices exhibit strain responses exceeding 10% upon switching. This atomic-layer structural modulation highlights a mechanism for functional oxide heterostructure development, offering pathways for advancements in nonvolatile memory, sensors, and energy-efficient transistors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Site-Selective Modification of Lanthanum Oxychloride to Modulate Halide-Ion Conduction

Design principles for solid-state halide-ion conduction remain poorly defined despite the increasing importance of halide ions as charge carriers in a variety of energy storage and electrochemical computing technologies. Here, we employ a siteselective modification strategy in which aliovalent cations are preferentially introduced at the La 3+ crystallographic site of LaOCl in the 2c Wyckoff position, enabling controlled generation of chloride vacancies and modification of lattice dynamics to enhance chloride-ion conductivity. Aliovalent substitution of La 3+ with Mg 2+ , Ca 2+ , and Sr 2+ generates charge-compensating Cl vacancies while preserving the matlockite crystal structure. X-ray excited optical luminescence measurements with Dy 3+ as a reporter chromophore evidence vacancy-derived midgap electronic states and an extended energy range of radiation-less Auger electron emission corresponding to substantial modification of electronic structure and local electrostatic potentials. Ca alloying at 8−10 at. % increases the chloride-ion conductivity by three- to 4 orders of magnitude as compared to unalloyed LaOCl, whereas comparable amounts of Sr- and Mg-alloying in LaOCl imbue less pronounced conductivity enhancements. Temperature-dependent Raman spectroscopy measurements reveal that Ca- and Sr-alloying substantially soften the La−Cl sublattice and yield a more compliant crystal lattice that can deform to accommodate Cl-ion migration. Structure solutions derived from Rietveld refinements to powder Xray diffraction reveal larger O−La−Cl bond-angle deviations and enhanced out-of-plane cation displacements for Ca- and Sr-alloyed compositions as compared to Mg-alloyed LaOCl. Such local distortions enhance chloride-ion mobility by reshaping and flattening vacancy migration energy landscapes and by modulating lattice dynamics governing anion conduction. We further illustrate that coalloying of Ca with Mg and Sr induces a nonmonotonic conductivity−defect stoichiometry relationship that can be rationalized based on cooperative interactions. Together, these results establish site-selective aliovalent alloying of LaOCl as an effective route to halide-ion solid electrolytes and provide broadly generalizable design principles for site-selective modification to induce vacancy formation and lattice softening to engender facile anion transport

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

Mitigating Band Tailing in Kesterite Solar Absorbers: Ab Initio Quantum Dynamics

Open-circuit voltage deficits are limiting factors in kesterite solar cells. Addressing this issue by suppressing band tailing and nonradiative charge recombination is essential for enhancing the performance. We employ ab initio nonadiabatic molecular dynamics to elucidate the origin of band tailing and charge losses and propose a mitigation strategy. The simulations show that Cu–Zn disorder, associated with antisite defect clusters [Cu Zn +Zn Cu ], is a significant source of band tailing in kesterites, as evidenced by the much larger Urbach energy in disordered than ordered kesterites. Cu–Zn disorder gives rise to new sulfur-centered coordination polyhedra, increases structural inhomogeneity, changes electrostatic potential at sulfur centers, and shifts the S(3p) orbital energy. Differences in the S(3p)/Cu(3d) and S(3p)/Sn(5s) hybridization strengths and the S(3p) orbital energy shift reduce the band gap by 0.37 eV. Furthermore, Cu–Zn disorder enhances vibrational motion of sulfur anions and surrounding cations, increasing band gap fluctuations by 15 meV. The stronger electron–phonon interactions reduce charge carrier lifetimes and limit the kesterite solar cell efficiency. Partial substitution of Zn with Cd facilitates structural ordering and significantly suppresses band tailing, particularly in disordered systems. The improvement can be attributed to the larger atomic radius and mass of Cd, which weakens bonding around the anion, suppresses S-related vibrations within the covalent tetrahedra, and reduces nonadiabatic coupling, thereby increasing charge carrier lifetimes. The reported results establish the key influence of cation disorder on band tailing and reduced charge carrier lifetimes in kesterites and highlight cation disorder engineering as a strategy to achieve high-efficiency kesterite solar cells.

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

Polarization-driven band topology evolution in twisted MoTe 2 and WSe 2

Motivated by recent experimental observations of opposite Chern numbers in R-type twisted MoTe 2 and WSe 2 homobilayers, we perform large-scale density-functional-theory calculations with machine learning force fields to investigate moiré band topology across a range of twist angles in both materials. We find that the Chern numbers of the moiré frontier bands change sign as a function of twist angle, and this change is driven by the competition between moiré ferroelectricity and piezoelectricity. Our large-scale calculations, enabled by machine learning methods, reveal crucial insights into interactions across different scales in twisted bilayer systems. The interplay between atomic-level relaxation effects and moiré-scale electrostatic potential variation opens new avenues for the design of intertwined topological and correlated states, including the possibility of mimicking higher Landau level physics in the absence of magnetic field.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Study of graphene p-n junctions formed by the electrostatic modification of the SiO2 substrate

Abstract We study the transport properties of mm-scale CVD graphene p-n junctions, which are formed in a single gated graphene field effect transistor configuration. Here, an electrical-stressing-voltage technique served to modify the electrostatic potential in the SiO 2 /Si substrate and create the p-n junction. We examine the transport characteristics about the Dirac points that are localized in the perturbed and unperturbed regions in the graphene channel and note the quantitative differences in the Hall effect between the perturbed and unperturbed regions. The results also show that the longitudinal resistance is highly sensitive to the external magnetic field when the Hall bar device operates as a p-n junction.

Science & Technology - Other Topics↗
DOE OSTI2024

Anomalous transport of multi-species edge plasma with the generalized Hasegawa–Wakatani model and the FLR effects

Anomalous transport of multi-species plasma is considered with the generalized Hasegawa–Wakatani model [A. R. Knyazev and S. I. Krasheninnikov, Phys. Plasmas 31, 012502 (2024)] further extended to incorporate the Finite Larmor Radius (FLR) effects. By introducing the “associated” enstrophy, it is shown that with no FLR effects (where anomalous transport of all ion species is described as a transport of passive scalars in the turbulent fields of the electrostatic potential and electron density fluctuations) the fluctuating densities of ion species converge to the state where they are linearly proportional to electron density and vorticity fluctuations, which confirm previous numerical findings of [A. R. Knyazev and S. I. Krasheninnikov, Phys. Plasmas 31, 012502 (2024)]. However, in contrast to the “cold” ion approximation, with the FLR effects included, both the plasma turbulence and the dynamics of all ion species become interconnected. Therefore, for simplicity, the FLR effects in this work were considered only for a small “trace” impurity fraction. It is found that for light (neon) “trace” impurity, the FLR effects reduce both anomalous flux and density fluctuations. However, for heavy (tungsten) “trace” impurity, the FLR effects exhibit non-monotonic impact on anomalous transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

On anomalous transport of multi-species plasma associated with the resistive ballooning and resistive drift waves driven turbulence

Anomalous transport of multi-species plasma related to the resistive ballooning and resistive drift wave turbulence is considered in a “cold” ion approximation. It is found that similar to the resistive drift wave turbulence [see A. R. Knyazev and S. I. Krasheninnikov, Phys. Plasmas 31, 012502 (2024); and S. I. Krasheninnikov and R. D. Smirnov, Phys. Plasmas (to be published)] the addition of the ballooning drive does not change the main features of anomalous transport of the multi-species plasma: (i) The transport of all ion species is described as a transport of the passive scalars in the turbulent field of the electrostatic potential and electron density perturbation; (ii) the density of ion species with a larger ratio of the mass to charge has the tendency to the accumulation/depletion in the vortices of plasma flow; and (iii) the cross-field transport of all plasma species (including electrons and ions) is described by the same anomalous transport coefficient.

Physics↗