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: Using density functional theory based first-principles calculations, the effect of surface hydrogenation on the structural, dynamical, electronic, and mechanical properties of monolayer washboard-gallenene (w-gallenene) is investigated. It is found that the dynamically stabilized strained monolayer of w-gallenene has a metallic nonmagnetic ground state. Both one-sided and two-sided hydrogenations of w-gallenene suppress its dynamical instability even when unstrained. Unlike one-sided hydrogenated monolayer w-gallenene (os-w-gallenene), two-sided hydrogenated monolayer w-gallenene (ts-w-gallenene) possesses the same crystal structure as w-gallenene. Electronic band structure calculations reveal that monolayers of hydrogenated derivatives of w-gallenene exhibit also metallic nonmagnetic ground state. Moreover, the linear-elastic constants, in-plane stiffness and Poisson ratio, are enhanced by hydrogenation, which is opposite to the behavior of other hydrogenated monolayer crystals. Furthermore, monolayer w-gallenene and ts-w-gallenene remain dynamically stable up to relatively higher biaxial strains as compared to borophene. With its enhanced dynamical stability, robust metallic character, and enhanced linear-elastic properties, hydrogenated monolayer w-gallenene is a potential candidate for nanodevice applications as a two-dimensional flexible metal.

Abstract: The in-plane stiffness (C), Poisson's ratio (nu), Young's modulus and ultimate strength (sigma) along two different crystallographic orientations are calculated for the single layer crystals: MoS2, HfS2 and TiS3 in 1H, 1T and monoclinic phases. We find that MoS2 and HfS2 have isotropic in-plane stiffnesses of 124.24 N m(-1) and 79.86 N m(-1), respectively. While for TiS3 the in-plane stiffness is highly anisotropic due to its monoclinic structure, with C-x = 83.33 N m(-1) and C-y = 133.56 N m(-1) (x and y are parallel to its longer and shorter in-plane lattice vectors.). HfS2 which is in the 1T phase has the smallest anisotropy in its ultimate strength, whereas TiS3 in the monoclinic phase has the largest. Along the armchair direction MoS2 has the largest sigma of 23.48 GPa, whereas along y TiS3 has the largest sigma of 18.32 GPa. We have further analyzed the band gap response of these materials under uniaxial tensile strain, and find that they exhibit different behavior. Along both armchair and zigzag directions, the band gap of MoS2 (HfS2) decreases (increases) as strain increases, and the response is almost isotropic. For TiS3, the band gap decreases when strain is along x, while if strain is along y, the band gap increases first and then decreases beyond a threshold strain value. The different characteristics observed in these sulphides with different structures shed light on the relationship between the structure and properties, which is useful for applications in nanotechnology.

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: Anatase nanosheets with exposed {001} facets
have gained increasing interest for photocatalytic applications. To
fully understand the structure-to-activity relation, combined
experimental and computational methods have been exploited.
Anatase nanosheets were prepared under hydrothermal conditions
in the presence of fluorine ions. High resolution scanning
transmission electron microscopy was used to fully characterize
the synthesized material, confirming the TiO2 nanosheet
morphology. Moreover, the surface structure and composition
of a single nanosheet could be determined by annular bright-field
scanning transmission electron microscopy (ABF-STEM) and
STEM electron energy loss spectroscopy (STEM-EELS). The photocatalytic activity was tested for the decomposition of organic
dyes rhodamine 6G and methyl orange and compared to a reference TiO2 anatase sample. The anatase nanosheets with exposed
{001} facets revealed a significantly lower photocatalytic activity compared to the reference. In order to understand the
mechanism for the catalytic performance, and to investigate the role of the presence of F−, light-induced electron paramagnetic
resonance (EPR) experiments were performed. The EPR results are in agreement with TEM, proving the presence of Ti3+
species close to the surface of the sample and allowing the analysis of the photoinduced formation of paramagnetic species.
Further, ab initio calculations of the anisotropic effective mass of electrons and electron holes in anatase show a very high effective
mass of electrons in the [001] direction, having a negative impact on the mobility of electrons toward the {001} surface and thus
the photocatalysis. Finally, motivated by the experimental results that indicate the presence of fluorine atoms at the surface, we
performed ab initio calculations to determine the position of the band edges in anatase slabs with different terminations of the
{001} surface. The presence of fluorine atoms near the surface is shown to strongly shift down the band edges, which indicates
another reason why it can be expected that the prepared samples with a large amount of {001} surface, but with fluorine atoms
near the surface, show only a low photocatalytic activity.

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: Plasma catalysis is gaining increasing interest for various environmental applications, but the crucial question is whether plasma can be created inside catalyst pores and under which conditions. In practice, various catalytic support materials are used, with various dielectric constants. We investigate here the influence of the dielectric constant on the plasma properties inside catalyst pores and in the sheath in front of the pores, for various pore sizes. The calculations are performed by a two-dimensional fluid model for an atmospheric pressure dielectric barrier discharge in helium. The electron impact ionization rate, electron temperature, electron and ion density, as well as the potential distribution and surface charge density, are analyzed for a better understanding of the discharge behavior inside catalyst pores. The results indicate that, in a 100 μm pore, the electron impact ionization in the pore, which is characteristic for the plasma generation inside the pore, is greatly enhanced for dielectric constants below 300. Smaller pore sizes only yield enhanced ionization for smaller dielectric constants, i.e., up to εr = 200, 150, and 50 for pore sizes of 50, 30, and 10 μm. Thus, the most common catalyst supports, i.e., Al2O3 and SiO2, which have dielectric constants around εr = 8−11 and 4.2, respectively, should allow more easily that microdischarges can be formed inside catalyst pores, even for smaller pore sizes. On the other hand, ferroelectric materials with dielectric constants above 300 never seem to yield plasma enhancement inside catalyst pores, not even for 100 μm pore sizes. Furthermore, it is clear that the dielectric constant of the material has a large effect on the extent of plasma enhancement inside the catalyst pores, especially in the range between εr = 4 and εr = 200. The obtained results are explained in detail based on the surface charge density at the pore walls,
and the potential distribution and electron temperature inside and above the pores. The results obtained with this model are
important for plasma catalysis, as the production plasma species in catalyst pores might affect the catalyst properties, and thus
improve the applications of plasma catalysis.

Abstract: We developed a microkinetic model to reveal the effects of plasma-generated radicals, intermediates, and vibrationally excited species on the catalytic hydrogenation of CO2 to CH3OH on a Cu(111) surface. As a benchmark, we first present the mechanisms of thermal catalytic CH3OH formation. Our model predicts that the reverse water-gas shift reaction followed by CO hydrogenation, together with the formate path, mainly contribute to CH3OH formation in thermal catalysis. Adding plasma-generated radicals and intermediates results in a higher CH3OH turnover frequency (TOF) by six to seven orders of magnitude, showing the potential of plasma-catalytic CO2 hydrogenation into CH3OH, in accordance with the literature. In addition, CO2 vibrational excitation further increases the CH3OH TOF, but the effect is limited due to relatively low vibrational temperatures under typical plasma catalysis conditions. The predicted increase in CH3OH formation by plasma catalysis is mainly attributed to the increased importance of the formate path. In addition, the conversion of plasma-generated CO to HCO* and subsequent HCOO* or H2CO* formation contribute to CH3OH formation. Both pathways bypass the HCOO* formation from CO2, which is the main bottleneck in the process. Hence, our model points toward the important role of CO, but also O, OH, and H radicals, as they influence the reactions that consume CO2 and CO. In addition, our model reveals that the H pressure should not be smaller than ca. half of the O pressure in the plasma as this would cause O* poisoning, which would result in very small product TOFs. Thus, plasma conditions should be targeted with a high CO and H content as this is favorable for CH3OH formation, while the O content should be minimized.

Abstract: Two-dimensional materials composed of transition metal carbides and nitrides (MXenes) are poised to revolutionize energy conversion and storage. In this work, we used density functional theory (DFT) to investigate the adsorption of Mg and Na adatoms on five M2CS2 monolayers (where M = Mo, Nb, Ti, V, and Zr) for battery applications. We assessed the stability of the adatom (i.e. Na and Mg)-monolayer systems by calculating adsorption and formation energies, as well as voltages as a function of surface coverage. For instance, we found that Mo2CS2 cannot support a full layer of Na nor even a single Mg atom. Na and Mg exhibit the strongest binding on Zr2CS2, followed by Ti2CS2, Nb2CS2 and V2CS2. Using the nudged elastic band method (NEB), we computed promising diffusion barriers for both dilute and nearly full ion surface coverage cases. In the dilute ion adsorption case, a single Mg and Na atom on Ti2CS2 experience similar to 0.47 eV and similar to 0.10 eV diffusion barriers between the lowest energy sites, respectively. For a nearly full surface coverage, a Na ion moving on Ti2CS2 experiences a similar to 0.33 eV energy barrier, implying a concentration-dependent diffusion barrier. Our molecular dynamics results indicate that the three (one) layers (layer) of the Mg (Na) ion on both surfaces of Ti2CS2 remain stable at T = 300 K. While, according to voltage calculations, Zr2CS2 can store Na up to three atomic layers, our MD simulations predict that the outermost layers detach from the Zr2CS2 monolayer due to the weak interaction between Na ions and the monolayer. This suggests that MD simulations are essential to confirm the stability of an ion-electrode system – an insight that is mostly absent in previous studies.

Abstract: We report on the ordering and dynamics of interacting colloidal particles confined by a parabolic potential. By means of Brownian dynamics simulations, we find that by varying the magnitude of the trap stiffness, it is possible to control the dimension of the system and, thus, explore both the structural transitions and the long-time self-diffusion coefficient as a function of the degree of confinement. We particularly study the structural ordering in the directions perpendicular and parallel to the confinement. Further analysis of the local distribution of the first-neighbors layer allows us to identify the different structural phases induced by the parabolic potential. These results are summarized in a structural state diagram that describes the way in which the colloidal suspension undergoes a structural re-ordering while increasing the confinement. To fully understand the particle dynamics, we take into account hydrodynamic interactions between colloids; the parabolic potential constricts the available space for the colloids, but it does not act on the solvent. Our findings show a non-linear behavior of the long-time self-diffusion coefficient that is associated to the structural transitions induced by the external field. (C) 2015 AIP Publishing LLC.

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: The effect of palladium- and ruthenium-based clusters on nanocrystalline tin dioxide interaction with oxygen was studied by temperature-programmed oxygen isotopic exchange with mass-spectrometry detection. The modification of aqueous sol-gel prepared SnO2 by palladium and, to a larger extent, by ruthenium, increases surface oxygen concentration on the materials. The revealed effects on oxygen exchange-lowering the threshold temperature, separation of surface oxygen contribution to the process, increase of heteroexchange rate and oxygen diffusion coefficient, decrease of activation energies of exchange and diffusion-were more intensive for Ru-modified SnO2 than in the case of SnO2/Pd. The superior promoting activity of ruthenium on tin dioxide interaction with oxygen was interpreted by favoring the dissociative O-2 adsorption and increasing the oxygen mobility, taking into account the structure and chemical composition of the modifier clusters.

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: In this paper, the splitting of CO2 in a pulsed plasma system, such as a dielectric barrier discharge (DBD), is evaluated from a chemical point of view by means of numerical modeling. For this purpose, a chemical reaction set of CO2 in an atmospheric pressure plasma is developed, including the vibrational states of CO2, O2, and CO. The simulated pulses are matched to the conditions of a filament (or microdischarge) and repeated with intervals of 1 μs. The influence of vibrationally excited CO2 as well as other neutral species, ions, and electrons on the CO2 splitting is discussed. Our calculations predict that the electrons have the largest contribution to the CO2 splitting at the conditions under study, by electron impact dissociation. The contribution of vibrationally excited CO2 levels in the splitting of CO2 is found be 6.4%, when only considering one microdischarge pulse and its afterglow, but it can be much higher for consecutive discharge pulses, as is typical for a filamentary DBD, when the interpulse time is short enough and accumulation effects in the vibrationally excited CO2 densities can occur.

Abstract: van der Waals (vdW) forces play an important role in the adsorption of molecules on the surface of solids. However, the choice of the most suitable vdW functional for different systems is an essential problem which must be addressed for different systems. The lack of a systematic study on the performance of the vdW functionals in the bulk and adsorption properties of metal-oxides motivated us to examine different vdW approaches and compute the bulk and molecular adsorption properties of alpha-Cr2O3, alpha-Fe2O3, and alpha-Al2O3. For the bulk properties, we compared our results for the heat of formation, cohesive energy, lattice parameters and bond distances between the different vdW functionals and available experimental data. Next we studied the adsorption of benzene and CH3 molecules on top of different oxide surfaces. We employed different approximations to exchange and correlation within DFT, namely, the Perdew-Burke-Ernzerhof (PBE) GGA, (PBE)+U, and vdW density functionals [ DFT(vdW-DF/DF2/optPBE/optB86b/optB88)+U] as well as DFT-D2/D3(+U) methods of Grimme for the bulk calculations and optB86b-vdW(+U) and DFT-D2(+U) for the adsorption energy calculations. Our results highlight the importance of vdW interactions not only in the adsorption of molecules, but importantly also for the bulk properties. Although the vdW contribution in the adsorption of CH3 (as a chemisorption interaction) is less important compared to the adsorption of benzene (as a physisorption interaction), this contribution is not negligible. Also adsorption of benzene on ferryl/chromyl terminated surfaces shows an important chemisorption contribution in which the vdW interactions become less significant.