1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
|
// This script was automatically generated!
// Edit dev/Source/Utils/write_atomic_data_cpp.py instead!
#ifndef WARPX_IONIZATION_TABLE_H_
#define WARPX_IONIZATION_TABLE_H_
#include <AMReX_AmrCore.H>
#include <map>
std::map<std::string, int> ion_map_ids = {
{"H", 0},
{"He", 1},
{"Li", 2},
{"Be", 3},
{"B", 4},
{"C", 5},
{"N", 6},
{"O", 7},
{"F", 8},
{"Ne", 9},
{"Na", 10},
{"Mg", 11},
{"Al", 12},
{"Si", 13},
{"P", 14},
{"S", 15},
{"Cl", 16},
{"Ar", 17},
{"Kr", 18},
{"Rb", 19},
{"Xe", 20},
{"Rn", 21} };
const int nelements = 22;
const int ion_atomic_numbers[nelements] = {
1, 2, 3, 4, 5, 6, 7, 8, 9, 10,
11, 12, 13, 14, 15, 16, 17, 18, 36, 37,
54, 86};
const int ion_energy_offsets[nelements] = {
0, 1, 3, 6, 10, 15, 21, 28, 36, 45,
55, 66, 78, 91, 105, 120, 136, 153, 171, 207,
244, 298};
const int energies_tab_length = 384;
const amrex::Real table_ionization_energies[energies_tab_length]{
// H
13.59843449,
// He
24.58738880, 54.4177650,
// Li
5.39171495, 75.6400964, 122.4543581,
// Be
9.322699, 18.21115, 153.896203, 217.7185843,
// B
8.298019, 25.15483, 37.93058, 259.3715, 340.226020,
// C
11.2602880, 24.383154, 47.88778, 64.49352, 392.090515, 489.993194,
// N
14.53413, 29.60125, 47.4453, 77.4735, 97.8901, 552.06732, 667.046116,
// O
13.618055, 35.12112, 54.93554, 77.41350, 113.8990, 138.1189, 739.32682,
871.40988,
// F
17.42282, 34.97081, 62.70798, 87.175, 114.249, 157.16311, 185.1868,
953.89804, 1103.11747,
// Ne
21.564540, 40.96297, 63.4233, 97.1900, 126.247, 157.934, 207.271,
239.0970, 1195.80783, 1362.19915,
// Na
5.1390769, 47.28636, 71.6200, 98.936, 138.404, 172.23, 208.504,
264.192, 299.856, 1465.13449, 1648.70218,
// Mg
7.646236, 15.035271, 80.1436, 109.2654, 141.33, 186.76, 225.02,
265.924, 327.99, 367.489, 1761.80487, 1962.66365,
// Al
5.985769, 18.82855, 28.447642, 119.9924, 153.8252, 190.49, 241.76,
284.64, 330.21, 398.65, 442.005, 2085.97700, 2304.14005,
// Si
8.15168, 16.34585, 33.49300, 45.14179, 166.767, 205.279, 246.57,
303.59, 351.28, 401.38, 476.273, 523.415, 2437.65813, 2673.17753,
// P
10.486686, 19.76949, 30.20264, 51.44387, 65.02511, 220.430, 263.57,
309.60, 372.31, 424.40, 479.44, 560.62, 611.741, 2816.90876,
3069.8415,
// S
10.36001, 23.33788, 34.86, 47.222, 72.5945, 88.0529, 280.954,
328.794, 379.84, 447.7, 504.55, 564.41, 651.96, 706.994,
3223.7807, 3494.1879,
// Cl
12.967632, 23.81364, 39.80, 53.24, 67.68, 96.94, 114.2013,
348.306, 400.851, 456.7, 530.0, 591.58, 656.30, 750.23,
809.198, 3658.3437, 3946.2909,
// Ar
15.7596117, 27.62967, 40.735, 59.58, 74.84, 91.290, 124.41,
143.4567, 422.60, 479.76, 540.4, 619.0, 685.5, 755.13,
855.5, 918.375, 4120.6656, 4426.2228,
// Kr
13.9996053, 24.35984, 35.838, 50.85, 64.69, 78.49, 109.13,
125.802, 233.0, 268, 308, 350, 391, 446,
492, 540, 591, 640, 785, 831.6, 882.8,
945, 999.0, 1042, 1155.0, 1205.23, 2928.9, 3072,
3228, 3380, 3584, 3752.0, 3971, 4109.083, 17296.420,
17936.209,
// Rb
4.1771280, 27.28954, 39.247, 52.20, 68.44, 82.9, 98.67,
132.79, 150.628, 277.12, 313.1, 356.0, 400, 443,
502, 550, 601, 654, 706.0, 857, 905.3,
958.9, 1024, 1080, 1125, 1242.5, 1294.57, 3133.3,
3281, 3443, 3600, 3815, 3988, 4214, 4356.865,
18305.884, 18965.516,
// Xe
12.1298436, 20.975, 31.05, 42.20, 54.1, 66.703, 91.6,
105.9778, 179.84, 202.0, 229.02, 255.0, 281, 314,
343, 374, 404, 434, 549, 582, 616,
650, 700, 736, 818, 857.0, 1493, 1571,
1653, 1742, 1826, 1919, 2023, 2113, 2209,
2300, 2556, 2637, 2726, 2811, 2975, 3068,
3243, 3333.8, 7660, 7889, 8144, 8382, 8971,
9243, 9581, 9810.37, 40271.724, 41299.71,
// Rn
10.74850, 21.4, 29.4, 36.9, 52.9, 64.0, 88.0,
102.0, 154.0, 173.9, 195.0, 218.0, 240, 264,
293, 317, 342, 367, 488, 520, 550,
580, 640, 680, 760, 800, 850, 920,
980, 1050, 1110, 1180, 1250, 1310, 1390,
1460, 1520, 1590, 1660, 1720, 2033, 2094,
2158, 2227, 2293, 2357, 2467, 2535, 2606,
2674, 2944, 3010, 3082, 3149, 3433, 3510,
3699, 3777, 6169, 6318, 6476, 6646, 6807,
6964, 7283, 7450, 7630, 7800, 8260, 8410,
8570, 8710, 9610, 9780, 10120, 10290, 21770,
22160, 22600, 22990, 26310, 26830, 27490, 27903.1,
110842.0, 112843.7
};
#endif // #ifndef WARPX_IONIZATION_TABLE_H_
|
