Clustering phenomena in adhesive hard sphere fluids. Theory of freezing

We develop a version of the physical cluster theory of the equation of state based on a definition of a “cluster” which depends both on particle positions and relative momenta, namely we choose to call a group of i molecules a physical cluster if their total energy (kinetic energy plus interactions) does not exceed the translation energy of the centre-of-mass of the whole group. This definition in addition with the “microcrystal” model of a cluster allows us to calculate the properties of individual clusters, to consider the clustering phenomena and to write down an equation of state for a dilute system as well as for a system at moderate and high densities. The interaction between molecular aggregates is taken into consideration by using the adhesive hard sphere model and the Percus-Yevick approximation. Based on the relations obtained we develop the theory of freezing and perform a numerical investigation of the equilibrium behavior of the system at various temperatures and densities. In particular, the theory enables us to describe the regular phase diagram which includes (in coordinates pressure-temperature) three branches intersecting at the triple point. It is interesting that the coexistence curve of the fluid-solid as well as liquid-gas phase transition has an end point.     

The evaporation/condensation transition of liquid droplets

The condensation of a supersaturated vapor enclosed in a finite system is considered. A phenomenological analysis reveals that the vapor is found to be stable at densities well above coexistence. The system size at which the supersaturated vapor condenses into a droplet is found to be governed by a typical length scale which depends on the coexistence densities, temperature and surface tension. When fluctuations are neglected, the chemical potential is seen to show a discontinuity at an effective spinodal point, where the inhomogeneous state becomes more stable than the homogeneous state. If fluctuations are taken into account, the transition is rounded, but the slope of the chemical potential versus density isotherm develops a discontinuity in the thermodynamic limit. In order to test the theoretical predictions, we perform a simulation study of droplet condensation for a Lennard-Jones fluid and obtain loops in the chemical potential versus density and pressure. By computing probability distributions for the cluster size, chemical potential, and internal energy, we confirm that the effective spinodal point may be identified with the occurrence of a first order phase transition, resulting in the condensation of a droplet. An accurate equation of state is employed in order to estimate the droplet size and the coexisting vapor density and good quantitative agreement with the simulation data is obtained. The results highlight the need of an accurate equation of state data for the Laplace equation to have predictive power.     

Phase behavior of attractive and repulsive ramp fluids: integral equation and computer simulation studies

Using computer simulations and a thermodynamically self-consistent integral equation we investigate the phase behavior and thermodynamic anomalies of a fluid composed of spherical particles interacting via a two-scale ramp potential (a hard core plus a repulsive and an attractive ramp) and the corresponding purely repulsive model. Both simulation and integral equation results predict a liquid-liquid demixing when attractive forces are present, in addition to a gas-liquid transition. Furthermore, a fluid-solid transition emerges in the neighborhood of the liquid-liquid transition region, leading to a phase diagram with a somewhat complicated topology. This solidification at moderate densities is also present in the repulsive ramp fluid, but in this case inhibits the fluid-fluid separation.     

球腔中氢键流体的吸附-解吸附转变:表面调控研究

集中讨论了球形微腔表面对腔中氢键流体相态结构的调控机制. 为了揭示微腔表面对腔中氢键流体相平衡的影响, 首先根据吸附-解吸附原理并利用经典流体的密度泛函理论计算了微腔中氢键流体的平衡密度分布, 进而通过吸附-解吸附等温线及巨势等温线绘制出体系的相图. 在此基础上, 重点考察了球腔尺寸、 表面作用强度和作用力程对氢键流体毛细凝聚及层化转变的影响. 结果表明, 这些因素可以有效地调控体系毛细凝聚和层化转变的临界约化温度、 临界密度和相区大小等特征, 从而阐明了表面调控的主要机制. 研究结果为设计相关吸附材料提供了理论参考.     

Sodium chloride in supercritical water as a function of density: potentials of mean force and an equation for the dissociation constant from 723 to 1073 K and from 0 to 0.9 g/cm(3)

The potential of mean force (PMF) of sodium chloride in water has been calculated by using the ab initio classical free-energy perturbation method at five state points: at 973 K with densities of 0.2796, 0.0935, and 0.0101 g/cm (3) and at 723 K with densities of 0.0897 and 0.0098 g/cm (3). The method is based on a QM-MM model in which Na-H 2O, Cl-H 2O, and Na-Cl interactions are calculated by ab initio methods. The water-water interactions are from the polarizable TIP4P-FQ model. The logarithm of the dissociation constant (log K c) has been calculated from the PMF. These predictions, together with experimental measurements, were used to derive an equation for log K c at densities from 0 to 0.9 g/cm (3) and temperatures from 723 to 1073 K, as well as from 600 to 1073 K for densities from 0.29 g/cm (3) to 0.9 g/cm (3). Extrapolation of the present equation below 723 K for densities less than 0.29 g/cm (3) does not fit the experimental results. This is attributed to long-range changes in the local dielectric constant due to the high compressibility. Comparisons with previous predictions and simulations are presented.     

Phase behavior of ionic fluids in slitlike pores: a density functional approach for the restricted primitive model

We present a density functional theory of nonuniform ionic fluids. This theory is based on the application of the electrostatic contribution to the free energy functional arising from mean spherical approximation for a bulk restricted primitive model and from the energy route bulk equation of state. In order to employ this functional we define a reference fluid and additional averaged densities, according to the approach introduced by Gillespie, Nonner and Eisenberg [J. Phys.: Condens. Matter 14, 12129 (2002)]. In the case of bulk systems the proposed theory reduces to the mean spherical approximation equation of state, arising from the energy route and thus it predicts the first-order phase transition. We use this theory to investigate the effects of confinement on the liquid-vapor equilibria. Two cases are considered, namely an electrolyte confined to the pore with uncharged walls and with charged walls. The dependence of the capillary evaporation diagrams on the pore width and on the electrostatic potential is determined.     

Equation of state for hard-spheres and excess function calculations of binary liquid mixtures

Using our previously proposed matrix method, an equation of state for hard spheres is presented, which can reproduce the exact values of the first-eight virial coefficients. This equation meets both the low density and the close-packed limits and can predicts the first order fluid–solid phase transition of hard spheres. The results obtained show that the new equation of state can correlate the simulation data of compressibility factor up to high densities better than other equations of state.The new equation of state is extended to mixtures of hard spheres and excess functions of various binary liquid mixtures are calculated using the perturbation theory of Leonard–Henderson–Barker. The results are compared with existing theoretical and experimental data and with those calculated by other hard-sphere equations of state.It is seen that the results obtained by the new equation of state is quite satisfactory compared to other equations of state for the hard spheres and mixture of hard spheres.     

Gas-liquid Transition in Charged Fluids

Abstract

The connexion between the equation of state of a classical fluid of non-polarizable ions, the character of the screening, and the appearance of long range oscillations in the chargetharge radial distribution function is examined. While considerations of stability lead to the usual inequalities for the inverse static dielectric function and the compressibility of the charged fluid, the square of the inverse screening length k ₃, does not need to be positive for thermodynamic stability. Through a study of an approximate equation of state for a twocomponent fluid of charged hard spheres, the regions of negative and positive k ² in the pressuredensity plane are related to a liquid phase and to an ionized-gas phase, respectively. The model fluid displays a gas-liquid critical point, above which the transition between the two types of screening is continuous. This behaviour of the charged-hard-spheres fluid is contrasted with the transition of a real ionic liquid to the molecular gaseous phase.     

Phase behavior of polymer/nanoparticle blends near a substrate

We use the recent fluids density functional theory of Tripathi and Chapman [Phys. Rev. Lett. 94, 087801 (2005); J. Chem. Phys. 122, 094506 (2005)] to investigate the phase behavior of athermal polymer/nanoparticle blends near a substrate. The blends are modeled as a mixture of hard spheres and freely jointed hard chains, near a hard wall. There is a first order phase transition present in these blends in which the nanoparticles expel the polymer from the surface to form a monolayer at a certain nanoparticle concentration. The nanoparticle transition density depends on the length of the polymer, the nanoparticle diameter, and the overall bulk density of the system. The phase transition is due to both packing entropy effects related to size asymmetry between the components and to the polymer configurational entropy, justifying the so-called "entropic push" observed in experiments. In addition, a layered state is found at higher densities which resembles that in colloidal crystals, in which the polymer and nanoparticles form alternating discrete layers. We show that this laminar state has nearly the same free energy as the homogeneously mixed fluid in the bulk and is nucleated by the surface.     

Equilibrium limit of diffusion equation for interacting particles: are bifurcations and multiplicity of correlation functions an artifact or real?

The structure of a fluid is analyzed by taking the equilibrium limit of a diffusion equation including the Giacomin-Lebowitz term for intermolecular interactions. This equation represents the differential mass balance in fluids with the Metropolis algorithm for fluxes; it allows a new qualitative yet analytical approximation for the direct correlation function over the entire range of fluid densities and temperatures. This approximation is analogous to a classical Ono-Kondo model for adsorption if the distribution of molecules around a central molecule is viewed as the adsorption of molecules on a central molecule. While this model qualitatively predicts known behavior for both gas and liquid phases, approaching a phase transition (e.g., condensation of gas into liquid) results in a bifurcation and multiplicity of the direct correlation function. The model predicts a sequence in the transformation of correlation functions from that of a gas to that of a liquid. This sequence starts with the appearance of an isolated loop in the direct correlation function, indicating states that are stable but cannot be achieved without perturbation of the system. However, the system seems to sense its proximity to a phase transition and reflects the distance to the phase boundary by the size and shape of this isolated loop in the direct correlation function. At lower temperatures, this loop merges with the gas-phase peak, indicating that clusters can form spontaneously. Then, these clusters grow into a high-density (liquid) phase. The notion of bifurcations and multiplicity in correlation functions is an unusual and controversial concept. Certainly, it is unexpected and raises important questions: (a) if such a behavior is not real, why does the diffusion equation predict such behavior, i.e., is it a mathematical artifact or is it due to conflicting physical assumptions? (b) if this behavior is real, how does one interpret it at a molecular level? Here, we present some interpretations, but they are open for discussion.     

Computational study of the melting-freezing transition in the quantum hard-sphere system for intermediate densities. I. Thermodynamic results

The points where the fluid-solid (face-centered-cubic) transition takes place in the quantum hard-sphere system, for reduced densities 0.85>rhoN*>0.5 (reduced de Broglie wavelengths lambdaB*相似文献   

Strange Kinetics of the C(2)H(6) + CN Reaction Explained

In this paper, we employ state of the art quantum chemical and transition state theory methods in making a priori kinetic predictions for the abstraction reaction of CN with ethane. This reaction, which has been studied experimentally over an exceptionally broad range of temperature (25-1140 K), exhibits an unusually strong minimum in the rate constant near 200 K. The present theoretical predictions, which are based on a careful consideration of the two distinct transition state regimes, quantitatively reproduce the measured rate constant over the full range of temperature, with no adjustable parameters. At low temperatures, the rate-determining step for such radical-molecule reactions involves the formation of a weakly bound van der Waals complex. At higher temperatures, the passage over a subthreshold saddle point on the potential energy surface, related to the formation and dissolution of chemical bonds, becomes the rate-determining step. The calculations illustrate the changing importance of the two transition states with increasing temperature and also clearly demonstrate the need for including accurate treatments of both transition states. The present two transition state model is an extension of that employed in our previous work on the C2H4 + OH reaction [J. Phys. Chem. A 2005, 109, 6031]. It incorporates direct ab initio evaluations of the potential in classical phase space integral based calculations of the fully coupled anharmonic transition state partition functions for both transition states. Comparisons with more standard rigid-rotor harmonic oscillator representations for the "inner" transition state illustrate the importance of variational, anharmonic, and nonrigid effects. The effects of tunneling through the "inner" saddle point and of dynamical correlations between the two transition states are also discussed. A study of the kinetic isotope effect provides a further test for the present two transition state model.     

Phase Behaviors of Soft-core Particle Systems

This paper reviews some of our recent works on phase behaviors of particulate systems with a soft-core interaction potential.The potential is purely repulsive and bounded,i.e.,it is finite even when two particles completely overlap.The one-sided linear spring(harmonic) potential is one of the representatives.This model system has been successively employed to study the jamming transition,i.e.,the formation of rigid and disordered packings of hard particles,and establish the jamming physics.This is actually based on the 4"hard,:aspect of the potential,because at low densities and when particle overlap is tiny the potential resembles the hard sphere limit.At high densities,the potential exhibits its"soft"aspect:with the increase of density,there are successive reentrant crystallizations with many types of solid phases.Taking advantage of the dual nature of the potential,we investigate the criticality of the jamming transition from different perspectives,extend the jamming scenario to high densities,reveal the novel density evolution of two-dimensional melting,and find unexpected formation of quasicrystals.It is surprising that such a simple potential can exhibit so rich and unexpected phenomena in phase transitions.The phase behaviors discussed in this paper are also highly regarded in polymer science,which may thus shed light on our understanding of polymeric systems or inspire new ideas in studies of polymers.     

Critical dynamics of dimers: implications for the glass transition

The Adam-Gibbs view of the glass transition relates the relaxation time to the configurational entropy, which goes continuously to zero at the so-called Kauzmann temperature. We examine this scenario in the context of a dimer model with an entropy-vanishing phase transition and stochastic loop dynamics. We propose a coarse-grained master equation for the order parameter dynamics which is used to compute the time-dependent autocorrelation function and the associated relaxation time. Using a combination of exact results, scaling arguments, and numerical diagonalizations of the master equation, we find nonexponential relaxation and a Vogel-Fulcher divergence of the relaxation time in the vicinity of the phase transition. Since in the dimer model the entropy stays finite all the way to the phase transition point and then jumps discontinuously to zero, we demonstrate a clear departure from the Adam-Gibbs scenario. Dimer coverings are the "inherent structures" of the canonical frustrated system, the triangular Ising antiferromagnet. Therefore, our results provide a new scenario for the glass transition in supercooled liquids in terms of inherent structure dynamics.     

Orientational ordering in hard rectangles: The role of three-body correlations

We investigate the effect of three-body correlations on the phase behavior of hard rectangle two-dimensional fluids. The third virial coefficient B3 is incorporated via an equation of state that recovers scaled particle theory for parallel hard rectangles. This coefficient, a functional of the orientational distribution function, is calculated by Monte Carlo integration, using an accurate parametrized distribution function, for various particle aspect ratios in the range of 1-25. A bifurcation analysis of the free energy calculated from the obtained equation of state is applied to find the isotropic (I)-uniaxial nematic (N(u)) and isotropic-tetratic nematic (N(t)) spinodals and to study the order of these phase transitions. We find that the relative stability of the N(t) phase with respect to the isotropic phase is enhanced by the introduction of B3. Finally, we have calculated the complete phase diagram using a variational procedure and compared the results with those obtained from scaled particle theory and with Monte Carlo simulations carried out for hard rectangles with various aspect ratios. The predictions of our proposed equation of state as regards the transition densities between the isotropic and orientationally ordered phases for small aspect ratios are in fair agreement with simulations. Also, the critical aspect ratio below which the N(t) phase becomes stable is predicted to increase due to three-body correlations, although the corresponding value is underestimated with respect to simulation.     

Mutual diffusion coefficient in binary mixtures in different aggregation states

An equation for the mutual diffusion coefficient for the components of a binary mixture occurring in various aggregation states has been derived. Situations of this type arise in narrow pores in which the adsorbate density sharply changes across the pore cross-section. The equation was derived with the assumption that molecules of the mixture components have a spherical shape and are similar in size. The constructed equations are based on the lattice-gas model applicable to any aggregation state. The migration of particles is described in terms of the transition state theory. At low mixture densities corresponding to an ideal gas phase, the equations are based on expressions of the rigorous kinetic theory of gases. The theory takes into account the change in the mechanism of particle migration in different phases: from pair collisions for the gas to the overcoming of the activation barrier by thermofluctuation for dense phases. __________ Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 8, pp. 1717–1725, August, 2005.     

Correlated chain dynamics in the nematic phase of 4-cyano-4'-hexyloxybiphenyl

We report in this paper measurements of the Zeeman and Quadrupolar spin-lattice relaxation times at two different deuteron Larmor frequencies for the nematic phase of the perdeuteriated nematogen 6OCB. A model of correlated internal motion is used to account for both the quadrupolar splittings and the spectral densities of motion. The nematic mean field is constructed using the additive potential method, while the conformational transitions among the allowed configurations are described by a master equation. For modelling the quadrupolar splittings, we used 729 conformations, although this seemed impossible when the spectral densities were fitted by minimizing the sum of squares of per cent errors. The pentane effect has, therefore, been used to limit the size of the transition rate matrix in the master equation. We found that the dynamic model for liquid crystals proposed by one of us (1991, Phys. Rev. A, 43, 4310) is essentially correct for 6OCB.     

Equation of state and structural properties of the Weeks-Chandler-Andersen fluid

Molecular dynamics simulations have been carried out for the equation of state and percolation properties of the Weeks-Chandler-Andersen (WCA) system in its fluid phase as functions of density and temperature. The compressibility factor Z collapses well for the various isotherms, using an effective particle diameter for the WCA particle which is (in the usual WCA reduced units) sigma(e)=2(16)(1+T)(16), where T is the temperature. A corresponding "effective" packing fraction is zeta(e)=pisigma(e) (3)N6V, for N particles in volume V, which therefore scales out the effects of temperature. Using zeta(e) the simulation derived Z can be fitted to a simple analytic form which is similar to the Carnahan-Starling hard sphere equation of state and which is valid at all temperatures and densities where the WCA fluid is thermodynamically stable. The data, however, are not scalable onto the hard sphere equation of state for the complete packing fraction range. We explored the continuum percolation behavior of the WCA fluids. The percolation distance sigma(p) for the various states collapses well onto a single curve when plotted as sigma(p)sigma(e) against zeta(e). The ratio sigma(p)sigma(e) exhibits a monotonic decrease with increasing zeta(e) between the percolation line for permeable spheres and the glass transition limit, where sigma(p)sigma(e) approximately 1. The percolation packing fraction was calculated as a function of effective packing fraction and fitted to an empirical expression. The local coordination number at the percolation threshold showed a transition between the soft core and hard core limits from ca. 2:74 to 1:5, as previously demonstrated in the literature for true hard spheres. A number of simple analytic expressions that represent quite well the percolation characteristics of the WCA system are proposed.