Generalized correlation for predicting the kinematic viscosity of liquid petroleum fractions

A two-parameter viscosity model proposed previously for pure liquids is extended to correlate the kinematic viscosity–temperature behavior for liquid petroleum fractions. The coefficients in the viscosity equation are related to the characterization properties of the petroleum fractions and a generalized kinematic viscosity–temperature correlation is then developed, which needs only specific gravity at 15.6°C and 50% boiling point as input parameters. The present method, when tested by predicting the experimental kinematic viscosities of 47 fractions from 15 world crude oils with total 250 data points, yielded reasonable results with an overall average absolute deviation of 4.2%.     

Semi-empirical relation between freezing point and critical point properties of a wide variety of molecular liquids

Based on experimental data for some 28 freezing liquids, we define a semi-empirical relation connecting surface tension and molar volume with the molar critical volume. It is shown that the semi-empirical equation deduced by means of this relation allows one to calculate the given critical property for various molecular liquids with high accuracy.     

An algorithm for mass matrix calculation of internally constrained molecular geometries

Dynamic models for molecular systems require the determination of corresponding mass matrix. For constrained geometries, these computations are often not trivial but need special considerations. Here, assembling the mass matrix of internally constrained molecular structures is formulated as an optimization problem. Analytical expressions are derived for the solution of the different possible cases depending on the rank of the constraint matrix. Geometrical interpretations are further used to enhance the solution concept. As an application, we evaluate the mass matrix for a constrained molecule undergoing an electron-transfer reaction. The preexponential factor for this reaction is computed based on the harmonic model.     

Systematically convergent correlation consistent basis sets for molecular core-valence correlation effects: the third-row atoms gallium through krypton

The family of correlation consistent polarized valence basis sets has been extended in order to account for core-core and core-valence correlation effects within the third-row, main group atoms gallium through krypton. Construction of the basis sets is similar to that of the atoms boron through argon, where either the difference between core-correlated and valence-only correlation energies were calculated via configuration interaction (CISD) computations on the ground electronic states of the atoms (named cc-pCVnZ) or the sets were optimized with respect to the core-valence correlation energy and a small weight of core-core correlation energy (cc-pwCVnZ). Due to the correlation of 3d orbitals, added shells of higher angular momentum exponents compared to the valence sets are necessary to describe the core region. The pattern of added core-correlating functions is (1s1p1d1f) for double-zeta, (2s2p2d2f1g) for triple-zeta, (3s3p3d3f2g1h) for quadruple-zeta, and (4s4p4d4f3g2h1i) for quintuple-zeta. Atomic and molecular results show good convergence to the CBS limit, with the cc-pwCVnZ sets showing improved convergence compared to the cc-pCVnZ ones for molecular core-valence correlation effects. After testing the basis sets on the homonuclear diatomics Ga2-Kr2 with coupled cluster wave functions, it is concluded that a treatment of core-valence correlation effects is essential for high-accuracy ab initio investigations of third-row-containing molecules. Though the basis sets are optimal for 3s3p3d correlation, preliminary atomic and molecular results show the basis sets to be efficient with respect to 3d-only correlation, and these potentially could be used with 3d-only correlation for more qualitative studies on larger species.     

A group-contribution correlation for predicting thermodynamic properties with the Perturbed-Soft-Chain Theory

The Perturbed-Soft-Chain theory (PSCT), an equation of state capable of predicting thermodynamic properties for a variety of compounds, is modified to be a group-contribution equation. The Group-PSCT (GPSCT) treats each pure compound as a sum of its constituent functional groups, and interactions among these groups are considered. The three pure-compound parameters, v*, T* and c, used in the original PSCT are calculated by using combining rules for five independent group parameters. The group parameters in this paper were determined from a correlation of the PSCT parameters.Using the GPSCT, calculated results show good agreement with experimental data. Liquid-densities and vapor pressures for twenty-six low molecular-weight compounds (alkanes, alkenes, and aromatics) have been calculated with average absolute errors less than 3.5% for liquid-densities and 7% for vapor-pressures. Though these values are somewhat larger than for the original PSCT, mixture predictions often are better for GPSCT than for PSCT, especially for mixtures involving polymers. This suggests that the GPSCT should prove useful for intermediate molecular weight (i.e., 200 – 1000) compounds for which little experimental data exists.     

An accurate empirical correlation for predicting natural gas viscosity

Natural gas viscosity is an important parameter in many gas and petroleum engineering calculations.This study presents a new empirical model for quickly calculating the natural gas viscosity.The model was derived from 4089 experimental viscosity data with varieties ranging from 0.01 to 21,and 1 to 3 of pseudo reduced pressure and temperature,respectively.The accuracy of this new empirical correlation has been compared with commonly used empirical models,including Lee et al.,Heidaryan et al.,Carr et al.,and Adel Elsharkawy correlations.The comparison indicates that this new empirical model can predict viscosity of natural gas with average absolute relative deviation percentage AARD (%) of 2.173.     

Density scaling of the transport properties of molecular and ionic liquids

Casalini and Roland [Phys. Rev. E 69, 062501 (2004); J. Non-Cryst. Solids 353, 3936 (2007)] and other authors have found that both the dielectric relaxation times and the viscosity, η, of liquids can be expressed solely as functions of the group (TV?(γ)), where T is the temperature, V is the molar volume, and γ a state-independent scaling exponent. Here we report scaling exponents γ, for the viscosities of 46 compounds, including 11 ionic liquids. A generalization of this thermodynamic scaling to other transport properties, namely, the self-diffusion coefficients for ionic and molecular liquids and the electrical conductivity for ionic liquids is examined. Scaling exponents, γ, for the electrical conductivities of six ionic liquids for which viscosity data are available, are found to be quite close to those obtained from viscosities. Using the scaling exponents obtained from viscosities it was possible to correlate molar conductivity over broad ranges of temperature and pressure. However, application of the same procedures to the self-diffusion coefficients, D, of six ionic and 13 molecular liquids leads to superpositioning of poorer quality, as the scaling yields different exponents from those obtained with viscosities and, in the case of the ionic liquids, slightly different values for the anion and the cation. This situation can be improved by using the ratio (D∕T), consistent with the Stokes-Einstein relation, yielding γ values closer to those of viscosity.     

Extended equation for the coexistence curve of molecular liquids in the vicinity of a critical point

An extended equation of state for substances along a phase interface, based on the Van der Waals model of a gas with fluctuations of the order parameter is confirmed by experimental data along the coexistence curve of a wide class of homogeneous and inhomogeneous molecular liquids under the Earth’s field of gravity. It is shown that the parameters of the extended equation for the coexistence curve of homogeneous and spatially inhomogeneous molecular liquids are linear functions of the compressibility factor. This allows us to predict the parameters of the equations of state of molecular liquids that are difficult to investigate experimentally.     

Prediction of melt rheological properties from GPC molecular weights

The SCLAIR^® solution polymerization platform produces a wide variety of ethylene-α-olefin copolymers and polyethylene homopolymers. Commercial products exhibit density and melt index values ranging from about 0.920 to 0.962 g/cm³ and 0.3–75 g/10 min respectively. Polymer molecular weight distributions can be tailored to meet a broad selection of end-use requirements. In this study, we have used a chemometric analysis approach using The Unscrambler^® software to demonstrate statistical correlations between rheological properties and fundamental structural parameters for thirty-three commercial SCLAIR polyethylenes. We demonstrate that molten rheological properties such as melt index, stress exponent, zero-shear viscosity, characteristic relaxation time, cross-over modulus and frequency show good non-linear correlations with molecular weight characteristics of SCLAIR products as determined by gel permeation chromatography (GPC). We also show that, with the use of Partial Least Squares (PLS) regression techniques, most melt rheological properties can be accurately predicted on the basis of GPC data.     

An accurate empirical correlation for predicting natural gas compressibility factors

The compressibility factor of natural gas is an important parameter in many gas and petroleum engineering calculations. This study presents a new empirical model for quick calculation of natural gas compressibility factors. The model was derived from 5844 experimental data of compressibility factors for a range of pseudo reduced pressures from 0.01 to 15 and pseudo reduced temperatures from 1 to 3. The accuracy of the new empirical correlation has been compared with commonly used existing methods. The comparison indicates the superiority of the new empirical model over the other methods used to calculate compressibility factor of natural gas with average absolute relative deviation percent (AARD%) of 0.6535.     

On the density scaling of pVT data and transport properties for molecular and ionic liquids

In this work, a general equation of state (EOS) recently derived by Grzybowski et al. [Phys. Rev. E 83, 041505 (2011)] is applied to 51 molecular and ionic liquids in order to perform density scaling of pVT data employing the scaling exponent γ(EOS). It is found that the scaling is excellent in most cases examined. γ(EOS) values range from 6.1 for ammonia to 13.3 for the ionic liquid [C(4)C(1)im][BF(4)]. These γ(EOS) values are compared with results recently reported by us [E. R. López, A. S. Pensado, M. J. P. Comu?as, A. A. H. Pádua, J. Fernández, and K. R. Harris, J. Chem. Phys. 134, 144507 (2011)] for the scaling exponent γ obtained for several different transport properties, namely, the viscosity, self-diffusion coefficient, and electrical conductivity. For the majority of the compounds examined, γ(EOS) > γ, but for hexane, heptane, octane, cyclopentane, cyclohexane, CCl(4), dimethyl carbonate, m-xylene, and decalin, γ(EOS) < γ. In addition, we find that the γ(EOS) values are very much higher than those of γ for alcohols, pentaerythritol esters, and ionic liquids. For viscosities and the self-diffusion coefficient-temperature ratio, we have tested the relation linking EOS and dynamic scaling parameters, proposed by Paluch et al. [J. Phys. Chem. Lett. 1, 987-992 (2010)] and Grzybowski et al. [J. Chem. Phys. 133, 161101 (2010); Phys. Rev. E 82, 013501 (2010)], that is, γ = (γ(EOS)/φ) + γ(G), where φ is the stretching parameter of the modified Avramov relation for the density scaling of a transport property, and γ(G) is the Gru?neisen constant. This relationship is based on data for structural relaxation times near the glass transition temperature for seven molecular liquids, including glass formers, and a single ionic liquid. For all the compounds examined in our much larger database the ratio (γ(EOS)/φ) is actually higher than γ, with the only exceptions of propylene carbonate and 1-methylnaphthalene. Therefore, it seems the relation proposed by Paluch et al. applies only in certain cases, and is really not generally applicable to liquid transport properties such as viscosities, self-diffusion coefficients or electrical conductivities when examined over broad ranges of temperature and pressure.     

Practical quantum mechanics-based fragment methods for predicting molecular crystal properties

Significant advances in fragment-based electronic structure methods have created a real alternative to force-field and density functional techniques in condensed-phase problems such as molecular crystals. This perspective article highlights some of the important challenges in modeling molecular crystals and discusses techniques for addressing them. First, we survey recent developments in fragment-based methods for molecular crystals. Second, we use examples from our own recent research on a fragment-based QM/MM method, the hybrid many-body interaction (HMBI) model, to analyze the physical requirements for a practical and effective molecular crystal model chemistry. We demonstrate that it is possible to predict molecular crystal lattice energies to within a couple kJ mol(-1) and lattice parameters to within a few percent in small-molecule crystals. Fragment methods provide a systematically improvable approach to making predictions in the condensed phase, which is critical to making robust predictions regarding the subtle energy differences found in molecular crystals.     

Correlation between rheological properties of zinc carboxylate liquids and molecular structure

Using Fourier transform infrared spectroscopy (FTIR), the viscosity of zinc 2-ethylhexanoate liquid has been found to correlate with the intensity of an asymmetric COO stretching resonance at 1632 cm-1. This is consistent with the presence of the zinc carboxylate polymer, catena-2-ethylhexanato-(O,O')-tri-micro-2-ethylhexanato-(O,O')-dizinc(II) as the origin of the viscosity, a conclusion that is further supported by theoretical predictions. Density functional theory has been used to assign the IR spectra of the zinc carboxylate dimer, catena-2-ethylhexanato-(O,O')-di-(tri-micro-2-ethylhexanato-(O,O')-dizinc(II)-formic) acid, and the model of the molecular liquid, micro-4-oxo-hexakis-(micro-2-ethylhexanato)-tetrazinc(II). The predicted spectra indicate that the decreased symmetry of the polymer relative to the zinc 2-ethylhexanoate liquid increases the intensity of the asymmetric carboxylate stretch at 1632 cm-1 and leads to the observed correlation.     

On the molecular theory of relaxation processes and the viscoelastic properties of magnetic liquids

The kinetic equations for one-and two-particle distribution functions including the contributions of spatial correlation of density and correlation of velocities were used to obtain equations of generalized hydrodynamics of magnetic liquids, whose transfer coefficients microscopically depended on the spatial and time scales. (The relaxing flows in these equations comprise kinetic and potential parts, which ensures the inclusion of translational and structural relaxation processes.) The system of generalized hydrodynamics equations obtained can be used to study transfer phenomena in magnetic liquids. Smoluchowski equations for binary density n ₂(q ₁, q ₂, t) and binary flow of particles I ₂(q ₁, q ₂, t) were obtained and their general solutions found to construct a closure of the initial kinetic equation for the one-particle distribution function. The asymptotic behavior of these solutions at low frequencies (ω → 0) was analyzed and found to coincide with the long-time asymptotes of autocorrelation functions. The viscoelastic properties of magnetic liquids were studied over a wide range of frequencies in the presence of an external magnetic field H. Microscopic equations for viscosity coefficients and elastic moduli were found and their asymptotic behavior in slow and fast processes considered.     

A simple and generalized model for predicting the density of ionic liquids

A simple and accurate model to predict the density of ionic liquids is presented. The proposed model is based on a generalized correlation that has been conveniently modified and experimental literature data have been used to fit the five model parameters, to finally propose an equation that allows predicting densities of any ionic liquid. The model uses the critical temperature, the critical volume, the normal boiling temperature and the molecular mass to estimate the density at temperatures commonly used in ionic liquid applications (270–360 K). A set of 602 density data for 146 ionic liquids has been used in the study. The results were compared with predictions of ten generalized corresponding states principle correlations available in the literature. These generalized correlations have not been applied to ionic liquids before so the appropriateness and accuracy of these models to ionic liquid density estimation are unknown until now. Results show that the new simple correlation gives low deviations and can be used with confidence in thermodynamic and engineering calculations.     

Atom-atom correlation order and its relation to the molecular properties

The atom-atom correlation order is defined and shown to have the special characteristics of the periodic change of its sign along the successive atom pairs and to be useful for the interpretation of some properties of aromatic molecules, that is, stabilities and reactivities. Moreover, the special appearances of the spectra of the alternant and non-alternant hydrocarbon isomers, e. g. naphthalene and azulene, are well understood by using this correlation order.

---

Zusammenfassung Die zwischenatomare Korrelation, deren Vorzeichen sich charakteristisch mit dem Abstand zweier Atome in Kette oder Ring ändert, wird definiert. Sie erweist sich als nützlich bei der Diskussion von Stabilität und Reaktivität von Aromaten sowie bei einem Vergleich der Spektren und Naphthalin und Azulen.

---

Résumé On définit la »corrélation interatomique« dont le signe change périodiquement avec la distance des deux atomes (en chaine ou cycle). Elle est utile à l'interprétation des stabilités et réactivités de molécules aromatiques, et reflète les caractéristiques des isomères alternant et non-alternant d'hydrocarbures, par exemple naphtalène et azulène.

     

Transport properties of room-temperature ionic liquids from classical molecular dynamics

Room-temperature ionic liquids (RTILs) have attracted much attention in the scientific community in the past decade due their novel and highly customizable properties. Nonetheless, their high viscosities pose serious limitations to the use of RTILs in practical applications. To elucidate some of the physical aspects behind transport properties of RTILs, extensive classical molecular dynamics calculations are reported. Here, in particular, bulk viscosities and ionic conductivities of butyl-methyl-imidazole based RTILs are presented over a wide range of temperatures. The dependence of the properties of the liquids on simulation parameters, e.g., system-size effects or the choice of the interaction potential, is analyzed in detail.     

Hydrogen bonds in protic ionic liquids and their correlation with physicochemical properties

The temperature dependence of the N-H proton chemical shift in protic ionic liquids (PILs) and FT-IR spectra of the N-H bonds indicated the presence of strong hydrogen bonds between the protonated cation and the anion, depending on the ΔpK(a) of the constituent acid and base, and their successive breaking with temperature, which may explain the characteristic properties of PILs such as relatively low ionicity and its decrease with temperature.     

Rheological properties of aqueous solutions of polyethylene glycols with various molecular weights

The dependences of the kinematic viscosity of dilute aqueous solutions of polyethylene glycols (PEGs) with various molecular weights in the temperature range of 293.15–323.15 K were investigated. The intrinsic viscosity, the Huggins constant, and the activation energy of a viscous flow were calculated for these solutions. Proposals regarding the structure of polymer macromolecules in solution are made. The constants of the Mark-Kuhn-Hauvink equation required to estimate the polymer molecular weights were determined for PEG-water systems at various temperatures.