******************************************************************************* * * * HFODD HFODD HFODD HFODD HFODD HFODD HFODD HFODD * * * ******************************************************************************* * * * SKYRME-HARTREE-FOCK CODE VERSION 1.75R * * * * ONE SYMMETRY-PLANE AND NO TIME-REVERSAL SYMMETRY * * * * DEFORMED CARTESIAN HARMONIC-OSCILLATOR BASIS * * * ******************************************************************************* * * * JACEK DOBACZEWSKI AND JERZY DUDEK * * * * INSTITUT DE RECHERCHES SUBATOMIQUES * * * * (FORMER CENTRE DE RECHERCHES NUCLEAIRES) * * * * STRASBOURG, 1993-2000 * * * ******************************************************************************* ******************************************************************************* * * * PARAMETER SET SKM*: T0= -2645.00 T1= 410.00 T2= -135.00 T3= 15595.00 * * * * POWER=0.1667 W=130 X0= 0.09000 X1= 0.00000 X2= 0.00000 X3= 0.00000 * * * ******************************************************************************* ******************************************************************************* * * * AUTOMATED HANDLING OF THE COULOMB FILE * * * ******************************************************************************* * * * THE COULOMB FILE WRITTEN ON THE DISC, FILE NAME (FILCOU): * * HFODD.COU * * * ******************************************************************************* * * * NUCLIDE: N = 86 Z = 66 NUMBER OF ITERATIONS = 50 CONTINUATION? = 0 * * * ******************************************************************************* * * * ITERATIONS WILL STOP WHEN THE STABILITY BECOMES LOWER THAN EPSITE=0.000100 * * FOR THE 5 CONSECUTIVE ITERATIONS * * * ******************************************************************************* * * * ITERATIONS WILL STOP WHEN THE PING-PONG DIVERGENCE OCCURS FOR EPSPNG=0.010 * * DURING 3 PAIRS OF CONSECUTIVE ITERATIONS * * * ******************************************************************************* * * * ITERATIONS WILL STOP WHEN THE CHAOTIC-DIVERGENCE CONDITION OCCURS 5 TIMES * * * ******************************************************************************* * * * SLOW-DOWN PARAMETERS: TIME-EVEN FIELDS = 0.50 TIME-ODD FIELDS = 0.50 * * * ******************************************************************************* * * * COULOMB PARAMETERS: NUMCOU = 80 NUMETA = 79 BOUCOU = 20.0 FURMAX =0.25 * * * ******************************************************************************* * * * MAXIMUM NUMBER OF MULTIPOLES CONSIDERED: FOR THE CONSTRAINTS, NMUCON = 2 * * FOR THE COULOMB FIELD, NMUCOU = 4 * * FOR THE OUTPUT INFO, NMUPRI = 4 * * * ******************************************************************************* * * * MULTIPOLE CONSTRAINTS: LAMBDA= 2 MIU= 0 MOMENT= 42.000 STIFFNESS=0.010 * * * ******************************************************************************* * * * PRINTING THE RESULTS FOR THE FOLLOWING ITERATIONS: THE FIRST (1) * * THE MIDDLE (0) * * AND/OR THE LAST (1) * * * ******************************************************************************* * * * NUMBERS OF LEVELS IN TIME-REVERSAL/CHARGE BLOCKS: TIME-UP TIME-DOWN * * NEUTRONS: 86 86 * * PROTONS: 66 66 * * * ******************************************************************************* * * * CALCULATIONS WITHOUT PAIRING * * * ******************************************************************************* * * * CALCULATIONS WITH PARITY/SIGNATURE SYMMETRY * * * ******************************************************************************* * * * CALCULATIONS WITH BROKEN TIME-REVERSAL SYMMETRY * * * ******************************************************************************* * * * INITIAL ROTATIONAL FREQUENCY OMEGAY = 0.500000 * * * ******************************************************************************* * * * LINEAR CONSTRAINT ON SPIN * * * *******************************************************************************
******************************************************************************* * * * PARITY/SIGNATURE CONFIGURATIONS: * * * * V A C U U M P A R T I C L E S H O L E S * * =========== ================= ========= * * (++) (+-) (-+) (--) (++) (+-) (-+) (--) (++) (+-) (-+) (--) * * * * NEUTRONS: 22 22 21 21 0 0 0 0 0 0 0 0 * * PROTONS : 16 16 17 17 0 0 0 0 0 0 0 0 * * * ******************************************************************************* ******************************************************************************* * * * CONVERGENCE REPORT * * * ******************************************************************************* * * * ITER ENERGY STABILITY Q20 GAMMA SPIN OMEGA HOW NICE * * * * 0 -558.029407 -761.683371 54.178 -0.069 98.777 0.500 0.422842 * * 1 -1163.372189 87.725963 48.776 0.047 52.033 0.500 1.081557 * * 2 -1195.980005 63.414957 47.722 -0.025 49.155 0.500 1.055992 * * 3 -1204.228467 38.471958 46.175 0.003 45.243 0.500 1.033002 * * 4 -1207.038466 24.979693 45.544 0.020 45.509 0.500 1.021132 * * 5 -1208.304258 15.158481 44.918 0.038 45.757 0.500 1.012705 * * 6 -1208.924216 9.559571 44.422 0.055 46.049 0.500 1.007971 * * 7 -1209.227387 6.030570 44.007 0.071 46.339 0.500 1.005012 * * 8 -1209.364247 3.884991 43.664 0.086 46.612 0.500 1.003223 * * 9 -1209.413569 2.543565 43.379 0.101 46.858 0.500 1.002108 * * 10 -1209.417195 1.696859 43.142 0.114 47.074 0.500 1.001405 * * 25 -1209.092516 0.016427 42.016 0.215 48.199 0.500 1.000014 * * 26 -1209.085246 0.012900 42.002 0.217 48.214 0.500 1.000011 * * 27 -1209.079025 0.010203 41.990 0.220 48.227 0.500 1.000008 * * 28 -1209.073706 0.008125 41.980 0.222 48.238 0.500 1.000007 * * 29 -1209.069161 0.006513 41.972 0.223 48.248 0.500 1.000005 * * 30 -1209.065279 0.005253 41.965 0.225 48.256 0.500 1.000004 * * 31 -1209.061965 0.004261 41.959 0.226 48.263 0.500 1.000004 * * 32 -1209.059136 0.003474 41.954 0.228 48.268 0.500 1.000003 * * 33 -1209.056722 0.002846 41.950 0.229 48.273 0.500 1.000002 * * 34 -1209.054662 0.002342 41.946 0.230 48.278 0.500 1.000002 * * 35 -1209.052904 0.001934 41.943 0.231 48.281 0.500 1.000002 * * 36 -1209.051404 0.001603 41.941 0.231 48.284 0.500 1.000001 * * 37 -1209.050124 0.001333 41.938 0.232 48.287 0.500 1.000001 * * 38 -1209.049031 0.001112 41.937 0.233 48.289 0.500 1.000001 * * 39 -1209.048098 0.000929 41.935 0.233 48.291 0.500 1.000001 * * 40 -1209.047301 0.000778 41.934 0.234 48.292 0.500 1.000001 * * 41 -1209.046621 0.000653 41.933 0.234 48.294 0.500 1.000001 * * 42 -1209.046040 0.000549 41.932 0.235 48.295 0.500 1.000000 * * 43 -1209.045544 0.000462 41.931 0.235 48.296 0.500 1.000000 * * 44 -1209.045120 0.000389 41.931 0.235 48.297 0.500 1.000000 * * 45 -1209.044757 0.000329 41.930 0.235 48.298 0.500 1.000000 * * 46 -1209.044447 0.000278 41.930 0.236 48.298 0.500 1.000000 * * 47 -1209.044182 0.000235 41.929 0.236 48.299 0.500 1.000000 * * 48 -1209.043955 0.000199 41.929 0.236 48.299 0.500 1.000000 * * 49 -1209.043761 0.000168 41.929 0.236 48.300 0.500 1.000000 * * * ******************************************************************************* ******************************************************************************* * * * DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * DRHO_ = 17.558301 8.908547 17.558301 0.258793 * * DRHOD = 12.481232 6.331119 12.481232 0.181006 * * DLPR_ = -3.718467 -1.885747 -3.718467 -0.053027 * * DTAU_ = 15.576998 7.972052 15.576998 0.367107 * * DSCU_ = 0.119021 0.063593 0.119021 0.008164 * * DDIV_ = 0.821038 0.417814 0.821038 0.014590 * * * * DSPI_ = 0.009448 0.005856 0.009448 0.002265 * * DSPID = 0.006720 0.004158 0.006720 0.001595 * * DLPS_ = -0.004751 -0.006086 -0.004751 -0.007421 * * DCUR_ = 0.058252 0.030270 0.058252 0.002287 * * DKIS_ = 0.014893 0.007474 0.014893 0.000055 * * DROT_ = 0.004337 0.002462 0.004337 0.000587 * * * ******************************************************************************* ******************************************************************************* * * * CONTRIBUTIONS TO ENERGY IN THE SKYRME FUNCTIONAL * * * ******************************************************************************* * * * TOTAL(T) SUM(S) ISOSCALAR(P) ISOVECTOR(M) * * -------- ------ ------------ ------------ * * ERHO_ = -24265.791099 6951.116415 -17415.639545 100.964861 * * ERHOD = 16220.401153 -4113.908512 12165.300865 -58.808223 * * ELPR_ = 317.231713 -64.527908 253.611067 -0.907262 * * ETAU_ = 1070.918626 -543.096073 540.327125 -12.504572 * * ESCU_ = 0.000000 4.332261 4.054159 0.278102 * * EDIV_ = -53.367500 -27.157931 -80.051250 -0.474181 * * ============ ============ ============ ============ * * SUM EVEN: -6710.607107 2206.758253 -4532.397579 28.548725 * * * * ESPI_ = -0.562245 3.872319 2.561340 0.748734 * * ESPID = 0.000000 -2.701576 -2.183390 -0.518186 * * ELPS_ = 0.000000 -0.208253 -0.081291 -0.126962 * * ECUR_ = -4.004811 2.062118 -2.020609 0.077916 * * EKIS_ = 0.000000 -0.509167 -0.507277 -0.001890 * * EROT_ = -0.281919 -0.160046 -0.422878 -0.019087 * * ============ ============ ============ ============ * * SUM ODD: -4.848975 2.355395 -2.654105 0.160525 * * * *******************************************************************************
******************************************************************************* * * * SINGLE-PARTICLE PROPERTIES: HARTREE-FOCK NEUTRONS * * * ******************************************************************************* * * * NO) ENERGY (++,+-,-+,--) | N,nz,lz,OMEG> <P> JY SY GFACT * * * * 75) -11.812 ( 0, 0,20, 0) | 7, 7, 0, 1/2> -100 3.125 0.328 0.105 * * 76) -11.744 ( 0, 0, 0,20) | 5, 3, 2, 3/2> -100 0.297 -0.063 -0.211 * * 77) -11.742 ( 0, 0,21, 0) | 5, 3, 2, 3/2> -100 1.567 0.110 0.071 * * 78) -11.527 (19, 0, 0, 0) | 4, 1, 1, 1/2> 100 0.130 -0.165 -1.268 * * 79) -11.320 ( 0,19, 0, 0) | 4, 1, 3, 5/2> 100 -0.152 -0.004 0.026 * * 80) -11.284 (20, 0, 0, 0) | 4, 1, 3, 5/2> 100 0.169 -0.109 -0.648 * * 81) -11.251 (21, 0, 0, 0) | 6, 5, 1, 1/2> 100 0.771 -0.167 -0.217 * * 82) -11.034 ( 0,20, 0, 0) | 4, 1, 1, 1/2> 100 0.522 -0.004 -0.008 * * 83) -10.899 ( 0,21, 0, 0) | 4, 1, 1, 1/2> 100 0.345 -0.016 -0.046 * * 84) -10.416 (22, 0, 0, 0) | 6, 4, 2, 5/2> 100 -0.014 -0.031 2.177 * * 85) -10.405 ( 0,22, 0, 0) | 6, 4, 2, 5/2> 100 0.074 -0.019 -0.263 * * 86) -9.561 ( 0, 0, 0,21) | 7, 7, 0, 1/2> -100 2.188 -0.036 -0.016 * * 87) -7.830 ( 0, 0,22, 0) | 5, 2, 1, 3/2> -100 0.839 0.124 0.148 * * 88) -7.741 ( 0, 0, 0,22) | 5, 2, 1, 3/2> -100 0.308 0.127 0.412 * * 89) -7.677 (23, 0, 0, 0) | 4, 0, 2, 5/2> 100 -0.282 0.087 -0.308 * * 90) -7.675 ( 0,23, 0, 0) | 4, 0, 2, 5/2> 100 -0.282 0.086 -0.304 * * 91) -7.521 ( 0, 0,23, 0) | 5, 2, 1, 3/2> -100 0.975 0.087 0.089 * * 92) -7.215 ( 0, 0, 0,23) | 7, 7, 0, 1/2> -100 0.026 -0.124 -4.839 * * 93) -7.208 ( 0, 0, 0,24) | 5, 1, 4, 9/2> -100 -0.334 0.008 -0.024 * * 94) -7.208 ( 0, 0,24, 0) | 5, 1, 4, 9/2> -100 -0.336 0.009 -0.026 * * 95) -6.936 ( 0,24, 0, 0) | 6, 4, 0, 1/2> 100 1.271 0.325 0.256 * * 96) -6.423 (24, 0, 0, 0) | 6, 4, 0, 1/2> 100 0.341 -0.415 -1.220 * * 97) -6.364 ( 0, 0, 0,25) | 5, 2, 1, 1/2> -100 0.258 -0.215 -0.832 * * 98) -6.263 (25, 0, 0, 0) | 6, 3, 3, 7/2> 100 -0.015 0.009 -0.602 * * 99) -6.263 ( 0,25, 0, 0) | 6, 3, 3, 7/2> 100 -0.009 0.019 -2.112 * * 100) -6.164 ( 0,26, 0, 0) | 4, 0, 0, 1/2> 100 0.141 0.448 3.168 * * 101) -5.819 ( 0, 0,25, 0) | 5, 2, 1, 1/2> -100 -0.045 -0.064 1.442 * * 102) -5.703 (26, 0, 0, 0) | 4, 0, 0, 1/2> 100 -0.575 -0.399 0.694 * * 103) -5.505 ( 0, 0,26, 0) | 5, 2, 3, 5/2> -100 0.095 -0.021 -0.223 * * 104) -5.503 ( 0, 0, 0,26) | 5, 2, 3, 5/2> -100 0.074 -0.012 -0.161 * * 105) -5.369 ( 0,27, 0, 0) | 6, 4, 2, 3/2> 100 0.350 0.004 0.011 * * 106) -5.366 (27, 0, 0, 0) | 6, 4, 2, 3/2> 100 0.608 -0.019 -0.032 * * 107) -5.172 ( 0,28, 0, 0) | 4, 0, 2, 3/2> 100 -0.519 -0.021 0.041 * * 108) -5.089 (28, 0, 0, 0) | 4, 0, 2, 3/2> 100 -0.445 -0.152 0.341 * * 109) -4.928 ( 0,29, 0, 0) | 4, 0, 4, 7/2> 100 -0.387 -0.031 0.080 * * 110) -4.926 (29, 0, 0, 0) | 4, 0, 4, 7/2> 100 -0.370 -0.040 0.109 * * 111) -4.793 ( 0,30, 0, 0) | 8, 8, 0, 1/2> 100 5.270 0.469 0.089 * * 112) -4.385 ( 0, 0, 0,27) | 7, 6, 1, 1/2> -100 1.575 -0.105 -0.067 * * 113) -4.348 ( 0, 0,27, 0) | 7, 7, 0, 1/2> -100 1.733 0.235 0.136 * * 114) -3.967 ( 0, 0,28, 0) | 7, 5, 2, 5/2> -100 0.234 -0.038 -0.161 * * 115) -3.927 ( 0, 0, 0,28) | 7, 5, 2, 5/2> -100 -0.157 -0.155 0.985 * * 116) -2.555 (30, 0, 0, 0) | 8, 7, 1, 3/2> 100 3.047 0.077 0.025 * * 117) -2.480 ( 0, 0, 0,29) | 5, 1, 2, 5/2> -100 -0.074 0.067 -0.915 * * 118) -2.478 ( 0, 0,29, 0) | 5, 1, 2, 5/2> -100 -0.067 0.069 -1.033 * * 119) -2.222 ( 0, 0, 0,30) | 5, 0, 5,11/2> -100 -0.602 -0.002 0.003 * * * ******************************************************************************* ******************************************************************************* * * * MULTIPOLE MOMENTS IN UNITS OF (10 FERMI)**LAMBDA TOTAL * * * ******************************************************************************* * * * Q00 =152.0000 ............. ............. ............. ............. * * * * Q10 = ZERO Q11 = ZERO ............. ............. ............. * * * * Q20 = 41.9283 Q21 = ZERO Q22 = 0.1730 ............. ............. * * * * Q30 = ZERO Q31 = ZERO Q32 = ZERO Q33 = ZERO ............. * * * * Q40 = 4.8357 Q41 = ZERO Q42 = 0.0074 Q43 = ZERO Q44 =-9.2E-04 * * * ******************************************************************************* ******************************************************************************* * * * ROOT-MEAN-SQUARE AND GEOMETRIC SIZES IN FERMIS TOTAL * * * * R_RMS = 5.5449 X_RMS = 2.3841 Y_RMS = 2.3703 Z_RMS = 4.4095 * * * * R_GEO = 7.1584 X_GEO = 5.3310 Y_GEO = 5.3001 Z_GEO = 9.8599 * * * ******************************************************************************* ******************************************************************************* * * * BOHR DEFORMATION PARAMETERS (FIRST-ORDER APPROXIMATION) TOTAL * * * ******************************************************************************* * * * B10 = ZERO B11 = ZERO ............. ............. ............. * * * * B20 = 0.7112 B21 = ZERO B22 = 0.0029 ............. ............. * * * * B30 = ZERO B31 = ZERO B32 = ZERO B33 = ZERO ............. * * * * B40 = 0.5075 B41 = ZERO B42 = 0.0011 B43 = ZERO B44 =-1.4E-04 * * * *******************************************************************************
******************************************************************************* * * * ANGULAR MOMENTA AND THE FIRST MOMENTS OF INERTIA FOR OMEGA = 0.500000 MEV * * * ******************************************************************************* * * * SPINS J(1) * * --------------------------- --------------------------- * * ORBITAL INTRINSIC TOTAL ORBITAL INTRINSIC TOTAL * * * * NEUTRONS 27.72295 0.90622 28.62917 55.44591 1.81243 57.25834 * * PROTONS 19.01301 0.65769 19.67069 38.02601 1.31537 39.34138 * * -------- * * TOTAL 46.73596 1.56390 48.29986 93.47192 3.12780 96.59972 * * * ******************************************************************************* ******************************************************************************* * * * NEUTRON CONFIGURATIONS * * ======================= * * P S 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 * * --- -------------------------------------------------------------- * * * * CONF: + + 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 * * VACC: + + 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 * * * * CONF: + - 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 * * VACC: + - 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - + 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - + 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - - 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * PROTON CONFIGURATIONS * * ======================= * * P S 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 * * --- -------------------------------------------------------------- * * * * CONF: + + 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + + 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: + - 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 * * VACC: + - 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - + 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 * * VACC: - + 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 * * * * CONF: - - 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 * * VACC: - - 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 * * * ******************************************************************************* * * * ENERGIES (MEV) * * * ******************************************************************************* * * * KINETIC: (NEU)= 1651.194254 (PRO)= 1107.180948 (TOT)= 2758.375202 * * SUM EPS: (NEU)= -2033.727026 (PRO)= -1132.746046 (TOT)= -3166.473072 * * PAIRING: (NEU)= 0.000000 (PRO)= 0.000000 (TOT)= 0.000000 * * * * COULOMB: (DIR)= 564.274457 (EXC)= -25.350821 (TOT)= 538.923637 * * * * CONSTR. (MULT)= 0.000051 SLOPE= 0.001433 CORR.= -0.030046 * * CONSTR. (SPIN)= -24.149930 SLOPE= 0.500000 CORR.= -12.074965 * * * * REARRANGEMENT ENERGY FROM THE SKYRME DENSITY-DEPENDENT TERMS= 1008.649255 * * ROUTHIAN (TOTAL ENERGY PLUS MULTIPOLE AND SPIN CONSTRAINTS)= -1233.193474 * * * * SPIN-ORB (EVE)= -80.525431 (ODD)= -0.441965 (TOT)= -80.967396 * * SKYRME: (EVE)= -4503.848854 (ODD)= -2.493580 (TOT)= -4506.342434 * * * * TOTAL: (STAB)= 0.000143 (SP)= -1209.043452 (FUN)= -1209.043595 * * * ******************************************************************************* ******************************************************************************* * * * REAL-CLOCK EXECUTION TIMES IN SUBROUTINES * * INCLUDING THEIR COMPLETE DOWN-CALLING TREES * * * ******************************************************************************* * * * 1023 => HFODD 0 => GEOMFC 21 => POWALL 0 => RECOUL * * * * 1 => NILSON 3 => INTKIN 180 => DIASIG 174 => DIAMAT * * * * 3 => AVPARI 96 => NILABS 64 => AVIMRE 37 => INTMUL * * * * 261 => DENSHF 64 => DENMAT 37 => MOMETS 50 => BEGINT * * * * 94 => INTCOU 16 => SKFILD 7 => RECORD 153 => INTEGH * * * * 24 => INTMAS 12 => INTCEN 45 => INTSOR 8 => ANGYSP * * * *******************************************************************************