rameterization consist in: i) a geometry optimization carried out in the target electronic state carried out specifically for each molecule; ii) a calculation of the Hessian matrix for each opti-mized geometry; iii) a number of relaxed torsional energy scans along each rotatable dihedral (vide infra) of the considered species, obtained minimizing the internal energy with respect to all ICs, but the scanned one. To demonstrate JOYCE3.0's capability to parameterize QMD-FFs for diverse chemical species and across different levels of theory used to build the QM reference database, Table summarizes the QM descriptors and methods employed in this work for the selected targets illustrated in Figure .