New Perspectives on the Charging Mechanisms of Supercapacitors

Mar 31, 2016 - Supercapacitors (or electric double-layer capacitors) are high-power energy storage devices that store charge at the interface between ...
0 downloads 13 Views 6MB Size
This is an open access article published under a Creative Commons Attribution (CC-BY) License, which permits unrestricted use, distribution and reproduction in any medium, provided the author and source are cited.

Perspective pubs.acs.org/JACS

New Perspectives on the Charging Mechanisms of Supercapacitors Alexander C. Forse,*,† Céline Merlet,† John M. Griffin,†,‡ and Clare P. Grey*,† †

Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K. Department of Chemistry, Lancaster University, Lancaster LA1 4YB, U.K.



This physical mechanism of charge storage gives rise to fast charge and discharge times and long cycle lives, characteristic properties that make supercapacitors attractive devices to complement batteries (which can store and deliver more energy but with slower charge and discharge times). Today, supercapacitors are used in a range of industrial, automotive, and electric utility applications including electric buses, trains, uninterruptible power supply systems, elevators, camera flashes, cranes, and engine starters.1,2 Their more widespread use could be facilitated by the development of new devices with improved energy densities, which retain the high power densities and long cycle lives that are characteristic of supercapacitors. Typical materials for supercapacitors are highlighted in Figure 1. Porous carbon electrode materials are generally prepared by the heat treatment and subsequent chemical activation of organic materials, such as coconut shells and wood,3 while a related class of materials, “carbide-derived carbons” (CDCs), are obtained from metal carbides by extracting the metal atoms.4 More exotic materials such as carbon nanotubes (CNTs) and graphenes are also being developed for supercapacitor application, but here we will focus our attention on disordered porous carbons (activated carbons) as they are well-studied and widely used in commercial devices due to their cheap price, facile synthesis, and sustainability. For the electrolyte, the most widely used systems are comprised of salts dissolved in organic solvents (e.g., tetraethylammonium tetrafluoroborate in acetonitrile solvent, NEt4−BF4/ACN). Such organic electrolytes offer a good balance of relatively large maximum operating voltages (∼2.5 V) and high ionic conductivities (∼20−60 mS·cm−1). The stored energy, E, is given by

ABSTRACT: Supercapacitors (or electric double-layer capacitors) are high-power energy storage devices that store charge at the interface between porous carbon electrodes and an electrolyte solution. These devices are already employed in heavy electric vehicles and electronic devices, and can complement batteries in a more sustainable future. Their widespread application could be facilitated by the development of devices that can store more energy, without compromising their fast charging and discharging times. In situ characterization methods and computational modeling techniques have recently been developed to study the molecular mechanisms of charge storage, with the hope that better devices can be rationally designed. In this Perspective, we bring together recent findings from a range of experimental and computational studies to give a detailed picture of the charging mechanisms of supercapacitors. Nuclear magnetic resonance experiments and molecular dynamics simulations have revealed that the electrode pores contain a considerable number of ions in the absence of an applied charging potential. Experiments and computer simulations have shown that different charging mechanisms can then operate when a potential is applied, going beyond the traditional view of charging by counter-ion adsorption. It is shown that charging almost always involves ion exchange (swapping of co-ions for counter-ions), and rarely occurs by counter-ion adsorption alone. We introduce a charging mechanism parameter that quantifies the mechanism and allows comparisons between different systems. The mechanism is found to depend strongly on the polarization of the electrode, and the choice of the electrolyte and electrode materials. In light of these advances we identify new directions for supercapacitor research. Further experimental and computational work is needed to explain the factors that control supercapacitor charging mechanisms, and to establish the links between mechanisms and performance. Increased understanding and control of charging mechanisms should lead to new strategies for developing next-generation supercapacitors with improved performances.

E=

(1)

where C is the cell capacitance, and V is the operating voltage. Thus, organic electrolytes are typically preferred to aqueous electrolytes that are limited to ∼1 V before water decomposes.5 Aqueous-based systems are, however, being studied for applications where cost is a critical parameter (e.g., on the electricity grid).6−8 Room-temperature ionic liquids are also emerging as alternative electrolytes for supercapacitors, with operating voltages as high as 4 V achievable.5,9−11 However, the increase in energy by eq 1 comes at a cost, as slower ionic transport (ionic conductivities are typically below 20 mS·cm−1) results in poorer device power performances. Most efforts to increase the energy density of supercapacitors have focused on the development of new carbon materials with increased capacitances, so that more energy can be stored, as

1. INTRODUCTION Supercapacitors (strictly, electric double-layer capacitors) store charge at the interface between porous carbon electrodes and an electrolyte solution (Figure 1). In contrast to batteries, charge storage in supercapacitors is non-faradaic and occurs by the physical adsorption and desorption of ions inside the pores of the carbon electrodes when an external voltage is applied. As electronic charge accumulates in an electrode, it is balanced at the interface by an equal and opposite ionic charge in the electrolyte. © 2016 American Chemical Society

1 2 CV 2

Received: February 25, 2016 Published: March 31, 2016 5731

DOI: 10.1021/jacs.6b02115 J. Am. Chem. Soc. 2016, 138, 5731−5744

Perspective

Journal of the American Chemical Society

Figure 1. Schematic view of a supercapacitor. Porous carbon materials with disordered structures are used as the electrodes, and the cell is soaked with an electrolyte that may be organic, aqueous or ionic liquid-based, with some typical electrolytes shown. Note, for simplicity the separator (which prevents short circuit), the binder that holds the electrode materials together and the current collectors are not shown. Schematic porous carbon structure adapted from ref 17 with permission from Springer.

observed as the carbon pore size was decreased below 1 nm. These results showed the importance of pore size, and challenged the previous view that pores smaller than the solvated electrolyte ions do not contribute significantly to the capacitance. It was hypothesized that in small pores ion desolvation allows a closer approach of charge centers at the electrode−electrolyte interface, which increases the capacitance.24,26 However, the capacitance increase was also observed in experimental,27 and theoretical,28−32 studies of ionic liquids, with the capacitance maximized when the pore size matched the ion size. It was even shown that the capacitance varies in an oscillatory fashion as the (slit-)pore width is varied.28−32 These findings indicate that factors beyond simple ion desolvation arguments are responsible for the anomalous capacitance increase. Interestingly, some experimental studies have reported that there is no correlation between pore size and capacitance,33,34 with the origin of these differences currently unclear. Needless to say, the observation of pore size effects led to a wave of interest in the mechanisms of charge storage in porous carbon electrodes. In the past 10 years, pioneering experimental and theoretical studies have led to an improved understanding of the molecular mechanisms of charge storage in supercapacitors. This perspective seeks to bring together recent findings from a range of studies to provide a coherent and detailed view of modern theories of supercapacitor charging mechanisms. Recent studies have shown that charging is not simply driven by adsorption of counter-ions into the electrode pores, as was previously believed (note, counter-ions are defined as having opposite charge to the electrode). NMR experiments and molecular dynamics (MD) simulations show that porous carbon electrodes typically contain a large number of ions in the absence of an applied potential, such that a range of different charging mechanisms, based on different amounts of ion adsorption and desorption, are then possible when a voltage is applied. New in situ experimental techniques (NMR spectroscopy, electrochemical quartz crystal microbalance (EQCM), infrared (IR) spectroscopy, and scattering approaches) have shown that the charging mechanism is sensitive to the electrode and electrolyte materials used, as well as the polarization of the electrode. Understanding and controlling the charge storage mechanism of supercapacitors may hold the key to developing next generation devices with

shown by eq 1. However, the success of such an approach requires an understanding of the carbon structure, and how this in turn affects the charge storage mechanism and the capacitance. This is a complex problem, as porous carbon materials lack long-range order, making their characterization challenging. These materials do, however, exhibit order on a local scale (at length scales up to 10 or 20 Å). While analysis of the broad Bragg peaks in diffraction experiments offers little information, pair distribution function (PDF) analysis and nuclear magnetic resonance (NMR) experiments show that typical activated carbons and CDCs consist of predominantly sp2-hybridized carbon atoms organized in a hexagonal arrangement.12−15 Recent studies have suggested that non-hexagonal rings, such as 5-membered16 and 7-membered rings,15 are also present, giving rise to curvature in the carbon sheets (see schematic structure in Figure 1).17,18 The curved and defective carbon sheets and fragments do not pack together well (structures are typically non-graphitizing, with graphite not formed even on heating to temperatures as high as 3000 °C),19,20 and nanometersized pores exist between the various carbon surfaces (see Figure 1). This porosity is typically characterized by pore size distributions obtained from the analysis of gas sorption isotherms. State of the art methods use classical density functional theory (DFT) to obtain the pore size distribution.21−23 The recently developed NLDFT (non-local DFT) and QSDFT (quenched solid DFT) methods both rely on the assumption that the porous material is a collection of pores with identical shapes (typically slit-shaped pores), and with different pore widths. Classical DFT is used to determine gas adsorption profiles inside the pores, and a set of partial isotherms corresponding to different pore sizes are generated. The pore size distribution is then obtained by fitting the experimental isotherm to a sum of partial isotherms. Owing to the relative ease of this technique, pore size distributions are the most frequently used metric in the characterization of carbon structure, and have been used as the primary tool in the search for relationships between carbon structure and capacitance. Seminal studies in 2006 showed that the carbon capacitance could be increased by optimizing the pore size of the carbon electrodes.24,25 An “anomalous” increase in capacitance was 5732

DOI: 10.1021/jacs.6b02115 J. Am. Chem. Soc. 2016, 138, 5731−5744

Perspective

Journal of the American Chemical Society

Figure 2. (a) Schematic illustrating the concept of filled and empty carbon pores at 0 V. The carbon (slit-)pore walls are represented by black rectangles. (b) Snapshot of an MD simulation showing the presence of in-pore ions and solvent molecules at 0 V. Red, blue, and green correspond to cations, anions, and solvent molecules, respectively, while the carbon surfaces are gray. See ref 40 for details of the MD study. (c) NMR (9.4 T) measurements of YP50F activated carbon soaked with a typical supercapacitor electrolyte, NEt4−BF4/dACN (1.5 M), recorded with magic angle spinning at a frequency of 5 kHz. Note, deuterated acetonitrile (dACN) was used here to allow a convenient separation of the signals from the cations and the solvent molecules. MD simulations and NMR experiments both reveal a significant number of in-pore ions in the absence of an applied potential.

sheet planes, inside curved hollows, and inside small pockets.40 BMI−PF6 dissolved in ACN to 1.5 M (i.e., an organic electrolyte) has also been studied under nanoconfinement, with a considerable number of in-pore ions again observed (Figure 2b). In this case the solvent molecules replace some of the in-pore ions, leading to a lower in-pore population relative to the neat ionic liquid study, with the small solvent molecules tending to occupy the most confined pore sites. Other theoretical studies based on idealized electrode geometries, such as CNTs41,42 and slit pores,28,43 have also shown spontaneous pore-filling behavior. Experimental results from NMR spectroscopy are in excellent accord with the findings from MD simulations, also revealing the predominance of filled pores in the absence of an applied potential (see Table 1). In NMR experiments, ions and solvent molecules inside carbon pores (“in-pore”) give rise to peaks which are distinct to those from species in bulk electrolyte between the carbon particles (“ex-pore”).14,44 By studying different NMR nuclei, the anions, cations, and solvent molecules can be studied separately. For example, for YP50F (a commercial activated carbon) soaked with the electrolyte NEt4−BF4 in ACN (1.5 M), clear resonances can be distinguished for in-pore

enhanced properties, and we identify new areas of research that can facilitate this process.

2. CHARACTERIZING THE ELECTRODE−ELECTROLYTE INTERFACE AT 0 V Before embarking on the study of supercapacitor charging mechanisms, it is crucial to have a detailed understanding of the structure of the electrode−electrolyte interface in the absence of an applied potential (i.e., at 0 V). As illustrated by Figure 2a, the carbon nanopores may be either filled with electrolyte ions, or they may contain no ions. These two possibilities have been described as “ionophilic” and “ionophobic” pores by Kornyshev and Kondrat (see later).35−37 Together, MD simulations and NMR spectroscopy experiments have shown that the pores of the carbon electrodes are generally filled at 0 V. Coarse-grained MD simulations of the ionic liquid butylmethylimidazolium hexafluorophosphate (BMI−PF6) showed a large number of ions inside the carbon nanopores of realistic carbon38 (CDC) structures at zero applied potential (Table 1).39 With the MD approach, subtly different local adsorption environments could be identified for confined ions, sitting at sheet edges, on top of 5733

DOI: 10.1021/jacs.6b02115 J. Am. Chem. Soc. 2016, 138, 5731−5744

Perspective

Journal of the American Chemical Society Table 1. In-Pore Ionic Populations at 0 V from MD Simulations and NMR Spectroscopya method

electrolyte

carbon

pore size (average)/nm

in-pore population at 0 V/mmol g−1 carbon

ref

NMR NMR MD MD MD MD MD

EMI−TFSI ionic liquid Pyr13−TFSI ionic liquid BMI−PF6 ionic liquid EMI−TFSI ionic liquid EMI−TFSI ionic liquid EMI−Cl ionic liquid EMI−BF4 ionic liquid

YP50F YP50F CDC slit pore slit pore MWCNT CNT

1.01 1.01 0.93 1.10 1.00 1.36 1.08

1.8 1.6 1.87 1.74 1.43 1.8b 0.7

53 53 39 43 28 41 42

NMR NMR MD NMR NMR NMR

EMI−TFSI/ACN (1.9 M) Pyr13−TFSI/ACN (1.8 M) BMI−PF6/ACN (1.5 M) PEt4−BF4/ACN (1.5 M) PEt4−BF4/ACN (0.75 M) PEt4−BF4/ACN (0.5 M)

YP50F YP50F CDC YP50F YP50F YP50F

1.01 1.01 0.93 1.01 1.01 1.01

1.2 1.1 0.68 0.86 0.47 0.31

53 53 40 52 52 52

NMR NMR

NaF aq (0.8 mol kg−1) NaF aq (0.8 mol kg−1)

PEEK AC PEEK AC

1.55 0.58

>0 0

55 55

a

In the NMR studies the commercially activated carbon YP50F was studied in all cases, except for in ref 55, where activated carbons derived from poly(ether-ether-ketone) (PEEK AC) were studied. In most of the MD simulations listed, idealized carbon geometries such as slit pores and multiwalled carbon nanotubes (MWCNTs) were studied, with the exception of the studies of model CDC electrodes in refs 39 and 40. bFor the calculation of the gravimetric in-pore population from ref 41, only the inner wall of the MWCNT was considered for the mass calculation. Note: Pyr13−TFSI is 1-methyl-1-propylpyrrolidinium bis(trifluoromethanesulfonyl)imide, and EMI−TFSI is 1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide.

anions (19F NMR), cations (1H NMR), and solvent (2H NMR) (Figure 2c). In each case the in-pore peak appears at a negative chemical shift (relative to the ex-pore) due to the local magnetic field originating from the circulation of the carbon’s delocalized π electrons in the applied magnetic field.14,45,46 These effects are largely independent of the choice of NMR nucleus, such that anions, cations, and solvent molecules all experience a similar change of chemical shift upon adsorption in a given carbon. Since NMR is quantitative, the number of ions inside the carbon pores can be readily determined by fitting and integration of the spectra. NMR studies have been carried out on carbons soaked with a range of organic electrolytes.14,44,47−52 For PEt4−BF4 in ACN (1.5 M) and the carbon YP50F, there was 0.86 mmol of in-pore anions and cations per gram of carbon, with this number decreasing at lower electrolyte concentrations (Table 1).52 NMR studies of ionic liquids and YP50F carbon have revealed larger numbers of in-pore ions at 0 V (Table 1).53 When the ionic liquid samples were diluted with ACN (to ∼1.8 M), the in-pore ion populations dropped as solvent molecules displaced some of the ions (Table 1).53 This is consistent with MD simulations on model CDC electrodes with and without a solvent present (Table 1). The NMR approach has also recently been applied to study the confinement of aqueous electrolytes, highlighting the versatility of the method.54,55 It was shown that for a carbon with a relatively small average pore width (0.58 nm) soaked with NaF (0.8 mol kg−1 aq), the ions were unable to enter the carbon pores in the absence of an applied potential, despite the presence of in-pore water.55 Measurements on a carbon with larger pores (average width 1.55 nm) showed there were a considerable number of in-pore ions, confirming that the absence of in-pore ions in the former carbon arose from steric factors. In the absence of any steric effects, it is currently unknown if it is possible to synthesize carbons with ionophobic pores (see later). We note that the study of the electrode− electrolyte interface in aqueous systems is highly complex,

as the chemical nature of the carbon surface will depend on the (de)protonation of the various functional groups, and OH− and H3O+ ions are present alongside the main electrolyte ions. Beyond NMR spectroscopy and MD simulations, small-angle X-ray scattering (SAXS) and small-angle neutron scattering (SANS) have also emerged as experimental methods to probe the wetting of carbon pores in the absence of an applied potential.56−60 Changes of the scattered neutron intensity after the addition of ACN to activated carbon fabric suggested that the solvent had wetted the carbon nanopores.56 Interestingly, this approach allows wetting to be studied as a function of pore size, and it was shown that the pore wetting of activated carbon fabrics was incomplete for the smallest nanopores.56 In contrast to NMR and MD simulations, however, these methods have not yet allowed an absolute quantification of the numbers of in-pore anions, cations, and solvent molecules.

3. STUDIES OF SUPERCAPACITOR CHARGING MECHANISMS 3.1. Possible Mechanisms of Charge Storage. Given the large number of ions inside the carbon pores at 0 V, a range of different charging mechanisms are possible (Figure 3).49 First, charge may be balanced by the adsorption of counter-ions. This is the traditional view of charging. A second possibility is that counter-ion adsorption is accompanied by simultaneous co-ion desorption from the pores, which we refer to as ion exchange (where co-ions are defined as having charge with the same sign as the electrode). A third possibility is that charging is driven purely by the desorption of co-ions (Figure 3). In each case the excess ionic charge inside the carbon pores is equal and opposite to the electronic charge stored in the carbon. In practice, charging may involve a combination of the different mechanisms shown in Figure 3. For example, ion-exchange and counter-ion adsorption could occur simultaneously. 5734

DOI: 10.1021/jacs.6b02115 J. Am. Chem. Soc. 2016, 138, 5731−5744

Perspective

Journal of the American Chemical Society

Figure 3. Different possible charging mechanisms for carbon pores that are initially filled with electrolyte: counter-ion adsorption, ion exchange, and co-ion desorption. The different charging mechanisms may be described by the charging mechanism parameter, X, as defined in eq 2. Example calculations of X (using eq 2) are shown for the three depicted charging mechanisms, with V0 taken as 0 V. A value of X = +1 is obtained for charging solely by counter-ion adsorption, while X = 0 is obtained for ion exchange, and X = −1 for co-ion desorption. As indicated by the scale on the right, X is a continuous variable (e.g., an X value intermediate between 1 and 0 would indicate that both io-exchange and counter-ion adsorption processes occur during charging). Part of this figure is adapted from ref 49 with permission from The Royal Society of Chemistry.

and experimental techniques, and the mechanism depends on the electrode and electrolyte materials studied. 3.2. Computational Methods. A schematic illustrating some of the main approaches that have been developed to study the charging mechanisms of supercapacitors is shown in Figure 4. Computational methods have been extremely successful in this area as they allow one to probe local length scales. The simulation of supercapacitors and their charging mechanisms presents a number of challenges which have been tackled through different approaches ranging from mean-field theories to MD simulations. One of the main difficulties is related to the complex structure of the carbon electrodes. In mean-field theories, electrodes are usually represented as single slit or cylindrical pores. In molecular simulations, it is also difficult to account for the complexity of the carbon structure as (i) an accurate description of the structure at the atomistic scale is not readily available and (ii) a realistic representation of the carbon structure requires the use of a large number of atoms in the simulations, which is time consuming. Nevertheless, MD simulations are valuable as they provide a precise description of the ion−ion correlations and packing effects that are important in the context of the concentrated electrolytes used in supercapacitors. Another challenging aspect is the electronic conductivity of the carbon electrodes. Experiments and simulations have shown that the materials and conditions used in the carbon synthesis have a large effect on the ordering of the carbon sheets.12,13,15,38,61−63 This ordering, as well as the presence of functional groups, impacts the electronic conductivity of the carbons and their performances as electrode materials.62,64,65 In molecular simulations, there are two main approaches to deal with the electrode charge: (i) constant charge and (ii) constant potential.66,67 In the constant charge approach, charges are assigned to the carbon atoms at the beginning of the simulations and kept constant throughout. As such, this approach neglects the existence of charge redistribution on the electrode. The converse approach is to consider the electrode as metallic and include the fluctuation of carbon atom charges in the simulations.68 While this approach

It is useful to describe the charging mechanism with a single mathematical quantity. Here we introduce the charging mechanism parameter, X(V,V0), defined as X(V , V0) =

N (V ) − N (V0) (Ncounter(V ) − Nco(V )) − (Ncounter(V0) − Nco(V0))

=

N (V ) − N (V0) ([|Q ionic(V )| − |Q ionic(V0)|]/e)

(2)

where N(V) is the total number of in-pore ions at a charging voltage V, N(V0) is the total number of in-pore ions at some initial voltage V0 (typically 0 V), and Ncounter(V) and Nco(V) are the number of in-pore counter-ions and co-ions, respectively, at a voltage V. Qionic(V) and Qionic(V0) are then the net in-pore ionic charges at the two voltages, and e is the elementary charge. Put simply, X gives a measure of the roles of counter-ion adsorption, ion exchange, and co-ion desorption in the charging mechanism. For the classical charging mechanism of counter-ion adsorption, X = +1. For ion exchange, X = 0, while for co-ion desorption, X = −1 (see example calculations in Figure 3). Intermediate values of X are also possible; for example, X = 0.3 would indicate that both ion exchange and counter-ion adsorption occur during charging, with ion exchange dominating (as 0.3 is closer to 0 than it is to 1). We stress that this parameter refers to charging; i.e., |V| must be greater than |V0|. It is also important to realize that X can take different values in the positive and negative electrodes, such that separate calculations must be performed for each, and that X may also depend on the studied voltage range. We note that for pores that are initially empty, charging must at least initially proceed by counter-ion adsorption (X = 1), as there are no in-pore co-ions available for desorption. For such ionophobic pores, charging may also involve the net movement of pairs of ions into the carbon pores. In this unusual case, X can take values greater than 1. As we will see in the next section, the different charging mechanisms can be probed with computational 5735

DOI: 10.1021/jacs.6b02115 J. Am. Chem. Soc. 2016, 138, 5731−5744

Perspective

Journal of the American Chemical Society

Figure 4. Schematic showing the different approaches for studying the charging mechanisms of supercapacitors, and their advantages and challenges.

above, this is very challenging, as it involves using a large number of atoms for the carbon porous structure and a constant potential approach, which are both computationally expensive. To date, only a few studies have included both of these characteristics.39,40,73,74 MD simulations have highlighted how the lack of “overscreening” effects in nanoporous electrodes contributes to the anomalous capacitance increase observed in nanometer-sized pores.39 In nanoporous electrodes, only a single layer (or a few layers) of adsorbed ions is present between the pore walls.39 This is in stark contrast to planar electrode surfaces, such as the graphite surface, where several layers of ions of alternating charge extend from the charged surface into the bulk.75,76 At planar surfaces, the first layer of adsorbed counterions carries a greater charge than the electrode surface (“overscreening” the surface charge), and subsequent layers of co-ions and counter-ions then balance the excess charge in the first layer. The lack of overscreening in nanoporous electrodes does not, however, account for all of the anomalous capacitance, and hence a more atomistic/molecular understanding of charge storage is required. To this end, MD simulations of a pure ionic liquid (BMI−PF6) and CDC electrodes have shown that charge storage does not occur by counter-ion adsorption alone (Table 2).39 In the positive electrode, charging occurred by a combination of ion exchange and counter-ion adsorption (X = 0.34), while in the negative electrode charging was mainly by ion exchange, with a small amount of co-ion desorption (X = −0.11). Thus, it appears that the PF6 anions are the more “active” species here, playing the dominant role in the charge storage process in both electrodes. Similar results were found for the same disordered carbon with the organic electrolyte (BMI−PF6/ACN), showing that the solvent does not significantly affect the capacitance or charging mechanism (Table 2).39,40 For a model supercapacitor system with slit-pore electrodes and

is likely to be more realistic, especially for carbons with high degrees of order,62 it is much more computationally expensive and has been used in only a limited number of studies. Further work is needed to assess the validity of these different approaches and also to develop quantum mechanical methods to treat the carbon charge. One way of circumventing the problem of computational cost is to explore the effects of confinement and polarization through analytical expressions. Indeed, while early theories considered ions as point charges and electrolytes as dilute, new advances in this field have allowed the inclusion of steric interactions and charge screening in the representation of the systems.31,69 Such mean-field theories predicted the “superionic effect” which states that the packing of ions of the same charge is easier in confined spaces because of an exponential screening of the electrostatic interactions by the charged pore wall. This charge screening has recently been explored in more detail by inserting atoms in gold nanotubes and CNTs using DFT calculations.70−72 The authors showed that the atoms were fully ionized upon insertion, turning the initially semi-conducting nanotubes into conductors. Both analytical theories and quantum chemistry simulations have provided valuable insights into the understanding of charging mechanisms by investigating the impact of charge screening. Nevertheless, the specific performances of different electrolytes, with or without solvent, can be understood only via techniques that explicitly account for ionic correlations, and the size and shape of the electrolyte molecules. Classical MD simulations are currently the method of choice here, although MD-DFT simulations will become increasingly important as computational power increases. Ideally, MD simulations of supercapacitors would include a realistic representation of the carbon electrode structure, as well as a description of its electronic conductivity. As mentioned 5736

DOI: 10.1021/jacs.6b02115 J. Am. Chem. Soc. 2016, 138, 5731−5744

Perspective

Journal of the American Chemical Society

where it was shown that charging mechanisms are different from slit pores due to local heterogeneities in the carbon structure.77,78 3.3. Electrochemical Quartz Crystal Microbalance (EQCM). Experiments using an EQCM were some of the first to investigate the molecular mechanisms of supercapacitor charging. Levi et al. demonstrated that by depositing porous carbon electrodes on a quartz crystal, it was possible to measure changes of electrode mass during the operation of a supercapacitor cell.79 In practice, variations of the resonance frequency of the quartz crystal are measured, and these are converted into mass changes using the Sauerbrey equation.80 These experiments have mostly been performed with dynamic charging conditions, with the voltage continuously swept between the voltage limits at a constant rate (i.e., cyclic voltammetry). The extracted mass changes are then typically compared to the predicted mass changes assuming the adsorption of a single counter-ion to balance each electronic charge in the electrode (i.e., assuming X = 1). While this assumption generally does not hold, deviations from the predicted mass changes can be used to infer the presence of ion-exchange, and qualitative studies of the supercapacitor charging mechanism are possible. Initial studies of activated carbon electrodes with an aqueous CsCl (0.5 M) electrolyte suggested that counter-ion adsorption processes dominated the charging mechanism. In contrast, further experiments on tetrabutylammonium chloride (NBu4Cl, 0.5 M aq), with the bulky NBu4 cation, suggested that charging was brought about predominantly by the migration of chloride anions (i.e., counter-ion adsorption in the positive electrode, and co-ion desorption in the negative electrode), with these experiments providing the first evidence that the charging mechanism depends on the electrolyte.79 In follow-on studies on organic electrolytes, different charging mechanisms could be discerned depending on the magnitude of the charge stored in the electrode (i.e., depending on the cell voltage), see Figure 5a.81 For low charge densities (region I), the measured mass change was smaller than that expected for pure counter-ion adsorption processes, indicating the presence of ion exchange. For larger (intermediate) charge densities (region II), the measured mass changes were consistent with a counter-ion adsorption mechanism. Finally, for high charge densities (region III), the measured mass changes were larger than those predicted for pure counter-ion adsorption, and it was suggested the source of the additional mass was extra solvent molecules that entered the electrode pores. This is surprising given that the pores are already densely packed with ions and solvent molecules, and we note that MD simulations have suggested that the in-pore solvent population does not significantly change during charging.74,82 Further studies are thus needed to investigate how the in-pore solvent population varies during charging. Interestingly, when electrolytes with different cations were studied with EQCM (Figure 5b), it was shown that the ion-exchange processes in the negative electrode became more significant for electrolytes with bulkier cations (e.g., NBu4−BF4), with the role of desorption of the smaller anions becoming more important.81 More recently, EQCM methods have been applied to devices with EMI−TFSI ionic liquid and CDC electrodes.83 Again, the charging mechanism was found to depend on the magnitude of the charge stored on the electrodes. In the positive electrode ion exchange was observed at low charge densities, with counter-ion adsorption processes then dominating at higher charge densities. In the negative electrode, on the other hand, counter-ion adsorption appeared to dominate over the full range of studied charge densities. Going further, the ionic liquids were diluted

Table 2. Details from MD Simulations at a Charging Voltage of 1 Va charging mechanism parameter, X carbon, electrolyte CDC, BMI−PF6 ionic liquid CDC, BMI−PF6/ACN (1.5 M) slit pore (1.1 nm), EMI−TFSI

Ncounter-ion − Nco‑ion positive negative (at 1 V)/mmol g−1 electrode electrode ref 1.12 1.03 0.33

0.34 0.57 −0.36

−0.11 −0.02 0.45

39 40 43

Ncounter-ion − Nco‑ion is proportional to the excess ionic charge (which is in turn directly proportional to gravimetric capacitance). Charging mechanism parameters, X, are also given.

a

the ionic liquid 1-ethyl-3-methylimidazolium− bis(trifluoromethanesulfonyl)imide (EMI−TFSI), charging in the positive electrode involved both ion exchange and co-ion desorption (X = −0.36), while charging in the negative electrode involved both ion exchange and counter-ion adsorption (X = 0.45); i.e., the EMI cations are more “active” in this system.43 These observations for the different MD studies indicate that ion packing effects and ion−ion interaction energies influence the charging mechanism that operates. We note that, for an identical potential difference and a similar pore size, the difference between the number of anions and cations in the charged electrodes (which is proportional to the capacitance) is 4 times larger in the disordered CDC electrodes compared to the single slit pores (Table 2). This is partly because both sides of the carbon sheets are accessible to ions in the disordered electrodes while only a single surface is accessible in the slit-pore system. Accounting for this, the capacitance is still twice as large for the disordered electrodes, indicating that the curved and defective nature of the carbon sheets facilitates charge separation. It has been shown that more confined “pocket-like” ion adsorption sites allow the more effective storage of charge, than “plane” sites,40 presumably because the curved surfaces maximize favorable Coulombic interactions between the carbon surfaces and the counter-ions (while also screening repulsive interactions between different counter-ions). These effects result in a much larger capacitance for disordered porous electrodes and highlight the importance of using complex porous carbon structures in simulations of supercapacitors. Besides providing a microscopic picture and fundamental understanding of the local processes occurring during supercapacitor charging, both mean-field theories and molecular simulations can be used to explore new ideas to improve energy storage efficiency in supercapacitors. An interesting example of these hypothetical explorations is the concept of ionophobicity and ionophilicity. Simulations can be used to investigate the effects of having initially empty (ionophobic) pores, on both the energy stored and the charging rate.35−37 Both mean-field theories and MD simulations for single slit pores have shown that ionophobic pores perform better in terms of charging rate,36 and suggest that, in some cases, the quantity of energy stored is higher for ionophobic pores.37 One interesting aspect is that the charging of initially filled pores usually leads to an overfilling of the porosity, corresponding to a temporary state where the density of ions is higher than the final density at the end of charging. This type of phenomena raises the issues of kinetic barriers and the difference between static and dynamic charging of supercapacitors. While the effects of ionophobicity and ionophilicity are now well studied in idealized slit-pore geometries, it would be interesting to study similar ideas in realistic three-dimensional porous geometries 5737

DOI: 10.1021/jacs.6b02115 J. Am. Chem. Soc. 2016, 138, 5731−5744

Perspective

Journal of the American Chemical Society

system. A detailed overview of the in situ NMR methodology can be found elsewhere.86 Figure 6a shows in situ NMR spectra for a supercapacitor cell with YP50F activated carbon electrodes and a PEt4−BF4/ACN (1.5 M electrolyte), with 31P and 19F NMR allowing the study of the PEt4 cations and BF4 anions, respectively.52 Changes of in-pore chemical shift during charging arise from the changes of the carbon electronic structure as electrons are added or removed from the electrodes.86 Crucially, changes of in-pore peak intensities relate directly to changes of the in-pore ion populations. By fitting the spectra the number of in-pore cations and anions could be determined, revealing that the charge storage mechanism is inherently different depending on the polarization of the electrode (see Figure 6b).52 In the positive electrode, charging in this system occurs mainly by ion exchange (X = 0), with simultaneous counter-ion adsorption and co-ion desorption. In contrast, charging in the positive electrode occurs purely by counter-ion adsorption (X = 1). Despite the different charging mechanisms, the excess ionic charge in the carbon pores balances the electronic charge in both electrodes (Figure 6c). Interestingly, the charging mechanism was also found to be invariant to the concentration of the electrolyte, with the same mechanisms observed for concentrations of 1.5, 0.75, and 0.5 M. This mirrors observations on neat and diluted ionic liquids studied by EQCM83 and MD simulations39,40 above, and further confirms the idea that the solvent concentration does not necessarily dictate which charging mechanism operates, but simply modulates the absolute in-pore populations. In our study, the findings from the NMR measurements were corroborated by EQCM measurements, where a small negative mass change was observed in the positive electrode, as the slightly heavier cations were desorbed from the pores while the lighter anions were adsorbed. Meanwhile, a positive mass change was observed in the negative electrode, in agreement with the counter-ion adsorption mechanism revealed by NMR. The mass increase in this electrode was greater than that expected due to cation adsorption, and the additional mass was ascribed to solvent molecules, with an estimated cation solvation number of 5.4. 2H NMR experiments were also performed to study the deuterated acetonitrile solvent, though the large peak line widths precluded the quantification of the in-pore solvent. Future in situ NMR studies may allow quantitative measurement of in-pore ion solvation numbers, allowing comparison with the results from EQCM analysis. Further NMR measurements on the same activated carbon (YP50F) and a range of different electrolytes have revealed that the charging mechanism can vary significantly when different electrolyte ions are studied. Ex situ NMR studies of supercapacitors with the ionic liquid Pyr13−TFSI showed that in the positive electrode charging took place by ion exchange and counter-ion adsorption (X = 0.3), while in the negative electrode ion exchange and co-ion desorption were observed (X = −0.4).53 In this system the TFSI anions are more active in the charge storage than the Pyr13 cations. This charging mechanism is markedly different to that observed for the PEt4−BF4/ACN system above, where the BF4 anions played no significant role in charging in the negative electrode (see schematic in Figure 7). Similarly, in situ 19F NMR measurements on YP50F with Li−TFSI/ACN (1.5 M) and Na−TFSI/ACN (1.5 M) electrolytes revealed significant co-ion desorption processes in the negative electrode, though the cations were not studied.49 The exact origin of these different charging mechanisms is currently unclear, but ion packing and ion−carbon interaction energies are expected to be important. Co-ion (anion) desorption was

Figure 5. EQCM measurements of the charge storage mechanism of supercapacitors. Γ gives the moles of ions per surface area of quartz crystal, obtained from the measured electrode mass changes with the assumption of charge storage solely by counter-ion adsorption processes. Dashed lines showing theoretical Γ values for charging by pure counter-ion adsorption. (a) Measurements for YP-17 activated carbon and NEt4−BF4 in propylene carbonate (PC) solvent at two different concentrations. (b) Further measurements for YP-17 carbon supercapacitors with different electrolyte salts (0.1 M in PC): TEA+, TBA+, and TOA+ are tetraethylammonium, tetrabutylammonium, and tetraoctylammonium (NOc4), respectively. Adapted with permission from ref 81. Copyright 2010 American Chemical Society.

with ACN. Interestingly, the charging mechanism was largely unaffected by the presence of solvent, though additional mass changes were detected in the negative electrode, assigned to the solvent molecules that the EMI cations carried into the pores. Following the approach developed by Levi for aqueous electrolytes,84 a solvation number of 3.7 could be estimated for the EMI cations by assuming a purely ion-adsorption-driven charging mechanism, and assuming that the unaccounted mass was due to ACN solvent molecules. Given the assumptions here, it would be beneficial if solvation numbers could be verified by a second method such as NMR spectroscopy (see later). EQCM studies have advanced our understanding of supercapacitor charging mechanisms. However, a significant limitation of these studies is that a single parameter is measured (the electrode mass), which depends on the number of cations, the number of anions and the number of solvent molecules. Thus, the populations of the various in-pore species cannot be determined. This makes it difficult to fully quantify the charging mechanism and obtain X values, such that information from other experimental techniques is crucial. 3.4. Nuclear Magnetic Resonance (NMR) Spectroscopy. As discussed in section 2, NMR experiments allow the absolute quantification of ions inside the carbon nanopores, with the separate study of anions, cations, and solvent possible. Along with co-workers, we have developed the in situ NMR approach for studying supercapacitors, such that in-pore ion populations can be tracked at different charging potentials in working supercapacitor cells.49,52,85,86 These measurements are typically performed in constant voltage mode; i.e., a fixed voltage is applied to the cell, and NMR spectra are acquired after equilibration of the 5738

DOI: 10.1021/jacs.6b02115 J. Am. Chem. Soc. 2016, 138, 5731−5744

Perspective

Journal of the American Chemical Society

Figure 6. (a) In situ NMR measurements of a supercapacitor with YP50F electrodes and PEt4−BF4/ACN (1.5 M) electrolyte. (b) In-pore ion populations at different charging voltages, obtained by fitting of the spectra in (a). (c) Ionic charge stored by in-pore ions, obtained from the in-pore ion populations in (b), as well as electronic charge, obtained from electrochemistry data. Reprinted with permission from ref 52. Copyright 2015 Macmillan Publishers Ltd.

Figure 7. Schematics showing charge storage mechanisms determined by NMR experiments. For the same activated carbon (YP50F) electrodes, the charging mechanism differs depending on the electrolyte: (a) PEt4−BF4/ACN (1.5 M) and (b) Pyr13−TFSI ionic liquid. Calculations of X are also given based on the experimental data. Note in part (b) the X value given for the negative electrode is the experimentally measured one, whereas based on the number of ions in the schematic one obtains an X value of −0.3. This simplification maintains clarity in the schematic. See refs 52 and53 for the original studies.

also observed in the positive electrode of a supercapacitor with the electrolyte NBu4−BF4, suggesting that the large size of the

NBu4 cations impeded their adsorption into the carbon pores, necessitating desorption of anions to store charge.49 This finding 5739

DOI: 10.1021/jacs.6b02115 J. Am. Chem. Soc. 2016, 138, 5731−5744

Perspective

Journal of the American Chemical Society

the carbon nanopores during charging, with experiments on nonporous onion-like carbons revealing no significant changes of absorbance during charging and supporting the hypothesis. These experiments demonstrated that ions entered the carbon nanopores during charging, though the charging mechanism was not fully quantified. Further experiments with EMI−TFSI and nanoporous carbon nanofiber (CNF) electrodes revealed similar losses of intensity for anions and cations in the positive electrode,88 though the intensity loss was more significant for the anions than the cations, explaining how an excess of ionic charge could develop inside the carbon pores. Experiments were also performed on modified CNFs, which had been activated with KOH at 800 °C and contained more oxygen-containing functional groups. In contrast to the measurements on the untreated CNFs, an increase in intensity was observed for the cation absorbances, indicating that cations were desorbed from the nanopores of the CNFs in the positive electrode. Together with the observation of the loss of intensity of the anion absorbance (as anions entered the pores of the CNFs), it was shown that an ion-exchange mechanism was operating. The measurements intriguingly showed how the charging mechanism of supercapacitors is dependent on the surface chemistry of the carbon materials.88 IR spectroscopy experiments have revealed new insights into the charging mechanisms of supercapacitors. However, a key limitation of the IR approach is that the in-pore ions cannot be directly detected, and instead one must rely on measurements of the bulk electrolyte surrounding the carbon particles. It is therefore challenging to make fully quantitative studies of the supercapacitor charging mechanism, and X values cannot be readily determined. 3.6. Scattering Approaches. Scattering-based methods have also been applied to the study of supercapacitor charging mechanisms. In situ SANS experiments utilize the different scattering properties of the various atoms in the electrolyte and electrode. By measuring the scattered neutron intensity from a single (working) electrode at different charging potentials, qualitative changes of the in-pore ion populations can be obtained. For example, in a study of aqueous H2SO4 electrolytes in activated carbon fabric electrodes, changes of scattered intensity were dominated by the migration of the hydrogen-containing ions and solvent molecules.57 On this basis it was inferred that ion-exchange processes bring about charging in activated carbon fabric electrodes, for example, with H3O+ replacing HSO4− and SO42− in the negative electrode. Similar findings were obtained for devices with the organic electrolyte NEt4−BF4 in dACN (1 M).56 Here, the increase in scattered intensity in the negative electrode indicated that the strong H-containing scatterers NEt4+ were adsorbed into the carbon pores, while neutron-absorbing 10BF4− were desorbed from the pores, i.e., an ion-exchange mechanism. Converse effects were observed in the positive electrode, again rationalized by an ion-exchange mechanism. An interesting feature of small-angle scattering measurements is that adsorption/ desorption processes can be monitored as a function of the carbon pore size. Boukhalfa et al. found that the largest changes of scattered intensity were observed for the smallest pores (i.e., at large scattering wavevectors). These measurements have the potential to complement MD simulations, where the relative ion populations in pores of different geometries and sizes can be tracked.40 However, in the SANS measurements full quantification of the in-pore ionic species at different potentials was not possible as the scattered neutron intensity (a single measured quantity) depends on the population of anions, cations, and solvent molecules.

is consistent with EQCM studies on activated carbons, where the role of the anions (co-ions) in the negative electrode became more significant for electrolytes with larger cations.81 It should, however, be kept in mind that the EQCM measurements were performed in a dynamic mode, with the frequency response of the quartz crystal measured during cyclic voltammetry measurements. Such dynamic measurements with EQCM are more likely to reveal kinetic effects (e.g., effects due to the different diffusion rates of the anions and cations) than the in situ NMR measurements that were performed at equilibrium after charging to a specific voltage. Other research groups have also applied NMR methods to the study of supercapacitor charging mechanisms. Ex situ NMR measurements by Deschamps et al. on NEt4−BF4/ACN showed that ion-exchange processes operated in both the positive and negative electrodes of two different activated carbons.48 More recently, Luo et al. have applied the in situ NMR methodology to supercapacitors with aqueous electrolytes.54,55 They studied an activated carbon material with small pores that were inaccessible to the ions (Na+ and F−) in the absence of an applied potential.55 However, when a potential above 0.4 V was applied, the F− anions were able to enter the pores of the positive electrode, with charging proceeding via counter-ion adsorption. Following discharge to 0 V, some ions remained in the carbon pores, suggesting the presence of a hysteresis in the charging mechanism. Interestingly, changes of in-pore chemical shift at ∼0.8 V indicated that the solvation number of the ions decreased above this voltage. The strength of the NMR approach arises from the ability to separately observe and quantify in-pore ions, as well as to separately detect anions and cations. This allows the absolute quantification of supercapacitor charging mechanisms (and X values can be determined). The technique, however, is limited to ions with NMR active nuclei, and if many experiments are to be performed at different voltages, ions containing sensitive nuclei such as 1H, 19F, 31P, 7Li, 11B, and 23Na must generally be studied. Another limitation is that in some cases the in-pore resonances cannot be clearly resolved (e.g., due to broad peaks, or small ring current shifts), and methods should thus be developed to improve spectral resolution for challenging carbon and electrolyte systems. For example, we have shown how crosspolarization experiments can be used to edit the NMR spectrum to reveal peaks solely due to ions at the carbon surface (i.e., in-pore ions).44 3.5. Infrared (IR) Spectroscopy. Beyond EQCM and NMR, IR spectroscopy has also been demonstrated as a useful probe of supercapacitor charging mechanisms.87,88 Here, changes in intensity of the absorbances from bond vibrations in the electrolyte anions and cations are monitored during charging, allowing the behavior of the two ions to be tracked separately. The IR radiation is directed on the working electrode of a supercapacitor cell (which is clamped onto the surface of a diamond attenuated total reflectance (ATR) crystal), and the reflected signal intensity is recorded to obtain an absorbance spectrum. In these measurements the IR radiation penetrates to ∼1 μm depth in the working electrode. In experiments on titanium carbide-derived carbon (TiC-CDC) supercapacitors with ionic liquid (EMI−TFSI) electrolyte, approximately equal losses of intensity were observed for both the anion and cation absorbances during charging, suggesting that both ions penetrated deeper into the carbon particles, beyond the depth to which the IR photons were able to penetrate.87 This was explained by postulating that both anions and cations had entered 5740

DOI: 10.1021/jacs.6b02115 J. Am. Chem. Soc. 2016, 138, 5731−5744

Perspective

Journal of the American Chemical Society Table 3. Selected Charging Mechanisms from in situ Characterization Methods and Computer Simulationsa charging mechanism method

carbon, electrolyte

cycling mode

NMR NMR and EQCM NMR

YP50F, PEt4−BF4/ACN (1.5 M) YP50F, PEt4−BF4/ACN (0.75 M)

static static

YP50F, PEt4−BF4/ACN (0.5 M)

static

NMR

YP50F, Pyr13−TFSI ionic liquid

MD

positive electrode

negative electrode

ref

ion exchange (X = 0)

counter-ion adsorption (X = 1)

52

static

ion exchange and counter-ion adsorption (X = 0.3)

ion exchange and co-ion desorption (X = −0.4)

53

CDC, BMI−PF6 ionic liquid

static

ion exchange and co-ion desorption (X = −0.11)

39

MD

CDC, BMI−PF6/ACN (1.5 M)

static

ion exchange and counter-ion adsorption (X = 0.34) ion exchange and counter-ion adsorption (X = 0.57)

ion exchange and co-ion desorption (X = −0.02)

40

MD

slit pore (1.1 nm), EMI−TFSI ionic liquid

static

ion exchange and co-ion desorption (X = −0.36)

ion exchange and counter-ion adsorption (X = 0.45)

43

NMR

ACs, NEt4−BF4/ACN (1 M)

static

ion exchange (X = 0)

ion exchange (X = 0)

48

XRT XRT XRT

YP-80, CsCl (aq) (1 M) YP-80, KCl (aq) (1 M) YP-80, NaCl (aq) (1 M)

dynamic dynamic dynamic

ion exchange (X = 0)

ion exchange (X = 0)

58

EQCM EQCM

CDC, EMI−TFSI ionic liquid CDC, EMI−TFSI/ACN (1.5 M)

dynamic dynamic

low V, ion exchange; high V, counter-ion adsorption

counter-ion adsorption

83

EQCM EQCM

YP17, NEt4−BF4 (0.1 M) YP17, NOc4−BF4 (0.1 M)

dynamic dynamic

low V, ion exchange; high V, counter-ion adsorption

low V, ion exchange; high V, counter-ion adsorption as above, except ion exchange over wider V range

81

IR IR

CNFs, EMI−TFSI ionic liquid modified CNFs, EMI−TFSI ionic liquid

dynamic dynamic

involves counter-ion adsorption involves ion exchange

not studied not studied

88

a

X values are given where possible.

A recent study by Prehal et al. developed in situ X-ray transmission (XRT) and SAXS methods to provide deeper insights into the charging mechanisms of supercapacitors with simple aqueous electrolytes (CsCl, KCl, and NaCl, all 1 M concentration) and the activated carbon YP-80.58 SAXS measurements on evacuated and electrolyte-soaked carbon electrodes confirmed that the pores were filled with electrolyte at 0 V. Variations of the XRT intensity during electrochemical cycling, taken together with the electrochemical data, then allowed semiquantitative measurements of the changes of ion concentrations at different charging potentials. The analysis relies on the different X-ray attenuation coefficients of the cations and anions, and it is assumed that changes of the amount of solvent (water) have no significant effects on the XRT signal. This analysis showed that an ion-exchange mechanism (X = 0) occurs for all three electrolytes studied, with the charging mechanism independent of the choice of cation. In situ SAXS measurements were consistent with these findings and also indicated that the counter-ions became more closely associated with the pore surfaces as the potential was increased.58 Further work with scattering-based approaches stands to further our understanding of supercapacitor charging mechanisms.

contain a considerable number of electrolyte ions in the absence of an applied potential. Larger in-pore populations are observed for more concentrated electrolytes, with the largest populations observed for ionic liquids. Given the large in-pore ion population in the absence of an applied potential, different charging mechanisms can operate when a potential is applied: counterion adsorption, co-ion desorption, and ion exchange (and combinations of these). We have introduced the charging mechanism parameter, X, to allow a convenient comparison of different charging mechanisms. In situ characterization experiments (NMR, EQCM, IR, and scattering methods) and simulations have shown that supercapacitor charging does not generally take place by counter-ion adsorption alone, as is the traditional view, and ion exchange plays an important role in the charge storage process. As shown by Table 3, in a wide range of experimental and computational studies, charging almost always involved some degree of ion exchange (X ≠ 1), and only rarely was charging driven purely by counter-ion adsorption (X = 1). Experiments have shown that the exact charging mechanism depends on the polarization of the electrode, the choice of electrolyte ions, and the choice of electrode material (Table 3). Surprisingly, the solvent concentration does not appear to significantly influence the charging mechanism that operates (Table 3), and it appears the role of the solvent is less significant than previously thought. This suggests that ion confinement effects,40 screening of ions by the carbon surfaces,31 as well as the lack of overscreening39 dominate the capacitance increase

4. SUMMARY AND OUTLOOK In summary, NMR spectroscopy measurements and MD simulations have shown that the pores of the carbon electrodes 5741

DOI: 10.1021/jacs.6b02115 J. Am. Chem. Soc. 2016, 138, 5731−5744

Perspective

Journal of the American Chemical Society reported for nanoporous carbons,24,25 rather than ion desolvation, while different in-pore ion populations at 0 V and fundamentally different charging mechanisms may also play a role. The wide range of charging mechanisms (X values) presented in Table 3 indicate that factors such as ion-packing energies and ion−carbon interaction energies dictate which charging mechanism is the thermodynamic one. Under dynamic (fast) charging conditions the charging mechanism may differ, though studies contrasting thermodynamic and kinetic charging mechanisms have not yet been carried out. When studying Table 3 it should be kept in mind that only NMR and MD simulations allow the absolute quantification of the different electrolyte species at different potentials, and that various assumptions must be made when inferring charging mechanisms from other in situ methods (e.g., with EQCM changes of in-pore anion, cation, and solvent populations have to be determined from a single measured parameter, the electrode mass). Further in situ characterization of supercapacitors with a wide range of different electrolytes and carbon materials will help to fully elucidate the factors which control supercapacitor charging mechanisms. As discussed above, each in situ characterization method comes with its limitations, and these methods must be further developed such that they can be applied to a wide range of systems in a quantitative manner, and on a range of charging time scales. Beyond experiments, advanced theoretical methods that can properly take into account the electronic structure of the ions and carbon electrodes could provide new insights. If one can properly understand the factors that influence the charging mechanism, it should be possible to tailor the mechanism by choice of the correct electrolyte-electrode combination. In particular, it would be interesting to study the relative magnitudes of ion−ion and ion−carbon interaction energies within the carbon pores. Additionally, one could investigate packing effects by studying ions with different shapes and sizes within carbon pores of different geometries. It is clear that systematic studies in which the electrolyte and carbon structure are varied independently will allow the understanding to progress here. This leads to an important question: How does the charging mechanism af fect the performance of a supercapacitor? The mechanism should have a significant effect on the power that supercapacitors can offer, and tailoring the mechanism may allow us to improve the power performance of devices. Indeed, theoretical work has indicated that initially empty pores (ionophobic pores) should charge more quickly than initially filled pores (ionophilic pores),36 though this effect is yet to be realized experimentally. Similarly, we expect that ion adsorption, ion exchange, and ion desorption mechanisms from initially filled pores should each result in different device power performances, and work must be done to establish which mechanism is optimal for fast charging. Purely counter-ion adsorption processes might be expected to allow fast charging in a front-like manner, with net migration of ions into the interior of the carbon nanoporosity, while ion exchange requires ionic migration in opposite directions. At the same time, these different mechanisms will bring about changes of in-pore ionic density, and therefore packing during charging, which will also affect the charging rate. For example, counter-ion adsorption mechanisms will increase the number of ions inside the carbon pores, and recent theoretical work has suggested that more densely packed pores result in slower ionic diffusion.89 As well as ion packing effects, interactions of the different ions with charged carbon surfaces will also affect in-pore transport processes.78 Clearly, experimental measurements

and simulations of diffusion and migration processes in charged carbon nanopores represents another exciting area of research. In principle, the charging mechanism will affect the capacitance, and therefore the energy density, that can be achieved in supercapacitors. Under thermodynamic conditions, the charge storage mechanism which operates is the one that minimizes the increase in free energy associated with charging, thus minimizing the voltage increase per unit charge (i.e., maximizing the capacitance as C = Q/V). Kondrat and Kornyshev point out that counter-ion adsorption is disfavored here, as there is an entropic penalty for an ion entering a pore, and there are also unfavorable electrostatic (enthalpic) terms associated with the packing of ions of the same charge inside the carbon pores (though this is alleviated to some extent by charge screening from the pore walls).37 In principle, the ion-exchange mechanism reduces the enthalpic penalty associated with denser ion packing because the total in-pore density remains essentially constant during charging, while the entropic penalty associated with charging will also be reduced. This may help to explain the prevalence of ion exchange mechanisms revealed by in situ characterization methods detailed in this article (Table 3). Charging by co-ion desorption should minimize the enthalpic penalty due to interactions between like charges, while also increasing entropy, and should thus maximize the capacitance.37 That said, charging by purely co-ion desorption (X = −1) has not yet been observed, indicating that other factors beyond these simple arguments are important. It is clear that additional work must be done to understand the interplay between supercapacitor charging mechanisms and capacitance. Under kinetic control, the charging mechanism will depend on the relative rates of in-pore motion of the anions and cations, and we again stress the need for further experimental and theoretical studies of these effects. If the diffusion rates of the different in-pore ions can be controlled, then it should be possible to control the kinetic charging mechanism and thus improve the capacitance. Finally, we stress that the structural complexity17,18 of porous carbon electrodes poses a significant challenge as we aim to design enhanced supercapacitors. Ideally the electrode structure would be modified in a controlled way to study its performance, though for this to happen new tools must be developed to characterize and model amorphous carbon structures.15 Beyond activated carbons, more ordered carbon materials based on nanotubes, graphenes, and templated materials may serve as model systems to probe structure−property relationships.



AUTHOR INFORMATION

Corresponding Authors

*[email protected] *[email protected] Notes

The authors declare no competing financial interest. This is a review article, and therefore all data underlying this study is cited in references.



ACKNOWLEDGMENTS The authors acknowledge the Sims Scholarship Cambridge (A.C.F.), the School of the Physical Sciences of the University of Cambridge (via an Oppenheimer Research Fellowship, C.M.), EPSRC (via the Supergen consortium, A.C.F. and J.M.G.), and the EU ERC (via an Advanced Fellowship to C.P.G.) for funding. We thank Nicole Trease, Andrew Ilott, Phoebe Allan, Elizabeth Humphreys, Paul Bayley, Hao Wang, Patrice Simon, Wan-Yu 5742

DOI: 10.1021/jacs.6b02115 J. Am. Chem. Soc. 2016, 138, 5731−5744

Perspective

Journal of the American Chemical Society

(37) Kondrat, S.; Kornyshev, A. A. Nanoscale Horiz. 2016, 1, 45. (38) Palmer, J. C.; Llobet, A.; Yeon, S.-H.; Fischer, J. E.; Shi, Y.; Gogotsi, Y.; Gubbins, K. E. Carbon 2010, 48, 1116. (39) Merlet, C.; Rotenberg, B.; Madden, P. A.; Taberna, P.-L.; Simon, P.; Gogotsi, Y.; Salanne, M. Nat. Mater. 2012, 11, 306. (40) Merlet, C.; Péan, C.; Rotenberg, B.; Madden, P. A.; Daffos, B.; Taberna, P.-L.; Simon, P.; Salanne, M. Nat. Commun. 2013, 4, 2701. (41) Chaban, V. V.; Prezhdo, O. V. ACS Nano 2014, 8, 8190. (42) Shim, Y.; Kim, H. J. ACS Nano 2010, 4, 2345. (43) Xing, L.; Vatamanu, J.; Borodin, O.; Bedrov, D. J. Phys. Chem. Lett. 2013, 4, 132. (44) Forse, A. C.; Griffin, J. M.; Wang, H.; Trease, N. M.; Presser, V.; Gogotsi, Y.; Simon, P.; Grey, C. P. Phys. Chem. Chem. Phys. 2013, 15, 7722. (45) Merlet, C.; Forse, A. C.; Griffin, J. M.; Frenkel, D.; Grey, C. P. J. Chem. Phys. 2015, 142, 094701. (46) Xing, Y.-Z.; Luo, Z.-X.; Kleinhammes, A.; Wu, Y. Carbon 2014, 77, 1132. (47) Borchardt, L.; Oschatz, M.; Paasch, S.; Kaskel, S.; Brunner, E. Phys. Chem. Chem. Phys. 2013, 15, 15177. (48) Deschamps, M.; Gilbert, E.; Azais, P.; Raymundo-Piñero, E.; Ammar, M. R.; Simon, P.; Massiot, D.; Béguin, F. Nat. Mater. 2013, 12, 351. (49) Griffin, J. M.; Forse, A. C.; Wang, H.; Trease, N. M.; Simon, P.; Grey, C. P. Faraday Discuss. 2014, 176, 49. (50) Urita, K.; Ide, N.; Isobe, K.; Furukawa, H.; Moriguchi, I. ACS Nano 2014, 8, 3614. (51) Kim, T.; Ideta, K.; Jung, D.; Saito, K.; Park, J.; Rhee, C. K.; Miyawaki, J.; Mochida, I.; Yoon, S.-H. RSC Adv. 2014, 4, 16726. (52) Griffin, J. M.; Forse, A. C.; Tsai, W.-Y.; Taberna, P.-L.; Simon, P.; Grey, C. P. Nat. Mater. 2015, 14, 812. (53) Forse, A. C.; Griffin, J. M.; Merlet, C.; Bayley, P. M.; Wang, H.; Simon, P.; Grey, C. P. J. Am. Chem. Soc. 2015, 137, 7231. (54) Luo, Z.-X.; Xing, Y.-Z.; Ling, Y.-C.; Kleinhammes, A.; Wu, Y. Nat. Commun. 2015, 6, 6358. (55) Luo, Z.-X.; Xing, Y.-Z.; Liu, S.; Ling, Y.-C.; Kleinhammes, A.; Wu, Y. J. Phys. Chem. Lett. 2015, 6, 5022. (56) Boukhalfa, S.; Gordon, D.; He, L.; Melnichenko, Y. B.; Nitta, N.; Magasinski, A.; Yushin, G. ACS Nano 2014, 8, 2495. (57) Boukhalfa, S.; He, L.; Melnichenko, Y. B.; Yushin, G. Angew. Chem., Int. Ed. 2013, 52, 4618. (58) Prehal, C.; Weingarth, D.; Perre, E.; Lechner, R. T.; Amenitsch, H.; Paris, O.; Presser, V. Energy Environ. Sci. 2015, 8, 1725. (59) Trognko, L.; Lecante, P.; Ratel-Ramond, N.; Rozier, P.; Daffos, B.; Taberna, P.-L.; Simon, P. Mater. Renew. Sustain. Energy 2015, 4, 17. (60) Bañuelos, J. L.; Feng, G.; Fulvio, P. F.; Li, S.; Rother, G.; Dai, S.; Cummings, P. T.; Wesolowski, D. J. Chem. Mater. 2014, 26, 1144. (61) Dash, R.; Chmiola, J.; Yushin, G.; Gogotsi, Y.; Laudisio, G.; Singer, J.; Fischer, J.; Kucheyev, S. Carbon 2006, 44, 2489. (62) Vora, P.; Gopu, P.; Rosario-Canales, M.; Pérez, C.; Gogotsi, Y.; Santiago-Avilés, J.; Kikkawa, J. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 84, 155114. (63) Takashiro, J.; Kudo, Y.; Kaneko, S.; Takai, K.; Ishii, T.; Kyotani, T.; Enoki, T.; Kiguchi, M. Phys. Chem. Chem. Phys. 2014, 16, 7280. (64) Dyatkin, B.; Gogotsi, Y. Faraday Discuss. 2014, 172, 139. (65) Weingarth, D.; Zeiger, M.; Jäckel, N.; Aslan, M.; Feng, G.; Presser, V. Adv. Energy Mater. 2014, 4, 1400316. (66) Merlet, C.; Rotenberg, B.; Madden, P. A.; Salanne, M. Phys. Chem. Chem. Phys. 2013, 15, 15781. (67) Merlet, C.; Péan, C.; Rotenberg, B.; Madden, P. A.; Simon, P.; Salanne, M. J. Phys. Chem. Lett. 2013, 4, 264. (68) Siepmann, J. I.; Sprik, M. J. Chem. Phys. 1995, 102, 511. (69) Kornyshev, A. A. J. Phys. Chem. B 2007, 111, 5545. (70) Goduljan, A.; Juarez, F.; Mohammadzadeh, L.; Quaino, P.; Santos, E.; Schmickler, W. Electrochem. Commun. 2014, 45, 48. (71) Mohammadzadeh, L.; Goduljan, A.; Juarez, F.; Quaino, P.; Santos, E.; Schmickler, W. Electrochim. Acta 2015, 162, 11. (72) Mohammadzadeh, L.; Goduljan, A.; Juarez, F.; Quaino, P.; Santos, E.; Schmickler, W. ChemPhysChem 2016, 17, 78.

Tsai, Yury Gogotsi, Mathieu Salanne, Benjamin Rotenberg, Alexei Kornyshev, Svyatoslav Kondrat, and Volker Presser for collaboration, and for stimulating discussions and insights into supercapacitors over the course of our research on this subject.



REFERENCES

(1) Conte, M. Fuel Cells 2010, 10, 806. (2) Miller, J. R.; Simon, P. Science 2008, 321, 651. (3) Marsh, H.; Rodriguez-Reinoso, F. Activated Carbon, 1st ed.; Elsevier Sci. Ltd.: Amsterdam, 2006. (4) Presser, V.; Heon, M.; Gogotsi, Y. Adv. Funct. Mater. 2011, 21, 810. (5) Zhong, C.; Deng, Y.; Hu, W.; Qiao, J.; Zhang, L.; Zhang, J. Chem. Soc. Rev. 2015, 44, 7484. (6) Fic, K.; Lota, G.; Meller, M.; Frackowiak, E. Energy Environ. Sci. 2012, 5, 5842. (7) Dyatkin, B.; Presser, V.; Heon, M.; Lukatskaya, M. R.; Beidaghi, M.; Gogotsi, Y. ChemSusChem 2013, 6, 2269. (8) Fu, C.; Grant, P. S. ACS Sustainable Chem. Eng. 2015, 3, 2831. (9) Béguin, F.; Presser, V.; Balducci, A.; Frackowiak, E. Adv. Mater. 2014, 26, 2219 DOI: 10.1002/adma.201304137. (10) Brandt, A.; Pohlmann, S.; Varzi, A.; Balducci, A.; Passerini, S. MRS Bull. 2013, 38, 554. (11) Lewandowski, A.; Olejniczak, A.; Galinski, M.; Stepniak, I. J. Power Sources 2010, 195, 5814. (12) Petkov, V.; Difrancesco, R. G.; Billinge, S. J. L.; Acharya, M.; Foley, H. C. Philos. Mag. B 1999, 79, 1519. (13) Acharya, M.; Strano, M. S.; Mathews, J. P.; Billinge, S. J. L.; Petkov, V.; Subramoney, S.; Foley, H. C. Philos. Mag. B 1999, 79, 1499. (14) Forse, A. C.; Griffin, J. M.; Presser, V.; Gogotsi, Y.; Grey, C. P. J. Phys. Chem. C 2014, 118, 7508. (15) Forse, A. C.; Merlet, C.; Allan, P. K.; Humphreys, E. K.; Griffin, J. M.; Aslan, M.; Zeiger, M.; Presser, V.; Gogotsi, Y.; Grey, C. P. Chem. Mater. 2015, 27, 6848. (16) Harris, P. J. F.; Liu, Z.; Suenaga, K. J. Phys.: Condens. Matter 2008, 20, 362201. (17) Harris, P. J. F. J. Mater. Sci. 2013, 48, 565. (18) Harris, P. Crit. Rev. Solid State Mater. Sci. 2005, 30, 235. (19) Franklin, R. E. Acta Crystallogr. 1951, 4, 253. (20) Franklin, R. E. Proc. R. Soc. London, Ser. A 1951, 209, 196. (21) Lastoskie, C.; Gubbins, K. E.; Quirke, N. J. Phys. Chem. 1993, 97, 4786. (22) Dombrowski, R. J.; Hyduke, D. R.; Lastoskie, C. M. Langmuir 2000, 16, 5041. (23) Neimark, A. V.; Lin, Y.; Ravikovitch, P. I.; Thommes, M. Carbon 2009, 47, 1617. (24) Chmiola, J.; Yushin, G.; Gogotsi, Y.; Portet, C.; Simon, P.; Taberna, P.-L. Science 2006, 313, 1760. (25) Raymundo-Piñero, E.; Kierzek, K.; Machnikowski, J.; Béguin, F. Carbon 2006, 44, 2498. (26) Chmiola, J.; Largeot, C.; Taberna, P.-L.; Simon, P.; Gogotsi, Y. Angew. Chem., Int. Ed. 2008, 47, 3392. (27) Largeot, C.; Portet, C.; Chmiola, J.; Taberna, P.-L.; Gogotsi, Y.; Simon, P. J. Am. Chem. Soc. 2008, 130, 2730. (28) Feng, G.; Cummings, P. T. J. Phys. Chem. Lett. 2011, 2, 2859. (29) Jiang, D.; Jin, Z.; Wu, J. Nano Lett. 2011, 11, 5373. (30) Wu, P.; Huang, J.; Meunier, V.; Sumpter, B. G.; Qiao, R. ACS Nano 2011, 5, 9044. (31) Kondrat, S.; Kornyshev, A. J. Phys.: Condens. Matter 2011, 23, 022201. (32) Kondrat, S.; Georgi, N.; Fedorov, M. V.; Kornyshev, A. A. Phys. Chem. Chem. Phys. 2011, 13, 11359. (33) Centeno, T. A.; Sereda, O.; Stoeckli, F. Phys. Chem. Chem. Phys. 2011, 13, 12403. (34) Garcia-Gomez, A.; Moreno-Fernandez, G.; Lobato, B.; Centeno, T. A. Phys. Chem. Chem. Phys. 2015, 17, 15687. (35) Kondrat, S.; Kornyshev, A. J. Phys. Chem. C 2013, 117, 12399. (36) Kondrat, S.; Wu, P.; Qiao, R.; Kornyshev, A. A. Nat. Mater. 2014, 13, 387. 5743

DOI: 10.1021/jacs.6b02115 J. Am. Chem. Soc. 2016, 138, 5731−5744

Perspective

Journal of the American Chemical Society (73) Varanasi, S. R.; Farmahini, A. H.; Bhatia, S. K. J. Phys. Chem. C 2015, 119, 28809. (74) Pean, C.; Rotenberg, B.; Simon, P.; Salanne, M. Electrochim. Acta 2016, DOI: 10.1016/j.electacta.2016.02.106. (75) Fedorov, M. V.; Kornyshev, A. A. Electrochim. Acta 2008, 53, 6835. (76) Merlet, C.; Salanne, M.; Rotenberg, B.; Madden, P. A. J. Phys. Chem. C 2011, 115, 16613. (77) Péan, C.; Merlet, C.; Rotenberg, B.; Madden, P. A.; Taberna, P.L.; Daffos, B.; Salanne, M.; Simon, P. ACS Nano 2014, 8, 1576. (78) Pean, C.; Daffos, B.; Rotenberg, B.; Levitz, P.; Haefele, M.; Taberna, P.-L.; Simon, P.; Salanne, M. J. Am. Chem. Soc. 2015, 137, 12627. (79) Levi, M. D.; Salitra, G.; Levy, N.; Aurbach, D.; Maier, J. Nat. Mater. 2009, 8, 872. (80) Sauerbrey, G. Eur. Phys. J. A 1959, 155, 206. (81) Levi, M. D.; Levy, N.; Sigalov, S.; Salitra, G.; Aurbach, D.; Maier, J. J. Am. Chem. Soc. 2010, 132, 13220. (82) Merlet, C. Ph.D. thesis, Université Pierre et Marie Curie, 2013; pp 127−131. (83) Tsai, W.-Y.; Taberna, P.-L.; Simon, P. J. Am. Chem. Soc. 2014, 136, 8722. (84) Levi, M. D.; Sigalov, S.; Salitra, G.; Elazari, R.; Aurbach, D. J. Phys. Chem. Lett. 2011, 2, 120. (85) Wang, H.; Köster, T. K.-J.; Trease, N. M.; Ségalini, J.; Taberna, P.L.; Simon, P.; Gogotsi, Y.; Grey, C. P. J. Am. Chem. Soc. 2011, 133, 19270. (86) Wang, H.; Forse, A. C.; Griffin, J. M.; Trease, N. M.; Trognko, L.; Taberna, P.-L.; Simon, P.; Grey, C. P. J. Am. Chem. Soc. 2013, 135, 18968. (87) Richey, F. W.; Dyatkin, B.; Gogotsi, Y.; Elabd, Y. A. J. Am. Chem. Soc. 2013, 135, 12818. (88) Richey, F. W.; Tran, C.; Kalra, V.; Elabd, Y. A. J. Phys. Chem. C 2014, 118, 21846. (89) He, Y.; Qiao, R.; Vatamanu, J.; Borodin, O.; Bedrov, D.; Huang, J.; Sumpter, B. G. J. Phys. Chem. Lett. 2016, 7, 36.

5744

DOI: 10.1021/jacs.6b02115 J. Am. Chem. Soc. 2016, 138, 5731−5744