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: We present here a new interatomic potential parameter set to predict the thermal conductivity of zirconium trisulfide monolayers. The generated Tersoff-type force field is parameterized using data collected with first-principles calculations. We use non-equilibrium molecular dynamics simulations to predict the thermal conductivity. The generated parameters result in very good agreement in structural, mechanical, and dynamical parameters. The room temperature lattice thermal conductivity ( kappa) of the considered crystal is predicted to be kappa x x = 25.69Wm – 1K – 1 and kappa y y = 42.38Wm – 1K – 1, which both agree well with their corresponding first-principles values with a discrepancy of less than 5%. Moreover, the calculated kappa variation with temperature (200 and 400 K) are comparable within the framework of the accuracy of both first-principles and molecular dynamics simulations.

Abstract: The rational design of two-dimensional (2D) piezoelectric materials has recently garnered great interest due to their increasing use in technological applications, including sensor technology, actuating devices, energy harvesting, and medical applications. Several materials possessing high piezoelectric response have been reported so far, but a high-throughput first-principles approach to estimate the piezoelectric potential of layered materials has not been performed yet. In this study, we systematically investigated the piezoelectric (e(11), d(11)) and elastic (C-11 and C-12) properties of 128 thermodynamically stable 2D semiconductor materials by employing first-principle methods. Our high-throughput approach demonstrates that the materials containing Group-V elements produce significantly high piezoelectric strain constants, d(11) > 40 pm V-1, and 49 of the materials considered have the e(11) coefficient higher than MoS2 insomuch as BrSSb has one of the largest d(11) with a value of 373.0 pm V-1. Moreover, we established a simple empirical model in order to estimate the d(11) coefficients by utilizing the relative ionic motion in the unit cell and the polarizability of the individual elements in the compounds.

Abstract: High-performance rechargeable batteries are becoming very important for high-end technologies with their ever increasing application areas. Hence, improving the performance of such batteries has become the main bottleneck to transferring high-end technologies to end users. In this study, we propose an argon intercalation strategy to enhance battery performance via engineering the interlayer spacing of honeycomb structures such as graphite, a common electrode material in lithium-ion batteries (LIBs). Herein, we systematically investigated the LIB performance of graphite and hexagonal boron nitride (h-BN) when argon atoms were sent into between their layers by using first-principles density-functional-theory calculations. Our results showed enhanced lithium binding for graphite and h-BN structures when argon atoms were intercalated. The increased interlayer space doubles the gravimetric lithium capacity for graphite, while the volumetric capacity also increased by around 20% even though the volume was also increased. The ab initio molecular dynamics simulations indicate the thermal stability of such graphite structures against any structural transformation and Li release. The nudged-elastic-band calculations showed that the migration energy barriers were drastically lowered, which promises fast charging capability for batteries containing graphite electrodes. Although a similar level of battery promise was not achieved for h-BN material, its enhanced battery capabilities by argon intercalation also support that the argon intercalation strategy can be a viable route to enhance such honeycomb battery electrodes.

Abstract: We report a high-throughput first-principles characterization of the structural, mechanical, electronic, and vibrational properties of tetragonal single-layer transition metal Xenes (t-TMXs). Our calculations revealed 22 dynamically, mechanically and chemically stable structures among the 96 possible free-standing layers present in the t-TMX family. As a fingerprint for their structural identification, we identified four characteristic Raman active phonon modes, namely three in-plane and one out-of-plane optical branches, with various intensities and frequencies depending on the material in question. Spin-polarized electronic calculations demonstrated that anti-ferromagnetic (AFM) metals, ferromagnetic (FM) metals, AFM semiconductors, and non-magnetic semiconductor materials exist within this family, evidencing the potential of t-TMXs for further use in multifunctional heterostructures.

Abstract: Two-dimensional heterostructures are an emerging class of materials for novel applications because of extensive engineering potential by tailoring intriguing properties of different layers as well as the ones arising from their interface. A systematic investigation of mechanical, electronic, and optical properties of possible heterostructures formed by bilayer structures graphenelike ZnO and MgO monolayers is presented. Different functionality of each layer makes these heterostructures very appealing for device applications. ZnO layer is convenient for electron transport in these structures, while MgO layer improves electron collection. At the outset, all of the four possible stacking configurations across the heterostructure are mechanically stable. In addition, stability analysis using phonon dispersion reveals that the AB stacking formed by placing the Mg atom on top of the O atom of the ZnO layer is also dynamically stable at zero temperature. Henceforth, we have investigated the optical properties of these stable heterostructures by applying many-body perturbation theory within the framework of GW approximation and solving the Bethe-Salpeter equation. It is demonstrated that strong excitonic effects reduce the optical band gap to the visible light spectrum range. These results show that this new two-dimensional form of ZnO/MgO heterostructures open an avenue for novel optoelectronic device applications.

Abstract: There is an increasing need to understand interfaces between two-dimensional materials to realize an energy efficient boundary with low contact resistance and small heat dissipation. In this respect, we investigated the impact of charge and substitutional atom doping on the electronic transport properties of the hybrid metallic-semiconducting lateral junctions, formed between metallic (1T and 1T(d)) and semiconducting (1H) phases of MoS2 by means of first-principles and non-equilibrium Green function formalism based calculations. Our results clearly revealed the strong influence of the type of interface and crystallographic orientation of the metallic phase on the transport properties of these systems. The Schottky barrier height, which is the dominant mechanism for contact resistance, was found to be as large as 0.63 eV and 1.19 eV for holes and electrons, respectively. We found that armchair interfaces are more conductive as compared to zigzag termination due to the presence of the metallic Mo zigzag chains that are directed along the transport direction. In order to manipulate these barrier heights we investigated the influence of electron doping of the metallic part (i.e. 1T(d) -MoS2). We observed that the Fermi level of the hybrid system moves towards the conduction band of semiconducting 1H-MoS2 due to filling of 4d-orbital of metallic MoS2, and thus the Schottky barrier for electrons decreases considerably. Besides electron doping, we also investigated the effect of substitutional doping of metallic MoS2 by replacing Mo atoms with either Re or Ta. Due to its valency, Re (Ta) behaves as a donor (acceptor) and reduces the Schottky barrier for electrons (holes). Since Re and Ta based transition metal dichalcogenides crystallize in either the 1T(d) or 1T phase, substitutional doping with these atom favors the stabilization of the 1T(d) phase of MoS2. Co-doping of hybrid structure results in an electronic structure, which facilities easy dissociation of excitons created in the 1H part.

Abstract: The surface termination of MXenes greatly determines the electrochemical properties and ion kinetics on their surfaces. So far, hydroxyl-, oxygen-, and fluorine-terminated MXenes have been widely studied for energy storage applications. Recently, sulfur-functionalized MXene structures, which possess low diffusion barriers, have been proposed as candidate materials to enhance battery performance. We performed first-principles calculations on the structural, stability, electrochemical, and ion dynamic properties of Li-adsorbed sulfur-functionalized groups 3B, 4B, 5B, and 6B transition-metal (M)-based MXenes (i.e., M2CS2 with M = Sc, Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, and W). We performed phonon calculations, which indicated that all of the above M2CS2 MXenes, except for Sc, are dynamically stable at T = 0 K. The ground-state structure of each M2CS2 monolayer depends on the type of M atom. For instance, while sulfur prefers to sit at the FCC site on Ti2CS2, it occupies the HCP site of Cr-based MXene. We determined the Li adsorption configurations at different concentrations using the cluster expansion method. The highest maximum open-circuit voltages were computed for the group 4B element (i.e., Ti, Zr, and Hf)-based M2CS2, which are larger than 2.1 V, while their average voltages are approximately 1 V. The maximum voltage for the group 6B element (i.e., Cr, Mo, W)-based M2CS2 is less than 1 V, and the average voltage is less than 0.71 V. We found that S functionalization is helpful for capacity improvements over the O-terminated MXenes. In this respect, the computed storage gravimetric capacity may reach up to 417.4 mAh/g for Ti2CS2 and 404.5 mAh/g for V2CS2. Ta-, Cr-, Mo-, and W-based M2CS2 MXenes show very low capacities, which are less than 100 mAh/g. The Li surface diffusion energy barriers for all of the considered MXenes are less than 0.22 eV, which is favorable for high charging and discharging rates. Finally, ab initio molecular dynamic simulations performed at 400 K and bond-length analysis with respect to Li concentration verify that selected promising systems are robust against thermally induced perturbations that may induce structural transformations or distortions and undesirable Li release.

Abstract: The elastic constant C-11 and piezoelectric stress constant e(1),(11) of two-dimensional (2D) dielectric materials comprising h-BN, 2H-MoS2, and other transition-metal dichalcogenides and dioxides are calculated using lattice dynamical theory. The results are compared with corresponding quantities obtained with ab initio calculations. We identify the difference between clamped-ion and relaxed-ion contributions with the dependence on inner strains which are due to the relative displacements of the ions in the unit cell. Lattice dynamics allows us to express the inner-strain contributions in terms of microscopic quantities such as effective ionic charges and optoacoustical couplings, which allows us to clarify differences in the piezoelectric behavior between h-BN and MoS2. Trends in the different microscopic quantities as functions of atomic composition are discussed.

Abstract: Thermal expansion behavior of two-dimensional (2D) nitrides and graphene were studied by ab initio molecular dynamics (MD) simulations as well as quasiharmonic approximation (QHA). Anharmonicity of the acoustic phonon modes are related to the unusual negative thermal expansion (NTE) behavior of the nitrides. Our results also hint that direct ab initio MD simulations are a more elaborate method to investigate thermal expansion behavior of 2D materials than the QHA. Nevertheless, giant NTE coefficients are found for h-GaN and h-AlN within the covered temperature range 100-600 K regardless of the chosen computational method. This unusual NTE of 2D nitrides is reasoned with the out-of-plane oscillations related to the rippling behavior of the monolayers.

Abstract: The effect of the number of stacking layers and the type of stacking on the electronic and optical properties of bilayer and trilayer black phosphorus are investigated by using first-principles calculations within the framework of density functional theory. We find that inclusion of many-body effects (i.e., electron-electron and electron-hole interactions) modifies strongly both the electronic and optical properties of black phosphorus. While trilayer black phosphorus with a particular stacking type is found to be a metal by using semilocal functionals, it is predicted to have an electronic band gap of 0.82 eV when many-body effects are taken into account within the G(0)W(0) scheme. Though different stacking types result in similar energetics, the size of the band gap and the optical response of bilayer and trilayer phosphorene are very sensitive to the number of layers and the stacking type. Regardless of the number of layers and the type of stacking, bilayer and trilayer black phosphorus are direct band gap semiconductors whose band gaps vary within a range of 0.3 eV. Stacking arrangements that are different from the ground state structure in both bilayer and trilayer black phosphorus exhibit significant modified valence bands along the zigzag direction and result in larger hole effective masses. The optical gap of bilayer (trilayer) black phosphorus varies by 0.4 (0.6) eV when changing the stacking type. The calculated binding energy of the bound exciton hardly changes with the type of stacking and is found to be 0.44 (0.30) eV for bilayer (trilayer) phosphorous.

Abstract: Direct experimental measurement of thermal expansion coefficient without substrate effects is a challenging task for two-dimensional (2D) materials, and its accurate estimation with large-scale ab initio molecular dynamics is computationally very expensive. Machine learning-based interatomic potentials trained with ab initio data have been successfully used in molecular dynamics simulations to decrease the computational cost without compromising the accuracy. In this study, we investigated using Gaussian approximation potentials to reproduce the density functional theory-level accuracy for graphene within both lattice dynamical and molecular dynamical methods, and to extend their applicability to larger length and time scales. Two such potentials are considered, GAP17 and GAP20. GAP17, which was trained with pristine graphene structures, is found to give closer results to density functional theory calculations at different scales. Further vibrational and structural analyses verify that the same conclusions can be deduced with density functional theory level in terms of the reasoning of the thermal expansion behavior, and the negative thermal expansion behavior is associated with long-range out-of-plane phonon vibrations. Thus, it is argued that the enabled larger system sizes by machine learning potentials may even enhance the accuracy compared to small-size-limited ab initio molecular dynamics.

Abstract: Piezoelectricity is a unique material property that allows one to convert mechanical energy into electrical one or vice versa. Transition metal dichalcogenides (TMDC) and transition metal dioxides (TMDO) are expected to have great potential for piezoelectric device applications due to their noncentrosymmetric and two-dimensional crystal structure. A detailed theoretical investigation of the piezoelectric stress (e 11 ) and piezoelectric strain (d 11 ) coefficients of single layer TMDCs and TMDOs with chemical formula MX 2 (where M= Cr, Mo, W, Ti, Zr, Hf, Sn and X = O, S, Se, Te) is presented by using first-principles calculations based on density func- tional theory. We predict that not only the Mo- and W-based members of this family but also the other materials with M= Cr, Ti, Zr and Sn exhibit highly promising piezoelectric properties. CrTe 2 has the largest e 11 and d 11 coefficients among the group VI elements (i.e., Cr, Mo, and W). In addition, the relaxed-ion e 11 and d 11 coefficients of SnS 2 are almost the same as those of CrTe 2 . Furthermore, TiO 2 and ZrO 2 pose comparable or even larger e 11 coefficients as compared to Mo- and W-based TMDCs and TMDOs. Our calculations reveal that TMDC and TMDO structures are strong candidates for future atomically thin piezoelectric applications such as transducers, sensors, and energy harvesting devices due to their piezoelectric coefficients that are comparable (even larger) to currently used bulk piezoelectric materials.

Abstract: Metallic contacts are critical components of electronic devices and the presence of a large Schottky barrier is detrimental for an optimal device operation. Here, we show by using first-principles calculations that a self-assembled monolayer (SAM) of polar molecules between the metal electrode and MoSe2 monolayer is able to convert the Schottky contact into an almost Ohmic contact. We choose -CH3 and -CF3 terminated short-chain alkylthiolate (i.e. SCH3 and fluorinated alkylthiolates (SCF3)) based SAMs to test our approach. We consider both high (Au) and low (Sc) work function metals in order to thoroughly elucidate the role of the metal work function. In the case of Sc, the Fermi level even moves into the conduction band of the MoSe2 monolayer upon SAM insertion between the metal surface and the MoSe2 monolayer, and hence possibly switches the contact type from Schottky to Ohmic. The usual Fermi level pinning at the metal-transition metal dichalcogenide (TMD) contact is shown to be completely removed upon the deposition of a SAM. Systematic analysis indicates that the work function of the metal surface and the energy level alignment between the metal electrode and the TMD monolayer can be tuned significantly by using SAMs as a buffer layer. These results clearly indicate the vast potential of the proposed interface engineering to modify the physical and chemical properties of MoSe2.

Abstract: MXenes are an emerging class of two-dimensional materials, which in their thinnest limit consist of a monolayer of carbon or nitrogen (X) sandwiched between two transition metal (M) layers. We have systematically searched for superconductivity among MXenes for a range of transition metal elements, based on a full first-principles characterization in combination with the Eliashberg formalism. Thus, we identified six superconducting MXenes: three carbides (Mo2C, W2C and Sc2C) and three nitrides (Mo2N, W2N and Ta2N). The highest critical temperature of similar to 16 K is found in Mo2N, for which a successful synthesis method has been established [Urbankowskiet al.,Nanoscale, 2017,9, 17722-17730]. Moreover, W2N presents a novel case of competing superconducting and charge density wave phases.