Updates to TRACK_TEST and TRACK_VISION Computer Programs

The computer programs TRACK_TEST and TRACK_VISION were previously developed to model profiles and optical appearances of tracks developed in solid-state nuclear track detectors. The programs were based on a track development model that involved the bulk etch rate Vb and the track etch rate Vt or the V function (i.e., Vt/Vb). The present work reported our work to update and modify these two programs. In the revised TRACK_TEST, two new V functions were added and enabled. Sample results for the CR-39 detector obtained using the three original and the two new V functions were compared. Discrepancies were within ~10% and <14% for incident alpha-particle energies of 1 MeV and >1 MeV, respectively. Another major revision of TRACK_TEST was to enable calculations for the Makrofol detector. In the revised TRACK_VISION, the two new V functions, as well as the option for the Makrofol detector, were also added. The experimental results on the Makrofol detectors were obtained (irradiated with 3.6-MeV alpha particles under normal incidence and then etched to achieve a removed detector thickness of 30 μm) for comparisons with the modeled results using the revised TRACK_VISION. The track diameters obtained from the experiment and model were 24.7 and 23.2 μm, respectively. Moreover, a bright area in the central parts, together with an outer dark ring, were present in both the simulated and experimental tracks. The track-opening diameters and the general optical appearances of the tracks were in good agreement.


Introduction
Solid-state nuclear track detectors (SSNTDs) have been widely used for radon measurements through counting and/or characterizing the tracks created by the incident alpha particles emitted form radon gas and its progeny and by the subsequent chemical etching (see, e.g., References [1][2][3][4][5][6][7]). When a heavy, charged particle passes through an SSNTD, a cylinder consisting of free chemical radicals and other chemical species is formed along the trajectory, and this cylinder is referred to as a "latent track". When the irradiated SSNTD is immersed into an aggressive solution called an etchant, the latent track will undergo more intensive chemical reactions compared to the unirradiated areas, leading to formation of a "track" that can be observed under an optical microscope.
Among the most widely used SSNTDs are polyallyldiglycol carbonate (PADC) films (commercially available as CR-39 detectors), cellulose nitrate films (commercially available as LR 115 detectors) and polycarbonate films (commercially available as Makrofol detectors). Different aspects of SSNTDs were reviewed by Nikezic and Yu [8]. Optical microscopes are commonly used for counting or characterizing alpha-particle tracks. It is also noted that optical features of the etched tracks, such as their optical appearances (bright and dark areas, etc.) can provide useful information on track parameters, such as track depths, measurements of which using common optical microscopes remain difficult.
To tackle these challenges, our group developed the computer program TRACK_TEST [9] in 2006 to theoretically model alpha-particle track profiles, which enabled the determination of parameters, including the lengths of the major and minor axes, as well as the depths of the tracks. Subsequently, based on TRACK_TEST, our group further developed the computer program TRACK_VISION [10] in 2008 to theoretically model the optical appearances of such tracks. These two computer programs were based on our proposed track development model, with some modifications [11]. The model involved two parameters-namely, the bulk etch rate V b and the track etch rate V t . For a homogeneous and isotropic detector material, the bulk etch rate V b is a constant for a specified etching condition (i.e., for the chosen etching chemical, chemical concentration, etching temperature, etc.). In contrast, V t generally varies along the trajectory of the incident ion that created the latent track and will vary according to the species and energy of the ion. The ratio V t /V b is referred to as the V function.
After the publication of the TRACK_TEST and TRACK_VISION programs [9,10], the idea was extended to studies on etched tracks developed in PADC films irradiated by protons [12]. The program TRACK_P calculated the track parameters and could also graphically provide two-dimensional and three-dimensional illustrations of proton tracks developed in the PADC films.

TRACK_TEST, TRACK_VISION and TRACK_P
TRACK_TEST enabled the calculations of track parameters in CR-39 and LR 115 SSNTDs, i.e., the lengths of the major and minor axes, as well as the depths of the alpha-particle tracks. Three and two built-in V functions were provided for the CR-39 SSNTD [13][14][15] and the LR 115 SSNTD [16,17], respectively. Based on the TRACK_TEST program, another program entitled TRACK_VISION was developed for visualizing the optical appearances of tracks developed in the chemically etched SSNTDs. Based on the geometrical optics, TRACK_VISION simulated light propagation through the tracks and computed the brightness of grid elements in the track wall [5]. A total of four built-in V functions were provided for the PADC detectors [13][14][15]18].
Since their publication, the programs TRACK_TEST and TRACK_VISION have been made available on the websites http://www.cityu.edu.hk/nru/Track_Test.htm and http: //www.cityu.edu.hk/nru/Track_Vision.htm (accessed on 13 February 2021), respectively. As of 30 January 2021, the programs TRACK_TEST and TRACK_VISION have been cited 119 and 45 times, according to Google Scholar. Although there have been advancements in the field, such as the publication of new V functions, as well as continuous requests from various research groups for the creation of modifications and new options, e.g., for modeling for the Makrofol detector, the two programs TRACK_TEST and TRACK_VISION have not yet been updated or modified since their first publications. This current work is devoted to such updates and modifications.

Programming and Simulation
Detailed information on programming and simulation can be found in Reference [9], but for completeness, the basic points are briefly reviewed here. TRACK_TEST first determines the particle range in the detector and the depth D along the particle trajectory to which the etchant penetrates. Coordinates of points on the two-dimensional track wall for normal incidence are then computed, which are utilized to generate the track profile for non-normal incidence through the translation and rotation of the coordinate system [11]. A three-dimensional track is subsequently rendered through the rotation of the points on the two-dimensional track wall around the particle trajectory, from which the track opening contour, as well as the major and minor axes, are determined.

Adding Two V Functions
In this revised version of TRACK_TEST, two new V functions (expressed through the residual range R , which is defined as the distance from a specific point on the particle track to its end point) were added and enabled. The first added V function was adopted from the work of Al-Jubbori [19] in the form with the constants a 1 = 0.098, a 2 = 1.86, a 3 = 37.78, a 4 = 36.98 and a 5 = 0.98 obtained by fitting the data provided by measuring the average lengths of the tracks (observed using an Optika B-193 microscope, Optika, Ponteranica, Italy) for different alpha-particle energies and etching times. In the revised TRACK_TEST, this function can be invoked by typing No. 4 when the program asks for the number of the V function. The second added V function was adopted from the work of Hermsdorf [20] in the form with constants a 1 = 390, a 2 = 2, a 3 = 1, a 4 = 5, a 5 = 80 and b 1 = 2.35. In the revised TRACK_TEST, this function can be invoked by typing No. 5 when the program asks for the number of the V function.
To summarize, together with the three built-in V functions originally provided in the previous version of TRACK_TEST for the CR-39 SSNTD, the five built-in V functions provided in the revised TRACK_TEST for the CR-39 SSNTD are: No. 1 from Reference [13], No. 2 from Reference [14], No. 3 from Reference [15], No. 4 from Reference [19] and No. 5 from Reference [20]. Figure 1 shows the relationship between the major-axis length (D major ) of the track opening in the CR-39 detector and the alpha-particle energy for the incident angle of 80 • . The results obtained using the five V functions are compared. Different results on D major were noticed. However, the discrepancies were not very significant. For example, for the incident alpha-particle energy of 1 MeV, all the results were in the range 12.2 ± 1.2 µm, which corresponded to discrepancies of only about 10%. For larger incident alpha-particle energies, the results were more scattered, in the range of 9.5 ± 1.3 µm, which, however, still corresponded to discrepancies of less than 14%.

Adding Makrofol Option
Another major revision of TRACK_TEST was to enable calculations for the Makrofol detector. In comparison to the CR-39 and LR115 detectors, there is, in general, a lack of data for the Makrofol detector. Even measurements on V b have not been common. One V function established for the Makrofol detector was given by Stajic et al. [21].
However, an unexpected behavior of this function was observed for low values of the residual range (i.e., the function was not equal to 1 when R → 0), so a slight modification was necessary to facilitate calculations by the computer programs. Furthermore, this V function led to numerical problems and incorrect plotting of the track profiles for overetched tracks. These difficulties prompted us to find a better V function for Makrofol detectors without such drawbacks. The data obtained using the method described by Stajic et al. [21] by fitting to the new function: Polymers 2021, 13, 560 4 of 9 gave fitting parameters a i , b i and c i , given in Table 1. The coefficient of determination of the fit was R 2 = 0.9998. In more detail, the data were calculated using the expressions involving the restricted energy loss [22,23], output data from SRIM-2013 [24] and the parameters provided by Somogyi et al. [23] for polycarbonate detectors. This new function was finally implemented in the modified programs in the present work.

Adding Makrofol Option
Another major revision of TRACK_TEST was to enable calculations for detector. In comparison to the CR-39 and LR115 detectors, there is, in gen data for the Makrofol detector. Even measurements on Vb have not been co function established for the Makrofol detector was given by Stajic et al. [21] However, an unexpected behavior of this function was observed for the residual range (i.e., the function was not equal to 1 when R′ → 0), so a cation was necessary to facilitate calculations by the computer programs. this V function led to numerical problems and incorrect plotting of the trac over-etched tracks. These difficulties prompted us to find a better V function detectors without such drawbacks. The data obtained using the method Stajic et al. [21] by fitting to the new function:  The bulk etch rate V b for the Makrofol detector was also determined by Stajic et al. [21]. For the PEW (potassium, ethanol and water) etchant (15-g KOH + 45-g H 2 O + 40-g ethyl alcohol), V b was found as 15 ± 1 µm/h. The relationships between the lengths of the major and minor axes of the track opening in the Makrofol detector and the alpha-particle energy are shown in Figure 2 for a removed layer of 30 µm and both incident angles of 60 • and 90 • . pha-particle energy are shown in Figure 2 for a removed layer of 30 μm and bo angles of 60° and 90°.

Programming and Simulation
Detailed information on the programming and simulation can be found i [10], but for completeness, the basic points are briefly reviewed here. TRAC simulated the appearance of tracks under the transparent mode of the op scope, where light rays came up from below the detector to strike the detec perpendicularly. The points on the three-dimensional track wall gen TRACK_TEST were linked to form "four-angle polygons", to which ray-tra dures were applied [25,26]. Any light ray exiting from the detector could th acterized by its direction and relative intensity. The relative intensities for all e the three-dimensional track wall were then determined.

Adding Two V Functions and Makrofol Option
The TRACK_VISION program was also modified to accommodate the functions for PADC films and to provide the new option for calculations for th detector. Both new V functions for PADC detectors described in Equations were programmed. A simulation of the optical appearance of tracks develo Makrofol detector was enabled. A sample image of the alpha-particle tracks on an irradiated and etched Makrofol detector captured using an optical mi shown in Figure 3. The detector was irradiated with alpha particles with an en MeV under normal incidence (using an electroplated 241 Am source and a colli

Programming and Simulation
Detailed information on the programming and simulation can be found in Reference [10], but for completeness, the basic points are briefly reviewed here. TRACK_VISION simulated the appearance of tracks under the transparent mode of the optical microscope, where light rays came up from below the detector to strike the detector surface perpendicularly. The points on the three-dimensional track wall generated by TRACK_TEST were linked to form "four-angle polygons", to which ray-tracing procedures were applied [25,26]. Any light ray exiting from the detector could then be characterized by its direction and relative intensity. The relative intensities for all elements on the three-dimensional track wall were then determined.

Adding Two V Functions and Makrofol Option
The TRACK_VISION program was also modified to accommodate the two new V functions for PADC films and to provide the new option for calculations for the Makrofol detector. Both new V functions for PADC detectors described in Equations (1) and (2) were programmed. A simulation of the optical appearance of tracks developed in the Makrofol detector was enabled. A sample image of the alpha-particle tracks developed on an irradiated and etched Makrofol detector captured using an optical microscope is shown in Figure 3. The detector was irradiated with alpha particles with an energy of 3.6 MeV under normal incidence (using an electroplated 241 Am source and a collimator with an appropriate thickness). The incident energy (3.6 MeV) was estimated based on the distance in the air between the source and the detector (i.e., the collimator thickness) and a function obtained by fitting SRIM-2013 output data for helium in the air. The detectors were etched for 2 h, and the removed detector thickness was determined as 30 ± 2 µm.  Corresponding to those alpha-particle tracks developed on an irradia Makrofol detector shown in Figure 3, the simulated results obtained b TRACK_VISION program are shown in Figure 4. The left panel of Figur track profile and track opening (the latter being a circle for normal incide using the revised TRACK_VISION, together with the position of the de before and after etching. On the other hand, the right panel of Figure 4 sh sponding optical track appearance. The diameter of the experimentally pha-particle tracks presented in Figure 3 was estimated to be 24.7 ± 1.4 µ m of the scale in Figure 3. The diameter given by the revised TRACK_VISION 23.2 μm, and the two values were deemed to have good agreement.  Corresponding to those alpha-particle tracks developed on an irradiated and etched Makrofol detector shown in Figure 3, the simulated results obtained by the revised TRACK_VISION program are shown in Figure 4. The left panel of Figure 4 shows the track profile and track opening (the latter being a circle for normal incidence) simulated using the revised TRACK_VISION, together with the position of the detector surface before and after etching. On the other hand, the right panel of Figure 4 shows the corresponding optical track appearance. The diameter of the experimentally obtained alpha-particle tracks presented in Figure 3 was estimated to be 24.7 ± 1.4 µm with the help of the scale in Figure 3. The diameter given by the revised TRACK_VISION program was 23.2 µm, and the two values were deemed to have good agreement. Corresponding to those alpha-particle tracks developed on an irradiated and etched Makrofol detector shown in Figure 3, the simulated results obtained by the revised TRACK_VISION program are shown in Figure 4. The left panel of Figure 4 shows the track profile and track opening (the latter being a circle for normal incidence) simulated using the revised TRACK_VISION, together with the position of the detector surface before and after etching. On the other hand, the right panel of Figure 4 shows the corresponding optical track appearance. The diameter of the experimentally obtained alpha-particle tracks presented in Figure 3 was estimated to be 24.7 ± 1.4 µ m with the help of the scale in Figure 3. The diameter given by the revised TRACK_VISION program was 23.2 μm, and the two values were deemed to have good agreement. As regards the optical appearance of the tracks, there were similarities and differences between the simulated and experimental results. Notably, a bright area was present in the central parts of both simulated and experimental tracks, which was explained by As regards the optical appearance of the tracks, there were similarities and differences between the simulated and experimental results. Notably, a bright area was present in the central parts of both simulated and experimental tracks, which was explained by the light passing through the tracks without a significant loss of intensity. Interestingly, however, the appearances of the tracks close to their perimeters were different for the simulated and experimental tracks. The simulated track shown in the right panel in Figure 4 showed an outer dark ring caused by a total reflection of the light coming from the bottom of the detector. In the experimental results shown in Figure 3, dark rings were also observed close to the perimeters of the tracks but were enclosed by bright interference fringes. It is remarked here that the interference of light was not taken into account in the present work, where only the geometrical optics were considered. The incorporation of wave optics will be carried out in our future works.

Notes on Installing and Running TRACK_TEST 2.0 and TRACK_VISION 2.0
The revised TRACK_TEST and TRACK_VISION programs are referred to as TRACK_TEST 2.0 and TRACK_VISION 2.0. These programs are freeware downloadable from the websites http://www.cityu.edu.hk/nru/Track_Test.htm and http://www.cityu.edu.hk/nru/ Track_Vision.htm (accessed on 13 February 2021), respectively, and are compressed in ZIP format. After downloading, the files for a chosen program (TRACK_TEST 2.0 or TRACK_VISION 2.0) should be extracted and stored in a directory created by the user. Note that different directories should be created for files from the two different programs. Further compilation is not needed. When a program is started, it will read input data from the file INPUT.DAT, which can be edited using any text editor. The programs enable calculations for three kinds of the most frequently used detectors-namely, the CR-39 detector and the LR 115 detector, as well as the Makrofol detector. The user can choose the detector by typing "C" for CR-39, "L" for LR 115 and "M" for Makrofol in INPUT.DAT. More detailed instructions on running the programs can be found on their respective webpages.
Upon successful running of the programs, the lengths of the major and minor axes of the track opening, as well as the track depth, are given. The following will also be shown on the computer screen: initial detector surface, detector surface after etching, track profile, contour of track opening and trajectory of the particle, as well as the scale for determining the real dimensions from the images. The aspect ratio of the screen should be set to 4:3 in order to show undistorted images. For TRACK_VISION 2.0, the simulated appearance of a track under the transparent mode of the optical microscope will also be displayed.

Conclusions
The current work provided updates and modifications to two previously developed computer programs: TRACK_TEST and TRACK_VISION.
The computer program TRACK_TEST was first developed in 2006 [9] to model alphaparticle track profiles that enable the determination of parameters including the lengths of major and minor axes, as well as the depths of the tracks. In the revised version of TRACK_TEST, two V functions were added, which were taken from the research works of Al-Jubbori [19] and Hermsdorf [20]. Moreover, the Makrofol option was added by adopting the V function from Stajic et al. [21], with some modifications. To differentiate between the original and revised versions of TRACK_TEST, the current revised version is entitled TRACK_TEST 2.0, and both the original version TRACK_TEST and the revised version TRACK_TEST 2.0 are available on the website http://www.cityu.edu.hk/nru/Track_Test. htm (accessed on 13 February 2021).
The computer program TRACK_VISION was developed in 2008 [10] to model the optical appearances of such tracks. In the revised version of TRACK_VISION, the two V functions for the CR-39 detector established by Al-Jubbori [19] and Hermsdorf [20] were also added, and the Makrofol option was also enabled. To differentiate between the original and revised versions of TRACK_VISION, the current revised version is entitled TRACK_VISION 2.0, and both the original version TRACK_VISION and the revised version TRACK_VISION 2.0 are available on the website http://www.cityu.edu.hk/nru/Track_ Vision.htm (accessed on 13 February 2021).