# Characterization of Energetic Protons Generated in the ShenGuang-II UP Petawatt Laser Interactions with Foil Targets

Authors: Huiya Liu; Anle Lei; Ning Kang; Honghai An; Zhiyong Xie; Yao Zhao; Shenlei Zhou; Mingying Sun; Baoqiang Zhu; Wei Wang; Jianqiang Zhu
Journal: Laser and Particle Beams
 (2021)

---

## Abstract

The characterization of energetic protons generated in the ShenGuang-II UP petawatt laser interactions with foil targets has been systematically studied. The proton energy spectra and angular distributions are measured with a radiochromic film stack. It shows that the proton energy spectra have a Boltzmann distribution with temperature of about 2.8 MeV and cutoff energy of about 20 MeV. The divergence angles of protons vary from 10° to 60°, dependent on the proton energy. The proton source size and location are investigated via the proton point-projection mesh imaging. The proton virtual sources are found to locate tens to hundreds of microns in front of the foil target, depending on the proton energies. A Monte Carlo simulation estimates the diameter of the virtual proton source to be about 12μm for the protons with energy of 16.8 MeV, which is much smaller than the laser focus size of about 50 μm. The spatial resolution of the 16.8 MeV proton imaging is quantified with the point spread function to be about 15 μm, which is consistent with the proton virtual source size. These results will be important for the users conducting experiments with the protons as a backlighting source on the ShenGuang-II UP petawatt laser.


---

## Body

## 1. Introduction

The energetic protons generated in short, intense laser interactions with foil targets have received much attention in recent years and have been widely studied both experimentally and theoretically [1–5]. The protons accelerated by the electric field at the rear side of the target via the target normal sheath acceleration (TNSA) mechanism have properties including short pulse duration and high brightness [6]. The protons have many potential applications, such as probing the fast dynamics of electric or magnetic fields in the laser produced plasmas [7–10], as a fast ignitor beam for laser-driven fusion [11], proton radiobiology [12–16], material science [17], and neutron production [18]. The characteristics of the laser-driven protons such as energy spectrum, divergence angle, and source size are essential for these applications.The energetic protons accelerated through TNSA are from the water or hydrocarbon contaminant layers on the target rear surface and can be accelerated to tens of MeV. Usually, in the point-projection imaging experiments, the size of the source decides the best resolution [19]. However, the transverse size of the proton source was measured to be much larger than the imaging resolution [20]. It was, therefore, considered that the protons were emitted in a quasilaminar fashion from a virtual source with a small size [11, 21].ShenGuang-II UP laser facility at the National Laboratory on High-Power Laser and Physics is a multifunctional experimental platform, which includes eight kilojoule-class nanosecond laser beams and one picosecond petawatt laser beam (called the SG-II UP petawatt laser, i.e., the SG-II UP PW laser) [22]. The SG-II UP PW laser is usually used for fast heating of an imploded dense core plasma or generating energetic protons as a backlight source to investigate the plasmas produced by the nanosecond lasers. In this article, the characterization of energetic protons generated in the SG-II UP PW laser interactions with foil targets was systematically presented. The proton point-projection imaging with Cu meshes in combination with a radiochromic film (RCF) stack was used to characterize the protons. The proton energy spectra, angular distribution, and the virtual source size and location were obtained. The spatial resolution of the proton imaging was also analysed and compared to the virtual source size. These results and findings will be important for the potential users to know the laser status and to understand the properties of the protons so as to design their experiments.


## 2. Experimental Setup

The experiment was conducted on the SG-II UP PW laser. The laser uses a hybrid technology combining optical parametric chirped-pulse amplifier and Nd:glass chirped-pulse amplifier to generate short pulse (1–10 ps) and large energy (up to 1000 J) at 1053 nm [22]. The experimental setup is shown in Figure 1. The laser was focused with an f/2.5 off-axis parabolic mirror to the 10 μm thickness Au foil targets, at an incidence angle of 21°. In some shots, a Cu mesh was put behind the foil target at a distance d=2 mm, as the object in the proton point-projection imaging to judge the quality of the proton imaging. The meshes (AG 200, Gilder Grids) were formed by 36 μm lines with 90 μm spacing. The protons emitting from the rear surface of the target passed through the mesh and were detected with RCF (HD-V2, GafChromic) stacks [23]. The front surface of the RCF stack was located L=32 mm behind the Cu mesh.Figure 1
Experimental setup.An X-ray pinhole camera (XPHC) with the theoretical resolution of about 13μm was used to monitor the laser focus spot [24, 25]. An imaging plate (BAS-SR, Fuji Film) was used as the X-ray detector. In the experiment, the incident laser energy and pulse duration were set at 130 J and 1 ps, respectively, which represents a routine status of the laser as a proton backlighter. Figure 2(a) shows a typical XPHC image of the laser focus. Because a lot of physical processes are involved in the X-ray emission, such as fast electron reflux, recirculation, and transport inside the target (collision, Ohmic heating, etc.), scientifically scaling the measured emitted X-ray profile with the real laser focus intensity on the target is difficult. In this article, we assume that the X-ray intensity obtained by XPHC scales linearly with the laser intensity. The cumulative laser energy fraction as functions of intensity and focus diameter can be obtained, as shown in Figure 2(b). It shows that about 50% of the laser energy is contained within a 60 μm diameter, resulting in laser intensity over 1.5 × 1018 W/cm2. The experimental results measured in each shot are shown in Figure 2(c), and the shot number 1 is chosen as a typical shot. The incident laser energies are stable at 130 J ± 15%, and the proton cutoff energies measured with the RCF stacks are all higher than 15 MeV.Figure 2
(a) A typical laser focus image recorded by an X-ray pinhole camera. PSL is the unit of the signal intensity on the imaging plate. (b) Curves of the cumulative laser energy fraction as functions of laser intensity and focus diameter obtained from (a). (c) Experimental results measured in each shot. The shot number 1 is chosen as the typical shot (data in the green box).
(a)(b)(c)


## 3. Experimental Results and Analysis

### 3.1. Proton Energy Spectra and Angular Distribution

The generated protons pass through the RCF layers, depositing their kinetic energy in the RCFs [26]. Typical proton images recorded on the RCFs scanned by using the EPSON V750 scanner for the same shot as in Figure 2(a) are shown in Figure 3(a). It is obvious that the lower-energy protons have higher intensity and larger angular divergence angle. The maximum proton energy reaches up to 18.9 MeV.Figure 3
(a) Proton images on RCF layers for the same shot as in Figure2(a). (b) Proton energy spectrum derived from (a), which was fitted with Boltzmann distribution. (c) Divergence angle of the protons. The results from VULCAN and 100TW-LULI lasers by Nürnberg et al. [26] are used for comparison.
(a)(b)(c)The energy spectrum of the protons can be obtained from the RCF images. Based on the calibration of the HD-V2-type RCF for protons with different doses in [23, 26, 27], the optical density of the RCFs can be converted into the deposited proton energy. The spatially integrated proton number on each RCF layer can, thus, be obtained from the deposited energy with SRIM simulations and the spectrum unfolding calculation [28, 29], as shown in Figure 3(b). The proton energy spectrum shows an exponential profile. We fit the spectrum with the Boltzmann distribution. The fitted energy spectrum is shown in Figure 3(b), with total proton number and the temperature kBT of about 1.2 × 1012 and 2.8 MeV, respectively. As we can see from the Figure 3(b), because of the saturation of the lowest-energy RCF layers, the measured proton numbers at low proton energies are less than the fitted result. Based on the fitting result, the conversion efficiency of laser energy to the proton beam can be calculated to be about 0.62%.With the given distances between the target and the RCF layers and sizes of the proton images, the divergence angles (referring the full beam aperture in this article) of the energy-resolved protons can be calculated, as shown in Figure3(c). The results from VULCAN and 100TW-LULI lasers by Nürnberg et al. [26] are also shown for comparison. In Figure 3(c), the proton energy E is scaled to the maximum proton energy Emax, here Emax=18.9 MeV, 29.7 MeV, and 16.2 MeV for the SG-II UP PW, VULCAN, and 100TW-LULI lasers, respectively. The divergence angles of protons in our experiment are from 10° to 60°, dependent on the proton energy, with the divergence angle decreasing with the increase in the proton energy. The divergence angles of energy-resolved protons generated by the SG-II UP PW laser and VULCAN laser are similar but relatively differ from the 100TW-LULI laser. Note that the laser energy for the SG-II UP PW and VULCAN shot is about 130 J and that for the 100TW-LULI shot is about 15 J. That may be the reason why their divergences of protons are different.


### 3.2. Proton Source Location and Size

It has been shown that the laser-driven protons are emitted from a virtual source in a quasilaminar fashion [20]. The location of the virtual source can be obtained through calculating the magnification rate M of the Cu mesh in the proton image. The magnification rate M of the point-projection imaging in our experiment is the ratio of the virtual source-to-RCF layer distance to the virtual source-to-mesh distance.(1)M=ν+d+L+lν+d,where ν is the distance from the virtual source to the foil target and l is the distance from the front surface of the RCF stack to a RCF layer with specific proton energy. In the experiment, the magnification rate M can be measured as M=dRCF/dmesh, where dRCF is the period of the mesh in the proton image and dmesh is the real period of the mesh. For the image of protons with energy of 16.8 MeV shown in Figure 3(a), dRCF=2.131 mm, dmesh=0.126 mm, l=2.2 mm, d=2 mm, and L=32 mm, we can then obtain ν=0.151 mm. Using the same method, the positions of the virtual sources for protons with different energies in front of the target are obtained with the uncertainty of about 7%–15%, as listed in Table 1.Table 1
Energy-resolved position of the virtual source in front of the targetν and virtual source size Svirtual (for VULCAN and 100TW-LULI lasers) or 2σ of ESF value (for SG-II UP PW) for three different laser systems. The results from VULCAN and 100TW-LULI lasers by Nürnberg et al. [26] are used for comparison.
Laser systemE (MeV)ν (μm)Svirtual or 2σ of ESF (μm)SG-II UP PW6.441819.27.939222.99.138113.710.630612.412.026810.213.323110.715.217911.016.815113.818.98917.9VULCAN6.237212.117.417311.0100TW-LULI4.7755.49.8311.413.3310.4To estimate the size of the virtual source, a series of Monte Carlo simulations based on Geant4 were carried out for the protons with various source sizes [11, 30]. In the simulations, spatial Gaussian density distribution in the proton source was used.(2)nx,y=n0exp−x2r2−y2r2,where n0 and r are the total proton number and the 1/e is the radius of the proton source, respectively. The spatially resolved optical density distribution in the RCF layer is calculated, as shown in Figure 4. By comparing the experimental optical density profile with the simulated ones, we find that the protons with a virtual source size of about 12 μm in 1/e diameter can well match the experimental optical density profile.Figure 4
(a) The optical density distribution of 16.8 MeV protons in Figure3(a). (b) The experimental optical density profile of 16.8 MeV protons in the white box region in (a), and the simulated optical density profile of 16.8 MeV protons with source sizes of 1 μm, 12 μm, and 50 μm in 1/e diameter, respectively.
(a)(b)


### 3.3. Spatial Resolution of the Proton Imaging

The spatial resolution of the proton imaging reflects the size of the proton virtual source. We, therefore, analysed the mesh images on RCFs to obtain the point spread function (PSF) so as to quantify the proton imaging spatial resolution. We first fitted the edge of the experimental optical density profile in Figure4(b) for X from 11 mm to 14 mm with an edge spread function (ESF) [31], which is(3)IESFx=I0⋅erfxσ+a0+a1x,where(4)erfxσ=2π∫0x/σe−t2dt,and I0, σ, a0, and a1 are fitting coefficients. The fitted ESF is shown in Figure 5(a), which gives σ=6.9 μm. Considering that the ESF is the integral response function of a Gaussian-type PSF, the PSF can be obtained as(5)IPSFx=I02πσexp−x2σ2+a1.Figure 5
(a) The experimental optical density profile in Figure4(b) for X from 11 mm to 14 mm was fitted by an edge spread function (ESF) with the coefficient2σ=13.8. M=dRCF/dmesh=16.9M=dRCF/dmesh=16.9 is the magnification rate of the imaging for the protons with energy of 16.8 MeV. (b) The MTF for the protons with different energies; the shaded area represents the range of the possible MTF.
(a)(b)So the full width at 1/e maximum of spatial resolution of proton imaging is 2σ=13.8μm. This is consistent with the size of the proton virtual source of about 12 μm obtained through Monte Carlo simulations. The energy-resolved 2σ of ESF values are listed in Table 1 and are compared with the virtual source sizes from VULCAN and 100TW-LULI lasers by Nürnberg et al. [26].In addition, the normalized Fourier transform of the PSF represents the modulation transfer function (MTF), which is(6)IMTFk=exp−k2σ24,where k=2π/λ and λ is the period of the modulation. Figure 5(b) shows the MTF for the protons with different energies. One sees that the spatial resolution of the 16.8 MeV proton imaging at 20% MTF is about 15 μm and that for other protons ranges from 12 to 25 μm.


## 3.1. Proton Energy Spectra and Angular Distribution

The generated protons pass through the RCF layers, depositing their kinetic energy in the RCFs [26]. Typical proton images recorded on the RCFs scanned by using the EPSON V750 scanner for the same shot as in Figure 2(a) are shown in Figure 3(a). It is obvious that the lower-energy protons have higher intensity and larger angular divergence angle. The maximum proton energy reaches up to 18.9 MeV.Figure 3
(a) Proton images on RCF layers for the same shot as in Figure2(a). (b) Proton energy spectrum derived from (a), which was fitted with Boltzmann distribution. (c) Divergence angle of the protons. The results from VULCAN and 100TW-LULI lasers by Nürnberg et al. [26] are used for comparison.
(a)(b)(c)The energy spectrum of the protons can be obtained from the RCF images. Based on the calibration of the HD-V2-type RCF for protons with different doses in [23, 26, 27], the optical density of the RCFs can be converted into the deposited proton energy. The spatially integrated proton number on each RCF layer can, thus, be obtained from the deposited energy with SRIM simulations and the spectrum unfolding calculation [28, 29], as shown in Figure 3(b). The proton energy spectrum shows an exponential profile. We fit the spectrum with the Boltzmann distribution. The fitted energy spectrum is shown in Figure 3(b), with total proton number and the temperature kBT of about 1.2 × 1012 and 2.8 MeV, respectively. As we can see from the Figure 3(b), because of the saturation of the lowest-energy RCF layers, the measured proton numbers at low proton energies are less than the fitted result. Based on the fitting result, the conversion efficiency of laser energy to the proton beam can be calculated to be about 0.62%.With the given distances between the target and the RCF layers and sizes of the proton images, the divergence angles (referring the full beam aperture in this article) of the energy-resolved protons can be calculated, as shown in Figure3(c). The results from VULCAN and 100TW-LULI lasers by Nürnberg et al. [26] are also shown for comparison. In Figure 3(c), the proton energy E is scaled to the maximum proton energy Emax, here Emax=18.9 MeV, 29.7 MeV, and 16.2 MeV for the SG-II UP PW, VULCAN, and 100TW-LULI lasers, respectively. The divergence angles of protons in our experiment are from 10° to 60°, dependent on the proton energy, with the divergence angle decreasing with the increase in the proton energy. The divergence angles of energy-resolved protons generated by the SG-II UP PW laser and VULCAN laser are similar but relatively differ from the 100TW-LULI laser. Note that the laser energy for the SG-II UP PW and VULCAN shot is about 130 J and that for the 100TW-LULI shot is about 15 J. That may be the reason why their divergences of protons are different.


## 3.2. Proton Source Location and Size

It has been shown that the laser-driven protons are emitted from a virtual source in a quasilaminar fashion [20]. The location of the virtual source can be obtained through calculating the magnification rate M of the Cu mesh in the proton image. The magnification rate M of the point-projection imaging in our experiment is the ratio of the virtual source-to-RCF layer distance to the virtual source-to-mesh distance.(1)M=ν+d+L+lν+d,where ν is the distance from the virtual source to the foil target and l is the distance from the front surface of the RCF stack to a RCF layer with specific proton energy. In the experiment, the magnification rate M can be measured as M=dRCF/dmesh, where dRCF is the period of the mesh in the proton image and dmesh is the real period of the mesh. For the image of protons with energy of 16.8 MeV shown in Figure 3(a), dRCF=2.131 mm, dmesh=0.126 mm, l=2.2 mm, d=2 mm, and L=32 mm, we can then obtain ν=0.151 mm. Using the same method, the positions of the virtual sources for protons with different energies in front of the target are obtained with the uncertainty of about 7%–15%, as listed in Table 1.Table 1
Energy-resolved position of the virtual source in front of the targetν and virtual source size Svirtual (for VULCAN and 100TW-LULI lasers) or 2σ of ESF value (for SG-II UP PW) for three different laser systems. The results from VULCAN and 100TW-LULI lasers by Nürnberg et al. [26] are used for comparison.
Laser systemE (MeV)ν (μm)Svirtual or 2σ of ESF (μm)SG-II UP PW6.441819.27.939222.99.138113.710.630612.412.026810.213.323110.715.217911.016.815113.818.98917.9VULCAN6.237212.117.417311.0100TW-LULI4.7755.49.8311.413.3310.4To estimate the size of the virtual source, a series of Monte Carlo simulations based on Geant4 were carried out for the protons with various source sizes [11, 30]. In the simulations, spatial Gaussian density distribution in the proton source was used.(2)nx,y=n0exp−x2r2−y2r2,where n0 and r are the total proton number and the 1/e is the radius of the proton source, respectively. The spatially resolved optical density distribution in the RCF layer is calculated, as shown in Figure 4. By comparing the experimental optical density profile with the simulated ones, we find that the protons with a virtual source size of about 12 μm in 1/e diameter can well match the experimental optical density profile.Figure 4
(a) The optical density distribution of 16.8 MeV protons in Figure3(a). (b) The experimental optical density profile of 16.8 MeV protons in the white box region in (a), and the simulated optical density profile of 16.8 MeV protons with source sizes of 1 μm, 12 μm, and 50 μm in 1/e diameter, respectively.
(a)(b)


## 3.3. Spatial Resolution of the Proton Imaging

The spatial resolution of the proton imaging reflects the size of the proton virtual source. We, therefore, analysed the mesh images on RCFs to obtain the point spread function (PSF) so as to quantify the proton imaging spatial resolution. We first fitted the edge of the experimental optical density profile in Figure4(b) for X from 11 mm to 14 mm with an edge spread function (ESF) [31], which is(3)IESFx=I0⋅erfxσ+a0+a1x,where(4)erfxσ=2π∫0x/σe−t2dt,and I0, σ, a0, and a1 are fitting coefficients. The fitted ESF is shown in Figure 5(a), which gives σ=6.9 μm. Considering that the ESF is the integral response function of a Gaussian-type PSF, the PSF can be obtained as(5)IPSFx=I02πσexp−x2σ2+a1.Figure 5
(a) The experimental optical density profile in Figure4(b) for X from 11 mm to 14 mm was fitted by an edge spread function (ESF) with the coefficient2σ=13.8. M=dRCF/dmesh=16.9M=dRCF/dmesh=16.9 is the magnification rate of the imaging for the protons with energy of 16.8 MeV. (b) The MTF for the protons with different energies; the shaded area represents the range of the possible MTF.
(a)(b)So the full width at 1/e maximum of spatial resolution of proton imaging is 2σ=13.8μm. This is consistent with the size of the proton virtual source of about 12 μm obtained through Monte Carlo simulations. The energy-resolved 2σ of ESF values are listed in Table 1 and are compared with the virtual source sizes from VULCAN and 100TW-LULI lasers by Nürnberg et al. [26].In addition, the normalized Fourier transform of the PSF represents the modulation transfer function (MTF), which is(6)IMTFk=exp−k2σ24,where k=2π/λ and λ is the period of the modulation. Figure 5(b) shows the MTF for the protons with different energies. One sees that the spatial resolution of the 16.8 MeV proton imaging at 20% MTF is about 15 μm and that for other protons ranges from 12 to 25 μm.


## 4. Conclusions

In summary, the characteristics of the energetic protons accelerated from the SG-II UP PW laser interactions with foil targets were investigated. The maximum proton energy obtained in the experiment was 18.9 MeV with a laser pulse of 130 J and 1 ps. The energy spectrum of protons measured by the RCF stack is fitted with the Boltzmann distribution with the temperature of about 2.8 MeV. The emission of the protons showed a decreasing divergence angle from 60° to 10° with the increase in the proton energy. Based on the proton point-projection mesh imaging, the proton virtual sources for protons with different energies were found to locate tens to hundreds of microns in front of the target, and the size was estimated with a Monte Carlo simulation, to be about 12μm in 1/e diameter for the protons with energy of 16.8 MeV. This size is consistent with the proton imaging spatial resolution of about 15 μm obtained from the point spread function. These results are of importance for experiments using the protons as a backlighting source on the SG-II UP PW laser.



---

*Source: 7205383-2021-08-18.xml*