To perform a purely data-driven study, we utilize the QM7-X dataset 37 containing 42 physicochemical properties of ≈ 4.2 M (equilibrium and non-equilibrium) organic molecules with up to seven heavy (non-hydrogen) atoms (including C, O, N, S and Cl), spanning a practically important subset of CCS. Accordingly, a subset of QM7-X considering only one equilibrium constitutional isomers and stereoisomers per unique molecular graph is selected for further analysis (≈ 13 k molecules). In QM7-X, the molecular structures were optimized using the third-order self-consistent charge density-functional tight binding method (DFTB3) supplemented with a treatment of many-body dispersion/van der Waals interactions via the MBD approach . However, for our studies concerning polarizability, α was computed directly, employing density functional perturbation theory (DFPT) by means of the PBE0 81 functional as implemented in the FHI-aims code (version 190205). The tight basis set settings were also applied for all atoms (C, H, N, O, S and Cl). To ensure the transferability of the values, we store the molecular (mean) polarizability (denoted as α and simply referred to as polarizability in other parts of the manuscript)