*******************************************************************************
* *
* HFODD HFODD HFODD HFODD HFODD HFODD HFODD HFODD *
* *
*******************************************************************************
* *
* SKYRME-HARTREE-FOCK-BOGOLYUBOV CODE VERSION: 2.08K *
* *
* 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 *
* *
*******************************************************************************
*******************************************************************************
* *
* CLASSICAL NUCLEAR SURFACE DEFINED FOR: N = 75 Z = 57 *
* *
*******************************************************************************
* *
* SPHERICAL SHAPE REQUESTED FOR THE SURFACE *
* *
*******************************************************************************
* *
* HOMEGA= 9.6629 FCHOM0= 1.2000 *
* *
* OSCILLATOR FREQUENCIES: HBAROX= 9.6629 HBAROY= 9.6629 HBAROZ= 9.6629 *
* *
* MOMENTS OF INERTIA: XMOMFC= 49.9200 YMOMFC= 49.9200 ZMOMFC= 49.9200 *
* *
* CENTERS OF MASS: CMSXFC= 0.0000 CMSYFC= 0.0000 CMSZFC= 0.0000 *
* *
*******************************************************************************
*******************************************************************************
* *
* 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) *
* *
*******************************************************************************
*******************************************************************************
* *
* COEFFICIENTS DEFINING THE SKYRME FUNCTIONAL *
* *
*******************************************************************************
* *
* TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) *
* -------- ------ ------------ ------------ *
* CRHO_ = -1382.012500 780.275000 -991.875000 390.137500 *
* CRHOD = 1299.583333 -649.791667 974.687500 -324.895833 *
* CLPR_ = -85.312500 34.218750 -68.203125 17.109375 *
* CTAU_ = 68.750000 -68.125000 34.687500 -34.062500 *
* CSCU_ = 0.000000 0.000000 0.000000 0.000000 *
* CDIV_ = -65.000000 -65.000000 -97.500000 -32.500000 *
* *
* CSPI_ = -59.512500 661.250000 271.112500 330.625000 *
* CSPID = 0.000000 -649.791667 -324.895833 -324.895833 *
* CLPS_ = 0.000000 34.218750 17.109375 17.109375 *
* CCUR_ = -68.750000 68.125000 -34.687500 34.062500 *
* CKIS_ = 0.000000 0.000000 0.000000 0.000000 *
* CROT_ = -65.000000 -65.000000 -97.500000 -32.500000 *
* *
*******************************************************************************
*******************************************************************************
* *
* PHYSICAL CONSTANTS: H_BARC=197.32891000 HBCOE2=137.03602000 *
* *
* XMASSN=938.90590000 XMASSP=938.27231000 *
* *
* HBMASS= 20.73000000 HBMRPA= 20.57295455 *
* *
* ECHAR2= 1.43997841 COULEX= -1.06350868 *
* *
*******************************************************************************
* *
* OSCILLATOR LENGTHS: X= 2.0716964 Y= 2.0716964 Z= 2.0716964 *
* *
* OSCILLATOR CONSTANTS: X= 0.4826962 Y= 0.4826962 Z= 0.4826962 *
* *
* OSCILLATOR FREQUENCIES: X= 9.6628920 Y= 9.6628920 Z= 9.6628920 *
* *
*******************************************************************************
* *
* BASIS CUT-OFF CONTROL PARAMETERS: NXMAXX = 10 NYMAXX = 10 NZMAXX = 10 *
* *
* OPTIMUM NUMBERS OF GAUSS POINTS: NXHERM = 22 NYHERM = 22 NZHERM = 22 *
* *
* NLIMIT= 286 LDBASE= 286 MCOUNT= 1331 *
* *
* ENECUT= 800.0000 ELIMIT= 111.1233 *
* *
*******************************************************************************
*******************************************************************************
* *
* SHAPE OF THE OSCILLATOR-BASIS DIAMOND (THREE CONSEQUTIVE PROJECTIONS) *
* *
*******************************************************************************
* *
* NX ===>>> 0 1 2 3 4 5 6 7 8 9 10 *
* MAX.NY => 10 9 8 7 6 5 4 3 2 1 0 *
* --------------------------------- *
* NY= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 *
* NY= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 *
* NY= 2 ( 8) | 8 7 6 5 4 3 2 1 0 *
* NY= 3 ( 7) | 7 6 5 4 3 2 1 0 *
* NY= 4 ( 6) | 6 5 4 3 2 1 0 *
* NY= 5 ( 5) | 5 4 3 2 1 0 *
* NY= 6 ( 4) | 4 3 2 1 0 *
* NY= 7 ( 3) | 3 2 1 0 *
* NY= 8 ( 2) | 2 1 0 *
* NY= 9 ( 1) | 1 0 *
* NY=10 ( 0) | 0 *
* *
* *
* NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 *
* MAX.NY => 10 9 8 7 6 5 4 3 2 1 0 *
* --------------------------------- *
* NY= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 *
* NY= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 *
* NY= 2 ( 8) | 8 7 6 5 4 3 2 1 0 *
* NY= 3 ( 7) | 7 6 5 4 3 2 1 0 *
* NY= 4 ( 6) | 6 5 4 3 2 1 0 *
* NY= 5 ( 5) | 5 4 3 2 1 0 *
* NY= 6 ( 4) | 4 3 2 1 0 *
* NY= 7 ( 3) | 3 2 1 0 *
* NY= 8 ( 2) | 2 1 0 *
* NY= 9 ( 1) | 1 0 *
* NY=10 ( 0) | 0 *
* *
* *
* NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 *
* MAX.NX => 10 9 8 7 6 5 4 3 2 1 0 *
* --------------------------------- *
* NX= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 *
* NX= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 *
* NX= 2 ( 8) | 8 7 6 5 4 3 2 1 0 *
* NX= 3 ( 7) | 7 6 5 4 3 2 1 0 *
* NX= 4 ( 6) | 6 5 4 3 2 1 0 *
* NX= 5 ( 5) | 5 4 3 2 1 0 *
* NX= 6 ( 4) | 4 3 2 1 0 *
* NX= 7 ( 3) | 3 2 1 0 *
* NX= 8 ( 2) | 2 1 0 *
* NX= 9 ( 1) | 1 0 *
* NX=10 ( 0) | 0 *
* *
*******************************************************************************
*******************************************************************************
* *
* NUCLIDE: N = 75 Z = 57 NUMBER OF ITERATIONS = 5000 CONTINUATION? = 0 *
* *
*******************************************************************************
* *
* ITERATIONS WILL STOP WHEN THE STABILITY BECOMES LOWER THAN EPSITE=0.000001 *
* FOR 5 CONSECUTIVE ITERATIONS *
* *
*******************************************************************************
* *
* SLOW-DOWN PARAMETERS: TIME-EVEN = 0.50 TIME-ODD = 0.50 PAIRING = 0.50 *
* *
*******************************************************************************
* *
* COULOMB PARAMETERS: NUMCOU = 80 NUMETA = 79 BOUCOU = 20.0 FURMAX =0.25 *
* *
*******************************************************************************
* *
* MAXIMUM NUMBERS OF MULTIPOLE MOMENTS USED FOR THE COULOMB FIELD, NMUCOU=4 *
* FOR THE CONSTRAINTS, NMUCON=2 *
* FOR THE OUTPUT INFO, NMUPRI=4 *
* *
*******************************************************************************
* *
* MULTIPOLE CONSTRAINTS: LAMBDA= 2 MIU= 0 MOMENT= 7.000 STIFFNESS= 0.500 *
* LAMBDA= 2 MIU= 2 MOMENT= 7.000 STIFFNESS= 0.500 *
* *
*******************************************************************************
* *
* MAXIMUM NUMBER OF MAGNETIC MOMENTS USED: FOR THE CONSTRAINTS, NMACON=0 *
* FOR THE OUTPUT INFO, NMAPRI=3 *
* *
*******************************************************************************
* *
* PRINTING THE RESULTS FOR THE FOLLOWING ITERATIONS: THE FIRST: YES *
* THE MIDDLE: NO *
* AND/OR THE LAST: YES *
* *
*******************************************************************************
* *
* PRINTING THE MOMENTS, RADII, AND DEFORMATIONS FOR: NEUTRONS: YES *
* PROTONS : YES *
* MASS : YES *
* *
*******************************************************************************
* *
* PRINTING THE MOMENTS, MOMENTA AND DEFORMATIONS IN THE INTRINSIC FRAME: YES *
* *
*******************************************************************************
* *
* NUMBERS OF LEVELS IN TIME-REVERSAL/CHARGE BLOCKS: TIME-UP TIME-DOWN *
* NEUTRONS: 75 75 *
* PROTONS: 57 57 *
* *
*******************************************************************************
* *
* CALCULATIONS WITHOUT PAIRING *
* *
*******************************************************************************
* *
* CALCULATIONS WITH PARITY/SIGNATURE SYMMETRY *
* *
*******************************************************************************
* *
* CALCULATIONS WITH TIME-ODD SYMMETRIES: T*SIMPLEX_X => CONSERVED *
* T*SIMPLEX_Y => NON CONSERVED *
* T*SIMPLEX_Z => CONSERVED *
* *
*******************************************************************************
* *
* CALCULATIONS WITH BROKEN TIME-REVERSAL SYMMETRY *
* *
*******************************************************************************
* *
* CALCULATIONS WITH ONE-DIMENSIONAL CRANKING FOR THE ROTATIONAL FREQUENCIES: *
* *
* OMEGAY = 0.000000 (ISOSCALAR) AND OMISOY = 0.000000 (ISOVECTOR) *
* *
*******************************************************************************
* *
* PARITY/SIGNATURE CONFIGURATIONS: *
* *
* V A C U U M P A R T I C L E S H O L E S *
* =========== ================= ========= *
* (++) (+-) (-+) (--) (++) (+-) (-+) (--) (++) (+-) (-+) (--) *
* *
* NEUTRONS: 19 19 19 18 0 0 0 0 0 0 0 0 *
* PROTONS : 15 15 14 13 0 0 0 0 0 0 0 0 *
* *
*******************************************************************************
*******************************************************************************
* *
* PARAMETERS OF THE NILSSON MODEL NEUTRONS *
* *
*******************************************************************************
* *
* ESCALE =70.809 KAPPA = 0.064 MIU = 0.600 *
* *
* HOMEGA = 9.663 CNILSS =-1.231 DNILSS =-0.369 *
* *
* OSCILLATOR FREQUENCIES: HBAROX = 9.663 HBAROY = 9.663 HBAROZ = 9.663 *
* *
*******************************************************************************
*******************************************************************************
* *
* SINGLE-PARTICLE PROPERTIES: NILSSON NEUTRONS *
* *
*******************************************************************************
* *
* NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEG>
GFACT *
* *
* 79) -11.263 (21, 0, 0, 0) | 6, 1, 5,11/2> 100 1.474 0.113 0.077 *
* 80) -11.263 ( 0,21, 0, 0) | 6, 1, 5,11/2> 100 -1.474 -0.113 0.077 *
* 81) -11.263 (22, 0, 0, 0) | 6, 2, 4, 9/2> 100 0.160 0.012 0.077 *
* 82) -11.263 ( 0,22, 0, 0) | 6, 2, 4, 9/2> 100 -0.160 -0.012 0.077 *
* 83) -11.263 (23, 0, 0, 0) | 6, 2, 4, 9/2> 100 -0.809 -0.062 0.077 *
* 84) -11.263 ( 0,23, 0, 0) | 6, 2, 4, 9/2> 100 0.809 0.062 0.077 *
* 85) -11.263 (24, 0, 0, 0) | 6, 0, 6,13/2> 100 -1.051 -0.081 0.077 *
* 86) -11.263 ( 0,24, 0, 0) | 6, 0, 6,13/2> 100 1.051 0.081 0.077 *
* 87) -11.263 (25, 0, 0, 0) | 6, 3, 1, 3/2> 100 -2.447 -0.188 0.077 *
* 88) -11.263 ( 0,25, 0, 0) | 6, 3, 1, 3/2> 100 2.447 0.188 0.077 *
* 89) -11.263 (26, 0, 0, 0) | 6, 1, 5,11/2> 100 0.452 0.035 0.077 *
* 90) -11.263 ( 0,26, 0, 0) | 6, 1, 5,11/2> 100 -0.452 -0.035 0.077 *
* 91) -11.263 (27, 0, 0, 0) | 6, 2, 4, 9/2> 100 -1.278 -0.098 0.077 *
* 92) -11.263 ( 0,27, 0, 0) | 6, 2, 4, 9/2> 100 1.278 0.098 0.077 *
* 93) -11.015 (28, 0, 0, 0) | 4, 1, 1, 1/2> 100 -0.020 0.004 -0.200 *
* 94) -11.015 ( 0,28, 0, 0) | 4, 1, 1, 1/2> 100 0.020 -0.004 -0.200 *
* 95) -11.015 (29, 0, 0, 0) | 4, 0, 2, 3/2> 100 1.020 -0.204 -0.200 *
* 96) -11.015 ( 0,29, 0, 0) | 4, 0, 2, 3/2> 100 -1.020 0.204 -0.200 *
* 97) -8.739 ( 0, 0, 0,20) | 5, 1, 2, 5/2> -100 0.572 0.082 0.143 *
* 98) -8.739 ( 0, 0, 0,21) | 5, 0, 3, 7/2> -100 -1.034 -0.148 0.143 *
* 99) -8.739 ( 0, 0,20, 0) | 5, 1, 2, 5/2> -100 -0.572 -0.082 0.143 *
* 100) -8.739 ( 0, 0,21, 0) | 5, 0, 3, 7/2> -100 1.034 0.148 0.143 *
* 101) -8.739 ( 0, 0, 0,22) | 5, 4, 1, 3/2> -100 -0.185 -0.026 0.143 *
* 102) -8.739 ( 0, 0,22, 0) | 5, 4, 1, 3/2> -100 0.185 0.026 0.143 *
* 103) -8.739 ( 0, 0, 0,23) | 5, 1, 0, 1/2> -100 -1.354 -0.193 0.143 *
* 104) -8.739 ( 0, 0,23, 0) | 5, 1, 0, 1/2> -100 1.354 0.193 0.143 *
* 105) -5.047 ( 0, 0, 0,24) | 7, 2, 5,11/2> -100 -0.531 -0.035 0.067 *
* 106) -5.047 ( 0, 0,24, 0) | 7, 2, 5,11/2> -100 0.531 0.035 0.067 *
* 107) -5.047 ( 0, 0, 0,25) | 7, 3, 4, 9/2> -100 -0.599 -0.040 0.067 *
* 108) -5.047 ( 0, 0,25, 0) | 7, 3, 4, 9/2> -100 0.599 0.040 0.067 *
* 109) -5.047 ( 0, 0, 0,26) | 7, 1, 6,13/2> -100 0.778 0.052 0.067 *
* 110) -5.047 ( 0, 0,26, 0) | 7, 1, 6,13/2> -100 -0.778 -0.052 0.067 *
* 111) -5.047 ( 0, 0, 0,27) | 7, 3, 4, 9/2> -100 0.236 0.016 0.067 *
* 112) -5.047 ( 0, 0,27, 0) | 7, 3, 4, 9/2> -100 -0.236 -0.016 0.067 *
* 113) -5.047 ( 0, 0, 0,28) | 7, 1, 6,13/2> -100 2.881 0.192 0.067 *
* 114) -5.047 ( 0, 0,28, 0) | 7, 1, 6,13/2> -100 -2.881 -0.192 0.067 *
* 115) -5.047 ( 0, 0, 0,29) | 7, 4, 1, 1/2> -100 -0.308 -0.021 0.067 *
* 116) -5.047 ( 0, 0,29, 0) | 7, 4, 1, 1/2> -100 0.308 0.021 0.067 *
* 117) -5.047 ( 0, 0, 0,30) | 7, 0, 7,15/2> -100 -1.324 -0.088 0.067 *
* 118) -5.047 ( 0, 0,30, 0) | 7, 0, 7,15/2> -100 1.324 0.088 0.067 *
* 119) -5.047 ( 0, 0, 0,31) | 7, 0, 7,15/2> -100 -5.133 -0.342 0.067 *
* 120) -5.047 ( 0, 0,31, 0) | 7, 0, 7,15/2> -100 5.133 0.342 0.067 *
* 121) -4.307 ( 0, 0, 0,32) | 5, 1, 4, 7/2> -100 -2.196 0.200 -0.091 *
* 122) -4.307 ( 0, 0,32, 0) | 5, 1, 4, 7/2> -100 2.196 -0.200 -0.091 *
* 123) -4.307 ( 0, 0, 0,33) | 5, 2, 3, 5/2> -100 1.723 -0.157 -0.091 *
* 124) -4.307 ( 0, 0,33, 0) | 5, 2, 3, 5/2> -100 -1.723 0.157 -0.091 *
* 125) -4.307 ( 0, 0, 0,34) | 5, 1, 4, 7/2> -100 -0.396 0.036 -0.091 *
* 126) -4.307 ( 0, 0,34, 0) | 5, 1, 4, 7/2> -100 0.396 -0.036 -0.091 *
* 127) -4.307 ( 0, 0, 0,35) | 5, 3, 2, 3/2> -100 1.431 -0.130 -0.091 *
* 128) -4.307 ( 0, 0,35, 0) | 5, 3, 2, 3/2> -100 -1.431 0.130 -0.091 *
* 129) -4.307 ( 0, 0, 0,36) | 5, 0, 5, 9/2> -100 1.937 -0.176 -0.091 *
* 130) -4.307 ( 0, 0,36, 0) | 5, 0, 5, 9/2> -100 -1.937 0.176 -0.091 *
* 131) -2.583 ( 0, 0, 0,37) | 5, 0, 1, 3/2> -100 -1.120 -0.373 0.333 *
* 132) -2.583 ( 0, 0,37, 0) | 5, 0, 1, 3/2> -100 1.120 0.373 0.333 *
* 133) -2.583 ( 0, 0, 0,38) | 5, 0, 1, 3/2> -100 0.120 0.040 0.333 *
* 134) -2.583 ( 0, 0,38, 0) | 5, 0, 1, 3/2> -100 -0.120 -0.040 0.333 *
* 135) -0.676 (30, 0, 0, 0) | 6, 1, 3, 7/2> 100 2.331 0.259 0.111 *
* 136) -0.676 ( 0,30, 0, 0) | 6, 1, 3, 7/2> 100 -2.331 -0.259 0.111 *
* 137) -0.676 (31, 0, 0, 0) | 6, 5, 1, 3/2> 100 -0.182 -0.020 0.111 *
* 138) -0.676 ( 0,31, 0, 0) | 6, 5, 1, 3/2> 100 0.182 0.020 0.111 *
* 139) -0.676 (32, 0, 0, 0) | 6, 0, 4, 9/2> 100 -0.577 -0.064 0.111 *
* 140) -0.676 ( 0,32, 0, 0) | 6, 0, 4, 9/2> 100 0.577 0.064 0.111 *
* 141) -0.676 (33, 0, 0, 0) | 6, 0, 4, 9/2> 100 -2.330 -0.259 0.111 *
* 142) -0.676 ( 0,33, 0, 0) | 6, 0, 4, 9/2> 100 2.330 0.259 0.111 *
* 143) -0.676 (34, 0, 0, 0) | 6, 5, 1, 3/2> 100 -1.742 -0.194 0.111 *
* 144) -0.676 ( 0,34, 0, 0) | 6, 5, 1, 3/2> 100 1.742 0.194 0.111 *
* 145) -0.121 ( 0, 0, 0,39) | 5, 1, 2, 3/2> -100 -0.612 0.087 -0.143 *
* 146) -0.121 ( 0, 0,39, 0) | 5, 1, 2, 3/2> -100 0.612 -0.087 -0.143 *
* 147) -0.121 ( 0, 0, 0,40) | 5, 0, 3, 5/2> -100 1.935 -0.276 -0.143 *
* 148) -0.121 ( 0, 0,40, 0) | 5, 0, 3, 5/2> -100 -1.935 0.276 -0.143 *
* 149) -0.121 ( 0, 0, 0,41) | 5, 0, 3, 5/2> -100 0.177 -0.025 -0.143 *
* 150) -0.121 ( 0, 0,41, 0) | 5, 0, 3, 5/2> -100 -0.177 0.025 -0.143 *
* *
*******************************************************************************
*******************************************************************************
* *
* MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] NEUTRONS *
* *
*******************************************************************************
* *
* Q00 = 75.0000 ............. ............. ............. ............. *
* *
* Q10 = ZERO Q11 = ZERO ............. ............. ............. *
* *
* Q20 = -0.0311 Q21 = -0.0276 Q22 = 0.0469 ............. ............. *
* *
* Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. *
* *
* Q40 = 0.0061 Q41 = 0.0047 Q42 = 9.8E-04 Q43 = -0.0018 Q44 = -0.0018 *
* *
*******************************************************************************
*******************************************************************************
* *
* MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] NEUTRONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* M10 = ZERO M1+1= ZERO ............. ............. ............. *
* M1-1= 0.1015 ............. ............. ............. *
* *
* M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. *
* M2-1= ZERO M2-2= ZERO ............. ............. *
* *
* M30 = ZERO M3+1= ZERO M3+2= ZERO M3+3= ZERO ............. *
* M3-1= -9.5621 M3-2= 2.3043 M3-3= 2.2079 ............. *
* *
*******************************************************************************
* *
* CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT NEUTRONS *
* *
* M1_X = 0.0000 M1_Y = 0.2937 M1_Z = 0.0000 M1 = 0.2937 *
* *
*******************************************************************************
*******************************************************************************
* *
* ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS NEUTRONS *
* *
* R_RMS = 4.4972 X_RMS = 2.6012 Y_RMS = 2.5943 Z_RMS = 2.5938 *
* *
* R_GEO = 5.8058 X_GEO = 5.8165 Y_GEO = 5.8010 Z_GEO = 5.7998 *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) NEUTRONS *
* *
*******************************************************************************
* *
* B10 = ZERO B11 = ZERO ............. ............. ............. *
* *
* B20 = -0.0016 B21 = 0.0050 B22 = 0.0025 ............. ............. *
* *
* B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. *
* *
* B40 = 0.0030 B41 = 0.0033 B42 = 6.8E-04 B43 = -0.0012 B44 = -0.0013 *
* *
*******************************************************************************
*******************************************************************************
* *
* PARAMETERS OF THE NILSSON MODEL PROTONS *
* *
*******************************************************************************
* *
* ESCALE =61.146 KAPPA = 0.064 MIU = 0.600 *
* *
* HOMEGA = 9.663 CNILSS =-1.231 DNILSS =-0.369 *
* *
* OSCILLATOR FREQUENCIES: HBAROX = 9.663 HBAROY = 9.663 HBAROZ = 9.663 *
* *
*******************************************************************************
*******************************************************************************
* *
* SINGLE-PARTICLE PROPERTIES: NILSSON PROTONS *
* *
*******************************************************************************
* *
* NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEG> GFACT *
* *
* 51) -8.186 ( 0, 0, 0,14) | 5, 3, 2, 5/2> -100 3.024 0.275 0.091 *
* 52) -8.186 ( 0, 0,14, 0) | 5, 3, 2, 5/2> -100 -3.024 -0.275 0.091 *
* 53) -8.186 ( 0, 0, 0,15) | 5, 1, 4, 9/2> -100 -1.815 -0.165 0.091 *
* 54) -8.186 ( 0, 0,15, 0) | 5, 1, 4, 9/2> -100 1.815 0.165 0.091 *
* 55) -8.186 ( 0, 0, 0,16) | 5, 2, 1, 3/2> -100 -3.500 -0.318 0.091 *
* 56) -8.186 ( 0, 0,16, 0) | 5, 2, 1, 3/2> -100 3.500 0.318 0.091 *
* 57) -8.186 ( 0, 0, 0,17) | 5, 0, 5,11/2> -100 -0.161 -0.015 0.091 *
* 58) -8.186 ( 0, 0,17, 0) | 5, 0, 5,11/2> -100 0.161 0.015 0.091 *
* 59) -8.186 ( 0, 0, 0,18) | 5, 3, 0, 1/2> -100 -1.400 -0.127 0.091 *
* 60) -8.186 ( 0, 0,18, 0) | 5, 3, 0, 1/2> -100 1.400 0.127 0.091 *
* 61) -8.186 ( 0, 0, 0,19) | 5, 2, 3, 7/2> -100 0.852 0.077 0.091 *
* 62) -8.186 ( 0, 0,19, 0) | 5, 2, 3, 7/2> -100 -0.852 -0.077 0.091 *
* 63) -7.508 (13, 0, 0, 0) | 4, 0, 2, 5/2> 100 0.300 0.060 0.200 *
* 64) -7.508 ( 0,13, 0, 0) | 4, 0, 2, 5/2> 100 -0.300 -0.060 0.200 *
* 65) -7.508 (14, 0, 0, 0) | 4, 0, 2, 5/2> 100 -1.181 -0.236 0.200 *
* 66) -7.508 ( 0,14, 0, 0) | 4, 0, 2, 5/2> 100 1.181 0.236 0.200 *
* 67) -7.508 (15, 0, 0, 0) | 4, 4, 0, 1/2> 100 -0.619 -0.124 0.200 *
* 68) -7.508 ( 0,15, 0, 0) | 4, 4, 0, 1/2> 100 0.619 0.124 0.200 *
* 69) -4.061 (16, 0, 0, 0) | 4, 1, 3, 5/2> 100 -2.313 0.257 -0.111 *
* 70) -4.061 ( 0,16, 0, 0) | 4, 1, 3, 5/2> 100 2.313 -0.257 -0.111 *
* 71) -4.061 (17, 0, 0, 0) | 4, 0, 4, 7/2> 100 1.787 -0.199 -0.111 *
* 72) -4.061 ( 0,17, 0, 0) | 4, 0, 4, 7/2> 100 -1.787 0.199 -0.111 *
* 73) -4.061 (18, 0, 0, 0) | 4, 1, 3, 5/2> 100 1.218 -0.135 -0.111 *
* 74) -4.061 ( 0,18, 0, 0) | 4, 1, 3, 5/2> 100 -1.218 0.135 -0.111 *
* 75) -4.061 (19, 0, 0, 0) | 4, 0, 4, 7/2> 100 1.308 -0.145 -0.111 *
* 76) -4.061 ( 0,19, 0, 0) | 4, 0, 4, 7/2> 100 -1.308 0.145 -0.111 *
* 77) -2.830 (20, 0, 0, 0) | 4, 0, 0, 1/2> 100 -0.500 -0.500 1.000 *
* 78) -2.830 ( 0,20, 0, 0) | 4, 0, 0, 1/2> 100 0.500 0.500 1.000 *
* 79) -1.600 (21, 0, 0, 0) | 6, 1, 5,11/2> 100 0.968 0.074 0.077 *
* 80) -1.600 ( 0,21, 0, 0) | 6, 1, 5,11/2> 100 -0.968 -0.074 0.077 *
* 81) -1.600 (22, 0, 0, 0) | 6, 3, 3, 7/2> 100 -1.737 -0.134 0.077 *
* 82) -1.600 ( 0,22, 0, 0) | 6, 3, 3, 7/2> 100 1.737 0.134 0.077 *
* 83) -1.600 (23, 0, 0, 0) | 6, 2, 4, 9/2> 100 -0.153 -0.012 0.077 *
* 84) -1.600 ( 0,23, 0, 0) | 6, 2, 4, 9/2> 100 0.153 0.012 0.077 *
* 85) -1.600 (24, 0, 0, 0) | 6, 3, 1, 3/2> 100 1.376 0.106 0.077 *
* 86) -1.600 ( 0,24, 0, 0) | 6, 3, 1, 3/2> 100 -1.376 -0.106 0.077 *
* 87) -1.600 (25, 0, 0, 0) | 6, 1, 5,11/2> 100 1.364 0.105 0.077 *
* 88) -1.600 ( 0,25, 0, 0) | 6, 1, 5,11/2> 100 -1.364 -0.105 0.077 *
* 89) -1.600 (26, 0, 0, 0) | 6, 3, 3, 7/2> 100 -5.380 -0.414 0.077 *
* 90) -1.600 ( 0,26, 0, 0) | 6, 3, 3, 7/2> 100 5.380 0.414 0.077 *
* 91) -1.600 (27, 0, 0, 0) | 6, 0, 6,13/2> 100 0.062 0.005 0.077 *
* 92) -1.600 ( 0,27, 0, 0) | 6, 0, 6,13/2> 100 -0.062 -0.005 0.077 *
* 93) -1.352 (28, 0, 0, 0) | 4, 1, 1, 1/2> 100 -0.018 0.004 -0.200 *
* 94) -1.352 ( 0,28, 0, 0) | 4, 1, 1, 1/2> 100 0.018 -0.004 -0.200 *
* 95) -1.352 (29, 0, 0, 0) | 4, 0, 2, 3/2> 100 1.018 -0.204 -0.200 *
* 96) -1.352 ( 0,29, 0, 0) | 4, 0, 2, 3/2> 100 -1.018 0.204 -0.200 *
* *
*******************************************************************************
*******************************************************************************
* *
* MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] PROTONS *
* *
*******************************************************************************
* *
* Q00 = 57.0000 ............. ............. ............. ............. *
* *
* Q10 = ZERO Q11 = ZERO ............. ............. ............. *
* *
* Q20 = 0.0474 Q21 = 0.0044 Q22 = 0.0647 ............. ............. *
* *
* Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. *
* *
* Q40 = -0.0040 Q41 = -0.0019 Q42 = -0.0022 Q43 = 1.3E-04 Q44 = -0.0042 *
* *
*******************************************************************************
*******************************************************************************
* *
* MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] PROTONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* M10 = ZERO M1+1= ZERO ............. ............. ............. *
* M1-1= -1.4804 ............. ............. ............. *
* *
* M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. *
* M2-1= ZERO M2-2= ZERO ............. ............. *
* *
* M30 = ZERO M3+1= ZERO M3+2= ZERO M3+3= ZERO ............. *
* M3-1= -1.1816 M3-2=-21.7951 M3-3= 9.0873 ............. *
* *
*******************************************************************************
* *
* CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT PROTONS *
* *
* M1_X = 0.0000 M1_Y = -4.2849 M1_Z = 0.0000 M1 = 4.2849 *
* *
*******************************************************************************
*******************************************************************************
* *
* ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS PROTONS *
* *
* R_RMS = 4.2902 X_RMS = 2.4808 Y_RMS = 2.4675 Z_RMS = 2.4825 *
* *
* R_GEO = 5.5386 X_GEO = 5.5471 Y_GEO = 5.5175 Z_GEO = 5.5511 *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) PROTONS *
* *
*******************************************************************************
* *
* B10 = ZERO B11 = ZERO ............. ............. ............. *
* *
* B20 = 0.0036 B21 = -0.0012 B22 = 0.0049 ............. ............. *
* *
* B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. *
* *
* B40 = -0.0031 B41 = -0.0021 B42 = -0.0025 B43 = 1.4E-04 B44 = -0.0047 *
* *
*******************************************************************************
*******************************************************************************
* *
* DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL *
* *
*******************************************************************************
* *
* TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) *
* -------- ------ ------------ ------------ *
* DRHO_ = 19.822156 10.058382 19.822156 0.294609 *
* DRHOD = 14.803733 7.505212 14.803733 0.206691 *
* DLPR_ = -4.210954 -2.170715 -4.210954 -0.130476 *
* DTAU_ = 22.039551 11.219684 22.039551 0.399817 *
* DSCU_ = 0.352122 0.179697 0.352122 0.007272 *
* DDIV_ = 1.458345 0.747412 1.458345 0.036479 *
* *
* DSPI_ = 0.004960 0.003588 0.004960 0.002217 *
* DSPID = 0.003401 0.002461 0.003401 0.001520 *
* DLPS_ = -0.017821 -0.013061 -0.017821 -0.008300 *
* DCUR_ = 0.001527 0.001059 0.001527 0.000591 *
* DKIS_ = -0.000529 -0.000466 -0.000529 -0.000403 *
* DROT_ = 0.000673 0.000493 0.000673 0.000313 *
* *
*******************************************************************************
*******************************************************************************
* *
* CONTRIBUTIONS TO ENERGY IN THE SKYRME FUNCTIONAL *
* *
*******************************************************************************
* *
* TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) *
* -------- ------ ------------ ------------ *
* ERHO_ = -27394.467363 7848.304309 -19661.100978 114.937924 *
* ERHOD = 19238.684061 -4876.824056 14429.013046 -67.153040 *
* ELPR_ = 359.247036 -74.279165 287.200241 -2.232370 *
* ETAU_ = 1515.219098 -764.340950 764.496909 -13.618761 *
* ESCU_ = 0.000000 0.000000 0.000000 0.000000 *
* EDIV_ = -94.792426 -48.581765 -142.188639 -1.185552 *
* ============ ============ ============ ============ *
* SUM EVEN: -6376.109593 2084.278373 -4322.579421 30.748202 *
* *
* ESPI_ = -0.295158 2.372791 1.344607 0.733026 *
* ESPID = 0.000000 -1.598974 -1.105088 -0.493886 *
* ELPS_ = 0.000000 -0.446923 -0.304909 -0.142015 *
* ECUR_ = -0.104948 0.072134 -0.052951 0.020137 *
* EKIS_ = 0.000000 0.000000 0.000000 0.000000 *
* EROT_ = -0.043761 -0.032044 -0.065641 -0.010163 *
* ============ ============ ============ ============ *
* SUM ODD: -0.443866 0.366984 -0.183982 0.107099 *
* *
*******************************************************************************
*******************************************************************************
* *
* EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL *
* *
* ALPHA =180.00000 BETA = 58.75884 GAMMA =180.00000 *
* *
*******************************************************************************
*******************************************************************************
* *
* ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL *
* *
* R_RMS = 4.4090 X_RMS = 2.5499 Y_RMS = 2.5403 Z_RMS = 2.5463 *
* *
* R_GEO = 5.6920 X_GEO = 5.7018 Y_GEO = 5.6803 Z_GEO = 5.6938 *
* *
*******************************************************************************
*******************************************************************************
* *
* MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL *
* *
*******************************************************************************
* *
* Q00 =132.0000 ............. ............. ............. ............. *
* *
* Q10 = ZERO Q11 = ZERO ............. ............. ............. *
* *
* Q20 = 0.0163 Q21 = -0.0231 Q22 = 0.1117 ............. ............. *
* *
* Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. *
* *
* Q40 = 0.0021 Q41 = 0.0028 Q42 = -0.0013 Q43 = -0.0017 Q44 = -0.0060 *
* *
*******************************************************************************
*******************************************************************************
* *
* MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] TOTAL *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* M10 = ZERO M1+1= ZERO ............. ............. ............. *
* M1-1= -1.3789 ............. ............. ............. *
* *
* M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. *
* M2-1= ZERO M2-2= ZERO ............. ............. *
* *
* M30 = ZERO M3+1= ZERO M3+2= ZERO M3+3= ZERO ............. *
* M3-1=-10.7438 M3-2=-19.4909 M3-3= 11.2952 ............. *
* *
*******************************************************************************
* *
* CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT TOTAL *
* *
* M1_X = 0.0000 M1_Y = -3.9912 M1_Z = 0.0000 M1 = 3.9912 *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) TOTAL *
* *
*******************************************************************************
* *
* B10 = ZERO B11 = ZERO ............. ............. ............. *
* *
* B20 = 5.0E-04 B21 = 0.0025 B22 = 0.0034 ............. ............. *
* *
* B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. *
* *
* B40 = 6.2E-04 B41 = 0.0012 B42 =-5.4E-04 B43 =-7.1E-04 B44 = -0.0026 *
* *
*******************************************************************************
*******************************************************************************
* *
* 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 = 0.1306 Q2+1= ZERO Q2+2= 0.0457 ............. ............. *
* 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.0078 Q4+1= -0.0021 Q4+2= -0.0016 Q4+3= -0.0029 Q4+4= -0.0022 *
* Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) [INTRINSIC FRAME] TOTAL *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* B10 = ZERO B1+1= ZERO ............. ............. ............. *
* B1-1= ZERO ............. ............. ............. *
* *
* B20 = 0.0040 B2+1= ZERO B2+2= 0.0014 ............. ............. *
* B2-1= ZERO B2-2= ZERO ............. ............. *
* *
* B30 = ZERO B3+1= ZERO B3+2= ZERO B3+3= ZERO ............. *
* B3-1= ZERO B3-2= ZERO B3-3= ZERO ............. *
* *
* B40 = -0.0024 B4+1=-8.8E-04 B4+2=-6.8E-04 B4+3= -0.0012 B4+4=-9.6E-04 *
* B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO *
* *
*******************************************************************************
*******************************************************************************
* *
* MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] [INTRINSIC FRAME] TOTAL *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* M10 = ZERO M1+1= ZERO ............. ............. ............. *
* M1-1= -1.3789 ............. ............. ............. *
* *
* M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. *
* M2-1= ZERO M2-2= ZERO ............. ............. *
* *
* M30 = ZERO M3+1= ZERO M3+2= ZERO M3+3= ZERO ............. *
* M3-1= 20.7339 M3-2= -4.6619 M3-3=-13.0873 ............. *
* *
*******************************************************************************
* *
* CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT [INTRINSIC FRAME] TOTAL *
* *
* M1_X = 0.0000 M1_Y = -3.9912 M1_Z = 0.0000 M1 = 3.9912 *
* *
*******************************************************************************
*******************************************************************************
* *
* ANGULAR MOMENTA AND THE FIRST MOMENTS OF INERTIA *
* *
*******************************************************************************
* *
* SPINS J(1) *
* --------------------------- --------------------------- *
* ORBITAL INTRINSIC TOTAL ORBITAL INTRINSIC TOTAL *
* *
* NEUTRONS -0.76757 -0.07676 -0.84433 *
* PROTONS -2.74921 -0.27492 -3.02413 *
* -------- *
* TOTAL -3.51678 -0.35168 -3.86846 *
* *
*******************************************************************************
* *
* 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.077 -0.768 -0.844 0.000 0.000 0.000 *
* PROTONS 0.000 0.000 0.000 -0.275 -2.749 -3.024 0.000 0.000 0.000 *
* -------- *
* TOTAL 0.000 0.000 0.000 -0.352 -3.517 -3.868 0.000 0.000 0.000 *
* *
*******************************************************************************
*******************************************************************************
* *
* NEUTRON CONFIGURATIONS *
* ======================= *
* P S 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 *
* --- -------------------------------------------------------------- *
* *
* CONF: + + 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 *
* VACC: + + * * * * * * * * * * * * * * * * * * * * * *
* *
* CONF: + - 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 *
* VACC: + - * * * * * * * * * * * * * * * * * * * * * *
* *
* CONF: - + 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 *
* VACC: - + * * * * * * * * * * * * * * * * * * * * * *
* *
* CONF: - - 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 *
* VACC: - - * * * * * * * * * * * * * * * * * * * * * *
* *
*******************************************************************************
* *
* PROTON CONFIGURATIONS *
* ======================= *
* P S 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 *
* --- -------------------------------------------------------------- *
* *
* CONF: + + 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 *
* VACC: + + * * * * * * * * * * * * * * * * * * * * * *
* *
* CONF: + - 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 *
* VACC: + - * * * * * * * * * * * * * * * * * * * * * *
* *
* CONF: - + 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 *
* VACC: - + * * * * * * * * * * * * * * * * * * * * * *
* *
* CONF: - - 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 *
* VACC: - - * * * * * * * * * * * * * * * * * * * * * *
* *
*******************************************************************************
* *
* ENERGIES (MEV) *
* *
*******************************************************************************
* *
* KINETIC: (NEU)= 1694.879561 (PRO)= 1172.274918 (TOT)= 2867.154480 *
* SUM EPS: (NEU)= -2040.091687 (PRO)= -1201.034464 (TOT)= -3241.126151 *
* 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)= -13.723524 (PRO)= -8.185527 AND AWAITS YOUR AD *
* *
* COULOMB: (DIR)= 514.411356 (EXC)= -24.076053 (TOT)= 490.335303 *
* *
* CONSTR. (MULT)= 48.110727 SLOPE= 13.872052 CORR.= -0.441453 *
* CONSTR. (SURF)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 *
* CONSTR. (SPIN)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 *
* *
* REARRANGEMENT ENERGY FROM THE SKYRME DENSITY-DEPENDENT TERMS= 1196.688419 *
* ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -886.307592 *
* *
* SPIN-ORB (EVE)= -143.374191 (ODD)= -0.075804 (TOT)= -143.449995 *
* SKYRME: (EVE)= -4291.831220 (ODD)= -0.076882 (TOT)= -4291.908102 *
* *
* TOTAL: (STAB)= -456.839833 (SP)= -1391.258152 (FUN)= -934.418319 *
* *
*******************************************************************************
*******************************************************************************
* *
* CONVERGENCE REPORT *
* *
*******************************************************************************
* *
* ITER ENERGY STABILITY Q20 GAMMA SPIN OMEGA ANGLE EPAIR *
* *
* 0 -934.418319 -456.839833 0.131 19.27 -3.87 0.000 0.000 0.00 *
* 1 -1051.109839 113.040704 4.257 37.89 5.30 0.000 0.000 0.00 *
* 2 -1080.306027 71.300163 3.567 45.29 8.14 0.000 0.000 0.00 *
* 3 -1087.837314 36.906723 5.073 40.82 5.50 0.000 0.000 0.00 *
* 4 -1090.707766 21.676060 5.264 43.47 5.58 0.000 0.000 0.00 *
* 5 -1092.361964 11.034720 5.785 43.12 5.62 0.000 0.000 0.00 *
* 6 -1093.245287 6.084546 6.030 43.76 5.62 0.000 0.000 0.00 *
* 7 -1093.767575 3.056759 6.267 43.92 5.60 0.000 0.000 0.00 *
* 8 -1094.067771 1.508854 6.423 44.19 5.59 0.000 0.000 0.00 *
* 9 -1094.243339 0.615325 6.549 44.37 5.58 0.000 0.000 0.00 *
* 10 -1094.346290 0.154723 6.640 44.53 5.58 0.000 0.000 0.00 *
* 11 -1094.407045 -0.085736 6.710 44.66 5.57 0.000 0.000 0.00 *
* 12 -1094.443088 -0.193763 6.763 44.77 5.57 0.000 0.000 0.00 *
* 13 -1094.464542 -0.232657 6.803 44.86 5.56 0.000 0.000 0.00 *
* 14 -1094.477349 -0.233703 6.834 44.93 5.56 0.000 0.000 0.00 *
* 15 -1094.485001 -0.216634 6.857 44.99 5.56 0.000 0.000 0.00 *
* 16 -1094.489571 -0.191761 6.875 45.04 5.56 0.000 0.000 0.00 *
* 17 -1094.492294 -0.164985 6.888 45.07 5.56 0.000 0.000 0.00 *
* 18 -1094.493908 -0.139250 6.899 45.10 5.56 0.000 0.000 0.00 *
* 19 -1094.494857 -0.115975 6.907 45.13 5.56 0.000 0.000 0.00 *
* 20 -1094.495408 -0.095664 6.914 45.14 5.56 0.000 0.000 0.00 *
* 21 -1094.495721 -0.078356 6.919 45.16 5.55 0.000 0.000 0.00 *
* 22 -1094.495894 -0.063842 6.923 45.17 5.55 0.000 0.000 0.00 *
* 23 -1094.495984 -0.051812 6.926 45.18 5.55 0.000 0.000 0.00 *
* 24 -1094.496026 -0.041926 6.928 45.19 5.55 0.000 0.000 0.00 *
* 25 -1094.496040 -0.033853 6.930 45.19 5.55 0.000 0.000 0.00 *
* 26 -1094.496040 -0.027293 6.931 45.19 5.55 0.000 0.000 0.00 *
* 27 -1094.496032 -0.021982 6.933 45.20 5.55 0.000 0.000 0.00 *
* 28 -1094.496021 -0.017694 6.934 45.20 5.55 0.000 0.000 0.00 *
* 29 -1094.496009 -0.014240 6.934 45.20 5.55 0.000 0.000 0.00 *
* 30 -1094.495996 -0.011461 6.935 45.20 5.55 0.000 0.000 0.00 *
* 31 -1094.495985 -0.009229 6.935 45.20 5.55 0.000 0.000 0.00 *
* 32 -1094.495974 -0.007436 6.936 45.20 5.55 0.000 0.000 0.00 *
* 33 -1094.495965 -0.005998 6.936 45.20 5.55 0.000 0.000 0.00 *
* 34 -1094.495956 -0.004844 6.937 45.20 5.55 0.000 0.000 0.00 *
* 35 -1094.495949 -0.003917 6.937 45.20 5.55 0.000 0.000 0.00 *
* 36 -1094.495942 -0.003173 6.937 45.20 5.55 0.000 0.000 0.00 *
* 37 -1094.495937 -0.002576 6.937 45.20 5.55 0.000 0.000 0.00 *
* 38 -1094.495932 -0.002095 6.937 45.20 5.55 0.000 0.000 0.00 *
* 39 -1094.495927 -0.001708 6.937 45.20 5.55 0.000 0.000 0.00 *
* 40 -1094.495924 -0.001396 6.937 45.20 5.55 0.000 0.000 0.00 *
* 41 -1094.495920 -0.001144 6.938 45.20 5.55 0.000 0.000 0.00 *
* 42 -1094.495917 -0.000940 6.938 45.20 5.55 0.000 0.000 0.00 *
* 43 -1094.495915 -0.000774 6.938 45.20 5.55 0.000 0.000 0.00 *
* 44 -1094.495913 -0.000640 6.938 45.20 5.55 0.000 0.000 0.00 *
* 45 -1094.495911 -0.000531 6.938 45.20 5.55 0.000 0.000 0.00 *
* 46 -1094.495909 -0.000442 6.938 45.20 5.55 0.000 0.000 0.00 *
* 47 -1094.495908 -0.000369 6.938 45.20 5.55 0.000 0.000 0.00 *
* 48 -1094.495906 -0.000309 6.938 45.20 5.55 0.000 0.000 0.00 *
* 49 -1094.495905 -0.000259 6.938 45.20 5.55 0.000 0.000 0.00 *
* 50 -1094.495904 -0.000219 6.938 45.20 5.55 0.000 0.000 0.00 *
* 51 -1094.495903 -0.000185 6.938 45.20 5.55 0.000 0.000 0.00 *
* 52 -1094.495902 -0.000157 6.938 45.20 5.55 0.000 0.000 0.00 *
* 53 -1094.495902 -0.000133 6.938 45.20 5.55 0.000 0.000 0.00 *
* 54 -1094.495901 -0.000114 6.938 45.20 5.55 0.000 0.000 0.00 *
* 55 -1094.495901 -0.000097 6.938 45.20 5.55 0.000 0.000 0.00 *
* 56 -1094.495900 -0.000083 6.938 45.20 5.55 0.000 0.000 0.00 *
* 57 -1094.495900 -0.000072 6.938 45.20 5.55 0.000 0.000 0.00 *
* 58 -1094.495899 -0.000062 6.938 45.20 5.55 0.000 0.000 0.00 *
* 59 -1094.495899 -0.000053 6.938 45.20 5.55 0.000 0.000 0.00 *
* 60 -1094.495899 -0.000046 6.938 45.20 5.55 0.000 0.000 0.00 *
* 61 -1094.495898 -0.000040 6.938 45.20 5.55 0.000 0.000 0.00 *
* 62 -1094.495898 -0.000035 6.938 45.20 5.55 0.000 0.000 0.00 *
* 63 -1094.495898 -0.000030 6.938 45.20 5.55 0.000 0.000 0.00 *
* 64 -1094.495898 -0.000026 6.938 45.20 5.55 0.000 0.000 0.00 *
* 65 -1094.495898 -0.000023 6.938 45.20 5.55 0.000 0.000 0.00 *
* 66 -1094.495898 -0.000020 6.938 45.20 5.55 0.000 0.000 0.00 *
* 67 -1094.495897 -0.000018 6.938 45.20 5.55 0.000 0.000 0.00 *
* 68 -1094.495897 -0.000015 6.938 45.20 5.55 0.000 0.000 0.00 *
* 69 -1094.495897 -0.000014 6.938 45.20 5.55 0.000 0.000 0.00 *
* 70 -1094.495897 -0.000012 6.938 45.20 5.55 0.000 0.000 0.00 *
* 71 -1094.495897 -0.000010 6.938 45.20 5.55 0.000 0.000 0.00 *
* 72 -1094.495897 -9.17E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 73 -1094.495897 -8.07E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 74 -1094.495897 -7.09E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 75 -1094.495897 -6.24E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 76 -1094.495897 -5.50E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 77 -1094.495897 -4.84E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 78 -1094.495897 -4.27E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 79 -1094.495897 -3.76E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 80 -1094.495897 -3.32E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 81 -1094.495897 -2.93E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 82 -1094.495897 -2.58E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 83 -1094.495897 -2.28E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 84 -1094.495897 -2.01E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 85 -1094.495897 -1.78E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 86 -1094.495897 -1.57E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 87 -1094.495897 -1.38E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 88 -1094.495897 -1.22E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 89 -1094.495897 -1.08E-06 6.938 45.20 5.55 0.000 0.000 0.00 *
* 90 -1094.495897 -9.54E-07 6.938 45.20 5.55 0.000 0.000 0.00 *
* 91 -1094.495897 -8.43E-07 6.938 45.20 5.55 0.000 0.000 0.00 *
* 92 -1094.495897 -7.44E-07 6.938 45.20 5.55 0.000 0.000 0.00 *
* 93 -1094.495897 -6.58E-07 6.938 45.20 5.55 0.000 0.000 0.00 *
* 94 -1094.495897 -5.81E-07 6.938 45.20 5.55 0.000 0.000 0.00 *
* *
*******************************************************************************
*******************************************************************************
* *
* SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS *
* *
*******************************************************************************
* *
* NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEG> GFACT *
* *
* 65) -11.694 ( 0,18, 0, 0) | 4, 1, 3, 5/2> 100 1.393 0.005 0.003 *
* 66) -11.653 (18, 0, 0, 0) | 4, 1, 3, 5/2> 100 -1.425 -0.014 0.010 *
* 67) -11.272 ( 0, 0,16, 0) | 5, 2, 3, 7/2> -100 3.111 0.259 0.083 *
* 68) -11.135 ( 0, 0, 0,16) | 5, 2, 3, 7/2> -100 -3.117 -0.251 0.080 *
* 69) -10.692 ( 0,19, 0, 0) | 4, 0, 4, 7/2> 100 0.698 -0.132 -0.188 *
* 70) -10.665 (19, 0, 0, 0) | 4, 0, 4, 7/2> 100 -0.769 0.125 -0.162 *
* 71) -10.598 ( 0, 0, 0,17) | 5, 2, 3, 7/2> -100 1.542 0.064 0.041 *
* 72) -10.525 ( 0, 0,17, 0) | 5, 2, 3, 7/2> -100 -1.222 -0.037 0.030 *
* 73) -9.913 ( 0, 0, 0,18) | 5, 1, 4, 9/2> -100 -0.282 -0.044 0.158 *
* 74) -9.841 ( 0, 0,18, 0) | 5, 1, 4, 9/2> -100 0.053 0.033 0.625 *
* 75) -8.687 (20, 0, 0, 0) | 4, 1, 1, 1/2> 100 -0.010 0.130 -1.E+01 *
* 76) -8.684 ( 0,20, 0, 0) | 4, 1, 1, 1/2> 100 0.004 -0.137 -3.E+01 *
* 77) -8.523 ( 0, 0,19, 0) | 5, 0, 5,11/2> -100 -0.098 0.010 -0.100 *
* 78) -8.483 ( 0, 0, 0,19) | 5, 0, 5,11/2> -100 -0.033 -0.018 0.541 *
* 79) -8.039 (21, 0, 0, 0) | 4, 0, 4, 7/2> 100 0.985 0.093 0.094 *
* 80) -8.028 ( 0,21, 0, 0) | 4, 0, 4, 7/2> 100 -0.981 -0.100 0.102 *
* 81) -7.576 ( 0, 0,20, 0) | 5, 2, 1, 1/2> -100 3.228 0.336 0.104 *
* 82) -7.472 ( 0, 0, 0,20) | 5, 2, 1, 1/2> -100 -3.220 -0.335 0.104 *
* 83) -6.230 ( 0, 0,21, 0) | 5, 5, 0, 1/2> -100 0.055 0.360 6.545 *
* 84) -6.198 ( 0, 0, 0,21) | 5, 5, 0, 1/2> -100 0.625 -0.347 -0.555 *
* 85) -6.130 ( 0,22, 0, 0) | 6, 4, 0, 1/2> 100 6.425 0.499 0.078 *
* 86) -5.856 ( 0, 0, 0,22) | 5, 3, 0, 1/2> -100 2.271 -0.297 -0.131 *
* 87) -5.742 ( 0, 0,22, 0) | 5, 3, 2, 3/2> -100 -2.280 0.302 -0.132 *
* 88) -5.569 (22, 0, 0, 0) | 4, 0, 2, 3/2> 100 -0.495 -0.302 0.611 *
* 89) -5.568 ( 0,23, 0, 0) | 4, 0, 2, 3/2> 100 0.494 0.299 0.606 *
* 90) -5.447 (23, 0, 0, 0) | 6, 4, 0, 1/2> 100 -6.425 -0.499 0.078 *
* 91) -5.079 ( 0, 0, 0,23) | 5, 1, 0, 1/2> -100 1.997 -0.320 -0.160 *
* 92) -5.002 ( 0, 0,23, 0) | 5, 2, 1, 3/2> -100 -2.456 0.300 -0.122 *
* 93) -3.968 ( 0, 0,24, 0) | 5, 0, 1, 3/2> -100 -1.051 0.247 -0.235 *
* 94) -3.965 (24, 0, 0, 0) | 6, 4, 2, 5/2> 100 5.361 0.301 0.056 *
* 95) -3.930 ( 0, 0, 0,24) | 5, 0, 1, 3/2> -100 0.945 -0.221 -0.234 *
* 96) -3.637 ( 0,24, 0, 0) | 6, 4, 2, 5/2> 100 -5.322 -0.297 0.056 *
* 97) -3.535 ( 0, 0, 0,25) | 5, 1, 2, 3/2> -100 1.709 0.072 0.042 *
* 98) -3.498 ( 0, 0,25, 0) | 5, 1, 2, 3/2> -100 -1.738 -0.095 0.055 *
* 99) -2.870 ( 0, 0,26, 0) | 5, 3, 2, 3/2> -100 0.672 -0.210 -0.313 *
* 100) -2.860 ( 0, 0, 0,26) | 5, 3, 2, 3/2> -100 -0.561 0.218 -0.389 *
* 101) -2.816 ( 0,25, 0, 0) | 6, 3, 3, 7/2> 100 4.183 0.283 0.068 *
* 102) -2.675 (25, 0, 0, 0) | 6, 3, 3, 7/2> 100 -4.191 -0.279 0.067 *
* 103) -2.113 ( 0, 0,27, 0) | 5, 2, 3, 5/2> -100 1.669 -0.169 -0.101 *
* 104) -2.068 ( 0, 0, 0,27) | 5, 2, 3, 5/2> -100 -1.640 0.168 -0.103 *
* 105) -1.978 (26, 0, 0, 0) | 6, 5, 1, 3/2> 100 2.887 0.141 0.049 *
* 106) -1.856 ( 0,26, 0, 0) | 6, 5, 1, 3/2> 100 -2.678 -0.127 0.048 *
* 107) -1.481 ( 0,27, 0, 0) | 6, 2, 4, 9/2> 100 1.153 0.086 0.075 *
* 108) -1.465 (27, 0, 0, 0) | 6, 2, 4, 9/2> 100 -1.116 -0.081 0.073 *
* 109) -1.367 ( 0, 0,28, 0) | 5, 0, 3, 5/2> -100 1.368 -0.279 -0.204 *
* 110) -1.343 ( 0, 0, 0,28) | 5, 0, 3, 7/2> -100 -1.584 0.272 -0.171 *
* 111) -0.662 ( 0,28, 0, 0) | 6, 1, 5,11/2> 100 -0.238 0.002 -0.010 *
* 112) -0.639 (28, 0, 0, 0) | 6, 1, 5,11/2> 100 0.068 -0.013 -0.188 *
* 113) -0.221 ( 0, 0, 0,29) | 5, 1, 0, 1/2> -100 -0.690 -0.106 0.153 *
* 114) -0.208 ( 0, 0,29, 0) | 5, 1, 0, 1/2> -100 0.735 0.103 0.140 *
* *
*******************************************************************************
*******************************************************************************
* *
* MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] NEUTRONS *
* *
*******************************************************************************
* *
* Q00 = 75.0000 ............. ............. ............. ............. *
* *
* Q10 = ZERO Q11 = ZERO ............. ............. ............. *
* *
* Q20 = 3.8497 Q21 = ZERO Q22 = 4.0563 ............. ............. *
* *
* Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. *
* *
* Q40 = 0.0442 Q41 = ZERO Q42 = 0.1208 Q43 = ZERO Q44 = 0.0524 *
* *
*******************************************************************************
*******************************************************************************
* *
* MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] NEUTRONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* M10 = ZERO M1+1= ZERO ............. ............. ............. *
* M1-1= -0.0837 ............. ............. ............. *
* *
* M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. *
* M2-1= ZERO M2-2= ZERO ............. ............. *
* *
* M30 = ZERO M3+1= ZERO M3+2= ZERO M3+3= ZERO ............. *
* M3-1= -3.0562 M3-2= ZERO M3-3= -6.3408 ............. *
* *
*******************************************************************************
* *
* CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT NEUTRONS *
* *
* M1_X = 0.0000 M1_Y = -0.2423 M1_Z = 0.0000 M1 = 0.2423 *
* *
*******************************************************************************
*******************************************************************************
* *
* ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS NEUTRONS *
* *
* R_RMS = 4.9005 X_RMS = 2.9514 Y_RMS = 2.3640 Z_RMS = 3.1171 *
* *
* R_GEO = 6.3266 X_GEO = 6.5996 Y_GEO = 5.2860 Z_GEO = 6.9700 *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) NEUTRONS *
* *
*******************************************************************************
* *
* B10 = ZERO B11 = ZERO ............. ............. ............. *
* *
* B20 = 0.1694 B21 = ZERO B22 = 0.1785 ............. ............. *
* *
* B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. *
* *
* B40 = 0.0154 B41 = ZERO B42 = 0.0596 B43 = ZERO B44 = 0.0259 *
* *
*******************************************************************************
*******************************************************************************
* *
* SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS *
* *
*******************************************************************************
* *
* NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEG> GFACT *
* *
* 37) -11.827 ( 0, 0, 0,11) | 3, 1, 0, 1/2> -100 -1.033 -0.331 0.320 *
* 38) -11.818 ( 0, 9, 0, 0) | 4, 2, 2, 5/2> 100 -3.408 -0.177 0.052 *
* 39) -11.794 ( 0, 0,11, 0) | 3, 1, 0, 1/2> -100 1.014 0.347 0.343 *
* 40) -11.707 ( 9, 0, 0, 0) | 4, 2, 2, 5/2> 100 3.393 0.117 0.035 *
* 41) -11.182 ( 0, 0,12, 0) | 3, 0, 1, 3/2> -100 -0.019 0.310 -2.E+01 *
* 42) -11.174 ( 0, 0, 0,12) | 3, 0, 3, 5/2> -100 0.040 -0.293 -7.358 *
* 43) -10.880 (10, 0, 0, 0) | 4, 2, 2, 5/2> 100 -2.221 -0.227 0.102 *
* 44) -10.691 ( 0,10, 0, 0) | 4, 2, 2, 5/2> 100 2.193 0.176 0.080 *
* 45) -10.103 (11, 0, 0, 0) | 4, 1, 3, 7/2> 100 0.785 -0.043 -0.055 *
* 46) -10.065 ( 0,11, 0, 0) | 4, 1, 3, 7/2> 100 -0.829 -0.019 0.023 *
* 47) -8.975 (12, 0, 0, 0) | 4, 0, 4, 9/2> 100 0.014 -0.053 -3.775 *
* 48) -8.930 ( 0,12, 0, 0) | 4, 0, 4, 9/2> 100 0.070 0.012 0.170 *
* 49) -7.131 ( 0, 0,13, 0) | 3, 0, 1, 1/2> -100 -0.499 -0.417 0.836 *
* 50) -7.120 ( 0, 0, 0,13) | 3, 0, 1, 1/2> -100 0.499 0.420 0.842 *
* 51) -6.998 (13, 0, 0, 0) | 4, 2, 2, 3/2> 100 3.260 -0.121 -0.037 *
* 52) -6.703 ( 0,13, 0, 0) | 4, 2, 2, 3/2> 100 -3.225 0.180 -0.056 *
* 53) -6.689 (14, 0, 0, 0) | 4, 1, 1, 1/2> 100 -2.270 -0.170 0.075 *
* 54) -6.422 ( 0,14, 0, 0) | 4, 4, 0, 1/2> 100 2.235 0.186 0.083 *
* 55) -5.553 ( 0, 0, 0,14) | 5, 3, 0, 1/2> -100 -5.450 -0.499 0.092 *
* 56) -5.512 (15, 0, 0, 0) | 4, 4, 0, 1/2> 100 0.621 -0.193 -0.311 *
* 57) -5.225 ( 0,15, 0, 0) | 4, 4, 0, 1/2> 100 -0.556 0.229 -0.412 *
* 58) -5.052 ( 0, 0,14, 0) | 5, 3, 0, 1/2> -100 5.454 0.499 0.092 *
* 59) -4.291 (16, 0, 0, 0) | 4, 0, 0, 1/2> 100 0.346 -0.276 -0.796 *
* 60) -4.013 ( 0,16, 0, 0) | 4, 0, 0, 1/2> 100 -0.414 0.301 -0.727 *
* 61) -3.618 ( 0, 0,15, 0) | 5, 3, 2, 5/2> -100 -4.423 -0.269 0.061 *
* 62) -3.435 ( 0, 0, 0,15) | 5, 3, 2, 5/2> -100 4.413 0.223 0.050 *
* 63) -3.243 ( 0,17, 0, 0) | 4, 2, 2, 3/2> 100 0.772 0.112 0.146 *
* 64) -3.203 (17, 0, 0, 0) | 4, 2, 2, 3/2> 100 -1.120 -0.090 0.080 *
* 65) -3.075 ( 0,18, 0, 0) | 4, 1, 3, 5/2> 100 0.257 0.094 0.368 *
* 66) -3.027 (18, 0, 0, 0) | 4, 1, 3, 5/2> 100 0.097 -0.072 -0.749 *
* 67) -2.556 ( 0, 0, 0,16) | 5, 2, 3, 7/2> -100 -3.311 -0.258 0.078 *
* 68) -2.354 ( 0, 0,16, 0) | 5, 2, 3, 7/2> -100 3.310 0.218 0.066 *
* 69) -2.027 ( 0,19, 0, 0) | 4, 0, 4, 7/2> 100 0.444 -0.005 -0.010 *
* 70) -1.953 (19, 0, 0, 0) | 4, 0, 4, 7/2> 100 -0.447 0.023 -0.052 *
* 71) -1.608 ( 0, 0,17, 0) | 5, 2, 3, 7/2> -100 -1.952 -0.105 0.054 *
* 72) -1.575 ( 0, 0, 0,17) | 5, 2, 3, 7/2> -100 1.766 0.043 0.024 *
* 73) -1.004 ( 0, 0, 0,18) | 5, 1, 4, 9/2> -100 -0.446 -0.073 0.163 *
* 74) -0.977 ( 0, 0,18, 0) | 5, 1, 4, 9/2> -100 0.572 0.036 0.064 *
* *
*******************************************************************************
*******************************************************************************
* *
* MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] PROTONS *
* *
*******************************************************************************
* *
* Q00 = 57.0000 ............. ............. ............. ............. *
* *
* Q10 = ZERO Q11 = ZERO ............. ............. ............. *
* *
* Q20 = 3.0884 Q21 = ZERO Q22 = 2.9301 ............. ............. *
* *
* Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. *
* *
* Q40 = 0.0440 Q41 = ZERO Q42 = 0.1061 Q43 = ZERO Q44 = 0.0462 *
* *
*******************************************************************************
*******************************************************************************
* *
* MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] PROTONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* M10 = ZERO M1+1= ZERO ............. ............. ............. *
* M1-1= 2.5350 ............. ............. ............. *
* *
* M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. *
* M2-1= ZERO M2-2= ZERO ............. ............. *
* *
* M30 = ZERO M3+1= ZERO M3+2= ZERO M3+3= ZERO ............. *
* M3-1= 74.8217 M3-2= ZERO M3-3= 62.4605 ............. *
* *
*******************************************************************************
* *
* CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT PROTONS *
* *
* M1_X = 0.0000 M1_Y = 7.3373 M1_Z = 0.0000 M1 = 7.3373 *
* *
*******************************************************************************
*******************************************************************************
* *
* ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS PROTONS *
* *
* R_RMS = 4.8195 X_RMS = 2.8850 Y_RMS = 2.3142 Z_RMS = 3.0901 *
* *
* R_GEO = 6.2219 X_GEO = 6.4511 Y_GEO = 5.1747 Z_GEO = 6.9096 *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) PROTONS *
* *
*******************************************************************************
* *
* B10 = ZERO B11 = ZERO ............. ............. ............. *
* *
* B20 = 0.1849 B21 = ZERO B22 = 0.1754 ............. ............. *
* *
* B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. *
* *
* B40 = 0.0216 B41 = ZERO B42 = 0.0736 B43 = ZERO B44 = 0.0320 *
* *
*******************************************************************************
*******************************************************************************
* *
* DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL *
* *
*******************************************************************************
* *
* TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) *
* -------- ------ ------------ ------------ *
* DRHO_ = 15.469618 7.855652 15.469618 0.241686 *
* DRHOD = 11.024977 5.597076 11.024977 0.169175 *
* DLPR_ = -3.283688 -1.667988 -3.283688 -0.052288 *
* DTAU_ = 13.887556 7.118566 13.887556 0.349576 *
* DSCU_ = 0.141200 0.082368 0.141200 0.023536 *
* DDIV_ = 0.833557 0.431065 0.833557 0.028572 *
* *
* DSPI_ = 0.003378 0.002684 0.003378 0.001991 *
* DSPID = 0.002313 0.001830 0.002313 0.001346 *
* DLPS_ = -0.010525 -0.008794 -0.010525 -0.007062 *
* DCUR_ = 0.001547 0.001433 0.001547 0.001319 *
* DKIS_ = 0.000171 -0.000168 0.000171 -0.000508 *
* DROT_ = 0.000502 0.000289 0.000502 0.000077 *
* *
*******************************************************************************
*******************************************************************************
* *
* CONTRIBUTIONS TO ENERGY IN THE SKYRME FUNCTIONAL *
* *
*******************************************************************************
* *
* TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) *
* -------- ------ ------------ ------------ *
* ERHO_ = -21379.205629 6129.568746 -15343.927485 94.290602 *
* ERHOD = 14327.876593 -3636.933451 10745.907445 -54.964302 *
* ELPR_ = 280.139612 -57.076460 223.957767 -0.894615 *
* ETAU_ = 954.769495 -484.952318 481.724609 -11.907432 *
* ESCU_ = 0.000000 0.000000 0.000000 0.000000 *
* EDIV_ = -54.181233 -28.019215 -81.271850 -0.928598 *
* ============ ============ ============ ============ *
* SUM EVEN: -5870.601162 1922.587303 -3973.609514 25.595655 *
* *
* ESPI_ = -0.201025 1.775002 0.915783 0.658194 *
* ESPID = 0.000000 -1.188852 -0.751593 -0.437259 *
* ELPS_ = 0.000000 -0.300913 -0.180078 -0.120835 *
* ECUR_ = -0.106340 0.097630 -0.053653 0.044944 *
* EKIS_ = 0.000000 0.000000 0.000000 0.000000 *
* EROT_ = -0.032607 -0.018797 -0.048911 -0.002493 *
* ============ ============ ============ ============ *
* SUM ODD: -0.339972 0.364070 -0.118453 0.142551 *
* *
*******************************************************************************
*******************************************************************************
* *
* EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL *
* *
* ALPHA =180.00000 BETA = 0.00000 GAMMA = 0.00000 *
* *
*******************************************************************************
*******************************************************************************
* *
* ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL *
* *
* R_RMS = 4.8657 X_RMS = 2.9229 Y_RMS = 2.3426 Z_RMS = 3.1054 *
* *
* R_GEO = 6.2816 X_GEO = 6.5359 Y_GEO = 5.2382 Z_GEO = 6.9440 *
* *
*******************************************************************************
*******************************************************************************
* *
* MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL *
* *
*******************************************************************************
* *
* Q00 =132.0000 ............. ............. ............. ............. *
* *
* Q10 = ZERO Q11 = ZERO ............. ............. ............. *
* *
* Q20 = 6.9381 Q21 = ZERO Q22 = 6.9864 ............. ............. *
* *
* Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. *
* *
* Q40 = 0.0882 Q41 = ZERO Q42 = 0.2269 Q43 = ZERO Q44 = 0.0986 *
* *
*******************************************************************************
*******************************************************************************
* *
* MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] TOTAL *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* M10 = ZERO M1+1= ZERO ............. ............. ............. *
* M1-1= 2.4513 ............. ............. ............. *
* *
* M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. *
* M2-1= ZERO M2-2= ZERO ............. ............. *
* *
* M30 = ZERO M3+1= ZERO M3+2= ZERO M3+3= ZERO ............. *
* M3-1= 71.7655 M3-2= ZERO M3-3= 56.1197 ............. *
* *
*******************************************************************************
* *
* CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT TOTAL *
* *
* M1_X = 0.0000 M1_Y = 7.0950 M1_Z = 0.0000 M1 = 7.0950 *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) TOTAL *
* *
*******************************************************************************
* *
* B10 = ZERO B11 = ZERO ............. ............. ............. *
* *
* B20 = 0.1760 B21 = ZERO B22 = 0.1772 ............. ............. *
* *
* B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. *
* *
* B40 = 0.0180 B41 = ZERO B42 = 0.0654 B43 = ZERO B44 = 0.0284 *
* *
*******************************************************************************
*******************************************************************************
* *
* 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.9381 Q2+1= ZERO Q2+2= 6.9864 ............. ............. *
* 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.0882 Q4+1= ZERO Q4+2= 0.2269 Q4+3= ZERO Q4+4= 0.0986 *
* Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) [INTRINSIC FRAME] TOTAL *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* B10 = ZERO B1+1= ZERO ............. ............. ............. *
* B1-1= ZERO ............. ............. ............. *
* *
* B20 = 0.1760 B2+1= ZERO B2+2= 0.1772 ............. ............. *
* B2-1= ZERO B2-2= ZERO ............. ............. *
* *
* B30 = ZERO B3+1= ZERO B3+2= ZERO B3+3= ZERO ............. *
* B3-1= ZERO B3-2= ZERO B3-3= ZERO ............. *
* *
* B40 = 0.0180 B4+1= ZERO B4+2= 0.0654 B4+3= ZERO B4+4= 0.0284 *
* B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO *
* *
*******************************************************************************
*******************************************************************************
* *
* MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] [INTRINSIC FRAME] TOTAL *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* M10 = ZERO M1+1= ZERO ............. ............. ............. *
* M1-1= -2.4513 ............. ............. ............. *
* *
* M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. *
* M2-1= ZERO M2-2= ZERO ............. ............. *
* *
* M30 = ZERO M3+1= ZERO M3+2= ZERO M3+3= ZERO ............. *
* M3-1=-71.7655 M3-2= ZERO M3-3=-56.1197 ............. *
* *
*******************************************************************************
* *
* CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT [INTRINSIC FRAME] TOTAL *
* *
* M1_X = 0.0000 M1_Y = -7.0950 M1_Z = 0.0000 M1 = 7.0950 *
* *
*******************************************************************************
*******************************************************************************
* *
* ANGULAR MOMENTA AND THE FIRST MOMENTS OF INERTIA *
* *
*******************************************************************************
* *
* SPINS J(1) *
* --------------------------- --------------------------- *
* ORBITAL INTRINSIC TOTAL ORBITAL INTRINSIC TOTAL *
* *
* NEUTRONS -0.02881 0.06333 0.03451 *
* PROTONS 5.12183 0.39661 5.51844 *
* -------- *
* TOTAL 5.09301 0.45994 5.55295 *
* *
*******************************************************************************
* *
* 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.063 0.029 -0.035 0.000 0.000 0.000 *
* PROTONS 0.000 0.000 0.000 -0.397 -5.122 -5.518 0.000 0.000 0.000 *
* -------- *
* TOTAL 0.000 0.000 0.000 -0.460 -5.093 -5.553 0.000 0.000 0.000 *
* *
*******************************************************************************
*******************************************************************************
* *
* NEUTRON CONFIGURATIONS *
* ======================= *
* P S 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 *
* --- -------------------------------------------------------------- *
* *
* 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 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 *
* *
* 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 0 0 0 0 0 0 0 0 0 0 0 *
* *
* CONF: - - 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 *
* VACC: - - 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 *
* *
*******************************************************************************
* *
* PROTON CONFIGURATIONS *
* ======================= *
* P S 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 *
* --- -------------------------------------------------------------- *
* *
* 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 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 0 0 0 0 0 0 0 0 0 0 0 *
* *
* CONF: - - 1 1 1 1 1 1 1 1 1 1 0 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)= 1459.889679 (PRO)= 957.681284 (TOT)= 2417.570963 *
* SUM EPS: (NEU)= -1763.878912 (PRO)= -1046.934101 (TOT)= -2810.813013 *
* 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)= -8.685433 (PRO)= -5.138303 AND AWAITS YOUR AD *
* *
* COULOMB: (DIR)= 457.896556 (EXC)= -21.973654 (TOT)= 435.922902 *
* *
* CONSTR. (MULT)= 0.002011 SLOPE= 0.075585 CORR.= -0.262536 *
* CONSTR. (SURF)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 *
* CONSTR. (SPIN)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 *
* *
* REARRANGEMENT ENERGY FROM THE SKYRME DENSITY-DEPENDENT TERMS= 890.812858 *
* ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -1094.493885 *
* *
* SPIN-ORB (EVE)= -82.200448 (ODD)= -0.051404 (TOT)= -82.251852 *
* SKYRME: (EVE)= -3948.013860 (ODD)= 0.024098 (TOT)= -3947.989762 *
* *
* TOTAL: (STAB)= -5.14E-07 (SP)= -1094.495897 (FUN)= -1094.495897 *
* *
*******************************************************************************
*******************************************************************************
* *
* CALCULATION OF THE DATA SET NO. 1 HAS BEEN SUCCESSFULLY COMPLETED *
* *
*******************************************************************************
*******************************************************************************
* *
* CLASSICAL NUCLEAR SURFACE DEFINED FOR: N = 75 Z = 57 *
* *
*******************************************************************************
* *
* SPHERICAL SHAPE REQUESTED FOR THE SURFACE *
* *
*******************************************************************************
* *
* HOMEGA= 9.6629 FCHOM0= 1.2000 *
* *
* OSCILLATOR FREQUENCIES: HBAROX= 9.6629 HBAROY= 9.6629 HBAROZ= 9.6629 *
* *
* MOMENTS OF INERTIA: XMOMFC= 49.9200 YMOMFC= 49.9200 ZMOMFC= 49.9200 *
* *
* CENTERS OF MASS: CMSXFC= 0.0000 CMSYFC= 0.0000 CMSZFC= 0.0000 *
* *
*******************************************************************************
*******************************************************************************
* *
* 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) *
* *
*******************************************************************************
*******************************************************************************
* *
* COEFFICIENTS DEFINING THE SKYRME FUNCTIONAL *
* *
*******************************************************************************
* *
* TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) *
* -------- ------ ------------ ------------ *
* CRHO_ = -1382.012500 780.275000 -991.875000 390.137500 *
* CRHOD = 1299.583333 -649.791667 974.687500 -324.895833 *
* CLPR_ = -85.312500 34.218750 -68.203125 17.109375 *
* CTAU_ = 68.750000 -68.125000 34.687500 -34.062500 *
* CSCU_ = 0.000000 0.000000 0.000000 0.000000 *
* CDIV_ = -65.000000 -65.000000 -97.500000 -32.500000 *
* *
* CSPI_ = -59.512500 661.250000 271.112500 330.625000 *
* CSPID = 0.000000 -649.791667 -324.895833 -324.895833 *
* CLPS_ = 0.000000 34.218750 17.109375 17.109375 *
* CCUR_ = -68.750000 68.125000 -34.687500 34.062500 *
* CKIS_ = 0.000000 0.000000 0.000000 0.000000 *
* CROT_ = -65.000000 -65.000000 -97.500000 -32.500000 *
* *
*******************************************************************************
*******************************************************************************
* *
* PHYSICAL CONSTANTS: H_BARC=197.32891000 HBCOE2=137.03602000 *
* *
* XMASSN=938.90590000 XMASSP=938.27231000 *
* *
* HBMASS= 20.73000000 HBMRPA= 20.57295455 *
* *
* ECHAR2= 1.43997841 COULEX= -1.06350868 *
* *
*******************************************************************************
* *
* OSCILLATOR LENGTHS: X= 2.0716964 Y= 2.0716964 Z= 2.0716964 *
* *
* OSCILLATOR CONSTANTS: X= 0.4826962 Y= 0.4826962 Z= 0.4826962 *
* *
* OSCILLATOR FREQUENCIES: X= 9.6628920 Y= 9.6628920 Z= 9.6628920 *
* *
*******************************************************************************
* *
* BASIS CUT-OFF CONTROL PARAMETERS: NXMAXX = 10 NYMAXX = 10 NZMAXX = 10 *
* *
* OPTIMUM NUMBERS OF GAUSS POINTS: NXHERM = 22 NYHERM = 22 NZHERM = 22 *
* *
* NLIMIT= 286 LDBASE= 286 MCOUNT= 1331 *
* *
* ENECUT= 800.0000 ELIMIT= 111.1233 *
* *
*******************************************************************************
*******************************************************************************
* *
* SHAPE OF THE OSCILLATOR-BASIS DIAMOND (THREE CONSEQUTIVE PROJECTIONS) *
* *
*******************************************************************************
* *
* NX ===>>> 0 1 2 3 4 5 6 7 8 9 10 *
* MAX.NY => 10 9 8 7 6 5 4 3 2 1 0 *
* --------------------------------- *
* NY= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 *
* NY= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 *
* NY= 2 ( 8) | 8 7 6 5 4 3 2 1 0 *
* NY= 3 ( 7) | 7 6 5 4 3 2 1 0 *
* NY= 4 ( 6) | 6 5 4 3 2 1 0 *
* NY= 5 ( 5) | 5 4 3 2 1 0 *
* NY= 6 ( 4) | 4 3 2 1 0 *
* NY= 7 ( 3) | 3 2 1 0 *
* NY= 8 ( 2) | 2 1 0 *
* NY= 9 ( 1) | 1 0 *
* NY=10 ( 0) | 0 *
* *
* *
* NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 *
* MAX.NY => 10 9 8 7 6 5 4 3 2 1 0 *
* --------------------------------- *
* NY= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 *
* NY= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 *
* NY= 2 ( 8) | 8 7 6 5 4 3 2 1 0 *
* NY= 3 ( 7) | 7 6 5 4 3 2 1 0 *
* NY= 4 ( 6) | 6 5 4 3 2 1 0 *
* NY= 5 ( 5) | 5 4 3 2 1 0 *
* NY= 6 ( 4) | 4 3 2 1 0 *
* NY= 7 ( 3) | 3 2 1 0 *
* NY= 8 ( 2) | 2 1 0 *
* NY= 9 ( 1) | 1 0 *
* NY=10 ( 0) | 0 *
* *
* *
* NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 *
* MAX.NX => 10 9 8 7 6 5 4 3 2 1 0 *
* --------------------------------- *
* NX= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 *
* NX= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 *
* NX= 2 ( 8) | 8 7 6 5 4 3 2 1 0 *
* NX= 3 ( 7) | 7 6 5 4 3 2 1 0 *
* NX= 4 ( 6) | 6 5 4 3 2 1 0 *
* NX= 5 ( 5) | 5 4 3 2 1 0 *
* NX= 6 ( 4) | 4 3 2 1 0 *
* NX= 7 ( 3) | 3 2 1 0 *
* NX= 8 ( 2) | 2 1 0 *
* NX= 9 ( 1) | 1 0 *
* NX=10 ( 0) | 0 *
* *
*******************************************************************************
*******************************************************************************
* *
* NUCLIDE: N = 75 Z = 57 NUMBER OF ITERATIONS = 5000 CONTINUATION? = 1 *
* *
*******************************************************************************
* *
* ITERATIONS WILL STOP WHEN THE STABILITY BECOMES LOWER THAN EPSITE=0.000001 *
* FOR 5 CONSECUTIVE ITERATIONS *
* *
*******************************************************************************
* *
* SLOW-DOWN PARAMETERS: TIME-EVEN = 0.50 TIME-ODD = 0.50 PAIRING = 0.50 *
* *
*******************************************************************************
* *
* COULOMB PARAMETERS: NUMCOU = 80 NUMETA = 79 BOUCOU = 20.0 FURMAX =0.25 *
* *
*******************************************************************************
* *
* MAXIMUM NUMBERS OF MULTIPOLE MOMENTS USED FOR THE COULOMB FIELD, NMUCOU=4 *
* FOR THE CONSTRAINTS, NMUCON=2 *
* FOR THE OUTPUT INFO, NMUPRI=4 *
* *
*******************************************************************************
* *
* MAXIMUM NUMBER OF MAGNETIC MOMENTS USED: FOR THE CONSTRAINTS, NMACON=0 *
* FOR THE OUTPUT INFO, NMAPRI=3 *
* *
*******************************************************************************
* *
* PRINTING THE RESULTS FOR THE FOLLOWING ITERATIONS: THE FIRST: YES *
* THE MIDDLE: NO *
* AND/OR THE LAST: YES *
* *
*******************************************************************************
* *
* PRINTING THE MOMENTS, RADII, AND DEFORMATIONS FOR: NEUTRONS: YES *
* PROTONS : YES *
* MASS : YES *
* *
*******************************************************************************
* *
* PRINTING THE MOMENTS, MOMENTA AND DEFORMATIONS IN THE INTRINSIC FRAME: YES *
* *
*******************************************************************************
* *
* NUMBERS OF LEVELS IN TIME-REVERSAL/CHARGE BLOCKS: TIME-UP TIME-DOWN *
* NEUTRONS: 75 75 *
* PROTONS: 57 57 *
* *
*******************************************************************************
* *
* CALCULATIONS WITHOUT PAIRING *
* *
*******************************************************************************
* *
* CALCULATIONS WITH PARITY SYMMETRY *
* *
*******************************************************************************
* *
* CALCULATIONS WITH TIME-ODD SYMMETRIES: T*SIMPLEX_X => NON CONSERVED *
* T*SIMPLEX_Y => NON CONSERVED *
* T*SIMPLEX_Z => NON CONSERVED *
* *
*******************************************************************************
* *
* CALCULATIONS WITH BROKEN TIME-REVERSAL SYMMETRY *
* *
*******************************************************************************
* *
* 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 *
* *
*******************************************************************************
*******************************************************************************
* *
* CALLING SUBROUTINE "RECORD" TO READ THE REPLAY FILE OBTAINED PREVIOUSLY *
* la132-a.rec *
* *
*******************************************************************************
*******************************************************************************
* *
* SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS *
* *
*******************************************************************************
* *
* NO) ENERGY (P=+,P=-) | N,nz,lz,OMEG> *
* *
* 64) -11.847 ( 0, 30) | 5, 2, 3, 7/2> -100 1.282 2.885 0.564 *
* 0.140 0.213 0.073 *
* 65) -11.825 ( 35, 0) | 4, 1, 3, 5/2> 100 -0.027 0.232 0.023 *
* -0.064 -0.019 -0.005 *
* 66) -11.569 ( 0, 31) | 5, 3, 2, 5/2> -100 1.813 -3.581 0.620 *
* 0.169 -0.209 0.056 *
* 67) -11.420 ( 36, 0) | 4, 1, 3, 5/2> 100 -0.125 -1.310 -0.184 *
* 0.068 -0.011 0.027 *
* 68) -11.093 ( 0, 32) | 5, 2, 3, 7/2> -100 1.478 0.246 0.580 *
* 0.133 0.032 0.056 *
* 69) -10.794 ( 37, 0) | 4, 0, 4, 7/2> 100 0.009 0.078 0.275 *
* 0.008 -0.113 -0.194 *
* 70) -10.698 ( 0, 33) | 5, 1, 4, 9/2> -100 -0.210 -0.919 0.200 *
* -0.032 -0.069 0.018 *
* 71) -10.488 ( 38, 0) | 4, 0, 4, 7/2> 100 -0.176 -1.054 -0.323 *
* -0.023 0.077 0.213 *
* 72) -10.308 ( 0, 34) | 5, 1, 4, 9/2> -100 -0.864 -0.682 1.423 *
* -0.098 -0.065 0.155 *
* 73) -9.916 ( 0, 35) | 5, 1, 4, 9/2> -100 -2.587 -0.502 -0.562 *
* -0.289 -0.020 -0.071 *
* 74) -9.197 ( 0, 36) | 5, 1, 4, 9/2> -100 -1.179 -0.599 -2.068 *
* -0.121 -0.041 -0.245 *
* 75) -9.131 ( 0, 37) | 5, 0, 5,11/2> -100 -0.670 -0.624 3.449 *
* -0.062 -0.033 0.331 *
* 76) -8.789 ( 39, 0) | 4, 1, 1, 1/2> 100 0.121 0.388 0.277 *
* -0.226 0.016 -0.166 *
* 77) -8.601 ( 40, 0) | 4, 1, 1, 1/2> 100 -0.277 0.250 -0.584 *
* 0.226 0.012 0.107 *
* 78) -8.168 ( 41, 0) | 4, 0, 4, 7/2> 100 -0.100 0.421 0.130 *
* -0.042 0.058 -0.122 *
* 79) -8.055 ( 0, 38) | 5, 2, 1, 1/2> -100 0.197 3.282 0.089 *
* 0.046 0.329 0.033 *
* 80) -7.823 ( 42, 0) | 4, 0, 4, 7/2> 100 -0.262 -1.068 -0.428 *
* 0.044 -0.102 0.184 *
* 81) -7.725 ( 0, 39) | 5, 0, 5,11/2> -100 -0.434 -0.520 -4.845 *
* -0.033 -0.030 -0.494 *
* 82) -7.101 ( 43, 0) | 6, 4, 0, 1/2> 100 0.507 6.406 0.317 *
* 0.053 0.495 0.038 *
* 83) -7.029 ( 0, 40) | 5, 5, 0, 1/2> -100 0.300 -3.059 0.106 *
* 0.072 -0.322 0.049 *
* 84) -6.523 ( 0, 41) | 5, 4, 1, 1/2> -100 0.409 3.130 0.042 *
* 0.017 -0.248 -0.022 *
* 85) -6.344 ( 0, 42) | 5, 5, 0, 1/2> -100 0.541 1.013 0.238 *
* 0.148 0.240 0.077 *
* 86) -6.024 ( 0, 43) | 5, 3, 0, 1/2> -100 -0.207 0.774 0.078 *
* -0.220 -0.242 -0.047 *
* 87) -5.700 ( 44, 0) | 4, 0, 2, 3/2> 100 0.007 0.219 0.278 *
* -0.187 0.130 -0.307 *
* 88) -5.470 ( 0, 44) | 5, 4, 1, 1/2> -100 0.231 -1.299 0.300 *
* -0.041 0.230 0.020 *
* 89) -5.378 ( 45, 0) | 4, 0, 2, 3/2> 100 -0.563 -0.225 -0.927 *
* 0.199 -0.134 0.316 *
* 90) -5.307 ( 0, 45) | 5, 1, 0, 1/2> -100 0.054 0.660 -0.040 *
* -0.002 -0.237 -0.090 *
* 91) -4.832 ( 46, 0) | 6, 3, 3, 7/2> 100 0.912 5.279 0.535 *
* 0.071 0.302 0.044 *
* 92) -4.646 ( 0, 46) | 5, 2, 1, 1/2> -100 0.127 -2.686 0.412 *
* -0.006 0.273 0.049 *
* 93) -4.587 ( 47, 0) | 6, 4, 0, 1/2> 100 0.748 -6.351 0.419 *
* 0.080 -0.489 0.051 *
* 94) -4.118 ( 0, 47) | 5, 0, 1, 3/2> -100 -0.016 0.248 0.737 *
* -0.012 -0.081 0.151 *
* 95) -3.852 ( 0, 48) | 5, 1, 2, 3/2> -100 0.029 0.732 0.408 *
* -0.003 0.165 0.068 *
* 96) -3.671 ( 0, 49) | 5, 0, 1, 3/2> -100 -0.007 -0.573 -0.571 *
* -0.036 0.001 -0.271 *
* 97) -3.546 ( 48, 0) | 6, 3, 3, 7/2> 100 1.294 4.048 0.574 *
* 0.118 0.257 0.061 *
* 98) -3.285 ( 0, 50) | 5, 4, 1, 1/2> -100 0.123 -0.825 -0.045 *
* -0.013 -0.150 -0.030 *
* 99) -3.052 ( 49, 0) | 6, 4, 2, 5/2> 100 1.915 -4.811 0.759 *
* 0.160 -0.279 0.064 *
* 100) -2.987 ( 0, 51) | 5, 2, 3, 5/2> -100 0.081 0.812 0.001 *
* 0.026 -0.149 -0.027 *
* 101) -2.803 ( 0, 52) | 5, 3, 2, 3/2> -100 0.070 -0.154 -0.036 *
* -0.005 0.231 -0.034 *
* 102) -2.625 ( 50, 0) | 6, 2, 4, 9/2> 100 2.038 2.115 0.930 *
* 0.158 0.127 0.075 *
* 103) -2.361 ( 0, 53) | 5, 1, 4, 7/2> -100 0.014 1.500 0.104 *
* -0.015 -0.172 0.016 *
* 104) -2.272 ( 51, 0) | 6, 2, 4, 9/2> 100 1.602 -2.125 0.641 *
* 0.132 -0.126 0.062 *
* 105) -1.891 ( 52, 0) | 6, 1, 5,11/2> 100 0.005 -0.593 0.280 *
* -0.013 -0.052 0.025 *
* 106) -1.859 ( 0, 54) | 5, 2, 3, 5/2> -100 0.047 -1.401 0.077 *
* 0.012 0.200 -0.007 *
* 107) -1.503 ( 0, 55) | 5, 0, 3, 5/2> -100 -0.098 0.338 0.137 *
* -0.026 -0.292 0.000 *
* 108) -1.424 ( 53, 0) | 6, 1, 5,11/2> 100 -0.738 -0.392 0.511 *
* -0.076 -0.028 0.038 *
* 109) -1.069 ( 54, 0) | 6, 1, 5,11/2> 100 -1.989 -0.890 2.151 *
* -0.185 -0.045 0.188 *
* 110) -1.062 ( 0, 56) | 5, 0, 3, 7/2> -100 0.003 -2.426 0.121 *
* -0.032 0.223 0.006 *
* 111) -0.676 ( 55, 0) | 6, 1, 5,11/2> 100 -3.075 -0.699 -2.052 *
* -0.286 -0.044 -0.198 *
* 112) -0.378 ( 0, 57) | 5, 2, 1, 1/2> -100 0.256 1.075 0.155 *
* -0.211 0.073 -0.019 *
* 113) -0.138 ( 0, 58) | 5, 1, 0, 1/2> -100 -0.159 0.033 -0.159 *
* 0.183 -0.057 -0.073 *
* 114) -0.008 ( 56, 0) | 6, 2, 0, 1/2> 100 0.207 4.013 0.392 *
* 0.044 0.362 0.057 *
* *
*******************************************************************************
*******************************************************************************
* *
* MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] NEUTRONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* Q00 = 75.0000 ............. ............. ............. ............. *
* ............. ............. ............. ............. *
* *
* Q10 = ZERO Q1+1= ZERO ............. ............. ............. *
* Q1-1= ZERO ............. ............. ............. *
* *
* Q20 = 3.8302 Q2+1= 0.0167 Q2+2= 4.0495 ............. ............. *
* Q2-1= -0.0099 Q2-2= 0.0105 ............. ............. *
* *
* Q30 = ZERO Q3+1= ZERO Q3+2= ZERO Q3+3= ZERO ............. *
* Q3-1= ZERO Q3-2= ZERO Q3-3= ZERO ............. *
* *
* Q40 = 0.0451 Q4+1= 0.0041 Q4+2= 0.1212 Q4+3= -0.0025 Q4+4= 0.0524 *
* Q4-1= -0.0011 Q4-2=-3.0E-04 Q4-3= 0.0014 Q4-4= 3.2E-04 *
* *
*******************************************************************************
*******************************************************************************
* *
* MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] NEUTRONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* M10 = -0.9321 M1+1= 0.0572 ............. ............. ............. *
* M1-1= -0.1041 ............. ............. ............. *
* *
* M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. *
* M2-1= ZERO M2-2= ZERO ............. ............. *
* *
* M30 = 43.5911 M3+1= -8.3644 M3+2=-15.8995 M3+3= -0.6233 ............. *
* M3-1= 0.3627 M3-2= 1.0374 M3-3= -1.7719 ............. *
* *
*******************************************************************************
* *
* CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT NEUTRONS *
* *
* M1_X = -0.1656 M1_Y = -0.3014 M1_Z = -1.9076 M1 = 1.9384 *
* *
*******************************************************************************
*******************************************************************************
* *
* ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS NEUTRONS *
* *
* R_RMS = 4.9004 X_RMS = 2.9516 Y_RMS = 2.3653 Z_RMS = 3.1156 *
* *
* R_GEO = 6.3264 X_GEO = 6.6001 Y_GEO = 5.2891 Z_GEO = 6.9667 *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) NEUTRONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* B10 = ZERO B1+1= ZERO ............. ............. ............. *
* B1-1= ZERO ............. ............. ............. *
* *
* B20 = 0.1686 B2+1= -0.0025 B2+2= 0.1782 ............. ............. *
* B2-1= 0.0015 B2-2= 4.6E-04 ............. ............. *
* *
* B30 = ZERO B3+1= ZERO B3+2= ZERO B3+3= ZERO ............. *
* B3-1= ZERO B3-2= ZERO B3-3= ZERO ............. *
* *
* B40 = 0.0157 B4+1= 0.0020 B4+2= 0.0598 B4+3= -0.0013 B4+4= 0.0258 *
* B4-1=-5.2E-04 B4-2=-1.5E-04 B4-3= 7.1E-04 B4-4= 1.6E-04 *
* *
*******************************************************************************
*******************************************************************************
* *
* SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS *
* *
*******************************************************************************
* *
* NO) ENERGY (P=+,P=-) | N,nz,lz,OMEG> *
* *
* 38) -11.945 ( 0, 21) | 3, 1, 0, 1/2> -100 -0.089 1.251 -0.094 *
* -0.021 0.342 -0.039 *
* 39) -11.689 ( 0, 22) | 3, 1, 0, 1/2> -100 -0.095 -0.604 -0.095 *
* 0.012 -0.329 0.019 *
* 40) -11.443 ( 18, 0) | 4, 2, 2, 5/2> 100 1.106 -3.034 0.615 *
* 0.110 -0.181 0.062 *
* 41) -11.168 ( 0, 23) | 3, 0, 1, 3/2> -100 -0.007 -0.309 -0.059 *
* -0.206 0.077 -0.269 *
* 42) -11.127 ( 0, 24) | 3, 0, 3, 5/2> -100 -0.164 -0.339 -0.188 *
* 0.203 -0.058 0.281 *
* 43) -11.125 ( 19, 0) | 4, 1, 3, 7/2> 100 1.046 1.812 0.418 *
* 0.141 0.114 0.074 *
* 44) -10.663 ( 20, 0) | 4, 1, 3, 7/2> 100 0.594 -1.243 0.452 *
* 0.074 -0.107 0.070 *
* 45) -10.222 ( 21, 0) | 4, 1, 3, 7/2> 100 -1.013 0.040 -0.110 *
* -0.156 -0.054 -0.041 *
* 46) -9.593 ( 22, 0) | 4, 1, 3, 7/2> 100 -1.088 -0.270 -0.947 *
* -0.146 -0.042 -0.140 *
* 47) -9.433 ( 23, 0) | 4, 0, 4, 9/2> 100 -1.312 -0.362 2.781 *
* -0.144 -0.038 0.340 *
* 48) -8.333 ( 24, 0) | 4, 0, 4, 9/2> 100 -0.478 -0.294 -3.888 *
* -0.047 -0.042 -0.485 *
* 49) -7.491 ( 25, 0) | 4, 2, 2, 3/2> 100 0.199 3.338 0.134 *
* -0.052 -0.118 -0.035 *
* 50) -7.187 ( 0, 25) | 3, 0, 1, 1/2> -100 -0.112 0.458 -0.089 *
* -0.101 0.386 -0.154 *
* 51) -7.042 ( 0, 26) | 3, 0, 1, 1/2> -100 -0.156 -0.458 -0.193 *
* 0.097 -0.383 0.151 *
* 52) -6.768 ( 26, 0) | 4, 1, 1, 1/2> 100 0.098 2.363 0.017 *
* 0.052 0.183 0.034 *
* 53) -6.405 ( 27, 0) | 4, 4, 0, 1/2> 100 0.396 -2.409 0.290 *
* -0.120 -0.093 -0.039 *
* 54) -6.187 ( 28, 0) | 4, 3, 1, 1/2> 100 -0.018 -3.056 -0.155 *
* 0.119 0.106 0.040 *
* 55) -5.888 ( 0, 27) | 5, 3, 0, 1/2> -100 0.516 5.422 0.362 *
* 0.070 0.491 0.055 *
* 56) -5.649 ( 29, 0) | 4, 1, 3, 5/2> 100 0.337 0.804 0.238 *
* -0.041 -0.178 -0.026 *
* 57) -5.175 ( 30, 0) | 4, 4, 0, 1/2> 100 0.210 -0.267 0.120 *
* 0.039 0.212 0.027 *
* 58) -4.816 ( 0, 28) | 5, 3, 0, 1/2> -100 0.556 -5.403 0.346 *
* 0.069 -0.492 0.049 *
* 59) -4.369 ( 31, 0) | 4, 0, 0, 1/2> 100 0.162 -0.004 0.466 *
* 0.056 -0.255 -0.003 *
* 60) -4.146 ( 0, 29) | 5, 3, 2, 5/2> -100 0.700 4.345 0.462 *
* 0.054 0.230 0.035 *
* 61) -3.967 ( 32, 0) | 4, 0, 0, 1/2> 100 -0.061 -0.739 0.502 *
* -0.066 0.277 -0.005 *
* 62) -3.503 ( 33, 0) | 4, 2, 2, 3/2> 100 -0.003 2.091 0.073 *
* -0.123 0.103 -0.112 *
* 63) -3.221 ( 34, 0) | 4, 2, 0, 1/2> 100 -0.249 1.193 -0.295 *
* 0.129 -0.047 0.108 *
* 64) -3.101 ( 0, 30) | 5, 3, 2, 5/2> -100 1.113 -4.124 0.697 *
* 0.100 -0.263 0.066 *
* 65) -3.002 ( 35, 0) | 4, 1, 3, 5/2> 100 0.038 -0.960 0.117 *
* -0.153 0.062 -0.143 *
* 66) -2.922 ( 0, 31) | 5, 2, 3, 7/2> -100 0.876 3.102 0.485 *
* 0.093 0.185 0.062 *
* 67) -2.785 ( 36, 0) | 4, 1, 3, 5/2> 100 -0.273 -1.638 -0.341 *
* 0.139 -0.063 0.141 *
* 68) -2.234 ( 0, 32) | 5, 2, 3, 7/2> -100 1.785 -1.921 0.558 *
* 0.177 -0.145 0.063 *
* 69) -2.094 ( 37, 0) | 4, 0, 4, 7/2> 100 0.029 -0.042 0.169 *
* -0.015 -0.004 -0.259 *
* 70) -1.959 ( 0, 33) | 5, 1, 4, 9/2> -100 0.648 0.526 0.170 *
* 0.050 -0.004 0.005 *
* 71) -1.813 ( 38, 0) | 4, 0, 4, 7/2> 100 -0.214 -0.679 -0.512 *
* 0.017 0.017 0.258 *
* 72) -1.469 ( 0, 34) | 5, 1, 4, 9/2> -100 -0.336 -0.490 0.967 *
* -0.052 -0.065 0.093 *
* 73) -0.974 ( 0, 35) | 5, 1, 4, 9/2> -100 -2.251 -0.300 -0.499 *
* -0.239 -0.032 -0.056 *
* 74) -0.394 ( 0, 36) | 5, 0, 5,11/2> -100 -0.792 -0.397 1.502 *
* -0.075 -0.037 0.136 *
* 75) -0.264 ( 0, 37) | 5, 0, 5,11/2> -100 -1.858 -0.395 0.023 *
* -0.188 -0.049 -0.019 *
* *
*******************************************************************************
*******************************************************************************
* *
* MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] PROTONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* Q00 = 57.0000 ............. ............. ............. ............. *
* ............. ............. ............. ............. *
* *
* Q10 = ZERO Q1+1= ZERO ............. ............. ............. *
* Q1-1= ZERO ............. ............. ............. *
* *
* Q20 = 3.0595 Q2+1= -0.0046 Q2+2= 2.9283 ............. ............. *
* Q2-1= 0.0076 Q2-2= -0.0286 ............. ............. *
* *
* Q30 = ZERO Q3+1= ZERO Q3+2= ZERO Q3+3= ZERO ............. *
* Q3-1= ZERO Q3-2= ZERO Q3-3= ZERO ............. *
* *
* Q40 = 0.0429 Q4+1=-4.7E-04 Q4+2= 0.1043 Q4+3= 9.4E-04 Q4+4= 0.0462 *
* Q4-1= -0.0018 Q4-2=-9.3E-04 Q4-3= -0.0014 Q4-4= -0.0022 *
* *
*******************************************************************************
*******************************************************************************
* *
* MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] PROTONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* M10 = 0.7970 M1+1= -0.8802 ............. ............. ............. *
* M1-1= 2.7620 ............. ............. ............. *
* *
* M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. *
* M2-1= ZERO M2-2= ZERO ............. ............. *
* *
* M30 = 27.1411 M3+1=-10.5805 M3+2= 11.1560 M3+3=-21.5096 ............. *
* M3-1= 68.0386 M3-2= -3.0177 M3-3= 64.6727 ............. *
* *
*******************************************************************************
* *
* CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT PROTONS *
* *
* M1_X = 2.5478 M1_Y = 7.9945 M1_Z = 1.6312 M1 = 8.5477 *
* *
*******************************************************************************
*******************************************************************************
* *
* ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS PROTONS *
* *
* R_RMS = 4.8190 X_RMS = 2.8860 Y_RMS = 2.3159 Z_RMS = 3.0871 *
* *
* R_GEO = 6.2213 X_GEO = 6.4534 Y_GEO = 5.1784 Z_GEO = 6.9029 *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) PROTONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* B10 = ZERO B1+1= ZERO ............. ............. ............. *
* B1-1= ZERO ............. ............. ............. *
* *
* B20 = 0.1832 B2+1= 9.6E-04 B2+2= 0.1754 ............. ............. *
* B2-1= -0.0016 B2-2= -0.0017 ............. ............. *
* *
* B30 = ZERO B3+1= ZERO B3+2= ZERO B3+3= ZERO ............. *
* B3-1= ZERO B3-2= ZERO B3-3= ZERO ............. *
* *
* B40 = 0.0211 B4+1=-3.3E-04 B4+2= 0.0724 B4+3= 6.5E-04 B4+4= 0.0321 *
* B4-1= -0.0012 B4-2=-6.4E-04 B4-3=-9.8E-04 B4-4= -0.0015 *
* *
*******************************************************************************
*******************************************************************************
* *
* DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL *
* *
*******************************************************************************
* *
* TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) *
* -------- ------ ------------ ------------ *
* DRHO_ = 15.466439 7.853944 15.466439 0.241450 *
* DRHOD = 11.022456 5.595730 11.022456 0.169004 *
* DLPR_ = -3.279619 -1.665748 -3.279619 -0.051877 *
* DTAU_ = 13.887252 7.118316 13.887252 0.349380 *
* DSCU_ = 0.141024 0.082131 0.141024 0.023239 *
* DDIV_ = 0.832573 0.430599 0.832573 0.028625 *
* *
* DSPI_ = 0.003215 0.002772 0.003215 0.002330 *
* DSPID = 0.002185 0.001881 0.002185 0.001577 *
* DLPS_ = -0.003116 -0.002846 -0.003116 -0.002577 *
* DCUR_ = 0.005919 0.004640 0.005919 0.003361 *
* DKIS_ = 0.003486 0.002875 0.003486 0.002264 *
* DROT_ = 0.000256 0.000224 0.000256 0.000192 *
* *
*******************************************************************************
*******************************************************************************
* *
* CONTRIBUTIONS TO ENERGY IN THE SKYRME FUNCTIONAL *
* *
*******************************************************************************
* *
* TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) *
* -------- ------ ------------ ------------ *
* ERHO_ = -21374.812452 6128.236527 -15340.774487 94.198562 *
* ERHOD = 14324.600436 -3636.058654 10743.450327 -54.908545 *
* ELPR_ = 279.792457 -56.999806 223.680233 -0.887582 *
* ETAU_ = 954.748563 -484.935260 481.714048 -11.900745 *
* ESCU_ = 0.000000 0.000000 0.000000 0.000000 *
* EDIV_ = -54.117256 -27.988938 -81.175885 -0.930310 *
* ============ ============ ============ ============ *
* SUM EVEN: -5869.788253 1922.253869 -3973.105764 25.571380 *
* *
* ESPI_ = -0.191308 1.833218 0.871515 0.770396 *
* ESPID = 0.000000 -1.222383 -0.709934 -0.512449 *
* ELPS_ = 0.000000 -0.097401 -0.053306 -0.044096 *
* ECUR_ = -0.406959 0.316102 -0.205330 0.114472 *
* EKIS_ = 0.000000 0.000000 0.000000 0.000000 *
* EROT_ = -0.016620 -0.014563 -0.024930 -0.006253 *
* ============ ============ ============ ============ *
* SUM ODD: -0.614887 0.814973 -0.121984 0.322070 *
* *
*******************************************************************************
*******************************************************************************
* *
* EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL *
* *
* ALPHA = -2.73332 BETA = 0.48182 GAMMA = 2.80793 *
* *
*******************************************************************************
*******************************************************************************
* *
* ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL *
* *
* R_RMS = 4.8654 X_RMS = 2.9235 Y_RMS = 2.3441 Z_RMS = 3.1033 *
* *
* R_GEO = 6.2812 X_GEO = 6.5371 Y_GEO = 5.2416 Z_GEO = 6.9392 *
* *
*******************************************************************************
*******************************************************************************
* *
* 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.8898 Q2+1= 0.0120 Q2+2= 6.9777 ............. ............. *
* Q2-1= -0.0022 Q2-2= -0.0181 ............. ............. *
* *
* Q30 = ZERO Q3+1= ZERO Q3+2= ZERO Q3+3= ZERO ............. *
* Q3-1= ZERO Q3-2= ZERO Q3-3= ZERO ............. *
* *
* Q40 = 0.0880 Q4+1= 0.0036 Q4+2= 0.2255 Q4+3= -0.0016 Q4+4= 0.0986 *
* Q4-1= -0.0028 Q4-2= -0.0012 Q4-3= 2.7E-05 Q4-4= -0.0019 *
* *
*******************************************************************************
*******************************************************************************
* *
* MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] TOTAL *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* M10 = -0.1351 M1+1= -0.8230 ............. ............. ............. *
* M1-1= 2.6579 ............. ............. ............. *
* *
* M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. *
* M2-1= ZERO M2-2= ZERO ............. ............. *
* *
* M30 = 70.7322 M3+1=-18.9449 M3+2= -4.7436 M3+3=-22.1329 ............. *
* M3-1= 68.4013 M3-2= -1.9803 M3-3= 62.9008 ............. *
* *
*******************************************************************************
* *
* CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT TOTAL *
* *
* M1_X = 2.3821 M1_Y = 7.6931 M1_Z = -0.2765 M1 = 8.0582 *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) TOTAL *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* B10 = ZERO B1+1= ZERO ............. ............. ............. *
* B1-1= ZERO ............. ............. ............. *
* *
* B20 = 0.1748 B2+1= -0.0011 B2+2= 0.1770 ............. ............. *
* B2-1= 2.0E-04 B2-2=-4.6E-04 ............. ............. *
* *
* B30 = ZERO B3+1= ZERO B3+2= ZERO B3+3= ZERO ............. *
* B3-1= ZERO B3-2= ZERO B3-3= ZERO ............. *
* *
* B40 = 0.0179 B4+1= 0.0011 B4+2= 0.0650 B4+3=-4.6E-04 B4+4= 0.0284 *
* B4-1=-8.2E-04 B4-2=-3.5E-04 B4-3= 7.7E-06 B4-4=-5.5E-04 *
* *
*******************************************************************************
*******************************************************************************
* *
* 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.8901 Q2+1= ZERO Q2+2= 6.9776 ............. ............. *
* 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.0879 Q4+1= 0.0013 Q4+2= 0.2256 Q4+3= 7.6E-04 Q4+4= 0.0986 *
* Q4-1= -0.0031 Q4-2=-5.9E-04 Q4-3=-2.1E-04 Q4-4= -0.0014 *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) [INTRINSIC FRAME] TOTAL *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* B10 = ZERO B1+1= ZERO ............. ............. ............. *
* B1-1= ZERO ............. ............. ............. *
* *
* B20 = 0.1748 B2+1= ZERO B2+2= 0.1770 ............. ............. *
* B2-1= ZERO B2-2= ZERO ............. ............. *
* *
* B30 = ZERO B3+1= ZERO B3+2= ZERO B3+3= ZERO ............. *
* B3-1= ZERO B3-2= ZERO B3-3= ZERO ............. *
* *
* B40 = 0.0179 B4+1= 3.7E-04 B4+2= 0.0650 B4+3= 2.2E-04 B4+4= 0.0284 *
* B4-1=-9.0E-04 B4-2=-1.7E-04 B4-3=-5.9E-05 B4-4=-4.0E-04 *
* *
*******************************************************************************
*******************************************************************************
* *
* MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] [INTRINSIC FRAME] TOTAL *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* M10 = -0.1238 M1+1= -0.8203 ............. ............. ............. *
* M1-1= 2.6590 ............. ............. ............. *
* *
* M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. *
* M2-1= ZERO M2-2= ZERO ............. ............. *
* *
* M30 = 71.3625 M3+1=-17.7582 M3+2= -4.6798 M3+3=-21.9344 ............. *
* M3-1= 68.3499 M3-2= -2.2758 M3-3= 63.0111 ............. *
* *
*******************************************************************************
* *
* CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT [INTRINSIC FRAME] TOTAL *
* *
* M1_X = 2.3743 M1_Y = 7.6963 M1_Z = -0.2534 M1 = 8.0582 *
* *
*******************************************************************************
*******************************************************************************
* *
* 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.043 2.145 2.188 0.079 0.485 0.564 0.499 6.111 6.609 *
* PROTONS 0.090 2.045 2.135 0.389 5.823 6.212 0.071 1.237 1.307 *
* -------- *
* TOTAL 0.133 4.190 4.323 0.468 6.308 6.776 0.569 7.347 7.916 *
* *
*******************************************************************************
* *
* ROTATIONAL FREQUENCIES IN THE THREE CARTESIAN DIRECTIONS: *
* *
* OMEGAX = 0.144338 OMEGAY = 0.144338 OMEGAZ = 0.144338 OMEGA = 0.250000 *
* *
*******************************************************************************
* *
* ANGLE BETWEEN ANGULAR MOMENTUM AND ANGULAR FREQUENCY VECTORS= 13.308 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.039 2.093 2.132 0.079 0.486 0.564 0.499 6.129 6.628 *
* PROTONS 0.089 2.027 2.116 0.389 5.825 6.214 0.072 1.256 1.328 *
* -------- *
* TOTAL 0.128 4.120 4.247 0.467 6.311 6.778 0.571 7.385 7.955 *
* *
*******************************************************************************
*******************************************************************************
* *
* 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)= 1459.861035 (PRO)= 957.859660 (TOT)= 2417.720696 *
* SUM EPS: (NEU)= -1764.608782 (PRO)= -1047.801431 (TOT)= -2812.410213 *
* 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)= -8.959750 (PRO)= -4.995772 AND AWAITS YOUR AD *
* *
* COULOMB: (DIR)= 457.912321 (EXC)= -21.972437 (TOT)= 435.939885 *
* *
* CONSTR. (MULT)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 *
* CONSTR. (SURF)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 *
* CONSTR. (SPIN)= -2.744611 SLOPE= 0.433013 CORR.= -1.372305 *
* *
* REARRANGEMENT ENERGY FROM THE SKYRME DENSITY-DEPENDENT TERMS= 890.609950 *
* ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -1096.418329 *
* *
* SPIN-ORB (EVE)= -82.106194 (ODD)= -0.031183 (TOT)= -82.137378 *
* SKYRME: (EVE)= -3947.534384 (ODD)= 0.200086 (TOT)= -3947.334298 *
* *
* TOTAL: (STAB)= -0.232831 (SP)= -1093.906549 (FUN)= -1093.673718 *
* *
*******************************************************************************
*******************************************************************************
* *
* ITERATION RESTARTED FROM THE DISC, FILE NAME (FILREP): *
* la132-a.rec *
* *
*******************************************************************************
* *
* 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 *
* 104 -1093.805889 0.019367 6.816 45.55 11.94 0.250 0.483 0.00 *
* 105 -1093.802823 0.016769 6.811 45.56 11.96 0.250 0.373 0.00 *
* 106 -1093.800539 0.014573 6.807 45.57 11.97 0.250 0.290 0.00 *
* 107 -1093.798822 0.012700 6.804 45.58 11.97 0.250 0.227 0.00 *
* 108 -1093.797520 0.011090 6.801 45.58 11.98 0.250 0.179 0.00 *
* 109 -1093.796523 0.009698 6.798 45.59 11.98 0.250 0.141 0.00 *
* 110 -1093.795751 0.008491 6.796 45.59 11.99 0.250 0.112 0.00 *
* 111 -1093.795147 0.007440 6.794 45.60 11.99 0.250 0.088 0.00 *
* 112 -1093.794670 0.006524 6.792 45.60 11.99 0.250 0.070 0.00 *
* 113 -1093.794290 0.005724 6.790 45.60 11.99 0.250 0.056 0.00 *
* 114 -1093.793984 0.005025 6.789 45.61 12.00 0.250 0.045 0.00 *
* 115 -1093.793736 0.004412 6.788 45.61 12.00 0.250 0.036 0.00 *
* 116 -1093.793534 0.003876 6.787 45.61 12.00 0.250 0.028 0.00 *
* 117 -1093.793366 0.003405 6.786 45.61 12.00 0.250 0.023 0.00 *
* 118 -1093.793228 0.002993 6.785 45.61 12.00 0.250 0.018 0.00 *
* 119 -1093.793112 0.002631 6.785 45.61 12.00 0.250 0.015 0.00 *
* 120 -1093.793015 0.002313 6.784 45.62 12.00 0.250 0.012 0.00 *
* 121 -1093.792933 0.002034 6.784 45.62 12.00 0.250 9.5E-03 0.00 *
* 122 -1093.792863 0.001789 6.783 45.62 12.00 0.250 7.7E-03 0.00 *
* 123 -1093.792804 0.001574 6.783 45.62 12.00 0.250 6.2E-03 0.00 *
* 124 -1093.792754 0.001385 6.782 45.62 12.00 0.250 5.0E-03 0.00 *
* 125 -1093.792711 0.001218 6.782 45.62 12.00 0.250 4.0E-03 0.00 *
* 126 -1093.792673 0.001072 6.782 45.62 12.00 0.250 3.3E-03 0.00 *
* 127 -1093.792641 0.000943 6.782 45.62 12.00 0.250 2.7E-03 0.00 *
* 128 -1093.792614 0.000830 6.781 45.62 12.00 0.250 2.2E-03 0.00 *
* 129 -1093.792590 0.000731 6.781 45.62 12.00 0.250 1.8E-03 0.00 *
* 130 -1093.792569 0.000643 6.781 45.62 12.00 0.250 1.4E-03 0.00 *
* 131 -1093.792551 0.000566 6.781 45.62 12.00 0.250 1.2E-03 0.00 *
* 132 -1093.792536 0.000498 6.781 45.62 12.00 0.250 9.7E-04 0.00 *
* 133 -1093.792522 0.000439 6.781 45.62 12.00 0.250 8.0E-04 0.00 *
* 134 -1093.792511 0.000386 6.781 45.62 12.00 0.250 6.5E-04 0.00 *
* 135 -1093.792501 0.000340 6.781 45.62 12.00 0.250 5.4E-04 0.00 *
* 136 -1093.792492 0.000300 6.780 45.62 12.00 0.250 4.5E-04 0.00 *
* 137 -1093.792484 0.000264 6.780 45.62 12.00 0.250 3.7E-04 0.00 *
* 138 -1093.792477 0.000232 6.780 45.62 12.00 0.250 3.1E-04 0.00 *
* 139 -1093.792471 0.000205 6.780 45.62 12.00 0.250 2.6E-04 0.00 *
* 140 -1093.792466 0.000180 6.780 45.62 12.00 0.250 2.2E-04 0.00 *
* 141 -1093.792462 0.000159 6.780 45.62 12.00 0.250 1.8E-04 0.00 *
* 142 -1093.792458 0.000140 6.780 45.62 12.00 0.250 1.5E-04 0.00 *
* 143 -1093.792454 0.000123 6.780 45.62 12.00 0.250 1.3E-04 0.00 *
* 144 -1093.792451 0.000109 6.780 45.62 12.00 0.250 1.1E-04 0.00 *
* 145 -1093.792449 0.000096 6.780 45.62 12.00 0.250 9.4E-05 0.00 *
* 146 -1093.792446 0.000084 6.780 45.62 12.00 0.250 8.1E-05 0.00 *
* 147 -1093.792444 0.000074 6.780 45.62 12.00 0.250 6.9E-05 0.00 *
* 148 -1093.792443 0.000065 6.780 45.62 12.00 0.250 6.0E-05 0.00 *
* 149 -1093.792441 0.000058 6.780 45.62 12.00 0.250 5.2E-05 0.00 *
* 150 -1093.792440 0.000051 6.780 45.62 12.00 0.250 4.5E-05 0.00 *
* 151 -1093.792439 0.000045 6.780 45.62 12.00 0.250 3.9E-05 0.00 *
* 152 -1093.792438 0.000039 6.780 45.62 12.00 0.250 3.4E-05 0.00 *
* 153 -1093.792437 0.000035 6.780 45.62 12.00 0.250 3.0E-05 0.00 *
* 154 -1093.792436 0.000031 6.780 45.62 12.00 0.250 2.7E-05 0.00 *
* 155 -1093.792435 0.000027 6.780 45.62 12.00 0.250 2.4E-05 0.00 *
* 156 -1093.792435 0.000024 6.780 45.62 12.00 0.250 2.1E-05 0.00 *
* 157 -1093.792434 0.000021 6.780 45.62 12.00 0.250 1.9E-05 0.00 *
* 158 -1093.792434 0.000018 6.780 45.62 12.00 0.250 1.7E-05 0.00 *
* 159 -1093.792433 0.000016 6.780 45.62 12.00 0.250 1.5E-05 0.00 *
* 160 -1093.792433 0.000014 6.780 45.62 12.00 0.250 1.3E-05 0.00 *
* 161 -1093.792433 0.000013 6.780 45.62 12.00 0.250 1.2E-05 0.00 *
* 162 -1093.792433 0.000011 6.780 45.62 12.00 0.250 1.1E-05 0.00 *
* 163 -1093.792432 9.82E-06 6.780 45.62 12.00 0.250 9.5E-06 0.00 *
* 164 -1093.792432 8.66E-06 6.780 45.62 12.00 0.250 8.5E-06 0.00 *
* 165 -1093.792432 7.63E-06 6.780 45.62 12.00 0.250 7.7E-06 0.00 *
* 166 -1093.792432 6.73E-06 6.780 45.62 12.00 0.250 6.9E-06 0.00 *
* 167 -1093.792432 5.93E-06 6.780 45.62 12.00 0.250 6.2E-06 0.00 *
* 168 -1093.792432 5.23E-06 6.780 45.62 12.00 0.250 5.6E-06 0.00 *
* 169 -1093.792431 4.61E-06 6.780 45.62 12.00 0.250 5.1E-06 0.00 *
* 170 -1093.792431 4.06E-06 6.780 45.62 12.00 0.250 4.4E-06 0.00 *
* 171 -1093.792431 3.58E-06 6.780 45.62 12.00 0.250 4.1E-06 0.00 *
* 172 -1093.792431 3.16E-06 6.780 45.62 12.00 0.250 3.7E-06 0.00 *
* 173 -1093.792431 2.78E-06 6.780 45.62 12.00 0.250 3.4E-06 0.00 *
* 174 -1093.792431 2.45E-06 6.780 45.62 12.00 0.250 3.0E-06 0.00 *
* 175 -1093.792431 2.16E-06 6.780 45.62 12.00 0.250 2.6E-06 0.00 *
* 176 -1093.792431 1.91E-06 6.780 45.62 12.00 0.250 2.4E-06 0.00 *
* 177 -1093.792431 1.68E-06 6.780 45.62 12.00 0.250 1.9E-06 0.00 *
* 178 -1093.792431 1.48E-06 6.780 45.62 12.00 0.250 1.9E-06 0.00 *
* 179 -1093.792431 1.31E-06 6.780 45.62 12.00 0.250 1.9E-06 0.00 *
* 180 -1093.792431 1.15E-06 6.780 45.62 12.00 0.250 1.7E-06 0.00 *
* 181 -1093.792431 1.01E-06 6.780 45.62 12.00 0.250 1.5E-06 0.00 *
* 182 -1093.792431 8.94E-07 6.780 45.62 12.00 0.250 1.7E-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 0.000 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> *
* *
* 64) -11.857 ( 35, 0) | 4, 1, 3, 5/2> 100 -0.001 0.762 -0.130 *
* 0.001 0.005 0.097 *
* 65) -11.762 ( 0, 30) | 5, 2, 3, 7/2> -100 0.013 3.191 0.663 *
* 0.001 0.246 0.086 *
* 66) -11.426 ( 0, 31) | 5, 3, 2, 5/2> -100 0.007 -4.118 0.835 *
* 0.001 -0.237 0.066 *
* 67) -11.395 ( 36, 0) | 4, 1, 3, 5/2> 100 -0.005 -1.072 -0.264 *
* -0.001 -0.012 -0.040 *
* 68) -11.047 ( 37, 0) | 4, 0, 4, 7/2> 100 0.005 -0.379 0.373 *
* -0.003 -0.076 -0.192 *
* 69) -10.940 ( 0, 32) | 5, 4, 1, 3/2> -100 0.015 0.980 0.977 *
* 0.001 0.064 0.083 *
* 70) -10.661 ( 0, 33) | 5, 1, 4, 9/2> -100 0.006 -1.888 0.596 *
* 0.001 -0.143 0.076 *
* 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 *
* 81) -7.669 ( 42, 0) | 4, 0, 4, 7/2> 100 -0.012 -0.743 -0.781 *
* 0.005 -0.063 0.357 *
* 82) -7.050 ( 43, 0) | 6, 4, 0, 1/2> 100 0.014 6.424 0.479 *
* 0.001 0.497 0.046 *
* 83) -6.975 ( 0, 40) | 5, 5, 0, 1/2> -100 -0.001 -3.116 0.163 *
* 0.001 -0.327 0.063 *
* 84) -6.439 ( 0, 41) | 5, 4, 1, 1/2> -100 0.005 3.522 0.089 *
* -0.001 -0.259 -0.052 *
* 85) -6.302 ( 0, 42) | 5, 5, 0, 1/2> -100 0.006 0.808 0.380 *
* 0.002 0.290 0.152 *
* 86) -6.028 ( 44, 0) | 4, 0, 0, 1/2> 100 0.004 0.118 0.269 *
* -0.005 0.064 -0.389 *
* 87) -5.997 ( 0, 43) | 5, 3, 0, 1/2> -100 0.002 0.853 0.051 *
* -0.002 -0.272 -0.130 *
* 88) -5.417 ( 0, 44) | 5, 4, 1, 1/2> -100 0.004 -0.991 0.380 *
* 0.000 0.167 0.002 *
* 89) -5.333 ( 0, 45) | 5, 1, 0, 1/2> -100 0.001 0.435 0.028 *
* -0.002 -0.178 -0.102 *
* 90) -5.195 ( 45, 0) | 4, 0, 2, 3/2> 100 -0.019 -0.126 -1.319 *
* 0.006 -0.064 0.422 *
* 91) -4.735 ( 46, 0) | 6, 3, 3, 7/2> 100 0.015 5.347 0.675 *
* 0.001 0.306 0.048 *
* 92) -4.609 ( 0, 46) | 5, 2, 1, 1/2> -100 0.005 -2.778 0.535 *
* 0.001 0.278 0.039 *
* 93) -4.485 ( 47, 0) | 6, 4, 0, 1/2> 100 0.001 -6.394 0.547 *
* 0.000 -0.494 0.059 *
* 94) -4.136 ( 0, 47) | 5, 0, 1, 3/2> -100 0.008 0.405 0.557 *
* 0.000 -0.182 0.019 *
* 95) -3.857 ( 0, 48) | 5, 1, 2, 3/2> -100 0.007 1.715 0.383 *
* 0.000 0.024 0.014 *
* 96) -3.811 ( 0, 49) | 5, 0, 1, 3/2> -100 -0.001 -1.854 0.085 *
* -0.001 0.231 -0.118 *
* 97) -3.460 ( 48, 0) | 6, 2, 4, 9/2> 100 0.015 4.237 0.716 *
* 0.001 0.275 0.072 *
* 98) -3.273 ( 0, 50) | 5, 4, 1, 1/2> -100 0.000 -0.977 0.072 *
* 0.000 -0.152 0.016 *
* 99) -2.980 ( 0, 51) | 5, 3, 2, 3/2> -100 0.002 0.736 0.053 *
* -0.001 -0.151 -0.094 *
* 100) -2.904 ( 49, 0) | 6, 3, 3, 7/2> 100 0.007 -5.191 0.885 *
* 0.001 -0.300 0.066 *
* 101) -2.821 ( 0, 52) | 5, 3, 2, 3/2> -100 0.000 0.137 -0.012 *
* -0.001 0.214 -0.062 *
* 102) -2.469 ( 50, 0) | 6, 2, 4, 9/2> 100 0.021 2.905 1.306 *
* 0.002 0.162 0.097 *
* 103) -2.377 ( 0, 53) | 5, 1, 4, 7/2> -100 0.004 1.555 0.134 *
* 0.001 -0.147 0.058 *
* 104) -2.117 ( 51, 0) | 6, 3, 3, 7/2> 100 0.015 -3.421 1.304 *
* 0.002 -0.212 0.128 *
* 105) -1.942 ( 0, 54) | 5, 0, 3, 5/2> -100 0.002 -1.178 0.204 *
* -0.003 0.159 -0.232 *
* 106) -1.812 ( 52, 0) | 6, 2, 4, 9/2> 100 0.002 0.628 0.127 *
* 0.000 0.025 0.008 *
* 107) -1.595 ( 0, 55) | 5, 0, 3, 7/2> -100 0.004 -0.128 0.318 *
* 0.001 -0.210 0.114 *
* 108) -1.468 ( 53, 0) | 6, 1, 5,11/2> 100 0.050 -1.468 3.700 *
* 0.004 -0.087 0.327 *
* 109) -1.353 ( 54, 0) | 6, 3, 3, 7/2> 100 -0.013 -0.387 -0.934 *
* -0.001 -0.013 -0.105 *
* 110) -1.033 ( 0, 56) | 5, 0, 3, 7/2> -100 0.000 -2.379 0.200 *
* 0.003 0.180 0.171 *
* 111) -0.750 ( 55, 0) | 6, 2, 4, 9/2> 100 -0.049 -0.667 -3.459 *
* -0.005 -0.042 -0.337 *
* 112) -0.470 ( 0, 57) | 5, 1, 4, 7/2> -100 0.006 1.340 0.300 *
* -0.001 0.079 -0.086 *
* 113) -0.195 ( 0, 58) | 5, 1, 0, 1/2> -100 -0.001 -0.254 -0.037 *
* -0.002 -0.085 -0.121 *
* 114) -0.193 ( 56, 0) | 6, 0, 6,13/2> 100 0.083 -0.707 5.970 *
* 0.007 -0.035 0.497 *
* *
*******************************************************************************
*******************************************************************************
* *
* MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] NEUTRONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* Q00 = 75.0000 ............. ............. ............. ............. *
* ............. ............. ............. ............. *
* *
* Q10 = ZERO Q1+1= ZERO ............. ............. ............. *
* Q1-1= ZERO ............. ............. ............. *
* *
* Q20 = 3.7668 Q2+1= 0.0101 Q2+2= 4.0154 ............. ............. *
* Q2-1= -0.0069 Q2-2= -0.0083 ............. ............. *
* *
* Q30 = ZERO Q3+1= ZERO Q3+2= ZERO Q3+3= ZERO ............. *
* Q3-1= ZERO Q3-2= ZERO Q3-3= ZERO ............. *
* *
* Q40 = 0.0424 Q4+1= 0.0022 Q4+2= 0.1177 Q4+3= -0.0020 Q4+4= 0.0525 *
* Q4-1= -0.0054 Q4-2=-4.5E-04 Q4-3= 0.0013 Q4-4=-2.0E-04 *
* *
*******************************************************************************
*******************************************************************************
* *
* MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] NEUTRONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* M10 = -0.7998 M1+1= 0.0080 ............. ............. ............. *
* M1-1= -0.0970 ............. ............. ............. *
* *
* M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. *
* M2-1= ZERO M2-2= ZERO ............. ............. *
* *
* M30 = 29.7412 M3+1= -1.0052 M3+2=-12.8551 M3+3= 0.2307 ............. *
* M3-1= -0.2873 M3-2= 0.0761 M3-3= -3.2046 ............. *
* *
*******************************************************************************
* *
* CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT NEUTRONS *
* *
* M1_X = -0.0232 M1_Y = -0.2809 M1_Z = -1.6369 M1 = 1.6610 *
* *
*******************************************************************************
*******************************************************************************
* *
* ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS NEUTRONS *
* *
* R_RMS = 4.8985 X_RMS = 2.9507 Y_RMS = 2.3697 Z_RMS = 3.1101 *
* *
* R_GEO = 6.3239 X_GEO = 6.5980 Y_GEO = 5.2989 Z_GEO = 6.9543 *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) NEUTRONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* B10 = ZERO B1+1= ZERO ............. ............. ............. *
* B1-1= ZERO ............. ............. ............. *
* *
* B20 = 0.1659 B2+1= -0.0015 B2+2= 0.1769 ............. ............. *
* B2-1= 0.0010 B2-2=-3.6E-04 ............. ............. *
* *
* B30 = ZERO B3+1= ZERO B3+2= ZERO B3+3= ZERO ............. *
* B3-1= ZERO B3-2= ZERO B3-3= ZERO ............. *
* *
* B40 = 0.0148 B4+1= 0.0011 B4+2= 0.0581 B4+3= -0.0010 B4+4= 0.0259 *
* B4-1= -0.0027 B4-2=-2.2E-04 B4-3= 6.3E-04 B4-4=-9.8E-05 *
* *
*******************************************************************************
*******************************************************************************
* *
* SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS *
* *
*******************************************************************************
* *
* NO) ENERGY (P=+,P=-) | N,nz,lz,OMEG> *
* *
* 38) -11.993 ( 0, 21) | 3, 1, 0, 1/2> -100 -0.001 1.215 -0.147 *
* -0.002 0.309 -0.176 *
* 39) -11.721 ( 0, 22) | 3, 1, 0, 1/2> -100 -0.003 -0.507 -0.185 *
* 0.001 -0.299 0.113 *
* 40) -11.416 ( 18, 0) | 4, 2, 2, 5/2> 100 0.011 -3.212 1.008 *
* 0.001 -0.195 0.102 *
* 41) -11.243 ( 0, 23) | 3, 0, 1, 3/2> -100 -0.002 -0.327 -0.145 *
* -0.006 0.035 -0.404 *
* 42) -11.108 ( 19, 0) | 4, 1, 3, 7/2> 100 0.017 2.128 1.039 *
* 0.002 0.126 0.168 *
* 43) -11.102 ( 0, 24) | 3, 0, 3, 5/2> -100 -0.005 -0.386 -0.301 *
* 0.006 -0.016 0.445 *
* 44) -10.629 ( 20, 0) | 4, 1, 3, 7/2> 100 0.007 -1.391 0.615 *
* 0.001 -0.105 0.100 *
* 45) -10.265 ( 21, 0) | 4, 2, 2, 5/2> 100 -0.008 0.218 -0.557 *
* -0.002 -0.029 -0.110 *
* 46) -10.072 ( 22, 0) | 4, 0, 4, 9/2> 100 0.046 -0.524 3.331 *
* 0.005 -0.060 0.395 *
* 47) -9.414 ( 23, 0) | 4, 1, 3, 7/2> 100 -0.038 -0.244 -2.706 *
* -0.005 -0.035 -0.363 *
* 48) -7.823 ( 24, 0) | 4, 0, 4, 9/2> 100 -0.055 -0.331 -3.892 *
* -0.007 -0.049 -0.488 *
* 49) -7.452 ( 25, 0) | 4, 2, 2, 3/2> 100 0.007 3.327 0.226 *
* -0.001 -0.120 -0.062 *
* 50) -7.284 ( 0, 25) | 3, 0, 1, 1/2> -100 -0.003 0.492 -0.258 *
* 0.000 0.416 -0.055 *
* 51) -7.121 ( 0, 26) | 3, 0, 1, 1/2> -100 -0.004 -0.496 -0.281 *
* 0.000 -0.410 0.047 *
* 52) -6.728 ( 26, 0) | 4, 1, 1, 1/2> 100 0.003 2.388 0.002 *
* 0.001 0.183 0.072 *
* 53) -6.391 ( 27, 0) | 4, 4, 0, 1/2> 100 0.004 -2.487 0.464 *
* -0.001 -0.074 -0.054 *
* 54) -6.123 ( 28, 0) | 4, 3, 1, 1/2> 100 -0.006 -3.036 -0.214 *
* 0.001 0.082 0.064 *
* 55) -5.864 ( 0, 27) | 5, 3, 0, 1/2> -100 0.014 5.436 0.562 *
* 0.002 0.492 0.083 *
* 56) -5.630 ( 29, 0) | 4, 1, 3, 5/2> 100 0.007 0.859 0.408 *
* -0.001 -0.175 -0.054 *
* 57) -5.150 ( 30, 0) | 4, 4, 0, 1/2> 100 0.002 -0.175 0.170 *
* 0.001 0.207 0.042 *
* 58) -4.745 ( 0, 28) | 5, 3, 0, 1/2> -100 0.002 -5.421 0.578 *
* 0.001 -0.492 0.081 *
* 59) -4.419 ( 31, 0) | 4, 0, 2, 3/2> 100 0.011 -0.043 0.793 *
* 0.000 -0.255 -0.005 *
* 60) -4.136 ( 0, 29) | 5, 2, 3, 7/2> -100 0.016 4.368 0.801 *
* 0.001 0.234 0.066 *
* 61) -4.035 ( 32, 0) | 4, 0, 2, 3/2> 100 0.012 -0.811 0.898 *
* 0.000 0.276 -0.025 *
* 62) -3.565 ( 33, 0) | 4, 2, 0, 1/2> 100 0.004 1.995 0.143 *
* -0.003 0.070 -0.201 *
* 63) -3.178 ( 34, 0) | 4, 2, 0, 1/2> 100 -0.006 1.326 -0.517 *
* 0.003 -0.016 0.217 *
* 64) -3.078 ( 35, 0) | 4, 1, 3, 5/2> 100 0.001 -1.020 0.171 *
* -0.003 0.045 -0.239 *
* 65) -3.066 ( 0, 30) | 5, 2, 3, 7/2> -100 0.011 -4.226 1.143 *
* 0.001 -0.271 0.109 *
* 66) -2.884 ( 0, 31) | 5, 1, 4, 9/2> -100 0.014 3.226 0.781 *
* 0.002 0.190 0.099 *
* 67) -2.733 ( 36, 0) | 4, 1, 3, 5/2> 100 -0.009 -1.528 -0.524 *
* 0.003 -0.049 0.222 *
* 68) -2.206 ( 37, 0) | 4, 0, 4, 7/2> 100 0.003 -0.146 0.204 *
* -0.005 -0.004 -0.360 *
* 69) -2.118 ( 0, 32) | 5, 1, 4, 9/2> -100 0.009 -3.045 0.850 *
* 0.001 -0.212 0.102 *
* 70) -1.982 ( 0, 33) | 5, 1, 4, 9/2> -100 0.024 1.813 1.563 *
* 0.002 0.089 0.143 *
* 71) -1.784 ( 38, 0) | 4, 0, 4, 7/2> 100 -0.009 -0.653 -0.613 *
* 0.005 0.021 0.349 *
* 72) -1.547 ( 0, 34) | 5, 3, 2, 5/2> -100 0.001 -0.587 0.109 *
* 0.000 -0.080 -0.001 *
* 73) -1.073 ( 0, 35) | 5, 0, 5,11/2> -100 0.046 -0.414 3.313 *
* 0.004 -0.040 0.317 *
* 74) -0.988 ( 0, 36) | 5, 2, 3, 7/2> -100 -0.008 -0.421 -0.526 *
* -0.001 -0.032 -0.067 *
* 75) -0.079 ( 0, 37) | 5, 1, 4, 9/2> -100 -0.053 -0.306 -3.748 *
* -0.006 -0.043 -0.414 *
* *
*******************************************************************************
*******************************************************************************
* *
* MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] PROTONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* Q00 = 57.0000 ............. ............. ............. ............. *
* ............. ............. ............. ............. *
* *
* Q10 = ZERO Q1+1= ZERO ............. ............. ............. *
* Q1-1= ZERO ............. ............. ............. *
* *
* Q20 = 3.0123 Q2+1= 0.0093 Q2+2= 2.9139 ............. ............. *
* Q2-1= 0.0046 Q2-2= -0.0065 ............. ............. *
* *
* Q30 = ZERO Q3+1= ZERO Q3+2= ZERO Q3+3= ZERO ............. *
* Q3-1= ZERO Q3-2= ZERO Q3-3= ZERO ............. *
* *
* Q40 = 0.0426 Q4+1= 0.0017 Q4+2= 0.1023 Q4+3= -0.0017 Q4+4= 0.0469 *
* Q4-1= -0.0047 Q4-2=-3.4E-04 Q4-3=-4.8E-04 Q4-4=-2.1E-04 *
* *
*******************************************************************************
*******************************************************************************
* *
* MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] PROTONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* M10 = 1.3488 M1+1= -0.0164 ............. ............. ............. *
* M1-1= 2.8124 ............. ............. ............. *
* *
* M20 = ZERO M2+1= ZERO M2+2= ZERO ............. ............. *
* M2-1= ZERO M2-2= ZERO ............. ............. *
* *
* M30 = 25.2153 M3+1= -0.3166 M3+2= 16.4334 M3+3= -0.4985 ............. *
* M3-1= 64.5581 M3-2= 0.2341 M3-3= 67.6302 ............. *
* *
*******************************************************************************
* *
* CARTESIAN COMPONENTS OF MAGNETIC DIPOLE MOMENT PROTONS *
* *
* M1_X = 0.0473 M1_Y = 8.1403 M1_Z = 2.7606 M1 = 8.5958 *
* *
*******************************************************************************
*******************************************************************************
* *
* ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS PROTONS *
* *
* R_RMS = 4.8172 X_RMS = 2.8862 Y_RMS = 2.3192 Z_RMS = 3.0817 *
* *
* R_GEO = 6.2190 X_GEO = 6.4537 Y_GEO = 5.1859 Z_GEO = 6.8909 *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) PROTONS *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* B10 = ZERO B1+1= ZERO ............. ............. ............. *
* B1-1= ZERO ............. ............. ............. *
* *
* B20 = 0.1805 B2+1= -0.0019 B2+2= 0.1746 ............. ............. *
* B2-1=-9.6E-04 B2-2=-3.9E-04 ............. ............. *
* *
* B30 = ZERO B3+1= ZERO B3+2= ZERO B3+3= ZERO ............. *
* B3-1= ZERO B3-2= ZERO B3-3= ZERO ............. *
* *
* B40 = 0.0209 B4+1= 0.0012 B4+2= 0.0710 B4+3= -0.0012 B4+4= 0.0326 *
* B4-1= -0.0032 B4-2=-2.4E-04 B4-3=-3.3E-04 B4-4=-1.5E-04 *
* *
*******************************************************************************
*******************************************************************************
* *
* DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL *
* *
*******************************************************************************
* *
* TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) *
* -------- ------ ------------ ------------ *
* DRHO_ = 15.470552 7.856029 15.470552 0.241506 *
* DRHOD = 11.025716 5.597383 11.025716 0.169051 *
* DLPR_ = -3.281263 -1.666589 -3.281263 -0.051916 *
* DTAU_ = 13.894330 7.121986 13.894330 0.349643 *
* DSCU_ = 0.141980 0.082666 0.141980 0.023351 *
* DDIV_ = 0.835435 0.432063 0.835435 0.028692 *
* *
* DSPI_ = 0.003283 0.002367 0.003283 0.001451 *
* DSPID = 0.002310 0.001653 0.002310 0.000996 *
* DLPS_ = -0.003385 -0.002590 -0.003385 -0.001794 *
* DCUR_ = 0.007178 0.005328 0.007178 0.003477 *
* DKIS_ = 0.003914 0.002640 0.003914 0.001366 *
* DROT_ = 0.001335 0.000797 0.001335 0.000260 *
* *
*******************************************************************************
*******************************************************************************
* *
* CONTRIBUTIONS TO ENERGY IN THE SKYRME FUNCTIONAL *
* *
*******************************************************************************
* *
* TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) *
* -------- ------ ------------ ------------ *
* ERHO_ = -21380.496528 6129.862975 -15344.853968 94.220415 *
* ERHOD = 14328.836772 -3637.133011 10746.627579 -54.923818 *
* ELPR_ = 279.932726 -57.028599 223.792372 -0.888244 *
* ETAU_ = 955.235161 -485.185312 481.959559 -11.909709 *
* ESCU_ = 0.000000 0.000000 0.000000 0.000000 *
* EDIV_ = -54.303245 -28.084109 -81.454868 -0.932487 *
* ============ ============ ============ ============ *
* SUM EVEN: -5870.795114 1922.431944 -3973.929326 25.566156 *
* *
* ESPI_ = -0.195402 1.565407 0.890165 0.479841 *
* ESPID = 0.000000 -1.074121 -0.750504 -0.323616 *
* ELPS_ = 0.000000 -0.088618 -0.057918 -0.030701 *
* ECUR_ = -0.493492 0.362950 -0.248989 0.118447 *
* EKIS_ = 0.000000 0.000000 0.000000 0.000000 *
* EROT_ = -0.086783 -0.051835 -0.130174 -0.008444 *
* ============ ============ ============ ============ *
* SUM ODD: -0.775677 0.713783 -0.297421 0.235526 *
* *
*******************************************************************************
*******************************************************************************
* *
* EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL *
* *
* ALPHA = -1.65924 BETA = 0.80185 GAMMA = 1.72034 *
* *
*******************************************************************************
*******************************************************************************
* *
* ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL *
* *
* R_RMS = 4.8635 X_RMS = 2.9230 Y_RMS = 2.3480 Z_RMS = 3.0978 *
* *
* R_GEO = 6.2788 X_GEO = 6.5361 Y_GEO = 5.2504 Z_GEO = 6.9270 *
* *
*******************************************************************************
*******************************************************************************
* *
* 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 *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) TOTAL *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* B10 = ZERO B1+1= ZERO ............. ............. ............. *
* B1-1= ZERO ............. ............. ............. *
* *
* B20 = 0.1721 B2+1= -0.0017 B2+2= 0.1759 ............. ............. *
* B2-1= 2.0E-04 B2-2=-3.7E-04 ............. ............. *
* *
* B30 = ZERO B3+1= ZERO B3+2= ZERO B3+3= ZERO ............. *
* B3-1= ZERO B3-2= ZERO B3-3= ZERO ............. *
* *
* B40 = 0.0173 B4+1= 0.0011 B4+2= 0.0635 B4+3= -0.0011 B4+4= 0.0287 *
* B4-1= -0.0029 B4-2=-2.3E-04 B4-3= 2.3E-04 B4-4=-1.2E-04 *
* *
*******************************************************************************
*******************************************************************************
* *
* 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 *
* *
*******************************************************************************
*******************************************************************************
* *
* BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) [INTRINSIC FRAME] TOTAL *
* *
* REAL PART FOR A NON-NEGATIVE PROJECTION *
* IMAGINARY PART FOR A NEGATIVE PROJECTION *
* *
*******************************************************************************
* *
* B10 = ZERO B1+1= ZERO ............. ............. ............. *
* B1-1= ZERO ............. ............. ............. *
* *
* B20 = 0.1721 B2+1= ZERO B2+2= 0.1759 ............. ............. *
* B2-1= ZERO B2-2= ZERO ............. ............. *
* *
* B30 = ZERO B3+1= ZERO B3+2= ZERO B3+3= ZERO ............. *
* B3-1= ZERO B3-2= ZERO B3-3= ZERO ............. *
* *
* B40 = 0.0173 B4+1= 2.7E-09 B4+2= 0.0635 B4+3=-1.8E-09 B4+4= 0.0287 *
* B4-1= -0.0030 B4-2= 1.9E-09 B4-3= 1.7E-04 B4-4=-3.8E-09 *
* *
*******************************************************************************
*******************************************************************************
* *
* 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 *
* *
*******************************************************************************
*******************************************************************************
* *
* CALCULATION OF THE DATA SET NO. 2 HAS BEEN SUCCESSFULLY COMPLETED *
* *
*******************************************************************************
*******************************************************************************
* *
* PROCESSOR EXECUTION TIMES IN SUBROUTINES *
* INCLUDING THEIR COMPLETE DOWN-CALLING TREES *
* *
*******************************************************************************
* *
* 1881 => HFODD 27 => POWALL 0 => NILSON 3 => INTKIN *
* *
* 38 => DIASIG 328 => DIAMAT 13 => AVANGY 30 => SPAVER *
* *
* 18 => AVOBSE 105 => INTMUL 2 => AVPARI 0 => AVSIMP *
* *
* 7 => NILABS 600 => DENSHF 194 => DENMAT 128 => MOMETS *
* *
* 0 => MOMSIF 0 => MOMVMU 10 => MAGMOM 7 => INTSPI *
* *
* 9 => SPIMOM 1 => LINAVR 1 => INTORB 248 => COUMAT *
* *
* 28 => BEGINT 247 => INTCOU 35 => SKFILD 0 => SHICOE *
* *
* 0 => SHIMOM 0 => SHISIF 0 => SHIMAG 0 => ROTMOM *
* *
* 13 => RECORD 201 => INTEGH 34 => INTMAS 12 => INTCEN *
* *
* 54 => INTSOR 294 => DIASIQ *
* *
*******************************************************************************