Geant4-11
G4AtomicShells_XDB_EADL.cc
Go to the documentation of this file.
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26//
28// //
29// File: G4AtomicShells_XDB_EADL.cc //
30// Author: D.H. Wright (SLAC) //
31// Date: 9 August 2018 //
32// //
33// Description: Class containing number of shells, electron configurations //
34// and binding energies for atoms from Z = 1 to Z = 120. //
35// Most entries are taken from the X-ray Data Book, with //
36// unmeasured values supplied by Kibedi. See documentation //
37// file in G4EMLOW7.3 or later, directory fluor/ //
38// //
40
42#include "G4SystemOfUnits.hh"
43
44const G4int
46{
47 1 , // nonexistent zero element
48
49 1, 1, 2, 2, 3, 3, 4, 4, 4, 4, // 1 - 10
50
51 5, 5, 6, 6, 7, 7, 7, 7, 8, 8, // 11 - 20
52
53 9, 9, 9, 10, 10, 10, 10, 10, 10, 10, // 21 - 30
54
55 11, 11, 12, 12, 12, 12, 13, 13, 14, 14, // 31 - 40
56
57 14, 15, 15, 15, 15, 14, 15, 15, 16, 16, // 41 - 50
58
59 17, 17, 17, 17, 18, 18, 19, 21, 21, 20, // 51 - 60
60
61 19, 20, 20, 21, 21, 20, 20, 20, 20, 20, // 61 - 70
62
63 21, 21, 21, 21, 22, 22, 22, 22, 22, 22, // 71 - 80
64
65 23, 23, 24, 24, 24, 24, 25, 25, 26, 26, // 81 - 90
66
67 27, 27, 27, 26, 27, 28, 28, 27, 27, 27, // 91 - 100
68
69 27, 27, 28, 28, 28, 28, 29, 29, 29, 29, // 101 - 110
70
71 29, 29, 30, 30, 31, 31, 31, 31, 32, 32
72};
73
74
75const G4int
77{
78 0, // nonexistent zero element
79
80 1, 2, 3, 5, 7, 10, 13, 17, 21, 25, // 1 - 10
81
82 29, 34, 39, 45, 51, 58, 65, 72, 79, 87, // 11 - 20
83
84 95, 104, 113, 122, 132, 142, 152, 162, 172, 182, // 21 - 30
85
86 192, 203, 214, 226, 238, 250, 262, 275, 288, 302, // 31 - 40
87
88 316, 330, 345, 360, 375, 390, 404, 419, 434, 450, // 41 - 50
89
90 466, 483, 500, 517, 534, 552, 570, 589, 610, 631, // 51 - 60
91
92 651, 670, 690, 710, 731, 752, 772, 792, 812, 832, // 61 - 70
93
94 852, 873, 894, 915, 936, 958, 980, 1002, 1024, 1046, // 71 - 80
95
96 1068, 1091, 1114, 1138, 1162, 1186, 1210, 1235, 1260, 1286, // 81 - 90
97
98 1312, 1339, 1366, 1393, 1419, 1446, 1474, 1502, 1529, 1556, // 91 - 100
99
100 1583, 1610, 1637, 1665, 1693, 1721, 1749, 1778, 1807, 1836, // 101 - 110
101
102 1865, 1894, 1923, 1953, 1983, 2014, 2045, 2076, 2107, 2139
103};
104
105
106const G4double
108{
109 // 0 ---------------------------------------------------------------------
110 0.0000, // nonexistent zero element
111
112 // H ---------------------------------------------------------------------
113 0.0136,
114
115 // He ---------------------------------------------------------------------
116 0.0246,
117
118 // Li ---------------------------------------------------------------------
119 0.0547, 0.0054,
120
121 // Be ---------------------------------------------------------------------
122 0.1115, 0.0093,
123
124 // B ---------------------------------------------------------------------
125 0.1880, 0.0126, 0.0047,
126
127 // C ---------------------------------------------------------------------
128 0.2842, 0.0181, 0.0090,
129
130 // N ---------------------------------------------------------------------
131 0.4099, 0.0373, 0.0089, 0.0097,
132
133 // O ---------------------------------------------------------------------
134 0.5431, 0.0416, 0.0068, 0.0074,
135
136 // F ---------------------------------------------------------------------
137 0.6967, 0.0340, 0.0084, 0.0087,
138
139 // Ne ---------------------------------------------------------------------
140 0.8702, 0.0485, 0.0217, 0.0216,
141
142 // Na ---------------------------------------------------------------------
143 1.0708, 0.0635, 0.03065, 0.03081, 0.0007,
144
145 // Mg ---------------------------------------------------------------------
146 1.3030, 0.0887, 0.04978, 0.0495, 0.0021,
147
148 // Al ---------------------------------------------------------------------
149 1.5596, 0.1178, 0.07295, 0.07255, 0.0007, 0.0060,
150
151 // Si ---------------------------------------------------------------------
152 1.8390, 0.1497, 0.09982, 0.09942, 0.0076, 0.0030,
153
154 // P ---------------------------------------------------------------------
155 2.1455, 0.1890, 0.1360, 0.1350, 0.0162, 0.0096, 0.0101,
156
157 // S ---------------------------------------------------------------------
158 2.4720, 0.2309, 0.1636, 0.1625, 0.0158, 0.0078, 0.0082,
159
160 // Cl ---------------------------------------------------------------------
161 2.8224, 0.2700, 0.2020, 0.2000, 0.0175, 0.0067, 0.0067,
162
163 // Ar ---------------------------------------------------------------------
164 3.2059, 0.3263, 0.2506, 0.2484, 0.0293, 0.0159, 0.0157,
165
166 // K ---------------------------------------------------------------------
167 3.6084, 0.3786, 0.2973, 0.2946, 0.0348, 0.0183, 0.0183,
168 0.0001,
169
170 // Ca ---------------------------------------------------------------------
171 4.0385, 0.4384, 0.3497, 0.3462, 0.0443, 0.0254, 0.0254,
172 0.0018,
173
174 // Sc ---------------------------------------------------------------------
175 4.4920, 0.4980, 0.4036, 0.3987, 0.0511, 0.0283, 0.0283,
176 0.0066, 0.0017,
177
178 // Ti ---------------------------------------------------------------------
179 4.9660, 0.5609, 0.4602, 0.4538, 0.0587, 0.0326, 0.0326,
180 0.0037, 0.0016,
181
182 // V ---------------------------------------------------------------------
183 5.4650, 0.6267, 0.5198, 0.5121, 0.0663, 0.0372, 0.0372,
184 0.0022, 0.0017,
185
186 // Cr ---------------------------------------------------------------------
187 5.9890, 0.6960, 0.5838, 0.5741, 0.0741, 0.0422, 0.0422,
188 0.0029, 0.0022, 0.0010,
189
190 // Mn ---------------------------------------------------------------------
191 6.5390, 0.7691, 0.6499, 0.6387, 0.0823, 0.0472, 0.0472,
192 0.0035, 0.0027, 0.0019,
193
194 // Fe ---------------------------------------------------------------------
195 7.1120, 0.8446, 0.7199, 0.7068, 0.0913, 0.0527, 0.0527,
196 0.0039, 0.0031, 0.0021,
197
198 // Co ---------------------------------------------------------------------
199 7.7090, 0.9251, 0.7932, 0.7781, 0.1010, 0.0589, 0.0599,
200 0.0033, 0.0027, 0.0019,
201
202 // Ni ---------------------------------------------------------------------
203 8.3330, 1.0086, 0.8700, 0.8527, 0.1108, 0.0680, 0.0662,
204 0.0039, 0.0033, 0.0022,
205
206 // Cu ---------------------------------------------------------------------
207 8.9790, 1.0967, 0.9523, 0.9327, 0.1225, 0.0773, 0.0751,
208 0.0018, 0.0015, 0.0012,
209
210 // Zn ---------------------------------------------------------------------
211 9.6590, 1.1962, 1.0449, 1.0218, 0.1398, 0.0914, 0.0886,
212 0.0102, 0.0101, 0.0013,
213
214 // Ga ---------------------------------------------------------------------
215 10.3670, 1.2990, 1.1432, 1.1164, 0.1595, 0.1035, 0.1000,
216 0.0187, 0.0187, 0.0056, 0.0008,
217
218 // Ge ---------------------------------------------------------------------
219 11.1030, 1.4146, 1.2481, 1.2170, 0.1801, 0.1249, 0.1208,
220 0.0298, 0.0292, 0.0090, 0.0023,
221
222 // As ---------------------------------------------------------------------
223 11.8670, 1.5270, 1.3591, 1.3236, 0.2047, 0.1462, 0.1412,
224 0.0417, 0.0417, 0.0125, 0.0025, 0.0025,
225
226 // Se ---------------------------------------------------------------------
227 12.6580, 1.6520, 1.4743, 1.4339, 0.2296, 0.1665, 0.1607,
228 0.0555, 0.0546, 0.0162, 0.0056, 0.0056,
229
230 // Br ---------------------------------------------------------------------
231 13.4740, 1.7820, 1.5960, 1.5500, 0.2570, 0.1890, 0.1820,
232 0.0700, 0.0690, 0.0273, 0.0052, 0.0046,
233
234 // Kr ---------------------------------------------------------------------
235 14.3260, 1.9210, 1.7309, 1.6784, 0.2928, 0.2222, 0.2144,
236 0.0950, 0.0938, 0.0275, 0.0141, 0.0141,
237
238 // Rb ---------------------------------------------------------------------
239 15.2000, 2.0650, 1.8640, 1.8040, 0.3267, 0.2487, 0.2391,
240 0.1130, 0.1120, 0.0305, 0.0163, 0.0153, 0.0040,
241
242 // Sr ---------------------------------------------------------------------
243 16.1050, 2.2160, 2.0070, 1.9400, 0.3587, 0.2803, 0.2700,
244 0.1360, 0.1342, 0.0389, 0.0213, 0.0201, 0.0050,
245
246 // Y ---------------------------------------------------------------------
247 17.0380, 2.3730, 2.1560, 2.0800, 0.3920, 0.3106, 0.2988,
248 0.1577, 0.1558, 0.0438, 0.0244, 0.0231, 0.0024, 0.0060,
249
250 // Zr ---------------------------------------------------------------------
251 17.9980, 2.5320, 2.3070, 2.2230, 0.4303, 0.3435, 0.3298,
252 0.1811, 0.1788, 0.0506, 0.0285, 0.0271, 0.0030, 0.0070,
253
254 // Nb ---------------------------------------------------------------------
255 18.9860, 2.6980, 2.4650, 2.3710, 0.4666, 0.3761, 0.3606,
256 0.2050, 0.2023, 0.0564, 0.0326, 0.0308, 0.0032, 0.0070,
257
258 // Mo ---------------------------------------------------------------------
259 20.0000, 2.8660, 2.6250, 2.5200, 0.5063, 0.4116, 0.3940,
260 0.2311, 0.2279, 0.0632, 0.0376, 0.0355, 0.0019, 0.0012,
261 0.0070,
262
263 // Tc ---------------------------------------------------------------------
264 21.0440, 3.0430, 2.7930, 2.6770, 0.5440, 0.4476, 0.4177,
265 0.2576, 0.2539, 0.0695, 0.0423, 0.0399, 0.0020, 0.0012,
266 0.0070,
267
268 // Ru ---------------------------------------------------------------------
269 22.1170, 3.2240, 2.9670, 2.8380, 0.5861, 0.4835, 0.4614,
270 0.2842, 0.2800, 0.0750, 0.0463, 0.0432, 0.0024, 0.0018,
271 0.0070,
272
273 // Rh ---------------------------------------------------------------------
274 23.2200, 3.4120, 3.1460, 3.0040, 0.6281, 0.5213, 0.4965,
275 0.3119, 0.3072, 0.0814, 0.0505, 0.0473, 0.0028, 0.0022,
276 0.0070,
277
278 // Pd ---------------------------------------------------------------------
279 24.3500, 3.6040, 3.3300, 3.1730, 0.6716, 0.5599, 0.5323,
280 0.3405, 0.3352, 0.0871, 0.0557, 0.0509, 0.0017, 0.0013,
281
282 // Ag ---------------------------------------------------------------------
283 25.5140, 3.8060, 3.5240, 3.3510, 0.7190, 0.6038, 0.5730,
284 0.3740, 0.3683, 0.0970, 0.0637, 0.0583, 0.0036, 0.0031,
285 0.0080,
286
287 // Cd ---------------------------------------------------------------------
288 26.7110, 4.0180, 3.7270, 3.5380, 0.7720, 0.6526, 0.6184,
289 0.4119, 0.4052, 0.1098, 0.0639, 0.0639, 0.0117, 0.0107,
290 0.0090,
291
292 // In ---------------------------------------------------------------------
293 27.9400, 4.2380, 3.9380, 3.7300, 0.8272, 0.7032, 0.6653,
294 0.4514, 0.4439, 0.1229, 0.0735, 0.0735, 0.0177, 0.0169,
295 0.0001, 0.0008,
296
297 // Sn ---------------------------------------------------------------------
298 29.2000, 4.4650, 4.1560, 3.9290, 0.8847, 0.7565, 0.7146,
299 0.4932, 0.4849, 0.1371, 0.0836, 0.0836, 0.0249, 0.0239,
300 0.0009, 0.0011,
301
302 // Sb ---------------------------------------------------------------------
303 30.4910, 4.6980, 4.3800, 4.1320, 0.9460, 0.8127, 0.7664,
304 0.5375, 0.5282, 0.1532, 0.0956, 0.0956, 0.0333, 0.0321,
305 0.0067, 0.0022, 0.0020,
306
307 // Te ---------------------------------------------------------------------
308 31.8140, 4.9390, 4.6120, 4.3410, 1.0060, 0.8708, 0.8200,
309 0.5834, 0.5730, 0.1694, 0.1033, 0.1033, 0.0419, 0.0404,
310 0.0116, 0.0026, 0.0020,
311
312 // I ---------------------------------------------------------------------
313 33.1690, 5.1880, 4.8520, 4.5570, 1.0720, 0.9310, 0.8750,
314 0.6308, 0.6193, 0.1860, 0.1230, 0.1230, 0.0506, 0.0489,
315 0.0136, 0.0038, 0.0029,
316
317 // Xe ---------------------------------------------------------------------
318 34.5610, 5.4530, 5.1070, 4.7860, 1.1487, 1.0021, 0.9406,
319 0.6890, 0.6764, 0.2132, 0.1467, 0.1455, 0.0695, 0.0675,
320 0.0233, 0.0134, 0.0121,
321
322 // Cs ---------------------------------------------------------------------
323 35.9850, 5.7140, 5.3590, 5.0120, 1.2110, 1.0710, 1.0030,
324 0.7405, 0.7266, 0.2323, 0.1724, 0.1613, 0.0798, 0.0775,
325 0.0227, 0.0142, 0.0121, 0.004,
326
327 // Ba ---------------------------------------------------------------------
328 37.4410, 5.9890, 5.6240, 5.2470, 1.2930, 1.1370, 1.0630,
329 0.7957, 0.7805, 0.2535, 0.1920, 0.1786, 0.0926, 0.0899,
330 0.0303, 0.0170, 0.0148, 0.0050,
331
332 // La ---------------------------------------------------------------------
333 38.9250, 6.2660, 5.8910, 5.4830, 1.3620, 1.2090, 1.1280,
334 0.8530, 0.8360, 0.2747, 0.2058, 0.1960, 0.1053, 0.1025,
335 0.0343, 0.0193, 0.0168, 0.0060, 0.0050,
336
337 // Ce ---------------------------------------------------------------------
338 40.4430, 6.5490, 6.1640, 5.7230, 1.4360, 1.2740, 1.1870,
339 0.9024, 0.8838, 0.2910, 0.2232, 0.2065, 0.1090, 0.1076,
340 0.0001, 0.0001, 0.0378, 0.0198, 0.0170, 0.0060, 0.0050,
341
342 // Pr ---------------------------------------------------------------------
343 41.9910, 6.8350, 6.4400, 5.9640, 1.5110, 1.3370, 1.2429,
344 0.9483, 0.9288, 0.3045, 0.2363, 0.2176, 0.1151, 0.1151,
345 0.0020, 0.0020, 0.0374, 0.0223, 0.0223, 0.0060, 0.0050,
346
347 // Nd ---------------------------------------------------------------------
348 43.5690, 7.1260, 6.7220, 6.2080, 1.5750, 1.4030, 1.2970,
349 1.0033, 0.9804, 0.3192, 0.2433, 0.2246, 0.1205, 0.1205,
350 0.0015, 0.0015, 0.0375, 0.0211, 0.0211, 0.0050,
351
352 // Pm ---------------------------------------------------------------------
353 45.1840, 7.4280, 7.0130, 6.4590, 1.6486, 1.4710, 1.3570,
354 1.0520, 1.0270, 0.3314, 0.2420, 0.2420, 0.1200, 0.1200,
355 0.0035, 0.0360, 0.0245, 0.0201, 0.0050,
356
357 // Sm ---------------------------------------------------------------------
358 46.8340, 7.7370, 7.3120, 6.7160, 1.7230, 1.5410, 1.4200,
359 1.1109, 1.0834, 0.3472, 0.2656, 0.2474, 0.1290, 0.1290,
360 0.0052, 0.0052, 0.0374, 0.0213, 0.0213, 0.0050,
361
362 // Eu ---------------------------------------------------------------------
363 48.5190, 8.0520, 7.6170, 6.9770, 1.8000, 1.6140, 1.4810,
364 1.1586, 1.1275, 0.3600, 0.2840, 0.2570, 0.1330, 0.1277,
365 0.0000, 0.0000, 0.0320, 0.0220, 0.0220, 0.0060,
366
367 // Gd ---------------------------------------------------------------------
368 50.2390, 8.3760, 7.9300, 7.2430, 1.8810, 1.6880, 1.5440,
369 1.2219, 1.1896, 0.3786, 0.2860, 0.2710, 0.1495, 0.1426,
370 0.0086, 0.0086, 0.0360, 0.0280, 0.0210, 0.0060, 0.0060,
371
372 // Tb ---------------------------------------------------------------------
373 51.9960, 8.7080, 8.2520, 7.5140, 1.9680, 1.7680, 1.6110,
374 1.2769, 1.2411, 0.3960, 0.3224, 0.2841, 0.1505, 0.1505,
375 0.0077, 0.0024, 0.0456, 0.0287, 0.0226, 0.0060, 0.0060,
376
377 // Dy ---------------------------------------------------------------------
378 53.7890, 9.0460, 8.5810, 7.7900, 2.0470, 1.8420, 1.6760,
379 1.3330, 1.2926, 0.4142, 0.3335, 0.2932, 0.1536, 0.1536,
380 0.0080, 0.0043, 0.0499, 0.0263, 0.0263, 0.0060,
381
382 // Ho ---------------------------------------------------------------------
383 55.6180, 9.3940, 8.9180, 8.0710, 2.1280, 1.9230, 1.7410,
384 1.3920, 1.3510, 0.4324, 0.3435, 0.3082, 0.1600, 0.1600,
385 0.0086, 0.0052, 0.0493, 0.0308, 0.0241, 0.0060,
386
387 // Er ---------------------------------------------------------------------
388 57.4860, 9.7510, 9.2640, 8.3580, 2.2070, 2.0060, 1.8120,
389 1.4530, 1.4090, 0.4498, 0.3662, 0.3202, 0.1676, 0.1676,
390 0.0053, 0.0047, 0.0506, 0.0314, 0.0247, 0.0060,
391
392 // Tm ---------------------------------------------------------------------
393 59.3900, 10.1160, 9.6170, 8.6480, 2.3070, 2.0900, 1.8850,
394 1.5150, 1.4680, 0.4709, 0.3859, 0.3326, 0.1755, 0.1755,
395 0.0062, 0.0046, 0.0547, 0.0318, 0.0250, 0.0060,
396
397 // Yb ---------------------------------------------------------------------
398 61.3320, 10.4860, 9.9780, 8.9440, 2.3980, 2.1730, 1.9500,
399 1.5760, 1.5280, 0.4805, 0.3887, 0.3397, 0.1912, 0.1824,
400 0.0025, 0.0013, 0.0520, 0.0303, 0.0241, 0.0060,
401
402 // Lu ---------------------------------------------------------------------
403 63.3140, 10.8700, 10.3490, 9.2440, 2.4910, 2.2640, 2.0240,
404 1.6390, 1.5890, 0.5068, 0.4124, 0.3592, 0.2061, 0.1963,
405 0.0089, 0.0075, 0.0573, 0.0336, 0.0267, 0.0046, 0.0070,
406
407 // Hf ---------------------------------------------------------------------
408 65.3510, 11.2710, 10.7390, 9.5610, 2.6010, 2.3650, 2.1080,
409 1.7160, 1.6620, 0.5380, 0.4382, 0.3807, 0.2200, 0.2115,
410 0.0159, 0.0142, 0.0642, 0.0380, 0.0299, 0.0066, 0.0070,
411
412 // Ta ---------------------------------------------------------------------
413 67.4160, 11.6820, 11.1360, 9.8810, 2.7080, 2.4690, 2.1940,
414 1.7930, 1.7350, 0.5634, 0.4634, 0.4009, 0.2379, 0.2264,
415 0.0235, 0.0216, 0.0697, 0.0422, 0.0327, 0.0057, 0.0080,
416
417 // W ---------------------------------------------------------------------
418 69.5250, 12.1000, 11.5440, 10.2070, 2.8200, 2.5750, 2.2810,
419 1.8720, 1.8090, 0.5941, 0.4904, 0.4236, 0.2559, 0.2435,
420 0.0336, 0.0314, 0.0756, 0.0453, 0.0368, 0.0061, 0.0080,
421
422 // Re ---------------------------------------------------------------------
423 71.6760, 12.5270, 11.9590, 10.5350, 2.9320, 2.6820, 2.3670,
424 1.9490, 1.8830, 0.6254, 0.5187, 0.4468, 0.2739, 0.2605,
425 0.0429, 0.0405, 0.0830, 0.0456, 0.0346, 0.0038, 0.0025,
426 0.0080,
427
428 // Os ---------------------------------------------------------------------
429 73.8710, 12.9680, 12.3850, 10.8710, 3.0490, 2.7920, 2.4570,
430 2.0310, 1.9600, 0.6582, 0.5491, 0.4707, 0.2931, 0.2785,
431 0.0534, 0.0507, 0.0840, 0.0580, 0.0445, 0.0009, 0.0004,
432 0.0080,
433
434 // Ir ---------------------------------------------------------------------
435 76.1110, 13.4190, 12.8240, 11.2150, 3.1740, 2.9090, 2.5510,
436 2.1160, 2.0400, 0.6911, 0.5778, 0.4958, 0.3119, 0.2963,
437 0.0638, 0.0608, 0.0952, 0.0630, 0.0480, 0.0042, 0.0032,
438 0.0090,
439
440 // Pt ---------------------------------------------------------------------
441 78.3950, 13.8800, 13.2730, 11.5640, 3.2960, 3.0270, 2.6450,
442 2.2020, 2.1220, 0.7254, 0.6091, 0.5194, 0.3316, 0.3146,
443 0.0745, 0.0712, 0.1017, 0.0653, 0.0517, 0.0028, 0.0014,
444 0.0090,
445
446 // Au ---------------------------------------------------------------------
447 80.7250, 14.3530, 13.7340, 11.9190, 3.4250, 3.1480, 2.7430,
448 2.2910, 2.2060, 0.7621, 0.6427, 0.5463, 0.3532, 0.3351,
449 0.0876, 0.0840, 0.1072, 0.0742, 0.0572, 0.0033, 0.0018,
450 0.0090,
451
452 // Hg ---------------------------------------------------------------------
453 83.1020, 14.8390, 14.2090, 12.2840, 3.5620, 3.2790, 2.8470,
454 2.3850, 2.2950, 0.8022, 0.6802, 0.5766, 0.3782, 0.3588,
455 0.1040, 0.0999, 0.1270, 0.0831, 0.0645, 0.0096, 0.0078,
456 0.0090,
457
458 // Tl --------------------------------------------------------------------
459 85.5300, 15.3470, 14.6980, 12.6580, 3.7040, 3.4160, 2.9570,
460 2.4850, 2.3890, 0.8462, 0.7205, 0.6095, 0.4057, 0.3850,
461 0.1222, 0.1178, 0.1360, 0.0946, 0.0735, 0.0147, 0.0125,
462 0.0060, 0.0008,
463
464 // Pb --------------------------------------------------------------------
465 88.0050, 15.8610, 15.200, 13.0350, 3.8510, 3.5540, 3.0660,
466 2.5860, 2.4840, 0.8918, 0.7619, 0.6435, 0.4343, 0.4122,
467 0.1417, 0.1369, 0.1470, 0.1064, 0.0833, 0.0207, 0.0181,
468 0.0031, 0.0007,
469
470 // Bi --------------------------------------------------------------------
471 90.5240, 16.3880, 15.7110, 13.4190, 3.9990, 3.6960, 3.1770,
472 2.6880, 2.5800, 0.9390, 0.8052, 0.6788, 0.4640, 0.4401,
473 0.1623, 0.1570, 0.1593, 0.1190, 0.0926, 0.0269, 0.0238,
474 0.0075, 0.0012, 0.0002,
475
476 // Po --------------------------------------------------------------------
477 93.1050, 16.9390, 16.2440, 13.8140, 4.1490, 3.8540, 3.3020,
478 2.7980, 2.6830, 0.9950, 0.8510, 0.7050, 0.5000, 0.4730,
479 0.1840, 0.1840, 0.1770, 0.1320, 0.1040, 0.0310, 0.0310,
480 0.0110, 0.0032, 0.0014,
481
482 // At --------------------------------------------------------------------
483 95.7300, 17.4930, 16.7850, 14.2140, 4.3170, 4.0080, 3.4260,
484 2.9090, 2.7870, 1.0420, 0.8860, 0.7400, 0.5330, 0.5070,
485 0.2100, 0.2100, 0.1950, 0.1480, 0.1150, 0.0400, 0.0400,
486 0.0150, 0.0057, 0.0028,
487
488 // Rn --------------------------------------------------------------------
489 98.4040, 18.0490, 17.3370, 14.6190, 4.4820, 4.1590, 3.5380,
490 3.0220, 2.8920, 1.0970, 0.9290, 0.7680, 0.5670, 0.5410,
491 0.2380, 0.2380, 0.2140, 0.1640, 0.1270, 0.0480, 0.0480,
492 0.0260, 0.0076, 0.0041,
493
494 // Fr --------------------------------------------------------------------
495 101.1370, 18.6390, 17.9070, 15.0310, 4.6520, 4.3270, 3.6630,
496 3.1360, 3.0000, 1.1530, 0.9800, 0.8100, 0.6030, 0.5770,
497 0.2680, 0.2680, 0.2340, 0.1820, 0.1400, 0.0580, 0.0580,
498 0.0340, 0.0150, 0.0150, 0.0040,
499
500 // Ra --------------------------------------------------------------------
501 103.9220, 19.2370, 18.4840, 15.4440, 4.8220, 4.4900, 3.7920,
502 3.2480, 3.1050, 1.2080, 1.0580, 0.8790, 0.6360, 0.6030,
503 0.2990, 0.2990, 0.2540, 0.2000, 0.1530, 0.0680, 0.0680,
504 0.0440, 0.0190, 0.0190, 0.0050,
505
506 // Ac --------------------------------------------------------------------
507 106.7550, 19.8400, 19.0830, 15.8710, 5.0020, 4.6560, 3.9090,
508 3.3700, 3.2190, 1.2690, 1.0800, 0.8900, 0.6750, 0.6390,
509 0.3190, 0.3190, 0.2720, 0.2150, 0.1670, 0.0800, 0.0800,
510 0.0398, 0.0241, 0.0170, 0.0060, 0.0060,
511
512 // Th --------------------------------------------------------------------
513 109.6510, 20.4720, 19.6930, 16.3000, 5.1820, 4.8300, 4.0460,
514 3.4910, 3.3320, 1.3300, 1.1680, 0.9664, 0.7121, 0.6752,
515 0.3424, 0.3331, 0.2900, 0.2290, 0.1820, 0.0925, 0.0854,
516 0.0414, 0.0245, 0.0166, 0.0060, 0.0060,
517
518 // Pa --------------------------------------------------------------------
519 112.6010, 21.1050, 20.3140, 16.7330, 5.3670, 5.0010, 4.1740,
520 3.6110, 3.4420, 1.3870, 1.2240, 1.0070, 0.7430, 0.7080,
521 0.3710, 0.3600, 0.3100, 0.2320, 0.2320, 0.0940, 0.0940,
522 0.0060, 0.0467, 0.0281, 0.0189, 0.0060, 0.0060,
523
524 // U --------------------------------------------------------------------
525 115.6060, 21.7570, 20.9480, 17.1660, 5.5480, 5.1820, 4.3030,
526 3.7280, 3.5520, 1.4390, 1.2710, 1.0430, 0.7783, 0.7362,
527 0.3882, 0.3774, 0.3210, 0.2570, 0.1920, 0.1028, 0.0942,
528 0.0060, 0.0439, 0.0268, 0.0168, 0.0060, 0.0060,
529
530 // Np --------------------------------------------------------------------
531 118.6690, 22.4281, 21.6013, 17.6106, 5.7396, 5.3667, 4.4334,
532 3.8498, 3.6652, 1.5001, 1.3277, 1.0860, 0.8161, 0.7708,
533 0.4143, 0.4034, 0.3398, 0.2834, 0.2062, 0.1093, 0.1012,
534 0.0060, 0.0500, 0.0293, 0.0175, 0.0060, 0.0060,
535
536 // Pu --------------------------------------------------------------------
537 121.7910, 23.0978, 22.2668, 18.0570, 5.9329, 5.5412, 4.5566,
538 3.9726, 3.7781, 1.5593, 1.3774, 1.1209, 0.8489, 0.8015,
539 0.4374, 0.4252, 0.3519, 0.2825, 0.2153, 0.1160, 0.1052,
540 0.0060, 0.0486, 0.0306, 0.0184, 0.0060,
541
542 // Am --------------------------------------------------------------------
543 124.9810, 23.8051, 22.9494, 18.5062, 6.1326, 5.7470, 4.7060,
544 4.0921, 3.8869, 1.6192, 1.4351, 1.1680, 0.8804, 0.8300,
545 0.4633, 0.4490, 0.3730, 0.3030, 0.2164, 0.1180, 0.1079,
546 0.0060, 0.0060, 0.0504, 0.0311, 0.0181, 0.0060,
547
548 // Cm --------------------------------------------------------------------
549 128.2000, 24.5200, 23.6600, 18.9700, 6.3350, 5.9400, 4.8500,
550 4.2300, 4.0150, 1.6850, 1.4950, 1.2050, 0.9250, 0.8700,
551 0.4950, 0.4850, 0.3970, 0.2350, 0.1350, 0.1200, 0.0920,
552 0.0060, 0.0060, 0.0670, 0.0500, 0.0060, 0.0060, 0.0970,
553
554 // Bk --------------------------------------------------------------------
555 131.6000, 25.2600, 24.3800, 19.4400, 6.5450, 6.1450, 4.9900,
556 4.3600, 4.1350, 1.7450, 1.5550, 1.2450, 0.9600, 0.9000,
557 0.5250, 0.5100, 0.4150, 0.3300, 0.2430, 0.1400, 0.1250,
558 0.0060, 0.0060, 0.0960, 0.0700, 0.0500, 0.0044, 0.0060,
559
560 // Cf --------------------------------------------------------------------
561 135.0000, 26.0100, 25.1200, 19.9100, 6.7550, 6.3500, 5.1350,
562 4.4950, 4.2600, 1.8100, 1.6150, 1.2850, 0.9950, 0.9300,
563 0.5550, 0.5400, 0.4300, 0.3450, 0.2500, 0.1450, 0.1300,
564 0.0060, 0.0060, 0.1000, 0.0700, 0.0500, 0.006,
565
566 // Es --------------------------------------------------------------------
567 138.4000, 26.7800, 25.8800, 20.3900, 6.9750, 6.5600, 5.2800,
568 4.6300, 4.3850, 1.8750, 1.6750, 1.3250, 1.0300, 0.9650,
569 0.5850, 0.5700, 0.4450, 0.3600, 0.2600, 0.1550, 0.1350,
570 0.0060, 0.0060, 0.1000, 0.0700, 0.0500, 0.0060,
571
572 // Fm --------------------------------------------------------------------
573 142.0000, 27.5700, 26.6600, 20.8700, 7.2000, 6.7800, 5.4250,
574 4.7650, 4.5100, 1.9400, 1.7400, 1.3700, 1.0650, 1.0000,
575 0.6150, 0.6000, 0.4600, 0.3750, 0.2700, 0.1650, 0.1400,
576 0.0060, 0.0060, 0.1000, 0.0800, 0.0500, 0.0060,
577
578 // Md --------------------------------------------------------------------
579 145.6000, 28.3800, 27.4600, 21.3600, 7.4300, 7.0000, 5.5750,
580 4.9050, 4.6330, 2.0050, 1.8050, 1.4150, 1.1050, 1.0350,
581 0.6450, 0.6300, 0.4800, 0.3900, 0.2800, 0.1750, 0.1450,
582 0.0070, 0.0070, 0.1100, 0.0800, 0.0500, 0.0060,
583
584 // No --------------------------------------------------------------------
585 149.3000, 29.2220, 28.2800, 21.8500, 7.6700, 7.2300, 5.7300,
586 5.0500, 4.7650, 2.0750, 1.8750, 1.4600, 1.1450, 1.0700,
587 0.6800, 0.6500, 0.5000, 0.4100, 0.2900, 0.1850, 0.1550,
588 0.0100, 0.0100, 0.1100, 0.0800, 0.0500, 0.0060,
589
590 // Lr --------------------------------------------------------------------
591 153.0000, 30.0800, 29.1300, 22.3500, 7.9200, 7.4700, 5.8900,
592 5.2000, 4.9000, 2.1500, 1.9500, 1.5050, 1.1900, 1.1100,
593 0.7200, 0.6900, 0.5200, 0.4300, 0.3000, 0.2000, 0.1600,
594 0.0150, 0.0150, 0.1200, 0.0900, 0.0600, 0.0070, 0.0070,
595
596 // Rf --------------------------------------------------------------------
597 156.9000, 30.9600, 30.0000, 22.8600, 8.1700, 7.7200, 6.0550,
598 5.3500, 5.0400, 2.2250, 2.0250, 1.5550, 1.2350, 1.1500,
599 0.7600, 0.7200, 0.5400, 0.4550, 0.3150, 0.2150, 0.1700,
600 0.0200, 0.0200, 0.1200, 0.0900, 0.0600, 0.0080, 0.0080,
601
602 // Db --------------------------------------------------------------------
603 160.8000, 31.8700, 30.9000, 23.3800, 8.4300, 7.9800, 6.6250,
604 5.5050, 5.1850, 2.3050, 2.1050, 1.6050, 1.2800, 1.1900,
605 0.8000, 0.7500, 0.5600, 0.4800, 0.3300, 0.2300, 0.1800,
606 0.0250, 0.0250, 0.1300, 0.1000, 0.0600, 0.0090, 0.0080,
607
608 // Sg --------------------------------------------------------------------
609 164.9000, 32.8000, 31.8200, 23.9100, 8.7000, 8.2400, 6.4000,
610 5.6650, 5.3300, 2.3850, 2.1900, 1.6600, 1.3250, 1.2350,
611 0.8400, 0.7800, 0.5850, 0.5050, 0.3450, 0.2450, 0.1900,
612 0.0300, 0.0300, 0.1400, 0.1100, 0.0600, 0.0090, 0.0080,
613
614 // Bh --------------------------------------------------------------------
615 169.0000, 33.7600, 32.7700, 24.4400, 8.9800, 8.5100, 6.5800,
616 5.8300, 5.4750, 2.4700, 2.2750, 1.7150, 1.3700, 1.2800,
617 0.8800, 0.8200, 0.6100, 0.5300, 0.3600, 0.2600, 0.2000,
618 0.0400, 0.0400, 0.1400, 0.1100, 0.0700, 0.0100, 0.0100,
619 0.0090,
620
621 // Hs --------------------------------------------------------------------
622 173.2000, 34.7500, 33.7400, 24.9200, 9.2700, 8.7900, 6.7600,
623 5.9950, 5.6200, 2.5550, 2.3650, 1.7700, 1.4200, 1.3250,
624 0.9200, 0.8600, 0.6350, 0.5550, 0.3750, 0.2750, 0.2100,
625 0.0500, 0.0500, 0.1500, 0.1200, 0.0700, 0.0100, 0.0100,
626 0.0090,
627
628 // Mt --------------------------------------------------------------------
629 177.6000, 35.7700, 34.7400, 25.5200, 9.5700, 9.0800, 6.9400,
630 6.1600, 5.7700, 2.6450, 2.4600, 1.8250, 1.4700, 1.3700,
631 0.9600, 0.9000, 0.6600, 0.5800, 0.3900, 0.2900, 0.2200,
632 0.0600, 0.0600, 0.1600, 0.1200, 0.0700, 0.0100, 0.0100,
633 0.0100,
634
635 // Ds --------------------------------------------------------------------
636 182.0000, 36.8100, 35.7700, 26.0700, 9.8700, 9.3700, 7.1200,
637 6.3300, 5.9200, 2.7350, 2.5600, 1.8800, 1.5200, 1.4150,
638 1.0000, 0.9400, 0.6900, 0.6100, 0.4050, 0.3050, 0.2300,
639 0.0700, 0.0700, 0.1600, 0.1300, 0.0800, 0.0100, 0.0100,
640 0.0100,
641
642 // Rg --------------------------------------------------------------------
643 186.1360, 37.8810, 36.7990, 26.6480, 10.1870, 9.6870, 7.2240,
644 6.4470, 6.0250, 2.9260, 2.6730, 1.9780, 1.5970, 1.4830,
645 0.9604, 0.9267, 0.7283, 0.6450, 0.2810, 0.2510, 0.2445,
646 0.0671, 0.0614, 0.1267, 0.0901, 0.0498, 0.0094, 0.0066,
647 0.0094,
648
649 // Cp --------------------------------------------------------------------
650 190.7170, 38.9840, 37.8960, 27.2130, 10.5080, 10.0020, 7.4060,
651 6.6180, 6.1780, 3.0330, 2.7760, 2.0390, 1.6520, 1.5320,
652 1.0016, 0.9662, 0.7624, 0.6770, 0.2960, 0.2650, 0.2577,
653 0.0753, 0.0690, 0.1349, 0.0969, 0.0528, 0.0104, 0.0073,
654 0.0098,
655
656 // Uut --------------------------------------------------------------------
657 195.4100, 40.1210, 39.0340, 27.7880, 10.8420, 10.3310, 7.5930,
658 6.7930, 6.3350, 3.1460, 2.8860, 2.1070, 1.7120, 1.5870,
659 1.0464, 1.0092, 0.8006, 0.7140, 0.3150, 0.2820, 0.2740,
660 0.0866, 0.0798, 0.1464, 0.1069, 0.0586, 0.0138, 0.0101,
661 0.0116, 0.0049,
662
663 // Fl --------------------------------------------------------------------
664 200.2140, 41.2940, 40.2070, 28.3700, 11.1840, 10.6720, 7.7830,
665 6.9720, 6.4950, 3.2640, 2.9990, 2.1740, 1.7740, 1.6420,
666 1.0923, 1.0533, 0.8409, 0.7530, 0.3350, 0.3000, 0.2910,
667 0.0984, 0.0911, 0.1583, 0.1174, 0.0646, 0.0173, 0.0130,
668 0.0135, 0.0057,
669
670 // Uup --------------------------------------------------------------------
671 205.1360, 42.5010, 41.4200, 28.9590, 11.5410, 11.0220, 7.9790,
672 7.1550, 6.6580, 3.3860, 3.1180, 2.2450, 1.8380, 1.7010,
673 1.1410, 1.1000, 0.8831, 0.7940, 0.3570, 0.3200, 0.3100,
674 0.1124, 0.1045, 0.1726, 0.1301, 0.0724, 0.0227, 0.0177,
675 0.0169, 0.0078, 0.0038,
676
677 // Lv --------------------------------------------------------------------
678 210.1770, 43.7500, 42.6740, 29.5560, 11.9070, 11.3870, 8.1760,
679 7.3410, 6.8240, 3.5130, 3.2410, 2.3160, 1.9030, 1.7590,
680 1.1907, 1.1477, 0.9275, 0.8380, 0.3790, 0.3400, 0.3296,
681 0.1269, 0.1184, 0.1874, 0.1434, 0.0805, 0.0282, 0.0226,
682 0.0203, 0.0101, 0.0044,
683
684 // Uus --------------------------------------------------------------------
685 215.3370, 44.0310, 43.9750, 30.1580, 12.2830, 11.7620, 8.3760,
686 7.5300, 6.9910, 3.6430, 3.4690, 2.3900, 1.9700, 1.8190,
687 1.2413, 1.1962, 0.9733, 0.8820, 0.4020, 0.3610, 0.3496,
688 0.1417, 0.1327, 0.2026, 0.1571, 0.0886, 0.0337, 0.0275,
689 0.0237, 0.0124, 0.0052,
690
691 // Uuo --------------------------------------------------------------------
692 220.6300, 46.3540, 45.3540, 30.7670, 12.6700, 12.1510, 8.5790,
693 7.7200, 7.1610, 3.7780, 3.5000, 2.4640, 2.0370, 1.8800,
694 1.2928, 1.2455, 1.0205, 0.9280, 0.4250, 0.3820, 0.3699,
695 0.1570, 0.1472, 0.2183, 0.1714, 0.0968, 0.0394, 0.0324,
696 0.0272, 0.0149, 0.0060,
697
698 // Uue --------------------------------------------------------------------
699 226.0630, 47.7240, 46.7140, 31.3860, 13.0740, 12.5560, 8.7890,
700 7.9180, 7.3360, 3.9200, 3.6390, 2.5430, 2.1090, 1.9440,
701 1.3476, 1.2982, 1.0720, 0.9800, 0.4520, 0.4060, 0.3931,
702 0.1751, 0.1646, 0.2371, 0.1887, 0.1075, 0.0476, 0.0399,
703 0.0331, 0.0197, 0.0084, 0.0036,
704
705 // Ubn --------------------------------------------------------------------
706 231.6340, 49.1390, 48.1580, 32.0130, 13.4890, 12.9750, 9.0000,
707 8.1170, 7.5130, 4.0660, 3.7840, 2.6220, 2.1820, 2.0100,
708 1.4040, 1.3523, 1.1256, 1.0330, 0.4790, 0.4310, 0.4173,
709 0.1941, 0.1830, 0.2570, 0.2071, 0.1190, 0.0565, 0.0479,
710 0.0396, 0.0251, 0.0113, 0.0042
711};
712
713const G4int
715{
716 0, // nonexistent zero element
717
718 // H ---------------------------------------------------------
719 1,
720
721 // He ---------------------------------------------------------
722 2,
723
724 // Li --------------------------------------------------------
725 2, 1,
726
727 // Be --------------------------------------------------------
728 2, 2,
729
730 // B --------------------------------------------------------
731 2, 2, 1,
732
733 // C ---------------------------------------------------------
734 2, 2, 2,
735
736 // N ---------------------------------------------------------
737 2, 2, 2, 1,
738
739 // O ---------------------------------------------------------
740 2, 2, 2, 2,
741
742 // F ---------------------------------------------------------
743 2, 2, 2, 3,
744
745 // Ne ---------------------------------------------------------
746 2, 2, 2, 4,
747
748 // Na ----------------------------------------------
749 2, 2, 2, 4, 1,
750
751 // Mg ----------------------------------------------
752 2, 2, 2, 4, 2,
753
754 // Al ----------------------------------------------
755 2, 2, 2, 4, 2, 1,
756
757 // Si ----------------------------------------------
758 2, 2, 2, 4, 2, 2,
759
760 // P ----------------------------------------------
761 2, 2, 2, 4, 2, 2, 1,
762
763 // S ----------------------------------------------
764 2, 2, 2, 4, 2, 2, 2,
765
766 // Cl ----------------------------------------------
767 2, 2, 2, 4, 2, 2, 3,
768
769 // Ar ----------------------------------------------
770 2, 2, 2, 4, 2, 2, 4,
771
772 // K ----------------------------------------------
773 2, 2, 2, 4, 2, 2, 4, 1,
774
775 // Ca ----------------------------------------------
776 2, 2, 2, 4, 2, 2, 4, 2,
777
778 // Sc ----------------------------------------------
779 2, 2, 2, 4, 2, 2, 4, 1, 2,
780
781 // Ti ----------------------------------------------
782 2, 2, 2, 4, 2, 2, 4, 2, 2,
783
784 // V ----------------------------------------------
785 2, 2, 2, 4, 2, 2, 4, 3, 2,
786
787 // Cr ----------------------------------------------
788 2, 2, 2, 4, 2, 2, 4, 4, 1, 1,
789
790 // Mn ----------------------------------------------
791 2, 2, 2, 4, 2, 2, 4, 4, 1, 2,
792
793 // Fe ----------------------------------------------
794 2, 2, 2, 4, 2, 2, 4, 4, 2, 2,
795
796 // Co ----------------------------------------------
797 2, 2, 2, 4, 2, 2, 4, 4, 3, 2,
798
799 // Ni ----------------------------------------------
800 2, 2, 2, 4, 2, 2, 4, 4, 4, 2,
801
802 // Cu ----------------------------------------------
803 2, 2, 2, 4, 2, 2, 4, 4, 6, 1,
804
805 // Zn ----------------------------------------------
806 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
807
808 // Ga ----------------------------------------------
809 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
810 1,
811
812 // Ge ----------------------------------------------
813 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
814 2,
815
816 // As ----------------------------------------------
817 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
818 2, 1,
819
820 // Se ----------------------------------------------
821 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
822 2, 2,
823
824 // Br ----------------------------------------------
825 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
826 2, 3,
827
828 // Kr ----------------------------------------------
829 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
830 2, 4,
831
832 // Rb ----------------------------------------------
833 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
834 2, 4, 1,
835
836 // Sr ----------------------------------------------
837 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
838 2, 4, 2,
839
840 // Y ----------------------------------------------
841 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
842 2, 4, 1, 2,
843
844 // Zr ----------------------------------------------
845 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
846 2, 4, 2, 2,
847
848 // Nb ----------------------------------------------
849 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
850 2, 4, 4, 1,
851
852 // Mo ----------------------------------------------
853 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
854 2, 4, 4, 1, 1,
855
856 // Tc ----------------------------------------------
857 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
858 2, 4, 4, 1, 2,
859
860 // Ru ----------------------------------------------
861 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
862 2, 4, 4, 3, 1,
863
864 // Rh ----------------------------------------------
865 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
866 2, 4, 4, 4, 1,
867
868 // Pd ----------------------------------------------
869 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
870 2, 4, 4, 6,
871
872 // Ag ----------------------------------------------
873 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
874 2, 4, 4, 6, 1,
875
876 // Cd ----------------------------------------------
877 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
878 2, 4, 4, 6, 2,
879
880 // In ----------------------------------------------
881 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
882 2, 4, 4, 6, 2, 1,
883
884 // Sn ----------------------------------------------
885 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
886 2, 4, 4, 6, 2, 2,
887
888 // Sb ----------------------------------------------
889 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
890 2, 4, 4, 6, 2, 2, 1,
891
892 // Te ----------------------------------------------
893 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
894 2, 4, 4, 6, 2, 2, 2,
895
896 // I ----------------------------------------------
897 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
898 2, 4, 4, 6, 2, 2, 3,
899
900 // Xe ----------------------------------------------
901 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
902 2, 4, 4, 6, 2, 2, 4,
903
904 // Cs ----------------------------------------------
905 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
906 2, 4, 4, 6, 2, 2, 4, 1,
907
908 // Ba ----------------------------------------------
909 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
910 2, 4, 4, 6, 2, 2, 4, 2,
911
912 // La ----------------------------------------------
913 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
914 2, 4, 4, 6, 2, 2, 4, 1, 2,
915
916 // Ce ----------------------------------------------
917 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
918 2, 4, 4, 6, 1, 1, 2, 2, 4, 1,
919 1,
920
921 // Pr ----------------------------------------------
922 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
923 2, 4, 4, 6, 2, 1, 2, 2, 4, 1,
924 1,
925
926 // Nd ----------------------------------------------
927 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
928 2, 4, 4, 6, 4, 1, 2, 2, 4, 1,
929
930 // Pm ----------------------------------------------
931 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
932 2, 4, 4, 6, 5, 2, 2, 4, 2,
933
934 // Sm ----------------------------------------------
935 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
936 2, 4, 4, 6, 6, 1, 2, 2, 4, 1,
937
938 // Eu ----------------------------------------------
939 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
940 2, 4, 4, 6, 6, 1, 2, 2, 4, 2,
941
942 // Gd ----------------------------------------------
943 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
944 2, 4, 4, 6, 6, 1, 2, 2, 4, 1,
945 2,
946
947 // Tb ----------------------------------------------
948 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
949 2, 4, 4, 6, 6, 2, 2, 2, 4, 1,
950 2,
951
952 // Dy ----------------------------------------------
953 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
954 2, 4, 4, 6, 6, 4, 2, 2, 4, 2,
955
956 // Ho ----------------------------------------------
957 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
958 2, 4, 4, 6, 6, 5, 2, 2, 4, 2,
959
960 // Er ----------------------------------------------
961 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
962 2, 4, 4, 6, 6, 6, 2, 2, 4, 2,
963
964 // Tm ----------------------------------------------
965 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
966 2, 4, 4, 6, 6, 7, 2, 2, 4, 2,
967
968 // Yb ----------------------------------------------
969 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
970 2, 4, 4, 6, 6, 8, 2, 2, 4, 2,
971
972 // Lu ----------------------------------------------
973 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
974 2, 4, 4, 6, 6, 8, 2, 2, 4, 1,
975 2,
976
977 // Hf ----------------------------------------------
978 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
979 2, 4, 4, 6, 6, 8, 2, 2, 4, 2,
980 2,
981
982 // Ta ----------------------------------------------
983 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
984 2, 4, 4, 6, 6, 8, 2, 2, 4, 3,
985 2,
986
987 // W ----------------------------------------------
988 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
989 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
990 2,
991
992 // Re ----------------------------------------------
993 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
994 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
995 1, 2,
996
997 // Os ----------------------------------------------
998 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
999 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1000 2, 2,
1001
1002 // Ir ----------------------------------------------
1003 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1004 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1005 3, 2,
1006
1007 // Pt ----------------------------------------------
1008 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1009 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1010 5, 1,
1011
1012 // Au ----------------------------------------------
1013 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1014 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1015 6, 1,
1016
1017 // Hg ----------------------------------------------
1018 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1019 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1020 6, 2,
1021
1022 // Tl ----------------------------------------------
1023 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1024 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1025 6, 2, 1,
1026
1027 // Pb ----------------------------------------------
1028 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1029 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1030 6, 2, 2,
1031
1032 // Bi ----------------------------------------------
1033 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1034 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1035 6, 2, 2, 1,
1036
1037 // Po ----------------------------------------------
1038 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1039 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1040 6, 2, 2, 2,
1041
1042 // At ----------------------------------------------
1043 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1044 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1045 6, 2, 2, 3,
1046
1047 // Rn ----------------------------------------------
1048 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1049 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1050 6, 2, 2, 4,
1051
1052 // Fr ----------------------------------------------
1053 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1054 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1055 6, 2, 2, 4, 1,
1056
1057 // Ra ----------------------------------------------
1058 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1059 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1060 6, 2, 2, 4, 2,
1061
1062 // Ac ----------------------------------------------
1063 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1064 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1065 6, 2, 2, 4, 1, 2,
1066
1067 // Th ----------------------------------------------
1068 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1069 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1070 6, 2, 2, 4, 2, 2,
1071
1072 // Pa ----------------------------------------------
1073 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1074 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1075 6, 2, 2, 2, 4, 1, 2,
1076
1077 // U ----------------------------------------------
1078 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1079 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1080 6, 3, 2, 2, 4, 1, 2,
1081
1082 // Np ----------------------------------------------
1083 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1084 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1085 6, 4, 2, 2, 4, 1, 2,
1086
1087 // Pu ----------------------------------------------
1088 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1089 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1090 6, 6, 2, 2, 4, 2,
1091
1092 // Am ----------------------------------------------
1093 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1094 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1095 6, 6, 1, 2, 2, 4, 2,
1096
1097 // Cm ----------------------------------------------
1098 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1099 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1100 6, 6, 1, 2, 2, 4, 1, 2,
1101
1102 // Bk ----------------------------------------------
1103 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1104 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1105 6, 6, 2, 2, 2, 4, 1, 2,
1106
1107 // Cf ----------------------------------------------
1108 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1109 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1110 6, 6, 4, 2, 2, 4, 2,
1111
1112 // Es ----------------------------------------------
1113 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1114 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1115 6, 6, 5, 2, 2, 4, 2,
1116
1117 // Fm ----------------------------------------------
1118 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1119 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1120 6, 6, 6, 2, 2, 4, 2,
1121
1122 // Md ----------------------------------------------
1123 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1124 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1125 6, 6, 7, 2, 2, 4, 2,
1126
1127 // No ----------------------------------------------
1128 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1129 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1130 6, 6, 8, 2, 2, 4, 2,
1131
1132 // Lr ----------------------------------------------
1133 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1134 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1135 6, 6, 8, 2, 2, 4, 1, 2,
1136
1137 // Rf ----------------------------------------------
1138 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1139 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1140 6, 6, 8, 2, 2, 4, 2, 2,
1141
1142 // Db ----------------------------------------------
1143 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1144 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1145 6, 6, 8, 2, 2, 4, 3, 2,
1146
1147 // Sg ----------------------------------------------
1148 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1149 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1150 6, 6, 8, 2, 2, 4, 4, 2,
1151
1152 // Bh ----------------------------------------------
1153 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1154 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1155 6, 6, 8, 2, 2, 4, 4, 1, 2,
1156
1157 // Hs ----------------------------------------------
1158 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1159 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1160 6, 6, 8, 2, 2, 4, 4, 2, 2,
1161
1162 // Mt ----------------------------------------------
1163 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1164 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1165 6, 6, 8, 2, 2, 4, 4, 3, 2,
1166
1167 // Ds ----------------------------------------------
1168 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1169 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1170 6, 6, 8, 2, 2, 4, 4, 5, 1,
1171
1172 // Rg ----------------------------------------------
1173 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1174 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1175 6, 6, 8, 2, 2, 4, 4, 6, 1,
1176
1177 // Cp ----------------------------------------------
1178 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1179 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1180 6, 6, 8, 2, 2, 4, 4, 6, 2,
1181
1182 // Uut ---------------------------------------------
1183 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1184 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1185 6, 6, 8, 2, 2, 4, 4, 6, 2, 1,
1186
1187 // Fl ----------------------------------------------
1188 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1189 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1190 6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
1191
1192 // Uup ---------------------------------------------
1193 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1194 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1195 6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
1196 1,
1197
1198 // Lv ----------------------------------------------
1199 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1200 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1201 6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
1202 2,
1203
1204 // Uus ---------------------------------------------
1205 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1206 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1207 6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
1208 3,
1209
1210 // Uuo ---------------------------------------------
1211 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1212 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1213 6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
1214 4,
1215
1216 // Uue ---------------------------------------------
1217 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1218 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1219 6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
1220 4, 1,
1221
1222 // Ubn ---------------------------------------------
1223 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1224 2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1225 6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
1226 4, 2
1227};
1228
1229
1230G4int
1232{
1233#ifdef G4VERBOSE
1234 if (Z < 0 || Z > 120) { Z = PrintErrorZ(Z, "GetNumberOfShells"); }
1235#endif
1236 return fNumberOfShells[Z];
1237}
1238
1239
1242{
1243#ifdef G4VERBOSE
1244 if (Z < 0 || Z > 120) { Z = PrintErrorZ(Z, "GetBindingEnergy"); }
1245 if(ShellNb < 0 || ShellNb >= fNumberOfShells[Z]) {
1246 ShellNb = PrintErrorShell(Z, ShellNb, "GetBindingEnergy");
1247 }
1248#endif
1249 return fBindingEnergies[fIndexOfShells[Z] + ShellNb]*CLHEP::keV;
1250}
1251
1252
1253G4int
1255{
1256#ifdef G4VERBOSE
1257 if (Z < 0 || Z > 120) { Z = PrintErrorZ(Z, "GetNumberOfElectrons"); }
1258 if(ShellNb < 0 || ShellNb >= fNumberOfShells[Z]) {
1259 ShellNb = PrintErrorShell(Z, ShellNb, "GetNumberOfElectrons");
1260 }
1261#endif
1262 return fNumberOfElectrons[fIndexOfShells[Z] + ShellNb];
1263}
1264
1265//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
1266
1267G4int
1269{
1270#ifdef G4VERBOSE
1271 if(Z<0 || Z>120) { Z = PrintErrorZ(Z, "GetNumberOfFreeElectrons"); }
1272#endif
1273 G4int idx = fIndexOfShells[Z];
1274 G4int idxmax = idx + fNumberOfShells[Z];
1275 G4int n = 0;
1276 for (G4int i=idx; i<idxmax; ++i) {
1277 if(fBindingEnergies[i]*CLHEP::keV <= th) { n += fNumberOfElectrons[i]; }
1278 }
1279 return n;
1280}
1281
1283{
1284#ifdef G4VERBOSE
1285 if ( Z < 0 || Z > 120) { Z = PrintErrorZ(Z, "GetTotalBindingEnergy"); }
1286#endif
1287 G4int idx = fIndexOfShells[Z];
1288 G4int idxmax = idx + fNumberOfShells[Z];
1289 G4double energy = 0.0;
1290 for (G4int i=idx; i<idxmax; ++i) {energy += fBindingEnergies[i];}
1291 return energy*CLHEP::keV;
1292}
1293
1294
1295#ifdef G4VERBOSE
1296
1297G4int G4AtomicShells_XDB_EADL::PrintErrorZ(G4int Z, const G4String& ss)
1298{
1299 G4String sss = "G4AtomicShells::"+ss+"()";
1301 ed << "Atomic number out of range Z= " << Z;
1302 G4Exception(sss,"mat060",FatalException,ed,"");
1303 return 1;
1304}
1305
1306
1307G4int
1308G4AtomicShells_XDB_EADL::PrintErrorShell(G4int Z, G4int n, const G4String& ss)
1309{
1310 G4String sss = "G4AtomicShells::"+ss+"()";
1312 ed << "Shell number out of range Nshell= " << n << " Z= " << Z;
1313 G4Exception(sss,"mat061",FatalException,ed,"");
1314 return 0;
1315}
1316#endif
1317
static const char sss[MAX_N_PAR+2]
Definition: Evaluator.cc:63
@ FatalException
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition: G4Exception.cc:35
std::ostringstream G4ExceptionDescription
Definition: G4Exception.hh:40
double G4double
Definition: G4Types.hh:83
int G4int
Definition: G4Types.hh:85
const G4int Z[17]
static G4int GetNumberOfShells(G4int Z)
static const G4int fNumberOfElectrons[2171]
static G4int GetNumberOfElectrons(G4int Z, G4int SubshellNb)
static const G4double fBindingEnergies[2171]
static G4double GetBindingEnergy(G4int Z, G4int SubshellNb)
static G4double GetTotalBindingEnergy(G4int Z)
static const G4int fNumberOfShells[121]
static G4int GetNumberOfFreeElectrons(G4int Z, G4double th)
static const G4int fIndexOfShells[121]
static constexpr double keV
G4double energy(const ThreeVector &p, const G4double m)