This Python tool calculates Hansen Solubility Parameters (HSP) using the GCVOL method for molar volume calculations. It extends and modifies the original Python code for HSP calculations from "Pencil and Paper Estimation of Hansen Solubility Parameters" (Mathieu, ACS Omega, 2018) by integrating the GCVOL method for improved molar volume estimation.
The tool provides calculations for:
- Molar volume using the GCVOL method
- Hansen Solubility Parameters (HSP)
- Dispersive component (δD)
- Polar component (δP)
- Hydrogen bonding component (δH)
- Python 3.x
- Required packages:
- rdkit
- pandas
- numpy
The tool requires two Excel files:
GCVOL_Group_Contributions.xlsx
: Contains GCVOL group contribution parameters (part of this repository)Database.xlsx
: Contains substance database
The tool requires Database.xlsx to contain your compounds data with the following columns:
Column A: Name/identifier of the compound Column B: SMILES notation of the molecular structure Column C: Temperature in Kelvin
Important Note: The identifier in Column A must match exactly with a compound name in the GCVOL table. The tool uses this identifier to look up the pre-calculated group contributions - it does not perform automatic group contribution analysis of the molecule. You need to manually determine the groups present in your molecule and ensure they match with an entry in the GCVOL table.
The molar volume is calculated using the GCVOL method, which uses group contributions and temperature-dependent parameters:
V_mol = Σ(ni * (Ai + Bi * T * 10-3 + Ci * T2 * 10-5))
where:
- ni: Number of groups of type i
- Ai, Bi, Ci: GCVOL parameters
- T: Temperature in Kelvin
The three HSP components are calculated using:
- Dispersive (δD): Based on molar refractivity and volume
- Polar (δP): Using group contribution method
- Hydrogen bonding (δH): Using group contribution method
This software is based on and extends two key methodologies:
- Originally described in: Mathieu, D. "Pencil and Paper Estimation of Hansen Solubility Parameters". ACS Omega, 2018, 3 (12)
- Core HSP calculation functions have been adapted from the original implementation
- The original code has been modified under academic fair use for research purposes
- Based on: Ihmels, E. C.; Gmehling, J. "Extension and Revision of the Group Contribution Method GCVOL for the Prediction of Pure Compound Liquid Densities". Ind. Eng. Chem. Res., 2002, 42 (2)
- Implementation of temperature-dependent molar volume calculations