The Electrostatic Screening Length in Concentrated Electrolytes

May 24, 2016 - This work was supported by The Leverhulme Trust (RPG-2015-328), the ERC (under Starting Grant LIQUISWITCH), and The John Fell Fund ...
0 downloads 0 Views 1MB Size
Letter pubs.acs.org/JPCL

The Electrostatic Screening Length in Concentrated Electrolytes Increases with Concentration Alexander M. Smith,*,† Alpha A. Lee,*,‡ and Susan Perkin*,† †

Department of Chemistry, Physical & Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States



S Supporting Information *

ABSTRACT: According to classical electrolyte theories interactions in dilute (low ion density) electrolytes decay exponentially with distance, with the Debye screening length the characteristic length scale. This decay length decreases monotonically with increasing ion concentration due to effective screening of charges over short distances. Thus, within the Debye model no long-range forces are expected in concentrated electrolytes. Here we reveal, using experimental detection of the interaction between two planar charged surfaces across a wide range of electrolytes, that beyond the dilute (Debye−Hückel) regime the screening length increases with increasing concentration. The screening lengths for all electrolytes studiedincluding aqueous NaCl solutions, ionic liquids diluted with propylene carbonate, and pure ionic liquidscollapse onto a single curve when scaled by the dielectric constant. This nonmonotonic variation of the screening length with concentration, and its generality across ionic liquids and aqueous salt solutions, demonstrates an important characteristic of concentrated electrolytes of substantial relevance from biology to energy storage.

E

densities allow the electric potential to decay over longer distances and so the electric field (gradient of potential) acts over greater distances from a charged particle or object. For example, a solution of NaCl in water at 25 °C and a concentration of 1 × 10−3 M has a Debye length of λD = 9.6 nm, corresponding to roughly 40 water diameters, whereas at a concentration of 0.1 M the Debye length is λD = 0.96 nm, equivalent to just four water diameters. At higher concentrations, when the Debye length becomes comparable to molecular dimensions, the Debye−Hückel theory is not expected to apply;5 ion size (steric or volume exclusion) effects and ion correlations come into play. Models to account for this include extensions based on the Poisson−Boltzmann equation,6,7 or analysis of the Ornstein−Zernike equation for binary ionic fluids,3,8 leading to the prediction of oscillatory (rather than monotonic) decay for concentrated electrolytes.9 Although it is well known within the theoretical community that the asymptotic decay length is determined by either charge−charge or hard core correlations and could exceed the Debye length in concentrated electrolytes,8,10 this has not been directly measured in electrolytes in the past, and its implications are not widely appreciated. We report here measurements of the interaction force between two smooth planar surfaces across electrolytes showing an exponentially decaying component to the interaction force extending over long-range even at high

lectrolytes are ubiquitous in nature and in technology: from the interior of cells to the oceans, from supercapacitors to nanoparticle dispersions, electrolytes act as both solvent and ion conduction medium. The “high concentration” end of the electrolyte spectrum is of particular relevance; examples of natural and technological electrolytes include human blood plasma (∼0.15 M), the Atlantic Ocean (∼0.6 M), supercapacitor electrolytes (1−2 M), and the Dead Sea (∼4.7 M). Despite this, the properties of concentrated electrolytes remain relatively uncharted territory, particularly from an experimental point of view, compared to the well-understood dilute electrolyte regime where Poisson−Boltzmann (meanfield) electrostatics applies and approximations such as in the Debye−Hückel theory are appropriate.1−4 Within the classical mean-field theories that explain screening and interactions in dilute electrolytes,1 the Debye length, λD, is the yardstick against which the reach of electrostatic interactions in solution is typically measured. The Debye length scales inversely with the square root of bulk ion density, ρ∞, according to ⎛ ε ε k T ⎞1/2 λD = ⎜⎜ r 0 B2 2 ⎟⎟ ⎝ ∑i ρ∞ i e zi ⎠

(1)

where ρ∞i is the number density of ion type i, εr is dielectric constant, ε0 is the permittivity of free space, T is the temperature, kB is the Boltzmann constant, and zi is the ion valency. Thus, within mean-field Poisson−Boltzmann models, a high concentration (or density) of ions leads to screening of an electrostatic potential over short distances, while low ion © 2016 American Chemical Society

Received: April 22, 2016 Accepted: May 24, 2016 Published: May 24, 2016 2157

DOI: 10.1021/acs.jpclett.6b00867 J. Phys. Chem. Lett. 2016, 7, 2157−2163

Letter

The Journal of Physical Chemistry Letters

Figure 1. Measured normal force, FN, normalized by radius of curvature, R, between crossed cylindrical mica surfaces across (A) pure [C4C1Pyrr][NTf2] ionic liquid and (B) 2 M NaCl(aq). Shown on log scale are the forces measured on approach of the surfaces, with exponential fits of the long-range force shown as solid lines. In each case the insets show, on linear scales, the combination of forces measured on approach (dots) and retraction of the surfaces (triangles), which reveals the nature of the oscillatory region at small distances.

we find that λexp increases again with increasing concentration, up to maximum values greater than the theoretical Debye length by a factor of ∼120. This trend is followed by all electrolytes studied including ionic liquids diluted with a polar solvent (propylene carbonate) and NaCl solutions in water. Scaling the measured decay length of interactions in each electrolyte or ionic liquid by its respective dielectric constant allows us to collapse the data onto a single curve, suggesting that these observations are independent of molecular architecture or details of molecular interactions and instead are a general feature of concentrated electrolytes with a common underlying mechanism. Measurements of interaction force between two atomically smooth and large-area mica sheets in crossed-cylinder configuration, across the electrolytes of interest, were carried

electrolyte concentrations; that is, the experimental decay length greatly exceeds the theoretical Debye length; λexp ≫ λD. This observation is in accordance with a recent observation of long-range screening in ionic liquids;11 we now provide measurements of a wide range of pure ionic liquids, ionic liquids dissolved in a polar solvent (propylene carbonate), and simple NaCl salt in water to demonstrate this to be a general observation of electrolytes at high concentration. We uncover a clear and strong nonmonotonic relationship between the exponential decay length of the observed longrange force, λexp, and the ion density (concentration), as follows. In dilute electrolytes the decay length decreases with concentration, as expected within the Debye−Hückel regime, reaching unmeasurably small λexp at ∼0.1M. Then, as the electrolyte concentration is increased further, beyond ∼0.5 M, 2158

DOI: 10.1021/acs.jpclett.6b00867 J. Phys. Chem. Lett. 2016, 7, 2157−2163

Letter

The Journal of Physical Chemistry Letters

Figure 2. Measured decay length, λexp, plotted as a function of c1/2, for aqueous NaCl solutions (filled circles), and solutions of an ionic liquid (IL), [C4C1Pyrr][NTf2], mixed with propylene carbonate (asterisk symbols) and compared with values in the literature13 (filled diamonds). The solid and dashed lines show how the theoretical Debye length varies with c1/2 for two values of dielectric constant, εr = 80.1 (solid line) and εr = 64 (dashed line). Details of the electrolytes, including the ionic liquid chemical structure, are in the Supporting Information.

out using a surface force balance (SFB).12 The apparatus and procedure is essentially similar to that used in early measurements of DLVO forces13 and experiments with pure ionic liquids.11,14 In brief: the back-silvered mica sheets create an interferometric cavity so that the absolute separation distance between them can be measured using white-light interferometry (measuring the fringes of equal chromatic order (FECO)). The electrolyte is held in a droplet between the crossedcylindrical lenses, such that the confined film is in equilibrium with a large excess of bulk liquid outside of the confined film. Interaction forces arising from the difference in pressure in the film and the bulk lead to the deflection of a spring upon which one of the mica sheets is mounted, and this is detected in the spectral fringes. All measurements were carried out at T = 294 K. Electrolytes (details in the Supporting Information) were injected between the two freshly cleaved mica sheets and forces, FN, were measured as a function of film thickness (or mica surface separation), D. In our previous works with ionic liquids we focused on the short-range oscillatory forces; here we show results of high-resolution measurements at long-range allowing us to observe both the oscillatory region and an exponentially decaying force that exists at distances beyond the oscillatory part of the interaction. Further details of the experimental procedure are in the Supporting Information. We begin by presenting and discussing the nature of the mica−mica interaction across concentrated electrolyte for two example cases: one pure ionic liquid (intrinsic concentration of the pure liquid = 3.3 M) and one solution of NaCl in water at

2.0 M. Figure 1A shows, for an example ionic liquid 1-butyl-1methylpyrrolidinium bis[(trifluoromethane)sulfonyl]imide, [C4C1Pyrr][NTf2], the total force law between negatively charged mica surfaces across the liquid as a function of the liquid film thickness. The data are plotted on a log−linear scale in the main panel, and the short-range part of the same data is plotted on a linear−linear scale in the inset. An oscillatory component to the force is apparent in the region 0−6 nm; this has been the focus of our previous studies,15 and we do not discuss it further here. Beyond the oscillatory region, from ∼6− 20 nm, we measure a monotonic and exponentially decaying repulsive force. The decay length of the exponential region is the “experimental decay length”, λexp, of the long-range force, in this case ∼8.4 nm, as discussed and compared with equivalent values for a range of electrolytes in the following discussion. Figure 1B shows a similar force law measured for a concentrated aqueous electrolyte: 2 M NaCl in water. Again, an oscillatory region is measured at short-range, in this case with an oscillatory period of ∼0.5 nm and extending to distances of 3 nm, followed by an exponentially decaying monotonic region from ∼3−6 nm with λexp = ∼ 1.1 nm. The forces within the oscillatory region, 0−3 nm, have been discussed in the past and are thought to be caused by squeezeout of sequential layers of hydrated (or partially hydrated) ions16 and will not be interpreted further here. We note, however, that detection of these oscillatory forces in NaCl is a crucial signature of cleanliness of the system; experiments where these were not observed were often found to be 2159

DOI: 10.1021/acs.jpclett.6b00867 J. Phys. Chem. Lett. 2016, 7, 2157−2163

Letter

The Journal of Physical Chemistry Letters irreproducible between different contact areas and were discarded due to contamination. The monotonic part, in the region 3−6 nm, has also been noted in studies of forces across other concentrated aqueous electrolytes in the past.17−19 We have conducted measurements similar to those in Figure 1 for a range of electrolytes at different concentrations, all measured at T = 294 K. In each case the asymptotic decay (manner in which the force decays toward zero at large D) was well represented by an exponentially decaying function; that is, for separations sufficiently far beyond the oscillatory regime F N /R ≈ exp(−D/λ exp ), where λ exp is a characteristic “experimental decay length” different for each electrolyte. In Figure 1 the values of λexp are the gradients of the black lines. We note here that the asymptotic decay length of the force between charged surfaces across a fluid is known to match the bulk asympototic decay in the same fluid,8 and so in this experiment we are probing in a systematic and high-resolution manner a property of bulk electrolyte. This decay length, λexp, can therefore be thought of as an experimentally determined screening length for interactions in the electrolyte. In the following sections we compare and discuss λexp for simple NaCl in water and an ionic liquid in propylene carbonate systematically varying the concentration in each caseand for a range of different pure ionic liquids. Figure 2 shows λexp plotted as a function of c1/2, for ease of comparison with the theoretical Debye length, λD, for ionic liquid diluted with propylene carbonate and NaCl salt dissolved in water. Also marked in Figure 2 are lines indicating the expected (theoretical) decay length, λD, for two values of εr; εr = 80.1 (appropriate for NaCl in water at low concentration) and εr = 64.0 (pure propylene carbonate). We note that here λD is calculated according to eq 1 with the assumption that all ions are free; that is, ρ∞ is the total concentration of ionic species; we do not propose that all ions are indeed “free” but simply that this is a suitable reference length scale with which to measure deviation from ideal (dilute electrolyte) behavior. It is clear by inspection of Figure 2 that λexp ≈ λD at low concentrations (c < 0.1 M; although this depends on dielectric constant as we show later), as has been demonstrated many times in the past,1 whereas at high concentrations there is strong deviation of λexp from λD. We make the following observations: (i) There appears to be a concentration of minimum λexp in the range cmin ≈ 0.1 to 0.5 M and (ii) above this concentration λexp increases with concentration, showing strong quantitative and qualitative deviation from the theoretical Debye length. The magnitude of the measured λexp is substantially lower for NaCl in water than for the propylene-carbonate electrolyte; this is accounted for by their different dielectric constants, as discussed later. We also studied pure ionic liquids and the resulting λexp values are shown in Figure 3. Variation in the intrinsic concentration for pure ionic liquids arises from differences in molecular volume and packing; details of the individual ionic liquids studied and their intrinsic concentrations are in the Supporting Information. In every case a long-range monotonic repulsion was observed with decay-length substantially greater than the theoretical value; that is, λexp ≫ λD (where λD ≈ 0.1 nm for these pure ionic liquids). Also plotted in Figure 3 are data points taken from other authors11,20,21 showing good agreement with respect to magnitudes of λexp. The scatter in values for pure ionic liquids in Figure 3 is at least partially attributed to their varying ion sizes and dielectric properties, as described later and shown in Figure 4.

Figure 3. Measured decay length, λexp, plotted as a function of molar ionic liquid concentration, c, for different pure ionic liquids studied here (circles) and compared with values in the literature11,20,21 (diamonds). The chemical details of each ionic liquid and the individual decay and concentration data corresponding to each are given in the Supporting Information.

When considering the origin of the longest range decay we arrive first at the question of the extent to which molecular (chemical) features of the ions play a role or whether the behavior can be explained purely in terms of the electrolyte concentration and dielectric constant. To illuminate this we now present, in Figure 4, dielectric-scaled values of λexp for aqueous NaCl, ionic liquid solutions in propylene carbonate, and pure ionic liquids (those where literature values of εr are available22,23). For aqueous NaCl solutions, εr depends strongly on concentration in the range 80.1 (pure water) to ∼42 (5 M NaCl).24,25 For mixtures of ionic liquid with propylene carbonate, εr was calculated using effective medium theory.26 It is apparent in Figure 4 that the measured decay lengths in all electrolytes studied appear to collapse onto a single curve, providing a unique relationship between the deviation from the Debye length, λexp/λD, and the dielectric-scaled ion density, as captured by the parameter d/λD (where d is the mean ion diameter in the electrolyte and so is constant for each electrolyte). The quantity d/λD can be interpreted as a dielectric-scaled concentration because 1/λD is proportional to c1/2/ε1/2. The dielectric constant depends on the background solvent and the polarizability of the ions and thus could be different even for ionic solutions with the same ion concentration and ion diameter. Therefore, for any particular salt and solvent combination d/λD depends mainly on concentration but in a way that is modulated by the dielectric constant, which also varies as the ion concentration increases. This useful dimensionless quantity allows comparison of measured length scales in electrolytes with very different dielectric constants. Additionally, there is an implicit prediction of how temperature should affect the results because it is accounted for in λD. Mean ion diameter, d, is estimated from the cube root of the volume per ion pair halved, which gives values in good agreement with approximately half of the measured spacing in the oscillatory region of the force profiles for ionic liquids. Our method agrees well if the force oscillations arise from combined cation and anion layers, as is currently accepted to be the structure in ionic liquids at charged surfaces. A similar value of d can also be inferred from the cation−anion correlation peak in X-ray scattering experiments with the same ionic liquid.27 This is likely to be the most accurate; however, there is a limited number of ionic liquids studied in the literature. We note that alternative methods could be used to 2160

DOI: 10.1021/acs.jpclett.6b00867 J. Phys. Chem. Lett. 2016, 7, 2157−2163

Letter

The Journal of Physical Chemistry Letters

Figure 4. Deviation of the experimental decay length from the Debye length, λexp/λD, plotted versus d/λD for pure ionic liquids (open diamonds), aqueous NaCl solutions (filled circles), and solutions of an ionic liquid, [C4C1Pyrr][NTf2], mixed with propylene carbonate (asterisk symbols). Literature values17,19 are also included for aqueous LiCl, KCl, and CsCl solutions (open circles), the details of which are in the Supporting Information. Data are plotted for ionic liquids for which values of εr are in the literature and are in the range 12−16. Values of εr for ionic liquid solutions in propylene carbonate were calculated using effective medium theory,26 while for aqueous NaCl solutions εr varies between 80.1 (pure water) and ∼42 (5 M NaCl).24 The inset shows the same data on a double logarithmic plot. The double logarithmic plot demonstrates the apparent power-law dependence in the range 1 < d/λD < 4. All quantitative data and molecular structures relating to this figure are detailed in the Supporting Information.

estimate or calculate the mean ion diameter, d, for example, taking into account the packing fraction of ions, and it is not immediately clear whether the “bare” or hydrated diameter should be taken for aqueous NaCl electrolytes. It will be important in future work to compare in detail the impact of different methods for assignment of d for both ionic liquids and simple electrolytes in a systematic way. All quantities used in preparing Figure 4 are provided in the Supporting Information. The fact that the screening length of ionic liquids and simple electrolytes appears to collapse onto a single curve (Figure 4) suggests a striking degree of similarity between ionic liquids and concentrated salt solutions; it appears that microphase separation and bulk nanostructuring that occurs in many ionic liquids due to interactions between nonpolar domains,28−33 or effects of ion size and shape anisotropy, do not alter the qualitative behavior of the interionic interactions. Although we do not intend to propose a theoretical model for the molecular origin of our observed trend in asymptotic decay with concentration here, we note that our observations are likely to be related to the works of Evans8 and Attard10 in the 1990s, where the mean spherical approximation (MSA) and hypernetted chain approximation (HNC) were used to study the transitions from monotonic to charge-oscillatory to densityoscillatory regimes with increasing electrolyte concentration. In

each of these works a qualitatively very similar nonmonotonic trend of λexp with d/λD was predicted, with the minimum in decay length occurring at the Kirkwood line, the point at which the charge−charge correlation function transitions from monotonic to oscillatory. The Kirkwood line is predicted to be at d/λD ≈ √2; close but not identical to the point where the electrostatic screening length begins to increase with concentration in Figure 4; however, we note substantial quantitative difference between the magnitude of the screening lengths measured here and the predictions of MSA and HNC. Our experiments are carried out using two atomically smooth and macroscopic mica sheets as the substrates for force measurement across electrolytes, and it is worth commenting upon the generality or otherwise of forces measured between mica sheets. Mica is negatively charged in electrolyte solutions,34 causing a perturbation to liquid structure in the near-surface region. This, in turn, causes variations in pressure in the thin film between the surfaces relative to the bulk liquid, as detected in our force measurement. The magnitude of the force and details of the oscillatory short-range part of the force are therefore likely to be altered in their details when the substrate is changed; however, the screening length of the monotonic part of the interaction, λexp, is expected to be a 2161

DOI: 10.1021/acs.jpclett.6b00867 J. Phys. Chem. Lett. 2016, 7, 2157−2163

The Journal of Physical Chemistry Letters



f unction of the bulk liquid only, as is the case for the Debye length, and thus is insensitive to the nature of the surface. The present experiments encourage comparison with several recent theoretical and experimental works. First, it is likely that the recently reported voltage-induced phase transition in ionic liquid−metal interfacial capacitance35 and maximum in differential capacitance as a function of ionic liquid concentration36 could be a part of a more general phenomenon also observable in concentrated electrolytes. Although this phase transition and anomaly in capacitance can be captured by positing a nonconvex bulk free energy37 or a charge-frustrated Ising model36 that models interactions specific to ionic liquids, it remains to be seen whether there is a more general theoretical description for this phenomenon. Second, the demonstration of general nonmonotonic dependence of decay length on electrolyte concentration should contribute toward interpreting the long-range forces detected in ionic liquids;11 a theory for long-range forces in ionic liquids must be consistent with our observed increase in screening length with concentration over a certain range. On the other hand, theories relating to the interplay between bulk and surface nanostructure38 need to keep in mind the observation that similar long-range forces occur in NaCl solutions where nanostructuring due to solvophobic interactions is not expected to occur; however, bulk oscillatory structuring is nonetheless expected in simple electrolytes, such as NaCl, at sufficiently high concentrations.39,40 Our direct measurements of the force between macroscopic and atomically smooth mica sheets across electrolytes have revealed a clear nonmonotonic dependence of decay length on concentration of the electrolyte. Within the Debye−Hückel region, below ca. 0.1 M, the screening length decreases with concentration according to the Debye theory, reaching a minimum when the ion diameter matches the Debye length; d/ λD ≈ 1. At higher concentrations, in the range 1 < d/λD < 5, the experimental decay length increases strongly with increasing concentration up to (λexp/λD) ≈ 120. Deviation of the experimental decay length from the theoretical Debye length, λexp/λD, is described by the single dimensionless parameter, d/ λD, for all electrolytes and ionic liquids studied, indicating a common mechanism for 1:1 electrolytes. Such substantial deviation from the canonical understanding of electrolytes is likely to be important for understanding electric fields and surface forces in concentrated electrolytes in a multitude of scenarios. These include intermolecular forces between interacting proteins and membranes in confined geometries, where ion concentration can be high, and in energy storage applications, such as batteries and supercapacitors, where highly concentrated electrolytes are confined within porous electrodes.



Letter

AUTHOR INFORMATION

Corresponding Authors

*A.M.S.: E-mail: [email protected]. *A.A.L.: E-mail: [email protected]. *S.P.: E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by The Leverhulme Trust (RPG2015-328), the ERC (under Starting Grant LIQUISWITCH), and The John Fell Fund (Oxford University). A.M.S. was supported by a Doctoral Prize from the EPSRC. A.A.L. was supported by a Fulbright Fellowship to Harvard University. We are grateful to Lorna Swift and Tom Welton for synthesis of a phosphonium ionic liquid. We thank Bob Evans for illuminating discussions. A.A.L. appreciates fruitful discussions with David Limmer.



REFERENCES

(1) Evans, D. F.; Wennerström, H. The Colloidal Domain: Where Physics, Chemistry, Biology, and Technology Meet; Wiley, 1999. (2) Levin, Y. Electrostatic Correlations: From Plasma to Biology. Rep. Prog. Phys. 2002, 65, 1577. (3) Jing, Y.; Jadhao, V.; Zwanikken, J. W.; Olvera de la Cruz, M. Ionic Structure in Liquids Confined by Dielectric Interfaces. J. Chem. Phys. 2015, 143, 194508. (4) Yochelis, A.; Singh, M. B.; Visoly-Fisher, I. Coupling Bulk and near-Electrode Interfacial Nanostructuring in Ionic Liquids. Chem. Mater. 2015, 27, 4169−4179. (5) Bockris, J. O. M.; Reddy, A. K. N. Modern Electrochemistry; Springer, 1998; Vol. 1 (6) Onsager, L. Theories of Concentrated Electrolytes. Chem. Rev. 1933, 13, 73−89. (7) Kjellander, R.; Mitchell, D. J. An Exact but Linear and Poisson Boltzmann-Like Theory for Electrolytes and Colloid Dispersions in the Primitive Model. Chem. Phys. Lett. 1992, 200, 76−82. (8) Leote de Carvalho, R. J. F.; Evans, R. The Decay of Correlations in Ionic Fluids. Mol. Phys. 1994, 83, 619−654. (9) Stillinger, F. H.; Lovett, R. Ion-Pair Theory of Concentrated Electrolytes. I. Basic Concepts. J. Chem. Phys. 1968, 48, 3858−3868. (10) Attard, P. Asymptotic Analysis of Primitive Model Electrolytes and the Electrical Double Layer. Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top. 1993, 48, 3604−3621. (11) Gebbie, M. A.; Dobbs, H. A.; Valtiner, M.; Israelachvili, J. N. Long-Range Electrostatic Screening in Ionic Liquids. Proc. Natl. Acad. Sci. U. S. A. 2015, 112, 7432−7437. (12) Smith, A. M.; Lovelock, K. R. J.; Perkin, S. Monolayer and Bilayer Structures in Ionic Liquids and Their Mixtures Confined to Nano-Films. Faraday Discuss. 2014, 167, 279−292. (13) Pashley, R. M. Dlvo and Hydration Forces between Mica Surfaces in Li+,Na+,K+, and Cs+ Electrolyte Solutions: - a Correlation of Double-Layer and Hydration Forces with Surface Cation-Exchange Properties. J. Colloid Interface Sci. 1981, 83, 531−546. (14) Perkin, S.; Crowhurst, L.; Niedermeyer, H.; Welton, T.; Smith, A. M.; Gosvami, N. N. Self-Assembly in the Electrical Double Layer of Ionic Liquids. Chem. Commun. 2011, 47, 6572−6574. (15) Smith, A. M.; Lovelock, K. R. J.; Gosvami, N. N.; Licence, P.; Dolan, A.; Welton, T.; Perkin, S. Monolayer to Bilayer Structural Transition in Confined Pyrrolidinium-Based Ionic Liquids. J. Phys. Chem. Lett. 2013, 4, 378−382. (16) Espinosa-Marzal, R. M.; Drobek, T.; Balmer, T.; Heuberger, M. P. Hydrated-Ion Ordering in Electrical Double Layers. Phys. Chem. Chem. Phys. 2012, 14, 6085−6093.

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpclett.6b00867. Details of the materials and methods; details of the ionic liquids used and the data points used to plot Figure 3; quantities used to perform scaling in Figure 4 and the final data points in numerical form used to plot Figure 4, including the literature values; discussion of alternative calculations of d and figures demonstrating how this impacts d/λD. (PDF) 2162

DOI: 10.1021/acs.jpclett.6b00867 J. Phys. Chem. Lett. 2016, 7, 2157−2163

Letter

The Journal of Physical Chemistry Letters

(37) Limmer, D. T. Interfacial Ordering and Accompanying Divergent Capacitance at Ionic Liquid-Metal Interfaces. Phys. Rev. Lett. 2015, 115, 256102. (38) Li, H.; Endres, F.; Atkin, R. Effect of Alkyl Chain Length and Anion Species on the Interfacial Nanostructure of Ionic Liquids at the Au(111)-Ionic Liquid Interface as a Function of Potential. Phys. Chem. Chem. Phys. 2013, 15, 14624−14633. (39) Kirkwood, J. G. Statistical Mechanics of Liquid Solutions. Chem. Rev. 1936, 19, 275−307. (40) Gavish, N.; Yochelis, A. Theory of Phase Separation and Polarization for Pure Ionic Liquids. J. Phys. Chem. Lett. 2016, 7, 1121− 1126.

(17) Pashley, R. M.; Israelachvili, J. N. Molecular Layering of Water in Thin Films between Mica Surfaces and Its Relation to Hydration Forces. J. Colloid Interface Sci. 1984, 101, 511−523. (18) Pashley, R. M. Forces between Mica Surfaces in La3+ and Cr3+ Electrolyte Solutions. J. Colloid Interface Sci. 1984, 102, 23−35. (19) Baimpos, T.; Shrestha, B. R.; Raman, S.; Valtiner, M. Effect of Interfacial Ion Structuring on Range and Magnitude of Electric Double Layer, Hydration, and Adhesive Interactions between Mica Surfaces in 0.05−3 M Li+ and Cs+ Electrolyte Solutions. Langmuir 2014, 30, 4322−4332. (20) Espinosa-Marzal, R. M.; Arcifa, A.; Rossi, A.; Spencer, N. D. Microslips to “Avalanches” in Confined, Molecular Layers of Ionic Liquids. J. Phys. Chem. Lett. 2014, 5, 179−184. (21) Espinosa-Marzal, R. M.; Arcifa, A.; Rossi, A.; Spencer, N. D. Ionic Liquids Confined in Hydrophilic Nanocontacts: Structure and Lubricity in the Presence of Water. J. Phys. Chem. C 2014, 118, 6491− 6503. (22) Wakai, C.; Oleinikova, A.; Ott, M.; Weingärtner, H. How Polar Are Ionic Liquids? Determination of the Static Dielectric Constant of an Imidazolium-Based Ionic Liquid by Microwave Dielectric Spectroscopy. J. Phys. Chem. B 2005, 109, 17028−17030. (23) Huang, M.-M.; Jiang, Y.; Sasisanker, P.; Driver, G. W.; Weingärtner, H. Static Relative Dielectric Permittivities of Ionic Liquids at 25 °C. J. Chem. Eng. Data 2011, 56, 1494−1499. (24) Buchner, R.; Hefter, G. T.; May, P. M. Dielectric Relaxation of Aqueous NaCl Solutions. J. Phys. Chem. A 1999, 103, 1−9. (25) Gavish, N.; Promislow, K. Dependence of the Dielectric Constant of Electrolyte Solutions on Ionic Concentration. arXiv preprint arXiv:1208.5169, 2012. (26) Bergman, D. J. The Dielectric Constant of a Composite Materiala Problem in Classical Physics. Phys. Rep. 1978, 43, 377− 407. (27) Santos, C. S.; Murthy, N. S.; Baker, G. A.; Castner, E. W. Communication: X-Ray Scattering from Ionic Liquids with Pyrrolidinium Cations. J. Chem. Phys. 2011, 134, 121101. (28) Urahata, S. M.; Ribeiro, M. C. C. Structure of Ionic Liquids of 1Alkyl-3-Methylimidazolium Cations: A Systematic Computer Simulation Study. J. Chem. Phys. 2004, 120, 1855−1863. (29) Wang, Y. T.; Voth, G. A. Unique Spatial Heterogeneity in Ionic Liquids. J. Am. Chem. Soc. 2005, 127, 12192−12193. (30) Canongia Lopes, J. N. A.; Pádua, A. A. H. Nanostructural Organization in Ionic Liquids. J. Phys. Chem. B 2006, 110, 3330−3335. (31) Triolo, A.; Russina, O.; Bleif, H.-J.; Di Cola, E. Nanoscale Segregation in Room Temperature Ionic Liquids. J. Phys. Chem. B 2007, 111, 4641−4644. (32) Russina, O.; Triolo, A.; Gontrani, L.; Caminiti, R.; Xiao, D.; Hines, L. G., Jr.; Bartsch, R. A.; Quitevis, E. L.; Pleckhova, N.; Seddon, K. R. Morphology and Intermolecular Dynamics of 1-Alkyl-3Methylimidazolium Bis{(Trifluoromethane)Sulfonyl}Amide Ionic Liquids: Structural and Dynamic Evidence of Nanoscale Segregation. J. Phys.: Condens. Matter 2009, 21, 424121. (33) Shimizu, K.; Costa Gomes, M. F.; Pádua, A. A.; Rebelo, L. P.; Canongia Lopes, J. N. Three Commentaries on the Nano-Segregated Structure of Ionic Liquids. J. Mol. Struct.: THEOCHEM 2010, 946, 70−76. (34) Perkin, S.; Albrecht, T.; Klein, J. Layering and Shear Properties of an Ionic Liquid, 1-Ethyl-3-Methylimidazolium Ethylsulfate, Confined to Nano-Films between Mica Surfaces. Phys. Chem. Chem. Phys. 2010, 12, 1243−1247. (35) Merlet, C.; Limmer, D. T.; Salanne, M.; van Roij, R.; Madden, P. A.; Chandler, D.; Rotenberg, B. The Electric Double Layer Has a Life of Its Own. J. Phys. Chem. C 2014, 118, 18291−18298. (36) Bozym, D. J.; Uralcan, B.; Limmer, D. T.; Pope, M. A.; Szamreta, N. J.; Debenedetti, P. G.; Aksay, I. A. Anomalous Capacitance Maximum of the Glassy Carbon−Ionic Liquid Interface through Dilution with Organic Solvents. J. Phys. Chem. Lett. 2015, 6, 2644−2648. 2163

DOI: 10.1021/acs.jpclett.6b00867 J. Phys. Chem. Lett. 2016, 7, 2157−2163