Calculation of properties

The COSMO-RS method allows to calculate the (pseudo-)chemical potential of a compound in the liquid phase, as well as in the gas phase, see the the COSMO-RS theory that was discussed before and Ref.[2]. In the ADF COSMO-RS implementation the following equations were used to calculate properties using these chemical potentials.

\[\begin{split}\sum_i x_i & = \sum_i y_i^{vapor} = \sum_i w_i = 1 \\ w_i & = x_i M_i^{pure}/M^{ave} \\ M^{ave} & = \sum_i x_i M_i^{pure} \\ p_i^{pure} & = exp \{ (\mu_i^{pure}-\mu_i^{gas})/RT \} \\ p_i^{vapor} & = x_i exp \{ (\mu_i^{solv}-\mu_i^{gas})/RT \} \\ p^{vapor} & = \sum_i p_i^{vapor} \\ y_i^{vapor} & = p_i^{vapor}/p^{vapor} \\ \gamma_i & = exp \{ (\mu_i^{solv}-\mu_i^{pure})/RT \} \\ a_i & = \gamma_i x_i \\ G^E & = \sum_i x_i (\mu_i^{solv}-\mu_i^{pure}) \\ H^E & = -RT^2 \partial \{ G^E/RT \}/ \partial T \\ G^E & = H^E - T S^E \\ \Delta_{vap} H & = RT^2/p^{vapor} \partial \{ p^{vapor} \}/ \partial T \\ k_H & = 1/V_{solvent} exp \{ (\mu_i^{gas}-\mu_i^{solv})/RT \} \\ k_H^{cc} & = 1/(k_H V_{solvent}) = \gamma_i p_i^{pure} \\ x_i^{SOL} & = 1/{\gamma_i} (T > T_m) \\ x_i^{SOL} & = 1/{\gamma_i} exp \{ \Delta H_{fus}(1/T_m-1/T)/R - \Delta C_p (ln(T_m/T) -T_m/T + 1)/R \} (T < T_m) \\ \Delta G_{solv}^{liq-solv} & = \mu_i^{solv}-\mu_i^{pure} \\ \Delta G_{solv}^{gas-solv} & = \mu_i^{solv}-\mu_i^{gas} + RT ln(V_{solvent}/V_{gas}) \\ log_{10} P_{solv1/solv2} & = 1/ln(10) (\mu_i^{solv2}-\mu_i^{solv1})/RT + log_{10} (V_{solv1}/V_{solv2}) \\ 1/LFL_{mix} & = \sum_i (y_i^{vapor}/LFL_i)\end{split}\]

The above equations are not always exact, some assume perfect gas behavior, for example.

The molar fraction \(x_i\) of each compound i of the solvent should add up to 1.

With the COSMO-RS method it is possible to predict vapor pressures. In the COSMO-RS model the free energy difference between the chemical potential in the gas phase (perfect gas with a reference state of 1 bar) and the chemical potential of the liquid phase has been defined in such a way that the equation: \(p_i = exp {(\mu_i^{pure} - \mu_i^{gas} )/RT}\), will give the pressure in units of bar. It is also possible to use experimental vapor pressures of pure compounds as input data for the calculation. This may increase the accuracy of the calculated vapor pressures in a mixture, for example.

In the COSMO-RS method the volume of 1 molecule in the liquid phase is approximated with the volume of the molecule shaped cavity, that is used in the COSMO calculations. In this way it is possible to calculate the volume of 1 mole of solvent molecules in the liquid phase. However, for properties that depend on such volumes, one can also use (related) experimental data as input data for the calculation.

The calculation of the boiling temperature of a solvent is performed with an iterative method. The temperature is varied until the calculated vapor pressure is within a certain threshold of the desired pressure.

Also the calculation of solubility of compound i is performed with an iterative method, since the activity coefficient \(\gamma_i\) depends on the molar fraction of this compound. The COSMO-RS method does not predict \(\Delta H_{fus}\), \(\Delta C_p\) , or \(T_m\) . These can be given as input data for the calculation of solubility calculations of solid compounds.

Starting from ADF2012 the Gibbs-Helmholtz equation is used to calculate the excess enthalpy of a mixture. Previously it was estimated using the misfit and hydrogen bonding energy of the mixture and its pure compounds.

Quantity Meaning
\(R\) Gas constant
\(T\) Temperature
\(x_i\) The molar fraction of compound i in a liquid solution
\(y_i\) vapor The molar fraction of compound i in the gas phase
\(w_i\) The mass fraction of compound i in a liquid solution
\(M_i\) pure The molar mass of the pure compound i
\(M^{ave}\) The average molar mass of the mixture
\(\gamma_i\) Activity coefficient of compound i in a liquid solution
\(a_i\) Activity of compound i in a liquid solution
\(p_i^{pure}\) The vapor pressure of the pure compound i
\(p_i^{vapor}\) The partial vapor pressure of compound i
\(p^{vapor}\) The total vapor pressure
\(\mu\)i gas The pseudochemical potential of the pure compound i in the gas phase
\(\mu\)i pure The pseudochemical potential of the pure compound i in the liquid phase
\(\mu\)i solv The pseudochemical potential of compound i in a liquid solution
\(G^E\) The excess Gibbs free energy
\(H^E\) The excess enthalpy, Gibbs-Helmholtz equation
\(\Delta_{vap} H\) The enthalpy of vaporization, Clausius-Clapeyron equation
\(E_i^{HB~pure}\) The hydrogen bond energy of the pure compound i in the liquid phase, see Ref.[2]
\(E_i^{HB}\) The partial hydrogen bond energy of compound i in a liquid solution
\(E_i^{misfit~pure}\) The misfit energy of the pure compound i in the liquid phase, see Ref.[2]
\(E_i^{misfit}\) The partial misfit energy of compound i in a liquid solution
\(\Delta G_{solv}^{liq-solv}\) The solvation Gibbs free energy from the pure compound liquid phase
  to the solvated phase
\(\Delta G_{solv}^{gas-solv}\) The solvation Gibbs free energy from the pure compound gas phase
  to the solvated phase, with a reference state of 1 mol/L in both phases
\(k_H\) Henry’s law constant: ratio between the liquid phase concentration of a compound
  and its partial vapor pressure in the gas phase
\(k_H^{cc}\) dimensionless Henry’s law constant: ratio between the liquid phase concentration
  of a compound and its gas phase concentration
\(k_{H~inv}^{px}\) Henry’s law constant, representing the volatility instead of the solubility,
  ratio between the partial vapor pressure of a compound in the gas phas
  and the molar fraction in the liquid phase”
\(V_{solvent}\) Volume of 1 mole of solvent molecules in the liquid phase
\(V_{gas}\) Volume of 1 mole of molecules in the gas phase (at 1 atm, perfect gas)
\(x_i^{SOL}\) Solubility of compound i in a solvent (molar fraction)
\(\Delta H_{fus}\) The enthalpy of fusion of compound i
\(\Delta C_p\) The \(\Delta\) heat capacity of fusion of compound i
\(T_m\) The melting temperature of compound i
\(log_{10} P_{solv1/solv2}\) The logarithm of the partition coefficient P, which is the ratio of the concentrations
  of a compound in two immiscible solvents, solvent 1 and solvent 2
\(LFL_i\) The flash point (lower flammable limit, LFL) of compound i
\(LFL_{mix}\) The flash point (lower flammable limit, LFL) of a mixture, Le Chatelier’s mixing rule

See also the COSMO-RS GUI tutorial for the calculation of the following properties:

  • solvent vapor pressure [1, 2]
  • boiling point of a solvent [1]
  • partition coefficients (log P) [1, 2], Octanol-Water partition coefficients (log POW ) [1]
  • activity coefficients [1, 2], solvation free energies [1], Henry’s law constants [1], \(pK_a\) values [1]
  • solubility [1, 2]
  • vapor-liquid diagram binary mixture (VLE/LLE) [1, 2]