For the benchmark datasets, we conduct hyperparameter optimizations using 20 itera- with QM features predicted from ML models, which we refer to as ml-QM-GNN models in the following sections. On the other hand, the D-MPNN models trained with the QM descriptors directly calculated by ωB97XD/def2-SVP are referred to as QM-GNN models.