Field theory for size- and charge-asymmetric primitive model of ionic systems: mean-field stability analysis and pretransitional effects

Mesoscopic Inhomogeneities in Concentrated Electrolytes

A mesoscopic theory for water-in-salt electrolytes combining density functional and field-theoretic methods is developed in order to explain the unexpectedly large period of the oscillatory decay of the disjoining pressure observed in recent experiments for the lithium bis(trifluoromethylsulfonyl)-imide (LiTFSI) salt [T. S. Groves et al., J. Phys. Chem. Lett. 2021, 12, 1702]. We assumed spherical ions with different diameters and implicit solvent, inducing strong, short-range attraction between ions of the same sign. For this highly simplified model, we calculated correlation functions. Our results indicate that mesoscopic inhomogeneities can occur when the sum of the Coulomb and the water-mediated interactions between like ions is attractive at short and repulsive at large distances. We adjusted the attractive part of the potential to the water-in-LiTFSI electrolyte and obtained both the period and the decay rate of the correlations, in semiquantitative agreement with the experiment. In particular, the decay length of the correlations increases nearly linearly with the volume fraction of ions.

Critical charge and density coupling in ionic spherical models

We investigate ionic criticality on the basis of a specially devised spherical model that accounts both for Coulomb and nonionic forces in binary systems. We show in detail here the consequences of the entanglement of density and charge correlation functions G_{NN} and G_{ZZ} on criticality and screening. We also show on this soluble model how, because of electroneutrality, the long-range Coulomb interactions do not change the universality class of criticality in the model driven primarily by sufficiently attractive nonionic interactions. Near criticality, G_{NN} and G_{ZZ} are fully decoupled in charge symmetric systems. However, in more realistic nonsymmetric models, charge and density fluctuations couple in leading order so that the charge and density correlation lengths diverge asymptotically in a similar way. Similarly, the Stillinger-Lovett sum rule, which characterizes a conducting fluid, is violated at criticality in nonsymmetric models when the critical-point density-decay exponent η vanishes. In addition, if quantum effects are accounted for semiclassically by incorporating algebraically decaying interactions, G_{ZZ} decays only as a power law in the whole phase space, contrary to the usually expected exponential Debye screening. We expect these results on this soluble toy model to be general and to reveal general mechanisms ruling ionic criticality.

Self-consistent theory for systems with mesoscopic fluctuations

We have developed a theory for inhomogeneous systems that allows for the incorporation of the effects of mesoscopic fluctuations. A hierarchy of equations relating the correlation and direct correlation functions for the local excess [Formula: see text] of the volume fraction of particles ζ has been obtained, and an approximation leading to a closed set of equations for the two-point functions has been introduced for the disordered inhomogeneous phase. We have numerically solved the self-consistent equations for one-dimensional (1D) and three-dimensional (3D) models with short-range attraction and long-range repulsion. Predictions for all of the qualitative properties of the 1D model agree with the exact results, but only semi-quantitative agreement is obtained in the simplest version of the theory. The effects of fluctuations in the two 3D models considered are significantly different, despite the very similar properties of these models in the mean-field approximation. In both cases we obtain the sequence of large-small-large compressibility for increasing ζ. The very small compressibility is accompanied by the oscillatory decay of correlations with correlation lengths that are orders of magnitude larger than the size of particles. In one of the two models considered, the small compressibility becomes very small and the large compressibility becomes very large with decreasing temperature, and eventually van der Waals loops appear. Further studies are necessary in order to determine the nature of the strongly inhomogeneous phase present for intermediate volume fractions in 3D.

Spatial inhomogeneities in ionic liquids, charged proteins, and charge stabilized colloids from collective variables theory

Effects of size and charge asymmetry between oppositely charged ions or particles on spatial inhomogeneities are studied for a large range of charge and size ratios. We perform a stability analysis of the primitive model of ionic systems with respect to periodic ordering using the collective variables-based theory. We extend previous studies [Ciach et al., Phys. Rev. E 75, 051505 (2007)] in several ways. First, we employ a nonlocal approximation for the reference hard-sphere fluid which leads to the Percus-Yevick pair direct correlation functions for the uniform case. Second, we use the Weeks-Chandler-Anderson regularization scheme for the Coulomb potential inside the hard core. We determine the relevant order parameter connected with the periodic ordering and analyze the character of the dominant fluctuations along the λ lines. We show that the above-mentioned modifications produce large quantitative and partly qualitative changes in the phase diagrams obtained previously. We discuss possible scenarios of the periodic ordering for the whole range of size and charge ratios of the two ionic species, covering electrolytes, ionic liquids, charged globular proteins or nanoparticles in aqueous solutions, and charge-stabilized colloids.

Thermodynamics and diffusion in size-symmetric and asymmetric dense electrolytes

MD simulation results for model size-symmetric and asymmetric electrolytes at high densities and temperatures (well outside the liquid-gas coexistence region) are generated and analyzed focusing on thermodynamic and diffusion properties. An extension of the mean spherical approximation for electrolytes originally derived for charged hard sphere fluids is adapted to these systems by exploiting the separation of short range and Coulomb interaction contributions intrinsic to these theoretical models and is found to perform well for predicting equation of state quantities. The diffusion coefficients of these electrolytes can also be reasonably well predicted using entropy scaling ideas suitably adapted to charged systems and mixtures. Thus, this approach may provide an avenue for studying dense electrolytes or complex molecular systems containing charged species at high pressures and temperatures.

Gas-liquid critical parameters of asymmetric models of ionic fluids

The effects of size and charge asymmetry on the gas-liquid critical parameters of a primitive model (PM) of ionic fluids are studied within the framework of the statistical field theory based on the collective variables method. Recently, this approach has enabled us to obtain the correct trends of the both critical parameters of the equisize charge-asymmetric PM without assuming ionic association. In this paper, we focus on the general case of an asymmetric PM characterized by the two parameters: hard-sphere diameter, lambda=sigma+/sigma-, and charge, z=q+/|q-|, ratios of the two ionic species. We derive an explicit expression for the chemical potential conjugate to the order parameter which includes the effects of correlations up to the third order. Based on this expression we consider the three versions of PM: a monovalent size-asymmetric PM (lambda not equal 1, z=1) , an equisize charge-asymmetric PM (lambda=1, z not equal 1) and a size- and charge-asymmetric PM (lambda not equal 1, z=2) . Similar to simulations, our theory predicts that the critical temperature and the critical density decrease with the increase in size asymmetry. Regarding the effects of charge asymmetry, we obtain the correct trend of the critical temperature with z , while the trend of the critical density obtained in this approximation is inconsistent with simulations, as well as with our previous results found in the higher-order approximation. We expect that the consideration of the higher-order correlations will lead to the correct trend of the critical density with charge asymmetry.

Correlations in the Structure and Dynamics of Ionic Liquids

A very great deal of the experimental work on room-temperature ionic liquids was done after high-performance computing became readily available for quantum-chemical or molecular-dynamic calculations. I explore the use of modern computational methods to guide or aid laboratory work, and the importance of ‘old-fashioned’ theory, from before the age of fast computers. Debye and Hückel published the first really important theoretical work on correlations between charged particles, and the Nernst–Einstein formula is still used to understand electrical conductivities. I assess the usefulness of all these theoretical methods and ideas, and discuss the particular difficulties presented by ionic liquids.

Universal sequence of ordered structures obtained from mesoscopic description of self-assembly

Mesoscopic theory for soft-matter systems that combines density functional and statistical field theory is derived by a systematic coarse-graining procedure. For particles interacting with spherically symmetric potentials of arbitrary form, the grand-thermodynamic potential consists of two terms. The first term is associated with microscopic length-scale fluctuations, and has the form of the standard density functional. The second term is associated with mesoscopic length-scale fluctuations, and has the form known from the statistical field theory. For the correlation function between density fluctuations in mesoscopic regions, a pair of equations similar to the Ornstein-Zernicke equation with a new closure is obtained. In the special case of weak ordering on the mesoscopic length scale, the theory takes a form similar to either the Landau-Ginzburg-Wilson (LGW) or the Landau-Brazovskii (LB) field theory, depending on the form of the interaction potential. Microscopic expressions for the phenomenological parameters that appear in the Landau-type theories are obtained. Within the framework of this theory, we obtain a universal sequence of phases: disordered, bcc, hexagonal, lamellar, inverted hexagonal, inverted bcc, disordered, for increasing density well below the close-packing density. The sequence of phases agrees with experimental observations and with simulations of many self-assembling systems. In addition to the above phases, more complex phases may appear depending on the interaction potentials. For a particular form of the short-range attraction long-range repulsion potential, we find the bicontinuous gyroid phase ( Ia3d symmetry) that may be related to a network-forming cluster of colloids in a mixture of colloids and nonadsorbing polymers.