Complementary Effects of Pore Accessibility and Decoordination on

Dec 7, 2015 - For a more comprehensive list of citations to this article, users are encouraged .... Agency, the U.S. Department of Justice, and the Mi...
0 downloads 0 Views 1MB Size
Subscriber access provided by CMU Libraries - http://library.cmich.edu

Article

Complementary Effects of Pore Accessibility and De-Coordination on the Capacitance of Nanoporous Carbon Electrochemical Supercapacitors Srinivasa Rao Varanasi, Amir Hajiahmadi Farmahini, and Suresh Kumar Bhatia J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.5b10712 • Publication Date (Web): 07 Dec 2015 Downloaded from http://pubs.acs.org on December 8, 2015

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Complementary Effects of Pore Accessibility and De-coordination on the Capacitance of Nanoporous Carbon Electrochemical Supercapacitors

Srinivasa Rao Varanasi, Amir H. Farmahini and Suresh K. Bhatia* School of Chemical Engineering, The University of Queensland, Brisbane, QLD 4072, Australia

Abstract We illustrate here the interplay of de-coordination and accessible pore volume in nanosized supercapacitors, using constant voltage Gibbs ensemble based grand canonical Monte Carlo simulations for three different microporous carbon electrodes of known atomistic structure, and 1-ethyl-3methylimidazolium boron tetrafluoride (EMI-BF4) as electrolyte. We demonstrate that the counter-ion coordination number decreases with pore size and this trend is similar for the electrodes considered, despite their different structures; suggesting that the pore shape is less important to this relation, at least for the carbons examined here. It is seen that ions with low coordination and/or completely de-coordinated ions induce maximum charge, while those with higher coordination induce less, in accordance with recent MD simulation results which demonstrate that ions in high degree of confinement (DOC) induce more charge than those in low DOC. Our results indicate that electrodes with different pore volumes can exhibit similar capacitances by balancing accessibility and de-coordination effects. Thus, similar capacitance may be obtained for electrodes having low pore volume, but which can adsorb a small amount of high charge inducers (de-coordinated ions) by virtue of having suitable pore size, and those having high pore volume and adsorbing many more low charge inducers (more highly coordinated ions). Keywords: Electrochemical supercapacitor; Ion de-coordination; Ionic liquid; Confinement effects; Nanoporous carbon. *To whom correspondence may be addressed. Tel.: +61 7 3365 4263. Email: [email protected].

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1. Introduction Carbon based electrochemical supercapacitors are an important class of charge storage systems that comprise two high surface area porous carbons as electrodes and an electrolyte, either solvent based or solvent free, commonly a room temperature ionic liquid (RTIL). These systems possess high cycle life and high power, since the storage mechanism is based on the formation of electric double layers (EDLs), involving physical changes in the electrode with no chemical reaction, making the charge/discharge mechanism highly reversible and rapid (having time scale of seconds or fractions of a second)1. Supercapacitors based on RTILs have received much attention in recent years, since they are environmentally friendly and support high operating voltages, despite showing lower capacitance and high equivalent series resistance (ESR) as compared to their solvent based counter-parts (usually having acetonitrile or propylene-carbonate as solvents). However, due to their non-volatile nature and higher thermal stability, RTILs are considered greener and safer options for electrochemical supercapacitor applications2-4.

Nanoporous carbide derived carbons (CDCs), such as that from titanium carbide (denoted TiC-DC) or silicon carbide (SiC-DC), and activated carbon (AC), have been demonstrated to be suitable as electrode materials in supercapacitor applications5-10. An efficient supercapacitor demands that an electrode material possess high surface area and pores that are commensurate with the size of the electrolyte ions11-14. The enhanced capacitance observed in experiments with electrodes having sub-nanometer pores is attributed to partial or complete desolvation of counter-ions adsorbed inside the pores11-12, 15. Recent studies have shown that room temperature ionic liquids (such as those based on imidazolium cations) used in combination with CDCs, carbon nanotubes and activated carbons exhibit good 2

ACS Paragon Plus Environment

Page 2 of 36

Page 3 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

electrochemical performance4, 10, 16-18. There are numerous studies in the literature reporting capacitances as high as 140 F/g based on the use of different types of porous carbons in combination with RTILs6, 15, 19-22. Greater capacitances have been achieved when the average pore size is commensurate with the size of the ion (desolvated ion in case of solvent electrolytes), indicating a possible correlation between ion size and pore size11-12,

23-24

.

However, there is a need to consider the pore size distribution (PSD), as broader PSDs may lead to invariance of capacitance over a wide range of mean pore sizes similar to that observed by Centeno et al25. To optimize materials for higher capacitance, it is important to understand the factors that are influencing the degree of ionic charge separation between electrodes and the mechanism leading to this. In situ NMR spectroscopic techniques have been developed to monitor the changes in the supercapacitor system during charge/discharge processes, promising to offer better insights into the microscopic mechanism of supercapacitor performance26-30.

Nevertheless, experimental techniques frequently have

limitations in uncovering microscopic mechanisms, as they provide information on average behavior but cannot yet track individual configurations or trajectories of ions and molecules, and theoretical analysis remains a powerful tool for understanding supercapacitor performance.

Molecular modelling techniques such as molecular dynamics (MD) and Monte Carlo (MC) simulations have emerged as promising tools to elucidate the molecular underpinnings of the electrochemical behavior of supercapacitor materials31-36. Theoretical and molecular simulation studies by Kornyshev and co-workers37-43 have revealed several interesting features of supercapacitors combining carbon electrodes with ionic liquid electrolytes. For instance, charging dynamics in ionophobic pores is relatively fast compared to ionophilic 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ones, thereby providing an opportunity to optimize the electrodes for enhanced energy as well as power densities38. These authors have shown that electrodes with monodispersed pore size distributions (PSDs) deliver higher energy densities than those with polydispersed PSDs40. It has also been demonstrated that pores with atomically rough surfaces exhibit substantial capacitance enhancement as compared to those with smooth surfaces44. MD simulation studies reveal that local structure of the carbon electrode has large influence on the charging of carbon atoms in the electrode, possibly by the distinct organization of electrolyte. Thus, contrasting electrochemical performance of different TiC-DC carbon structures having similar pore size distribution has been predicted31. The poor performance of TiC-DC prepared at 1200 °C as compared to TiC-DC prepared at 950 °C is due to the presence of small graphitic domains that facilitate over- screening effects, such as those occurring in graphitic pores31. More disordered carbon structures deliver better electrochemical performance than the less disordered ones possessing graphitic domains45. MD simulation studies by Merlet et al46 demonstrate that the ions adsorbed in high degree confinements induce counter-charge on the electrode atoms more effectively than those in the low degree confinements, and such ions lose most of their solvation/coordination shell neighbors, leading to partial/complete desolvation or de-coordination. Hence, it is important to consider the microstructure of the electrode pores in addition to surface area, average pore size and pore size distribution, while optimizing electrodes for superior electrochemical performance.

In the present study, we perform Gibbs ensemble based grand canonical Monte Carlo simulations for three distinct microporous carbon structures, activated carbon (ACF-15), SiCDC and TiC-CDC, as electrodes and a RTIL, 1-ethyl-3methylimidazolium boron tetrafluoride (EMI-BF4), as electrolyte. Our simulations and analysis of the microscopic structure and 4

ACS Paragon Plus Environment

Page 4 of 36

Page 5 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

energetics of the ions in the carbons reveal comparable values of capacitance in these distinct structures, despite significant differences in pore volume, and demonstrate this behavior to arise from to two complementary factors, namely, pore accessibility and de-coordination in nano-confinement. Our results explain how these two factors complement each other in order to deliver similar capacitances in the electrodes examined.

2. Computational details: Electrode models: We have considered three different microporous carbon structures independently as electrodes, namely, activated carbon fibre (ACF-15), SiC-DC and TiC-DC. Atomistic models of these materials have been developed using hybrid reverse Monte Carlo (HRMC) simulation as reported elsewhere47-49. We have taken three unit cells of ACF-15 along the zaxis (total of 3486 atoms), one unit cell of SiC-DC (2998 atoms) and one unit cell of TiC-DC (2570 atoms) as electrodes in their respective systems. The choice of three unit cells has been made for ACF-15 so as to ensure roughly equal number of carbon atoms (i.e. approximately equal mass) for all electrodes. The electrodes have been kept 10 nm apart along the z-axis while periodic boundary conditions have been imposed along x and y axes.

Electrolyte model and intermolecular potential: The ionic liquid has been modelled using a coarse grained four site model proposed by Roy et al50, in which each imidazolium cation is represented by three connected sites whereas BF4 anions are represented as single sites. Each site in the electrolyte ions is treated as a LennardJones (LJ) sphere with a point charge placed on it. The potential parameters have been taken from the work of Merlet et al51 (see supplementary information). The charge on each carbon 5

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 36

atom of the electrodes is modelled as a Gaussian distribution centered on the atom, following Vatamanu et al.52 and earlier studies53-54. The total interaction energy of the system, comprising two carbon electrodes and ionic liquid electrolyte is given by52:

Φ to ta l = Φ L J + Φ E

(1)

Where Φ LJ ΦE =



i , j ,i ≠ j

 σ ij = ∑ 4 ε ij    r  i , j ,i ≠ j   ij qi q j 4 π ε r ε 0 rij



1 2



  

i , j ,i ≠ j

    q i q j e rfc ( γ rij ) + 4 π ε r ε 0 rij

12

 σ ij −  rij 

  

6

(2)

∑ i

q i2 γδ ( i , G )

π

(3)

Here rij is distance between two atoms (or ions) i, j = 1, 2, qi is partial charge on atom (or ion) i, σij, εij are Lennard-Jones (LJ) parameters for the i-j interaction, εo is permittivity of vacuum, εr is relative permittivity and 1/γ is width of the Gaussian distribution. The first term in eqn. (1) is the contribution from all the LJ interactions (eqn. (2)), excluding those between electrode atoms. The second term (eqn. (3)) is the contribution from electrostatic interactions. The second and third terms in eqn. (3) arise due to the Gaussian nature of the electrode charges. The delta function in eqn. (3) is zero for point charge and one for Gaussian charge distribution. The first term in eqn. (3) has been calculated using the modified Wolf method5558

. The convergence parameter in the Wolf method has been adjusted between 1-2 nm-1,

depending upon the simulation cut-off to ensure the electrostatic energies to converge to those obtained by Ewald summation approach. The cut-off length for all the interactions has been set to be 14.75 Å, 20 Å and 17.5 Å for ACF-15, SiC-DC and TiC-DC respectively, corresponding to half of the smallest dimension of the simulation cell. The specific width of the Gaussians (1/γ) is chosen to be 0.2 Å, ensuring the short range approximation to the second term in eqn. (3) is valid within the cut-off.52. LJ interactions between the atoms of the electrodes have not been considered, as the intramolecular degrees of freedom of the porous 6

ACS Paragon Plus Environment

Page 7 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

electrode structures are assumed to be frozen. Hence there are only electrostatic interactions between the fluctuating electrode charges. While the electrochemical dilatometry studies on carbon electrodes have shown that there can be charge-induced microscopic structural changes in the electrode due to pumping of electrons in and out, and also due to insertion of electrolyte ions inside smaller pores, inclusion of such effects is beyond the scope of the present work. Image interactions have not been treated explicitly in our simulations. The cross interaction parameters between dissimilar atoms have been calculated using standard Lorentz-Berthelot mixing rules. All the simulations have been performed in vacuum (ϵr = 1). Simulation procedure Pore characterization of electrode materials: Pore size distributions (PSDs) of the carbon models have been determined using the spherical probe geometric approximation technique proposed by Gelb and Gubbins59,60. The method makes use of the Metropolis Monte Carlo algorithm to insert random points into the unit cell. The center to center distance of a successfully inserted probe molecule with any atom of the network must be always larger than ߪ௦௢௟௜ௗି௣௥௢௕௘ to ensure no overlap occurs. Once statistically adequate sampling of the phase space was achieved, the porosity is calculated based on the number of accepted moves divided by total number of insertion attempts. For the PSD, the acceptance of a successful trial insertion is followed by a random walk away from its original point which is repeated in many different random directions. This way the maximum diameter of a spherical probe which simultaneously encompasses the original insertion point and does not overlap any network atom is determined. The pore sizes are calculated from the Connolly surface which is defined to be the surface traced by the exterior of the probe particle while rolling over the pore wall. 7

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 36

Constant voltage Gibbs ensemble based GCMC simulations: After characterizing the pore structure of the electrode materials, we have performed Gibbs ensemble based constant voltage grand canonical Monte Carlo simulations on porous carbon electrodes, following the recent algorithm of Punnathanam35, as

implemented by us

elsewhere36. The grand canonical partition function at constant voltage is given by:

Ξ (µ , Φ ,V , T ) =





∑ σ

∑ eβ µ

A = −∞

[ N + Φ θt ]

Q ( N ,θ t ,V , T )

(4)

N =0

where Q ( N ,θ t ,V , T ) =

N



N1 A = 0

N



Q AC ( N 1 A , N1C , N , θ t , V , T )

(5)

N1 C = 0

where QAC ( N1 A , N1C , N ,θt ,V , T ) =

(qV )2 N ..... e− βU ds 2 N N1A ! N1C ! N2 A ! N2C ! ∫ ∫

(6)

and the probability density distribution is given by:

f (s2N , N1A ,N1C ,N2A ,N2C ,θt ) ∝

(qV)2N βµN βΦθt −βU e e e N1A!N1C!N2A!N2C!

(7)

Here, µ, Φ, V and T are chemical potential, applied voltage, volume and temperature respectively. θt is the total charge on the electrode, β is (kBT)-1 and q is eβµ. The number of cations and anions in each of the electrodes are N1C, N2C and N1A, N2A respectively, where 1, 2 are the indices of the electrodes. Imposing the following condition ensures charge neutrality in the overall system:

N1 A + N 2 A = N1C + N 2C = N

(8)

There are six different Monte Carlo moves (formulated using Eqn. 6) performed during the simulation, namely, translation, insertion, deletion, charge plus ion transfer (CPIT), ion 8

ACS Paragon Plus Environment

Page 9 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

transfer (IT) and charge transfer (CT). In a translation move, an ion is chosen randomly from a randomly selected electrode and moved by a small distance (determined by the maximum allowed acceptance, 35%), and rotated (in case of EMI+) by a small angle with respect to a randomly selected axis passing through the center of mass of the selected ion. In an insertion (or deletion) move, an ion pair is created in a randomly chosen electrode at random positions (similarly an ion pair is selected for deletion in a deletion move). In a charge plus ion transfer move (CPIT), an ion (cation or anion) is chosen randomly from a randomly selected electrode and transferred to the other electrode, and

a charge equivalent to the ion’s charge is

added/subtracted to/from the origin/destination electrode respectively. By doing this, the CPIT move preserves the charge neutrality of each of the electrodes. In an ion transfer (IT) move, an ion (cation or anion) is selected randomly from a randomly chosen electrode and transferred to the other electrode. Unlike the CPIT move, no charge compensation on the respective electrodes is performed, and therefore this move does not preserve charge neutrality of the individual electrodes, although overall charge neutrality is preserved following eqn (8). In a charge transfer move (CT), a small amount of charge is transferred between electrodes. The transferred charge is shared only between randomly chosen 500 or half the number of atoms, whichever is smaller, in the destination electrode, which ensures non-uniform charge distribution among the atoms in the destination electrode. All these moves are accepted or rejected according to Metropolis scheme, and the acceptance probability for any move is:

 f new  acc PMC  move = min 1,  f old 

9

ACS Paragon Plus Environment

(9)

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

where fnew and fold represent probability density distributions (see Eqn. (6)) for the states after and before the moves respectively. Further details of the Monte Carlo moves and their acceptance criterion can be found elsewhere35-36. The simulations are initiated starting with empty electrodes with a zero charge on each carbon atom. The voltage has been set at 2V. The system is equilibrated for 20,000 Monte Carlo (MC) cycles, comprising about 5 million MC moves. The number of Monte Carlo moves in each cycle is equal to the total number of ions in the system, while the minimum number of MC moves in each MC cycle is set to be 100. A production run of about 60,000 MC cycles is performed during which the coordinates of the ions and charges on the electrode atoms are stored at the end of each MC cycle for subsequent analysis. All the simulations have been performed at a temperature of 300 K and a bulk activity of 0.001nm-3. At this activity, bulk simulations on the ionic liquid model, provided a number density of 3.61 ion pairs/nm3, which is very close to what is expected from the density values given by Merlet et al51 (3.86 ion pairs/nm3). All Monte Carlo moves have been given equal priority. 3. Results and discussions Pore characterization: Prior to performing Gibbs ensemble based grand canonical Monte Carlo simulations, we have characterized the pore structure of the three porous carbon electrodes using the spherical probe geometric approximation technique, with three different hard sphere probes of diameters 2.64 Å, 5.84 Å and 4.51 Å representing He, EMI+ and BF4respectively. Figures 1a, 1b and 1c illustrate the resulting pore size distributions. The PSDs indicate that the three structures possess a broad distribution between 0-8Å, for all the three probes. It is seen that the specific pore volume of ACF-15 is relatively higher than that of other two structures, up to the maximum pore size (10.5 Å). In this pore range (0-10.5 Å), the 10

ACS Paragon Plus Environment

Page 10 of 36

Page 11 of 36

specific pore volume of SiC-DC is somewhat lesser than that of ACF-15. However, it possesses pores larger than 10.5 Å and up to 13.5 Å, absent in ACF-15 and TiC-DC. TiC-DC possesses two narrow peaks between 8.6 and 10.4 Å, lying at 9 Å and 10 Å respectively. It is seen that the specific pore volume is smaller for TiC-DC in comparison to the other two carbons, though the latter two structures may compete with each other depending on the probe size. The values of specific pore volume for the three structures with different probes have been tabulated in Table 1. The specific pore volume of ACF-15 is slightly higher than that of SiC-DC. The difference between these two is 0.055, 0.011 and 0.04 cc/g/Å respectively with Helium, EMI+ and BF4- probes. Possibly, this variation between different probes is due to the fact that He and BF4- could probe the smaller pores around 5 Å and EMI+ could not. Most of the pore volume in SiC-DC is contributed from the pores in the range 9-13 Å while in ACF-15 0.20

it is mostly from

(a) He

ACF-15 SiC-DC TiC-DC

pores between 8-

0.15

PSD (cc/g/Å)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

10.5 Å. 0.10

0.05

0.00 1

2

3

4

5

6

7

8

9

10

11

12

13

pore size (Å)

11

ACS Paragon Plus Environment

14

15

The Journal of Physical Chemistry

0.20

(b) EMI+

ACF-15 SiC-DC TiC-DC

PSD (cc/g/Å)

0.15

0.10

0.05

0.00 4

5

6

7

8

9

10

11

12

13

14

15

pore size (Å) 0.20

(c) BF4-

ACF-15 SiC-DC TiC-DC

0.15

PSD (cc/g/Å)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 36

0.10

0.05

0.00 4

5

6

7

8

9

10

11

12

13

14

15

pore size (Å)

Figure 1: Pore size distributions of ACF-15, SiC-DC and TiC-DC with hard sphere probes of size comparable to (a) He, (b) EMI+ and (c) BF4-. We also note that the cation is represented as a single sphere of effective size of 5.84Å, though it is a combination of three connected sites in the actual coarse grained model (see supplementary information). Hence, the cation may see a slightly different PSD than that shown in Figure 1b, and the difference in the accessibility of ACF-15 and SiC-DC in the constant voltage simulations may be slightly different from that predicted from the PSD and specific pore volumes. However, the specific pore volumes (Table 1) using three different 12

ACS Paragon Plus Environment

Page 13 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

probes indicate that the apparent accessibility of TiC-DC is much smaller than that of ACF15 and SiC-DC, and this order may not alter during the constant voltage simulations.

Table 1: Specific pore volume of ACF-15, SiC-DC and TiC-DC calculated using spherical probe approximation technique with three different probes.

Carbon

Specific pore volume (cc/g) using different size probes He EMI+ BF4(2.64 Å) (5.84 Å) (4.51 Å)

ACF-15

0.615

0.459

0.54

SiC-DC

0.560

0.448

0.50

TiC-DC

0.335

0.190

0.256

Charge storage capability and mechanism: Following characterization of the pore structure of the three electrode materials, we have conducted Gibbs ensemble based grand canonical ensemble simulations on three different systems comprising ACF-15, SiC-DC and TiC-DC as electrodes with EMI-BF4 electrolyte. These three carbon structures are distinctly different with respect to carbon density, PSD and specific pore volume (c.f. Figure 1 and Table 1). The gravimetric and volumetric charge densities per electrode are tabulated in Table 2. The integral capacitance has been calculated from: Cint = σ / ∆V

(10)

where σ is the charge density accumulated on one of the electrodes, which can be either gravimetric or volumetric. We note that it is important to consider both volume and mass of the electrodes while optimizing the electrode materials for supercapacitor applications.

13

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 36

Hence, the gravimetric and volumetric charge densities are calculated using mass and volume of the dry electrodes, with ‘< >’ representing the ensemble average. The denominator, ∆V is the applied potential difference between the two electrodes. The gravimetric charge density values for ACF-15 and SiC-DC are close to each other while TiC-DC is not significantly different. The gravimetric and volumetric capacitances calculated from eqn. (10) are also given in Table 2. The capacitance values (particularly volumetric) for the activated carbon are in agreement with those reported in the literature15 (100 F/g, 50 F/cc), while those for CDCs are lower than values for similar carbons

6, 20-21

(160 F/g, 85 F/cc) , possibly due to the

narrow pore structure of our electrodes. The values of gravimetric and volumetric capacitances are quite similar to each other in all the three electrodes, as their densities are close to 1.0 g/cm3. The competing results obtained for TiC-DC, despite its lower specific volume, are counter-intuitive and quite intriguing, requiring a more detailed probe into the charging mechanism at the microscopic level.

Table 2: Values of the gravimetric charge density and capacitance, and volumetric charge density and capacitance for electrodes based on ACF-15, SiC-DC and TiC-DC. Carbon

Gravimetric charge density (C/g)

Gravimetric capacitance (F/g)

Volumetric charge density (C/cc)

Volumetric capacitance (F/cc)

ACF-15

117.63

58.815

103.16

51.58

SiC-DC

117.80

58.90

110.496

55.248

TiC-DC

101.38

50.69

120.64

60.32

We note that the charging of a porous electrode is not just filling the electrode with the electrolyte but it is mainly the (ionic) charge separation between the electrodes when an external voltage is applied. The amount of charge induced in an electrode is proportional to 14

ACS Paragon Plus Environment

Page 15 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

the difference between the number of counter-ions and co-ions in the electrode and does not depend upon the total amount of the electrolyte electro-sorbed into the electrode. The values of total number density and the difference between counter-ion to co-ion number densities are tabulated in Table 3. The number density values indicate that relatively high amount of electrolyte has been electro-sorbed into SiC-DC electrodes as compared to other two electrode materials; in agreement with the comparison of specific pore volumes estimated using neutral probes (c.f. Table 1). Figures 2a, 2b and 2c illustrate snapshots of cathodes of ACF-15, SiC-DC and TiC-DC filled with EMI-BF4 electrolyte respectively. Visual inspection of the snapshots confirms that TiC-DC is filled with relatively lower amount of electrolyte as compared to the other two. Despite lower total number density, the difference in number densities of counter-ions and co-ions in the electrodes of TiC-DC are significantly higher than in the other two, indicating greater charge separation, leading to values of charge density and capacitance comparable to those of ACF-15 and SiC-DC.

Table 3: Values of total number density (number of ions per unit volume) and the difference between the number densities of counter-ion and co-ion in the electrodes, for ACF-15, SiCDC, TiC-DC electrodes at 2V. Carbon

Total number density, ρ (nm-3)

ρcounter-ion - ρco-ion (nm-3)

ACF-15

2.348

0.728

SiC-DC

2.625

0.797

TiC-DC

1.796

0.840

Having determined the relative accessibility of the electrodes from the specific pore volumes estimated by the spherical probe geometric approximation technique using neutral probes, and the overall number density of the electrolyte electro-sorbed into the electrodes during constant voltage simulations, we have also analyzed the wetting behavior of the electrodes with respect to the counter-ion as well as the co-ion. 15

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

In order to compare the wettability of porous carbon structures at the set voltage, we have defined a parameter called wetting length, which is the distance between the ion and carbon atom, within which the number of carbon atoms around a given ion type are enumerated. The fraction of carbon atoms within the wetting length from the given ion has been evaluated and its variation with wetting length is illustrated in Figure 3, for the counter-ion in cathode (cation). Figure 4 depicts that for co-ion in cathode (anion). It is seen that about 75% carbon atoms in SiC-DC have been exposed to counter-ions, while it is 70 and 60% for ACF-15 and TiC-DC respectively, at a wetting length of 0.5 nm.

Figure 2: Snapshots of ion-electrode configurations in (a) ACF-15 (b) SiC-DC and (c) TiCDC cathodes. Yellow and blue spheres represent EMI+ (three sites) cation and BF4 (one site) anions respectively. The electrodes are shown in cyan coloured bonded representation. Some of the inaccessible regions have been encircled. The fraction of carbon atoms exposed to counter-ion is always greater for SiC-DC than the other two, at all wetting lengths between 0.35 and 0.5 nm. This indicates that the carbon 16

ACS Paragon Plus Environment

Page 16 of 36

Page 17 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

atoms in TiC-DC are less accessible to counter-ions of the RTIL than those in ACF-15 and SiC-DC, as evident from the PSDs depicted in Figure 1. The increase in fraction of carbons with wetting length (c.f. Figures 3 and 4) is less steep in case of TiC-DC, indicating that there are more bulky regions in TiC-DC which are not accessible to the electrolyte (marked by circles in Figure 2c), though there exist less bulkier inaccessible regions in ACF-15 (see Figure 2a) and SiC-DC (see Figure 2b). The co-ions in SiC-DC could approach carbons more closely than in TiC-DC, which might happen in larger pores of size between 1 nm and 1.3 nm (Figure 4). The fraction of carbons within 0.5 nm of wetting length from co-ions in the cathode is about 5% and 4% in SiC-DC and ACF-15 respectively, whereas it is only 2% in TiC-DC. The co-ions will induce counter charge on the nearby carbon atoms, leading to reduced charge density. From these observations, we can infer that the larger pores, particularly those above 1 nm, will accommodate more co-ions besides counter-ions, leading to reduced charge separation in ACF-15 and SiC-DC electrodes. The fraction of carbons at the shortest approachable distance (see insets in Figs. 3 and 4) is comparatively higher in TiC-DC than in the other two, indicating tighter confinement in TiC-DC and this will be reflected by the ion-electrode energy, particularly the short-range part (LJ), which will be discussed later. Besides the relative number of counter-ions and co-ions inside the electrode, the structure of coordination shells of co-ions around counter-ions will provide further insights into the charge storage mechanism. The coordination shell structure is completely modified when the ionic liquid is subjected to confinement, depending upon pore size, shape and pore size distribution.

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1.0 ACF-15 SiC-DC TiC-DC

0.8

cation in cathode

0.6

0.4

fraction of carbons

fraction of carbons

0.0040 0.0035 0.0030 0.0025 0.0020 0.0015 0.30

0.2

wetting length (nm)

0.0 (a)

0.30

0.35

0.40

0.45

0.50

wetting length (nm)

Figure 3: Fraction of carbon atoms wetted by the cations in the ACF-15, SiC-DC and TiCDC cathodes, as a function of wetting length.

0.06 ACF-15 SiC-DC TiC-DC

0.05

anion in cathode

0.04

0.03

0.02

fraction of carbons

fraction of carbons

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

4e-7

3e-7

2e-7

1e-7 0.33995

0.34000

0.01

0.34005

0.34010

0.34015

0.34020

wetting length (nm )

0.00

0.36

0.38

0.40

0.42

0.44

0.46

0.48

0.50

0.52

wetting length (nm)

Figure 4: Fraction of carbon atoms wetted by the anions in the ACF-15, SiC-DC and TiCDC cathodes, as a function of wetting length.

18

ACS Paragon Plus Environment

Page 18 of 36

Page 19 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

In bulk, each cation is surrounded by about 7 anions in the first coordination shell (coordination number). However, the coordination number will be less than 7 when the ions are subjected to confinement31, 46. The coordination number obtained by integrating the first peak of the radial distribution function up to the first minimum will be an average quantity, but does not provide adequate information about the counter-ion and co-ion coordination inside the electrodes. Hence, we have calculated the distribution of coordination number of co-ions around each counter-ion within 0.63 nm (the position of first minimum of RDF between cation and anion, see Figure S5 of supplementary information) for all the three structures and in bulk as well. Figure 5 depicts the distribution of coordination number for anions around cations in cathodes of all the three porous structures at 2V, in TiC-DC with no charge on carbon atoms and in bulk. The coordination number exhibits a sharp peak at 7 in the bulk and has a broad distribution between 0 and 6 in ACF-15 and SiC-DC. It is evident that the coordination shell is effectively modified depending on the level of confinement. Unpaired counter-ions (coordination number zero) will polarize the electrode more than the paired ones. The most probable coordination number in ACF-15 and SiC-DC is 2 while it is 1 for TiC-CDC. The maximum possible coordination number is 5 and 6 in ACF-15 and SiCDC respectively, while it is 3 in TiC-DC at 2V. For SiC-DC, the lower coordination numbers (0 to 2) are marginally less probable compared to ACF-15; in contrast the higher coordination numbers (3 to 6) are marginally less probable in ACF-15. This explains why ACF-15 attains charge density nearly equal to that of SiC-DC, despite possessing lower specific pore volume compared to SiC-DC. The coordination number distribution in neutral TiC-DC is quite distinct from that in bulk, shifted towards lower coordination numbers (Figure 5), showing that the cations lose most of their coordination shell partners inside the nanoporous TiC-DC structure even in the absence of external field. This indicates that the nanoporous TiC-DC 19

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

structure modifies the coordination shell significantly and facilitates effective charge separation process between the TiC-DC electrodes. The completely de-coordinated (coordination number = 0) and single coordinated (coordination number = 1) counter-ions are more probable in TiC-DC than in the other two carbons. These observations from our Monte Carlo simulations are consistent with what has previously been reported from MD simulation studies on RTILs in TiC-DC electrodes46. Counter-ions that are completely de-coordinated will polarize the electrodes more effectively than the others. Complete de-coordination will occur only when the ion size is commensurate with the pore size, and such pores seem to be more probable in TiC-DC than in ACF-15 and SiC-DC. The larger coordination shells can appear only in the larger pores, particularly those larger than 1 nm, such as those present particularly in SiC-DC. Figure 6 demonstrates the average coordination number of counterions as a function of pore size, determined by the spherical probe geometric approximation technique (see simulation procedure). It shows that the average coordination number is between 1 and 2 for pores of size between 0.7 and 0.9 nm, and 3 in 1 nm pores in all the three structures. Higher coordination numbers are found in SiC-DC in pores larger than 1 nm. The variation of coordination number appears to be similar for all three carbons, implying that the pore shape is less important to this relation at least for the carbons considered here. However, this cannot be considered generic until a systematic study is performed on a variety nanoporous carbons with varying complexity in the pore structure. Besides probing the significant de-coordination in sub-nanometer pores, it is also important to examine the amount of counter-charge induced in the electrode atoms locally by these differently coordinated counter-ions. Figure 7 illustrates the variation of charge induced by each ion on the electrode atoms with coordination number in both cathode and anode of the three carbon structures. Maximum charge, up to 1.0e, is induced by completely de20

ACS Paragon Plus Environment

Page 20 of 36

Page 21 of 36

coordinated ions, whereas the ions corresponding to higher coordination numbers induce very little or almost no charge on the electrode atoms.

0.6

ACF-15 SiC-DC TiC-DC TiC-DC (no charge) Bulk

0.5

fraction of ions

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

0.4

cation in cathode

0.3

0.2

0.1

0.0 -1

0

1

2

3

4

5

6

7

8

9

10

coordination number

Figure 5: Distribution of anion coordination number around cations in the ACF-15, SiC-DC and TiC-DC cathodes at 2V, in neutral TiC-DC and in the bulk. This analysis provides a basis to explain why these three electrodes exhibit almost similar capacitances, despite having distinct pore volumes and pore accessibilities. TiC-DC attains most of its capacitance by electro-sorbing high charge inducers (de-coordinated ions for instance with coordination number 0 or 1) whereas ACF-15 and SiC-DC attain the similar capacitance by electro-sorbing more of somewhat highly coordinated ions including the decoordinated ones. Thus, the shortcoming of weaker de-coordination effects in ACF-15 and SiC-DC is overcome by higher accessibility, and vice-versa in TiC-DC. The asymmetry in the charge induced by counter-ions in cathode (cation) and anode (anion) stems from the

21

ACS Paragon Plus Environment

The Journal of Physical Chemistry

asymmetry in the structure of cation and anion (see Figure S4 of supplementary information for coordination number distribution of anions in anode).

6

ACF-15 SiC-DC TiC-DC

5

coordination number

cation in cathode

4

3

2

1

0 6

7

8

9

10

11

12

13

14

pore size (Å)

Figure 6: Average coordination number of cation in cathode as a function of the size of its local pore. 1.5 ACF-15 (+) SiC-DC (+) TiC-DC (+) ACF-15 (-) SiC-DC (-) TiC-DC (-)

1.0

charge induced per ion (e)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 36

anode

0.5

0.0

-0.5

cathode -1.0

-1.5 0

2

4

6

8

coordination number of counter-ion

Figure 7: Average charge induced by adsorbed counter-ions on their local electrode pore atoms as a function of their coordination number. Here ‘+’ and ‘-’ represent cathode and anode respectively. 22

ACS Paragon Plus Environment

Page 23 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

After analyzing the structure of the porous electrodes and coordination shells of counter-ions electro-sorbed inside the electrode, we have also analyzed energies contributed from LJ as well as electrostatic interactions. The distribution of LJ energy per ion contributed only from the ions of the RTIL inside the cathode pores for all the three porous materials are shown in Figure 8. The most probable ion-ion LJ energy per ion in TiC-CDC is around 4.0 kJ/mol, while it is -6.8 kJ/mole and -7.3 kJ/mol in ACF-15 and SiC-DC respectively. These values are much smaller in magnitude than that in the bulk (-19.5 kJ/mol),. The ion-ion LJ energy depends upon the number of immediate neighbors, essentially the coordination shells of counter-ions in this context. It indicates that the extent of clustering of ions (number of nearest neighbors) is much higher in ACF-15 and SiC-DC than in TiC-DC. The extent of clustering in ACF-15 is somewhat lesser than in SiC-DC. This difference can be attributed to the presence of larger pores, for instance, those larger than 1 nm in ACF-15 and SiC-DC. The ion-ion LJ energy in neutral (i.e. uncharged) TiC-DC is close to that in TiC-DC at 2V, and more repulsive than that of bulk, indicating that the neutral structure significantly reduces the number of coordination shell neighbors even in the absence of an external field. However, it is also crucial to understand how these ionic clusters interact with the porous structures under which they are subjected to confinement. The distributions of LJ interaction energy per ion, between ions and the electrode atoms are shown in Figure 9, in all the three structures. It shows that the LJ energy between ions and electrode atoms in TiC-DC is more negative than that in ACF-15 and SiC-DC which indicates that the ions in TiC-DC are subjected to stronger confinement than in ACF-15 and SiC-DC, also indicated by the relatively higher fraction of carbons found in TiC-DC at the shortest wetting length (see insets of Figure 3 and 4). The confinement in ACF-15 appears to be somewhat stronger than in SiC-DC, which can again be 23

ACS Paragon Plus Environment

The Journal of Physical Chemistry

attributed to their distinct pore size distributions. The ion-electrode LJ energy in TiC-DC cathode at 2V is somewhat more attractive than in neutral TiC-DC because of the changes in the electrolyte configurations and charge separation due to the external field.

ACF-15 SiC-DC TiC-DC TiC-DC (no charge) Bulk

normalised distribution function

2.8 2.4 2.0 1.6 1.2 0.8 0.4 0.0 -21

-20

-19

-9

-8

-7

-6

-5

-4

-3

-2

-1

ELJ (kJ/mol)

Figure 8: Distribution of ion-ion LJ energy per ion in the ACF-15, SiC-DC and TiC-DC electrodes at 2V, in uncharged TiC-DC and in the bulk.

1.8 ACF-15 SiC-DC TiC-DC TiC-DC (no charge)

1.6

normalised distribution function

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 36

1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 -35

-34

-33

-32

-31

-30

-29

-28

-27

-26

-25

ELJIL-ELEC (kJ/mol)

Figure 9: Distribution of ion-electrode LJ energy per ion in the ACF-15, SiC-DC, TiC-DC electrodes at 2V, and in uncharged TiC-CDC.

24

ACS Paragon Plus Environment

Page 25 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

In supercapacitor systems, electrostatic energy is the main component of energy which can indicate the extent of ionic separation between the electrodes. Figure 10 illustrates the ion-ion electrostatic energy per ion for the three electrode materials. It is evident that the electrostatic energy per ion is more repulsive under confinement than in bulk, and is more repulsive in TiC-DC than in ACF-15 and SiC-DC, and is somewhat repulsive in ACF-15 than in SiC-DC. More repulsive ion-ion electrostatic energy indicates the domination of completely decoordinated or less coordinated counter-ions and increased level of confinement, indicated by rapid increase in the radial distribution function between counter-ion and co-ion before the first peak (illustrated in Figure S5 of supplementary information). The ion-ion electrostatic energy per ion in neutral TiC-DC is more repulsive than that in the bulk (Figure 10 (inset)), indicating a significant reduction of number of coordination partners of a given ion (Figure 5). The repulsive nature of electrostatic energy in TiC-DC cathode at 2V with respect to that in neutral TiC-DC is attributed to the charge separation between the electrodes (electrodes attracting more ions of counter charge).

However, the neutral pores of TiC-DC are

instrumental in breaking the coordination shell of the ions, aiding the external field to induce charge density more effectively than in other structures.

25

ACS Paragon Plus Environment

1.2

normalised distribution function

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

distribution function

The Journal of Physical Chemistry

1.0

0.8

Page 26 of 36

2.0 TiC-DC (no charge) Bulk

1.5

1.0

0.5

0.0

0.6

-122 -120 -118

0.4

ACF-15 SiC-DC TiC-DC

-72 -69

ECoul (kJ/mol)

0.2

-66

-63

0.0 -48

-44

-40

-28

-26

-24

-22

-20

-18

ECoul (kJ/mol)

Figure 10: Distribution of ion-ion Coulomb energy per ion in the ACF-15, SiC-DC, TiC-DC electrodes at 2V, uncharged TiC-DC and bulk (inset).

4. Conclusions We have performed Gibbs ensemble based grand canonical ensemble Monte Carlo simulations on ACF-15, SiC-DC and TiC-DC microporous carbons as electrodes in combination with an RTIL. Our results demonstrate that all the three structures exhibit similar gravimetric and volumetric capacitances, despite wide contrast in their pore structure. The counter-ion coordination with co-ions is greatly reduced under confinement (decoordination effect), similar to desolvation of ions in solvent based electrolytes in high confinements46.

The distribution of coordination number of counter-ions (EMI+ in the

cathode) has a broad distribution in the range 0-6 in ACF-15 and SiC-DC (compared to average coordination number of 7 in the bulk), whereas this distribution is quite narrow in TiC-DC. TiC-DC exhibits capacitance close to that of ACF-15 and SiC-DC, despite 26

ACS Paragon Plus Environment

Page 27 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

possessing almost half the amount of specific pore volume as compared to ACF-15 and SiCDC. The Lennard-Jones part of the energy between electrolyte ions is lower and that between electrolyte and electrode is higher in TiC-DC than in ACF-15 and SiC-DC, despite relatively small amount of electro-sorption in TiC-DC; this is indicative of stronger confinement experienced by the electrolyte in TiC-DC. The contribution of Coulomb energy to ion-ion interaction energy is more repulsive in TiC-DC which indicates that the counter-ions are mostly de-coordinated in TiC-DC than in ACF-15 and SiC-DC. The coordination numbers of cations are greatly reduced even in the neutral pores of uncharged TiC-DC, facilitating stronger de-coordination effects when an external voltage is applied. Completely decoordinated counter-ions induce higher counter-charge in the local electrode atoms than those with one or more co-ion neighbors, consistent with the degree of confinement effect on counter-ion coordination number and subsequent effect on local carbon atoms, observed using MD simulations by Merlet et al46. TiC-DC delivers similar capacitance as ACF-15 and SiC-DC, complementing the shortage of accessibility by electro-sorbing high charge inducing counter-ions such as those with coordination number 0 or 1 while ACF-15 and SiC-DC complement the shortage of de-coordinated counter-ions by providing accessibility to many low charge inducers (somewhat high coordinated counter-ions). Our results imply that for excellent electrochemical performance, the electrode should possess a porous structure that has maximum accessibility to the electrolyte as well as pores of size that facilitate decoordination of counter-ions. 5. Acknowledgements This research has been supported by a grant (DP150101824) from the Australian Research Council under the Discovery Scheme. We thank Dr. J. C. Palmer, University of Houston, for providing the HRMC-derived structure of TiC-DC. 27

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ASSOCIATED CONTENT Supporting Information TiC-CDC 800 model, interaction model of electrolyte, coordination number distribution of anion in anode, and radial distribution function between cation-anion in cathode. The Supporting Information is available free of charge on the ACS Publications website at DOI: .

28

ACS Paragon Plus Environment

Page 28 of 36

Page 29 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

References 1.

Conway, B. E., Electrochemical Supercapacitors: Scientific Fundamentals and Technological Applications. Plenum Press: New York, 1999.

2.

Beguin, F.; Presser, V.; Balducci, A.; Frackowiak, E., Carbons and Electrolytes for Advanced Supercapacitors. Adv. Mater. 2014, 26 (14), 2219-2251.

3.

Ghandi, K., A Review of Ionic Liquids, Their Limits and Applications. Green and Sustainable Chemistry 2014, 04 (01), 44-53.

4.

Wei, D.; Ng, T. W., Application of Novel Room Temperature Ionic Liquids In Flexible Supercapacitors. Electrochem. Commun. 2009, 11 (10), 1996-1999.

5.

Su, D. S.; Schlogl, R., Nanostructured Carbon and Carbon Nanocomposites for Electrochemical Energy Storage Applications. Chem. Sus. Chem. 2010, 3 (2), 136-68.

6.

Presser, V.; Heon, M.; Gogotsi, Y., Carbide-Derived Carbons - From Porous Networks to Nanotubes and Graphene. Adv. Funct. Mater. 2011, 21 (5), 810-833.

7.

Li, X.; Wei, B., Supercapacitors Based on Nanostructured Carbon. Nano Energy 2013, 2 (2), 159-173.

8.

Ji, H.; Zhao, X.; Qiao, Z.; Jung, J.; Zhu, Y.; Lu, Y.; Zhang, L. L.; MacDonald, A. H.; Ruoff, R. S., Capacitance of Carbon-Based Electrical Double-Layer Capacitors. Nat. Commun. 2014, 5, 3317.

9.

Guo, Y.-G.; Hu, J.-S.; Wan, L.-J., Nanostructured Materials for Electrochemical Energy Conversion and Storage Devices. Adv. Mater. 2008, 20 (15), 2878-2887.

10.

Candelaria, S. L.; Shao, Y.; Zhou, W.; Li, X.; Xiao, J.; Zhang, J.-G.; Wang, Y.; Liu, J.; Li, J.; Cao, G., Nanostructured Carbon for Energy Storage and Conversion. Nano Energy 2012, 1 (2), 195-220.

29

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

11.

Chmiola, J.; Yushin, G.; Gogotsi, Y.; Portet, C.; Simon, P.; Taberna, P. L., Anomalous Increase in Carbon Capacitance at Pore Sizes Less Than 1 Nanometer. Science 2006, 313 (5794), 1760-1763.

12.

Chmiola, J.; Largeot, C.; Taberna, P. L.; Simon, P.; Gogotsi, Y., Desolvation of Ions in Subnanometer Pores and its Effect on Capacitance And Double-Layer Theory. Angew. Chem. Int. Edit. 2008, 47 (18), 3392-3395.

13.

Huang, J.; Sumpter, B. G.; Meunier, V., A Universal Model for Nanoporous Carbon Supercapacitors Applicable to Diverse Pore Regimes, Carbon Materials, And Electrolytes. Chemistry 2008, 14 (22), 6614-26.

14.

Jiang, D. E.; Jin, Z.; Wu, J., Oscillation of Capacitance Inside Nanopores. Nano Lett. 2011, 11 (12), 5373-5377.

15.

Largeot, C.; Portet, C.; Chmiola, J.; Taberna, P. L.; Gogotsi, Y.; Simon, P., Relation Between the Ion Size and Pore Size for an Electric Double-Layer Capacitor. J. Am. Chem. Soc. 2008, 130 (9), 2730-2731.

16.

Brandt, A.; Pohlmann, S.; Varzi, A.; Balducci, A.; Passerini, S., Ionic Liquids in Supercapacitors. MRS Bull. 2013, 38 (7), 554-559.

17.

Feng, G.; Li, S.; Presser, V.; Cummings, P. T., Molecular Insights into Carbon Supercapacitors Based on Room-Temperature Ionic Liquids. J. Phys. Chem. Lett. 2013, 4 (19), 3367-3376.

18.

Timperman, L.; Beguin, F.; Frackowiak, E.; Anouti, M., Comparative Study of Two Protic Ionic Liquids as Electrolyte for Electrical Double-Layer Capacitors. J. Electrochem. Soc. 2013, 161 (3), A228-A238.

19.

Simon, P.; Gogotsi, Y., Capacitive Energy Storage in Nanostructured CarbonElectrolyte Systems. Acc. Chem. Res. 2013, 46 (5), 1094-1103. 30

ACS Paragon Plus Environment

Page 30 of 36

Page 31 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

20.

Simon, P.; Gogotsi, Y., Materials for Electrochemical Capacitors. Nat. Mater. 2008, 7 (11), 845-854.

21.

Perez, C. R.; Yeon, S. H.; Segalini, J.; Presser, V.; Taberna, P. L.; Simon, P.; Gogotsi, Y., Structure and Electrochemical Performance of Carbide-Derived Carbon Nanopowders. Adv. Funct. Mater. 2013, 23 (8), 1081-1089.

22.

Tsai, W. Y.; Gao, P. C.; Daffos, B.; Taberna, P. L.; Perez, C. R.; Gogotsi, Y.; Favier, F.; Simon, P., Ordered Mesoporous Silicon Carbide-Derived Carbon for High-Power Supercapacitors. Electrochem. Commun. 2013, 34, 109-112.

23.

Feng, G.; Cummings, P. T., Supercapacitor Capacitance Exhibits Oscillatory Behavior as a Function of Nanopore Size. J. Phys. Chem. Lett. 2011, 2 (22), 2859-2864.

24.

Kalluri, R. K.; Biener, M. M.; Suss, M. E.; Merrill, M. D.; Stadermann, M.; Santiago, J. G.; Baumann, T. F.; Biener, J.; Striolo, A., Unraveling the Potential and Pore-Size Dependent Capacitance of Slit-Shaped Graphitic Carbon Pores in Aqueous Electrolytes. Phys. Chem. Chem. Phys. 2013, 15 (7), 2309-2320.

25.

Centeno, T. A.; Sereda, O.; Stoeckli, F., Capacitance in Carbon Pores of 0.7 to 15 nm: a Regular Pattern. Phys. Chem. Chem. Phys. 2011, 13 (27), 12403-12406.

26.

Wang, H.; Koster, T. K. J.; Trease, N. M.; Segalini, J.; Taberna, P. L.; Simon, P.; Gogotsi, Y.; Grey, C. P., Real-Time NMR Studies of Electrochemical Double-Layer Capacitors. J. Am. Chem. Soc. 2011, 133 (48), 19270-19273.

27.

Wang, H.; Forse, A. C.; Griffin, J. M.; Trease, N. M.; Trognko, L.; Taberna, P. L.; Simon, P.; Grey, C. P., In Situ NMR Spectroscopy of Supercapacitors: Insight into the Charge Storage Mechanism. J. Am. Chem. Soc. 2013, 135 (50), 18968-18980.

31

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

28.

Griffin, J. M.; Forse, A. C.; Wang, H.; Trease, N. M.; Taberna, P. L.; Simon, P.; Grey, C. P., Ion Counting in Supercapacitor Electrodes Using NMR Spectroscopy. Farad. Disc. 2014, 176, 49-68.

29.

Forse, A. C.; Griffin, J. M.; Merlet, C.; Bayley, P. M.; Wang, H.; Simon, P.; Grey, C. P., NMR Study of Ion Dynamics and Charge Storage in Ionic Liquid Supercapacitors. J. Am. Chem. Soc. 2015, 137 (22), 7231-7242.

30.

Forse, A. C.; Griffin, J. M.; Wang, H.; Trease, N. M.; Presser, V.; Gogotsi, Y.; Simon, P.; Grey, C. P., Nuclear Magnetic Resonance Study of Ion Adsorption on Microporous Carbide-Derived Carbon. Phys. Chem. Chem. Phys. 2013, 15 (20), 7722-7730.

31.

Merlet, C.; Rotenberg, B.; Madden, P. A.; Taberna, P. L.; Simon, P.; Gogotsi, Y.; Salanne, M., On the Molecular Origin of Supercapacitance in Nanoporous Carbon Electrodes. Nat. Mater. 2012, 11 (4), 306-310.

32.

Vatamanu, J.; Hu, Z.; Bedrov, D.; Perez, C.; Gogotsi, Y., Increasing Energy Storage in Electrochemical Capacitors with Ionic Liquid Electrolytes and Nanostructured Carbon Electrodes. J. Phys. Chem. Lett. 2013, 4 (17), 2829-2837.

33.

Wu, P.; Huang, J. S.; Meunier, V.; Sumpter, B. G.; Qiao, R., Complex Capacitance Scaling in Ionic Liquids-Filled Nanopores. ACS Nano 2011, 5 (11), 9044-9051.

34.

Gorniak, R.; Lamperski, S., Investigation of the Electrical Double Layer with a Graphene Electrode by the Grand Canonical Monte Carlo Simulation. J. Phys. Chem. C 2014, 118 (6), 3156-3161.

35.

Punnathanam, S. N., A Gibbs-Ensemble Based Technique for Monte Carlo Simulation of Electric Double Layer Capacitors (EDLC) at Constant Voltage. J. Chem. Phys. 2014, 140,174110.

32

ACS Paragon Plus Environment

Page 32 of 36

Page 33 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

36.

Varanasi, S. R.; Bhatia, S. K., Capacitance Optimization in Nanoscale Electrochemical Supercapacitors. J. Phys. Chem. C 2015, 119 (31), 17573-17584.

37.

Fedorov, M. V.; Kornyshev, A. A., Ionic Liquids at Electrified Interfaces. Chem. Rev. 2014, 114 (5), 2978-3036.

38.

Kondrat, S.; Wu, P.; Qiao, R.; Kornyshev, A. A., Accelerating Charging Dynamics in Subnanometre Pores (vol 13, pg 387, 2014). Nat. Mater. 2014, 13, 387-393.

39.

Lee, A. A.; Kondrat, S.; Oshanin, G.; Kornyshev, A. A., Charging Dynamics of Supercapacitors with Narrow Cylindrical Nanopores. Nano Tech. 2014, 25, 315401.

40.

Kondrat, S.; Perez, C. R.; Presser, V.; Gogotsi, Y.; Kornyshev, A. A., Effect of Pore Size and its Dispersity on the Energy Storage in Nanoporous Supercapacitors. Ener. Environ. Sci. 2012, 5 (4), 6474-6479.

41.

Kondrat, S.; Kornyshev, A., Superionic State in Double-Layer Capacitors with Nanoporous Electrodes J. Phys. Cond. Mat. 2011, 23, 022201.

42.

Kondrat, S.; Georgi, N.; Fedorov, M. V.; Kornyshev, A. A., A Superionic State in Nano-Porous Double-Layer Capacitors: Insights from Monte Carlo Simulations. Phys. Chem. Chem. Phys. 2011, 13 (23), 11359-11366.

43.

Lee, A. A.; Kondrat, S.; Kornyshev, A. A., Single-File Charge Storage in Conducting Nanopores. Phys. Rev. Lett. 2014, 113, 048701.

44.

Vatamanu, J.; Vatamanu, M.; Bedrov, D., Non-Faradaic Energy Storage by Room Temperature Ionic Liquids in Nanoporous Electrodes. ACS Nano 2015, 9 (6), 59996017.

45.

Deschamps, M.; Gilbert, E.; Azais, P.; Raymundo-Pinero, E.; Ammar, M. R.; Simon, P.; Massiot, D.; Beguin, F., Exploring Electrolyte Organization in Supercapacitor Electrodes with Solid-State NMR. Nat. Mater. 2013, 12 (4), 351-358. 33

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

46.

Merlet, C.; Pean, C.; Rotenberg, B.; Madden, P. A.; Daffos, B.; Taberna, P. L.; Simon, P.; Salanne, M., Highly Confined Ions Store Charge More Efficiently in Supercapacitors. Nat. Commun. 2013, 4, 2701.

47.

Nguyen, T. X.; Cohaut, N.; Bae, J. S.; Bhatia, S. K., New Method for Atomistic Modeling of the Microstructure of Activated Carbons using Hybrid Reverse Monte Carlo Simulation. Langmuir 2008, 24 (15), 7912-7922.

48.

Farmahini, A. H.; Opletal, G.; Bhatia, S. K., Structural Modelling of Silicon CarbideDerived Nanoporous Carbon by Hybrid Reverse Monte Carlo Simulation. J. Phys. Chem. C 2013, 117 (27), 14081-14094.

49.

Palmer, J. C.; Jain, S. J.; Gubbins, K. E.; Cohaut, N.; Fischer, J. E.; Dash, R. K.; Gogotsi, Y., Hybrid Reverse Monte Carlo Simulations of Microporous Carbons. In Characterisation of Porous Solids VIII, Royal Society of Chemistry Special Publications,: 2009, 318, 56-63.

50.

Roy, D.; Maroncelli, M., An Improved Four-Site Ionic Liquid Model. J. Phys. Chem. B 2010, 114 (39), 12629-12631.

51.

Merlet, C.; Salanne, M.; Rotenberg, B., New Coarse-Grained Models of Imidazolium Ionic Liquids for Bulk and Interfacial Molecular Simulations. J. Phys. Chem. C 2012, 116 (14), 7687-7693.

52.

Vatamanu, J.; Borodin, O.; Smith, G. D., Molecular Dynamics Simulations of Atomically Flat and Nanoporous Electrodes with A Molten Salt Electrolyte. Phys. Chem. Chem. Phys. 2010, 12 (1), 170-182.

53.

Reed, S. K.; Lanning, O. J.; Madden, P. A., Electrochemical Interface Between an Ionic Liquid and A Model Metallic Electrode. J. Chem. Phys. 2007, 126 (8).

34

ACS Paragon Plus Environment

Page 34 of 36

Page 35 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

54.

Siepmann, J. I.; Sprik, M., Influence of Surface-Topology and Electrostatic Potential on Water Electrode Systems. J. Chem. Phys. 1995, 102 (1), 511-524.

55.

Wolf, D.; Keblinski, P.; Phillpot, S. R.; Eggebrecht, J., Exact Method for the Simulation of Coulombic Systems by Spherically Truncated, Pairwise R(-1) Summation. J. Chem. Phys. 1999, 110 (17), 8254-8282.

56.

Ma, Y.; Garofalini, S. H., Molecular Dynamics Simulations of Beta-SiC Using Both Fixed Charge and Variable Charge Models. J. Chem. Phys. 2008, 128, 084505.

57.

Fennell, C. J.; Gezelter, J. D., Is the Ewald Summation Still Necessary? Pairwise Alternatives to the Accepted Standard for Long-Range Electrostatics. J. Chem. Phys. 2006, 124, 234104.

58.

Ma, Y.; Garofalini, S. H., Application of the Wolf Damped Coulomb Method to Simulations of SiC. J. Chem. Phys. 2005, 122, 094508.

59.

Gelb, L. D.; Gubbins, K. E., Characterization of Porous Glasses: Simulation Models, Adsorption Isotherms, and the Brunauer-Emmett-Teller Analysis Method. Langmuir 1998, 14 (8), 2097-2111.

60.

Gelb, L. D.; Gubbins, K. E., Pore Size Distributions in Porous Glasses: A Computer Simulation Study. Langmuir 1999, 15 (2), 305-308.

35

ACS Paragon Plus Environment

The Journal of Physical Chemistry

TOC Graphic 6

ACF-15 SiC-DC TiC-DC

5

coordination number

cation in cathode

4

3

2

1

0 6

7

8

9

10

11

12

13

14

pore size (Å)

1.5 ACF-15 (+) SiC-DC (+) TiC-DC (+) ACF-15 (-) SiC-DC (-) TiC-DC (-)

1.0

charge induced per ion (e)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 36 of 36

anode

0.5

0.0

-0.5

cathode -1.0

-1.5 0

2

4

coordination number of counter-ion

36

ACS Paragon Plus Environment

6

8