Calculation of properties¶
The COSMORS 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 COSMORS theory that was discussed before and Ref. 1. In the ADF COSMORS implementation the following equations were used to calculate properties using these chemical potentials.
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 COSMORS method it is possible to predict vapor pressures. In the COSMORS 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 COSMORS 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 COSMORS 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 GibbsHelmholtz 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, GibbsHelmholtz equation 
\(G^{mix}\) 
The Gibbs energy of mixing 
\(\Delta_{vap} H\) 
The enthalpy of vaporization, ClausiusClapeyron equation 
\(E_i^{HB~pure}\) 
The hydrogen bond energy of the pure compound i in the liquid phase, see Ref. 1 
\(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. 1 
\(E_i^{misfit}\) 
The partial misfit energy of compound i in a liquid solution 
\(\Delta G_{solv}^{liqsolv}\) 
The solvation Gibbs free energy from the pure compound liquid phase 
to the solvated phase 

\(\Delta G_{solv}^{gassolv}\) 
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 COSMORS GUI tutorial for the calculation of the following properties:
boiling point of a solvent [1]
partition coefficients (log P) [1, 2], OctanolWater partition coefficients (log P_{OW} ) [1]
activity coefficients [1, 2], solvation free energies [1], Henry’s law constants [1], \(pK_a\) values [1]
Ionic liquids in COSMORS 2020
The activity coefficient of a compound i solvated in an ionic liquid is an important thermodynamic property. In COSMORS 2020 one can treat the ionic liquid as one compound, which means that the value of the activity coefficient is calculated in the standard way most applications report them. In particular, in COSMORS 2020 one can treat the ionic liquid as one compound, which only has the dissociated form.
Ionic liquids in COSMORS <=2019
The activity coefficient of a compound i solvated in an ionic liquid is an important thermodynamic property. The cation and anion, which have been treated separately, will be used in equal amounts to ensure an electroneutral mixture in the COSMORS calculation.
In other applications cationanion pair have been treated as one molecule, however, below we will treat the cation and anion as two separate molecules, which is needed in older versions of COSMORS <=2019. This has consequences for the value of the activity coefficient.
For example, for a 1:1 IL (i.e., [A]^{+} [B]^{} ), the activity coefficient at a finite concentration of solute i in the binary mixture (IL + solute) can be calculated by
 γ_{i} ^{bin} = (γ_{i} ^{tern} x_{i} ^{tern} )/x_{i} ^{bin} = γ_{i} ^{tern} /(1+x_{IL} ^{bin} )
where the superscript “tern” represents the hypothetical ternary system comprising cation, anion and solute i, with
 x_{cation} ^{tern} = x_{anion} ^{tern}
 x_{cation} ^{tern} + x_{anion} ^{tern} + x_{i} ^{tern} = 1
and the superscript “bin” represents the binary mixture comprising solute and IL, with
 x_{IL} ^{bin} + x_{i} ^{bin} = 1
Accordingly, the activity coefficient of a solute i in the binary mixture (IL + solute) at infinite dilution is simplified as
 γ_{i} ^{bin} = 0.5 γ_{i} ^{tern} (at infinite dilution)
Thus in this case we should scale the activity coefficient at infinite dilution γ_{i} ^{tern} , which is directly obtained from the COSMORS calculation, with a factor of 0.5.
Similarly, for a ternary system comprising component i, component j and an ionic liquid, the activity coefficient at finite concentration of component i can be calculated by
 γ_{i} ^{tern} = γ_{i} ^{quart} /(1+x_{IL} ^{tern} )
where the superscript “quart” represents the hypothetical quaternary system comprised of cation, anion, solute i and solute j, with:
 x_{cation} ^{quart} = x_{anion} ^{quart}
 x_{cation} ^{quart} + x_{anion} ^{quart} + x_{i} ^{quart} + x_{j} ^{quart} = 1
and the superscript “tern” represents the ternary mixture comprising solute i, j, and IL, with
 x_{IL} ^{tern} + x_{i} ^{tern} + x_{j} ^{tern} = 1
References
 1(1,2,3)
A. Klamt, V. Jonas, T. Bürger and J.C. Lohrenz, Refinement and Parametrization of COSMORS. J. Phys. Chem. A 102, 5074 (1998)