Beichen
Liu
,
Wenxiao
Guo
,
Seth R.
Anderson
,
Samuel G.
Johnstone
,
Siqi
Wu
,
Megan C.
Herrington
and
Matthew A.
Gebbie
*
Department of Chemical and Biological Engineering, University of Wisconsin–Madison, Madison, Wisconsin 53706, USA. E-mail: gebbie@wisc.edu
First published on 1st December 2023
Electric double layers are crucial to energy storage and electrocatalytic device performance. While double layer formation originates in electrostatic interactions, electric double layer properties are governed by a balance of both electrostatic and entropic driving forces. Favorable ion–surface electrostatic interactions attract counterions to charged surfaces to compensate, or “screen,” potentials, but the confinement of these same ions from a bulk reservoir to the interface incurs an entropic penalty. Here, we use a dicationic imidazolium ionic liquid and its monovalent analogue to explore how cation valence and entropy influence double layer formation and electrochemical reactivity using CO2 electroreduction as a model reaction. We find that divalent and monovalent cations display similar CO2 reduction kinetics but differ vastly in steady-state reactivity due to rapid electrochemically induced precipitation of insulating dicationic (bi)carbonate films. Using in situ surface-enhanced Raman scattering spectroscopy, we find that potential-dependent cation reorientation occurs at similar potentials between the two ionic liquids, but the introduction of a covalent link in the divalent cation imparts a more ordered double layer structure that favors (bi)carbonate precipitation. In mixed monovalent-divalent electrolytes, we find that the divalent cations dominate interfacial properties by preferentially accumulating at surfaces even at very low relative concentrations. Our findings confirm that ion entropy plays a key role in modulating local electrochemical environments. Furthermore, we highlight how double layer properties are sensitive to the properties of counterions that pay the lowest entropic penalty to accumulate at interfaces. Overall, we illustrate that ion entropy provides a new knob to tune reaction microenvironments and unveil how entropy plays a major role in modulating electrochemical reactivity in mixed ion electrolytes.
To date, the electrocatalysis community has demonstrated that CO2 can be reduced to a variety of potentially useful products by judiciously tuning catalyst materials12–15 and electrolyte properties.12,16–18 The formation of CO is typically the first reduction step for the production of more complex C2+ products.19,20 However, reducing CO2 to CO requires significant energy input due to the relative stability of the CO2 double bonds,21,22 making electrocatalytic reduction of CO2 to CO an especially relevant model reaction for fundamental studies.
In most cases, CO2 reduction to CO proceeds through a two-electron, two-proton reduction process (eqn (1)), where the intermediate has been reported to be a negatively-charged radical, CO−2.21,22 The formation of this intermediate, as well as the transformation of CO2 from a nonpolar molecule to the polar CO molecule, means that this reaction is particularly sensitive to potential gradients and therefore presents a unique avenue for the study of electric double layers.
CO2 + 2H+ + 2e− → CO + H2O | (1) |
Notably, CO2 electroreduction is carried out at high cathodic potentials, often exceeding −1 V below open circuit potential.23 These intense electrode polarizations induce the formation of nanostructured electric double layers,24–29 where ion–ion correlations and excluded volume effects play defining roles in determining local potential gradients and reaction microenvironments.30,31 Hence, electric double layer formation at electrocatalytic interfaces often hinges on collective ion self-assembly, which is not considered in classical electric double layer theories such as the Gouy–Chapman–Stern model, where ion size, ion shape, and ion–ion interactions are neglected.32–35
In this work, we aim to expand our understanding of how cation entropy influences ion assembly, electric double layer formation, and ultimately local reaction microenvironments under conditions of high applied bias that are inherent to CO2 electrocatalytic reduction. We selected ionic liquid-mediated CO2 reduction to CO on Ag electrodes as a model system, as ionic liquid molecular structures can be systematically modulated via synthetic protocols.36
Ionic liquids are salts that are often composed of weakly coordinating organic cations and anions and have been heavily studied as excellent promoters of CO2 reduction in both aqueous37–40 and non-aqueous electrolytes.30,41,42 Imidazolium-based salts are the prototypical class of ionic liquids used for mediating CO2 electroreduction due to their propensity to donate protons and coordinate CO2.37 Increasingly, however, studies are highlighting how other organic cations that cannot donate protons or form coordination interactions can also effectively promote CO2 electroreduction,43,44 leaving many open questions surrounding how and why organic cations can facilitate electrocatalytic reactions.
We hypothesize that cation entropy provides an important and underexplored element in explaining why ionic liquids can facilitate efficient CO2 electroreduction. Electric double layer formation hinges on both entropic and electrostatic interactions,45 where the large size and irregular shape of ionic liquid ions can lead to additional dielectric fluctuations and responses which can contribute to electrostatic screening and enhance interfacial potential gradients.32,46,47 An ion confined to the double layer not only experiences electrostatic interactions, but also incurs entropic penalties.48 However, the degree to which the entropic penalty influences double layer formation, and subsequently electrostatic screening, remains poorly understood.
Here, we systematically isolate the influence of entropy on electric double layer formation by comparing electric double layer properties and electrocatalytic CO2 reduction rates in two ionic liquids that are structurally identical except for a covalent linkage that forms a dimer analogue of a common monovalent ionic liquid. Specifically, we examine the role of entropy in double layer formation by comparing the interfacial behavior and reactivity of 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide ([BMIm][TFSI]) and its dicationic analogue, 1,8-bis(3-methylimidazolium-1-yl)octane bis(trifluoromethanesulfonyl)imide ([BisMImO][TFSI]2) (Fig. 1).
Ions confined to the interface must pay an entropic penalty due to the loss of translational motion of , where kB and T represent the Boltzmann constant (J/K) and temperature (K), respectively. As this entropic penalty pertains to any single molecule regardless of size or shape, [BisMImO]2+ pays an entropic penalty of while [BMIm]+ must pay an entropic penalty of to confine two positive charges to the interface. Therefore, it is more thermodynamically favorable to confine the [BisMImO]2+ cation to the double layer than 2 [BMIm]+ ions for the same amount of charge screening.
We hypothesized that salts with lower entropic penalties per charge carried, such as [BisMImO][TFSI]2, will preferentially accumulate at interfaces, generally leading to thinner and more efficiently-screening double layers. However, the effects of a thinner double layer may be tempered by the additional electrostatic interactions between the imidazolium rings of the [BisMImO]2+ cation, as well as the imposed constraint on the number of positions the imidazoliums can adopt relative to each other in [BisMImO]2+ when compared to [BMIm]+. We evaluated the electrochemical reactivity and interfacial structure of [BisMImO][TFSI]2 and [BMIm][TFSI] in acetonitrile to avoid the competing hydrogen evolution reaction49 and on Ag working electrodes for its selectivity for CO.50
Our findings demonstrate that entropic contributions to double layer formation can dominate interfacial ion assembly and govern reaction rates. Using cyclic voltammetry (CV), we show that [BisMImO][TFSI]2 generally exhibits lower current densities than [BMIm][TFSI], but that current densities begin to converge in more concentrated electrolytes. We then use chronoamperometry (CA) and product analyses such as gas chromatography (GC) to show that CO2 electroreduction in [BisMImO][TFSI]2-containing electrolytes lead to noticeable precipitation of the bicarbonate salt at the electrode surface, which compromises steady state reactivity, while [BMIm][TFSI] does not.
To better understand the structure within the double layer and how it is related to the reaction rate, we performed in situ electrochemical surface-enhanced Raman scattering (SERS) spectroscopy and found that [BisMImO][TFSI]2 more readily forms ordered double layers, especially under inert gas-sparging conditions. Furthermore, in a mixture of the two ionic liquids, we provide evidence that electrocatalytic performance and interfacial structuring are dominated by the ionic liquid with the lower entropic penalty of de-mixing from the bulk electrolyte.
The dicationic ionic liquid, 1,8-bis(3-methylimidazolium-1-yl)octane bis(trifluoromethanesulfonyl)imide ([BisMImO][TFSI]2) was synthesized according to established literature procedures.51,52 Briefly, 1-methylimidazole (99%, Thermo Scientific) and 1,8-dibromooctane (98%, Sigma-Aldrich) were added to a round-bottom flask in a 2:1 molar ratio and allowed to react for 24 hours at room temperature to make 1,8-bis(3-methylimidazolium-1-yl)octane dibromide ([BisMImO][Br]2). The [BisMImO][Br]2 was then washed with acetone and dried, resulting in a dry and white powder.
To conduct the anion exchange, [BisMImO][Br]2 and [Li][TFSI] (>98%, TCI) in a 1:2 molar ratio were added to a round bottom flask and fully dissolved in MilliQ H2O. The mixture was stirred at room temperature for 3 hours. The water layer was decanted off and the resulting ionic liquid was washed with MilliQ H2O. After several washes, the water layer was tested with a silver nitrate halide test to confirm the removal of halides.
The crude [BisMImO][TFSI]2 was then purified by first dissolving it in acetonitrile and then stirring it with activated charcoal (Sigma-Aldrich) for 2–3 days. After filtering out the activated charcoal, the solution was dried using a rotary evaporator (BUCHI R-300M Rotavapor) followed by a vacuum drying oven (Yamato ADP200C) at 80 °C for 2 days. The resulting [BisMImO][TFSI]2 was a viscous and colorless oil that slowly began to solidify upon storage. A more detailed method is presented in ESI.†
The water content of HPLC grade acetonitrile (>99.9%, Sigma-Aldrich) was adjusted to be 900–1000 ppm according to our previous protocol to ensure that protons were available for the reaction.30 Potassium bromide was purchased from MP Biomedicals and silver nitrate was from Alfa Aesar.
For cyclic voltammetry and steady-state electrolysis, we used a 3 mm polycrystalline Ag electrode (geometric surface area of 0.071 cm2) purchased from BASi as the working electrode, which was polished to a mirror finish using 15 μm, 3 μm, and 1 μm diamond polish followed by 0.05 μm alumina (BASi) before ultrasonication in water. All electrolytes were purged with either Ar or CO2 prior to electrochemical measurements. CV measurements were conducted in one-compartment borosilicate electrochemical cells with a 3D-printed cap to ensure consistent electrode placement. All presented CV measurements were carried out at a scan rate of 100 mV s−1 and the results are presented in their raw form, as IR-corrections do not change the conclusions of this study, since the current passed is low enough to negate the minor changes in solution resistance that occur as the ion concentration is altered in the reported electrolytes (Tables S1, S2 and Fig. S1, S2, ESI†).
The cells were continuously stirred with Teflon stir bars (VWR) and CO2 was continuously purged through the cathode compartment at 10 standard cubic centimeters per minute (sccm). The resulting gas stream outlet fed directly into a gas chromatograph (Multiple Gas Analyzer #5, SRI Instruments), which was equipped with a flame-ionization detector (FID) coupled with a methanizer and a thermal conductivity detector (TCD). The FID was connected with a HayeSep D and a Molesieve 5A column (Restek) while the TCD was connected to a HayeSep D column. The carrier gas was ultrahigh purity grade Ar (Airgas).
Product analysis was conducted every 6 minutes and faradaic efficiencies were calculated using the following equation:56
(2) |
We also used NMR spectroscopy to identify products in the catholyte after electrolysis. The NMR (400 MHz, Bruker Avance) was equipped with a BBFO probe.
Briefly, the electrode was submerged in an aqueous solution of 0.1 M [K][Br] under continuous Ar purging. The electrode was then cycled between −0.3 V and +0.3 V vs. open-circuit potential (OCP) for 4 cycles. The electrode was then held at −0.3 V vs. OCP for 10 seconds to desorb bromide from the silver surface. The electrode was thoroughly rinsed with water after roughening and dried prior to use. A mechanically polished and coiled Pt wire was the counter electrode and a non-aqueous Ag/Ag+ (0.01 M [Ag][NO3]) reference electrode was used.
The electrolytes were purged with either Ar or CO2 for 20 minutes at 10 sccm prior to measurement. SERS spectra were collected using a Thermo-Fisher Scientific DXRxi Raman Imaging Microscope with 532 nm excitation laser at 10× magnification. The laser power was 8.0 mW, the acquisition time was 0.125 seconds, and 25 scans were collected for each spectrum. The potential was swept from −0.5 V to −2.5 V vs. Ag/Ag+ and then back to −0.7 V vs. Ag/Ag+ at intervals of 0.2 V. Only the forward scans are reported in the main text as the trends were reversible.
We found that the [BisMImO][TFSI]2 electrolytes display a maximum current density at similar salt concentrations as [BMIm][TFSI] but that the current densities of [BisMImO][TFSI]2-containing electrolytes are lower than those containing [BMIm][TFSI] at all concentrations (Fig. 2(A)). More interestingly, after accounting for the total number of charge-carrying imidazoliums in each cation, we observe that the discrepancies in current densities between the two ionic liquids are larger at lower concentrations (<1 M) but converge to similar current densities as concentration increases (Fig. 2(B)).
We attribute the discrepancy in CO2 reduction rates at lower imidazolium concentrations to the reorientation of the cations in the double layer and their corresponding role in the reaction. The imidazolium cation is thought to stabilize the CO−2 intermediate through either complexing with the C2 position,37,58 or through hydrogen bonding between CO−2 and the C4 and C5 protons.42,59 The covalent link between two imidazolium rings in the [BisMImO]2+ cation serves to further constrain the imidazolium ring orientation relative to the CO−2 intermediate. As such, the imidazoliums in [BisMImO]2+ have fewer degrees of freedom and hence fewer stabilizing conformations compared to the [BMIm]+ cation. Additionally, at low concentrations, the [BMIm]+ cations remain fully solvated and “free” while the imidazoliums in the [BisMImO]2+ cations experience enhanced interactions due to the covalent linkage, which could promote the formation of electrode-blocking precipitates containing (bi)carbonate anions. The imposed constraint and corresponding entropic limitations are likely the reason behind the attenuated reactivity of the [BisMImO][TFSI]2 at concentrations lower than 1 M of imidazolium cation.
As concentration increases past 1 M, however, the emergence of bulk ion–ion correlations and the formation of multi-ion aggregates makes screening less efficient than the predictions of classical double layer theory.60 At these concentrations, ion–ion interactions between [BMIm]+ cations also decrease screening efficiency such that correlations within [BMIm]+ electrolytes become similar to those in [BisMImO]2+ electrolytes. As a result, increasing concentration towards the neat limit yields converging current densities for the two ionic liquids as double layer screening efficiency decreases and the double layer thickness increases.
We then conducted chronoamperometry (CA) studies combined with product analysis techniques to evaluate the two ionic liquids for steady-state electrolysis performance. CA experiments were conducted using polycrystalline Ag working electrodes in a well-mixed H-cell with continuous CO2 sparging at 10 sccm. The outlet gas was directed into the sample loop of a gas chromatograph to identify gaseous products.
We notice that within 10–15 minutes of applying a potential of −2.5 V vs. Ag/Ag+, the current density and faradaic efficiency of CO2 conversion to CO for all [BisMImO][TFSI]2 electrolytes drop and remain at a low level for the remainder of the experiment (Fig. 3(A)). In contrast, the [BMIm][TFSI] electrolytes maintain a relatively constant current density and faradaic efficiency (Fig. 3(B)). Upon inspection of the working electrode after electrolysis with the dicationic analogue, we observed that a solid precipitate had formed on the surface and completely occluded the working electrode. This precipitate was not observed after electrolysis in [BMIm][TFSI], suggesting that covalently joining the two imidazoliums together has the effect of inducing precipitate formation.
Next, we performed CA at incremented potentials following the reaction onset (Fig. S7, ESI†) and observed this precipitate formation at potentials near the onset potential of CO2 reduction but did not observe the precipitate under Ar-sparged conditions (Fig. S8, ESI†). Further, even low concentrations (0.0625 M) of [BisMImO][TFSI]2 displayed this behavior, as a decrease in current density was observed even between different cycles of CV (Fig. S9, ESI†).
The formation of the obstructing precipitate on the electrode under CO2-sparged and its absence under Ar-sparged conditions indicates that the formation of byproducts during CO2 reduction is the likely cause of the precipitate. Indeed, (bi)carbonate formation and subsequent precipitation of (bi)carbonate salts (Fig. S10, ESI†), which is caused by the reaction between interface-generated [OH]− and dissolved CO2, is a common issue in CO2 electrocatalysis.10,61 The precipitate that is formed lowers the current density and the subsequent concentration of gaseous products, which can approach or fall below the detection limits of our instrumentation and makes rigorous analysis of product distributions challenging. Analysis of the carbonate precipitate suggests that the precipitate can trap products like formate at the electrode–electrolyte interface, which can further decrease the faradaic efficiency for CO (Fig. S11, ESI†). Carbonate salt precipitate formation is also observed in inorganic electrolytes and remains a common source of diminishing reactor performance in steady state CO2 electrolysis.62,63 Our observations are also consistent with our prior work,30,31 where we previously reported the formation of a solid precipitate in electrolytes with 1-ethyl-2,3-dimethylimidazolium tetrafluoroborate, in which the bicarbonate anion paired with the imidazolium cation was insoluble in acetonitrile.
Hence, simply joining two imidazolium rings together dramatically alters the steady-state electrocatalytic behavior. It is difficult to overstate the importance of ion molecular structure and how it impacts both double layer formation and reactivity. Prior works have demonstrated significant cation effects in both aqueous64 and non-aqueous systems,65 which often focus on how changing the cation structure or size can enhance local concentration gradients66 or the availability of protons.67 Here, we show that cation entropy can have similarly drastic effects on CO2 electroreduction by preserving the inherent structure of the imidazolium cation but introducing an entropic constraint. To further investigate the double layer that is formed when the cation is entropically-constrained, we used surface-enhanced Raman scattering (SERS) spectroscopy to examine the structure of the double layer at various applied potentials.
In the alkyl substituent region, we observed the rise of a peak at 1003 cm−1, which is particularly pronounced in [BisMImO][TFSI]2 electrolytes and is observed primarily in Ar-sparged conditions (Fig. 4). Upon the reverse potential scan in Ar-sparged electrolytes, the intensity of the 1003 cm−1 peak diminishes. Interestingly, we found that there was a marked difference in the growth of the 1003 cm−1 peak when the electrolyte had been sparged with Ar compared to CO2 (Fig. 4). Additionally, while we do observe peaks that may be assigned to Raman-active modes of carbonate, these peaks overlap with those associated with the alkyl chain region of the cation (Fig. S12, ESI†) and definitive assignment of these peaks in CO2-sparged solutions is therefore challenging.
To analyze the growth of the 1003 cm−1 peak, we took the ratio of the peak intensity of the 1003 cm−1 peak relative to the peak intensity of the 1388 cm−1 peak, which is associated with the stretching mode of the imidazolium ring structure and remains largely unchanged with applied potential. We find that the peak intensity of the 1003 cm−1 peak surpasses that of the 1388 cm−1 peak at −1.9 V vs. Ag/Ag+ (Fig. S13, ESI†). After reaching this potential, the intensity of the 1003 cm−1 peak becomes several times larger than the 1388 cm−1 peak in Ar-sparged solutions. A similar trend is observed in [BMIm][TFSI] solutions, but the growth of the 1003 cm−1 peak in Ar-sparged solutions is much more subdued (Fig. S14, ESI†).
The 1003 cm−1 peak is associated with the alkyl substituents on the imidazolium and has been previously shown to correspond to C–C stretching.68,69 The pronounced 1003 cm−1 peaks likely arise from interfacial segregation of the alkyl chains from the imidazolium ring structures given that the alkyl chains, especially for [BisMImO]2+, are less polar than the positively charged imidazolium ring. Additionally, the formation of segregated alkyl chain domains have previously been observed in monovalent imidazolium ionic liquid systems.28,70 There was also a greater enhancement of the 1003 cm−1 signal in [BisMImO]2+, suggesting that the connecting alkyl chain imparts more organization to the double layer. Notably, the peak at 1003 cm−1 disappeared altogether after −1.9 V vs. Ag/Ag+ under CO2-sparged conditions, which is approximately the onset potential for CO2 electroreduction.
The CO2 reduction process requires the transport of molecules across the double layer and forms primarily CO but can also generate unwanted byproducts such as [OH]− (eqn (3)) and, by recombination of the [OH]− with CO2, (bi)carbonate.31,61
CO2 + H2O + 2e− → CO + 2OH− | (3) |
Transporting CO2 towards the electrode and transporting products and byproducts away from the electrode can disrupt the structure that would have formed in the absence of a reaction in the double layer. Furthermore, the formation of negatively charged products like the CO−2 intermediate or (bi)carbonates can electrostatically balance the charge of the cations already in the double layer, therefore negating the screening effect of the cation. As a result, more cations are recruited to the double layer as the neutral clusters of ions diffuse away from the surface, thus creating a more dynamic and disrupted double layer structure.
We find that the increase in the intensity of the 1003 cm−1 peak under Ar-sparged conditions was also observed in N2-sparged electrolytes (Fig. S15, ESI†). Combined with the CO2-sparged results, this suggests that the double layer under inert conditions is fundamentally different from double layer under reactive conditions, in which the transport of reactants and (by)products creates a more dynamic interface. This finding has far-reaching implications, the most significant of which is that double layer models developed for inert conditions, or capacitive double layers, are not directly applicable to dynamic reactive double layers.
The imidazolium ring region of 1344 to 1580 cm−1 can be used to qualitatively understand the reorientation of the imidazolium ring due to the development of an electric field as a potential is applied to the electrode. In SERS spectroscopy, the vibration modes with tensors perpendicular to the electrode plane are enhanced so trends in signal intensities can be used to gauge the strength of the Raman mode. We observed the red-shifting of the 1580 cm−1 peak to 1561 cm−1, which is indicative of increased proximity between the C4 and C5 protons with the surface72–74 (Fig. 5). We also observed a decrease in the intensity of the 1415 cm−1 peak, which is associated with the in-plane stretching mode of the imidazolium,74,75 and a corresponding increase in the intensity of the 1344 cm−1 peak, associated with CH2(N) stretching,74–76 as the potential becomes more cathodic for both [BisMImO][TFSI]2 and [BMIm][TFSI] electrolytes. Taken together, these trends in ring Raman modes as a function of applied potential suggests that the imidazolium ring begins to adopt an orientation parallel to the electrode surface as applied potential increases,73,75,76 likely with the alkyl substituents roughly also parallel to the electrode.
We used the ratio of the intensity of the 1344 and 1415 cm−1 peaks as a metric for assessing the strength of the electric field, as stronger electric fields can reorient an imidazolium ring at less cathodic potentials than weaker electric fields (Fig. 5(B)). Although we observe the same general trend of diminishing 1415 cm−1 peak intensities between the two ionic liquids, the potentials at which the ratio of the peak intensities become greater than 1 are quite different at various concentrations.
For example, the intensity inversion for 0.125 M [BisMImO][TFSI]2 and 0.250 M [BMIm][TFSI] are very similar and occur at approximately −1.0 V vs. Ag/Ag+. At higher concentrations, however, we observed a delayed inversion potential of −1.6 V vs. Ag/Ag+ for the 1.65 M [BMIm][TFSI] electrolyte. The inversion potential for 0.825 M [BisMImO][TFSI]2 occurs at approximately −1.0 V vs. Ag/Ag+, indicating little difference between the low and intermediate concentrations of [BisMImO][TFSI]2.
The imidazolium ring reorientation inversion potential allows us to gauge how bulk ion–ion correlations impact double layer formation. For the lower concentrations of [BisMImO][TFSI]2 and [BMIm][TFSI], the intensity inversion occurs at similar potentials, as ion–ion correlations have not yet begun to dominate. In the 1.65 M [BMIm][TFSI] electrolyte, we observed a delayed inversion potential, which can be explained by ion–ion interactions, as well as the formation of clusters and aggregates in the bulk electrolyte that decrease screening efficiency and thus lower the strength of the interfacial electric field. Although there are the same number of charges, 0.825 M [BisMImO][TFSI]2 did not have such a delay and the potential of peak intensity inversion was more similar to lower concentration electrolytes, indicating a stronger electric field in 0.825 M [BisMImO][TFSI]2 than 1.65 M [BMIm][TFSI].
One explanation for this is that the linkage between the two imidazoliums in [BisMImO]2+ can help with the potential-dependent reorientation through a decrease in accessible orientations of imidazolium rings relative to [BMIm]+, hence imparting a “pre-structure” to the double layer. The ratio of the 1344 cm−1 and 1415 cm−1 peaks, however, was smaller than that of the lower concentrations especially at potentials more cathodic than −1.5 V vs. Ag/Ag+. Thus, the SERS results indicate that the extent of this reorientation is smaller in more concentrated electrolytes and is likely hindered by ion–ion interactions.
First, we find that by adding 5 mM [BisMImO][TFSI]2 to 120 mM [BMIm][TFSI], the steady-state current density drops to 25% of the current density observed for the 125 mM [BMIm][TFSI] electrolyte. As the ratio of [BisMImO][TFSI]2 increases, the steady state current density decreases until it becomes smaller than −2.5 mA cm−2 in a mixture of 25 mM [BisMImO][TFSI]2 and 100 mM [BMIm][TFSI]. Notably, faradaic efficiency for CO decreases gradually until it reaches approximately 50% at 125 mM [BisMImO][TFSI]2 (Fig. 6).
We observed a pronounced decrease in current density with the addition of relatively small amounts of [BisMImO][TFSI]2 as well as precipitate formation, which together indicates that a distinct [BisMImO]-enriched phase forms on the electrode surface. Interestingly, the potential applied to the electrode surface dictated the physical appearance of the precipitate in solutions containing only [BisMImO][TFSI]2. For instance, a solid white crystal nucleated and propagated on the surface when a potential of −2.5 V vs. Ag/Ag+ was held under CO2-sparged conditions (Fig. S3, ESI†). If the potential was quickly stepped to −2.7 V vs. Ag/Ag+ and held, however, a dynamic precipitate formation process was observed that occurred in two steps. Evidence of liquid–liquid phase separation was observed as islands of a liquid phase formed on the electrode surface and began to coalesce. Within these liquid droplets, solid precipitates began to form, ultimately resulting in the formation of an orange gel-like film spanning the entire electrode surface, likely composed of the (bi)carbonate salt and decomposed imidazolium (Fig. S16, Supplemental Video, ESI†). These same trends were observed in mixtures of [BisMImO][TFSI]2 and [BMIm][TFSI].
Hence, we find that decreasing the entropic penalty of constraining cations within electric double layers appears to amplify driving forces for electrochemically induced phase separation. Electrochemically induced phase separation has been previously demonstrated by Leibler and coworkers to be possible in binary mixtures of silicone oil and paraffin oil,77 which was followed up with a theoretical model for electric field gradient-induced phase separation in ion-containing mixtures.78 In experiments by Lahiri and coworkers, the mechanism for potential-dependent phase separation in ionic liquid-containing electrolytes was through the formation of unstable bonds, resulting in the eventual decomposition of the electrolyte into two phases.79
Even more interestingly, electrolytes containing only 10 mM [BisMImO][TFSI]2 (with 115 mM [BMIm][TFSI]) formed a clear liquid phase on the electrode surface that was distinguishable from the bulk electrolyte under CO2-sparged conditions but not in Ar-sparged conditions. We hypothesize that this is due to the formation of the [BisMImO]2+ (bi)carbonate salt, where the concentration of the salt is not sufficient to induce precipitation. Additionally, the presence of [BMIm]+ could disrupt the formation of a solid crystal since the (bi)carbonate salt of [BMIm]+ does not precipitate out of solution, leading to the metastable liquid–liquid phase separation process.
To the best of our knowledge, liquid–liquid phase separation due to the active electrocatalytic processes in ionic liquid electrolytes has not been reported in literature. This separate liquid phase was fluid and moved when the cell was lightly disturbed before reverting to its original state (Supplemental Video, ESI†). A white crystal began to nucleate and grow after more than 30 minutes of applied potential of −2.7 V vs. Ag/Ag+. Further increases in the ratio of [BisMImO][TFSI]2 to [BMIm][TFSI] in the electrolyte mixtures began to yield precipitates that were indistinguishable from that of electrolytes containing only [BisMImO][TFSI]2. Hence, we conclude that the formation of solid precipitates in these mixed electrolytes appears to pass through a non-classical multistep nucleation and growth process, as explained by De Yoreo and colleagues for carbonate precipitation in biological systems.80
The low concentration required to observe this phase separation indicates that the [BisMImO][TFSI]2 preferentially accumulates in the double layer but crystal formation is suppressed or delayed due to the presence of [BMIm]+ cations with lower propensity to form solid (bi)carbonate crystals. Although the full mechanism is the subject of ongoing work, the potential-dependent liquid–liquid phase separation we observe in these mixed ionic liquids appears to have intriguing connections to the process of coacervate formation in biology, where changes in pH and ionic strength drive formation of distinct liquid phases.80 We therefore expect that the interface would contain higher concentrations of [BisMImO][TFSI]2 than the bulk.
To investigate the composition of the double layer, we collected SERS spectra of the electrolyte mixtures, focusing on the electrolyte compositions with low ratios of [BisMImO][TFSI]2:[BMIm][TFSI]. Given that even mixtures of low concentrations of [BisMImO][TFSI]2 exhibited electrochemical behavior more like that of [BisMImO][TFSI]2-only electrolytes, we expected that the SERS spectra of the interface would reflect an enhanced abundance of interfacial [BisMImO][TFSI]2 in these mixed electrolytes.
Under Ar-sparged conditions, we indeed found that the electrolytes containing small amounts of the [BisMImO][TFSI]2 behaved more like the electrolyte containing primarily [BisMImO][TFSI]2. In short, the 1003 cm−1 peak outgrew the 1388 cm−1 peak in intensity in all solutions containing [BisMImO][TFSI]2 regardless of whether [BMIm][TFSI] was also present (Fig. S17, ESI†). We observed the appearance of the 1003 cm−1 peak under Ar-sparged conditions in a mixture of 10 mM [BisMImO][TFSI]2 and 115 mM [BMIm][TFSI] (Fig. 6(B)). This suggests there is organization of the alkyl chain domains of the cations and an enhanced presence of the [BisMImO]+ cation at the interface. In other mixtures, the ratio of the peak intensities in the mixed electrolytes began to approach that of the 125 mM [BisMImO][TFSI]2 solution as the concentration of [BisMImO][TFSI]2 increased (Fig. S18, ESI†), suggesting that the [BisMImO]2+ concentrates within the double layer. Otherwise, the low ratio of [BisMImO][TFSI]2 to [BMIm][TFSI] should yield SERS spectra more like an electrolyte containing only [BMIm][TFSI].
These findings, combined with the electrolysis data, point to the preferential accumulation of the lower entropic penalty ionic liquid, [BisMImO][TFSI]2, at the electrode/electrolyte interface. In short, the electrochemical behavior of a mixed electrolyte is dictated by the salt with the lowest entropic penalty for double layer formation (Fig. 7). Intriguingly, our results on entropic considerations of double layer formation may also help explain observations made in aqueous studies of mixed alkali metal salts, where the reactivity of mixed electrolytes containing mixtures of Li+ and larger alkali metal cations is dominated by the reactivity of the larger alkali metal cation,81 or the enhancement of CO2 electroreduction in the presence of multivalent ions.61
When confined within the double layer, the [BMIm]+ ion has more orientational degrees of freedom. The imidazoliums in the [BisMImO]2+ ion, on the other hand, have a reduced number of available orientations available relative to one another. The covalent bond in [BisMImO]2+ also serves to impose an additional layer of ion correlations, resulting in smaller current densities at lower concentrations due to the formation of precipitates. These effects become less pronounced at high concentrations, however, as the emergence of ion–ion correlations in [BMIm]+ attenuates the difference between CO2 reduction in [BisMImO]2+ and [BMIm]+, leading to converging reaction rates in the highly concentrated ionic liquid limit.
Analysis of SERS results in [BisMImO][TFSI]2 electrolytes also point towards a fundamental difference in the structure of the double layer under inert or reactive conditions. We found that under inert conditions, the long alkyl chains of the [BisMImO]2+ formed segregated domains or local nanostructures. When CO2 is present, however, the double layer structure becomes disrupted. This can be attributed to the dynamic nature of a reactive double layer, in which the transport of reactants and products disturbs the interfacial structure, and formation of charged byproducts disrupts the local electric field. This difference in double layer structure under two different gas sparging conditions highlights the complex nature of double layers and suggests that conventional models of double layers, which assume no faradaic reactions, may not be directly applicable to reactive double layers.
In [BisMImO][TFSI]2 electrolytes, we observed solid precipitate formation, which passivated the electrode surface and prevented further CO2 electroreduction from occurring. These phase separations only occurred under CO2-sparged conditions, suggesting that byproducts of the CO2 reduction reaction facilitate the phase separation. Additionally, the overall appearance of the solid precipitate differs at different applied potentials, indicating that the presence of other electrochemical products, such as those made by the deactivation of imidazolium, can also modify for surface passivation processes.
The inclusion of an ion like [BMIm]+ in mixed electrolytes disrupted the solid [BisMImO]–(bi)carbonate precipitate crystal structure, thus leading to the liquid–liquid phase separation observed in the 10 mM [BisMImO][TFSI]2 115 mM [BMIm][TFSI] solutions. Precipitate formation poses a serious hurdle in the adoption of CO2 electrolysis cells, as formation of solid (bi)carbonate salts creates a passivating layer on the electrode surface, ultimately leading to electrolytic cell failure.10,61 Our mixed electrolyte experiments show that this precipitation process occurs via a multistep mechanism, in which liquid–liquid phase separation first occurs and is followed by precipitate nucleation. Multistep nucleation processes have been demonstrated to have much lower energy barriers than classical nucleation,80 which may explain the prevalence of the precipitate problem and also offer methods for prevention. We have shown that mixing two salts can have a stabilizing effect on the metastable liquid formed at the electrode surface, but higher compositional complexity may be a more suitable approach for completely suppressing nucleation.
Our findings highlight the key role entropy plays in determining interfacial potential gradients and electrochemical reactivity. We also demonstrate that entropy–reactivity relationships offer a new approach to tuning electrocatalytic reactions, in which the design of an electrochemical system must account for both ion entropy of confinement to the double layer, as well as the ability of ions to reorient at the interface. Double layer models must take both entropic considerations and reactive conditions into account to understand the dynamic nature of the reactions occurring in double layers. Ions at the interface can facilitate or hinder a reaction, as evidenced by the onset of liquid–liquid phase separation followed by solid precipitate formation in [BisMImO][TFSI]2 electrolytes.
Our observation of multi-step nucleation processes of solid precipitates in mixed electrolytes also warrants reexamination of electrochemically induced phase separation and its pertinence for electrolyte engineering. Multistep nucleation mechanisms, akin to those described by De Yoreo and coworkers for biological carbonate precipitation,80 have much lower energy barriers than classical nucleation. It is possible that (bi)carbonate salts of [BMIm]+ undergo classical nucleation with high energy barriers, while [BisMImO]2+ salts tend to undergo multi-step nucleation pathways and are therefore more prone to precipitate formation. Thus, understanding how entropy governs double layer formation can offer fundamentally new ways to leverage entropy as an additional parameter for double layer engineering and can offer an exciting opportunity to reexamine precipitate nucleation and growth mechanisms.
Footnote |
† Electronic supplementary information (ESI) available: Additional experimental data. See DOI: https://doi.org/10.1039/d3sm01302b |
This journal is © The Royal Society of Chemistry 2024 |