Unexplored Pathways to Charge Storage in Supercapacitors - The

Dec 9, 2018 - The observations unearth the limitations associated with the previous theoretical models, which are routinely employed to analyze the ...
0 downloads 0 Views 5MB Size
Subscriber access provided by University of Winnipeg Library

C: Energy Conversion and Storage; Energy and Charge Transport

Unexplored Pathways to Charge Storage in Supercapacitors Anvesh Dixit, Sagnik Middya, Shirsendu Mitra, Surjendu Maity, Mitradip Bhattacharjee, and Dipankar Bandyopadhyay J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b10326 • Publication Date (Web): 09 Dec 2018 Downloaded from http://pubs.acs.org on December 9, 2018

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 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 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.

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 57 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

Unexplored Pathways to Charge Storage in Supercapacitors Anvesh Dixit,a Sagnik Middya,b Shirsendu Mitra,a Surjendu Maity,b Mitradip Bhattacharjee,b Dipankar Bandyopadhyay a,b* aDepartment

bCentre

of Chemical Engineering, Indian Institute of Technology Guwahati, India for Nanotechnology, Indian Institute of Technology Guwahati, India

*Corresponding author. Email: [email protected]

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 2 of 57

Abstract The study reveals various unexplored pathways to energy storage in the parallel and curved plate supercapacitors (SCs). The spatiotemporal variations in the electric field intensity of such SCs were found to have a significant influence on their performance. The observations unearth the limitations associated with the previous theoretical models, which are routinely employed to analyze the performance of SCs by considering electrical double layers (EDLs) as capacitors near the electrodes. The time-dependent electrochemical behaviors of SCs obtained from the Nyquist and Bode diagrams of electrochemical impedance spectroscopy showed, (i) electrode polarization at the higher frequency sweeps, (ii) immobile Helmholtz layer formation at the midfrequency zone, and (iii) formation of diffuse layer of EDL at low-frequency-regime. The results suggest that charge storage of SCs heavily depend upon electrode geometry, type of electrolyte, electrolyte concentration, electrode separation, separator type, and dielectric relaxation of the electrolyte. A theoretical model composed of Poisson-Nernst-Plank equations for the electric field in electrolyte and Laplace equation for the electric field in electrodes were coupled with NavierStokes equations for the fluid flow was numerically solved with appropriate boundary conditions to uncover the pathways to supercapacitance during the experiments. The experimental and theoretical studies together reveal that the use of the potential drop across the EDL originating

ACS Paragon Plus Environment

Page 3 of 57 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

from the opposing electric fields due to electrode polarization and EDL formation could provide more accurate pathways to supercapacitance of such SCs.

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 4 of 57

Introduction The vision of a greener world1–4 contemplate the usage of high quality ubiquitous renewable energies replacing the popular fossil or nuclear fuels.5,6 In this direction, of late, extensive efforts have been put forward to obtain grid-scale7 energy from the wind,8 solar,9 solar-thermal,10 geothermal,11 tidal,12

or ocean-thermal13 resources. Further, the technologies have been

miniaturized to improve efficacy through process intensification and VLSI.14,15 However, most of these renewable resources often lead to a power having quality less than expected owing to their temporal fluctuations due to various environmental reasons.15 In order to resolve this issue, the renewable power generators are often integrated with the energy storage devices before use. Recent studies indicate that the supercapacitors (SCs) can be one of the better alternatives to store energy because they bridge the gap between an energy source and user.16 A number of recent works report that the integration of SCs with energy generators can achieve higher device efficiency,17,18 suitable for futuristic portable19 or large scale20 applications. In particular, one of the major focus of research related to SCs has been directed towards the invention materials for electrode preparation21,22, which possesses superior surface-to-volume ratio and electrical conductivity than the existing ones.16,23–28 In this direction, among the other materials, graphene and its derivatives21 have been extensively used as the active material

ACS Paragon Plus Environment

Page 5 of 57 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

because they have the capacity to achieve a theoretical intrinsic capacitance up to 21 µF/cm2.26 Further, extensive research efforts have been observed in optimizing the electrical conductivity and ion transport properties of the electrolyte between the SC electrodes.29 The design of separators between the electrodes for more efficient ion transport is found to be another important aspect in this regard.30 Importantly, efforts have also been made to uncover the underlying physics associated with the charging and discharging of the SCs, which may enable further improvements in their performance.31 For example, according to the current understanding of charge storage mechanism of the SCs, the energy is stored at the electrode-electrolyte interface due to the accumulation of opposite charges in the electrode surface and electrical double layer (EDL).32 In this situation, the capacitance per unit square area of the charge separation at the electrodeelectrolyte interface is defined as the intrinsic capacitance of an SC. The performance of the SCs are experimentally characterized with the help of the electrochemical impedance spectroscopy (EIS), which helps in evaluating the equivalent series resistances, charge transfer resistance, Warburg resistance16,33 alongside providing the real and imaginary capacitances.34

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 57

Figure 1. Image (a) shows the profile for the net electric field potential (ψnet) between a pair of electrodes representing a Helmholtz electrical double layer (EDL). Image (b) illustrates the equivalent circuit derived from the system shown in the image (a) using two capacitor model (TCM). Image (c) shows the distribution of ions in a diffuse layer following the Gouy-Chapman model. Image (d) represents the Gouy-Chapman-Stern (GCS) model consisting of an immobile Stern layer and mobile diffuse layer in the EDL.

Traditionally, the charge separations at the electrodes are theoretically analyzed by the ‘two capacitor model’ (TCM), as schematically shown in Figures 1(a) and in the equivalent circuit diagram in the Figure 1(b). The subsequent variation in the electric field potential (ψ) across the SC is also depicted in the figure. The TCM model considers two immobile layers of opposite charges at the electrode-electrolyte interfaces, analogous to a solid state capacitor. The TCM also assumes the pristine Helmholtz model of EDL to be valid.32 Further, the TCM model shows that the theoretical capacitance of a SC can only be the functions of electrode surface area and Helmholtz layer thickness. However, experimental EIS studies in the prior-art16,23,33 suggest that

ACS Paragon Plus Environment

Page 7 of 57 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 capacitance of the SCs can also vary with scan rate, frequency of electric field intensity, electric field potential, type of electrolyte, current discharge rates, and pressure between electrodes, among other parameters. Certainly, a better physical model interpreting the data obtained from the EIS during the charging and discharging cycles of SCs is the need of the hour.37 In this direction, the initial attempt was to use the Gouy-Chapman model for the EDL in the TCM having only one mobile layer of ionic distribution following the Maxwell-Boltzmann statistics near the electrodes, as shown in the Figure 1(c). The major drawback of this model was found to be the lack of ionic interactions, which is significant for the SCs using higher electrolyte concentration.38,39 The present understanding of the formation of EDL22-25,40 around the electrodes is most comprehensively explained by the Gouy-Chapman-Stern (GCS) model,41 as shown in the Figure 1(d). This model combines the Helmholtz and Gouy-Chapman models and suggests that the EDL consists of an immobile Stern layer near the electrode and a diffuse layer. Thus, the most recent works explain all the EIS data related to the SCs employing the details of the EDL through the GCS model while considering interactions between the solvated ions inside the electrolyte.41,42 However, most of these studies employ a parameter ‘differential capacitance’43 to account for most of the non-idealities associated with capacitance measurement. Clearly, a quantitative SC performance based on the time-dependent parameters is missing in the prior-art

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 57

owing to the complexities associated with the dynamics of the ion transport near the electrodeelectrolyte interface during charge-discharge cycles. A recent seminal theoretical work on the evaluation of the electrode, electrolyte and diffuse layer resistances from the Nyquist plot of EIS is directed toward this end.44 Further, significant efforts have been made to study the characteristics of a metal-electrolyte interfaces in a SC and batteries using EIS.45,46 These studies underlines the importance of time-dependent analysis to evaluate the performance of electrochemical devices. In view of this background art, the present study uncovers the physical mechanisms of ion transport in parallel plate (PPS) and curved supercapacitors (CPS) with the help of EIS spectroscopy supplemented by a set of analytical and numerical models to describe these experimental behaviours. Initially, the EIS has been employed to uncover the charge storage mechanisms at different time-scales or frequencies47 for different electrode geometries, electrode separation distances, electrolyte concentrations, types of electrolytes, and separators. The Nyquist and Bode diagrams obtained from EIS uncovered that the electrode polarization for SCs took place at smaller time scales, which was followed by the formation of the immobile Stern layer of the EDL at the intermediate time scales. Finally, the diffuse layer was formed at larger time scales. The time dependent ionic movements during the charge separation and subsequent

ACS Paragon Plus Environment

Page 9 of 57 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

dielectric relaxation of the electrolyte is found to play a significant role in the performance of the SC. The proposed analytical and numerical models evaluate the spatiotemporal escalation in the capacitance of the SCs based on the difference between the electric fields due to, (i) polarization at the electrodes under the applied field and (ii) subsequent generation of the electric field in the electrolyte during the formation of the EDL during dielectric relaxation. The theoretical models proposed here is found explain the experimental observations in a more comprehensive manner than previously reported ones.

Methods Preparation and characterization of SC: The chemicals used in the experiments have been outlined in the Section A of the Electronic Supporting Information (ESI). Graphene oxide (GO) was synthesized using modified Hummer’s method.48,49 The detailed methodology of thermally reducing graphene oxide to produce reduced graphene oxide (rGO) and subsequently using it for preparing the active material, has been described in Section B of ESI. The active materials were coated on processed Cu plates, acting as the current collectors, to fabricate functional electrodes for parallel plate (PPS) and curved plate (CPS) supercapacitors. A step-wise demonstration of the same is given in Section B of the ESI. During the study, separators were soaked with an electrolyte, and then the electrode-separator-electrode configuration led to the formation of the

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 10 of 57

proposed PPS. The edges of a PPS were bent against the graph paper to maintain a constant radius of curvature (κ) for the CPS. The active material coated electrode was characterized using Raman spectroscopy (LABram, Horriba, Japan) and the details are presented in Section C of the ESI. Design of Separators: Whatman filter paper (FP), Kimtech tissue paper (TP), and PVDF were used as the separator materials between the electrodes of PPS and CPS. PVDF nanofibers were prepared using electrospinning device (ESPIN NANOTECH - Super ES2) with 16% wt./vol. PVDF polymer in DMF and Acetone solvent (2:3 ratio). The operating voltage and distance between the ejection tip and collector drum were maintained at 20 kV and 15 cm respectively. The FESEM images of the FP, TP and PVDF surfaces are shown in the images (a) – (d) of Figure S3 in the ESI. Table S1 of the ESI summarizes the porosity of these surfaces. The thicknesses of FP, TP, and PVDF surface after wetting with the electrolyte were found to be ~100 µm, ~45 µm, and 100 µm respectively, which were measured under optical microscope (Leica DM 2500). The variations in separation distance between the electrodes of the PPS and CPS were achieved by inserting different numbers (1-5) of wet FP layers. In order to keep the distance between the electrodes (de) similar for the PPS with different separators, we employed two layers of TP (de ~90 µm) for a PPS when compared with the PPS with one layer of FP or PVDF (de ~ 100 µm). Although, in most

ACS Paragon Plus Environment

Page 11 of 57 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

of the experiments we employed single layer of a separator or multilayers of the same separator, we also performed experiments with composite separators. For example, a comparative study was carried out with two TP layers or layers of FP, PFP, and PVDF as the composite separator. Electrochemical analysis: The electrochemical studies were performed using a potentiostat (CHInstruments, 600C). Cyclic Voltammetry analysis was performed for the PPS and CPS with KOH electrolyte in the range of electric field potential (φ) of, 0 V to -0.3 V. The CV analysis was not performed in the high voltage range because there was a possibility of Cu-corrosion in the basic medium, which could lead to and excess current during charging cycle and an asymmetric CV curve. The electrochemical impedance spectroscopy was performed in a frequency range from 100 kHz to 0.1 Hz with an electric field potential amplitude of 0.01 V.

Results and discussion 1. Electrode geometry and distance: In the beginning of this section, we discuss the effects of the electrode geometry and distance on the capacitance of the PPS and CPS. Towards that end, a series of CV analysis was performed for the supercapacitors (SCs) with different radius of curvatures (κ ), ranging from 0 m-1 (PPS) to 0.2 m-1 (CPS) with same mass loading (2 mg/cm2) of active material and electrode separation distance, de = 100 µm. Figures 2(Ia) and 2(Ib) display the typical CV characteristics, that is the

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 12 of 57

variations in specific current, 𝐼𝑆 = 𝐼0/𝑚, with the applied electric potential, 𝜑. The results were obtained at two separate scan rates, 𝜈 = 50 mV/s and 100 mV/s, respectively, for the different PPS and CPS configurations. A smaller area under the discharge current curve in the positive yaxis confirmed a significant decrease in capacitance as κ increased from 0 m-1 (PPS) to 0.2 m-1 (CPS). This observation was true for both a number of scan rates in which two are reported in this manuscript.

Figure 2. Cyclic voltammetry (CV) measurements associated with the different supercapacitor arrangements such as PPS and CPS, as shown in the previous figure. Images (Ia) and (Ib) show the variations in the specific current, IS, with the applied electric field potential, φ, when the electrode separation distance was, de = 100 µm. Here, IS = I0/m, where I0 is the measured current from CV and m is the mass of the active material. The images (Ia) and (Ib) correspond to the scan rates, ν = 50 mV/s and 100 mV/s, respectively. The radius of curvature (κ) were measured, as shown in the inset of the image (Ia) where κ = 0 m-1 represents PPS while κ = 0.035 m-1, 0.123 m-1, and 0.2 m-1 show different CPS configurations. Images (IIa) and (IIb) show the CV when κ =

ACS Paragon Plus Environment

Page 13 of 57 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.035 for de = 100 µm and 500 µm, and scan rates, ν = 50 mV/s and 100 mV/s, respectively. The separator used for all the experiments was filter paper (FP) soaked in 4M KOH.

It may be noted here that previously employed two-capacitor model (TCM) is in general unable to predict this reduction in capacitance as the curvature of the SCs increase from PPS to CPS. This is because, in TCM, the EDLs near the electrodes are modelled like a pair of infinitely long parallel plate capacitors. In this regard, the areas and thicknesses of the EDLs correspond to the surface areas of electrodes and the thicknesses of the Helmholtz layers, respectively.50 This assumption also holds good in the case of CPS, since the Helmholtz layer thickness is of the order of micro to nanometer51 while the electrode dimensions are in centimeter. Subsequently, the TCM model assume the same surface area and Helmholtz layer thickness for both PPS and CPS, which is the major reason behind its limitation to differentiate the capacitance of these configurations. However, the experiments shown in the Figure 2 uncovered that the capacitance in the PPS and CPS were different. A quantitative variation in capacitance and its comparison with capacitance theoretically calculated using TCM is shown in Figure S4 of the ESI. We anticipated that the electric field in the different SC configurations varied with the change in the geometry of the electrodes which in turn influenced the variation in the capacitance. Later, in the section 3.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

Page 14 of 57

we propose a more comprehensive and accurate pathway based on electric field of the SCs to evaluate their capacitance. In order to identify the other limitations of the TCM model, we further explored the effects of de for the CPS. As anticipated, the Figures 2(IIa) and 2(IIb) uncovered that for all scan rates the capacitance of the CPS decreased when de was increased. Again, the reduction in the capacitance of the CPS with increase in de could not be explained by the TCM because it did not account for de in the calculations of capacitance. The experiments corroborated that the TCM model was insufficient to provide a comprehensive understanding of the charge storage of the SCs considered. In the PPS configuration, the electric field and its intensity were not expected to vary significantly across the electrodes owing to their planar geometry and constant separation distance. In contrast, for the CPS, the electric field intensity in between curved electrodes were expected to vary for a given applied electric field potential, φ (inset of Figure S6 of ESI). Thus, the theoretical model proposed later consider the electric field and its intensity to evaluate the capacitance of the SCs. 2. Frequency dependent PPS performance: In order to prove our hypothesis further, we attempted to understand electrochemical performance of the SCs more comprehensively. For this purpose, the key mechanisms at the

ACS Paragon Plus Environment

Page 15 of 57 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

molecular scale during the charging process were explored in detail. In this regard, we anticipated that the polarization of the current collector electrodes and the formation of the EDL near the electrode-electrolyte interface were the two important events during charge storage process in the SC. Since these events happen at different time scales, it was also important to analyze them with respect to time. Thus, we employed EIS, which emulates a pair of sinusoidal chargingdischarging cycles at different frequencies (time-scale) to uncover these spatiotemporal characteristics during charge storage.

Figure 3. Electrochemical impedance spectroscopy (EIS) measurements associated with the PPS having different electrolytes and separator materials. Image (a) shows the high-frequency region (> 10 kHz) of Nyquist plots (real (Z’) versus imaginary impedance (-Z”)) obtained from EIS analysis of PPSs with two different electrolytes 4M KOH (circular symbols) and 1M H2SO4 (square symbols). Image (b) illustrates the same for different separators, namely, filter paper (FP), tissue paper (TP), perforated filter paper (PFP), and poly-vinylidene fluoride (PVDF). In all the experiments, de = 100 µm.

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 16 of 57

A typical Nyquist plot obtained from the EIS of SCs consists of a semicircle at high-frequencies followed by a line with nearly 45° slope at the mid-frequency level and a vertical line at the lower frequencies.26 In general, the EIS data obtained for various SCs are employed to evaluate the C1,

C2, R1, R2, and R3, as shown in the Figure 1b after fitting the same with an equivalent circuit obtained from TCM.27,47 However, as we argued in this study that the TCM has limitations in predicting the capacitance in the situations where the electric field intensity varies between the electrodes. The major reasons behind these variations in the electric field can be the polarization of the electrodes and subsequent formation of the EDL at molecular level during the charge storage cycle. The frequency dependent EIS analyses shown in the following section gives a more accurate interpretation of those parameters. 2.1. High-frequency EIS (>10 kHz): An EIS study of otherwise identical SCs having 1M H2SO4 and 4M KOH as electrolytes is shown in Figure 3(a). In this regard, the x-intercept, the real part of impedance (Z’) at highfrequency, is traditionally termed as the equivalent series resistance (ESR) of the SCs.16 Figure 3(b) shows that the x-intercept was similar for both the cases. However, it may be noted here that the ESR includes all the resistances due to the current collector, active material, binder, electrolyte, and separator.24,26,35,50 Importantly, since we used two different electrolytes of different

ACS Paragon Plus Environment

Page 17 of 57 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

electrical resistivity, as shown in the Table S1 of the ESI, the ESR must be different for the two cases analyzed. A similar observation was made when different separator materials with different porosity were used keeping electrode separation, electrolyte, and other parameters same, as shown in Figure 3(b). Again, the figure suggests that the ESR was independent of the separator, which was quite unlikely owing to the use of different materials as separators. Surely, these experimental observations contradict the traditional method of estimating the ESR, which has been frequently employed to evaluate the electrochemical performance of SCs at high-frequency. In this direction, we made an attempt identify the origin of the SC parameters based on the high-frequency results of EIS after a careful analysis on the dynamics of the charging and discharging processes. It is well known that during the charging process, the observed current is a consequence of the accumulation of free charges on the metallic electrodes. Subsequently, an electric field is almost instantaneously established between the electrodes. However, the formation of EDL due to the drifting of ions of the electrolyte in response to this electric field is a much slower process. For example, the typical characteristic time-scale of charge relaxation for an electrolyte can be evaluated from, 𝜆𝐷𝑑𝑒/𝐷, where 𝜆𝐷 is Debye length, 𝑑𝑒 is electrode separation distance, and D is diffusivity.54 Thus, for the experiments reported in the present work, the theoretical time scales for the aqueous electrolytes (𝜆𝐷 ~ 1 – 100 nm, D ~ 103 µm2/s)54 and 𝑑𝑒 ~

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 18 of 57

100 µm we obtained the time scales in the range of, 10-4 to 10-2 s. The frequency at which the EIS intersect Z’-axis was 105 Hz equivalent to 10-5 s, which was at least one order of magnitude smaller than the theoretically predicted ones. Thus, a smaller experimental time scale in the high-frequency EIS analysis as compared to the theoretically predicted charge relaxation time highlighted the following mechanism. It is well known that the electrolyte resistance arises owing to the restriction in the movement of ions of an electrolyte under an electric field. At smaller time scales (i.e. at higher frequencies) these movements are minimal and contribute marginally to the effective electrical resistance. At the same time, the weakness of the ionic movement also rules out the ion-separator interaction. This again leads to the separator’s meagre contribution towards the net electrical resistance. Subsequently, we obtain identical electrical resistances at high-frequency EIS when the separator or the electrolytes were varied, as shown in the Figure 3. In fact, the variations in the electrolyte strength (e.g. 0.5M to 4M KOH) led to marginal variations in the EIS response, as shown in the Figure 4(Ia). However, the Figure 4(Ib) shows that when de was increased while keeping other parameters same, the x-intercept of Z’ in the EIS kept shifting towards higher electrical resistances. These experiments highlighted that indeed the electric field intensity between the electrodes had major

ACS Paragon Plus Environment

Page 19 of 57 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

influence in determining the real impedance (Z’) or the electrical resistance of the PPS. In this situation, the increase in the electrical resistance of the PPS were attributed to the polarization of the electrodes due to lower electric field intensities. The EIS at the high-frequency zone corroborated that the electrical resistance of the PPS would depend on electrode resistance44 and their polarization, which would vary with the applied electric field intensity. In fact, the midfrequency region of EIS also showed the importance of electrode polarization charge storage mechanism, as discussed in the following section. In a Nyquist plot, the semi-circular region in the mid-frequency zone indicates the concurrent dissipation (Z’) and storage (-Z”) of electrical energy.55 An equivalent circuit of such a system could be envisaged as a resistor and a capacitor connected in parallel.56 While the resistor dissipates the electrical energy, the capacitors stores the same. For such a system, the magnitude of real impedance, ΔZ’, equaling the horizontal spread (or diameter) of the semi-circular region (Figure 4(Ib)), is in general termed as the charge transfer resistance.24,29,57 Previous studies associated ΔZ’ with the charge transfer resistance only when there is a redox reaction at the electrode-electrolyte interface. For example, in the cases of pseudo-capacitors ΔZ’ provides the information on the amount of electrical energy is dissipated during the redox reactions.58

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 20 of 57

Figure 4. Images (Ia) and (IIa) represent the Nyquist and Bode (imaginary impedance (-Z”) versus frequency (f)) plots obtained from the EIS of a PPSs with different KOH concentrations (4M, 2M, 1M, and 0.5M) in which de = 100 µm. Images (Ib) and (IIb) represent the Nyquist and Bode plots for de = 100 µm - 500 µm when the electrolyte was 4M KOH. Images (Ic) and (IIc) represent the Bode plots for the PPS when 4M KOH and 1M H2SO4 were used as electrolytes and de =100 µm. In all the experiments FP was used as the separator material.

2.2 Mid-frequency EIS (10 Hz-10 kHz): Importantly, the charge storage in most of the SC configurations, such as the PPS shown here, is non-Faradaic type.47 Thus, we expect that there is no redox reaction near the electrodes of a PPS, except during the dielectric breakdown or electrode corrosion. However, the experiments here uncover that ΔZ’ of the EIS spectra in the mid-frequency zone can also be an important parameter for the non-Faradaic SCs even in absence of any redox reactions between the electrodes and electrolyte. Figure 4(1b) shows that, although the electrolyte concentration and

ACS Paragon Plus Environment

Page 21 of 57 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

over potential were kept similar, the width of semi-circular regions increased with de for a PPS. The experiments suggested that even in absence of a redox reaction between the electrodes and electrolyte the magnitude of ΔZ’ increased. Certainly, the variation in the electrical resistance with

de in the mid-frequency zone of the EIS had some other physical significance at ionic level. In order to reveal the underlying mechanism of this phenomenon, we put up the following set of hypotheses: (i) during the electrode polarization, marginal ionic movements were observed in the PPS, which was previously established with the help of higher frequency results of EIS; (ii) the counter-ions in the electrolyte ‘drifted’ towards electrodes with the variation in the electric field to cause counter-charge accumulation near the polarized electrodes to initiate the formation of EDL with the progressive reduction in the frequency of the EIS; (iii) in the mid-frequency zone, a stronger ionic drift as compared to ionic diffusion ensured the formation of the Stern or immobile Helmholtz layer of EDL near the electrodes; (iv) subsequently, an electric field (Eh) was generated due to the formation of EDL near the electrodes, which was opposite to the one established due to electrode polarization (Ep); (iv) in the process, the bound charges and the ions present in the electrolyte underwent a frequency dependent dielectric relaxation, which led to a temporal increase in the frequency dependent dielectric constant of the electrolyte.

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 22 of 57

In order to establish these arguments, we studied the resistive (Z’ vs. f) and capacitive (-Z” vs.

f) Bode plots of the EIS for the PPS considered in the Figure 4a. Figures 4(Ic) and 4(IIc) show that in high-frequency region, as anticipated there was no ionic movement and the electrochemical responses were similar for both electrolytes. However, the responses started bifurcating for the different electrolytes in the mid-frequency range of 100-1000 Hz. The plots show that different ionic mobility and charge densities of 4M KOH and 1M H2SO4 led to different resistive and capacitive properties, respectively. In fact, the Bode plots uncovered that the experimentally observed dielectric relaxation starts in the time range of 10-3-10-2s, which was of the same order of the theoretically calculated characteristic time-scale for charge relaxation, which was discussed in Figure 4a.54 Importantly, the bifurcation of the resistive and capacitive Bode plots also indicated a variation in the intrinsic property of the electrolytes after ~10-3-10-2s. For example, as the counter-ions move towards the electrodes, the dissipation and storage of electrical energy initiated owing to the movement of ions and subsequently the formation of the EDL. These two processes can be thought of as a resistor and capacitor in parallel, which resulted in the semi-circular nature of the Nyquist plot shown in the Figure 5(Ib).55,56 As a result of this ionic movement, Figure 4(Ic) shows

ACS Paragon Plus Environment

Page 23 of 57 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

an increase in resistance due to energy dissipation while Figure 4(IIc) illustrates the energy storage due to the EDL formation. The plots in the Figure 4(Ib) also indicate that, at the initial stages of the charge storage of the PPS, the concentration gradient of the counter ions near the unshielded electrode surface was rather insignificant. For example, when de = 100 μm the plot shows, Z’ = 0.5 Ω while -Z” ~ 0. With the progress in time, the applied field drifted the counter-ions in the electrolyte more towards polarized electrodes. At this stage, the drift of the ions was expected to be stronger than the diffusion, which enabled the formation of a compact and immobile counter-ionic layer on the electrodes. The end of the semi-circular region of the Nyquist plot signified the completion of the immobile Helmholtz or Stern layer formation due to ion drift. For example, when de = 100 μm the plot shows that the formation of the immobile layer of EDL was complete at 686.6 Hz (0.0014 s) when Z’ = 0.75 Ω while -Z” ~ 0.3 Ω. Following this, the formation of the diffuse layer of the EDL started, which are analyzed by the Warburg element in EIS plots at the low-frequency range where diffusional ionic movement are significant, as discussed later.16,59 Importantly, before analyzing the low-frequency data of EIS, we investigated the effects of applied electric field intensity (Ep) and electrolyte concentration on the dynamics of EDL formation, in the mid-frequency range of EIS. Previously, Figure 4(Ib) pointed out that the width of semi-

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 24 of 57

circular region increased when the electric field intensity between the electrode reduced. This observation is similar to the case for Li-ion batteries where diffusional resistance of electrolyte increases with electrode separation distance.46 The frequencies of the inflection points (fc) at which semi-circular region completed are indicated by arrows on the figure. The plot suggests that fc reduced from 686.6 Hz for de = 100 µm to about 175.8 Hz for de = 500 µm. As completion of semi-circular region marks the end of the Stern layer formation, decreasing frequency of inflection point suggested that the requirement of more time for the formation of Stern layer at a smaller Ep. This observation was in accordance with the predictions from the proposed theoretical model discussed later,54,59,60 which considered the flux of drifting ions to be the function of the applied electrostatic field. In a way ∆Z’ could be associated with the ionic resistance faced during the formation of Stern layer which in turn would be the functions of electrolyte resistance and Ep. Furthermore, the experimental observations in the Bode plot on the imaginary impedance (-Z”) in Figure 4(IIb) suggests steady state concentration of counter-ions in the Stern layer did not vary with Ep.59 In contrast, Figure 4(IIa) suggests that the steady state concentration of counter-ions in the Stern layer increased with the concentration of KOH. Importantly, the Figures 4(Ia) and (IIa) show that although the electrochemical responses were same for all electrolyte concentration at the higher frequency EIS zone, the response varied with the molarity of the electrolyte KOH in the

ACS Paragon Plus Environment

Page 25 of 57 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

mid-frequency region. Thus, the capacitance increased with the concentration of KOH owing to the accumulation of the larger concentration of the counter-ions in the Stern layer. Concisely, the response of PPS in the mid-frequency EIS could uncover the dynamics of the formation of the Stern Layer in an EDL with the variations in the Ep and electrolyte strength. 2.3. Low-frequency EIS (