Computer simulation of the structure of the electrochemical double layer

We analyze and compare the structure of the electrochemical double layer obtained from molecular dynamics simulations of concentrated aqueous NaCl and CsF solutions near a model electrode. The electrode is modeled as a corrugated external potential in conjunction with the image charge model. Calculations are performed for uncharged electrodes and for electrodes carrying positive or negative surface charges.     

Electrochemistry in nanometer-wide electrochemical cells

The electrochemical properties of an electrochemical cell defined by two concentric spherical electrodes, separated by a 1 to 20-nm-wide gap, and a freely diffusing electrochemically active molecule (e.g., ferrocene) have been investigated by coupling of Brownian dynamics simulations with long-range electron-transfer probability values. The simulation creates a trajectory of a single molecule and calculates the likelihood that the molecule undergoes a redox reaction during each time interval based on a probability-distance function derived from literature first-order kinetic data for a surface-bound ferrocene. Steady-state voltammograms for the single-molecule concentric spherical electrochemical cell are computed and are used to extract a heterogeneous electron-transfer rate for the freely diffusing molecule redox reaction. The Brownian dynamics simulations also indicate that long-range electron transfer, between the redox molecule and electrode, leads to nonsigmoidal-shaped i-E characteristics when the distance between electrodes approaches the characteristic redox tunneling decay length. The long-range electron transfer generates a "tunneling depletion layer" that results in a potential-dependent diffusion-limited current.     

Molecular dynamics simulations of electrolyte solutions at the (100) goethite surface

Molecular dynamics simulations of electrolyte solutions in contact with a neutral (100) goethite (alpha-FeOOH) surface were used to probe the structure of the mineral-water interface and gain insight into the adsorption properties of monovalent ions. Three electrolyte solutions were considered: NaCl, CsCl, and CsF. The electrolyte ions were chosen to cover a range of ionic sizes and affinities for the aqueous phase. The molecular dynamics simulations indicate the presence of a structured interfacial region resulting from the strong interaction of water with the mineral surface. The specific arrangement and preferred orientation of water that arise from this interaction create adsorption sites in the interfacial region, i.e., as far as 15 A away from the surface, and hence give rise to a strong correlation between the water and ion distributions. The structure of the hydrated ion, its effect on the water arrangement at the interface, and the strength of the ion-water bond are found to be key factors that determine the location and extent of ion adsorption at the interface. Additionally, in all simulations, we find a build up of positive charges near the surface due to cation adsorption, which is compensated by an accumulation of anions in the next few angstr?ms. This creates an excess of negative charges, which is in turn compensated by an excess of positive charges, and so on. As we modeled a neutral surface, the structure of the electrolyte distribution arises from the complex interplay of the interactions between the surface, water, and the electrolyte ions rather than from the need to neutralize a surface charge. In addition, our simulations indicate that the electrolyte distribution does not resemble that of a classical electrical double layer. Indeed, our calculations predict the presence of several condensed layers and oscillations in the net charge away from the surface.     

A molecular simulation study on the role of ion sizes and dielectric images in near-surface ion distribution far from the strong coupling limit

A series of Monte Carlo simulations of the planar electric double layers are carried out in the primitive model for two electrolyte mixtures next to a smooth and uniformly charged hard wall representing an ideal biological interface with low and moderate surface charge densities. The structural information of the double layers is applied to reveal charge inversion and overcharging through the addition of multivalent electrolyte at a certain physiological concentration. Various values for the radius of the ions are taken into account to capture the impact of short-range correlations. Meanwhile, the influence of image charges on ion distribution is analyzed, which stems from dielectric discontinuity between the interior and exterior of the membrane matrix. It is clearly shown that depending on the amount of foreign salt, the large size of charged species regardless of its polarity plays a positive role in promoting charge inversion. Moreover, our findings indicate that charge inversion do not signify the reversal of the electrophoretic mobility, in consistent with the recent theoretical predictions by Horno and co-workers [J. Colloid Interface Sci. 356, 325 (2011)]. In addition, the depletion effect triggered by repulsive image forces which are intertwined with the excluded volume correlations gives rise to an anomalous overcharging for low surface charged surface in the high concentrations of trivalent salt. Overall, the ion distribution in a double layer is exclusively governed by entropic and electrostatic contributions but with preferentially leading status for different magnitudes of surface charge.     

1,3-二氧戊环基LiCF3SO3电解液对锂硫电池正极材料单质硫的电化学性能影响

测试了二元和多元溶剂组分的1,3-二氧戊环基LiCF3SO3电解液的粘度、离子电导率和单质硫的溶解度. 研究结果表明, 由较强的给电子能力溶剂组成的低粘度电解液较容易提高单质硫的氧化还原反应活性和可逆性能, 有利于提高单质硫在2.10 V附近的低放电平台电位和放电比容量. DOL-DME LiCF3SO3电解液能够较好地改善单质硫电极的表面钝化层结构, 促进电活性物质离子扩散和降低界面电荷传递阻抗, 从而表现出很好的放电倍率特性. 在室温下充放电流密度分别为0.1和0.2 mA/cm2时, 单质硫的首次放电比容量为792 mA·h/g, 第29次放电比容量达到412 mA·h/g.     

Effects of image charges, interfacial charge discreteness, and surface roughness on the zeta potential of spherical electric double layers

We investigate the effects of image charges, interfacial charge discreteness, and surface roughness on spherical electric double layer structures in electrolyte solutions with divalent counterions in the setting of the primitive model. By using Monte Carlo simulations and the image charge method, the zeta potential profile and the integrated charge distribution function are computed for varying surface charge strengths and salt concentrations. Systematic comparisons were carried out between three distinct models for interfacial charges: (1) SURF1 with uniform surface charges, (2) SURF2 with discrete point charges on the interface, and (3) SURF3 with discrete interfacial charges and finite excluded volume. By comparing the integrated charge distribution function and the zeta potential profile, we argue that the potential at the distance of one ion diameter from the macroion surface is a suitable location to define the zeta potential. In SURF2 model, we find that image charge effects strongly enhance charge inversion for monovalent interfacial charges, and strongly suppress charge inversion for multivalent interfacial charges. For SURF3, the image charge effect becomes much smaller. Finally, with image charges in action, we find that excluded volumes (in SURF3) suppress charge inversion for monovalent interfacial charges and enhance charge inversion for multivalent interfacial charges. Overall, our results demonstrate that all these aspects, i.e., image charges, interfacial charge discreteness, their excluding volumes, have significant impacts on zeta potentials of electric double layers.     

The role of supporting electrolyte in heterogeneous electron transfer

The effects of supporting electrolyte on the kinetics of the elementary step of electron transfer are considered as unavoidable interplay of interfacial phenomena and ionic equilibria in solution. For the former, the problems to separate contributions of electrostatic electrode-reactant interactions and specific adsorption are addressed, and various aspects of the traditional Frumkin correction (“psi-prime effect”) are discussed. The construction of corrected Tafel plots is shown to be a procedure containing the internal contradiction resulting in an uncertainty. This uncertainty can be eliminated by combining the principles of traditional analysis of the “double layer” effects with physical theory instead of phenomenological approaches. Specific manifestations of parallel electron transfer to an ensemble of reacting species are presented in the context of “mean reactant charge in solution bulk.” The approach to account for non-spherical shape and inhomogeneous charge distribution in reacting species is considered in terms of “molecular psi-prime effect.” Finally, some comments are given on analogy of “double layer” effects at metal/solution interface and interfacial phenomena specific for more complex and highly relevant electrochemical systems.     

Quantum dynamics simulations of interfacial electron transfer in sensitized TiO2 semiconductors

Ab initio DFT molecular dynamics simulations are combined with quantum dynamics calculations of electronic relaxation to investigate the interfacial electron transfer in catechol/TiO(2)-anatase nanostructures under vacuum conditions. It is found that the primary process in the interfacial electron-transfer dynamics involves an ultrafast (tau(1) approximately 6 fs) electron-injection event that localizes the charge in the Ti(4+) surface ions next to the catechol adsorbate. The primary event is followed by charge delocalization (i.e., carrier diffusion) through the TiO(2)-anatase crystal, an anisotropic diffusional process that can be up to an order of magnitude slower along the [-101] direction than carrier relaxation along the [010] and [101] directions in the anatase crystal. It is shown that both the mechanism of electron injection and the time scales for interfacial electron transfer are quite sensitive to the symmetry of the electronic state initially populated in the adsorbate molecule. The results are particularly relevant to the understanding of surface charge separation in efficient mechanisms of molecular-based photovoltaic devices.     

Hierarchical multiscale simulation of electrokinetic transport in silica nanochannels at the point of zero charge

Effects of nanoscale confinement and partial charges that stem from quantum calculations are investigated in silica slit channels filled with 1 M KCl at the point of zero charge by using a hierarchical multiscale simulation methodology. Partial charges of both bulk and surface atoms from ab initio quantum calculations that take into account bond polarization and electronegativity are used in molecular dynamics (MD) simulations to obtain ion and water concentration profiles for channel widths of 1.1, 2.1, 2.75, and 4.1 nm. The interfacial electron density profiles of simulations matched well with that of recent X-ray reflectivity experiments. By simulating corresponding channels with no partial charges, it was observed that the partial charges affect the concentration profiles and transport properties such as diffusion coefficients and mobilities up to a distance of about 3 sigma(O)(-)(O) from the surface. Both in uncharged and partially charged cases, oscillations in concentration profiles of K(+) and Cl(-) ions give rise to an electro-osmotic flow in the presence of an external electric field, indicating the presence of an electric double layer at net zero surface charge, contrary to the expectations from classical continuum theory. I-V curves in a channel-bath system using ionic mobilities from MD simulations were significantly different for channels with and without partial charges for channel widths less than 4.1 nm.     

Dynamics of charge transfer states on metal surfaces: the competition between reactivity and quenching

The dynamics of excited states of adsorbates on surfaces caused by charge transfer is studied. Both negative and positive charge transfer processes are possible. In particular we are interested in positive charge transfer from a metal surface to molecular or atomic oxygen adsorbed on the surface. Once the negatively charged oxygen on the surface loses an electron it becomes chemically activated. The ability of this species to react depends on the quenching time or back transfer. The analysis of these processes is based on a set of diabatic potential energy surfaces each representing a different charged oxygen species. The dynamics is followed by solving the multichannel time-dependent Schr?dinger equation or Liouville von Neumann equation. Due to the nonadiabatic character of these reactions large isotope effects are predicted.     

Time-domain ab initio study of charge relaxation and recombination in dye-sensitized TiO2

In order to investigate the electron dynamics at the alizarin/I2-/TiO2 interface this study uses a novel state-of-the-art quantum-classical approach that combines time-dependent density functional theory with surface hopping in the Kohn-Sham basis. Representing the dye-sensitized semiconductor Gr?tzel cell with the I-/I3- mediator, the system addresses the problems of an organic/inorganic, molecule/bulk interface that are commonly encountered in molecular electronics, photovoltaics, and photoelectrochemistry. The processes studied include the relaxation of the injected electron inside the TiO2 conduction band (CB), the back electron transfer (ET) from TiO2 to alizarin, the ET from the surface to the electrolyte, and the regeneration of the neutral chromophore by ET from the electrolyte to alizarin. Developing a theoretical understanding of these processes is crucial for improving solar cell design and optimizing photovoltaic current and voltage. The simulations carried out for the entire system that contains many electronic states reproduce the experimental time scales and provide detailed insights into the ET dynamics. In particular, they demonstrate the differences between the optimized geometric and electronic structure of the system at 0 K and the experimentally relevant structure at ambient temperature. The relaxation of the injected electron inside the TiO2 CB, which affects the solar cell voltage, is shown to occur on a 100 fs time scale and occurs simultaneously with the electron delocalization into the semiconductor bulk. The transfer of the electron trapped at the surface to the ground state of alizarin proceeds on a 1 ps time scale and is facilitated by vibrational modes localized on alizarin. If the electrolyte mediator is capable of approaching the semiconductor surface, it can form a stable complex and short-circuit the cell by accepting the photoexcited electron on a subpicosecond time scale. The ET from TiO2 to both alizarin and the electrolyte diminishes the solar cell current. Finally, the simulations show that the electrolyte can efficiently regenerate the neutral chromophore. This is true even though the two species do not form a chemical bond and, therefore, the electronic coupling between them is weaker than in the TiO2-chromophore and TiO2-electrolyte donor-acceptor pairs. The chromophore-electrolyte coupling can occur both directly through space and indirectly through bonding to the semiconductor surface. The ET events involving the electrolyte are promoted primarily by the electrolyte vibrational modes.     

Electrochemistry,ion adsorption and dynamics in the double layer: a study of NaCl(aq) on graphite

Graphite and related sp² carbons are ubiquitous electrode materials with particular promise for use in e.g., energy storage and desalination devices, but very little is known about the properties of the carbon–electrolyte double layer at technologically relevant concentrations. Here, the (electrified) graphite–NaCl(aq) interface was examined using constant chemical potential molecular dynamics (CμMD) simulations; this approach avoids ion depletion (due to surface adsorption) and maintains a constant concentration, electroneutral bulk solution beyond the surface. Specific Na⁺ adsorption at the graphite basal surface causes charging of the interface in the absence of an applied potential. At moderate bulk concentrations, this leads to accumulation of counter-ions in a diffuse layer to balance the effective surface charge, consistent with established models of the electrical double layer. Beyond ∼0.6 M, however, a combination of over-screening and ion crowding in the double layer results in alternating compact layers of charge density perpendicular to the interface. The transition to this regime is marked by an increasing double layer size and anomalous negative shifts to the potential of zero charge with incremental changes to the bulk concentration. Our observations are supported by changes to the position of the differential capacitance minimum measured by electrochemical impedance spectroscopy, and are explained in terms of the screening behaviour and asymmetric ion adsorption. Furthermore, a striking level of agreement between the differential capacitance from solution evaluated in simulations and measured in experiments allows us to critically assess electrochemical capacitance measurements which have previously been considered to report simply on the density of states of the graphite material at the potential of zero charge. Our work shows that the solution side of the double layer provides the more dominant contribution to the overall measured capacitance. Finally, ion crowding at the highest concentrations (beyond ∼5 M) leads to the formation of liquid-like NaCl clusters confined to highly non-ideal regions of the double layer, where ion diffusion is up to five times slower than in the bulk. The implications of changes to the speciation of ions on reactive events in the double layer are discussed.

CμMD reveals multi-layer electrolyte screening in the double layer beyond 0.6 M, which affects ion activities, speciation and mobility; asymmetric charge screening explains concentration dependent changes to electrochemical properties.     

Computer simulation of electron transfer processes across the electrode|electrolyte interface: a treatment of solvent and electrode polarizability

A polarizable molecular dynamics model for adiabatic electron transfer across the electrode|electrolyte interface is presented. The electronic polarizability of the water and of the metal electrode is accounted for by a dynamical fluctuating charge algorithm, image charges, and the Ewald summation adapted for a conducting interface. The effects of the solvent electronic polarizability are studied by computing the diabatic and adiabatic free energy curves for both polarizable and non-polarizable water models. This represents the first effort to compute the adiabatic free energy curves from simulation for a fully polarizable electrochemical system.     

Brownian dynamics simulation of protein association

Summary The Brownian Dynamics (BD) method is applied to study the diffusive dynamics and interaction of two proteins, cytochrome c (CYTC) and cytochrome c peroxidase (CYP). We examine the role of protein electrostatic charge distribution in the facilitation of protein-protein docking prior to the electron transfer step, assessing the influence of individual charged amino acid residues. Accurate interaction potentials are computed by iterating the linearized Poisson-Boltzmann (PB) equation around the larger protein CYP. The low dielectric constant inside proteins, electrolyte screening effects and irregular protein surface topography are taken into account. We observe a large ensemble of electrostatically stable encounter complexes seemingly with acceptable geometric requirements for electron transfer rather than a single dominant complex. Stabilities of the large variety of docking complexes are rationalized in terms of generalized charged residue complementarities. However, it is found that the electrostatic interactions giving rise to complex stabilities are somewhat nonspecific in nature. A large series of additional simulations are performed in which individual charged residues on CYTC have been chemically modified. Resulting perturbations of the association rate are significant and qualitatively similar to results observed in comparable kinetics experiments. We therefore demonstrate the potential of the Brownian dynamics method to estimate the effects of site-directed mutagenesis on protein-protein and protein-ligand diffusional association rates.     

Electrode polarization effects on interfacial kinetics of ionic liquid at graphite surface: An extended lagrangian-based constant potential molecular dynamics simulation study

Computational models including electrode polarization can be essential to study electrode/electrolyte interfacial phenomena more realistically. We present here a constant-potential classical molecular dynamics simulation method based on the extended Lagrangian formulation where the fluctuating electrode atomic charges are treated as independent dynamical variables. The method is applied to a graphite/ionic liquid system for the validation and the interfacial kinetics study. While the correct adiabatic dynamics is achieved with a sufficiently small fictitious mass of charge, static properties have been shown to be almost insensitive to the fictitious mass. As for the kinetics study, electrical double layer (EDL) relaxation and ion desorption from the electrode surface are considered. We found that the polarization slows EDL relaxation greatly whereas it has little impact on the ion desorption kinetics. The findings suggest that the polarization is essential to estimate the kinetics in nonequilibrium processes, not in equilibrium. © 2019 Wiley Periodicals, Inc.     

Modeling of ionic relaxation around a biomembrane disk

A simplified Brownian dynamics model and the corresponding software implementation have been developed for the simulation of electrolyte dynamics on the mesoscopic scale. In addition to direct control simulations, the model system has been verified by a quantitative comparison with the Debye-Hückel theory. As a first application, the model was used to simulate ionic relaxation processes following abrupt intramembrane charge rearrangements in the case of a disk shaped membrane. In addition to its general implications, the obtained properties of the relaxation kinetics confirm the assumptions of the theory of the so-called suspension method, a technique capable of tracing molecular charge motions of membrane proteins in three dimensions.     

Influence of thermal fluctuations on interfacial electron transfer in functionalized TiO2 semiconductors

The influence of thermal fluctuations on the dynamics of interfacial electron transfer in sensitized TiO2-anatase semiconductors is investigated by combining ab initio DFT molecular dynamics simulations and quantum dynamics propagation of transient electronic excitations. It is shown that thermal nuclear fluctuations speed up the underlying interfacial electron transfer dynamics by introducing nonadiabatic transitions between electron acceptor states, localized in the vicinity of the photoexcited adsorbate, and delocalized states extended throughout the semiconductor material, creating additional relaxation pathways for carrier diffusion. Furthermore, it is shown that room-temperature thermal fluctuations reduce the anisotropic character of charge diffusion along different directions in the anatase crystal and make similar the rates for electron injection from adsorbate states of different character. The reported results are particularly relevant to the understanding of temperature effects on surface charge separation mechanisms in molecular-based photo-optic devices.