The critical behavior of the refractive index near liquid-liquid critical points

The nature of the critical behavior in the refractive index n is revisited in the framework of the complete scaling formulation. A comparison is made with the critical behavior of n as derived from the Lorentz-Lorenz equation. Analogue anomalies to those predicted for the dielectric constant ε, namely, a leading |t|(2β) singularity in the coexistence-curve diameter in the two-phase region and a |t|(1-α) along the critical isopleth in the one phase region, are expected in both cases. However, significant differences as regards the amplitudes of both singularities are obtained from the two approaches. Analysis of some literature data along coexistence in the two-phase region and along the critical isopleth in the one-phase region provide evidence of an intrinsic effect, independent of the density, in the critical anomalies of n. This effect is governed by the shift of the critical temperature with an electric field, which is supposed to take smaller values at optical frequencies than at low frequencies in the Hz to MHz range.     

Phase behavior of model confined fluids. Influence of substrate-fluid interaction strength

We examine the relationship between the macroscopic phase behavior of nanoconfined fluids and the nature of microscopic interactions between a confining substrate and fluid. Two model slit-pore systems are explored using grand canonical transition-matrix Monte Carlo simulation. One system consists of a square-well fluid interacting with a square-well substrate, and the other contains an embedded point charge model of lysozyme interacting with a mica surface. Fluid phase diagrams are constructed for a broad range of substrate conditions. Our results indicate that one observes a maximum in the critical temperature of the fluid phase envelope upon variation of substrate strength for a given slit width. Both systems studied exhibit such maxima at intermediate wall strength. The physical rationale for this observation suggests that this behavior should be generally expected. We introduce two metrics that enable one to predict conditions that produce maxima in critical temperature. The first is related to the contact angle a fluid develops at a single confining substrate. The second is based upon virial coefficient information and requires knowledge of the substrate-fluid and fluid-fluid interaction potentials only.     

Application of a renormalization-group treatment to the statistical associating fluid theory for potentials of variable range (SAFT-VR)

An accurate prediction of phase behavior at conditions far and close to criticality cannot be accomplished by mean-field based theories that do not incorporate long-range density fluctuations. A treatment based on renormalization-group (RG) theory as developed by White and co-workers has proven to be very successful in improving the predictions of the critical region with different equations of state. The basis of the method is an iterative procedure to account for contributions to the free energy of density fluctuations of increasing wavelengths. The RG method has been combined with a number of versions of the statistical associating fluid theory (SAFT), by implementing White's earliest ideas with the improvements of Prausnitz and co-workers. Typically, this treatment involves two adjustable parameters: a cutoff wavelength L for density fluctuations and an average gradient of the wavelet function Φ. In this work, the SAFT-VR (variable range) equation of state is extended with a similar crossover treatment which, however, follows closely the most recent improvements introduced by White. The interpretation of White's latter developments allows us to establish a straightforward method which enables Φ to be evaluated; only the cutoff wavelength L then needs to be adjusted. The approach used here begins with an initial free energy incorporating only contributions from short-wavelength fluctuations, which are treated locally. The contribution from long-wavelength fluctuations is incorporated through an iterative procedure based on attractive interactions which incorporate the structure of the fluid following the ideas of perturbation theories and using a mapping that allows integration of the radial distribution function. Good agreement close and far from the critical region is obtained using a unique fitted parameter L that can be easily related to the range of the potential. In this way the thermodynamic properties of a square-well (SW) fluid are given by the same number of independent intermolecular model parameters as in the classical equation. Far from the critical region the approach provides the correct limiting behavior reducing to the classical equation (SAFT-VR). In the critical region the β critical exponent is calculated and is found to take values close to the universal value. In SAFT-VR the free energy of an associating chain fluid is obtained following the thermodynamic perturbation theory of Wertheim from the knowledge of the free energy and radial distribution function of a reference monomer fluid. By determining L for SW fluids of varying well width a unique equation of state is obtained for chain and associating systems without further adjustment of critical parameters. We use computer simulation data of the phase behavior of chain and associating SW fluids to test the accuracy of the new equation.     

A Thermodynamic Model for Square-well Chain Fluid: Theory and Monte Carlo Simulation

A thermodynamic model for the freely jointed square-well chain fluids was developed based on the thermodynamic perturbation theory of Barker-Henderson, Zhang and Wertheim. In this derivation Zhang's expressions for square-well monomers improved from Barker-Henderson compressibility approximation were adopted as the reference fluid, and Wertheim's polymerization method was used to obtain the free energy term due to the bond connectivity. An analytic expression for the Helmholtz free energy of the square-well chain fluids was obtained. The expression without adjustable parameters leads to the thermodynamic consistent predictions of the compressibility factors, residual internal energy and constant-volume heat capacity for dimer, 4-mer, 8-mer and 16-mer square-well fluids. The results are in good agreement with the Monte Carlo simulation. To obtain the MC data of residual internal energy and the constant-volume heat capacity needed, NVT MC simulations were performed for these square-well chain fluids.     

Phase behavior of colloids and proteins in aqueous suspensions: Theory and computer simulations

The fluid phase behavior of colloidal suspensions with short-range attractive interactions is studied by means of Monte Carlo computer simulations and two theoretical approximations, namely, the discrete perturbation theory and the so-called self-consistent Ornstein-Zernike approximation. The suspensions are modeled as hard-core attractive Yukawa (HCAY) and Asakura-Oosawa (AO) fluids. A detailed comparison of the liquid-vapor phase diagrams obtained through different routes is presented. We confirm Noro-Frenkel's extended law of scaling according to which the properties of a short-ranged fluid at a given temperature and density are independent of the detailed form of the interaction, but just depend on the value of the second virial coefficient. By mapping the HCAY and AO fluids onto an equivalent square-well fluid of appropriate range at the critical point we show that the critical temperature as a function of the effective range is independent of the interaction potential, i.e., all curves fall in a master curve. Our findings are corroborated with recent experimental data for lysozyme proteins.     

Determination of critical indices by "slow" spectroscopy: NMR shifts by statistical thermodynamics and density functional theory calculations

The temperature dependencies of NMR shifts in the critical region of two coexisting phases have been simulated using statistical thermodynamics and graph-theory consideration of equilibrium processes of molecular association. Microparameters of magnetic screening of various water and water/pyridine structures used in the statistical averaging have been evaluated by density functional theory calculations (PBE1PBE and B3PW91 functionals in the 6-311++G** basis set). The gauge-including atomic orbital (GIAO) approach has been applied to ensure gauge invariance of the results. Solvent effects were taken into account by a polarized continuum model (PCM). NMR shifts "order parameters" (Deltadelta = |delta+ - delta-|) and "diameters" (phidelta = |(delta+ + delta-)/2 - deltaC|, where delta+, delta-, and deltaC are the chemical shifts of coexisting phases and at the critical point respectively) have been calculated in each case close to the lower critical solution point (TL) and processed using linear regression analysis of Deltadelta approximately |T - TL| and phidelta approximately |T - TL| in the log-log plot. It has been shown that the critical index beta can be evaluated with high precision from the slope of Deltadelta = f(T - TL) at any realistic set of model input parameters. The slope of diameter has been found to depend on both input beta and alpha values. The obtained phidelta slopes (0.58-0.63) are very close to 2beta values. The results are discussed within the concept of complete scaling. Results of simulation are compared and supported by experimental NMR data for water/2,6-lutidine, acetic anhydride/n-heptane, and acetic anhydride/cyclohexane systems.     

Role of solvent for globular proteins in solution

The properties of the solvent affect the behavior of the solution. We propose a model that accounts for the contribution of the solvent free energy to the free energy of globular proteins in solution. For the case of an attractive square-well potential, we obtain an exact mapping of the phase diagram of this model without solvent to the model that includes the solute-solvent contribution. In particular we find for appropriate choices of parameters upper critical points, lower critical points, and even closed loops with both upper and lower critical points similar to those found before [Macromolecules 36, 5843 (2003)]. In the general case of systems whose interactions are not attractive square wells, this mapping procedure can be a first approximation to understand the phase diagram in the presence of solvent. We also present simulation results for both the square-well model and a modified Lennard-Jones model.     

Collective and self-diffusion coefficients in an ionic critical mixture: 3-methylpyridine+water+NaBr

The dynamics of concentration fluctuations of three critical samples of the 3-methylpyridine (3MP)+water+NaBr system have been measured by photon correlation spectroscopy. The collective-diffusion coefficient D shows the usual Ising behavior near the critical temperature T(c). However, as |T-T(c)| increases, the dynamic correlation length calculated from D, xi, takes values higher than the correlation length of the critical fluctuations calculated from static light scattering, xi(s). At the largest |T-T(c)| measured, xi approaches the value, xi(0,d) approximately equal to 1.13 nm, while the amplitude of xi(s) is xi(0,s)=0.38 nm. Pulsed-gradient NMR spectroscopy points out the existence of two dynamic contributions. One of them is consistent with the existence of molecular entities of hydrodynamic radius 0.31 nm, while the other one indicates the existence of aggregates rich in 3MP of radius 1.16 nm. The existence of the aggregates may explain the apparent anomalous behavior of the dynamic light scattering experiments for this system far from the critical point.     

Calculating thermodynamic properties from perturbation theory: II. An analytic representation for the square-well chain fluid

An analytic representation of thermodynamic properties of the freely jointed square-well chain fluid is developed based on the thermodynamic perturbation theory of Barker–Henderson, Zhang and Weitheim. By using a real function expression for the radial distribution function and incorporating structural information for square-well monomer of TPT1 model, an analytic expression for the Helmholtz energy of square-well chain fluid is expanded from Zhang’s analytic expressions for thermodynamic properties of square-well monomer. The expression leads to good predictions of the compressibility factor, residual internal energy and constant-volume heat capacity for 4-mer, 8-mer and 16-mer square-well fluids when compared with the Monte Carlo (MC) simulation results. The incorporating structural information for square-well dimer of TPT-D model is also calculated. To obtain the constant-volume heat capacity needed, NVT MC simulations were performed.     

Molecular association of heteronuclear vibrating square-well dumbbells in liquid-vapor phase equilibrium

Molecular aggregates are formed by heteronuclear vibrating square-well dumbbells. In a recent article [G. A. Chapela and J. Alejandre, J. Chem. Phys., 132(10), 104704 (2010)], it is shown that heteronuclear vibrating square-well dumbbells with a diameter ratio between particles of 1/2 and interacting potential ratio of 4 form micelles of different sizes and shapes which manifest themselves in both the liquid and vapor phases, up to and above the critical point. This means that micellization and phase separation are present simultaneously in this simple model. These systems present a maximum in the critical temperature when plotted against the potential well depth of the second particle ε(2). In the same publication, it was speculated that the formation of micelles was responsible for the appearance of the maximum. A thorough study on this phenomena is presented here and it is found that there is a threshold on the size of the second particle and its corresponding depth of interaction potential, where the micelles are formed. If the diameter and well depth of the second particle are small enough for the first and deep enough for the second, micelles are formed. For σ(2)/σ(1) between 0.25 and 0.65 and ε(2)/ε(1) larger than 5.7, micelles are formed up to and above the critical temperature. Outside these ranges micelles appear only at temperatures lower than the critical point. There is a strong temperature dependence on the formation and persistence of the aggregates. For the deepest wells and large enough second particles, a gel interconnected aggregate is obtained. In this work, the micelles are formed at temperatures as low as the triple point and as high as the critical point and, in some cases, persist well above it. The presence of these maxima in critical temperatures T(c) when plotted against ε(2) as follows. At lower values of ε(2), an increase of T(c) is obtained as is expected by the increase of the attractive volume as indicated by the principle of corresponding states. As ε(2) increases further, the formation of molecular aggregates produce a saturation effect of the deepening of the potential well by encapsulating the particles of the second kind inside the micelles, so the resulting T(c) represents a new poly disperse system of molecular aggregates and not the original heteronuclear vibrating square-well dumbbells. The surface tension is also analyzed for these systems, and it is shown that decreases with increasing attraction due to the formation of molecular aggregates.     

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.     

Parametric studies of interaction strengths in polymer/CO2 systems: discontinuous molecular dynamics simulations

Discontinuous molecular dynamics simulations are performed on homopolymer/solvent and surfactant/solvent systems. The homopolymer and surfactant molecules are modeled as freely jointed square-well chains. Solvent molecules are modeled as both hard spheres and square-well spheres. We explore how the various interaction parameters affect the types of phase behavior and micellization observed in the homopolymer/solvent and surfactant/solvent systems. Increasing the packing fraction of homopolymers in both hard-sphere solvents and square-well solvents increases the solvent's ability to dissolve homopolymers only when the segment-solvent interaction strength exceeds a critical value. Although only upper critical solution temperature (UCST) behavior is observed for homopolymers in hard-sphere solvents, both UCST and lower critical solution temperature (LCST) behavior are observed for homopolymers in square-well solvents, depending upon the interaction strengths and chain length. This indicates that it is necessary to account for the solvent-solvent attraction to model LCST behavior in supercritical CO2. Our simulation results on surfactants in hard-sphere solvents show that it is necessary to account for the interactions experienced by both the head and tail blocks in order to capture the essential features of surfactant/supercritical CO2 systems.     

Vapor-liquid equilibrium and critical behavior of the square-well fluid of variable range: a theoretical study

The vapor-liquid phase behavior and the critical behavior of the square-well (SW) fluid are investigated as a function of the interaction range, lambdain [1.25, 3], by means of the self-consistent Ornstein-Zernike approximation (SCOZA) and analytical equations of state based on a perturbation theory [A. L. Benavides and F. del Rio, Mol. Phys. 68, 983 (1989); A. Gil-Villegas, F. del Rio, and A. L. Benavides, Fluid Phase Equilib. 119, 97 (1996)]. For this purpose the SCOZA, which has been restricted up to now to a few model systems, has been generalized to hard-core systems with arbitrary interaction potentials requiring a fully numerical solution of an integro-partial differential equation. Both approaches, in general, describe well the liquid-vapor phase diagram of the square-well fluid when compared with simulation data. SCOZA yields very precise predictions for the coexistence curves in the case of long ranged SW interaction (lambda>1.5), and the perturbation theory is able to predict the binodal curves and the saturated pressures, for all interaction ranges considered if one stays away from the critical region. In all cases, the SCOZA gives very good predictions for the critical temperatures and the critical pressures, while the perturbation theory approach tends to slightly overestimate these quantities. Furthermore, we propose analytical expressions for the critical temperatures and pressures as a function of the square-well range.     

A completely analytic equation of state for the square-well chain fluid of variable well width

A completely analytic perturbation theory equation of state for the freely-jointed square-well chain fluid of variable well width (1 ≤ λ ≤ 2) is developed and tested against Monte Carlo simulation data. The equation of state is based on second-order Barker and Henderson perturbation theory to calculate the thermodynamic properties of the reference monomer fluid, and on first-order Wertheim thermodynamic perturbation theory to account for the connectivity of monomers to form chains. By using a recently developed real function expression for the radial distribution function of hard spheres in perturbation theory, we obtain analytic, closed form expressions for the Helmholtz free energy and the radial distribution function of square-well monomers of any well width. This information is used as the reference fluid in the perturbation theory of Wertheim to obtain an analytic equation of state, without adjustable parameters, that leads to good predictions of the compressibility factors and residual internal energies for 4-mer, 8-mer and 16-mer square-well fluids when compared with the simulation results. Further, very good results are obtained when this equation of state with temperature-independent parameters is used to correlate the vapor pressures and critical points of the linear alkanes from methane to n-decane.     

The effect of density on the properties of short chain fluids

We incorporate density dependence into continuum Born-Green-Yvon (BGY) theory through calculation of the end-to-end intramolecular correlation function. Whereas in previous studies we had only performed this calculation for the case of an isolated (zero-density) square-well chain of m segments (3相似文献   

Theoretical prediction of multiple fluid-fluid transitions in monocomponent fluids

The authors use the analytical equation of state obtained by the discrete perturbation theory [A. L. Benavides and A. Gil-Villegas, Mol. Phys. 97, 1225 (1999)] to study the phase diagram of fluids with discrete spherical potentials formed by a repulsive square-shoulder plus an attractive square-well interaction (SS+SW). This interaction is characterized by the usual energy and size parameters plus three dimensionless parameters: two of them measuring the widths of the SS and the SW and the third the relative height of the SS. The matter of interest is that, for certain values of the interaction parameters, the SS+SW systems exhibit more than one first-order fluid-fluid transition. The evidence that several real substances (such as water, phosphorus, carbon, and silica, among others) exhibit an extra liquid-liquid transition has drawn interest into the study of interactions responsible for this behavior. The simple SS+SW fluid is one of the systems that, in spite of being spherically symmetric, shows multiple fluid-fluid transitions. In this work the authors investigate systematically the effect on the phase diagram of varying the interaction parameters. The use of an analytical free-energy equation gives a clear thermodynamic picture of the emergence of different types of critical points, throwing new light on the phase behavior of these fluids and thus clarifying previous results obtained by other techniques. The interplay of attractive and repulsive forces with several scale lengths produces very rich phase diagrams, including cases with three critical points. The region of the interaction-parameter space where multiple critical points appear is mapped for various families of interactions.     

Molecular dynamical calculations on the viscosity of a square-well fluid

The coefficient of viscosity for a square-well fluid is calculated by molecular dynamics as a function of the well-depth for densities up to the region of the fluid-solid phase transition. The inclusion of an attractive contribution in the intermolecular potential has a profound influence on the behaviour of the viscosity coeffient and is also responsible for the qualitative correspondence with real systems which has been found for densities above the critical one.     

A crossover SAFT-VR equation of state for pure fluids: preliminary results for light hydrocarbons

SAFT is perhaps the most versatile, fundamentally, based engineering equation of state in use today. However, in common with all analytic equations of state, SAFT exhibits classical behavior in the critical region rather than the non-analytical, singular behavior seen in real fluids. Recently, so-called crossover equations of state have been developed which solve this shortcoming by incorporating the scaling laws valid asymptotically close to the critical point while reducing to the original classical equation of state far from the critical point. We have combined the SAFT-VR equation of state with an analytical crossover technique to obtain the SAFT-VRX equation of state. The SAFT-VRX approach combines the accurate low temperature behavior of SAFT-VR with a precise representation of the critical region. Preliminary results are presented for hydrocarbon systems which illustrate the accuracy of the SAFT-VRX approach over the entire fluid phase region.     

Deuterium isotope effect on the atomic alignment dependence in the reaction of oriented Ar (3P2) with (CH3CN)2 and (CD3CN)2 dimers

The effect of atomic alignment on CN (B2Sigma+) formation has been studied in the reaction of oriented Ar (3P2) with (CX3CN)2 (X = H, D). The reaction cross-section for each magnetic M'(J) substate in the collision frame sigma|M'(J)|(H(D),d) relative to the cross-section sigma0(H,m) in the CH(3)CN reaction was determined to be sigma0(H,d)/sigma|1|(H,d)/sigma|2|(H,d)/sigma0(D,d)/sigma|1|(D,d):/sigma|2|(D,d)= 0.87/1.00/0.98/1.58/1.93/1.78. A notable deuterium isotope effect was observed. In contrast with the monomer reactions, a significant decrease of sigma0(H(D),d) relative to the other cross-sections of sigma|M'(J)|(H(D),d) was observed.     

Equation of state for square-well chain molecules with variable range II. Extension to mixtures

An equation of state (EOS) developed in our previous work for square-well chain molecules with variable range is further extended to the mixtures of non-associating fluids. The volumetric properties of binary mixtures for small molecules as well as polymer blends can well be predicted without using adjustable parameter. With one temperature-independent binary interaction parameter, satisfactory correlations for experimental vapor–liquid equilibria (VLE) data of binary normal fluid mixtures at low and elevated pressures are obtained. In addition, VLE of n-alkane mixtures and nitrogen + n-alkane mixtures at high pressures are well predicted using this EOS. The phase behavior calculations on polymer mixture solutions are also investigated using one-fluid mixing rule. The equilibrium pressure and solubility of gas in polymer are evaluated with a single adjustable parameter and good results are obtained. The calculated results for gas + polymer systems are compared with those from other equations of state.