It has been noted that, to describe all of the basic atom type combinations in GAFF, some 20,000 angle parameters would be required. Such large parameter libraries are commonplace with OPLS3 containing 15,236 angle-bending parameters, with a continuing effort to expand this list as new chemistries are encountered. To generate these parameters, general FFs have to use a wide range of reference data combining experiment and QM. QM data actually already play a role in the derivation of the majority of the transferable parameters in these FFs due to the lack of experimental data available for unique chemical species and the ease of generating accurate QM data on-the-fly. While many of the equilibrium terms are collected from x-ray crystallography and NMR studies of small molecules, some have to be determined from QM predicted minimum energy structures. Force constants are then manually fit in an iterative process which aims to recreate the QM vibrational frequencies using an initial guess for the other required parameters as described in the development of CGENFF 5 and AMBER. While this method is effective, it does create interdependencies in the FF parameters as the force constants are dependent on the rest of the original parameter set, meaning that ideally all parameters should be continually updated in a self-consistent fashion until convergence is reached. Instead, we have adopted the modified Seminario method for deriving bond and angle force field parameters. The standard Seminario method derives force constants directly from the QM Hessian matrix and has been incorporated into specialized FF fitting tools for metal complexes such as the VFFDT plugin, or in the MCPB.py program which is part of AmberTools. This method estimates force constants by projecting the decomposed forces felt by an atom due to the displacement of a neighboring atom onto their mutual bond vector. However, this method results in undesirably stiff force constants due to the double counting of angle bending contributions in larger molecules. The modified method, however, accounts for an atom's chemical environment and has been shown to recreate QM vibrational frequencies with a low average error of 6.3% across all vibrational modes for a wide range of molecules. The ability to accurately derive the bonded parameters directly from the QM Hessian matrix without the need for initial parameter guesses simplifies the procedure for non-expert users by removing sources of human error and also speeds up the process making it suitable for automation. We shall also show that the derived force constants retain a low percentage error in recreating QM vibrational frequencies when combined with the rest of the QUBE FF.