  We describe the new version 4.0 of the code hfbtho that solves the nuclear
Hartree-Fock-Bogoliubov problem by using the deformed harmonic oscillator basis
in cylindrical coordinates. In the new version, we have implemented the
restoration of rotational, particle number, and reflection symmetry. The
restoration of rotational symmetry does not require using bases closed under
rotation. Furthermore, the SeaLL1 functional was added and the calculation of
the Coulomb potential was improved. Finally, the code has been refactored to
facilitate maintenance and future developments.
