Effect of a nearby charge-ordered phase on correlation functions in ionic systems

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.

Correct scaling of the correlation length from a theory for concentrated electrolytes

Self-consistent theory for concentrated electrolytes is developed. Oscillatory decay of the charge-charge correlation function with the decay length that shows perfect agreement with the experimentally discovered and so far unexplained scaling is obtained. For the density-density correlations, monotonic asymptotic decay with the decay length comparable with the decay length of the charge correlations is found. We show that the correlation lengths in concentrated electrolytes depend crucially on the local variance of the charge density.

Structural crossover in a model fluid exhibiting two length scales: Repercussions for quasicrystal formation

We investigate the liquid state structure of the two-dimensional model introduced by Barkan et al. [Phys. Rev. Lett. 113, 098304 (2014)10.1103/PhysRevLett.113.098304], which exhibits quasicrystalline and other unusual solid phases, focusing on the radial distribution function g(r) and its asymptotic decay r→∞. For this particular model system, we find that as the density is increased there is a structural crossover from damped oscillatory asymptotic decay with one wavelength to damped oscillatory asymptotic decay with another distinct wavelength. The ratio of these wavelengths is ≈1.932. Following the locus in the phase diagram of this structural crossover leads directly to the region where quasicrystals are found. We argue that identifying and following such a crossover line in the phase diagram towards higher densities where the solid phase(s) occur is a good strategy for finding quasicrystals in a wide variety of systems. We also show how the pole analysis of the asymptotic decay of equilibrium fluid correlations is intimately connected with the nonequilibrium growth or decay of small-amplitude density fluctuations in a bulk fluid.

Density functional theory for systems with mesoscopic inhomogeneities

We study the effects of fluctuations on the mesoscopic length scale on systems with mesoscopic inhomogeneities. Equations for the correlation function and for the average volume fraction are derived in the self-consistent Gaussian approximation. The equations are further simplified by postulating the expression for the structure factor consistent with scattering experiments for self-assembling systems. Predictions of the approximate theory are verified by a comparison with the exact results obtained earlier for the one-dimensional lattice model with first-neighbor attraction and third-neighbor repulsion. We find qualitative agreement for the correlation function, the equation of state and the dependence of the chemical potential μ on the volume fraction ζ. Our results confirm also that strong inhomogeneities in the disordered phase are found only in the case of strong repulsion. The inhomogeneities are reflected in an oscillatory decay of the correlation function with a very large correlation length, three inflection points in the [Formula: see text] curve and a compressibility that for increasing ζ takes very large, very small and again very large values.

Modeling the structure of liquids and crystals using one- and two-component modified phase-field crystal models

A modified phase-field crystal model in which the free energy may be minimized by an order parameter profile having isolated bumps is investigated. The phase diagram is calculated in one and two dimensions and we locate the regions where modulated and uniform phases are formed and also regions where localized states are formed. We investigate the effectiveness of the phase-field crystal model for describing fluids and crystals with defects. We further consider a two-component model and elucidate how the structure transforms from hexagonal crystalline ordering to square ordering as the concentration changes. Our conclusion contains a discussion of possible interpretations of the order parameter field.

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.

Two-dimensional fluid with competing interactions exhibiting microphase separation: theory for bulk and interfacial properties

Colloidal particles that are confined to an interface such as the air-water interface are an example of a two-dimensional fluid. Such dispersions have been observed to spontaneously form cluster and stripe morphologies in certain systems with isotropic pair potentials between the particles, due to the fact that the pair interaction between the colloids has competing attraction and repulsion over different length scales. Here we present a simple density functional theory for a model of such a two-dimensional fluid. The theory predicts a bulk phase diagram exhibiting cluster, stripe, and bubble modulated phases, in addition to homogeneous fluid phases. Comparing with simulation results for this model from the literature, we find that the theory is qualitatively reliable. The model allows for a detailed investigation of the structure of the fluid and we are able to obtain simple approximate expressions for the static structure factor and for the length scale characterizing the modulations in the microphase separated phases. We also investigate the behavior of the system under confinement between two parallel hard walls. We find that the confined fluid phase behavior can be rather complex.

Charge fluctuations and correlation lengths in finite electrolytes

Fluctuations of the charge Q_{Lambda} inside a subdomain Lambda embedded in an electrolyte contained in a finite cubical box of dimensions LxLxL with periodic boundary conditions are investigated. When Lambda is an LxL "slab" of width W , asymptotically exact expressions for the mean-square fluctuation Q_{Lambda};{2} are obtained in terms of the Lebowitz length xi_{L}(T,rho) and of the "true" or asymptotic screening/decay length xi_{Z,infinity}(T,rho) together, when the charge correlation decay is oscillatory, with the characteristic wavelength lambda_{Z}(T,rho) . In finite systems, the normalized charge fluctuations exhibit threefold scaling behavior in the ratios xi_{Z,infinity}W , Wlambda_{Z} , and WL . This enables one to estimate all the correlation lengths away from criticality quite precisely from finite-size grand canonical Monte Carlo simulations. The results for xi_{Z,infinity},lambda_{Z} , and xi_{L} are presented for the restricted primitive model or hard-sphere 1:1 electrolyte for densities rho less, similar1.3rho_{c} and T greater, similar4T_{c} . The fitted values compare favorably with the expectations of generalized Debye-Hückel theory [Lee and Fisher, Europhys. Lett. 39, 611 (1997)]; specifically, if xi_{D} proportional, variant(Trho);{12} is the Debye length, we find xi_{Z,infinity}<xi_{D}<xi_{L} , although, since ion-pairing is neglected, the charge oscillations set in only for densities approximately 1.9 times larger than predicted.

Phase behavior of a fluid with competing attractive and repulsive interactions

Fluids in which the interparticle potential has a hard core, is attractive at moderate separations, and repulsive at large separations are known to exhibit unusual phase behavior, including stable inhomogeneous phases. Here we report a joint simulation and theoretical study of such a fluid, focusing on the relationship between the liquid-vapor transition line and any new phases. The phase diagram is studied as a function of the amplitude of the attraction for a certain fixed amplitude of the long ranged repulsion. We find that the effect of the repulsion is to substitute the liquid-vapor critical point and a portion of the associated liquid-vapor transition line, by two first-order transitions. One of these transitions separates the vapor from a fluid of spherical liquidlike clusters; the other separates the liquid from a fluid of spherical voids. At low temperature, the two transition lines intersect one another and a vapor-liquid transition line at a triple point. While most integral equation theories are unable to describe the new phase transitions, the Percus-Yevick approximation does succeed in capturing the vapor-cluster transition, as well as aspects of the structure of the cluster fluid, in reasonable agreement with the simulation results.

Why do ultrasoft repulsive particles cluster and crystallize? Analytical results from density-functional theory

We demonstrate the accuracy of the hypernetted chain closure and of the mean-field approximation for the calculation of the fluid-state properties of systems interacting by means of bounded and positive pair potentials with oscillating Fourier transforms. Subsequently, we prove the validity of a bilinear, random-phase density functional for arbitrary inhomogeneous phases of the same systems. On the basis of this functional, we calculate analytically the freezing parameters of the latter. We demonstrate explicitly that the stable crystals feature a lattice constant that is independent of density and whose value is dictated by the position of the negative minimum of the Fourier transform of the pair potential. This property is equivalent with the existence of clusters, whose population scales proportionally to the density. We establish that regardless of the form of the interaction potential and of the location on the freezing line, all cluster crystals have a universal Lindemann ratio Lf=0.189 at freezing. We further make an explicit link between the aforementioned density functional and the harmonic theory of crystals. This allows us to establish an equivalence between the emergence of clusters and the existence of negative Fourier components of the interaction potential. Finally, we make a connection between the class of models at hand and the system of infinite-dimensional hard spheres, when the limits of interaction steepness and space dimension are both taken to infinity in a particularly described fashion.

Nonasymptotic Critical Behavior of a Ternary Ionic System

Refractive indices n and salt concentrations ms of coexisting phases of the ternary system 1,4-dioxane + water + potassium chloride were measured along the liquid-liquid-solid coexistence curve near the liquid-liquid critical end point. Refractive index measurements were carried out in the range 0.689 x 10-3 < t = (T - Tc)/Tc < 0.118 while salt concentrations were determined for the temperature range 1.84 x 10-3 < t < 8.07 x 10-2. From these experimental results, compositions fD (mass fraction of dioxane on a salt-free basis) and densities rho of coexisting phases were obtained. The shape of the coexistence curve was analyzed using alternatively n, ms, fD, and rho as order parameters. In all cases, the obtained coexistence curve displays, asymptotically, Ising behavior. Outside the asymptotic critical domain, n, ms, and rho show significant deviations of the effective critical exponent from its Ising value, while the concentration variable fD requires no corrections to simple scaling. On the basis of the present results, we conclude that this system shows no indication of multicritical behavior.

Model colloidal fluid with competing interactions: Bulk and interfacial properties

Using a simple mean field density functional theory (DFT), the authors investigate the structure and phase behavior of a model colloidal fluid composed of particles interacting via a pair potential which has a hard core of diameter sigma, is attractive Yukawa at intermediate separations, and is repulsive Yukawa at large separations. The authors analyze the form of the asymptotic decay of the bulk fluid correlation functions, comparing results from DFT with those from the self-consistent Ornstein-Zernike approximation (SCOZA). In both theories the authors find rich crossover behavior, whereby the ultimate decay of correlation functions changes from monotonic to long wavelength damped oscillatory decay on crossing certain lines in the phase diagram or sometimes from oscillatory to oscillatory with a longer wavelength. For some choices of potential parameters the authors find, within the DFT, a lambda line at which the fluid becomes unstable with respect to periodic density fluctuations. SCOZA fails to yield solutions for state points near such a lambda line. The propensity towards clustering of particles, which is reflected by the presence of a long wavelength (>>sigma) slowly decaying oscillatory pair correlation function, and a structure factor that exhibits a very sharp maximum at small but nonzero wave numbers, is enhanced in states near the lambda line. The authors present density profiles for the planar liquid-gas interface and for fluids adsorbed at a planar hard wall. The presence of a nearby lambda transition gives rise to pronounced long wavelength oscillations in the one-body density profiles at both types of interface.

Field-theoretic description of ionic crystallization in the restricted primitive model

Effects of charge-density fluctuations on a phase behavior of the restricted primitive model are studied within a field-theoretic formalism. We focus on a lambda line of continuous transitions between charge-ordered and charge-disordered phases that is observed in several mean-field theories, but is absent in simulation results. In our study the RPM is reduced to a phi(6) theory, and a fluctuation contribution to a grand thermodynamic potential is obtained by generalizing the Brazovskii approach. We find that in a presence of fluctuations the lambda -line disappears. Instead, a fluctuation-induced first-order transition to a charge-ordered phase appears in the same region of a phase diagram, where the liquid-ionic-crystal transition is obtained in simulations. Our results indicate that the charge-ordered phase should be identified with an ionic crystal.

A molecular dynamics computer simulation study of room-temperature ionic liquids. I. Equilibrium solvation structure and free energetics

Solvation in 1-ethyl-3-methylmidazolium chloride and in 1-ethyl-3-methylimidazolium hexafluorophosphate near equilibrium is investigated via molecular dynamics computer simulations with diatomic and benzenelike molecules employed as probe solutes. It is found that electrostriction plays an important role in both solvation structure and free energetics. The angular and radial distributions of cations and anions become more structured and their densities near the solute become enhanced as the solute charge separation grows. Due to the enhancement in structural rigidity induced by electrostriction, the force constant associated with solvent configuration fluctuations relevant to charge shift and transfer processes is also found to increase. The effective polarity and reorganization free energies of these ionic liquids are analyzed and compared with those of highly polar acetonitrile. Their screening behavior of electric charges is also investigated.