Abstract: Two-dimensional Transition Metal Dichalcogenides (TMDs) possessing extraordinary physical properties at reduced dimensionality have attracted interest due to their promise in electronic and optical device applications. However, TMD monolayers can show a broad range of different properties depending on their crystal phase; for example, H phases are usually semiconductors, while the T phases are metallic. Thus, controlling phase transitions has become critical for device applications. In this study, the energetically low-lying crystal structures of pristine and Janus TMDs are investigated by using ab initio Nudged Elastic Band and molecular dynamics simulations to provide a general explanation for their phase stability and transition properties. Across all materials investigated, the T phase is found to be the least stable and the H phase is the most stable except for WTe2, while the T' and T '' phases change places according to the TMD material. The transition energy barriers are found to be large enough to hint that even the higher energy phases are unlikely to undergo a phase transition to a more stable phase if they can be achieved except for the least stable T phase, which has zero barrier towards the T ' phase. Indeed, in molecular dynamics simulations the thermodynamically least stable T phase transformed into the T ' phase spontaneously while in general no other phase transition was observed up to 2100 K for the other three phases. Thus, the examined T ', T '' and H phases were shown to be mostly stable and do not readily transform into another phase. Furthermore, so-called mixed phase calculations considered in our study explain the experimentally observed lateral hybrid structures and point out that the coexistence of different phases is strongly stable against phase transitions. Indeed, stable complex structures such as metal-semiconductor-metal architectures, which have immense potential to be used in future device applications, are also possible based on our investigation.

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: 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: We have characterized theoretically the work-function modifications of the (111) surfaces of gold and silver upon deposition of self-assembled monolayers based on methanethiol and trifluoromethanethiol. A comparative analysis is made between the experimental results and those obtained from two widely used approaches based on density functional theory. The contributions to the total work-function modifications are estimated on the basis of two decomposition schemes of the thiol molecules that have been proposed in the literature. The contributions are found to differ significantly between the two approaches, as do the corresponding adsorption energies.

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: 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: Two-dimensional (2D) materials have attracted growing interest due to their potential use in the next generation of nanoelectronic and optoelectronic applications. On the basis of first-principles calculations based on density functional theory, we first investigate the electronic and mechanical properties of single layer boron phosphide (h-BP). Our calculations show that h-BP is a mechanically stable 2D material with a direct band gap of 0.9 eV at the K-point, promising for both electronic and optoelectronic applications. We next investigate the electron transport properties of a p-n junction constructed from single layer boron phosphide (h-BP) using the non-equilibrium Green's function formalism. The n-and p-type doping of BP are achieved by substitutional doping of B with C and P with Si, respectively. C(Si) substitutional doping creates donor (acceptor) states close to the conduction (valence) band edge of BP, which are essential to construct an efficient p-n junction. By modifying the structure and doping concentration, it is possible to tune the electronic and transport properties of the p-n junction which exhibits not only diode characteristics with a large current rectification but also negative differential resistance (NDR). The degree of NDR can be easily tuned via device engineering.

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: Copper deposited by thermal evaporation onto pristine and oxygen plasma treated carbon nanotubes (CNTs) diffuse over the CNT surface, coalescing and forming crystalline islands. The nucleation sites of the islands are preferentially defects, and more homogeneous island dispersion was observed at the CNT oxygen functionalized surface. The presence of weakly bound oxygen atoms at the CNT surface induces the formation of CuO bonds at the Cu/CNT interface, as described through density functional calculations. Exposure to air allows further oxidation to facetted crystalline Cu2O. Oxygen plasma pre-treatment represents a promising route for homogenous disperse Cu2O nanoparticle decoration of CNTs.

Abstract: The interaction between evaporated Pt and pristine or oxygen-plasma-treated multiwall carbon nanotubes (CNTs) is investigated. Pt is found to nucleate at defect sites, whether initially present or introduced by oxygen plasma treatment. The plasma treatment induces a uniform dispersion of Pt nanoparticles at the CNT surface. The absence of additional features in the C 1s core level spectrum indicates that no mixed PtC phase is formed. The formation of COPt bonds at the cluster-CNT interface is suggested to reduce the electronic interaction between Pt nanoparticles and the CNT surface.

Abstract: Chalcopyrite semiconductors are of considerable interest for application as absorber layers in thin-film photovoltaic cells. When growing films of these compounds, however, they are often found to contain CuAu-like domains, a metastable phase of chalcopyrite. It has been reported that for CuInS2, the presence of the CuAu-like phase improves the short circuit current of the chalcopyrite-based photovoltaic cell. We investigate the thermodynamic stability of both phases for a selected list of I-III-VI2 materials using a first-principles density functional theory approach. For the CuIn-VI2 compounds, the difference in formation energy between the chalcopyrite and CuAu-like phase is found to be close to 2 meV per atom, indicating a high likelihood of the presence of CuAu-like domains. Next, we calculate the spectroscopic limited maximum efficiency (SLME) of the CuAu-like phase and compare the results with those of the corresponding chalcopyrite phase. We identify several candidates with a high efficiency, such as CuAu-like CuInS2, for which we obtain an SLME of 29% at a thickness of 500 nm. We observe that the SLME can have values above the Shockley-Queisser (SQ) limit, and show that this can occur because the SQ limit assumes the absorptivity to be a step function, thus overestimating the radiative recombination in the detailed balance approach. This means that it is possible to find higher theoretical efficiencies within this framework simply by calculating the J-V characteristic with an absorption spectrum. Finally, we expand our SLME analysis to indirect band gap absorbers by studying silicon, and find that the SLME quickly overestimates the reverse saturation current of indirect band gap materials, drastically lowering their calculated efficiency.

Abstract: We have performed a first-principles study of the p- and n-type conductivity in CuIn1−xGaxSe2 due to native point defects, based on the HSE06 hybrid functional. Band alignment shows that the band gap becomes larger with x due to the increasing conduction band minimum, rendering it hard to establish n-type conductivity in CuGaSe2. From the defect formation energies, we find that In/GaCu is a shallow donor, while VCu, VIn/Ga and CuIn/Ga act as shallow acceptors. Using the total charge neutrality of ionized defects and intrinsic charge carriers to determine the Fermi level, we show that under In-rich growth conditions InCu causes strongly n-type conductivity in CuInSe2. Under increasingly In-poor growth conditions, the conductivity type in CuInSe2 alters to p-type and compensation of the acceptors by InCu reduces, as also observed in photoluminescence experiments. In CuGaSe2, the native acceptors pin the Fermi level far away from the conduction band minimum, thus inhibiting n-type conductivity. On the other hand, CuGaSe2 shows strong p-type conductivity under a wide range of Ga-poor growth conditions. Maximal p-type conductivity in CuIn1−xGaxSe2 is reached under In/Ga-poor growth conditions, in agreement with charge concentration measurements on samples with In/Ga-poor stoichiometry, and is primarily due to the dominant acceptor CuIn/Ga.

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: 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.