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 quantum electrostatic heterostructure (QEH) model allows for efficient computation of the dielectric screening properties of layered van der Waals (vdW)-bonded heterostructures in terms of the dielectric functions of the individual two-dimensional (2D) layers. Here, we extend the QEH model by including (1) contributions to the dielectric function from infrared active phonons in the 2D layers, (2) screening from homogeneous bulk substrates, and (3) intraband screening from free carriers in doped 2D semiconductor layers. We demonstrate the potential of the extended QEH model by calculating the dispersion of coupled phonons in multilayer stacks of hexagonal boron-nitride (hBN), the strong hybridization of plasmons and optical phonons in graphene/hBN heterostructures, the effect of substrate screening on the exciton series of monolayer MoS2, and the properties of hyperbolic plasmons in a doped phosphorene sheet. The new QEH code is distributed as a Python package with a simple command line interface and a comprehensive library of dielectric building blocks for the most common 2D materials, providing an efficient open platform for dielectric modeling of realistic vdW heterostructures.