A study of the phase behavior of a simple model of chiral molecules and enantiomeric mixtures

We present a study of the solid-fluid and solid-solid phase equilibrium for molecular models representative of chiral molecules and enantiomeric mixtures. The models consist of four hard sphere interaction sites of different diameters in a tetrahedral arrangement with the fifth hard sphere interaction site at the center of the tetrahedron. The volumetric properties and free energies of the pure enantiomers and binary mixtures were calculated in both fluid and solid phases using isobaric Monte Carlo simulations. The models exhibit essentially ideal solution behavior in the fluid phase with little chiral discrimination. In the solid phase the effects of chirality are much greater. Solid-fluid phase behavior involving the pure enantiomer solids and also racemic compounds was calculated. The calculations indicate that, depending on the relative sizes of the hard sphere interaction sites, packing effects alone can be sufficient to stabilize a racemic compound with respect to the pure enantiomer solids.     

Effect of charge asymmetry and charge screening on structure of superlattices formed by oppositely charged colloidal particles

Colloidal suspensions made up of oppositely charged particles have been shown to self-assemble into substitutionally ordered superlattices. For a given colloidal suspension, the structure of the superlattice formed from self-assembly depends on its composition, charges on the particles, and charge screening. In this study we have computed the pressure-composition phase diagrams of colloidal suspensions made up of binary mixtures of equal sized and oppositely charged particles interacting via hard core Yukawa potential for varying values of charge screening and charge asymmetry. The systems are studied under conditions where the thermal energy is equal or greater in magnitude to the contact energy of the particles and the Debye screening length is smaller than the size of the particles. Our studies show that charge asymmetry has a significant effect on the ability of colloidal suspensions to form substitutionally ordered superlattices. Slight deviations of the charges from the stoichiometric ratio are found to drastically reduce the thermodynamic stability of substitutionally ordered superlattices. These studies also show that for equal-sized particles, there is an optimum amount of charge screening that favors the formation of substitutionally ordered superlattices.     

Solid-fluid and solid-solid phase equilibrium in a model of n-alkane mixtures

Solid-fluid and solid-solid phase equilibrium for binary mixtures of hard sphere chains modeling n-hexane, n-heptane, and n-octane has been calculated using Monte Carlo computer simulations. Thermodynamic integration was used to calculate the Gibbs free energy and chemical potentials in the solid and fluid phases from pure component reference values. A multiple stage free energy perturbation method was used to calculate the composition derivative of the Gibbs free energy. Equation of state and free energy data for the fluid phase indicate ideal solution behavior. Nonideality is much more significant in the solid phase with only partial solubility of shorter chains in the longer chains and essentially no solubility at the other end of the composition range. The miscibility decreases with increasing chain length difference between the components. For the model of n-hexane/n-octane mixtures solid--solid phase separation has been observed directly in some of the simulations, with the components segregating between the layers of the solid structure. The behavior is similar to that seen in some binary n-alkane mixtures with longer chain lengths but comparable chain length ratios between the components. Such phase separation, although indicated thermodynamically, is not seen directly in the simulations of the n-heptane/n-octane mixture due to the difference in the pure component crystal structures.     

Phase diagram of hard snowman-shaped particles

We present the phase diagram of hard snowman-shaped particles calculated using Monte Carlo simulations and free energy calculations. The snowman particles consist of two hard spheres rigidly attached at their surfaces. We find a rich phase behavior with isotropic, plastic crystal, and aperiodic crystal phases. The crystalline phases found to be stable for a given sphere diameter ratio correspond mostly to the close packed structures predicted for equimolar binary hard-sphere mixtures of the same diameter ratio. However, our results also show several crystal-crystal phase transitions, with structures with a higher degree of degeneracy found to be stable at lower densities, while those with the best packing are found to be stable at higher densities.     

Effects of Polydispersity on the Phase Behavior of Additive Hard Spheres in Solution

The ability to separate enzymes, nucleic acids, cells, and viruses is an important asset in life sciences. This can be realised by using their spontaneous asymmetric partitioning over two macromolecular aqueous phases in equilibrium with one another. Such phases can already form while mixing two different types of macromolecules in water. We investigate the effect of polydispersity of the macromolecules on the two-phase formation. We study theoretically the phase behavior of a model polydisperse system: an asymmetric binary mixture of hard spheres, of which the smaller component is monodisperse and the larger component is polydisperse. The interactions are modelled in terms of the second virial coefficient and are assumed to be additive hard sphere interactions. The polydisperse component is subdivided into sub-components and has an average size ten times the size of the monodisperse component. We calculate the theoretical liquid–liquid phase separation boundary (the binodal), the critical point, and the spinodal. We vary the distribution of the polydisperse component in terms of skewness, modality, polydispersity, and number of sub-components. We compare the phase behavior of the polydisperse mixtures with their concomittant monodisperse mixtures. We find that the largest species in the larger (polydisperse) component causes the largest shift in the position of the phase boundary, critical point, and spinodal compared to the binary monodisperse binary mixtures. The polydisperse component also shows fractionation. The smaller species of the polydisperse component favor the phase enriched in the smaller component. This phase also has a higher-volume fraction compared to the monodisperse mixture.     

Analysis of retrograde behavior and the cross-over effect in supercritical fluids

In this paper, results based upon thermodynamic stability theory are developed which lead to a set of both necessary and sufficient conditions for the existence of retrograde behavior in a multicomponent solute system dissolved in a pure supercritical fluid. While experimental evidence of retrograde behavior in single solute systems has been known for some time, recently data have been obtained showing the retrograde effect in binary solute systems dissolved in a pure supercritical fluid.In such systems, cross-over regions may be defined. These are pressure—temperature regions where the solubility of one solute increases while that of the other decreases with a change in temperature at constant pressure. The existence of cross-over regions in multicomponent mixtures can have implications for a new separation technique using pure supercritical fluids. In conjunction with an equation of state, the results derived here allow cross-over regions to be predicted, thus enabling one to identify candidate systems and thermodynamic conditions for the cross-over process. For this work a variation of a perturbed hard sphere model equation of state was used for the calculations.     

Generic pathways to stability in concentrated protein mixtures

We present a series of experimental results that disclose the crucial role of ionic strength and partial volume fractions in the control of the phase behaviour of binary protein mixtures. Our findings can be understood as that the ionic strength determines the relative contribution of the entropy of the protein counter-ions to the overall thermodynamics of the system. Associative phase separation and crystallization observed at, respectively, low and high ionic strength are suppressed at intermediate salt concentrations, where the entropy gain upon releasing the counter-ions from the double layer of the proteins is negligible and the entropy loss upon confining the counter-ions within the protein crystal phase significant. Moreover, we find that the partial volume fraction of the protein prone to crystallize determines the crystallization boundary and that the presence of other proteins strongly delays crystallization, leading to temporarily stable mixtures. These findings suggest that stability in more complex protein mixtures, such as the cytosol, relates to the ionic strength and protein composition rather than to protein specific properties.     

Polymorphism in AB(13) nanoparticle superlattices: an example of semiconductor-metal metamaterials

Colloidal crystallization of nanoparticles with different functionalities into multicomponent assemblies provides a route to new classes of ordered nanocomposites with novel properties tunable by the choice of the constituent building blocks. While theories based on hard sphere approximation predict crystallization of only a few stable binary phases (NaCl-, AlB(2)- and NaZn(13)-type), we find that additional polymorphs of lower packing density are possible. We demonstrate that PbSe and Pd nanoparticles can be reproducibly crystallized into two polymorphs with AB(13) stoichiometry. One polymorph is isostructural with the intermetallic compound NaZn(13) and is consistent with dense packing of hard spheres driven by entropy. The second unanticipated polymorph is of lower packing density. This observation underscores the shortcomings of applying simple space-filling principles to the crystallization of organically passivated nanocrystals and further motivates the development of models that incorporate combinations of hard-sphere, van der Waals, dipolar, and hydrophobic forces. This work demonstrates that ordered periodic structures with lower packing density are achievable and provides the first example of a binary semiconductor-metal superlattice using a combination of PbSe-Pd nanocrystals.     

Phase field theory of interfaces and crystal nucleation in a eutectic system of fcc structure: I. Transitions in the one-phase liquid region

The phase field theory (PFT) has been applied to predict equilibrium interfacial properties and nucleation barrier in the binary eutectic system Ag-Cu using double well and interpolation functions deduced from a Ginzburg-Landau expansion that considers fcc (face centered cubic) crystal symmetries. The temperature and composition dependent free energies of the liquid and solid phases are taken from CALculation of PHAse Diagrams-type calculations. The model parameters of PFT are fixed so as to recover an interface thickness of approximately 1 nm from molecular dynamics simulations and the interfacial free energies from the experimental dihedral angles available for the pure components. A nontrivial temperature and composition dependence for the equilibrium interfacial free energy is observed. Mapping the possible nucleation pathways, we find that the Ag and Cu rich critical fluctuations compete against each other in the neighborhood of the eutectic composition. The Tolman length is positive and shows a maximum as a function of undercooling. The PFT predictions for the critical undercooling are found to be consistent with experimental results. These results support the view that heterogeneous nucleation took place in the undercooling experiments available at present. We also present calculations using the classical droplet model [classical nucleation theory (CNT)] and a phenomenological diffuse interface theory (DIT). While the predictions of the CNT with a purely entropic interfacial free energy underestimate the critical undercooling, the DIT results appear to be in a reasonable agreement with the PFT predictions.     

The effects of shape and flexibility on bio-engineered fd-virus suspensions

We present a theoretical model to describe binary mixtures of semi-flexible rods, applied here to fd-virus suspensions. We investigate the effects of rod stiffness on both monodisperse and binary systems, studying thick-thin and long-short mixtures. For monodisperse systems, we find that fd-virus particles have to be made extremely stiff to even approach the behavior of rigid rods. For thick-thin mixtures, we find increasingly rich phase behavior as the rods are either made more flexible or if their diameter ratio is increased. For long-short rod mixtures we find that the phase behavior is controlled by the relative stiffness of the rods, with increasing the stiffness of the long rods or decreasing that of the short rods resulting in richer phase behavior. We also calculate the state point dependent effective shape of the rods. The flexible rods studied here always behave as shorter, thicker rigid rods, but with an effective shape that varies widely throughout the phase diagrams, and plays a key role in determining phase behavior.     

Structure, phase behavior, and inhomogeneous fluid properties of binary dendrimer mixtures

The effective pair potentials between different kinds of dendrimers in solution can be well approximated by appropriate Gaussian functions. We find that in binary dendrimer mixtures the range and strength of the effective interactions depend strongly upon the specific dendrimer architecture. We consider two different types of dendrimer mixtures, employing the Gaussian effective pair potentials, to determine the bulk fluid structure and phase behavior. Using a simple mean field density functional theory (DFT) we find good agreement between theory and simulation results for the bulk fluid structure. Depending on the mixture, we find bulk fluid-fluid phase separation (macrophase separation) or microphase separation, i.e., a transition to a state characterized by undamped periodic concentration fluctuations. We also determine the inhomogeneous fluid structure for confinement in spherical cavities. Again, we find good agreement between the DFT and simulation results. For the dendrimer mixture exhibiting microphase separation, we observe a rather striking pattern formation under confinement.     

Density-functional study of model bidisperse ferrocolloids in an external magnetic field

We present phase diagrams of a model bidisperse ferrocolloid consisting of a binary mixture of dipolar hard spheres (DHSs) under the influence of an external magnetic field. The dipole moments of the particles are chosen proportional to the particle volume to mimic real ferrocolloids, and we focus on dipole-dominated systems where isotropic attractive interactions are absent. Our results are based on density-functional theory in the modified mean-field (MMF) approximation. For one-component DHS fluids in external fields, and for corresponding mixtures dominated by one of the components, MMF theory predicts the tricritical point of the transition between an isotropic gas and a ferromagnetic liquid occurring at zero field to be changed into a critical point separating two magnetically ordered phases of different density. The corresponding critical temperature displays a nonmonotonic dependence on the field strength. Completely different behavior is found for the critical temperature related to the demixing phase transitions appearing in strongly asymmetric mixtures [G. M. Range and S. H. L. Klapp, Phys. Rev. E 70, 061407 (2004)]. For such systems, we find a monotonic decrease of the demixing critical temperature with increasing field. The field strength dependence of the critical temperature can therefore be tuned between nonmonotonic and monotonic behaviors just by changing the composition of the mixture--e.g., by adjusting the chemical potentials. This allows us to efficiently control the influence of external magnetic fields on the phase behavior over a large temperature interval.     

Integral equation theory for symmetric nonadditive hard sphere mixtures

An integral equation theory is presented for the pair correlation functions and phase behavior of symmetric nonadditive hard sphere mixtures with hard sphere diameters given by sigma(A)(A)() = sigma(BB) = lambdad and sigma(AB) = d. This mixture exhibits a fluid-fluid phase separation into an A-rich phase and a B-rich phase at high densities. The theory incorporates, into the closure approximation, all terms that can be calculated exactly in the density expansion of the direct correlation functions. We find that the closure approximation developed in this work is accurate for the structure and phase behavior over the entire range of lambda, when compared to computer simulations, and is significantly more accurate than the previous theories.     

Configurational entropy of binary hard-disk glasses: nonexistence of an ideal glass transition

We study the thermodynamics of a binary hard-disk mixture in which the ratio of disk diameters is kappa=1.4. We use a recently developed molecular dynamics algorithm to calculate the free-volume entropy of glassy configurations and obtain the configurational entropy (degeneracy) of the supercompressed liquid as a function of density. We find that the configurational entropy of the glasses near the kinetic glass transition is very close to the mixing entropy, suggesting that the degeneracy is zero only for the phase-separated crystal. We explicitly construct an exponential number of jammed packings with densities spanning the spectrum from the accepted "amorphous" glassy state to the phase-separated crystal, thus showing that there is no ideal glass transition in binary hard-disk mixtures. This construction also demonstrates that the ideal glass, defined as having zero configurational entropy, is not amorphous, but instead is nothing more than a phase-separated crystal. This critique of the presumed existence of an ideal glass parallels our previous critique of the idea that there is a most-dense random (close) packing for hard spheres [Torquato et al., Phys. Rev. Lett. 84, 2064 (2000)]. We also perform free-energy calculations to determine the equilibrium phase behavior of the system. The calculations predict a first-order freezing transition at a density below the kinetic glass transition. However, this transition appears to be strongly kinetically suppressed and is not observed directly. New simulation techniques are needed in order to gain a more complete understanding of the thermodynamic and kinetic behavior of the binary disk mixture and, in particular, of the demixing process during crystallization.     

On the influence of a patterned substrate on crystallization in suspensions of hard spheres

We present a computer simulation study on crystal nucleation and growth in supersaturated suspensions of mono-disperse hard spheres induced by a triangular lattice substrate. The main result is that compressed substrates are wet by the crystalline phase (the crystalline phase directly appears without any induction time), while for stretched substrates we observe heterogeneous nucleation. The shapes of the nucleated crystallites fluctuate strongly. In the case of homogeneous nucleation amorphous precursors have been observed [T. Schilling et al., Phys. Rev. Lett. 105(2), 025701 (2010)]. For heterogeneous nucleation we do not find such precursors. The fluid is directly transformed into highly ordered crystallites.     

Shape factors in equations of state. Part II. Repulsion phenomena in multicomponent chain fluids

An equation of state for the multicomponent fluid phase of nonattracting rigid particles of arbitrary shape is presented. The equation is a generalization of a previously presented equation of state for pure fluids of rigid particles; the approach describes the volumetric properties of a pure fluid in terms of a shape factor, zeta, which can be back calculated by scaling the volumetric properties of pure fluids to that of a hard sphere. The performance of the proposed equation is tested against mixtures of chain fluids immersed in a "monomeric" solvent of hard spheres of equal and different sizes. Extensive new Monte Carlo simulation data are presented for 19 binary mixtures of hard homonuclear tangent freely-jointed hard sphere chains (pearl-necklace) of various lengths (three to five segments), with spheres of several size ratios and at various compositions. The performance of the proposed equation is compared to the hard-sphere SAFT approach and found to be of comparable accuracy. The equation proposed is further tested for mixtures of spheres with spherocylinders. In all cases, the equation proved to be accurate and simple to use.     

Phase behavior and interfacial properties of nonadditive mixtures of Onsager rods

Within a second virial theory, we study bulk phase diagrams as well as the free planar isotropic-nematic interface of binary mixtures of nonadditive thin and thick hard rods. For species of the same type, the excluded volume is determined only by the dimensions of the particles, whereas for dissimilar ones it is taken to be larger or smaller than that, giving rise to a nonadditivity that can be positive or negative. We argue that such a nonadditivity can result from modeling of soft interactions as effective hard-core interactions. The nonadditivity enhances or reduces the fractionation at isotropic-nematic (IN) coexistence and may induce or suppress a demixing of the high-density nematic phase into two nematic phases of different composition (N(1) and N(2)), depending on whether the nonadditivity is positive or negative. The interfacial tension between coexisting isotropic and nematic phases shows an increase with increasing fractionation at the IN interface, and complete wetting of the IN(2) interface by the N(1) phase upon approach of the triple-point coexistence. In all explored cases bulk and interfacial properties of the nonadditive mixtures exhibit a striking and quite unexpected similarity with the properties of additive mixtures of different diameter ratio.     

PHYSICAL-CHEMICAL BEHAVIOR OF CONFECTIONERY FATS

Confectionery coatings employ hard butter fat systems made from both lauric and non-lauric source oils. These oils are routinely modified by a combination of treatments including hydrogenation, fractionation and interesterification to achieve desired physical properties. Such processing methods create heterogeneous triglyceride mixtures consisting of a variety of compositional and positional isomers. Published phase diagrams of “simple” binary triglyceride mixtures of closely related molecules are complex, and suggest that innumberable unique liquid and solid phases may co-exist at any given temperature ( and pressure ) in vastly more complex triglyceride mixtures such as confectionery hard butters. Thus we may view confectionery fat systems as multiphasic mixtures (liquid, solid and compositional) with a propensity to undergo liquid content fluctuations and crystal size/morphology changes in response to slight changes in temperature. A true equilibrium among all phases may indeed never be attained, and a potential for movement of certain components in response to temperature change is probably constant. Surface growth of long needle-like fat crystals, “fat bloom”, most likely results from this non-equilibrium condition and serves to reduce the system's free energy. The ever present, ever changing liquid phase(s) is viewed as the vehicle for free energy minimization via triglyceride migration and ongoing crystal growth.     

Polymer induced changes of the crystallization scenario in suspensions of hard sphere like microgel particles

We investigated the crystallization scenario of highly cross linked polystyrene particles dispersed in the good solvent 2-ethylnaphtalene and their mixtures with non-adsorbing low molecular weight polysterene polymer using time resolved static light scattering. The samples were prepared slightly below the melting volume fraction of the polymer free system. For the polymer free samples, we obtained polycrystalline solids via crystallization scenario known from hard sphere suspensions with little competition of wall crystal formation. Addition of non-adsorbing low molecular weight polystyrene polymer leads to a considerably slowing down of the bulk crystallization kinetics. We observed a delay of the precursor to crystal conversion for the bulk crystallization while the induction times for the wall nucleation are reduced. The increased polymer concentration thus shifts the balance between the two competing crystallization pathways giving the possibility to tune the relative amount of wall based crystals.     

New phase for one-component hard spheres

A completely new phase for one-component hard spheres is reported in an unexpected region of the phase diagram. The new phase is observed at compressibility factors intermediate between the solid and the metastable branches. It can be obtained from either Monte Carlo simulations alone or a combination of Monte Carlo and molecular dynamics calculations. An analysis of the intermediate scattering function data shows that the new phase is in a stable equilibrium. Radial distribution function data, configurational snapshots, bond order parameters, and translational order parameters obtained from molecular simulations indicate that the new phase is significantly different from the isotropic liquid, metastable, or crystalline phases traditionally observed in hard sphere systems. This result significantly changes our previous understanding of the behavior of hard spheres.