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: Ultrathin bismuth titanate films (Bi2Ti2O7, 5-25 nm) are deposited onto SiO2/Si substrates by aqueous chemical solution deposition and their evolution during annealing is studied. The films crystallize into a preferentially oriented, pure pyrochlore phase between 500 and 700 degrees C, depending on the film thickness and the total thermal budget. Crystallization causes a strong increase of surface roughness compared to amorphous films. An increase of the interfacial layer thickness is observed after anneal at 600 degrees C, together with intermixing of bismuth with the substrate as shown by TEM-EDX. The band gap was determined to be similar to 3 eV from photoconductivity measurements and high dielectric constants between 30 and 130 were determined from capacitance voltage measurements, depending on the processing conditions. (C) 2012 Elsevier Ltd. All rights reserved.

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: Simulation results of electrostatics in Si cylindrical junctionless nanowire transistors with a homogenous channel are presented. Junctionless transistors including strain and arbitrary crystallographic orientations are studied. Size quantization effects are simulated by self-consistent solutions of the Poisson and Schrodinger equations. The 6 x 6 k.p method is employed for the calculation of the valence subband structure in a junctionless nanowire pFET. The influence of stress/strain and crystallographic channel orientation on to the electrostatics in terms of subband structure, charge density, and C-V curve is systematically studied. (C) 2011 Elsevier Ltd. All rights reserved.

Abstract: The electronic properties of two-dimensional honeycomb structures of molybdenum disulfide (MoS(2)) subjected to biaxial strain have been investigated using first-principles calculations based on density functional theory. On applying compressive or tensile bi-axial strain on bi-layer and mono-layer MoS(2), the electronic properties are predicted to change from semiconducting to metallic. These changes present very interesting possibilities for engineering the electronic properties of two-dimensional structures of MoS(2).

Abstract: Using reactive molecular dynamics simulations by means of the ReaxFF potential, we studied the growth mechanism of ultrathin silica (SiO2) layers during hyperthermal oxidation at room temperature. Oxidation of Si(100){2 × 1} surfaces by both atomic and molecular oxygen was investigated in the energy range 15 eV. The oxidation mechanism, which differs from thermal oxidation, is discussed. In the case of oxidation by molecular O2, silica is quickly formed and the thickness of the formed layers remains limited compared to oxidation by atomic oxygen. The Si/SiO2 interfaces are analyzed in terms of partial charges and angle distributions. The obtained structures of the ultrathin SiO2 films are amorphous, including some intrinsic defects. This study is important for the fabrication of silica-based devices in the micro- and nanoelectronics industry, and more specifically for the fabrication of metal oxide semiconductor devices.

Abstract: The growth mechanism and chirality formation of a single-walled carbon nanotube (SWNT) on a surface-bound nickel nanocluster are investigated by hybrid reactive molecular dynamics/force-biased Monte Carlo simulations. The validity of the interatomic potential used, the so-called ReaxFF potential, for simulating catalytic SWNT growth is demonstrated. The SWNT growth process was found to be in agreement with previous studies and observed to proceed through a number of distinct steps, viz., the dissolution of carbon in the metallic particle, the surface segregation of carbon with the formation of aggregated carbon clusters on the surface, the formation of graphitic islands that grow into SWNT caps, and finally continued growth of the SWNT. Moreover, it is clearly illustrated in the present study that during the growth process, the carbon network is continuously restructured by a metal-mediated process, thereby healing many topological defects. It is also found that a cap can nucleate and disappear again, which was not observed in previous simulations. Encapsulation of the nanoparticle is observed to be prevented by the carbon network migrating as a whole over the cluster surface. Finally, for the first time, the chirality of the growing SWNT cap is observed to change from (11,0) over (9,3) to (7,7). It is demonstrated that this change in chirality is due to the metal-mediated restructuring process.

Abstract: Germanium combined with high-κ dielectrics is investigated for the next generations of CMOS devices. Therefore, we study reaction mechanisms for Al2O3 atomic layer deposition on sulfur passivated Ge using calculations based on density functional theory and total reflection X-ray fluorescence (TXRF). TXRF indicates 6 S/nm2 and 4 Al/nm2 after the first TMA/H2O reaction cycle, and growth inhibition from the second reaction cycle on. Calculations are performed on molecular clusters representing −GeSH surface sites. The calculations confirm that the TMA reaction does not affect the S content. On fully SH-terminated Ge, TMA favorably reacts with up to three −GeSH sites, resulting in a near tetrahedral Al coordination. Electron deficient structures with a GeS site shared between two Al atoms are proposed. The impact of the cluster size on the structures and reaction energetics is systematically 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: Nonlinear electron plasma waves driven by a finite-charged particle pulse or rigid object moving at relativistic speeds are investigated. Quasi-stationary smooth and spiky wake waves comoving with the object are found. Localized soliton-like solutions are also shown to exist. Relativistic effects tend to prevent their formation because of the electron mass increase. The application of the very-large-amplitude wake density waves as a source of ultrahigh-energy cosmic-ray events is discussed.

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: The electronic properties of hydrogenated silicene and germanene, so called silicane and germanane, respectively, are investigated using first-principles calculations based on density functional theory. Two different atomic configurations are found to be stable and energetically degenerate. Upon the adsorption of hydrogen, an energy gap opens in silicene and germanene. Their energy gaps are next computed using the HSE hybrid functional as well as the G(0)W(0) many-body perturbation method. These materials are found to be wide band-gap semiconductors, the type of gap in silicane (direct or indirect) depending on its atomic configuration. Germanane is predicted to be a direct-gap material, independent of its atomic configuration, with an average energy gap of about 3.2 eV, this material thus being potentially interesting for optoelectronic applications in the blue/violet spectral range. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3595682]

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.