Nuclear resonance fluorescence drug inspection

There is an increasing challenge to prevent illicit drug smuggling across borders and seaports. However, the existing techniques in-and-of-themselves are not sufficient to identify the illicit drugs rapidly and accurately. In the present study, combining nuclear resonance fluorescence (NRF) spectroscopy and the element (or isotope) ratio approach, we present a novel inspection method that can simultaneously reveal the elemental (or isotopic) composition of the illicit drugs, such as widely abused methamphetamine, cocaine, heroin, ketamine and morphine. In the NRF spectroscopy, the nuclei are excited by the induced photon beam, and measurement of the characteristic energies of the emitted γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document} rays from the distinct energy levels in the excited nuclei provides “fingerprints” of the interested elements in the illicit drugs. The element ratio approach is further used to identify drug elemental composition in principle. Monte Carlo simulations show that four NRF peaks from the nuclei 12\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{12}$$\end{document}C, 14\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{14}$$\end{document}N and 16\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{16}$$\end{document}O can be detected with high significance of 7−24σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} using an induced photon beam flux of 1011\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{11}$$\end{document}. The ratio of 14N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{14}N$$\end{document}/12C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{12}C$$\end{document} and/or 16O\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{16}O$$\end{document}/12C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{12}C$$\end{document} for illicit drugs inspected are then extracted using the element ratio approach. It is found that the present results of simulations are in good agreement with the theoretical calculations. The feasibility to detect the illicit drugs, inside the 15-mm-thick iron shielding, or surrounded by thin benign materials, is also discussed. It is indicated that, using the state-of-the-art γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document}-ray source of high intensity and energy-tunability, the proposed method has a great potential for identifying drugs and explosives in a realistic measurement time.

The smuggling of contraband is one of the biggest threats to public health, human welfare and national security. Currently, illicit drug trafficking is becoming more and more rampant. For instance, the global cocaine production and seizures reach record highs in 2017, reaching 1976 tons and 1275 tons, respectively 1 . There are large quantities of unseized drugs still passing readily across national borders, which presents a continual challenge to the law enforcement officials. The fight against drug smuggling at ports of entry still relies much on experienced narcotic agents and well trained sniffer dogs. So far, Raman spectroscopy 2,3 is widely used on the forensic analysis of drugs. However, this technique is only limited to identify drug trafficking patterns and distribution networks, thus leaving behind the problem of large quantities of unseized drugs. It seems impossible to identify smuggled drugs by checking on every single parcel, air courier and maritime container. An efficient nondestructive inspection method for concealed illicit drugs is indispensable for reducing this combined danger.
Nuclear resonance fluorescence 4 (NRF) describes the X(γ,γ')X reaction in which a photon γ is resonantly absorbed by the nucleus X and then re-emitted as the excited nucleus subsequently transits to its lower states, particularly the ground state. The energy of NRF γ rays depends merely on the distribution of nuclear levels, and the spectrometry of NRF γ rays can then be used for nondestructive detection of heavily shielded materials, considering deep penetration of energetic γ rays. This nondestructive detection could be beneficial in industrial applications [5][6][7] as well as nuclear safeguards 8,9 . Bremsstrahlung photons 10 and neutrons 11,12 have been studied as probes for measuring composition of materials that are concealed by heavy shields. However, bremsstrahlung photons generally trigger the NRF processes in low-lying excited states (within 3 MeV), because of decreasing intensity with energy and high dose to the shielded objects, and the exact determination of neutron energy is difficult because of the multiple neutron scattering inside the objects. Recently, laser-Compton scattering (LCS) γ rays have been proposed for nondestructive detection of special nuclear materials (SNM) 13,14 (such as 208 Pb and 238 U) and chemical compounds 15,16 . Moreover, the LCS γ-ray source can generate energy-tunable, www.nature.com/scientificreports/ quasi-monochromatic γ rays 17,18 and have been used for studying nuclear physics 19,20 and nuclear astrophysics 21 , industrial and medical applications [22][23][24][25] .
In this study, we investigate through Monte Carlo (MC) simulations nondestructive inspection of illicit drugs by exploiting NRF interaction in light elements with an intense LCS γ-ray source (see Fig. 1). It uses isotopespecific data signatures, which can be used to obtain a unique fingerprint of the object, to identify elemental components ( 12 C, 14 N and 16 O) of illicit drugs. Unlike fast neutrons 11 , which does not have this isotope-specificity for both low and high Z nuclei, the backscatter NRF signals can be made highly specific and sensitive to the presence and abundance of individual isotopes. Also unlike bremsstrahlung radiations, the LCS γ rays can reduce the background flux at the low-energy side produced from atomic processes in the object 26 .
There are five widely abused drugs to be considered, including cocaine, heroin, ketamine, methamphetamine and morphine. To identify these abused drugs by their NRF signatures, we propose a novel approach of element ratio, 14 N/ 12 C and 16 O/ 12 C . Results show that four NRF signals from drug components can be readily differentiated by spectrometry of NRF γ rays, even in the presence of metal shielding. The significance of the resulting NRF peaks reaches 7−24σ with LCS γ-ray flux of 10 11 . The element ratios of 14 N/ 12 C and/or 16 O/ 12 C simulated for five abused drugs are in good agreement with the theoretical values. It hence indicates that this combined detection method can identify the illicit drugs with high confidence in a practical measurement time.

NRF yields in backscatter inspection.
In this section, we use NRF cross section, σ NRF (E) , and angular distribution, W(θ) (see Methods), to construct an semi-analytical expression for the expected NRF counts. This expression can be further used to validate the simulation algorithm (see Methods). For a collimated LCS photon beam of incident flux, I(E) interacting with the target, a small part of photon flux near the resonant energy E NRF will undergo resonant (NRF) and non-resonant (atomic) interactions inside the target. The resulting NRF yield then produces a double-differential rate of NRF interactions in the infinitesimal solid angle d 27 , where E is the energy of primary LCS photons. E ′ = E − E rec is the energy of photons induced by NRF. θ is the emission angle of the NRF photons relative to the beam propagation, L is the thickness of the irradiated target, ǫ(E ′ ) is the intrinsic detection efficiency for photopeak, µ NRF (E) = Nσ NRF (E) denotes the linear attenuation coefficient with N being the target number density, and µ eff (E, E ′ , θ) is the effective attenuation coefficient, which is given by Here µ nr (E) and µ nr (E ′ ) are the non-resonant attenuation coefficient for the primary LCS photons and secondary photons, respectively. If the energy of the γ-ray beam, which is produced by Compton backscattering of a laser pulse from high-energy electron beam, is identical with the excitation energy of the nucleus of interest, it can be effectively absorbed in the nucleus, which de-excites subsequently by emission of NRF photons. (b) Schematic view of the NRF-based inspection system using a LCS γ-ray source. After collimation, the LCS γ-ray beam penetrates through the iron shield and then impinges on the drug target, producing NRF γ rays. In order to realize backscatter NRF inspection, four high-purity germanium (HPGe) detectors are assembled on a detection plane with large angle of 135 • with respect to the beam propagation. The lead plates are installed in front of HPGe detectors to improve the ability to operate detectors by reducing the intensity of the low-energy portion of the scattered beam. www.nature.com/scientificreports/ In this study, the backscatter inspection setup is employed (see Fig. 1). Rather than forward-scatter ones, the backscatter setup can be used to suppress significantly exponential decaying backgrounds induced by nonresonant processes, and hence to improve the signal to noise (S/N) ratio 28 . We run simulations on NRF γ-ray emission and detection of 238 U at energies below 2.5 MeV. The irradiation of 238 U target is considered as a representative scenario since the isotope 238 U has several visible NRF peaks which have been observed in previous measurements 29,30 . The photon fluence I(E) in each simulation is composed of 10 11 photons with cut-off energy of 2.5 MeV, and is directed orthogonally at an isotopically-pure cylinder target without iron shielding. Note that such flux can be readily delivered by the state-of-the-art LCS γ-ray source within a few of seconds 31,32 . Figure 2a shows partial energy spectrum recorded by four HPGe detectors considering an achievable energy resolution of 0.05%. Seven NRF peaks are seen for 238 U, from which these peak counts were extracted accordingly.

NRF signatures of illicit drugs.
Similarly, NRF signatures are obtained for five drug materials, cocaine, morphine, methamphetamine, ketamine and heroin. In the simulations, the cut-off energy of LCS γ-ray beam is shifted from 2.5 MeV to ∼ 7.2 MeV, which is obtained by the Compton backscattering of a 532 nm laser pulse from a 452 MeV electron beam with an energy spread of 0.05%. Regarding the MC simulation models and conditions, more detailed information is given in section of Methods. Figure 3 shows the exemplary energy spectrum for cocaine material, which is merged with four separate HPGe detectors to improve the S/N ratio. One can see four NRF signals can be clearly detected, where 12 C contributes the 4438 keV (0→2→ 0) peak, 14 N the 7029 keV (1→2→ 1) peak, and 16 O the 6917 keV (0→2→ 0) and 7116 keV (0→1→ 0) peaks. Three NRF signal regions near 4438 keV 12 C, 7029 keV 14 N and 7116 keV 16 O are then fitted separately with Gaussian function, on top of a quasi-linear decaying continuum background. The fitting function for each of these NRF peaks is written as where c 1 and c 2 describe the shape of the background, and a, E NRF and σ E are the area, mean and standard deviation (SD) fit parameters of the peak. The fitting function can describe well the simulation curve (Fig. 3).
Once the 5-parameter fit (and associated parameter uncertainties) for each NRF peak is computed using Eq. (3), the detected NRF rate in counts were extracted as simply Y = a/�E , where the division by the spectrum bin width E enforces proper dimensions and normalization 9 . The peak uncertainty is calculated by δY = √ Y + 2B with B being the background yields. The NRF peak and background yields, and associated peak significance for five illicit drugs are summarized in Supplementary Section S1. Although the integration cross section is as low as 6.32 eV · b, the NRF signal at 4438 keV for 12 C can be detected with high confidence of 16−20σ because of high mass fraction. For cocaine, heroin and morphine, the NRF lines at 6917 and 7116 keV of 16 O exhibit comparable confidence (13−24σ ) to the 12 C. The NRF peak at 7029 keV for 14 N has less production yield and then relatively low confidence ( ∼ 7 σ ). This is mainly attributed to very small mass fraction while the integration cross section for the 14 N is almost two times higher than the 12 C. Moreover, the background level is almost one order of magnitude lower than the NRF signal. This is due to the non-resonant effect contributes visibly to the low energy region of the spectrum, given by that the LCS photon beam has sufficient intensity at the cutoff energy (7200 keV), which matches well to the 7116 keV 16 Counts per 2 keV  1)), while the points correspond to NRF yields extracted from the energy spectrum.
Here error bars denote ±1σ statistical uncertainties. www.nature.com/scientificreports/ Element ratio of illicit drugs. To further distinguish the types of illicit drugs, two-dimensional element ratio approach is proposed. After extracting the NRF peak yields, the element ratio 14 N/ 12 C can be obtained by where the Y14 N /Y12 C is the ratio between the 14 N and 12 C peak yields (see Supplementary Section S1). The intensity ratio I(E12 C )/I(E14 N ) and efficiency ratio ǫ(E12 C )/ǫ(E14 N ) were simulated with the code MCLCSS 33 , together with the Geant4 toolkit 34 . The ratios σ12 C (E)dE/ σ14 N (E)dE and W12 C (θ)/W14 N (θ) are given by Eqs. (5) and (8), respectively. According to the rules of error propagation, the uncertainties associated with the predicted element ratios can be obtained by 14 Fig. 4, together with their theoretical ratios as comparison. According to Eq. (4), the relative deviations of both the ratios 16 O/ 12 C and 14 N/ 12 C were calculated to be 0.4−11.6 % . One can see that ketamine and morphine can be well distinguished from three other drugs when considering either the ratio 16 O/ 12 C or the 14 N/ 12 C . Methamphetamine can be identified though the ratio of 14 N/ 12 C with the NRF peak at 7029 keV, since the ratio value is visibly higher than others. However, the ratio 16 O/ 12 C is not applicable for methamphetamine considering it has no oxygen content. For cocaine and heroin, they are almost identical regarding the ratio 16 O/ 12 C and have insignificant discrepancy (0.059:0.048) in the ratio 14 N/ 12 C , which suggests that accurate discrimination of cocaine and heroin is challenging when the statistic is not sufficient.
It is also necessary to check the element ratio for benign materials that are not dissimilar to the various narcotics. Caffeine and related chemicals in coffee can be regarded as good candidates of benign materials since Photon energy (MeV) www.nature.com/scientificreports/ coffee is generally used by smugglers to overwhelm drug-smelling dogs 35 . We then performed such simulation by replacing the drug material inspected with caffeine (C 8 H 10 N 4 O 2 ). The element ratios, 16 O/ 12 C and 14 N/ 12 C , are calculated accordingly and compared to those of the five illicit drugs (see Fig. 4). The comparison shows the inspected drugs have significant difference in element ratios relative to the caffeine.

Inspection of concealed cocaine.
We have presented inspection of illicit drugs without any blocking. A more realistic scenario is to smuggle drugs in heavy concealment. Thus, it is of great importance to discuss the inspection feasibility on concealed drugs (see Fig. 1). Here we used 15-mm-thick iron to enclose drugs to be inspected such as cocaine. Such shielding reduces the photon intensity from 1.0 (relative) to 0.75 at ∼ 7 MeV, so it is for the NRF γ rays emitted. Meanwhile, due to the γ-ray attenuation, a large number of backgrounds were produced. Figure 5 shows that the background photons produced in total are almost two orders of magnitude higher than without iron shielding (see Fig. 3). For cocaine, the peak significance at 4438, 6917, 7029 and 7116 keV are reduced from 16.7, 23.4, 7.7 and 15.6σ to 3.5, 13.9, 3.6 and 9.1σ , respectively. In order to relieve the background effect on the signal detection, we introduced a 1-mm-thick lead filter in front of each HPGe detector. This reduces significantly the background for peak fitting at 4438 keV such that the NRF peak at 4438 keV 12 C can be observed more clearly (see Fig. 5). This also provides a potential approach to enhance the S/N in the presence of metal shielding. In order to demonstrate the feasibility to inspect drugs even when the illicit drugs are surrounded by benign materials containing C, H, O and N, we simulated the inspection of cocaine concealed in a wrapping material consist of caffeine (C 8 H 10 N 4 O 2 ) with different areal densities. According to the simulated spectra, the NRF peak yields are extracted, and then the element ratios, ( 16 O/ 12 C) 6917 , ( 16 O/ 12 C) 7116 and ( 14 N/ 12 C) 7029 , and their associated uncertainties are obtained with the same manner described in the previous subsection. It is shown that the predicted element ratios, ( 16 O/ 12 C) 6917 and ( 16 O/ 12 C) 7116 , stay almost the same as the areal density of caffeine grows (see Table 1). This is because caffeine has an element ratio ( 16 O/ 12 C) of 0.25, which is very close  www.nature.com/scientificreports/ to that of cocaine. One can see that when the areal density of caffeine is less than 0.30 g/cm 2 , all predicted ratios are consistent with the theoretical values within the statistical uncertainty. As a result, the element ratio approach is still valid to identify drugs wrapped in thin materials containing C, H, O and N. However, for caffeine denser than 0.30 g/cm 2 , the predicted ratio ( 14 N/ 12 C) 7029 increases rapidly due to the fact that the ratio of ( 14 N/ 12 C) for caffeine is visibly larger than that of the cocaine.

Discussion
We have proposed a novel inspection method for achieving drug components and composition, which can be used to obtain a unique fingerprint of the drug. The proposed method is based on NRF spectroscopy plus element ratio approach, benefiting from an intense LCS γ-ray source. MC simulations show that NRF signatures of the drug components, 12 C , 14 N and 16 O, can be detected with a high confidence of 7−24σ with a γ-ray flux of 10 11 . The element ratios, 16 O/ 12 C and 14 N/ 12 C , for five widely abused drugs are then identical with the exact ratios within relative deviation of 11.6 %.
We further evaluate systematic uncertainty in element ratios predicted with Eq. (4), in view of uncertainty of NRF cross sections, instability of LCS γ-ray beam and unavailability of angular correlation W(θ) (see Section S3 of Supplementary materials). First, the predicted elemental ratios, ( 16 O/ 12 C) 6917 , ( 16 O/ 12 C) 7116 and ( 14 N/ 12 C) 7029 , have respectively uncertainties of 6.9%, 7.8% and 13.5% due to the uncertainty in NRF cross sections. Second, the instability of the LCS beam has invisible impacts on the element ratio predictions because of the intrinsic synchronous fluctuation of photon intensity at different photon energies. Third, the NRF lines at 4438 keV ( 12 C), and 6917 keV ( 16 O) and 7116 keV ( 16 O) have a purely quadrupole-quadrupole or dipole-dipole transition and the resulting W(θ) have exact expressions, which will not lead to observable uncertainties in element ratios. However, an exact expression of W(θ) is not available for the NRF line at 7029 keV ( 14 N), due to the unavailability of mixing ratio. Currently it is very hard to evaluate the effect of uncertainty in W(θ) on element ratio. Here we call for both experimental and theoretical efforts towards the accurate determination of real values of W(θ) . Moreover, the element ratio approach can hardly be used to identify benign materials with very closed element concentrations to the drug materials inspected, considering that the predicted element ratios would be the same as those of the drug materials. In this situation, the inspector should be alerted to take further actions, such as chemical identification with the Raman spectroscopy.
We should note that the γ-ray beam can probe different parts of inspected material due to the spatial nonuniformity, which makes NRF γ rays emitting from different positions of the target and then induces uncertainty to the simulated NRF yields. In our simulations, such uncertainty has been included in the relative deviations of the element ratios predicted. When the spot size of the γ-ray beam incident upon the target (see Fig. 6) is far smaller than the target dimension (see the caption of Fig. 3), it can be regarded as a kind of point-like γ-ray  www.nature.com/scientificreports/ source with good directionality. Then the spatial non-uniformity of such a beam would not impact on the fidelity of the element ratio predictions. In a realistic setup, in case the beam size is comparable to the target dimension, one should address this issue carefully. The detectable drug mass can reach gram level, given by that the target size can be shrank from 5 cm (see caption of Fig. 3) to a few mm, which is comparable to the size of LCS γ-ray beam produced by hundreds of MeV electron beam 25 . Note that such small target would not reduce the NRF yields but benefit the NRF γ-ray transportation to the HPGe detectors. These results suggest that the NRF spectroscopy combined with the element ratio approach will enable us to identify drugs in a portable package within a practical measurement time. Our proposed method can still be extended to screen containers and identify SNM relevant to weapon of mass destruction and hazardous chemical compounds like explosives. Recently, compact X/γ-ray source generators have been developed [36][37][38][39][40] , which may offer great opportunity to adopt the proposed method to seaports and national borders.

Methods
NRF γ rays. NRF refers to nuclear resonant absorption of a γ photon followed by the de-excitation with emission of γ rays. The energies of the states excited by NRF is inherent in the atomic number and the mass of the nucleus of interest (see Fig. 1b). The NRF cross section for absorption via the resonant energy level E r is given by Breit-Wigner distribution 41 : where Ŵ is the width of the level at E r , Ŵ 0 is the partial width for transitions between E r and the ground state, is the Plank constant, c is speed of light.
In practise, the NRF cross section should be calculated with taking Doppler broadening into account. Approximating the true Voigt profile to a Gaussian profile, Eq. (5) is then transformed into 42 : Here = E r k B T/Mc 2 is the Doppler width, k B is the Boltzman constant, T is the absolute temperature, and M is the mass number of the nucleus. For low-Z isotopes of interest, their fundamental widths are typically ∼ 10 meV. The effective width of the cross section after Doppler broadening increases to ∼ 10 eV. However, considering that in our case the irradiation targets exist in a molecular form, M should be replaced with the mass number of the molecule. As a consequence, the Doppler widths decrease significantly by a factor of ∼ 5. Imperfect detector resolution further broadens the detectable NRF resolution to widths of about 2-3.5 keV (see Fig. 3). The NRF cross sections of 12 C, 14 N and 16 O are discussed in Supplementary Section S2.
Due to conservation of energy and momentum, a free nucleus undergoing NRF will recoil with kinetic energy E rec determined by the Compton-like formula: where θ is photon scattering angle relative to its incident direction. It is worth to note that the recoil energy for the drug molecules is ∼ 200 eV, while the width of the Doppler-broaden NRF interaction is ∼ 2 eV. Therein, the energy of the emitting photons are too small to trigger another NRF interaction of the same transition.
NRF is generally considered to occur only between states that differ by two or fewer units of angular momentum. The angular distribution of NRF γ rays is analogous to that of γ-ray cascades. For an NRF interaction of transitions J a (L 1 )J b (L 2 )J c , with L 1 and L 2 being the multipole orders of the excitation and de-excitation, respectively, the angular distribution W(θ) can be written as 43 : where P 2n (cosθ) is the Legendre polynomial expansion, and A 2n is given by: where F 2n (L 1 J a J b ) and F 2n (L 2 J c J b ) are two constants depending on spin states of transitions and photon multipolarities 44 . The angular correlations of NRF γ rays from 12 C, 14 N and 16 O are discussed in Supplementary Section S3.
MC simulations. In this study, LCS γ-ray beam was produced by Compton backscattering of 532-nm-wavelength laser pulse with high-energy electron beam. The generation of the LCS γ-ray beam was simulated with the code MCLCSS 33 . After generation, the delivery, transportation and/or collimation of the LCS γ rays were performed with Geant4 toolkit.
To model NRF interactions, we developed a new class G4NRF into the Geant4 toolkit 34,45 . The pure virtual method G4VUserPhysicslist::ConstrcutProcess() was implemented in the simulation, and the method AddDis-creteProcess() was used for registration of the newly developed physics process. Customizing NRF process into the simulation requires the implementation of two features: first, to provide the cross sections for interaction W(θ) = 1 + A 2 P 2 (cosθ) + · · · + A 2n P 2n (cosθ), www.nature.com/scientificreports/ and second, to determine the final state resulting from the interaction. The calculation of a series of NRF cross sections was given by Eq. (6). The information of final states is provided by Eq. (8).
In all simulations, the root-mean-square (RMS) value of the electron energy spread of electron is fixed to 0.05%. The normalized electron emittance is ǫ n = 0.2 mm · mrad and the RMS beam size at the interaction point is σ e = 0.03 mm. In the scenario of 238 U irradiation, a 266-MeV electron beam was employed to generate the LCS γ-ray beam with cut-off energy of 2.5 MeV. Such un-collimated γ-ray beam was then incident upon a cylindrical 238 U target with a radius of 5.0 cm. In the scenario of drug inspection, the energy of electron beam was adjusted to 452 MeV, such that the γ-ray beam generated reaches a cut-off energy of ∼7.2 MeV. In addition, a collimator with collimation angle of 1.5 mrad was employed to avoid the drug material irradiated by LCS photons at low energies, and the beam spot size on the target surface is several mm. In this scenario, the cocaine, heroin, ketamine, methamphetamine and morphine targets are also cylindrical ones with radius of 5.0 cm and areal density of 9.0 g/cm 2 . Their compositions are specified to be C 17 H 21 NO 4 , C 21 H 23 NO 5 , C 13 H 16 ONCl, C 10 H 15 N and C 17 H 19 NO 3 , respectively. The spectral and spatial distributions of the LCS γ-ray bemas used for these target irradiations are shown in Fig. 6.
The backscatter NRF photons resulted from irradiation targets are recorded by four HPGe detectors assembled at 135 • relative to the LCS beam direction (see Fig. 1). It is supposed that the HPGe detectors have an energy resolution of 0.05% (in RMS), which is achievable with the present detector technology 46,47 . The Ge crystals have a 10 cm diameter and a 10 cm length. The full-energy peak efficiency of each HPGe detector, ǫ(E) , is simulated with the Geant4 toolkit. Accordingly, the ǫ(E) at 4438, 7029, 6917 and 7116 keV are extracted in order to perform element ratio predictions. Four HPGe detectors are shielded by a 1.0-mm-thick lead plate, avoiding a large part of low-energy photons entering into the HPGe detector.