Computational evaluation of a novel beta radiation probe design using integrated circuits

Researchers at Texas A&M University (TAMU) have designed the radiation integrated circuit (RIC) for deployment as a new radiation detection system. Most integrated circuits are susceptible to radiation-induced failures, and decades of research have gone into solving this problem. Research at TAMU has led to a novel integrated circuit design that utilizes both radiation-hardened areas (RHAs) and radiation-sensitive areas (RSAs) to take advantage of these failures. The RSAs are susceptible to charged particle interactions, allowing the RIC to detect alpha and beta particles. However, beta particles are more penetrating compared to alpha particles, resulting in a lower interaction probability for beta particles incident on a bare RIC. In any material, the higher the beta energy, the deeper the beta particle can penetrate; therefore, the use of a wedge-shaped attenuator for beta particle detection not only increases interaction probability, but also provides the capability to perform maximum beta energy discrimination in the field. The objective of this research was to optimize the design of the RIC. Monte Carlo N-particle radiation transport code (MCNP) simulations assessed the beta particle detection and maximum energy discrimination performance of plate glass, borosilicate (Pyrex®) glass, acrylic (Lucite®), and natural rubber attenuators. In this proof-of-concept analysis, natural rubber was observed to be the optimal attenuating material for the beta probe with respect to maximum energy discrimination capability and weight, but all materials considered proved to be good candidates. The results of this study are promising and indicate the potential to achieve maximum beta particle energy discrimination of 50 keV using a wedged, natural rubber attenuator on the RIC.


Background
With the miniaturization of semiconductor technology, circuits on chips have become increasingly sensitive to radiation damage and, therefore, more likely not to function properly in a high radiation field. Researchers at Texas A&M University (TAMU) have been able to convert this failure into a feature with the radiation integrated circuit (RIC). The RIC contains both radiation-hardened areas (RHAs) and radiation-sensitive areas (RSAs) [1]. Figure 1 provides a simple overview of the RIC design. The RSAs contain electrical circuits which are sensitive to charged particle interactions [2,3]. A single RHA, designed to be highly resistant to radiation effects, monitors multiple RSAs for changes in voltage, indicative of a charged particle interaction [2,3]. If an RHA registers no voltage change in a particular RSA, the RHA moves on to observe the next RSA in the chain. This cycle continues until a voltage change is detected in an RSA. An RSA that suffers a charged particle interaction will recover to its normal state within several nanoseconds, allowing an RHA to monitor about 40 RSAs, given a worstcase internal clock frequency of 1 GHz.
Hundreds of thousands (~ 10 5 ) of RSAs can be fabricated on a 1 × 1 cm RIC, comprising about 75% of the chip area [1,3]. Tens of thousands (~ 10 4 ) of RHAs consume about 15% of the chip area, and the remainder of the RIC is used to implement signal processing and clocking circuits. The distribution of the total active area into many distinct, active regions renders the likelihood of multiple particle strikes in a single region between resets very low, especially in the low activity fields a survey probe of this type would be deployed in. As such, the system dead time is expected to be negligible. Survey probes utilizing the RIC can offer several advantages over current systems: reduced power consumption, reduced weight, relatively low cost when produced in bulk, and significantly smaller packaging because detection and signal processing are handled on a single chip [1][2][3].

Working principle of the probe
With RSAs sensitive to charged particles, RICs can be employed to detect beta particles. However, beta particles are far more penetrating than alpha particles, reducing the probability of a beta particle interaction with a bare RIC. Utilizing a wedge-shaped attenuator (Fig. 2) for beta particle detection not only increases the interaction probability, but also provides the capability to perform maximum beta energy (E β,max ) discrimination in the field. Beta particle penetration (range) in matter is a function of the particle's energy and the material's properties. For a survey probe comprised of RICs, the attenuator thickness over each sector of the chip will be known, allowing the determination of E β,max for beta contamination in the field. When the probe is fabricated, RICs will be placed on the inclined surface of the wedge such that beta particles will penetrate varying lengths of attenuator before arriving at the RIC surface. Figure 2 exhibits the working principle used to determine E β,max . In the figure, a beta particle source is depicted below the probe base, shown along with the wedge-shaped attenuator. A beta particle of known energy has a known maximum penetration depth (MPD) in any given material. In the figure, an arc is shown which represents the MPD in the material. The arrowed lines indicate the tracks of five identical beta particles, labeled 1 through 5, that strike the probe at different angles. For beta particles 1 and 5, the attenuator thickness is greater than the MPD. These particles will not reach the chip. For beta particle 3, the attenuator thickness is less than the MPD. This particle, which will penetrate the wedge, can either pass through the chip without affecting an RSA or interact in the chip. For beta particles 2 and 4, the attenuator thickness is slightly less than the MPD in the material. These particles will be near the end of their track and have minimal energy near the RIC surface, providing the highest probability of interaction. For this research, areas on the RIC where the attenuator thickness is slightly less than the MPD are jointly referred to as the region of interest (ROI). This ROI is a function of the beta particle source position along the probe, E β,max , and the attenuating material. Analyses of different attenuator materials were focused on interactions in the ROI, as explained in the next section.
It is important to note that the scope of this research is only to introduce a novel beta particle detector design and to explore through simulation how it could potentially be used to indicate the energy of the beta particles it registers. This research does not address the challenges of manufacturing or operational use in the field such as more complex source geometries, the presence of multiple source isotopes, and the added complexity of energy discrimination given the continuous energy spectrum of beta decay. The simulation work presented in this paper was used as a basis for creating the RIC which has since been fabricated. The challenges of manufacturing  Illustration of maximum beta particle penetration distance in a wedge-shaped attenuator, the phenomenon by which the RIC can perform maximum energy discrimination a finalized design and field use in complex beta radiation fields will be addressed in follow-on research.

Materials and geometry
Simulations of the RIC were carried out using Monte Carlo N-particle radiation transport code, MCNP6.1 [4]. In these simulations, the attenuator material for the beta survey probe was modeled as a wedge. Four low atomic number (low-Z) materials were considered in this research: borosilicate (Pyrex®) glass, plate glass, acrylic (Lucite®), and natural rubber. The isotopics of the materials were acquired from Pacific Northwest National Laboratory's Compendium of Material Composition Data for Radiation Transport Modeling, henceforth referred to as The Compendium [5]. Isotope ZAID (atomic number, Z, and atomic mass, A, identification) values, weight fractions, and the density of each material are summarized in Table 1. Several simulations were carried out for each attenuator material to assess its energy discrimination capability by varying the source position and E β,max .
As a handheld survey instrument, a probe comprised of RICs must balance energy discrimination capability, physical size, and weight. Low-Z materials were selected due to their superior beta particle stopping power in the target energy regime and low density. For low incident energies (less than about 20 MeV) like those considered in this study, beta particles lose energy more rapidly in low-Z materials due to their higher collision rate compared to high-Z materials [6]. It is only at higher incident energies that the superior energy reduction per collision of high-Z materials dominates; therefore, beta particles have shorter range in low-Z materials within the energy regime considered. Probe weight was an especially significant concern which is why only very lightweight, organic materials were ultimately selected. Organic materials are known to degrade (yellowing, loss of elasticity, etc.) when irradiated to sufficient doses. For example, many polymers and elastomers for the medical field are considered radiation compatible up to doses on the order of tens of kGy [7]. This is several orders of magnitude greater than the annual dose limit to radiation workers in the United States of 5000 mrem or 0.05 Gy [8]. Radiation fields intense enough to accumulate to degrading doses within a timeframe that would bring the probe's utility into question would also prevent human entry. Additionally, field instruments of the type envisioned are characterized at regular intervals to check for performance changes. As such, material degradation was not considered in this investigation. Table 2 lists the modeled dimensions for each attenuator material. The w, l, and h values represent the width, length, and height (maximum thickness) of the different attenuators as shown in Fig. 3. For example, the plate glass attenuator was 5 cm long by 5 cm wide, with a thickness that increased linearly along the length from 0 to 0.5 cm. For each material, the attenuator dimensions were selected to allow E β,max discrimination of beta particles with energies in excess of 2.1 MeV. For each attenuator material, the linear stopping power for beta particles varies, so the attenuators' dimensions were selected accordingly assuming linear penetration, which is a rough approximation used in working problems [6].

MCNP simulations
A surface flux detector (tally type F2 per MCNP code nomenclature) was created along the inclined surface of the wedge, centered around the Y-axis. Tally segments (the FS card in MCNP) were utilized to split the detector (tally) into many bins of interest. Each bin had a width along the X-axis of 0.1 cm and a length along the Y-axis of 0.1 cm, resulting in a bin area slightly greater than 1 mm 2 on the inclined surface. Such an area on the RIC chip contains several thousand RSAs. These bins lay flush to one another such that a 5 cm long probe model has 50 tally bins along its inclined surface. Simulations were carried out to evaluate the E β,max discrimination capability of each attenuator material. For this assessment, the initial energy of beta particles was varied, and the position of the source was varied accordingly. For each material, an isotropic, monoenergetic beta particle point source was simulated with energies ranging from 1.95 to 2.10 MeV in increments of 0.01 MeV.
In each simulation, the point source was placed 1 cm below the attenuator (Z = − 1 cm) and centered along the X-axis (X = 0 cm). The source was also positioned along the Y-axis such that the attenuator thickness directly above the source equaled the material's MPD for the beta particle energy under consideration. MPD was calculated with the empirical formula given by Eq. 1 [6]. This empirical calculation of the MPD (Turner [6] uses the nomenclature, "range") assumes the continuous slowing down approximation (CSDA) and "straight-ahead travel." The result is an overestimation of electron penetration, which will shift the peak electron fluence rate along the inclined surface in the positive-Y direction where the attenuator is less thick. Since the purpose of the study is to observe separation of the peaks at different energies, not to locate peaks precisely, these assumptions are acceptable for this research so long as the shift of expected peak positions is consistent across the energy range considered. The monoenergetic beta energy in the above simulations represents the maximum energy of all beta particles emitted. In reality, however, beta decay is a three-body problem, resulting in a continuous energy spectrum for beta particles. This spectrum can be predicted by the Fermi theory [9].
On a manufactured probe utilizing RICs, energy will be indicated by the position of the particle fluence rate peak along the detector surface, so energy discrimination can be measured by the spacing of peaks of different energies along the surface. The location of the fluence rate peak for each energy in the 1.95 to 2.10 MeV range was located, and the spatial energy resolution of the detector was calculated as the negative ratio of the difference of maximum and minimum energies to the difference of maximum and minimum peak locations along the Y-axis as depicted in Eq. (2). The negative ratio is required because the peak position of the maximum energy is less than that of the minimum energy in the chosen coordinate system. By this formulation, a lower value of spatial energy resolution indicates better energy discrimination capability because there is greater distance between the peaks. (1)

Results and discussion
In MCNP6.1, plate glass, borosilicate (Pyrex®) glass, acrylic (Lucite®), and natural rubber were simulated as wedgeshaped attenuators to characterize the maximum beta energy discrimination performance of a RIC-based survey instrument. In each simulation, the number of beta particles traversing the attenuator and eventually entering RSAs of the RIC was tallied. The analysis of the tallies focused on the ROIs for each attenuator.

Plate glass
The fluence rate tallies as a function of scoring bin position along the Y-axis are presented in Fig. 4. Vertical error bars are displayed on the data points, but most are not visible because they are very small (< 1%). For clarity, only the results for initial energies of 1.95, 2.00, 2.05, and 2.10 MeV are plotted. No tally multipliers were applied, so the results are representative of a 1 Bq source, hence the low fluence rate magnitudes. The relative errors of most data points were well under 10%, satisfying the statistical check for a surface detector (F2 tally) in MCNP. Some tallies on the thickest portion of the detector (0 cm ≤ Y ≤ 0.5 cm) had relative errors in excess of 10%. The fluence rates in this region are very low, so the error bars are not visible.
However, this region is isolated from the nearest energy peak studied in this research by about 1 cm, so the slightly high relative errors are not of concern. A connecting line is shown to highlight the trend of the results. The trend in the data is dominated by the attenuation of beta particles in the material: higher energy peaks occur on the left where the material is thicker, and lower energy peaks occur on the right where the material is thinner. To the left of the peaks, most particles do not penetrate all the way through the thicker portion of the material, and as a result, the fluence rate values are much lower. To the immediate right of the peaks, the beta particles can penetrate the attenuator and be tallied, but further to the right the total transport distance increases, reducing the fluence rate.
To the far right of the peaks, the MPD becomes increasingly greater than the attenuation thickness, so the particles pass all the way through the material, resulting in increased fluence rate. This general behavior is illustrated in Fig. 5, which depicts the relationship between MPD and the linear travel distance of a particle. The empirical MPD formulation of Eq. 1 was used to calculate the MPD plotted in Fig. 5, slightly overestimating the MPD. The increasing fluence rate to the far right of the peaks is a result of simulation setup and would not be observed in a probe. The surface detector (MCNP F2 tally) was implemented to allow simplified problem geometry and reduce computation time. However, the F2 tally also scores higher energy beta particles passing through the thinner portions of the material rather than only lower energy particles expected to interact with the RIC. The result is an overestimation of particle interactions to the far right of the peaks. This follows from the fact that the thickness of attenuator material traversed by particles in this region is much less than the Fig. 4 A plot of the beta particle tallies of all scoring regions along the length of the probe with a plate glass attenuator Fig. 5 The attenuation thickness that a beta particle traverses between the source and the inclined surface of the plate glass attenuator as a function of bin position on the Y-axis. This plot is for a 2 MeV source placed at (X, Y, Z) = (0, 1.046, −1) cm MPD as calculated from the CSDA. Since the focus of this research is the separation of peaks, this is an acceptable discrepancy but must be taken into account if mathematical models are to be developed describing the detector response.
The MPD for a 2 MeV beta particle in plate glass is 0.395 cm. From 1.05 cm along the probe length all the way to 5.00 cm, the attenuator thickness for a beta particle emitted from a source at (X, Y, Z) = (0, 1.046, − 1) cm is less than the MPD (Fig. 5). However, between 1.05 cm and 2.35 cm along the probe length, the attenuator thickness is only slightly less than the MPD of a 2 MeV beta particle. After traversing this distance, the beta particles are at the end of their track and within an energy range that can more readily interact in the RIC RSA. Similar trends were observed for all beta energies simulated.
For the plate glass attenuator, the ROI is between 1.05 and 2.35 cm. The defined probe region encompasses fluence rate maxima for all initial beta energies from 1.95 MeV to 2.10 MeV. Figure 6 displays the same data plotted in Fig. 4, focusing on the ROI (1.05-2.35 cm). In the ROI, the trend for the attenuator thickness was convex upward as a function of the probe length (Fig. 5). In this region, there is a local minimum attenuator thickness which corresponds to the scoring bin containing the maximum number of counts. For the 2 MeV beta particles, the peak count was estimated across the tally area at 1.75 cm into the probe length. The closer look provided by Fig. 6 reveals that each energy peak occurred within a distinct tally bin, critical for energy discrimination. Fig. 6 A plot of the beta particle tallies of all scoring regions in the ROI of the probe with a plate glass attenuator Fig. 7 A plot of the beta particle tallies of all scoring regions in the ROI of the probe with a borosilicate glass attenuator Fig. 8 A plot of the beta particle tallies of all scoring regions in the ROI of the probe with an acrylic attenuator Fig. 9 A plot of the beta particle tallies of all scoring regions in the ROI of the probe with a natural rubber attenuator

Discrimination analysis
Figures 7, 8, and 9 display the number of beta particles tallied in the ROIs of borosilicate (Pyrex®) glass, acrylic (Lucite®), and natural rubber, respectively. The ROIs (with respect to the 2 MeV reference) were 0.75-2.05 cm for Pyrex® glass, 1.55-2.85 cm for Lucite®, and 0.45-1.85 cm for natural rubber. Analogous to Fig. 6 for plate glass, the plotted values are the beta particle fluence rate tallied in each bin [cm −2 s −1 ], presented as functions of the initial beta particle energy and distance along the probe length. The trends in these figures are similar to the trend observed in Fig. 6 with the key difference being the MPD of the beta particles in each material. The MPD for 2 MeV beta particles is 0.426 cm, 0.798 cm, and 1.032 cm in Pyrex® glass, Lucite®, and natural rubber, respectively.
The energy peaks in the plate glass (Fig. 6) were all separated by 0.1 cm, a single tally bin in the chosen discretization. The energy peaks in the borosilicate glass (Fig. 7) and acrylic (Fig. 8) were separated by 0.1-0.2 cm (1 or 2 bins). The energy peaks in the natural rubber ( Fig. 9) were all separated by 0.2 cm (2 bins). Table 3 summarizes the energy discrimination performance and mass for each attenuator material. All materials investigated performed comparably, but natural rubber was observed to show the best spatial energy resolution (lowest energy change per unit distance, 0.25 MeV cm −1 ). The better energy resolution for natural rubber can be explained by its low effective atomic number (5.41) and effective mass number (10.71) compared to the other attenuator materials. For low initial particle energies (less than 20 MeV) like those considered in this investigation, beta particles lose energy more rapidly in materials of low atomic number due to the higher collision rate compared to high-Z materials [6]. The substantial hydrogen content of natural rubber reduces the effective atomic number and enhances the attenuation of beta particles in the material. Additionally, with a range of 13.9-20.8 g, all attenuators had similarly low masses. Natural rubber performed the best in these simulations, but all four materials investigated are good candidates for a RIC survey probe. Since the beta particle count drives the signal generated in the RSA, a large difference in counts from one position to the next along an attenuator makes energy discrimination easier to perform.
This calculated energy resolution serves as an estimator of hardware performance. The manufactured RIC features 70 distinct detection elements arranged parallel to one another, 52 of which can be monitored individually (Fig. 10). Hardware limitations and design priorities prevent all 70 elements from being monitored individually. For simplicity, if the 52 monitorable elements are presumed to be spaced evenly along the 1 cm width of the RIC, the RIC has a resolution of about 4.8 keV per element as given by the spatial energy resolution of 250 keV cm −1 for a natural rubber attenuator. If two tally bins (about 10 detection elements) of spacing between peaks are required for reliable peak discrimination-the worst-case scenario in the simulations-then E β,max discrimination to about 50 keV may be possible using a natural rubber wedged attenuator.
The natural rubber wedge design was 7 cm in length and 5 cm in width and the height increased linearly from 0 to 1.1 cm (Table 2). This design can measure E β,max in excess of 2.1 MeV. The range of E β,max detection can be extended by simply increasing the length and height. Conversely, an extremely compact design option is achievable but will reduce the range of measurable E β,max values. For example, to detect E β,max in the range of 0.5 to 2.1 MeV, the wedge height (thickness) would necessarily range from 0.18 to 1.1 cm and the length would be reduced to 5.87 cm.

Conclusions
The RIC RSAs are sensitive to charged particles, allowing them to detect alpha particles and beta particles. Beta particle detection with the RIC is more challenging because they are more penetrating than alpha particles and are less likely to interact on the bare integrated circuit. In order to optimize beta particle detection, attenuator materials can be placed on a RIC to reduce incident energy so that the probability of interaction is increased. Another advantage of using an attenuator is that the maximum depth a beta particle can penetrate is a function of its incident energy and the material properties. In this research, MCNP6.1 simulations were carried out to investigate the use of wedge-shaped attenuators made of plate glass, borosilicate (Pyrex®) glass, acrylic (Lucite®), and natural rubber for E β,max discrimination. If energy discrimination is achievable, survey probes made of RICs will have considerable advantages over existing in situ beta particle detection technology.
To evaluate the energy discrimination capability of each material, a spatial energy resolution was calculated. This ratio reflects the spatial separation of energy. The finer the spatial energy resolution, the better the E β,max discrimination capability. The natural rubber wedged attenuator had the best spatial energy resolution of 0.25 MeV cm −1 compared to 0.375, 0.375, and 0.5 MeV cm −1 for Pyrex®, Lucite®, and plate glass, respectively. However, these values are comparable. Likewise, the masses of the attenuators made of each material were comparable, ranging from 13.9-20.8 g. As such, each material investigated is a good candidate for a RIC survey probe. The calculated spatial energy resolution serves as an estimator of energy discrimination capability of the RIC. Thereby, there exists a potential to perform E β,max discrimination within 50 keV using a natural rubber wedged attenuator.
For this research, only a wedge-shaped attenuator was used to investigate the possibility of E β,max discrimination. The authors acknowledge that other attenuator geometries may be more efficient and such possibilities will be investigated in future work.

Compliance with ethical standards
Conflict of interest The authors have no conflict of interest to declare.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.