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: Carbon nanotubes (CNTs) are nowadays routinely grown in a thermal CVD setup. State-of-the-art plasma-enhanced CVD (PECVD) growth, however, offers advantages over thermal CVD. A lower growth temperature and the growth of aligned freestanding single-walled CNTs (SWNTs) makes the technique very attractive. The atomic scale growth mechanisms of PECVD CNT growth, however, remain currently entirely unexplored. In this contribution, we employed molecular dynamics simulations to focus on the effect of applying an electric field on the SWNT growth process, as one of the effects coming into play in PECVD. Using sufficiently strong fields results in (a) alignment of the growing SWNTs, (b) a better ordering of the carbon network, and (c) a higher growth rate relative to thermal growth rate. We suggest that these effects are due to the small charge transfer occurring in the Ni/C system. These simulations constitute the first study of PECVD growth of SWNTs on the atomic level.

Abstract: A hybrid model, called the hybrid plasma equipment model, was used to study inductively coupled SF6 plasmas used for Si etching applications. The plasma properties such as number densities of electrons, positive and negative ions, and neutrals are calculated under typical etching conditions. The electron kinetics is analysed by means of the electron energy probability function. The plasma chemistry taking place in pure SF6 and in an Ar/SF6 mixture is also discussed, and finally the effect of the argon fraction on the plasma properties is investigated.

Abstract: An inductively coupled plasma (ICP) is analyzed by means of experiments and numerical simulation. Important plasma properties are analyzed, namely, the effective temperature inside the central channel and the mean flow velocity inside the plasma. Furthermore, the effect of torches with different injector diameters is studied by the model. The temperature inside the central channel is determined from the end-on collected line-to-background ratio in dependence of the injector gas flow rates. Within the limits of 3% deviation, the results of the simulation and the experiments are in good agreement in the range of flow rates relevant for the analysis of relatively large droplets, i.e., 50 μm. The deviation increases for higher gas flow rates but stays below 6% for all flow rates studied. The velocity of the gas inside the coil region was determined by side-on analyte emission measurements with single monodisperse droplet introduction and by the analysis of the injector gas path lines in the simulation. In the downstream region significantly higher velocities were found than in the upstream region in both the simulation and the experiment. The quantitative values show good agreement in the downstream region. In the upstream region, deviations were found in the absolute values which can be attributed to the flow conditions in that region and because the methods used for velocity determination are not fully consistent. Eddy structures are found in the simulated flow lines. These affect strongly the way taken by the path lines of the injector gas and they can explain the very long analytical signals found in the experiments at low flow rates. Simulations were performed for different injector diameters in order to find conditions where good analyte transport and optimum signals can be expected. The results clearly show the existence of a transition flow rate which marks the lower limit for effective analyte transport conditions through the plasma. A rule-of-thumb equation was extracted from the results from which the transition flow rate can be estimated for different injector diameters and different injector gas compositions.

Abstract: In this paper, simulations are performed to gain a better insight into the properties of a Cl2/Ar plasma, with and without O2, during plasma etching of Si. Both plasma and surface properties are calculated in a self-consistent manner. Special attention is paid to the behavior of etch products coming from the wafer or the walls, and how the chamber walls can affect the plasma and the resulting etch process. Two modeling cases are considered. In the first case, the reactor walls are defined as clean (Al2O3), whereas in the second case a SiO2 coating is introduced on the reactor walls before the etching process, so that oxygen will be sputtered from the walls and introduced into the plasma. For this reason, a detailed reaction set is presented for a Cl2/O2/Ar plasma containing etched species, as well as an extensive reaction set for surface processes, including physical and chemical sputtering, chemical etching and deposition processes. Density and flux profiles of various species are presented for a better understanding of the bulk plasma during the etching process. Detailed information is also given on the composition of the surfaces at various locations of the reactor, on the etch products in the plasma and on the surface loss probabilities of the plasma species at the walls, with different compositions. It is found that in the clean chamber, walls are mostly chlorinated (Al2Cl3), with a thin layer of etch products residing on the wall. In the coated chamber, an oxy-chloride layer is grown on the walls for a few nanometers during the etching process. The Cl atom wall loss probability is found to decrease significantly in the coated chamber, hence increasing the etch rate. SiCl2, SiCl4 and SiCl3 are found to be the main etch products in the plasma, with the fraction of SiCl2 being always slightly higher. The simulation results compare well with experimental data available from the literature.

Abstract: A model for an atmospheric pressure inductively coupled plasma (ICP) is developed which allows rather easy extension to a variable number of species and ionisation degrees. This encompasses an easy calculation of transport parameters for mixtures, ionisation and heat capacity. The ICP is modeled in an axisymmetric geometry, taking into account the gas streaming into a flowing ambient gas. A mixture of argon and helium is applied in the injector gas stream as it is often done in laser ablation ICP spectrometry. The results show a strong influence of the added helium on the center of the ICP, which is important for chemical analysis. The length of the central channel is significantly increased and the temperature inside is significantly higher than in the case of pure argon. This means that higher gas volume flow rates can be applied by addition of helium compared to the use of pure argon. This has the advantage that the gas velocity in the transport system towards the ICP can be increased, which allows shorter washout-times. Consequently, shorter measurement times can be achieved, e.g. for spatial mapping analyses in laser ablation ICP spectrometry. Furthermore, the higher temperature and the longer effective plasma length will increase the maximum size of droplets or particles injected into the ICP that are completely evaporated at the detection site. Thus, we expect an increase of the analytical performance of the ICP by helium addition to the injector gas.

Abstract: In this review paper, an overview is given of different modelling efforts for plasmas used for the formation and growth of nanostructured materials. This includes both the plasma chemistry, providing information on the precursors for nanostructure formation, as well as the growth processes itself. We limit ourselves to carbon (and silicon) nanostructures. Examples of the plasma modelling comprise nanoparticle formation in silane and hydrocarbon plasmas, as well as the plasma chemistry giving rise to carbon nanostructure formation, such as (ultra)nanocrystalline diamond ((U)NCD) and carbon nanotubes (CNTs). The second part of the paper deals with the simulation of the (plasma-based) growth mechanisms of the same carbon nanostructures, i.e. (U)NCD and CNTs, both by mechanistic modelling and detailed atomistic simulations.

Abstract: The conversion of methane to value-added chemicals and fuels is considered to be one of the challenges of the 21st century. In this paper we study, by means of fluid modeling, the conversion of methane to higher hydrocarbons or oxygenates by partial oxidation with CO2 or O2 in a dielectric barrier discharge. Sixty-nine different plasma species (electrons, ions, molecules, radicals) are included in the model, as well as a comprehensive set of chemical reactions. The calculation results presented in this paper include the conversion of the reactants and the yields of the reaction products as a function of residence time in the reactor, for different gas mixing ratios. Syngas (i.e. H2 + CO) and higher hydrocarbons (C2Hx) are typically found to be important reaction products.

Abstract: In this article, surface processes occurring during the etching of Si with a Cl2/O2/Ar plasma are investigated by means of experiments and modeling. Cl2-based plasmas are commonly used to etch silicon, while a small fraction of O2 is added to protect the sidewalls from lateral etching during the shallow trench isolation process. When the oxygen fraction exceeds a critical value, the wafer surface process changes from an etching regime to a deposition regime, drastically reducing the etch rate. This effect is commonly referred to as the etch stop phenomenon. To gain better understanding of this mechanism, the oxygen fraction is varied in the gas mixture and special attention is paid to the effects of oxygen and of the redeposition of non-volatile etched species on the overall etch/deposition process. It is found that, when the O2 flow is increased, the etch process changes from successful etching to the formation of a rough surface, and eventually to the actual growth of an oxide layer which completely blocks the etching of the underlying Si. The size of this etch stop island was found to increase as a function of oxygen flow, while its thickness was dependent on the amount of Si etched. This suggests that the growth of the oxide layer mainly depends on the redeposition of non-volatile etch products. The abrupt change in the etch rate as a function of oxygen fraction was not found back in the oxygen content of the plasma, suggesting the competitive nature between oxidation and chlorination at the wafer. Finally, the wafer and reactor wall compositions were investigated by modeling and it was found that the surface rapidly consisted mainly of SiO2 when the O2 flow was increased above about 15 sccm.

Abstract: A Monte Carlo model has been developed for investigating the electron behavior in a dual-magnetron sputter deposition system. To describe the three-dimensional (3D) geometry, different reference frames, i.e. a local and a global coordinate system, were used. In this study, the influence of both closed and mirror magnetic field configurations on the plasma properties is investigated. In the case of a closed magnetic field configuration, the calculated electron trajectories show that if an electron is emitted in (or near) the center of the cathode, where the influence of the magnetic field is low, it is able to travel from one magnetron to the other. On the other hand, when an electron is created at the race track area, it is more or less trapped in the strong magnetic field and cannot easily escape to the second magnetron region. In the case of a mirror magnetic field configuration, irrespective of where the electron is emitted from the cathode, it cannot travel from one magnetron to the other because the magnetic field lines guide the electron to the substrate. Moreover, the electron density and electron impact ionization rate have been calculated and studied in detail for both configurations.

Abstract: In this paper, the effects of prolonged bias-enhanced nucleation (prolonged BEN) on the growth mechanisms of diamond are investigated by molecular dynamics (MD) and combined MD-Metropolis Monte Carlo (MD-MMC) simulations. First, cumulative impacts of CxHy+ and Hx+ on an a-C:H/nanodiamond composite were simulated; second, nonconsecutive impacts of the dominant ions were simulated in order to understand the observed phenomena in more detail. As stated in the existing literature, the growth of diamond structures during prolonged BEN is a process that takes place below the surface of the growing film. The investigation of the penetration behavior of CxHy+ and Hx+ species shows that the carbon-containing ions remain trapped within this amorphous phase where they dominate mechanisms like precipitation of sp3 carbon clusters. The H+ ions, however, penetrate into the crystalline phase at high bias voltages (>100 V), destroying the perfect diamond structure. The experimentally measured reduction of grain sizes at high bias voltage, reported in the literature, might thus be related to penetrating H+ ions. Furthermore, the CxHy+ ions are found to be the most efficient sputtering agents, preventing the build up of defective material.

Abstract: In this paper, an O2 inductively coupled plasma used for plasma enhanced atomic layer deposition of Al2O3 thin films is investigated by means of modeling. This work intends to provide more information about basic plasma properties such as species densities and species fluxes to the substrate as a function of power and pressure, which might be hard to measure experimentally. For this purpose, a hybrid model developed by Kushner et al is applied to calculate the plasma characteristics in the reactor volume for different chamber pressures ranging from 1 to 10 mTorr and different coil powers ranging from 50 to 500 W. Density profiles of the various oxygen containing plasma species are reported as well as fluxes to the substrate under various operating conditions. Furthermore, different orientations of the substrate, which can be placed vertically or horizontally in the reactor, are taken into account. In addition, special attention is paid to the recombination process of atomic oxygen on the different reactor walls under the stated operating conditions. From this work it can be concluded that the plasma properties change significantly in different locations of the reactor. The plasma density near the cylindrical coil is high, while it is almost negligible in the neighborhood of the substrate. Ion and excited species fluxes to the substrate are found to be very low and negligible. Finally, the orientation of the substrate has a minor effect on the flux of O2, while it has a significant effect on the flux of O. In the horizontal configuration, the flux of atomic oxygen can be up to one order of magnitude lower than in the vertical configuration.

Abstract: Using self-consistent KohnSham density-functional theory molecular dynamics simulations, we demonstrate the theoretical possibility to synthesize NiC60, the incarfullerene Ni@C60 and the heterofullerene C59Ni in an ion implantation setup. The corresponding formation mechanisms of all three complexes are elucidated as a function of the ion implantation energy and impact location, suggesting possible routes for selectively synthesizing these complexes.

Abstract: In this theoretical study, the electron emission from a copper surface under ultrashort pulsed laser irradiation is investigated using a one dimensional particle in cell model. Thermionic emission as well as multi-photon photoelectron emission were taken into account. The emitted electrons create a negative space charge above the target, consequently the generated electric field reduces the electron emission by several orders of magnitude. The simulations indicate that the space charge effect should be considered when investigating electron emission related phenomena in materials under ultrashort pulsed laser irradiation of metals.the word abstract, but do replace the rest of this text. ©2010 American Institute of Physics

Abstract: The behavior of hydrocarbon species at step edges of diamond terraces is investigated by means of combined molecular dynamics−Metropolis Monte Carlo simulations. The results show that the formation of ballas-like diamond films (like UNCD) and well-faceted diamond films (like NCD) can be related to the gas phase concentrations of CxHy in a new manner: Species that have high concentrations above the growing UNCD films suppress the extension of step edges through defect formation. The species that are present above the growing NCD film, however, enhance the extension of diamond terraces, which is believed to result in well-faceted diamond films. Furthermore, it is shown that, during UNCD growth, CxHy species with x ≥ 2 play an important role, in contrast to the currently adopted CVD diamond growth mechanism. Finally, the probabilities for the extension of the diamond (100) terrace are much higher than those for the diamond (111) terrace, which is in full agreement with the experimental observation that diamond (100) facets are more favored than diamond (111) facets during CVD diamond growth.