Abstract: We report on a first-principles study of the structural and electronic properties of a Sigma3 (112) grain boundary model in CuInSe2. The study focuses on a coherent, stoichiometry preserving, cation–Se terminated grain boundary, addressing the properties of the grain boundary as such, as well as the effect
of well known defects in CuInSe2. We show that in spite of its apparent simplicity, such a grain boundary exhibits a very rich phenomenology, providing an explanation for several of the experimentally observed properties of grain boundaries in CuInSe2 thin films. In particular, we show that the combined effect of Cu vacancies and cation antisites can result in the observed Cu depletion with no In enrichment at the grain boundaries. Furthermore, Cu vacancies are unlikely to produce a hole barrier at the grain boundaries, but Na may indeed have such an effect. We find that Na-on-Cu defects will tend to form abundantly at
the grain boundaries, and can provide a mechanism for the carrier depletion and/or type inversion experimentally reported.

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: 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: 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: 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: Using first principles calculations we investigate the structural and electronic properties of interfaces between fluorographane and MoS2. Unsymmetrical functionalization of graphene with H and F results in an intrinsic dipole moment perpendicular to the plane of the buckled graphene skeleton. Depending on the orientation of this dipole moment, the electronic properties of a physically absorbed MoS2 monolayer can be switched from n-to p-type or vice versa. We show that one can realize vanishing n-type/p-type Schottky barrier heights when contacting MoS2 to fluorographane. By applying a perpendicular electric field, the size of the Schottky barrier and the degree of doping can be tuned. Our calculations indicate that a fluorographane monolayer is a promising candidate for bipolar doping of MoS2, which is vital in the design of novel technological applications based on two-dimensional materials.

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: Perovskite oxides are active room-temperature bifunctional oxygen electrocatalysts in alkaline media, capable of performing the oxygen reduction reaction (ORR) and oxygen evolution reaction (OER) with lower combined overpotentials relative to their precious metal counterparts. However, their semiconducting nature necessitates the use of activated carbons as conductive supports to generate applicably relevant current densities. In efforts to advance the performance and theory of oxide electrocatalysts, the chemical and physical properties of the oxide material often take precedence over contributions from the conductive additive. In this work, we find that carbon plays an important synergistic role in improving the performance of La1-xSrxCoO3- (0 x 1) electrocatalysts through the activation of O-2 and spillover of radical oxygen intermediates, HO2- and O-2(-), which is further reduced through chemical decomposition of HO2- on the perovskite surface. Through a combination of thin-film rotating disk electrochemical characterization of the hydrogen peroxide intermediate reactions (hydrogen peroxide reduction reaction (HPRR), hydrogen peroxide oxidation reaction (HPOR)) and oxygen reduction reaction (ORR), surface chemical analysis, HR-TEM, and microkinetic modeling on La1-xSrxCoO3- (0 x 1)/carbon (with nitrogen and non-nitrogen doped carbons) composite electrocatalysts, we deconvolute the mechanistic aspects and contributions to reactivity of the oxide and carbon support.

Abstract: Cold atmospheric plasma in contact with solutions has many applications, but its chemistry contains many unknowns such as the undescribed reactions with solutes. By combining experiments and modelling, we report the first direct demonstration of the reaction of chloride with oxygen atoms in aqueous solutions exposed to cold plasma.

Abstract: Modeling the oxidation process of silicon nanowires through reactive force field based molecular dynamics simulations suggests that the formation of Si epoxide defects occurs both at the Si/SiOx interface and at the nanowire surface, whereas for flat surfaces, this defect is experimentally observed to occur only at the interface as a result of stress. In this paper, we argue that the increasing curvature stabilizes the defect at the nanowire surface, as suggested by our density functional theory calculations. The latter can have important consequences for the opto-electronic properties of thin silicon nanowires, since the epoxide induces an electronic state within the band gap. Removing the epoxide defect by hydrogenation is expected to be possible but becomes increasingly difficult with a reduction of the diameter of the nanowires.

Abstract: The thermodynamic aspects of various 2D materials are explored using Density Functional Theory (DFT). Various metal chalcogenides (MX2, M = metal, chalcogen X = S, Se, Te) are investigated with respect to their interaction and stability under different ambient conditions met in the integration process of a transistor device. Their interaction with high- dielectrics is also addressed, in order to assess their possible integration in Complementary Metal Oxide Semiconductor (CMOS) field effect transistors. 2D materials show promise for high performance nanoelectronic devices, but the presence of defects (vacancies, grain boundaries,...) can significantly impact their electronic properties. To assess the impact of defects, their enthalpies of formation and their signature levels in the density of states have been studied. We find, consistently with literature reports, that chalcogen vacancies are the most likely source of defects. It is shown that while pristine 2D materials are in general stable whenever set in contact with different ambient atmospheres, the presence of defective sites affects the electronic properties of the 2D materials to varying degrees. We observe that all the 2D materials studied in the present work show strong reactivity towards radical oxygen plasma treatments while reactivity towards other common gas phase chemical such as O-2 and H2O and groups present at the high- surface varies significantly between species. While energy band-gaps, effective masses and contact resistivities are key criteria in selection of 2D materials for scaled CMOS and tunneling based devices, the phase and ambient stabilities might also play a very important role in the development of reliable nanoelectronic applications.

Abstract: In this study, it is predicted by density functional theory calculations that graphene-like novel ultra-thin phases of manganese fluoride crystals, that have nonlayered structures in their bulk form, can be stabilized by fluorination of manganese dichalcogenide crystals. First, it is shown that substitution of fluorine atoms with chalcogens in the manganese dichalcogenide host lattice is favorable. Among possible crystal formations, three stable ultra-thin structures of manganese fluoride, 1H-MnF2, 1T-MnF2 and MnF3, are found to be stable by total energy optimization calculations. In addition, phonon calculations and Raman activity analysis reveal that predicted novel single-layers are dynamically stable crystal structures displaying distinctive characteristic peaks in their vibrational spectrum enabling experimental determination of the corresponding phases. Differing from 1H-MnF2 antiferromagnetic (AFM) large gap semiconductor, 1T-MnF2 and MnF3 single-layers are semiconductors with ferromagnetic (FM) ground state.

Abstract: While spectroscopic data on small hydrocarbons in interstellar media in combination with crossed molecular beam (CMB) experiments have provided a wealth of information on astrochemically relevant species, much of the underlying mechanistic pathways of their formation remain elusive. Therefore, in this work, the chemical reaction mechanisms of C(³P_J) + C₆H₆and C⁺(²P) + C₆H₆systems using the quantum mechanical molecular dynamics (QMMD) technique at the PBE0-D3(BJ) level of theory is investigated, mimicking a CMB experiment. Both the dynamics of the reactions as well as the electronic structure for the purpose of the reaction network are evaluated. The method is validated for the first reaction by comparison to the available experimental data. The reaction scheme for the C(³P_J) + C₆H₆system covers the literature data,<italic>e.g.</italic>the major products are the 1,2-didehydrocycloheptatrienyl radical (C₇H₅) and benzocyclopropenyl radical (C₆H₅–CH), and it reveals the existence of less common pathways for the first time. The chemistry of the C⁺(²P_J) + C₆H₆system is found to be much richer, and we have found that this is because of more exothermic reactions in this system in comparison to those in the C(³P_J) + C₆H₆system. Moreover, using the QMMD simulation, a number of reaction paths have been revealed that produce three distinct classes of reaction products with different ring sizes. All in all, at all the collision energies and orientations, the major product is the heptagon molecular ion for the ionic system. It is also revealed that the collision orientation has a dominant effect on the reaction products in both systems, while the collision energy mostly affects the charged system. These simulations both prove the applicability of this approach to simulate crossed molecular beams, and provide fundamental information on reactions relevant for the interstellar medium.

Abstract: Research progress on single layer group III monochalcogenides has been increasing rapidly owing to their interesting physics. Herein, we investigate the dynamically stable single layer forms of XBi (X = Ge, Si or Sn) using density functional theory calculations. Phonon band dispersion calculations and ab initio molecular dynamics simulations reveal the dynamical and thermal stability of the considered monolayers. Raman spectra calculations indicate the existence of 5 Raman active phonon modes, 3 of which are prominent and can be observed in possible Raman measurements. The electronic band structures of the XBi single layers were investigated with and without the effects of spin-orbit coupling (SOC). Our results show that XBi single layers show semiconducting properties with narrow band gap values without SOC. However, only single layer SiBi is an indirect band gap semiconductor, while GeBi and SnBi exhibit metallic behaviors when adding spin-orbit coupling effects. In addition, the calculated linear elastic parameters indicate the soft nature of the predicted monolayers. Moreover, our predictions for the thermoelectric properties of single layer XBi reveal that SiBi is a good thermoelectric material with increasing temperature. Overall, it is proposed that single layer XBi structures can be alternative, stable 2D single layers with varying electronic and thermoelectric properties.