Abstract: Recent years have seen a tremendous rise of two-dimensional (2D) magnetic materials, several of which were verified experimentally. However, most of the theoretical predictions to date rely on ab initio methods, at zero temperature and fluctuation-free, while one certainly expects detrimental quantum fluctuations at finite temperatures. Here, we present the solution of the quantum Heisenberg model for honeycomb/hexagonal lattices with anisotropic exchange interaction up to third nearest neighbors and in an applied field in arbitrary direction, which answers the question whether long-range magnetization can indeed survive in the ultrathin limit of materials, up to which temperature, and what the characteristic excitation (magnon) frequencies are, all essential to envisaged applications of magnetic 2D materials. We find that long-range magnetic order persists at finite temperature for materials with overall easy-axis anisotropy. We validate the calculations on the examples of monolayers CrI3, CrBr3, and MnSe2. Moreover, we provide an easy-to-use tool to calculate Curie temperatures of new 2D computational materials.

Abstract: We investigate experimentally the terahertz (THz) optoelectronic properties of monolayer (ML) tungsten disulfide (WS2) placed on different substrates using THz time-domain spectroscopy (TDS). We find that the THz optical response of n-type ML WS2 depends sensitively on the choice of the substrate. This dependence is found to be a consequence of substrate induced charge transfer, extra scattering centers, and electronic localization. Through fitting the experimental results with the Drude-Smith formula, we can determine the key sample parameters (e.g., the electronic relaxation time, electron density, and electronic localization factor) of ML WS2 on different substrates. The temperature dependence of these parameters is examined. Our results show that the THz TDS technique is an efficient non-contact method that can be utilized to characterize and investigate the optoelectronic properties of nano-devices based on ML WS2.

Abstract: The two-dimensional (2D) transition metal trihalide (TMX3, X = Cl, Br, I) family has attracted considerable attention in recent years due to the realization of CrCl3, CrBr3, and CrI3 monolayers. Up to now, the main focus of the theoretically predicted TMX3 monolayers has been on the Chern insulator states, which can realize the quantum anomalous Hall effect. Here, using first-principles calculations, we theoretically demonstrate that the stable OsCl3 monolayer has a ferromagnetic ground state and a spin-polarized Dirac point without spin-orbit coupling (SOC), which disappears in the band structure of a Janus OsBr1.5Cl1.5 monolayer. We find that OsCl3 exhibits in-plane magnetization when SOC is included. By manipulating the magnetization direction along the C-2 symmetry axis of the OsCl3 structure, a gapless half-Dirac semimetal state with SOC can be achieved, which is different from the gapped Chern insulator state. Both semimetal states of OsCl3 monolayer without and with SOC exhibit a linear half-Dirac point (twofold degenerate) with high Fermi velocities. The achievement of the 2D semimetal state with SOC is expected to be found in other TMX3 monolayers, and we confirm it in a TiCl3 monolayer. This provides a different perspective to study the band structure with SOC of the 2D TMX3 family.

Abstract: Recently, interlayer and intralayer excitons in transition metal dichalcogenide heterobilayers have been studied both experimentally and theoretically. In spite of a growing interest, these layer-resolved excitons in the presence of external stimuli, such as strain, remain not fully understood. Here, using density-functional theory calculations with many-body effects, we explore the excitonic properties of vertically stacked MoSe2/WSe2 heterobilayer in the presence of in-plane biaxial strain of up to 5%. We calculate the strain dependence of exciton absorption spectrum, oscillator strength, wave function, and binding energy by solving the Bethe-Salpeter equation on top of the standard GW approach. We identify the interlayer and intralayer excitons by analyzing their electron-hole weights and spatial wave functions. We show that with the increase in strain magnitude, the absorption spectrum of the interlayer and intralayer excitons is red-shifted and re-ordered, and the binding energies of these layer-resolved excitons decrease monotonically and almost linearly. We derive the sensitivity of exciton binding energy to the applied strain and find that the intralayer excitons are more sensitive to strain than the interlayer excitons. For instance, a sensitivity of -7.9 meV/% is derived for the intra-MoSe2-layer excitons, which is followed by -7.4 meV/% for the intra-WSe2-layer excitons, and by -4.2 meV/% for the interlayer excitons. Our results indicate that interlayer and intralayer excitons in vertically stacked MoSe2/WSe2 heterobilayer are efficiently tunable by in-plane biaxial strain.

Abstract: A systematic investigation is undertaken for studying the optoelectronic properties of single crystal diamond (SCD) grown by microwave plasma chemical vapor deposition (MPCVD). It is indicated that, without intentional doping and surface treatment during the sample growth, the terahertz (THz) optical conduction in SCD is mainly affected by surface H-terminations, -OH-, O- and N-based functional groups. By using THz time-domain spectroscopy (TDS), we measure the transmittance, the complex dielectric constant and optical conductivity σ(ω) of SCD. We find that SCD does not show typical semiconductor characteristics in THz regime, where σ(ω) cannot be described rightly by the conventional Drude formula. Via fitting the real and imaginary parts of σ(ω) to the Drude-Smith formula, the ratio of the average carrier density to the effective electron mass γ = ne/m*, the electronic relaxation time τ and the electronic backscattering or localization factor can be determined optically. The temperature dependence of these parameters is examined. From the temperature dependence of γ, a metallic to semiconductor transition is observed at about T = 10 K. The temperature dependence of τ is mainly induced by electron coupling with acoustic-phonons and there is a significant effect of photon-induced electron backscattering or localization in SCD. This work demonstrates that THz TDS is a powerful technique in studying SCD which contains H-, N- and O-based bonds and has low electron density and high dc resistivity. The results obtained from this study can benefit us to gain an in-depth understanding of SCD and may provide new guidance for the application of SCD as electronic, optical and optoelectronic materials.

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: 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: 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: 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: 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: 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: 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: Motivated by a recent experiment that reported the synthesis of a new 2D material nitrogenated holey graphene (C2N) [Mahmood et al., Nat. Commun., 2015, 6, 6486], the electronic, magnetic, and mechanical properties of nitrogenated (C2N), phosphorated (C2P) and arsenicated (C2As) monolayer holey graphene structures are investigated using first-principles calculations. Our total energy calculations indicate that, similar to the C2N monolayer, the formation of the other two holey structures are also energetically feasible. Calculated cohesive energies for each monolayer show a decreasing trend going from the C2N to C2As structure. Remarkably, all the holey monolayers considered are direct band gap semiconductors. Regarding the mechanical properties (in-plane stiffness and Poisson ratio), we find that C2N has the highest in-plane stiffness and the largest Poisson ratio among the three monolayers. In addition, our calculations reveal that for the C2N, C2P and C2As monolayers, creation of N and P defects changes the semiconducting behavior to a metallic ground state while the inclusion of double H impurities in all holey structures results in magnetic ground states. As an alternative to the experimentally synthesized C2N, C2P and C2As are mechanically stable and flexible semiconductors which are important for potential applications in optoelectronics.

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