Skip to content

Commit

Permalink
add K,Ca,Sc,Zn,Ga,Ge,As,Se,Br to ionization tables (#3646)
Browse files Browse the repository at this point in the history
  • Loading branch information
@lucafedeli88
lucafedeli88 authored Jan 30, 2023
1 parent 8a9f9c4 commit 42661be
Show file tree
Hide file tree
Showing 2 changed files with 593 additions and 12 deletions.
197 changes: 185 additions & 12 deletions Source/Utils/Physics/IonizationEnergiesTable.H
View file
Original file line number Diff line number Diff line change
Expand Up @@ -30,26 +30,43 @@ namespace utils::physics
{"S", 15},
{"Cl", 16},
{"Ar", 17},
{"Ti", 18},
{"Cu", 19},
{"Kr", 20},
{"Rb", 21},
{"Xe", 22},
{"Rn", 23} };
{"K", 18},
{"Ca", 19},
{"Sc", 20},
{"Ti", 21},
{"V", 22},
{"Cr", 23},
{"Mn", 24},
{"Fe", 25},
{"Co", 26},
{"Ni", 27},
{"Cu", 28},
{"Zn", 29},
{"Ga", 30},
{"Ge", 31},
{"As", 32},
{"Se", 33},
{"Br", 34},
{"Kr", 35},
{"Rb", 36},
{"Xe", 37},
{"Rn", 38} };

constexpr int nelements = 24;
constexpr int nelements = 39;

constexpr int ion_atomic_numbers[nelements] = {
1, 2, 3, 4, 5, 6, 7, 8, 9, 10,
11, 12, 13, 14, 15, 16, 17, 18, 22, 29,
36, 37, 54, 86};
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, 54, 86};

constexpr int ion_energy_offsets[nelements] = {
0, 1, 3, 6, 10, 15, 21, 28, 36, 45,
55, 66, 78, 91, 105, 120, 136, 153, 171, 193,
222, 258, 295, 349};
55, 66, 78, 91, 105, 120, 136, 153, 171, 190,
210, 231, 253, 276, 300, 325, 351, 378, 406, 435,
465, 496, 528, 561, 595, 630, 666, 703, 757};

constexpr int energies_tab_length = 435;
constexpr int energies_tab_length = 843;

constexpr amrex::Real table_ionization_energies[energies_tab_length]{
// H
Expand Down Expand Up @@ -133,6 +150,30 @@ namespace utils::physics
amrex::Real(479.76), amrex::Real(540.4), amrex::Real(619.0),
amrex::Real(685.5), amrex::Real(755.13), amrex::Real(855.5),
amrex::Real(918.375), amrex::Real(4120.6656), amrex::Real(4426.2228),
// K
amrex::Real(4.34066373), amrex::Real(31.62500), amrex::Real(45.8031),
amrex::Real(60.917), amrex::Real(82.66), amrex::Real(99.44),
amrex::Real(117.56), amrex::Real(154.87), amrex::Real(175.8174),
amrex::Real(503.67), amrex::Real(565.6), amrex::Real(631.1),
amrex::Real(714.7), amrex::Real(786.3), amrex::Real(860.92),
amrex::Real(967.7), amrex::Real(1034.542), amrex::Real(4610.87018),
amrex::Real(4934.04979),
// Ca
amrex::Real(6.11315547), amrex::Real(11.871719), amrex::Real(50.91316),
amrex::Real(67.2732), amrex::Real(84.34), amrex::Real(108.78),
amrex::Real(127.21), amrex::Real(147.24), amrex::Real(188.54),
amrex::Real(211.275), amrex::Real(591.60), amrex::Real(658.2),
amrex::Real(728.6), amrex::Real(817.2), amrex::Real(894.0),
amrex::Real(973.7), amrex::Real(1086.8), amrex::Real(1157.726),
amrex::Real(5128.8578), amrex::Real(5469.86358),
// Sc
amrex::Real(6.56149), amrex::Real(12.79977), amrex::Real(24.756839),
amrex::Real(73.4894), amrex::Real(91.95), amrex::Real(110.68),
amrex::Real(137.99), amrex::Real(158.08), amrex::Real(180.03),
amrex::Real(225.18), amrex::Real(249.798), amrex::Real(687.36),
amrex::Real(757.7), amrex::Real(833.2), amrex::Real(926.5),
amrex::Real(1008.6), amrex::Real(1093.5), amrex::Real(1213.1),
amrex::Real(1287.957), amrex::Real(5674.9037), amrex::Real(6033.75643),
// Ti
amrex::Real(6.828120), amrex::Real(13.5755), amrex::Real(27.49171),
amrex::Real(43.26717), amrex::Real(99.299), amrex::Real(119.533),
Expand All @@ -142,6 +183,65 @@ namespace utils::physics
amrex::Real(1042.5), amrex::Real(1130.2), amrex::Real(1220.3),
amrex::Real(1346.3), amrex::Real(1425.257), amrex::Real(6249.0226),
amrex::Real(6625.81023),
// V
amrex::Real(6.746187), amrex::Real(14.634), amrex::Real(29.3111),
amrex::Real(46.709), amrex::Real(65.28165), amrex::Real(128.125),
amrex::Real(150.72), amrex::Real(173.55), amrex::Real(206.0),
amrex::Real(230.5), amrex::Real(254.8), amrex::Real(308.5),
amrex::Real(336.274), amrex::Real(896.0), amrex::Real(977.2),
amrex::Real(1062.9), amrex::Real(1165.2), amrex::Real(1258.9),
amrex::Real(1354.2), amrex::Real(1486.7), amrex::Real(1569.656),
amrex::Real(6851.3109), amrex::Real(7246.12624),
// Cr
amrex::Real(6.76651), amrex::Real(16.486305), amrex::Real(30.959),
amrex::Real(49.16), amrex::Real(69.46), amrex::Real(90.6349),
amrex::Real(160.29), amrex::Real(184.76), amrex::Real(209.5),
amrex::Real(244.5), amrex::Real(270.8), amrex::Real(296.7),
amrex::Real(354.7), amrex::Real(384.163), amrex::Real(1011.6),
amrex::Real(1097.2), amrex::Real(1188.0), amrex::Real(1294.8),
amrex::Real(1394.5), amrex::Real(1495.1), amrex::Real(1634.1),
amrex::Real(1721.183), amrex::Real(7481.8628), amrex::Real(7894.80289),
// Mn
amrex::Real(7.4340380), amrex::Real(15.63999), amrex::Real(33.668),
amrex::Real(51.21), amrex::Real(72.41), amrex::Real(95.604),
amrex::Real(119.203), amrex::Real(195.5), amrex::Real(221.89),
amrex::Real(248.6), amrex::Real(286.1), amrex::Real(314.4),
amrex::Real(343.6), amrex::Real(402.95), amrex::Real(435.172),
amrex::Real(1133.7), amrex::Real(1224.1), amrex::Real(1320.3),
amrex::Real(1430.9), amrex::Real(1537.2), amrex::Real(1643.2),
amrex::Real(1788.7), amrex::Real(1879.873), amrex::Real(8140.7872),
amrex::Real(8571.95438),
// Fe
amrex::Real(7.9024681), amrex::Real(16.19921), amrex::Real(30.651),
amrex::Real(54.91), amrex::Real(75.00), amrex::Real(98.985),
amrex::Real(124.9671), amrex::Real(151.060), amrex::Real(233.6),
amrex::Real(262.10), amrex::Real(290.9), amrex::Real(330.8),
amrex::Real(361.0), amrex::Real(392.2), amrex::Real(456.2),
amrex::Real(489.312), amrex::Real(1262.7), amrex::Real(1357.8),
amrex::Real(1460), amrex::Real(1575.6), amrex::Real(1687.0),
amrex::Real(1798.4), amrex::Real(1950.4), amrex::Real(2045.759),
amrex::Real(8828.1879), amrex::Real(9277.6886),
// Co
amrex::Real(7.88101), amrex::Real(17.0844), amrex::Real(33.50),
amrex::Real(51.27), amrex::Real(79.50), amrex::Real(102.00),
amrex::Real(128.9), amrex::Real(157.8), amrex::Real(186.14),
amrex::Real(275.4), amrex::Real(305.32), amrex::Real(336.1),
amrex::Real(378.5), amrex::Real(410.0), amrex::Real(441.1),
amrex::Real(511.96), amrex::Real(546.588), amrex::Real(1397.2),
amrex::Real(1504.5), amrex::Real(1606), amrex::Real(1724),
amrex::Real(1844.0), amrex::Real(1960.8), amrex::Real(2119.4),
amrex::Real(2218.876), amrex::Real(9544.1833), amrex::Real(10012.1297),
// Ni
amrex::Real(7.639878), amrex::Real(18.168838), amrex::Real(35.187),
amrex::Real(54.92), amrex::Real(76.06), amrex::Real(108.0),
amrex::Real(132.0), amrex::Real(162.0), amrex::Real(193.2),
amrex::Real(224.7), amrex::Real(319.5), amrex::Real(351.6),
amrex::Real(384.5), amrex::Real(429.3), amrex::Real(462.8),
amrex::Real(495.4), amrex::Real(571.07), amrex::Real(607.020),
amrex::Real(1540.1), amrex::Real(1646), amrex::Real(1758),
amrex::Real(1880), amrex::Real(2008.1), amrex::Real(2130.5),
amrex::Real(2295.6), amrex::Real(2399.259), amrex::Real(10288.8862),
amrex::Real(10775.3948),
// Cu
amrex::Real(7.726380), amrex::Real(20.29239), amrex::Real(36.841),
amrex::Real(57.38), amrex::Real(79.8), amrex::Real(103.0),
Expand All @@ -153,6 +253,79 @@ namespace utils::physics
amrex::Real(1918), amrex::Real(2044), amrex::Real(2179.4),
amrex::Real(2307.3), amrex::Real(2479.1), amrex::Real(2586.954),
amrex::Real(11062.4313), amrex::Real(11567.613),
// Zn
amrex::Real(9.394197), amrex::Real(17.96439), amrex::Real(39.72330),
amrex::Real(59.573), amrex::Real(82.6), amrex::Real(108.0),
amrex::Real(133.9), amrex::Real(173.9), amrex::Real(203.0),
amrex::Real(238.0), amrex::Real(274.4), amrex::Real(310.8),
amrex::Real(417.6), amrex::Real(453.4), amrex::Real(490.6),
amrex::Real(540.0), amrex::Real(577.8), amrex::Real(613.3),
amrex::Real(697.5), amrex::Real(737.366), amrex::Real(1846.8),
amrex::Real(1961), amrex::Real(2085), amrex::Real(2214),
amrex::Real(2358.0), amrex::Real(2491.5), amrex::Real(2669.9),
amrex::Real(2781.996), amrex::Real(11864.9399), amrex::Real(12388.9427),
// Ga
amrex::Real(5.9993020), amrex::Real(20.51514), amrex::Real(30.72576),
amrex::Real(63.241), amrex::Real(86.01), amrex::Real(112.7),
amrex::Real(140.8), amrex::Real(169.9), amrex::Real(211.0),
amrex::Real(244.0), amrex::Real(280.0), amrex::Real(319.0),
amrex::Real(356), amrex::Real(471.2), amrex::Real(508.8),
amrex::Real(548.3), amrex::Real(599.8), amrex::Real(640.0),
amrex::Real(677), amrex::Real(765.7), amrex::Real(807.308),
amrex::Real(2010.0), amrex::Real(2129), amrex::Real(2258),
amrex::Real(2391), amrex::Real(2543.9), amrex::Real(2683.0),
amrex::Real(2868), amrex::Real(2984.426), amrex::Real(12696.5575),
amrex::Real(13239.5029),
// Ge
amrex::Real(7.899435), amrex::Real(15.934610), amrex::Real(34.0576),
amrex::Real(45.7155), amrex::Real(90.500), amrex::Real(115.90),
amrex::Real(144.9), amrex::Real(176.4), amrex::Real(212.5),
amrex::Real(252.1), amrex::Real(286), amrex::Real(326.0),
amrex::Real(367), amrex::Real(407), amrex::Real(527.9),
amrex::Real(567.3), amrex::Real(609.1), amrex::Real(662.8),
amrex::Real(706.7), amrex::Real(744), amrex::Real(837.1),
amrex::Real(880.44), amrex::Real(2178.2), amrex::Real(2304),
amrex::Real(2439), amrex::Real(2575), amrex::Real(2737.1),
amrex::Real(2881.9), amrex::Real(3074), amrex::Real(3194.293),
amrex::Real(13557.4208), amrex::Real(14119.4457),
// As
amrex::Real(9.78855), amrex::Real(18.5892), amrex::Real(28.349),
amrex::Real(50.15), amrex::Real(62.77), amrex::Real(121.19),
amrex::Real(147.0), amrex::Real(180.0), amrex::Real(213.0),
amrex::Real(247.0), amrex::Real(296.0), amrex::Real(333.0),
amrex::Real(375), amrex::Real(418), amrex::Real(460),
amrex::Real(587.6), amrex::Real(628.8), amrex::Real(672.9),
amrex::Real(728.9), amrex::Real(774.0), amrex::Real(814),
amrex::Real(911.7), amrex::Real(956.79), amrex::Real(2356.9),
amrex::Real(2486), amrex::Real(2626), amrex::Real(2766),
amrex::Real(2938), amrex::Real(3088.1), amrex::Real(3287),
amrex::Real(3411.643), amrex::Real(14447.678), amrex::Real(15028.9251),
// Se
amrex::Real(9.752368), amrex::Real(21.196), amrex::Real(31.697),
amrex::Real(42.947), amrex::Real(68.30), amrex::Real(81.83),
amrex::Real(155.327), amrex::Real(184.0), amrex::Real(219.0),
amrex::Real(255.0), amrex::Real(291.0), amrex::Real(342.9),
amrex::Real(383), amrex::Real(426), amrex::Real(473),
amrex::Real(517), amrex::Real(650.5), amrex::Real(693.4),
amrex::Real(739.8), amrex::Real(798), amrex::Real(845.8),
amrex::Real(887), amrex::Real(989.6), amrex::Real(1036.36),
amrex::Real(2540.7), amrex::Real(2674), amrex::Real(2820),
amrex::Real(2964), amrex::Real(3146), amrex::Real(3301.8),
amrex::Real(3507), amrex::Real(3636.526), amrex::Real(15367.491),
amrex::Real(15968.1075),
// Br
amrex::Real(11.81381), amrex::Real(21.591), amrex::Real(34.871),
amrex::Real(47.782), amrex::Real(59.595), amrex::Real(87.390),
amrex::Real(103.03), amrex::Real(192.61), amrex::Real(224.0),
amrex::Real(261.0), amrex::Real(301.0), amrex::Real(338.0),
amrex::Real(393), amrex::Real(436), amrex::Real(481),
amrex::Real(530), amrex::Real(577), amrex::Real(716.3),
amrex::Real(761.0), amrex::Real(809.8), amrex::Real(870),
amrex::Real(920.8), amrex::Real(963), amrex::Real(1070.6),
amrex::Real(1119.17), amrex::Real(2731.4), amrex::Real(2869),
amrex::Real(3021), amrex::Real(3169), amrex::Real(3361),
amrex::Real(3523.1), amrex::Real(3735), amrex::Real(3868.986),
amrex::Real(16317.011), amrex::Real(16937.1497),
// Kr
amrex::Real(13.9996053), amrex::Real(24.35984), amrex::Real(35.838),
amrex::Real(50.85), amrex::Real(64.69), amrex::Real(78.49),
Expand Down
Loading

0 comments on commit 42661be

Please sign in to comment.