Nucleation of self-associating fluids: free versus activated association

We use the density functional theory of statistical mechanics in a square gradient approximation to analyze the structure, size, and work of formation of critical nuclei in self-associating fluids where association reduces the strength of the interactions between bonded particles. This effect is expected in systems of strongly dipolar particles that associate into chains. In this work we analyze the nucleation behavior of two types of self-associating fluids: a system comprised of particles that can freely associate, and a system in which the association process involves a thermally activated initiation step. For the first case, we explore the properties of critical nuclei in fluids that exhibit a metastable critical point between a vapor phase and a highly associated liquid phase. In fluids where the association dynamics involves an initiation step, we investigate the nucleation behavior in the vicinity of the polymerization transition. In both cases critical nuclei undergo a structural transition that shares many of the features of the coil-globule transition reported in Monte Carlo simulations of strongly dipolar Stockmayer fluids. Our results suggest that the sharp structural transition observed in these simulations is evidence of the existence of a second-order or nearly second-order association transition in these model fluids.     

Phase separations in liquid crystal-colloid mixtures

We present a mean-field theory to describe phase separations in mixtures of a nematic liquid crystal and a colloidal particle. The theory takes into account an orientational ordering of liquid crystals and a crystalline ordering of colloidal particles. We calculate phase diagrams on the temperature-concentration plane, depending on interactions between a liquid crystal and a colloidal surface and a coupling between nematic and crystalline ordering. We find various phase separation processes, such as a nematic-crystal phase separation and nematic-isotropic-crystal triple point. Inside binodal curves, we find new unstable and metastable regions which are important in phase ordering dynamics. We also find a stable nematic-crystalline (NC) phase, where colloidal particles dispersed in a nematic phase can form a crystalline structure. The coexistence between two NC phases with different concentrations can be appear though the coupling between nematic and crystalline ordering.     

Metastable extension of the liquid-vapor phase equilibrium curve and surface tension

The method of molecular dynamics has been used to calculate the parameters of liquid-vapor phase equilibrium and the surface tension in a two-phase system of 4096 Lennard-Jones particles. Calculations have been made in a range from the triple point to near-critical temperature and also at temperatures below the triple point corresponding to the metastable equilibrium of a supercooled liquid and supersaturated vapor. To determine the surface tension, along with a mechanical approach a thermodynamic one has been used as well. The latter was based on calculation of the excess internal energy of an interfacial layer. It has been shown that in accuracy the thermodynamic approach is as good as the more sophisticated mechanical one. Low-temperature asymptotics of the phase-equilibrium curve and also of liquid and vapor spinodals have been considered in the Lennard-Jones and the van der Waals models. The behavior of the surface tension and the excess internal energy of an interfacial layer at T-->0 is discussed.     

Thermodynamics and kinetics of bubble nucleation: Simulation methodology

The simulation of homogeneous liquid to vapor nucleation is investigated using three rare-event algorithms, boxed molecular dynamics, hybrid umbrella sampling Monte Carlo, and forward flux sampling. Using novel implementations of these methods for efficient use in the isothermal-isobaric ensemble, the free energy barrier to nucleation and the kinetic rate are obtained for a Lennard-Jones fluid at stretched and at superheated conditions. From the free energy surface mapped as a function of two order parameters, the global density and largest bubble volume, we find that the free energy barrier height is larger when projected over bubble volume. Using a regression analysis of forward flux sampling results, we show that bubble volume is a more ideal reaction coordinate than global density to quantify the progression of the metastable liquid toward the stable vapor phase and the intervening free energy barrier. Contrary to the assumptions of theoretical approaches, we find that the bubble takes on cohesive non-spherical shapes with irregular and (sometimes highly) undulating surfaces. Overall, the resulting free energy barriers and rates agree well between the methods, providing a set of complementary algorithms useful for studies of different types of nucleation events.     

Phase behavior of density-dependent pair potentials

Phase diagram is calculated by a recently proposed third-order thermodynamic perturbation theory (TPT) for fluid phase and a recently proposed first-order TPT for solid phases; the underlying interparticle potential consists of a hard sphere repulsion and a perturbation tail of an attractive inverse power law type or Yukawa type whose range varies with bulk densities. It is found that besides usual phase transitions associated with density-independent potentials, the density dependence of the perturbation tail evokes some additional novel phase transitions including isostructural solid-solid transition and liquid-liquid transition. Novel triple points are also exhibited which includes stable fluid (vapor or liquid)-face-centered cubic(fcc)-fcc and liquid-liquid-fcc, metastable liquid-body-centered cubic(bcc)-bcc. It also is found that the phase diagram sensitively depends on the density dependence and the concrete mathematical form of the underlying potentials. Some of the disclosed novel transitions has been observed experimentally in complex fluids and molecular liquids, while others still remain to be experimentally verified.     

Phase behavior of C60 by computer simulation using ab‐initio interaction potential

A first‐principles intermolecular potential recently proposed by Pacheco and Ramalho [Phys Rev Lett 1997, 79, 3873–3876] has been used with the Gibbs ensemble and Gibbs–Duhem integration Monte Carlo methods to simulate the vapor–liquid and fluid–solid coexistence properties of C₆₀. The critical properties were calculated by fitting the results to the laws of rectilinear diameters and order parameter scaling. The triple‐point properties were determined from the limiting behavior of the Gibbs ensemble vapor–liquid simulations at the lowest temperature range. A stable liquid phase is predicted for temperatures between 1570±20 and 2006±27 K and densities between 0.444±0.003 and 1.05±0.01 nm^?3. The estimated critical and triple‐point pressures are, respectively, 35±6 and 5±16 bars. We show for the first time, to our knowledge, that it is possible, strictly by computer simulation, to estimate a triple point for C₆₀ in accordance with the predictions of theoretical methods and the basic concepts of thermodynamics. The liquid and fluid radial distribution functions indicate the presence of solid or glasslike features. This may support the suggestion of a more cooperative interaction of clusters in C₆₀. A comparison of our results with the data obtained by other authors is presented and discussed. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem 84: 375–387, 2001     

Free-energy landscape of nucleation with an intermediate metastable phase studied using capillarity approximation

Capillarity approximation is used to study the free-energy landscape of nucleation when an intermediate metastable phase exists. The critical nucleus that corresponds to the saddle point of the free-energy landscape as well as the whole free-energy landscape can be studied using this capillarity approximation, and various scenarios of nucleation and growth can be elucidated. In this study, we consider a model in which a stable solid phase nucleates within a metastable vapor phase when an intermediate metastable liquid phase exists. We predict that a composite critical nucleus that consists of a solid core and a liquid wetting layer as well as pure liquid and pure solid critical nuclei can exist depending not only on the supersaturation of the liquid phase relative to that of the vapor phase but also on the wetting behavior of the liquid surrounding the solid. The existence of liquid critical nucleus indicates that the phase transformation from metastable vapor to stable solid occurs via the intermediate metastable liquid phase, which is quite similar to the scenario of nucleation observed in proteins and colloidal systems. By studying the minimum-free-energy path on the free-energy landscape, we can study the evolution of the composition of solid and liquid within nuclei which is not limited to the critical nucleus.     

受限于纳米狭缝中的四缔合Lennard-Lones流体的相平衡和界面张力研究

分别用改进的基础测量理论和平均球近似理论表达短程作用和长程作用对四缔合Lennard-Jones流体的过剩自由能的贡献. 在密度函泛理论的框架下, 研究了平均密度等温线, 密度分布, 未缔合分子在平衡汽相和液相中的分布, 相平衡以及平衡时的界面张力等热力学性质. 分析了缔合能量, 流体-固体作用和孔宽对受限于纳米狭缝中的四缔合Lennard-Jones流体相行为的影响.     

Re-entrant phase behavior for systems with competition between phase separation and self-assembly

In patchy particle systems where there is a competition between the self-assembly of finite clusters and liquid-vapor phase separation, re-entrant phase behavior can be observed, with the system passing from a monomeric vapor phase to a region of liquid-vapor phase coexistence and then to a vapor phase of clusters as the temperature is decreased at constant density. Here, we present a classical statistical mechanical approach to the determination of the complete phase diagram of such a system. We model the system as a van der Waals fluid, but one where the monomers can assemble into monodisperse clusters that have no attractive interactions with any of the other species. The resulting phase diagrams show a clear region of re-entrance. However, for the most physically reasonable parameter values of the model, this behavior is restricted to a certain range of density, with phase separation still persisting at high densities.     

The intrinsic structure of the water surface

An operational procedure to obtain the intrinsic structure of liquid surfaces is applied here to a molecular dynamics simulation of water, with a model of point charges for the molecular interactions. The method, which had been recently proposed and used for simple fluids, is successfully extended to a molecular liquid with the complex bond structure of water. The elimination of the capillary wave fluctuations, in the intrinsic density and orientation profiles, gives a new overall view of the water surface, at the sharpest molecular level, and without the size-dependent broadening observed in the mean profiles. The molecules belonging to the outer liquid layer are clearly identified, and we find that only these molecules exhibit a clear preferential orientation to lie flat on the surface. Moreover, there is a strong correlation between the dipolar structure and the local curvatures of the intrinsic surface, so that at the extrusions of the intrinsic surface the molecular dipoles point preferentially toward the vapor side of the interface. Finally, we have found an intrinsic density layering structure, although the inner structure is strongly damped beyond the second layer.     

Experimental phase diagram of symmetric binary colloidal mixtures with opposite charges

The phase behavior of equimolar mixtures of oppositely charged colloidal systems with similar absolute charges is studied experimentally as a function of the salt concentration in the system and the colloid volume fraction. As the salt concentration increases, fluids of irreversible clusters, gels, liquid-gas coexistence, and finally, homogeneous fluids, are observed. Previous simulations of similar mixtures of Derjaguin-Landau-Verwey-Overbeek (DLVO) particles indeed showed the transition from homogeneous fluids to liquid-gas separation, but also predicted a reentrant fluid phase at low salt concentrations, which is not found in the experiments. Possibly, the fluid of clusters could be caused by a nonergodicity transition responsible for the gel phase in the reentrant fluid phase. Liquid-gas separation takes a delay time after the sample is prepared, whereas gels collapse from the beginning. The density of the liquid in coexistence with a vapor phase depends linearly on the overall colloid density of the system. The vapor, on the other hand, is comprised of equilibrium clusters, as expected from the simulations.     

Stability of nematic and smectic phases in rod-like mesogens with orientation-dependent attractive interactions

The stability of isotropic (I), nematic (N), smectic A (Sm A), and hexatic (Hex) liquid crystalline phases is studied for a fluid of molecules with a rod-like shape and dispersive interactions dependent on orientation. The fluid is modeled with the spherocylindrical Gay-Berne-Kihara interaction potential proposed in a recent work, with parameters favoring parallel pair orientations. The liquid crystal phase diagram is characterized for different molecular aspect ratios by means of Monte Carlo simulations in the isobaric-isothermal ensemble. Three types of triple points are observed, namely, I-Sm A-Hex, I-N-Sm A, and N-Sm A-Hex, leading to island-shape domains for the smectic A phase. The resulting phase diagrams are compared with those derived previously for prolate fluids of ellipsoidal and spherocylindrical symmetry. It is concluded that the stability of the layered phases with respect to the nematic phase is enhanced in the spherocylindrical fluids due to geometrical constraints. Furthermore, the anisotropy of the dispersive interactions induces a stronger dependence of the overall phase diagram on temperature and aids in the energetic stabilization of the hexatic crystalline phase with respect to the fluid smectic A phase.     

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.     

Direct determination of phase behavior of square-well fluids

We have combined Gibbs ensemble Monte Carlo simulations with the aggregation volume-biased method in conjunction with the Gibbs-Duhem method to provide the first direct estimates for the vapor-solid, vapor-liquid, and liquid-solid phase coexistences of square-well fluids with three different ranges of attraction. Our results are consistent with the previous simulations and verify the notion that the vapor-liquid coexistence behavior becomes metastable for cases where the attraction well becomes smaller than 1.25 times the particle diameter. In these cases no triple point is found.     

Extended Veytsman statistics for the solid–fluid transition in the framework of lattice fluid

Lattice fluid can describe a vapor–liquid transition but not a solid–fluid transition. In this work, we propose a simple and analytic term which yields a solid–fluid transition when coupled with a lattice based equation of state (EOS). The proposed term is derived based on the two assumptions that (1) solid can be considered as highly associated phase affected by strong attractive force and (2) this force is distinct from the conventional attractive forces yielding a vapor–liquid transition. To formulate these assumptions, we extend Veytsman statistics by modifying its density dependency. The derived term was combined with a quasi-chemical nonrandom lattice fluid theory (QLF) developed by the authors. The combined model was found to require only two parameters besides 3 QLF parameters for physical properties calculation of three phases. When tested against equilibrium properties of 8 components, the combined model was found to closely reproduce melting pressure, sublimation pressure, and vapor pressure, but underestimate solid density as well as heat of melting at the triple point temperature. It was found that the present approach can yield a solid–liquid transition at all temperatures.     

Vapor-liquid and vapor-solid phase equilibria for united-atom benzene models near their triple points: the importance of quadrupolar interactions

Gibbs ensemble Monte Carlo simulations were used to calculate the vapor-liquid and vapor-solid coexistence curves for benzene using two simple united-atom models. An extension of the Gibbs ensemble method that makes use of an elongated box containing a slab of the condensed phase with a vapor phase along one axis was employed for the simulations of the vapor-solid equilibria and the vapor-liquid equilibria at very low reduced temperatures. Configurational-bias and aggregation-volume-bias Monte Carlo techniques were applied to improve the sampling of particle transfers between the two simulation boxes and between the vapor and condensed-phase regions of the elongated box. An isotropic united-atom representation with six Lennard-Jones sites at the positions of the carbon atoms was used for both force fields, but one model contained three additional out-of-plane partial charge sites to explicitly represent benzene's quadrupolar interactions. Both models were fitted to reproduce the critical temperature and density of benzene and yield a fair representation of the vapor-liquid coexistence curve. In contrast, differences between the models are very large for the vapor-solid coexistence curve. In particular, the lack of explicit quadrupolar interactions for the 6-site model greatly reduces the energetic differences between liquid and solid phases, and this model yields a triple point temperature that is about a factor of 2 too low. In contrast, the 9-site model predicts a triple point of benzene at T = 253 +/- 6 K and p = 2.3 +/- 0.8 kPa in satisfactory agreement with the experimental data (T = 278.7 K and p = 4.785 kPa).     

On the universal behavior of some thermodynamic properties along the whole liquid-vapor coexistence curve

When thermodynamic properties of a pure substance are transformed to reduced form by using both critical- and triple-point values, the corresponding experimental data along the whole liquid-vapor coexistence curve can be correlated with a very simple analytical expression that interpolates between the behavior near the triple and the critical points. The leading terms of this expression contain only two parameters: the critical exponent and the slope at the triple point. For a given thermodynamic property, the critical exponent has a universal character but the slope at the triple point can vary significantly from one substance to another. However, for certain thermodynamic properties including the difference of coexisting densities, the enthalpy of vaporization, and the surface tension of the saturated liquid, one finds that the slope at the triple point also has a nearly universal value for a wide class of fluids. These thermodynamic properties thus show a corresponding apparently universal behavior along the whole coexistence curve.     

The interaction of the theophylline metastable phase with water vapor

The conditions of hydration of the stable and metastable theophylline phases were determined. Two-phase metastable phase/monohydrate and stable phase/monohydrate equilibrium pressures were measured at 25, 30, and 35°C. The metastable phase began to react with water vapor at lower relative humidities than the stable phase. Processes that occurred with the metastable and stable theophylline phases over various water pressure ranges were considered. The metastable phase exhibited an unusual behavior at 25°C and relative humidity 47%. At constant water vapor pressure and temperature, theophylline was initially hydrated and then lost water and again became anhydrous. Two consecutive processes occurred in the system, the formation of theophylline monohydrate from the metastable phase and its decomposition to the stable phase. The ratio between the rates of these processes determined the content of the monohydrate at the given time moment.