TEST RUN INPUT
============================================================================
| This file (la132-a.dat) contains the input data for the code HFODD. |
| It defines two consecutive runs that lead to a planar solution in 132La. |
============================================================================
------------ Nuclide -------------
NUCLIDE
75 57
---------- Convergence -----------
ITERATIONS
5000
ITERAT_EPS
0.000001
MAXANTIOSC
5
------- Skyrme interaction -------
SKYRME-STD
1 0 0 0 0
--------- Configurations ---------
VACSIG_NEU PPSP PPSM PMSP PMSM
19 19 19 18
VACSIG_PRO PPSP PPSM PMSP PMSM
15 15 14 13
------------ HO basis ------------
BASIS_SIZE
10 286 800.
SURFAC_PAR
75 57 1.230
SURFAC_DEF
2 0 0.000
SURFAC_DEF
4 0 0.000
---------- Output data -----------
MAX_MAGNET
0 3
----------- File names -----------
REVIEWFILE
la132-a.rev
REPLAYFILE
la132-a.rec
RECORDFILE
la132-a.rec
---------- Constraints -----------
OMEGAY
0.000
MULTCONSTR
2 0 0.5 7.0 1
MULTCONSTR
2 2 0.5 7.0 1
--------- Starting point ---------
RESTART
0
------------ Calculate -----------
EXECUTE
----------- Symmetries -----------
SIMPLEXY
0
SIGNATUREY
0
PARITY
1
TSIMPLEX3D
0 0 0
--------- Configurations ---------
VACPAR_NEU PARP PARM
38 37
VACPAR_PRO PARP PARM
30 27
---------- Constraints -----------
OMEGA_XYZ
0.1443375673 0.1443375673 0.1443375673 1
OMEGA_TURN
1
MULTCONSTR
2 0 0.5 7.0 0
MULTCONSTR
2 2 0.5 7.0 0
--------- Starting point ---------
RESTART
1
CONT_OMEGA
0
------------ Calculate -----------
EXECUTE
------------ Terminate -----------
ALL_DONE
TEST RUN OUTPUT
******************************************************************************* * * * HFODD HFODD HFODD HFODD HFODD HFODD HFODD HFODD * * * ******************************************************************************* * * * SKYRME-HARTREE-FOCK-BOGOLYUBOV CODE VERSION: 2.07F * * * * NO SYMMETRY-PLANES AND NO TIME-REVERSAL SYMMETRY * * * * DEFORMED CARTESIAN HARMONIC-OSCILLATOR BASIS * * * ******************************************************************************* * * * J. DOBACZEWSKI, J. DUDEK, AND P. OLBRATOWSKI * * * * INSTITUT DE RECHERCHES SUBATOMIQUES, STRASBOURG * * * * INSTYTUT FIZYKI TEORETYCZNEJ, WARSZAWA * * * * 1993-2004 * * * ******************************************************************************* ******************************************************************************* * * * SKYRME FORCE DEFINITION * * * ******************************************************************************* * * * PARAMETER SET SKM*: T0= -2645.00 T1= 410.00 T2= -135.00 T3= 15595.00 * * * * POWER=1/6 W=130.000 X0= 0.09000 X1= 0.00000 X2= 0.00000 X3= 0.00000 * * * * ETA_J=0 ETA_W=0 ETACM=0 ETA_M=1 HBM= 20.73000 (SKYRME-FORCE SPIN-ORBIT) * * * ******************************************************************************* ******************************************************************************* * * * CALCULATIONS WITH ADJUSTABLE ISOSCALAR ROTATIONAL FREQUENCY VECTOR, WHICH * * IN EACH ITERATION IS TURNED TO BE PARALLEL TO THE ANGULAR MOMENTUM VECTOR * * * * DIRECTION OF THE ROTATIONAL FREQUENCY VECTOR THAT IS GIVEN BELOW WILL BE * * USED ONLY IN THE FIRST ITERATION, AND LATER WILL BE IGNORED * * * ******************************************************************************* * * * CALCULATIONS WITH THE TILTED-AXIS CRANKING * * * * FOR THE ISOSCALAR ROTATIONAL FREQUENCIES OF: * * * * OMEGAX = 0.144338 OMEGAY = 0.144338 OMEGAZ = 0.144338 OMEGA = 0.250000 * * * ******************************************************************************* * * * PARITY CONFIGURATIONS: * * * * V A C U U M P A R T I C L E S H O L E S * * =========== ================= =========== * * ( + ) ( - ) ( + ) ( - ) ( + ) ( - ) * * * * NEUTRONS: 38 37 0 0 0 0 * * PROTONS : 30 27 0 0 0 0 * * * ******************************************************************************* ******************************************************************************* * * * CONVERGENCE REPORT * * * ******************************************************************************* * * * ITER ENERGY STABILITY Q20 GAMMA SPIN OMEGA ANGLE EPAIR * * * * 96 -1093.673718 -0.232831 6.890 45.36 11.28 0.250 13.308 0.00 * * 97 -1093.850272 0.025915 6.876 45.40 11.40 0.250 8.286 0.00 * * 98 -1093.854682 0.044356 6.863 45.43 11.59 0.250 4.645 0.00 * * 99 -1093.844776 0.041488 6.852 45.46 11.71 0.250 2.825 0.00 * * 100 -1093.832812 0.036000 6.843 45.49 11.79 0.250 1.773 0.00 * * 101 -1093.823120 0.030728 6.834 45.51 11.85 0.250 1.207 0.00 * * 102 -1093.815613 0.026215 6.827 45.52 11.89 0.250 0.857 0.00 * * 103 -1093.810028 0.022473 6.821 45.54 11.92 0.250 0.637 0.00 * * 179 -1093.792431 1.31E-06 6.780 45.62 12.00 0.250 1.7E-06 0.00 * * 180 -1093.792431 1.15E-06 6.780 45.62 12.00 0.250 1.5E-06 0.00 * * 181 -1093.792431 1.01E-06 6.780 45.62 12.00 0.250 8.5E-07 0.00 * * 182 -1093.792431 8.94E-07 6.780 45.62 12.00 0.250 1.5E-06 0.00 * * 183 -1093.792431 7.88E-07 6.780 45.62 12.00 0.250 8.5E-07 0.00 * * 184 -1093.792431 6.95E-07 6.780 45.62 12.00 0.250 1.5E-06 0.00 * * 185 -1093.792431 6.13E-07 6.780 45.62 12.00 0.250 8.5E-07 0.00 * * 186 -1093.792431 5.40E-07 6.780 45.62 12.00 0.250 0.000 0.00 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY (P=+,P=-) | N,nz,lz,OMEG> <P> <JX> <JY> <JZ> * * * * 71) -10.460 ( 0, 34) | 5, 1, 4, 9/2> -100 0.029 -0.189 2.078 * * 0.003 -0.022 0.223 * * 72) -10.338 ( 38, 0) | 4, 0, 2, 3/2> 100 -0.007 -0.811 -0.436 * * 0.005 0.036 0.372 * * 73) -10.042 ( 0, 35) | 5, 2, 3, 7/2> -100 -0.034 -0.579 -2.392 * * -0.004 -0.027 -0.285 * * 74) -9.318 ( 0, 36) | 5, 0, 5,11/2> -100 0.070 -0.630 5.021 * * 0.007 -0.034 0.495 * * 75) -9.145 ( 0, 37) | 5, 1, 4, 9/2> -100 -0.053 -0.580 -3.765 * * -0.006 -0.040 -0.429 * * 76) -8.966 ( 39, 0) | 4, 1, 1, 1/2> 100 0.005 0.589 0.287 * * -0.003 -0.011 -0.240 * * 77) -8.637 ( 40, 0) | 4, 0, 2, 5/2> 100 -0.004 -0.152 -0.282 * * -0.002 0.045 -0.168 * * 78) -8.315 ( 41, 0) | 4, 1, 1, 1/2> 100 -0.003 0.295 -0.247 * * 0.001 0.013 0.065 * * 79) -8.018 ( 0, 38) | 5, 2, 1, 1/2> -100 0.006 3.305 0.150 * * 0.001 0.332 0.034 * * 80) -7.731 ( 0, 39) | 5, 0, 5,11/2> -100 -0.068 -0.541 -4.847 * * -0.007 -0.029 -0.496 * * * ******************************************************************************* ******************************************************************************* * * * EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL * * * * ALPHA = -1.65924 BETA = 0.80185 GAMMA = 1.72034 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 =132.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 6.7791 Q2+1= 0.0194 Q2+2= 6.9293 ............. ............. * * Q2-1= -0.0022 Q2-2= -0.0147 ............. ............. * * * * Q30 = ZERO Q3+1= ZERO Q3+2= ZERO Q3+3= ZERO ............. * * Q3-1= ZERO Q3-2= ZERO Q3-3= ZERO ............. * * * * Q40 = 0.0849 Q4+1= 0.0039 Q4+2= 0.2199 Q4+3= -0.0038 Q4+4= 0.0994 * * Q4-1= -0.0101 Q4-2=-7.9E-04 Q4-3= 7.9E-04 Q4-4=-4.1E-04 * * * ******************************************************************************* ******************************************************************************* * * * MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * M10 = 0.5491 M1+1= -0.0083 ............. ............. ............. * * M1-1= 2.7154 ............. ............. ............. * * * * M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. * * M2-1= ZERO M2-2= ZERO ............. ............. * * * * M30 = 54.9564 M3+1= -1.3218 M3+2= 3.5784 M3+3= -0.2678 ............. * * M3-1= 64.2708 M3-2= 0.3102 M3-3= 64.4256 ............. * * * ******************************************************************************* * * * CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT TOTAL * * * * M1_X = 0.0241 M1_Y = 7.8594 M1_Z = 1.1237 M1 = 7.9394 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] [INTRINSIC FRAME] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 =132.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 6.7799 Q2+1= ZERO Q2+2= 6.9289 ............. ............. * * Q2-1= ZERO Q2-2= ZERO ............. ............. * * * * Q30 = ZERO Q3+1= ZERO Q3+2= ZERO Q3+3= ZERO ............. * * Q3-1= ZERO Q3-2= ZERO Q3-3= ZERO ............. * * * * Q40 = 0.0848 Q4+1= 9.3E-09 Q4+2= 0.2200 Q4+3=-6.3E-09 Q4+4= 0.0994 * * Q4-1= -0.0103 Q4-2= 6.7E-09 Q4-3= 5.7E-04 Q4-4=-1.3E-08 * * * ******************************************************************************* ******************************************************************************* * * * MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] [INTRINSIC FRAME] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * M10 = 0.5507 M1+1= 1.4E-07 ............. ............. ............. * * M1-1= 2.7152 ............. ............. ............. * * * * M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. * * M2-1= ZERO M2-2= ZERO ............. ............. * * * * M30 = 55.0803 M3+1=-4.0E-06 M3+2= 3.6402 M3+3= 2.3E-06 ............. * * M3-1= 64.2333 M3-2= 4.5E-07 M3-3= 64.4215 ............. * * * ******************************************************************************* * * * CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT [INTRINSIC FRAME] TOTAL * * * * M1_X = 0.0000 M1_Y = 7.8590 M1_Z = 1.1271 M1 = 7.9394 * * * ******************************************************************************* ******************************************************************************* * * * ANGULAR MOMENTA IN THE THREE CARTESIAN DIRECTIONS * * * ******************************************************************************* * * * SPIN IN X-DIRECTION SPIN IN Y-DIRECTION SPIN IN Z-DIRECTION * * ------------------- ------------------- ------------------- * * INTRI ORBIT TOTAL INTRI ORBIT TOTAL INTRI ORBIT TOTAL * * * * NEUTRONS 0.006 0.101 0.107 0.073 0.518 0.591 0.428 7.209 7.637 * * PROTONS 0.002 0.035 0.037 0.392 5.949 6.341 0.131 2.030 2.160 * * -------- * * TOTAL 0.008 0.136 0.144 0.466 6.466 6.932 0.559 9.238 9.797 * * * ******************************************************************************* * * * ROTATIONAL FREQUENCIES IN THE THREE CARTESIAN DIRECTIONS: * * * * OMEGAX = 0.003009 OMEGAY = 0.144390 OMEGAZ = 0.204065 OMEGA = 0.250000 * * * ******************************************************************************* * * * ANGLE BETWEEN ANGULAR MOMENTUM AND ANGULAR FREQUENCY VECTORS= 8.5E-07 DEG * * * ******************************************************************************* * * * ANGULAR MOMENTA IN THE INTRINSIC REFERENCE FRAME * * * ******************************************************************************* * * * SPIN IN X-DIRECTION SPIN IN Y-DIRECTION SPIN IN Z-DIRECTION * * ------------------- ------------------- ------------------- * * INTRI ORBIT TOTAL INTRI ORBIT TOTAL INTRI ORBIT TOTAL * * * * NEUTRONS 0.000 0.000 0.000 0.073 0.515 0.588 0.428 7.210 7.638 * * PROTONS 0.000 0.000 0.000 0.392 5.948 6.340 0.131 2.032 2.163 * * -------- * * TOTAL 0.000 0.000 0.000 0.466 6.463 6.928 0.559 9.242 9.801 * * * ******************************************************************************* ******************************************************************************* * * * NEUTRON CONFIGURATIONS * * ======================= * * PAR 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 * * --- -------------------------------------------------------------- * * * * CONF: (+) 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 * * VACC: (+) 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 * * * * CONF: (-) 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 * * VACC: (-) 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * PROTON CONFIGURATIONS * * ======================= * * PAR 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 * * --- -------------------------------------------------------------- * * * * CONF: (+) 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 * * VACC: (+) 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 * * * * CONF: (-) 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 * * VACC: (-) 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * ENERGIES (MEV) * * * ******************************************************************************* * * * KINETIC: (NEU)= 1460.491929 (PRO)= 958.114770 (TOT)= 2418.606700 * * SUM EPS: (NEU)= -1764.973911 (PRO)= -1047.797062 (TOT)= -2812.770974 * * PAIRING: (NEU)= 0.000000 (PRO)= 0.000000 (TOT)= 0.000000 * * P-REARR: (NEU)= 0.000000 (PRO)= 0.000000 (TOT)= 0.000000 * * * * PAIRGAP: (NEU)= 0.000000 (PRO)= 0.000000 THIS SPACE IS EMPTY * * E-FERMI: (NEU)= -9.055110 (PRO)= -4.947826 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 458.000345 (EXC)= -21.974411 (TOT)= 436.025934 * * * * CONSTR. (MULT)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 * * CONSTR. (SURF)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 * * CONSTR. (SPIN)= -3.000627 SLOPE= 0.351464 CORR.= -1.500314 * * * * REARRANGEMENT ENERGY FROM THE SKYRME DENSITY-DEPENDENT TERMS= 890.885803 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -1096.793058 * * * * SPIN-ORB (EVE)= -82.387355 (ODD)= -0.138618 (TOT)= -82.525973 * * SKYRME: (EVE)= -3948.363170 (ODD)= -0.061894 (TOT)= -3948.425064 * * * * TOTAL: (STAB)= 4.76E-07 (SP)= -1093.792430 (FUN)= -1093.792431 * * * *******************************************************************************