Density functional theory in Firefly is implemented for restricted, unrestricted, and restricted-open wavefunctions. For each wavefunction, energies and analytical gradients are available. Analytical Hessians are currently not available. Instead, Hessians can be calculated in a semi-numerical fashion. Double hybrid functionals are an exception – for these functionals, only energies are currentl...