*******************************************************************************
* *
* 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 *
* *
*******************************************************************************