Hipace
IonizationEnergiesTable.H
Go to the documentation of this file.
1 /* Copyright 2021
2  *
3  * This file is part of HiPACE++.
4  *
5  * Authors: AlexanderSinn
6  * License: BSD-3-Clause-LBNL
7  */
8 // This script was automatically generated!
9 // Edit src/utils/write_atomic_data_cpp.py instead!
10 #ifndef HIPACE_IONIZATION_TABLE_H_
11 #define HIPACE_IONIZATION_TABLE_H_
12 
13 #include <AMReX_AmrCore.H>
14 #include <AMReX_REAL.H>
15 #include <map>
16 
17 // Reference:
18 // Kramida, A., Ralchenko, Yu., Reader, J., and NIST ASD Team (2014).
19 // NIST Atomic Spectra Database (ver. 5.2), [Online].
20 // Available: http://physics.nist.gov/asd [2017, March 3].
21 //
22 // The Data written below is a reformatting of the data referenced form NIST.
23 
24 std::map<std::string, int> ion_map_ids = {
25  {"H", 0},
26  {"He", 1},
27  {"Li", 2},
28  {"Be", 3},
29  {"B", 4},
30  {"C", 5},
31  {"N", 6},
32  {"O", 7},
33  {"F", 8},
34  {"Ne", 9},
35  {"Na", 10},
36  {"Mg", 11},
37  {"Al", 12},
38  {"Si", 13},
39  {"P", 14},
40  {"S", 15},
41  {"Cl", 16},
42  {"Ar", 17},
43  {"Cu", 18},
44  {"Kr", 19},
45  {"Rb", 20},
46  {"Xe", 21},
47  {"Rn", 22} };
48 
49 constexpr int nelements = 23;
50 
51 constexpr int ion_atomic_numbers[nelements] = {
52  1, 2, 3, 4, 5, 6, 7, 8, 9, 10,
53  11, 12, 13, 14, 15, 16, 17, 18, 29, 36,
54  37, 54, 86};
55 
56 constexpr int ion_energy_offsets[nelements] = {
57  0, 1, 3, 6, 10, 15, 21, 28, 36, 45,
58  55, 66, 78, 91, 105, 120, 136, 153, 171, 200,
59  236, 273, 327};
60 
61 constexpr int energies_tab_length = 413;
62 
64  // H
65  amrex::Real(13.59843449),
66  // He
67  amrex::Real(24.58738880), amrex::Real(54.4177650),
68  // Li
69  amrex::Real(5.39171495), amrex::Real(75.6400964), amrex::Real(122.4543581),
70  // Be
71  amrex::Real(9.322699), amrex::Real(18.21115), amrex::Real(153.896203),
72  amrex::Real(217.7185843),
73  // B
74  amrex::Real(8.298019), amrex::Real(25.15483), amrex::Real(37.93058),
75  amrex::Real(259.3715), amrex::Real(340.226020),
76  // C
77  amrex::Real(11.2602880), amrex::Real(24.383154), amrex::Real(47.88778),
78  amrex::Real(64.49352), amrex::Real(392.090515), amrex::Real(489.993194),
79  // N
80  amrex::Real(14.53413), amrex::Real(29.60125), amrex::Real(47.4453),
81  amrex::Real(77.4735), amrex::Real(97.8901), amrex::Real(552.06732),
82  amrex::Real(667.046116),
83  // O
84  amrex::Real(13.618055), amrex::Real(35.12112), amrex::Real(54.93554),
85  amrex::Real(77.41350), amrex::Real(113.8990), amrex::Real(138.1189),
86  amrex::Real(739.32682), amrex::Real(871.40988),
87  // F
88  amrex::Real(17.42282), amrex::Real(34.97081), amrex::Real(62.70798),
89  amrex::Real(87.175), amrex::Real(114.249), amrex::Real(157.16311),
90  amrex::Real(185.1868), amrex::Real(953.89804), amrex::Real(1103.11747),
91  // Ne
92  amrex::Real(21.564540), amrex::Real(40.96297), amrex::Real(63.4233),
93  amrex::Real(97.1900), amrex::Real(126.247), amrex::Real(157.934),
94  amrex::Real(207.271), amrex::Real(239.0970), amrex::Real(1195.80783),
95  amrex::Real(1362.19915),
96  // Na
97  amrex::Real(5.1390769), amrex::Real(47.28636), amrex::Real(71.6200),
98  amrex::Real(98.936), amrex::Real(138.404), amrex::Real(172.23),
99  amrex::Real(208.504), amrex::Real(264.192), amrex::Real(299.856),
100  amrex::Real(1465.13449), amrex::Real(1648.70218),
101  // Mg
102  amrex::Real(7.646236), amrex::Real(15.035271), amrex::Real(80.1436),
103  amrex::Real(109.2654), amrex::Real(141.33), amrex::Real(186.76),
104  amrex::Real(225.02), amrex::Real(265.924), amrex::Real(327.99),
105  amrex::Real(367.489), amrex::Real(1761.80487), amrex::Real(1962.66365),
106  // Al
107  amrex::Real(5.985769), amrex::Real(18.82855), amrex::Real(28.447642),
108  amrex::Real(119.9924), amrex::Real(153.8252), amrex::Real(190.49),
109  amrex::Real(241.76), amrex::Real(284.64), amrex::Real(330.21),
110  amrex::Real(398.65), amrex::Real(442.005), amrex::Real(2085.97700),
111  amrex::Real(2304.14005),
112  // Si
113  amrex::Real(8.15168), amrex::Real(16.34585), amrex::Real(33.49300),
114  amrex::Real(45.14179), amrex::Real(166.767), amrex::Real(205.279),
115  amrex::Real(246.57), amrex::Real(303.59), amrex::Real(351.28),
116  amrex::Real(401.38), amrex::Real(476.273), amrex::Real(523.415),
117  amrex::Real(2437.65813), amrex::Real(2673.17753),
118  // P
119  amrex::Real(10.486686), amrex::Real(19.76949), amrex::Real(30.20264),
120  amrex::Real(51.44387), amrex::Real(65.02511), amrex::Real(220.430),
121  amrex::Real(263.57), amrex::Real(309.60), amrex::Real(372.31),
122  amrex::Real(424.40), amrex::Real(479.44), amrex::Real(560.62),
123  amrex::Real(611.741), amrex::Real(2816.90876), amrex::Real(3069.8415),
124  // S
125  amrex::Real(10.36001), amrex::Real(23.33788), amrex::Real(34.86),
126  amrex::Real(47.222), amrex::Real(72.5945), amrex::Real(88.0529),
127  amrex::Real(280.954), amrex::Real(328.794), amrex::Real(379.84),
128  amrex::Real(447.7), amrex::Real(504.55), amrex::Real(564.41),
129  amrex::Real(651.96), amrex::Real(706.994), amrex::Real(3223.7807),
130  amrex::Real(3494.1879),
131  // Cl
132  amrex::Real(12.967632), amrex::Real(23.81364), amrex::Real(39.80),
133  amrex::Real(53.24), amrex::Real(67.68), amrex::Real(96.94),
134  amrex::Real(114.2013), amrex::Real(348.306), amrex::Real(400.851),
135  amrex::Real(456.7), amrex::Real(530.0), amrex::Real(591.58),
136  amrex::Real(656.30), amrex::Real(750.23), amrex::Real(809.198),
137  amrex::Real(3658.3437), amrex::Real(3946.2909),
138  // Ar
139  amrex::Real(15.7596117), amrex::Real(27.62967), amrex::Real(40.735),
140  amrex::Real(59.58), amrex::Real(74.84), amrex::Real(91.290),
141  amrex::Real(124.41), amrex::Real(143.4567), amrex::Real(422.60),
142  amrex::Real(479.76), amrex::Real(540.4), amrex::Real(619.0),
143  amrex::Real(685.5), amrex::Real(755.13), amrex::Real(855.5),
144  amrex::Real(918.375), amrex::Real(4120.6656), amrex::Real(4426.2228),
145  // Cu
146  amrex::Real(7.726380), amrex::Real(20.29239), amrex::Real(36.841),
147  amrex::Real(57.38), amrex::Real(79.8), amrex::Real(103.0),
148  amrex::Real(139.0), amrex::Real(166.0), amrex::Real(198.0),
149  amrex::Real(232.2), amrex::Real(265.33), amrex::Real(367.0),
150  amrex::Real(401.0), amrex::Real(436.0), amrex::Real(483.1),
151  amrex::Real(518.7), amrex::Real(552.8), amrex::Real(632.5),
152  amrex::Real(670.608), amrex::Real(1690.5), amrex::Real(1800),
153  amrex::Real(1918), amrex::Real(2044), amrex::Real(2179.4),
154  amrex::Real(2307.3), amrex::Real(2479.1), amrex::Real(2586.954),
155  amrex::Real(11062.4313), amrex::Real(11567.613),
156  // Kr
157  amrex::Real(13.9996053), amrex::Real(24.35984), amrex::Real(35.838),
158  amrex::Real(50.85), amrex::Real(64.69), amrex::Real(78.49),
159  amrex::Real(109.13), amrex::Real(125.802), amrex::Real(233.0),
160  amrex::Real(268), amrex::Real(308), amrex::Real(350),
161  amrex::Real(391), amrex::Real(446), amrex::Real(492),
162  amrex::Real(540), amrex::Real(591), amrex::Real(640),
163  amrex::Real(785), amrex::Real(831.6), amrex::Real(882.8),
164  amrex::Real(945), amrex::Real(999.0), amrex::Real(1042),
165  amrex::Real(1155.0), amrex::Real(1205.23), amrex::Real(2928.9),
166  amrex::Real(3072), amrex::Real(3228), amrex::Real(3380),
167  amrex::Real(3584), amrex::Real(3752.0), amrex::Real(3971),
168  amrex::Real(4109.083), amrex::Real(17296.420), amrex::Real(17936.209),
169  // Rb
170  amrex::Real(4.1771280), amrex::Real(27.28954), amrex::Real(39.247),
171  amrex::Real(52.20), amrex::Real(68.44), amrex::Real(82.9),
172  amrex::Real(98.67), amrex::Real(132.79), amrex::Real(150.628),
173  amrex::Real(277.12), amrex::Real(313.1), amrex::Real(356.0),
174  amrex::Real(400), amrex::Real(443), amrex::Real(502),
175  amrex::Real(550), amrex::Real(601), amrex::Real(654),
176  amrex::Real(706.0), amrex::Real(857), amrex::Real(905.3),
177  amrex::Real(958.9), amrex::Real(1024), amrex::Real(1080),
178  amrex::Real(1125), amrex::Real(1242.5), amrex::Real(1294.57),
179  amrex::Real(3133.3), amrex::Real(3281), amrex::Real(3443),
180  amrex::Real(3600), amrex::Real(3815), amrex::Real(3988),
181  amrex::Real(4214), amrex::Real(4356.865), amrex::Real(18305.884),
182  amrex::Real(18965.516),
183  // Xe
184  amrex::Real(12.1298436), amrex::Real(20.975), amrex::Real(31.05),
185  amrex::Real(42.20), amrex::Real(54.1), amrex::Real(66.703),
186  amrex::Real(91.6), amrex::Real(105.9778), amrex::Real(179.84),
187  amrex::Real(202.0), amrex::Real(229.02), amrex::Real(255.0),
188  amrex::Real(281), amrex::Real(314), amrex::Real(343),
189  amrex::Real(374), amrex::Real(404), amrex::Real(434),
190  amrex::Real(549), amrex::Real(582), amrex::Real(616),
191  amrex::Real(650), amrex::Real(700), amrex::Real(736),
192  amrex::Real(818), amrex::Real(857.0), amrex::Real(1493),
193  amrex::Real(1571), amrex::Real(1653), amrex::Real(1742),
194  amrex::Real(1826), amrex::Real(1919), amrex::Real(2023),
195  amrex::Real(2113), amrex::Real(2209), amrex::Real(2300),
196  amrex::Real(2556), amrex::Real(2637), amrex::Real(2726),
197  amrex::Real(2811), amrex::Real(2975), amrex::Real(3068),
198  amrex::Real(3243), amrex::Real(3333.8), amrex::Real(7660),
199  amrex::Real(7889), amrex::Real(8144), amrex::Real(8382),
200  amrex::Real(8971), amrex::Real(9243), amrex::Real(9581),
201  amrex::Real(9810.37), amrex::Real(40271.724), amrex::Real(41299.71),
202  // Rn
203  amrex::Real(10.74850), amrex::Real(21.4), amrex::Real(29.4),
204  amrex::Real(36.9), amrex::Real(52.9), amrex::Real(64.0),
205  amrex::Real(88.0), amrex::Real(102.0), amrex::Real(154.0),
206  amrex::Real(173.9), amrex::Real(195.0), amrex::Real(218.0),
207  amrex::Real(240), amrex::Real(264), amrex::Real(293),
208  amrex::Real(317), amrex::Real(342), amrex::Real(367),
209  amrex::Real(488), amrex::Real(520), amrex::Real(550),
210  amrex::Real(580), amrex::Real(640), amrex::Real(680),
211  amrex::Real(760), amrex::Real(800), amrex::Real(850),
212  amrex::Real(920), amrex::Real(980), amrex::Real(1050),
213  amrex::Real(1110), amrex::Real(1180), amrex::Real(1250),
214  amrex::Real(1310), amrex::Real(1390), amrex::Real(1460),
215  amrex::Real(1520), amrex::Real(1590), amrex::Real(1660),
216  amrex::Real(1720), amrex::Real(2033), amrex::Real(2094),
217  amrex::Real(2158), amrex::Real(2227), amrex::Real(2293),
218  amrex::Real(2357), amrex::Real(2467), amrex::Real(2535),
219  amrex::Real(2606), amrex::Real(2674), amrex::Real(2944),
220  amrex::Real(3010), amrex::Real(3082), amrex::Real(3149),
221  amrex::Real(3433), amrex::Real(3510), amrex::Real(3699),
222  amrex::Real(3777), amrex::Real(6169), amrex::Real(6318),
223  amrex::Real(6476), amrex::Real(6646), amrex::Real(6807),
224  amrex::Real(6964), amrex::Real(7283), amrex::Real(7450),
225  amrex::Real(7630), amrex::Real(7800), amrex::Real(8260),
226  amrex::Real(8410), amrex::Real(8570), amrex::Real(8710),
227  amrex::Real(9610), amrex::Real(9780), amrex::Real(10120),
228  amrex::Real(10290), amrex::Real(21770), amrex::Real(22160),
229  amrex::Real(22600), amrex::Real(22990), amrex::Real(26310),
230  amrex::Real(26830), amrex::Real(27490), amrex::Real(27903.1),
231  amrex::Real(110842.0), amrex::Real(112843.7)
232 };
233 
234 #endif // #ifndef HIPACE_IONIZATION_TABLE_H_
std::map< std::string, int > ion_map_ids
Definition: IonizationEnergiesTable.H:24
constexpr int ion_atomic_numbers[nelements]
Definition: IonizationEnergiesTable.H:51
constexpr int ion_energy_offsets[nelements]
Definition: IonizationEnergiesTable.H:56
constexpr int energies_tab_length
Definition: IonizationEnergiesTable.H:61
constexpr int nelements
Definition: IonizationEnergiesTable.H:49
constexpr amrex::Real table_ionization_energies[energies_tab_length]
Definition: IonizationEnergiesTable.H:63