Fluid-solid equilibrium of a charged hard-sphere model

Monte Carlo study of molten salt with charge asymmetry near the electrode surface

Results of the Monte Carlo simulation of the electrode | molten salt or ionic liquid interface are reported. The system investigated is approximated by the primitive model of electrolyte being in contact with a charged hard wall. Ions differ in charges, namely anions are divalent and cations are monovalent but they are of the same diameter d = 400 pm. The temperature analysis of heat capacity at a constant volume Cv and the anion radial distribution function, g2-/2-, allowed the choice of temperature of the study, which is T = 2800 K and corresponds to T(*) = 0.34 (definition of reduced temperature T(*) in text). The differential capacitance curve of the interface with the molten salt or ionic liquid at c = 5.79 M has a distorted bell shape. It is shown that with increasing electrolyte concentration from c = 0.4 to 5 M the differential capacitance curves undergo transition from U shape to bell shape.

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.

On the isobaric thermal expansivity of liquids

The temperature and pressure dependence of isobaric thermal expansivity, α(p), in liquids is discussed in this paper. Reported literature data allow general trends in this property that are consistent with experimental evidence to be established. Thus, a negative pressure dependence is to be expected except around the critical point. On the other hand, α(p) exhibits broad regions of negative and positive temperature dependence in the (T, p) plane depending on the nature of the particular liquid. These trends are rationalized here in terms of various molecular-based equations of state. The analysis of the Lennard-Jones, hard sphere square well and restricted primitive model equations allows understanding the differences in the α(p) behavior between liquids of diverse chemical nature (polar, nonpolar, and ionic): broader regions of negative temperature and positive pressure dependencies are obtained for liquids characterized by larger ranges of the interparticle potential. Also, using the statistical associating fluid theory (SAFT) allowed the behavior of more complex systems (basically, those potentially involving chain and association effects) to be described. The effect of chain length is rather simple: increasing it is apparently equivalent to raise the interaction range. By contrast, association presents a quite complex effect on α(p), which comes from a balance between the dispersive and associative parts of the interaction potential. Thus, if SAFT parameters are adjusted to obtain low association ability, α(p) is affected by each mechanism at clearly separate regions, one at low temperature, due to association, and the other to dispersive forces, which has its origin in fluctuations related with vapor-liquid transition.

Complete phase behavior of the symmetrical colloidal electrolyte

We computed the complete phase diagram of the symmetrical colloidal electrolyte by means of Monte Carlo simulations. Thermodynamic integration, together with the Einstein-crystal method, and Gibbs-Duhem integration were used to calculate the equilibrium phase behavior. The system was modeled via the linear screening theory, where the electrostatic interactions are screened by the presence of salt in the medium, characterized by the inverse Debye length, kappa (in this work kappasigma=6). Our results show that at high temperature, the hard-sphere picture is recovered, i.e., the liquid crystallizes into a fcc crystal that does not exhibit charge ordering. In the low temperature region, the liquid freezes into a CsCl structure because charge correlations enhance the pairing between oppositely charged colloids, making the liquid-gas transition metastable with respect to crystallization. Upon increasing density, the CsCl solid transforms into a CuAu-like crystal and this one, in turn, transforms into a tetragonal ordered crystal near close packing. Finally, we have studied the ordered-disordered transitions finding three triple points where the phases in coexistence are liquid-CsCl-disordered fcc, CsCl-CuAu-disordered fcc, and CuAu-tetragonal-disordered fcc.

Solid phase thermodynamic perturbation theory: Test and application to multiple solid phases

A simple procedure for the determination of hard sphere (HS) solid phase radial distribution function (rdf) is proposed, which, thanks to its physical foundation, allows for extension to other crystal structures besides the fcc structure. The validity of the procedure is confirmed by comparing (1) the predicted HS solid phase rdf's with corresponding simulation data and (2) the predicted non-HS solid phase Helmholtz free energy by the present solid phase first-order thermodynamic perturbation theory (TPT) whose numerical implementation depends on the HS solid phase rdf's as input, with the corresponding predictions also by the first-order TPT but the required HS solid phase rdf is given by an "exact" empirical simulation-fitted formula. The present solid phase first-order TPT predicts isostructural fcc-fcc transition of a hard core attractive Yukawa fluid, in very satisfactory agreement with the corresponding simulation data and is far more accurate than a recent thermodynamically consistent density functional perturbation theory. The present solid phase first-order TPT is employed to investigate multiple solid phases. It is found that a short-ranged potential, even if it is continuous and differentiable or is superimposed over a long-ranged potential, is sufficient to induce the multiple solid phases. When the potential range is short enough, not only isostructural fcc-fcc transition but also isostructural bcc-bcc transition, simple cubic (sc)-sc transition, or even fcc-bcc, fcc-sc, and bcc-sc transitions can be induced. Even triple point involving three solid phases becomes possible. The multiple solid phases can be stable or metastable depending on the potential parameters.

The melting point of ice Ih for common water models calculated from direct coexistence of the solid-liquid interface

In this work we present an implementation for the calculation of the melting point of ice I(h) from direct coexistence of the solid-liquid interface. We use molecular dynamics simulations of boxes containing liquid water and ice in contact. The implementation is based on the analysis of the evolution of the total energy along NpT simulations at different temperatures. We report the calculation of the melting point of ice I(h) at 1 bar for seven water models: SPC/E, TIP4P, TIP4P-Ew, TIP4P/ice, TIP4P/2005, TIP5P, and TIP5P-E. The results for the melting temperature from the direct coexistence simulations of this work are in agreement (within the statistical uncertainty) with those obtained previously by us from free energy calculations. By taking into account the results of this work and those of our free energy calculations, recommended values of the melting point of ice I(h) at 1 bar for the above mentioned water models are provided.

Gas-liquid phase separation in oppositely charged colloids: Stability and interfacial tension

We study the phase behavior and the interfacial tension of the screened Coulomb (Yukawa) restricted primitive model (YRPM) of oppositely charged hard spheres with diameter sigma using Monte Carlo simulations. We determine the gas-liquid and gas-solid phase transitions using free energy calculations and grand-canonical Monte Carlo simulations for varying inverse Debye screening length kappa. We find that the gas-liquid phase separation is stable for kappasigma<or=4, and that the critical temperature decreases upon increasing the screening of the interaction (decreasing the range of the interaction). In addition, we determine the gas-liquid interfacial tension using grand-canonical Monte Carlo simulations. The interfacial tension decreases upon increasing the range of the interaction. In particular, we find that simple scaling can be used to relate the interfacial tension of the YRPM to that of the restricted primitive model, where particles interact with bare Coulomb interactions.

Phase behavior of the lattice restricted primitive model with nearest neighbor exclusion

The global phase behavior of the lattice restricted primitive model with nearest neighbor exclusion has been studied by grand canonical Monte Carlo simulations. The phase diagram is dominated by a fluid (or charge-disordered solid) to charge-ordered solid transition that terminates at the maximum density rho*(max)= sqrt 2 and reduced temperature T* approximately equal to 0.29. At that point, there is a first-order phase transition between two phases of the same density, one charge-ordered, and the other charge-disordered. The liquid-vapor transition for the model is metastable, lying entirely within the fluid-solid phase envelope.

CuAu structure in the restricted primitive model and oppositely charged colloids

We study the phase behavior of oppositely charged equal-size hard spheres both theoretically and experimentally, using Monte Carlo simulations and confocal microscopy. In the simulations, two systems are considered: the restricted primitive model (RPM) and a system of screened Coulomb particles. We construct the phase diagrams of both systems by computer simulations and predict a novel solid phase that has the CuAu structure. In addition, the CuAu structure is observed experimentally in a system of oppositely charged colloids. The qualitative agreement between the RPM, the screened Coulomb system, and the experiments shows that colloids form a suitable model system to study phase behavior in ionic systems.

Three-dimensional binary superlattices of oppositely charged colloids

We report the equilibrium self-assembly of binary crystals of oppositely charged colloidal microspheres at high density. By varying the magnitude of the charge on near equal-sized spheres we show that the structure of the binary crystal may be switched between face-centered cubic, cesium chloride, and sodium chloride. We interpret these transformations in terms of a competition between entropic and Coulombic forces.

Molecular dynamics simulations of the surface tension of ionic liquids

We report molecular dynamics computer simulations of the surface tension and interfacial thickness of ionic liquid-vapor interfaces modeled with a soft core primitive model potential. We find that the surface tension shows an anomalous oscillatory behavior with interfacial area. This observation is discussed in terms of finite size effects introduced by the periodic boundary conditions employed in computer simulations. Otherwise we show that the thickness of the liquid-vapor interface increases with surface area as predicted by the capillary wave theory. Data on the surface tension of size-asymmetric ionic liquids are reported and compared with experimental data of molten salts. Our data suggest that the surface tensions of size-asymmetric ionic liquids do not follow a corresponding states law.

Characterization of the order-disorder transition of a charged hard-sphere model

Monte Carlo simulations at constant pressure are used to characterize the structure of the restricted primitive model in the tetragonal-ordered solid phase. A method to estimate the location of the order-disorder transition and the densities of the coexistence phases is discussed. The results support the weakly first-order character of the transition.

Surface tension of the restrictive primitive model for ionic liquids

Hybrid molecular dynamics and Monte Carlo simulations are performed to study the liquid-vapor interface of the restricted primitive model (RPM) of ionic fluids. We report for the first time simulation results of the surface tension associated to this interface. The RPM accurately predicts experimental surface tensions of ionic salts and good agreement with theoretical predictions that include the idea of ion association is found. The simulation results indicate that the structure of an ionic liquid-vapor interface is rather rough. This is reflected in the interfacial thickness, larger than that observed in simple fluids and water.

Critical behavior of ionic solids

Phase transitions of lattice models of ionic crystals are studied by computer simulation. The nature of order-disorder transitions on different crystal structures is established and compared with the behavior of related Ising models. It is found that for both, continuous and first order transitions the basic features seem to be similar to those of Ising systems.

Order-disorder transition in the solid phase of a charged hard sphere model

We investigate the solid phases of the restricted primitive model (RPM). Monte Carlo simulations show the existence of an order-disorder transition from a substitutionally disordered face centered cubic lattice (fcc) to a new ordered fcc structure which is proposed as the ground state of the RPM at the close packing density. Our results suggest that the new phase might turn out in a new triple point in the RPM phase diagram involving three solid phases: CsCl, fcc ordered and fcc disordered structures. The order-disorder transition is also studied using the cell theory. The theory shows good agreement with the simulation results and suggests that the transition is weakly first order.