SEARCH · Engineering Papers
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.
We establish a solid/water interface model which treats ionic defects in solids and ions in water uniformly, allowing one to study the space charge layer and the electric double layer profiles under doping and pH effects.
In photoelectrochemical cells, promising devices for directly converting solar energy into storable chemical fuels, the spatial variation of the electrostatic potential across the semiconductor–electrolyte junction is the key parameter that determines the cell performance. In principle, electric field induced second harmonic generation (EFISH) provides a contactless in situ spectroscopic tool to measure the spatial variation of electrostatic potential. However, the total second harmonic generation (SHG) signal contains the contributions of the EFISH signals of semiconductor space charge layer and the electric double layer, in addition to the SHG signal of the electrode surface. The interference of these complex quantities hinders their analysis. In this work, to understand and deconvolute their contributions to the total SHG signals, bias-dependent SHG measurements are performed on the rutile TiO2(100)–electrolyte junction as a function of light polarization and crystal azimuthal angle (angle of the incident plane relative to the crystal [001] axis). A quadratic response between SHG intensity and the applied potential is observed in both the accumulation and depletion regions of TiO2. The relative phase difference and amplitude ratio are extracted at selected azimuthal angles and light polarizations. At 0° azimuthal angle and s-in–p-out polarization, the SHG intensity minimum has the best match with the TiO2 flatband potential due to the orthogonal relative phase difference between bias-dependent and bias-independent SHG terms. We further measure the pH-dependent flatband potential and probe the photovoltage under open circuit conditions using the EFISH technique, demonstrating the capability of this contactless method for measuring electrostatic potential at semiconductor–electrolyte junctions.
Electric double layer capacitors (EDLCs) are prominent energy storage systems that constitute the foundation of more reliable and sustainable energy infrastructures. Modern EDLCs often incorporate ionic liquids (ILs) as a key component in their electrolytes, leveraging the high electrochemical stability of ILs to enhance device performance. The performance and functionality of these capacitors also hinge on the interfacial behavior of ILs at electrode surfaces, which remain insufficiently understood. Here, we performed synchrotron infrared nanospectroscopy (SINS) in combination with density functional theory (DFT) calculations to investigate the electric double layers (EDLs) of three imidazolium-based ILs in a custom-designed graphene liquid cell. This approach revealed new insights into the dynamics of IL EDLs and the underpinning factors originating from the IL structures. Variations in anion size and structure were found to tune the ILs’ ability to form EDLs with compact and closely-correlated ion arrangements, which are critical for enhancing their capacitive performance. These findings highlight the intricate interactions between ions governed by their structures and charge behaviors, which underscores the opportunities for targeted design of IL-based electrolytes to optimize EDLC functionality.
Electric double layers are ubiquitous in science and engineering and are of current interest, owing to their applications in the stabilization of colloidal suspensions and as supercapacitors. While the structure and properties of electric double layers in electrolyte solutions near a charged surface are well characterized, there are subtleties in calculating thermodynamic properties from the free energy of a system with charged surfaces. These subtleties arise from the difference in the free energy between systems with constant surface charge and constant surface potential. In this work, we present a systematic, pedagogical framework to properly account for the different specifications on charged bodies in electrolyte solutions. Our approach is fully variational—that is, all free energies, boundary conditions, relevant electrostatic equations, and thermodynamic quantities are systematically derived using variational principles of thermodynamics. We illustrate our approach by considering a simple electrolyte solution between two charged surfaces using the Poisson–Boltzmann theory. Our results highlight the importance of using the proper thermodynamic potential and provide a general framework for calculating thermodynamic properties of electrolyte solutions near charged surfaces. Specifically, we present the calculation of the pressure and the surface tension between two charged surfaces for different boundary conditions, including mixed boundary conditions.
Not provided.
Significance For centuries the double layer at the solid/electrolyte interface has been a central concept in electrochemistry. Today, it is still crucial for virtually all renewable energy storage and conversion technologies. Here, the double-layer formation is probed by THz spectroscopy with ultrabright synchrotron light as a source. Our results capture the molecular details of double-layer formation at positively/negatively charged Au electrodes for an NaCl electrolyte. We reveal a contrasting response applying positive versus negative bias, which is dictated by the interfacial water network and rationalized by accompanying molecular dynamics simulations and electronic-structure calculations. While Na + is directly attracted toward the negatively charged electrode, stripping of the Cl − hydration shell is observed only at larger potential values.
For many slow solar-fuel-forming reactions, the accumulation of photogenerated minority carriers on the photoelectrode surface leads to light-induced band edge unpinning, affecting the junction properties by decreasing band bending in the semiconductor space charge layer and increasing the driving force of surface reactions in the electric double layer. In this study, we demonstrate a contactless operando electric field-induced second harmonic generation (EFISH) method for measuring the band bending change (δΔΦ SCR L ) on photoelectrodes upon photoexcitation. For ndoped rutile TiO 2 water oxidation photoanodes at pH 7, δΔΦ SCR L increases at more positive potentials or higher illumination power density until it reaches saturation values. We show that under fast mass transport conditions, δΔΦ SCR L is exclusively attributed to the accumulated charged rate-determining species that can be regarded as temporary surface states, and the relationship between the photocurrent and δΔΦ SCR L can be well modeled by assuming that hole trap states function as the reaction center. Kinetic isotope experiments identify proton-coupled electron transfer as the rate-determining step and suggest a possible chemical nature of the key intermediate. We demonstrate that light-induced band edge unpinning is a beneficial feature under high illumination conditions for oxygen evolution reaction on TiO 2 because it maintains the photon-to-current conversion efficiency by enhancing the surface reaction driving force, shedding light on the actual device application.
Complex electrical conductivity is increasingly used to monitor subsurface processes associated with microbial activities because microbial cells mostly have surface charges and thus electrical double layers. Although highly saline environments are frequently encountered in coastal and marine sediments, there are limited data available on the complex conductivity associated with microbial activities under a high-salinity condition. Therefore, we have developed the spectral responses of complex conductivity of sand associated with bacterial growth and biofilm formation under a highly saline condition of approximately 1% salinity and approximately 2 S/m pore water conductivity with an emphasis on relaxation behavior. A column test is performed, in which the model bacteria Shewanella oneidensis MR-1 are stimulated for cell growth and biofilm formation in a sand pack, whereas the complex conductivity is monitored from 0.01 Hz to 10 kHz. The test results indicate that the real conductivity increases in the early stage due to the microbial metabolites and the increased surface conduction with cell growth but soon begin to decrease because of the reduction of charge passages due to bioclogging. However, the imaginary conductivity significantly increases with time, and clear bell-shaped relaxation behaviors are observed with the peak frequency of 0.1–1 Hz, associated with the double-layer polarization of cells and electrically conductive pili and biofilms. The Cole-Cole relaxation model appears to capture such relaxation behaviors well, and the modeling results indicate gradual increases in normalized chargeability and decreases in relaxation time during bacterial growth and biofilm formation in the highly saline condition. Comparison with previous literature confirms that the high-salinity condition further increases the normalized chargeability, whereas it suppresses the phase shift and thus the imaginary conductivity. Our results suggest that the complex conductivity can effectively capture microbial biomass formation in sands under a highly saline condition.
We report the results of constant-potential molecular dynamics simulations of the double layer interface between molten 2LiF–BeF 2 (FLiBe) and 23LiF–6NaF–21KF (FLiNaK) fluoride mixtures and idealized solid electrodes. Employing methods similar to those used in studies of chloride double layers, we compute the structure and differential capacitance of molten fluoride electric double layers as a function of applied voltage. The role of molten salt structure is probed through comparisons between FLiBe and FLiNaK, which serve as models for strong and weak associate-forming salts, respectively. In FLiBe, screening involves changes in Be–F–Be angles and alignment of the oligomers parallel to the electrode, while in FLiNaK, the electric field is screened mainly by rearrangement of individual ions, predominantly the polarizable potassium cation.
Oxide-water interfaces underpin ion separation, catalysis, and electrochemical energy technologies, where the electrical double layer (EDL) controls adsorption, transport, and reactivity. Yet, the molecular-scale link between pH-dependent surface protonation, hydration-layer structure, and counter-ion adsorption remains poorly defined. Here, we combine in situ crystal truncation rod (CTR) and resonant anomalous X-ray reflectivity (RAXR) with streaming potential measurements and ab initio molecular dynamics (AIMD) simulations to resolve the chemical and structural evolution of the EDL at the single-crystal alumina (012)-water interface in 10 mM Rb+ over pH 3-12. CTR measurements reveal two distinct adsorbed water layers at ~2.2 and ~3.5 Å above the surface that each shift toward the substrate at transition pHs near 6.5 and 10.6, respectively, directly reflecting changes in primary hydration layer structure in response to the deprotonation of bridging and terminal aluminol groups. RAXR shows a 10-fold increase in Rb+ coverage and a decrease in mean adsorption height from ~3.5 to ~2.7 Å with increasing pH, indicating enhanced counter-ion binding accompanied by Stern layer contraction. Streaming potential measurements demonstrate that the zeta potential, i.e., potential at the hydrodynamic shear plane, is positive at pH 3 and becomes negative at pH ≥3.5. This negative charge magnitude increases with increasing pH, consistent with progressive surface deprotonation at higher pH. AIMD identifies inner- and outer-sphere Rb+ complexes whose adsorption heights and coordination geometries depend sensitively on the protonation state of surface oxygens, providing atomistic support for the experimentally inferred trends. These measurements establish two discrete, site-specific pH transitions in hydration-layer structure that track aluminol (de)protonation and quantitatively link them to a pH-driven contraction of the Stern layer (increasing Rb+ coverage and decreasing adsorption height). This provides a direct structural basis for connecting surface acid-base chemistry to ion binding distances at an oxide-water interface.
The solvation environments of Li + in conventional nonaqueous battery electrolytes, such as LiPF 6 in mixtures of ethylene carbaronate (EC) and ethyl methyl carbonate (EMC), are often used to rationalize transport properties and solid electrolyte interphase (SEI) formation. Solvation environments in the compact electrical double layer (EDL) next to the electrode, also known as the Helmholtz layer, determine (partially) what species can react to form the SEI, with bulk solvation environments often being used as a proxy. Here, we develop and test a theory of cation solvation in the Helmholtz layer of nonaqueous Li-ion battery electrolytes. First, we validate the theory against bulk and diffuse EDL atomistic molecular dynamics (MD) simulations of LiPF 6 EC/EMC mixtures as a function of surface charge, where we find the theory can qualitatively capture the solvation environments. Next, we turn to the Helmholtz layer, where we find the main effect of the solvation structures next to the electrode is an apparent reduction in the number of binding sites between Li + and the solvents, again where we find reasonable agreement with our developed theory. Finally, by solving a simplified version of the theory, we find that the probability of Li + binding to each solvent remains equal to the bulk probability, suggesting that the bulk solvation environments are a reasonable place to start when understanding battery electrolytes. Our developed formalism can be parametrized from bulk MD simulations and used to predict the solvation environments in the Helmholtz layer through reducing the number of available coordination sites, which can be used to determine what could react and form the SEI.
Energy generation through nanofluidics is a topic of great nanotechnological relevance. In this work, we conduct all-atom molecular dynamics (MD) simulations of the transport of water and ions in a pressure-driven flow in nanochannels grafted with charged polyelectrolyte (PE) brushes and discover the possibility of simultaneous electrokinetic energy generation and flow enhancement (henceforth denoted as the electroslippage effect). Such PE-brush-functionalized nanochannels have been recently shown to demonstrate an overscreening (OS) effect (characterized by the presence of a greater number of screening counterions within the PE brush layer than needed to screen the PE brush charges), a consequent presence of excess coions within the PE brush-free bulk, and a coion-driven electroosmotic (EOS) transport in the presence of small to moderate applied axial electric fields. Here, however, we find that the streaming current, which represents the current generated by the flow-driven downstream advection of the charge imbalance present within the electric double layer (EDL) that screen the PE brush charges, is governed by the migration of the counterions. This stems from the fact that the highest contribution to the overall streaming current arises from the region near the PE brush-water interface (where there is an excess of counterions), while the brush-free bulk yields a hitherto unreported, but small, coion-dictated streaming current. This downstream advection of the charge imbalance (and the resultant counterion-driven streaming current) eventually leads to the development of an electric field (streaming electric field) in the direction that is opposite to the direction of the counterion-driven streaming current. The streaming current and the streaming electric field interact to generate the electrokinetic energy. Equally importantly, this streaming electric field induces an EOS transport, which becomes coion-driven, due to the presence of excess coions in the brush-free bulk. For the case of nanochannels grafted with negatively charged PE brushes, the streaming electric field will be in a direction that is opposite to the pressure-driven transport, and hence the coion (or anion) driven EOS flow, will be in the same direction as the pressure-driven transport. On the other hand, for the case of nanochannels grafted with positively charged PE brushes, the streaming electric field will be in the same direction as the pressure-driven flow, and hence the coion (or cation) driven EOS flow, will again be in the same direction as the pressure-driven flow. Therefore, whenever there occurs the presence of the OS and the resulting coion-driven EOS transport in PE brush grafted nanochannels, regardless of the sign of the charges of the PE brushes, this EOS transport will always aid the pressure-driven transport and will cause the most fascinating increase in the net volume flow rate across the nanochannel cross section, which is the electroslippage effect.
Knowledge Gap: The aggregation of clay minerals—layered silicate nanoparticles—strongly impacts fluid flow, solute migration, and solid mechanics in soils, sediments, and sedimentary rocks. Experimental and computational characterization of clay aggregation is inhibited by the delicate water-mediated nature of clay colloidal interactions and by the range of spatial scales involved, from 1 nm thick platelets to flocs with dimensions up to micrometers or more. Simulations: Using a new coarse-grained molecular dynamics (CGMD) approach, we predicted the microstructure, dynamics, and rheology of hydrated smectite (more precisely, montmorillonite) clay gels containing up to 2,000 clay platelets on length scales up to 0.1 μm. Further, simulations investigated the impact of simulation time, platelet diameters (6 to 25nm), and the ratio of Na to Ca exchangeable cations on the assembly of tactoids (i.e., stacks of parallel clay platelets) and larger aggregates (i.e., assemblages of tactoids). We analyzed structural features including tactoid size and size distribution, basal spacing, counterion distribution in the electrical double layer, clay association modes, and the rheological properties of smectite gels. Findings: Our results demonstrate new potential to characterize and understand clay aggregation in dilute suspensions and gels on a scale of thousands of particles with explicit representation of counterion clouds and with accuracy approaching that of all-atom molecular dynamics (MD) simulations. For example, our simulations predict the strong impact of Na/Ca ratio on clay tactoid formation and the shear-thinning rheology of clay gels.
Classical electric double layer (EDL) models have been widely used to describe ion distributions at charged solid-water interfaces in dilute electrolytes. However, the chemistry of EDLs remains poorly constrained at high ionic strength where ion-ion correlations control non-classical behavior such as overcharging, i. e., the accumulation of counter-ions in amounts exceeding the substrate's surface charge. Here, we provide direct experimental observations of correlated cation and anion distributions adsorbed at the muscovite (001)-aqueous electrolyte interface as a function of dissolved RbBr concentration ([RbBr]=0.01–5.8 M) using resonant anomalous X-ray reflectivity. Our results show alternating cation-anion layers in the EDL when [RbBr]≳100 mM, whose spatial extension (i. e., ~20 Å from the surface) far exceeds the dimension of the classical Stern layer. Comparison to RbCl and RbI electrolytes indicates that these behaviors are sensitive to the choice of co-ion. This new in-depth molecular-scale understanding of the EDL structure during transition from classical to non-classical regimes supports the development of realistic EDL models for technologies operating at high salinity such as water purification applications or modern electrochemical storage.
Significance Ion distributions at charged solid–water interfaces, referred to as the electrical double layer (EDL), are poorly understood at high ion concentrations, in part due to the lack of molecular-scale descriptions of the interactions between adsorbed hydrated ions. Here, direct visualization of the salinity-dependent evolution of EDL structure reveals molecular origins of nonclassical transformation of the EDL, in which charge overscreening and heterogeneous nucleation are driven by ion–ion correlations at the interfaces. This manifestation of the atomistic basis of nonclassical behaviors provides a much-needed understanding of the impact of ion cooperativity at charged interfaces for the development of predictive models for element transport in natural environments and advanced technologies for material growth and synthesis in saline environments.
Accurately predicting solid electrolyte interphase (SEI) formation requires explicitly resolving the electric double layer (EDL) structure, which deviates significantly from that of the bulk electrolyte. Although an established molecular dynamics (MD) and Density Functional Theory (DFT) framework can model SEI formation by evaluating reduction reactions of local clusters in the EDL, it suffers from a combinatorial computational bottleneck. To overcome this limitation, we introduce a machine-learning-accelerated simulation workflow (MD–DFT–ML), integrating a gradient-boosted regression model trained on EDL composition data to efficiently predict reduction potentials. We apply this framework to seven fluorinated electrolytes comprising fluorinated anions, a fluorinated ester solvent, two types of diluent (ion-solvating ester vs. non-solvating ether), and an FEC additive. The analysis shows that the EDL selectively accumulates cation-binding species; consequently, the non–cation-binding ether diluent rarely enters the EDL and makes minimal contributions to SEI formation. DFT calculations on statistically representative EDL clusters provide reduction potentials and fluorine-release pathways, while the ML model, which substantially reduces the DFT workload, predicts cluster reduction energies with a mean absolute error of 0.1 eV. The combined MD–DFT–ML approach also quantifies contributions from different sources to LiF formation in the SEI. This methodology establishes a generalizable route for multiscale modeling electrolyte and interphase design for next-generation electrochemical energy-storage systems.
Recent one-dimensional and two-dimensional numerical simulations of double layers (DLs) in the electric fields of the auroral plasma are reviewed, with reference to observational data. It is found that two-dimensional DLs driven by current sheets of finite thickness have different characteristics, depending on whether the layer thickness is less than or much greater than the ion gyroradius: When thickness is less than ion gyroradius, V-shaped DLs form with nearly equal parallel and perpendicular potential drops; when layer thickness is much greater than ion gyroradius the major parallel potential drop occurs outside the current sheet and the perpendicular electric fields are localized at the edges of the current sheet. It is shown that some features of the simulated fields, such as the amplitudes and scale lengths, are qualitatively similar to those observed in space.
The Earth's auroral zone contains dynamic processes occurring on scales from the length of an auroral zone field line which characterizes Alfven wave propagation to the scale of microscopic processes which occur over a few Debye lengths. These processes interact in a time-dependent fashion since the current carried by the Alfven waves can excite microscopic turbulence which can in turn provide dissipation of the Alfven wave energy. This review will first describe the dynamic aspects of auroral current structures with emphasis on consequences for models of microscopic turbulence. A number of models of microscopic turbulence will be introduced into a large-scale model of Alfven wave propagation to determine the effect of various models on the overall structure of auroral currents. In particular, the effects of a double layer electric field which scales with the plasma temperature and Debye length is compared with the effect of anomalous resistivity due to electrostatic ion cyclotron turbulence in which the electric field scales with the magnetic field strength. It is found that the double layer model is less diffusive than in the resistive model leading to the possibility of narrow, intense current structures.