Molecular dynamics simulations of Si etching in Cl- and Br-based plasmas: Cl+ and Br+ ion incidence in the presence of Cl and Br neutrals

Abstract

Classical molecular dynamics (MD) simulations have been performed for Cl+ and Br+ ions incident on Si(100) surfaces with Cl and Br neutrals, respectively, to gain a better understanding of the ion-enhanced surface reaction kinetics during Si etching in Cl- and Br-based plasmas. The ions were incident normally on surfaces with translational energies in the range Ei = 20–500 eV, and low-energy neutrals of En = 0.01 eV were also incident normally thereon with the neutral-to-ion flux ratio in the range Γn0/Γi0 = 0–100, where an improved Stillinger--Weber potential form was employed for the interatomic potential concerned. The etch yields and thresholds presently simulated were in agreement with the experimental results previously reported for Si etching in Cl2 and Br2 plasmas as well as in Cl+, Cl2+, and Br+ beams, and the product stoichiometry simulated was consistent with that observed during Ar+ beam incidence on Si in Cl2. Moreover, the surface coverage of halogen atoms, halogenated layer thickness, surface stoichiometry, and depth profile of surface products simulated for Γn0/Γi0 = 100 were in excellent agreement with the observations depending on Ei reported for Si etching in Cl2 plasmas. The MD also indicated that the yield, coverage, and surface layer thickness are smaller in Si/Br than in Si/Cl system, while the percentage of higher halogenated species in product and surface stoichiometries is larger in Si/Br. The MD further indicated that in both systems, the translational energy distributions of products and halogen adsorbates desorbed from surfaces are approximated by two Maxwellians of temperature T1 ≈ 2500 K and T2 ≈ 7000–40 000 K. These energy distributions are discussed in terms of the desorption or evaporation from hot spots formed through chemically enhanced physical sputtering and physically enhanced chemical sputtering, which have so far been speculated to both occur in the ion-enhanced surface reaction kinetics of plasma etching.

Molecular Dynamics Simulations for Plasma‐Surface Interactions

“…As a result, there is a rich literature on MD‐simulations of plasma etching, addressing topics such as the influence of the etchant gas, the energy of the plasma species, the role of concurrent ion bombardment, the surface temperature, etc. A typical example is given in ref . where MD simulations are applied to investigate the effect of Cl + and Br + ion impacts on the etching process of silicon in the presence of Cl and Br.…”

Different Etching Mechanisms of Diamond by Oxygen and Hydrogen Plasma: a Reactive Molecular Dynamics Study

“…Molecular dynamics (MD) simulation can be used to study the etching mechanisms at the atomic scale. Several researchers have already applied the MD simulations to the etching problems. , However, the traditional MD simulation cannot handle the chemical reactions, and the mechanisms of the etching process are still unclear. Additionally, researchers have used reactive MD methods such as ReaxFF − and tight-binding quantum chemical molecular dynamics (TB-QCMD) − to provide deeper insights into the detailed chemical reactions.…”

“…irradiation such as atomic O,20,21 atomic Cl and Br,16 SF 6 /O 2 , and CF 2 /CF 3 radicals 23,24. Thus, based on the above discussion, the etchant species used in our simulation are set as neutral atoms.…”

Molecular dynamic simulation study of plasma etching L1 FePt media in embedded mask patterning (EMP) process

“…Atoms in the next 6 layers upwards are "thermostated" to maintain the system at the set temperature 300K. 12 The 10 MLs above the thermostat layer are "mobile" layers that are only constrained through interatomic metallic bonds. The source of etch ion particles is confined in a 3D "slab" about 5A above the top surface.…”

“…Classical Molecular Dynamics (MD) simulation has been used to understand the ion sputtering and surface reaction kinetics during etching process. 9,[12][13][14] In MD model of etch process, substrate atoms, etchant particles (ions) and neutrals (molecules) are constructed in one ensemble with defined interatomic potentials. The time evolution of these particles are tracked and updated via solution of Newton's equation of motions.…”