Paul A. Madden

On the molecular origin of supercapacitance in nanoporous carbon electrodes

Lightweight, low-cost supercapacitors with the capability of rapidly storing a large amount of electrical energy can contribute to meeting continuous energy demands and effectively levelling the cyclic nature of renewable energy sources. The excellent electrochemical performance of supercapacitors is due to a reversible ion adsorption in porous carbon electrodes. Recently, it was demonstrated that ions from the electrolyte could enter sub nanometre pores, greatly increasing the capacitance. However, the molecular mechanism of this enhancement remains poorly understood. Here we provide the first quantitative picture of the structure of an ionic liquid adsorbed inside realistically modelled microporous carbon electrodes. We show how the separation of the positive and negative ions occurs inside the porous disordered carbons, yielding much higher capacitance values (125 F g(-1)) than with simpler electrode geometries. The proposed mechanism opens the door for the design of materials with improved energy storage capabilities. It also sheds new light on situations where ion adsorption in porous structures or membranes plays a role.

Molecular dynamics without effective potentials via the Car-Parrinello approach

An introduction to the Car-Parrinello molecular-dynamics method is given, written from the viewpoint of a molecular physicist. This scheme makes it possible to simulate atomic and molecular motion in cluster or bulk systems on an energy surface that is determined ‘on the fly’, by finding an adiabatic electronic state appropriate to the instantaneous nuclear positions. The crucial stability of the method is discussed at length. The development of computationally tractable expressions is considered in detail, using (local) density-functional techniques to describe the electronic energy surface. An extended account of these techniques is given, suitable for those unfamiliar with them.

Highly confined ions store charge more efficiently in supercapacitors

Liquids exhibit specific properties when they are adsorbed in nanoporous structures. This is particularly true in the context of supercapacitors, for which an anomalous increase in performance has been observed for nanoporous electrodes. This enhancement has been traditionally attributed in experimental studies to the effect of confinement of the ions from the electrolyte inside sub-nanometre pores, which is accompanied by their partial desolvation. Here we perform molecular dynamics simulations of realistic supercapacitors and show that this picture is correct at the microscopic scale. We provide a detailed analysis of the various environments experienced by the ions. We pick out four different adsorption types, and we, respectively, label them as edge, planar, hollow and pocket sites upon increase of the coordination of the molecular species by carbon atoms from the electrode. We show that both the desolvation and the local charge stored on the electrode increase with the degree of confinement.

Electrochemical interface between an ionic liquid and a model metallic electrode.

A molecular dynamics simulation model for an electroactive interface in which a metallic electrode is maintained at a preset electrical potential is described. The model, based on earlier work of Siepmann and Sprik [J. Chem. Phys. 102, 511 (1995)], uses variable charges whose magnitudes are adjusted on-the-fly according to a variational procedure to maintain the constant potential condition. As such, the model also allows for the polarization of the electrode by the electrolyte, sometimes described by the introduction of image charges. The model has been implemented in a description of an electrochemical cell as a pair of parallel planar electrodes separated by the electrolyte using a two-dimensional Ewald summation method. The method has been applied to examine the interfacial structure in two ionic liquids, consisting of binary mixtures of molten salts, chosen to exemplify the influences of dissimilar cation size and charge. The stronger coordination of the smaller and more highly charged cations by the anions prevents them from approaching even the negatively charged electrode closely. This has consequences for the capacitance of the electrode and will also have an impact on the rates of electron transfer processes. The calculated capacitances exhibit qualitatively the same dependence on the applied potential as has been observed in experimental studies.

Properties of liquid CS2 from the allowed light scattering spectra

Measurements of the anisotropic Rayleigh and the allowed Raman spectra of carbon disulphide over the normal liquid range are described. The Raman reorientation times agree well with N.M.R. values and are qualitatively well predicted by the Stokes-Einstein-Debye formula with slip boundary conditions, although significant departures are found. The Hubbard relationship is obeyed. A careful comparison of the Rayleigh and anisotropic Raman intensities indicates that the orientational correlation parameter (g 2) decreases with increasing temperature. Using values of g 2 obtained in this way the ratio of the Rayleigh and Raman reorientation times is shown to be consistent with a small value for the correlation in angular velocity between molecules. The vibrational dephasing time is found to be roughly constant with temperature in contrast with the predictions of the theories for this quantity. An attempt is made to correlate its behaviour with that of the collision induced Raman spectra of CS2.

Nonlinear counterion screening in colloidal suspensions

A new ‘‘ab initio’’ method is presented which is designed to simulate highly asymmetric systems of charged particles such as micellar solutions and charge‐stabilized colloidal suspensions. The hybrid description considers the macroion degrees of freedom explicitly, while the microscopic counterions are treated within the framework of density functional theory. The counterion density profile is treated as a dynamical variable which is coupled to the macroion positions; the corresponding equation of motions are derived from a Lagrangian which contains a fictitious kinetic energy term associated with the inhomogeneous counterion density, with a fictitious mass chosen so that the counterions stay as close as possible to the surface of lowest free energy (adiabatic condition). The discontinuous behavior of the counterion density profile at the macroion surfaces is suppressed by the use of a classical pseudopotential scheme without spoiling the rapid variation of the counterion density profile outside the macro...

Hydration of metal surfaces can be dynamically heterogeneous and hydrophobic

We present a study of the solvation properties of model aqueous electrode interfaces. The exposed electrodes we study strongly bind water and have closed packed crystalline surfaces, which template an ordered water adlayer adjacent to the interface. We find that these ordered water structures facilitate collective responses in the presence of solutes that are correlated over large lengthscales and across long timescales. Specifically, we show that the liquid water adjacent to the ordered adlayers forms a soft, liquid-vapor-like interface with concomitant manifestations of hydrophobicity. Temporal defects in the adlayer configurations create a dynamic heterogeneity in the degree to which different regions of the interface attract hydrophobic species. The structure and heterogeneous dynamics of the adlayer defects depend upon the geometry of the underlying ordered metal surface. For both 100 and 111 surfaces, the dynamical heterogeneity relaxes on times longer than nanoseconds. Along with analyzing time scales associated with these effects, we highlight implications for electrolysis and the particular catalytic efficiency of platinum.We have applied molecular dynamics and methods of importance sampling to study structure and dynamics of liquid water in contact with metal surfaces. The specific surfaces considered resemble the 100 and 111 faces of platinum. Several results emerge that should apply generally, not just to platinum. These results are generic consequences of water molecules binding strongly to surfaces that are incommensurate with favorable hydrogen-bonding patterns. We show that adlayers of water under these conditions have frustrated structures that interact unfavorably with adjacent liquid water. We elucidate dynamical processes of water in these cases that extend over a broad range of timescales, from less than picoseconds to more than nanoseconds. Associated spatial correlations extend over nanometers. We show that adlayer reorganization occurs intermittently, and each reorganization event correlates motions of several molecules. We show that soft liquid interfaces form adjacent to the adlayer, as is generally characteristic of liquid water adjacent to a hydrophobic surface. The infrequent adlayer reorganization produces a hydrophobic heterogeneity that we characterize by studying the degrees by which different regions of the adlayers attract small hydrophobic particles. Consequences for electrochemistry are discussed in the context of hydronium ions being attracted from the liquid to the metal–adlayer surface.

Simulating Supercapacitors: Can We Model Electrodes As Constant Charge Surfaces?

Supercapacitors based on an ionic liquid electrolyte and graphite or nanoporous carbon electrodes are simulated using molecular dynamics. We compare a simplified electrode model in which a constant, uniform charge is assigned to each carbon atom with a realistic model in which a constant potential is applied between the electrodes (the carbon charges are allowed to fluctuate). We show that the simulations performed with the simplified model do not provide a correct description of the properties of the system. First, the structure of the adsorbed electrolyte is partly modified. Second, dramatic differences are observed for the dynamics of the system during transient regimes. In particular, upon application of a constant applied potential difference, the increase in the temperature, due to the Joule effect, associated with the creation of an electric current across the cell follows Ohms law, while unphysically high temperatures are rapidly observed when constant charges are assigned to each carbon atom.

Polarization effects in ionic systems from first principles

A classical application of the Car-Parrinello method in the computer simulation of ionic systems is demonstrated. The induction effects in the interionic interactions are included in addition to the short-range repulsion and dispersion effects described by effective pair potentials. By representing the induced dipoles as a pair of dynamically variable charges fixed on the ends of a rod and extending the Lagrangian accordingly, the self-consistent induced dipoles at each time step are generated from the values at the previous time step, without the need for explicit minimization. Coulomb-field-induced and overlap-induced polarization effects are included and these are parameterized by ab initio electronic structure calculations. For simple ionic systems it is shown that the neglect of the overlap-induced dipoles leads to a poor representation of real systems.

The electric double layer has a life of its own

Using molecular dynamics simulations with recently developed importance sampling methods, we show that the differential capacitance of a model ionic liquid based double-layer capacitor exhibits an anomalous dependence on the applied electrical potential. Such behavior is qualitatively incompatible with standard mean-field theories of the electrical double layer, but is consistent with observations made in experiment. The anomalous response results from structural changes induced in the interfacial region of the ionic liquid as it develops a charge density to screen the charge induced on the electrode surface. These structural changes are strongly influenced by the out-of-plane layering of the electrolyte and are multifaceted, including an abrupt local ordering of the ions adsorbed in the plane of the electrode surface, reorientation of molecular ions, and the spontaneous exchange of ions between different layers of the electrolyte close to the electrode surface. The local ordering exhibits signatures of a first-order phase transition, which would indicate a singular charge-density transition in a macroscopic limit.