Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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 OSTI2021

Computational study of the electrostatic potential and charges of multivalent ionic liquid molecules

The electrostatic interactions in ionic liquids (ILs) can be difficult to quantify, but this information is very valuable for improving intermolecular potentials for performing molecular-level simulations. Although polarizable forcefields can provide improved accuracy, the computational efficiency of fixed charge molecular models is often preferred for larger systems. Here, we propose a multi-scale screening approach for analyzing the charge scaling behavior and other electrostatic properties of ILs, which is more attainable than previous analysis methods that rely on the IL crystalline structures. Based on isolated molecules and small clusters, we study the electrostatic potential properties and charges of multivalent IL molecules using density-functional theory calculations and ab initio molecular dynamics simulations. The charges of the IL molecules are calculated and compared using several different charge estimation techniques. Although there are some differences among the estimated partial charges on the atoms, the total anion and cation charges are very similar. In conclusion, the charge scaling factors are substantial, and they are only slightly affected by the IL cluster size model.

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

When do short-range atomistic machine-learning models fall short?

In this work, we explore the role of long-range interactions in atomistic machine-learning models by analyzing the effects on fitting accuracy, isolated cluster properties, and bulk thermodynamic properties. Such models have become increasingly popular in molecular simulations given their ability to learn highly complex and multi-dimensional interactions within a local environment; however, many of them fundamentally lack a description of explicit long-range interactions. In order to provide a well-defined benchmark system with precisely known pairwise interactions, we chose as the reference model a flexible version of the Extended Simple Point Charge (SPC/E) water model. Our analysis shows that while local representations are sufficient for predictions of the condensed liquid phase, the short-range nature of machine-learning models falls short in representing cluster and vapor phase properties. These findings provide an improved understanding of the role of long-range interactions in machine learning models and the regimes where they are necessary.

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

Periodic GFN1-xTB Tight Binding: A Generalized Ewald Partitioning Scheme for the Klopman–Ohno Function

A novel formulation is presented for the treatment of electrostatics in the periodic GFN1-xTB tight-binding model. Periodic GFN1-xTB is hindered by the functional form of the second-order electrostatics, which only recovers Coulombic behavior at large interatomic distances and lacks a closed-form solution for its Fourier transform. We address this by introducing a binomial expansion of the Klopman–Ohno function to partition short- and long-range interactions, enabling the use of a generalized Ewald summation for the solution of the electrostatic energy. This approach is general and is applicable to any damped potential of the form |R n + c| –m . Benchmarks on the X23 molecular crystal dataset and a range of prototypical bulk semiconductors demonstrate that this systematic treatment of the electrostatics eliminates unphysical behavior in the equation of state curves. In the bulk systems studied, we observe a mean absolute error in total energy of 35 meV/atom, comparable to the machine-learned universal force field, M3GNet, and sufficiently precise for structure relaxation. These results highlight the promising potential of GFN1-xTB as a universal tight-binding parametrization.

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

Role of water model on ion dissociation at ambient conditions

In this work, we study ion pair dissociation in water at ambient conditions using a combination of classical and ab initio approaches. The goal of this study is to disentangle the sources of discrepancy observed in computed potentials of mean force. In particular, we aim to understand why some models favor the stability of solvent-separated ion pairs vs contact ion pairs. We found that some observed differences can be explained by non-converged simulation parameters. However, we also unveil that for some models, small changes in the solution density can have significant effects on modifying the equilibrium balance between the two configurations. We conclude that the thermodynamic stability of contact and solvent-separated ion pairs is very sensitive to the dielectric properties of the underlying simulation model. In general, classical models are very robust in providing a similar estimation of the contact ion pair stability, while this is much more variable in density functional theory-based models. The barrier to transition from the solvent-separated to contact ion pair is fundamentally dependent on the balance between electrostatic potential energy and entropy. This reflects the importance of water intra- and inter-molecular polarizability in obtaining an accurate description of the screened ion–ion interactions.

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

Charge fluctuations from molecular simulations in the constant-potential ensemble

We revisit the statistical mechanics of charge fluctuations in capacitors. In constant-potential classical molecular simulations, the atomic charges of electrode atoms are treated as additional degrees of freedom which evolve in time so as to satisfy the constraint of fixed electrostatic potential for each configuration of the electrolyte. The present work clarifies the role of the overall electroneutrality constraint, as well as the link between the averages computed within the Born–Oppenheimer approximation and that of the full constant-potential ensemble. This allows us in particular to derive a complete fluctuation–dissipation relation for the differential capacitance, that includes a contribution from the charge fluctuations (around the charges satisfying the constant-potential and electroneutrality constraints) also present in the absence of an electrolyte. We provide a simple expression for this contribution from the elements of the inverse of the matrix defining the quadratic form of the fluctuating charges in the energy. We then illustrate numerically the validity of our results, and recover the expected continuum result for an empty capacitor with structureless electrodes at large inter-electrode distances. By considering a variety of liquids between graphite electrodes, we confirm that this contribution to the total differential capacitance is small compared to that induced by the thermal fluctuations of the electrolyte.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Structure, dynamics, and electrochemistry of choline chloride/ethylene glycol eutectic solvents at an electrode surface explored by molecular dynamics simulations

Choline chloride and ethylene glycol mixtures with 1:2, 1:4, and 1:6 molar ratios on the surfaces of graphite and gold electrodes were studied using classical molecular dynamics simulations. Both neutral and charged electrodes were considered. The liquid composition, solvation structure, molecular orientation, and dynamics at the electrode surface are significantly different from those of the bulk liquid. These properties strongly depend on the electrode material and charge density, whereas they are less sensitive to the overall solvent composition. The effect of the electrode on the composition, structure, and orientation of the liquid fades beyond ∼10 Å from the surface of the electrode. This distance corresponds to about two layers of the structured solvent, despite the fact that the layered structure extends to at least five layers or about 25 Å. However, the electrode influences solvent dynamics over a longer distance. The electrochemical properties of the eutectic solvent at both electrode surfaces were also studied. The simulations captured the experimental differential capacitance shapes for both electrode systems, although the magnitudes and exact shapes differ. The simulations further revealed that the solvent in the first solvation layer can both overscreen and underscreen the electrode charges depending on the electrode material and electrode potential.

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

Driving force and pathway in polyelectrolyte complex coacervation

There is notable discrepancy between experiments and coarse-grained model studies regarding the thermodynamic driving force in polyelectrolyte complex coacervation: experiments find the free energy change to be dominated by entropy, while simulations using coarse-grained models with implicit solvent usually report a large, even dominant energetic contribution in systems with weak to intermediate electrostatic strength. Here, using coarse-grained, implicit-solvent molecular dynamics simulation combined with thermodynamic analysis, we study the potential of mean force (PMF) in the two key stages on the coacervation pathway for symmetric polyelectrolyte mixtures: polycation–polyanion complexation and polyion pair–pair condensation. We show that the temperature dependence in the dielectric constant of water gives rise to a substantial entropic contribution in the electrostatic interaction. By accounting for this electrostatic entropy, which is due to solvent reorganization, we find that under common conditions (monovalent ions, room temperature) for aqueous systems, both stages are strongly entropy-driven with negligible or even unfavorable energetic contributions, consistent with experimental results. Furthermore, for weak to intermediate electrostatic strengths, this electrostatic entropy, rather than the counterion-release entropy, is the primary entropy contribution. From the calculated PMF, we find that the supernatant phase consists predominantly of polyion pairs with vanishingly small concentration of bare polyelectrolytes, and we provide an estimate of the spinodal of the supernatant phase. Finally, we show that prior to contact, two neutral polyion pairs weakly attract each other by mutually induced polarization, providing the initial driving force for the fusion of the pairs.

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

Mesoporous peptide frameworks engineered from crystallizable collagen-mimetic peptide amphiphiles

The rational design of porous frameworks with tunable pore dimensions and chemical functionalities is a critical step toward their implementation in diverse applications. While traditional porous materials are typically constructed from abiotic components, there is increasing interest in employing biologically derived building blocks (e.g., peptides and proteins) that offer unmatched structural and functional diversity. Here, we report the construction of crystalline mesoporous frameworks that are self-assembled from amphiphilic collagen-mimetic peptides. Comprehensive structural characterization via microscopy, spectroscopy, and computational techniques provides insights into the assembly packing model, in which hexagonally packed channels are interconnected by antiparallel-aligned collagen triple helices via hydrophobic and electrostatic interactions. Lastly, we demonstrate the functional potential of aCMP frameworks through the encapsulation of various molecular guests, including doxorubicin, an anti-cancer drug. Overall, this work establishes a class of mesoporous frameworks, derived from synthetically engineerable peptide conjugates, marking a significant step forward in broadening the architectural scope and application potential of peptide-based materials.

Perez, Anthony R↗
DOE OSTI2023

Electrical Manipulation of Quantum Coherence in a Two-Level Molecular System

We report the manipulation of ultrafast quantum coherence of a two-level single hydrogen molecular system by employing static electric field from the sample bias in a femtosecond terahertz scanning tunneling microscope. A H 2 molecule adsorbed on the polar Cu 2⁢ N surface develops an electric dipole and exhibits a giant Stark effect. An avoided crossing of the quantum state energy levels is derived from the resonant frequency of the single H 2 two levels in a double-well potential. The dephasing time of the initial wave packet can also be changed by applying the electric field. The electrical manipulation for different tunneling gaps in three dimensions allows quantification of the surface electrostatic fields at the atomic scale. Further, our work demonstrated the potential application of molecules as controllable two-level molecular systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Predicting Outcomes of Nanoparticle Attachment by Connecting Atomistic, Interfacial, Particle, and Aggregate Scales

Predicting nanoparticle aggregation and attachment phenomena requires a rigorous understanding of the interplay among crystal structure, particle morphology, surface chemistry, solution conditions, and interparticle forces, yet no comprehensive picture exists. We used an integrated suite of experimental, theoretical, and simulation methods to resolve the effect of solution pH on the aggregation of boehmite nanoplatelets, a case study with important implications for the environmental management of legacy nuclear waste. Real-time observations showed that the particles attach preferentially along the (010) planes at pH 8.5 and the (101) planes at pH 11. To rationalize these results, we established the connection between key physicochemical phenomena across the relevant length scales. Starting from molecular-scale simulations of surface hydroxyl reactivity, we developed an interfacial-scale model of the corresponding electrostatic potentials, with subsequent particle-scale calculations of the resulting driving forces allowing successful prediction of the attachment modes. Finally, we scaled these phenomena to understand the collective structure at the aggregate-scale . Here, our results indicate that facet-specific differences in surface chemistry produce heterogeneous surface charge distributions that are coupled to particle anisotropy and shape-dependent hydrodynamic forces, to play a key role in controlling aggregation behavior.

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

Molecular bridge-mediated ultralow-power gas sensing

Abstract We report the electrical detection of captured gases through measurement of the quantum tunneling characteristics of gas-mediated molecular junctions formed across nanogaps. The gas-sensing nanogap device consists of a pair of vertically stacked gold electrodes separated by an insulating 6 nm spacer (~1.5 nm of sputtered α-Si and ~4.5 nm ALD SiO 2 ), which is notched ~10 nm into the stack between the gold electrodes. The exposed gold surface is functionalized with a self-assembled monolayer (SAM) of conjugated thiol linker molecules. When the device is exposed to a target gas (1,5-diaminopentane), the SAM layer electrostatically captures the target gas molecules, forming a molecular bridge across the nanogap. The gas capture lowers the barrier potential for electron tunneling across the notched edge region, from ~5 eV to ~0.9 eV and establishes additional conducting paths for charge transport between the gold electrodes, leading to a substantial decrease in junction resistance. We demonstrated an output resistance change of >10 8 times upon exposure to 80 ppm diamine target gas as well as ultralow standby power consumption of <15 pW, confirming electron tunneling through molecular bridges for ultralow-power gas sensing.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Ion density waves driving the formation of filamentary dust structures

The PlasmaKristall-4 (PK-4) experiment on the International Space Station allows for the study of the 3-dimensional interaction between plasma and dust particles. Previous simulations of the PK-4 environment have discovered fast moving ionization waves in the dc discharge. These ionization waves vary the plasma parameters by up to an order of magnitude, which may affect the mechanisms responsible for the self-organization of chains seen in the PK-4 experiment. Here, we adapt a molecular dynamics simulation to employ temporally varying plasma conditions in order to investigate the effect on the dust charging and electrostatic potential. In order to describe the differences between the average of the plasma conditions and the time-varying plasma condition, we present a model to reproduce the interaction that takes into account the negative potential from the dust grain and the positive potential from the ion wake.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Electrostatic Potential of Functional Cations as a Predictor of Hydroxide Diffusion Pathways in Nanoconfined Environments of Anion Exchange Membranes

Nanoconfined anion exchange membranes (AEMs) play a vital role in emerging electrochemical technologies. The ability to control dominant hydroxide diffusion pathways is an important goal in the design of nanoconfined AEMs. Such control can shorten hydroxide transport pathways between electrodes, reduce transport resistance, and enhance device performance. In this work, we propose an electrostatic potential (ESP) approach to explore the effect of the polymer electrolyte cation spacing on hydroxide diffusion pathways from a molecular perspective. By exploring cation ESP energy surfaces and validating outcomes through prior ab initio molecular dynamics simulations of nanoconfined AEMs, we find that we can achieve control over preferred hydroxide diffusion pathways by adjusting the cation spacing. The results presented in this work provide a unique and straightforward approach to predict preferential hydroxide diffusion pathways, enabling efficient design of highly conductive nanoconfined AEM materials for electrochemical technologies.

25 ENERGY STORAGE↗
DOE OSTI2026

How Does Water Dissociation Work in Bipolar Membranes?

Bipolar membranes (BPMs) create counteracting spatial gradients of pH and electrostatic potential in electrochemical systems, enabling applications in pH regulation, electrocatalysis, and separations. At the polarized junction of a BPM the water dissociation (WD, 2H2O ⇌ H3O+ + OH-) reaction can be driven, but it remains poorly understood. In this Perspective, we integrate molecular insights from bulk-water autoionization and the associated field effects with continuum descriptions of BPM electrostatics and experimental WD kinetic analyses to describe possible mechanisms of voltage-driven WD. Pristine BPM junctions highlight both the limits of primarily electric-field-driven WD and the practical challenges of junction stability at extreme reverse bias. Introducing heterogeneous catalyst layers, commonly metal oxides and graphene oxides, accelerates WD by orders of magnitude through hypothesized coupled effects in which surface acid-base functionality and high-density hydroxyl sites mediate proton-transfer steps, and catalyst mobile electronic/ionic charges redistribute the junction electric potential drop to shape the local electric fields and reactive microenvironments. Kinetic analyses suggest two regimes of heterogeneous WD mechanism, including field-driven ordering of interfacial water and a Second-Wien-Effect dissociation-barrier lowering. We conclude by defining the key unknown variables (local pH, electrostatic potential, catalyst charge state and relationships among mechanisms) and outlining experimental and multiscale modeling strategies needed for predictive WD catalysis and for controlling related ion-transfer reactions.

Wu, Yifan↗
DOE OSTI2020

Modeling Adsorption in Silica Pores via Minkowski Functionals and Molecular Electrostatic Moments

Capillary condensation phenomena are important in various technological and environmental processes. Using molecular simulations, we study the confined phase behavior of fluids relevant to carbon sequestration and shale gas production. As a first step toward translating information from the molecular to the pore scale, we express the thermodynamic potential and excess adsorption of methane, nitrogen, carbon dioxide, and water in terms of the pore’s geometric properties via Minkowski functionals. This mathematical reconstruction agrees very well with molecular simulations data. Our results show that the fluid molecular electrostatic moments are positively correlated with the number of adsorption layers in the pore. Moreover, stronger electrostatic moments lead to adsorption at lower pressures. These findings can be applied to improve pore-scale thermodynamic and transport models.

42 ENGINEERING↗
DOE OSTI2025

Influence of counterion substitution on the properties of imidazolium-based ionic liquid clusters

Due to their unique physiochemical properties that may be tailored for specific purposes, ionic liquids (ILs) have been investigated for various applications, including chemical separations, catalysis, energy storage, and space propulsion. The different cations and anions comprising ILs may be selected to optimize a range of desired properties, such as thermal stability, ionic conductivity, and volatility, leading to the designation of certain ILs as designer “green” solvents. The effect of counterions on the properties of ILs is of both fundamental scientific interest and technological importance. Herein, we report a systematic experimental and theoretical investigation of the size, charge, stability toward dissociation, and geometric/electronic structure of 1-ethyl-3-methyl imidazolium (EMIM)-based IL clusters containing two different atomic counterions (i.e., bromide [Br − ] and iodide [I − ]). This work extends our studies of EMIM + cations with atomic chloride (Cl − ) and molecular tetrafluoroborate (BF 4 − ) anions reported previously by Baxter et al. [Chem. Mater. 34, 2612 (2022)] and Zhang et al . [J. Phys. Chem. Lett. 11, 6844 (2020)], respectively. Distributions of anionic IL clusters were generated in the gas phase using electrospray ionization and characterized by high mass resolution mass spectrometry, energy-resolved collision-induced dissociation, and negative ion photoelectron spectroscopy experiments. The experimental results reveal anion-dependent trends in the size distribution, relative abundance, ionic charge state, stability toward dissociation, and electron binding energies of the IL clusters. Complementary global optimization theory provides molecular-level insights into the bonding and electronic structure of a selected subset of clusters, including their low energy structures and electrostatic potential maps, and how these fundamental characteristics are influenced by anion substitution. Collectively, our findings demonstrate how the fundamental properties of ILs, which determine their suitability for many applications, may be tuned by substituting counterions. These observations are critical in the sub-nanometer cluster size regime where phenomena do not scale predictably to the bulk phase, and each atom counts toward determining behavior.

cluster↗
DOE OSTI2021

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↗