Abstract: We present here some of our modeling efforts for PECVD growth of nanostructured carbon materials with focus on amorphous hydrogenated carbon. Experimental data from an expanding thermal plasma setup were used as input for the simulations. Attention was focused both on the film growth mechanism, as well as on the hydrocarbon reaction mechanisms during growth of the films. 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. The film growth results are in correspondence with the experiment. Furthermore, it is found that thin a-C:H films can be densified using an additional H-flux towards the substrate.

Abstract: Understanding the fundamental dynamics in carbon nanotube (CNT) catalysts is of primary importance to understand CNT nucleation. This Letter reports on calculated bond switching (BS) rates in pure and carbon containing nickel nanoclusters. The rates are analyzed in terms of their temperature dependent spatial distribution and the mobility of the cluster atoms. The BS mechanism is found to change from vibrational to diffusional at around 900 K, with a corresponding strong increase in activation energy. Furthermore, the BS activation energy is observed to decrease as the carbon content in the cluster increases, resulting in an effective liquification of the cluster.

Abstract: In this letter we compare the effect of a radio-frequency sine, a low frequency sine, a rectangular and a pulsed dc voltage profile on the calculated electron production and power consumption in the dielectric barrier discharge. We also demonstrate using calculated potential distribution profiles of high time and space resolution how the pulsed dc discharge generates a secondary discharge pulse by deactivating the power supply.

Abstract: Mesoporous carbon, which can be replicated from mesoporous silica and whose surface is hydrophobic, can be an ideal catalyst for the esterification of fatty acids. Here we report an easy and low cost way to prepare sulfonic acid group-functionalized mesoporous carbon. A sample of calcined mesoporous silica SBA-15 was added to an aqueous sucrose solution followed by drying and calcination at different temperatures. In contrast to existing procedures, the obtained hybrid Si/C material was then first sulfonated in H2SO4, before the final removal of the silica template in order to stabilize the porous structure towards the liquid phase sulfonation treatment. Thus the silicacarbon composites, instead of the mesoporous carbon, were successfully sulfonated to introduce SO3H groups, while keeping the ordered mesoporous structure intact. The influence of carbonization temperature was investigated, suggesting an optimum temperature of 873 K. The SO3H group-functionalized mesoporous carbon, denoted as CMK-3-873-SO3H, was characterized by means of XRD, N2 physisorption, SEM, FT-IR, elemental analysis and TEM. It followed that a uniform mesoporous carbon was obtained with an average pore size of 3.89 nm, a specific surface of 807 m2/g and a SO3H group loading of 0.39 meq/g of dry material. Compared with other solid acid catalysts, the resulting material shows enhanced activity in the acid-catalyzed esterification of oleic acid with methanol, and can be used repeatedly. The increased catalytic performance is attributed to the hydrophobic surface and larger pore size of the new catalyst. It can effectively accommodate long chain fatty acids and reject formed water, making the active sites easily accessible.

Abstract: In this paper, an overview is given of modeling activities going on in our research group, for describing the plasma chemistry and plasmasurface interactions in reactive plasmas. The plasma chemistry is calculated by a fluid approach or by hybrid Monte Carlo (MC)fluid modeling. An example of both is illustrated in the first part of the paper. The example of fluid modeling is given for a dielectric barrier discharge (DBD) in CH4/O2, to describe the partial oxidation of CH4 into value-added chemicals. The example of hybrid MCfluid modeling concerns an inductively coupled plasma (ICP) etch reactor in Ar/Cl2/O2, including also the description of the etch process. The second part of the paper deals with the treatment of plasmasurface interactions on the atomic level, with molecular dynamics (MD) simulations or a combination of MD and MC simulations.

Abstract: Different zeolites were tested for selective removal of methyl elaidate (trans isomer) from an equimolar mixture with methyl oleate (cis isomer). Sorption experiments of the geometric isomers show that only ZSM-5 samples with reduced Al content in the framework are able to discriminate among the bent cis and the linear trans fatty acid methyl esters. Hydrogenation experiments of equimolar methyl oleate and elaidate mixtures at low temperature (65 °C) and high hydrogen pressure (6.0 MPa), using Pt catalysts, confirm this result. Only with a Pt/NaZSM-5 catalyst outspoken selectivity for the hydrogenation of the trans isomer is obtained. In order to prepare a selective Pt/ZSM-5 catalyst, the influence of Pt addition (impregnation, ion-exchange and competitive ion-exchange) and Pt activation (different calcination and reduction temperatures) on the Pt-distribution and Pt particle size was investigated using SEM, bright-field and HR TEM, EDX, electron tomography, CO-chemisorption, XPS, XRD, and UVvis measurements. The best result in terms of hydrogenation activity and selectivity is obtained with a Pt/ZSM-5 catalyst, which is prepared via competitive ion-exchange, followed by slow calcination up to 350 °C under high O2 flow and a reduction up to 500 °C under H2. This preparation method leads to a Pt/ZSM-5 catalyst with the best Pt distribution and the smallest Pt clusters occluded in the zeolite structure. Finally, the influence of zeolite crystal size, morphology, and elemental composition of ZSM-5 on hydrogenation activity and selectivity was investigated in detail.

Abstract: Molecular dynamics simulations using improved TersoffBrenner potential parameters were performed to investigate Cl+ etching of a {2×1} reconstructed Si(100) surface. Steady-state Si etching accompanying the Cl coverage of the surface is observed. Furthermore, a steady-state chlorinated reaction layer is formed. The thickness of this reaction layer is found to increase with increasing energy. The stoichiometry of SiClx species in the reaction layer is found to be SiCl:SiCl2:SiCl3 = 1.0:0.14:0.008 at 50 eV. These results are in excellent agreement with available experimental data. While elemental Si products are created by physical sputtering, most SiClx (0<x<4) etch products are produced by chemical-enhanced physical sputtering.

Abstract: In this paper, we present the results of combined molecular dynamics−Metropolis Monte Carlo (MD-MMC) simulations of hydrocarbon species at flat diamond (100)2 × 1 and (111)1 × 1 surfaces. The investigated species are considered to be the most important growth species for (ultra)nanocrystalline diamond ((U)NCD) growth. When applying the MMC algorithm to stuck species at monoradical sites, bonding changes are only seen for CH2. The sequence of the bond breaking and formation as put forward by the MMC simulations mimics the insertion of CH2 into a surface dimer as proposed in the standard growth model of diamond. For hydrocarbon species attached to two adjacent radical (biradical) sites, the MMC simulations give rise to significant changes in the bonding structure. For UNCD, the combinations of C3 and C3H2, and C3 and C4H2 (at diamond (100)2 × 1) and C and C2H2 (at diamond (111)1 × 1) are the most successful in nucleating new crystal layers. For NCD, the following combinations pursue the diamond structure the best: C2H2 and C3H2 (at diamond (100)2 × 1) and CH2 and C2H2 (at diamond (111)1 × 1). The different behaviors of the hydrocarbon species at the two diamond surfaces are related to the different sterical hindrances at the diamond surfaces.

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.