No abstract

The self-consistent Ornstein-Zernike approximation (SCOZA), the generalized mean spherical approximation (GMSA), the modified hypernetted chain (MHNC) approximation, and the hierarchical reference theory (HRT) are applied to the determination of thermodynamic and structural properties, and the phase diagram of the hard-core Yukawa fluid (HCYF). We investigate different Yukawa-tail screening lengths lambda, ranging from lambda=1.8 (a value appropriate to approximate the shape of the Lennard-Jones potential) to lambda=9 (suitable for a simple one-body modelization of complex fluids like colloidal suspensions and globular protein solutions). The comparison of the results obtained with computer simulation data shows that at relatively low lambda's all the theories are fairly accurate in the prediction of thermodynamic and structural properties; as far as the phase diagram is concerned, the SCOZA and HRT are able to predict the binodal line and the critical parameters in a quantitative manner. At lambda=4 some discrepancies begin to emerge in the performances of the different theoretical approaches: the MHNC remains, on the whole, reasonably accurate in predicting the energy and the contact value of the radial distribution function; the SCOZA predicts well the equation of state up to the highest lambda values investigated. The GMSA and the MHNC underestimate and overestimate, respectively, the liquid coexisting density, while the SCOZA and HRT yield liquid branches that fall between the two former theoretical predictions, although both appear to overestimate the critical temperature somewhat. At higher lambda's the GMSA and MHNC binodals further worsen, while the SCOZA appears to remain usefully predictive. In general, the predictions of all the theories tend to slightly worsen at low temperatures and high density. The determination of the freezing line, performed by means of a one-phase "freezing criterion" (due to other authors) is not particularly satisfactory within either the SCOZA or the MHNC; the GMSA prediction for the freezing line at lambda=7 and 9 is instead able to follow in a qualitative manner the pattern of the solid-vapor coexistence line as determined through computer simulation studies. The necessity of further assessments of the freezing predictions is also discussed. Finally, versions of the GMSA, SCOZA, and HRT that can be expected to be more accurate for interactions with extremely short-ranged attractions are identified.

An integral-equation approach combined with molecular dynamics simulations based on the Girifalco spherical intermolecular potential has been used to predict the phase diagram for rigid C60 molecules. The boundary of the liquid-vapor coexistence region and the location of freezing and melting lines have been sketched out. The liquid phase is only observed in a very narrow temperature range compared with atomic systems (e.g., the rare gases). Unfortunately, the dense fluid is predicted to exist above -1800 K, which is sufficiently high that the C6o molecule may be unstable.

The free energy of the solid and fluid phases of the Girifalco C60 model are determined through extensive Monte Carlo simulations. In this model the molecules interact through a spherical pair potential, characterized by a narrow and attractive well, adjacent to a harshly repulsive core. We have used the Widom test particle method and a mapping from an Einstein crystal, in order to estimate the absolute free energy in the fluid and solid phases, respectively; we have then determined the free energy along several isotherms, and the whole phase diagram, by means of standard thermodynamic integrations. The dependence of the simulation's results on the size of the sample is also monitored in a number of cases.We highlight how the interplay between the liquid-vapor and the liquid-solid coexistence conditions determines the existence of a narrow liquid pocket in the phase diagram, whose stability is assessed and confirmed in agreement with previous studies. In particular, the critical temperature follows closely an extended corresponding-states rule recently outlined by Noro and Frenkel [J. Chem. Phys. 113, 2941 (2000)].We discuss the emerging "energetic" properties of the system, which drive the phase behavior in systems interacting through short-range forces [A. A. Louis, Phil. Trans. R. Soc. A 359, 939 (2001)], in order to explain the discrepancy between the predictions of several structural indicators and the results of full free energy calculations, to locate the fluid phase boundaries.More generally, we aim to provide extended reference data for calculations of the free energy of the C60 fullerite in the low temperature regime, as for the determination of the phase diagram of higher order Cn>60 fullerenes and other fullerene-related materials, whose description is based on the same model adopted in this work.

Articles you may be interested inPhase diagrams of model C 60 and C 70 fullerenes from short-range attractive potentialsa)Computer simulation study of the global phase behavior of linear rigid Lennard-Jones chain molecules: Comparison with flexible models

Recent studies have confirmed that model C60 has a stable liquid phase confined to a rather narrow temperature range; in this letter we determine the phase diagram of model C70 and show that also this fullerene possesses a liquid phase whose existence range is wider than that of C60. A "corresponding state" behaviour of C60 and C70, emerging from a description in terms of reduced quantities of both the interaction potential and the phase diagrams' parameters, is then extended to higher-order fullerenes. Critical-and triple-point densities and temperatures of C76, C84 and C92 are thereby predicted. The decreasing accuracy of such a representation of Cn>70 phase behaviour with increasing n is discussed in relation to the model potential features.

We numerically investigate colloidal dimers with asymmetric interaction strengths to study how the interplay between molecular geometry, excluded volume effects and attractive forces determines the overall phase behavior of such systems. Specifically, our model is constituted by two rigidly-connected tangent hard spheres interacting with other particles in the first instance via identical square-well attractions. Then, one of the square-well interactions is progressively weakened, until only the corresponding bare hard-core repulsion survives, giving rise to a "Janus dumbbell" model. We investigate structure, thermodynamics and phase behavior of the model by means of successive umbrella sampling and Monte Carlo simulations. In most of the cases, the system behaves as a standard simple fluid, characterized by a gas-liquid phase separation, for sufficiently low temperatures. In these conditions we observe a remarkable linear scaling of the critical temperature as a function of the interaction strength. But, as the interaction potential approaches the Janus dumbbell limit, we observe the spontaneous formation of self-assembled lamellar structures, preempting the gas-liquid phase separation. Comparison with previous studies allows us to pinpoint the role of the interaction range in controlling the onset of ordered structures and the competition between the formation of these structures and gas-liquid condensation.

The modified-hypernetted-chain theory is applied to the determination of the phase diagram of the Lennard-Jones (LJ) fiuid, and of a model of C«& previously investigated [Phys. Rev. Lett. 71, 1200 (1993)] through molecular-dynamics (MD) simulation and a different theoretical approach. In the LJ case the agreement with available MD data is quantitative and superior to other theories. For C60, the phase diagram obtained is in quite good agreement with previous MD results: in particular, the theory confirms the existence of a liquid phase between 1600 and 1920 K, the estimated triple point and critical temperature, respectively.