Abstract: The interaction of silicene, the silicon counterpart of graphene, with (0001) ZnS surfaces is investigated theoretically, using first-principles simulations. The charge transfer occurring at the silicene/(0001) ZnS interface leads to the opening of an indirect energy band gap of about 0.7 eV in silicene. Remarkably, the nature (indirect or direct) and magnitude of the energy band gap of silicene can be controlled by an external electric field: the energy gap is predicted to become direct for electric fields larger than about 0.5 V angstrom(-1), and the direct energy gap decreases approximately linearly with the applied electric field. The predicted electric field tunable energy band gap of the silicene/(0001) ZnS interface is very promising for its potential use in nanoelectronic devices.

Abstract: Monte Carlo simulations of coronene molecules in single-walled carbon nanotubes (SWCNTs) and dicoronylene molecules in SWCNTs are performed. Depending on the diameter D of the encapsulating SWCNT, regimes favoring the formation of ordered, one-dimensional (1D) stacks of tilted molecules (D <= 1.7 nm for coronene@SWCNT, 1.5 nm <= D <= 1.7 nm for dicoronylene@SWCNT) and regimes with disordered molecular arrangements and increased translational mobilities enabling the thermally induced polymerization of neighboring molecules resulting in the formation of graphene nanoribbons (GNRs) are observed. The results show that the diameter of the encapsulating nanotube is a crucial parameter for the controlled synthesis of either highly ordered 1D structures or GNR precursors.

Abstract: Using fixed-node diffusion quantum Monte Carlo (FN-DMC) simulation we investigate the electron correlation in all-metal aromatic clusters MAl4- (with M = Li, Na, K, Rb, Cu, Ag and Au). The electron detachment energies and electron affinities of the clusters are obtained. The vertical electron detachment energies obtained from the FN-DMC calculations are in very good agreement with the available experimental results. Calculations are also performed within the Hartree-Fock approximation, density-functional theory (DFT), and the couple-cluster (CCSD(T)) method. From the obtained results, we analyse the impact of the electron correlation effects in these bimetallic clusters and find that the correlation of the valence electrons contributes significantly to the detachment energies and electron affinities, varying between 20% and 50% of their total values. Furthermore, we discuss the electron correlation effects on the stability of the clusters as well as the accuracy of the DFT and CCSD(T) calculations in the present systems.

Abstract: The absence of a direct-to-indirect band gap transition in ReS2 when going from the monolayer to bulk makes it special among the other semiconducting transition metal dichalcogenides. The functionalization of this promising layered material emerges as a necessity for the next generation technological applications. Here, the structural, electronic, and magnetic properties of substitutionally doped ReS2 monolayers at either the S or Re site were systematically studied by using first principles density functional calculations. We found that substitutional doping of ReS2 depends sensitively on the growth conditions of ReS2. Among the large number of non-metallic atoms, namely H, B, C, Se, Te, F, Br, Cl, As, P. and N, we identified the most promising candidates for n-type and p-type doping of ReS2. While Cl is an ideal candidate for n-type doping, P appears to be the most promising candidate for p-type doping of the ReS2 monolayer. We also investigated the doping of ReS2 with metal atoms, namely Mo, W, Ti, V. Cr, Co, Fe, Mn, Ni, Cu, Nb, Zn, Ru, Os and Pt. Mo, Nb, Ti, and V atoms are found to be easily incorporated in a single layer of ReS2 as substitutional impurities at the Re site for all growth conditions considered in this work. Tuning chemical potentials of dopant atoms energetically makes it possible to dope ReS2 with Fe, Co, Cr, Mn, W, Ru, and Os at the Re site. We observe a robust trend for the magnetic moments when substituting a Re atom with metal atoms such that depending on the electronic configuration of dopant atoms, the net magnetic moment of the doped ReS2 becomes either 0 or 1 mu(B). Among the metallic dopants, Mo is the best candidate for p-type doping of ReS2 owing to its favorable energetics and promising electronic properties.

Abstract: Past research efforts aiming at obtaining stable p-type ZnO have been based on complexes involving nitrogen doping. A recent experiment by (J. G. Reynolds et al., Appl. Phys. Lett., 2013, 102, 152114) demonstrated a significant ([similar]1018 cm−3) p-type behavior in N-doped ZnO films after appropriate annealing. The p-type conductivity was attributed to a VZnNOH shallow acceptor complex, formed by a Zn vacancy (VZn), N substituting O (NO), and H interstitial (Hi). We present here a first-principles hybrid functional study of this complex compared to the one without hydrogen. Our results confirm that the VZnNOH complex acts as an acceptor in ZnO. We find that H plays an important role, because it lowers the formation energy of the complex with respect to VZnNO, a complex known to exhibit (unstable) p-type behavior. However, this additional H atom also occupies the hole level at the origin of the shallow behavior of VZnNO, leaving only two states empty higher in the band gap and making the VZnNOH complex a deep acceptor. Therefore, we conclude that the cause of the observed p-type conductivity in experiment is not the presence of the VZnNOH complex, but probably the formation of the VZnNO complex during the annealing process.

Abstract: Ionic diffusion through the novel (AlxMg1-2xLix)Al2O4 spinel electrolyte is investigated using first-principles calculations, combined with the Kinetic Monte Carlo algorithm. We observe that the ionic diffusion increases with the lithium content x. Furthermore, the structural parameters, formation enthalpies and electronic structures of (AlxMg1-2xLix)Al2O4 are calculated for various stoichiometries. The overall results indicate the (AlxMg1-2xLix)Al2O4 stoichiometries x = 0.2...0.3 as most promising. The (AlxMg1-2xLix)Al2O4 electrolyte is a potential candidate for the all-spinel solid-state battery stack, with the material epitaxially grown between well-known spinel electrodes, such as LiyMn2O4 and Li4+3yTi5O12 (y = 0...1). Due to their identical crystal structure, a good electrolyte-electrode interface is expected.

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: Gold nanocluster properties exhibit unique size-dependence. In this contribution, we employ reactive molecular dynamics simulations to calculate the size- and temperature-dependent surface energies, strain energies and atomic displacements for icosahedral, cuboctahedral, truncated octahedral and decahedral Au-nanoclusters. The calculations demonstrate that the surface energy decreases with increasing cluster size at 0 K but increases with size at higher temperatures. The calculated melting curves as a function of cluster size demonstrate the Gibbs-Thomson effect. Atomic displacements and strain are found to strongly depend on the cluster size and both are found to increase with increasing cluster size. These results are of importance for understanding the size-and temperature-dependent surface processes on gold nanoclusters.

Abstract: Hydrogen plays an essential role during the in situ assembly of tailored catalytic materials, and serves as key ingredient in multifarious chemical reactions promoted by these catalysts. Despite intensive debate for several decades, the existence and nature of hydrogen-involved mechanisms – such as hydrogen-spillover, surface migration – have not been unambiguously proven and elucidated up to date. Here, Pt-Ga alloy formation is used as a probe reaction to study the behavior and atomic transport of H and Ga, starting from Pt nanoparticles on hydrotalcite-derived Mg(Ga)(Al)Ox supports. In situ XANES spectroscopy, time-resolved TAP kinetic experiments, HAADF-STEM imaging and EDX mapping are combined to probe Pt, Ga and H in a series of H2 reduction experiments up to 650 degrees C. Mg(Ga)(Al)Ox by itself dissociates hydrogen, but these dissociated hydrogen species do not induce significant reduction of Ga3+ cations in the support. Only in the presence of Pt, partial reduction of Ga3+ into Gadelta+ is observed, suggesting that different reaction mechanisms dominate for Pt- and Mg(Ga)(Al)Ox-dissociated hydrogen species. This partial reduction of Ga3+ is made possible by Pt-dissociated H species which spillover onto non-reducible Mg(Al)Ox or partially reducible Mg(Ga)(Al)Ox and undergo long-range transport over the support surface. Moderately mobile Gadelta+Ox migrates towards Pt clusters, where Gadelta+ is only fully reduced to Ga0 on condition of immediate stabilization inside Pt-Ga alloyed nanoparticles.

Abstract: Recently penta-graphene was proposed as a stable two-dimensional carbon allotrope consisting of a single layer of interconnected carbon pentagons [Zhang et al., PNAS, 2015, 112, 2372]. Its silicon counterpart, penta-silicene, however, is not stable. In this work, we show that multilayers of penta-silicene form stable materials with semiconducting or metallic properties, depending on the stacking mode. We demonstrate their dynamic stability through their phonon spectrum and using molecular dynamics. A particular type of bilayer penta-silicene is found to have lower energy than all of the known hexagonal silicene bilayers and forms therefore the most stable bilayer silicon material predicted so far. The electronic and mechanical properties of these new silicon allotropes are studied in detail and their behavior under strain is investigated. We demonstrate that strain can be used to tune its band gap.

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.

Abstract: Despite the numerous studies on the adsorption of different proteins onto mesoporous titanium dioxide and indications on the important role of buffer solutions in bioactivity, a systematic study on the impact of the buffer on the protein incorporation into porous substrates is still lacking. We here studied the interaction between a commercial mesoporous TiO2 and three of the most used buffers for protein incorporation, i.e. HEPES, Tris and phosphate buffer. In addition, this paper analyzes the adsorption of horse heart myoglobin (hhMb) onto commercial mesoporous TiO2 as a model system to test the influence of buffers on the protein incorporation behavior in mesoporous TiO2. N2 sorption analysis, FT-IR and TGA/DTG measurements were used to evaluate the interaction between the buffers and the TiO2 surface, and the effect of such an interaction on hhMb adsorption. Cyclic voltammetry (CV) and electron paramagnetic resonance (EPR) were used to detect changes in the microenvironment surrounding the heme. The three buffers show a completely different interaction with the TiO2 surface, which drastically affects the adsorption of myoglobin as well as its structure and electrochemical activity. Therefore, special attention is required while choosing the buffer medium to avoid misguided evaluation of protein adsorption on mesoporous TiO2.

Abstract: In this study, we investigated the diffusion of H-atoms to the subsurface and their further diffusion into the bulk of a Ni(111) crystal by means of density functional theory calculations in the context of thermal and plasma-assisted catalysis. The H-atoms at the surface can originate from the dissociative adsorption of H2 or CH4 molecules, determining the surface H-coverage. When a threshold H-coverage is passed, corresponding to 1.00 ML for the crystalline Ni(111) surface, the surface-bound H-atoms start to diffuse to the subsurface. A similar threshold coverage is observed for the interstitial H-coverage. Once the interstitial sites are filled up with a coverage above 1.00 ML of H, dissolution of interstitial H-atoms to the layer below the interstitial sites will be initiated. Hence, by applying a high pressure or inducing a reactive plasma and high temperature, increasing the H-flux to the surface, a large amount of hydrogen can diffuse in a crystalline metal like Ni and can be absorbed. The formation of metal hydride may modify the entire reaction kinetics of the system. Equivalently, the H-atoms in the bulk can easily go back to the surface and release a large amount of heat. In a plasma process, H-atoms are formed in the plasma, and therefore the energy barrier for dissociative adsorption is dismissed, thus allowing achievement of the threshold coverage without applying a high pressure as in a thermal process. As a result, depending on the crystal plane and type of metal, a large number of H-atoms can be dissolved (absorbed) in the metal catalyst, explaining the high efficiency of plasma-assisted catalytic reactions. Here, the mechanism of H-dissolution is established as a chemical identifier for the investigation of the reaction kinetics of a chemical process.

Abstract: A combined experimental and first-principles study is performed to study the origin of conductivity in
ZnO:Al nanoparticles synthesized under controlled conditions via a reflux route using benzylamine as a
solvent. The experimental characterization of the samples by Raman, nuclear magnetic resonance (NMR)
and conductivity measurements indicates that upon annealing in nitrogen, the Al atoms at interstitial
positions migrate to the substitutional positions, creating at the same time Zn interstitials. We provide
evidence for the fact that the formed complex of AlZn and Zni corresponds to the origin of the Knight
shifted peak (KS) we observe in 27Al NMR. As far as we know, the role of this complex has not been
discussed in the literature to date. However, our first-principles calculations show that such a complex is
indeed energetically favoured over the isolated Al interstitial positions. In our calculations we also
address the charge state of the Al interstitials. Further, Zn interstitials can migrate from Al_Zn and possibly
also form Zn clusters, leading to the observed increased conductivity.

Abstract: The vast biomedical potential of cold atmospheric pressure plasmas (CAPs) is governed by the formation of reactive species. These biologically active species are formed upon the interaction of CAPs with the surroundings. In biological milieu, water plays an essential role. The development of biomedical CAPs thus requires understanding of the sources of the reactive species in aqueous media exposed to the plasma. This is especially important in case of the COST RF plasma jet, which is developed as a reference microplasma system. In this work, we investigated the formation of the OH radicals, H atoms and H2O2 in aqueous solutions exposed to the COST plasma jet. This was done by combining experimental and modelling approaches. The liquid phase species were analysed using UV-Vis spectroscopy and spin trapping with hydrogen isotopes and electron paramagnetic resonance (EPR) spectroscopy. The discrimination between the species formed from the liquid phase and the gas phase molecules was performed by EPR and 1H-NMR analyses of the liquid samples. The concentrations of the reactive species in the gas phase plasma were obtained using a zero-dimensional (0D) chemical kinetics computational model. A three-dimensional (3D) fluid dynamics model was developed to provide information on the induced humidity in the plasma effluent. The comparison of the experimentally obtained trends for the formation of the species as a function of the feed gas and effluent humidity with the modelling results suggest that all reactive species detected in our system are mostly formed in the gas phase plasma inside the COST jet, with minor amounts arising from the plasma effluent humidity.