1. Introduction
Compact neutron generators based on 2H(d,n)3He (D-D) and 3H(d,n)4He (D-T) fusion reactions are among the most important neutron sources. With advantages over isotropic neutron sources in the radiation safety and the neutron output adjustability, they have important applications in scientific research and neutron application technologies, such as neutron imaging [1] and neutron activation analysis [2-4]. A compact D-D neutron generator (CDDNG) with a peak neutron output of up to 2.48 × 108 n/s is being developed at Lanzhou University for real-time neutron activation analysis [5].
The scheme of the CDDNG is shown in Fig. 1. The outline of the CDDNG is a cylinder with a length of 984 mm and a diameter of 234 mm. The CDDNG consists of an ion source, an extraction-accelerating electrode, a target assembly, a high voltage insulator assembly, a vacuum vessel, and a vacuum system. Deuterium ions are produced from a duoplasmatron ion source. To simplify the mode of power supply and ensure the safety and reliability during operation, a negative high voltage is applied to the target, and the ion source is kept at ground potential. The gap between the extraction-accelerating electrode and the plasma electrode of the ion source is 55 mm. Negative high voltage is applied to the extraction-accelerating electrode with the aid of the high voltage insulator assembly and a high voltage cable. The electric field is generated between the ion source and the extraction-accelerating electrode to extract and accelerate the deuterium ions to the target where neutrons are produced by D-D nuclear reactions. To achieve a longer target lifetime, a self-loading solid molybdenum target (50 mm in diameter and 1 mm in thickness) is used as the neutron production target. The surface of the target is 45° to the beam direction. A high vacuum degree (~10-4 Pa) in the vacuum vessel is provided by a set of vacuum pump systems, which consist of a mechanical pump and a turbo molecular pump. The target is cooled by fluorine containing coolant (3M Fluorinert FC-3283), which is pumped to the back of the target through the pipe. To achieve an excellent high-voltage insulation performance, the ceramic insulator tube is filled with transformer oil.
-201905/1001-8042-30-05-016/alternativeImage/1001-8042-30-05-016-F001.jpg)
Because the extraction-accelerating electrode and the target are at negative potential and the ion source is at ground potential, secondary electrons released from the target surface bombarded with deuterium ions[6] will be accelerated back to the ion source, producing a secondary electron current. This effect will increase the load of the high voltage power supply [7] and will produce frequent high voltage breakdown [8] and significant bremsstrahlung X-ray emission [9,10]. These effects affect both the accurate measurement of the deuterium ion current on the target [11,12] and the stable operation of the CDDNG [8], as well as posing a potential radiation safety hazard [9]. Therefore, it is imperative to find a suitable way to suppress the secondary electrons.
In this study, an electric field method and a magnetic field method [13-15] for suppressing the secondary electrons are designed and experimentally investigated. In addition, the secondary electron emission coefficient for deuterium ion impacting on molybdenum is estimated.
2. Experimental methods and results
2.1 Electric field method
The scheme of the designed electric field method for suppressing the secondary electrons is shown in Fig. 2. A ceramic ring is used to isolate the target from the extraction-accelerating electrode. A resistance (R) is used to connect the extraction-accelerating electrode to the target. When the deuterium current bombards the target and flows through the resistance to the extraction-accelerating electrode, a bias voltage is generated between the target and the extraction-accelerating electrode. The bias voltage (Vb) of the target relative to the extraction-accelerating electrode can be estimated with the following equation,
-201905/1001-8042-30-05-016/alternativeImage/1001-8042-30-05-016-F002.jpg)
where, Id is the deuterium ion current on the target. A current of secondary electrons, which is formed by the secondary electrons hitting the extraction-accelerating electrode, also affects the calculation of Vb through the resistance R to the target. The aim of this experiment is to obtain the threshold suppression voltage at which the secondary electrons are effectively suppressed. In this case, the secondary electron current can be negligible.
According to Eq.(1), when Id is constant, Vb increases linearly with increasing R. The secondary electrons should be completely suppressed when Vb is large enough. To study the performance of the secondary electron suppression of the electric field method shown in Fig. 2, experimental measurements of the high voltage load current (Iload) with the accelerating voltage (Vex) are taken under the same arc current (Iarc) of the ion source and various resistors (R = 0 kΩ, R = 10.2 kΩ, R = 51.1 kΩ, R = 102 kΩ, R = 133.2 kΩ, R = 200 kΩ, R = 400 kΩ, and R = 500 kΩ). The high voltage load current Iload is the read back of the high voltage power supply to the given high voltage. These resistance values are obtained by a combination of different metal film resistors. The rated power of a single metal film resistor is 2 W. The current Iload should be the sum of the extraction deuterium current (Id), the secondary electron current (Ie), and the leakage current of high voltage (Ileak):
The experimental results of Ileak, with Id = 0 mA, show that Ileak will not be greater than 0.1 mA. even at accelerating voltage of 150 kV, thus, Ileak can be ignored. Therefore, Iload can be approximately calculated as the sum of Id and Ie:
Fig. 3 shows the experimental results of Iload as a function of the accelerating voltage (Vex) for various resistances (R) under the same arc current (Iarc) of the ion source. As shown in Fig. 3, for R = 0 kΩ (Vb = 0 V), Iload increases rapidly with increasing Vex and, when Vex = 120 kV, Iload is up to 9.2 mA. As R increases, the increase in Iload becomes slower with increasing Vex and, at the same Vex, Iload decreases rapidly. This means that more and more secondary electrons are suppressed and, when R = 0 kΩ (Vb = 0 V), a large part of Ie is in part of the Iload. When R ≥102 kΩ and Vex ≥ 80 kV, Iload becomes a constant of 2 mA. This suggest that the secondary electrons have been completely suppressed, and Iload is equal to Id. In other words, the effective secondary electrons suppression can be achieved when the bias voltage between the target and the extraction-accelerating electrode is greater than or equal to 204 V (the product of R = 102 kΩ and Id = 2 mA from Eq. (1)).
-201905/1001-8042-30-05-016/alternativeImage/1001-8042-30-05-016-F003.jpg)
2.2 Magnetic field method
The designed magnetic field method for suppressing secondary electrons is shown in Fig. 4a. A pair of magnets are aligned with opposite poles facing each other to produce a magnetic field parallel to the target surface. The target is connected to the extraction-accelerating electrode through a wire. The Sm2Co17 permanent magnet is selected because of its high remanent flux density and thermostability. Under three different permanent magnet sizes (Φ 20 × 5 mm, Φ 20 × 10 mm, and Φ 20 × 12 mm), the magnetic field distributions along the center axis of the magnets (the red line in Fig. 4a) are measured with a Hall magnetometer. The results are shown in Fig. 4b. The magnetic field strength at the center of the target surface (X = 0 mm) is approximately 76 Gauss for Φ 20 × 5 mm, 131 Gauss for Φ 20 × 10 mm, and 150 Gauss for Φ 20 × 12 mm.
-201905/1001-8042-30-05-016/alternativeImage/1001-8042-30-05-016-F004.jpg)
To study the performance of secondary electron suppression of the magnetic field method shown in Fig. 4a, the load current (Iload) as a function of the accelerating voltage (Vex) were measured with and without a magnet (three different sizes) under the same arc current (Iarc) of the ion source. The experimental results are shown in Fig. 5. It can be seen that Iload is significantly lower with the Sm2Co17 magnets. The secondary electrons can be suppressed by the designed magnetic field method. When the magnet size is Φ 20 × 5 mm, the secondary electrons are not completely suppressed. The effects of secondary electron suppression with the magnets of size Φ 20 × 10 mm and Φ 20 × 12 mm are superior than those with the Φ 20 × 5 mm magnets.
-201905/1001-8042-30-05-016/alternativeImage/1001-8042-30-05-016-F005.jpg)
2.3 Comparison of the electrostatic and magnetic suppressions
Both the electric field and magnetic field methods have affected the suppression of secondary electrons. To compare their performance in secondary electrons suppression, the dependence of Iload on Vex for different arc currents (Iarc) of the ion source were measured with the electric field method with R = 500 kΩ and the magnetic field method with the magnet of size Φ 20 × 12 mm. During the experiment, the accelerating voltage (Vex) was increased over time from 15 kV to 120 kV, in steps of 5kV, and the arc current (Iarc) was increased over time from 0.4 A to 0.8 A, in steps of 0.2 A. The experimental results are shown in Fig. 6. As shown in the figure, for different Iarc, the Iload gradually reaches a stable value with increasing Vex in both methods. This confirms that both methods can effectively suppress secondary electrons. However, the Iload with the magnetic field method is significantly higher than that with the electric field method, while the Iarc of the ion source is higher in the former. A higher arc current (Iarc) means a bigger deuterium ion beam, that is, the suppression effect of the magnetic field will become worse as the deuterium ion beam increases. There are two possible reasons for this: one is the space charge effect, which can occur in the vicinity of the target. and the other is the demagnetization of the permanent magnet due to ion bombardment and heat conduction from the target.
-201905/1001-8042-30-05-016/alternativeImage/1001-8042-30-05-016-F006.jpg)
In summary, the designed electric field method is superior to the magnetic field method for secondary electrons suppression. Therefore, the electric field method will be employed in our CDDNG. For a better suppression effect in a lower deuterium ion current, the resistance should be 500 kΩ. The resistance is located in the transformer oil of the high voltage insulator assembly to achieve good cooling. After the implementation of the designed electric field method, the CDDNG was continuously operated for more than 3 hours without high voltage breakdown. When measures to suppress secondary electrons are not taken, high voltage breakdown occurs frequently, especially when the load current is close to 10 mA. Furthermore, the dose of X rays is reduced by about 10 times.
2.4 Estimation of the secondary electron emission coefficient
In this section, the secondary electron emission coefficient for deuterium ions incident on molybdenum is estimated based on the experimental data of the electrostatic suppression shown in Fig. 3. The secondary electron emission coefficient should be the ratio of the number of secondary electrons ejected from the target to the number of incident deuterium ions. It is assumed that the secondary electrons are not suppressed at all with R = 0 kΩ, and are completely suppressed with R = 500 kΩ and Vex> 80 kV. According to the experimental data in Fig. 3, the secondary electron emission coefficient λ of deuterium beam incident on molybdenum can be estimated by the following equation:
where, Ed, which is approximately equal to the product of Vex and the charge of deuterium ion, is the energy of the deuterium ion, Iload(0,
However, some secondary electrons can be suppressed even if Vb =0 V (i.e. R= 0 kΩ), owing to the target being installed inside the cavity of the cylindrical extraction-accelerating electrode (see Fig. 2). In other words, Iload(0,
It is quite difficult to measure the value of η. Therefore, a simulation investigation on the transport of the secondary electrons was carried out by a PIC code [16]. In the simulation, we skip the complex physical process of the secondary electrons produced by the deuterium ion interaction with the target [17], because electrons are produced as long as the deuterium ions are incident on the target surface. It is assumed that a deuterium ion produces three secondary electrons. The angular distribution of electrons ejected is set to be a cosine distribution. After contact with the target, the deuterium ions are 'killed', which means they disappear in the subsequent simulation, irrespective of the cascade effect of the secondary electrons, which means that the secondary electrons disappear directly after hitting any electrode and there is no secondary particle production. The current of deuterium ions is 2 mA.
A typical simulation image of the space distribution of the secondary electrons is shown in Fig. 7a, under Vex = 120 kV and R = 0 kΩ. Fig. 7b shows the potential distribution between the extraction-accelerating electrode and the target. The value of η can be calculated by counting the number of secondary electrons on each electrode. The calculated results of η are shown in Fig. 7c, with Vex increasing from 80 kV to 120 kV in steps of 10 kV. The value of η as a function of Vex shows a good linear relationship. The linear fitting values of η are used to calculate the coefficient λ. Fig. 8 shows the calculated results of the coefficient λ according to equation (5) with the correction of the ratio η. It can be seen that λ increases with increasing deuterium beam energy. In the energy range of 80 to 120 keV, the value of λ is between 5–5.5 with a deuterium ion glancing incidence of 45°.
-201905/1001-8042-30-05-016/alternativeImage/1001-8042-30-05-016-F007.jpg)
-201905/1001-8042-30-05-016/alternativeImage/1001-8042-30-05-016-F008.jpg)
The secondary electron emission coefficient on the metal target bombarded by deuterium ion has been widely investigated in the past decades [18-20]. However, most of these studies were carried out under the condition of normal incidence. According to references [21,22], the secondary electron emission coefficient with glancing incidence of 45° (λ45°) can be calculated by the following equation,
where λ⊥ is the secondary electron emission coefficient under normal incidence. For comparison with the previous data, we convert λ45° into λ⊥ according to equation (6). Fig. 8 shows the calculated results. The value of coefficient λ⊥ is between 3.5–3.9 at energies of 80–120 keV. The experimental data from reference [20] for the deuterium ions (D1+) normally incident on molybdenum are also shown in Fig. 8. It can be seen that λ45° is bigger than λ⊥. This is because the glancing incidence leads to a shallower penetration and it makes it easier for the generated carriers to exit the target surface. Furthermore, our estimated results (λ⊥) are significantly greater than the experimental data from reference [17]. The reason is that the mixed deuterium ion beam (D1+, D2+, and D3+) is used in our experiment. The previous experimental results show that the coefficient (λ⊥) for D2+ ions is about two times that of equal-energy D1+ ions [23]. According to our previous experimental results [24], the proportion of D1+ ions in the mixed deuterium ion beam (D1+, D2+, and D3+) is less than 10% under the condition of low arc current (such as Iarc< 1 A) for the duoplasmatron ion source.
3. Conclusion
An experimental study of two methods for secondary electrons suppression in CDDNG was carried out. The methods are the electrostatic suppression and the magnetic suppression. Both methods can effectively suppress secondary electrons. By comparison, the designed electric field method is superior to the designed magnetic field method. Therefore, the electric field method will be applied in our CDDNG. The resistance will be taken as 500 kΩ for a better suppression effect in a lower deuterium ion current.
The secondary electron emission coefficient for the mixed deuterium ions incident on molybdenum is estimated based on the experimental data of the electrostatic suppression. In the deuterium energy range of 80 to 120 keV, the estimated secondary electron emission coefficients are approximately 5 to 5.5 for mixed deuterium ion with glancing incidence of 45° and 3.5 to 3.9 for mixed deuterium ion with normal incidence. The above results may help to analyze the secondary electron problem in neutron tubes.
Neutron Generators for Analytical Purposes, IAEA Radiation Technology Reports No. 1
,An neutron generator-based NIPGA system for on-site analysis
. Nucl. Sci. Tech., 21:Instrumental neutron activation analysis of beachrock samples of South East Coast of Tamilnadu, India
. Nucl. Sci. Tech., 18:204-211(2007). doi: 10.1016/S1001-8042(07)60047-5Development of a compact DD neutron generator
, J. Instrum. 13(01), P01013 (2018). doi: 10.1088/1748-0221/13/01/P01013The emission of secondary electrons under high energy positive ion bombardment
, Phys. Rev. 55(5), 463 (1939). doi: 10.1103/PhysRev.55.463Magnetic suppression of secondary electrons in plasma immersion ion implantation
, Appl. Phys. Lett. 86(2), 023509 (2005). doi: 10.1063/1.1852704Role of Positive Ions in High-Voltage Breakdown in Vacuum
, J. Appl. Phys. 26(5), 596-599 (1955). doi: 10.1063/1.1722047Development and characterization of a D–D fast neutron generator for imaging applications
, Appl. Radiat. Isotopes. 96, 114-121 (2015). doi: 10.1016/j.apradiso.2014.11.017An algorithm for Monte Carlo simulation of bremsstrahlung emission by electrons
. Nucl. Sci. Tech., 28: 71(2017). doi: 10.1007/s41365-017-0218-7A neutron tube with constant output (1010 n/sec) for activation analysis and reactor applications
, Nucl. Instrum. Methods. 33(2), 283-288 (1965). doi: 10.1016/0029-554X(65)90055-8Low Voltage 14 MeV Neutron Source
, Rev. Sci. Instrum. 24(6), 424-427 (1953). doi: 10.1063/1.1770737Neutron generation in small sealed accelerating tubes, Russ
. Phys. J 28(8), 609-632 (1985). doi: 10.1007/BF00895162Development of D–D neutron generator
, Nucl. Instrum. Meth. B 241(1-4), 917-920 (2005). doi: 10.1016/j.nimb.2005.07.170Beam-induced back-streaming electron suppression analysis for an accelerator type neutron generator designed for 40Ar/39Ar geochronology
, Appl. Radiat. Isotopes. 125, 124-128 (2017). doi: 10.1016/j.apradiso.2017.04.017User-configurable MAGIC for electromagnetic PIC calculations, Comput
. Phys. Commun. 87(1-2), 54-86(1995).Calculations of Secondary Electron Yield of Graphene Coated Copper for Vacuum Electronic Applications
, AIP Advances 8(1),015325 (2018). doi: 10.1063/1.5019360Secondary electron emission from scandium, erbium, scandium deuteride, and erbium deuteride under deuteron bombardment
, J. Appl. Phys. 43(8),3318-3324 (1972). doi: 10.1063/1.1661714Deuteron induced secondary electron emission from titanium deuteride surface
, Nucl. Instrum. Meth. B 280,1-4 (2012). doi: 10.1016/j.nimb.2012.02.033Data compendium for plasma-surface interactions
, Nucl. Fusion. 24(S1), S9(1984). doi: 10.1088/0029-5515/24/S1/001Angular dependence of the ion-induced secondary-electron yield from solids
, Phys. Rev. B 24(7), 3749, (1981). doi: 10.1103/PhysRevB.24.3749High-voltage glow discharges in D2 gas. I. Diagnostic measurements
, Phys. Rev. 124(4), 969 (1961). doi: 10.1103/PhysRev.124.969