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: ZnM2O4 (M = Co, Rh, Ir) spinels are considered as a class of potential p-type transparent conducting oxides (TCOs). We report the formation energy of acceptor-like defects using first principles calculations with an advanced hybrid exchange-correlation functional (HSE06) within density functional theory (DFT). Due to the discrepancies between the theoretically obtained band gaps with this hybrid functional and the – scattered – experimental results, we also perform GW calculations to support the validity of the description of these spinels with the HSE06 functional. The considered defects are the cation vacancy and antisite defects, which are supposed to be the leading source of disorder in the spinel structures. We also discuss the band alignments in these spinels. The calculated formation energies indicate that the antisite defects ZnM (Zn replacing M, M = Co, Rh, Ir) and VZn act as shallow acceptors in ZnCo2O4, ZnRh2O4 and ZnIr2O4, which explains the experimentally observed p-type conductivity in those systems. Moreover, our systematic study indicates that the ZnIr antisite defect has the lowest formation energy in the group and it corroborates the highest p-type conductivity reported for ZnIr2O4 among the group of ZnM2O4 spinels. To gain further insight into factors affecting the p-type conductivity, we have also investigated the formation of localized small polarons by calculating the self-trapping energy of the holes.

Abstract: Nitrogen-doped graphene was recently synthesized and was reported to be a catalyst for hydrogen dissociative adsorption under a perpendicular applied electric field (F). In this work, the diffusion of H atoms on N-doped graphene, in the presence and absence of an applied perpendicular electric field, is studied using density functional theory. We demonstrate that the applied field can significantly facilitate the binding of hydrogen molecules on N-doped graphene through dissociative adsorption and diffusion on the surface. By removing the applied field the absorbed H atoms can be released efficiently. Our theoretical calculation indicates that N-doped graphene is a promising hydrogen storage material with reversible hydrogen adsorption/desorption where the applied electric field can act as a switch for the uptake/release processes.

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: Motivated by the recent successful synthesis of highly crystalline ultrathin BiTeCl and BiTeBr layered sheets [Debarati Hajra et al., ACS Nano, 2020, 14, 15626], herein for the first time, we carry out a comprehensive study on the structural and electronic properties of BiTeCl and BiTeBr Janus monolayers using density functional theory (DFT) calculations. Different structural and electronic parameters including the lattice constant, bond lengths, layer thickness in the z-direction, different interatomic angles, work function, charge density difference, cohesive energy and Rashba coefficients are determined to acquire a deep understanding of these monolayers. The calculations show good stability of the studied single layers. BiTeCl and BiTeBr monolayers are semiconductors with electronic bandgaps of 0.83 and 0.80 eV, respectively. The results also show that the semiconductor-metal transformation can be induced by increasing the number of layers. In addition, the engineering of the electronic structure is also studied by applying an electric field, and mechanical uniaxial and biaxial strain. The results show a significant change of the bandgaps and that an indirect-direct band-gap transition can be induced. This study highlights the positive prospect for the application of BiTeCl and BiTeBr layered sheets in novel electronic and energy conversion systems.

Abstract: In the present work, the effect of various embedded atom impurities on tuning electronic and magnetic properties of C3N4 and C4N3 nanosheets have been studied using first-principles calculations. Our calculations show that C3N4 is a semiconductor and it exhibits extraordinary electronic properties such as dilute-magnetic semiconductor (with H, F, Cl, Be, V, Fe and Co); metal (with N, P, Mg and Ca), half-metal (with Li, Na, K, Al, Sc, Cr, Mn, and Cu) and semiconductor (with O, S, B, C, Si, Ti, Ni and Zn) with the band gaps in the range of 0.3-2.0 eV depending on the species of embedded atom. The calculated electronic properties reveal that C4N3 is a half-metal and it retains half-metallic character with embedded H, O, S, F, B, N, P, Be, Mg, Al, Sc, V, Fe, Ni and Zn atoms. The substitution of Cl, C, Cr and Mn atoms create ferromagnetic-metal character in the C4N3 nanosheet, embedded Co and Cu atoms exhibit a dilute-magnetic semiconductor nature, and embedded Ti atoms result in the system becoming a semiconductor. Therefore, our results reveal the fact that the band gap and magnetism can be modified or induced by various atom impurities, thus, offering effective possibilities to tune the electronic and magnetic properties of C3N4 and C4N3 nanosheets.

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: 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: 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: 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: 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 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: Recently, Sr-doped LaCrO3 has been experimentally introduced as a new p-type transparent conducting oxide. It is demonstrated that substituting Sr for La results in inducing p-type conductivity in LaCrO3. Performing first principles calculations we study the electronic structure and formation energy of various point defects in LaCrO3. Our results for the formation energies show that in addition to Sr, two more divalent defects, Ca and Ba, substituting for La in LaCrO3, behave as shallow acceptors in line with previous experimental reports. We further demonstrate that under oxygen-poor growth conditions, these shallow acceptors will be compensated by intrinsic donor-like defects (an oxygen vacancy and Cr on an oxygen site), but in the oxygen-rich growth regime the shallow acceptors have the lowest formation energies between all considered defects and will lead to p-type conductivity.

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.