Variation after Particle-Number Projection for the HFB Method with the   Skyrme Energy Density Functional
  Variation after particle-number restoration is incorporated for the first time into the Hartree-Fock-Bogoliubov framework employing the Skyrme energy density functional with zero-range pairing. The resulting projected HFB equations can be expressed in terms of the local gauge-angle-dependent densities. Results of projected calculations are compared with those obtained within the Lipkin-Nogami method in the standard version and with the Lipkin-Nogami method followed by exact particle-number projection. 
M.V. Stoitsov, J. Dobaczewski, R. Kirchner, W. Nazarewicz, J. Terasaki