Abstract: Modeling results are presented to gain a better insight in the properties of a SiH4/O2/Ar inductively coupled plasma (ICP) and how it interacts with a silicon substrate (wafer), as applied in the microelectronics industry for the fabrication of electronic devices. The SiH4/O2/Ar ICP is used for the filling of microtrenches with isolating material (SiO2), as applied in shallow trench isolation (STI). In this article, a detailed reaction set that describes the plasma chemistry of SiH4/O2/Ar discharges as well as surface processes, such as sputtering, oxidation, and deposition, is presented. Results are presented on the plasma properties during the plasma enhanced chemical vapor deposition process (PECVD) for different gas ratios, as well as on the shape of the filled trenches and the surface compositions of the deposited layers. For the operating conditions under study it is found that the most important species accounting for deposition are SiH2, SiH3O, SiH3 and SiH2O, while SiH+2, SiH+3, O+2 and Ar+ are the dominant species for sputtering of the surface. By diluting the precursor gas (SiH4) in the mixture, the deposition rate versus sputtering rate can be controlled for a desired trench filling process. From the calculation results it is clear that a high deposition rate will result in undesired void formation during the trench filling, while a small deposition rate will result in undesired trench bottom and mask damage by sputtering. By varying the SiH4/O2 ratio, the chemical composition of the deposited layer will be influenced. However, even at the highest SiH4/O2 ratio investigated (i.e., 3.2:1; low oxygen content), the bulk deposited layer consists mainly of SiO2, suggesting that low-volatile silane species deposit first and subsequently become oxidized instead of being oxidized first in the plasma before deposition. Finally, it was found that the top surface of the deposited layer contained less oxygen due to preferential sputtering of O atoms, making the top layer more Si-rich. However, this effect is negligible at a SiH4/O2 ratio of 2:1 or lower.

Abstract: The reaction behavior of species that may affect the growth of ultrananocrystal line and nanocrystalline diamond ((U)NCD) films is investigated by means of molecular dynamics simulations. Impacts of CHx (x = 0 – 4), C2Hx (x=0-6), C3Hx (x=0-2), C4Hx (x = 0 – 2), H, and H-2 on clean and hydrogenated diamond (100)2 x 1 and (111) 1 x 1 surfaces at two different substrate temperatures are simulated. We find that the different bonding structures of the two surfaces cause different temperature effects on the sticking efficiency. These results predict a temperature-dependent ratio of diamond (100) and (111) growth. Furthermore, predictions of which are the most important hydrocarbon species for (U)NCD growth are made.

Abstract: A 2d3v Particle-in-cell/Monte Carlo collisions (PIC/MCC) model was constructed for an Ar/N2 reactive gas mixture in a magnetron discharge. A titanium target was used, in order to study the sputter deposition of a TiNx thin film. Cathode currents and voltages were calculated self-consistently and compared with experiments. Also, ion fluxes to the cathode were calculated, which cause sputtering of the target. The sputtered atom fluxes from the target, and to the substrate were calculated, in order to visualize the deposition of the TiNx film.

Abstract: In this review paper, we present several examples of reactive molecular dynamics simulations, which contribute to a better understanding of the underlying mechanisms in plasma medicine on the atomic scale. This includes the interaction of important reactive oxygen plasma species with the outer cell wall of both gram-positive and gram-negative bacteria, and with lipids present in human skin. Moreover, as most biomolecules are surrounded by a liquid biofilm, the behavior of these plasma species in a liquid (water) layer is presented as well. Finally, a perspective for future atomic scale modeling studies is given, in the field of plasma medicine in general, and for cancer treatment in particular.

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: We have characterized theoretically the work-function modifications of the (111) surfaces of gold and silver upon deposition of self-assembled monolayers based on methanethiol and trifluoromethanethiol. A comparative analysis is made between the experimental results and those obtained from two widely used approaches based on density functional theory. The contributions to the total work-function modifications are estimated on the basis of two decomposition schemes of the thiol molecules that have been proposed in the literature. The contributions are found to differ significantly between the two approaches, as do the corresponding adsorption energies.

Abstract: The catalytic oxidative dehydrogenation of hydrocarbons is an industrially important process, in which selectivity is a key issue. We here investigate the conversion of methanol to formaldehyde on a vanadia surface employing long timescale simulations, reaching a time scale of seconds. In particular, we compare the thermal process to the case where an additional external electric field is applied, as would be the case in a direct plasma-catalysis setup. We find that the electric field influences the retention time of the molecules at the catalyst surface. These simulations provide an atomic scale insight in the thermal catalytic oxidative dehydrogenation process, and in how an external electric field may affect this process.

Abstract: A hybrid Monte Carlofluid modeling network is developed for an argon-hydrogen mixture, to predict the effect of small amounts of hydrogen added to a dc argon glow discharge. The species considered in the model include the Ar gas atoms, electrons, Ar+ ions and fast Ar atoms, ArH+, H+, H+2 and H+3 ions, and H atoms and H2 molecules, as well as Ar metastable atoms, sputtered Cu atoms, and the corresponding Cu+ ions. Sixty-five reactions between these species are incorporated in the model. The effect of hydrogen on various calculation results is investigated, such as the species densities, the relative role of different production and loss processes for the various species, the cathode sputtering rate and contributions by different bombarding species, and the dissociation degree of H2 and the ionization degree of Ar and Cu. The calculation results are presented and discussed for 1% H2 addition, and comparison is also made with a pure argon discharge and with only 0.1% H2 addition.

Abstract: The electron distribution function (EDF) in an electron cyclotron resonance (ECR) discharge is far from Maxwellian. The self-consistent simulation of ECR discharges requires a calculation of the EDF on every magnetic line for various ion density profiles. The straightforward self-consistent simulation of ECR discharges using the Monte Carlo technique for the EDF calculation is very computer time expensive, since the electron and ion time scales are very different. An electron Boltzmann kinetic equation averaged over the fast electron bouncing and pitch-angle scattering was derived in order to develop an effective and operative tool for the fast modeling (FM) of low-pressure ECR discharges. An analytical solution for the EDF in a loss cone was derived. To check the validity of the FM, one-dimensional (in coordinate) and two-dimensional (in velocity) Monte Carlo simulation codes were developed. The validity of the fast modeling method is proved by comparison with the Monte Carlo simulations. The complete system of equations for FM is presented and ready for use in a comprehensive study of ECR discharges. The variations of plasma density and of wall and sheath potentials are analyzed by solving a self-consistent set of equations for the EDF.

Abstract: Electron acceleration from the interaction of an intense short-pulse laser with low density plasma is considered. The relation between direct electron acceleration within the laser pulse and that in the wake is investigated analytically. The magnitude and location of the ponderomotive-force-caused charge separation field with respect to that of the pulse determine the relative effectiveness of the two acceleration mechanisms. It is shown that there is an optimum condition for acceleration in the wake. Electron acceleration within the pulse dominates as the pulse becomes sufficiently short, and the latter directly drives and even traps the electrons. The latter can reach ultrahigh energies and can be extracted by impinging the pulse on a solid target. (C) 2003 American Institute of Physics.