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

Abstract: The theoretical investigation on structural, vibrational, and electronic properties of zinc-blende (ZB) AgI were carried out employing first principles density functional theory calculations. Thermoelectric properties then were predicted through semi-classical Boltzmann transport equations within the constant relaxation time approximation. Equilibrium lattice parameter, bulk modulus, elastic constants, and vibrational properties were calculated by using generalized gradient approximation. Calculated properties are in good agreement with available experimental values. Electronic and thermoelectric properties were investigated both with and without considering spin-orbit coupling (SOC) effect which is found to have a strong influence on p-type Seebeck coefficient as well as the power factor of the ZB-AgI. By inclusion of SOC, a reduction of the band-gap and p-type Seebeck coefficients as well as the power factor was found which is the indication of that spin-orbit interaction cannot be ignored for p-type thermoelectric properties of the ZB-AgI. By using deformation potential theory for electronic relaxation time and experimentally predicted lattice thermal conductivity, we obtained aZTvalue 1.69 (0.89) at 400 K for n-type (p-type) carrier concentration of 1.5 x 10(18)(4.6 x10(19)) cm(-3)that makes ZB-AgI as a promising room temperature thermoelectric material.

Abstract: The adsorption and diffusion of Li, Na, K and Ca atoms on a Mo2C monolayer are systematically investigated by using first principles methods. We found that the considered metal atoms are strongly bound to the Mo2C monolayer. However, the adsorption energies of these alkali and earth alkali elements decrease as the coverage increases due to the enhanced repulsion between the metal ions. We predict a significant charge transfer from the ad-atoms to the Mo2C monolayer, which indicates clearly the cationic state of the metal atoms. The metallic character of both pristine and doped Mo2C ensures a good electronic conduction that is essential for an optimal anode material. Low migration energy barriers are predicted as small as 43 meV for Li, 19 meV for Na and 15 meV for K, which result in the very fast diffusion of these atoms on Mo2C. For Mo2C, we found a storage capacity larger than 400 mA h g(-1) by the inclusion of multilayer adsorption. Mo2C expands slightly upon deposition of Li and Na even at high concentrations, which ensures the good cyclic stability of the atomic layer. The calculated average voltage of 0.68 V for Li and 0.30 V for Na ions makes Mo2C attractive for low charging voltage applications.

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: 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: 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: 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: Efficient thermal transport control is a fundamental issue for electronic device applications such as information, communication, and energy storage technologies in modern electronics in order to achieve desired thermal conditions. Structural defects in materials provide a mechanism to adjust the thermal transport properties of these materials on demand. In this context, the effect of structural defects on lattice thermal conductivities of two-dimensional hexagonal binary group-III nitride (XN, X = B, Al, and Ga) semiconductors is systematically investigated by means of classical molecular dynamics simulations performed with recently developed transferable inter-atomic potentials accurately describing defect energies. Here, two different Green-Kubo based approaches and another approach based on non-equilibrium molecular dynamics are compared in order to get an overall understanding. Our investigation clearly shows that defect concentrations of 3% decrease the thermal conductivity of systems containing these nitrites up to 95%. Results hint that structural defects can be used as effective adjustment parameters in controlling thermal transport properties in device applications associated with these materials. Published under an exclusive license by AIP Publishing.

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