Abstract: A hybrid model is developed for describing the effects of oxygen addition to argon glow discharges. The species taken into account in the model include Ar atoms in the ground state and the metastable level, O2 gas molecules in the ground state and two metastable levels, O atoms in the ground state and one metastable level, O3 molecules, Ar+, O+, O2+ and O− ions, as well as the electrons. The hybrid model consists of a Monte Carlo model for electrons and fluid models for the other plasma species. In total, 87 different reactions between the various plasma species are taken into account. Calculation results include the species densities and the importance of their production and loss processes, as well as the dissociation degree of oxygen. The effect of different O2 additions on these calculation results, as well as on the sputtering rates, is discussed.

Abstract: The purpose of this Brief Report is to point out the mistake in a formula for anisotropic electron scattering, previously published in Phys. Rev. A 41, 1112 (1990), which is widely used in Monte Carlo models of gas discharges. Anisotropic electron scattering is investigated based on the screened Coulomb potential between electrons and neutral atoms. The approach is also applied for electron scattering by nonpolar neutral molecules. Differential cross sections for electron scattering by Ar, N2, and CH4 are constructed on the basis of momentum and integrated cross sections. The formula derived in this paper is useful for Monte Carlo simulations of gas discharges.

Abstract: In this paper, we investigate the etch process of so-called low-k organic material by means of a N2/H2 capacitively coupled plasma, as applied in the micro-electronics industry for the manufacturing of computer chips. In recent years, such an organic material has emerged as a possible alternative for replacing bulk SiO2 as a dielectric material in the back-end-of-line, because of the smaller parasitic capacity between adjacent conducting lines, and thus a faster propagation of the electrical signals throughout the chip. Numerical simulations with a hybrid plasma model, using an extensive plasma and surface chemistry set, as well as experiments are performed, focusing on the plasma properties as well as the actual etching process, to obtain a better insight into the underlying mechanisms. Furthermore, the effects of gas pressure, applied power and gas composition are investigated to try to optimize the etch process. In general, the plasma density reaches a maximum near the wafer edge due to the so-called 'edge effect'. As a result, the etch rate is not uniform but will also reach its maximum near the wafer edge. The pressure seems not to have a big effect. A higher power increases the etch rate, but the uniformity becomes (slightly) worse. The gas mixing ratio has no significant effect on the etch process, except when a pure H2 or N2 plasma is used, illustrating the synergistic effects of a N2/H2 plasma. In fact, our calculations reveal that the N2/H2 plasma entails an ion-enhanced etch process. The simulation results are in reasonable agreement with the experimental values. The microscopic etch profile shows the desired anisotropic shape under all conditions under study.

Abstract: A one-dimensional fluid model for a dielectric barrier discharge in methane, used as a chemical reactor for gas conversion, is developed. The model describes the gas phase chemistry governing the conversion process of methane to higher hydrocarbons. The spatially averaged densities of the various plasma species as a function of time are discussed. Besides, the conversion of methane and the yields of the reaction products as a function of the residence time in the reactor are shown and compared with experimental data. Higher hydrocarbons (C2Hy and C3Hy) and hydrogen gas are typically found to be important reaction products. Furthermore, the main underlying reaction pathways are determined.

Abstract: Computer simulations are performed for an argon inductively coupled plasma (ICP) with a capacitive radio-frequency bias power, to investigate the bias effect on the discharge mode transition and on the plasma characteristics at various ICP currents, bias voltages, and bias frequencies. When the bias frequency is fixed at 13.56 MHz and the ICP current is low, e.g., 6A, the spatiotemporal averaged plasma density increases monotonically with bias voltage, and the bias effect is already prominent at a bias voltage of 90 V. The maximum of the ionization rate moves toward the bottom electrode, which indicates clearly the discharge mode transition in inductive/capacitive discharges. At higher ICP currents, i.e., 11 and 13 A, the plasma density decreases first and then increases with bias voltage, due to the competing mechanisms between the ion acceleration power dissipation and the capacitive power deposition. At 11 A, the bias effect is still important, but it is noticeable only at higher bias voltages. At 13 A, the ionization rate is characterized by a maximum at the reactor center near the dielectric window at all selected bias voltages, which indicates that the ICP power, instead of the bias power, plays a dominant role under this condition, and no mode transition is observed. Indeed, the ratio of the bias power to the total power is lower than 0.4 over a wide range of bias voltages, i.e., 0300V. Besides the effect of ICP current, also the effect of various bias frequencies is investigated. It is found that the modulation of the bias power to the spatiotemporal distributions of the ionization rate at 2MHz is strikingly different from the behavior observed at higher bias frequencies. Furthermore, the minimum of the plasma density appears at different bias voltages, i.e., 120V at 2MHz and 90V at 27.12 MHz.

Abstract: A two-dimensional fluid model, including the full set of Maxwell equations, has been developed and applied to investigate the effect of a phase shift between two power sources on the radial uniformity of several plasma characteristics in a hydrogen capacitively coupled plasma. This study was carried out at various frequencies in the range 13.56200 MHz. When the frequency is low, at 13.56 MHz, the plasma density is characterized by an off-axis peak when both power sources are in-phase (phgr = 0), and the best radial uniformity is obtained at phgr = π. This trend can be explained because the radial nonuniformity caused by the electrostatic edge effect can be effectively suppressed by the phase-shift effect at a phase difference equal to π. When the frequency rises to 60 MHz, the plasma density profiles shift smoothly from edge-peaked over uniform to centre-peaked as the phase difference increases, due to the pronounced standing-wave effect, and the best radial uniformity is reached at phgr = 0.3π. At a frequency of 100 MHz, a similar behaviour is observed, except that the maximum of the plasma density moves again towards the radial edge at the reverse-phase case (phgr = π), because of the dominant skin effect. When the frequency is 200 MHz, the bulk plasma density increases significantly with increasing phase-shift values, and a better uniformity is obtained at phgr = 0.4π. This is because the density in the centre increases faster than at the radial edge as the phase difference rises, due to the increasing power deposition Pz in the centre and the decreasing power density Pr at the radial edge. As the phase difference increases to π, the maximum near the radial edge becomes obvious again. This is because the skin effect has a predominant influence on the plasma density under this condition, resulting in a high density at the radial edge. Moreover, the axial ion flux increases monotonically with phase difference, and exhibits similar profiles to the plasma density. The calculation results illustrate that the radial uniformity of the various plasma characteristics is strongly dependent on the applied frequency and the phase shift between both power sources, which is important to realize, for controlling the uniformity of the plasma etch and deposition processes.

Abstract: A fluid model is self-consistently established to investigate the harmonic effects in an inductively coupled plasma, where the electromagnetic field is solved by the finite difference time domain technique. The spatiotemporal distribution of harmonic current density, harmonic potential, and other plasma quantities, such as radio frequency power deposition, plasma density, and electron temperature, have been investigated. Distinct differences in current density have been observed when calculated with and without Lorentz force, which indicates that the nonlinear Lorentz force plays an important role in the harmonic effects, especially at low frequencies. Moreover, the even harmonics are larger than the odd harmonics both in the current density and the potential. Finally, the dependence of various plasma quantities with and without the Lorentz force on various driving frequencies is also examined. It is shown that the deposited power density decreases and the depth of penetration increases slightly because of the Lorentz force. The electron density increases distinctly while the electron temperature remains almost the same when the Lorentz force is taken into account.

Abstract: PECVD of a nanoscale SiO2 capping layer using low pressure SiCl4/O-2/Ar plasmas is numerically investigated. The purpose of this capping layer is to restore photoresist profiles with improved line edge roughness. A 2D plasma and Monte Carlo feature profile model are applied for this purpose. The deposited films are calculated for various operating conditions to obtain a layer with desired shape. An increase in pressure results in more isotropic deposition with a higher deposition rate, while a higher power creates a more anisotropic process. Dilution of the gas mixture with Ar does not result in an identical capping layer shape with a thickness linearly correlated to the dilution. Finally, a substrate bias seems to allow proper control of the vertical deposition rate versus sidewall deposition as desired.

Abstract: We investigate the theoretical possibility of single layer graphene formation on a nickel surface at different substrate temperatures under far-from-equilibrium high precursor flux conditions, employing state-of-the-art hybrid reactive molecular dynamics/uniform acceptance force bias Monte Carlo simulations. It is predicted that under these conditions, the formation of a single layer graphene-like film may proceed through a combined depositionsegregation mechanism on a nickel substrate, rather than by pure surface segregation as is typically observed for metals with high carbon solubility. At 900 K and above, nearly continuous graphene layers are obtained. These simulations suggest that single layer graphene deposition is theoretically possible on Ni under high flux conditions.

Abstract: In this work, a hybrid model is used to investigate the effect of different gas ratios on the Si etching and polymer film deposition characteristics in an Ar/CF4 inductively coupled plasma. The influence of the surface processes on the bulk plasma properties is studied, and also the spatial characteristics of important gas phase and etched species. The densities of F and CF2 decrease when the surface module is included in the simulations, due to the species consumption caused by etching and polymer deposition. The influence of the surface processes on the bulk plasma depends on the Ar/CF4 gas ratio. The deposited polymer becomes thicker at high CF4 content because of more abundant CFx radicals. As a result of the competition between the polymer thickness and the F flux, the etch rate first increases and then decreases upon increasing the CF4 content. The electron properties, more specifically the electron density profile, affect the Si etch characteristics substantially by determining the radical density and flux profiles. In fact, the radial profile of the etch rate is more uniform at low CF4 content since the electron density has a smooth distribution. At high CF4 content, the etch rate is less uniform with a minimum halfway along the wafer radius, because the electron density distribution is more localized. Therefore, our calculations predict that it is better to work at relatively high Ar/CF4 gas ratios, in order to obtain high etch rate and good profile uniformity for etch applications. This, in fact, corresponds to the typical experimental etch conditions in Ar/CF4 gas mixtures as found in the literature, where Ar is typically present at a much higher concentration than CF4.