Abstract: The recent advantages of the fabrication of artificial nanochannels enabled new research on the molecular transport, permeance, and selectivity of various gases and molecules. However, the physisorption/chemisorption of the unwanted molecules (usually hydrocarbons) inside nanochannels results in the alteration of the functionality of the nanochannels. We investigated contamination due to hydrocarbon molecules, nanochannels made of graphene, hexagonal boron nitride, BC2N, and molybdenum disulfide using molecular dynamics simulations. We found that for a certain size of nanochannel (i.e., h = 0.7 nm), as a result of the anomalous hydrophilic nature of nanochannels made of graphene, the hydrocarbons are fully adsorbed in the nanochannel, giving rise to full uptake. An increasing temperature plays an important role in unclogging, while pressure does not have a significant role. The results of our pioneering work contribute to a better understanding and highlight the important factors in alleviating the contamination and unclogging of nanochannels, which are in good agreement with the results of recent experiments.

Abstract: CsPbI3 is a promising perovskite material for photovoltaic applications in its photoactive perovskite or black phase. However, the material degrades to a photovoltaically inactive or yellow phase at room temperature. Various mitigation strategies are currently being developed to increase the lifetime of the black phase, many of which rely on inducing strains in the material that hinder the black-to-yellow phase transition. Physical insight into how these strategies exactly induce strain as well as knowledge of the spatial extent over which these strains impact the material is crucial to optimize these approaches but is still lacking. Herein, we combine machine learning potential-based molecular dynamics simulations with our in silico strain engineering approach to accurately quantify strained large-scale atomic structures on a nanosecond time scale. To this end, we first model the strain fields introduced by atomic substitutions as they form the most elementary strain sources. We demonstrate that the magnitude of the induced strain fields decays exponentially with the distance from the strain source, following a decay rate that is largely independent of the specific substitution. Second, we show that the total strain field induced by multiple strain sources can be predicted to an excellent approximation by summing the strain fields of each individual source. Finally, through a case study, we illustrate how this additive character allows us to explain how complex strain fields, induced by spatially extended strain sources, can be predicted by adequately combining the strain fields caused by local strain sources. Hence, the strain additivity proposed here can be adopted to further our insight into the complex strain behavior in perovskites and to design strain from the atomic level onward to enhance their sought-after phase stability.

Abstract: Molecular dynamics simulations are essential for a better understanding of dissociative chemisorption on metal surfaces, which is often the rate-controlling step in heterogeneous and plasma catalysis. The workhorse quasi-classical trajectory approach ubiquitous in molecular dynamics is able to accurately predict reactivity only for high translational and low vibrational energies. In contrast, catalytically relevant conditions generally involve low translational and elevated vibrational energies. Existing quantum dynamics approaches are intractable or approximate as a result of the large number of degrees of freedom present in molecule−metal surface reactions. Here, we extend a ring polymer molecular dynamics approach to fully include, for the first time, the degrees of freedom of a moving metal surface. With this approach, experimental sticking probabilities for the dissociative chemisorption of methane on Pt(111) are reproduced for a large range of translational and vibrational energies by including nuclear quantum effects and employing full-dimensional simulations.

Abstract: The resistance of a TiO2 anatase (001) surface to coke formation was studied in the context of dry reforming of methane using density functional theory (DFT) calculations. As carbon atoms act as precursors for coke formation, the resistance to coke formation can be measured by the carbon coverage of the surface. This is related to the stability of different CHx (x = 0−3) species and their rate of hydrogenation and dehydrogenation on the TiO2 surface. Therefore, we studied the reaction mechanisms and their corresponding rates as a function of the temperature for the dehydrogenation of the species on the surface. We found that the stabilities of C and CH are significantly lower than those of CH3 and CH2. The hydrogenation rates of the different species are significantly higher than the dehydrogenation rates in a temperature range of 300−1000 K. Furthermore, we found that dehydrogenation of CH3, CH2, and CH will only occur at appreciable rates starting from 600, 900, and 900 K, respectively. On the basis of these results, it is clear that the anatase (001) surface has a high coke resistance, and it is thus not likely that the surface will become poisoned by coke during dry reforming of methane. As the rate limiting step in dry reforming is the dissociative adsorption of CH4, we studied an alternative approach to thermal catalysis. We found that the temperature threshold for dry reforming is at least 700 K. This threshold temperature may be lowered by the use of plasma-catalysis, where the appreciable rates of adsorption of plasma-generated CHx radicals result in bypassing the rate limiting step of the reaction.

Abstract: A multiscale modeling and simulation approach, including first-principles calculations, ab initio molecular dynamics simulations, and a tight binding approach, is employed to study band flattening of the electronic band structure of oxidized monolayer graphene. The width offlat bands can be tuned by strain, the external electric field, and the density of functional groups and their distribution. A transition to a conducting state is found for monolayer graphene with impurities when it is subjected to an electric field of similar to 1.0 V/angstrom. Several parallel impurity-induced flat bands appear in the low-energy spectrum of monolayer graphene when the number of epoxy groups is changed. The width of the flat band decreases with an increase in tensile strain but is independent of the electric field strength. Here an alternative and easy route for obtaining band flattening in thermodynamically stable functionalized monolayer graphene is introduced. Our work discloses a new avenue for research on band flattening in monolayer graphene.

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: Grain boundaries between 60 degrees rotated and twinned crystals constitute the dominant type of extended line defects in two-dimensional transition metal dichalcogenides (2D MX2) when grown on a single crystalline template through van der Waals epitaxy. The two most common 60 degrees grain boundaries in MX2 layers, i.e., beta- and gamma-boundaries, introduce distinct distortion and strain into the 2D lattice. They impart a localized tensile or compressive strain on the subsequent layer, respectively, due to van der Waals coupling in bilayer MX2 as determined by combining atomic resolution electron microscopy, geometric phase analysis, and density functional theory. Based on these observations, an alternate route to strain engineering through controlling intrinsic van der Waals forces in homobilayer MX2 is proposed. In contrast to the commonly used external means, this approach enables the localized application of strain to tune the electronic properties of the 2D semiconducting channel in ultra-scaled nanoelectronic applications.

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: 3D characterization of assemblies of nanoparticles is of great importance to determine their structure-property connection. Such investigations become increasingly more challenging when the assemblies become larger and more compact. In this paper, we propose an optimized approach for electron tomography to minimize artefacts related to beam broadening in High Angle Annular Dark-Field Scanning Transmission Electron Microscopy mode. These artefacts are typically present at one side of the reconstructed 3D data set for thick nanoparticle assemblies. To overcome this problem, we propose a procedure in which two tomographic tilt series of the same sample are acquired. After acquiring the first series, the sample is flipped over 180o, and a second tilt series is acquired. By merging the two reconstructions, blurring in the reconstructed volume is minimized. Next, this approach is combined with an advanced three-dimensional reconstruction algorithm yielding quantitative structural information. Here, the approach is applied to a thick and compact assembly of spherical Au nanoparticles, but the methodology can we used to investigate a broad range of samples.

Abstract: To gain insight into the nature of the adhesion mechanism between hydroxyapatite (HA) and rutile (rTiO(2)), the mutual affinity between their surfaces was systematically studied using density functional theory (DFT). We calculated both bulk and surface properties of HA and rTiO(2), and explored the interfacial bonding mechanism of amorphous HA (aHA) surface onto amorphous as well as stoichiometric and nonstoichiometric crystalline rTiO(2). Formation energies of bridging and subbridging oxygen vacancies considered in the rTiO(2)(110) surface were evaluated and compared with other theoretical and experimental results. The interfacial interaction was evaluated through the work of adhesion. For the aHA/rTiO(2)(110) interfaces, the work of adhesion is found to depend strongly on the chemical environment of the rTiO(2)(110) surface. Electronic analysis indicates that the charge transfer is very small in the case of interface formation between aHA and crystalline rTiO(2)(110). In contrast, significant charge transfer occurs between aHA and amorphous rTiO(2) (aTiO(2)) slabs during the formation of the interface. Charge density difference (CDD) analysis indicates that the dominant interactions in the interface have significant covalent character, and in particular the Ti-O and Ca-O bonds. Thus, the obtained results reveal that the aHA/aTiO(2) interface shows a more preferable interaction and is thermodynamically more stable than other interfaces. These results are particularly important for improving the long-term stability of HA-based implants.

Abstract: Supersonic flows provide a high thermodynamic
nonequilibrium, which is crucial for energy-efficient conversion of
CO 2 in microwave plasmas and are therefore of great interest.
However, the effect of the flow on the chemical reactions is poorly
understood. In this work, we present a combined flow and plasma
chemical kinetics model of a microwave CO 2 plasma in a Laval
nozzle setup. The effects of the flow field on the different dissociation
and recombination mechanisms, the vibrational distribution, and the
vibrational transfer mechanism are discussed. In addition, the effect
of experimental parameters, like position of power deposition, outlet
pressure, and specific energy input, on the CO 2 conversion and
energy efficiency is examined. The short residence time of the gas in
the plasma region, the shockwave, and the maximum critical heat,
and thus power, that can be added to the flow to avoid thermal
choking are the main obstacles to reaching high energy efficiencies.

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.