Free energy and phase equilibria for the restricted primitive model of ionic fluids from Monte Carlo simulations

Theory and Monte Carlo simulation of the ideal gas with shell particles in the canonical, isothermal-isobaric, grand canonical, and Gibbs ensembles

Theories of small systems play an important role in the fundamental understanding of finite size effects in statistical mechanics, as well as the validation of molecular simulation results as no computer can simulate fluids in the thermodynamic limit. Previously, a shell particle was included in the isothermal-isobaric ensemble in order to resolve an ambiguity in the resulting partition function. The shell particle removed either redundant volume states or redundant translational degrees of freedom of the system and yielded quantitative differences from traditional simulations in this ensemble. In this work, we investigate the effect of including a shell particle in the canonical, grand canonical, and Gibbs ensembles. For systems comprised of a pure component ideal gas, analytical expressions for various thermodynamic properties are obtained. We also derive the Metropolis Monte Carlo simulation acceptance criteria for these ensembles with shell particles, and the results of the simulations of an ideal gas are in excellent agreement with the theoretical predictions. The system size dependence of various important ensemble averages is also analyzed.

Anomalous Underscreening in the Restricted Primitive Model

Underscreening is a collective term for charge correlations in electrolytes decaying slower than the Debye length. Anomalous underscreening refers to phenomenology that cannot be attributed alone to steric interactions. Experiments with concentrated electrolytes and ionic fluids report anomalous underscreening, which so far has not been observed in simulation. We present Molecular Dynamics simulation results exhibiting anomalous underscreening that can be connected to cluster formation. A theory that accounts for ion pairing confirms the trend. Our results challenge the classic understanding of dense electrolytes impacting the design of technologies for energy storage and conversion.

Self-Consistent Description of Vapor-Liquid Interface in Ionic Fluids

Inhomogeneity of ion correlation widely exists in many physicochemical, soft matter, and biological systems. Here, we apply the modified Gaussian renormalized fluctuation theory to study the classic example of the vapor-liquid interface of ionic fluids. The ion correlation is decomposed into a short-range contribution associated with the local electrostatic environment and a long-range contribution accounting for the spatially varying ionic strength and dielectric permittivity. For symmetric salt, both the coexistence curve and the interfacial tension predicted by our theory are in quantitative agreement with simulation data reported in the literature. Furthermore, we provide the first theoretical prediction of interfacial structure for asymmetric salt, highlighting the importance of capturing local charge separation.

Coexistence of Transient Liquid Droplets and Amorphous Solid Particles in Nonclassical Crystallization of Cerium Oxalate

Crystallization from solution often occurs via "nonclassical" routes; that is, it involves transient, non-crystalline states like reactant-rich liquid droplets and amorphous particles. However, in mineral crystals, the well-defined thermodynamic character of liquid droplets and whether they convert─or not─into amorphous phases have remained unassessed. Here, by combining cryo-transmission electron microscopy and X-ray scattering down to a 250 ms reaction time, we unveil that crystallization of cerium oxalate involves a metastable chemical equilibrium between transient liquid droplets and solid amorphous particles: contrary to the usual expectation, reactant-rich droplets do not evolve into amorphous solids. Instead, at concentrations above 2.5 to 10 mmol L^-1, both amorphous and reactant-rich liquid phases coexist for several tens of seconds and their molar fractions remain constant and follow the lever rule in a multicomponent phase diagram. Such a metastable chemical equilibrium between solid and liquid precursors has been so far overlooked in multistep nucleation theories and highlights the interest of rationalizing phase transformations using multicomponent phase diagrams not only when designing and recycling rare earths materials but also more generally when describing nonclassical crystallization.

The primitive model in classical density functional theory: beyond the standard mean-field approximation

The primitive model describes ions by point charges with an additional hard-core interaction. In classical density-functional theory (DFT) the mean-field electrostatic contribution can be obtained from the first order of a functional perturbation of the pair potential for an uncharged reference system of hard spheres. This mean-field electrostatic term particularly contributes at particle separations that are forbidden due to hard-core overlap. In this work we modify the mean-field contribution such that the pair potential is constant for distances smaller than the contact distance of the ions. We motivate our modification by the underlying splitting of the potential, which is similar to the splitting of the Weeks-Chandler-Andersen potential and leads to higher-order terms in the respective expansion of the functional around the reference system. The resulting formalism involves weighted densities similar to the ones found in fundamental measure theory. To test our modifications, we analyze and compare density profiles, direct and total correlation functions, and the thermodynamic consistency of the functional via a widely established sum rule and the virial pressure formula for our modified functional, for established functionals, and for data from computer simulations. We found that our modifications clearly show improvements compared to the standard mean-field functional, especially when predicting layering effects and direct correlation functions in high concentration scenarios; for the latter we also find improved consistency when calculated via different thermodynamic routes. In conclusion, we demonstrate how modifications toward higher order corrections beyond mean-field functionals can be made and how they perform, by this providing a basis for systematic future improvements in classical DFT for the description of electrostatic interactions.

Thermoresponsive Ionic Liquid/Water Mixtures: From Nanostructuring to Phase Separation

The thermodynamics, structures, and applications of thermoresponsive systems, consisting primarily of water solutions of organic salts, are reviewed. The focus is on organic salts of low melting temperatures, belonging to the ionic liquid (IL) family. The thermo-responsiveness is represented by a temperature driven transition between a homogeneous liquid state and a biphasic state, comprising an IL-rich phase and a solvent-rich phase, divided by a relatively sharp interface. Demixing occurs either with decreasing temperatures, developing from an upper critical solution temperature (UCST), or, less often, with increasing temperatures, arising from a lower critical solution temperature (LCST). In the former case, the enthalpy and entropy of mixing are both positive, and enthalpy prevails at low T. In the latter case, the enthalpy and entropy of mixing are both negative, and entropy drives the demixing with increasing T. Experiments and computer simulations highlight the contiguity of these phase separations with the nanoscale inhomogeneity (nanostructuring), displayed by several ILs and IL solutions. Current applications in extraction, separation, and catalysis are briefly reviewed. Moreover, future applications in forward osmosis desalination, low-enthalpy thermal storage, and water harvesting from the atmosphere are discussed in more detail.

The differential capacitance as a probe for the electric double layer structure and the electrolyte bulk composition

In this work, we theoretically study the differential capacitance of an aqueous electrolyte in contact with a planar electrode, using classical density functional theory, and show how this measurable quantity can be used as a probe to better understand the structure and composition of the electric double layer at play. Specifically, we show how small trace amounts of divalent ions can influence the differential capacitance greatly and also how small ions dominate its behavior for high electrode potentials. In this study, we consider primitive model electrolytes and not only use the standard definition of the differential capacitance but also derive a new expression from mechanical equilibrium in a planar geometry. This expression reveals explicitly that the first layer of ions near the charged surface is key to its understanding. Our insights might be used as a guide in experiments to better understand the electrolyte-electrode interface as well as the (composition of the) bulk electrolyte.

Primitive model electrolytes in the near and far field: Decay lengths from DFT and simulations

Inspired by recent experimental observations of anomalously large decay lengths in concentrated electrolytes, we revisit the Restricted Primitive Model (RPM) for an aqueous electrolyte. We investigate the asymptotic decay lengths of the one-body ionic density profiles for the RPM in contact with a planar electrode using classical Density Functional Theory (DFT) and compare these with the decay lengths of the corresponding two-body correlation functions in bulk systems, obtained in previous Integral Equation Theory (IET) studies. Extensive Molecular Dynamics (MD) simulations are employed to complement the DFT and IET predictions. Our DFT calculations incorporate electrostatic interactions between the ions using three different (existing) approaches: one is based on the simplest mean-field treatment of Coulomb interactions (MFC), while the other two employ the Mean Spherical Approximation (MSA). The MSAc invokes only the MSA bulk direct correlation function, whereas the MSAu also incorporates the MSA bulk internal energy. Although MSAu yields profiles that are in excellent agreement with MD simulations in the near field, in the far field, we observe that the decay lengths are consistent between IET, MSAc, and MD simulations, whereas those from MFC and MSAu deviate significantly. Using DFT, we calculated the solvation force, which relates directly to surface force experiments. We find that its decay length is neither qualitatively nor quantitatively close to the large decay lengths measured in experiments and conclude that the latter cannot be accounted for by the primitive model. The anomalously large decay lengths found in surface force measurements require an explanation that lies beyond primitive models.

Improving the efficiency of Monte Carlo simulations of ions using expanded grand canonical ensembles

While ionic liquids have promising applications as industrial solvents, predicting their fluid phase properties and coexistence remains a challenge. Grand canonical Monte Carlo simulation is an effective method for such predictions, but equilibration is hampered by the apparent requirement to insert and delete neutral sets of ions simultaneously in order to maintain charge neutrality. For relatively high densities and low temperatures, previously developed methods have been shown to be essential in improving equilibration by gradual insertion and deletion of these neutral sets of ions. We introduce an expanded ensemble approach which may be used in conjunction with these existing methods to further improve efficiency. Individual ions are inserted or deleted in one Monte Carlo trial rather than simultaneous insertion/deletion of neutral sets. We show how charge neutrality is maintained and show rigorous quantitative agreement between the conventional and the proposed expanded ensemble approaches, but with up to an order of magnitude increase in efficiency at high densities. The expanded ensemble approach is also more straightforward to implement than simultaneous insertion/deletion of neutral sets, and its implementation is demonstrated within open source software.

Designing Electrostatic Interactions via Polyelectrolyte Monomer Sequence

Charged polymers are ubiquitous in biological systems because electrostatic interactions can drive complicated structure formation and respond to environmental parameters such as ionic strength and pH. In these systems, function emerges from sophisticated molecular design; for example, intrinsically disordered proteins leverage specific sequences of monomeric charges to control the formation and function of intracellular compartments known as membraneless organelles. The role of a charged monomer sequence in dictating the strength of electrostatic interactions remains poorly understood despite extensive evidence that sequence is a powerful tool biology uses to tune soft materials. In this article, we use a combination of theory, experiment, and simulation to establish the physical principles governing sequence-driven control of electrostatic interactions. We predict how arbitrary sequences of charge give rise to drastic changes in electrostatic interactions and correspondingly phase behavior. We generalize a transfer matrix formalism that describes a phase separation phenomenon known as "complex coacervation" and provide a theoretical framework to predict the phase behavior of charge sequences. This work thus provides insights into both how charge sequence is used in biology and how it could be used to engineer properties of synthetic polymer systems.

Understanding the Reactive Adsorption of H2 S and CO2 in Sodium-Exchanged Zeolites

Purifying sour natural gas streams containing hydrogen sulfide and carbon dioxide has been a long-standing environmental and economic challenge. In the presence of cation-exchanged zeolites, these two acid gases can react to form carbonyl sulfide and water (H₂ S+CO₂ ⇌H₂ O+COS), but this reaction is rarely accounted for. In this work, we carry out reactive first-principles Monte Carlo (RxFPMC) simulations for mixtures of H₂ S and CO₂ in all-silica and Na-exchanged forms of zeolite beta to understand the governing principles driving the enhanced conversion. The RxFPMC simulations show that the presence of Na⁺ cations can change the equilibrium constant by several orders of magnitude compared to the gas phase or in all-silica beta. The shift in the reaction equilibrium is caused by very strong interactions of H₂ O with Na⁺ that reduce the reaction enthalpy by about 20 kJ mol^-1 . The simulations also demonstrate that the siting of Al atoms in the framework plays an important role. The RxFPMC method presented here is applicable to any chemical conversion in any confined environment, where strong interactions of guest molecules with the host framework and high activation energies limit the use of other computational approaches to study reaction equilibria.

Structure of electric double layers in capacitive systems and to what extent (classical) density functional theory describes it

Ongoing scientific interest is aimed at the properties and structure of electric double layers (EDLs), which are crucial for capacitive energy storage, water treatment, and energy harvesting technologies like supercapacitors, desalination devices, blue engines, and thermocapacitive heat-to-current converters. A promising tool to describe their physics on a microscopic level is (classical) density functional theory (DFT), which can be applied in order to analyze pair correlations and charge ordering in the primitive model of charged hard spheres. This simple model captures the main properties of ionic liquids and solutions and it predicts many of the phenomena that occur in EDLs. The latter often lead to anomalous response in the differential capacitance of EDLs. This work constructively reviews the powerful theoretical framework of DFT and its recent developments regarding the description of EDLs. It explains to what extent current approaches in DFT describe structural ordering and in-plane transitions in EDLs, which occur when the corresponding electrodes are charged. Further, the review briefly summarizes the history of modeling EDLs, presents applications, and points out limitations and strengths in present theoretical approaches. It concludes that DFT as a sophisticated microscopic theory for ionic systems is expecting a challenging but promising future in both fundamental research and applications in supercapacitive technologies.

Effective short-range Coulomb correction to model the aggregation behavior of ionic surfactants

We present a short-range correction to the Coulomb potential to investigate the aggregation of amphiphilic molecules in aqueous solutions. The proposed modification allows to quantitatively reproduce the distribution of counterions above the critical micelle concentration (CMC) or, equivalently, the degree of ionization, α, of the micellar clusters. In particular, our theoretical framework has been applied to unveil the behavior of the cationic surfactant C24H49N2O2 (+) CH3SO4 (-), which offers a wide range of applications in the thriving and growing personal care market. A reliable and unambiguous estimation of α is essential to correctly understand many crucial features of the micellar solutions, such as their viscoelastic behavior and transport properties, in order to provide sound formulations for the above mentioned personal care solutions. We have validated our theory by performing extensive lattice Monte Carlo simulations, which show an excellent agreement with experimental observations. More specifically, our coarse-grained model is able to reproduce and predict the complex morphology of the micelles observed at equilibrium. Additionally, our simulation results disclose the existence of a transition from a monodisperse to a bidisperse size distribution of aggregates, unveiling the intriguing existence of a second CMC.

Modelling aqueous solubility of sodium chloride in clays at thermodynamic conditions of hydraulic fracturing by molecular simulations

To address the high salinity of flow-back water during hydraulic fracturing, we have studied the equilibrium partitioning of NaCl and water between the bulk phase and clay pores.

Charge ordering and scattering pre-peaks in ionic liquids and alcohols

The structural properties of ionic liquids and alcohols are viewed under the charge ordering process as a common basis to explain the peculiarity of their radiation scattering properties, namely the presence, or absence, of a scattering pre-peak.

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.

Method for sampling compact configurations for semistiff polymers

The sampling of compact configurations is crucial when investigating structural properties of semistiff polymers, like proteins and DNA, using Monte Carlo methods. A sampling scheme for a continuous model based on configuration biasing is introduced, tested, and compared with conventional methods. The proposed configuration biased Monte Carlo method, used together with the Wang-Landau sampling scheme, enables us to obtain any thermodynamic property within the statistical ensemble in use. Using the proposed method, it is possible to collect statistical data of interest for a wide range of compactions (from stretched up to several toroid loops) in a single computer experiment. A second-order-like stretched-toroid phase transition is observed for a semistiff polymer, and the critical temperature is estimated.

Computer simulations of the restricted primitive model at very low temperature and density

The problem of successfully simulating ionic fluids at low temperature and low density states is well known in the simulation literature: using conventional methods, the system is not able to equilibrate rapidly due to the presence of strongly associated cation-anion pairs. In this paper we present a numerical method for speeding up computer simulations of the restricted primitive model (RPM) at low temperatures (around the critical temperature) and at very low densities (down to 10(-10)σ(-3), where σ is the ion diameter). Experimentally, this regime corresponds to typical concentrations of electrolytes in nonaqueous solvents. As far as we are aware, this is the first time that the RPM has been equilibrated at such extremely low concentrations. More generally, this method could be used to equilibrate other systems that form aggregates at low concentrations.

A computational study of electrolyte adsorption in a simple model for intercalated clays

A pillared interlayered clay is represented by a two-dimensional quenched charged disordered medium, in which the pillar configuration is produced by the quench of a two-dimensional electrolyte and the subsequent removal of the anions (that act as a template). The cation charge is counterbalanced by a neutralizing background that is an ideal representation of the layer's negative charge in the experimental system. In this paper we investigate the adsorption of electrolyte particles in this charged disordered medium resorting both to the use of the replica Ornstein-Zernike equation in the hypernetted chain approximation and grand canonical Monte Carlo simulations. The theoretical approach qualitatively reproduces the simulated behavior of the adsorbed fluids. Theoretical estimates of the material porosities obtained for various types of pillar distributions are in good agreement with the simulation. We investigate the influence of the matrix on correlation functions and adsorption isotherms.

Modeling the anisotropic self-assembly of spherical polymer-grafted nanoparticles

Recent experimental results demonstrated that polymer grafted nanoparticles in solvents display self-assembly behavior similar to the microphase separation of block copolymers and other amphiphiles. We present a mean-field theory and complementary computer simulations to shed light on the parametric underpinnings of the experimental observations. Our theory suggests that such self-assembled structures occur most readily when the nanoparticle size is comparable to the radius of gyration of the polymer brush chains. Much smaller particle sizes are predicted to yield uniform particle dispersions, while larger particles are expected to agglomerate due to phase separation from the solvent. Selected aspects of our theoretical predictions are corroborated by computer simulations.

Charge correlation effects on ionization of weak polyelectrolytes

Ionization curves of weak polyelectrolytes were obtained as a function of the charge coupling strength from Monte Carlo simulations. In contrast to many earlier studies, the present work treats counterions explicitly, thus allowing the investigation of charge correlation effects at strong couplings. For conditions representing typical weak polyelectrolytes in water near room temperature, ionization is suppressed because of interactions between nearby dissociated groups, as also seen in prior work. A novel finding here is that, for stronger couplings, relevant for non-aqueous environments in the absence of added salt, the opposite behavior is observed-ionization is enhanced relative to the behavior of the isolated groups due to ion-counterion correlation effects. The fraction of dissociated groups as a function of position along the chain also behaves non-monotonically. Dissociation is highest near the ends of the chains for aqueous polyelectrolytes and highest at the chain middle segments for non-aqueous environments. At intermediate coupling strengths, dissociable groups appear to behave in a nearly ideal fashion, even though chain dimensions still show strong expansion effects due to ionization. These findings provide physical insights on the impact of competition between acid/base chemical equilibrium and electrostatic attractions in ionizable systems.

Stability mechanisms for plate-like nanoparticles immersed in a macroion dispersion

An integral equation theory and Monte Carlo simulations are applied to study a model macroion solution confined between two parallel plates immersed in a 1:1 electrolyte and the macroions' counterions. We analyze the cases in which plates are: (a) uncharged; (b) when they are like-charged to the macroions; (c) when they are oppositely charged to the macroions. For all cases a long range oscillatory behavior of the induced charge density between the plates is found (implying an overcompensation/undercompensation of the plates' charge density) and a correlation between the confined and outside fluids. The behavior of the force is discussed in terms of the macroion and ion structure inside and outside the plates. A good agreement is found between theoretical and simulation results.

Liquid-liquid phase separation in solutions of ionic liquids: phase diagrams, corresponding state analysis and comparison with simulations of the primitive model

Phase diagrams of ionic solutions of the ionic liquid C(18)mim(+)NTF(2)(-) (1-n-octadecyl-3-methyl imidazolium bistrifluormethylsulfonylimide) in decalin, cyclohexane and methylcyclohexane are reported and compared with that of solutions of other imidazolium ionic liquids with the anions NTF(2)(-), Cl(-) and BF4(-) in arenes, CCl(4), alcohols and water. The phase diagrams are analysed presuming Ising criticality and taking into account the asymmetry of the phase diagrams. The resulting parameters are compared with simulation results for equal-sized charged hard spheres in a dielectric continuum, the restricted primitive model (RPM) and the primitive model (PM) that allows for ions of different size. In the RPM temperature scale the critical temperatures vary almost linearly with the dielectric permittivity of the solvent. The RPM critical temperatures of the solutions in non-polar solvents are very similar, somewhat below the RPM value. Correlations with the boiling temperatures of the solvents and a dependence on the length of the side chain of the imidazolium cations show that dispersion interactions modify the phase transition, which is mainly determined by Coulomb forces. Critical concentrations, widths of the phase diagrams and the slopes of the diameter are different for the solutions in protic and aprotic solvents. The phase diagrams of the solutions in alcohols and water get a lower critical solution point when represented in RPM variables.

A cluster algorithm for Monte Carlo simulation at constant pressure

We propose an efficient algorithm to sample the volume in Monte Carlo simulations in the isobaric-isothermal ensemble. The method is designed to be applied in the simulation of hard-core models at high density. The algorithm is based in the generation of clusters of particles. At the volume change step, the distances between pairs of particles belonging to the same cluster do not change. This is done by rescaling the positions of the center of mass of each cluster instead of the position of each individual particle. We have tested the performance of the algorithm by simulating fluid and solid phases of hard spheres, finding that in both cases the algorithm is much more efficient than the standard procedure. Moreover, the efficiency of the method measured in terms of correlation "time" does not depend on the system size in contrast with the standard method, in which the sampling becomes rapidly inefficient as the system size increases. We have used the procedure to compute with high precision the equation of state of the face-centered-cubic phase of the hard sphere system for different system sizes. Using these results we have estimated the equation of state at the thermodynamic limit. The results are compared to different equations of state proposed in literature.