3.1. Electrolyte Solvation Structure Analysis
cMD simulations were carried out to investigate the specified solvation structure caused by the asymmetric structure on the LiDFTFSI-EC/EMC. The solvation structure of an electrolyte is determined by its molecular structure, which can be distinguished into three categories [
32,
42]: (1) solvents surrounding Li
+ (SSL), (2) Li
+-anion single pairs (LASPs) and (3) Li
+-anion clusters (LACs). Strong Li
+-solvent interaction should be avoided, but it inhibits desolvation and triggers solvent cointercalation during Li
+ reduction while charging [
43]. The comparisons of sovlation structures between LiDFTFSI and LiTFSI are as shown in
Figure 1a and
Figure S3a. Based on Li
+ coordination with oxygen, the proportions of SSL, LASP, and LAC in the LiDFTFSI system are 35%, 27.5%, and 37.5%, respectively, while in the LiTFSI system, these are 47.5%, 27.5%, and 25%, respectively. Thus, the proportion of SSL in the LiDFTFSI system is much lower than that in the LiTFSI system. The significant increase in the proportion of LAC solvation structures suggests a shift in the solvation structure from SSL to LASP and LAC.
In the LiDFTFSI system, a greater number of anions penetrate the solvation sheath and contribute to the first solvation sheath of Li
+. To understand such an effect, DFT calculations were used to obtain the electrostatic potential (ESP) energy surface of TFSI
– and DFTFSI
–. In
Figure 1b, the (-CF
2H) substituent of DFTFSI
– differs significantly from the (-CF
3) substituent of TFSI
–, resulting in a less negative charge and enhancing anion-Li
+ interaction. Meanwhile, DFTFSI
– facilitates hydrogen bonds with solvent molecules, in turn, enhancing ionic conductivity. A considerable amount of DFTFSI
– enters the first solvation sheath, displacing solvent molecules and enriching the solvation structure with DFTFSI
–.
Table S1 provides additional statistics on Li
+ coordination. Compared to the LiTFSI system, the LiDFTFSI system exhibits a higher average anion coordination number and a lower average solvent coordination number. This anion-rich solvation structure can facilitate salt preferential decomposition, resulting in the formation of an inorganic-rich SEI interface.
The radial distribution function (RDF) provides further insights into the detailed structure features of anions and solvent molecules within the solvation sheath. As shown in
Figure 1c, Li
+-O(DFTFSI
–), Li
+-O(EC) and Li
+-O(EMC) peaks at 0.208, 0.216 and 0.212 nm, respectively, indicating that DFTFSI
– has strong coordination with Li
+, followed by EMC and EC, further confirming that DFTFSI
– penetrates into the first solvation sheath. In LiTFSI, Li
+-O(TFSI
–), Li
+-O(EC) and Li
+-O(EMC) peaks at 0.212, 0.212 and 0.210 nm, respectively. Compared to LiDFTFSI, the anion peak exhibit a systematical right-shift, and EMC has larger coordination than TFSI
– and EC (
Figure S3b). The coordination number of Li
+-O(DFTFSI
–) increases compared to that of Li
+-O(TFSI
–), and the coordination number distribution plot also shows a decrease in the coordination number of Li
+-O(EC) in the LiDFTFSI system. These observations indicate that DFTFSI
– weakens the interaction between EC and Li
+ and enters the solvation sheath, thereby modifying the solvation structure, consistent with previous report [
44].
3.2. Underlying mechanism of electrolyte reduction and SEI formation
Hybrid molecular dynamics simulations were carried out to investigate the decomposition and the initial SEI formation of LiDFTFSI-EC/EMC electrolyte. Electrolyte reduction is crucial for SEI formation at the lithium metal anode, a 2.86 ns HAIR simulation using the model shown in
Figure S2b to reveal the detailed SEI reaction. As shown in
Figure 2, DFTFSI
– decomposed within 7 ps AIMD simulation. The detailed reactions are as follows: the N‒S bond on ‒CF
2H side breaks at 0.19 ps (
Figure 2a), indicating that the asymmetric substituent increases the reactivity. Subsequently, the C‒F bonds on ‒CF
2H quickly break at 1.21 ps due to the formation of LiF as reduced by Li
0. As shown in
Figure 2a, the S‒‒O bond breaks at 1.51 ps, along with Li
2O formation. After extended simulation, the dehydrogenation reaction occurs at 72.60 ps, leading to LiH formation, consistent with experimental report. These results demonstrate that introducing hydrogen atoms to the asymmetric group increases the salt reactivity, and produces ionic conductive LiH and mechanically stable LiF.
On the other hand, the ‒CF
3 group detached from the salt at 3.20 ps through C‒S bond cleavage, similar to the previous report [
37]. Deep reduction of DFTFSI
–, followed by S‒O, N‒S, and C‒F cleavages, contributed to the inorganic components of the SEI (
Figure 2b). The breaking of the N‒S bond occurred at 5.26 ps, and all F atoms were released from DFTFSI
– at 6.88 ps (
Figure 2b). This detailed mechanism demonstrated that asymmetric DFTFSI
– can rapidly react with Li
0 and generate products, including LiF, Li
2O, LiH, and Li
3N, which can potentially improve the ionic conductivity, regulate the mechanical properties of the SEI and inhibit dendrite formation [
6]. The decomposition mechanism of DFTFSI
– revealed by AIMD simulation supports the experimental hypothesis, as illustrated in
Figure 2c and
Figure S4.
The solvent decomposition mechanism is crucial. AIMD simulations indicate that EC is reduced during charging, contributing to SEI formation. EC starts to decompose within 50 ps, and three different initial pathways can be distinguished based on the reactive trajectories. At 2.51 ps, in a Li
+-rich environment, EC undergoes a one-electron reduction, breaking its C‒O bond and transforming from a ring to a chain structure, forming OC
2H
4OCO
–. Subsequently, OC
2H
4OCO
– follows two pathways: it accepts another electron to form OC
2H
4O
2–, releasing one CO at 4.00 ps (
Figure 3a), and releases CO
2 while forming the free radical anion •C
2H
4O
–, consistenting with previous DFT [
45]. At 61.21 ps, •C
2H
4O
– is further reduced by Li
0 and forms Li
2O along with one electron transfer (
Figure 3b). In addition to one-electron reduction, two-electron reduction was also observed in the simulation. EC underwent a two-electron reduction to form CO
32– and C
2H
4, as previously observed in experiment. Further reductions occurred at 50.60 and 56.50 ps, resulting in the formation of Li
2O from CO
32– (
Figure 3c).
The decomposition of EMC is slower than EC and generally falls into two categories. As shown in
Figure 3d, the C‒O bond breaks with Li
2O at 50.6 ps, and then the C‒O bond breaks, leading to the formation of CO and CH
3O
– at 55.69 ps. The other reaction path involves the initial detachment of the ‒OCH
3 group from the carbon chain via one-electron reduction, followed by the breaking of the C‒O bond in the remaining ‒OCOCH
2CH
3 species, yielding CO and C
2H
5O
– (
Figure 3e). EMC decomposition did not occur in the first 50 ps of the AIMD simulation but began in the early stages of the HAIR simulation, revealing similar reaction mechanism.
After initial decomposition, deeper reactions continue, forming the main products of the SEI. The products and the primitive SEI formed after 2.86 ns of HAIR simulation are shown in
Figures S5 and S6. The inner SEI layer, close to the anode surface, predominantly contains inorganic compounds like Li
2O, LiF, etc. [
46,
47]. These inorganic compounds improve mechanical properties, inhibit lithium dendrite growth, decrease electrolyte consumption, and promote a uniform and stable SEI, thereby potentially enhancing battery cycle performance. Instead, the outer SEI layer, close to the electrolyte, consists mainly of carbon-oxygen organic compounds from solvent decomposition, consistent with previous experiment [
35,
36]. To further analyze the products and atomic structure of the SEI, the XPS spectra of C 1 s, O 1 s and F 1 s are simulated, as shown in
Figure 4. C atoms exhibit a complex chemical environment characterized by various bond types, including C‒O, C‒C, C‒H, and C=O bonds (
Figure 4a), primarily arising from the decomposition of solvents into RC‒O (e.g., CH
3OLi, C
2H
5OLi). In
Figure 4b, O atoms are either bonded to C, forming C‒O and C=O bonds, which confirm the presence of carbon-oxygen groups and correlate with the C 1 s XPS spectrum, or exist as Li
2O. The schematic diagram of DFTFSI
– and solvent decomposition shows that F atoms bond with Li
+ to form LiF (
Figure 4c). The theoretical XPS predictions well agree with the experimental results, demonstrating the reliability of the simulations.
3.3. DFTFSI–Passivated Al
In addition to regulating the anode interphase, LiDFTFSI, as the reactive active component, can remarkably suppress the dissolution of the Al
0 current collector at high potentials (>4.2 V versus Li/Li
+) in comparison with LiTFSI. The Al
0 current collector maintains good stability within the battery, because Al foil naturally forms a protective Al
2O
3 layer [
52]. However, when the voltage exceeds the threshold, Al
2O
3 can be destroyed in the organic solution, exposing fresh Al
0. This leads to progressive corrosion of the Al
0 current collector. The asymmetric structure of LiDFTFSI and the high reactivity of the ‒CF
2H group rapidly dissociate to form AlF
3, which then protects the anode. To test such hypothesis, we examined the defluorination process on the ‒CF
2 group on the Al metal surface, as shown in
Figure 5. DFT calculations show that the Gibbs free energy change (Δ
G) for this process is -1.64 eV, compared to -0.94 eV for LiTFSI. The results indicate that LiDFTFSI defluorinates more readily on the Al electrode, facilitating the subsequent formation of AlF
3, well validating our hypothesis. Therefore, the protective mechanism of LiDFTFSI on Al foil can be described as follows: LiDFTFSI diffuses and adsorbs onto the Al foil, decomposes, and the resulting F ions react with Al
3+ to form an AlF
3 passivation layer. This layer prevents further exposure of the Al metal and suppresses interfacial side reactions.