105 Computational and Experimental Research in Materials and Renewable Energy (CERiMRE) Volume 5, Issue 2, page 105-113 eISSN : 2747-173X Submitted : June 29, 2022 Accepted : September 8, 2022 Online : November 11, 2022 DOI: 10.19184/cerimre.v5i2.32069 Innovative Compact Molten Salt Reactor (ICMSR) Analysis for Mo-99 Production Iza Shafera Hardiyanti1,a, A Suparmi1 and Andang Widi Harto2 1Physics Department, Universitas Sebelas Maret, Surakarta, Indonesia 2Nuclear Physics Department and Engineering Physics Department, Universitas Gadjah Mada, Yogyakarta, Indonesia aEmail: izashafera@student.uns.ac.id Abstract. Mo-99 isotope production calculation in the ICMSR (Innovative Compact Molten Salt Reactor) with the computer code MCNP6 has been carried out. ICMSR is a conceptual design of the MSR type reactor that uses NaF-ThF4-UF4 fuel with an enrichment of 235U of 19.75%. This reactor operates in thermal neutron spectrum with a graphite moderator. ICMSR is a power reactor that produces Mo-99 as a by- product. Calculations carried out for 12 days of operation show that the reactor condition is still critical so that there will be no intervention from refueling. The total Mo-99 produced until the 12th day is 9.118 x 106 Ci. Mo-99 can be extracted from the reactor as long as the power reactor is operating so it will be economically advantageous. Keywords: 99Mo, isotope production, ICMSR, MCNP, criticality Introduction Almost all hospitals in the world use radioisotopes in medicine. Most of it is for diagnosis. Technetium-99m (𝑇1/2 = 6.02 β„Ž) is the most widely used radioisotope [1]. More than 80% of radio- diagnostic procedures in the world use this isotope [2, 3, 4]. Technetium-99m can be generated from the decay of the parent nuclide Molybdenum-99 (𝑇1/2 = 66 β„Ž) [5]. The Mo-99 will decay to Tc-99m which then decays to its ground state Tc-99 by emitting 143 keV 𝛾 βˆ’rays which are known to have the same wavelength as diagnostic x-ray devices. The Tc-99m is ideal for diagnostic because of its short half-life so it is safe for patients [6]. However, this also makes it difficult to supply to medical centers. To solve the problem of sending Tc-99m to end user, a generator is used. A generator is a system that contains a radioactive parent nuclide with relatively long half- life that decays into short-lived daughter nuclides [7]. The generator used here refers to the Mo- 99 generator. The five main supply countries of Mo-99 are Canada, Belgium, France, Netherlands and South Africa. Most isotopes are generated from research reactors [8]. In 2008, shortages of Tc-99m began to emerge as the two Mo-99 reactors that supply 66% of the world’s supply was shut down for maintenance. That reactor is National Research Universal (NRU) in Canada and Petten Nuclear Reactor (HFR) in Netherlands [9, 10]. Discontinuation of isotope supply from both reactors led to many delays and cancellations of medical procedures. The situation is getting worse: firstly, the reactor is old, secondly, using highly enriched uranium (HEU) which is not compatible with nonproliferation [11]. In order to meet the need of the domestic market, several research and development were carried out. One of the parties authorized to produce Mo-99 is PT INUKI. However, the production is not 106 Computational and Experimental Research in Materials and Renewable Energy (CERiMRE) Volume 5, Issue 2, page 105-113 eISSN : 2747-173X Submitted : June 29, 2022 Accepted : September 8, 2022 Online : November 11, 2022 DOI: 10.19184/cerimre.v5i2.32069 sufficient to meet national needs [12]. From the results of the research reactor (RSG GAS), the production capacity of Mo-99 is 300 Ci per batch (week) [13]. The development of reactors to produce radioisotopes is also being carried out, i.e SAMOP and CAMOLYP [12, 14, 15]. An alternative technology that can support the fulfillment of the growing demand for Mo-99 is Molten Salt Reactor (MSR) [16]. The big advantage of using MSR for the production of Mo-99 is the possibility of online extraction [17]. In this study, we propose Innovative Compact Molten Salt Reactor (ICMSR). ICMSR is a power reactor, not a special reactor to produce Mo-99. However, during its operation, this reactor produces by-products, one of which is Mo-99. So this research was conducted to calculate and find out how much amount of Mo-99 can be produced from reactor operation for support reactors such as SAMOP and CAMOLYP to supply this isotope. Theoretical Background The production of Mo-99 in general can be carried out in two ways, namely: fission fragments from U-235 and neutron capture reactions from Mo-98. When the U-235 is hit by thermal neutrons, the atomic nucleus to become unstable and lose its shape, then split into smaller elements while releasing heat energy and releasing 2-3 new neutrons. One of those elements is Mo-99. In addition, neutron capture by M-98 can also produce Mo-99 accompanied by 𝛾 βˆ’ray emission. However, the production of Mo-99 from capture reaction is not used for commercial scale production because it is inefficient [18]. Thus, the production of this isotope depends on the fission reactions in nuclear reactors. Molybdenum is considered a noble metal in molten salt so it is easily reduced in molten [19]. The reactor that uses molten salt is the Molten Salt Reactor (MSR). The ability of MSR to produce Mo- 99 has been carried out and shows good results [20]. The accumulation of each isotope in a nuclear reactor depends on fission, decay, activation and extraction [17]. This can apply to any isotope, not just Mo-99. Firstly is fission. The rate equation for the fission reaction can be written as: �̂�𝑓 Γ— π‘Œπ‘– = 𝐹𝑃 π‘”π‘’π‘›π‘’π‘Ÿπ‘Žπ‘‘π‘–π‘œπ‘› π‘Ÿπ‘Žπ‘‘π‘’ (1) Where, �̂�𝑓 = π‘π‘œπ‘€π‘’π‘Ÿ 𝑑𝑒𝑛𝑠𝑖𝑑𝑦 π‘’π‘›π‘’π‘Ÿπ‘”π‘¦ π‘Ÿπ‘’π‘™π‘’π‘Žπ‘ π‘’/π‘“π‘–π‘ π‘ π‘–π‘œπ‘› , π‘Œπ‘– is the fission product yield per fission. The subscript 𝑖 refers to the isotope of interest. Secondly is decay. The decay equation can be divided into two, namely the decay of the parent isotope and the decay of the isotope of interest. The generation of isotopes from parent isotopes can be written as: βˆ‘ π›½π‘πœ†π‘π‘π‘ = π‘‘π‘’π‘π‘Žπ‘¦ π‘œπ‘“ π‘π‘Žπ‘Ÿπ‘’π‘›π‘‘ π‘–π‘ π‘œπ‘‘π‘œπ‘π‘’π‘ π‘ (2) Where, 𝛽 is the branch fraction, πœ† is the decay constant and 𝑁 is the number density. The subscript 𝑝 refers to the parent isotopes. This equation is especially important for isotopes which have short half-lives. This isotope will remove and decay into its daughter isotopes. The decay of 107 Computational and Experimental Research in Materials and Renewable Energy (CERiMRE) Volume 5, Issue 2, page 105-113 eISSN : 2747-173X Submitted : June 29, 2022 Accepted : September 8, 2022 Online : November 11, 2022 DOI: 10.19184/cerimre.v5i2.32069 the isotope of interest is also needed if isotope unstable and decay. This equation can be written as: βˆ’πœ†π‘– 𝑁𝑖 = π‘‘π‘’π‘π‘Žπ‘¦ π‘œπ‘“ π‘‘β„Žπ‘’ π‘–π‘ π‘œπ‘‘π‘œπ‘π‘’ π‘œπ‘“ π‘–π‘›π‘‘π‘’π‘Ÿπ‘’π‘ π‘‘ (3) Where, πœ† is the decay constant and 𝑁 is the number density. The subscript 𝑖 refers to the isotope of interest. Thirdly is activation. Activation occurs when a neutron is absorbed by an isotope. Activation of a nuclear reaction depends on the type of isotope, the energy of the neutron and the isotopic number density. In reactors where the fuel is considered homogeneous and has uniform neutron energy such as MSR, the activation equation can be approximated by the following equation �̂�𝑖 π‘₯ = 𝑁𝑖 πœŽπ‘– π‘₯ πœ™ = π‘Ÿπ‘’π‘Žπ‘π‘‘π‘–π‘œπ‘› π‘Ÿπ‘Žπ‘‘π‘’ (4) Where, οΏ½Μ‚οΏ½ is reaction rate, 𝜎 is microscopic cross section, πœ™ is neutron flux. The superscript π‘₯ refers to nuclear reaction of interest. The neutron flux in this equation can be eliminated by looking at the relationship between the reaction rates of two different isotopes by assuming that the flux is mono energetic. �̂�𝑖 π‘₯ π‘π‘–πœŽπ‘– π‘₯ = οΏ½Μ‚οΏ½ 𝑗 𝑦 π‘π‘—πœŽπ‘— 𝑦 = πœ™ (5) So the activation equation can be rewritten as follows: �̂�𝑖 π‘₯ = π‘π‘–πœŽπ‘– π‘₯ π‘π‘—πœŽπ‘— 𝑦 �̂�𝑗 𝑦 (6) The nuclear reaction we want to observe on activation is an isotope capture reaction in a fission fuel, so the equation becomes: �̂�𝑖 𝑐 = π‘π‘–πœŽπ‘– 𝑐 π‘π‘“π‘’π‘’π‘™πœŽπ‘“π‘’π‘’π‘™ 𝑓 οΏ½Μ‚οΏ½ 𝑓 (7) Where, superscript 𝑐 refers to capture reaction and superscript 𝑓 refers to fission reaction. Activation can do two things: generate isotopes and remove isotopes. Equation (7) refers to removal of the isotope due to activation. Meanwhile, the generation of isotopes produced from the activation product precursors can be written by replacing the subscript 𝑖 with the subscript π‘Ž. The total activation equation is as follows: βˆ’ π‘π‘–πœŽπ‘– 𝑐 π‘π‘“π‘’π‘’π‘™πœŽπ‘“π‘’π‘’π‘™ 𝑓 οΏ½Μ‚οΏ½ 𝑓 + π‘π‘ŽπœŽπ‘Ž 𝑐 π‘π‘“π‘’π‘’π‘™πœŽπ‘“π‘’π‘’π‘™ 𝑓 οΏ½Μ‚οΏ½ 𝑓 = π‘‘π‘œπ‘‘π‘Žπ‘™ π‘Žπ‘π‘‘π‘–π‘£π‘Žπ‘‘π‘–π‘œπ‘› 𝑑𝑒𝑛𝑠𝑖𝑑𝑦 (8) Later the isotope of interest will be extracted, the extraction of isotopes such as Mo-99 which is classified as a noble metal fulfills the equation 108 Computational and Experimental Research in Materials and Renewable Energy (CERiMRE) Volume 5, Issue 2, page 105-113 eISSN : 2747-173X Submitted : June 29, 2022 Accepted : September 8, 2022 Online : November 11, 2022 DOI: 10.19184/cerimre.v5i2.32069 βˆ’π‘˜π‘ 𝐴𝑁𝑖 = 𝑒π‘₯π‘‘π‘Ÿπ‘Žπ‘π‘‘π‘–π‘œπ‘› (9) Where, π‘˜π‘ is the mass transfer coefficient, 𝐴 is the surface area available for extraction and N is number density. So that the desired isotope accumulation equation in the reactor can be written: 𝑑𝑁𝑖 𝑑𝑑 = �̂�𝑓 π‘Œπ‘– βˆ’ πœ†π‘– 𝑁𝑖 + βˆ‘ π›½π‘πœ†π‘π‘π‘π‘ βˆ’ π‘π‘–πœŽπ‘– 𝑐 π‘π‘“π‘’π‘’π‘™πœŽπ‘“π‘’π‘’π‘™ 𝑓 οΏ½Μ‚οΏ½ 𝑓 + π‘π‘ŽπœŽπ‘Ž 𝑐 π‘π‘“π‘’π‘’π‘™πœŽπ‘“π‘’π‘’π‘™ 𝑓 οΏ½Μ‚οΏ½ 𝑓 βˆ’ π‘˜π‘ 𝐴𝑁𝑖 (10) Materials and Methods The type of MSR proposed in this study is the conceptual design of ICMSR (Innovative Compact Molten Salt Reactor). ICMSR is a power reactor with NaF – ThF4 – UF4 fuel. ICMSR uses 19.75% uranium enrichment to avoid proliferation. The core of this reactor is a HeteType core, which is a heterogeneous graphite fuel and moderator layout, allowing molten salt to flow through many graphite blocks. The HeteType core design has the advantage that it requires the smallest initial fuel loading and the lowest radiotoxicity of the remaining spent fuel in long term operation [16]. In the ICMSR design, the heat exchanger (HE) is integrated with the core, the aim is to reduce fuel inventory in the core and reduce the size of the cooling system. ICMSR is equipped with a noble and volatile gas extraction system. Noble gases can be removed by purging an inert gas into the molten salt. With the extraction system, the utility of the neutron is higher because neutron poisons such as Xe-135 no longer exist. Higher neutron utility means more collisions with uranium-235, resulting in higher production of fission products such as Mo-99. ICMSR parameter and geometry is show in Table 1 and Figure 1. Table 1. Input parameter of ICMSR Parameter Value (unit) Power 75 MWe; 187.5 MWth Temperature (during operation) 900K Fuel NaF-ThF4-UF4 U-235 enrichment 19.75% Core diameter; high 220 cm; 246 cm Moderator Graphite Moderator diameter; high 28 cm; 150 cm Intermediate coolant NaF-KF Vessel reactor Hastelloy-N Vessel diameter; high 400 cm; 10.51 m 109 Computational and Experimental Research in Materials and Renewable Energy (CERiMRE) Volume 5, Issue 2, page 105-113 eISSN : 2747-173X Submitted : June 29, 2022 Accepted : September 8, 2022 Online : November 11, 2022 DOI: 10.19184/cerimre.v5i2.32069 Figure 1. Configuration of modular ICMSR system The fuel depletion analysis was carried out for 12 days. The scenario of fuel depletion by removing gas fission products is carried out with the OMIT command. The analysis will only carried out at a power of 75 MWe (thermodynamic efficiency = 40%) where this value corresponds to the operating power of the ICMSR. Although in this reactor operation all data on isotope depletion information is available, this paper will only focus on the production of Mo-99. The flowchart of this study is shown in Figure 2. Figure 2. Flowchart of the study 110 Computational and Experimental Research in Materials and Renewable Energy (CERiMRE) Volume 5, Issue 2, page 105-113 eISSN : 2747-173X Submitted : June 29, 2022 Accepted : September 8, 2022 Online : November 11, 2022 DOI: 10.19184/cerimre.v5i2.32069 The method used for calculation of the Mo-99 isotope is a computation method. This study has used MCNP6. The MCNP (Monte Carlo N Particle) is a computer program that developed since 1963 at Los Alamos National Laboratory (LANL), United States. MCNP6 is the latest version of MCNP which has been successfully developed by LANL. MCNP6 is the result of merging MCNP5 and MCNPX into a single computer code program that contains the features of both. In addition, MCNP6 also includes new calculation capabilities such as calculation of reactor kinetic parameters, unstructured Tally FMESH, and sensitivity analysis card [21]. The monte carlo method is a calculation method that theoretically imitates a statistical or random process. This method is usually used to solve complex problems, which can no longer be solved by deterministic calculation methods. This method is implemented by simulating every single probabilistic event that occurs in a process, one by one in sequence. The probability distribution that occurs in each event is calculated randomly. In this method, it takes a lot of repetition, so that the whole phenomenon being simulated can be fully and realistically described. Multiplication factor calculation in MCNP can be done with several estimators [22]: a. Collision estimator π‘˜π‘’π‘“π‘“ 𝐢 = 1 𝑁 𝛴𝑖 π‘Šπ‘– [ π›΄π‘˜π‘“π‘˜οΏ½Μ…οΏ½π‘˜πœŽπ‘“π‘˜ π›΄π‘˜π‘“π‘˜πœŽπ‘‡π‘˜ ] (11) Where, 𝑖 is summed over all collision in a cycle where fission is possible; π‘˜ is summed over all nuclides of the material involved in the π‘–π‘‘β„Ž collision; πœŽπ‘‡π‘˜ is total microscopic cross section; πœŽπ‘“π‘˜ is microscopic fission cross section; οΏ½Μ…οΏ½π‘˜ is average number of prompt or total neutrons produced per fission by the collision nuclide at the incident energy; π‘“π‘˜ is atomic fraction for nuclide k; N is nominal source size for cycle; π‘Šπ‘– is weight of particle entering collision. b. Absorption estimator π‘˜π‘’π‘“π‘“ 𝐴 = 1 𝑁 𝛴𝑖 π‘Šπ‘– οΏ½Μ…οΏ½π‘˜ πœŽπ‘“π‘˜ πœŽπ‘π‘˜+ πœŽπ‘“π‘˜ (12) Where, 𝑖 is summed over each analog absorption event in the π‘˜π‘‘β„Ž nuclide. c. Track length estimators π‘˜π‘’π‘“π‘“ 𝑇𝐿 = 1 𝑁 𝛴𝑖 π‘Šπ‘– 𝜌 𝑑 π›΄π‘˜ π‘“π‘˜ οΏ½Μ…οΏ½π‘˜ πœŽπ‘“π‘˜ (13) Where; 𝑖 is summed over all neutron trajectories; 𝜌 is atomic density in the cell; 𝑑 is the trajectory track length from the last event. The combination average of the three estimators is: π‘˜π‘’π‘“π‘“ 𝑏𝑒𝑠𝑑 = βˆ‘ π‘˜π‘’π‘“π‘“ 𝑖 βˆ—π‘–,π‘–β‰ π‘—β‰ π‘˜ (πΆπ‘—π‘—βˆ—πΆπ‘˜π‘˜βˆ’πΆπ‘—π‘—βˆ—πΆπ‘–π‘˜βˆ’πΆπ‘˜π‘˜βˆ—πΆπ‘–π‘—+πΆπ‘—π‘˜βˆ—(𝐢𝑖𝑗+πΆπ‘–π‘˜βˆ’πΆπ‘—π‘˜)) βˆ‘ πΆπ‘—π‘—βˆ—πΆπ‘˜π‘˜βˆ’πΆπ‘—π‘—βˆ—πΆπ‘–π‘˜βˆ’πΆπ‘˜π‘˜βˆ—πΆπ‘–π‘—+πΆπ‘—π‘˜βˆ—(𝐢𝑖𝑗+πΆπ‘–π‘˜βˆ’πΆπ‘—π‘˜)𝑖,π‘–β‰ π‘—β‰ π‘˜ (14) Where: 𝑖 = C, A, and TL 111 Computational and Experimental Research in Materials and Renewable Energy (CERiMRE) Volume 5, Issue 2, page 105-113 eISSN : 2747-173X Submitted : June 29, 2022 Accepted : September 8, 2022 Online : November 11, 2022 DOI: 10.19184/cerimre.v5i2.32069 𝐢𝑖𝑗 = covariance = 1 𝑀 βˆ‘ π‘˜π‘’π‘“π‘“π‘š 𝑖 π‘˜π‘’π‘“π‘“π‘š 𝑗 βˆ’ ( 1 𝑀 βˆ‘ π‘˜π‘’π‘“π‘“π‘š 𝑖 π‘š ) ( 1 𝑀 βˆ‘ π‘˜π‘’π‘“π‘“π‘š 𝑗 π‘š )π‘š m is the index for the mth active cycle and M is the number of active cycles. Results and Discussion The burnup calculation with the extraction of noble gas for 12 days has been done. The results showed that during 12 days of operation, the reactor was still in critical condition. This indicates that there is no need to add fuel online at that time. This condition can be achieved with a mole fraction of UF4 of 5.6% with an enrichment of 19.75%. Basically, the higher the two values, the higher the criticality value. But keep in mind, that the reactor will be safe when the criticality value is 1 or the reactivity value is zero and continues to be rated constant. The results show that keff value is slightly above 1. This supercritical condition is still relatively reasonable where the excess reactivity can be used to increase power if necessary. The results of the criticality value during operation can be seen in Figure 3. Figure 3. The criticality levels per day of ICMSR reactor operation Mo-99 is one of the major radionuclides in FPs. Table 2 shows the accumulation of Mo-99 in the ICMSR for 12 days. On day 0, Mo-99 has not yet formed; meaning there is only fuel in the reactor core. The highest increase in Mo-99 production was on the first day at 2.136 x 106 Ci. This production is the result of the fission of fuels. On six days, the resulting average is around 1.245 x 106 Ci, while on the 6th to 12th day the average generated is 0.275 x 106 Ci. The total Mo-99 produced until the 12th day is 9.118 x 106 Ci. In ICMSR, Mo-99 can be extracted during power reactor operation through fission product extraction. So it is very profitable from an economic point of view. 112 Computational and Experimental Research in Materials and Renewable Energy (CERiMRE) Volume 5, Issue 2, page 105-113 eISSN : 2747-173X Submitted : June 29, 2022 Accepted : September 8, 2022 Online : November 11, 2022 DOI: 10.19184/cerimre.v5i2.32069 Table 2. Specific production of Mo-99 as by product of ICMSR power reactor Day Specific production (Ci) Day Specific production (Ci) 1 2.136 x 106 7 7.946 x 106 2 3.794 x 106 8 8.306 x 106 3 5.083 x 106 9 8.597 x 106 4 6.085 x 106 10 8.822 x 106 5 6.863 x 106 11 8.893 x 106 6 7.472 x 106 12 9.118 x 106 Conclusions The criticality of the ICMSR reactor for 12 days of operation shows the rector is still in a critical condition so that there will be no intervention for refueling or power if Mo-99 extraction is carried out in that time span. Mo-99 production for six days of operation showed a high average of 1.245 x 106 Ci. Meanwhile, on the 6th to 12th day the average generated is 0.275 x 106 Ci. The total Mo- 99 produced until the 12th day is 9.118 x 106 Ci. So it can be state that the ICMSR power reactor produces an isotope that is useful for medical purposes, namely Mo-99 as a by product. References [1] IAEA, 1999, Production Technologies for Molybdenum-99 and Technetium-99m, IAEA-TECDOC-1605. [2] M. R. A. Pillai, A. Dash, and F. F. R. K. Jr, 2015, Diversification of 99Mo/99mTc separation: non-fission reactor production of 99Mo as a strategy for enhancing 99mTc availability, J Nucl Med, 56 (1), 159-161. [3] IAEA, 2015, Feasibility of Producing Molybdenum-99 on a Small Scale Using Fission of Low Enriched Uranium or Neutron Activation of Natural Molybdenum, Technical Report Series No. 478, STI/DOC/010/478 ISBN 978-92-0-114713-4. [4] I. N. Goldman, N. Ramamoorthy, and P, 2008, presented at RERTR Int. Meet. [5] WNA, 2013, Radioisotopes in medicine. [6] M. M. Khalil, 2011, Basic Sciences of Nuclear Medicine, Springer, United Kingdom [7] J. J. P. De lima, 2011, Nuclear Medicine Physics, Taylor and Francis Group, US [8] S. Cherry, J. Sorenson, and M. Phelps, 2012, Physics in Nuclear Medicine, Elsevier, Amsterdam, Netherlands [9] G. S. Thomas and J. Maddahi, 2010, The technetium shortage, J Nucl Cardiol, 17 (6), 993 – 998. [10] T. J. Ruth, 2014, The Medical Isotope Crisis: How We Got Here and Where We Are Going, Journal of Nuclear Medicine Technology, 42 (4), 245 – 248. [11] Y. Nagai, and Y. Hatsukawa, 2009, Production of 99Mo for Nuclear Medicine by 100Mo (n, 2n)99Mo, Journal of the Physical Society of Japan, 78 (3). 113 Computational and Experimental Research in Materials and Renewable Energy (CERiMRE) Volume 5, Issue 2, page 105-113 eISSN : 2747-173X Submitted : June 29, 2022 Accepted : September 8, 2022 Online : November 11, 2022 DOI: 10.19184/cerimre.v5i2.32069 [12] Syarip, T. Sutondo, E. T. Budisantoso, and E. Susiantini, 2018, Design and Development of Subcritical Reactor by Using Aqueous Fuel for 99Mo Production, A. Physical and Computational Sciences, 55 (1), 21 -26. [13] S. Kuntjoro, 2016, Analysis of Mo-99 Isotope Activation in the RSG-GAS Reactor (in Indonesian), Sigma Epsilon, 20, 13 -20. [14] M. I. W. Farezza, and Syarip, 2018, Mo-99 Isotope Production Calculation of SAMOP Reactor Experimental Facility, International Conference on Computation in Science and Engineering. [15] D. Bartolomeus and S. Syarip, 2019, Analysis of Mo-99 production as function of CAMOLYP reactor power, 4th Padjajaran International Physics Symposium. [16] R. J. Sheu, C. C. Chao, O. Feynberg, amd Y. W. H. Liu, 2014, A fuel depletion analysis of the MSRE and three conceptual small molten-salt reactors for Mo-99 production, Annals of nUclear Energy, 71. [17] M. N. Stoddart, J. N. Harb, and M. J. Memmott, 2019, Numerical analysis of isotope production in molten salt reactors: A case study for Molybdenum-99 production, Annals of Nuclear Energy, 129, 56 – 61. [18] D. A. Pakholik, O. Y. Kochnov, V. V. Kolesov, and V. V. Fomichev, 2021, Increasing the production of the Mo-99 isotope by modernizing the design of targets irradiated in the experimental channels of the VVR-c reactor, Nuclear Energy and Technology, 7 (4), 291 – 295. [19] W. Grimes, 1970, Molten-salt reactor chemistry, Nucl Apply Techn, 8 (2), 137 – 155. [20] D. Y. Chuvilin and V. A. Zafryadskii, 2009, New method of producing 99Mo in molten salt fluoride fuel, Atomic Energy, 107 (3), 185 – 193. [21] C. J. Werner, J. Armstrong, F. B. Brown, J. S. Bull, L. Casswell, L. J. Cox,…, B. C. Kiedrowski, 2017, MCNP User’s Manual Code Version 6.2, Los Alamos National Laboratory, 746. [22] X-5 Monte Carlo Team, 2008, MCNP-A General Monte Carlo N – Particle Transport Code, Version 5 Volume I: Overview and Theory, Los Alamos National Laboratory.