Cheminformatics is an interdisciplinary field that combines principles of chemistry, computer science, and information technology to process, store, analyze, and interpret chemical data. One area of cheminformatics is quantitative structure-property relationship (QSPR) modeling which is a computational approach that correlates the structural attributes of chemical compounds with their physical, chemical, or biological properties to predict the behavior and characteristics of new or untested compounds. Structure descriptors deliver contemporary mathematical tools required for QSPR modeling. One of a significant class of such descriptors is graph-based descriptors known as graphical descriptors. A degree-based graphical descriptor/invariant ( ) a -vertex graph = ( , ) has a general structure = Sigma is an element of deg , deg , where is bivariate symmetric map, and deg is the degree of vertex is an element of . For is an element of R {0}, if = (deg x deg ) (resp. = (deg + deg ) , then is called the general product-connectivity (resp. sum-connectivity ) index of . Moreover, the general Sombor index has the structure = (deg2 x deg2 ) .By choosing the heat capacity and the entropy as representatives of thermodynamic properties, we in this paper find optimal value(s) of which deliver the strongest potential of the predictors is an element of { , , } predicting and of benzenoid hydrocarbons. In order to achieve this, we employ tools such as discrete optimization and multivariate regression analysis. This, in turn, study completely solves two open problems proposed in the literature.