Abstract
High-energy, high-dose, microfocus X-ray computed tomography (HHM CT) is one of the most effective methods for high-resolution X-ray radiography inspection of high-density samples with fine structures. Minimizing the effective focal spot size of the X-ray source can significantly improve the spatial resolution and the quality of the sample images, which is critical and important for the performance of HHM CT. The objective of this study is to present a 9 MeV HHM CT prototype based on a high-average-current photo-injector in which X-rays with about 70μm focal spot size are produced via using tightly focused electron beams with 65/66μm beam size to hit an optimized tungsten target. In digital radiography (DR) experiment using this HHM CT, clear imaging of a standard 0.1 mm lead DR resolution phantom reveals a resolution of 6 lp/mm (line pairs per mm), while a 5 lp/mm resolution is obtained in CT mode using another resolution phantom made of 10 mm ferrum. Moreover, comparing with the common CT systems, a better turbine blade prototype image was obtained with this HHM CT system, which also indicates the promising application potentials of HHM CT in non-destructive inspection or testing for high-density fine-structure samples.
Keywords
Introduction
For half a century, X-ray computed tomography (CT) has been attracting growing interest in the manufacturing industry, clinical application and other non-destructive testing (NDTs), since it allows both the external geometrical and internal features of an object to be measured [1–4]. In many scenarios, such as the precise and fast detection of a large workpiece with high density and fine structures, high spatial resolution, high penetrating power, and high-dose of an X-ray CT are required simultaneously. Unfortunately, such an ideal X-ray CT with all these merits has not been created to date.
Several methods have been explored to improve the overall performance of X-ray CT. Among them, low-energy microfocus X-ray CT [5, 6] can achieve excellent resolution due to small focal spot size, even less than 1μm [7, 8], but the penetrability is limited by low beam energy (equipped with X-ray tubes with a few tens of to hundreds of kV). Therefore, it is mainly used to detect low density or small size objects such as biological bodies [9] or circuit chips. The high-energy X-ray CT, with a MeV X-ray source generated by the linear accelerators (LINAC), has much more industrial applications to detect high-density objects [10]. The biggest challenge of high-energy industrial CT is to get high-resolution for the large and high-density objects. At present, the typical focal spot of high energy X-ray CT is 1.5–2.0 mm [11], which limits the spatial resolution to 1.5–2 lp/mm and hinders its applications. The spatial resolution is determined by the X-ray source focal spot size, the detector pixel size, and the magnification factor of the imaging system [12]. However, high-energy X-ray CT radiation detectors are limited by detect quantum efficiency, machining performances, etc., which means its pixel size is difficult to improve significantly [13, 14]. Therefore, high-energy, high-dose, microfocus X-ray computed tomography (HHM CT) becomes one of the most effective methods for high-resolution X-ray radiography inspection of high-density samples with fine structures.
In addition, new attempts to minimize the focal spot size of X-ray source have been carried out recently, including minimizing the bremsstrahlung target to about 100μm [15, 16], using bremsstrahlung or all-optical inverse Compton scattering driven by laser-accelerated electron beams [17–20] with ultra-small focal spot size, and so on. These efforts improve the resolution of high-energy CT significantly, but none of them could provide a high dose rate enough. To achieve a high dose rate, bremsstrahlung with high average current electron beams is always preferred. Moreover, the electron source should have low emittance, low energy spread, and low peak current to be strongly focused to micron-scale to generate microfocus X-rays. Thus, Rhodotron [21–23] and photo-injector with high-average-current [24] are suitable for the electron source of HHM CT. In particular, the photo-injector used to drive a high average power free electron laser can generate the highest quality electron beam [25–28].
In this paper, we explore the feasibility of a 9 MeV HHM CT based on a novel high-average-current photo-cathode injector, which holds great promise to generate dream X-rays with high dose and micron-scale focal spot size. In Sec. II, we report the detailed design and implementation of a 9 MeV HHM CT prototype based on the beamline of Chengdu CAEP Terahertz Free Electron Laser (CTFEL) facility [29]. In Sec. III, we present the proof-of-principle imaging experiments carried out with this HHM CT. Discussion of these results and future plan of HHM CT system are given in Sec. IV.
Design and implement of high-energy high-dose microfocus X-ray CT
In a CT system, the cutoff frequency f co [lp/mm] = (BW [mm]) -1, where BW is the equivalent beam width, which represents the spatial resolution and can be written as [30, 31]:
where d is the pixel size of the detector, a is the focal size of X-ray source, and M is the geometric magnification. M is defined as M = D
SD
/D
SO
, where D
SD
and D
SO
represent the distance between the source and the detector and the distance between the source and the object, respectively. The optimized geometric magnification corresponding to the best resolution is:
By taking Equation (2) into Equation (1), the best resolution BW with optimized magnification factor can be derived as:
When d ⪡ a, M| BW min ∼ 1, BW ∼ d, meaning that the resolution is determined by the detector’s pixel size. For most low-energy micro-focus X-ray CT, d ⪢ a, M| BW min ∼ ∞ , BW ∼ a, the resolution is mainly determined by the X-ray source focal size. In normal situations of common high energy X-ray CT machines, it is difficult to satisfy neither d ⪡ a nor d ⪢ a, then the spatial resolution of CT system is jointly determined by both two factors of the detector’s pixel size and the X-ray source focal size. A straightforward conclusion can be reached from Equation (3) that the smaller the X-ray focal spot size or the detector’s pixel size is, the better the spatial resolution can be obtained. For the high-energy X-ray CT system, the detector’s pixel size is difficult to reduce significantly. Microfocus with much smaller X-ray source focal size will be a good solution to reach a better spatial resolution.
As mentioned above, the development of high-energy, high-dose, microfocus X-ray sources has imposed stringent requirements on electron sources, such as high-flux, low emittance, high energy, and low energy spread. The well-developed high-average-current photo-injectors can fulfill these requirements and are suitable to drive HHM CT. Therefore, we design and implement an HHM CT on the base of the CTFEL photo-injector. The whole layout of HHM CT is present in Fig. 1, which is mainly composed of four sections, including photo-injector section, beam diagnostic section, beam matching section as well as X-ray imaging section.

Schematic diagram of the HHM CT driven by CTFEL photo-injector. Note: not to scale.
The high-average-current photo-injector composes a high-voltage photocathode DC gun [32], solenoids, a normal-conducting RF buncher and a 2×4 cell superconducting linac [33]. Electron beams are generated by illuminating a drive laser pulse on the GaAs photocathode and accelerator to 320 keV in the DC gun. Solenoid at the exit of the DC gun is for emittance compensation as well as envelope control. The buncher is to compress the pulse duration of electron bunches to avoid large energy spread in the main accelerator. The super-conducting linac is to boost beam energy to 9 MeV. Parameters of the electron beams produced by the photo-injector are listed in Table 1.
Parameters of electron beams in CTFEL
The diagnostic section uses the combination of a YAG scintillator and a CCD camera to detect the transverse profile while uses a dipole magnet to measure the energy and energy spread of the electron beams. The beam matching section is made up of three quadrupole magnets (i.e., triplet) to achieve a strong focusing of the electron beam. A removable YAG screen measures the final focus of the electron beam (on the target position) with 10μm resolution. An optimized tungsten target is to produce high-energy, high-flux X-rays via bremsstrahlung. The object to be tested is mounted on a rotating platform, 476 mm downstream of the tungsten target. A flat panel detector with 2048×2048 pixels, 0.2 mm pixel size, and DRZ + scintillators is employed to collect the transmitted X-rays. The X-ray detector is about 2000 mm downstream of the target, forming a magnification factor of 4.2. Design considerations of microfocus electron beams, efficient X-ray conversion targets, and other limiting factors are present in the following subsections.
Behavior of an electron beam focused by a strong focusing element can be analogous to the focusing of light by a strong convex lens, where larger initial beam size and shorter focal length are helpful for a smaller final beam size achievable. However, the existence of space charge forces makes focusing of electron beams more complicated. A good estimation of beams’ transverse dimension evolution after the focusing triplet can be realized by the well-known envelope equation [34]:
where, R is the envelope, I b is the beam current, I A ≈ 17KA is the Alfven current, γ is the Lorentz factor, β is the speed in the unit of light speed and ɛ _ n is the normalized transverse emittance. Physically, the second and third term in Equation (4) represent the envelope expansion induced by space charge effect and transverse emittance, respectively.
Knowing the initial value of R and R′, the above differential equation can be solved numerically. By employing thin-lens approximation, R0 is the beam size at entrance of the triplet and does not change within the triplet, while R0′ ≈ R0/f, with f being the focal length. Then we can calculate the beam envelope evolution between the triplet and the tungsten target using electron beam parameters listed in Table 1. Simulation results are shown in Fig. 2(a), where larger initial beam size and short focal length is beneficial for obtaining < 100μm electron spot size on target.

(a) Simulated beam envelope after focusing triplet with varying initial beam size and focal length. (b) Measured transverse profile of focused electron beam on target place with < 70μm lateral size (x: horizontal direction, y: vertical direction).
In practical experiments, a strong focusing triplet is located about 25 cm ahead of the tungsten target. By finely tuning the initial beam size R0 and the focusing strength of three quadrupoles, a minimum beam size of 70μm (full-width-half-maximum, FWHM) can be obtained, as shown in Fig. 2(b). Since the X-ray source dimension is partly determined by the electron beam size on the target, obtaining of < 100μm microfocus electron beam paved the way for generating microfocus X-rays.
With a specified electron beam energy, the conversion efficiency of the X-ray is determined by the thickness of the tungsten. With too thick target, the X-ray conversion efficiency will be reduced due to self-absorption of bremsstrahlung photons by the target, while the conversion efficiency is still low for too thin target, with less energy deposited by electron beam power inside the target. Therefore, the thickness of the tungsten target should be carefully optimized.
Following deductions in Ref [35], where considered energy loss of electrons and self-absorption of produced X-ray, a best thickness of the target d0 can be estimated by:

Simulated X-ray photon yield (a) and energy spectrum (b) of 9 MeV electrons with tungsten target of different thicknesses.
Also, the dose and dose rate of HHM CT should be taken into consideration. The maximum exposure at one meter position after the target can be estimated using following empirical formula [36, 37]:
The absorbed dose D in the air is:
The corresponding dose rate within one minute can be written as:
The calculation results of dose rate of 9 MeV electron beam with tungsten target are listed in
Once the microfocus electron beam has been obtained, focal spot size widening caused by scattering of electrons by the thick target needs to consider. The scattering and the focal spot size widening of the 9 MeV electron beam are simulated by CASINO [38], as shown in Fig. 4. Simulated results showed that, compared with the initial electron distribution, widening of focal spot size of 9 MeV electrons by 3 mm thick tungsten target is not obvious. Only a small portion of electrons are scattered to outer region of initial lateral size of electrons. The statistical data analysis showed that the focal spot size is increased from 100μm to about 110μm, indicating that a near 100μm X-ray source can be obtained in our HHM CT system.

Focal spot widening of electron beams by target scattering.
To investigate the effectiveness of the high-energy high-dose microfocus X-ray source described in above section, DR images of a standard X-ray resolution test phantom are acquired to evaluate the resolving power. The spatial resolution is calibrated by identifying the smallest distinguishable line pair in the original images. Next, the imaging mode is shifted to CT mode of the resolution phantom, where to measure the spatial resolution of our HHM CT system. Last, a representative test object in an industrial CT-turbine blade prototype is imaged in CT mode to evaluate the practical performance of our HHM CT. By comparing imaging results with common high energy industrial CT and low energy microfocus CT, the superiority of such an HHM CT is perfectly showcased.
DR experiment for spatial resolution calibration
Following the standard procedure in the X-ray imaging community, DR experiment of a resolution phantom is carried out to evaluate the spatial resolution of the HHM CT. Here, a test phantom made of 0.1 mm lead and containing eight groups of line-pairs is chosen. It can generate following spatial frequencies (in lp/mm): 3.4, 3.7, 4.0, 4.3, 4.6, 5.0, 5.5, 6.0, where the spatial frequency is determined by the line-pair slit width ω as f = 1/2ω. Once the DR image has been acquired, the modulation transfer function (MTF) can be measured by calculating the line profile in terms of gray-scale level along the line perpendicular to the line-pair region. The spatial resolution can be defined as the smallest achieved line-pair, which corresponds to MTF value of 10%.
In this DR experiment, using microfocus X-rays as listed in Table 2, a clear image of the resolution phantom is obtained with a total exposure time of 64 seconds (i.e. 64 frames average, one second per frame), as shown in Fig. 5(a). The measured line profile in Fig. 5(b) showed that the MTF value at 5 lp/mm is 24.6%, which is much better than the criteria of 10%. Though the MTF value drops as spatial frequency increases, this value at the highest 6 lp/mm is still 17.1%, thus indicating a lower limit of 6 lp/mm for spatial resolution in DR mode.
The calculation results of the dose rate with 3 mm tungsten target and 9 MeV electrons
The calculation results of the dose rate with 3 mm tungsten target and 9 MeV electrons

DR image of the resolution phantom: (a) original DR image; (b) line profile of the left four line-pairs in (a), indicated by the red line.
Apart from the DR experiment, a CT experiment of another resolution phantom is also carried out with the HHM CT. This resolution phantom consists of four groups of line-pairs, including 2 lp/mm, 3 lp/mm, 4 lp/mm and 5 lp/mm. The material of this resolution phantom is ferrum, while the thickness is 10 mm. The CT resolution phantom is scanned in 1000 views over 360 degree, within a whole scan time of 1000 seconds. The CT experiment results are shown in Fig. 6. As can be seen from Fig.6(a), the CT slice imaging result shows that the line-pair of 5 lp/mm is clearly distinguished. The MTF curve of the resolution phantom is shown in Fig. 6(c), in which the MTF of 5 lp/mm line-pair is as high as 30%. These results indicate that a conservative estimation of 5 lp/mm spatial resolution can be achieved in the HHM CT, in well agreement with the DR imaging results. To our best knowledge, this is the first time that spatial resolution of 100μm has been demonstrated when using a MeV high-energy, high-dose, microfocuse X-ray CT.

CT image of the resolution phantom. (a) the CT slice imaging result; (b) the CT 3D result; (c) the gray distribution curve of the CT image and the measured MTF.
Finally, the practical performance of the HHM CT is checked with the CT experiment of a turbine blade prototype [39]. The main system parameters of this HHM CT, 450 kV CT and conventional high energy CT are listed in Table 3. The turbine blade is scanned in 1000 views over 360 degrees, with a scanning period of only 16.67 minutes. Reconstruction of the turbine blade is realized using the Feldkamp-Davis-Kress (FDK) [40] algorithm. Here, the resolving power of the HHM CT on this turbine blade is calibrated by measuring the edge response function (ERF) and the corresponding point spread function (PSF) of the inner arc region of the blade, which is a key part that determines the image quality of the turbine blade prototype, as indicated by the red short line in Fig.7(a). With this HHM CT, boundaries of the blade are clearly identified in CT image, while the calibrated PSF is 261μm. In comparison, though the CT image of the blade using a 450 kV microfocus X-ray CT is measured to hold a spatial resolution of 364μm, the blade cannot be fully penetrated, leaving certain parts indistinguishable, see Fig. 7(c). Moreover, CT image of the turbine blade is also obtained with a conventional high-energy industrial X-ray CT, as shown in Fig. 7(e). Obviously, the boundary of the blade is severely blurred due to the large focal spot size of the X-ray, and the calculated PSF indicates a much poorer resolution of 857μm. In this sense, considerable advantages, including high spatial resolution as well as high penetration power provided by this HHM CT can be taken of, to significantly promote the detection capabilities of such sophisticated workpieces.
The main parameters of HHM CT, 450 kV CT and conventional high energy CT
The main parameters of HHM CT, 450 kV CT and conventional high energy CT

Turbine blade CT images: (a) CT image of the HHM CT, (b) the edge response function (ERF) (blue line) and the PSF (orange line) of the HHM CT image in the red line, (c) image of the 450 kV CT, (d) the ERF (blue line) and the PSF (orange line) of the 450 kV CT image in the red line, (e) image of the common high energy CT, (f) the ERF (blue line) and the PSF (orange line) of the common high energy CT image in the red line.
This study provides a comprehensive design, implement and experimental verification of a HHM CT prototype on the base of CTFEL high-average-current photo-injector. Realization of microfocus X-ray source is achieved by a strongly focused < 70μm electron beam spot on the tungsten target. The target thickness is carefully optimized, with an overall consideration of photon yielding, energy spectrum and spot size widening effect. Experimental investigations of resolution phantoms are performed with the HHM CT, in both DR and CT imaging mode. Imaging results of resolution phantoms show that a DR spatial resolution of up to 6 lp/mm and a CT spatial resolution better than 5 lp/mm can be readily obtained. Moreover, the practical performance of the HHM CT is checked with the imaging of a turbine blade prototype, where the blade is fully penetrated and clearly imaged with a quite good resolution of 261μm. By comparing imaging results with conventional low energy microfocus X-ray CT and high energy industrial X-ray CT, the superiority of such an HHM CT with simultaneously obtainable high resolution power and high penetration power is well established. Success of this HHM CT prototype should pave the way to standard HHM CT machine in the future.
Conclusion
In summary, this study paves the way to the development of HHM CT driven by a high-average-current photo-injector. This technology may become an important means of high-energy metrology CT. More importantly, successful imaging of the turbine blade using this HHM CT should have a major impact on non-destructive testing and non-destructive evaluation of high-density fine-structured components, such as those used in aviation, aerospace, ships, automobiles, etc.
Footnotes
Acknowledgments
This work is supported by the National Natural Science Foundation of China (NSFC Grants No.11975218, 11905210, 11805192, 12005211 and 11605190), and by the Innovation Foundation of CAEP with grant (No.CX2019036, CX2019037).
Data availability
The data that support the findings of this study are available from the corresponding authors upon reasonable request.
