******************************************************************************* * * * 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.23 AT 14:01:58.614 * * * ******************************************************************************* ******************************************************************************* * * * CLASSICAL NUCLEAR SURFACE DEFINED FOR: N = 20 Z = 20 * * * ******************************************************************************* * * * SPHERICAL SHAPE REQUESTED FOR THE SURFACE * * * ******************************************************************************* * * * HOMEGA= 14.3862 FCHOM0= 1.2000 * * * * OSCILLATOR FREQUENCIES: HBAROX= 14.3862 HBAROY= 14.3862 HBAROZ= 14.3862 * * * * MOMENTS OF INERTIA: XMOMFC= 6.8244 YMOMFC= 6.8244 ZMOMFC= 6.8244 * * * * CENTERS OF MASS: CMSXFC= 0.0000 CMSYFC= 0.0000 CMSZFC= 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * SKYRME FORCE DEFINITION * * * ******************************************************************************* * * * PARAMETER SET SV : T0= -1248.29 T1= 970.56 T2= 107.22 T3= 0.00 * * * * POWER=1/1 W=150.000 X0= -0.17000 X1= 0.00000 X2= 0.00000 X3= 1.00000 * * * * ETA_J=0 ETA_W=0 ETACM=0 ETA_M=1 HBM= 20.73000 (SKYRME-FORCE SPIN-ORBIT) * * * ******************************************************************************* ******************************************************************************* * * * NON-STANDARD TREATMENT OF THE SKYRME FORCE: ETA_J=1 ETA_W=0 ETACM=0 ETA_M=1 * * * * THIS INCLUDES: * * * * NON-STANDARD TREATMENT OF THE TENSOR TERMS: THESE TERMS HAVE BEEN INCLUDED * * * ******************************************************************************* ******************************************************************************* * * * COUPLING CONSTANTS DEFINING THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * CRHO_ = -571.092675 205.967850 -468.108750 102.983925 * * CRHOD = 0.000000 0.000000 0.000000 0.000000 * * CLPR_ = -175.278750 94.340625 -128.108438 47.170312 * * CTAU_ = 269.445000 -107.917500 215.486250 -53.958750 * * CSCU_ = 0.000000 107.917500 53.958750 53.958750 * * CDIV_ = -75.000000 -75.000000 -112.500000 -37.500000 * * * * CSPI_ = 53.052325 312.072500 209.088575 156.036250 * * CSPID = 0.000000 0.000000 0.000000 0.000000 * * CLPS_ = 0.000000 94.340625 47.170312 47.170312 * * CCUR_ = -269.445000 107.917500 -215.486250 53.958750 * * CKIS_ = 0.000000 -107.917500 -53.958750 -53.958750 * * CROT_ = -75.000000 -75.000000 -112.500000 -37.500000 * * * ******************************************************************************* ******************************************************************************* * * * PHYSICAL CONSTANTS: H_BARC=197.32891000 HBCOE2=137.03602000 * * * * XMASSN=938.90590000 XMASSP=938.27231000 * * * * HBMASS= 20.73000000 HBMRPA= 20.21175000 * * * * ECHAR2= 1.43997841 COULEX= -1.06350868 * * * ******************************************************************************* * * * OSCILLATOR LENGTHS: X= 1.6978805 Y= 1.6978805 Z= 1.6978805 * * * * OSCILLATOR CONSTANTS: X= 0.5889696 Y= 0.5889696 Z= 0.5889696 * * * * OSCILLATOR FREQUENCIES: X=14.3861673 Y=14.3861673 Z=14.3861673 * * * ******************************************************************************* * * * BASIS CUT-OFF CONTROL PARAMETERS: NXMAXX = 10 NYMAXX = 10 NZMAXX = 10 * * * * OPTIMUM NUMBERS OF GAUSS POINTS: NXHERM = 22 NYHERM = 22 NZHERM = 22 * * * * NLIMIT= 286 LDBASE= 286 MCOUNT= 1331 * * * * ENECUT= 800.0000 ELIMIT= 165.4409 * * * ******************************************************************************* ******************************************************************************* * * * SHAPE OF THE OSCILLATOR-BASIS DIAMOND (THREE CONSECUTIVE PROJECTIONS) * * * ******************************************************************************* * * * NX ===>>> 0 1 2 3 4 5 6 7 8 9 10 * * MAX.NY => 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------- * * NY= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NY= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NY= 2 ( 8) | 8 7 6 5 4 3 2 1 0 * * NY= 3 ( 7) | 7 6 5 4 3 2 1 0 * * NY= 4 ( 6) | 6 5 4 3 2 1 0 * * NY= 5 ( 5) | 5 4 3 2 1 0 * * NY= 6 ( 4) | 4 3 2 1 0 * * NY= 7 ( 3) | 3 2 1 0 * * NY= 8 ( 2) | 2 1 0 * * NY= 9 ( 1) | 1 0 * * NY=10 ( 0) | 0 * * * * * * NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 * * MAX.NY => 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------- * * NY= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NY= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NY= 2 ( 8) | 8 7 6 5 4 3 2 1 0 * * NY= 3 ( 7) | 7 6 5 4 3 2 1 0 * * NY= 4 ( 6) | 6 5 4 3 2 1 0 * * NY= 5 ( 5) | 5 4 3 2 1 0 * * NY= 6 ( 4) | 4 3 2 1 0 * * NY= 7 ( 3) | 3 2 1 0 * * NY= 8 ( 2) | 2 1 0 * * NY= 9 ( 1) | 1 0 * * NY=10 ( 0) | 0 * * * * * * NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 * * MAX.NX => 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------- * * NX= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NX= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NX= 2 ( 8) | 8 7 6 5 4 3 2 1 0 * * NX= 3 ( 7) | 7 6 5 4 3 2 1 0 * * NX= 4 ( 6) | 6 5 4 3 2 1 0 * * NX= 5 ( 5) | 5 4 3 2 1 0 * * NX= 6 ( 4) | 4 3 2 1 0 * * NX= 7 ( 3) | 3 2 1 0 * * NX= 8 ( 2) | 2 1 0 * * NX= 9 ( 1) | 1 0 * * NX=10 ( 0) | 0 * * * ******************************************************************************* ******************************************************************************* * * * AUTOMATED HANDLING OF THE COULOMB FILE * * * ******************************************************************************* * * * THE COULOMB FILE WRITTEN ON THE DISC, FILE NAME (FILCOU): * * ca40_iso.cou * * * ******************************************************************************* * * * NUCLIDE: N = 20 Z = 20 NUMBER OF ITERATIONS = 100 CONTINUATION? = 0 * * * ******************************************************************************* * * * ITERATIONS WILL STOP WHEN THE STABILITY BECOMES LOWER THAN EPSITE= 1.0E-09 * * FOR 3 CONSECUTIVE ITERATIONS * * * ******************************************************************************* * * * SLOW-DOWN PARAMETERS: * * TIME-EVEN = 0.50 TIME-ODD = 0.50 PAIRING = 0.50 LIPKIN-NOGAMI = 0.50 * * * ******************************************************************************* * * * BROYDEN METHOD IS: ON * * * * TRIGGERED ONLY WHEN STABILITY IS LOWER THAN : 10000.000 MEV * * INITIAL SLOWING FACTOR (BEFORE TRIGGER) : 0.50 (=SLOWEV) * * BROYDEN SLOWING FACTOR (AFTER TRIGGER) : 0.70 (=1-ALPHAM) * * NUMBER OF ITERATIONS RETAINED IN MEMORY : 7 * * * ******************************************************************************* * * * DIRECT COULOMB TERMS SOLVED BY GREEN FUNCTION AND FOURIER TRANSFORM METHOD * * COULOMB PARAMETERS: NUMCOU = 80 NUMETA = 79 BOUCOU = 20.0 FURMAX =0.25 * * * ******************************************************************************* * * * EXCHANGE COULOMB TERMS TREATED BY THE SLATER APPROXIMATION * * * ******************************************************************************* * * * MAXIMUM NUMBERS OF MULTIPOLE MOMENTS USED FOR THE CONSTRAINTS, NMUCON=2 * * FOR THE COULOMB FIELD, NMUCOU=4 * * FOR THE OUTPUT INFO, NMUPRI=4 * * * ******************************************************************************* * * * MULTIPOLE CONSTRAINTS: LAMBDA= 2 MIU= 0 MOMENT= 0.200 STIFFNESS= 0.250 * * LAMBDA= 2 MIU= 2 MOMENT= 0.000 STIFFNESS= 0.250 * * * ******************************************************************************* * * * PRINTING THE RESULTS FOR THE FOLLOWING ITERATIONS: THE FIRST: NO * * 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: 20 20 * * 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: 20 0 0 * * PROTONS : 20 0 0 * * * ******************************************************************************* ******************************************************************************* * * * PARAMETERS OF THE NILSSON MODEL NEUTRONS * * * ******************************************************************************* * * * ESCALE =58.352 KAPPA = 0.080 MIU = 0.000 * * * * HOMEGA =14.386 CNILSS =-2.302 DNILSS = 0.000 * * * * OSCILLATOR FREQUENCIES: X=14.386167272 Y=14.386167272 Z=14.386167272 * * * ******************************************************************************* ******************************************************************************* * * * PARAMETERS OF THE NILSSON MODEL PROTONS * * * ******************************************************************************* * * * ESCALE =58.352 KAPPA = 0.080 MIU = 0.000 * * * * HOMEGA =14.386 CNILSS =-2.302 DNILSS = 0.000 * * * * OSCILLATOR FREQUENCIES: X=14.386167272 Y=14.386167272 Z=14.386167272 * * * ******************************************************************************* ******************************************************************************* * * * CONVERGENCE REPORT * * * ******************************************************************************* * * * ITER ENERGY STABILITY Q_2 GAMMA SPIN OMEGA ANGLE EPAIR * * * * 0 -215.815784 329.853764 0.000 9.03 0.00 0.000 0.000 0.00 * * 1 -297.791858 397.282073 0.007 0.00 0.00 0.000 0.000 0.00 * * 2 -313.977917 357.094710 0.005 0.00 0.00 0.000 0.000 0.00 * * 3 -339.544980 55.764975 0.003 0.00 0.00 0.000 0.000 0.00 * * 4 -340.288984 62.133242 0.003 0.00 0.00 0.000 0.000 0.00 * * 5 -341.861369 28.881163 0.003 0.00 0.00 0.000 0.000 0.00 * * 6 -342.032158 25.207110 0.003 0.00 0.00 0.000 0.000 0.00 * * 7 -342.489394 -6.723878 0.003 0.00 0.00 0.000 0.000 0.00 * * 8 -342.498474 -3.807345 0.003 0.00 0.00 0.000 0.000 0.00 * * 9 -342.504544 -1.117880 0.003 0.00 0.00 0.000 0.000 0.00 * * 10 -342.505235 -1.408908 0.003 0.00 0.00 0.000 0.000 0.00 * * 11 -342.506082 0.141367 0.003 0.00 0.00 0.000 0.000 0.00 * * 12 -342.506161 -0.070233 0.003 0.00 0.00 0.000 0.000 0.00 * * 13 -342.506167 -0.002417 0.003 0.00 0.00 0.000 0.000 0.00 * * 14 -342.506168 0.009814 0.003 0.00 0.00 0.000 0.000 0.00 * * 15 -342.506168 0.001018 0.003 0.00 0.00 0.000 0.000 0.00 * * 16 -342.506168 0.000533 0.003 0.00 0.00 0.000 0.000 0.00 * * 17 -342.506169 0.001275 0.003 0.00 0.00 0.000 0.000 0.00 * * 18 -342.506169 -0.000459 0.003 0.00 0.00 0.000 0.000 0.00 * * 19 -342.506169 0.000310 0.003 0.00 0.00 0.000 0.000 0.00 * * 20 -342.506169 -0.000526 0.003 0.00 0.00 0.000 0.000 0.00 * * 21 -342.506169 -0.000100 0.003 0.00 0.00 0.000 0.000 0.00 * * 22 -342.506169 0.000055 0.003 0.00 0.00 0.000 0.000 0.00 * * 23 -342.506169 -0.000053 0.003 0.00 0.00 0.000 0.000 0.00 * * 24 -342.506169 -3.04E-06 0.003 0.00 0.00 0.000 0.000 0.00 * * 25 -342.506169 -0.000010 0.003 0.00 0.00 0.000 0.000 0.00 * * 26 -342.506169 -1.96E-06 0.003 0.00 0.00 0.000 0.000 0.00 * * 27 -342.506169 3.45E-06 0.003 0.00 0.00 0.000 0.000 0.00 * * 28 -342.506169 -4.25E-06 0.003 0.00 0.00 0.000 0.000 0.00 * * 29 -342.506169 4.76E-07 0.003 0.00 0.00 0.000 0.000 0.00 * * 30 -342.506169 -7.15E-07 0.003 0.00 0.00 0.000 0.000 0.00 * * 31 -342.506169 2.30E-07 0.003 0.00 0.00 0.000 0.000 0.00 * * 32 -342.506169 2.80E-07 0.003 0.00 0.00 0.000 0.000 0.00 * * 33 -342.506169 3.74E-08 0.003 0.00 0.00 0.000 0.000 0.00 * * 34 -342.506169 7.56E-09 0.003 0.00 0.00 0.000 0.000 0.00 * * 35 -342.506169 -3.80E-08 0.003 0.00 0.00 0.000 0.000 0.00 * * 36 -342.506169 1.33E-08 0.003 0.00 0.00 0.000 0.000 0.00 * * 37 -342.506169 6.16E-09 0.003 0.00 0.00 0.000 0.000 0.00 * * 38 -342.506169 3.28E-09 0.003 0.00 0.00 0.000 0.000 0.00 * * 39 -342.506169 2.10E-09 0.003 0.00 0.00 0.000 0.000 0.00 * * 40 -342.506169 8.26E-10 0.003 0.00 0.00 0.000 0.000 0.00 * * 41 -342.506169 5.86E-10 0.003 0.00 0.00 0.000 0.000 0.00 * * 42 -342.506169 -2.88E-11 0.003 0.00 0.00 0.000 0.000 0.00 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY | N,nz,lz,OMEz>

* * * * 9) -27.609 -0.8803 | 2, 2, 0, 1/2> 100 0.701 1.320 -0.040 * * 0.140 0.264 -0.008 * * 10) -27.609 0.8803 | 2, 2, 0, 1/2> 100 -0.701 -1.320 0.040 * * -0.140 -0.264 0.008 * * 11) -27.601 0.9180 | 2, 1, 1, 3/2> 100 0.000 0.000 -0.286 * * 0.000 0.000 -0.057 * * 12) -27.601 -0.9180 | 2, 1, 1, 3/2> 100 0.000 0.000 0.286 * * 0.000 0.000 0.057 * * 13) -27.585 -0.8488 | 2, 0, 2, 5/2> 100 0.000 0.000 -1.016 * * 0.000 0.000 -0.203 * * 14) -27.585 0.8488 | 2, 0, 2, 5/2> 100 0.000 0.000 1.016 * * 0.000 0.000 0.203 * * 15) -19.134 -0.9997 | 2, 0, 0, 1/2> 100 0.013 0.494 -0.002 * * 0.013 0.499 -0.002 * * 16) -19.134 0.9997 | 2, 0, 0, 1/2> 100 -0.013 -0.494 0.002 * * -0.013 -0.499 0.002 * * 17) -18.926 0.9986 | 2, 1, 1, 1/2> 100 -0.050 0.993 0.008 * * 0.010 -0.201 -0.002 * * 18) -18.926 -0.9986 | 2, 1, 1, 1/2> 100 0.050 -0.993 -0.008 * * -0.010 0.201 0.002 * * 19) -18.906 0.9844 | 2, 0, 2, 3/2> 100 0.000 0.000 -0.146 * * 0.000 0.000 0.029 * * 20) -18.906 -0.9844 | 2, 0, 2, 3/2> 100 0.000 0.000 0.146 * * 0.000 0.000 -0.029 * * 21) -7.399 -0.9998 | 3, 2, 1, 1/2> -100 -0.025 2.000 -0.006 * * -0.004 0.286 -0.001 * * 22) -7.399 0.9998 | 3, 2, 1, 1/2> -100 0.025 -2.000 0.006 * * 0.004 -0.286 0.001 * * 23) -7.393 -0.9678 | 3, 2, 1, 3/2> -100 0.000 0.000 -0.171 * * 0.000 0.000 -0.024 * * 24) -7.393 0.9678 | 3, 2, 1, 3/2> -100 0.000 0.000 0.171 * * 0.000 0.000 0.024 * * 25) -7.383 -0.9777 | 3, 1, 2, 5/2> -100 0.000 0.000 0.086 * * 0.000 0.000 0.012 * * 26) -7.383 0.9777 | 3, 1, 2, 5/2> -100 0.000 0.000 -0.086 * * 0.000 0.000 -0.012 * * 27) -7.367 0.5746 | 3, 0, 3, 7/2> -100 0.000 0.000 -2.788 * * 0.000 0.000 -0.398 * * 28) -7.367 -0.5746 | 3, 0, 3, 7/2> -100 0.000 0.000 2.788 * * 0.000 0.000 0.398 * * 29) -1.400 0.9915 | 3, 3, 0, 1/2> -100 -0.054 -0.991 0.059 * * -0.018 -0.334 0.021 * * 30) -1.400 -0.9915 | 3, 3, 0, 1/2> -100 0.054 0.991 -0.059 * * 0.018 0.334 -0.021 * * 31) -1.375 -0.2763 | 3, 0, 1, 3/2> -100 0.000 0.000 0.872 * * 0.000 0.000 0.291 * * 32) -1.375 0.2763 | 3, 0, 1, 3/2> -100 0.000 0.000 -0.872 * * 0.000 0.000 -0.291 * * 33) 0.780 0.1731 | 3, 0, 1, 1/2> -100 -0.028 0.087 0.492 * * 0.009 -0.028 -0.171 * * 34) 0.780 -0.1731 | 3, 0, 1, 1/2> -100 0.028 -0.087 -0.492 * * -0.009 0.028 0.171 * * 35) 2.183 0.0124 | 3, 2, 1, 1/2> -100 0.922 0.019 -0.394 * * -0.132 -0.003 0.057 * * 36) 2.183 -0.0124 | 3, 2, 1, 1/2> -100 -0.922 -0.019 0.394 * * 0.132 0.003 -0.057 * * 37) 2.193 0.8828 | 3, 1, 2, 3/2> -100 0.000 0.000 -0.649 * * 0.000 0.000 0.093 * * 38) 2.193 -0.8828 | 3, 1, 2, 3/2> -100 0.000 0.000 0.649 * * 0.000 0.000 -0.093 * * 39) 2.214 0.5453 | 3, 0, 3, 5/2> -100 0.000 0.000 -0.697 * * 0.000 0.000 0.100 * * 40) 2.214 -0.5453 | 3, 0, 3, 5/2> -100 0.000 0.000 0.697 * * 0.000 0.000 -0.100 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] NEUTRONS * * * * 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 = 0.0013 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 = 4.1E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4= 7.8E-10 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS NEUTRONS * * * * R_RMS = 3.3240 X_RMS = 1.9188 Y_RMS = 1.9188 Z_RMS = 1.9197 * * * * R_GEO = 4.2912 X_GEO = 4.2906 Y_GEO = 4.2906 Z_GEO = 4.2925 * * * ******************************************************************************* ******************************************************************************* * * * 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 = 4.6E-04 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.6E-07 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4= 6.8E-09 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.0604 FM^(-3) BOHR RADIUS = 4.2912 FM * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS * * * ******************************************************************************* * * * NO) ENERGY | N,nz,lz,OMEz>

* * * * 9) -20.268 0.9991 | 2, 2, 0, 1/2> 100 0.059 -1.499 -0.007 * * 0.012 -0.300 -0.001 * * 10) -20.268 -0.9991 | 2, 2, 0, 1/2> 100 -0.059 1.499 0.007 * * -0.012 0.300 0.001 * * 11) -20.260 -0.1244 | 2, 1, 1, 3/2> 100 0.000 0.000 -1.380 * * 0.000 0.000 -0.277 * * 12) -20.260 0.1244 | 2, 1, 1, 3/2> 100 0.000 0.000 1.380 * * 0.000 0.000 0.277 * * 13) -20.244 0.9812 | 2, 0, 2, 5/2> 100 0.000 0.000 -0.192 * * 0.000 0.000 -0.038 * * 14) -20.244 -0.9812 | 2, 0, 2, 5/2> 100 0.000 0.000 0.192 * * 0.000 0.000 0.038 * * 15) -11.978 -1.0000 | 2, 0, 0, 1/2> 100 0.004 0.490 0.002 * * 0.004 0.498 0.002 * * 16) -11.978 1.0000 | 2, 0, 0, 1/2> 100 -0.004 -0.490 -0.002 * * -0.004 -0.498 -0.002 * * 17) -11.824 -0.9999 | 2, 1, 1, 1/2> 100 0.003 -0.989 -0.005 * * -0.001 0.202 0.001 * * 18) -11.824 0.9999 | 2, 1, 1, 1/2> 100 -0.003 0.989 0.005 * * 0.001 -0.202 -0.001 * * 19) -11.804 0.9939 | 2, 0, 2, 3/2> 100 0.000 0.000 -0.153 * * 0.000 0.000 0.031 * * 20) -11.804 -0.9939 | 2, 0, 2, 3/2> 100 0.000 0.000 0.153 * * 0.000 0.000 -0.031 * * 21) -0.805 -0.7104 | 3, 2, 1, 1/2> -100 0.161 1.421 -0.350 * * 0.023 0.203 -0.050 * * 22) -0.805 0.7104 | 3, 2, 1, 1/2> -100 -0.161 -1.421 0.350 * * -0.023 -0.203 0.050 * * 23) -0.800 -0.1755 | 3, 2, 1, 3/2> -100 0.000 0.000 1.401 * * 0.000 0.000 0.201 * * 24) -0.800 0.1755 | 3, 2, 1, 3/2> -100 0.000 0.000 -1.401 * * 0.000 0.000 -0.201 * * 25) -0.789 0.9968 | 3, 1, 2, 5/2> -100 0.000 0.000 -0.119 * * 0.000 0.000 -0.017 * * 26) -0.789 -0.9968 | 3, 1, 2, 5/2> -100 0.000 0.000 0.119 * * 0.000 0.000 0.017 * * 27) -0.773 0.7043 | 3, 0, 3, 7/2> -100 0.000 0.000 0.736 * * 0.000 0.000 0.105 * * 28) -0.773 -0.7043 | 3, 0, 3, 7/2> -100 0.000 0.000 -0.736 * * 0.000 0.000 -0.105 * * 29) 4.740 -0.4278 | 3, 3, 0, 1/2> -100 0.178 0.428 -0.443 * * 0.060 0.144 -0.155 * * 30) 4.740 0.4278 | 3, 3, 0, 1/2> -100 -0.178 -0.428 0.443 * * -0.060 -0.144 0.155 * * 31) 4.764 -0.9736 | 3, 0, 1, 3/2> -100 0.000 0.000 0.146 * * 0.000 0.000 0.049 * * 32) 4.764 0.9736 | 3, 0, 1, 3/2> -100 0.000 0.000 -0.146 * * 0.000 0.000 -0.049 * * 33) 6.706 0.9442 | 3, 0, 1, 1/2> -100 -0.085 0.472 -0.141 * * 0.028 -0.153 0.049 * * 34) 6.706 -0.9442 | 3, 0, 1, 1/2> -100 0.085 -0.472 0.141 * * -0.028 0.153 -0.049 * * 35) 8.265 0.0283 | 3, 2, 1, 1/2> -100 -0.143 0.042 -0.498 * * 0.020 -0.006 0.071 * * 36) 8.265 -0.0283 | 3, 2, 1, 1/2> -100 0.143 -0.042 0.498 * * -0.020 0.006 -0.071 * * 37) 8.275 -0.5886 | 3, 1, 2, 3/2> -100 0.000 0.000 1.210 * * 0.000 0.000 -0.173 * * 38) 8.275 0.5886 | 3, 1, 2, 3/2> -100 0.000 0.000 -1.210 * * 0.000 0.000 0.173 * * 39) 8.296 -0.1704 | 3, 0, 3, 5/2> -100 0.000 0.000 -2.145 * * 0.000 0.000 0.307 * * 40) 8.296 0.1704 | 3, 0, 3, 5/2> -100 0.000 0.000 2.145 * * 0.000 0.000 -0.307 * * * ******************************************************************************* ******************************************************************************* * * * 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 = 0.0014 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.9E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4=-1.2E-08 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= ZERO * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS PROTONS * * * * R_RMS = 3.3696 X_RMS = 1.9451 Y_RMS = 1.9451 Z_RMS = 1.9460 * * * * R_GEO = 4.3501 X_GEO = 4.3495 Y_GEO = 4.3495 Z_GEO = 4.3515 * * * ******************************************************************************* ******************************************************************************* * * * 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 = 4.9E-04 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.7E-07 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4=-9.6E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.0580 FM^(-3) BOHR RADIUS = 4.3501 FM * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 4.245202 2.123090 4.245202 0.000977 * * DRHOD = 0.556645 0.278389 0.556645 0.000133 * * DLPR_ = -1.546566 -0.773821 -1.546566 -0.001076 * * DTAU_ = 3.342010 1.671568 3.342010 0.001126 * * DSCU_ = 0.000248 0.000124 0.000248 0.000001 * * DDIV_ = 0.013944 0.006959 0.013944 -0.000027 * * * * 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_ = -2424.403968 437.288222 -1987.216367 0.100621 * * ERHOD = 0.000000 0.000000 0.000000 0.000000 * * ELPR_ = 271.080181 -73.002763 198.128172 -0.050755 * * ETAU_ = 900.487818 -180.391415 720.157149 -0.060746 * * ESCU_ = 0.000000 0.013416 0.013387 0.000030 * * EDIV_ = -1.045781 -0.521895 -1.568671 0.000995 * * ============ ============ ============ ============ * * SUM EVEN: -1253.881750 183.385565 -1070.486330 -0.009855 * * * * 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 = 20.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 0.0013 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 = 4.1E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4= 7.8E-10 * * 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 = 4.6E-04 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.6E-07 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4= 6.8E-09 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.0604 FM^(-3) BOHR RADIUS = 4.2912 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 = 0.0014 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.9E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4=-1.2E-08 * * 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 = 4.9E-04 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.7E-07 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4=-9.6E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.0580 FM^(-3) BOHR RADIUS = 4.3501 FM * * * ******************************************************************************* ******************************************************************************* * * * EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL * * * * ALPHA =129.14888 BETA = 0.00000 GAMMA = 0.00000 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 3.3469 X_RMS = 1.9320 Y_RMS = 1.9320 Z_RMS = 1.9329 * * * * R_GEO = 4.3208 X_GEO = 4.3201 Y_GEO = 4.3201 Z_GEO = 4.3221 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 40.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 0.0027 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.0E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4=-1.1E-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 = 4.7E-04 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-07 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4=-4.6E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1184 FM^(-3) BOHR RADIUS = 4.3208 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 = 40.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 0.0027 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.0E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4=-1.1E-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 = 4.7E-04 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-07 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4=-4.6E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= ZERO * * * ******************************************************************************* * * * BOHR DENSITY = 0.1184 FM^(-3) BOHR RADIUS = 4.3208 FM * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] [INTRINSIC FRAME] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 40.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 0.0027 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.0E-08 Q4+1= ZERO Q4+2= ZERO Q4+3= ZERO Q4+4= 9.9E-09 * * Q4-1= ZERO Q4-2= ZERO Q4-3= ZERO Q4-4= 4.3E-09 * * * ******************************************************************************* ******************************************************************************* * * * 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 = 4.7E-04 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-07 B4+1= ZERO B4+2= ZERO B4+3= ZERO B4+4= 4.2E-08 * * B4-1= ZERO B4-2= ZERO B4-3= ZERO B4-4= 1.8E-08 * * * ******************************************************************************* * * * BOHR DENSITY = 0.1184 FM^(-3) BOHR RADIUS = 4.3208 FM * * * ******************************************************************************* ******************************************************************************* * * * LINEAR MOMENTA IN FERMI^(-1) TOTAL * * * * P_X = 0.0000 P_Y = 0.0000 P_Z = 0.0000 P_R = 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * ANGULAR MOMENTA IN THE THREE CARTESIAN DIRECTIONS * * * ******************************************************************************* * * * SPIN IN X-DIRECTION SPIN IN Y-DIRECTION SPIN IN Z-DIRECTION * * ------------------- ------------------- ------------------- * * INTRI ORBIT TOTAL INTRI ORBIT TOTAL INTRI ORBIT TOTAL * * * * NEUTRONS 0.000 0.000 0.000 0.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 * * ======================= * * 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 * * * ******************************************************************************* * * * 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)= 331.803322 (PRO)= 323.081109 (TOT)= 654.884431 * * SUM EPS: (NEU)= -743.109718 (PRO)= -591.723990 (TOT)= -1334.833708 * * PAIRING: (NEU)= 0.000000 (PRO)= 0.000000 (TOT)= 0.000000 * * P-REARR: (NEU)= 0.000000 (PRO)= 0.000000 (TOT)= 0.000000 * * * * PAIRGAP: (NEU)= 0.000000 (PRO)= 0.000000 THIS SPACE IS EMPTY * * E-FERMI: (NEU)= -13.152236 (PRO)= -6.304733 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 80.700573 (EXC)= -7.594987 (TOT)= 73.105585 * * * * CONSTR. (MULT)= 0.009734 SLOPE= 0.098659 CORR.= -0.000132 * * 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= 0.000000 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -342.496435 * * * * SPIN-ORB (EVE)= -1.567675 (ODD)= 0.000000 (TOT)= -1.567675 * * SKYRME: (EVE)= -1070.496185 (ODD)= 0.000000 (TOT)= -1070.496185 * * * * TOTAL: (STAB)= -1.14E-10 (SP)= -342.506169 (FUN)= -342.506169 * * * ******************************************************************************* ******************************************************************************* * * * CALLING SUBROUTINE "RECORD" TO SAVE THE RECORD FILE CONTAINING RESULTS * * ca40_iso.rec * * THE MEAN FIELDS OBTAINED IN THE LAST-BUT-ONE ITERATION HAVE BEEN STORED * * * * CALLING SUBROUTINE "RECFIL" TO SAVE THE FIELDS FILE CONTAINING RESULTS * * ca40_iso.fil * * * ******************************************************************************* ******************************************************************************* * * * CALCULATION OF THE DATA SET NO. 1 HAS BEEN SUCCESSFULLY COMPLETED * * * ******************************************************************************* ******************************************************************************* * * * CLASSICAL NUCLEAR SURFACE DEFINED FOR: N = 20 Z = 20 * * * ******************************************************************************* * * * SPHERICAL SHAPE REQUESTED FOR THE SURFACE * * * ******************************************************************************* * * * HOMEGA= 14.3862 FCHOM0= 1.2000 * * * * OSCILLATOR FREQUENCIES: HBAROX= 14.3862 HBAROY= 14.3862 HBAROZ= 14.3862 * * * * MOMENTS OF INERTIA: XMOMFC= 6.8244 YMOMFC= 6.8244 ZMOMFC= 6.8244 * * * * CENTERS OF MASS: CMSXFC= 0.0000 CMSYFC= 0.0000 CMSZFC= 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * SKYRME FORCE DEFINITION * * * ******************************************************************************* * * * PARAMETER SET SV : T0= -1248.29 T1= 970.56 T2= 107.22 T3= 0.00 * * * * POWER=1/1 W=150.000 X0= -0.17000 X1= 0.00000 X2= 0.00000 X3= 1.00000 * * * * ETA_J=0 ETA_W=0 ETACM=0 ETA_M=1 HBM= 20.73000 (SKYRME-FORCE SPIN-ORBIT) * * * ******************************************************************************* ******************************************************************************* * * * NON-STANDARD TREATMENT OF THE SKYRME FORCE: ETA_J=1 ETA_W=0 ETACM=0 ETA_M=1 * * * * THIS INCLUDES: * * * * NON-STANDARD TREATMENT OF THE TENSOR TERMS: THESE TERMS HAVE BEEN INCLUDED * * * ******************************************************************************* ******************************************************************************* * * * COUPLING CONSTANTS DEFINING THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * CRHO_ = -571.092675 205.967850 -468.108750 102.983925 * * CRHOD = 0.000000 0.000000 0.000000 0.000000 * * CLPR_ = -175.278750 94.340625 -128.108438 47.170312 * * CTAU_ = 269.445000 -107.917500 215.486250 -53.958750 * * CSCU_ = 0.000000 107.917500 53.958750 53.958750 * * CDIV_ = -75.000000 -75.000000 -112.500000 -37.500000 * * * * CSPI_ = 53.052325 312.072500 209.088575 156.036250 * * CSPID = 0.000000 0.000000 0.000000 0.000000 * * CLPS_ = 0.000000 94.340625 47.170312 47.170312 * * CCUR_ = -269.445000 107.917500 -215.486250 53.958750 * * CKIS_ = 0.000000 -107.917500 -53.958750 -53.958750 * * CROT_ = -75.000000 -75.000000 -112.500000 -37.500000 * * * ******************************************************************************* ******************************************************************************* * * * PHYSICAL CONSTANTS: H_BARC=197.32891000 HBCOE2=137.03602000 * * * * XMASSN=938.90590000 XMASSP=938.27231000 * * * * HBMASS= 20.73000000 HBMRPA= 20.21175000 * * * * ECHAR2= 1.43997841 COULEX= 0.00000000 * * * ******************************************************************************* * * * OSCILLATOR LENGTHS: X= 1.6978805 Y= 1.6978805 Z= 1.6978805 * * * * OSCILLATOR CONSTANTS: X= 0.5889696 Y= 0.5889696 Z= 0.5889696 * * * * OSCILLATOR FREQUENCIES: X=14.3861673 Y=14.3861673 Z=14.3861673 * * * ******************************************************************************* * * * BASIS CUT-OFF CONTROL PARAMETERS: NXMAXX = 10 NYMAXX = 10 NZMAXX = 10 * * * * OPTIMUM NUMBERS OF GAUSS POINTS: NXHERM = 22 NYHERM = 22 NZHERM = 22 * * * * NLIMIT= 286 LDBASE= 286 MCOUNT= 1331 * * * * ENECUT= 800.0000 ELIMIT= 165.4409 * * * ******************************************************************************* ******************************************************************************* * * * SHAPE OF THE OSCILLATOR-BASIS DIAMOND (THREE CONSECUTIVE PROJECTIONS) * * * ******************************************************************************* * * * NX ===>>> 0 1 2 3 4 5 6 7 8 9 10 * * MAX.NY => 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------- * * NY= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NY= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NY= 2 ( 8) | 8 7 6 5 4 3 2 1 0 * * NY= 3 ( 7) | 7 6 5 4 3 2 1 0 * * NY= 4 ( 6) | 6 5 4 3 2 1 0 * * NY= 5 ( 5) | 5 4 3 2 1 0 * * NY= 6 ( 4) | 4 3 2 1 0 * * NY= 7 ( 3) | 3 2 1 0 * * NY= 8 ( 2) | 2 1 0 * * NY= 9 ( 1) | 1 0 * * NY=10 ( 0) | 0 * * * * * * NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 * * MAX.NY => 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------- * * NY= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NY= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NY= 2 ( 8) | 8 7 6 5 4 3 2 1 0 * * NY= 3 ( 7) | 7 6 5 4 3 2 1 0 * * NY= 4 ( 6) | 6 5 4 3 2 1 0 * * NY= 5 ( 5) | 5 4 3 2 1 0 * * NY= 6 ( 4) | 4 3 2 1 0 * * NY= 7 ( 3) | 3 2 1 0 * * NY= 8 ( 2) | 2 1 0 * * NY= 9 ( 1) | 1 0 * * NY=10 ( 0) | 0 * * * * * * NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 * * MAX.NX => 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------- * * NX= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NX= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NX= 2 ( 8) | 8 7 6 5 4 3 2 1 0 * * NX= 3 ( 7) | 7 6 5 4 3 2 1 0 * * NX= 4 ( 6) | 6 5 4 3 2 1 0 * * NX= 5 ( 5) | 5 4 3 2 1 0 * * NX= 6 ( 4) | 4 3 2 1 0 * * NX= 7 ( 3) | 3 2 1 0 * * NX= 8 ( 2) | 2 1 0 * * NX= 9 ( 1) | 1 0 * * NX=10 ( 0) | 0 * * * ******************************************************************************* ******************************************************************************* * * * AUTOMATED HANDLING OF THE COULOMB FILE * * * ******************************************************************************* * * * THE COULOMB FILE READ FROM THE DISC, FILE NAME (FILCOU): * * ca40_iso.cou * * * ******************************************************************************* * * * NUCLIDE: N = 20 Z = 20 NUMBER OF ITERATIONS = 30 CONTINUATION? = 1 * * * ******************************************************************************* * * * ITERATIONS WILL STOP WHEN THE STABILITY BECOMES LOWER THAN EPSITE= 1.0E-09 * * FOR 3 CONSECUTIVE ITERATIONS * * * ******************************************************************************* * * * 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 EXPANDING POTENTIAL 1/R INTO 5 GAUSSIANS * * EXPANSION PARAMETERS CALCULATED BY USING THE GAUSS-LEGENDRE INTEGRATION * * * ******************************************************************************* * * * EXCHANGE COULOMB TERMS SOLVED BY EXPANDING POTENTIAL 1/R INTO 5 GAUSSIANS * * EXPANSION PARAMETERS CALCULATED BY USING THE GAUSS-LEGENDRE INTEGRATION * * * ******************************************************************************* * * * MEAN-FIELD MATRIX ELEMENTS WILL BE READ FROM DISC TO CONTINUE ITERATIONS * * * ******************************************************************************* * * * MAXIMUM NUMBERS OF MULTIPOLE MOMENTS USED FOR THE CONSTRAINTS, NMUCON=2 * * FOR THE COULOMB FIELD, NMUCOU=4 * * FOR THE OUTPUT INFO, NMUPRI=4 * * * ******************************************************************************* * * * PRINTING THE RESULTS FOR THE FOLLOWING ITERATIONS: THE FIRST: NO * * 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: 20 20 * * 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: 20 0 0 * * PROTONS : 20 0 0 * * * ******************************************************************************* ******************************************************************************* * * * CALLING SUBROUTINE "RECORD" TO READ THE REPLAY FILE OBTAINED PREVIOUSLY * * ca40_iso.rec * * * * CALLING SUBROUTINE "RECFIL" TO READ THE FIELDS FILE OBTAINED PREVIOUSLY * * ca40_iso.fil * * * ******************************************************************************* ******************************************************************************* * * * ITERATION RESTARTED FROM THE DISC, THE FOLLOWING FILE(S) READ: * * ca40_iso.rec * * ca40_iso.fil * * * ******************************************************************************* * * * CONVERGENCE REPORT * * * ******************************************************************************* * * * ITER ENERGY STABILITY Q_2 GAMMA SPIN OMEGA ANGLE EPAIR * * * * 44 -342.858837 2.884199 0.003 0.00 0.00 0.000 0.000 0.00 * * 45 -342.858082 -2.267045 0.001 0.00 0.00 0.000 0.000 0.00 * * 46 -342.858925 -1.203311 0.001 0.00 0.00 0.000 0.000 0.00 * * 47 -342.859170 -0.660290 0.000 0.00 0.00 0.000 0.000 0.00 * * 48 -342.859253 -0.368246 0.000 0.00 0.00 0.000 0.000 0.00 * * 49 -342.859285 -0.209253 0.000 0.00 0.00 0.000 0.000 0.00 * * 50 -342.859298 -0.121367 0.000 0.00 0.00 0.000 0.000 0.00 * * 51 -342.859304 -0.071929 0.000 0.00 0.00 0.000 0.000 0.00 * * 52 -342.859306 -0.043573 0.000 0.00 0.00 0.000 0.000 0.00 * * 53 -342.859308 -0.026965 0.000 0.00 0.00 0.000 0.000 0.00 * * 54 -342.859308 -0.017025 0.000 0.00 0.00 0.000 0.000 0.00 * * 55 -342.859309 -0.010946 0.000 0.00 0.00 0.000 0.000 0.00 * * 56 -342.859309 -0.007151 0.000 0.00 0.00 0.000 0.000 0.00 * * 57 -342.859309 -0.004736 0.000 0.00 0.00 0.000 0.000 0.00 * * 58 -342.859309 -0.003172 0.000 0.00 0.00 0.000 0.000 0.00 * * 59 -342.859309 -0.002145 0.000 0.00 0.00 0.000 0.000 0.00 * * 60 -342.859309 -0.001461 0.000 0.00 0.00 0.000 0.000 0.00 * * 61 -342.859309 -0.001002 0.000 0.00 0.00 0.000 0.000 0.00 * * 62 -342.859309 -0.000690 0.000 0.00 0.00 0.000 0.000 0.00 * * 63 -342.859309 -0.000477 0.000 0.00 0.00 0.000 0.000 0.00 * * 64 -342.859309 -0.000331 0.000 0.00 0.00 0.000 0.000 0.00 * * 65 -342.859309 -0.000230 0.000 0.00 0.00 0.000 0.000 0.00 * * 66 -342.859309 -0.000160 0.000 0.00 0.00 0.000 0.000 0.00 * * 67 -342.859309 -0.000112 0.000 0.00 0.00 0.000 0.000 0.00 * * 68 -342.859309 -0.000078 0.000 0.00 0.00 0.000 0.000 0.00 * * 69 -342.859309 -0.000055 0.000 0.00 0.00 0.000 0.000 0.00 * * 70 -342.859309 -0.000038 0.000 0.00 0.00 0.000 0.000 0.00 * * 71 -342.859309 -0.000027 0.000 0.00 0.00 0.000 0.000 0.00 * * 72 -342.859309 -0.000019 0.000 0.00 0.00 0.000 0.000 0.00 * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY | N,nz,lz,OMEz>

* * * * 9) -27.600 0.9993 | 2, 2, 0, 1/2> 100 -0.006 -1.499 -0.019 * * -0.001 -0.300 -0.004 * * 10) -27.600 -0.9993 | 2, 2, 0, 1/2> 100 0.006 1.499 0.019 * * 0.001 0.300 0.004 * * 11) -27.600 -0.9962 | 2, 1, 1, 3/2> 100 0.000 0.000 0.033 * * 0.000 0.000 0.007 * * 12) -27.600 0.9962 | 2, 1, 1, 3/2> 100 0.000 0.000 -0.033 * * 0.000 0.000 -0.007 * * 13) -27.600 -0.9829 | 2, 0, 2, 5/2> 100 0.000 0.000 0.014 * * 0.000 0.000 0.003 * * 14) -27.600 0.9829 | 2, 0, 2, 5/2> 100 0.000 0.000 -0.014 * * 0.000 0.000 -0.003 * * 15) -19.132 1.0000 | 2, 0, 0, 1/2> 100 -0.001 -0.500 -0.001 * * -0.001 -0.500 -0.001 * * 16) -19.132 -1.0000 | 2, 0, 0, 1/2> 100 0.001 0.500 0.001 * * 0.001 0.500 0.001 * * 17) -18.917 1.0000 | 2, 1, 1, 1/2> 100 -0.002 1.000 0.001 * * 0.000 -0.200 0.000 * * 18) -18.917 -1.0000 | 2, 1, 1, 1/2> 100 0.002 -1.000 -0.001 * * 0.000 0.200 0.000 * * 19) -18.917 -0.9544 | 2, 0, 2, 3/2> 100 0.000 0.000 -0.021 * * 0.000 0.000 0.004 * * 20) -18.917 0.9544 | 2, 0, 2, 3/2> 100 0.000 0.000 0.021 * * 0.000 0.000 -0.004 * * 21) -7.386 0.7212 | 3, 2, 1, 1/2> -100 0.137 -1.442 0.345 * * 0.020 -0.206 0.049 * * 22) -7.386 -0.7212 | 3, 2, 1, 1/2> -100 -0.137 1.442 -0.345 * * -0.020 0.206 -0.049 * * 23) -7.386 0.0750 | 3, 2, 1, 3/2> -100 0.000 0.000 1.335 * * 0.000 0.000 0.191 * * 24) -7.386 -0.0750 | 3, 2, 1, 3/2> -100 0.000 0.000 -1.335 * * 0.000 0.000 -0.191 * * 25) -7.386 -0.0549 | 3, 1, 2, 5/2> -100 0.000 0.000 0.177 * * 0.000 0.000 0.025 * * 26) -7.386 0.0549 | 3, 1, 2, 5/2> -100 0.000 0.000 -0.177 * * 0.000 0.000 -0.025 * * 27) -7.386 0.8055 | 3, 0, 3, 7/2> -100 0.000 0.000 -0.948 * * 0.000 0.000 -0.135 * * 28) -7.386 -0.8055 | 3, 0, 3, 7/2> -100 0.000 0.000 0.948 * * 0.000 0.000 0.135 * * 29) -1.387 0.9517 | 3, 3, 0, 1/2> -100 -0.058 -0.952 -0.151 * * -0.019 -0.317 -0.050 * * 30) -1.387 -0.9517 | 3, 3, 0, 1/2> -100 0.058 0.952 0.151 * * 0.019 0.317 0.050 * * 31) -1.387 -0.8875 | 3, 0, 1, 3/2> -100 0.000 0.000 0.238 * * 0.000 0.000 0.079 * * 32) -1.387 0.8875 | 3, 0, 1, 3/2> -100 0.000 0.000 -0.238 * * 0.000 0.000 -0.079 * * 33) 0.781 -0.7017 | 3, 0, 1, 1/2> -100 0.241 -0.351 -0.262 * * -0.080 0.117 0.087 * * 34) 0.781 0.7017 | 3, 0, 1, 1/2> -100 -0.241 0.351 0.262 * * 0.080 -0.117 -0.087 * * 35) 2.197 -0.9242 | 3, 2, 1, 1/2> -100 -0.096 -1.386 0.188 * * 0.014 0.198 -0.027 * * 36) 2.197 0.9242 | 3, 2, 1, 1/2> -100 0.096 1.386 -0.188 * * -0.014 -0.198 0.027 * * 37) 2.197 0.6222 | 3, 1, 2, 3/2> -100 0.000 0.000 0.964 * * 0.000 0.000 -0.138 * * 38) 2.197 -0.6222 | 3, 1, 2, 3/2> -100 0.000 0.000 -0.964 * * 0.000 0.000 0.138 * * 39) 2.197 -0.1320 | 3, 0, 3, 5/2> -100 0.000 0.000 -0.525 * * 0.000 0.000 0.075 * * 40) 2.197 0.1320 | 3, 0, 3, 5/2> -100 0.000 0.000 0.525 * * 0.000 0.000 -0.075 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] NEUTRONS * * * * 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 = 6.0E-08 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 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS NEUTRONS * * * * R_RMS = 3.3240 X_RMS = 1.9191 Y_RMS = 1.9191 Z_RMS = 1.9191 * * * * R_GEO = 4.2913 X_GEO = 4.2913 Y_GEO = 4.2913 Z_GEO = 4.2913 * * * ******************************************************************************* ******************************************************************************* * * * 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 = 2.1E-08 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.0604 FM^(-3) BOHR RADIUS = 4.2913 FM * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS * * * ******************************************************************************* * * * NO) ENERGY | N,nz,lz,OMEz>

* * * * 9) -20.467 0.9980 | 2, 2, 0, 1/2> 100 -0.041 -1.497 -0.028 * * -0.008 -0.299 -0.006 * * 10) -20.467 -0.9980 | 2, 2, 0, 1/2> 100 0.041 1.497 0.028 * * 0.008 0.299 0.006 * * 11) -20.467 -0.9988 | 2, 1, 1, 3/2> 100 0.000 0.000 0.052 * * 0.000 0.000 0.010 * * 12) -20.467 0.9988 | 2, 1, 1, 3/2> 100 0.000 0.000 -0.052 * * 0.000 0.000 -0.010 * * 13) -20.467 -0.9988 | 2, 0, 2, 5/2> 100 0.000 0.000 0.041 * * 0.000 0.000 0.008 * * 14) -20.467 0.9988 | 2, 0, 2, 5/2> 100 0.000 0.000 -0.041 * * 0.000 0.000 -0.008 * * 15) -12.166 0.9851 | 2, 0, 0, 1/2> 100 0.085 -0.493 -0.012 * * 0.085 -0.493 -0.012 * * 16) -12.166 -0.9851 | 2, 0, 0, 1/2> 100 -0.085 0.493 0.012 * * -0.085 0.493 0.012 * * 17) -12.018 1.0000 | 2, 1, 1, 1/2> 100 0.000 1.000 0.002 * * 0.000 -0.200 0.000 * * 18) -12.018 -1.0000 | 2, 1, 1, 1/2> 100 0.000 -1.000 -0.002 * * 0.000 0.200 0.000 * * 19) -12.018 0.9848 | 2, 0, 2, 3/2> 100 0.000 0.000 0.015 * * 0.000 0.000 -0.003 * * 20) -12.018 -0.9848 | 2, 0, 2, 3/2> 100 0.000 0.000 -0.015 * * 0.000 0.000 0.003 * * 21) -0.637 0.9434 | 3, 2, 1, 1/2> -100 -0.247 -1.887 0.154 * * -0.035 -0.270 0.022 * * 22) -0.637 -0.9434 | 3, 2, 1, 1/2> -100 0.247 1.887 -0.154 * * 0.035 0.270 -0.022 * * 23) -0.637 -0.6713 | 3, 2, 1, 3/2> -100 0.000 0.000 -1.102 * * 0.000 0.000 -0.157 * * 24) -0.637 0.6713 | 3, 2, 1, 3/2> -100 0.000 0.000 1.102 * * 0.000 0.000 0.157 * * 25) -0.637 0.6364 | 3, 1, 2, 5/2> -100 0.000 0.000 0.949 * * 0.000 0.000 0.136 * * 26) -0.637 -0.6364 | 3, 1, 2, 5/2> -100 0.000 0.000 -0.949 * * 0.000 0.000 -0.136 * * 27) -0.637 0.6119 | 3, 0, 3, 7/2> -100 0.000 0.000 0.994 * * 0.000 0.000 0.142 * * 28) -0.637 -0.6119 | 3, 0, 3, 7/2> -100 0.000 0.000 -0.994 * * 0.000 0.000 -0.142 * * 29) 4.893 0.9785 | 3, 3, 0, 1/2> -100 -0.125 -0.978 0.082 * * -0.042 -0.326 0.027 * * 30) 4.893 -0.9785 | 3, 3, 0, 1/2> -100 0.125 0.978 -0.082 * * 0.042 0.326 -0.027 * * 31) 4.893 0.7788 | 3, 0, 1, 3/2> -100 0.000 0.000 -0.442 * * 0.000 0.000 -0.147 * * 32) 4.893 -0.7788 | 3, 0, 1, 3/2> -100 0.000 0.000 0.442 * * 0.000 0.000 0.147 * * 33) 6.836 -0.9913 | 3, 0, 1, 1/2> -100 -0.062 -0.496 0.023 * * 0.021 0.165 -0.008 * * 34) 6.836 0.9913 | 3, 0, 1, 1/2> -100 0.062 0.496 -0.023 * * -0.021 -0.165 0.008 * * 35) 8.419 0.9178 | 3, 2, 1, 1/2> -100 -0.516 1.377 -0.099 * * 0.074 -0.197 0.014 * * 36) 8.419 -0.9178 | 3, 2, 1, 1/2> -100 0.516 -1.377 0.099 * * -0.074 0.197 -0.014 * * 37) 8.419 -0.7175 | 3, 1, 2, 3/2> -100 0.000 0.000 0.890 * * 0.000 0.000 -0.127 * * 38) 8.419 0.7175 | 3, 1, 2, 3/2> -100 0.000 0.000 -0.890 * * 0.000 0.000 0.127 * * 39) 8.419 0.6363 | 3, 0, 3, 5/2> -100 0.000 0.000 0.272 * * 0.000 0.000 -0.039 * * 40) 8.419 -0.6363 | 3, 0, 3, 5/2> -100 0.000 0.000 -0.272 * * 0.000 0.000 0.039 * * * ******************************************************************************* ******************************************************************************* * * * 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 = 6.1E-08 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 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS PROTONS * * * * R_RMS = 3.3693 X_RMS = 1.9452 Y_RMS = 1.9452 Z_RMS = 1.9452 * * * * R_GEO = 4.3497 X_GEO = 4.3497 Y_GEO = 4.3497 Z_GEO = 4.3497 * * * ******************************************************************************* ******************************************************************************* * * * 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 = 2.1E-08 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.0580 FM^(-3) BOHR RADIUS = 4.3497 FM * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 4.245751 2.123358 4.245751 0.000964 * * DRHOD = 0.556756 0.278443 0.556756 0.000131 * * DLPR_ = -1.547172 -0.774112 -1.547172 -0.001051 * * DTAU_ = 3.342615 1.671862 3.342615 0.001110 * * DSCU_ = 0.000247 0.000124 0.000247 0.000001 * * DDIV_ = 0.013928 0.006951 0.013928 -0.000025 * * * * 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_ = -2424.717561 437.343382 -1987.473410 0.099232 * * ERHOD = 0.000000 0.000000 0.000000 0.000000 * * ELPR_ = 271.186351 -73.030177 198.205771 -0.049596 * * ETAU_ = 900.650954 -180.423214 720.287615 -0.059876 * * ESCU_ = 0.000000 0.013373 0.013345 0.000028 * * EDIV_ = -1.044596 -0.521347 -1.566894 0.000951 * * ============ ============ ============ ============ * * SUM EVEN: -1253.924852 183.382017 -1070.533574 -0.009261 * * * * 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 = 20.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 6.0E-08 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) [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 = 2.1E-08 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.0604 FM^(-3) BOHR RADIUS = 4.2913 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 = 6.1E-08 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) [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 = 2.1E-08 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.0580 FM^(-3) BOHR RADIUS = 4.3497 FM * * * ******************************************************************************* ******************************************************************************* * * * EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL * * * * ALPHA =259.60687 BETA = 0.00000 GAMMA = 0.00000 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 3.3467 X_RMS = 1.9322 Y_RMS = 1.9322 Z_RMS = 1.9322 * * * * R_GEO = 4.3206 X_GEO = 4.3206 Y_GEO = 4.3206 Z_GEO = 4.3206 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 40.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 1.2E-07 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) TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = 2.1E-08 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.1184 FM^(-3) BOHR RADIUS = 4.3206 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 = 40.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 1.2E-07 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) [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 = 2.1E-08 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.1184 FM^(-3) BOHR RADIUS = 4.3206 FM * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] [INTRINSIC FRAME] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 40.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 1.2E-07 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 = 2.1E-08 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.1184 FM^(-3) BOHR RADIUS = 4.3206 FM * * * ******************************************************************************* ******************************************************************************* * * * LINEAR MOMENTA IN FERMI^(-1) TOTAL * * * * P_X = 0.0000 P_Y = 0.0000 P_Z = 0.0000 P_R = 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * ANGULAR MOMENTA IN THE THREE CARTESIAN DIRECTIONS * * * ******************************************************************************* * * * SPIN IN X-DIRECTION SPIN IN Y-DIRECTION SPIN IN Z-DIRECTION * * ------------------- ------------------- ------------------- * * INTRI ORBIT TOTAL INTRI ORBIT TOTAL INTRI ORBIT TOTAL * * * * NEUTRONS 0.000 0.000 0.000 0.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 * * ======================= * * 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 * * * ******************************************************************************* * * * 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)= 331.798302 (PRO)= 323.126285 (TOT)= 654.924587 * * SUM EPS: (NEU)= -743.128900 (PRO)= -597.514331 (TOT)= -1340.643231 * * PAIRING: (NEU)= 0.000000 (PRO)= 0.000000 (TOT)= 0.000000 * * P-REARR: (NEU)= 0.000000 (PRO)= 0.000000 (TOT)= 0.000000 * * * * PAIRGAP: (NEU)= 0.000000 (PRO)= 0.000000 THIS SPACE IS EMPTY * * E-FERMI: (NEU)= -13.151954 (PRO)= -6.327659 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 80.716917 (EXC)= -7.957977 (TOT)= 72.758940 * * COULEXC: (SCA)= -7.955261 (VEC)= -0.002716 (TOT)= -7.957977 * * * * 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= 0.000000 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -342.859309 * * * * SPIN-ORB (EVE)= -1.565943 (ODD)= 0.000000 (TOT)= -1.565943 * * SKYRME: (EVE)= -1070.542835 (ODD)= 0.000000 (TOT)= -1070.542835 * * * * TOTAL: (STAB)= -0.000013 (SP)= -342.859322 (FUN)= -342.859309 * * * ******************************************************************************* ******************************************************************************* * * * CALLING SUBROUTINE "RECORD" TO SAVE THE RECORD FILE CONTAINING RESULTS * * ca40_iso.rec * * THE MEAN FIELDS OBTAINED IN THE LAST-BUT-ONE ITERATION HAVE BEEN STORED * * * * CALLING SUBROUTINE "RECFIL" TO SAVE THE FIELDS FILE CONTAINING RESULTS * * ca40_iso.fil * * * ******************************************************************************* ******************************************************************************* * * * WARNING! THE STABILITY REQUIRED WAS NOT ACHIEVED. NO-OF-ITERATIONS EXCEEDED * * WARNING! THE STABILITY REQUIRED WAS NOT ACHIEVED. NO-OF-ITERATIONS EXCEEDED * * * ******************************************************************************* * * ******************************************************************************* * * * CALCULATION OF THE DATA SET NO. 2 HAS BEEN SUCCESSFULLY COMPLETED * * * ******************************************************************************* ******************************************************************************* * * * CLASSICAL NUCLEAR SURFACE DEFINED FOR: N = 20 Z = 20 * * * ******************************************************************************* * * * SPHERICAL SHAPE REQUESTED FOR THE SURFACE * * * ******************************************************************************* * * * HOMEGA= 14.3862 FCHOM0= 1.2000 * * * * OSCILLATOR FREQUENCIES: HBAROX= 14.3862 HBAROY= 14.3862 HBAROZ= 14.3862 * * * * MOMENTS OF INERTIA: XMOMFC= 6.8244 YMOMFC= 6.8244 ZMOMFC= 6.8244 * * * * CENTERS OF MASS: CMSXFC= 0.0000 CMSYFC= 0.0000 CMSZFC= 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * SKYRME FORCE DEFINITION * * * ******************************************************************************* * * * PARAMETER SET SV : T0= -1248.29 T1= 970.56 T2= 107.22 T3= 0.00 * * * * POWER=1/1 W=150.000 X0= -0.17000 X1= 0.00000 X2= 0.00000 X3= 1.00000 * * * * ETA_J=0 ETA_W=0 ETACM=0 ETA_M=1 HBM= 20.73000 (SKYRME-FORCE SPIN-ORBIT) * * * ******************************************************************************* ******************************************************************************* * * * NON-STANDARD TREATMENT OF THE SKYRME FORCE: ETA_J=1 ETA_W=0 ETACM=0 ETA_M=1 * * * * THIS INCLUDES: * * * * NON-STANDARD TREATMENT OF THE TENSOR TERMS: THESE TERMS HAVE BEEN INCLUDED * * * ******************************************************************************* ******************************************************************************* * * * COUPLING CONSTANTS DEFINING THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * CRHO_ = -571.092675 205.967850 -468.108750 102.983925 * * CRHOD = 0.000000 0.000000 0.000000 0.000000 * * CLPR_ = -175.278750 94.340625 -128.108438 47.170312 * * CTAU_ = 269.445000 -107.917500 215.486250 -53.958750 * * CSCU_ = 0.000000 107.917500 53.958750 53.958750 * * CDIV_ = -75.000000 -75.000000 -112.500000 -37.500000 * * * * CSPI_ = 53.052325 312.072500 209.088575 156.036250 * * CSPID = 0.000000 0.000000 0.000000 0.000000 * * CLPS_ = 0.000000 94.340625 47.170312 47.170312 * * CCUR_ = -269.445000 107.917500 -215.486250 53.958750 * * CKIS_ = 0.000000 -107.917500 -53.958750 -53.958750 * * CROT_ = -75.000000 -75.000000 -112.500000 -37.500000 * * * ******************************************************************************* ******************************************************************************* * * * PHYSICAL CONSTANTS: H_BARC=197.32891000 HBCOE2=137.03602000 * * * * XMASSN=938.90590000 XMASSP=938.27231000 * * * * HBMASS= 20.73000000 HBMRPA= 20.21175000 * * * * ECHAR2= 1.43997841 COULEX= 0.00000000 * * * ******************************************************************************* * * * OSCILLATOR LENGTHS: X= 1.6978805 Y= 1.6978805 Z= 1.6978805 * * * * OSCILLATOR CONSTANTS: X= 0.5889696 Y= 0.5889696 Z= 0.5889696 * * * * OSCILLATOR FREQUENCIES: X=14.3861673 Y=14.3861673 Z=14.3861673 * * * ******************************************************************************* * * * BASIS CUT-OFF CONTROL PARAMETERS: NXMAXX = 10 NYMAXX = 10 NZMAXX = 10 * * * * OPTIMUM NUMBERS OF GAUSS POINTS: NXHERM = 22 NYHERM = 22 NZHERM = 22 * * * * NLIMIT= 286 LDBASE= 286 MCOUNT= 1331 * * * * ENECUT= 800.0000 ELIMIT= 165.4409 * * * ******************************************************************************* ******************************************************************************* * * * SHAPE OF THE OSCILLATOR-BASIS DIAMOND (THREE CONSECUTIVE PROJECTIONS) * * * ******************************************************************************* * * * NX ===>>> 0 1 2 3 4 5 6 7 8 9 10 * * MAX.NY => 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------- * * NY= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NY= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NY= 2 ( 8) | 8 7 6 5 4 3 2 1 0 * * NY= 3 ( 7) | 7 6 5 4 3 2 1 0 * * NY= 4 ( 6) | 6 5 4 3 2 1 0 * * NY= 5 ( 5) | 5 4 3 2 1 0 * * NY= 6 ( 4) | 4 3 2 1 0 * * NY= 7 ( 3) | 3 2 1 0 * * NY= 8 ( 2) | 2 1 0 * * NY= 9 ( 1) | 1 0 * * NY=10 ( 0) | 0 * * * * * * NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 * * MAX.NY => 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------- * * NY= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NY= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NY= 2 ( 8) | 8 7 6 5 4 3 2 1 0 * * NY= 3 ( 7) | 7 6 5 4 3 2 1 0 * * NY= 4 ( 6) | 6 5 4 3 2 1 0 * * NY= 5 ( 5) | 5 4 3 2 1 0 * * NY= 6 ( 4) | 4 3 2 1 0 * * NY= 7 ( 3) | 3 2 1 0 * * NY= 8 ( 2) | 2 1 0 * * NY= 9 ( 1) | 1 0 * * NY=10 ( 0) | 0 * * * * * * NZ ===>>> 0 1 2 3 4 5 6 7 8 9 10 * * MAX.NX => 10 9 8 7 6 5 4 3 2 1 0 * * --------------------------------- * * NX= 0 (10) | 10 9 8 7 6 5 4 3 2 1 0 * * NX= 1 ( 9) | 9 8 7 6 5 4 3 2 1 0 * * NX= 2 ( 8) | 8 7 6 5 4 3 2 1 0 * * NX= 3 ( 7) | 7 6 5 4 3 2 1 0 * * NX= 4 ( 6) | 6 5 4 3 2 1 0 * * NX= 5 ( 5) | 5 4 3 2 1 0 * * NX= 6 ( 4) | 4 3 2 1 0 * * NX= 7 ( 3) | 3 2 1 0 * * NX= 8 ( 2) | 2 1 0 * * NX= 9 ( 1) | 1 0 * * NX=10 ( 0) | 0 * * * ******************************************************************************* ******************************************************************************* * * * AUTOMATED HANDLING OF THE COULOMB FILE * * * ******************************************************************************* * * * THE COULOMB FILE READ FROM THE DISC, FILE NAME (FILCOU): * * ca40_iso.cou * * * ******************************************************************************* * * * NUCLIDE: N = 20 Z = 20 NUMBER OF ITERATIONS = 1 CONTINUATION? = 1 * * * ******************************************************************************* * * * ITERATIONS WILL STOP WHEN THE STABILITY BECOMES LOWER THAN EPSITE= 1.0E-09 * * FOR 3 CONSECUTIVE ITERATIONS * * * ******************************************************************************* * * * SLOW-DOWN PARAMETERS: * * TIME-EVEN = 1.00 TIME-ODD = 1.00 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 EXPANDING POTENTIAL 1/R INTO 5 GAUSSIANS * * EXPANSION PARAMETERS CALCULATED BY USING THE GAUSS-LEGENDRE INTEGRATION * * * ******************************************************************************* * * * EXCHANGE COULOMB TERMS SOLVED BY EXPANDING POTENTIAL 1/R INTO 5 GAUSSIANS * * EXPANSION PARAMETERS CALCULATED BY USING THE GAUSS-LEGENDRE INTEGRATION * * * ******************************************************************************* * * * MEAN-FIELD MATRIX ELEMENTS WILL BE READ FROM DISC TO CONTINUE ITERATIONS * * * ******************************************************************************* * * * MAXIMUM NUMBERS OF MULTIPOLE MOMENTS USED FOR THE CONSTRAINTS, NMUCON=2 * * FOR THE COULOMB FIELD, NMUCOU=4 * * FOR THE OUTPUT INFO, NMUPRI=4 * * * ******************************************************************************* * * * PRINTING THE RESULTS FOR THE FOLLOWING ITERATIONS: THE FIRST: NO * * 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: 20 20 * * PROTONS: 20 20 * * * ******************************************************************************* * * * REDUCED KERNELS AND MATRIX ELEMENTS: * * FOR ELECTRIC MULTIPOLES CALCULATED UP TO = 2, PRINTED IF LARGER THAN 0.000 * * FOR MAGNETIC MULTIPOLES CALCULATED UP TO = 1, PRINTED IF LARGER THAN 0.000 * * FOR SURFACE MULTIPOLES CALCULATED UP TO = 0, PRINTED IF LARGER THAN 0.000 * * * ******************************************************************************* * * * PRINTING THE RESULTS FOR ANGULAR-MOMENTUM-PROJECTED STATES: 0 < 2*I < 999 * * AVERAGE ENERGIES OF PROJECTED STATES: YES * * KERNELS BETWEEN THE PROJECTED STATES: NO * * SUM RULES COMPARED TO THE H-F VALUES: YES * * K-MIXED ENERGIES OF PROJECTED STATES: YES * * REDUCED K E R N E L S BETWEEN THE PROJECTED STATES: NO * * REDUCED MATRIX ELEMENTS BETWEEN THE K-MIXED STATES: NO * * * ******************************************************************************* * * * THE GCM KERNELS IN FUNCTION OF THE EULER ANGLES WILL BE READ AND/OR STORED * * * ******************************************************************************* * * * CALCULATIONS WITHOUT PAIRING * * * ******************************************************************************* * * * CALCULATIONS WITHOUT SIMPLEX, SIGNATURE, OR PARITY SYMMETRIES * * * ******************************************************************************* * * * CALCULATIONS WITH TIME-ODD SYMMETRIES: T*SIMPLEX_X => NON CONSERVED * * T*SIMPLEX_Y => NON CONSERVED * * T*SIMPLEX_Z => NON CONSERVED * * * ******************************************************************************* * * * CALCULATIONS WITH BROKEN TIME-REVERSAL SYMMETRY * * * ******************************************************************************* * * * CALCULATIONS WITH ONE-DIMENSIONAL CRANKING FOR THE ROTATIONAL FREQUENCIES: * * * * OMEGAY = 0.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: 20 0 0 * * PROTONS : 20 0 0 * * * ******************************************************************************* ******************************************************************************* * * * CALLING SUBROUTINE "RECORD" TO READ THE REPLAY FILE OBTAINED PREVIOUSLY * * ca40_iso.rec * * * * CALLING SUBROUTINE "RECFIL" TO READ THE FIELDS FILE OBTAINED PREVIOUSLY * * ca40_iso.fil * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY | N,nz,lz,OMEz>

* * * * 9) -27.600 0.9993 | 2, 2, 0, 1/2> 100 -0.006 -1.499 -0.019 * * -0.001 -0.300 -0.004 * * 10) -27.600 -0.9993 | 2, 2, 0, 1/2> 100 0.006 1.499 0.019 * * 0.001 0.300 0.004 * * 11) -27.600 -0.9962 | 2, 1, 1, 3/2> 100 0.000 0.000 0.033 * * 0.000 0.000 0.007 * * 12) -27.600 0.9962 | 2, 1, 1, 3/2> 100 0.000 0.000 -0.033 * * 0.000 0.000 -0.007 * * 13) -27.600 -0.9829 | 2, 0, 2, 5/2> 100 0.000 0.000 0.014 * * 0.000 0.000 0.003 * * 14) -27.600 0.9829 | 2, 0, 2, 5/2> 100 0.000 0.000 -0.014 * * 0.000 0.000 -0.003 * * 15) -19.132 1.0000 | 2, 0, 0, 1/2> 100 -0.001 -0.500 -0.001 * * -0.001 -0.500 -0.001 * * 16) -19.132 -1.0000 | 2, 0, 0, 1/2> 100 0.001 0.500 0.001 * * 0.001 0.500 0.001 * * 17) -18.917 1.0000 | 2, 1, 1, 1/2> 100 -0.002 1.000 0.001 * * 0.000 -0.200 0.000 * * 18) -18.917 -1.0000 | 2, 1, 1, 1/2> 100 0.002 -1.000 -0.001 * * 0.000 0.200 0.000 * * 19) -18.917 -0.9544 | 2, 0, 2, 3/2> 100 0.000 0.000 -0.021 * * 0.000 0.000 0.004 * * 20) -18.917 0.9544 | 2, 0, 2, 3/2> 100 0.000 0.000 0.021 * * 0.000 0.000 -0.004 * * 21) -7.386 0.7212 | 3, 2, 1, 1/2> -100 0.137 -1.442 0.345 * * 0.020 -0.206 0.049 * * 22) -7.386 -0.7212 | 3, 2, 1, 1/2> -100 -0.137 1.442 -0.345 * * -0.020 0.206 -0.049 * * 23) -7.386 0.0750 | 3, 2, 1, 3/2> -100 0.000 0.000 1.335 * * 0.000 0.000 0.191 * * 24) -7.386 -0.0750 | 3, 2, 1, 3/2> -100 0.000 0.000 -1.335 * * 0.000 0.000 -0.191 * * 25) -7.386 -0.0549 | 3, 1, 2, 5/2> -100 0.000 0.000 0.177 * * 0.000 0.000 0.025 * * 26) -7.386 0.0549 | 3, 1, 2, 5/2> -100 0.000 0.000 -0.177 * * 0.000 0.000 -0.025 * * 27) -7.386 0.8055 | 3, 0, 3, 7/2> -100 0.000 0.000 -0.948 * * 0.000 0.000 -0.135 * * 28) -7.386 -0.8055 | 3, 0, 3, 7/2> -100 0.000 0.000 0.948 * * 0.000 0.000 0.135 * * 29) -1.387 0.9517 | 3, 3, 0, 1/2> -100 -0.058 -0.952 -0.151 * * -0.019 -0.317 -0.050 * * 30) -1.387 -0.9517 | 3, 3, 0, 1/2> -100 0.058 0.952 0.151 * * 0.019 0.317 0.050 * * 31) -1.387 -0.8875 | 3, 0, 1, 3/2> -100 0.000 0.000 0.238 * * 0.000 0.000 0.079 * * 32) -1.387 0.8875 | 3, 0, 1, 3/2> -100 0.000 0.000 -0.238 * * 0.000 0.000 -0.079 * * 33) 0.781 -0.7017 | 3, 0, 1, 1/2> -100 0.241 -0.351 -0.262 * * -0.080 0.117 0.087 * * 34) 0.781 0.7017 | 3, 0, 1, 1/2> -100 -0.241 0.351 0.262 * * 0.080 -0.117 -0.087 * * 35) 2.197 -0.9242 | 3, 2, 1, 1/2> -100 -0.096 -1.386 0.188 * * 0.014 0.198 -0.027 * * 36) 2.197 0.9242 | 3, 2, 1, 1/2> -100 0.096 1.386 -0.188 * * -0.014 -0.198 0.027 * * 37) 2.197 0.6222 | 3, 1, 2, 3/2> -100 0.000 0.000 0.964 * * 0.000 0.000 -0.138 * * 38) 2.197 -0.6222 | 3, 1, 2, 3/2> -100 0.000 0.000 -0.964 * * 0.000 0.000 0.138 * * 39) 2.197 -0.1320 | 3, 0, 3, 5/2> -100 0.000 0.000 -0.525 * * 0.000 0.000 0.075 * * 40) 2.197 0.1320 | 3, 0, 3, 5/2> -100 0.000 0.000 0.525 * * 0.000 0.000 -0.075 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] NEUTRONS * * * * 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 = 6.0E-08 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 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS NEUTRONS * * * * R_RMS = 3.3240 X_RMS = 1.9191 Y_RMS = 1.9191 Z_RMS = 1.9191 * * * * R_GEO = 4.2913 X_GEO = 4.2913 Y_GEO = 4.2913 Z_GEO = 4.2913 * * * ******************************************************************************* ******************************************************************************* * * * 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 = 2.1E-08 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.0604 FM^(-3) BOHR RADIUS = 4.2913 FM * * * ******************************************************************************* ******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK PROTONS * * * ******************************************************************************* * * * NO) ENERGY | N,nz,lz,OMEz>

* * * * 9) -20.467 0.9980 | 2, 2, 0, 1/2> 100 -0.041 -1.497 -0.028 * * -0.008 -0.299 -0.006 * * 10) -20.467 -0.9980 | 2, 2, 0, 1/2> 100 0.041 1.497 0.028 * * 0.008 0.299 0.006 * * 11) -20.467 -0.9988 | 2, 1, 1, 3/2> 100 0.000 0.000 0.052 * * 0.000 0.000 0.010 * * 12) -20.467 0.9988 | 2, 1, 1, 3/2> 100 0.000 0.000 -0.052 * * 0.000 0.000 -0.010 * * 13) -20.467 -0.9988 | 2, 0, 2, 5/2> 100 0.000 0.000 0.041 * * 0.000 0.000 0.008 * * 14) -20.467 0.9988 | 2, 0, 2, 5/2> 100 0.000 0.000 -0.041 * * 0.000 0.000 -0.008 * * 15) -12.166 0.9851 | 2, 0, 0, 1/2> 100 0.085 -0.493 -0.012 * * 0.085 -0.493 -0.012 * * 16) -12.166 -0.9851 | 2, 0, 0, 1/2> 100 -0.085 0.493 0.012 * * -0.085 0.493 0.012 * * 17) -12.018 1.0000 | 2, 1, 1, 1/2> 100 0.000 1.000 0.002 * * 0.000 -0.200 0.000 * * 18) -12.018 -1.0000 | 2, 1, 1, 1/2> 100 0.000 -1.000 -0.002 * * 0.000 0.200 0.000 * * 19) -12.018 0.9848 | 2, 0, 2, 3/2> 100 0.000 0.000 0.015 * * 0.000 0.000 -0.003 * * 20) -12.018 -0.9848 | 2, 0, 2, 3/2> 100 0.000 0.000 -0.015 * * 0.000 0.000 0.003 * * 21) -0.637 0.9434 | 3, 2, 1, 1/2> -100 -0.247 -1.887 0.154 * * -0.035 -0.270 0.022 * * 22) -0.637 -0.9434 | 3, 2, 1, 1/2> -100 0.247 1.887 -0.154 * * 0.035 0.270 -0.022 * * 23) -0.637 -0.6713 | 3, 2, 1, 3/2> -100 0.000 0.000 -1.102 * * 0.000 0.000 -0.157 * * 24) -0.637 0.6713 | 3, 2, 1, 3/2> -100 0.000 0.000 1.102 * * 0.000 0.000 0.157 * * 25) -0.637 0.6364 | 3, 1, 2, 5/2> -100 0.000 0.000 0.949 * * 0.000 0.000 0.136 * * 26) -0.637 -0.6364 | 3, 1, 2, 5/2> -100 0.000 0.000 -0.949 * * 0.000 0.000 -0.136 * * 27) -0.637 0.6119 | 3, 0, 3, 7/2> -100 0.000 0.000 0.994 * * 0.000 0.000 0.142 * * 28) -0.637 -0.6119 | 3, 0, 3, 7/2> -100 0.000 0.000 -0.994 * * 0.000 0.000 -0.142 * * 29) 4.893 0.9785 | 3, 3, 0, 1/2> -100 -0.125 -0.978 0.082 * * -0.042 -0.326 0.027 * * 30) 4.893 -0.9785 | 3, 3, 0, 1/2> -100 0.125 0.978 -0.082 * * 0.042 0.326 -0.027 * * 31) 4.893 0.7788 | 3, 0, 1, 3/2> -100 0.000 0.000 -0.442 * * 0.000 0.000 -0.147 * * 32) 4.893 -0.7788 | 3, 0, 1, 3/2> -100 0.000 0.000 0.442 * * 0.000 0.000 0.147 * * 33) 6.836 -0.9913 | 3, 0, 1, 1/2> -100 -0.062 -0.496 0.023 * * 0.021 0.165 -0.008 * * 34) 6.836 0.9913 | 3, 0, 1, 1/2> -100 0.062 0.496 -0.023 * * -0.021 -0.165 0.008 * * 35) 8.419 0.9178 | 3, 2, 1, 1/2> -100 -0.516 1.377 -0.099 * * 0.074 -0.197 0.014 * * 36) 8.419 -0.9178 | 3, 2, 1, 1/2> -100 0.516 -1.377 0.099 * * -0.074 0.197 -0.014 * * 37) 8.419 -0.7175 | 3, 1, 2, 3/2> -100 0.000 0.000 0.890 * * 0.000 0.000 -0.127 * * 38) 8.419 0.7175 | 3, 1, 2, 3/2> -100 0.000 0.000 -0.890 * * 0.000 0.000 0.127 * * 39) 8.419 0.6363 | 3, 0, 3, 5/2> -100 0.000 0.000 0.272 * * 0.000 0.000 -0.039 * * 40) 8.419 -0.6363 | 3, 0, 3, 5/2> -100 0.000 0.000 -0.272 * * 0.000 0.000 0.039 * * * ******************************************************************************* ******************************************************************************* * * * 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 = 6.1E-08 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 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS PROTONS * * * * R_RMS = 3.3693 X_RMS = 1.9452 Y_RMS = 1.9452 Z_RMS = 1.9452 * * * * R_GEO = 4.3497 X_GEO = 4.3497 Y_GEO = 4.3497 Z_GEO = 4.3497 * * * ******************************************************************************* ******************************************************************************* * * * 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 = 2.1E-08 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.0580 FM^(-3) BOHR RADIUS = 4.3497 FM * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 4.245751 2.123358 4.245751 0.000964 * * DRHOD = 0.556756 0.278443 0.556756 0.000131 * * DLPR_ = -1.547172 -0.774112 -1.547172 -0.001051 * * DTAU_ = 3.342615 1.671862 3.342615 0.001110 * * DSCU_ = 0.000247 0.000124 0.000247 0.000001 * * DDIV_ = 0.013928 0.006951 0.013928 -0.000025 * * * * 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_ = -2424.717561 437.343382 -1987.473410 0.099232 * * ERHOD = 0.000000 0.000000 0.000000 0.000000 * * ELPR_ = 271.186351 -73.030177 198.205771 -0.049596 * * ETAU_ = 900.650954 -180.423214 720.287615 -0.059876 * * ESCU_ = 0.000000 0.013373 0.013345 0.000028 * * EDIV_ = -1.044596 -0.521347 -1.566894 0.000951 * * ============ ============ ============ ============ * * SUM EVEN: -1253.924852 183.382017 -1070.533574 -0.009261 * * * * 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 = 20.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 6.0E-08 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) [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 = 2.1E-08 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.0604 FM^(-3) BOHR RADIUS = 4.2913 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 = 6.1E-08 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) [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 = 2.1E-08 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.0580 FM^(-3) BOHR RADIUS = 4.3497 FM * * * ******************************************************************************* ******************************************************************************* * * * EULER ANGLES OF THE PRINCIPAL-AXES FRAME IN DEGREES TOTAL * * * * ALPHA =259.60687 BETA = 0.00000 GAMMA = 0.00000 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 3.3467 X_RMS = 1.9322 Y_RMS = 1.9322 Z_RMS = 1.9322 * * * * R_GEO = 4.3206 X_GEO = 4.3206 Y_GEO = 4.3206 Z_GEO = 4.3206 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 40.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 1.2E-07 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) TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * B10 = ZERO B1+1= ZERO ............. ............. ............. * * B1-1= ZERO ............. ............. ............. * * * * B20 = 2.1E-08 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.1184 FM^(-3) BOHR RADIUS = 4.3206 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 = 40.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 1.2E-07 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) [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 = 2.1E-08 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.1184 FM^(-3) BOHR RADIUS = 4.3206 FM * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS [UNITS: (10 FERMI)^LAMBDA] [INTRINSIC FRAME] TOTAL * * * * REAL PART FOR A NON-NEGATIVE PROJECTION * * IMAGINARY PART FOR A NEGATIVE PROJECTION * * * ******************************************************************************* * * * Q00 = 40.0000 ............. ............. ............. ............. * * ............. ............. ............. ............. * * * * Q10 = ZERO Q1+1= ZERO ............. ............. ............. * * Q1-1= ZERO ............. ............. ............. * * * * Q20 = 1.2E-07 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 = 2.1E-08 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.1184 FM^(-3) BOHR RADIUS = 4.3206 FM * * * ******************************************************************************* ******************************************************************************* * * * LINEAR MOMENTA IN FERMI^(-1) TOTAL * * * * P_X = 0.0000 P_Y = 0.0000 P_Z = 0.0000 P_R = 0.0000 * * * ******************************************************************************* ******************************************************************************* * * * ANGULAR MOMENTA IN THE THREE CARTESIAN DIRECTIONS * * * ******************************************************************************* * * * SPIN IN X-DIRECTION SPIN IN Y-DIRECTION SPIN IN Z-DIRECTION * * ------------------- ------------------- ------------------- * * INTRI ORBIT TOTAL INTRI ORBIT TOTAL INTRI ORBIT TOTAL * * * * NEUTRONS 0.000 0.000 0.000 0.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 * * ======================= * * 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 * * * ******************************************************************************* * * * 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)= 331.798302 (PRO)= 323.126285 (TOT)= 654.924587 * * SUM EPS: (NEU)= -743.128900 (PRO)= -597.514331 (TOT)= -1340.643231 * * PAIRING: (NEU)= 0.000000 (PRO)= 0.000000 (TOT)= 0.000000 * * P-REARR: (NEU)= 0.000000 (PRO)= 0.000000 (TOT)= 0.000000 * * * * PAIRGAP: (NEU)= 0.000000 (PRO)= 0.000000 THIS SPACE IS EMPTY * * E-FERMI: (NEU)= -13.151954 (PRO)= -6.327659 AND AWAITS YOUR AD * * * * COULOMB: (DIR)= 80.716917 (EXC)= -7.957977 (TOT)= 72.758940 * * COULEXC: (SCA)= -7.955261 (VEC)= -0.002716 (TOT)= -7.957977 * * * * 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= 0.000000 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -342.859309 * * * * SPIN-ORB (EVE)= -1.565943 (ODD)= 0.000000 (TOT)= -1.565943 * * SKYRME: (EVE)= -1070.542835 (ODD)= 0.000000 (TOT)= -1070.542835 * * * * TOTAL: (STAB)= -0.000013 (SP)= -342.859322 (FUN)= -342.859309 * * * ******************************************************************************* ******************************************************************************* * * * ONLY THE ISOSPIN PROJECTION IS PERFORMED * * 8 GAUSS-LEGENDRE KNOTS IS USED * * * ******************************************************************************* ******************************************************************************* * * * THE MATRIX OF KERNELS THAT HAVE ALREADY BEEN CALCULATED IS AS FOLLOWS: * * * ******************************************************************************* * * * GAMMA * * A - * * L |0| * * P |0| * * H |1| * * A - * * 1|0| * * - * * * ******************************************************************************* * * * THE GCM KERNELS IN FUNCTION OF THE EULER ANGLES ARE SAVED; FILE (FILKER): * * ca40_iso.ker * * * ******************************************************************************* * * * N001-L001-R001 1 <= I <= 1 AND 1 <= K <= 1; TOT 1 PNTS ARE SAVD * * * ******************************************************************************* * * * KERNELS AND AVERAGE VALUES OBTAINED FOR ANGULAR-MOMENTUM PROJECTED STATES * * * ******************************************************************************* * * * REAL(KERNEL) IMAG(KERNEL) AVERAGE ENERGY * * * * T= 0 NORM = 0.993940797784 0.000000000000 -342.634058 * * T= 1 NORM = 0.006047755176 0.000000000000 -300.197315 * * T= 2 NORM = 0.000011434975 0.000000000000 -257.946328 * * T= 3 NORM = 0.000000012057 0.000000000000 -215.737825 * * T= 4 NORM = 0.000000000008 0.000000000000 -171.548853 * * T= 5 NORM = 0.000000000000 0.000000000000 -2129.641926 * * * ******************************************************************************* * * * SUM RULES OBTAINED FOR ANGULAR-MOMENTUM PROJECTED STATES VS THE HF VALUES * * * ******************************************************************************* * * * REAL(SUM RULE) IMAG(SUM RULE) HF VALUE * * * * NORM= 1.000000000000 0.000000000000 * * SKYRME= -1070.542835 0.000000 -1070.542835 * * EKIN_T= 654.924587 0.000000 654.924587 * * COUL_D= 80.716917 0.000000 80.716917 * * COUL_E= -7.957977 0.000000 -7.957977 * * COUL_T= 72.758940 0.000000 72.758940 * * * * MULTIPOLE L= 0 M=-5 QMUL_P= 0.000000 0.000000 0.000000 * * L= 0 M=-4 QMUL_P= 0.000000 0.000000 0.000000 * * L= 0 M=-3 QMUL_P= 0.000000 0.000000 0.000000 * * L= 0 M=-2 QMUL_P= 0.000000 0.000000 0.000000 * * L= 0 M=-1 QMUL_P= 0.000000 0.000000 0.000000 * * L= 0 M= 0 QMUL_P= 0.000000 0.000000 20.000000 * * L= 0 M= 1 QMUL_P= 0.000000 0.000000 0.000000 * * L= 0 M= 2 QMUL_P= 0.000000 0.000000 0.000000 * * L= 0 M= 3 QMUL_P= 0.000000 0.000000 0.000000 * * L= 0 M= 4 QMUL_P= 0.000000 0.000000 0.000000 * * L= 0 M= 5 QMUL_P= 0.000000 0.000000 0.000000 * * L= 1 M=-5 QMUL_P= 0.000000 0.000000 0.000000 * * L= 1 M=-4 QMUL_P= 0.000000 0.000000 0.000000 * * L= 1 M=-3 QMUL_P= 0.000000 0.000000 0.000000 * * L= 1 M=-2 QMUL_P= 0.000000 0.000000 0.000000 * * L= 1 M=-1 QMUL_P= 0.000000 0.000000 0.000000 * * L= 1 M= 0 QMUL_P= 0.000000 0.000000 0.000000 * * L= 1 M= 1 QMUL_P= 0.000000 0.000000 0.000000 * * L= 1 M= 2 QMUL_P= 0.000000 0.000000 0.000000 * * L= 1 M= 3 QMUL_P= 0.000000 0.000000 0.000000 * * L= 1 M= 4 QMUL_P= 0.000000 0.000000 0.000000 * * L= 1 M= 5 QMUL_P= 0.000000 0.000000 0.000000 * * L= 2 M=-5 QMUL_P= 0.000000 0.000000 0.000000 * * L= 2 M=-4 QMUL_P= 0.000000 0.000000 0.000000 * * L= 2 M=-3 QMUL_P= 0.000000 0.000000 0.000000 * * L= 2 M=-2 QMUL_P= 0.000000 0.000000 0.000000 * * L= 2 M=-1 QMUL_P= 0.000000 0.000000 0.000000 * * L= 2 M= 0 QMUL_P= 0.000000 0.000000 0.000000 * * L= 2 M= 1 QMUL_P= 0.000000 0.000000 0.000000 * * L= 2 M= 2 QMUL_P= 0.000000 0.000000 0.000000 * * L= 2 M= 3 QMUL_P= 0.000000 0.000000 0.000000 * * L= 2 M= 4 QMUL_P= 0.000000 0.000000 0.000000 * * L= 2 M= 5 QMUL_P= 0.000000 0.000000 0.000000 * * * ******************************************************************************* * * * CUT-OFF "EPSISO" = 0.00000000010000 ==> GOOD-T BASIS OF DIM = 4 * * * ******************************************************************************* * * * ISOSPIN-MIXED EIGENSTATES * * ------------------------- * * * * N EIGENENERGY i T |C_i|^2 Re[C_i] Im[C_i] * * * * 1 -342.860677 1 0 0.994703 0.997348 0.000000 * * 2 1 0.005296 0.072772 0.000000 * * 3 2 0.000002 0.001316 0.000000 * * 4 3 0.000000 -0.000061 0.000000 * * DOMINANT AMPLITUDE SQUARED EQUALS: 0.9947025220 AND CORRESPONDS TO T = 0 * * COULOMB MIXING IN THIS STATE IS: 0.0052974780 [ 0.529748%] * * * * * * 2 -300.253660 1 0 0.005276 0.072636 0.000000 * * 2 1 0.988062 -0.994013 0.000000 * * 3 2 0.006659 -0.081605 0.000000 * * 4 3 0.000003 -0.001655 0.000000 * * DOMINANT AMPLITUDE SQUARED EQUALS: 0.9880618341 AND CORRESPONDS TO T = 1 * * COULOMB MIXING IN THIS STATE IS: 0.0119381659 [ 1.193817%] * * * * * * 3 -258.016938 1 0 0.000021 -0.004628 0.000000 * * 2 1 0.006609 0.081296 0.000000 * * 3 2 0.985082 -0.992513 0.000000 * * 4 3 0.008288 -0.091039 0.000000 * * DOMINANT AMPLITUDE SQUARED EQUALS: 0.9850815887 AND CORRESPONDS TO T = 2 * * COULOMB MIXING IN THIS STATE IS: 0.0149184113 [ 1.491841%] * * * * * * 4 -215.384252 1 0 0.000000 0.000241 0.000000 * * 2 1 0.000033 -0.005784 0.000000 * * 3 2 0.008257 0.090869 0.000000 * * 4 3 0.991709 -0.995846 0.000000 * * DOMINANT AMPLITUDE SQUARED EQUALS: 0.9917092280 AND CORRESPONDS TO T = 3 * * COULOMB MIXING IN THIS STATE IS: 0.0082907720 [ 0.829077%] * * * * * ******************************************************************************* ******************************************************************************* * * * CALLING SUBROUTINE "RECORD" TO SAVE THE RECORD FILE CONTAINING RESULTS * * ca40_iso.rec * * THE MEAN FIELDS OBTAINED IN THE LAST-BUT-ONE ITERATION HAVE BEEN STORED * * * * CALLING SUBROUTINE "RECFIL" TO SAVE THE FIELDS FILE CONTAINING RESULTS * * ca40_iso.fil * * * ******************************************************************************* ******************************************************************************* * * * WARNING! THE STABILITY REQUIRED WAS NOT ACHIEVED. NO-OF-ITERATIONS EXCEEDED * * WARNING! THE STABILITY REQUIRED WAS NOT ACHIEVED. NO-OF-ITERATIONS EXCEEDED * * * ******************************************************************************* * * ******************************************************************************* * * * CALCULATION OF THE DATA SET NO. 3 HAS BEEN SUCCESSFULLY COMPLETED * * * ******************************************************************************* ******************************************************************************* * * * PROCESSOR EXECUTION TIMES IN SUBROUTINES * * INCLUDING THEIR COMPLETE DOWN-CALLING TREES * * * ******************************************************************************* * * * 5565 => HFODD 5043 => COUENE 4983 => GAUOPT 3685 => PROANG * * * * 299 => DENSHF 89 => DIASIZ 88 => DIAMAT 51 => DEFGOG * * * * 35 => COUMAT 35 => INTCOU 33 => INTEGH 13 => MOMETS * * * * 10 => INTSOR 9 => INTMUL 8 => DENMAC 7 => POWALL * * * * 7 => INTMAS 5 => BEGINT 5 => DOBROY 3 => INTCEN * * * * 2 => SPAVER 1 => SKFILD 1 => AVOBSE 1 => RECORD * * * * 1 => AVANGY 1 => MAGMOM 0 => SPIMOM 0 => LINAVR * * * * 0 => ROTAVR 0 => INTSPI 0 => INTKIN 0 => NILABS * * * * 0 => INTORB 0 => BEGINC 0 => NILSON 0 => RECOUL * * * * 0 => AVPARI 0 => SHICOE 0 => PAXMOM 0 => MOMSIF * * * * 0 => SHIMOM 0 => MOMVMU 0 => SHISIF 0 => SHIMAG * * * * 0 => AVSIMP 0 => DIAPRO * * * ******************************************************************************* * * * NUMBERS OF CALLS TO SUBROUTINES * * * ******************************************************************************* * * * 3934 => INTMUL 2880 => BEGINT 1665 => DEFGOG 1110 => GAUOPT * * * * 576 => INTMAS 576 => INTCEN 576 => INTSOR 328 => INTKIN * * * * 237 => DIAMAT 231 => PAXMOM 182 => DENSHF 182 => DENMAC * * * * 150 => DIASIZ 150 => AVPARI 150 => MOMETS 150 => MOMSIF * * * * 150 => MOMVMU 150 => MAGMOM 150 => INTSPI 150 => SPIMOM * * * * 150 => LINAVR 150 => ROTAVR 150 => SHICOE 144 => INTEGH * * * * 111 => COUENE 77 => RECORD 75 => INTORB 75 => SKFILD * * * * 75 => SHIMOM 75 => SHISIF 75 => SHIMAG 66 => SPAVER * * * * 44 => COUMAT 44 => BEGINC 44 => INTCOU 43 => DOBROY * * * * 6 => AVANGY 6 => AVOBSE 6 => AVSIMP 6 => NILABS * * * * 3 => RECOUL 2 => NILSON 1 => HFODD 1 => POWALL * * * * 1 => PROANG 1 => DIAPRO * * * ******************************************************************************* * * * EXECUTION ENDS ON 2011.02.23 AT 15:34:57.564 * * * ******************************************************************************* ALL DONE