Signatures of critical phenomena in rare atom clusters

We study the occurrence of a critical behaviour related to a liquid-gas phase transition in an atomic cluster. In particular, we analyse some experimental possibilities to investigate the role of finite size effects on the phase behaviour. In mass distributions, conditional moments, as well as scaled factorial moments of fragment mass distributions we find evidence for the presence of critical behaviour. It can be connected, by the use of Fisher's droplet model and Campi analysis, to a liquid-gas phase transition.     

Effect of confinement on the liquid-liquid phase transition of supercooled water

We report on an observation of the phase transition between two liquid phases of supercooled confined water in simulations. The temperature of the liquid-liquid transition of water at zero pressure slightly decreases due to confinement in the hydrophobic pore. The hydrophilic confinement affects this temperature in the opposite direction and shifts the critical point of the liquid-liquid transition to a higher pressure. As a result, in a strongly hydrophilic pore the liquid-liquid phase transition becomes continuous at zero pressure, indicating the shift of its critical point from negative to a positive pressure. These findings indicate that experimental studies of water confined in the pores of various hydrophobicity/hydrophilicity may clarify the location of the liquid-liquid critical point of bulk water.     

Phase separation in mixtures of Yukawa and charged Yukawa particles from Gibbs ensemble Monte Carlo simulations and the mean spherical approximation

The phase equilibrium of mixtures of Yukawa and charged Yukawa particles is studied by means of Gibbs ensemble Monte Carlo (GEMC) simulation method and the mean spherical approximation (MSA). The strength of the Coulomb energy compared to that of the Yukawa attraction is characterized by a coupling constant. For low coupling constants a classical vapor--liquid phase separation appears with a good agreement between GEMC and the MSA. For high coupling constant, a phase separation between a salt poor and a salt rich phase occurs that resembles the phase equilibrium behavior of the solvent primitive model.     

Liquid-liquid phase transitions in supercooled water studied by computer simulations of various water models

Liquid-liquid and liquid-vapor coexistence regions of various water models were determined by Monte Carlo (MC) simulations of isotherms of density fluctuation-restricted systems and by Gibbs ensemble MC simulations. All studied water models show multiple liquid-liquid phase transitions in the supercooled region: we observe two transitions of the TIP4P, TIP5P, and SPCE models and three transitions of the ST2 model. The location of these phase transitions with respect to the liquid-vapor coexistence curve and the glass temperature is highly sensitive to the water model and its implementation. We suggest that the apparent thermodynamic singularity of real liquid water in the supercooled region at about 228 K is caused by an approach to the spinodal of the first (lowest density) liquid-liquid phase transition. The well-known density maximum of liquid water at 277 K is related to the second liquid-liquid phase transition, which is located at positive pressures with a critical point close to the maximum. A possible order parameter and the universality class of liquid-liquid phase transitions in one-component fluids are discussed.     

The thermodynamic properties of Al-Ce melts

Thermodynamic modeling of Al-Ce melts was performed at a total pressure p = 10⁵ Pa in an argon atmosphere over the temperature and concentration ranges 1773–2373 K and 0 ≤ x _Ce ≤ 1, respectively, corresponding to the liquid-phase region of the phase diagram of the system. The concentration and temperature dependences of the thermodynamic characteristics of the melts and the contents of melt components and gas phase components over the melts were calculated, the temperatures, enthalpies, and entropies of the liquid-gas transition were determined, and the phase diagram of the liquid-liquid + gas-gas transition was obtained.     

Phase relations in the lanthanum nitrate (copper nitrate)-poly(vinylpyrrolidone)-water systems

Room-temperature isothermal sections of the phase diagrams of lanthanum nitrate-poly(vinylpyrrolidone)-water and copper nitrate-poly(vinylpyrrolidone)-water systems were studied. The following features were found: a wide region of homogeneous water-polymer solutions, liquid-liquid phase separation field, and a three-phase region in which two liquids coexist with salt crystals. In the lanthanum nitrate system, liquid-liquid phase separation has a lower critical solution point (polythermal sections were studied); in the copper nitrate system, it has an upper critical solution point. The type of diagram for unstudied systems is predicted based on the analysis of polymer-salt phase diagrams.     

Liquid—gas transition for hard spheres with attractive Yukawa tail interactions

We describe the liquid-gas transition in the hard sphere system with Yukawa tail interactions in the mean spherical approximation. The dependence of critical temperature and density on the range of the interaction is shown and the spinodal curve for a short-ranged potential and a long—ranged potential is presented. The compressibility, energy and virial pressures are presented for a long-ranged potential. Liquid phase pressures are calculated by integrating round the coexistence region, rather than through it.     

Vapor-liquid equilibrium of ethanol by molecular dynamics simulation and Voronoi tessellation

Explicit atom simulations of ethanol were performed by molecular dynamics using the OPLS-AA potential. The phase densities were determined self-consistently by comparing the distribution of Voronoi volumes from two-phase and single-phase simulations. This is the first demonstration of the use of Voronoi tessellation in two-phase molecular dynamics simulation of polyatomic fluids. This technique removes all arbitrary determination of the phase diagram by using single-phase simulations to self-consistently validate the probability distribution of Voronoi volumes of the liquid and vapor phases extracted from the two-phase molecular dynamics simulations. Properties from the two phase simulations include critical temperature, critical density, critical pressure, phase diagram, surface tension, and molecule orientation at the interface. The simulations were performed from 375 to 472 K. Also investigated were the vapor pressure and hydrogen bonding along the two phase envelope. The phase envelope agrees extremely well with literature values from GEMC at lower temperatures. The combined use of two-phase molecular dynamics simulation and Voronoi tessellation allows us to extend the phase diagram toward the critical point.     

Thermal properties of ionic systems near the liquid-liquid critical point

Isobaric heat capacity per unit volume, C(p), and excess molar enthalpy, h(E), were determined in the vicinity of the critical point for a set of binary systems formed by an ionic liquid and a molecular solvent. Moreover, and, since critical composition had to be accurately determined, liquid-liquid equilibrium curves were also obtained using a calorimetric method. The systems were selected with a view on representing, near room temperature, examples from clearly solvophobic to clearly coulombic behavior, which traditionally was related with the electric permittivity of the solvent. The chosen molecular compounds are: ethanol, 1-butanol, 1-hexanol, 1,3-dichloropropane, and diethylcarbonate, whereas ionic liquids are formed by imidazolium-based cations and tetrafluoroborate or bis-(trifluromethylsulfonyl)amide anions. The results reveal that solvophobic critical behavior-systems with molecular solvents of high dielectric permittivity-is very similar to that found for molecular binary systems. However, coulombic systems-those with low permittivity molecular solvents-show strong deviations from the results usually found for these magnitudes near the liquid-liquid phase transition. They present an extremely small critical anomaly in C(p)-several orders of magnitude lower than those typically obtained for binary mixtures-and extremely low h(E)-for one system even negative, fact not observed, up to date, for any liquid-liquid transition in the nearness of an upper critical solution temperature.     

Behavior of a polymer chain immersed in a binary mixture of solvents

The behavior of a polymer chain immersed in a binary solvent mixture is investigated via a single-polymer simulation using an effective Hamiltonian, where the solvent effects are taken into account through a density-functional theory for polymer-solvent admixtures. The liquid-liquid phase separation of the binary solvent mixture is modeled as that of a Lennard-Jones binary fluid mixture with weakly attractive interactions between the different components. Two types of energetic preferences of the polymer chain for the better solvent-(A) no preferential solvophilicity and (B) strong preferential solvophilicity-are employed as polymer-solvent interaction models. The radius of gyration and the polymer-solvent radial distribution functions are determined from the simulations of various molar fractions along an isotherm slightly above the critical temperature of the liquid-liquid phase separation. These quantities near the critical point conspicuously depend on the strength of the preferential solvophilicity. In the case where the polymer exhibits a strong preferential solvophilicity, a remarkable expansion of the polymer chain is observed near the critical point. On the other hand, in the case where the polymer has no preferential solvophilicity, no characteristic variation of the polymer conformation is observed even near the critical point. These results indicate that the expansion of a polymer chain enhances the local phase separation around it, acting as a nucleus of demixing in the vicinity of the critical point. This phenomenon in binary solvents near the liquid-liquid critical point is similar to the expansion of the polymer chain in one-component supercritical solvents near the liquid-vapor critical point, which we have reported [T. Sumi and H. Sekino J. Chem. Phys. 122, 194910 (2005)].     

Polymer-induced phase separation and crystallization in immunoglobulin G solutions

We study the effects of the size of polymer additives and ionic strength on the phase behavior of a nonglobular protein-immunoglobulin G (IgG)-by using a simple four-site model to mimic the shape of IgG. The interaction potential between the protein molecules consists of a Derjaguin-Landau-Verwey-Overbeek-type colloidal potential and an Asakura-Oosawa depletion potential arising from the addition of polymer. Liquid-liquid equilibria and fluid-solid equilibria are calculated by using the Gibbs ensemble Monte Carlo technique and the Gibbs-Duhem integration (GDI) method, respectively. Absolute Helmholtz energy is also calculated to get an initial coexisting point as required by GDI. The results reveal a nonmonotonic dependence of the critical polymer concentration rho(PEG) (*) (i.e., the minimum polymer concentration needed to induce liquid-liquid phase separation) on the polymer-to-protein size ratio q (equivalently, the range of the polymer-induced depletion interaction potential). We have developed a simple equation for estimating the minimum amount of polymer needed to induce the liquid-liquid phase separation and show that rho(PEG) (*) approximately [q(1+q)(3)]. The results also show that the liquid-liquid phase separation is metastable for low-molecular weight polymers (q=0.2) but stable at large molecular weights (q=1.0), thereby indicating that small sizes of polymer are required for protein crystallization. The simulation results provide practical guidelines for the selection of polymer size and ionic strength for protein phase separation and crystallization.     

Crystalline-amorphous interaction in relation to the phase diagrams of binary polymer blends containing a crystalline constituent

The present article describes an equilibrium theory for determining binary phase diagrams of various crystalline-amorphous polymer blends by taking into account the contributions from both liquid-liquid phase separation between the constituents and solid-liquid phase transition of the crystalline component. An analytical expression for determining a crystal-amorphous interaction parameter is deduced based on the solid-liquid transition, involving the solidus and liquidus lines in conjunction with the coexistence curve of an upper critical solution temperature type. Of particular importance is that the crystalline-amorphous interaction parameter can be determined directly from the melting point depression data. The present analysis is therefore different from the conventional Flory-Huggins interaction parameter, which is associated with the liquid-liquid phase separation. The validity of the present theory is tested with the experimental phase diagrams of blends of poly(ethylene oxide)/diacrylate and poly(vinyl alcohol)/cellulose.     

Corresponding-states laws for protein solutions

The solvent around protein molecules in solutions is structured and this structuring introduces a repulsion in the intermolecular interaction potential at intermediate separations. We use Monte Carlo simulations with isotropic, pair-additive systems interacting with such potentials. We test if the liquid-liquid and liquid-solid phase lines in model protein solutions can be predicted from universal curves and a pair of experimentally determined parameters, as done for atomic and colloid materials using several laws of corresponding states. As predictors, we test three properties at the critical point for liquid-liquid separation: temperature, as in the original van der Waals law, the second virial coefficient, and a modified second virial coefficient, all paired with the critical volume fraction. We find that the van der Waals law is best obeyed and appears more general than its original formulation: A single universal curve describes all tested nonconformal isotropic pair-additive systems. Published experimental data for the liquid-liquid equilibrium for several proteins at various conditions follow a single van der Waals curve. For the solid-liquid equilibrium, we find that no single system property serves as its predictor. We go beyond corresponding-states correlations and put forth semiempirical laws, which allow prediction of the critical temperature and volume fraction solely based on the range of attraction of the intermolecular interaction potential.     

Monte Carlo simulation of polarizable systems: Early rejection scheme for improving the performance of adiabatic nuclear and electronic sampling Monte Carlo simulations

An early rejection scheme for trial moves in adiabatic nuclear and electronic sampling Monte Carlo simulation (ANES-MC) of polarizable intermolecular potential models is presented. The proposed algorithm is based on Swendsen–Wang filter functions for prediction of success or failure of trial moves in Monte Carlo simulations. The goal was to reduce the amount of calculations involved in ANES-MC electronic moves, by foreseeing the success of an attempt before making those moves. The new method was employed in Gibbs ensemble Monte Carlo (GEMC) simulations of the polarizable simple point charge-fluctuating charge (SPC-FQ) model of water. The overall improvement in GEMC depends on the number of swap attempts (transfer molecules between phases) in one Monte Carlo cycle. The proposed method allows this number to increase, enhancing the chemical potential equalization. For a system with 300 SPC-FQ water molecules, for example, the fractions of early rejected transfers were about 0.9998 and 0.9994 at 373 and 423 K, respectively. This means that the transfer moves consume only a very small part of the overall computing effort, making GEMC almost equivalent to a simulation in the canonical ensemble.     

GEMC和GDI方法计算流体气液相平衡的比较

考察采用TraPPE联合原子和OPLS全原子力场两种分子力场, Gibbs系综蒙特卡罗(GEMC)方法和Gibbs-Duhem积分(GDI)方法计算流体气液相平衡的适用性、计算速度、计算精度等问题. 结果表明, 在采用全原子力场情况下, GDI方法比GEMC方法极大地节省了计算时间. 从计算结果来看, 两种方法各有适用范围, 在使用时可互为补充. 在给定力场的前提下, 两种方法所得到的液相密度、蒸发焓、临界温度和临界密度相差不大, 而当力场中的缺陷导致蒸发焓的计算不够准确时, 两种计算方法得到的气体的压力和密度明显不同,进而导致预测的临界压力也明显不同.     

Liquid-gas separation in colloidal electrolytes

The liquid-gas transition of an electroneutral mixture of oppositely charged colloids, studied by Monte Carlo simulations, is found in the low-temperature-low-density region. The critical temperature shows a nonmonotonous behavior as a function of the interaction range, kappa(-1), with a maximum at kappasigma approximately 10, implying an island of coexistence in the kappa-rho plane. The system is arranged in such a way that each particle is surrounded by shells of particles with alternating charge. In contrast with the electrolyte primitive model, both neutral and charged clusters are obtained in the vapor phase.     

On the long wavelength thermodynamic limit of a neutron diffraction experiment in the vicinity of a liquid-liquid critical point. Application to the concentration fluctuations in the Li-ND3 system

A quantitative comparison is made between the thermodynamics as obtained from the long wavelength limit of a small angle neutron scattering experiment in the vicinity of a liquid-liquid critical point for the Li-ND₃ system and the data obtained from vapour pressure measurements. The agreement is fair. It is shown how the comparison always implies an underlying model of the interacting species leading to the liquid-liquid phase separation.     

On chip steady liquid-gas phase separation for flexible generation of dissolved gas concentration gradient

In this study, steady liquid-gas phase separation is realized by applying a hydrophobic small microchannel array (SMA) to bridge two large microchannels, one for liquid phase and one for gas phase. In this structure, a capillary pressure difference between that in the SMA and the larger channel results in a steady liquid-gas interface. The generated liquid-gas interface allows for fast gas dissolving speed. By coupling the liquid-gas interface with a one directional fluidic field, a steady dissolved gas concentration gradient (DgCG) is generated. The DgCG distribution is easily designable for linear or exponential modes, providing improved flexibility for gas participated processes on chip. To demonstrate its applicability, a CO(2) DgCG chip is fabricated and applied for screening CaCO(3) crystal growth conditions in the DgCG chip. Crystals with transitional structures are successfully fabricated, which is consistent with the CO(2) DgCG distribution.     

Thermodynamic, dynamic, structural, and excess entropy anomalies for core-softened potentials

Using molecular dynamic simulations, we study three families of continuous core-softened potentials consisting of two length scales: a shoulder scale and an attractive scale. All the families have the same slope between the two length scales but exhibit different potential energy gap between them. For each family three shoulder depths are analyzed. We show that all these systems exhibit a liquid-liquid phase transition between a high density liquid phase and a low density liquid phase ending at a critical point. The critical temperature is the same for all cases suggesting that the critical temperature is only dependent on the slope between the two scales. The critical pressure decreases with the decrease of the potential energy gap between the two scales suggesting that the pressure is responsible for forming the high density liquid. We also show, using the radial distribution function and the excess entropy analysis, that the density, the diffusion, and the structural anomalies are present if particles move from the attractive scale to the shoulder scale with the increase of the temperature indicating that the anomalous behavior depends only in what happens up to the second coordination shell.