Abstract: Using reactive molecular dynamics simulations, the melting behavior of nickelcarbon nanoclusters is examined. The phase diagrams of icosahedral and Wulff polyhedron clusters are determined using both the Lindemann index and the potential energy. Formulae are derived for calculating the equilibrium constants and the solid and liquid fractions during a phase transition, allowing more rational determination of the melting temperature with respect to the arbitrary Lindemann value. These results give more insight into the properties of nickelcarbon nanoclusters in general and can specifically be very useful for a better understanding of the synthesis of carbon nanotubes using the catalytic chemical vapor deposition method.

Abstract: In this computational study, a gliding arc plasma reactor with a reverse-vortex flow stabilization is modelled for the first time by a fluid plasma description. The plasma reactor operates with argon gas at atmospheric pressure. The gas flow is simulated using the k-ε Reynolds-averaged Navier–Stokes turbulent model. A quasi-neutral fluid plasma model is used for computing the plasma properties. The plasma arc movement in the reactor is observed, and the results for the gas flow, electrical characteristics, plasma density, electron temperature, and gas temperature are analyzed.

Abstract: We study the properties of a capacitive 13.56 MHz discharge with a mixture of Ar/C₂H₂ taking into account the plasmochemistry and growth of heavy hydrocarbons. A hybrid model was developed to combine the kinetic description for electron motion and the fluid approach for negative and positive ion transports and plasmochemical processes. A significant change in plasma parameters related to injection of 5.8% portion of acetylene in argon was observed and analyzed. We found that the electronegativity of the mixture is about 30%. The densities of negatively and positively charged heavy hydrocarbons are sufficiently large to be precursors for the formation of nanoparticles in the discharge volume.

Abstract: An inductively coupled plasma, connected to the sampling cone of a mass spectrometer, is computationally investigated. The occurrence of rotational motion of the auxiliary and carrier gas flows is studied. The effects of operating parameters, i.e., applied power and gas flow rates, as well as geometrical parameters, i.e., sampler orifice diameter and injector inlet diameter, are investigated. Our calculations predict that at higher applied power the auxiliary and carrier gas flows inside the torch move more forward to the sampling cone, which is validated experimentally for the auxiliary gas flow, by means of an Elan 6000 ICP-MS. Furthermore, an increase of the gas flow rates can also modify the occurrence of rotational motion. This is especially true for the carrier gas flow rate, which has a more pronounced effect to reduce the backward motion than the flow rates of the auxiliary and cooling gas. Moreover, a larger sampler orifice (e.g., 2 mm instead of 1 mm) reduces the backward flow of the auxiliary gas path lines. Finally, according to our model, an injector inlet of 2 mm diameter causes more rotations in the carrier gas flow than an injector inlet diameter of 1.5 mm, which can be avoided again by changing the operating parameters.

Abstract: We studied the transport of copper droplets through an inductively coupled plasma, connected to the sampling cone of a mass spectrometer, by means of a computational model. The sample droplets are followed until they become evaporated. They are inserted as liquid particles from the central inlet and the effects of injection position (i.e. “on” and “off” axis), droplet diameter, as well as mass loading flow rate are investigated. It is shown that more “on-axis” injection of the droplets leads to a more straight path line, so that the droplets move less in the radial direction and are evaporated more on the central axis, enabling a better sample transfer efficiency to the sampler cone. Furthermore, there are optimum ranges of diameters and flow rates, which guarantee the proper position of evaporation along the torch, i.e. not too early, so that the sample can get lost in the torch, and not too late, which reduces the chance of becoming ionized before reaching the sampler.

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: Molecular dynamics simulations using the Brenner potential have been performed to investigate reaction mechanisms of various hydrocarbon radicals with low kinetic energies on amorphous hydrogenated carbon (a-C:H) surfaces and to simulate thin a-C:H film growth. Experimental data from an expanding thermal plasma setup were used as input for the simulations. The hydrocarbon reaction mechanisms were studied both during growth of the films and on a set of surface sites specific for a-C:H surfaces. Thin film growth was studied using experimentally detected growth species. It is found that the reaction mechanisms and sticking coefficients are dependent on the specific surface sites, and the structural properties of the growth radicals. Furthermore, it is found that thin a-C:H films can be densified using an additional H-flux towards the substrate.

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 multiphoton 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.