Solid-solid and liquid-solid phase equilibria for the restricted primitive 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.

Anomalously large capacitance of an ionic liquid described by the restricted primitive model

We use Monte Carlo simulations to examine the simplest model of a room-temperature ionic liquid (RTIL), called the "restricted primitive model," at a metal surface. We find that at moderately low temperatures the capacitance of the metal-RTIL interface is so large that the effective thickness of the electrostatic double layer is up to three times smaller than the ion radius. To interpret these results we suggest an approach which is based on the interaction between discrete ions and their image charges in the metal surface and which therefore goes beyond the mean-field approximation. When a voltage is applied across the interface, the strong image attraction causes counterions to condense onto the metal surface to form compact ion-image dipoles. These dipoles repel each other to form a correlated liquid. When the surface density of these dipoles is low, the insertion of an additional dipole does not require much energy. This leads to a large capacitance C that decreases monotonically with voltage V, producing a "bell-shaped" curve C(V). We also consider what happens when the electrode is made from a semimetal rather than a perfect metal. In this case, the finite screening radius of the electrode shifts the reflection plane for image charges to the interior of the electrode, and we arrive at a "camel-shaped" C(V). These predictions seem to be in qualitative agreement with experiment.

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.

Oppositely charged colloidal binary mixtures: A colloidal analog of the restricted primitive model

The equilibrium phase diagram of a colloidal system composed of 1:1 mixture of positive and negative particles with equal charge is studied by means of Monte Carlo simulations. The system is the colloidal analog of the restricted primitive model (RPM) for ionic fluids. A liquid-gas transition is found in the low-temperature-low-density region, similar to the liquid-gas transition in the RPM. The fluid-crystal transition is also studied, and the liquid phase is shown to be stable in a narrow range of temperatures. In the liquid, the pair distribution function shows alternating layers of particles with opposite sign of charge surrounding every particle. In the vapor phase, clusters of particles are observed, again in agreement with the RPM. However, a decreasing distribution of clusters is obtained, instead of the discrimination between charged and neutral clusters found in the RPM.

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.

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.