CC BY
26
Ionization Efficiency in a Hot Flat Disc-Shaped Cavity
M. Turek
Institute of Physics, Maria Curie-Sklodowska University in Lublin, pl. M. Curie-Sklodowskiej 1, 20-031 Lublin, Poland
Received 14.02.2020
Accepted for publication 20.04.2020
Abstract
Hot cavity ion sources of different kinds are widely used in nuclear and mass spectroscopy, especially in on-line isotope separation devices attracting attention of scientists and engineers looking for high ionization efficiency, robustness and beam purity. In the paper a new type of hot ionizer cavity is proposed: namely cavity having the shape of a flat disc, which may be especially suitable for short-lived nuclides to be ionized.
A numerical model of the ion source is presented in the paper. The particle tracking code takes into account ionization at hot surfaces and enables modeling of both flat disc cavity and standard elongated cavity ionizers. The code enables calculation of total ionization efficiency and is suitable for stable and long-lived nuclides.
Influence of the flat disc cavity geometry (thickness and radius) and its temperature on total ionization efficiency was considered - it was shown that the efficiency increases with cavity radius due to the growing number of particle-wall collisions. This effect may be important in the case of the hard-to-ionize nuclides.
The optimal ionizer geometry is characterized by 90 % efficiency, even for substances with rather low ionization coefficient (of order 0.05). The role played by the size of the extraction opening is explained -it is demonstrated that the ionization efficiency increases due to the opening radius reduction. It is also proven that extraction voltage of 1-2 kV is sufficient to maintain optimal ionizer efficiency.
Keywords: ionizer, ionization efficiency, cavity, extraction opening.
DOI: 10.21122/2220-9506-2020-11-2-132-139
Адрес для переписки:
M. Turek
Institute of Physics, Maria Curie-Sklodowska University in Lublin, pl. M. Curie-Sklodowskiej 1, 20-031 Lublin, Poland e-mail: [email protected]
Address for correspondence:
M. Turek
Institute of Physics, Maria Curie-Sklodowska University in Lublin, pl. M. Curie-Sklodowskiej 1, 20-031 Lublin, Poland e-mail: [email protected]
Для цитирования:
M. Turek.
Ionization Efficiency in a Hot Flat Disc-Shaped Cavity. Приборы и методы измерений. 2020. - Т. 11, № 2. - С. 132-139. DOI: 10.21122/2220-9506-2020-11-2-132-139
For citation: M. Turek.
Ionization Efficiency in a Hot Flat Disc-Shaped Cavity. Devices and Methods of Measurements. 2020, vol. 11, no. 2, pp. 132-139. DOI: 10.21122/2220-9506-2020-11-2-132-139
Эффективность ионизации в горячей плоской
полости
М. Турек
Институт физики, Университет имени Марии Кюри-Склодовской в Люблине, пл. имени М. Кюри-Складовской, 1, г. Люблин 20-031, Польша
Поступила 14.02.2020 Принята к печати 20.04.2020
Ионизаторы с горячей полостью различного типа находят широкое применение в атомной и масс-спектроскопии, в частности, в устройствах для разделения изотопов в режиме онлайн, представляют большой интерес для учёных и инженеров вследствие высокой эффективности ионизации, надёжности и чистоты луча. В работе предложен новый тип горячей ионизационной полости, а именно полости в форме плоского диска, особенно эффективной для ионизации короткоживущих нуклидов.
Представлена численная модель ионного источника. Модель отслеживания частиц учитывает ионизацию на горячих поверхностях и позволяет моделировать как полость плоского диска, так и стандартные ионизаторы с удлинённой полостью. Модель позволяет выполнять расчёт общей эффективности ионизации и применима к стабильным и долгоживущим нуклидам.
Рассмотрено влияние геометрии полости плоского диска (толщина и радиус) и его температуры на общую эффективность ионизации - показано, что эффективность увеличивается с радиусом полости из-за растущего числа столкновений частиц со стенками. Данный эффект может оказаться важен для трудноионизируемых нуклидов.
Оптимальная геометрия ионизатора характеризуется эффективностью 90 % даже для трудноионизируемых веществ с коэффициентом ионизации порядка 0,05. Объясняется роль, которую играет размер экстракционного отверстия - показано, что эффективность ионизации увеличивается из-за уменьшения радиуса отверстия. Также доказано, что выходное напряжение 1-2 кВ достаточно для поддержания оптимальной эффективности.
Ключевые слова: ионизатор, эффективность ионизации, полость, экстракционное отверстие.
DOI: 10.21122/2220-9506-2020-11-2-132-139
Адрес для переписки:
M. Turek
Institute of Physics, Maria Curie-Sklodowska University in Lublin, pl. M. Curie-Sklodowskiej 1, 20-031 Lublin, Poland e-mail: [email protected]
Address for correspondence:
M. Turek
Institute of Physics, Maria Curie-Sklodowska University in Lublin, pl. M. Curie-Sklodowskiej 1, 20-031 Lublin, Poland e-mail: [email protected]
Для цитирования:
M. Turek.
Ionization Efficiency in a Hot Flat Disc-Shaped Cavity. Приборы и методы измерений. 2020. - Т. 11, № 2. - С. 132-139. DOI: 10.21122/2220-9506-2020-11-2-132-139
For citation: M. Turek.
Ionization Efficiency in a Hot Flat Disc-Shaped Cavity. Devices and Methods of Measurements. 2020, vol. 11, no. 2, pp. 132-139. DOI: 10.21122/2220-9506-2020-11-2-132-139
Introduction
Surface ionization ion sources equipped with hot cavities proved to be useful tools e.g. in nuclear spectroscopy [1], mass spectroscopy [2, 3] and are widely used in on-line isotope separation (ISOL) projects [4-7]. Hot cavity ion sources in its pure form were invented almost five decades ago [8, 9] and evolved into the resonant ionization laser ion sources (RILIS) [10-12] based on an excitation of the valence electron to the continuum by tunable lasers. It is worth mentioning that the evolution has not stopped and next hybrid ion sources for ISOL purposes still appear [13]. Numerous advantages of hot cavity ion sources like: compact design, robustness and reliability, very low energy spread of the ion beam and its high purity and, probably most of all, high efficiency and very small amounts of ionized substance needed to obtain a good quality ion beam made these devices attractive to physicist and engineers.
There is a variety of hot cavity ion sources, but in their classical form the basic part is a semi-opened ionizer made of refractory metal or ceramics, heated to a high working temperature (» 2500 K or more) either by electron beams [14] or ohmically [15, 16] or even inductively [17]. Most often the ionizer is connected to a target irradiated in order to produce isotopes by a kind of transfer line [18], but in some solutions the ionizer itself could serve as the target and the produced isotopes are released immediately into the hot cavity [14]. As it has been already said, the cylindrical ionizer is the most common solution, but spherical ones are also used in some devices [19, 20].
In a series of previous papers a numerical models of a variety of hot cavity ion source were presented, including those concerning tubular [21-23], spherical or hemispherical [24] and also conical ionizers [25]. Results presented by other groups also showed that ionizer of complex shapes are characterized by very high efficiencies [26, 27]. At the early version of the model ionization of stable isotopes was considered, later on the model was upgraded in order to take into account effects of radioactive decay and delays of the particle emission due to diffusion and effusion in the ionizer [28]. The ionization model is based on the assumption that multiple collisions of particles with hot ionizer walls enhance total ionization efficiency. It should be mentioned here that e.g. electron impact ionization could also provide very important contribution, especially in the case of hard-to-ionize elements [22].
As it has been said before, geometry of the ionizer cavity may be crucial for obtaining high ionization efficiency. It was e. g. demonstrated that spherical ionizers with small extraction openings achieve higher efficiency than tubular semi-opened ionizers in the case stable nuclides, while in the case of short-lived isotopes either hemispherical or conical ionizers are superb. In the current paper another shape of the ionizer is postulated. It is a flat disc cavity ionizer, which could be understood as a limiting case of cylindrical ionizer, but characterized by shorter length compared to its radius. The flat shape of the cavity should result in a large number of collisions of the particle traveling to the extraction opening. Simultaneously, small length of the cavity makes the penetration of the extraction field easier, leading to fast and effective ion evacuation. It should be also mentioned that manufacturing of the flat disc cavity seems to be easier task than in the case of spherical or even conical cavities, especially in the case of hard-to-machine materials. The paper contains brief description of the numerical model suitable for stable and very long-lived isotopes. Dependency of the ionization efficiency on the cavity shape (its elongation an radius) is studied. The influence of the extraction opening size on the efficiency is under consideration. The current-voltage curves obtained for different values of ionization coefficient are also presented and discussed. The changes of the ion yield due the ionizer temperature are investigated.
Numerical model
The numerical code used for calculation of ionization efficiency in a disc cavity is a trajectory tracking code similar to that considered in previous papers [21-25, 28]. The code follows trajectories of both ions and neutral particles that are confined in a hot cavity until they reach the extraction opening. A schematic cross-section through the simulated cavity-extraction electrode system is presented in Figure 1.
It should be stressed once again that standard cylindrical ionizer is characterized by r¡ << L condition, whilst for a flat disc cavity L should be much smaller or at least comparable to r. The results presented in the paper were obtained for L of order of 1 mm. The simulation system is discretized using by 100^400x400 rectangular numerical mesh with cell sizes Ax = Ay = Az = 0.05 mm. The electrostatic potential is worked out by numerical solving of Laplace equation (with boundary condition imposed
by electrodes) employing the successive overrelaxation method, as in [29, 30]. The classical equations of motion are integrated using 4th order Runge-Kutta algorithm. The forces used by the push subroutine are found by linear interpolation of the electric field values at six nodes nearest to the particle position. Starting velocities of particles depend on the ionizer temperature, while their initial directions are chosen randomly. A neutral particle or ion that hits the ionizer's hot internal surface could undergo ionization/neutralization with a probability (also called ionization coefficient) P, which is related to ionization degree by formula:
P = a/(1+a). (1)
The latter magnitude (a) is usually defined as the ratio of the numbers of ions and neutral atoms detaching from the hot surface and could be estimated using the Saha-Langmuir formula:
a = Gexp(-(V - k) / kT). (2)
In the above formula V and ^ are the ionization potential and the work function of the ionizer material, respectively, while the G prefactor depends e.g. on the atom-surface reflection coefficient. As it was already mentioned, the code follows the particles until they exit the cavity through the extraction hole and registers the total numbers of ions (N+) and of neutrals (No). The total ion source ionization efficiency is calculated as:
в =
N,
N + N
(3)
Figure 1 - Schematic view of the simulated system
One should have in mind that particles inside the cavity undergo many collisions with the ionizer hot walls. As the collision number reaches hundreds or even thousands, and charged particles are caught by the extraction field, the P. could be by orders
of magnitude larger than ionization probability during a single collision, predicted by the Saha-Langmuir formula (2).
Simulation results
The influence of ionizer cavity shape on ionization efficiency was investigated as the first. Namely, r¡ radius was changed starting from 0.8 mm up to 7 mm while the ionizer length L was kept constant (L = 1 mm). The flat extraction electrode on potential Vext = -2 kV) was set at the distance d = 1 mm from the extraction hole of radius re = 0.5 mm. The extraction channel length was chosen as l = 0.2, unless otherwise stated. The simulation timestep was chosen as 10-2 |is. Simulation runs were performed using 20000 test particles of 150 a.m.u. mass. The ionizer temperature was set to kT = 0.3 eV, unless otherwise stated. Calculation results are shown in Figure 2. As one can see, P coefficient changed over two and a half decades. Obtained total efficiency increases very fast with P. As one may expect ionization efficiency increases with the cavity radius and this effect is especially important in the case of hard-to-ionize substances (the case of P of order 0.01 or less). The total ionization efficiency could be several tens times larger than efficiency predicted by Saha-Langmuir law. For example, for P = 0.01 the total efficiency reaches 0.5 in the case of flat enough cavities (r¿ = 4 mm or more). The 90 % efficiency, marked as a dotted horizontal line, could be achieved for P larger than 0.06, which is significantly better than efficiencies achieved in the case of cylindrical [21] and hemispherical ionizers and comparable to that reached in spherical ionizers [24]. As it was already mentioned, the efficiency gain is reached due to the increasing number of particle-wall collisions in the case of flat cavities. As one can see in Figure 2. the average number of collisions increase dramatically with r-. For small values of P it reaches 700 in the case of very flat cavity while it is - 30 for the compact shape. It could be seen that <nhit> decreases with P - the easier is an atom to be ionized, the faster is caught by the extraction field.
On the other hand, in a flat disc shaped cavity the volume that is far away from the extraction field increases very fast with the ionizer radius. Decreasing probability of ion extraction leads to the saturation of ionization efficiency growth with r-. As it could be seen in Figure 3 an excessive increase
of r - above 7 mm in the considered case - is pointless, as it will not result in efficiency improvement.
Figure 2 - Total ionization efficiency as a function of ionization coefficient p for different values of ionizer radius ri
Figure 3 - Average number of particle-wall collisions calculated for different values of ionizer radius ri
The influence of extraction hole radius was also under investigation. Calculations were done for a disc cavity of diameter rt = 5 mm and L = 1 mm. Results for re changing in the range from 0.25 mm up to 4 mm are presented in Figure 4. One can see that total ionization efficiency increases with re. Once again that effect is more important for smaller p. The explanation is very similar to that given in the previous case: the smaller is the extraction hole, the more collision with the hot wall each particle undergoes, which enlarges the total ionization probability. The average number of collisions as function of p is shown in Figure 5. The increase of the average number of collisions by three orders of magnitude when the extraction hole radius is reduced to 0.2 mm is meaningful. One should, however, have in mind that the presented results were obtained for stable isotopes. More particle-wall collisions result in longer total time a particle stays in the ionizer, which may play a crucial role in the case of short-lived isotopes.
Figure 4 - Total ionization efficiency as a function of ionization coefficient p for different values of extraction opening radius re
Figure 5 - Average number of particle-wall collisions calculated for different values of extraction opening radius re
The role played by the disc cavity thickness L was also checked out. Simulations were performed for rt = 2 mm and re = 0.05 mm, all other parameters were as in the previous cases. The results are shown in Figure 6. One can see slight reduction of ionization efficiency for a very thin cavity (L = 0.5 mm). This could be probably due to the fact that the penetration of extraction field is limited in the case of such a flat cavity. For larger values of L the changes of efficiency are rather subtle, independently on the ionization coefficient.
1.0 1.5
I, mm
Figure 6 - Influence of disc-shaped cavity thickness L on the total ionization efficiency
Ionizer temperature may play an important role in achieving large ionization efficiencies. Figure 7 presents results of simulations of P. evolution with kT for a flat disk cavity. Simulations were performed for L = 1 mm, f = 5 mm and re = 0.5 mm. Two cases were considered: ionization coefficient in the lowest considered temperature kT = 0.2 eV was set either to 0.01 or 0.001. Ionization coefficient for higher temperatures was calculated using Saha-Langmuir formula assuming Vr^e = 0.9 eV. Changes of ionization coefficient with kT are shown in Figure 7b. It could bee seen that P increases several times in the considered temperature range. This is not the only magnitude that changes with ionizer temperature. Hollow symbols show the dependence of the desorbing particle initial temperature on kT - that change was also taken into account by the code. As it is shown in Figure 7a, total ionization efficiency increases with kT. The increase is nearly linear for P(0.2 eV) = 0.001, while a kind of saturation could be expected in the cases of higher initial ionization coefficient p.
0.8
0.6 -
со.
0.4 -
0.2 -
-■- (3, ( kT) for p ( 0.2 eV ) = 0.01 -•- ps (kT) for P ( 0.2 eV ) = 0.001
0.20 0.25 0.30
kT, eV
0.35
b
Figure 7 - Total ionization efficiency as a function of ionizer temperature kT (a) and dependence of ionization coefficient P and initial velocity of the desorbing particle on ionizer temperature kT (b)
As for other shapes of the ionizer, simulations of the current-voltage curves were done (see Figure 8).
Figure 8 - Current-voltage curves calculated for different values of ionization coefficient P
The case of L = 1 mm, f = 5 mm and re = 0.5 mm was considered (the same set of parameters as previously). Extraction voltage was increased up to 20 kV. One can observe very rapid increase of ion current with Vext rising up to ~ 1 kV. Above 1-2 kV a saturation of current-voltage curve is observed. This happens when the extraction field is able to catch and extract ions almost immediately after they are created. It should be kept in mind that excessive (more than 1-2 keV in the considered case) increase of the Vext is pointless -it does not lead to any significant increase of P. and may result in electrical breakdown. It is worth noting that the shape of the calculated curves is close to that obtained for tubular [23, 28] ionizers and to experimental curves measured for the tubular ionizers [31].
a
Conclusion
Stable isotope ionization in a novel kind of hot cavity, namely flat disc-shaped one, was considered in the paper. The brief description of the numerical model was given for completeness. Influence of the ionizer geometry on ionization efficiency was considered. It was shown that the efficiency increases with the disc ionizer radius reaching values higher than these obtained for elongated cylindrical or hemispherical ionizer. Ionization efficiency of 90 % could be reached even in the cases of the ionization coefficient of order 0.05. This effect is mostly to the large number of particle-ionizer wall collisions (reaching several hundreds) in the case of flat cavities.
Numerical simulation performed using the model showed also that:
- ionization efficiency of stable isotopes could be significantly enhanced by the reduction of the extraction opening radius;
- excessive reduction of that thickness (below 0.5 mm) may lower the total ionization efficiency;
- extraction voltage of 1-2 kV is sufficient to achieve almost maximal efficiencies at given geometry.
Simulations of the influence of the ionizer temperature on the total ionization efficiency showed that the efficiency increases with kT, mostly due to the fact that ionization probability usually increases with the surface temperature. Thus, it is crucial to maintain possible high ionizer temperature (for typical ionization conditions) and reduce chances of forming "cold spots'' that could reduce efficiency to a large extent.
Acknowledgments
The research was cofunded by th Plenipotent Representative of the Government of the Republic of Poland at JINR Dubna in the frame of the project № 75/07/2020.
References
1. Studer D., Maske L., Windpassinger P., Wendt K. Laser spectroscopy of the 1001-nm ground-state transition in dysprosium. Phys. Rev. A, 2018, vol. 98, pp. 042504. DOI: https://doi.org/10.1103/PhysRevA.98.042504
2. Duan Y., Danen R.E., Yan X., Steiner R., Cuadrado J., Wayne D., Majidi V., Olivares J.A. Characterization of an improved thermal ionization cavity source for mass spectrometry. Journal of the American Society for
Mass Spectrometry, 1999, vol. 10, pp. 917-1052. DOI: 10.1016/S1044-0305(99)00065-3
3. Maden C., Trinquier A., Fauré A.-L., Hubert A., Pointurier F., Rickli J., Bourdon B. Design of a prototype thermal ionization cavity source intended for isotope ratio analysis. International Journal of Mass Spectrometry, 2018, vol. 434, pp. 70-80.
DOI: 10.1016/j.ijms.2018.09.006
4. Babcock C., Day Goodacre T., Gottberg A. Target and Ion Source Development for Better Beams in the ARIEL Era. IOP Conf. Series: Journal of Physics: Conf. Series, 2018, vol. 1067, pp. 052019.
DOI: 10.1088/1742-6596/1067/5/052019
5. Alton G.D., Liu Y., Stracener D.W. High-efficiency target ion sources for radioactive ion beam generation. Rev. Sci. Instrum., 2006, vol. 77, pp. 03A711. DOI: 10.1063/1.2173968
6. Köster U., Arndt O., Bouquerel E., Fedo-seyev V.N., Franberg H., Joinet A., Jost C., Kerki-nes I.S.K., Kirchner R. The TARGISOL Collaboration, Progress in ISOL target-ion-source-system. Nucl. Instrum Meth. B, 2008, vol. 266, pp. 4229-4239.
DOI: 10.1016/j.nimb.2008.05.152
7. Woo H.J., Kang B.H., Tshoo K., Seo C.S., Hwang W., Park Y.H., Yoon J.W., Yoo S.H., Kim Y.K., Jang D.Y. Overview of the ISOL facility for the RISP. Journal of the Korean Physical Society, 2015, vol. 66, pp. 443-448. DOI: 10.3938/jkps.66.443
8. Beyer G.J., Herrmann E., Piotrowski A., Rai-ko V.I., Tyroff H. A new method for rare-earth isotope separation. Nucl. Instrum. Meth. 1971, vol. 96, pp. 437439. DOI: 10.1016/0029-554X(71)90613-6
9. Johnson P.G., Bolson A., Henderson C.M. A high temperature ion source for isotope separators. Nucl. Instrum. Meth, 1973, vol. 106, pp. 83-87.
DOI: 10.1016/0029-554X(73)90049-9
10. Liu Y., Jost C.U., Mendez II A.J., Stracener D.W., Williams C.L., Gross C.J., Grzywacz R.K., Ma-durga M., Miernik K., Miller D. On-line commissioning of the HRIBF resonant ionization laser ion source. Nucl. Instrum. and Meth. B, 2013, vol. 298, pp. 5-12.
DOI: 10.1016/j.nimb.2012.12.041
11. Lecesne N. Laser ion sources for radioactive beams. Rev. Sci. Instrum., 2012, vol. 83, pp. 02A916. DOI: 10.1063/1.3681148
12. Henares J.L., Lecesne N., Hijazi L., Bastin B., Kron T., Lassen J., Le Blanc F., Leroy R., Osmond B., Raeder S., Schneider F., Wendt K. Hot-cavity studies for the Resonance Ionization Laser Ion Source. Nucl Instr. Meth. B, 2016, vol. 830, pp. 520-525.
DOI: 10.1016/j.nima.2015.10.061
13. Day Goodacre T., Billowes J., Catherall R., Cocolios T.E., Crepieux B., Fedorov D.V., Fedosse-ev V.N., Gaffney L.P., Giles T., Gottberg A., Lynch K.M.,
Marsh B.A., Mendonça T.M., Ramos J.P., Rossel R.E., Rothe S., Sels S., Sotty C., Stora T., Van Beveren C., Veinhard M. Blurring the boundaries between ion sources: The application of the RILIS inside a FEBIAD type ion source at ISOLDE. Nucl Instr. Meth. B, 2016, vol. 376, pp. 39-45. DOI: 10.1016/j.nimb.2016.03.005
14. Kalinnikov V.G., Gromov K.Ya., Janicki M., Yushkevich Yu.V., Potempa A.W., Egorov V.G., Byst-rov V.A., Kotovsky N.Yu., Evtisov S.V. Experimental complex to study nuclei far from the beta-stability line -ISOL-facility YASNAPP-2. Nucl. Instr. and Meth. B, 1992, vol. 70, pp. 62-68.
DOI: 10.1016/0168-583X(92)95910-J
15. Zhai L., Deng H., Wei G., Li Z., Wang C., Li X., Zhou G., Su Y., Zhang Z. A new, ohmic-heating based thermal ionization cavity source for mass spectrometry. International Journal of Mass Spectrometry, 2011, vol. 305, pp. 45-49. DOI: 10.1016/j.ijms.2011.05.015
16. Eléon C., Jardin P., Gaubert G., Saintlaurent M., Alcantaranunez J., Alvesconde R. Development of a surface ionization source for the production of radioactive alkali ion beams in SPIRAL. Nucl. Instr. and Meth. B, 2008, vol. 266, pp. 4362-4367.
DOI: 10.1016/j.nimb.2008.05.067
17. Reponen M., Moore I.D., Pohjalainen I., Rothe S., Savonen M., Sonnenschein V., Voss A. An inductively heated hot cavity catcher laser ion source. Rev Sci Instrum. 2015, vol. 86, pp. 123501.
DOI: 10.1063/1.4936569
18. Alton G.D., Zhang Y. A fast effusive-flow vapor-transport system for ISOL-based radioactive ion beam facilities. Nucl. Instrum. and Meth. A, 2005, vol. 539, pp. 540-546. DOI: 10.1016/j.nima.2004.11.027
19. Alton G.D., Liu Y., Zaim H., Murray S.N. An efficient negative surface ionization source for RIB generation. Nucl. Instrum. and Meth. B, 2003, vol. 211, pp. 425-435. DOI: 10.1016/S0168-583X(03)01365-X
20. Hausladen P.A., Weisser D.C., Lobanov N.R., Fifield L.K., Wallace H.J. Simple concepts for ion source improvement. Nucl. Instrum. and Meth. B, 2002, vol. 190, pp. 402-404.
DOI: 10.1016/S0168-583X(01)01307-6
21. Turek M., Pyszniak K., Drozdziel A., Sie-lanko J. Ionization efficiency calculations for cavity thermoionization ion source. Vacuum, 2008, vol. 82, pp. 1103-1106. DOI: 10.1016/j.vacuum.2008.01.025
22. Turek M., Pyszniak K., Drozdziel A. Influence of electron impact ionization on the efficiency of thermoemission ion source. Vacuum, 2009, vol. 83, pp. S260-S263. DOI: 10.1016/j.vacuum.2009.01.077
23. Turek M., Drozdziel A., Pyszniak K., Maczka D., Slowinski B. Simulations of ionization in a hot cavity surface ion source. Rev. Sci. Instrum, 2012, vol. 83, pp. 023303. DOI: 10.1063/1.3685247
24. Turek M. Modeling of Ionization in a Spherical Surface Ionizer. Acta Phys. Pol. A, 2011, vol. 120, pp. 188-191. DOI: 10.12693/APhysPolA.120.188
25. Turek M. Ionisation Efficiency in Conical Hot Cavities. Acta Phys. Pol. A, 2017, vol. 132, pp. 259-263. DOI: 10.12693/APhysPolA.132.259
26. Maden C., Baur H., Fauré A.-L., Hubert A., Pointurier F., Bourdon B. Determination of ionization efficiencies of thermal ionization cavity sources by numerical simulation of charged particle trajectories including space charge. Int. J. Mass Spectr, 2016, vol. 405, pp. 39-49. DOI: 10.1016/j.ijms.2016.05.013
27. Liu Y., Batchelder J.C., Galindo-Uribarri A., Chu R., Fan S., Romero-Romero E., Stracener D.W. Ion source development for ultratrace detection of uranium and thorium. Nucl. Instrum. and Meth. B, 2015, vol. 361, pp. 267-272. DOI: 10.1016/j.nimb.2015.04.081
28. Turek M. Ionization of short-lived isotopes in a hot cavity - Numerical simulations. Vacuum, 2014, vol. 104, pp. 1-12. DOI: 10.1016/j.vacuum.2013.12.016
29. Hadjidimos A. Successive overrelaxation (SOR) and related methods. Journal of Computational and Applied Mathematics, 2000, vol. 123, pp. 177-199. DOI: 10.1016/S0377-0427(00)00403-9
30. Press W.H., Teukolsky S.A., Vetterling W.T., Flannery B.P. Numerical recipes in FORTRAN (2nd ed.): The art of scientific computing, 1992, Cambridge University Press New York.
31. Latuszynski A., Maczka D. High temperature cavity thermo-ionizer. Vacuum, 1998, vol. 51, pp. 109112. DOI: 10.1016/S0042-207X(98)00142-0
