Abstract: The investigation of thermal transport properties of novel two-dimensional materials is crucially important in order to assess their potential to be used in future technological applications, such as thermoelectric power generation. In this respect, the lattice thermal transport properties of the monolayer structures of group VA elements (P, As, Sb, Bi, PAs, PSb, PBi, AsSb, AsBi, SbBi, P3As1, P3Sb1, P1As3, and As3Sb1) with a black phosphorus like puckered structure were systematically investigated by first-principles calculations and an iterative solution of the phonon Boltzmann transport equation. Phosphorene was found to have the highest lattice thermal conductivity, , due to its low average atomic mass and strong interatomic bonding character. As a matter of course, anisotropic was obtained for all the considered materials, owing to anisotropy in frequency values and phonon group velocities calculated for these structures. However, the determined linear correlation between the anisotropy in the values of P, As, and Sb is significant. The results corresponding to the studied compound structures clearly point out that thermal (electronic) conductivity of pristine monolayers might be suppressed (improved) by alloying them with the same group elements. For instance, the room temperature of PBi along the armchair direction was predicted to be as low as 1.5 W m(-1) K-1, whereas that of P was predicted to be 21 W m(-1) K-1. In spite of the apparent differences in structural and vibrational properties, we peculiarly revealed an intriguing correlation between the values of all the considered materials as = c(1) + c(2)/m(2), in particular along the zigzag direction. Furthermore, our calculations on compound structures clearly showed that the thermoelectric potential of these materials can be improved by suppressing their thermal properties. The presence of ultra-low values and high electrical conductivity (especially along the armchair direction) makes this class of monolayers promising candidates for thermoelectric applications.

Abstract: Using first-principles calculations, we demonstrated that high lithium storage capacity and fast kinetics are achieved for Ti3C2O2 by preintercalating organic molecules. As a proof-of-concept, two different quinone molecules, namely 1,4-benzoquinone (C6H4O2) and tetrafluoro-1,4-benzoquinone (C6F4O2) were selected as the molecular linkers to demonstrate the feasibility of this interlayer engineering strategy for energy storage. As compared to Ti3C2O2 bilayer without linker molecules, our pillared structures facilitate a much faster ion transport, promising a higher charge/discharge rate for Li. For example, while the diffusion barrier of a single Li ion within pristine Ti3C2O2 bilayer is at least 1.0 eV, it becomes 0.3 eV in pillared structures, which is comparable and even lower than that of commercial materials. At high Li concentrations, the calculated diffusion barriers are as low as 0.4 eV. Out-of-plane migration of Li ions is hindered due to large barrier energy with a value of around 1-1.35 eV. Concerning storage capacity, we can only intercalate one monolayer of Li within pristine Ti3C2O2 bilayer. In contrast, pillared structures offer significantly higher storage capacity. Our calculations showed that at least two layers of Li can be intercalated between Ti3C2O2 layers without forming bulk Li and losing the pillared structure upon Li loading/unloading. A small change in the in-plane lattice parameters (<0.5%) and volume (<1.0%) and ab initio molecular dynamics simulations prove the stability of the pillared structures against Li intercalation and thermal effects. Intercalated molecules avoid the large contraction/expansion of the whole structure, which is one of the key problems in electrochemical energy storage. Pillared structures allow us to realize electrodes with high capacity and fast kinetics. Our results open new research paths for improving the performance of not only MXenes but also other layered materials for supercapacitor and battery applications.

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: Molybdenum carbide has various applications for which studying the material using classical molecular dynamics simulations would be valuable. Here, we develop an empirical potential within the Tersoff formalism using particle swarm optimization for the orthorhombic phase of Mo2C. The developed potential is shown to predict lattice constants, elastic properties, and equation of state results that are consistent with current and previously reported results from experiments and first principles calculations. We demonstrate the potential with simulations of indentation using multiple indenter sizes that load and unload in three different directions relative to the crystallographic lattice of orthorhombic Mo2C. Direction-dependent force-displacement trends are analyzed and explained in terms of the spatial distributions of stress and strain within the material during indentation. This study reveals the anisotropic elasticity of orthorhombic Mo2C and, more generally, provides researchers with a new empirical potential that can be used to explore the properties and behavior of the material going forward.

Abstract: Recently, an extremely low lattice thermal conductivity value has been reported for the alkali-based telluride material BaIn2Te4. The value is comparable with low-thermal conductivity metal chalcogenides, and the glass limit is highly intriguing. Therefore, to shed light on this issue, we performed first-principles phonon thermal transport calculations. We predicted highly anisotropic lattice thermal conductivity along different directions via the solution of the linearized phonon Boltzmann transport equation. More importantly, we determined several different factors as the main sources of the predicted ultralow lattice thermal conductivity of this crystal, such as the strong interactions between low-frequency optical phonons and acoustic phonons, small phonon group velocities, and lattice anharmonicity indicated by large negative mode Gruneisen parameters. Along with thermal transport calculations, we also investigated the electronic transport properties by accurately calculating the scattering mechanisms, namely the acoustic deformation potential, ionized impurity, and polar optical scatterings. The inclusion of spin-orbit coupling (SOC) for electronic structure is found to strongly affect the p-type Seebeck coefficients. Finally, we calculated the thermoelectric properties accurately, and the optimal ZT value of p-type doping, which originated from high Seebeck coefficients, was predicted to exceed unity after 700 K and have a direction averaged value of 1.63 (1.76 in the y-direction) at 1000 K around 2 x 1020 cm-3 hole concentration. For n-type doping, a ZT around 3.2 x 1019 cm-3 concentration was predicted to be a direction-averaged value of 1.40 (1.76 in the z-direction) at 1000 K, mostly originating from its high electron mobility. With the experimental evidence of high thermal stability, we showed that the BaIn2Te4 compound has the potential to be a promising mid- to high-temperature thermoelectric material for both p-type and n-type systems with appropriate doping.

Abstract: Single-layer, large-scale two-dimensional material growth is still a challenge for their wide-range usage. Therefore, we carried out a comprehensive study of monolayer MoS2 growth by CVD investigating the influence of growth zone configuration and precursors ratio. We first compared the two commonly used approaches regarding the relative substrate and precursor positions, namely, horizontal and face-down configurations where facedown approach is found to be more favorable to obtain larger flakes under identical growth conditions. Secondly, we used different types of substrate holders to investigate the influence of the Mo and S vapor confinement on the resulting diffusion environment. We suggest that local changes of the S to Mo vapor ratio in the growth zone is a key factor for the change of shape, size and uniformity of the resulting MoS2 formations, which is also confirmed by performing depositions under different precursor ratios. Therefore, to obtain continuous monolayer films, the S to Mo vapor ratio is needed to be kept within a certain range throughout the substrate. As a conclusion, we obtained monolayer triangles with a side length of 90 mu m and circles with a diameter of 500 mu m and continuous films with an area of 85 0 mu m x 1 cm when the S-to-Mo vapor ratio is optimized. (C) 2017 The Japan Society of Applied Physics

Abstract: Electrenes, an atomically thin form of layered electrides, are very recent members of the 2D materials family. In this work, we employed first principle calculations to determine stable, exfoliatable, and application-promising 2D electrene materials among possible M2X compounds, where M is a group II-A metal and X is a nonmetal element (C, N, P, As, and Sb). The promise of stable electrene compounds for battery applications is assessed via their exfoliation energy, adsorption properties, and migration energy barriers toward relevant Li, Na, K, and Ca atoms. Our calculations revealed five new stable electrene candidates in addition to previously known Ca2N and Sr2N. Among these seven dynamically stable electrenes, Ba2As, Ba2P, Ba2Sb, Ca2N, Sr2N, and Sr2P are found to be very promising for either K or Na ion batteries due to their extremely low migration energy barriers (5-16 meV), which roughly demonstrates 105 times higher mobility than graphene and two to four times higher mobility than other promising 2D materials such as MXene (Mo2C).

Abstract: Recently, a group of MAX phases, (Mo2/3Y1/3)(2)AlC, (Mo2/3Sc1/3)(2)AlC, (W2/3Sc1/3)(2)AlC,(W2/3Y1/3)(2)AlC, and (V-2/3 Zr-1/3)(2)AlC, with in-plane ordered double transition metals, named iMAX phases, have been synthesized. Experimentally, some of these MAX phases can be chemically exfoliated into two-dimensional (2D) single- or multilayered transition metal carbides, so-called MXenes. Accordingly, the 2D nanostructures derived from iMAX phases are named iMXenes. Here we investigate the structural stabilities and electronic structures of the experimentally discovered iMAX phases and their possible iMXene derivatives. We show that the iMAX phases and their pristine, F, or OH-terminated iMXenes are metallic. However, upon 0 termination, (Mo2/3Y1/3)(2)C, (Mo2/3Sc1/3)(2)C, (W2/3Y1/3)(2)C, and (W2/3Sc1/3)(2)C iMXenes turn into semiconductors. Owing to the absence of centrosymmetry, the semiconducting iMXenes may find applications in piezoelectricity. Our calculations reveal that the semiconducting iMXenes possess giant piezoelectric coefficients as large as 45 x 10(-)(10) C/m.

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: 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 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: The atomically laminated layered ternary transition-metal borides (the MAB phases) have demonstrated outstanding properties and have been applied in various fields. Understanding their thermal and mechanical properties is critical to determining their applicability in various fields such as high-temperature applications. To achieve this, we conducted first-principles calculations based on density-functional theory and the quasiharmonic approximation to determine the thermal expansion coefficients, Gruneisen parameters, bulk moduli, hardness, thermal conductivity, electron-phonon coupling parameters, and the structural and vibrational properties of MoAlB, WAlB, Cr2AlB2, and Ti2InB2. We found varying degrees of anisotropy in the thermal expansion and mechanical properties in spite of similarities in their crystal structures. MoAlB has a mild degree of anisotropy in its thermal expansion coefficient (TEC), while Cr2AlB2 and WAlB display the highest level of TEC anisotropy. We assessed various empirical models to calculate hardness and thermal conductivity, and correlated the calculated values with the material properties such as elastic moduli, Gruneisen parameter, Debye temperature, and type of bonding. Owing to their higher Gruneisen parameters, implying a greater degree of anharmonicity in lattice vibrations and lower phonon group velocities, MoAlB and WAlB have significantly lower lattice thermal conductivity values than those of Cr2AlB2 and Ti2InB2. The hardness and lattice thermal conductivity of MAB phases can be predicted with high accuracy if one utilizes an appropriate model.

Abstract: Two-dimensional materials (2DMs) continue to attract a lot of attention, particularly for their extreme flexibility and superior thermal properties. Molecular dynamics simulations are among the most powerful methods for computing these properties, but their reliability depends on the accuracy of interatomic interactions. While first principles approaches provide the most accurate description of interatomic forces, they are computationally expensive. In contrast, classical force fields are computationally efficient, but have limited accuracy in interatomic force description. Machine learning interatomic potentials, such as Gaussian Approximation Potentials, trained on density functional theory (DFT) calculations offer a compromise by providing both accurate estimation and computational efficiency. In this work, we present a systematic procedure to develop Gaussian approximation potentials for selected 2DMs, graphene, buckled silicene, and h-XN (X = B, Al, and Ga, as binary compounds) structures. We validate our approach through calculations that require various levels of accuracy in interatomic interactions. The calculated phonon dispersion curves and lattice thermal conductivity, obtained through harmonic and anharmonic force constants (including fourth order) are in excellent agreement with DFT results. HIPHIVE calculations, in which the generated GAP potentials were used to compute higher-order force constants instead of DFT, demonstrated the first-principles level accuracy of the potentials for interatomic force description. Molecular dynamics simulations based on phonon density of states calculations, which agree closely with DFT-based calculations, also show the success of the generated potentials in high-temperature simulations.

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