******************************************************************************* * * * HFODD HFODD HFODD HFODD HFODD HFODD HFODD HFODD * * * ******************************************************************************* * * * SKYRME-HARTREE-FOCK-BOGOLYUBOV CODE VERSION: 2.38J * * * * NO SYMMETRY-PLANES AND NO TIME-REVERSAL SYMMETRY * * * * DEFORMED CARTESIAN HARMONIC-OSCILLATOR BASIS * * * ******************************************************************************* * * * J. DOBACZEWSKI, B.G. CARLSSON, J. DUDEK, J. ENGEL * * P. OLBRATOWSKI, P. POWALOWSKI, M. SADZIAK, W. SATULA * * N. SCHUNCK, A. STASZCZAK, M. STOITSOV, M. ZALEWSKI * * AND H. ZDUNCZUK * * * * INSTYTUT FIZYKI TEORETYCZNEJ, WARSZAWA * * * * 1993-2009 * * * ******************************************************************************* ******************************************************************************* * * * CODE COMPILED WITH THE FOLLOWING ARRAY DIMENSIONS AND SWITCHES: * * * ******************************************************************************* * * * NDBASE = 680 NDSTAT = 181 NDXHRM = 30 NDYHRM = 30 NDZHRM = 30 * * * * NDMAIN = 14 NDMULT = 9 NDMULR = 4 NDLAMB = 9 NDITER = 5000 * * * * NDAKNO = 50 NDBKNO = 50 NDPROI = 50 NDCOUL = 80 NDPOLS = 25 * * * * IPARAL = 0 I_CRAY = 0 IZHPEV = 3 * * * ******************************************************************************* * * * EXECUTION BEGINS ON 2009.03.05 AT 09:20:29.528 * * * ******************************************************************************* ******************************************************************************* * * * CLASSICAL NUCLEAR SURFACE DEFINED FOR: N = 32 Z = 32 * * * ******************************************************************************* * * * SPHERICAL SHAPE REQUESTED FOR THE SURFACE * * * ******************************************************************************* * * * HOMEGA= 12.3000 FCHOM0= 1.2000 * * * * OSCILLATOR FREQUENCIES: HBAROX= 12.3000 HBAROY= 12.3000 HBAROZ= 12.3000 * * * * MOMENTS OF INERTIA: XMOMFC= 14.9370 YMOMFC= 14.9370 ZMOMFC= 14.9370 * * * * CENTERS OF MASS: CMSXFC= 0.0000 CMSYFC= 0.0000 CMSZFC= 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * SKYRME FORCE DEFINITION * * * ******************************************************************************* * * * PARAMETER SET SIII: T0= -1128.75 T1= 395.00 T2= -95.00 T3= 14000.00 * * * * POWER=1/1 W=120.000 X0= 0.45000 X1= 0.00000 X2= 0.00000 X3= 1.00000 * * * * ETA_J=0 ETA_W=0 ETACM=0 ETA_M=1 HBM= 20.73000 (SKYRME-FORCE SPIN-ORBIT) * * * ******************************************************************************* ******************************************************************************* * * * COUPLING CONSTANTS DEFINING THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * CRHO_ = -691.359375 536.156250 -423.281250 268.078125 * * CRHOD = 1750.000000 -1750.000000 875.000000 -875.000000 * * CLPR_ = -80.000000 34.062500 -62.968750 17.031250 * * CTAU_ = 75.000000 -61.250000 44.375000 -30.625000 * * CSCU_ = 0.000000 0.000000 0.000000 0.000000 * * CDIV_ = -60.000000 -60.000000 -90.000000 -30.000000 * * * * CSPI_ = -126.984375 282.187500 14.109375 141.093750 * * CSPID = 583.333333 -583.333333 291.666667 -291.666667 * * CLPS_ = 0.000000 34.062500 17.031250 17.031250 * * CCUR_ = -75.000000 61.250000 -44.375000 30.625000 * * CKIS_ = 0.000000 0.000000 0.000000 0.000000 * * CROT_ = -60.000000 -60.000000 -90.000000 -30.000000 * * * ******************************************************************************* ******************************************************************************* * * * PHYSICAL CONSTANTS: H_BARC=197.32891000 HBCOE2=137.03602000 * * * * XMASSN=938.90590000 XMASSP=938.27231000 * * * * HBMASS= 20.73000000 HBMRPA= 20.40609375 * * * * ECHAR2= 1.43997841 COULEX= -1.06350868 * * * ******************************************************************************* * * * OSCILLATOR LENGTHS: X= 1.8362302 Y= 1.8362302 Z= 1.8362302 * * * * OSCILLATOR CONSTANTS: X= 0.5445940 Y= 0.5445940 Z= 0.5445940 * * * * OSCILLATOR FREQUENCIES: X=12.3000000 Y=12.3000000 Z=12.3000000 * * * ******************************************************************************* * * * BASIS CUT-OFF CONTROL PARAMETERS: NXMAXX = 14 NYMAXX = 14 NZMAXX = 14 * * * * OPTIMUM NUMBERS OF GAUSS POINTS: NXHERM = 30 NYHERM = 30 NZHERM = 30 * * * * NLIMIT= 680 LDBASE= 680 MCOUNT= 3375 * * * * ENECUT= 800.0000 ELIMIT= 190.6500 * * * ******************************************************************************* ******************************************************************************* * * * SHAPE OF THE OSCILLATOR-BASIS DIAMOND (THREE CONSECUTIVE PROJECTIONS) * * * ******************************************************************************* * * * NX ===>>> 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 * * MAX.NY => 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------------------- * * NY= 0 (14) | 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 1 (13) | 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 2 (12) | 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 3 (11) | 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 4 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NY= 5 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NY= 6 ( 8) | 8 7 6 5 4 3 2 1 0 * * NY= 7 ( 7) | 7 6 5 4 3 2 1 0 * * NY= 8 ( 6) | 6 5 4 3 2 1 0 * * NY= 9 ( 5) | 5 4 3 2 1 0 * * NY=10 ( 4) | 4 3 2 1 0 * * NY=11 ( 3) | 3 2 1 0 * * NY=12 ( 2) | 2 1 0 * * NY=13 ( 1) | 1 0 * * NY=14 ( 0) | 0 * * * * * * NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 * * MAX.NY => 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------------------- * * NY= 0 (14) | 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 1 (13) | 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 2 (12) | 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 3 (11) | 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 4 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NY= 5 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NY= 6 ( 8) | 8 7 6 5 4 3 2 1 0 * * NY= 7 ( 7) | 7 6 5 4 3 2 1 0 * * NY= 8 ( 6) | 6 5 4 3 2 1 0 * * NY= 9 ( 5) | 5 4 3 2 1 0 * * NY=10 ( 4) | 4 3 2 1 0 * * NY=11 ( 3) | 3 2 1 0 * * NY=12 ( 2) | 2 1 0 * * NY=13 ( 1) | 1 0 * * NY=14 ( 0) | 0 * * * * * * NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 * * MAX.NX => 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------------------- * * NX= 0 (14) | 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 1 (13) | 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 2 (12) | 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 3 (11) | 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 4 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NX= 5 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NX= 6 ( 8) | 8 7 6 5 4 3 2 1 0 * * NX= 7 ( 7) | 7 6 5 4 3 2 1 0 * * NX= 8 ( 6) | 6 5 4 3 2 1 0 * * NX= 9 ( 5) | 5 4 3 2 1 0 * * NX=10 ( 4) | 4 3 2 1 0 * * NX=11 ( 3) | 3 2 1 0 * * NX=12 ( 2) | 2 1 0 * * NX=13 ( 1) | 1 0 * * NX=14 ( 0) | 0 * * * ******************************************************************************* ******************************************************************************* * * * AUTOMATED HANDLING OF THE COULOMB FILE * * * ******************************************************************************* * * * THE COULOMB FILE WRITTEN ON THE DISC, FILE NAME (FILCOU): * * ge064-a.cou * * * ******************************************************************************* * * * NUCLIDE: N = 32 Z = 32 NUMBER OF ITERATIONS = 20 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 LIPKIN-NOGAMI = 0.50 * * * ******************************************************************************* * * * DIRECT COULOMB TERMS SOLVED BY GREEN FUNCTION AND FOURIER TRANSFORM METHOD * * COULOMB PARAMETERS: NUMCOU = 80 NUMETA = 79 BOUCOU = 20.0 FURMAX =0.25 * * * ******************************************************************************* * * * EXCHANGE COULOMB TERMS TREATED BY THE SLATER APPROXIMATION * * * ******************************************************************************* * * * MAXIMUM NUMBERS OF MULTIPOLE MOMENTS USED FOR THE CONSTRAINTS, NMUCON=2 * * FOR THE COULOMB FIELD, NMUCOU=4 * * FOR THE OUTPUT INFO, NMUPRI=6 * * * ******************************************************************************* * * * MULTIPOLE CONSTRAINTS: LAMBDA= 2 MIU= 0 MOMENT= 2.700 STIFFNESS= 1.000 * * LAMBDA= 2 MIU= 2 MOMENT= -1.300 STIFFNESS= 1.000 * * * ******************************************************************************* * * * 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: NO * * PROTONS : NO * * MASS : YES * * * ******************************************************************************* * * * PRINTING THE MOMENTS, MOMENTA AND DEFORMATIONS IN THE INTRINSIC FRAME: NO * * * ******************************************************************************* * * * NILSSON LABELS OF SINGLE-PARTICLE STATES SHOWN WITH RESPECT TO THE: Z-AXIS * * * ******************************************************************************* * * * NUMBERS OF LEVELS IN TIME-REVERSAL/CHARGE BLOCKS: TIME-UP TIME-DOWN * * NEUTRONS: 32 32 * * PROTONS: 32 32 * * * ******************************************************************************* * * * CALCULATIONS WITHOUT PAIRING * * * ******************************************************************************* * * * CALCULATIONS WITH PARITY/SIGNATURE SYMMETRY * * * ******************************************************************************* * * * CALCULATIONS WITH TIME-ODD SYMMETRIES: T*SIMPLEX_X => CONSERVED * * T*SIMPLEX_Y => CONSERVED * * T*SIMPLEX_Z => CONSERVED * * * ******************************************************************************* * * * CALCULATIONS WITHOUT ROTATION * * * ******************************************************************************* ******************************************************************************* * * * PARITY/SIGNATURE CONFIGURATIONS: * * * * V A C U U M P A R T I C L E S H O L E S * * =========== ================= ========= * * (++) (+-) (-+) (--) (++) (+-) (-+) (--) (++) (+-) (-+) (--) * * * * NEUTRONS: 7 7 9 9 0 0 0 0 0 0 0 0 * * PROTONS : 7 7 9 9 0 0 0 0 0 0 0 0 * * * ******************************************************************************* ******************************************************************************* * * * PARAMETERS OF THE NILSSON MODEL NEUTRONS * * * ******************************************************************************* * * * ESCALE =63.350 KAPPA = 0.080 MIU = 0.000 * * * * HOMEGA =12.300 CNILSS =-1.968 DNILSS = 0.000 * * * * OSCILLATOR FREQUENCIES: HBAROX =12.300 HBAROY =12.300 HBAROZ =12.300 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: NILSSON NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 2) -34.568 ( 0 ----- 1) | 1, 1, 0, 1/2> -100 -1.000 -0.333 0.333 * * 3) -34.568 ( 0 ----- 2) | 1, 0, 1, 3/2> -100 0.000 0.000 0.333 * * 4) -28.664 ( 0 ----- 3) | 1, 0, 1, 1/2> -100 0.500 -0.167 -0.333 * * 5) -24.236 ( 2 ----- 0) | 2, 2, 0, 1/2> 100 -1.212 -0.242 0.200 * * 6) -24.236 ( 3 ----- 0) | 2, 0, 2, 5/2> 100 -0.069 -0.014 0.200 * * 7) -24.236 ( 4 ----- 0) | 2, 1, 1, 3/2> 100 -0.220 -0.044 0.200 * * 8) -20.300 ( 5 ----- 0) | 2, 0, 0, 1/2> 100 -0.500 -0.500 1.000 * * 9) -14.396 ( 6 ----- 0) | 2, 0, 2, 3/2> 100 -0.348 0.070 -0.200 * * 10) -14.396 ( 7 ----- 0) | 2, 0, 2, 3/2> 100 1.348 -0.270 -0.200 * * 11) -13.904 ( 0 ----- 4) | 3, 1, 2, 5/2> -100 0.177 0.025 0.143 * * 12) -13.904 ( 0 ----- 5) | 3, 2, 1, 3/2> -100 -1.014 -0.145 0.143 * * 13) -13.904 ( 0 ----- 6) | 3, 0, 3, 7/2> -100 0.575 0.082 0.143 * * 14) -13.904 ( 0 ----- 7) | 3, 1, 0, 1/2> -100 -1.739 -0.248 0.143 * * 15) -9.968 ( 0 ----- 8) | 3, 3, 0, 1/2> -100 -1.003 -0.334 0.333 * * 16) -9.968 ( 0 ----- 9) | 3, 0, 1, 3/2> -100 0.003 0.001 0.333 * * 17) -4.064 ( 0 ----- 10) | 3, 0, 1, 1/2> -100 0.500 -0.167 -0.333 * * 18) -3.572 ( 8 ----- 0) | 4, 0, 4, 9/2> 100 -1.505 -0.167 0.111 * * 19) -3.572 ( 9 ----- 0) | 4, 0, 4, 9/2> 100 -0.363 -0.040 0.111 * * 20) -3.572 (10 ----- 0) | 4, 2, 2, 5/2> 100 1.216 0.135 0.111 * * 21) -3.572 (11 ----- 0) | 4, 2, 2, 5/2> 100 -2.233 -0.248 0.111 * * 22) -3.572 (12 ----- 0) | 4, 1, 3, 7/2> 100 0.385 0.043 0.111 * * 23) -0.128 ( 0 ----- 11) | 3, 2, 1, 1/2> -100 0.219 -0.031 -0.143 * * 24) -0.128 ( 0 ----- 12) | 3, 1, 2, 3/2> -100 2.288 -0.327 -0.143 * * 25) -0.128 ( 0 ----- 13) | 3, 0, 3, 5/2> -100 -1.007 0.144 -0.143 * * * ******************************************************************************* ******************************************************************************* * * * PARAMETERS OF THE NILSSON MODEL PROTONS * * * ******************************************************************************* * * * ESCALE =63.350 KAPPA = 0.080 MIU = 0.000 * * * * HOMEGA =12.300 CNILSS =-1.968 DNILSS = 0.000 * * * * OSCILLATOR FREQUENCIES: HBAROX =12.300 HBAROY =12.300 HBAROZ =12.300 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: NILSSON PROTONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 2) -34.568 ( 0 ----- 1) | 1, 1, 0, 1/2> -100 -1.000 -0.333 0.333 * * 3) -34.568 ( 0 ----- 2) | 1, 0, 1, 3/2> -100 0.000 0.000 0.333 * * 4) -28.664 ( 0 ----- 3) | 1, 0, 1, 1/2> -100 0.500 -0.167 -0.333 * * 5) -24.236 ( 2 ----- 0) | 2, 2, 0, 1/2> 100 -1.212 -0.242 0.200 * * 6) -24.236 ( 3 ----- 0) | 2, 0, 2, 5/2> 100 -0.069 -0.014 0.200 * * 7) -24.236 ( 4 ----- 0) | 2, 1, 1, 3/2> 100 -0.220 -0.044 0.200 * * 8) -20.300 ( 5 ----- 0) | 2, 0, 0, 1/2> 100 -0.500 -0.500 1.000 * * 9) -14.396 ( 6 ----- 0) | 2, 0, 2, 3/2> 100 -0.348 0.070 -0.200 * * 10) -14.396 ( 7 ----- 0) | 2, 0, 2, 3/2> 100 1.348 -0.270 -0.200 * * 11) -13.904 ( 0 ----- 4) | 3, 1, 2, 5/2> -100 0.177 0.025 0.143 * * 12) -13.904 ( 0 ----- 5) | 3, 2, 1, 3/2> -100 -1.014 -0.145 0.143 * * 13) -13.904 ( 0 ----- 6) | 3, 0, 3, 7/2> -100 0.575 0.082 0.143 * * 14) -13.904 ( 0 ----- 7) | 3, 1, 0, 1/2> -100 -1.739 -0.248 0.143 * * 15) -9.968 ( 0 ----- 8) | 3, 3, 0, 1/2> -100 -1.003 -0.334 0.333 * * 16) -9.968 ( 0 ----- 9) | 3, 0, 1, 3/2> -100 0.003 0.001 0.333 * * 17) -4.064 ( 0 ----- 10) | 3, 0, 1, 1/2> -100 0.500 -0.167 -0.333 * * 18) -3.572 ( 8 ----- 0) | 4, 0, 4, 9/2> 100 -1.505 -0.167 0.111 * * 19) -3.572 ( 9 ----- 0) | 4, 0, 4, 9/2> 100 -0.363 -0.040 0.111 * * 20) -3.572 (10 ----- 0) | 4, 2, 2, 5/2> 100 1.216 0.135 0.111 * * 21) -3.572 (11 ----- 0) | 4, 2, 2, 5/2> 100 -2.233 -0.248 0.111 * * 22) -3.572 (12 ----- 0) | 4, 1, 3, 7/2> 100 0.385 0.043 0.111 * * 23) -0.128 ( 0 ----- 11) | 3, 2, 1, 1/2> -100 0.219 -0.031 -0.143 * * 24) -0.128 ( 0 ----- 12) | 3, 1, 2, 3/2> -100 2.288 -0.327 -0.143 * * 25) -0.128 ( 0 ----- 13) | 3, 0, 3, 5/2> -100 -1.007 0.144 -0.143 * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 9.803645 4.901823 9.803645 0.000000 * * DRHOD = 1.919966 0.959983 1.919966 0.000000 * * DLPR_ = -3.120921 -1.560461 -3.120921 0.000000 * * DTAU_ = 11.048155 5.524077 11.048155 0.000000 * * DSCU_ = 0.290743 0.145372 0.290743 0.000000 * * DDIV_ = 1.163843 0.581921 1.163843 0.000000 * * * * DSPI_ = 0.000000 0.000000 0.000000 0.000000 * * DSPID = 0.000000 0.000000 0.000000 0.000000 * * DLPS_ = 0.000000 0.000000 0.000000 0.000000 * * DCUR_ = 0.000000 0.000000 0.000000 0.000000 * * DKIS_ = 0.000000 0.000000 0.000000 0.000000 * * DROT_ = 0.000000 0.000000 0.000000 0.000000 * * * ******************************************************************************* ******************************************************************************* ******************************************************************************* * * * CONTRIBUTIONS TO ENERGY IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * ERHO_ = -6777.842174 2628.142884 -4149.699290 0.000000 * * ERHOD = 3359.941157 -1679.970579 1679.970579 0.000000 * * ELPR_ = 249.673689 -53.153188 196.520501 0.000000 * * ETAU_ = 828.611610 -338.349741 490.261870 0.000000 * * ESCU_ = 0.000000 0.000000 0.000000 0.000000 * * EDIV_ = -69.830572 -34.915286 -104.745858 0.000000 * * ============ ============ ============ ============ * * SUM EVEN: -2409.446290 521.754091 -1887.692199 0.000000 * * * * ESPI_ = 0.000000 0.000000 0.000000 0.000000 * * ESPID = 0.000000 0.000000 0.000000 0.000000 * * ELPS_ = 0.000000 0.000000 0.000000 0.000000 * * ECUR_ = 0.000000 0.000000 0.000000 0.000000 * * EKIS_ = 0.000000 0.000000 0.000000 0.000000 * * EROT_ = 0.000000 0.000000 0.000000 0.000000 * * ============ ============ ============ ============ * * SUM ODD: 0.000000 0.000000 0.000000 0.000000 * * * ******************************************************************************* ******************************************************************************* * * * EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL * * * * ALPHA = 0.00000 BETA = 86.81452 GAMMA = 90.00000 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 3.4653 X_RMS = 2.0007 Y_RMS = 2.0007 Z_RMS = 2.0007 * * * * R_GEO = 4.4737 X_GEO = 4.4737 Y_GEO = 4.4737 Z_GEO = 4.4737 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * ******************************************************************************* * * * Q00 = 64.0000 ............. ............. ............. ............. * * * * Q10 = ZERO Q11 = ZERO ............. ............. ............. * * * * Q20 = ZERO Q21 = ZERO Q22 = ZERO ............. ............. * * * * Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. * * * * Q40 = ZERO Q41 = ZERO Q42 = ZERO Q43 = ZERO Q44 = ZERO * * * * Q50 = ZERO Q51 = ZERO Q52 = ZERO Q53 = ZERO Q54 = ZERO * * * * Q55 = ZERO * * * * Q60 = ZERO Q61 = ZERO Q62 = ZERO Q63 = ZERO Q64 = ZERO * * * * Q65 = ZERO Q66 = ZERO * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (EXACT MULTIPOLE MOMENTS) TOTAL * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = ZERO B21 = ZERO B22 = ZERO ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 = ZERO B41 = ZERO B42 = ZERO B43 = ZERO B44 = ZERO * * * * B50 = ZERO B51 = ZERO B52 = ZERO B53 = ZERO B54 = ZERO * * * * B55 = ZERO * * * * B60 = ZERO B61 = ZERO B62 = ZERO B63 = ZERO B64 = ZERO * * * * B65 = ZERO B66 = ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1706 FM^(-3) BOHR RADIUS = 4.4737 FM * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) TOTAL * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = ZERO B21 = ZERO B22 = ZERO ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 = ZERO B41 = ZERO B42 = ZERO B43 = ZERO B44 = ZERO * * * * B50 = ZERO B51 = ZERO B52 = ZERO B53 = ZERO B54 = ZERO * * * * B55 = ZERO * * * * B60 = ZERO B61 = ZERO B62 = ZERO B63 = ZERO B64 = ZERO * * * * B65 = ZERO B66 = ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1706 FM^(-3) BOHR RADIUS = 4.4737 FM * * * ******************************************************************************* ******************************************************************************* * * * NEUTRON CONFIGURATIONS * * ======================= * * P S 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 * * --- -------------------------------------------------------------- * * * * CONF: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * PROTON CONFIGURATIONS * * ======================= * * P S 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 * * --- -------------------------------------------------------------- * * * * CONF: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * ENERGIES (MEV) * * * ******************************************************************************* * * * KINETIC: (NEU)= 690.145312 (PRO)= 690.145312 (TOT)= 1380.290625 * * SUM EPS: (NEU)= -680.104000 (PRO)= -680.104000 (TOT)= -1360.208000 * * 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)= -7.016000 (PRO)= -7.016000 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 201.316506 (EXC)= -13.777065 (TOT)= 187.539441 * * * * CONSTR. (MULT)= 8.980000 SLOPE= 2.800000 CORR.= 0.000000 * * 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= 839.985289 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -310.882133 * * * * SPIN-ORB (EVE)= -104.745858 (ODD)= 0.000000 (TOT)= -104.745858 * * SKYRME: (EVE)= -1887.692199 (ODD)= 0.000000 (TOT)= -1887.692199 * * * * TOTAL: (STAB)= -514.674199 (SP)= -834.536332 (FUN)= -319.862133 * * * ******************************************************************************* ******************************************************************************* * * * CONVERGENCE REPORT * * * ******************************************************************************* * * * ITER ENERGY STABILITY Q_2 GAMMA SPIN OMEGA ANGLE EPAIR * * * * 0 -319.862133 -514.674199 0.000 41.34 0.00 0.000 0.000 0.00 * * 1 -370.533448 73.232244 2.039 18.79 0.00 0.000 0.000 0.00 * * 2 -485.150276 158.398451 1.194 18.38 0.00 0.000 0.000 0.00 * * 3 -504.665151 97.236471 1.669 15.29 0.00 0.000 0.000 0.00 * * 4 -517.027347 56.748154 1.821 16.05 0.00 0.000 0.000 0.00 * * 5 -526.452482 32.524813 2.077 17.33 0.00 0.000 0.000 0.00 * * 6 -533.127076 15.918449 2.296 19.24 0.00 0.000 0.000 0.00 * * 7 -536.949965 3.884326 2.450 20.38 0.00 0.000 0.000 0.00 * * 8 -538.766756 -2.110647 2.566 21.27 0.00 0.000 0.000 0.00 * * 9 -539.578809 -4.043832 2.656 21.96 0.00 0.000 0.000 0.00 * * 10 -539.953187 -3.912343 2.728 22.51 0.00 0.000 0.000 0.00 * * 11 -540.139284 -3.136869 2.786 22.97 0.00 0.000 0.000 0.00 * * 12 -540.240434 -2.283317 2.833 23.35 0.00 0.000 0.000 0.00 * * 13 -540.300015 -1.588341 2.871 23.67 0.00 0.000 0.000 0.00 * * 14 -540.337381 -1.071092 2.901 23.95 0.00 0.000 0.000 0.00 * * 15 -540.361848 -0.714356 2.925 24.18 0.00 0.000 0.000 0.00 * * 16 -540.378345 -0.473085 2.945 24.39 0.00 0.000 0.000 0.00 * * 17 -540.389687 -0.315025 2.961 24.56 0.00 0.000 0.000 0.00 * * 18 -540.397597 -0.211369 2.974 24.71 0.00 0.000 0.000 0.00 * * 19 -540.403172 -0.144237 2.985 24.84 0.00 0.000 0.000 0.00 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 3) -39.302 ( 0 ----- 2) | 1, 0, 1, 3/2> -100 0.137 0.295 2.144 * * 4) -36.870 ( 0 ----- 3) | 1, 0, 1, 1/2> -100 0.045 -0.362 -7.962 * * 5) -31.038 ( 2 ----- 0) | 2, 2, 0, 1/2> 100 -0.440 -0.155 0.353 * * 6) -29.978 ( 3 ----- 0) | 2, 1, 1, 3/2> 100 -0.748 -0.050 0.067 * * 7) -28.059 ( 4 ----- 0) | 2, 0, 2, 5/2> 100 -0.031 -0.077 2.451 * * 8) -26.499 ( 5 ----- 0) | 2, 1, 1, 1/2> 100 0.433 -0.311 -0.718 * * 9) -24.082 ( 6 ----- 0) | 2, 0, 2, 3/2> 100 -0.280 0.012 -0.044 * * 10) -21.262 ( 7 ----- 0) | 2, 0, 0, 1/2> 100 0.072 -0.420 -5.830 * * 11) -20.932 ( 0 ----- 4) | 3, 2, 1, 1/2> -100 -0.062 -0.005 0.080 * * 12) -19.735 ( 0 ----- 5) | 3, 2, 1, 3/2> -100 -1.534 -0.220 0.144 * * 13) -18.718 ( 0 ----- 6) | 3, 1, 2, 5/2> -100 -0.153 -0.057 0.371 * * 14) -16.551 ( 0 ----- 7) | 3, 0, 3, 7/2> -100 -0.005 0.011 -2.115 * * 15) -14.815 ( 0 ----- 8) | 3, 2, 1, 1/2> -100 0.404 -0.363 -0.898 * * 16) -13.687 ( 0 ----- 9) | 3, 1, 2, 3/2> -100 -0.388 0.063 -0.162 * * 17) -11.051 ( 0 ----- 10) | 3, 1, 0, 1/2> -100 0.220 -0.371 -1.691 * * 18) -10.785 ( 0 ----- 11) | 3, 0, 1, 3/2> -100 0.453 0.277 0.611 * * 19) -10.593 ( 8 ----- 0) | 4, 3, 1, 1/2> 100 0.047 0.011 0.242 * * 20) -9.239 ( 9 ----- 0) | 4, 3, 1, 3/2> 100 -1.101 -0.136 0.124 * * 21) -9.048 ( 0 ----- 12) | 3, 0, 3, 5/2> -100 -0.314 -0.005 0.016 * * 22) -8.537 (10 ----- 0) | 4, 2, 2, 5/2> 100 -1.125 -0.148 0.132 * * 23) -7.191 (11 ----- 0) | 4, 1, 3, 7/2> 100 -0.041 0.005 -0.115 * * 24) -6.708 ( 0 ----- 13) | 3, 0, 1, 1/2> -100 0.378 -0.330 -0.872 * * 25) -4.742 (12 ----- 0) | 4, 0, 4, 9/2> 100 0.000 -0.003 7.718 * * 26) -3.916 (13 ----- 0) | 4, 4, 0, 1/2> 100 -0.228 -0.365 1.601 * * 27) -2.892 (14 ----- 0) | 4, 3, 1, 1/2> 100 0.187 0.245 1.312 * * 28) -1.199 (15 ----- 0) | 4, 1, 1, 3/2> 100 -0.129 -0.041 0.318 * * 29) -0.604 (16 ----- 0) | 4, 2, 2, 3/2> 100 0.248 -0.371 -1.495 * * 30) -0.053 ( 0 ----- 14) | 5, 3, 0, 1/2> -100 0.020 0.002 0.126 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 1) -38.419 ( 1 ----- 0) | 0, 0, 0, 1/2> 100 -0.498 -0.500 1.004 * * 2) -30.196 ( 0 ----- 1) | 1, 1, 0, 1/2> -100 -0.675 -0.436 0.645 * * 3) -28.761 ( 0 ----- 2) | 1, 0, 1, 3/2> -100 0.152 0.304 2.003 * * 4) -26.324 ( 0 ----- 3) | 1, 0, 1, 1/2> -100 0.029 -0.368 -1.E+01 * * 5) -20.775 ( 2 ----- 0) | 2, 2, 0, 1/2> 100 -0.405 -0.142 0.351 * * 6) -19.715 ( 3 ----- 0) | 2, 1, 1, 3/2> 100 -0.778 -0.057 0.074 * * 7) -17.761 ( 4 ----- 0) | 2, 0, 2, 5/2> 100 -0.042 -0.076 1.809 * * 8) -16.258 ( 5 ----- 0) | 2, 1, 1, 1/2> 100 0.449 -0.311 -0.692 * * 9) -13.812 ( 6 ----- 0) | 2, 0, 2, 3/2> 100 -0.257 0.014 -0.055 * * 10) -10.983 ( 7 ----- 0) | 2, 0, 0, 1/2> 100 0.039 -0.428 -1.E+01 * * 11) -10.939 ( 0 ----- 4) | 3, 2, 1, 1/2> -100 -0.042 0.000 0.003 * * 12) -9.732 ( 0 ----- 5) | 3, 2, 1, 3/2> -100 -1.513 -0.218 0.144 * * 13) -8.730 ( 0 ----- 6) | 3, 1, 2, 5/2> -100 -0.201 -0.062 0.309 * * 14) -6.546 ( 0 ----- 7) | 3, 0, 3, 7/2> -100 -0.005 0.011 -2.216 * * 15) -4.866 ( 0 ----- 8) | 3, 2, 1, 1/2> -100 0.442 -0.352 -0.797 * * 16) -3.779 ( 0 ----- 9) | 3, 1, 2, 3/2> -100 -0.410 0.050 -0.121 * * 17) -1.202 ( 0 ----- 10) | 3, 1, 2, 3/2> -100 0.251 -0.379 -1.514 * * 18) -0.936 ( 0 ----- 11) | 3, 0, 1, 3/2> -100 0.438 0.288 0.659 * * 19) -0.904 ( 8 ----- 0) | 4, 3, 1, 1/2> 100 0.051 0.012 0.238 * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 6.957660 3.479917 6.957660 0.002174 * * DRHOD = 0.890601 0.445440 0.890601 0.000278 * * DLPR_ = -2.217658 -1.109739 -2.217658 -0.001821 * * DTAU_ = 5.678356 2.840471 5.678356 0.002587 * * DSCU_ = 0.061889 0.030971 0.061889 0.000052 * * DDIV_ = 0.387166 0.193731 0.387166 0.000296 * * * * DSPI_ = 0.000000 0.000000 0.000000 0.000000 * * DSPID = 0.000000 0.000000 0.000000 0.000000 * * DLPS_ = 0.000000 0.000000 0.000000 0.000000 * * DCUR_ = 0.000000 0.000000 0.000000 0.000000 * * DKIS_ = 0.000000 0.000000 0.000000 0.000000 * * DROT_ = 0.000000 0.000000 0.000000 0.000000 * * * ******************************************************************************* ******************************************************************************* ******************************************************************************* * * * CONTRIBUTIONS TO ENERGY IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * ERHO_ = -4810.243404 1865.779261 -2945.046982 0.582839 * * ERHOD = 1558.551814 -779.519438 779.275907 -0.243532 * * ELPR_ = 177.412615 -37.800493 139.643132 -0.031010 * * ETAU_ = 425.876679 -173.978859 251.977035 -0.079215 * * ESCU_ = 0.000000 0.000000 0.000000 0.000000 * * EDIV_ = -23.229930 -11.623848 -34.844895 -0.008883 * * ============ ============ ============ ============ * * SUM EVEN: -2671.632226 862.856623 -1808.995803 0.220199 * * * * ESPI_ = 0.000000 0.000000 0.000000 0.000000 * * ESPID = 0.000000 0.000000 0.000000 0.000000 * * ELPS_ = 0.000000 0.000000 0.000000 0.000000 * * ECUR_ = 0.000000 0.000000 0.000000 0.000000 * * EKIS_ = 0.000000 0.000000 0.000000 0.000000 * * EROT_ = 0.000000 0.000000 0.000000 0.000000 * * ============ ============ ============ ============ * * SUM ODD: 0.000000 0.000000 0.000000 0.000000 * * * ******************************************************************************* ******************************************************************************* * * * EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL * * * * ALPHA =270.00000 BETA = 0.00000 GAMMA = 0.00000 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 3.9235 X_RMS = 1.9634 Y_RMS = 2.2347 Z_RMS = 2.5583 * * * * R_GEO = 5.0653 X_GEO = 4.3904 Y_GEO = 4.9970 Z_GEO = 5.7206 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * ******************************************************************************* * * * Q00 = 64.0000 ............. ............. ............. ............. * * * * Q10 = ZERO Q11 = ZERO ............. ............. ............. * * * * Q20 = 2.7142 Q21 = ZERO Q22 = -1.2626 ............. ............. * * * * Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. * * * * Q40 = -0.0046 Q41 = ZERO Q42 = -0.0408 Q43 = ZERO Q44 = 0.0017 * * * * Q50 = ZERO Q51 = ZERO Q52 = ZERO Q53 = ZERO Q54 = ZERO * * * * Q55 = ZERO * * * * Q60 = -0.0032 Q61 = ZERO Q62 = 3.8E-04 Q63 = ZERO Q64 =-4.4E-04 * * * * Q65 = ZERO Q66 = 2.2E-05 * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (EXACT MULTIPOLE MOMENTS) TOTAL * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = 0.2167 B21 = ZERO B22 = -0.1186 ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 = -0.0273 B41 = ZERO B42 = -0.0358 B43 = ZERO B44 = -0.0045 * * * * B50 = ZERO B51 = ZERO B52 = ZERO B53 = ZERO B54 = ZERO * * * * B55 = ZERO * * * * B60 = -0.0067 B61 = ZERO B62 = 0.0155 B63 = ZERO B64 = -0.0104 * * * * B65 = ZERO B66 =-3.0E-04 * * * ******************************************************************************* * * * BOHR DENSITY = 0.1220 FM^(-3) BOHR RADIUS = 4.9785 FM * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) TOTAL * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = 0.2184 B21 = ZERO B22 = -0.1016 ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 = -0.0046 B41 = ZERO B42 = -0.0573 B43 = ZERO B44 = 0.0025 * * * * B50 = ZERO B51 = ZERO B52 = ZERO B53 = ZERO B54 = ZERO * * * * B55 = ZERO * * * * B60 = -0.0124 B61 = ZERO B62 = 0.0021 B63 = ZERO B64 = -0.0024 * * * * B65 = ZERO B66 = 1.2E-04 * * * ******************************************************************************* * * * BOHR DENSITY = 0.1176 FM^(-3) BOHR RADIUS = 5.0653 FM * * * ******************************************************************************* ******************************************************************************* * * * NEUTRON CONFIGURATIONS * * ======================= * * P S 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 * * --- -------------------------------------------------------------- * * * * CONF: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * PROTON CONFIGURATIONS * * ======================= * * P S 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 * * --- -------------------------------------------------------------- * * * * CONF: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * ENERGIES (MEV) * * * ******************************************************************************* * * * KINETIC: (NEU)= 560.348378 (PRO)= 544.035786 (TOT)= 1104.384164 * * SUM EPS: (NEU)= -862.957676 (PRO)= -535.192516 (TOT)= -1398.150192 * * 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)= -12.368983 (PRO)= -2.490693 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 176.334811 (EXC)= -12.350511 (TOT)= 163.984300 * * * * CONSTR. (MULT)= 0.001604 SLOPE= -0.103337 CORR.= -0.008673 * * 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= 389.516188 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -540.405535 * * * * SPIN-ORB (EVE)= -34.853778 (ODD)= 0.000000 (TOT)= -34.853778 * * SKYRME: (EVE)= -1808.775604 (ODD)= 0.000000 (TOT)= -1808.775604 * * * * TOTAL: (STAB)= -0.100226 (SP)= -540.507365 (FUN)= -540.407139 * * * ******************************************************************************* ******************************************************************************* * * * CALLING SUBROUTINE "RECORD" TO SAVE THE RECORD FILE CONTAINING RESULTS * * ge064-a.rec * * * ******************************************************************************* ******************************************************************************* * * * WARNING! THE STABILITY REQUIRED WAS NOT ACHIEVED. NO-OF-ITERATIONS EXCEEDED * * WARNING! THE STABILITY REQUIRED WAS NOT ACHIEVED. NO-OF-ITERATIONS EXCEEDED * * * ******************************************************************************* * * ******************************************************************************* * * * CALCULATION OF THE DATA SET NO. 1 HAS BEEN SUCCESSFULLY COMPLETED * * * ******************************************************************************* ******************************************************************************* * * * CLASSICAL NUCLEAR SURFACE DEFINED FOR: N = 32 Z = 32 * * * ******************************************************************************* * * * SPHERICAL SHAPE REQUESTED FOR THE SURFACE * * * ******************************************************************************* * * * HOMEGA= 12.3000 FCHOM0= 1.2000 * * * * OSCILLATOR FREQUENCIES: HBAROX= 12.3000 HBAROY= 12.3000 HBAROZ= 12.3000 * * * * MOMENTS OF INERTIA: XMOMFC= 14.9370 YMOMFC= 14.9370 ZMOMFC= 14.9370 * * * * CENTERS OF MASS: CMSXFC= 0.0000 CMSYFC= 0.0000 CMSZFC= 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * SKYRME FORCE DEFINITION * * * ******************************************************************************* * * * PARAMETER SET SIII: T0= -1128.75 T1= 395.00 T2= -95.00 T3= 14000.00 * * * * POWER=1/1 W=120.000 X0= 0.45000 X1= 0.00000 X2= 0.00000 X3= 1.00000 * * * * ETA_J=0 ETA_W=0 ETACM=0 ETA_M=1 HBM= 20.73000 (SKYRME-FORCE SPIN-ORBIT) * * * ******************************************************************************* ******************************************************************************* * * * COUPLING CONSTANTS DEFINING THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * CRHO_ = -691.359375 536.156250 -423.281250 268.078125 * * CRHOD = 1750.000000 -1750.000000 875.000000 -875.000000 * * CLPR_ = -80.000000 34.062500 -62.968750 17.031250 * * CTAU_ = 75.000000 -61.250000 44.375000 -30.625000 * * CSCU_ = 0.000000 0.000000 0.000000 0.000000 * * CDIV_ = -60.000000 -60.000000 -90.000000 -30.000000 * * * * CSPI_ = -126.984375 282.187500 14.109375 141.093750 * * CSPID = 583.333333 -583.333333 291.666667 -291.666667 * * CLPS_ = 0.000000 34.062500 17.031250 17.031250 * * CCUR_ = -75.000000 61.250000 -44.375000 30.625000 * * CKIS_ = 0.000000 0.000000 0.000000 0.000000 * * CROT_ = -60.000000 -60.000000 -90.000000 -30.000000 * * * ******************************************************************************* ******************************************************************************* * * * PHYSICAL CONSTANTS: H_BARC=197.32891000 HBCOE2=137.03602000 * * * * XMASSN=938.90590000 XMASSP=938.27231000 * * * * HBMASS= 20.73000000 HBMRPA= 20.40609375 * * * * ECHAR2= 1.43997841 COULEX= -1.06350868 * * * ******************************************************************************* * * * OSCILLATOR LENGTHS: X= 1.8362302 Y= 1.8362302 Z= 1.8362302 * * * * OSCILLATOR CONSTANTS: X= 0.5445940 Y= 0.5445940 Z= 0.5445940 * * * * OSCILLATOR FREQUENCIES: X=12.3000000 Y=12.3000000 Z=12.3000000 * * * ******************************************************************************* * * * BASIS CUT-OFF CONTROL PARAMETERS: NXMAXX = 14 NYMAXX = 14 NZMAXX = 14 * * * * OPTIMUM NUMBERS OF GAUSS POINTS: NXHERM = 30 NYHERM = 30 NZHERM = 30 * * * * NLIMIT= 680 LDBASE= 680 MCOUNT= 3375 * * * * ENECUT= 800.0000 ELIMIT= 190.6500 * * * ******************************************************************************* ******************************************************************************* * * * SHAPE OF THE OSCILLATOR-BASIS DIAMOND (THREE CONSECUTIVE PROJECTIONS) * * * ******************************************************************************* * * * NX ===>>> 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 * * MAX.NY => 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------------------- * * NY= 0 (14) | 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 1 (13) | 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 2 (12) | 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 3 (11) | 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 4 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NY= 5 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NY= 6 ( 8) | 8 7 6 5 4 3 2 1 0 * * NY= 7 ( 7) | 7 6 5 4 3 2 1 0 * * NY= 8 ( 6) | 6 5 4 3 2 1 0 * * NY= 9 ( 5) | 5 4 3 2 1 0 * * NY=10 ( 4) | 4 3 2 1 0 * * NY=11 ( 3) | 3 2 1 0 * * NY=12 ( 2) | 2 1 0 * * NY=13 ( 1) | 1 0 * * NY=14 ( 0) | 0 * * * * * * NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 * * MAX.NY => 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------------------- * * NY= 0 (14) | 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 1 (13) | 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 2 (12) | 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 3 (11) | 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 4 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NY= 5 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NY= 6 ( 8) | 8 7 6 5 4 3 2 1 0 * * NY= 7 ( 7) | 7 6 5 4 3 2 1 0 * * NY= 8 ( 6) | 6 5 4 3 2 1 0 * * NY= 9 ( 5) | 5 4 3 2 1 0 * * NY=10 ( 4) | 4 3 2 1 0 * * NY=11 ( 3) | 3 2 1 0 * * NY=12 ( 2) | 2 1 0 * * NY=13 ( 1) | 1 0 * * NY=14 ( 0) | 0 * * * * * * NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 * * MAX.NX => 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------------------- * * NX= 0 (14) | 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 1 (13) | 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 2 (12) | 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 3 (11) | 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 4 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NX= 5 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NX= 6 ( 8) | 8 7 6 5 4 3 2 1 0 * * NX= 7 ( 7) | 7 6 5 4 3 2 1 0 * * NX= 8 ( 6) | 6 5 4 3 2 1 0 * * NX= 9 ( 5) | 5 4 3 2 1 0 * * NX=10 ( 4) | 4 3 2 1 0 * * NX=11 ( 3) | 3 2 1 0 * * NX=12 ( 2) | 2 1 0 * * NX=13 ( 1) | 1 0 * * NX=14 ( 0) | 0 * * * ******************************************************************************* ******************************************************************************* * * * AUTOMATED HANDLING OF THE COULOMB FILE * * * ******************************************************************************* * * * THE COULOMB FILE READ FROM THE DISC, FILE NAME (FILCOU): * * ge064-a.cou * * * ******************************************************************************* * * * NUCLIDE: N = 32 Z = 32 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 LIPKIN-NOGAMI = 0.50 * * * ******************************************************************************* * * * BROYDEN ACCELERATION WILL START WHEN STABILITY IS LOWER THAN: 1000.000 MEV * * BROYDEN MIXING PARAMETER: ALPHAM = 0.80, ITERATIONS IN MEMORY: N_ITER = 7 * * * ******************************************************************************* * * * DIRECT COULOMB TERMS SOLVED BY GREEN FUNCTION AND FOURIER TRANSFORM METHOD * * COULOMB PARAMETERS: NUMCOU = 80 NUMETA = 79 BOUCOU = 20.0 FURMAX =0.25 * * * ******************************************************************************* * * * EXCHANGE COULOMB TERMS TREATED BY THE SLATER APPROXIMATION * * * ******************************************************************************* * * * MAXIMUM NUMBERS OF MULTIPOLE MOMENTS USED FOR THE CONSTRAINTS, NMUCON=2 * * FOR THE COULOMB FIELD, NMUCOU=4 * * FOR THE OUTPUT INFO, NMUPRI=6 * * * ******************************************************************************* * * * 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: NO * * PROTONS : NO * * MASS : YES * * * ******************************************************************************* * * * PRINTING THE MOMENTS, MOMENTA AND DEFORMATIONS IN THE INTRINSIC FRAME: NO * * * ******************************************************************************* * * * NILSSON LABELS OF SINGLE-PARTICLE STATES SHOWN WITH RESPECT TO THE: Z-AXIS * * * ******************************************************************************* * * * NUMBERS OF LEVELS IN TIME-REVERSAL/CHARGE BLOCKS: TIME-UP TIME-DOWN * * NEUTRONS: 32 32 * * PROTONS: 32 32 * * * ******************************************************************************* * * * CALCULATIONS WITHOUT PAIRING * * * ******************************************************************************* * * * CALCULATIONS WITH PARITY/SIGNATURE SYMMETRY * * * ******************************************************************************* * * * CALCULATIONS WITH TIME-ODD SYMMETRIES: T*SIMPLEX_X => CONSERVED * * T*SIMPLEX_Y => CONSERVED * * T*SIMPLEX_Z => CONSERVED * * * ******************************************************************************* * * * CALCULATIONS WITHOUT ROTATION * * * ******************************************************************************* ******************************************************************************* * * * PARITY/SIGNATURE CONFIGURATIONS: * * * * V A C U U M P A R T I C L E S H O L E S * * =========== ================= ========= * * (++) (+-) (-+) (--) (++) (+-) (-+) (--) (++) (+-) (-+) (--) * * * * NEUTRONS: 7 7 9 9 0 0 0 0 0 0 0 0 * * PROTONS : 7 7 9 9 0 0 0 0 0 0 0 0 * * * ******************************************************************************* ******************************************************************************* * * * CALLING SUBROUTINE "RECORD" TO READ THE REPLAY FILE OBTAINED PREVIOUSLY * * ge064-a.rec * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 3) -39.297 ( 0 ----- 2) | 1, 0, 1, 3/2> -100 0.141 0.296 2.108 * * 4) -36.855 ( 0 ----- 3) | 1, 0, 1, 1/2> -100 0.043 -0.363 -8.396 * * 5) -31.044 ( 2 ----- 0) | 2, 2, 0, 1/2> 100 -0.443 -0.156 0.353 * * 6) -29.981 ( 3 ----- 0) | 2, 1, 1, 3/2> 100 -0.743 -0.048 0.065 * * 7) -28.049 ( 4 ----- 0) | 2, 0, 2, 5/2> 100 -0.032 -0.078 2.409 * * 8) -26.501 ( 5 ----- 0) | 2, 1, 1, 1/2> 100 0.431 -0.311 -0.721 * * 9) -24.072 ( 6 ----- 0) | 2, 0, 2, 3/2> 100 -0.279 0.013 -0.046 * * 10) -21.241 ( 7 ----- 0) | 2, 0, 0, 1/2> 100 0.072 -0.419 -5.855 * * 11) -20.938 ( 0 ----- 4) | 3, 2, 1, 1/2> -100 -0.063 -0.005 0.080 * * 12) -19.741 ( 0 ----- 5) | 3, 2, 1, 3/2> -100 -1.534 -0.220 0.144 * * 13) -18.718 ( 0 ----- 6) | 3, 1, 2, 5/2> -100 -0.151 -0.057 0.375 * * 14) -16.538 ( 0 ----- 7) | 3, 0, 3, 7/2> -100 -0.005 0.011 -2.242 * * 15) -14.824 ( 0 ----- 8) | 3, 2, 1, 1/2> -100 0.403 -0.363 -0.900 * * 16) -13.690 ( 0 ----- 9) | 3, 1, 2, 3/2> -100 -0.389 0.063 -0.162 * * 17) -11.047 ( 0 ----- 10) | 3, 1, 0, 1/2> -100 0.224 -0.372 -1.662 * * 18) -10.773 ( 0 ----- 11) | 3, 0, 1, 3/2> -100 0.451 0.278 0.616 * * 19) -10.600 ( 8 ----- 0) | 4, 3, 1, 1/2> 100 0.047 0.011 0.242 * * 20) -9.247 ( 9 ----- 0) | 4, 3, 1, 3/2> 100 -1.107 -0.137 0.124 * * 21) -9.031 ( 0 ----- 12) | 3, 0, 3, 5/2> -100 -0.314 -0.005 0.015 * * 22) -8.543 (10 ----- 0) | 4, 2, 2, 5/2> 100 -1.119 -0.148 0.132 * * 23) -7.190 (11 ----- 0) | 4, 1, 3, 7/2> 100 -0.040 0.005 -0.120 * * 24) -6.688 ( 0 ----- 13) | 3, 0, 1, 1/2> -100 0.375 -0.330 -0.879 * * 25) -4.725 (12 ----- 0) | 4, 0, 4, 9/2> 100 0.000 -0.003 7.391 * * 26) -3.927 (13 ----- 0) | 4, 4, 0, 1/2> 100 -0.228 -0.367 1.605 * * 27) -2.899 (14 ----- 0) | 4, 3, 1, 1/2> 100 0.185 0.246 1.325 * * 28) -1.198 (15 ----- 0) | 4, 1, 1, 3/2> 100 -0.128 -0.041 0.324 * * 29) -0.610 (16 ----- 0) | 4, 2, 2, 3/2> 100 0.251 -0.368 -1.465 * * 30) -0.060 ( 0 ----- 14) | 5, 3, 0, 1/2> -100 0.020 0.003 0.127 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 1) -38.415 ( 1 ----- 0) | 0, 0, 0, 1/2> 100 -0.498 -0.500 1.004 * * 2) -30.201 ( 0 ----- 1) | 1, 1, 0, 1/2> -100 -0.676 -0.437 0.646 * * 3) -28.756 ( 0 ----- 2) | 1, 0, 1, 3/2> -100 0.155 0.306 1.974 * * 4) -26.310 ( 0 ----- 3) | 1, 0, 1, 1/2> -100 0.027 -0.369 -1.E+01 * * 5) -20.781 ( 2 ----- 0) | 2, 2, 0, 1/2> 100 -0.408 -0.143 0.352 * * 6) -19.719 ( 3 ----- 0) | 2, 1, 1, 3/2> 100 -0.774 -0.056 0.072 * * 7) -17.752 ( 4 ----- 0) | 2, 0, 2, 5/2> 100 -0.043 -0.077 1.790 * * 8) -16.260 ( 5 ----- 0) | 2, 1, 1, 1/2> 100 0.448 -0.311 -0.695 * * 9) -13.803 ( 6 ----- 0) | 2, 0, 2, 3/2> 100 -0.257 0.015 -0.058 * * 10) -10.963 ( 7 ----- 0) | 2, 0, 0, 1/2> 100 0.038 -0.428 -1.E+01 * * 11) -10.946 ( 0 ----- 4) | 3, 2, 1, 1/2> -100 -0.043 0.000 0.004 * * 12) -9.739 ( 0 ----- 5) | 3, 2, 1, 3/2> -100 -1.514 -0.218 0.144 * * 13) -8.731 ( 0 ----- 6) | 3, 1, 2, 5/2> -100 -0.199 -0.062 0.311 * * 14) -6.534 ( 0 ----- 7) | 3, 0, 3, 7/2> -100 -0.005 0.011 -2.362 * * 15) -4.875 ( 0 ----- 8) | 3, 2, 1, 1/2> -100 0.441 -0.352 -0.798 * * 16) -3.782 ( 0 ----- 9) | 3, 1, 2, 3/2> -100 -0.410 0.050 -0.121 * * 17) -1.201 ( 0 ----- 10) | 3, 1, 2, 3/2> -100 0.254 -0.380 -1.495 * * 18) -0.923 ( 0 ----- 11) | 3, 0, 1, 3/2> -100 0.436 0.289 0.661 * * 19) -0.911 ( 8 ----- 0) | 4, 3, 1, 1/2> 100 0.051 0.012 0.238 * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 6.957427 3.479801 6.957427 0.002174 * * DRHOD = 0.890549 0.445413 0.890549 0.000278 * * DLPR_ = -2.217743 -1.109782 -2.217743 -0.001821 * * DTAU_ = 5.677858 2.840222 5.677858 0.002586 * * DSCU_ = 0.061804 0.030928 0.061804 0.000052 * * DDIV_ = 0.386887 0.193591 0.386887 0.000295 * * * * DSPI_ = 0.000000 0.000000 0.000000 0.000000 * * DSPID = 0.000000 0.000000 0.000000 0.000000 * * DLPS_ = 0.000000 0.000000 0.000000 0.000000 * * DCUR_ = 0.000000 0.000000 0.000000 0.000000 * * DKIS_ = 0.000000 0.000000 0.000000 0.000000 * * DROT_ = 0.000000 0.000000 0.000000 0.000000 * * * ******************************************************************************* ******************************************************************************* ******************************************************************************* * * * CONTRIBUTIONS TO ENERGY IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * ERHO_ = -4810.082655 1865.716908 -2944.948565 0.582817 * * ERHOD = 1558.459889 -779.473459 779.229945 -0.243515 * * ELPR_ = 177.419455 -37.801949 139.648516 -0.031011 * * ETAU_ = 425.839313 -173.963591 251.954927 -0.079205 * * ESCU_ = 0.000000 0.000000 0.000000 0.000000 * * EDIV_ = -23.213209 -11.615469 -34.819813 -0.008865 * * ============ ============ ============ ============ * * SUM EVEN: -2671.577207 862.862439 -1808.934990 0.220222 * * * * ESPI_ = 0.000000 0.000000 0.000000 0.000000 * * ESPID = 0.000000 0.000000 0.000000 0.000000 * * ELPS_ = 0.000000 0.000000 0.000000 0.000000 * * ECUR_ = 0.000000 0.000000 0.000000 0.000000 * * EKIS_ = 0.000000 0.000000 0.000000 0.000000 * * EROT_ = 0.000000 0.000000 0.000000 0.000000 * * ============ ============ ============ ============ * * SUM ODD: 0.000000 0.000000 0.000000 0.000000 * * * ******************************************************************************* ******************************************************************************* * * * EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL * * * * ALPHA = 90.00000 BETA = 0.00000 GAMMA = 0.00000 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 3.9238 X_RMS = 1.9629 Y_RMS = 2.2347 Z_RMS = 2.5592 * * * * R_GEO = 5.0656 X_GEO = 4.3891 Y_GEO = 4.9969 Z_GEO = 5.7226 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * ******************************************************************************* * * * Q00 = 64.0000 ............. ............. ............. ............. * * * * Q10 = ZERO Q11 = ZERO ............. ............. ............. * * * * Q20 = 2.7218 Q21 = ZERO Q22 = -1.2647 ............. ............. * * * * Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. * * * * Q40 = -0.0047 Q41 = ZERO Q42 = -0.0407 Q43 = ZERO Q44 = 0.0017 * * * * Q50 = ZERO Q51 = ZERO Q52 = ZERO Q53 = ZERO Q54 = ZERO * * * * Q55 = ZERO * * * * Q60 = -0.0033 Q61 = ZERO Q62 = 3.9E-04 Q63 = ZERO Q64 =-4.6E-04 * * * * Q65 = ZERO Q66 = 2.2E-05 * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (EXACT MULTIPOLE MOMENTS) TOTAL * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = 0.2173 B21 = ZERO B22 = -0.1189 ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 = -0.0275 B41 = ZERO B42 = -0.0357 B43 = ZERO B44 = -0.0046 * * * * B50 = ZERO B51 = ZERO B52 = ZERO B53 = ZERO B54 = ZERO * * * * B55 = ZERO * * * * B60 = -0.0068 B61 = ZERO B62 = 0.0156 B63 = ZERO B64 = -0.0105 * * * * B65 = ZERO B66 =-3.0E-04 * * * ******************************************************************************* * * * BOHR DENSITY = 0.1220 FM^(-3) BOHR RADIUS = 4.9784 FM * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) TOTAL * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = 0.2189 B21 = ZERO B22 = -0.1017 ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 = -0.0047 B41 = ZERO B42 = -0.0572 B43 = ZERO B44 = 0.0024 * * * * B50 = ZERO B51 = ZERO B52 = ZERO B53 = ZERO B54 = ZERO * * * * B55 = ZERO * * * * B60 = -0.0126 B61 = ZERO B62 = 0.0022 B63 = ZERO B64 = -0.0025 * * * * B65 = ZERO B66 = 1.2E-04 * * * ******************************************************************************* * * * BOHR DENSITY = 0.1175 FM^(-3) BOHR RADIUS = 5.0656 FM * * * ******************************************************************************* ******************************************************************************* * * * NEUTRON CONFIGURATIONS * * ======================= * * P S 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 * * --- -------------------------------------------------------------- * * * * CONF: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * PROTON CONFIGURATIONS * * ======================= * * P S 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 * * --- -------------------------------------------------------------- * * * * CONF: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * ENERGIES (MEV) * * * ******************************************************************************* * * * KINETIC: (NEU)= 560.318123 (PRO)= 544.009289 (TOT)= 1104.327412 * * SUM EPS: (NEU)= -862.880886 (PRO)= -535.130276 (TOT)= -1398.011162 * * 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)= -12.368916 (PRO)= -2.491586 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 176.328216 (EXC)= -12.350372 (TOT)= 163.977844 * * * * CONSTR. (MULT)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 * * 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= 389.493215 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -540.409513 * * * * SPIN-ORB (EVE)= -34.828678 (ODD)= 0.000000 (TOT)= -34.828678 * * SKYRME: (EVE)= -1808.714769 (ODD)= 0.000000 (TOT)= -1808.714769 * * * * TOTAL: (STAB)= -0.042368 (SP)= -540.451881 (FUN)= -540.409513 * * * ******************************************************************************* ******************************************************************************* * * * ITERATION RESTARTED FROM THE DISC, THE FOLLOWING FILE(S) READ: * * ge064-a.rec * * * ******************************************************************************* * * * CONVERGENCE REPORT * * * ******************************************************************************* * * * ITER ENERGY STABILITY Q_2 GAMMA SPIN OMEGA ANGLE EPAIR * * * * 21 -540.409513 -0.042368 3.001 24.92 0.00 0.000 0.000 0.00 * * 22 -540.411664 -0.028724 3.008 24.98 0.00 0.000 0.000 0.00 * * 23 -540.417976 0.042960 3.043 25.29 0.00 0.000 0.000 0.00 * * 24 -540.418453 0.017936 3.051 25.38 0.00 0.000 0.000 0.00 * * 25 -540.418547 0.005002 3.051 25.41 0.00 0.000 0.000 0.00 * * 26 -540.418684 -0.005928 3.051 25.48 0.00 0.000 0.000 0.00 * * 27 -540.418746 -0.009570 3.052 25.53 0.00 0.000 0.000 0.00 * * 28 -540.418767 -0.007681 3.052 25.58 0.00 0.000 0.000 0.00 * * 29 -540.418768 -0.003235 3.052 25.58 0.00 0.000 0.000 0.00 * * 30 -540.418768 -0.000712 3.052 25.58 0.00 0.000 0.000 0.00 * * 31 -540.418768 0.000073 3.052 25.58 0.00 0.000 0.000 0.00 * * 32 -540.418768 0.000332 3.052 25.58 0.00 0.000 0.000 0.00 * * 33 -540.418768 0.000357 3.052 25.58 0.00 0.000 0.000 0.00 * * 34 -540.418768 0.000107 3.052 25.58 0.00 0.000 0.000 0.00 * * 35 -540.418768 -0.000013 3.052 25.58 0.00 0.000 0.000 0.00 * * 36 -540.418768 -0.000031 3.052 25.58 0.00 0.000 0.000 0.00 * * 37 -540.418768 -0.000013 3.052 25.58 0.00 0.000 0.000 0.00 * * 38 -540.418768 -3.49E-06 3.052 25.58 0.00 0.000 0.000 0.00 * * 39 -540.418768 2.61E-07 3.052 25.58 0.00 0.000 0.000 0.00 * * 40 -540.418768 1.32E-06 3.052 25.58 0.00 0.000 0.000 0.00 * * 41 -540.418768 4.39E-07 3.052 25.58 0.00 0.000 0.000 0.00 * * 42 -540.418768 6.85E-08 3.052 25.58 0.00 0.000 0.000 0.00 * * 43 -540.418768 1.74E-08 3.052 25.58 0.00 0.000 0.000 0.00 * * 44 -540.418768 1.29E-08 3.052 25.58 0.00 0.000 0.000 0.00 * * 45 -540.418768 1.55E-08 3.052 25.58 0.00 0.000 0.000 0.00 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 3) -39.301 ( 0 ----- 2) | 1, 0, 1, 3/2> -100 0.155 0.309 1.994 * * 4) -36.748 ( 0 ----- 3) | 1, 0, 1, 1/2> -100 0.021 -0.375 -2.E+01 * * 5) -31.088 ( 2 ----- 0) | 2, 2, 0, 1/2> 100 -0.412 -0.147 0.357 * * 6) -30.001 ( 3 ----- 0) | 2, 1, 1, 3/2> 100 -0.756 -0.050 0.067 * * 7) -28.035 ( 4 ----- 0) | 2, 0, 2, 5/2> 100 -0.043 -0.090 2.085 * * 8) -26.499 ( 5 ----- 0) | 2, 1, 1, 1/2> 100 0.417 -0.318 -0.762 * * 9) -24.064 ( 6 ----- 0) | 2, 0, 2, 3/2> 100 -0.268 0.026 -0.096 * * 10) -21.072 ( 7 ----- 0) | 2, 0, 0, 1/2> 100 0.068 -0.420 -6.202 * * 11) -21.004 ( 0 ----- 4) | 3, 2, 1, 1/2> -100 -0.041 0.000 -0.001 * * 12) -19.759 ( 0 ----- 5) | 3, 2, 1, 3/2> -100 -1.525 -0.221 0.145 * * 13) -18.738 ( 0 ----- 6) | 3, 1, 2, 5/2> -100 -0.173 -0.063 0.366 * * 14) -16.515 ( 0 ----- 7) | 3, 0, 3, 7/2> -100 -0.004 0.014 -3.212 * * 15) -14.860 ( 0 ----- 8) | 3, 2, 1, 1/2> -100 0.398 -0.364 -0.913 * * 16) -13.744 ( 0 ----- 9) | 3, 1, 2, 3/2> -100 -0.387 0.068 -0.175 * * 17) -10.986 ( 0 ----- 10) | 3, 1, 0, 1/2> -100 0.217 -0.366 -1.683 * * 18) -10.806 ( 0 ----- 11) | 3, 0, 1, 1/2> -100 0.474 0.258 0.544 * * 19) -10.681 ( 8 ----- 0) | 4, 3, 1, 1/2> 100 0.051 0.012 0.239 * * 20) -9.278 ( 9 ----- 0) | 4, 3, 1, 3/2> 100 -1.010 -0.126 0.124 * * 21) -8.941 ( 0 ----- 12) | 3, 0, 3, 5/2> -100 -0.317 0.012 -0.038 * * 22) -8.566 (10 ----- 0) | 4, 2, 2, 5/2> 100 -1.204 -0.159 0.132 * * 23) -7.214 (11 ----- 0) | 4, 1, 3, 7/2> 100 -0.046 0.006 -0.121 * * 24) -6.485 ( 0 ----- 13) | 3, 0, 1, 1/2> -100 0.357 -0.339 -0.948 * * 25) -4.697 (12 ----- 0) | 4, 0, 4, 9/2> 100 0.000 -0.005 2.E+01 * * 26) -3.974 (13 ----- 0) | 4, 4, 0, 1/2> 100 -0.221 -0.362 1.636 * * 27) -2.960 (14 ----- 0) | 4, 3, 1, 1/2> 100 0.177 0.241 1.362 * * 28) -1.255 (15 ----- 0) | 4, 1, 1, 3/2> 100 -0.102 -0.025 0.242 * * 29) -0.613 (16 ----- 0) | 4, 2, 2, 3/2> 100 0.252 -0.386 -1.528 * * 30) -0.152 ( 0 ----- 14) | 5, 3, 0, 1/2> -100 0.018 0.002 0.124 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 1) -38.394 ( 1 ----- 0) | 0, 0, 0, 1/2> 100 -0.498 -0.500 1.005 * * 2) -30.227 ( 0 ----- 1) | 1, 1, 0, 1/2> -100 -0.668 -0.437 0.653 * * 3) -28.766 ( 0 ----- 2) | 1, 0, 1, 3/2> -100 0.169 0.317 1.882 * * 4) -26.207 ( 0 ----- 3) | 1, 0, 1, 1/2> -100 0.006 -0.381 -7.E+01 * * 5) -20.828 ( 2 ----- 0) | 2, 2, 0, 1/2> 100 -0.377 -0.134 0.356 * * 6) -19.744 ( 3 ----- 0) | 2, 1, 1, 3/2> 100 -0.786 -0.058 0.074 * * 7) -17.741 ( 4 ----- 0) | 2, 0, 2, 5/2> 100 -0.055 -0.089 1.618 * * 8) -16.262 ( 5 ----- 0) | 2, 1, 1, 1/2> 100 0.434 -0.318 -0.731 * * 9) -13.799 ( 6 ----- 0) | 2, 0, 2, 3/2> 100 -0.247 0.027 -0.110 * * 10) -11.014 ( 0 ----- 4) | 3, 2, 1, 1/2> -100 -0.023 0.005 -0.196 * * 11) -10.802 ( 7 ----- 0) | 2, 0, 0, 1/2> 100 0.036 -0.428 -1.E+01 * * 12) -9.761 ( 0 ----- 5) | 3, 2, 1, 3/2> -100 -1.499 -0.218 0.145 * * 13) -8.753 ( 0 ----- 6) | 3, 1, 2, 5/2> -100 -0.226 -0.069 0.306 * * 14) -6.514 ( 0 ----- 7) | 3, 0, 3, 7/2> -100 -0.004 0.014 -3.594 * * 15) -4.914 ( 0 ----- 8) | 3, 2, 1, 1/2> -100 0.435 -0.353 -0.811 * * 16) -3.838 ( 0 ----- 9) | 3, 1, 2, 3/2> -100 -0.407 0.054 -0.132 * * 17) -1.147 ( 0 ----- 10) | 3, 1, 0, 1/2> -100 0.256 -0.380 -1.482 * * 18) -0.994 ( 8 ----- 0) | 4, 3, 1, 1/2> 100 0.054 0.013 0.236 * * 19) -0.960 ( 0 ----- 11) | 3, 0, 1, 1/2> -100 0.449 0.275 0.612 * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 6.955028 3.478600 6.955028 0.002173 * * DRHOD = 0.890012 0.445145 0.890012 0.000278 * * DLPR_ = -2.217661 -1.109740 -2.217661 -0.001820 * * DTAU_ = 5.673465 2.838024 5.673465 0.002583 * * DSCU_ = 0.060947 0.030499 0.060947 0.000051 * * DDIV_ = 0.384268 0.192280 0.384268 0.000292 * * * * DSPI_ = 0.000000 0.000000 0.000000 0.000000 * * DSPID = 0.000000 0.000000 0.000000 0.000000 * * DLPS_ = 0.000000 0.000000 0.000000 0.000000 * * DCUR_ = 0.000000 0.000000 0.000000 0.000000 * * DKIS_ = 0.000000 0.000000 0.000000 0.000000 * * DROT_ = 0.000000 0.000000 0.000000 0.000000 * * * ******************************************************************************* ******************************************************************************* ******************************************************************************* * * * CONTRIBUTIONS TO ENERGY IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * ERHO_ = -4808.423767 1865.073257 -2943.932919 0.582409 * * ERHOD = 1557.521208 -779.003872 778.760604 -0.243268 * * ELPR_ = 177.412894 -37.800530 139.643352 -0.030988 * * ETAU_ = 425.509865 -173.828966 251.760003 -0.079105 * * ESCU_ = 0.000000 0.000000 0.000000 0.000000 * * EDIV_ = -23.056052 -11.536777 -34.584077 -0.008752 * * ============ ============ ============ ============ * * SUM EVEN: -2671.035851 862.903111 -1808.353036 0.220296 * * * * ESPI_ = 0.000000 0.000000 0.000000 0.000000 * * ESPID = 0.000000 0.000000 0.000000 0.000000 * * ELPS_ = 0.000000 0.000000 0.000000 0.000000 * * ECUR_ = 0.000000 0.000000 0.000000 0.000000 * * EKIS_ = 0.000000 0.000000 0.000000 0.000000 * * EROT_ = 0.000000 0.000000 0.000000 0.000000 * * ============ ============ ============ ============ * * SUM ODD: 0.000000 0.000000 0.000000 0.000000 * * * ******************************************************************************* ******************************************************************************* * * * EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL * * * * ALPHA =270.00000 BETA = 0.00000 GAMMA = 0.00000 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 3.9260 X_RMS = 1.9562 Y_RMS = 2.2395 Z_RMS = 2.5635 * * * * R_GEO = 5.0685 X_GEO = 4.3742 Y_GEO = 5.0077 Z_GEO = 5.7322 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * ******************************************************************************* * * * Q00 = 64.0000 ............. ............. ............. ............. * * * * Q10 = ZERO Q11 = ZERO ............. ............. ............. * * * * Q20 = 2.7527 Q21 = ZERO Q22 = -1.3178 ............. ............. * * * * Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. * * * * Q40 = -0.0058 Q41 = ZERO Q42 = -0.0405 Q43 = ZERO Q44 = 0.0026 * * * * Q50 = ZERO Q51 = ZERO Q52 = ZERO Q53 = ZERO Q54 = ZERO * * * * Q55 = ZERO * * * * Q60 = -0.0035 Q61 = ZERO Q62 = 5.5E-04 Q63 = ZERO Q64 =-5.3E-04 * * * * Q65 = ZERO Q66 =-2.5E-05 * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (EXACT MULTIPOLE MOMENTS) TOTAL * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = 0.2203 B21 = ZERO B22 = -0.1245 ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 = -0.0290 B41 = ZERO B42 = -0.0341 B43 = ZERO B44 = -0.0037 * * * * B50 = ZERO B51 = ZERO B52 = ZERO B53 = ZERO B54 = ZERO * * * * B55 = ZERO * * * * B60 = -0.0072 B61 = ZERO B62 = 0.0161 B63 = ZERO B64 = -0.0112 * * * * B65 = ZERO B66 =-4.9E-04 * * * ******************************************************************************* * * * BOHR DENSITY = 0.1219 FM^(-3) BOHR RADIUS = 4.9779 FM * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) TOTAL * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = 0.2212 B21 = ZERO B22 = -0.1059 ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 = -0.0058 B41 = ZERO B42 = -0.0568 B43 = ZERO B44 = 0.0036 * * * * B50 = ZERO B51 = ZERO B52 = ZERO B53 = ZERO B54 = ZERO * * * * B55 = ZERO * * * * B60 = -0.0135 B61 = ZERO B62 = 0.0030 B63 = ZERO B64 = -0.0029 * * * * B65 = ZERO B66 =-1.4E-04 * * * ******************************************************************************* * * * BOHR DENSITY = 0.1173 FM^(-3) BOHR RADIUS = 5.0685 FM * * * ******************************************************************************* ******************************************************************************* * * * NEUTRON CONFIGURATIONS * * ======================= * * P S 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 * * --- -------------------------------------------------------------- * * * * CONF: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * PROTON CONFIGURATIONS * * ======================= * * P S 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 * * --- -------------------------------------------------------------- * * * * CONF: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + + 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 2 2 2 2 2 2 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * ENERGIES (MEV) * * * ******************************************************************************* * * * KINETIC: (NEU)= 560.039045 (PRO)= 543.749611 (TOT)= 1103.788656 * * SUM EPS: (NEU)= -862.750333 (PRO)= -535.125986 (TOT)= -1397.876319 * * 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)= -12.364786 (PRO)= -2.492330 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 176.274120 (EXC)= -12.348805 (TOT)= 163.925316 * * * * CONSTR. (MULT)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 * * 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= 389.258668 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -540.418768 * * * * SPIN-ORB (EVE)= -34.592829 (ODD)= 0.000000 (TOT)= -34.592829 * * SKYRME: (EVE)= -1808.132740 (ODD)= 0.000000 (TOT)= -1808.132740 * * * * TOTAL: (STAB)= -4.08E-08 (SP)= -540.418768 (FUN)= -540.418768 * * * ******************************************************************************* ******************************************************************************* * * * CALLING SUBROUTINE "RECORD" TO SAVE THE RECORD FILE CONTAINING RESULTS * * ge064-a.rec * * * ******************************************************************************* ******************************************************************************* * * * CALCULATION OF THE DATA SET NO. 2 HAS BEEN SUCCESSFULLY COMPLETED * * * ******************************************************************************* ******************************************************************************* * * * CLASSICAL NUCLEAR SURFACE DEFINED FOR: N = 32 Z = 32 * * * ******************************************************************************* * * * SPHERICAL SHAPE REQUESTED FOR THE SURFACE * * * ******************************************************************************* * * * HOMEGA= 12.3000 FCHOM0= 1.2000 * * * * OSCILLATOR FREQUENCIES: HBAROX= 12.3000 HBAROY= 12.3000 HBAROZ= 12.3000 * * * * MOMENTS OF INERTIA: XMOMFC= 14.9370 YMOMFC= 14.9370 ZMOMFC= 14.9370 * * * * CENTERS OF MASS: CMSXFC= 0.0000 CMSYFC= 0.0000 CMSZFC= 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * SKYRME FORCE DEFINITION * * * ******************************************************************************* * * * PARAMETER SET SIII: T0= -1128.75 T1= 395.00 T2= -95.00 T3= 14000.00 * * * * POWER=1/1 W=120.000 X0= 0.45000 X1= 0.00000 X2= 0.00000 X3= 1.00000 * * * * ETA_J=0 ETA_W=0 ETACM=0 ETA_M=1 HBM= 20.73000 (SKYRME-FORCE SPIN-ORBIT) * * * ******************************************************************************* ******************************************************************************* * * * COUPLING CONSTANTS DEFINING THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * CRHO_ = -691.359375 536.156250 -423.281250 268.078125 * * CRHOD = 1750.000000 -1750.000000 875.000000 -875.000000 * * CLPR_ = -80.000000 34.062500 -62.968750 17.031250 * * CTAU_ = 75.000000 -61.250000 44.375000 -30.625000 * * CSCU_ = 0.000000 0.000000 0.000000 0.000000 * * CDIV_ = -60.000000 -60.000000 -90.000000 -30.000000 * * * * CSPI_ = -126.984375 282.187500 14.109375 141.093750 * * CSPID = 583.333333 -583.333333 291.666667 -291.666667 * * CLPS_ = 0.000000 34.062500 17.031250 17.031250 * * CCUR_ = -75.000000 61.250000 -44.375000 30.625000 * * CKIS_ = 0.000000 0.000000 0.000000 0.000000 * * CROT_ = -60.000000 -60.000000 -90.000000 -30.000000 * * * ******************************************************************************* ******************************************************************************* * * * PHYSICAL CONSTANTS: H_BARC=197.32891000 HBCOE2=137.03602000 * * * * XMASSN=938.90590000 XMASSP=938.27231000 * * * * HBMASS= 20.73000000 HBMRPA= 20.40609375 * * * * ECHAR2= 1.43997841 COULEX= -1.06350868 * * * ******************************************************************************* * * * OSCILLATOR LENGTHS: X= 1.8362302 Y= 1.8362302 Z= 1.8362302 * * * * OSCILLATOR CONSTANTS: X= 0.5445940 Y= 0.5445940 Z= 0.5445940 * * * * OSCILLATOR FREQUENCIES: X=12.3000000 Y=12.3000000 Z=12.3000000 * * * ******************************************************************************* * * * BASIS CUT-OFF CONTROL PARAMETERS: NXMAXX = 14 NYMAXX = 14 NZMAXX = 14 * * * * OPTIMUM NUMBERS OF GAUSS POINTS: NXHERM = 30 NYHERM = 30 NZHERM = 30 * * * * NLIMIT= 680 LDBASE= 680 MCOUNT= 3375 * * * * ENECUT= 800.0000 ELIMIT= 190.6500 * * * ******************************************************************************* ******************************************************************************* * * * SHAPE OF THE OSCILLATOR-BASIS DIAMOND (THREE CONSECUTIVE PROJECTIONS) * * * ******************************************************************************* * * * NX ===>>> 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 * * MAX.NY => 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------------------- * * NY= 0 (14) | 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 1 (13) | 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 2 (12) | 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 3 (11) | 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 4 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NY= 5 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NY= 6 ( 8) | 8 7 6 5 4 3 2 1 0 * * NY= 7 ( 7) | 7 6 5 4 3 2 1 0 * * NY= 8 ( 6) | 6 5 4 3 2 1 0 * * NY= 9 ( 5) | 5 4 3 2 1 0 * * NY=10 ( 4) | 4 3 2 1 0 * * NY=11 ( 3) | 3 2 1 0 * * NY=12 ( 2) | 2 1 0 * * NY=13 ( 1) | 1 0 * * NY=14 ( 0) | 0 * * * * * * NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 * * MAX.NY => 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------------------- * * NY= 0 (14) | 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 1 (13) | 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 2 (12) | 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 3 (11) | 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 4 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NY= 5 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NY= 6 ( 8) | 8 7 6 5 4 3 2 1 0 * * NY= 7 ( 7) | 7 6 5 4 3 2 1 0 * * NY= 8 ( 6) | 6 5 4 3 2 1 0 * * NY= 9 ( 5) | 5 4 3 2 1 0 * * NY=10 ( 4) | 4 3 2 1 0 * * NY=11 ( 3) | 3 2 1 0 * * NY=12 ( 2) | 2 1 0 * * NY=13 ( 1) | 1 0 * * NY=14 ( 0) | 0 * * * * * * NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 * * MAX.NX => 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------------------- * * NX= 0 (14) | 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 1 (13) | 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 2 (12) | 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 3 (11) | 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 4 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NX= 5 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NX= 6 ( 8) | 8 7 6 5 4 3 2 1 0 * * NX= 7 ( 7) | 7 6 5 4 3 2 1 0 * * NX= 8 ( 6) | 6 5 4 3 2 1 0 * * NX= 9 ( 5) | 5 4 3 2 1 0 * * NX=10 ( 4) | 4 3 2 1 0 * * NX=11 ( 3) | 3 2 1 0 * * NX=12 ( 2) | 2 1 0 * * NX=13 ( 1) | 1 0 * * NX=14 ( 0) | 0 * * * ******************************************************************************* ******************************************************************************* * * * AUTOMATED HANDLING OF THE COULOMB FILE * * * ******************************************************************************* * * * THE COULOMB FILE READ FROM THE DISC, FILE NAME (FILCOU): * * ge064-a.cou * * * ******************************************************************************* * * * NUCLIDE: N = 32 Z = 32 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 LIPKIN-NOGAMI = 0.50 * * * ******************************************************************************* * * * BROYDEN ACCELERATION WILL START WHEN STABILITY IS LOWER THAN: 1000.000 MEV * * BROYDEN MIXING PARAMETER: ALPHAM = 0.80, ITERATIONS IN MEMORY: N_ITER = 7 * * * ******************************************************************************* * * * DIRECT COULOMB TERMS SOLVED BY GREEN FUNCTION AND FOURIER TRANSFORM METHOD * * COULOMB PARAMETERS: NUMCOU = 80 NUMETA = 79 BOUCOU = 20.0 FURMAX =0.25 * * * ******************************************************************************* * * * EXCHANGE COULOMB TERMS TREATED BY THE SLATER APPROXIMATION * * * ******************************************************************************* * * * MAXIMUM NUMBERS OF MULTIPOLE MOMENTS USED FOR THE CONSTRAINTS, NMUCON=2 * * FOR THE COULOMB FIELD, NMUCOU=4 * * FOR THE OUTPUT INFO, NMUPRI=6 * * * ******************************************************************************* * * * 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: NO * * PROTONS : NO * * MASS : YES * * * ******************************************************************************* * * * PRINTING THE MOMENTS, MOMENTA AND DEFORMATIONS IN THE INTRINSIC FRAME: NO * * * ******************************************************************************* * * * NILSSON LABELS OF SINGLE-PARTICLE STATES SHOWN WITH RESPECT TO THE: Z-AXIS * * * ******************************************************************************* * * * NUMBERS OF LEVELS IN TIME-REVERSAL/CHARGE BLOCKS: TIME-UP TIME-DOWN * * NEUTRONS: 32 32 * * PROTONS: 32 32 * * * ******************************************************************************* * * * 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.575000 (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: 7 7 9 9 0 0 0 0 0 0 0 0 * * PROTONS : 7 7 9 9 0 0 0 0 0 0 0 0 * * * ******************************************************************************* ******************************************************************************* * * * CALLING SUBROUTINE "RECORD" TO READ THE REPLAY FILE OBTAINED PREVIOUSLY * * ge064-a.rec * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 5) -39.316 ( 0, 0, 0, 2) | 1, 0, 1, 3/2> -100 -0.085 0.191 -2.250 * * 6) -39.154 ( 0, 0, 2, 0) | 1, 0, 1, 3/2> -100 -0.334 -0.327 0.980 * * 7) -36.680 ( 0, 0, 0, 3) | 1, 0, 1, 1/2> -100 -0.277 -0.416 1.503 * * 8) -36.678 ( 0, 0, 3, 0) | 1, 0, 1, 1/2> -100 -0.206 0.341 -1.656 * * 9) -31.798 ( 0, 2, 0, 0) | 2, 2, 0, 1/2> 100 1.798 0.452 0.251 * * 10) -31.260 ( 2, 0, 0, 0) | 2, 1, 1, 3/2> 100 0.767 0.173 0.226 * * 11) -30.185 ( 0, 3, 0, 0) | 2, 1, 1, 3/2> 100 0.099 -0.159 -1.614 * * 12) -29.425 ( 3, 0, 0, 0) | 2, 2, 0, 1/2> 100 -0.974 -0.275 0.282 * * 13) -27.921 ( 0, 4, 0, 0) | 2, 0, 2, 5/2> 100 -0.412 0.039 -0.096 * * 14) -27.873 ( 4, 0, 0, 0) | 2, 0, 2, 5/2> 100 -0.454 -0.126 0.277 * * 15) -26.720 ( 5, 0, 0, 0) | 2, 1, 1, 1/2> 100 0.248 -0.364 -1.463 * * 16) -26.294 ( 0, 5, 0, 0) | 2, 1, 1, 1/2> 100 -0.279 0.294 -1.054 * * 17) -24.195 ( 0, 6, 0, 0) | 2, 0, 2, 3/2> 100 0.180 -0.030 -0.167 * * 18) -23.891 ( 6, 0, 0, 0) | 2, 0, 2, 3/2> 100 -0.330 0.027 -0.082 * * 19) -21.827 ( 0, 0, 4, 0) | 3, 3, 0, 1/2> -100 2.547 0.429 0.168 * * 20) -21.499 ( 0, 0, 0, 4) | 3, 2, 1, 3/2> -100 1.486 0.267 0.180 * * 21) -21.022 ( 7, 0, 0, 0) | 2, 0, 0, 1/2> 100 -0.237 -0.436 1.843 * * 22) -20.944 ( 0, 7, 0, 0) | 2, 0, 0, 1/2> 100 -0.375 0.404 -1.078 * * 23) -20.510 ( 0, 0, 5, 0) | 3, 2, 1, 3/2> -100 0.967 0.071 0.074 * * 24) -19.389 ( 0, 0, 0, 5) | 3, 1, 2, 5/2> -100 -0.018 -0.089 5.044 * * 25) -18.469 ( 0, 0, 6, 0) | 3, 1, 2, 5/2> -100 -0.738 -0.124 0.168 * * 26) -17.949 ( 0, 0, 0, 6) | 3, 3, 0, 1/2> -100 -1.985 -0.352 0.177 * * 27) -16.298 ( 0, 0, 7, 0) | 3, 0, 3, 7/2> -100 -0.734 -0.073 0.099 * * 28) -16.282 ( 0, 0, 0, 7) | 3, 0, 3, 7/2> -100 -0.803 -0.062 0.077 * * 29) -15.249 ( 0, 0, 0, 8) | 3, 2, 1, 1/2> -100 0.884 -0.369 -0.417 * * 30) -14.814 ( 0, 0, 8, 0) | 3, 3, 0, 1/2> -100 0.230 0.373 1.622 * * 31) -14.027 ( 0, 0, 9, 0) | 3, 1, 2, 3/2> -100 0.583 -0.086 -0.148 * * 32) -13.589 ( 0, 0, 0, 9) | 3, 1, 2, 3/2> -100 -0.154 0.069 -0.450 * * 33) -11.599 ( 0, 8, 0, 0) | 4, 3, 1, 1/2> 100 3.199 0.411 0.128 * * 34) -11.414 ( 8, 0, 0, 0) | 4, 3, 1, 3/2> 100 2.206 0.299 0.136 * * 35) -11.168 ( 0, 0, 0,10) | 3, 1, 2, 3/2> -100 0.399 -0.158 -0.396 * * 36) -10.952 ( 0, 0, 0,11) | 3, 1, 0, 1/2> -100 0.067 0.012 0.177 * * 37) -10.845 ( 0, 0,10, 0) | 3, 0, 1, 3/2> -100 -0.305 0.217 -0.711 * * 38) -10.553 ( 0, 9, 0, 0) | 4, 3, 1, 3/2> 100 1.957 0.210 0.107 * * 39) -10.476 ( 0, 0,11, 0) | 3, 0, 1, 1/2> -100 -0.634 -0.139 0.219 * * 40) -9.544 ( 9, 0, 0, 0) | 4, 2, 2, 5/2> 100 0.838 0.044 0.053 * * 41) -9.035 ( 0, 0,12, 0) | 3, 0, 3, 5/2> -100 0.017 -0.017 -0.982 * * 42) -8.667 ( 0, 0, 0,12) | 3, 0, 3, 5/2> -100 -0.633 0.015 -0.024 * * 43) -8.501 ( 0,10, 0, 0) | 4, 1, 3, 7/2> 100 -0.200 -0.075 0.372 * * 44) -7.513 (10, 0, 0, 0) | 4, 1, 3, 7/2> 100 -1.134 -0.169 0.149 * * 45) -6.734 ( 0,11, 0, 0) | 4, 1, 3, 7/2> 100 -1.433 -0.188 0.131 * * 46) -6.524 ( 0, 0, 0,13) | 3, 0, 1, 1/2> -100 -0.220 -0.334 1.518 * * 47) -6.393 (11, 0, 0, 0) | 4, 1, 3, 7/2> 100 -2.689 -0.352 0.131 * * 48) -6.120 ( 0, 0,13, 0) | 3, 0, 1, 1/2> -100 -0.903 0.349 -0.387 * * 49) -4.466 ( 0,12, 0, 0) | 4, 0, 4, 9/2> 100 -0.385 0.056 -0.146 * * 50) -4.440 (12, 0, 0, 0) | 4, 0, 4, 9/2> 100 -0.890 -0.062 0.069 * * 51) -4.343 ( 0,13, 0, 0) | 4, 4, 0, 1/2> 100 0.615 0.313 0.508 * * 52) -4.124 (13, 0, 0, 0) | 4, 4, 0, 1/2> 100 0.720 -0.348 -0.483 * * 53) -3.276 (14, 0, 0, 0) | 4, 3, 1, 1/2> 100 0.865 0.242 0.280 * * 54) -3.075 ( 0,14, 0, 0) | 4, 2, 2, 3/2> 100 0.555 -0.264 -0.475 * * 55) -1.397 ( 0,15, 0, 0) | 4, 1, 1, 3/2> 100 0.342 -0.041 -0.120 * * 56) -1.328 (15, 0, 0, 0) | 4, 2, 0, 1/2> 100 0.373 -0.064 -0.172 * * 57) -1.211 ( 0, 0,14, 0) | 5, 3, 0, 1/2> -100 3.815 0.401 0.105 * * 58) -1.098 ( 0, 0, 0,14) | 5, 4, 1, 3/2> -100 2.953 0.322 0.109 * * 59) -0.809 (16, 0, 0, 0) | 4, 2, 2, 3/2> 100 0.454 -0.370 -0.817 * * 60) -0.532 ( 0,16, 0, 0) | 4, 1, 1, 1/2> 100 -0.006 0.398 -7.E+01 * * 61) -0.357 ( 0, 0,15, 0) | 5, 3, 2, 5/2> -100 2.961 0.287 0.097 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 1) -38.680 ( 0, 1, 0, 0) | 0, 0, 0, 1/2> 100 0.499 0.500 1.001 * * 2) -38.108 ( 1, 0, 0, 0) | 0, 0, 0, 1/2> 100 -0.496 -0.500 1.008 * * 3) -30.727 ( 0, 0, 1, 0) | 1, 1, 0, 1/2> -100 1.032 0.486 0.471 * * 4) -29.993 ( 0, 0, 0, 1) | 1, 1, 0, 1/2> -100 -0.140 -0.286 2.047 * * 5) -28.795 ( 0, 0, 0, 2) | 1, 0, 1, 3/2> -100 -0.054 0.206 -3.795 * * 6) -28.615 ( 0, 0, 2, 0) | 1, 0, 1, 3/2> -100 -0.336 -0.334 0.992 * * 7) -26.145 ( 0, 0, 3, 0) | 1, 0, 1, 1/2> -100 -0.194 0.348 -1.789 * * 8) -26.130 ( 0, 0, 0, 3) | 1, 0, 1, 1/2> -100 -0.293 -0.421 1.437 * * 9) -21.525 ( 0, 2, 0, 0) | 2, 2, 0, 1/2> 100 1.791 0.450 0.251 * * 10) -21.011 ( 2, 0, 0, 0) | 2, 1, 1, 3/2> 100 0.772 0.180 0.233 * * 11) -19.937 ( 0, 3, 0, 0) | 2, 1, 1, 3/2> 100 0.098 -0.165 -1.678 * * 12) -19.159 ( 3, 0, 0, 0) | 2, 2, 0, 1/2> 100 -0.973 -0.278 0.286 * * 13) -17.635 ( 0, 4, 0, 0) | 2, 0, 2, 5/2> 100 -0.394 0.041 -0.103 * * 14) -17.575 ( 4, 0, 0, 0) | 2, 0, 2, 5/2> 100 -0.449 -0.121 0.269 * * 15) -16.490 ( 5, 0, 0, 0) | 2, 1, 1, 1/2> 100 0.247 -0.367 -1.486 * * 16) -16.047 ( 0, 5, 0, 0) | 2, 1, 1, 1/2> 100 -0.298 0.291 -0.977 * * 17) -13.916 ( 0, 6, 0, 0) | 2, 0, 2, 3/2> 100 0.157 -0.028 -0.181 * * 18) -13.637 ( 6, 0, 0, 0) | 2, 0, 2, 3/2> 100 -0.314 0.029 -0.093 * * 19) -11.820 ( 0, 0, 4, 0) | 3, 3, 0, 1/2> -100 2.532 0.426 0.168 * * 20) -11.513 ( 0, 0, 0, 4) | 3, 2, 1, 3/2> -100 1.491 0.269 0.181 * * 21) -10.736 ( 7, 0, 0, 0) | 2, 0, 0, 1/2> 100 -0.262 -0.443 1.695 * * 22) -10.693 ( 0, 7, 0, 0) | 2, 0, 0, 1/2> 100 -0.342 0.411 -1.203 * * 23) -10.526 ( 0, 0, 5, 0) | 3, 2, 1, 3/2> -100 0.992 0.075 0.075 * * 24) -9.416 ( 0, 0, 0, 5) | 3, 1, 2, 5/2> -100 -0.006 -0.088 1.E+01 * * 25) -8.489 ( 0, 0, 6, 0) | 3, 1, 2, 5/2> -100 -0.747 -0.128 0.171 * * 26) -7.930 ( 0, 0, 0, 6) | 3, 3, 0, 1/2> -100 -2.006 -0.353 0.176 * * 27) -6.298 ( 0, 0, 7, 0) | 3, 0, 3, 7/2> -100 -0.731 -0.071 0.097 * * 28) -6.281 ( 0, 0, 0, 7) | 3, 0, 3, 7/2> -100 -0.805 -0.061 0.076 * * 29) -5.328 ( 0, 0, 0, 8) | 3, 2, 1, 1/2> -100 0.933 -0.365 -0.392 * * 30) -4.853 ( 0, 0, 8, 0) | 3, 3, 0, 1/2> -100 0.212 0.353 1.662 * * 31) -4.136 ( 0, 0, 9, 0) | 3, 1, 2, 3/2> -100 0.612 -0.064 -0.104 * * 32) -3.675 ( 0, 0, 0, 9) | 3, 1, 2, 3/2> -100 -0.163 0.063 -0.385 * * 33) -1.908 ( 0, 8, 0, 0) | 4, 3, 1, 1/2> 100 3.196 0.411 0.129 * * 34) -1.732 ( 8, 0, 0, 0) | 4, 3, 1, 3/2> 100 2.221 0.302 0.136 * * 35) -1.338 ( 0, 0, 0,10) | 3, 1, 2, 3/2> -100 0.390 -0.169 -0.432 * * 36) -1.109 ( 0, 0, 0,11) | 3, 1, 0, 1/2> -100 0.100 0.026 0.260 * * 37) -0.989 ( 0, 0,10, 0) | 3, 0, 1, 3/2> -100 -0.324 0.247 -0.764 * * 38) -0.881 ( 0, 9, 0, 0) | 4, 2, 2, 5/2> 100 1.976 0.213 0.108 * * 39) -0.642 ( 0, 0,11, 0) | 3, 0, 1, 1/2> -100 -0.619 -0.170 0.275 * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 6.949899 3.476037 6.949899 0.002175 * * DRHOD = 0.888969 0.444623 0.888969 0.000278 * * DLPR_ = -2.210555 -1.106184 -2.210555 -0.001813 * * DTAU_ = 5.674011 2.838302 5.674011 0.002593 * * DSCU_ = 0.058657 0.029353 0.058657 0.000049 * * DDIV_ = 0.380791 0.190540 0.380791 0.000289 * * * * DSPI_ = 0.002437 0.001220 0.002437 0.000003 * * DSPID = 0.000325 0.000163 0.000325 0.000000 * * DLPS_ = -0.006307 -0.003159 -0.006307 -0.000012 * * DCUR_ = 0.003582 0.001792 0.003582 0.000002 * * DKIS_ = 0.000548 0.000274 0.000548 0.000000 * * DROT_ = 0.000693 0.000347 0.000693 0.000000 * * * ******************************************************************************* ******************************************************************************* ******************************************************************************* * * * CONTRIBUTIONS TO ENERGY IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * ERHO_ = -4804.877709 1863.698856 -2941.761863 0.583010 * * ERHOD = 1555.694932 -778.091117 777.847466 -0.243651 * * ELPR_ = 176.844416 -37.679392 139.195898 -0.030873 * * ETAU_ = 425.550856 -173.846001 251.784257 -0.079401 * * ESCU_ = 0.000000 0.000000 0.000000 0.000000 * * EDIV_ = -22.847443 -11.432386 -34.271164 -0.008664 * * ============ ============ ============ ============ * * SUM EVEN: -2669.634947 862.649961 -1807.205406 0.220420 * * * * ESPI_ = -0.309437 0.344275 0.034382 0.000456 * * ESPID = 0.189505 -0.094869 0.094753 -0.000116 * * ELPS_ = 0.000000 -0.107614 -0.107412 -0.000202 * * ECUR_ = -0.268674 0.109774 -0.158965 0.000065 * * EKIS_ = 0.000000 0.000000 0.000000 0.000000 * * EROT_ = -0.041593 -0.020808 -0.062390 -0.000011 * * ============ ============ ============ ============ * * SUM ODD: -0.430199 0.230758 -0.199633 0.000192 * * * ******************************************************************************* ******************************************************************************* * * * EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL * * * * ALPHA =270.00000 BETA = 0.00000 GAMMA = 0.00000 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 3.9261 X_RMS = 1.9595 Y_RMS = 2.2400 Z_RMS = 2.5606 * * * * R_GEO = 5.0686 X_GEO = 4.3815 Y_GEO = 5.0088 Z_GEO = 5.7258 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * ******************************************************************************* * * * Q00 = 64.0000 ............. ............. ............. ............. * * * * Q10 = ZERO Q11 = ZERO ............. ............. ............. * * * * Q20 = 2.7241 Q21 = ZERO Q22 = -1.3060 ............. ............. * * * * Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. * * * * Q40 = -0.0073 Q41 = ZERO Q42 = -0.0400 Q43 = ZERO Q44 = 0.0024 * * * * Q50 = ZERO Q51 = ZERO Q52 = ZERO Q53 = ZERO Q54 = ZERO * * * * Q55 = ZERO * * * * Q60 = -0.0035 Q61 = ZERO Q62 = 5.8E-04 Q63 = ZERO Q64 =-5.4E-04 * * * * Q65 = ZERO Q66 =-5.2E-05 * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (EXACT MULTIPOLE MOMENTS) TOTAL * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = 0.2182 B21 = ZERO B22 = -0.1232 ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 = -0.0300 B41 = ZERO B42 = -0.0336 B43 = ZERO B44 = -0.0040 * * * * B50 = ZERO B51 = ZERO B52 = ZERO B53 = ZERO B54 = ZERO * * * * B55 = ZERO * * * * B60 = -0.0068 B61 = ZERO B62 = 0.0157 B63 = ZERO B64 = -0.0110 * * * * B65 = ZERO B66 =-8.1E-04 * * * ******************************************************************************* * * * BOHR DENSITY = 0.1218 FM^(-3) BOHR RADIUS = 4.9797 FM * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) TOTAL * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = 0.2189 B21 = ZERO B22 = -0.1049 ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 = -0.0072 B41 = ZERO B42 = -0.0561 B43 = ZERO B44 = 0.0033 * * * * B50 = ZERO B51 = ZERO B52 = ZERO B53 = ZERO B54 = ZERO * * * * B55 = ZERO * * * * B60 = -0.0133 B61 = ZERO B62 = 0.0032 B63 = ZERO B64 = -0.0029 * * * * B65 = ZERO B66 =-2.9E-04 * * * ******************************************************************************* * * * BOHR DENSITY = 0.1173 FM^(-3) BOHR RADIUS = 5.0686 FM * * * ******************************************************************************* ******************************************************************************* * * * ANGULAR MOMENTA AND THE FIRST MOMENTS OF INERTIA * * * ******************************************************************************* * * * SPINS J(1) * * --------------------------- --------------------------- * * ORBITAL INTRINSIC TOTAL ORBITAL INTRINSIC TOTAL * * * * NEUTRONS 2.26052 0.05416 2.31468 3.93133 0.09419 4.02552 * * PROTONS 2.31218 0.05460 2.36678 4.02118 0.09495 4.11613 * * -------- * * TOTAL 4.57269 0.10876 4.68145 7.95251 0.18915 8.14166 * * * ******************************************************************************* ******************************************************************************* * * * NEUTRON CONFIGURATIONS * * ======================= * * P S 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 * * --- -------------------------------------------------------------- * * * * CONF: + + 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + + 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: + - 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + - 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - + 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - + 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - - 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * PROTON CONFIGURATIONS * * ======================= * * P S 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 * * --- -------------------------------------------------------------- * * * * CONF: + + 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + + 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: + - 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + - 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - + 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - + 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - - 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * ENERGIES (MEV) * * * ******************************************************************************* * * * KINETIC: (NEU)= 560.156254 (PRO)= 543.823448 (TOT)= 1103.979702 * * SUM EPS: (NEU)= -863.426775 (PRO)= -535.818138 (TOT)= -1399.244913 * * 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)= -12.593997 (PRO)= -2.791129 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 176.252661 (EXC)= -12.344494 (TOT)= 163.908167 * * * * CONSTR. (MULT)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 * * CONSTR. (SURF)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 * * CONSTR. (SPIN)= -2.691836 SLOPE= 0.575000 CORR.= -1.345918 * * * * REARRANGEMENT ENERGY FROM THE SKYRME DENSITY-DEPENDENT TERMS= 388.849226 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -541.988394 * * * * SPIN-ORB (EVE)= -34.279829 (ODD)= -0.062401 (TOT)= -34.342229 * * SKYRME: (EVE)= -1806.984987 (ODD)= -0.199441 (TOT)= -1807.184427 * * * * TOTAL: (STAB)= 0.045813 (SP)= -539.250745 (FUN)= -539.296558 * * * ******************************************************************************* ******************************************************************************* * * * ITERATION RESTARTED FROM THE DISC, THE FOLLOWING FILE(S) READ: * * ge064-a.rec * * * ******************************************************************************* * * * CONVERGENCE REPORT * * * ******************************************************************************* * * * ITER ENERGY STABILITY Q_2 GAMMA SPIN OMEGA ANGLE EPAIR * * * * 47 -539.296558 0.045813 3.021 25.61 4.68 0.575 0.000 0.00 * * 48 -539.194158 0.128136 3.009 25.60 5.01 0.575 0.000 0.00 * * 49 -538.979910 -0.010588 2.986 25.61 5.50 0.575 0.000 0.00 * * 50 -538.856772 0.014367 2.967 25.58 5.74 0.575 0.000 0.00 * * 51 -538.784952 -0.000621 2.948 25.55 5.88 0.575 0.000 0.00 * * 52 -538.753712 -0.001574 2.930 25.50 5.94 0.575 0.000 0.00 * * 53 -538.746304 -0.005310 2.911 25.41 5.96 0.575 0.000 0.00 * * 54 -538.744729 -0.008603 2.899 25.32 5.97 0.575 0.000 0.00 * * 55 -538.736179 -0.005238 2.890 25.22 5.98 0.575 0.000 0.00 * * 56 -538.728987 -0.001724 2.885 25.15 5.99 0.575 0.000 0.00 * * 57 -538.727068 0.004782 2.884 25.11 6.00 0.575 0.000 0.00 * * 58 -538.727490 0.004508 2.884 25.08 6.00 0.575 0.000 0.00 * * 59 -538.728065 0.001703 2.884 25.05 6.00 0.575 0.000 0.00 * * 60 -538.728231 -0.000074 2.884 25.04 6.00 0.575 0.000 0.00 * * 61 -538.727673 -0.000693 2.884 25.04 6.00 0.575 0.000 0.00 * * 62 -538.727187 -0.000225 2.884 25.03 6.00 0.575 0.000 0.00 * * 63 -538.726924 0.000109 2.884 25.03 6.00 0.575 0.000 0.00 * * 64 -538.726870 0.000106 2.884 25.03 6.00 0.575 0.000 0.00 * * 65 -538.726889 0.000017 2.884 25.03 6.00 0.575 0.000 0.00 * * 66 -538.726889 -0.000043 2.884 25.03 6.00 0.575 0.000 0.00 * * 67 -538.726880 -0.000027 2.884 25.03 6.00 0.575 0.000 0.00 * * 68 -538.726865 3.62E-06 2.884 25.03 6.00 0.575 0.000 0.00 * * 69 -538.726857 0.000015 2.884 25.03 6.00 0.575 0.000 0.00 * * 70 -538.726854 9.34E-06 2.884 25.03 6.00 0.575 0.000 0.00 * * 71 -538.726853 -1.29E-06 2.884 25.03 6.00 0.575 0.000 0.00 * * 72 -538.726852 -4.38E-06 2.884 25.03 6.00 0.575 0.000 0.00 * * 73 -538.726850 -8.91E-07 2.884 25.03 6.00 0.575 0.000 0.00 * * 74 -538.726848 2.13E-06 2.884 25.03 6.00 0.575 0.000 0.00 * * 75 -538.726847 2.23E-06 2.884 25.03 6.00 0.575 0.000 0.00 * * 76 -538.726847 5.26E-07 2.884 25.03 6.00 0.575 0.000 0.00 * * 77 -538.726847 -5.00E-07 2.884 25.03 6.00 0.575 0.000 0.00 * * 78 -538.726847 -4.26E-07 2.884 25.03 6.00 0.575 0.000 0.00 * * 79 -538.726847 -8.92E-08 2.884 25.03 6.00 0.575 0.000 0.00 * * 80 -538.726847 9.62E-08 2.884 25.03 6.00 0.575 0.000 0.00 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 5) -39.297 ( 0, 0, 0, 2) | 1, 0, 1, 3/2> -100 -0.108 0.145 -1.338 * * 6) -39.113 ( 0, 0, 2, 0) | 1, 0, 1, 3/2> -100 -0.359 -0.316 0.881 * * 7) -36.812 ( 0, 0, 0, 3) | 1, 0, 1, 1/2> -100 -0.422 -0.423 1.002 * * 8) -36.589 ( 0, 0, 3, 0) | 1, 0, 1, 1/2> -100 -0.277 0.328 -1.181 * * 9) -31.990 ( 0, 2, 0, 0) | 2, 2, 0, 1/2> 100 1.968 0.469 0.238 * * 10) -31.324 ( 2, 0, 0, 0) | 2, 1, 1, 3/2> 100 0.899 0.184 0.205 * * 11) -30.142 ( 0, 3, 0, 0) | 2, 1, 1, 3/2> 100 0.086 -0.171 -1.983 * * 12) -29.253 ( 3, 0, 0, 0) | 2, 2, 0, 1/2> 100 -0.788 -0.267 0.339 * * 13) -27.897 ( 4, 0, 0, 0) | 2, 0, 2, 5/2> 100 -0.389 -0.129 0.332 * * 14) -27.892 ( 0, 4, 0, 0) | 2, 0, 2, 5/2> 100 -0.492 0.019 -0.039 * * 15) -26.861 ( 5, 0, 0, 0) | 2, 1, 1, 1/2> 100 -0.045 -0.379 8.334 * * 16) -26.069 ( 0, 5, 0, 0) | 2, 1, 1, 1/2> 100 -0.202 0.281 -1.389 * * 17) -24.299 ( 0, 6, 0, 0) | 2, 0, 2, 3/2> 100 0.160 0.010 0.065 * * 18) -23.781 ( 6, 0, 0, 0) | 2, 0, 2, 3/2> 100 -0.367 0.028 -0.075 * * 19) -21.922 ( 0, 0, 4, 0) | 3, 3, 0, 1/2> -100 2.728 0.446 0.163 * * 20) -21.530 ( 0, 0, 0, 4) | 3, 2, 1, 3/2> -100 1.641 0.280 0.171 * * 21) -21.253 ( 7, 0, 0, 0) | 2, 0, 0, 1/2> 100 -0.282 -0.436 1.547 * * 22) -21.198 ( 0, 7, 0, 0) | 2, 0, 0, 1/2> 100 -0.501 0.392 -0.781 * * 23) -20.498 ( 0, 0, 5, 0) | 3, 2, 1, 3/2> -100 0.968 0.064 0.066 * * 24) -19.321 ( 0, 0, 0, 5) | 3, 1, 2, 5/2> -100 0.025 -0.085 -3.388 * * 25) -18.353 ( 0, 0, 6, 0) | 3, 1, 2, 5/2> -100 -0.770 -0.137 0.177 * * 26) -17.733 ( 0, 0, 0, 6) | 3, 3, 0, 1/2> -100 -2.044 -0.361 0.177 * * 27) -16.318 ( 0, 0, 7, 0) | 3, 0, 3, 7/2> -100 -0.822 -0.078 0.095 * * 28) -16.300 ( 0, 0, 0, 7) | 3, 0, 3, 7/2> -100 -0.924 -0.081 0.088 * * 29) -15.313 ( 0, 0, 0, 8) | 3, 2, 1, 1/2> -100 1.075 -0.365 -0.340 * * 30) -14.763 ( 0, 0, 8, 0) | 3, 3, 0, 1/2> -100 0.455 0.392 0.863 * * 31) -14.103 ( 0, 0, 9, 0) | 3, 1, 2, 3/2> -100 0.642 -0.089 -0.138 * * 32) -13.493 ( 0, 0, 0, 9) | 3, 1, 0, 1/2> -100 -0.067 0.081 -1.203 * * 33) -11.584 ( 0, 8, 0, 0) | 4, 3, 1, 1/2> 100 3.393 0.427 0.126 * * 34) -11.371 ( 8, 0, 0, 0) | 4, 3, 1, 3/2> 100 2.357 0.311 0.132 * * 35) -11.309 ( 0, 0, 0,10) | 3, 1, 2, 3/2> -100 0.261 -0.258 -0.989 * * 36) -10.989 ( 0, 0, 0,11) | 3, 0, 1, 1/2> -100 0.127 0.110 0.863 * * 37) -10.836 ( 0, 0,10, 0) | 3, 0, 1, 3/2> -100 -0.414 0.208 -0.501 * * 38) -10.501 ( 0, 9, 0, 0) | 4, 2, 2, 5/2> 100 1.948 0.205 0.105 * * 39) -10.396 ( 0, 0,11, 0) | 3, 0, 1, 1/2> -100 -0.679 -0.164 0.241 * * 40) -9.478 ( 9, 0, 0, 0) | 4, 2, 2, 5/2> 100 0.855 0.043 0.050 * * 41) -9.200 ( 0, 0,12, 0) | 3, 0, 3, 5/2> -100 -0.044 0.007 -0.163 * * 42) -8.750 ( 0, 0, 0,12) | 3, 0, 3, 5/2> -100 -0.725 -0.012 0.016 * * 43) -8.411 ( 0,10, 0, 0) | 4, 1, 3, 7/2> 100 -0.170 -0.074 0.436 * * 44) -7.394 (10, 0, 0, 0) | 4, 1, 3, 7/2> 100 -1.024 -0.153 0.149 * * 45) -6.908 ( 0, 0, 0,13) | 3, 0, 1, 1/2> -100 -0.340 -0.309 0.909 * * 46) -6.627 ( 0,11, 0, 0) | 4, 1, 3, 7/2> 100 -1.554 -0.203 0.130 * * 47) -6.260 ( 0, 0,13, 0) | 3, 0, 1, 1/2> -100 -1.031 0.350 -0.340 * * 48) -6.210 (11, 0, 0, 0) | 4, 4, 0, 1/2> 100 -2.951 -0.381 0.129 * * 49) -4.514 ( 0,12, 0, 0) | 4, 0, 4, 9/2> 100 -0.406 0.051 -0.126 * * 50) -4.500 (12, 0, 0, 0) | 4, 0, 4, 9/2> 100 -0.967 -0.066 0.068 * * 51) -4.449 ( 0,13, 0, 0) | 4, 4, 0, 1/2> 100 0.862 0.328 0.381 * * 52) -4.089 (13, 0, 0, 0) | 4, 4, 0, 1/2> 100 0.874 -0.318 -0.364 * * 53) -3.298 (14, 0, 0, 0) | 4, 3, 1, 1/2> 100 1.035 0.222 0.214 * * 54) -2.972 ( 0,14, 0, 0) | 4, 2, 2, 3/2> 100 0.692 -0.271 -0.391 * * 55) -1.433 ( 0,15, 0, 0) | 4, 1, 1, 3/2> 100 0.387 -0.047 -0.120 * * 56) -1.288 (15, 0, 0, 0) | 4, 2, 0, 1/2> 100 0.526 -0.100 -0.190 * * 57) -1.116 ( 0, 0,14, 0) | 5, 3, 0, 1/2> -100 4.030 0.416 0.103 * * 58) -0.987 ( 0, 0, 0,14) | 5, 4, 1, 3/2> -100 3.110 0.332 0.107 * * 59) -0.759 (16, 0, 0, 0) | 4, 2, 2, 3/2> 100 0.690 -0.336 -0.487 * * 60) -0.495 ( 0,16, 0, 0) | 4, 1, 1, 1/2> 100 -0.130 0.390 -3.004 * * 61) -0.257 ( 0, 0,15, 0) | 5, 3, 2, 5/2> -100 2.942 0.282 0.096 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 1) -38.819 ( 0, 1, 0, 0) | 0, 0, 0, 1/2> 100 0.500 0.500 0.999 * * 2) -38.017 ( 1, 0, 0, 0) | 0, 0, 0, 1/2> 100 -0.496 -0.500 1.009 * * 3) -30.876 ( 0, 0, 1, 0) | 1, 1, 0, 1/2> -100 1.135 0.489 0.430 * * 4) -29.946 ( 0, 0, 0, 1) | 1, 1, 0, 1/2> -100 0.026 -0.238 -9.145 * * 5) -28.778 ( 0, 0, 0, 2) | 1, 0, 1, 3/2> -100 -0.073 0.166 -2.261 * * 6) -28.569 ( 0, 0, 2, 0) | 1, 0, 1, 3/2> -100 -0.360 -0.321 0.893 * * 7) -26.251 ( 0, 0, 0, 3) | 1, 0, 1, 1/2> -100 -0.436 -0.427 0.981 * * 8) -26.065 ( 0, 0, 3, 0) | 1, 0, 1, 1/2> -100 -0.269 0.332 -1.237 * * 9) -21.701 ( 0, 2, 0, 0) | 2, 2, 0, 1/2> 100 1.963 0.467 0.238 * * 10) -21.064 ( 2, 0, 0, 0) | 2, 1, 1, 3/2> 100 0.898 0.188 0.210 * * 11) -19.894 ( 0, 3, 0, 0) | 2, 1, 1, 3/2> 100 0.082 -0.177 -2.172 * * 12) -18.991 ( 3, 0, 0, 0) | 2, 2, 0, 1/2> 100 -0.787 -0.269 0.341 * * 13) -17.605 ( 0, 4, 0, 0) | 2, 0, 2, 5/2> 100 -0.472 0.022 -0.048 * * 14) -17.597 ( 4, 0, 0, 0) | 2, 0, 2, 5/2> 100 -0.375 -0.123 0.327 * * 15) -16.621 ( 5, 0, 0, 0) | 2, 1, 1, 1/2> 100 -0.048 -0.381 7.995 * * 16) -15.827 ( 0, 5, 0, 0) | 2, 1, 1, 1/2> 100 -0.224 0.277 -1.237 * * 17) -14.007 ( 0, 6, 0, 0) | 2, 0, 2, 3/2> 100 0.137 0.011 0.082 * * 18) -13.534 ( 6, 0, 0, 0) | 2, 0, 2, 3/2> 100 -0.352 0.030 -0.084 * * 19) -11.903 ( 0, 0, 4, 0) | 3, 3, 0, 1/2> -100 2.717 0.444 0.163 * * 20) -11.535 ( 0, 0, 0, 4) | 3, 2, 1, 3/2> -100 1.642 0.282 0.172 * * 21) -10.949 ( 0, 7, 0, 0) | 2, 0, 0, 1/2> 100 -0.465 0.399 -0.859 * * 22) -10.946 ( 7, 0, 0, 0) | 2, 0, 0, 1/2> 100 -0.309 -0.445 1.441 * * 23) -10.512 ( 0, 0, 5, 0) | 3, 2, 1, 3/2> -100 0.985 0.067 0.068 * * 24) -9.349 ( 0, 0, 0, 5) | 3, 1, 2, 5/2> -100 0.031 -0.086 -2.768 * * 25) -8.375 ( 0, 0, 6, 0) | 3, 1, 2, 5/2> -100 -0.777 -0.140 0.180 * * 26) -7.718 ( 0, 0, 0, 6) | 3, 3, 0, 1/2> -100 -2.061 -0.361 0.175 * * 27) -6.315 ( 0, 0, 7, 0) | 3, 0, 3, 7/2> -100 -0.819 -0.076 0.093 * * 28) -6.296 ( 0, 0, 0, 7) | 3, 0, 3, 7/2> -100 -0.930 -0.080 0.086 * * 29) -5.383 ( 0, 0, 0, 8) | 3, 2, 1, 1/2> -100 1.125 -0.362 -0.322 * * 30) -4.799 ( 0, 0, 8, 0) | 3, 3, 0, 1/2> -100 0.432 0.372 0.861 * * 31) -4.207 ( 0, 0, 9, 0) | 3, 1, 2, 3/2> -100 0.674 -0.065 -0.097 * * 32) -3.579 ( 0, 0, 0, 9) | 3, 1, 0, 1/2> -100 -0.081 0.075 -0.927 * * 33) -1.886 ( 0, 8, 0, 0) | 4, 3, 1, 1/2> 100 3.390 0.427 0.126 * * 34) -1.683 ( 8, 0, 0, 0) | 4, 3, 1, 3/2> 100 2.369 0.313 0.132 * * 35) -1.467 ( 0, 0, 0,10) | 3, 1, 2, 3/2> -100 0.270 -0.260 -0.960 * * 36) -1.142 ( 0, 0, 0,11) | 3, 0, 1, 1/2> -100 0.156 0.116 0.739 * * 37) -0.979 ( 0, 0,10, 0) | 3, 0, 1, 3/2> -100 -0.435 0.245 -0.562 * * 38) -0.826 ( 0, 9, 0, 0) | 4, 2, 2, 5/2> 100 1.963 0.207 0.106 * * 39) -0.564 ( 0, 0,11, 0) | 3, 0, 1, 1/2> -100 -0.656 -0.201 0.306 * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 6.951690 3.476935 6.951690 0.002179 * * DRHOD = 0.889465 0.444872 0.889465 0.000279 * * DLPR_ = -2.205403 -1.103605 -2.205403 -0.001808 * * DTAU_ = 5.686796 2.844702 5.686796 0.002609 * * DSCU_ = 0.059089 0.029569 0.059089 0.000049 * * DDIV_ = 0.385001 0.192648 0.385001 0.000294 * * * * DSPI_ = 0.005459 0.002732 0.005459 0.000005 * * DSPID = 0.000747 0.000374 0.000747 0.000001 * * DLPS_ = -0.013795 -0.006907 -0.013795 -0.000019 * * DCUR_ = 0.007976 0.003990 0.007976 0.000005 * * DKIS_ = 0.001997 0.000999 0.001997 0.000000 * * DROT_ = 0.004163 0.002083 0.004163 0.000003 * * * ******************************************************************************* ******************************************************************************* ******************************************************************************* * * * CONTRIBUTIONS TO ENERGY IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * ERHO_ = -4806.116086 1864.180234 -2942.520053 0.584200 * * ERHOD = 1556.564306 -778.526631 778.282153 -0.244478 * * ELPR_ = 176.432206 -37.591557 138.871443 -0.030794 * * ETAU_ = 426.509687 -174.238021 252.351565 -0.079899 * * ESCU_ = 0.000000 0.000000 0.000000 0.000000 * * EDIV_ = -23.100076 -11.558870 -34.650113 -0.008832 * * ============ ============ ============ ============ * * SUM EVEN: -2669.709962 862.265155 -1807.665005 0.220198 * * * * ESPI_ = -0.693172 0.770832 0.077019 0.000641 * * ESPID = 0.435886 -0.218122 0.217943 -0.000179 * * ELPS_ = 0.000000 -0.235269 -0.234938 -0.000331 * * ECUR_ = -0.598189 0.244403 -0.353928 0.000142 * * EKIS_ = 0.000000 0.000000 0.000000 0.000000 * * EROT_ = -0.249794 -0.124983 -0.374692 -0.000086 * * ============ ============ ============ ============ * * SUM ODD: -1.105269 0.436860 -0.668596 0.000187 * * * ******************************************************************************* ******************************************************************************* * * * EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL * * * * ALPHA =270.00000 BETA = 0.00000 GAMMA = 0.00000 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 3.9215 X_RMS = 1.9736 Y_RMS = 2.2351 Z_RMS = 2.5469 * * * * R_GEO = 5.0626 X_GEO = 4.4132 Y_GEO = 4.9979 Z_GEO = 5.6951 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * ******************************************************************************* * * * Q00 = 64.0000 ............. ............. ............. ............. * * * * Q10 = ZERO Q11 = ZERO ............. ............. ............. * * * * Q20 = 2.6128 Q21 = ZERO Q22 = -1.2201 ............. ............. * * * * Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. * * * * Q40 = -0.0099 Q41 = ZERO Q42 = -0.0394 Q43 = ZERO Q44 = 0.0013 * * * * Q50 = ZERO Q51 = ZERO Q52 = ZERO Q53 = ZERO Q54 = ZERO * * * * Q55 = ZERO * * * * Q60 = -0.0032 Q61 = ZERO Q62 = 4.9E-04 Q63 = ZERO Q64 =-4.6E-04 * * * * Q65 = ZERO Q66 =-1.7E-05 * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (EXACT MULTIPOLE MOMENTS) TOTAL * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = 0.2094 B21 = ZERO B22 = -0.1138 ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 = -0.0306 B41 = ZERO B42 = -0.0348 B43 = ZERO B44 = -0.0049 * * * * B50 = ZERO B51 = ZERO B52 = ZERO B53 = ZERO B54 = ZERO * * * * B55 = ZERO * * * * B60 = -0.0054 B61 = ZERO B62 = 0.0147 B63 = ZERO B64 = -0.0098 * * * * B65 = ZERO B66 =-7.7E-04 * * * ******************************************************************************* * * * BOHR DENSITY = 0.1219 FM^(-3) BOHR RADIUS = 4.9815 FM * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) TOTAL * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = 0.2104 B21 = ZERO B22 = -0.0983 ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 = -0.0098 B41 = ZERO B42 = -0.0556 B43 = ZERO B44 = 0.0018 * * * * B50 = ZERO B51 = ZERO B52 = ZERO B53 = ZERO B54 = ZERO * * * * B55 = ZERO * * * * B60 = -0.0123 B61 = ZERO B62 = 0.0027 B63 = ZERO B64 = -0.0025 * * * * B65 = ZERO B66 =-9.6E-05 * * * ******************************************************************************* * * * BOHR DENSITY = 0.1178 FM^(-3) BOHR RADIUS = 5.0626 FM * * * ******************************************************************************* ******************************************************************************* * * * ANGULAR MOMENTA AND THE FIRST MOMENTS OF INERTIA * * * ******************************************************************************* * * * SPINS J(1) * * --------------------------- --------------------------- * * ORBITAL INTRINSIC TOTAL ORBITAL INTRINSIC TOTAL * * * * NEUTRONS 2.91368 0.06726 2.98094 5.06727 0.11697 5.18424 * * PROTONS 2.94839 0.06862 3.01702 5.12764 0.11934 5.24698 * * -------- * * TOTAL 5.86208 0.13588 5.99795 10.19491 0.23631 10.43123 * * * ******************************************************************************* ******************************************************************************* * * * NEUTRON CONFIGURATIONS * * ======================= * * P S 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 * * --- -------------------------------------------------------------- * * * * CONF: + + 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + + 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: + - 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + - 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - + 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - + 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - - 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * PROTON CONFIGURATIONS * * ======================= * * P S 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 * * --- -------------------------------------------------------------- * * * * CONF: + + 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + + 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: + - 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + - 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - + 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - + 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - - 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * ENERGIES (MEV) * * * ******************************************************************************* * * * KINETIC: (NEU)= 560.882688 (PRO)= 544.495493 (TOT)= 1105.378180 * * SUM EPS: (NEU)= -863.767069 (PRO)= -536.028071 (TOT)= -1399.795140 * * 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)= -12.538534 (PRO)= -2.732136 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 176.353366 (EXC)= -12.345178 (TOT)= 164.008189 * * * * CONSTR. (MULT)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 * * CONSTR. (SURF)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 * * CONSTR. (SPIN)= -3.448824 SLOPE= 0.575000 CORR.= -1.724412 * * * * REARRANGEMENT ENERGY FROM THE SKYRME DENSITY-DEPENDENT TERMS= 389.127720 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -542.175671 * * * * SPIN-ORB (EVE)= -34.658945 (ODD)= -0.374778 (TOT)= -35.033723 * * SKYRME: (EVE)= -1807.444807 (ODD)= -0.668409 (TOT)= -1808.113216 * * * * TOTAL: (STAB)= 9.89E-08 (SP)= -538.726847 (FUN)= -538.726847 * * * ******************************************************************************* ******************************************************************************* * * * CALLING SUBROUTINE "RECORD" TO SAVE THE RECORD FILE CONTAINING RESULTS * * ge064-a.rec * * * ******************************************************************************* ******************************************************************************* * * * CALCULATION OF THE DATA SET NO. 3 HAS BEEN SUCCESSFULLY COMPLETED * * * ******************************************************************************* ******************************************************************************* * * * CLASSICAL NUCLEAR SURFACE DEFINED FOR: N = 32 Z = 32 * * * ******************************************************************************* * * * SPHERICAL SHAPE REQUESTED FOR THE SURFACE * * * ******************************************************************************* * * * HOMEGA= 12.3000 FCHOM0= 1.2000 * * * * OSCILLATOR FREQUENCIES: HBAROX= 12.3000 HBAROY= 12.3000 HBAROZ= 12.3000 * * * * MOMENTS OF INERTIA: XMOMFC= 14.9370 YMOMFC= 14.9370 ZMOMFC= 14.9370 * * * * CENTERS OF MASS: CMSXFC= 0.0000 CMSYFC= 0.0000 CMSZFC= 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * SKYRME FORCE DEFINITION * * * ******************************************************************************* * * * PARAMETER SET SIII: T0= -1128.75 T1= 395.00 T2= -95.00 T3= 14000.00 * * * * POWER=1/1 W=120.000 X0= 0.45000 X1= 0.00000 X2= 0.00000 X3= 1.00000 * * * * ETA_J=0 ETA_W=0 ETACM=0 ETA_M=1 HBM= 20.73000 (SKYRME-FORCE SPIN-ORBIT) * * * ******************************************************************************* ******************************************************************************* * * * COUPLING CONSTANTS DEFINING THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * CRHO_ = -691.359375 536.156250 -423.281250 268.078125 * * CRHOD = 1750.000000 -1750.000000 875.000000 -875.000000 * * CLPR_ = -80.000000 34.062500 -62.968750 17.031250 * * CTAU_ = 75.000000 -61.250000 44.375000 -30.625000 * * CSCU_ = 0.000000 0.000000 0.000000 0.000000 * * CDIV_ = -60.000000 -60.000000 -90.000000 -30.000000 * * * * CSPI_ = -126.984375 282.187500 14.109375 141.093750 * * CSPID = 583.333333 -583.333333 291.666667 -291.666667 * * CLPS_ = 0.000000 34.062500 17.031250 17.031250 * * CCUR_ = -75.000000 61.250000 -44.375000 30.625000 * * CKIS_ = 0.000000 0.000000 0.000000 0.000000 * * CROT_ = -60.000000 -60.000000 -90.000000 -30.000000 * * * ******************************************************************************* ******************************************************************************* * * * PHYSICAL CONSTANTS: H_BARC=197.32891000 HBCOE2=137.03602000 * * * * XMASSN=938.90590000 XMASSP=938.27231000 * * * * HBMASS= 20.73000000 HBMRPA= 20.40609375 * * * * ECHAR2= 1.43997841 COULEX= -1.06350868 * * * ******************************************************************************* * * * OSCILLATOR LENGTHS: X= 1.8362302 Y= 1.8362302 Z= 1.8362302 * * * * OSCILLATOR CONSTANTS: X= 0.5445940 Y= 0.5445940 Z= 0.5445940 * * * * OSCILLATOR FREQUENCIES: X=12.3000000 Y=12.3000000 Z=12.3000000 * * * ******************************************************************************* * * * BASIS CUT-OFF CONTROL PARAMETERS: NXMAXX = 14 NYMAXX = 14 NZMAXX = 14 * * * * OPTIMUM NUMBERS OF GAUSS POINTS: NXHERM = 30 NYHERM = 30 NZHERM = 30 * * * * NLIMIT= 680 LDBASE= 680 MCOUNT= 3375 * * * * ENECUT= 800.0000 ELIMIT= 190.6500 * * * ******************************************************************************* ******************************************************************************* * * * SHAPE OF THE OSCILLATOR-BASIS DIAMOND (THREE CONSECUTIVE PROJECTIONS) * * * ******************************************************************************* * * * NX ===>>> 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 * * MAX.NY => 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------------------- * * NY= 0 (14) | 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 1 (13) | 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 2 (12) | 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 3 (11) | 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 4 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NY= 5 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NY= 6 ( 8) | 8 7 6 5 4 3 2 1 0 * * NY= 7 ( 7) | 7 6 5 4 3 2 1 0 * * NY= 8 ( 6) | 6 5 4 3 2 1 0 * * NY= 9 ( 5) | 5 4 3 2 1 0 * * NY=10 ( 4) | 4 3 2 1 0 * * NY=11 ( 3) | 3 2 1 0 * * NY=12 ( 2) | 2 1 0 * * NY=13 ( 1) | 1 0 * * NY=14 ( 0) | 0 * * * * * * NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 * * MAX.NY => 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------------------- * * NY= 0 (14) | 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 1 (13) | 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 2 (12) | 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 3 (11) | 11 10 9 8 7 6 5 4 3 2 1 0 * * NY= 4 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NY= 5 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NY= 6 ( 8) | 8 7 6 5 4 3 2 1 0 * * NY= 7 ( 7) | 7 6 5 4 3 2 1 0 * * NY= 8 ( 6) | 6 5 4 3 2 1 0 * * NY= 9 ( 5) | 5 4 3 2 1 0 * * NY=10 ( 4) | 4 3 2 1 0 * * NY=11 ( 3) | 3 2 1 0 * * NY=12 ( 2) | 2 1 0 * * NY=13 ( 1) | 1 0 * * NY=14 ( 0) | 0 * * * * * * NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 * * MAX.NX => 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------------------- * * NX= 0 (14) | 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 1 (13) | 13 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 2 (12) | 12 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 3 (11) | 11 10 9 8 7 6 5 4 3 2 1 0 * * NX= 4 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NX= 5 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NX= 6 ( 8) | 8 7 6 5 4 3 2 1 0 * * NX= 7 ( 7) | 7 6 5 4 3 2 1 0 * * NX= 8 ( 6) | 6 5 4 3 2 1 0 * * NX= 9 ( 5) | 5 4 3 2 1 0 * * NX=10 ( 4) | 4 3 2 1 0 * * NX=11 ( 3) | 3 2 1 0 * * NX=12 ( 2) | 2 1 0 * * NX=13 ( 1) | 1 0 * * NX=14 ( 0) | 0 * * * ******************************************************************************* ******************************************************************************* * * * AUTOMATED HANDLING OF THE COULOMB FILE * * * ******************************************************************************* * * * THE COULOMB FILE READ FROM THE DISC, FILE NAME (FILCOU): * * ge064-a.cou * * * //////////////////// INCORRECT NUMBER OF MULTIPOLES //////////////////// //////////////////// FOUND ON THE COULOMB FILE: //////////////////// //////////////////// READ SHOULD BE //////////////////// //////////////////// NUPOLS= 9 25 //////////////////// ******************************************************************************* * * * THE COULOMB DATA HAVE BEEN RECALCULATED * * * ******************************************************************************* * * * THE COULOMB FILE WRITTEN ON THE DISC, FILE NAME (FILCOU): * * ge064-a.cou * * * ******************************************************************************* * * * NUCLIDE: N = 32 Z = 32 NUMBER OF ITERATIONS = 1 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 LIPKIN-NOGAMI = 0.50 * * * ******************************************************************************* * * * DIRECT COULOMB TERMS SOLVED BY GREEN FUNCTION AND FOURIER TRANSFORM METHOD * * COULOMB PARAMETERS: NUMCOU = 80 NUMETA = 79 BOUCOU = 20.0 FURMAX =0.25 * * * ******************************************************************************* * * * EXCHANGE COULOMB TERMS TREATED BY THE SLATER APPROXIMATION * * * ******************************************************************************* * * * MAXIMUM NUMBERS OF MULTIPOLE MOMENTS USED FOR THE CONSTRAINTS, NMUCON=2 * * FOR THE COULOMB FIELD, NMUCOU=4 * * FOR THE OUTPUT INFO, NMUPRI=6 * * * ******************************************************************************* * * * 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: NO * * PROTONS : NO * * MASS : YES * * * ******************************************************************************* * * * PRINTING THE MOMENTS, MOMENTA AND DEFORMATIONS IN THE INTRINSIC FRAME: NO * * * ******************************************************************************* * * * NILSSON LABELS OF SINGLE-PARTICLE STATES SHOWN WITH RESPECT TO THE: Z-AXIS * * * ******************************************************************************* * * * NUMBERS OF LEVELS IN TIME-REVERSAL/CHARGE BLOCKS: TIME-UP TIME-DOWN * * NEUTRONS: 32 32 * * PROTONS: 32 32 * * * ******************************************************************************* * * * THE ANGULAR MOMENTUM PROJECTION IS PERFORMED AFTER THE ITERATIONS ARE DONE * * * * PROJECTION ON VALUES OF ANGULAR MOMENTA BETWEEN J = 0 AND J = 5 HBAR * * * * NUMBERS OF INTEGRATION POINTS IN EULER ANGLES: 10 (ALPHA/GAMMA), 10 (BETA) * * * ******************************************************************************* * * * REDUCED KERNELS AND MATRIX ELEMENTS: * * FOR ELECTRIC MULTIPOLES CALCULATED UP TO = 2, PRINTED IF LARGER THAN 0.001 * * FOR MAGNETIC MULTIPOLES CALCULATED UP TO = 1, PRINTED IF LARGER THAN 0.001 * * FOR SURFACE MULTIPOLES CALCULATED UP TO = 0, PRINTED IF LARGER THAN 0.001 * * * ******************************************************************************* * * * PRINTING THE RESULTS FOR ANGULAR-MOMENTUM-PROJECTED STATES: 0 < 2*I < 10 * * AVERAGE ENERGIES OF PROJECTED STATES: YES * * KERNELS BETWEEN THE PROJECTED STATES: NO * * SUM RULES COMPARED TO THE H-F VALUES: YES * * K-MIXED ENERGIES OF PROJECTED STATES: YES * * REDUCED K E R N E L S BETWEEN THE PROJECTED STATES: YES * * REDUCED MATRIX ELEMENTS BETWEEN THE K-MIXED STATES: YES * * * ******************************************************************************* * * * THE GCM KERNELS IN FUNCTION OF THE EULER ANGLES WILL BE READ AND/OR STORED * * * ******************************************************************************* * * * CALCULATIONS WITHOUT PAIRING * * * ******************************************************************************* * * * CALCULATIONS WITHOUT SIMPLEX, SIGNATURE, OR PARITY SYMMETRIES * * * ******************************************************************************* * * * 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 ONE-DIMENSIONAL CRANKING FOR THE ROTATIONAL FREQUENCIES: * * * * OMEGAY = 0.575000 (ISOSCALAR) AND OMISOY = 0.000000 (ISOVECTOR) * * * ******************************************************************************* * * * NO-SYMMETRY CONFIGURATIONS: * * * * V A C U U M P A R T I C L E S H O L E S * * =========== ================= =========== * * * * NEUTRONS: 32 0 0 * * PROTONS : 32 0 0 * * * ******************************************************************************* ******************************************************************************* * * * THE RECORD FILE WILL NOT BE SAVED; AND THIS RUN CANNOT BE LATER CONTINUED * * * ******************************************************************************* ******************************************************************************* * * * CALLING SUBROUTINE "RECORD" TO READ THE REPLAY FILE OBTAINED PREVIOUSLY * * ge064-a.rec * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY | N,nz,lz,OMEz>

* * * * 5) -39.297 1.0000 | 1, 0, 1, 3/2> -100 0.000 -0.108 0.000 * * 0.000 0.145 0.000 * * 6) -39.113 -1.0000 | 1, 0, 1, 3/2> -100 0.000 -0.359 0.000 * * 0.000 -0.316 0.000 * * 7) -36.812 1.0000 | 1, 0, 1, 1/2> -100 0.000 -0.422 0.000 * * 0.000 -0.423 0.000 * * 8) -36.589 -1.0000 | 1, 0, 1, 1/2> -100 0.000 -0.277 0.000 * * 0.000 0.328 0.000 * * 9) -31.990 -1.0000 | 2, 2, 0, 1/2> 100 0.000 1.968 0.000 * * 0.000 0.469 0.000 * * 10) -31.324 1.0000 | 2, 1, 1, 3/2> 100 0.000 0.899 0.000 * * 0.000 0.184 0.000 * * 11) -30.142 -1.0000 | 2, 1, 1, 3/2> 100 0.000 0.086 0.000 * * 0.000 -0.171 0.000 * * 12) -29.253 1.0000 | 2, 2, 0, 1/2> 100 0.000 -0.788 0.000 * * 0.000 -0.267 0.000 * * 13) -27.897 1.0000 | 2, 0, 2, 5/2> 100 0.000 -0.389 0.000 * * 0.000 -0.129 0.000 * * 14) -27.892 -1.0000 | 2, 0, 2, 5/2> 100 0.000 -0.492 0.000 * * 0.000 0.019 0.000 * * 15) -26.861 1.0000 | 2, 1, 1, 1/2> 100 0.000 -0.045 0.000 * * 0.000 -0.379 0.000 * * 16) -26.069 -1.0000 | 2, 1, 1, 1/2> 100 0.000 -0.202 0.000 * * 0.000 0.281 0.000 * * 17) -24.299 -1.0000 | 2, 0, 2, 3/2> 100 0.000 0.160 0.000 * * 0.000 0.010 0.000 * * 18) -23.781 1.0000 | 2, 0, 2, 3/2> 100 0.000 -0.367 0.000 * * 0.000 0.028 0.000 * * 19) -21.922 -1.0000 | 3, 3, 0, 1/2> -100 0.000 2.728 0.000 * * 0.000 0.446 0.000 * * 20) -21.530 1.0000 | 3, 2, 1, 3/2> -100 0.000 1.641 0.000 * * 0.000 0.280 0.000 * * 21) -21.253 1.0000 | 2, 0, 0, 1/2> 100 0.000 -0.282 0.000 * * 0.000 -0.436 0.000 * * 22) -21.198 -1.0000 | 2, 0, 0, 1/2> 100 0.000 -0.501 0.000 * * 0.000 0.392 0.000 * * 23) -20.498 -1.0000 | 3, 2, 1, 3/2> -100 0.000 0.968 0.000 * * 0.000 0.064 0.000 * * 24) -19.321 1.0000 | 3, 1, 2, 5/2> -100 0.000 0.025 0.000 * * 0.000 -0.085 0.000 * * 25) -18.353 -1.0000 | 3, 1, 2, 5/2> -100 0.000 -0.770 0.000 * * 0.000 -0.137 0.000 * * 26) -17.733 1.0000 | 3, 3, 0, 1/2> -100 0.000 -2.044 0.000 * * 0.000 -0.361 0.000 * * 27) -16.318 -1.0000 | 3, 0, 3, 7/2> -100 0.000 -0.822 0.000 * * 0.000 -0.078 0.000 * * 28) -16.300 1.0000 | 3, 0, 3, 7/2> -100 0.000 -0.924 0.000 * * 0.000 -0.081 0.000 * * 29) -15.313 1.0000 | 3, 2, 1, 1/2> -100 0.000 1.075 0.000 * * 0.000 -0.365 0.000 * * 30) -14.763 -1.0000 | 3, 3, 0, 1/2> -100 0.000 0.455 0.000 * * 0.000 0.392 0.000 * * 31) -14.103 -1.0000 | 3, 1, 2, 3/2> -100 0.000 0.642 0.000 * * 0.000 -0.089 0.000 * * 32) -13.493 1.0000 | 3, 1, 0, 1/2> -100 0.000 -0.067 0.000 * * 0.000 0.081 0.000 * * 33) -11.584 -1.0000 | 4, 3, 1, 1/2> 100 0.000 3.393 0.000 * * 0.000 0.427 0.000 * * 34) -11.371 1.0000 | 4, 3, 1, 3/2> 100 0.000 2.357 0.000 * * 0.000 0.311 0.000 * * 35) -11.309 1.0000 | 3, 1, 2, 3/2> -100 0.000 0.261 0.000 * * 0.000 -0.258 0.000 * * 36) -10.989 1.0000 | 3, 0, 1, 1/2> -100 0.000 0.127 0.000 * * 0.000 0.110 0.000 * * 37) -10.836 -1.0000 | 3, 0, 1, 3/2> -100 0.000 -0.414 0.000 * * 0.000 0.208 0.000 * * 38) -10.501 -1.0000 | 4, 2, 2, 5/2> 100 0.000 1.948 0.000 * * 0.000 0.205 0.000 * * 39) -10.396 -1.0000 | 3, 0, 1, 1/2> -100 0.000 -0.679 0.000 * * 0.000 -0.164 0.000 * * 40) -9.478 1.0000 | 4, 2, 2, 5/2> 100 0.000 0.855 0.000 * * 0.000 0.043 0.000 * * 41) -9.200 -1.0000 | 3, 0, 3, 5/2> -100 0.000 -0.044 0.000 * * 0.000 0.007 0.000 * * 42) -8.750 1.0000 | 3, 0, 3, 5/2> -100 0.000 -0.725 0.000 * * 0.000 -0.012 0.000 * * 43) -8.411 -1.0000 | 4, 1, 3, 7/2> 100 0.000 -0.170 0.000 * * 0.000 -0.074 0.000 * * 44) -7.394 1.0000 | 4, 1, 3, 7/2> 100 0.000 -1.024 0.000 * * 0.000 -0.153 0.000 * * 45) -6.908 1.0000 | 3, 0, 1, 1/2> -100 0.000 -0.340 0.000 * * 0.000 -0.309 0.000 * * 46) -6.627 -1.0000 | 4, 1, 3, 7/2> 100 0.000 -1.554 0.000 * * 0.000 -0.203 0.000 * * 47) -6.260 -1.0000 | 3, 0, 1, 1/2> -100 0.000 -1.031 0.000 * * 0.000 0.350 0.000 * * 48) -6.210 1.0000 | 4, 4, 0, 1/2> 100 0.000 -2.951 0.000 * * 0.000 -0.381 0.000 * * 49) -4.514 -1.0000 | 4, 0, 4, 9/2> 100 0.000 -0.406 0.000 * * 0.000 0.051 0.000 * * 50) -4.500 1.0000 | 4, 0, 4, 9/2> 100 0.000 -0.967 0.000 * * 0.000 -0.066 0.000 * * 51) -4.449 -1.0000 | 4, 4, 0, 1/2> 100 0.000 0.862 0.000 * * 0.000 0.328 0.000 * * 52) -4.089 1.0000 | 4, 4, 0, 1/2> 100 0.000 0.874 0.000 * * 0.000 -0.318 0.000 * * 53) -3.298 1.0000 | 4, 3, 1, 1/2> 100 0.000 1.035 0.000 * * 0.000 0.222 0.000 * * 54) -2.972 -1.0000 | 4, 2, 2, 3/2> 100 0.000 0.692 0.000 * * 0.000 -0.271 0.000 * * 55) -1.433 -1.0000 | 4, 1, 1, 3/2> 100 0.000 0.387 0.000 * * 0.000 -0.047 0.000 * * 56) -1.288 1.0000 | 4, 2, 0, 1/2> 100 0.000 0.526 0.000 * * 0.000 -0.100 0.000 * * 57) -1.116 -1.0000 | 5, 3, 0, 1/2> -100 0.000 4.030 0.000 * * 0.000 0.416 0.000 * * 58) -0.987 1.0000 | 5, 4, 1, 3/2> -100 0.000 3.110 0.000 * * 0.000 0.332 0.000 * * 59) -0.759 1.0000 | 4, 2, 2, 3/2> 100 0.000 0.690 0.000 * * 0.000 -0.336 0.000 * * 60) -0.495 -1.0000 | 4, 1, 1, 1/2> 100 0.000 -0.130 0.000 * * 0.000 0.390 0.000 * * 61) -0.257 -1.0000 | 5, 3, 2, 5/2> -100 0.000 2.942 0.000 * * 0.000 0.282 0.000 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS * * * ******************************************************************************* * * * NO) ENERGY | N,nz,lz,OMEz>

* * * * 1) -38.819 -1.0000 | 0, 0, 0, 1/2> 100 0.000 0.500 0.000 * * 0.000 0.500 0.000 * * 2) -38.017 1.0000 | 0, 0, 0, 1/2> 100 0.000 -0.496 0.000 * * 0.000 -0.500 0.000 * * 3) -30.876 -1.0000 | 1, 1, 0, 1/2> -100 0.000 1.135 0.000 * * 0.000 0.489 0.000 * * 4) -29.946 1.0000 | 1, 1, 0, 1/2> -100 0.000 0.026 0.000 * * 0.000 -0.238 0.000 * * 5) -28.778 1.0000 | 1, 0, 1, 3/2> -100 0.000 -0.073 0.000 * * 0.000 0.166 0.000 * * 6) -28.569 -1.0000 | 1, 0, 1, 3/2> -100 0.000 -0.360 0.000 * * 0.000 -0.321 0.000 * * 7) -26.251 1.0000 | 1, 0, 1, 1/2> -100 0.000 -0.436 0.000 * * 0.000 -0.427 0.000 * * 8) -26.065 -1.0000 | 1, 0, 1, 1/2> -100 0.000 -0.269 0.000 * * 0.000 0.332 0.000 * * 9) -21.701 -1.0000 | 2, 2, 0, 1/2> 100 0.000 1.963 0.000 * * 0.000 0.467 0.000 * * 10) -21.064 1.0000 | 2, 1, 1, 3/2> 100 0.000 0.898 0.000 * * 0.000 0.188 0.000 * * 11) -19.894 -1.0000 | 2, 1, 1, 3/2> 100 0.000 0.082 0.000 * * 0.000 -0.177 0.000 * * 12) -18.991 1.0000 | 2, 2, 0, 1/2> 100 0.000 -0.787 0.000 * * 0.000 -0.269 0.000 * * 13) -17.605 -1.0000 | 2, 0, 2, 5/2> 100 0.000 -0.472 0.000 * * 0.000 0.022 0.000 * * 14) -17.597 1.0000 | 2, 0, 2, 5/2> 100 0.000 -0.375 0.000 * * 0.000 -0.123 0.000 * * 15) -16.621 1.0000 | 2, 1, 1, 1/2> 100 0.000 -0.048 0.000 * * 0.000 -0.381 0.000 * * 16) -15.827 -1.0000 | 2, 1, 1, 1/2> 100 0.000 -0.224 0.000 * * 0.000 0.277 0.000 * * 17) -14.007 -1.0000 | 2, 0, 2, 3/2> 100 0.000 0.137 0.000 * * 0.000 0.011 0.000 * * 18) -13.534 1.0000 | 2, 0, 2, 3/2> 100 0.000 -0.352 0.000 * * 0.000 0.030 0.000 * * 19) -11.903 -1.0000 | 3, 3, 0, 1/2> -100 0.000 2.717 0.000 * * 0.000 0.444 0.000 * * 20) -11.535 1.0000 | 3, 2, 1, 3/2> -100 0.000 1.642 0.000 * * 0.000 0.282 0.000 * * 21) -10.949 -1.0000 | 2, 0, 0, 1/2> 100 0.000 -0.465 0.000 * * 0.000 0.399 0.000 * * 22) -10.946 1.0000 | 2, 0, 0, 1/2> 100 0.000 -0.309 0.000 * * 0.000 -0.445 0.000 * * 23) -10.512 -1.0000 | 3, 2, 1, 3/2> -100 0.000 0.985 0.000 * * 0.000 0.067 0.000 * * 24) -9.349 1.0000 | 3, 1, 2, 5/2> -100 0.000 0.031 0.000 * * 0.000 -0.086 0.000 * * 25) -8.375 -1.0000 | 3, 1, 2, 5/2> -100 0.000 -0.777 0.000 * * 0.000 -0.140 0.000 * * 26) -7.718 1.0000 | 3, 3, 0, 1/2> -100 0.000 -2.061 0.000 * * 0.000 -0.361 0.000 * * 27) -6.315 -1.0000 | 3, 0, 3, 7/2> -100 0.000 -0.819 0.000 * * 0.000 -0.076 0.000 * * 28) -6.296 1.0000 | 3, 0, 3, 7/2> -100 0.000 -0.930 0.000 * * 0.000 -0.080 0.000 * * 29) -5.383 1.0000 | 3, 2, 1, 1/2> -100 0.000 1.125 0.000 * * 0.000 -0.362 0.000 * * 30) -4.799 -1.0000 | 3, 3, 0, 1/2> -100 0.000 0.432 0.000 * * 0.000 0.372 0.000 * * 31) -4.207 -1.0000 | 3, 1, 2, 3/2> -100 0.000 0.674 0.000 * * 0.000 -0.065 0.000 * * 32) -3.579 1.0000 | 3, 1, 0, 1/2> -100 0.000 -0.081 0.000 * * 0.000 0.075 0.000 * * 33) -1.886 -1.0000 | 4, 3, 1, 1/2> 100 0.000 3.390 0.000 * * 0.000 0.427 0.000 * * 34) -1.683 1.0000 | 4, 3, 1, 3/2> 100 0.000 2.369 0.000 * * 0.000 0.313 0.000 * * 35) -1.467 1.0000 | 3, 1, 2, 3/2> -100 0.000 0.270 0.000 * * 0.000 -0.260 0.000 * * 36) -1.142 1.0000 | 3, 0, 1, 1/2> -100 0.000 0.156 0.000 * * 0.000 0.116 0.000 * * 37) -0.979 -1.0000 | 3, 0, 1, 3/2> -100 0.000 -0.435 0.000 * * 0.000 0.245 0.000 * * 38) -0.826 -1.0000 | 4, 2, 2, 5/2> 100 0.000 1.963 0.000 * * 0.000 0.207 0.000 * * 39) -0.564 -1.0000 | 3, 0, 1, 1/2> -100 0.000 -0.656 0.000 * * 0.000 -0.201 0.000 * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 6.951690 3.476935 6.951690 0.002179 * * DRHOD = 0.889465 0.444872 0.889465 0.000279 * * DLPR_ = -2.205403 -1.103605 -2.205403 -0.001808 * * DTAU_ = 5.686796 2.844702 5.686796 0.002609 * * DSCU_ = 0.059089 0.029569 0.059089 0.000049 * * DDIV_ = 0.385001 0.192648 0.385001 0.000294 * * * * DSPI_ = 0.005459 0.002732 0.005459 0.000005 * * DSPID = 0.000747 0.000374 0.000747 0.000001 * * DLPS_ = -0.013795 -0.006907 -0.013795 -0.000019 * * DCUR_ = 0.007976 0.003990 0.007976 0.000005 * * DKIS_ = 0.001997 0.000999 0.001997 0.000000 * * DROT_ = 0.004163 0.002083 0.004163 0.000003 * * * ******************************************************************************* ******************************************************************************* ******************************************************************************* * * * CONTRIBUTIONS TO ENERGY IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * ERHO_ = -4806.116087 1864.180234 -2942.520053 0.584200 * * ERHOD = 1556.564307 -778.526631 778.282153 -0.244478 * * ELPR_ = 176.432206 -37.591557 138.871443 -0.030794 * * ETAU_ = 426.509687 -174.238021 252.351565 -0.079899 * * ESCU_ = 0.000000 0.000000 0.000000 0.000000 * * EDIV_ = -23.100076 -11.558870 -34.650113 -0.008832 * * ============ ============ ============ ============ * * SUM EVEN: -2669.709962 862.265155 -1807.665005 0.220198 * * * * ESPI_ = -0.693172 0.770832 0.077019 0.000641 * * ESPID = 0.435886 -0.218122 0.217943 -0.000179 * * ELPS_ = 0.000000 -0.235269 -0.234938 -0.000331 * * ECUR_ = -0.598189 0.244403 -0.353928 0.000142 * * EKIS_ = 0.000000 0.000000 0.000000 0.000000 * * EROT_ = -0.249794 -0.124983 -0.374692 -0.000086 * * ============ ============ ============ ============ * * SUM ODD: -1.105269 0.436860 -0.668596 0.000187 * * * ******************************************************************************* ******************************************************************************* * * * CENTER-OF-MASS COORDINATES IN FERMIS TOTAL * * * * R_CMS = 0.0000 X_CMS = 0.0000 Y_CMS = 0.0000 Z_CMS = 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL * * * * ALPHA =270.00000 BETA = 0.00000 GAMMA = 0.00000 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 3.9215 X_RMS = 1.9736 Y_RMS = 2.2351 Z_RMS = 2.5469 * * * * R_GEO = 5.0626 X_GEO = 4.4132 Y_GEO = 4.9979 Z_GEO = 5.6951 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 64.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 2.6128 Q2+1= ZERO Q2+2= -1.2201 ............. ............. * * 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.0099 Q4+1= ZERO Q4+2= -0.0394 Q4+3= ZERO Q4+4= 0.0013 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * * Q50 = ZERO Q5+1= ZERO Q5+2= ZERO Q5+3= ZERO Q5+4= ZERO * * Q5-1= ZERO Q5-2= ZERO Q5-3= ZERO Q5-4= ZERO * * * * Q5+5= ZERO * * Q5-5= ZERO * * * * Q60 = -0.0032 Q6+1= ZERO Q6+2= 4.9E-04 Q6+3= ZERO Q6+4=-4.6E-04 * * Q6-1= ZERO Q6-2= ZERO Q6-3= ZERO Q6-4= ZERO * * * * Q6+5= ZERO Q6+6=-1.7E-05 * * Q6-5= ZERO Q6-6= ZERO * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (EXACT MULTIPOLE MOMENTS) TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = 0.2094 B2+1= ZERO B2+2= -0.1138 ............. ............. * * 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.0306 B4+1= ZERO B4+2= -0.0348 B4+3= ZERO B4+4= -0.0049 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * * B50 = ZERO B5+1= ZERO B5+2= ZERO B5+3= ZERO B5+4= ZERO * * B5-1= ZERO B5-2= ZERO B5-3= ZERO B5-4= ZERO * * * * B5+5= ZERO * * B5-5= ZERO * * * * B60 = -0.0054 B6+1= ZERO B6+2= 0.0147 B6+3= ZERO B6+4= -0.0098 * * B6-1= ZERO B6-2= ZERO B6-3= ZERO B6-4= ZERO * * * * B6+5= ZERO B6+6=-7.7E-04 * * B6-5= ZERO B6-6= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1219 FM^(-3) BOHR RADIUS = 4.9815 FM * * * ******************************************************************************* ******************************************************************************* * * * 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.2104 B2+1= ZERO B2+2= -0.0983 ............. ............. * * 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.0098 B4+1= ZERO B4+2= -0.0556 B4+3= ZERO B4+4= 0.0018 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * * B50 = ZERO B5+1= ZERO B5+2= ZERO B5+3= ZERO B5+4= ZERO * * B5-1= ZERO B5-2= ZERO B5-3= ZERO B5-4= ZERO * * * * B5+5= ZERO * * B5-5= ZERO * * * * B60 = -0.0123 B6+1= ZERO B6+2= 0.0027 B6+3= ZERO B6+4= -0.0025 * * B6-1= ZERO B6-2= ZERO B6-3= ZERO B6-4= ZERO * * * * B6+5= ZERO B6+6=-9.6E-05 * * B6-5= ZERO B6-6= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1178 FM^(-3) BOHR RADIUS = 5.0626 FM * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] [CENT-MASS FRAME] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 64.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 2.6128 Q2+1= ZERO Q2+2= -1.2201 ............. ............. * * 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.0099 Q4+1= ZERO Q4+2= -0.0394 Q4+3= ZERO Q4+4= 0.0013 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * * Q50 = ZERO Q5+1= ZERO Q5+2= ZERO Q5+3= ZERO Q5+4= ZERO * * Q5-1= ZERO Q5-2= ZERO Q5-3= ZERO Q5-4= ZERO * * * * Q5+5= ZERO * * Q5-5= ZERO * * * * Q60 = -0.0032 Q6+1= ZERO Q6+2= 4.9E-04 Q6+3= ZERO Q6+4=-4.6E-04 * * Q6-1= ZERO Q6-2= ZERO Q6-3= ZERO Q6-4= ZERO * * * * Q6+5= ZERO Q6+6=-1.7E-05 * * Q6-5= ZERO Q6-6= ZERO * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (EXACT MULTIPOLE MOMENTS) [CENT-MASS 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.2094 B2+1= ZERO B2+2= -0.1138 ............. ............. * * 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.0306 B4+1= ZERO B4+2= -0.0348 B4+3= ZERO B4+4= -0.0049 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * * B50 = ZERO B5+1= ZERO B5+2= ZERO B5+3= ZERO B5+4= ZERO * * B5-1= ZERO B5-2= ZERO B5-3= ZERO B5-4= ZERO * * * * B5+5= ZERO * * B5-5= ZERO * * * * B60 = -0.0054 B6+1= ZERO B6+2= 0.0147 B6+3= ZERO B6+4= -0.0098 * * B6-1= ZERO B6-2= ZERO B6-3= ZERO B6-4= ZERO * * * * B6+5= ZERO B6+6=-7.7E-04 * * B6-5= ZERO B6-6= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1219 FM^(-3) BOHR RADIUS = 4.9815 FM * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) [CENT-MASS 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.2104 B2+1= ZERO B2+2= -0.0983 ............. ............. * * 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.0098 B4+1= ZERO B4+2= -0.0556 B4+3= ZERO B4+4= 0.0018 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * * B50 = ZERO B5+1= ZERO B5+2= ZERO B5+3= ZERO B5+4= ZERO * * B5-1= ZERO B5-2= ZERO B5-3= ZERO B5-4= ZERO * * * * B5+5= ZERO * * B5-5= ZERO * * * * B60 = -0.0123 B6+1= ZERO B6+2= 0.0027 B6+3= ZERO B6+4= -0.0025 * * B6-1= ZERO B6-2= ZERO B6-3= ZERO B6-4= ZERO * * * * B6+5= ZERO B6+6=-9.6E-05 * * B6-5= ZERO B6-6= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1178 FM^(-3) BOHR RADIUS = 5.0626 FM * * * ******************************************************************************* ******************************************************************************* * * * LINEAR MOMENTA IN FERMI^(-1) TOTAL * * * * P_X = 0.0000 P_Y = 0.0000 P_Z = 0.0000 P_R = 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * 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.000 0.000 0.000 0.067 2.914 2.981 0.000 0.000 0.000 * * PROTONS 0.000 0.000 0.000 0.069 2.948 3.017 0.000 0.000 0.000 * * -------- * * TOTAL 0.000 0.000 0.000 0.136 5.862 5.998 0.000 0.000 0.000 * * * ******************************************************************************* * * * ANGLE BETWEEN ANGULAR MOMENTUM AND ANGULAR FREQUENCY VECTORS= 0.000 DEG * * * ******************************************************************************* * * * ANGULAR MOMENTA IN THE CENTER-OF-MASS 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.067 2.914 2.981 0.000 0.000 0.000 * * PROTONS 0.000 0.000 0.000 0.069 2.948 3.017 0.000 0.000 0.000 * * -------- * * TOTAL 0.000 0.000 0.000 0.136 5.862 5.998 0.000 0.000 0.000 * * * ******************************************************************************* ******************************************************************************* * * * NEUTRON CONFIGURATIONS * * ======================= * * 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 * * -------------------------------------------------------------------- * * * * CONF: 1 1 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 1 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * PROTON CONFIGURATIONS * * ======================= * * 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 * * -------------------------------------------------------------------- * * * * CONF: 1 1 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 1 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * ENERGIES (MEV) * * * ******************************************************************************* * * * KINETIC: (NEU)= 560.882688 (PRO)= 544.495493 (TOT)= 1105.378181 * * SUM EPS: (NEU)= -863.767069 (PRO)= -536.028071 (TOT)= -1399.795140 * * 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)= -12.538534 (PRO)= -2.732136 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 176.353367 (EXC)= -12.345178 (TOT)= 164.008189 * * * * CONSTR. (MULT)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 * * CONSTR. (SURF)= 0.000000 SLOPE= 0.000000 CORR.= 0.000000 * * CONSTR. (SPIN)= -3.448824 SLOPE= 0.575000 CORR.= -1.724412 * * * * REARRANGEMENT ENERGY FROM THE SKYRME DENSITY-DEPENDENT TERMS= 389.127720 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -542.175671 * * * * SPIN-ORB (EVE)= -34.658945 (ODD)= -0.374778 (TOT)= -35.033723 * * SKYRME: (EVE)= -1807.444807 (ODD)= -0.668409 (TOT)= -1808.113216 * * * * TOTAL: (STAB)= 2.41E-09 (SP)= -538.726847 (FUN)= -538.726847 * * * ******************************************************************************* ******************************************************************************* * * * THE GCM KERNELS IN FUNCTION OF THE EULER ANGLES ARE SAVED; FILE (FILKER): * * ge064-a.ker * * * ******************************************************************************* * * * N001-L001-R001 1 <= I <= 10 AND 1 <= K <= 10; TOT 100 PNTS ARE SAVD * * * ******************************************************************************* ******************************************************************************* * * * KERNELS AND AVERAGE VALUES OBTAINED FOR ANGULAR-MOMENTUM PROJECTED STATES * * * ******************************************************************************* * * * REAL(KERNEL) IMAG(KERNEL) AVERAGE ENERGY * * * * I,K= 0 0 NORM = 0.008035213920 0.000000000000 -542.137728 * * I,K= 1 -1 NORM = 0.000867253993 0.000000000000 -509.681810 * * I,K= 1 0 NORM = 0.000000000000 0.000000000000 -1424.000000 * * I,K= 1 1 NORM = 0.000867253993 0.000000000000 -509.681810 * * I,K= 2 -2 NORM = 0.014229647523 0.000000000000 -541.456596 * * I,K= 2 -1 NORM = 0.003145537179 0.000000000000 -530.120651 * * I,K= 2 0 NORM = 0.034483061433 0.000000000000 -541.966975 * * I,K= 2 1 NORM = 0.003145537179 0.000000000000 -530.120651 * * I,K= 2 2 NORM = 0.014229647523 0.000000000000 -541.456596 * * I,K= 3 -3 NORM = 0.001860645585 0.000000000000 -522.392590 * * I,K= 3 -2 NORM = 0.016828660109 0.000000000000 -541.247773 * * I,K= 3 -1 NORM = 0.002278617993 0.000000000000 -525.077250 * * I,K= 3 0 NORM = 0.000000000000 0.000000000000 -1726.769231 * * I,K= 3 1 NORM = 0.002278617993 0.000000000000 -525.077250 * * I,K= 3 2 NORM = 0.016828660109 0.000000000000 -541.247773 * * I,K= 3 3 NORM = 0.001860645585 0.000000000000 -522.392590 * * I,K= 4 -4 NORM = 0.011187959146 0.000000000000 -538.708768 * * I,K= 4 -3 NORM = 0.003162059172 0.000000000000 -531.632290 * * I,K= 4 -2 NORM = 0.018931571592 0.000000000000 -541.275028 * * I,K= 4 -1 NORM = 0.009960784427 0.000000000000 -539.256396 * * I,K= 4 0 NORM = 0.044802398411 0.000000000000 -541.701436 * * I,K= 4 1 NORM = 0.009960784427 0.000000000000 -539.256396 * * I,K= 4 2 NORM = 0.018931571592 0.000000000000 -541.275028 * * I,K= 4 3 NORM = 0.003162059172 0.000000000000 -531.632290 * * I,K= 4 4 NORM = 0.011187959146 0.000000000000 -538.708768 * * I,K= 5 -5 NORM = 0.002068177721 0.000000000000 -525.538628 * * I,K= 5 -4 NORM = 0.010587976562 0.000000000000 -538.450837 * * I,K= 5 -3 NORM = 0.004288497776 0.000000000000 -535.933341 * * I,K= 5 -2 NORM = 0.016252757944 0.000000000000 -540.734127 * * I,K= 5 -1 NORM = 0.003102066068 0.000000000000 -534.723124 * * I,K= 5 0 NORM = 0.000000000000 0.000000000000 -1920.000000 * * I,K= 5 1 NORM = 0.003102066068 0.000000000000 -534.723124 * * I,K= 5 2 NORM = 0.016252757944 0.000000000000 -540.734127 * * I,K= 5 3 NORM = 0.004288497776 0.000000000000 -535.933341 * * I,K= 5 4 NORM = 0.010587976562 0.000000000000 -538.450837 * * I,K= 5 5 NORM = 0.002068177721 0.000000000000 -525.538628 * * * ******************************************************************************* * * * SUM RULES OBTAINED FOR ANGULAR-MOMENTUM PROJECTED STATES VS THE HF VALUES * * * ******************************************************************************* * * * REAL(SUM RULE) IMAG(SUM RULE) HF VALUE * * * * NORM= 0.324825099346 0.000000000000 * * SKYRME= -587.075787 0.000000 -1808.113216 * * EKIN_T= 358.576805 0.000000 1105.378181 * * COUL_D= 57.290175 0.000000 176.353367 * * COUL_E= -4.008438 0.000000 -12.345178 * * * * MULTIPOLE L= 0 M= 0 QMUL_P= 10.394403 0.000000 32.000000 * * L= 1 M=-1 QMUL_P= 0.000000 0.000000 0.000000 * * L= 1 M= 0 QMUL_P= 0.000000 0.000000 0.000000 * * L= 1 M= 1 QMUL_P= 0.000000 0.000000 0.000000 * * L= 2 M=-2 QMUL_P= -0.206203 0.000000 0.000000 * * L= 2 M=-1 QMUL_P= 0.000000 0.032421 0.000000 * * L= 2 M= 0 QMUL_P= 0.432164 0.000000 1.333827 * * L= 2 M= 1 QMUL_P= 0.000000 -0.032421 0.000000 * * L= 2 M= 2 QMUL_P= -0.206203 0.000000 -0.621629 * * * ******************************************************************************* * * * REDUCED KERNELS OF PROTON MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] * * ONLY THE RESULTS WITH ABSOLUTE VALUES LARGER THAN 0.001 ARE BELOW PRINTED * * * ******************************************************************************* * * * * IL KL L IR KR <||Q||> * IL KL L IR KR <||Q||> * * * * * 0 0 2 2 -2 -0.008 * 2 -2 2 2 -2 0.025 * * 2 0 2 2 -2 -0.017 * 2 2 2 2 -2 0.025 * * 1 -1 2 2 -1 0.003 * 1 1 2 2 -1 0.003 * * 0 0 2 2 0 0.024 * 2 -2 2 2 0 -0.017 * * 2 0 2 2 0 -0.055 * 2 2 2 2 0 -0.017 * * 1 -1 2 2 1 -0.003 * 1 1 2 2 1 -0.003 * * 0 0 2 2 2 -0.008 * 2 -2 2 2 2 0.025 * * 2 0 2 2 2 -0.017 * 2 2 2 2 2 0.025 * * 1 -1 2 3 -3 -0.001 * 1 1 2 3 -3 -0.001 * * 3 -3 2 3 -3 0.004 * 3 3 2 3 -3 0.004 * * 2 -2 2 3 -2 0.034 * 2 0 2 3 -2 0.022 * * 2 2 2 3 -2 0.034 * 1 -1 2 3 -1 0.002 * * 1 1 2 3 -1 0.002 * 2 -1 2 3 -1 0.002 * * 2 1 2 3 -1 -0.002 * 3 -1 2 3 -1 -0.003 * * 3 1 2 3 -1 -0.003 * 1 -1 2 3 1 0.002 * * 1 1 2 3 1 0.002 * 2 -1 2 3 1 0.002 * * 2 1 2 3 1 -0.002 * 3 -1 2 3 1 -0.003 * * 3 1 2 3 1 -0.003 * 2 -2 2 3 2 -0.034 * * 2 0 2 3 2 -0.022 * 2 2 2 3 2 -0.034 * * 1 -1 2 3 3 -0.001 * 1 1 2 3 3 -0.001 * * 3 -3 2 3 3 0.004 * 3 3 2 3 3 0.004 * * 2 -2 2 4 -4 -0.015 * 2 0 2 4 -4 0.003 * * 2 2 2 4 -4 -0.015 * 3 -2 2 4 -4 0.008 * * 3 2 2 4 -4 -0.008 * 4 -4 2 4 -4 0.034 * * 4 -2 2 4 -4 -0.009 * 4 0 2 4 -4 0.004 * * 4 2 2 4 -4 -0.009 * 4 4 2 4 -4 0.034 * * 2 -1 2 4 -3 -0.004 * 2 1 2 4 -3 0.004 * * 3 -3 2 4 -3 0.005 * 3 -1 2 4 -3 -0.002 * * 3 1 2 4 -3 -0.002 * 3 3 2 4 -3 0.005 * * 4 -3 2 4 -3 0.003 * 4 -1 2 4 -3 -0.003 * * 4 1 2 4 -3 0.003 * 4 3 2 4 -3 -0.003 * * 2 -2 2 4 -2 0.026 * 2 0 2 4 -2 -0.030 * * 2 2 2 4 -2 0.026 * 3 -2 2 4 -2 0.040 * * 3 2 2 4 -2 -0.040 * 4 -4 2 4 -2 -0.009 * * 4 -2 2 4 -2 -0.007 * 4 0 2 4 -2 -0.018 * * 4 2 2 4 -2 -0.007 * 4 4 2 4 -2 -0.009 * * 2 -1 2 4 -1 0.013 * 2 1 2 4 -1 -0.013 * * 3 -3 2 4 -1 -0.005 * 3 -1 2 4 -1 0.008 * * 3 1 2 4 -1 0.008 * 3 3 2 4 -1 -0.005 * * 4 -3 2 4 -1 -0.003 * 4 -1 2 4 -1 -0.007 * * 4 1 2 4 -1 0.007 * 4 3 2 4 -1 0.003 * * 2 -2 2 4 0 -0.014 * 2 0 2 4 0 0.093 * * 2 2 2 4 0 -0.014 * 3 -2 2 4 0 -0.028 * * 3 2 2 4 0 0.028 * 4 -4 2 4 0 0.004 * * 4 -2 2 4 0 -0.018 * 4 0 2 4 0 -0.078 * * 4 2 2 4 0 -0.018 * 4 4 2 4 0 0.004 * * 2 -1 2 4 1 -0.013 * 2 1 2 4 1 0.013 * * 3 -3 2 4 1 0.005 * 3 -1 2 4 1 -0.008 * * 3 1 2 4 1 -0.008 * 3 3 2 4 1 0.005 * * 4 -3 2 4 1 0.003 * 4 -1 2 4 1 0.007 * * 4 1 2 4 1 -0.007 * 4 3 2 4 1 -0.003 * * 2 -2 2 4 2 0.026 * 2 0 2 4 2 -0.030 * * 2 2 2 4 2 0.026 * 3 -2 2 4 2 0.040 * * 3 2 2 4 2 -0.040 * 4 -4 2 4 2 -0.009 * * 4 -2 2 4 2 -0.007 * 4 0 2 4 2 -0.018 * * 4 2 2 4 2 -0.007 * 4 4 2 4 2 -0.009 * * 2 -1 2 4 3 0.004 * 2 1 2 4 3 -0.004 * * 3 -3 2 4 3 -0.005 * 3 -1 2 4 3 0.002 * * 3 1 2 4 3 0.002 * 3 3 2 4 3 -0.005 * * 4 -3 2 4 3 -0.003 * 4 -1 2 4 3 0.003 * * 4 1 2 4 3 -0.003 * 4 3 2 4 3 0.003 * * 2 -2 2 4 4 -0.015 * 2 0 2 4 4 0.003 * * 2 2 2 4 4 -0.015 * 3 -2 2 4 4 0.008 * * 3 2 2 4 4 -0.008 * 4 -4 2 4 4 0.034 * * 4 -2 2 4 4 -0.009 * 4 0 2 4 4 0.004 * * 4 2 2 4 4 -0.009 * 4 4 2 4 4 0.034 * * 3 -3 2 5 -5 -0.003 * 3 -1 2 5 -5 0.001 * * 3 1 2 5 -5 0.001 * 3 3 2 5 -5 -0.003 * * 4 -1 2 5 -5 -0.001 * 4 1 2 5 -5 0.001 * * 5 -5 2 5 -5 0.007 * 5 -3 2 5 -5 -0.004 * * 5 -1 2 5 -5 0.001 * 5 1 2 5 -5 0.001 * * 5 3 2 5 -5 -0.004 * 5 5 2 5 -5 0.007 * * 3 -2 2 5 -4 -0.017 * 3 2 2 5 -4 0.017 * * 4 -4 2 5 -4 0.028 * 4 -2 2 5 -4 0.007 * * 4 0 2 5 -4 -0.006 * 4 2 2 5 -4 0.007 * * 4 4 2 5 -4 0.028 * 5 -4 2 5 -4 0.017 * * 5 -2 2 5 -4 -0.008 * 5 2 2 5 -4 0.008 * * 5 4 2 5 -4 -0.017 * 3 -3 2 5 -3 0.005 * * 3 -1 2 5 -3 -0.005 * 3 1 2 5 -3 -0.005 * * 3 3 2 5 -3 0.005 * 4 -3 2 5 -3 0.007 * * 4 -1 2 5 -3 0.002 * 4 1 2 5 -3 -0.002 * * 4 3 2 5 -3 -0.007 * 5 -5 2 5 -3 -0.004 * * 5 -3 2 5 -3 0.001 * 5 -1 2 5 -3 0.001 * * 5 1 2 5 -3 0.001 * 5 3 2 5 -3 0.001 * * 5 5 2 5 -3 -0.004 * 3 -2 2 5 -2 0.036 * * 3 2 2 5 -2 -0.036 * 4 -4 2 5 -2 -0.009 * * 4 -2 2 5 -2 0.030 * 4 0 2 5 -2 0.023 * * 4 2 2 5 -2 0.030 * 4 4 2 5 -2 -0.009 * * 5 -4 2 5 -2 -0.008 * 5 -2 2 5 -2 -0.020 * * 5 2 2 5 -2 0.020 * 5 4 2 5 -2 0.008 * * 3 -3 2 5 -1 -0.003 * 3 -1 2 5 -1 0.006 * * 3 1 2 5 -1 0.006 * 3 3 2 5 -1 -0.003 * * 4 -3 2 5 -1 -0.005 * 4 -1 2 5 -1 0.001 * * 4 1 2 5 -1 -0.001 * 4 3 2 5 -1 0.005 * * 5 -5 2 5 -1 0.001 * 5 -3 2 5 -1 0.001 * * 5 -1 2 5 -1 -0.006 * 5 1 2 5 -1 -0.006 * * 5 3 2 5 -1 0.001 * 5 5 2 5 -1 0.001 * * 3 -3 2 5 1 -0.003 * 3 -1 2 5 1 0.006 * * 3 1 2 5 1 0.006 * 3 3 2 5 1 -0.003 * * 4 -3 2 5 1 -0.005 * 4 -1 2 5 1 0.001 * * 4 1 2 5 1 -0.001 * 4 3 2 5 1 0.005 * * 5 -5 2 5 1 0.001 * 5 -3 2 5 1 0.001 * * 5 -1 2 5 1 -0.006 * 5 1 2 5 1 -0.006 * * 5 3 2 5 1 0.001 * 5 5 2 5 1 0.001 * * 3 -2 2 5 2 -0.036 * 3 2 2 5 2 0.036 * * 4 -4 2 5 2 0.009 * 4 -2 2 5 2 -0.030 * * 4 0 2 5 2 -0.023 * 4 2 2 5 2 -0.030 * * 4 4 2 5 2 0.009 * 5 -4 2 5 2 0.008 * * 5 -2 2 5 2 0.020 * 5 2 2 5 2 -0.020 * * 5 4 2 5 2 -0.008 * 3 -3 2 5 3 0.005 * * 3 -1 2 5 3 -0.005 * 3 1 2 5 3 -0.005 * * 3 3 2 5 3 0.005 * 4 -3 2 5 3 0.007 * * 4 -1 2 5 3 0.002 * 4 1 2 5 3 -0.002 * * 4 3 2 5 3 -0.007 * 5 -5 2 5 3 -0.004 * * 5 -3 2 5 3 0.001 * 5 -1 2 5 3 0.001 * * 5 1 2 5 3 0.001 * 5 3 2 5 3 0.001 * * 5 5 2 5 3 -0.004 * 3 -2 2 5 4 0.017 * * 3 2 2 5 4 -0.017 * 4 -4 2 5 4 -0.028 * * 4 -2 2 5 4 -0.007 * 4 0 2 5 4 0.006 * * 4 2 2 5 4 -0.007 * 4 4 2 5 4 -0.028 * * 5 -4 2 5 4 -0.017 * 5 -2 2 5 4 0.008 * * 5 2 2 5 4 -0.008 * 5 4 2 5 4 0.017 * * 3 -3 2 5 5 -0.003 * 3 -1 2 5 5 0.001 * * 3 1 2 5 5 0.001 * 3 3 2 5 5 -0.003 * * 4 -1 2 5 5 -0.001 * 4 1 2 5 5 0.001 * * 5 -5 2 5 5 0.007 * 5 -3 2 5 5 -0.004 * * 5 -1 2 5 5 0.001 * 5 1 2 5 5 0.001 * * 5 3 2 5 5 -0.004 * 5 5 2 5 5 0.007 * * * ******************************************************************************* * * * REDUCED KERNELS OF PROTON MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] * * ONLY THE RESULTS WITH ABSOLUTE VALUES LARGER THAN 0.001 ARE BELOW PRINTED * * * ******************************************************************************* * * * * IL KL L IR KR <||M||> * IL KL L IR KR <||M||> * * * * * 2 -2 1 2 -2 0.019 * 2 0 1 2 -2 -0.003 * * 2 2 1 2 -2 0.019 * 2 -1 1 2 -1 0.006 * * 2 1 1 2 -1 -0.006 * 2 -2 1 2 0 -0.003 * * 2 0 1 2 0 0.052 * 2 2 1 2 0 -0.003 * * 2 -1 1 2 1 -0.006 * 2 1 1 2 1 0.006 * * 2 -2 1 2 2 0.019 * 2 0 1 2 2 -0.003 * * 2 2 1 2 2 0.019 * 3 -3 1 3 -3 0.004 * * 3 -1 1 3 -3 -0.002 * 3 1 1 3 -3 -0.002 * * 3 3 1 3 -3 0.004 * 2 -2 1 3 -2 -0.001 * * 2 2 1 3 -2 -0.001 * 3 -2 1 3 -2 0.041 * * 3 2 1 3 -2 -0.041 * 3 -3 1 3 -1 -0.002 * * 3 -1 1 3 -1 0.007 * 3 1 1 3 -1 0.007 * * 3 3 1 3 -1 -0.002 * 3 -3 1 3 1 -0.002 * * 3 -1 1 3 1 0.007 * 3 1 1 3 1 0.007 * * 3 3 1 3 1 -0.002 * 2 -2 1 3 2 0.001 * * 2 2 1 3 2 0.001 * 3 -2 1 3 2 -0.041 * * 3 2 1 3 2 0.041 * 3 -3 1 3 3 0.004 * * 3 -1 1 3 3 -0.002 * 3 1 1 3 3 -0.002 * * 3 3 1 3 3 0.004 * 4 -4 1 4 -4 0.026 * * 4 -2 1 4 -4 -0.006 * 4 2 1 4 -4 -0.006 * * 4 4 1 4 -4 0.026 * 3 -3 1 4 -3 -0.001 * * 3 -1 1 4 -3 0.001 * 3 1 1 4 -3 0.001 * * 3 3 1 4 -3 -0.001 * 4 -3 1 4 -3 0.012 * * 4 -1 1 4 -3 -0.010 * 4 1 1 4 -3 0.010 * * 4 3 1 4 -3 -0.012 * 3 -2 1 4 -2 -0.003 * * 3 2 1 4 -2 0.003 * 4 -4 1 4 -2 -0.006 * * 4 -2 1 4 -2 0.069 * 4 0 1 4 -2 -0.042 * * 4 2 1 4 -2 0.069 * 4 4 1 4 -2 -0.006 * * 4 -3 1 4 -1 -0.010 * 4 -1 1 4 -1 0.044 * * 4 1 1 4 -1 -0.044 * 4 3 1 4 -1 0.010 * * 4 -2 1 4 0 -0.042 * 4 0 1 4 0 0.173 * * 4 2 1 4 0 -0.042 * 4 -3 1 4 1 0.010 * * 4 -1 1 4 1 -0.044 * 4 1 1 4 1 0.044 * * 4 3 1 4 1 -0.010 * 3 -2 1 4 2 -0.003 * * 3 2 1 4 2 0.003 * 4 -4 1 4 2 -0.006 * * 4 -2 1 4 2 0.069 * 4 0 1 4 2 -0.042 * * 4 2 1 4 2 0.069 * 4 4 1 4 2 -0.006 * * 3 -3 1 4 3 0.001 * 3 -1 1 4 3 -0.001 * * 3 1 1 4 3 -0.001 * 3 3 1 4 3 0.001 * * 4 -3 1 4 3 -0.012 * 4 -1 1 4 3 0.010 * * 4 1 1 4 3 -0.010 * 4 3 1 4 3 0.012 * * 4 -4 1 4 4 0.026 * 4 -2 1 4 4 -0.006 * * 4 2 1 4 4 -0.006 * 4 4 1 4 4 0.026 * * 4 -1 1 5 -5 0.003 * 4 1 1 5 -5 -0.003 * * 5 -5 1 5 -5 0.002 * 5 -3 1 5 -5 -0.003 * * 5 3 1 5 -5 -0.003 * 5 5 1 5 -5 0.002 * * 4 -4 1 5 -4 -0.008 * 4 -2 1 5 -4 0.001 * * 4 0 1 5 -4 0.003 * 4 2 1 5 -4 0.001 * * 4 4 1 5 -4 -0.008 * 5 -4 1 5 -4 0.042 * * 5 -2 1 5 -4 -0.018 * 5 2 1 5 -4 0.018 * * 5 4 1 5 -4 -0.042 * 4 -3 1 5 -3 -0.002 * * 4 -1 1 5 -3 0.002 * 4 1 1 5 -3 -0.002 * * 4 3 1 5 -3 0.002 * 5 -5 1 5 -3 -0.003 * * 5 -3 1 5 -3 0.025 * 5 -1 1 5 -3 -0.017 * * 5 1 1 5 -3 -0.017 * 5 3 1 5 -3 0.025 * * 5 5 1 5 -3 -0.003 * 4 -2 1 5 -2 -0.003 * * 4 0 1 5 -2 0.006 * 4 2 1 5 -2 -0.003 * * 5 -4 1 5 -2 -0.018 * 5 -2 1 5 -2 0.088 * * 5 2 1 5 -2 -0.088 * 5 4 1 5 -2 0.018 * * 4 -3 1 5 -1 -0.001 * 4 -1 1 5 -1 -0.004 * * 4 1 1 5 -1 0.004 * 4 3 1 5 -1 0.001 * * 5 -3 1 5 -1 -0.017 * 5 -1 1 5 -1 0.022 * * 5 1 1 5 -1 0.022 * 5 3 1 5 -1 -0.017 * * 4 -3 1 5 1 -0.001 * 4 -1 1 5 1 -0.004 * * 4 1 1 5 1 0.004 * 4 3 1 5 1 0.001 * * 5 -3 1 5 1 -0.017 * 5 -1 1 5 1 0.022 * * 5 1 1 5 1 0.022 * 5 3 1 5 1 -0.017 * * 4 -2 1 5 2 0.003 * 4 0 1 5 2 -0.006 * * 4 2 1 5 2 0.003 * 5 -4 1 5 2 0.018 * * 5 -2 1 5 2 -0.088 * 5 2 1 5 2 0.088 * * 5 4 1 5 2 -0.018 * 4 -3 1 5 3 -0.002 * * 4 -1 1 5 3 0.002 * 4 1 1 5 3 -0.002 * * 4 3 1 5 3 0.002 * 5 -5 1 5 3 -0.003 * * 5 -3 1 5 3 0.025 * 5 -1 1 5 3 -0.017 * * 5 1 1 5 3 -0.017 * 5 3 1 5 3 0.025 * * 5 5 1 5 3 -0.003 * 4 -4 1 5 4 0.008 * * 4 -2 1 5 4 -0.001 * 4 0 1 5 4 -0.003 * * 4 2 1 5 4 -0.001 * 4 4 1 5 4 0.008 * * 5 -4 1 5 4 -0.042 * 5 -2 1 5 4 0.018 * * 5 2 1 5 4 -0.018 * 5 4 1 5 4 0.042 * * 4 -1 1 5 5 0.003 * 4 1 1 5 5 -0.003 * * 5 -5 1 5 5 0.002 * 5 -3 1 5 5 -0.003 * * 5 3 1 5 5 -0.003 * 5 5 1 5 5 0.002 * * * ******************************************************************************* * * * RESULTS OF THE K-MIXING CALCULATION FOR THE FOLLOWING OVERLAP CUTOFF DATA * * ICUTOV=1, CUTOVE= 1.0E-02, CUTOVF= 0.0E+00 * * * ******************************************************************************* * | * * I N OVERLAP ENERGY | I N OVERLAP ENERGY * * | * *--------------------------------------|--------------------------------------* * 2 1 4.02E-02 -541.817486 | 3 1 3.78E-02 -540.997226 * * 2 2 2.71E-02 -541.344809 | * *--------------------------------------|--------------------------------------* * 4 1 7.90E-02 -541.775089 | 5 1 5.06E-02 -540.367910 * * 4 2 3.04E-02 -540.196720 | 5 2 2.04E-02 -537.665199 * * 4 3 2.06E-02 -538.486591 | * *--------------------------------------|--------------------------------------* ******************************************************************************* * * * REDUCED MAT.ELEMS OF PROTON MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] * * ONLY THE RESULTS WITH ABSOLUTE VALUES LARGER THAN 0.001 ARE BELOW PRINTED * * * ******************************************************************************* * * * * IL NL L IR NR <||Q||> * IL NL L IR NR <||Q||> * * * * * 2 1 2 2 1 -1.142 * 2 2 2 2 1 -1.419 * * 2 1 2 2 2 -1.419 * 2 2 2 2 2 1.142 * * 3 1 2 3 1 -0.007 * 2 1 2 4 1 -2.353 * * 2 2 2 4 1 0.543 * 4 1 2 4 1 -0.518 * * 4 2 2 4 1 1.023 * 4 3 2 4 1 -0.417 * * 2 1 2 4 2 0.554 * 2 2 2 4 2 1.264 * * 4 1 2 4 2 1.023 * 4 2 2 4 2 -2.302 * * 4 3 2 4 2 -0.763 * 2 1 2 4 3 0.159 * * 2 2 2 4 3 -0.831 * 4 1 2 4 3 -0.417 * * 4 2 2 4 3 -0.763 * 4 3 2 4 3 2.832 * * 3 1 2 5 1 -2.170 * 5 1 2 5 1 -0.937 * * 5 2 2 5 1 1.296 * 3 1 2 5 2 -0.381 * * 5 1 2 5 2 1.296 * 5 2 2 5 2 0.982 * * * ******************************************************************************* * * * REDUCED MAT.ELEMS OF PROTON MAGNETIC MOMENTS [MAGNETON*FERMI^(LAMBDA-1)] * * ONLY THE RESULTS WITH ABSOLUTE VALUES LARGER THAN 0.001 ARE BELOW PRINTED * * * ******************************************************************************* * * * * IL NL L IR NR <||M||> * IL NL L IR NR <||M||> * * * * * 2 1 1 2 1 1.509 * 2 2 1 2 1 0.015 * * 2 1 1 2 2 0.015 * 2 2 1 2 2 1.382 * * 3 1 1 3 1 2.433 * 4 1 1 4 1 3.642 * * 4 2 1 4 1 -0.026 * 4 3 1 4 1 0.070 * * 4 1 1 4 2 -0.026 * 4 2 1 4 2 3.579 * * 4 3 1 4 2 0.057 * 4 1 1 4 3 0.070 * * 4 2 1 4 3 0.057 * 4 3 1 4 3 2.375 * * 5 1 1 5 1 4.843 * 5 2 1 5 1 -0.168 * * 5 1 1 5 2 -0.168 * 5 2 1 5 2 3.700 * * * ******************************************************************************* ******************************************************************************* * * * WARNING! THE STABILITY REQUIRED WAS NOT ACHIEVED. NO-OF-ITERATIONS EXCEEDED * * WARNING! THE STABILITY REQUIRED WAS NOT ACHIEVED. NO-OF-ITERATIONS EXCEEDED * * * ******************************************************************************* * * ******************************************************************************* * * * CALCULATION OF THE DATA SET NO. 4 HAS BEEN SUCCESSFULLY COMPLETED * * * ******************************************************************************* ******************************************************************************* * * * PROCESSOR EXECUTION TIMES IN SUBROUTINES * * INCLUDING THEIR COMPLETE DOWN-CALLING TREES * * * ******************************************************************************* * * * 15137 => HFODD 14338 => PROANG 8151 => DENSHF 2994 => DENMAC * * * * 1263 => MOMETS 1010 => COUMAT 995 => INTCOU 982 => INTMUL * * * * 468 => MAGMOM 424 => ROTWAV 194 => INTSPI 185 => INTEGH * * * * 147 => INTORB 106 => DIAMAT 89 => DIASIG 58 => INTSOR * * * * 53 => SPAVER 37 => AVOBSE 23 => INTKIN 23 => DOBROY * * * * 21 => INTMAS 19 => DIASIZ 18 => AVANGY 16 => BEGINT * * * * 16 => SKFILD 15 => POWALL 14 => BEGINC 14 => INTCEN * * * * 7 => NILABS 6 => RECORD 3 => LINAVR 2 => SPIMOM * * * * 0 => NILSON 0 => AVPARI 0 => RECOUL 0 => MOMSIF * * * * 0 => AVSIMP 0 => MOMVMU 0 => ROTMOM 0 => SHISIF * * * * 0 => SHICOE 0 => SHIMAG 0 => REDPRO 0 => SHIMOM * * * * 0 => DIAPRO * * * ******************************************************************************* * * * NUMBERS OF CALLS TO SUBROUTINES * * * ******************************************************************************* * * * 26766 => INTMUL 4002 => INTSPI 4001 => INTORB 2452 => BEGINT * * * * 2332 => INTKIN 2166 => DENSHF 2166 => DENMAC 2000 => ROTWAV * * * * 1166 => MOMETS 1166 => MOMSIF 1166 => MAGMOM 1086 => COUMAT * * * * 1086 => BEGINC 1086 => INTCOU 564 => INTSOR 563 => DIAMAT * * * * 380 => INTMAS 380 => INTCEN 255 => ROTMOM 166 => AVPARI * * * * 166 => MOMVMU 166 => SPIMOM 166 => LINAVR 164 => DIASIG * * * * 164 => INTEGH 106 => SPAVER 85 => RECORD 83 => SKFILD * * * * 83 => SHICOE 83 => SHIMOM 83 => SHISIF 83 => SHIMAG * * * * 61 => DOBROY 14 => AVANGY 14 => AVOBSE 14 => AVSIMP * * * * 14 => NILABS 5 => RECOUL 2 => POWALL 2 => NILSON * * * * 2 => DIASIZ 1 => HFODD 1 => PROANG 1 => DIAPRO * * * * 1 => REDPRO * * * ******************************************************************************* * * * EXECUTION ENDS ON 2009.03.05 AT 22:31:21.642 * * * *******************************************************************************