Predicting mixture phase equilibria and critical behavior using the SAFT-VRX approach

The SAFT-VRX equation of state combines the SAFT-VR equation with a crossover function that smoothly transforms the classical equation into a nonanalytical form close to the critical point. By a combinination of the accuracy of the SAFT-VR approach away from the critical region with the asymptotic scaling behavior seen at the critical point of real fluids, the SAFT-VRX equation can accurately describe the global fluid phase diagram. In previous work, we demonstrated that the SAFT-VRX equation very accurately describes the pvT and phase behavior of both nonassociating and associating pure fluids, with a minimum of fitting to experimental data. Here, we present a generalized SAFT-VRX equation of state for binary mixtures that is found to accurately predict the vapor-liquid equilibrium and pvT behavior of the systems studied. In particular, we examine binary mixtures of n-alkanes and carbon dioxide + n-alkanes. The SAFT-VRX equation accurately describes not only the gas-liquid critical locus for these systems but also the vapor-liquid equilibrium phase diagrams and thermal properties in single-phase regions.     

The Ornstein-Zernike equation, critical phenomena, and the equation of state for liquids and gases

It is shown that the Ornstein-Zernike (OZ) equation has two solutions: the standard one, which depends explicitly on the interaction potential, and a second universal one, resulting from the infinity point of the partition function. It is stressed that there are two pressure components: the standard one and a universal one that is valid over the whole of the phase plane. It is concluded that the universal solution parameters depend in general on definite integrals of functions dependent on the interaction potential. In the vicinity of the critical point, however, the dependence on the interaction potential vanishes; i.e., the solution becomes fully universal. It is shown that in this range of the phase diagram, all results of the theory of critical phenomena (scaling theory) follow from the OZ equation.     

Theory of critical phenomena in liquids

A theory of critical phenomena in liquids is constructed based on postulates of classical statistical mechanics (i.e., without invoking hypotheses underlying scaling theory). This theory allows for both long-and short-wave fluctuations, arising in the vicinity of a critical point. Equations for evaluating the parameters of the critical point (critical density, critical temperature, and critical pressure) from the given interaction potential are derived. An equation of state is obtained for the critical isotherm. From this equation it follows that in the vicinity of the critical point, the critical indices have classical values; scaling indices arise away from the critical point.     

Theoretical analysis of internal fluctuations and bistability in CO oxidation on nanoscale surfaces

The bistable CO oxidation on a nanoscale surface is characterized by a limited number of reacting molecules on the catalytic area. Internal fluctuations due to finite-size effects are studied by the master equation with a Langmuir-Hinshelwood mechanism for CO oxidation. Analytical solutions can be found in a reduced one-component model after the adiabatic elimination of one variable which in our case is the oxygen coverage. It is shown that near the critical point, with decreasing surface area, one cannot distinguish between two macroscopically stable stationary states. This is a consequence of the large fluctuations in the coverage which occur on a fast time scale. Under these conditions, the transition times between the macroscopic states also are no longer separated from the short-time scale of the coverage fluctuations as is the case for large surface areas and far away from the critical point. The corresponding stationary solutions of the probability distribution and the mean first passage times calculated in the reduced model are supported by numerics of the full two-component model.     

Calculation of the vibration frequencies of alpha-quartz: the effect of Hamiltonian and basis set

The central-zone vibrational spectrum of alpha-quartz (SiO2) is calculated by building the Hessian matrix numerically from the analytical gradients of the energy with respect to the atomic coordinates. The nonanalytical part is obtained with a finite field supercell approach for the high-frequency dielectric constant and a Wannier function scheme for the evaluation of Born charges. The results obtained with four different Hamiltonians, namely Hartree-Fock, DFT in its local (LDA) and nonlocal gradient corrected (PBE) approximation, and hybrid B3LYP, are discussed, showing that B3LYP performs far better than LDA and PBE, which in turn provide better results than HF, as the mean absolute difference from experimental frequencies is 6, 18, 21, and 44 cm(-1), respectively, when a split valence basis set containing two sets of polarization functions is used. For the LDA results, comparison is possible with previous calculations based on the Density Functional Perturbation Theory and usage of a plane-wave basis set. The effects associated with the use of basis sets of increasing size are also investigated. It turns out that a split valence plus a single set of d polarization functions provides frequencies that differ from the ones obtained with a double set of d functions and a set of f functions on all atoms by on average less than 5 cm(-1).     

Scaled equations for the coexistence curve,the capillary constant and the surface tension of n-alkanes

A review of experimental data of several fluids shows that their coexistence curve follows a power law in reduced temperature at the approach of the critical point, with an universal exponent equal to 0.325, their capillary constant a power law with an universal exponent equal to 0.925 and their surface tension a power law with an universal exponent equal to 1.26. In the critical region, the concept of two-scale-factor universality was used to predict the density difference amplitude, the capillary constant amplitude, and the surface tension amplitude between near critical vapor and liquid phases. A comparison with amplitudes determined from experimental data is given. In order to extend this universality all along the liquid–gas coexistence curve from the triple point to the critical point for n-alkanes, a mean field approximation was used far away from T_C. We show that the density difference, the capillary constant and the surface tension can be calculated with a reasonable accuracy by generalized scaled equations adding only two empirical constants. A comparison between calculated and experimental data is presented.     

On the universal behavior of some thermodynamic properties along the whole liquid-vapor coexistence curve

When thermodynamic properties of a pure substance are transformed to reduced form by using both critical- and triple-point values, the corresponding experimental data along the whole liquid-vapor coexistence curve can be correlated with a very simple analytical expression that interpolates between the behavior near the triple and the critical points. The leading terms of this expression contain only two parameters: the critical exponent and the slope at the triple point. For a given thermodynamic property, the critical exponent has a universal character but the slope at the triple point can vary significantly from one substance to another. However, for certain thermodynamic properties including the difference of coexisting densities, the enthalpy of vaporization, and the surface tension of the saturated liquid, one finds that the slope at the triple point also has a nearly universal value for a wide class of fluids. These thermodynamic properties thus show a corresponding apparently universal behavior along the whole coexistence curve.     

On the nature of the Møller-Plesset critical point

It has been suggested [F. H. Stillinger, J. Chem. Phys. 112, 9711 (2000)] that the convergence or divergence of M?ller-Plesset perturbation theory is determined by a critical point at a negative value of the perturbation parameter z at which an electron cluster dissociates from the nuclei. This conjecture is examined using configuration-interaction computations as a function of z and using a quadratic approximant analysis of the high-order perturbation series. Results are presented for the He, Ne, and Ar atoms and the hydrogen fluoride molecule. The original theoretical analysis used the true Hamiltonian without the approximation of a finite basis set. In practice, the singularity structure depends strongly on the choice of basis set. Standard basis sets cannot model dissociation to an electron cluster, but if the basis includes diffuse functions then it can model another critical point corresponding to complete dissociation of all the valence electrons. This point is farther from the origin of the z plane than is the critical point for the electron cluster, but it is still close enough to cause divergence of the perturbation series. For the hydrogen fluoride molecule a critical point is present even without diffuse functions. The basis functions centered on the H atom are far enough from the F atom to model the escape of electrons away from the fluorine end of the molecule. For the Ar atom a critical point for a one-electron ionization, which was not previously predicted, seems to be present at a positive value of the perturbation parameter. Implications of the existence of critical points for quantum-chemical applications are discussed.     

Concentration dependence of the interfacial tension for aqueous two-phase polymer solutions of dextran and polyethylene glycol

We studied the interfacial tension between coexisting phases of aqueous solutions of dextran and polyethylene glycol. First, we characterized the phase diagram of the system and located the binodal. Second, the tie lines between the coexisting phases were determined using a method that only requires measuring the density of the coexisting phases. The interfacial tension was then measured by a spinning drop tensiometer over a broad range of polymer concentrations close to and above the critical point. In this range, the interfacial tension increases by 4 orders of magnitude with increasing polymer concentration. The scaling exponents of the interfacial tension, the correlation length, and order parameters were evaluated and showed a crossover behavior depending on the distance to the critical concentration. The scaling exponent of the interfacial tension attains the value 1.50 ± 0.01 further away from the critical point, in good agreement with mean field theory, but the increased value 1.67 ± 0.10 closer to this point, which disagrees with the Ising value 1.26. We discuss possible reasons for this discrepancy. The composition and density differences between the two coexisting phases, which may be taken as two possible order parameters, showed the expected crossover from mean field behavior to Ising model behavior as the critical point is approached. The crossover behavior of aqueous two-phase polymer solutions with increasing concentration is similar to that of polymer solutions undergoing phase separation induced by lowering the temperature.     

Crossover Peng-Robinson equation of state with introduction of a new expression for the crossover function

In this research, we use the original Peng-Robinson (PR) equation of state (EOS) for pure fluids and develop a crossover cubic equation of state which incorporates the scaling laws asymptotically close to the critical point and it is transformed into the original cubic equation of state far away from the critical point. The modified EOS is transformed to ideal gas EOS in the limit of zero density. A new formulation for the crossover function is introduced in this work. The new crossover function ensures more accurate change from the singular behavior of fluids inside the regular classical behavior outside the critical region. The crossover PR (CPR) EOS is applied to describe thermodynamic properties of pure fluids (normal alkanes from methane to n-hexane, carbon dioxide, hydrogen sulfide and R125). It is shown that over wide ranges of state, the CPR EOS yields the thermodynamic properties of fluids with much more accuracy than the original PR EOS. The CPR EOS is then used for mixtures by introducing mixing rules for the pure component parameters. Higher accuracy is observed in comparison with the classical PR EOS in the mixture critical region.     

Inhomogeneous polymer solutions: Theory of the interfacial free energy and the composition profile near the consolute point

The Flory–Huggins formulation of the combinatorial entropy, supplemented with residual free energy, is applied locally to obtain the interfacial free energy and the concentration profile of polymer in the interface between two demixed polymer solution phases. Two choices were investigated for the residual free energy: a “regular solution” formulation and an empirical formulation of Koningsveld for polystyrene in cyclohexane. Asymptotic, analytical solutions of the equations near the critical solution point and solutions obtained by numerical calculations are given as a function of temperature for several molecular weights. At temperatures farther below the critical temperature the equations have no solutions. The reason for this is not entirely clear. The local formulation of the free energy used here is an improved version of a previous one, which gave wrong results for asymmetric systems (polymer in a low molecular weight solvent). This newer version is consistent with our theory of critical opalescence and gives a relation between the interface “thickness” and the correlation range of the concentration fluctuations. The calculated correlation ranges were in good accord with those found experimentally by Debye, Chu, and Woerman. That the newer version of our equations for an interface gives no acceptable solutions at lower temperatures could be caused by a “collapse” of a diffuse to a sharp interface as suggested by Nose.     

Modelling aqueous solutions near the critical point of water

We review the thermodynamic properties of dilute solution near the critical point of the solvent. Two examples are discussed, a solution of a non-electrolyte and a solution of an electrolyte. The limiting behavior of the electrolyte solutions is modelled with a Debye-Huckel term in the Helmholtz free energy. The partial molar properties, in particular the volume and isobaric thermal expansion are examined in detail. The derivation of these properties is introduced by considering the geometry of the thermodynamic surfaces near to and far from the critical point of the solvent. We conclude that the properties of solutions near the solvent critical point are dominated by that feature; solution properties cannot be adequately modelled without including the functional forms associated with the critical point.     

Critical phenomena in liquids (theory)

It is shown that the Ornstein-Zernike equation has two solutions, classic (analytic) and critical (not analytic). Closure equations of the HNC, PY, etc. types well known in the theory of liquids correspond to the first solution. The second solution presupposes that the bridge functional of a system has the form of the sum B = B ^rg + B ^cr, where B ^rg is a regular (analytic) function of density and B ^cr is a critical (not analytic) function. It is shown the classic solution determines the coordinates of the critical point, critical amplitudes, and the other parameters of critical phenomena depending on the individual liquid characteristics. The critical solution determines the critical indices and the equations relating them. The regions in which classic and critical solutions are valid are separated on the phase plane by a line of singular points, at which the second derivatives of pressure experience discontinuity.     

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.     

Refining of nonionic surfactant micellization theory based on the law of mass action

Two approaches to determining critical micelle concentration (CMC) are assessed, i.e., from the inflection point in the curve for the concentration dependence of the degree of micellization and as K^1/(1–n), where K is the constant of the law of mass action and n is the aggregation number. The latter approach makes the theory simpler, while the former explicitly expresses the critical degree of micellization via the aggregation number. The concentrations of monomers and micelles are analyzed as functions of the overall concentration of a surfactant in a micellar solution. These functions look much simpler in the graphical form as compared with their complex exact analytical representation. This has resulted in derivation of simple analytical approximations for these functions, with these approximations being useful for calculations. The concentration dependence of the surfactant diffusion coefficient has been considered based on these approximations. It turned out that this dependence not only provides the known method for determining the diffusion coefficient of micelles, but also gives the possibility in principle to determine the aggregation number from the slope of the dependence of the diffusion coefficient on the inverse concentration (counted from the CMC in the CMC units). This new method for determining the aggregation number has been tested using the literature data on the diffusion coefficient of penta(ethylene glycol)-1-hexyl ether in an aqueous solution.     

Spatial hydration maps and dynamics of naphthalene in ambient and supercritical water

The hydration structures and dynamics of naphthalene in aqueous solution are examined using molecular-dynamics simulations. The simulations are performed at several state points along the coexistence curve of water up to the critical point, and above the critical point with the density fixed at 0.3 g/cm(3). Spatial maps of local atomic pair-density are presented which show a detailed picture of the hydration shell around a bicyclic aromatic structure. The self-diffusion coefficient of naphthalene is also calculated. It is shown that water molecules tend to form pi-type complexes with the two aromatic regions of naphthalene, where water acts as the H-bond donor. At ambient conditions, the hydration shell of naphthalene is comprised, on average, of about 39 water molecules. Within this shell, two water molecules can be identified as pi-coordinating, forming close to one H-bond to the aromatic rings. With increasing temperature, the hydration of naphthalene changes dramatically, leading to the disappearance of the pi-coordination near the critical point.     

Theory of the potential of a donor ion in silicon

In this paper a variational principle has been formulated for obtaining an approximate solution of Poisson's equation for the potential of a donor ion embedded in Si. A comparison of a two parameter approximate analytical potential with the exact numerical potential shows excellent agreement for several electron concentrations. It is also found that the present theory, using the concept of a “spatially variable dielectric constant,” leads to a donor-ion potential which deviates significantly from the potential of Dingle, which has been the traditional starting point of calculations of electron mobility limited by ionized-impurity scattering.     

Anwendung einiger Methoden der automatischen Zeichenerkennung auf die interpretation von Massenspektren

Some pattern recognition methods are briefly discussed from the geometric point of view: classification by distance measurements to prototypes, classification by decision-planes and classification by the nearest-neighbor-method. The last two methods give good results with regard to the automatic determination of molecular structures from low resolution mass spectra. It is emphasized, that pattern recognition methods may be useful for interpreting various types of data pools in physical and analytical chemistry.     

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.