We here use the QUBEKit code and protocol to parametrize a bespoke force field for the CV + chromophore using ab initio data obtained with the rangeseparated hybrid CAM-B3LYP functional and 6-31G(d) basis set. QUBEKit interfaces with the Gaussian electronic structure code to obtain the necessary optimized structures, frequencies, torsion scans of the potential energy surface, and other properties. The ground state optimized structures and frequencies were obtained in vacuum, along with the torsion scan of the potential to fit the dihedral force constants.