Mechanistic Insight into Ion Transport in Layered Materials for Sodium-Ion Rechargeable Batteries

Abstract

As a feasible alternative to Li-ion batteries (LIBs), Na-ion batteries (SIBs) are an emerging energy storage solution for electric vehicles, electronic devices, and large grid level energy storage systems. The TMD (transition metal dichalcogenides) component MoS2 and the NASICON structured-type NaxTi2(PO4)3 are two materials of interest as next-generation SIB materials. MoS2 is a 2D layered material with a graphene-like structure that can be used as a high-capacity anode or Na-ion host material as well as an electrocatalyst. Our experiments on the Al-doped 1T phase of MoS2 on reduced graphene oxide (rGO) nanosheets support demonstrated good electrochemical performance, such as 450, 400, etc. mAh/g at current densities of 0.05, 0.1 A/g. Experiments further revealed that Na-ion intercalation supports increased interlayer distance in 1T Al-MoS2@rGO, resulting in a stable phase structure. Density functional theory (DFT)-based simulations are utilized to study the intercalation properties of Na+ ions, in addition to GITT, EIS, and CV analysis. In this study, the Na-ion diffusion mechanism for 1T Al-MoS2@rGO interface and Al-MoS2-MoS2 interlayer host structure is explored. According to the experiments, the host structure is doped with 5% Al at a specific stable site, and Na-ions are introduced to the interface and interlayer (an interlayer distance of 9.3 Å is chosen, as estimated by the experiments). The Na-ion is hopped in both structures, and activation energies (Eact) are calculated to be 57.9 kJ/mol and 40.3 kJ/mol in the 1T Al-MoS2@rGO interface at migration distances of 0-2 Å and 2-4 Å, respectively. The related Na-ion diffusion coefficients are computed using an Arrhenius type fit to be in the order of 10-13 cm2/s and 10-10 cm2/s. Similarly, the Na-ion activation energy barriers in the 1T Al-MoS2-MoS2 interlayer are estimated to be between 49.8 and 59.6 kJ/mol, with related diffusivity ranging between 10-11 and 10-13 cm2/s. The outcomes of the DFT simulations are consistent with those of our experiments. The Na-ion transport properties in the NaTi2(PO4)3 bulk structure via the Na-ion vacancy assisted mechanism are investigated using MD simulations and DFT. The self-diffusion coefficients (DNa+) of bulk are examined, and then the lattice strain effect is proposed to alter its transport properties by introducing the lattice strain in uniaxial (x, y, z) and biaxial (xy, xz, yz) directions, and the changes are compared for strained and unstrained structure. In comparison to the unstrained structure at 323K, a significant improvement by 1-2 order change in DNa+ values (10-8 cm2/s - 10-10 cm2/s) is measured. DFT is also used to investigate the diffusion mechanisms along the various planes (line [001], [100], and [010]) of the bulk NaTi2(PO4)3 structure and Eact is determined along those directions, respectively. The values obtained are consistent with the experimental value ranges obtained for similar materials. In addition, a mixed ion-system is represented by inserting Li-ions at the surface of a NaTi2(PO4)3 structure, yielding the formula LixNay-xTi2(PO4)3. The sluggish kinetics of Li-ion in the NaTi2(PO4)3 structure were proposed by MD analysis. This corresponds well with experimental data that the NaTi2(PO4)3 framework exhibits preferred Na-ion insertion when cycling in mixed-ion electrolytes of varied concentrations, implying slow diffusion of Li-ions in the bulk Li+/Na+ mixed-ion crystal.