Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electrostatic”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Structure–Capacitance Relationships of Graphene/Ionic Liquid Electrolyte Double Layers

The differential capacitance profile of electrochemical interfaces reflects the physical properties of the double layer. For carbon electrodes and ionic liquid based electrolytes, these capacitance profiles are not fully understood. In this work, we utilize constant voltage molecular dynamics simulations to compute differential capacitance profiles of ionic liquids [BMIm + ][BF 4 - ] and [BMIm + ][TFSI - ] mixed with acetonitrile and 1,2-dichloroethane, at model graphene electrodes. We find that both pure and 10% mole fraction ionic liquid electrolytes exhibit camel shaped capacitance profiles with two peaks on either side of a minimum centered at the potential of zero charge. This profile shape results from electric field induced rearrangement of ion structure within the inner layer closest to the electrode interface. At low potential, the ionic liquid inner layer is concentrated with non-polar tri-fluoromethyl and butyl functional groups of the anions and cations, corresponding to the minimum of the capacitance profiles. With increasing voltage, electrostatic interactions of polar/charged functional groups with the electrode surface compete with these non-polar interactions, leading to ion rearrangement that increases the inner layer charge density and results in higher capacitance. After the ion restructuring is complete, the response saturates and capacitance diminishes. The presence of organic solvent significantly changes the composition of the inner layer. For example, strong non-polar interactions between dichloroethane molecules and the graphene surface substantially block ion/electrode contact at moderate potentials. Overall, our simulations highlight the dynamic nature of the inner region of organic electrolyte double layers, and the sensitive dependence on electrolyte composition and applied voltage.

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

Preferential Direction of Electron Transfers at a Dye–Metal Oxide Interface with an Insulating Fluorinated Self-Assembled Monolayer and MgO

Slowing non-desirable electron transfer reactions at metal oxide-dye interfaces is important for many technologies. Specifically, after an interfacial photoinduced charge separation event at a metal oxide-dye interface, it is critically important to limit the rate of electron transfer reactions back to sensitizers and to limit electron transfer reactions between the electrolyte and the metal oxide. Ruthenium-based dyes at metal oxide interfaces are widely used in many fields; however, these dyes often have poor surface insulation resulting in fast recombination kinetics with transition metal-based redox shuttles in an electrolyte. This work explores two semiconductor surface modification strategies designed to minimize recombination events of electrons in TiO 2 with oxidized redox shuttles using a fluorinated siloxane insulator (PFTS) and a metal oxide insulator (MgO) with a well-known Ru dye, B11. Additionally, the influence of these treatments on the rate and duration of photoinduced interfacial charge separation at the TiO 2 -dye interface was examined. The TiO 2 -dye-redox shuttle systems were studied via dye-sensitized solar cell current-voltage curve, incident photon-to-current conversion efficiency, small-modulated photovoltage transient, time-correlated single photon counting, and transient absorption spectroscopy measurements. MgO was found to decrease the rate of the electron transfer reaction from the metal oxide to the electrolyte, decrease the rate of the electron transfer reaction to the oxidized dye from TiO 2 , increase the electron transfer reaction rate from a reduced redox shuttle to an oxidized dye, and decrease the electron injection rate from the photoexcited dye to TiO2. Interestingly, PFTS was found to desirably improve on these rates relative to MgO or untreated TiO 2 . Here, a model based on electrostatic interactions is presented to explain the exceptional behavior of PFTS with DFT computational analysis of PFTS supporting this model.

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

Adsorption and Catalytic Activity of Gold Nanoparticles in Mesoporous Silica: Effect of Pore Size and Dispersion Salinity

The assembled state of nanoparticles (NPs) within porous matrices plays a governing role in directing their biological, electronic, and catalytic properties. However, the effects of the spatial confinement and environmental factors, such as salinity, on the NP assemblies within the pores are poorly understood. In this study, we use adsorption isotherms, spectrophotometry, and small-angle neutron scattering to develop a better understanding of the effect of spatial confinement on the assembled state and catalytic performance of gold (Au) NPs in propylamine-functionalized SBA-15 and MCM-41 mesoporous silica materials (mSiO 2 ). We carry out a detailed investigation of the effect of pore diameter and ionic strength on the packing and spatial distribution of AuNPs within mSiO 2 to get a comprehensive insight into the structure, functioning, and activity of these NPs. We demonstrate the ability of the adsorbed AuNPs to withstand aggregation under high salinity conditions. We attribute the observed preservation of the adsorbed state of AuNPs to the strong electrostatic attraction between oppositely charged pore walls and AuNPs. The preservation of the structure allows the AuNPs to retain their catalytic activity for a model reaction in high salinity aqueous solution, here, the reduction of p-nitrophenol to p-aminophenol, which otherwise is significantly diminished due to bulk aggregation of the AuNPs. This fundamental study demonstrates the critical role of confinement and dispersion salinity on the adsorption and catalytic performance of NPs.

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

Random Walks and Sticky Surfaces: Single-Molecule Measurements of Solute Diffusion in Ethanol/Water-Filled Anodic Alumina Nanopores

Nanoporous anodic aluminum oxide (AAO) membranes are now being explored for use in advanced chemical separations, including in the dehydration of biofuels such as ethanol. Optimization of membrane performance requires an in-depth understanding of how solvent mixtures and solutes behave under nanoconfinement. In this work, the diffusion of rhodamine B (RhB) dye through 10 nm and 20 nm AAO nanopores filled with a series of ethanol/water mixtures (0 - 33% water) is explored by fluorescence correlation spectroscopy (FCS). RhB was found to diffuse through the pores by two distinct mechanisms with mean diffusion coefficients, $D_f$ and $D_s$, reflecting fast and slow diffusive motions, respectively, with values that differ by nearly two orders of magnitude. Further, both $D_f$ and $D_s$ increased with pore size and were significantly smaller than $D_b$, the RhB diffusion coefficient in bulk liquid. Mean $D_f$ values follow a composition-dependent trend that closely mimics the viscosity dependence of $D_b$. Additional slowing of fast RhB diffusion is attributed to both hydrodynamic drag and electrostatic interactions with the nanopore surface. The mean $D_s$ values exhibit a different trend with increasing water content, revealing an increase in $D_s$ and a decrease in the contributions of slow diffusion to the observed dynamics. The fluorescence time transient data used in the analysis show that the slow diffusion process is strongly hindered and likely involves frequent adsorption of RhB to the pore surfaces. These results provide new insights into the detailed molecular-level mechanisms of mass transport in nanoporous AAO membranes.

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

Electric Field Effects on Water and Ion Structure and Diffusion at the Orthoclase (001)–Water Interface

Understanding the electrochemical properties of mineral–water interfaces tends to rely upon electrical double layer (EDL) models, but these models are based on the assumption that electrostatic equilibrium is constantly maintained. In reality, interfacial reactions, ion diffusion, and their electrochemical signatures are based in nonequilibrium conditions of locally or globally imbalanced electrical fields where current EDL models have limited purview. In this work, we performed molecular dynamics (MD) simulations of the orthoclase (001) surface in contact with a 1 M NaCl aqueous solution under various electric fields, to explore the interplay between EDL structure and dynamics when perturbed by electric fields of different direction and strength, by confinement, and by different distributions of structural surface charge. The simulations showed that confinement between two opposing (001) surfaces led to the development of an induced field when the applied field was perpendicular to the surfaces and, as a result, to ionic diffusion coefficients that were independent of electric field strength. In contrast, when the applied field was parallel to the surfaces, confinement resulted in ionic diffusion coefficients that were more strongly dependent on the magnitude of the electric field than in bulk water. Differences in the density and distribution of aluminol groups on the two surfaces had a significant impact on how the interfacial structure and dynamics varied in the presence of an electric field. Notably, these differences resulted in an electro-osmotic flow with opposite directions at the two surfaces under parallel applied electric field. Overall, the MD simulations highlighted the importance of considering atomic-level structure and heterogeneities when developing models of the electrochemical properties of mineral–water interfaces.

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

Heavy Versus Light Lanthanide Selectivity for Graphene Oxide Films is Concentration Dependent

Rare earths are important materials in various technologies such as catalysis and optoelectronics. Graphene oxide (GO) is a promising material for separation applications, including the isolation of lanthanides from complex mixtures. Previous works using fatty acid monolayers have demonstrated preferential heavy versus light lanthanide adsorption, which has been attributed to differences in lanthanide ion size. In this work, we used interfacial X-ray fluorescence measurements to reveal that GO thin films at the air/water interface have no lanthanide selectivity for dilute subphases. However, at high subphase concentrations, ~8 times more Lu is adsorbed than La. By comparing the GO results with an ideal monolayer with a carboxylic acid headgroup, arachidic acid (AA), we demonstrate that the number of Lu ions adsorbed to GO is significantly higher than the number expected to compensate for the surface charge. Vibrational sum frequency generation (SFG) spectroscopy results on both GO thin films and AA monolayers reveal a red-shifted SFG signal in the OH region, which we attribute to partial dehydration of the adsorbed ions and carboxylic acid headgroups. Liquid surface X-ray reflectivity data show that the GO thin film structure does not significantly change between the very dilute and concentrated subphases. We speculate that the functional groups of both GO and AA facilitate cation dehydration, which is essential for ion adsorption. Heavy lanthanide Lu has stronger ion–ion correlations that can overcome the electrostatic repulsion between cations at higher concentrations compared to light lanthanide La, meaning GO and AA can exhibit apparent overcharge with Lu. Lastly, the layered structure of the GO films and reactive chemical nature of GO itself can accommodate ion adsorption.

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

Ion Solvation-Driven Liquid–Liquid Phase Separation in Divalent Electrolytes with Miscible Organic Solvents

Liquid–liquid phase separation (LLPS) is a common phenomenon, but LLPS of electrolytes prepared in miscible organic solvents is rarely documented. Here we report four cases of LLPS that occur in MgTFSI 2 or ZnTFSI 2 electrolytes in mixed organic solvents. The conditions for the formation of this LLPS share four common features: cations with high charge density; bulky anions with low charge density; a strongly coordinating solvent with high compressibility; and a relatively weak coordinating solvent with low relative permittivity and high mobility. With these conditions, the cations tend to draw the strongly coordinating solvents together to form a densely-packed, energetically-favorable ion solvation region, while the cosolvents with low relative permittivity tend to repel ions to form an ion-depleted dilute phase. Furthermore, in the dilute upper layer with more than 300 solvent molecules per cation, ion pairs and large ion aggregates are clearly evidenced, due to the incomplete screening of electrostatic interactions by the weaker cosolvent. This LLPS driven by ion solvation may be more common in multivalent electrolytes and is a design consideration with mixed solvents that should not be overlooked.

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

Effect of Flexibility in Molecular Simulations of Carbon Dioxide Adsorption and Diffusion in a Cuprous Triazolate Framework

Using fixed atom force fields to model gas adsorption in flexible metal organic frameworks (MOFs) is known to pose difficulties in accurately reproducing and predicting experimental results. This paper studies the difference in accuracy between flexible and fixed atom force fields in reproducing CO 2 adsorption measurements in MAF-2 ([Cu(etz)]∞ (MAF-2, Hetz) 3,5-diethyl-1,2,4-triazole), an NbO-type triazolate scaffold with a bcu cavity system and attached ethyl groups. The flexible force field used to run the hybrid molecular dynamics and grand canonical Monte Carlo calculations were generated using the QuickFF software incorporating van der Waals parameters from the Universal Force Field (UFF) and density derived electrostatic and chemical (DDEC) charges. The fixed atom force field used was composed of UFF van der Waals parameters together with DDEC charges. The calculations were run at 298 K and at pressures of 0.1, 0.3, 0.5, 0.8, and 1 bar. It was observed that for this MOF the rigid force field overpredicted gas adsorption, whereas the flexible force field values closely matched experimental data. In the flexible structure, the freely moving ethyl groups of MAF-2 hindered adsorption, reducing the interaction energy between CO 2 and the N atoms of the triazolate framework as well as reducing the size of the largest cavity diameter. Here, the combination of these factors led to improved prediction of adsorption values with the flexible forcefield as compared to the rigid forcefield, demonstrating the need for modeling MOFs in a way more indicative of their behavior.

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

Understanding Anomalous Gas-Phase Peak Shifts in Dip-and-Pull Ambient Pressure XPS Experiments

Dip-and-pull ambient pressure X-ray photoelectron spectroscopy (AP-XPS) holds promise to uncover elementary processes of (photo)electrochemistry. We show, however, that the sample for such experiments should preferably be nonporous and the potential on the surface homogeneous. We carried out dip-and-pull AP-XPS experiments on a hematite thin film sample under the photoelectrochemical oxygen evolution reaction (OER) and find unexpected O 1s core level shifts. Upon electrochemical biasing under simulated solar light illumination, the gas-phase water peak shifted more than the electrolyte peak. To uncover the origin of the unexpected larger shift of the gas-phase peak, we performed electrostatic simulations using COMSOL, to map the potential field in the relevant volume between the sample and the first aperture of the XPS spectrometer. A number of geometric models were considered. We find that when the potential on the sample surface is inhomogeneous, e.g., with ionically isolated electrolyte patches, the gas-phase peak of the spectrum can shift more than the peak due to the electrolyte film. This suggests that at the measured sample position, the local potential was not as set by the potentiostat. Despite this, we find reversible consumption and replenishment of hydroxide in the spectra, which, due to OH– being the reactant of the OER in alkaline electrolyte, makes sense chemically. We propose that this is linked to OH– diffusion across the measured sample position, driven by the potential dependent consumption and replenishment of OH– at the nearby well-connected surface regions.

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

Unveiling the Role of Surface Ir-Oxo Species in O 2 Evolution at IrO 2 Electrocatalysts via Embedded Cluster Multireference Calculations

Understanding the mechanisms driving the oxygen evolution reaction (OER) on iridium oxide (IrO 2 )-based catalysts is essential to improving their performance and enabling an actual scale-up of water-splitting photoelectrochemical cells. The mechanistic pathways at IrO 2 interfaces have been extensively investigated computationally using density functional theory (DFT), which predicts a high-energy barrier for the last step of the OER of molecular oxygen detachment and release from the catalyst surface. Nevertheless, surface O 2 over- and under-binding results by standard generalized gradient approximation and hybrid density functionals, respectively, call for further analysis of this crucial step via multireference methods. Aiming at unveiling the nature of such a barrier, we hereby address the formation of O 2 from the most-stable IrO 2 (110) surface with both periodic DFT and an electrostatic embedded cluster approach at the n-electron valence-state perturbation theory. With this multireference approach, we find a value for the aforementioned energy barrier that is much closer to experimental indications than DFT ones. Here, an in-depth analysis of the involved molecular orbitals suggests that the origin of this barrier is related to the breaking of a π interaction between O 2 and Ir surface atom and to a significant additional O 2 interaction with adjacent electrophilic Ir-oxo species, which is present under experimental operating conditions. Besides shedding light on the mechanism of the OER on IrO 2 , these findings point out the importance of multireference methods for dissecting complex reactions at electrocatalytic interfaces and pave the route for further investigations with effective embedding approaches.

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

Insights into the Oxidation Mechanism of Vivianite to Metavivianite from First-Principles Calculations

Vivianite, a hydrous ferrous iron-phosphate mineral (Fe 3 (PO 4 ) 2 ·8H 2 O), readily oxidizes in contact with air yielding the less hydrous mixed-valent iron-phosphate mineral metavivianite. This topotactic transformation nominally occurs by oxidative dehydrogenation in which outgoing electrons from the iron sublattice are charge compensated by hydrolysis of structural water. However, the details of this internal charge balancing mechanism that allows the structure to remain electrostatically stable remain unknown. Here, in this study, we use density functional theory (DFT) calculations and ab initio thermodynamics (AIT) to evaluate the energetics of this process in terms of hydrogen release as a function of environmental variables such as partial pressure and temperature. The results show a thermodynamic driving force for vivianite phase transformation as oxidation progresses that is triggered by its rigid structure that has a limited accommodation for hydrogen vacancies. In contrast, metavivianite has a more flexible lattice and hydrogen bond network that stabilizes these hydrogen defects. Metavivianite is shown to be a stable intermediate for 66% residual Fe 2+ down to 33% Fe 2+ , below which other phases/structures such as santabarbarite should be more thermodynamically favorable. Our study provides a basis for experimental tests of our mechanistic findings and helps fill a basic knowledge gap about the solid-state process that defines how vivianite interacts with its surrounding environment.

Sassi, Michel [Pacific Northwest National Laborato↗
DOE OSTI2020

Nanoparticle-Enhanced Plasma Discharge Using Nanosecond High-Voltage Pulses

By discharging nanosecond high-voltage (5 kV) pulses across an insulating substrate containing Au, Pt, or Cu nanoparticles, a 3 order of magnitude (1000×) enhancement in the generation of plasma can be achieved through local field enhancement on the surface of the nanoparticles. The low-temperature nature of this transient plasma is crucial to maintaining the structural integrity of these delicate nanoparticles. These nanoparticles provide up to a 1000-fold enhancement in the generation of the plasma, which is localized to the surface of the nanoparticles where it is potentially useful (e.g., for catalysis). Additionally, we performed both time-domain and frequency-domain calculations of the electromagnetic response of the nanoparticles based on high-resolution transmission electron microscope (HRTEM) images, which show local field enhancement of the nanosecond high-voltage pulse on the order of 3×. Since the plasma initiation depends exponentially on the peak electric field strength, this 3-fold increase in the local electric field can result in a several orders of magnitude increases in the generation of plasma at a given applied external field strength. In order to rule out plasmon-resonance enhancement, which is often associated with small metal nanoparticles, we performed finite difference time domain (FDTD) simulations in the optical frequency domain, which show that the effect of plasmon resonance is negligible for Pt nanoparticles. We therefore attribute the nanoparticle-based enhancement to the generation of plasma (an electrostatic effect) rather than enhanced coupling of light from the near field to the far field via the plasmon resonance phenomenon (an optical effect).

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ion Pairing and Redissociaton in Low-Permittivity Electrolytes for Multivalent Battery Applications

Detailed speciation of electrolytes as a function of chemical system and concentration provides the foundation for understanding bulk transport as well as possible decomposition mechanisms. In particular, multivalent electrolytes have shown a strong coupling between anodic stability and solvation structure. Furthermore, solvents that are found to exhibit reasonable stability against alkaline-earth metals generally exhibit low permittivity, which typically increases the complexity of the electrolyte species. To improve our understanding of ionic population and associated transport in these important classes of electrolytes, the speciation of Mg(TFSI)2 in monoglyme and diglyme systems is studied via a multiscale thermodynamic model using first-principles calculations for ion association and molecular dynamics simulations for dielectric properties. The results are then compared to Raman and dielectric relaxation spectroscopies, which independently confirm the modeling insights. We find that the significant presence of free ions in the low-permittivity glymes in the concentration range from 0.02 to 0.6 M is well-explained by the low-permittivity redissociation hypothesis. Here, salt speciation is largely dictated by long-range electrostatics, which includes permittivity increases due to polar contact ion pairs. The present results suggest that other low-permittivity multivalent electrolytes may also reach high conductivities as a result of redissociation.

25 ENERGY STORAGE↗
DOE OSTI2020

Spatially Resolved Investigation of Mixed Valence and Insulator-to-Metal Transition in an Organic Salt

Using scanning tunneling microscopy/spectroscopy (STM/STS), we investigate the evolution of electronic structures across the boundaries of 7,7,8,8-tetracyanoquinodimethane (TCNQ) and K-TCNQ assemblies on a weakly interacting substrate. Despite the semiconducting/insulating nature of TCNQ (TCNQ 0 ) and K-TCNQ (TCNQ -1 ), a continuum metallic-like density of states extending deep (~1.5 nm) into the TCNQ assembly is observed near domain boundary. We attribute the formation of these states to the abrupt change of molecular valence, which perturbs the electrostatics of the junction and creates local electric fields as evidenced by the band bending near domain boundary. To the best of our knowledge, this study provides the first microscopic understanding of the crucial physics occurring near domain boundaries of mixed valence in K-TCNQ, or broadly speaking charge transfer complexes, which highlights these boundaries as potential “weak” points to initiate the electric field induced insulator-to-metal transition.

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

Connecting Hindered Transport in Porous Media across Length Scales: From Single-Pore to Macroscopic

Hindered mass transport is widely observed in various porous media, including filtration membranes and porous rocks. However, there is no universal model capable of predicting transport in porous media due to the heterogeneity of porous structures and the complexity of the underlying microscopic mechanisms, including electrostatic and hydrodynamic effects. Here, we used a highly ordered porous medium (i.e. a silica inverse opal film) as a model system to directly explore the effects of geometric parameters (i.e. pore size, pore throat size and tracer particle size) and microscopic interaction parameters (e.g. controlled by ionic strength) on nanoparticle transport in porous environments using three-dimensional single-particle tracking. We found a linear scaling relation between the macroscopic diffusion coefficient and microscopic diffusion behaviors and a combination of parameters associated with microscopic features and phenomena, including both geometric effects and particle-level interactions. The proportionality coefficient relating micro and macro behaviors was complex and related to the connectivity of the matrix and pore-size variation, which could lead to tortuous diffusion pathways, hindering macroscopic transport. Furthermore, these findings provide important insights into nanoparticles and macromolecules transport in various porous media.

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

Unique Orientation of the Solid–Solid Interface at the Janus Particle Boundary Induced by Ionic Liquids

This study reveals the unique role on Janus particles of the solid–solid interface at the boundary in determining particle interactions and assembly. In an aqueous ionic liquid (IL) solution, Janus spheres adopt intriguing orientations with their boundaries pinned on the glass substrate. In addition, it was further discovered that the orientation was affected by the particle amphiphilicity as well as the chemical structure and concentration of the IL. Further characterization suggests that the adsorption on the hydrophilic side is due to both an electrostatic interaction and hydrogen bonding, while adsorption on the hydrophobic side is due to hydrophobic attraction. Through the concerted interplay of all these interactions, the amphiphilic boundary may attract an excessive amount of IL cations, which guide the unique orientations of the Janus spheres. The results highlight the importance of the Janus boundary that has not been recognized previously. Adsorption at the solid–solid interfaces may inspire new applications in areas such as separation and catalysis.

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

Size-Dependent Onset of Nitric Acid Dissociation in Cs + ·(HNO 3 )(H 2 O) n =0–11 Clusters at 20 K

We report the water-mediated charge separation of nitric acid upon incorporation into size-selected Cs + ∙(HNO 3 )(H 2 O) n=0-11 clusters at 20 K. Dramatic spectral changes are observed in the range n=7-9 that are traced to the formation of many isomeric structures associated with intermediate transfer of the acidic proton to the water network. This transfer is complete by n=10, which exhibits much simpler vibrational band patterns consistent with those expected for a tri-coordinated hydronium ion (the Eigen motif) along with the NO stretching bands predicted for a hydrated NO 3 – anion that is directly complexed to the Cs + cation. Theoretical analysis of the n=10 spectrum indicates that the dissociated ions adopt a solvent-separated ion-pair configuration such that the Cs + and H 3 O + cations flank the NO 3 – anion in a microhydrated salt bridge. In conclusion, this charge separation motif is evidently assisted by the electrostatic stabilization of the product NO 3 – /H 3 O + ion pair by the proximal metal ion.

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

Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data

Ionic solvation phenomena in liquids involve intense interactions in the inner solvation shell. For interactions beyond the first shell, the ion–solvent interaction energies result from the sum of many smaller-magnitude contributions that can still include polarization effects. Deep neural network (DNN) methods have recently found wide application in developing efficient molecular models that maintain near-quantum accuracy. In this work, we extend the DeePMD-kit code to produce accurate molecular multipole moments in the bulk and near interfaces. The new method is validated by comparing the DNN moments with those generated by ab initio simulations. The moments are used to compute the electrostatic potential at the center of a molecular-sized hydrophobic cavity in water. The results show that the fields produced by the DNN models are in quantitative agreement with the AIMD-derived values. These efficient methods will open the door to more accurate solvation models for large solutes such as proteins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗