Next Article in Journal
Impedance Spectroscopy Sensing Material Properties for Self-Tuning Ratio Control in Pharmaceutical Industry
Next Article in Special Issue
A Novel Method for Spectrometry Based on Imaging Systems
Previous Article in Journal
Boundary-Aware Hashing for Hamming Space Retrieval
Previous Article in Special Issue
Study of Surface Emissions of 220Rn (Thoron) at Two Sites in the Campi Flegrei Caldera (Italy) during Volcanic Unrest in the Period 2011–2017
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study of 222−220Rn Measurement Systems Based on Electrostatic Collection by Using Geant4+COMSOL Simulation

1
Department of Mathematics and Physics, University of Campania “L. Vanvitelli”, Viale Lincoln 5, 81100 Caserta, Italy
2
Department of Agricultural Sciences, University of Naples “Federico II”, Portici, Italy
3
National Institute of Nuclear Physics, Via Cinthia 21, 80126 Napoli, Italy
*
Author to whom correspondence should be addressed.
Submission received: 21 October 2021 / Revised: 20 December 2021 / Accepted: 1 January 2022 / Published: 5 January 2022

Abstract

:
Using Monte Carlo (with Geant4) and COMSOL simulations, the authors have defined a useful tool to reproduce the alpha spectroscopy of 222Rn, 220Rn and their ionized daughters by measurement systems based on electrostatic collection on a silicon detector, inside a metallic chamber. Several applications have been performed: (i) simulating commercial devices worldwide used, and comparing them with experimental theoretical results; (ii) studying of realization of new measurement systems through investigation of the detection efficiency versus different chamber geometries. New considerations and steps forward have been drawn. The present work is a novelty in the literature concerning this research framework.

1. Introduction

Radon (Rn) is a natural radioactive gas present in the environment. Two are the main isotopes: 222Rn (radon) with a decay half-life of 3.823 days, and 220Rn (thoron) with a half-life of 55.8 s. Both are carefully studied worldwide by the scientific community especially for radioprotection purposes [1,2] and for geophysical phenomena tracking [3,4]. These goals require measurement systems suitable for monitoring the gas isotopes. Various methods for Rn detection have been developed throughout the years, those based on solid-state nuclear track detection and electrostatic collection are the most widespread and used nowadays [5]. The first method utilizes plastic devices sensitive to α-radiation emitted by Rn and decay products, for mean long-term monitoring [6]. In the second method, naturally ionized Rn daughter, inside a metallic chamber put to a high voltage, are electrostatically driven on to the surface of a silicon detector set to ground potential, enabling the collection of full alpha energies and obtaining high resolution α-spectrometry. The electrostatic collection systems have optimal detection efficiencies; they allow the continuous detection of the Rn by distinguishing both gas isotopes with their progeny [7]. Currently, three are the main commercial devices for measuring environmental Rn concentrations based on electrostatic collection (by diffusion or pumping the gas in the detection chamber): RAD7 (Durrige Company, Billerica, MA, USA) [8]; RaMonA (University of Naples and University of Campania-Caserta, Italy) [9]; Radim3A (Tesla Company, Pruhonice, Czech Republic) [10]. Rn measurements and the development of new detection systems remain popular topics [11]. Using of tools for the simulation is a required and useful part for the complete evaluation and realization of such kind of Rn measurement systems. In particular, it is necessary to reproduce the geometry of the system, and simulate the α-spectrum of the Rn with progeny and the correspondent detection efficiency for each involved radionuclides. This leads to being able to study the efficiency versus the different geometries of the system, in order to find the best experimental solution to be realized. Different Rn measurement systems can, hence, be developed for several purposes. Tools for the simulation of the particles passage through matter are based on Monte Carlo method, the major ones are FLUKA and Geant4 [12]. Both can handle complex 3D-geometries using combinatorial geometry algorithms.
In literature, papers focused on Rn monitoring only use commercial devices for their purposes, without worrying about simulations [13,14]; while papers focused on investigation of new Rn measurement systems only deal with the realization and the testing of them without handling the simulation [15,16,17]. Within the research framework of simulation of Rn detection systems, in particular those based on electrostatic collection, very few studies are published, e.g., a low-cost monitor for the thoron [18,19] and radon [20] measurement was designed. These studies obtained results by using FLUKA or Geant4 with restrictions due to the gaps of both simulation tools. In [19] only non-ionized 220Rn atoms are simulated with FLUKA, while in [18,20] the efficiency of 220Rn and the spatial distribution of 214Po ions were determined using Geant4 with a constant electric field inside the metallic chamber. In fact, FLUKA does not allow to add the electric field, contrary to Geant4 where it is possible but constant or as a function of spatial coordinates. For these reasons, in the present work we introduce the COMSOL software, a powerful tool for the resolution of problems at partial differential equations, which simulates an electric field as gradient of the applied potential difference within metallic chamber and creates the correspondent map useful for Geant4 input [21]. The aim of our work is to provide a useful and standard methodology for the study of radon and thoron measurement systems based on electrostatic collection by using Geant4+COMSOL simulation tools. The main commercial devices will be simulated and the results compared with the experimental ones. Due to the short 220Rn half-life, in works [18,19] emerge the importance of the realization of systems for the direct detection of the α-particles produced by the non-ionized 220Rn atoms, to avoid their energy loss before they reach the detector. Therefore, a simulation study of the detection efficiency of 222−220Rn and their progeny versus different geometries of measurements system is also performed.

2. Software

Geant4 (Geometry and Tracking, 4th generation) is the used software for the simulation of the electrostatic detection of 222−220Rn and their progeny on a silicon detector inside a high voltage metallic chamber. Geant4, a publicly available Monte Carlo toolkit developed at CERN, is used to perform accurate simulations of the interaction and transport of particles through matter [22]. Its fields of application include high energy, nuclear and accelerator physics, as well as studies in medical and space science. Unlike the other computational Monte Carlo codes, such as FLUKA, Geant4 is not an executable program but rather a set of C++ class libraries for running, tracking, stepping and collecting information. The version 10.7 of Geant4 for Windows platform, released in December 2020, has been used in this work. COMSOL Multiphysics (or simply COMSOL) is used to simulate the electric field generated within the metallic chamber walls of the Rn measurement system for the electrostatic collection of the gas isotopes daughters. COMSOL is a finite element analysis and free solver software package for various physics and engineering applications [21]. It includes an environment for modeling any physical phenomenon that can be described using ordinary (ODEs) or partial differential equations (PDEs). Electrostatic module of COMSOL has been used to determine the values of the electric fields, as potential gradient, to be used in Geant4. COMSOL version 5.5 for Windows platform, released in February 2020, has been used in this work. Geant4 provides several widget toolkits for creating graphical user interfaces and visualization system of results. Qt application has been chosen being fast, smart and open-source. Moreover, a Matlab code has been developed and used to plot the spectrum data obtained with Geant4+COMSOL simulation and for the statistical analysis of the results.

3. Methods and Simulation

3.1. Decay Process and Electrostatic Collection

222Rn and 220Rn decay by emitting α-particles with energy of 5.49 MeV and 6.30 MeV respectively. 222Rn α-emitter daughters are 218Po (6.02 MeV) and 214Po (7.68 MeV); while those of 220Rn are 216Po (6.77 MeV), 212Bi (6.09 MeV) and 212Po (8.77 MeV). The α-particles in air lose some or all of their energy before to hit the detector, except those emitted nearby that can be detected with full energy [23]. It is well known that the daughter nuclei of 222Rn and 220Rn become positively charged in air: at atmospheric pressure and temperature, with a standard relative humidity of 45%, the fraction of ionized daughters of radon e thoron is 88%, the remaining 12% are neutral [24]. Therefore, it is possible to electrostatically collect the 88% part that will decay with full α-energy on detector, inside a metallic chamber. To achieve the electrostatic collection, a high potential difference is applied between chamber and detector: an electric field is formed inside the detection chamber and the positive charged ions can be collected on the detector surface grounded. Instead, 222Rn, 220Rn and about their 12% daughters are neutrally charged in air; a direct detection of them is possible and the emitted α-particles arrive with energy depending on the sizes of the chamber [18,24].
From the simulation point of view, to reproduce the physical process above described in Geant4, the pre-package physics list QGSP_BERT_HP, recommended for low energy dosimetric application, is used, which contains the process G4RadioactiveDecay to simulate radioactive decay processes. The radionuclides production is made by G4VUserPrimaryGeneratorAction class, where through the use of the GeneratePrimaries() method they are defined by specifying the atomic number, mass number, ion charge and excitation energy. Each radionuclide is simulated separately, distinguishing between ionized and non-ionized and respecting the percentage proportions of the decay chain to which it belongs. For our purposes, 222−220Rn and their progeny are generated randomly with a uniform probability density function (PDF) within the whole chamber volume. This is performed with the random number generator by G4UniformRand(). Emitted α-particles produced during decay are generated by choosing a random emission direction in the whole 4π solid angle.

3.2. Geometry of Measurement System and α-Detector

The measurement system consists of a silicon detector that is an active element and a metallic chamber, of with different shapes and dimensions, in which the detector is mounted. Since high voltage is applied to the chamber walls, the detector and the chamber are electrically isolated by means of a thin polyethylene layer wrapping the detector non active surfaces [20]. Silicon detectors are known to trap α-radiation and to measure the effect of incident charged particles. They provide high-resolution; the electrostatic transport allows high-efficiency alpha spectroscopy and very low background.
From the simulation point of view, to design the system geometry and define the materials in Geant4, the Construct() method is used by the G4RunManager class. Materials, the ‘world’ environment and all components must be defined. Firstly, the creation of the ‘world’ environment is required as physical volume that represents the experimental area and contains all other components. Inside the ‘world’, the shape of the system with the detector can be defined, by G4VSolid class. This class has several preset shapes already implemented, such as parallelepiped, hemisphere and cylinder with G4Box, G4Sphere and G4Tubs respectively. The volume inside the system, where the investigated radionuclides move, is filled by air. The G4Material Data-Base allows choosing the material to be assigned to the system components. The silicon detector uses G4_Si material; the air volume uses G4_AIR, which is surrounded by a metal grid made of G4_Fe material. The insulator, made up by G4_POLYETHYLENE, is placed between detector and metal grid. The 3D-coordinates of the system within the ‘world’ are defined by G4VPhysicalVolume. The silicon detector is made sensitive to α-radiation by SensitiveDetector+HitsCollection tools.

3.3. Calculation of the Electric Field with COMSOL

In order to electrostatically drive the 222−220Rn progeny ions on the silicon detector, a specific high voltage (V) is applied to the metallic walls of the chamber to generate an electric field ( E ) between them and the detector surface, according to the following PDEs:
E = V
· ( ε 0   ε r   E ) = ρ   [ Gauss s   law ]
where ε0 and εr are respectively the dielectric constants in vacuum and relative to the material, and ρ is the volume charge density [25]. In our case, high voltage equipment, space charges due to ions are negligible (ρ = 0) and therefore the equation to be solved for the dielectric media is:
· [   ε 0   ε r   ( V ) ] = 0 .
The program solves (3) for voltage (potential) distribution V over user defined domain with sources and boundary condition.
From the simulation point of view, in COMSOL it is necessary to reproduce the same geometry and materials of the measurement system defined in Geant4, so that the dimensions coincide. This is essential in order to avoid code errors of the interaction between the two simulation tools. The COMSOL panel control is well organized and particularly intuitive: with the ‘Model Tree’ command it is possible to edit each step of the whole simulation process, i.e., defining 3D-geometry and materials of the system, specifying the PDEs that describe the physical phenomenon, solving the PDEs and then obtaining the solutions. By using Geometry Tool of ‘Model Tree’ command it is possible to introduce several geometrical shapes (such as blocks, cones, cylinders, spheres) for the realization of the system components. The materials of the objects previously created is performed using the Materiel Library provided by the software. A defined high voltage is applied to the metallic walls of the chamber, and the optimal bias voltage (to ground) is applied to the detector (25 V). The involved PDEs, above described, for the generation of the electric field are selected by Electric Potential Tool. Then, it is possible to start the simulation by Run command. The data that the user can save and plot are stored in the Results section of the software. COMSOL simulates the flow rate of air sampling assembly, electric field map and particle trajectories in the 3D-geometry. Using the 3D Plot Group, the simulated potential distribution and the electric field spatial distribution can be obtained respectively. The electric field map, to be used in Geant4, is exported to a text-file by using the command Copy Plot Data to Clipboard. The field map contains the coordinates and the value of the electric field in the chamber volume.

3.4. Simulation and Spectrum Formation

Once the electric field map is created (Section 3.3) and fed to Geant4 simulation (Section 3.1 and Section 3.2), it is possible to launch the application by using/run/beamOn command with 108 primary starting particles source (enough to make small the error). The developed simulation is time-independent; this means: (i) Rn diffusion processes in the system are not taken into account and, hence, each radionuclide is originated in the whole volume; (ii) the equilibrium between parent and daughter radionuclides is assumed (after ~4 h for 222Rn, and after ~4 days for 220Rn). The simulation results are used to generate an energy spectrum. The G4AnalysisManager class is used to create and fill a 1D-histogram (the spectrum) of the energy deposited per event for sensitive detector volume.
These data are stored in a CSV-file for each individual radionuclide; then, they are combined together into a single file using specially developed Matlab code.

4. Results and Discussion

Worldwide-used commercial devices for Rn detection, based on electrostatic collection, have been simulated with our developed procedure by using Geant4+COMSOL. The results will be compared with the experimental ones provided in literature and/or stated by the manufacturer company. The simulated devices are:
  • RAD7 (Durrige Company, USA), likely the most known and used, consisting of a hollow hemispherical metallic chamber of 700 ± 5 cm3 volume, with a silicon surface barrier detector, having a 300 ± 1 mm2 circular active area, centered on the bottom [26].
  • RaMonA (University of Naples and University of Campania-Caserta, Italy), consisting of a 785 ± 8 cm3 hollow cylindrical metallic chamber, 10.0 ± 0.1 cm high and 10.0 ± 0.1 cm diameter, with a silicon surface barrier detector, having a 50.0 ± 0.4 mm2 circular active area, centered on the top [27].
  • Radim3A (Tesla company, Czech Republic), consisting of a hollow hemispherical metallic chamber of 830 ± 5 cm3 volume, and a silicon surface barrier detector, having a 450 ± 2 mm2 circular active area, centered 1 cm below on the bottom [10].
The sampling of the measuring air takes place by diffusion of 222−220Rn inside the chamber in Radim3A; inside RaMonA the air is pumped; while RAD7 works in both ways. The high voltage power circuit of the three systems charges the metallic chamber to a potential of 2500 V for RAD7, 2000 V for Radim3A, and 3500 V for RaMonA, creating an electric field throughout the volume that propels 222−220Rn positively charged particles onto the detector [8,9,10].
In all systems, the detector is insulated in a plastic holder. The systems have an experimentally estimated detection efficiency of 218Po (about the same assumed for 216Po [23]) of: 0.007 ± 0.001 cpm·m3·Bq−1 (0.11 ± 0.03 cps·L·Bq−1) for RAD7, 0.8 ± 0.1 count·m3·Bq−1·h−1 (0.22 ± 0.02 cps·L·Bq−1) for Radim3A, and 0.047 ± 0.005 cps·L·Bq−1 for RaMonA [26,27]. Considering the volume of the three chambers, the systems have 218−216Po detection efficiency of: 16 ± 4% for RAD7, 27 ± 3% for Radim3A, and 6 ± 1% for RaMonA.
Geant4+COMSOL simulation reproduce the geometries of the three investigated measurement systems, with the specific features of the electric field generated inside the chamber. The thicknesses of the metallic chamber, the detector insulator holder and the detector sensitive surface are set at 1 mm, 2 mm and 500 μm, respectively (not affecting the simulation). In Figure 1 are reported the 3D-model geometry of the three Rn measurement systems in Geant4, together with the correspondent simulated electric potential spatial distributions and electric field flow lines inside the metallic chamber by COMSOL.
An example of spectrum of mixed atmospheres of 222Rn + 220Rn from Geant4+COMSOL simulation in the RaMonA device is reported in Figure 2, compared with an experimental one [23,28].
The simulated efficiency, obtained by the spectrum from Geant4+COMSOL, is computed as the ratio between the alpha particles detected on the Silicon detector and the number of alpha primary starting particles generated in the detection volume of chamber. The efficiency uncertainty comes from the propagation error of both elements of the ratio [19], where the numerator error is the square root of the counts sum. The efficiencies of the three measurement systems obtained by Geant4+COMSOL are the following: 16 ± 2% for RAD7, 5 ± 2% for RaMonA, and 28 ± 3% for Radim3A. These results are perfectly consistent, considering the uncertainties, with the experimental ones provided in literature and/or stated by the manufacturer, which are 16 ± 4% for RAD7, 6 ± 1% for RaMonA, and 27 ± 3% for Radim3A.
A further study has been performed to simulate, by Geant4+COMSOL, two different versions of RaMonA system, called Ramonino and Ramonello, used in diffusion mode respectively as the reference monitor in the Rn exposure facility of Universities of Naples and Campania, and referred to the ENEA (National Metrology Institute of Ionizing Radiation Metrology in Italy) respectively [29,30,31,32,33]. Ramonino and Ramonello consist of 190 ± 4 cm3 and 300 ± 5 cm3 hollow hemispherical metallic chambers with silicon surface barrier detectors, having a 50.0 ± 0.4 mm2 and 100 ± 1 mm2 circular active area, on the bottom, respectively [28]. The basis of the hemispherical chamber of both systems has a hole for the detector (insulated in a plastic holder). A high voltage of 1500 V is applied to the metallic chamber of both systems, to electrostatically collect 222−220Rn ionized daughters on the detector. Experimentally, the measurement systems have a 218−216Po detection efficiency of: 0.0108 ± 0.0003 cps·L·Bq−1 for Ramonino, and 0.0275 ± 0.0008 cps·L·Bq−1 for Ramonello [28]. Considering the volume of the two measurement chambers, Ramonino and Ramonello have 6.0 ± 0.2% and 9.0 ± 0.3% of 218−216Po detection efficiencies. In addition, in this case, the simulated efficiencies of Ramonimo and Ramonello, i.e., 6 ± 1% and 8 ± 2%, by Geant4+COMSOL agree within uncertainty with the experimental ones. No figures of the Ramonino and Ramonello systems are shown because they are in [28] and, also, have geometries similar to those of RAD7 and Radim3A and, therefore, from the simulation point of view they would have been similar and brought only an overload of repetitive figures.
Based on all previous considerations, a new Rn monitoring device (Ramonetto) has been designed and realized to be used in diffusion mode as new improved reference monitor in the Rn exposure facility of Universities of Naples and Campania. It consists of a 1072 ± 3 cm3 hollow hemispherical metallic mesh chamber with a silicon PIN photodiode detector, having a 784 ± 2 mm2 square active area, on the bottom (insulated in a plastic holder). A voltage of 2500 V is applied to the metallic chamber to electrostatically collect 222−220Rn ionized daughters on the detector. The measurement system has a high 218−216Po detection efficiency of: 0.25 ± 0.03 cps·L·Bq−1, and 24 ± 3% considering the volume of the measurement chamber. Again, the efficiency of the new device is consistent with the simulated one of 23 ± 4% by Geant4+COMSOL. The experimental efficiency of Ramonino, Ramonello and Ramonetto were obtained using 222Rn + 220Rn mixed atmospheres realized by calibrated sources of a pylon RNC-RN-1025 226Ra for 222Rn, and 232Th salts for 220Rn [28]. Figure 3a shows Ramonetto measurement system. It is worth mentioning that to choose the right high voltage value to be applied to the chamber walls of Ramonetto, experimental and simulated tests have been carried out on the 218Po counts per seconds (cps) versus the applied voltage value. As seen in Figure 3(a.1), the value of 2500 V has been chosen since the maximum cps saturation is reached. As Ramonetto geometry shape is equal to those of RAD7/Radim3A, no figures of Geant4 and COMSOL are shown being the simulations similar.
Referring to the work [18] that presents a new reference measurement system (Sabot chamber) mainly focused on the 220Rn (and also its ionized progeny) direct detection by diffusion, small enough to reduce as much as possible the energy loss of the non-ionized 220Rn alpha particles before they reach the detector (due to 220Rn short half-life of 55.8 s), we reproduced experimentally and by simulation that chamber. The Sabot chamber is a 1 cm high metallic chamber with a surface composed of a silicon detector of diameter 4 cm: the active volume is 13.0 ± 0.1 cm3. With 1000 V potential applied to the chamber, the 220Rn efficiency value is 35 ± 2%, by using a Pylon 228Th source to create a 220Rn atmosphere [18]. We have made a similar chamber (RadonSmall, Figure 3b) using a hallow parallelepiped metallic mesh chamber, having a volume of 7840 ± 8 mm3. A silicon PIN photodiode detector, having a square active area of 784 ± 2 mm2 and insulated in a plastic holder, acts as a basis. By using 220Rn atmosphere realized with a calibrated source of 232Th salts, and by applying a voltage of 200 V on the metallic walls of the chamber, the 220Rn detection efficiency is 24 ± 2%. In addition, for RadonSmall, to choose the best voltage value to be applied to the chamber walls, experimental and simulated tests have been carried out on the 216Po counts per second (cps) as the applied voltage varied, to reach the maximum cps saturation. The 220Rn efficiency by Geant4+COMSOL simulations (Figure 3(b.1)) gives the values of 23 ± 4%, in agreement with the experimental one. It is important to point out that in Sabot chamber and in RadonSmall, the non-ionized α-particles peak of 220Rn is obtained with the procedure described in [18,19]: subtracting from the total spectrum (220Rn + 216Po + 212Bi + 212Po) the spectrum collected after eliminating the 220Rn and consisting of 212Bi + 212Po that is achieved simply by removing the 228−232Th source.
Strengthened by all previous results, a final investigation has been conducted with the aim of simulating the detection efficiency of 222−220Rn and their progeny (218−216Po) versus different sizes of the three main geometries of the electrostatic collection measurements system (i.e., parallelepiped, cylinder and hemisphere), also varying the applied voltages. The silicon PIN photodiode of 784 ± 2 mm2 square active area is fixed as detector in each simulation, always insulated in a plastic holder. Table 1 shows the obtained results. The main conclusions are: (i) detection efficiency is strongly dependent on the size of the chambers and on the voltage value applied to obtain electrostatic collection; (ii) detection efficiency decreases as chamber’s height and/or radius increases, in agreement with the theoretical result that the detection efficiency of a radioactive source decrease as the distance source-detector increases; (iii) 222−220Rn efficiency depends on the height of the chambers due to the fact that both isotopes do not ionize, contrary to 218−216Po ions that are electrostatic collected, therefore, small chambers in height ensure the collection of the α-particles of non-ionized 220−222Rn nuclei; (iv) the electric field inside the chamber is not strong enough to collect all the ions with low voltage value, while increasing the voltage the electric field intensifies and most of the ions are collected till the efficiency saturation is reached.
The overall results of the present work, prove the proper and correct working of the developed Geant4+COMSOL procedure in simulating measurement systems for the detection of 222−220Rn and their ionized progeny by electrostatic collection on a silicon detector, inside a metallic chamber.

5. Conclusions

The present paper provides, for the first time in literature, a useful tool to design and study, by simulation, the Radon and Thoron measurement systems based on the electrostatic collection of the α-particle progeny. The proposed tool is a simulation procedure that uses a Monte Carlo method (Geant4) combined with a mathematical modeling software for solving partial differential equations (COMSOL) to recreate the entire process of Rn detection. Geant4 is for the interaction and transport of α-particles through matter, while COMSOL determines the values of the electric fields applied to the collection chamber. To validate the tool, in this work we simulated worldwide commercial devices and compared our results with known experimental ones; then, we studied new measurement systems to optimize the chamber geometry and to obtain the best collection efficiency. The following achievements have been obtained: (1) the simulation tool gives results in agreement with those from the experiments, both in efficiency and in measurement spectrum; (2) based on the simulations, new devices were realized having efficiency greater than 20% and high counting rate for measuring directly Rn isotopes and their progeny; (3) a functional datasheet that compare the detection efficiency of 222−220Rn and their progeny (218−216Po), the different collection chamber geometries, and the dependence from the applied voltages.

Author Contributions

Conceptualization, F.A. and C.S.; methodology, F.A. and L.R.; software, F.A. and L.R.; validation, F.A., C.S. and V.R.; formal analysis, F.A. and L.R.; investigation, F.A., L.R., C.S., V.R. and A.D.; resources, F.A., C.S., V.R. and A.D.; data curation, F.A., L.R., C.S., V.R. and A.D.; writing—original draft preparation, F.A. and L.R.; writing—review and editing, F.A., L.R., C.S., V.R. and A.D.; visualization, F.A., C.S. and V.R.; supervision, F.A. and C.S.; project administration, F.A., C.S. and A. D.; funding acquisition, C.S., V.R. and A.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Siročić, A.P.; Stanko, D.; Sakač, N.; Dogančić, D.; Trojko, T. Short-term measurement of indoor radon concentration in northern Croatia. Appl. Sci. 2020, 10, 2341. [Google Scholar] [CrossRef] [Green Version]
  2. Ambrosino, F.; Sabbarese, C.; Giudicepietro, F.; De Cesare, W.; Pugliese, M.; Roca, V. Study of surface emissions of 220Rn (thoron) at two sites in the Campi Flegrei caldera (Italy) during volcanic unrest in the period 2011–2017. Appl. Sci. 2021, 11, 5809. [Google Scholar] [CrossRef]
  3. Woith, H. Radon earthquake precursor: A short review. Eur. Phys. J. Spec. Top. 2015, 224, 611–627. [Google Scholar] [CrossRef]
  4. Ambrosino, F.; Sabbarese, C.; Roca, V.; Giudicepietro, F.; De Cesare, W. Connection between 222Rn emission and geophysical-geochemical parameters recorded during the volcanic unrest at Campi Flegrei caldera (2011–2017). Appl. Radiat. Isot. 2020, 166, 109385. [Google Scholar] [CrossRef] [PubMed]
  5. Teterev, Y.G.; Tretyakova, S.P.; Golovchenko, A.N.; Ilić, R.; Skvarč, J. Determination of the absolute radon concentration in air using an electrostatic method and CR-39 detector. Radiat. Meas. 1995, 25, 645–646. [Google Scholar] [CrossRef]
  6. Sabbarese, C.; Ambrosino, F.; Roca, V. Analysis by scanner of tracks produced by radon alpha particles in CR-39 detectors. Radiat. Prot. Dosim. 2020, 191, 154–159. [Google Scholar] [CrossRef]
  7. Lane-Smith, D.; Sims, K.W.W. The effect of CO2 on the measurement of 220Rn and 222Rn with instruments utilising electrostatic precipitation. Acta Geophys. 2013, 61, 822–830. [Google Scholar] [CrossRef]
  8. Salih, N.F.; Jafri, Z.M.; Aswood, M.S. Measurement of radon concentration in blood and urine samples collected from female cancer patients using RAD7. J. Radiat. Res. Appl. Sci. 2016, 9, 332–336. [Google Scholar] [CrossRef] [Green Version]
  9. Ambrosino, F.; Sabbarese, C.; Roca, V.; Giudicepietro, F.; Chiodini, G. Analysis of 7-years Radon time series at Campi Flegrei area (Naples, Italy) using artificial neural network method. Appl. Radiat. Isot. 2020, 163, 109239. [Google Scholar] [CrossRef] [PubMed]
  10. Theodoulou, G.; Parpottas, Y.; Tsertos, H. Systematic grid-wise radon concentration measurements and first radon map in Cyprus. Radiat. Meas. 2012, 47, 451–460. [Google Scholar] [CrossRef] [Green Version]
  11. Sohrabi, M.; Ghahremani, M. Novel panorama megasize environmental radon monitor. Radiat. Phys. Chem. 2021, 181, 109325. [Google Scholar] [CrossRef]
  12. Song, Y.; Shin, J.; Cho, S.; Yoo, S.; Cho, I.; Kim, E.; Jung, W.; Choi, S.; Oh, K. A comparison study of the ridge filter parameter by using FLUKA and GEANT4 simulation codes. J. Korean Phys. Soc. 2015, 67, 96–102. [Google Scholar] [CrossRef]
  13. Pantelić, G.; Čeliković, I.; Živanović, M.; Vukanac, I.; Nikolić, J.K.; Cinelli, G.; Gruber, V. Qualitative overview of indoor radon surveys in Europe. J. Environ. Radioact. 2019, 204, 163–174. [Google Scholar] [CrossRef] [PubMed]
  14. Ambrosino, F.; Thinová, L.; Briestenský, M.; Šebela, S.; Sabbarese, C. Detecting time series anomalies using hybrid methods applied to Radon signals recorded in caves for possible correlation with earthquakes. Acta Geod. Geophys. 2020, 55, 405–420. [Google Scholar] [CrossRef]
  15. Nikolaev, V.A. Progress in the development of track radiometers for radon measurements. Radiochemistry 2019, 61, 396–407. [Google Scholar] [CrossRef]
  16. Martín Sánchez, A.; dela Torre Pérez, J. Estimating retrospective indoor radon concentrations with a new device. Appl. Radiat. Isot. 2012, 70, 2742–2745. [Google Scholar] [CrossRef]
  17. Alvarellos, A.; Chao, A.L.; Rabuñal, J.R.; García-Vidaurrázaga, M.D.; Pazos, A. Development of an automatic low-cost air quality control system: A radon application. Appl. Sci. 2021, 11, 2169. [Google Scholar] [CrossRef]
  18. Sabot, B.; Pierre, S.; Michielsen, N.; Bondiguel, S.; Cassette, P. A new thoron atmosphere reference measurement system. Appl. Radiat. Isot. 2016, 109, 205–209. [Google Scholar] [CrossRef]
  19. Ambrosino, F.; Sabbarese, C.; Buompane, R.; Roca, V. Development and calibration of a method for direct measurement of 220Rn (thoron) activity concentration. Appl. Radiat. Isot. 2020, 166, 109310. [Google Scholar] [CrossRef]
  20. Barlas, E.; Bayrak, A.; Emirhan, E.; Haciomeroglu, S.; Ozben, C.S. Spatial distribution of 214Po ions in the electrostatic collection. Appl. Radiat. Isot. 2013, 80, 23–26. [Google Scholar] [CrossRef]
  21. Bird, M.B.; Butler, S.L.; Hawkes, C.D.; Kotzer, T. Numerical modeling of fluid and electrical currents through geometries based on synchrotron X-ray tomographic images of reservoir rocks using Avizo and COMSOL. Comput. Geosci. 2014, 73, 6–16. [Google Scholar] [CrossRef]
  22. Agostinelli, S.; Allison, J.; Amako, K.; Apostolakis, J.; Araujo, H.; Arce, P.; Asai, M.; Axen, D.; Banerjee, S.; Barrand, G.; et al. Geant4-a simulation toolkit. Nucl. Instrum. Methods Phys. Res. A 2003, 506, 250–303. [Google Scholar] [CrossRef] [Green Version]
  23. Sabbarese, C.; Ambrosino, F.; Buompane, R.; Pugliese, M.; Roca, V. Analysis of alpha particles spectra of the Radon and Thoron progenies generated by an electrostatic collection detector using new software. Appl. Radiat. Isot. 2017, 122, 180–185. [Google Scholar] [CrossRef]
  24. Hopke, P.K. The initial atmospheric behaviour of radon decay products. J. Radioanal. Nucl. Chem. 1996, 203, 353–375. [Google Scholar] [CrossRef]
  25. Gu, S.; Ghosh, S. A finite element model for coupled 3D transient electromagnetic and structural dynamics problems. Comput. Mech. 2014, 54, 407–424. [Google Scholar] [CrossRef]
  26. Tan, Y.; Xiao, D.; Tang, Q.; Shan, J.; Zhou, Q.; Liu, X. Research on the lower detection efficiency of the RAD7 for 220Rn than for 222Rn. J. Instrum. 2014, 9, T06001. [Google Scholar] [CrossRef]
  27. Ambrosino, F.; Thinová, L.; Briestenský, M.; Giudicepietro, F.; Roca, V.; Sabbarese, C. Analysis of geophysical and meteorological parameters influencing 222Rn activity concentration in Mladeč caves (Czech Republic) and in soils of Phlegrean Fields caldera (Italy). Appl. Radiat. Isot. 2020, 160, 109140. [Google Scholar] [CrossRef]
  28. Ambrosino, F.; Buompane, R.; Pugliese, M.; Roca, V.; Sabbarese, C. RaMonA system for radon and thoron measurement. Nuovo Cimento C 2018, 41, 222. [Google Scholar]
  29. Berlier, F.; Cardellini, F.; Chiaberto, E.; Garlati, L.; Giuffrida, D.; Faure Ragani, M.; Leonardi, F.; Magnoni, M.; Minchillo, G.; Prandstatter, A.; et al. Main results of the second AIRP international radon-in-field intercomparison for passive measurement devices. Radiat. Meas. 2019, 128, 106177. [Google Scholar] [CrossRef]
  30. Sabbarese, C.; Ambrosino, F.; D’Onofrio, A.; Roca, V. Radiological characterization of natural building materials from the Campania region (Southern Italy). Constr. Build. Mater. 2021, 268, 121087. [Google Scholar] [CrossRef]
  31. Ambrosino, F.; Stellato, L.; Sabbarese, C. A case study on possible radiological contamination in the Lo Uttaro landfill site (Caserta, Italy). J. Phys. Conf. Ser. 2020, 1548, 012001. [Google Scholar] [CrossRef]
  32. Bochicchio, F.; Campos-Venuti, G.; Piermattei, S.; Nuccetelli, C.; Risica, S.; Tommasino, L.; Torri, G.; Magnoni, M.; Agnesod, G.; Sgorbati, G.; et al. Annual average and seasonal variations of residential radon concentration for all the Italian Regions. Radiat. Meas. 2005, 40, 686–694. [Google Scholar] [CrossRef]
  33. Sabbarese, C.; Ambrosino, F.; D’Onofrio, A.; Roca, V.; Pugliese, M. Natural radioactivity in soils and materials of the Campania region (Italy). Nuovo Cimento C 2020, 43, 152. [Google Scholar]
Figure 1. Schematic view of the Geant4 and COMSOL simulations of the three main worldwide commercial 222−220Rn measurement system RAD7 (a), RaMonA (b) and Radim3A (c). For Geant4 are shown the 3D geometry and simulation of the primary starting particles source; for COMSOL are shown the simulated electric potential spatial distributions and electric field flow lines.
Figure 1. Schematic view of the Geant4 and COMSOL simulations of the three main worldwide commercial 222−220Rn measurement system RAD7 (a), RaMonA (b) and Radim3A (c). For Geant4 are shown the 3D geometry and simulation of the primary starting particles source; for COMSOL are shown the simulated electric potential spatial distributions and electric field flow lines.
Applsci 12 00507 g001
Figure 2. Experimental (a) and Geant4+COMSOL simulated (b) spectrum of 222Rn (6344 ± 541 Bq) + 220Rn (3258 ± 356 Bq) mixed atmospheres in the RaMonA measurement system [23].
Figure 2. Experimental (a) and Geant4+COMSOL simulated (b) spectrum of 222Rn (6344 ± 541 Bq) + 220Rn (3258 ± 356 Bq) mixed atmospheres in the RaMonA measurement system [23].
Applsci 12 00507 g002
Figure 3. Ramonetto measurement chamber (a) and the graph of the experimental and simulated tests for choosing the best HV to apply (a.1). RadonSmall chamber, based on the prototype in [18] (b), with the correspondent Geant4+COMSOL simulations (b.1): Geant4 geometry and primary starting particles source (top), COMSOL electric potential and electric field (bottom).
Figure 3. Ramonetto measurement chamber (a) and the graph of the experimental and simulated tests for choosing the best HV to apply (a.1). RadonSmall chamber, based on the prototype in [18] (b), with the correspondent Geant4+COMSOL simulations (b.1): Geant4 geometry and primary starting particles source (top), COMSOL electric potential and electric field (bottom).
Applsci 12 00507 g003
Table 1. Efficiency results of 222−220Rn and 218−216Po from the Geant4+COMSOL simulations of measurement systems based on electrostatic collection of α-particles, by varying the size of the collection chambers and the high voltage applied to the metallic walls of the chambers.
Table 1. Efficiency results of 222−220Rn and 218−216Po from the Geant4+COMSOL simulations of measurement systems based on electrostatic collection of α-particles, by varying the size of the collection chambers and the high voltage applied to the metallic walls of the chambers.
Chamber (r = radius, h = height, s = side) [mm]High Voltage [V]
Uncertainty = ±5 V
218−216Po Efficiency [%]
Mean Uncertainty = ±2%
222−220Rn Efficiency (%)
Mean Uncertainty = ±2%
Parallelepiped—h = 10, s = 2870120200350500263031323123
Parallelepiped—h = 40, s = 28350500700900100016192324247
Hemisphere—h = 20, r = 2050100200350500182530353528
Hemisphere—h = 30, r = 30100200400600800131823272713
Cylinder—h = 10, r = 2070120200350500192629292825
Cylinder—h = 40, r = 203507008009501100791514146
High Voltage [V]
Uncertainty = ±5 V
Chamber (r = radius, h = height, s = side) [mm]218−216Po Efficiency [%]
Mean Uncertainty = ±2%
222−220Rn Efficiency (%)
Mean Uncertainty = ±2%
350Parallelepiped s = 28474326211344371496
h = 1h = 5h = 20h = 30h = 45
350Cylinder r = 2038332214634291384
h = 1h = 5h = 20h = 30h = 45
500Cylinder h = 1039322316113630242015
r = 10r = 15r = 25r = 30r = 35
500Hemisphere4237312013383124105
r = 10r = 15r = 25r = 35r = 40
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rinaldi, L.; Ambrosino, F.; Roca, V.; D’Onofrio, A.; Sabbarese, C. Study of 222−220Rn Measurement Systems Based on Electrostatic Collection by Using Geant4+COMSOL Simulation. Appl. Sci. 2022, 12, 507. https://0-doi-org.brum.beds.ac.uk/10.3390/app12010507

AMA Style

Rinaldi L, Ambrosino F, Roca V, D’Onofrio A, Sabbarese C. Study of 222−220Rn Measurement Systems Based on Electrostatic Collection by Using Geant4+COMSOL Simulation. Applied Sciences. 2022; 12(1):507. https://0-doi-org.brum.beds.ac.uk/10.3390/app12010507

Chicago/Turabian Style

Rinaldi, Luigi, Fabrizio Ambrosino, Vincenzo Roca, Antonio D’Onofrio, and Carlo Sabbarese. 2022. "Study of 222−220Rn Measurement Systems Based on Electrostatic Collection by Using Geant4+COMSOL Simulation" Applied Sciences 12, no. 1: 507. https://0-doi-org.brum.beds.ac.uk/10.3390/app12010507

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop