Criticality and phase behavior in the restricted-primitive model electrolyte: Description of ion association

Response to "Comment on 'Binding Debye-Hückel theory for associative electrolyte solutions'" [J. Chem. Phys. 159, 154503 (2023)]

This Response addresses critiques raised about the Binding Debye-Hückel (BiDH) theory [Naseri Boroujeni et al., J. Chem. Phys. 159, 154503 (2023)] by Simonin and Bernard [J.-P. Simonin and O. Bernard, J. Chem. Phys. (2024)]. The critiques questioned the foundational framework of the Debye-Hückel (DH) theory, the relevance of ion pairing in primitive model fluids, and the accuracy of the BiDH model compared to mean spherical approximation model. Through a systematic rebuttal, supported by extensive literature review and comparison with Monte Carlo simulation data, this Response addresses these concerns. It demonstrates the efficacy of DH theory in describing real electrolyte solutions, validates the relevance of ion pairing in primitive model fluids, and establishes the BiDH model's accuracy in describing electrolyte properties.

Structure and Screening in Confined Electrolytes: The Role of Ion Association

The effect of ionic association on the structure and property of confined electrolytes is investigated using the classical density functional theory. We find that ionic association strongly affects the ion distribution, surface force, and screening behavior of confined electrolytes. The decay length ξ, which can describe the screening effect of high-concentration electrolytes, satisfies a scaling relationship ξ/λD ∼ (σ/λD)n, with λD being the Debye length and σ representing the ion diameter. We find that n = 1.5 in the nonassociation model, which is contributed by the charge correlation, but n = 3 in the association model, which is contributed by the density correlation. The ion association changes the concentration-dependent characteristics of the screening length by promoting the shift of the decay behavior from the charge-dominated regime to the density-dominated regime. Our result reveals the importance of ion association for electrolyte structure and screening behaviors.

Wettability of ultra-small pores of carbon electrodes by size-asymmetric ionic fluids

Recently, we studied the phase behavior of ionic fluids under confinement using the classical density functional theory within the framework of the restricted primitive model. The theoretical results indicate that narrowing the pore size may lead to a drastic reduction in the electric double layer capacitance, while increasing the surface electrical potential would improve the ionic accessibility of micropores. In this work, we extend the theoretical investigation to systems containing size-asymmetric electrolytes that may exhibit a vapor-liquid like phase transition in the bulk phase. The effects of pore size and surface electric potential on the phase diagram and microscopic structures of the confined electrolytes were studied over a broad range of parameters. We found that decreasing the pore size or increasing the surface potential could destabilize the liquid phase in micropores, and capillary evaporation could occur regardless of the size asymmetry between cations and anions. Compared to that in a symmetric ionic system, the vapor-liquid phase separation is more likely to take place as the size asymmetry becomes more pronounced. The phase transition would alter the "accessibility" of ions to micropores and lead to coexisting micropores with different surface charge densities as identified by Monte Carlo simulation.

Does capillary evaporation limit the accessibility of nonaqueous electrolytes to the ultrasmall pores of carbon electrodes?

Porous carbons have been widely utilized as electrode materials for capacitive energy storage. Whereas the importance of pore size and geometry on the device performance has been well recognized, little guidance is available for identification of carbon materials with ideal porous structures. In this work, we study the phase behavior of ionic fluids in slit pores using the classical density functional theory. Within the framework of the restricted primitive model for nonaqueous electrolytes, we demonstrate that the accessibility of micropores depends not only on the ionic diameters (or desolvation) but also on their wetting behavior intrinsically related to the vapor-liquid or liquid-liquid phase separation of the bulk ionic systems. Narrowing the pore size from several tens of nanometers to subnanometers may lead to a drastic reduction in the capacitance due to capillary evaporation. The wettability of micropores deteriorates as the pore size is reduced but can be noticeably improved by raising the surface electrical potential. The theoretical results provide fresh insights into the properties of confined ionic systems beyond electric double layer models commonly employed for rational design/selection of electrolytes and electrode materials.

Numerical study of density functional theory with mean spherical approximation for ionic condensation in highly charged confined electrolytes

We investigate numerically a density functional theory (DFT) for strongly confined ionic solutions in the canonical ensemble by comparing predictions of ionic concentration profiles and pressure for the double-layer configuration to those obtained with Monte Carlo (MC) simulations and the simpler Poisson-Boltzmann (PB) approach. The DFT consists of a bulk (ion-ion) and an ion-solid part. The bulk part includes nonideal terms accounting for long-range electrostatic and short-range steric correlations between ions and is evaluated with the mean spherical approximation and the local density approximation. The ion-solid part treats the ion-solid interactions at the mean-field level through the solution of a Poisson problem. The main findings are that ionic concentration profiles are generally better described by PB than by DFT, although DFT captures the nonmonotone co-ion profile missed by PB. Instead, DFT yields more accurate pressure predictions than PB, showing in particular that nonideal effects are important to describe highly confined ionic solutions. Finally, we present a numerical methodology capable of handling nonconvex minimization problems so as to explore DFT predictions when the reduced temperature falls below the critical temperature.

Multiply associating electrolytes in the binding mean spherical approximation: thermodynamic properties and speciation

Ionic solutions exhibiting multiple association are described within the binding mean spherical approximation (BiMSA). This model is based on the Wertheim formalism, in the framework of the primitive model at the McMillan-Mayer level. The cation and the anion form the various complexes according to stepwise complexation-equilibria. Analytic expressions for the Helmholtz energy, the internal energy, the speciation, and for the osmotic and activity coefficients are given considering a binary solution with an arbitrary number of association sites on one type of ion (polyion) and one site on the ions of opposite sign (counterions). As an alternative, mean field expressions, as developed in SAFT-type theories, are also presented. The result obtained from the latter approximate method exhibits a reasonable agreement with those from BiMSA for the speciation, and a remarkable one for the osmotic coefficient.

Monte Carlo simulation of the nonadditive restricted primitive model of ionic fluids: phase diagram and clustering

We report an accurate Monte Carlo calculation of the phase diagram and clustering properties of the restricted primitive model with nonadditive hard-sphere diameters. At high density the positively nonadditive fluid shows more clustering than in the additive model and the negatively nonadditive fluid shows less clustering than in the additive model; at low density the reverse scenario appears. A negative nonadditivity tends to favor the formation of neutrally charged clusters starting from the dipole. A positive nonadditivity favors the pairing of like ions at high density. The critical point of the gas-liquid phase transition moves at higher temperatures and higher densities for a negative nonadditivity and at lower temperatures and lower densities for a positive nonadditivity. The law of corresponding states does not seem to hold strictly. Our results can be used to interpret recent experimental works on room temperature ionic liquids.

The one-dimensional Coulomb lattice fluid capacitor

The one-dimensional Coulomb lattice fluid in a capacitor configuration is studied. The model is formally exactly soluble via a transfer operator method within a field theoretic representation of the model. The only interactions present in the model are the one-dimensional Coulomb interaction between cations and anions and the steric interaction imposed by restricting the maximal occupancy at any lattice site to one particle. Despite the simplicity of the model, a wide range of intriguing physical phenomena arise, some of which are strongly reminiscent of those seen in experiments and numerical simulations of three-dimensional ionic liquid based capacitors. Notably, we find regimes where over-screening and density oscillations are seen near the capacitor plates. The capacitance is also shown to exhibit strong oscillations as a function of applied voltage. It is also shown that the corresponding mean-field theory misses most of these effects. The analytical results are confirmed by extensive numerical simulations.

Microscopic approach to the kinetics of pattern formation of charged molecules on surfaces

A microscopic formalism based on computing many-particle densities is applied to the analysis of the diffusion-controlled kinetics of pattern formation in oppositely charged molecules on surfaces or adsorbed at interfaces with competing long-range Coulomb and short-range Lennard-Jones interactions. Particular attention is paid to the proper molecular treatment of energetic interactions driving pattern formation in inhomogeneous systems. The reverse Monte Carlo method is used to visualize the spatial molecular distribution based on the calculated radial distribution functions (joint correlation functions). We show the formation of charge domains for certain combinations of temperature and dynamical interaction parameters. The charge segregation evolves into quasicrystalline clusters of charges, due to the competing long- and short-range interactions. The clusters initially co-exist with a gas phase of charges that eventually add to the clusters, generating "fingers" or line of charges of the same sign, very different than the nanopatterns expected by molecular dynamics in systems with competing interactions in two dimensions, such as strain or dipolar versus van der Waals interactions.

Temperature dependence of the double layer capacitance for the restricted primitive model of an electrolyte solution from a density functional approach

We apply a different version of the density functional theory, given by Pizio, Patrykiejew, and Sokolowski [J. Chem. Phys. 121, 11957 (2004)], for a nonuniform restricted primitive model of an electrolyte solution to evaluate the temperature dependence of the capacitance of an electric double layer. We show that this theory is capable of reproducing the computer simulation data at a quantitative level. In particular, the reversal of the temperature dependence of the capacitance at low temperatures is predicted. This phenomenon has been difficult to predict from theory. Further, this theory also leads to an accurate description of the double layer structure.

Size-asymmetric primitive model at low temperature: description of ion pairing and location of the critical point

We argue that Bjerrum's approach to ion pairing is inappropriate for the size-asymmetric primitive model in the neighborhood of its critical point, and propose a new approach based on the Stillinger-Lovett pairing procedure. The new approach recursively scales up the ion size until linear approximations are suitable for analyzing such a model. To locate the critical point, a residual van der Waals interaction between pairs is added, with an energy cutoff adjusted to match the critical temperature of the restricted primitive model. The locations and downward trends of T(c) and rho(c) with asymmetry are found to compare favorably with simulations.

Criticality in charge-asymmetric hard-sphere ionic fluids

Phase separation and criticality are analyzed in z:1 charge-asymmetric ionic fluids of equisized hard spheres by generalizing the Debye-Hückel approach combined with ionic association, cluster solvation by charged ions, and hard-core interactions, following lines developed by Fisher and Levin for the 1:1 case (i.e., the restricted primitive model). Explicit analytical calculations for 2:1 and 3:1 systems account for ionic association into dimers, trimers, and tetramers and subsequent multipolar cluster solvation. The reduced critical temperatures, Tc* (normalized by z), decrease with charge asymmetry, while the critical densities increase rapidly with . The results compare favorably with simulations and represent a distinct improvement over all current theories such as the mean spherical approximation, symmetric Poisson-Boltzmann theory, etc. For z not equal to 1, the interphase Galvani (or absolute electrostatic) potential difference, Deltaphi(T), between coexisting liquid and vapor phases is calculated and found to vanish as absolute value (T-Tc) beta when T-->Tc-with, since our approximations are classical, beta = (1/2). Above Tc, the compressibility maxima and so-called k-inflection loci (which aid the fast and accurate determination of the critical parameters) are found to exhibit a strong z dependence.

Phase coexistence in a polydisperse charged hard-sphere fluid: Polymer mean spherical approximation

We have reconsidered the phase behavior of a polydisperse mixture of charged hard spheres (CHSs) introducing the concept of minimal size neutral clusters. We thus take into account ionic association effects observed in charged systems close to the phase boundary where the properties of the system are dominated by the presence of neutral clusters while the amount of free ions or charged clusters is negligible. With this concept we clearly pass beyond the simple level of the mean spherical approximation (MSA) that we have presented in our recent study of a polydisperse mixture of CHS [Yu. V. Kalyuzhnyi, G. Kahl, and P. T. Cummings, J. Chem. Phys. 120, 10133 (2004)]. Restricting ourselves to a 1:1 and possibly size-asymmetric model we treat the resulting polydisperse mixture of neutral, polar dimers within the framework of the polymer MSA, i.e., a concept that--similar as the MSA--readily can be generalized from the case of a mixture with a finite number of components to the polydisperse case: again, the model belongs to the class of truncatable free-energy models so that we can map the formally infinitely many coexistence equations onto a finite set of coupled, nonlinear equations in the generalized moments of the distribution function that characterizes the system. This allows us to determine the full phase diagram (in terms of binodals as well as cloud and shadow curves), we can study fractionation effects on the level of the distribution functions of the coexisting daughter phases, and we propose estimates on how the location of the critical point might vary in a polydisperse mixture with an increasing size asymmetry and polydispersity.

Phase behavior of the restricted primitive model of ionic fluids with association in slitlike pores. Density-functional approach

We present results of investigations of condensation of restricted primitive model of electrolyte solutions with association between oppositely charged ions confined to slitlike pores. The associative interaction leads to the formation of ionic pairs. It is accounted for by incorporating the first-order thermodynamic perturbation theory into the free energy functional. In order to elucidate the role of association, the phase diagrams are compared with those obtained by us recently [O. Pizio et al., J. Chem. Phys. 121, 11957 (2004)] for the restricted primitive model. The inclusion of the association into the theory leads to lowering the critical temperature for the fluid confined to pores with uncharged and with charged walls. We have observed that the average fraction of bonded ions is high along the coexistence envelope.

Phase behavior of ionic fluids in slitlike pores: A density functional approach for the restricted primitive model

We present a density functional theory of nonuniform ionic fluids. This theory is based on the application of the electrostatic contribution to the free energy functional arising from mean spherical approximation for a bulk restricted primitive model and from the energy route bulk equation of state. In order to employ this functional we define a reference fluid and additional averaged densities, according to the approach introduced by Gillespie, Nonner and Eisenberg [J. Phys.: Condens. Matter 14, 12129 (2002)]. In the case of bulk systems the proposed theory reduces to the mean spherical approximation equation of state, arising from the energy route and thus it predicts the first-order phase transition. We use this theory to investigate the effects of confinement on the liquid-vapor equilibria. Two cases are considered, namely an electrolyte confined to the pore with uncharged walls and with charged walls. The dependence of the capillary evaporation diagrams on the pore width and on the electrostatic potential is determined.

Phase behavior and critical properties of size-asymmetric, primitive-model electrolytes

The theory of J. Jiang et al. [J. Chem. Phys. 116, 7977 (2002)] for size-symmetric electrolytes is extended to size-asymmetric electrolytes. When compared to molecular-simulation results, this extension gives the correct trend of critical properties with size asymmetry.

Phase diagram of charged dumbbells: a random phase approximation approach

The phase diagram of the charged hard dumbbell system (hard spheres of opposite unit charge fixed at contact) is obtained with the use of the random phase approximation (RPA). The effect of the impenetrability of charged spheres on charge-charge fluctuations is described by introduction of a modified electrostatic potential. The correlations of ions in a pair are included via a correlation function in the RPA. The coexistence curve is in good agreement with Monte Carlo simulations. The relevance of the theory to the restricted primitive model is discussed.

Phase coexistence in polydisperse charged hard-sphere fluids: Mean spherical approximation

Taking advantage of the availability of the analytic solution of the mean spherical approximation for a mixture of charged hard spheres with an arbitrary number of components we show that the polydisperse fluid mixture of charged hard spheres belongs to the class of truncatable free energy models, i.e., to those systems where the thermodynamic properties can be represented by a finite number of (generalized) moments of the distribution function that characterizes the mixture. Thus, the formally infinitely many equations that determine the parameters of the two coexisting phases can be mapped onto a system of coupled nonlinear equations in these moments. We present the formalism and demonstrate the power of this approach for two systems; we calculate the full phase diagram in terms of cloud and shadow curves as well as binodals and discuss the distribution functions of the coexisting daughter phases and their charge distributions.

Precipitation of oppositely charged polyelectrolytes in salt solutions

We study phase separation in symmetric solutions of weakly charged flexible chains of opposite sign. Precipitation is caused by effective attractions due to charge fluctuations and by short-range attractions between monomers. The contribution from charge fluctuations is computed within the random phase approximation (RPA), which takes into account the connectivity of charges in the polyions. The impenetrability of the ions is accounted for by using a modified Coulomb potential in the RPA. In good solvent conditions the precipitate monotonically swells and eventually dissolves upon addition of salt. However, near the theta-solvent condition, but still in the good solvent, the precipitate can be stable at any salt concentration. Moreover, the density of the precipitate after initial decrease can increase with addition of salt. This effect is a result of redistribution of salt between the precipitate and the supernatant, which is due to an interplay of electrostatic and hardcore interactions. For not too weakly charged polyions the precipitate properties become strongly dependent on temperature even in good solvent conditions.

Thermodynamics of electrolytes on anisotropic lattices

The phase behavior of ionic fluids on simple cubic and tetragonal (anisotropic) lattices has been studied by grand canonical Monte Carlo simulations. Systems with both the true lattice Coulombic potential and continuous-space 1/r electrostatic interactions have been investigated. At all degrees of anisotropy, only coexistence between a disordered low-density phase and an ordered high-density phase with the structure similar to ionic crystal was found, in contrast to recent theoretical predictions. Tricritical parameters were determined to be increasing functions of anisotropy parameters which is consistent with theoretical calculations based on the Debye-Hückel approach. At large anisotropies a two-dimensional-like behavior is observed, from which we estimated the dimensionless tricritical temperature and density for the two-dimensional square lattice electrolyte to be T*(tri)=0.14 and rho*(tri)=0.70.

Dipolar origin of the gas-liquid coexistence of the hard-core 1:1 electrolyte model

We present a systematic study of the effect of the ion pairing on the gas-liquid phase transition of hard-core 1:1 electrolyte models. We study a class of dipolar dimer models that depend on a parameter R(c), the maximum separation between the ions that compose the dimer. This parameter can vary from sigma(+/-) that corresponds to the tightly tethered dipolar dimer model to R(c)--> infinity that corresponds to the Stillinger-Lovett description of the free ion system. The coexistence curve and critical point parameters are obtained as a function of R(c) by grand-canonical Monte Carlo techniques. Our results show that this dependence is smooth but nonmonotonic and converges asymptotically towards the free ion case for relatively small values of R(c). This fact allows us to describe the gas-liquid transition in the free ion model as a transition between two dimerized fluid phases. The role of the unpaired ions can be considered as a perturbation of this picture.