Phase diagram and surface tension of the hard-core attractive Yukawa model of variable range: Monte Carlo simulations

The liquid-vapor phase diagram and surface tension for hard-core Yukawa potential with 4相似文献   

Phase diagram of a binary symmetric hard-core Yukawa mixture

We assess the accuracy of the self-consistent Ornstein-Zernike approximation for a binary symmetric hard-core Yukawa mixture by comparison with Monte Carlo simulations of the phase diagrams obtained for different choices of the ratio alpha of the unlike-to-like interactions. In particular, from the results obtained at alpha=0.75 we find evidence for a critical endpoint in contrast to recent studies based on integral equation and hierarchical reference theories. The variation of the phase diagrams with range of the Yukawa potential is investigated.     

Application of the optimized Baxter model to the hard-core attractive Yukawa system

We perform Monte Carlo simulations on the hard-core attractive Yukawa system to test the optimized Baxter model that was introduced by Prinsen and Odijk [J. Chem. Phys. 121, 6525 (2004)] to study a fluid phase of spherical particles interacting through a short-range pair potential. We compare the chemical potentials and pressures from the simulations with analytical predictions from the optimized Baxter model. We show that the model is accurate to within 10% over a range of volume fractions from 0.1 to 0.4, interaction strengths up to three times the thermal energy, and interaction ranges from 6% to 20% of the particle diameter, and performs even better in most cases. We furthermore establish the consistency of the model by showing that the thermodynamic properties of the Yukawa fluid computed via simulations may be understood on the basis of one similarity variable, the stickiness parameter defined within the optimized Baxter model. Finally, we show that the optimized Baxter model works significantly better than an often used, naive method determining the stickiness parameter by equating the respective second virial coefficients based on the attractive Yukawa and Baxter potentials.     

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.     

Sedimentation equilibrium of colloidal suspensions in a planar pore based on density functional theory and the hard-core attractive Yukawa model

The sedimentation equilibrium of colloidal suspensions modeled by hard-core attractive Yukawa (HCAY) fluids in a planar pore is studied. The density profile of the HCAY fluid in a gravitational field and its distribution between the pore and uniform phases are investigated by a density functional theory (DFT) approach, which results from employing a recently proposed parameter-free version of the Lagrangian theorem-based density functional approximation (Zhou, S. Phys. Lett. A 2003, 319, 279) for hard-sphere fluids to the hard-core part of the HCAY fluid, and the second-order functional perturbation expansion approximation to the tail part as was done in a recent partitioned density functional approximation (Zhou, S. Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top. 2003, 68, 061201). The resultant DFT approach is, thus, the first adjustable parameter-free DFT for HCAY fluids. The validity of the present DFT for HCAY fluids of reduced range parameter z(red) = 1.8 under various external potentials is established in the first of the papers cited previously. The present DFT for HCAY fluids can predict the radial distribution function for the bulk HCAY fluid accurately in the colloidal limit (large value of z(red)), and in the hard-sphere limit, its prediction for the density profile of the hard-sphere fluid in a gravitational field is in very good agreement with the existing simulation data. The dependence of the density profile and distribution coefficient on the magnitude of the interparticle attraction, gravitational field, and degree of confinement is investigated in detail by the present DFT approach. Intuitive and qualitative analyses are also compared with the quantitative DFT calculational results.     

Interfacial and coexistence properties of soft spheres with a short-range attractive Yukawa fluid: molecular dynamics simulations

Molecular dynamics simulations have been carried out to obtain the interfacial and coexistence properties of soft-sphere attractive Yukawa (SAY) fluids with short attraction range, κ = 10, 9, 8, 7, 6, and 5. All our simulation results are new. These data are also compared with the recently reported results in the literature of hard-core attractive Yukawa (HAY) fluids. We show that the interfacial and coexistence properties of both potentials are different. For the surveyed systems, here we show that all coexistence curves collapse into a master curve when we rescale with their respective critical points and the surface tension curves form a single master curve when we plot γ* vs. T/T(c).     

Polymorph selection in the crystallization of hard-core Yukawa system

Colloid-colloid interactions in charge-stabilized dispersions can to some extent be represented by the hard-core Yukawa model. The crystallization process and polymorph selection of hard-core Yukawa model are studied by means of smart Monte Carlo simulations in the region of face-centered-cubic (fcc) phase. The contact value of hard-core Yukawa potential and the volume fraction of the colloids are fixed, while the Debye screening length can be varied. In the early stage of the crystallization, the precursors with relatively ordered liquid structure have been observed. Although the crystal structure of thermodynamically stable phase is fcc, the system crystallizes into a mixture of fcc and hexagonal close-packed (hcp) structures under small Debye screening length since the colloidal particles act as effective hard spheres. In the intermediate range of Debye screening length, the system crystallizes into a mixture of fcc, hcp, and body-centered-cubic (bcc). The existence of metastable hcp and bcc structures can be interpreted as a manifestation of the Ostwald’s step rule. Until the Debye screening length is large enough, the crystal structure obtained is almost a complete fcc suggesting the system eventually reaches to a thermodynamically stable state.     

Structure and thermodynamics of hard-core Yukawa fluids: thermodynamic perturbation approaches

The thermodynamic perturbation theories, which are based on the power series of a coupling constant (λ-expansion), have been proposed for studying the structural and thermodynamic properties of a hard-core Yukawa (HCY) fluid: one (A1-approximation) is the perturbation theory based on the hard-sphere repulsion as a reference system. The other (A2-approximation) is the perturbation theory based on the reference system which incorporates both the repulsive and short-range attractive interactions. The first-order mean-spherical approximation (FMSA) provided by Tang and Lu [J. Chem. Phys. 99, 9828 (1993)] has been employed for investigating the thermodynamic properties of a HCY fluid using the alternative method via the direct correlation function. The calculated results show that (i) the A1 and A2 approximations are in excellent agreements with previous computer simulation results in the literature and compare with the semi-empirical works of Shukla including the higher-order free energy terms, (ii) the A1 and A2 approximations are better than the FMSA and the mean-spherical approximation, (iii) the A2-approximation compares with the A1-approximation, even though the perturbation effect of an A2-approximation is much smaller than that of an A1-approximation, and that (iv) the FMSA study is particularly of advantage in providing the structure and thermodynamics in a simple and analytic manner.     

Close-packed structures and phase diagram of soft spheres in cylindrical pores

It is shown for a model system consisting of spherical particles confined in cylindrical pores that the first ten close-packed phases are in one-to-one correspondence with the first ten ways of folding a triangular lattice, each being characterized by a roll-up vector like the single-walled carbon nanotube. Phase diagrams in pressure-diameter and temperature-diameter planes are obtained by inherent-structure calculation and molecular dynamics simulation. The phase boundaries dividing two adjacent phases are infinitely sharp in the low-temperature limit but are blurred as temperature is increased. Existence of such phase boundaries explains rich, diameter-sensitive phase behavior unique for cylindrically confined systems.     

Thermodynamic and structural properties of repulsive hard-core Yukawa fluid: integral equation theory, perturbation theory and Monte Carlo simulations

The thermodynamic and structural properties of purely repulsive hard-core Yukawa particles in the fluid state are determined through Monte Carlo simulation and modeled using perturbation theory and integral equation theory in the mean spherical approximation (MSA). Systems of particles with Yukawa screening lengths of 1.8, 3.0, and 5.0 are examined with results compared to variations of MSA and perturbation theory. Thermodynamic properties were predicted well by both theories in the fluid region up to the fluid-solid phase boundary. Further, we found that a simplified exponential version of the MSA is the most accurate at predicting radial distribution function at contact. Radial distribution function of repulsive hard-core Yukawa particles are also reported. The results show that methods based on MSA and perturbation theory that are typically applied to the attractive hard-core Yukawa potential can also be extended to the purely repulsive hard-core Yukawa potential.     

Structures and adsorption of binary hard-core Yukawa mixtures in a slitlike pore: grand canonical Monte Carlo simulation and density-functional study

The grand canonical ensemble Monte Carlo simulation and density-functional theory are applied to calculate the structures, local mole fractions, and adsorption isotherms of binary hard-core Yukawa mixtures in a slitlike pore as well as the radial distribution functions of bulk mixtures. The excess Helmholtz energy functional is a combination of the modified fundamental measure theory of Yu and Wu [J. Chem. Phys. 117, 10156 (2002)] for the hard-core contribution and a corrected mean-field theory for the attractive contribution. A comparison of the theoretical results with the results from the Monte Carlo simulations shows that the corrected theory improves the density profiles of binary hard-core Yukawa mixtures in the vicinity of contact over the original mean-field theory. Both the present corrected theory and the simulations suggest that depletion and desorption occur at low temperature, and the local segregation can be observed in most cases. For binary mixtures in the hard slitlike pore, the present corrected theory predicts more accurate surface excesses than the original one does, while in the case of the attractive pore, no improvement is found in the prediction of a surface excess of the smaller molecule.     

Freezing lines of colloidal Yukawa spheres. I. A Rogers-Young integral equation study

Using the Rogers-Young (RY) integral equation scheme for the static structure factor combined with the one-phase Hansen-Verlet (HV) freezing rule, we study the equilibrium structure and two-parameter freezing lines of colloidal particles with Yukawa-type pair interactions representing charge-stabilized silica spheres suspended in dimethylformamide (DMF). Results are presented for a vast range of concentrations, salinities and effective charges covering particles with masked excluded-volume interactions. The freezing lines were obtained for the low-charge and high-charge solutions of the static structure factor, for various two-parameter sets of experimentally accessible system parameters. All RY-HV based freezing lines can be mapped on a universal fluid-solid coexistence line in good agreement with computer simulation predictions. The RY-HV calculations extend the freezing lines obtained in earlier simulations to a broader parameter range. The experimentally observed fluid-bcc-fluid reentrant transition of charged silica spheres in DMF can be explained using the freezing lines obtained in this work.     

Structure of inhomogeneous attractive and repulsive hard-core yukawa fluid: grand canonical Monte Carlo simulation and density functional theory study

A density functional theory is proposed for an inhomogeneous hard-core Yukawa (HCY) fluid based on Rosenfeld's perturbative method. The excess Helmholtz energy functional is derived from a modified fundamental measure theory for the hard-core repulsion and a quadratic functional Taylor expansion for the long-ranged attractive or repulsive interactions. To test the established theory, grand canonical ensemble Monte Carlo simulations are carried out to simulate the density profiles of attractive and repulsive HCY fluid near a wall. Comparison with the results from the Monte Carlo simulations shows that the present density functional theory gives accurate density profiles for both attractive and repulsive HCY fluid near a wall. Both the present theory and simulations suggest that there is depletion for attractive HCY fluid at low temperature, but no depletion is found for repulsive HCY fluid. The calculated results indicate that the present density functional theory is better than those of the modified version of the Lovett-Mou-Buff-Wertheim and other density functional theories. The present theory is simple in form and computationally efficient. It predicts accurate radial distribution functions of both attractive and repulsive HCY fluid except for the repulsive case at high density, where the theory overestimates the radial distribution function in the vicinity of contact.     

The low-pressure phase diagram of sulfur

The experimental values of the (p, T) phase diagram of sulfur have been assessed through empirical equations commonly used. This treatment shows insufficiencies, namely those leading to significant uncertainties in the coordinates of the low-pressure triple points. An alternative way to try and find reliable equations for the phase equilibrium curves is provided by the exact integration of Clapeyron equation established and developed in a previous work. No fitting to (p, T) experimental data is then necessary. According to this method equations have been found to all the sulfur low-pressure equilibrium lines, and at the same time providing the enthalpies of transition between the phases as functions of temperature, and reliable values for the coordinates of the low-pressure triple points. The details of all these calculations and of the corresponding results are the object of the present paper. The (p, T) low-pressure phase diagram of sulfur is reconstructed on the basis of the results of the calculations carried out following these lines.     

Structural-mechanical phase diagram of isotactic polypropylene

The polymorphic transformations associated with the plastic deformation of isotactic polypropylene samples of different stereoregularity, prepared with different metallocene catalysts, have been studied. Crystals of alpha or gamma forms, present in the unstretched samples, transform into the mesomorphic form by stretching. The formation of the mesophase facilitates successive further deformation up to very high strains and produces development of outstanding unusual properties of high flexibility and elasticity. The phase diagram of isotactic polypropylene, where the regions of stability of the different polymorphic forms are defined as a function of stereoregularity and degree of deformation, is reported.     

The phase diagram of ammonium nitrate

The pressure-temperature (P-T) phase diagram of ammonium nitrate (AN) [NH(4)NO(3)] has been determined using synchrotron x-ray diffraction (XRD) and Raman spectroscopy measurements. Phase boundaries were established by characterizing phase transitions to the high temperature polymorphs during multiple P-T measurements using both XRD and Raman spectroscopy measurements. At room temperature, the ambient pressure orthorhombic (Pmmn) AN-IV phase was stable up to 45 GPa and no phase transitions were observed. AN-IV phase was also observed to be stable in a large P-T phase space. The phase boundaries are steep with a small phase stability regime for high temperature phases. A P-V-T equation of state based on a high temperature Birch-Murnaghan formalism was obtained by simultaneously fitting the P-V isotherms at 298, 325, 446, and 467 K, thermal expansion data at 1 bar, and volumes from P-T ramping experiments. Anomalous thermal expansion behavior of AN was observed at high pressure with a modest negative thermal expansion in the 3-11 GPa range for temperatures up to 467 K. The role of vibrational anharmonicity in this anomalous thermal expansion behavior has been established using high P-T Raman spectroscopy.