******************************************************************************* * * * S I N G L E - C O R E V E R S I O N * * * ******************************************************************************* * * * HFODD HFODD HFODD HFODD HFODD HFODD HFODD HFODD * * * ******************************************************************************* * * * SKYRME-HARTREE-FOCK-BOGOLYUBOV CODE VERSION: 2.49T * * * * NO SYMMETRY-PLANES AND NO TIME-REVERSAL SYMMETRY * * * * DEFORMED CARTESIAN HARMONIC-OSCILLATOR BASIS * * * ******************************************************************************* * * * J. DOBACZEWSKI, B.G. CARLSSON, J. DUDEK, J. ENGEL * * J. MCDONNELL, P. OLBRATOWSKI, P. POWALOWSKI, M. SADZIAK * * J. SARICH, W. SATULA, N. SCHUNCK, J.A. SHEIKH * * A. STASZCZAK, M. STOITSOV, P. TOIVANEN, M. ZALEWSKI * * AND H. ZDUNCZUK * * * * INSTYTUT FIZYKI TEORETYCZNEJ, WARSZAWA * * LAWRENCE LIVERMORE NATIONAL LABORATORY, USA * * * * 1993-2011 * * * ******************************************************************************* ******************************************************************************* * * * CODE COMPILED WITH THE FOLLOWING ARRAY DIMENSIONS AND SWITCHES: * * * ******************************************************************************* * * * NDBASE = 680 NDSTAT = 680 NDXHRM = 40 NDYHRM = 40 NDZHRM = 40 * * * * NDMAIN = 16 NDMULT = 9 NDMULR = 4 NDLAMB = 9 NDITER = 5000 * * * * NDAKNO = 1 NDBKNO = 1 NDPROI = 20 NDCOUL = 80 NDPOLS = 25 * * * * NDPROT = 10 NDBTKN = 10 * * * * IPARAL = 0 I_CRAY = 0 * * * * * * PRE-PROCESSOR OPTIONS: * * * * switch_port = 1 switch_diag = 3 switch_cray = 0 * * * * switch_nagl = 0 switch_quad = 0 switch_vect = 1 * * * ******************************************************************************* * * * EXECUTION BEGINS ON 2011.02.24 AT 01:03:12.820 * * * ******************************************************************************* ******************************************************************************* * * * CLASSICAL NUCLEAR SURFACE DEFINED FOR: N = 28 Z = 20 * * * ******************************************************************************* * * * SPHERICAL SHAPE REQUESTED FOR THE SURFACE * * * ******************************************************************************* * * * HOMEGA= 13.5379 FCHOM0= 1.2000 * * * * OSCILLATOR FREQUENCIES: HBAROX= 13.5379 HBAROY= 13.5379 HBAROZ= 13.5379 * * * * MOMENTS OF INERTIA: XMOMFC= 9.2482 YMOMFC= 9.2482 ZMOMFC= 9.2482 * * * * CENTERS OF MASS: CMSXFC= 0.0000 CMSYFC= 0.0000 CMSZFC= 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * SKYRME FORCE DEFINITION * * * ******************************************************************************* * * * PARAMETER SET SLY4: T0= -2488.91 T1= 486.82 T2= -546.39 T3= 13777.00 * * * * POWER=1/6 W=123.000 X0= 0.83400 X1= -0.34400 X2= -1.00000 X3= 1.35400 * * * * ETA_J=0 ETA_W=0 ETACM=0 ETA_M=1 HBM= 20.73553 (SKYRME-FORCE SPIN-ORBIT) * * * ******************************************************************************* ******************************************************************************* * * * COUPLING CONSTANTS DEFINING THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * CRHO_ = -1763.394861 1660.104971 -933.342375 830.052485 * * CRHOD = 1925.335750 -2128.546500 861.062500 -1064.273250 * * CLPR_ = -92.653338 31.314270 -76.996203 15.657135 * * CTAU_ = 32.471951 49.313473 57.128687 24.656736 * * CSCU_ = 0.000000 0.000000 0.000000 0.000000 * * CDIV_ = -61.500000 -61.500000 -92.250000 -30.750000 * * * * CSPI_ = -518.938360 622.228250 -207.824235 311.114125 * * CSPID = 777.252417 -574.041667 490.231583 -287.020833 * * CLPS_ = 32.774724 28.564344 47.056896 14.282172 * * CCUR_ = -32.471951 -49.313473 -57.128687 -24.656736 * * CKIS_ = 0.000000 0.000000 0.000000 0.000000 * * CROT_ = -61.500000 -61.500000 -92.250000 -30.750000 * * * ******************************************************************************* ******************************************************************************* * * * PHYSICAL CONSTANTS: H_BARC=197.32891000 HBCOE2=137.03602000 * * * * XMASSN=938.90590000 XMASSP=938.27231000 * * * * HBMASS= 20.73553000 HBMRPA= 20.30353979 * * * * ECHAR2= 1.43997841 COULEX= -1.06350868 * * * ******************************************************************************* * * * OSCILLATOR LENGTHS: X= 1.7502658 Y= 1.7502658 Z= 1.7502658 * * * * OSCILLATOR CONSTANTS: X= 0.5713418 Y= 0.5713418 Z= 0.5713418 * * * * OSCILLATOR FREQUENCIES: X=13.5379017 Y=13.5379017 Z=13.5379017 * * * ******************************************************************************* * * * 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= 209.8375 * * * ******************************************************************************* ******************************************************************************* * * * 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): * * ca48_lip.cou * * * ******************************************************************************* * * * NUCLIDE: N = 28 Z = 20 NUMBER OF ITERATIONS = 200 CONTINUATION? = 0 * * * ******************************************************************************* * * * ITERATIONS WILL STOP WHEN THE STABILITY BECOMES LOWER THAN EPSITE= 1.0E-07 * * FOR 5 CONSECUTIVE ITERATIONS * * * ******************************************************************************* * * * ITERATIONS WILL STOP WHEN THE PING-PONG DIVERGENCE OCCURS FOR EPSPNG=0.010 * * DURING 3 PAIRS OF CONSECUTIVE ITERATIONS * * * ******************************************************************************* * * * ITERATIONS WILL STOP WHEN THE CHAOTIC-DIVERGENCE CONDITION OCCURS 5 TIMES * * * ******************************************************************************* * * * SLOW-DOWN PARAMETERS: * * TIME-EVEN = 0.50 TIME-ODD = 0.50 PAIRING = 0.50 LIPKIN-NOGAMI = 0.50 * * * ******************************************************************************* * * * BROYDEN METHOD IS: OFF * * * * ITERATIONS WILL PROCEED BY STANDARD LINEAR MIXING WITH SLOWING FACTORS * * DISPLAYED EARLIER * * * ******************************************************************************* * * * 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=3 * * FOR THE COULOMB FIELD, NMUCOU=4 * * FOR THE OUTPUT INFO, NMUPRI=4 * * * ******************************************************************************* * * * PRINTING THE RESULTS FOR THE FOLLOWING ITERATIONS: THE FIRST: YES * * THE MIDDLE: NO * * AND/OR THE LAST: YES * * * ******************************************************************************* * * * PRINTING THE MOMENTS, RADII, AND DEFORMATIONS FOR: NEUTRONS: YES * * PROTONS : YES * * MASS : YES * * * ******************************************************************************* * * * PRINTING THE MOMENTS, MOMENTA AND DEFORMATIONS IN THE INTRINSIC FRAME: YES * * * ******************************************************************************* * * * 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: 28 28 * * PROTONS: 20 20 * * * ******************************************************************************* * * * 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 7 7 0 0 0 0 0 0 0 0 * * PROTONS : 7 7 3 3 0 0 0 0 0 0 0 0 * * * ******************************************************************************* ******************************************************************************* * * * PARAMETERS OF THE NILSSON MODEL NEUTRONS * * * ******************************************************************************* * * * ESCALE =68.921 KAPPA = 0.080 MIU = 0.000 * * * * HOMEGA =13.538 CNILSS =-2.166 DNILSS = 0.000 * * * * OSCILLATOR FREQUENCIES: X=13.537901720 Y=13.537901720 Z=13.537901720 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: NILSSON NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 1) -48.614 ( 1 ----- 0) | 0, 0, 0, 1/2> 100 -0.500 -0.500 1.000 * * 2) -37.242 ( 0 ----- 1) | 1, 1, 0, 1/2> -100 -1.000 -0.333 0.333 * * 3) -37.242 ( 0 ----- 2) | 1, 0, 1, 3/2> -100 0.000 0.000 0.333 * * 4) -30.744 ( 0 ----- 3) | 1, 0, 1, 1/2> -100 0.500 -0.167 -0.333 * * 5) -25.870 ( 2 ----- 0) | 2, 1, 1, 3/2> 100 -1.041 -0.208 0.200 * * 6) -25.870 ( 3 ----- 0) | 2, 0, 2, 5/2> 100 0.588 0.118 0.200 * * 7) -25.870 ( 4 ----- 0) | 2, 1, 1, 3/2> 100 -1.047 -0.209 0.200 * * 8) -21.538 ( 5 ----- 0) | 2, 0, 0, 1/2> 100 -0.500 -0.500 1.000 * * 9) -15.040 ( 6 ----- 0) | 2, 1, 1, 1/2> 100 1.277 -0.255 -0.200 * * 10) -15.040 ( 7 ----- 0) | 2, 0, 2, 3/2> 100 -0.277 0.055 -0.200 * * 11) -14.498 ( 0 ----- 4) | 3, 2, 1, 3/2> -100 0.156 0.022 0.143 * * 12) -14.498 ( 0 ----- 5) | 3, 1, 2, 5/2> -100 -0.039 -0.006 0.143 * * 13) -14.498 ( 0 ----- 6) | 3, 1, 0, 1/2> -100 -1.876 -0.268 0.143 * * 14) -14.498 ( 0 ----- 7) | 3, 0, 3, 7/2> -100 -0.241 -0.034 0.143 * * 15) -10.166 ( 0 ----- 8) | 3, 0, 1, 3/2> -100 -0.002 -0.001 0.333 * * 16) -10.166 ( 0 ----- 9) | 3, 3, 0, 1/2> -100 -0.998 -0.333 0.333 * * 17) -3.668 ( 0 ----- 10) | 3, 0, 1, 1/2> -100 0.500 -0.167 -0.333 * * 18) -3.126 ( 8 ----- 0) | 4, 3, 1, 3/2> 100 -1.375 -0.153 0.111 * * 19) -3.126 ( 9 ----- 0) | 4, 1, 3, 7/2> 100 0.958 0.106 0.111 * * 20) -3.126 (10 ----- 0) | 4, 0, 4, 9/2> 100 -1.940 -0.216 0.111 * * 21) -3.126 (11 ----- 0) | 4, 2, 2, 5/2> 100 -0.848 -0.094 0.111 * * 22) -3.126 (12 ----- 0) | 4, 0, 4, 9/2> 100 0.705 0.078 0.111 * * 23) 0.664 ( 0 ----- 11) | 3, 1, 2, 3/2> -100 -0.474 0.068 -0.143 * * 24) 0.664 ( 0 ----- 12) | 3, 0, 3, 5/2> -100 0.106 -0.015 -0.143 * * 25) 0.664 ( 0 ----- 13) | 3, 2, 1, 1/2> -100 1.868 -0.267 -0.143 * * 26) 1.206 (13 ----- 0) | 4, 1, 1, 3/2> 100 -0.712 -0.142 0.200 * * 27) 1.206 (14 ----- 0) | 4, 4, 0, 1/2> 100 -1.013 -0.203 0.200 * * 28) 1.206 (15 ----- 0) | 4, 0, 2, 5/2> 100 0.225 0.045 0.200 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] NEUTRONS * * * ******************************************************************************* * * * Q00 = 28.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 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS NEUTRONS * * * * R_RMS = 3.2408 X_RMS = 1.8711 Y_RMS = 1.8711 Z_RMS = 1.8711 * * * * R_GEO = 4.1839 X_GEO = 4.1839 Y_GEO = 4.1839 Z_GEO = 4.1839 * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) NEUTRONS * * * ******************************************************************************* * * * 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 * * * ******************************************************************************* * * * BOHR DENSITY = 0.0913 FM^(-3) BOHR RADIUS = 4.1839 FM * * * ******************************************************************************* ******************************************************************************* * * * PARAMETERS OF THE NILSSON MODEL PROTONS * * * ******************************************************************************* * * * ESCALE =55.383 KAPPA = 0.080 MIU = 0.000 * * * * HOMEGA =13.538 CNILSS =-2.166 DNILSS = 0.000 * * * * OSCILLATOR FREQUENCIES: X=13.537901720 Y=13.537901720 Z=13.537901720 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: NILSSON PROTONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 1) -35.076 ( 1 ----- 0) | 0, 0, 0, 1/2> 100 -0.500 -0.500 1.000 * * 2) -23.704 ( 0 ----- 1) | 1, 1, 0, 1/2> -100 -1.000 -0.333 0.333 * * 3) -23.704 ( 0 ----- 2) | 1, 0, 1, 3/2> -100 0.000 0.000 0.333 * * 4) -17.206 ( 0 ----- 3) | 1, 0, 1, 1/2> -100 0.500 -0.167 -0.333 * * 5) -12.332 ( 2 ----- 0) | 2, 2, 0, 1/2> 100 -0.010 -0.002 0.200 * * 6) -12.332 ( 3 ----- 0) | 2, 1, 1, 3/2> 100 -0.805 -0.161 0.200 * * 7) -12.332 ( 4 ----- 0) | 2, 0, 2, 5/2> 100 -0.685 -0.137 0.200 * * 8) -8.000 ( 5 ----- 0) | 2, 0, 0, 1/2> 100 -0.500 -0.500 1.000 * * 9) -1.502 ( 6 ----- 0) | 2, 1, 1, 1/2> 100 0.153 -0.031 -0.200 * * 10) -1.502 ( 7 ----- 0) | 2, 0, 2, 3/2> 100 0.847 -0.169 -0.200 * * 11) -0.960 ( 0 ----- 4) | 3, 1, 2, 5/2> -100 0.099 0.014 0.143 * * 12) -0.960 ( 0 ----- 5) | 3, 2, 1, 3/2> -100 0.810 0.116 0.143 * * 13) -0.960 ( 0 ----- 6) | 3, 0, 3, 7/2> -100 -0.996 -0.142 0.143 * * 14) -0.960 ( 0 ----- 7) | 3, 0, 3, 7/2> -100 -1.913 -0.273 0.143 * * 15) 3.372 ( 0 ----- 8) | 3, 3, 0, 1/2> -100 -0.914 -0.305 0.333 * * 16) 3.372 ( 0 ----- 9) | 3, 0, 1, 3/2> -100 -0.086 -0.029 0.333 * * 17) 9.870 ( 0 ----- 10) | 3, 0, 1, 1/2> -100 0.500 -0.167 -0.333 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] PROTONS * * * ******************************************************************************* * * * Q00 = 20.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 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS PROTONS * * * * R_RMS = 3.0315 X_RMS = 1.7502 Y_RMS = 1.7502 Z_RMS = 1.7502 * * * * R_GEO = 3.9136 X_GEO = 3.9136 Y_GEO = 3.9136 Z_GEO = 3.9136 * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) PROTONS * * * ******************************************************************************* * * * 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 * * * ******************************************************************************* * * * BOHR DENSITY = 0.0797 FM^(-3) BOHR RADIUS = 3.9136 FM * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 7.312763 3.728959 7.312763 0.145154 * * DRHOD = 5.495046 2.798037 5.495046 0.101028 * * DLPR_ = -3.035498 -1.566956 -3.035498 -0.098414 * * DTAU_ = 7.776399 3.987980 7.776399 0.199561 * * DSCU_ = 0.065609 0.065609 0.065609 0.065609 * * DDIV_ = 0.566955 0.294413 0.566955 0.021870 * * * * 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_ = -12895.289160 6190.462699 -6825.311835 120.485374 * * ERHOD = 10579.809393 -5955.752229 4731.578441 -107.521276 * * ELPR_ = 281.249049 -49.068083 233.721842 -1.540876 * * ETAU_ = 252.514832 196.661134 444.255441 4.920524 * * ESCU_ = 0.000000 0.000000 0.000000 0.000000 * * EDIV_ = -34.867763 -18.106375 -52.301645 -0.672493 * * ============ ============ ============ ============ * * SUM EVEN: -1816.583648 364.197146 -1468.057755 15.671253 * * * * 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 =180.00000 BETA = 1.36436 GAMMA = 90.00000 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 3.1553 X_RMS = 1.8217 Y_RMS = 1.8217 Z_RMS = 1.8217 * * * * R_GEO = 4.0735 X_GEO = 4.0735 Y_GEO = 4.0735 Z_GEO = 4.0735 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * ******************************************************************************* * * * Q00 = 48.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 * * * ******************************************************************************* ******************************************************************************* * * * 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 * * * ******************************************************************************* * * * BOHR DENSITY = 0.1695 FM^(-3) BOHR RADIUS = 4.0735 FM * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] [INTRINSIC FRAME] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 48.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = ZERO Q2+1= ZERO Q2+2= ZERO ............. ............. * * 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 = ZERO Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4= ZERO * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) [INTRINSIC FRAME] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 = ZERO B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4= ZERO * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1695 FM^(-3) BOHR RADIUS = 4.0735 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 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 * * * ******************************************************************************* * * * 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 N * * * * CONF: - - 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 N * * * ******************************************************************************* * * * ENERGIES (MEV) * * * ******************************************************************************* * * * KINETIC: (NEU)= 636.281381 (PRO)= 397.675863 (TOT)= 1033.957244 * * SUM EPS: (NEU)= -682.122594 (PRO)= -295.379017 (TOT)= -977.501611 * * 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.332129 (PRO)= -1.231049 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 89.953646 (EXC)= -8.426092 (TOT)= 81.527555 * * * * 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= 385.338097 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -336.901703 * * * * SPIN-ORB (EVE)= -52.974138 (ODD)= 0.000000 (TOT)= -52.974138 * * SKYRME: (EVE)= -1452.386502 (ODD)= 0.000000 (TOT)= -1452.386502 * * * * TOTAL: (STAB)= -23.017275 (SP)= -359.918978 (FUN)= -336.901703 * * * ******************************************************************************* ******************************************************************************* * * * CONVERGENCE REPORT * * * ******************************************************************************* * * * ITER ENERGY STABILITY Q_2 GAMMA SPIN OMEGA ANGLE EPAIR * * * * 0 -336.901703 -23.017275 0.000 15.80 0.00 0.000 0.000 0.00 * * 1 -298.218758 384.125957 0.000 56.37 0.00 0.000 0.000 0.00 * * 2 -384.458438 273.000566 0.000 46.22 0.00 0.000 0.000 0.00 * * 3 -405.037221 159.491431 0.000 46.64 0.00 0.000 0.000 0.00 * * 4 -412.387979 88.169208 0.000 46.93 0.00 0.000 0.000 0.00 * * 5 -415.410697 47.853564 0.000 47.05 0.00 0.000 0.000 0.00 * * 6 -416.714994 25.644843 0.000 47.12 0.00 0.000 0.000 0.00 * * 7 -417.285457 13.641512 0.000 47.15 0.00 0.000 0.000 0.00 * * 8 -417.536476 7.221303 0.000 47.18 0.00 0.000 0.000 0.00 * * 9 -417.647245 3.811385 0.000 47.20 0.00 0.000 0.000 0.00 * * 10 -417.696209 2.008426 0.000 47.21 0.00 0.000 0.000 0.00 * * 11 -417.717882 1.057771 0.000 47.21 0.00 0.000 0.000 0.00 * * 12 -417.727486 0.557243 0.000 47.22 0.00 0.000 0.000 0.00 * * 13 -417.731748 0.293827 0.000 47.22 0.00 0.000 0.000 0.00 * * 14 -417.733640 0.155145 0.000 47.23 0.00 0.000 0.000 0.00 * * 15 -417.734483 0.082058 0.000 47.23 0.00 0.000 0.000 0.00 * * 16 -417.734858 0.043482 0.000 47.23 0.00 0.000 0.000 0.00 * * 17 -417.735025 0.023083 0.000 47.23 0.00 0.000 0.000 0.00 * * 18 -417.735100 0.012273 0.000 47.23 0.00 0.000 0.000 0.00 * * 19 -417.735134 0.006533 0.000 47.23 0.00 0.000 0.000 0.00 * * 20 -417.735149 0.003479 0.000 47.23 0.00 0.000 0.000 0.00 * * 21 -417.735156 0.001851 0.000 47.23 0.00 0.000 0.000 0.00 * * 22 -417.735159 0.000983 0.000 47.23 0.00 0.000 0.000 0.00 * * 23 -417.735160 0.000519 0.000 47.23 0.00 0.000 0.000 0.00 * * 24 -417.735161 0.000272 0.000 47.23 0.00 0.000 0.000 0.00 * * 25 -417.735161 0.000141 0.000 47.23 0.00 0.000 0.000 0.00 * * 26 -417.735162 0.000072 0.000 47.23 0.00 0.000 0.000 0.00 * * 27 -417.735162 0.000035 0.000 47.23 0.00 0.000 0.000 0.00 * * 28 -417.735162 0.000017 0.000 47.23 0.00 0.000 0.000 0.00 * * 29 -417.735162 7.17E-06 0.000 47.23 0.00 0.000 0.000 0.00 * * 30 -417.735162 2.54E-06 0.000 47.23 0.00 0.000 0.000 0.00 * * 31 -417.735162 4.09E-07 0.000 47.23 0.00 0.000 0.000 0.00 * * 32 -417.735162 -4.69E-07 0.000 47.23 0.00 0.000 0.000 0.00 * * 33 -417.735162 -7.45E-07 0.000 47.23 0.00 0.000 0.000 0.00 * * 34 -417.735162 -7.49E-07 0.000 47.23 0.00 0.000 0.000 0.00 * * 35 -417.735162 -6.49E-07 0.000 47.23 0.00 0.000 0.000 0.00 * * 36 -417.735162 -5.21E-07 0.000 47.23 0.00 0.000 0.000 0.00 * * 37 -417.735162 -4.00E-07 0.000 47.23 0.00 0.000 0.000 0.00 * * 38 -417.735162 -2.98E-07 0.000 47.23 0.00 0.000 0.000 0.00 * * 39 -417.735162 -2.18E-07 0.000 47.23 0.00 0.000 0.000 0.00 * * 40 -417.735162 -1.57E-07 0.000 47.23 0.00 0.000 0.000 0.00 * * 41 -417.735162 -1.12E-07 0.000 47.23 0.00 0.000 0.000 0.00 * * 42 -417.735162 -7.86E-08 0.000 47.23 0.00 0.000 0.000 0.00 * * 43 -417.735162 -5.50E-08 0.000 47.23 0.00 0.000 0.000 0.00 * * 44 -417.735162 -3.81E-08 0.000 47.23 0.00 0.000 0.000 0.00 * * 45 -417.735162 -2.63E-08 0.000 47.23 0.00 0.000 0.000 0.00 * * 46 -417.735162 -1.83E-08 0.000 47.23 0.00 0.000 0.000 0.00 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 1) -49.784 ( 1 ----- 0) | 0, 0, 0, 1/2> 100 -0.500 -0.500 1.000 * * 2) -36.183 ( 0 ----- 1) | 1, 1, 0, 1/2> -100 0.179 0.060 0.333 * * 3) -36.183 ( 0 ----- 2) | 1, 0, 1, 3/2> -100 -1.179 -0.393 0.333 * * 4) -33.338 ( 0 ----- 3) | 1, 0, 1, 1/2> -100 0.500 -0.167 -0.333 * * 5) -22.599 ( 2 ----- 0) | 2, 1, 1, 3/2> 100 0.834 0.167 0.200 * * 6) -22.599 ( 3 ----- 0) | 2, 0, 2, 5/2> 100 -0.715 -0.143 0.200 * * 7) -22.599 ( 4 ----- 0) | 2, 0, 2, 5/2> 100 -1.619 -0.324 0.200 * * 8) -17.603 ( 5 ----- 0) | 2, 0, 0, 1/2> 100 -0.500 -0.500 1.000 * * 9) -16.547 ( 6 ----- 0) | 2, 1, 1, 1/2> 100 -0.179 0.036 -0.200 * * 10) -16.547 ( 7 ----- 0) | 2, 0, 2, 3/2> 100 1.179 -0.236 -0.200 * * 11) -9.793 ( 0 ----- 4) | 3, 1, 2, 5/2> -100 1.499 0.214 0.143 * * 12) -9.793 ( 0 ----- 5) | 3, 2, 1, 3/2> -100 -0.706 -0.101 0.143 * * 13) -9.793 ( 0 ----- 6) | 3, 0, 3, 7/2> -100 -1.613 -0.230 0.143 * * 14) -9.793 ( 0 ----- 7) | 3, 0, 3, 7/2> -100 -1.179 -0.168 0.143 * * 15) -5.535 ( 0 ----- 8) | 3, 0, 1, 3/2> -100 0.179 0.060 0.333 * * 16) -5.535 ( 0 ----- 9) | 3, 0, 1, 3/2> -100 -1.179 -0.393 0.333 * * 17) -3.539 ( 0 ----- 10) | 3, 0, 1, 1/2> -100 0.500 -0.167 -0.333 * * 18) -1.334 ( 0 ----- 11) | 3, 1, 2, 3/2> -100 -0.835 0.119 -0.143 * * 19) -1.334 ( 0 ----- 12) | 3, 2, 1, 1/2> -100 0.715 -0.102 -0.143 * * 20) -1.334 ( 0 ----- 13) | 3, 0, 3, 5/2> -100 1.619 -0.231 -0.143 * * 21) 2.125 ( 8 ----- 0) | 4, 2, 2, 5/2> 100 -0.646 -0.072 0.111 * * 22) 2.125 ( 9 ----- 0) | 4, 1, 3, 7/2> 100 2.461 0.273 0.111 * * 23) 2.125 (10 ----- 0) | 4, 0, 4, 9/2> 100 -0.130 -0.014 0.111 * * 24) 2.125 (11 ----- 0) | 4, 0, 4, 9/2> 100 -2.342 -0.260 0.111 * * 25) 2.125 (12 ----- 0) | 4, 0, 4, 9/2> 100 -1.843 -0.205 0.111 * * 26) 3.802 (13 ----- 0) | 4, 0, 0, 1/2> 100 -0.500 -0.500 1.000 * * 27) 4.024 (14 ----- 0) | 4, 1, 1, 3/2> 100 0.834 0.167 0.200 * * 28) 4.024 (15 ----- 0) | 4, 0, 2, 5/2> 100 -0.713 -0.143 0.200 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] NEUTRONS * * * ******************************************************************************* * * * Q00 = 28.0000 ............. ............. ............. ............. * * * * Q10 = ZERO Q11 = ZERO ............. ............. ............. * * * * Q20 = ZERO Q21 = ZERO Q22 =-3.0E-10 ............. ............. * * * * Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. * * * * Q40 =-1.4E-08 Q41 = ZERO Q42 = ZERO Q43 = ZERO Q44 =-8.5E-09 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS NEUTRONS * * * * R_RMS = 3.6062 X_RMS = 2.0821 Y_RMS = 2.0821 Z_RMS = 2.0821 * * * * R_GEO = 4.6556 X_GEO = 4.6556 Y_GEO = 4.6556 Z_GEO = 4.6556 * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) NEUTRONS * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = ZERO B21 = ZERO B22 = ZERO ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 =-4.5E-08 B41 = ZERO B42 = ZERO B43 = ZERO B44 =-3.8E-08 * * * ******************************************************************************* * * * BOHR DENSITY = 0.0662 FM^(-3) BOHR RADIUS = 4.6556 FM * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEz>

GFACT * * * * 1) -46.032 ( 1 ----- 0) | 0, 0, 0, 1/2> 100 -0.500 -0.500 1.000 * * 2) -34.173 ( 0 ----- 1) | 1, 1, 0, 1/2> -100 0.179 0.060 0.333 * * 3) -34.173 ( 0 ----- 2) | 1, 0, 1, 3/2> -100 -1.179 -0.393 0.333 * * 4) -31.023 ( 0 ----- 3) | 1, 0, 1, 1/2> -100 0.500 -0.167 -0.333 * * 5) -22.103 ( 2 ----- 0) | 2, 1, 1, 3/2> 100 0.835 0.167 0.200 * * 6) -22.103 ( 3 ----- 0) | 2, 0, 2, 5/2> 100 -0.715 -0.143 0.200 * * 7) -22.103 ( 4 ----- 0) | 2, 0, 2, 5/2> 100 -1.620 -0.324 0.200 * * 8) -16.814 ( 5 ----- 0) | 2, 0, 0, 1/2> 100 -0.500 -0.500 1.000 * * 9) -16.048 ( 6 ----- 0) | 2, 1, 1, 1/2> 100 -0.179 0.036 -0.200 * * 10) -16.048 ( 7 ----- 0) | 2, 0, 2, 3/2> 100 1.179 -0.236 -0.200 * * 11) -10.226 ( 0 ----- 4) | 3, 1, 2, 5/2> -100 1.499 0.214 0.143 * * 12) -10.226 ( 0 ----- 5) | 3, 2, 1, 3/2> -100 -0.708 -0.101 0.143 * * 13) -10.226 ( 0 ----- 6) | 3, 0, 3, 7/2> -100 -1.617 -0.231 0.143 * * 14) -10.226 ( 0 ----- 7) | 3, 0, 3, 7/2> -100 -1.174 -0.168 0.143 * * 15) -4.813 ( 0 ----- 8) | 3, 0, 1, 3/2> -100 0.179 0.060 0.333 * * 16) -4.813 ( 0 ----- 9) | 3, 0, 1, 3/2> -100 -1.179 -0.393 0.333 * * 17) -2.807 ( 0 ----- 10) | 3, 0, 1, 1/2> -100 0.500 -0.167 -0.333 * * 18) -2.024 ( 0 ----- 11) | 3, 1, 2, 3/2> -100 -0.835 0.119 -0.143 * * 19) -2.024 ( 0 ----- 12) | 3, 2, 1, 1/2> -100 0.715 -0.102 -0.143 * * 20) -2.024 ( 0 ----- 13) | 3, 0, 3, 5/2> -100 1.619 -0.231 -0.143 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] PROTONS * * * ******************************************************************************* * * * Q00 = 20.0000 ............. ............. ............. ............. * * * * Q10 = ZERO Q11 = ZERO ............. ............. ............. * * * * Q20 = ZERO Q21 = ZERO Q22 =-1.5E-10 ............. ............. * * * * Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. * * * * Q40 =-7.8E-09 Q41 = ZERO Q42 = ZERO Q43 = ZERO Q44 =-4.7E-09 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS PROTONS * * * * R_RMS = 3.4541 X_RMS = 1.9942 Y_RMS = 1.9942 Z_RMS = 1.9942 * * * * R_GEO = 4.4593 X_GEO = 4.4593 Y_GEO = 4.4593 Z_GEO = 4.4593 * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) PROTONS * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = ZERO B21 = ZERO B22 = ZERO ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 =-4.1E-08 B41 = ZERO B42 = ZERO B43 = ZERO B44 =-3.5E-08 * * * ******************************************************************************* * * * BOHR DENSITY = 0.0538 FM^(-3) BOHR RADIUS = 4.4593 FM * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 5.144197 2.626106 5.144197 0.108015 * * DRHOD = 3.633662 1.853599 3.633662 0.073536 * * DLPR_ = -1.749188 -0.896397 -1.749188 -0.043605 * * DTAU_ = 4.229355 2.185379 4.229355 0.141403 * * DSCU_ = 0.040516 0.039138 0.040516 0.037761 * * DDIV_ = 0.339868 0.185889 0.339868 0.031909 * * * * 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_ = -9071.249938 4359.611277 -4801.296721 89.658060 * * ERHOD = 6996.018768 -3945.471112 3128.809825 -78.262168 * * ELPR_ = 162.068141 -28.070002 134.680863 -0.682724 * * ETAU_ = 137.335397 107.768629 241.617480 3.486546 * * ESCU_ = 0.000000 0.000000 0.000000 0.000000 * * EDIV_ = -20.901880 -11.432143 -31.352820 -0.981203 * * ============ ============ ============ ============ * * SUM EVEN: -1796.729511 482.406650 -1327.541373 13.218512 * * * * 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 = 90.00000 GAMMA =179.99993 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 3.5437 X_RMS = 2.0459 Y_RMS = 2.0459 Z_RMS = 2.0459 * * * * R_GEO = 4.5748 X_GEO = 4.5748 Y_GEO = 4.5748 Z_GEO = 4.5748 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * ******************************************************************************* * * * Q00 = 48.0000 ............. ............. ............. ............. * * * * Q10 = ZERO Q11 = ZERO ............. ............. ............. * * * * Q20 = 1.4E-10 Q21 = ZERO Q22 =-4.5E-10 ............. ............. * * * * Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. * * * * Q40 =-2.2E-08 Q41 = ZERO Q42 = ZERO Q43 = ZERO Q44 =-1.3E-08 * * * ******************************************************************************* ******************************************************************************* * * * 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 =-4.4E-08 B41 = ZERO B42 = ZERO B43 = ZERO B44 =-3.7E-08 * * * ******************************************************************************* * * * BOHR DENSITY = 0.1197 FM^(-3) BOHR RADIUS = 4.5748 FM * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] [INTRINSIC FRAME] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 48.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 3.2E-10 Q2+1= ZERO Q2+2= 3.5E-10 ............. ............. * * 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 =-2.2E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4=-1.3E-08 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) [INTRINSIC FRAME] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 =-4.4E-08 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4=-3.7E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1197 FM^(-3) BOHR RADIUS = 4.5748 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 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 * * * ******************************************************************************* * * * 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 N * * * * CONF: - - 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 2 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 N * * * ******************************************************************************* * * * ENERGIES (MEV) * * * ******************************************************************************* * * * KINETIC: (NEU)= 516.414796 (PRO)= 309.027310 (TOT)= 825.442105 * * SUM EPS: (NEU)= -626.310864 (PRO)= -521.234772 (TOT)= -1147.545636 * * 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.664247 (PRO)= -13.136820 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 78.558869 (EXC)= -7.413275 (TOT)= 71.145594 * * * * 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= 254.212305 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -417.735162 * * * * SPIN-ORB (EVE)= -32.334023 (ODD)= 0.000000 (TOT)= -32.334023 * * SKYRME: (EVE)= -1314.322862 (ODD)= 0.000000 (TOT)= -1314.322862 * * * * TOTAL: (STAB)= -1.26E-08 (SP)= -417.735162 (FUN)= -417.735162 * * * ******************************************************************************* ******************************************************************************* * * * CALLING SUBROUTINE "RECORD" TO SAVE THE RECORD FILE CONTAINING RESULTS * * ca48_lip.rec * * THE MEAN FIELDS OBTAINED IN THE LAST-BUT-ONE ITERATION HAVE BEEN STORED * * * * CALLING SUBROUTINE "RECFIL" TO SAVE THE FIELDS FILE CONTAINING RESULTS * * ca48_lip.fic * * * ******************************************************************************* ******************************************************************************* * * * CALCULATION OF THE DATA SET NO. 1 HAS BEEN SUCCESSFULLY COMPLETED * * * ******************************************************************************* ******************************************************************************* * * * CLASSICAL NUCLEAR SURFACE DEFINED FOR: N = 28 Z = 20 * * * ******************************************************************************* * * * SPHERICAL SHAPE REQUESTED FOR THE SURFACE * * * ******************************************************************************* * * * HOMEGA= 13.5379 FCHOM0= 1.2000 * * * * OSCILLATOR FREQUENCIES: HBAROX= 13.5379 HBAROY= 13.5379 HBAROZ= 13.5379 * * * * MOMENTS OF INERTIA: XMOMFC= 9.2482 YMOMFC= 9.2482 ZMOMFC= 9.2482 * * * * CENTERS OF MASS: CMSXFC= 0.0000 CMSYFC= 0.0000 CMSZFC= 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * SKYRME FORCE DEFINITION * * * ******************************************************************************* * * * PARAMETER SET SLY4: T0= -2488.91 T1= 486.82 T2= -546.39 T3= 13777.00 * * * * POWER=1/6 W=123.000 X0= 0.83400 X1= -0.34400 X2= -1.00000 X3= 1.35400 * * * * ETA_J=0 ETA_W=0 ETACM=0 ETA_M=1 HBM= 20.73553 (SKYRME-FORCE SPIN-ORBIT) * * * ******************************************************************************* ******************************************************************************* * * * NON-STANDARD TREATMENT OF THE SKYRME FORCE: ETA_J=0 ETA_W=0 ETACM=2 ETA_M=1 * * * * THIS INCLUDES: * * * * NON-STANDARD TREATMENT OF THE CENTER-OF-MASS CORRECTION FOR THIS FORCE: * * METHOD BASED ON THE TWO-BODY OPERATOR INCLUDED BEFORE VARIATION IS NOW USED * * * ******************************************************************************* ******************************************************************************* * * * COUPLING CONSTANTS DEFINING THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * CRHO_ = -1763.394861 1660.104971 -933.342375 830.052485 * * CRHOD = 1925.335750 -2128.546500 861.062500 -1064.273250 * * CLPR_ = -92.653338 31.314270 -76.996203 15.657135 * * CTAU_ = 32.471951 49.313473 57.128687 24.656736 * * CSCU_ = 0.000000 0.000000 0.000000 0.000000 * * CDIV_ = -61.500000 -61.500000 -92.250000 -30.750000 * * * * CSPI_ = -518.938360 622.228250 -207.824235 311.114125 * * CSPID = 777.252417 -574.041667 490.231583 -287.020833 * * CLPS_ = 32.774724 28.564344 47.056896 14.282172 * * CCUR_ = -32.471951 -49.313473 -57.128687 -24.656736 * * CKIS_ = 0.000000 0.000000 0.000000 0.000000 * * CROT_ = -61.500000 -61.500000 -92.250000 -30.750000 * * * ******************************************************************************* ******************************************************************************* * * * PHYSICAL CONSTANTS: H_BARC=197.32891000 HBCOE2=137.03602000 * * * * XMASSN=938.90590000 XMASSP=938.27231000 * * * * HBMASS= 20.73553000 HBMRPA= 20.73553000 * * * * ECHAR2= 1.43997841 COULEX= -1.06350868 * * * ******************************************************************************* * * * OSCILLATOR LENGTHS: X= 1.7502658 Y= 1.7502658 Z= 1.7502658 * * * * OSCILLATOR CONSTANTS: X= 0.5713418 Y= 0.5713418 Z= 0.5713418 * * * * OSCILLATOR FREQUENCIES: X=13.5379017 Y=13.5379017 Z=13.5379017 * * * ******************************************************************************* * * * 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= 209.8375 * * * ******************************************************************************* ******************************************************************************* * * * 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): * * ca48_lip.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): * * ca48_lip.cou * * * ******************************************************************************* * * * NUCLIDE: N = 28 Z = 20 NUMBER OF ITERATIONS = 200 CONTINUATION? = 1 * * * ******************************************************************************* * * * ITERATIONS WILL STOP WHEN THE STABILITY BECOMES LOWER THAN EPSITE= 1.0E-07 * * FOR 5 CONSECUTIVE ITERATIONS * * * ******************************************************************************* * * * ITERATIONS WILL STOP WHEN THE PING-PONG DIVERGENCE OCCURS FOR EPSPNG=0.010 * * DURING 3 PAIRS OF CONSECUTIVE ITERATIONS * * * ******************************************************************************* * * * ITERATIONS WILL STOP WHEN THE CHAOTIC-DIVERGENCE CONDITION OCCURS 5 TIMES * * * ******************************************************************************* * * * SLOW-DOWN PARAMETERS: * * TIME-EVEN = 0.50 TIME-ODD = 0.50 PAIRING = 0.50 LIPKIN-NOGAMI = 0.50 * * * ******************************************************************************* * * * BROYDEN METHOD IS: OFF * * * * ITERATIONS WILL PROCEED BY STANDARD LINEAR MIXING WITH SLOWING FACTORS * * DISPLAYED EARLIER * * * ******************************************************************************* * * * 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 * * * ******************************************************************************* * * * MEAN-FIELD MATRIX ELEMENTS WILL BE READ FROM DISC TO CONTINUE ITERATIONS * * * ******************************************************************************* * * * MAXIMUM NUMBERS OF MULTIPOLE MOMENTS USED FOR THE CONSTRAINTS, NMUCON=3 * * FOR THE COULOMB FIELD, NMUCOU=4 * * FOR THE OUTPUT INFO, NMUPRI=4 * * * ******************************************************************************* * * * NUCLEON MASS RENORMALIZED IRENMA= 1 DISTAX=0.00 DISTAY=0.00 DISTAZ=2.00 * * * ******************************************************************************* * * * NUCLEON MASS RENORMALIZATION INITIALIZED WITH FACTORS X=1.00 Y=1.00 Z=1.00 * * * ******************************************************************************* * * * PRINTING THE RESULTS FOR THE FOLLOWING ITERATIONS: THE FIRST: YES * * THE MIDDLE: NO * * AND/OR THE LAST: YES * * * ******************************************************************************* * * * PRINTING THE MOMENTS, RADII, AND DEFORMATIONS FOR: NEUTRONS: YES * * PROTONS : YES * * MASS : YES * * * ******************************************************************************* * * * PRINTING THE MOMENTS, MOMENTA AND DEFORMATIONS IN THE INTRINSIC FRAME: YES * * * ******************************************************************************* * * * 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: 28 28 * * PROTONS: 20 20 * * * ******************************************************************************* * * * 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.000000 (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: 28 0 0 * * PROTONS : 20 0 0 * * * ******************************************************************************* ******************************************************************************* * * * CALLING SUBROUTINE "RECORD" TO READ THE REPLAY FILE OBTAINED PREVIOUSLY * * ca48_lip.rec * * * * CALLING SUBROUTINE "RECFIL" TO READ THE FIELDS FILE OBTAINED PREVIOUSLY * * ca48_lip.fic * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY | N,nz,lz,OMEz>

* * * * 1) -49.784 1.0000 | 0, 0, 0, 1/2> 100 0.000 -0.500 0.000 * * 0.000 -0.500 0.000 * * 2) -49.784 -1.0000 | 0, 0, 0, 1/2> 100 0.000 0.500 0.000 * * 0.000 0.500 0.000 * * 3) -36.183 1.0000 | 1, 1, 0, 1/2> -100 0.000 0.179 0.000 * * 0.000 0.060 0.000 * * 4) -36.183 -1.0000 | 1, 1, 0, 1/2> -100 0.000 -0.179 0.000 * * 0.000 -0.060 0.000 * * 5) -36.183 -1.0000 | 1, 0, 1, 3/2> -100 0.000 1.179 0.000 * * 0.000 0.393 0.000 * * 6) -36.183 1.0000 | 1, 0, 1, 3/2> -100 0.000 -1.179 0.000 * * 0.000 -0.393 0.000 * * 7) -33.338 -1.0000 | 1, 0, 1, 1/2> -100 0.000 -0.500 0.000 * * 0.000 0.167 0.000 * * 8) -33.338 1.0000 | 1, 0, 1, 1/2> -100 0.000 0.500 0.000 * * 0.000 -0.167 0.000 * * 9) -22.599 -1.0000 | 2, 1, 1, 3/2> 100 0.000 -0.834 0.000 * * 0.000 -0.167 0.000 * * 10) -22.599 1.0000 | 2, 1, 1, 3/2> 100 0.000 0.834 0.000 * * 0.000 0.167 0.000 * * 11) -22.599 -1.0000 | 2, 0, 2, 5/2> 100 0.000 0.715 0.000 * * 0.000 0.143 0.000 * * 12) -22.599 1.0000 | 2, 0, 2, 5/2> 100 0.000 -0.715 0.000 * * 0.000 -0.143 0.000 * * 13) -22.599 1.0000 | 2, 0, 2, 5/2> 100 0.000 -1.619 0.000 * * 0.000 -0.324 0.000 * * 14) -22.599 -1.0000 | 2, 0, 2, 5/2> 100 0.000 1.619 0.000 * * 0.000 0.324 0.000 * * 15) -17.603 1.0000 | 2, 0, 0, 1/2> 100 0.000 -0.500 0.000 * * 0.000 -0.500 0.000 * * 16) -17.603 -1.0000 | 2, 0, 0, 1/2> 100 0.000 0.500 0.000 * * 0.000 0.500 0.000 * * 17) -16.547 -1.0000 | 2, 1, 1, 1/2> 100 0.000 0.179 0.000 * * 0.000 -0.036 0.000 * * 18) -16.547 1.0000 | 2, 1, 1, 1/2> 100 0.000 -0.179 0.000 * * 0.000 0.036 0.000 * * 19) -16.547 1.0000 | 2, 0, 2, 3/2> 100 0.000 1.179 0.000 * * 0.000 -0.236 0.000 * * 20) -16.547 -1.0000 | 2, 0, 2, 3/2> 100 0.000 -1.179 0.000 * * 0.000 0.236 0.000 * * 21) -9.793 1.0000 | 3, 1, 2, 5/2> -100 0.000 1.499 0.000 * * 0.000 0.214 0.000 * * 22) -9.793 -1.0000 | 3, 1, 2, 5/2> -100 0.000 -1.499 0.000 * * 0.000 -0.214 0.000 * * 23) -9.793 -1.0000 | 3, 2, 1, 3/2> -100 0.000 0.706 0.000 * * 0.000 0.101 0.000 * * 24) -9.793 1.0000 | 3, 2, 1, 3/2> -100 0.000 -0.706 0.000 * * 0.000 -0.101 0.000 * * 25) -9.793 -1.0000 | 3, 0, 3, 7/2> -100 0.000 1.613 0.000 * * 0.000 0.230 0.000 * * 26) -9.793 1.0000 | 3, 0, 3, 7/2> -100 0.000 -1.613 0.000 * * 0.000 -0.230 0.000 * * 27) -9.793 1.0000 | 3, 0, 3, 7/2> -100 0.000 -1.179 0.000 * * 0.000 -0.168 0.000 * * 28) -9.793 -1.0000 | 3, 0, 3, 7/2> -100 0.000 1.179 0.000 * * 0.000 0.168 0.000 * * 29) -5.535 1.0000 | 3, 0, 1, 3/2> -100 0.000 0.179 0.000 * * 0.000 0.060 0.000 * * 30) -5.535 -1.0000 | 3, 0, 1, 3/2> -100 0.000 -0.179 0.000 * * 0.000 -0.060 0.000 * * 31) -5.535 1.0000 | 3, 0, 1, 3/2> -100 0.000 -1.179 0.000 * * 0.000 -0.393 0.000 * * 32) -5.535 -1.0000 | 3, 0, 1, 3/2> -100 0.000 1.179 0.000 * * 0.000 0.393 0.000 * * 33) -3.539 1.0000 | 3, 0, 1, 1/2> -100 0.000 0.500 0.000 * * 0.000 -0.167 0.000 * * 34) -3.539 -1.0000 | 3, 0, 1, 1/2> -100 0.000 -0.500 0.000 * * 0.000 0.167 0.000 * * 35) -1.334 -1.0000 | 3, 1, 2, 3/2> -100 0.000 0.835 0.000 * * 0.000 -0.119 0.000 * * 36) -1.334 1.0000 | 3, 1, 2, 3/2> -100 0.000 -0.835 0.000 * * 0.000 0.119 0.000 * * 37) -1.334 -1.0000 | 3, 2, 1, 1/2> -100 0.000 -0.715 0.000 * * 0.000 0.102 0.000 * * 38) -1.334 1.0000 | 3, 2, 1, 1/2> -100 0.000 0.715 0.000 * * 0.000 -0.102 0.000 * * 39) -1.334 -1.0000 | 3, 0, 3, 5/2> -100 0.000 -1.619 0.000 * * 0.000 0.231 0.000 * * 40) -1.334 1.0000 | 3, 0, 3, 5/2> -100 0.000 1.619 0.000 * * 0.000 -0.231 0.000 * * 41) 2.125 -1.0000 | 4, 2, 2, 5/2> 100 0.000 0.646 0.000 * * 0.000 0.072 0.000 * * 42) 2.125 1.0000 | 4, 2, 2, 5/2> 100 0.000 -0.646 0.000 * * 0.000 -0.072 0.000 * * 43) 2.125 1.0000 | 4, 1, 3, 7/2> 100 0.000 2.462 0.000 * * 0.000 0.274 0.000 * * 44) 2.125 -1.0000 | 4, 1, 3, 7/2> 100 0.000 -2.461 0.000 * * 0.000 -0.273 0.000 * * 45) 2.125 1.0000 | 4, 0, 4, 9/2> 100 0.000 -0.130 0.000 * * 0.000 -0.014 0.000 * * 46) 2.125 -1.0000 | 4, 0, 4, 9/2> 100 0.000 0.130 0.000 * * 0.000 0.014 0.000 * * 47) 2.125 1.0000 | 4, 0, 4, 9/2> 100 0.000 -2.342 0.000 * * 0.000 -0.260 0.000 * * 48) 2.125 -1.0000 | 4, 0, 4, 9/2> 100 0.000 2.342 0.000 * * 0.000 0.260 0.000 * * 49) 2.125 1.0000 | 4, 0, 4, 9/2> 100 0.000 -1.843 0.000 * * 0.000 -0.205 0.000 * * 50) 2.125 -1.0000 | 4, 0, 4, 9/2> 100 0.000 1.843 0.000 * * 0.000 0.205 0.000 * * 51) 3.802 -1.0000 | 4, 0, 0, 1/2> 100 0.000 0.500 0.000 * * 0.000 0.500 0.000 * * 52) 3.802 1.0000 | 4, 0, 0, 1/2> 100 0.000 -0.500 0.000 * * 0.000 -0.500 0.000 * * 53) 4.024 1.0000 | 4, 1, 1, 3/2> 100 0.000 0.834 0.000 * * 0.000 0.167 0.000 * * 54) 4.024 -1.0000 | 4, 1, 1, 3/2> 100 0.000 -0.834 0.000 * * 0.000 -0.167 0.000 * * 55) 4.024 -1.0000 | 4, 0, 2, 5/2> 100 0.000 0.713 0.000 * * 0.000 0.143 0.000 * * 56) 4.024 1.0000 | 4, 0, 2, 5/2> 100 0.000 -0.713 0.000 * * 0.000 -0.143 0.000 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] NEUTRONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 28.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = ZERO Q2+1= ZERO Q2+2=-3.0E-10 ............. ............. * * 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 =-1.4E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4=-8.5E-09 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS NEUTRONS * * * * R_RMS = 3.6062 X_RMS = 2.0821 Y_RMS = 2.0821 Z_RMS = 2.0821 * * * * R_GEO = 4.6556 X_GEO = 4.6556 Y_GEO = 4.6556 Z_GEO = 4.6556 * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) NEUTRONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 =-4.5E-08 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4=-3.8E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.0662 FM^(-3) BOHR RADIUS = 4.6556 FM * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS * * * ******************************************************************************* * * * NO) ENERGY | N,nz,lz,OMEz>

* * * * 1) -46.032 1.0000 | 0, 0, 0, 1/2> 100 0.000 -0.500 0.000 * * 0.000 -0.500 0.000 * * 2) -46.032 -1.0000 | 0, 0, 0, 1/2> 100 0.000 0.500 0.000 * * 0.000 0.500 0.000 * * 3) -34.173 1.0000 | 1, 1, 0, 1/2> -100 0.000 0.179 0.000 * * 0.000 0.060 0.000 * * 4) -34.173 -1.0000 | 1, 1, 0, 1/2> -100 0.000 -0.179 0.000 * * 0.000 -0.060 0.000 * * 5) -34.173 -1.0000 | 1, 0, 1, 3/2> -100 0.000 1.179 0.000 * * 0.000 0.393 0.000 * * 6) -34.173 1.0000 | 1, 0, 1, 3/2> -100 0.000 -1.179 0.000 * * 0.000 -0.393 0.000 * * 7) -31.023 1.0000 | 1, 0, 1, 1/2> -100 0.000 0.500 0.000 * * 0.000 -0.167 0.000 * * 8) -31.023 -1.0000 | 1, 0, 1, 1/2> -100 0.000 -0.500 0.000 * * 0.000 0.167 0.000 * * 9) -22.103 1.0000 | 2, 1, 1, 3/2> 100 0.000 0.835 0.000 * * 0.000 0.167 0.000 * * 10) -22.103 -1.0000 | 2, 1, 1, 3/2> 100 0.000 -0.835 0.000 * * 0.000 -0.167 0.000 * * 11) -22.103 1.0000 | 2, 0, 2, 5/2> 100 0.000 -0.715 0.000 * * 0.000 -0.143 0.000 * * 12) -22.103 -1.0000 | 2, 0, 2, 5/2> 100 0.000 0.715 0.000 * * 0.000 0.143 0.000 * * 13) -22.103 -1.0000 | 2, 0, 2, 5/2> 100 0.000 1.620 0.000 * * 0.000 0.324 0.000 * * 14) -22.103 1.0000 | 2, 0, 2, 5/2> 100 0.000 -1.620 0.000 * * 0.000 -0.324 0.000 * * 15) -16.814 -1.0000 | 2, 0, 0, 1/2> 100 0.000 0.500 0.000 * * 0.000 0.500 0.000 * * 16) -16.814 1.0000 | 2, 0, 0, 1/2> 100 0.000 -0.500 0.000 * * 0.000 -0.500 0.000 * * 17) -16.048 -1.0000 | 2, 1, 1, 1/2> 100 0.000 0.179 0.000 * * 0.000 -0.036 0.000 * * 18) -16.048 1.0000 | 2, 1, 1, 1/2> 100 0.000 -0.179 0.000 * * 0.000 0.036 0.000 * * 19) -16.048 1.0000 | 2, 0, 2, 3/2> 100 0.000 1.179 0.000 * * 0.000 -0.236 0.000 * * 20) -16.048 -1.0000 | 2, 0, 2, 3/2> 100 0.000 -1.179 0.000 * * 0.000 0.236 0.000 * * 21) -10.226 1.0000 | 3, 1, 2, 5/2> -100 0.000 1.499 0.000 * * 0.000 0.214 0.000 * * 22) -10.226 -1.0000 | 3, 1, 2, 5/2> -100 0.000 -1.499 0.000 * * 0.000 -0.214 0.000 * * 23) -10.226 1.0000 | 3, 2, 1, 3/2> -100 0.000 -0.708 0.000 * * 0.000 -0.101 0.000 * * 24) -10.226 -1.0000 | 3, 2, 1, 3/2> -100 0.000 0.708 0.000 * * 0.000 0.101 0.000 * * 25) -10.226 1.0000 | 3, 0, 3, 7/2> -100 0.000 -1.616 0.000 * * 0.000 -0.231 0.000 * * 26) -10.226 -1.0000 | 3, 0, 3, 7/2> -100 0.000 1.616 0.000 * * 0.000 0.231 0.000 * * 27) -10.226 1.0000 | 3, 0, 3, 7/2> -100 0.000 -1.174 0.000 * * 0.000 -0.168 0.000 * * 28) -10.226 -1.0000 | 3, 0, 3, 7/2> -100 0.000 1.174 0.000 * * 0.000 0.168 0.000 * * 29) -4.813 -1.0000 | 3, 0, 1, 3/2> -100 0.000 -0.179 0.000 * * 0.000 -0.060 0.000 * * 30) -4.813 1.0000 | 3, 0, 1, 3/2> -100 0.000 0.179 0.000 * * 0.000 0.060 0.000 * * 31) -4.813 -1.0000 | 3, 0, 1, 3/2> -100 0.000 1.179 0.000 * * 0.000 0.393 0.000 * * 32) -4.813 1.0000 | 3, 0, 1, 3/2> -100 0.000 -1.179 0.000 * * 0.000 -0.393 0.000 * * 33) -2.807 -1.0000 | 3, 0, 1, 1/2> -100 0.000 -0.500 0.000 * * 0.000 0.167 0.000 * * 34) -2.807 1.0000 | 3, 0, 1, 1/2> -100 0.000 0.500 0.000 * * 0.000 -0.167 0.000 * * 35) -2.024 1.0000 | 3, 1, 2, 3/2> -100 0.000 -0.835 0.000 * * 0.000 0.119 0.000 * * 36) -2.024 -1.0000 | 3, 1, 2, 3/2> -100 0.000 0.835 0.000 * * 0.000 -0.119 0.000 * * 37) -2.024 1.0000 | 3, 2, 1, 1/2> -100 0.000 0.715 0.000 * * 0.000 -0.102 0.000 * * 38) -2.024 -1.0000 | 3, 2, 1, 1/2> -100 0.000 -0.715 0.000 * * 0.000 0.102 0.000 * * 39) -2.024 -1.0000 | 3, 0, 3, 5/2> -100 0.000 -1.619 0.000 * * 0.000 0.231 0.000 * * 40) -2.024 1.0000 | 3, 0, 3, 5/2> -100 0.000 1.619 0.000 * * 0.000 -0.231 0.000 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] PROTONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 20.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = ZERO Q2+1= ZERO Q2+2=-1.5E-10 ............. ............. * * 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 =-7.8E-09 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4=-4.7E-09 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS PROTONS * * * * R_RMS = 3.4541 X_RMS = 1.9942 Y_RMS = 1.9942 Z_RMS = 1.9942 * * * * R_GEO = 4.4593 X_GEO = 4.4593 Y_GEO = 4.4593 Z_GEO = 4.4593 * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) PROTONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 =-4.1E-08 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4=-3.5E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.0538 FM^(-3) BOHR RADIUS = 4.4593 FM * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 5.144197 2.626106 5.144197 0.108015 * * DRHOD = 3.633662 1.853599 3.633662 0.073536 * * DLPR_ = -1.749188 -0.896397 -1.749188 -0.043605 * * DTAU_ = 4.229355 2.185379 4.229355 0.141403 * * DSCU_ = 0.040516 0.039138 0.040516 0.037761 * * DDIV_ = 0.339868 0.185889 0.339868 0.031909 * * * * 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_ = -9071.249938 4359.611277 -4801.296721 89.658060 * * ERHOD = 6996.018768 -3945.471112 3128.809825 -78.262168 * * ELPR_ = 162.068141 -28.070002 134.680863 -0.682724 * * ETAU_ = 137.335397 107.768629 241.617480 3.486546 * * ESCU_ = 0.000000 0.000000 0.000000 0.000000 * * EDIV_ = -20.901880 -11.432143 -31.352820 -0.981203 * * ============ ============ ============ ============ * * SUM EVEN: -1796.729511 482.406650 -1327.541373 13.218512 * * * * 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 * * * ******************************************************************************* ******************************************************************************* * * * CENTER-OF-MASS COORDINATES IN FERMIS TOTAL * * * * R_CMS = 0.0000 X_CMS = 0.0000 Y_CMS = 0.0000 Z_CMS = 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] [CENT-MASS FRAME] NEUTRONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 28.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = ZERO Q2+1= ZERO Q2+2=-3.0E-10 ............. ............. * * 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 =-1.4E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4=-8.5E-09 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) [CENT-MASS FRAME] NEUTRONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 =-4.5E-08 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4=-3.8E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.0662 FM^(-3) BOHR RADIUS = 4.6556 FM * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] [CENT-MASS FRAME] PROTONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 20.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = ZERO Q2+1= ZERO Q2+2=-1.5E-10 ............. ............. * * 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 =-7.8E-09 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4=-4.7E-09 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) [CENT-MASS FRAME] PROTONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 =-4.1E-08 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4=-3.5E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.0538 FM^(-3) BOHR RADIUS = 4.4593 FM * * * ******************************************************************************* ******************************************************************************* * * * EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL * * * * ALPHA = 90.00000 BETA = 90.00000 GAMMA =179.99992 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 3.5437 X_RMS = 2.0459 Y_RMS = 2.0459 Z_RMS = 2.0459 * * * * R_GEO = 4.5748 X_GEO = 4.5748 Y_GEO = 4.5748 Z_GEO = 4.5748 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 48.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 1.4E-10 Q2+1= ZERO Q2+2=-4.5E-10 ............. ............. * * 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 =-2.2E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4=-1.3E-08 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * 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 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 =-4.4E-08 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4=-3.7E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1197 FM^(-3) BOHR RADIUS = 4.5748 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 = 48.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 1.4E-10 Q2+1= ZERO Q2+2=-4.5E-10 ............. ............. * * 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 =-2.2E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4=-1.3E-08 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * 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 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 =-4.4E-08 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4=-3.7E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1197 FM^(-3) BOHR RADIUS = 4.5748 FM * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] [INTRINSIC FRAME] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 48.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 3.2E-10 Q2+1= ZERO Q2+2= 3.5E-10 ............. ............. * * 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 =-2.2E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4=-1.3E-08 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) [INTRINSIC FRAME] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 =-4.4E-08 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4=-3.7E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1197 FM^(-3) BOHR RADIUS = 4.5748 FM * * * ******************************************************************************* ******************************************************************************* * * * LINEAR MOMENTA IN FERMI^(-1) TOTAL * * * * P_X = 0.0000 P_Y = 0.0000 P_Z = 0.0000 P_R = 0.0000 * * * * AVERAGE LINEAR MOMENTUM SQUARED IN FERMI^(-2) = 18.7222 * * * ******************************************************************************* ******************************************************************************* * * * 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.000 0.000 0.000 0.000 0.000 0.000 * * PROTONS 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 * * -------- * * TOTAL 0.000 0.000 0.000 0.000 0.000 0.000 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.000 0.000 0.000 0.000 0.000 0.000 * * PROTONS 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 * * -------- * * TOTAL 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 * * * ******************************************************************************* * * * ANGULAR MOMENTA IN THE INTRINSIC REFERENCE FRAME * * * ******************************************************************************* * * * SPIN IN X-DIRECTION SPIN IN Y-DIRECTION SPIN IN Z-DIRECTION * * ------------------- ------------------- ------------------- * * INTRI ORBIT TOTAL INTRI ORBIT TOTAL INTRI ORBIT TOTAL * * * * NEUTRONS 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 * * PROTONS 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 * * -------- * * TOTAL 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 * * * ******************************************************************************* ******************************************************************************* * * * NEUTRON CONFIGURATIONS * * ======================= * * 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 * * -------------------------------------------------------------------- * * * * 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 * * ======================= * * 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 * * -------------------------------------------------------------------- * * * * 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)= 527.402345 (PRO)= 315.602359 (TOT)= 843.004703 * * SUM EPS: (NEU)= -626.310864 (PRO)= -521.234772 (TOT)= -1147.545636 * * 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.664247 (PRO)= -13.136820 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 78.558869 (EXC)= -7.413275 (TOT)= 71.145594 * * * * 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= 254.212305 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -408.260382 * * * * TWO-BODY CENTER OF MASS CORRECTION INCLUDED BEFORE VARIATION= -8.087818 * * * * SPIN-ORB (EVE)= -32.334023 (ODD)= 0.000000 (TOT)= -32.334023 * * SKYRME: (EVE)= -1314.322862 (ODD)= 0.000000 (TOT)= -1314.322862 * * * * TOTAL: (STAB)= -9.474780 (SP)= -417.735162 (FUN)= -408.260382 * * * ******************************************************************************* ******************************************************************************* * * * RENORMALIZED MASSES IN MEV*FERMI^2 AND RATIOS W.R. TO M*A TOTAL * * * * M_X = 0.440255 M_Y = 0.440255 M_Z = 0.440255 M_R = 0.4403 * * * * F_X = 0.981228 F_Y = 0.981228 F_Z = 0.981228 F_R = 0.9812 * * * ******************************************************************************* * * * GAUSSIAN-OVERLAP MASSES IN MEV*FERMI^2 AND RATIOS W.R. TO M*A TOTAL * * * * M_X = 0.437108 M_Y = 0.437108 M_Z = 0.437108 M_R = 0.4371 * * * * F_X = 0.988291 F_Y = 0.988291 F_Z = 0.988291 F_R = 0.9883 * * * ******************************************************************************* ******************************************************************************* * * * ITERATION RESTARTED FROM THE DISC, THE FOLLOWING FILE(S) READ: * * ca48_lip.rec * * ca48_lip.fic * * * ******************************************************************************* * * * CONVERGENCE REPORT * * * ******************************************************************************* * * * ITER ENERGY STABILITY Q_2 GAMMA SPIN OMEGA ANGLE EPAIR * * * * 48 -408.260382 -9.474780 0.000 47.23 0.00 0.000 0.000 0.00 * * 49 -408.480547 0.191390 0.000 47.23 0.00 0.000 0.000 0.00 * * 50 -408.480286 0.044967 0.000 47.23 0.00 0.000 0.000 0.00 * * 51 -408.485704 0.035814 0.000 47.23 0.00 0.000 0.000 0.00 * * 52 -408.486271 0.021917 0.000 47.23 0.00 0.000 0.000 0.00 * * 53 -408.487435 0.014045 0.000 47.23 0.00 0.000 0.000 0.00 * * 54 -408.487817 0.008814 0.000 47.23 0.00 0.000 0.000 0.00 * * 55 -408.488151 0.005544 0.000 47.24 0.00 0.000 0.000 0.00 * * 56 -408.488321 0.003461 0.000 47.24 0.00 0.000 0.000 0.00 * * 57 -408.488438 0.002157 0.000 47.23 0.00 0.000 0.000 0.00 * * 58 -408.488508 0.001339 0.000 47.23 0.00 0.000 0.000 0.00 * * 59 -408.488553 0.000831 0.000 47.23 0.00 0.000 0.000 0.00 * * 60 -408.488581 0.000514 0.000 47.24 0.00 0.000 0.000 0.00 * * 61 -408.488599 0.000318 0.000 47.25 0.00 0.000 0.000 0.00 * * 62 -408.488610 0.000197 0.000 47.23 0.00 0.000 0.000 0.00 * * 63 -408.488617 0.000122 0.000 47.24 0.00 0.000 0.000 0.00 * * 64 -408.488622 0.000075 0.000 47.23 0.00 0.000 0.000 0.00 * * 65 -408.488625 0.000047 0.000 47.23 0.00 0.000 0.000 0.00 * * 66 -408.488627 0.000029 0.000 47.26 0.00 0.000 0.000 0.00 * * 67 -408.488628 0.000018 0.000 47.28 0.00 0.000 0.000 0.00 * * 68 -408.488629 0.000011 0.000 47.26 0.00 0.000 0.000 0.00 * * 69 -408.488630 6.88E-06 0.000 47.25 0.00 0.000 0.000 0.00 * * 70 -408.488630 4.27E-06 0.000 47.21 0.00 0.000 0.000 0.00 * * 71 -408.488630 2.66E-06 0.000 47.25 0.00 0.000 0.000 0.00 * * 72 -408.488630 1.65E-06 0.000 47.30 0.00 0.000 0.000 0.00 * * 73 -408.488631 1.03E-06 0.000 47.23 0.00 0.000 0.000 0.00 * * 74 -408.488631 6.42E-07 0.000 47.30 0.00 0.000 0.000 0.00 * * 75 -408.488631 4.01E-07 0.000 47.28 0.00 0.000 0.000 0.00 * * 76 -408.488631 2.50E-07 0.000 47.35 0.00 0.000 0.000 0.00 * * 77 -408.488631 1.57E-07 0.000 47.40 0.00 0.000 0.000 0.00 * * 78 -408.488631 9.84E-08 0.000 47.36 0.00 0.000 0.000 0.00 * * 79 -408.488631 6.14E-08 0.000 47.46 0.00 0.000 0.000 0.00 * * 80 -408.488631 3.86E-08 0.000 47.46 0.00 0.000 0.000 0.00 * * 81 -408.488631 2.40E-08 0.000 47.36 0.00 0.000 0.000 0.00 * * 82 -408.488631 1.51E-08 0.000 47.24 0.00 0.000 0.000 0.00 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY | N,nz,lz,OMEz>

* * * * 1) -49.329 0.9982 | 0, 0, 0, 1/2> 100 0.025 -0.499 -0.016 * * 0.025 -0.499 -0.016 * * 2) -49.329 -0.9982 | 0, 0, 0, 1/2> 100 -0.025 0.499 0.016 * * -0.025 0.499 0.016 * * 3) -35.557 0.9960 | 1, 1, 0, 1/2> -100 -0.072 0.145 -0.012 * * -0.024 0.048 -0.004 * * 4) -35.557 -0.9960 | 1, 1, 0, 1/2> -100 0.071 -0.173 0.019 * * 0.024 -0.058 0.006 * * 5) -35.557 0.9972 | 1, 0, 1, 3/2> -100 0.006 -1.142 -0.065 * * 0.002 -0.381 -0.022 * * 6) -35.557 -0.9972 | 1, 0, 1, 3/2> -100 -0.006 1.171 0.058 * * -0.002 0.390 0.019 * * 7) -32.773 0.9944 | 1, 0, 1, 1/2> -100 0.027 0.497 -0.046 * * -0.009 -0.166 0.015 * * 8) -32.773 -0.9944 | 1, 0, 1, 1/2> -100 -0.027 -0.497 0.046 * * 0.009 0.166 -0.015 * * 9) -21.901 0.2755 | 2, 1, 1, 3/2> 100 -0.604 0.230 0.527 * * -0.121 0.046 0.105 * * 10) -21.901 -0.2755 | 2, 1, 1, 3/2> 100 0.604 -0.230 -0.527 * * 0.121 -0.046 -0.105 * * 11) -21.901 -0.9881 | 2, 0, 2, 5/2> 100 0.265 0.712 0.047 * * 0.053 0.142 0.009 * * 12) -21.901 0.9881 | 2, 0, 2, 5/2> 100 -0.265 -0.716 -0.047 * * -0.053 -0.143 -0.009 * * 13) -21.901 -1.0000 | 2, 0, 2, 5/2> 100 0.003 1.612 -0.008 * * 0.001 0.322 -0.002 * * 14) -21.901 1.0000 | 2, 0, 2, 5/2> 100 -0.003 -1.608 0.008 * * -0.001 -0.322 0.002 * * 15) -17.372 0.9997 | 2, 0, 0, 1/2> 100 0.012 -0.500 -0.004 * * 0.012 -0.500 -0.004 * * 16) -17.372 -0.9997 | 2, 0, 0, 1/2> 100 -0.012 0.500 0.004 * * -0.012 0.500 0.004 * * 17) -16.433 -0.9876 | 2, 1, 1, 1/2> 100 -0.195 0.172 -0.017 * * 0.039 -0.034 0.003 * * 18) -16.433 0.9876 | 2, 1, 1, 1/2> 100 0.194 -0.187 0.017 * * -0.039 0.037 -0.003 * * 19) -16.433 0.9999 | 2, 0, 2, 3/2> 100 -0.001 1.189 0.008 * * 0.000 -0.238 -0.002 * * 20) -16.433 -0.9999 | 2, 0, 2, 3/2> 100 0.001 -1.174 -0.008 * * 0.000 0.235 0.002 * * 21) -9.476 -0.9510 | 3, 0, 3, 7/2> -100 0.094 1.109 -0.348 * * 0.013 0.158 -0.050 * * 22) -9.476 0.9510 | 3, 0, 3, 7/2> -100 -0.094 -1.109 0.348 * * -0.013 -0.158 0.050 * * 23) -9.476 -0.9994 | 3, 2, 1, 3/2> -100 -0.024 0.723 0.029 * * -0.003 0.103 0.004 * * 24) -9.476 0.9993 | 3, 2, 1, 3/2> -100 0.025 -0.707 -0.039 * * 0.004 -0.101 -0.006 * * 25) -9.476 0.9160 | 3, 0, 3, 7/2> -100 0.083 -1.490 0.502 * * 0.012 -0.213 0.072 * * 26) -9.476 -0.9160 | 3, 0, 3, 7/2> -100 -0.084 1.474 -0.492 * * -0.012 0.211 -0.070 * * 27) -9.476 -0.4541 | 3, 1, 2, 5/2> -100 1.322 -0.681 -0.198 * * 0.189 -0.097 -0.028 * * 28) -9.476 0.4541 | 3, 1, 2, 5/2> -100 -1.322 0.681 0.198 * * -0.189 0.097 0.028 * * 29) -5.273 -0.0585 | 3, 0, 1, 3/2> -100 0.078 -0.029 0.176 * * 0.026 -0.010 0.059 * * 30) -5.273 0.0587 | 3, 0, 1, 3/2> -100 -0.077 0.002 -0.194 * * -0.026 0.001 -0.065 * * 31) -5.273 -0.2336 | 3, 0, 1, 3/2> -100 -0.153 0.278 0.759 * * -0.051 0.093 0.253 * * 32) -5.273 0.2334 | 3, 0, 1, 3/2> -100 0.153 -0.251 -0.742 * * 0.051 -0.084 -0.247 * * 33) -3.315 -0.2984 | 3, 0, 1, 1/2> -100 -0.423 -0.149 -0.222 * * 0.141 0.050 0.074 * * 34) -3.315 0.2984 | 3, 0, 1, 1/2> -100 0.423 0.149 0.222 * * -0.141 -0.050 -0.074 * * 35) -1.122 0.6244 | 3, 2, 1, 1/2> -100 0.403 0.443 -0.724 * * -0.058 -0.063 0.103 * * 36) -1.122 -0.6244 | 3, 2, 1, 1/2> -100 -0.403 -0.450 0.722 * * 0.058 0.064 -0.103 * * 37) -1.122 0.1960 | 3, 0, 3, 5/2> -100 -0.162 0.321 1.272 * * 0.023 -0.046 -0.182 * * 38) -1.122 -0.1960 | 3, 0, 3, 5/2> -100 0.162 -0.314 -1.270 * * -0.023 0.045 0.181 * * 39) -1.122 0.1575 | 3, 1, 2, 3/2> -100 -0.072 -0.131 0.820 * * 0.010 0.019 -0.117 * * 40) -1.122 -0.1575 | 3, 1, 2, 3/2> -100 0.072 0.131 -0.820 * * -0.010 -0.019 0.117 * * 41) 2.524 -1.0000 | 4, 0, 4, 9/2> 100 0.001 1.833 0.005 * * 0.000 0.204 0.001 * * 42) 2.524 1.0000 | 4, 0, 4, 9/2> 100 -0.001 -1.833 -0.005 * * 0.000 -0.204 -0.001 * * 43) 2.524 -0.0902 | 4, 3, 1, 3/2> 100 0.523 -0.045 -0.092 * * 0.058 -0.005 -0.010 * * 44) 2.524 0.0902 | 4, 3, 1, 3/2> 100 -0.523 0.058 0.092 * * -0.058 0.006 0.010 * * 45) 2.524 0.9400 | 4, 0, 4, 9/2> 100 -0.276 -1.609 0.224 * * -0.031 -0.179 0.025 * * 46) 2.524 -0.9400 | 4, 0, 4, 9/2> 100 0.276 1.596 -0.224 * * 0.031 0.177 -0.025 * * 47) 2.524 0.9893 | 4, 2, 2, 5/2> 100 -0.358 -1.418 0.002 * * -0.040 -0.158 0.000 * * 48) 2.524 -0.9893 | 4, 2, 2, 5/2> 100 0.360 1.329 -0.009 * * 0.040 0.148 -0.001 * * 49) 2.524 0.9900 | 4, 1, 3, 7/2> 100 -0.257 1.884 0.053 * * -0.029 0.209 0.006 * * 50) 2.524 -0.9900 | 4, 1, 3, 7/2> 100 0.255 -1.796 -0.046 * * 0.028 -0.200 -0.005 * * 51) 3.777 0.9999 | 4, 0, 0, 1/2> 100 0.008 -0.500 0.001 * * 0.008 -0.500 0.001 * * 52) 3.777 -0.9999 | 4, 0, 0, 1/2> 100 -0.008 0.500 -0.001 * * -0.008 0.500 -0.001 * * 53) 4.075 0.9074 | 4, 1, 1, 3/2> 100 0.339 0.756 0.086 * * 0.068 0.151 0.017 * * 54) 4.075 -0.9074 | 4, 1, 1, 3/2> 100 -0.339 -0.756 -0.086 * * -0.068 -0.151 -0.017 * * 55) 4.075 0.9947 | 4, 4, 0, 1/2> 100 -0.025 -0.997 0.076 * * -0.005 -0.199 0.015 * * 56) 4.075 -0.9947 | 4, 4, 0, 1/2> 100 0.025 1.006 -0.076 * * 0.005 0.201 -0.015 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] NEUTRONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 28.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = ZERO Q2+1= ZERO Q2+2= ZERO ............. ............. * * 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 = 7.9E-09 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4= 4.7E-09 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS NEUTRONS * * * * R_RMS = 3.6167 X_RMS = 2.0881 Y_RMS = 2.0881 Z_RMS = 2.0881 * * * * R_GEO = 4.6691 X_GEO = 4.6691 Y_GEO = 4.6691 Z_GEO = 4.6691 * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) NEUTRONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 = 2.5E-08 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4= 2.1E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.0657 FM^(-3) BOHR RADIUS = 4.6691 FM * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS * * * ******************************************************************************* * * * NO) ENERGY | N,nz,lz,OMEz>

* * * * 1) -45.635 -0.9998 | 0, 0, 0, 1/2> 100 0.008 0.500 0.006 * * 0.008 0.500 0.006 * * 2) -45.635 0.9998 | 0, 0, 0, 1/2> 100 -0.008 -0.500 -0.006 * * -0.008 -0.500 -0.006 * * 3) -33.603 0.9353 | 1, 1, 0, 1/2> -100 -0.428 0.152 -0.033 * * -0.143 0.051 -0.011 * * 4) -33.603 -0.9352 | 1, 1, 0, 1/2> -100 0.422 -0.181 0.049 * * 0.141 -0.060 0.016 * * 5) -33.603 0.6896 | 1, 0, 1, 3/2> -100 0.303 -0.799 -0.293 * * 0.101 -0.266 -0.098 * * 6) -33.603 -0.6897 | 1, 0, 1, 3/2> -100 -0.297 0.828 0.276 * * -0.099 0.276 0.092 * * 7) -30.516 -0.7207 | 1, 0, 1, 1/2> -100 -0.242 -0.360 -0.248 * * 0.081 0.120 0.083 * * 8) -30.516 0.7207 | 1, 0, 1, 1/2> -100 0.242 0.360 0.248 * * -0.081 -0.120 -0.083 * * 9) -21.894 0.9656 | 2, 0, 2, 5/2> 100 0.161 -0.683 0.236 * * 0.032 -0.137 0.047 * * 10) -21.894 -0.9656 | 2, 0, 2, 5/2> 100 -0.161 0.699 -0.237 * * -0.032 0.140 -0.047 * * 11) -21.894 0.9992 | 2, 0, 2, 5/2> 100 0.019 -1.624 0.025 * * 0.004 -0.325 0.005 * * 12) -21.894 -0.9992 | 2, 0, 2, 5/2> 100 -0.019 1.608 -0.024 * * -0.004 0.322 -0.005 * * 13) -21.894 0.9912 | 2, 1, 1, 3/2> 100 -0.110 0.826 -0.006 * * -0.022 0.165 -0.001 * * 14) -21.894 -0.9912 | 2, 1, 1, 3/2> 100 0.110 -0.826 0.006 * * 0.022 -0.165 0.001 * * 15) -16.591 0.9992 | 2, 0, 0, 1/2> 100 0.020 -0.500 -0.006 * * 0.020 -0.500 -0.006 * * 16) -16.591 -0.9992 | 2, 0, 0, 1/2> 100 -0.020 0.500 0.006 * * -0.020 0.500 0.006 * * 17) -15.940 -0.9946 | 2, 1, 1, 1/2> 100 -0.111 0.173 -0.015 * * 0.022 -0.035 0.003 * * 18) -15.940 0.9946 | 2, 1, 1, 1/2> 100 0.111 -0.183 0.015 * * -0.022 0.037 -0.003 * * 19) -15.940 0.9996 | 2, 0, 2, 3/2> 100 0.014 1.184 0.000 * * -0.003 -0.237 0.000 * * 20) -15.940 -0.9996 | 2, 0, 2, 3/2> 100 -0.014 -1.174 0.000 * * 0.003 0.235 0.000 * * 21) -9.925 -0.9995 | 3, 0, 3, 7/2> -100 0.028 1.166 -0.026 * * 0.004 0.167 -0.004 * * 22) -9.925 0.9995 | 3, 0, 3, 7/2> -100 -0.028 -1.166 0.026 * * -0.004 -0.167 0.004 * * 23) -9.925 -0.9893 | 3, 2, 1, 3/2> -100 0.098 0.708 0.131 * * 0.014 0.101 0.019 * * 24) -9.925 0.9893 | 3, 2, 1, 3/2> -100 -0.098 -0.698 -0.132 * * -0.014 -0.100 -0.019 * * 25) -9.925 -0.9987 | 3, 0, 3, 7/2> -100 -0.005 1.616 -0.068 * * -0.001 0.231 -0.010 * * 26) -9.925 0.9987 | 3, 0, 3, 7/2> -100 0.005 -1.626 0.068 * * 0.001 -0.232 0.010 * * 27) -9.925 -0.5782 | 3, 1, 2, 5/2> -100 -0.649 -0.867 1.037 * * -0.093 -0.124 0.148 * * 28) -9.925 0.5782 | 3, 1, 2, 5/2> -100 0.649 0.867 -1.037 * * 0.093 0.124 -0.148 * * 29) -4.547 -0.7989 | 3, 0, 1, 3/2> -100 -0.511 -0.144 -0.105 * * -0.170 -0.048 -0.035 * * 30) -4.547 0.7989 | 3, 0, 1, 3/2> -100 0.514 0.137 0.092 * * 0.171 0.046 0.031 * * 31) -4.547 -0.8254 | 3, 0, 1, 3/2> -100 -0.205 0.974 0.300 * * -0.068 0.325 0.100 * * 32) -4.547 0.8254 | 3, 0, 1, 3/2> -100 0.203 -0.967 -0.287 * * 0.068 -0.322 -0.096 * * 33) -2.590 -0.5122 | 3, 0, 1, 1/2> -100 0.395 -0.256 -0.169 * * -0.132 0.085 0.056 * * 34) -2.590 0.5122 | 3, 0, 1, 1/2> -100 -0.395 0.256 0.169 * * 0.132 -0.085 -0.056 * * 35) -1.817 -0.6386 | 3, 2, 1, 1/2> -100 -0.034 -0.460 -0.747 * * 0.005 0.066 0.107 * * 36) -1.817 0.6386 | 3, 2, 1, 1/2> -100 0.035 0.456 0.744 * * -0.005 -0.065 -0.106 * * 37) -1.817 0.7947 | 3, 0, 3, 5/2> -100 0.116 1.289 0.772 * * -0.017 -0.184 -0.110 * * 38) -1.817 -0.7947 | 3, 0, 3, 5/2> -100 -0.117 -1.285 -0.770 * * 0.017 0.184 0.110 * * 39) -1.817 -0.6861 | 3, 1, 2, 3/2> -100 -0.091 0.572 -0.599 * * 0.013 -0.082 0.086 * * 40) -1.817 0.6861 | 3, 1, 2, 3/2> -100 0.091 -0.572 0.599 * * -0.013 0.082 -0.086 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] PROTONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 20.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = ZERO Q2+1= ZERO Q2+2= ZERO ............. ............. * * 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 = 2.8E-09 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4= 1.7E-09 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS PROTONS * * * * R_RMS = 3.4627 X_RMS = 1.9992 Y_RMS = 1.9992 Z_RMS = 1.9992 * * * * R_GEO = 4.4704 X_GEO = 4.4704 Y_GEO = 4.4704 Z_GEO = 4.4704 * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) PROTONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 = 1.5E-08 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4= 1.2E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.0534 FM^(-3) BOHR RADIUS = 4.4704 FM * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 5.101959 2.604474 5.101959 0.106989 * * DRHOD = 3.598420 1.835563 3.598420 0.072707 * * DLPR_ = -1.732322 -0.887806 -1.732322 -0.043290 * * DTAU_ = 4.170343 2.154646 4.170343 0.138949 * * DSCU_ = 0.039990 0.038642 0.039990 0.037293 * * DDIV_ = 0.336008 0.183646 0.336008 0.031284 * * * * 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_ = -8996.767510 4323.700029 -4761.874123 88.806642 * * ERHOD = 6928.166104 -3907.081906 3098.464268 -77.380071 * * ELPR_ = 160.505415 -27.800996 133.382216 -0.677797 * * ETAU_ = 135.419169 106.253080 238.246214 3.426035 * * ESCU_ = 0.000000 0.000000 0.000000 0.000000 * * EDIV_ = -20.664462 -11.294216 -30.996693 -0.961985 * * ============ ============ ============ ============ * * SUM EVEN: -1793.341283 483.775991 -1322.778118 13.212825 * * * * 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 * * * ******************************************************************************* ******************************************************************************* * * * CENTER-OF-MASS COORDINATES IN FERMIS TOTAL * * * * R_CMS = 0.0000 X_CMS = 0.0000 Y_CMS = 0.0000 Z_CMS = 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] [CENT-MASS FRAME] NEUTRONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 28.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = ZERO Q2+1= ZERO Q2+2= ZERO ............. ............. * * 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 = 7.9E-09 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4= 4.7E-09 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) [CENT-MASS FRAME] NEUTRONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 = 2.5E-08 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4= 2.1E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.0657 FM^(-3) BOHR RADIUS = 4.6691 FM * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] [CENT-MASS FRAME] PROTONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 20.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = ZERO Q2+1= ZERO Q2+2= ZERO ............. ............. * * 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 = 2.8E-09 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4= 1.7E-09 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) [CENT-MASS FRAME] PROTONS * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 = 1.5E-08 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4= 1.2E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.0534 FM^(-3) BOHR RADIUS = 4.4704 FM * * * ******************************************************************************* ******************************************************************************* * * * EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL * * * * ALPHA =270.00087 BETA = 89.86863 GAMMA = 0.00655 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 3.5533 X_RMS = 2.0515 Y_RMS = 2.0515 Z_RMS = 2.0515 * * * * R_GEO = 4.5873 X_GEO = 4.5873 Y_GEO = 4.5873 Z_GEO = 4.5873 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 48.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = ZERO Q2+1= ZERO Q2+2= ZERO ............. ............. * * 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 = 1.1E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4= 6.4E-09 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * 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 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 = 2.1E-08 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4= 1.8E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1187 FM^(-3) BOHR RADIUS = 4.5873 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 = 48.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = ZERO Q2+1= ZERO Q2+2= ZERO ............. ............. * * 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 = 1.1E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4= 6.4E-09 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * 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 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 = 2.1E-08 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4= 1.8E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1187 FM^(-3) BOHR RADIUS = 4.5873 FM * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] [INTRINSIC FRAME] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 48.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = ZERO Q2+1= ZERO Q2+2= ZERO ............. ............. * * 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 = 1.1E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4= 6.4E-09 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATIONS (FIRST-ORDER APPROXIMATION) [INTRINSIC FRAME] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = ZERO B2+1= ZERO B2+2= ZERO ............. ............. * * 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 = 2.1E-08 B4+1=-1.5E-10 B4+2= ZERO B4+3= ZERO B4+4= 1.8E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1187 FM^(-3) BOHR RADIUS = 4.5873 FM * * * ******************************************************************************* ******************************************************************************* * * * LINEAR MOMENTA IN FERMI^(-1) TOTAL * * * * P_X = 0.0000 P_Y = 0.0000 P_Z = 0.0000 P_R = 0.0000 * * * * AVERAGE LINEAR MOMENTUM SQUARED IN FERMI^(-2) = 18.8317 * * * ******************************************************************************* ******************************************************************************* * * * 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.000 0.000 0.000 0.000 0.000 0.000 * * PROTONS 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 * * -------- * * TOTAL 0.000 0.000 0.000 0.000 0.000 0.000 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.000 0.000 0.000 0.000 0.000 0.000 * * PROTONS 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 * * -------- * * TOTAL 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 * * * ******************************************************************************* * * * ANGULAR MOMENTA IN THE INTRINSIC REFERENCE FRAME * * * ******************************************************************************* * * * SPIN IN X-DIRECTION SPIN IN Y-DIRECTION SPIN IN Z-DIRECTION * * ------------------- ------------------- ------------------- * * INTRI ORBIT TOTAL INTRI ORBIT TOTAL INTRI ORBIT TOTAL * * * * NEUTRONS 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 * * PROTONS 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 * * -------- * * TOTAL 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 * * * ******************************************************************************* ******************************************************************************* * * * NEUTRON CONFIGURATIONS * * ======================= * * 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 * * -------------------------------------------------------------------- * * * * 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 * * ======================= * * 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 * * -------------------------------------------------------------------- * * * * 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)= 524.341653 (PRO)= 314.072461 (TOT)= 838.414114 * * SUM EPS: (NEU)= -614.121955 (PRO)= -515.019657 (TOT)= -1129.141612 * * 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.374194 (PRO)= -12.932373 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 78.351019 (EXC)= -7.395864 (TOT)= 70.955156 * * * * 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= 251.757016 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -408.488631 * * * * TWO-BODY CENTER OF MASS CORRECTION INCLUDED BEFORE VARIATION= -8.292608 * * * * SPIN-ORB (EVE)= -31.958678 (ODD)= 0.000000 (TOT)= -31.958678 * * SKYRME: (EVE)= -1309.565292 (ODD)= 0.000000 (TOT)= -1309.565292 * * * * TOTAL: (STAB)= 9.24E-09 (SP)= -408.488631 (FUN)= -408.488631 * * * ******************************************************************************* ******************************************************************************* * * * RENORMALIZED MASSES IN MEV*FERMI^2 AND RATIOS W.R. TO M*A TOTAL * * * * M_X = 0.440354 M_Y = 0.440354 M_Z = 0.440354 M_R = 0.4404 * * * * F_X = 0.981006 F_Y = 0.981006 F_Z = 0.981006 F_R = 0.9810 * * * ******************************************************************************* * * * GAUSSIAN-OVERLAP MASSES IN MEV*FERMI^2 AND RATIOS W.R. TO M*A TOTAL * * * * M_X = 0.437209 M_Y = 0.437209 M_Z = 0.437209 M_R = 0.4372 * * * * F_X = 0.988063 F_Y = 0.988063 F_Z = 0.988063 F_R = 0.9881 * * * ******************************************************************************* ******************************************************************************* * * * CALLING SUBROUTINE "RECORD" TO SAVE THE RECORD FILE CONTAINING RESULTS * * ca48_ren.rec * * THE MEAN FIELDS OBTAINED IN THE LAST-BUT-ONE ITERATION HAVE BEEN STORED * * * * CALLING SUBROUTINE "RECFIL" TO SAVE THE FIELDS FILE CONTAINING RESULTS * * ca48_ren.fic * * * ******************************************************************************* ******************************************************************************* * * * CALCULATION OF THE DATA SET NO. 2 HAS BEEN SUCCESSFULLY COMPLETED * * * ******************************************************************************* ******************************************************************************* * * * PROCESSOR EXECUTION TIMES IN SUBROUTINES * * INCLUDING THEIR COMPLETE DOWN-CALLING TREES * * * ******************************************************************************* * * * 18743 => HFODD 16774 => LINAVR 7408 => RENMAS 747 => DENSHF * * * * 602 => DIAMAT 581 => DIASIZ 151 => INTEGH 137 => COUMAT * * * * 135 => INTCOU 112 => MOMETS 78 => DENMAC 77 => INTMUL * * * * 63 => INTSOR 29 => INTMAS 24 => DIASIG 19 => BEGINT * * * * 16 => SPAVER 13 => INTCEN 10 => POWALL 10 => AVOBSE * * * * 6 => AVANGY 5 => SKFILD 5 => SHIWAV 3 => SPIMOM * * * * 3 => RECORD 2 => MAGMOM 2 => ROTAVR 2 => INTKIN * * * * 2 => NILABS 2 => BEGINC 1 => INTSPI 0 => INTORB * * * * 0 => NILSON 0 => AVPARI 0 => RECOUL 0 => SHICOE * * * * 0 => PAXMOM 0 => AVSIMP 0 => MOMSIF 0 => MOMVMU * * * * 0 => SHIMOM 0 => SHISIF 0 => SHIMAG * * * ******************************************************************************* * * * NUMBERS OF CALLS TO SUBROUTINES * * * ******************************************************************************* * * * 4480 => INTMUL 2434 => BEGINT 562 => INTSOR 478 => INTKIN * * * * 374 => INTMAS 374 => INTCEN 348 => DIAMAT 312 => DENSHF * * * * 312 => DENMAC 312 => LINAVR 256 => PAXMOM 240 => MOMETS * * * * 168 => AVPARI 168 => MOMSIF 168 => MOMVMU 168 => MAGMOM * * * * 168 => SPIMOM 168 => ROTAVR 168 => SHICOE 164 => INTEGH * * * * 157 => COUMAT 157 => BEGINC 157 => INTCOU 144 => SHIWAV * * * * 96 => DIASIG 85 => RECORD 84 => SKFILD 84 => SHIMOM * * * * 84 => SHISIF 84 => SHIMAG 72 => SPAVER 72 => DIASIZ * * * * 72 => INTSPI 36 => INTORB 36 => RENMAS 8 => AVANGY * * * * 8 => AVOBSE 8 => AVSIMP 8 => NILABS 3 => RECOUL * * * * 2 => POWALL 2 => NILSON 1 => HFODD * * * ******************************************************************************* * * * EXECUTION ENDS ON 2011.02.24 AT 06:16:26.514 * * * ******************************************************************************* ALL DONE