Image-charge effects on ion adsorption near aqueous interfaces

GCMe: Efficient Implementation of the Gaussian Core Model with Smeared Electrostatic Interactions for Molecular Dynamics Simulations of Soft Matter Systems

In recent years, molecular dynamics (MD) simulations have emerged as an essential tool for understanding the structure, dynamics, and phase behavior of charged soft matter systems. To explore phenomena across greater length and time scales in MD simulations, molecules are often coarse-grained for better computational performance. However, commonly used force fields represent particles as hard-core interaction centers with point charges, which often overemphasizes the packing effect and short-range electrostatics, especially in systems with bulky deformable organic molecules and systems with strong coarse-graining. This underscores the need for an efficient soft-core model to physically capture the effective interactions between coarse-grained particles. To this end, we implement a soft-core model uniting the Gaussian core model with smeared electrostatic interactions that is phenomenologically equivalent to recent theoretical models. We first parametrize it generically using water as the model solvent. Then, we benchmark its performance in the OpenMM toolkit for different boundary conditions to highlight a computational speedup of up to 34 × compared to commonly used force fields and existing implementations. Finally, we demonstrate its utility by investigating how boundary polarizability affects the adsorption behavior of a polyelectrolyte solution on perfectly conducting and nonmetal boundaries.

Accounting for the Quantum Capacitance of Graphite in Constant Potential Molecular Dynamics Simulations

Molecular dynamics (MD) simulations at a constant electric potential are an essential tool to study electrochemical processes, providing microscopic information on the structural, thermodynamic, and dynamical properties. Despite the numerous advances in the simulation of electrodes, they fail to accurately represent the electronic structure of materials such as graphite. In this work, a simple parameterization method that allows to tune the metallicity of the electrode based on a quantum chemistry calculation of the density of states (DOS) is introduced. As a first illustration, the interface between graphite electrodes and two different liquid electrolytes, an aqueous solution of NaCl and a pure ionic liquid, at different applied potentials are studied. It is shown that the simulations reproduce qualitatively the experimentally-measured capacitance; in particular, they yield a minimum of capacitance at the point of zero charge (PZC), which is due to the quantum capacitance (QC) contribution. An analysis of the structure of the adsorbed liquids allows to understand why the ionic liquid displays a lower capacitance despite its large ionic concentration. In addition to its relevance for the important class of carbonaceous electrodes, this method can be applied to any electrode materials (e.g. 2D materials, conducting polymers, etc), thus enabling molecular simulation studies of complex electrochemical devices in the future.

Excitation generated preferential binding sites for ethane on porous carbon-copper porphyrin sorbents: ethane/ethylene adsorptive separation improved by light

Energy-efficient separation of C2H6/C2H4 is a great challenge, for which adsorptive separation is very promising. C2H6-selective adsorption has big implications, while the design of C2H6-sorbents with ideal adsorption capability, particularly with the C2H6/C2H4-selectivity exceeded 2.0, is still challenging. Instead of the current strategies such as chemical modification or pore space modulation, we propose a new methodology for the design of C2H6-sorbents. With a Cu-TCPP [TCPP = 5,10,15,20-tetrakis(4-carboxyphenyl)porphyrin] framework dispersed onto a microporous carbon and a hierarchical-pore carbon, two composite sorbents are fabricated. The composite sorbents exhibit enhanced C2H6-selective adsorption capabilities with visible light, particularly the composite sorbent based on the hierarchical-pore carbon, whose C2H6 and C2H4 adsorption capacities (0 °C, 1 bar) are targetedly increased by 27% and only 1.8% with visible light, and therefore, an C2H6-selectivity (C2H6/C2H4 = 10/90, v/v) of 4.8 can be realized. With visible light, the adsorption force of the C2H6 molecule can be asymmetrically enhanced by the excitation enriched electron density over the adsorption sites formed via the close interaction between the Cu-TCPP and the carbon layer, whereas that of the C2H4 molecule is symmetrically altered and the forces cancelled each other out. This strategy may open up a new route for energy-efficient adsorptive separation of C2H6/C2H4 with light.

Impedance of nanocapacitors from molecular simulations to understand the dynamics of confined electrolytes

Nanoelectrochemical devices have become a promising candidate technology across various applications, including sensing and energy storage, and provide new platforms for studying fundamental properties of electrode/electrolyte interfaces. In this work, we employ constant-potential molecular dynamics simulations to investigate the impedance of gold-aqueous electrolyte nanocapacitors, exploiting a recently introduced fluctuation-dissipation relation. In particular, we relate the frequency-dependent impedance of these nanocapacitors to the complex conductivity of the bulk electrolyte in different regimes, and use this connection to design simple but accurate equivalent circuit models. We show that the electrode/electrolyte interfacial contribution is essentially capacitive and that the electrolyte response is bulk-like even when the interelectrode distance is only a few nanometers, provided that the latter is sufficiently large compared to the Debye screening length. We extensively compare our simulation results with spectroscopy experiments and predictions from analytical theories. In contrast to experiments, direct access in simulations to the ionic and solvent contributions to the polarization allows us to highlight their significant and persistent anticorrelation and to investigate the microscopic origin of the timescales observed in the impedance spectrum. This work opens avenues for the molecular interpretation of impedance measurements, and offers valuable contributions for future developments of accurate coarse-grained representations of confined electrolytes.

Generalized Helmholtz model describes capacitance profiles of ionic liquids and concentrated aqueous electrolytes

In this work, we propose and validate a generalization of the Helmholtz model that can account for both "bell-shaped" and "camel-shaped" differential capacitance profiles of concentrated electrolytes, the latter being characteristic of ionic liquids. The generalization is based on introducing voltage dependence of both the dielectric constant "ϵr(V)" and thickness "L(V)" of the inner Helmholtz layer, as validated by molecular dynamics (MD) simulations. We utilize MD simulations to study the capacitance profiles of three different electrochemical interfaces: (1) graphite/[BMIm+][BF4-] ionic liquid interface; (2) Au(100)/[BMIm+][BF4-] ionic liquid interface; (3) Au(100)/1M [Na+][Cl-] aqueous interface. We compute the voltage dependence of ϵr(V) and L(V) and demonstrate that the generalized Helmholtz model qualitatively describes both camel-shaped and bell-shaped differential capacitance profiles of ionic liquids and concentrated aqueous electrolytes (in lieu of specific ion adsorption). In particular, the camel-shaped capacitance profile that is characteristic of ionic liquid electrolytes arises simply from combination of the voltage-dependent trends of ϵr(V) and L(V). Furthermore, explicit analysis of the inner layer charge density for both concentrated aqueous and ionic liquid double layers reveal similarities, with these charge distributions typically exhibiting a dipolar region closest to the electrode followed by a monopolar peak at larger distances. It is appealing that a generalized Helmholtz model can provide a unified description of the inner layer structure and capacitance profile for seemingly disparate aqueous and ionic liquid electrolytes.

Electron Spin-Dependent Electrocatalysis for the Oxygen Reduction Reaction in a Chiro-Self-Assembled Iron Phthalocyanine Device

The chiral-induced spin selectivity effect (CISS) is a breakthrough phenomenon that has revolutionized the field of electrocatalysis. We report the first study on the electron spin-dependent electrocatalysis for the oxygen reduction reaction, ORR, using iron phthalocyanine, FePc, a well-known molecular catalyst for this reaction. The FePc complex belongs to the non-precious catalysts group, whose active site, FeN4, emulates catalytic centers of biocatalysts such as Cytochrome c. This study presents an experimental platform involving FePc self-assembled to a gold electrode surface using chiral peptides (L and D enantiomers), i.e., chiro-self-assembled FePc systems (CSAFePc). The chiral peptides behave as spin filters axial ligands of the FePc. One of the main findings is that the peptides' handedness and length in CSAFePc can optimize the kinetics and thermodynamic factors governing ORR. Moreover, the D-enantiomer promotes the highest electrocatalytic activity of FePc for ORR, shifting the onset potential up to 1.01 V vs. RHE in an alkaline medium, a potential close to the reversible potential of the O2 /H2 O couple. Therefore, this work has exciting implications for developing highly efficient and bioinspired catalysts, considering that, in biological organisms, biocatalysts that promote O2 reduction to water comprise L-enantiomers.

Water Electric Field Induced Modulation of the Wetting of Hexagonal Boron Nitride: Insights from Multiscale Modeling of Many-Body Polarization

Understanding the behavior of water contacting two-dimensional materials, such as hexagonal boron nitride (hBN), is important in practical applications, including seawater desalination and energy harvesting. Water, being a polar solvent, can strongly polarize the hBN surface via the electric fields that it generates. However, there is a lack of molecular-level understanding about the role of polarization effects at the hBN/water interface, including its effect on the wetting properties of water. In this study, we develop a theoretical framework that introduces an all-atomistic polarizable force field to accurately model the interactions of water molecules with hBN surfaces. The force field is then utilized to self-consistently describe the water-induced polarization of hBN using the classical Drude oscillator model, including predicting the hBN-water binding energies which are found to be in excellent agreement with diffusion Monte Carlo (DMC) predictions. By carrying out molecular dynamics (MD) simulations, we demonstrate that the polarizable force field yields a water contact angle on multilayered hBN which is in close agreement with the recent experimentally reported values. Conversely, an implicit modeling of the hBN-water polarization energy utilizing a Lennard-Jones (LJ) potential, a commonly utilized approximation in previous MD simulation studies, leads to a considerably lower water contact angle. This difference in the predicted contact angles is attributed to the significant energy-entropy compensation resulting from the incorporation of polarization effects at the hBN-water interface. Our work highlights the importance of self-consistently modeling the hBN-water polarization energy and offers insights into the wetting-related interfacial phenomena of water on polarizable materials.

Multiscale Modeling of Aqueous Electric Double Layers

From the stability of colloidal suspensions to the charging of electrodes, electric double layers play a pivotal role in aqueous systems. The interactions between interfaces, water molecules, ions and other solutes making up the electrical double layer span length scales from Ångströms to micrometers and are notoriously complex. Therefore, explaining experimental observations in terms of the double layer's molecular structure has been a long-standing challenge in physical chemistry, yet recent advances in simulations techniques and computational power have led to tremendous progress. In particular, the past decades have seen the development of a multiscale theoretical framework based on the combination of quantum density functional theory, force-field based simulations and continuum theory. In this Review, we discuss these theoretical developments and make quantitative comparisons to experimental results from, among other techniques, sum-frequency generation, atomic-force microscopy, and electrokinetics. Starting from the vapor/water interface, we treat a range of qualitatively different types of surfaces, varying from soft to solid, from hydrophilic to hydrophobic, and from charged to uncharged.

Regulating the electrical double layer to prevent water electrolysis for wet ionic liquids with cheap salts

Hydrophobic ionic liquids (ILs), broadly utilized as electrolytes, face limitations in practical applications due to their hygroscopicity, which narrows their electrochemical windows via water electrolysis. Herein, we scrutinized the impact of incorporating cheap salts on the electrochemical stability of wet hydrophobic ILs. We observed that alkali ions effectively manipulate the solvation structure of water and regulate the electrical double layer (EDL) structure by subtly adjusting the free energy distribution of water in wet ILs. Specifically, alkali ions significantly disrupted the hydrogen bond network, reducing free water, strengthening the O-H bond, and lowering water activity in bulk electrolytes. This effect was particularly pronounced in EDL regions, where most water molecules were repelled from both the cathode and anode with the disappearance of the H-bond network connectivity along the EDL. The residual interfacial water underwent reorientation, inhibiting water electrolysis and thus enhancing the electrochemical window of wet hydrophobic ILs. This theoretical proposition was confirmed by cyclic voltammetry measurements, demonstrating a 45% enhancement in the electrochemical windows for salt-in-wet ILs, approximating the dry one. This work offers feasible strategies for tuning the EDL and managing interfacial water activity, expanding the comprehension of interface engineering for advanced electrochemical systems.

Surface polarization enhances ionic transport and correlations in electrolyte solutions nanoconfined by conductors

Layered materials that perform mixed electron and ion transport are promising for energy harvesting, water desalination, and bioinspired functionalities. These functionalities depend on the interaction between ionic and electronic charges on the surface of materials. Here we investigate ion transport by an external electric field in an electrolyte solution confined in slit-like channels formed by two surfaces separated by distances that fit only a few water layers. We study different electrolyte solutions containing monovalent, divalent, and trivalent cations, and we consider walls made of non-polarizable surfaces and conductors. We show that considering the surface polarization of the confining surfaces can result in a significant increase in ionic conduction. The ionic conductivity is increased because the conductors' screening of electrostatic interactions enhances ionic correlations, leading to faster collective transport within the slit. While important, the change in water's dielectric constant in confinement is not enough to explain the enhancement of ion transport in polarizable slit-like channels.

Phase transitions of ionic fluids in nanoporous electrodes

In this work, we utilise grand canonical Metropolis Monte Carlo simulations, to establish pore-induced freezing of restricted primitive model fluids. A planar pore model is utilised, with walls that are initially neutral, and either non-conducting or perfectly conducting. The phase of the confined electrolyte (solid/fluid) displays an oscillatory dependence on surface separation, in narrow pores. Conditions are chosen so that the bulk is composed of a stable fluid electrolyte. The tendency for the electrolyte to freeze in narrow pores is somewhat stronger in systems with non-conducting walls. We also demonstrate that an applied potential will, above a threshold value, melt a frozen electrolyte. In these cases, the capacitance, as measured by the average surface charge density divided by the applied potential, will be almost vanishing if the applied potential is below this threshold value. We do not see any evidence for a "superionic fluid", which has been hypothesised to generate a strong capacitance in narrow pores, due to an efficient screening of like-charge repulsions by image charges.

Ion and water adsorption to graphene and graphene oxide surfaces

Graphene and graphene oxide (GO) are two particularly promising nanomaterials for a range of applications including energy storage, catalysis, and separations. Understanding the nanoscale interactions between ions and water near graphene and GO surfaces is critical for advancing our fundamental knowledge of these systems and downstream application success. This minireview highlights the necessity of using surface-specific experimental probes and computational techniques to fully characterize these interfaces, including the nanomaterial, surrounding water, and any adsorbed ions, if present. Key experimental and simulation studies considering water and ion structures near both graphene and GO are discussed. The major findings are: water forms 1-3 hydration layers near graphene; ions adsorb electrostatically to graphene under an applied potential; the chemical and physical properties of GO vary considerably depending on the synthesis route; and these variations influence water and ion adsorption to GO. Lastly, we offer outlooks and perspectives for these research areas.

Synergistic Effect of Salt and Anionic Surfactants on Interfacial Tension Reduction: Insights from Molecular Dynamics Simulations

Surfactants are commonly utilized in chemical flooding processes alongside salt to effectively decrease interfacial tension (IFT). However, the underlying microscopic mechanism for the synergistic effect of salt and surfactants on oil displacement remains ambiguous. Herein, the structure and properties of the interface between water and n-dodecane are studied by means of molecular dynamics simulations, considering three types of anionic surfactants and two types of salts. As the salt concentration (ρsalt) increases, the IFT first decreases to a minimum value, followed by a subsequent increase to higher values. The salt ions reduce the IFT only at low ρsalt due to the salt screening effect and ion bridging effect, both of which contribute to a decrease in the nearest head-to-head distance of surfactants. By incorporating salt doping, the IFTs can be reduced by at most 5%. Notably, the IFTs of different surfactants are mainly determined by the hydrogen bond interactions between oxygen atoms in the headgroup and water molecules. The presence of a greater number of oxygen atoms corresponds to lower IFT values. Specifically, for alkyl ethoxylate sulfate, the ethoxy groups play a crucial role in reducing the IFTs. This study provides valuable insights into formulating anionic surfactants that are applicable to oil recovery processes in petroleum reservoirs using saline water.

Molecular dynamics simulations of electrified interfaces including the metal polarisation

Understanding electrified interfaces requires an accurate description of the electric double layer which also takes into account the metal polarisation. Here we present a simple approach to the molecular dynamics simulation of electrified interfaces which combines fixed charges and a core-shell model for the description of the polarisable electron density on the metal electrode. The approach has been applied to the Au(111) surface in contact with a NaCl aqueous electrolyte solution in order to calculate the differential capacitance and to gain a detailed picture of the charging mechanism. Metal polarisation enhances the interfacial capacitance with a difference between the cathode and anode. In particular, we find that the influence of the metal polarisation on the electric double layer depends on the ion's solvation shell structure and, for the investigated interface, is more important at the cathode, where it modifies the sodium ion distribution.

PiNNwall: Heterogeneous Electrode Models from Integrating Machine Learning and Atomistic Simulation

Electrochemical energy storage always involves the capacitive process. The prevailing electrode model used in the molecular simulation of polarizable electrode-electrolyte systems is the Siepmann-Sprik model developed for perfect metal electrodes. This model has been recently extended to study the metallicity in the electrode by including the Thomas-Fermi screening length. Nevertheless, a further extension to heterogeneous electrode models requires introducing chemical specificity, which does not have any analytical recipes. Here, we address this challenge by integrating the atomistic machine learning code (PiNN) for generating the base charge and response kernel and the classical molecular dynamics code (MetalWalls) dedicated to the modeling of electrochemical systems, and this leads to the development of the PiNNwall interface. Apart from the cases of chemically doped graphene and graphene oxide electrodes as shown in this study, the PiNNwall interface also allows us to probe polarized oxide surfaces in which both the proton charge and the electronic charge can coexist. Therefore, this work opens the door for modeling heterogeneous and complex electrode materials often found in energy storage systems.

Fluids and Electrolytes under Confinement in Single-Digit Nanopores

Confined fluids and electrolyte solutions in nanopores exhibit rich and surprising physics and chemistry that impact the mass transport and energy efficiency in many important natural systems and industrial applications. Existing theories often fail to predict the exotic effects observed in the narrowest of such pores, called single-digit nanopores (SDNs), which have diameters or conduit widths of less than 10 nm, and have only recently become accessible for experimental measurements. What SDNs reveal has been surprising, including a rapidly increasing number of examples such as extraordinarily fast water transport, distorted fluid-phase boundaries, strong ion-correlation and quantum effects, and dielectric anomalies that are not observed in larger pores. Exploiting these effects presents myriad opportunities in both basic and applied research that stand to impact a host of new technologies at the water-energy nexus, from new membranes for precise separations and water purification to new gas permeable materials for water electrolyzers and energy-storage devices. SDNs also present unique opportunities to achieve ultrasensitive and selective chemical sensing at the single-ion and single-molecule limit. In this review article, we summarize the progress on nanofluidics of SDNs, with a focus on the confinement effects that arise in these extremely narrow nanopores. The recent development of precision model systems, transformative experimental tools, and multiscale theories that have played enabling roles in advancing this frontier are reviewed. We also identify new knowledge gaps in our understanding of nanofluidic transport and provide an outlook for the future challenges and opportunities at this rapidly advancing frontier.

Coupled Interactions at the Ionic Graphene-Water Interface

We compute ionic free energy adsorption profiles at an aqueous graphene interface by developing a self-consistent approach. To do so, we design a microscopic model for water and put the liquid on an equal footing with the graphene described by its electronic band structure. By evaluating progressively the electronic and dipolar coupled electrostatic interactions, we show that the coupling level including mutual graphene and water screening permits one to recover remarkably the precision of extensive quantum simulations. We further derive the potential of mean force evolution of several alkali cations.

Predictive Molecular Models for Charged Materials Systems: From Energy Materials to Biomacromolecules

Electrostatic interactions play a dominant role in charged materials systems. Understanding the complex correlation between macroscopic properties with microscopic structures is of critical importance to develop rational design strategies for advanced materials. But the complexity of this challenging task is augmented by interfaces present in the charged materials systems, such as electrode-electrolyte interfaces or biological membranes. Over the last decades, predictive molecular simulations that are founded in fundamental physics and optimized for charged interfacial systems have proven their value in providing molecular understanding of physicochemical properties and functional mechanisms for diverse materials. Novel design strategies utilizing predictive models have been suggested as promising route for the rational design of materials with tailored properties. Here, an overview of recent advances in the understanding of charged interfacial systems aided by predictive molecular simulations is presented. Focusing on three types of charged interfaces found in energy materials and biomacromolecules, how the molecular models characterize ion structure, charge transport, morphology relation to the environment, and the thermodynamics/kinetics of molecular binding at the interfaces is discussed. The critical analysis brings two prominent field of energy materials and biological science under common perspective, to stimulate crossover in both research field that have been largely separated.

Measuring anion binding at biomembrane interfaces

The quantification of anion binding by molecular receptors within lipid bilayers remains challenging. Here we measure anion binding in lipid bilayers by creating a fluorescent macrocycle featuring a strong sulfate affinity. We find the determinants of anion binding in lipid bilayers to be different from those expected that govern anion binding in solution. Charge-dense anions H2PO4- and Cl- that prevail in dimethyl sulfoxide fail to bind to the macrocycle in lipids. In stark contrast, ClO4- and I- that hardly bind in dimethyl sulfoxide show surprisingly significant affinities for the macrocycle in lipids. We reveal a lipid bilayer anion binding principle that depends on anion polarisability and bilayer penetration depth of complexes leading to unexpected advantages of charge-diffuse anions. These insights enhance our understanding of how biological systems select anions and guide the design of functional molecular systems operating at biomembrane interfaces.

Microscopic Simulations of Electrochemical Double-Layer Capacitors

Electrochemical double-layer capacitors (EDLCs) are devices allowing the storage or production of electricity. They function through the adsorption of ions from an electrolyte on high-surface-area electrodes and are characterized by short charging/discharging times and long cycle-life compared to batteries. Microscopic simulations are now widely used to characterize the structural, dynamical, and adsorption properties of these devices, complementing electrochemical experiments and in situ spectroscopic analyses. In this review, we discuss the main families of simulation methods that have been developed and their application to the main family of EDLCs, which include nanoporous carbon electrodes. We focus on the adsorption of organic ions for electricity storage applications as well as aqueous systems in the context of blue energy harvesting and desalination. We finally provide perspectives for further improvement of the predictive power of simulations, in particular for future devices with complex electrode compositions.

Understanding the Electric Double-Layer Structure, Capacitance, and Charging Dynamics

Significant progress has been made in recent years in theoretical modeling of the electric double layer (EDL), a key concept in electrochemistry important for energy storage, electrocatalysis, and multitudes of other technological applications. However, major challenges remain in understanding the microscopic details of the electrochemical interface and charging mechanisms under realistic conditions. This review delves into theoretical methods to describe the equilibrium and dynamic responses of the EDL structure and capacitance for electrochemical systems commonly deployed for capacitive energy storage. Special emphasis is given to recent advances that intend to capture the nonclassical EDL behavior such as oscillatory ion distributions, polarization of nonmetallic electrodes, charge transfer, and various forms of phase transitions in the micropores of electrodes interfacing with an organic electrolyte or ionic liquid. This comprehensive analysis highlights theoretical insights into predictable relationships between materials characteristics and electrochemical performance and offers a perspective on opportunities for further development toward rational design and optimization of electrochemical systems.

A coarse-grained model of room-temperature ionic liquids between metal electrodes: a molecular dynamics study

Recent mean-field theories predict that room-temperature ionic liquid (RTIL) electric double-layer capacitors (EDLCs) undergo a spontaneous surface charge separation (SSCS) with no applied potential. In this study, we construct a coarse-grained molecular model that corresponds to the mean-field models to directly simulate the behavior of RTILs without invoking mean-field approximations. In addition to observing the SSCS transition, we highlight the importance of the image charge interactions and explore the enhanced in-plane ordering on the electrodes, two effects not accounted for by the mean-field theories. By calculating and comparing the differential capacitance for RTILs confined between perfectly conducting and non-metal electrodes, we show that the image charge interactions drastically improve the energy storage properties of RTIL EDLCs.

Interactions between conducting surfaces in salt solutions

In this work, we simulate interactions between two perfectly conducting surfaces, immersed in a salt solution. We demonstrate that these forces are quantitatively different from those between (equally charged) non-conducting surfaces. There is, for instance, a significant repulsion between net neutral surfaces. On the other hand, there are also qualitative similarities, with behaviours found with non-conducting surfaces. For instance, there is a non-monotonic dependence of the free energy barrier height, on the salt concentration, and the minimum essentially coincides with a flat profile of the apparent surface charge density (i.e. the effective net surface charge density, some distance away from the surface, when accounting for ion neutralization), outside the so-called Stern layer. These conditions can be described as "perfect surface charge neutralization". Despite observed quantitative differences, we demonstrate that it might be possible to mimic a dispersion containing charged colloidal metal particles by a simpler model system with charged non-conducting particles, using modified particle-ion interactions.

Microscopic origin of the effect of substrate metallicity on interfacial free energies

We investigate the effect of the metallic character of solid substrates on solid-liquid interfacial thermodynamics using molecular simulations. Building on the recent development of a semiclassical Thomas-Fermi model to tune the metallicity in classical molecular dynamics simulations, we introduce a thermodynamic integration framework to compute the evolution of the interfacial free energy as a function of the Thomas-Fermi screening length. We validate this approach against analytical results for empty capacitors and by comparing the predictions in the presence of an electrolyte with values determined from the contact angle of droplets on the surface. The general expression derived in this work highlights the role of the charge distribution within the metal. We further propose a simple model to interpret the evolution of the interfacial free energy with voltage and Thomas-Fermi length, which allows us to identify the charge correlations within the metal as the microscopic origin of the evolution of the interfacial free energy with the metallic character of the substrate. This methodology opens the door to the molecular-scale study of the effect of the metallic character of the substrate on confinement-induced transitions in ionic systems, as reported in recent atomic force microscopy and surface force apparatus experiments.

A molecular perspective on induced charges on a metallic surface

Understanding the response of the surface of metallic solids to external electric field sources is crucial to characterize electrode-electrolyte interfaces. Continuum electrostatics offer a simple description of the induced charge density at the electrode surface. However, such a simple description does not take into account features related to the atomic structure of the solid and to the molecular nature of the solvent and of the dissolved ions. In order to illustrate such effects and assess the ability of continuum electrostatics to describe the induced charge distribution, we investigate the behavior of a gold electrode interacting with sodium or chloride ions fixed at various positions, in a vacuum or in water, using all-atom constant-potential classical molecular dynamics simulations. Our analysis highlights important similarities between the two approaches, especially under vacuum conditions and when the ion is sufficiently far from the surface, as well as some limitations of the continuum description, namely, neglecting the charges induced by the adsorbed solvent molecules and the screening effect of the solvent when the ion is close to the surface. While the detailed features of the charge distribution are system-specific, we expect some of our generic conclusions on the induced charge density to hold for other ions, solvents, and electrode surfaces. Beyond this particular case, the present study also illustrates the relevance of such molecular simulations to serve as a reference for the design of improved implicit solvent models of electrode-electrolyte interfaces.

Conformation and Ionization Behavior of Charge-Regulating Polyelectrolyte Brushes in a Poor Solvent

Understanding the structural response of weak polyelectrolyte brushes upon external stimuli is crucial for their applications ranging from modifying surface properties to the development of smart and intelligent materials. In this work, coarse-grained molecular dynamics simulations were carried out to investigate the conformation and ionization behavior of charge-regulating polyelectrolyte brushes under poor solvent conditions, using an implicit solvent model. The results show that, while the thickness of a sparse polyelectrolyte brush shows a similar behavior to that of a single chain, namely, a monotonic change as a function of solvent quality (modeled by an effective segment-segment attraction strength parameter) and solution pH, a dense polyelectrolyte brush exhibits more complex behavior. An unexpected reexpansion is observed when the effective segment-segment attraction strength is further increased, especially in the case of a high pH. In the latter case, strong attraction in polymer segments promotes the formation of large, interchain, cylindrical aggregates, leading to an increase in brush thickness.

Effect of the metallicity on the capacitance of gold-aqueous sodium chloride interfaces

Electrochemistry experiments have established that the capacitance of electrode-electrolyte interfaces is much larger for good metals, such as gold and platinum, than for carbon-based materials. Despite the development of elaborate electrode interaction potentials, to date molecular dynamics simulations are not able to capture this effect. Here, we show that changing the width of the Gaussian charge distribution used to represent the atomic charges in gold is an effective way to tune its metallicity. Larger Gaussian widths lead to a capacitance of aqueous solutions (pure water and 1 M NaCl) in good agreement with recent ab initio molecular dynamics results. For pure water, the increase in the capacitance is not accompanied by structural changes, while in the presence of salt, the Na⁺ cations tend to adsorb significantly on the surface. For a strongly metallic gold electrode, these ions can even form inner sphere complexes on hollow sites of the surface.

Effects of Confinement and Ion Adsorption in Ionic Liquid Supercapacitors with Nanoporous Electrodes

We investigate the effects of pore size and ion adsorption on the room-temperature ionic liquid capacitor with nanoporous electrodes, with a focus on optimizing the capacitance and energy storage. Using a recently developed modified BSK model accounting for both ion correlations and nonelectrostatic interactions, we find that ion crowding proximate to the electrode surface induced by the spontaneous charge separation due to strong ion correlations is responsible for the anomalous increase in the capacitance with decreasing pore sizes observed in experiments. Reducing the strength of ion correlations increases the capacitance and suppresses the anomalous size dependence. For a given pore size, the capacitance peak diverges when the ion correlation strength α reaches a critical value, α_sc,L. The capacitance peak shifts to smaller pore size as α decreases because of rapid decrease of α_sc,L with decreasing pore size. Asymmetric preferential ion adsorption is shown to lead to significantly enhanced energy storage close to the transition point for any pore sizes. For a given correlation strength, the energy storage is optimal at a pore size where α = α_sc,L.