Abstract: We investigate the role of transition metal atoms of group V-b (V, Nb, Ta) and VI-b (Cr, Mo, W) as n- or p-type dopants in anatase TiO$_2$ using thermodynamic
principles and density functional theory with the Heyd-Scuseria-Ernzerhof HSE06 hybrid functional. The HSE06 functional provides a realistic value for the band gap, which ensures a correct classification of dopants as shallow or deep donors or acceptors. Defect formation energies and thermodynamic transition levels are calculated taking into account the constraints imposed by the stability of TiO$_2$ and the solubility limit of the impurities.
Nb, Ta, W and Mo are identified as shallow donors. Although W provides two electrons, Nb and Ta show a considerable lower formation energy, in particular under O-poor conditions. Mo donates in principle one electron, but under specific conditions can turn into a double donor. V impurities are deep donors and Cr
shows up as an amphoteric defect, thereby acting as an electron trapping center in n-type TiO$_2$ especially under O-rich conditions. A comparison with the available experimental data yields excellent agreement.

Abstract: Strain engineering is an effective way to tune the intrinsic properties of a material. Here, we show by using first-principles calculations that both uniaxial and biaxial tensile strain applied to monolayer TiS3 are able to significantly modify its intrinsic mobility. From the elastic modulus and the phonon dispersion relation we determine the tensile strain range where structure dynamical stability of the monolayer is guaranteed. Within this region, we find more than one order of enhancement of the acoustic phonon limited mobility at 300 K (100 K), i.e. from 1.71 x 10(4) (5.13 x 10(4)) cm(2) V-1 s(-1) to 5.53 x 10(6) (1.66 x 10(6)) cm(2) V-1 s(-1). The degree of anisotropy in both mobility and effective mass can be tuned by using tensile strain. Furthermore, we can either increase or decrease the band gap of TiS3 monolayer by applying strain along different crystal directions. This property allows us to use TiS3 not only in electronic but also in optical applications.

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: 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: <script type='text/javascript'>document.write(unpmarked('To realize spintronic devices based on topological insulators (TIs), well-defined interfaces between magnetic metals and TIs are required. Here, we characterize atomically precisely the interface between the 3d transition metal Fe and the TI Bi2Te3 at different stages of its formation. Using photoelectron diffraction and holography, we show that after deposition of up to 3 monolayers Fe on Bi2Te3 at room temperature, the Fe atoms are ordered at the interface despite the surface disorder revealed by our scanning-tunneling microscopy images. We find that Fe occupies two different sites: a hollow adatom deeply relaxed into the Bi2Te3 quintuple layers and an interstitial atom between the third (Te) and fourth (Bi) atomic layers. For both sites, our core-level photoemission spectra and density-functional theory calculations demonstrate simultaneous chemical bonding of Fe to both Te and Bi atoms. We further show that upon deposition of Fe up to a thickness of 20 nm, the Fe atoms penetrate deeper into the bulk forming a 2-5 nm interface layer containing FeTe. In addition, excessive Bi is pushed down into the bulk of Bi2Te3 leading to the formation of septuple layers of Bi3Te4 within a distance of similar to 25 nm from the interface. Controlling the magnetic properties of the complex interface structures revealed by our work will be of critical importance when optimizing the efficiency of spin injection in TI-based devices.'));

Abstract: A number of recent computational material design studies based on density functional theory (DFT) calculations have put forward a new class of materials with electrically switchable chemical characteristics that can be exploited in the development of tunable gas storage and electrocatalytic applications. We find systematic flaws in almost every computational study of gas adsorption on polarized or charged surfaces, stemming from an improper and unreproducible treatment of periodicity, leading to very large errors of up to 3 eV in some cases. Two simple corrective procedures that lead to consistent results are proposed, constituting a crucial course correction to the research in the field.

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: 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: When confined between graphene layers, water behaves differently from the bulk and exhibits unusual properties such as fast water flow and ordering into a crystal. The hydrogen-bonded network is affected by the limited space and by the characteristics of the confining walls. The presence of an extraordinary number of hydronium and hydroxide ions in narrow channels has the following effects: (i) they affect water permeation through the channel, (ii) they may interact with functional groups on the graphene oxide surface and on the edges, and (iii) they change the thermochemistry of water, which are fundamentally important to understand, especially when confined water is subjected to an external electric field. Here we study the physical properties of water when confined between two graphene sheets and containing hydronium and hydroxide. We found that: (i) there is a disruption in the solvation structure of the ions, which is also affected by the layered structure of confined water, (ii) hydronium and hydroxide occupy specific regions inside the nanochannel, with a prevalence of hydronium (hydroxide) ions at the edges (interior), and (iii) ions recombine more slowly in confined systems than in bulk water, with the recombination process depending on the channel height and commensurability between the size of the molecules and the nanochannel height – a decay of 20% (40%) in the number of ions in 8 ps is observed for a channel height of h = 7 angstrom (bulk water). Our work reveals distinctive properties of water confined in a nanocapillary in the presence of additional hydronium and hydroxide ions.

Abstract: Defects are inevitably present in materials, and their existence in a material strongly affects its fundamental physical properties. We have systematically investigated the effects of surface adsorption, substitutional impurities, defect engineering, an electric field and strain engineering on the structural, electronic and magnetic properties of antimonene nanosheets, using spin-polarized density functional calculations based on first-principles. The adsorption or substitution of atoms can locally modify the atomic and electronic structures as well as induce a variety of electronic behaviors including metal, half-metal, ferromagnetic metal, dilute magnetic semiconductor and spin-glass semiconductor. Our calculations show that the presence of typical defects (vacancies and Stone-Wales defect) in antimonene affects the geometrical symmetry as well as the band gap in the electronic band structure and induces magnetism to antimonene. Moreover, by applying an external electric field and strain (uniaxial and biaxial), the electronic structure of antimonene can be easily modified. The calculation results presented in this paper provide a fundamental insight into the tunable nature of the electronic properties of antimonene, supporting its promise for use in future applications.

Abstract: The mechanisms of plasma in medicine are broadly attributed to plasma-derived reactive oxygen and nitrogen species (RONS). In order to exert any intracellular effects, these plasma-derived RONS must first traverse a major barrier in the cell membrane. The cell membrane lipid composition, and thereby the magnitude of this barrier, is highly variable between cells depending on type and state (e.g. it is widely accepted that healthy and cancerous cells have different membrane lipid compositions). In this study, we investigate how plasma-derived RONS interactions with lipid membrane components can potentially be exploited in the future for treatment of diseases. We couple phospholipid vesicle experiments, used as simple cell models, with molecular dynamics (MD) simulations of the lipid membrane to provide new insights into how the interplay between phospholipids and cholesterol may influence the response of healthy and diseased cell membranes to plasma-derived RONS. We focus on the (i) lipid tail saturation degree, (ii) lipid head group type, and (iii) membrane cholesterol fraction. Using encapsulated molecular probes, we study the influence of the above membrane components on the ingress of RONS into the vesicles, and subsequent DNA damage. Our results indicate that all of the above membrane components can enhance or suppress RONS uptake, depending on their relative concentration within the membrane. Further, we show that higher RONS uptake into the vesicles does not always correlate with increased DNA damage, which is attributed to ROS reactivity and lifetime. The MD simulations indicate the multifactorial chemical and physical processes at play, including (i) lipid oxidation, (ii) lipid packing, and (iii) lipid rafts formation. The methods and findings presented here provide a platform of knowledge that could be leveraged in the development of therapies relying on the action of plasma, in which the cell membrane and oxidative stress response in cells is targeted.

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.

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.