Instrumentation and application of the ion beam analysis line of the in situ ion beam system

Instrumentation and application of the ion beam analysis line of the in situ ion beam system 1 first-author 1 1 2 1 corresp Corresponding author, djfu@whu.edu.cn Corresponding author, djfu@whu.edu.cn djfu@whu.edu.cn Accelerator Laboratory, Key Laboratory of Artificial Micro- and Nano-Materials of Ministry of Education of China, School of Physics and Technology, Wuhan University, Wuhan 430072, China Accelerator Laboratory, Key Laboratory of Artificial Micro- and Nano-Materials of Ministry of Education of China, School of Physics and Technology, Wuhan University, Wuhan 430072, China State Key Laboratory of Advanced Welding and Joining, Harbin Institute of Technology, Harbin 150001, China State Key Laboratory of Advanced Welding and Joining, Harbin Institute of Technology, Harbin 150001, China An ion beam analysis system was established on a 1.7 MV tandem accelerator, enabling Rutherford backscattering (RBS), elastic recoil detection (ERD), nuclear reaction analysis (NRA) and channeling measurements. The system was tested by performing qualitative and quantitative analysis of Si, Ni/Si, BiFeO 3 :La/Si, MoC/Mo/Si and TiBN/Si samples. RBS of a BiFeO 3 :La film was used as system calibration. Tested by ion beam channeling, a Si(100) is of good crystallinity ( χ min = 3.01%). For thin film samples, the measured thickness agrees well with simulation results by SIMNRA. In particular, composition of a MoC/Mo/Si and TiBN film samples were analyzed by RBS and non-Rutherford elastic backscattering. Ion beam analysis Channeling Rutherford backscattering Non-Rutherford elastic backscattering I. INTRODUCTION 1 Based on interactions of energetic ions with target materials, ion beam analysis techniques are widely used in characterization of thin films and multilayer structures. They provide quantitative information of elemental components of the solid from several atomic layers to micrometer depth from the surface, or depth profiling of impurities. Therefore, it attracted much attention in the ion beam community and many ion beam systems have been developed and upgraded. These include conventional ion beam analysis, microbeams [ 1 1 - 4 4 1-4 ], and combination of ion beam systems, such as the 5 MV tandem accelerator interfaced with the 200 MeV proton synchrotron in the Wakasa Wan Energy Research Center [ 5 5 ] and the multi-ion irradiation platform JANNUS (Joint Accelerators for Nano Sciences and Nuclear Simulation) at Saclay where the first triple beam irradiation using Fe, He and H ion beams was carried out in 2010 [ 6 6 ]. Ion beam analysis methods are frequently used due to their efficiency in providing a broad spectrum of qualitative and quantitative information on elemental composition and crystal structure of solids. While Rutherford backscattering (RBS) provides composition and depth profiles of heavy elements in a light matrix, elastic recoil detection analysis (ERDA) is unique in profiling light elements (C, B, Li, H etc.) in a heavy matrix. Nuclear reaction analysis (NRA) provides specific elemental profiling which can resolve isotopes. Proton induced X-ray emission (PIXE) enables simultaneous detection of elements of intermediate atomic numbers. Ion beam channeling, compared with random ion beam scattering spectroscopy, gives high precision measurement of crystallinity of single-crystalline wafers or epitaxial layers [ 7 7 ]. Ion beam analysis methods have a common advantage of non-destructiveness in the sen se that the corresponding ion-solid interactions do not consume the target material, since the beam current is at analytical level, i.e., a few n A. Damage may result only, for sensitive samples, from ionization or heating by the ion beams, sputtering or nuclear stopping and in the case of channeling studies in single crystals, from defect production inducing dechanneling [ 8 8 ]. The accelerator system in Wuhan University is multifunction equipment, namely, the 1.7 MV tandem accelerator is interfaced with a 200 kV ion implanter and a transmission electron microscope (TEM). The ability of ion beam technique is extended to emerging fields of physics and materials science. In the present work, we have designed RBS/PIXE/C system aiming at characterization of solids with a broad range of elements in the same chamber and in situ measurement of the implanted sample in a dual beam chamber. The advantage of low-energy RBS technique lies in its capability of quantitative analysis of major and minor constituents lying in the first 0.5–2.0 μM of a material depending on the sample structure and composition. And the detection limits for some light elements may further be improved by using resonant scattering. In this paper we report the tests of basic functions of the system based on RBS and channeling, and demonstrate quantitative characterization of elemental profiling, and interface analysis, and crystal quality of wafers and thin film materials. II. SYSTEM DESIGN 1 Figure 1 Figure 1 shows a schematic diagram of the in situ ion beam system. Ion beams extracted from the negative ion source by cesium sputtering (SNICS) are transported to the low energy cluster ion beam chamber for low energy ion implantation or to the RBS/PIXE/C chamber after transportation lines of the tandem accelerator. The dual beam chamber is used for in situ implantation with ions from the tandetron or 200 kV implanter. It will also be possible to conduct in situ ion beam analysis of the implanted samples with the establishment of the RBS/PIXE/C system in the dual beam chamber. The ion beam analysis line thus consists of the SNICS, tandetron and analysis chamber. A schematic diagram of the measurement system is shown in Fig. 2 Fig. 2 . Conventional ion beam scattering, recoil detection and nuclear reaction are performed with Si:Au surface barrier detectors (SBD) connected through a charge sensitive preamplifier, an amplifier and a multichannel analyzer (MCA). The total number of ions is recorded by a charge integrator. For channeling measurement, a single channel analyzer (SCA) and a linear rate-meter are used. For PIXE measurement, we have an Ortec SLP-1080 detector with liquid nitrogen cooling and a miniature detector XR-100CR with thermal electric cooling. The SLP-1080 detector is put into the chamber through a cinlyndrical port. The XR-100CR detector is placed in the chamber, and the electronics for the RBS measurement can be used for signal amplification and recording. Ion beams extracted from the SNICS are focused and deflected before they were accelerated to 1–3 MeV, depending on energy requirement for different measurements. The accelerated ion beams are focused again by the electrostatic quadrupole lens and deflected into the analysis chamber by a magnetic switcher. The samples are placed on a ϕ 70 mm copper holder mounted on a four-axis goniometer, which is computer controlled through a stepping motor controller. During the measurement, the stepping motor controller controls movement of the goniometer, providing three orthogonal rotations (two tilts and one azimuth) and one translation ( y -axis). The z -axis is aligned with the incident ion beam direction. The y -axis is defined by orthogonality with both the x - and z -axis, which complete a right-hand triad. Specification of the sample orientation is described by three angles, γ, η and ζ, corresponding to rotation of the x -, y - and z -axes, respectively. The goniometer enables the samples to move along the y -axis or rotate around the x -, y - and z -axes. The spatial relationship between the sample and incident ion beam is described by the angle θ = ( η 2 + ζ 2 ) 1/2 and the azimuthal angle φ. In RBS analysis, the ion beam was collimated to 0.3 mm size by an aperture slot. The SBD, with an effective area of 50 mm 2 , is positioned at a distance of 10–15 cm to the sample and 150–170 to the beam incidence. The signals, after preamplifier and the spectroscopic amplifier, were sent to the MCA installed in the computer. The current integrator records the charge collected during the experiment. The RBS/C measurement has three more steps. First, signals from the amplifier are sent to a gate interface of the MCA through a single channel analyzer (SCA), simultaneously, an energy window is selected through coincidence by regulating the upper and lower level discriminators of the SCA. Next, the gate interface is disconnected and signals of the selected energy range from the SCA are sent to the linear rate-meters. Finally, tilt or rotation of the crystal sample relative to the beam incidence is controlled by the stepping motor controller, and the reduction of scattering yield is displayed on the linear rate-meters to decide channeling effect of the incident beams [ 9 9 ]. Detailed information of how to determine the channel of crystal sample by computerized control program was described in the literature [ 10 10 ]. III. ANALYSIS AND DISCUSSION 1 A. System calibration and determination of crystal quality 2 The RBS system was calibrated by measuring the Au peak position and Si edge of a Au-deposited Si sample. RBS measurement of a BiFeO 3 :La layer sample was performed. BiFeO 3 is a typical multiferroic material with ferroelectric and antiferromagnetic properties, has and is of broad application prospects in spintronics and multibit memory devices. Doping it with rare earth elements may significantly improve its ferroelectric and ferromagnetic properties. It is important to measure each element accurately. A thin film of La-doped BiFeO 3 was prepared by sol gel method on a Si(100) substrate. Fig. 3(a) Fig. 3(a) shows RBS spectrum of the sample using 2.9 MeV 7 Li 2+ beams with a backscattering angle of 160. The leading edge of Si is at channel 147 (at half-height) and the Bi, La, Fe and O peak positions are at channel numbers of 386, 362, 265, and 62, respectively. They correspond to the kinematic factors of 0.8778 for Bi, 0.8219 for La, 0.6126 for Fe, 0.3713 for Si, and 0.1607 for O. This indicates that Li ion beam has excellent mass resolution for many elements of film analysis. Assuming that the deposited La-BiFeO 3 thin layer is distributed uniformly on the substrate, the film thickness is estimated at 135 nm. Compared with other methods, the RBS analysis is very accurate, with an error of just 2%. Figure 3(b) Figure 3(b) shows RBS/C measurement results of a single crystal Si(100) sample using 2.75 MeV Li ion beams. Compare to the random spectrum, the backscattering yield greatly reduced in the aligned spectrum, the flat shape of which demonstrates that the incident ions go deeply into the crystal channels without dechanneling. At the high-energy edge, a small peak was observed, revealing an amorphous layer in the surface owning to the cutting during sample preparation. The minimum yield, χ min , defined as the backscattered yield ratio of the aligned spectrum to the random spectrum in the same channel region, i.e., the same depth, represents structural quality of the whole lattice of a monatomic crystal. For a perfect crystal, χ min could be 1–2% [ 9 9 ]. The measured χ min for the Si(100) sample is 3.01%, indicating a high quality crystal structure. B. Quantitative analysis of thin films 2 Thickness of a Ni film prepared by electron beam evaporation on SiO 2 :Si substrate was measured by RBS with 2.9 MeV 7 Li 2+ in normal incidence at a backscattering angle of 160 ( Fig. 4 Fig. 4 ). The yield of backscattered Li ions was detected by a Si:Au surface barrier detector with an effective detection area of 50 mm 2 , positioned at 112.6 mm from the sample, with a solid angle of Ω =3.94 msr. The measured spectrum overlaps with the theoretical spectrum simulated using the iterative analytical code of SIMNRA. The simulated areal density of the Ni film is 2.21×10 18 atoms/cm 2 , corresponding to a thickness of 242 nm. From the measured spectrum, the Ni peak area is A = 510000 cts, and the areal density Nt = 2.20×10 18 atoms/cm 2 , corresponding to a thickness of 240 nm, using the equation of Nt = A / σΩQ , where the scattering cross section of Ni being σ =1.25 bar from d σ /d Ω = [Z 1 Z 2 e 2 /(4 E )]2(4/ sin 2 θ ), and the number of incident ions being Q = 4.7×10 13 (7.5 μC). C. Composition analysis of multi-element and multilayer samples 2 A MoC thin film was prepared by multi-arc ion plating on Si substrate with Mo film as a transition layer. The Mo and C contents, and thickness of the MoC film, were determined by RBS analysis using a 1.75 MeV proton beam and the backscattering angle of θ = 170°. The RBS spectra are shown in Fig. 5(a) Fig. 5(a) together with a fit using the SIMNRA code. The Mo signal is approximately a rectangular peak with the front edge rising sharply. It indicates that molybdenum in the MoC layer is quite uniformly distributed and the hump at the back edge is the signal of Mo from transition layer. The carbon signal is a sharp peak rather than a hillock. Because non-Rutherford elastic backscattering occurred for carbon and the non-Rutherford backscattering cross section is 60 times larger than that of Rutherford backscattering at 1.74 MeV. [ 7 7 ] With 1.75 MeV protons, non-Rutherford elastic backscattering of carbon occurred at depth near the surface. A raised Si edge can also be observed due to the enhanced elastic scattering cross-section for Si atoms at 1.75 MeV of protons. According to the simulation result, the content of Mo is 0.71 and the content of C is 0.29 in the compound film layer with an area density of 1.52×10 19 atom/cm 2 . The Mo transition layer is 296.4 nm in thickness. Fig. 5(b) Fig. 5(b) shows RBS spectra of the sample measured under energies of 1.75–1.79 MeV. It can be seen that the spectrum broadened and the C peaks moved to the low energy region with increasing incident beam energy, due to that occurrence of non-Rutherford elastic backscattering of carbon deepened with the proton energy, and the energy resolution worsened with increasing proton energy. D. The detection and analysis of light elements in compound film 2 It is quite difficult to detect films containing light elements on a substrate of higher atomic mass by backscattering techniques, in which the yield of substrate is too high to distinguish signals of the light elements. In this case, one can use resonant scattering, which is of greater scattering cross section than Rutherford scattering. Appropriate higher energy should be applied so that resonance nuclear reaction occurs for certain elements. A sample of TiBN film deposited on Si substrate by multi-arc ion plating method was analyzed by 2 MeV proton with a backscattering angle of θ = 170°. Fig. 6 Fig. 6 is the energy spectra of the sample, with the fitting line and spectrum of the TiBN film after subtraction of Si substrate. The strong non-Rutherford resonances in 1 H- 14 N and 1 H- 11 B cross sections lead to sufficient enhancement of the N and B signals and areal density of the film can be calculated with the peak areas. The net integrated counts of the Ti, N and B peaks in the blue spectrum are A Ti = 186921 cts, A N = 118432 cts and A B = 47248 cts, respectively. A large uncertainty may be caused by subtraction of the substantial signal due to the Si substrate. Accuracy of the N and B yields is determined by the error of intercept and slope of the straight regression line y = y 0 + kx , where y is an approximation of the yield Y in channel x , y 0 is the intercept and k is the slope, best fitted into the experimental data. With 2 MeV protons, the cross section of Ti is Rutherford and the cross sections of N and B are non-Rutherford. The 1 H- 14 N cross-section enhancement factor is ( σ / σ R ) N = 5 and the 1 H- 11 B factor is ( σ / σ B ) N = 7 [ 7 7 ]. The cross section of Rutherford backscattering of Ti, B and N can be calculated as <math display="block" id="M1"><mrow><msub><mi>σ</mi><mrow><mtext>R,Ti</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>E</mi><mn>0</mn></msub><msup><mrow><mn>,170</mn><mo stretchy="false">)</mo><mo>=</mo><mo stretchy="false">(</mo><mn>0.6359</mn><mo>/</mo><mn>2</mn></mrow><mn>2</mn></msup><mo stretchy="false">)</mo><mo>×</mo><mi>d</mi><mo>−</mo><mn>24</mn><mo>=</mo><mo stretchy="false">[</mo><mi>p</mi><mi>e</mi><mi>r</mi><mo>−</mo><mi>m</mi><mi>o</mi><mi>d</mi><mi>e</mi><mo>=</mo><mi>s</mi><mi>y</mi><mi>m</mi><mi>b</mi><mi>o</mi><mi>l</mi><mo stretchy="false">]</mo><mn>1.590</mn><mi>d</mi><mo>−</mo><mn>25</mn><mi>s</mi><mi>r</mi><mn>,</mn></mrow></math> <math display="block" id="M2"><mrow><msub><mi>σ</mi><mrow><mtext>R,N</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>E</mi><mn>0</mn></msub><msup><mrow><mn>,170</mn><mo stretchy="false">)</mo><mo>=</mo><mo stretchy="false">(</mo><mn>0.0638</mn><mo>/</mo><mn>2</mn></mrow><mn>2</mn></msup><mo stretchy="false">)</mo><mo>×</mo><mi>d</mi><mo>−</mo><mn>24</mn><mo>=</mo><mo stretchy="false">[</mo><mi>p</mi><mi>e</mi><mi>r</mi><mo>−</mo><mi>m</mi><mi>o</mi><mi>d</mi><mi>e</mi><mo>=</mo><mi>s</mi><mi>y</mi><mi>m</mi><mi>b</mi><mi>o</mi><mi>l</mi><mo stretchy="false">]</mo><mn>1.595</mn><mi>d</mi><mo>−</mo><mn>26</mn><mi>s</mi><mi>r</mi><mn>,</mn></mrow></math> <math display="block" id="M3"><mrow><msub><mi>σ</mi><mrow><mtext>R,B</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>E</mi><mn>0</mn></msub><msup><mrow><mn>,170</mn><mo stretchy="false">)</mo><mo>=</mo><mo stretchy="false">(</mo><mn>0.0323</mn><mo>/</mo><mn>2</mn></mrow><mn>2</mn></msup><mo stretchy="false">)</mo><mo>×</mo><mi>d</mi><mo>−</mo><mn>24</mn><mo>=</mo><mo stretchy="false">[</mo><mi>p</mi><mi>e</mi><mi>r</mi><mo>−</mo><mi>m</mi><mi>o</mi><mi>d</mi><mi>e</mi><mo>=</mo><mi>s</mi><mi>y</mi><mi>m</mi><mi>b</mi><mi>o</mi><mi>l</mi><mo stretchy="false">]</mo><mn>8.078</mn><mi>d</mi><mo>−</mo><mn>27</mn><mi>s</mi><mi>r</mi><mn>.</mn></mrow></math> The number of backscattered particles detected using the surface energy approximation for normal incidence of beam can be expressed as A 0 = σ ( E 0 ) ΩQNt , where σ ( E 0 ) is scattering cross-section at incident proton energy E 0 , Ω is solid angle, Q is incident proton flux, N is number density of target atoms, and t is film thickness. Then, for Ti: N and Ti: B ratios from the backscattered spectrum, we have <math display="block" id="M4"><mrow><msub><mi>N</mi><mrow><mtext>Ti</mtext></mrow></msub><mo>/</mo><msub><mi>N</mi><mtext>N</mtext></msub><mo>=</mo><mo stretchy="false">(</mo><msub><mi>A</mi><mrow><mtext>Ti</mtext></mrow></msub><mo>/</mo><msub><mi>A</mi><mtext>N</mtext></msub><mo stretchy="false">)</mo><mo stretchy="false">[</mo><msub><mi>σ</mi><mrow><mtext>R,N</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>E</mi><mn>0</mn></msub><mn>,170</mn><mo stretchy="false">)</mo><mo>/</mo><msub><mi>σ</mi><mrow><mtext>R,Ti</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>E</mi><mn>0</mn></msub><mn>,170</mn><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo stretchy="false">(</mo><mi>σ</mi><mo>/</mo><msub><mi>σ</mi><mtext>R</mtext></msub><msub><mo stretchy="false">)</mo><mtext>N</mtext></msub><mn>,</mn></mrow></math> <math display="block" id="M5"><mrow><msub><mi>N</mi><mrow><mtext>Ti</mtext></mrow></msub><mo>/</mo><msub><mi>N</mi><mtext>B</mtext></msub><mo>=</mo><mo stretchy="false">(</mo><msub><mi>A</mi><mrow><mtext>Ti</mtext></mrow></msub><mo>/</mo><msub><mi>A</mi><mtext>B</mtext></msub><mo stretchy="false">)</mo><mo stretchy="false">[</mo><msub><mi>σ</mi><mrow><mtext>R,B</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>E</mi><mn>0</mn></msub><mn>,170</mn><mo stretchy="false">)</mo><mo>/</mo><msub><mi>σ</mi><mrow><mtext>R,Ti</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>E</mi><mn>0</mn></msub><mn>,170</mn><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo stretchy="false">(</mo><mi>σ</mi><mo>/</mo><msub><mi>σ</mi><mtext>R</mtext></msub><msub><mo stretchy="false">)</mo><mtext>B</mtext></msub><mn>.</mn></mrow></math> The average stoichiometric ratio for this film was calculated as Ti:N:B=1.4:1.75:1. IV. CONCLUSION 1 Measurements were carried out by ion beam analysis establishment on the 1.7 MV tandem accelerator at Wuhan University. The RBS spectrum calibration using a BiFeO 3 :La/Si sample indicates that the system works. Crystal quality of Si(100) was of χ min =3.01%, tested by RBS/C. A Ni film thickness was measured at 240 nm. A MoC film on Si substrate with a Mo transition layer was analyzed, the film content ratio being Mo:C= 0.71:0.29, with the Mo transition layer of 296.4 nm thick. Light element analysis of TiBN film was performed by non-Rutherford elastic backscattering, the stoichiometric ratio being Ti/N/B = 1.4:1.75:1. References C Jeynes , N P Barradas , M J Blewett , et al . Improved ion beam analysis facilities at the University of Surrey . Nucl Instrum Meth B , 1998 , 136-138 : 1229 - 1234 . DOI: 10.1016/S0168-583X(97)00818-5 http://doi.org/10.1016/S0168-583X(97)00818-5 Improved ion beam analysis facilities at the University of Surrey J G Lopez , F J Ager , M B Rank , et al . CNA: The first accelerator-based IBA facility in Spain . Nucl Instrum Meth B , 2000 , 161-163 : 1137 - 1142 . DOI: 10.1016/S0168-583X(99)00702-8 http://doi.org/10.1016/S0168-583X(99)00702-8 CNA: The first accelerator-based IBA facility in Spain M Uhrmacher , H Hofsass . Ion accelerator facilities at the University of Gottingen . Nucl Instrum Meth B , 2005 , 240 : 48 - 54 . DOI: 10.1016/j.nimb.2005.06.087 http://doi.org/10.1016/j.nimb.2005.06.087 Ion accelerator facilities at the University of Gottingen M Roumie , B Nsouli , K Zahraman , et al . First accelerator based ion beam analysis facility in Lebanon: development and applications . Nucl Instrum Meth B , 2004 , 219-220 : 389 - 393 . DOI: 10.1016/j.nimb.2004.01.088 http://doi.org/10.1016/j.nimb.2004.01.088 First accelerator based ion beam analysis facility in Lebanon: development and applications S Hatori , T Kurita , Y Hayashi , et al . Developments and applications of accelerator system at the Wakasa Wan Energy Research Center . Nucl Instrum Meth B , 2005 , 241 : 862 - 869 . DOI: 10.1016/j.nimb.2005.07.176 http://doi.org/10.1016/j.nimb.2005.07.176 Developments and applications of accelerator system at the Wakasa Wan Energy Research Center S Pellegrino , P Trocellier , S Miro , et al . The JANNUS Saclay facility: A new platform for materials irradiation,implantation and ion beam analysis . Nucl Instrum Meth B , 2012 , 273 : 213 - 217 . DOI: 10.1016/j.nimb.2011.07.078 http://doi.org/10.1016/j.nimb.2011.07.078 The JANNUS Saclay facility: A new platform for materials irradiation,implantation and ion beam analysis Y Q Wang and M A Nastasi . Handbook of modern ion beam materials analysis , Material Research Society , 2009 . Handbook of modern ion beam materials analysis G Amsel , G Battistig . The impact on materials science of ion beam analysis with electrostatic accelerators . Nucl Instrum Meth B , 2005 , 240 : 1 - 12 . DOI: 10.1016/j.nimb.2005.06.078 http://doi.org/10.1016/j.nimb.2005.06.078 The impact on materials science of ion beam analysis with electrostatic accelerators W K Chu , J W Mayer , M A Nicolet . Backscattering spectrometry , New York : Academic Press , 1978 . Backscattering spectrometry J He , J C Lee , M Li , et al . Computerized control and operation of rutherford backscattering/channeling for an in situ ion beam system and its application for measurement of Si(001) and ZnO(001) . Chin Phys Lett , 2011 , 28 : 012901 . DOI: 10.1088/0256-307X/28/1/012901 http://doi.org/10.1088/0256-307X/28/1/012901 Computerized control and operation of rutherford backscattering/channeling for an in situ ion beam system and its application for measurement of Si(001) and ZnO(001) 10.13538/j.1001-8042/nst.26.010202 Supported by the National Natural Science Foundation of China (No. 11405117) and the State Key Lab of Advanced Welding and Joining of Harbin Institute of Technology (No. AWJ-M13-03) 2014-03-12 2014-04-09 2015-02-02 2015-02-20 2021-04-11 2015 Nuclear Science and Techniques 1 26