Abstract: We evaluate the importance of quantum effects on the sticking of H2 on Al(110) for conditions that are close to those of molecular beam experiments that have been done on this system. Calculations with the quasi-classical trajectory (QCT) method and with quantum dynamics (QD) are performed using a model in which only motion in the six molecular degrees of freedom is allowed. The potential energy surface used has a minimum barrier height close to the value recently obtained with the quantum Monte Carlo method. Monte Carlo averaging over the initial rovibrational states allowed the QD calculations to be done with an order of magnitude smaller computational expense. The sticking probability curve computed with QD is shifted to lower energies relative to the QCT curve by 0.21 to 0.05 kcal/mol, with the highest shift obtained for the lowest incidence energy. Quantum effects are therefore expected to play a small role in calculations that would evaluate the accuracy of electronic structure methods for determining the minimum barrier height to dissociative chemisorption for H2 + Al(110) on the basis of the standard procedure for comparing results of theory with molecular beam experiments.

Abstract: The gliding arc plasmatron (GAP) is a highly efficient atmospheric plasma source, which is very promising for CO2 conversion applications. To understand its operation principles and to improve its application, we present here comprehensive modeling results, obtained by means of computational fluid dynamics simulations and plasma modeling. Because of the complexity of the CO2 plasma, a full 3D plasma model would be computationally impractical. Therefore, we combine a 3D turbulent gas flow model with a 2D plasma and gas heating model in order to calculate the plasma parameters and CO2 conversion characteristics. In addition, a complete 3D gas flow and plasma model with simplified argon chemistry is used to evaluate the gliding arc evolution in space and time. The calculated values are compared with experimental data from literature as much as possible in order to validate the model. The insights obtained in this study are very helpful for improving the application of CO2 conversion, as they allow us to identify the limiting factors in the performance, based on which solutions can be provided on how to further improve the capabilities of CO2 conversion in the GAP.

Abstract: We present a first-principles investigation on the stability, electronic structure, and mechanical response of ultrathin heterostructures composed of single layers of InSe and SiGe. First, by performing total energy optimization and phonon calculations, we show that single layers of InSe and SiGe can form dynamically stable heterostructures in 12 different stacking types. Valence and conduction band edges of the heterobilayers form a type-I heterojunction having a tiny band gap ranging between 0.09 and 0.48 eV. Calculations on elastic-stiffness tensor reveal that two mechanically soft single layers form a heterostructure which is stiffer than the constituent layers because of relatively strong interlayer interaction. Moreover, phonon analysis shows that the bilayer heterostructure has highly Raman active modes at 205.3 and 43.7 cm(-1), stemming from the out-of-plane interlayer mode and layer breathing mode, respectively. Our results show that, as a stable type-I heterojunction, ultrathin heterobilayer of InSe/SiGe holds promise for nanoscale device applications.

Abstract: The phase transition processes induced by ultrashort, 100 fs pulsed laser irradiation of Au, Cu, and Ni are studied by means of a combined atomistic-continuum approach. A moderately low absorbed laser fluence range, from 200 to 600 J/m2 is considered to study phase transitions by means of a local and a nonlocal order parameter. At low laser fluences, the occurrence of layer-by-layer evaporation has been observed, which suggests a direct solid to vapor transition. The calculated amount of molten material remains very limited under the conditions studied, especially for Ni. Therefore, our results show that a kinetic equation that describes a direct solid to vapor transition might be the best approach to model laser-induced phase transitions by continuum models. Furthermore, the results provide more insight into the applicability of analytical superheating theories that were implemented in continuum models and help the understanding of nonequilibrium phase transitions.

Abstract: Morphological characterization of nanoparticle assemblies and hybrid nanomaterials is critical in determining their structure-property relationships as well as in the development of structures with desired properties. Electron tomography has become a widely utilized technique for the three-dimensional characterization of nanoparticle assemblies. However, the extraction of quantitative morphological parameters from the reconstructed volume can be a complex and labor-intensive task. In this study, we aim to overcome this challenge by automating the volumetric segmentation process applied to three-dimensional reconstructions of nanoparticle assemblies. The key to enabling automated characterization is to assess the performance of different volumetric segmentation methods in accurately extracting predefined quantitative descriptors for morphological characterization. In our methodology, we compare the quantitative descriptors obtained through manual segmentation with those obtained through automated segmentation methods, to evaluate their accuracy and effectiveness. To show generality, our study focuses on the characterization of assemblies of CdSe/CdS quantum dots, gold nanospheres and CdSe/CdS encapsulated in polymeric micelles, and silica-coated gold nanorods decorated with both CdSe/CdS or PbS quantum dots. We use two unsupervised segmentation algorithms: the watershed transform and the spherical Hough transform. Our results demonstrate that the choice of automated segmentation method is crucial for accurately extracting the predefined quantitative descriptors. Specifically, the spherical Hough transform exhibits superior performance in accurately extracting quantitative descriptors, such as particle size and interparticle distance, thereby allowing for an objective, efficient, and reliable volumetric segmentation of complex nanoparticle assemblies.

Abstract: Zinc(II)-phthalocyanine tetrasulfonic acid (Zn-PcS) molecules physisorbed on bilayer epitaxial graphene on silicon carbide (SiC(0001)) were studied by using scanning tunneling microscopy/spectroscopy (STM/STS) and density functional theory (DFT). Two different methods were used to deposit Zn-PcS molecules and regardless of the method being used, the surface coverage stayed very low indicating the weakness of surface-molecule interaction. STS measurements revealed that derivative of tunneling current with respect to voltage (dI/dV) measured on Zn-PcS molecules did not exhibit the characteristic dip observed on dI/dV curves of pristine bilayer epitaxial graphene. DFT calculations show that the energy of the lowest unoccupied molecular orbital (LUMO) of the Zn-PcS molecule is below the Dirac point of graphene which enhances local density of states (LDOS). We attribute the disappearance of the dip in the dI/dV curves measured on the Zn-PcS/bilayer system to the LUMO of Zn-PcS. Charge density calculations along Zn-PcS/graphene interface reveal that there is a small charge transfer from graphene to the molecule. Calculated adsorption energy (3.13 eV) of the molecule is notably low and is consistent with the observed low surface coverage at room temperature.

Abstract: Printed chalcopyrite thin films have attracted considerable attention in recent years due to their potential in the high-throughput production of photovoltaic devices. To improve the homogeneity of printed CuInSe2 (CISe) layers, chemical additives such as binder can be added to the precursor ink. In this contribution, we investigate the influence of the dicyandiamide (DCDA) content, used as a binder in the precursor ink, on the physical and electrical properties of printed CISe solar cells. It is shown that the use of the binder leads to a dense absorber, composed of large CISe grains close to the surface, while the bulk of the layer consists of CISe crystallites embedded in a CuxS particle based matrix, resulting from the limited sintering of the precursor in this region. The expected additional carbon contamination of the CISe layer due to the addition of the binder appears to be limited, and the optical properties of the CISe layer are similar to the reference sample without additive. The electrical characterization of the corresponding CISe/CdS solar cells shows a degradation of the efficiency of the devices, due to a modification in the predominant recombination mechanisms and a limitation of the space charge region width when using the binder; both effects could be explained by the inhomogeneity of the bulk of the CISe absorber and high defect density at the CISe/CuxS-based matrix interface.

Abstract: We study the mechanisms of CO2 conversion in a nanosecond repetitively pulsed (NRP) discharge, by means of a chemical kinetics model. The calculated conversions and energy efficiencies are in reasonable agreement with experimental results over a wide range of specific energy input values, and the same applies to the evolution of gas temperature and CO2 conversion as a function of time in the afterglow, indicating that our model provides a realistic picture of the underlying mechanisms in the NRP discharge and can be used to identify its limitations and thus to suggest further improvements. Our model predicts that vibrational excitation is very important in the NRP discharge, explaining why this type of plasma yields energy-efficient CO2 conversion. A significant part of the CO2 dissociation occurs by electronic excitation from the lower vibrational levels toward repulsive electronic states, thus resulting in dissociation. However, vibration−translation (VT) relaxation (depopulating the higher vibrational levels) and CO + O recombination (CO + O + M → CO2 + M), as well as mixing of the converted gas with fresh gas entering the plasma in between the pulses, are limiting factors for the conversion and energy efficiency. Our model predicts that extra cooling, slowing down the rate of VT relaxation and of the above recombination reaction, thus enhancing the contribution of the highest vibrational levels to the overall CO2 dissociation, can further improve the performance of the NRP discharge for energy-efficient CO2 conversion.

Abstract: Molecular dynamics (MD) simulations were performed to provide atomic scale insight in the initial interaction between hydroxyl radicals (OH) and peptide systems in solution. These OH radicals are representative reactive oxygen species produced by cold atmospheric plasmas. The use of plasma for biomedical applications is gaining increasing interest, but the fundamental mechanisms behind the plasma modifications still remain largely elusive. This study helps to gain more insight in the underlying mechanisms of plasma medicine but is also more generally applicable to peptide oxidation, of interest for other applications. Combining both reactive and nonreactive MD simulations, we are able to elucidate the reactivity of the amino acids inside the peptide systems and their effect on their structure up to 1 μs. Additionally, experiments were performed, treating the simulated peptides with a plasma jet. The computational results presented here correlate well with the obtained experimental data and highlight the importance of the chemical environment for the reactivity of the individual amino acids, so that specific amino acids are attacked in higher numbers than expected. Furthermore, the long time scale simulations suggest that a single oxidation has an effect on the 3D conformation due to an increase in hydrophilicity and intra- and intermolecular interactions.

Abstract: We study the oxygen thermal stability of two binary
systems. The larger particles are magnetic amorphous Co (7.2 nm) or
Fe3O4 (7.5 nm) nanocrystals, whereas the smaller ones (3.7 nm) are
Au nanocrystals. The nanocrystal ordering as well as the choice of the
magnetic nanoparticles very much influence the stability of the binary
system. A perfect crystalline structure is obtained with the Fe3O4/Au
binary supracrystals. For the Co/Au binary system, oxidation of Co
results in the chemical transformation from Co to CoO, where the size
of the amorphous Co nanoparticles increases from 7.2 to 9.8 nm in
diameter. During the volume expansion of the Co nanoparticles, Au
nanoparticles within the binary assemblies coalesce and are at the
origin of the instability of the binary nanoparticle supracrystals. On the
other hand, for the Fe3O4/Au binary system, the oxidation of Fe3O4 to
γ-Fe2O3 does not lead to a size change of the nanoparticles, which
maintains the stability of the binary nanoparticle supracrystals. A similar behavior is observed for an AlB2-type Co−Ag binary
system: The crystalline structure is maintained, whereas in disordered assemblies, coalescence of Ag nanocrystals is observed.

Abstract: Because of its excellent properties in nanotechnology applications, atomic layer deposition of ruthenium (Ru) has been the subject of numerous experimental studies. Recently, two different Ru precursors were compared for plasma-enhanced atomic layer deposition (PEALD) of Ru, and their reactivity was found to be different. Inhibition was observed for bis(ethylcyclopentadienyl)ruthenium (Ru(EtCp)(2)), while nearly linear growth behavior was observed for (methylcyclopentadienyl-pyrrolyl)ruthenium (Ru(MeCp)Py). To understand this difference in reactivity, we investigate the adsorption of RuCp, and RuCpPy (i.e., without substituents) on a TiN surface using calculations based on periodic boundary conditions density functional theory (DFT) combined with experiments based on Rutherford backscattering spectroscopy (RBS). The calculations demonstrate that the RuCpPy precursor chemisorbs on the TiN(100) surface while the RuCp2 precursor only physisorbs. We propose a reaction mechanism for the chemisorption of RuCpPy. The area density of the calculated RuCpPy surface species is compared with the experimental values from RBS. The impact of a H-plasma is also investigated. The DFT calculations and experimental results from RBS provide insight into the adsorption processes of the RuCpPy and RuCp2 precursors on the TiN(100) surface.

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: Atomic, layer deposition,(ALD of ruthenium using two ruthenium precursors, i.e., Ru(C5H5)(2) (RuCp2) and Ru(C5H5)(C4H4N) (RuCpPy), is studied using density functional theory. By investigating the reaction mechanisms On bare ruthenium surfaces, i.e., (001), (101), and (100), and H-terminated surfaces, an atomistic insight in the Ru ALD is provided. The calculated results show that on the Ru surfaces both RuCp2 and RuCpPy an undergo dehydrogenation and ligand dissociation reactions. RuCpPy is more reactive than RuCp2. By forming a, strong, bond between N of Py and Ru of the surface, RuCpPy can easily chemisorb on the surfaces. The reactions of RuCp2,On the Surfaces are less favorable the adsorption is not strong enough This could be a,factor contributing to the higher growth-per-cycle of Ru using RuCpPy, as observed experimentally. By Studying, the adsorption on H-terminated Ru surfaces, We showed that H Can prevent the adsorption of the precursors, thus inhibiting the growth of Ru. Our calculations indicate that the H content on the surface can have an impact on the growth-per-cycle. Finally, our simulations also demonstrate large impacts of the surface structure on the reaction mechanisms. Of the three surfaces, the (100) surface, which is the less stable and has a zigzag surface structure, is also the most reactive one.

Abstract: Cyrogenic etching of silicon is envisaged to enable better control over plasma processing in the microelectronics industry, albeit little is known about the fundamental differences compared to the room temperature process. We here present molecular dynamics simulations carried out to obtain sticking probabilities, thermal desorption rates, surface diffusion speeds, and sputter yields of F, F2, Si, SiF, SiF2, SiF3, SiF4, and the corresponding ions on Si(100) and on SiF13 surfaces, both at cryogenic and near room temperature. The different surface behavior during conventional etching and cryoetching is discussed. F2 is found to be relatively reactive compared to other species like SiF03. Thermal desorption occurs at a significantly lower rate under cryogenic conditions, which results in an accumulation of physisorbed species. Moreover, ion incorporation is often observed for ions with energies of 30400 eV, which results in a relatively low net sputter yield. The obtained results suggest that the actual etching of Si, under both cryogenic and near room temperature conditions, is based on the complete conversion of the Si surface to physisorbed SiF4, followed by subsequent sputtering of these molecules, instead of direct sputtering of the SiF03 surface.

Abstract: The adsorption of C and CHx radicals on anatase (001) was studied using DFT within the generalized gradient approximation using the Perde-Burke-Ernzerhof (PBE) functional. We have studied the influence of oxygen vacancies in and at the surface on the adsorption properties of the radicals. For the oxygen vacancies in anatase (001), the most stable vacancy is located at the surface. For this vacancy, the maximal adsorption strength of C and CH decreases compared to the adsorption on the stoichiometric surface, but it increases for CH2 and CH3. If an oxygen vacancy is present in the first subsurface layer, the maximal adsorption strength increases for C, CH, CH2, and CH3. When the vacancy is present in the next subsurface layer, we find that only the CH3 adsorption is enhanced, while the maximal adsorption energies for the other radical species decrease. Not only does the precise location of the oxygen vacancy determine the maximal adsorption interaction, it also influences the adsorption strengths of the radicals at different surface configurations. This determines the probability of finding a certain adsorption configuration at the surface, which in turn influences the possible surface reactions. We find that C preferentially adsorbs far away from the oxygen vacancy, while CH2 and CH3 adsorb preferentially at the oxygen vacancy site. A fraction of CH partially adsorbs at the oxygen vacancy, and another fraction adsorbs further away from the vacancy.