The one-dimensional Coulomb lattice fluid capacitor

Theory and Simulations of Ionic Liquids in Nanoconfinement

Room-temperature ionic liquids (RTILs) have exciting properties such as nonvolatility, large electrochemical windows, and remarkable variety, drawing much interest in energy storage, gating, electrocatalysis, tunable lubrication, and other applications. Confined RTILs appear in various situations, for instance, in pores of nanostructured electrodes of supercapacitors and batteries, as such electrodes increase the contact area with RTILs and enhance the total capacitance and stored energy, between crossed cylinders in surface force balance experiments, between a tip and a sample in atomic force microscopy, and between sliding surfaces in tribology experiments, where RTILs act as lubricants. The properties and functioning of RTILs in confinement, especially nanoconfinement, result in fascinating structural and dynamic phenomena, including layering, overscreening and crowding, nanoscale capillary freezing, quantized and electrotunable friction, and superionic state. This review offers a comprehensive analysis of the fundamental physical phenomena controlling the properties of such systems and the current state-of-the-art theoretical and simulation approaches developed for their description. We discuss these approaches sequentially by increasing atomistic complexity, paying particular attention to new physical phenomena emerging in nanoscale confinement. This review covers theoretical models, most of which are based on mapping the problems on pertinent statistical mechanics models with exact analytical solutions, allowing systematic analysis and new physical insights to develop more easily. We also describe a classical density functional theory, which offers a reliable and computationally inexpensive tool to account for some microscopic details and correlations that simplified models often fail to consider. Molecular simulations play a vital role in studying confined ionic liquids, enabling deep microscopic insights otherwise unavailable to researchers. We describe the basics of various simulation approaches and discuss their challenges and applicability to specific problems, focusing on RTIL structure in cylindrical and slit confinement and how it relates to friction and capacitive and dynamic properties of confined ions.

Capacitive energy storage in single-file pores: Exactly solvable models and simulations

Understanding charge storage in low-dimensional electrodes is crucial for developing novel ecologically friendly devices for capacitive energy storage and conversion and water desalination. Exactly solvable models allow in-depth analyses and essential physical insights into the charging mechanisms. So far, however, such analytical approaches have been mainly limited to lattice models. Herein, we develop a versatile, exactly solvable, one-dimensional off-lattice model for charging single-file pores. Unlike the lattice model, this model shows an excellent quantitative agreement with three-dimensional Monte Carlo simulations. With analytical calculations and simulations, we show that the differential capacitance can be bell-shaped (one peak), camel-shaped (two peaks), or have four peaks. Transformations between these capacitance shapes can be induced by changing pore ionophilicity, by changing cation-anion size asymmetry, or by adding solvent. We find that the camel-shaped capacitance, characteristic of dilute electrolytes, appears for strongly ionophilic pores with high ion densities, which we relate to charging mechanisms specific to narrow pores. We also derive a large-voltage asymptotic expression for the capacitance, showing that the capacitance decays to zero as the inverse square of the voltage, C ∼ u^-2. This dependence follows from hard-core interactions and is not captured by the lattice model.

One-dimensional colloidal model with dielectric inhomogeneity

We consider a one-dimensional model allowing analytical derivation of the effective interactions between two charged colloids. We evaluate exactly the partition function for an electroneutral salt-free suspension with dielectric jumps at the colloids' position. We derive a contact relation with the pressure that shows there is like-charge attraction, whether or not the counterions are confined between the colloids. In contrast to the homogeneous dielectric case, there is the possibility for the colloids to attract despite the number of counterions (N) being even. The results are shown to recover the mean-field prediction in the limit N→∞.

Superionic Liquids in Conducting Nanoslits: Insights from Theory and Simulations

Mapping the theory of charging supercapacitors with nanostructured electrodes on known lattice models of statistical physics is an interesting task, aimed at revealing generic features of capacitive energy storage in such systems. The main advantage of this approach is the possibility to obtain analytical solutions that allow new physical insights to be more easily developed. But how general the predictions of such theories could be? How sensitive are they to the choice of the lattice? Herein, we address these questions in relation to our previous description of such systems using the Bethe-lattice approach and Monte Carlo simulations. Remarkably, we find a surprisingly good agreement between the analytical theory and simulations. In addition, we reveal a striking correlation between the ability to store energy and ion ordering inside a pore, suggesting that such ordering can be beneficial for energy storage.

Dynamics of Ion Channels via Non-Hermitian Quantum Mechanics

We study dynamics and thermodynamics of ion transport in narrow, water-filled channels, considered as effective 1D Coulomb systems. The long range nature of the inter-ion interactions comes about due to the dielectric constants mismatch between the water and the surrounding medium, confining the electric filed to stay mostly within the water-filled channel. Statistical mechanics of such Coulomb systems is dominated by entropic effects which may be accurately accounted for by mapping onto an effective quantum mechanics. In presence of multivalent ions the corresponding quantum mechanics appears to be non-Hermitian. In this review we discuss a framework for semiclassical calculations for the effective non-Hermitian Hamiltonians. Non-Hermiticity elevates WKB action integrals from the real line to closed cycles on a complex Riemann surfaces where direct calculations are not attainable. We circumvent this issue by applying tools from algebraic topology, such as the Picard-Fuchs equation. We discuss how its solutions relate to the thermodynamics and correlation functions of multivalent solutions within narrow, water-filled channels.

Charge oscillations in ionic liquids: A microscopic cluster model

In spite of their enormous applications as alternative energy storage devices and lubricants, room-temperature ionic liquids (ILs) still pose many challenges from a pure scientific viewpoint. We develop an IL microscopic theory in terms of ionic clusters, which describes the IL behavior close to charged interfaces. The full structure factor of finite-size clusters is considered and allows us to retain fine and essential details of the system as a whole. Beside the reduction in the screening, it is shown that ionic clusters cause the charge density to oscillate near charged boundaries, with alternating ion-size thick layers, in agreement with experiments. We distinguish between short-range oscillations that persist for a few ionic layers close to the boundary, as opposed to long-range damped oscillations that hold throughout the bulk. The former can be captured by finite-size ion pairs, while the latter is associated with larger clusters with a pronounced quadrupole (or higher) moment. The long-wavelength limit of our theory recovers the well-known Bazant-Storey-Kornyshev (BSK) equation in the linear regime, and elucidates the microscopic origin of the BSK phenomenological parameters.

Superionic liquids in conducting nanoslits: A variety of phase transitions and ensuing charging behavior

We develop a theory of charge storage in ultranarrow slitlike pores of nanostructured electrodes. Our analysis is based on the Blume-Capel model in an external field, which we solve analytically on a Bethe lattice. The obtained solutions allow us to explore the complete phase diagram of confined ionic liquids in terms of the key parameters characterizing the system, such as pore ionophilicity, interionic interaction energy, and voltage. The phase diagram includes the lines of first- and second-order, direct and re-entrant phase transitions, which are manifested by singularities in the corresponding capacitance-voltage plots. Testing our predictions experimentally requires monodisperse, conducting ultranarrow slit pores, to permit only one layer of ions, and thick pore walls, to prevent interionic interactions across the pore walls. However, some qualitative features, which distinguish the behavior of ionophilic and ionophobic pores and their underlying physics, may emerge in future experimental studies of more complex electrode structures.

Ionic Coulomb blockade as a fractional Wien effect

Recent advances in nanofluidics have allowed the exploration of ion transport down to molecular-scale confinement, yet artificial porins are still far from reaching the advanced functionalities of biological ion machinery. Achieving single ion transport that is tunable by an external gate-the ionic analogue of electronic Coulomb blockade-would open new avenues in this quest. However, an understanding of ionic Coulomb blockade beyond the electronic analogy is still lacking. Here, we show that the many-body dynamics of ions in a charged nanochannel result in quantized and strongly nonlinear ionic transport, in full agreement with molecular simulations. We find that ionic Coulomb blockade occurs when, upon sufficient confinement, oppositely charged ions form 'Bjerrum pairs', and the conduction proceeds through a mechanism reminiscent of Onsager's Wien effect. Our findings open the way to novel nanofluidic functionalities, such as an ion pump based on ionic Coulomb blockade, inspired by its electronic counterpart.

Phase behaviour and structure of a superionic liquid in nonpolarized nanoconfinement

The ion-ion interactions become exponentially screened for ions confined in ultranarrow metallic pores. To study the phase behaviour of an assembly of such ions, called a superionic liquid, we develop a statistical theory formulated on bipartite lattices, which allows an analytical solution within the Bethe-lattice approach. Our solution predicts the existence of ordered and disordered phases in which ions form a crystal-like structure and a homogeneous mixture, respectively. The transition between these two phases can potentially be first or second order, depending on the ion diameter, degree of confinement and pore ionophobicity. We supplement our analytical results by three-dimensional off-lattice Monte Carlo simulations of an ionic liquid in slit nanopores. The simulations predict formation of ionic clusters and ordered snake-like patterns, leading to characteristic close-standing peaks in the cation-cation and anion-anion radial distribution functions.

Coarse-grained simulations of an ionic liquid-based capacitor: I. Density, ion size, and valency effects

We introduce a hierarchy of generic coarse-grained models of ionic liquids of increasing complexity. We use them in molecular dynamics simulations to study the differential capacitance of a capacitor consisting of an ionic liquid between two planar electrodes. The primary goal is to explain the complex dependence of the differential capacitance Cd on the electrode potential U in simple terms, e.g. in terms of the size and valency of the ions. For this purpose we introduce the symmetric model A, which qualitatively reproduces the Cd(U) dependence predicted by the mean-field theory but also reveals strong quantitative deviations. We further introduce size asymmetry in model A by increasing the cation size. In model B we vary the cation valency, keeping the sizes of both ions constant. We show that simultaneous increases in size and valency may compensate for each other, leading to a Cd(U) very similar to that for the symmetric case. We interpret distinct features in Cd(U) on the basis of the density profiles of the ions and charge density profiles. We focus on the first two ion layers at the electrode, and demonstrate that the polarization of the ionic liquid proceeds through replacement of one ion type by the other, in contrast to the simple increase in ion concentrations typical for dilute systems. The understanding gained for the simple models serves as a reference for interpretation of complex effects of ion size, valency and shape. This is carried through in part II (a separate article) where we show how the planar shape of ions in model C brings new features to the Cd(U) curve and also to the polarization mechanism.