Complete phase behavior of the symmetrical colloidal electrolyte

Calculation Methods of Solution Chemical Potential and Application in Emulsion Microencapsulation

Several new biased sampling methods were summarized for solution chemical potential calculation methods in the field of emulsion microencapsulation. The principles, features, and calculation efficiencies of various biased Widom insertion sampling methods were introduced, including volume detection bias, simulation ensemble bias, and particle insertion bias. The proper matches between various types of solution in emulsion and biased Widom methods were suggested, following detailed analyses on the biased insertion techniques. The volume detection bias methods effectively improved the accuracy of the data and the calculation efficiency by inserting detection particles and were suggested to be used for the calculation of solvent chemical potential for the homogeneous aqueous phase of the emulsion. The chemical potential of water, argon, and fluorobenzene (a typical solvent of the oil phase in double emulsion) was calculated by a new, optimized volume detection bias proposed by this work. The recently developed Well-Tempered(WT)-Metadynamics method skillfully constructed low-density regions for particle insertion and dynamically adjusted the system configuration according to the potential energy around the detection point, and hence, could be used for the oil-polymer mixtures of microencapsulation emulsion. For the macromolecule solutes in the oil or aqueous phase of the emulsion, the particle insertion bias could be applied to greatly increase the success rate of Widom insertions. Readers were expected to choose appropriate biased Widom methods to carry out their calculations on chemical potential, fugacity, and solubility of solutions based on the system molecular properties, inspired by this paper.

Parasitic crystallization of colloidal electrolytes: growing a metastable crystal from the nucleus of a stable phase

Colloidal particles have been extensively used to comprehend the main principles governing liquid-crystal nucleation. Multiple mechanisms and frameworks have been proposed, through either experiments or computational approaches, to rationalise the ubiquitous formation of colloidal crystals. In this work, we elucidate the nucleation scenario behind the crystallization of oppositely charged colloids. By performing molecular dynamics simulations of colloidal electrolytes in combination with the Seeding technique, we evaluate the fundamental factors, such as the nucleation rate, free energy barrier, surface tension and kinetic pre-factor, that determine the liquid-to-solid transition of several crystalline polymorphs. Our results show that at a high packing fraction, there is a cross-over between the nucleation of the CsCl structure and that of a substitutionally disordered fcc phase, despite the CuAu crystal being the most stable phase. We demonstrate that the crucial factor in determining which phase nucleates the fastest is the free energy cost of the cluster formation rather than their kinetic ability to grow from the liquid. While at a low packing fraction, the stable phase, CsCl, is the one that nucleates and subsequently grows, we show how at moderate and high packing fractions, a disordered fcc phase subsequently grows regardless of the nature of the nucleating phase, termed parasitic crystallization. Taken together, our results provide a panoramic perspective of the complex nucleation scenario of oppositely charged colloids at moderate temperature and rationalise the different thermodynamic and kinetic aspects behind it.

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.

Computing the free energy of molecular solids by the Einstein molecule approach: ices XIII and XIV, hard-dumbbells and a patchy model of proteins

The recently proposed Einstein molecule approach is extended to compute the free energy of molecular solids. This method is a variant of the Einstein crystal method of Frenkel and Ladd [J. Chem. Phys. 81, 3188 (1984)]. In order to show its applicability, we have computed the free energy of a hard-dumbbell solid, of two recently discovered solid phases of water, namely, ice XIII and ice XIV, where the interactions between water molecules are described by the rigid nonpolarizable TIP4P/2005 model potential, and of several solid phases that are thermodynamically stable for an anisotropic patchy model with octahedral symmetry which mimics proteins. Our calculations show that both the Einstein crystal method and the Einstein molecule approach yield the same results within statistical uncertainty. In addition, we have studied in detail some subtle issues concerning the calculation of the free energy of molecular solids. First, for solids with noncubic symmetry, we have studied the effect of the shape of the simulation box on the free energy. Our results show that the equilibrium shape of the simulation box must be used to compute the free energy in order to avoid the appearance of artificial stress in the system that will result in an increase in the free energy. In complex solids, such as the solid phases of water, another difficulty is related to the choice of the reference structure. As in some cases there is no obvious orientation of the molecules; it is not clear how to generate the reference structure. Our results will show that, as long as the structure is not too far from the equilibrium structure, the calculated free energy is invariant to the reference structure used in the free energy calculations. Finally, the strong size dependence of the free energy of solids is also studied.