Tailoring Nitrogen Terminals on MXene Enables Fast Charging and Stable Cycling Na-Ion Batteries at Low Temperature

An interlayer confined strategy to realize the substitution of nitrogen terminals between Ti3C2 MXene layers is proposed. The targeted Ti3C2-Nfunct anode exhibits fast-charging ability and great superiority in cycle life at − 25 °C. Ti3C2-Nfunct not only possesses a lower Na-ion diffusion energy barrier and charge transfer activation energy, but also exhibits Na+-solvent co-intercalation behavior at low temperature. An interlayer confined strategy to realize the substitution of nitrogen terminals between Ti3C2 MXene layers is proposed. The targeted Ti3C2-Nfunct anode exhibits fast-charging ability and great superiority in cycle life at − 25 °C. Ti3C2-Nfunct not only possesses a lower Na-ion diffusion energy barrier and charge transfer activation energy, but also exhibits Na+-solvent co-intercalation behavior at low temperature. Sodium-ion batteries stand a chance of enabling fast charging ability and long lifespan while operating at low temperature (low-T). However, sluggish kinetics and aggravated dendrites present two major challenges for anodes to achieve the goal at low-T. Herein, we propose an interlayer confined strategy for tailoring nitrogen terminals on Ti3C2 MXene (Ti3C2-Nfunct) to address these issues. The introduction of nitrogen terminals endows Ti3C2-Nfunct with large interlayer space and charge redistribution, improved conductivity and sufficient adsorption sites for Na+, which improves the possibility of Ti3C2 for accommodating more Na atoms, further enhancing the Na+ storage capability of Ti3C2. As revealed, Ti3C2-Nfunct not only possesses a lower Na-ion diffusion energy barrier and charge transfer activation energy, but also exhibits Na+-solvent co-intercalation behavior to circumvent a high de-solvation energy barrier at low-T. Besides, the solid electrolyte interface dominated by inorganic compounds is more beneficial for the Na+ transfer at the electrode/electrolyte interface. Compared with of the unmodified sample, Ti3C2-Nfunct exhibits a twofold capacity (201 mAh g−1), fast-charging ability (18 min at 80% capacity retention), and great superiority in cycle life (80.9%@5000 cycles) at − 25 °C. When coupling with Na3V2(PO4)2F3 cathode, the Ti3C2-Nfunct//NVPF exhibits high energy density and cycle stability at − 25 °C.


Introduction
The ever-increasing demand such as space exploration, military defence, and electric vehicles impels researchers to focus on developing high-performance energy storage devices at low temperature (low-T) [1][2][3]. Although great efforts have been made to improve the low-T performance, lithium-ion batteries (LIBs) cannot be charged at high currents at low-T and are accompanied by safety risks due to serious dendrite problems [4]. And these issues will be more prominent if commercial graphite is used as the anode. In this regard, researchers turn their attention to sodium-ion batteries (SIBs) owing to the low cost and superior temperature tolerance compared to LIBs [5]. In addition, the de-solvation energy of Na + is about 25-30% smaller than Li + , meaning a lower activation barrier of Na (de) insertion [6], which is promising to achieve faster charging and higher battery performance at low-T according to the Arrhenius formula [7]. However, their low-T performance is still limited due to the subzero-temperature induced sluggish ionic transport in electrodes, decreased ions conductivity in electrolytes, and increased impedance between electrode and electrolyte [8,9].
To solve these issues, external warming devices and selfheating systems are applied successively (Scheme 1a), which realize the low-T operation with minimal battery performance loss [10]. However, an extra high current is required for initial activation. Besides, the self-heating devices need to change the structure of the battery and lead to uneven heat distribution; thus there are certain safety risks [11]. Electrolytes are considered as a key reason for battery performance loss, due to the increased viscosity and poor ionic transport at low-T. Great efforts, such as employing low-melting-point solvents, high ions conductive electrolytes, novel salt additives and designing the solvation structure of electrolytes, are made to develop electrolytes applicable to low-T conditions [12][13][14]. Though several electrolytes (e.g., the mixture of PC and EC, ether-based electrolytes) have been found suitable for low-T operation, low-T batteries are still hindered by the slow charging process (Scheme 1b), if carbonbased anodes are employed [15,16]. Serious polarization at high current density under low-T leads to the fast arrival of cut-off voltages, resulting in minimal accessible capacity and notorious dendrites formation. As proposed by the Sand time model, the time (τ) of the appearance of dendrites is inversely proportional to the square of the current density (J) and directly proportional to ionic diffusion (D) [17][18][19].
in which μ c and μ a are the mobilities of cation and anion, and C 0 is the initial concentration of the electrolyte. Therefore, it is vital to develop novel anode materials that possess fast interfacial kinetics to conquer these challenges at low-T. Ti 3 C 2 , a typical member of the MXenes family, possessing high conductivity, large interlayer spacing, low Na + diffusion barrier, high theoretical capacity and appropriate operation voltage, is prospective anode material of SIBs at low-T [20,21]. However, the surface terminations, such as -OH and -F, highly affect the electronic conductivity and adsorb ability of cations [22,23]. Restacking of nanosheets causes poor electrolyte wettability and ion accessibility, leading to sluggish ion diffusion and limited reversible capacity [24,25]. In general, the performance of MXenes has been improved at low-T by constructing 3D architectures to enhance ionic accessibility, interlayer pillaring and surficial groups regulation [9,[26][27][28][29][30][31]. Tailoring surficial terminals on MXenes is a method to improve the charge transfer process by directly regulating electronic structure [23,32]. On the other hand, surficial terminals would promote interfacial chemical bonding between ions and substrates to suppress the formation of dendrites [33]. However, compared to the conventional gas doping methods, it is still a great challenge to find a facile and efficient strategy to tailor terminals between Ti 3 C 2 MXene layers. Moreover, the ions storage mechanism of MXenes under low temperature is still rarely reported. Herein, to lower the charge transfer barrier and boost ions diffusion kinetics, an interlayer confined strategy for tailoring nitrogen terminals on Ti 3 C 2 has been proposed, which realizes the fast charging ability at − 25 °C (Scheme 1c). The pre-intercalated cetyltrimethylammonium bromide (CTAB) not only introduces nitrogen source into Ti 3 C 2 MXene layers for designing nitrogen terminals (Ti 3 C 2 -N funct ), but also supports the layer structure during the annealing process. The corresponding structural transformation and formation mechanism are demonstrated by variable temperature in situ X-ray diffraction (XRD). As revealed by density functional theory (DFT) calculations, tailoring nitrogen terminals would lead to charge redistribution on Ti 3 C 2 layer, narrowing the bandgap, endowing the sodiophilicity, and reducing the diffusion energy barrier. The comprehensive X-ray photoelectron spectroscopy (XPS) analyses indicate that the solid electrolyte interface (SEI) formed on the Ti 3 C 2 and Ti 3 C 2 -N funct electrode shows different composition changes with more inorganic compounds in the interior when formed on Ti 3 C 2 -N funct electrode, which is essential for the Na + transfer at the electrode/electrolyte interface. In addition, the multi-scale physical characterizations show that Ti 3 C 2 -N funct might possess Na + -solvent co-intercalation behavior to circumvent the de-solvation process to achieve fast kinetics at low-T. Therefore, the designed Ti 3 C 2 -N funct anodes deliver high reversible capacity, fast-charging ability, and superior cycling stability (80.9% after 5000 cycles) at − 25 °C. When coupling with Na 3 V 2 (PO 4 ) 2 F 3 cathode, the full cells yield high energy density and high capacity retention at − 25 °C.

Preparation of Ti 3 C 2 -N funct MXene
As depicted in Fig. S1, the Ti 3 C 2 -N funct is realized by the in situ thermal decomposition of the pre-intercalated CTAB molecules between the interlayers of Ti 3 C 2 (Ti 3 C 2 -CT confined ). First, the intercalation of CTAB and its effect on the layer structure of Ti 3 C 2 is proved and analysed. The electrostatic interaction between negative-charged surface of Ti 3 C 2 and the positive-charged CTAB impels Ti 3 C 2 sheets to expand to accommodate CTAB, which guarantees the intercalation of nitrogen source in layers of Ti 3 C 2 [34][35][36][37][38]. Fourier transform infrared spectroscopy (FTIR) spectra demonstrate that the peaks corresponding to R 4 N + and -CH 2 stretching modes of cetyltrimethylammonium (CTA + ) shift to a lower wave in Ti 3 C 2 -CT confined , implying the electrostatic interaction between CTAB and Ti 3 C 2 (Fig. 1a). XRD patterns exhibit that the (002) peak of Ti 3 C 2 -CT confined become wider and shift to a lower angle compared to Ti 3 C 2 (Fig. S2). This is completely different from the simple overlaying of (002) peak and CTAB peak observed in Ti 3 C 2 -CT mix (Ti 3 C 2 physically mixed with CTAB powders), further verifying that CTAB has successfully intercalated in the layers of Ti 3 C 2 . Such a phenomenon is consistent with high-resolution transmission electron microscopy (HRTEM) images, wherein the interlayer spacing of Ti 3 C 2 enlarges to 1.45 nm from 1.13 nm after CTAB pre-intercalation (Fig. 1e, f). Second, the interlayered confined CTAB will decompose and part of the released N atoms will be captured by Ti 3 C 2 to form nitrogen terminals during the annealing process. Temperature is a crucial factor that strongly influence the formation of surface terminations on MXenes. To determine the optimal tailoring temperature, in situ high-temperature XRD and thermogravimetric analysis (TGA) were performed to track the thermal structural transformation of Ti 3 C 2 -CT confined as demonstrated in Figs. 1b and S3. When the temperature rises from 25 to 200 °C, the (002) peak of Ti 3 C 2 -CT confined slightly shifts to a higher angle, and a weight loss of 2.2% is observed both in Ti 3 C 2 and Ti 3 C 2 -CT confined , which could be ascribed to the loss of water adsorbed in the layers of Ti 3 C 2 [39]. The salient deviation of TGA curves from 200 to 250 °C between Ti 3 C 2 and Ti 3 C 2 -CT confined can be attributed to the decomposition of CTAB, in which N atoms are released from CTAB (Fig. S3). With the temperature increasing above 300 °C, the removal of surficial terminals will take place as in previous reports [29,40]. When the temperature rises to 350 °C, there is an obvious drop in the TGA curves of Ti 3 C 2 -CT confined compared to Ti 3 C 2 and pure CTAB powder. Such a phenomenon could be attributed to the low carbonization ratio of CTAB and the confined effect of MXene on CTAB. A new peak is gradually formed at 25.2° (Fig. 1b), which belongs to anatase TiO 2 (JCPDS No. 71-1166) [41]. However, the (002) and (110) peaks belonging to Ti 3 C 2 are still maintained, implying the layer structure could be kept at 350 °C. When temperature keeps rising, the intensity of (002) peak decreases and more characteristic peaks of anatase TiO 2 appear, confirming that Ti 3 C 2 would inevitably transform to anatase TiO 2 at a higher temperature. To confirm the tailoring temperature, the electrochemical performance of Ti 3 C 2 -CT confined annealing between 300 and 400 °C were compared, as shown in Fig. S4. Due to the low tailoring efficiency under 300 °C and serious oxidation of Ti 3 C 2 under 400 °C, 350 °C is considered as an appropriate temperature for tailoring nitrogen terminals on Ti 3 C 2 . Fig. 1 a FTIR spectra of CTAB powder, Ti 3 C 2 and Ti 3 C 2 -CT confined . b In situ high-temperature XRD patterns, enlarged view within 5°-8° of Ti 3 C 2 -CT confined within 25-600 °C, and enlarged view of XRD pattern at 350 °C. c XRD patterns of Ti 3 C 2 , Ti 3 C 2 -350, Ti 3 C 2 -N mix and Ti 3 C 2 -N funct . d Illustration of Ti 3 C 2 during the tailoring process. HRTEM images of e Ti 3 C 2 , f Ti 3 C 2 -CT confined and g Ti 3 C 2 -Nfunct. h Ti, O, N, and C elements mappings of Ti 3 C 2 -N funct 1 3 To illustrate the important role of CTAB pre-intercalation in the preparation of Ti 3 C 2 -N funct , XRD patterns of Ti 3 C 2 -N funct , Ti 3 C 2 -N mix (product of Ti 3 C 2 -CT mix after annealing at 350 °C), Ti 3 C 2 -350 (Ti 3 C 2 powder after annealing at 350 °C), and Ti 3 C 2 are compared (Fig. 1c). The main diffraction peaks of Ti 3 C 2 -N mix are similar to that of Ti 3 C 2 and the peak of anatase TiO 2 is also observed in Ti 3 C 2 -350. Moreover, a thick carbon layer on Ti 3 C 2 -N mix also indicated that CTAB decomposed outside of Ti 3 C 2 layers (Fig. S5). Since the introduction of CTAB, more TiO 2 peaks are observed in Ti 3 C 2 -N funct than in other samples. The split-up of (002) peak in Ti 3 C 2 -N funct further indicates that Ti 3 C 2 -N funct still keeps large interlayer spacing after the decomposition of pillared CTAB. This is consistent well with the TEM images shown in Fig. 1e-g, wherein Ti 3 C 2 shows a layer distance of 1.13 nm and expands to 1.45 nm after the CTAB intercalation, then reduces to 1.35 nm after thermal treatment. Figure 1d illustrates the structure evolution of Ti 3 C 2 during the preparation of Ti 3 C 2 -N funct and the phenomenon of (002) peak splitting. Moreover, the N element is distributed uniformly on Ti 3 C 2 -N funct (Fig. 1h). The crinkled nanosheets remain after thermal treatment (Fig.  S6), ensuring electrolyte wettability and ion accessibility. These results demonstrate that CTAB pre-intercalation and thermal treatment can achieve uniform N doping and make Ti 3 C 2 -N funct maintain well-layered structure with a large lattice distance, while the physical mixing of CTAB cannot achieve this effect.

Electronic Structure Analysis
XPS was applied to record the difference of the valence states of Ti/O elements and the existence of N between Ti 3 C 2 -CT confined and Ti 3 C 2 -N funct . The Ti 2p spectra of Ti 3 C 2 -CT confined , Ti 3 C 2 -N funct , and Ti 3 C 2 are provided in Figs. 2a and S7a, a new peak appears at 456.5 eV, which is attributed to Ti-N in Ti 3 C 2 -N funct , indicating that N atoms bond with Ti atoms on Ti 3 C 2 [42][43][44]. Moreover, the O 1s spectra reflect the change of surficial composition on Ti 3 C 2 ( Fig. 2b) [45,46]. The peak at 533.8 eV, corresponding to -ON derived from the electrostatic interaction between N atoms in CTAB and surficial O-terminals of Ti 3 C 2 , disappears in Ti 3 C 2 -N funct . In addition, the C-Ti peak of Ti 3 C 2 , Ti 3 C 2 -CT confined and Ti 3 C 2 -N funct also reflect the electrostatic interaction between CTAB and Ti 3 C 2 (Fig. S7).
Different from the existence of lattice substitution for C atoms (396.0 eV) in previous reports, the N 1s peak can only be split into two peaks located at 399.8 and 402.0 eV in this work, corresponding to the surficial N-terminals and adsorbed N atoms, labelled as N funct and N ads (Fig. 2c) [32,42]. It is clear that, after the thermal treatment, the proportion of N funct increases sharply, further verifying the successful transformation of adsorbed N atoms in CTAB to the bonding N atoms in Ti 3 C 2 -N funct . The concentrations of nitrogen element in Ti 3 C 2 -CT confined and Ti 3 C 2 -N funct are 3.63 and 1.05 at%, as exhibited in Table S1.
DFT calculation was conducted to figure out the effect of tailoring N-terminals on the electronic structure of Ti 3 C 2 . A 3 × 2 × 1 supercell of Ti 3 C 2 O 2 is constructed as the initial structure as demonstrated in Fig. S8a. An oxygen atom is substituted by a nitrogen atom in the supercell to approximate the N atoms concentration in Ti 3 C 2 -N funct (Fig. S8b), corresponding to a chemical stoichiometry of Ti 3 C 2 O 1.83 N 0.17 . Charge accumulation around N atoms and depletion near Ti atoms are verified by electron density difference (EDD) in Fig. 2d. Based on Eq. 2, the plane-averaged EDD Δρ and the charge transfer ΔQ are delivered to quantitatively compare the charge discrepancy after N substitution (Fig. 2d right) [8]. It could be observed that more charge will transfer from N atoms to Ti atoms compared to that of O atoms. The Mulliken charge distribution demonstrates that Ti atoms bond with N atoms and inner Ti atoms obtain more electrons with a reduced charge (Fig. 2d). These results are consistent with Ti 2p XPS spectra in Fig. 2a, in which the peaks of Ti-C and Ti 2+ -C slightly shift to low binding energy due to more electrons provided by N atoms. The above observations suggest that the N-terminals could induce a charge redistribution, leading to the formation of more active sites of redox reactions [44].
The density of states (DOS) of Ti 3 C 2 O 2 and Ti 3 C 2 O 1.83 N 0.17 are shown in Fig. 2e. As observed, the conduction band (CB) of Ti 3 C 2 O 2 and Ti 3 C 2 O 1.83 N 0.17 is mainly contributed by Ti 3d, while the valence band (VB) originates from the hybridization of Ti 3d, C 2p and O 2p [32]. The bandgap is reduced from 0.15 to 0.13 eV after N atom substitution, manifesting that electrons could easily migrate from VB to CB, implying the enhanced electronic conductivity of Ti 3 C 2 O 1.83 N 0.17 . The DOS shape of Ti 3 C 2 O 1.83 N 0.17 below the Fermi level is sharper, indicating that the substitution of N atoms makes electrons more localized [47], which is beneficial for the chemical bonding with Na + to increase the driving force for nucleation to inhibit the formation of metallic Na nuclei [33]. Such improvement is particularly important at low temperatures and can effectively inhibit the formation of sodium dendrites to achieve fast charging ability. This speculation will be confirmed in the next part.
Furthermore, the adsorption energy of Na atom on Ti 3 C 2 O 2 and Ti 3 C 2 O 1.83 N 0.17 is investigated to demonstrate the sodiophilicity of N-groups on Ti 3 C 2 . The optimized geometric structures of a Na atom on monolayer Ti 3 C 2 O 2 and Ti 3 C 2 O 1.83 N 0.17 are exhibited in Figs. S9 and S10. As exhibited in Fig. S11, after tailoring N-terminals on Ti 3 C 2 , the adsorption energy becomes more negative, confirming the adsorption stability of Na. It is worthy to notice that the adsorption sites near the N atoms possess the lowest adsorption energy among types of C top and Ti top, suggesting that the surficial N-terminals would induce Na + to form a Na-N-Ti interaction to keep the lowest energy state. Such a phenomenon indicates that the existence of N-terminal on Ti 3 C 2 could regulate the Na deposition behavior. Moreover,  Fig. S12, wherein the Ti 3 C 2 O 1.83 N 0.17 possesses larger interlayer spacing, in accordance with HRTEM results (Fig. 1 g). The negative formation energy suggests that Na + prefers to intercalate into the layers of Ti 3 C 2 . The result suggests that the surficial nitrogen terminals regulation could effectively improve the electronic conductivity, enhancing the sodiophilicity of Ti 3 C 2 , and improve the possibility of Ti 3 C 2 for accommodating more Na atoms, further enhancing the Na + storage capability of Ti 3 C 2 .

Electrochemical Performance at − 25 °C
Based on the well-tailored Ti 3 C 2 -N funct anode, half-cells are assembled with metallic Na to evaluate the electrochemical performance at − 25 °C. Ti 3 C 2 -N funct delivers a much higher capacity, which is nearly 2 times that of Ti 3 C 2 at 25 and − 25 °C (Fig. 3a). Besides, the capacity loss caused by the temperature drop is greatly reduced, and the reversible capacity of Ti 3 C 2 -N funct at − 25 °C is about 77% of room temperature. The first charging/discharging curves of Ti 3 C 2 and Ti 3 C 2 -N funct demonstrate that SEI forms in the first discharge process, which leads to the unpleasant initial Coulombic efficiency (Fig. S13). Figure 3b is the rate capability of the two anodes at − 25 °C, the specific capacities of Ti 3 C 2 -N funct are 201, 182, 172, 160, 143, 126, and 90 mAh g −1 at 0.05, 0.1, 0.2, 0.5, 1.0, 2.0, and 5.0 A g −1 , while about half of the capacity is obtained in Ti 3 C 2 . The Ti 3 C 2 -N funct anode delivers fast-charging/discharging ability at low-T, which could recharge to 80% capacity (160 mAh g −1 ) within 18 min and 72% capacity (144 mAh g −1 ) in 8 min (Fig. S14). After 5000 cycles at 1.0 A g −1 at − 25 °C, the Ti 3 C 2 -N funct electrode maintains capacity retention of 80.9% (Fig. 3d), superior to most reported low-T batteries as shown in Table S2 and Fig. 3c [48][49][50].
In addition, Na 3 V 2 (PO 4 ) 2 F 3 cathode, due to its high capacity retention and reversibility of Na + at low-T [51], was selected as the basis of eventual full-cell construction (Fig. 3g). The electrochemical performance of Na 3 V 2 (PO 4 ) 2 F 3 cathodes are shown in Fig. S15. The , which keeps high capacity retention (95%) compared to room temperature (Figs. 3e and S16). The designed Ti 3 C 2 -N funct //NVPF SIB displays excellent cycle stability at − 25 °C (Fig. 3f), which keeps capacity retention of 88.1% after 200 cycles at 0.5C. Moreover, as shown in Fig. S17, the Ti 3 C 2 -N funct //NVPF full cell achieves a maximum energy density and power density of 414.14 Wh kg −1 and 763.98 W kg −1 at − 25 °C, respectively (based on the cathode material). Hence, the interlayer confined strategy for tailoring N-terminals on Ti 3 C 2 is effective for improving the Na + storage performance at low-T of MXenes including temperature adaptability, fast-charging ability, and ultra-long lifespan.

Analysis of Kinetics of Ti 3 C 2 -N funct Anode at Low Temperature
Given a deeper insight into the effects of tailoring N-terminals on Ti 3 C 2 , DFT calculations and various experimental analyses were delivered to investigate the kinetics of Ti 3 C 2 and Ti 3 C 2 -N funct . To clarify the fast diffusion process of Na + in Ti 3 C 2 -N funct , DFT calculations of energy barrier of possible Na + diffusion paths in bilayer framework for Ti 3 C 2 O 2 and Ti 3 C 2 O 1.83 N 0.17 were carried out. As previous results in Fig. S11, the site on top of C atoms is preferred Na adsorption sites. Therefore, two possible pathways of Na + between the two nearest neighboring C sites in the bilayer of Ti 3 C 2 O 2 are explored (Fig. 4a). For O-Path C-C , the Na + is hopping directly to the nearest C site in a one-step path; for O-Path C-Ti-C , the Na + is migrating along the pathway from the top of C atoms to the top of Ti atoms and then to the nearest C atoms [22,52]. The calculated energy barrier of two paths are 288.7 and 132.2 meV, indicating that Na + tends to migrate along the path C-Ti-C . When an oxygen atom is substituted by a nitrogen atom in Ti 3 C 2 O 1.83 N 0.17 , C atoms with different chemical environment are emerging, labelled as C A , C B , C C , and C D . Therefore, three diffusion directions of Na + in the bilayer of Ti 3 C 2 O 1.83 N 0.17 are explored as shown in Figs. 4b, c and S18-S19. It is noteworthy that the diffusion energy barrier of path C-Ti-C in Ti 3 C 2 O 1.83 N 0.17 is lower than that of Ti 3 C 2 O 2 and the energy barrier for Na + hopping over the Ti-N bonds is lower than that of Ti-O bonds, demonstrating that tailoring N-terminals on Ti 3 C 2 is efficient to facilitate Na + diffusion. Especially, the N-Path C A -Ti 1 -C B for Na + migration is presented in Fig. 4b, with a calculated energy barrier as low as 69.9 meV, which is much lower than that of typical Na-ion insertion materials as previous reported, such as TiO 2 (2.20 eV), TiS 2 (1.20 eV), and TiSe 2 (0.5 eV), confirming the fast kinetics of Ti 3 C 2 -N funct [53][54][55]. The Na + storage behavior of Ti 3 C 2 -N funct electrodes at − 25 °C is verified by the cyclic voltammetry (CV) curves at various scan rates from 0.2 to 2.0 mV s −1 (Fig. S20a). According to Eq. 3: in which the Na + storage mechanisms can be defined as faradic ion intercalation (b = 0.5) or capacitive response (b = 1.0) [56]. The calculated b values for peak C 1 , C 2 , A 1 , and A 2 are 0.89, 0.95, 0.88, and 0.94, respectively, demonstrating that the Ti 3 C 2 -N funct electrodes possess fast kinetics at − 25 °C (Fig. S20b). The Na + diffusion coefficients and polarization of Ti 3 C 2 and Ti 3 C 2 -N funct electrodes at − 25 °C were identified by the galvanostatic intermittent titration technique (GITT) (Fig. 4d, details in Fig. S21). The Ti 3 C 2 -N funct delivers much lower polarization during the charge/discharge process than that of Ti 3 C 2 at − 25 °C, which might suppress the nucleation rate of Na dendrites as shown in Fig. 4e [18]. The calculated Na + diffusion coefficients of Ti 3 C 2 and Ti 3 C 2 -N funct electrodes at 25 and − 25 °C are demonstrated in Fig. 4f. At room temperature, the Na + diffusion coefficients of Ti 3 C 2 and Ti 3 C 2 -N funct are very close. However, while transferred to − 25 °C, Ti 3 C 2 -N funct could nearly maintain the same Na + diffusion coefficients at room temperature, which is much higher than that of Ti 3 C 2 , indicating the better adaptability of Ti 3 C 2 -N funct at low-T. This discrepancy is significant, as proposed by the Sand time model [17], the nearly 5.6 times difference in Na + diffusion coefficients is expected to prolong the time (τ) of dendrites appearing, which could guarantee the cycle stability at low-T. Therefore, beneficial from the interfacial chemical bonding of Na + and the low nucleation rate of Na dendrites, the Ti 3 C 2 -N funct delivers excellent cycle stability.
The interfacial charge transfer dominates battery performance at low-T as previous studies reported, the activation energy (E a ) holds the key to determining the interfacial kinetics [7]. Based on the Butler-Volmer equation and Arrhenius equation, the relationship between charge transfer resistance (R ct ) and E a is generalized as: in which A is a constant, R is the gas constant, T is the temperature. As shown in Fig. 4g-h, the temperature-dependent EIS was performed to measure the value of R ct from − 25 to 25 °C for Ti 3 C 2 and Ti 3 C 2 -N funct electrodes. The equivalent circuit model of Ti 3 C 2 and Ti 3 C 2 -N funct electrodes is exhibited in Fig. S22. The fitted R ct values of Ti 3 C 2 and Ti 3 C 2 -N funct electrodes at various temperatures are exhibited in Table S3. It is noteworthy that the R ct values of two electrodes are low, indicating the fast charge transfer process at − 25 °C. Based on Eq. 4, the calculated values of E a for Ti 3 C 2 and Ti 3 C 2 -N funct electrodes are 200.8 and 182.2 meV (Fig. 4i), respectively. This result demonstrates that Ti 3 C 2 -N funct could serve a nearly 10% reduction in activation energy during the interfacial charge transfer process than that of Ti 3 C 2 . All these observations further confirm that tailoring N-terminals on Ti 3 C 2 not only accelerate interfacial kinetics and Na + diffusion but also lower the charge transfer energy barrier, which is essential for achieving fastcharging ability at low-T.

Analysis of Electrode/Electrolyte Interface
SEI is a crucial composition on electrode surface to affect the ion transport, which determines the electrochemical performance at low-T. To get deep insight of the SEI composition, high-resolution XPS of C 1s, O 1s, F 1s, and Na 1s spectra were collected from the electrode surface and in-depth of 10 nm with Ar + sputtering for both Ti 3 C 2 and Ti 3 C 2 -N funct anodes (after discharging to 0.01 V at − 25 °C), as shown in Figs. 5a-c and S23. The C 1s spectra are fitted using peaks with binding energies of 284.2 (C-Ti), 284.8 (C-C), 285.6 (C-O), 286.3 (CH 2 -CF 2 ), 287.4 (C=O), 288.3 (RCH 2 -F), 289.5 (O-C=O), and 290.8 (-CF 2 -) eV, which are consistent with Ti 3 C 2 , the main reduction products of diglyme solvent, Na 2 CO 3 and PVDF [57]. These are also corresponding to the peaks located at 530.6 (C-Ti-O x ), 531.0 (Na 2 CO 3 ), 531.7 (C-Ti-(OH) x ), 532.6 (R-O-Na), and 533.8 (C=O) eV (Fig. 5b). Due to the reduction of TiO 2 , the peak at 529.6 eV (Na-O-Ti) belonging to Na x TiO 2 could be observed in O 1s spectra of Ti 3 C 2 -N funct [58]. The peak located at 685.0 eV in F 1s spectra representing Na-F bond and the Na 1s spectra shows a peak at 1072.2 eV (Na-F and Na-O), demonstrating the existence of sodium compounds (Figs. 5c and S23) [51,59]. Note that, the peak at 686.0 eV in F 1s spectra is belonging to the complex fluorosulfate, which is decomposed from NaCF 3 SO 3 salt. Such a phenomenon could also be observed in the S 2p spectra (Figs. 5d and S24). Therefore, it can be inferred that the SEI layers formed on Ti 3 C 2 and Ti 3 C 2 -N funct anodes are composed of both organic compounds (RCH 2 ONa and ROCO 2 Na) and inorganic compounds (NaF and Na 2 CO 3 ). The complex fluorosulfate is distributed near the surface of Ti 3 C 2 electrode. It is noteworthy that there is a difference between SEI formed on the two electrodes, where the organic compounds are the main components in the SEI of Ti 3 C 2 electrode and the inorganic compounds are dominated in the SEI of Ti 3 C 2 -N funct electrode.
The compositions of SEI layer vary during the charge-discharge process. Figures 5e, f and S25 demonstrate the elemental compositions of SEI at different states of charge (SOC) and depths of Ti 3 C 2 and Ti 3 C 2 -N funct anodes. For both electrodes, the C and F concentrations decrease during the discharge process and increase when charging back to 2.5 V, while the trends of the O and Na concentrations are opposite. It could be ascribed to the formation and reconstruction of SEI layer during the discharge and charge process. Moreover, there is a difference between Ti 3 C 2 and Ti 3 C 2 -N funct electrodes when discharging to 0.01 V. For the Ti 3 C 2 -N funct electrode, the C concentration decreases sharply with sputtering depths increasing, while the Na concentration increases and the O concentration decreases slightly. This phenomenon indicates that there is a thin organic compounds layer distributed near the surface and the inorganic compounds dominate the interior of the SEI. In contrast, there is no obvious change of the C and Na concentrations in the SEI of the Ti 3 C 2 electrode. This means the organic compounds are distributed homogeneously in the SEI on the Ti 3 C 2 electrode. Beneficial from the low charge transfer resistance of inorganic compounds as reported [58], the Ti 3 C 2 -N funct electrode delivers lower activation energy during the interfacial charge transfer process as illustrated in Fig. 5g, which is consistent with the above dynamic analysis.
Moreover, the Ti 2p spectra of Ti 3 C 2 and Ti 3 C 2 -N funct at different SOC after Ar + etching to a depth of 10 nm are shown in Fig. 5d. The Ti 2p spectra of Ti 3 C 2 -N funct have greatly changed after discharged to 0.01 V, wherein the ratio of Ti 2+ -C increases and the peaks of Ti 2+ -C, Ti 3+ -C, and Ti-N shift to the lower binding energy, indicating the reduced valence of Ti and the formation of Ti-O-Na and Na-N-Ti interaction [60]. After the extraction of Na + , the partial Ti 2+ is oxidized to Ti 3+ when charged back to 2.5 V, demonstrating the reversible redox reactions of the Ti 3+ / Ti 2+ couple. However, as for the Ti 3 C 2 electrode, the ratio of Ti 2+ -C and Ti 3+ -C in Ti 2p spectra nearly maintained during the discharging and charging process. Such a phenomenon demonstrates that tailoring N-terminals on Ti 3 C 2 provides more active sites for the redox reaction and induces more charge transfer during the charge/discharge process, leading to higher capacity. Overall, due to the inorganic component of SEI layer and more active sites for redox reaction, it is easy to understand the reason why Ti 3 C 2 -N funct delivers better sodium storage performance than that of Ti 3 C 2 as shown in Fig. 3.

Na + Storage Mechanism of Ti 3 C 2 -N funct at Low Temperature
How do Ti 3 C 2 -N funct participate in the electrochemical reactions at low-T? Cycling performance is used to compare the Na + storage capability of Ti 3 C 2 -N funct at different conditions (Fig. 6c). Interestingly, Ti 3 C 2 -N funct electrode delivers stable cycling performance at − 25 °C than that at room temperature. To facilitate a deep understanding of the sodiation behavior of Ti 3 C 2 -N funct , in situ and ex situ XRD measurements were made to demonstrate the structure reversibility during the charging-discharging process under different conditions. In situ XRD patterns with 2θ ranging from 6.0° to 9.0° and 22.5°-27.5° of Ti 3 C 2 -N funct operating at room temperature are displayed in Fig. 6a. According to the discharge curve, the sodiation process can be divided into surficial adsorption and intercalation. Although the (002) peak belonging to Ti 3 C 2 fixes at 8.3° during the continuous discharging-charging process, the (002′) peak shifts to a higher angle during discharging, and then could not return when get the fully charged state. Moreover, the (008') peak disappears during discharging and reemerges at the charged state, which could be attributed to the intercalation of Na + . These intercalated Na + will attract the MXene layers to each other, leading to a reduced interlayer space [61,62]. The obvious variation for the layer structure of Ti 3 C 2 -N funct during the charging-discharging process would lead to unstable cycling performance at room temperature. The reflections located between 22.5° and 27.5° are attributed to the sodiation process of TiO 2 , in which TiO 2 is transformed to sodium titanates (JCPDS No. 78-1590), titanium suboxide (JCPDS No. 72-2101) [58]. However, the TiO 2 particles forming during the tailoring process didn't contribute to the increased capacity of Ti 3 C 2 -N funct as shown in Fig. S26. In addition, TiO 2 -N mix (deriving from the annealing product of TiO 2 mixed with CTAB powders) exhibits an unpleasant performance than pure TiO 2 , further proving that the improved performance of Ti 3 C 2 -N funct exactly stemmed from tailoring surficial N-terminals rather than the contribution of TiO 2 (Fig. S27).
To explore the Na + storage mechanism at low-T, ex situ XRD patterns of Ti 3 C 2 -N funct electrodes operating at − 25 °C were recorded (Fig. 6b). Interestingly, after aging for 1 h in Na-ion half-cell at − 25 °C, the potential of Ti 3 C 2 -N funct reduced 0.05 V (Fig. S28a). It could be attributed to the chemical pre-sodiation of Ti 3 C 2 -N funct during the aging process, leading to the broaden of (002) peak, the shift of (004) peak of Ti 3 C 2 -N funct , and the disappearance of peaks of anatase TiO 2 (Fig. S28b). It has been reported that the anatase TiO 2 would transform into amorphous structure when sodiated at low temperature [63]. The widened (002) peak shifts to a lower angle when firstly discharged to 0.8 V, and then it returns to the higher angle at the fully discharged state. Two new lattice peaks are detected at 28.5° and 31.4°, corresponding to Na 2 CO 3 and NaO 2 . Different from operating at room temperature, the interlayer spacing of Ti 3 C 2 -N funct is enlarged before discharging to 0.8 V, and then it shrinks during subsequent discharging process at low-T. Such phenomenon could be correlated to the intercalation of the solvent molecules [61,64]. Moreover, the slight the variation of (002) peak indicates the stable structure of Fig. 6 a In situ XRD patterns of Ti 3 C 2 -N funct operating at room temperature with 2θ ranging from 6.0° to 9.0° and 22.5°-27.5°, the first two charge/discharge curves corresponding to the in situ XRD patterns. b Ex situ XRD patterns of Ti 3 C 2 -N funct electrodes operating at − 25 °C. c Comparison of cycling performance for Ti 3 C 2 -N funct operating at room temperature and − 25 °C. d Ex situ FTIR spectra of Ti 3 C 2 -N funct electrodes. e Illustration of Na + -solvent co-intercalation behavior in Ti 3 C 2 -N funct at low-T Ti 3 C 2 -N funct , leading to outstanding cycling performance at low-T. The elemental mapping of Ti 3 C 2 -N funct at a deepdischarged state is exhibited in Fig. S29 further prove the Na + intercalation at low-T.
It is believed that the desolvation process is the rate-determining step at low-T due to the high desolvation energy barrier [65,66]. As previous studies have indicated that there is an ion-solvent co-intercalation behavior in graphitic structure when operated in the ether-based electrolytes with high solvation energy [67,68], we speculate that the graphitelike MXenes possessing larger interlayer spacing might deliver the similar ion-solvent co-intercalation behavior to circumvent the desolvation process to achieve fast kinetics at low-T. To confirm the hypothesis, the ex situ FTIR spectra of Ti 3 C 2 -N funct electrodes during sodiation and desodiation are collected. The FTIR spectra of the electrolyte (1 M NaCF 3 SO 3 in diglyme) and pure solvent (diglyme) are compared in Fig. S30, demonstrating that the peak at 858 cm −1 belonging to the solvated Na + [69]. The characteristic peaks of the solvated Na + appear at the full discharged state and disappear when charged back to 2.5 V in the ex situ FTIR spectra of Ti 3 C 2 -N funct and Ti 3 C 2 (Figs. 6d and S31). Our previous work shows that if the solvent is simply adsorbed on the surface of the material, its signal is independent of the state of charge and could be detected during the whole charging-discharging process [70]. In Ti 3 C 2 -N funct , the characteristic peaks of diglyme solvent are strongly associated with SOC, suggesting that Na + and solvent might possess similar migration behavior. Therefore, integrating these evidences into account, we deduced that Ti 3 C 2 -N funct possesses Na + -solvent co-intercalation behavior during the charge transfer process as illustrated in Fig. 6e, which could avoid the high desolvation energy barrier to realize the fastcharging ability at low-T.

Conclusions
In summary, we propose and demonstrate that tailoring nitrogen terminals on Ti 3 C 2 through the interlayer confined strategy is crucial to enable high-performance SIBs at low temperature. N atoms derived from the decomposition of confined CTAB molecules directly substitute the surface terminals, which tailor the in-plane structure of Ti 3 C 2 . The interfacial kinetics and energy storage mechanism at − 25 °C of Ti 3 C 2 -N funct are investigated. It is found that tailoring nitrogen terminals could boost Na + diffusion kinetics and lower charge transfer barrier by the synergistic effects of large interlayer spacing, charge redistribution, and strong adsorption, empowering Ti 3 C 2 -N funct with higher Na + diffusion coefficient and a 10% reduction in activation energy at low-T. The inorganic compounds in the SEI on Ti 3 C 2 -N funct are beneficial for Na + transfer. Moreover, the ion-solvent co-intercalation behavior endows Ti 3 C 2 -N funct with fastcharging ability at low-T. As expected, the Ti 3 C 2 -N funct anodes deliver high capacity retention, fast-charging ability (charging 80% capacity within 18 min), and ultra-long lifespan (5000 cycles with a capacity retention of 80.9%) at − 25 °C, far exceeding that of pristine Ti 3 C 2 . The assembled Ti 3 C 2 -N funct //NVPF full cells also deliver high energy density and cycling stability at − 25 °C. This work opens avenues for the development of other 2D materials in constructing high-energy storage systems at low temperatures.