Together with the FORTRAN source code in the file hfodd.f, an example of the output file is provided in la132-a.out. Selected lines from this file are presented in section TEST RUN OUTPUT below. This output file corresponds to the input file la132-a.dat reproduced in section TEST RUN INPUT below.
The output file contains two consecutive runs of the code HFODD.
The first run starts from the Nilsson model, see Sec. I-3.10, and in
95 iterations finds a triaxial solution in La for
=0. The triaxial shape is forced by using
appropriate constraints for moments
and
. Such a solution is obtained with parity
and signature symmetries conserved, and with the time-reversal broken
in order to describe odd numbers of neutrons and protons. Therefore,
even for
=0 one obtains solution with a non-zero
value of the angular momentum, which is oriented along the
axis.
Since this first run is performed within conditions that have already
been available in the previous version (v1.75r), in section TEST RUN
OUTPUT we skip the corresponding output lines.
The second run starts from the solution obtained in the first run and
performs calculations with all symmetries broken except parity, and
for the initial
angular momentum vector of
MeV.
Since for such a low value of the angular frequency, only a planar
solution exists [20], the nucleus has to turn in space in
such way that the angular momentum vector coincides with the symmetry
plane spanned by the short and long axes of the mass distribution.
Standard iteration, with IMOVAX=0, see Secs. 2.3 and
3.5, requires several thousand iterations to achieve the
alignment of the angular momentum and angular frequency vectors
below the angle of 0.01
. On the other hand, iteration with
IMOVAX=1, as shown in the example of file la132-a.dat,
achieves the alignment of 0.000001
within 92 iterations.
Section SKYRME FORCE DEFINITION lists the name and parameters of the Skyrme force together with parameters KETA_J, KETA_W, KETACM, and KETA_M that define the way the given force should be used (see Sec. 3.1).
Section CALCULATIONS WITH THE TILTED-AXIS CRANKING gives
values of components of the angular frequency vector and its length.
For switches IMOVAX=1 or 1 and IOCONT=0,
these vales are used only in the first iteration, and later ignored,
because the angular frequency vector is in each iteration readjusted
to be aligned or anti-aligned with the angular momentum vector.
Section PARITY CONFIGURATIONS gives numbers of neutron and proton states in the two parity blocks, see Sec. 3.3.
Section CONVERGENCE REPORT gives information on the performed iterations. One line per one iteration lists values described in Sec. II-4, except that version (v2.07f) prints the multipole moments in the intrinsic frame, the last before last column gives the angle between the angular frequency and angular momentum vectors, and the last column gives the total pairing energy.
Section SINGLE-PARTICLE PROPERTIES gives information on the s.p. states. For broken simplex symmetry, three projections of the total angular momentum and intrinsic spin are printed in the first and second line, respectively, for each s.p. state.
Section EULER ANGLES OF THE PRINCIPAL-AXES FRAME gives the
standard Euler angles ,
, and
[8]
in degrees, see Sec. 2.4, which define the orientation of
the principal-axes (intrinsic) frame of reference with respect to the
original frame.
Section MULTIPOLE MOMENTS [UNITS: (10 FERMI)^ LAMBDA] gives values of multipole moments with respect to the original frame of reference. Similarly, section MULTIPOLE MOMENTS [UNITS: (10 FERMI)^ LAMBDA] [INTRINSIC FRAME] gives analogous values with respect to the intrinsic frame of reference. Whenever the parity symmetry is broken, yet another analogous section gives information on values of multipole moments with respect to the center-of-mass reference frame.
Sections MAGNETIC MOMENTS [MAGNETON*FERMI^ (LAMBDA-1)] and MAGNETIC MOMENTS [MAGNETON*FERMI^ (LAMBDA-1)] [INTRINSIC FRAME] give values of magnetic moments with respect to the original and intrinsic reference frame, respectively.
Section ANGULAR MOMENTA IN THE THREE CARTESIAN DIRECTIONS gives values of neutron, proton, and total projections of the total angular momentum and intrinsic spin on the three Cartesian axes.
Sections NEUTRON CONFIGURATIONS and PROTON CONFIGURATIONS give visual representation of states occupied in the parity blocks, in analogy to those described in Sec. I-4.
Section ENERGIES gives summary of energies calculated in the last iteration. Apart from entries described in Sec. I-4, it gives in addition the pairing rearrangement energies P-REARR, pairing gaps PAIRGAP, Fermi energies E-FERMI, and surface multipole constraint energies CONSTR. (SURF).
Alternatively, one can perform the second run with the
-simplex
conserved (ISIMTX=1) and the initial angular
frequency vector of
MeV and
=0. In this way, the angular momentum is confined
within the symmetry plane
-
spanned by the short and long axes
of the mass distribution. The corresponding input and output data
files are provided in files la132-b.dat and
la132-b.out, respectively. Since the planar solution in
question corresponds to projections of the angular momentum on the
short and long axes that are not equal to one another, the
nucleus has now turn in space about the
axis. For IMOVAX=0,
this process is also very slow, while for IMOVAX=1, as in file
la132-b.dat it takes again only 92 iterations. Apart from
different orientations in space, solutions obtained by using input
data files la132-a.dat and la132-b.dat are entirely
equivalent. Output file la132-b.out is not shown in section
TEST RUN OUTPUT below.