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: 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 structural stability of 2D transition metal dichalcogenide (TMD) formations is of particular importance for their reliable device performance in nano-electronics and opto-electronics. Recent observations show that the CVD-grown TMD monolayers are likely to encounter stability problems such as cracking or fracturing when they are kept under ambient conditions. Here, two different growth configurations are investigated and a favorable growth geometry is proposed, which also sheds light onto the growth mechanism and provides a solution for the stability and fracture formation issues for TMDs specifically for MoS2 monolayers. It is shown that 18 months naturally and thermally aged MoS2 monolayer flakes grown using specifically developed conditions, retain their stability. To understand the mechanism of the structural deterioration, two possible effective mechanisms, S vacancy defects and growth-induced tensile stress, are assessed by the first principle calculations where the role of S vacancy defects in obtaining oxidation resistant MoS2 monolayer flakes is revealed to be rather more critical. Hence, these simulations, time-dependent observations and thermal aging experiments show that durability and stability of 2D MoS2 flakes can be controlled by CVD growth configuration.

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.