Abstract: The interface heat transfer of two layers induced by van der Waals (vdW) contacts is theoretically investigated, based on first-principles calculations at low temperatures. The results suggest that out-of-plane acoustic phonons with low frequencies dominate the interface thermal transport due to the vdW interaction. The interface thermal conductivity is proportional to the cubic of temperature at very low temperatures, but becomes linearly proportional to temperature as temperature increases. We show that manipulating the strain alters vdW coupling, leading to increased interfacial thermal conductivity at the interface. Our findings provide valuable insights into the interface heat transport in vdW heterostructures and support further design and optimization of electronic and optoelectronic nanodevices based on vdW contacts. The heat transfer induced by van der Waals contacts is dominated by ZA phonons. The interface thermal conductivity is proportional to the cubic of temperature, but becomes linearly proportional to temperature as temperature increases.

Abstract: Topological semimetal phases in two-dimensional (2D) materials have gained widespread interest due to their potential applications in novel nanoscale devices. Despite the growing number of studies on 2D topological nodal lines (NLs), candidates with significant topological features that combine nontrivial topological semimetal phase with superconductivity are still rare. Herein, we predict Al2B2 and AlB4 monolayers as new 2D nonmagnetic Dirac nodal line semimetals with several novel features. Our extensive electronic structure calculations combined with analytical studies reveal that, in addition to multiple Dirac points, these 2D configurations host various highly dispersed NLs around the Fermi level, all of which are semimetal states protected by time-reversal and in-plane mirror symmetries. The most intriguing NL in Al2B2 encloses the K point and crosses the Fermi level, showing a considerable dispersion and thus providing a fresh playground to explore exotic properties in dispersive Dirac nodal lines. More strikingly, for the AlB4 monolayer, we provide the first evidence for a set of 2D nonmagnetic open type-II NLs coexisting with superconductivity at a rather high transition temperature. The coexistence of superconductivity and nontrivial band topology in AlB4 not only makes it a promising material to exhibit novel topological superconducting phases, but also a rather large energy dispersion of type-II nodal lines in this configuration may offer a platform for the realization of novel topological features in the 2D limit.

Abstract: The interface-dependent electronic, vibrational, piezoelectric, and optical properties of van der Waals heterobilayers, formed by buckled GeO (b-GeO) and Janus MoSO structures, are investigated by means of first-principles calculations. The electronic band dispersions show that O/Ge and S/O interface formations result in a type-II band alignment with direct and indirect band gaps, respectively. In contrast, O/O and S/Ge interfaces give rise to the formation of a type-I band alignment with an indirect band gap. By considering the Bethe-Salpeter equation (BSE) on top of G(0)W(0) approximation, it is shown that different interfaces can be distinguished from each other by means of the optical absorption spectra as a consequence of the band alignments. Additionally, the low- and high-frequency regimes of the Raman spectra are also different for each interface type. The alignment of the individual dipoles, which is interface-dependent, either weakens or strengthens the net dipole of the heterobilayers and results in tunable piezoelectric coefficients. The results indicate that the possible heterobilayers of b-GeO/MoSO asymmetric structures possess various electronic, optical, and piezoelectric properties arising from the different interface formations and can be distinguished by means of various spectroscopic techniques.

Abstract: In the present work, the identification of two hexagonal phases of germanium oxides (namely GeO2 and GeO) through the vibrational and optical properties is reported using density functional theory calculations. While structural optimizations show that single-layer GeO2 and GeO crystallize in 1T and buckled phases, phonon band dispersions reveal the dynamical stability of each structure. First-order off-resonant Raman spectral predictions demonstrate that each free-standing single-layer possesses characteristic peaks that are representative for the identification of the germanium oxide phase. On the other hand, electronic band dispersion analysis shows the insulating and large-gap semiconducting nature of single-layer GeO2 and GeO, respectively. Moreover, optical absorption, reflectance, and transmittance spectra obtained by means of G(0)W(0)-BSE calculations reveal the existence of tightly bound excitons in each phase, displaying strong optical absorption. Furthermore, the excitonic gaps are found to be at deep UV and visible portions of the spectrum, for GeO2 and GeO crystals, with energies of 6.24 and 3.10 eV, respectively. In addition, at the prominent excitonic resonances, single-layers display high reflectivity with a zero transmittance, which is another indication of the strong light-matter interaction inside the crystal medium.

Abstract: Pure hydrogen production via water splitting is an ideal strategy for producing clean and sustainable energy. Two-dimensional (2D) cadmium chalcogenide single-layers with a tetragonal crystal structure, namely Tetra-CdX (X = S, Se, and Te) monolayers, are theoretically predicted by means of density functional theory (DFT). Their structural stability and electronic and optical properties are investigated. We find that Tetra-CdX single-layers are thermodynamically stable. Their stability decreases as we go down the 6A group in the periodic table, i.e., from X = S to Se, and Te which also means that the electronegativity decreases. All considered novel monolayers are indirect band gap semiconductors. Using the HSE06 functional the electronic band gaps of CdS, CdSe, and CdTe monolayers are predicted to be 3.10 eV, 2.97 eV, and 2.90 eV, respectively. The impact of mechanical strain on the physical properties was studied, which indicates that compressive strain increases the band gap and tensile strain decreases the band gap. The optical properties of the Tetra-CdX monolayers show the ability of these monolayers to absorb visible light. Due to the suitable band gaps and band edge positions of Tetra-CdX, these newly discovered 2D materials are promising for photocatalytic water splitting.

Abstract: High-performance thermoelectric materials are critical in recuperating the thermal losses in various machinery and promising in renewable energy applications. In this respect, the search for novel thermoelectric materials has attracted considerable attention. In particular, low dimensional materials have been proposed as potential candidates due to their unique and controllable thermal and electronic transport properties. The considerable potential of several two-dimensional materials as thermoelectric devices has already been uncovered and many new candidates that merit further research have been suggested. In this regard, we comprehensively investigate the thermoelectric coefficients and electronic fitness function (EFF) of a large family of structurally isotropic and anisotropic two-dimensional layered materials using density functional theory combined with semi-classical Boltzmann transport theory. With this high-throughput screening, we bring to light additional 2D crystals that haven't been previously classified as favorable TE materials. We predict that Pb2Se2, GeS2, As-2, NiS2, Hf2O6, Zr2O6, AsBrS, ISbTe, ISbSe, AsISe, and AsITe are promising isotropic thermoelectric materials due to their considerably high EFF values. In addition to these materials, Hf2Br4, Zr2Br4, Hf2Cl4, Zr2Cl4, Hf2O6, Zr(2)O(6)and Os(2)O(4)exhibit strong anisotropy and possess prominently high EFF values.

Abstract: Using the simplified linear combination of atomic orbitals (LCAO) method in combination with ab initio calculations, we construct a tight-binding (TB) model for two different crystal structures of monolayer gallium: a(100)- and b(010)-Gallenene. The analytical expression for the Hamiltonian and numerical results for the overlap matrix elements between different orbitals of the Ga atoms and for the Slater and Koster (SK) integrals are obtained. We find that the compaction of different structures affects significantly the formation of the orbitals. The results for a(100)-Gallenene can be very well explained with an orthogonal basis set, while for b(010)-Gallenene we have to assume a non-orthogonal basis set in order to construct the TB model. Moreover, the transmission properties of nanoribbons of both monolayers oriented along the AC and ZZ directions are also investigated and it is shown that both AC- and ZZ-b(010)-Gallenene nanoribbons exhibit semiconducting behavior with zero transmission while those of a(100)-Gallenene nanoribbons are metallic.

Abstract: MXenes are the newest class of two-dimensional (2D) materials, and they offer great potential in a wide range of applications including electronic devices, sensors, and thermoelectric and energy storage materials. In this work, we combined the outstanding electrical conductivity, that is essential for battery applications, of graphene with MXene monolayers (M2CX2 where M = Sc, Ti, V and X = OH, O) to explore its potential in Li battery applications. Through first principles calculations, we determined the stable stacking configurations of M2CX2/graphene bilayer heterostructures and their Li atom intercalation by calculating the Li binding energy, diffusion barrier and voltage. We found that: (1) for the ground state stacking, the interlayer binding is strong, yet the interlayer friction is small; (2) Li binds more strongly to the O-terminated monolayer, bilayer and heterostructure MXene systems when compared with the OHterminated MXenes due to the H+ induced repulsion to the Li atoms. The binding energy of Li decreases as the Li concentration increases due to enhanced repulsive interaction between the positively charged Li ions; (3) Ti2CO2/graphene and V2CO2/graphene heterostructures exhibit large Li atom binding energies making them the most promising candidates for battery applications. When fully loaded with Li atoms, the binding energy is -1.43 eV per Li atom and -1.78 eV per Li atom for Ti2CO2/graphene and V2CO2/graphene, respectively. These two heterostructures exhibit a nice compromise between storage capacity and kinetics. For example, the diffusion barrier of Li in Ti2CO2/graphene is around 0.3 eV which is comparable to that of graphite. Additionally, the calculated average voltages are 1.49 V and 1.93 V for Ti2CO2/graphene and V2CO2/graphene structures, respectively; (4) a small change in the in-plane lattice parameters (<1%), interatomic bond lengths and interlayer distances (<0.5 angstrom) proves the stability of the heterostructures against Li intercalation, and the impending phase separation into constituent layers and capacity fading during charge-discharge cycles in real battery applications; (5) as compared to bare M2CX2 bilayers, M2CX2/graphene heterostructures have lower molecular mass, offering high storage capacity; (6) the presence of graphene ensures good electrical conductivity that is essential for battery applications. Given these advantages, Ti2CO2/graphene and V2CO2/graphene heterostructures are predicted to be promising for lithium-ion battery applications.

Abstract: The newest members of a two-dimensional material family, involving transition metal carbides and nitrides (called MXenes), have garnered increasing attention due to their tunable electronic and thermal properties depending on the chemical composition and functionalization. This flexibility can be exploited to fabricate efficient electrochemical energy storage (batteries) and energy conversion (thermoelectric) devices. In this study, we calculated the Seebeck coefficients and lattice thermal conductivity values of oxygen terminated M2CO2 (where M = Ti, Zr, Hf, Sc) monolayer MXene crystals in two different functionalization configurations (model-II (MD-II) and model-III (MD-III)), using density functional theory and Boltzmann transport theory. We estimated the thermoelectric figure-of-merit, zT, of these materials by two different approaches, as well. First of all, we found that the structural model (i.e. adsorption site of oxygen atom on the surface of MXene) has a paramount impact on the electronic and thermoelectric properties of MXene crystals, which can be exploited to engineer the thermoelectric properties of these materials. The lattice thermal conductivity kappa(l), Seebeck coefficient and zT values may vary by 40% depending on the structural model. The MD-III configuration always has the larger band gap, Seebeck coefficient and zT, and smaller kappa(l) as compared to the MD-II structure due to a larger band gap, highly flat valence band and reduced crystal symmetry in the former. The MD-III configuration of Ti2CO2 and Zr2CO2 has the lowest kappa(l) as compared to the same configuration of Hf2CO2 and Sc2CO2. Among all the considered structures, the MD-II configuration of Hf2CO2 has the highest kappa(l), and Ti2CO2 and Zr2CO2 in the MD-III configuration have the lowest kappa(l). For instance, while the band gap of the MD-II configuration of Ti2CO2 is 0.26 eV, it becomes 0.69 eV in MD-III. The zT(max) value may reach up to 1.1 depending on the structural model of MXene.

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: First-principles density functional theory (DFT) calculations were carried out to investigate the structural and electronic properties of beryllium (Be) doped and Be and boron (B) co-doped graphene systems. We observed that not only the concentration of impurity atoms is important to tune the band-gap to some desired level, but also the specific substitution sites play a key role. In our system, which consists of 32 atoms, a maximum of 4Be and, in the co-doped state, 2Be and 3B atom substitutions are investigated. Both dopants are electron deficient relative to C atoms and cause the Fermi level to shift downward (p-type doping). A maximum band gap of 1.44 eV can be achieved on incorporation of 4Be atoms. The introduction of Be is more sensitive in terms of geometry and stability than B. However, in opening the energy gap, Be is more effective than B and N (nitrogen). Our results offer the possibility to modify the band-gap of graphene sufficiently for utilization in diverse electronic device applications.

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: We present a high-throughput computing scheme based on density functional theory (DFT) to generate a class of oxides and screen them with the aim of identifying those that might be electronically appropriate for transparent conducting oxide (TCO) applications. The screening criteria used are a minimum band gap to ensure sufficient transparency, a band edge alignment consistent with easy n- or p-type dopability, and a minimum thermodynamic phase stability to be experimentally synthesizable. Following this scheme we screened 23 binary and 1518 ternary bixbyite oxides in order to identify promising candidates, which can then be a subject of an in-depth study. The results for the known TCOs are in good agreement with the reported data in the literature. We suggest a list of several new potential TCOs, including both n- and p-type compounds.

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