The code HFBTHO [1] solves the self-consistent HFB equations by using the axial (2D) Transformed Harmonic Oscillator (THO) basis [2], which allows for a correct treatment of the single-quasiparticle wave function asymptotics. As discussed recently [3], the THO technique is a method of choice for performing massive nuclear structure calculations including weakly bound systems. In order to fully test the formalism, in the present study we present results obtained with the axial (2D) HFBTHO (v1.64) code compared to those obtained with two other codes: HFODD (v2.08i) [4], which uses a Cartesian (3D) Harmonic Oscillator (HO) basis and spherical (1D) HFBRAD [5], which uses a lattice of points in the radial coordinate.
In Ref. [3] we have published the first complete mass table of even-even nuclei obtained by using the THO method for the SLy4 Skyrme force [6]. Here we discuss one specific improvement of the method, and also present new results obtained with the SkP Skyrme force [7]. More details, including downloadable tables of ground-state properties can be found at http://www.fuw.edu.pl/~ dobaczew/thodri/thodri.html.