Monte Carlo method for determining the response of portable gamma detector for in-situ measurement of terrestrial gamma ray field

Monte Carlo method for determining the response of portable gamma detector for in-situ measurement of terrestrial gamma ray field a b corresp first-author Corresponding author: Email address: askrib@gmail.com (Boubaker Askri) Corresponding author: Email address: askrib@gmail.com (Boubaker Askri) askrib@gmail.com Department of Physics, College of Science, Qassim University, P.O.Box: 6666-Buraidah, K.S.A Department of Physics, College of Science, Qassim University, P.O.Box: 6666-Buraidah, K.S.A Research Unit of Nuclear Physics and High Energy, Faculty of Sciences of Tunis, University of Tunis El Manar 2092 Tunis, Tunisia Research Unit of Nuclear Physics and High Energy, Faculty of Sciences of Tunis, University of Tunis El Manar 2092 Tunis, Tunisia A Monte Carlo method to study the response of portable detector to terrestrial gamma ray is proposed. This method is based on two-stage Monte Carlo simulation. First, the probability distributions of the phase space coordinates of the events that are most likely to be detected are reconstructed at the phase space shell level. The phase space shell is a closed surface enclosing the detector. The detector response to events originating from the phase space shell is then studied. The full absorption spectra as well as the partial absorption spectra are obtained for natural radionuclides uniformly distributed in the ground. For validation, this method is applied to a HpGe portable detector previously studied. The previous study is based on a semi-empirical model. Good agreement is achieved when we compare the full energy peak efficiencies and the total in-situ spectra obtained by the two methods. As an application, the effective depth of the activity of the 137 Cs artificial radionuclide in the soil is determined from the low-energy part of the total in-situ spectrum. Gamma ray spectroscopy portable detector Monte Carlo Geant4 Radioactivity 1- Introduction 1 Primordial radionuclides have been present in the crust of the Earth since its formation. The natural decay-chains of the uranium and thorium series and of 40 K radionuclide are usually accompanied by the emission of gamma ray radiation gathering information about the primary source. The release of anthropogenic radionuclides from atmospheric testing of nuclear explosive devices and releases from nuclear reactors presents another source of environmental gamma ray radiation. Depending on gamma ray interaction mechanisms with different materials, many gamma radiation detectors have been conceived and studied [ 1 1 , 2 2 ]. These detectors are used for several purposes, ranging from astrophysical observations [ 3 3 ] to gamma ray field measurement in radiation processing and radioprotection [ 4 4 - 9 9 4-9 ]. In the present work, we are concerned specifically with two types of detectors; the semiconductor-based detectors and the scintillation detectors. Their principles are described in ref. [ 10 10 ] and [ 6 6 ]. An advantage of such detectors is the possibility of manufacturing light and portable prototypes. In the present work, a semi-conductor based portable detector has been studied, but the same method could be applied to a scintillation portable detector. Portable detectors are widely used for in-situ measurements to collect information about the source distribution in the ground. The response of a portable detector to incoming photons is measured as a pulse height distribution, which essentially represents the spectral distribution of energy depositions in its active volume. Full energy depositions of primary photons lead to full energy peaks. But a substantial fraction of the recorded events corresponds to photons for which only a partial deposition of their incoming energy took place, due to the escape of Compton-scattered photons and of secondary radiation (fluorescence photons, annihilation gamma-rays, bremsstrahlung and electrons) from the active volume, as well as to energy loss after scattering from the non-active parts (housing and cryostat) of the detector. In conventional applications of in-situ gamma-ray spectrometry with portable detectors, only counts of full energy peaks are determined in the measured detector response. The photopeaks coupled with careful and precise energy calibration, are used for the determination of the nuclides in a sample. The net area of the peak is directly related to the activity of each isotope after correcting for the efficiency of the detector, the branching ratio of the source and its half life. But it is of considerable interest to obtain from the detector response the full spectral distribution of the incident photon fluence, including the continuum component coming from photons that have scattered in the environment. This component contains valuable information about the source distribution [ 11 11 ]. For example, given the energy dependence of the attenuation of gamma rays in the soil and the ratio of the scattered radiation to the unscattered one, it is possible to extract information about the depth distribution of the source [ 11 11 ]. In-situ Gamma ray field measurement requires a prior knowledge of the detector response. However, the large scale geometry of the soil and air media involved and the reduced dimensions of the detector make it difficult to perform a field calibration of the portable detectors. To deal with this problem, a semi-empirical method has been widely applied in previous works [ 6 6 - 9 9 6-9 , 12 12 , 13 13 ]. In the semi-empirical method, the response of a portable gamma detector is described by a calibration factor. The calibration factor for in-situ measurements can be expanded into a term corresponding to the intrinsic efficiency, which deals with the partial absorption of the energy of the incident events, and a term corresponding to the angular response of the detector. According to the literature, these two terms are obtained experimentally in the laboratory by studying the response of the portable detector to the gamma ray flux emitted by point-like radioactive sources of known activity and known emission energy distant from the portable detector and positioned at different angles. In certain works, the same experience has been simulated by numerical techniques [ 9 9 ]. This protocol attempts to mimic the semi-infinite geometry which the detector will be exposed to in real measurement conditions. Details of this method can be found in ref. [ 10 10 ] and in ref. [ 13 13 ]. Despite being effective, the semi-empirical method does not take into account the real distribution of the gamma ray field in the air emitted by sources in the ground. This represents a source of uncertainty, especially when one tries to identify sources of very low activity. Nowadays, Monte Carlo methods are widely used to solve complex physical and mathematical problems. In Monte Carlo simulation of radiation transport, the history of a particle is viewed as a random sequence of free flights that end with an interactive event. The Monte Carlo model is viewed as a good reproduction of reality. Better description of the distribution of the gamma ray field could be achieved by applying a convenient Monte Carlo simulation. Thereby the response of a portable detector could be determined as accurately as possible, avoiding the uncertainty hat arises from the application of the semi-empirical method. It is worth to mention that in the literature, it was not possible for authors to use only Monte Carlo simulation for the determination of the entire response of portable detectors for an open field geometry such as in-situ measurements. Rather Monte Carlo technique was a support for a semi-empirical protocol of calibration. In the present work, a Monte Carlo based approach is proposed to determine the portable gamma detector response for in-situ measurement. This approach is based on the assumption that each photon track reaching the detector should originate from an event starting at a spherical surface enclosing this detector. We denote this surface by the detector phase space shell. The approach consists in two-stage Monte Carlo simulation. In the first stage, an optimized geometry of the soil-air medium is applied to determine the probability distribution of the phase space coordinates of the events that will originate from the detector phase space shell. In the second stage, events starting from the phase space shell are propagated and the detector response to those events is studied.The full energy absorption spectrum and the partial energy spectrum are obtained. This approach permits to generate events in the detector phase space shell that contribute effectively to the gamma ray field distribution around the detector volume. It also permits to reduce considerably the amount of computation time needed to calculate the portable detector response in the open soil-air geometry. For validation, this method is applied to the Canbera p type HpGe detector used in ref. [ 13 13 ]. The simulated response is compared to the one obtained in the previous study based on semi-empirical method. As an application of the method, the effective depth of 137 Cs artificial radionuclide is calculated in the same location where the experimental spectrum is measured. 2- Material and Methods 1 2-1 Method description 2 The transport of a gamma photon that starts from the source and propagates in the soil-air medium consists of a dynamic evolution in a 7-dimentional phase space where the position of the photon is described by the three coordinates x, y and z, its direction is described by the three components of its momentum direction p x , p y and p z and its energy by one variable. In this description, the detector is represented by a geometric region in space. This region serves at scoring the useful events. For the case of in-situ gamma ray field measurement, the dimensions of the detection region are typically small compared to the whole geometry of propagation. Monte Carlo methods are based in determining the detector response by simulating all relevant physical processes taking place along the path of a photon emitted by the source. The main disadvantage of this technique is that a large number of histories must be simulated to obtain a low statistical uncertainty. Therefore, Monte Carlo methods need long computing times. To circumvent this drawback, the detector response should be obtained indirectly. The basic idea behind the actual work is that every event which corresponds to a photon track reaching the detector region is considered to be starting at a closed surface surrounding the detector. Let’s call this surface the detector phase space shell. The detector response to gamma ray photons emitted by radionuclides distributed in the ground could be obtained by tracking only the photons that emerge from the shell events. Every event is equivalent to a set of seven coordinates sampled from predetermined probability distributions. It is possible to reconstruct the probability distribution of the phase space coordinates at the shell level by tracking every photon emitted by a source in the ground which has a high probability to reach the center of the detector region. This is achieved by simulating the soil-air medium by an optimized geometry allowing for the tracking of the photons most likely to be detected. The optimization of the soil-air geometry is based on the fact that the probability for a radiation to reach the detector is strongly correlated to the thickness of matter separating the emitter from the detector. The contribution of the cross-sections of the physical interactions encountered by the photons is ensured throw the linear attenuation coefficient. Since the gamma ray is strongly attenuated in the soil within a nominal depth d equivalent to few tens of centimetres of soil thickness, a probability cut-off could be defined for an emitted radiation to reach the detector. If the detector is at distance h from the ground level, the attenuation factor for a reference point source at the nominal depth d on the vertical axe containing the detector is <math id="M1"><mrow><mo> </mo><msub><mi>μ</mi><mi>s</mi></msub><mi>d</mi><mo>+</mo><msub><mi>μ</mi><mi>a</mi></msub><mi>h</mi></mrow></math> . <math id="M2"><mrow><msub><mi>μ</mi><mi>s</mi></msub></mrow></math> and <math id="M3"><mrow><msub><mi>μ</mi><mi>a</mi></msub></mrow></math> are the linear attenuation coefficients in soil and air respectively. Thus, for an emitted radiation by a point source within the earth distant of <math id="M4"><mrow><msub><mi>x</mi><mi>s</mi></msub><mo>+</mo><msub><mi>x</mi><mi>a</mi></msub><mtext> </mtext></mrow></math> from the detector the total attenuation factor verifies <math display="block" id="M5"><mrow><msub><mi>μ</mi><mi>s</mi></msub><msub><mi>x</mi><mi>s</mi></msub><mo>+</mo><msub><mi>μ</mi><mi>a</mi></msub><msub><mi>x</mi><mi>a</mi></msub><mo>≤</mo><msub><mi>μ</mi><mi>s</mi></msub><mi>d</mi><mo>+</mo><msub><mi>μ</mi><mi>a</mi></msub><mi>h</mi></mrow></math> It was demonstrated in ref. [ 8 8 ] that the set of the useful sources emitting this radiation are contained in a soil volume bordered by the surface defined by the equation <math display="block" id="M6"><mrow><mi>r</mi><mo stretchy="false">(</mo><mi>z</mi><mo stretchy="false">)</mo><mo>=</mo><msup><mrow><mrow><mo>(</mo><mrow><msup><mrow><mo stretchy="false">(</mo><mover accent="true"><mi>μ</mi><mo>^</mo></mover><mi>h</mi><mo>+</mo><mi>d</mi><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mi>h</mi><mo>−</mo><mi>z</mi></mrow><mrow><mover accent="true"><mi>μ</mi><mo>^</mo></mover><mi>h</mi><mo>−</mo><mi>z</mi></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><mo>−</mo><msup><mrow><mo stretchy="false">(</mo><mi>h</mi><mo>−</mo><mi>z</mi><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow><mo>)</mo></mrow></mrow><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></mrow></msup></mrow></math> where r represents the radius at soil depth z , <math id="M7"><mrow><mover accent="true"><mi>μ</mi><mo>^</mo></mover><mo>=</mo><mfrac bevelled="true"><mrow><msub><mi>μ</mi><mi>a</mi></msub></mrow><mrow><msub><mi>μ</mi><mi>s</mi></msub></mrow></mfrac></mrow></math> is the ratio of linear attenuation coefficients. Following the same arguments, the air is supposed to be composed of virtual sources reemitting scattered radiation to the detector centre. The emitted radiation beyond a reference point in air distant of <math id="M8"><mrow><msub><mi>x</mi><mn>0</mn></msub></mrow></math> from the detector is supposed to have a probability to reach the detector lower than a given cut-off. This cut-off corresponds to an attenuation factor <math id="M9"><mrow><msub><mi>μ</mi><mi>a</mi></msub><msub><mi>x</mi><mn>0</mn></msub></mrow></math> in the air material considered. The corresponding air geometry is a spherical volume centred on the detector and of radius <math id="M10"><mrow><mo> </mo><msub><mi>x</mi><mn>0</mn></msub></mrow></math> . For convenience we choose <math id="M11"><mrow><msub><mtext>x</mtext><mtext>0</mtext></msub><mo>=</mo><msqrt><mrow><msup><mi>r</mi><mn>2</mn></msup><mo stretchy="false">(</mo><mi>z</mi><mo>=</mo><mn>0</mn><mo stretchy="false">)</mo><mo>+</mo><msup><mi>h</mi><mn>2</mn></msup></mrow></msqrt></mrow></math> . This reference point corresponds to the furthest real source from the detector which lies on the ground. In this way, the suitable soil-air geometry depends on the nominal depth of soil as sole parameter. Detailed properties of the optimized geometry and its application to the gamma ray field measurement are found in [ 14 14 ] and [ 15 15 ]. Based on this methodology, the simulation of the response of a portable gamma detector is fulfilled throw two independent stages. In the first stage, the probability distributions of the phase space coordinates of the events at the shell level are reconstructed. In this stage, the soil-air medium is simulated by the optimized geometry. It is worth to mention that this stage is performed only once. The resulting distributions could be used to simulate the response of any kind of portable detector. In the second stage, the response of the studied portable detector to events generated according to the distributions reconstructed in the first stage is determined. A diagram explaining the methodology is shown in Fig. 1 Fig. 1 . For computer simulation, we use the Geant4 toolkit. Geant4 is a Monte Carlo transport code that is based on object-oriented technology. It was developed using the C++ computer language and provides transparency of physical implementation [ 16 16 , 17 17 , 18 18 ]. The low-energy electromagnetic physics package of Geant4 was used in this simulation. This package includes improved models for the simulation of photon, electron, hadron and ion interactions. It also includes models for the simulation of atomic relaxation and other atomic shell effects as well as unique models describing the interactions of radiation with biological systems [ 19 19 ]. The package is capable of describing interactions at energies down to 100 eV and can perform very detailed simulations of particle transport in a medium. Geant4 geometry package allows a set of geometric shapes along with some Boolean operations, such as intersection and union, to build more complex geometries. 2-2 Probability distributions of phase space coordinates of the events at the detector phase space shell level 2 The soil-air medium was simulated by an optimized geometry of nominal depth equal to 50 cm. This choice is adequate since 90% of gamma radiation with the maximal energy of 3 MeV is attenuated within this soil thickness. The air radius is thus equal to 622 m. Soil and air compositions are shown in Table 1 Table 1 . 1 1 Elements 1 1 Soil (density=1.6 g/cm 3 ) Weight fraction 1 1 Air(density=1.29 mg/cm 3 ) Weight fraction 1 1 H 1 1 2.2 1 1 - 1 1 O 1 1 57.5 1 1 23.2 1 1 Al 1 1 8.5 1 1 - 1 1 Si 1 1 26.2 1 1 - 1 1 Fe 1 1 5.6 1 1 - 1 1 N 1 1 - 1 1 75.5 1 1 Ar 1 1 - 1 1 1.3 The detector phase space shell consists of a sphere of 80 cm radius, centered on the detector. This highly symmetric closed shape would enable an unbiased simulation. Four photon-interactions processes — photoelectric absorption, Compton scattering, Rayleigh scattering and pair production — were considered. The simulation was performed for primary energies corresponding to the natural radioactive series of 238 U, 232 Th and the 40 K radionuclide. In total 49 energies of the natural radioactive spectrum were simulated. The data was taken from [ 20 20 ]. These energies correspond to the most intense radionuclides. A total number of 10 8 of photons were generated uniformly in the soil per primary energy. For each inward crossing of the shell sphere by an incoming photon, its energy, position and momentum direction were recorded. The resulting distributions of the seven coordinates were normalized to the total number of records. In this manner we obtained the probability distribution of each coordinate. Some of the normalized distributions are shown in Fig. 2 Fig. 2 . As clearly seen, the distributions of position and momentum coordinates in the XY plane are almost flat. This indicates the rotational symmetry of the problem around the Z axis. It is worth to mention that for a soil volume simulated using the optimized volume; we achieved a variance reduction of two orders of magnitude compared with classical Monte Carlo simulation with a considerable reduction of computation time. An average time per event of 4×10 -5 s on a Core Duo processor was recorded. 3- Simulation of the response of p-type coaxial High Purity Ge (Li) portable detector 1 We simulated the response of a Canberra p-type coaxial High purity Ge (Li) detector (HpGe) of Φ 58 mm×56 mm with a relative photo-peak efficiency of 30% at 1332 keV. This is the same detector used in [ 13 13 ]. Its layout as supplied by the manufacturer is shown in Fig. 3 Fig. 3 . It has been shown in [ 13 13 ] that the detector dead layer value should be 1.84 mm instead of 0.84 mm provided by the manufacturer. We toke this detector dead layer value throughout the simulation. 10 8 of events were generated at the shell level per energy of the natural radioactive spectrum. Each event corresponds to a photon with sampled energy, direction and position from the relevant distributions already prepared in the first stage. In GEANT4, the following electro-magnetic processes were included: Compton scattering, photo-electric effect, Rayleigh effect, pair production, multiple scattering, fluorescence and Auger effect, bremsstrahlung and ionization. The full energy peak efficiency is defined as the ratio of the number of recorded events corresponding to full absorption in the detector crystal of incident primary photons having the energy of the source, to the total number of the incident primary photons on the detector. In gamma ray spectroscopy, the full energy peak efficiency allows transforming total peak counts to concentration of the correspondent radionuclide in the ground. The detector full energy peak efficiency is function of the primary energy. Fig. 4 Fig. 4 shows the full energy peak efficiency for different radionuclides as deduced by our method, compared to the efficiency curve deduced from [ 13 13 ]. The values obtained by our method are in agreement with those calculated in [ 13 13 ]. The observed discrepancy at very low energy could be explained by the angular response of the detector considered to be uniform in the previous work, a hypothesis which is not true for the low energy gamma photons which experience severe scattering in the air medium. In our method, by generating events at the shell level according to the obtained distributions in the first sta ge, we take into account the real distribution of the gamma ray field in the detector proximity. Furthermore, in ref. [ 13 13 ], the efficiency curve has been extrapolated to higher energies. This could be the source of the observed discrepancy between the two methods at energy values higher than 2000 keV. For the energies between 600 keV and 2000 keV, efficiencies calculated by Monte Carlo method exceed the ones obtained from [ 13 13 ] by less than 5%. It is worth to mention that in ref. [ 13 13 ], the efficiencies were calculated with 5 % accuracy. We conclude that the observed discrepancy is not relevant in the energy range of interest. In Fig. 5 Fig. 5 , the spectra of partial absorption and full absorption are shown for different radionuclides of the natural decay series. It can be noted that a large amount of the partial absorption events occurred at the low energy part of the spectra. We also notice the appearance of some secondary peaks at this low part. These peaks are due to the escape and to the creation of X-ray and Auger electrons in the detector housing materials. Since the most important radioactive peaks are located above 20 keV, a cut-off applied at this energy permits to cancel the underlying undesirable background. For radionuclides with higher energies, the partial absorption events became more and more frequent due to the escape of high energy photons. This photon escape is due to the small volume of the Germanium crystal. Thus, the detector used in this work is more sensitive to radionuclides emitting low energy gamma rays. In a real detector, if the incident gamma energy E γ is above 1022 keV, pair production events result in the production of two 511 keV annihilation gamma-rays. If only one of these gamma-rays escapes while the other is completely absorbed in the detector, a separate peak in the spectrum appears at the energy E γ -511 keV. This peak is called the single escape peak. If both annihilation gamma-rays escape this gives rise to the double escape peak at E γ -1022 keV. The double escape peaks are shown in Fig. 5 Fig. 5 for 40 K and 208 Tl radionuclide's spectra. These peaks are very pronounced because pair production is the most relevant interaction of high energy gamma photons with the detector crystal. Furthermore, in the 208 Tl spectrum, there is a second peak composed of full absorption events of 511 keV annihilation photons. The intensity of the single escape peak and double escape peak depends usually on the size of the Ge crystal. For a small detector compared with the mean free path of the secondary gamma radiation, all annihilation photons escape from the detector. The spectrum will show a double escape peak but no single escape peak. This could explain the absence of single escape peak in 40 K and 208 Tl radionuclide's spectra. The pronounced double escape peaks of high energy spectra represent a significant background proving the fact that the actual detector is suitable for identifying low energy gamma sources. Fig. 6 Fig. 6 shows the gamma ray simulated spectrum for the natural decay chain series after removing the partial absorption spectrum. This spectrum is normalized to its experimentally observed full-energy peak. In the same figure, the gamma ray experimental spectrum recorded in-situ is shown. A matrix stripping method explained in detail in [ 13 13 ] was applied to the experimental spectrum to remove the partial absorption events. The shapes of the two spectra are close. The observed discrepancy could be attributable to the accuracy of the stripping method described in [ 13 13 ]. The discrepancy could also be attributable to the background component due to the electronics and contamination of the detector material which is not taken into account in our model. This background component is dominant at the high energy part of the spectrum where the detector sensitivity to the emitted energies is degraded. As a result, the discrepancy at this part of the spectrum is the most relevant. The discrepancy could also be attributable to how accurate is the assumption that the different radionuclides are homogeneously distributed in the soil and possibly to the different chemical composition and moisture of the soil where the measure was taken from the data used in the simulation ( Table 1 Table 1 ). 4- Application 1 As a consequence of nuclear probes and nuclear accidents a certain amount of the fission product 137 Cs was ejected into the atmosphere, from where it was deposited on the soil. 137 Cs is a gamma emitter with energy of 661.6 keV and with half-life of 30.17 years. Several models have been developed to describe 137 Cs migration in soils and to explain its vertical distribution [ 22 22 , 23 23 , 24 24 , 25 25 ]. After the nuclear tests in the sixties and the Chernobyl accident it was initially possible to describe the depth profile of the 137 Cs isotope in the soil by an exponential distribution. It was also shown [ 25 25 ] that for any exponential distribution of 137 Cs in the soil, an equivalent planar distribution of 137 Cs radionuclides can be found at an effective depth, which yields the same photon fluence energy and angular distribution as the original one. In the present work, the effective depth was determined in tow independent steps. In the first step, we applied a cut-off at the bin energy containing the 137 Cs peak to the experimental and to the simulated in-situ spectra shown in Fig.6 Fig.6 . We conserved only the low-energy part of the two spectra. It is worth to mention that the simulated spectrum is due to natural radionuclides uniformly distributed in the ground. The simulated spectrum was then stripped from the experimental one, leaving one with the residual spectrum belonging to the 137 Cs isotope only. This spectrum was utilized to obtain information on the 137 Cs depth distribution. In the second step, we assumed a plane distribution of 137 Cs activity at a certain depth beneath the surface. For each location of the plane source, a Monte Carlo simulation was conducted. According to the optimization applied to the soil geometry, any plane source is reduced to a disc of radius r(z) given by Eq. (2). The effective depth was determined by comparing the simulated 137 Cs spectrum due to the plane source and normalized to the experimental full-energy peak value, with the residual experimental spectrum. The location of the 137 Cs plane beneath the soil surface is varied until a good match between the two spectra is obtained. The optimal effective depth was found equal to 40 cm as inferred from Fig. 7 Fig. 7 . For an effective depth of 40 cm, the planar source is reduced to a disc of only 145 cm radius. It is worth to mention that the experimental in-situ measurement was taken in the soil region located in the vicinity of the site studied in ref. [ 26 26 ]. In this experimental study, the authors found a maximal activity of 137 Cs located between 30 cm and 40 cm which is in accordance with our prediction. 5- Conclusion 1 A Monte Carlo based approach to study the response of a portable gamma detector to terrestrial gamma ray was proposed. This approach consists in two-stage Monte Carlo simulation. In the first stage, the probability distributions of the phase space coordinates at the detector phase space shell were reconstructed for natural radionuclides uniformly distributed in the ground. For this purpose, an optimized geometry of the soil-air medium was implemented to ensure the tracking of only the events most probable to be detected. These distributions serve at studying the response of any type of portable detector for in-situ measurement. The response to events originating at the phase space shell level according to the already prepared distributions of a p type HpGe detector was determined. This method allowed to separate the background signal, consisting of the partial absorption energy spectrum, from the desired signal, consisting of the full absorption energy spectrum, for any natural radionuclide uniformly distributed in the ground. Comparison with a previous study which uses a model based on semi-empirical method has shown a good agreement for the full energy peak efficiency curve and the total in-situ spectrum. As an application of this method, the effective depth of 137 Cs artificial radionuclide in the ground at the location where the in-situ spectrum was carried out was determined with good precision. This approach is among the first attempts to determine the entire response of portable HpGe detectors for in-situ measurement of gamma ray field by Monte Carlo simulation. We believe that we could improve this approach by taking into account other parameters such as soil and air real properties, the detector resolution [ 27 27 ] and the detector electronic background to enhance its predictive power. References NIU Li-Bo , LI Yu-Lan , ZHANG Lan et al . Performance simulation and structure design of Binode CdZnTe gamma-ray detector . Nucl Sci Tech , 2014 , 25 ( 1 ): 010406 . DOI: 10.13538/j.1001-8042/nst.25.010406 http://doi.org/10.13538/j.1001-8042/nst.25.010406 Performance simulation and structure design of Binode CdZnTe gamma-ray detector XU Wenzhen , LIU Yanfen , TAN Zongquan et al . Geant4 simulation of plastic scintillators for a prototype μSR spectrometer . Nucl Sci Tech , 2013 , 24 ( 4 ): 040401 . DOI: 10.13538/j.1001-8042/nst.2013.04.011 http://doi.org/10.13538/j.1001-8042/nst.2013.04.011 Geant4 simulation of plastic scintillators for a prototype μSR spectrometer AHARONIAN F , BUCKLEY, J , KIFUNE, T et al . High energy astrophysics with ground-based gamma ray detectors . Reports on Progress in Physics , 2008 , 71 : 096901 . High energy astrophysics with ground-based gamma ray detectors Mejri A , Farah K , Eleuch H et al . Application of commercial glass in gamma radiation processing . Radiation Measurements , 2008 , 43 : 1372 - 1376 . DOI: 10.1016/j.radmeas.2008.05.010 http://doi.org/10.1016/j.radmeas.2008.05.010 Application of commercial glass in gamma radiation processing Mejri A , Farah K , Hosni F et al . An Empirical Model for Predicting Harwell Dose Variation During Gamma Process Interruption Using Experimental Design . M_APAN-Journal of Metrology Society of India , 2015 , 30 : 85 - 90 . DOI: 10.1007/s12647-014-0129-6 http://doi.org/10.1007/s12647-014-0129-6 An Empirical Model for Predicting Harwell Dose Variation During Gamma Process Interruption Using Experimental Design Glenn F K . Radiation Detection and Measurement . New York : John Wiley and Sons , 1989 . Radiation Detection and Measurement Boson J , Plamboeck A , Rameba H et al . Evaluation of Monte Carlo-based calibrations of HPGe detectors for in situ gamma-ray spectrometry . Journal of Environmental Radioactivity , 2009 , 100 : 935 - 940 DOI: 10.1016/j.jenvrad.2009.06.006 http://doi.org/10.1016/j.jenvrad.2009.06.006 Evaluation of Monte Carlo-based calibrations of HPGe detectors for in situ gamma-ray spectrometry Fehrenbacher G , Meckbach R , Jacob P . Unfolding the response of a Ge detector used for in-situ gamma-ray spectrometry . Nucl Instr Meth A , 1996 , 383 : 454 - 462 . DOI: 10.1016/S0168-9002(96)00824-8 http://doi.org/10.1016/S0168-9002(96)00824-8 Unfolding the response of a Ge detector used for in-situ gamma-ray spectrometry Hurtado S , Garcia-Leon M , Garcia-Tenorio R . Monte Carlo simulation of the response of a germanium detector for low-level spectrometry measurements using GEANT4 . Applied Radiation and Isotopes , 2004 , 61 : 139 - 143 . DOI: 10.1016/j.apradiso.2004.03.035 http://doi.org/10.1016/j.apradiso.2004.03.035 Monte Carlo simulation of the response of a germanium detector for low-level spectrometry measurements using GEANT4 Miller K M , Beck H L . Indoor Gamma and Cosmic Ray Exposure Rate Measurements Using a GE Spectrometer and Pressurised Ionisation Chamber . Radiat Protect Dosim , 1984 , 7 : 185 - 189 . Indoor Gamma and Cosmic Ray Exposure Rate Measurements Using a GE Spectrometer and Pressurised Ionisation Chamber Korun M , Martin R , Pucelj B . In-situ measurements of the radioactive fallout deposit . Nucl Instr Meth A , 1991 , 300 : 611 - 615 . DOI: 10.1016/0168-9002(91)90398-A http://doi.org/10.1016/0168-9002(91)90398-A In-situ measurements of the radioactive fallout deposit Beck H L , De Campo J , Gogolak C . In-situ Ge(Li) and Na(Tl) gamma-ray Spectrometry . New York : US Department of Energy . 258 , 1972 . In-situ Ge(Li) and Na(Tl) gamma-ray Spectrometry Askri B , Trabelsi A , Baccari B . Method for converting in-situ gamma ray spectra of a portable Ge detector to an incident photon flux energy distribution based on Monte Carlo simulation . Nucl Sci Tech , 2008 , 19 : 358 - 364 . DOI: 10.1016/S1001-8042(09)60019-1 http://doi.org/10.1016/S1001-8042(09)60019-1 Method for converting in-situ gamma ray spectra of a portable Ge detector to an incident photon flux energy distribution based on Monte Carlo simulation Askri B , Manai K , Trabelsi A et al . Optimised geometry to calculate dose rate conversion coefficient for external exposure to photons . Radiat. Prot. Dosim , 2008 , 128 : 279 - 288 DOI: 10.1093/rpd/ncm426 http://doi.org/10.1093/rpd/ncm426 Optimised geometry to calculate dose rate conversion coefficient for external exposure to photons Askri B . Application of optimized geometry for the Monte Carlo simulation of a gamma-ray field in air created by sources distributed in the ground . Radiation Measurements , 2015 , 72 : 1 - 11 . DOI: 10.1016/j.radmeas.2014.11.006 http://doi.org/10.1016/j.radmeas.2014.11.006 Application of optimized geometry for the Monte Carlo simulation of a gamma-ray field in air created by sources distributed in the ground Agostinelli S , Allison J , Amako K et al . Geant4—a simulation toolkit . Nucl Instr Meth A , 2003 , 506 : 250 - 303 . Geant4—a simulation toolkit XU Qi , YU Gang-Lin , WANG Kan et al . Research on GPU-accelerated algorithm in 3D finite difference neutron diffusion calculation method . Nucl Sci Tech , 2014 , 25 ( 1 ): 010501 . DOI: 10.13538/j.1001-8042/nst.25.010501 http://doi.org/10.13538/j.1001-8042/nst.25.010501 Research on GPU-accelerated algorithm in 3D finite difference neutron diffusion calculation method Askri B . Optimization of a photoneutron source based on 10 MeV electron beam using Geant4 Monte Carlo code . Nucl Instr Meth B , 2015 , 360 : 1 - 8 DOI: 10.1016/j.nimb.2015.07.132 http://doi.org/10.1016/j.nimb.2015.07.132 Optimization of a photoneutron source based on 10 MeV electron beam using Geant4 Monte Carlo code Chauvie S et al . Nuclear Science Symposium Conference Record . vol. 3 . IEEE 1881 , 2004 Nuclear Science Symposium Conference Record Chu, S. , Ekstrm, L. , Firestone, R. , 1999 . The lund/lbnl nuclear data search . http://nucleardata.nuclear.lu.se/nucleardata/toi Version 2.0 http://nucleardata.nuclear.lu.se/nucleardata/toi Version 2.0 The lund/lbnl nuclear data search ROOT: An object-oriented data analysis framework . http://root.cern.ch/ http://root.cern.ch/ ROOT: An object-oriented data analysis framework Szerbin P , Koblinger-Bokori E , Koblinger L et al . Caesium-137 migration in Hungarian soils . Sci Total Environ , 1999 , 227 : 215 - 227 . DOI: 10.1016/S0048-9697(99)00017-0 http://doi.org/10.1016/S0048-9697(99)00017-0 Caesium-137 migration in Hungarian soils Bunzl, K , Migration of fallout-radionuclides in the soil: effect of non-uniformity of the sorption properties on the activity-depth profiles . Radiat Environ Biophys , 2001 , 40 : 237 - 241 . DOI: 10.1007/s004110100102 http://doi.org/10.1007/s004110100102 Migration of fallout-radionuclides in the soil: effect of non-uniformity of the sorption properties on the activity-depth profiles Likar A , Vidmar T , Lipoglavsek M et al . Monte Carlo calculation of entire in situ gamma-ray spectra . J Env Rad , 2004 , 72 : 163 - 168 . DOI: 10.1016/S0265-931X(03)00198-X http://doi.org/10.1016/S0265-931X(03)00198-X Monte Carlo calculation of entire in situ gamma-ray spectra Likar A , Omahen G , Vidmar T et al . Method to determine the depth of Cs-137 in soil from in-situ gamma-ray spectrometry . J Phys D: Appl Phys , 2000 , 33 : 2825 - 2830 . DOI: 10.1088/0022-3727/33/21/323 http://doi.org/10.1088/0022-3727/33/21/323 Method to determine the depth of Cs-137 in soil from in-situ gamma-ray spectrometry Trabelsi Y , Gharbi F , El-Ghali A et al . Recent sedimentation rates in Garaet El Ichkeul Lake, NW Tunisia, as affected by the construction of dams and a regulatory sluice . J Soils Sediments , 2012 , 12 : 784 - 796 . DOI: 10.1007/s11368-012-0496-y http://doi.org/10.1007/s11368-012-0496-y Recent sedimentation rates in Garaet El Ichkeul Lake, NW Tunisia, as affected by the construction of dams and a regulatory sluice WANG Yiming , WEI Yixiang . Implicit FWHM calibration for gamma-ray spectra . Nucl Sci Tech , 2013 , 24 ( 2 ): 020403 . DOI: 10.13538/j.1001-8042/nst.2013.02.005 http://doi.org/10.13538/j.1001-8042/nst.2013.02.005 Implicit FWHM calibration for gamma-ray spectra 10.1007/s41365-016-0089-3 2015-08-22 2015-10-13 2015-10-22 2016-07-08 2016-08-20 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2016 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-016-0089-3 http://dx.doi.org/10.1007/s41365-016-0089-3 . 2016 Nuclear Science and Techniques 04 27