In the code HFBTHO (v1.66p), the HFB energy (9) can be minimized
under the constraint of a fixed quadrupole moment. This option
should be used if one is interested in the potential energy surface
of a nucleus along the quadrupole collective coordinate.
The quadrupole constraint is assumed in the standard
quadratic form [40]: