Various contributions to the osmotic second virial coefficient in protein-water-cosolvent solutions

An analysis of the cosolvent concentration dependence of the osmotic second virial coefficient (OSVC) in water-protein-cosolvent mixtures is developed. The Kirkwood-Buff fluctuation theory for ternary mixtures is used as the main theoretical tool. On its basis, the OSVC is expressed in terms of the thermodynamic properties of infinitely dilute (with respect to the protein) water-protein-cosolvent mixtures. These properties can be divided into two groups: (1) those of infinitely dilute protein solutions (such as the partial molar volume of a protein at infinite dilution and the derivatives of the protein activity coefficient with respect to the protein and water molar fractions) and (2) those of the protein-free water-cosolvent mixture (such as its concentrations, the isothermal compressibility, the partial molar volumes, and the derivative of the water activity coefficient with respect to the water molar fraction). Expressions are derived for the OSVC of ideal mixtures and for a mixture in which only the binary mixed solvent is ideal. The latter expression contains three contributions: (1) one due to the protein-solvent interactions B2(p-s), which is connected to the preferential binding parameter, (2) another one due to protein/protein interactions (B2(p-p)), and (3) a third one representing an ideal mixture contribution (B2(id)). The cosolvent composition dependencies of these three contributions were examined for several water-protein-cosolvent mixtures using experimental data regarding the OSVC and the preferential binding parameter. For the water-lysozyme-arginine mixture, it was found that OSVC exhibits the behavior of an ideal mixture and that B2(id) provides the main contribution to the OSVC. For the other mixtures considered (water-Hm MalDH-NaCl, water-Hm MalDH-(NH4)2SO4, and water-lysozyme-NaCl mixtures), it was found that the contribution of the protein-solvent interactions B2(p-s) is responsible for the composition dependence of the OSVC on the cosolvent concentration, whereas the two remaining contributions (B2(p-p)) and B2(id)) are almost composition independent.     

Hydrophobic interactions and osmotic second virial coefficients for methanol in water

A recent theory of the hydrophobic effect together with a simple model for an alcohol molecule is used to calculate the osmotic (McMillan-Mayer) second virial coefficientB ₂ for methanol dissolved in water. We use this calculation to study the validity of common arguments that try to draw microscopic structural information from experimental virial coefficient data. In disagreement with many workers, we find that the hydrophobic interaction between hard spheres in water is attractive and that its strength diminishes as temperature is raised. Models that have come to the opposite conclusions have neglected complications inherent to real solutes such as the role of the hydroxy groups in affecting the correlations between the apolar portions of neighboring alcohols. The calculations reported here indicate that this neglect is a poor approximation for methanol. Our calculations also show that osmotic virial coefficients are sensitive to subtle details in the potentials of mean force. Therefore, slowly varying (e.g., dispersion) interactions may also contribute significantly to the values of these coefficients without significantly changing the solvent structure near the solute molecules.     

Generalised virial equation of state for natural gas systems

In this work we present a generalised virial equation of state for natural gas systems under custody transfer conditions. The model is based on corresponding states expressions for the second and third virial coefficients with argon as the reference fluid. These functional forms involve 12 adjustable coefficients. For the extension to mixtures we propose a one-fluid mixture model with binary interaction parameters in the combining rules for the mixture critical temperature and density. We obtained overall average absolute deviations (AAD) of 0.04 and 0.08% in pure-fluid compression factors and speeds of sound; AADs of 0.07 and 0.19% in compression factors and speeds of sound, respectively, of binary mixtures and AADs of 0.047, and 0.13% in natural gas compression factors and speeds of sound, respectively. These results compare favourably with equivalent calculations with other generalised virial coefficient models.     

An application of the modified Poisson- Boltzmann equation in studies of osmotic properties of micellar solutions

The spherical cell model of colloidal solutions is applied in calculations of the osmotic pressure of micellar systems. The predictions of the nonlinear Poisson-Boltzmann equation (MPB) and of the Modified Poisson-Boltzmann equation (MPB) containing the leading terms of the fluctuation potential and the exclusion volume corrections to the mean potential acting on simple ions are compared with the results of recent computer simulations. Both PB and MPB seem satisfactory for solutions with monovalent counterions while the MPB is preferable for studies of the solutions containing divalent countenons.On leave from the University of Ljubljana, Ljubljana, Yugoslavia     

Interdependence of enthalpic and entropic contributions to the second osmotic virial coefficient

The interdependence of the enthalpic contribution A_{2, H} and the entropic contribution A_{2, s} to the second osmotic virial coefficient for a given polymer-solvent system has been investigated from the experimental and the theoretical point of view. Experimentally, the following common facts were observed for various systems at temperatures and pressures below the critical values for the solvent. Both the isobaric and isothermal dependences can be approximated over relatively wide ranges of A_{2, H} by linear relations with a slope deviating only slightly, but in a characteristic manner from a value of ?1. When the temperature is increased at constant pressure one moves along an isobar towards higher A_{2, H}; when the pressure is increased at constant temperature, one moves along an isotherm in the opposite direction, i.e., towards lower A₂, _H. Theoretically this behavior can be described in a qualitative manner, starting from a relation derived by Patterson and Delmas on the basis of the Prigogine corresponding-states theory. The reasons for the lack of quantitative agreement are discussed.     

Ab initio virial equation of state for argon using a new nonadditive three-body potential

An ab initio nonadditive three-body potential for argon has been developed using quantum-chemical calculations at the CCSD(T) and CCSDT levels of theory. Applying this potential together with a recent ab initio pair potential from the literature, the third and fourth to seventh pressure virial coefficients of argon were computed by standard numerical integration and the Mayer-sampling Monte Carlo method, respectively, for a wide temperature range. All calculated virial coefficients were fitted separately as polynomials in temperature. The results for the third virial coefficient agree with values evaluated directly from experimental data and with those computed for other nonadditive three-body potentials. We also redetermined the second and third virial coefficients from the best experimental pρT data utilizing the computed higher virial coefficients as constraints. Thus, a significantly closer agreement of the calculated third virial coefficients with the experimental data was achieved. For different orders of the virial expansion, pρT data have been calculated and compared with results from high quality measurements in the gaseous and supercritical region. The theoretically predicted pressures are within the very small experimental errors of ±0.02% for p ≤ 12 MPa in the supercritical region near room temperature, whereas for subcritical temperatures the deviations increase up to +0.3%. The computed pressure at the critical density and temperature is about 1.3% below the experimental value. At pressures between 200 MPa and 1000 MPa and at 373 K, the calculated values deviate by 1% to 9% from the experimental results.     

On the suitability of the virial equation for modeling the solubility of solids in supercritical fluids

Five model systems, the van der Waals fluid, the Soave-Redlich-Kwong fluid, the Peng-Robinson fluid, the hard-sphere fluid, and the square-well fluid, are used to examine the performance of the truncated virial expansion in describing the fugacity of a solute at infinite dilution in a solvent. It is demonstrated that the virial fugacity results deteriorate at significantly lower densities as the solute becomes larger. This has consequences for attempts to describe the solubility of solids in supercritical fluids, where the virial expansion, truncated after the third virial coefficient, has been considered as a modeling option. The results of this work suggest that, for the densities and solute-to-solvent size ratios commonly encountered in supercritical extraction, the truncated virial expansion should not be expected to describe correctly the solute fugacity, and therefore any success it has in fitting solubility data should be viewed with caution.     

Osmotic pressure,atomic pressure and the virial equation of state of polymer solutions: Monte Carlo simulations of a bead-spring model

A recently introduced coarse-grained model of polymer chains is studied analyzing various contributions to the pressure as obtained from the virial theorem as a function of chain length N, temperature T and density ϕ. The off-lattice model of the polymer chains has anharmonic springs between the beads, but of finite extensibility, and the Morse-type interaction between beads is repulsive at very short distances and attractive at intermediate distances. Solvent molecules are not explicitly included. It is found that the covalent forces along the chain (modelled by the spring potentials) contribute a negative term to the pressure, irrespective of temperature, which vanishes linearly in ϕ as ϕ → 0. In contrast, both contributions to the pressure due to intrachain nonbonded forces and due to forces between different chains change sign from high temperatures (T ≫ θ, θ the theta-temperature) where they are positive, to low temperature where both parts of the pressure become negative. It is shown that the total pressure has the expected behavior with temperature near the θ-temperature, i.e., Δp ≡ p_tot − k_B · T_p ∼ (T − θ). We study also the concentration and chainlength dependence of the various contributions to the pressure in the good solvent regime and interpret them with scaling predictions.     

Pressure dependence of the second virial coefficient of dilute polystyrene solutions

The pressure derivatives of the second virial coefficients [dA₂/dP; 0.1 ≤ P (MPa) ≤ 35.0] for dilute polystyrene (PS) solutions in good, θ, and poor solvents were measured with static light scattering. The solvent quality improved (dA₂/dP > 0) in the good and poor solvents that we investigated (toluene, chloroform; and methylcyclohexane) but deteriorated (dA₂/dP < 0) in θ solvents (cyclohexane and 50‐50 cis,trans‐decalin). The effects of temperature [22 < T (°C) < 45] and molecular weight [25 × 10³ < weight‐average molecular weight (amu mol^?1) < 900 × 10³] on dA₂/dP for PS/cyclohexane solutions were examined. © 2003 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 41: 3070–3076, 2003     

Inverse power potentials: virial coefficients and a general equation of state

Virial coefficients up to the seventh are calculated for pair potentials depending on inverse powers of separation, for inverse powers from 5 to 80. Unlike the limiting (infinite inverse power) hard-sphere potential, some virial coefficients for finite inverse power potentials are found to be negative. This makes resummation of the virial series for general inverse power potentials more difficult than that for hard spheres, and some alternative resummation methods are presented and compared. A general equation of state is proposed for fluids of particles interacting through inverse power pair potentials, for inverse powers greater than about 10. This includes the "molecular" inverse power of 12, for which the current results support and extend the results of previous studies.     

The second and the third virial coefficients of athermal polymer solutions

The second and the third virial coefficients in the lattice model of athermal mixtures of molecules of different sizes are calculated. All computations have been done for two- and three-dimensional simple square and simple cubic lattices.     

牛血清蛋白在NaCl水溶液中渗透压的研究

将NaCl水溶液中带电蛋白质分子间的静电排斥作用用Yukawa位能函数描述,蛋白质分子间的色散作用用Lennard-Jones位能函数描述,代替了经典DLVO理论。基于McMillan-Mayer渗透压统计理论,将Duh和Mier-Y-Teran的Yukawa状态方程与Cotterman等的Lennard-Jones微扰公式相结合,建立了一个新的状态方程,并采用该状态方程研究了含有NaCl的牛血清蛋白(简称BSA)水溶液的渗透压,依据BSA在0.15mol/LNaCl水溶液中的回归参数可预测其在1～5mol/LNaCl水溶液中的渗透压。该状态方程采用了两个回归参数(硬球直径和LJ参数),具有一定的关联和预测效果,并与己有的其它理论模型进行了比较。     

多元电解质溶液表面张力的新型预测方程

将半理想溶液理论和Butler方程相结合建立了预测多元电解质溶液表面张力的新型线性预测方程.新方程可由二元系数据预测多元系的表面张力数据,而不涉及任何多元交互作用参数.利用不同温度下24个混合电解质溶液的表面张力数据对新方程进行了系统检验.结果表明新方程可利用298.15K时二元系的渗透压系数和不同温度下二元系的表面张力数据预测不同温度下高浓度的多元系的表面张力数据,且预测结果与实验数据符合得很好,并且预测结果普遍优于基于Pitzer方程的表面张力模型.     

A coupled Ramani equation: multiple soliton solutions

In this work, a sixth-order coupled Ramani equation is investigated. The simplified form of the Hirota’s method is used for analytic treatment of this equation. The constraint condition between coefficients of the spatial and the temporal variables is treated. Multiple soliton solutions and multiple singular soliton solutions are formally derived for this model.     

On theAbraham equation for salt activities in concentrated solutions

Summary It is shown that the equations derived byVoigt andAbraham for the calculation of salt activities in concentrated solutions are equivalent and that they satisfy theGibbs-Duhem equation.

---

ZurAbraham-Gleichung zur Berechnung der Aktivitäten von Salzen in konzentrierten Lösungen (Kurze Mitt.)

Zusammenfassung Es wird gezeigt, daß die vonVoigt undAbraham zur Berechnung der Aktivitäten von Salzen in konzentrierten Lösungen hergeleiteten Gleichungen äquivalent sind und dieGibbs-Duhem-Beziehung erfüllen.

     

An equation of state for hydrochloric acid solutions

Gibbons, R.M. and Laughton, A.P., 1984. An equation of state for hydrochloric acid solutions. Fluid Phase Equilibria, 18:61–68.A modified Redlich-Kwong-Soave (RKS) equation is presented for hydrochloric acid solutions. It is shown that the equation with one binary parameter predicts the azeotropic behaviour of the system accurately and the heats of solution within 10%. These results suggest that the new equation can be applied to a wide range of polar systems.     

A new cubic equation of state for prediction of VLE of polymer solutions

A new cubic equation of state (CEOS) is proposed based on temperature–pressure superposition principle. A generic CEOS form, with the Peng–Robinson parameters, is used and a temperature-dependent attractive term a(T)$a(T)$ is developed, allowing an easy calculation of thermodynamic properties and vapor–liquid equilibrium. The new equation is applied to pure polymer and polymer solutions and its results are compared with those of two others equations of state. For polymer solutions, two mixing rules without binary interaction parameters were used. The vapor–liquid equilibrium (VLE) predictions showed good agreement with experimental data as well as pressure–volume–temperature (PVT) behavior of polymer liquids, attesting the appropriate form of the new equation proposed.     

Connection between the virial equation of state and physical clusters in a low density vapor

We carry out Monte Carlo simulations of physical Lennard-Jones and water clusters and show that the number of physical clusters in vapor is directly related to the virial equation of state. This relation holds at temperatures clearly below the critical temperatures, in other words, as long as the cluster-cluster interactions can be neglected--a typical assumption used in theories of nucleation. Above a certain threshold cluster size depending on temperature and interaction potential, the change in cluster work of formation can be calculated analytically with the recently proposed scaling law. The breakdown of the scaling law below the threshold sizes is accurately modeled with the low order virial coefficients. Our results indicate that high order virial coefficients can be analytically calculated from the lower order coefficients when the scaling law for cluster work of formation is valid. The scaling law also allows the calculation of the surface tension and equilibrium vapor density with computationally efficient simulations of physical clusters. Our calculated values are in good agreement with those obtained with other methods. We also present our results for the curvature dependent surface tension of water clusters.