Diagnostics and two-dimensional simulation of low-frequency inductively coupled plasmas with neutral gas heating and electron heat fluxes

This article presents the results on the diagnostics and numerical modeling of low-frequency (∼460 KHz) inductively coupled plasmas generated in a cylindrical metal chamber by an external flat spiral coil. Experimental data on the electron number densities and temperatures, electron energy distribution functions, and optical emission intensities of the abundant plasma species in low/intermediate pressure argon discharges are included. The spatial profiles of the plasma density, electron temperature, and excited argon species are computed, for different rf powers and working gas pressures, using the two-dimensional fluid approach. The model allows one to achieve a reasonable agreement between the computed and experimental data. The effect of the neutral gas temperature on the plasma parameters is also investigated. It is shown that neutral gas heating (at rf powers≥0.55kW) is one of the key factors that control the electron number density and temperature. The dependence of the average rf power loss, per electron-ion pair created, on the working gas pressure shows that the electron heat flux to the walls appears to be a critical factor in the total power loss in the discharge.

Surface fluxes of Si and C adatoms at initial growth stages of SiC quantum dots

Self-assembly of highly stoichiometric SiC quantum dots still remains a major challenge for the gas/plasma-based nanodot synthesis. By means of a multiscale hybrid numerical simulation of the initial stage (0.1-2.5 s into the process) of deposition of SiCSi (100) quantum dot nuclei, it is shown that equal Si and kst atom deposition fluxes result in strong nonstoichiometric nanodot composition due to very different surface fluxes of Si and C adatoms to the quantum dots. At this stage, the surface fluxes of Si and C adatoms to SiC nanodots can be effectively controlled by manipulating the SiC atom influx ratio and the Si (100) surface temperature. It is demonstrated that at a surface temperature of 800 K the surface fluxes can be equalized after only 0.05 s into the process; however, it takes more then 1 s at a surface temperature of 600 K. Based on the results of this study, effective strategies to maintain a stoichiometric ([Si] [C] =1:1) elemental ratio during the initial stages of deposition of SiCSi (100) quantum dot nuclei in a neutral/ionized gas-based process are proposed.

Plasma-deposited Ge nanoisland films on Si : is Stranski–Krastanow fragmentation unavoidable?

The formation of Ge quantum dot arrays by deposition from a low-temperature plasma environment is investigated by kinetic Monte Carlo numerical simulation. It is demonstrated that balancing of the Ge influx from the plasma against surface diffusion provides an effective control of the surface processes and can result in the formation of very small densely packed quantum dots. In the supply-controlled mode, a continuous layer is formed which is then followed by the usual Stranski-Krastanow fragmentation with a nanocluster size of 10 nm. In the diffusion-controlled mode, with the oversupply relative to the surface diffusion rate, nanoclusters with a characteristic size of 3 nm are formed. Higher temperatures change the mode to supply controlled and thus encourage formation of the continuous layer that then fragments into an array of large size. The use of a high rate of deposition, easily accessible using plasma techniques, changes the mode to diffusion controlled and thus encourages formation of a dense array of small nanoislands.

Tailoring the composition of self-assembled Si1−xCx quantum dots : simulation of plasma/ion-related controls

Precise control of composition and internal structure is essential for a variety of novel technological applications which require highly tailored binary quantum dots (QDs) with predictable optoelectronic and mechanical properties. The delicate balancing act between incoming flux and substrate temperature required for the growth of compositionally graded (Si1-xC x; x varies throughout the internal structure), core-multishell (discrete shells of Si and C or combinations thereof) and selected composition (x set) QDs on low-temperature plasma/ion-flux-exposed Si(100) surfaces is investigated via a hybrid numerical simulation. Incident Si and C ions lead to localized substrate heating and a reduction in surface diffusion activation energy. It is shown that by incorporating ions in the influx, a steady-state composition is reached more quickly (for selected composition QDs) and the composition gradient of a Si1-xCx QD may be fine tuned; additionally (with other deposition conditions remaining the same), larger QDs are obtained on average. It is suggested that ionizing a portion of the influx is another way to control the average size of the QDs, and ultimately, their internal structure. Advantages that can be gained by utilizing plasma/ion-related controls to facilitate the growth of highly tailored, compositionally controlled quantum dots are discussed as well.

Nanosphere monolayer-templated, ion-assisted nanofeature etching in dielectric materials : a numerical simulation of nanoscale ion flux topography

The results of numerical simulations of nanometer precision distributions of microscopic ion fluxes in ion-assisted etching of nanoscale features on the surfaces of dielectric materials using a self-assembled monolayer of spherical nanoparticles as a mask are presented. It is shown that the ion fluxes to the substrate and nanosphere surfaces can be effectively controlled by the plasma parameters and the external bias applied to the substrate. By proper adjustment of these parameters, the ion flux can be focused onto the areas uncovered by the nanospheres. Under certain conditions, the ion flux distributions feature sophisticated hexagonal patterns, which may lead to very different nanofeature etching profiles. The results presented are generic and suggest viable ways to overcome some of the limitations of the existing plasma-assisted nanolithography.

Multi-scale hybrid numerical simulation of the growth of high-aspect-ratio nanostructures

Recently, a variety high-aspect-ratio nanostructures have been grown and profiled for various applications ranging from field emission transistors to gene/drug delivery devices. However, fabricating and processing arrays of these structures and determining how changing certain physical parameters affects the final outcome is quite challenging. We have developed several modules that can be used to simulate the processes of various physical vapour deposition systems from precursor interaction in the gas phase to gas-surface interactions and surface processes. In this paper, multi-scale hybrid numerical simulations are used to study how low-temperature non-equilibrium plasmas can be employed in the processing of high-aspect-ratio structures such that the resulting nanostructures have properties suitable for their eventual device application. We show that whilst using plasma techniques is beneficial in many nanofabrication processes, it is especially useful in making dense arrays of high-aspect-ratio nanostructures.

Nonlinear geometric and material computational technique : higher-order element with refined plastic hinge approach

This paper addresses of the advanced computational technique of steel structures for both simulation capacities simultaneously; specifically, they are the higher-order element formulation with element load effect (geometric nonlinearities) as well as the refined plastic hinge method (material nonlinearities). This advanced computational technique can capture the real behaviour of a whole second-order inelastic structure, which in turn ensures the structural safety and adequacy of the structure. Therefore, the emphasis of this paper is to advocate that the advanced computational technique can replace the traditional empirical design approach. In the meantime, the practitioner should be educated how to make use of the advanced computational technique on the second-order inelastic design of a structure, as this approach is the future structural engineering design. It means the future engineer should understand the computational technique clearly; realize the behaviour of a structure with respect to the numerical analysis thoroughly; justify the numerical result correctly; especially the fool-proof ultimate finite element is yet to come, of which is competent in modelling behaviour, user-friendly in numerical modelling and versatile for all structural forms and various materials. Hence the high-quality engineer is required, who can confidently manipulate the advanced computational technique for the design of a complex structure but not vice versa.

Simulation of gas-phase nanoparticle dynamics in the plasma-enhanced chemical vapor deposition of carbon nanostructures

The results of 1D simulation of nanoparticle dynamics in the areas adjacent to nanostructured carbon-based films exposed to chemically active complex plasma of CH4 + H2 + Ar gas mixtures are presented. The nanoparticle-loaded near-substrate (including sheath and presheath) areas of a low-frequency (0.5 MHz) inductively coupled plasma facility for the PECVD growth of the ordered carbon-based nanotip structures are considered. The conditions allowing one to predict the size of particles that can pass through the plasma sheath and softly land onto the surface are formulated. The possibility of soft nano-cluster deposition without any additional acceleration common for some existing nano-cluster deposition schemes is demonstrated. The effect of the substrate heating power and the average atomic mass of neutral species is studied numerically and verified experimentally.

Two-dimensional simulation of nanoassembly precursor species in Ar+H2+C2H2 reactive plasmas

The results of two-dimensional fluid simulation of number densities and fluxes of the main building blocks and surface preparation species involved in nanoassembly of carbon-based nanopatterns in Ar+H2+C2H2 reactive plasmas are reported. It is shown that the process parameters and non-uniformity of surface fluxes of each particular species may affect the targeted nanopattern quality. The results can be used to improve predictability of plasma-aided nanofabrication processes and optimize the parameters of plasma nanotools.KGaA, Weinheim.

Two-dimensional simulation of nanoparticle deposition from high-density plasmas on microstructured surfaces

Selective and controlled deposition of plasma-grown nanoparticles is one of the pressing problems of plasma-aided nanofabrication. The results of advanced numerical simulations of motion of charge-variable nanoparticles in the plasma presheath and sheath areas and in localized microscopic electric fields created by surface microstructures are reported. Conditions for site-selective deposition of such nanoparticles onto individual microstructures and open surface areas within a periodic micropattern are formulated. The effects of plasma parameters, surface potential, and micropattern features on nanoparticle deposition are investigated and explained using particle charging and plasma force arguments. The results are generic and applicable to a broad range of nanoparticle-generating plasmas and practical problems ranging from management of nanoparticle contamination in microelectronics to site-selective nanoparticle deposition into specified device locations, and synthesis of advanced microporous materials and nanoparticle superlattices. © 2007 American Institute of Physics.

Simulation of ion flux distribution in conductive and nonconductive nanotip patterns

The distribution of flux of carbon-bearing cations over nanopatterned surfaces with conductive nanotips and nonconductive nanoislands is simulated using the Monte-Carlo technique. It is shown that the ion current is focused to nanotip surfaces when the negative substrate bias is low and only slightly perturbed at higher substrate biases. In the low-bias case, the mean horizontal ion displacement caused by the nanotip electric field exceeds 10 nm. However, at higher substrate biases, this value reduces down to 2 nm. In the nonconductive nanopattern case, the ion current distribution is highly nonuniform, with distinctive zones of depleted current density around the nanoislands. The simulation results suggest the efficient means to control ion fluxes in plasma-aided nanofabrication of ordered nanopatterns, such as nanotip microemitter structures and quantum dot or nanoparticle arrays. © World Scientific Publishing Company.

Numerical simulation of self-organized nano-islands in plasma-based assembly of quantum dot arrays

This work presents the details of the numerical model used in simulation of self-organization of nano-islands on solid surfaces in plasma-assisted assembly of quantum dot structures. The model includes the near-substrate non-neutral layer (plasma sheath) and a nanostructured solid deposition surface and accounts for the incoming flux of and energy of ions from the plasma, surface temperature-controlled adatom migration about the surface, adatom collisions with other adatoms and nano-islands, adatom inflow to the growing nano-islands from the plasma and from the two-dimensional vapour on the surface, and particle evaporation to the ambient space and the two-dimensional vapour. The differences in surface concentrations of adatoms in different areas within the quantum dot pattern significantly affect the self-organization of the nano-islands. The model allows one to formulate the conditions when certain islands grow, and certain ones shrink or even dissolve and relate them to the process control parameters. Surface coverage by selforganized quantum dots obtained from numerical simulation appears to be in reasonable agreement with the available experimental results.

Numerical simulation of nanoparticle-generating electronegative plasmas in the PECVD of nanostructured silicon film

The results of numerical simulation of the equilibrium parameters of a low pressure nanopowder-generating discharge in silane for the plasma enhanced chemical vapor deposition (PECVD) of nanostructured silicon-based films are presented. It is shown that a low electron temperature and a low density of negative SiH3 - ions are favorable for the PECVD process. This opens a possibility to predict the main parameters of the reactive plasma and plasma-nucleated nanoparticles, and hence, to control the quality of silicon nanofilms.

A novel approach for numerical simulation of plant tissue shrinkage during drying

Drying is a key processing techniques used in food engineering which demands continual developments on advanced analysis techniques in order to optimize the product and the process. In this regard, plant based materials are a frequent subject of interest where microstructural studies can provide a clearer understanding on the fundamental physical mechanisms involved. In this context, considering numerous challenges of using conventional numerical grid-based modelling techniques, a meshfree particle based model was developed to simulate extreme deformations of plant microstructure during drying. The proposed technique is based on a particle based meshfree method: Smoothed Particle Hydrodynamics (SPH) and a Discrete Element Method (DEM). A tissue model was developed by aggrading individual cells modelled with SPH-DEM coupled approach by initializing the cells as hexagons and aggregating them to form a tissue. The model also involves a middle lamella resembling real tissues. Using the model, different dried tissue states were simulated with different moisture content, the turgor pressure, and cell wall contraction effects. Compared to the state of the art grid-based microscale plant tissue drying models, the proposed model is capable of simulating plant tissues at lower moisture contents which results in excessive shrinkage and cell wall wrinkling. Model predictions were compared with experimental findings and a fairly good agreement was observed both qualitatively and quantitatively.