Hostname: page-component-7c8c6479df-xxrs7 Total loading time: 0 Render date: 2024-03-28T10:45:36.584Z Has data issue: false hasContentIssue false

The application of proton spectrometers at the SG-III facility for ICF implosion areal density diagnostics

Published online by Cambridge University Press:  29 October 2015

Xing Zhang
Affiliation:
Laser Fusion Research Center, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900, China
Jianhua Zheng
Affiliation:
Laser Fusion Research Center, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900, China
Ji Yan
Affiliation:
Laser Fusion Research Center, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900, China
Zhenghua Yang
Affiliation:
Laser Fusion Research Center, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900, China
Ming Su
Affiliation:
Laser Fusion Research Center, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900, China
Yudong Pu
Affiliation:
Laser Fusion Research Center, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900, China
Pin Yang
Affiliation:
Laser Fusion Research Center, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900, China
Xufei Xie
Affiliation:
Laser Fusion Research Center, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900, China
Li Chen
Affiliation:
Laser Fusion Research Center, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900, China
Ming Chen
Affiliation:
Laser Fusion Research Center, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900, China
Tianxuan Huang
Affiliation:
Laser Fusion Research Center, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900, China
Shao’en Jiang
Affiliation:
Laser Fusion Research Center, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900, China
Shenye Liu*
Affiliation:
Laser Fusion Research Center, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900, China
Jiamin Yang
Affiliation:
Laser Fusion Research Center, Chinese Academy of Engineering Physics, Mianyang, Sichuan 621900, China
*
Correspondence to:  S. Liu, PO Box 919-986, No. 64 Mianshan Road, Mianyang, Sichuan province 621900, China. Email: lsye1029@163.com

Abstract

Charged particle diagnostics is one of the required techniques for implosion areal density diagnostics at the SG-III facility. Several proton spectrometers are under development, and some preliminary areal density diagnostics have been carried out. The response of the key detector, CR39, to charged particles was investigated in detail. A new track profile simulation code based on a semi-empirical model was developed. The energy response of the CR39 detector was calibrated with the accelerator protons and alphas from a 241Am source. A proton spectrometer based on the filtered CR39 detector was developed, and D–D primary proton measurements were implemented. A step range filter spectrometer was developed, and preliminary areal density diagnostics was carried out. A wedged range filter spectrometer array made of Si with a higher resolution was designed and developed at the SG-III facility. A particle response simulation code by the Monte Carlo method and a spectra unfolding code were developed. The capability was evaluated in detail by simulations.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
The online version of this article is published within an Open Access environment subject to the conditions of the Creative Commons Attribution licence .
Copyright
© The Author(s) 2015

1. Introduction

The SG-III facility has 48 laser beams with a maximum output energy of 180 kJ. It is a new platform for inertial confinement fusion (ICF) experiments. A high enough areal density (${\it\rho}$R) is necessary for ignition, and how to realize this is one of the most important issues in implosion physics research. Measurement of the energy downshift of the fusion protons is an effective technique in passive ${\it\rho}$R diagnostics. In D–D fuel implosion, protons are emitted in a primary D–D reaction and a secondary D–3He reaction with energies of 3 and 12–17 MeV, respectively. In D–3He fuel implosion, protons from primary D–3He, D–D and 3He–3He reactions are dominant. When the protons, generated in the core, pass through the fuel and shell plasmas, there will be an energy loss, which is related to ${\it\rho}$R, due to the interaction with the plasma. The stopping power of the protons in dense plasma, quite different from that in non-ionized matter, can be calculated with some well-defined models, such as the Li–Petrasso (L–P) model[Reference Li and Petrasso1] and the Brown–Preston–Singleton (BPS) model[Reference Browna, Preston and Singleton2]. This technique has been widely used in ${\it\rho}$R diagnostics at NIF[Reference Zylstra, Frenje, Séguin, Rosenberg, Rinderknecht, Gatu Johnson, Casey, Sinenian, Manuel, Waugh, Sio, Li, Petrasso, Friedrich, Knittel, Bionta, McKernan, Callahan, Collins, Dewald, Döppner, Edwards, Glenzer, Hicks, Landen, London, Mackinnon, Meezan, Prasad, Ralph, Richardson, Rygg, Sepke, Weber, Zacharias, Moses, Kilkenny, Nikroo, Sangster, Glebov, Stoeckl, Olson, Leeper, Kline, Kyrala and Wilson3], OMEGA[Reference Séguin, Li, Frenje, Hicks, Green, Kurebayashi, Petrasso, Soures, Meyerhofer, Glebov, Radha, Stoeckl, Roberts, Sorce, Sangster, Cable, Fletcher and Padalino4], GEKKO-XII[Reference Kitagawa, Tanaka, Nakai, Yamanaka, Nishihara, Azechi, Miyanaga, Norimatsu, Kanabe, Chen, Richard, Sato, Furukawa and Nakai5] and some other ICF facilities. Moreover, the ${\it\rho}$R evolution in D–3He implosion can be diagnosed by the different energy downshifts at the shock bang time and compression bang time[Reference Frenje, Li, Séguin, Deciantis, Kurebayashi, Rygg, Petrasso, Delettrez, Glebov, Stoeckl, Marshall, Meyerhofer, Sangster, Smalyuk and Soures6, Reference Li, Séguin, Frenje, Rosenberg, Petrasso, Amendt, Koch, Landen, Park, Robey, Town, Casner, Philippe, Betti, Knauer, Meyerhofer, Back, Kilkenny and Nikroo7]. Besides, the fuel areal density in D–D implosion can be diagnosed by the yield ratio of secondary protons and primary neutrons, and the residual CH ablator areal density in D–T implosion can be diagnosed by the neutron recoil proton spectra. Through ${\it\rho}$R measurements at multiple view angles by several proton spectrometers, the compression asymmetry can be diagnosed[Reference Séguin, Li, Frenje, Hicks, Green, Kurebayashi, Petrasso, Soures, Meyerhofer, Glebov, Radha, Stoeckl, Roberts, Sorce, Sangster, Cable, Fletcher and Padalino4]. Furthermore, the mix effect can be studied by fusion proton measurement in some special targets[Reference Li, Séguin, Frenje, Rosenberg, Petrasso, Amendt, Koch, Landen, Park, Robey, Town, Casner, Philippe, Betti, Knauer, Meyerhofer, Back, Kilkenny and Nikroo7Reference Casey, Woods, Smalyuk, Hurricane, Glebov, Stoeckl, Theobald, Wallace, Nikroo, Schoff, Wu, Frenje, Landen, Remington and Glendinning10]. The measured 3He–3He proton spectra provide a new approach for laboratory astrophysics research.

Several kinds of proton spectrometer have been developed at NIF, OMEGA, NOVA, GEKKO-XII and some other ICF facilities, such as the step range filter (SRF) spectrometer[Reference Séguin, Frenje, Li, Hicks, Kurebayashi, Rygg, Schwartz, Petrasso, Roberts, Soures, Meyerhofer, Sangster, Knauer, Sorce, Glebov, Stoeckl, Phillips, Leeper, Fletcher and Padalino11, Reference Rosenberg, Zylstra, Frenje, Rinderknecht, Gatu Johnson, Waugh, Séguin, Sio, Sinenian, Li, Petrasso, Glebov, Hohenberger, Stoeckl, Sangster, Yeamans, LePape, Mackinnon, Bionta, Talison, Casey, Landen, Moran, Zacharias, Kilkenny and Nikroo12], the wedged range filter (WRF) spectrometer[Reference Séguin, Frenje, Li, Hicks, Kurebayashi, Rygg, Schwartz, Petrasso, Roberts, Soures, Meyerhofer, Sangster, Knauer, Sorce, Glebov, Stoeckl, Phillips, Leeper, Fletcher and Padalino11, Reference Séguin, Sinenian, Rosenberg, Zylstra, Manuel, Sio, Waugh, Rinderknecht, Gatu Johnson, Frenje, Li, Petrasso, Sangster and Roberts13] and the magnet-based spectrometer[Reference Kitagawa, Tanaka, Nakai, Yamanaka, Nishihara, Azechi, Miyanaga, Norimatsu, Kanabe, Chen, Richard, Sato, Furukawa and Nakai5, Reference Séguin, Frenje, Li, Hicks, Kurebayashi, Rygg, Schwartz, Petrasso, Roberts, Soures, Meyerhofer, Sangster, Knauer, Sorce, Glebov, Stoeckl, Phillips, Leeper, Fletcher and Padalino11]. The first two are compact spectrometers with high detection efficiencies and perform well in D–3He reaction proton measurement. The last one has the highest energy resolution, while the space occupation and the cost are higher. For all of the spectrometers, the CR39 detector is the key detector because of the high sensitivity to proton, alpha and some other ionizing particles and quite low sensitivity to electron, x-ray and ${\it\gamma}$-ray. In this work, the development of the proton spectrometers for ${\it\rho}$R diagnostics at the SG-III facility is introduced. The detailed charged particle response of the CR39 detector is studied in Section 2. A new track profile simulation code is presented, as well as an experimental study of the energy response. A proton spectrometer based on the filtered CR39 detector and the D–D primary proton measurements are described in Section 3. Section 4 describes the setup of a SRF spectrometer and preliminary ${\it\rho}$R diagnostics. Section 5 describes the setup of a Si-WRF spectrometer array at the SG-III facility and its capability evaluation. Some concluding remarks are listed in Section 6.

2. Energy response of the CR39 detector to charged particles

2.1. Processing of the the CR39 detector

The CR39 is a kind of solid-state nuclear track detector. It is a transparent plastic with a chemical composition of $\text{C}_{12}\text{H}_{18}\text{O}_{7}$. When a charged particle passes through a CR39 detector, it deposits energy along the trail, mostly by ionization. A narrow damaged trail in the CR39, on the scale of nm, is formed. The amount of damage along the trail is mainly determined by the local electronic stopping power of the incident particle[Reference Fleischer, Price and Walker14]. When the irradiated CR39 is exposed in an etchant, the surface is etched at a bulk velocity $v_{\text{B}}$, while the damaged trail is etched at a track velocity $v_{\text{T}}$. A visible track, on the scale of ${\rm\mu}\text{m}$, is formed after etching. The appearance parameters of the tracks are related to the species and energy of the incident particles. The CR39 detector has been widely used in many fields, such as ICF proton diagnostics, dosimetry for dose assessment and particle detection in space[Reference Séguin, Frenje, Li, Hicks, Kurebayashi, Rygg, Schwartz, Petrasso, Roberts, Soures, Meyerhofer, Sangster, Knauer, Sorce, Glebov, Stoeckl, Phillips, Leeper, Fletcher and Padalino11, Reference Caresana, Ferrarini, Fuerstner and Mayer15Reference Zhou, Semones, Guetersloh, Zapp, Weyland and Benton17].

In this work, the etchant is $6~\text{mol}/\text{L}$ NaOH solution at a temperature of $80\,^{\circ }\text{C}$. The typical etching time is 6 h, and it can be set in a range of 2–9 h to make sure that the tracks are visible and there is no saturation. After etching, the CR39 detectors are rinsed with dilute acetic acid and highly purified water, and then dried. The appearance parameters of the particle track profiles are obtained from photos of the etched detector’s surfaces, by a digital microscope with objective magnifications of $20\times$ and $50\times$. The spatial resolution of the track diameter is approximately 1 and $0.4~{\rm\mu}\text{m}$ for the two magnifications, respectively.

2.2. Simulation of the track profile induced by a charged particle

Three track parameters are mostly used in the particle analysis: diameter, contrast and eccentricity. The diameter is always used to deduce the particle energy and species, while the contrast and eccentricity are used to reject the intrinsic tracks and neutron induced tracks[Reference Séguin, Frenje, Li, Hicks, Kurebayashi, Rygg, Schwartz, Petrasso, Roberts, Soures, Meyerhofer, Sangster, Knauer, Sorce, Glebov, Stoeckl, Phillips, Leeper, Fletcher and Padalino11]. The track diameter has a nonlinear response to the particle energy, and it varies with the etching conditions. In the measurement with a specific etching condition (etchant, temperature and time), a well-defined empirical $V$-function ($V=v_{\text{T}}/v_{\text{B}}$) based etching kinematics model can describe the response well, and it has been studied in many works[Reference Caresana, Ferrarini, Fuerstner and Mayer15, Reference Zhou, Semones, Guetersloh, Zapp, Weyland and Benton17, Reference Nikezic and Yu18]. In the application of ICF proton measurement, the etching time varies for different fusion yields and some other conditions. A semi-empirical model[Reference Hicks19] of the $V$-function, independent of the etching time, has been proposed, as shown in Equation (1), where $\left(\frac{dE}{dx}\right)_{elec}$ is the local electronic stopping power along the trail, and $k$ and $n$ are the free parameters. In this section, a new etching kinematics simulation code based on a semi-empirical model is developed. The electronic stopping power is calculated by the TRIM code[Reference Ziegler, Ziegler and Biersack20]. The track profiles of perpendicular incident particles and oblique incident particles are simulated.

(1)$$\begin{eqnarray}V=\frac{V_{\text{T}}}{V_{\text{B}}}=1+k\left(\frac{dE}{dx}\right)_{elec}^{n}.\end{eqnarray}$$

In the case of perpendicular incidence, the track profile can be calculated in two steps, as shown in Figure 1(a). The damaged trail is first etched to the point $A^{\prime }(y^{\prime },z^{\prime })$ in a time of $t^{\prime }$ at a velocity of $v_{\text{T}}(z)$. Then, the wall of the nano-hole is etched at a velocity of $v_{\text{B}}$ in the direction determined by the angle of ${\it\xi}$. In the time of ($t{-}t^{\prime }$), the track profile reaches the point $A(y_{A},z_{A})$. The two-step process can be described by Equations (2)–(4). The bulk etching depth can be calculated by Equation (5). With Equations (2)–(5), the track profile is obtained, shown as the solid line in Figure 1(a). The diameter, $D$, of the track on the surface is calculated using Equation (6). If the etchant has reached the particle range position R during the etching time, while the bulk etching depth is still less than the range (overetched case), as shown in Figure 1(b), the etchant first etches the trail to the range position at $v_{\text{T}}(z)$, and then etches around the wall at $v_{\text{B}}$ to form a spherical bottom. The solution of the track profile is similar to that in the non-overetched case. If the bulk etching depth exceeds the range of the particle (fully overetched case), the track profile is spherical, and the track image will have a low contrast and will usually be rejected.

(2)$$\begin{eqnarray}\displaystyle & z^{\prime }=\int _{0}^{t^{\prime }}v_{\text{T}}(z)dz, & \displaystyle\end{eqnarray}$$
(3)$$\begin{eqnarray}\displaystyle & \cos {\it\xi}=v_{\text{B}}/v_{\text{T}}, & \displaystyle\end{eqnarray}$$
(4)$$\begin{eqnarray}\displaystyle & z_{A}-z^{\prime }=v_{\text{B}}(t-t^{\prime })\cos {\it\xi}, & \displaystyle\end{eqnarray}$$
(5)$$\begin{eqnarray}\displaystyle & z_{A}=v_{\text{B}}t, & \displaystyle\end{eqnarray}$$
(6)$$\begin{eqnarray}\displaystyle & D=2y_{A}=2v_{\text{B}}(t-t^{\prime })\sin {\it\xi}. & \displaystyle\end{eqnarray}$$

Figure 1. The track profiles in (a) the non-overetched case and (b) the overetched case for perpendicular incident particles.

In the case of oblique incidence, the track profile is solved by Equations (2)–(4) in $y^{\prime }$$z^{\prime }$ coordinates and Equation (5) in $y$$z$ coordinates, as shown in Figure 2. The shape of the track on the surface of the etched CR39 is no longer circular. If the bulk etching depth is not very large, for example surface $a$ or $b$, the surface track shape will be an ellipse. If the etching depth is very large, for example surface $c$, the shape will be an egg shape. This case often occurs in the tracks induced by neutron recoil protons.

Figure 2. The track profiles for oblique incident particles.

Figure 3. The measured and calculated energy responses of CR39 to proton, deuteron, triton and alpha.

Figure 4. (a) The schematic setup of the D–D primary proton spectrometer at the SG-III prototype facility and (b) the measured proton spectra in the cool gas target (solid histogram) and the cryogenic $\text{D}_{2}$ target (dashed histogram) experiments with Gaussian fitting lines (blue lines).

2.3. Energy response calibration

The CR39 detector in this work was manufactured by Tastrak in UK. The energy responses of the track diameter to protons and alphas were measured at some energies, as the square and diamond points shown in Figure 3. For protons, the responses to 16 energies in the range of 0.7–3 MeV and three energies in the range of 4–10 MeV were measured at a $2\times 1.7~\text{MV}$ tandem accelerator and a $2\times 6~\text{MV}$ tandem accelerator, respectively, at Peking University. The uncertainties in the proton energy were less than 1%, and the incident proton beams were perpendicular to the CR39 surface. A timed aperture with a diameter of 1 cm was used to control the exposure time to avoid track saturation on the CR39. For alphas, the response to the 5.5 MeV alpha particles from a 241Am source was measured. To limit the incident angle and the number of the alpha particles, the timed aperture was placed in between, 20 cm away from the source and 10 cm away from the CR39 plastic. The area of the CR39 plastic used for track readout was approximately 1 cm2. The irradiated CR39 plastic was etched in $6~\text{mol}/\text{L}$ NaOH solution at a temperature of $80\,^{\circ }\text{C}$ for 6 h. By the mass loss of the unirradiated CR39 plastic before and after etching, the bulk etching velocity was measured as $\text{v}_{\text{B}}=1.793\pm 0.015~{\rm\mu}\text{m}/\text{h}$. By iterative fitting of the track diameters to the particle energies, the free parameters in Equation (1) were obtained as $k=3.3\times 10^{-4}$ and $n=2.4$. The calculated energy responses of the diameter to proton, deuteron, triton and alpha are also shown in Figure 3. With increase of the particle energy, the diameter decreases and the diameter decrease rate also decreases. Considering the track image spatial resolution and requirement of particle energy resolution, in routine proton measurements at the SG-III facility, readout proton tracks with diameters just in the range of 5–15 ${\rm\mu}\text{m}$ are used. Reference [Reference Sinenian, Rosenberg, Manuel, McDuffee, Casey, Zylstra, Rinderknecht, Gatu Johnson, Séguin, Frenje, Li and Petrasso21] reported that the piece-to-piece variation of the proton track diameter was approximately 19%. Considering the same etching conditions in this work as those in Ref. [Reference Sinenian, Rosenberg, Manuel, McDuffee, Casey, Zylstra, Rinderknecht, Gatu Johnson, Séguin, Frenje, Li and Petrasso21], in the following simulation work, a Gaussian distributed track diameter with the mean value of Figure 3 and a deviation of 19% was adopted. The detection efficiency was calculated as the probability that the track diameter would fall in the range of 5–$15~{\rm\mu}\text{m}$.

3. D–D fusion primary proton measurement by a filtered CR39 detector

In low yield $\text{D}_{2}$ fuel implosion experiments, the energy downshift of the primary proton spectra is a possible way to diagnose the areal density, which is less than $30~\text{mg}/\text{cm}^{2}$. The mean energy of the protons emitted in the hot-spot region is 3 MeV, with a Doppler broadening related to the ion temperature. When the areal density is less than 10 mg/cm2, the energy downshift is approximately linear with the areal density. In laser directly driven implosion, the disturbance comes from the ablator protons and the radial acceleration effect. Most ablator protons have low energy (${<}1~\text{MeV}$) and can be shielded with an appropriate filter. When the proton emission time falls in the laser pulse duration, the radial electric field in the corona region will accelerate the protons and contribute an energy upshift to the measured spectra, which will present a lower diagnosed ${\it\rho}$R. Distinct energy upshifts for different experimental conditions are necessary. However, the relatively different proton energy downshifts can be used to infer the variation of ${\it\rho}$R in similar shots. The typical energy upshift of D–D reaction protons at OMEGA is approximately 0.5 MeV[Reference Hicks, Li, Séguin, Ram, Frenje, Petrasso, Soures, Glebov, Meyerhofer, Roberts, Sorce, Stockl, Sangster and Phillips22].

In laser directly driven cool gas target and cryogenic $\text{D}_{2}$ target implosion experiments at the SG-III prototype facility, the D–D fusion primary proton spectra were also measured. The cool gas target has an inner radius of $381~{\rm\mu}\text{m}$ and a CH shell thickness of $9~{\rm\mu}\text{m}$, filled with $\text{D}_{2}$ gas (30 K, 200 kPa). The cryogenic target has an inner radius of $389~{\rm\mu}\text{m}$, a CH shell thickness of $9~{\rm\mu}\text{m}$ and an inner $\text{D}_{2}$ ice thickness of $40~{\rm\mu}\text{m}$. There were eight laser beams (1 ns pulse duration and 6.4 kJ total energy) driving the targets symmetrically. Due to the low fusion yield (${\sim}10^{7}$), an Al filtered CR39 detector was adopted to increase the detection efficiency. The schematic experimental setup is shown in Figure 4(a). The detector was set up below the cryogenic shielding, facing the laser entrance hole (LEH). The distance between the target and the CR39 was approximately 25 cm. Al filters with thicknesses of 20 and $25~{\rm\mu}\text{m}$ were used for the cool gas target and the cryogenic target, respectively. The CR39 recorded proton spectra were obtained based on the proton track histograms. With the TRIM calculated Al filter induced proton energy loss matrices in the energy range of 1–5 MeV, the target emission proton spectra were obtained as shown in Figure 4(b). It can be seen that the proton energy downshift of the cryogenic target was higher than that of the cool gas target, which represented a higher implosion areal density. Some further study of the energy upshift due to radial electric field acceleration will be carried out, to give the proton energy downshift in the fuel and shell and quantify the areal density.

Figure 5. (a) The schematic diagram of the SRF proton spectrometer and (b) the simulated proton energy response matrix.

4. D–3He reaction proton measurement by SRF spectrometer

4.1. Design and response simulation of the SRF spectrometer

In ICF implosions, D–3He reaction protons are emitted by the primary reaction in the D–3He fuel implosion and the secondary reaction in the D–D fuel implosion. The proton yield and the energy downshift can be used to deduce the areal density[Reference Séguin, Li, Frenje, Hicks, Green, Kurebayashi, Petrasso, Soures, Meyerhofer, Glebov, Radha, Stoeckl, Roberts, Sorce, Sangster, Cable, Fletcher and Padalino4, Reference Hicks19]. The high average proton energy (15 MeV) makes it possible to diagnose the areal density up to $100~\text{mg}/\text{cm}^{2}$. As a result, the energy dynamic range of the spectrometer is required to be wide, 8–18 MeV for example.

The SRF proton spectrometer uses a range of stepped filters to measure protons in different energy ranges[Reference Rosenberg, Zylstra, Frenje, Rinderknecht, Gatu Johnson, Waugh, Séguin, Sio, Sinenian, Li, Petrasso, Glebov, Hohenberger, Stoeckl, Sangster, Yeamans, LePape, Mackinnon, Bionta, Talison, Casey, Landen, Moran, Zacharias, Kilkenny and Nikroo12]. An SRF proton spectrometer has been built at the SG-III facility for D–3He reaction proton measurement. The schematic diagram is shown in Figure 5(a). A set of 12 Al filters, with thicknesses from 500 to $1800~{\rm\mu}\text{m}$ with an interval of 100 ${\rm\mu}\text{m}$, was adopted in the SRF spectrometer. The CR39 detector was divided into 12 parts to record the transmitted protons. Proton tracks with diameters in the range of 5–15 ${\rm\mu}\text{m}$ were used. The low energy protons were ranged out by the filter to each cell, while the high energy protons were not counted by the CR39. Every cell in the spectrometer looked like an energy band-pass filter to the protons. In the response function simulation, the track diameters of the transmitted protons were calculated according to a Gaussian distribution with the mean values of Figure 3 and a deviation of 19%. The response function of each cell to proton energy is shown in Figure 5(b). With the help of the energy response matrix, a spectrum unfolding code based on a weighted averaging algorithm was also developed.

Figure 6. The measured secondary proton spectra in the shot with a target diameter of $800~{\rm\mu}\text{m}$.

4.2. D–D implosion secondary proton measurement on the SG-III facility

The SRF proton spectrometer was first used in a $\text{D}_{2}$ exploding pusher target experiment on the SG-III facility in 2014. The spectrometer was installed on a general diagnostics instrument manipulator (DIM) and accessed the chamber in the middle plane. The distance between the target and the spectrometer was approximately 16 cm. The D–D secondary proton spectra in three shots were measured. The three targets had the same shell of $2~{\rm\mu}\text{m}$ thick $\text{SiO}_{2}$, while they had different diameters of 800, 620 and $450~{\rm\mu}\text{m}$. They were filled with 10 atm $\text{D}_{2}$ fuel. There were 32 laser beams directly driving the target, with a total energy of 30 kJ. In each shot, the hit positions on the target surface of the laser beams were shifted a little to decrease the root mean square (RMS) of the laser illumination. The measured proton spectrum for the $800~{\rm\mu}\text{m}$ diameter target is shown in Figure 6. The spectra in the other two shots were similar. The emitted protons were distributed mainly in the energy range of 12–18 MeV. There were no significant energy downshifts. This was consistent with expectations. Because the areal density of the exploding pusher target was very low, expected to be less than $1~\text{mg}/\text{cm}^{2}$, the energy downshift of secondary protons would be approximately 0.1 MeV, which was too small compared with the spectrometer energy resolution. However, the results show that the SRF spectrometer works well and is able to diagnose the areal density when it is higher than $10~\text{mg}/\text{cm}^{2}$. Moreover, from the spectra, the yields of the secondary protons were measured as $4.0\times 10^{6}$, $1.9\times 10^{6}$ and $9\times 10^{5}$, respectively. The neutron yields in the three shots were $3.1\times 10^{6}$, $1.9\times 10^{6}$ and $1.7\times 10^{6}$, respectively. A Monte Carlo code was developed to calculate the secondary proton emission, including the uniform hot-spot model and the 3He+ stopping power data calculated by the L–P model. According to the measured yield ratios of the secondary proton to primary neutron[Reference Cable and Hatchett23], the estimated areal densities of the three shots were 0.8, 0.7 and $0.4~\text{mg}/\text{cm}^{2}$, respectively.

5. A Si-WRF proton spectrometer for high energy proton diagnostics

5.1. Design and response simulation of the Si-WRF spectrometer

To overcome the disadvantage of low energy resolution in the SRF spectrometer, the WRF spectrometer was developed, which uses a wedged filter to replace the filter package. As illustrated in Figure 7(b), for incident mono-energetic protons, in the thicker end of the filter the protons will be ranged out, while in the thinner end the tracks of the transmitted protons will not be counted because of a track diameter smaller than $5~{\rm\mu}\text{m}$. The readout proton tracks are distributed in a narrow band on the surface of the CR39. By unfolding the proton track position spectra with the help of a well-defined energy response matrix, a high resolution incident proton spectrum can be obtained. There are already four types of WRF spectrometer, made from Al or $\text{ZrO}_{2}$ filters, developed at the OMEGA and NIF facilities[Reference Séguin, Sinenian, Rosenberg, Zylstra, Manuel, Sio, Waugh, Rinderknecht, Gatu Johnson, Frenje, Li, Petrasso, Sangster and Roberts13]. A WRF spectrometer, made from a Si filter, has been set up at the SG-III facility for 6–18 MeV proton measurement. The Si filter has similar characteristics to the Al filter in the particle range, but the deformation is much smaller during storage. The structure of a single spectrometer is illustrated in Figure 7(b). A WRF spectrometer array consisting of five spectrometers was developed, as shown in Figure 7(a). The thickness range of the Si filter is 120–$2000~{\rm\mu}\text{m}$. Both the width and length of the filter are 4 mm. The WRF array will be installed on a DIM and will access the chamber in the middle plane. There will be three view angles in the poloidal direction and three in the azimuthal direction, which may be useful in implosion asymmetry diagnostics. The distance between the target and the array can be adjusted in the range of 16–50 cm. The maximum poloidal and azimuthal view angles are $\pm 20^{\circ }$.

Figure 7. (a) The alignment of the Si-WRF spectrometer array at the SG-III facility (b) and a schematic illustration of a single WRF spectrometer[Reference Séguin, Frenje, Li, Hicks, Kurebayashi, Rygg, Schwartz, Petrasso, Roberts, Soures, Meyerhofer, Sangster, Knauer, Sorce, Glebov, Stoeckl, Phillips, Leeper, Fletcher and Padalino11].

A Monte Carlo simulation code based on the Geant4 package[Reference Agostinelli, Allisonas, Amako, Apostolakis, Araujo, Arce, Asai, Axen, Banerjee, Barrand, Behner, Bellagamba, Boudreau, Broglia, Brunengo, Burkhardt, Chauvie, Chuma, Chytracek, Cooperman, Cosmo, Degtyarenko, DellAcqua, Depaola, Dietrich, Enami, Feliciello, Ferguson, Fesefeldt, Folger, Foppiano, Forti, Garelli, Giani, Giannitrapani, Gibin, Cadenas, Gonzalez, Abril, Greeniaus, Greiner, Grichine, Grossheim, Guatelli, Gumplinger, Hamatsu, Hashimoto, Hasui, Heikkinen, Howard, Ivanchenko, Johnson, Jones, Kallenbach, Kanaya, Kawabata, Kawabata, Kawaguti, Kelner, Kent, Kimura, Kodama, Kokoulin, Kossov, Kurashige, Lamanna, Lampen, Lara, Lefebure, Lei, Liendl, Lockman, Longo, Magni, Maire, Medernach, Minamimoto, Freitas, Morita, Murakami, Nagamatu, Nartallo, Nieminen, Nishimura, Ohtsubo, Okamura, O’Neale, Oohata, Paech, Perl, Pfeiffer, Pia, Ranjard, Rybin, Sadilov, Di Salvo, Santin, Sasaki, Savvas, Sawada, Scherer, Sei, Sirotenko, Smith, Starkov, Stoecker, Sulkimo, Takahata, Tanaka, Tcherniaev, Tehrani, Tropeano, Truscott, Uno, Urban, Urban, Verderi, Walkden, Wander, Weber, Wellisch, Wenaus, Williams, Wright, Yamada, Yoshida and Zschiesche24] was developed to study the spectrometer response and evaluate its capability. The simulated proton position spectra for 9 and 15 MeV mono-energetic protons are shown in Figure 8(a). The 1–4 MeV transmitted proton tracks are distributed in a 3 mm wide band. This indicates that the bin width of the recorded proton position spectra should be no more than 1 mm. The proton energy responses in the range of 4–19 MeV were also simulated with an energy interval of 0.1 MeV, as shown in Figure 8(b). This shows that the designed WRF spectrometer is able to work well in the range of 6–18 MeV. It can also be used as a well-defined energy response matrix in the following spectra unfolding code.

Figure 8. (a) The energy response to 9 and 15 MeV mono-energetic protons and (b) the energy response matrix in the range of 4–19 MeV.

5.2. Capability evaluation

First, the capability of energy spectra reconstruction was evaluated. The test ‘experimental’ proton position spectra were simulated in the cases of different incident proton energy distributions. The spectra were unfolded by the MAXED code[Reference Reginatto, Goldhagen and Neumann25], which uses the maximum entropy algorithm. The energy response matrix was the simulated one in Figure 8(b). In the cases of mono-energetic proton incidence, the broadenings of the unfolded energy spectra were approximately 0.2 MeV, which represented the energy resolution of the WRF spectrometer. In the case of D–3He implosion primary proton incidence, two input Gaussian spectra were adopted with mean energies of 15 and 12 MeV and broadenings of 405 and 780 keV, respectively, as shown in Figure 9(a). The one with higher energy represented the proton emission spectrum at 5 keV ion temperature without energy downshift, while the other one represented the energy-downshifted proton spectrum in an areal density of 60–$70~\text{mg}/\text{cm}^{2}$. In the simulated ‘experiment’ with 106 proton incidence, the unfolded spectra are also shown in Figure 9(a). They are consistent with the input Gaussian spectra. In the case of D–D implosion secondary proton incidence, a rectangular spectrum of 12–17 MeV was adopted as the simulation input. The unfolded spectrum and the input spectrum are consistent, as shown in Figure 9(b).

Figure 9. (a) Comparisons of the unfolded spectra and the simulation input spectra in the cases of D–3He primary proton and (b) D–D secondary proton.

Figure 10. The unfolded energy spectra under different incident proton counts in the cases of (a) D–3He primary protons and (b) D–D secondary protons.

Second, the effect of the proton counts was investigated. In the case of D–3He implosion, with an incident proton count of 105–108, which represents a yield range of 108–1011, the unfolded spectra are shown in Figure 10(a). This shows that when the D–3He primary proton yield exceeds 108, the WRF spectrometer can work well. In the case of D–D implosion, with an incident secondary count of 104–107, which represents a primary yield range of 1010–1013, the unfolded spectra are shown in Figure 10(b). This shows that the necessary yield for WRF spectrometer measurement is approximately 1010.

6. Conclusions

In this work, three kinds of CR39-based proton spectrometers were applied at the SG-III facility for ICF implosion areal density diagnostics. The energy response of the CR39 detector to charged particles was studied by simulations and experiments. A Geant4 particle transport code including the track profile calculation code will be developed to study the background track profiles of the scattered particles and neutron recoil protons. A filtered CR39 detector has been applied in the measurement of D–D primary protons in low yield shots. A preliminary analysis of the areal density variation was implemented at the SG-III prototype facility. For the yields and spectra measurements of the D–3He reaction protons at the SG-III facility, a SRF spectrometer has already been applied in D–D fuel implosion areal density diagnostics. A Si-WRF spectrometer array has been designed and developed at the SG-III facility. Its capability was evaluated by Geant4 simulations in detail. In the following D–3He fuel implosion experiments, it will be applied in the diagnostics of areal density and asymmetry, mix effect studies and some plasma nuclear physics studies. By installing a CH foil in front, it can be upgraded to a compact nWRF spectrometer[Reference Zylstra, Gatu Johnson, Frenje, Séguin, Rinderknecht, Rosenberg, Sio, Li, Petrasso, McCluskey, Mastrosimone, Glebov, Forrest, Stoeckl and Sangster26] for D–T reaction neutron spectra measurement. Moreover, for high resolution measurement of the charged particle spectra, a magnet-based charged particle spectrometer is also under development at the SG-III facility.

Acknowledgements

This work was supported by the Open Project of the State Key Laboratory of Nuclear Physics and Technology at Peking University (No. 2014-02). The authors are very grateful to Professor Xiaotang Ren and Professor Hongji Ma in Peking University for the operation of the two accelerators. The authors are also very grateful to Professor Tieshuan Fan and his team in Peking University for their help with the spectra unfolding code.

References

Li, C. K. and Petrasso, R. D. Phys. Rev. Lett. 70, 3059 (1993).Google Scholar
Browna, L. S. Preston, D. L. and Singleton, R. L. Jr. Phys. Rep. 410, 237 (2005).Google Scholar
Zylstra, A. B. Frenje, J. A. Séguin, F. H. Rosenberg, M. J. Rinderknecht, H. G. Gatu Johnson, M. Casey, D. T. Sinenian, N. Manuel, M. J.-E. Waugh, C. J. Sio, H. W. Li, C. K. Petrasso, R. D. Friedrich, S. Knittel, K. Bionta, R. McKernan, M. Callahan, D. Collins, G. W. Dewald, E. Döppner, T. Edwards, M. J. Glenzer, S. Hicks, D. G. Landen, O. L. London, R. Mackinnon, A. Meezan, N. Prasad, R. R. Ralph, J. Richardson, M. Rygg, J. R. Sepke, S. Weber, S. Zacharias, R. Moses, E. Kilkenny, J. Nikroo, A. Sangster, T. C. Glebov, V. Stoeckl, C. Olson, R. Leeper, R. J. Kline, J. Kyrala, G. and Wilson, D. Rev. Sci. Instrum. 83, 10D901 (2012).Google Scholar
Séguin, F. H. Li, C. K. Frenje, J. A. Hicks, D. G. Green, K. M. Kurebayashi, S. Petrasso, R. D. Soures, J. M. Meyerhofer, D. D. Glebov, V. Yu. Radha, P. B. Stoeckl, C. Roberts, S. Sorce, C. Sangster, T. C. Cable, M. D. Fletcher, K. and Padalino, S. Phys. Plasmas 9, 2725 (2002).Google Scholar
Kitagawa, Y. Tanaka, K. A. Nakai, M. Yamanaka, T. Nishihara, K. Azechi, H. Miyanaga, N. Norimatsu, T. Kanabe, T. Chen, C. Richard, A. Sato, M. Furukawa, H. and Nakai, S. Phys. Rev. Lett. 75, 3130 (1995).Google Scholar
Frenje, J. A. Li, C. K. Séguin, F. H. Deciantis, J. Kurebayashi, S. Rygg, J. R. Petrasso, R. D. Delettrez, J. Glebov, V. Yu. Stoeckl, C. Marshall, F. J. Meyerhofer, D. D. Sangster, T. C. Smalyuk, V. A. and Soures, J. M. Phys. Plasmas 11, 2798 (2004).Google Scholar
Li, C. K. Séguin, F. H. Frenje, J. A. Rosenberg, M. Petrasso, R. D. Amendt, P. A. Koch, J. A. Landen, O. L. Park, H. S. Robey, H. F. Town, R. P. J. Casner, A. Philippe, F. Betti, R. Knauer, J. P. Meyerhofer, D. D. Back, C. A. Kilkenny, J. D. and Nikroo, A. Science 327, 1231 (2010).Google Scholar
Rygg, J. R. Frenje, J. A. Li, C. K. Séguin, F. H. Petrasso, R. D. Glebov, V. Y. Meyerhofer, D. D. Sangster, T. C. and Stoeckl, C. Phys. Rev. Lett. 98, 215002 (2007).Google Scholar
Rinderknecht, H. G. Sio, H. Li, C. K. Zylstra, A. B. Rosenberg, M. J. Amendt, P. Delettrez, J. Bellei, C. Frenje, J. A. Gatu Johnson, M. Séguin, F. H. Petrasso, R. D. Betti, R. Glebov, V. Yu. Meyerhofer, D. D. Sangster, T. C. Stoeckl, C. Landen, O. Smalyuk, V. A. Wilks, S. Greenwood, A. and Nikroo, A. Phys. Rev. Lett. 112, 135001 (2014).Google Scholar
Casey, D. T. Woods, D. T. Smalyuk, V. A. Hurricane, O. A. Glebov, V. Y. Stoeckl, C. Theobald, W. Wallace, R. Nikroo, A. Schoff, M. Wu, K. J. Frenje, J. A. Landen, O. L. Remington, B. A. and Glendinning, G. Phys. Rev. Lett. 114, 205002 (2015).Google Scholar
Séguin, F. H. Frenje, J. A. Li, C. K. Hicks, D. G. Kurebayashi, S. Rygg, J. R. Schwartz, B.-E. Petrasso, R. D. Roberts, S. Soures, J. M. Meyerhofer, D. D. Sangster, T. C. Knauer, J. P. Sorce, C. Glebov, V. Yu. Stoeckl, C. Phillips, T. W. Leeper, R. J. Fletcher, K. and Padalino, S. Rev. Sci. Instrum. 74, 975 (2003).Google Scholar
Rosenberg, M. J. Zylstra, A. B. Frenje, J. A. Rinderknecht, H. G. Gatu Johnson, M. Waugh, C. J. Séguin, F. H. Sio, H. Sinenian, N. Li, C. K. Petrasso, R. D. Glebov, V. Yu. Hohenberger, M. Stoeckl, C. Sangster, T. C. Yeamans, C. B. LePape, S. Mackinnon, A. J. Bionta, R. M. Talison, B. Casey, D. T. Landen, O. L. Moran, M. J. Zacharias, R. A. Kilkenny, J. D. and Nikroo, A. Rev. Sci. Instrum. 85, 103504 (2014).Google Scholar
Séguin, F. H. Sinenian, N. Rosenberg, M. Zylstra, A. Manuel, M. J.-E. Sio, H. Waugh, C. Rinderknecht, H. G. Gatu Johnson, M. Frenje, J. Li, C. K. Petrasso, R. Sangster, T. C. and Roberts, S. Rev. Sci. Instrum. 83, 10D908 (2012).Google Scholar
Fleischer, R. L. Price, P. B. and Walker, R. M. J. Appl. Phys. 36, 3645 (1965).Google Scholar
Caresana, M. Ferrarini, M. Fuerstner, M. and Mayer, S. Nucl. Instrum. Methods Phys. Res. A 683, 8 (2012).Google Scholar
Immè, G. Morelli, D. Aranzulla, M. Catalano, R. and Mangano, G. Radiat. Meas. 50, 253 (2013).Google Scholar
Zhou, D. Semones, E. Guetersloh, S. Zapp, N. Weyland, M. and Benton, E. R. Radiat. Meas. 45, 916 (2010).Google Scholar
Nikezic, D. and Yu, K. N. Comput. Phys. Comm. 174, 160 (2006).Google Scholar
Hicks, D. G. “Charged-particle spectroscopy: a new window on inertial confinement fusion”, PhD. Thesis (Massachusetts Institute of Technology, 1999).Google Scholar
Ziegler, J. F. Ziegler, M. D. and Biersack, J. P. Nucl. Instrum. Methods Phys. Res. B 268, 1818 (2010).Google Scholar
Sinenian, N. Rosenberg, M. J. Manuel, M. McDuffee, S. C. Casey, D. T. Zylstra, A. B. Rinderknecht, H. G. Gatu Johnson, M. Séguin, F. H. Frenje, J. A. Li, C. K. and Petrasso, R. D. Rev. Sci. Instrum. 82, 103303 (2011).Google Scholar
Hicks, D. G. Li, C. K. Séguin, F. H. Ram, A. K. Frenje, J. A. Petrasso, R. D. Soures, J. M. Glebov, V. Yu. Meyerhofer, D. D. Roberts, S. Sorce, C. Stockl, C. Sangster, T. C. and Phillips, T. W. Phys. Plasmas 7, 5106 (2000).Google Scholar
Cable, M. D. and Hatchett, S. P. J. Appl. Phys. 62, 2233 (1987).Google Scholar
Agostinelli, S. Allisonas, J. Amako, K. Apostolakis, J. Araujo, H. Arce, P. Asai, M. Axen, D. Banerjee, S. Barrand, G. Behner, F. Bellagamba, L. Boudreau, J. Broglia, L. Brunengo, A. Burkhardt, H. Chauvie, S. Chuma, J. Chytracek, R. Cooperman, G. Cosmo, G. Degtyarenko, P. DellAcqua, A. Depaola, G. Dietrich, D. Enami, R. Feliciello, A. Ferguson, C. Fesefeldt, H. Folger, G. Foppiano, F. Forti, A. Garelli, S. Giani, S. Giannitrapani, R. Gibin, D. Cadenas, J. J. G. Gonzalez, I. Abril, G. G. Greeniaus, G. Greiner, W. Grichine, V. Grossheim, A. Guatelli, S. Gumplinger, P. Hamatsu, R. Hashimoto, K. Hasui, H. Heikkinen, A. Howard, A. Ivanchenko, V. Johnson, A. Jones, F. W. Kallenbach, J. Kanaya, N. Kawabata, M. Kawabata, Y. Kawaguti, M. Kelner, S. Kent, P. Kimura, A. Kodama, T. Kokoulin, R. Kossov, M. Kurashige, H. Lamanna, E. Lampen, T. Lara, V. Lefebure, V. Lei, F. Liendl, M. Lockman, W. Longo, F. Magni, S. Maire, M. Medernach, E. Minamimoto, K. Freitas, P. M. Morita, Y. Murakami, K. Nagamatu, M. Nartallo, R. Nieminen, P. Nishimura, T. Ohtsubo, K. Okamura, M. O’Neale, S. Oohata, Y. Paech, K. Perl, J. Pfeiffer, A. Pia, M. G. Ranjard, F. Rybin, A. Sadilov, S. Di Salvo, E. Santin, G. Sasaki, T. Savvas, N. Sawada, Y. Scherer, S. Sei, S. Sirotenko, V. Smith, D. Starkov, N. Stoecker, H. Sulkimo, J. Takahata, M. Tanaka, S. Tcherniaev, E. Tehrani, E. S. Tropeano, M. Truscott, P. Uno, H. Urban, L. Urban, P. Verderi, M. Walkden, A. Wander, W. Weber, H. Wellisch, J. P. Wenaus, T. Williams, D. C. Wright, D. Yamada, T. Yoshida, H. and Zschiesche, D. Nucl. Instrum. Methods Phys. Res. A 506, 250 (2003).Google Scholar
Reginatto, M. Goldhagen, P. and Neumann, S. Nucl. Instrum. Methods Phys. Res. A 476, 242 (2002).Google Scholar
Zylstra, A. B. Gatu Johnson, M. Frenje, J. A. Séguin, F. H. Rinderknecht, H. G. Rosenberg, M. J. Sio, H. W. Li, C. K. Petrasso, R. D. McCluskey, M. Mastrosimone, D. Glebov, V. Y. Forrest, C. Stoeckl, C. and Sangster, T. C. Rev. Sci. Instrum. 85, 063502 (2014).Google Scholar
Figure 0

Figure 1. The track profiles in (a) the non-overetched case and (b) the overetched case for perpendicular incident particles.

Figure 1

Figure 2. The track profiles for oblique incident particles.

Figure 2

Figure 3. The measured and calculated energy responses of CR39 to proton, deuteron, triton and alpha.

Figure 3

Figure 4. (a) The schematic setup of the D–D primary proton spectrometer at the SG-III prototype facility and (b) the measured proton spectra in the cool gas target (solid histogram) and the cryogenic $\text{D}_{2}$ target (dashed histogram) experiments with Gaussian fitting lines (blue lines).

Figure 4

Figure 5. (a) The schematic diagram of the SRF proton spectrometer and (b) the simulated proton energy response matrix.

Figure 5

Figure 6. The measured secondary proton spectra in the shot with a target diameter of $800~{\rm\mu}\text{m}$.

Figure 6

Figure 7. (a) The alignment of the Si-WRF spectrometer array at the SG-III facility (b) and a schematic illustration of a single WRF spectrometer[11].

Figure 7

Figure 8. (a) The energy response to 9 and 15 MeV mono-energetic protons and (b) the energy response matrix in the range of 4–19 MeV.

Figure 8

Figure 9. (a) Comparisons of the unfolded spectra and the simulation input spectra in the cases of D–3He primary proton and (b) D–D secondary proton.

Figure 9

Figure 10. The unfolded energy spectra under different incident proton counts in the cases of (a) D–3He primary protons and (b) D–D secondary protons.