5.16 ME (OP,GOSI)M<strong>and</strong>atory suboption of OP,GOSIThis suboption of OP,GOSI is used both to input <strong>and</strong> to catalog the starting set of matrix elements aswell as to set constraints on the variation of these matrix elements for the least-squares search procedure.The suboption LEVE must immediately precede ME since the catalog of the matrix elements is performedusing the level scheme <strong>and</strong> state indices assigned using the LEVE comm<strong>and</strong>. Up to 500 matrix elements areallowed.Although similar in many respects, the input to the OP,GOSI version of ME is more extensive thanthat required by the OP,COUL version of ME (Section 5-15). The OP,GOSI version of ME differs from theOP,COUL version in the following respects:(a)Restrictions are placed on the range over which each matrix element is allowed to vary during theleast-squares search procedure. These restrictions are used to prevent the code from finding unphysicalsolutions. Moreover, these restrictions limit the range of coupling coefficients ζ over which the ζ-dependence of the q-parameters is fitted.(b) Specifications are given defining which matrix elements are to be treated as free variables, <strong>and</strong>which are to be kept fixed or varied conserving a preset coupling with other matrix elements. Thisallows a reduction in the number of unknowns by using other knowledge such as lifetime, branchingratio or multiple mixing-ratio <strong>data</strong> to restrict the number of free parameters used in the least-squaressearch. Note that lifetimes, branching ratios, multiple mixing ratios <strong>and</strong> E2 matrix elements alsocan be included explicitly in the <strong>data</strong> set used for the least-squares search (See OP,YIEL, Section5-29). Restrictions on matrix elements can be overridden by the OP,RE,A , OP,RE,C <strong>and</strong> OP,RE,Foptions or, conversely, additional restrictions can be imposed using the FIX or LCK comm<strong>and</strong> ofthe suboption CONT without changing the ME input file. Input from the OP,GOSI version of MErequires a five-entry record.MEA summary of the input format is as follows:λ 1 0, 0, 0, 0Specify multipolarity.INDEX1, ±INDEX2, ME, R 1 , R 2upper (R 2 ) limits.List of matrix elements for multipolarity λ 1 plus lower (R 1 )<strong>and</strong>··λ 2 , 0, 0, 0, 0Specify multipolarity.INDEX1, ±INDEX2, ME, R 1 , R 2 List matrix elements for multipolarity λ 2 .··0, 0, 0, 0, 0 - Terminates input.Matrix elements for each multipolarity are input as a set preceded by a single record defining the multipolarity,i.e., λ, 0, 0, 0, 0 where λ= 1through 6 from E1 through E6 respectively, while λ =7for M1 <strong>and</strong>λ =8for M2.Matrix elements are input as:INDEX1, ±INDEX2, ME, R 1 , R 2input.INDEX# is the user-given state number as definedintheLEVE90
ME =< INDEX2||E(M)λ||INDEX1 > The multipole matrix element defined by equation 2.16. Itis given in units of e.barns 1/2 for Eλ matrix elements <strong>and</strong> μ N .barns (1−1)/2 for Mλ matrix elements.The sign assigned INDEX2 plays no role in the definition of the matrix element. A negative signsignifies coupled matrix elements as discussed below.The lower <strong>and</strong> upper limits,respectively, between which the given matrix element ME is allowedto vary. Obviously R 2 >R 1 . Equality of R 1 <strong>and</strong> R 2 implies that this given matrix element is kept fixed atthe value ME. Note that in this case, i.e.,R 1 = R 2 , R1 need not be equal to the ME. For example record1, 2, 0.5, 2, 2 is equivalent to 1, 2, 0.5, 0.5, 0.5Fixing matrix elements as in the first example is recommended because of two features of the code. First,when an OP,RE comm<strong>and</strong> is used the limits of fixed matrix elements are set as R 2 = |R 2 | <strong>and</strong> R 1 = −|R 1 |.If fixing is done as in the second example the matrix element will be allowed to vary only in one direction,i.e., between ±0.5. Secondly, when fitting the q-parameters, the ζ-ranges are set according to the actuallimits R 1 ,R 2 .Iffixed matrix elements are released at a later stage during the <strong>analysis</strong>, then there is less riskof incorrect extrapolation if the approach used in the first example is employed rather than later extendingthe limits without recalculating the q-parameter maps. Note that neither R 1 nor R 2 should be exactly zero,use a small number instead.A negative sign assigned to INDEX2 specifies that the matrix element defined by the pair of indicesINDEX1, INDEX2 is not a free variable but is a coupled one. In this case R 1 <strong>and</strong> R 2 are no longer upper<strong>and</strong> lower limits but the pair of indices of the matrix element to which this matrix element is coupled. Thecode automatically assigns upper <strong>and</strong> lower limits to the coupled matrix elements using the upper <strong>and</strong> lowerlimits given to the one to which it is coupled while preserving the ratio of the initial values of the coupledmatrix elements. These upper <strong>and</strong> lower limits calculated by the code are used by the code if this couplingis subsequently released. The ratio of coupled matrix elements set by the initial values is preserved in theleast squares search if the coupling is not released.]As an example consider the pair of matrix elements:1, 2, 2.0, −4, 42, −3, 1., 1, 2The second statement signifies that the matrix element connecting state 2 to 3 is coupled to the matrixelement connecting states 1 <strong>and</strong> 2 in the ratio:M(2 + 3)M(1 + 2) =+0.5This ratio will be preserved in the least squares search if this coupling is not released. Limits (R 1 ,R 2 )of−2, +2 will be assigned by the code to the matrix element connecting states 2 <strong>and</strong> 3. These limits are usefulif the coupling of these matrix elements is subsequently released.The convention of coupling matrix elements already described is valid only within a single multipolarity.The coupling of matrix elements belonging to different multipolarities is performed by using 100λ +R 2 asinput for the index R 2 of the “slave“ matrix element where λ is the multipolarity of the master matrixelement. The convention that λ= 7for M1 <strong>and</strong> λ= 8for M2 is still valid.An important restriction is that there can be only one “master“ matrix element with a number of slavescoupled to it if a group of matrix elements are specified to be mutually related. For example, a valid sequenceis1, 2, 2., −4, 4 E2 set of matrix elements2, −3, 2., 1, 2·2, −3, 0.5, 1, 202 M1 set of matrix elements91
- Page 1:
COULOMB EXCITATION DATA ANALYSIS CO
- Page 4 and 5:
10 MINIMIZATION BY SIMULATED ANNEAL
- Page 6 and 7:
1 INTRODUCTION1.1 Gosia suite of Co
- Page 8 and 9:
104 Ru, 110 Pd, 165 Ho, 166 Er, 186
- Page 13 and 14:
Figure 1: Coordinate system used to
- Page 15 and 16:
Cλ E =1.116547 · (13.889122) λ (
- Page 17 and 18:
Figure 2: The orbital integrals R 2
- Page 19 and 20:
2.2 Gamma Decay Following Electroma
- Page 21 and 22:
where :d 2 σ= σ R (θ p ) X R kχ
- Page 23 and 24:
Formula 2.49 is valid only for t mu
- Page 25 and 26:
à XK(α) =exp−iτ i (E γ )x i (
- Page 27 and 28:
important to have an accurate knowl
- Page 29 and 30:
3 APPROXIMATE EVALUATION OF EXCITAT
- Page 31 and 32:
with the reduced matrix element M c
- Page 33 and 34:
q (20)s (0 + → 2 + ) · M 1 ζ (2
- Page 35 and 36:
esults of minimization and error ru
- Page 37 and 38:
adjustment of the stepsize accordin
- Page 39 and 40: approximation reliability improves
- Page 41 and 42: Zd 2 σ(I → I f )Y (I → I f )=s
- Page 43 and 44: 4.5 MinimizationThe minimization, i
- Page 45 and 46: X(CC k Yk c − Yk e ) 2 /σ 2 k =m
- Page 47 and 48: However, estimation of the stepsize
- Page 49 and 50: It can be shown that as long as the
- Page 51 and 52: een exceeded; third, the user-given
- Page 53 and 54: where f k stands for the functional
- Page 55 and 56: x i + δx i Rx iexp ¡ − 1 2 χ2
- Page 57 and 58: method used for the minimization, i
- Page 59 and 60: OP,ERRO (ERRORS) (5.6):Activates th
- Page 61 and 62: -----OP,SIXJ (SIX-j SYMBOL) (5.25):
- Page 63 and 64: 5.3 CONT (CONTROL)This suboption of
- Page 65 and 66: I,I1 Ranges of matrix elements to b
- Page 67 and 68: CODE DEFAULT OTHER CONSEQUENCES OF
- Page 69 and 70: 5.4 OP,CORR (CORRECT )This executio
- Page 71 and 72: 5.6 OP,ERRO (ERRORS)ThemoduleofGOSI
- Page 73 and 74: 5.7 OP,EXIT (EXIT)This option cause
- Page 75 and 76: M AControls the number of magnetic
- Page 77 and 78: 5.10 OP,GDET (GE DETECTORS)This opt
- Page 79 and 80: 5.12 OP,INTG (INTEGRATE)This comman
- Page 81 and 82: ¡ dE¢dx1 ..¡ dEdx¢Stopping powe
- Page 83 and 84: NI1, NI2 Number of subdivisions of
- Page 85 and 86: 5.13 LEVE (LEVELS)Mandatory subopti
- Page 87 and 88: 5.15 ME (OP,COUL)Mandatory suboptio
- Page 89: Figure 10: Model system having 4 st
- Page 93 and 94: When entering matrix elements in th
- Page 95 and 96: There are no restrictions concernin
- Page 97 and 98: 5.18 OP,POIN (POINT CALCULATION)Thi
- Page 99 and 100: 5.20 OP,RAW (RAW UNCORRECTED γ YIE
- Page 101 and 102: 5.21 OP,RE,A (RELEASE,A)This option
- Page 103 and 104: 5.25 OP,SIXJ (SIXJ SYMBOL)This stan
- Page 105 and 106: 5.27 OP,THEO (COLLECTIVE MODEL ME)C
- Page 107 and 108: 2,5,1,-2,23,5,1,-2,23,6,1,-2,2Matri
- Page 109 and 110: 5.29 OP,TROU (TROUBLE)This troubles
- Page 111 and 112: to that of the previous experiment,
- Page 113 and 114: To reduce the unnecessary input, on
- Page 115 and 116: OP,STAR or OP,POIN under OP,GOSI. N
- Page 117 and 118: 5.31 INPUT OF EXPERIMENTAL γ-RAY Y
- Page 119 and 120: 6 QUADRUPOLE ROTATION INVARIANTS -
- Page 121 and 122: *½P 5 (J) = s(E2 × E2) J ׯh¾
- Page 123 and 124: The expectation value of cos3δ can
- Page 125 and 126: where ē is an arbitratry vector. D
- Page 127 and 128: achieved using “mixed“ calculat
- Page 129 and 130: TAPE9 Contains the parameters neede
- Page 131 and 132: TAPE18 Input file, containing the i
- Page 133 and 134: 7.4.4 CALCULATION OF THE INTEGRATED
- Page 135 and 136: OP,EXITInput: TAPE4,TAPE7,TAPE9Outp
- Page 137 and 138: OP,ERRO0,MS,MEND,1,0,RMAXand the fi
- Page 139 and 140: 8 SIMULTANEOUS COULOMB EXCITATION:
- Page 141 and 142:
4, 3, 1kr88.corKr corrected yields
- Page 143 and 144:
0 Correction for in-flight decay ch
- Page 145 and 146:
OP, ERRO Estimation of errors of fi
- Page 147 and 148:
9 COULOMB EXCITATION OF ISOMERIC ST
- Page 149 and 150:
configurations with a probability e
- Page 151 and 152:
The average range covered by each m
- Page 153 and 154:
SFX,NTOTI1(1),I2(1),RSIGN(1)I1(2),I
- Page 155 and 156:
11.2 LearningtoWriteGosiaInputsThe
- Page 157 and 158:
(1.6 MeV)1.1 MeV0.75 MeV0.4 MeV0.08
- Page 159 and 160:
Define the germaniumdetector geomet
- Page 161 and 162:
Figure 15: Flow diagram for Gosia m
- Page 163 and 164:
gosia < 2-make-correction-factors.i
- Page 165 and 166:
Issue the commandgosia < 9-diag-err
- Page 167 and 168:
At this point, it is suggested to c
- Page 169 and 170:
calculation.) In this case, a copy
- Page 171 and 172:
4,-4, -3.705, 3,44,5, 4.626, 3.,7.5
- Page 173 and 174:
90145901459014590145901459014590145
- Page 175 and 176:
.10.028921.10.026031.10.023431.10.0
- Page 177 and 178:
5,5,634,650,82.000,84.000634,638,64
- Page 179 and 180:
***********************************
- Page 181 and 182:
*** CHISQ= 0.134003E+01 ***MATRIX E
- Page 183 and 184:
CALCULATED AND EXPERIMENTAL YIELDS
- Page 185 and 186:
11.7 Annotated excerpt from a Coulo
- Page 187 and 188:
11.8 Accuracy and speed of calculat
- Page 189 and 190:
18,10.056,0.068,0.082,0.1,0.12,0.15
- Page 191 and 192:
line 152 Eu 182 Tanumber (keV) (keV
- Page 193 and 194:
1.6 Normalization between data sets
- Page 195 and 196:
13 GOSIA 2007 RELEASE NOTESThese no
- Page 197 and 198:
Matrix elements 500(April 1990, T.
- Page 199 and 200:
14 GOSIA Manual UpdatesDATE UPDATE2
- Page 201 and 202:
[KIB08]T.Kibédi,T.W.Burrows,M.B.Tr