P.
Jankowski
*abc,
W.
Wieczorek
ac and
P.
Johansson
bc
aFaculty of Chemistry, Warsaw University of Technology, ul. Noakowskiego 3, 00-664 Warsaw, Poland. E-mail: piotr.jankowski@chalmers.se
bDepartment of Physics, Chalmers University of Technology, SE-412 96, Gothenburg, Sweden
cALISTORE-ERI European Research Institute, 33 rue Saint Leu, 80039, Amiens, France
First published on 19th May 2016
A new anion design concept, based on combining a boron atom as the central atom and conjugated systems as ligands, is presented as a route for finding alternative Li-salts for lithium-ion batteries. The properties of a wide range of novel anions designed in this way have been evaluated by DFT calculations focusing on three different fundamental success factors/measures: the strength of the cation–anion interaction, ultimately determining both the solubility and the ionic conductivity, the oxidation limit, determining their possible use vs. high voltage cathodes, and the reduction stability, revealing a possible role of the anion in the SEI-formation at the anode. For a few anions superior properties vs. today's existing or suggested anions are predicted, especially the very low cation–anion interaction strengths are promising features. The design route itself is shown to be versatile in determining the correlation between different choices of ligands and the resulting overall properties – where the most striking feature is the decreased lithium cation interaction energy upon using the (1Z,3Z)-buta-1,3-diene-1,2,3,4-tetracarbonitrile ligands. This also opens avenues for the further design of novel anions beyond those with a boron central atom.
Luckily, many of the above requisite basic properties can be successfully estimated by standard computational chemistry methods, enabling us to search for new promising Li-salts/anions by in silico screening methods and thereby saving time and money by complementing/replacing costly exploratory synthesis.7–10 For instance, the salt dissociation degree can be evaluated as the interaction energy between the lithium cation and the anion, the ion-pair dissociation energy, which strongly depends on the charge distribution in the anion.11 To achieve the lowest energy of dissociation by the most extensive delocalization of the negative charge within the anion, strongly electron withdrawing groups are applied to the anion structures; e.g. –F, –CF3, –CN, >CO, and >SO2.12
Also the ESWs can be evaluated,10,13 and while the oxidation potential limit of the anion unquestionably should be as high as possible to accommodate high voltage cathode active materials, there is no similar strict demand for the reduction potential limit. Even the most used solvents, the organic carbonates, undergo reduction processes at potentials of ca. 1 V vs. Li+/Li0. Therefore, it is essential to provide a kinetic stability against further solvent reduction by forming a stable solid electrolyte interphase (SEI) layer on the anode surface. A salt with the ability to decompose in a controlled manner and create a SEI layer is therefore highly desirable14 – exemplified by lithium salts with anions based on boron central atoms such as lithium bis(malonato)borate (LiBMB) and lithium bis(oxalato)borate (LiBOB).15 Also the lithium bis(trifluoromethanesulfonylimide) (LiTFSI) salt forms adequate SEI layers, but only at very high concentrations.16 Moreover, even the thermal stability can be assessed, including the prediction of decomposition products, and has provided clear explanations for the thermal decomposition mechanism of LiPF6-based electrolytes.17 Recently, COSMO18 approaches have been used for electrolyte studies to predict properties such as viscosities, solubilities and melting points.19
Previously, many variations of the imide (N(SO2CnF2n+1)2−),12 including e.g. the TFSI and the bis(fluorosulfonyl)imide (FSI) anions, and the Hückel (CnN5−n(CF3)m(CN)n−m−),20,21 including e.g. the 4,5-dicyano-1,2,3-triazolate (DCTA) and the 4,5-dicyano-2-(trifluoromethyl)imidazolide (TDI) anions, families of anions have been in silico screened as PF6− alternatives. So far, some experiments suggest that LiPF6 indeed could be replaced by one of these three main concepts, as represented by, for example, LiTFSI,22 LiBOB23 or LiTDI.24 They all exhibit higher thermal stabilities, are resistant to hydrolysis, and provide sufficient conductivities of their electrolytes. Application of cyano rather than fluorinated ligands, as for TDI, also makes for less environmentally demanding salts – as we avoid fluorine synthesis. Each of these salts, however, has some fundamental drawbacks and common to all is a higher cost of production. Furthermore, LiTFSI does not properly passivate the aluminum current collector, causing its corrosion and degradation during the entire life of a LIB,25 while LiBOB shows low solubility in organic carbonate solvents, a critical parameter especially at low operation temperatures.26 LiTDI in turn, has problems in creating a stable SEI.27
Herein, we apply in silico screening to a new anion design concept, combining two current leading ideas: a tetrahedral boron atom as the anion core28–31 and a conjugated system as the anion ligand;24,32 here composed of parts of imidazole and pyrrole rings of Hückel anions. All the resulting anions are extensively compared with a wide range of other relevant anions: (i) tetrahedral boron-centered anions: BF4−, B(CH3)4−, BPh4−, BISON [B(CN)4−], BPFPB [(C2(CF3)4O2)2B−],23 BBB [(1,2-C6H4O2)2B−],33 4F-BBB [(1,2-C6F4O2)2B−],34 BCB [(C5O5)2B−],35 BNB [(2,3-C10H6O2)2B−],34 BBPB [(2,2′-C12H6O2)2B−],36 BSB [(C7H4O3)2B−],37 BOB38 and BMB;23 (ii) conjugated anions: DCTA,39 TIM (C6N5−),40 TCP (C8N5−),40 PCPI [C3(CN)5−],40 TDI,41 TDBI (C9N4H2CF3−)42 and TDPI (C7N6CF3−)43 and (iii) standard LIB electrolyte anions: PF6−, AsF6−, and TFSI. All anion data are compared by three different measures: the cation–anion interaction strength, stability vs. oxidation, and stability vs. reduction.
The electronic dissociation energy (ΔEd) moving from a Li+–anion ion-pair to separated ions was calculated for all anions using an approach outlined before.44,45 Initially the geometry of each anion was optimized and subsequently a lithium cation was manually inserted at a wide range of possible coordination sites including mono-, bi- and tri-dentate coordination, creating up to 15 different ion-pairs for each anion, and the geometries of all the resulting ion-pairs were optimized, all using the B3LYP functional46 and the 6-311+G(d) basis set. The basis set superposition error (BSSE) was neglected in the ion-pair calculations as it was previously found to be typically only a few kJ mol−1.44 The obtained electron densities were used for an analysis of the charge distribution using the atoms-in-molecules (AIM) theory.47
To assess the oxidative electrochemical stability limit (Eox) the vertical ionization potentials (IPs) (ΔEv's) were calculated as the differences between the anion and the neutral radical electronic energies, without any geometry relaxation of the radical i.e. change from the anion geometry (the Franck–Condon approximation). For a proper assessment of the reduction limit (Ered) via the electron affinity (EA) it proved necessary to apply a more accurate approach, where the reduced species was relaxed and also the solvent influence was considered by applying the C-PCM SCRF solvation model with all solvent parameters for water, foremost ε = 78. Water is chosen as a simple proxy for a large dielectric solvent. Due to the possibility of ion-pairs being reduced rather than the anions, the reduction potentials were also calculated for the most stable contact ion-pairs. All potentials were adjusted towards the Li+/Li0 reference by −1.46 V.48 For all electrochemical stability limit calculations the M06-2X functional49 was used, as it was previously found to provide the best predictive power.50 This functional was also used to calculate the chemical hardness (η) to evaluate the reactivity of the new species. As defined by Pearson,51 the chemical hardness is the average of the vertical IP and the vertical EA: η = (IP − EA)/2. All computations were made using the Gaussian09 program.52
Anion | ΔEd [kJ mol−1] B3LYP | ΔEv [V vs. Li+/Li0] M06-2X |
η [eV]
M06-2X |
E
red![]() |
E
red![]() |
---|---|---|---|---|---|
BDB | 467 | 5.80 | 4.48 | 2.00 | 1.84 |
DFDB | 530 | 5.00 | 4.58 | 1.85 | 2.17 |
DCDB | 478 | 5.87 | 4.61 | 2.02 | 2.36 |
BTMB | 503 | 5.73 | 4.97 | 1.44 | 1.99 |
DFTMB | 550 | 5.19 | 5.11 | 1.40 | 1.72 |
DCTMB | 498 | 6.03 | 5.09 | 1.59 | 1.98 |
BFB | 551 | 5.12 | 5.34 | 0.66 | 1.01 |
DFFB | 578 | 4.63 | 5.50 | 0.59 | 0.98 |
DCFB | 518 | 5.54 | 5.48 | 0.83 | 0.87 |
BTCB | 454 | 4.88 | 3.66 | 1.56 | 1.67 |
DFTCB | 479 | 3.90 | 3.76 | 1.63 | 1.75 |
DCTCB | 474 | 4.29 | 3.66 | 1.87 | 1.92 |
BBTCB | 426 | 5.44 | 3.75 | 2.67 | 2.82 |
DFBTCB | 477 | 4.93 | 3.79 | 2.48 | 2.62 |
DCBTCB | 463 | 5.44 | 3.75 | 2.62 | 2.79 |
BF4− | 602 | 6.85 | 6.94 | −1.32 | −0.84 |
B(CH3)4− | 613 | 2.83 | 4.27 | −1.03 | −0.74 |
BPh4− | 554 | 3.28 | 4.35 | — | — |
BISON | 479 | 7.04 | 6.29 | −0.13 | 0.23 |
BPFPB | 533 | 5.94 | 5.48 | — | — |
BBB | 580 | 2.64 | 3.81 | −0.64 | −0.46 |
4F-BBB | 497 | 3.84 | 4.84 | −0.57 | 0.51 |
BCB | 480 | 4.78 | 3.80 | 1.74 | 2.35 |
BNB | 551 | 1.89 | 3.97 | 0.05 | 0.27 |
BBPB | 593 | 3.08 | 3.95 | −0.29 | 0.03 |
BSB | 593 | 3.83 | 4.26 | 0.13 | 0.73 |
BOB | 520 | 5.67 | 5.37 | 1.28 | 1.72 |
BMB | 525 | 5.51 | 5.35 | −0.74 | 0.50 |
DCTA | 541 | 4.26 | 4.99 | 0.34 | 0.51 |
TIM | 493 | 4.02 | 4.65 | 0.41 | 0.44 |
TCP | 479 | 4.23 | 4.33 | 0.81 | 0.97 |
PCPI | 525 | 3.59 | 3.54 | 1.89 | 2.18 |
TDI | 516 | 4.05 | 4.74 | 0.35 | 0.59 |
TDBI | 516 | 4.09 | 4.37 | 0.64 | 0.82 |
TDPI | 505 | 4.64 | 4.21 | 1.53 | 1.74 |
PF6− | 567 | 7.85 | 7.36 | −1.36 | −0.57 |
AsF6− | 558 | 8.24 | 7.16 | 1.32 | 1.67 |
TFSI | 592 | 5.79 | 5.15 | 1.22 | 1.73 |
The introduction of further modifications, such as adding –CN groups to nitrogen atoms to form [(Z)-1,2-dicyanoethene-1,2-diyl]diazanidecarbonitrile ligands or even a higher substitution degree by replacing these nitrogen atoms by carbon atoms, results in even lower ΔEd values: 454 and 426 kJ mol−1, respectively, for BTCB and BBTCB. The very low cation–anion interaction strength obtained for the latter anion is in fact lower than for any other LIB electrolyte relevant ion-pair that we know of. The comparable lithium salts from the literature with anions based on a boron central atom, except the said BISON, reaches ΔEd's of 520–525 kJ mol−1 for BOB and BMB, thus similar but higher, but for many others ΔEd is closer to 600 kJ mol−1 such as for BF4− and B(CH3)4−. Also the Hückel salts are not able to compete; their ΔEd's range from 479 kJ mol−1 for the pyrolates e.g. TCP, via 493–516 kJ mol−1 for imidazolates e.g. TIM and TDI, to 541 kJ mol−1 for triazolates e.g. DCTA. As for a comparison with the standard PF6− anion and its ion-pair with Li+, ΔEd is 567 kJ mol−1.
As the lithium ion is quite small, thus rendering steric hindrance less likely, the nature of the significant decrease in the interaction energy can primarily be addressed by the anion charge distribution. One crucial parameter, the ligand electron withdrawing force, here measured by the resulting partial charge on the central boron atom, reveals the connecting atom to be the more important design feature, the effect largely stops after the first atom. All the resulting charge distributions for the novel anions are listed in Table 2. A fluorine ligand, a very strong electron withdrawing moiety, can cause an increase in the boron atom charge from +1.57 for BH4− to +2.36 for BF4−, while the oxygen connecting ligands result in ca. +2.30, as seen for BOB and BMB, and more medium effects are observed for carbon connecting ligands such as for B(CH3)4− resulting in a boron atom charge of +1.83. However, there are large differences between normal carbon connecting ligands such as –CH3 and –Ph (ca. +1.8) and –CN ligands (ca. +2.0), probably due to the very strong carbon–nitrogen triple-bond.
Anion | Core | Ligands | |||||
---|---|---|---|---|---|---|---|
B | –CN | –F | Conjugated | ||||
Backbone | –CN | –F | –CF3 | ||||
a Nitrile groups close to the boron atom. b Nitrile groups far from the boron atom. | |||||||
BDB | +2.12 | — | — | −1.03 | −0.26 | — | — |
DFDB | +2.25 | — | −0.82 | −1.02 | −0.29 | — | — |
DCDB | +2.06 | −0.75 | — | −1.03 | −0.26 | — | — |
BTMB | +2.12 | — | — | −1.15 | — | — | −0.21 |
DFTMB | +2.25 | — | −0.82 | −1.14 | — | — | −0.23 |
DCTMB | +2.06 | −0.76 | — | −1.16 | — | — | −0.20 |
BFB | +2.13 | — | — | −0.31 | — | −0.63 | — |
DFFB | +2.25 | — | −0.83 | −0.32 | — | −0.63 | — |
DCFB | +2.07 | −0.77 | — | −0.27 | — | −0.63 | — |
BTCB | +2.18 | — | — | −1.77 | +0.28a | — | — |
−0.19b | |||||||
DFTCB | +2.29 | — | −0.82 | −1.67 | +0.25a | — | — |
−0.23b | |||||||
DCTCB | +2.10 | −0.76 | — | −1.70 | +0.27a | — | — |
−0.21b | |||||||
BBTCB | +1.84 | — | — | −0.46 | −0.25a | — | — |
−0.22b | |||||||
DFBTCB | +2.14 | — | −0.81 | −0.40 | −0.29a | — | — |
−0.27b | |||||||
DCBTCB | +1.92 | −0.76 | — | −0.42 | −0.25a | — | — |
−0.24b |
The novel anions with a ligand coordination of the boron atom by four nitrogen atoms, as for BDB, BTMB, and BFB, provide an average electron withdrawing force resulting in a boron atom charge of ca. +2.12 (Table 2). This level appears to be the most appropriate, ensuring a fair charge distribution between the core and the ligands, balancing the need for proper delocalization and the risk of creating negatively charged parts/fragments attractive for lithium ion interaction. The change to BTCB by adding –CN groups connected to nitrogen atoms only slightly increases the observed effect (+2.18). Further replacement of these nitrogen atoms by carbon atoms, as for BBTCB, results in a decrease to a typical value for carbon connected ligands, +1.84.
As stated above, the part of the ligand not connecting the central boron atom seems to have no significant direct influence on the core charge. Thus we as a first approximation consider the core and ligand charge distributions independently. As the observed cation–anion interaction energies result from delocalization of the negative charge, comparing BFB, BTMB and BDB, the –F ligand at BFB has the strongest electron withdrawing force with a ligand charge of −0.63, and will thus possibly strongly interact with the lithium cation via this fluorine ligand. In turn, the –CF3 group of BTMB has a much weaker effect, with the charge of the trifluoromethyl ligand being equal to −0.21. However, inside the latter ligand there is a significant gradient of charge density: a highly positive carbon atom (+1.61) is surrounded by three strongly electron withdrawing fluorine atoms (−0.61). Thus also the –CF3 ligand can be expected to coordinate strongly to a lithium cation, especially if being able to use two adjacent fluorine atoms in a bi-dentate manner. In contrast, the partial charges of the –CN groups in BDB are −0.26 and this provides the lowest ΔEd, while the large charge density gradients, +0.91 at C and −1.17 at N, seem to be of less importance and at the same time the lithium cation coordination is mono-dentate.
The exchange of an aromatic system for a corresponding conjugated backbone connected to the tetrahedral boron atom seems to be advantageous for oxidation resistance. Application of the imidazole-derived ligand from the anions TDI and TIM, anions which exhibit ΔEv's of 4.05–4.09 V vs. Li+/Li0, to the boron-centered anions causes an increased stability; 5.80 V vs. Li+/Li0 for the BDB anion. Likewise, changing from a typical Hückel anion to a boron-centered anion using a pyrrole ring derived ligand results in an enhancement of ΔEv: 4.23 vs. 5.44 V vs. Li+/Li0, for TCP and BBTCB, respectively. Lithium salts with such high oxidation stabilities are very desirable as they allow for application of cathode materials with higher voltages, and thus improved energy density of the battery. However it should be noted that experimental values of the oxidation potential limit (Eox) for Hückel anions (TCP, TIM) are usually higher than predicted with the ΔEv approach, above 5 V vs. Li+/Li0.54 On the other hand, ΔEv's of boron centered anions (BOB, MOB, B(CH3)4−)38,55,56 overestimate Eox's, which are usually ca. 1 V lower.50 Despite these difficulties in correlation, the calculated ΔEv's for anions like BDB, DCDB, BTMB and DCTMB are high enough; 5.73–6.03 V vs. Li+/Li0, to consider their application to the cathode material LiNi0.5Mn1.5O4, requiring stability up to 4.7 V vs. Li+/Li0.53 On the other hand, also the salts with the lowest stabilities vs. oxidation, DFTCB and DCTCB with ΔEv's of 3.90–4.29 V vs. Li+/Li0, might still be able to work with standard cathode materials such as LiFePO4 with a mere need for stability up to 3.45 V vs. Li+/Li0.57
![]() | ||
Fig. 3 Reduction potentials for anions (circles) and ion-pairs (squares) of the novel lithium salts. |
The most stable anions/ion-pairs, i.e. the lowest reduction potentials were obtained for anions containing the (1,2-difluoroethane-1,2-diylidene)diazanide ligand, such as BFB, DFFB and DCFB, 1.01, 0.98 and 0.87 V vs. Li+/Li0 for their ion-pairs. Changing to conjugated ligands with –CN groups results in larger increases in the reduction ability, 1.84, 2.17, and 2.36 V vs. Li+/Li0 for the ion-pairs of BDB, BFBD and DCBD, respectively. Application of another –CN group to nitrogen atoms causes a little drop in Ered; 1.67 vs. Li+/Li0 for BTCB. The least stable anions are those containing the pyrrole-derived ligand – for BBTCB, DFTCB and DCBTCB, the reduction potentials all exceed 2.5 V vs. Li+/Li0, seemingly preventing any use in LIBs.
For the series of halide anions F−, Cl− and Br− their η values of 6.02, 4.89 and 4.37 eV50 correlate well with both their ΔEd's of 778, 645 and 612 kJ mol−1,45 and their ΔEv's of 1.71, 2.23 and 2.03 V vs. Li+/Li0.50 Our results show that the chemical hardness is in very good agreement with the degree of charge delocalization within the anion (Table 1). The application of fluorine atoms, both as the boron center and as part of the conjugated ligand, results in higher η values; 4.58, 5.34 and 5.50 eV for DFDB, BFB and DFFB, respectively, a little lower value is observed for ligands with the –CF3 group; 4.97 and 5.11 eV for BTMB and DDFTMB. The use of –CN ligands results in lower chemical hardness: down to 4.48 eV for BDB, but the application of up to four –CN groups per ligand results in very low chemical hardness, below 4 eV as shown for BTCB and BBTCB, 3.66 eV and 3.75 eV, respectively. The very low cation–anion interaction energies for the BBTCB and BTCB anions can be explained by this and in general such low levels of chemical hardness are rarely observed, even for completely conjugated and aromatic Hückel anions, only the PCPI pseudo-aromatic anion reaches as low as 3.59 eV. The common LIB anions, such as PF6− and TFSI, are above 7 and 5 eV, respectively.
This journal is © the Owner Societies 2016 |