• Journal of the Korean Vacuum Society
• /
• v.12 no.3
• /
• pp.162-167
• /
• 2003

To find the concentration according to the depth-direction of ions implanted in the sample, with sputtering of the sample surface, one needs the depth profiling of ion implanted in the sample. On measuring of depth profiling, the sputtering rate to affect depth direction, is calculated by SRIM simulation. When ion is implanted in the sample, the atomic density of the sample rises up a little, and it alters sputtering yield. This alteration then causes differences of sputtering rate to affect depth-direction, on measuring of depth profiling. With the usage of SRIM Monte Carlo simulation code, one calculates sputtering rate, with sputtering yield by the alteration of atomic density of the sample through ion implantation. As a result, it goes to prove that its difference affects depth distribution, on measuring of depth profiling.

• Nuclear Engineering and Technology
• /
• v.53 no.10
• /
• pp.3158-3163
• /
• 2021

Chromium material is commonly used for fusion plasma facing applications because of the low neutron activation property. The Monte Carlo method is one of the useful ways to investigate the ion-target interactions. In this study, Chromium target irradiated by protons was investigated using Monte Carlo based simulation tools. In this context, the calculations of radiation damage on Chromium material irradiated with protons at 17.9 and 22.3 MeV energies were carried out using GEANT4 and SRIM codes. Besides, the cross sections for proton interaction with Chromium target were calculated by the TALYS 1.9 code using CTM + FGM, BSFGM, and GSFM level densities. As a result, GEANT4, SRIM and TALYS 1.9 codes provide a suitable tool for the predictions of radiation damage and cross cross section with proton irradiation.

• Nuclear Engineering and Technology
• /
• v.52 no.4
• /
• pp.827-834
• /
• 2020

Damage of tungsten due to helium ions of a PF device was studied. The tungsten was analyzed by SEM and AFM after irradiation. SEM revealed fine bubbles of helium atoms with diameters of a few nanometers, which join and form larger bubbles and blisters on the surface of tungsten. This observation confirmed the results of molecular dynamics simulation. SEM analysis after etching of the irradiated surface indicated cavities with depth range of 35-85 nm. The average fluence of helium ion of the PF device was calculated about 5.2 × 10¹⁵ cm^-2 per shot, using Lee code. Energy spectrum of helium ions was estimated using a Thomson parabola spectrometer as a function of dN/dE ∝ E^-2.8 in the energy range of 10-200 keV. The characteristics of helium ion beam was imported to SRIM code. SRIM revealed that the maximum DPA and maximum helium concentration occur in the depth range of 20-50 nm. SRIM also showed that at depth of 30 nm, all of the tungsten atoms are displaced after 20 shots, while at depth of higher than 85 nm the destruction is insignificant. There is a close match between SRIM results and the measured depths of cavities in SEM images of tungsten after etching.

• Journal of the Korean Institute of Electrical and Electronic Material Engineers
• /
• v.17 no.11
• /
• pp.1156-1166
• /
• 2004

Antimony profiles by MeV implantation are measured by secondary ion mass spectrometry (SIMS) and spreading resistance (SR). The moments of SIMS and simulated profiles are calculated and compared for the exact range in MeV energy. SRIM, DUPEX, ICECREM, and TSUPREM4 simulation programs are used for the calculation of range 1D, 2D. SRIM is a Monte Carlo simulation program and different inter-atomic potentials can be used for the calculation of nuclear stopping power cross-section (Sn) and range moments. Nevertheless, the range parameters were not influenced from nuclear stopping power in MeV. Through the modification of electronic stopping power cross-section (Se), the results of simulation are remarkably improved and matched very well with SIMS data. The values of electronic stopping power are optimized for Sb high energy implantation. For the electrical activation, Sb implanted samples are annealed under $N_2$ and $O_2$ ambient. Finally, Oxidation retard diffusion(ORD) effect of Sb implanted sample are demonstrated by SR measurements and ICECREM simulation.

• Journal of the Korean Institute of Electrical and Electronic Material Engineers
• /
• v.29 no.12
• /
• pp.750-758
• /
• 2016

For the channel doping of shallow junction and retrograde well formation in CMOS, indium can be implanted in silicon. The retrograde doping profiles can serve the needs of channel engineering in deep MOS devices for punch-through suppression and threshold voltage control. Indium is heavier element than B, $BF_2$ and Ga ions. It also has low coefficient of diffusion at high temperatures. Indium ions can be cause the erode of wafer surface during the implantation process due to sputtering. For the ultra shallow junction, indium ions can be implanted for p-doping in silicon. UT-MARLOWE and SRIM as Monte carlo ion-implant models have been developed for indium implantation into single crystal and amorphous silicon, respectively. An analytical tool was used to carry out for the annealing process from the extracted simulation data. For the 1D (one-dimensional) and 2D (two-dimensional) diffused profiles, the analytical model is also developed a simulation program with $C^{{+}{+}}$ code. It is very useful to simulate the indium profiles in implanted and annealed silicon autonomously. The fundamental ion-solid interactions and sputtering effects of ion implantation are discussed and explained using SRIM and T-dyn programs. The exact control of indium doping profiles can be suggested as a future technology for the extreme shallow junction in the fabrication process of integrated circuits.

• Nuclear Engineering and Technology
• /
• v.55 no.4
• /
• pp.1287-1299
• /
• 2023

A new Monte-Carlo-based computer program (RDS-BASIC) is developed to simulate the transport of energetic ions in pure matter. This computer program is utilizing an algorithm that uses detailed numerical solutions for the classical scattering integral for evaluating the outcomes of the binary collision processes. This approach is adopted by several prominent similar simulation programs and is known to provide results with higher accuracy compared to other approaches that use approximations to shorten the simulation time. Furthermore, RDS-BASIC simulation program contains special methods to reduce the displacement energy threshold of surface atoms. This implementation is found essential for accurate simulation results for sputtering yield in the case of very low energy ions irradiation (near sputtering energy threshold) and also successfully solve the problem of simultaneously obtaining an acceptable number of atomic displacements per incident ions. Results of our simulation for several irradiation systems are presented and compared with their respective TRIM (SRIM-2013) and the state-of-the-art SDTrimSP simulation results. Our sputtering simulation results were also compared with available experimental data. The simulation execution time for these different simulation programs has also been compared.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2012.02a
• /
• pp.526-526
• /
• 2012

DC magnetron discharges were studied using three-dimensional self-consistent particle-in-cell and Monte Carlo collisional (PIC-MCC) simulation codes. Two rectangular sputter sources (120 mm * 250 mm and 380 mm * 200 mm target sizes) were used in the simulation modeling. The number of incident ions to the Cu target as a function of position and simulation time was obtained. The target erosion profile was calculated by using the incident ions and the sputtering yields of the Cu target calculated with SRIM codes. The maximum ion density of the ion density distribution in the discharge was about $10^{10}cm^{-3}$ due to the calculation speed limit. The result may be less than one or two order of magnitude smaller than the real maximum ion density. However, the target erosion profiles of the two sputter sources were in good agreement with the measured target erosion profiles except for the erosion profile near the target surface, in which which the measured erosion width was broader than the simulation erosion width.

• Journal of the Korean Institute of Electrical and Electronic Material Engineers
• /
• v.27 no.11
• /
• pp.679-685
• /
• 2014

For the measurements of surface shape milled using FIB (focused ion beam), the silicon bulk, $Si_3N_4/Si$, and Al/Si samples are used and observed the shapes milled from different sputtering rates, incident angles of $Ga^+$ ions bombardment, beam current, and target material. These conditions also can be influenced the sputtering rate, raster image, and milled shape. The fundamental ion-solid interactions of FIB milling are discussed and explained using TRIM programs (SRIM, TC, and T-dyn). The damaged layers caused by bombarding of $Ga^+$ ions were observed on the surface of target materials. The simulated results were shown a little bit deviation with the experimental data due to relatively small sputtering rate on the sample surface. The simulation results showed about 10.6% tolerance from the measured data at 200 pA. On the other hand, the improved analytical model of damaged layer was matched well with experimental XTEM (cross-sectional transmission electron microscopy) data.

• Journal of the Korean Institute of Electrical and Electronic Material Engineers
• /
• v.31 no.6
• /
• pp.377-385
• /
• 2018

For the investigation of dopant profiles in implanted $Si_{1-x}Ge_x$, the implanted B and As profiles are measured using SIMS (secondary ion mass spectrometry). The fundamental ion-solid interactions of implantation in $Si_{1-x}Ge_x$ are discussed and explained using SRIM, UT-marlowe, and T-dyn programs. The annealed simulation profiles are also analyzed and compared with experimental data. In comparison with the SIMS data, the boron simulation results show 8% deviations of $R_p$ and 1.8% deviations of ${\Delta}R_p$ owing to relatively small lattice strain and relaxation on the sample surface. In comparison with the SIMS data, the simulation results show 4.7% deviations of $R_p$ and 8.1% deviations of ${\Delta}R_p$ in the arsenic implanted $Si_{0.2}Ge_{0.8}$ layer and 8.5% deviations of $R_p$ and 38% deviations of ${\Delta}R_p$ in the $Si_{0.5}Ge_{0.5}$ layer. An analytical method for obtaining the dopant profile is proposed and also compared with experimental and simulation data herein. For the high-speed CMOSFET (complementary metal oxide semiconductor field effect transistor) and HBT (heterojunction bipolar transistor), the study of dopant profiles in the $Si_{1-x}Ge_x$ layer becomes more important for accurate device scaling and fabrication technologies.

• Journal of the Korean Institute of Electrical and Electronic Material Engineers
• /
• v.24 no.8
• /
• pp.609-615
• /
• 2011

This paper is studied for the improvement of the characteristics of gate oxide with 3-nm-thick gate oxide by deuterium ion implantation methode. Deuterium ions were implanted to account for the topography of the overlaying layers and placing the D peak at the top of gate oxide. A short anneal at forming gas to nitrogen was performed to remove the damage of D-implantation. We simulated the deuterium ion implantation to find the optimum condition by SRIM (stopping and range of ions in matter) tool. We got the optimum condition by the results of simulation. We compare the electrical characteristics of the optimum condition with others terms. We also analyzed the electrical characteristics to change the annealing conditions after deuterium ion implantation. The results of the analysis, the breakdown time of the gate oxide was prolonged in the optimum condition. And a variety of annealing, we realized the dielectric property that annealing is good at longer time. However, the high temperature is bad because of thermal stress.